1
|
Loreti A, Freixas VM, Avagliano D, Segatta F, Song H, Tretiak S, Mukamel S, Garavelli M, Govind N, Nenov A. WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation. J Chem Theory Comput 2024; 20:4804-4819. [PMID: 38828948 DOI: 10.1021/acs.jctc.4c00310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.
Collapse
Affiliation(s)
- Alessandro Loreti
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Via Piero Gobetti 85, Bologna 40129, Italy
| | - Victor Manuel Freixas
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, California 92697, United States
| | - Davide Avagliano
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Via Piero Gobetti 85, Bologna 40129, Italy
| | - Francesco Segatta
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Via Piero Gobetti 85, Bologna 40129, Italy
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Shaul Mukamel
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, California 92697, United States
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Via Piero Gobetti 85, Bologna 40129, Italy
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Via Piero Gobetti 85, Bologna 40129, Italy
| |
Collapse
|
2
|
Shu Y, Truhlar DG. Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients. J Chem Theory Comput 2024; 20:4396-4426. [PMID: 38819014 DOI: 10.1021/acs.jctc.4c00424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
3
|
Taylor JT, Tozer DJ, Curchod BFE. On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2). J Chem Phys 2023; 159:214115. [PMID: 38059547 DOI: 10.1063/5.0176140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Accepted: 11/14/2023] [Indexed: 12/08/2023] Open
Abstract
Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order [ADC(2)]-approaches often used in excited-state molecular dynamics simulations-to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine-highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.
Collapse
Affiliation(s)
- Jack T Taylor
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - David J Tozer
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Basile F E Curchod
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
4
|
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
|
5
|
Levine BG, Durden AS, Esch MP, Liang F, Shu Y. CAS without SCF-Why to use CASCI and where to get the orbitals. J Chem Phys 2021; 154:090902. [PMID: 33685182 DOI: 10.1063/5.0042147] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree-Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Finally, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.
Collapse
Affiliation(s)
- Benjamin G Levine
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Andrew S Durden
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Michael P Esch
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Fangchun Liang
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
6
|
Levine BG, Esch MP, Fales BS, Hardwick DT, Peng WT, Shu Y. Conical Intersections at the Nanoscale: Molecular Ideas for Materials. Annu Rev Phys Chem 2019; 70:21-43. [DOI: 10.1146/annurev-physchem-042018-052425] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The ability to predict and describe nonradiative processes in molecules via the identification and characterization of conical intersections is one of the greatest recent successes of theoretical chemistry. Only recently, however, has this concept been extended to materials science, where nonradiative recombination limits the efficiencies of materials for various optoelectronic applications. In this review, we present recent advances in the theoretical study of conical intersections in semiconductor nanomaterials. After briefly introducing conical intersections, we argue that specific defects in materials can induce conical intersections between the ground and first excited electronic states, thus introducing pathways for nonradiative recombination. We present recent developments in theoretical methods, computational tools, and chemical intuition for the prediction of such defect-induced conical intersections. Through examples in various nanomaterials, we illustrate the significance of conical intersections for nanoscience. We also discuss challenges facing research in this area and opportunities for progress.
Collapse
Affiliation(s)
- Benjamin G. Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Michael P. Esch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - B. Scott Fales
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Dylan T. Hardwick
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Wei-Tao Peng
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
7
|
Levine BG, Peng WT, Esch MP. Locality of conical intersections in semiconductor nanomaterials. Phys Chem Chem Phys 2019; 21:10870-10878. [PMID: 31106323 DOI: 10.1039/c9cp01584a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A predictive theory connecting atomic structure to the rate of recombination would enable the rational design of semiconductor nanomaterials for optoelectronic applications. Recently our group has demonstrated that the theoretical study of conical intersections can serve this purpose. Here we review recent work in this area, focusing on the thesis that low-energy conical intersections in nanomaterials share a common feature: locality. We define a conical intersection as local if (a) the intersecting states differ by the excitation of an electron between spatially local orbitals, and (b) the intersection is accessed when the energies of these orbitals are tuned by local distortions of the geometry. After illustrating the locality of the conical intersection responsible for recombination at dangling bond defects in silicon, we demonstrate the locality of low-energy conical intersections in cases where locality may be a surprise. First, we demonstrate the locality of low-energy self-trapped conical intersections in a pristine silicon nanocrystal, which has no defects that one would expect to serve as the center of a local intersection. Second, we demonstrate that the lowest energy intersection in a silicon system with two neighboring dangling bond defects localizes to a single defect site. We discuss the profound implications of locality for predicting the rate of recombination and suggest that the locality of intersections could be exploited in the experimental study of recombination, where spectroscopic studies of molecular models of defects could provide new insights.
Collapse
Affiliation(s)
- Benjamin G Levine
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA.
| | - Wei-Tao Peng
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA.
| | - Michael P Esch
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA.
| |
Collapse
|
8
|
Ghosh S, Verma P, Cramer CJ, Gagliardi L, Truhlar DG. Combining Wave Function Methods with Density Functional Theory for Excited States. Chem Rev 2018; 118:7249-7292. [PMID: 30044618 DOI: 10.1021/acs.chemrev.8b00193] [Citation(s) in RCA: 128] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.
Collapse
Affiliation(s)
- Soumen Ghosh
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Christopher J Cramer
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| |
Collapse
|
9
|
Sand AM, Hoyer CE, Truhlar DG, Gagliardi L. State-interaction pair-density functional theory. J Chem Phys 2018; 149:024106. [DOI: 10.1063/1.5036727] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Affiliation(s)
- Andrew M. Sand
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute, The University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Chad E. Hoyer
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute, The University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute, The University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and The Minnesota Supercomputing Institute, The University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|