1
|
García-González F, Otero JC, Ávila Ferrer FJ, Santoro F, Aranda D. Linear Vibronic Coupling Approach for Surface-Enhanced Raman Scattering: Quantifying the Charge-Transfer Enhancement Mechanism. J Chem Theory Comput 2024; 20:3850-3863. [PMID: 38687961 PMCID: PMC11099975 DOI: 10.1021/acs.jctc.4c00061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Revised: 04/12/2024] [Accepted: 04/15/2024] [Indexed: 05/02/2024]
Abstract
The outstanding amplification observed in surface-enhanced Raman scattering (SERS) is due to several enhancement mechanisms, and standing out among them are the plasmonic (PL) and charge-transfer (CT) mechanisms. The theoretical estimation of the enhancement factors of the CT mechanism is challenging because the excited-state coupling between bright plasmons and dark CT states must be properly introduced into the model to obtain reliable intensities. In this work, we aim at simulating electrochemical SERS spectra, considering models of pyridine on silver clusters subjected to an external electric field E⃗ that represents the effect of an electrode potential Vel. The method adopts quantum dynamical propagations of nuclear wavepackets on the coupled PL and CT states described with linear vibronic coupling models parametrized for each E⃗ through a fragment-based maximum-overlap diabatization. By presenting results at different values of E⃗, we show that indeed there is a relation between the population transferred to the CT states and the total scattered intensity. The tuning and detuning processes of the CT states with the bright PLs as a function of the electric field are in good agreement with those observed in experiments. Finally, our estimations for the CT enhancement factors predict values in the order of 105 to 106, meaning that when the CT and PL states are both in resonance with the excitation wavelength, the CT and PL enhancements are comparable, and vibrational bands whose intensity is amplified by different mechanisms can be observed together, in agreement with what was measured by typical experiments on silver electrodes.
Collapse
Affiliation(s)
- Francisco García-González
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| | - Juan Carlos Otero
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| | - Francisco J. Ávila Ferrer
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| | - Fabrizio Santoro
- Istituto
di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca
del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Daniel Aranda
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| |
Collapse
|
2
|
Malpathak S, Church MS, Ananth N. A Semiclassical Framework for Mixed Quantum Classical Dynamics. J Phys Chem A 2022; 126:6359-6375. [PMID: 36070472 DOI: 10.1021/acs.jpca.2c03467] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Semiclassical (SC) approximations for quantum dynamic simulations in complex chemical systems range from rigorously accurate methods that are computationally expensive to methods that exhibit near-classical scaling with system size but are limited in their ability to describe quantum effects. In practical studies of high-dimensional reactions, neither extreme is the best choice: frequently a high-level quantum mechanical description is only required for a handful of modes, while the majority of environment modes that do not play a key role in the reactive event of interest are well served with a lower level of theory. In this feature, we introduce modified Filinov filtration as a powerful tool to construct mixed quantum-classical SC theories where different subsystems can be quantized to different extents without introducing ad hoc intersubsystem interaction terms. We demonstrate that these Filinov-based SC methods can systematically tune between quantum and classical limit SC behavior, offering a practical way forward to accurate and computationally efficient simulations of high-dimensional quantum processes.
Collapse
Affiliation(s)
- Shreyas Malpathak
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University Ithaca, New York 14853, United States
| | - Matthew S Church
- Department of Chemistry, Brown University, Providence, Rhode Island 02906, United States
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University Ithaca, New York 14853, United States
| |
Collapse
|
3
|
McKeon CA, Hamed SM, Bruneval F, Neaton JB. An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules. J Chem Phys 2022; 157:074103. [DOI: 10.1063/5.0097582] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.
Collapse
Affiliation(s)
- Caroline A. McKeon
- Department of Physics, University of California, Berkeley, California 94720, USA
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Natural Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Samia M. Hamed
- Department of Physics, University of California, Berkeley, California 94720, USA
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Natural Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Fabien Bruneval
- Université Paris-Saclay, CEA, Service de Recherches de Métallurgie Physique, 91191 Gif-sur-Yvette, France
| | - Jeffrey B. Neaton
- Department of Physics, University of California, Berkeley, California 94720, USA
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Natural Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Kavli ENSI, University of California, Berkeley, California 94720, USA
| |
Collapse
|
4
|
Segalina A, Aranda D, Green JA, Cristino V, Caramori S, Prampolini G, Pastore M, Santoro F. How the Interplay among Conformational Disorder, Solvation, Local, and Charge-Transfer Excitations Affects the Absorption Spectrum and Photoinduced Dynamics of Perylene Diimide Dimers: A Molecular Dynamics/Quantum Vibronic Approach. J Chem Theory Comput 2022; 18:3718-3736. [PMID: 35377648 PMCID: PMC9202308 DOI: 10.1021/acs.jctc.2c00063] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Indexed: 12/12/2022]
Abstract
In this contribution we present a mixed quantum-classical dynamical approach for the computation of vibronic absorption spectra of molecular aggregates and their nonadiabatic dynamics, taking into account the coupling between local excitations (LE) and charge-transfer (CT) states. The approach is based on an adiabatic (Ad) separation between the soft degrees of freedom (DoFs) of the system and the stiff vibrations, which are described by the quantum dynamics (QD) of wave packets (WPs) moving on the coupled potential energy surfaces (PESs) of the LE and CT states. These PESs are described with a linear vibronic coupling (LVC) Hamiltonian, parameterized by an overlap-based diabatization on the grounds of time-dependent density functional theory computations. The WPs time evolution is computed with the multiconfiguration time-dependent Hartree method, using effective modes defined through a hierarchical representation of the LVC Hamiltonian. The soft DoFs are sampled with classical molecular dynamics (MD), and the coupling between the slow and fast DoFs is included by recomputing the key parameters of the LVC Hamiltonians, specifically for each MD configuration. This method, named Ad-MD|gLVC, is applied to a perylene diimide (PDI) dimer in acetonitrile and water solutions, and it is shown to accurately reproduce the change in the vibronic features of the absorption spectrum upon aggregation. Moreover, the microscopic insight offered by the MD trajectories allows for a detailed understanding of the role played by the fluctuation of the aggregate structure on the shape of the vibronic spectrum and on the population of LE and CT states. The nonadiabatic QD predicts an extremely fast (∼50 fs) energy transfer between the two LEs. CT states have only a moderate effect on the absorption spectrum, despite the fact that after photoexcitation they are shown to acquire a fast and non-negligible population, highlighting their relevance in dictating the charge separation and transport in PDI-based optical devices.
Collapse
Affiliation(s)
- Alekos Segalina
- Université
de Lorraine and CNRS, LPCT, UMR 7019, F-54000 Nancy, France
| | - Daniel Aranda
- Instituto
de Ciencia Molecular (ICMol), Universidad
de Valencia, Catedrático J. Beltrán 2, 46980 Paterna, Valencia, Spain
| | - James A. Green
- Consiglio
Nazionale delle Ricerche, Istituto di Biostrutture
e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy
| | - Vito Cristino
- Dipartimento
di Scienze Chimiche, Farmaceutiche ed Agrarie, Via Fossato di Mortara 17, 44121 Ferrara, Italy
| | - Stefano Caramori
- Dipartimento
di Scienze Chimiche, Farmaceutiche ed Agrarie, Via Fossato di Mortara 17, 44121 Ferrara, Italy
| | - Giacomo Prampolini
- Istituto
di Chimica dei Composti Organo Metallici, Consiglio Nazionale delle Ricerche, (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi
1, I-56124 Pisa, Italy
| | | | - Fabrizio Santoro
- Istituto
di Chimica dei Composti Organo Metallici, Consiglio Nazionale delle Ricerche, (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi
1, I-56124 Pisa, Italy
| |
Collapse
|
5
|
Jang SJ, Burghardt I, Hsu CP, Bardeen CJ. Excitons: Energetics and spatiotemporal dynamics. J Chem Phys 2021; 155:200401. [PMID: 34852498 DOI: 10.1063/5.0075292] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Seogjoo J Jang
- Department of Chemistry and Biochemistry, Queens College, City University of New York, 65-30 Kissena Boulevard, Queens, New York 11367, USA and PhD Programs in Chemistry and Physics, and Initiative for the Theoretical Sciences, Graduate Center, City University of New York, 365 Fifth Avenue, New York, New York 10016, USA
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Frankfurt, Germany
| | - Chao-Ping Hsu
- Institute of Chemistry, Academia Sinica, Taipei 115, Taiwan and Physics Division, National Center for Theoretical Sciences, Taipei 106, Taiwan
| | - Christopher J Bardeen
- Department of Chemistry, University of California, Riverside, Riverside, California 92521, USA
| |
Collapse
|
6
|
Provazza J, Tempelaar R, Coker DF. Analytic and numerical vibronic spectra from quasi-classical trajectory ensembles. J Chem Phys 2021; 155:014108. [PMID: 34241392 DOI: 10.1063/5.0053735] [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/14/2022] Open
Abstract
The truncated Wigner approximation to quantum dynamics in phase space is explored in the context of computing vibronic line shapes for monomer linear optical spectra. We consider multiple model potential forms including a shifted harmonic oscillator with both equal and unequal frequencies on the ground and excited state potentials as well as a shifted Morse potential model. For the equal-frequency shifted harmonic oscillator model, we derive an analytic expression for the exact vibronic line shape that emphasizes the importance of using a quantum mechanical distribution of phase space initial conditions. For the unequal-frequency shifted harmonic oscillator model, we are no longer able to obtain an exact expression for the vibronic line shape in terms of independent deterministic classical trajectories. We show how one can rigorously account for corrections to the truncated Wigner approximation through nonlinear responses of the line shape function to momentum fluctuations along a classical trajectory and demonstrate the qualitative improvement in the resulting spectrum when the leading-order quantum correction is included. Finally, we numerically simulate absorption spectra of a highly anharmonic shifted Morse potential model. We find that, while finite quantization and the dissociation limit are captured with reasonable accuracy, there is a qualitative breakdown of the quasi-classical trajectory ensemble's ability to describe the vibronic line shape when the relative shift in Morse potentials becomes large. The work presented here provides clarity on the origin of unphysical negative features known to contaminate absorption spectra computed with quasi-classical trajectory ensembles.
Collapse
Affiliation(s)
- Justin Provazza
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - Roel Tempelaar
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - David F Coker
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| |
Collapse
|