Finite-range pairing in nuclear density functional theory
Pith reviewed 2026-05-22 12:05 UTC · model grok-4.3
The pith
Folding pair densities with Gaussians of radius near 1 fm stabilizes pairing calculations in nuclear DFT.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Folding the pair density with a Gaussian of finite range yields a pairing functional whose results become independent of the size of the single-particle basis once the folding radius is chosen near 1 fm, thereby removing the continuum divergences that plague zero-range functionals while leaving the quality of bound-state predictions essentially unchanged.
What carries the argument
Gaussian folding of the pair density, which converts a zero-range pairing interaction into a finite-range one that regularizes the pairing field against continuum divergences.
If this is right
- Pairing calculations become stable and basis-size independent even in very large spaces containing many continuum states.
- Gradient-dependent pairing functionals such as Fayans can be used reliably without strong cutoff dependence.
- Global surveys of nuclear properties gain robustness because results no longer hinge on arbitrary choices for marginally occupied states.
- The same regularization can be applied to other observables that involve pair densities without introducing new parameters.
Where Pith is reading between the lines
- The 1 fm scale may connect to the effective range of the underlying nucleon-nucleon force and could be tested by varying the folding function shape.
- The approach might extend naturally to time-dependent density functional calculations or to pairing in neutron-star crust matter where continuum effects are large.
- Comparison with ab initio pairing results in medium-mass nuclei could check whether the regularization preserves microscopic accuracy.
Load-bearing premise
The Gaussian folding at a radius near 1 fm keeps the essential short-range pairing physics intact while removing the divergences.
What would settle it
Numerical results with the folded functional that deviate markedly from converged zero-range results in small model spaces or from experimental pairing gaps would indicate the regularization has altered the physics.
Figures
read the original abstract
Pairing correlations are ubiquitous in low-energy states of atomic nuclei. To incorporate them within nuclear density functional theory, often used for global computations of nuclear properties, pairing functionals that generate nucleonic pair densities and pairing fields are introduced. Many pairing functionals currently used can be traced back to zero-range nucleon-nucleon interactions. Unfortunately, such functionals are plagued by deficiencies that become apparent in large model spaces that contain unbound single-particle (continuum) states. In particular, the underlying computational schemes diverge as the single-particle space increases, and the results depend on how marginally occupied states are incorporated. These problems become more pronounced for pairing functionals that contain gradient-density dependence, such as in the Fayans functional. To remedy this, finite-range pairing functionals are introduced. In this study, this is done by folding the pair density with Gaussians. We show that a folding radius of about 1\,fm offers the best compromise between quality and stability, and substantially reduces the pathological behavior in different numerical applications.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes finite-range pairing functionals for nuclear density functional theory by folding the pair density with a Gaussian to address divergences and numerical instabilities that arise with zero-range interactions in large model spaces containing continuum states. It focuses on mitigating problems that are exacerbated in gradient-dependent functionals such as the Fayans form and concludes that a folding radius of approximately 1 fm provides the optimal compromise between quality and numerical stability while substantially reducing pathological behaviors across applications.
Significance. If the numerical demonstrations confirm that bound-state observables remain essentially unchanged while stability improves, the approach would offer a low-overhead practical fix for continuum-related issues in global DFT calculations, adding only a single parameter without requiring a full microscopic finite-range force. This could enhance reliability of pairing predictions in extended bases and reduce sensitivity to cutoff choices.
major comments (2)
- [Abstract] Abstract: the assertion that a folding radius of about 1 fm 'offers the best compromise between quality and stability' and 'substantially reduces the pathological behavior' is presented without any quantitative tables, error estimates, or explicit comparisons to zero-range results, leaving the central claim without visible supporting data in the summary of findings.
- [Folding procedure and numerical tests] The folding construction (direct Gaussian smearing of the pair density rather than derivation from a finite-range force): because only the radius is added and no compensating readjustment of the pairing strength is introduced, any alteration of the effective pairing field or its coupling to density gradients must be shown not to shift gaps or odd-even staggering even in small, fully convergent bases; explicit verification of equivalence for bound-state observables is required to support the claim that essential physics is preserved.
minor comments (1)
- Clarify in the text how the Gaussian folding is implemented numerically (e.g., whether it is applied before or after discretization) and ensure all convergence plots explicitly label the model-space sizes and radii tested.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major comment below and indicate the revisions we will make.
read point-by-point responses
-
Referee: [Abstract] Abstract: the assertion that a folding radius of about 1 fm 'offers the best compromise between quality and stability' and 'substantially reduces the pathological behavior' is presented without any quantitative tables, error estimates, or explicit comparisons to zero-range results, leaving the central claim without visible supporting data in the summary of findings.
Authors: The abstract is a concise summary of conclusions supported by the detailed numerical tests and comparisons in the body of the manuscript, including systematic variations of folding radii, model-space sizes, and direct contrasts with zero-range results shown in the figures and tables of the results section. We agree that the abstract would benefit from a brief quantitative anchor. In the revised manuscript we will update the abstract to reference the observed reduction in pathological cutoff dependence while preserving the overall length and readability. revision: yes
-
Referee: [Folding procedure and numerical tests] The folding construction (direct Gaussian smearing of the pair density rather than derivation from a finite-range force): because only the radius is added and no compensating readjustment of the pairing strength is introduced, any alteration of the effective pairing field or its coupling to density gradients must be shown not to shift gaps or odd-even staggering even in small, fully convergent bases; explicit verification of equivalence for bound-state observables is required to support the claim that essential physics is preserved.
Authors: We acknowledge that the Gaussian folding is a direct regularization of the pair density rather than a microscopic finite-range force. In the manuscript we already demonstrate that the integrated pairing strength is preserved and that, for a folding radius of 1 fm, the changes to pairing gaps and odd-even staggering remain within numerical precision in fully convergent small bases. To make this verification more explicit we will add a short dedicated paragraph (or small table) in the revised version that directly compares bound-state observables between the original zero-range functional and the folded version in those small spaces. revision: yes
Circularity Check
No significant circularity in derivation chain
full rationale
The paper introduces finite-range pairing functionals by folding the pair density with a Gaussian of chosen radius ~1 fm as a direct modification to address continuum divergences in zero-range functionals. The radius is selected numerically for best compromise between quality and stability, with results demonstrated in applications. This construction is presented as an independent ansatz and numerical fix rather than a derivation that reduces by definition or self-citation to the same data or inputs it aims to predict. No load-bearing steps, self-definitional equations, or fitted inputs called predictions are evident that would make central claims equivalent to their inputs by construction. The approach remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- Gaussian folding radius
axioms (1)
- domain assumption Standard nuclear density functional theory framework remains valid when the pairing interaction is made finite-range via Gaussian folding.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
finite-range pairing functionals are introduced. In this study, this is done by folding the pair density with Gaussians. We show that a folding radius of about 1 fm offers the best compromise between quality and stability
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The pairing density ˘ρF(r) ... ˆFφ_α ... Gaussian folding: F(r) ∝ exp(−r²/˜R_F²)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
The local densities Nuclear energy-density functionals are usually ex- pressed in terms of a local densities and currents. As we consider here only ground states of even-even nuclei, we can limit ourselves to the time-even densities: local densityρ, kinetic-energy densityτ, spin-orbit densityJ, and pair density ˘ρ. They read: ρq(r) = X α∈q f(cut) α v2 α φ...
-
[2]
space The pairing occupation amplitudev α decreases slowly with s.p
Soft cutoff in s.p. space The pairing occupation amplitudev α decreases slowly with s.p. energyε α. Thus, the summations in (A1) do not converge if all s.p. states are considered. A possi- ble way to remedy this problem is to apply a soft cutoff weight depending on s.p. energyε α [12]: f(cut) α = 1 1 + exp [εα −(ϵ F −δϵ)]/(δϵ/10) ,(A3) whereδϵis the energ...
-
[3]
A density dependence in the Skyrme energy is included only in the volume term∝ρ 2 T
The Skyrme energy functional The Skyrme functional for time-even energy density can be written as [6] E= Z d3r{E kin +E Sk +E Sk,pair}+E Coul +E cm,(A4) 10 where Ekin = ℏ2 2m τ, ESk,T = 1X T=0 h C ρ T ρ2 T +C ∆ρ T ρT ∆ρT +C τ T ρT τT +C J T J2 T +C ∇J T ρT ∇·J T i , ESk,pair = 1 4 1− ρ ρpair X q∈{p,n} V0,q ˘ρ2 q , ECoul =e 2 1 2 ZZ d3r d3r′ ρp(r)ρp(r′) |r...
-
[4]
The Fayans energy functional Alternatively, we also discuss results obtained with the Fayans EDF [26]. The functional used in this work has a form where each of the three terms in the pairing func- tional is split into separate proton and neutron terms, each with its own strength coefficient. We abbreviate that as Fy(IVP3) to denote isospin variable pairi...
-
[5]
J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Micro- scopic theory of superconductivity, Phys. Rev.106, 162 (1957)
work page 1957
-
[6]
A. Bohr, B. R. Mottelson, and D. Pines, Possible anal- ogy between the excitation spectra of nuclei and those of the superconducting metallic state, Phys. Rev.110, 936 (1958)
work page 1958
-
[7]
P. Ring and P. Schuck,The Nuclear Many-Body Problem (Springer–Verl., New York, Heidelberg, Berlin, 1980)
work page 1980
-
[8]
D. J. Dean and M. Hjorth-Jensen, Pairing in nuclear sys- tems: from neutron stars to finite nuclei, Rev. Mod. Phys. 75, 607 (2003)
work page 2003
-
[9]
R. A. Broglia and V. Zelevinsky,Fifty Years of Nuclear BCS : Pairing in Finite Systems(World Scientific, 2013)
work page 2013
-
[10]
M. Bender, P.-H. Heenen, and P.-G. Reinhard, Self- consistent mean-field models for nuclear structure, Rev. Mod. Phys.75, 121 (2003)
work page 2003
-
[11]
Schunck, ed.,Energy Density Functional Methods for Atomic Nuclei, 2053-2563 (IOP Publishing, 2019)
N. Schunck, ed.,Energy Density Functional Methods for Atomic Nuclei, 2053-2563 (IOP Publishing, 2019)
work page 2053
-
[12]
E. Sapershte˘ ın and M. Troitski˘ ı, Sov. J. Nucl. Phys.1, 284 (1965)
work page 1965
-
[13]
Z. Bochnacki, I. Holban, and I. Mikhailov, Residual in- teraction in nuclei, Nucl. Phys. A97, 33 (1967)
work page 1967
-
[14]
R. R. Chasman, Density-dependent delta interactions and actinide pairing matrix elements, Phys. Rev. C14, 1935 (1976)
work page 1935
-
[15]
S. G. Kadmenski˘ ı, Y. L. Ratis, K. Rybak, and V. Furman, Choice of density-dependent effective interaction andα decay of heavy spherical nuclei, Sov. J. Nucl. Phys.27, 481 (1978)
work page 1978
-
[16]
S. J. Krieger, P. Bonche, H. Flocard, P. Quentin, and M. S. Weiss, An improved pairing interaction for mean– field calculations using skyrme potentials, Nucl. Phys. A 517, 275 (1990)
work page 1990
- [17]
-
[18]
S. Fayans and D. Zawischa, Towards a better parametrization of the nuclear pairing force: density dependence with gradient term, Phys. Lett. B383, 19 (1996)
work page 1996
-
[19]
E. Garrido, P. Sarriguren, E. Moya de Guerra, and P. Schuck, Effective density-dependent pairing forces in theT= 1 andT= 0 channels, Phys. Rev. C60, 064312 (1999)
work page 1999
-
[20]
J. Dobaczewski, W. Nazarewicz, and P.-G. Rein- hard, Pairing interaction and self-consistent densities in neutron-rich nuclei, Nucl. Phys. A693, 361 (2001)
work page 2001
- [21]
-
[22]
A. Bulgac and Y. Yu, Renormalization of the Hartree- Fock-Bogoliubov equations in the case of a zero range pairing interaction, Phys. Rev. Lett.88, 042504 (2002)
work page 2002
-
[23]
J. Dobaczewski, W. Nazarewicz, T. R. Werner, J. F. Berger, C. R. Chinn, and J. Decharg´ e, Mean-field de- scription of ground-state properties of drip-line nuclei: Pairing and continuum effects, Phys. Rev. C53, 2809 (1996)
work page 1996
-
[24]
T. Nikˇ si´ c, P. Ring, and D. Vretenar, Renormalized rela- tivistic Hartree-Bogoliubov equations with a zero-range pairing interaction, Phys. Rev. C71, 044320 (2005)
work page 2005
- [25]
-
[26]
P. J. Borycki, J. Dobaczewski, W. Nazarewicz, and M. V. Stoitsov, Pairing renormalization and regulariza- tion within the local density approximation, Phys. Rev. C73, 044319 (2006)
work page 2006
-
[27]
J. C. Pei, A. T. Kruppa, and W. Nazarewicz, Quasi- particle continuum and resonances in the hartree-fock- bogoliubov theory, Phys. Rev. C84, 024311 (2011)
work page 2011
-
[28]
J. C. Pei, N. Fei, Y. N. Zhang, and P. Schuck, Generalized second-order Thomas-Fermi method for superfluid Fermi systems , Phys. Rev. C92, 064316 (2015)
work page 2015
-
[29]
J. Dobaczewski and W. Nazarewicz, Hartree-Fock- Bogoliubov solution of the pairing Hamiltonian in finite nuclei, inFifty Years of Nuclear BCS, edited by R. A. Broglia and V. Zelevinsky (World Scientific, 2013) p. 40
work page 2013
-
[30]
P.-G. Reinhard and W. Nazarewicz, Toward a global de- scription of nuclear charge radii: Exploring the Fayans energy density functional, Phys. Rev. C95, 064328 (2017)
work page 2017
-
[31]
J. C. Pei, M. V. Stoitsov, G. I. Fann, W. Nazarewicz, 12 N. Schunck, and F. R. Xu, Deformed coordinate-space hartree-fock-bogoliubov approach to weakly bound nuclei and large deformations, Phys. Rev. C78, 064306 (2008)
work page 2008
-
[32]
M. Chen, T. Li, B. Schuetrumpf, P.-G. Reinhard, and W. Nazarewicz, Three-dimensional Skyrme Hartree- Fock-Bogoliubov solver in coordinate-space representa- tion, Comput. Phys. Commun.276, 108344 (2022)
work page 2022
- [33]
- [34]
- [35]
-
[36]
T. Duguet, Bare vs effective pairing forces: A microscopic finite-range interaction for Hartree-Fock-Bogolyubov cal- culations in coordinate space, Phys. Rev. C69, 054317 (2004)
work page 2004
-
[37]
Y. Tian, Z.-y. Ma, and P. Ring, Separable pairing force for relativistic quasiparticle random-phase approxima- tion, Phys. Rev. C79, 064301 (2009)
work page 2009
-
[38]
Y. Tian, Z.-y. Ma, and P. Ring, A finite range pairing force for density functional theory in superfluid nuclei, Phys. Lett. B676, 44–50 (2009)
work page 2009
-
[39]
Z.-y. Ma, Y. Tian, and P. Ring, Density functional theory with a separable pairing force in finite nuclei, Nucl. Phys. A834, 50c (2010)
work page 2010
-
[40]
P. Vesel´ y, J. Dobaczewski, N. Michel, and J. Toivanen, Finite-range separable pairing interaction within the new N3LO DFT approach, J. Phys.: Conf. Ser.267, 012027 (2011)
work page 2011
-
[41]
S. Teeti and A. V. Afanasjev, Global study of separable pairing interaction in covariant density functional theory, Phys. Rev. C103, 034310 (2021)
work page 2021
-
[42]
A. Migdal,Theory of Finite Fermi Systems and Appli- cations to Atomic Nuclei, Interscience monographs and texts in physics and astronomy (Interscience, New York, 1967)
work page 1967
- [43]
-
[44]
P. G. Reinhard, M. Bender, K. Rutz, and J. A. Maruhn, An HFB scheme in natural orbitals, Z. Phys. A358, 277 (1997)
work page 1997
-
[45]
N. Tajima, Canonical-basis solution of the Hartree-Fock- Bogoliubov equation on a three-dimensional Cartesian mesh, Phys. Rev. C69, 034305 (2004)
work page 2004
-
[46]
Y. N. Zhang, J. C. Pei, and F. R. Xu, Hartree-fock- bogoliubov descriptions of deformed weakly bound nu- clei in large coordinate spaces, Phys. Rev. C88, 054305 (2013)
work page 2013
-
[47]
T. Mitchell, J. Bollinger, X.-P. Huang, W. Itano, J. Tan, B. Jelenkovi´ c, and D. Wineland, Crystalline order in strongly coupled ion plasmas, inAIP Conf. Proc., Vol. 498 (1999) pp. 353–366
work page 1999
-
[48]
A. Piel, A. Homann, M. Klindworth, A. Melzer, C. Zafiu, V. Nosenko, and J. Goree, Waves and oscillations in plasma crystals, J. Phys. B36, 533 (2003)
work page 2003
-
[49]
P. Kl¨ upfel, P.-G. Reinhard, T. J. B¨ urvenich, and J. A. Maruhn, Variations on a theme by Skyrme: A systematic study of adjustments of model parameters, Phys. Rev. C 79, 034310 (2009)
work page 2009
-
[50]
A. J. Miller, K. Minamisono, A. Klose, D. Garand, C. Kujawa, J. D. Lantis, Y. Liu, B. Maaß, P. F. Mantica, W. Nazarewicz, W. N¨ ortersh¨ auser, S. V. Pineda, P. G. . Reinhard, D. M. Rossi, F. Sommer, C. Sumithrarachchi, A. Teigelh¨ ofer, and J. Watkins, Proton superfluidity and charge radii in proton-rich calcium isotopes, Nat. Phys. 15, 432 (2019)
work page 2019
-
[51]
P.-G. Reinhard, J. O’Neal, S. M. Wild, and W. Nazarewicz, Extended Fayans energy density functional: optimization and analysis, J. Phys. G51, 105101 (2024)
work page 2024
-
[52]
J. Karthein, C. M. Ricketts, R. F. Garcia Ruiz, J. Bil- lowes, C. L. Binnersley, T. E. Cocolios, J. Dobaczewski, G. J. Farooq-Smith, K. T. Flanagan, G. Georgiev, W. Gins, R. P. de Groote, F. P. Gustafsson, J. D. Holt, A. Kanellakopoulos, A. Koszor´ us, D. Leimbach, K. M. Lynch, T. Miyagi, W. Nazarewicz, G. Neyens, P.-G. Reinhard, B. K. Sahoo, A. R. Vernon...
work page 2024
-
[53]
J. Dobaczewski, H. Flocard, and J. Treiner, Hartree-fock- bogolyubov description of nuclei near the neutron-drip line, Nucl. Phys. A422, 103 (1984)
work page 1984
-
[54]
P.-G. Reinhard, B. Schuetrumpf, and J. Maruhn, The ax- ial hartree–fock + bcs code skyax, Comp. Phys. Comm. 258, 107603 (2021)
work page 2021
-
[55]
A. Bulgac and Y.-l. Yu, Local density approximation for pairing correlations in nuclei, (2001), arXiv:nucl- th/0109083
- [56]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.