1
|
Bai S, Zhang P, Beratan DN. Using Adiabatic Energy Splitting To Compute Dexter Energy Transfer Couplings. J Phys Chem A 2024; 128:5721-5729. [PMID: 38968620 DOI: 10.1021/acs.jpca.3c08146] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/07/2024]
Abstract
Dexter energy transfer and transport (DET) are of broad interest in energy science, and DET rates depend on electronic couplings between donor and acceptor species. DET couplings are challenging to compute since they originate from both one- and two-particle interactions, and the strength of this interaction drops approximately exponentially with donor-acceptor distances. Using adiabatic energy splitting to compute DET couplings has advantages because adiabatic states can be calculated directly using conventional quantum chemical methods. We describe a minimum energy splitting method to compute the DET coupling by altering molecular geometries to drive the systems into a T1/T2 energy quasi-degenerate-activated DA complex. We explore the accuracy of various quantum chemical approaches to calculate the Dexter couplings.
Collapse
Affiliation(s)
- Shuming Bai
- 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, China
- University of Chinese Academy of Sciences, Beijing 100049, China
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Peng Zhang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - David N Beratan
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
- Department of Biochemistry, Duke University, Durham, North Carolina 27710, United States
- Department of Physics, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
2
|
Humeniuk A. Approximate Functionals for Multistate Density Functional Theory. J Chem Theory Comput 2024; 20:5497-5509. [PMID: 38905701 DOI: 10.1021/acs.jctc.4c00330] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/23/2024]
Abstract
Recently Lu and Gao [J. Phys. Chem. Lett. 2022, 13 (33), 7762-7769] published a new, rigorous density functional theory for excited states and proved that the projection of the kinetic and electron-repulsion operators into the subspace of the lowest electronic states are a universal functional of the matrix density D(r). This is the first attempt to find an approximation to the multistate universal functional F [ D ( r ) ] . It is shown that F (i) does not explicitly depend on the number of electronic states and (ii) is an analytic matrix functional. The Thomas-Fermi-Dirac-von Weizsäcker model and the correlation energy of the homogeneous electron gas are turned into matrix functionals guided by two principles: that each matrix functional should transform properly under basis set transformations and that the ground state functional should be recovered for a single electronic state. Lieb-Oxford-like bounds on the average kinetic and electron-repulsion energies in the subspace are given. When evaluated on the numerically exact matrix density of LiF, this simple approximation reproduces the matrix elements of the electron-repulsion operator in the basis of the exact eigenstates accurately for all bond lengths. In particular the off-diagonal elements of the effective Hamiltonian that come from the interactions of different electronic states can be calculated with the same or better accuracy than the diagonal elements. Unsurprisingly, the largest error comes from the kinetic energy functional. More exact conditions that constrain the functional form of F are needed to go beyond the local density approximation.
Collapse
Affiliation(s)
- Alexander Humeniuk
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-Ku, Kyoto 606-8502, Japan
| |
Collapse
|
3
|
Qiu T, Bian X, Tao Z, Subotnik JE. A fast and smooth one-electron approach for investigating charge transfer states and D1-D0 crossings for systems with odd numbers of electrons. J Chem Phys 2024; 160:214115. [PMID: 38832731 DOI: 10.1063/5.0209036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2024] [Accepted: 05/13/2024] [Indexed: 06/05/2024] Open
Abstract
We propose an efficient version of ensemble Hartree-Fock/density functional theory to calculate a set of two charge-transfer states for systems with odd-numbers of electrons. The approach does require definitions of donor/acceptor fragments; however, the approach is not very sensitive to such definitions-even in the limit of very strong electronic coupling. The key ansatz is that, by mandating that the vector space spanned by the active orbitals projects equally onto the donor and acceptor fragments, such a constraint eliminates all intra-molecular local excitations and makes it far easier to generate potential energy surfaces that are smooth over a wide region of configuration space. The method is fast, working with only two electron configurations, and should be useful for ab initio non-adiabatic dynamics in the near future.
Collapse
Affiliation(s)
- Tian Qiu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Zhen Tao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
4
|
Zhong J, Soudackov AV, Hammes-Schiffer S. Probing Nonadiabaticity of Proton-Coupled Electron Transfer in Ribonucleotide Reductase. J Phys Chem Lett 2024; 15:1686-1693. [PMID: 38315651 DOI: 10.1021/acs.jpclett.3c03552] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2024]
Abstract
The enzyme ribonucleotide reductase, which is essential for DNA synthesis, initiates the conversion of ribonucleotides to deoxyribonucleotides via radical transfer over a 32 Å pathway composed of proton-coupled electron transfer (PCET) reactions. Previously, the first three PCET reactions in the α subunit were investigated with hybrid quantum mechanical/molecular mechanical (QM/MM) free energy simulations. Herein, the fourth PCET reaction in this subunit between C439 and guanosine diphosphate (GDP) is simulated and found to be slightly exoergic with a relatively high free energy barrier. To further elucidate the mechanisms of all four PCET reactions, we analyzed the vibronic and electron-proton nonadiabaticities. This analysis suggests that interfacial PCET between Y356 and Y731 is vibronically and electronically nonadiabatic, whereas PCET between Y731 and Y730 and between C439 and GDP is fully adiabatic and PCET between Y730 and C439 is in the intermediate regime. These insights provide guidance for selecting suitable rate constant expressions for these PCET reactions.
Collapse
Affiliation(s)
- Jiayun Zhong
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Alexander V Soudackov
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
5
|
Wu X, Cao C, Zhou C, Wu W. Hybrid Density Functional Valence Bond Method with Multistate Treatment. J Chem Theory Comput 2024. [PMID: 38279919 DOI: 10.1021/acs.jctc.3c01170] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2024]
Abstract
Recently, a hybrid density functional valence bond (VB) method, λ-DFVB(U), has been proposed and shown to give accuracy that is comparable to that of CASPT2 in calculations of atomization energies, atomic excitation energies, and reaction barriers, while its computational cost is approximately the same as the valence bond self-consistent-field (VBSCF) method. However, the interaction between electronic states is not included in λ-DFVB(U) since the last step of λ-DFVB(U) is not a diagonalization of the Hamiltonian matrix on the electronic state basis. Therefore, λ-DFVB(U) gives the wrong topology of the potential energy surfaces (PESs) near the conical intersection region. In the present paper, we propose a novel hybrid density functional VB method with multistate treatment, named λ-DFVB(MS), in which an effective Hamiltonian matrix is constructed on the basis of the diabatic states obtained by the valence-bond-based compression approach for the diabatization scheme, and the interaction between electronic states can be included through the diagonalization of the effective Hamiltonian matrix. Test calculations show that λ-DFVB(MS) gives the correct topology of the PESs near the conical intersection region. We also show that the VBSCF wave function with selected VB structures can be applied as a reference in λ-DFVB(MS).
Collapse
Affiliation(s)
- Xun Wu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen 361005, China
| | - Chan Cao
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen 361005, China
| | - Chen Zhou
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen 361005, China
| | - Wei Wu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen 361005, China
| |
Collapse
|
6
|
Rapacioli M, Buey MY, Spiegelman F. Addressing electronic and dynamical evolution of molecules and molecular clusters: DFTB simulations of energy relaxation in polycyclic aromatic hydrocarbons. Phys Chem Chem Phys 2024; 26:1499-1515. [PMID: 37933901 PMCID: PMC10793726 DOI: 10.1039/d3cp02852f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2023] [Accepted: 10/17/2023] [Indexed: 11/08/2023]
Abstract
We present a review of the capabilities of the density functional based Tight Binding (DFTB) scheme to address the electronic relaxation and dynamical evolution of molecules and molecular clusters following energy deposition via either collision or photoabsorption. The basics and extensions of DFTB for addressing these systems and in particular their electronic states and their dynamical evolution are reviewed. Applications to PAH molecules and clusters, carbonaceous systems of major interest in astrochemical/astrophysical context, are reported. A variety of processes are examined and discussed such as collisional hydrogenation, fast collisional processes and induced electronic and charge dynamics, collision-induced fragmentation, photo-induced fragmentation, relaxation in high electronic states, electronic-to-vibrational energy conversion and statistical versus non-statistical fragmentation. This review illustrates how simulations may help to unravel different relaxation mechanisms depending on various factors such as the system size, specific electronic structure or excitation conditions, in close connection with experiments.
Collapse
Affiliation(s)
- Mathias Rapacioli
- Laboratoire de Chimie et Physique Quantique (LCPQ/FERMI), UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France.
| | - Maysa Yusef Buey
- Laboratoire de Chimie et Physique Quantique (LCPQ/FERMI), UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France.
| | - Fernand Spiegelman
- Laboratoire de Chimie et Physique Quantique (LCPQ/FERMI), UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France.
| |
Collapse
|
7
|
Melčák M, Šebesta F, Heyda J, Gray HB, Záliš S, Vlček A. Tryptophan to Tryptophan Hole Hopping in an Azurin Construct. J Phys Chem B 2024; 128:96-108. [PMID: 38145895 PMCID: PMC10788906 DOI: 10.1021/acs.jpcb.3c06568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2023] [Revised: 11/18/2023] [Accepted: 11/30/2023] [Indexed: 12/27/2023]
Abstract
Electron transfer (ET) between neutral and cationic tryptophan residues in the azurin construct [ReI(H126)(CO)3(dmp)](W124)(W122)CuI (dmp = 4,7-Me2-1,10-phenanthroline) was investigated by Born-Oppenheimer quantum-mechanics/molecular mechanics/molecular dynamics (QM/MM/MD) simulations. We focused on W124•+ ← W122 ET, which is the middle step of the photochemical hole-hopping process *ReII(CO)3(dmp•-) ← W124 ← W122 ← CuI, where sequential hopping amounts to nearly 10,000-fold acceleration over single-step tunneling (ACS Cent. Sci. 2019, 5, 192-200). In accordance with experiments, UKS-DFT QM/MM/MD simulations identified forward and reverse steps of W124•+ ↔ W122 ET equilibrium, as well as back ET ReI(CO)3(dmp•-) → W124•+ that restores *ReII(CO)3(dmp•-). Strong electronic coupling between the two indoles (≥40 meV in the crossing region) makes the productive W124•+ ← W122 ET adiabatic. Energies of the two redox states are driven to degeneracy by fluctuations of the electrostatic potential at the two indoles, mainly caused by water solvation, with contributions from the protein dynamics in the W122 vicinity. ET probability depends on the orientation of Re(CO)3(dmp) relative to W124 and its rotation diminishes the hopping yield. Comparison with hole hopping in natural systems reveals structural and dynamics factors that are important for designing efficient hole-hopping processes.
Collapse
Affiliation(s)
- Martin Melčák
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, CZ-182 23 Prague, Czech Republic
- Department
of Physical Chemistry, University of Chemistry
and Technology Prague, Technická 5, CZ-166 28 Prague, Czech Republic
| | - Filip Šebesta
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, CZ-182 23 Prague, Czech Republic
- Department
of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-121 16 Prague, Czech Republic
| | - Jan Heyda
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, CZ-182 23 Prague, Czech Republic
- Department
of Physical Chemistry, University of Chemistry
and Technology Prague, Technická 5, CZ-166 28 Prague, Czech Republic
| | - Harry B. Gray
- Beckman
Institute, California Institute of Technology, Pasadena, California 91125, United States
| | - Stanislav Záliš
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, CZ-182 23 Prague, Czech Republic
| | - Antonín Vlček
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, CZ-182 23 Prague, Czech Republic
- Department
of Chemistry, Queen Mary University of London, London E1 4NS, U.K.
| |
Collapse
|
8
|
Kempfer-Robertson EM, Mahler AD, Haase MN, Roe P, Thompson LM. Nonorthogonal Active Space Decomposition of Wave Functions with Multiple Correlation Mechanisms. J Phys Chem Lett 2022; 13:12041-12048. [PMID: 36541869 DOI: 10.1021/acs.jpclett.2c03349] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The nonorthogonal active space decomposition (NO-ASD) methodology is proposed for describing systems containing multiple correlation mechanisms. NO-ASD partitions the wave function by a correlation mechanism, such that the interactions between different correlation mechanisms are treated with an effective Hamiltonian approach, while interactions between correlated orbitals in the same correlation mechanism are treated explicitly. As a result, the determinant expansion scales polynomially with the number of correlation mechanisms rather than exponentially, which significantly reduces the factorial scaling associated with the size of the correlated orbital space. Despite the nonorthogonal framework of NO-ASD, the approach can take advantage of computational efficient matrix element evaluation when performing nonorthogonal coupling of orthogonal determinant expansions. In this work, we introduce and examine the NO-ASD approach in comparison to complete active space methods to establish how the NO-ASD approach reduces the problem dimensionality and the extent to which it affects the amount of correlation energy recovered. Calculations are performed on ozone, nickel-acetylene, and isomers of μ-oxo dicopper ammonia.
Collapse
Affiliation(s)
| | - Andrew D Mahler
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| | - Meagan N Haase
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| | - Piper Roe
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| | - Lee M Thompson
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| |
Collapse
|
9
|
Athavale V, Teh HH, Shao Y, Subotnik J. Analytical gradients and derivative couplings for the TDDFT-1D method. J Chem Phys 2022; 157:244110. [PMID: 36586994 DOI: 10.1063/5.0130404] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
Abstract
We derive and implement analytic gradients and derivative couplings for time-dependent density functional theory plus one double (TDDFT-1D) which is a semiempirical configuration interaction method whereby the Hamiltonian is diagonalized in a basis of all singly excited configurations and one doubly excited configuration as constructed from a set of reference Kohn-Sham orbitals. We validate the implementation by comparing against finite difference values. Furthermore, we show that our implementation can locate both optimized geometries and minimum-energy crossing points along conical seams of S1/S0 surfaces for a set of test cases.
Collapse
Affiliation(s)
- Vishikh Athavale
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Hung-Hsuan Teh
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Joseph Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
10
|
Niemeyer N, Eschenbach P, Bensberg M, Tölle J, Hellmann L, Lampe L, Massolle A, Rikus A, Schnieders D, Unsleber JP, Neugebauer J. The subsystem quantum chemistry program
Serenity. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1647] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Niklas Niemeyer
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Patrick Eschenbach
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Moritz Bensberg
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Johannes Tölle
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Lars Hellmann
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Lukas Lampe
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Anja Massolle
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Anton Rikus
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - David Schnieders
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Jan P. Unsleber
- Laboratorium für Physikalische Chemie ETH Zürich Zürich Switzerland
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| |
Collapse
|
11
|
Eschenbach P, Neugebauer J. Subsystem density-functional theory: A reliable tool for spin-density based properties. J Chem Phys 2022; 157:130902. [PMID: 36209003 DOI: 10.1063/5.0103091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Subsystem density-functional theory compiles a set of features that allow for efficiently calculating properties of very large open-shell radical systems such as organic radical crystals, proteins, or deoxyribonucleic acid stacks. It is computationally less costly than correlated ab initio wave function approaches and can pragmatically avoid the overdelocalization problem of Kohn-Sham density-functional theory without employing hard constraints on the electron-density. Additionally, subsystem density-functional theory calculations commonly start from isolated fragment electron densities, pragmatically preserving a priori specified subsystem spin-patterns throughout the calculation. Methods based on subsystem density-functional theory have seen a rapid development over the past years and have become important tools for describing open-shell properties. In this Perspective, we address open questions and possible developments toward challenging future applications in connection with subsystem density-functional theory for spin-dependent properties.
Collapse
Affiliation(s)
- Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
12
|
Glebov IO, Poddubnyy VV, Khokhlov DV. Perturbative Expansion of Nonorthogonal Product Approach for Charge Transfer States. J Phys Chem A 2022; 126:5800-5813. [PMID: 35998639 DOI: 10.1021/acs.jpca.2c02612] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Modeling of the excited states of multichromophoric systems is crucial for the understanding of photosynthesis functioning. The excitonic Hamiltonian method is widely used for such calculations. Excited states of the combined system are constructed from the wave functions of individual chromophores while their interactions are described by excitonic couplings. In the current study, we enhance a previously proposed nonorthogonal product approach to incorporate dynamic correlation effects accounted for by the multireference perturbation theory. We discuss the problems of constructing the excitonic Hamiltonian including charge transfer states for the molecular systems where the overlap contribution to the excitonic couplings is non-negligible. The benchmark calculations were performed for a model system. It was shown that the overlap component of the excitonic coupling is of great importance. The enhanced method provides an accurate description of the excited state energies and other properties.
Collapse
Affiliation(s)
- Ilya O Glebov
- Chemistry Department, Moscow State University, Leninskie gory 1-3, Moscow 119991, Russia
| | - Vladimir V Poddubnyy
- Chemistry Department, Moscow State University, Leninskie gory 1-3, Moscow 119991, Russia
| | - Daniil V Khokhlov
- Chemistry Department, Moscow State University, Leninskie gory 1-3, Moscow 119991, Russia
| |
Collapse
|
13
|
Ahart CS, Rosso KM, Blumberger J. Implementation and Validation of Constrained Density Functional Theory Forces in the CP2K Package. J Chem Theory Comput 2022; 18:4438-4446. [PMID: 35700315 PMCID: PMC9281399 DOI: 10.1021/acs.jctc.2c00284] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Constrained density
functional theory (CDFT) is a powerful tool
for the prediction of electron transfer parameters in condensed phase
simulations at a reasonable computational cost. In this work we present
an extension to CDFT in the popular mixed Gaussian/plane wave electronic
structure package CP2K, implementing the additional force terms arising
from a constraint based on Hirshfeld charge partitioning. This improves
upon the existing Becke partitioning scheme, which is prone to give
unphysical atomic charges. We verify this implementation for a variety
of systems: electron transfer in (H2O)2+ in a vacuum, electron tunnelling
between oxygen vacancy centers in solid MgO, and electron self-exchange
in aqueous Ru2+–Ru3+. We find good agreement
with previous plane-wave CDFT results for the same systems, but at
a significantly lower computational cost, and we discuss the general
reliability of condensed phase CDFT calculations.
Collapse
Affiliation(s)
- Christian S Ahart
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Kevin M Rosso
- Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
14
|
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
|
15
|
Grofe A, Li X. Relativistic nonorthogonal configuration interaction: application to L 2,3-edge X-ray spectroscopy. Phys Chem Chem Phys 2022; 24:10745-10756. [PMID: 35451435 DOI: 10.1039/d2cp01127a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
In this article, we develop a relativistic exact-two-component nonorthogonal configuration interaction (X2C-NOCI) for computing L-edge X-ray spectra. This article to our knowledge is the first time NOCI has been used for relativistic wave functions. A set of molecular complexes, including SF6, SiCl4 and [FeCl6]3-, are used to demonstrate the accuracy and computational scaling of the X2C-NOCI method. Our results suggest that X2C-NOCI is able to satisfactorily capture the main features of the L2,3-edge X-ray absorption spectra. Excitations from the core require a large amount of orbital relaxation to yield reasonable energies and X2C-NOCI allows us to treat orbital optimization explicitly. However, the cost of computing the nonorthogonal coupling is higher than in conventional CI. Here, we propose an improved integral screening using overlap-scaled density combined with a continuous measure of the generalized Slater-Condon rules that allows us to estimate if an element is zero before attempting a two-electron integral contraction.
Collapse
Affiliation(s)
- Adam Grofe
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA.
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA.
| |
Collapse
|
16
|
Nishio S, Kurashige Y. Importance of dynamical electron correlation in diabatic couplings of electron-exchange processes. J Chem Phys 2022; 156:114107. [PMID: 35317578 DOI: 10.1063/5.0075978] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We demonstrate the importance of the dynamical electron correlation effect in diabatic couplings of electron-exchange processes in molecular aggregates. To perform a multireference perturbation theory with large active space of molecular aggregates, an efficient low-rank approximation is applied to the complete active space self-consistent field reference functions. It is known that kinetic rates of electron-exchange processes, such as singlet fission, triplet-triplet annihilation, and triplet exciton transfer, are not sufficiently explained by the direct term of the diabatic couplings but efficiently mediated by the low-lying charge transfer states if the two molecules are in close proximity. It is presented in this paper, however, that regardless of the distance of the molecules, the direct term is considerably underestimated by up to three orders of magnitude without the dynamical electron correlation, i.e., the diabatic states expressed in the active space are not adequate to quantitatively reproduce the electron-exchange processes.
Collapse
Affiliation(s)
- Soichiro Nishio
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
17
|
Warburton RE, Soudackov AV, Hammes-Schiffer S. Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer. Chem Rev 2022; 122:10599-10650. [PMID: 35230812 DOI: 10.1021/acs.chemrev.1c00929] [Citation(s) in RCA: 42] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. A vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential-dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.
Collapse
Affiliation(s)
- Robert E Warburton
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Alexander V Soudackov
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
18
|
Ren M, Zhang L, Jiao Y, Chen Z, Wu W. Extended Mulliken-Hush Method with Applications to the Theoretical Study of Electron Transfer. J Chem Theory Comput 2021; 17:6861-6875. [PMID: 34605634 DOI: 10.1021/acs.jctc.1c00603] [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 novel adiabatic-to-diabatic (ATD) transformation strategy, namely, the extended Mulliken-Hush (XMH) method, is proposed to evaluate diabatic properties including electronic couplings, potential energy surfaces, and their crossings. The XMH method is developed by adopting our recently proposed ATD transformation formula of a general vectorial physical observable, in which a useful ATD transformation is further determined by using an auxiliary dipole between localized frontier orbitals as a simple approximation of the diabatic transition dipole. The XMH method is simple and practical that provides a flexible way to construct diabatic states. To some extent, it can be regarded as an extension of the generalized Mulliken-Hush (GMH) method since the latter takes a stronger approximation, in which the diabatic transition dipole is assumed to be vanishing. Test calculations on the HeH2+ system show that the electronic couplings predicted by the XMH method are closer to the ones calculated by the valence bond block-diagonalization approach than the GMH ones since the XMH method takes into account both the magnitude and direction of the diabatic transition dipole, which is consistent with the properties of this molecule. In the study of electron transfer in the two kinds of donor-bridge-acceptor systems, the XMH method maintains the simplicity of the GMH method and gives reasonable results even when the latter fails, wherein the diabatic transition dipole is nearly perpendicular to the difference of the initial and final adiabatic dipoles. More importantly, the XMH method can be easily combined with high-level electronic structure methods, in which the properties of the ground and excited states may be more accurately calculated, and hence, one may expect that further development of the XMH method would result in a general computational model for studying electron transfer reactions.
Collapse
Affiliation(s)
- Mingxing Ren
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Lina Zhang
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Yang Jiao
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Zhenhua Chen
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Wei Wu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
19
|
Herbert JM, Paul SK. Interaction Energy Analysis of Monovalent Inorganic Anions in Bulk Water Versus Air/Water Interface. Molecules 2021; 26:6719. [PMID: 34771133 PMCID: PMC8588468 DOI: 10.3390/molecules26216719] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2021] [Revised: 10/29/2021] [Accepted: 11/03/2021] [Indexed: 11/16/2022] Open
Abstract
Soft anions exhibit surface activity at the air/water interface that can be probed using surface-sensitive vibrational spectroscopy, but the structural implications of this surface activity remain a matter of debate. Here, we examine the nature of anion-water interactions at the air/water interface using a combination of molecular dynamics simulations and quantum-mechanical energy decomposition analysis based on symmetry-adapted perturbation theory. Results are presented for a set of monovalent anions, including Cl-, Br-, I-, CN-, OCN-, SCN-, NO2-, NO3-, and ClOn- (n=1,2,3,4), several of which are archetypal examples of surface-active species. In all cases, we find that average anion-water interaction energies are systematically larger in bulk water although the difference (with respect to the same quantity computed in the interfacial environment) is well within the magnitude of the instantaneous fluctuations. Specifically for the surface-active species Br-(aq), I-(aq), ClO4-(aq), and SCN-(aq), and also for ClO-(aq), the charge-transfer (CT) energy is found to be larger at the interface than it is in bulk water, by an amount that is greater than the standard deviation of the fluctuations. The Cl-(aq) ion has a slightly larger CT energy at the interface, but NO3-(aq) does not; these two species are borderline cases where consensus is lacking regarding their surface activity. However, CT stabilization amounts to <20% of the total induction energy for each of the ions considered here, and CT-free polarization energies are systematically larger in bulk water in all cases. As such, the role of these effects in the surface activity of soft anions remains unclear. This analysis complements our recent work suggesting that the short-range solvation structure around these ions is scarcely different at the air/water interface from what it is in bulk water. Together, these observations suggest that changes in first-shell hydration structure around soft anions cannot explain observed surface activities.
Collapse
Affiliation(s)
- John M. Herbert
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, OH 43210, USA;
| | | |
Collapse
|
20
|
Practical treatment of singlet oxygen with density-functional theory and the multiplet-sum method. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02852-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
21
|
Zhao R, Grofe A, Wang Z, Bao P, Chen X, Liu W, Gao J. Dynamic-then-Static Approach for Core Excitations of Open-Shell Molecules. J Phys Chem Lett 2021; 12:7409-7417. [PMID: 34328742 DOI: 10.1021/acs.jpclett.1c02039] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Delta self-consistent-field methods are widely used in studies of electronically excited states. However, the nonaufbau determinants are generally spin-contaminated. Here, we describe a general approach for spin-coupling interactions of open-shell molecules, making use of multistate density functional theory (MSDFT). In particular, the effective exchange integrals that determine spin coupling are obtained by enforcing the multiplet degeneracy of the S+1 state in the MS = S manifold. Consequently, they are consistent with the energy of the high-spin state that is adequately treated by Kohn-Sham density functional theory (DFT) and, thereby, free of double counting of correlation. The method was applied to core excitations of open-shell molecules and compared with those by spin-adapted time-dependent DFT. An excellent agreement with experiment was found employing the BLYP functional and aug-cc-pCVQZ basis set. Overall, MSDFT provides an effective combination of the strengths of DFT and wave function theory to achieve efficiency and accuracy.
Collapse
Affiliation(s)
- Ruoqi Zhao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
- Institute of Theoretical Chemistry, Jilin University Changchun 130023, Jilin, China
| | - Adam Grofe
- Institute of Theoretical Chemistry, Jilin University Changchun 130023, Jilin, China
| | - Zikuan Wang
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, Shandong, China
| | - 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, China
| | - Xin Chen
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, Shandong, China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis 55455, Minnesota, United States
| |
Collapse
|
22
|
Matsika S. Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections. Chem Rev 2021; 121:9407-9449. [PMID: 34156838 DOI: 10.1021/acs.chemrev.1c00074] [Citation(s) in RCA: 59] [Impact Index Per Article: 19.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.
Collapse
Affiliation(s)
- Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
23
|
Bourne Worster S, Feighan O, Manby FR. Reliable transition properties from excited-state mean-field calculations. J Chem Phys 2021; 154:124106. [DOI: 10.1063/5.0041233] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Susannah Bourne Worster
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Oliver Feighan
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Frederick R. Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
24
|
Herbert JM, Carter-Fenk K. Electrostatics, Charge Transfer, and the Nature of the Halide-Water Hydrogen Bond. J Phys Chem A 2021; 125:1243-1256. [PMID: 33502859 DOI: 10.1021/acs.jpca.0c11356] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Binary halide-water complexes X-(H2O) are examined by means of symmetry-adapted perturbation theory, using charge-constrained promolecular reference densities to extract a meaningful charge-transfer component from the induction energy. As is known, the X-(H2O) potential energy surface (for X = F, Cl, Br, or I) is characterized by symmetric left and right hydrogen bonds separated by a C2v-symmetric saddle point, with a tunneling barrier height that is <2 kcal/mol except in the case of F-(H2O). Our analysis demonstrates that the charge-transfer energy is correspondingly small (<2 kcal/mol except for X = F), considerably smaller than the electrostatic interaction energy. Nevertheless, charge transfer plays a crucial role determining the conformational preferences of X-(H2O) and provides a driving force for the formation of quasi-linear X··· H-O hydrogen bonds. Charge-transfer energies correlate well with measured O-H vibrational redshifts for the halide-water complexes and also for OH-(H2O) and NO2-(H2O), providing some indication of a general mechanism.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
25
|
Zheng P, Ji C, Ying F, Su P, Wu W. A Valence-Bond-Based Multiconfigurational Density Functional Theory: The λ-DFVB Method Revisited. Molecules 2021; 26:521. [PMID: 33498268 PMCID: PMC7863953 DOI: 10.3390/molecules26030521] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Revised: 01/17/2021] [Accepted: 01/18/2021] [Indexed: 11/29/2022] Open
Abstract
A recently developed valence-bond-based multireference density functional theory, named λ-DFVB, is revisited in this paper. λ-DFVB remedies the double-counting error of electron correlation by decomposing the electron-electron interactions into the wave function term and density functional term with a variable parameter λ. The λ value is defined as a function of the free valence index in our previous scheme, denoted as λ-DFVB(K) in this paper. Here we revisit the λ-DFVB method and present a new scheme based on natural orbital occupation numbers (NOONs) for parameter λ, named λ-DFVB(IS), to simplify the process of λ-DFVB calculation. In λ-DFVB(IS), the parameter λ is defined as a function of NOONs, which are straightforwardly determined from the many-electron wave function of the molecule. Furthermore, λ-DFVB(IS) does not involve further self-consistent field calculation after performing the valence bond self-consistent field (VBSCF) calculation, and thus, the computational effort in λ-DFVB(IS) is approximately the same as the VBSCF method, greatly reduced from λ-DFVB(K). The performance of λ-DFVB(IS) was investigated on a broader range of molecular properties, including equilibrium bond lengths and dissociation energies, atomization energies, atomic excitation energies, and chemical reaction barriers. The computational results show that λ-DFVB(IS) is more robust without losing accuracy and comparable in accuracy to high-level multireference wave function methods, such as CASPT2.
Collapse
Affiliation(s)
| | | | | | - Peifeng Su
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, The State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China; (P.Z.); (C.J.); (F.Y.)
| | - Wei Wu
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, The State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China; (P.Z.); (C.J.); (F.Y.)
| |
Collapse
|
26
|
Barragan AM, Soudackov AV, Luthey-Schulten Z, Hammes-Schiffer S, Schulten K, Solov'yov IA. Theoretical Description of the Primary Proton-Coupled Electron Transfer Reaction in the Cytochrome bc1 Complex. J Am Chem Soc 2021; 143:715-723. [PMID: 33397104 DOI: 10.1021/jacs.0c07799] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The cytochrome bc1 complex is a transmembrane enzymatic protein complex that plays a central role in cellular energy production and is present in both photosynthetic and respiratory chain organelles. Its reaction mechanism is initiated by the binding of a quinol molecule to an active site, followed by a series of charge transfer reactions between the quinol and protein subunits. Previous work hypothesized that the primary reaction was a concerted proton-coupled electron transfer (PCET) reaction because of the apparent absence of intermediate states associated with single proton or electron transfer reactions. In the present study, the kinetics of the primary bc1 complex PCET reaction is investigated with a vibronically nonadiabatic PCET theory in conjunction with all-atom molecular dynamics simulations and electronic structure calculations. The computed rate constants and relatively high kinetic isotope effects are consistent with experimental measurements on related biomimetic systems. The analysis implicates a concerted PCET mechanism with significant hydrogen tunneling and nonadiabatic effects in the bc1 complex. Moreover, the employed theoretical framework is shown to serve as a general strategy for describing PCET reactions in bioenergetic systems.
Collapse
Affiliation(s)
- Angela M Barragan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 North Mathews Avenue, Urbana, Illinois 61801, United States.,Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States
| | - Alexander V Soudackov
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Zaida Luthey-Schulten
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 North Mathews Avenue, Urbana, Illinois 61801, United States.,Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States.,Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Klaus Schulten
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 North Mathews Avenue, Urbana, Illinois 61801, United States.,Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States
| | - Ilia A Solov'yov
- Department of Physics, Carl von Ossietzky University, Carl-von-Ossietzky-Strasse 9-11, Oldenburg 26111, Germany
| |
Collapse
|
27
|
Coleman CN, Tapping PC, Huxley MT, Kee TW, Huang DM, Doonan CJ, Sumby CJ. Structural modulation of the photophysical and electronic properties of pyrene-based 3D metal–organic frameworks derived from s-block metals. CrystEngComm 2021. [DOI: 10.1039/d0ce01505a] [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
Materials in which charge delocalization and migration can be tuned are critical for electronic applications.
Collapse
Affiliation(s)
- Christopher N. Coleman
- Department of Chemistry and Centre for Advanced Nanomaterials
- The University of Adelaide
- Adelaide
- Australia
| | - Patrick C. Tapping
- Department of Chemistry and Centre for Advanced Nanomaterials
- The University of Adelaide
- Adelaide
- Australia
| | - Michael T. Huxley
- Department of Chemistry and Centre for Advanced Nanomaterials
- The University of Adelaide
- Adelaide
- Australia
| | - Tak W. Kee
- Department of Chemistry and Centre for Advanced Nanomaterials
- The University of Adelaide
- Adelaide
- Australia
| | - David M. Huang
- Department of Chemistry and Centre for Advanced Nanomaterials
- The University of Adelaide
- Adelaide
- Australia
| | - Christian J. Doonan
- Department of Chemistry and Centre for Advanced Nanomaterials
- The University of Adelaide
- Adelaide
- Australia
| | - Christopher J. Sumby
- Department of Chemistry and Centre for Advanced Nanomaterials
- The University of Adelaide
- Adelaide
- Australia
| |
Collapse
|
28
|
Bao P, Hettich CP, Shi Q, Gao J. Block-Localized Excitation for Excimer Complex and Diabatic Coupling. J Chem Theory Comput 2020; 17:240-254. [PMID: 33370101 DOI: 10.1021/acs.jctc.0c01015] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We describe a block-localized excitation (BLE) method to carry out constrained optimization of block-localized orbitals for constructing valence bond-like, diabatic excited configurations using multistate density functional theory (MSDFT). The method is an extension of the previous block-localized wave function method through a fragment-based ΔSCF approach to optimize excited determinants within a molecular complex. In BLE, both the number of electrons and the electronic spin of different fragments in a whole system can be constrained, whereas electrostatic, exchange, and polarization interactions among different blocks can be fully taken into account of. To avoid optimization collapse to unwanted states, a ΔSCF projection scheme and a maximum overlap of wave function approach have been presented. The method is illustrated by the excimer complex of two naphthalene molecules. With a minimum of eight spin-adapted configurational state functions, it was found that the inversion of La- and Lb- states near the optimal structure of the excimer complex is correctly produced, which is in quantitative agreement with DMRG-CASPT2 calculations and experiments. Trends in the computed transfer integrals associated with excited-state energy transfer both in the singlet and triplet states are discussed. The results suggest that MSDFT may be used as an efficient approach to treat intermolecular interactions in excited states with a minimal active space (MAS) for interpretation of the results and for dynamic simulations, although the selection of a small active space is often system dependent.
Collapse
Affiliation(s)
- 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, China
| | - Christian P Hettich
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Qiang Shi
- 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, China.,University of the Chinese Academy of Sciences, Beijing 100049, China
| | - Jiali Gao
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States.,Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China.,Beijing University Shenzhen Graduate School, Shenzhen 518055, China
| |
Collapse
|
29
|
Mao Y, Montoya-Castillo A, Markland TE. Excited state diabatization on the cheap using DFT: Photoinduced electron and hole transfer. J Chem Phys 2020; 153:244111. [PMID: 33380087 DOI: 10.1063/5.0035593] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Excited state electron and hole transfer underpin fundamental steps in processes such as exciton dissociation at photovoltaic heterojunctions, photoinduced charge transfer at electrodes, and electron transfer in photosynthetic reaction centers. Diabatic states corresponding to charge or excitation localized species, such as locally excited and charge transfer states, provide a physically intuitive framework to simulate and understand these processes. However, obtaining accurate diabatic states and their couplings from adiabatic electronic states generally leads to inaccurate results when combined with low-tier electronic structure methods, such as time-dependent density functional theory, and exorbitant computational cost when combined with high-level wavefunction-based methods. Here, we introduce a density functional theory (DFT)-based diabatization scheme that directly constructs the diabatic states using absolutely localized molecular orbitals (ALMOs), which we denote as Δ-ALMO(MSDFT2). We demonstrate that our method, which combines ALMO calculations with the ΔSCF technique to construct electronically excited diabatic states and obtains their couplings with charge-transfer states using our MSDFT2 scheme, gives accurate results for excited state electron and hole transfer in both charged and uncharged systems that underlie DNA repair, charge separation in donor-acceptor dyads, chromophore-to-solvent electron transfer, and singlet fission. This framework for the accurate and efficient construction of excited state diabats and evaluation of their couplings directly from DFT thus offers a route to simulate and elucidate photoinduced electron and hole transfer in large disordered systems, such as those encountered in the condensed phase.
Collapse
Affiliation(s)
- Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | | | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
30
|
Al‐Nu'airat J, Oluwoye I, Zeinali N, Altarawneh M, Dlugogorski BZ. Review of Chemical Reactivity of Singlet Oxygen with Organic Fuels and Contaminants. CHEM REC 2020; 21:315-342. [DOI: 10.1002/tcr.202000143] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Revised: 11/26/2020] [Indexed: 01/03/2023]
Affiliation(s)
- Jomana Al‐Nu'airat
- Murdoch University Discipline of Chemistry and Physics, College of Science, Health, Engineering and Education 90 South Street Murdoch WA 6150 Australia
| | - Ibukun Oluwoye
- Murdoch University Discipline of Chemistry and Physics, College of Science, Health, Engineering and Education 90 South Street Murdoch WA 6150 Australia
| | - Nassim Zeinali
- Murdoch University Discipline of Chemistry and Physics, College of Science, Health, Engineering and Education 90 South Street Murdoch WA 6150 Australia
| | - Mohammednoor Altarawneh
- United Arab Emirates University Chemical and Petroleum Engineering Department Sheikh Khalifa bin Zayed St Al-Ain 15551 United Arab Emirates
| | - Bogdan Z. Dlugogorski
- Charles Darwin University Energy and Resources Institute, Ellengowan Drive Darwin NT 0909 Australia
| |
Collapse
|
31
|
Zarotiadis RA, Burton HGA, Thom AJW. Towards a Holomorphic Density Functional Theory. J Chem Theory Comput 2020; 16:7400-7412. [PMID: 33211476 DOI: 10.1021/acs.jctc.0c00822] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Self-consistent-field (SCF) approximations formulated using Hartree-Fock (HF) or Kohn-Sham density-functional theory (KS-DFT) have the potential to yield multiple solutions. However, the formal relationship between multiple solutions identified using HF or KS-DFT remains generally unknown. We investigate the connection between multiple SCF solutions for HF or KS-DFT by introducing a parameterized functional that scales between the two representations. Using the hydrogen molecule and a model of electron transfer, we continuously map multiple solutions from the HF potential to a KS-DFT description. We discover that multiple solutions can coalesce and vanish as the functional changes, forming a direct analogy with the disappearance of real HF solutions along a change in molecular structure. To overcome this disappearance of solutions, we develop a complex-analytic extension of DFT-the "holomorphic DFT" approach-that allows every SCF stationary state to be analytically continued across all molecular structures and exchange-correlation functionals.
Collapse
Affiliation(s)
- Rhiannon A Zarotiadis
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Hugh G A Burton
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K.,Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ
| | - Alex J W Thom
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| |
Collapse
|
32
|
Dontot L, Spiegelman F, Zamith S, Rapacioli M. Dependence upon charge of the vibrational spectra of small Polycyclic Aromatic Hydrocarbon clusters: the example of pyrene. THE EUROPEAN PHYSICAL JOURNAL. D, ATOMIC, MOLECULAR, AND OPTICAL PHYSICS 2020; 74:216. [PMID: 33597829 PMCID: PMC7116754 DOI: 10.1140/epjd/e2020-10081-0] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/13/2020] [Revised: 07/18/2020] [Accepted: 09/01/2020] [Indexed: 06/12/2023]
Abstract
Infrared spectra are computed for neutral and cationic clusters of Polycyclic Aromatic Hydrocarbon clusters, namely( C 16 H 10 ) n = 1 , 4 ( 0 / + ) , using the Density Functional based Tight Binding scheme combined with a Configuration Interaction scheme (DFTB-CI) in the double harmonic approximation. Cross-comparison is carried out with DFT and simple DFTB. Similarly to the monomer cation, the IR spectra of cluster cations are characterized by a depletion of the intensity of the CH stretch modes around 3000 cm-1, with a weak revival for n = 3 and 4. The in-plane CCC modes in the region 1400-2000 cm-1 are enhanced while the CH bending modes in the range 700-1000 cm-1 are significantly weakened with respect to the monomer cation, in particular for n = 2. Finally, soft modes corresponding to diedral fluctuations of the monomers within the central stack of the ion structure, possibly mixed with monomer folding, are also observed in the region 70-120 cm-1.
Collapse
Affiliation(s)
- Léo Dontot
- Laboratoire de Chimie et Physique Quantique (LCPQ/IRSAMC), UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Fernand Spiegelman
- Laboratoire de Chimie et Physique Quantique (LCPQ/IRSAMC), UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Sébastien Zamith
- Laboratoire Collisions Agrégats Réactivité (LCAR/IRSAMC), UMR5589, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Mathias Rapacioli
- Laboratoire de Chimie et Physique Quantique (LCPQ/IRSAMC), UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| |
Collapse
|
33
|
Marzolf DR, McKenzie AM, O’Malley MC, Ponomarenko NS, Swaim CM, Brittain TJ, Simmons NL, Pokkuluri PR, Mulfort KL, Tiede DM, Kokhan O. Mimicking Natural Photosynthesis: Designing Ultrafast Photosensitized Electron Transfer into Multiheme Cytochrome Protein Nanowires. NANOMATERIALS (BASEL, SWITZERLAND) 2020; 10:E2143. [PMID: 33126541 PMCID: PMC7693585 DOI: 10.3390/nano10112143] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Revised: 10/23/2020] [Accepted: 10/26/2020] [Indexed: 05/02/2023]
Abstract
Efficient nanomaterials for artificial photosynthesis require fast and robust unidirectional electron transfer (ET) from photosensitizers through charge-separation and accumulation units to redox-active catalytic sites. We explored the ultrafast time-scale limits of photo-induced charge transfer between a Ru(II)tris(bipyridine) derivative photosensitizer and PpcA, a 3-heme c-type cytochrome serving as a nanoscale biological wire. Four covalent attachment sites (K28C, K29C, K52C, and G53C) were engineered in PpcA enabling site-specific covalent labeling with expected donor-acceptor (DA) distances of 4-8 Å. X-ray scattering results demonstrated that mutations and chemical labeling did not disrupt the structure of the proteins. Time-resolved spectroscopy revealed three orders of magnitude difference in charge transfer rates for the systems with otherwise similar DA distances and the same number of covalent bonds separating donors and acceptors. All-atom molecular dynamics simulations provided additional insight into the structure-function requirements for ultrafast charge transfer and the requirement of van der Waals contact between aromatic atoms of photosensitizers and hemes in order to observe sub-nanosecond ET. This work demonstrates opportunities to utilize multi-heme c-cytochromes as frameworks for designing ultrafast light-driven ET into charge-accumulating biohybrid model systems, and ultimately for mimicking the photosynthetic paradigm of efficiently coupling ultrafast, light-driven electron transfer chemistry to multi-step catalysis within small, experimentally versatile photosynthetic biohybrid assemblies.
Collapse
Affiliation(s)
- Daniel R. Marzolf
- Department of Chemistry and Biochemistry, James Madison University, Harrisonburg, VA 22807, USA; (D.R.M.); (A.M.M.); (C.M.S.); (T.J.B.)
| | - Aidan M. McKenzie
- Department of Chemistry and Biochemistry, James Madison University, Harrisonburg, VA 22807, USA; (D.R.M.); (A.M.M.); (C.M.S.); (T.J.B.)
| | - Matthew C. O’Malley
- Department of Chemistry and Biochemistry, James Madison University, Harrisonburg, VA 22807, USA; (D.R.M.); (A.M.M.); (C.M.S.); (T.J.B.)
| | - Nina S. Ponomarenko
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, IL 60439, USA; (N.S.P.); (K.L.M.); (D.M.T.)
| | - Coleman M. Swaim
- Department of Chemistry and Biochemistry, James Madison University, Harrisonburg, VA 22807, USA; (D.R.M.); (A.M.M.); (C.M.S.); (T.J.B.)
| | - Tyler J. Brittain
- Department of Chemistry and Biochemistry, James Madison University, Harrisonburg, VA 22807, USA; (D.R.M.); (A.M.M.); (C.M.S.); (T.J.B.)
| | - Natalie L. Simmons
- Department of Biology, James Madison University, Harrisonburg, VA 22807, USA;
| | | | - Karen L. Mulfort
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, IL 60439, USA; (N.S.P.); (K.L.M.); (D.M.T.)
| | - David M. Tiede
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, IL 60439, USA; (N.S.P.); (K.L.M.); (D.M.T.)
| | - Oleksandr Kokhan
- Department of Chemistry and Biochemistry, James Madison University, Harrisonburg, VA 22807, USA; (D.R.M.); (A.M.M.); (C.M.S.); (T.J.B.)
| |
Collapse
|
34
|
Dong X, Mahler AD, Kempfer-Robertson EM, Thompson LM. Global Elucidation of Self-Consistent Field Solution Space Using Basin Hopping. J Chem Theory Comput 2020; 16:5635-5644. [DOI: 10.1021/acs.jctc.0c00488] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xinju Dong
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40205, United States
| | - Andrew D. Mahler
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40205, United States
| | | | - Lee M. Thompson
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40205, United States
| |
Collapse
|
35
|
Abstract
This work explores the reactivity of singlet oxygen with respect to two typical reactions: cycloaddition to anthracene and excitation energy transfer (EET) to a carotenoid using diabatic states with multistate density functional theory (MSDFT). Noticeably, the degenerate state 1Δg has distinct open-shell (OS) and closed-shell (CS) components, and the closed-shell component showed more reactivity than the open-shell one due to the strong diabatic couplings to the product diabatic states. The diabatic perspective presented in this work could also apply to general singlet fission processes.
Collapse
Affiliation(s)
- Zexing Qu
- Institute of Theoretical Chemistry, Jilin University, Changchun, 130023, China.
| |
Collapse
|
36
|
Kühne TD, Iannuzzi M, Del Ben M, Rybkin VV, Seewald P, Stein F, Laino T, Khaliullin RZ, Schütt O, Schiffmann F, Golze D, Wilhelm J, Chulkov S, Bani-Hashemian MH, Weber V, Borštnik U, Taillefumier M, Jakobovits AS, Lazzaro A, Pabst H, Müller T, Schade R, Guidon M, Andermatt S, Holmberg N, Schenter GK, Hehn A, Bussy A, Belleflamme F, Tabacchi G, Glöß A, Lass M, Bethune I, Mundy CJ, Plessl C, Watkins M, VandeVondele J, Krack M, Hutter J. CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations. J Chem Phys 2020; 152:194103. [PMID: 33687235 DOI: 10.1063/5.0007045] [Citation(s) in RCA: 924] [Impact Index Per Article: 231.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.
Collapse
Affiliation(s)
- Thomas D Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Mauro Del Ben
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Vladimir V Rybkin
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Patrick Seewald
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Frederick Stein
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Teodoro Laino
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | - Rustam Z Khaliullin
- Department of Chemistry, McGill University, CH-801 Sherbrooke St. West, Montreal, Quebec H3A 0B8, Canada
| | - Ole Schütt
- Department of Materials, ETH Zürich, CH-8092 Zürich, Switzerland
| | | | - Dorothea Golze
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| | - Jan Wilhelm
- Institute of Theoretical Physics, University of Regensburg, Universitätsstraße 31, D-93053 Regensburg, Germany
| | - Sergey Chulkov
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | | | - Valéry Weber
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | | | | | | | | | - Hans Pabst
- Intel Extreme Computing, Software and Systems, Zürich, Switzerland
| | - Tiziano Müller
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Robert Schade
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Manuel Guidon
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Samuel Andermatt
- Integrated Systems Laboratory, ETH Zürich, CH-8092 Zürich, Switzerland
| | - Nico Holmberg
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076 Aalto, Finland
| | - Gregory K Schenter
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Anna Hehn
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Augustin Bussy
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Fabian Belleflamme
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Gloria Tabacchi
- Department of Science and High Technology, University of Insubria and INSTM, via Valleggio 9, I-22100 Como, Italy
| | - Andreas Glöß
- BASF SE, Carl-Bosch-Straße 38, D-67056 Ludwigshafen am Rhein, Germany
| | - Michael Lass
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Iain Bethune
- Hartree Centre, Science and Technology Facilities Council, Sci-Tech Daresbury, Warrington WA4 4AD, United Kingdom
| | - Christopher J Mundy
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Christian Plessl
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Matt Watkins
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | - Joost VandeVondele
- Swiss National Supercomputing Centre (CSCS), ETH Zürich, Zürich, Switzerland
| | - Matthias Krack
- Laboratory for Scientific Computing and Modelling, Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
37
|
Prentice JCA, Aarons J, Womack JC, Allen AEA, Andrinopoulos L, Anton L, Bell RA, Bhandari A, Bramley GA, Charlton RJ, Clements RJ, Cole DJ, Constantinescu G, Corsetti F, Dubois SMM, Duff KKB, Escartín JM, Greco A, Hill Q, Lee LP, Linscott E, O'Regan DD, Phipps MJS, Ratcliff LE, Serrano ÁR, Tait EW, Teobaldi G, Vitale V, Yeung N, Zuehlsdorff TJ, Dziedzic J, Haynes PD, Hine NDM, Mostofi AA, Payne MC, Skylaris CK. The ONETEP linear-scaling density functional theory program. J Chem Phys 2020; 152:174111. [PMID: 32384832 DOI: 10.1063/5.0004445] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange-correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications.
Collapse
Affiliation(s)
- Joseph C A Prentice
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Jolyon Aarons
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - James C Womack
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Alice E A Allen
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Lampros Andrinopoulos
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Lucian Anton
- UKAEA, Culham Science Centre, Abingdon OX14 3DB, United Kingdom
| | - Robert A Bell
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Arihant Bhandari
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Gabriel A Bramley
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Robert J Charlton
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Rebecca J Clements
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Daniel J Cole
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| | - Gabriel Constantinescu
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Fabiano Corsetti
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Simon M-M Dubois
- Institute of Condensed Matter and Nanosciences, Université Catholique de Louvain, Louvain-la-Neuve, Belgium
| | - Kevin K B Duff
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - José María Escartín
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Andrea Greco
- Department of Physics, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Quintin Hill
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Louis P Lee
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Edward Linscott
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - David D O'Regan
- School of Physics, AMBER, and CRANN Institute, Trinity College Dublin, The University of Dublin, Dublin 2, Ireland
| | - Maximillian J S Phipps
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Laura E Ratcliff
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Álvaro Ruiz Serrano
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Edward W Tait
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Gilberto Teobaldi
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Valerio Vitale
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nelson Yeung
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Peter D Haynes
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Nicholas D M Hine
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Arash A Mostofi
- Department of Materials, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Mike C Payne
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
38
|
Spiegelman F, Tarrat N, Cuny J, Dontot L, Posenitskiy E, Martí C, Simon A, Rapacioli M. Density-functional tight-binding: basic concepts and applications to molecules and clusters. ADVANCES IN PHYSICS: X 2020; 5:1710252. [PMID: 33154977 PMCID: PMC7116320 DOI: 10.1080/23746149.2019.1710252] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2019] [Accepted: 12/19/2019] [Indexed: 06/10/2023] Open
Abstract
The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.
Collapse
Affiliation(s)
- Fernand Spiegelman
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Nathalie Tarrat
- CEMES, Université de Toulouse (UPS), CNRS, UPR8011, Toulouse, Toulouse, France
| | - Jérôme Cuny
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Leo Dontot
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Evgeny Posenitskiy
- Laboratoire Collisions Agrégats et Réactivité LCAR/IRSAMC, UMR5589, Université de Toulouse (UPS) and CNRS, Toulouse, France
| | - Carles Martí
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
- Laboratoire de Chimie, UMR5182, Ecole Normale Supérieure de Lyon, Université de Lyon and CNRS, Lyon, France
| | - Aude Simon
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Mathias Rapacioli
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| |
Collapse
|
39
|
Charge reduction in ions in the ionic liquid 1-ethy-2,3-dimethylimidazolium bis(trifluoromethanesulfonyl)imide on the Au(111) surface. Theor Chem Acc 2020. [DOI: 10.1007/s00214-019-2527-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
40
|
Multistate density functional theory applied with 3 unpaired electrons in 3 orbitals: The singdoublet and tripdoublet states of the ethylene cation. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.136803] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
41
|
Fornari RP, de Silva P. Unexpectedly Large Couplings Between Orthogonal Units in Anthraquinone Polymers. Chemistry 2019; 25:14651-14658. [PMID: 31509312 DOI: 10.1002/chem.201903227] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2019] [Indexed: 01/06/2023]
Abstract
The unusual electronic properties of directly linked 1,4-polyanthraquinones (14PAQs) are investigated. The dihedral angle between the molecular planes of anthraquinones (AQs) is found to be close to 90°. Contrary to the prevailing notion that the interaction between orthogonal units is negligible due to broken π-electron conjugation, the coupling between neighboring AQ units does not have a minimum at 90° and is much larger than that expected. The unexpectedly large electronic coupling between orthogonal AQ units is explained by the interaction between the lone pairs of the carbonyl oxygen and the π system of the neighboring unit, which allows favorable overlap between frontier molecular orbitals at the orthogonal geometry. It is shown that this effect, which is described computationally for the first time, can be strengthened by adding more quinone units. The effect of thermal fluctuations on the couplings is assessed through ab initio molecular dynamics simulations. The distributions of the couplings reveal that electron transport is resilient to dynamic disorder in all systems considered, whereas the hole couplings are much more sensitive to disorder. Lone pair-π interactions are described, as a previously largely overlooked conjugation mechanism, for incorporation into a new class of disorder-resilient semiconducting redox polymers.
Collapse
Affiliation(s)
- Rocco P Fornari
- Department of Energy Conversion and Storage, Technical University of Denmark, Fysikvej 309, 2800 Kongens, Lyngby, Denmark
| | - Piotr de Silva
- Department of Energy Conversion and Storage, Technical University of Denmark, Fysikvej 309, 2800 Kongens, Lyngby, Denmark
| |
Collapse
|
42
|
Dontot L, Spiegelman F, Rapacioli M. Structures and Energetics of Neutral and Cationic Pyrene Clusters. J Phys Chem A 2019; 123:9531-9543. [PMID: 31589446 PMCID: PMC6917508 DOI: 10.1021/acs.jpca.9b07007] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The low energy structures of neutral and cationic pyrene clusters containing up to seven molecules are searched through a global exploration scheme combining parallel tempering Monte Carlo algorithm and local quenches. The potential energies are computed at the density functional based tight binding level for neutrals and configuration interaction density functional based tight binding for cations in order to treat properly the charge resonance. New simplified versions of these schemes are also presented and used during the global exploration. Neutral clusters are shown to be made of compact assemblies of sub-blocs containing up to three units whereas cations present a charged dimer or trimer core surrounded by neutral units. The structural features of the clusters are analyzed and correlated for the cation with the charge distribution. The stability of clusters is also discussed in terms of cohesive and evaporation energies. Adiabatic and vertical ionization potentials are also discussed.
Collapse
Affiliation(s)
- Léo Dontot
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Fernand Spiegelman
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Mathias Rapacioli
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| |
Collapse
|
43
|
Mao Y, Montoya-Castillo A, Markland TE. Accurate and efficient DFT-based diabatization for hole and electron transfer using absolutely localized molecular orbitals. J Chem Phys 2019; 151:164114. [DOI: 10.1063/1.5125275] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Affiliation(s)
- Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | | | - Thomas E. Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
44
|
Mester D, Kállay M. Combined Density Functional and Algebraic-Diagrammatic Construction Approach for Accurate Excitation Energies and Transition Moments. J Chem Theory Comput 2019; 15:4440-4453. [PMID: 31265275 DOI: 10.1021/acs.jctc.9b00391] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
A composite of time-dependent density functional theory (TDDFT) and the second-order algebraic-diagrammatic construction [ADC(2)] approach is presented for efficient calculation of spectral properties of molecules. Our method can be regarded as a new excited-state double-hybrid (DH) approach or a dressed TDDFT scheme, but it can also be interpreted as an empirically tuned ADC(2) model. Several combinations of exchange-correlation functionals and spin-scaling schemes are explored. Our best-performing method includes the Perdew, Burke, and Ernzerhof exchange and Perdew's 1986 correlation functional and employs the scaled-opposite-spin approximation for the higher-order terms. The computation time of the new method scales as the fourth power of the system size, and an efficient cost-reduction approach is also presented, which further speeds up the calculations. Our benchmark calculations show that the proposed model outperforms not only the existing DH approaches and ADC(2) variants but also the considerably more expensive coupled-cluster methods.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science , Budapest University of Technology and Economics , P.O. Box 91, H-1521 Budapest , Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science , Budapest University of Technology and Economics , P.O. Box 91, H-1521 Budapest , Hungary
| |
Collapse
|
45
|
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
|
46
|
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
|
47
|
Alagna N, Han J, Wollscheid N, Perez Lustres JL, Herz J, Hahn S, Koser S, Paulus F, Bunz UHF, Dreuw A, Buckup T, Motzkus M. Tailoring Ultrafast Singlet Fission by the Chemical Modification of Phenazinothiadiazoles. J Am Chem Soc 2019; 141:8834-8845. [DOI: 10.1021/jacs.9b01079] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Nicolò Alagna
- Centre for Advanced Materials
, Im Neuenheimer Feld 225, D-69120 Heidelberg, Germany
| | | | - Nikolaus Wollscheid
- Centre for Advanced Materials
, Im Neuenheimer Feld 225, D-69120 Heidelberg, Germany
| | - J. Luis Perez Lustres
- Centre for Advanced Materials
, Im Neuenheimer Feld 225, D-69120 Heidelberg, Germany
| | | | | | | | | | - Uwe H. F. Bunz
- Centre for Advanced Materials
, Im Neuenheimer Feld 225, D-69120 Heidelberg, Germany
| | - Andreas Dreuw
- Centre for Advanced Materials
, Im Neuenheimer Feld 225, D-69120 Heidelberg, Germany
| | - Tiago Buckup
- Centre for Advanced Materials
, Im Neuenheimer Feld 225, D-69120 Heidelberg, Germany
| | - Marcus Motzkus
- Centre for Advanced Materials
, Im Neuenheimer Feld 225, D-69120 Heidelberg, Germany
| |
Collapse
|
48
|
Ying F, Zhou C, Zheng P, Luan J, Su P, Wu W. λ-Density Functional Valence Bond: A Valence Bond-Based Multiconfigurational Density Functional Theory With a Single Variable Hybrid Parameter. Front Chem 2019; 7:225. [PMID: 31041304 PMCID: PMC6476929 DOI: 10.3389/fchem.2019.00225] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2019] [Accepted: 03/22/2019] [Indexed: 11/13/2022] Open
Abstract
A new valence bond (VB)-based multireference density functional theory (MRDFT) method, named λ-DFVB, is presented in this paper. The method follows the idea of the hybrid multireference density functional method theory proposed by Sharkas et al. (2012). λ-DFVB combines the valence bond self-consistent field (VBSCF) method with Kohn–Sham density functional theory (KS-DFT) by decomposing the electron–electron interactions with a hybrid parameter λ. Different from the Toulouse's scheme, the hybrid parameter λ in λ-DFVB is variable, defined as a function of a multireference character of a molecular system. Furthermore, the EC correlation energy of a leading determinant is introduced to ensure size consistency at the dissociation limit. Satisfactory results of test calculations, including potential energy surfaces, bond dissociation energies, reaction barriers, and singlet–triplet energy gaps, show the potential capability of λ-DFVB for molecular systems with strong correlation.
Collapse
Affiliation(s)
- Fuming Ying
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen University, Xiamen, China.,The State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, China.,College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Chen Zhou
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen University, Xiamen, China.,The State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, China.,College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Peikun Zheng
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen University, Xiamen, China.,The State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, China.,College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Jiamin Luan
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen University, Xiamen, China.,The State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, China.,College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Peifeng Su
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen University, Xiamen, China.,The State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, China.,College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| | - Wei Wu
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen University, Xiamen, China.,The State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, China.,College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China
| |
Collapse
|
49
|
Glebov IO, Kozlov MI, Poddubnyy VV. Comparison of the Coulomb and non-orthogonal approaches to the construction of the exciton Hamiltonian. COMPUT THEOR CHEM 2019. [DOI: 10.1016/j.comptc.2019.02.010] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
50
|
Bai S, Barbatti M. Mechanism of Spin-Exchange Internal Conversion: Practical Proxies for Diabatic and Nonadiabatic Couplings. J Chem Theory Comput 2019; 15:1503-1513. [PMID: 30735372 DOI: 10.1021/acs.jctc.8b00923] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Spin-exchange internal conversion (SEIC) is a general class of reactions having singlet fission and triplet fusion as particular cases. Based on a charge transfer (CT) mediated mechanism and analytical derivation with a model Hamiltonian, we propose proxies for estimating the coupling strength in both diabatic and adiabatic pictures for general SEIC reactions. In the diabatic picture, we demonstrated the existence of a bilinear relationship between the coupling strength and molecular orbital overlap, which provides a practical way to predict diabatic couplings. In the adiabatic picture, we showed that nonadiabatic couplings can be approximated by simple functions of the wave function CT coefficients. These approaches were verified through the investigation of singlet oxygen photosensitization, where both 1Δg and 1Σg oxygen states can be competitively generated by a triplet fusion reaction. The interplay between the CT-mediated mechanism, the spatial factors of the bimolecular complex, and the electronic structure of the oxygen molecule during the reaction explains the curiously small coupling to the 1Σg state along specific incidence directions. The results from both the diabatic and adiabatic pictures provide a comprehensive understanding of the reaction mechanism, which applies to general SEIC problems.
Collapse
Affiliation(s)
- Shuming Bai
- Aix-Marseille Univ., CNRS, ICR , 13397 Marseille , France
| | - Mario Barbatti
- Aix-Marseille Univ., CNRS, ICR , 13397 Marseille , France
| |
Collapse
|