Symmetrized correlation function for liquid para-hydrogen using complex-time pair-product propagators

We present a simple and efficient method for calculating symmetrized time correlation functions of neat quantum fluids. Using the pair-product approximation to each complex-time quantum mechanical propagator, symmetrized correlation functions are written in terms of a double integral for each degree of freedom with a purely positive integrand. At moderate temperatures and densities, where the pair-product approximation to the Boltzmann operator is sufficiently accurate, the method leads to quantitative results for the early time part of the correlation function. The method is tested extensively on liquid para-hydrogen at 25 K and used to obtain accurate quantum mechanical results for the initial 0.2 ps segment of the symmetrized velocity autocorrelation function of this system, as well as the incoherent dynamic structure factor at certain momentum transfer values.     

Coherent exciton transport in dendrimers and continuous-time quantum walks

We model coherent exciton transport in dendrimers by continuous-time quantum walks. For dendrimers up to the second generation the coherent transport shows perfect recurrences when the initial excitation starts at the central node. For larger dendrimers, the recurrence ceases to be perfect, a fact which resembles results for discrete quantum carpets. Moreover, depending on the initial excitation site, we find that the coherent transport to certain nodes of the dendrimer has a very low probability. When the initial excitation starts from the central node, the problem can be mapped onto a line which simplifies the computational effort. Furthermore, the long time average of the quantum mechanical transition probabilities between pairs of nodes shows characteristic patterns and allows us to classify the nodes into clusters with identical limiting probabilities. For the (space) average of the quantum mechanical probability to be still or to be again at the initial site, we obtain, based on the Cauchy-Schwarz inequality, a simple lower bound which depends only on the eigenvalue spectrum of the Hamiltonian.     

Interpreting nonlinear vibrational spectroscopy with the classical mechanical analogs of double-sided Feynman diagrams

Observables in coherent, multiple-pulse infrared spectroscopy may be computed from a vibrational nonlinear response function. This response function is conventionally calculated quantum-mechanically, but the challenges in applying quantum mechanics to large, anharmonic systems motivate the examination of classical mechanical vibrational nonlinear response functions. We present an approximate formulation of the classical mechanical third-order vibrational response function for an anharmonic solute oscillator interacting with a harmonic solvent, which establishes a clear connection between classical and quantum mechanical treatments. This formalism permits the identification of the classical mechanical analog of the pure dephasing of a quantum mechanical degree of freedom, and suggests the construction of classical mechanical analogs of the double-sided Feynman diagrams of quantum mechanics, which are widely applied to nonlinear spectroscopy. Application of a rotating wave approximation permits the analytic extraction of signals obeying particular spatial phase matching conditions from a classical-mechanical response function. Calculations of the third-order response function for an anharmonic oscillator coupled to a harmonic solvent are compared to numerically correct classical mechanical results.     

Time-dependent quantum dynamical simulations of C2 condensation under extreme conditions

We report theoretical studies of the initial phase of bulk C(2) condensation into carbon nano-structures by means of Born-Oppenheimer and time-dependent quantum mechanical Liouville-von Neumann molecular dynamics based on the density-functional tight-binding (DFTB) framework for electrons. We observe that the time-dependent quantum mechanical approach leads to faster formation of carbon nanostructures than analogous Born-Oppenheimer simulations. Our results suggest that the condensation of bulk carbon is nonadiabatic in nature, with the critical role of electronic stopping as in ion-irradiation of materials. Contrary to time-dependent quantum mechanical simulations, Born-Oppenheimer dynamics incorrectly predict that the short carbon chains obtained from initial reactive collisions between C(2) quickly evaporate, leading to much lower probability of secondary collisions and condensation. We also discuss some deficiencies in Born-Oppenheimer dynamics that lead to unphysical charge polarization and electron transfer.     

A semiclassical study of the thermal conductivity of low temperature liquids

The conventional classical energy current auto-correlation function has been extended into a quantum mechanical version and then approximated by the linearized semiclassical initial value representation approach. Comparison of the thermal conductivity to simulation results shows that about 15% quantum correction to the classical molecular dynamics results for liquid neon are quantitatively predicted. For liquid para-hydrogen the quantum effects are sufficiently large that the linearized semiclassical approach is only 20% accurate, while for both liquid He(4) and He(3) the thermal conductivity disagrees by a factor of 2, although exchange effects appear to play a minor role.     

On the correspondence between classical and quantum mechanics in macromolecular systems: Too much classical chaos

In a recent study we found the classical dynamics of a polyethylene (PE) chain to exhibit low dimensional chaos at temperatures as low as a few Kelvin. These results strongly suggest that classical molecular dynamic simulations in polymer systems can grossly overestimate vibrational motion, which consequently results in disordered structures. In contrast, quantum mechanical calculations using Internal Coordinate Quantum Monte Carlo (an improved method with an initial conjecture for the correct wave function) indicate that the quantum ground state for a three-dimensional model PE chain is far more rigid than determined from molecular dynamics (MD) simulations, even at energies as low as a small fraction of the ground state energy. This result casts uncertainty on the reliability of MD estimates of dynamical or structural quantities relevant to the study of some macromolecular systems.     

Optimal choice of dividing surface for the computation of quantum reaction rates

We consider the calculation of quantum mechanical rate constants for chemical reactions via algorithms that utilize short-time values of the symmetrized flux-flux correlation function. We argue that the dividing surface that makes optimal use of the short-time quantum information is the surface that minimizes the value at the origin of the symmetrized flux-flux correlation function. We also demonstrate that, in the classical limit, this quantum variational criterion produces the same dividing surface as Wigner's variational principle. Finally, we argue that the quantum variational criterion behaves in a nearly optimal fashion with respect to the minimization of the extent of re-crossing flux.     

Calculations of vibrational energy transfer using semi-classicals matrix theory

We utilize an extension of Miller's semi-classical S matrix theory to calculate resonant and non-resonant vibrational energy transfer probabilities. The collisions under study are collinear D₂D₂ interactions at energies below and above the classical dynamic transition threshold and collinear H₂D₂ interactions above the dynamic threshold. Below threshold we employ an initial angle representation. We find that the computed probabilities are in substantial agreement with exact quantum mechanical computations and represent a major improvement over quasi-classical results. At energies above threshold we apply the first order and the classical semi-classical versions of the theory. The results indicate fair agreement with quantum mechanical calculations, but no significant improvement over quasi-classical results.     

Classical and quantum mechanical calculations of conformational probability density distributions. Two-rotor molecules

In order to study the dynamical structure of a two-rotor molecule, such as acetone, as a function of temperature, conformational probability density distributions are computed by using three different approaches: the so-called current approach, the classical approach, and the quantum mechanical oscillator approach. It is found that the three procedures yield comparable results, at least at normal temperature (25°C), although the current and, especially, the classical approaches give rise to too sharp distributions when compared with the quantum mechanical results. Owing to its simplicity, the current approach may be used advantageously, and it is easily extendible to many-rotor systems. Finally, it is verified that deuteration does not affect appreciably the conformation.     

Quantum crystallography applied to crystalline maleic anhydride

Quantum crystallography (QCr) is a term that concerns techniques for using crystallographic information to enhance quantum mechanical calculations and the information derived from them. In our approach to QCr, we use molecular orbitals and a single‐determinant density matrix formalism to develop a quantum mechanical model. Our initial application to a test material, crystalline maleic anhydride, involved the adjustment of the elements in the density (projector) matrix and some others in the quantum mechanical model. The purpose was to optimize the fit between the experimental structure factor magnitudes and the values of those magnitudes obtained from the quantum mechanical model. The adjustment of the projector matrix preserved the idempotency and normalization properties of the matrix. In this application, it was also found that it was necessary to correct the X‐ray diffraction data for systematic errors. An effective statistical method for doing this was developed from quantum mechanical theory. There were a number of special features of this investigation that emerged as it progressed. The mirror plane in maleic anhydride, for example, was quite useful because, in the absence of significant interactions between the molecules in the crystal, charge distributions on both sides of the mirror plane should be essentially the same. Deviations raised questions that resulted in improved procedures. The quality of theoretical results as a function of basis set and mode of calculation is also part of this investigation. One result of the information obtained from various aspects of this study is the potential for greater efficiency in the procedures and calculations. The calculations for maleic anhydride based on its structure concern the number of electrons per atom, various energies, and electron density contours. Related theoretical calculations based on geometry optimization were also made. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 439–450, 1999     

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

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

Thermodynamic integration from classical to quantum mechanics

We present a new method for calculating quantum mechanical corrections to classical free energies, based on thermodynamic integration from classical to quantum mechanics. In contrast to previous methods, our method is numerically stable even in the presence of strong quantum delocalization. We first illustrate the method and its relationship to a well-established method with an analysis of a one-dimensional harmonic oscillator. We then show that our method can be used to calculate the quantum mechanical contributions to the free energies of ice and water for a flexible water model, a problem for which the established method is unstable.     

Effects of anharmonicity on nonadiabatic electron transfer: a model

The effect of anharmonicity in the intramolecular modes of a model system for exothermic intramolecular nonadiabatic electron transfer is probed by examining the dependence of the transition probability on the exoergicity. The Franck-Condon factor for the Morse potential is written in terms of the Gauss hypergeometric function both for a ground initial state and for the general case, and comparisons are made between the first-order perturbation theory results for transition probability for harmonic and Morse oscillators. These results are verified with quantum dynamical simulations using wave-packet propagations on a numerical grid. The transition-probability expression incorporating a high-frequency quantum mode and low-frequency medium mode is compared for Morse and harmonic oscillators in different temperature ranges and with various coarse-graining treatments of the delta function from the Fermi golden rule expression. We find that significant deviations from the harmonic approximation are expected for even moderately anharmonic quantum modes at large values of exoergicity. The addition of a second quantum mode of opposite displacement negates the anharmonic effect at small energy change, but in the inverted regime a significantly flatter dependence on exoergicity is predicted for anharmonic modes.     

Probabilistic evolution approach to the expectation value dynamics of quantum mechanical operators,part I: integral representation of Kronecker power series and multivariate Hausdorff moment problems

This is the first one of two companion papers focusing on the establishment of a new path for the expectation value dynamics of the quantum mechanical operators. The main goal of these studies is to do quantum mechanics without explicitly solving Schrödinger wave equation, in other words, without using wave functions except their initially given forms. This goal is achieved by using Ehrenfest theorem and utilizing probabilistic evolution approach (PEA). PEA, first introduced by Metin Demiralp, is a method providing solutions to the nonlinear ordinary differential equations by transforming them to a set of linear ODEs at the cost of denumerably infinite dimensionality. It is recently shown that this method produces analytic solutions, if the initial conditions are given appropriately at some special cases. However, generalization of these conditions to the quantum mechanical applications is not straightforward due to the dispersion of the quantum mechanical systems. For this purpose, multivariate moment problems for the integral representation of the Kronecker power series are introduced and then solved yielding to more specific and precise convergence analysis for the quantum mechanical applications.     

Impacts of quantization on the properties of liquid water

The results of classical and quantum simulations of liquid water over a wide range of temperatures are compared to probe the impact of quantization on the properties of liquid water. We show that, when treated quantum mechanically, water molecules have an enhanced probability of accessing nontetrahedral coordination in the local three-dimensional structure. We discuss how this enhanced probability, also called "effective tunneling", is related to the dynamics of the hydrogen-bond breaking and molecular diffusion in the liquid. We explore in detail how local molecular environments affect the manifestation of quantum effects and identify a previously unreported and apparently unique behavior of the quantum mechanical uncertainty of the water molecule as a function of temperature. The nonmonotonic behavior of the quantum mechanical uncertainty with temperature is shown to be due to the notable strength of the water-water interaction in the condensed phase and becomes further evidence of the importance of the water structure in the properties of this ubiquitous liquid.     

IR spectra of N-methylacetamide in water predicted by combined quantum mechanical/molecular mechanical molecular dynamics simulations

We applied the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulation method in assessing IR spectra of N-methylacetamide and its deuterated form in aqueous solutions. The model peptide is treated at the Austin Model 1 (AM1) level and the induced dipole effects by the solvent are incorporated in fluctuating solute dipole moments, which are calculated using partial charges from Mulliken population analyses without resorting to any available high-level ab initio dipole moment data. Fourier transform of the solute dipole autocorrelation function produces in silico IR spectra, in which the relative peak intensities and bandwidths of major amide bands are quantitatively compatible with experimental results only when both geometric and electronic polarizations of the peptide by the solvent are dealt with at the same quantum-mechanical level. We cast light on the importance of addressing dynamic charge fluctuations of the solute in calculating IR spectra by comparing classical and QM/MM MD simulation results. We propose the adjustable scaling factors for each amide mode to be directly compared with experimental data.     

The flux-flux correlation function for anharmonic barriers

The flux-flux correlation function formalism is a standard and widely used approach for the computation of reaction rates. In this paper we introduce a method to compute the classical and quantum flux-flux correlation functions for anharmonic barriers essentially analytically through the use of the classical and quantum normal forms. In the quantum case we show that for a general f degree-of-freedom system having an index one saddle the quantum normal form reduces the computation of the flux-flux correlation function to that of an effective one-dimensional anharmonic barrier. The example of the computation of the quantum flux-flux correlation function for a fourth order anharmonic barrier is worked out in detail, and we present an analytical expression for the quantum mechanical microcanonical flux-flux correlation function. We then give a discussion of the short-time and harmonic limits.