Py3BR: A software for computing atomic three-body recombination rates.

The three-body recombination reaction, or ternary association, is a termolecular reaction leading to a molecule after a three-body encounter that plays a vital role in many relevant scenarios in chemical physics. Here, we introduce the Python 3-Body Recombination program, which is dedicated to the computation of atomic three-body recombination rate coefficients. The software is based on a classical trajectory approach in hyperspherical coordinates after mapping the three-body problem as a single particle in a higher-dimensional space. This theoretical approach is fully general and applicable to any ion-atom-atom or atom-atom-atom three-body process. The predictive power of the methodology has been tested in several different experimental scenarios, reaching a good description of every system. The code structure is presented alongside examples and tests to ensure the software's capacity. In addition, the performance of the software after parallelization is shown.

On the role of non-additive interactions in three-body recombination.

Non-additive forces are a cornerstone of molecular spectroscopy and reaction dynamics. However, the relevance of non-additive forces in three-body recombination remains largely unexplored. In this work, we present a global study on the impact of non-additive interactions in three-body recombination: atom-atom-atom and ion-atom-atom. Our study explores these reactions in a wide range of energies, from the cold to the hyperthermal regime, finding no effect of non-additive interactions. Therefore, pair-wise interactions are enough to describe the three-body recombination dynamics adequately.

Ion-atom-atom three-body recombination: From the cold to the thermal regime.

We present a study on ion-atom-atom reaction A + A + B+ in a wide range of systems and collision energies ranging from 100 µK to 105 K, analyzing two possible products: molecules and molecular ions. The dynamics is performed via a direct three-body formalism based on a classical trajectory method in hyperspherical coordinates developed in Pérez-Ríos et al. [J. Chem. Phys. 140, 044307 (2014)]. Our chief finding is that the dissociation energy of the molecular ion product acts as a threshold energy, separating the low- and high-energy regimes. In the low-energy regime, the long-range tail of the three-body potential dictates the fate of the reaction and the main reaction product. On the contrary, in the high-energy regime, the short-range of atom-atom and atom-ion interaction potential dominate the dynamics, enhancing molecular formation.

Ozone Formation in Ternary Collisions: Theory and Experiment Reconciled.

The present Letter shows that the formation of ozone in ternary collisions O+O_{2}+M-the primary mechanism of ozone formation in the stratosphere-at temperatures below 200 K (for M=Ar) proceeds through a formation of a temporary complex MO_{2}, while at temperatures above ∼700 K, the reaction proceeds mainly through a formation of long-lived vibrational resonances of O_{3}^{*}. At intermediate temperatures 200-700 K, the process cannot be viewed as a two-step mechanism, often used to simplify and approximate collisions of three atoms or molecules. The developed theoretical approach is applied to the reaction O+O_{2}+Ar because of extensive experimental data available. The rate coefficients for the formation of O_{3} in ternary collisions O+O_{2}+Ar without using two-step approximations were computed for the first time as a function of collision energy. Thermally averaged coefficients were derived for temperatures 5-900 K. It is found that the majority of O_{3} molecules formed initially are weakly bound. Accounting for the process of vibrational quenching of the nascent population, a good agreement with available experimental data for temperatures 100-900 K is obtained.

On the formation of van der Waals complexes through three-body recombination.

In this work, we show that van der Waals molecules X-RG (where RG is the rare gas atom) may be created through direct three-body recombination collisions, i.e., X + RG + RG → X-RG + RG. In particular, the three-body recombination rate at temperatures relevant for buffer gas cell experiments is calculated via a classical trajectory method in hyperspherical coordinates [Pérez-Ríos et al., J. Chem. Phys. 140, 044307 (2014)]. As a result, it is found that the formation of van der Waals molecules in buffer gas cells (1 K ≲ T ≲ 10 K) is dominated by the long-range tail (distances larger than the LeRoy radius) of the X-RG interaction. For higher temperatures, the short-range region of the potential becomes more significant. Moreover, we notice that the rate of formation of van der Walls molecules is of the same order of the magnitude independent of the chemical properties of X. As a consequence, almost any X-RG molecule may be created and observed in a buffer gas cell under proper conditions.

Classical threshold law for the formation of van der Waals molecules.

We study the role of pairwise long-range interactions in the formation of van der Waals molecules through direct three-body recombination processes A + B + B → AB + B, based on a classical trajectory method in hyperspherical coordinates developed in our earlier works [J. Pérez-Ríos et al., J. Chem. Phys. 140, 044307 (2014); M. Mirahmadi and J. Pérez-Ríos, J. Chem. Phys. 154, 034305 (2021)]. In particular, we find the effective long-range potential in hyperspherical coordinates with an exact expression in terms of dispersion coefficients of pairwise potentials. Exploiting this relation, we derive a classical threshold law for the total cross section and the three-body recombination rate yielding an analytical expression for the three-body recombination rate as a function of the pairwise long-range coefficients of the involved partners.

Dynamics of polar polarizable rotors acted upon by unipolar electromagnetic pulses: From the sudden to the adiabatic regime.

We study, analytically as well as numerically, the dynamics that arises from the interaction of a polar polarizable rigid rotor with single unipolar electromagnetic pulses of varying length, Δτ, with respect to the rotational period of the rotor, τ r . In the sudden, non-adiabatic limit, Δτ ≪ τ r , we derive analytic expressions for the rotor's wavefunctions, kinetic energies, and field-free evolution of orientation and alignment. We verify the analytic results by solving the corresponding time-dependent Schrödinger equation numerically and extend the temporal range of the interactions considered all the way to the adiabatic limit, Δτ > τ r , where general analytic solutions beyond the field-free case are no longer available. The effects of the orienting and aligning interactions as well as of their combination on the post-pulse populations of the rotational states are visualized as functions of the orienting and aligning kick strengths in terms of population quilts. Quantum carpets that encapsulate the evolution of the rotational wavepackets provide the space-time portraits of the resulting dynamics. The population quilts and quantum carpets reveal that purely orienting, purely aligning, or even-break combined interactions each exhibit sui generis dynamics. In the intermediate temporal regime, we find that the wavepackets as functions of the orienting and aligning kick strengths show resonances that correspond to diminished kinetic energies at particular values of the pulse duration.