1
|
Romanelli M, Corni S. Identifying Differences between Semiclassical and Full-Quantum Descriptions of Plexcitons. J Phys Chem Lett 2024:9326-9334. [PMID: 39236151 DOI: 10.1021/acs.jpclett.4c01519] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2024]
Abstract
Strong light-matter coupling between molecules and plasmonic nanoparticles gives rise to new hybrid eigenstates of the coupled system, commonly referred to as polaritons or, more precisely, plexcitons. Over the past decade, it has been amply shown that molecular electron dynamics and photophysics can be drastically affected by such interactions, thus paving the way for light-induced control of molecular excited state properties and reactivity. Here, by combining the ab initio molecular description and classical or quantum modeling of arbitrarily shaped plasmonic nanostructures within the stochastic Schrödinger equation, we present two approaches, one semiclassical and one full-quantum, to follow in real time the electronic dynamics of plexcitons while realistically taking plasmonic dissipative losses into account. The full-quantum theory is compared with the semiclassical analogue under different interaction regimes, showing (numerically and theoretically) that even in the weak-field and weak-coupling limit a small-yet-observable difference arises.
Collapse
Affiliation(s)
- Marco Romanelli
- Department of Chemical Sciences, University of Padova, via Marzolo 1, 35131 Padova, Italy
| | - Stefano Corni
- Department of Chemical Sciences, University of Padova, via Marzolo 1, 35131 Padova, Italy
- CNR Institute of Nanoscience, via Campi 213/A, 41125 Modena, Italy
- Padua Quantum Technologies Research Center, University of Padova, 35131 Padova, Italy
| |
Collapse
|
2
|
Capone M, Romanelli M, Castaldo D, Parolin G, Bello A, Gil G, Vanzan M. A Vision for the Future of Multiscale Modeling. ACS PHYSICAL CHEMISTRY AU 2024; 4:202-225. [PMID: 38800726 PMCID: PMC11117712 DOI: 10.1021/acsphyschemau.3c00080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/30/2023] [Revised: 01/31/2024] [Accepted: 02/01/2024] [Indexed: 05/29/2024]
Abstract
The rise of modern computer science enabled physical chemistry to make enormous progresses in understanding and harnessing natural and artificial phenomena. Nevertheless, despite the advances achieved over past decades, computational resources are still insufficient to thoroughly simulate extended systems from first principles. Indeed, countless biological, catalytic and photophysical processes require ab initio treatments to be properly described, but the breadth of length and time scales involved makes it practically unfeasible. A way to address these issues is to couple theories and algorithms working at different scales by dividing the system into domains treated at different levels of approximation, ranging from quantum mechanics to classical molecular dynamics, even including continuum electrodynamics. This approach is known as multiscale modeling and its use over the past 60 years has led to remarkable results. Considering the rapid advances in theory, algorithm design, and computing power, we believe multiscale modeling will massively grow into a dominant research methodology in the forthcoming years. Hereby we describe the main approaches developed within its realm, highlighting their achievements and current drawbacks, eventually proposing a plausible direction for future developments considering also the emergence of new computational techniques such as machine learning and quantum computing. We then discuss how advanced multiscale modeling methods could be exploited to address critical scientific challenges, focusing on the simulation of complex light-harvesting processes, such as natural photosynthesis. While doing so, we suggest a cutting-edge computational paradigm consisting in performing simultaneous multiscale calculations on a system allowing the various domains, treated with appropriate accuracy, to move and extend while they properly interact with each other. Although this vision is very ambitious, we believe the quick development of computer science will lead to both massive improvements and widespread use of these techniques, resulting in enormous progresses in physical chemistry and, eventually, in our society.
Collapse
Affiliation(s)
- Matteo Capone
- Department
of Physical and Chemical Sciences, University
of L’Aquila, L’Aquila 67010, Italy
| | - Marco Romanelli
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
| | - Davide Castaldo
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
| | - Giovanni Parolin
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
| | - Alessandro Bello
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
- Department
of Physics, Informatics and Mathematics, University of Modena and Reggio Emilia, Modena 41125, Italy
| | - Gabriel Gil
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
- Instituto
de Cibernética, Matemática y Física (ICIMAF), La Habana 10400, Cuba
| | - Mirko Vanzan
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
- Department
of Physics, University of Milano, Milano 20133, Italy
| |
Collapse
|
3
|
Yamada A. Molecular Simulation Study of Surface-Enhanced Raman Scattering of Liquid Water. J Phys Chem A 2024; 128:51-61. [PMID: 38127813 DOI: 10.1021/acs.jpca.3c05027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
We developed in our previous study [J. Chem. Phys., 2021, 155, 174118, J. Phys. Chem. A, 2022, 126, 4762] a classical electronic and molecular dynamics simulation method to describe the optical response of metal material in solution based on an atomistic model by incorporating the classical equation of motion for free electrons under an applied electric field. To show further usefulness of the method, in the present study, we apply it to surface-enhanced Raman scattering of liquid water to examine the signal enhancement of the solution system caused by plasmon resonance effects of a silver nanoparticle.
Collapse
Affiliation(s)
- Atsushi Yamada
- Department of Applied Chemistry, National Defense Academy, 1-10-20 Hashirimizu, Yokosuka, Kanagawa 239-8686, Japan
| |
Collapse
|
4
|
Dall’Osto G, Corni S. Time Resolved Raman Scattering of Molecules: A Quantum Mechanics Approach with Stochastic Schroedinger Equation. J Phys Chem A 2022; 126:8088-8100. [PMID: 36278928 PMCID: PMC9639147 DOI: 10.1021/acs.jpca.2c05245] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Raman scattering is a very powerful tool employed to
characterize
molecular systems. Here we propose a novel theoretical strategy to
calculate the Raman cross-section in time domain, by computing the
cumulative Raman signal emitted during the molecular evolution in
time. Our model is based on a numerical propagation of the vibronic
wave function under the effect of a light pulse of arbitrary shape.
This approach can therefore tackle a variety of experimental setups.
Both resonance and nonresonance Raman scattering can be retrieved,
and also the time-dependent fluorescence emission is computed. The
model has been applied to porphyrin considering both resonance and
nonresonance conditions and varying the incident field duration. Moreover
the effect of the vibrational relaxation, which should be taken into
account when its time scale is similar to that of the Raman emission,
has been included through the stochastic Schroedinger equation approach.
Collapse
Affiliation(s)
- Giulia Dall’Osto
- Department of Chemical Sciences, University of Padova, via Marzolo 1, Padova, 35131, Italy
| | - Stefano Corni
- Department of Chemical Sciences, University of Padova, via Marzolo 1, Padova, 35131, Italy
- CNR Institute of Nanoscience, via Campi 213/A, Modena, 41125, Italy
| |
Collapse
|
5
|
Yamada A. Computational Analyses of Plasmonics of a Silver Nanoparticle in a Vacuum and in a Water Solution by Classical Electronic and Molecular Dynamics Simulations. J Phys Chem A 2022; 126:4762-4771. [PMID: 35853175 DOI: 10.1021/acs.jpca.2c02811] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present basic optical responses of a silver nanoparticle (Ag309) in a vacuum and a water solution obtained by classical electronic and molecular dynamics (CEMD) calculations, where the CEMD is our previously developed force-field based molecular dynamics simulation method that incorporates the classical equation of motion for free electrons in metal and an interaction with the applied oscillating electric field (Yamada, A. J. Chem. Phys. 2021, 155, 174118)). The present work is the follow-up of the previous study. Calculated absorption spectra in the visible region with various conditions are reported together with parameter determination for realistic analyses as well as the verification of the method. Time-domain optical responses of Ag309 in water solvent are as well analyzed in terms of the plasmon resonance excitation under a subpicosecond light pulse and its thermal relaxation.
Collapse
Affiliation(s)
- Atsushi Yamada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| |
Collapse
|
6
|
Yamada A. Classical electronic and molecular dynamics simulation for optical response of metal system. J Chem Phys 2021; 155:174118. [PMID: 34742192 DOI: 10.1063/5.0067144] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
An extended molecular dynamics simulation that incorporates classical free electron dynamics in the framework of the force-field model has been developed to enable us to describe the optical response of metal materials under the visible light electric field. In the simulation, dynamical atomic point charges follow equations of motion of classical free electrons that include Coulomb interactions with the oscillating field and surrounding atomic sites and collision effects from nearby electrons and ions. This scheme allows us to simulate an interacting system of metals with molecules using an ordinary polarizable force-field and preserves energy conservation in the case without applying an external electric field. As the first applications, we show that the presented simulation accurately reproduces (i) the classical image potential in a metal-charge interaction system and (ii) the dielectric function of bulk metal. We also demonstrate (iii) calculations of absorption spectra of metal nano-particles with and without a water solvent at room temperature, showing reasonable red-shift by the solvent effect, and (iv) plasmon resonant excitation of the metal nano-particle in solution under the visible light pulse and succeeding energy relaxation of the absorbed light energy from electrons to atoms on the metal and to the water solvent. Our attempt thus opens the possibility to expand the force-field based molecular dynamics simulation to an alternative tool for optical-related fields.
Collapse
Affiliation(s)
- Atsushi Yamada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| |
Collapse
|
7
|
Grobas Illobre P, Marsili M, Corni S, Stener M, Toffoli D, Coccia E. Time-Resolved Excited-State Analysis of Molecular Electron Dynamics by TDDFT and Bethe-Salpeter Equation Formalisms. J Chem Theory Comput 2021; 17:6314-6329. [PMID: 34486881 PMCID: PMC8515806 DOI: 10.1021/acs.jctc.1c00211] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2021] [Indexed: 12/16/2022]
Abstract
In this work, a theoretical and computational set of tools to study and analyze time-resolved electron dynamics in molecules, under the influence of one or more external pulses, is presented. By coupling electronic-structure methods with the resolution of the time-dependent Schrödinger equation, we developed and implemented the time-resolved induced density of the electronic wavepacket, the time-resolved formulation of the differential projection density of states (ΔPDOS), and of transition contribution map (TCM) to look at the single-electron orbital occupation and localization change in time. Moreover, to further quantify the possible charge transfer, we also defined the energy-integrated ΔPDOS and the fragment-projected TCM. We have used time-dependent density-functional theory (TDDFT), as implemented in ADF software, and the Bethe-Salpeter equation, as provided by MolGW package, for the description of the electronic excited states. This suite of postprocessing tools also provides the time evolution of the electronic states of the system of interest. To illustrate the usefulness of these postprocessing tools, excited-state populations have been computed for HBDI (the chromophore of GFP) and DNQDI molecules interacting with a sequence of two pulses. Time-resolved descriptors have been applied to study the time-resolved electron dynamics of HBDI, DNQDI, LiCN (being a model system for dipole switching upon highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) electronic excitation), and Ag22. The computational analysis tools presented in this article can be employed to help the interpretation of fast and ultrafast spectroscopies on molecular, supramolecular, and composite systems.
Collapse
Affiliation(s)
- P. Grobas Illobre
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| | - M. Marsili
- Dipartimento
di Scienze Chimiche, Universitá di
Padova, via Marzolo 1, Padova 35131, Italy
| | - S. Corni
- Dipartimento
di Scienze Chimiche, Universitá di
Padova, via Marzolo 1, Padova 35131, Italy
- CNR
Istituto di Nanoscienze, via Campi 213/A, Modena 41125, Italy
| | - M. Stener
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| | - D. Toffoli
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| | - E. Coccia
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| |
Collapse
|
8
|
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: 28] [Impact Index Per Article: 9.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
|
9
|
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
|
10
|
Rosa M, Gil G, Corni S, Cammi R. Quantum optimal control theory for solvated systems. J Chem Phys 2019; 151:194109. [PMID: 31757146 DOI: 10.1063/1.5125184] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
In this work, we generalize the quantum optimal control theory (QOCT) of molecules subject to ultrashort laser pulses to the case of solvated systems, explicitly including the solvent dielectric properties in the system's quantum Hamiltonian. A reliable description of the solvent polarization is accounted for within the polarizable continuum model (PCM). The electron dynamics for the molecules in solution is coupled with the dynamics of the surrounding polarizable environment, which affects the features of the optimized laser pulse. To illustrate such effects, numerical applications of the developed method to the study of optimal population of selected excited states of two molecular solvated systems are presented and discussed.
Collapse
Affiliation(s)
- Marta Rosa
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Padova, Italy
| | - Gabriel Gil
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Padova, Italy
| | - Stefano Corni
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Padova, Italy
| | - Roberto Cammi
- Dipartimento di Scienze Chimiche, della Vita e della Sostenibilità Ambientale, Università di Parma, Parma, Italy
| |
Collapse
|
11
|
Coccia E, Troiani F, Corni S. Probing quantum coherence in ultrafast molecular processes: Anab initioapproach to open quantum systems. J Chem Phys 2018; 148:204112. [DOI: 10.1063/1.5022976] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Emanuele Coccia
- Department of Chemical Sciences, University of Padova, via Marzolo 1, Padova, Italy
- CNR Institute of Nanoscience, via Campi 213/A, Modena, Italy
| | - Filippo Troiani
- CNR Institute of Nanoscience, via Campi 213/A, Modena, Italy
| | - Stefano Corni
- Department of Chemical Sciences, University of Padova, via Marzolo 1, Padova, Italy
- CNR Institute of Nanoscience, via Campi 213/A, Modena, Italy
| |
Collapse
|
12
|
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
|
13
|
Goings JJ, Lestrange PJ, Li X. Real‐time time‐dependent electronic structure theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1341] [Citation(s) in RCA: 84] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
| | | | - Xiaosong Li
- Department of ChemistryUniversity of Washington Seattle WA USA
| |
Collapse
|