1
|
Solov’yov AV, Verkhovtsev AV, Mason NJ, Amos RA, Bald I, Baldacchino G, Dromey B, Falk M, Fedor J, Gerhards L, Hausmann M, Hildenbrand G, Hrabovský M, Kadlec S, Kočišek J, Lépine F, Ming S, Nisbet A, Ricketts K, Sala L, Schlathölter T, Wheatley AEH, Solov’yov IA. Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment. Chem Rev 2024; 124:8014-8129. [PMID: 38842266 PMCID: PMC11240271 DOI: 10.1021/acs.chemrev.3c00902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2023] [Revised: 05/02/2024] [Accepted: 05/10/2024] [Indexed: 06/07/2024]
Abstract
This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.
Collapse
Affiliation(s)
| | | | - Nigel J. Mason
- School
of Physics and Astronomy, University of
Kent, Canterbury CT2 7NH, United
Kingdom
| | - Richard A. Amos
- Department
of Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, U.K.
| | - Ilko Bald
- Institute
of Chemistry, University of Potsdam, Karl-Liebknecht-Str. 24-25, 14476 Potsdam, Germany
| | - Gérard Baldacchino
- Université
Paris-Saclay, CEA, LIDYL, 91191 Gif-sur-Yvette, France
- CY Cergy Paris Université,
CEA, LIDYL, 91191 Gif-sur-Yvette, France
| | - Brendan Dromey
- Centre
for Light Matter Interactions, School of Mathematics and Physics, Queen’s University Belfast, Belfast BT7 1NN, United Kingdom
| | - Martin Falk
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, 61200 Brno, Czech Republic
- Kirchhoff-Institute
for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
| | - Juraj Fedor
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Luca Gerhards
- Institute
of Physics, Carl von Ossietzky University, Carl-von-Ossietzky-Str. 9-11, 26129 Oldenburg, Germany
| | - Michael Hausmann
- Kirchhoff-Institute
for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
| | - Georg Hildenbrand
- Kirchhoff-Institute
for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
- Faculty
of Engineering, University of Applied Sciences
Aschaffenburg, Würzburger
Str. 45, 63743 Aschaffenburg, Germany
| | | | - Stanislav Kadlec
- Eaton European
Innovation Center, Bořivojova
2380, 25263 Roztoky, Czech Republic
| | - Jaroslav Kočišek
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Franck Lépine
- Université
Claude Bernard Lyon 1, CNRS, Institut Lumière
Matière, F-69622, Villeurbanne, France
| | - Siyi Ming
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield
Road, Cambridge CB2 1EW, United Kingdom
| | - Andrew Nisbet
- Department
of Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, U.K.
| | - Kate Ricketts
- Department
of Targeted Intervention, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Leo Sala
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Thomas Schlathölter
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh
4, 9747 AG Groningen, The Netherlands
- University
College Groningen, University of Groningen, Hoendiepskade 23/24, 9718 BG Groningen, The Netherlands
| | - Andrew E. H. Wheatley
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield
Road, Cambridge CB2 1EW, United Kingdom
| | - Ilia A. Solov’yov
- Institute
of Physics, Carl von Ossietzky University, Carl-von-Ossietzky-Str. 9-11, 26129 Oldenburg, Germany
| |
Collapse
|
2
|
Huang X, Zhang W, Liang W. Time-dependent Kohn-Sham electron dynamics coupled with nonequilibrium plasmonic response via atomistic electromagnetic model. J Chem Phys 2024; 160:214106. [PMID: 38828813 DOI: 10.1063/5.0205845] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2024] [Accepted: 05/20/2024] [Indexed: 06/05/2024] Open
Abstract
Computational modeling of plasmon-mediated molecular photophysical and photochemical behaviors can help us better understand and tune the bound molecular properties and reactivity and make better decisions to design and control nanostructures. However, computational investigations of coupled plasmon-molecule systems are challenging due to the lack of accurate and efficient protocols to simulate these systems. Here, we present a hybrid scheme by combining the real-time time-dependent density functional theory (RT-TDDFT) approach with the time-domain frequency dependent fluctuating charge (TD-ωFQ) model. At first, we transform ωFQ in the frequency-domain, an atomistic electromagnetic model for the plasmonic response of plasmonic metal nanoparticles (PMNPs), into the time-domain and derive its equation-of-motion formulation. The TD-ωFQ introduces the nonequilibrium plasmonic response of PMNPs and atomistic interactions to the electronic excitation of the quantum mechanical (QM) region. Then, we combine TD-ωFQ with RT-TDDFT. The derived RT-TDDFT/TD-ωFQ scheme allows us to effectively simulate the plasmon-mediated "real-time" electronic dynamics and even the coupled electron-nuclear dynamics by combining them with the nuclear dynamics approaches. As a first application of the RT-TDDFT/TD-ωFQ method, we study the nonradiative decay rate and plasmon-enhanced absorption spectra of two small molecules in the proximity of sodium MNPs. Thanks to the atomistic nature of the ωFQ model, the edge effect of MNP on absorption enhancement has also been investigated and unveiled.
Collapse
Affiliation(s)
- Xunkun Huang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Wenshu Zhang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - WanZhen Liang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| |
Collapse
|
3
|
Bondanza M, Nottoli T, Nottoli M, Cupellini L, Lipparini F, Mennucci B. The OpenMMPol library for polarizable QM/MM calculations of properties and dynamics. J Chem Phys 2024; 160:134106. [PMID: 38557842 DOI: 10.1063/5.0198251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Accepted: 03/13/2024] [Indexed: 04/04/2024] Open
Abstract
We present a new library designed to provide a simple and straightforward way to implement QM/AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) and other polarizable QM/MM (Molecular Mechanics) methods based on induced point dipoles. The library, herein referred to as OpenMMPol, is free and open-sourced and is engineered to address the increasing demand for accurate and efficient QM/MM simulations. OpenMMPol is specifically designed to allow polarizable QM/MM calculations of ground state energies and gradients and excitation properties. Key features of OpenMMPol include a modular architecture facilitating extensibility, parallel computing capabilities for enhanced performance on modern cluster architectures, a user-friendly interface for intuitive implementation, and a simple and flexible structure for providing input data. To show the capabilities offered by the library, we present an interface with PySCF to perform QM/AMOEBA molecular dynamics, geometry optimization, and excited-state calculation based on (time-dependent) density functional theory.
Collapse
Affiliation(s)
- Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Tommaso Nottoli
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Michele Nottoli
- Institute of Applied Analysis and Numerical Simulation, Universität Stuttgart, Pfaffenwaldring 57, D-70569 Stuttgart, Germany
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
4
|
Salvadori G, Mazzeo P, Accomasso D, Cupellini L, Mennucci B. Deciphering Photoreceptors Through Atomistic Modeling from Light Absorption to Conformational Response. J Mol Biol 2024; 436:168358. [PMID: 37944793 DOI: 10.1016/j.jmb.2023.168358] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Revised: 10/28/2023] [Accepted: 11/02/2023] [Indexed: 11/12/2023]
Abstract
In this review, we discuss the successes and challenges of the atomistic modeling of photoreceptors. Throughout our presentation, we integrate explanations of the primary methodological approaches, ranging from quantum mechanical descriptions to classical enhanced sampling methods, all while providing illustrative examples of their practical application to specific systems. To enhance the effectiveness of our analysis, our primary focus has been directed towards the examination of applications across three distinct photoreceptors. These include an example of Blue Light-Using Flavin (BLUF) domains, a bacteriophytochrome, and the orange carotenoid protein (OCP) employed by cyanobacteria for photoprotection. Particular emphasis will be placed on the pivotal role played by the protein matrix in fine-tuning the initial photochemical event within the embedded chromophore. Furthermore, we will investigate how this localized perturbation initiates a cascade of events propagating from the binding pocket throughout the entire protein structure, thanks to the intricate network of interactions between the chromophore and the protein.
Collapse
Affiliation(s)
- Giacomo Salvadori
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Patrizia Mazzeo
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Davide Accomasso
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Lorenzo Cupellini
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| |
Collapse
|
5
|
Nicoli L, Giovannini T, Cappelli C. Assessing the quality of QM/MM approaches to describe vacuo-to-water solvatochromic shifts. J Chem Phys 2022; 157:214101. [PMID: 36511555 DOI: 10.1063/5.0118664] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The performance of different quantum mechanics/molecular mechanics embedding models to compute vacuo-to-water solvatochromic shifts is investigated. In particular, both nonpolarizable and polarizable approaches are analyzed and computed results are compared to reference experimental data. We show that none of the approaches outperform the others and that errors strongly depend on the nature of the molecular transition to be described. Thus, we prove that the best choice of embedding model highly depends on the molecular system and that the use of a specific approach as a black box can lead to significant errors and, sometimes, totally wrong predictions.
Collapse
Affiliation(s)
- Luca Nicoli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
6
|
Kumari B, Huwaidi A, Robert G, Cloutier P, Bass AD, Sanche L, Wagner JR. Shape Resonances in DNA: Nucleobase Release, Reduction, and Dideoxynucleoside Products Induced by 1.3 to 2.3 eV Electrons. J Phys Chem B 2022; 126:5175-5184. [PMID: 35793462 DOI: 10.1021/acs.jpcb.2c01851] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding the details of DNA damage caused by high-energy particles or photons is complicated by the multitude of reactive species, arising from the ionization and dissociation of H2O, DNA, and protein. In this work, oligonucleotides (ODNs) are irradiated with a beam of low-energy electrons of 1.3 to 2.3 eV, which can only induce damage via the decay of shape resonances into various dissociative electron attachment channels. Using LC-MS/MS analysis, the major products are the release of nonmodified nucleobases (NB; Cyt ≫ Thy ∼ Ade > Gua). Additional damage includes 5,6-dihydropyrimidines (dHT > dHU) and eight nucleosides with modified sugar moieties consisting of 2',3'- and 2',5'-dideoxynucleosides (ddG > ddA ∼ ddC > ddT). The distribution of products is remarkably different in a 16-mer ODN compared to that observed previously with thymidylyl-(3'-5')-thymidine. This difference is explained by electron delocalization occurring within a sufficiently long strand, the DEA theory of O'Malley, and recent time-dependent density functional theory calculations.
Collapse
Affiliation(s)
- Bhavini Kumari
- Département de Médecine Nucléaire et Radiobiologie, Faculté de Médecine et des Sciences de la Santé, Université de Sherbrooke, 3001 12e Avenue Nord, Québec J1H 5N4, Canada
| | - Alaa Huwaidi
- Département de Médecine Nucléaire et Radiobiologie, Faculté de Médecine et des Sciences de la Santé, Université de Sherbrooke, 3001 12e Avenue Nord, Québec J1H 5N4, Canada
| | - Gabriel Robert
- Département de Médecine Nucléaire et Radiobiologie, Faculté de Médecine et des Sciences de la Santé, Université de Sherbrooke, 3001 12e Avenue Nord, Québec J1H 5N4, Canada
| | - Pierre Cloutier
- Département de Médecine Nucléaire et Radiobiologie, Faculté de Médecine et des Sciences de la Santé, Université de Sherbrooke, 3001 12e Avenue Nord, Québec J1H 5N4, Canada
| | - Andrew D Bass
- Département de Médecine Nucléaire et Radiobiologie, Faculté de Médecine et des Sciences de la Santé, Université de Sherbrooke, 3001 12e Avenue Nord, Québec J1H 5N4, Canada
| | - Léon Sanche
- Département de Médecine Nucléaire et Radiobiologie, Faculté de Médecine et des Sciences de la Santé, Université de Sherbrooke, 3001 12e Avenue Nord, Québec J1H 5N4, Canada
| | - J Richard Wagner
- Département de Médecine Nucléaire et Radiobiologie, Faculté de Médecine et des Sciences de la Santé, Université de Sherbrooke, 3001 12e Avenue Nord, Québec J1H 5N4, Canada
| |
Collapse
|
7
|
Suzuki K, Maeda S. Multistructural microiteration combined with QM/MM-ONIOM electrostatic embedding. Phys Chem Chem Phys 2022; 24:16762-16773. [PMID: 35775395 DOI: 10.1039/d2cp02270b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Multistructural microiteration (MSM) is a method to take account of contributions of multiple surrounding structures in a geometrical optimization or reaction path calculation using the quantum mechanics/molecular mechanics (QM/MM) ONIOM method. In this study, we combined MSM with the electrostatic embedding (EE) scheme of the QM/MM-ONIOM method by extending its original formulation for mechanical embedding (ME). MSM-EE takes account of the polarization in the QM region induced by point charges assigned to atoms in the multiple surrounding structures, where the point charges are scaled by the weight factor of each surrounding structure determined through MSM. The performance of MSM-EE was compared with that of the other methods, i.e., ONIOM-ME, ONIOM-EE, and MSM-ME, by applying them to three chemical processes: (1) chorismate-to-prephenate transformation in aqueous solution, (2) the same transformation as (1) in an enzyme, and (3) hydroxylation in p-hydroxybenzoate hydroxylase. These numerical tests of MSM-EE yielded barriers and reaction energies close to experimental values with computational costs comparable to those of the other three methods.
Collapse
Affiliation(s)
- Kimichi Suzuki
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo 001-0021, Japan. .,Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan.,JST, ERATO Maeda Artificial Intelligence for Chemical Reaction Design and Discovery Project, Sapporo 060-0810, Japan
| | - Satoshi Maeda
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo 001-0021, Japan. .,Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan.,JST, ERATO Maeda Artificial Intelligence for Chemical Reaction Design and Discovery Project, Sapporo 060-0810, Japan.,Research and Services Division of Materials Data and Integrated System (MaDIS), National Institute for Materials Science (NIMS), Tsukuba 305-0044, Japan
| |
Collapse
|
8
|
Korsaye FA, de la Lande A, Ciofini I. Following the density evolution using real time density functional theory and density based indexes: Application to model push-pull molecules. J Comput Chem 2022; 43:1464-1473. [PMID: 35766295 DOI: 10.1002/jcc.26932] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Revised: 04/19/2022] [Accepted: 04/19/2022] [Indexed: 02/01/2023]
Abstract
Considering as test case a family of organic rod like push-pull molecules, we derived and applied density based index enabling the description and diagnostic of the electronic density evolution in real time-time dependent density functional theory (RT-TDDFT) simulations. In particular, both the charge transfer (CT) distance and a diagnostic index, the DCT and MAC RT respectively, were computed on the fly from the density distribution obtained at a given time and the reference ground state density and their mean values were compared with what obtained at Linear Response-TDDFT level. Besides giving a way of analyzing the density redistribution occurring in time, these tools allowed to show how RT-TDDFT, which is definitely a powerful method to model the evolution of the density in CT or charge separation processes, can be affected by the same artifacts known for LR-TDDFT approaches and, particularly, to those related to the use of approximate exchange correlation functionals. The analysis here performed allowed to identify and discard on fly the electronic configurations corresponding to spurious situations.
Collapse
Affiliation(s)
- Feven Alemu Korsaye
- PSL University, CNRS, Chimie ParisTech-PSL, Institute of Chemistry for Life and Health Sciences (i-CLeHS), Theoretical Chemistry and Modelling Group (CTM), Paris, France.,Institut de Chimie Physique, Université Paris Saclay, CNRS, UMR 8000, Orsay, France
| | - Aurélien de la Lande
- Institut de Chimie Physique, Université Paris Saclay, CNRS, UMR 8000, Orsay, France
| | - Ilaria Ciofini
- PSL University, CNRS, Chimie ParisTech-PSL, Institute of Chemistry for Life and Health Sciences (i-CLeHS), Theoretical Chemistry and Modelling Group (CTM), Paris, France
| |
Collapse
|
9
|
Salahub DR. Multiscale molecular modelling: from electronic structure to dynamics of nanosystems and beyond. Phys Chem Chem Phys 2022; 24:9051-9081. [PMID: 35389399 DOI: 10.1039/d1cp05928a] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Important contemporary biological and materials problems often depend on interactions that span orders of magnitude differences in spatial and temporal dimensions. This Tutorial Review attempts to provide an introduction to such fascinating problems through a series of case studies, aimed at beginning researchers, graduate students, postdocs and more senior colleagues who are changing direction to focus on multiscale aspects of their research. The choice of specific examples is highly personal, with examples either chosen from our own work or outstanding multiscale efforts from the literature. I start with various embedding schemes, as exemplified by polarizable continuum models, 3-D RISM, molecular DFT and frozen-density embedding. Next, QM/MM (quantum mechanical/molecular mechanical) techniques are the workhorse of pm-to-nm/ps-to-ns simulations; examples are drawn from enzymes and from nanocatalysis for oil-sands upgrading. Using polarizable force-fields in the QM/MM framework represents a burgeoning subfield; with examples from ion channels and electron dynamics in molecules subject to strong external fields, probing the atto-second dynamics of the electrons with RT-TDDFT (real-time - time-dependent density functional theory) eventually coupled with nuclear motion through the Ehrenfest approximation. This is followed by a section on coarse graining, bridging dimensions from atoms to cells. The penultimate chapter gives a quick overview of multiscale approaches that extend into the meso- and macro-scales, building on atomistic and coarse-grained techniques to enter the world of materials engineering, on the one hand, and cell biology, on the other. A final chapter gives just a glimpse of the burgeoning impact of machine learning on the structure-dynamics front. I aim to capture the excitement of contemporary leading-edge breakthroughs in the description of physico-chemical systems and processes in complex environments, with only enough historical content to provide context and aid the next generation of methodological development. While I aim also for a clear description of the essence of methodological breakthroughs, equations are kept to a minimum and detailed formalism and implementation details are left to the references. My approach is very selective (case studies) rather than exhaustive. I think that these case studies should provide fodder to build as complete a reference tree on multiscale modelling as the reader may wish, through forward and backward citation analysis. I hope that my choices of cases will excite interest in newcomers and help to fuel the growth of multiscale modelling in general.
Collapse
Affiliation(s)
- Dennis R Salahub
- Department of Chemistry, Department of Physics and Astronomy, CMS-Centre for Molecular Simulation, IQST-Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, Calgary, Alberta, T2N 1N4, Canada.
| |
Collapse
|
10
|
Wozniak AP, Przybytek M, Lewenstein M, Moszynski R. Effects of electronic correlation on the high harmonic generation in helium: a time-dependent configuration interaction singles vs time-dependent full configuration interaction study. J Chem Phys 2022; 156:174106. [DOI: 10.1063/5.0087384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we investigate the effects of full electronic correlation on the high harmonic generation in the helium atom subjected to laser pulses of extremely high intensity. To do this, we perform real-time propagations of the helium atom wavefunction using quantum chemistry methods coupled to Gaussian basis sets. The calculations are done within the real-time time-dependent configuration interaction framework, at two levels of theory: time-dependent configuration interation with single excitations (TD-CIS, uncorrelated method) and time-dependent full configuration interaction (TD-FCI, fully correlated method). The electronic wavefunction is expanded in Dunning basis sets supplemented with functions adapted to describing highly excited and continuum states. We also compare the TD-CI results with grid-based propagations of the helium atom within the single-active-electron approximation. Our results show that when including the dynamical electron correlation, a noticeable improvement to the description of HHG can be achieved, in terms of e.g. a more constant intensity in the lower energy part of the harmonic plateau. However, such effects can be captured only if the basis set used suffices to reproduce the most basic features, such as the HHG cutoff position, at the uncorrelated level of theory.
Collapse
|
11
|
Yusef Buey M, Mineva T, Rapacioli M. Coupling density functional based tight binding with class 1 force fields in a hybrid QM/MM scheme. Theor Chem Acc 2022. [DOI: 10.1007/s00214-022-02878-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
12
|
De Santis M, Vallet V, Gomes ASP. Environment Effects on X-Ray Absorption Spectra With Quantum Embedded Real-Time Time-Dependent Density Functional Theory Approaches. Front Chem 2022; 10:823246. [PMID: 35295974 PMCID: PMC8919347 DOI: 10.3389/fchem.2022.823246] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2021] [Accepted: 02/04/2022] [Indexed: 11/13/2022] Open
Abstract
In this work we implement the real-time time-dependent block-orthogonalized Manby-Miller embedding (rt-BOMME) approach alongside our previously developed real-time frozen density embedding time-dependent density functional theory (rt-TDDFT-in-DFT FDE) code, and investigate these methods’ performance in reproducing X-ray absorption spectra (XAS) obtained with standard rt-TDDFT simulations, for model systems comprised of solvated fluoride and chloride ions ([X@(H2O)8−, X = F, Cl). We observe that for ground-state quantities such as core orbital energies, the BOMME approach shows significantly better agreement with supermolecular results than FDE for the strongly interacting fluoride system, while for chloride the two embedding approaches show more similar results. For the excited states, we see that while FDE (constrained not to have the environment densities relaxed in the ground state) is in good agreement with the reference calculations for the region around the K and L1 edges, and is capable of reproducing the splitting of the 1s1 (n + 1)p1 final states (n + 1 being the lowest virtual p orbital of the halides), it by and large fails to properly reproduce the 1s1 (n + 2)p1 states and misses the electronic states arising from excitation to orbitals with important contributions from the solvent. The BOMME results, on the other hand, provide a faithful qualitative representation of the spectra in all energy regions considered, though its intrinsic approximation of employing a lower-accuracy exchange-correlation functional for the environment induces non-negligible shifts in peak positions for the excitations from the halide to the environment. Our results thus confirm that QM/QM embedding approaches are viable alternatives to standard real-time simulations of X-ray absorption spectra of species in complex or confined environments.
Collapse
|
13
|
Coccia E, Luppi E. Time-dependent ab initioapproaches for high-harmonic generation spectroscopy. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 34:073001. [PMID: 34731835 DOI: 10.1088/1361-648x/ac3608] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Accepted: 11/03/2021] [Indexed: 06/13/2023]
Abstract
High-harmonic generation (HHG) is a nonlinear physical process used for the production of ultrashort pulses in XUV region, which are then used for investigating ultrafast phenomena in time-resolved spectroscopies. Moreover, HHG signal itself encodes information on electronic structure and dynamics of the target, possibly coupled to the nuclear degrees of freedom. Investigating HHG signal leads to HHG spectroscopy, which is applied to atoms, molecules, solids and recently also to liquids. Analysing the number of generated harmonics, their intensity and shape gives a detailed insight of, e.g., ionisation and recombination channels occurring in the strong-field dynamics. A number of valuable theoretical models has been developed over the years to explain and interpret HHG features, with the three-step model being the most known one. Originally, these models neglect the complexity of the propagating electronic wavefunction, by only using an approximated formulation of ground and continuum states. Many effects unravelled by HHG spectroscopy are instead due to electron correlation effects, quantum interference, and Rydberg-state contributions, which are all properly captured by anab initioelectronic-structure approach. In this review we have collected recent advances in modelling HHG by means ofab initiotime-dependent approaches relying on the propagation of the time-dependent Schrödinger equation (or derived equations) in presence of a very intense electromagnetic field. We limit ourselves to gas-phase atomic and molecular targets, and to solids. We focus on the various levels of theory employed for describing the electronic structure of the target, coupled with strong-field dynamics and ionisation approaches, and on the basis used to represent electronic states. Selected applications and perspectives for future developments are also given.
Collapse
Affiliation(s)
- Emanuele Coccia
- Dipartimento di Scienze Chimiche e Farmaceutiche, University of Trieste, via Giorgieri 1, 34127 Trieste, Italy
| | - Eleonora Luppi
- Sorbonne Universités, UPMC Univ Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005 Paris, France
- CNRS, UMR 7616, Laboratoire de Chimie Théorique, F-75005 Paris, France
| |
Collapse
|
14
|
de la Lande A, Denisov S, Mostafavi M. The mystery of sub-picosecond charge transfer following irradiation of hydrated uridine monophosphate. Phys Chem Chem Phys 2021; 23:21148-21162. [PMID: 34528029 DOI: 10.1039/d0cp06482c] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The early mechanisms by which ionizing rays damage biological structures by so-called direct effects are largely elusive. In a recent picosecond pulse radiolysis study of concentrated uridine monophosphate solutions [J. Ma, S. A. Denisov, J.-L. Marignier, P. Pernot, A. Adhikary, S. Seki and M. Mostafavi, J. Phys. Chem. Lett., 2018, 9, 5105], unexpected results were found regarding the oxidation of the nucleobase. The signature of the oxidized nucleobase could not be detected 5 ps after the electron pulse, but only the oxidized phosphate, raising intriguing questions about the identity of charge-transfer mechanisms that could explain the absence of U+. We address here this question by means of advanced first-principles atomistic simulations of solvated uridine monophosphate, combining Density Functional Theory (DFT) with polarizable embedding schemes. We contrast three very distinct mechanisms of charge transfer covering the atto-, femto- and pico-second timescales. We first investigate the ionization mechanism and subsequent hole/charge migrations on a timescale of attoseconds to a few femtoseconds under the frozen nuclei approximation. We then consider a nuclear-driven phosphate-to-oxidized-nucleobase electron transfer, showing that it is an uncompetitive reaction channel on the sub-picosecond timescale, despite its high exothermicity and significant electronic coupling. Finally, we show that non-adiabatic charge transfer is enabled by femtosecond nuclear relaxation after ionization. We show that electronic decoherence and the electronic coupling strength are the key parameters that determine the hopping probabilities. Our results provide important insight into the interplay between electronics and nuclear motions in the early stages of the multiscale responses of biological matter subjected to ionizing radiation.
Collapse
Affiliation(s)
- Aurélien de la Lande
- Institut de Chimie Physique, CNRS, Université Paris Saclay (UMR 8000), 15 Avenue Jean Perrin, 91405, France.
| | - Sergey Denisov
- Institut de Chimie Physique, CNRS, Université Paris Saclay (UMR 8000), 15 Avenue Jean Perrin, 91405, France.
| | - Mehran Mostafavi
- Institut de Chimie Physique, CNRS, Université Paris Saclay (UMR 8000), 15 Avenue Jean Perrin, 91405, France.
| |
Collapse
|
15
|
Semmeq A, Badawi M, Dziurla MA, Ouaskit S, Monari A. Nucleic Acids under Stress: Understanding and Simulating Nucleobase Fragmentation Pathways. Chempluschem 2021; 86:1426-1435. [PMID: 34637193 DOI: 10.1002/cplu.202100323] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2021] [Revised: 09/24/2021] [Indexed: 11/10/2022]
Abstract
The effects of radiations on nucleic acids and their constituents is widely studied across several research fields using different experimental and theoretical protocols. While a large number of studies were performed in this context, many fundamental physical and chemical effects are still being investigated, particularly involving the effect of the biological environment. As an example, the interpretation of experimental nucleic acid bases mass spectra, and hence inferring their reactivity in complex environment still poses great challenge. This Minireview summarizes recent theoretical advancements aiming to predict and interpret the reactivity of nucleic acid bases. We focus not only on the understanding of the inherent fragmentation pathways of isolated nucleobases but also on the modeling of a realistic nano-environments highlighting the importance of molecular dynamics simulations and the non-innocent role of the environment and also the possibility to open novel fragmentation pathways.
Collapse
Affiliation(s)
| | - Michael Badawi
- Université de Lorraine and CNRS, UMR 7019 LPCT, 54000, Nancy, France
| | | | - Said Ouaskit
- Laboratoire de Physique de la Matière Condensée, Faculté de Sciences Ben M'sick, University Hassan II of Casablanca, Morocco
| | - Antonio Monari
- Université de Lorraine and CNRS, UMR 7019 LPCT, 54000, Nancy, France
- Université de Paris and CNRS, ITODYS, 75006, Paris, France
| |
Collapse
|
16
|
Reliability and performances of real-time time-dependent auxiliary density functional theory. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02819-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
17
|
Femtosecond responses of hydrated DNA irradiated by ionizing rays focus on the sugar-phosphate part. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02778-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
18
|
Nottoli M, Bondanza M, Lipparini F, Mennucci B. An enhanced sampling QM/AMOEBA approach: The case of the excited state intramolecular proton transfer in solvated 3-hydroxyflavone. J Chem Phys 2021; 154:184107. [PMID: 34241028 DOI: 10.1063/5.0046844] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
We present an extension of the polarizable quantum mechanical (QM)/AMOEBA approach to enhanced sampling techniques. This is achieved by connecting the enhanced sampling PLUMED library to the machinery based on the interface of Gaussian and Tinker to perform QM/AMOEBA molecular dynamics. As an application, we study the excited state intramolecular proton transfer of 3-hydroxyflavone in two solvents: methanol and methylcyclohexane. By using a combination of molecular dynamics and umbrella sampling, we find an ultrafast component of the transfer, which is common to the two solvents, and a much slower component, which is active in the protic solvent only. The mechanisms of the two components are explained in terms of intramolecular vibrational redistribution and intermolecular hydrogen-bonding, respectively. Ground and excited state free energies along an effective reaction coordinate are finally obtained allowing for a detailed analysis of the solvent mediated mechanism.
Collapse
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
19
|
Nottoli M, Cupellini L, Lipparini F, Granucci G, Mennucci B. Multiscale Models for Light-Driven Processes. Annu Rev Phys Chem 2021; 72:489-513. [PMID: 33561359 DOI: 10.1146/annurev-physchem-090419-104031] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Multiscale models combining quantum mechanical and classical descriptions are a very popular strategy to simulate properties and processes of complex systems. Many alternative formulations have been developed, and they are now available in all of the most widely used quantum chemistry packages. Their application to the study of light-driven processes, however, is more recent, and some methodological and numerical problems have yet to be solved. This is especially the case for the polarizable formulation of these models, the recent advances in which we review here. Specifically, we identify and describe the most important specificities that the polarizable formulation introduces into both the simulation of excited-state dynamics and the modeling of excitation energy and electron transfer processes.
Collapse
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Giovanni Granucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| |
Collapse
|
20
|
Omar KA, Hasnaoui K, de la Lande A. First-Principles Simulations of Biological Molecules Subjected to Ionizing Radiation. Annu Rev Phys Chem 2021; 72:445-465. [PMID: 33878897 DOI: 10.1146/annurev-physchem-101419-013639] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Ionizing rays cause damage to genomes, proteins, and signaling pathways that normally regulate cell activity, with harmful consequences such as accelerated aging, tumors, and cancers but also with beneficial effects in the context of radiotherapies. While the great pace of research in the twentieth century led to the identification of the molecular mechanisms for chemical lesions on the building blocks of biomacromolecules, the last two decades have brought renewed questions, for example, regarding the formation of clustered damage or the rich chemistry involving the secondary electrons produced by radiolysis. Radiation chemistry is now meeting attosecond science, providing extraordinary opportunities to unravel the very first stages of biological matter radiolysis. This review provides an overview of the recent progress made in this direction, focusing mainly on the atto- to femto- to picosecond timescales. We review promising applications of time-dependent density functional theory in this context.
Collapse
Affiliation(s)
- Karwan Ali Omar
- Institut de Chimie Physique, CNRS UMR 8000, Université Paris-Saclay, 91405 Orsay, France; .,Department of Chemistry, College of Education, University of Sulaimani, 41005 Kurdistan, Iraq
| | - Karim Hasnaoui
- High Performance Computing User Support Team, Institut du Développement et des Ressources en Informatique Scientifique (IDRIS), 91403 Orsay, France.,Maison de la Simulation, CNRS, Commissariat à l'Energie Atomique et aux Énergies Alternatives (CEA), Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Aurélien de la Lande
- Institut de Chimie Physique, CNRS UMR 8000, Université Paris-Saclay, 91405 Orsay, France;
| |
Collapse
|
21
|
Santoro F, Green JA, Martinez-Fernandez L, Cerezo J, Improta R. Quantum and semiclassical dynamical studies of nonadiabatic processes in solution: achievements and perspectives. Phys Chem Chem Phys 2021; 23:8181-8199. [PMID: 33875988 DOI: 10.1039/d0cp05907b] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
We concisely review the main methodological approaches to model nonadiabatic dynamics in isotropic solutions and their applications. Three general classes of models are identified as the most used to include solvent effects in the simulations. The first model describes the solvent as a set of harmonic collective modes coupled to the solute degrees of freedom, and the second as a continuum, while the third explicitly includes solvent molecules in the calculations. The issues related to the use of these models in semiclassical and quantum dynamical simulations are discussed, as well as the main limitations and perspectives of each approach.
Collapse
Affiliation(s)
- Fabrizio Santoro
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy.
| | - James A Green
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy.
| | - Lara Martinez-Fernandez
- Departamento de Química, Facultad de Ciencias and Institute for Advanced Research in Chemistry (IADCHEM), Universidad Autónoma de Madrid, Campus de Excelencia UAM-CSIC, 28049 Madrid, Spain
| | - Javier Cerezo
- Departamento de Química, Facultad de Ciencias and Institute for Advanced Research in Chemistry (IADCHEM), Universidad Autónoma de Madrid, Campus de Excelencia UAM-CSIC, 28049 Madrid, Spain
| | - Roberto Improta
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy.
| |
Collapse
|
22
|
Pauletti CF, Coccia E, Luppi E. Role of exchange and correlation in high-harmonic generation spectra of H 2, N 2, and CO 2: Real-time time-dependent electronic-structure approaches. J Chem Phys 2021; 154:014101. [PMID: 33412879 DOI: 10.1063/5.0033072] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
This study arises from the attempt to answer the following question: how different descriptions of electronic exchange and correlation affect the high-harmonic generation (HHG) spectroscopy of H2, N2, and CO2 molecules? We compare HHG spectra for H2, N2, and CO2 with different ab initio electronic structure methods: real-time time-dependent configuration interaction and real-time time-dependent density functional theory (RT-TDDFT) using truncated basis sets composed of correlated wave functions expanded on Gaussian basis sets. In the framework of RT-TDDFT, we employ Perdew-Burke-Ernzerhof (PBE) and long-range corrected Perdew-Burke-Ernzerhof (LC-ωPBE) functionals. We study HHG spectroscopy by disentangling the effect of electronic exchange and correlation. We first analyze the electronic exchange alone, and in the case of RT-TDDFT with LC-ωPBE, we use ω = 0.3 and ω = 0.4 to tune the percentage of long-range Hartree-Fock exchange and short-range exchange PBE. Then, we added the correlation as described by the PBE functional. All the methods give very similar HHG spectra, and they seem not to be particularly sensitive to the different description of exchange and correlation or to the correct asymptotic behavior of the Coulomb potential. Despite this general trend, some differences are found in the region connecting the cutoff and the background. Here, the harmonics can be resolved with different accuracy depending on the theoretical schemes used. We believe that the investigation of the molecular continuum and its coupling with strong fields merits further theoretical investigations in the near future.
Collapse
Affiliation(s)
| | - Emanuele Coccia
- Dipartimento di Scienze Chimiche e Farmaceutiche, Via Giorgieri 1, Trieste Italy
| | - Eleonora Luppi
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, F-75005 Paris, France
| |
Collapse
|
23
|
Nottoli M, Lipparini F. General formulation of polarizable embedding models and of their coupling. J Chem Phys 2020; 153:224108. [PMID: 33317291 DOI: 10.1063/5.0035165] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Univeristà di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Univeristà di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
24
|
Müller C, Sharma M, Sierka M. Real-time time-dependent density functional theory using density fitting and the continuous fast multipole method. J Comput Chem 2020; 41:2573-2582. [PMID: 33464600 DOI: 10.1002/jcc.26412] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Revised: 08/17/2020] [Accepted: 08/18/2020] [Indexed: 11/07/2022]
Abstract
An implementation of real-time time-dependent density functional theory (RT-TDDFT) within the TURBOMOLE program package is reported using Gaussian-type orbitals as basis functions, second and fourth order Magnus propagator, and the self-consistent field as well as the predictor-corrector time integration schemes. The Coulomb contribution to the Kohn-Sham matrix is calculated combining density fitting approximation and the continuous fast multipole method. Performance of the implementation is benchmarked for molecular systems with different sizes and dimensionalities. For linear alkane chains, the wall time for density matrix time propagation step is comparable to the Kohn-Sham (KS) matrix construction. However, for larger two- and three-dimensional molecules, with up to about 5,000 basis functions, the computational effort of RT-TDDFT calculations is dominated by the KS matrix evaluation. In addition, the maximum time step is evaluated using a set of small molecules of different polarities. The photoabsorption spectra of several molecular systems calculated using RT-TDDFT are compared to those obtained using linear response time-dependent density functional theory and coupled cluster methods.
Collapse
Affiliation(s)
- Carolin Müller
- Otto Schott Institute of Materials Research, Friedrich Schiller University of Jena, Jena, Germany
| | - Manas Sharma
- Otto Schott Institute of Materials Research, Friedrich Schiller University of Jena, Jena, Germany
| | - Marek Sierka
- Otto Schott Institute of Materials Research, Friedrich Schiller University of Jena, Jena, Germany
| |
Collapse
|
25
|
Nottoli M, Mennucci B, Lipparini F. Excited state Born-Oppenheimer molecular dynamics through coupling between time dependent DFT and AMOEBA. Phys Chem Chem Phys 2020; 22:19532-19541. [PMID: 32844823 DOI: 10.1039/d0cp03688a] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We present the implementation of excited state Born-Oppenheimer molecular dynamics (BOMD) using a polarizable QM/MM approach based on a time-dependent density functional theory (TDDFT) formulation and the AMOEBA force field. The implementation relies on an interface between Tinker and Gaussian software and it uses an algorithm for the calculation of QM/MM energy and forces which scales linearly with the number of MM atoms. The resulting code can perform TDDFT/AMOEBA BOMD simulations on real-life systems with standard computational resources. As a test case, the method is applied to the study of the mechanism of locally-excited to charge-transfer conversion in dimethylaminobenzonitrile in a polar solvent. Our simulations confirm that such a conversion is governed by the twisting of the dimethylamino group which is accompanied by an important reorientation of solvent molecules.
Collapse
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy.
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy.
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy.
| |
Collapse
|
26
|
De Santis M, Belpassi L, Jacob CR, Severo Pereira Gomes A, Tarantelli F, Visscher L, Storchi L. Environmental Effects with Frozen-Density Embedding in Real-Time Time-Dependent Density Functional Theory Using Localized Basis Functions. J Chem Theory Comput 2020; 16:5695-5711. [PMID: 32786918 PMCID: PMC8009524 DOI: 10.1021/acs.jctc.0c00603] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Indexed: 12/14/2022]
Abstract
Frozen-density embedding (FDE) represents a versatile embedding scheme to describe the environmental effect on electron dynamics in molecular systems. The extension of the general theory of FDE to the real-time time-dependent Kohn-Sham method has previously been presented and implemented in plane waves and periodic boundary conditions [Pavanello, M.; J. Chem. Phys. 2015, 142, 154116]. In the current paper, we extend our recent formulation of the real-time time-dependent Kohn-Sham method based on localized basis set functions and developed within the Psi4NumPy framework to the FDE scheme. The latter has been implemented in its "uncoupled" flavor (in which the time evolution is only carried out for the active subsystem, while the environment subsystems remain at their ground state), using and adapting the FDE implementation already available in the PyEmbed module of the scripting framework PyADF. The implementation was facilitated by the fact that both Psi4NumPy and PyADF, being native Python API, provided an ideal framework of development using the Python advantages in terms of code readability and reusability. We employed this new implementation to investigate the stability of the time-propagation procedure, which is based on an efficient predictor/corrector second-order midpoint Magnus propagator employing an exact diagonalization, in combination with the FDE scheme. We demonstrate that the inclusion of the FDE potential does not introduce any numerical instability in time propagation of the density matrix of the active subsystem, and in the limit of the weak external field, the numerical results for low-lying transition energies are consistent with those obtained using the reference FDE calculations based on the linear-response TDDFT. The method is found to give stable numerical results also in the presence of a strong external field inducing nonlinear effects. Preliminary results are reported for high harmonic generation (HHG) of a water molecule embedded in a small water cluster. The effect of the embedding potential is evident in the HHG spectrum reducing the number of the well-resolved high harmonics at high energy with respect to the free water. This is consistent with a shift toward lower ionization energy passing from an isolated water molecule to a small water cluster. The computational burden for the propagation step increases approximately linearly with the size of the surrounding frozen environment. Furthermore, we have also shown that the updating frequency of the embedding potential may be significantly reduced, much less than one per time step, without jeopardizing the accuracy of the transition energies.
Collapse
Affiliation(s)
- Matteo De Santis
- Dipartimento di
Chimica, Biologia e Biotecnologie, Università
degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Istituto di Scienze
e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche
c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Leonardo Belpassi
- Istituto di Scienze
e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche
c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Christoph R. Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstr. 17, 38106 Braunschweig, Germany
| | | | - Francesco Tarantelli
- Dipartimento di
Chimica, Biologia e Biotecnologie, Università
degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Lucas Visscher
- Theoretical Chemistry, Faculty of Science, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Loriano Storchi
- Istituto di Scienze
e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche
c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Dipartimento di Farmacia, Università
degli Studi ‘G. D’Annunzio’, Via dei Vestini 31, 66100 Chieti, Italy
| |
Collapse
|
27
|
Li X, Govind N, Isborn C, DePrince AE, Lopata K. Real-Time Time-Dependent Electronic Structure Theory. Chem Rev 2020; 120:9951-9993. [DOI: 10.1021/acs.chemrev.0c00223] [Citation(s) in RCA: 64] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christine Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
28
|
De Santis M, Storchi L, Belpassi L, Quiney HM, Tarantelli F. PyBERTHART: A Relativistic Real-Time Four-Component TDDFT Implementation Using Prototyping Techniques Based on Python. J Chem Theory Comput 2020; 16:2410-2429. [PMID: 32101419 DOI: 10.1021/acs.jctc.0c00053] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Matteo De Santis
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Istituto di Scienze e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Loriano Storchi
- Istituto di Scienze e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Dipartimento di Farmacia, Università degli Studi ‘G. D’Annunzio’, Via dei Vestini 31, 66100 Chieti, Italy
| | - Leonardo Belpassi
- Istituto di Scienze e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Harry M. Quiney
- ARC Centre of Excellence for Advanced Molecular Imaging, School of Physics, The University of Melbourne, 3010 Victoria, Australia
| | - Francesco Tarantelli
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Istituto di Scienze e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| |
Collapse
|
29
|
Lingerfelt DB, Ganesh P, Jakowski J, Sumpter BG. Understanding Beam-Induced Electronic Excitations in Materials. J Chem Theory Comput 2020; 16:1200-1214. [DOI: 10.1021/acs.jctc.9b00792] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David B. Lingerfelt
- Nanomaterials Theory Institute, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Panchapakesan Ganesh
- Nanomaterials Theory Institute, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Jacek Jakowski
- Nanomaterials Theory Institute, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Bobby G. Sumpter
- Nanomaterials Theory Institute, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| |
Collapse
|
30
|
Bondanza M, Nottoli M, Cupellini L, Lipparini F, Mennucci B. Polarizable embedding QM/MM: the future gold standard for complex (bio)systems? Phys Chem Chem Phys 2020; 22:14433-14448. [DOI: 10.1039/d0cp02119a] [Citation(s) in RCA: 59] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We provide a perspective of the induced dipole formulation of polarizable QM/MM, showing how efficient implementations will enable their application to the modeling of dynamics, spectroscopy, and reactivity in complex biosystems.
Collapse
Affiliation(s)
- Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| |
Collapse
|
31
|
Alvarez-Ibarra A, Parise A, Hasnaoui K, de la Lande A. The physical stage of radiolysis of solvated DNA by high-energy-transfer particles: insights from new first principles simulations. Phys Chem Chem Phys 2020; 22:7747-7758. [DOI: 10.1039/d0cp00165a] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Electron dynamics simulations based on density functional theory are carried out on nanometric molecular systems to decipher the primary processes following irradiation of bio-macromolecules by high energy transfer charged particles.
Collapse
Affiliation(s)
| | - Angela Parise
- Université Paris-Saclay
- CNRS
- Institut de Chimie Physique UMR8000
- Orsay
- France
| | - Karim Hasnaoui
- Institut du Développement et des Ressources en Informatique Scientifique
- Rue John von Neumann
- Orsay
- France
- Maison de la Simulation
| | | |
Collapse
|
32
|
de la Lande A, Alvarez-Ibarra A, Hasnaoui K, Cailliez F, Wu X, Mineva T, Cuny J, Calaminici P, López-Sosa L, Geudtner G, Navizet I, Garcia Iriepa C, Salahub DR, Köster AM. Molecular Simulations with in-deMon2k QM/MM, a Tutorial-Review. Molecules 2019; 24:molecules24091653. [PMID: 31035516 PMCID: PMC6539060 DOI: 10.3390/molecules24091653] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2019] [Revised: 04/11/2019] [Accepted: 04/12/2019] [Indexed: 12/18/2022] Open
Abstract
deMon2k is a readily available program specialized in Density Functional Theory (DFT) simulations within the framework of Auxiliary DFT. This article is intended as a tutorial-review of the capabilities of the program for molecular simulations involving ground and excited electronic states. The program implements an additive QM/MM (quantum mechanics/molecular mechanics) module relying either on non-polarizable or polarizable force fields. QM/MM methodologies available in deMon2k include ground-state geometry optimizations, ground-state Born-Oppenheimer molecular dynamics simulations, Ehrenfest non-adiabatic molecular dynamics simulations, and attosecond electron dynamics. In addition several electric and magnetic properties can be computed with QM/MM. We review the framework implemented in the program, including the most recently implemented options (link atoms, implicit continuum for remote environments, metadynamics, etc.), together with six applicative examples. The applications involve (i) a reactivity study of a cyclic organic molecule in water; (ii) the establishment of free-energy profiles for nucleophilic-substitution reactions by the umbrella sampling method; (iii) the construction of two-dimensional free energy maps by metadynamics simulations; (iv) the simulation of UV-visible absorption spectra of a solvated chromophore molecule; (v) the simulation of a free energy profile for an electron transfer reaction within Marcus theory; and (vi) the simulation of fragmentation of a peptide after collision with a high-energy proton.
Collapse
Affiliation(s)
- Aurélien de la Lande
- Laboratoire de Chimie Physique, CNRS, Université Paris Sud, Université Paris Saclay, 15 avenue Jean Perrin, 91405 Orsay, France.
| | - Aurelio Alvarez-Ibarra
- Laboratoire de Chimie Physique, CNRS, Université Paris Sud, Université Paris Saclay, 15 avenue Jean Perrin, 91405 Orsay, France.
| | - Karim Hasnaoui
- Laboratoire de Chimie Physique, CNRS, Université Paris Sud, Université Paris Saclay, 15 avenue Jean Perrin, 91405 Orsay, France.
| | - Fabien Cailliez
- Laboratoire de Chimie Physique, CNRS, Université Paris Sud, Université Paris Saclay, 15 avenue Jean Perrin, 91405 Orsay, France.
| | - Xiaojing Wu
- Laboratoire de Chimie Physique, CNRS, Université Paris Sud, Université Paris Saclay, 15 avenue Jean Perrin, 91405 Orsay, France.
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique, PSL University, 75005 Paris, France.
| | - Tzonka Mineva
- Matériaux Avancés pour la Catalyse et la Santé, UMR 5253 CNRS/UM/ENSCM, Institut Charles Gerhardt de Montpellier (ICGM) Montpellier CEDEX 5, 34090 Montpellier, France.
| | - Jérôme Cuny
- Laboratoire de Chimie et Physique Quantiques, IRSAMC, Université Paul Sabatier, 118 Route de Narbonne, 31062 Toulouse CEDEX 4, France.
| | - Patrizia Calaminici
- Programa de Doctorado en Nanociencias y Nanotecnología, CINVESTAV, Av. Instituto Politécnico Nacional, 2508, A.P. 14-740, Ciudad de México 07000, Mexico.
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional, 2508, A.P. 14-740, Ciudad de México 07000, México.
| | - Luis López-Sosa
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional, 2508, A.P. 14-740, Ciudad de México 07000, México.
| | - Gerald Geudtner
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional, 2508, A.P. 14-740, Ciudad de México 07000, México.
| | - Isabelle Navizet
- Laboratoire Modélisation et Simulation Multi Échelle, Université Paris-Est, MSME, UMR 8208 CNRS, UPEM, 5 bd Descartes, 77454 Marne-la-Vallée, France.
| | - Cristina Garcia Iriepa
- Laboratoire Modélisation et Simulation Multi Échelle, Université Paris-Est, MSME, UMR 8208 CNRS, UPEM, 5 bd Descartes, 77454 Marne-la-Vallée, France.
| | - Dennis R Salahub
- Department of Chemistry, Centre for Molecular Simulation, Institute for Quantum Science and Technology and Quantum Alberta, University of Calgary, 2500 University Drive N.W., Calgary, AB T2N 1N4, Canada.
- College of Chemistry and Chemical Engineering, Henan University of Technology, No. 100, Lian Hua Street, High-Tech Development Zone, Zhengzhou 450001, China.
| | - Andreas M Köster
- Programa de Doctorado en Nanociencias y Nanotecnología, CINVESTAV, Av. Instituto Politécnico Nacional, 2508, A.P. 14-740, Ciudad de México 07000, Mexico.
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional, 2508, A.P. 14-740, Ciudad de México 07000, México.
| |
Collapse
|
33
|
Wildman A, Donati G, Lipparini F, Mennucci B, Li X. Nonequilibrium Environment Dynamics in a Frequency-Dependent Polarizable Embedding Model. J Chem Theory Comput 2018; 15:43-51. [PMID: 30512961 DOI: 10.1021/acs.jctc.8b00836] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Hybrid quantum mechanical/molecular mechanical (QM/MM) models are some of the most powerful and computationally feasible approaches to account for solvent effects or more general environmental perturbations on quantum chemical systems. In their more recent formulations (known as polarizable embedding) they can account for electrostatic and mutual polarization effects between the QM and the MM subsystems. In this paper, a polarizable embedding scheme based on induced dipoles that is able both to describe electron evolution of the embedded QM system in an efficient manner as well as to capture the frequency dependent behavior of the solvent is proposed, namely, ωMMPol. The effects of this frequency-dependent solvent on a time-dependent model system-the Rabi oscillations of H2+ in a resonant field-are considered. The solvent is shown to introduce only mild perturbations when the excitation frequencies of the solvent and the solute are off-resonant. However, the dynamics of the H2+ are fundamentally changed in the presence of a near-resonant excitation solvent. The effectiveness of ωMMPol to simulating realistic chemical systems is demonstrated by capturing charge transfer dynamics within a solvated system.
Collapse
Affiliation(s)
- Andrew Wildman
- Department of Chemistry , University of Washington , Seattle , Washington 98195 , United States
| | - Greta Donati
- Department of Chemistry , University of Washington , Seattle , Washington 98195 , United States
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale , Universitá di Pisa , Via Risorgimento 35 , 56126 Pisa , Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale , Universitá di Pisa , Via Risorgimento 35 , 56126 Pisa , Italy
| | - Xiaosong Li
- Department of Chemistry , University of Washington , Seattle , Washington 98195 , United States
| |
Collapse
|
34
|
Morzan UN, Alonso de Armiño DJ, Foglia NO, Ramírez F, González Lebrero MC, Scherlis DA, Estrin DA. Spectroscopy in Complex Environments from QM–MM Simulations. Chem Rev 2018; 118:4071-4113. [DOI: 10.1021/acs.chemrev.8b00026] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Uriel N. Morzan
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Diego J. Alonso de Armiño
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Nicolás O. Foglia
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Francisco Ramírez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Mariano C. González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Damián A. Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Darío A. Estrin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| |
Collapse
|
35
|
Parise A, Alvarez-Ibarra A, Wu X, Zhao X, Pilmé J, Lande ADL. Quantum Chemical Topology of the Electron Localization Function in the Field of Attosecond Electron Dynamics. J Phys Chem Lett 2018; 9:844-850. [PMID: 29384381 DOI: 10.1021/acs.jpclett.7b03379] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We report original analyses of attosecond electron dynamics of molecules subject to collisions by high energy charged particles based on Real-Time Time-Dependent-Density-Functional-Theory simulations coupled to Topological Analyses of the Electron Localization Function (TA-TD-ELF). We investigate irradiation of water and guanine. TA-TD-ELF enables qualitative and quantitative characterizations of bond breaking and formation, of charge migration within topological basins, or of electron attachment to the colliding particle. Whereas the Lewis-VSEPR structure of gas phase water is blown out within a few attoseconds after collision, that of guanine is far more robust and reconstitutes rapidly after impact even though the molecule remains electronically excited. This difference is accounted by the presence of the electron bath surrounding the impact point which enables energy relaxation within the molecule. Our approach should stimulate future studies to unravel the early steps following irradiation of various types of systems (isolated molecules, biomolecules, nanoclusters, solids, etc.) and is also readily applicable to irradiation by photons of various energies.
Collapse
Affiliation(s)
- Angela Parise
- Laboratoire de Chimie Physique, Université Paris Sud, CNRS , Université Paris Saclay. 15 avenue Jean Perrin, F91405 Orsay, France
| | - Aurelio Alvarez-Ibarra
- Laboratoire de Chimie Physique, Université Paris Sud, CNRS , Université Paris Saclay. 15 avenue Jean Perrin, F91405 Orsay, France
| | - Xiaojing Wu
- Laboratoire de Chimie Physique, Université Paris Sud, CNRS , Université Paris Saclay. 15 avenue Jean Perrin, F91405 Orsay, France
| | - Xiaodong Zhao
- Laboratoire de Chimie Physique, Université Paris Sud, CNRS , Université Paris Saclay. 15 avenue Jean Perrin, F91405 Orsay, France
| | - Julien Pilmé
- Laboratoire de Chimie Théorique, Sorbonne Universités, Université Pierre et Marie Curie, CNRS , F75005 Paris, France
| | - Aurélien de la Lande
- Laboratoire de Chimie Physique, Université Paris Sud, CNRS , Université Paris Saclay. 15 avenue Jean Perrin, F91405 Orsay, France
| |
Collapse
|
36
|
Donati G, Wildman A, Caprasecca S, Lingerfelt DB, Lipparini F, Mennucci B, Li X. Coupling Real-Time Time-Dependent Density Functional Theory with Polarizable Force Field. J Phys Chem Lett 2017; 8:5283-5289. [PMID: 28994290 DOI: 10.1021/acs.jpclett.7b02320] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Real-time time-dependent density functional theory (RT-TDDFT) is a powerful tool for obtaining spectroscopic observables and understanding complex, time-dependent properties. Currently, performing RT-TDDFT calculations on large, fully quantum mechanical systems is not computationally feasible. Previously, polarizable mixed quantum mechanical and molecular mechanical (QM/MMPol) models have been successful in providing accurate, yet efficient, approximations to a fully quantum mechanical system. Here we develop a coupling scheme between induced dipole based QM/MMPol and RT-TDDFT. Our approach is validated by comparing calculated spectra with both real-time and linear-response TDDFT calculations. The model developed within provides an accurate method for performing RT-TDDFT calculations on extended systems while accounting for mutual polarization between the quantum mechanical and molecular mechanical regions.
Collapse
Affiliation(s)
- Greta Donati
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - Andrew Wildman
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - Stefano Caprasecca
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa , Via G. Moruzzi 13, 56124 Pisa, Italy
| | - David B Lingerfelt
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa , Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa , Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Xiaosong Li
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| |
Collapse
|