{"id":"86a9923b-2adc-49d0-b4d0-7f92469a8e2e","arxiv_id":"2505.09540","paper_version":3,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Valence orbitals are largely insensitive to the form of soft confinement, a new exponential potential localizes numerical atomic orbitals faster than polynomial potentials, and FHI-aims singular-potential defaults yield basis-set truncation errors ranging from about 1 microhartree to 10…","lead":"This paper compares four families of soft confinement potentials used to build numerical atomic orbital basis sets, and proposes a new exponential confinement potential that localizes orbitals faster than existing polynomial forms. It maps how orbital shapes and basis-set truncation errors depend on confinement parameters across the periodic table, which matters for the accuracy and cost of DFT codes such as FHI-aims, SIESTA, and GPAW.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (21) is not a basis-set truncation error: it measures confinement-induced deformation in an exact FEM calculation, so the BSTE analysis and the fixed-truncation-radius proposal in Figs. 16-18 are unsupported.","rationale":"The reader's weakest assumption concerned spherical averaging and the hard-wall FEM boundary. Those are not the most load-bearing weak points: Mg and Ca are closed-shell atoms, so spherical averaging is exact for them, and the hard-wall truncation is introduced with an explicit 1 μEh convergence criterion. The central qualitative claims about orbital insensitivity and faster exponential-potential localization are well supported by the variational FEM results and the systematic tables and figures. The genuine soft spot is the BSTE analysis in Section IVC3. Eq. (21) subtracts the confinement energy from the exact confined total energy and then subtracts the unconfined total energy, yielding the deformation energy of the confined density in the unconfined Hamiltonian. This quantity does not involve a finite NAO basis, so calling it a basis-set truncation error is a category error. The paper's proposed method for selecting r_i by fixing this quantity to a target value may still be a useful heuristic, but the claim that it controls NAO truncation error is not established. The correct test is to perform an actual unconfined atomic calculation in a finite NAO basis generated from the confined orbitals and compare the resulting energy error with Eq. (21). Because the BSTE claim is one of the abstract's three advertised contributions, the verdict should be conditional on correcting or reinterpreting this analysis, while the other central claims can remain accepted.","tokens_in":33661,"tokens_out":28182,"duration_ms":294226,"concrete_test":"For a representative set of atoms (e.g., Mg, Ne, Ar), construct a finite NAO basis from the radial functions generated by Eq. (12) at selected r_i values, then recompute the unconfined atom variationally in that basis. Define the true basis-set truncation error as E_unconfined^NAO − E_unconfined and compare it with ΔE(ri) from Eq. (21) for the same r_i. If the true error is not systematically related to Eq. (21), then Eq. (21) is not a basis-set truncation error and Figs. 16-18 should be reinterpreted as deformation-energy measures.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's final contribution, 'we assess NAO basis set truncation errors' (abstract; Section IVC3), rests on Eq. (21): ΔE(ri) = E_confined(ri) − E_confinement(ri) − E_unconfined. With E_confined the exact FEM total energy under the singular confinement and E_confinement the expectation value of V_c, this expression equals ⟨ψ_c|H0|ψ_c⟩ − E_unconfined, i.e., the energy of the unconfined atom evaluated at the frozen confined density minus the exact unconfined energy. This is a confinement-induced deformation energy, not a basis-set truncation error. No finite NAO basis is constructed or used in Eq. (21); the FEM calculation is converged to the complete-basis limit. Consequently, Figs. 16-18 measure how strongly confinement deforms each atom, and the proposal to choose r_i by 'fixing the BSTE' does not control any truncation error of a generated NAO basis. The qualitative insensitivity and exponential-potential localization claims are unaffected, but the basis-set truncation-error interpretation and the associated computational-savings conclusions are not supported.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reviews soft confinement potentials used in numerical atomic orbital (NAO) generation and studies four families of such potentials for the Mg and Ca atoms with fully numerical DFT calculations in the HelFEM finite-element program. The potentials considered are the finite-barrier potential of Eq. (6), the polynomial potential of Eq. (9), a newly introduced exponential potential of Eq. (14), and the singular potentials of Eqs. (10)-(12) with n = 1, 2, 3. The authors report that the valence orbital shapes are qualitatively insensitive to the form of the confinement potential, that all soft potentials approach the hard-wall limit smoothly and systematically, and that the new exponential potential yields smaller truncation radii than the polynomial potential at the same parameters. The paper also presents an analysis in Section IVC3 that is described as an assessment of NAO basis-set truncation errors (BSTEs) for the singular potentials, and it proposes fixing the confinement parameters by setting this quantity to a target value for atoms H-Xe.","tokens_in":33996,"tokens_out":6941,"duration_ms":72700,"significance":"The review and systematic numerical comparison of confinement potentials are valuable for the NAO community, and the proposed exponential confinement potential of Eq. (14) is a useful addition that appears to provide genuinely improved localization in the tested cases. The calculations are variational, converged to a stated 1 microhartree criterion, systematically scanned over parameters, and performed with the open-source HelFEM code, which strengthens confidence in the qualitative findings. The orbital-contraction and hard-wall-limit results for Mg and Ca are well supported by the figures and tables. However, the final BSTE analysis in Section IVC3 rests on a quantity that is not a basis-set truncation error, so the associated conclusions about computational savings are not supported as written.","major_comments":[{"comment":"The quantity defined in Eq. (21) is not a basis-set truncation error. Since the FEM calculations are converged to the complete-basis limit (Section III), no finite NAO basis is constructed or truncated, and Eq. (21) reduces to <psi_c|H_0|psi_c> - E_unconfined, i.e., the energy of the unconfined atom evaluated at the confined ground-state density. This is a confinement-induced deformation energy, not an error caused by truncating an NAO basis. Therefore the values labelled BSTE in Figs. 16-18 and the proposal to fix this quantity when choosing r_i do not control a basis-set truncation error of generated NAOs, and the r_c^6 computational-savings conclusion in Section IVC3 is unsupported. Please rename this quantity and revise the abstract, Section IVC3, and Section V accordingly, or alternatively carry out an actual calculation with a finite NAO basis to quantify the true truncation error.","section":"Section IVC3, Eq. (21)"},{"comment":"Because the plotted quantities are deformation energies rather than basis-set truncation errors, the proposed fixed-BSTE parameter-selection scheme is not a validated method for controlling NAO basis quality. The periodic trends in r_i shown in Figs. 17 and 18 reflect how strongly each atom responds to confinement, not how incomplete a generated NAO basis would be. The concluding statement in Section V that considerable computational savings may be achievable from this scheme therefore needs to be re-evaluated once the quantity is correctly interpreted.","section":"Section IVC3, Figs. 17-18; Section V"}],"minor_comments":[{"comment":"The radial Schrodinger equation as written is inconsistent with Eq. (4): the standard form should read [-1/2(d^2/dr^2 + (2/r)d/dr) + l(l+1)/(2r^2) + V(r)]R(r) = E R(r). The signs of the centrifugal and potential terms and the missing factor 1/2 do not affect the subsequent continuity argument about R''(r), but they should be corrected.","section":"Section III, Eq. (19)"},{"comment":"The phrase 'the the finite element method' contains a duplicated definite article and should read 'the finite element method'.","section":"Abstract and Section III"},{"comment":"The word 'adressed' is a typo and should be 'addressed'.","section":"Section V"},{"comment":"The table caption says 'the 3s orbital of the Mg atom in the finite-barrier potential at r = 4 a0'; this should be 'with r0 = 4 a0' to avoid confusion with a radial coordinate value.","section":"Table II"}],"recommendation":"major_revision","confidential_remarks":"The paper contains a solid and useful review plus well-executed numerical experiments for Mg and Ca, and the exponential-potential result is likely to be of interest to the NAO community. The central problem is that the BSTE section does not measure what its name claims; this can be fixed by reframing the analysis as confinement-induced deformation and removing the basis-set-truncation and computational-savings claims, or by adding genuine finite-basis NAO calculations. I do not see grounds for rejection, but the final contribution needs substantial revision before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a useful, mostly sound methods paper that deserves referee time, but the last contribution is mislabeled. The review of confinement potentials is thorough and the new exponential potential in Eq. (14) is genuine. The FEM calculations are variational, converged to 1 µEh, and the qualitative claim that orbital contraction is insensitive to the potential form is well supported for Mg and Ca. The hard-wall approach analysis and Tables III–IV are solid. The asymptotic appendix is nice.\n\nThe soft spot is Section IVC3. The stress-test concern lands. Eq. (21) does not measure a basis-set truncation error: E_confined − E_confinement is the energy of the unconfined Hamiltonian at the confined density, so the difference from E_unconfined is a confinement-induced deformation energy in an exact FEM calculation. No finite NAO basis is constructed or used. That means Figs. 16–18 and the 'fixed BSTE' proposal are not what they claim to be. They are maps of how much each atom's density is deformed by the singular potential at a given ri, and choosing ri to make that deformation small is a plausible heuristic, but it is not controlling the truncation error of a generated NAO basis. The computational-savings conclusions need to be rewritten accordingly.\n\nMinor: the calculations use spherically averaged densities, which is standard in NAO studies but could matter for open-shell anisotropic cases; the paper acknowledges this. Some of the fixed-radius periodic table data would be more convincing with a second atom or a few open-shell examples, but that is a refinement.\n\nWho is it for: people building or using NAO basis sets (FHI-aims, SIESTA, GPAW) and anyone doing confined-atom benchmarks. The review half alone has value. I would cite the exponential potential and the asymptotic analysis, not the BSTE numbers.\n\nRecommendation: send to peer review. The right referee will ask for a careful rewrite of Section IVC3 and its abstract claim about BSTEs, but the core material is sound and the new functional forms are worth publishing.","headline":"A useful and mostly sound systematic comparison of NAO confinement potentials with a genuinely new exponential potential, but the final BSTE analysis mislabels a confinement deformation energy as a basis-set truncation error.","tokens_in":34503,"tokens_out":3099,"would_cite":true,"duration_ms":33258,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The valence orbitals of Mg and Ca are nearly independent of which soft confinement potential is used to generate them.","keywords":["numerical atomic orbitals","confinement potentials","density functional theory","finite element method","exponential confinement potential","basis-set truncation errors","hard-wall limit","Mg and Ca atoms"],"falsifier":"Run the same confinement comparison for an open-shell atom such as O or Fe with unrestricted, non-spherically averaged DFT (or unrestricted Hartree-Fock) and compare the 2p or 3d radial orbitals under polynomial, exponential, and finite-barrier confinement at matched truncation error; if the orbital shapes differ by more than a small tolerance, the central insensitivity claim fails.","tokens_in":33501,"feed_emoji":"⚛️","tokens_out":5585,"duration_ms":52898,"temperature":0.7,"pith_summary":"This paper asks whether the specific soft potential used to confine an atom matters for the orbitals that come out. Working on Mg and Ca, it finds that the valence 3s and 4s orbitals are nearly the same for finite-barrier, polynomial, exponential, and singular confinement potentials, as long as the potential switches on outside the core. The paper also introduces a new exponential confinement potential that makes orbitals decay faster than the standard polynomial potential with the same exponent, and shows that any soft potential approaches the hard-wall limit smoothly when made steeper. The practical payoff is that numerical atomic orbital (NAO) basis sets can be built with whatever smooth potential is most convenient for integral evaluation, then cut off with a hard-wall boundary at negligible energy cost.","feed_headline":"Orbital shapes barely change across confinement potentials","feed_subtitle":"Mg and Ca study says NAO builders can pick any smooth wall, plus a new exponential option.","key_machinery":"The machinery is a set of radial confinement potentials appended to the atomic Kohn-Sham Hamiltonian, solved with a high-order finite element method in which the radial grid ends at a practical infinity $r_\\infty$ that acts as a hard wall. The new exponential potential, $V_c(r) = N![\\exp(r/r_0) - \\sum_{k=0}^{N-1}(r/r_0)^k/k!]$, has the same small-$r$ Taylor behavior as the polynomial potential $(r/r_0)^N$ but grows exponentially at large $r$, which is what gives faster decay. Hermite interpolating polynomials let the calculation force the orbital derivative to zero at $r_\\infty$, making the truncation check parameter-free. The key comparison tool is the norm distance between a soft-confined orbital and the best-fit hard-wall orbital, which quantifies how smoothly each potential approaches the hard-wall limit.","core_discovery":"The central claim is that ground-state atomic orbital shapes under soft confinement are largely independent of the confinement potential's functional form. For the Mg 3s and Ca 4s orbitals, the radial functions obtained with the four potential families coincide to the eye once the confinement onset radius is large enough that the core is untouched; differences appear only for strong, small-cavity confinement. The paper demonstrates this by direct fully numerical Kohn-Sham calculations, by showing that all potentials approach the hard-wall orbital systematically as their steepness grows, and by fitting hard-wall orbitals to soft-confined orbitals with small L2 norms. As a novel alternative, the exponential potential of eq. (14) is shown to localize the 3s orbital at smaller truncation radii than the polynomial potential with the same $N$, improving basis sparsity; and fixing basis-set truncation errors to a target value yields periodic-table-wide confinement radii that fluctuate periodically with atomic number $Z$.","pith_inferences":["If the spherically averaged density premise holds, the insensitivity result should extend to other closed-shell atoms and to open-shell atoms treated with fractional occupations; a natural next test is unrestricted non-spherical calculations.","The exponential potential's faster decay could make it useful outside NAO generation, for example in models of quantum dots or pressure-confined atoms where exponential localization is physically motivated.","The fixed-truncation-error calibration scheme suggests a path to automated, functional-by-functional NAO basis generation across the periodic table, reducing the manual parameter tuning currently needed in codes that use singular potentials.","The failed convergence of the $n=4$ singular potential, despite its promising asymptotic faster decay, points to a numerical-stabilization problem worth solving, since tighter localization would directly cut polyatomic integral costs."],"forward_implications":["NAO generators can choose confinement potentials for smoothness, integral ease, or strict support without changing the atom's essential valence shape.","The new exponential potential yields smaller truncation radii than polynomial confinement at equal $N$, so basis sets made from it should give sparser Hamiltonian matrices and cheaper polyatomic calculations.","Adding a hard-wall cutoff at a radius where the soft-confined orbital is already negligible (1 $\\mu E_h$ energy error) gives strictly localized, smooth NAO basis functions and removes numerical noise.","Fixing the basis-set truncation error to a target ($10^{-2}$, $10^{-3}$, or $10^{-4}$ $E_h$) produces consistent confinement radii across the periodic table, with periodic fluctuations that suggest where computational savings are largest.","For singular potentials with denominator exponent $n=1,2,3$, the exponent $n=3$ gives the smallest radii and $n=2$ the largest, though all differences stay below about 0.5 $\\AA$."],"supporting_citations":[{"why":"The authors' previous hard-wall confinement study, whose methodology and baseline the present work extends.","marker":"[81]"},{"why":"Introduced the singular confinement potential of eq. (10), one of the four potential families studied here.","marker":"[15]"},{"why":"Introduced the $n=2$ singular potential of eq. (11) and the default parameters used in the basis-set truncation error analysis.","marker":"[27]"},{"why":"Proposed combining soft confinement with a hard-wall boundary, the localization strategy the paper validates and systematizes.","marker":"[33]"},{"why":"Supplies the Hermite interpolating polynomial basis that allows the derivative to be forced to zero at the truncation radius.","marker":"[93]"},{"why":"Describes the finite element method implementation used for all fully numerical calculations in the paper.","marker":"[85]"}],"fun_headline_variants":["Orbitals ignore smooth-wall shape, study finds","Smooth walls barely bend Mg and Ca orbitals","Exponential wall cuts NAO radius, same orbital","Confinement shape moot for ground-state orbitals","Hard-wall limit reached with any steep potential"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"All results rest on DFT calculations with spherically averaged electron densities; if that averaging distorts valence orbital shapes for open-shell or nonspherical atoms, the insensitivity claim and the reported truncation radii may not carry over to real NAO generation.","fun_headline_variants_meta":{"raw":{"variants":["Orbitals ignore smooth-wall shape, study finds","Smooth walls barely bend Mg and Ca orbitals","Exponential wall cuts NAO radius, same orbital","Confinement shape moot for ground-state orbitals","Hard-wall limit reached with any steep potential"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000134,"raw_usage":{"total_tokens":1183,"prompt_tokens":1032,"completion_tokens":151,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":648,"completion_tokens_details":{"reasoning_tokens":79}},"tokens_in":648,"tokens_out":151,"duration_ms":2381,"temperature":1.0,"reasoning_tokens":79,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:30:23.116190+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same confinement comparison for an open-shell atom such as O or Fe with unrestricted, non-spherically averaged DFT (or unrestricted Hartree-Fock) and compare the 2p or 3d radial orbitals under polynomial, exponential, and finite-barrier confinement at matched truncation error; if the orbital shapes differ by more than a small tolerance, the central insensitivity claim fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Proposed combining soft confinement with a hard-wall boundary, the localization strategy the paper validates and systematizes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Hermite interpolating polynomial basis that allows the derivative to be forced to zero at the truncation radius."},{"cited_title":"Lewis , author K","cited_arxiv_id":null,"evidence_quote":"Describes the finite element method implementation used for all fully numerical calculations in the paper."}],"review_version":1}