{"id":"e987d206-0415-43ae-b79a-905b4508956c","arxiv_id":"1908.06175","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"A single-Gaussian variational calculation suggests the Born-Oppenheimer approximation is reliable for H2 but not for molecules where muons replace electrons.","lead":"The paper models hydrogen molecules where one or both electrons are replaced by muons, using one Gaussian wave function per particle. It suggests that the standard fixed-nucleus approximation works for ordinary H2 but not for muon-containing versions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Upper-bound variational energies cannot establish that mono-muonated H2 is unbound, and the text itself later predicts weak resonance binding; the BO-validity-only-for-H2 claim rests on that unsupported no-binding result.","rationale":"The reader's weakest-assumption analysis correctly identifies the central logical gap: a variational upper bound cannot prove the absence of binding. My reading of the full text reinforces this and adds a decisive internal contradiction: the paper explicitly says a resonance will 'very likely' lead to weak binding in the same system for which it earlier claimed 'there is no binding energy.' Since the BO-validity-only-for-H2 suggestion in Section 5 is tied to the mono-muonated system's lack of a chemical bond, this unsupported and self-contradictory no-binding result is the most load-bearing weakness. I do not see a reason to move beyond the reader's CONDITIONAL verdict, because the central suggestion could still survive if a converged non-BO calculation supports it; but the present evidence does not establish it.","tokens_in":9140,"tokens_out":13060,"duration_ms":150413,"concrete_test":"Repeat the mono-muonated four-body calculation with an explicitly correlated Gaussian basis of at least 1000 terms, after exact center-of-mass projection, and compute the converged energy as a function of the interatomic separation. Compare it with the converged H + (p mu) dissociation threshold obtained in the same basis. If the finite-separation minimum lies below the threshold, the no-binding claim is refuted; if it remains above with a quantified convergence error, the claim would then have actual variational support. The same calculation should also report the non-BO vs clamped-nucleus energy difference for all three molecules to test the BO-validity-only-for-H2 statement directly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4's no-binding conclusion for mono-muonated di-hydrogen is not licensed by the Rayleigh-Ritz method. The optimized trial energy at finite c lying above the c = infinity energy only gives an upper bound on the exact ground-state energy; it cannot exclude an exact state below the dissociation threshold. The ansatz is far too crude for this to be a reliable null result: the separated-atom limit is built from single Gaussians per particle, so the variational threshold is not converged and the basis error is large on the scale of the tiny binding margin being probed. The authors themselves state that 'significant quantitative meaning cannot be attributed to' the mono-muonated curves for c < 0.06, the region where the resonance is discussed. Most directly, the text contradicts itself: after the center-of-mass correction it says 'It seems very likely that the resonance effect will lead to a weak binding between the non-muonic hydrogen atom and the muonic hydrogen atom', whereas Section 4 asserts 'there is no binding energy.' The Section 5 suggestion that BO is valid only for di-hydrogen depends on this unsupported no-bond result, and the orbital-radii argument provides no quantitative estimate of the BO error. The central claim therefore remains unverified by the calculation as presented.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a four-fermion model of two protons and two negatively charged leptons, with masses varying to represent H2, mono-muonated dihydrogen (one electron replaced by a muon), and di-muonated dihydrogen. The trial wave function is a symmetrized product of four Gaussian orbitals centered at ±c/2, one for each particle; the Rayleigh quotients for the kinetic and Coulomb terms are evaluated analytically and minimized variationally over the orbital exponents and the center separation c. The authors report binding for H2 and di-muonated H2 and no binding for mono-muonated H2, add a center-of-mass correction, and use the resulting orbital radii to suggest that the Born-Oppenheimer approximation is valid only for ordinary H2 among the three systems considered.","tokens_in":9350,"tokens_out":4565,"duration_ms":48701,"significance":"If its conclusions were established, the paper would provide a compact, fully analytic variational illustration of the Born-Oppenheimer approximation and of the effect of lepton mass on molecular binding. The derivation of the Rayleigh quotients is clean and transparent, and the explicit treatment of the center-of-mass correction is a useful pedagogical feature. The paper also deserves credit for stating its own limitation that the curves have no significant quantitative meaning for small c. However, the main physical claims are not supported by the calculation as presented: the no-binding conclusion for mono-muonated H2 rests on a variational upper-bound argument, and the claim that the Born-Oppenheimer approximation is valid only for di-hydrogen is asserted from qualitative orbital radii rather than demonstrated. These issues affect the central message rather than only the presentation.","major_comments":[{"comment":"The conclusion that mono-muonated dihydrogen has 'no binding energy' is not licensed by the calculation. Rayleigh-Ritz variational energies are upper bounds to the exact ground-state energy, so finding E_trial(c) > E_trial(∞) does not exclude an exact state below the dissociation threshold. In fact, the quoted difference is only about 0.0051 a.u. (between -49.6078 at c=0 and -49.6129 at c=∞), while the separated-atom threshold itself is built from single Gaussians per particle; the variational error in that threshold is plausibly of the same order as the claimed non-binding margin. The paper therefore proves at most that the chosen trial ansatz does not bind the mono-muonated system in this model, not that the physical molecule is unbound.","section":"Section 4, Figure 3"},{"comment":"The text directly contradicts itself: after discussing the center-of-mass-corrected curves, it states that 'It seems very likely that the resonance effect will lead to a weak binding between the non-muonic hydrogen atom and the muonic hydrogen atom,' whereas Section 4 asserts that no chemical bond exists because there is no binding energy. The paper also warns that 'significant quantitative meaning cannot be attributed to the curves in Figures 1, 2 and 3' for c < 0.06, which is precisely the region where the resonance and possible weak binding are discussed. This internal inconsistency must be resolved, and any conclusion about binding or non-binding must be stated as a property of the trial model rather than of the physical system.","section":"Center of mass correction subsection"},{"comment":"The claim that the Born-Oppenheimer approximation 'is valid only in the case of the unsubstituted di-hydrogen molecule' is not established. The argument uses only the qualitative radii of the particle orbitals relative to their centers, without any quantitative measure of the error introduced by the BO approximation. A defensible claim would require comparing BO and non-BO energies for all three systems, or at least estimating the omitted kinetic-energy contribution of the heavy particles. Because the no-binding result for mono-muonated H2 is itself unsupported (see the first comment), the BO-validity conclusion that leans on it inherits the same fragility.","section":"Section 5"}],"minor_comments":[{"comment":"There are several typographical errors: 'Fot c = ∞' in the Figure 3 caption, 'inf ty' in the Figure 5 caption, 'aniversary' in Section 5, and 'Borm-Openheimer' in references [9]–[11]. These should be corrected.","section":"Throughout"},{"comment":"The operator S is described as 'the symmetrizing operator,' but it only exchanges particles within the lepton pair and within the proton pair rather than fully symmetrizing all four particles. The intended permutation symmetry should be stated explicitly to avoid confusion.","section":"Section 2"},{"comment":"The integral I_eP_V_EE and related expressions contain the factor (αβ' − γβ) in the denominator and an error function of the same quantity; in the H2 limit where α=γ and β=β', this combination vanishes and the limiting form should be supplied.","section":"Appendix"},{"comment":"The bullet list gives binding energies and orbital radii without error bars or convergence information; since the trial basis is a single Gaussian per particle, at least a brief statement about expected basis-set error would help the reader interpret the quantitative claims.","section":"Section 4, first paragraph"}],"recommendation":"major_revision","confidential_remarks":"The analytic machinery is sound and the paper is short, but the variational-principle error is fundamental to the main no-binding claim and the BO-validity claim is much stronger than the evidence. A revision that reframes all conclusions as statements about the trial model, reconciles the 'no binding' and 'weak binding' passages, and softens the Section 5 claim could make the paper publishable. The authors may also wish to compare with existing muonic-molecule calculations to calibrate the basis error; as it stands, the central physical message is not verified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a compact, readable variational calculation with analytic Gaussian integrals, but the central suggestion—that BO is valid only for ordinary H2 because mono-muonated dihydrogen has no chemical bond—is not established by the evidence presented. The Rayleigh–Ritz bound cuts against that conclusion, and the text itself later floats weak resonance binding.\n\nWhat the paper does well: for three mass regimes it sets up symmetrized four-particle Gaussian wavefunctions, evaluates the matrix elements analytically, varies exponents and center separation, and handles the center-of-mass contamination explicitly. That is reproducible from the appendix and works as a nice pedagogical demonstration that small Coulombic molecules can be treated beyond fixed nuclei without heavy machinery. The H2 and di-muonated H2 numbers are sensible order-of-magnitude results for such a crude basis, and the self-citations to [6,7] are appropriate because those papers use the same model.\n\nThe soft spot is the one your reader flagged, and it is load-bearing. Section 4 says there is no binding energy for mono-muonated dihydrogen because the optimized trial energy at finite c lies above the c→∞ energy. But that trial energy is only an upper bound to the exact ground-state energy; a higher value cannot exclude a deeper exact state. The basis is one Gaussian per particle, and the energy difference at stake is ~0.005 a.u. in a total of ~−49.6 a.u., so basis error is likely much larger than the margin being probed. The authors themselves warn that the mono-muonated curves carry no quantitative meaning for c<0.06, which covers the region where the alleged resonance would appear. Then, in the center-of-mass section, they say it seems very likely the resonance effect will lead to weak binding between the non-muonic and muonic hydrogen atoms—directly contradicting Section 4. The Section 5 claim that BO holds only for unsubstituted H2 inherits this fragility, because it relies on the orbital-radii picture rather than on a quantitative estimate of the BO error.\n\nDespite that, the analytic derivation is sound and the variational framework itself is not flawed; the problem is over-interpretation. I think this deserves a serious referee, with the expectation of major revision: the mono-muonated result should be reframed as an upper-bound illustration, and the BO-validity claim should be either removed or heavily qualified. I would not cite the no-bond conclusion in my own work, but I might bring this short paper to a reading group to discuss where variational bounds can and cannot support physical conclusions.","headline":"A clean variational exercise with analytic Gaussian integrals, but the no-bond claim for mono-muonated H2 overreaches what an upper-bound calculation can show, and the paper contradicts itself on weak resonance binding.","tokens_in":9915,"tokens_out":3057,"would_cite":false,"duration_ms":34122,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["36.10.-k","31.15.vj"],"model":"deepseek-v4-flash","headline":"The paper argues that a mass swap from electron to muon changes whether a hydrogen-like molecule binds: ordinary di-hydrogen and di-muonated di-hydrogen bind, while the mixed mono-muonated molecule does not, and the Born-Oppenheimer…","keywords":["Born-Oppenheimer approximation","muonated hydrogen","muon-catalyzed fusion","variational calculation","Gaussian orbitals","four-fermion Coulomb system","molecular binding","non-Born-Oppenheimer effects"],"falsifier":"Run a high-accuracy non-Born-Oppenheimer calculation on the system with two protons, one electron, and one muon, using a many-Gaussian or explicitly correlated basis. If the optimized ground-state energy falls below the separated-atom limit, the paper's no-bond conclusion is falsified; the same calculation's one-particle densities would show whether the muon and proton orbital radii remain comparable, directly testing the radii-based Born-Oppenheimer-validity criterion.","tokens_in":8895,"feed_emoji":"⚛️","tokens_out":10943,"duration_ms":100650,"temperature":0.7,"pith_summary":"To mark the 90th anniversary of the Born-Oppenheimer approximation, the paper asks when a four-fermion Coulomb molecule still binds and when the approximation applies. Using a symmetrized product of one Gaussian orbital per particle and varying particle masses, it reports that ordinary di-hydrogen binds with 0.126 a.u. and that di-muonated di-hydrogen binds much more strongly at a bond distance of 0.014 a.u. The mixed molecule, with one electron and one muon, shows no binding: the optimized energy at finite separation sits above the separated-atom value because the muon breaks the degeneracy among permuted configurations, suppressing the resonance that lowers the energy in the symmetric molecules. From the computed orbital radii, the paper concludes that the Born-Oppenheimer approximation is valid only for unsubstituted di-hydrogen, where the proton is much heavier than the electron. The question matters because muon-substituted hydrogen is central to muon-catalyzed fusion, and a failure of the Born-Oppenheimer approximation would change how such systems must be modeled.","feed_headline":"Muon swap dissolves hydrogen's chemical bond","feed_subtitle":"A four-particle variational model says the Born-Oppenheimer approximation fails once muons replace electrons.","key_machinery":"The machinery is a variational wavefunction made of four Gaussian orbitals, one per particle, grouped in pairs around two atomic centers separated by a vector $c$; each orbital has its own width, and the product is symmetrized over exchanges of identical particles. All Coulomb and kinetic expectation values reduce to closed-form Gaussian integrals, several involving erf functions, so the energy can be minimized as a function of the widths and $c$ with no numerical quadrature. Translational contamination is removed by subtracting the center-of-mass kinetic energy, defining $H_{\\rm int}=H-P^2/(2M_{\\rm tot})$. The decisive mechanism is degeneracy: when the two leptons have equal mass, the four permuted configurations share one energy and can resonate to bind the molecule, but in the mixed electron-muon system that degeneracy is broken and the resonance effect disappears; the resulting orbital radii provide the paper's criterion for whether the Born-Oppenheimer approximation is trustworthy.","core_discovery":"On its own terms, the central claim is that the binding of these four-particle molecules is controlled by a mass-symmetry switch. In H2 and in fully muon-substituted H2, the four wavefunctions obtained by exchanging identical leptons and identical protons are degenerate, so their linear combinations produce a resonance that lowers the energy; the model yields binding energies of 0.125949 a.u. and 7.6337 a.u. (7.38 a.u. after center-of-mass correction) at bond distances of 1.626 a.u. and 0.0142 a.u., respectively. In the mixed system, the muon and electron configurations have different energies, the resonance is suppressed, and the finite-separation variational energy never falls below the separated-atom energy, so the paper concludes there is no chemical bond. The same wavefunction's orbital radii are then used as a diagnostic: the electron orbit radius in H2 is about six times the proton orbit radius, while in muonated species the heavy and light radii are comparable, so fixing the heavy particles as delta-function sources is no longer justified. The proposed conclusion is that the Born-Oppenheimer approximation is valid only for ordinary di-hydrogen.","pith_inferences":["Because a variational upper bound cannot rule out a deeper exact energy, the no-bond conclusion for mono-muonated di-hydrogen is not settled: a multi-Gaussian or explicitly correlated calculation could still find a shallow bound state, and the paper's result should be read as a model-based prediction rather than a proof.","The radii-based BO-validity criterion could be made quantitative by sweeping the lepton mass continuously from electron to muon and locating where the proton and lepton orbital radii cross or where binding disappears; the resulting threshold would be a testable prediction for other exotic molecules.","If the criterion generalizes, then any exotic atom whose lightest charged particle has a mass approaching its nuclear partner's mass, such as other heavy-lepton or antimatter systems, would lie outside the Born-Oppenheimer regime and require nonadiabatic treatment."],"forward_implications":["If the Born-Oppenheimer approximation is invalid for mono- and di-muonated hydrogen, then muon-catalyzed fusion models that rely on potential-energy surfaces should be replaced by a full four-particle dynamical treatment.","The di-muonated molecule's equilibrium distance of about 0.014 a.u. puts the two protons roughly 100 times closer than in H2, and its tiny moment of inertia implies no rotational spectrum, consistent with the nuclear-compression picture of muon catalysis.","The mixed molecule, even if not chemically bound, is said to be near a resonance between configurations with almost equal energies, with the protons very close together; the paper suggests this resonance may provide a weak binding relevant to catalysis.","The closed-form single-Gaussian variational treatment, including the center-of-mass correction, provides a reusable recipe for other four-fermion Coulombic systems such as dipositronium."],"supporting_citations":[{"why":"Defines the Born-Oppenheimer approximation whose applicability to the three molecules is the paper's subject.","marker":"[1]"},{"why":"Supplies the physical motivation: replacing an electron by a muon brings nuclei about 200 times closer, enabling muon-catalyzed fusion.","marker":"[3]"},{"why":"Provides the variational strategy for four-particle Gaussian orbitals and the observation that the proton-to-muon mass ratio is only about 10.","marker":"[6]"},{"why":"Earlier conclusion that muonated molecules are not reliably described by the Born-Oppenheimer approximation and that the chemical bond is a resonance effect absent in the mono-muonated case.","marker":"[7]"},{"why":"Shows how to remove translational motion from nuclear-orbital wavefunctions, informing the center-of-mass correction used here.","marker":"[11]"},{"why":"Gives the erf integral identity used to evaluate all the Gaussian Coulomb integrals in closed form.","marker":"[12]"}],"fun_headline_variants":["Muon swap unbinds hydrogen molecule","Muonated H2 loses chemical bond","Born-Oppenheimer fails for muon H2","Muon for electron breaks H2 bond","Mass symmetry controls hydrogen bonding"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on trusting the single-Gaussian trial wavefunction to represent the true ground state: because it is a variational upper bound, a more flexible wavefunction could still find a bound mixed molecule, and the same model supplies the orbital radii that motivate the Born-Oppenheimer-validity claim.","fun_headline_variants_meta":{"raw":{"variants":["Muon swap unbinds hydrogen molecule","Muonated H2 loses chemical bond","Born-Oppenheimer fails for muon H2","Muon for electron breaks H2 bond","Mass symmetry controls hydrogen bonding"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00081,"raw_usage":{"total_tokens":3517,"prompt_tokens":873,"completion_tokens":2644,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":489,"completion_tokens_details":{"reasoning_tokens":2580}},"tokens_in":489,"tokens_out":2644,"duration_ms":18449,"temperature":1.0,"reasoning_tokens":2580,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:54:14.440269+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a high-accuracy non-Born-Oppenheimer calculation on the system with two protons, one electron, and one muon, using a many-Gaussian or explicitly correlated basis. If the optimized ground-state energy falls below the separated-atom limit, the paper's no-bond conclusion is falsified; the same calculation's one-particle densities would show whether the muon and proton orbital radii remain comparable, directly testing the radii-based Born-Oppenheimer-validity criterion.","supporting_citations":[{"cited_title":"Zur quantentheorie der moleke n","cited_arxiv_id":null,"evidence_quote":"Defines the Born-Oppenheimer approximation whose applicability to the three molecules is the paper's subject."},{"cited_title":"µ-mesonic molecules. 2. Molecular-Ion For- mation and Nuclear Catalysis","cited_arxiv_id":null,"evidence_quote":"Supplies the physical motivation: replacing an electron by a muon brings nuclei about 200 times closer, enabling muon-catalyzed fusion."},{"cited_title":"Or- bitals of the dipositronium","cited_arxiv_id":null,"evidence_quote":"Provides the variational strategy for four-particle Gaussian orbitals and the observation that the proton-to-muon mass ratio is only about 10."},{"cited_title":"On dipositronium and molecular hydrogen: similarities and differences","cited_arxiv_id":null,"evidence_quote":"Earlier conclusion that muonated molecules are not reliably described by the Born-Oppenheimer approximation and that the chemical bond is a resonance effect absent in the mono-muonated case."},{"cited_title":"Non-Borm Openheimer effects predicted by translation free nuclear orbitals plus molecular orbital method","cited_arxiv_id":null,"evidence_quote":"Shows how to remove translational motion from nuclear-orbital wavefunctions, informing the center-of-mass correction used here."},{"cited_title":"Abramowitz and I.A","cited_arxiv_id":null,"evidence_quote":"Gives the erf integral identity used to evaluate all the Gaussian Coulomb integrals in closed form."}],"review_version":1}