Search

A Smolyak algorithm adapted to system-bath separation is proposed for rigorous quantum simulations. This technique combines a sparse grid method with the system-bath concept in a specific configuration without limitations on the form of the Hamiltonian, thus achieving a highly efficient convergence of the excitation transitions for the "system" part. Our approach provides a general way to overcome the perennial convergence problem for the standard Smolyak algorithm and enables the simulation of floppy molecules with more than a hundred degrees of freedom.The efficiency of the present method is illustrated on the simulation of H$_2$ caged in an sII clathrate hydrate including two kinds of cage modes. The transition energies are converged by increasing the number of normal modes of water molecules. Our results confirm the triplet splittings of both translational and rotational ($j=1$) transitions of the H$_2$ molecule. Furthermore, they show a slight increase of the translational transitions with respect to the ones in a rigid cage., Comment: 5 pages, 3 figures

The reaction between D$^+$ and H$_2$ plays an important role in astrochemistry at low temperatures and also serves as a prototype for simple ion-molecule reaction. Its ground $\tilde{X}~^1A^{\prime}$ state has a very small thermodynamic barrier (up to 1.8$\times 10^{-2}$ eV) and the reaction proceeds through the formation of an intermediate complex lying within the potential well of depth of at least 0.2 eV thus representing a challenge for dynamical studies. In the present work, we analyze the title reaction within the temperature range of 20 $-$ 100 K by means of ring polymer molecular dynamics (RPMD) and quasi-classical trajectory (QCT) methods over the full-dimensional global potential energy surface developed by Aguado et al. [A. Aguado, O. Roncero, C. Tablero, C. Sanz, and M. Paniagua, J. Chem. Phys., 2000, 112, 1240]. The computed thermal RPMD and QCT rate coefficients are found to be almost independent of temperature and fall within the range of 1.34 $-$ 2.01$\times$10$^{-9}$ cm$^3$ s$^{-1}$. They are also in a very good agreement with the previous time-independent quantum mechanical and statistical quantum method calculations. Furthermore, we observe that the choice of asymptotic separation distance between the reactants can markedly alter the rate coefficient in the low temperature regime (20 $-$ 50 K). Therefore it is of utmost importance to correctly assign the value of this parameter for dynamical studies, particularly at very low temperatures of astrochemical importance. We finally conclude that experimental rate measurements for the title reaction are highly desirable in future.

The HDO/H2O ratio is a powerful diagnostic to understand the evolution of water from the first stages of star formation to the formation of planets and comets. Our aim is to determine precisely the abundance distribution of HDO towards the low-mass protostar IRAS16293-2422 and learn more about the water formation mechanisms by determining the HDO/H2O abundance ratio. A spectral survey of the source IRAS16293-2422 was carried out in the framework of the CHESS Herschel Key program with the HIFI instrument, allowing the detection of numerous HDO lines. Other transitions have been observed previously with ground-based telescopes. The spherical Monte Carlo radiative transfer code RATRAN was used to reproduce the observed line profiles of HDO by assuming an abundance jump. To determine the H2O abundance throughout the envelope, a similar study was made of the H2-18O observed lines, as the H2O main isotope lines are contaminated by the outflows. We derive an inner HDO abundance of 1.7e-7 and an outer HDO abundance of 8e-11. To reproduce the HDO absorption lines, it is necessary to add an absorbing layer in front of the envelope. It may correspond to a water-rich layer created by the photodesorption of the ices at the edges of the molecular cloud. The HDO/H2O ratio is ~1.4-5.8% in the hot corino whereas it is ~0.2-2.2% in the outer envelope. It is estimated at ~4.8% in the added absorbing layer. Although it is clearly higher than the cosmic D/H abundance, the HDO/H2O ratio remains lower than the D/H ratio derived for other deuterated molecules observed in the same source. The similarity of the ratios derived in the hot corino and in the added absorbing layer suggests that water formed before the gravitational collapse of the protostar, contrary to formaldehyde and methanol, which formed later once the CO molecules had depleted on the grains., Comment: 12 pages, 11 figures