pith. sign in

arxiv: 2606.27631 · v1 · pith:POECSFOMnew · submitted 2026-06-26 · ⚛️ nucl-th · nucl-ex

Oblate-prolate shape mixing and E0 transition in 28Si

Yasutaka Taniguchi , Masaaki Kimura This is my paper

Pith reviewed 2026-06-29 00:45 UTC · model grok-4.3

classification ⚛️ nucl-th nucl-ex
keywords shape coexistenceoblate prolate mixing28SiE0 transitionnuclear structureAMD GCMGogny force
0
0 comments X

The pith

In silicon-28 the ground state is dominated by the oblate shape with the prolate component limited to under 20 percent.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constrains the oblate-prolate shape mixing in the ground and low-lying states of 28Si using AMD and GCM calculations. By fitting to the charge radius, quadrupole moment, and E2 transition strengths while varying the Gogny force, it determines that the oblate component dominates the ground state. The prolate admixture is limited to less than 20 percent, and an upper bound is placed on the E0 transition strength. This helps clarify the degree of shape coexistence, which has been discussed for decades but remained poorly quantified.

Core claim

Oblate and prolate 0+ and 2+ configurations are obtained by antisymmetrized molecular dynamics combined with the generator coordinate method. Using these as basis states, the mixing amplitudes are constrained by simultaneously reproducing the measured charge radius, the quadrupole moment of the 2_1+ state, and the in-band and inter-band B(E2) values, with variation of the density-dependent term in the Gogny interaction. In the ground state the oblate component is dominant and the prolate component is limited to less than about 20%. For the 2_1+ state the allowed prolate component is smaller. An upper limit of rho^2(E0;0_3+ to 0_1+) less than or equal to 0.206 is obtained.

What carries the argument

Mixing amplitudes of oblate and prolate configurations from AMD+GCM constrained by experimental E2 strengths, charge radius and quadrupole moment.

If this is right

  • The prolate component in the 2_1+ state is smaller than in the ground state.
  • The low-lying 0+ states may exhibit substantial oblate-prolate mixing.
  • The E0 transition strength has an upper limit but is not tightly constrained.
  • A measurement of the inter-band E0 transition strength would provide a quantitative determination of the mixing amplitude.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This limited mixing suggests shape coexistence in 28Si is less pronounced than in some other nuclei with similar mass.
  • The method of constraining mixing via multiple observables could be applied to other sd-shell nuclei to test the robustness of shape coexistence predictions.
  • Future E0 measurements could directly test the predicted mixing amplitudes.

Load-bearing premise

The pure oblate and prolate configurations generated by AMD+GCM form an adequate basis whose mixing can be constrained solely by the chosen experimental observables and the variation of the Gogny density-dependent term.

What would settle it

A measurement showing the prolate component in the ground state exceeds 20 percent or finding rho^2(E0;0_3+ to 0_1+) larger than 0.206 would challenge the constrained mixing amplitudes.

Figures

Figures reproduced from arXiv: 2606.27631 by Masaaki Kimura, Yasutaka Taniguchi.

Figure 1
Figure 1. Figure 1: FIG. 1. Theoretical and experimental level scheme of low [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Proton radius (upper), Q-moments (middle), and [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Overlap regions of the constraints in the [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Constraints in the [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 2
Figure 2. Figure 2: The oblate and prolate basis states have different [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Predicted dimensionless E0 transition strength of the [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
read the original abstract

Background: oblate-prolate shape coexistence in $^{28}$Si has been discussed for decades, but the degree of shape mixing between these configurations remains poorly constrained. Purpose: We constrain the oblate-prolate mixing amplitudes in $^{28}$Si using available experimental information and discuss the inter-band E0 transition strength. Methods: Oblate and prolate $0^+$ and $2^+$ configurations are obtained by antisymmetrized molecular dynamics combined with the generator coordinate method. Using these configurations as the basis states, we constrain the mixing amplitudes by simultaneously reproducing the measured charge radius, the quadrupole moment of the $2_1^+$ state, and the in-band and inter-band $B(\mathrm{E2})$ values. The strength of the density-dependent term in the Gogny interaction is also varied within a reasonable range. Results: In the ground state, the oblate component is dominant, and the prolate component in the ground state is limited to less than about $20\%$. For the $2_1^+$ state, the allowed prolate component is smaller than that in the ground state. The present analysis does not tightly constrain the corresponding E0 transition strength, but an upper limit of $\rho^2(\mathrm{E0};0_3^+\rightarrow0_1^+) \lesssim 0.206$ is obtained. Conclusions: The low-lying $0^+$ states of $^{28}$Si may exhibit substantial oblate-prolate mixing. A measurement of the inter-band E0 transition strength would provide a quantitative determination of the mixing amplitude.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript uses antisymmetrized molecular dynamics combined with the generator coordinate method (AMD+GCM) to generate pure oblate and prolate 0+ and 2+ configurations in 28Si. These serve as a two-configuration basis whose mixing amplitudes are constrained by simultaneously reproducing the experimental charge radius, Q(2_1+), and selected in-band and inter-band B(E2) values, while also varying the strength of the density-dependent term in the Gogny interaction within a reasonable range. The central results are that the ground state is dominantly oblate with a prolate component limited to less than about 20%, the prolate admixture in the 2_1+ state is even smaller, and an upper limit ρ²(E0; 0_3+ → 0_1+) ≲ 0.206 is obtained from the allowed mixing range.

Significance. If the central claim holds, the work supplies quantitative bounds on long-discussed oblate-prolate shape coexistence in 28Si and yields a concrete, falsifiable upper limit on an inter-band E0 strength that can be tested experimentally. The explicit variation of one Gogny parameter within a stated range provides a limited robustness check on the mixing amplitudes, which is a positive feature of the analysis.

major comments (2)
  1. [Methods] Methods paragraph: the procedure used to determine the mixing amplitudes (e.g., whether a χ² minimization, grid search, or manual adjustment is employed, and how the three observables plus the Gogny variation are weighted) is not specified. This detail is load-bearing for the reported <20% prolate bound, because different fitting choices could alter the allowed mixing range.
  2. [Results] Results paragraph: the upper limit ρ²(E0;0_3+→0_1+) ≲ 0.206 is stated to follow from the allowed mixing amplitudes after the fit, yet no explicit expression or numerical steps connecting the mixing coefficients to the E0 matrix element are provided. Without this, it is impossible to verify whether the numerical value is robust or sensitive to the precise definition of the E0 operator within the two-configuration basis.
minor comments (2)
  1. [Abstract] The abstract lists the fitted observables but does not state the numerical range explored for the Gogny density-dependent strength; this range should be given explicitly.
  2. Notation for the E0 strength (ρ²) is introduced without a reference to the conventional definition used in the field; a brief reminder or citation would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation of our work and the constructive comments. We address the two major comments below and will revise the manuscript to incorporate the requested clarifications.

read point-by-point responses

  1. Referee: [Methods] Methods paragraph: the procedure used to determine the mixing amplitudes (e.g., whether a χ² minimization, grid search, or manual adjustment is employed, and how the three observables plus the Gogny variation are weighted) is not specified. This detail is load-bearing for the reported <20% prolate bound, because different fitting choices could alter the allowed mixing range.

    Authors: We agree that the fitting procedure requires explicit description. The mixing amplitudes were determined via a systematic grid search over the oblate-prolate mixing coefficients for the 0+ and 2+ states (in 5% steps), while stepping the density-dependent Gogny parameter over its stated range. Acceptable solutions were those simultaneously reproducing the experimental charge radius, Q(2_1+), and selected B(E2) values within their uncertainties (or 10% when uncertainties were not quoted). We will add a dedicated paragraph in the Methods section detailing this grid-search approach and the acceptance criteria, which directly supports the robustness of the <20% prolate limit. revision: yes

  2. Referee: [Results] Results paragraph: the upper limit ρ²(E0;0_3+→0_1+) ≲ 0.206 is stated to follow from the allowed mixing amplitudes after the fit, yet no explicit expression or numerical steps connecting the mixing coefficients to the E0 matrix element are provided. Without this, it is impossible to verify whether the numerical value is robust or sensitive to the precise definition of the E0 operator within the two-configuration basis.

    Authors: We acknowledge the omission. Within the two-configuration basis the E0 matrix element is ρ(E0;0_3+→0_1+) = |α_o β_p <O|E0|P> + α_p β_o <P|E0|O>| where α,β are the mixing amplitudes for the oblate (o) and prolate (p) components of the 0_1+ and 0_3+ states; the squared strength is then normalized by the standard factor 4π/Z²R⁴. The upper limit 0.206 arises from the maximum allowed off-diagonal mixing consistent with the fitted range. We will insert the explicit formula together with a short numerical example tracing the largest admissible mixing coefficients to the quoted bound, and we will note the definition of the E0 operator employed. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper generates oblate and prolate basis states via AMD+GCM, then constrains their mixing amplitudes by fitting to a set of experimental inputs (charge radius, Q(2_1^+), selected B(E2) values) while varying one Gogny parameter. From the resulting allowed mixing range it reports an upper bound on the E0 matrix element, which is not among the fitted observables. This is a standard constrained-parameter analysis whose output (the E0 limit) is independent of the inputs by construction; no equation reduces to a tautology, no fitted quantity is relabeled as a prediction, and no load-bearing self-citation chain is present. The central claim therefore remains non-circular.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that AMD+GCM generates sufficiently pure oblate and prolate basis states and that varying one Gogny parameter plus fitting mixing amplitudes to three observables is sufficient to constrain the physics.

free parameters (2)
  • oblate-prolate mixing amplitudes
    Determined by simultaneous fit to charge radius, Q(2_1+), and B(E2) values
  • strength of density-dependent term in Gogny interaction
    Varied within a reasonable range to explore model dependence
axioms (1)
  • domain assumption AMD+GCM configurations represent distinct pure oblate and prolate shapes suitable as basis states
    Invoked in Methods when constructing the mixing basis

pith-pipeline@v0.9.1-grok · 5819 in / 1402 out tokens · 28531 ms · 2026-06-29T00:45:12.551646+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

21 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

    Bohr, Rotational motion in nuclei, Reviews of Modern Physics 48, 365 (1976)

    A. Bohr, Rotational motion in nuclei, Reviews of Modern Physics 48, 365 (1976)

    1976

  2. [2]

    Bender, P.-H

    M. Bender, P.-H. Heenen, and P.-G. Reinhard, Self- consistent mean-field models for nuclear structure, Re- views of Modern Physics 75, 121 (2003)

    2003

  3. [3]

    Prochniak and S

    L. Prochniak and S. G. Rohozinski, Quadrupole collective states within the Bohr collective Hamiltonian, Journal of Physics G: Nuclear and Particle Physics 36, 123101 (2009)

    2009

  4. [4]

    Nikšić, D

    T. Nikšić, D. Vretenar, and P. Ring, Relativistic nu- clear energy density functionals: Mean-field and beyond, Progress in Particle and Nuclear Physics 66, 519 (2011)

    2011

  5. [5]

    Heyde and J

    K. Heyde and J. L. Wood, Shape coexistence in atomic nuclei, Reviews of Modern Physics 83, 1467 (2011)

    2011

  6. [6]

    P. E. Garrett, M. Zielińska, and E. Clément, An exper- imental view on shape coexistence in nuclei, Progress in Particle and Nuclear Physics 124, 103931 (2022)

    2022

  7. [7]

    Bernier and M

    J.-P. Bernier and M. Harvey, The 28si nucleus in the su3 model, Nuclear Physics A 94, 593 (1967)

    1967

  8. [8]

    S. D. Gupta and M. Harvey, The 28si nucleus in the hartree–fock model, Nuclear Physics A 94, 602 (1967)

    1967

  9. [9]

    G. D. Dang, Mixing of hartree–fock states in 28si, Nuclear Physics A 157, 231 (1970)

    1970

  10. [10]

    Kanada-En’yo, Deformations in n = 14 isotones, Physical Review C 71, 014303 (2005)

    Y. Kanada-En’yo, Deformations in n = 14 isotones, Physical Review C 71, 014303 (2005)

    2005

  11. [11]

    Taniguchi, Y

    Y. Taniguchi, Y. Kanada-En’yo, and M. Kimura, Cluster structures and superdeformation in 28Si, Physical Review C 80, 044316 (2009). 6

    2009

  12. [12]

    Chiba, Y

    Y. Chiba, Y. Taniguchi, and M. Kimura, Inversion dou- blets of reflection-asymmetric clustering in 28Si and their isoscalar monopole and dipole transitions, Physical Re- view C 95, 044328 (2017)

    2017

  13. [13]

    Frycz, J

    D. Frycz, J. Menéndez, A. Rios, B. Bally, T. R. Rodríguez, and A. M. Romero, Shell-model study of 28si: Coexistence of oblate, prolate, and superdeformed shapes, Physical Review C 110, 054326 (2024)

    2024

  14. [14]

    N. J. Stone, Table of nuclear electric quadrupole mo- ments, Atomic Data and Nuclear Data Tables 111–112, 1 (2016)

    2016

  15. [15]

    M. S. Basunia, Nuclear data sheets for a = 28 , Nuclear Data Sheets 114, 1189 (2013)

    2013

  16. [16]

    Deformed Base Antisymmetrized Molecular Dynamics and its Application to ^{20}Ne

    M. Kimura, Deformed-basis antisymmetrized molecular dynamics and its application to 20Ne, Physical Review C 69, 044319 (2004), arXiv:nucl-th/0311062

    work page internal anchor Pith review Pith/arXiv arXiv 2004

  17. [17]

    Dechargé and D

    J. Dechargé and D. Gogny, Hartree–fock–bogolyubov cal- culations with the D1 effective interaction on spherical nuclei, Physical Review C 21, 1568 (1980)

    1980

  18. [18]

    J. F. Berger, M. Girod, and D. Gogny, Time-dependent quantum collective dynamics applied to nuclear fission, Computer Physics Communications 63, 365 (1991)

    1991

  19. [19]

    Angeli and K

    I. Angeli and K. P. Marinova, Table of experimental nu- clear ground state charge radii: An update, Atomic Data and Nuclear Data Tables 99, 69 (2013)

    2013

  20. [20]

    Delaroche, M

    J.-P. Delaroche, M. Girod, J. Libert, H. Goutte, S. Hi- laire, S. Péru, N. Pillet, and G. F. Bertsch, Structure of even–even nuclei using a mapped collective Bohr hamil- tonian and the D1S Gogny interaction, Physical Review C 81, 014303 (2010)

    2010

  21. [21]

    Foris and C

    B. Foris and C. N. Papanicolas, Electron scattering and nuclear structure, Annu. Rev. Nucl. Part. Sci. 37, 133 (1987)

    1987