{"id":"87042189-1fca-4237-b182-21cb6e3f6782","arxiv_id":"2505.05455","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":2.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"A variational calculation with a one-parameter effective charge reproduces lithium and beryllium ground state energies to 2.5% and 3.1% of experiment.","lead":"The paper calculates the ground state energies of lithium and beryllium atoms using a trial wavefunction with a single adjustable effective charge. It is a textbook-level variational calculation, and the energies it finds are within about 3% of experiment.","discovery_kind":"incremental","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the central variational energies are correct upper bounds; the single-λ restriction is a limitation, and the printed prefactor/integral typos do not alter the final values.","rationale":"The reader's conditional verdict is appropriate, but my stress-test does not find a new load-bearing objection. I verified the one- and two-body matrix elements independently: the kinetic and nuclear averages follow from hydrogenic expectations with the actual nuclear charge, and the electron-electron terms match Slater-Condon rules when the exchange integral is 16/729. The final fractions in Eqs. (21) and (32) are algebraically consistent with Eqs. (19) and (30). Thus the central claim—that these are valid variational upper bounds agreeing with experiment to 2.5–3.1%—is supported. The single-λ assumption is the weakest physical modeling choice, but because the variational principle only requires an upper bound, it does not create a correctness risk. The real defects are typographical: the wrong I3 denominator and inconsistent prefactors in the displayed derivation. These are serious for reproducibility but not for the numerical claim. Therefore no verdict change is warranted; the paper should still be conditional on cleaning up the presentation.","tokens_in":7028,"tokens_out":23652,"duration_ms":255876,"concrete_test":"Recompute V_ee for the lithium determinant with the Slater-Condon rules: V_ee = I1 + 2I2 - I3 = 5/8 + 34/81 - 16/729 = 5965/5832 in units of λα²mc². If this coefficient reproduces Eq. (19)'s -33401/5832 after adding the nuclear term -27/4, the final energy claim is sound and the printed errors are confined to typos; if not, the numerical claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Independent re-derivation confirms the central numbers. For Li, the Slater determinant gives T=9/8 λ², V_nuc=-27/4 λ, and V_ee=J(1s,1s)+2J(1s,2s)-K(1s,2s)=5/8+34/81-16/729=5965/5832 λ, which leads exactly to λ0=33401/13122 and -198.345 eV. For Be, T=5/4 λ², V_nuc=-10λ, and V_ee=I1+I4+4I2-2I3=586373/373248λ, giving -386.663 eV. The one-parameter ansatz is a modeling simplification; since the result is a rigorous variational upper bound, this is not a correctness flaw. The load-bearing weakness is therefore not physical but presentational: Eq. (18)/A8 print I3=16/272 instead of 16/729, and Eqs. (15)/(28) display prefactors ('3×', '6/4!') that are inconsistent with the bracket expressions yet are not used in the final coefficients. These typos make the derivation non-reproducible as printed but do not invalidate the stated energies.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents variational calculations of the ground-state energies of lithium and beryllium. The trial wavefunctions are Slater determinants built from hydrogenic 1s and 2s orbitals that share a single effective charge λ, chosen to minimize the energy expectation. Using the virial theorem for the one-electron terms and direct/exchange Coulomb integrals for the electron-electron repulsion, the authors obtain -198.345 eV for Li and -386.663 eV for Be, compared with experimental values of -203.48 eV and -399.14 eV (relative errors 2.5% and 3.1%). The manuscript concludes that the variational method, together with the antisymmetry requirement and atomic shell structure, gives good agreement with experiment.","tokens_in":7263,"tokens_out":16732,"duration_ms":145164,"significance":"If the printed derivation is corrected, the central results are valid variational upper bounds and the agreement with experiment is a genuine, falsifiable test of the one-parameter screened-hydrogenic model. The paper's main strength is that the final energies are independently reproducible from the stated integrals (once the misprints in Eqs. (15), (18), and (28) are fixed), and the calculation is self-contained rather than fitted to the target energies. The physical content is, however, elementary: the single-λ ansatz is a known approximation, and the paper does not go beyond standard textbook material. Its value for a research journal is therefore primarily pedagogical.","major_comments":[{"comment":"The exchange integral I3 is printed as 16/272 × λ/a_B, but the correct value is 16/729 × λ/a_B. This value is load-bearing: the coefficient 5965/5832 in Eq. (15) and Eq. (19) equals 5/8 + 2×(17/81) − 16/729, so the printed 16/272 is not the value used. As printed, a reader who follows Eq. (18) obtains a different electron-electron energy and hence a different minimized energy. The denominator and the intermediate algebra in Eq. (A8) must be corrected.","section":"II, Eq. (18) and Appendix A, Eq. (A8)"},{"comment":"The prefactors are internally inconsistent. In Eq. (15), the bracket expression integrates to the total electron-electron energy J_1s1s + 2J_1s2s − K_1s2s, so the leading factor '3×' would triple-count the interaction; the final coefficient 5965/5832 corresponds to the bracket without the '3×'. In Eq. (28), the factor 6αℏc×(1/4!) times the bracket equals the correct total J_1s1s + J_2s2s + 4J_1s2s − 2K, so the initial '3×' is again spurious. These prefactors do not enter the final coefficients, but their presence makes the derivation non-reproducible as printed.","section":"II, Eq. (15) and III, Eq. (28)"},{"comment":"The trial wavefunction assumes that the 1s and 2s electrons share a single effective charge λ. This imposes equal screening on core and valence electrons and prevents the wavefunction from describing different effective charges for the two shells. The paper should explicitly state this restriction and discuss its likely effect on the variational energy (allowing independent λ_1s and λ_2s would lower the upper bound). The absence of this discussion is a gap, though not a correctness error, because the variational principle guarantees the computed value is an upper bound.","section":"II, Eq. (6) and III, Eq. (22)"}],"minor_comments":[{"comment":"The line '(4312)^2 + 4312(−1342 − 4213 + 1243)' appears twice, and the permutation 4213 is missing from its expected place in the list; this should be corrected to give the complete 24-term expansion.","section":"Appendix B"},{"comment":"The angle-bracket notation in Eqs. (15) and (28) is undefined; please clarify whether ⟨e²/(4πε0 r12)⟩ denotes the expectation in Ψ and whether the subsequent integrals are single-pair or total pair contributions.","section":"II, Eq. (15) and III, Eq. (28)"},{"comment":"The experimental values −203.48 eV and −399.14 eV are cited without a reference; please provide a source for these values.","section":"II and III"},{"comment":"The reference list consists almost entirely of Chinese-language teaching journals; adding a standard textbook or review reference for the variational treatment of many-electron atoms would help readers assess the novelty and context.","section":"References"},{"comment":"The notation a'_B = 2a_B in Eq. (6) is introduced without explanation; a brief comment that this is the scaled Bohr radius for the hydrogenic 2s orbital with effective charge λ would improve readability.","section":"II, Eq. (6)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a straightforward, textbook-level calculation. Its central numbers are correct variational upper bounds, and the only substantive issues are misprints in the printed derivation. The manuscript might be more suitable for a pedagogical journal than for a research journal in atomic physics; however, if the journal publishes such teaching-oriented contributions, it could be acceptable after a thorough correction of the equations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short answer: this is a standard variational calculation for lithium and beryllium, and the central numbers are right. I re-derived the coefficients: for Li, the Slater determinant gives V_ee = 5965/5832 λ, λ0 = 33401/13122, and -198.345 eV; for Be, V_ee = 586373/373248 λ, λ0 = 3146107/933120, and -386.663 eV. Those agree with experiment to about 2.5% and 3.1%, and they are rigorous upper bounds because the trial wavefunction is normalized and the energy is minimized variationally.\n\nWhat the paper does well: it is a self-contained, explicit treatment of antisymmetry and shell structure, with the six-dimensional integrals worked out in an appendix. It does not oversell the method, and the comparison to experiment is honest. It is not new — refs. [11-17] already cover variational calculations for these atoms — but as a pedagogical derivation it is clear and reproducible in principle.\n\nSoft spots, in decreasing order of importance. First, Eq. (18) and Eq. (A8) print the exchange integral I3 as 16/272; the correct value, used later in Eq. (15), is 16/729. As printed, the derivation cannot be followed, even though the final energy is unaffected. Second, Eqs. (15) and (28) have prefactor displays (\"3×\", \"6/4!\") that do not match the bracketed integrands; again, the final coefficients are fine, but the printed steps are sloppy. Third, the trial wavefunction uses a single effective charge λ for both 1s and 2s electrons. That is a real modeling restriction — the 1s and 2s shells would generally screen differently — and it biases the answer. It is not a correctness flaw, because the variational principle still gives a valid upper bound, but the choice is motivated by simplicity rather than physics.\n\nThe citation pattern is unremarkable and there is no self-citation. The paper makes no inflated claims. For a teaching-oriented journal, this is a suitable manuscript after the typos are corrected. I would not cite it in a research paper, but a referee can verify the calculation quickly and the errors are mechanical, not conceptual. Send it to peer review if the venue values worked pedagogical examples.","headline":"Routine but correct variational calculation for Li and Be; the final energies check out, and the printed exchange-integral typo is the main thing to fix.","tokens_in":7769,"tokens_out":3381,"would_cite":false,"duration_ms":36832,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that a Slater-determinant trial wavefunction with a single effective charge reproduces the ground-state energies of lithium and beryllium to within 2.5% and 3.1% of experiment.","keywords":["variational method","ground state energy","lithium atom","beryllium atom","antisymmetric wavefunction","Slater determinant","atomic shell structure","effective charge"],"falsifier":"Repeat the minimization with two independent variational parameters, $\\lambda_{1s}$ for the 1s electrons and $\\lambda_{2s}$ for the 2s electrons. If the energy falls substantially below $-198.345$ eV for lithium and $-386.663$ eV for beryllium, the shared-$\\lambda$ assumption is the limiting simplification; if it barely moves, the one-parameter ansatz is exonerated.","tokens_in":6812,"feed_emoji":"⚛️","tokens_out":10530,"duration_ms":101957,"temperature":0.7,"pith_summary":"The paper sets out to show that a deliberately small variational calculation, forced to respect the indistinguishability of electrons and the atomic shell structure, reproduces the ground-state energies of lithium and beryllium to within a few percent. The trial wavefunction is a Slater determinant built from hydrogenic 1s and 2s orbitals with one shared screening charge $\\lambda$, and minimizing the energy over $\\lambda$ gives $-198.345$ eV for lithium and $-386.663$ eV for beryllium, against experimental values of $-203.48$ eV and $-399.14$ eV. The significance of the claim is that the two qualitative facts about many-electron atoms, antisymmetry and shell structure, do most of the work even when every other detail is kept intentionally simple.","feed_headline":"One screening charge puts Li and Be energies within 3.1 percent.","feed_subtitle":"The variational ground-state energies come out at -198.345 eV for lithium and -386.663 eV for beryllium.","key_machinery":"The central object is a one-parameter antisymmetrized trial wavefunction: a Slater determinant of hydrogenic 1s and 2s spin-orbitals, $\\psi_{100} \\propto e^{-\\lambda r/a_B}$ and $\\psi_{200} \\propto (1 - \\lambda r/(2a_B)) e^{-\\lambda r/(2a_B)}$, all sharing the same effective charge $\\lambda$. The argument is carried by two tools: the virial theorem fixes the kinetic-energy and electron-nucleus expectation values directly in terms of $\\lambda$, and the electron-electron repulsion is reduced to six-dimensional integrals evaluated by expanding $1/r_{12}$ in Legendre polynomials and integrating over angles. Because the resulting energy function is the quadratic $E(\\lambda) = a\\lambda^2 - b\\lambda$, the minimization is closed form: $\\lambda_0 = b/(2a)$ and $E_{\\min} = -b^2/(4a)$.","core_discovery":"On the paper's own terms, the central discovery is that the antisymmetrized trial wavefunction $\\Psi = (1/\\sqrt{N!})\\sum_P \\delta_P P[\\varphi_a(1)\\varphi_b(2)\\cdots]$, with $\\varphi_a = \\psi_{100}\\chi_+$, $\\varphi_b = \\psi_{100}\\chi_-$, $\\varphi_c = \\psi_{200}\\chi_+$, and (for beryllium) $\\varphi_d = \\psi_{200}\\chi_-$, leads to a closed-form energy function $E(\\lambda) = a\\lambda^2 - b\\lambda$ in units of $\\alpha^2 mc^2$. For lithium, $a = 9/8$ and $b = 33401/5832$, giving $\\lambda_0 = 33401/13122 \\approx 2.5454$ and $E \\approx -198.345$ eV; for beryllium, $a = 5/4$ and $b = 3146107/373248$, giving $\\lambda_0 = 3146107/933120 \\approx 3.3716$ and $E \\approx -386.663$ eV. The antisymmetry contributes an explicit exchange integral to the electron-electron interaction, and the shell structure decides which orbitals enter the determinant; the paper reports the resulting errors as 2.5% and 3.1% relative to experiment.","pith_inferences":["A testable extension would replace the shared $\\lambda$ with independent variational charges $\\lambda_{1s}$ and $\\lambda_{2s}$; if the energy minimum drops significantly below the quoted values, the single-parameter form is the main source of the residual error.","Applying the same construction along the lithium and beryllium isoelectronic sequences would produce predicted curves $\\lambda_0(Z)$ and error trends that the paper does not report, offering a cheap check of whether the ansatz remains adequate as the nuclear charge grows.","Because the Hamiltonian in Eq. (4) contains only nonrelativistic Coulomb terms while the comparison values are experimental, the reported agreement may involve a partial cancellation between missing correlation energy and missing relativistic corrections; disentangling the two is beyond the paper."],"forward_implications":["For any atom whose ground configuration is two filled $s$ shells, the same determinant construction gives a closed-form quadratic energy function, so the calculation remains fully analytic.","By the variational theorem, the quoted energies are upper bounds to the exact nonrelativistic energies, so part of the apparent closeness to experiment reflects a cancellation between neglected correlation and neglected relativistic or QED shifts.","The optimal effective charges, $\\lambda_0 \\approx 2.55$ for lithium and $3.37$ for beryllium, quantify how strongly the outer electrons screen the nuclear charge.","The method needs only a handful of Coulomb integrals and no numerical diagonalization, so it can be repeated by hand in a graduate quantum mechanics course and extended to lithium-like and beryllium-like ions by changing $Z$."],"supporting_citations":[{"why":"States the variational theorem used throughout: the energy expectation in any normalized trial state is an upper bound on the ground-state energy.","marker":"[9]"},{"why":"Textbook variational treatment of the helium atom that supplies the pattern this paper extends to lithium and beryllium.","marker":"[10]"},{"why":"Provides the Legendre expansion of the inverse electron-electron distance used in the appendix to evaluate the Coulomb integrals.","marker":"[20]"}],"fun_headline_variants":["Antisymmetric wavefunction predicts Li and Be ground states to 3%","Variational trick: Li and Be ground states within 3.1 percent","Screening charge plus exchange nails Li and Be ground-state energies","Closed-form solution puts Li and Be within 3 percent of experiment"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that every electron, in the 1s shell and in the 2s shell, feels the same effective nuclear charge, so one number $\\lambda$ controls all screening; if the inner and outer electrons screen the nucleus differently, the trial wavefunction is too rigid to capture that difference.","fun_headline_variants_meta":{"raw":{"variants":["Antisymmetric wavefunction predicts Li and Be ground states to 3%","Variational trick: Li and Be ground states within 3.1 percent","Screening charge plus exchange nails Li and Be ground-state energies","Closed-form solution puts Li and Be within 3 percent of experiment"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001065,"raw_usage":{"total_tokens":4426,"prompt_tokens":865,"completion_tokens":3561,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":481,"completion_tokens_details":{"reasoning_tokens":3483}},"tokens_in":481,"tokens_out":3561,"duration_ms":26338,"temperature":1.0,"reasoning_tokens":3483,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:03:56.195212+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the minimization with two independent variational parameters, $\\lambda_{1s}$ for the 1s electrons and $\\lambda_{2s}$ for the 2s electrons. If the energy falls substantially below $-198.345$ eV for lithium and $-386.663$ eV for beryllium, the shared-$\\lambda$ assumption is the limiting simplification; if it barely moves, the one-parameter ansatz is exonerated.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"States the variational theorem used throughout: the energy expectation in any normalized trial state is an upper bound on the ground-state energy."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Textbook variational treatment of the helium atom that supplies the pattern this paper extends to lithium and beryllium."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Legendre expansion of the inverse electron-electron distance used in the appendix to evaluate the Coulomb integrals."}],"review_version":1}