{"id":"7329697b-d87d-4a10-ab81-e66a70db66ab","arxiv_id":"1908.05016","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"At the critical nuclear charge where an N-electron atom is on the verge of unbinding, the ground state decays faster than exp(-C sqrt(r)) in the K outermost electron distances, assuming Z_c < N-K.","lead":"This paper proves that the ground state of a many-electron atom at the critical nuclear charge, the moment it just barely binds, decays rapidly in the positions of the outermost electrons. The proof gives a new decay estimate at threshold without assuming an energy gap, extending the authors' companion method from helium to general atoms.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 2.3 omits the hypotheses used in its final positivity step: it never states Z=Z_c, K maximal, or Z_c<N-K, yet A1 requires W2=(N-K)/(1+δ)-Z>0 and A2 requires E_Z^(N-K-1)-E_Z^(N)>0.","rationale":"The reader's weakest assumption names exactly the same issue, so I agree. I cannot find a more severe flaw: the localization/IMS machinery in Section 3 is standard, Lemma 3.1 supplies the needed potential lower bounds, and the final positivity could be made to work once the missing hypotheses are supplied. The statement gap is nevertheless load-bearing because the theorem is the paper's central claim: it is the only result that justifies the abstract's decay statement, and as written it covers couplings for which W2 is negative and the proof's core inequality cannot be positive. The issue is not merely cosmetic wording; without Z=Z_c and K maximal, A2 may have no spectral gap to dominate the error terms. Because the concern is fixable by adding assumptions already present in Theorem 2.1/abstract, the verdict should remain CONDITIONAL rather than moving to REJECT. The alternative existence proof in Remark 2.2 is deferred to [7], so I did not count it as independent confirmation.","tokens_in":10453,"tokens_out":8340,"duration_ms":86441,"concrete_test":"Compute the decisive A1 term at a noncritical coupling. Take Z=N-K and δ=1/2; then W2=(N-K)/(3/2)-(N-K)=-(N-K)/3<0, so the leading term of A1 is negative for every choice of c_k and the proof's positivity step fails. If the intended domain is only the critical case, amend Theorem 2.3 to state explicitly: Z=Z_c, K is the largest integer with E_{Z_c}^{(N)}=E_{Z_c}^{(N-K)}, and Z_c<N-K, and verify that with these hypotheses W2>0 (choosing δ small) and E_{Z_c}^{(N-K-1)}-E_{Z_c}^{(N)}>0. This check determines whether the flaw is a missing hypothesis or a genuine counterexample to the statement.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing problem is that Theorem 2.3 is stated for an arbitrary normalized eigenfunction ψ_Z of H_Z^(N), with K appearing in G but never defined in the theorem. In the proof, K is silently taken from Theorem 2.1 and two quantitative conditions are used as if they were hypotheses. First, the positivity of A1 (p. 12) needs W2=(N-K)/(1+δ)-Z>0. This is exactly the abstract's assumption Z_c<N-K, but the theorem never says Z=Z_c, and if Z is larger, W2<0 and A1 is not positive. Second, positivity of A2 uses 'the assumption that E_Z^(N-K-1)-E_Z^(N)>0'. At Z=Z_c with K chosen as the largest integer with E_{Z_c}^{(N)}=E_{Z_c}^{(N-K)} this gap follows from maximality; for arbitrary Z it is not guaranteed and can fail. Since the whole method reduces to showing the operators A1,A2,A3 in (8) are positive, the conclusion e^Gψ∈L² is not established for the stated domain. The fix is to add the missing hypotheses to Theorem 2.3; the proof outline then appears consistent. The companion-paper deferral in Remark 2.2 is a separate limitation, but not the main correctness risk.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the decay of ground states of N-electron atoms at the critical nuclear charge Z_c, i.e., at the edge of the essential spectrum where the ground-state energy coincides with the ionization threshold. The main result (Theorem 2.3) asserts that if ψ_Z is a normalized ground state of the Hamiltonian H_Z^{(N)} at energy E_Z^{(N)}, then e^G ψ_Z is square-integrable, where G is a weighted sum of the square roots of the K largest electron distances in the region where the K outer electrons are well separated from the inner ones, and a weighted linear sum otherwise. The proof uses a weighted IMS localization method developed in the companion paper [7], exploiting the repulsive part of the interelectronic potential to obtain positivity of certain quadratic forms without assuming a spectral gap. Existence of the threshold ground state is imported from [3]; the paper also sketches an alternative existence proof via tightness. The theorem is stated for an arbitrary ground state of H_Z^{(N)} without specifying the critical value Z=Z_c or the definition of K.","tokens_in":10753,"tokens_out":18358,"duration_ms":167850,"significance":"If the stated result is correct, it gives a quantitative decay estimate for a bound state exactly at the threshold of a many-body Coulomb system, with no spectral gap and no dependence on the Born-Oppenheimer approximation (though the main theorem as written assumes an infinitely heavy nucleus; the finite-mass case is only addressed in Remark 2.5). This is a genuinely interesting and nontrivial question, and the method of using repulsive interactions to compensate for the absence of a gap is a useful complement to Agmon-type techniques. The estimates in Lemma 3.1 are plausible and supported by an iterative computation, and the argument is presented in a self-contained way apart from the deferred companion-paper details. A strength is that the constants in the theorem have an explicit form, and the argument is adaptable to bosons or distinguishable particles. However, the main theorem statement omits essential hypotheses used in the proof, which is a serious correctness issue that must be addressed.","major_comments":[{"comment":"Theorem 2.3 is not correctly stated. The integer K appears in the definition of G in Eq. (3) and throughout the proof, but it is never defined in the theorem. The proof silently takes K from Theorem 2.1—the largest integer with E_{Z_c}^{(N)}=E_{Z_c}^{(N-K)}—and also assumes Z=Z_c. The final positivity step on p. 12 additionally requires W2=(N-K)/(1+δ)-Z>0 and the strict gap E_Z^{(N-K-1)}-E_Z^{(N)}>0, neither of which appears among the theorem's hypotheses. As written, the theorem asserts a decay statement for an arbitrary ground state of H_Z^{(N)} for arbitrary Z and an arbitrary natural number K, for which the operators A1 and A2 in (8) need not be positive. The theorem should be restated with Z=Z_c, with K the specified maximal integer, and with the assumptions Z_c<N-K (or equivalently W2>0 for a suitably small δ) and E_{Z_c}^{(N-K-1)}>E_{Z_c}^{(N)} made explicit.","section":"Theorem 2.3 (Section 2), p. 5"},{"comment":"The positivity argument for A2 is incomplete as written. The lower bound for |∇F|^2 in Lemma 3.2 on the relevant region contains a constant term d_k^2 that does not vanish as R→∞. The sentence 'the assumption that E_Z^{(N-K-1)}-E_Z^{(N)}>0 and the fact that we take sufficiently large R' does not account for this constant negative contribution; only the terms W1/|x|_k and e_k/sqrt(|x|_k) are small for large R. The proof needs to state explicitly that the constants C_m and K_m in the definition of F (Eq. (6)) are chosen small enough so that sum_k d_k^2 < E_Z^{(N-K-1)}-E_Z^{(N)}, which is possible but is not mentioned. Additionally, the text at this point says 'with N-K > Z', but the condition needed for A1>0 is (N-K)/(1+δ) > Z (i.e., W2>0); the weaker condition N-K>Z is insufficient because of the factor 1/(1+δ) in the definition of W2.","section":"Proof of Theorem 2.3, final step (p. 12)"}],"minor_comments":[{"comment":"There are grammatical errors ('atoms undergoes', 'we derive upper bound for the bound state', 'Our method do not require') that should be corrected.","section":"Abstract"},{"comment":"The statement 'We make the standard assumptions that the nucleus is infinitely heavy and at the origin, i.e. Born-Oppenheimer approximation' contradicts the abstract's claim that the method does not require the Born-Oppenheimer approximation. The abstract should be qualified, since the main theorem as stated assumes an infinitely heavy nucleus; the finite-mass extension appears only in Remark 2.5 and with an additional mass condition.","section":"Section 2, first paragraph"},{"comment":"The constants C_m and K_m in Eq. (3) are left unspecified; the statement should mention that they are chosen as in the proof, or that they are sufficiently small positive constants.","section":"Theorem 2.3"},{"comment":"The first line 'Let U = -Z/|x_j| + sum_{j≠k} 1/|x_j-x_k|' should read 'U = sum_{j=1}^N -Z/|x_j| + sum_{j≠k} 1/|x_j-x_k|'.","section":"Lemma 3.1"},{"comment":"The computation '|∇|x|_m|^2 = |sum_j ∂_j |x|_m|^2 = |∂_k |x_k||^2 = 1' is written in nonstandard notation; it should say that on Ω, the gradient of |x|_m is the unit vector pointing in the direction of the corresponding coordinate, so its squared norm is 1.","section":"Proof of Lemma 3.2"},{"comment":"The letter K is used both for the number of outer electrons and as a generic constant in inequalities (e.g., '≤K' on pp. 10-12). This is confusing; a different symbol (e.g., C) for the constants would help.","section":"Throughout the proof"},{"comment":"The sentence 'it remains to take the limit ε→0 in (8)' should also address the role of R; after (8), a sentence explaining that for a fixed large R the tail is controlled and the inner region is harmless because e^F is bounded on {|x|_outer≤R} would be useful.","section":"End of proof of Theorem 2.3"},{"comment":"The reference to [7] as 'will appear on arXiv shortly' is not a stable reference; the authors should update it if available, or state the necessary facts from [7] explicitly when they are used.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The central claim is plausible, and the paper is worth publishing after the theorem is restated with the missing hypotheses. The proof relies on the companion paper [7] for some technical details, which may slow refereeing but is not fatal if [7] is available. The manuscript would also benefit from a careful proofreading pass. I would like to see the revised version before making a final decision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a genuinely new result — sqrt-exponential decay for N-electron atoms at critical charge — and the method is interesting. The main theorem as stated, however, is missing its own hypotheses. The proof outline is plausible and the flaw is fixable; this needs a serious revision, not a rejection.\n\nWhat is actually new: Theorem 2.3 gives, for an atom at Z = Z_c, decay of the ground state like exp(-C Σ sqrt(|x|_k)) over the K outermost electron coordinates, with no spectral gap and without Agmon-type assumptions. Earlier threshold existence results ([3], [5]) did not provide decay rates. The trick of using the repulsive inter-electron potential to gain positivity without a gap is real and worth taking seriously. The generalized IMS localization formula for piecewise C^1 cutoffs in the appendix is also a useful piece of work.\n\nWhere it is soft: the stress-test note is correct. Theorem 2.3 says 'let H_Z^(N) be given and let ψ_Z be a normalized eigenfunction,' defines G with K, and never states Z = Z_c, never says K is the maximal integer from Theorem 2.1, never states Z < N-K. In the proof, positivity of A1 uses W2 = (N-K)/(1+δ) - Z > 0, and positivity of A2 uses E_Z^(N-K-1) - E_Z^(N) > 0. Both enter as silent assumptions. Take the theorem literally and it is false: for Z > N-K, W2 is negative, and the conclusion is not established. The fix is simple — add the standing hypotheses from Theorem 2.1 to the statement of Theorem 2.3 — and the proof outline becomes consistent.\n\nSmaller issues: K is never defined in the theorem; key technical estimates are deferred to the companion paper [7], including the alternative existence proof; the abstract says the method needs no Born-Oppenheimer approximation, but Section 2 assumes an infinitely heavy nucleus and Remark 2.5 says finite mass requires extra work; there are typos (Lemma 3.1's U is missing a sum, and the proof of Theorem A.2 cites Theorem A.2 itself).\n\nNet: the core idea is sound, the result is new, and once the theorem is restated correctly it is a solid contribution. It deserves a referee. I'd want the referee to insist on the hypothesis fix and to flag the companion-paper dependency before publication. I would not cite it in its current form.","headline":"New threshold decay result for atoms that is probably right, but Theorem 2.3 is missing the hypotheses its own proof needs; fixable, worth refereeing.","tokens_in":11298,"tokens_out":4717,"would_cite":false,"duration_ms":43243,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81Q10","81V45"],"pacs":[],"model":"deepseek-v4-flash","headline":"At the critical nuclear charge $Z_c$, every ground state of an $N$-electron atom decays at least as fast as $\\exp(-C\\sum\\sqrt{|x_k|})$ over the $K$ outermost electrons, with no spectral gap or Born-Oppenheimer approximation.","keywords":["threshold bound states","critical nuclear charge","exponential decay of eigenfunctions","Coulomb many-body systems","Schrödinger operators","no spectral gap","Born-Oppenheimer approximation","tightness"],"falsifier":"For the two-electron atom at its critical charge $Z_c\\approx 0.91$ (a case where a threshold ground state is known to exist), compute the ground state numerically and test whether $\\int_{|x_1|,|x_2|>R}|e^{C(\\sqrt{|x_1|}+\\sqrt{|x_2|})}\\psi(x)|^2\\,dx$ stays bounded as $R\\to\\infty$; if it diverges for every $C>0$, Theorem 2.3 would be false for that system.","tokens_in":10239,"feed_emoji":"⚛️","tokens_out":16114,"duration_ms":138166,"temperature":0.7,"pith_summary":"An atom with $N$ electrons stops binding when the nuclear charge is lowered to a critical value $Z_c$; at exactly that coupling the ground-state energy lies at the edge of the essential spectrum. This paper claims that the ground state still decays at least as fast as $\\exp(-C\\sum_{k=1}^{K}\\sqrt{|x_k|})$, where the sum runs over the $K$ largest electron–nucleus distances and $K$ is the number of electrons whose removal leaves the ground-state energy unchanged, provided $Z_c<N-K$. The proof needs no spectral gap and no Born-Oppenheimer approximation; the repulsive electron–electron interaction is what supplies the missing positivity. If the claim is right, threshold atoms decay super-polynomially in their outermost coordinates, and the same estimate gives a new, tightness-based route to proving that such ground states exist.","feed_headline":"Critical atoms decay like exp(-C√r), no spectral gap needed","feed_subtitle":"New bound on the K outermost electrons needs no spectral gap and no Born-Oppenheimer approximation.","key_machinery":"The engine of the proof is a weighted $L^2$ estimate built from a piecewise-$C^1$ version of the IMS localization formula. Configuration space is split according to how many inner electrons lie within a $\\delta$-fraction of the outermost distance, and an upper bound $F=\\sum_{m=N+1-K}^{N}(C_m\\sqrt{|x|_m}+K_m|x|_m\\chi^{\\perp}_{0,2\\delta})$ is introduced for the weight $G$. After subtracting localization errors of the form $|\\nabla F|^2$, the task reduces to proving positivity of three quadratic forms $A_1,A_2,A_3$. Positivity is driven by $Z_c<N-K$, which makes $W_2=(N-K)/(1+\\delta)-Z$ positive, so the outer electrons feel net repulsion large enough to dominate $|\\nabla F|^2$; the square-root decay comes from $|\\nabla \\sqrt{|x|}|^2\\approx 1/(4|x|)$ while the repulsive term is also of order $1/|x|$.","core_discovery":"At $Z=Z_c$, the paper's central result (Theorem 2.3) gives an $L^2$ decay bound: for the Hamiltonian $H_Z^{(N)}=\\sum_{j=1}^{N}(-\\Delta_j-Z/|x_j|)+\\sum_{j\\neq k}1/|x_j-x_k|$ on the antisymmetric subspace $L^2_a(\\mathbb{R}^{3N})$, every normalized eigenfunction $\\psi_Z$ with $H_Z^{(N)}\\psi_Z=E_Z^{(N)}\\psi_Z$ satisfies $e^G\\psi_Z\\in L^2_a(\\mathbb{R}^{3N})$. The weight $G$ is a sum of square roots, $G=\\sum_{m=N+1-K}^{N}C_m\\sqrt{|x|_m}$, whenever the $K$ outermost electrons are well separated from the $N-K$ inner ones, and a linear sum of the same outer distances otherwise. Thus, although the eigenvalue sits exactly at the threshold of the essential spectrum and no gap is available, the wavefunction decays no slower than $\\exp(-C\\sqrt{r})$ in the outer coordinates. The proof localizes the eigenvalue equation with cutoff functions and uses the repulsive Coulomb terms to make the relevant quadratic forms positive; no Born-Oppenheimer approximation is used.","pith_inferences":["One natural extension is to ask whether the same repulsion-driven sqrt-decay appears at thresholds set by cluster breakup rather than single-electron loss; the positivity trick may carry over to Coulomb systems with several nuclei.","The strict inequality $Z_c<N-K$ suggests a transition as $Z_c$ approaches $N-K$ from below: the constants in the decay exponent should deteriorate, possibly marking the change from a threshold bound state to no bound state. A quantitative version of that blow-up is not in the paper.","The $L^2$ decay statement could likely be upgraded to pointwise bounds via standard elliptic estimates once $e^G\\psi\\in L^2$ is known; the paper does not attempt this.","Since the method avoids the Born-Oppenheimer approximation, it may extend to molecular Hamiltonians with finite nuclear masses and give decay rates at autoionization thresholds, where the relevant outer coordinates are electron distances from the center of mass."],"forward_implications":["At the critical charge, the ground-state wavefunction is square-integrable against $\\exp(2C\\sum\\sqrt{|x|_m})$, so the electron density in the $K$ outermost channels decays super-polynomially in the outer coordinates.","Because the decay estimate is uniform for every $Z>Z_c$, it combines with tightness criteria to give an alternative proof that a threshold ground state exists.","The argument does not use the fermionic sign, so the same decay bound applies to bosons and to distinguishable particles.","For a finite-mass nucleus whose mass is at least the total electron mass, the same estimate holds after straightforward modifications.","When the outer and inner electrons are not well separated, the theorem yields only linear exponential decay in the outer coordinates, reflecting stronger screening in that configuration."],"supporting_citations":[{"why":"Supplies the existence of the threshold ground state whose decay the paper then proves.","marker":"[3]"},{"why":"The companion paper whose repulsion-driven weighted-estimate method is adapted here.","marker":"[7]"},{"why":"Provides the IMS localization formula used to convert the eigenvalue equation into quadratic-form inequalities.","marker":"[4]"},{"why":"The tightness criteria used to turn the uniform decay estimate into an alternative existence proof.","marker":"[8]"}],"fun_headline_variants":["Critical atoms decay like exp(-C√r) without a spectral gap","No spectral gap, no Born-Oppenheimer: √r decay at critical charge","At the brink: bound states decay exponentially in √r","Threshold eigenstates: exponential decay in outer coordinates"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof relies on $Z_c<N-K$, with $K$ the largest number of electrons whose removal leaves the ground-state energy unchanged, so that the repulsive Coulomb terms dominate the localization error; if that inequality fails, the final positivity step that produces the decay bound is not available.","fun_headline_variants_meta":{"raw":{"variants":["Critical atoms decay like exp(-C√r) without a spectral gap","No spectral gap, no Born-Oppenheimer: √r decay at critical charge","At the brink: bound states decay exponentially in √r","Threshold eigenstates: exponential decay in outer coordinates"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000207,"raw_usage":{"total_tokens":1407,"prompt_tokens":962,"completion_tokens":445,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":578,"completion_tokens_details":{"reasoning_tokens":372}},"tokens_in":578,"tokens_out":445,"duration_ms":5320,"temperature":1.0,"reasoning_tokens":372,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:26:50.428088+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For the two-electron atom at its critical charge $Z_c\\approx 0.91$ (a case where a threshold ground state is known to exist), compute the ground state numerically and test whether $\\int_{|x_1|,|x_2|>R}|e^{C(\\sqrt{|x_1|}+\\sqrt{|x_2|})}\\psi(x)|^2\\,dx$ stays bounded as $R\\to\\infty$; if it diverges for every $C>0$, Theorem 2.3 would be false for that system.","supporting_citations":[{"cited_title":"Frank, Elliott H","cited_arxiv_id":null,"evidence_quote":"Supplies the existence of the threshold ground state whose decay the paper then proves."},{"cited_title":"existence and decay rates of bound states at thres holds; helium , will appear on arXiv shortly (2019), 25","cited_arxiv_id":null,"evidence_quote":"The companion paper whose repulsion-driven weighted-estimate method is adapted here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the IMS localization formula used to convert the eigenvalue equation into quadratic-form inequalities."},{"cited_title":"Nonlinear Sci","cited_arxiv_id":null,"evidence_quote":"The tightness criteria used to turn the uniform decay estimate into an alternative existence proof."}],"review_version":1}