1
|
Lee IS, Filatov M, Min SK. Formulation of transition dipole gradients for non-adiabatic dynamics with polaritonic states. J Chem Phys 2024; 160:154103. [PMID: 38624116 DOI: 10.1063/5.0202095] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Accepted: 03/31/2024] [Indexed: 04/17/2024] Open
Abstract
A general formulation of the strong coupling between photons confined in a cavity and molecular electronic states is developed for the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham method. The light-matter interaction is included in the Jaynes-Cummings model, which requires the derivation and implementation of the analytical derivatives of the transition dipole moments between the molecular electronic states. The developed formalism is tested in the simulations of the nonadiabatic dynamics in the polaritonic states resulting from the strong coupling between the cavity photon mode and the ground and excited states of the penta-2,4-dieniminium cation, also known as PSB3. Comparison with the field-free simulations of the excited-state decay dynamics in PSB3 reveals that the light-matter coupling can considerably alter the decay dynamics by increasing the excited state lifetime and hindering photochemically induced torsion about the C=C double bonds of PSB3. The necessity of obtaining analytical transition dipole gradients for the accurate propagation of the dynamics is underlined.
Collapse
Affiliation(s)
- In Seong Lee
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
| | - Michael Filatov
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
| | - Seung Kyu Min
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, Republic of Korea
| |
Collapse
|
2
|
Wodyński A, Kaupp M. Local Hybrid Functional Applicable to Weakly and Strongly Correlated Systems. J Chem Theory Comput 2022; 18:6111-6123. [PMID: 36170626 DOI: 10.1021/acs.jctc.2c00795] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The recent idea (Wodyński, A.; Arbuznikov, A. V.; Kaupp M. J. Chem. Phys. 2021, 155, 144101) to augment local hybrid functionals by a strong-correlation (sc) factor obtained from the adiabatic connection in the spirit of the KP16 model has been extended and applied to generate the accurate sc-corrected local hybrid functional scLH22t. By damping small values of the ratio between nondynamical and dynamical correlation entering the correction factor, it has become possible to avoid double counting of nondynamical correlation for weakly correlated situations and thereby preserve the excellent accuracy of the underlying LH20t local hybrid for such cases almost perfectly. On the other hand, scLH22t improves substantially over LH20t in reducing fractional-spin errors (FSEs), in providing improved spin-restricted bond dissociation curves, and in treating some typical systems with multireference character. The obtained FSEs are similar to those of the KP16/B13 model and slightly larger than for B13, but performance for weakly correlated systems is better than for these two related methods, which are also difficult to use self-consistently. The recent DM21 functional based on the training of a deep neural network still performs somewhat better than scLH22t but allows no physical insights into the origins of reduced FSEs. Examination of local mixing functions (LMFs) for the corrected scLH22t and uncorrected LH20t functionals provides further insights: in weakly correlated situations, the LMF remains essentially unchanged. Strong-correlation effects manifest in a reduction of the LMF values in certain regions of space, even to the extent of producing negative LMF values. It is suggested that this is the mechanism by which also DM21, which may be viewed as a range-separated local hybrid, is able to reduce FSEs.
Collapse
Affiliation(s)
- Artur Wodyński
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
3
|
Baerends EJ. Chemical potential, derivative discontinuity, fractional electrons, jump of the Kohn-Sham potential, atoms as thermodynamic open systems, and other (mis)conceptions of the density functional theory of electrons in molecules. Phys Chem Chem Phys 2022; 24:12745-12766. [PMID: 35593143 DOI: 10.1039/d2cp01585d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Many references exist in the density functional theory (DFT) literature to the chemical potential of the electrons in an atom or a molecule. The origin of this notion has been the identification of the Lagrange multiplier μ = ∂E/∂N in the Euler-Lagrange variational equation for the ground state density as the chemical potential of the electrons. We first discuss why the Lagrange multiplier in this case is an arbitrary constant and therefore cannot be a physical characteristic of an atom or molecule. The switching of the energy derivative ("chemical potential") from -I to -A when the electron number crosses the integer, called integer discontinuity or derivative discontinuity, is not physical but only occurs when the nonphysical noninteger electron systems and the corresponding energy and derivative ∂E/∂N are chosen in a specific discontinuous way. The question is discussed whether in fact the thermodynamical concept of a chemical potential can be defined for the electrons in such few-electron systems as atoms and molecules. The conclusion is that such systems lack important characteristics of thermodynamic systems and do not afford the definition of a chemical potential. They also cannot be considered as analogues of the open systems of thermodynamics that can exchange particles with an environment (a particle bath or other members of a Gibbsian ensemble). Thermodynamical (statistical mechanical) concepts like chemical potential, open systems, grand canonical ensemble etc. are not applicable to a few electron system like an atom or molecule. A number of topics in DFT are critically reviewed in light of these findings: jumps in the Kohn-Sham potential when crossing an integer number of electrons, the band gap problem, the deviation-from-straight-lines error, and the role of ensembles in DFT.
Collapse
|
4
|
Elayan IA, Gupta R, Hollett JW. ΔNO and the complexities of electron correlation in simple hydrogen clusters. J Chem Phys 2022; 156:094102. [DOI: 10.1063/5.0073227] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
5
|
Lee S, Park W, Nakata H, Filatov M, Choi CH. Recent advances in ensemble density functional theory and linear response theory for strong correlation. B KOREAN CHEM SOC 2021. [DOI: 10.1002/bkcs.12429] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Seunghoon Lee
- Division of Chemistry and Chemical Engineering California Institute of Technology Pasadena California USA
| | - Woojin Park
- Department of Chemistry Kyungpook National University Daegu South Korea
| | - Hiroya Nakata
- Department of Chemistry Kyungpook National University Daegu South Korea
| | - Michael Filatov
- Department of Chemistry Kyungpook National University Daegu South Korea
| | - Cheol Ho Choi
- Department of Chemistry Kyungpook National University Daegu South Korea
| |
Collapse
|
6
|
Filatov M, Lee S, Choi CH. Description of Sudden Polarization in the Excited Electronic States with an Ensemble Density Functional Theory Method. J Chem Theory Comput 2021; 17:5123-5139. [PMID: 34319730 DOI: 10.1021/acs.jctc.1c00479] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Sudden polarization (SP) is one of the manifestations of electron transfer in the electronically excited states of molecules. Proposed initially to explain the unusual reactivity of photoexcited olefins, SP often occurs in the excited states of molecules possessing strongly correlated diradical ground state. Theoretical description of SP involves mixing between the singly excited and the doubly excited zwitterionic states, which makes it inaccessible with the use of the popular linear-response time-dependent density functional theory methods. In this work, an extended variant of the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS, or SSR) method is applied to study SP in a number of organic diradical systems. To this end, the analytical derivative formalism is derived and implemented for the SSR(3,2) method (see the main text for explanation of the acronym), which enables the automatic geometry optimization and obtains the relaxed density matrices as well as the electron binding energies and respective Dyson's orbitals. Application of the new method to SP in the lowest singlet excited state of ethylene agrees with the results obtained previously with the use of multireference methods of wavefunction theory. A number of interesting manifestations of SP are observed, such as the charge transfer in photoexcited tetramethyleneethene (TME) diradical mediated by the vibrational motion and conductivity switching in the excited state of a donor-acceptor dyad placed in an external electric field.
Collapse
Affiliation(s)
- Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| |
Collapse
|
7
|
Desmarais JK, Erba A, Flament JP, Kirtman B. Perturbation Theory Treatment of Spin-Orbit Coupling II: A Coupled Perturbed Kohn-Sham Method. J Chem Theory Comput 2021; 17:4712-4732. [PMID: 34286577 DOI: 10.1021/acs.jctc.1c00460] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A noncanonical coupled perturbed Kohn-Sham density functional theory (KS-DFT)/Hartree-Fock (HF) treatment of spin-orbit coupling (SOC) is provided. We take the scalar-relativistic KS-DFT/HF solution, obtained with a relativistic effective core potential, as the zeroth-order approximation. Explicit expressions are given for the total energy through the 4th order, which satisfy the 2n + 1 rule. Second-order expressions are provided for orbital energies and density variables of spin-current DFT. Test calculations are carried out on the halogen homonuclear diatomic and hydride molecules, including 6p and 7p elements, as well as open-shell negative ions. The computed properties through second or third order match well with those from reference two-component self-consistent field calculations for total and orbital energies as well as spin-current densities. In only one case (At2-) did a significant deviation occur for the remaining density variables. Our coupled perturbation theory approach provides an efficient way of adding the effect of SOC to a scalar-relativistic single-reference KS-DFT/HF treatment, in particular because it does not require diagonalization in the two-component spinor basis, leading to saving factors on the number of required floating-point operations that may exceed one order of magnitude.
Collapse
Affiliation(s)
- Jacques K Desmarais
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy.,Université de Pau et des Pays de l'Adour, E2S UPPA, CNRS, IPREM, 64000 Pau, France
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Bernard Kirtman
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
8
|
Desmarais JK, Komorovsky S, Flament JP, Erba A. Spin–orbit coupling from a two-component self-consistent approach. II. Non-collinear density functional theories. J Chem Phys 2021; 154:204110. [DOI: 10.1063/5.0051447] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Jacques K. Desmarais
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
- Université de Pau et des Pays de l’Adour, E2S UPPA, CNRS, IPREM, Pau, France
| | - Stanislav Komorovsky
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Dúbravská cesta 9, SK-84536 Bratislava, Slovakia
| | - Jean-Pierre Flament
- Université de Lille, CNRS, UMR 8523—PhLAM—Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| |
Collapse
|
9
|
Filatov M, Lee S, Nakata H, Choi CH. Signatures of Conical Intersection Dynamics in the Time-Resolved Photoelectron Spectrum of Furan: Theoretical Modeling with an Ensemble Density Functional Theory Method. Int J Mol Sci 2021; 22:4276. [PMID: 33924097 PMCID: PMC8074317 DOI: 10.3390/ijms22084276] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Revised: 04/13/2021] [Accepted: 04/15/2021] [Indexed: 12/13/2022] Open
Abstract
The non-adiabatic dynamics of furan excited in the ππ* state (S2 in the Franck-Condon geometry) was studied using non-adiabatic molecular dynamics simulations in connection with an ensemble density functional method. The time-resolved photoelectron spectra were theoretically simulated in a wide range of electron binding energies that covered the valence as well as the core electrons. The dynamics of the decay (rise) of the photoelectron signal were compared with the excited-state population dynamics. It was observed that the photoelectron signal decay parameters at certain electron binding energies displayed a good correlation with the events occurring during the excited-state dynamics. Thus, the time profile of the photoelectron intensity of the K-shell electrons of oxygen (decay constant of 34 ± 3 fs) showed a reasonable correlation with the time of passage through conical intersections with the ground state (47 ± 2 fs). The ground-state recovery constant of the photoelectron signal (121 ± 30 fs) was in good agreement with the theoretically obtained excited-state lifetime (93 ± 9 fs), as well as with the experimentally estimated recovery time constant (ca. 110 fs). Hence, it is proposed to complement the traditional TRPES observations with the trXPS (or trNEXAFS) measurements to obtain more reliable estimates of the most mechanistically important events during the excited-state dynamics.
Collapse
Affiliation(s)
- Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 702-701, Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA 91125, USA;
| | - Hiroya Nakata
- R & D Center Kagoshima, Kyocera, 1-4 Kokubu Yamashita-cho, Kirishima-shi, Kagoshima 899-4312, Japan;
| | - Cheol-Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 702-701, Korea
| |
Collapse
|
10
|
Filatov M, Lee S, Nakata H, Choi CH. Computation of Molecular Electron Affinities Using an Ensemble Density Functional Theory Method. J Phys Chem A 2020; 124:7795-7804. [PMID: 32900199 DOI: 10.1021/acs.jpca.0c06976] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
The computation of electron attachment energies (electron affinities) was implemented in connection with an ensemble density functional theory method, the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method. With the use of the extended Koopmans' theorem, the electron affinities and the respective Dyson orbitals are obtained directly for the neutral molecule, thus avoiding the necessity to compute the ionized system. Together with the EKT-SSR (extended Koopmans' theorem-SSR) method for ionization potentials, which was developed earlier, EKT-SSR for electron affinities completes the implementation of the EKT-SSR formalism, which can now be used for obtaining electron detachment as well as the electron attachment energies of molecules in the ground and excited electronic states. The extended EKT-SSR method was tested in the calculation of several closed-shell molecules. For the molecules in the ground states, the EKT-SSR energies of Dyson's orbitals are virtually identical to the energies of the unoccupied orbitals in the usual single-reference spin-restricted Kohn-Sham calculations. For the molecules in the excited states, EKT-SSR predicts an increase of the most positive electron affinity by approximately the amount of the vertical excitation energy. The electron affinities of a number of diradicals were calculated with EKT-SSR and compared with the available experimental data. With the use of a standard density functional (BH&HLYP), the EKT-SSR electron affinities deviate on average by ca. 0.2 eV from the experimental data. It is expected that the agreement with the experiment can be improved by designing density functionals parametrized for ionization energies.
Collapse
Affiliation(s)
- Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Hiroya Nakata
- R&D Center Kagoshima, Kyocera, 1-4 Kokubu Yamashita-cho, Kirishima-shi, Kagoshima 899-4312, Japan
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| |
Collapse
|
11
|
Filatov M, Lee S, Choi CH. Computation of Molecular Ionization Energies Using an Ensemble Density Functional Theory Method. J Chem Theory Comput 2020; 16:4489-4504. [PMID: 32421323 DOI: 10.1021/acs.jctc.0c00218] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Computation of the ionization energies and of the respective Dyson orbitals based on the use of the extended Koopmans theorem (EKT) is implemented in connection with an ensemble density functional theory (eDFT) method, the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method. The new methodology enables fast computation of the ionization energies and evaluation of the respective Dyson orbitals, the square norms of which are related with the ionization probabilities, in the ground and excited electronic states of molecules. As the application of EKT recycles the intermediate quantities from the SSR analytical energy gradient, evaluation of the ionization energies and probabilities can be carried out on-the-fly during the nonadiabatic molecular dynamics simulations. This opens up a perspective for fast theoretical simulation of the time-resolved photoelectron spectroscopy observations. In the present work, the new methodology is tested in the computation of the ionization energies and Dyson orbitals of several molecules in the ground and excited electronic states, including strongly correlated species, such as the ozone molecule, dissociating chemical bonds, and conical intersections.
Collapse
Affiliation(s)
- Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| |
Collapse
|
12
|
Takayanagi T. Quantum dynamics analysis of transition-state spectrum for the SH + H 2S → H 2S + SH reaction. Phys Chem Chem Phys 2020; 22:19845-19854. [DOI: 10.1039/d0cp03072d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Quantum dynamics calculations were performed to understand transition-state spectroscopy of the SH + H2S hydrogen atom transfer reaction.
Collapse
|
13
|
Filatov M, Liu F, Martínez TJ. Analytical derivatives of the individual state energies in ensemble density functional theory method. I. General formalism. J Chem Phys 2018; 147:034113. [PMID: 28734302 DOI: 10.1063/1.4994542] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The state-averaged (SA) spin restricted ensemble referenced Kohn-Sham (REKS) method and its state interaction (SI) extension, SI-SA-REKS, enable one to describe correctly the shape of the ground and excited potential energy surfaces of molecules undergoing bond breaking/bond formation reactions including features such as conical intersections crucial for theoretical modeling of non-adiabatic reactions. Until recently, application of the SA-REKS and SI-SA-REKS methods to modeling the dynamics of such reactions was obstructed due to the lack of the analytical energy derivatives. In this work, the analytical derivatives of the individual SA-REKS and SI-SA-REKS energies are derived. The final analytic gradient expressions are formulated entirely in terms of traces of matrix products and are presented in the form convenient for implementation in the traditional quantum chemical codes employing basis set expansions of the molecular orbitals. The implementation and benchmarking of the derived formalism will be described in a subsequent article of this series.
Collapse
Affiliation(s)
- Michael Filatov
- Department of Chemistry, School of Natural Sciences, Ulsan National Institute of Science and Technology (UNIST), Ulsan 689-798, South Korea
| | - Fang Liu
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| |
Collapse
|
14
|
Filatov M, Martínez TJ, Kim KS. Description of ground and excited electronic states by ensemble density functional method with extended active space. J Chem Phys 2018; 147:064104. [PMID: 28810777 DOI: 10.1063/1.4996873] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An extended variant of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, the REKS(4,4) method, designed to describe the ground electronic states of strongly multireference systems is modified to enable calculation of excited states within the time-independent variational formalism. The new method, the state-interaction state-averaged REKS(4,4), i.e., SI-SA-REKS(4,4), is capable of describing several excited states of a molecule involving double bond cleavage, polyradical character, or multiple chromophoric units. We demonstrate that the new method correctly describes the ground and the lowest singlet excited states of a molecule (ethylene) undergoing double bond cleavage. The applicability of the new method for excitonic states is illustrated with π stacked ethylene and tetracene dimers. We conclude that the new method can describe a wide range of multireference phenomena.
Collapse
Affiliation(s)
- Michael Filatov
- Department of Chemistry, School of Natural Sciences, Ulsan National Institute of Science and Technology (UNIST), Ulsan 44919, South Korea
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Kwang S Kim
- Department of Chemistry, School of Natural Sciences, Ulsan National Institute of Science and Technology (UNIST), Ulsan 44919, South Korea
| |
Collapse
|
15
|
van Leeuwen R, Neugebauer J, Visscher L, Bickelhaupt FM. Editorial for PCCP themed issue "Developments in Density Functional Theory". Phys Chem Chem Phys 2018; 18:20864-7. [PMID: 27406292 DOI: 10.1039/c6cp90143c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This issue provides an overview of the state-of-the-art of DFT, ranging from mathematical and software developments, via topics in chemical bonding theory, to all kinds of molecular and material properties. Through this issue, we also celebrate the enormous contributions that Evert Jan Baerends has made to this field.
Collapse
Affiliation(s)
- Robert van Leeuwen
- Department of Physics, University of Jyväskylä, Survontie 9 40014 Jyväskylä, Finland.
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation (CMTC), Westfälische Wilhelms-Universität Münster, Corrensstrasse 40, 48149 Münster, Germany.
| | - Lucas Visscher
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands.
| | - F Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands. and Institute of Molecules and Materials (IMM), Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands.
| |
Collapse
|
16
|
Mahler A, Janesko BG, Moncho S, Brothers EN. Why are GGAs so accurate for reaction kinetics on surfaces? Systematic comparison of hybrid vs. nonhybrid DFT for representative reactions. J Chem Phys 2017. [DOI: 10.1063/1.4986404] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
17
|
Gao J, Grofe A, Ren H, Bao P. Beyond Kohn-Sham Approximation: Hybrid Multistate Wave Function and Density Functional Theory. J Phys Chem Lett 2016; 7:5143-5149. [PMID: 27973892 PMCID: PMC5793875 DOI: 10.1021/acs.jpclett.6b02455] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
A multistate density functional theory (MSDFT) is presented in which the energies and densities for the ground and excited states are treated on the same footing using multiconfigurational approaches. The method can be applied to systems with strong correlation and to correctly describe the dimensionality of the conical intersections between strongly coupled dissociative potential energy surfaces. A dynamic-then-static framework for treating electron correlation is developed to first incorporate dynamic correlation into contracted state functions through block-localized Kohn-Sham density functional theory (KSDFT), followed by diagonalization of the effective Hamiltonian to include static correlation. MSDFT can be regarded as a hybrid of wave function and density functional theory. The method is built on and makes use of the current approximate density functional developed in KSDFT, yet it retains its computational efficiency to treat strongly correlated systems that are problematic for KSDFT but too large for accurate WFT. The results presented in this work show that MSDFT can be applied to photochemical processes involving conical intersections.
Collapse
Affiliation(s)
- Jiali Gao
- Theoretical Chemistry Institute, Jilin University, Changchun 130023, P. R. China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Adam Grofe
- Theoretical Chemistry Institute, Jilin University, Changchun 130023, P. R. China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Haisheng Ren
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Peng Bao
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
| |
Collapse
|
18
|
Filatov M, Martínez TJ, Kim KS. Using the GVB Ansatz to develop ensemble DFT method for describing multiple strongly correlated electron pairs. Phys Chem Chem Phys 2016; 18:21040-50. [PMID: 26947515 DOI: 10.1039/c6cp00236f] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Ensemble density functional theory (DFT) furnishes a rigorous theoretical framework for describing the non-dynamic electron correlation arising from (near) degeneracy of several electronic configurations. Ensemble DFT naturally leads to fractional occupation numbers (FONs) for several Kohn-Sham (KS) orbitals, which thereby become variational parameters of the methodology. The currently available implementation of ensemble DFT in the form of the spin-restricted ensemble-referenced KS (REKS) method was originally designed for systems with only two fractionally occupied KS orbitals, which was sufficient to accurately describe dissociation of a single chemical bond or the singlet ground state of biradicaloid species. To extend applicability of the method to systems with several dissociating bonds or to polyradical species, more fractionally occupied orbitals must be included in the ensemble description. Here we investigate a possibility of developing the extended REKS methodology with the help of the generalized valence bond (GVB) wavefunction theory. The use of GVB enables one to derive a simple and physically transparent energy expression depending explicitly on the FONs of several KS orbitals. In this way, a version of the REKS method with four electrons in four fractionally occupied orbitals is derived and its accuracy in the calculation of various types of strongly correlated molecules is investigated. We propose a possible scheme to ameliorate the partial size-inconsistency that results from perfect spin-pairing. We conjecture that perfect pairing natural orbital (NO) functionals of reduced density matrix functional theory (RDMFT) should also display partial size-inconsistency.
Collapse
Affiliation(s)
- Michael Filatov
- Department of Chemistry, School of Natural Sciences, Ulsan National Institute of Science and Technology (UNIST), Ulsan 689-798, Korea.
| | | | | |
Collapse
|
19
|
Nikiforov A, Gamez JA, Thiel W, Huix-Rotllant M, Filatov M. Assessment of approximate computational methods for conical intersections and branching plane vectors in organic molecules. J Chem Phys 2015; 141:124122. [PMID: 25273427 DOI: 10.1063/1.4896372] [Citation(s) in RCA: 59] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 Å or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems.
Collapse
Affiliation(s)
- Alexander Nikiforov
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Jose A Gamez
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Miquel Huix-Rotllant
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, D-60438 Frankfurt am Main, Germany
| | - Michael Filatov
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
20
|
Filatov M, Huix-Rotllant M, Burghardt I. Ensemble density functional theory method correctly describes bond dissociation, excited state electron transfer, and double excitations. J Chem Phys 2015; 142:184104. [DOI: 10.1063/1.4919773] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Michael Filatov
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, USA
| | - Miquel Huix-Rotllant
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, D-60438 Frankfurt am Main, Germany
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, D-60438 Frankfurt am Main, Germany
| |
Collapse
|
21
|
Becke AD. Perspective: Fifty years of density-functional theory in chemical physics. J Chem Phys 2015; 140:18A301. [PMID: 24832308 DOI: 10.1063/1.4869598] [Citation(s) in RCA: 648] [Impact Index Per Article: 72.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Since its formal inception in 1964-1965, Kohn-Sham density-functional theory (KS-DFT) has become the most popular electronic structure method in computational physics and chemistry. Its popularity stems from its beautifully simple conceptual framework and computational elegance. The rise of KS-DFT in chemical physics began in earnest in the mid 1980s, when crucial developments in its exchange-correlation term gave the theory predictive power competitive with well-developed wave-function methods. Today KS-DFT finds itself under increasing pressure to deliver higher and higher accuracy and to adapt to ever more challenging problems. If we are not mindful, however, these pressures may submerge the theory in the wave-function sea. KS-DFT might be lost. I am hopeful the Kohn-Sham philosophical, theoretical, and computational framework can be preserved. This Perspective outlines the history, basic concepts, and present status of KS-DFT in chemical physics, and offers suggestions for its future development.
Collapse
Affiliation(s)
- Axel D Becke
- Department of Chemistry, Dalhousie University, 6274 Coburg Rd., P.O. Box 15000, Halifax, Nova Scotia B3H 4R2, Canada
| |
Collapse
|
22
|
Huix-Rotllant M, Nikiforov A, Thiel W, Filatov M. Description of Conical Intersections with Density Functional Methods. Top Curr Chem (Cham) 2015; 368:445-76. [PMID: 25896441 DOI: 10.1007/128_2015_631] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Conical intersections are perhaps the most significant mechanistic features of chemical reactions occurring through excited states. By providing funnels for efficient non-adiabatic population transfer, conical intersections govern the branching ratio of products of such reactions, similar to what the transition states do for ground-state reactivity. In this regard, intersections between the ground and the lowest excited states play a special role, and the correct description of the potential energy surfaces in their vicinity is crucial for understanding the mechanism and dynamics of excited-state reactions. The methods of density functional theory, such as time-dependent density functional theory, are widely used to describe the excited states of large molecules. However, are these methods suitable for describing the conical intersections or do they lead to artifacts and, consequently, to erroneous description of reaction dynamics? Here we address the first part of this question and analyze the ability of several density functional approaches, including the linear-response time-dependent approach as well as the spin-flip and ensemble formalisms, to provide the correct description of conical intersections and the potential energy surfaces in their vicinity. It is demonstrated that the commonly used linear-response time-dependent theory does not yield a proper description of these features and that one should instead use alternative computational approaches.
Collapse
Affiliation(s)
- Miquel Huix-Rotllant
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438, Frankfurt am Main, Germany
| | | | | | | |
Collapse
|
23
|
Filatov M. Ensemble DFT Approach to Excited States of Strongly Correlated Molecular Systems. Top Curr Chem (Cham) 2015; 368:97-124. [PMID: 25906417 DOI: 10.1007/128_2015_630] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Ensemble density functional theory (DFT) is a novel time-independent formalism for obtaining excitation energies of many-body fermionic systems. A considerable advantage of ensemble DFT over the more common Kohn-Sham (KS) DFT and time-dependent DFT formalisms is that it enables one to account for strong non-dynamic electron correlation in the ground and excited states of molecular systems in a transparent and accurate fashion. Despite its positive aspects, ensemble DFT has not so far found its way into the repertoire of methods of modern computational chemistry, probably because of the perceived lack of practically affordable implementations of the theory. The spin-restricted ensemble-referenced KS (REKS) method is perhaps the first computationally feasible implementation of the ideas behind ensemble DFT which enables one to describe accurately electronic transitions in a wide class of molecular systems, including strongly correlated molecules (biradicals, molecules undergoing bond breaking/formation), extended π-conjugated systems, donor-acceptor charge transfer adducts, etc.
Collapse
Affiliation(s)
- Michael Filatov
- Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstr. 4, 53115, Bonn, Germany,
| |
Collapse
|
24
|
Becke AD. Fractional Kohn–Sham Occupancies from a Strong-Correlation Density Functional. Top Curr Chem (Cham) 2014. [DOI: 10.1007/128_2014_581] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
|
25
|
Filatov M. Description of electron transfer in the ground and excited states of organic donor–acceptor systems by single-reference and multi-reference density functional methods. J Chem Phys 2014; 141:124123. [DOI: 10.1063/1.4896455] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michael Filatov
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
26
|
Filatov M, Huix-Rotllant M. Assessment of density functional theory based ΔSCF (self-consistent field) and linear response methods for longest wavelength excited states of extended π-conjugated molecular systems. J Chem Phys 2014; 141:024112. [DOI: 10.1063/1.4887087] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Michael Filatov
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Miquel Huix-Rotllant
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, D-60438 Frankfurt am Main, Germany
| |
Collapse
|
27
|
Filatov M, Olivucci M. Designing Conical Intersections for Light-Driven Single Molecule Rotary Motors: From Precessional to Axial Motion. J Org Chem 2014; 79:3587-600. [DOI: 10.1021/jo5004289] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Michael Filatov
- Institut
für Physikalische und Theoretische Chemie, Universität Bonn, Beringstrasse 4, D-53115 Bonn, Germany
| | - Massimo Olivucci
- Dipartimento
di Chimica, Università di Siena, Siena, Italy
- Department
of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| |
Collapse
|
28
|
Varsano D, Barborini M, Guidoni L. Kohn-Sham orbitals and potentials from quantum Monte Carlo molecular densities. J Chem Phys 2014; 140:054102. [DOI: 10.1063/1.4863213] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
29
|
Pribram-Jones A, Pittalis S, Gross EKU, Burke K. Thermal Density Functional Theory in Context. LECTURE NOTES IN COMPUTATIONAL SCIENCE AND ENGINEERING 2014. [DOI: 10.1007/978-3-319-04912-0_2] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
30
|
Filatov M. Assessment of Density Functional Methods for Obtaining Geometries at Conical Intersections in Organic Molecules. J Chem Theory Comput 2013; 9:4526-4541. [PMID: 24124402 PMCID: PMC3794441 DOI: 10.1021/ct400598b] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2013] [Indexed: 11/28/2022]
Abstract
![]()
A number of commonly
available density functionals have been tested
for their ability to describe the energetics and the geometry at conical
intersections in connection with the spin-restricted ensemble referenced
Kohn–Sham (REKS) method. The minimum energy conical intersections
have been optimized for several molecular systems, which are widely
used as paradigmatic models of photochemical rearrangements and models
of biological chromophores. The results of the calculations are analyzed
using the sign-change theorem of Longuet-Higgins and a method of elementary
reaction coordinates of Haas et al. The latter approach helps to elucidate
the differences between the geometries at conical intersections as
predicted by the multireference wave function ab initio methods and by the density functional methods. Overall, the BH&HLYP
density functional yields the best results for the conical intersection
geometries and energetics.
Collapse
Affiliation(s)
- Michael Filatov
- Institut für Physikalische und Theoretische Chemie, Universität Bonn , Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
31
|
Huix-Rotllant M, Filatov M, Gozem S, Schapiro I, Olivucci M, Ferré N. Assessment of Density Functional Theory for Describing the Correlation Effects on the Ground and Excited State Potential Energy Surfaces of a Retinal Chromophore Model. J Chem Theory Comput 2013; 9:3917-32. [PMID: 26592387 DOI: 10.1021/ct4003465] [Citation(s) in RCA: 78] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
In the quest for a cost-effective level of theory able to describe a large portion of the ground and excited potential energy surfaces of large chromophores, promising approaches are rooted in various approximations to the exact density functional theory (DFT). In the present work, we investigate how generalized Kohn-Sham DFT (GKS-DFT), time-dependent DFT (TDDFT), and spin-restricted ensemble-DFT (REKS) methods perform along three important paths characterizing a model retinal chromophore (the penta-2,4-dieniminium cation) in a region of near-degeneracy (close to a conical intersection) with respect to reference high-level multiconfigurational wave function methods. If GKS-DFT correctly describes the closed-shell charge transfer state, only TDDFT and REKS approaches give access to the open-shell diradical, one which sometimes corresponds to the electronic ground state. It is demonstrated that the main drawback of the usual DFT-based methods lies in the absence of interactions between the charge transfer and the diradicaloid configurations. Hence, we test a new computational scheme based on the State-averaged REKS (SA-REKS) approach, which explicitly includes these interactions into account. The State-Interaction SA-REKS (SI-SA-REKS) method significantly improves on the REKS and the SA-REKS results for the target system. The similarities and differences between DFT and wave function-based approaches are analyzed according to (1) the active space dimensions of the wave function-based methods and (2) the relative electronegativities of the allyl and protonated Schiff base moieties.
Collapse
Affiliation(s)
- Miquel Huix-Rotllant
- Aix-Marseille Université, CNRS, Institut de Chimie Radicalaire , 13397 Marseille Cedex 20, France
| | - Michael Filatov
- Institut für Physicalische und Theoretische Chemie, Universität Bonn , Beringstr. 4, 53115 Bonn, Germany
| | - Samer Gozem
- Department of Chemistry, Bowling Green State University , Bowling Green, Ohio 43402, United States
| | - Igor Schapiro
- Department of Chemistry, Bowling Green State University , Bowling Green, Ohio 43402, United States.,Max Planck Institute for Chemical Energy Conversion , Stiftstr. 34 - 36, Mülheim an der Ruhr, Germany
| | - Massimo Olivucci
- Department of Chemistry, Bowling Green State University , Bowling Green, Ohio 43402, United States.,Dipartimento di Chimica, Università di Siena , Siena, Italy
| | - Nicolas Ferré
- Aix-Marseille Université, CNRS, Institut de Chimie Radicalaire , 13397 Marseille Cedex 20, France
| |
Collapse
|
32
|
Becke AD. Communication: Two-determinant mixing with a strong-correlation density functional. J Chem Phys 2013; 139:021104. [DOI: 10.1063/1.4812926] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
33
|
Klüpfel S, Klüpfel P, Jónsson H. The effect of the Perdew-Zunger self-interaction correction to density functionals on the energetics of small molecules. J Chem Phys 2013; 137:124102. [PMID: 23020319 DOI: 10.1063/1.4752229] [Citation(s) in RCA: 69] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Self-consistent calculations using the Perdew-Zunger self-interaction correction (PZ-SIC) to local density and gradient dependent energy functionals are presented for the binding energy and equilibrium geometry of small molecules as well as energy barriers of reactions. The effect of the correction is to reduce binding energy and bond lengths and increase activation energy barriers when bond breaking is involved. The accuracy of the corrected functionals varies strongly, the correction to the binding energy being too weak for the local density approximation but too strong for the gradient dependent functionals considered. For the Perdew, Burke, and Ernzerhof (PBE) functional, a scaling of the PZ-SIC by one half gives improved results on average for both binding energy and bond lengths. The PZ-SIC does not necessarily give more accurate total energy, but it can result in a better cancellation of errors. An essential aspect of these calculations is the use of complex orbitals. A restriction to real orbitals leads to less accurate results as was recently shown for atoms [S. Klüpfel, P. Klüpfel, and H. Jónsson, Phys. Rev. A 84, 050501 (2011)]. The molecular geometry of radicals can be strongly affected by PZ-SIC. An incorrect, non-linear structure of the C(2)H radical predicted by PBE is corrected by PZ-SIC. The CH(3) radical is correctly predicted to be planar when complex orbitals are used, while it is non-planar when the PZ-SIC calculation is restricted to real orbitals.
Collapse
Affiliation(s)
- Simon Klüpfel
- Science Institute of the University of Iceland, Reykjavík, Iceland
| | | | | |
Collapse
|
34
|
Ess DH, Cook TC. Unrestricted prescriptions for open-shell singlet diradicals: using economical ab initio and density functional theory to calculate singlet-triplet gaps and bond dissociation curves. J Phys Chem A 2012; 116:4922-9. [PMID: 22578025 DOI: 10.1021/jp300633j] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Here we present and test several computational prescriptions for calculating singlet-triplet (ST) gap energies and bond dissociation curves for open-shell singlet diradicals using economical unrestricted single reference type calculations. For ST gap energies from Slipchenko and Krylov's atom and molecule test set (C, O, Si, NH, NF, OH(+), O(2), CH(2), and NH(2)(+)) spin unrestricted Hartree-Fock and MP2 energies result in errors greater than 15 kcal/mol. However, spin-projected (SP) Hartree-Fock theory in combination with spin-component-scaled (SCS) or scaled-opposite-spin (SOS) second-order perturbation theory gives ST gap energies with a mean unsigned error (MUE) of less than 2 kcal/mol. Density functionals generally give poor results for unrestricted energies and only the ωB97X-D, the M06, and the M06-2X functionals provide reasonable accuracy after spin-projection with MUE values of 4.7, 4.3, and 3.0 kcal/mol, respectively, with the 6-311++G(2d,2p) basis set. We also present a new one parameter hybrid density functional, diradical-1 (DR-1), based on Adamo and Barone's modified PW exchange functional with the PW91 correlation functional. This DR-1 method gives a mean error (ME) of 0.0 kcal/mol and a MUE value of 1.3 kcal/mol for ST gap energies. As another test of unrestricted methods the bond dissociation curves for methane (CH(4)) and hydrofluoric acid (H-F) were calculated with the M06-2X, DR-1, and ωB97X-D density functionals. All three of these functionals give reasonable results for the methane C-H bond but result in errors greater than 50 kcal/mol for the H-F bond dissociation. Spin-projection is found to significantly degrade bond dissociation curves past ~2.2 Å. Although unrestricted Hartree-Fock theory provides a very poor description of H-F bond dissociation, unrestricted SCS-MP2 and SOS-MP2 methods give accurate results.
Collapse
Affiliation(s)
- Daniel H Ess
- Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602, USA.
| | | |
Collapse
|
35
|
Abstract
The classic density-functional theory (DFT) formalism introduced by Hohenberg, Kohn, and Sham in the mid-1960s is based on the idea that the complicated N-electron wave function can be replaced with the mathematically simpler 1-electron charge density in electronic structure calculations of the ground stationary state. As such, ordinary DFT cannot treat time-dependent (TD) problems nor describe excited electronic states. In 1984, Runge and Gross proved a theorem making TD-DFT formally exact. Information about electronic excited states may be obtained from this theory through the linear response (LR) theory formalism. Beginning in the mid-1990s, LR-TD-DFT became increasingly popular for calculating absorption and other spectra of medium- and large-sized molecules. Its ease of use and relatively good accuracy has now brought LR-TD-DFT to the forefront for this type of application. As the number and the diversity of applications of TD-DFT have grown, so too has our understanding of the strengths and weaknesses of the approximate functionals commonly used for TD-DFT. The objective of this article is to continue where a previous review of TD-DFT in Volume 55 of the Annual Review of Physical Chemistry left off and highlight some of the problems and solutions from the point of view of applied physical chemistry. Because doubly-excited states have a particularly important role to play in bond dissociation and formation in both thermal and photochemistry, particular emphasis is placed on the problem of going beyond or around the TD-DFT adiabatic approximation, which limits TD-DFT calculations to nominally singly-excited states.
Collapse
Affiliation(s)
- M E Casida
- Laboratoire de Chimie Théorique, Département de Chimie Moléculaire, Institut de Chimie Moléculaire de Grenoble, Université Joseph Fourier, France.
| | | |
Collapse
|
36
|
Ess DH, Johnson ER, Hu X, Yang W. Singlet−Triplet Energy Gaps for Diradicals from Fractional-Spin Density-Functional Theory. J Phys Chem A 2010; 115:76-83. [DOI: 10.1021/jp109280y] [Citation(s) in RCA: 98] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- Daniel H. Ess
- Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602, United States, School of Natural Sciences, University of California—Merced, 5200 North Lake Road, Merced, California 95343, United States, and Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Erin R. Johnson
- Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602, United States, School of Natural Sciences, University of California—Merced, 5200 North Lake Road, Merced, California 95343, United States, and Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Xiangqian Hu
- Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602, United States, School of Natural Sciences, University of California—Merced, 5200 North Lake Road, Merced, California 95343, United States, and Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Weitao Yang
- Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602, United States, School of Natural Sciences, University of California—Merced, 5200 North Lake Road, Merced, California 95343, United States, and Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
37
|
Hermann A, Schwerdtfeger P. Complete basis set limit second-order Møller–Plesset calculations for the fcc lattices of neon, argon, krypton, and xenon. J Chem Phys 2009; 131:244508. [DOI: 10.1063/1.3279303] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
38
|
|
39
|
Janesko BG, Scuseria GE. Hartree-Fock orbitals significantly improve the reaction barrier heights predicted by semilocal density functionals. J Chem Phys 2008; 128:244112. [PMID: 18601322 DOI: 10.1063/1.2940738] [Citation(s) in RCA: 75] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Semilocal density functional theory predictions for the barrier heights of representative hydrogen transfer, heavy-atom transfer, and nucleophilic substitution reactions are significantly improved in non-self-consistent calculations using Hartree-Fock orbitals. Orbitals from hybrid calculations yield related improvements. These results provide insight into compensating for one-electron self-interaction error in semilocal density functional theory.
Collapse
|
40
|
Bahmann H, Ernzerhof M. Generalized-gradient exchange-correlation hole obtained from a correlation factor ansatz. J Chem Phys 2008; 128:234104. [DOI: 10.1063/1.2937447] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Wu Q, Cheng CL, Van Voorhis T. Configuration interaction based on constrained density functional theory: a multireference method. J Chem Phys 2008; 127:164119. [PMID: 17979331 DOI: 10.1063/1.2800022] [Citation(s) in RCA: 122] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Existing density functional theory (DFT) methods are typically very effective in capturing dynamic correlation, but run into difficulty treating near-degenerate systems where static correlation becomes important. In this work, we propose a configuration interaction (CI) method that allows one to use a multireference approach to treat static correlation but incorporates DFT's efficacy for the dynamic part as well. The new technique uses localized charge or spin states built by a constrained DFT approach to construct an active space in which the effective Hamiltonian matrix is built. These local configurations have significantly less static correlation compared to their delocalized counterparts and possess an essentially constant amount of self-interaction error. Thus their energies can be reliably calculated by DFT with existing functionals. Using a small number of local configurations as different references in the active space, a simple CI step is then able to recover the static correlation missing from the localized states. Practical issues of choosing configurations and adjusting constraint values are discussed, employing as examples the ground state dissociation curves of H(2) (+), H(2), and LiF. Excellent results are obtained for these curves at all interatomic distances, which is a strong indication that this method can be used to accurately describe bond breaking and forming processes.
Collapse
Affiliation(s)
- Qin Wu
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA.
| | | | | |
Collapse
|
42
|
Sheehan SM, Parsons BF, Zhou J, Garand E, Yen TA, Moore DT, Neumark DM. Characterization of cyclic and linear C3H− and C3H via anion photoelectron spectroscopy. J Chem Phys 2008; 128:034301. [DOI: 10.1063/1.2812561] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
|
43
|
Bochevarov AD, Friesner RA. The densities produced by the density functional theory: Comparison to full configuration interaction. J Chem Phys 2008; 128:034102. [DOI: 10.1063/1.2821123] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
|
44
|
Lawson Daku LM, Linares J, Boillot ML. Ab Initio Static and Molecular Dynamics Study of 4-Styrylpyridine. Chemphyschem 2007; 8:1402-16. [PMID: 17530730 DOI: 10.1002/cphc.200700117] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
We report an in-depth theoretical study of 4-styrylpyridine in its singlet S(0) ground state. The geometries and the relative stabilities of the trans and cis isomers were investigated within density functional theory (DFT) as well as within Hartree-Fock (HF), second-order Møller-Plesset (MP2), and coupled cluster (CC) theories. The DFT calculations were performed using the B3LYP and PBE functionals, with basis sets of different qualities, and gave results that are very consistent with each other. The molecular structure is thus predicted to be planar at the energy minimum, which is associated with the trans conformation, and to become markedly twisted at the minimum of higher energy, which is associated with the cis conformation. The results of the calculations performed with the post-HF methods approach those obtained with the DFT methods, provided that the level of treatment of the electronic correlation is high enough and that sufficiently flexible basis sets are used. Calculations carried out within DFT also allowed the determination of the geometry and the energy of the molecule at the biradicaloid transition state associated with the thermal cis=trans isomerization and at the transition states associated with the enantiomerization of the cis isomer and with the rotations of the pyridinyl and phenyl groups in the trans and cis isomers. Car-Parrinello molecular dynamics simulations were also performed at 50, 150, and 300 K using the PBE functional. The studies allowed us to evidence the highly flexible nature of the molecule in both conformations. In particular, the trans isomer was found to exist mainly in a nonplanar form at finite temperatures, while the rotation of the pyridinyl ring in the cis isomer was incidentally observed to take place within approximately 1 ps during the simulation carried out at 150 K on this isomer.
Collapse
Affiliation(s)
- Latévi M Lawson Daku
- Département de chimie physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland.
| | | | | |
Collapse
|
45
|
Philipsen PHT, Baerends EJ. Role of the Fermi Surface in Adsorbate−Metal Interactions: An Energy Decomposition Analysis. J Phys Chem B 2006; 110:12470-9. [PMID: 16800574 DOI: 10.1021/jp060886e] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We present the result of a fragment-based energy decomposition analysis on some molecule-surface interactions. The analysis allows us to quantify the Pauli repulsion, its relief, and the attractive orbital interaction energy. In a metal, the existence of incompletely occupied energy bands causes significant relief of the Pauli repulsion due to escape of antibonding electrons to unoccupied states at the Fermi energy. This is the key electronic structure feature of metals that causes metal-molecule bond energies to be stronger and dissociation barriers of chemisorbed molecules to be much lower than those in comparable systems with no or one metal atom. As examples, we discuss the energy decomposition for the activated dissociation of hydrogen on the Cu surface and its unactivated dissociation on Pd, and for the (activated) chemisorption of N2 on W. We show that in all cases the relief of Pauli repulsion is of crucial importance for the chemisorption energy and for the low (or nonexistent) dissociation barriers. The barrier to the chemisorption well for nitrogen on tungsten is clearly related to a late relief of the Pauli repulsion. The relief of Pauli repulsion is important in lowering the barrier to dissociation of H2 on both Cu and Pd, but the difference in barrier heights for Cu and Pd appears to not be due to stronger relief of Pauli repulsion on Pd but primarily to the Pauli repulsion itself being stronger on Cu than on Pd, the relief energy being quite comparable on the two metals.
Collapse
Affiliation(s)
- P H T Philipsen
- Department of Chemistry, Section Theoretical Chemistry, Vrije Universiteit, Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | | |
Collapse
|
46
|
Fuchs M, Niquet YM, Gonze X, Burke K. Describing static correlation in bond dissociation by Kohn–Sham density functional theory. J Chem Phys 2005; 122:094116. [PMID: 15836121 DOI: 10.1063/1.1858371] [Citation(s) in RCA: 135] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We show that density functional theory within the RPA (random phase approximation for the exchange-correlation energy) provides a correct description of bond dissociation in H(2) in a spin-restricted Kohn-Sham formalism, i.e., without artificial symmetry breaking. We present accurate adiabatic connection curves both at equilibrium and beyond the Coulson-Fisher point. The strong curvature at large bond length implies important static (left-right) correlation, justifying modern hybrid functional constructions but also demonstrating their limitations. Although exact at infinite separation and accurate near the equilibrium bond length, the RPA dissociation curve displays unphysical repulsion at larger but finite bond lengths. Going beyond the RPA by including the exact exchange kernel (RPA+X), we find a similar repulsion. We argue that this deficiency is due to the absence of double excitations in adiabatic linear response theory. Further analyzing the H(2) dissociation limit we show that the RPA+X is not size consistent, in contrast to the RPA.
Collapse
Affiliation(s)
- M Fuchs
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | | | | | | |
Collapse
|
47
|
Hernández * MI, Campos-Martínez J, Van Caillie C, Di Domenico † D. Wave packet calculations for H2+H2collisions: isotopic substitution effects. Mol Phys 2004. [DOI: 10.1080/00268970412331294793] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
48
|
Andersson S, Grüning M. Performance of Density Functionals for Calculating Barrier Heights of Chemical Reactions Relevant to Astrophysics. J Phys Chem A 2004. [DOI: 10.1021/jp040448c] [Citation(s) in RCA: 77] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Stefan Andersson
- Leiden Observatory, P.O. Box 9513, 2300 RA Leiden, The Netherlands, Theoretical Chemistry, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands, and Section Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Myrta Grüning
- Leiden Observatory, P.O. Box 9513, 2300 RA Leiden, The Netherlands, Theoretical Chemistry, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands, and Section Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
49
|
Grüning M, Gritsenko OV, Baerends EJ. Improved Description of Chemical Barriers with Generalized Gradient Approximations (GGAs) and Meta-GGAs. J Phys Chem A 2004. [DOI: 10.1021/jp031227u] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- M. Grüning
- Section Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - O. V. Gritsenko
- Section Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - E. J. Baerends
- Section Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
50
|
Khait YG, Hoffmann MR. Multireference spin-adapted variant of density functional theory. J Chem Phys 2004; 120:5005-16. [PMID: 15267366 DOI: 10.1063/1.1646352] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new Kohn-Sham formalism is developed for studying the lowest molecular electronic states of given space and spin symmetry whose densities are represented by weighted sums of several reference configurations. Unlike standard spin-density functional theory, the new formalism uses total spin conserving spin-density operators and spin-invariant density matrices so that the method is fully spin-adapted and solves the so-called spin-symmetry dilemma. The formalism permits the use of an arbitrary set of reference (noninteracting) configurations with any number of open shells. It is shown that the requirement of degeneracy of the total noninteracting energies of the reference configurations (or configuration state functions) is equivalent to the stationary condition of the exact energy relative to the weights of the configurations (or configuration state functions). Consequently, at any molecular geometry, the weights can be determined by minimization of the energy, and, for given reference weights, the Kohn-Sham orbitals can be determined. From this viewpoint, the developed theory can be interpreted as an analog of the multiconfiguration self-consistent field approach within density functional theory.
Collapse
Affiliation(s)
- Yuriy G Khait
- Chemistry Department, University of North Dakota, Grand Forks, North Dakota 58202-9024, USA
| | | |
Collapse
|