Long-range proton and hydroxide ion transfer dynamics at the water/CeO2 interface in the nanosecond regime: reactive molecular dynamics simulations and kinetic analysis.

The structural properties, dynamical behaviors, and ion transport phenomena at the interface between water and cerium oxide are investigated by reactive molecular dynamics (MD) simulations employing neural network potentials (NNPs). The NNPs are trained to reproduce density functional theory (DFT) results, and DFT-based MD (DFT-MD) simulations with enhanced sampling techniques and refinement schemes are employed to efficiently and systematically acquire training data that include diverse hydrogen-bonding configurations caused by proton hopping events. The water interfaces with two low-index surfaces of (111) and (110) are explored with these NNPs, and the structure and long-range proton and hydroxide ion transfer dynamics are examined with unprecedented system sizes and long simulation times. Various types of proton hopping events at the interface are categorized and analyzed in detail. Furthermore, in order to decipher the proton and hydroxide ion transport phenomena along the surface, a counting analysis based on the semi-Markov process is formulated and applied to the MD trajectories to obtain reaction rates by considering the transport as stochastic jump processes. Through this model, the coupling between hopping events, vibrational motions, and hydrogen bond networks at the interface are quantitatively examined, and the high activity and ion transport phenomena at the water/CeO2 interface are unequivocally revealed in the nanosecond regime.

On-the-fly kinetic Monte Carlo simulations with neural network potentials for surface diffusion and reaction.

We develop an adaptive scheme in the kinetic Monte Carlo simulations, where the adsorption and activation energies of all elementary steps, including the effects of other adsorbates, are evaluated "on-the-fly" by employing the neural network potentials. The configurations and energies evaluated during the simulations are stored for reuse when the same configurations are sampled in a later step. The present scheme is applied to hydrogen adsorption and diffusion on the Pd(111) and Pt(111) surfaces and the CO oxidation reaction on the Pt(111) surface. The effects of interactions between adsorbates, i.e., adsorbate-adsorbate lateral interactions, are examined in detail by comparing the simulations without considering lateral interactions. This study demonstrates the importance of lateral interactions in surface diffusion and reactions and the potential of our scheme for applications in a wide variety of heterogeneous catalytic reactions.

Catalytic ammonia synthesis on HY-zeolite-supported angstrom-size molybdenum cluster.

The development of new catalysts with high N2 activation ability is an effective approach for low-temperature ammonia synthesis. Herein, we report a novel angstrom-size molybdenum metal cluster catalyst for efficient ammonia synthesis. This catalyst is prepared by the impregnation of a molybdenum halide cluster complex with an octahedral Mo6 metal core on HY zeolite, followed by the removal of all the halide ligands by activation with hydrogen. In this activation, the size of the Mo6 cluster (ca. 7 Å) is almost retained. The resulting angstrom-size cluster shows catalytic activity for ammonia synthesis from N2 and H2, and the reaction proceeds continuously even at 200 °C under 5.0 MPa. DFT calculations suggest that N[triple bond, length as m-dash]N bond cleavage is promoted by the cooperation of the multiple molybdenum sites.

Collective bath coordinate mapping of "hierarchy" in hierarchical equations of motion.

The theory of hierarchical equations of motion (HEOM) is one of the standard methods to give exact evaluations of the dynamics as coupled to harmonic oscillator environments. However, the theory is numerically demanding due to its hierarchy, which is the set of auxiliary elements introduced to capture the non-Markovian and non-perturbative effects of environments. When system-bath coupling becomes relatively strong, the required computational resources and precision move beyond the regime that can be currently handled. This article presents a new representation of HEOM theory in which the hierarchy is mapped into a continuous space of a collective bath coordinate and several auxiliary coordinates as the form of the quantum Fokker-Planck equation. This representation gives a rigorous time evolution of the bath coordinate distribution and is more stable and efficient than the original HEOM theory, particularly when there is a strong system-bath coupling. We demonstrate the suitability of this approach to treat vibronic system models coupled to environments.

Generalization of the hierarchical equations of motion theory for efficient calculations with arbitrary correlation functions.

The hierarchical equations of motion (HEOM) theory is one of the standard methods to rigorously describe open quantum dynamics coupled to harmonic environments. Such a model is used to capture non-Markovian and non-perturbative effects of environments appearing in ultrafast phenomena. In the regular framework of the HEOM theory, the environment correlation functions are restricted to linear combinations of exponential functions. In this article, we present a new formulation of the HEOM theory including treatment of non-exponential correlation functions, which enables us to describe general environmental effects more efficiently and stably than the original theory and other generalizations. The library and its Python binding we developed to perform simulations based on our approach, named LibHEOM and PyHEOM, respectively, are provided as the supplementary material.

Modeling and analyzing a photo-driven molecular motor system: Ratchet dynamics and non-linear optical spectra.

A light-driven molecular motor system is investigated using a multi-state Brownian ratchet model described by a single effective coordinate with multiple electronic states in a dissipative environment. The rotational motion of the motor system is investigated on the basis of wavepacket dynamics. A current determined from the interplay between a fast photochemical isomerization (photoisomerization) process triggered by pulses and a slow thermal isomerization (thermalization) process arising from an overdamped environment is numerically evaluated. For this purpose, we employ the multi-state low-temperature quantum Smoluchowski equations that allow us to simulate the fast quantum electronic dynamics in the overdamped environment, where conventional approaches, such as the Zusman equation approach, fail to apply due to the positivity problem. We analyze the motor efficiency by numerically integrating the equations of motion for a rotator system driven by repeatedly impulsive excitations. When the time scales of the pulse repetition, photoisomerization, and thermalization processes are separated, the average rotational speed of the motor is determined by the time scale of thermalization. In this regime, the average rotational current can be described by a simple equation derived from a rate equation for the thermalization process. When laser pulses are applied repeatedly and the time scales of the photoisomerization and pulse repetition are close, the details of the photoisomerization process become important to analyze the entire rotational process. We examine the possibility of observing the photoisomerization and the thermalization processes associated with stationary rotating dynamics of the motor system by spectroscopic means, e.g., pump-probe, transient absorption, and two-dimensional electronic spectroscopy techniques.

Low-Temperature Quantum Fokker-Planck and Smoluchowski Equations and Their Extension to Multistate Systems.

Simulating electron-nucleus coupled dynamics poses a nontrivial challenge and an important problem in the investigation of ultrafast processes involving coupled electronic and vibrational dynamics. Because irreversibility of the system dynamics results from thermal activation and dissipation caused by the environment, in dynamical studies, it is necessary to include heat bath degrees of freedom in the total system. When the system dynamics involves high-energy electronic transitions, the environment is regarded to be in a low-temperature regime and we must treat it quantum mechanically. In this Article, we present rigorous and versatile approaches for investigating the dynamics of open systems with coupled electronic and vibrational degrees of freedom within a fully quantum mechanical framework. These approaches are based on a quantum Fokker-Planck equation and a quantum Smoluchowski equation employing a heat bath with an Ohmic spectral density, with non-Markovian low-temperature correction terms, and extensions of these equations to the case of multistate systems. The accuracy of these equations was numerically examined for a single-state Brownian system, while their applicability was examined for multistate double-well systems by comparing their results with those of the fewest-switch surface hopping and Ehrenfest methods with a classical Markovian Langevin force. Comparison of the transient absorption spectra obtained using these methods clearly reveals the importance of the quantum low-temperature correction terms. These equations allow us to treat nonadiabatic dynamics in an efficient way, while maintaining numerical accuracy. The C++ source codes that we developed, which allow for the treatment of the phase and coordinate space dynamics with any single-state or multistate potential forms, are provided as Supporting Information.

Probing photoisomerization processes by means of multi-dimensional electronic spectroscopy: The multi-state quantum hierarchical Fokker-Planck equation approach.

Photoisomerization in a system with multiple electronic states and anharmonic potential surfaces in a dissipative environment is investigated using a rigorous numerical method employing quantum hierarchical Fokker-Planck equations (QHFPEs) for multi-state systems. We have developed a computer code incorporating QHFPE for general-purpose computing on graphics processing units that can treat multi-state systems in phase space with any strength of diabatic coupling of electronic states under non-perturbative and non-Markovian system-bath interactions. This approach facilitates the calculation of both linear and nonlinear spectra. We computed Wigner distributions for excited, ground, and coherent states. We then investigated excited state dynamics involving transitions among these states by analyzing linear absorption and transient absorption processes and multi-dimensional electronic spectra with various values of heat bath parameters. Our results provide predictions for spectroscopic measurements of photoisomerization dynamics. The motion of excitation and ground state wavepackets and their coherence involved in the photoisomerization were observed as the profiles of positive and negative peaks of two-dimensional spectra.

Analysis of 2D THz-Raman spectroscopy using a non-Markovian Brownian oscillator model with nonlinear system-bath interactions.

We explore and describe the roles of inter-molecular vibrations employing a Brownian oscillator (BO) model with linear-linear (LL) and square-linear (SL) system-bath interactions, which we use to analyze two-dimensional (2D) THz-Raman spectra obtained by means of molecular dynamics (MD) simulations. In addition to linear infrared absorption (1D IR), we calculated 2D Raman-THz-THz, THz-Raman-THz, and THz-THz-Raman signals for liquid formamide, water, and methanol using an equilibrium non-equilibrium hybrid MD simulation. The calculated 1D IR and 2D THz-Raman signals are compared with results obtained from the LL+SL BO model applied through use of hierarchal Fokker-Planck equations with non-perturbative and non-Markovian noise. We find that all of the qualitative features of the 2D profiles of the signals obtained from the MD simulations are reproduced with the LL+SL BO model, indicating that this model captures the essential features of the inter-molecular motion. We analyze the fitted 2D profiles in terms of anharmonicity, nonlinear polarizability, and dephasing time. The origins of the echo peaks of the librational motion and the elongated peaks parallel to the probe direction are elucidated using optical Liouville paths.