{"id":"4431d488-5974-4a90-ae5c-a96885ac5302","arxiv_id":"2412.18747","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A muon bound to a noble-gas atom undergoes orbital collapse at a critical angular momentum, limiting which states can be populated by capture.","lead":"Muonic argon atoms show a sudden orbital collapse: muon states with angular momentum above a critical value lc no longer overlap the electrons that would capture them. The paper predicts lc with a simple mass-scaled formula and tests the resulting upper limit with three computational methods.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The hard upper limit on captured muon angular momentum rests on the single-active-electron/static-potential approximation; multi-electron or dynamic-screening channels could populate l>37, so the cutoff needs a multi-electron check.","rationale":"After reading the manuscript, the central claim is precisely the hard cutoff at lc=37. The DFT energy-structure result is direct, but the dynamical confirmation (TDSE/CTMC) uses the same single-active-electron model potentials that enter the prediction, so it cannot independently rule out contributions from other electrons or from dynamic screening. The paper itself lists these as limitations. The reader's CONDITIONAL verdict already captures this; my stress test does not identify a different, more severe flaw. The remaining uncertainty is whether multi-electron channels populate l>37; a two-active-electron or TDDFT calculation, or a quantitative overlap calculation with all occupied orbitals, would settle it. If the multi-electron capture probability remains negligible, the claim stands; if not, the verdict should move toward rejection or unverified. Therefore I recommend UNCHANGED (conditional remains appropriate).","tokens_in":8974,"tokens_out":16681,"duration_ms":180266,"concrete_test":"Run the same muon-Ar collision at 10 eV with two active electrons (3s and 3p) in the TDSE, or with a time-dependent DFT treatment of all 18 electrons, and compare the total capture probability into l>37 states with the 3p-only result of Fig. 4(b). Also compute the squared overlaps S_l = sum_{occupied e} |<psi_mu(n,l)|phi_e>|^2 for l=37 and 38; if P(l>37) or S_38 is more than ~1% of the peak capture probability, the hard upper limit is violated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that lc=37 is a hard upper limit on the angular momentum of muon-captured states in Ar—is established only inside the single-active-particle, static-model-potential framework. The paper states 'We consider only the 3p electron as an active electron' and the TDSE/CTMC simulations evolve this one electron and the muon in a potential fitted to the neutral-atom ground state. The authors concede that other orbitals (3s, 2p) contribute (with 'lower probabilities') and that the model potential 'does not consider dynamic changes of the electron wavefunction when the heavy particle locates in different orbitals.' Since the upper limit is a statement about all capture channels, a nonzero multi-electron or polarization-induced contribution to l>37 would invalidate it. The qualitative 'almost no overlap' argument from Fig. 1 is not quantitative: the final outer muon wavefunction has a small but nonzero tail into the electron cloud, and the capture amplitude is a scattering matrix element, not just a bound-state overlap. Thus the dynamical confirmation does not close the channel; it assumes it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a DFT study of muonic Ar atoms in which the muon orbital is found to undergo a collapse at a critical angular momentum lc (about 36–37 for Ar): for circular states with lµ > lc the muon wavefunction is pushed outside the electron cloud, while for lµ = lc it penetrates inside. The authors interpret lc as an upper limit on the angular momentum of muon-captured states and confirm this with TDSE and CTMC simulations using a single-active-electron model with model potentials. They also propose a semiclassical scaling relation, lc(lc+1) = mµ Lc^2, where Lc is the local maximum of the classical circular-orbit angular momentum in an atomic model potential, and use it to estimate lc for muonic and antiprotonic noble-gas atoms, with values in Table I compared against DFT.","tokens_in":9227,"tokens_out":9693,"duration_ms":83831,"significance":"If correct, this is a qualitatively new prediction: orbital collapse of a heavy negative particle in a low-Z atom, giving a sharp upper bound on the angular momentum of initially captured states. Because the initial capture state determines the subsequent x-ray cascade, the proposed upper limit is directly relevant to recent and planned muonic-atom experiments (e.g., muonic Ar and Ne). The paper offers a simple, potentially widely applicable scaling estimate for lc, and it cross-checks the dynamics with two independent methods (TDSE and CTMC). The strengths are the direct DFT wavefunction evidence, the explicit cross-method dynamical confirmation, and the authors' candid acknowledgment of the single-active-electron and static-potential limitations. The main weaknesses are that the dynamical confirmation is carried out entirely within those approximations, and that the scaling method's agreement with DFT is partly a consistency check because the model potential is fitted to the same DFT potential.","major_comments":[{"comment":"The values in the 'lc (µ: scaled)' column do not follow from Eq. (6) with the stated rule of rounding to the nearest integer. For Ne, Lc = 1.98 gives lc(lc+1) ≈ 1.98^2 × 207 ≈ 811, whose nearest integer is lc = 28, not 27; for Ar, Lc = 2.65 gives lc(lc+1) ≈ 1453, and the nearest integer is lc = 38, not 37. Please specify the exact value of mµ and the rounding convention used, and correct the table or the text accordingly; as printed, the agreement between the scaled and DFT columns is not reproducible.","section":"Table I / Eq. (6)"},{"comment":"The claim that lc provides an upper limit for muon-captured states is established only within the single-active-electron approximation with static model potentials, as the authors themselves state ('We consider only the 3p electron as an active electron' and the potential 'does not consider dynamic changes of the electron wavefunction when the heavy particle locates in different orbitals'). Since the abstract and conclusion phrase this as a general property of muon-captured states, please either qualify the upper limit as model-dependent or provide a quantitative estimate of the contribution to l > lc from other active orbitals (3s, 2p) and from dynamic response of the electron cloud. Without such a test, the statement that capture into lµ > lc is 'negligible' goes beyond what the presented calculations can demonstrate.","section":"Abstract / §2 (single-active-particle)"},{"comment":"The 'simple method' of Eq. (6) is validated in Table I against DFT calculations, but model potential 1 is fitted to the same self-interaction-free DFT potential (Ref. [34]); the agreement is therefore a consistency check rather than an independent prediction. Please report Lc and the resulting lc obtained with model potential 2 (fitted to a Hartree-Fock potential) as well, or otherwise show that the scaling prediction is insensitive to the choice of model potential and to the fitting target.","section":"§4 / Table I"}],"minor_comments":[{"comment":"Equation (2) as printed is ambiguous: the placement of parentheses is unclear, and the current typesetting appears to show −(Z−1)[(η/ξ)(e^{ξr}−1)+1]^{−1}+1 all divided by r, which is not dimensionally consistent. Please rewrite the formula unambiguously.","section":"Eq. (2)"},{"comment":"The capture probabilities in Fig. 4 are normalized to the largest value; to support the claim that states with lµ ≥ lc are 'negligibly' populated, please report absolute capture probabilities or at least a numerical upper bound relative to the peak for the lµ ≥ lc region.","section":"Fig. 4"},{"comment":"The caption of Fig. 1 is difficult to parse; please rewrite it to state clearly which line style corresponds to which muon state and which to the 3p electron.","section":"Fig. 1 caption"},{"comment":"The text uses 'principle quantum number' in the introduction; the correct spelling is 'principal quantum number'.","section":"Introduction"},{"comment":"Please specify in the text whether mµ in Eq. (6) is the muon mass in atomic units or the reduced mass of the muonic atom, and use the same convention throughout the paper.","section":"Eq. (6)"}],"recommendation":"major_revision","confidential_remarks":"The core phenomenon—a DFT-predicted orbital collapse of a muon in Ar with lc near 37—appears plausible and interesting, and the dynamical cross-check with TDSE and CTMC is a valuable addition. However, the numerical inconsistency in Table I with Eq. (6) must be fixed, and the paper should more carefully delineate the scope of the 'upper limit' claim given the single-active-electron approximation. The paper is within scope for a physics journal. I recommend major revision rather than rejection because the central idea is sound and the issues are addressable within the manuscript's framework."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the prediction that a muon in a noble-gas atom undergoes orbital collapse at a critical angular momentum (l_c ~ 37 for Ar), so that capture into states with l > l_c is negligible. The scaling rule l_c(l_c+1) = m_mu L_c^2 is simple, non-obvious, and testable. The paper earns credit for that.\n\nWhat it does well: the DFT wavefunctions directly show the collapse at the circular-state level; the classical L(r) barrier argument is clean; and both the TDSE and CTMC calculations produce negligible capture probability for l > l_c, which is a real dynamical check. Table I, comparing l_c from the scaling formula with DFT results across Ne, Ar, Kr, Xe, is encouraging. The CTMCm remapping is post hoc, but the authors say why it is needed and it brings the two methods into reasonable agreement. Citation practice is fine; electron-collapse literature is there.\n\nWhere it is softer: the phrase \"upper limit\" is stronger than the evidence. The whole calculation uses a single-active-electron approximation (only the 3p electron) and a static model potential fitted to ground-state DFT or HF. The TDSE and CTMC dynamics then use those same model potentials, so the dynamical confirmation is not independent of the potential that creates the barrier. Figure 1's \"almost no overlap\" is qualitative; capture is a scattering matrix element, not just a bound-state overlap, and a small tail in the muon wavefunction can still contribute. The authors concede that 3s and 2p channels contribute with lower probability and that the model potential ignores dynamic changes of the electron cloud. Those channels could populate l > l_c at some small level. So the cutoff is probably right for the dominant capture path, but as a strict upper bound on all channels it is not proven. Agreement with the muonic Ne experiment is interpretive, not a direct test. Data being available only on request is a minor but unnecessary obstacle.\n\nOverall, the central mechanism is plausible and the paper is worth engaging. The stress-test concern about multi-electron effects is legitimate but not fatal; a multi-electron calculation or a softened claim would address it. I would send this to peer review with a request to revise the upper-limit language and deposit the potentials and raw capture probabilities. The paper is useful for the muonic-atom and antiprotonic-atom community, and I would cite the scaling rule as a proposed mechanism rather than an established bound.","headline":"A genuinely new and testable prediction of a mass-scaled orbital collapse that caps muon capture angular momentum in noble gases, even though the 'hard upper limit' is softer than the calculations can strictly support.","tokens_in":747,"tokens_out":1884,"would_cite":true,"duration_ms":30027,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["36.10.-k"],"model":"deepseek-v4-flash","headline":"This paper argues that a heavy negative particle captured into an atom can undergo orbital collapse in a light atom, something previously associated only with f- or d-orbitals in heavy elements.","keywords":["muonic atoms","orbital collapse","muon capture","exotic atoms","angular momentum ceiling","time-dependent Schrödinger equation","classical trajectory Monte Carlo","density functional theory"],"falsifier":"Measure the initial-state distribution of muonic Ar—for instance from the electronic K x-ray spectrum or from muonic cascade lines—and look for any population with $l_\\mu\\ge37$; finding a non-negligible high-$l$ component would falsify the ceiling. A numerical counterpart is a TDSE calculation with the 3s and 2p electrons active, checking whether any flux lands in $l\\ge37$ states.","tokens_in":8758,"feed_emoji":"⚛️","tokens_out":10053,"duration_ms":90380,"temperature":0.7,"pith_summary":"This paper argues that a heavy negative particle captured into an atom can undergo orbital collapse in a light atom, something previously associated only with f- or d-orbitals in heavy elements. The central claim is that in muonic argon the muon orbital collapses at a critical angular momentum $l_c\\approx37$: the circular state $(37,36)$ lies outside the 3p electron cloud with almost no overlap, while $(37,35)$ and lower states shrink into the inner well, so states with $l_\\mu\\ge l_c$ are essentially never populated in capture. The authors confirm the ceiling by solving the time-dependent Schrödinger equation and by a classical trajectory Monte Carlo method in the single-active-particle approximation, and they propose a simple estimate $l_c(l_c+1)=m_\\mu L_c^2$ from the local maximum of the classical circular-orbit angular momentum in an atomic model potential. If correct, the result fixes the initial angular momentum distribution of muonic and antiprotonic noble-gas atoms, which is what subsequent x-ray cascades begin from.","feed_headline":"Muon orbit in argon collapses at angular momentum 37","feed_subtitle":"The collapse caps which states a captured muon can occupy, reshaping predicted x-ray cascades.","key_machinery":"The load-bearing identity is the scaling relation $l_c(l_c+1)=m_\\mu L_c^2$, where $L_c$ is the local maximum of the classical circular-orbit angular momentum $L(r)=\\sqrt{-2mr^2V(r)}$ for a particle of mass $m$ in an atomic model potential $V(r)$. It converts a purely classical feature—the jump of the stable circular orbit from an inner radius $R_c$ to an outer radius $R_0$ as $L$ passes $L_c$—into a quantum angular momentum at which the muon's centrifugal barrier no longer allows an inner bound state. The accompanying quantum object is the effective potential $V_{\\rm eff}(r)=L_e^2/(2r^2)+V(r)=l_\\mu(l_\\mu+1)/(2m_\\mu r^2)+V(r)$, whose barrier height for Ar corresponds to $l_e\\approx2.15$, a non-integer for electrons but an allowed state, $l_\\mu=37$, for a muon. The dynamical calculations are built on the single-active-particle approximation, with the 3p electron and the muon sharing one model potential, and on the CTMCm remapping of Eq. (8).","core_discovery":"The central discovery is that the muon orbital changes character discontinuously at $l_c$: in muonic Ar, the $(37,36)$ circular orbital is localized outside the 3p electron orbital and has almost no overlap with it, whereas the $(37,35)$ orbital is confined inside and overlaps strongly. The mechanism is the effective radial potential $V_{\\rm eff}(r)=l_\\mu(l_\\mu+1)/(2m_\\mu r^2)+V(r)$, which develops a centrifugal barrier separating an inner and an outer well. For $l_\\mu=37$ the inner well cannot hold a bound state, while for $l_\\mu=36$ it can; the same barrier corresponds to an electron angular momentum $l_e\\approx2.15$, which is forbidden for argon, explaining why electrons do not show this collapse. Consequently the state-specified muon capture probability into states with $l_\\mu\\ge l_c=37$ is negligibly small in both TDSE and CTMC simulations; the two methods agree reasonably after CTMC's classical binding energies are remapped to quantum states using the eigenenergies of the radial Schrödinger equation. The authors show the phenomenon is general, giving $l_c\\approx27,37,43,50$ for muonic Ne, Ar, Kr, Xe and larger values for antiprotonic atoms.","pith_inferences":["A multi-active-electron TDSE calculation that includes the 3s and 2p electrons would test whether those channels populate states above $l_c$; if they do, the ceiling becomes a statement about the dominant channel rather than a strict upper bound.","The small discrepancies between the scaled and DFT values of $l_c$ in Table I suggest that allowing the electron cloud to respond dynamically as the muon moves inward could shift $l_c$ by one or two units.","The same recipe—locate $L_c$ in any model potential and apply $l_c(l_c+1)=m L_c^2$—could predict capture ceilings for pions, kaons, or antihydrogen formation, yielding a testable mass-scaling sequence.","A direct experimental falsification would be to infer the initial $l_\\mu$ distribution of muonic Ar from the electronic K x-ray spectrum; the paper's ceiling predicts no population at or above $l_\\mu=37$."],"forward_implications":["Muonic Ar capture populations are confined below $l_c$; cascade calculations can safely omit high-$l$ initial states.","Since $l_c$ scales with the square root of the particle mass, antiprotonic noble gases have ceilings near 80–150, making the effect a general feature of heavy-particle capture.","The remapped CTMC method (CTMCm) gives state-specified capture probabilities in reasonable agreement with TDSE, providing a practical computational tool for other exotic atoms.","Orbital collapse appears already in muonic Ne, and the trend across Ne, Ar, Kr, and Xe matches the implication of earlier muonic Ne experiments that the initial $n_\\mu$ exceeds 20."],"supporting_citations":[{"why":"Measured electronic K x rays from muonic Ar show initial capture states higher than the mass-scaling value, motivating and constraining the present calculation.","marker":"[29]"},{"why":"Supplies the self-interaction-free DFT potential and muon orbital densities used to identify the collapse at $l_c$.","marker":"[32]"},{"why":"Defines model potential 1, the six-parameter form whose circular-orbit maximum gives $L_c$ and hence $l_c$.","marker":"[33]"},{"why":"Provides the Hartree-Fock parameters for model potential 2, the two-parameter potential used in the CTMC and TDSE simulations.","marker":"[37]"},{"why":"Establishes the TDSE approach for Coulomb three-body capture that is extended here to muon capture in argon.","marker":"[24]"},{"why":"Gives the TDSE method and muon–H capture benchmark from which the present muon-capture TDSE is adapted.","marker":"[38]"},{"why":"Originates the classical trajectory Monte Carlo method used for the capture-probability calculations.","marker":"[25]"},{"why":"Shows the binding-energy to quantum-state mapping for antiproton capture that the authors modify into CTMCm.","marker":"[39]"},{"why":"Fermi's Thomas-Fermi orbital-collapse prediction is the phenomenon the paper generalizes to exotic-atom particles.","marker":"[1]"}],"fun_headline_variants":["Muon orbit collapse in argon caps captured states","Argon muon orbital collapse at l=37 limits capture","Orbital collapse in exotic atoms alters muon dynamics","Muonic argon: l=37 collapse affects x-ray cascades","Exotic atom orbital collapse sets muon capture limit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes argon can be represented by one active 3p electron and a muon moving in a single, fixed model potential, so that capture through the 3s or 2p electrons or through dynamic response of the electron cloud is too weak to populate states above $l_c$; if that assumption fails, the claimed hard ceiling softens.","fun_headline_variants_meta":{"raw":{"variants":["Muon orbit collapse in argon caps captured states","Argon muon orbital collapse at l=37 limits capture","Orbital collapse in exotic atoms alters muon dynamics","Muonic argon: l=37 collapse affects x-ray cascades","Exotic atom orbital collapse sets muon capture limit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000372,"raw_usage":{"total_tokens":2000,"prompt_tokens":966,"completion_tokens":1034,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":582,"completion_tokens_details":{"reasoning_tokens":953}},"tokens_in":582,"tokens_out":1034,"duration_ms":8773,"temperature":1.0,"reasoning_tokens":953,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:30:58.022039+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the initial-state distribution of muonic Ar—for instance from the electronic K x-ray spectrum or from muonic cascade lines—and look for any population with $l_\\mu\\ge37$; finding a non-negligible high-$l$ component would falsify the ceiling. A numerical counterpart is a TDSE calculation with the 3s and 2p electrons active, checking whether any flux lands in $l\\ge37$ states.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the self-interaction-free DFT potential and muon orbital densities used to identify the collapse at $l_c$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Hartree-Fock parameters for model potential 2, the two-parameter potential used in the CTMC and TDSE simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the TDSE approach for Coulomb three-body capture that is extended here to muon capture in argon."},{"cited_title":"Green, An analytic independent particle model for atoms, in Advances in Quantum Chemistry, Vol","cited_arxiv_id":null,"evidence_quote":"Gives the TDSE method and muon–H capture benchmark from which the present muon-capture TDSE is adapted."},{"cited_title":"Abrines and I","cited_arxiv_id":null,"evidence_quote":"Originates the classical trajectory Monte Carlo method used for the capture-probability calculations."},{"cited_title":"Tokesi, B","cited_arxiv_id":null,"evidence_quote":"Shows the binding-energy to quantum-state mapping for antiproton capture that the authors modify into CTMCm."}],"review_version":1}