{"id":"af7ceceb-7fef-42cc-91b9-e97869911737","arxiv_id":"2507.20275","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Genuine gamma and beta vibrational bands are identified in 166Er and 162Dy above the gamma band using a phase-specified T-plot of Monte Carlo Shell Model wavefunctions.","lead":"Large-scale shell-model calculations identify genuine vibrational states above the gamma band in the heavy deformed nuclei 166Er and 162Dy. The new phase-specified T-plot analysis distinguishes these vibrations from rotations and predicts B(E2) values that future experiments can test.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"PT-plot phase convention and orthonormalization stability are not demonstrated; the sign pattern used to assign gamma/beta vibrations may be an expansion artifact rather than a state property.","rationale":"The paper's PT-plot is a creative diagnostic and the qualitative comparison with the harmonic oscillator is suggestive, but the central assignments depend on a sign pattern whose definition carries a phase-gauge freedom. The reader's weakest_assumption already points to the PT-plot interpretation; I sharpen that concern by identifying the specific non-invariance of the plotted overlap under allowed phase/unitary transformations of the norm-matrix eigenvectors and of the MCSM eigenstates. The paper does not provide a quantitative node-counting measure or a demonstration that the sign pattern is stable under these transformations, so the assignments 0+3 as gamma vibration in 166Er and 0+6 as gamma vibration in 162Dy are not yet established. The cited B(E2) values and the wider PES support for the 162Dy beta vibration are real independent evidence, which is why the paper should not be rejected outright. The concrete numerical test proposed here would settle whether the sign pattern is a property of the state or of the expansion convention, and therefore whether the conditional acceptance should become a full acceptance.","tokens_in":9102,"tokens_out":12689,"duration_ms":161137,"concrete_test":"Recompute the PT-plot for 166Er 0+3 with (a) the norm-matrix eigenvectors multiplied by independent fixed signs (e.g., (-1)^i) and (b) the global phase of the 0+3 eigenvector fixed by Re<0+1|0+3> > 0. If the red/yellow boundary in Fig. 2(c) remains at the same gamma bin (up to an overall color swap), the node is robust; if the boundary shifts or washes out, the vibrational assignment is an artifact of the expansion. As a cross-check, repeat the analysis with the 0+2 state as reference and verify that a physical node is reference-independent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests entirely on reading a sign change in the PT-plot as a nodal excitation in gamma or beta_2. In the End Matter construction, the color of each circle is the real part of the overlap <tilde phi_i^(R)|tilde phi_i^(k)> between K-integrated states formed after diagonalizing the norm matrix (Eq. 3). This overlap is not an observable: its phase is sensitive to the arbitrary global phase of each MCSM eigenvector and to the unitary freedom of the norm-matrix eigenvectors when eigenvalues are near-degenerate. Footnote [61] only asserts that the orthogonalization mixes states with very close deformation parameters, so the beta_2/gamma labels of the circles are preserved; it does not address the stability of the relative phases that define the red/yellow pattern. If near-degenerate eigenvalues occur, the |tilde phi_i^(k)> can rotate and the sign pattern of Fig. 2(c) can change without changing any energy or B(E2). Moreover, the histograms of Fig. 3 are expansion coefficients in a non-orthogonal discrete basis, not a collective wave function; a red/yellow separation is necessary but not sufficient for a one-node vibration, because orthogonality to a nodeless 0+1 ground state can itself force sign oscillations in such an expansion. Without a phase-convention-independent definition of the node and a quantitative node count, the identification of 0+3, 2+5, and 0+6 as vibrational states is not established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a large-scale QVSM/MCSM study of 166Er and 162Dy and introduces an extension of the T-plot, termed the phase-specified T-plot (PT-plot), in which the real part of the overlap between K-integrated basis states of a reference and a target eigenstate is used to color circles at each basis vector's (β2, γ) position. Based on a red/yellow sign pattern and a harmonic-oscillator analogy, the authors identify the 0+3 state of 166Er as a K^P = 0+ gamma vibration built on the ground state, the 2+5 state as a gamma vibration built on the 2+2 state, the 0+2 state of 162Dy as a beta vibration, and the 0+6 state as a gamma vibration; shape-coexistence bands are also discussed. Energies and selected B(E2) values are compared with experiment.","tokens_in":9434,"tokens_out":7741,"duration_ms":87127,"significance":"The result is potentially significant for the nuclear structure of heavy deformed nuclei: it argues that genuine vibrational modes appear above the rotational gamma band, and that the traditional double-gamma-phonon 0+ interpretation may be a single gamma-vibrational phonon in a triaxial context. The underlying QVSM/MCSM calculations are state of the art (Hilbert spaces of dimension about 10^33), and the energies and B(E2) values are genuine outputs of the shell-model Hamiltonian, not fitted to the vibrational observables; some predictions have experimental counterparts, e.g., B(E2; 4+γγ → 2+2) ≈ 8.4 W.u. in 162Dy and the predicted strong decay of the 0+3 state in 166Er to the 2+2 state. The novelty is high if the identification holds. However, the central classification currently rests on a qualitative visual pattern, so the significance is conditional on making that pattern quantitative and phase-convention independent.","major_comments":[{"comment":"The PT-plot color is defined as the real part of ⟨tilde φ_i^(R)|tilde φ_i^(k)⟩, and this quantity is not invariant under the arbitrary global phase that each energy eigenvector |Ψ^k⟩ carries, nor under the residual unitary freedom of the norm-matrix eigenvectors when eigenvalues are near-degenerate. Footnote [61] asserts only that the orthogonalization mixes states with very close deformation parameters, so that β2 and γ labels are preserved; it does not address the stability of the relative phases that define the red/yellow pattern. If the phase of |Ψ^k⟩ is flipped, every circle changes color; if near-degenerate norm eigenvalues are present, the tilde basis can rotate and the pattern of Fig. 2(c) can change without changing any energy or B(E2). The authors should demonstrate, e.g., by diagonalizing the norm matrix with random degenerate-subspace rotations and by testing the effect of global sign conventions, that the reported sign change is a robust property of the state rather than a gauge artifact.","section":"End Matter, Eq. (3); Fig. 2(c); footnote [61]"},{"comment":"The red/yellow histograms are expansion coefficients in a non-orthogonal, discrete MCSM basis, not a collective wavefunction in γ or β2. A sign change in these coefficients is necessary but not sufficient for a one-node vibration: orthogonality to a nodeless 0+1 ground state can itself force sign oscillations in such an expansion, and the binning/normalization procedure adds additional convention dependence. The harmonic-oscillator analogy is qualitative; the paper reports no quantitative measure of the node position, no count of sign changes, no significance test against a null pattern (e.g., random signs with the same ground-state overlap), and no check that the pattern survives a different basis truncation. The authors should provide an invariant node definition, for example by constructing a collective wavefunction in a suitably orthonormalized coordinate projection and counting its zeros, or by computing the overlap of the candidate state with a one-phonon excitation operator acting on the reference state.","section":"Fig. 3(b)–(e) and the paragraph comparing them with Fig. 3(a)"},{"comment":"The assignment of the 2+5 state as a gamma vibration 'on top of the 2+2 state' uses the 2+2 state as reference because 2+1 is a member of the ground band. But the same analysis is not applied to the rotational members of the proposed vibrational band, and no B(E2) or E0 transition connecting the 2+5 state to the 2+2 band is given to corroborate the phonon picture. Since the reference-state choice changes which states are compared, the authors should show that the node pattern is stable when the reference state is varied within the bands (e.g., using the 4+ member of the ground band as an alternative reference) and that the proposed vibrational band members share the same PT-plot signature.","section":"Fig. 2(e) and the discussion of the 2+5 state"}],"minor_comments":[{"comment":"The acronym 'PT-plot' is introduced, but 'T-plot' is used alone in several places (e.g., End Matter, 'The T-plot displays...'), which is confusing; please use the full name consistently or state explicitly when the standard T-plot is meant.","section":"Throughout"},{"comment":"The superscript (k) in |tilde φ_i^(k)⟩ and the eigenstate index k in |Ψ^k⟩ are the same symbol but refer to different objects; this dual use should be disambiguated, for instance by using a different label for the K-integrated basis states.","section":"End Matter, Eq. (3)"},{"comment":"'one-and-half HO shell' should read 'one-and-a-half HO shell'; also 'a` la A. Bohr' should be typeset as 'à la A. Bohr'.","section":"Main text, first paragraph"},{"comment":"The criterion for omitting experimental and theoretical levels is not stated; since the comparison of band structures depends on which levels are included, the selection rule should be specified.","section":"Fig. 1 and Fig. 4 captions"},{"comment":"Footnote [61] is placed at the end of the T-plot definition but is not referenced in the main text; consider inserting the reference at the first use of the T-plot, as it is directly relevant to the interpretation of the circle positions.","section":"Footnote [61]"}],"recommendation":"major_revision","confidential_remarks":"The underlying MCSM calculations are credible and the paper addresses a timely question, but the central claim is currently not established in a falsifiable way from the presented material. I would encourage the editor to request a revision that adds a quantitative, phase-invariant node analysis and a stability test of the PT-plot sign pattern. If the authors can supply those, the paper would be a strong addition to the literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper introduces the Phase-specified T-plot (PT-plot), a genuinely new diagnostic that adds amplitude phase information to the usual T-plot, and uses it to argue that specific low-lying states in 166Er and 162Dy are one-node vibrational excitations built on triaxial rotors. The underlying QVSM-MCSM calculations are state of the art, and the B(E2) predictions give experimental hooks. But the central assignment depends on reading sign changes in a non-orthogonal basis expansion as nodal structure, and that interpretive step is not yet backed by a phase-invariant or quantitative measure.\n\nWhat is new and good: The PT-plot is a real step forward: showing phase on the deformation plane can expose patterns that magnitude-only T-plots miss. The calculations are large-scale and the paper compares energies and B(E2) values to experiment, with some nice agreements, e.g., the 4+ to 2+ transition in 162Dy. The link to the earlier triaxial-rotor interpretation of the gamma band is coherent, and the paper is honest about where it is predicting rather than explaining.\n\nThe soft spot is exactly the one the stress-test note raises. The overlap used to color PT-plot circles, Eq. (3), is not an observable. It depends on the global phase of each MCSM eigenvector and on the unitary freedom in diagonalizing the norm matrix when eigenvalues are near-degenerate. Footnote 61 only says the orthogonalization mixes states with close deformation parameters; it does not show that the relative phases defining red versus yellow are stable. The sign pattern in Fig. 2(c) could in principle change without any energy or B(E2) changing. Moreover, the histograms in Fig. 3 are expansion coefficients in a non-orthogonal discrete basis, not a collective wavefunction. Orthogonality to the nodeless ground state can itself force sign oscillations, so red/yellow separation is necessary but not sufficient for a one-node vibration. A quantitative node count, or a phase-convention-independent observable, is missing. The harmonic-oscillator analogy is heuristic.\n\nThat said, the energies and B(E2) values are genuine outputs, not fitted to the target observables. The assignments are plausible, and the paper does not overclaim beyond what the pattern suggests. The selection of displayed states could bias the visual impression, but the appendix at least shows some mixed cases.\n\nWho is this for: nuclear structure physicists interested in collective modes and large-scale shell-model calculations. It deserves a serious referee: the method is new, the calculations are heavy, and the predictions are testable. A referee should push for a quantitative definition of the vibrational node and robustness tests of the PT-plot sign pattern. I would engage with it, expecting the interpretive claim to need major revision or careful strengthening.","headline":"New PT-plot diagnostic plus state-of-the-art MCSM calculations make a plausible but not yet solid case for gamma and beta vibrations in 166Er and 162Dy; the vibrational assignments rest on phase patterns whose gauge invariance is not demonstrated.","tokens_in":9941,"tokens_out":3977,"would_cite":false,"duration_ms":40068,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Genuine beta and gamma vibrations are identified above the rotational gamma band in 166Er and 162Dy.","keywords":["nuclear vibrations","gamma vibration","beta vibration","shape coexistence","phase-specified T-plot","Monte Carlo Shell Model","166Er","162Dy"],"falsifier":"Measure the gamma-decay branch from the $0^+_3$ state in $^{166}$Er: the calculation predicts $B(E2;0^+_3\\rightarrow2^+_2)=9.8$ W.u., so a much weaker experimental transition would contradict the assignment; equally, recompute the PT-plot with a different reference state or with an enlarged basis and check whether the red-to-yellow sign boundary in $\\gamma$ stays at the same deformation.","tokens_in":8931,"feed_emoji":"⚛️","tokens_out":8567,"duration_ms":95828,"temperature":0.7,"pith_summary":"Building on the finding that low-lying gamma bands in heavy deformed nuclei are rotational $K^P=2^+$ excitations of triaxial states, this paper asks where genuine vibrational excitations live. Analysing $^{166}$Er and $^{162}$Dy with ultra-large configuration-interaction calculations, it identifies the $0^+_3$ state in $^{166}$Er as the $K^P=0^+$ gamma-vibrational band head and the $2^+_5$ state as a gamma vibration built on the triaxial $2^+_2$ state; in $^{162}$Dy the $0^+_2$ state is a $\\beta$ vibration and the $0^+_6$ state a gamma vibration. The identification uses a new phase-specified T-plot in which the sign of each basis-vector amplitude is plotted against a reference state, so that a sign change across the $\\gamma$ or $\\beta_2$ coordinate signals a node, exactly as in a harmonic-oscillator wave function. If the picture holds, genuine vibrational modes do exist in heavy deformed nuclei, but they sit above the rotational gamma band, and shape-coexistence states can be interleaved at lower energies.","feed_headline":"Genuine beta and gamma vibrations found above the gamma band","feed_subtitle":"A sign-sensitive T-plot locates one-node waves in 166Er and 162Dy, replacing old labels with testable predictions.","key_machinery":"The key instrument is the phase-specified T-plot (PT-plot). A T-plot places a circle at the $(\\beta_2,\\gamma)$ deformation of each basis vector, with circle area proportional to the squared amplitude in the eigenstate; the PT-plot additionally colours each circle by the sign of the real part of the amplitude relative to a chosen reference state, computed through overlaps of the orthonormalised basis vectors. Because the amplitudes are predominantly real for the states of interest, the red/yellow sign pattern is meaningful. The work of the PT-plot is to display a node: when a state's amplitudes switch sign along $\\gamma$ or $\\beta_2$, the state is read as the first excited wave function of a harmonic oscillator in that coordinate, i.e. a one-phonon vibrational excitation. The reference state is chosen state by state; for example, the $2^+_5$ state is plotted against the $2^+_2$ state, not the $2^+_1$ state, because the $2^+_2$ state is the rotational partner of the ground state and the $2^+_1$ is not.","core_discovery":"On the paper's own terms, the central discovery is that $\\beta$- and gamma-vibrational modes are real, identifiable excitations in strongly deformed heavy nuclei, and that the traditional labelling missed them because the low-lying gamma band is actually a rotational excitation of a triaxial intrinsic state. In $^{166}$Er the $0^+_3$ state shows a one-node amplitude pattern in the $\\gamma$ coordinate relative to the $0^+_1$ ground state, making it the $K^P=0^+$ gamma-vibrational band head; the $2^+_5$ state shows the same pattern relative to the $2^+_2$ triaxial rotational state, making it a gamma vibration on top of that state. In $^{162}$Dy the $0^+_2$ state shows a one-node pattern in the $\\beta_2$ coordinate (a $\\beta$ vibration), and the $0^+_6$ state shows a gamma vibration. The authors further find that the $K^P=0^+$ vibrational band lies below the $K^P=2^+$ vibrational band because of the $K$-splitting of the triaxial intrinsic state, and that the $0^+$ gamma-vibrational head decays strongly to the $2^+_\\gamma$ state, a decay that could be misread as a double-gamma-phonon signature.","pith_inferences":["If the PT-plot sign-change criterion is a true nodal diagnostic, the same method could be applied across the deformed rare-earth and actinide regions; one testable prediction is that a low beta band will appear only where the potential energy surface is wide along $\\beta_2$, and will be absent where it is narrow.","A natural generalisation is to require the node location to be stable when several different reference states are used for the same eigenstate; that would turn the current diagnostic into a quantitative test of anharmonicity.","If these assignments are right, the traditional beta- and gamma-phonon counting in deformed nuclei needs revision: measured $B(E2)$ branchings that have been attributed to double-phonon states may in some cases be single-phonon decays from a $K^P=0^+$ gamma vibration."],"forward_implications":["The $0^+_3$ state in $^{166}$Er should be re-labelled as a $K^P=0^+$ gamma vibration rather than a double-gamma-phonon candidate; its strong decay to the $2^+_2$ state follows from single-phonon vibrational character.","The $0^+_2$ state in $^{162}$Dy is a beta vibration whose rotational band includes the $2^+_3$ state, while the $0^+_6$ state is a gamma vibration on top of the ground state.","Genuine vibrational band heads in these nuclei typically appear above the rotational gamma band, with the notable exception of the $^{162}$Dy beta band, which sits low.","Shape-coexistence states, such as the $0^+_2$ state in $^{166}$Er, can appear at even lower energies than vibrational band heads, so spectroscopy in this region must separate vibrational from coexisting deformed configurations.","In $^{166}$Er the $2^+_5$ state is a prediction of a gamma vibration on top of the triaxial $2^+_2$ band, giving a concrete experimental target."],"supporting_citations":[{"why":"Establishes that the low-lying gamma band is a rotational $K^P=2^+$ excitation of a triaxial state, the premise that forces a search for true vibrations.","marker":"[12]"},{"why":"Provides the self-organization mechanism and systematic context for heavy deformed nuclei used to interpret low-lying coexisting states.","marker":"[13]"},{"why":"Presents the quasi-particle vacua shell model used to solve the many-body problem in the ultra-large configuration space.","marker":"[14]"},{"why":"The Monte Carlo Shell Model framework whose basis vectors and projection are extended with additional vectors to target vibrational states.","marker":"[15, 16]"},{"why":"Supplies the proton-proton and neutron-neutron parts of the effective interaction used in the calculation.","marker":"[22]"},{"why":"Supplies the $V_{MU}$ interaction completing the proton-neutron channel.","marker":"[23]"},{"why":"Compiled experimental levels used as the comparison data in the energy-level figures.","marker":"[30]"},{"why":"Introduces the orthonormalized version of the T-plot that the present PT-plot extends.","marker":"[52]"},{"why":"Recent measured $B(E2;4^+\\rightarrow2^+)$ in $^{162}$Dy that the calculation reproduces, supporting the PT-plot assignments.","marker":"[56]"}],"fun_headline_variants":["Real vibrations found above the gamma band","True beta and gamma modes found in deformed nuclei","T-plot reveals genuine nuclear vibrations","Vibrational modes above the gamma band are real"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a sign change in the PT-plot, measured against a chosen reference state, truly marks a node of the collective wave function in the $\\gamma$ or $\\beta_2$ coordinate; if the reference-state convention or the orthogonalisation is unstable, the vibrational assignment would not be established.","fun_headline_variants_meta":{"raw":{"variants":["Real vibrations found above the gamma band","True beta and gamma modes found in deformed nuclei","T-plot reveals genuine nuclear vibrations","Vibrational modes above the gamma band are real"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000564,"raw_usage":{"total_tokens":2728,"prompt_tokens":1053,"completion_tokens":1675,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":669,"completion_tokens_details":{"reasoning_tokens":1619}},"tokens_in":669,"tokens_out":1675,"duration_ms":14869,"temperature":1.0,"reasoning_tokens":1619,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T13:40:40.377996+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the gamma-decay branch from the $0^+_3$ state in $^{166}$Er: the calculation predicts $B(E2;0^+_3\\rightarrow2^+_2)=9.8$ W.u., so a much weaker experimental transition would contradict the assignment; equally, recompute the PT-plot with a different reference state or with an enlarged basis and check whether the red-to-yellow sign boundary in $\\gamma$ stays at the same deformation.","supporting_citations":[{"cited_title":"Phase-specified T-plot (PT- plot)","cited_arxiv_id":null,"evidence_quote":"Establishes that the low-lying gamma band is a rotational $K^P=2^+$ excitation of a triaxial state, the premise that forces a search for true vibrations."},{"cited_title":"Otsuka, Y","cited_arxiv_id":null,"evidence_quote":"Provides the self-organization mechanism and systematic context for heavy deformed nuclei used to interpret low-lying coexisting states."},{"cited_title":"Otsuka, Y","cited_arxiv_id":null,"evidence_quote":"Presents the quasi-particle vacua shell model used to solve the many-body problem in the ultra-large configuration space."},{"cited_title":"Otsuka, T","cited_arxiv_id":null,"evidence_quote":"Supplies the $V_{MU}$ interaction completing the proton-neutron channel."},{"cited_title":"Selset al., Phys","cited_arxiv_id":null,"evidence_quote":"Compiled experimental levels used as the comparison data in the energy-level figures."},{"cited_title":"Otsuka and Y","cited_arxiv_id":null,"evidence_quote":"Introduces the orthonormalized version of the T-plot that the present PT-plot extends."}],"review_version":1}