No Arabic abstract
Fractionation of isotopes among distinct molecules or phases is a quantum effect which is often exploited to obtain insights on reaction mechanisms, biochemical, geochemical and atmospheric phenomena. Accurate evaluation of isotope ratios in atomistic simulations is challenging, because one needs to perform a thermodynamic integration with respect to the isotope mass, along with time-consuming path integral calculations. By re-formulating the problem as a particle exchange in the ring polymer partition function, we derive new estimators giving direct access to the differential partitioning of isotopes, which can simplify the calculations by avoiding thermodynamic integration. We demonstrate the efficiency of these estimators by applying them to investigate the isotope fractionation ratios in the gas-phase Zundel cation, and in a few simple hydrocarbons.
At room temperature, the quantum contribution to the kinetic energy of a water molecule exceeds the classical contribution by an order of magnitude. The quantum kinetic energy (QKE) of a water molecule is modulated by its local chemical environment and leads to uneven partitioning of isotopes between different phases in thermal equilibrium, which would not occur if the nuclei behaved classically. In this work, we use ab initio path integral simulations to show that QKEs of the water molecules and the equilibrium isotope fractionation ratios of the oxygen and hydrogen isotopes are sensitive probes of the hydrogen bonding structures in aqueous ionic solutions. In particular, we demonstrate how the QKE of water molecules in path integral simulations can be decomposed into translational, rotational and vibrational degrees of freedom, and use them to determine the impact of solvation on different molecular motions. By analyzing the QKEs and isotope fractionation ratios, we show how the addition of the Na$^+$, Cl$^-$ and HPO$_4^{2-}$ ions perturbs the competition between quantum effects in liquid water and impacts their local solvation structures.
Path reweighting is a principally exact method to estimate dynamic properties from biased simulations - provided that the path probability ratio matches the stochastic integrator used in the simulation. Previously reported path probability ratios match the Euler-Maruyama scheme for overdamped Langevin dynamics. Since MD simulations use Langevin dynamics rather than overdamped Langevin dynamics, this severely impedes the application of path reweighting methods. Here, we derive the path probability ratio $M_L$ for Langevin dynamics propagated by a variant of the Langevin Leapfrog integrator. This new path probability ratio allows for exact reweighting of Langevin dynamics propagated by this integrator. We also show that a previously derived approximate path probability ratio $M_{mathrm{approx}}$ differs from the exact $M_L$ only by $mathcal{O}(xi^4Delta t^4)$, and thus yields highly accurate dynamic reweighting results. ($Delta t$ is the integration time step, $xi$ is the collision rate.) The results are tested and the efficiency of path-reweighting is explored using butane as an example.
Revelation of chlorine and bromine isotope fractionation of halogenated organic compounds (HOCs) in electron ionization mass spectrometry (EI-MS) is crucial for compound-specific chlorine/bromine isotope analysis (CSIA-Cl/Br) using gas chromatography EI-MS (GC-EI-MS). This study systematically investigated chlorine/bromine isotope fractionation in EI-MS of HOCs including 12 organochlorines and 5 organobromines using GC-double focus magnetic-sector high resolution MS (GC-DFS-HRMS). Chlorine/bromine isotope fractionation behaviors of the HOCs in EI-MS showed varied isotope fractionation patterns and extents depending on compounds. Besides, isotope fractionation patterns and extents varied at different EI energies, demonstrating potential impacts of EI energy on the chlorine/bromine isotope fractionation. Hypotheses of inter-ion and intra-ion isotope fractionations were applied to interpreting the isotope fractionation behaviors. The inter-ion and intra-ion isotope fractionations counteractively contributed to the apparent isotope ratio for a certain dehalogenated product ion. The isotope fractionation mechanisms were tentatively elucidated on basis of the quasi-equilibrium theory. In the light of the findings of this study, isotope ratio evaluation scheme using complete molecular ions and the EI source with sufficient stable EI energies may be helpful to achieve optimal precision and accuracy of CSIA-Cl/Br data. The method and results of this study can help to predict isotope fractionation of HOCs during dehalogenation processes and further to reveal the dehalogenation pathways.
Path integral Monte Carlo approach is used to study the coupled quantum dynamics of the electron and nuclei in hydrogen molecule ion. The coupling effects are demonstrated by comparing differences in adiabatic Born--Oppenheimer and non-adiabatic simulations, and inspecting projections of the full three-body dynamics onto adiabatic Born--Oppenheimer approximation. Coupling of electron and nuclear quantum dynamics is clearly seen. Nuclear pair correlation function is found to broaden by 0.040 a_0 and average bond length is larger by 0.056 a_0. Also, non-adiabatic correction to the binding energy is found. Electronic distribution is affected less, and therefore, we could say that the adiabatic approximation is better for the electron than for the nuclei.
We investigate the gas-phase and grain-surface chemistry in the inner 30 AU of a typical protoplanetary disk using a new model which calculates the gas temperature by solving the gas heating and cooling balance and which has an improved treatment of the UV radiation field. We discuss inner-disk chemistry in general, obtaining excellent agreement with recent observations which have probed the material in the inner regions of protoplanetary disks. We also apply our model to study the isotopic fractionation of carbon. Results show that the fractionation ratio, 12C/13C, of the system varies with radius and height in the disk. Different behaviour is seen in the fractionation of different species. We compare our results with 12C/13C ratios in the Solar System comets, and find a stark contrast, indicative of reprocessing.