No Arabic abstract
High-temperature reactions widely exist in nature. However, they are difficult to be characterized either experimentally or computationally. The routinely used minimum energy path (MEP) model in computational modeling of chemical reactions is not justified to describe high-temperature reactions since high-energy structures are actively involved there. In this study, using CH4 decomposition on the Cu(111) surface as an example, we systematically compare MEP results with those obtained by explicitly sampling all relevant structures via ab initio molecular dynamics (AIMD) simulations at different temperatures. Interestingly, we find that, for reactions protected by a strong steric hindrance effect, the MEP is still effectively followed even at a temperature close to the Cu melting point. In contrast, without such a protection, the flexibility of surface Cu atoms can lead to a significant free energy barrier reduction at a high temperature. Accordingly, some conclusions about graphene growth mechanisms based on MEP calculations should be revisited. Physical insights provided by this study can deepen our understanding on high-temperature surface reactions.
The free energy profile of a reaction can be estimated in a molecular-dynamics approach by imposing a mechanical constraint along a reaction coordinate (RC). Many recent studies have shown that the temperature can greatly influence the path followed by the reactants. Here, we propose a practical way to construct the minimum energy path directly on the free energy surface (FES) at a given temperature. First, we follow the blue-moon ensemble method to derive the expression of the free energy gradient for a given RC. These derivatives are then used to find the actual minimum energy reaction path at finite temperature, in a way similar to the Intrinsic Reaction Path of Fukui on the potential energy surface [K Fukui J. Phys. Chem. 74, 4161 (1970)]. Once the path is know, one can calculate the free energy profile using thermodynamic integration. We also show that the mass-metric correction cancels for many types of constraints, making the procedure easy to use. Finally, the minimum free energy path at 300 K for the addition of the 1,1-dichlorocarbene to ethylene is compared with a path based on a simple one-dimensional reaction coordinate. A comparison is also given with the reaction path at 0 K.
We investigate the influence of a stochastically fluctuating step-barrier potential on bimolecular reaction rates by exact analytical theory and stochastic simulations. We demonstrate that the system exhibits a new resonant reaction behavior with rate enhancement if an appropriately defined fluctuation decay length is of the order of the system size. Importantly, we find that in the proximity of resonance the standard reciprocal additivity law for diffusion and surface reaction rates is violated due to the dynamical coupling of multiple kinetic processes. Together, these findings may have important repercussions on the correct interpretation of various kinetic reaction problems in complex systems, as, e.g., in biomolecular association or catalysis.
We study the elastic scattering of slow electrons by two-atomic molecule in the frame of non-overlapping atomic potentials model. The molecular continuum wave function is represented as a combination of a plane wave and two spherical s-waves, generated by the centers of atomic spheres. The asymptotic of this function determines in closed form the amplitude of elastic electron scattering. We show that this amplitude cannot be represented as a series of spherical functions. Therefore, it is impossible to use straightly the usual S-matrix methods to determine the scattering phases for non-spherical targets. We show that far from molecule the continuum wave function can be presented as an expansion in other than spherical orthonormal functions. The coefficients of this expansion determine the molecular scattering phases for non-spherical molecular systems. In such an approach, we calculate the Wigner times delay for slow electron scattered by two-atomic target.
The uniaxial negative thermal expansion in pentacene crystals along $a$ is a particularity in the series of the oligoacenes, and exeptionally large for a crystalline solid. Full x-ray structure analysis from 120 K to 413 K reveals that the dominant thermal motion is a libration of the rigid molecules about their long axes, modifying the intermolecular angle which describes the herringbone packing within the layers. This herringbone angle increases with temperature (by 0.3 -- 0.6$^{circ}$ per 100 K), and causes an anisotropic rearrangement of the molecules within the layers, i.e. an expansion in the $b$ direction, and a distinct contraction along $a$. Additionally, a larger herringbone angle improves the cofacial overlap between adjacent, parallel molecules, and thus enhances the attractive van der Waals forces.
There have existed for a long time a paradigm that TiO phases at ambient conditions are stable only if structural vacancies are available. Using an evolutionary algorithm, we perform an ab initio search of possible zero-temperature polymorphs of TiO in wide pressure interval. We obtain the Gibbs energy of the competing phases taking into account entropy via quasiharmonic approximation and build the pressure-temperature diagram of the system. We reveal that two vacancy-free hexagonal phases are the most stable at relatively low temperatures in a wide range of pressures. The transition between these phases takes place at 28 GPa. Only above 1290 K at ambient pressure the phases with vacancies (B1-derived) become stable. In particular, the high-pressure hexagonal phase is shown to have unusual electronic properties, with a pronounced pseudo-gap in the electronic spectrum. The comparison of DFT-GGA and GW calculations demonstrates that the account for many-body corrections significantly changes the electronic spectrum near the Fermi energy.