{"id":"93dea55a-7e24-4ae7-9de0-3c19794c0ff2","arxiv_id":"2506.19618","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A full 3D multichannel calculation of two dipolar atoms in a cigar trap shows that standard quasi-1D pseudopotentials miss short-range corrections, including a residual interaction at the magic angle and energy-dependent shifts of confinement-induced resonances.","lead":"Ultracold dipolar atoms squeezed into a cigar-shaped trap do not behave like a simple one-dimensional gas: a full three-dimensional calculation shows that the standard 1D model misses short-range corrections. These corrections shift confinement-induced resonances and leave a residual interaction even at the 'magic angle' where dipolar forces supposedly vanish.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 1D effective potential's renormalization is asserted but never demonstrated; because the bare adiabatic potential deepens with channel number, the '3D correction' and the fitted agreement in Fig. 4 are not uniquely defined without a convergence check.","rationale":"The paper's qualitative message—that the standard Gaussian transverse-mode pseudopotential misses short-distance coupling to excited transverse modes—is well motivated and supported by the adiabatic-potential plots. However, the quantitative claim that a regularized effective 1D potential reproduces the 3D spectrum depends on a renormalization that is only named, not derived or tested. The paper explicitly notes that the lowest adiabatic potential deepens with channel number and asserts that the bound-state energy is the physical parameter, but gives neither the renormalization equation nor a convergence study. This is the same weakest point the Reader identified. In addition, the agreement shown in Fig. 4 is obtained by fitting the modified resonance position and strength to the same 3D energies used for validation, so it cannot independently establish accuracy. The single-channel QDT treatment is an acknowledged simplification and is less central to the claim. The recommended conditional verdict is unchanged; the requested convergence and independent-validation tests would turn the conditional into an accept or reject.","tokens_in":11533,"tokens_out":10738,"duration_ms":113818,"concrete_test":"For θF=0 and a fixed physical scattering length (e.g., near the resonance in Fig. 2), recompute the 3D spectrum with N=50, 100, 200, and 400 transverse channels. For each N, build the lowest adiabatic potential, renormalize the contact term so the 1D model reproduces one reference bound-state energy from the 3D calculation, and then predict the next five energy levels and the g1D(a3D) curve. If those predictions shift by more than 5% of the level spacing near the CIR as N doubles, the effective potential is not unique. As an independent check, solve the same 3D Hamiltonian with a cylindrical-coordinate B-spline solver and compare with the 1D model using g1D fixed by that reference state, without refitting to the spectrum.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is that a regularized effective 1D potential reproduces the full 3D spectrum (Sec. IV). The construction in Sec. III diagonalizes Vdd in a truncated transverse harmonic basis. As the paper states, the lowest adiabatic potential deepens as the number of channels is increased, because a Dirac delta appears with the same coefficient in every matrix element; for a matrix with uniform off-diagonal contact terms the lowest eigenvalue scales with the channel number N. The paper adds that 'the physical parameter of this model is the bound state energy,' but no renormalization procedure or N-convergence study is supplied. Without this, the '3D correction' is cutoff-dependent and the extracted g1D(a3D) mapping is not a well-defined object. The validation in Fig. 4 is also partly circular: the same 3D spectrum is used to fix the modified resonance position and strength in Eq. (13). This matters most near CIRs and at the magic angle, where the paper predicts a quantitative failure of Eq. (10).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the dimensional reduction of two-body dipolar interactions in a quasi-one-dimensional cigar-shaped harmonic trap. The authors solve the full 3D two-body problem with van der Waals and dipole-dipole potentials using multichannel Numerov propagation with QDT boundary conditions, then compare the resulting spectrum with effective 1D models. They show that the standard single-transverse-mode Gaussian reduction, Eq. (10), fails at short distances, that the lowest adiabatic potential in a truncated harmonic-oscillator basis deviates from the Gaussian result, and that the mapping between the 3D scattering length and the 1D coupling constant is modified and depends on the specific Feshbach resonance. They also report residual interaction at the magic angle, where the standard 1D dipolar term vanishes.","tokens_in":11774,"tokens_out":3217,"duration_ms":36166,"significance":"If the quantitative claims hold, the paper would provide an important correction to the widely used 1D pseudopotential for dipolar gases, with direct consequences for interpreting experiments on excited states, confinement-induced resonances, and many-body phases in reduced dimension. The numerical approach, combining multichannel calculations with QDT boundary conditions, is credible, and the explicit wavefunction cross-sections in Fig. 3 give physical support to the qualitative conclusion that transverse-mode coupling matters at short distances. The paper also honestly identifies nonuniversal, branch-dependent behavior in the 3D-to-1D mapping. However, the quantitative validation is weakened by the absence of convergence tests for the channel-truncated potential and by a validation procedure that fits the 1D model to the same 3D spectrum it then compares against. These issues must be addressed before the central claim of a quantitatively accurate effective 1D potential can be accepted.","major_comments":[{"comment":"The lowest adiabatic potential is explicitly channel-number dependent: the text states that 'the lowest adiabatic potential gets deeper with the increasing number of included channels.' No renormalization procedure or N-convergence study is supplied. Since the effective 1D potential is derived from this channel-truncated potential, the claimed '3D correction' is cutoff-dependent and the fitted g1D(a3D) mapping is not a uniquely defined object unless a well-defined limit or renormalization condition is established.","section":"Section III, Eq. (12) and Fig. 1F"},{"comment":"The validation is partly circular. The full 3D spectrum is used to fix the modified resonance position and strength in Eq. (13), and the same spectrum is then used to demonstrate agreement with the corrected 1D model. This does not independently establish that the effective 1D model 'accurately reproduces' 3D results. Please provide an out-of-sample test, for example predicting a spectrum at a different transverse confinement, a different dipole angle, or a different energy window without refitting the short-range parameters.","section":"Section IV, Fig. 4"},{"comment":"The statement that 'each branch of the energy spectrum leads to a slightly different fit' means that the extracted g1D is not a single parameter of the model but depends on the branch/energy. The paper does not report the fit residuals per branch or a quantitative criterion for how many branches the 1D model must match. Without such numbers, the visual agreement in Fig. 4 cannot be assessed and the energy-dependent correction remains a fitting statement rather than a predictive result.","section":"Section III, Fig. 4 (lower) and text after Eq. (13)"},{"comment":"The mapping in Eq. (13) is labeled g(1D)_dd even though the text immediately says 'In the absence of dipolar interactions (g(3D)_dd = 0)' and the formula is the standard Olshanii result for the contact coupling g1D. This notation error obscures the role of the dipolar contribution and should be corrected, since Eq. (13) is central to the claimed modified mapping.","section":"Section II.A, Eq. (13)"}],"minor_comments":[{"comment":"The caption says 'ω∥ = 118kHz and ω⊥ = 10ω⊥'; the second relation should presumably read ω⊥ = 10ω∥ to match the text in Section II.A.","section":"Fig. 1 caption"},{"comment":"The text refers to 'Fig. 3A' and 'Fig. 3B' when discussing agreement between 1D and 3D spectra and branch-dependent fits, but Fig. 3 shows wavefunction cross sections. The spectral comparisons appear in Fig. 4; please correct the cross-references.","section":"Section III and Fig. 3"},{"comment":"Equation (10) is split across two displayed lines with a dangling '∝' and the argument of the exponential/erfc is not fully specified. Please define u = x/l⊥ and state the precise regularization of the δ-function term.","section":"Eq. (10)"},{"comment":"The conclusion that the corrected 1D model works 'outside of confinement-induced resonance (CIR) regimes' is stated without a precise definition of what 'outside' means. Please specify the range of |a3D|/l⊥ or the proximity to the CIR where the model is expected to be quantitatively accurate.","section":"Section IV"},{"comment":"The paper mentions 'Dysprosium at θF = 0' in Fig. 1E and 'Erbium particles' in Fig. 3, but does not clearly state the atomic species and parameters used for each figure. Please add explicit parameter tables or state in the captions which species, trap frequencies, and dipole strengths are used.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a relevant and timely problem, and the qualitative message is likely correct: single-mode Gaussian reduction is unreliable at short distances for dipolar interactions. The main barrier is that the quantitative effective 1D potential is not uniquely defined without a channel-convergence or renormalization analysis, and the central validation is self-referential. The lack of any convergence test or residual statistics is surprising for a numerical methods paper and should be addressed before publication. I would not recommend reject at this stage, but the required additions are substantive and go beyond cosmetic revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nYou should know two things about this one. First, it does a real multichannel 3D calculation with a realistic van der Waals plus dipolar potential, QDT boundary conditions, and uses it to test the standard single-transverse-mode 1D pseudopotential. Second, the main quantitative claim—that a renormalized 1D model reproduces the full 3D spectrum—rests partly on a fit to the same spectrum, and the renormalization itself is never shown to converge.\n\nWhat's new: prior work on quasi-1D dipolar scattering used model potentials or generalized pseudopotentials. This paper does the full few-body numerics for Er-like and Dy-like parameters and finds two things worth remembering. The Gaussian reduction fails at distances below about 0.5 l⊥, where coupling to excited transverse modes is strong, and the resulting transverse wavefunctions are not harmonic ground states (Fig. 3). Even at the magic angle, where the effective dipolar potential vanishes in the standard formula, there is a residual short-range interaction decaying as x^-6 that shifts a1D. That is a concrete, testable prediction for experiments that use the magic angle as a contact-only reference.\n\nThe soft spots are where the paper moves from qualitative to quantitative. The adiabatic potential is built by diagonalizing the dipole interaction in a truncated transverse basis. The text states plainly that the lowest adiabatic potential deepens as the number of channels is increased, because a delta-function contact appears with the same coefficient in every matrix element. For a matrix with uniform off-diagonal contact terms, the lowest eigenvalue scales with the number of channels. The paper then says \"the physical parameter of this model is the bound state energy,\" but there is no renormalization calculation and no N-convergence study. Without that, the \"3D correction\" in Fig. 4 is cutoff-dependent, and the extracted g1D(a3D) mapping is not yet a well-defined object. On top of that, the validation is partly circular: the modified resonance position and strength in Eq. (13) are fitted to the same 3D spectrum that the corrected 1D model is then compared against. The authors acknowledge the branch-dependence of the fit but do not address the circularity directly. No code or data is shipped, which makes independent verification harder, though that alone is not disqualifying.\n\nThese problems are real but not fatal. The qualitative conclusions—single-mode Gaussian is quantitatively wrong near confinement-induced resonances and at the magic angle—are supported by the structure of the adiabatic potentials and wavefunction cross-sections, independently of the fitting. What's missing is a demonstration that the effective potential stabilizes as channels are added, and a clearer statement of what is fitted versus predicted. A serious referee should ask for a convergence test and a cleaner separation between fit and validation. The paper deserves referee time; it's an important correction to a widely used approximation in dipolar quantum gas experiments.","headline":"Realistic 3D multichannel calculation shows the standard 1D dipolar pseudopotential is quantitatively unreliable at short distances, but the paper's quantitative mapping rests on an unproven renormalization and a partly circular validation.","tokens_in":12298,"tokens_out":4097,"would_cite":true,"duration_ms":37235,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["67.85.-d","34.50.-s"],"model":"deepseek-v4-flash","headline":"The standard one-dimensional pseudopotential for dipolar gases in tight traps is quantitatively wrong; the paper derives a regularized 3D-corrected effective potential that matches full three-dimensional calculations.","keywords":["quasi-one-dimensional dipolar gas","effective 1D interaction","confinement-induced resonance","dipolar pseudopotential","dimensional reduction","quantum defect theory","magic angle","ultracold dipolar atoms"],"falsifier":"Measure the two-body energy spectrum or confinement-induced resonance position of dysprosium or erbium atoms in a cigar trap with dipole angle tuned through the magic angle, for a fixed 3D scattering length. If the spectrum follows the Gaussian-integrated pseudopotential of Eq. (10), with a vanishing dipolar correction at the magic angle and an unshifted Olshanii resonance, the paper's central claim is wrong; agreement with the renormalized potential would support it.","tokens_in":11348,"feed_emoji":"🧲","tokens_out":5040,"duration_ms":48487,"temperature":0.7,"pith_summary":"This paper argues that the widely used one-dimensional pseudopotential for dipolar atoms in cigar-shaped traps, obtained by integrating out the transverse harmonic-oscillator ground state, misses real three-dimensional physics. Solving the full 3D two-body problem with realistic van der Waals and dipolar potentials shows that coupling to excited transverse modes deepens the effective interaction at short distances and shifts confinement-induced resonances. The authors construct a regularized effective 1D potential, renormalized by the bound-state energy, that reproduces the full 3D spectrum. If correct, every experiment that uses the magic angle as a purely contact-interaction reference, and every many-body calculation built on the old pseudopotential, needs a species- and resonance-specific 3D-to-1D mapping.","feed_headline":"Standard 1D dipolar interaction is wrong in tight traps","feed_subtitle":"Full 3D scattering shows coupling to excited transverse modes reshapes interactions, even at the magic angle.","key_machinery":"The load-bearing object is the renormalized effective one-dimensional potential built from the lowest adiabatic curve obtained by diagonalizing the dipole-dipole interaction in the basis of transverse harmonic-oscillator modes, with the short-range part fixed by the molecular bound-state energy rather than by the bare 3D parameters. This potential is then used with the exact 3D-to-1D mapping for the coupling constant, which follows the Olshanii formula $g^{(1D)} \\propto a_{3D}/(1 - C a_{3D}/l_\\perp)$ but with modified resonance position and strength and with energy-dependent corrections. The machinery also includes multichannel quantum-defect theory with a short-range phase and Numerov integration to obtain the trapped 3D spectrum that the 1D model must reproduce.","core_discovery":"The central claim is that the effective interaction in a quasi-one-dimensional dipolar gas is determined by the full 3D scattering problem, not by projecting onto the lowest transverse mode. The paper shows that the regularized dipole-dipole potential obtained by integrating out the transverse ground state (Eq. 10) deviates from the true lowest adiabatic potential already at distances around half the transverse oscillator length, and that the deviation grows as more channels are included. The physical input is the molecular bound-state energy; with that input, a renormalized 1D potential matches the 3D energy levels. The paper further finds that the 3D-to-1D mapping between the 3D scattering length and the 1D coupling constant follows the Olshanii form but with shifted resonance position and strength, with distinct branches for s-wave- and d-wave-dominated resonances, and that residual dipolar corrections persist even at the magic angle where the old pseudopotential predicts no dipolar interaction at all.","pith_inferences":["The energy-dependence of the fitted $g_{1D}$ suggests that the correct low-energy description of q1D dipolar gases may be a two-channel or frequency-dependent effective interaction, rather than a single static potential.","Since the physical parameter is the bound-state energy, the renormalized potential should transfer across species and trap geometries once that energy is matched, offering a testable prediction for molecular dipolar gases now becoming available.","The residual $x^{-6}$ tail at the magic angle, though short-ranged and weak, may still matter for many-body phases sensitive to small scattering-length shifts, such as droplet stability boundaries.","A direct extension would be to compute the same corrected potential for tilted fields $\\theta_F$ between 0 and $\\pi/2$ and to tabulate the resonance-shift function, which would let experimental groups interpolate rather than rerun 3D calculations."],"forward_implications":["Magic-angle experiments cannot be treated as purely contact-interacting reference systems; a residual short-range dipolar potential shifts the 1D scattering length.","Confinement-induced resonances shift relative to the single-mode prediction, so locating them experimentally requires solving the 3D problem or using the corrected mapping.","Each Feshbach resonance used to tune $a_{3D}$ produces a different 3D-to-1D mapping, meaning 1D parameters are nonuniversal across resonances.","Excited-state and quench experiments access energies where the effective range of the interaction matters; a constant $g_{1D}$ fitted to low-lying levels will fail.","Quantum-droplet and extended-Hubbard-model simulations that take Eq. (10) as input need the renormalized potential to get phase boundaries right."],"supporting_citations":[{"why":"Supplies the analytic 3D-to-1D coupling mapping $g^{(1D)} \\propto a_{3D}/(1 - C a_{3D}/l_\\perp)$ that the paper generalizes.","marker":"[57]"},{"why":"Source of the regularized single-mode dipolar potential (Eq. 10) that the paper shows is insufficient.","marker":"[60]"},{"why":"Erratum to the same, providing the corrected form of the regularized potential used as baseline.","marker":"[61]"},{"why":"Previous q1D dipolar treatment assuming ground-state transverse Gaussians, the approach the paper upgrades.","marker":"[64]"},{"why":"Earlier study of repulsive dipole orientation showing coupling to higher transverse modes matters.","marker":"[67]"},{"why":"Dipolar confinement-induced resonances in waveguides, context for the CIR shifts studied here.","marker":"[65]"},{"why":"Van der Waals quantum-defect solutions for the short-range phase used in the boundary conditions.","marker":"[73]"},{"why":"Quasi-universal dipolar scattering, basis for the claim of nonuniversal resonance-dependent branches.","marker":"[71]"}],"fun_headline_variants":["3D scattering reshapes dipolar interactions in 1D traps","Tight traps expose flaws in 1D dipolar gas model","Magic angle still feels dipolar forces: 3D matters","1D pseudopotential fails for dipolar gases in tight traps","Full 3D treatment needed for dipolar gas interactions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the claim that truncating the adiabatic channel basis and renormalizing by the bound-state energy gives a unique effective 1D potential; if different channel sets or different bound states yield different corrected potentials, the '3D correction' is not a well-defined quantity.","fun_headline_variants_meta":{"raw":{"variants":["3D scattering reshapes dipolar interactions in 1D traps","Tight traps expose flaws in 1D dipolar gas model","Magic angle still feels dipolar forces: 3D matters","1D pseudopotential fails for dipolar gases in tight traps","Full 3D treatment needed for dipolar gas interactions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000568,"raw_usage":{"total_tokens":2631,"prompt_tokens":825,"completion_tokens":1806,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":441,"completion_tokens_details":{"reasoning_tokens":1732}},"tokens_in":441,"tokens_out":1806,"duration_ms":12707,"temperature":1.0,"reasoning_tokens":1732,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T23:06:40.163980+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the two-body energy spectrum or confinement-induced resonance position of dysprosium or erbium atoms in a cigar trap with dipole angle tuned through the magic angle, for a fixed 3D scattering length. If the spectrum follows the Gaussian-integrated pseudopotential of Eq. (10), with a vanishing dipolar correction at the magic angle and an unshifted Olshanii resonance, the paper's central claim is wrong; agreement with the renormalized potential would support it.","supporting_citations":[{"cited_title":"Deuretzbacher, J","cited_arxiv_id":null,"evidence_quote":"Source of the regularized single-mode dipolar potential (Eq. 10) that the paper shows is insufficient."},{"cited_title":"Deuretzbacher, J","cited_arxiv_id":null,"evidence_quote":"Erratum to the same, providing the corrected form of the regularized potential used as baseline."},{"cited_title":"Sinha and L","cited_arxiv_id":null,"evidence_quote":"Previous q1D dipolar treatment assuming ground-state transverse Gaussians, the approach the paper upgrades."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier study of repulsive dipole orientation showing coupling to higher transverse modes matters."},{"cited_title":"Giannakeas, V","cited_arxiv_id":null,"evidence_quote":"Dipolar confinement-induced resonances in waveguides, context for the CIR shifts studied here."},{"cited_title":"Gao, Solutions of the schrödinger equation for an at- tractive 1/r6 potential, Phys","cited_arxiv_id":null,"evidence_quote":"Van der Waals quantum-defect solutions for the short-range phase used in the boundary conditions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Quasi-universal dipolar scattering, basis for the claim of nonuniversal resonance-dependent branches."}],"review_version":1}