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.

Molecular dynamics-driven global potential energy surfaces: Application to the AlF dimer.

In this work, we present a full-dimensional potential energy surface for AlF-AlF. We apply a general machine learning approach for full-dimensional potential energy surfaces, employing an active learning scheme trained on ab initio points, whose size grows based on the accuracy required. The training points are selected based on molecular dynamics simulations, choosing the most suitable configurations for different collision energy and mapping the most relevant part of the potential energy landscape of the system. The present approach does not require long-range information and is entirely general. As a result, it is possible to provide the full-dimensional AlF-AlF potential energy surface, requiring ≲0.01% of the configurations to be calculated ab initio. Furthermore, we analyze the general properties of the AlF-AlF system, finding critical differences with other reported results on CaF or bi-alkali dimers.

Barrier Height Prediction by Machine Learning Correction of Semiempirical Calculations.

Different machine learning (ML) models are proposed in the present work to predict density functional theory-quality barrier heights (BHs) from semiempirical quantum mechanical (SQM) calculations. The ML models include a multitask deep neural network, gradient-boosted trees by means of the XGBoost interface, and Gaussian process regression. The obtained mean absolute errors are similar to those of previous models considering the same number of data points. The ML corrections proposed in this paper could be useful for rapid screening of the large reaction networks that appear in combustion chemistry or in astrochemistry. Finally, our results show that 70% of the features with the highest impact on model output are bespoke predictors. This custom-made set of predictors could be employed by future Δ-ML models to improve the quantitative prediction of other reaction properties.

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.

The performance of CCSD(T) for the calculation of dipole moments in diatomics.

This work analyzes the accuracy of the coupled cluster with single, double, and perturbative triple excitation [CCSD(T)] method for predicting dipole moments. In particular, we benchmark CCSD(T) predictions for the equilibrium bond length, vibrational frequency, and dipole moment versus accurate experimental data. As a result, we find that CCSD(T) leads to accurate dipole moments. However, in some cases, it disagrees with the experimental values, and the disagreement can not be satisfactorily explained via relativistic or multi-reference effects. Therefore, our results indicate that benchmark studies for energy and geometry properties do not accurately describe other electron density magnitudes.

AlF-AlF Reaction Dynamics between 200 K and 1000 K: Reaction Mechanisms and Intermediate Complex Characterization.

AlF is a relevant molecule in astrochemistry as a tracer of F-bearing molecules. Additionally, AlF presents diagonal Franck-Condon factors and can be created very efficiently in the lab, which makes it a prototypical molecular for laser cooling. However, very little is known about the reaction dynamics of AlF. In this work, we report on the reaction dynamics of AlF-AlF between 200 and 1000 K using ab initio molecular dynamics and a highly efficient active learning approach for the potential energy surface, including all degrees of freedom. As a result, we identify the main reaction mechanisms and the lifetime of the intermediate complex AlF-AlF relevant to astrochemistry environments and regions in buffer gas cells.

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.

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.

Complex Reaction Network Thermodynamic and Kinetic Autoconstruction Based on Ab Initio Statistical Mechanics: A Case Study of O2 Activation on Ag4 Clusters.

An approach based on ab initio statistical mechanics is demonstrated for autoconstructing complex reaction networks. Ab initio molecular dynamics combined with Markov state models are employed to study relevant transitions and corresponding thermodynamic and kinetic properties of a reaction. To explore the capability and flexibility of this approach, we present a study of oxygen activation on Ag4 as a model reaction. Specifically, with the same sampled trajectories, it is possible to study the structural effects and the reaction rate of the cited reaction. The results show that this approach is suitable for automatized construction of reaction networks, especially for non-well-studied reactions, which can benefit from this ab initio molecular dynamics based approach to construct comprehensive reaction networks with Markov state models without prior knowledge about the potential energy landscape.

Predicting second virial coefficients of organic and inorganic compounds using Gaussian process regression.

We show that by using intuitive and accessible molecular features it is possible to predict the temperature-dependent second virial coefficient of organic and inorganic compounds with Gaussian process regression. In particular, we built a low dimensional representation of features based on intrinsic molecular properties, topology and physical properties relevant for the characterization of molecule-molecule interactions. The featurization was used to predict second virial coefficients in the interpolative regime with a relative error ⪅1% and to extrapolate the prediction to temperatures outside of the training range for each compound in the dataset with a relative error of 2.1%. Additionally, the model's predictive abilities were extended to organic molecules unseen in the training process, yielding a prediction with a relative error of 2.7%. Test molecules must be well-represented in the training set by instances of their families, which are high in variety. The method shows a generally better performance when compared to several semi-empirical procedures employed in the prediction of the quantity. Therefore, apart from being robust, the present Gaussian process regression model is extensible to a variety of organic and inorganic compounds.

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.

On the relationship between spectroscopic constants of diatomic molecules: a machine learning approach.

Through a machine learning approach, we show that the equilibrium distance, harmonic vibrational frequency and binding energy of diatomic molecules are related, independently of the nature of the bond of a molecule; they depend solely on the group and period of the constituent atoms. As a result, we show that by employing the group and period of the atoms that form a molecule, the spectroscopic constants are predicted with an accuracy of <5%, whereas for the A-excited electronic state it is needed to include other atomic properties leading to an accuracy of <11%.

A data-driven approach to determine dipole moments of diatomic molecules.

We present a data-driven approach for the prediction of the electric dipole moment of diatomic molecules, which is one of the most relevant molecular properties. In particular, we apply Gaussian process regression to a novel dataset to show that dipole moments of diatomic molecules can be learned, and hence predicted, with a relative error ⪅5%. The dataset contains the dipole moment of 162 diatomic molecules, the most exhaustive and unbiased dataset of dipole moments up to date. Our findings show that the dipole moment of diatomic molecules depends on atomic properties of the constituents atoms: electron affinity and ionization potential, as well as on (a feature related to) the first derivative of the electronic kinetic energy at the equilibrium distance.

The diatomic molecular spectroscopy database.

MOTIVATION: The spectroscopy of diatomic molecules is an important research area in chemical physics due to its relevance in astrochemistry, combustion chemistry, and ultracold physics. However, there is currently no database where the user can easily retrieve, in a useful format, the spectroscopic constants of a given molecule. A similar situation appears concerning the vibrational Franck-Condon factors for diatomic molecules, a crucial parameter to infer laser cooling prospects for molecules. To address this problem, and inspired by the idea that data should be open and freely accessible, we have developed a user-friendly website (https://rios.mp.fhi.mpg.de) where the user can retrieve spectroscopic constants and Franck-Condon factors in useful formats. IMPLEMENTATION: In this database, the spectroscopic constants of the ground states and first excited states of the diatomic molecules are accessible from the website and can be retrieved in readable formats. The website is implemented within the LAMP web service stacks. In particular, using Linux as the operative system, Apache as the HTTP Server, MySQL as the database management system, and PHP as the programming language for the web. Furthermore, the user can register and upload new data. This project is licensed under the Free-Libre/Open Source Software (FLOSS) license Apache License 2.0 which allows free and open access to the codes as well as efficient collaboration in the maintenance of the software. CONCLUSIONS AND IMPACT: The present data-driven website presents essential information in a user-friendly manner and may help the chemical physics community to identify molecules that should be explored through spectroscopic techniques.

Observation of Quantum Interference and Coherent Control in a Photochemical Reaction.

Coherent control of reactants remains a long-standing challenge in quantum chemistry. In particular, we have studied laser-induced molecular formation (photoassociation) in a Raman-dressed spin-orbit-coupled ^{87}Rb Bose-Einstein condensate, whose spin quantum state is a superposition of multiple bare spin components. In contrast to the notably different photoassociation-induced fractional atom losses observed for the bare spin components of a statistical mixture, a superposition state with a comparable spin composition displays the same fractional loss on every spin component. We interpret this as the superposition state itself undergoing photoassociation. For superposition states induced by a large Raman coupling and zero Raman detuning, we observe a nearly complete suppression of the photoassociation rate. This suppression is consistent with a model based upon quantum destructive interference between two photoassociation pathways for colliding atoms with different spin combinations. This model also explains the measured dependence of the photoassociation rate on the Raman detuning at a moderate Raman coupling. Our work thus suggests that preparing atoms in quantum superpositions may represent a powerful new technique to coherently control photochemical reactions.

Classical Fractals and Quantum Chaos in Ultracold Dipolar Collisions.

We examine a dipolar-gas model to address fundamental issues regarding the correspondence between classical chaos and quantum observations in ultracold dipolar collisions. The theoretical model consists of a short-range Lennard-Jones potential well with an anisotropic, long-range dipole-dipole interaction between two atoms. Both the classical and quantum dynamics are explored for the same Hamiltonian of the system. The classical chaotic scattering is revealed by the fractals developed in the scattering function (defined as the final atom separation as a function of initial conditions), while the quantum chaotic features lead to the repulsion of the eigenphases from the corresponding quantum S matrix. The nearest-eigenphase-spacing statistics have an intermediate behavior between the Poisson and the Wigner-Dyson distributions. The character of the distribution can be controlled by changing an effective Planck constant or the dipole moment. The degree of quantum chaos shows a good correspondence with the overall average of the classical scattering function. The results presented here also provide helpful insights for understanding the role of the inherent dipole-dipole interaction in the currently ongoing experiments on ultracold collisions of highly magnetic atoms.

Observation of pendular butterfly Rydberg molecules.

Engineering molecules with a tunable bond length and defined quantum states lies at the heart of quantum chemistry. The unconventional binding mechanism of Rydberg molecules makes them a promising candidate to implement such tunable molecules. A very peculiar type of Rydberg molecules are the so-called butterfly molecules, which are bound by a shape resonance in the electron-perturber scattering. Here we report the observation of these exotic molecules and employ their exceptional properties to engineer their bond length, vibrational state, angular momentum and orientation in a small electric field. Combining the variable bond length with their giant dipole moment of several hundred Debye, we observe counter-intuitive molecules which locate the average electron position beyond the internuclear distance.

Effective Atom-Molecule Conversions Using Radio Frequency Fields.

The present study is inspired by the Wieman group experiment [Phys. Rev. Lett. 2005, 95, 190404], in which they use a slow modulated magnetic field to effectively transfer rubidium atoms into cold molecules near a Feshbach resonance. We develop a time-dependent collision theory based on two channel model potentials to study the atom-molecule population transfer induced by a single-color radio frequency field in an ultracold 87 Rb gas. Wave-packet dynamical simulations allow an investigation of both bound-bound transitions and free-bound transitions. The effects of temperature, detuning and the RF amplitude on the population transfer are discussed in detail. Some of our simulations suggest that oscillatory atom-molecule conversion could originate from the long coherence time of the wave packet. This coherence time is unusually long in ultracold gases because the collision energy is typically quite well-defined.

Energy Scaling of Cold Atom-Atom-Ion Three-Body Recombination.

We study three-body recombination of Ba^{+}+Rb+Rb in the mK regime where a single ^{138}Ba^{+} ion in a Paul trap is immersed into a cloud of ultracold ^{87}Rb atoms. We measure the energy dependence of the three-body rate coefficient k_{3} and compare the results to the theoretical prediction, k_{3}∝E_{col}^{-3/4}, where E_{col} is the collision energy. We find agreement if we assume that the nonthermal ion energy distribution is determined by at least two different micromotion induced energy scales. Furthermore, using classical trajectory calculations we predict how the median binding energy of the formed molecules scales with the collision energy. Our studies give new insights into the kinetics of an ion immersed in an ultracold atom cloud and yield important prospects for atom-ion experiments targeting the s-wave regime.