Search

We investigate the coarse-graining of host-guest systems under the perspective of the local distribution of pore occupancies, along with the physical meaning and actual computability of the coarse-interaction terms. We show that the widely accepted approach, in which the contributions to the free energy given by the molecules located in two neighboring pores are estimated through Monte Carlo simulations where the two pores are kept separated from the rest of the system, leads to inaccurate results at high sorbate densities. In the coarse-graining strategy that we propose, which is based on the Bethe-Peierls approximation, density-independent interaction terms are instead computed according to local effective potentials that take into account the correlations between the pore pair and its surroundings by means of mean-field correction terms without the need for simulating the pore pair separately. Use of the interaction parameters obtained this way allows the coarse-grained system to reproduce more closely the equilibrium properties of the original one. Results are shown for lattice-gases where the local free energy can be computed exactly and for a system of Lennard-Jones particles under the effect of a static confining field. [ABSTRACT FROM AUTHOR]

A very simple equation, FAppij=[(∂2(Va i-Va j)/∂Q2)/(Va i-Va j)]1/2/2, giving a reliable magnitude of non-adiabatic coupling terms (NACTs, Fij's) based on adiabatic potential energies only (Va i and Va j) was discovered, and its reliability was tested for several prototypes of same-symmetry interstate crossings in LiF, C2, NH3Cl, and C6H5SH molecules. Our theoretical derivation starts from the analysis of the relationship between the Lorentzian dependence of NACTs along a diabatization coordinate and the well-established linear vibronic coupling scheme. This analysis results in a very simple equation, α=2κ/Δc, enabling the evaluation of the Lorentz function αα parameter in terms of the coupling constant κκ and the energy gap Δc (Δc= Va i-Va jQc) between adiabatic states at the crossing point Qc. Subsequently, it was shown that Qc corresponds to the point where FAppij exhibit maximum values if we set the coupling parameter as κ=[(Va i-Va j)⋅(∂2(Va i-Va j)/∂Q2)]1/2Qc/2. Finally, we conjectured that this relation could give reasonable values of NACTs not only at the crossing point but also at other geometries near Qc. In this final approximation, the pre-defined crossing point Qc is not required. The results of our test demonstrate that the approximation works much better than initially expected. The present new method does not depend on the selection of an ab initio method for adiabatic electronic states but is currently limited to local non-adiabatic regions where only two electronic states are dominantly involved within a nuclear degree of freedom. [ABSTRACT FROM AUTHOR]

We formulate mixed semiclassical-classical (SC-Cl) propagators by adding a further approximation to the phase-space SC propagators, which have been formulated in our previous paper [S. Koda, J. Chem. Phys. 143, 244110 (2015)]. We first show that the stationary phase approximation over the operation of the phase-space van Vleck propagator on initial distribution functions results in the classical mechanical time propagation. Then, after dividing the degrees of freedom (DOFs) of the total system into the semiclassical DOFs and the classical DOFs, the SC-Cl van Vleck propagator and the SC-Cl Herman-Kluk (HK) propagator are derived by performing the stationary phase approximation only with respect to the classical DOFs. These SC-Cl propagators are naturally decomposed to products of the phase-space SC propagators and the classical mechanical propagators when the system does not have any interaction between the semiclassical and the classical DOFs. In addition, we also numerically compare the original phase-space HK (full HK) propagator and the SC-Cl HK propagator in terms of accuracy and efficiency to find that the accuracy of the SC-Cl HK propagator can be comparable to that of the full HK propagator although the latter is more accurate than the former in general. On the other hand, we confirm that the convergence speed of the SC-Cl HK propagator is faster than that of the full HK propagator. The present numerical tests indicate that the SC-Cl HK propagator can be more accurate than the full HK propagator when they use a same and finite number of classical trajectories due to the balance of the accuracy and the efficiency. [ABSTRACT FROM AUTHOR]

A mixed quantal-semiquantal theory is presented in which the semiquantal squeezed-state wave packet describes the heavy degrees of freedom. Starting from the mean-field equations of motion that are naturally derived from the time-dependent variational principle, we introduce the stochastic particle description for both the quantal and semiquantal parts in an aim to take into account the interparticle correlation, in particular the "quantum backreaction" beyond the mean-field approximation. A numerical application on a model of O2 scattering from a Pt surface demonstrates that the proposed scheme gives correct asymptotic behavior of the scattering probability, with improvement over the mixed quantum-classical scheme with Bohmian particles, which is comprehended by comparing the Bohmian and the stochastic trajectories. [ABSTRACT FROM AUTHOR]

Here, I present a practical recipe for obtaining contravariant vibration-rotation metric tensors, and thus the kinetic energy operators, when some degrees of freedom are constrained rigidly. An element of the contravariant metric tensor is obtained as a sum of dot products of contravariant measuring vectors, which are obtained from their unconstrained counterparts by adding a frozen mode correction. The present method applies in principle for any choice of shape coordinates and a body-frame for which the contravariant measuring vectors can be evaluated. In contrast to the existing methods, the present method does not involve evaluation of covariant metric tensors, matrix inversions, chain rules of derivation, or numerical differentiation. It is applied in the sequel paper [L. Partanen, J. Pesonen, E. Sjöholm, and L. Halonen, J. Chem. Phys. 139, 144311 (2013)] to study the effects of several different approximations to the kinetic energy operator, when the two large-amplitude OH-torsional motions in H2SO4 are of interest. [ABSTRACT FROM AUTHOR]

Expressions for reaction rate constants in multidimensional chemical systems are derived by applying semiclassical approximation to the quantum path integrals of the ImF formulation of reaction rate theory. First, the transverse degrees of freedom orthogonal to the reaction coordinate are treated within the steepest descent approximation, after which the semiclassical approximation is applied to the remaining reaction coordinate. Thus derived, the semiclassical expressions account for the multidimensional nature of quantum effects and accurately incorporate nuclear quantum effects such as multidimensional tunneling and zero point energies. The obtained expressions are applicable in the broad temperature range from the deep tunneling to high-temperature regimes. The present paper provides derivation of the semiclassical instanton expressions proposed by Kryvohuz [J. Chem. Phys. 134, 114103 (2011)]. [ABSTRACT FROM AUTHOR]

Mixed quantum/classical (MQC) simulations treat the majority of a system classically and reserve quantum mechanics only for a few degrees of freedom that actively participate in the chemical process(es) of interest. In MQC calculations, the quantum and classical degrees of freedom are coupled together using pseudopotentials. Although most pseudopotentials are developed empirically, there are methods for deriving pseudopotentials using the results of quantum chemistry calculations, which guarantee that the explicitly-treated valence electron wave functions remain orthogonal to the implicitly-treated core electron orbitals. Whether empirical or analytically derived in nature, to date all such pseudopotentials have been subject to the frozen core approximation (FCA) that ignores how changes in the nuclear coordinates alter the core orbitals, which in turn affects the wave function of the valence electrons. In this paper, we present a way to go beyond the FCA by developing pseudopotentials that respond to these changes. In other words, we show how to derive an analytic expression for a pseudopotential that is an explicit function of nuclear coordinates, thus accounting for the polarization effects experienced by atomic cores in different chemical environments. We then use this formalism to develop a coordinate-dependent pseudopotential for the bonding electron of the sodium dimer cation molecule and we show how the analytic representation of this potential can be used in one-electron MQC simulations that provide the accuracy of a fully quantum mechanical Hartree-Fock (HF) calculation at all internuclear separations. We also show that one-electron MQC simulations of Na2+ using our coordinate-dependent pseudopotential provide a significant advantage in accuracy compared to frozen core potentials with no additional computational expense. This is because use of a frozen core potential produces a charge density for the bonding electron of Na2+ that is too localized on the molecule, leading to significant overbinding of the valence electron. This means that FCA calculations are subject to inaccuracies of order ∼10% in the calculated bond length and vibrational frequency of the molecule relative to a full HF calculation; these errors are fully corrected by using our coordinate-dependent pseudopotential. Overall, our findings indicate that even for molecules like Na2+, which have a simple electronic structure that might be expected to be well-treated within the FCA, the importance of including the effects of the changing core molecular orbitals on the bonding electrons cannot be overlooked. [ABSTRACT FROM AUTHOR]

The microscopic structure of fluids of simple spheres is well known. However, the constituents of most real-life fluids are non-spherical, leading to a coupling between the rotational and translational degrees of freedom. The structure of simple dense fluids of spheroids - ellipsoids of revolution - was only recently determined by direct experimental techniques [A. P. Cohen, E. Janai, E. Mogilko, A. B. Schofield, and E. Sloutskin, Phys. Rev. Lett. 107, 238301 (2011)]. Using confocal microscopy, it was demonstrated that the structure of these simple fluids cannot be described by hard particle models based on the widely used Percus-Yevick approximation. In this paper, we describe a new protocol for determining the shape of the experimental spheroids, which allows us to expand our previous microscopy measurements of these fluids. To avoid the approximations in the theoretical approach, we have also used molecular dynamics simulations to reproduce the experimental radial distribution functions g(r) and estimate the contribution of charge effects to the interactions. Accounting for these charge effects within the Percus-Yevick framework leads to similar agreement with the experiment. [ABSTRACT FROM AUTHOR]

The vibrational coupled cluster (VCC) equations are analyzed in terms of vibrational Mo\ller-Plesset perturbation theory aiming specifically at the importance of four-mode couplings. Based on this analysis, new VCC methods are derived for the calculation of anharmonic vibrational energies and vibrational spectra using vibrational coupled cluster response theory. It is shown how the effect of four-mode coupling and excitations can be efficiently and accurately described using approximations for their inclusion. Two closely related approaches are suggested. The computational scaling of the so-called VCC[3pt4F] method is not higher than the fifth power in the number of vibrational degrees of freedom when up to four-mode coupling terms are present in the Hamiltonian and only fourth order when only up to three-mode couplings are present. With a further approximation, one obtains the VCC[3pt4] model which is shown to scale with at most the fourth power in the number of vibrational degrees of freedom for Hamiltonians with both three- and four-mode coupling levels, while sharing the most important characteristics with VCC[3pt4F]. Sample calculations reported for selected tetra-atomic molecules as well as the larger dioxirane and ethylene oxide molecules support that the new models are accurate and useful. [ABSTRACT FROM AUTHOR]

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 Schrö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. [ABSTRACT FROM AUTHOR]

We present a new approach to calculate real-time quantum dynamics in complex systems. The formalism is based on the partitioning of a system's environment into 'core' and 'reservoir' modes with the former to be treated quantum mechanically and the latter classically. The presented method only requires the calculation of the system's reduced density matrix averaged over the quantum core degrees of freedom which is then coupled to a classically evolved reservoir to treat the remaining modes. We demonstrate our approach by applying it to the spin-boson problem using the noninteracting blip approximation to treat the system and core, and Ehrenfest dynamics to treat the reservoir. The resulting hybrid methodology is accurate for both fast and slow baths, since it naturally reduces to its composite methods in their respective regimes of validity. In addition, our combined method is shown to yield good results in intermediate regimes where neither approximation alone is accurate and to perform equally well for both strong and weak system-bath coupling. Our approach therefore provides an accurate and efficient methodology for calculating quantum dynamics in complex systems. [ABSTRACT FROM AUTHOR]

We derive a many-site version of the non-Markovian time-convolutionless polaron master equation [Jang et al., J. Chem Phys. 129, 101104 (2008)] to describe electronic excitation dynamics in multichromophoric systems. By treating electronic and vibrational degrees of freedom in a combined frame (polaron frame), this theory is capable of interpolating between weak and strong exciton-phonon coupling and is able to account for initial non-equilibrium bath states and spatially correlated environments. Besides outlining a general expression for the expected value of any electronic system observable in the original frame, we also discuss implications of the Markovian and Secular approximations highlighting that they need not hold in the untransformed frame despite being strictly satisfied in the polaron frame. The key features of the theory are illustrated using as an example a four-site subsystem of the Fenna-Mathews-Olson light-harvesting complex. For a spectral density including a localised mode, we show that oscillations of site populations may only be observed when non-equilibrium bath effects are taken into account. Furthermore, we illustrate how this formalism allows us to identify the electronic and vibrational components of the oscillatory dynamics. [ABSTRACT FROM AUTHOR]

The solution of the time-dependent Schrödinger equation for systems of interacting electrons is generally a prohibitive task, for which approximate methods are necessary. Popular approaches, such as the time-dependent Hartree-Fock (TDHF) approximation and time-dependent density functional theory (TDDFT), are essentially single-configurational schemes. TDHF is by construction incapable of fully accounting for the excited character of the electronic states involved in many physical processes of interest; TDDFT, although exact in principle, is limited by the currently available exchange-correlation functionals. On the other hand, multiconfigurational methods, such as the multiconfigurational time-dependent Hartree-Fock (MCTDHF) approach, provide an accurate description of the excited states and can be systematically improved. However, the computational cost becomes prohibitive as the number of degrees of freedom increases, and thus, at present, the MCTDHF method is only practical for few-electron systems. In this work, we propose an alternative approach which effectively establishes a compromise between efficiency and accuracy, by retaining the smallest possible number of configurations that catches the essential features of the electronic wavefunction. Based on a time-dependent variational principle, we derive the MCTDHF working equation for a multiconfigurational expansion with fixed coefficients and specialise to the case of general open-shell states, which are relevant for many physical processes of interest. [ABSTRACT FROM AUTHOR]

Internal conversion is an inherently quantum mechanical process. To date, 'on the fly' computation of internal conversion rates is limited to harmonic approximations, which would seem to be especially unsuitable, given that the typical transition to the ground electronic state occurs at energies which are far from the harmonic limit. It is thus of interest to study the applicability of the semiclassial initial value representation (SCIVR) approach which is in principle amenable to on the fly studies even with 'many' degrees of freedom. In this paper we study the applicability of the Herman-Kluk (HK) SCIVR to a model system with two coupled and anharmonic degrees of freedom. We find that (a) the HK SCIVR is a good approximation to the exact quantum dynamics; (b) computation of the first order correction to the HK-SCIVR approximation corroborates the accuracy; (c) by studying a large parameter range, we find that the harmonic approximation is mostly unsatisfactory; and (d) for the specific model used, the coupling between the modes was found to be relatively unimportant. These results imply that the HK-SCIVR methodology is a good candidate for on the fly studies of internal conversion processes of 'large' molecules. [ABSTRACT FROM AUTHOR]

We investigate the structure formation of amphiphilic molecules at planar walls using density functional theory. The molecules are modeled as (hard) spheres composed of a hydrophilic and hydrophobic part. The orientation of the resulting Janus particles is described as a vector representing an internal degree of freedom. Our density functional approach involves fundamental measure theory combined with a mean-field approximation for the anisotropic interaction. Considering neutral, hydrophilic, and hydrophobic walls, we study the adsorption of the particles, focusing on the competition between the surface field and the interaction-induced ordering phenomena. Finally, we consider systems confined between two planar walls. It is shown that the anisotropic Janus interaction yields pronounced frustration effects at low temperatures. [ABSTRACT FROM AUTHOR]

The frozen Gaussian approximation to the quantum propagator may be a viable method for obtaining 'on the fly' quantum dynamical information on systems with many degrees of freedom. However, it has two severe limitations, it rapidly loses normalization and one needs to know the Gaussian averaged potential, hence it is not a purely local theory in the force field. These limitations are in principle remedied by using the Herman-Kluk (HK) form for the semiclassical propagator. The HK propagator approximately conserves unitarity for relatively long times and depends only locally on the bare potential and its second derivatives. However, the HK propagator involves a much more expensive computation due to the need for evaluating the monodromy matrix elements. In this paper, we (a) derive a new formula for the normalization integral based on a prefactor free HK propagator which is amenable to 'on the fly' computations; (b) show that a frozen Gaussian version of the normalization integral is not readily computable 'on the fly'; (c) provide a new insight into how the HK prefactor leads to approximate unitarity; and (d) how one may construct a prefactor free approximation which combines the advantages of the frozen Gaussian and the HK propagators. The theoretical developments are backed by numerical examples on a Morse oscillator and a quartic double well potential. [ABSTRACT FROM AUTHOR]

We introduce a new formula for the acceleration weight factor in the hyperdynamics simulation method, the use of which correctly provides an exact simulation of the true dynamics of a system. This new form of hyperdynamics is valid and applicable where the transition state theory (TST) is applicable and also where the TST is not applicable. To illustrate this new formulation, we perform hyperdynamics simulations for four systems ranging from one degree of freedom to 591 degrees of freedom: (1) We first analyze free diffusion having one degree of freedom. This system does not have a transition state. The TST and the original form of hyperdynamics are not applicable. Using the new form of hyperdynamics, we compute mean square displacement for a range of time. The results obtained agree perfectly with the analytical formula. (2) Then we examine the classical Kramers escape rate problem. The rate computed is in perfect agreement with the Kramers formula over a broad range of temperature. (3) We also study another classical problem: Computing the rate of effusion out of a cubic box through a tiny hole. This problem does not involve an energy barrier. Thus, the original form of hyperdynamics excludes the possibility of using a nonzero bias and is inappropriate. However, with the new weight factor formula, our new form of hyperdynamics can be easily implemented and it produces the exact results. (4) To illustrate applicability to systems of many degrees of freedom, we analyze diffusion of an atom adsorbed on the (001) surface of an fcc crystal. The system is modeled by an atom on top of a slab of six atomic layers. Each layer has 49 atoms. With the bottom two layers of atoms fixed, this system has 591 degrees of freedom. With very modest computing effort, we are able to characterize its diffusion pathways in the exchange-with-the-substrate and hop-over-the-bridge mechanisms. [ABSTRACT FROM AUTHOR]

The forward-backward (FB) approximation as applied to semiclassical initial value representations (SCIVR’s) has enabled the practical application of the SCIVR methodology to systems with many degrees of freedom. However, to date a systematic representation of the exact quantum dynamics in terms of the FB-SCIVR has proven elusive. In this paper, we provide a new derivation of a forward-backward phase space SCIVR expression (FBPS-SCIVR) derived previously by Thompson and Makri [Phys. Rev. E 59, R4729 (1999)]. This enables us to represent quantum correlation functions exactly in terms of a series whose leading order term is the FBPS-SCIVR expression. Numerical examples for systems with over 50 degrees of freedom are presented for the spin boson problem. Comparison of the FBPS-SCIVR with the numerically exact results of Wang [J. Chem. Phys. 113, 9948 (2000)] obtained using a multiconfigurational time dependent method shows that the leading order FBPS-SCIVR term already provides an excellent approximation. [ABSTRACT FROM AUTHOR]

A semiclassical (SC) approximation to the quantum mechanical propagator for nonadiabatic systems is derived. Our derivation starts with an exact path integral expression that uses canonical coherent states for the nuclear degrees of freedom and spin coherent states for the electronic degrees of freedom. A stationary path approximation (SPA) is then applied to the path integral to obtain the SC approximation. The SPA results in complex classical trajectories of both nuclear and electronic degrees of freedom and a double ended boundary condition. The root search problem is solved using the previously proposed “real trajectory local search” algorithm. The SC approximation is tested on three simple one dimensional two-state systems proposed by Tully [J. Chem. Phys. 93, 1061 (1990)], and the SC results are compared to Ehrenfest and surface hopping predictions. Excellent agreement with quantum results is reached when the SC trajectory is far away from caustics. We discuss the origin of caustics in this SC formalism and the strengths and weaknesses of this approach. [ABSTRACT FROM AUTHOR]

The field-free orientation induced by a single-cycle terahertz (THz) laser pulse is studied for two 'heavy' molecules, NaI and IBr. Two methods are used and compared in the calculations: One is to solve the exact time-dependent Schrödinger equation (ETDSE) considering the full-rovibrational degrees of freedom, and the other is to invoke the rigid-rotor approximation (RRA). Calculations are performed for the central frequency varying from 0.05 to 1.0 THz and for the peak intensity taken to be 5 × 107, 2 × 108, and 5 × 108 W/cm2, respectively. The degree of field-free orientation, , is strongly dependent on the central frequency and the peak intensity of the single-cycle THz pulse. The maximum degree of field-free orientation is determined to be 0.84 for NaI and 0.63 for IBr in these given ranges of frequency and intensity. The molecular orientation obtained by the RRA calculations is in good agreement with that obtained by the ETDSE in the given parameter region. [ABSTRACT FROM AUTHOR]

Powerful approximate methods for propagating the density matrix of complex systems that are conveniently described in terms of electronic subsystem states and nuclear degrees of freedom have recently been developed that involve linearizing the density matrix propagator in the difference between the forward and backward paths of the nuclear degrees of freedom while keeping the interference effects between the different forward and backward paths of the electronic subsystem described in terms of the mapping Hamiltonian formalism and semi-classical mechanics. Here we demonstrate that different approaches to developing the linearized approximation to the density matrix propagator can yield a mean-field like approximate propagator in which the nuclear variables evolve classically subject to Ehrenfest-like forces that involve an average over quantum subsystem states, and by adopting an alternative approach to linearizing we obtain an algorithm that involves classical like nuclear dynamics influenced by a quantum subsystem state dependent force reminiscent of trajectory surface hopping methods. We show how these different short time approximations can be implemented iteratively to achieve accurate, stable long time propagation and explore their implementation in different representations. The merits of the different approximate quantum dynamics methods that are thus consistently derived from the density matrix propagator starting point and different partial linearization approximations are explored in various model system studies of multi-state scattering problems and dissipative non-adiabatic relaxation in condensed phase environments that demonstrate the capabilities of these different types of approximations for treating non-adiabatic electronic relaxation, bifurcation of nuclear distributions, and the passage from nonequilibrium coherent dynamics at short times to long time thermal equilibration in the presence of a model dissipative environment. [ABSTRACT FROM AUTHOR]

The semiclassical instanton approach discussed by Kryvohuz [J. Chem. Phys. 134, 114103 (2011)] is applied to calculate kinetic H/D isotope effect (KIE) of intramolecular hydrogen transfer in cis-1,3-pentadiene. All 33 vibrational degrees of freedom are treated quantum mechanically with semiclassical approximation. Nuclear quantum effects such as tunneling under the barrier and zero-point energy are automatically incorporated in the theory, and are shown to be responsible for the observed appreciable kinetic isotope effect in cis-1,3-pentadiene. Over the barrier passage is also automatically included. Numerical calculations are performed on an empirical valence bond potential energy surface and compared with the previous experimental and theoretical studies. An estimation of heavy-atom 12C/13C KIE in the same system is also provided and the factors contributing to it are discussed. [ABSTRACT FROM AUTHOR]

The approach of defining quantum corrections on nuclear dynamics of molecular systems incorporated approximately into selected degrees of freedom, is described. The approach is based on the Madelung-de-Broglie-Bohm formulation of time-dependent quantum mechanics which represents a wavefunction in terms of an ensemble of trajectories. The trajectories follow classical laws of motion except that the quantum potential, dependent on the wavefunction amplitude and its derivatives, is added to the external, classical potential. In this framework the quantum potential, determined approximately for practical reasons, is included only into the 'quantum' degrees of freedom describing light particles such as protons, while neglecting with the quantum force for the heavy, nearly classical nuclei. The entire system comprised of light and heavy particles is described by a single wavefunction of full dimensionality. The coordinate space of heavy particles is divided into spatial domains or subspaces. The quantum force acting on the light particles is determined for each domain of similar configurations of the heavy nuclei. This approach effectively introduces parametric dependence of the reduced dimensionality quantum force, on classical degrees of freedom. This strategy improves accuracy of the quantum force and does not restrict interaction between the domains. The concept is illustrated for two-dimensional scattering systems, where the quantum force is required to reproduce vibrational energy of the quantum degree of freedom. [ABSTRACT FROM AUTHOR]