Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFirst-Principles Semiclassical Initial Value Representation Molecular Dynamics
128
0
0.0
(
0
)
Added by
Al\\'an Aspuru-Guzik
Publication date
2009
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
A method for carrying out semiclassical initial value representation calculations using first-principles molecular dynamics (FP-SC-IVR) is presented. This method can extract the full vibrational power spectrum of carbon dioxide from a single trajectory providing numerical results that agree with experiment even for Fermi resonant states. The computational demands of the method are comparable to those of classical single-trajectory calculations, while describing uniquely quantum features such as the zero-point energy and Fermi resonances. By propagating the nuclear degrees of freedom using first-principles Born-Oppenheimer molecular dynamics, the stability of the method presented is improved considerably when compared to dynamics carried out using fitted potential energy surfaces and numerical derivatives.
rate research
Read More
We extend the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), a semiclassical method for computing real-time correlation functions, to electronically nonadiabatic systems using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian to treat electronic and nuclear degrees of freedom (dofs) within a consistent dynamic framework. We introduce an efficient symplectic integration scheme, the MInt algorithm, for numerical time-evolution of the nuclear and electronic phase space variables as well as the Monodromy matrix, under the non-separable MMST Hamiltonian. We then calculate the probability of transmission through a curve-crossing in model two-level systems and show that in the quantum limit MQC-IVR is in good agreement with the exact quantum results, whereas in the classical limit the method yields results in keeping with mean-field approaches like the Linearized Semiclassical IVR. Finally, exploiting the ability of MQC-IVR to quantize different dofs to different extents, we present a detailed study of the extents to which quantizing the nuclear and electronic dofs improves numerical convergence properties without significant loss of accuracy.
We present in detail and validate an effective Monte Carlo approach for the calculation of the nuclear vibrational densities via integration of molecular eigenfunctions that we have preliminary employed to calculate the densities of the ground and the excited OH stretch vibrational states in protonated glycine molecule [C. Aieta et. al. Nat. Commun. 11, 4348 (2020)]. Here, we first validate and discuss in detail the features of the method on a benchmark water molecule. Then, we apply it to calculate on-the-fly the ab initio anharmonic nuclear densities in correspondence of the fundamental transitions of NH and CH stretches in protonated glycine. We show how we can gain both qualitative and quantitative physical insight by inspection of different one-nucleus densities and assign a character to spectroscopic absorption peaks using the expansion of vibrational states in terms of harmonic basis functions. The visualization of the nuclear vibrations in a purely quantum picture allows us to observe and quantify the effects of anharmonicity on the molecular structure, and to exploit the effect of IR excitations on specific bonds or functional groups, beyond the harmonic approximation. We also calculate the quantum probability distribution of bond-lengths, angles and dihedrals of the molecule. Notably, we observe how in the case of one type of fundamental NH stretching the typical harmonic nodal pattern is absent in the anharmonic distribution.
Chemisorbed molecules at a fuel cell electrode are a very sensitive probe of the surrounding electrochemical environment, and one that can be accurately monitored with different spectroscopic techniques. We develop a comprehensive electrochemical model to study molecular chemisorption at either constant charge or fixed applied voltage, and calculate from first principles the voltage dependence of vibrational frequencies -- the vibrational Stark effect -- for CO adsorbed on close-packed platinum electrodes. The predicted vibrational Stark slopes are found to be in very good agreement with experimental electrochemical spectroscopy data, thereby resolving previous controversies in the quantitative interpretation of in-situ experiments and elucidating the relation between canonical and grand-canonicaldescriptions of vibrational surface phenomena.
We discuss the key steps that have to be followed to calculate coherent quantum transport in molecular and atomic-scale systems, making emphasis on the ab-initio Gaussian Embedded Cluster Method recently developed by the authors. We present various results on a simple system such as a clean Au nanocontact and the same nanocontact in the presence of hydrogen that illustrate the applicability of this method in the study and interpretation of a large range of experiments in the field of molecular electronics.
The electrostatic continuum solvent model developed by Fattebert and Gygi is combined with a first-principles formulation of the cavitation energy based on a natural quantum-mechanical definition for the surface of a solute. Despite its simplicity, the cavitation contribution calculated by this approach is found to be in remarkable agreement with that obtained by more complex algorithms relying on a large set of parameters. Our model allows for very efficient Car-Parrinello simulations of finite or extended systems in solution, and demonstrates a level of accuracy as good as that of established quantum-chemistry continuum solvent methods. We apply this approach to the study of tetracyanoethylene dimers in dichloromethane, providing valuable structural and dynamical insights on the dimerization phenomenon.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
