1
|
Worth GA, Robb MA. Controlling Electronic Coherences and the Curvature Induced by the Derivative Coupling at a Conical Intersection: A Quantum Ehrenfest (QuEh) Protocol for Reaction Path Following Application to "Channel 3" Benzene Photochemistry. J Phys Chem A 2024. [PMID: 38917388 DOI: 10.1021/acs.jpca.4c02449] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/27/2024]
Abstract
We report a protocol for the implementation of "reaction path following" from a transition state through a conical intersection, including both the path curvature induced by the derivative coupling and the corresponding induced electronic coherences. This protocol focuses on the "central" Gaussian wavepacket (initially unexcited) in the quantum Ehrenfest (QuEh) method. Like the reaction path following, the normal mode corresponding to the imaginary frequency at the transition state is given an initial momentum. The protocol is applied to the "channel 3" radiationless decay of benzene. We also demonstrate that one can enhance the effect of the derivative coupling and the electronic coherence with an IR pulse.
Collapse
Affiliation(s)
- Graham A Worth
- Department of Chemistry, University College London, 20, Gordon Street, WC1H 0AJ London, U.K
| | - Michael A Robb
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, W12 0BZ London, U.K
| |
Collapse
|
2
|
Solov'yov AV, Verkhovtsev AV, Mason NJ, Amos RA, Bald I, Baldacchino G, Dromey B, Falk M, Fedor J, Gerhards L, Hausmann M, Hildenbrand G, Hrabovský M, Kadlec S, Kočišek J, Lépine F, Ming S, Nisbet A, Ricketts K, Sala L, Schlathölter T, Wheatley AEH, Solov'yov IA. Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment. Chem Rev 2024. [PMID: 38842266 DOI: 10.1021/acs.chemrev.3c00902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2024]
Abstract
This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.
Collapse
Affiliation(s)
- Andrey V Solov'yov
- MBN Research Center, Altenhöferallee 3, 60438 Frankfurt am Main, Germany
| | | | - Nigel J Mason
- School of Physics and Astronomy, University of Kent, Canterbury CT2 7NH, United Kingdom
| | - Richard A Amos
- Department of Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, U.K
| | - Ilko Bald
- Institute of Chemistry, University of Potsdam, Karl-Liebknecht-Str. 24-25, 14476 Potsdam, Germany
| | - Gérard Baldacchino
- Université Paris-Saclay, CEA, LIDYL, 91191 Gif-sur-Yvette, France
- CY Cergy Paris Université, CEA, LIDYL, 91191 Gif-sur-Yvette, France
| | - Brendan Dromey
- Centre for Light Matter Interactions, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, United Kingdom
| | - Martin Falk
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 61200 Brno, Czech Republic
- Kirchhoff-Institute for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
| | - Juraj Fedor
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Luca Gerhards
- Institute of Physics, Carl von Ossietzky University, Carl-von-Ossietzky-Str. 9-11, 26129 Oldenburg, Germany
| | - Michael Hausmann
- Kirchhoff-Institute for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
| | - Georg Hildenbrand
- Kirchhoff-Institute for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
- Faculty of Engineering, University of Applied Sciences Aschaffenburg, Würzburger Str. 45, 63743 Aschaffenburg, Germany
| | | | - Stanislav Kadlec
- Eaton European Innovation Center, Bořivojova 2380, 25263 Roztoky, Czech Republic
| | - Jaroslav Kočišek
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Franck Lépine
- Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeurbanne, France
| | - Siyi Ming
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Andrew Nisbet
- Department of Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, U.K
| | - Kate Ricketts
- Department of Targeted Intervention, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Leo Sala
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Thomas Schlathölter
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
- University College Groningen, University of Groningen, Hoendiepskade 23/24, 9718 BG Groningen, The Netherlands
| | - Andrew E H Wheatley
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Ilia A Solov'yov
- Institute of Physics, Carl von Ossietzky University, Carl-von-Ossietzky-Str. 9-11, 26129 Oldenburg, Germany
| |
Collapse
|
3
|
Guo Z, Zhang Z, Deng Y, Wang J, Ye D, Liu J, Liu Y. Probing H_{2} Double Ionization with Bicircular Laser Fields. PHYSICAL REVIEW LETTERS 2024; 132:143201. [PMID: 38640361 DOI: 10.1103/physrevlett.132.143201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/04/2023] [Revised: 11/30/2023] [Accepted: 03/14/2024] [Indexed: 04/21/2024]
Abstract
We present a kinematically complete study on strong-field double ionization of H_{2} molecules in two-color bicircular laser fields. The releasing times of electrons and protons are recorded with the double-hand attoclock. We observe the relative emission angles of two electrons oscillate with the kinetic energy release of protons, indicating the internal concerted four-body fragmentation. Using a three-dimensional molecular semiclassical ensemble model, we have disentangled the attosecond correlated electron emission in H_{2} double ionization. This work reveals the strong electron-nuclear coupling in the molecular bond breaking and may open up a new approach to experimentally accessing the intramolecular electron and bond dynamics with bicircular fields.
Collapse
Affiliation(s)
- Zhenning Guo
- State Key Laboratory for Mesoscopic Physics and Frontiers Science Center for Nano-optoelectronics, School of Physics, Peking University, Beijing 100871, China
| | - Zhihe Zhang
- Center for Applied Physics and Technology, HEDPS, Peking University, Beijing 100871, China
| | - Yongkai Deng
- State Key Laboratory for Mesoscopic Physics and Frontiers Science Center for Nano-optoelectronics, School of Physics, Peking University, Beijing 100871, China
| | - Jiguo Wang
- State Key Laboratory for Mesoscopic Physics and Frontiers Science Center for Nano-optoelectronics, School of Physics, Peking University, Beijing 100871, China
| | - Difa Ye
- National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
| | - Jie Liu
- Center for Applied Physics and Technology, HEDPS, Peking University, Beijing 100871, China
- Graduate School, China Academy of Engineering Physics, Beijing 100193, China
| | - Yunquan Liu
- State Key Laboratory for Mesoscopic Physics and Frontiers Science Center for Nano-optoelectronics, School of Physics, Peking University, Beijing 100871, China
- Center for Applied Physics and Technology, HEDPS, Peking University, Beijing 100871, China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| |
Collapse
|
4
|
Villaseco Arribas E, Maitra NT, Agostini F. Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states. J Chem Phys 2024; 160:054102. [PMID: 38310471 DOI: 10.1063/5.0186984] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2023] [Accepted: 01/11/2024] [Indexed: 02/05/2024] Open
Abstract
Advances in coherent light sources and development of pump-probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron-nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. This work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.
Collapse
Affiliation(s)
- Evaristo Villaseco Arribas
- Department of Physics, Rutgers University, Newark, New Jersey 07102, USA
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark, New Jersey 07102, USA
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
5
|
Ofstad BS, Wibowo-Teale M, Kristiansen HE, Aurbakken E, Kitsaras MP, Schøyen ØS, Hauge E, Irons TJP, Kvaal S, Stopkowicz S, Wibowo-Teale AM, Pedersen TB. Magnetic optical rotation from real-time simulations in finite magnetic fields. J Chem Phys 2023; 159:204109. [PMID: 38018753 DOI: 10.1063/5.0171927] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Accepted: 11/01/2023] [Indexed: 11/30/2023] Open
Abstract
We present a numerical approach to magnetic optical rotation based on real-time time-dependent electronic-structure theory. Not relying on perturbation expansions in the magnetic field strength, the formulation allows us to test the range of validity of the linear relation between the rotation angle per unit path length and the magnetic field strength that was established empirically by Verdet 160 years ago. Results obtained from time-dependent coupled-cluster and time-dependent current density-functional theory are presented for the closed-shell molecules H2, HF, and CO in magnetic fields up to 55 kT at standard temperature and pressure conditions. We find that Verdet's linearity remains valid up to roughly 10-20 kT, above which significant deviations from linearity are observed. Among the three current density-functional approximations tested in this work, the current-dependent Tao-Perdew-Staroverov-Scuseria hybrid functional performs the best in comparison with time-dependent coupled-cluster singles and doubles results for the magnetic optical rotation.
Collapse
Affiliation(s)
- Benedicte Sverdrup Ofstad
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| | - Meilani Wibowo-Teale
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - Håkon Emil Kristiansen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| | - Einar Aurbakken
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| | - Marios Petros Kitsaras
- Physical and Theoretical Chemistry, Saarland University, Campus B2.2, 66123 Saarbruecken, Germany
| | | | - Eirill Hauge
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
- Department of Numerical Analysis and Scientific Computing, Simula Research Laboratory, 0164 Oslo, Norway
| | - Tom J P Irons
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - Simen Kvaal
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| | - Stella Stopkowicz
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
- Physical and Theoretical Chemistry, Saarland University, Campus B2.2, 66123 Saarbruecken, Germany
| | - Andrew M Wibowo-Teale
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - Thomas Bondo Pedersen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| |
Collapse
|
6
|
Yang YN, Chen SQ, Zhang ZH, Jiang H, Chen M, Li Y, He F. Harmonic Suppression Induced by Three-Electron Dynamics of Li in Strong Laser Fields. PHYSICAL REVIEW LETTERS 2023; 131:183201. [PMID: 37977615 DOI: 10.1103/physrevlett.131.183201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/24/2022] [Revised: 02/20/2023] [Accepted: 10/03/2023] [Indexed: 11/19/2023]
Abstract
We build a model to elucidate the high harmonic generation in combined EUV and midinfrared laser fields by embodying the spin-resolved three-electron dynamics. The EUV pulse ionizes an inner-shell electron, and the midinfrared laser drives the photoelectron and steers the electron-ion rescattering. Depending on the spin of the photoelectron, the residual ion including two bound electrons can be either in a single spin configuration or in a coherent superposition of different spin configurations. In the latter case, the two electrons in the ion swap their orbits, leading to a deep valley in the harmonic spectrum. The model results agree with the time-dependent Schrödinger equation simulations including three active electrons. The intriguing picture explored in this work is fundamentally distinguished from all reported scenarios relied on spin-orbit coupling, but originates from the exchanges asymmetry of two-electron wave functions.
Collapse
Affiliation(s)
- Yu-Ning Yang
- Key Laboratory for Laser Plasmas (Ministry of Education) and School of Physics and Astronomy, Collaborative innovation center for IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China
| | - Su-Qi Chen
- Key Laboratory for Laser Plasmas (Ministry of Education) and School of Physics and Astronomy, Collaborative innovation center for IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China
| | - Zhao-Han Zhang
- Key Laboratory for Laser Plasmas (Ministry of Education) and School of Physics and Astronomy, Collaborative innovation center for IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China
| | - Hui Jiang
- Key Laboratory for Laser Plasmas (Ministry of Education) and School of Physics and Astronomy, Collaborative innovation center for IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China
| | - Min Chen
- Key Laboratory for Laser Plasmas (Ministry of Education) and School of Physics and Astronomy, Collaborative innovation center for IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China
| | - Yang Li
- Key Laboratory for Laser Plasmas (Ministry of Education) and School of Physics and Astronomy, Collaborative innovation center for IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China
| | - Feng He
- Key Laboratory for Laser Plasmas (Ministry of Education) and School of Physics and Astronomy, Collaborative innovation center for IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China
- CAS Center for Excellence in Ultra-intense Laser Science, Shanghai 201800, China
| |
Collapse
|
7
|
Di Felice R, Mayes ML, Richard RM, Williams-Young DB, Chan GKL, de Jong WA, Govind N, Head-Gordon M, Hermes MR, Kowalski K, Li X, Lischka H, Mueller KT, Mutlu E, Niklasson AMN, Pederson MR, Peng B, Shepard R, Valeev EF, van Schilfgaarde M, Vlaisavljevich B, Windus TL, Xantheas SS, Zhang X, Zimmerman PM. A Perspective on Sustainable Computational Chemistry Software Development and Integration. J Chem Theory Comput 2023; 19:7056-7076. [PMID: 37769271 PMCID: PMC10601486 DOI: 10.1021/acs.jctc.3c00419] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Indexed: 09/30/2023]
Abstract
The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.
Collapse
Affiliation(s)
- Rosa Di Felice
- Departments
of Physics and Astronomy and Quantitative and Computational Biology, University of Southern California, Los Angeles, California 90089, United States
- CNR-NANO
Modena, Modena 41125, Italy
| | - Maricris L. Mayes
- Department
of Chemistry and Biochemistry, University
of Massachusetts Dartmouth, North Dartmouth, Massachusetts 02747, United States
| | | | | | - Garnet Kin-Lic Chan
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Wibe A. de Jong
- Lawrence
Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Niranjan Govind
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Martin Head-Gordon
- Pitzer Center
for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Karol Kowalski
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Xiaosong Li
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409, United States
| | - Karl T. Mueller
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Erdal Mutlu
- Advanced
Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Anders M. N. Niklasson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Mark R. Pederson
- Department
of Physics, The University of Texas at El
Paso, El Paso, Texas 79968, United States
| | - Bo Peng
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, United States
| | - Edward F. Valeev
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | | | - Bess Vlaisavljevich
- Department
of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Theresa L. Windus
- Department
of Chemistry, Iowa State University and
Ames Laboratory, Ames, Iowa 50011, United States
| | - Sotiris S. Xantheas
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
- Advanced
Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Xing Zhang
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Paul M. Zimmerman
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
8
|
Herbert JM, Zhu Y, Alam B, Ojha AK. Time-Dependent Density Functional Theory for X-ray Absorption Spectra: Comparing the Real-Time Approach to Linear Response. J Chem Theory Comput 2023; 19:6745-6760. [PMID: 37708349 DOI: 10.1021/acs.jctc.3c00673] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/16/2023]
Abstract
We simulate X-ray absorption spectra at elemental K-edges using time-dependent density functional theory (TDDFT) in both its conventional linear-response implementation and its explicitly time-dependent or "real-time" formulation. Real-time TDDFT simulations enable broadband spectra calculations without the need to invoke frozen occupied orbitals ("core/valence separation"), but we find that these spectra are often contaminated by transitions to the continuum that originate from lower-energy core and semicore orbitals. This problem becomes acute in triple-ζ basis sets, although it is sometimes sidestepped in double-ζ basis sets. Transitions to the continuum acquire surprisingly large dipole oscillator strengths, leading to spectra that are difficult to interpret. Meaningful spectra can be recovered by means of a filtering technique that decomposes the spectrum into contributions from individual occupied orbitals, and the same procedure can be used to separate L- and K-edge spectra arising from different elements within a given molecule. In contrast, conventional linear-response TDDFT requires core/valence separation but is free of these artifacts. It is also significantly more efficient than the real-time approach, even when hundreds of individual states are needed to reproduce near-edge absorption features and even when Padé approximants are used to reduce the real-time simulations to just 2-4 fs of time propagation. Despite the cost, the real-time approach may be useful to examine the validity of the core/valence separation approximation.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
- Chemical Physics Graduate Program, The Ohio State University, Columbus, Ohio 43210, United States
| | - Ying Zhu
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
- Chemical Physics Graduate Program, The Ohio State University, Columbus, Ohio 43210, United States
| | - Bushra Alam
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Avik Kumar Ojha
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
9
|
Calegari F, Martin F. Open questions in attochemistry. Commun Chem 2023; 6:184. [PMID: 37666969 PMCID: PMC10477171 DOI: 10.1038/s42004-023-00989-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2023] [Accepted: 08/17/2023] [Indexed: 09/06/2023] Open
Affiliation(s)
- Francesca Calegari
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607, Hamburg, Germany.
- The Hamburg Centre for Ultrafast Imaging, Universität Hamburg, Luruper Chaussee 149, 22761, Hamburg, Germany.
- Institut für Experimentalphysik, Universität Hamburg, Luruper Chaussee 149, 22761, Hamburg, Germany.
| | - Fernando Martin
- Instituto Madrileño de Estudios Avanzados en Nanociencia, Cantoblanco, 28049, Madrid, Spain.
- Departamento de Química, Universidad Autónoma de Madrid, 28049, Madrid, Spain.
| |
Collapse
|
10
|
Schlegel HB. Charge Migration in HCCI Cations Probed by Strong Field Ionization: Time-Dependent Configuration Interaction and Vibrational Wavepacket Simulations. J Phys Chem A 2023; 127:6040-6050. [PMID: 37459461 DOI: 10.1021/acs.jpca.3c02667] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
Strong field ionization of neutral iodoacetylene (HCCI) can produce a coherent superposition of the X and A cations and results in charge migration between the CC π orbital and the iodine π-type lone pair. This charge migration causes oscillations in the rate of strong field ionization of the cation to the dication that can be monitored using intense few-cycle probe pulses. The dynamics and strong field ionization of the coherent superposition the X and A states of HCCI+ have been modeled by time-dependent configuration interaction (TDCI) simulations. When the nuclei are allowed to move, the electronic wavefunctions need to be multiplied by vibrational wavefunctions. Nuclear motion has been modeled by vibrational packets moving on quadratic approximations to the potential energy surfaces for the X and A states of the cation. The overlap of the vibrational wavepackets decays in about 10-15 fs. Consequently, the oscillations in the strong field ionization decay on the same time scale. A revival of the vibrational overlap and in the oscillations of the strong field ionization is seen at 60-110 fs. TDCI simulations show that the decay and revival of the charge migration can be monitored by strong field ionization with intense 2- and 4-cycle linearly polarized 800 nm pulses. The revival is also seen with 7-cycle pulses.
Collapse
Affiliation(s)
- H Bernhard Schlegel
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| |
Collapse
|
11
|
Danilov D, Jenkins AJ, Bearpark MJ, Worth GA, Robb MA. Coherent Mixing of Singlet and Triplet States in Acrolein and Ketene: A Computational Strategy for Simulating the Electron-Nuclear Dynamics of Intersystem Crossing. J Phys Chem Lett 2023:6127-6134. [PMID: 37364275 DOI: 10.1021/acs.jpclett.3c01187] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/28/2023]
Abstract
We present a theoretical study of intersystem crossing (ISC) in acrolein and ketene with the Ehrenfest method that can describe a superposition of singlet and triplet states. Our simulations illustrate a new mechanistic effect of ISC, namely, that a superposition of singlets and triplets yields nonadiabatic dynamics characteristic of that superposition rather than the constituent state potential energy surfaces. This effect is particularly significant in ketene, where mixing of singlet and triplet states along the approach to a singlet/singlet conical intersection occurs, with the spin-orbit coupling (SOC) remaining small throughout. In both cases, the effects require many recrossings of the singlet/triplet state crossing seam, consistent with the textbook treatment of ISC.
Collapse
Affiliation(s)
- Don Danilov
- Department of Chemistry, Imperial College London, Molecular Sciences Research Hub, 82 Wood Lane, W12 0BZ London, United Kingdom
| | - Andrew J Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Michael J Bearpark
- Department of Chemistry, Imperial College London, Molecular Sciences Research Hub, 82 Wood Lane, W12 0BZ London, United Kingdom
| | - Graham A Worth
- Department of Chemistry, University College London, 20 Gordon St., WC1H 0AJ London, United Kingdom
| | - Michael A Robb
- Department of Chemistry, Imperial College London, Molecular Sciences Research Hub, 82 Wood Lane, W12 0BZ London, United Kingdom
| |
Collapse
|
12
|
Koushki AM. High-order harmonic generation from aligned HCN molecules under orthogonally and linearly polarized two-color laser fields. J Mol Model 2023; 29:137. [PMID: 37052706 DOI: 10.1007/s00894-023-05544-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2022] [Accepted: 04/03/2023] [Indexed: 04/14/2023]
Abstract
Molecular high-order harmonic generation and molecular orbital ionization probabilities are calculated under orthogonally and linearly polarized two-color laser fields. When a second-harmonic field is applied, the high-order harmonics generated under the linearly polarized two-color laser fields in the antiparallel case are stronger than those generated in the orthogonal polarization case and even stronger than those of the parallel polarization case. The results show that ionization probabilities of various orbitals and harmonic orders are dependent on spatial symmetry of molecular orbitals. It is found that the ionization of low-lying Kohn-Sham molecular orbitals contributes significantly to the ionization and molecular high-order harmonic generation processes. The ionization probability maximum occurs when molecular orbital densities are maximum in the direction of laser field polarization. Furthermore, we show that the degeneracy of π orbitals is broken when the laser-molecule alignment angle deviates from the field axis. Accordingly, we indicated one component of the π orbital is effectively contributed to the ionization and high-order harmonic generation processes. Finally, to confirm the recollision model in the high-order harmonic generation, the quantum time-frequency analysis is used to extract electron paths information on subcycle time scales.
Collapse
Affiliation(s)
- A M Koushki
- Photonics and Quantum Technologies Research School, Nuclear Science and Technology Research Institute, North Kargar Avenue, P. O. Box 11365-8489, Tehran, Iran.
| |
Collapse
|
13
|
Yao HB, Qu QW, Zhang ZH, Wang JW, Gao J, Hu CX, Li H, Wu J, He F. Multiphoton Ionization Reduction of Atoms in Two-Color Femtosecond Laser Fields. PHYSICAL REVIEW LETTERS 2023; 130:113201. [PMID: 37001077 DOI: 10.1103/physrevlett.130.113201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/04/2022] [Accepted: 02/23/2023] [Indexed: 06/19/2023]
Abstract
We report the ionization reduction of atoms in two-color femtosecond laser fields in this joint theoretical-experimental study. For the multiphoton ionization of atoms using a 400 nm laser pulse, the ionization probability is reduced if another relatively weak 800 nm laser pulse is overlapped. Such ionization reduction consistently occurs regardless of the relative phase between the two pulses. The time-dependent Schrödinger equation simulation results indicate that with the assisted 800 nm photons the electron can be launched to Rydberg states with large angular quantum numbers, which stand off the nuclei and thus are hard to be freed in the multiphoton regime. This mechanism works for hydrogen, helium, and probably some other atoms if two-color laser fields are properly tuned.
Collapse
Affiliation(s)
- Hong-Bin Yao
- Key Laboratory for Laser Plasmas (Ministry of Education) and School of Physics and Astronomy, Collaborative innovation center for IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China
- Key Laboratory of New Energy and Materials Research of Xinjiang Education Department, Xinjiang Institute of Engineering, Urumqi 830091, China
| | - Qi-Wen Qu
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
| | - Zhao-Han Zhang
- Key Laboratory for Laser Plasmas (Ministry of Education) and School of Physics and Astronomy, Collaborative innovation center for IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China
| | - Jia-Wei Wang
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
| | - Jian Gao
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
- Chongqing Key Laboratory of Precision Optics, Chongqing Institute of East China Normal University, Chongqing 401121, China
| | - Chen-Xi Hu
- Key Laboratory for Laser Plasmas (Ministry of Education) and School of Physics and Astronomy, Collaborative innovation center for IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China
| | - Hui Li
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
| | - Jian Wu
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
- Chongqing Key Laboratory of Precision Optics, Chongqing Institute of East China Normal University, Chongqing 401121, China
- CAS Center for Excellence in Ultra-intense Laser Science, Shanghai 201800, China
| | - Feng He
- Key Laboratory for Laser Plasmas (Ministry of Education) and School of Physics and Astronomy, Collaborative innovation center for IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China
- CAS Center for Excellence in Ultra-intense Laser Science, Shanghai 201800, China
| |
Collapse
|
14
|
Gelin MF, Chen L, Domcke W. Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals. Chem Rev 2022; 122:17339-17396. [PMID: 36278801 DOI: 10.1021/acs.chemrev.2c00329] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Lipeng Chen
- Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Strasse 38, D-01187 Dresden, Germany
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, D-85747 Garching,Germany
| |
Collapse
|
15
|
Zhu Y, Herbert JM. High harmonic spectra computed using time-dependent Kohn-Sham theory with Gaussian orbitals and a complex absorbing potential. J Chem Phys 2022; 156:204123. [PMID: 35649850 DOI: 10.1063/5.0079910] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
High harmonic spectra for H2 and H2 + are simulated by solving the time-dependent Kohn-Sham equation in the presence of a strong laser field using an atom-centered Gaussian representation of the density and a complex absorbing potential. The latter serves to mitigate artifacts associated with the finite extent of the basis functions, including spurious reflection of the outgoing electronic wave packet. Interference between the outgoing and reflected waves manifests as peak broadening in the spectrum as well as the appearance of spurious high-energy peaks after the harmonic progression has terminated. We demonstrate that well-resolved spectra can be obtained through the use of an atom-centered absorbing potential. As compared to grid-based algorithms, the present approach is more readily extensible to larger molecules.
Collapse
Affiliation(s)
- Ying Zhu
- Department of Chemistry and Biochemistry, and Chemical Physics Graduate Program, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M Herbert
- Department of Chemistry and Biochemistry, and Chemical Physics Graduate Program, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
16
|
Schlegel HB, Hoerner P, Li W. Ionization of HCCI Neutral and Cations by Strong Laser Fields Simulated With Time Dependent Configuration Interaction. Front Chem 2022; 10:866137. [PMID: 35548678 PMCID: PMC9081608 DOI: 10.3389/fchem.2022.866137] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2022] [Accepted: 03/24/2022] [Indexed: 12/05/2022] Open
Abstract
Strong field ionization of neutral iodoacetylene (HCCI) can produce a coherent superposition of the X and A cations. This superposition results in charge migration between the CC π orbital and the iodine π-type lone pair which can be monitored by strong field ionization with short, intense probe pulses. Strong field ionization of the X and A states of HCCI cation was simulated with time-dependent configuration interaction using singly ionized configurations and singly excited, singly ionized configurations (TD-CISD-IP) and an absorbing boundary. Studies with static fields were used to obtain the 3-dimensional angular dependence of instantaneous ionization rates by strong fields and the orbitals involved in producing the cations and dications. The frequency of charge oscillation is determined by the energy separation of the X and A states; this separation can change depending on the direction and strength of the field. Furthermore, fields along the molecular axis can cause extensive mixing between the field-free X and A configurations. For coherent superpositions of the X and A states, the charge oscillations are characterized by two frequencies-the driving frequency of the laser field of the probe pulse and the intrinsic frequency due to the energy separation between the X and A states. For linear and circularly polarized pulses, the ionization rates show marked differences that depend on the polarization direction of the pulse, the carrier envelope phase and initial phase of the superposition. Varying the initial phase of the superposition at the beginning of the probe pulse is analogous to changing the delay between the pump and probe pulses. The charge oscillation in the coherent superposition of the X and A states results in maxima and minima in the ionization yield as a function of the superposition phase.
Collapse
|
17
|
Coccia E, Luppi E. Time-dependent ab initioapproaches for high-harmonic generation spectroscopy. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 34:073001. [PMID: 34731835 DOI: 10.1088/1361-648x/ac3608] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Accepted: 11/03/2021] [Indexed: 06/13/2023]
Abstract
High-harmonic generation (HHG) is a nonlinear physical process used for the production of ultrashort pulses in XUV region, which are then used for investigating ultrafast phenomena in time-resolved spectroscopies. Moreover, HHG signal itself encodes information on electronic structure and dynamics of the target, possibly coupled to the nuclear degrees of freedom. Investigating HHG signal leads to HHG spectroscopy, which is applied to atoms, molecules, solids and recently also to liquids. Analysing the number of generated harmonics, their intensity and shape gives a detailed insight of, e.g., ionisation and recombination channels occurring in the strong-field dynamics. A number of valuable theoretical models has been developed over the years to explain and interpret HHG features, with the three-step model being the most known one. Originally, these models neglect the complexity of the propagating electronic wavefunction, by only using an approximated formulation of ground and continuum states. Many effects unravelled by HHG spectroscopy are instead due to electron correlation effects, quantum interference, and Rydberg-state contributions, which are all properly captured by anab initioelectronic-structure approach. In this review we have collected recent advances in modelling HHG by means ofab initiotime-dependent approaches relying on the propagation of the time-dependent Schrödinger equation (or derived equations) in presence of a very intense electromagnetic field. We limit ourselves to gas-phase atomic and molecular targets, and to solids. We focus on the various levels of theory employed for describing the electronic structure of the target, coupled with strong-field dynamics and ionisation approaches, and on the basis used to represent electronic states. Selected applications and perspectives for future developments are also given.
Collapse
Affiliation(s)
- Emanuele Coccia
- Dipartimento di Scienze Chimiche e Farmaceutiche, University of Trieste, via Giorgieri 1, 34127 Trieste, Italy
| | - Eleonora Luppi
- Sorbonne Universités, UPMC Univ Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005 Paris, France
- CNRS, UMR 7616, Laboratoire de Chimie Théorique, F-75005 Paris, France
| |
Collapse
|
18
|
Ghomashi B, Reiff R, Becker A. Coherence in macroscopic high harmonic generation for spatial focal phase distributions of monochromatic and broadband Gaussian laser pulses. OPTICS EXPRESS 2021; 29:40146-40160. [PMID: 34809362 DOI: 10.1364/oe.444317] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/26/2021] [Accepted: 11/08/2021] [Indexed: 06/13/2023]
Abstract
Efficient application of ultrafast laser sources from high harmonic generation requires an understanding of how the spectrum can be controlled - the extent of the highest harmonics and the strength and cleanness of the harmonic lines. We study one important aspect in the coherent build-up of macroscopic high-order harmonic generation, namely the impact of different phase distributions in the focal area on the features of the generated radiation. Specifically, we compare the high harmonic signals for the commonly-used Gouy distribution of a monochromatic beam with those for the phase distribution of a short broadband Gaussian pulse. To this end, we apply a theoretical model in which the microscopic yields are obtained via interpolation of results of the time-dependent Schrödinger equation, which are then used in an individual-emitter approach to determine the macroscopic signals. Regions of poor and good coherent build-up as a function of the position of the gas jet are identified using measures for the strength of the harmonic lines and for the impact of off-harmonic radiation. While the largest extent of the spectra as well as the strongest contribution of off-harmonic radiation is found for positioning the gas jet after the focus for both distributions, the relative strength of the harmonics is overall weaker for the short Gaussian pulse distribution and the spectra differ for a gas jet positioned at the focus. These differences are mainly caused by the additional dependence of the focal phase in the transverse direction for the short Gaussian pulse distribution.
Collapse
|
19
|
Wang C, Ulusoy IS, Aebersold LE, Wilson AK. Multi-configuration electron-nuclear dynamics: An open-shell approach. J Chem Phys 2021; 155:154103. [PMID: 34686063 DOI: 10.1063/5.0063478] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The multi-configuration electron-nuclear dynamics for open shell systems with a spin-unrestricted formalism is described. The mean fields are evaluated using second-order reduced density matrices for electronic and nuclear degrees of freedom. Applications to light-element diatomics including equilibrium geometries, electronic energies, dipole moments, and absorption spectra are presented. The von Neumann entropies for different spin states of a LiH molecule are compared.
Collapse
Affiliation(s)
- Cong Wang
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Inga S Ulusoy
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Lucas E Aebersold
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| |
Collapse
|
20
|
Abstract
The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this "on-the-fly" approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. This approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH2O → H2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.
Collapse
Affiliation(s)
- H. Bernhard Schlegel
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| |
Collapse
|
21
|
Grobas Illobre P, Marsili M, Corni S, Stener M, Toffoli D, Coccia E. Time-Resolved Excited-State Analysis of Molecular Electron Dynamics by TDDFT and Bethe-Salpeter Equation Formalisms. J Chem Theory Comput 2021; 17:6314-6329. [PMID: 34486881 PMCID: PMC8515806 DOI: 10.1021/acs.jctc.1c00211] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2021] [Indexed: 12/16/2022]
Abstract
In this work, a theoretical and computational set of tools to study and analyze time-resolved electron dynamics in molecules, under the influence of one or more external pulses, is presented. By coupling electronic-structure methods with the resolution of the time-dependent Schrödinger equation, we developed and implemented the time-resolved induced density of the electronic wavepacket, the time-resolved formulation of the differential projection density of states (ΔPDOS), and of transition contribution map (TCM) to look at the single-electron orbital occupation and localization change in time. Moreover, to further quantify the possible charge transfer, we also defined the energy-integrated ΔPDOS and the fragment-projected TCM. We have used time-dependent density-functional theory (TDDFT), as implemented in ADF software, and the Bethe-Salpeter equation, as provided by MolGW package, for the description of the electronic excited states. This suite of postprocessing tools also provides the time evolution of the electronic states of the system of interest. To illustrate the usefulness of these postprocessing tools, excited-state populations have been computed for HBDI (the chromophore of GFP) and DNQDI molecules interacting with a sequence of two pulses. Time-resolved descriptors have been applied to study the time-resolved electron dynamics of HBDI, DNQDI, LiCN (being a model system for dipole switching upon highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) electronic excitation), and Ag22. The computational analysis tools presented in this article can be employed to help the interpretation of fast and ultrafast spectroscopies on molecular, supramolecular, and composite systems.
Collapse
Affiliation(s)
- P. Grobas Illobre
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| | - M. Marsili
- Dipartimento
di Scienze Chimiche, Universitá di
Padova, via Marzolo 1, Padova 35131, Italy
| | - S. Corni
- Dipartimento
di Scienze Chimiche, Universitá di
Padova, via Marzolo 1, Padova 35131, Italy
- CNR
Istituto di Nanoscienze, via Campi 213/A, Modena 41125, Italy
| | - M. Stener
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| | - D. Toffoli
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| | - E. Coccia
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| |
Collapse
|
22
|
Aebersold LE, Ulusoy IS, Wilson AK. Electron-nuclear quantum dynamics of diatomic molecules: nonadiabatic signatures in molecular spectra. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1988743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Lucas E. Aebersold
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| | - Inga S. Ulusoy
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
- Scientific Software Center, Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany
| | - Angela K. Wilson
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| |
Collapse
|
23
|
Hoerner P, Li W, Schlegel HB. Sequential double ionization of molecules by strong laser fields simulated with time-dependent configuration interaction. J Chem Phys 2021; 155:114103. [PMID: 34551546 DOI: 10.1063/5.0060365] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new time-dependent configuration interaction method has been developed for simulating strong field sequential double ionization of molecular systems. Ionization of the neutral is simulated by time-dependent configuration interaction with single excitations (TD-CIS) and an absorbing boundary. At each time step, the ionized part of the wavefunction from the TD-CIS calculation is transferred to a second time-dependent configuration interaction simulation for ionization of the cation to the dication. The second simulation uses a CISD-IP wavefunction that consists of singly ionized configurations and singly excited, singly ionized configurations (TD-CISD-IP). The transfer between the TD-CIS and TD-CISD-IP simulations is accomplished by partitioning the first ionization rate into contributions from individual orbitals or by singular value decomposition of the absorbed wavefunction. Sequential double ionization simulations have been carried out for HBr in five cycle 800 nm linearly polarized pulses and HI (with spin-orbit coupling) in four cycle 800 nm circularly polarized pulses, with intensities chosen so that the population of the neutral was depleted by the mid-pulse. The singular value decomposition of the cation produced by the first ionization is dominated by a single component for the two orientations considered. The population of the cation rises and then falls as it is ionized to the dication. Depending on the pulse shape and field strength, the ionization of the cation to the dication can continue for several half cycles. For HI with circularly polarized light, the rates for both the first and second ionization peak when the electric field is aligned with the pπ orbital.
Collapse
Affiliation(s)
- Paul Hoerner
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| | - Wen Li
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| | - H Bernhard Schlegel
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| |
Collapse
|
24
|
Merritt ICD, Jacquemin D, Vacher M. Attochemistry: Is Controlling Electrons the Future of Photochemistry? J Phys Chem Lett 2021; 12:8404-8415. [PMID: 34436903 DOI: 10.1021/acs.jpclett.1c02016] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Controlling matter with light has always been a great challenge, leading to the ever-expanding field of photochemistry. In addition, since the first generation of light pulses of attosecond (1 as = 10-18 s) duration, a great deal of effort has been devoted to observing and controlling electrons on their intrinsic time scale. Because of their short duration, attosecond pulses have a large spectral bandwidth populating several electronically excited states in a coherent manner, i.e., an electronic wavepacket. Because of interference, such a wavepacket has a new electronic distribution implying a potentially different and totally new reactivity as compared to traditional photochemistry, leading to the novel concept of "attochemistry". This nascent field requires the support of theory right from the start. In this Perspective, we discuss the opportunities offered by attochemistry, the related challenges, and the current and future state-of-the-art developments in theoretical chemistry needed to model it accurately.
Collapse
Affiliation(s)
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM, UMR 6230, F-44000 Nantes, France
| | - Morgane Vacher
- Université de Nantes, CNRS, CEISAM, UMR 6230, F-44000 Nantes, France
| |
Collapse
|
25
|
Baiardi A. Electron Dynamics with the Time-Dependent Density Matrix Renormalization Group. J Chem Theory Comput 2021; 17:3320-3334. [PMID: 34043347 DOI: 10.1021/acs.jctc.0c01048] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
In this work, we simulate the electron dynamics in molecular systems with the time-dependent density matrix renormalization group (TD-DMRG) algorithm. We leverage the generality of the so-called tangent-space TD-DMRG formulation and design a computational framework in which the dynamics is driven by the exact nonrelativistic electronic Hamiltonian. We show that by parametrizing the wave function as a matrix product state, we can accurately simulate the dynamics of systems including up to 20 electrons and 32 orbitals. We apply the TD-DMRG algorithm to three problems that are hardly targeted by time-independent methods: the calculation of molecular (hyper)polarizabilities, the simulation of electronic absorption spectra, and the study of ultrafast ionization dynamics.
Collapse
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, Zürich 8093, Switzerland
| |
Collapse
|
26
|
Lee MK, Hoerner P, Li W, Schlegel HB. Effect of spin-orbit coupling on strong field ionization simulated with time-dependent configuration interaction. J Chem Phys 2020; 153:244109. [PMID: 33380070 DOI: 10.1063/5.0034807] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Time-dependent configuration interaction with a complex absorbing potential has been used to simulate strong field ionization by intense laser fields. Because spin-orbit coupling changes the energies of the ground and excited states, it can affect the strong field ionization rate for molecules containing heavy atoms. Configuration interaction with single excitations (CIS) has been employed for strong field ionization of closed shell systems. Single and double excitation configuration interaction with ionization (CISD-IP) has been used to treat ionization of degenerate states of cations on an equal footing. The CISD-IP wavefunction consists of ionizing single (one hole) and double (two hole/one particle) excitations from the neutral atom. Spin-orbit coupling has been implemented using an effective one electron spin-orbit coupling operator. The effective nuclear charge in the spin-orbit coupling operator has been optimized for Ar+, Kr+, Xe+, HX+ (X = Cl, Br, and I). Spin-orbit effects on angular dependence of the strong field ionization have been studied for HX and HX+. The effects of spin-orbit coupling are largest for ionization from the π orbitals of HX+. In a static field, oscillations are seen between the 2Π3/2 and 2Π1/2 states of HX+. For ionization of HX+ by a two cycle circularly polarized pulse, a single peak is seen when the maximum in the carrier envelope is perpendicular to the molecular axis and two peaks are seen when it is parallel to the axis. This is the result of the greater ionization rate for the π orbitals than for the σ orbitals.
Collapse
Affiliation(s)
- Mi Kyung Lee
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| | - Paul Hoerner
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| | - Wen Li
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| | - H Bernhard Schlegel
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| |
Collapse
|
27
|
Jahnke T, Hergenhahn U, Winter B, Dörner R, Frühling U, Demekhin PV, Gokhberg K, Cederbaum LS, Ehresmann A, Knie A, Dreuw A. Interatomic and Intermolecular Coulombic Decay. Chem Rev 2020; 120:11295-11369. [PMID: 33035051 PMCID: PMC7596762 DOI: 10.1021/acs.chemrev.0c00106] [Citation(s) in RCA: 65] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Interatomic or intermolecular Coulombic decay (ICD) is a nonlocal electronic decay mechanism occurring in weakly bound matter. In an ICD process, energy released by electronic relaxation of an excited atom or molecule leads to ionization of a neighboring one via Coulombic electron interactions. ICD has been predicted theoretically in the mid nineties of the last century, and its existence has been confirmed experimentally approximately ten years later. Since then, a number of fundamental and applied aspects have been studied in this quickly growing field of research. This review provides an introduction to ICD and draws the connection to related energy transfer and ionization processes. The theoretical approaches for the description of ICD as well as the experimental techniques developed and employed for its investigation are described. The existing body of literature on experimental and theoretical studies of ICD processes in different atomic and molecular systems is reviewed.
Collapse
Affiliation(s)
- Till Jahnke
- Institut für Kernphysik, Goethe Universität, Max-von-Laue-Str. 1, 60438 Frankfurt, Germany
| | - Uwe Hergenhahn
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany.,Max Planck Institute for Plasma Physics, Wendelsteinstr. 1, 17491 Greifswald, Germany.,Leibniz Institute of Surface Engineering (IOM), 04318 Leipzig, Germany
| | - Bernd Winter
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Reinhard Dörner
- Institut für Kernphysik, Goethe Universität, Max-von-Laue-Str. 1, 60438 Frankfurt, Germany
| | - Ulrike Frühling
- Institut für Experimentalphysik and Center for Ultrafast Imaging, Universität Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Philipp V Demekhin
- Institut für Physik und CINSaT, Universität Kassel, Heinrich-Plett-Str. 40, 34132 Kassel, Germany
| | - Kirill Gokhberg
- Physical-Chemistry Institute, Ruprecht-Karls University, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Lorenz S Cederbaum
- Physical-Chemistry Institute, Ruprecht-Karls University, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Arno Ehresmann
- Institut für Physik und CINSaT, Universität Kassel, Heinrich-Plett-Str. 40, 34132 Kassel, Germany
| | - André Knie
- Institut für Physik und CINSaT, Universität Kassel, Heinrich-Plett-Str. 40, 34132 Kassel, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
28
|
Quantum state holography to reconstruct the molecular wave packet using an attosecond XUV-XUV pump-probe technique. Sci Rep 2020; 10:12981. [PMID: 32737413 PMCID: PMC7395139 DOI: 10.1038/s41598-020-69733-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2020] [Accepted: 07/15/2020] [Indexed: 11/08/2022] Open
Abstract
An attosecond molecular interferometer is proposed by using a XUV-XUV pump-probe scheme. The interferograms resulting in the photoelectron distributions enable the full reconstruction of the molecular wave packet associated to excited states using a quantum state holographic approach that, to our knowledge, has only been proposed for simple atomic targets combining attosecond XUV pulses with IR light. In contrast with existing works, we investigate schemes where one- and two-photon absorption paths contribute to ionize the hydrogen molecule and show that it is possible to retrieve the excitation dynamics even when imprinted in a minority channel. Furthermore, we provide a systematic analysis of the time-frequency maps that reveal the distinct, but tightly coupled, motion of electrons and nuclei.
Collapse
|
29
|
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
|
30
|
Cabral Tenorio BN, Chaer Nascimento MA, Rocha AB, Coriani S. Lanczos-based equation-of-motion coupled-cluster singles-and-doubles approach to the total photoionization cross section of valence excited states. J Chem Phys 2019; 151:184106. [DOI: 10.1063/1.5125125] [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)
- Bruno Nunes Cabral Tenorio
- Universidade Federal do Rio de Janeiro, UFRJ, Instituto de Química, Av. Athos da Silveira Ramos, 149, Rio de Janeiro, RJ 21941-909, Brazil
| | - Marco Antonio Chaer Nascimento
- Universidade Federal do Rio de Janeiro, UFRJ, Instituto de Química, Av. Athos da Silveira Ramos, 149, Rio de Janeiro, RJ 21941-909, Brazil
| | - Alexandre Braga Rocha
- Universidade Federal do Rio de Janeiro, UFRJ, Instituto de Química, Av. Athos da Silveira Ramos, 149, Rio de Janeiro, RJ 21941-909, Brazil
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet, Building 207, DK-2800 Kongens Lyngby, Denmark
| |
Collapse
|