skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Model space diabatization for quantum photochemistry

Abstract

Diabatization is a procedure that transforms multiple adiabatic electronic states to a new representation in which the potential energy surfaces and the couplings between states due to the electronic Hamiltonian operator are smooth, and the couplings due to nuclear momentum are negligible. In this work, we propose a simple and general diabatization strategy, called model space diabatization, that is applicable to multi-configuration quasidegenerate perturbation theory (MC-QDPT) or its extended version (XMC-QDPT). An advantage over previous diabatization schemes is that dynamical correlation calculations are based on standard post-multi-configurational self-consistent field (MCSCF) multi-state methods even though the diabatization is based on state-averaged MCSCF results. The strategy is illustrated here by applications to LiH, LiF, and thioanisole, with the fourfold-way diabatization and XMC-QDPT, and the results illustrate its validity.

Authors:
;  [1]; ;  [2]
  1. Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455 (United States)
  2. Department of Chemistry, Iowa State University, Ames, Iowa 50011 (United States)
Publication Date:
OSTI Identifier:
22416093
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Chemical Physics; Journal Volume: 142; Journal Issue: 6; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; ANISOLE; CORRELATIONS; COUPLING; HAMILTONIANS; LITHIUM FLUORIDES; LITHIUM HYDRIDES; ORGANIC SULFUR COMPOUNDS; PERTURBATION THEORY; PHOTOCHEMISTRY; POTENTIAL ENERGY; SELF-CONSISTENT FIELD; STANDARDS; SURFACES

Citation Formats

Li, Shaohong L., Truhlar, Donald G., E-mail: truhlar@umn.edu, Schmidt, Michael W., and Gordon, Mark S.. Model space diabatization for quantum photochemistry. United States: N. p., 2015. Web. doi:10.1063/1.4907038.
Li, Shaohong L., Truhlar, Donald G., E-mail: truhlar@umn.edu, Schmidt, Michael W., & Gordon, Mark S.. Model space diabatization for quantum photochemistry. United States. doi:10.1063/1.4907038.
Li, Shaohong L., Truhlar, Donald G., E-mail: truhlar@umn.edu, Schmidt, Michael W., and Gordon, Mark S.. Sat . "Model space diabatization for quantum photochemistry". United States. doi:10.1063/1.4907038.
@article{osti_22416093,
title = {Model space diabatization for quantum photochemistry},
author = {Li, Shaohong L. and Truhlar, Donald G., E-mail: truhlar@umn.edu and Schmidt, Michael W. and Gordon, Mark S.},
abstractNote = {Diabatization is a procedure that transforms multiple adiabatic electronic states to a new representation in which the potential energy surfaces and the couplings between states due to the electronic Hamiltonian operator are smooth, and the couplings due to nuclear momentum are negligible. In this work, we propose a simple and general diabatization strategy, called model space diabatization, that is applicable to multi-configuration quasidegenerate perturbation theory (MC-QDPT) or its extended version (XMC-QDPT). An advantage over previous diabatization schemes is that dynamical correlation calculations are based on standard post-multi-configurational self-consistent field (MCSCF) multi-state methods even though the diabatization is based on state-averaged MCSCF results. The strategy is illustrated here by applications to LiH, LiF, and thioanisole, with the fourfold-way diabatization and XMC-QDPT, and the results illustrate its validity.},
doi = {10.1063/1.4907038},
journal = {Journal of Chemical Physics},
number = 6,
volume = 142,
place = {United States},
year = {Sat Feb 14 00:00:00 EST 2015},
month = {Sat Feb 14 00:00:00 EST 2015}
}
  • Cited by 17
  • We present a comprehensive photochemistry model for exploration of the chemical composition of terrestrial exoplanet atmospheres. The photochemistry model is designed from the ground up to have the capacity to treat all types of terrestrial planet atmospheres, ranging from oxidizing through reducing, which makes the code suitable for applications for the wide range of anticipated terrestrial exoplanet compositions. The one-dimensional chemical transport model treats up to 800 chemical reactions, photochemical processes, dry and wet deposition, surface emission, and thermal escape of O, H, C, N, and S bearing species, as well as formation and deposition of elemental sulfur and sulfuricmore » acid aerosols. We validate the model by computing the atmospheric composition of current Earth and Mars and find agreement with observations of major trace gases in Earth's and Mars' atmospheres. We simulate several plausible atmospheric scenarios of terrestrial exoplanets and choose three benchmark cases for atmospheres from reducing to oxidizing. The most interesting finding is that atomic hydrogen is always a more abundant reactive radical than the hydroxyl radical in anoxic atmospheres. Whether atomic hydrogen is the most important removal path for a molecule of interest also depends on the relevant reaction rates. We also find that volcanic carbon compounds (i.e., CH{sub 4} and CO{sub 2}) are chemically long-lived and tend to be well mixed in both reducing and oxidizing atmospheres, and their dry deposition velocities to the surface control the atmospheric oxidation states. Furthermore, we revisit whether photochemically produced oxygen can cause false positives for detecting oxygenic photosynthesis, and find that in 1 bar CO{sub 2}-rich atmospheres oxygen and ozone may build up to levels that have conventionally been accepted as signatures of life, if there is no surface emission of reducing gases. The atmospheric scenarios presented in this paper can serve as the benchmark atmospheres for quickly assessing the lifetime of trace gases in reducing, weakly oxidizing, and highly oxidizing atmospheres on terrestrial exoplanets for the exploration of possible biosignature gases.« less
  • We model the triplet-triplet energy-transfer experiments from the Closs group [Closs, G. L.; et al. J. Am. Chem. Soc. 1988, 110, 2652.] using a combination of Marcus theory and either Boys or Edmiston-Ruedenberg localized diabatization, and we show that relative and absolute rates of electronic excitation transfer may be computed successfully. For the case where both the donor and acceptor occupy equatorial positions on a rigid cyclohexane bridge, we find β calc = 2.8 per C-C bond, compared with the experimental value β exp = 2.6. This work highlights the power of using localized diabatization methods as a tool formore » modeling nonequilibrium processes.« less
  • We model the triplet-triplet energy-transfer experiments from the Closs group [Closs, G. L.; et al. J. Am. Chem. Soc. 1988, 110, 2652.] using a combination of Marcus theory and either Boys or Edmiston-Ruedenberg localized diabatization, and we show that relative and absolute rates of electronic excitation transfer may be computed successfully. For the case where both the donor and acceptor occupy equatorial positions on a rigid cyclohexane bridge, we find βcalc = 2.8 per C-C bond, compared with the experimental value βexp = 2.6. This work highlights the power of using localized diabatization methods as a tool for modeling nonequilibriummore » processes.« less
  • For a system which undergoes electron or energy transfer in a polar solvent, we define the diabatic states to be the initial and final states of the system, before and after the nonequilibrium transfer process. We consider two models for the system-solvent interactions: A solvent which is linearly polarized in space and a solvent which responds linearly to the system. From these models, we derive two new schemes for obtaining diabatic states from ab initio calculations of the isolated system in the absence of solvent. These algorithms resemble standard approaches for orbital localization, namely, the Boys and Edmiston-Ruedenberg (ER) formalisms.more » We show that Boys localization is appropriate for describing electron transfer [Subotnik et al., J. Chem. Phys. 129, 244101 (2008)] while ER describes both electron and energy transfer. Neither the Boys nor the ER methods require definitions of donor or acceptor fragments and both are computationally inexpensive. We investigate one chemical example, the case of oligomethylphenyl-3, and we provide attachment/detachment plots whereby the ER diabatic states are seen to have localized electron-hole pairs.« less