Shape transitions and ground-state properties of tungsten isotopes in covariant density functional theory
Pith reviewed 2026-05-10 17:30 UTC · model grok-4.3
The pith
Covariant density functional theory maps shape transitions across tungsten isotopes up to the neutron drip line at N=184
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The calculations reveal a dynamic shape evolution: spherical configurations dominate at the magic neutron numbers N=82 and N=126, prolate shapes prevail through the intermediate region, and shape coexistence occurs in isotopes such as 158W, 160W, 194W, 196W, 206W and near 244-248W. A potential subshell closure at N=118 is indicated by anomalies in two-neutron separation energies together with vanishing pairing energies. The neutron drip line is located at N=184 where the shape returns to spherical symmetry. The results remain consistent with available experimental data and with independent calculations that use Skyrme, finite-range droplet, and other relativistic mean-field approaches.
What carries the argument
Quadrupole deformation parameter beta_2 and self-consistent potential energy curves generated by relativistic mean-field calculations with pairing included, using the DD-ME1, DD-ME2, DD-PC1 and DD-PCX functionals.
If this is right
- Shape coexistence in the listed isotopes implies low-lying excited 0+ states whose energies and transition strengths can be measured directly.
- The vanishing pairing at N=118 produces a local minimum in separation energies that should appear in precise mass measurements.
- The return to sphericity at N=184 sets an upper limit on the isotopes that can participate in r-process nucleosynthesis paths.
- Strong agreement with Skyrme SLy4 and NL3 results indicates that the predicted deformation trends are not strongly functional-dependent.
Where Pith is reading between the lines
- The same pattern of shape evolution and the N=118 feature may appear in neighboring isotopic chains such as osmium or hafnium.
- If the drip-line prediction holds, it constrains the fission barriers and beta-decay rates needed for superheavy-element synthesis models.
- Extending the calculations to odd-A tungsten isotopes would test whether the single-particle spectrum remains consistent with the even-even results.
Load-bearing premise
The four relativistic density functionals remain accurate when extrapolated far beyond the mass region where they were originally fitted, especially near the neutron drip line.
What would settle it
An experimental measurement that places the neutron drip line for tungsten before N=184 or that finds no evidence of shape coexistence in 194W or 196W would contradict the calculated location and deformation patterns.
Figures
read the original abstract
This study investigates the structural evolution of even-even tungsten isotopes ($^{154\text{--}264}$W) using covariant density functional theory (CDFT) with four relativistic functionals: DD-ME1, DD-ME2, DD-PC1, and DD-PCX. Key nuclear properties, including binding energies, quadrupole deformation parameters, two-neutron separation energies, neutron pairing energies, nuclear radii, and potential energy curves, are analyzed to explore shape transitions and stability from neutron-deficient to neutron-rich isotopes up to the drip line. The results reveal a dynamic shape evolution, with spherical configurations at $N = 82$ and $N = 126$, prolate dominance in intermediate regions, and shape coexistence in isotopes such as $^{158}$W, $^{160}$W, $^{194}$W, $^{196}$W, $^{206}$W, and near $^{244\text{--}248}$W. A potential subshell closure at $N = 118$ is identified, supported by anomalies in separation energies and vanishing pairing energies. The neutron drip line is predicted at $N = 184$, marked by a return to spherical symmetry. Comparisons with experimental data and other theoretical models, including the deformed Hartree-Fock-Bogoliubov method with the Skyrme SLy4 interaction, the Finite Range Droplet Model, and the Relativistic Mean Field model with NL3, show strong agreement, validating the robustness of CDFT. These findings enhance our understanding of nuclear structure in the medium-to-heavy mass region and provide insights relevant to r-process nucleosynthesis, thereby guiding future experimental studies at radioactive ion beam facilities.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper applies covariant density functional theory with four relativistic functionals (DD-ME1, DD-ME2, DD-PC1, DD-PCX) to even-even tungsten isotopes ranging from ^{154}W to ^{264}W. It examines binding energies, quadrupole deformations, two-neutron separation energies, pairing energies, nuclear radii, and potential energy curves to map out shape transitions, identify a possible subshell closure at N=118, and locate the neutron drip line at N=184 where spherical symmetry returns. The results are compared to experimental data and other models such as Skyrme SLy4, FRDM, and NL3.
Significance. The study provides a comprehensive view of shape evolution across a wide isotopic chain in tungsten, highlighting regions of shape coexistence and potential magic numbers. The use of multiple functionals and cross-comparisons with other theoretical approaches adds credibility to the findings on shape transitions in the stable region. However, the predictions for the neutron-rich side, particularly the drip line, carry larger uncertainties due to the extrapolation involved. If the subshell closure and drip-line location are confirmed, they would be relevant for understanding r-process nucleosynthesis paths.
major comments (2)
- [Neutron drip line and shape near N=184] The central claim of the neutron drip line at N=184 with a return to spherical symmetry is load-bearing for the paper's conclusions on stability limits. However, this relies on the four functionals' performance at extreme neutron excess (N-Z ~ 100 for ^{264}W), where they were not originally fitted. The manuscript should provide a quantitative estimate of the uncertainty in the two-neutron separation energies or perform a sensitivity study by varying the isovector parameters to show how much the drip line can shift.
- [Identification of subshell closure at N=118] The anomalies in two-neutron separation energies and vanishing pairing energies at N=118 are used to suggest a subshell closure. To establish this as robust rather than functional-specific, the paper should show the single-particle energy spectra or the evolution of the Fermi surface for at least two functionals in that region, as small changes in pairing or deformation could mimic such anomalies.
minor comments (3)
- [Abstract] The abstract claims 'strong agreement' with experiment and other models but provides no quantitative measures such as root-mean-square deviations for binding energies or deformation parameters. Including such metrics would improve clarity.
- [Figures] The potential energy curves and deformation plots would benefit from indicating the energy scale and labeling the functionals clearly in each panel to facilitate direct comparison.
- [Methods] A brief discussion of the fitting procedure or the range of nuclei used to calibrate the four functionals would help readers assess the reliability of the extrapolation to the drip line.
Simulated Author's Rebuttal
We thank the referee for the thorough review and constructive feedback on our manuscript. We have carefully considered the major comments and provide point-by-point responses below. Where appropriate, we have revised the manuscript to address the concerns.
read point-by-point responses
-
Referee: [Neutron drip line and shape near N=184] The central claim of the neutron drip line at N=184 with a return to spherical symmetry is load-bearing for the paper's conclusions on stability limits. However, this relies on the four functionals' performance at extreme neutron excess (N-Z ~ 100 for ^{264}W), where they were not originally fitted. The manuscript should provide a quantitative estimate of the uncertainty in the two-neutron separation energies or perform a sensitivity study by varying the isovector parameters to show how much the drip line can shift.
Authors: We acknowledge that predictions at the neutron drip line involve extrapolation and thus uncertainties. The manuscript already employs four different covariant density functionals, which were fitted to different datasets, providing a built-in measure of model dependence. In the revised version, we have added a quantitative estimate of the uncertainty by reporting the standard deviation of the two-neutron separation energies across the four functionals near N=184, which is on the order of 0.5-1 MeV. This indicates that the location of the drip line at N=184 is stable within the model variations. A dedicated sensitivity analysis by varying isovector parameters would require a separate study, but the current multi-functional approach addresses the robustness to a significant extent. revision: partial
-
Referee: [Identification of subshell closure at N=118] The anomalies in two-neutron separation energies and vanishing pairing energies at N=118 are used to suggest a subshell closure. To establish this as robust rather than functional-specific, the paper should show the single-particle energy spectra or the evolution of the Fermi surface for at least two functionals in that region, as small changes in pairing or deformation could mimic such anomalies.
Authors: We agree that displaying the single-particle spectra would strengthen the claim. In the revised manuscript, we have included the single-particle neutron energy levels for the DD-ME2 and DD-PC1 functionals around N=118, demonstrating a clear gap at the Fermi surface consistent with the observed anomalies in separation energies and pairing. This supports that the subshell closure is a feature of the mean-field structure rather than an artifact of pairing or deformation. revision: yes
Circularity Check
No circularity: standard application of pre-existing functionals
full rationale
The paper applies four established covariant density functionals (DD-ME1, DD-ME2, DD-PC1, DD-PCX) to even-even tungsten isotopes from N=82 to the drip line. These functionals predate the study and were fitted on other nuclei; the tungsten calculations are direct numerical outputs of the standard CDFT equations. Reported features (shape evolution, N=118 subshell signature, N=184 drip line) emerge from the computed binding energies, deformations, separation energies and pairing gaps. External comparisons to experiment, Skyrme HFB, FRDM and NL3 RMF provide independent benchmarks. No self-definitional redefinitions, fitted-input predictions, or load-bearing self-citations appear in the derivation chain.
Axiom & Free-Parameter Ledger
free parameters (1)
- Parameters of DD-ME1, DD-ME2, DD-PC1, DD-PCX functionals
axioms (1)
- domain assumption Covariant density functional theory with the chosen functionals accurately describes ground-state properties and deformations across the tungsten isotopic chain.
Reference graph
Works this paper leans on
-
[1]
Meson-Exchange Model The meson-exchange model describes the nucleus as a system of Dirac nucleons interacting via finite-mass meson exchange [ 37, 38, 42, 57, 64–67]. The principal mesons mediating the interaction are the isoscalar-scalar σ meson, the isoscalar-vector ω meson, and the isovector- vector ρ meson. The corresponding Lagrangian density reads [...
-
[2]
Point-Coupling Model The point-coupling formulation replaces meson ex- change with contact interactions among nucleons [ 76–78]. The Lagrangian reads: L = ¯ψ(iγµ∂µ − m)ψ − 1 2 αS(ˆρ)( ¯ψψ)2 − 1 2 αV (ˆρ)( ¯ψγ µψ)2 − 1 2 αT V (ˆρ)( ¯ψ⃗ τ γµψ)2 − 1 2 δS(∂ν ¯ψψ)(∂ν ¯ψψ) − e ¯ψγ µAµ (1 − τ3) 2 ψ. (7) We adopt the DD-PC1 [ 58] and DD-PCX [ 59] param- eterizati...
-
[3]
= −Gδ(R − R′)P (r)P (r′), (11) where R = 1√ 2 (r1 + r2) and r = 1√ 2 (r1 − r2), with P (r) taken as a Gaussian: P (r) = 1 (4πa2)3/2 exp ( − r2 2a2 ) . (12) The pairing strength is taken as G = 728 MeV·fm3 with a range parameter a = 0 .644 fm, calibrated to re- produce the pairing gap of the Gogny D1S interaction in symmetric nuclear matter [ 84]. C. Compu...
- [4]
-
[5]
G G Kiss and Z Podolyák European Physical Journal A 60 175 (2024)
work page 2024
-
[6]
Physical Review C 72 047303 (2005) 13
P D Stevenson et al. Physical Review C 72 047303 (2005) 13
work page 2005
-
[7]
U Rahaman International Journal of Modern Physics E 34 2550043 (2025)
work page 2025
-
[8]
U Rahaman, M Ikram, I A Rather and A A Usmani Brazilian Journal of Physics 54 105 (2024)
work page 2024
-
[9]
Physical Review C 80 064308 (2009)
N Alkhomashi et al. Physical Review C 80 064308 (2009)
work page 2009
-
[10]
Physical Review C 63 011304 (2000)
C Wheldon et al. Physical Review C 63 011304 (2000)
work page 2000
-
[11]
Physical Review C 90 021301 (2014)
P R John et al. Physical Review C 90 021301 (2014)
work page 2014
-
[12]
European Physical Journal A 53 1 (2017)
A Rohilla et al. European Physical Journal A 53 1 (2017)
work page 2017
-
[13]
Acta Physica Polonica B 48 337 (2017)
S K Chamoli et al. Acta Physica Polonica B 48 337 (2017)
work page 2017
-
[14]
P Cejnar, J Jolie and R F Casten Reviews of Modern Physics 82 2155 (2010)
work page 2010
-
[15]
L M Robledo, R Rodríguez-Guzmán and P Sarriguren Journal of Physics G: Nuclear and Particle Physics 36 115104 (2009)
work page 2009
-
[16]
P Sarriguren, R Rodríguez-Guzmán and L M Robledo Physical Review C 77 064322 (2008)
work page 2008
-
[17]
K Nomura, T Otsuka, R Rodríguez-Guzmán, L M Rob- ledo, P Sarriguren, P H Regan, P D Stevenson and Z Podolyák Physical Review C 83 054303 (2011)
work page 2011
-
[18]
K Nomura, T Otsuka, R Rodríguez-Guzmán, L M Rob- ledo and P Sarriguren Physical Review C 84 054316 (2011)
work page 2011
-
[19]
J E García-Ramos, K Heyde, L M Robledo and R Rodríguez-Guzmán Physical Review C 89 034313 (2014)
work page 2014
-
[20]
Physical Review C 111 L031301 (2025)
K Sharma, P Chowdhury, A M Rogers et al. Physical Review C 111 L031301 (2025)
work page 2025
-
[21]
A Dadwal and H M Mittal European Physical Journal A 53 132 (2017)
work page 2017
-
[22]
A M Khalaf, M M Sirag and K E Abdelmageed Chinese Journal of Physics 54 329 (2016)
work page 2016
-
[23]
P Jain, Y Kumar, V Kumar and P Bangotra Universe 9 362 (2023)
work page 2023
-
[24]
N Sharma and H M Mittal International Journal of Modern Physics E 22 1350053 (2013)
work page 2013
-
[25]
Physical Review C 61 024607 (2000)
B Djerroud et al. Physical Review C 61 024607 (2000)
work page 2000
-
[26]
Physical Review C 79 031304 (2009)
J Barea, R Bijker, A Frank et al. Physical Review C 79 031304 (2009)
work page 2009
-
[27]
U Rahaman, M Ikram, M Imran and A A Usmani In- ternational Journal of Modern Physics E 29 2050076 (2020)
work page 2020
-
[28]
A A Usmani, S A Abbas, U Rahaman, M Ikram and F H Bhat International Journal of Modern Physics E 27 1850060 (2018)
work page 2018
-
[29]
H Morinaga Physical Review 101 254 (1956)
work page 1956
-
[30]
D Bonatsos, A Martinou, S K Peroulis, T J Mertzimekis and N Minkov Atoms 11 117 (2023)
work page 2023
-
[31]
R V F Janssens and T L Khoo Annual Review of Nuclear and Particle Science 41 321 (1991)
work page 1991
-
[32]
Physical Review C 102 034311 (2020)
K E Karakatsanis et al. Physical Review C 102 034311 (2020)
work page 2020
-
[33]
P D Stevenson, M P Brine, Zs Podolyák, P H Regan, P M Walker and J Rikovska Stone Physical Review C 72 047303 (2005)
work page 2005
-
[34]
Physical Review C 88 044301 (2013)
P J R Mason et al. Physical Review C 88 044301 (2013)
work page 2013
-
[35]
Physical Review C 82 051304 (2010)
G J Lane et al. Physical Review C 82 051304 (2010)
work page 2010
-
[36]
Proceedings of the National Academy of Sciences USA 118 e2020680118 (2021)
N Nakanishi et al. Proceedings of the National Academy of Sciences USA 118 e2020680118 (2021)
work page 2021
- [37]
-
[38]
Earth and Planetary Science Letters 651 119139 (2025)
D C Rubie et al. Earth and Planetary Science Letters 651 119139 (2025)
work page 2025
-
[39]
Proceedings of the National Academy of Sciences USA 113 2886 (2016)
G Budde et al. Proceedings of the National Academy of Sciences USA 113 2886 (2016)
work page 2016
-
[40]
T Nikšić, N Paar, D Vretenar and P Ring Computer Physics Communications 185 1808 (2014)
work page 2014
-
[41]
G A Lalazissis, T Nikšić, D Vretenar and P Ring Phys- ical Review C 71 024312 (2005)
work page 2005
-
[42]
X Roca-Maza, X Viñas, M Centelles, P Ring and P Schuck Physical Review C 84 054309 (2011)
work page 2011
-
[43]
T Nikšić, D Vretenar, P Finelli and P Ring Physical Review C 66 024306 (2002)
work page 2002
-
[44]
S E Agbemava, A V Afanasjev, D Ray and P Ring Phys- ical Review C 89 054320 (2014)
work page 2014
-
[45]
P Ring Progress in Particle and Nuclear Physics 37 193 (1996)
work page 1996
-
[46]
J Meng, H Toki, S G Zhou, S Q Zhang, W H Long and L S Geng Progress in Particle and Nuclear Physics 57 470 (2006)
work page 2006
-
[47]
H Liang, J Meng and S G Zhou Physics Reports 570 1 (2015)
work page 2015
-
[48]
S Shen, H Liang, W H Long, J Meng and P Ring Progress in Particle and Nuclear Physics 109 103713 (2019)
work page 2019
-
[49]
Frontiers of Physics 9 529 (2014)
Q S Zhang, Z M Niu, Z P Li et al. Frontiers of Physics 9 529 (2014)
work page 2014
-
[50]
P W Zhao and Z P Li International Journal of Modern Physics E 27 1830007 (2018)
work page 2018
-
[51]
Z X Ren and P W Zhao Physical Review C 102 021301 (2020)
work page 2020
-
[52]
10 (World Scientific, Singapore) (2016)
J Meng Relativistic Density Functional for Nuclear Structure International Review of Nuclear Physics Vol. 10 (World Scientific, Singapore) (2016)
work page 2016
-
[53]
H Abusara, A V Afanasjev and P Ring Physical Review C 85 024314 (2012)
work page 2012
-
[54]
S E Agbemava, A V Afanasjev, T Nakatsukasa and P Ring Physical Review C 92 054310 (2015)
work page 2015
-
[55]
S E Agbemava, A V Afanasjev, D Ray and P Ring Phys- ical Review C 95 054324 (2017)
work page 2017
-
[56]
J Meng and S G Zhou Journal of Physics G: Nuclear and Particle Physics 42 093101 (2015)
work page 2015
-
[57]
M Matev, A V Afanasjev, J Dobaczewski, G A Lalazissis and W Nazarewicz Physical Review C 76 034304 (2007)
work page 2007
-
[58]
A V Afanasjev Physical Review C 78 054303 (2008)
work page 2008
-
[59]
A V Afanasjev and S E Agbemava Physical Review C 93 054310 (2016)
work page 2016
-
[60]
T Nikšić, D Vretenar and P Ring Physical Review C 66 064302 (2002)
work page 2002
-
[61]
T Nikšić, D Vretenar and P Ring Physical Review C 78 034318 (2008)
work page 2008
-
[62]
E Yüksel, T Marketin and N Paar Physical Review C 99 034318 (2019)
work page 2019
-
[63]
J Dobaczewski, H Flocard and J Treiner Nuclear Physics A 422 103 (1984)
work page 1984
-
[64]
E Chabanat, P Bonche, P Haensel, J Meyer and R Scha- effer Nuclear Physics A 635 231 (1998)
work page 1998
-
[65]
J Erler, N Birge, M Kortelainen, W Nazarewicz, E Olsen, A M Perhac and M Stoitsov Nature 486 509 (2012)
work page 2012
-
[66]
M L Cescato and P Ring Physical Review C 57 134 (1998)
work page 1998
-
[67]
D Vretenar, A V Afanasjev, G A Lalazissis and P Ring Physics Reports 409 101 (2005)
work page 2005
-
[68]
S Typel and H H Wolter Nuclear Physics A 656 331 (1999)
work page 1999
-
[69]
R Brockmann and H Toki Physical Review Letters 68 3408 (1992) 14
work page 1992
-
[70]
F Hofmann, C M Keil and H Lenske Physical Review C 64 034314 (2001)
work page 2001
-
[71]
Y K Gambhir, P Ring and A Thimet Annals of Physics 198 132 (1990)
work page 1990
-
[72]
B D Serot and J D Walecka Advances in Nuclear Physics 16 1 (1986)
work page 1986
-
[73]
B D Serot and J D Walecka International Journal of Modern Physics E 6 515 (1997)
work page 1997
-
[74]
J D Walecka Annals of Physics 83 491 (1974)
work page 1974
-
[75]
J Boguta and A R Bodmer Nuclear Physics A 292 413 (1977)
work page 1977
-
[76]
P G Reinhard, M Rufa, J Maruhn, W Greiner and J Friedrich Zeitschrift für Physik A 323 13 (1986)
work page 1986
-
[77]
G A Lalazissis, J König and P Ring Physical Review C 55 540 (1997)
work page 1997
-
[78]
B G Todd-Rutel and J Piekarewicz Physical Review Let- ters 95 122501 (2005)
work page 2005
-
[79]
B A Nikolaus, T Hoch and D G Madland Physical Re- view C 46 1757 (1992)
work page 1992
-
[80]
J J Rusnak and R J Furnstahl Nuclear Physics A 627 495 (1997)
work page 1997
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.