No Arabic abstract
Clathrate hydrates hold considerable promise as safe and economical materials for hydrogen storage. Here we present a quantum mechanical study of H$_2$ and D$_2$ diffusion through a hexagonal face shared by two large cages of clathrate hydrates over a wide range of temperatures. Path integral molecular dynamics simulations are used to compute the free-energy profiles for the diffusion of H$_2$ and D$_2$ as a function of temperature. Ring polymer molecular dynamics rate theory, incorporating both exact quantum statistics and approximate quantum dynamical effects, is utilized in the calculations of the H$_2$ and D$_2$ diffusion rates in a broad temperature interval. We find that the shape of the quantum free-energy profiles and their height relative to the classical free energy barriers at a given temperature, as well as the rate of diffusion, are profoundly affected by competing quantum effects: above 25 K, zero-point energy (ZPE) perpendicular to the reaction path for diffusion between cavities decreases the quantum rate compared to the classical rate, whereas at lower temperatures tunneling outcompetes the ZPE and as result the quantum rate is greater than the classical rate.
Water management in a hydrogen polymer electrolyte membrane (PEM) fuel cell is critical for performance. The impact of thermal conductivity and water vapor diffusion coefficients in a gas diffusion layer (GDL) has been studied by a mathematical model. The fraction of product water that is removed in the vapour phase through the GDL as a function of GDL properties and operating conditions has been calculated and discussed. Furthermore, the current model enables identification of conditions when condensation occurs in each GDL component and calculation of temperature gradient across the interface between different layers, providing insight into the overall mechanism of water transport in a given cell design. Water transport mode and condensation conditions in the GDL components depend on the combination of water vapor diffusion coefficients and thermal conductivities of the GDL components. Different types of GDL and water removal scenarios have been identified and related to experimentally-determined GDL properties.
A fundamental question in the study of chemical reactions is how reactions proceed at a collision energy close to absolute zero. This question is no longer hypothetical: quantum degenerate gases of atoms and molecules can now be created at temperatures lower than a few tens of nanoKelvin. In this work we consider the benchmark ultracold reaction between, the most-celebrated ultracold molecule, KRb and K. For the first time we map out an accurate ab initio ground state potential energy surface of the KRbK complex in full dimensionality and report numerically exact quantum-mechanical reaction dynamics. The distribution of rotationally resolved rates is shown to be Poissonian. An analysis of the hyperspherical adiabatic potential curves explains this statistical character revealing a chaotic distribution for the short-range collision complex that plays a key role in governing the reaction outcome. We compare this with a lighter system with a smaller density of states (here the LiYbLi trimer) which displays random, and not chaotic, behavior.
We introduce a new optical tool - a two-dimensional optical centrifuge, capable of aligning molecules in extreme rotational states. Unlike the conventional centrifuge, which confines the molecules in the plane of their rotation, its two-dimensional version aligns the molecules along a well-defined axis, similarly to the effect of a single linearly polarized laser pulse, but at a much higher level of rotational excitation. The increased robustness of ultra-high rotational states with respect to collisions results in a longer life time of the created alignment in dense media, offering new possibilities for studying and utilizing aligned molecular ensembles under ambient conditions.
We study the formation of C$_{18}$H and C$_{18}$H$_2$ by irradiating a cyclo[$18$]carbon molecule with atomic and molecular hydrogen at impact energy, $E$, in the range of 0.5-25 eV. We utilize the density-functional tight-binding method to perform molecular dynamics simulations to emulate the interaction of a carbon ring when colliding with atomic or molecular hydrogen. From our results, the formation of the C$_{18}$H molecules is likely to occur upon irradiating by H atoms at $E < 10$ eV and by H$_2$ molecules at $2 < E < 15$ eV center of mass energy. Formation of C$_{18}$H$_2$ molecules is only observed at around $E = 2$ eV. Our results show that the absorption of hydrogen is more prone in atomic than in molecular hydrogen atmosphere. Thus, we find that the probability of physio-absorption reaches up to 80 % for atomic projectiles with $E < 5$ eV but only up to 10 % for the molecular ones. Our analysis shows that the deformation of the carbon ring due to the hydrogen bonding produces transition from $sp$ to $sp^2$ hybridization. The angle between the carbon atoms at the locations near to the H bond in the resulting ring is not 120$^o$ but instead 110$^o$ degrees. No molecular fragmentation of the C$_{18}$ ring is observed.
We calculate the linear and non-linear susceptibilities of periodic longitudinal chains of hydrogen dimers with different bond-length alternations using a diffusion quantum Monte Carlo approach. These quantities are derived from the changes in electronic polarization as a function of applied finite electric field - an approach we recently introduced and made possible by the use of a Berry-phase, many-body electric-enthalpy functional. Calculated susceptibilities and hyper-susceptibilities are found to be in excellent agreement with the best estimates available from quantum chemistry - usually extrapolations to the infinite-chain limit of calculations for chains of finite length. It is found that while exchange effects dominate the proper description of the susceptibilities, second hyper-susceptibilities are greatly affected by electronic correlations. We also assess how different approximations to the nodal surface of the many-body wavefunction affect the accuracy of the calculated susceptibilities.