{"id":"d0592ecf-8a87-4efb-8667-5f2926255e30","arxiv_id":"1909.02141","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Triatomic butterfly Rydberg molecules with two ground-state atoms are predicted to bind in two parity classes with distinct geometries and vibrational spectra.","lead":"Ultralong-range Rydberg molecules can bind a second atom, forming triatomic butterfly molecules. The paper maps their shapes, vibrational spectra, and a design rule for where they sit.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No convergence study supports the asserted n=29-31/lmin=2 hybrid-basis truncation; the predicted trimer wells and vibrational spacings could shift if more manifolds or partial waves are included.","rationale":"The paper is a coherent theoretical extension of known butterfly dimers to trimers. The machinery—Fermi-Omont pseudopotential, hybrid dimer-orbital/quantum-defect basis, Born-Oppenheimer PES, and DVR/finite-difference vibrational calculation—is standard in this subfield, and the analytic treatment of the odd-parity class provides independent support. The parity decoupling and bosonic symmetry constraints are internally consistent. The central claim does not depend on agreement with external consensus; it is a new prediction. The weakest point is exactly the asserted but undemonstrated convergence of the basis truncation, which is also what the reader identified. I considered whether neglect of spin structure or the L=0 projection might be more serious, but the paper acknowledges spin as a future quantitative refinement, and L=0 projection is standard for vibrational spectra of triatomic Rydberg molecules. I also considered whether the qualitative 'building principle' for even trimers is too phenomenological; it is presented as an interpretive guide, not as a premise of the binding claim. Thus no fatal flaw is apparent. The conditionality is appropriate: the paper should supply a convergence study or ship code/parameter files before the quantitative spectra are taken as reliable. Because the reader's verdict is already CONDITIONAL, no change is needed.","tokens_in":11181,"tokens_out":4252,"duration_ms":48434,"concrete_test":"Recompute the representative PES cuts of Figs. 2-5 with an enlarged hybrid basis, e.g., M=5 manifolds (n=28-32) and lmin=3, and compare (i) odd-trimer minima at theta=pi, (ii) even-trimer minima at theta=pi and theta=0.35pi, and (iii) the five lowest vibrational energies for each parity class. If all minima positions shift by less than about 10 a0 and energies by less than about 1 GHz, the n=29-31/lmin=2 truncation is adequate and the central claim stands; if larger shifts occur, the reported spectra and even the binding energies are unconverged.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that resonant p-wave scattering binds a second ground-state atom to a Rydberg core, forming triatomic butterfly molecules with two parity classes and computable vibrational spectra. For this claim to hold quantitatively, the potential energy surfaces must be converged with respect to the hybrid-basis truncation. Section 2 states: 'In our present calculations we use the n = 29, 30, and 31 Rydberg manifolds and include quantum defects for s, p, and d waves (lmin = 2). These parameters give adequately converged potential energy surfaces.' No convergence study, table, or error estimate supports this assertion. This is load-bearing for two reasons. First, the butterfly PES is made finite only through coupling to neighboring Rydberg manifolds and to the low-l quantum-defect states; the p-wave scattering volume diverges at resonance, so the truncation is not a benign high-l cutoff but part of the regularization. Second, the paper's quantitative outputs—well depths, minima positions, and vibrational spacings such as 30 MHz, 400 MHz, 1.2 GHz, and 2.5 GHz—are all read off these surfaces. If adding n=28 or n=32 shifts a well by more than its depth, the existence or localizability of the predicted trimer states is not established by the calculation as presented. The manuscript itself flags the omission only indirectly by asserting adequacy, and no code or parameter files are shipped to allow direct reproduction. This is an addressable numerical-convergence concern, not a conceptual flaw; it does not by itself falsify the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a theoretical study of triatomic 'butterfly' Rydberg molecules, in which two ground-state atoms are bound to a Rydberg core by resonant p-wave scattering. The authors construct Born-Oppenheimer potential energy surfaces from the Fermi-Omont pseudopotential using a hybrid basis of dimer orbitals and quantum-defect states for n=29, 30, and 31. They divide the trimers into odd and even parity classes, give analytic expressions for the odd surfaces in Eqs. (12)-(14), and analyze the even surfaces using correlations with dimer orbital gradients. Nuclear eigenstates are computed with a finite-difference/DVR approach, yielding vibrational spectra with spacings around 30 MHz, 400 MHz, 1.2 GHz, and 2.5 GHz. The paper concludes that butterfly trimers exist in two classes with qualitatively different potential landscapes and vibrational dynamics.","tokens_in":11502,"tokens_out":4003,"duration_ms":49984,"significance":"If the predicted states are real, this work significantly extends ultra-long-range Rydberg molecules to a new triatomic class and provides a useful symmetry-based framework for future studies. The analytic form of the odd-trimer surfaces and the careful separation into parity classes are valuable and go beyond prior dimer work. The main quantitative predictions, however, rest on a convergence assertion that is not demonstrated; the claimed well depths and vibrational frequencies should therefore be treated as provisional.","major_comments":[{"comment":"The statement 'These parameters give adequately converged potential energy surfaces' is not backed by any convergence study. The hybrid basis truncates to n=29-31 and lmin=2, and the paper's quantitative outputs, including minima positions in Figs. 3-6 and the vibrational spacings in Section 5, are read off these surfaces. Because the p-wave scattering volume diverges and is regularized only through coupling to neighboring Rydberg manifolds and quantum-defect states, this truncation is not a benign high-l cutoff. Please provide convergence data for, e.g., n=28/32 and lmin=3, and quantify how well depths and vibrational frequencies shift under these changes.","section":"Section 2, paragraph after Eq. (11)"},{"comment":"The Hamiltonian in Eq. (1) neglects all spin degrees of freedom and assumes scattering occurs only in the triplet channel. The conclusion acknowledges that including spin is 'necessary for quantitative predictions,' yet Section 5 reports vibrational frequencies such as 30 MHz, 400 MHz, 1.2 GHz, and 2.5 GHz without any estimate of the spin-induced corrections. The authors should either estimate the effect of spin on the lowest vibrational states or explicitly label the quoted spectra as qualitative rather than quantitative.","section":"Section 2, Eq. (1) and Section 6"},{"comment":"The hybrid-basis method is central to the calculation, but the manuscript does not state the numerical parameters of the diagonalization, such as grid sizes, number of dimer orbitals per manifold, or the tolerance used in the generalized eigenproblem. Without these details, an independent check of the potential energy surfaces is difficult. Please add a brief numerical-convergence statement covering these parameters.","section":"Section 2, Eq. (7)"}],"minor_comments":[{"comment":"'Rb 3 potential surface' should be 'Rb trimer potential surface' to avoid ambiguity with a three-body system vs. a state label.","section":"Figure 2 caption"},{"comment":"In the sentence 'very similiar to trilobite trimers', 'similiar' should be 'similar'.","section":"Section 5, last paragraph"},{"comment":"The overlap integral in Eq. (3) uses the label n for the principal quantum number, but the dependence on n is not explicitly retained in the notation; please clarify that the sums are over l>lmin and m for a fixed n.","section":"Section 2, Eq. (3)"}],"recommendation":"major_revision","confidential_remarks":"The missing convergence study is the main blocker; everything else is addressable in revision. The paper fits the journal's scope and the topic is timely, but the quantitative claims should not be accepted without numerical evidence. A short reproducibility statement would also strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid extension of the Rydberg trimer work to p-wave-bound butterfly molecules, and it deserves refereeing. The thing to know: the central qualitative claim—that resonant p-wave scattering can bind a second ground-state atom and give two parity classes—holds up. The quantitative numbers (well depths, the 30 MHz / 400 MHz / 1.2 GHz / 2.5 GHz spacings) rest on a convergence assertion that the paper does not back with data.\n\nWhat is actually new: first calculation of the full 3D potential energy surface for triatomic butterfly molecules, the odd/even reflection-parity classification, analytic expressions for the odd-parity surfaces (Eqs. 12–14), and a gradient-based building principle for the even trimers. The vibrational mode ordering (bending stiffer than stretch) differs from s-wave and trilobite trimers, which gives experimentalists a concrete handle. The hybrid-basis machinery is used correctly, and the scattering volumes and quantum defects are fixed inputs from prior work, not fit to the target result. No circularity problem.\n\nSoft spots, in proportion. The load-bearing one: Section 2 states that the n=29–31, lmin=2 parameters give 'adequately converged potential energy surfaces,' but no convergence study is shown. This is not cosmetic. The butterfly potential is made finite only through coupling to neighboring Rydberg manifolds and low-l defect states; the p-wave scattering volume diverges at resonance, so the truncation is part of the regularization. If adding n=28 or n=32 shifts a well by more than its depth, the predicted vibrational states are not established by the paper as written. Addressable and probably fixable, but a real gap. Second, spin is neglected; the authors acknowledge this limits quantitative predictions. Acceptable for a first survey. Third, no code or parameter files are shipped, so the surfaces are not independently checkable without reimplementation. Minor: no error estimates anywhere.\n\nWho this is for: Rydberg molecule theorists and atomic physics experimenters looking for new ultralong-range molecular targets. This paper merits a serious referee; the referee should ask for a convergence check and some sensitivity analysis around the scattering model. I would not desk-reject it.","headline":"First full treatment of triatomic butterfly molecules, qualitatively sound but quantitatively provisional because the basis truncation convergence is asserted, not demonstrated.","tokens_in":12024,"tokens_out":2748,"would_cite":true,"duration_ms":25788,"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 paper claims that resonant p-wave scattering can bind a second ground-state atom to a Rydberg core, forming triatomic butterfly molecules, and that these come in two reflection-parity classes with distinct potential surfaces and…","keywords":["Rydberg molecules","butterfly molecules","p-wave scattering resonance","triatomic molecules","Born-Oppenheimer potential energy surfaces","vibrational spectra","Fermi pseudopotential","ultracold rubidium"],"falsifier":"Compute the same potential surfaces with n up to 32–35 (or with f-wave quantum defects included): if the deep odd-trimer minima near R1=R2 approximately 316 a0 and 340 a0 move by more than a few GHz or disappear, the truncation assumption fails. Alternatively, a microwave spectrum of Rb n approximately 30 dimers plus a second atom that shows no vibrational resonances near the predicted 400 MHz, 1.2 GHz, and 2.5 GHz spacings would refute the quantitative claim.","tokens_in":10981,"feed_emoji":"🦋","tokens_out":7184,"duration_ms":62868,"temperature":0.7,"pith_summary":"The paper sets out to show that resonant p-wave scattering between a Rydberg electron and a neutral atom—the mechanism behind diatomic butterfly molecules—can bind two ground-state atoms to one Rydberg core, creating triatomic butterfly molecules. It computes the three-dimensional Born-Oppenheimer potential energy surfaces and vibrational eigenstates of these trimers, and it separates them into two classes by the parity of the electronic wave function under reflection through the molecular plane. Odd trimers are predicted to be stable only near collinear geometries, with nearly independent stretching and bending modes and well-separated vibrational ladders; even trimers have many more potential minima, but only a few deep and isolated enough to localize vibrational states. The paper further claims a building principle: even-trimer minima sit where the second ground-state atom finds a steep gradient (rather than a maximum) of the diatomic butterfly orbital. If correct, this predicts a new class of ultra-long-range polyatomic Rydberg molecules with specific geometries, vibrational spacings, and large dipole moments.","feed_headline":"Two atoms bind to a Rydberg core as butterfly trimers","feed_subtitle":"Odd-parity trimers collapse to collinear shapes; even ones offer many minima but few wells deep enough to trap states.","key_machinery":"The central machinery is the Fermi-Omont pseudopotential Hamiltonian for a Rydberg electron interacting with two ground-state atoms, reduced to a compact hybrid basis of dimer orbitals—the trilobite (s-wave) and three butterfly orbitals (R-, $\\theta$-, phi-) per atom—plus quantum-defect states. The argument is carried by the overlap integrals between these orbitals: their nodal and gradient structure determines where trimer minima appear, and the reflection parity of the phi orbital splits the problem into odd and even trimers. For the odd trimers an explicit two-surface formula shows the cross-coupling vanishes like 1/$sin^{4}$($\\theta$), explaining why trimer features appear only near collinear geometry.","core_discovery":"Restricted to the n=30 butterfly manifold, the paper's central discovery is that the p-wave shape resonance can bind two ground-state 87Rb atoms to a Rydberg core, and that the resulting trimers split cleanly by reflection parity. Odd trimers (built from phi-dimer orbitals, which are odd under reflection through the molecular plane) have smooth potential surfaces away from collinearity; at $\\theta$=pi the cross-term coupling is large enough to create deep minima, and the paper finds vibrational states localized near equal bond lengths R1=R2 approximately 316 a0, 340 a0, and 270 a0 with harmonic-like ladders whose asymmetric-stretch, symmetric-stretch, and bending spacings are approximately 400 MHz, 1.2 GHz, and 2.5 GHz. Even trimers (built from R- and $\\theta$-dimer orbitals) have oscillatory surfaces over the full angular range and a plethora of minima; most are too shallow or poorly isolated to support bound states, but a subset at collinear and non-collinear geometries does. The even minima can be understood by overlaying four derivative overlap surfaces, with minima typically at maxima of at least one orbital gradient, often coincident with a node of the underlying wave function. The paper emphasizes that the order of vibrational spacings in the odd trimer differs from both low-l triatomic Rydberg molecules and trilobite trimers.","pith_inferences":["A direct test of the parity classification would be to measure the electronic wave function's nodal structure via Stark or field-ionization imaging; the paper does not propose an experiment.","The building principle—minima track gradient maxima of dimer orbitals—suggests that p-wave Rydberg aggregates beyond trimers could be designed by placing additional atoms at steep-gradient positions of the cluster's own electronic wave function, though the paper stops at trimers.","The reported 400 MHz, 1.2 GHz, and 2.5 GHz ladders imply that microwave double-resonance spectroscopy on a cold Rb gas near n=30 might resolve these states; the paper does not estimate signal strengths or lifetimes.","The neglect of spin structure means quantitative agreement with experiment may require the full spin-coupled treatment the paper lists as future work; the qualitative geometry and parity predictions are likely robust."],"forward_implications":["Odd butterfly trimers should be observable as collinear, equal-bond-length states near n=30, with a predicted ground state at R1=R2 approximately 316 a0 and a nearby symmetric-stretch partner at approximately 340 a0.","The distinct vibrational hierarchy (bending spacing near 2.5 GHz exceeding symmetric stretch near 1.2 GHz) gives a spectroscopic fingerprint separating butterfly trimers from trilobite and low-l trimers.","Even butterfly trimers offer multiple equilibrium geometries, including R1 not equal to R2 and non-collinear arrangements; only the deepest, best-isolated wells should support long-lived vibrational states.","Because many even-trimer equilibria mix R- and theta-butterfly orbitals with very different dipole moments, those states should have large, geometry-dependent electric dipole moments and respond to weak external fields."],"supporting_citations":[{"why":"Establishes the butterfly dimer potential from the p-wave resonance, providing the single-atom binding mechanism the trimer extends.","marker":"[11]"},{"why":"Gives the butterfly dimer electronic structure and scattering-volume description used for the single-atom case.","marker":"[13]"},{"why":"Introduces the hybrid dimer-orbital basis and overlap formalism the trimer calculation builds on.","marker":"[19]"},{"why":"Provides the triatomic Rydberg molecule framework and vibrational Hamiltonian for low-l states.","marker":"[20]"},{"why":"Provides the trilobite-trimer treatment and the bosonic vibrational method whose results are extended here.","marker":"[23]"},{"why":"Supplies the multi-manifold hybrid basis including quantum-defect states used for convergence of butterfly surfaces.","marker":"[30]"},{"why":"Defines the Fermi-Omont pseudopotential with arbitrary partial waves underlying the Hamiltonian.","marker":"[25]"},{"why":"Reports experimental observation of butterfly molecules and their pendular response, motivating the trimer extension.","marker":"[14]"},{"why":"Provides the energy-dependent p-wave scattering volume that sets the allowed bond-length range.","marker":"[10]"}],"fun_headline_variants":["Butterfly trimers: odd smooth, even rugged, parity decides","Odd butterfly trimers collapse, even ones offer rare traps","Two atoms bind to a Rydberg core, forming butterfly trimers","Building principle for butterfly trimers from dimer orbitals","Smooth odd, rugged even: the two faces of butterfly trimers"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that including only the n=29, 30, and 31 Rydberg manifolds with quantum defects for s, p, and d waves gives converged potential energy surfaces; the paper states this without presenting a convergence study, so a shift or loss of minima at higher n would undo the quantitative spectra.","fun_headline_variants_meta":{"raw":{"variants":["Butterfly trimers: odd smooth, even rugged, parity decides","Odd butterfly trimers collapse, even ones offer rare traps","Two atoms bind to a Rydberg core, forming butterfly trimers","Building principle for butterfly trimers from dimer orbitals","Smooth odd, rugged even: the two faces of butterfly trimers"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000442,"raw_usage":{"total_tokens":2245,"prompt_tokens":959,"completion_tokens":1286,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":575,"completion_tokens_details":{"reasoning_tokens":1198}},"tokens_in":575,"tokens_out":1286,"duration_ms":13523,"temperature":1.0,"reasoning_tokens":1198,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:59:01.225829+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same potential surfaces with n up to 32–35 (or with f-wave quantum defects included): if the deep odd-trimer minima near R1=R2 approximately 316 a0 and 340 a0 move by more than a few GHz or disappear, the truncation assumption fails. Alternatively, a microwave spectrum of Rb n approximately 30 dimers plus a second atom that shows no vibrational resonances near the predicted 400 MHz, 1.2 GHz, and 2.5 GHz spacings would refute the quantitative claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the butterfly dimer potential from the p-wave resonance, providing the single-atom binding mechanism the trimer extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the butterfly dimer electronic structure and scattering-volume description used for the single-atom case."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the hybrid dimer-orbital basis and overlap formalism the trimer calculation builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the triatomic Rydberg molecule framework and vibrational Hamiltonian for low-l states."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the trilobite-trimer treatment and the bosonic vibrational method whose results are extended here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the multi-manifold hybrid basis including quantum-defect states used for convergence of butterfly surfaces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Fermi-Omont pseudopotential with arbitrary partial waves underlying the Hamiltonian."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports experimental observation of butterfly molecules and their pendular response, motivating the trimer extension."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the energy-dependent p-wave scattering volume that sets the allowed bond-length range."}],"review_version":1}