On the role of van der Waals interaction in chemical reactions at low temperatures

Exploring HCl-HCl interactions: QZVPP calculations, improved Lennard-Jones potential, and second virial coefficient analysis for thermodynamics and industrial applications

In this paper, we present a comprehensive analysis of HCl-HCl interactions, including QZVPP calculations, energy fitting, conformation validation, and the determination of the second virial coefficient B using improved Lennard-Jones (ILJ) potential parameters. To acquire accurate interaction energies, initial QZVPP calculations are performed on approximately 1851 randomly generated HCl-HCl conformations. Then, these energies are used to fit an improved Lennard-Jones potential energy surface, allowing for a robust description of HCl-HCl interactions. The ILJ potential parameters are then used to validate particular HCl dimer conformations, ensuring their stability and consistency with experimental observations. The correlation between calculated and experimental conformations strengthens the validity of the ILJ potential parameters. In addition, the second viral coefficient B is calculated at various temperatures using the ILJ potential. The obtained B values are compared to experimental data, demonstrating close agreement, and validating the ILJ potential's ability to accurately capture the intermolecular interactions and gas-phase behavior of the HCl-HCl system. The results of this study demonstrate the effective implementation of QZVPP calculations, energy fitting, and ILJ potential parameters in validating HCl-HCl conformations and accurately determining the second virial coefficient B. The high degree of concordance between calculated B values and experimental data demonstrates the validity of the ILJ potential and its suitability for modeling HCl-HCl interactions. This research contributes to a greater comprehension of HCl-HCl interactions and their implications for numerous chemical and atmospheric processes. The validated conformations, energy fitting method, and calculated second virial coefficients provide valuable instruments for future research and pave the way for more accurate modeling and simulations of HCl-HCl systems.

Quantum-Controlled Collisions of H2 Molecules

The amount of information that can be obtained from a scattering experiment depends upon the precision with which the quantum states are defined in the incoming channel. By precisely defining the incoming states and measuring the outgoing states in a scattering experiment, we set up the boundary condition for experimentally solving the Schrödinger equation. In this Perspective we discuss cold inelastic scattering experiments using the most theoretically tractable H₂ and its isotopologues as the target. We prepare the target in a precisely defined rovibrational (v, j, m) quantum state using a special coherent optical technique called the Stark-induced adiabatic Raman passage (SARP). v and j represent the quantum numbers of the vibrational and rotational energy levels, and m refers to the projection of the rotational angular momentum vector j on a suitable quantization axis in the laboratory frame. Selection of the m quantum numbers defines the alignment of the molecular frame, which is necessary to probe the anisotropic interactions. For us to achieve the collision temperature in the range of a few degrees Kelvin, we co-expand the colliding partners in a mixed supersonic beam that is collimated to define a direction for the collision velocity. When the bond axis is aligned with respect to a well-defined collision velocity, SARP achieves stereodynamic control at the quantum scale. Through various examples of rotationally inelastic cold scattering experiments, we show how SARP coherently controls the dynamics of anisotropic interactions by preparing quantum superpositions of the orientational m states within a single rovibrational (v, j) energy state. A partial wave analysis, which has been developed for the cold scattering experiments, shows dominance of a resonant orbital that leaves its mark in the scattering angular distribution. These highly controlled cold collision experiments at the single partial wave limit allow the most direct comparison with the results of theoretical computations, necessary for accurate modeling of the molecular interaction potential.

Evidence for the association of triatomic molecules in ultracold 23Na40K + 40K mixtures

Ultracold assembly of diatomic molecules has enabled great advances in controlled chemistry, ultracold chemical physics and quantum simulation with molecules^1-3. Extending the ultracold association to triatomic molecules will offer many new research opportunities and challenges in these fields. A possible approach is to form triatomic molecules in a mixture of ultracold atoms and diatomic molecules by using a Feshbach resonance between them^4,5. Although ultracold atom-diatomic-molecule Feshbach resonances have been observed recently^6,7, using these resonances to form triatomic molecules remains challenging. Here we report on evidence of the association of triatomic molecules near the Feshbach resonance between ²³Na⁴⁰K molecules in the rovibrational ground state and ⁴⁰K atoms. We apply a radio-frequency pulse to drive the free-bound transition in ultracold mixtures of ²³Na⁴⁰K and ⁴⁰K and monitor the loss of ²³Na⁴⁰K molecules. The association of triatomic molecules manifests itself as an additional loss feature in the radio-frequency spectra, which can be distinguished from the atomic loss feature. The observation that the distance between the association feature and the atomic transition changes with the magnetic field provides strong evidence for the formation of triatomic molecules. The binding energy of the triatomic molecules is estimated from the measurements. Our work contributes to the understanding of the complex ultracold atom-molecule Feshbach resonances and may open up an avenue towards the preparation and control of ultracold triatomic molecules.

Feshbach resonances in D + HD(v = 1, j = 0) reaction at low collision energies

Feshbach resonances in D + HD(v = 1, j = 0) reaction are studied by using the time-independent quantum method. The integral cross section (ICS) results present three Feshbach resonance peaks, which are different from H + HD(v = 1, j = 0) reaction dominated by only one peak. These resonances are attributed to coupling with adiabatic effective potentials of D + HD(v = 1, j = 1) reaction, and the most obvious peak is contributed by J = 1 at 83.16 cm^-1 collision energy. For J = 0 and 2, the resonances are related with the same L partial wave and present a double-peak structure in total ICS. The characteristics of product angular distribution show that the resonance of J = 1 is long-lived, while the lifetimes are relatively shorter for the resonance of J = 0 and 2.

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.

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.

Cold molecules: Progress in quantum engineering of chemistry and quantum matter

Cooling atoms to ultralow temperatures has produced a wealth of opportunities in fundamental physics, precision metrology, and quantum science. The more recent application of sophisticated cooling techniques to molecules, which has been more challenging to implement owing to the complexity of molecular structures, has now opened the door to the longstanding goal of precisely controlling molecular internal and external degrees of freedom and the resulting interaction processes. This line of research can leverage fundamental insights into how molecules interact and evolve to enable the control of reaction chemistry and the design and realization of a range of advanced quantum materials.

Orientation effects in Cl + H2 inelastic collisions: characterization of the mechanisms

Based on quantum mechanical scattering (QM) calculations, we have analyzed the polarization of the product hydrogen molecule in Cl + H(2) (v = 0, j = 0) inelastic collisions. The spatial arrangements adopted by the rotational angular momentum and internuclear axis of the departing molecule have been characterized and used to prove that two distinct mechanisms, corresponding to different dynamical regimes, are responsible for the inelastic collisions. Such mechanisms, named as low-b and high-b, correlate with well defined ranges of impact parameter values, add in an essentially incoherent way, and can be clearly differentiated through the quantum mechanical polarization moment that measures the orientation of the products rotational angular momentum with respect to the scattering plane. Other directional effects turn out to fail when it comes to distinguishing the mechanisms. Quasiclassical trajectories (QCT) calculations have been used as a supplement to the purely quantum mechanical analysis. By combining QM and QCT results, which are in very good agreement, we have succeeded in obtaining a clear and meaningful picture of how the two types of collisions take place.

Dynamical regimes on the Cl + H2 collisions: inelastic rainbow scattering

While Cl + H(2) reactive collisions have been a subject of numerous experimental and theoretical studies, inelastic collisions leading to rotational energy transfer and/or vibrational excitation have been largely ignored. In this work, extensive quantum mechanical calculations covering the 0.5-1.5 eV total energy range and various initial rovibrational states have been carried out and used to perform a joint study of inelastic and reactive Cl + H(2) collisions. Quasiclassical trajectories calculations complement the quantum mechanical results. The analysis of the inelastic transition probabilities has revealed the existence of two distinct dynamical regimes that correlate with low and high impact parameters, b, and are neatly separated by glory scattering. It has been found that while high-b collisions are mainly responsible for |Δj| = 2 transitions which dominate the inelastic scattering, they are very inefficient in promoting higher |Δj| transitions. The effectiveness of this type of collision also drops with rotational excitation of H(2). In contrast, reactive scattering, that competes with |Δj| > 2 inelastic transitions, is exclusively caused by low-b collisions, and it is greatly favored when the reactants get rotationally excited. Previous studies focusing on the reactivity of the Cl + H(2) system established that the van der Waals well located in the entrance channel play a key role in determining the mechanism of the collisions. Our results prove this to be also a case for inelastic processes, where the origin of the double dynamical regime can be traced back to the influence exerted by this well that shapes the topology of the entrance channel of the Cl-H(2) system.

Infrared Spectroscopy of Prereactive Aluminum−, Gallium−, and Indium−HCN Entrance Channel Complexes Solvated in Helium Nanodroplets

Prereactive metal atom-HCN entrance channel complexes [M-HCN (M=Al, Ga, In)] have been stabilized in helium nanodroplets. Rotationally resolved infrared spectra are reported for the CH stretching vibration of the linear nitrogen-bound HCN-Ga and HCN-In complexes that show significant perturbation due to spin-orbit coupling of the 2Pi1/2 ground state with the 2Sigma1/2 state which are degenerate at long range. Six unresolved bands are also observed and assigned to the linear hydrogen-bound isomers of Al-HCN, Ga-HCN, and In-HCN corresponding to the fundamental CH stretching vibration and a combination band involving the CH stretch plus intermolecular stretch for each isomer. A nitrogen-bound HCN-Al complex is not observed, which is attributed to reaction, even at 0.37 K. This conclusion is supported by the observation of a weakly bound complex containing two HCN's and one Al atom which, from the analysis of its rotationally resolved zero-field and Stark spectra is assigned to a weakly bound complex of a HCNAl reaction product and a second HCN molecule. Theoretical calculations are presented to elucidate the reaction mechanisms and energetics of these metal atom reactions with HCN.

Dynamics of the O(3P) + H2 reaction at low temperatures: comparison of quasiclassical trajectory with quantum scattering calculations

Quasiclassical trajectory and quantum-mechanical scattering calculations are reported for the O((3)P) + H(2)(X (1)Sigma(g) (+);upsilon = 1-3,j = 0)-->OH(X (2)Pi) + H((2)S) reaction at energies close to the reaction threshold. The dynamics of the reaction have been investigated for zero total angular momentum using the lowest (3)A" potential-energy surface developed by Rogers et al. [J. Phys. Chem. A 104, 2308 (2000)] and its recent extensions by Brandao et al. [J. Chem. Phys. 121, 8861 (2004)] which provide an improved description of the van der Waals interaction. Good agreement is observed for this system between quasiclassical and quantal results for incident kinetic energies above the tunneling regime. Quantum-mechanical calculations also confirm recent theoretical predictions of a strong collision-energy dependence of the OH(v(') = 0)OH(v' = 1) product branching ratio in the O((3)P) + H(2)(v = 1) reaction, which explains the differences observed in OH vibrational populations between experiments using different O((3)P) sources.

A high-resolution infrared spectroscopic investigation of the halogen atom–HCN entrance channel complexes solvated in superfluid helium droplets

Rotationally resolved infrared spectra are reported for the X-HCN (X = Cl, Br, I) binary complexes solvated in helium nanodroplets. These results are directly compared with those obtained previously for the corresponding X-HF complexes [J. M. Merritt, J. Küpper and R. E. Miller, Phys. Chem. Chem. Phys., 2005, 7, 67]. For bromine and iodine atoms complexed with HCN, two linear structures are observed and assigned to the (2)Sigma(1/2) and (2)Pi(3/2) ground electronic states of the nitrogen and hydrogen bound geometries, respectively. Experiments for HCN + chlorine atoms give rise to only a single band which is attributed to the nitrogen bound isomer. That the hydrogen bound isomer is not stabilized is rationalized in terms of a lowering of the isomerization barrier by spin-orbit coupling. Theoretical calculations with and without spin-orbit coupling have also been performed and are compared with our experimental results. The possibility of stabilizing high-energy structures containing multiple radicals is discussed, motivated by preliminary spectroscopic evidence for the di-radical Br-HCCCN-Br complex. Spectra for the corresponding molecular halogen HCN-X(2) complexes are also presented.

Infrared laser spectroscopy of the CH3–HCN radical complex stabilized in helium nanodroplets

The CH3-HCN and CD3-HCN radical complexes have been formed in helium nanodroplets by sequential pickup of a CH3 (CD3) radical and a HCN molecule and have been studied by high-resolution infrared laser spectroscopy. The complexes have a hydrogen-bonded structure with C3v symmetry, as inferred from the analysis of their rotationally resolved nu = 1 <-- 0 H-CN vibrational bands. The A rotational constants of the complexes are found to change significantly upon vibrational excitation of the C-H stretch of HCN within the complex, DeltaA = A'-A" = -0.04 cm(-1) (for CH3-HCN), whereas the B rotational constants are found to be 2.9 times smaller than that predicted by theory. The reduction in B can be attributed to the effects of helium solvation, whereas the large DeltaA is found to be a sensitive probe of the vibrational averaging dynamics of such weakly bound systems. The complex has a permanent electric dipole moment of 3.1 +/- 0.2 D, as measured by Stark spectroscopy. A vibration-vibration resonance is observed to couple the excited C-H stretching vibration of HCN within the complex to the lower-frequency C-H stretches of the methyl radical. Deuteration of the methyl radical was used to detune these levels from resonance, increasing the lifetime of the complex by a factor of 2. Ab initio calculations for the energies and molecular parameters of the stationary points on the CN+CH4 --> HCN+CH3 potential-energy surface are also presented.

Reactivity enhancement of ultracold O(P3)+H2 collisions by van der Waals interactions

The role of van der Waals forces in O((3)P)+H(2)(upsilon=1,j=0) collisions is investigated theoretically at low and ultralow temperatures. Quantum scattering calculations have been performed for zero total angular momentum using the lowest London-Eyring-Polanyi-Sato double-polynomial (3)A(") potential-energy surface reported by [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)] and its recent BMS1 and BMS2 extensions developed by [Brandao et al., J. Chem. Phys. 121, 8861 (2004)] which provide a more accurate treatment of the van der Waals interaction. Our calculations show that van der Waals forces strongly influence chemical reactivity at ultracold translational energies. The presence of a zero-energy resonance for the BMS1 surface is found to enhance reactivity in the ultracold regime and shift the Wigner threshold to lower temperatures.

Quantum dynamics of the Li+HF→H+LiF reaction at ultralow temperatures

Quantum-mechanical calculations are reported for the Li+HF(v=0,1,j=0)-->H+LiF(v',j') bimolecular scattering process at low and ultralow temperatures. Calculations have been performed for zero total angular momentum using a recent high-accuracy potential-energy surface for the X2A' electronic ground state. For Li+HF(v=0,j=0), the reaction is dominated by resonances due to the decay of metastable states of the Li cdots,...F-H van der Waals complex. Assignment of these resonances has been carried out by calculating the eigenenergies of the quasibound states. We also find that while chemical reactivity is greatly enhanced by vibrational excitation, the resonances get mostly washed out in the reaction of vibrationally excited HF with Li atoms. In addition, we find that at low energies, the reaction is significantly suppressed due to the less-efficient tunneling of the relatively heavy fluorine atom.

Quantum mechanical investigation of rovibrational relaxation of H2 and D2 by collisions with Ar atoms

We report quantum mechanical calculations of cross sections and rate coefficients for rovibrational relaxation of H2 and D2 by collisions with Ar atoms over a wide range of temperatures including the ultracold limit. Limiting values of the rate coefficients for vibrational and rotational quenching at zero temperature were computed and sensitivity of the results to the choice of the interaction potential is investigated. We also demonstrate dramatic change in the behavior of the rate coefficients at low temperatures when the van der Waals potential supports a quasibound level very close to the dissociation threshold.

Theoretical Study of H2...I- van der Waals Anion Complex

The ab initio potential energy surface (PES) for the weak interaction of hydrogen molecule with iodine anion is presented. The surface was obtained by the supermolecular method at the MP4(SDTQ) level of theory. Our calculations indicate the van der Waals (vdW) system for the linear configuration at rH-H = 0.752 Å and R = 3.76 Å with a well depth of De = 2096 μEh. The presented PES reveals also a transition state for the perpendicular arrangement at rH-H = 0.7416 Å and R = 4.63 Å with an interaction energy of -113 μEh. The physical origin of stability of the vdW H2...I- structure with respect to the H2...X- (X = F, Cl, Br) one was analysed by the symmetry adapted perturbation theory (SAPT) based on the single determinant HF wave function. The separation of the interaction energy shows that the dispersion forces play a much more important role for the systems with Cl, Br and I than for H2...F- and their importance slightly increases in the order Cl < Br < I. The global importance of the electrostatic and the induction energies decreases in the order F > Cl > Br > I.