Exact quantum master equation via the calculus on path integrals

An exact quantum master equation formalism is constructed for the efficient evaluation of quantum non-Markovian dissipation beyond the weak system-bath interaction regime in the presence of time-dependent external field. A novel truncation scheme is further proposed and compared with other approaches to close the resulting hierarchically coupled equations of motion. The interplay between system-bath interaction strength, non-Markovian property, and required level of hierarchy is also demonstrated with the aid of simple spin-boson systems.     

Non-Born-Oppenheimer path in anti-Hermitian dynamics for nonadiabatic transitions

A serious difficulty in the semiclassical Ehrenfest theory for nonadiabatic transitions is that a path passing across the avoided crossing is forced to run on a potential averaged over comprising adiabatic potential surfaces that commit the avoided crossing. Therefore once a path passes through the crossing region, it immediately becomes incompatible with the standard view of "classical trajectory" running on an adiabatic surface. This casts a fundamental question to the theoretical structure of chemical dynamics. In this paper, we propose a non-Born-Oppenheimer path that is generated by an anti-Hermitian Hamiltonian, whose complex-valued eigenenergies can cross in their real parts and avoid crossing in the imaginary parts in the vicinity of the nonadiabatic transition region. We discuss the properties of this non-Born-Oppenheimer path and thereby show its compatibility with the Born-Oppenheimer classical trajectories. This theory not only allows the geometrical branching of the paths but gives the nonadiabatic transition amplitudes and quantum phases along the generated paths.     

A surface hopping algorithm for nonadiabatic minimum energy path calculations

The article introduces a robust algorithm for the computation of minimum energy paths transiting along regions of near‐to or degeneracy of adiabatic states. The method facilitates studies of excited state reactivity involving weakly avoided crossings and conical intersections. Based on the analysis of the change in the multiconfigurational wave function the algorithm takes the decision whether the optimization should continue following the same electronic state or switch to a different state. This algorithm helps to overcome convergence difficulties near degeneracies. The implementation in the MOLCAS quantum chemistry package is discussed. To demonstrate the utility of the proposed procedure four examples of application are provided: thymine, asulam, 1,2‐dioxetane, and a three‐double‐bond model of the 11‐cis‐retinal protonated Schiff base. © 2015 Wiley Periodicals, Inc.     

Phase-equilibrium calculations for n-alkane + alkanol systems using continuous thermodynamics

A new association model based on continuous thermodynamics is introduced and applied to six systems of the type n-alkane (n-hexane, n-heptane, n-octane) + alkanol (methanol, ethanol). The alkanol is considered to be a mixture of chain associates with the composition described by a continuous distribution function. This distribution function is derived as an analytical expression from the mass action law applied to the association equilibrium. To consider the entropic contribution originating from the size differences of the molecules (associates) activity coefficients based on Flory–Huggins model are included in the mass action law. Unlike the molecular-mass distribution of a polymer the chain-length distribution of the associates depends on the temperature and on the mole fraction of the alkanol. The treatment of vapor–liquid equilibrium and liquid–liquid equilibrium is similar to that of an oil system or of a polymer solution using continuous thermodynamics. Different to other chemical models of association there is no additive split into a physical and a chemical contribution. The equilibrium constants of association were fitted to vapor-pressure data of methanol and ethanol. The model needs only one interaction parameter being independent of temperature and taking the same value for all systems studied. Considering the simplicity of the model, both the liquid–liquid equilibrium of the three methanol systems and the vapor–liquid equilibrium of all six systems are predicted with reasonable accuracy.     

On the properties of a primitive semiclassical surface hopping propagator for nonadiabatic quantum dynamics

A previously developed nonadiabatic semiclassical surface hopping propagator [M. F. Herman J. Chem. Phys. 103, 8081 (1995)] is further studied. The propagator has been shown to satisfy the time-dependent Schrodinger equation (TDSE) through order h, and the O(h2) terms are treated as small errors, consistent with standard semiclassical analysis. Energy is conserved at each hopping point and the change in momentum accompanying each hop is parallel to the direction of the nonadiabatic coupling vector resulting in both transmission and reflection types of hops. Quantum mechanical analysis and numerical calculations presented in this paper show that the h2 terms involving the interstate coupling functions have significant effects on the quantum transition probabilities. Motivated by these data, the h2 terms are analyzed for the nonadiabatic semiclassical propagator. It is shown that the propagator can satisfy the TDSE for multidimensional systems by including another type of nonclassical trajectories that reflect on the same surfaces. This h2 analysis gives three conditions for these three types of trajectories so that their coefficients are uniquely determined. Besides the nonadiabatic semiclassical propagator, a numerically useful quantum propagator in the adiabatic representation is developed to describe nonadiabatic transitions.     

Exact on-event expressions for discrete potential systems

The properties of systems composed of atoms interacting though discrete potentials are dictated by a series of events which occur between pairs of atoms. There are only four basic event types for pairwise discrete potentials and the square-well/shoulder systems studied here exhibit them all. Closed analytical expressions are derived for the on-event kinetic energy distribution functions for an atom, which are distinct from the Maxwell-Boltzmann distribution function. Exact expressions are derived that directly relate the pressure and temperature of equilibrium discrete potential systems to the rates of each type of event. The pressure can be determined from knowledge of only the rate of core and bounce events. The temperature is given by the ratio of the number of bounce events to the number of disassociation/association events. All these expressions are validated with event-driven molecular dynamics simulations and agree with the data within the statistical precision of the simulations.     

A quantum generalization of intrinsic reaction coordinate using path integral centroid coordinates

We propose a generalization of the intrinsic reaction coordinate (IRC) for quantum many-body systems described in terms of the mass-weighted ring polymer centroids in the imaginary-time path integral theory. This novel kind of reaction coordinate, which may be called the "centroid IRC," corresponds to the minimum free energy path connecting reactant and product states with a least amount of reversible work applied to the center of masses of the quantum nuclei, i.e., the centroids. We provide a numerical procedure to obtain the centroid IRC based on first principles by combining ab initio path integral simulation with the string method. This approach is applied to NH(3) molecule and N(2)H(5) (-) ion as well as their deuterated isotopomers to study the importance of nuclear quantum effects in the intramolecular and intermolecular proton transfer reactions. We find that, in the intramolecular proton transfer (inversion) of NH(3), the free energy barrier for the centroid variables decreases with an amount of about 20% compared to the classical one at the room temperature. In the intermolecular proton transfer of N(2)H(5) (-), the centroid IRC is largely deviated from the "classical" IRC, and the free energy barrier is reduced by the quantum effects even more drastically.     

Open path atmospheric spectroscopy using room temperature operated pulsed quantum cascade laser

We report the application of a distributed feedback quantum cascade laser for 5.8 km long open path spectroscopic monitoring of ozone, water vapor and CO(2). The thermal chirp during a 140 or 200 ns long excitation pulse is used for fast wavelength scanning. The fast wavelength scanning has the advantage of the measured spectra not being affected by atmospheric turbulence, which is essential for long open path measurements. An almost linear tuning of about 0.6 and 1.2 cm(-1) is achieved, respectively. Lines from the nu(3) vibrational band of the ozone spectra centered at 1,031 and 1,049 cm(-1) is used for ozone detection by differential absorption. The lowest column densities (LCD) for ozone of the order of 0.3 ppmm retrieved from the absorption spectra for averaging times less than 20s are better then the LCD value of 2 ppmm measured with UV DOAS systems. The intrinsic haze immunity of mid-IR laser sources is an additional important advantage of mid-IR open path spectroscopy, compared with standard UV-vis DOAS. The third major advantage of the method is the possibility to measure more inorganic and organic atmospheric species compared to the UV-vis DOAS.     

Exact dynamics of dissipative electronic systems and quantum transport: Hierarchical equations of motion approach

A generalized quantum master equation theory that governs the exact, nonperturbative quantum dissipation and quantum transport is formulated in terms of hierarchically coupled equations of motion for an arbitrary electronic system in contact with electrodes under either a stationary or a nonstationary electrochemical potential bias. The theoretical construction starts with the influence functional in path integral, in which the electron creation and annihilation operators are Grassmann variables. Time derivatives on the influence functionals are then performed in a hierarchical manner. Both the multiple-frequency dispersion and the non-Markovian reservoir parametrization schemes are considered for the desired hierarchy construction. The resulting hierarchical equations of motion formalism is in principle exact and applicable to arbitrary electronic systems, including Coulomb interactions, under the influence of arbitrary time-dependent applied bias voltage and external fields. Both the conventional quantum master equation and the real-time diagrammatic formalism of Schon and co-workers can be readily obtained at well defined limits of the present theory. We also show that for a noninteracting electron system, the present hierarchical equations of motion formalism terminates at the second tier exactly, and the Landuer-Buttiker transport current expression is recovered. The present theory renders an exact and numerically tractable tool to evaluate various transient and stationary quantum transport properties of many-electron systems, together with the involving nonperturbative dissipative dynamics.     

Reaction path potential for complex systems derived from combined ab initio quantum mechanical and molecular mechanical calculations

Combined ab initio quantum mechanical and molecular mechanical calculations have been widely used for modeling chemical reactions in complex systems such as enzymes, with most applications being based on the determination of a minimum energy path connecting the reactant through the transition state to the product in the enzyme environment. However, statistical mechanics sampling and reaction dynamics calculations with a combined ab initio quantum mechanical (QM) and molecular mechanical (MM) potential are still not feasible because of the computational costs associated mainly with the ab initio quantum mechanical calculations for the QM subsystem. To address this issue, a reaction path potential energy surface is developed here for statistical mechanics and dynamics simulation of chemical reactions in enzymes and other complex systems. The reaction path potential follows the ideas from the reaction path Hamiltonian of Miller, Handy and Adams for gas phase chemical reactions but is designed specifically for large systems that are described with combined ab initio quantum mechanical and molecular mechanical methods. The reaction path potential is an analytical energy expression of the combined quantum mechanical and molecular mechanical potential energy along the minimum energy path. An expansion around the minimum energy path is made in both the nuclear and the electronic degrees of freedom for the QM subsystem internal energy, while the energy of the subsystem described with MM remains unchanged from that in the combined quantum mechanical and molecular mechanical expression and the electrostatic interaction between the QM and MM subsystems is described as the interaction of the MM charges with the QM charges. The QM charges are polarizable in response to the changes in both the MM and the QM degrees of freedom through a new response kernel developed in the present work. The input data for constructing the reaction path potential are energies, vibrational frequencies, and electron density response properties of the QM subsystem along the minimum energy path, all of which can be obtained from the combined quantum mechanical and molecular mechanical calculations. Once constructed, it costs much less for its evaluation. Thus, the reaction path potential provides a potential energy surface for rigorous statistical mechanics and reaction dynamics calculations of complex systems. As an example, the method is applied to the statistical mechanical calculations for the potential of mean force of the chemical reaction in triosephosphate isomerase.     

General laser interaction theory in atom-diatom systems for both adiabatic and nonadiabatic cases

This paper develops the general theory for laser fields interacting with bimolecular systems. In this study, we choose to use the multipolar gauge on the basis of gauge invariance. We consider both the adiabatic and nonadiabatic cases and find they produce similar interaction pictures. As an application of this theory, we present the study of rovibrational energy transfer in Ar + CO collisions in the presence of an intense laser field.     

Single molecule photon counting statistics for quantum mechanical chromophore dynamics

We extend the generating function technique for calculation of single molecule photon emission statistics (Zheng, Y.; Brown, F. L. H. Phys. Rev. Lett. 2003, 90, 238305) to systems governed by multi-level quantum dynamics. This opens up the possibility to study phenomena that are outside the realm of purely stochastic and mixed quantum-stochastic models. In particular, the present methodology allows for calculation of photon statistics that are spectrally resolved and subject to quantum coherence. Several model calculations illustrate the generality of the technique and highlight quantitative and qualitative differences between quantum mechanical models and related stochastic approximations when they arise. Calculations suggest that studying photon statistics as a function of photon frequency has the potential to reveal more about system dynamics than the usual broadband detection schemes.     

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.     

Exact solution of multidimensional hyper‐radial Schrödinger equation for many‐electron quantum systems

In quantum theory, solving Schrödinger equation analytically for larger atomic and molecular systems with cluster of electrons and nuclei persists to be a tortuous challenge. Here, we consider, Schrödinger equation in arbitrary N‐dimensional space corresponding to inverse‐power law potential function originating from a multitude of interactions participating in a many‐electron quantum system for exact solution within the framework of Frobenius method via the formulation of an ansatz to the hyper‐radial wave function. Analytical expressions for energy spectra, and hyper‐radial wave functions in terms of known coefficients of inverse‐power potential function, and wave function parameters have been obtained. A generalized two‐term recurrence relation for power series expansion coefficients has been established. © 2016 Wiley Periodicals, Inc.     

SWKB formalism for confined quantum systems

A formalism is developed to study spatially confined one‐dimensional quantum mechanical systems in the framework of the supersymmetric Wentzel–Kramers–Brillouin (SWKB) method. The approximation technique is applied to estimate the energy eigenvalues of two confined potentials—the harmonic oscillator V(x)=x² and the screened Coulomb potential V(x)=−V₀sech²x. The results thus obtained are found to be in better agreement with the exact numerical values than are those from the ordinary WKB approach. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 79: 267–273, 2000     

On identifiability for chemical systems from measurable variables

The dynamics of the composition of chemical species in reacting systems can be characterized by a set of autonomous differential equations derived from mass conservation principles and some elementary hypothesis related to chemical reactivity. These sets of ordinary differential equations (ODEs) are basically non-linear, their complexity grows as much increases the number of substances present in the reacting media and can be characterized by a set of phenomenological constants (kinetic rate constants) which contains all the relevant information about the physical system. The determination of these kinetic constants is critical for the design or control of chemical systems from a technological point of view but the non-linear nature of the ODEs implies that there are hidden correlations between the parameters which maybe can be revealed with a identifiability analysis.     

Feynman path integral representation for many-fermion interacting systems

A many-fermion interacting system is investigated within the scenario of the Feynman path integral representation of quantum mechanics. Short-time propagator algorithms and a basis set, closely related to the coherent states, are used to obtain the many-body analytic propagator. A second-quantized Hamiltonian involving a restricted set of two-body interactions and the whole set of Coulomb interactions are separately and shown to lead to an exact and an approximate propagator, respectively. In the latter case, use of a grand canonical ensemble allows the grand partition function and the density operator matrix to be readily obtained. No further approximations are required in the calculation of the trace of the evolution operator involved in the evaluation of statistical expectation values. © 1994 John Wiley & Sons, Inc.