On the Bayesian approach to calculating time correlation functions in quantum systems; reaction dynamics and spectroscopy

An open-chain imaginary-time path-integral sampling approach to the calculation of approximate symmetrized quantum time correlation functions

We introduce a scheme for approximating quantum time correlation functions numerically within the Feynman path integral formulation. Starting with the symmetrized version of the correlation function expressed as a discretized path integral, we introduce a change of integration variables often used in the derivation of trajectory-based semiclassical methods. In particular, we transform to sum and difference variables between forward and backward complex-time propagation paths. Once the transformation is performed, the potential energy is expanded in powers of the difference variables, which allows us to perform the integrals over these variables analytically. The manner in which this procedure is carried out results in an open-chain path integral (in the remaining sum variables) with a modified potential that is evaluated using imaginary-time path-integral sampling rather than requiring the generation of a large ensemble of trajectories. Consequently, any number of path integral sampling schemes can be employed to compute the remaining path integral, including Monte Carlo, path-integral molecular dynamics, or enhanced path-integral molecular dynamics. We believe that this approach constitutes a different perspective in semiclassical-type approximations to quantum time correlation functions. Importantly, we argue that our approximation can be systematically improved within a cumulant expansion formalism. We test this approximation on a set of one-dimensional problems that are commonly used to benchmark approximate quantum dynamical schemes. We show that the method is at least as accurate as the popular ring-polymer molecular dynamics technique and linearized semiclassical initial value representation for correlation functions of linear operators in most of these examples and improves the accuracy of correlation functions of nonlinear operators.

Conformational Transition Pathways of Epidermal Growth Factor Receptor Kinase Domain from Multiple Molecular Dynamics Simulations and Bayesian Clustering

The epidermal growth factor receptor (EGFR) is aberrantly activated in various cancer cells and an important target for cancer treatment. Deep understanding of EGFR conformational changes between the active and inactive states is of pharmaceutical interest. Here we present a strategy combining multiply targeted molecular dynamics simulations, unbiased molecular dynamics simulations, and Bayesian clustering to investigate transition pathways during the activation/inactivation process of EGFR kinase domain. Two distinct pathways between the active and inactive forms are designed, explored, and compared. Based on Bayesian clustering and rough two-dimensional free energy surfaces, the energy-favorable pathway is recognized, though DFG-flip happens in both pathways. In addition, another pathway with different intermediate states appears in our simulations. Comparison of distinct pathways also indicates that disruption of the Lys745-Glu762 interaction is critically important in DFG-flip while movement of the A-loop significantly facilitates the conformational change. Our simulations yield new insights into EGFR conformational transitions. Moreover, our results verify that this approach is valid and efficient in sampling of protein conformational changes and comparison of distinct pathways.

A Bayesian statistics approach to multiscale coarse graining

Coarse-grained (CG) modeling provides a promising way to investigate many important physical and biological phenomena over large spatial and temporal scales. The multiscale coarse-graining (MS-CG) method has been proven to be a thermodynamically consistent way to systematically derive a CG model from atomistic force information, as shown in a variety of systems, ranging from simple liquids to proteins embedded in lipid bilayers. In the present work, Bayes' theorem, an advanced statistical tool widely used in signal processing and pattern recognition, is adopted to further improve the MS-CG force field obtained from the CG modeling. This approach can regularize the linear equation resulting from the underlying force-matching methodology, therefore substantially improving the quality of the MS-CG force field, especially for the regions with limited sampling. Moreover, this Bayesian approach can naturally provide an error estimation for each force field parameter, from which one can know the extent the results can be trusted. The robustness and accuracy of the Bayesian MS-CG algorithm is demonstrated for three different systems, including simple liquid methanol, polyalanine peptide solvated in explicit water, and a much more complicated peptide assembly with 32 NNQQNY hexapeptides.

On the short-time limit of ring polymer molecular dynamics

We examine the short-time accuracy of a class of approximate quantum dynamical techniques that includes the centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD) methods. Both of these methods are based on the path integral molecular dynamics (PIMD) technique for calculating the exact static equilibrium properties of quantum mechanical systems. For Kubo-transformed real-time correlation functions involving operators that are linear functions of positions or momenta, the RPMD and (adiabatic) CMD approximations differ only in the choice of the artificial mass matrix of the system of ring polymer beads that is employed in PIMD. The obvious ansatz for a general method of this type is therefore to regard the elements of the PIMD (or Parrinello-Rahman) mass matrix as an adjustable set of parameters that can be chosen to improve the accuracy of the resulting approximation. We show here that this ansatz leads uniquely to the RPMD approximation when the criterion that is used to select the mass matrix is the short-time accuracy of the Kubo-transformed correlation function. In particular, we show that the leading error in the RPMD position autocorrelation function is O(t(8)) and the error in the velocity autocorrelation function is O(t(6)), for a general anharmonic potential. The corresponding errors in the CMD approximation are O(t(6)) and O(t(4)), respectively.

Evaluation of quantum correlation functions from classical data: Anharmonic models

The previously introduced method of evaluating quantum mechanical time correlation functions using as input only classical simulation data is generalized and applied to two anharmonic model systems, as a further test. The quantum correction approach utilizes the relation between a general quantum correlation function and its classical analog. For the tested models, we obtain numerical results of nonlinear correlation functions with comparable accuracy to that of the centroid molecular dynamics method, although the present method is much simpler to implement and not limited to real valued quantum correlation functions.

Hybrid quantum/classical molecular dynamics for a proton transfer reaction coupled to a dissipative bath

A hybrid quantum/classical molecular dynamics approach is applied to a proton transfer reaction represented by a symmetric double well system coupled to a dissipative bath. In this approach, the proton is treated quantum mechanically and all bath modes are treated classically. The transition state theory rate constant is obtained from the potential of mean force, which is generated along a collective reaction coordinate with umbrella sampling techniques. The transmission coefficient, which accounts for dynamical recrossings of the dividing surface, is calculated with a reactive flux approach combined with the molecular dynamics with quantum transitions surface hopping method. The hybrid quantum/classical results agree well with numerically exact results in the spatial-diffusion-controlled regime, which is most relevant for proton transfer in proteins. This hybrid quantum/classical approach has already been shown to be computationally practical for studying proton transfer in large biological systems. These results have important implications for future applications to hydrogen transfer reactions in solution and proteins.

Time-dependent perturbation theory for vibrational energy relaxation and dephasing in peptides and proteins

Without invoking the Markov approximation, we derive formulas for vibrational energy relaxation (VER) and dephasing for an anharmonic system oscillator using a time-dependent perturbation theory. The system-bath Hamiltonian contains more than the third order coupling terms since we take a normal mode picture as a zeroth order approximation. When we invoke the Markov approximation, our theory reduces to the Maradudin-Fein formula which is used to describe the VER properties of glass and proteins. When the system anharmonicity and the renormalization effect due to the environment vanishes, our formulas reduce to those derived by and Mikami and Okazaki [J. Chem. Phys. 121, 10052 (2004)] invoking the path-integral influence functional method with the second order cumulant expansion. We apply our formulas to VER of the amide I mode of a small amino-acid like molecule, N-methylacetamide, in heavy water.

Semiclassical calculation of nonadiabatic thermal rate constants: Application to condensed phase reactions

The nonadiabatic transition state theory proposed recently by Zhao et al. [J. Chem. Phys. 121, 8854 (2004)] is extended to calculate rate constants of complex systems by using the Monte Carlo and umbrella sampling methods. Surface hopping molecular dynamics technique is incorporated to take into account the dynamic recrossing effect. A nontrivial benchmark model of the nonadiabatic reaction in the condensed phase is used for the numerical test. It is found that our semiclassical results agree well with those produced by the rigorous quantum mechanical method. Comparing with available analytical approaches, we find that the simple statistical theory proposed by Straub and Berne [J. Chem. Phys. 87, 6111 (1987)] is applicable for a wide friction region although their formula is obtained using Landau-Zener [Phys. Z. Sowjetunion 2, 46 (1932); Proc. R. Soc. London, Ser. A 137, 696 (1932)] nonadiabatic transition probability along a one-dimensional diffusive coordinate. We also investigate how the nuclear tunneling events affect the dependence of the rate constant on the friction.

Transport properties of normal liquid helium: Comparison of various methodologies

We revisit the problem of self-diffusion in normal liquid helium above the lambda transition. Several different methods are applied to compute the velocity autocorrelation function. Since it is still impossible to determine the exact result for the velocity autocorrelation function from simulation, we appeal to the computation of short-time moments to determine the accuracy of the different approaches at short times. The main conclusion reached from our study is that both the quantum mode-coupling theory and the numerical analytic continuation approach must be regarded as a viable and competitive methods for the computation of dynamical properties of quantum systems.

QUANTUM MODE-COUPLING THEORY: Formulation and Applications to Normal and Supercooled Quantum Liquids

▪ Abstract  We review our recent efforts to formulate and study a mode-coupling approach to real-time dynamic fluctuations in quantum liquids. Comparison is made between the theory and recent neutron scattering experiments performed on liquid ortho-deuterium and para-hydrogen. We discuss extensions of the theory to supercooled and glassy states where quantum fluctuations compete with thermal fluctuations. Experimental scenarios for quantum glassy liquids are briefly discussed.

Quantum mode-coupling theory for binary mixtures

We extend the quantum mode-coupling theory of neat liquids to the case of binary mixtures, in order to study supercooled liquids where quantum fluctuations may compete with thermal fluctuations. We apply the theory to a generic model of a binary mixture of Lennard-Jones particles. Our treatment may be used to study quantum aging and exotic glass melting scenarios in structural supercooled quantum liquids.

A fully self-consistent treatment of collective fluctuations in quantum liquids

The problem of calculating collective density fluctuations in quantum liquids is revisited. A fully quantum mechanical self-consistent treatment based on a quantum mode-coupling theory is presented. The theory is compared with the maximum entropy analytic continuation approach and with available experimental results. The quantum mode-coupling theory provides semiquantitative results for both short and long time dynamics. The proper description of long time phenomena is important in future study of problems related to the physics of glassy quantum systems, and to the study of collective fluctuations in Bose fluids.

The calculation of transport properties in quantum liquids using the maximum entropy numerical analytic continuation method: application to liquid para-hydrogen

We present a method based on augmenting an exact relation between a frequency-dependent diffusion constant and the imaginary time velocity autocorrelation function, combined with the maximum entropy numerical analytic continuation approach to study transport properties in quantum liquids. The method is applied to the case of liquid para-hydrogen at two thermodynamic state points: a liquid near the triple point and a high-temperature liquid. Good agreement for the self-diffusion constant and for the real-time velocity autocorrelation function is obtained in comparison to experimental measurements and other theoretical predictions. Improvement of the methodology and future applications are discussed.