Ricardo Pérez de Tudela

Stress Test for Quantum Dynamics Approximations: Deep Tunneling in the Muonium Exchange Reaction D + HMu → DMu + H

Quantum effects play a crucial role in chemical reactions involving light atoms at low temperatures, especially when a light particle is exchanged between two heavier partners. Different theoretical methodologies have been developed in the last decades attempting to describe zero-point energy and tunneling effects without abandoning a classical or semiclassical framework. In this work, we have chosen the D + HMu → DMu + H reaction as a stress test system for three well-established methods: two representative versions of transition state theory (TST), canonical variational theory and semiclassical instanton, and ring polymer molecular dynamics (RPMD). These calculations will be compared with accurate quantum mechanical results. Despite its apparent simplicity, the exchange of the extremely light muonium atom (0.114 u) becomes a most challenging reaction for conventional methods. The main result of this work is that RPMD provides an overall better performance than TST-based methods for such a demanding reaction. RPMD might well turn out to be a useful tool beyond TST applicability.

Quantum features of a barely bound molecular dopant: Cs2(3Σu) in bosonic helium droplets of variable size.

We present in this work the study of small (4)He(N)-Cs(2)((3)Σ(u)) aggregates (2 ≤ N ≤ 30) through combined variational, diffusion Monte Carlo (DMC), and path integral Monte Carlo (PIMC) calculations. The full surface is modeled as an addition of He-Cs(2) interactions and He-He potentials. Given the negligible strength and large range of the He-Cs(2) interaction as compared with the one for He-He, a propensity of the helium atoms to pack themselves together, leaving outside the molecular dopant is to be expected. DMC calculations determine the onset of helium gathering at N = 3. To analyze energetic and structural properties as a function of N, PIMC calculations with no bosonic exchange, i.e., Boltzmann statistics, at low temperatures are carried out. At T = 0.1 K, although acceptable one-particle He-Cs(2) distributions are obtained, two-particle He-He distributions are not well described, indicating that the proper symmetry should be taken into account. PIMC distributions at T = 1 K already compare well with DMC ones and show minor exchange effects, although binding energies are still far from having converged in terms of the number of quantum beads. As N increases, the He-He PIMC pair correlation function shows a clear tendency to coincide with the experimental boson-liquid helium one at that temperature. It supports the picture of a helium droplet which carries the molecular impurity on its surface, as found earlier for other triplet dimers.

Energetics and solvation structure of a dihalogen dopant (I2) in 4He clusters

The energetics and structure of small HeNI2 clusters are analyzed as the size of the system changes, with N up to 38. The full interaction between the I2 molecule and the He atoms is based on analytical ab initio He-I2 potentials plus the He-He interaction, obtained from first-principle calculations. The most stable structures, as a function of the number of solvent He atoms, are obtained by employing an evolutionary algorithm and compared with CCSD(T) and MP2 ab initio computations. Further, the classical description is completed by explicitly including thermal corrections and quantum features, such as zero-point-energy values and spatial delocalization. From quantum PIMC calculations, the binding energies and radial/angular probability density distributions of the thermal equilibrium state for selected-size clusters are computed at a low temperature. The sequential formation of regular shell structures is analyzed and discussed for both classical and quantum treatments.

Examination of the Feynman–Hibbs Approach in the Study of NeN-Coronene Clusters at Low Temperatures

Feynman-Hibbs (FH) effective potentials constitute an appealing approach for investigations of many-body systems at thermal equilibrium since they allow us to easily include quantum corrections within standard classical simulations. In this work we apply the FH formulation to the study of NeN-coronene clusters (N = 1-4, 14) in the 2-14 K temperature range. Quadratic (FH2) and quartic (FH4) contributions to the effective potentials are built upon Ne-Ne and Ne-coronene analytical potentials. In particular, a new corrected expression for the FH4 effective potential is reported. FH2 and FH4 cluster energies and structures-obtained from energy optimization through a basin-hopping algorithm as well as classical Monte Carlo simulations-are reported and compared with reference path integral Monte Carlo calculations. For temperatures T > 4 K, both FH2 and FH4 potentials are able to correct the purely classical calculations in a consistent way. However, the FH approach fails at lower temperatures, especially the quartic correction. It is thus crucial to assess the range of applicability of this formulation and, in particular, to apply the FH4 potentials with great caution. A simple model of N isotropic harmonic oscillators allows us to propose a means of estimating the cutoff temperature for the validity of the method, which is found to increase with the number of atoms adsorbed on the coronene molecule.

Publisher’s Note: Acid Dissociation in HCl-Water Clusters is Temperature Dependent and Cannot be Detected Based on Dipole Moments [Phys. Rev. Lett. 119 , 223001 (2017)]

This corrects the article DOI: 10.1103/PhysRevLett.119.223001.