{"id":"ca2ebb97-54ea-438d-ba62-0dc9d43be09c","arxiv_id":"2506.13354","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Finite-temperature relativistic QRPA calculations predict that low-energy E1 strength in nickel isotopes grows up to 2.5 times when heated to 2 MeV, with tables provided for experimental comparison.","lead":"What happens to the weak 'pygmy' electric dipole resonance in nickel nuclei when they are hot? This paper calculates it for eight nickel isotopes up to 2 MeV and predicts the hot strength grows up to 2.5 times, giving experimental teams concrete numbers to test.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Headline numbers depend on ad hoc Lorentzian width and no continuum treatment; sensitivity to Γ and basis is undocumented, so the 2.5× enhancement and 62Ni benchmark are not yet robust.","rationale":"The reader identified the neglect of beyond-RPA correlations and damping as the weakest assumption. I sharpen this to the concrete and testable issue of the fixed Lorentzian width and the missing continuum treatment, which directly govern the headline integrated strengths. The sign error in Eq. (6) compounds the verification problem. Because these are addressable with additional numerical checks and do not necessarily invalidate the underlying thermal-unblocking mechanism, the conditional verdict stands: the quantitative benchmarks should be accompanied by sensitivity analysis before being treated as validated.","tokens_in":14682,"tokens_out":10335,"duration_ms":101017,"concrete_test":"Recompute SEWS for 62Ni and 66Ni at T=0, 1.6, and 2 MeV using Γ=0.5, 1.0, and 2.0 MeV, and also with a larger basis or a continuum-QRPA (e.g., finite-temperature continuum QRPA of Refs. [47,48]). If the 0–12 MeV SEWS at T=2 MeV changes by more than ~20% across Γ or basis, or if the 62Ni T=1.6 MeV 8–12 MeV SEWS leaves the 3.5–4.5% range for reasonable Γ, the factor-2.5 claim and the \"excellent agreement\" are not robust. Fixing the Eq. (6) denominator and reproducing Fig. 1 would also confirm the calculations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Load-bearing concern: the quantitative outputs — the up-to-2.5× enhancement of 0–12 MeV E1 strength at T=2 MeV and the SEWS=4.32% \"excellent agreement\" for 62Ni at T=1.6 MeV — are calculated from a discrete FT-RQRPA spectrum smoothed with a fixed Lorentzian width Γ=1.0 MeV (Eq. 6). The integrated SEWS is sensitive to this Γ because the window boundaries cut Lorentzian tails, and the 0–12 MeV window extends well above the neutron separation energy in neutron-rich Ni (S_n≈6–8 MeV), where a discrete-basis RQRPA misses continuum and energy-dependent spreading effects. No Γ-sensitivity study or continuum-QRPA comparison (e.g., Refs. [47,48]) is provided, and the 62Ni benchmark is quoted for a post hoc 8–12 MeV window with the same Γ, so the agreement with the preliminary ~4% experimental value may be partly an artifact of the chosen smoothing and window. In addition, Eq. (6) as printed has a sign error in the denominator (minus instead of plus), so the published expressions cannot reproduce the plotted strength functions, leaving the numerical results unverifiable from the paper alone.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports finite-temperature relativistic quasiparticle random-phase approximation (FT-RQRPA) calculations of the isovector E1 strength in the nickel isotopes 56–70Ni for temperatures T = 0–2 MeV, using the DD-PCX point-coupling energy density functional. The central claims are that the integrated pygmy dipole strength in the E = 0–12 MeV window grows with temperature, reaching up to 2.5 times its zero-temperature value at T = 2 MeV in neutron-rich isotopes; that nuclei near N ≈ Z acquire low-energy E1 strength only at finite temperature; and that the calculated value SEWS = 4.32% for 62Ni at T = 1.6 MeV in the E = 8–12 MeV interval agrees with a preliminary experimental value of about 4%. The paper also provides a microscopic analysis of the participating two-quasiparticle configurations and, as supplementary material, tabulated SEWS and cumulative B(E1) values intended for direct comparison with experiments.","tokens_in":14941,"tokens_out":5192,"duration_ms":49899,"significance":"If the predicted thermal enhancement is confirmed, it would establish hot pygmy dipole strength as a sizeable, temperature-driven effect in medium-mass nuclei, with direct implications for the interpretation of ongoing fusion-evaporation experiments and for radiative neutron-capture rates in stellar environments. The paper has clear strengths: the T = 0 limit is benchmarked against existing 68,70Ni data, the FT-RQRPA calculations are self-consistent and not fitted to the hot-PDS data, and the supplementary tables provide quantitative, falsifiable benchmarks. However, the robustness of the headline numbers currently hinges on the adopted Lorentzian smoothing width and on a discrete-basis treatment of an energy region that extends above the neutron separation energy, so the quantitative claims are conditional on additional sensitivity checks.","major_comments":[{"comment":"The Lorentzian averaging formula as printed has a sign error in the denominator: it reads (E - E_w)^2 - Γ^2/4 instead of (E - E_w)^2 + Γ^2/4. The printed expression is not a normalized Lorentzian and can change sign, so it cannot reproduce the strength functions shown in Figs. 1–3 and 6. Because these figures and the derived SEWS values are the basis of the central quantitative claims, the formula must be corrected and the normalization stated explicitly.","section":"Theoretical framework, Eq. (6)"},{"comment":"The claimed 'excellent agreement' for 62Ni is not robust as presented: SEWS = 4.32% is quoted for the E = 8–12 MeV window, which is selected after the preliminary experimental value (approximately 4%, with no uncertainty stated) is known, while the preceding analysis uses E = 0–12 MeV. Since the integrated SEWS is obtained from Lorentzian-smoothed discrete states with a fixed width Γ = 1.0 MeV, the result depends on Γ through the window boundaries. The manuscript provides no Γ-sensitivity study and no basis-convergence test, so the agreement may be partly an artifact of the chosen smoothing and window. Please add such a study and also quote the 0–12 MeV SEWS value for 62Ni.","section":"Results and discussion, 62Ni benchmark and Fig. 3"},{"comment":"For neutron-rich nickel isotopes the neutron separation energy lies in the range S_n ≈ 6–8 MeV, so a substantial part of the E = 0–12 MeV interval lies in the continuum. The FT-RQRPA as applied here uses a discrete basis and omits coupling to the continuum as well as energy-dependent spreading effects. The paper cites continuum finite-temperature QRPA implementations (Refs. [47,48]) but does not compare with them. The magnitude of the reported up-to-2.5× enhancement and of the 62Ni benchmark could be affected by these omitted contributions, and the manuscript should at least quantify the expected effect, for example by comparing with a continuum-QRPA result for one representative isotope.","section":"Results and discussion, continuum effects"}],"minor_comments":[{"comment":"The caption states that panel (a) shows SEWS as a function of neutron number N and panel (b) as a function of temperature T, whereas the text describing Fig. 3(a) and 3(b) states the opposite ordering. Please correct the inconsistency.","section":"Fig. 3 caption"},{"comment":"The caption contains a duplicated unit: 'E = 0–6 MeV and 6 < E ≤ 12 MeV MeV'. Please fix the typographical error.","section":"Fig. 6 caption"},{"comment":"The notation B(EJ, E_w) = |Σ_{cd}(b^π_cd(E_w) + b^ν_cd(E_w))|² is ambiguous because B(EJ, E_w) already denotes the reduced transition probability from Eq. (4). Please clarify that the quantities b^π and b^ν are the partial amplitudes and that their coherent sum gives the total E1 transition amplitude for the selected state.","section":"Eq. (5)"},{"comment":"The statement that the strength 'increases by up to a factor of 2.5' would be more informative if the specific isotope(s) and the exact energy/temperature conditions were stated explicitly; consider adding this detail.","section":"Abstract and Conclusion"},{"comment":"Reference [6] appears incomplete, and some conference contributions (Refs. [38–41]) may benefit from additional bibliographic information (e.g., page numbers or DOIs where available).","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a timely topic and provides useful benchmark tables, but the validation of the central numbers is currently overstated relative to the uncertainties introduced by the Lorentzian smoothing, the post hoc energy window, and the discrete-basis treatment above the neutron separation energy. I do not see grounds for rejection; the concerns are addressable with additional sensitivity studies and a corrected Eq. (6). The authors should also temper the phrase 'excellent agreement' when comparing to a preliminary experimental value with no quoted uncertainty."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one for the systematic FT-RQRPA survey, not for the headline numbers. The paper gives the first consistent look at low-energy E1 strength across 56-70Ni from T=0 to 2 MeV, with tabulated SEWS and cumulative B(E1) as functions of temperature and energy window. That is genuinely useful for the ongoing fusion-evaporation campaigns, and the zero-temperature check against 68,70Ni data is a reasonable sanity check. The mechanism—thermal unblocking of single-particle states, mostly neutron configurations, with pairing collapse at Tc sharpening the effect—is clearly presented and consistent with earlier FT-RPA/QRPA work. Credit where due: the calculation is not fitted to any HPDS data point, so the central prediction is independent.\n\nThe soft spots are real but mostly fixable. Eq. (6) as printed has a minus sign in the Lorentzian denominator, so the strength functions cannot be reproduced from the paper. That's a typo, but it has to be fixed. More important: all quantitative outputs use a fixed smoothing width Γ=1.0 MeV, and there is no sensitivity study. Integrated SEWS over windows like 0-12 MeV or 8-12 MeV cuts the tails of the Lorentzians, so the factor of 2.5 and the 4.32% for 62Ni could shift with Γ. The 62Ni comparison is also a bit too kind: the 8-12 MeV window was adopted after the preliminary experimental value was known, and the agreement is with a preliminary number. I would soften 'excellent agreement' accordingly. The discrete-basis RQRPA also has no continuum treatment, so above the neutron separation energy the predicted strength functions should be viewed as semi-quantitative; for integrated B(E1) in the low-energy region this is probably acceptable, but it deserves a caveat. Table II shows 70Ni underpredicted by roughly 40%, which is a reminder that the model's neutron-rich performance is not perfect.\n\nNone of this breaks the qualitative message: thermally occupied states do open new E1 channels and the pygmy strength does grow with temperature. But the quantitative benchmarks are not yet robust enough to be called validated. This deserves a serious referee, not a desk reject. If the authors fix the sign, add a Γ-sensitivity check or at least a discussion of how the integrated results depend on the smoothing width, and recalibrate the language about experimental agreement, it becomes a solid reference for the HPDS experiments.","headline":"A useful benchmark survey of hot pygmy dipole strength in Ni, but the quantitative claims need a corrected Lorentzian formula and a smoothing-width sensitivity study before the 'excellent agreement' with experiment should be taken at face value.","tokens_in":15471,"tokens_out":3124,"would_cite":true,"duration_ms":31841,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V35"],"pacs":["21.60.Jz","24.30.Cz"],"model":"deepseek-v4-flash","headline":"This paper predicts that heating nickel isotopes to 2 MeV multiplies low-energy electric dipole strength by up to 2.5 times, and provides experiment-ready benchmarks for the hot pygmy dipole strength.","keywords":["hot pygmy dipole strength","finite-temperature relativistic quasiparticle RPA","nickel isotopes","E1 strength","thermal unblocking","pygmy dipole resonance","nuclear astrophysics","fusion-evaporation experiments"],"falsifier":"Measure the gamma-ray strength function of nickel-62 or nickel-66 from fusion-evaporation reactions at temperatures of 1.5 to 2 MeV: if the integrated E1 strength below 12 MeV does not rise toward 2.5 times the zero-temperature value, or if the 8-12 MeV SEWS in nickel-62 at $T=1.6$ MeV is far from 4.32%, the thermal-enhancement claim would be falsified.","tokens_in":14454,"feed_emoji":"🔥","tokens_out":12303,"duration_ms":104485,"temperature":0.7,"pith_summary":"At temperatures up to 2 MeV, the low-energy electric dipole response of nickel isotopes is predicted to change dramatically: neutron-rich isotopes develop a hot pygmy dipole strength that can reach 2.5 times its zero-temperature value, while near-symmetric nuclei acquire low-energy dipole strength that is absent at zero temperature. The mechanism is thermal unblocking: heating smears the Fermi surface, partially populating higher single-particle orbitals, and opens new two-quasiparticle transition channels. The paper provides tables of energy-weighted E1 strength and cumulative B(E1) values for nickel-56 through nickel-70 over temperature and energy windows, intended as direct benchmarks for ongoing fusion-evaporation experiments. These benchmarks matter because low-energy dipole strength feeds radiative neutron-capture rates in astrophysical environments and probes the symmetry energy.","feed_headline":"Heating nickel to 2 MeV multiplies low-energy gamma strength by 2.5","feed_subtitle":"New predictions give fusion-evaporation experiments and stellar rate models concrete benchmarks to test.","key_machinery":"The central object is the finite-temperature relativistic quasiparticle random phase approximation (FT-RQRPA), a thermally averaged extension of RPA built on finite-temperature Hartree-Bardeen-Cooper-Schrieffer occupation probabilities that include Fermi-Dirac factors, with the QRPA matrix diagonalized self-consistently using the DD-PCX point-coupling energy density functional. This machinery generates temperature-dependent electric dipole strength distributions, transition probabilities, and single-particle occupation numbers, and it is what lets the authors trace how thermal unblocking opens new low-energy excitation channels.","core_discovery":"Using the finite-temperature relativistic quasiparticle random phase approximation (FT-RQRPA) with the DD-PCX relativistic point-coupling interaction, the authors compute the isovector electric dipole response of nickel isotopes from mass 56 to 70 at temperatures from $T=0$ to 2 MeV. They find that in neutron-rich isotopes the integrated E1 strength in the $E=0$ to 12 MeV region increases by up to a factor of 2.5 at $T=2$ MeV compared with the zero-temperature case, and that in lighter isotopes such as nickel-56 a low-energy pygmy component appears only once the temperature rises. The calculations attribute this to thermal unblocking: occupation probabilities of orbitals above the Fermi level grow, and new quasiparticle configurations such as neutron $1g_{9/2}$ to $1h_{11/2}$ and proton $2p_{3/2}$ to $2d_{5/2}$ transitions contribute, while some zero-temperature transitions are suppressed. For nickel-62 at $T=1.6$ MeV, the predicted SEWS in the 8-12 MeV window is 4.32%, matching a preliminary experimental value of about 4% from fusion-evaporation studies.","pith_inferences":["If thermal unblocking is the dominant mechanism, similar enhancement should appear in other open-shell medium-mass nuclei with many single-particle levels near the Fermi surface, not just nickel.","A decisive test would be to extract gamma-ray strength from the same nickel isotope produced at several beam energies corresponding to different internal temperatures; the predicted monotonic growth of low-energy strength with temperature would confirm the mechanism.","If confirmed, the temperature dependence argues for including T-dependent E1 strength functions in statistical-model codes for stellar neutron-capture rates instead of scaling zero-temperature strength ad hoc.","Beyond 2 MeV the enhancement may saturate or be affected by width effects, so extending the calculation to higher temperatures and comparing with data would map the thermal limit."],"forward_implications":["The predicted SEWS and cumulative B(E1) tables for nickel-56 through nickel-70 give ongoing fusion-evaporation experiments direct quantitative benchmarks for hot pygmy dipole strength.","In neutron-rich nickel isotopes, low-energy E1 strength at $T=2$ MeV should be up to 2.5 times the zero-temperature value, a signal within reach of gamma-ray spectroscopy.","Isotopes near neutron-proton symmetry, which show no pygmy strength at zero temperature, should develop visible low-energy dipole strength at temperatures around 1.5 to 2 MeV.","A substantial part of the E1 strength moves from the high-energy region down to the 0-12 MeV region as temperature rises, so the hot pygmy strength grows partly at the expense of higher-energy response.","Temperature-dependent low-energy E1 strength of this kind can affect radiative neutron-capture rates used in models of heavy-element synthesis in stars."],"supporting_citations":[{"why":"Establishes the FT-RQRPA framework whose equations produce the temperature-dependent E1 spectra.","marker":"[50]"},{"why":"Companion paper developing the finite-temperature quasiparticle RPA formulation used here.","marker":"[51]"},{"why":"Supplies the DD-PCX interaction with realistic symmetry energy that governs the E1 transition strength.","marker":"[52]"},{"why":"Provides the temperature-dependent matrix elements and B(E1) transition expressions used in the calculation.","marker":"[46]"},{"why":"Documents thermal unblocking as the source of enhanced low-energy dipole strength at finite temperature.","marker":"[34]"},{"why":"Zero-temperature experimental SEWS for nickel-68 used to benchmark the T=0 limit.","marker":"[31]"},{"why":"Second nickel-68 measurement of low-lying E1 strength used as a benchmark.","marker":"[32]"},{"why":"Nickel-70 experimental SEWS values for neutron-rich benchmark comparison.","marker":"[33]"},{"why":"Reports the preliminary 62Ni hot-pygmy measurement, about 4% at T around 1.6 MeV, that the calculation matches.","marker":"[38]"}],"fun_headline_variants":["Hot nickel: 2.5x low-energy dipole strength at 2 MeV","Thermal unblocking: Ni isotopes gain pygmy dipole at high T","Nickel-62 matches experiment: SEWS 4.32% at 1.6 MeV","Pygmy dipole in nickel emerges as temperature rises","Hot pygmy dipole: predictions for Ni-56 to Ni-70"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central premise is that a thermally occupied mean field plus quasiparticle RPA, without beyond-RPA damping, anharmonicity, or pairing fluctuations, captures the temperature evolution of the low-energy E1 response.","fun_headline_variants_meta":{"raw":{"variants":["Hot nickel: 2.5x low-energy dipole strength at 2 MeV","Thermal unblocking: Ni isotopes gain pygmy dipole at high T","Nickel-62 matches experiment: SEWS 4.32% at 1.6 MeV","Pygmy dipole in nickel emerges as temperature rises","Hot pygmy dipole: predictions for Ni-56 to Ni-70"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000834,"raw_usage":{"total_tokens":3676,"prompt_tokens":1020,"completion_tokens":2656,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":636,"completion_tokens_details":{"reasoning_tokens":2552}},"tokens_in":636,"tokens_out":2656,"duration_ms":18432,"temperature":1.0,"reasoning_tokens":2552,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:04:37.097026+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the gamma-ray strength function of nickel-62 or nickel-66 from fusion-evaporation reactions at temperatures of 1.5 to 2 MeV: if the integrated E1 strength below 12 MeV does not rise toward 2.5 times the zero-temperature value, or if the 8-12 MeV SEWS in nickel-62 at $T=1.6$ MeV is far from 4.32%, the thermal-enhancement claim would be falsified.","supporting_citations":[{"cited_title":"Wibowo and E","cited_arxiv_id":null,"evidence_quote":"Establishes the FT-RQRPA framework whose equations produce the temperature-dependent E1 spectra."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Companion paper developing the finite-temperature quasiparticle RPA formulation used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the DD-PCX interaction with realistic symmetry energy that governs the E1 transition strength."},{"cited_title":"Y¨ uksel, G","cited_arxiv_id":null,"evidence_quote":"Provides the temperature-dependent matrix elements and B(E1) transition expressions used in the calculation."},{"cited_title":"Wieland, A","cited_arxiv_id":null,"evidence_quote":"Documents thermal unblocking as the source of enhanced low-energy dipole strength at finite temperature."},{"cited_title":"Bertulani, The European Physical Journal A 55, 240 (2019)","cited_arxiv_id":null,"evidence_quote":"Zero-temperature experimental SEWS for nickel-68 used to benchmark the T=0 limit."},{"cited_title":"Wieland, A","cited_arxiv_id":null,"evidence_quote":"Second nickel-68 measurement of low-lying E1 strength used as a benchmark."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Nickel-70 experimental SEWS values for neutron-rich benchmark comparison."},{"cited_title":"Wieland, A","cited_arxiv_id":null,"evidence_quote":"Reports the preliminary 62Ni hot-pygmy measurement, about 4% at T around 1.6 MeV, that the calculation matches."}],"review_version":2}