Bond making and bond breaking in molecular dynamics

The theoretical concept of a potential energy surface, so important for the accepted pictures of molecular structure and bonding, plays a key role in most current methods of molecular reaction dynamics. However, the lack of quality potential energy surfaces for anything but the smallest of systems and the great expense and difficulty in producing several surfaces and their nonadiabatic coupling terms with sufficient accuracy is becoming a hindrance to accurate reaction dynamics. The electron nuclear dynamics (END) theory provides a new way to address this problem by circumventing the use of potential energy surfaces while accounting fully for nonadiabatic couplings. Preliminary results for the prototypical low-energy reactive collision between the hydrogen molecular ion and the hydrogen molecule are given for both electron transfer and chemical exchange. © 1996 John Wiley & Sons, Inc.     

利用基于直接动力学的轨线面跳跃方法处理非绝热过程

势能面交叉引起的非绝热过程广泛存在于光化学和光物理中。对这一过程进行描述是理论化学的重要挑战之一。非绝热过程涉及原子核与电子之间的耦合运动,因此量子化学的基本假设之一"玻恩-奥本海默"近似被打破,所以对其进行描述需要发展新的动力学理论方法。在这些方法中,Tully发展的最少轨线面跳跃方法凭借易于程序化、便于计算等优点已经发展成为处理非绝热问题的主要动力学方法之一。其中原子核以经典的方式在单一势能面上进行演化,电子以量子的方式沿着同一轨线进行演化。在整个演化过程中,非绝热跃迁通过轨线在不同势能面间的跃迁来描述,其中跳跃发生的几率与电子的演化有关。如果将该方法与从头算直接动力学相结合,可以在全原子水平上研究实际分子体系的非绝热动力学,给出其激发态寿命、非绝热动力学中分子的主要运动方式、反应通道以及分支比等重要信息。本文旨在讨论最少面跳跃直接动力学方法研究非绝热问题的一些进展,包括动力学基本理论,特别关注将最少面跳跃方法和直接动力学结合的数值实现细节,同时讨论该方法在研究实际体系当中的一些应用,并对轨线面跳跃方法下一步发展的一些方向进行合理的展望。     

Role of Rydberg states in the photochemical dynamics of ethylene

We use the ab initio multiple spawning method with potential energy surfaces and nonadiabatic coupling vectors computed from multistate multireference perturbation theory (MSPT2) to follow the dynamics of ethylene after photoexcitation. We introduce an analytic formulation for the nonadiabatic coupling vector in the context of MSPT2 calculations. We explicitly include the low-lying 3s Rydberg state which has been neglected in previous ab initio molecular dynamics studies of this process. We find that although the 3s Rydberg state lies below the optically bright ππ* state, little population gets trapped on this state. Instead, the 3s Rydberg state is largely a spectator in the photodynamics, with little effect on the quenching mechanism or excited state lifetime. We predict the time-resolved photoelectron spectrum for ethylene and point out the signature of Rydberg state involvement that should be easily observed.     

First-Principle Calculations of the Adsorption, Dissociation and Diffusion of Hydrogen on the Mg(0001) Surface

First-principle pseudopotential plane wave calculations and the Nudged Elastic Band method based on density functional theory (DFT) have been used in this article to study the dissociation of molecular hydrogen on a Mg(0001) surface and the subsequent diffusion of atomic hydrogen into the magnesium substrate. First, the dissociation pathway of H2 and the relative barrier were investigated. It was shown that physical adsorption rather than chemisorption of molecular hydrogen was observed in the calculation of the dissociation process of molecular hydrogen. Also, the diffusion process of atomic hydrogen on Mg(0001) was presented. The surface effect, which affected the diffusion of hydrogen obviously, was observed. Finally, comparing the values of the activation energies for the steps of dissociation, diffusion, and desorption, our calculation further showed that the dissociation of H2 and the desorption of hydride were the rate-limiting steps.     

H在Mg(0001)表面吸附、解离和扩散的第一性原理研究

采用基于密度泛函的第一性原理方法, 同时结合Nudged Elastic Band方法, 系统研究了H2分子和H原子在Mg(0001)表面的吸附过程. 给出了H2分子的解离路径和势垒, 结果表明H2分子的吸附过程中仅存在物理吸附; 在给出H原子在Mg(0001)表面的吸附势能面的基础上, 进一步研究了H原子在Mg(0001)表面及体内的扩散过程. 计算发现, Mg(0001) slab存在表面效应, 且对H原子的表面扩散影响较明显. 在此基础上, 通过比较解离、扩散和放氢环节的激活能数据, 为H2分子的解离和氢化物的放氢过程是速控步骤这一结论提供了理论支持.     

Quantum and dynamical effects of proton donor-acceptor vibrational motion in nonadiabatic proton-coupled electron transfer reactions

This paper presents a general theoretical formulation for proton-coupled electron transfer (PCET) reactions. The solute is represented by a multistate valence bond model, and the active electrons and transferring proton(s) are treated quantum mechanically. This formulation enables the classical or quantum mechanical treatment of the proton donor-acceptor vibrational mode, as well as the dynamical treatment of the proton donor-acceptor mode and the solvent. Nonadiabatic rate expressions are presented for PCET reactions in a number of well-defined limits for both dielectric continuum and molecular representations of the environment. The dynamical rate expressions account for correlations between the fluctuations of the proton donor-acceptor distance and the nonadiabatic PCET coupling. The quantities in the rate expressions can be calculated with a dielectric continuum model or a molecular dynamics simulation of the full system. The significance of the quantum and dynamical effects of the proton donor-acceptor mode is illustrated with applications to model PCET systems.     

Dynamics of charge transfer states on metal surfaces: the competition between reactivity and quenching

The dynamics of excited states of adsorbates on surfaces caused by charge transfer is studied. Both negative and positive charge transfer processes are possible. In particular we are interested in positive charge transfer from a metal surface to molecular or atomic oxygen adsorbed on the surface. Once the negatively charged oxygen on the surface loses an electron it becomes chemically activated. The ability of this species to react depends on the quenching time or back transfer. The analysis of these processes is based on a set of diabatic potential energy surfaces each representing a different charged oxygen species. The dynamics is followed by solving the multichannel time-dependent Schr?dinger equation or Liouville von Neumann equation. Due to the nonadiabatic character of these reactions large isotope effects are predicted.     

Trajectory-based solution of the nonadiabatic quantum dynamics equations: an on-the-fly approach for molecular dynamics simulations

The non-relativistic quantum dynamics of nuclei and electrons is solved within the framework of quantum hydrodynamics using the adiabatic representation of the electronic states. An on-the-fly trajectory-based nonadiabatic molecular dynamics algorithm is derived, which is also able to capture nuclear quantum effects that are missing in the traditional trajectory surface hopping approach based on the independent trajectory approximation. The use of correlated trajectories produces quantum dynamics, which is in principle exact and computationally very efficient. The method is first tested on a series of model potentials and then applied to study the molecular collision of H with H(2) using on-the-fly TDDFT potential energy surfaces and nonadiabatic coupling vectors.     

A chemical route to control molecular mobility on graphene

Controlling the motion of nanoscale building blocks on chemically contaminated or modified substrates is a bottleneck in bottom-up approaches to develop high performance Nanoscale Electro-Mechanical Systems (NEMS). Nevertheless, how such modification of a substrate surface affects the mobility of an admolecule has not been well understood. Here, we employ molecular dynamics (MD) simulations to study the surface diffusion of a C(60) admolecule on both pure and hydrogenated graphene. By changing temperature and hydrogen coverage, we obtain a diagram which shows evidently the existence of three distinct regimes of the surface diffusion of a C(60) admolecule, that is, superdiffusion, normal diffusion, and subdiffusion. Surprisingly, our simulations also show that minute hydrogenation on graphene leads to a giant reduction in molecular mobility. A theoretical model which takes into account the effects of both random traps and barriers is developed to predict the relation between the diffusion coefficient and the temperature and hydrogen coverage. The model predictions are in good agreement with our molecular dynamics simulation results.     

Semiclassical theory of electronically nonadiabatic chemical dynamics: incorporation of the Zhu-Nakamura theory into the frozen Gaussian propagation method

The title theory is developed by combining the Herman-Kluk semiclassical theory for adiabatic propagation on single potential-energy surface and the semiclassical Zhu-Nakamura theory for nonadiabatic transition. The formulation with use of natural mathematical principles leads to a quite simple expression for the propagator based on classical trajectories and simple formulas are derived for overall adiabatic and nonadiabatic processes. The theory is applied to electronically nonadiabatic photodissociation processes: a one-dimensional problem of H2+ in a cw (continuous wave) laser field and a two-dimensional model problem of H2O in a cw laser field. The theory is found to work well for the propagation duration of several molecular vibrational periods and wide energy range. Although the formulation is made for the case of laser induced nonadiabatic processes, it is straightforwardly applicable to ordinary electronically nonadiabatic chemical dynamics.     

Influence of thermal fluctuations on interfacial electron transfer in functionalized TiO2 semiconductors

The influence of thermal fluctuations on the dynamics of interfacial electron transfer in sensitized TiO2-anatase semiconductors is investigated by combining ab initio DFT molecular dynamics simulations and quantum dynamics propagation of transient electronic excitations. It is shown that thermal nuclear fluctuations speed up the underlying interfacial electron transfer dynamics by introducing nonadiabatic transitions between electron acceptor states, localized in the vicinity of the photoexcited adsorbate, and delocalized states extended throughout the semiconductor material, creating additional relaxation pathways for carrier diffusion. Furthermore, it is shown that room-temperature thermal fluctuations reduce the anisotropic character of charge diffusion along different directions in the anatase crystal and make similar the rates for electron injection from adsorbate states of different character. The reported results are particularly relevant to the understanding of temperature effects on surface charge separation mechanisms in molecular-based photo-optic devices.     

Pneumatochemical impedance spectroscopy. 2. Dynamics of hydrogen sorption by metals

In this paper, pneumatochemical impedance spectroscopy is used to analyze multistep reaction mechanisms such as those observed in solid solution domains of LaNi5-H2(g) systems. It is shown that hydrogen sorption is a two-step mechanism including (i) dissociative surface chemisorption of molecular hydrogen and (ii) atomic hydrogen bulk transport by diffusion. Data fitting of experimental transfer functions with model equations yields the value of the kinetic parameter associated with each individual reaction step, i.e., surface sorption resistances and hydrogen bulk diffusion coefficients. The technique is used to follow the activation procedure of the sample as well as the degradation of sorption properties in oxygen-containing hydrogen atmospheres. A decrease in sorption kinetics is attributed to surface oxidation, whereas bulk properties remain unchanged. The perspectives offered by the technique which potentially can be used to optimize surface and bulk composition of IMC for increased sorption rates are discussed.     

Measuring nonadiabaticity of molecular quantum dynamics with quantum fidelity and with its efficient semiclassical approximation

We propose to measure nonadiabaticity of molecular quantum dynamics rigorously with the quantum fidelity between the Born-Oppenheimer and fully nonadiabatic dynamics. It is shown that this measure of nonadiabaticity applies in situations where other criteria, such as the energy gap criterion or the extent of population transfer, fail. We further propose to estimate this quantum fidelity efficiently with a generalization of the dephasing representation to multiple surfaces. Two variants of the multiple-surface dephasing representation (MSDR) are introduced, in which the nuclei are propagated either with the fewest-switches surface hopping or with the locally mean field dynamics (LMFD). The LMFD can be interpreted as the Ehrenfest dynamics of an ensemble of nuclear trajectories, and has been used previously in the nonadiabatic semiclassical initial value representation. In addition to propagating an ensemble of classical trajectories, the MSDR requires evaluating nonadiabatic couplings and solving the Schro?dinger (or more generally, the quantum Liouville-von Neumann) equation for a single discrete degree of freedom. The MSDR can be also used in the diabatic basis to measure the importance of the diabatic couplings. The method is tested on three model problems introduced by Tully and on a two-surface model of dissociation of NaI.     

Ab Initio Nonadiabatic Dynamics of Semiconductor Materials via Surface Hopping Method

Photoinduced carrier dynamic processes are without doubt the main driving force responsible for the efficient performance of semiconductor nanomaterials in applications like photoconversion and photonics. Nevertheless, establishing theoretical insights into these processes is computationally challenging owing to the multiple factors involved in the processes, namely reaction rate, material surface area, material composition etc. Modelling of photoinduced carrier dynamic processes can be performed via nonadiabatic molecular dynamics (NA-MD) methods, which are methods specifically designed to solve the time-dependent Schrodinger equation with the inclusion of nonadiabatic couplings. Among NA-MD methods, surface hopping methods have been proven to be a mighty tool to mimic the competitive nonadiabatic processes in semiconductor nanomaterials, a worth noticing feature is its exceptional balance between accuracy and computational cost. Consequently, surface hopping is the method of choice for modelling ultrafast dynamics and more complex phenomena like charge separation in Janus transition metal dichalcogenides-based van der Waals heterojunction materials. Covering latest stateof-the-art numerical simulations along with experimental results in the field, this review aims to provide a basic understanding of the tight relation between semiconductor nanomaterials and the proper simulation of their properties via surface hopping methods. Special stress is put on emerging state-ot-the-art techniques. By highlighting the challenge imposed by new materials, we depict emerging creative approaches, including high-level electronic structure methods and NA-MD methods to model nonadiabatic systems with high complexity.     

Exploring the role of decoherence in condensed-phase nonadiabatic dynamics: a comparison of different mixed quantum/classical simulation algorithms for the excited hydrated electron

Mixed quantum/classical (MQC) molecular dynamics simulation has become the method of choice for simulating the dynamics of quantum mechanical objects that interact with condensed-phase systems. There are many MQC algorithms available, however, and in cases where nonadiabatic coupling is important, different algorithms may lead to different results. Thus, it has been difficult to reach definitive conclusions about relaxation dynamics using nonadiabatic MQC methods because one is never certain whether any given algorithm includes enough of the necessary physics. In this paper, we explore the physics underlying different nonadiabatic MQC algorithms by comparing and contrasting the excited-state relaxation dynamics of the prototypical condensed-phase MQC system, the hydrated electron, calculated using different algorithms, including: fewest-switches surface hopping, stationary-phase surface hopping, and mean-field dynamics with surface hopping. We also describe in detail how a new nonadiabatic algorithm, mean-field dynamics with stochastic decoherence (MF-SD), is to be implemented for condensed-phase problems, and we apply MF-SD to the excited-state relaxation of the hydrated electron. Our discussion emphasizes the different ways quantum decoherence is treated in each algorithm and the resulting implications for hydrated-electron relaxation dynamics. We find that for three MQC methods that use Tully's fewest-switches criterion to determine surface hopping probabilities, the excited-state lifetime of the electron is the same. Moreover, the nonequilibrium solvent response function of the excited hydrated electron is the same with all of the nonadiabatic MQC algorithms discussed here, so that all of the algorithms would produce similar agreement with experiment. Despite the identical solvent response predicted by each MQC algorithm, we find that MF-SD allows much more mixing of multiple basis states into the quantum wave function than do other methods. This leads to an excited-state lifetime that is longer with MF-SD than with any method that incorporates nonadiabatic effects with the fewest-switches surface hopping criterion.     

“Superpermeability” and “pumping” of atomic hydrogen through palladium membranes

Permeation of atomic as well as molecular hydrogen through palladium membranes has been investigated experimentally in the temperature range from room temperature to 200 °C and at a higher incident flux of hydrogen atoms on palladium surface than in previous studies. The results demonstrate that phenomena of ‘superpermeability’ and ‘pumping’ of atomic gases through metal membranes are of a common nature. A theoretical model based on chemical thermodynamics and diffusion theory adequately describes the quantitative relationships observed in experiments. It was found that permeability of atomic hydrogen depends strongly on the magnitude of surface incident flux and membrane temperature.     

An analysis of the accuracy of an initial value representation surface hopping wave function in the interaction and asymptotic regions

The behavior of an initial value representation surface hopping wave function is examined. Since this method is an initial value representation for the semiclassical solution of the time independent Schrodinger equation for nonadiabatic problems, it has computational advantages over the primitive surface hopping wave function. The primitive wave function has been shown to provide transition probabilities that accurately compare with quantum results for model problems. The analysis presented in this work shows that the multistate initial value representation surface hopping wave function should approach the primitive result in asymptotic regions and provide transition probabilities with the same level of accuracy for scattering problems as the primitive method.     

Ab initio nonadiabatic dynamics study of ultrafast radiationless decay over conical intersections illustrated on the Na3F cluster

We present a theoretical approach for the ultrafast nonadiabatic dynamics based on the ab initio molecular dynamics carried out "on the fly" in the framework of the configuration interaction method combined with Tully's surface hopping algorithm for nonadiabatic transitions. This approach combined with our Wigner distribution approach allows us to perform accurate simulations of femtosecond pump-probe spectra in the systems where radiationless transitions among electronic states take place. In this paper we illustrate this by theoretical simulation of ultrafast processes and nonradiative relaxation in the Na(3)F cluster, involving three excited states and the ground electronic state. Furthermore, we show that our accurate simulation of the photoionization pump-probe spectrum is in full agreement with the experimental signal. Based on the nonadiabatic dynamics at high level of accuracy and taking into account all degrees of freedom, the nonradiative lifetime for the 1 (1)B(1) excited state of Na(3)F has been determined to be approximately 900 fs.     

Nonadiabatic surface hopping Herman-Kluk semiclassical initial value representation method revisited: applications to Tully's three model systems

The nonadiabatic surface hopping Herman-Kluk (HK) semiclassical initial value representation (SC-IVR) method for nonadiabatic problems is reformulated. The method has the same spirit as Tully's surface hopping technique [J. Chem. Phys. 93, 1061 (1990)] and almost keeps the same structure as the original single-surface HK SC-IVR method except that trajectories can hop to other surfaces according to the hopping probabilities and phases, which can be easily integrated along the paths. The method is based on a rather general nonadiabatic semiclassical surface hopping theory developed by Herman [J. Chem. Phys. 103, 8081 (1995)], which has been shown to be accurate to the first order in h and through all the orders of the nonadiabatic coupling amplitude. Our simulation studies on the three model systems suggested by Tully demonstrate that this method is practical and capable of describing nonadiabatic quantum dynamics for various coupling situations in very good agreement with benchmark calculations.