Pith. sign in

REVIEW 4 major objections 4 minor 72 references

This paper claims that a regularized Born series solves coupled-channel quantum scattering in O(N²) operations with guaranteed convergence.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-03 12:53 UTC pith:YVJBTTKW

load-bearing objection A credible O(N^2) single-column coupled-channel method with genuinely useful closed-channel handling and a convincing He+CO benchmark; the 'guaranteed convergence' label currently outruns the spectral verification in the paper. the 4 major comments →

arxiv 2601.01159 v2 pith:YVJBTTKW submitted 2026-01-03 physics.chem-ph physics.atom-phphysics.comp-ph

A quadratic-scaling algorithm with guaranteed convergence for quantum coupled-channel calculations

Hubert J. J\'o\'zwiak , Md Muktadir Rahman , Timur V. Tscherbul This is my paper
classification physics.chem-ph physics.atom-phphysics.comp-ph PACS 34.50.-s34.50.Ez
keywords coupled-channel scatteringLippmann-Schwinger equationBorn seriesWeinberg eigenvaluesspectral regularizationquadratic scalingclosed channelsFeshbach resonances
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper aims to break the cubic-scaling barrier of standard coupled-channel scattering calculations by solving the Lippmann-Schwinger equation for a single column of the S-matrix instead of propagating a full wavefunction matrix. Its WISE algorithm identifies the Weinberg eigenvalues of the scattering kernel that lie outside the unit circle, subtracts their subspace, and iterates the remaining regularized kernel, whose spectral radius is below one. This makes closed channels and Feshbach resonances tractable and reduces the cost from O(N³) to O(N²) per calculation. The authors test the method on He+CO and CO+N₂ collisions, reporting agreement with benchmark coupled-channel results and a measured quadratic scaling of total CPU time.

Core claim

The paper's central claim is that the multichannel Lippmann-Schchinger equation for one column of the S-matrix can be solved in O(N²) operations with guaranteed convergence, regardless of interaction strength or anisotropy. The divergence of the naive Born series is traced to the presence of Weinberg eigenvalues with magnitude at least one. The WISE algorithm splits the kernel into a divergent part built from those eigenvectors and a regularized remainder with spectral radius strictly below one, solves the divergent part exactly through a small linear system, and iterates the convergent series for the rest. Closed channels are incorporated through a stable ratio-propagation construction of t

What carries the argument

The load-bearing object is the spectrum of the discretized multichannel Lippmann-Schwinger kernel K = G U_off, and in particular its Weinberg eigenvalues η_k, defined as eigenvalues of that kernel. The WISE decomposition K = K_R + K_D splits off eigenvectors with |η_k| ≥ 1 into a separable divergent operator K_D, leaving a regularized remainder K_R whose spectral radius is less than one. The full solution is then the convergent Born series for K_R plus a finite correction involving the divergent eigenvectors, with coefficients fixed by a small linear system. A stable ratio-propagation method constructs the Green's function without explicit evaluation of exponentially growing closed-channel w

Load-bearing premise

The guarantee stands only if the eigensolver finds every Weinberg eigenvalue with magnitude at least one and the kernel is diagonalizable on that divergent subspace; if an eigenvalue is missed or defective, the regularized kernel may still have spectral radius at least one and the series can diverge, and the numerical agreement additionally assumes the uniform radial grid is fine enough to represent the Green's operator faithfully.

What would settle it

For a system with a closed channel, compute the full spectrum of the regularized kernel K_R after regularization: if any eigenvalue has magnitude at least one, the claimed guaranteed convergence fails. On the numerical side, rerun the CO+N₂ benchmark with a halved radial step, e.g., Δr = 0.05 a₀, and check that the WISE cross-sections move onto the reference coupled-channel curve; if they do not, the discretized kernel, not the regularization, is controlling the error.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • State-to-state scattering for molecules prepared in a specific incoming quantum state can be computed in O(N²) time, a quadratic reduction from the O(N³) cost of standard propagation methods.
  • Closed channels can be included robustly, so Feshbach resonances and threshold effects are captured rather than lost in open-channel-only approximations.
  • Highly anisotropic molecule-molecule systems such as CO+N₂ become feasible for systematic energy-resolved calculations, including the scans needed for atmospheric line-shape modeling.
  • The regularization supplies a formal convergence guarantee for the iterative series, removing the need for ad-hoc acceleration schemes in a broad class of coupled-channel problems.
  • Weinberg eigenvalue spectra provide a quantitative diagnostic of how strongly perturbative an interaction potential is, which can guide basis-set and method choices.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • A cheap, auditable robustness check is left implicit: after forming K_R, one could directly verify that its spectral radius is below one by testing K_R v_k ≈ 0 for each removed eigenvector; adding this residual check would turn the guarantee from an assumption into a numerical fact.
  • Because the benchmarks show the eigenvalue search dominating total cost, future gains for very large channel counts likely depend on improving the eigensolver step more than on the Born iterations themselves.
  • The residual differences visible in the CO+N₂ results could plausibly be removed by adaptive radial grids at no extra asymptotic cost, which the paper notes but does not implement.
  • The Weinberg spectrum could be repurposed as a principled channel-selection tool: channels that push eigenvalues outside the unit circle are precisely the ones responsible for strong coupling, so they indicate which closed channels matter most.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The manuscript introduces the Weinberg-regularized Iterative Series Expansion (WISE) algorithm for quantum coupled-channel scattering. The method recasts the coupled-channel equations as a multichannel Lippmann-Schwinger equation for a single column of the S-matrix, discretizes the radial coordinate, and solves the resulting linear system by Born iteration. To cure divergence of the Born series, the kernel is decomposed into a divergent part K_D built from Weinberg eigenvalues with |η|≥1 and a regularized part K_R claimed to have spectral radius <1; the divergent part is treated by a small linear system. Closed channels are included through a ratio-propagation construction of the Green's function. The algorithm is validated against standard CC results for He+CO and CO+N2 collisions, and timing benchmarks are reported for channel counts up to 343.

Significance. If the central claims hold, this is a substantial methodological advance: quadratic scaling for a single S-matrix column, explicit inclusion of closed channels, and removal of the Born-series divergence would extend rigorous CC calculations to systems currently considered intractable. The He+CO benchmark is strong and clean: it reproduces independent CC reference results, including Feshbach resonances, without fitted parameters. The CO+N2 application to a molecule-molecule system with strong anisotropy is also a meaningful stress test. The matrix-free implementation and the stable ratio-propagation construction of the Green's function are well-designed engineering. However, the paper's main guarantee is asserted rather than proved, and several numerical details are under-specified; the work is promising but needs revision before the central claims are fully supported.

major comments (4)
  1. [Spectral Regularization and Guaranteed Convergence, Eq. (16)] The statement after Eq. (16) that ρ(K_R)<1 'by construction' is the load-bearing step for the guaranteed-convergence claim. The construction using right and left eigenvectors of K is valid only if the discretized kernel is diagonalizable on the subspace spanned by all eigenvalues with |η|≥1 and if the Arnoldi/ARPACK run returns every such eigenvalue. For a non-symmetric matrix, a defective Jordan block or a missed eigenvalue leaves ρ(K_R)≥1; then the Born series for u0^(R) and v_k^(R) in Eqs. (17)-(18) need not converge, and Eq. (19) need not represent the solution of Eq. (4). Please provide a theorem with explicit assumptions, or add numerical verification: residuals ||K v_k - η_k v_k|| and ||K† w_k - η_k* w_k||, a completeness check for the selected divergent set, and an estimate of ρ(K_R).
  2. [Application, Fig. 3(b), and Methods (grid)] The abstract and text describe the CO+N2 results as 'exact quantum results,' but Fig. 3(b) shows visible residual differences between WISE and CC. The Methods state that Δr=0.1 'ensures strict convergence (<1%) of benchmark CC calculations' and then attribute the WISE residuals to the uniform radial discretization. This conflates two sources of error: the trapezoid quadrature in Eq. (10) is not identical to the propagator used in the benchmark, and the same grid step can give different quadrature errors in the two formulations. Please quantify the residuals and report a convergence study in Δr, or an adaptive quadrature, before using the word 'exact'.
  3. [Methods, Numerical Implementation of the Iterative Algorithm] The fixed-point iteration is stopped when the L2 norm of the update ||u^(n)-u^(n-1)||_2 falls below 0.1. This is an absolute tolerance, not a relative one; for wavefunction components of order unity it is a very loose stopping criterion. Since the paper claims numerically exact agreement with CC, please justify this tolerance or perform a convergence test with ε decreased by factors of 10 and report the effect on the S-matrix elements.
  4. [Scaling and Summary/Outlook] The quadratic-scaling claim is supported by timing data only up to N=343 and by visual fits in Fig. 4. The total time includes an Arnoldi eigensolver whose iteration count is not bounded as N grows, and the Outlook estimates n_D ∝ N (up to 10^5), in which case a dense solve of the n_D×n_D M-matrix would be O(n_D^3). To support 'quadratic scaling guaranteed,' please give a complexity argument for the eigensolver and regularization steps, or report scaling of each stage with N over a wider range.
minor comments (4)
  1. [Abstract] The abstract as supplied mentions a 'contour-integral construction' that does not appear anywhere in the body; the convergence guarantee is instead based on spectral regularization. Please reconcile the abstract with the actual method.
  2. [Summary of the Algorithm] The statement that K_R v_k ≈0 should be an equality if K_D is constructed from the complete divergent eigenspace; the approximate equality signals an unquantified Arnoldi error. Please explain or remove the approximation.
  3. [Methods, Numerical Implementation] The M-matrix in Eq. (21) is described as sparse, but the implementation uses the dense LAPACK driver ZGESV. Clarify when a sparse solver would be needed and how the sparse structure would be exploited.
  4. [Methods, ARPACK] Please report ARPACK restart parameters, the number of eigenvalues requested, and residual tolerances so that the calculation is reproducible and the claimed convergence of the divergent subspace can be assessed.

Circularity Check

0 steps flagged

No significant circularity: the central derivation is benchmarked against independent CC results, and the self-citations are context, not load-bearing premises.

full rationale

The paper's derivation chain is not circular. The WISE algorithm solves the multichannel Lippmann-Schrödinger equation by Born iteration after spectrally removing the Weinberg subspace with |η_k| ≥ 1; the statement that ρ(K_R) < 1 'by construction' is a legitimate definitional construction (the divergent modes are removed), not a fitted parameter renamed as a prediction. The benchmark S-matrix elements and cross-sections are compared with independent standard CC calculations on published ab initio PESs, and no parameter is tuned to force agreement with the target observables. The authors' self-citations (refs. 14, 16, 28, 43, 62, 63) supply potential energy surfaces, line-shape context, and an estimate of conventional CC capabilities; none is invoked as the unverified basis for the convergence guarantee. The main weakness—that the 'guaranteed convergence' claim assumes the Arnoldi iterations find every eigenvalue with |η| ≥ 1 and that K is diagonalizable on that subspace—is an unproven mathematical/computational premise and therefore a correctness risk, not circularity, because the premise is not equivalent to the conclusions it is used to support.

Axiom & Free-Parameter Ledger

4 free parameters · 6 axioms · 0 invented entities

The method uses external PESs and standard scattering theory; no target-fitting parameters. The main unproven input is the spectral assumption on K needed for the regularization guarantee, plus numerical parameters (grid, tolerance, basis) that affect the claimed agreement. No invented physical entities are introduced.

free parameters (4)
  • Born-series convergence tolerance ε = 0.1
    Chosen by hand to stop fixed-point iteration based on the L2 norm of the update; affects reported accuracy but is not fitted to target observables.
  • Radial grid step Δr and grid size = 0.05 a0 / Ng=933 (He+CO); 0.1 a0 / Ng=458 (CO+N2)
    Chosen for the initial implementation; the CO+N2 value is explicitly blamed for residual differences in Fig. 3(b).
  • Basis truncation for CO+N2 = jCO=0-14, jN2=4,6,8
    Chosen to converge the target cross-section to within 20%; a numerical convergence parameter rather than a fitted constant.
  • Arnoldi subspace/restart parameters = not specified
    The number of Arnoldi vectors and restart settings are not given, yet the convergence guarantee depends on capturing all divergent eigenvalues.
axioms (6)
  • standard math The multichannel LS equation (4) with Green's function (9) is equivalent to the CC equations (3) with the proper boundary conditions.
    Invoked at Eqs. (4)-(9) as standard scattering theory; needed for the entire iterative formulation.
  • domain assumption The non-symmetric kernel K has a complete set of left and right eigenvectors for its divergent subspace.
    Needed for the spectral split in Eqs. (16)-(19); no diagonalizability or completeness theorem is provided.
  • ad hoc to paper After removing all eigenvalues with |η|≥1, the regularized kernel K_R has spectral radius <1, and Arnoldi iteration finds every such eigenvalue.
    The 'guaranteed convergence' claim rests on this; no residual check or error bound is supplied to verify that no divergent eigenvalue was missed.
  • domain assumption The closed-channel Green's function can be constructed stably from solution ratios via Eqs. (13)-(15).
    The ratio-propagation formulas are derived from asymptotic Bessel forms but are not proven to be uniformly stable for all closed channels.
  • domain assumption The ab initio PESs from refs. [71] and [54] are accurate enough for benchmark comparisons.
    External inputs from prior literature; the benchmarking value of the paper depends on them.
  • domain assumption The rigid-rotor approximation with frozen bond lengths for CO and N2.
    Used to define the internal Hamiltonian H_int in the Methods; standard for this type of scattering calculation.

pith-pipeline@v1.3.0-alltime-deepseek · 18850 in / 20811 out tokens · 206153 ms · 2026-08-03T12:53:40.417459+00:00 · methodology

0 comments
read the original abstract

Rigorous quantum dynamics calculations provide essential insights into complex scattering phenomena across atomic and molecular physics, chemical reaction dynamics, and astrochemistry. However, the application of the gold-standard quantum coupled-channel (CC) method has been fundamentally constrained by a steep cubic scaling of computational cost $[{O}(N^3)]$. Here, we develop a general, rigorous, and robust method for solving the time-independent Schr\"odinger equation for a single column of the scattering S-matrix with quadratic scaling $[{O}(N^2)]$ in the number of channels. The Weinberg-regularized Iterative Series Expansion (WISE) algorithm resolves the divergence issues affecting iterative techniques by applying a regularization procedure to the kernel of the multichannel Lippmann-Schwinger integral equation. The method also explicitly incorporates closed-channel effects, including those responsible for multichannel Feshbach resonances. We demonstrate the power of this approach by performing rigorous calculations on He + CO and CO + N$_2$ collisions, achieving exact quantum results with quadratic scaling guaranteed by a contour-integral construction. Our results establish a highly scalable computational paradigm, enabling state-to-state quantum scattering computations for complex molecular systems.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

72 extracted references · 40 canonical work pages

  1. [1]

    Phys.82, 1225–1286 (2010), doi:10.1103/RevModPhys.82.1225,https://link.aps.org/ doi/10.1103/RevModPhys.82.1225

    C.Chin,R.Grimm,P.Julienne,E.Tiesinga,Feshbachresonancesinultracoldgases.Rev.Mod. Phys.82, 1225–1286 (2010), doi:10.1103/RevModPhys.82.1225,https://link.aps.org/ doi/10.1103/RevModPhys.82.1225

    doi:10.1103/revmodphys.82.1225 2010

  2. [2]

    C.H.Greene,P.Giannakeas,J.Pérez-Ríos,Universalfew-bodyphysicsandclusterformation. Rev. Mod. Phys.89, 035006 (2017), doi:10.1103/RevModPhys.89.035006,https://link. aps.org/doi/10.1103/RevModPhys.89.035006

    work page doi:10.1103/revmodphys.89.035006 2017

  3. [3]

    R.Krems,MoleculesinElectromagneticFields:FromUltracoldPhysicstoControlledChem- istry(Wiley-VCH) (2019)

    2019

  4. [4]

    S. C. Althorpe, D. C. Clary, Quantum Scattering Calculations on Chemical Reactions.Annu. Rev. Phys. Chem.54(1), 493–529 (2003), doi:10.1146/annurev.physchem.54.011002.103750, http://dx.doi.org/10.1146/annurev.physchem.54.011002.103750

    arXiv 2003

  5. [5]

    D. C. Clary, Quantum Dynamics of Chemical Reactions.Science321(5890), 789 (2008), doi:10.1126/science.1157718

    work page doi:10.1126/science.1157718 2008

  6. [6]

    R. V. Krems, Cold Controlled Chemistry.Phys. Chem. Chem. Phys.10, 4079–4092 (2008)

    2008

  7. [7]

    Balakrishnan, Perspective: Ultracold molecules and the dawn of cold controlled chemistry

    N. Balakrishnan, Perspective: Ultracold molecules and the dawn of cold controlled chemistry. J. Chem. Phys.145(15), 150901 (2016), doi:10.1063/1.4964096

    doi:10.1063/1.4964096 2016

  8. [8]

    Roueff, F

    E. Roueff, F. Lique, Molecular Excitation in the Interstellar Medium: Recent Advances in Collisional, Radiative, and Chemical Processes.Chem. Rev.113(12), 8906–8938 (2013), doi:10.1021/cr400145a,http://dx.doi.org/10.1021/cr400145a

    doi:10.1021/cr400145a 2013

  9. [9]

    E. F. van Dishoeck, E. Herbst, D. A. Neufeld, Interstellar Water Chemistry: From Laboratory to Observations.Chem. Rev.113(12), 9043–9085 (2013), doi:10.1021/cr4003177,https: //doi.org/10.1021/cr4003177

    doi:10.1021/cr4003177 2013

  10. [10]

    Hagino, K

    K. Hagino, K. Ogata, A. M. Moro, Coupled-channels calculations for nuclear reac- tions: From exotic nuclei to superheavy elements.Prog. Part. Nucl. Phys.125, 103951 24 (2022), doi:https://doi.org/10.1016/j.ppnp.2022.103951,https://www.sciencedirect. com/science/article/pii/S0146641022000126

    arXiv 2022

  11. [11]

    J.Hartmann,C.Boulet,D.Robert,CollisionalEffectsonMolecularSpectra(Elsevier)(2008), doi:10.1016/b978-0-444-52017-3.x0001-5

    work page doi:10.1016/b978-0-444-52017-3.x0001-5 2008

  12. [12]

    Hartmann,et al., Recent advances in collisional effects on spectra of molecular gases and their practical consequences.J

    J.-M. Hartmann,et al., Recent advances in collisional effects on spectra of molecular gases and their practical consequences.J. Quant. Spectrosc. Radiat. Transf.213, 178 – 227 (2018), doi:https://doi.org/10.1016/j.jqsrt.2018.03.016

    work page doi:10.1016/j.jqsrt.2018.03.016 2018

  13. [13]

    P. J. Rayer,Pressure Broadening of Spectral Lines: The Theory of Line Shape in Atmospheric Physics(Cambridge University Press) (2020)

    2020

  14. [14]

    H.Jóźwiak,F.Thibault,H.Cybulski,P.Wcisło,AbinitioinvestigationoftheCO–N 2 quantum scattering: The collisional perturbation of the pure rotational R(0) line in CO.J. Chem. Phys. 154(5) (2021), doi:10.1063/5.0040438,http://dx.doi.org/10.1063/5.0040438

    work page doi:10.1063/5.0040438 2021

  15. [15]

    A. Olejnik,et al., Ab initio quantum scattering calculations and a new potential energy sur- face for the HCl(𝑋1Σ+)-O2(𝑋 3Σ− 𝑔) system: Collision-induced line shape parameters for O2- perturbed R(0) 0–0 line in H35Cl.J. Chem. Phys.159(13) (2023), doi:10.1063/5.0169968, http://dx.doi.org/10.1063/5.0169968

    work page doi:10.1063/5.0169968 2023

  16. [16]

    Astron.Astrophys.687,A69(2024),doi:10.1051/0004-6361/202449889,https://doi.org/ 10.1051/0004-6361/202449889

    H.Jóźwiak,etal.,AccuratereferencespectraofHDinanH 2–Hebathforplanetaryapplications. Astron.Astrophys.687,A69(2024),doi:10.1051/0004-6361/202449889,https://doi.org/ 10.1051/0004-6361/202449889

    work page doi:10.1051/0004-6361/202449889 2024

  17. [17]

    Tennyson, S

    J. Tennyson, S. N. Yurchenko, Laboratory spectra of hot molecules: Data needs for hot super- Earth exoplanets.Mol. Astrophys.8, 1–18 (2017), doi:10.1016/j.molap.2017.05.002,http: //dx.doi.org/10.1016/j.molap.2017.05.002

    work page doi:10.1016/j.molap.2017.05.002 2017

  18. [18]

    Niraula,et al., The impending opacity challenge in exoplanet atmospheric characterization

    P. Niraula,et al., The impending opacity challenge in exoplanet atmospheric characterization. Nat. Astron.6(11), 1287–1295 (2022), doi:10.1038/s41550-022-01773-1,http://dx.doi. org/10.1038/s41550-022-01773-1. 25

    doi:10.1038/s41550-022-01773-1 2022

  19. [19]

    K. L. Chubb,et al., Data availability and requirements relevant for the Ariel space mission and other exoplanet atmosphere applications.RAS Techniques and Instruments3(1), 636–690 (2024), doi:10.1093/rasti/rzae039,http://dx.doi.org/10.1093/rasti/rzae039

    doi:10.1093/rasti/rzae039 2024

  20. [20]

    Wiesenfeld,et al., Ab Initio Quantum Dynamics as a Scalable Solution to the Exoplanet Opacity Challenge: A Case Study of CO2 in a Hydrogen Atmosphere.Astrophys

    L. Wiesenfeld,et al., Ab Initio Quantum Dynamics as a Scalable Solution to the Exoplanet Opacity Challenge: A Case Study of CO2 in a Hydrogen Atmosphere.Astrophys. J.981(2), 148 (2025), doi:10.3847/1538-4357/adb02e,http://dx.doi.org/10.3847/1538-4357/ adb02e

    doi:10.3847/1538-4357/adb02e 2025

  21. [21]

    J. L. Bohn, A. M. Rey, J. Ye, Cold molecules: Progress in quantum engineering of chemistry and quantum matter.Science357, 1002–1010 (2017)

    2017

  22. [22]

    S. N. Vogels,et al., Imaging resonances in low-energy NO-He inelastic collisions.Science 350(6262), 787 (2015)

    2015

  23. [23]

    Bause, A

    R. Bause, A. Christianen, A. Schindewolf, I. Bloch, X.-Y. Luo, Ultracold Sticky Collisions: TheoreticalandExperimentalStatus.J.Phys.Chem.A127(3),729–741(2023),doi:10.1021/ acs.jpca.2c08095,http://dx.doi.org/10.1021/acs.jpca.2c08095

    work page doi:10.1021/acs.jpca.2c08095 2023

  24. [24]

    Margulis,et al., Tomography of Feshbach resonance states.Science380(6640), 77–81 (2023), doi:10.1126/science.adf9888,https://doi.org/10.1126/science.adf9888

    B. Margulis,et al., Tomography of Feshbach resonance states.Science380(6640), 77–81 (2023), doi:10.1126/science.adf9888,https://doi.org/10.1126/science.adf9888

    work page doi:10.1126/science.adf9888 2023

  25. [25]

    K.P.Horn,etal.,FeshbachResonancesinColdCollisions:BenchmarkingState-of-the-ArtAb Initio Potential Energy Surfaces.J. Phys. Chem. Lett.16, 7862–7867 (2025), doi:10.1021/acs. jpclett.5c01581,https://pubs.acs.org/doi/abs/10.1021/acs.jpclett.5c01581

    doi:10.1021/acs 2025

  26. [26]

    Tang,et al., Quantum state–resolved molecular dipolar collisions over four decades of energy.Science379(6636),1031–1036(2023),doi:10.1126/science.adf9836,https://doi

    G. Tang,et al., Quantum state–resolved molecular dipolar collisions over four decades of energy.Science379(6636),1031–1036(2023),doi:10.1126/science.adf9836,https://doi. org/10.1126/science.adf9836

    work page doi:10.1126/science.adf9836 2023

  27. [27]

    J. J. Park,et al., Spectrum of Feshbach Resonances inNaLi+NaCollisions.Phys. Rev. X 13, 031018 (2023), doi:10.1103/PhysRevX.13.031018,https://link.aps.org/doi/10. 1103/PhysRevX.13.031018. 26

    work page doi:10.1103/physrevx.13.031018 2023

  28. [28]

    Morita, M

    M. Morita, M. B. Kosicki, P. S. Żuchowski, P. Brumer, T. V. Tscherbul, Magnetic Fesh- bach resonances in ultracold atom-molecule collisions.Phys. Rev. A110, L021301 (2024), doi:10.1103/PhysRevA.110.L021301,https://link.aps.org/doi/10.1103/PhysRevA. 110.L021301

    work page doi:10.1103/physreva.110.l021301 2024

  29. [29]

    Y.-X.Liu,etal.,Hyperfine-to-rotationalenergytransferinultracoldatom–moleculecollisions ofRbandKRb.Nat.Chem.17(5),688–694(2025),doi:10.1038/s41557-025-01778-z,https: //doi.org/10.1038/s41557-025-01778-z

    doi:10.1038/s41557-025-01778-z 2025

  30. [30]

    M. B. Soley, E. J. Heller, Ultracold Molecular Collisions: Quasiclassical, Semiclassical, and Classical Approaches in the Quantum Regime.Chem. Rev.125(14), 6609–6652 (2025), doi: 10.1021/acs.chemrev.4c01014,https://doi.org/10.1021/acs.chemrev.4c01014

    work page doi:10.1021/acs.chemrev.4c01014 2025

  31. [31]

    J. P. D’Incao, Few-body physics in resonantly interacting ultracold quantum gases.J. Phys. B51(4), 043001 (2018), doi:10.1088/1361-6455/aaa116,http://dx.doi.org/10.1088/ 1361-6455/aaa116

    work page doi:10.1088/1361-6455/aaa116 2018

  32. [32]

    Haze,et al., Controlling few-body reaction pathways using a Feshbach resonance.Nat

    S. Haze,et al., Controlling few-body reaction pathways using a Feshbach resonance.Nat. Phys.21(2), 228–232 (2025), doi:10.1038/s41567-024-02726-3,https://doi.org/10. 1038/s41567-024-02726-3

    doi:10.1038/s41567-024-02726-3 2025

  33. [33]

    M. D. Higgins, C. H. Greene, Five-body recombination of identical bosons.Proc. Natl. Acad. Sci.USA122(18),e2503390122(2025),doi:10.1073/pnas.2503390122,https://doi.org/ 10.1073/pnas.2503390122

    work page doi:10.1073/pnas.2503390122 2025

  34. [34]

    Skouteris, J

    D. Skouteris, J. Castillo, D. Manolopoulos, ABC: a quantum reactive scattering program. Comput. Phys. Commun.133(1), 128–135 (2000), doi:10.1016/s0010-4655(00)00167-3, http://dx.doi.org/10.1016/s0010-4655(00)00167-3

    work page doi:10.1016/s0010-4655(00)00167-3 2000

  35. [35]

    Bostan, B

    D. Bostan, B. Mandal, D. Babikov, Mixed quantum/classical theory for rotationally inelastic scattering of identical collision partners revised.Phys. Chem. Chem. Phys.26(43), 27567– 27582 (2024), doi:10.1039/D4CP02806F,http://dx.doi.org/10.1039/D4CP02806F. 27

    work page doi:10.1039/d4cp02806f 2024

  36. [36]

    L. D. Thomas, Solution of the coupled equations of inelastic atom–molecule scattering for a single initial state.J. Chem. Phys.70(6), 2979–2985 (1979), doi:10.1063/1.437836,https: //doi.org/10.1063/1.437836

    work page doi:10.1063/1.437836 1979

  37. [37]

    L. D. Thomas, Solution of the coupled equations of inelastic atom–molecule scattering for a single initial state. II. Use of nondiagonal matrix Green functions.J. Chem. Phys.76(10), 4925–4931 (1982), doi:10.1063/1.442838,https://doi.org/10.1063/1.442838

    work page doi:10.1063/1.442838 1982

  38. [38]

    Rhoades-Brown, M

    M. Rhoades-Brown, M. H. Macfarlane, S. C. Pieper, Techniques for heavy-ion coupled- channels calculations. II. Iterative solution of the coupled radial equations.Phys. Rev. C 21, 2436–2446 (1980), doi:10.1103/PhysRevC.21.2436,https://link.aps.org/doi/10. 1103/PhysRevC.21.2436

    work page doi:10.1103/physrevc.21.2436 1980

  39. [39]

    A. M. Arthurs, A. Dalgarno, The theory of scattering by a rigid rotator.Proc. R. Soc. Lond. Ser. A256(1287), 540–551 (1960), doi:10.1098/rspa.1960.0125,http://dx.doi.org/10. 1098/rspa.1960.0125

    arXiv 1960

  40. [40]

    Johnson, The multichannel log-derivative method for scattering calculations.J

    B. Johnson, The multichannel log-derivative method for scattering calculations.J. Comput. Phys.13(3), 445–449 (1973), doi:10.1016/0021-9991(73)90049-1,http://dx.doi.org/ 10.1016/0021-9991(73)90049-1

    work page doi:10.1016/0021-9991(73)90049-1 1973

  41. [41]

    B. R. Johnson, The renormalized Numerov method applied to calculating bound states of the coupled-channel Schroedinger equation.J. Chem. Phys.69(10), 4678–4688 (1978), doi: 10.1063/1.436421,http://dx.doi.org/10.1063/1.436421

    doi:10.1063/1.436421 1978

  42. [42]

    D. E. Manolopoulos, An improved log derivative method for inelastic scattering.J. Chem. Phys.85(11), 6425–6429 (1986), doi:10.1063/1.451472,http://dx.doi.org/10.1063/ 1.451472

    work page doi:10.1063/1.451472 1986

  43. [43]

    Y. V. Suleimanov, T. V. Tscherbul, R. V. Krems, Efficient method for quantum calculations of molecule-molecule scattering properties in a magnetic field.J. Chem. Phys.137(2) (2012), doi:10.1063/1.4733288,http://dx.doi.org/10.1063/1.4733288

    work page doi:10.1063/1.4733288 2012

  44. [44]

    R. G. Newton,Scattering theory of waves and particles(Springer Science & Business Media) (2013). 28

    2013

  45. [45]

    Child,Molecular Collision Theory(Wiley) (1974)

    M. Child,Molecular Collision Theory(Wiley) (1974)

    1974

  46. [46]

    Theory.J

    R.T.Pack,G.A.Parker,Quantumreactivescatteringinthreedimensionsusinghyperspherical (APH) coordinates. Theory.J. Chem. Phys.87(7), 3888–3921 (1987), doi:10.1063/1.452944, https://doi.org/10.1063/1.452944

    work page doi:10.1063/1.452944 1987

  47. [47]

    D. E. Manolopoulos, R. E. Wyatt, D. C. Clary, Iterative solution in quantum scattering theory. The log derivative Kohn approach.J. Chem. Soc., Faraday Trans.86(10), 1641–1648 (1990), doi:10.1039/FT9908601641,http://dx.doi.org/10.1039/FT9908601641

    work page doi:10.1039/ft9908601641 1990

  48. [48]

    Weinberg, Quasiparticles and the Born Series.Phys

    S. Weinberg, Quasiparticles and the Born Series.Phys. Rev.131, 440–460 (1963), doi:10. 1103/PhysRev.131.440,https://link.aps.org/doi/10.1103/PhysRev.131.440

    doi:10.1103/physrev.131.440 1963

  49. [49]

    Hoppe, C

    J. Hoppe, C. Drischler, R. J. Furnstahl, K. Hebeler, A. Schwenk, Weinberg eigenvalues for chiral nucleon-nucleon interactions.Phys. Rev. C96, 054002 (2017), doi:10.1103/PhysRevC. 96.054002,https://link.aps.org/doi/10.1103/PhysRevC.96.054002

    doi:10.1103/physrevc 2017

  50. [50]

    Balakrishnan, A

    N. Balakrishnan, A. Dalgarno, R. C. Forrey, Vibrational relaxation of CO by collisions with 4Heatultracoldtemperatures.J.Chem.Phys.113(2),621–627(2000),doi:10.1063/1.481838, https://doi.org/10.1063/1.481838

    work page doi:10.1063/1.481838 2000

  51. [51]

    Cecchi-Pestellini, E

    C. Cecchi-Pestellini, E. Bodo, N. Balakrishnan, A. Dalgarno, Rotational and Vibrational Ex- citation of CO Molecules by Collisions with4He Atoms.Astrophys. J.571(2), 1015 (2002), doi:10.1086/340020,https://doi.org/10.1086/340020

    work page doi:10.1086/340020 2002

  52. [52]

    P. M. Florian, M. Hoster, R. C. Forrey, Rotational relaxation in ultracoldCO+Hecollisions. Phys. Rev. A70, 032709 (2004), doi:10.1103/PhysRevA.70.032709,https://link.aps. org/doi/10.1103/PhysRevA.70.032709

    work page doi:10.1103/physreva.70.032709 2004

  53. [53]

    Chem.7(4), 349–353 (2015), doi: 10.1038/nchem.2204,https://doi.org/10.1038/nchem.2204

    A.Bergeat,J.Onvlee,C.Naulin,A.vanderAvoird,M.Costes,Quantumdynamicalresonances in low-energy CO(𝑗=0) + He inelastic collisions.Nat. Chem.7(4), 349–353 (2015), doi: 10.1038/nchem.2204,https://doi.org/10.1038/nchem.2204

    work page doi:10.1038/nchem.2204 2015

  54. [54]

    Cybulski,et al., Ab initio study of the CO–N2 complex: a new highly accurate inter- molecular potential energy surface and rovibrational spectrum.Phys

    H. Cybulski,et al., Ab initio study of the CO–N2 complex: a new highly accurate inter- molecular potential energy surface and rovibrational spectrum.Phys. Chem. Chem. Phys. 29 20(18), 12624–12636 (2018), doi:10.1039/c8cp01373j,http://dx.doi.org/10.1039/ c8cp01373j

    work page doi:10.1039/c8cp01373j 2018

  55. [55]

    R. D. Stephens, S. H. Cadle, Remote Sensing Measurements of Carbon Monoxide Emissions from On-Road Vehicles.J. Air Waste Manag. Assoc.41(1), 39–46 (1991), doi:10.1080/ 10473289.1991.10466823

    arXiv 1991

  56. [56]

    W. W. McMillan,et al., An observational and modeling strategy to investigate the impact of remote sources on local air quality: A Houston, Texas, case study from the Second Texas Air Quality Study (TexAQS II).J. Geophys. Res.115(D1) (2010), doi:10.1029/2009jd011973

    work page doi:10.1029/2009jd011973 2010

  57. [57]

    Shan,et al., Regional CO emission estimated from ground-based remote sensing at Hefei site, China.Atmos

    C. Shan,et al., Regional CO emission estimated from ground-based remote sensing at Hefei site, China.Atmos. Res.222, 25–35 (2019), doi:10.1016/j.atmosres.2019.02.005

    work page doi:10.1016/j.atmosres.2019.02.005 2019

  58. [58]

    J. C. W. W. W. Carole, N. Beale,Global climate change linkages: acid rain, air quality, and stratospheric ozone(Springer Science & Business Media) (1989)

    1989

  59. [59]

    J. S. Daniel, S. Solomon, On the climate forcing of carbon monoxide.J. Geophys. Res.: Atmospheres103(D11), 13249–13260 (1998), doi:10.1029/98jd00822

    work page doi:10.1029/98jd00822 1998

  60. [60]

    A. Nishiyama,et al., Line-shape study of CO perturbed by N 2 with mid- infrared frequency comb-based Fourier-transform spectroscopy.Measurement227, 114273 (2024), doi:10.1016/j.measurement.2024.114273,http://dx.doi.org/10.1016/ j.measurement.2024.114273

    arXiv 2024

  61. [61]

    D. J. Jacob,Introduction to atmospheric chemistry(Princeton University Press) (1999)

    1999

  62. [62]

    Paredes-Roibás, R

    D. Paredes-Roibás, R. Z. Martínez, H. Jóźwiak, F. Thibault, Collisional line broadening and mixing in the Raman spectrum of CO perturbed by N2: Experimental measurements and theoreticalcalculations.J.Quant.Spectrosc.Radiat.Transf.275,107868(2021),doi:10.1016/ j.jqsrt.2021.107868,http://dx.doi.org/10.1016/j.jqsrt.2021.107868

    arXiv 2021

  63. [63]

    A.Zadrożny,H.Jóźwiak,E.Quintas-Sánchez,R.Dawes,P.Wcisło,Abinitioquantumscatter- ing calculations for the CO–O2 system and a new CO–O2 potential energy surface: O2 and air 30 broadening of the R(0) line in CO.J. Chem. Phys.157(17) (2022), doi:10.1063/5.0115654, http://dx.doi.org/10.1063/5.0115654

    work page doi:10.1063/5.0115654 2022

  64. [64]

    R. B. Lehoucq, D. C. Sorensen, C. Yang,ARPACK users’ guide: solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods(SIAM) (1998)

    1998

  65. [65]

    Lique, A

    F. Lique, A. Faure, eds.,Gas-Phase Chemistry in Space, 2514-3433 (IOP Publishing) (2019), doi:10.1088/2514-3433/aae1b5,https://doi.org/10.1088/2514-3433/aae1b5

    work page doi:10.1088/2514-3433/aae1b5 2019

  66. [66]

    𝑗𝑧 -conserving coupled states approximation.J

    P.McGuire,D.J.Kouri,Quantummechanicalclosecouplingapproachtomolecularcollisions. 𝑗𝑧 -conserving coupled states approximation.J. Chem. Phys.60(6), 2488–2499 (1974), doi: 10.1063/1.1681388,https://doi.org/10.1063/1.1681388

    work page doi:10.1063/1.1681388 1974

  67. [67]

    C. Joy, B. Mandal, D. Bostan, M.-L. Dubernet, D. Babikov, Mixed quantum/classical theory (MQCT) approach to the dynamics of molecule–molecule collisions in complex systems. FaradayDiscuss.251,225–248(2024),doi:10.1039/d3fd00166k,http://dx.doi.org/10. 1039/d3fd00166k

    work page doi:10.1039/d3fd00166k 2024

  68. [68]

    D. G. Truhlar, Long-standing themes in computational chemical dynamics.Comput. Phys. Commun.84(1–3), 78–90 (1994), doi:10.1016/0010-4655(94)90204-6,http://dx.doi. org/10.1016/0010-4655(94)90204-6

    work page doi:10.1016/0010-4655(94)90204-6 1994

  69. [69]

    Leforestier, K

    C. Leforestier, K. Szalewicz, A. van der Avoird, Spectra of water dimer from a newab initio potential with flexible monomers.J. Chem. Phys.137(1), 014305 (2012), doi:10.1063/1. 4722338,https://doi.org/10.1063/1.4722338

    doi:10.1063/1 2012

  70. [70]

    X.-G. Wang, S. Yang, J. Carrington, Tucker, D. H. Zhang, A numerically exact calculation of vibration–rotation–tunnelinglevelsofwaterdimeronanewaccuratepotentialenergysurface: Achieving sub-cm−1 accuracy from the terahertz to the infrared.J. Chem. Phys.163(14), 144308 (2025), doi:10.1063/5.0294206,https://doi.org/10.1063/5.0294206

    doi:10.1063/5.0294206 2025

  71. [71]

    K. A. Peterson, G. C. McBane, A hierarchical family of three-dimensional potential energy surfacesforHe-CO.J.Chem.Phys.123(8)(2005),doi:10.1063/1.1947194,http://dx.doi. org/10.1063/1.1947194. 31

    work page doi:10.1063/1.1947194 2005

  72. [72]

    K. P. Huber, G. Herzberg,Constants of diatomic molecules(Springer US, Boston, MA), pp. 8–689 (1979), doi:10.1007/978-1-4757-0961-2_2,https://doi.org/10.1007/ 978-1-4757-0961-2_2. Acknowledgments We thank George McBane and Jeremy Hutson for stimulating discussions in the early stages of this work and Nathan Prins for valuable comments. H. J. J. was support...

    doi:10.1007/978-1-4757-0961-2_2 1979