1
|
Herring C, Montemore MM. Recent Advances in Real-Time Time-Dependent Density Functional Theory Simulations of Plasmonic Nanostructures and Plasmonic Photocatalysis. ACS NANOSCIENCE AU 2023; 3:269-279. [PMID: 37601917 PMCID: PMC10436373 DOI: 10.1021/acsnanoscienceau.2c00061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/12/2022] [Revised: 05/02/2023] [Accepted: 05/08/2023] [Indexed: 08/22/2023]
Abstract
Plasmonic catalysis provides a possible means for driving chemical reactions under relatively mild conditions. Rational design of these systems is impeded by the difficulty in understanding the electron dynamics and their interplay with reactions. Real-time, time-dependent density functional theory (RT-TDDFT) can provide dynamic information on excited states in plasmonic systems, including those relevant to plasmonic catalysis, at time scales and length scales that are otherwise out of reach of many experimental techniques. Here, we discuss previous RT-TDDFT studies of plasmonic systems, focusing on recent work that gains insight into plasmonic catalysis. These studies provide insight into plasmon dynamics, including size effects and the role of specific electronic states. Further, these studies provide significant insight into mechanisms underlying plasmonic catalysis, showing the importance of charge transfer between metal and adsorbate states, as well as local field enhancement, in different systems.
Collapse
Affiliation(s)
- Connor
J. Herring
- Department of Chemical and Biomolecular
Engineering, Tulane University, New Orleans, Louisiana 70115, United States
| | - Matthew M. Montemore
- Department of Chemical and Biomolecular
Engineering, Tulane University, New Orleans, Louisiana 70115, United States
| |
Collapse
|
2
|
Yuwono SH, Cooper BC, Zhang T, Li X, DePrince AE. Time-dependent equation-of-motion coupled-cluster simulations with a defective Hamiltonian. J Chem Phys 2023; 159:044113. [PMID: 37497820 DOI: 10.1063/5.0157852] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2023] [Accepted: 07/11/2023] [Indexed: 07/28/2023] Open
Abstract
Simulations of laser-induced electron dynamics in a molecular system are performed using time-dependent (TD) equation-of-motion (EOM) coupled-cluster (CC) theory. The target system has been chosen to highlight potential shortcomings of truncated TD-EOM-CC methods [represented in this work by TD-EOM-CC with single and double excitations (TD-EOM-CCSD)], where unphysical spectroscopic features can emerge. Specifically, we explore driven resonant electronic excitations in magnesium fluoride in the proximity of an avoided crossing. Near the avoided crossing, the CCSD similarity-transformed Hamiltonian is defective, meaning that it has complex eigenvalues, and oscillator strengths may take on negative values. When an external field is applied to drive transitions to states exhibiting these traits, unphysical dynamics are observed. For example, the stationary states that make up the time-dependent state acquire populations that can be negative, exceed one, or even complex-valued.
Collapse
Affiliation(s)
- Stephen H Yuwono
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| | - Brandon C Cooper
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| | - Tianyuan Zhang
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
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
|
6
|
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
|
7
|
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
|
8
|
Yang M, Sissay A, Chen M, Lopata K. Intruder Peak-Free Transient Inner-Shell Spectra Using Real-Time Simulations. J Chem Theory Comput 2022; 18:992-1002. [PMID: 35025498 DOI: 10.1021/acs.jctc.1c00079] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Real-time methods are convenient for simulating core-level absorption spectra but suffer from nonphysical intruder peaks when using atom-centered basis sets. In transient absorption spectra, these peaks exhibit highly nonphysical time-dependent modulations in their energies and oscillator strengths. In this paper, we address the origins of these intruder peaks and propose a straightforward and effective solution based on a filtered dipole operator. In combination with real-time time-dependent density functional theory (RT-TDDFT), we demonstrate how to compute intruder-free attosecond transient X-ray absorption spectra for the aminophenol (C6H7NO) oxygen and nitrogen K-edges and the α-quartz (SiO2) silicon L-edge. Without filtering, the computed spectra are qualitatively wrong. This procedure is suitable for both static and transient inner-shell spectroscopy studies and can easily be implemented in a range of real-time methodologies.
Collapse
Affiliation(s)
- Mengqi Yang
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Adonay Sissay
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Min Chen
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States.,Center for Computation and Technology, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
9
|
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
|
10
|
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
|
11
|
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
|
12
|
Levine BG, Durden AS, Esch MP, Liang F, Shu Y. CAS without SCF-Why to use CASCI and where to get the orbitals. J Chem Phys 2021; 154:090902. [PMID: 33685182 DOI: 10.1063/5.0042147] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree-Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Finally, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.
Collapse
Affiliation(s)
- Benjamin G Levine
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Andrew S Durden
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Michael P Esch
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Fangchun Liang
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
13
|
Mosquera MA, Jones LO, Kang G, Ratner MA, Schatz GC. Second Linear Response Theory and the Analytic Calculation of Excited-State Properties. J Phys Chem A 2021; 125:1093-1102. [PMID: 33497573 DOI: 10.1021/acs.jpca.0c10152] [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/28/2022]
Abstract
We present a method based on second linear response time-dependent density functional theory (TDDFT) to calculate permanent and transition multipoles of excited states, which are required to compute excited-state absorption/emission spectra and multiphoton optical processes, among others. In previous work, we examined computations based on second linear response theory in which linear response TDDFT was employed twice. In contrast, the present methodology requires information from only a single linear response calculation to compute the excited-state properties. These are evaluated analytically through various algebraic operations involving electron repulsion integrals and excitation vectors. The present derivation focuses on full many-body wave functions instead of single orbitals, as in our previous approach. We test the proposed method by applying it to several diatomic and triatomic molecules. This shows that the computed excited-state dipoles are consistent with respect to reference equation-of-motion coupled-cluster calculations.
Collapse
Affiliation(s)
- Martín A Mosquera
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Department of Chemistry and Biochemistry, Montana State University, Bozeman, Montana 59717, United States
| | - Leighton O Jones
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Gyeongwon Kang
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Mark A Ratner
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - George C Schatz
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| |
Collapse
|
14
|
Müller C, Sharma M, Sierka M. Real-time time-dependent density functional theory using density fitting and the continuous fast multipole method. J Comput Chem 2020; 41:2573-2582. [PMID: 33464600 DOI: 10.1002/jcc.26412] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Revised: 08/17/2020] [Accepted: 08/18/2020] [Indexed: 11/07/2022]
Abstract
An implementation of real-time time-dependent density functional theory (RT-TDDFT) within the TURBOMOLE program package is reported using Gaussian-type orbitals as basis functions, second and fourth order Magnus propagator, and the self-consistent field as well as the predictor-corrector time integration schemes. The Coulomb contribution to the Kohn-Sham matrix is calculated combining density fitting approximation and the continuous fast multipole method. Performance of the implementation is benchmarked for molecular systems with different sizes and dimensionalities. For linear alkane chains, the wall time for density matrix time propagation step is comparable to the Kohn-Sham (KS) matrix construction. However, for larger two- and three-dimensional molecules, with up to about 5,000 basis functions, the computational effort of RT-TDDFT calculations is dominated by the KS matrix evaluation. In addition, the maximum time step is evaluated using a set of small molecules of different polarities. The photoabsorption spectra of several molecular systems calculated using RT-TDDFT are compared to those obtained using linear response time-dependent density functional theory and coupled cluster methods.
Collapse
Affiliation(s)
- Carolin Müller
- Otto Schott Institute of Materials Research, Friedrich Schiller University of Jena, Jena, Germany
| | - Manas Sharma
- Otto Schott Institute of Materials Research, Friedrich Schiller University of Jena, Jena, Germany
| | - Marek Sierka
- Otto Schott Institute of Materials Research, Friedrich Schiller University of Jena, Jena, Germany
| |
Collapse
|
15
|
Li X, Govind N, Isborn C, DePrince AE, Lopata K. Real-Time Time-Dependent Electronic Structure Theory. Chem Rev 2020; 120:9951-9993. [DOI: 10.1021/acs.chemrev.0c00223] [Citation(s) in RCA: 64] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christine Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
16
|
Krumland J, Valencia AM, Pittalis S, Rozzi CA, Cocchi C. Understanding real-time time-dependent density-functional theory simulations of ultrafast laser-induced dynamics in organic molecules. J Chem Phys 2020; 153:054106. [PMID: 32770886 DOI: 10.1063/5.0008194] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Real-time time-dependent density functional theory, in conjunction with the Ehrenfest molecular dynamics scheme, is becoming a popular methodology to investigate ultrafast phenomena on the nanoscale. Thanks to recent developments, it is also possible to explicitly include in the simulations a time-dependent laser pulse, thereby accessing the transient excitation regime. However, the complexity entailed in these calculations calls for in-depth analysis of the accessible and yet approximate (either "dressed" or "bare") quantities in order to evaluate their ability to provide us with a realistic picture of the simulated processes. In this work, we analyze the ultrafast dynamics of three small molecules (ethylene, benzene, and thiophene) excited by a resonant laser pulse in the framework of the adiabatic local-density approximation. The electronic response to the laser perturbation in terms of induced dipole moment and excited-state population is compared to the results given by an exactly solvable two-level model. In this way, we can interpret the charge-carrier dynamics in terms of simple estimators, such as the number of excited electrons. From the computed transient absorption spectra, we unravel the appearance of nonlinear effects such as excited-state absorption and vibronic coupling. In this way, we observe that the laser excitation affects the vibrational spectrum by enhancing the anharmonicities therein, while the coherent vibrational motion contributes to stabilizing the electronic excitation already within a few tens of femtoseconds.
Collapse
Affiliation(s)
- Jannis Krumland
- Humboldt-Universität zu Berlin, Physics Department and IRIS Adlershof, 12489 Berlin, Germany
| | - Ana M Valencia
- Humboldt-Universität zu Berlin, Physics Department and IRIS Adlershof, 12489 Berlin, Germany
| | | | | | - Caterina Cocchi
- Humboldt-Universität zu Berlin, Physics Department and IRIS Adlershof, 12489 Berlin, Germany
| |
Collapse
|
17
|
Chen M, Lopata K. First-Principles Simulations of X-ray Transient Absorption for Probing Attosecond Electron Dynamics. J Chem Theory Comput 2020; 16:4470-4478. [PMID: 32470295 PMCID: PMC7467644 DOI: 10.1021/acs.jctc.0c00122] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
X-ray transient absorption spectroscopy (XTAS) is a promising technique for measuring electron dynamics in molecules and solids with attosecond time resolutions. In XTAS, the elemental specificity and spatial locality of core-to-valence X-ray absorption is exploited to relate modulations in the time-resolved absorption spectra to local electron density variations around particular atoms. However, interpreting these absorption modulations and frequency shifts as a function of the time delay in terms of dynamics can be challenging. In this paper, we present a first-principles study of attosecond XTAS in a selection of simple molecules based on real-time time-dependent density functional theory (RT-TDDFT) with constrained DFT to emulate the state of the system following the interaction with a ultraviolet pump laser. In general, there is a decrease in the optical density and a blue shift in the frequency with increasing electron density around the absorbing atom. In carbon monoxide (CO), modulations in the O K-edge occur at the frequency of the valence electron dynamics, while for dioxygen (O2) they occur at twice the frequency, due to the indistinguishability of the oxygen atoms. In 4-aminophenol (H2NC6H4OH), likewise, there is a decrease in the optical density and a blue shift in the frequency for the oxygen and nitrogen K-edges with increasing charge density on the O and N, respectively. Similar effects are observed in the nitrogen K-edge for a long-range charge-transfer excitation in a benzene (C6H6)-tetracyanoethylene (C6N4; TCNE) dimer but with weaker modulations due to the delocalization of the charge across the entire TCNE molecule. Additionally, in all cases, there are pre-edge features corresponding to core transitions to depopulated orbitals. These potentially offer a background-free signal that only appears in pumped molecules.
Collapse
Affiliation(s)
- Min Chen
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States.,Center for Computation and Technology, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
18
|
Bonafé FP, Aradi B, Hourahine B, Medrano CR, Hernández FJ, Frauenheim T, Sánchez CG. A Real-Time Time-Dependent Density Functional Tight-Binding Implementation for Semiclassical Excited State Electron–Nuclear Dynamics and Pump–Probe Spectroscopy Simulations. J Chem Theory Comput 2020; 16:4454-4469. [DOI: 10.1021/acs.jctc.9b01217] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Franco P. Bonafé
- Max Planck Institute for the Structure and Dynamics of Matter, Hamburg, Germany
- Universidad Nacional de Córdoba, Facultad de Ciencias Químicas, Departamento de Química Teórica y Computacional, Córdoba, Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba, INFIQC (CONICET - Universidad Nacional de Córdoba), Córdoba, Argentina
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, Universitát Bremen, Bremen, Germany
| | - Ben Hourahine
- SUPA, Department of Physics, John Anderson Building, The University of Strathclyde, 107 Rottenrow, Glasgow G15 6QN, United Kingdom
| | - Carlos R. Medrano
- Universidad Nacional de Córdoba, Facultad de Ciencias Químicas, Departamento de Química Teórica y Computacional, Córdoba, Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba, INFIQC (CONICET - Universidad Nacional de Córdoba), Córdoba, Argentina
| | - Federico J. Hernández
- Universidad Nacional de Córdoba, Facultad de Ciencias Químicas, Departamento de Química Teórica y Computacional, Córdoba, Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba, INFIQC (CONICET - Universidad Nacional de Córdoba), Córdoba, Argentina
- Department of Physics, Universidad de Santiago de Chile, Av. Ecuador 3493, Santiago, Chile
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, Universitát Bremen, Bremen, Germany
- Computational Science Research Center (CSRC) Beijing and Computational Science and Applied Research (CSAR) Institute, Shenzhen, China
| | - Cristián G. Sánchez
- Instituto Interdisciplinario de Ciencias Básicas, Universidad Nacional de Cuyo, CONICET, Facultad de Ciencias Exactas y Naturales, Mendoza, Argentina
| |
Collapse
|
19
|
Usukura J, Hiyama M, Kurata M, Hazama Y, Qiu XP, Winnik FM, Akiyama H, Koga N. Theoretical Study of the Wavelength Selection for the Photocleavage of Coumarin-caged D-luciferin. Photochem Photobiol 2020; 96:805-814. [PMID: 31907932 DOI: 10.1111/php.13212] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2019] [Accepted: 12/31/2019] [Indexed: 11/29/2022]
Abstract
The equilibrium structures and optical properties of the photolabile caged luciferin, (7-diethylaminocoumarin-4-yl)methyl caged D-luciferin (DEACM-caged D-luciferin), in aqueous solution were investigated via quantum chemical calculations. The probable conformers of DEACM-caged D-luciferin were determined by potential energy curve scans and structural optimizations. We identified 40 possible conformers of DEACM-caged D-luciferin in water by comparing the Gibbs free energy of the optimized structures. Despite the difference in their structures, the conformers were similar in terms of assignments, oscillator strengths and energies of the three low-lying excited states. From the concentrations of the conformers and their oscillator strengths, we obtained a theoretical UV/Vis spectrum of DEACM-caged D-luciferin that has two main bands of shape nearly identical to the experimental UV/Vis spectrum. The absorption bands with maxima ~ 384 and 339 nm were attributed to the electronic excitations of the caged group and the luciferin moiety, respectively, by analysis of the theoretical UV/Vis spectrum. Furthermore, the analysis showed that DEACM-caged D-luciferin is excited in the caged group only by light of wavelength ranging within 400-430 nm, which is in the long-wavelength tail of the 384 nm band. This should be tested to lower damage upon photocleavage.
Collapse
Affiliation(s)
- Junko Usukura
- Institute for Solid Physics, University of Tokyo, Kashiwa, Japan
| | - Miyabi Hiyama
- Graduate School of Science and Technology, Gunma University, Kiryu, Japan
| | - Maki Kurata
- Institute for Solid Physics, University of Tokyo, Kashiwa, Japan
| | - Yuji Hazama
- Institute for Solid Physics, University of Tokyo, Kashiwa, Japan.,AIST-UTokyo Advanced Operando-Measurement Technology Open Innovation Laboratory (OPERANDO-OIL), Kashiwa, Japan
| | - Xing-Ping Qiu
- Département de Chimie, Université de Montréal, Montréal, Canada
| | - Francoise M Winnik
- Department of Chemistry, University of Helsinki, Helsinki, Finland.,International Center for Materials Nanoarchitectonics (MANA), National Institute for Materials Science (NIMS), Tsukuba, Japan
| | - Hidefumi Akiyama
- Institute for Solid Physics, University of Tokyo, Kashiwa, Japan.,AIST-UTokyo Advanced Operando-Measurement Technology Open Innovation Laboratory (OPERANDO-OIL), Kashiwa, Japan
| | - Nobuaki Koga
- Graduate School of Informatics, Nagoya University, Nagoya, Japan
| |
Collapse
|
20
|
Nenov A, Segatta F, Bruner A, Mukamel S, Garavelli M. X-ray linear and non-linear spectroscopy of the ESCA molecule. J Chem Phys 2019; 151:114110. [DOI: 10.1063/1.5116699] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Artur Nenov
- Dipartimento di Chimica Industriale “Toso Montanari”, Università degli studi di Bologna, Viale del Risorgimento 4,
40136 Bologna, Italy
| | - Francesco Segatta
- Dipartimento di Chimica Industriale “Toso Montanari”, Università degli studi di Bologna, Viale del Risorgimento 4,
40136 Bologna, Italy
| | - Adam Bruner
- Department of Chemistry and Physics and Astronomy, University of California, Irvine, California 92697,
USA
| | - Shaul Mukamel
- Department of Chemistry and Physics and Astronomy, University of California, Irvine, California 92697,
USA
| | - Marco Garavelli
- Dipartimento di Chimica Industriale “Toso Montanari”, Università degli studi di Bologna, Viale del Risorgimento 4,
40136 Bologna, Italy
| |
Collapse
|
21
|
Ghosh S, Asher JC, Gagliardi L, Cramer CJ, Govind N. A semiempirical effective Hamiltonian based approach for analyzing excited state wave functions and computing excited state absorption spectra using real-time dynamics. J Chem Phys 2019; 150:104103. [DOI: 10.1063/1.5061746] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Soumen Ghosh
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, USA
| | - Jason C. Asher
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, USA
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, USA
| | - Christopher J. Cramer
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, USA
| | - Niranjan Govind
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99338, USA
| |
Collapse
|
22
|
Hernández FJ, Bonafé FP, Aradi B, Frauenheim T, Sánchez CG. Simulation of Impulsive Vibrational Spectroscopy. J Phys Chem A 2019; 123:2065-2072. [DOI: 10.1021/acs.jpca.9b00307] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Federico J. Hernández
- Universidad Nacional de Córdoba. Facultad de Ciencias Quı́micas, Departamento de Quı́mica Teórica y Computacional, Córdoba Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba, INFIQC (CONICET-Universidad Nacional de Córdoba), Córdoba 5000, Argentina
| | - Franco P. Bonafé
- Universidad Nacional de Córdoba. Facultad de Ciencias Quı́micas, Departamento de Quı́mica Teórica y Computacional, Córdoba Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba, INFIQC (CONICET-Universidad Nacional de Córdoba), Córdoba 5000, Argentina
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, Universität Bremen, Bremen 28359, Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, Universität Bremen, Bremen 28359, Germany
| | - Cristián G. Sánchez
- Universidad Nacional de Córdoba. Facultad de Ciencias Quı́micas, Departamento de Quı́mica Teórica y Computacional, Córdoba Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba, INFIQC (CONICET-Universidad Nacional de Córdoba), Córdoba 5000, Argentina
| |
Collapse
|
23
|
Liu H, Jenkins AJ, Wildman A, Frisch MJ, Lipparini F, Mennucci B, Li X. Time-Dependent Complete Active Space Embedded in a Polarizable Force Field. J Chem Theory Comput 2019; 15:1633-1641. [DOI: 10.1021/acs.jctc.8b01152] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Hongbin Liu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Andrew J. Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Andrew Wildman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Michael J. Frisch
- Gaussian Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, 56126 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, 56126 Pisa, Italy
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
24
|
Comparison and convergence of optical absorption spectra of noble metal nanoparticles computed using linear-response and real-time time-dependent density functional theories. COMPUT THEOR CHEM 2018. [DOI: 10.1016/j.comptc.2018.11.005] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
25
|
Mattiat J, Luber S. Efficient calculation of (resonance) Raman spectra and excitation profiles with real-time propagation. J Chem Phys 2018; 149:174108. [DOI: 10.1063/1.5051250] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Johann Mattiat
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, Zurich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, Zurich, Switzerland
| |
Collapse
|
26
|
Lian C, Hu SQ, Guan MX, Meng S. Momentum-resolved TDDFT algorithm in atomic basis for real time tracking of electronic excitation. J Chem Phys 2018; 149:154104. [PMID: 30342439 DOI: 10.1063/1.5036543] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ultrafast electronic dynamics in solids lies at the core of modern condensed matter and materials physics. To build up a practical ab initio method for studying solids under photoexcitation, we develop a momentum-resolved real-time time dependent density functional theory (rt-TDDFT) algorithm using numerical atomic basis, together with the implementation of both the length and vector gauge of the electromagnetic field. When applied to simulate elementary excitations in two-dimensional materials such as graphene, different excitation modes, only distinguishable in momentum space, are observed. The momentum-resolved rt-TDDFT is important and computationally efficient for the study of ultrafast dynamics in extended systems.
Collapse
Affiliation(s)
- Chao Lian
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
| | - Shi-Qi Hu
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
| | - Meng-Xue Guan
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
| | - Sheng Meng
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
| |
Collapse
|
27
|
Peng WT, Fales BS, Levine BG. Simulating Electron Dynamics of Complex Molecules with Time-Dependent Complete Active Space Configuration Interaction. J Chem Theory Comput 2018; 14:4129-4138. [PMID: 29986143 DOI: 10.1021/acs.jctc.8b00381] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Time-dependent electronic structure methods are growing in popularity as tools for modeling ultrafast and/or nonlinear processes, for computing spectra, and as the electronic structure component of mean-field molecular dynamics simulations. Time-dependent configuration interaction (TD-CI) offers several advantages over the widely used real-time time-dependent density functional theory: namely, that it correctly models Rabi oscillations; it offers a spin-pure description of open-shell systems; and a hierarchy of TD-CI methods can be defined that systematically approach the exact solution of the time-dependent Schrodinger equation (TDSE). In this work, we present a novel TD-CI approach that extends TD-CI to large complete active-space configuration expansions. Such extension is enabled by use of a direct configuration interaction approach that eliminates the need to explicitly build, store, or diagonalize the Hamiltonian matrix. Graphics processing unit (GPU) acceleration enables fast solution of the TDSE even for large active spaces-up to 12 electrons in 12 orbitals (853776 determinants) in this work. A symplectic split operator propagator yields long-time norm conservation. We demonstrate the applicability of our approach by computing the response of a large molecule with a strongly correlated ground state, decacene (C42H24), to various pulses (δ-function, transform limited, chirped). Our simulations predict that chirped pulses can be used to induce dipole-forbidden transitions. Simulations of decacene using the 6-31G(d) basis set and a 12 electrons/12 orbitals active space took 20.1 h to propagate for 100 fs with a 1 attosecond time step on a single NVIDIA K40 GPU. Convergence with respect to time step is found to depend on the property being computed and the chosen active space.
Collapse
Affiliation(s)
- Wei-Tao Peng
- Department of Chemistry , Michigan State University , East Lansing , Michigan 48824 , United States
| | - B Scott Fales
- Department of Chemistry and the PULSE Institute , Stanford University , Stanford , California 94305 , United States.,SLAC National Accelerator Laboratory , Menlo Park , California 94025 , United States
| | - Benjamin G Levine
- Department of Chemistry , Michigan State University , East Lansing , Michigan 48824 , United States
| |
Collapse
|
28
|
Bonafé FP, Hernández FJ, Aradi B, Frauenheim T, Sánchez CG. Fully Atomistic Real-Time Simulations of Transient Absorption Spectroscopy. J Phys Chem Lett 2018; 9:4355-4359. [PMID: 30024765 DOI: 10.1021/acs.jpclett.8b01659] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We have implemented an electron-nuclear real-time propagation scheme for the calculation of transient absorption spectra. When this technique is applied to the study of ultrafast dynamics of Soret-excited zinc(II) tetraphenylporphyrin in the subpicosecond time scale, quantum beats in the transient absorption caused by impulsively excited molecular vibrations are observed. The launching mechanism of such vibrations can be regarded as a displacive excitation of the zinc-pyrrole and pyrrole C-C bonds.
Collapse
Affiliation(s)
- Franco P Bonafé
- Departamento de Química Teórica y Computacional, Facultad de Ciencias Químicas , Universidad Nacional de Córdoba , Córdoba , Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba , INFIQC (CONICET - Universidad Nacional de Córdoba) , Córdoba , Argentina
| | - Federico J Hernández
- Departamento de Química Teórica y Computacional, Facultad de Ciencias Químicas , Universidad Nacional de Córdoba , Córdoba , Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba , INFIQC (CONICET - Universidad Nacional de Córdoba) , Córdoba , Argentina
| | - Bálint Aradi
- Bremen Center for Computational Materials Science , Universität Bremen , Bremen , Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science , Universität Bremen , Bremen , Germany
| | - Cristián G Sánchez
- Departamento de Química Teórica y Computacional, Facultad de Ciencias Químicas , Universidad Nacional de Córdoba , Córdoba , Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba , INFIQC (CONICET - Universidad Nacional de Córdoba) , Córdoba , Argentina
| |
Collapse
|
29
|
Chalabala J, Uhlig F, Slavíček P. Assessment of Real-Time Time-Dependent Density Functional Theory (RT-TDDFT) in Radiation Chemistry: Ionized Water Dimer. J Phys Chem A 2018. [PMID: 29513531 DOI: 10.1021/acs.jpca.8b01259] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Ionization in the condensed phase and molecular clusters leads to a complicated chain of processes with coupled electron-nuclear dynamics. It is difficult to describe such dynamics with conventional nonadiabatic molecular dynamics schemes since the number of states swiftly increases as the molecular system grows. It is therefore attractive to use a direct electron and nuclear propagation such as the real-time time-dependent density functional theory (RT-TDDFT). Here we report a RT-TDDFT benchmark study on simulations of singly and doubly ionized states of a water monomer and dimer as a prototype for more complex processes in a condensed phase. We employed the RT-TDDFT based Ehrenfest molecular dynamics with a generalized gradient approximate (GGA) functional and compared it with wave-function-based surface hopping (SH) simulations. We found that the initial dynamics of a singly HOMO ionized water dimer is similar for both the RT-TDDFT/GGA and the SH simulations but leads to completely different reaction channels on a longer time scale. This failure is attributed to the self-interaction error in the GGA functionals and it can be avoided by using hybrid functionals with large fraction of exact exchange (represented here by the BHandHLYP functional). The simulations of doubly ionized states are reasonably described already at the GGA level. This suggests that the RT-TDDFT/GGA method could describe processes following the autoionization processes such as Auger emission, while its applicability to more complex processes such as intermolecular Coulombic decay remains limited.
Collapse
Affiliation(s)
- Jan Chalabala
- Department of Physical Chemistry , University of Chemistry and Technology , Technická 5 , 16628 Prague , Czech Republic
| | - Frank Uhlig
- Department of Physical Chemistry , University of Chemistry and Technology , Technická 5 , 16628 Prague , Czech Republic.,Institute for Computational Physics , University of Stuttgart , Allmandring 3 , 70569 Stuttgart , Germany
| | - Petr Slavíček
- Department of Physical Chemistry , University of Chemistry and Technology , Technická 5 , 16628 Prague , Czech Republic.,Jaroslav Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic , Dolejškova 3 , 18200 Prague , Czech Republic
| |
Collapse
|
30
|
Kretchmer JS, Chan GKL. A real-time extension of density matrix embedding theory for non-equilibrium electron dynamics. J Chem Phys 2018; 148:054108. [DOI: 10.1063/1.5012766] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Joshua S. Kretchmer
- Divsion of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Garnet Kin-Lic Chan
- Divsion of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
31
|
Bowman DN, Asher JC, Fischer SA, Cramer CJ, Govind N. Excited-state absorption in tetrapyridyl porphyrins: comparing real-time and quadratic-response time-dependent density functional theory. Phys Chem Chem Phys 2018; 19:27452-27462. [PMID: 28975162 DOI: 10.1039/c7cp04567k] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Three meso-substituted tetrapyridyl porphyrins (free base, Ni(ii), and Cu(ii)) were investigated for their optical limiting (OL) capabilities using real-time (RT-), linear-response (LR-), and quadratic-response (QR-) time-dependent density functional theory (TDDFT) methods. These species are experimentally known to display a prominent reverse saturable absorption feature between the Q and B bands of the ground-state absorption (GSA), which has been attributed to increased excited-state absorption (ESA) relative to GSA. A recently developed RT-TDDFT based method for calculating ESA from a LR-TDDFT density was utilized with eight exchange-correlation functionals (BLYP, PBE, B3LYP, CAM-B3LYP, PBE0, M06, BHLYP, and BHandH) and contrasted with calculations of ESA using QR-TDDFT with five exchange-correlation functionals (BLYP, B3LYP, CAM-B3LYP, BHLYP, and BHandH). This allowed for comparison between functionals with varying amounts of exact exchange as well as between the ability of RT-TDDFT and QR-TDDFT to reproduce OL behavior in porphyrin systems. The absorption peak positions and intensities for GSA and ESA are significantly impacted by the choice of DFT functional, with the most critical factor identified as the amount of exact exchange in the functional form. Calculating ESA with QR-TDDFT is found to be significantly more sensitive to the amount of exact exchange than GSA and ESA with RT-TDDFT, as well as GSA with LR-TDDFT. An analogous behavior is also demonstrated for the polycyclic aromatic hydrocarbon coronene. This is problematic when using the same approximate functional for calculation of both GSA and ESA, as the LR- and QR-TDDFT excitation energies will not have similar errors. Overall, the RT-TDDFT method with hybrid functionals reproduces the OL features for the porphyrin systems studied here and is a viable computational approach for efficient screening of molecular complexes for OL properties.
Collapse
Affiliation(s)
- David N Bowman
- Department of Chemistry, Supercomputing Institute and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota 55455, USA.
| | | | | | | | | |
Collapse
|
32
|
Lestrange PJ, Hoffmann MR, Li X. Time-Dependent Configuration Interaction Using the Graphical Unitary Group Approach: Nonlinear Electric Properties. NOVEL ELECTRONIC STRUCTURE THEORY: GENERAL INNOVATIONS AND STRONGLY CORRELATED SYSTEMS 2018. [DOI: 10.1016/bs.aiq.2017.06.003] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
33
|
Goings JJ, Lestrange PJ, Li X. Real‐time time‐dependent electronic structure theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1341] [Citation(s) in RCA: 84] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
| | | | - Xiaosong Li
- Department of ChemistryUniversity of Washington Seattle WA USA
| |
Collapse
|
34
|
Rossi TP, Kuisma M, Puska MJ, Nieminen RM, Erhart P. Kohn–Sham Decomposition in Real-Time Time-Dependent Density-Functional Theory: An Efficient Tool for Analyzing Plasmonic Excitations. J Chem Theory Comput 2017; 13:4779-4790. [DOI: 10.1021/acs.jctc.7b00589] [Citation(s) in RCA: 72] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Tuomas P. Rossi
- COMP
Centre of Excellence, Department of Applied Physics, Aalto University, P.O.
Box 11100, FI-00076 Aalto, Finland
| | - Mikael Kuisma
- Department
of Physics, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
- Department
of Chemistry, Nanoscience Center, University of Jyväskylä, FI-40014 Jyväskylä, Finland
| | - Martti J. Puska
- COMP
Centre of Excellence, Department of Applied Physics, Aalto University, P.O.
Box 11100, FI-00076 Aalto, Finland
| | - Risto M. Nieminen
- COMP
Centre of Excellence, Department of Applied Physics, Aalto University, P.O.
Box 11100, FI-00076 Aalto, Finland
| | - Paul Erhart
- Department
of Physics, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
| |
Collapse
|
35
|
Bruner A, Hernandez S, Mauger F, Abanador PM, LaMaster DJ, Gaarde MB, Schafer KJ, Lopata K. Attosecond Charge Migration with TDDFT: Accurate Dynamics from a Well-Defined Initial State. J Phys Chem Lett 2017; 8:3991-3996. [PMID: 28792225 DOI: 10.1021/acs.jpclett.7b01652] [Citation(s) in RCA: 46] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
We investigate the ability of time-dependent density functional theory (TDDFT) to capture attosecond valence electron dynamics resulting from sudden X-ray ionization of a core electron. In this special case the initial state can be constructed unambiguously, allowing for a simple test of the accuracy of the dynamics. The response following nitrogen K-edge ionization in nitrosobenzene shows excellent agreement with fourth-order algebraic diagrammatic construction (ADC(4)) results, suggesting that a properly chosen initial state allows TDDFT to adequately capture attosecond charge migration. Visualizing hole motion using an electron localization picture (ELF), we provide an intuitive chemical interpretation of the charge migration as a superposition of Lewis dot resonance structures.
Collapse
Affiliation(s)
- Adam Bruner
- Department of Chemistry, Louisiana State University , Baton Rouge, Louisiana 70803, United States
| | - Samuel Hernandez
- Department of Chemistry, Louisiana State University , Baton Rouge, Louisiana 70803, United States
| | - François Mauger
- Department of Physics and Astronomy, Louisiana State University , Baton Rouge, Louisiana 70803, United States
| | - Paul M Abanador
- Department of Physics and Astronomy, Louisiana State University , Baton Rouge, Louisiana 70803, United States
| | - Daniel J LaMaster
- Department of Chemistry, Louisiana State University , Baton Rouge, Louisiana 70803, United States
| | - Mette B Gaarde
- Department of Physics and Astronomy, Louisiana State University , Baton Rouge, Louisiana 70803, United States
| | - Kenneth J Schafer
- Department of Physics and Astronomy, Louisiana State University , Baton Rouge, Louisiana 70803, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University , Baton Rouge, Louisiana 70803, United States
- Center for Computation & Technology, Louisiana State University , Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
36
|
Mosquera MA, Jackson NE, Fauvell TJ, Kelley MS, Chen LX, Schatz GC, Ratner MA. Exciton Absorption Spectra by Linear Response Methods: Application to Conjugated Polymers. J Am Chem Soc 2017; 139:3728-3735. [DOI: 10.1021/jacs.6b12405] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Martín A. Mosquera
- Department
of Chemistry and the Materials Research Center, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Nicholas E. Jackson
- Department
of Chemistry and the Materials Research Center, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
- Chemical
Sciences and Engineering Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439, United States
| | - Thomas J. Fauvell
- Department
of Chemistry and the Materials Research Center, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew S. Kelley
- Department
of Chemistry and the Materials Research Center, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Lin X. Chen
- Department
of Chemistry and the Materials Research Center, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
- Chemical
Sciences and Engineering Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439, United States
| | - George C. Schatz
- Department
of Chemistry and the Materials Research Center, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Mark A. Ratner
- Department
of Chemistry and the Materials Research Center, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| |
Collapse
|
37
|
Sissay A, Abanador P, Mauger F, Gaarde M, Schafer KJ, Lopata K. Angle-dependent strong-field molecular ionization rates with tuned range-separated time-dependent density functional theory. J Chem Phys 2016; 145:094105. [DOI: 10.1063/1.4961731] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Adonay Sissay
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Paul Abanador
- Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - François Mauger
- Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Mette 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
|
38
|
Luo K, Fuks JI, Maitra NT. Studies of spuriously shifting resonances in time-dependent density functional theory. J Chem Phys 2016; 145:044101. [DOI: 10.1063/1.4955447] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Kai Luo
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
- Department of Chemistry and Chemical Biology, Cornell University, 259 East Ave., Ithaca, New York 14850, USA
| | - Johanna I. Fuks
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| | - Neepa T. Maitra
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| |
Collapse
|
39
|
Zhang GP, Bai YH, George TF. Ultrafast reduction of exchange splitting in ferromagnetic nickel. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:236004. [PMID: 27160931 DOI: 10.1088/0953-8984/28/23/236004] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
A decade ago Rhie et al (2003 Phys. Rev. Lett. 90 247201) reported that when ferromagnetic nickel is subject to an intense ultrashort laser pulse, its exchange splitting is reduced quickly. But to simulate such reduction remains a big challenge. The popular rigid band approximation (RBA), where both the band structure and the exchange splitting are held fixed before and after laser excitation, is unsuitable for this purpose, while the time-dependent density functional theory could be time-consuming. To overcome these difficulties, we propose a time-dependent Liouville and density functional theory (TDLDFT) that integrates the time-dependent Liouville equation into the density functional theory. As a result, the excited charge density is reiterated back into the Kohn-Sham equation, and the band structure is allowed to change dynamically. Even with the ground-state density functional, a larger demagnetization than RBA is found; after we expand Ortenzi's spin scaling method into an excited-state (laser) density functional, we find that the exchange splitting is indeed strongly reduced, as seen in the experiment. Both the majority and minority bands are shifted toward the Fermi level, but the majority shifts a lot more. The ultrafast reduction in exchange splitting occurs concomitantly with demagnetization. While our current theory is still unable to yield the same percentage loss in the spin moment as observed in the experiment, it predicts a correct trend that agrees with the experiments. With a better functional, we believe that our results can be further improved.
Collapse
Affiliation(s)
- G P Zhang
- Department of Physics, Indiana State University, Terre Haute, IN 47809, USA
| | | | | |
Collapse
|