1
|
Wang YS, Zhong Manis JX, Rohan MC, Orlando TM, Kretchmer JS. Modeling Intermolecular Coulombic Decay with Non-Hermitian Real-Time Time-Dependent Density Functional Theory. J Phys Chem Lett 2024:7806-7813. [PMID: 39052307 DOI: 10.1021/acs.jpclett.4c01146] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2024]
Abstract
In this work, we investigate the capability of using real-time time-dependent density functional theory (RT-TDDFT) in conjunction with a complex absorbing potential (CAP) to simulate the intermolecular Coulombic decay (ICD) processes following the ionization of an inner-valence electron. We examine the ICD dynamics in a series of noncovalent bonded dimer systems, including hydrogen-bonded and purely van der Waals (VdW)-bonded systems. In comparison to previous work, we show that RT-TDDFT simulations with a CAP correctly capture the ICD phenomenon in systems exhibiting a stronger binding energy. The calculated time scales for ICD of the studied systems are in the range of 5-50 fs, in agreement with previous studies. However, there is a breakdown in the accuracy of the methodology for the pure VdW-bonded systems. Overall, the presented RT-TDDFT/CAP methodology provides a powerful tool for differentiating between competing electronic relaxation pathways following inner-valence or core ionization without necessitating any a priori assumptions.
Collapse
Affiliation(s)
- Yi-Siang Wang
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - James X Zhong Manis
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Matthew C Rohan
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Thomas M Orlando
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Joshua S Kretchmer
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| |
Collapse
|
2
|
Shepard C, Zhou R, Bost J, Carney TE, Yao Y, Kanai Y. Efficient exact exchange using Wannier functions and other related developments in planewave-pseudopotential implementation of RT-TDDFT. J Chem Phys 2024; 161:024111. [PMID: 38984957 DOI: 10.1063/5.0211238] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2024] [Accepted: 06/19/2024] [Indexed: 07/11/2024] Open
Abstract
The plane-wave pseudopotential (PW-PP) formalism is widely used for the first-principles electronic structure calculation of extended periodic systems. The PW-PP approach has also been adapted for real-time time-dependent density functional theory (RT-TDDFT) to investigate time-dependent electronic dynamical phenomena. In this work, we detail recent advances in the PW-PP formalism for RT-TDDFT, particularly how maximally localized Wannier functions (MLWFs) are used to accelerate simulations using the exact exchange. We also discuss several related developments, including an anti-Hermitian correction for the time-dependent MLWFs (TD-MLWFs) when a time-dependent electric field is applied, the refinement procedure for TD-MLWFs, comparison of the velocity and length gauge approaches for applying an electric field, and elimination of long-range electrostatic interaction, as well as usage of a complex absorbing potential for modeling isolated systems when using the PW-PP formalism.
Collapse
Affiliation(s)
- Christopher Shepard
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Ruiyi Zhou
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - John Bost
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Thomas E Carney
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
- Department of Physics and Astronomy, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| |
Collapse
|
3
|
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
|
4
|
Schreder L, Luber S. Propagated (fragment) Pipek-Mezey Wannier functions in real-time time-dependent density functional theory. J Chem Phys 2024; 160:214117. [PMID: 38832736 DOI: 10.1063/5.0203442] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2024] [Accepted: 05/14/2024] [Indexed: 06/05/2024] Open
Abstract
Localization procedures are an important tool for analysis of complex systems in quantum chemistry, since canonical molecular orbitals are delocalized and can, therefore, be difficult to align with chemical intuition and obscure information at the local level of the system. This especially applies to calculations obeying periodic boundary conditions. The most commonly used approach to localization is Foster-Boys Wannier functions, which use a unitary transformation to jointly minimize the second moment of the orbitals. This procedure has proven to be robust and fast but has a side effect of often mixing σ- and π-type orbitals. σ/π-separation is achieved by the Pipek-Mezey Wannier function (PMWF) approach [Lehtola and Jónsson, J. Chem. Theory Comput. 10, 642 (2014) and Jónsson et al., J. Chem. Theory Comput. 13, 460 (2017)], which defines the spread functional in terms of partial charges instead. We have implemented a PMWF algorithm in the CP2K software package using the Cardoso-Souloumiac algorithm to enable their application to real-time time-dependent density functional theory. The method is demonstrated on stacked CO2 molecules, linear acetylenic carbon, boron and nitrogen co-doped graphene, and nitrogen-vacancy doped diamond. Finally, we discuss its computational scaling and recent efforts to improve it with fragment approaches.
Collapse
Affiliation(s)
- Lukas Schreder
- University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Sandra Luber
- University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| |
Collapse
|
5
|
Korsaye FA, Perrella F, Petrone A, Adamo C, Rega N, Ciofini I. Monitoring Density Redistribution at the Excited State in a Dual Emitting Molecule: An Analysis Based on Real-Time Density Functional Theory and Density Descriptors. J Phys Chem A 2024; 128:4324-4334. [PMID: 38758031 DOI: 10.1021/acs.jpca.4c00934] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2024]
Abstract
In this work, we computed and analyzed, by means of density-based descriptors, the real-time evolution of both the locally excited (LE) and charge-transfer (CT) excited states for the planar and twisted conformations of the DMABN (4-(N,N-dimethylamino)benzonitrile) molecule using real-time time-dependent density functional theory (DFT) and three different exchange-correlation energy functionals (EXC) belonging to the same family (the PBE one). Our results based on the analysis of density-based descriptors show that the underlying EXC modifies the evolution in time of the density. In particular, comparing the frequency of density reorganization computed with the three functionals (PBE, PBE0, and LC-PBE), we found that the frequency of electronic interconversion of the individual determinants involved during the dynamics increases from PBE to PBE0 and to LC-PBE. This allows us to show that there is a correlation between the delocalization of the electronic density and the frequency of reorganization. In particular, the greater the mean hole-electron distance during the dynamics, the lower is the frequency of density reorganization.
Collapse
Affiliation(s)
- Feven-Alemu Korsaye
- Chimie ParisTech, PSL University, CNRS, Institute of Chemistry for Life and Health Sciences, F-75005 Paris, France
| | - Fulvio Perrella
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italia
| | - Alessio Petrone
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italia
- Dipartimento di Scienze Chimiche, Università degli Studi di Napoli Federico II, via Cintia 21, I-80126 Napoli, Italia
- Istituto Nazionale di Fisica Nucleare, Unità di Napoli, via Cintia 21, I-80126 Napoli, Italia
| | - Carlo Adamo
- Chimie ParisTech, PSL University, CNRS, Institute of Chemistry for Life and Health Sciences, F-75005 Paris, France
| | - Nadia Rega
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italia
- Dipartimento di Scienze Chimiche, Università degli Studi di Napoli Federico II, via Cintia 21, I-80126 Napoli, Italia
- Istituto Nazionale di Fisica Nucleare, Unità di Napoli, via Cintia 21, I-80126 Napoli, Italia
| | - Ilaria Ciofini
- Chimie ParisTech, PSL University, CNRS, Institute of Chemistry for Life and Health Sciences, F-75005 Paris, France
| |
Collapse
|
6
|
Schreder L, Luber S. Chiral Spectroscopy of Bulk Systems with Propagated Localized Orbitals. J Chem Theory Comput 2024; 20:3894-3910. [PMID: 38661175 DOI: 10.1021/acs.jctc.4c00235] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
We present approaches for the simulation of electronic circular dichroism, Raman, and Raman optical activity (ROA) spectra for isolated and periodic systems as well as subsystem analysis thereof. The method is based on the use of time-dependent maximally localized Wannier functions in the CP2K package and accounts for origin dependencies inherent to the Gaussian and plane wave with pseudopotentials approach as well as the origin dependence of the magnetic dipole and electric quadrupole operators. Tests on the H-bonded enantiomers of alanine by harmonic normal-mode analysis and on an aqueous solution of l-alanine by ab initio molecular dynamics obeying periodic boundary conditions (PBCs) are presented as total and subsystem-resolved spectra. To our knowledge, this is the first instance of an ROA spectrum derived from real-time propagation obeying PBCs and the first ROA simulation considering off-, pre-, and on-resonance effects within PBCs.
Collapse
Affiliation(s)
- Lukas Schreder
- University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Sandra Luber
- University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| |
Collapse
|
7
|
Woźniak AP, Moszyński R. Modeling of High-Harmonic Generation in the C 60 Fullerene Using Ab Initio, DFT-Based, and Semiempirical Methods. J Phys Chem A 2024; 128:2683-2702. [PMID: 38534023 PMCID: PMC11017253 DOI: 10.1021/acs.jpca.3c07865] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Revised: 03/04/2024] [Accepted: 03/13/2024] [Indexed: 03/28/2024]
Abstract
We report calculations of the high-harmonic generation spectra of the C60 fullerene molecule carried out by employing a diverse set of real-time time-dependent quantum chemical methods. All methodologies involve expanding the propagated electronic wave function in bases consisting of the ground and singly excited time-independent eigenstates obtained through the solution of the corresponding linear-response equations. We identify the correlation and exchange effect in the spectra by comparing the results from methods relying on the Hartree-Fock reference determinant with those obtained using approaches based on the density functional theory with different exchange-correlation functionals. The effect of the full random-phase approximation treatment of the excited electronic states is also analyzed and compared with the configuration interaction singles and the Tamm-Dancoff approximation. We also showcase the fact that the real-time extension of the semiempirical method INDO/S can be effectively applied for an approximate description of laser-driven dynamics in large systems.
Collapse
Affiliation(s)
| | - Robert Moszyński
- Faculty of Chemistry, University
of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| |
Collapse
|
8
|
Xu J, Carney TE, Zhou R, Shepard C, Kanai Y. Real-Time Time-Dependent Density Functional Theory for Simulating Nonequilibrium Electron Dynamics. J Am Chem Soc 2024; 146:5011-5029. [PMID: 38362887 DOI: 10.1021/jacs.3c08226] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2024]
Abstract
The explicit real-time propagation approach for time-dependent density functional theory (RT-TDDFT) has increasingly become a popular first-principles computational method for modeling various time-dependent electronic properties of complex chemical systems. In this Perspective, we provide a nontechnical discussion of how this first-principles simulation approach has been used to gain novel physical insights into nonequilibrium electron dynamics phenomena in recent years. Following a concise overview of the RT-TDDFT methodology from a practical standpoint, we discuss our recent studies on the electronic stopping of DNA in water and the Floquet topological phase as examples. Our discussion focuses on how RT-TDDFT simulations played a unique role in deriving new scientific understandings. We then discuss existing challenges and some new advances at the frontier of RT-TDDFT method development for studying increasingly complex dynamic phenomena and systems.
Collapse
Affiliation(s)
- Jianhang Xu
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Thomas E Carney
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Ruiyi Zhou
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Christopher Shepard
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
- Department of Physics and Astronomy, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
9
|
Crisci L, Coppola F, Petrone A, Rega N. Tuning ultrafast time-evolution of photo-induced charge-transfer states: A real-time electronic dynamics study in substituted indenotetracene derivatives. J Comput Chem 2024; 45:210-221. [PMID: 37706600 DOI: 10.1002/jcc.27231] [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/20/2023] [Revised: 08/31/2023] [Accepted: 09/05/2023] [Indexed: 09/15/2023]
Abstract
Photo-induced charge transfer (CT) states are pivotal in many technological and biological processes. A deeper knowledge of such states is mandatory for modeling the charge migration dynamics. Real-time time-dependent density functional theory (RT-TD-DFT) electronic dynamics simulations are employed to explicitly observe the electronic density time-evolution upon photo-excitation. Asymmetrically substituted indenotetracene molecules, given their potential application as n-type semiconductors in organic photovoltaic materials, are here investigated. Effects of substituents with different electron-donating characters are analyzed in terms of the overall electronic energy spacing and resulting ultrafast CT dynamics through linear response (LR-)TD-DFT and RT-TD-DFT based approaches. The combination of the computational techniques here employed provided direct access to the electronic density reorganization in time and to its spatial and rational representation in terms of molecular orbital occupation time evolution. Such results can be exploited to design peculiar directional charge dynamics, crucial when photoactive materials are used for light-harvesting applications.
Collapse
Affiliation(s)
- Luigi Crisci
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Normale Superiore di Pisa, Pisa, Italy
| | | | - Alessio Petrone
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Superiore Meridionale, Naples, Italy
- Istituto Nazionale Di Fisica Nucleare, Sezione di Napoli, Complesso Universitario di M.S. Angelo ed. 6, Naples, Italy
| | - Nadia Rega
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Superiore Meridionale, Naples, Italy
- Istituto Nazionale Di Fisica Nucleare, Sezione di Napoli, Complesso Universitario di M.S. Angelo ed. 6, Naples, Italy
| |
Collapse
|
10
|
Buttarazzi E, Perrella F, Rega N, Petrone A. Watching the Interplay between Photoinduced Ultrafast Charge Dynamics and Nuclear Vibrations. J Chem Theory Comput 2023; 19:8751-8766. [PMID: 37991892 PMCID: PMC10720350 DOI: 10.1021/acs.jctc.3c00855] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2023] [Revised: 11/02/2023] [Accepted: 11/03/2023] [Indexed: 11/24/2023]
Abstract
Here is presented the ultrafast hole-electron dynamics of photoinduced metal to ligand charge-transfer (MLCT) states in a Ru(II) complex, [Ru(dcbpy)2(NCS)2]4- (dcbpy = 4,4'-dicarboxy-2,2'-bipyridine), a photoactive molecule employed in dye sensitized solar cells. Via cutting-edge computational techniques, a tailored computational protocol is here presented and developed to provide a detailed analysis of the electronic manifold coupled with nuclear vibrations to better understand the nonradiative pathways and the resulting overall dye performances in light-harvesting processes (electron injection). Thus, the effects of different vibrational modes were investigated on both the electronic levels and charge transfer dynamics through a theoretical-computational approach. First, the linear response time-dependent density functional (LR-TDDFT) formalism was employed to characterize excitation energies and spacing among electronic levels (the electronic layouts). Then, to understand the ultrafast (femtosecond) charge dynamics on the molecular scale, we relied on the nonperturbative mean-field quantum electronic dynamics via real-time (RT-) TDDFT. Three vibrational modes were selected, representative for collective nuclear movements that can have a significant influence on the electronic structure: two involving NCS- ligands and one involving dcbpy ligands. As main results, we observed that such MLCT states, under vibrational distortions, are strongly affected and a faster interligand electron transfer mechanism is observed along with an increasing MLCT character of the adiabatic electronic states approaching closer in energy due to the vibrations. Such findings can help both in providing a molecular picture of multidimensional vibro-electronic spectroscopic techniques, used to characterize ultrafast coherent and noncoherent dynamics of complex systems, and to improve dye performances with particular attention to the study of energy or charge transport processes and vibronic couplings.
Collapse
Affiliation(s)
- Edoardo Buttarazzi
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di Monte S. Angelo, Via Cintia 21, I-80126 Napoli, Italy
| | - Fulvio Perrella
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
| | - Nadia Rega
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di Monte S. Angelo, Via Cintia 21, I-80126 Napoli, Italy
- Istituto
Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario
di Monte S. Angelo ed. 6, Via Cintia, I-80126 Napoli, Italy
| | - Alessio Petrone
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di Monte S. Angelo, Via Cintia 21, I-80126 Napoli, Italy
- Istituto
Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario
di Monte S. Angelo ed. 6, Via Cintia, I-80126 Napoli, Italy
| |
Collapse
|
11
|
Peyton BG, Stewart ZJ, Weidman JD, Wilson AK. Tailoring light-induced charge transfer and intersystem crossing in FeCO using time-dependent spin-orbit configuration interaction. J Chem Phys 2023; 159:204108. [PMID: 38014783 DOI: 10.1063/5.0173529] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Accepted: 11/02/2023] [Indexed: 11/29/2023] Open
Abstract
Real-time (RT) electronic structure methods provide a natural framework for describing light-matter interactions in arbitrary time-dependent electromagnetic fields (EMF). Optically induced excited state transitions are of particular interest, which require tuned EMF to drive population transfer to and from the specific state(s) of interest. Intersystem crossing, or spin-flip, may be driven through shaped EMF or laser pulses. These transitions can result in long-lived "spin-trapped" excited states, which are especially useful for materials requiring charge separation or protracted excited state lifetimes. Time-dependent configuration interaction (TDCI) is unique among RT methods in that it may be implemented in a basis of eigenstates, allowing for rapid propagation of the time-dependent Schrödinger equation. The recent spin-orbit TDCI (TD-SOCI) enables a real-time description of spin-flip dynamics in an arbitrary EMF and, therefore, provides an ideal framework for rational pulse design. The present study explores the mechanism of multiple spin-flip pathways for a model transition metal complex, FeCO, using shaped pulses designed to drive controlled intersystem crossing and charge transfer. These results show that extremely tunable excited state dynamics can be achieved by considering the dipole transition matrix elements between the states of interest.
Collapse
Affiliation(s)
- Benjamin G Peyton
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Zachary J Stewart
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Jared D Weidman
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
12
|
Dyck O, Almutlaq J, Lingerfelt D, Swett JL, Oxley MP, Huang B, Lupini AR, Englund D, Jesse S. Direct imaging of electron density with a scanning transmission electron microscope. Nat Commun 2023; 14:7550. [PMID: 37985658 PMCID: PMC10662251 DOI: 10.1038/s41467-023-42256-9] [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/15/2022] [Accepted: 07/27/2023] [Indexed: 11/22/2023] Open
Abstract
Recent studies of secondary electron (SE) emission in scanning transmission electron microscopes suggest that material's properties such as electrical conductivity, connectivity, and work function can be probed with atomic scale resolution using a technique known as secondary electron e-beam-induced current (SEEBIC). Here, we apply the SEEBIC imaging technique to a stacked 2D heterostructure device to reveal the spatially resolved electron density of an encapsulated WSe2 layer. We find that the double Se lattice site shows higher emission than the W site, which is at odds with first-principles modelling of valence ionization of an isolated WSe2 cluster. These results illustrate that atomic level SEEBIC contrast within a single material is possible and that an enhanced understanding of atomic scale SE emission is required to account for the observed contrast. In turn, this suggests that, in the future, subtle information about interlayer bonding and the effect on electron orbitals could be directly revealed with this technique.
Collapse
Affiliation(s)
- Ondrej Dyck
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, TN, USA.
| | | | - David Lingerfelt
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, TN, USA
| | - Jacob L Swett
- Biodesign Institute, Arizona State University, Tempe, 87287, AZ, USA
| | - Mark P Oxley
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, TN, USA
| | - Bevin Huang
- Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Andrew R Lupini
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, TN, USA
| | - Dirk Englund
- Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Stephen Jesse
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, TN, USA
| |
Collapse
|
13
|
Hauge E, Kristiansen HE, Konecny L, Kadek M, Repisky M, Pedersen TB. Cost-Efficient High-Resolution Linear Absorption Spectra through Extrapolating the Dipole Moment from Real-Time Time-Dependent Electronic-Structure Theory. J Chem Theory Comput 2023; 19:7764-7775. [PMID: 37874968 PMCID: PMC10653104 DOI: 10.1021/acs.jctc.3c00727] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Revised: 10/02/2023] [Accepted: 10/03/2023] [Indexed: 10/26/2023]
Abstract
We present a novel function fitting method for approximating the propagation of the time-dependent electric dipole moment from real-time electronic structure calculations. Real-time calculations of the electronic absorption spectrum require discrete Fourier transforms of the electric dipole moment. The spectral resolution is determined by the total propagation time, i.e., the trajectory length of the dipole moment, causing a high computational cost. Our developed method uses function fitting on shorter trajectories of the dipole moment, achieving arbitrary spectral resolution through extrapolation. Numerical testing shows that the fitting method can reproduce high-resolution spectra by using short dipole trajectories. The method converges with as little as 100 a.u. dipole trajectories for some systems, though the difficulty converging increases with the spectral density. We also introduce an error estimate of the fit, reliably assessing its convergence and hence the quality of the approximated spectrum.
Collapse
Affiliation(s)
- Eirill Hauge
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
- Department
of Numerical Analysis and Scientific Computing, Simula Research Laboratory, Kristian Augusts Gate 23, 0164 Oslo, Norway
| | - Håkon Emil Kristiansen
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| | - Lukas Konecny
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø—The Arctic University
of Norway, N-9037 Tromsø, Norway
- Center
for Free Electron Laser, Max Planck Institute
for the Structure and Dynamics of Matter Science, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Marius Kadek
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø—The Arctic University
of Norway, N-9037 Tromsø, Norway
- Department
of Physics, Northeastern University, Boston, Massachusetts 02115, United States
| | - Michal Repisky
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø—The Arctic University
of Norway, N-9037 Tromsø, Norway
- Department
of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University, Ilkovicova 6, SK-84215 Bratislava, Slovakia
| | - Thomas Bondo Pedersen
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| |
Collapse
|
14
|
Mokkath JH. Plasmon induced hot carrier generation in a pyridine@Au20 composite. Phys Chem Chem Phys 2023; 25:28750-28760. [PMID: 37850351 DOI: 10.1039/d3cp02907g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2023]
Abstract
Using time-dependent density functional theory calculations, we have investigated the generation of hot carriers (HCs) in a system comprising a pyridine molecule and a tetrahedral Au20 plasmonic cluster. Our findings indicate that the decay of the localized surface plasmon resonance (LSPR) induced in the pyridine@Au20 system by a laser pulse facilitates the direct transfer of hot electrons from the occupied states of the Au20 cluster to the unoccupied molecular orbitals of pyridine. Notably, we have identified that the interparticle gap distance between the Au20 cluster and the pyridine molecule plays a critical role in controlling the generation of HCs. By precisely controlling the interaction between the plasmonic cluster and the molecule, we can effectively manipulate the energy distribution of the generated HCs. These insights have the potential to drive advancements in the development of more efficient systems for plasmonic catalysis.
Collapse
Affiliation(s)
- Junais Habeeb Mokkath
- Quantum Nanophotonics Simulations Lab, Department of Physics, Kuwait College of Science and Technology, Doha Area, 7th Ring Road, P.O. Box 27235, Kuwait.
| |
Collapse
|
15
|
Peyton BG, Wang Z, Crawford TD. Reduced Scaling Real-Time Coupled Cluster Theory. J Phys Chem A 2023; 127:8486-8499. [PMID: 37782945 DOI: 10.1021/acs.jpca.3c05151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/04/2023]
Abstract
Real-time coupled cluster (CC) methods have several advantages over their frequency-domain counterparts, namely, response and equation of motion CC theories. Broadband spectra, strong fields, and pulse manipulation allow for the simulation of complex spectroscopies that are unreachable using frequency-domain approaches. Due to the high-order polynomial scaling, the required numerical time propagation of the CC residual expressions is a computationally demanding process. This scaling may be reduced by local correlation schemes, which aim to reduce the size of the (virtual) orbital space by truncation according to user-defined parameters. We present the first application of local correlation to real-time CC. As in previous studies of locally correlated frequency-domain CC, traditional local correlation schemes are of limited utility for field-dependent properties; however, a perturbation-aware scheme proves promising. A detailed analysis of the amplitude dynamics suggests that the main challenge is a strong time dependence of the wave function sparsity.
Collapse
Affiliation(s)
- Benjamin G Peyton
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Zhe Wang
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - T Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
16
|
Herbert JM, Zhu Y, Alam B, Ojha AK. Time-Dependent Density Functional Theory for X-ray Absorption Spectra: Comparing the Real-Time Approach to Linear Response. J Chem Theory Comput 2023; 19:6745-6760. [PMID: 37708349 DOI: 10.1021/acs.jctc.3c00673] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/16/2023]
Abstract
We simulate X-ray absorption spectra at elemental K-edges using time-dependent density functional theory (TDDFT) in both its conventional linear-response implementation and its explicitly time-dependent or "real-time" formulation. Real-time TDDFT simulations enable broadband spectra calculations without the need to invoke frozen occupied orbitals ("core/valence separation"), but we find that these spectra are often contaminated by transitions to the continuum that originate from lower-energy core and semicore orbitals. This problem becomes acute in triple-ζ basis sets, although it is sometimes sidestepped in double-ζ basis sets. Transitions to the continuum acquire surprisingly large dipole oscillator strengths, leading to spectra that are difficult to interpret. Meaningful spectra can be recovered by means of a filtering technique that decomposes the spectrum into contributions from individual occupied orbitals, and the same procedure can be used to separate L- and K-edge spectra arising from different elements within a given molecule. In contrast, conventional linear-response TDDFT requires core/valence separation but is free of these artifacts. It is also significantly more efficient than the real-time approach, even when hundreds of individual states are needed to reproduce near-edge absorption features and even when Padé approximants are used to reduce the real-time simulations to just 2-4 fs of time propagation. Despite the cost, the real-time approach may be useful to examine the validity of the core/valence separation approximation.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
- Chemical Physics Graduate Program, The Ohio State University, Columbus, Ohio 43210, United States
| | - Ying Zhu
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
- Chemical Physics Graduate Program, The Ohio State University, Columbus, Ohio 43210, United States
| | - Bushra Alam
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Avik Kumar Ojha
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
17
|
Richter R, Aschebrock T, Schelter I, Kümmel S. Meta-generalized gradient approximations in time dependent generalized Kohn-Sham theory: Importance of the current density correction. J Chem Phys 2023; 159:124117. [PMID: 38127400 DOI: 10.1063/5.0167972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2023] [Accepted: 09/07/2023] [Indexed: 12/23/2023] Open
Abstract
We revisit the use of Meta-Generalized Gradient Approximations (mGGAs) in time-dependent density functional theory, reviewing conceptual questions and solving the generalized Kohn-Sham equations by real-time propagation. After discussing the technical aspects of using mGGAs in combination with pseudopotentials and comparing real-space and basis set results, we focus on investigating the importance of the current-density based gauge invariance correction. For the two modern mGGAs that we investigate in this work, TASK and r2SCAN, we observe that for some systems, the current density correction leads to negligible changes, but for others, it changes excitation energies by up to 40% and more than 0.8 eV. In the cases that we study, the agreement with the reference data is improved by the current density correction.
Collapse
Affiliation(s)
- Rian Richter
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| | - Thilo Aschebrock
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| | - Ingo Schelter
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| | - Stephan Kümmel
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| |
Collapse
|
18
|
Hawks A, Altman D, Faddis R, Wagner EM, Bell KJJ, Charland-Martin A, Collier GS. Relating Design and Optoelectronic Properties of 1,4-Dihydropyrrolo[3,2- b]pyrroles Bearing Biphenyl Substituents. J Phys Chem B 2023; 127:7352-7360. [PMID: 37561612 PMCID: PMC10461294 DOI: 10.1021/acs.jpcb.3c03061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Revised: 07/17/2023] [Indexed: 08/12/2023]
Abstract
Understanding the influence of peripheral functionality on optoelectronic properties of conjugated materials is an important task for the continued development of chromophores for myriad applications. Here, π-extended 1,4-dihydropyrrolo[3,2-b]pyrrole (DHPP) chromophores with varying electron-donating or electron-withdrawing capabilities were synthesized via Suzuki cross-coupling reactions, and the influence of functionality on optoelectronic properties was elucidated. First, chromophores display distinct differences in the UV-vis absorbance spectra measured via UV-vis absorbance spectroscopy in addition to changes in the onset of oxidation measured with cyclic voltammetry and differential pulse voltammetry. Solution oxidation studies found that variations in the electron-donating and -withdrawing capabilities result in different absorbance profiles of the radical cations that correspond to quantifiably different colors. In addition to fundamental insights into the molecular design of DHPP chromophores and their optoelectronic properties, two chromophores display high-contrast electrochromism, which makes them potentially compelling in electronic devices. Overall, this study represents the ability to fine-tune the optoelectronic properties of DHPP chromophores in their neutral and oxidized states and expands the understanding of structure-property relationships that will guide the continued development of DHPP-based materials.
Collapse
Affiliation(s)
- Allison
M. Hawks
- Department of Chemistry and
Biochemistry, Kennesaw State University, Kennesaw, Georgia 30144, United States
| | - Drake Altman
- Department of Chemistry and
Biochemistry, Kennesaw State University, Kennesaw, Georgia 30144, United States
| | - Ryan Faddis
- Department of Chemistry and
Biochemistry, Kennesaw State University, Kennesaw, Georgia 30144, United States
| | - Ethan M. Wagner
- Department of Chemistry and
Biochemistry, Kennesaw State University, Kennesaw, Georgia 30144, United States
| | - Kenneth-John J. Bell
- Department of Chemistry and
Biochemistry, Kennesaw State University, Kennesaw, Georgia 30144, United States
| | - Ariane Charland-Martin
- Department of Chemistry and
Biochemistry, Kennesaw State University, Kennesaw, Georgia 30144, United States
| | - Graham S. Collier
- Department of Chemistry and
Biochemistry, Kennesaw State University, Kennesaw, Georgia 30144, United States
| |
Collapse
|
19
|
Ranka K, Isborn CM. Size-dependent errors in real-time electron density propagation. J Chem Phys 2023; 158:2887545. [PMID: 37125706 DOI: 10.1063/5.0142515] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2023] [Accepted: 04/14/2023] [Indexed: 05/02/2023] Open
Abstract
Real-time (RT) electron density propagation with time-dependent density functional theory (TDDFT) or Hartree-Fock (TDHF) is one of the most popular methods to model the charge transfer in molecules and materials. However, both RT-TDHF and RT-TDDFT within the adiabatic approximation are known to produce inaccurate evolution of the electron density away from the ground state in model systems, leading to large errors in charge transfer and erroneous shifting of peaks in absorption spectra. Given the poor performance of these methods with small model systems and the widespread use of the methods with larger molecular and material systems, here we bridge the gap in our understanding of these methods and examine the size-dependence of errors in RT density propagation. We analyze the performance of RT density propagation for systems of increasing size during the application of a continuous resonant field to induce Rabi-like oscillations, during charge-transfer dynamics, and for peak shifting in simulated absorption spectra. We find that the errors in the electron dynamics are indeed size dependent for these phenomena, with the largest system producing the results most aligned with those expected from linear response theory. The results suggest that although the RT-TDHF and RT-TDDFT methods may produce severe errors for model systems, the errors in charge transfer and resonantly driven electron dynamics may be much less significant for more realistic, large-scale molecules and materials.
Collapse
Affiliation(s)
- Karnamohit Ranka
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| | - Christine M Isborn
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
20
|
Bhat V, Callaway CP, Risko C. Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials. Chem Rev 2023. [PMID: 37141497 DOI: 10.1021/acs.chemrev.2c00704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.
Collapse
Affiliation(s)
- Vinayak Bhat
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| | - Connor P Callaway
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| | - Chad Risko
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| |
Collapse
|
21
|
Folorunso AS, Mauger F, Hamer KA, Jayasinghe DD, Wahyutama IS, Ragains JR, Jones RR, DiMauro LF, Gaarde MB, Schafer KJ, Lopata K. Attochemistry Regulation of Charge Migration. J Phys Chem A 2023; 127:1894-1900. [PMID: 36791088 PMCID: PMC9986869 DOI: 10.1021/acs.jpca.3c00568] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2023]
Abstract
Charge migration (CM) is a coherent attosecond process that involves the movement of localized holes across a molecule. To determine the relationship between a molecule's structure and the CM dynamics it exhibits, we perform systematic studies of para-functionalized bromobenzene molecules (X-C6H4-R) using real-time time-dependent density functional theory. We initiate valence-electron dynamics by emulating rapid strong-field ionization leading to a localized hole on the bromine atom. The resulting CM, which takes on the order of 1 fs, occurs via an X localized → C6H4 delocalized → R localized mechanism. Interestingly, the hole contrast on the acceptor functional group increases with increasing electron-donating strength. This trend is well-described by the Hammett σ value of the group, which is a commonly used metric for quantifying the effect of functionalization on the chemical reactivity of benzene derivatives. These results suggest that simple attochemistry principles and a density-based picture can be used to predict and understand CM.
Collapse
Affiliation(s)
| | | | | | | | | | | | - Robert R Jones
- Department of Physics, University of Virginia, Charlottesville, Virginia 22904, United States
| | - Louis F DiMauro
- Department of Physics, The Ohio State University, Columbus, Ohio 43210, United States
| | | | | | | |
Collapse
|
22
|
Moitra T, Konecny L, Kadek M, Rubio A, Repisky M. Accurate Relativistic Real-Time Time-Dependent Density Functional Theory for Valence and Core Attosecond Transient Absorption Spectroscopy. J Phys Chem Lett 2023; 14:1714-1724. [PMID: 36757216 PMCID: PMC9940299 DOI: 10.1021/acs.jpclett.2c03599] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2022] [Accepted: 02/03/2023] [Indexed: 06/18/2023]
Abstract
First principles theoretical modeling of out-of-equilibrium processes observed in attosecond pump-probe transient absorption spectroscopy (TAS) triggering pure electron dynamics remains a challenging task, especially for heavy elements and/or core excitations containing fingerprints of scalar and spin-orbit relativistic effects. To address this, we formulate a methodology for simulating TAS within the relativistic real-time, time-dependent density functional theory (RT-TDDFT) framework, for both the valence and core energy regimes. Especially for TAS, full four-component (4c) RT simulations are feasible but computationally demanding. Therefore, in addition to the 4c approach, we also introduce the atomic mean-field exact two-component (amfX2C) Hamiltonian accounting for one- and two-electron picture-change corrections within RT-TDDFT. amfX2C preserves the accuracy of the parent 4c method at a fraction of its computational cost. Finally, we apply the methodology to study valence and near-L2,3-edge TAS processes of experimentally relevant systems and provide additional physical insights using relativistic nonequilibrium response theory.
Collapse
Affiliation(s)
- Torsha Moitra
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, 9037 Tromsø, Norway
| | - Lukas Konecny
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, 9037 Tromsø, Norway
- Max
Planck Institute for the Structure and Dynamics of Matter, Center for Free Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Marius Kadek
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, 9037 Tromsø, Norway
- Department
of Physics, Northeastern University, Boston, Massachusetts 02115, United States
- Algorithmiq
Ltd., Kanavakatu 3C, FI-00160 Helsinki, Finland
| | - Angel Rubio
- Max
Planck Institute for the Structure and Dynamics of Matter, Center for Free Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
- Center
for Computational Quantum Physics (CCQ), The Flatiron Institute, 162 Fifth Avenue, New York New York 10010, United States
- Nano-Bio
Spectroscopy Group, Departamento de Física de Materiales, Universidad del País Vasco, 20018 San Sebastian, Spain
| | - Michal Repisky
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, 9037 Tromsø, Norway
- Department
of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University, 84104 Bratislava, Slovakia
| |
Collapse
|
23
|
Kanungo B, Rufus ND, Gavini V. Efficient All-Electron Time-Dependent Density Functional Theory Calculations Using an Enriched Finite Element Basis. J Chem Theory Comput 2023; 19:978-991. [PMID: 36656153 DOI: 10.1021/acs.jctc.2c01052] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
We present an efficient and systematically convergent approach to all-electron real-time time-dependent density functional theory (TDDFT) calculations using a mixed basis, termed as enriched finite element (EFE) basis. The EFE basis augments the classical finite element basis (CFE) with a compactly supported numerical atom-centered basis, obtained from atomic ground-state DFT calculations. Particularly, we orthogonalize the enrichment functions with respect to the classical finite element basis to ensure good conditioning of the resultant basis. We employ the second-order Magnus propagator in conjunction with an adaptive Krylov subspace method for efficient time evolution of the Kohn-Sham orbitals. We rely on a priori error estimates to guide our choice of an adaptive finite element mesh as well as the time step to be used in the TDDFT calculations. We observe close to optimal rates of convergence of the dipole moment with respect to spatial and temporal discretizations. Notably, we attain a 50-100 times speedup for the EFE basis over the CFE basis. We also demonstrate the efficacy of the EFE basis for both linear and nonlinear responses by studying the absorption spectra in sodium clusters, the linear to nonlinear response transition in the green fluorescence protein chromophore, and the higher harmonic generation in the magnesium dimer. Lastly, we attain good parallel scalability of our numerical implementation of the EFE basis for up to ∼1000 processors, using a benchmark system of a 50-atom sodium nanocluster.
Collapse
Affiliation(s)
- Bikash Kanungo
- Department of Mechanical Engineering, University of Michigan, Ann Arbor, Michigan48109, United States
| | - Nelson D Rufus
- Department of Mechanical Engineering, University of Michigan, Ann Arbor, Michigan48109, United States
| | - Vikram Gavini
- Department of Mechanical Engineering, University of Michigan, Ann Arbor, Michigan48109, United States.,Department of Materials Science and Engineering, University of Michigan, Ann Arbor, Michigan48109, United States
| |
Collapse
|
24
|
Perrella F, Petrone A, Rega N. Understanding Charge Dynamics in Dense Electronic Manifolds in Complex Environments. J Chem Theory Comput 2023; 19:626-639. [PMID: 36602443 PMCID: PMC9878732 DOI: 10.1021/acs.jctc.2c00794] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Indexed: 01/06/2023]
Abstract
Photoinduced charge transfer (CT) excited states and their relaxation mechanisms can be highly interdependent on the environment effects and the consequent changes in the electronic density. Providing a molecular interpretation of the ultrafast (subpicosecond) interplay between initial photoexcited states in such dense electronic manifolds in condensed phase is crucial for improving and understanding such phenomena. Real-time time-dependent density functional theory is here the method of choice to observe the charge density, explicitly propagated in an ultrafast time domain, along with all time-dependent properties that can be easily extracted from it. A designed protocol of analysis for real-time electronic dynamics to be applied to time evolving electronic density related properties to characterize both in time and in space CT dynamics in complex systems is here introduced and validated, proposing easy to be read cross-correlation maps. As case studies to test such tools, we present the photoinduced charge-transfer electronic dynamics of 5-benzyluracil, a mimic of nucleic acid/protein interactions, and the metal-to-ligand charge-transfer electronic dynamics in water solution of [Ru(dcbpy)2(NCS)2]4-, dcbpy = (4,4'-dicarboxy-2,2'-bipyridine), or "N34-", a dye sensitizer for solar cells. Electrostatic and explicit ab initio treatment of solvent molecules have been compared in the latter case, revealing the importance of the accurate modeling of mutual solute-solvent polarization on CT kinetics. We observed that explicit quantum mechanical treatment of solvent slowed down the charge carriers mobilities with respect to the gas-phase. When all water molecules were modeled instead as simpler embedded point charges, the electronic dynamics appeared enhanced, with a reduced hole-electron distance and higher mean velocities due to the close fixed charges and an artificially increased polarization effect. Such analysis tools and the presented case studies can help to unveil the influence of the electronic manifold, as well as of the finite temperature-induced structural distortions and the environment on the ultrafast charge motions.
Collapse
Affiliation(s)
- Fulvio Perrella
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di
M.S. Angelo, via Cintia 21, I-80126, Napoli, Italy
| | - Alessio Petrone
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di
M.S. Angelo, via Cintia 21, I-80126, Napoli, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138, Napoli, Italy
- Istituto
Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario di Monte S. Angelo ed. 6, via Cintia, 80126, Napoli, Italia
| | - Nadia Rega
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di
M.S. Angelo, via Cintia 21, I-80126, Napoli, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138, Napoli, Italy
- Istituto
Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario di Monte S. Angelo ed. 6, via Cintia, 80126, Napoli, Italia
- CRIB,
Centro Interdipartimentale di Ricerca sui Biomateriali, Piazzale Tecchio 80, I-80125, Napoli, Italy
| |
Collapse
|
25
|
Fyta M. Functionalized electrodes embedded in nanopores: read-out enhancement? Chem Asian J 2023; 18:e202200916. [PMID: 36372991 PMCID: PMC10107472 DOI: 10.1002/asia.202200916] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Revised: 11/12/2022] [Accepted: 11/12/2022] [Indexed: 11/16/2022]
Abstract
In this review, functionalized nanogaps embedded in nanopores are discussed in view of their high biosensitivity in detecting biomolecules, their length, type, and sequence. Specific focus is given on nanoelectrodes functionalized with tiny nanometer-sized diamond-like particles offering vast functionalization possibilities for gold junction electrodes. This choice of the functionalization, in turn, offers nucleotide-specific binding possibilities improving the detection signals arising from such functionalized electrodes potentially embedded in a nanopore. The review sheds light onto the use and enhancement of the tunnelling recognition in functionalized nanogaps towards sensing DNA nucleotides and mutation detection, providing important input for a practical realization.
Collapse
Affiliation(s)
- Maria Fyta
- Computational Biotechnology, RWTH-Aachen University, Worringerweg 3, 52072, Aachen, Germany
| |
Collapse
|
26
|
Woźniak AP, Lewenstein M, Moszyński R. Exploring the attosecond laser-driven electron dynamics in the hydrogen molecule with different real-time time-dependent configuration interaction approaches. ADVANCES IN QUANTUM CHEMISTRY 2023. [DOI: 10.1016/bs.aiq.2023.01.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/05/2023]
|
27
|
Classical Nuclear Motion: Comparison to Approaches with Quantum Mechanical Nuclear Motion. HYDROGEN 2022. [DOI: 10.3390/hydrogen4010002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
Ab initio molecular dynamics combines a classical description of nuclear motion with a density-functional description of the electronic cloud. This approach nicely describes chemical reactions. A possible conclusion is that a quantum mechanical description of nuclear motion is not needed. Using Occam’s razor, this means that, being the simpler approach, classical nuclear motion is preferable. In this paper, it is claimed that nuclear motion is classical, and this hypothesis will be tested in comparison to methods with quantum mechanical nuclear motion. In particular, we apply ab initio molecular dynamics to two photoreactions involving hydrogen. Hydrogen, as the lightest element, is often assumed to show quantum mechanical tunneling. We will see that the classical picture is fully sufficient. The quantum mechanical view leads to phenomena that are difficult to understand, such as the entanglement of nuclear motion. In contrast, it is easy to understand the simple classical picture which assumes that nuclear motion is steady and uniform unless a force is acting. Of course, such a hypothesis must be verified for many systems and phenomena, and this paper is one more step in this direction.
Collapse
|
28
|
Lingerfelt DB, Yoshimura A, Jakowski J, Ganesh P, Sumpter BG. Extracting Inelastic Scattering Cross Sections for Finite and Aperiodic Materials from Electronic Dynamics Simulations. J Chem Theory Comput 2022; 18:7093-7107. [PMID: 36375179 DOI: 10.1021/acs.jctc.2c00882] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Explicit time-dependent electronic structure theory methods are increasingly prevalent in the areas of condensed matter physics and quantum chemistry, with the broad-band optical absorptivity of molecular and small condensed-phase systems nowadays routinely studied with such approaches. In this paper, it is demonstrated that electronic dynamics simulations can similarly be employed to study cross sections for the scattering-induced electronic excitations probed in nonresonant inelastic X-ray scattering and momentum-resolved electron energy loss spectroscopies. A method is put forth for evaluating the electronic dynamic structure factor, which involves the application of a momentum boost-type perturbation and transformation of the resulting reciprocal space density fluctuations into the frequency domain. Good agreement is first demonstrated between the dynamic structure factor extracted from these electronic dynamics simulations and the corresponding transition matrix elements from linear response theory. The method is then applied to some extended (quasi)one-dimensional systems, for which the wave vector becomes a good quantum number in the thermodynamic limit. Finally, the dispersion of many-body excitations in a series of hydrogen-terminated graphene flakes (and twisted bilayers thereof) is investigated to highlight the utility of the presented approach for capturing morphology-dependent effects in the inelastic scattering cross sections of nanostructured and/or noncrystalline materials.
Collapse
Affiliation(s)
- David B Lingerfelt
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee37831, United States
| | - Anthony Yoshimura
- Lawrence Livermore National Laboratory, Livermore, California94550, United States
| | - Jacek Jakowski
- Computing and Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee37831, United States
| | - Panchapakesan Ganesh
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee37831, United States
| | - Bobby G Sumpter
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee37831, United States
| |
Collapse
|
29
|
Trepl T, Schelter I, Kümmel S. Analyzing Excitation-Energy Transfer Based on the Time-Dependent Density Functional Theory in Real Time. J Chem Theory Comput 2022; 18:6577-6587. [PMID: 36268773 DOI: 10.1021/acs.jctc.2c00600] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Excitation-energy transfer is a key step in processes such as photosynthesis that convert light into other forms of energy. Time-dependent density functional theory (DFT) in real time is ideal for the first-principles simulation of such processes due to its computational efficiency. We here demonstrate how real-time DFT can be used for analyzing excitation-energy transfer from first-principles. We discuss several measures of energy transfer that are based solely on the time-dependent density, are well founded in the DFT framework, allow for intuitive understanding and visualization, and reproduce important limiting cases of an analytical model. We demonstrate their usefulness in calculations for model systems, both with static nuclei and in the context of DFT-based Ehrenfest dynamics.
Collapse
Affiliation(s)
- T Trepl
- Theoretical Physics IV, University of Bayreuth, Bayreuth95440, Germany
| | - I Schelter
- Theoretical Physics IV, University of Bayreuth, Bayreuth95440, Germany
| | - S Kümmel
- Theoretical Physics IV, University of Bayreuth, Bayreuth95440, Germany
| |
Collapse
|
30
|
Mattiat J, Luber S. Comparison of Length, Velocity, and Symmetric Gauges for the Calculation of Absorption and Electric Circular Dichroism Spectra with Real-Time Time-Dependent Density Functional Theory. J Chem Theory Comput 2022; 18:5513-5526. [PMID: 36041170 DOI: 10.1021/acs.jctc.2c00644] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A velocity and symmetric gauge implementation for real-time time-dependent density functional theory (RT-TDDFT) in the CP2K package using a Gaussian and plane wave approach is presented, including the explicit gauge-transformed contributions due to the nonlocal part of pseudopotentials. Absorption spectra of gas-phase α-pinene are calculated in length and velocity gauges in the long-wavelength approximation for the application of a δ pulse in linear and full order. The velocity gauge implementation is also applied to a solvated uracil molecule to showcase its use within periodic boundary conditions (PBC). For the calculation of the expectation value of the electric dipole moment in PBC, both the velocity representation and the modern theory of polarization give equivalent absorption spectra if a distributed reference point is used for the nonlocal term of the velocity operator. The discussion of linear response theory takes place in a unified framework in terms of linear response functions in propagator notation, distinguishing the parts of the linear response functions associated with perturbation and response. To further investigate gauge dependence, electric circular dichroism (ECD) spectra of α-pinene were calculated either as magnetic response to an electric field perturbation, in length or velocity gauge, or as electric response to a magnetic field perturbation in the symmetric gauge. Both approaches, electric and magnetic perturbations, have been found to yield equivalent ECD spectra.
Collapse
Affiliation(s)
- Johann Mattiat
- Department of Chemistry, University of Zurich, Zurich 8057, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zurich, Zurich 8057, Switzerland
| |
Collapse
|
31
|
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
|
32
|
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
|
33
|
Lambropoulos K, Alvertis AM, Morphis A, Simserides C. Cyclo[18]carbon including zero-point motion: ground state, first singlet and triplet excitations, and hole transfer. Phys Chem Chem Phys 2022; 24:7779-7787. [PMID: 35293921 DOI: 10.1039/d2cp00343k] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Recent synthesis of cyclo[18]carbon has spurred increasing interest in carbon rings. We focus on a comparative inspection of ground and excited states, as well as of hole transfer properties of cumulenic and polyynic cyclo[18]carbon via Density Functional Theory (DFT), time-dependent DFT (TD-DFT) and real-time time-dependent DFT (RT-TDDFT). Zero-point vibrations are also accounted for, using a Monte Carlo sampling technique and a less exact, yet mode-resolved, quadratic approximation. The inclusion of zero-point vibrations leads to a red-shift on the HOMO-LUMO gap and the first singlet and triplet excitation energies of both conformations, correcting the values of the 'static' configurations by 9% to 24%. Next, we oxidize the molecule, creating a hole at one carbon atom. Hole transfer along polyynic cyclo[18]carbon is decreased in magnitude compared to its cumulenic counterpart and lacks the symmetric features the latter displays. Contributions by each mode to energy changes and hole transfer between diametrically opposed atoms vary, with specific bond-stretching modes being dominant.
Collapse
Affiliation(s)
- Konstantinos Lambropoulos
- Department of Physics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos GR-15784, Athens, Greece.
| | - Antonios M Alvertis
- Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, UK.,Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA. .,Department of Physics, University of California Berkeley, Berkeley, California 94720, USA
| | - Andreas Morphis
- Department of Physics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos GR-15784, Athens, Greece.
| | - Constantinos Simserides
- Department of Physics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos GR-15784, Athens, Greece.
| |
Collapse
|
34
|
Ghosal A, Roy AK. A real-time TDDFT scheme for strong-field interaction in Cartesian coordinate grid. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2022.139562] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
35
|
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
|
36
|
Durden AS, Levine BG. Floquet Time-Dependent Configuration Interaction for Modeling Ultrafast Electron Dynamics. J Chem Theory Comput 2022; 18:795-806. [DOI: 10.1021/acs.jctc.1c01009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Andrew S. Durden
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, United States
| | - Benjamin G. Levine
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, United States
| |
Collapse
|
37
|
Mattiat J, Luber S. Recent Progress in the Simulation of Chiral Systems with Real Time Propagation Methods. Helv Chim Acta 2021. [DOI: 10.1002/hlca.202100154] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Johann Mattiat
- Department of Chemistry University of Zurich Winterthurerstrasse 190 CH-8057 Zurich Switzerland
| | - Sandra Luber
- Department of Chemistry University of Zurich Winterthurerstrasse 190 CH-8057 Zurich Switzerland
| |
Collapse
|
38
|
Hekele J, Yao Y, Kanai Y, Blum V, Kratzer P. All-electron real-time and imaginary-time time-dependent density functional theory within a numeric atom-centered basis function framework. J Chem Phys 2021; 155:154801. [PMID: 34686041 DOI: 10.1063/5.0066753] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Real-time time-dependent density functional theory (RT-TDDFT) is an attractive tool to model quantum dynamics by real-time propagation without the linear response approximation. Sharing the same technical framework of RT-TDDFT, imaginary-time time-dependent density functional theory (it-TDDFT) is a recently developed robust-convergence ground state method. Presented here are high-precision all-electron RT-TDDFT and it-TDDFT implementations within a numerical atom-centered orbital (NAO) basis function framework in the FHI-aims code. We discuss the theoretical background and technical choices in our implementation. First, RT-TDDFT results are validated against linear-response TDDFT results. Specifically, we analyze the NAO basis sets' convergence for Thiel's test set of small molecules and confirm the importance of the augmentation basis functions for adequate convergence. Adopting a velocity-gauge formalism, we next demonstrate applications for systems with periodic boundary conditions. Taking advantage of the all-electron full-potential implementation, we present applications for core level spectra. For it-TDDFT, we confirm that within the all-electron NAO formalism, it-TDDFT can successfully converge systems that are difficult to converge in the standard self-consistent field method. We finally benchmark our implementation for systems up to ∼500 atoms. The implementation exhibits almost linear weak and strong scaling behavior.
Collapse
Affiliation(s)
- Joscha Hekele
- Faculty of Physics, University of Duisburg-Essen, Lotharstr. 1, 47057 Duisburg, Germany
| | - Yi Yao
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | - Volker Blum
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - Peter Kratzer
- Faculty of Physics, University of Duisburg-Essen, Lotharstr. 1, 47057 Duisburg, Germany
| |
Collapse
|
39
|
Khalili F, Vafaee M, Shokri B. Attosecond charge migration following oxygen K-shell ionization in DNA bases and base pairs. Phys Chem Chem Phys 2021; 23:23005-23013. [PMID: 34611693 DOI: 10.1039/d1cp02920g] [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
Core ionization of DNA begins a cascade of events which could lead to cellular inactivation or death. The created core-hole following an impulse inner-shell ionization of molecules naturally decays in the auger timescale. We simulated charge migration (CM) phenomena following an impulsive core ionization of individual DNA bases at the oxygen K-edge which occurs before Auger decay of the oxygen. Our approach is based on real-time time dependent density functional theory (RT-TDDFT). It is shown that the pronounced hole fluctuation observed around bonds of the initial core-hole results in various valence orbital migrations. Also, the same photo-core-ionized dynamics is studied for the related base pairs. We investigate the role of base pairing and H-bonding interactions in the attosecond CM dynamics. In particular, the creation of a core-hole in the oxygen involved in H-bonding leads to an enhancement of charge migration relative to the respective single bases. Importantly, the hole oscillation of the adenine-thymine base pair upon creation of a core-hole at the oxygen, which does not contribute to the donor-acceptor interactions (not H-bonded), decreases compared to the single thymine base. Understanding the detailed dynamics of the localized core-hole initiating CM process would open the way for chemically controlling DNA damage/repair in the future.
Collapse
Affiliation(s)
- Fatemeh Khalili
- Department of Physics, Shahid Beheshti University, Velenjak, Tehran 19839, Iran.
| | - Mohsen Vafaee
- Department of Chemistry, Tarbiat Modares University, P. O. Box 14115-175, Tehran, Iran.
| | - Babak Shokri
- Department of Physics, Shahid Beheshti University, Velenjak, Tehran 19839, Iran. .,Laser-Plasma Research Institute, Shahid Beheshti University, Velenjak, Tehran 19839, Iran
| |
Collapse
|
40
|
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
|
41
|
Schreder L, Luber S. Local approaches for electric dipole moments in periodic systems and their application to real-time time-dependent density functional theory. J Chem Phys 2021; 155:134116. [PMID: 34624999 DOI: 10.1063/5.0058793] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Within periodic boundary conditions, the traditional quantum mechanical position operator is ill-defined, necessitating the use of alternative methods, most commonly the Berry phase formulation in the modern theory of polarization. Since any information about local properties is lost in this change of framework, the Berry phase formulation can only determine the total electric polarization of a system. Previous approaches toward recovering local electric dipole moments have been based on applying the conventional dipole moment operator to the centers of maximally localized Wannier functions (MLWFs). Recently, another approach to local electric dipole moments has been demonstrated in the field of subsystem density functional theory (DFT) embedding. We demonstrate in this work that this approach, aside from its use in ground state DFT-based molecular dynamics, can also be applied to obtain electric dipole moments during real-time propagated time-dependent DFT (RT-TDDFT). Moreover, we present an analogous approach to obtain local electric dipole moments from MLWFs, which enables subsystem analysis in cases where DFT embedding is not applicable. The techniques were implemented in the quantum chemistry software CP2K for the mixed Gaussian and plane wave method and applied to cis-diimide and water in the gas phase, cis-diimide in aqueous solution, and a liquid mixture of dimethyl carbonate and ethylene carbonate to obtain absorption and infrared spectra decomposed into localized subsystem contributions.
Collapse
Affiliation(s)
- Lukas Schreder
- University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Sandra Luber
- University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| |
Collapse
|
42
|
Shepard C, Zhou R, Yost DC, Yao Y, Kanai Y. Simulating electronic excitation and dynamics with real-time propagation approach to TDDFT within plane-wave pseudopotential formulation. J Chem Phys 2021; 155:100901. [PMID: 34525811 DOI: 10.1063/5.0057587] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We give a perspective on simulating electronic excitation and dynamics using the real-time propagation approach to time-dependent density functional theory (RT-TDDFT) in the plane-wave pseudopotential formulation. RT-TDDFT is implemented in various numerical formalisms in recent years, and its practical application often dictates the most appropriate implementation of the theory. We discuss recent developments and challenges, emphasizing numerical aspects of studying real systems. Several applications of RT-TDDFT simulation are discussed to highlight how the approach is used to study interesting electronic excitation and dynamics phenomena in recent years.
Collapse
Affiliation(s)
- Christopher Shepard
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| | - Ruiyi Zhou
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| | - Dillon C Yost
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02142, USA
| | - Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| |
Collapse
|
43
|
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]
|
44
|
Sun J, Ding Z, Yu Y, Liang W. Nonlinear features of Fano resonance: a QM/EM study. Phys Chem Chem Phys 2021; 23:15994-16004. [PMID: 34318831 DOI: 10.1039/d1cp02459k] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The nonlinear Fano effects on the absorption of hybrid systems composed of a silver nanosphere and an indoline dye molecule have been systematically investigated by the hybrid approach, which combines the quantum mechanics method (QM) with the computational electromagnetic method (EM). The absorption spectra of the dye molecule in the proximity of an Ag nanoparticle have been calculated by changing the incident field intensity, the phenomenological dephasing of molecular excitation, and the enhancement ratio of the near field. The contribution of molecular nonlinear response properties and the quantum interferences of the incident and scattered fields and of resonant plasmon-molecular excitations to the spectra has been identified. It is in no doubt that Fano resonance due to the plasmon-molecular interaction can appear in both the weak and strong field regimes; however, the Fano effect is more pronounced in the strong field regime where quantum interference leads to a nonlinear Fano effect controlled by a complex field-dependent Fano factor. When the incident field is strong enough, the resonance antisymmetry structure is spectrally resolved, and it changes with the change of the field intensity. As the field intensity varies from weak to strong, the Fano lineshape's asymmetry increases with increasing intensity in the beginning, and then decreases with a further increase of the field intensity attributed to the increase of the detuning energy induced by the integrated energy shift upon field dressing during the excitation. Decreasing the enhancement ratio of the near field or the dephasing of molecular excitation can also control the spectral lineshape transformation from an asymmetric profile to a symmetric Lorentzian lineshape. These findings are consistent with previous experimental and theoretical observations arisen by quantum interferences and are expected to stimulate further work toward exploring the plasmon-molecular interplay and the applications of Fano resonance in optical switching and sensing.
Collapse
Affiliation(s)
- Jin Sun
- School of Physics and Optoelectronics Engineering, Energy Materials and Devices Key Lab of Anhui Province for Photoelectric Conversion, Anhui University, Hefei 230601, People's Republic of China
| | | | | | | |
Collapse
|
45
|
Baiardi A. Electron Dynamics with the Time-Dependent Density Matrix Renormalization Group. J Chem Theory Comput 2021; 17:3320-3334. [PMID: 34043347 DOI: 10.1021/acs.jctc.0c01048] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
In this work, we simulate the electron dynamics in molecular systems with the time-dependent density matrix renormalization group (TD-DMRG) algorithm. We leverage the generality of the so-called tangent-space TD-DMRG formulation and design a computational framework in which the dynamics is driven by the exact nonrelativistic electronic Hamiltonian. We show that by parametrizing the wave function as a matrix product state, we can accurately simulate the dynamics of systems including up to 20 electrons and 32 orbitals. We apply the TD-DMRG algorithm to three problems that are hardly targeted by time-independent methods: the calculation of molecular (hyper)polarizabilities, the simulation of electronic absorption spectra, and the study of ultrafast ionization dynamics.
Collapse
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, Zürich 8093, Switzerland
| |
Collapse
|
46
|
Moreno Carrascosa A, Yang M, Yong H, Ma L, Kirrander A, Weber PM, Lopata K. Mapping static core-holes and ring-currents with X-ray scattering. Faraday Discuss 2021; 228:60-81. [PMID: 33605956 DOI: 10.1039/d0fd00124d] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Measuring the attosecond movement of electrons in molecules is challenging due to the high temporal and spatial resolutions required. X-ray scattering-based methods are promising, but many questions remain concerning the sensitivity of the scattering signals to changes in density, as well as the means of reconstructing the dynamics from these signals. In this paper, we present simulations of stationary core-holes and electron dynamics following inner-shell ionization of the oxazole molecule. Using a combination of time-dependent density functional theory simulations along with X-ray scattering theory, we demonstrate that the sudden core-hole ionization produces a significant change in the X-ray scattering response and how the electron currents across the molecule should manifest as measurable modulations to the time dependent X-ray scattering signal. This suggests that X-ray scattering is a viable probe for measuring electronic processes at time scales faster than nuclear motion.
Collapse
Affiliation(s)
| | - Mengqi Yang
- Department of Chemistry, 232 Choppin Hall, Baton Rouge, Louisiana 70803, USA
| | - Haiwang Yong
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Lingyu Ma
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Adam Kirrander
- EaStCHEM, School of Chemistry, University of Edinburgh, David Brewster Road, EH9 3FJ Edinburgh, UK
| | - Peter M Weber
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Kenneth Lopata
- Department of Chemistry, 232 Choppin Hall, Baton Rouge, Louisiana 70803, USA and Center for Computation and Technology, Louisiana State University, Baton Roug, Louisiana 70803, USA.
| |
Collapse
|
47
|
Mantela M, Morphis A, Lambropoulos K, Simserides C, Di Felice R. Effects of Structural Dynamics on Charge Carrier Transfer in B-DNA: A Combined MD and RT-TDDFT Study. J Phys Chem B 2021; 125:3986-4003. [PMID: 33857373 DOI: 10.1021/acs.jpcb.0c11489] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Hole transfer along the axis of duplex DNA has been the focus of physical chemistry research for decades, with implications in diverse fields, from nanotechnology to cell oxidative damage. Computational approaches are particularly amenable for this problem, to complement experimental data for interpretation of transfer mechanisms. To be predictive, computational results need to account for the inherent mobility of biological molecules during the time frame of experimental measurements. Here, we address the structural variability of B-DNA and its effects on hole transfer in a combined molecular dynamics (MD) and real-time time-dependent density functional theory (RT-TDDFT) study. Our results show that quantities that characterize the charge transfer process, such as the time-dependent dipole moment and hole population at a specific site, are sensitive to structural changes that occur on the nanosecond time scale. We extend the range of physical properties for which such a correlation has been observed, further establishing the fact that quantitative computational data on charge transfer properties should include statistical averages. Furthermore, we use the RT-TDDFT results to assess an efficient tight-binding method suitable for high-throughput predictions. We demonstrate that charge transfer, although affected by structural variability, on average, remains strong in AA and GG dimers.
Collapse
Affiliation(s)
- Marilena Mantela
- Department of Physics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos GR-15784, Athens, Greece
| | - Andreas Morphis
- Department of Physics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos GR-15784, Athens, Greece
| | - Konstantinos Lambropoulos
- Department of Physics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos GR-15784, Athens, Greece
| | - Constantinos Simserides
- Department of Physics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos GR-15784, Athens, Greece
| | | |
Collapse
|
48
|
Wibowo M, Irons TJP, Teale AM. Modeling Ultrafast Electron Dynamics in Strong Magnetic Fields Using Real-Time Time-Dependent Electronic Structure Methods. J Chem Theory Comput 2021; 17:2137-2165. [PMID: 33724806 PMCID: PMC8047917 DOI: 10.1021/acs.jctc.0c01269] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
![]()
An implementation
of real-time time-dependent Hartree–Fock
(RT-TDHF) and current density functional theory (RT-TDCDFT) for molecules
in strong uniform magnetic fields is presented. In contrast to earlier
implementations, the present work enables the use of the RT-TDCDFT
formalism, which explicitly includes field-dependent terms in the
exchange–correlation functional. A range of current-dependent
exchange–correlation functionals based on the TPSS functional
are considered, including a range-separated variant, which is particularly
suitable for application to excited state calculations. The performance
of a wide range of propagator algorithms for real-time methods is
investigated in this context. A recently proposed molecular orbital
pair decomposition analysis allows for assignment of electronic transitions,
providing detailed information about which molecular orbitals are
involved in each excitation. The application of these methods is demonstrated
for the electronic absorption spectra of N2 and H2O both in the absence and in the presence of a magnetic field. The
dependence of electronic spectra on the magnetic field strength and
its orientation relative to the molecule is studied. The complex evolution
of the absorption spectra with magnetic field is rationalized using
the molecular orbital pair decomposition analysis, which provides
crucial insight in strong fields where the spectra are radically different
from their zero-field counterparts.
Collapse
Affiliation(s)
- Meilani Wibowo
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - Tom J P Irons
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - Andrew M Teale
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom.,Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, N-0315 Oslo, Norway
| |
Collapse
|
49
|
Folorunso AS, Bruner A, Mauger F, Hamer KA, Hernandez S, Jones RR, DiMauro LF, Gaarde MB, Schafer KJ, Lopata K. Molecular Modes of Attosecond Charge Migration. PHYSICAL REVIEW LETTERS 2021; 126:133002. [PMID: 33861123 DOI: 10.1103/physrevlett.126.133002] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Accepted: 02/24/2021] [Indexed: 06/12/2023]
Abstract
First-principles calculations are employed to elucidate the modes of attosecond charge migration (CM) in halogenated hydrocarbon chains. We use constrained density functional theory (DFT) to emulate the creation of a localized hole on the halogen and follow the subsequent dynamics via time-dependent DFT. We find low-frequency CM modes (∼1 eV) that propagate across the molecule and study their dependence on length, bond order, and halogenation. We observe that the CM speed (∼4 Å/fs) is largely independent of molecule length, but is lower for triple-bonded versus double-bonded molecules. Additionally, as the halogen mass increases, the hole travels in a more particlelike manner as it moves across the molecule. These heuristics will be useful in identifying molecules and optimal CM detection methods for future experiments, especially for halogenated hydrocarbons which are promising targets for ionization-triggered CM.
Collapse
Affiliation(s)
- Aderonke S Folorunso
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Adam Bruner
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - François Mauger
- Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Kyle A Hamer
- Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Samuel Hernandez
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Robert R Jones
- Department of Physics, University of Virginia, Charlottesville, Virginia 22904, USA
| | - Louis F DiMauro
- Department of Physics, The Ohio State University, Columbus, Ohio 43210, USA
| | - Mette B Gaarde
- Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Kenneth J Schafer
- Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
- Center for Computation and Technology, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| |
Collapse
|
50
|
Woźniak AP, Lesiuk M, Przybytek M, Efimov DK, Prauzner-Bechcicki JS, Mandrysz M, Ciappina M, Pisanty E, Zakrzewski J, Lewenstein M, Moszyński R. A systematic construction of Gaussian basis sets for the description of laser field ionization and high-harmonic generation. J Chem Phys 2021; 154:094111. [PMID: 33685145 DOI: 10.1063/5.0040879] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A precise understanding of mechanisms governing the dynamics of electrons in atoms and molecules subjected to intense laser fields has a key importance for the description of attosecond processes such as the high-harmonic generation and ionization. From the theoretical point of view, this is still a challenging task, as new approaches to solve the time-dependent Schrödinger equation with both good accuracy and efficiency are still emerging. Until recently, the purely numerical methods of real-time propagation of the wavefunction using finite grids have been frequently and successfully used to capture the electron dynamics in small one- or two-electron systems. However, as the main focus of attoscience shifts toward many-electron systems, such techniques are no longer effective and need to be replaced by more approximate but computationally efficient ones. In this paper, we explore the increasingly popular method of expanding the wavefunction of the examined system into a linear combination of atomic orbitals and present a novel systematic scheme for constructing an optimal Gaussian basis set suitable for the description of excited and continuum atomic or molecular states. We analyze the performance of the proposed basis sets by carrying out a series of time-dependent configuration interaction calculations for the hydrogen atom in fields of intensity varying from 5 × 1013 W/cm2 to 5 × 1014 W/cm2. We also compare the results with the data obtained using Gaussian basis sets proposed previously by other authors.
Collapse
Affiliation(s)
| | - Michał Lesiuk
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Michał Przybytek
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Dmitry K Efimov
- Institute of Theoretical Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Jakub S Prauzner-Bechcicki
- Marian Smoluchowski Institute of Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Michał Mandrysz
- Institute of Theoretical Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Marcelo Ciappina
- ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, Av. Carl Friedrich Gauss 3, 08860, Castelldefels, Barcelona, Spain
| | - Emilio Pisanty
- ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, Av. Carl Friedrich Gauss 3, 08860, Castelldefels, Barcelona, Spain
| | - Jakub Zakrzewski
- Institute of Theoretical Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Maciej Lewenstein
- ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, Av. Carl Friedrich Gauss 3, 08860, Castelldefels, Barcelona, Spain
| | - Robert Moszyński
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|