arXiv:2610.20882

A tutorial on quantum scattering theory and coupled-channel calculations of bimolecular collisions

Karl P. Horn, Christiane P. Koch

quant-phphysics.atom-phphysics.chem-phSubmitted Oct 14, 2026• Announced Oct 16, 2026

Abstract

Collisions between atoms and molecules govern how energy is exchanged and redistributed in the gas phase, underpinning phenomena from spectral line broadening and interstellar molecular excitation to the cooling and trapping of molecules. A quantitative account of these collisions requires connecting the theoretical description of the interaction between the collision partners to measurable quantities such as cross sections. However, scattering theory as presented in quantum mechanics textbooks does not readily translate into numerical methods. Here, we formulate the scattering problem in Dirac notation, as customary in present-day quantum mechanics courses, and show how to build a practical numerical framework from it. We derive the "coupled channel" representation of the scattering problem with correct asymptotic behaviour and introduce Numerov's method to obtain time-independent wavefunctions, from which we calculate cross sections, the observables measured in scattering experiments. Starting from single channel scattering, we build up to the realistic collision of a ground-state atom with a diatomic molecule.

Related papers

Ultracold collisions of molecules

Goulven Quéméner

This paper deals with the theory of collisions between two ultracold particles with a special focus on molecules. It describes the general features of the scattering theory of two particles with internal structure, using a time-independent quantum formalism. It starts from the Schrödinger equation and introduces the experimental observables such as the differential or integral cross sections, and rate coefficients. Using a partial-wave expansion of the scattering wavefunction, the radial motion of the collision is described through a linear system of coupled equations, which is solved numerically. Using a matching procedure of the scattering wavefunction with its asymptotic form, the observables such as cross sections and rate coefficients are obtained from the extraction of the reactance, scattering and transition matrices. The example of the collision of two dipolar molecules in the presence of an electric field is presented, showing how dipolar interactions and collisions can be controlled.

arXiv:1703.09174Mar 27, 20174 citations

CHAPTER 12. Ultracold Collisions of Molecules

Goulven Quéméner

This chapter deals with the theory of collisions between two ultracold particles with a special focus on molecules. It describes the general features of the scattering theory of two particles with internal structure, using a time-independent quantum formalism. It starts from the Schrodinger equation and introduces experimental observables such as the differential or integral cross-sections, and rate coefficients. Using a partial-wave expansion of the scattering wave function, the radial motion of the collision is described through a linear system of coupled equations, which is solved numerically. Using a matching procedure of the scattering wavefunction with its asymptotic form, the observables such as cross sections and rate coefficients are obtained from the extraction of the reactance, scattering and transition matrices. The example of the collision of two dipolar molecules in the presence of an electric field is presented, showing how dipolar interactions and collisions can be controlled.

RSC theoretical and computational chemistry seriesJan 1, 20173 citations

Modern scattering theory

Manuel Valiente, Nikolaj Thomas Zinner

Almost any of the great textbooks available on scattering theory continues to explain how a collision experiment measures scattering cross sections for different events, and relates the relevant S-matrix elements, or the scattering amplitudes (onshell T-matrix), to those cross sections. We highly recommend Taylor’s book in this respect. In this volume and the volume that will follow, however, our motivation does not lie in (although discussions can be extended to) the theoretical description of collision experiments. We are mostly interested in applications of scattering theory in few-body physics and effective interactions in few- and many-body physics. Therefore, we consider a method that is having huge success in the theoretical and numerical description of collisions in nuclear physics and quantum chromodynamics, as well as ultracold atomic and condensed matter physics recently. The method allows for an estimation of low-energy scattering parameters and binding energies from exact diagonalisation, in the low-energy sector, of the system’s Hamiltonian, and is originally due to Lüscher [2]. Lacking another term for it, we have decided to simply label this and related methods ‘modern scattering theory’.

Institute of Physics eBooksSep 1, 20230 citations

A Hybrid Quantum-Classical Algorithm for Multichannel Quantum Scattering of Atoms and Molecules

Xiaodong Xing, Alejandro Gomez Cadavid, Artur F. Izmaylov +1 more

We propose a hybrid quantum-classical algorithm for solving the time-independent Schrödinger equation for atomic and molecular collisions. The algorithm is based on the S -matrix version of the Kohn variational principle, which computes the fundamental scattering S -matrix by inverting the Hamiltonian matrix expressed in the basis of square-integrable functions. The computational bottleneck of the classical algorithm─symmetric matrix inversion─is addressed here using the variational quantum linear solver (VQLS), a recently developed noisy intermediate-scale quantum (NISQ) algorithm for solving systems of linear equations. We apply our algorithm to single- and multichannel quantum scattering problems, obtaining accurate vibrational relaxation probabilities in collinear atom–molecule collisions. We also show how the algorithm could be scaled up to simulate collisions of large polyatomic molecules. Our results demonstrate that it is possible to calculate scattering cross sections and rates for complex molecular collisions on NISQ quantum processors, opening up the possibility of scalable digital quantum computation of gas-phase bimolecular collisions and reactions of relevance to astrochemistry and ultracold chemistry.

arXiv:2304.06089The Journal of Physical Chemistry LettersJun 30, 202312 citations

Time-Dependent Wave Packet Dynamics Calculations of Cross Sections for Ultracold Scattering of Molecules

Jiayu Huang, Shu Liu, Dong Hui Zhang +1 more

Because the de Broglie wavelength of ultracold molecules is very large, the cross sections for collisions of molecules at ultracold temperatures are always computed by the time-independent quantum scattering approach. Here, we report the first accurate time-dependent wave packet dynamics calculation for reactive scattering of ultracold molecules. Wave packet dynamics calculations can be applied to molecular systems with more dimensions and provide real-time information on the process of bond rearrangement and/or energy exchange in molecular collisions. Our work thus makes possible the extension of rigorous quantum calculations of ultracold reaction properties to polyatomic molecules and adds a new powerful tool for the study of ultracold chemistry.

arXiv:1710.09057Physical Review LettersApr 4, 201844 citations

Cold collisions of metastable neon atoms

Christian Cop

In this thesis we investigate theoretically the cold collision properties of metastable neon (Ne*). The collisions of metastable rare gases (Rg*) exhibit a plethora of interesting phenomena: Due to the high internal excitation energies of the metastable atoms ionization processes are very likely. The produced ions can be detected with single-ion efficiency in the experiment. In experiments with atomic Rg* gases at cold temperatures T~1 mK, this allows for in-situ monitoring of the real-time dynamics of the atomic gas and for example for the study of two-particle correlations. Furthermore, in cold Rg* gases one can examine the prospects of creating a Bose-Einstein condensate or a degenerate Fermi gas of Rg*. Ne* is a promising candidate for Bose-Einstein condensation and in the group of G. Birkl the properties of cold gases of Ne* are investigated experimentally. Two-body loss rate coefficients due to the ionization processes and elastic cross sections have been measured. It is useful to complement the experimental measurements with a theoretical study in order to obtain a better understanding of the cold Ne* gases and the collision physics of Ne*. At the low temperatures of the atomic gases, these gases are dilute with particle number densities of n~10^9 cm^{-3}. For these low densities, the behavior of the gas is determined mainly by two-body collision physics which is solved by quantum scattering theory in terms of the S matrix. The collision physics is given by short-range interactions and long-range interactions of the atoms. The short-range and long-range interaction potentials of Ne* have been calculated in the molecular basis of Ne2. In this thesis we first calculate short-range and long-range interactions of atoms and molecular basis states of diatomic molecules with a simpler electronic structure than Ne* and Ne2 in order to introduce the notation and to discuss the characteristics of the calculated molecular interaction potentials and in order to employ them in the scattering calculations for Ne*. Instead of solving the full scattering equations by taking into account all the molecular interaction potentials, we demonstrate that the collision physics of Ne* can be described in terms of a coupled two-channel model with a single interaction potential only which describes the elastic scattering of Ne*. In this model, ionization in Ne* collisions is described by the transition of the upper elastic interaction channel to the lower channel, representing the loss or ionization channel. We introduce two versions of the two-channel model in this work. In the first version, the two channels are given by square-well potentials. We can solve the scattering equations analytically and for complex wave numbers k and study the solutions in the complex k plane. With an expansion of the S matrix in terms of its poles in the complex k plane we find a parametrization of the two-body loss rate coefficients. We show that these coefficients describe ionizing collisions of Ne* by comparing them to the experimental measurements and to the numerical results obtained by the two-channel model in the second version. In this second version of the two-channel model, the elastic scattering channel is given by a realistic interaction potential of Ne*, consisting of a calculated short-range and long-range molecular potential, and the loss channel by a model ionization potential. For this model we calculate the two-body loss rate coefficients and the elastic cross sections for the isotope mixtures of Ne* which have been measured experimentally and demonstrate that the free potential parameters can be optimized to the experimental data to achieve very good agreement of the numerical results with the experimental measurements. We discuss the validity of the two-channel model by comparing the results to existing models for cold Rg* collisions and find that the two-channel model of this work is a useful extension to these models.

Jan 1, 20170 citations