1
|
Neufeld O, Tancogne-Dejean N, Rubio A. Benchmarking Functionals for Strong-Field Light-Matter Interactions in Adiabatic Time-Dependent Density Functional Theory. J Phys Chem Lett 2024; 15:7254-7264. [PMID: 38976844 PMCID: PMC11261632 DOI: 10.1021/acs.jpclett.4c01383] [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/11/2024] [Revised: 06/29/2024] [Accepted: 07/03/2024] [Indexed: 07/10/2024]
Abstract
In recent years, time-dependent density functional theory (TDDFT) has been extensively employed for highly nonlinear optics in molecules and solids, including high harmonic generation (HHG), photoemission, and more. TDDFT exhibits a relatively low numerical cost while still describing both light-matter and electron-electron interactions ab initio, making it highly appealing. However, the majority of implementations of the theory utilize the simplest possible approximations for the exchange-correlation (XC) functional-either the local density or generalized gradient approximations, which are traditionally considered to have rather poor chemical accuracy. We present the first systematic study of the XC functional effect on molecular HHG, testing various levels of theory. Our numerical results suggest justification for using simpler approximations for the XC functional, showing that hybrid and meta functionals (as well as Hartree-Fock) can, at times, lead to poor and unphysical results. The specific source of the failure in more elaborate functionals should be topic of future work, but we hypothesize that its origin might be connected to the adiabatic approximation of TDDFT.
Collapse
Affiliation(s)
- Ofer Neufeld
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
| | - Nicolas Tancogne-Dejean
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
| | - Angel Rubio
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
- Center
for Computational Quantum Physics (CCQ), The Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|
2
|
de las Heras A, Bonafé FP, Hernández-García C, Rubio A, Neufeld O. Tunable Tesla-Scale Magnetic Attosecond Pulses through Ring-Current Gating. J Phys Chem Lett 2023; 14:11160-11167. [PMID: 38054653 PMCID: PMC10726360 DOI: 10.1021/acs.jpclett.3c02899] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2023] [Revised: 11/28/2023] [Accepted: 11/30/2023] [Indexed: 12/07/2023]
Abstract
Coherent control over electron dynamics in atoms and molecules using high-intensity circularly polarized laser pulses gives rise to current loops, resulting in the emission of magnetic fields. We propose, and demonstrate with ab initio calculations, "current-gating" schemes to generate direct or alternating-current magnetic pulses in the infrared spectral region, with highly tunable waveform and frequency, and showing femtosecond-to-attosecond pulse duration. In optimal conditions, the magnetic pulse can be highly isolated from the driving laser and exhibits a high flux density (∼1 T at a few hundred nanometers from the source, with a pulse duration of 787 attoseconds) for application in forefront experiments of ultrafast spectroscopy. Our work paves the way toward the generation of attosecond magnetic fields to probe ultrafast magnetization, chiral responses, and spin dynamics.
Collapse
Affiliation(s)
- Alba de las Heras
- Grupo
de Investigación en Aplicaciones del Láser y Fotónica,
Departamento de Física Aplicada, Universidad de Salamanca, Salamanca 37008, Spain
| | - Franco P. Bonafé
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
| | - Carlos Hernández-García
- Grupo
de Investigación en Aplicaciones del Láser y Fotónica,
Departamento de Física Aplicada, Universidad de Salamanca, Salamanca 37008, Spain
| | - Angel Rubio
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
- Center
for Computational Quantum Physics, The Flatiron
Institute, New York 10010, United States
- Nano-Bio
Spectroscopy Group, Departamento de Física de Materiales, Universidad del País Vasco, San Sebastían 20018, Spain
| | - Ofer Neufeld
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
| |
Collapse
|
3
|
Hellgren M, Baguet L. Strengths and limitations of the adiabatic exact-exchange kernel for total energy calculations. J Chem Phys 2023; 158:2889488. [PMID: 37158324 DOI: 10.1063/5.0146423] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Accepted: 04/24/2023] [Indexed: 05/10/2023] Open
Abstract
We investigate the adiabatic approximation to the exact-exchange kernel for calculating correlation energies within the adiabatic-connection fluctuation-dissipation framework of time-dependent density functional theory. A numerical study is performed on a set of systems having bonds of different character (H2 and N2 molecules, H-chain, H2-dimer, solid-Ar, and the H2O-dimer). We find that the adiabatic kernel can be sufficient in strongly bound covalent systems, yielding similar bond lengths and binding energies. However, for non-covalent systems, the adiabatic kernel introduces significant errors around equilibrium geometry, systematically overestimating the interaction energy. The origin of this behavior is investigated by studying a model dimer composed of one-dimensional, closed-shell atoms, interacting via soft-Coulomb potentials. The kernel is shown to exhibit a strong frequency dependence at small to intermediate atomic separation that affects both the low-energy spectrum and the exchange-correlation hole obtained from the corresponding diagonal of the two-particle density matrix.
Collapse
Affiliation(s)
- Maria Hellgren
- Sorbonne Université, MNHN, UMR CNRS 7590, IMPMC, 4 place Jussieu, 75005 Paris, France
| | - Lucas Baguet
- CEA, DAM, DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière en Conditions Extrêmes, 91680 Bruyères-le-Châtel, France
| |
Collapse
|
4
|
Electronic excited states in deep variational Monte Carlo. Nat Commun 2023; 14:274. [PMID: 36650151 PMCID: PMC9845370 DOI: 10.1038/s41467-022-35534-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Accepted: 12/08/2022] [Indexed: 01/19/2023] Open
Abstract
Obtaining accurate ground and low-lying excited states of electronic systems is crucial in a multitude of important applications. One ab initio method for solving the Schrödinger equation that scales favorably for large systems is variational quantum Monte Carlo (QMC). The recently introduced deep QMC approach uses ansatzes represented by deep neural networks and generates nearly exact ground-state solutions for molecules containing up to a few dozen electrons, with the potential to scale to much larger systems where other highly accurate methods are not feasible. In this paper, we extend one such ansatz (PauliNet) to compute electronic excited states. We demonstrate our method on various small atoms and molecules and consistently achieve high accuracy for low-lying states. To highlight the method's potential, we compute the first excited state of the much larger benzene molecule, as well as the conical intersection of ethylene, with PauliNet matching results of more expensive high-level methods.
Collapse
|
5
|
Albareda G, Lively K, Sato SA, Kelly A, Rubio A. Conditional Wave Function Theory: A Unified Treatment of Molecular Structure and Nonadiabatic Dynamics. J Chem Theory Comput 2021; 17:7321-7340. [PMID: 34752108 PMCID: PMC8675140 DOI: 10.1021/acs.jctc.1c00772] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Indexed: 11/28/2022]
Abstract
We demonstrate that a conditional wave function theory enables a unified and efficient treatment of the equilibrium structure and nonadiabatic dynamics of correlated electron-ion systems. The conditional decomposition of the many-body wave function formally recasts the full interacting wave function of a closed system as a set of lower-dimensional (conditional) coupled "slices". We formulate a variational wave function ansatz based on a set of conditional wave function slices and demonstrate its accuracy by determining the structural and time-dependent response properties of the hydrogen molecule. We then extend this approach to include time-dependent conditional wave functions and address paradigmatic nonequilibrium processes including strong-field molecular ionization, laser-driven proton transfer, and nuclear quantum effects induced by a conical intersection. This work paves the road for the application of conditional wave function theory in equilibrium and out-of-equilibrium ab initio molecular simulations of finite and extended systems.
Collapse
Affiliation(s)
- Guillermo Albareda
- Nano-Bio
Spectroscopy Group and European Theoretical Spectroscopy Facility
(ETSF), Universidad del País Vasco
(UPV/EHU), Av. Tolosa
72, 20018 San Sebastian, Spain
- Institute
of Theoretical and Computational Chemistry, University of Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Kevin Lively
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
- The
Hamburg Centre for Ultrafast Imaging, University
of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Shunsuke A. Sato
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
- Center
for Computational Sciences, University of
Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| | - Aaron Kelly
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
- The
Hamburg Centre for Ultrafast Imaging, University
of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
- Department
of Chemistry, Dalhousie University, Halifax, Nova Scotia B3H 4R2, Canada
| | - Angel Rubio
- Nano-Bio
Spectroscopy Group and European Theoretical Spectroscopy Facility
(ETSF), Universidad del País Vasco
(UPV/EHU), Av. Tolosa
72, 20018 San Sebastian, Spain
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
- The
Hamburg Centre for Ultrafast Imaging, University
of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
- Center
for Computational Quantum Physics (CCQ), Flatiron Institute, 162 Fifth Avenue, New York, New York 10010, United
States
| |
Collapse
|
6
|
Direct treatment of interaction between laser-field and electrons for simulating laser processing of metals. Sci Rep 2021; 11:14626. [PMID: 34272433 PMCID: PMC8285478 DOI: 10.1038/s41598-021-94036-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Accepted: 07/06/2021] [Indexed: 11/08/2022] Open
Abstract
Laser ablation is often simulated by the two-temperature model in which electrons are assumed to be thermalized by laser irradiation, while an explicit representation of interaction between laser-field and electrons is challenging but beneficial as being free from any adjustable parameters. Here, an ab initio method based on the time-dependent density functional theory (TDDFT) in which electron-ion dynamics under a laser field are numerically simulated is examined as a tool for simulating femtosecond laser processing of metals. Laser-induced volume expansion in surface normal directions of Cu(111) and Ni(111) surfaces are simulated by using repeating slab models. The amount of simulated volume expansion is compared between Cu(111) and Ni(111) slabs for the same laser pulse conditions, and the Ni slab is found to expand more than the Cu slab despite the smaller thermal expansion coefficient of Ni compared with Cu. The analyzed electronic excitation and lattice motion were compared to those in the two-temperature model. The threshold fluence to release surface Cu atom deduced from current TDDFT approach is found to be comparable to those of Cu ablation reported experimentally.
Collapse
|
7
|
Yao X, Lee Y, Ceresoli D, Cho K. First-Principles Study on Electron-Induced Excitations of Atomic Layer Deposition Precursors: Inelastic Electron Wave Packet Scattering with Cobalt Tricarbonyl Nitrosyl Co(CO) 3NO Using Time-Dependent Density Functional Theory. J Phys Chem A 2021; 125:4524-4533. [PMID: 34019398 DOI: 10.1021/acs.jpca.0c11309] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A quantitative study on inelastic electron scattering with a molecule is of significant importance for understanding the essential mechanisms of electron-induced gas-phase and surface chemical reactions in their excited electronic states. A key issue to be addressed is the quantitatively detailed inelastic electron collision processes with a realistic molecular target, associated with electron excitation that leads to potential ionization and dissociation reactions of the molecule. Using the real-time time-dependent density functional theory (TDDFT) modeling, we present quantitative findings on the energy transfers and internal excitations for the low energy (up to 270 eV) electron wave packet impact with the molecular target cobalt tricarbonyl nitrosyl (CTN, Co(CO)3NO) that is used as a precursor in electron-enhanced atomic layer deposition (EE-ALD) growth of Co films. Our modeling shows the quantitative dependence of the wave packet sizes, target molecule orientations, and impact parameters on the energy transfer in this inelastic electron scattering process. It is found that the wave packet sizes have little effect on the overall profile of the internal multiple excited states, whereas different target orientations can cause significantly different internal excited states. To evaluate the quantitative prediction capability, the inelastic scattering cross-section of a hydrogen atom is calculated and compared with the experimental data, leading to a constant scaling factor over the whole energy range. The present study demonstrates the remarkable potential of TDDFT for simulating the inelastic electron scattering process, which provides critical information for future exploration of electronic excitations in a wide range of electron-induced chemical reactions in current technological applications.
Collapse
Affiliation(s)
- Xiaolong Yao
- Department of Materials Science and Engineering, University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Yeonghun Lee
- Department of Materials Science and Engineering, University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Davide Ceresoli
- CNR-SCITEC, Istituto di Scienze e Tecnologie Chimiche "G. Natta", via Golgi 19, 20133 Milan, Italy
| | - Kyeongjae Cho
- Department of Materials Science and Engineering, University of Texas at Dallas, Richardson, Texas 75080, United States
| |
Collapse
|
8
|
Bhan L, Covington C, Rivas J, Varga K. Simulation of photo-electron spectrum and electron scattering by dual time propagation. J Chem Phys 2021; 154:114110. [PMID: 33752384 DOI: 10.1063/5.0045591] [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/15/2022] Open
Abstract
A dual time propagation approach is introduced to describe electron scattering and ionization. The space is divided into two regions, a central region with a full time-dependent Hamiltonian and an outer region where the kinetic operator and the laser field dominate. The two regions are connected by a source term. Time-dependent density functional theory calculations of wave packet scattering on molecules and photoelectron spectrum due to circularly polarized laser are presented to illustrate the efficiency and applicability of the approach.
Collapse
Affiliation(s)
- Luke Bhan
- Department of Physics and Astronomy, Vanderbilt University, Nashville, Tennessee 37235, USA
| | - Cody Covington
- Department of Chemistry, Austin Peay State University, Clarksville, Tennessee 37044, USA
| | - Jason Rivas
- Department of Chemistry, Austin Peay State University, Clarksville, Tennessee 37044, USA
| | - Kálmán Varga
- Department of Physics and Astronomy, Vanderbilt University, Nashville, Tennessee 37235, USA
| |
Collapse
|
9
|
Lacombe L, Maitra NT. Developing new and understanding old approximations in TDDFT. Faraday Discuss 2020; 224:382-401. [PMID: 32926040 DOI: 10.1039/d0fd00049c] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
When a system has evolved far from a ground-state, the adiabatic approximations commonly used in time-dependent density functional theory calculations completely fail in some applications, while giving qualitatively good predictions in others, and sometimes even quantitative predictions. It is not clearly understood why this is so, and developing practical approximations going beyond the adiabatic approximation remains a challenge. This paper explores three different lines of investigation. First, an expression for the exact time-dependent exchange-correlation potential suggests that the accuracy of an adiabatic approximation is intimately related to the deviation between the natural orbital occupation numbers of the physical system and those of the Kohn-Sham system, and we explore this on some exactly-solvable model systems. The exact expression further suggests a path to go beyond the adiabatic approximations, and in the second part we discuss a newly proposed class of memory-dependent approximations developed in this way. Finally, we derive a new expression for the exact exchange-correlation potential from a coupling-constant path integration.
Collapse
Affiliation(s)
- Lionel Lacombe
- Department of Physics, Rutgers University, Newark, New Jersey 07102, USA.
| | | |
Collapse
|
10
|
Lee Y, Yao X, Fischetti MV, Cho K. Real-time ab initio simulation of inelastic electron scattering using the exact, density functional, and alternative approaches. Phys Chem Chem Phys 2020; 22:8616-8624. [DOI: 10.1039/c9cp06376e] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Inelastic electron scattering phenomena in chemical/physical/materials interests: electron radiation damage in materials; DNA damaged by electron scattering; electron therapy; electron microscope; electron-beam-induced deposition for nanofabrication.
Collapse
Affiliation(s)
- Yeonghun Lee
- Department of Materials Science and Engineering
- University of Texas at Dallas
- Richardson
- USA
| | - Xiaolong Yao
- Department of Materials Science and Engineering
- University of Texas at Dallas
- Richardson
- USA
| | - Massimo V. Fischetti
- Department of Materials Science and Engineering
- University of Texas at Dallas
- Richardson
- USA
| | - Kyeongjae Cho
- Department of Materials Science and Engineering
- University of Texas at Dallas
- Richardson
- USA
| |
Collapse
|
11
|
Lacombe L, Maitra NT. Density-Matrix Coupled Time-Dependent Exchange-Correlation Functional Approximations. J Chem Theory Comput 2019; 15:1672-1678. [DOI: 10.1021/acs.jctc.8b01159] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Lionel Lacombe
- 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, United States
| | - 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, United States
| |
Collapse
|
12
|
Suzuki Y, Hagiwara S, Watanabe K. Time-Dependent Multicomponent Density Functional Theory for Coupled Electron-Positron Dynamics. PHYSICAL REVIEW LETTERS 2018; 121:133001. [PMID: 30312034 DOI: 10.1103/physrevlett.121.133001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/12/2018] [Indexed: 06/08/2023]
Abstract
Electron-positron interactions have been utilized in various fields of science. Here we develop time-dependent multicomponent density functional theory to study the coupled electron-positron dynamics from first principles. We prove that there are coupled time-dependent single-particle equations that can provide the electron and positron density dynamics, and derive the formally exact expression for their effective potentials. Introducing the adiabatic local density approximation to time-dependent electron-positron correlation, we apply the theory to the dynamics of a positronic lithium hydride molecule under a laser field. We demonstrate the significance of the coupling between electronic and positronic motion by revealing the complex positron detachment mechanism and the suppression of electronic resonant excitation by the screening effect of the positron.
Collapse
Affiliation(s)
- Yasumitsu Suzuki
- Department of Physics, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| | - Satoshi Hagiwara
- Department of Physics, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| | - Kazuyuki Watanabe
- Department of Physics, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| |
Collapse
|
13
|
Fuks JI, Lacombe L, Nielsen SEB, Maitra NT. Exploring non-adiabatic approximations to the exchange–correlation functional of TDDFT. Phys Chem Chem Phys 2018; 20:26145-26160. [DOI: 10.1039/c8cp03957g] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Decomposition of the exact time-dependent exchange–correlation potential offers a new starting point to build approximations with memory.
Collapse
Affiliation(s)
- Johanna I. Fuks
- Departamento de Física and IFIBA
- FCEN
- Universidad de Buenos Aires
- Ciudad Universitaria
- C1428EHA Ciudad de Buenos Aires
| | - Lionel Lacombe
- Department of Physics and Astronomy
- Hunter College and the Graduate Center of the City University of New York
- New York
- USA
| | - Søren E. B. Nielsen
- Max Planck Institute for the Structure and Dynamics of Matter and Center for Free-Electron Laser Science
- 22761 Hamburg
- Germany
| | - Neepa T. Maitra
- Department of Physics and Astronomy
- Hunter College and the Graduate Center of the City University of New York
- New York
- USA
| |
Collapse
|