Adiabatic-Molecular Dynamics Generalized Vertical Hessian Approach: A Mixed Quantum Classical Method to Compute Electronic Spectra of Flexible Molecules in the Condensed Phase
Entity
UAM. Departamento de QuímicaPublisher
American Chemical SocietyDate
2020-02-11Citation
Journal of Chemical Theory and Computation 16.2 (2020): 1215-1231ISSN
1549-9618 (print); 1549-9626 (online)Funded by
This work has received funding from the European Union’s Horizon 2020 research and innovation programme MSCA-ITN under grant agreement no. 765266 (LightDyNAmics). Computational resources provided by the Centro de Cálculo Científico at Universidad Autónoma de Madrid (CCC-UAM) and by SCBI (Supercomputing and Bioinformatics) center of Universidad de Málaga are also acknowledged. D.A. and F.A. acknowledge financial support from Spanish “Ministerio de Economía y Competitividad” (project CTQ2015-65816-R). D.A. acknowledges Fundación Ramón Areces (Spain) for funding his postdoctoral stay at ICCOM-CNR Pisa. J.C. and D.A. acknowledge the Pisa Unit of ICCOM-CNR for hospitality.Project
info:eu-repo/H2020-MSCAITN-2019-765266; Gobierno de España. CTQ2015-65816-REditor's Version
https://doi.org/10.1021/acs.jctc.9b01009Subjects
Time-Dependent Density Functional Theory; Simulation; Normal Modes; QuímicaRights
© 2019 American Chemical SocietyAbstract
We present a general mixed quantum classical method that couples classical molecular dynamics (MD) and vibronic models to compute the shape of electronic spectra of flexible molecules in the condensed phase without, in principle, any phenomenological broadening. It is based on a partition of the nuclear motions of the solute + solvent system in "soft" and "stiff" vibrational modes and an adiabatic hypothesis that assumes that stiff modes are much faster than soft ones. In this framework, the spectrum is rigorously expressed as a conformational integral of quantum vibronic spectra along the stiff coordinates only. Soft modes enter at the classical level through the conformational distribution that is sampled with classical MD runs. In each configuration, reduced-dimensionality quadratic Hamiltonians are built in the space of the stiff coordinates only, thanks to a generalization of the Vertical Hessian harmonic model and an iterative application of projectors in internal coordinates to remove soft modes. Quantum vibronic spectra, specific for each sampled configuration of the soft coordinates, are then computed at the desired temperature with efficient time-dependent techniques, and the global spectrum simply arises from their average. For consistency of the whole procedure, classical MD runs are performed with quantum-mechanically derived force fields, parameterized at the same level of theory selected for generating the quadratic Hamiltonians along the stiff coordinates. Application to N-methyl-6-oxyquinolinium betaine in water, dithiophene in ethanol, and cyanidine in water is presented to show the performance of the method
Files in this item
Google Scholar:Cerezo Bastida, Javier
-
Aranda, Daniel
-
Ávila Ferrer, Francisco José
-
Prampolini, Giacomo
-
Santoro, Fabrizio
This item appears in the following Collection(s)
Related items
Showing items related by title, author, creator and subject.