1
|
Zhang Q, Shao X, Li W, Mi W, Pavanello M, Akimov AV. Nonadiabatic molecular dynamics with subsystem density functional theory: application to crystalline pentacene. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2024; 36:385901. [PMID: 38866023 DOI: 10.1088/1361-648x/ad577d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2024] [Accepted: 06/12/2024] [Indexed: 06/14/2024]
Abstract
In this work, we report the development and assessment of the nonadiabatic molecular dynamics approach with the electronic structure calculations based on the linearly scaling subsystem density functional method. The approach is implemented in an open-source embedded Quantum Espresso/Libra software specially designed for nonadiabatic dynamics simulations in extended systems. As proof of the applicability of this method to large condensed-matter systems, we examine the dynamics of nonradiative relaxation of excess excitation energy in pentacene crystals with the simulation supercells containing more than 600 atoms. We find that increased structural disorder observed in larger supercell models induces larger nonadiabatic couplings of electronic states and accelerates the relaxation dynamics of excited states. We conduct a comparative analysis of several quantum-classical trajectory surface hopping schemes, including two new methods proposed in this work (revised decoherence-induced surface hopping and instantaneous decoherence at frustrated hops). Most of the tested schemes suggest fast energy relaxation occurring with the timescales in the 0.7-2.0 ps range, but they significantly overestimate the ground state recovery rates. Only the modified simplified decay of mixing approach yields a notably slower relaxation timescales of 8-14 ps, with a significantly inhibited ground state recovery.
Collapse
Affiliation(s)
- Qingxin Zhang
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, NY 14260, United States of America
| | - Xuecheng Shao
- Department of Physics, Rutgers University, The State University of New Jersey, Newark, NJ 07102, United States of America
| | - Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, People's Republic of China
| | - Wenhui Mi
- Key Laboratory of Material Simulation Methods & Software of Ministry of Education, College of Physics, Jilin University, Changchun 130012, People's Republic of China
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, People's Republic of China
| | - Michele Pavanello
- Department of Physics, Rutgers University, The State University of New Jersey, Newark, NJ 07102, United States of America
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, NY 14260, United States of America
| |
Collapse
|
2
|
Selenius E, Sigurdarson AE, Schmerwitz YLA, Levi G. Orbital-Optimized Versus Time-Dependent Density Functional Calculations of Intramolecular Charge Transfer Excited States. J Chem Theory Comput 2024; 20:3809-3822. [PMID: 38695313 DOI: 10.1021/acs.jctc.3c01319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2024]
Abstract
The performance of time-independent, orbital-optimized calculations of excited states is assessed with respect to charge transfer excitations in organic molecules in comparison to the linear-response time-dependent density functional theory (TD-DFT) approach. A direct optimization method to converge on saddle points of the electronic energy surface is used to carry out calculations with the local density approximation (LDA) and the generalized gradient approximation (GGA) functionals PBE and BLYP for a set of 27 excitations in 15 molecules. The time-independent approach is fully variational and provides a relaxed excited state electron density from which the extent of charge transfer is quantified. The TD-DFT calculations are generally found to provide larger charge transfer distances compared to the orbital-optimized calculations, even when including orbital relaxation effects with the Z-vector method. While the error on the excitation energy relative to theoretical best estimates is found to increase with the extent of charge transfer up to ca. -2 eV for TD-DFT, no correlation is observed for the orbital-optimized approach. The orbital-optimized calculations with the LDA and the GGA functionals provide a mean absolute error of ∼0.7 eV, outperforming TD-DFT with both local and global hybrid functionals for excitations with a long-range charge transfer character. Orbital-optimized calculations with the global hybrid functional B3LYP and the range-separated hybrid functional CAM-B3LYP on a selection of states with short- and long-range charge transfer indicate that inclusion of exact exchange has a small effect on the charge transfer distance, while it significantly improves the excitation energy, with the best-performing functional CAM-B3LYP providing an absolute error typically around 0.15 eV.
Collapse
Affiliation(s)
- Elli Selenius
- Science Institute of the University of Iceland, Reykjavík 107, Iceland
| | | | | | - Gianluca Levi
- Science Institute of the University of Iceland, Reykjavík 107, Iceland
| |
Collapse
|
3
|
Wang X, Krause P, Kirschbaum T, Palczynski K, Dzubiella J, Bande A. Photo-excited charge transfer from adamantane to electronic bound states in water. Phys Chem Chem Phys 2024; 26:8158-8176. [PMID: 38380443 DOI: 10.1039/d3cp04602h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2024]
Abstract
Aqueous nanodiamonds illuminated by UV light produce free solvated electrons, which may drive high-energy reduction reactions in water. However, the influence of water conformations on the excited-state electron-transfer mechanism are still under debate. In this work, we offer a theoretical study of charge-transfer states in adamantane-water structures obtained by linear-response time-dependent density-functional theory. Small water clusters with broken hydrogen bonds are found to efficiently bind the electron from adamantane. A distinction is made with respect to the nature of the water clusters: some bind the electron in a water cavity, others along a strong permanent total dipole. These two types of bound states are more strongly binding, the higher their electron affinity and their positive electrostatic potential, the latter being dominated by the energy of the lowest unoccupied molecular orbital of the isolated water clusters. Structural sampling in a thermal equilibrium at room temperature via molecular dynamics snapshots confirms under which conditions the underlying waters clusters can occur and verifies that broken hydrogen bonds in the water network close to adamantane can create traps for the solvated electron.
Collapse
Affiliation(s)
- Xiangfei Wang
- Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, 14109 Berlin, Germany.
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Arnimallee 22, 14195 Berlin, Germany
| | - Pascal Krause
- Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, 14109 Berlin, Germany.
| | - Thorren Kirschbaum
- Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, 14109 Berlin, Germany.
- Institute of Mathematics, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Karol Palczynski
- Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, 14109 Berlin, Germany.
| | - Joachim Dzubiella
- Applied Theoretical Physics - Computational Physics, Albert-Ludwigs-Universität Freiburg, Hermann-Herder Straße 3, 79104 Freiburg, Germany.
| | - Annika Bande
- Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, 14109 Berlin, Germany.
- Institute of Inorganic Chemistry, Leibniz University Hannover, Callinstr. 9, 30167 Hannover, Germany
| |
Collapse
|
4
|
Sigurdarson AE, Schmerwitz YLA, Tveiten DKV, Levi G, Jónsson H. Orbital-optimized density functional calculations of molecular Rydberg excited states with real space grid representation and self-interaction correction. J Chem Phys 2023; 159:214109. [PMID: 38047508 DOI: 10.1063/5.0179271] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2023] [Accepted: 10/30/2023] [Indexed: 12/05/2023] Open
Abstract
Density functional calculations of Rydberg excited states up to high energy are carried out for several molecules using an approach where the orbitals are variationally optimized by converging on saddle points on the electronic energy surface within a real space grid representation. Remarkably good agreement with experimental estimates of the excitation energy is obtained using the generalized gradient approximation (GGA) functional of Perdew, Burke, and Ernzerhof (PBE) when Perdew-Zunger self-interaction correction is applied in combination with complex-valued orbitals. Even without the correction, the PBE functional gives quite good results despite the fact that corresponding Rydberg virtual orbitals have positive energy in the ground state calculation. Results obtained using the Tao, Perdew, Staroverov, and Scuseria (TPSS) and r2SCAN meta-GGA functionals are also presented, but they do not provide a systematic improvement over the results from the uncorrected PBE functional. The grid representation combined with the projector augmented-wave approach gives a simpler and better representation of diffuse Rydberg orbitals than a linear combination of atomic orbitals with commonly used basis sets, the latter leading to an overestimation of the excitation energy due to confinement of the excited states.
Collapse
Affiliation(s)
- Alec E Sigurdarson
- Science Institute and Faculty of Physical Sciences, University of Iceland, Reykjavík, Iceland
| | - Yorick L A Schmerwitz
- Science Institute and Faculty of Physical Sciences, University of Iceland, Reykjavík, Iceland
| | - Dagrún K V Tveiten
- Science Institute and Faculty of Physical Sciences, University of Iceland, Reykjavík, Iceland
| | - Gianluca Levi
- Science Institute and Faculty of Physical Sciences, University of Iceland, Reykjavík, Iceland
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, Reykjavík, Iceland
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
5
|
Janesko BG. Core-Projected Hybrids Fix Systematic Errors in Time-Dependent Density Functional Theory Predicted Core-Electron Excitations. J Chem Theory Comput 2023. [PMID: 37437304 DOI: 10.1021/acs.jctc.3c00312] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/14/2023]
Abstract
Linear response time-dependent density functional theory (TDDFT) is widely applied to valence, Rydberg, and charge-transfer excitations but, in its current form, makes large errors for core-electron excitations. This work demonstrates that the admixture of nonlocal exact exchange in atomic core regions significantly improves TDDFT-predicted core excitations. Exact exchange admixture is accomplished using projected hybrid density functional theory [ J. Chem. Theory Comput. 2023, 19, 837-847]. Scalar relativistic TDDFT calculations using core-projected B3LYP accurately model core excitations of second-period elements C-F and third-period elements Si-Cl, without sacrificing performance for the relative shifts of core excitation energies. Predicted K-edge X-ray near absorption edge structure (XANES) of a series of sulfur standards highlight the value of this approach. Core-projected hybrids appear to be a practical solution to TDDFT's limitations for core excitations, in the way that long-range-corrected hybrids are a practical solution to TDDFT's limitations for Rydberg and charge-transfer excitations.
Collapse
Affiliation(s)
- Benjamin G Janesko
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, Texas 76129, United States
| |
Collapse
|
6
|
Katayama T, Choi TK, Khakhulin D, Dohn AO, Milne CJ, Vankó G, Németh Z, Lima FA, Szlachetko J, Sato T, Nozawa S, Adachi SI, Yabashi M, Penfold TJ, Gawelda W, Levi G. Atomic-scale observation of solvent reorganization influencing photoinduced structural dynamics in a copper complex photosensitizer. Chem Sci 2023; 14:2572-2584. [PMID: 36908966 PMCID: PMC9993854 DOI: 10.1039/d2sc06600a] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2022] [Accepted: 01/30/2023] [Indexed: 02/04/2023] Open
Abstract
Photochemical reactions in solution are governed by a complex interplay between transient intramolecular electronic and nuclear structural changes and accompanying solvent rearrangements. State-of-the-art time-resolved X-ray solution scattering has emerged in the last decade as a powerful technique to observe solute and solvent motions in real time. However, disentangling solute and solvent dynamics and how they mutually influence each other remains challenging. Here, we simultaneously measure femtosecond X-ray emission and scattering to track both the intramolecular and solvation structural dynamics following photoexcitation of a solvated copper photosensitizer. Quantitative analysis assisted by molecular dynamics simulations reveals a two-step ligand flattening strongly coupled to the solvent reorganization, which conventional optical methods could not discern. First, a ballistic flattening triggers coherent motions of surrounding acetonitrile molecules. In turn, the approach of acetonitrile molecules to the copper atom mediates the decay of intramolecular coherent vibrations and induces a further ligand flattening. These direct structural insights reveal that photoinduced solute and solvent motions can be intimately intertwined, explaining how the key initial steps of light harvesting are affected by the solvent on the atomic time and length scale. Ultimately, this work takes a step forward in understanding the microscopic mechanisms of the bidirectional influence between transient solvent reorganization and photoinduced solute structural dynamics.
Collapse
Affiliation(s)
- Tetsuo Katayama
- Japan Synchrotron Radiation Research Institute Kouto 1-1-1, Sayo Hyogo 679-5198 Japan.,RIKEN SPring-8 Center 1-1-1 Kouto, Sayo Hyogo 679-5148 Japan
| | - Tae-Kyu Choi
- XFEL Division, Pohang Accelerator Laboratory Jigok-ro 127-80 Pohang 37673 Republic of Korea
| | | | - Asmus O Dohn
- Science Institute, University of Iceland 107 Reykjavík Iceland .,DTU Physics, Technical University of Denmark Kongens Lyngby Denmark
| | | | - György Vankó
- Wigner Research Centre for Physics, Hungarian Academy of Sciences H-1525 Budapest Hungary
| | - Zoltán Németh
- Wigner Research Centre for Physics, Hungarian Academy of Sciences H-1525 Budapest Hungary
| | | | - Jakub Szlachetko
- SOLARIS National Synchrotron Radiation Centre, Jagiellonian University PL-30392 Kraków Poland
| | - Tokushi Sato
- European XFEL Holzkoppel 4, Schenefeld 22869 Germany
| | - Shunsuke Nozawa
- Institute of Materials Structure Science, High Energy Accelerator Research Organization (KEK) 1-1 Oho Tsukuba Ibaraki 305-0801 Japan.,Department of Materials Structure Science, School of High Energy Accelerator Science, The Graduate University for Advanced Studies 1-1 Oho Tsukuba Ibaraki 305-0801 Japan
| | - Shin-Ichi Adachi
- Institute of Materials Structure Science, High Energy Accelerator Research Organization (KEK) 1-1 Oho Tsukuba Ibaraki 305-0801 Japan.,Department of Materials Structure Science, School of High Energy Accelerator Science, The Graduate University for Advanced Studies 1-1 Oho Tsukuba Ibaraki 305-0801 Japan
| | - Makina Yabashi
- RIKEN SPring-8 Center 1-1-1 Kouto, Sayo Hyogo 679-5148 Japan
| | - Thomas J Penfold
- Chemistry-School of Natural and Environmental Sciences, Newcastle University Newcastle Upon-Tyne NE1 7RU UK
| | - Wojciech Gawelda
- Departamento de Química, Universidad Autónoma de Madrid, Campus Cantoblanco 28047 Madrid Spain.,IMDEA-Nanociencia, Campus Cantoblanco C/Faraday 9 28049 Madrid Spain.,Faculty of Physics, Adam Mickiewicz University 61-614 Poznań Poland
| | - Gianluca Levi
- Science Institute, University of Iceland 107 Reykjavík Iceland
| |
Collapse
|
7
|
Nottoli M, Mazzeo P, Lipparini F, Cupellini L, Mennucci B. A ΔSCF model for excited states within a polarisable embedding. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2089605] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Pisa, Italy
| | - Patrizia Mazzeo
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Pisa, Italy
- Scuola Normale Superiore, Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Pisa, Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Pisa, Italy
| |
Collapse
|
8
|
Mališ M, Vandaele E, Luber S. Spin-Orbit Couplings for Nonadiabatic Molecular Dynamics at the ΔSCF Level. J Chem Theory Comput 2022; 18:4082-4094. [PMID: 35666703 DOI: 10.1021/acs.jctc.1c01046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A procedure for the calculation of spin-orbit coupling (SOC) at the delta self-consistent field (ΔSCF) level of theory is presented. Singlet and triplet excited electronic states obtained with the ΔSCF method are expanded into a linear combination of singly excited Slater determinants composed of ground electronic state Kohn-Sham orbitals. This alleviates the nonorthogonality between excited and ground electronic states and introduces a framework, similar to the auxiliary wave function at the time-dependent density functional theory (TD-DFT) level, for the calculation of observables. The ΔSCF observables of the formaldehyde system were compared to reference TD-DFT values. Our procedure gives all components (energies, gradients, nonadiabatic couplings, and SOC terms) at the ΔSCF level of theory for conducting efficient, full-atomistic nonadiabatic molecular dynamics with intersystem crossing, particularly in condensed phase systems.
Collapse
Affiliation(s)
- Momir Mališ
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Eva Vandaele
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| |
Collapse
|
9
|
Schmerwitz YLA, Ivanov AV, Jónsson EÖ, Jónsson H, Levi G. Variational Density Functional Calculations of Excited States: Conical Intersection and Avoided Crossing in Ethylene Bond Twisting. J Phys Chem Lett 2022; 13:3990-3999. [PMID: 35481754 DOI: 10.1021/acs.jpclett.2c00741] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Theoretical studies of photochemical processes require a description of the energy surfaces of excited electronic states, especially near degeneracies, where transitions between states are most likely. Systems relevant to photochemical applications are typically too large for high-level multireference methods, and while time-dependent density functional theory (TDDFT) is efficient, it can fail to provide the required accuracy. A variational, time-independent density functional approach is applied to the twisting of the double bond and pyramidal distortion in ethylene, the quintessential model for photochemical studies. By allowing for symmetry breaking, the calculated energy surfaces exhibit the correct topology around the twisted-pyramidalized conical intersection even when using a semilocal functional approximation, and by including explicit self-interaction correction, the torsional energy curves are in close agreement with published multireference results. The findings of the present work point to the possibility of using a single determinant time-independent density functional approach to simulate nonadiabatic dynamics, even for large systems where multireference methods are impractical and TDDFT is often not accurate enough.
Collapse
Affiliation(s)
| | - Aleksei V Ivanov
- Science Institute of the University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Elvar Ö Jónsson
- Science Institute of the University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Hannes Jónsson
- Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
- Department of Applied Physics, Aalto University, FI-00076 Espoo, Finland
| | - Gianluca Levi
- Science Institute of the University of Iceland, VR-III, 107 Reykjavík, Iceland
| |
Collapse
|
10
|
Kumar C, Luber S. Robust ΔSCF calculations with direct energy functional minimization methods and STEP for molecules and materials. J Chem Phys 2022; 156:154104. [PMID: 35459303 DOI: 10.1063/5.0075927] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The direct energy functional minimization method using the orbital transformation (OT) scheme in the program package CP2K has been employed for Δ self-consistent field (ΔSCF) calculations. The OT method for non-uniform molecular orbitals occupations allows us to apply the ΔSCF method for various kinds of molecules and periodic systems. Vertical excitation energies of heteroaromatic molecules and condensed phase systems, such as solvated ethylene and solvated uracil obeying periodic boundary conditions, are reported using the ΔSCF method. In addition, a Re-phosphate molecule attached to the surface of anatase (TiO2) has been investigated. Additionally, we have implemented a recently proposed state-targeted energy projection ΔSCF algorithm [K. Carter-Fenk and J. M. Herbert, J. Chem. Theory Comput. 16(8), 5067-5082 (2020)] for diagonalization based SCF in CP2K. It is found that the OT scheme provides a smooth and robust SCF convergence for all investigated excitation energies and (non-)periodic systems.
Collapse
Affiliation(s)
- Chandan Kumar
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| |
Collapse
|