Time‐dependent schrödinger equation with Markovian outgoing wave boundary conditions: Applications to quantum tunneling dynamics and photoionization

Markovian outgoing wave boundary conditions are introduced as an approximate method to reduce the size of the computational grid for time integration of the time‐dependent Schrödinger equation. The ratio and polynomial methods developed as open boundary conditions are applied to the wave function at the boundaries of the computational grid. This computational method is used to study the wave packet dynamics for a metastable well, a double well, and strong‐field ionization of a model atom. Accurate results demonstrate that this method can significantly reduce the number of grid points required in a dynamical calculation for quantum dynamical problems. © 2012 Wiley Periodicals, Inc.     

Quantum dynamics of dissociative adsorption of H2 and D2: The time-dependent wave packet method

A time-dependent quantum wave packet method was used to study the dynamics of dissociative adsorption of H₂ and D₂ on a flat and static surface. The molecule-surface interaction is described using a modified London-Eyring-Polanyi-Sato (LEPS) type potential for the H₂/Ni(100) system. The three-dimensional (3-D) dissociation probabilities were calculated for different initial rovibrational states as a function of initial incident energies. Our results show that the dissociation of the diatomic rotational states whose quantum numbers satisfyj+m = odd is forbidden at low energies for the homonuclear H_z and D₂ due to the selection rule. The effect of the rotational orientation of diatoms on adsorption predicts that the in-plane rotation (m = j) is more favorable for dissociation than the out-of-plane rotation (m = 0). Enhanced dissociation for vibrationally excited molecules and the significant enhancement of the dissociation probability of H2 when compared to D₂ were explained reasonably in terms of quantum mechanical zero-point energies, the tunneling effect and the reflection from an activation barrier. Project supported by the National Natural Science Foundation of China (Grant No. 19694033) and partially by the Science Foundation for Overseas Chinese Scholars and Students, administered by the State Education Commission of China (Grant No. 1992), by the State Key Laboratory of Theoretical and Computational Chemistry of Jilin University at Changchun (Grant No. 98011, and by the Natural Science Foundation of Shandong Province (Grant No. Y96B03022)     

Quantum dynamics of dissociative adsorption of H_2 and D_2:The time-dependent wave packet method

Ｒｅｃｅｎｔｐｒｏｇｒｅｓｓｉｎｓｕｒｆａｃｅｓｐｅｃｔｒｏｓｃｏｐｙａｎｄｍｏｌｅｃｕｌａｒｂｅａｍｓｃａｔｔｅｒｉｎｇａｎｄｄｅｔｅｃｔｉｏｎｔｅｃｈｎｉｑｕｅｓｍａｋｅｓｉｔｎｏｗｐｏｓｓｉｂｌｅａｔｔｈｅｍｉｃｒｏｓｃｏｐｉｃｌｅｖｅｌｔｏａｃｃｕｒａｔｅｌｙｍｅａｓｕｒｅｓｕｃｈｔｈｉｎｇｓａｓｄｉｓｓｏｃｉａｔｉｏｎｒａｔｅｓ,ａｄｓｏｒｂａｔｅｂｉｎｄｉｎｇａｎｄｇｅｏｍｅｔｒｙ,ａｎｄｍｏｂｉｌｉｔｉｅｓｏｆａｄｓｏｒｂｅｄｓｐｅｃｉｅｓｏｎｓｕｒｆａｃｅｓ.Ｔｈｉｓｋｉｎｄｏｆｅ…     

Quantum dynamics of inelastic scattering with a moving grid

The time‐dependent discrete variable representation (TDDVR) of a wave function with grid points defined by the Hermite part of the Gauss–Hermite (G‐H) basis set introduces quantum corrections to classical mechanics. The grid points in this method follow classical trajectory and the approach converges to the exact quantum formulation with sufficient trajectories (TDDVR points) but just with a single grid point; only classical mechanics performs the dynamics. This newly formulated approach (developed for handling time‐dependent molecular quantum dynamics) has been explored to calculate vibrational transitions in the inelastic scattering processes. Traditional quantum mechanical results exhibit an excellent agreement with TDDVR profiles during the entire propagation when enough grid points are included in the quantum‐classical dynamics. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Time-dependent Fourier grid Hamiltonian method for modelling real-time quantum dynamics: Theoretical models and applications

A local grid method for modelling real-time quantum dynamical events is formulated. The formulation is straightforward for 1-D systems. For more than one dimension, appeal has to be made to mean-field approximation of the appropriate kind. The simplest mean-field model results in time-dependent Hartree-Fourier grid method. The relationship of the proposed method with some other methods available in the literature is examined. A few examples of numerical applications dealing with (i) the dynamics of dissociation and ionization processes in diatoms and atoms respectively and (ii) tunnelling dynamics in the intramolecular H-atom transfer phenomenon are presented.     

An extension of the grid empowered molecular simulator to quantum reactive scattering

Within the activities of the D37 COST Action, we have further developed the quantum dynamics framework of the grid empowered molecular simulator (GEMS) implemented on the segment of the European grid available to the COMPCHEM (computational chemistry) virtual organization. GEMS does now include in a full ab initio approach, the evaluation of the detailed quantum (both time dependent and time independent) dynamics of small systems starting from the calculation of the electronic structure properties as well as the direct calculation of thermalized properties. Illustrative, full dimensional applications of the extended simulator to the H + H(2) , N + N(2) , and O + O(2) systems are presented.     

Time-dependent quantum study of the kinetics of the H(2S) + FO(2II) OH(2II) + F(2P) reaction

Quantum mechanical wave packet calculations are carried out for the H((2)S) + FO((2)II) --> OH((2)II) + F((2)P) reaction on the adiabatic potential energy surface of the ground (3)A' triplet state. The state-to-state and state-to-all reaction probabilities for total angular momentum J = 0 have been calculated. The probabilities for J > 0 have been estimated from the J = 0 results by using J-shifting approximation based on a capture model. Then, the integral cross sections and initial state-selected rate constants have been calculated. The calculations show that the initial state-selected reaction probabilities are dominated by many sharp peaks. The reaction cross section does not manifest any sharp oscillations and the initial state-selected rate constants are sensitive to the temperature.     

Code interoperability and standard data formats in quantum chemistry and quantum dynamics: The Q5/D5Cost data model

Code interoperability and the search for domain‐specific standard data formats represent critical issues in many areas of computational science. The advent of novel computing infrastructures such as computational grids and clouds make these issues even more urgent. The design and implementation of a common data format for quantum chemistry (QC) and quantum dynamics (QD) computer programs is discussed with reference to the research performed in the course of two Collaboration in Science and Technology Actions. The specific data models adopted, Q5Cost and D5Cost, are shown to work for a number of interoperating codes, regardless of the type and amount of information (small or large datasets) to be exchanged. The codes are either interfaced directly, or transfer data by means of wrappers; both types of data exchange are supported by the Q5/D5Cost library. Further, the exchange of data between QC and QD codes is addressed. As a proof of concept, the H + H₂ reaction is discussed. The proposed scheme is shown to provide an excellent basis for cooperative code development, even across domain boundaries. Moreover, the scheme presented is found to be useful also as a production tool in the grid distributed computing environment. © 2013 Wiley Periodicals, Inc.     

Efficient evaluation of the accuracy of molecular quantum dynamics on an approximate analytical or interpolated ab initio potential energy surface

Ab initio electronic structure methods have reached a satisfactory accuracy for the calculation of static properties, but remain too expensive for quantum dynamical calculations. Recently, an efficient semiclassical method was proposed to evaluate the accuracy of quantum dynamics on an approximate potential without having to perform the expensive quantum dynamics on the accurate potential. Here, this method is applied for the first time to evaluate the accuracy of quantum dynamics on an approximate analytical or interpolated potential in comparison to the quantum dynamics on an accurate potential obtained by an ab initio electronic structure method. Specifically, the vibrational dynamics of H₂ on a Morse potential is compared with that on the full CI potential, and the photodissociation dynamics of CO₂ on a LEPS potential with that on the excited ¹Π surface computed at the EOM‐CCSD/aug‐cc‐pVDZ level of theory. Finally, the effect of discretization of a potential energy surface on the quantum dynamics is evaluated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2426–2435, 2010     

Stochastic dynamics simulation of alanine dipeptide: Including solvation interaction determined by boundary element method

A merger of the Poisson–Boltzmann equation and stochastic dynamics simulation is examined using illustrative calculations of alanine dipeptide. The boundary element method (BEM) is used to calculate the hydration forces acting on the solute molecule based on the surroundings. Computational efficiency is achieved by the use of a simple hydration model and a coarse boundary element. Nonetheless, the conformational distribution obtained from this new method is reasonable compared with other theoretical and computational results. Detailed analysis has been accomplished in terms of the hydration interactions and solvation energies. The results indicate that the new simulation method provides an obvious improvement over the conventional stochastic dynamics simulation technique. The further improvement of the hydration model and future application to large molecules are also discussed. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1440–1449, 1997     

Two‐dimensional reactive scattering with transmitted quantum trajectories

A simple and easy‐to‐implement method is presented for the study of time‐dependent reaction dynamics by propagating an ensemble of transmitted quantum trajectories. During the trajectory evolution, reflected trajectories are gradually removed and all the remaining trajectories represent the transmitted subensemble. The removal process of reflected trajectories avoids numerical instabilities arising from node formation in the reactant region, and allows stable long‐time propagation of transmitted trajectories. This method is applied to a two‐dimensional model chemical reaction. Excellent computational results are obtained for the time‐dependent reaction probabilities evaluated by the time integration of the probability flux. © 2014 Wiley Periodicals, Inc.     

Differential diffusion quantum Monte Carlo method: determination of potential energy surfaces of molecules

A differential approach for self-optimizing diffusion Monte Carlo calculation was proposed in this paper, which is a new algorithm combining three techniques such as optimizing, diffusion and correlation sampling. This method can be used to directly compute the energy differential between two systems in the diffusion process, making the statistical error of calculation be reduced to order of 10-5 hartree, and recover about more than 80% of the correlation energy. We employed this approach to set up a potential energy surface of a molecule, used a "rigid move" model, and utilized Jacobi transformation to make energy calculation for two configurations of a molecule having good positive correlation. So, an accurate energy differential could be obtained, and the potential energy surface with good quality can be depicted. In calculation, a technique called "post-equilibrium remaining sample was set up firstly, which can save about 50% of computation expense. This novel algorithm was used to study the potenti     

The quantum dynamics of the reactions N+H2(HD,D2) and their vibrational excitation effect

The N(⁴S)+H₂ reaction and its isotopic variants have been investigated by means of time‐dependent quantum wave packet with split operator method on the ground state potential energy surface (Zhai and Han, J. Chem. Phys. 2011, 135, 104314). The reaction probabilities, integral cross sections, branching ratio of the integral cross sections, and effect of vibrational excitation of H₂, HD, and D₂ diatomic molecules are presented and discussed. The results reveal that the intramolecular isotopic effect is greater than the intermolecular one, and that the vibrational excitation of the diatomic molecules can promote the progress of this reaction. In addition, a limited number of rigorous Coriolis coupling calculations of the integral cross sections of the N(⁴S)+H₂ reaction have been carried out. Also shown is that since the Coriolis coupling plays a small part in this accurate quantum calculation, the cheaper centrifugal sudden calculations here reported are effective for this reactive system. © 2014 Wiley Periodicals, Inc.     

Time‐dependent Fourier grid Hamiltonian method for studying the curve crossing dynamics: Dynamics of predissociation of NaI

A time‐dependent Fourier route is proposed to study quantum dynamics in the presence of curve crossing. The method is tested on a model problem and is then used to study the predissociation dynamics of NaI by a femtosecond laser pulse. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Applications of automatic differentiation to the time‐dependent Schrödinger equation

A new approach based upon the Taylor series method is proposed for propagating solutions of the time‐dependent Schrödinger equation. Replacing the spatial derivative of the wave function with finite difference formulas, we derive a recursive formula for the evaluation of Taylor coefficients. The automatic differentiation technique is used to recursively calculate the required Taylor coefficients. We also develop an implicit scheme for the recursive evaluation of these coefficients. We then advance the solution in time using a Taylor series expansion. Excellent computational results are obtained when this method is applied to a one‐dimensional reflectionless potential and a two‐dimensional barrier transmission problem. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Quantum dynamics of relaxation of a pair of coupled Morse oscillators: Effects of mass and electrical asymmetries

A multiconfiguration time‐dependent Hartree method based recipe is used to study the role of mass and electrical asymmetry in controlling the quantum dynamics of relaxation of a locally excited O? H bond in a water molecule, modeled by a pair of interacting Morse oscillators. The fast periodic energy transfer between the two equivalent O? H bonds in H? O? H is replaced by a rather slow process when one of the H atom is replaced by a deuterium atom. Application of static electric field along the O? D bond in HOD molecule is seen to either enhance or damp the relaxation rate, depending on the strength of the applied field. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Rapid QM/MM approach for biomolecular systems under periodic boundary conditions: Combination of the density‐functional tight‐binding theory and particle mesh Ewald method

A quantum mechanical/molecular mechanical (QM/MM) approach based on the density‐functional tight‐binding (DFTB) theory is a useful tool for analyzing chemical reaction systems in detail. In this study, an efficient QM/MM method is developed by the combination of the DFTB/MM and particle mesh Ewald (PME) methods. Because the Fock matrix, which is required in the DFTB calculation, is analytically obtained by the PME method, the Coulomb energy is accurately and rapidly computed. For assessing the performance of this method, DFTB/MM calculations and molecular dynamics simulation are conducted for a system consisting of two amyloid‐β(1‐16) peptides and a zinc ion in explicit water under periodic boundary conditions. As compared with that of the conventional Ewald summation method, the computational cost of the Coulomb energy by utilizing the present approach is drastically reduced, i.e., 166.5 times faster. Furthermore, the deviation of the electronic energy is less than . © 2016 Wiley Periodicals, Inc.     

Three pillars for achieving quantum mechanical molecular dynamics simulations of huge systems: Divide‐and‐conquer,density‐functional tight‐binding,and massively parallel computation

The linear‐scaling divide‐and‐conquer (DC) quantum chemical methodology is applied to the density‐functional tight‐binding (DFTB) theory to develop a massively parallel program that achieves on‐the‐fly molecular reaction dynamics simulations of huge systems from scratch. The functions to perform large scale geometry optimization and molecular dynamics with DC‐DFTB potential energy surface are implemented to the program called DC‐DFTB‐K. A novel interpolation‐based algorithm is developed for parallelizing the determination of the Fermi level in the DC method. The performance of the DC‐DFTB‐K program is assessed using a laboratory computer and the K computer. Numerical tests show the high efficiency of the DC‐DFTB‐K program, a single‐point energy gradient calculation of a one‐million‐atom system is completed within 60 s using 7290 nodes of the K computer. © 2016 Wiley Periodicals, Inc.