1
|
Wang Y, Havenridge S, Aikens CM. Implementation of energy and gradient for the TDDFT-approximate auxiliary function (aas) method. J Chem Phys 2024; 161:024101. [PMID: 38980092 DOI: 10.1063/5.0213587] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2024] [Accepted: 06/15/2024] [Indexed: 07/10/2024] Open
Abstract
In this work, we have implemented the time-dependent density functional theory approximate auxiliary s function (TDDFT-aas) method, which is an approximate TDDFT method. Instead of calculating the exact two-center electron integrals in the K coupling matrix when solving the Casida equation, we approximate the integrals, thereby reducing the computational cost. In contrast to the related TDDFT plus tight-binding (TDDFT+TB) method, a new type of gamma function is used in the coupling matrix that does not depend on the tight-binding parameters. The calculated absorption spectra of silver and gold nanoparticles using TDDFT-aas show good agreement with TDDFT and TDDFT+TB results. In addition, we have implemented the analytical excited-state gradients for the TDDFT-aas method, which makes it possible to calculate the emission energy of molecular systems.
Collapse
Affiliation(s)
- Yuchen Wang
- Department of Chemistry, Kansas State University, Manhattan, Kansas 66506, USA
| | - Shana Havenridge
- Department of Chemistry, Kansas State University, Manhattan, Kansas 66506, USA
| | - Christine M Aikens
- Department of Chemistry, Kansas State University, Manhattan, Kansas 66506, USA
| |
Collapse
|
2
|
Lee IS, Filatov M, Min SK. Formulation of transition dipole gradients for non-adiabatic dynamics with polaritonic states. J Chem Phys 2024; 160:154103. [PMID: 38624116 DOI: 10.1063/5.0202095] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Accepted: 03/31/2024] [Indexed: 04/17/2024] Open
Abstract
A general formulation of the strong coupling between photons confined in a cavity and molecular electronic states is developed for the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham method. The light-matter interaction is included in the Jaynes-Cummings model, which requires the derivation and implementation of the analytical derivatives of the transition dipole moments between the molecular electronic states. The developed formalism is tested in the simulations of the nonadiabatic dynamics in the polaritonic states resulting from the strong coupling between the cavity photon mode and the ground and excited states of the penta-2,4-dieniminium cation, also known as PSB3. Comparison with the field-free simulations of the excited-state decay dynamics in PSB3 reveals that the light-matter coupling can considerably alter the decay dynamics by increasing the excited state lifetime and hindering photochemically induced torsion about the C=C double bonds of PSB3. The necessity of obtaining analytical transition dipole gradients for the accurate propagation of the dynamics is underlined.
Collapse
Affiliation(s)
- In Seong Lee
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
| | - Michael Filatov
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
| | - Seung Kyu Min
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, Republic of Korea
| |
Collapse
|
3
|
Zhou Z, Della Sala F, Parker SM. Minimal Auxiliary Basis Set Approach for the Electronic Excitation Spectra of Organic Molecules. J Phys Chem Lett 2023; 14:1968-1976. [PMID: 36787711 DOI: 10.1021/acs.jpclett.2c03698] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
We report a minimal auxiliary basis model for time-dependent density functional theory (TDDFT) with hybrid density functionals that can accurately reproduce excitation energies and absorption spectra from TDDFT while reducing cost by about 2 orders of magnitude. Our method, dubbed TDDFT-ris, employs the resolution-of-the-identity technique with just one s-type auxiliary basis function per atom for the linear response operator, where the Gaussian exponents are parametrized across the periodic table using tabulated atomic radii with a single global scaling factor. By tuning on a small test set, we determine a single functional-independent scale factor that balances errors in excitation energies and absorption spectra. Benchmarked on organic molecules and compared to standard TDDFT, TDDFT-ris has an average energy error of only 0.06 eV and yields absorption spectra in close agreement with TDDFT. Thus, TDDFT-ris enables simulation of realistic absorption spectra in large molecules that would be inaccessible from standard TDDFT.
Collapse
Affiliation(s)
- Zehao Zhou
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, Ohio 44106, United States
| | - Fabio Della Sala
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, 73010 Arnesano (LE), Italy
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
| | - Shane M Parker
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, Ohio 44106, United States
| |
Collapse
|
4
|
Bertoni AI, Sánchez CG. Data-driven approach for benchmarking DFTB-approximate excited state methods. Phys Chem Chem Phys 2023; 25:3789-3798. [PMID: 36645084 DOI: 10.1039/d2cp04979a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
In this work we propose a chemically-informed data-driven approach to benchmark the approximate density-functional tight-binding (DFTB) excited state (ES) methods that are currently available within the DFTB+ suite. By taking advantage of the large volume of low-detail ES-data in the machine learning (ML) dataset, QM8, we were able to extract valuable insights regarding the limitations of the benchmarked methods in terms of the approximations made to the parent formalism, density-functional theory (DFT), while providing recommendations on how to overcome them. For this benchmark, we compared the first singlet-singlet vertical excitation energies (E1) predicted by the DFTB-approximate methods with predictions of less approximate methods from the reference ML-dataset. For the nearly 21800 organic molecules in the GDB-8 chemical space, we were able to identify clear trends in the E1 prediction error distributions, with respect to second-order approximate coupled cluster (CC2), showing a strong dependence on chemical identity.
Collapse
Affiliation(s)
- Andrés I Bertoni
- Instituto Interdisciplinario de Ciencias Básicas (ICB-CONICET), Universidad Nacional de Cuyo, Padre Jorge Contreras 1300, Mendoza 5502, Argentina.
| | - Cristián G Sánchez
- Instituto Interdisciplinario de Ciencias Básicas (ICB-CONICET), Universidad Nacional de Cuyo, Padre Jorge Contreras 1300, Mendoza 5502, Argentina.
| |
Collapse
|
5
|
Kim TI, Lee IS, Kim H, Min SK. Calculation of exciton couplings based on density functional tight-binding coupled to state-interaction state-averaged ensemble-referenced Kohn-Sham approach. J Chem Phys 2023; 158:044106. [PMID: 36725518 DOI: 10.1063/5.0132361] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
We introduce the combination of the density functional tight binding (DFTB) approach, including onsite correction (OC) and long-range corrected (LC) functional and the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method with extended active space involving four electrons and four orbitals [LC-OC-DFTB/SSR(4,4)], to investigate exciton couplings in multichromophoric systems, such as organic crystals and molecular aggregates. We employ the LC-OC-DFTB/SSR(4,4) method to calculate the excitonic coupling in anthracene and tetracene. As a result, the LC-OC-DFTB/SSR(4,4) method provides a reliable description of the locally excited (LE) state in a single chromophore and the excitonic couplings between chromophores with reasonable accuracy compared to the experiment and the conventional SSR(4,4) method. In addition, the thermal fluctuation of excitonic couplings from dynamic nuclear motion in an anthracene crystal with LC-OC-DFTB/SSR(4,4) shows a similar fluctuation of excitonic coupling and spectral density with those of first-principle calculations. We conclude that LC-OC-DFTB/SSR(4,4) is capable of providing reasonable features related to LE states, such as Frenkel exciton with efficient computational cost.
Collapse
Affiliation(s)
- Tae In Kim
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - In Seong Lee
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Hwon Kim
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Seung Kyu Min
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| |
Collapse
|
6
|
Wu X, Wen S, Song H, Frauenheim T, Tretiak S, Yam C, Zhang Y. Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method. J Chem Phys 2022; 157:084114. [DOI: 10.1063/5.0100339] [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
Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. The nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.
Collapse
Affiliation(s)
- Xiaoyan Wu
- Shenzhen JL Computational Science and Applied Research Institute, Shenzhen JL Computational Science and Applied Research Institute, China
| | | | - Huajing Song
- Los Alamos National Laboratory, United States of America
| | | | - Sergei Tretiak
- Theoretical Division, T-1, Los Alamos National Laboratory, United States of America
| | - ChiYung Yam
- Beijing Computational Science Research Center, Beijing Computational Science Research Center, China
| | - Yu Zhang
- Theoretical Division, Los Alamos National Laboratory, United States of America
| |
Collapse
|
7
|
Mukherjee S, Pinheiro M, Demoulin B, Barbatti M. Simulations of molecular photodynamics in long timescales. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200382. [PMID: 35341303 PMCID: PMC8958277 DOI: 10.1098/rsta.2020.0382] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Accepted: 08/12/2021] [Indexed: 05/04/2023]
Abstract
Nonadiabatic dynamics simulations in the long timescale (much longer than 10 ps) are the next challenge in computational photochemistry. This paper delimits the scope of what we expect from methods to run such simulations: they should work in full nuclear dimensionality, be general enough to tackle any type of molecule and not require unrealistic computational resources. We examine the main methodological challenges we should venture to advance the field, including the computational costs of the electronic structure calculations, stability of the integration methods, accuracy of the nonadiabatic dynamics algorithms and software optimization. Based on simulations designed to shed light on each of these issues, we show how machine learning may be a crucial element for long time-scale dynamics, either as a surrogate for electronic structure calculations or aiding the parameterization of model Hamiltonians. We show that conventional methods for integrating classical equations should be adequate to extended simulations up to 1 ns and that surface hopping agrees semiquantitatively with wave packet propagation in the weak-coupling regime. We also describe our optimization of the Newton-X program to reduce computational overheads in data processing and storage. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
Affiliation(s)
| | - Max Pinheiro
- Aix Marseille University, CNRS, ICR, Marseille, France
| | | | - Mario Barbatti
- Aix Marseille University, CNRS, ICR, Marseille, France
- Institut Universitaire de France, 75231 Paris, France
| |
Collapse
|
8
|
Lee IS, Min SK. Generalized Formulation of the Density Functional Tight Binding-Based Restricted Ensemble Kohn-Sham Method with Onsite Correction to Long-Range Correction. J Chem Theory Comput 2022; 18:3391-3409. [PMID: 35549266 DOI: 10.1021/acs.jctc.2c00037] [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
We present a generalized formulation for the combination of the density functional tight binding (DFTB) approach and the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method by considering onsite correction (OC) as well as the long-range corrected (LC) functional. The OC contribution provides more accurate energies and analytic gradients for individual microstates, while the multireference character of the SSR provides the correct description for conical intersections. We benchmark the LC-OC-DFTB/SSR method against various DFTB calculation methods for excitation energies and conical intersection structures with π/π* or n/π* characters. Furthermore, we perform excited-state molecular dynamics simulations with a molecular rotary motor with variations of LC-OC-DFTB/SSR approaches. We show that the OC contribution to the LC functional is crucial to obtain the correct geometry of conical intersections.
Collapse
Affiliation(s)
- In Seong Lee
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, South Korea
| | - Seung Kyu Min
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, South Korea
| |
Collapse
|
9
|
Sarkar R, Loos PF, Boggio-Pasqua M, Jacquemin D. Assessing the Performances of CASPT2 and NEVPT2 for Vertical Excitation Energies. J Chem Theory Comput 2022; 18:2418-2436. [PMID: 35333060 DOI: 10.1021/acs.jctc.1c01197] [Citation(s) in RCA: 25] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events but also to provide reference values for vertical transition energies, hence allowing benchmarking of lower-order models. In this category, both the complete-active-space second-order perturbation theory (CASPT2) and the N-electron valence state second-order perturbation theory (NEVPT2) are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches, the present work relies on highly accurate (±0.03 eV) aug-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 n → π*, 119 π → π*, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially contracted (PC) and strongly contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of 0.11 and 0.13 eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focused on single-reference wave function methods ( J. Chem. Theory Comput. 2018, 14, 4360; J. Chem. Theory Comput. 2020, 16, 1711), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character.
Collapse
Affiliation(s)
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, CNRS, UPS, Université de Toulouse, Toulouse 31062, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques, CNRS, UPS, Université de Toulouse, Toulouse 31062, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
10
|
Hutama AS, Marlina LA, Chou CP, Irle S, Hofer TS. Development of Density-Functional Tight-Binding Parameters for the Molecular Dynamics Simulation of Zirconia, Yttria, and Yttria-Stabilized Zirconia. ACS OMEGA 2021; 6:20530-20548. [PMID: 34395999 PMCID: PMC8359130 DOI: 10.1021/acsomega.1c02411] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/17/2021] [Accepted: 07/15/2021] [Indexed: 06/13/2023]
Abstract
In this work, a set of density-functional tight-binding (DFTB) parameters for the Zr-Zr, Zr-O, Y-Y, Y-O, and Zr-Y interactions was developed for bulk and surface simulations of ZrO2 (zirconia), Y2O3 (yttria), and yttria-stabilized zirconia (YSZ) materials. The parameterization lays the ground work for realistic simulations of zirconia-, yttria-, and YSZ-based electrolytes in solid oxide fuel cells and YSZ-based catalysts on long timescales and relevant size scales. The parameterization was validated for the zirconia and yttria polymorphs observed under standard conditions based on density functional theory calculations and experimental data. Additionally, we performed DFTB-based molecular dynamics (MD) simulations to compute structural and vibrational properties of these materials. The results show that the parameters can give a qualitatively correct phase ordering of zirconia, where the tetragonal phase is more stable than the cubic phase at a lower temperature. The lattice parameters are only slightly overestimated by 0.05-0.1 Å (2% error), still within the typical accuracy of first-principles methods. Additionally, the MD results confirm that zirconia and yttria phases are stable against transformations under standard conditions. The parameterization also predicts that vibrational spectra are within the range of 100-1000 cm-1 for zirconia and 100-800 cm-1 for yttria, which is in good agreement with predictions both from full quantum mechanics and a recently developed classical force field. To further demonstrate the advantage of the developed DFTB parameters in terms of computational resources, we conducted DFTB/MD simulations of the YSZ4 and YS12 models containing approximately 750 atoms.
Collapse
Affiliation(s)
- Aulia Sukma Hutama
- Department
of Chemistry, Faculty of Mathematics and Natural Sciences, Universitas Gadjah Mada, Yogyakarta 55281, Indonesia
| | - Lala Adetia Marlina
- Department
of Chemistry, Faculty of Mathematics and Natural Sciences, Universitas Gadjah Mada, Yogyakarta 55281, Indonesia
| | - Chien-Pin Chou
- Department
of Applied Chemistry, National Chiao Tung
University, Hsinchu 30010, Taiwan
| | - Stephan Irle
- Computational
Sciences and Engineering Division & Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, United States
| | - Thomas S. Hofer
- Theoretical
Chemistry Division, Institute of General, Inorganic and Theoretical
Chemistry, Center for Chemistry and Biomedicine, University of Innsbruck, Innsbruck A-6020, Austria
| |
Collapse
|
11
|
Lee IS, Ha JK, Han D, Kim TI, Moon SW, Min SK. PyUNIxMD: A Python-based excited state molecular dynamics package. J Comput Chem 2021; 42:1755-1766. [PMID: 34197646 PMCID: PMC8362049 DOI: 10.1002/jcc.26711] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 06/10/2021] [Accepted: 06/16/2021] [Indexed: 01/17/2023]
Abstract
Theoretical/computational description of excited state molecular dynamics is nowadays a crucial tool for understanding light-matter interactions in many materials. Here we present an open-source Python-based nonadiabatic molecular dynamics program package, namely PyUNIxMD, to deal with mixed quantum-classical dynamics for correlated electron-nuclear propagation. The PyUNIxMD provides many interfaces for quantum chemical calculation methods with commercial and noncommercial ab initio and semiempirical quantum chemistry programs. In addition, the PyUNIxMD offers many nonadiabatic molecular dynamics algorithms such as fewest-switch surface hopping and its derivatives as well as decoherence-induced surface hopping based on the exact factorization (DISH-XF) and coupled-trajectory mixed quantum-classical dynamics (CTMQC) for general purposes. Detailed structures and flows of PyUNIxMD are explained for the further implementations by developers. We perform a nonadiabatic molecular dynamics simulation for a molecular motor system as a simple demonstration.
Collapse
Affiliation(s)
- In Seong Lee
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Jong-Kwon Ha
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Daeho Han
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Tae In Kim
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Sung Wook Moon
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Seung Kyu Min
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| |
Collapse
|
12
|
de Wergifosse M, Grimme S. Perspective on Simplified Quantum Chemistry Methods for Excited States and Response Properties. J Phys Chem A 2021; 125:3841-3851. [PMID: 33928774 DOI: 10.1021/acs.jpca.1c02362] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
We review recent developments in the framework of simplified quantum chemistry for excited state and optical response properties (sTD-DFT) and present future challenges for new method developments to improve accuracy and extend the range of application. In recent years, the scope of sTD-DFT was extended to molecular response calculations of the polarizability, optical rotation, first hyperpolarizability, two-photon absorption (2PA), and excited-state absorption for large systems with hundreds to thousands of atoms. The recently introduced spin-flip simplified time-dependent density functional theory (SF-sTD-DFT) variant enables an ultrafast treatment for diradicals and related strongly correlated systems. A few drawbacks were also identified, specifically for the computation of 2PA cross sections. We propose solutions to this problem and how to generally improve the accuracy of simplified schemes. New possible simplified schemes are also introduced for strongly correlated systems, e.g., with a second-order perturbative correlation correction. Interpretation tools that can extract chemical structure-property relationships from excited state or response calculations are also discussed. In particular, the recently introduced method-agnostic RespA approach based on natural response orbitals (NROs) as the key concept is employed.
Collapse
Affiliation(s)
- Marc de Wergifosse
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstrasse 4, D-53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstrasse 4, D-53115 Bonn, Germany
| |
Collapse
|
13
|
Darghouth AAMHM, Casida ME, Zhu X, Natarajan B, Su H, Humeniuk A, Titov E, Miao X, Mitrić R. Effect of varying the TD-lc-DFTB range-separation parameter on charge and energy transfer in a model pentacene/buckminsterfullerene heterojunction. J Chem Phys 2021; 154:054102. [PMID: 33557554 DOI: 10.1063/5.0024559] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Atomistic modeling of energy and charge transfer at the heterojunction of organic solar cells is an active field with many remaining outstanding questions owing, in part, to the difficulties in performing reliable photodynamics calculations on very large systems. One approach to being able to overcome these difficulties is to design and apply an appropriate simplified method. Density-functional tight binding (DFTB) has become a popular form of approximate density-functional theory based on a minimal valence basis set and neglect of all but two center integrals. We report the results of our tests of a recent long-range correction (lc) [A. Humeniuk and R. Mitrić, J. Chem. Phys. 143, 134120 (2015)] for time-dependent (TD) lc-DFTB by carrying out TD-lc-DFTB fewest switches surface hopping calculations of energy and charge transfer times using the relatively new DFTBABY [A. Humeniuk and R. Mitrić, Comput. Phys. Commun. 221, 174 (2017)] program. An advantage of this method is the ability to run enough trajectories to get meaningful ensemble averages. Our interest in the present work is less in determining exact energy and charge transfer rates than in understanding how the results of these calculations vary with the value of the range-separation parameter (Rlc = 1/μ) for a model organic solar cell heterojunction consisting of a gas-phase van der Waals complex P/F made up of a single pentacene (P) molecule together with a single buckminsterfullerene (F) molecule. The default value of Rlc = 3.03 a0 is found to be much too small as neither energy nor charge transfer is observed until Rlc ≈ 10 a0. Tests at a single geometry show that the best agreement with high-quality ab initio spectra is obtained in the limit of no lc (i.e., very large Rlc). A plot of energy and charge transfer rates as a function of Rlc is provided, which suggests that a value of Rlc ≈ 15 a0 yields the typical literature (condensed-phase) charge transfer time of about 100 fs. However, energy and charge transfer times become as high as ∼300 fs for Rlc ≈ 25 a0. A closer examination of the charge transfer process P*/F → P+/F- shows that the initial electron transfer is accompanied by a partial delocalization of the P hole onto F, which then relocalizes back onto P, consistent with a polaron-like picture in which the nuclei relax to stabilize the resultant redistribution of charges.
Collapse
Affiliation(s)
| | - Mark E Casida
- Laboratoire de Spectrométrie, Interactions et Chimie Théorique (SITh), Département de Chimie Moléculaire (DCM), Institut de Chimie Moléculaire de Grenoble (ICMG), Université Grenoble-Alpes, 301 rue de la Chimie, CS 40700, 38058 Grenoble Cedex 9, France
| | - Xi Zhu
- Institute of Advanced Studies, Nanyang Technological University, 60 Nanyang View, 639673, Singapore
| | - Bhaarathi Natarajan
- Institute of Advanced Studies, Nanyang Technological University, 60 Nanyang View, 639673, Singapore
| | - Haibin Su
- Institute of Advanced Studies, Nanyang Technological University, 60 Nanyang View, 639673, Singapore
| | - Alexander Humeniuk
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| | - Evgenii Titov
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| | - Xincheng Miao
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| | - Roland Mitrić
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| |
Collapse
|
14
|
Wagstaffe M, Wenthaus L, Dominguez-Castro A, Chung S, Lana Semione GD, Palutke S, Mercurio G, Dziarzhytski S, Redlin H, Klemke N, Yang Y, Frauenheim T, Dominguez A, Kärtner F, Rubio A, Wurth W, Stierle A, Noei H. Ultrafast Real-Time Dynamics of CO Oxidation over an Oxide Photocatalyst. ACS Catal 2020. [DOI: 10.1021/acscatal.0c04098] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Lukas Wenthaus
- Deutsches Elektronen-Synchrotron, Hamburg D-22607, Germany
- Center for Free-Electron Laser Science, Hamburg D-22761, Germany
| | | | - Simon Chung
- Deutsches Elektronen-Synchrotron, Hamburg D-22607, Germany
- Fachbereich Physik Universität Hamburg, Hamburg D-20355, Germany
| | | | | | | | | | - Harald Redlin
- Deutsches Elektronen-Synchrotron, Hamburg D-22607, Germany
| | - Nicolai Klemke
- Center for Free-Electron Laser Science, Hamburg D-22761, Germany
| | - Yudong Yang
- Center for Free-Electron Laser Science, Hamburg D-22761, Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Material Science (BCCMS), Bremen D-28359, Germany
- Computational Science and Applied Research Institute (CSAR), Shenzhen 518110, China
- Beijing Computational Science Research Center (CSRC), Beijing 100193, China
| | - Adriel Dominguez
- Bremen Center for Computational Material Science (BCCMS), Bremen D-28359, Germany
- Computational Science and Applied Research Institute (CSAR), Shenzhen 518110, China
- Beijing Computational Science Research Center (CSRC), Beijing 100193, China
- Nano-Bio Spectroscopy Group, Departamento de Fisica de Materiales, Universidad del País Vasco, UPV/EHU, San Sebastián 20018, Spain
| | - Franz Kärtner
- Center for Free-Electron Laser Science, Hamburg D-22761, Germany
| | - Angel Rubio
- Center for Free-Electron Laser Science, Hamburg D-22761, Germany
- Nano-Bio Spectroscopy Group, Departamento de Fisica de Materiales, Universidad del País Vasco, UPV/EHU, San Sebastián 20018, Spain
- Max Planck Institute for the Structure and Dynamics of Matter, Hamburg D-22761, Germany
- Center for Computational Quantum Physics, Flatiron Institute, New York 10010, United States
| | - Wilfried Wurth
- Deutsches Elektronen-Synchrotron, Hamburg D-22607, Germany
- Center for Free-Electron Laser Science, Hamburg D-22761, Germany
- Fachbereich Physik Universität Hamburg, Hamburg D-20355, Germany
| | - Andreas Stierle
- Deutsches Elektronen-Synchrotron, Hamburg D-22607, Germany
- Fachbereich Physik Universität Hamburg, Hamburg D-20355, Germany
| | - Heshmat Noei
- Deutsches Elektronen-Synchrotron, Hamburg D-22607, Germany
| |
Collapse
|
15
|
Maity S, Bold BM, Prajapati JD, Sokolov M, Kubař T, Elstner M, Kleinekathöfer U. DFTB/MM Molecular Dynamics Simulations of the FMO Light-Harvesting Complex. J Phys Chem Lett 2020; 11:8660-8667. [PMID: 32991176 DOI: 10.1021/acs.jpclett.0c02526] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Because of the size of light-harvesting complexes and the involvement of electronic degrees of freedom, computationally these systems need to be treated with a combined quantum-classical description. To this end, Born-Oppenheimer molecular dynamics simulations have been employed in a quantum mechanics/molecular mechanics (QM/MM) fashion for the ground state followed by excitation energy calculations again in a QM/MM scheme for the Fenna-Matthews-Olson (FMO) complex. The self-consistent-charge density functional tight-binding (DFTB) method electrostatically coupled to a classical description of the environment was applied to perform the ground-state dynamics. Subsequently, long-range-corrected time-dependent DFTB calculations were performed to determine the excitation energy fluctuations of the individual bacteriochlorophyll a molecules. The spectral densities obtained using this approach show an excellent agreement with experimental findings. In addition, the fluctuating site energies and couplings were used to estimate the exciton transfer dynamics.
Collapse
Affiliation(s)
- Sayan Maity
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Beatrix M Bold
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | | | - Monja Sokolov
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Tomáš Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
- Institute of Biological Interfaces, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
16
|
Giannone G, Della Sala F. Minimal auxiliary basis set for time-dependent density functional theory and comparison with tight-binding approximations: Application to silver nanoparticles. J Chem Phys 2020; 153:084110. [DOI: 10.1063/5.0020545] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Giulia Giannone
- Center for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
- Department of Mathematics and Physics “E. De Giorgi,” University of Salento, Via Arnesano, Lecce, Italy
| | - Fabio Della Sala
- Center for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
| |
Collapse
|
17
|
Bonafé FP, Aradi B, Hourahine B, Medrano CR, Hernández FJ, Frauenheim T, Sánchez CG. A Real-Time Time-Dependent Density Functional Tight-Binding Implementation for Semiclassical Excited State Electron–Nuclear Dynamics and Pump–Probe Spectroscopy Simulations. J Chem Theory Comput 2020; 16:4454-4469. [DOI: 10.1021/acs.jctc.9b01217] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Franco P. Bonafé
- Max Planck Institute for the Structure and Dynamics of Matter, Hamburg, Germany
- Universidad Nacional de Córdoba, Facultad de Ciencias Químicas, Departamento de Química Teórica y Computacional, Córdoba, Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba, INFIQC (CONICET - Universidad Nacional de Córdoba), Córdoba, Argentina
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, Universitát Bremen, Bremen, Germany
| | - Ben Hourahine
- SUPA, Department of Physics, John Anderson Building, The University of Strathclyde, 107 Rottenrow, Glasgow G15 6QN, United Kingdom
| | - Carlos R. Medrano
- Universidad Nacional de Córdoba, Facultad de Ciencias Químicas, Departamento de Química Teórica y Computacional, Córdoba, Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba, INFIQC (CONICET - Universidad Nacional de Córdoba), Córdoba, Argentina
| | - Federico J. Hernández
- Universidad Nacional de Córdoba, Facultad de Ciencias Químicas, Departamento de Química Teórica y Computacional, Córdoba, Argentina
- Instituto de Investigaciones en Fisicoquímica de Córdoba, INFIQC (CONICET - Universidad Nacional de Córdoba), Córdoba, Argentina
- Department of Physics, Universidad de Santiago de Chile, Av. Ecuador 3493, Santiago, Chile
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, Universitát Bremen, Bremen, Germany
- Computational Science Research Center (CSRC) Beijing and Computational Science and Applied Research (CSAR) Institute, Shenzhen, China
| | - Cristián G. Sánchez
- Instituto Interdisciplinario de Ciencias Básicas, Universidad Nacional de Cuyo, CONICET, Facultad de Ciencias Exactas y Naturales, Mendoza, Argentina
| |
Collapse
|
18
|
Hourahine B, Aradi B, Blum V, Bonafé F, Buccheri A, Camacho C, Cevallos C, Deshaye MY, Dumitrică T, Dominguez A, Ehlert S, Elstner M, van der Heide T, Hermann J, Irle S, Kranz JJ, Köhler C, Kowalczyk T, Kubař T, Lee IS, Lutsker V, Maurer RJ, Min SK, Mitchell I, Negre C, Niehaus TA, Niklasson AMN, Page AJ, Pecchia A, Penazzi G, Persson MP, Řezáč J, Sánchez CG, Sternberg M, Stöhr M, Stuckenberg F, Tkatchenko A, Yu VWZ, Frauenheim T. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations. J Chem Phys 2020; 152:124101. [PMID: 32241125 DOI: 10.1063/1.5143190] [Citation(s) in RCA: 409] [Impact Index Per Article: 102.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green's functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.
Collapse
Affiliation(s)
- B Hourahine
- SUPA, Department of Physics, The University of Strathclyde, Glasgow G4 0NG, United Kingdom
| | - B Aradi
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - V Blum
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - F Bonafé
- Max Planck Institute for the Structure and Dynamics of Matter, Hamburg, Germany
| | - A Buccheri
- School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - C Camacho
- School of Chemistry, University of Costa Rica, San José 11501-2060, Costa Rica
| | - C Cevallos
- School of Chemistry, University of Costa Rica, San José 11501-2060, Costa Rica
| | - M Y Deshaye
- Department of Chemistry and Advanced Materials Science and Engineering Center, Western Washington University, Bellingham, Washington 98225, USA
| | - T Dumitrică
- Department of Mechanical Engineering, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - A Dominguez
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - S Ehlert
- University of Bonn, Bonn, Germany
| | - M Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - T van der Heide
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - J Hermann
- Freie Universität Berlin, Berlin, Germany
| | - S Irle
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - J J Kranz
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - C Köhler
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - T Kowalczyk
- Department of Chemistry and Advanced Materials Science and Engineering Center, Western Washington University, Bellingham, Washington 98225, USA
| | - T Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - I S Lee
- Department of Chemistry, Ulsan National Institute of Science and Technology, Ulsan, South Korea
| | - V Lutsker
- Institut I - Theoretische Physik, University of Regensburg, Regensburg, Germany
| | - R J Maurer
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - S K Min
- Department of Chemistry, Ulsan National Institute of Science and Technology, Ulsan, South Korea
| | - I Mitchell
- Center for Multidimensional Carbon Materials, Institute for Basic Science (IBS), Ulsan 44919, South Korea
| | - C Negre
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - T A Niehaus
- Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - A M N Niklasson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - A J Page
- School of Environmental and Life Sciences, University of Newcastle, Callaghan, Australia
| | - A Pecchia
- CNR-ISMN, Via Salaria km 29.300, 00015 Monterotondo Stazione, Rome, Italy
| | - G Penazzi
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - M P Persson
- Dassault Systemes, Cambridge, United Kingdom
| | - J Řezáč
- Institute of Organic Chemistry and Biochemistry AS CR, Prague, Czech Republic
| | - C G Sánchez
- Instituto Interdisciplinario de Ciencias Básicas, Universidad Nacional de Cuyo, CONICET, Facultad de Ciencias Exactas y Naturales, Mendoza, Argentina
| | - M Sternberg
- Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - M Stöhr
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - F Stuckenberg
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - A Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - V W-Z Yu
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - T Frauenheim
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| |
Collapse
|
19
|
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
|
20
|
Loos PF, Lipparini F, Boggio-Pasqua M, Scemama A, Jacquemin D. A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Medium Sized Molecules. J Chem Theory Comput 2020; 16:1711-1741. [PMID: 31986042 DOI: 10.1021/acs.jctc.9b01216] [Citation(s) in RCA: 107] [Impact Index Per Article: 26.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Following our previous work focusing on compounds containing up to 3 non-hydrogen atoms [J. Chem. Theory Comput. 2018, 14, 4360-4379], we present here highly accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms: acetone, acrolein, benzene, butadiene, cyanoacetylene, cyanoformaldehyde, cyanogen, cyclopentadiene, cyclopropenone, cyclopropenethione, diacetylene, furan, glyoxal, imidazole, isobutene, methylenecyclopropene, propynal, pyrazine, pyridazine, pyridine, pyrimidine, pyrrole, tetrazine, thioacetone, thiophene, thiopropynal, and triazine. To obtain these energies, we use equation-of-motion/linear-response coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ) and selected configuration interaction (SCI) calculations (with tens of millions of determinants in the reference space), as well as the multiconfigurational n-electron valence state perturbation theory (NEVPT2) method. All these approaches are applied in combination with diffuse-containing atomic basis sets. For all transitions, we report at least CC3/aug-cc-pVQZ vertical excitation energies as well as CC3/aug-cc-pVTZ oscillator strengths for each dipole-allowed transition. We show that CC3 almost systematically delivers transition energies in agreement with higher-level methods with a typical deviation of ±0.04 eV, except for transitions with a dominant double excitation character where the error is much larger. The present contribution gathers a large, diverse, and accurate set of more than 200 highly accurate transition energies for states of various natures (valence, Rydberg, singlet, triplet, n → π*, π → π*, ...). We use this series of theoretical best estimates to benchmark a series of popular methods for excited state calculations: CIS(D), ADC(2), CC2, STEOM-CCSD, EOM-CCSD, CCSDR(3), CCSDT-3, CC3, and NEVPT2. The results of these benchmarks are compared to the available literature data.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Denis Jacquemin
- CEISAM Lab, UMR 6230, Université de Nantes, CNRS, F-44000 Nantes, France
| |
Collapse
|
21
|
Fihey A, Jacquemin D. Performances of Density Functional Tight-Binding Methods for Describing Ground and Excited State Geometries of Organic Molecules. J Chem Theory Comput 2019; 15:6267-6276. [DOI: 10.1021/acs.jctc.9b00688] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Arnaud Fihey
- Institut des Sciences Chimiques de Rennes, UMR 6226 CNRS, Université de Rennes 1, 263 Av. du Général Leclerc, 35042 Cedex Rennes, France
| | - Denis Jacquemin
- Laboratoire CEISAM - UMR CNRS 6230, Université de Nantes, 2 Rue de la Houssinière, BP 92208, 44322 Cedex 3 Nantes, France
| |
Collapse
|
22
|
Hanasaki K, Kanno M, Niehaus TA, Kono H. An efficient approximate algorithm for nonadiabatic molecular dynamics. J Chem Phys 2019; 149:244117. [PMID: 30599729 DOI: 10.1063/1.5046757] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
We propose a modification to the nonadiabatic surface hopping calculation method formulated in a paper by Yu et al. [Phys. Chem. Chem. Phys. 16, 25883 (2014)], which is a multidimensional extension of the Zhu-Nakamura theory with a practical diabatic gradient estimation algorithm. In our modification, their diabatic gradient estimation algorithm, which is based on a simple interpolation of the adiabatic potential energy surfaces, is replaced by an algorithm using the numerical derivatives of the adiabatic gradients. We then apply the algorithm to several models of nonadiabatic dynamics, both analytic and ab initio models, to numerically demonstrate that our method indeed widens the applicability and robustness of their method. We also discuss the validity and limitations of our new nonadiabatic surface hopping method while considering in mind potential applications to excited-state dynamics of biomolecules or unconventional nonadiabatic dynamics such as radiation decay processes in ultraintense X-ray fields.
Collapse
Affiliation(s)
- Kota Hanasaki
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| | - Manabu Kanno
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| | - Thomas A Niehaus
- University of Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeubanne, France
| | - Hirohiko Kono
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| |
Collapse
|
23
|
Tran T, Prlj A, Lin KH, Hollas D, Corminboeuf C. Mechanisms of fluorescence quenching in prototypical aggregation-induced emission systems: excited state dynamics with TD-DFTB. Phys Chem Chem Phys 2019; 21:9026-9035. [DOI: 10.1039/c9cp00691e] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
A recent implementation of time-dependent tight-binding density functional theory is employed in excited state molecular dynamics for the investigation of the fluorescence quenching mechanism in 3 prototypical aggregation-induced emission systems.
Collapse
Affiliation(s)
- Thierry Tran
- Laboratory for Computational Molecular Design
- Institute of Chemical Sciences and Engineering
- Ecole Polytechnique Federale de Lausanne (EPFL)
- CH-1015 Lausanne
- Switzerland
| | - Antonio Prlj
- Laboratory for Computational Molecular Design
- Institute of Chemical Sciences and Engineering
- Ecole Polytechnique Federale de Lausanne (EPFL)
- CH-1015 Lausanne
- Switzerland
| | - Kun-Han Lin
- Laboratory for Computational Molecular Design
- Institute of Chemical Sciences and Engineering
- Ecole Polytechnique Federale de Lausanne (EPFL)
- CH-1015 Lausanne
- Switzerland
| | - Daniel Hollas
- Laboratory for Computational Molecular Design
- Institute of Chemical Sciences and Engineering
- Ecole Polytechnique Federale de Lausanne (EPFL)
- CH-1015 Lausanne
- Switzerland
| | - Clémence Corminboeuf
- Laboratory for Computational Molecular Design
- Institute of Chemical Sciences and Engineering
- Ecole Polytechnique Federale de Lausanne (EPFL)
- CH-1015 Lausanne
- Switzerland
| |
Collapse
|
24
|
Darghouth AAMHM, Correa GC, Juillard S, Casida ME, Humeniuk A, Mitrić R. Davydov-type excitonic effects on the absorption spectra of parallel-stacked and herringbone aggregates of pentacene: Time-dependent density-functional theory and time-dependent density-functional tight binding. J Chem Phys 2018; 149:134111. [PMID: 30292200 DOI: 10.1063/1.5025624] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
Exciton formation leads to J-bands in solid pentacene. Describing these exciton bands represents a challenge for both time-dependent (TD) density-functional theory (DFT) and for its semi-empirical analog, namely, for TD density-functional tight binding (DFTB) for three reasons: (i) solid pentacene and pentacene aggregates are bound only by van der Waals forces which are notoriously difficult to describe with DFT and DFTB, (ii) the proper description of the long-range coupling between molecules, needed to describe Davydov splitting, is not easy to include in TD-DFT with traditional functionals and in TD-DFTB, and (iii) mixing may occur between local and charge transfer excitons, which may, in turn, require special functionals. We assess how far TD-DFTB has progressed toward a correct description of this type of exciton by including both a dispersion correction for the ground state and a range-separated hybrid functional for the excited state and comparing the results against corresponding TD-CAM-B3LYP/CAM-B3LYP+D3 results. Analytic results for parallel-stacked ethylene are derived which go beyond Kasha's exciton model [M. Kasha, H. R. Rawls, and A. El-Bayoumi, Pure Appl. Chem. 11, 371 (1965)] in that we are able to make a clear distinction between charge transfer and energy transfer excitons. This is further confirmed when it is shown that range-separated hybrids have a markedly greater effect on charge-transfer excitons than on energy-transfer excitons in the case of parallel-stacked pentacenes. TD-DFT calculations with the CAM-B3LYP functional and TD-lc-DFT calculations lead to negligible excitonic corrections for the herringbone crystal structure, possibly because of an overcorrection of charge-transfer effects (CAM refers to Coulomb attenuated method). In this case, TD-DFT calculations with the B3LYP functional or TD-DFTB calculations parameterized to B3LYP give the best results for excitonic corrections for the herringbone crystal structure as judged from comparison with experimental spectra and with Bethe-Salpeter equation calculations from the literature.
Collapse
Affiliation(s)
- Ala Aldin M H M Darghouth
- Laboratoire de Chimie Théorique, Département de Chimie Moléculaire (DCM), Institut de Chimie Moléculaire de Grenoble (ICMG), Université Grenoble-Alpes, 301 Rue de la Chimie, CS 40700, 38058 Grenoble Cedex 9, France
| | - Gabriela Calinao Correa
- Laboratoire de Chimie Théorique, Département de Chimie Moléculaire (DCM), Institut de Chimie Moléculaire de Grenoble (ICMG), Université Grenoble-Alpes, 301 Rue de la Chimie, CS 40700, 38058 Grenoble Cedex 9, France
| | - Sacha Juillard
- Laboratoire de Chimie Théorique, Département de Chimie Moléculaire (DCM), Institut de Chimie Moléculaire de Grenoble (ICMG), Université Grenoble-Alpes, 301 Rue de la Chimie, CS 40700, 38058 Grenoble Cedex 9, France
| | - Mark E Casida
- Laboratoire de Chimie Théorique, Département de Chimie Moléculaire (DCM), Institut de Chimie Moléculaire de Grenoble (ICMG), Université Grenoble-Alpes, 301 Rue de la Chimie, CS 40700, 38058 Grenoble Cedex 9, France
| | - Alexander Humeniuk
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| | - Roland Mitrić
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| |
Collapse
|
25
|
Loos PF, Scemama A, Blondel A, Garniron Y, Caffarel M, Jacquemin D. A Mountaineering Strategy to Excited States: Highly Accurate Reference Energies and Benchmarks. J Chem Theory Comput 2018; 14:4360-4379. [DOI: 10.1021/acs.jctc.8b00406] [Citation(s) in RCA: 152] [Impact Index Per Article: 25.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31013 Toulouse Cedex 6, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31013 Toulouse Cedex 6, France
| | - Aymeric Blondel
- Laboratoire CEISAM - UMR CNRS 6230, Université de Nantes, 2 Rue de la Houssinière, BP 92208, 44322 Nantes Cedex 3, France
| | - Yann Garniron
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31013 Toulouse Cedex 6, France
| | - Michel Caffarel
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31013 Toulouse Cedex 6, France
| | - Denis Jacquemin
- Laboratoire CEISAM - UMR CNRS 6230, Université de Nantes, 2 Rue de la Houssinière, BP 92208, 44322 Nantes Cedex 3, France
| |
Collapse
|
26
|
Crespo-Otero R, Barbatti M. Recent Advances and Perspectives on Nonadiabatic Mixed Quantum–Classical Dynamics. Chem Rev 2018; 118:7026-7068. [DOI: 10.1021/acs.chemrev.7b00577] [Citation(s) in RCA: 301] [Impact Index Per Article: 50.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rachel Crespo-Otero
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | | |
Collapse
|
27
|
Ullah N, Chen S, Zhang R. Excited state dynamics study of the self-trapped exciton formation in silicon nanosheets. Phys Chem Chem Phys 2018; 20:29299-29305. [DOI: 10.1039/c8cp04806a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
After excitation to S1 (1), the exciton takes ∼450–850 femtoseconds to relax into the self-trapped (ST) state (2) with the occurrence of strong localization and a large Stokes shift, due to the significant stretching of the Si–Si bonds.
Collapse
Affiliation(s)
- Naeem Ullah
- Department of Physics, City University of Hong Kong
- Hong Kong SAR
- China
| | - Shunwei Chen
- Department of Physics, City University of Hong Kong
- Hong Kong SAR
- China
| | - Ruiqin Zhang
- Department of Physics, City University of Hong Kong
- Hong Kong SAR
- China
- Beijing Computational Science Research Center
- Beijing 100193
| |
Collapse
|
28
|
Stojanović L, Aziz SG, Hilal RH, Plasser F, Niehaus TA, Barbatti M. Nonadiabatic Dynamics of Cycloparaphenylenes with TD-DFTB Surface Hopping. J Chem Theory Comput 2017; 13:5846-5860. [DOI: 10.1021/acs.jctc.7b01000] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - Saadullah G. Aziz
- Chemistry
Department, Faculty of Science, King Abdulaziz University, Jeddah B.O.
208203, Saudi Arabia
| | - Rifaat H. Hilal
- Chemistry
Department, Faculty of Science, King Abdulaziz University, Jeddah B.O.
208203, Saudi Arabia
- Chemistry
Department, Faculty of Science, Cairo University, Giza, Egypt
| | - Felix Plasser
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Thomas A. Niehaus
- Univ Lyon, Université
Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Lyon, France
| | | |
Collapse
|
29
|
Yang Y, Dominguez A, Zhang D, Lutsker V, Niehaus TA, Frauenheim T, Yang W. Charge transfer excitations from particle-particle random phase approximation—Opportunities and challenges arising from two-electron deficient systems. J Chem Phys 2017; 146:124104. [DOI: 10.1063/1.4977928] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yang Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Adriel Dominguez
- Max Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
| | - Du Zhang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Vitalij Lutsker
- Department of Theoretical Physics, University of Regensburg, 93040 Regensburg, Germany
| | - Thomas A. Niehaus
- Université Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, Universität Bremen, Am Fallturm, 28359 Bremen, Germany
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
- Department of Physics, Duke University, Durham, North Carolina 27708, USA
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry and Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
30
|
Kranz JJ, Elstner M, Aradi B, Frauenheim T, Lutsker V, Garcia AD, Niehaus TA. Time-Dependent Extension of the Long-Range Corrected Density Functional Based Tight-Binding Method. J Chem Theory Comput 2017; 13:1737-1747. [PMID: 28272887 DOI: 10.1021/acs.jctc.6b01243] [Citation(s) in RCA: 60] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We present a consistent linear response formulation of the density functional based tight-binding method for long-range corrected exchange-correlation functionals (LC-DFTB). Besides a detailed account of derivation and implementation of the method, we also test the new scheme on a variety of systems considered to be problematic for conventional local/semilocal time-dependent density functional theory (TD-DFT). To this class belong the optical properties of polyacenes and nucleobases, as well as charge transfer excited states in molecular dimers. We find that the approximate LC-DFTB method exhibits the same general trends and similar accuracy as range-separated DFT methods at significantly reduced computational cost. The scheme should be especially useful in the determination of the electronic excited states of very large molecules, for which conventional TD-DFT is supposed to fail due to a multitude of artificial low energy states.
Collapse
Affiliation(s)
- Julian J Kranz
- Institute of Physical Chemistry and Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology , Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry and Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology , Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Bálint Aradi
- BCCMS, University of Bremen , 28359 Bremen, Germany
| | | | - Vitalij Lutsker
- Department of Theoretical Physics, University of Regensburg , 93040 Regensburg, Germany
| | - Adriel Dominguez Garcia
- Max Planck Institute for the Structure and Dynamics of Matter and Center for Free-Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Thomas A Niehaus
- Université Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Lyon, France
| |
Collapse
|
31
|
Rüger R, van Lenthe E, Heine T, Visscher L. Tight-binding approximations to time-dependent density functional theory - A fast approach for the calculation of electronically excited states. J Chem Phys 2017; 144:184103. [PMID: 27179467 DOI: 10.1063/1.4948647] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
We propose a new method of calculating electronically excited states that combines a density functional theory based ground state calculation with a linear response treatment that employs approximations used in the time-dependent density functional based tight binding (TD-DFTB) approach. The new method termed time-dependent density functional theory TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive TD-DFT calculations. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB.
Collapse
Affiliation(s)
- Robert Rüger
- Scientific Computing & Modelling NV, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Erik van Lenthe
- Scientific Computing & Modelling NV, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Thomas Heine
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Linnéstr. 2, 04103 Leipzig, Germany
| | - Lucas Visscher
- Department of Theoretical Chemistry, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
32
|
Rüger R, Niehaus T, van Lenthe E, Heine T, Visscher L. Vibrationally resolved UV/Vis spectroscopy with time-dependent density functional based tight binding. J Chem Phys 2016; 145:184102. [DOI: 10.1063/1.4966918] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Robert Rüger
- Scientific Computing and Modelling NV, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
- Department of Theoretical Chemistry, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Linnéstr. 2, 04103 Leipzig, Germany
| | - Thomas Niehaus
- Univ Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Erik van Lenthe
- Scientific Computing and Modelling NV, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Thomas Heine
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Linnéstr. 2, 04103 Leipzig, Germany
| | - Lucas Visscher
- Department of Theoretical Chemistry, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
33
|
Lutsker V, Aradi B, Niehaus TA. Implementation and benchmark of a long-range corrected functional in the density functional based tight-binding method. J Chem Phys 2016; 143:184107. [PMID: 26567646 DOI: 10.1063/1.4935095] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Bridging the gap between first principles methods and empirical schemes, the density functional based tight-binding method (DFTB) has become a versatile tool in predictive atomistic simulations over the past years. One of the major restrictions of this method is the limitation to local or gradient corrected exchange-correlation functionals. This excludes the important class of hybrid or long-range corrected functionals, which are advantageous in thermochemistry, as well as in the computation of vibrational, photoelectron, and optical spectra. The present work provides a detailed account of the implementation of DFTB for a long-range corrected functional in generalized Kohn-Sham theory. We apply the method to a set of organic molecules and compare ionization potentials and electron affinities with the original DFTB method and higher level theory. The new scheme cures the significant overpolarization in electric fields found for local DFTB, which parallels the functional dependence in first principles density functional theory (DFT). At the same time, the computational savings with respect to full DFT calculations are not compromised as evidenced by numerical benchmark data.
Collapse
Affiliation(s)
- V Lutsker
- Department of Theoretical Physics, University of Regensburg, 93040 Regensburg, Germany
| | - B Aradi
- BCCMS, University of Bremen, 28359 Bremen, Germany
| | - T A Niehaus
- Department of Theoretical Physics, University of Regensburg, 93040 Regensburg, Germany
| |
Collapse
|
34
|
Plötz PA, Polyutov SP, Ivanov SD, Fennel F, Wolter S, Niehaus T, Xie Z, Lochbrunner S, Würthner F, Kühn O. Biphasic aggregation of a perylene bisimide dye identified by exciton-vibrational spectra. Phys Chem Chem Phys 2016; 18:25110-25119. [DOI: 10.1039/c6cp04898f] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The quantum efficiency of light emission supramolecular aggregates strongly depends on the intermolecular coupling. We present a molecule which demonstrates two different aggregated structures with high and low quantum efficiency. The spectral signatures can be understood by simulating the aggregated structures and the corresponding exciton-vibrational spectra.
Collapse
Affiliation(s)
- P.-A. Plötz
- Institut für Physik
- Universität Rostock
- 18059 Rostock
- Germany
| | - S. P. Polyutov
- Institut für Physik
- Universität Rostock
- 18059 Rostock
- Germany
- Laboratory for Nonlinear Optics and Spectroscopy
| | - S. D. Ivanov
- Institut für Physik
- Universität Rostock
- 18059 Rostock
- Germany
| | - F. Fennel
- Institut für Physik
- Universität Rostock
- 18059 Rostock
- Germany
| | - S. Wolter
- Institut für Physik
- Universität Rostock
- 18059 Rostock
- Germany
| | - T. Niehaus
- Université Claude Bernard Lyon 1
- CNRS
- Institut Lumière Matière
- F-69622
- France
| | - Z. Xie
- Institut für Organische Chemie & Center for Nanosystems Chemistry
- Universität Würzburg
- 97074 Würzburg
- Germany
| | - S. Lochbrunner
- Institut für Physik
- Universität Rostock
- 18059 Rostock
- Germany
| | - F. Würthner
- Institut für Organische Chemie & Center for Nanosystems Chemistry
- Universität Würzburg
- 97074 Würzburg
- Germany
| | - O. Kühn
- Institut für Physik
- Universität Rostock
- 18059 Rostock
- Germany
| |
Collapse
|
35
|
Kowalczyk T, Le K, Irle S. Self-Consistent Optimization of Excited States within Density-Functional Tight-Binding. J Chem Theory Comput 2015; 12:313-23. [DOI: 10.1021/acs.jctc.5b00734] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Tim Kowalczyk
- Institute
of Transformative Bio-Molecules (WPI-ITbM) and Department of Chemistry,
Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
- Department
of Chemistry, Advanced Materials Science and Engineering Center, and
Institute for Energy Studies, Western Washington University, Bellingham, Washington 98225, United States
| | - Khoa Le
- Department
of Chemistry, Advanced Materials Science and Engineering Center, and
Institute for Energy Studies, Western Washington University, Bellingham, Washington 98225, United States
| | - Stephan Irle
- Institute
of Transformative Bio-Molecules (WPI-ITbM) and Department of Chemistry,
Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
| |
Collapse
|
36
|
Humeniuk A, Mitrić R. Long-range correction for tight-binding TD-DFT. J Chem Phys 2015; 143:134120. [DOI: 10.1063/1.4931179] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Alexander Humeniuk
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians Universität Würzburg, Emil-Fischer-Straße 42, 97074 Würzburg, Germany
| | - Roland Mitrić
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians Universität Würzburg, Emil-Fischer-Straße 42, 97074 Würzburg, Germany
| |
Collapse
|
37
|
Time-dependent density-functional tight-binding method with the third-order expansion of electron density. J Chem Phys 2015; 143:094108. [DOI: 10.1063/1.4929926] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
38
|
Gaus M, Jin H, Demapan D, Christensen AS, Goyal P, Elstner M, Cui Q. DFTB3 Parametrization for Copper: The Importance of Orbital Angular Momentum Dependence of Hubbard Parameters. J Chem Theory Comput 2015; 11:4205-19. [PMID: 26575916 PMCID: PMC4827604 DOI: 10.1021/acs.jctc.5b00600] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We report the parametrization of a density functional tight binding method (DFTB3) for copper in a spin-polarized formulation. The parametrization is consistent with the framework of 3OB for main group elements (ONCHPS) and can be readily used for biological applications that involve copper proteins/peptides. The key to our parametrization is to introduce orbital angular momentum dependence of the Hubbard parameter and its charge derivative, thus allowing the 3d and 4s orbitals to adopt different sizes and responses to the change of charge state. The parametrization has been tested by applying to a fairly broad set of molecules of biological relevance, and the properties of interest include optimized geometries, ligand binding energies, and ligand proton affinities. Compared to the reference QM level (B3LYP/aug-cc-pVTZ, which is shown here to be similar to the B97-1 and CCSD(T) results, in terms of many properties of interest for a set of small copper containing molecules), our parametrization generally gives reliable structural properties for both Cu(I) and Cu(II) compounds, although several exceptions are also noted. For energetics, the results are more accurate for neutral ligands than for charged ligands, likely reflecting the minimal basis limitation of DFTB3; the results generally outperform NDDO based methods such as PM6 and even PBE with the 6-31+G(d,p) basis. For all ligand types, single-point B3LYP calculations at DFTB3 geometries give results very close (∼1-2 kcal/mol) to the reference B3LYP values, highlighting the consistency between DFTB3 and B3LYP structures. Possible further developments of the DFTB3 model for a better treatment of transition-metal ions are also discussed. In the current form, our first generation of DFTB3 copper model is expected to be particularly valuable as a method that drives sampling in systems that feature a dynamical copper binding site.
Collapse
Affiliation(s)
- Michael Gaus
- Department of Chemistry, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Haiyun Jin
- Department of Chemistry, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Darren Demapan
- Department of Chemistry, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Anders S Christensen
- Department of Chemistry, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Puja Goyal
- Department of Chemistry, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology , Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - Qiang Cui
- Department of Chemistry, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
39
|
Fihey A, Hettich C, Touzeau J, Maurel F, Perrier A, Köhler C, Aradi B, Frauenheim T. SCC-DFTB parameters for simulating hybrid gold-thiolates compounds. J Comput Chem 2015; 36:2075-87. [PMID: 26280464 DOI: 10.1002/jcc.24046] [Citation(s) in RCA: 60] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2015] [Revised: 07/09/2015] [Accepted: 07/10/2015] [Indexed: 01/11/2023]
Abstract
We present a parametrization of a self-consistent charge density functional-based tight-binding scheme (SCC-DFTB) to describe gold-organic hybrid systems by adding new Au-X (X = Au, H, C, S, N, O) parameters to a previous set designed for organic molecules. With the aim of describing gold-thiolates systems within the DFTB framework, the resulting parameters are successively compared with density functional theory (DFT) data for the description of Au bulk, Aun gold clusters (n = 2, 4, 8, 20), and Aun SCH3 (n = 3 and 25) molecular-sized models. The geometrical, energetic, and electronic parameters obtained at the SCC-DFTB level for the small Au3 SCH3 gold-thiolate compound compare very well with DFT results, and prove that the different binding situations of the sulfur atom on gold are correctly described with the current parameters. For a larger gold-thiolate model, Au25 SCH3 , the electronic density of states and the potential energy surfaces resulting from the chemisorption of the molecule on the gold aggregate obtained with the new SCC-DFTB parameters are also in good agreement with DFT results.
Collapse
Affiliation(s)
- Arnaud Fihey
- Laboratoire Interfaces, Traitements, Organisation Et Dynamique Des Systèmes (ITODYS), CNRS UMR 7086, Université Paris Diderot Sorbonne Paris Cité, Bâtiment Lavoisier, 15 Rue Jean Antoine De Baïf, Paris Cedex 13, 75205, France
| | - Christian Hettich
- Bremen Center for Computational Materials Science (BCCMS), Universität Bremen, Am Fallturm 1, Bremen, 28359, Germany
| | - Jérémy Touzeau
- Laboratoire Interfaces, Traitements, Organisation Et Dynamique Des Systèmes (ITODYS), CNRS UMR 7086, Université Paris Diderot Sorbonne Paris Cité, Bâtiment Lavoisier, 15 Rue Jean Antoine De Baïf, Paris Cedex 13, 75205, France
| | - François Maurel
- Laboratoire Interfaces, Traitements, Organisation Et Dynamique Des Systèmes (ITODYS), CNRS UMR 7086, Université Paris Diderot Sorbonne Paris Cité, Bâtiment Lavoisier, 15 Rue Jean Antoine De Baïf, Paris Cedex 13, 75205, France
| | - Aurélie Perrier
- Laboratoire Interfaces, Traitements, Organisation Et Dynamique Des Systèmes (ITODYS), CNRS UMR 7086, Université Paris Diderot Sorbonne Paris Cité, Bâtiment Lavoisier, 15 Rue Jean Antoine De Baïf, Paris Cedex 13, 75205, France
| | - Christof Köhler
- Bremen Center for Computational Materials Science (BCCMS), Universität Bremen, Am Fallturm 1, Bremen, 28359, Germany
| | - Bálint Aradi
- Bremen Center for Computational Materials Science (BCCMS), Universität Bremen, Am Fallturm 1, Bremen, 28359, Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science (BCCMS), Universität Bremen, Am Fallturm 1, Bremen, 28359, Germany
| |
Collapse
|
40
|
Nurbawono A, Liu S, Zhang C. Modeling optical properties of silicon clusters by first principles: From a few atoms to large nanocrystals. J Chem Phys 2015; 142:154705. [DOI: 10.1063/1.4918588] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Argo Nurbawono
- Department of Physics and the Centre for Advanced 2D Materials, National University of Singapore, 2 Science Drive 3, Singapore
| | - Shuanglong Liu
- Department of Physics and the Centre for Advanced 2D Materials, National University of Singapore, 2 Science Drive 3, Singapore
| | - Chun Zhang
- Department of Physics and the Centre for Advanced 2D Materials, National University of Singapore, 2 Science Drive 3, Singapore
- Department of Chemistry, National University of Singapore, 3 Science Drive 3, Singapore
| |
Collapse
|
41
|
Domínguez A, Niehaus TA, Frauenheim T. Accurate hydrogen bond energies within the density functional tight binding method. J Phys Chem A 2015; 119:3535-44. [PMID: 25763597 DOI: 10.1021/acs.jpca.5b01732] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The density-functional-based tight-binding (DFTB) approach has been recently extended by incorporating one-center exchange-like terms in the expansion of the multicenter integrals. This goes beyond the Mulliken approximation and leads to a scheme which treats in a self-consistent way the fluctuations of the whole dual density matrix and not only its diagonal elements (Mulliken charges). To date, only the performance of this new formalism to reproduce excited-state properties has been assessed (Domínguez et al. J. Chem. Theory Comput., 2013, 9, 4901-4914). Here we study the effect of our corrections on the computation of hydrogen bond energies for water clusters and water-containing systems. The limitations of traditional DFTB to reproduce hydrogen bonds has been acknowledged often. We compare our results for a set of 22 small water clusters and water-containing systems as well as for five water hexadecamers to those obtained with the DFTB3 method. Additionally, we combine our extension with a third-order energy expansion in the charge fluctuations. Our results show that the new formalisms significantly improve upon original DFTB.
Collapse
Affiliation(s)
- A Domínguez
- †Bremen Center for Computational Materials Science, Universität Bremen, Am Fallturm 1, 28359 Bremen, Germany
| | - T A Niehaus
- ‡Department of Theoretical Physics, University of Regensburg, 93040 Regensburg, Germany
| | - T Frauenheim
- †Bremen Center for Computational Materials Science, Universität Bremen, Am Fallturm 1, 28359 Bremen, Germany
| |
Collapse
|
42
|
Niehaus TA, Hofbeck T, Yersin H. Charge-transfer excited states in phosphorescent organo-transition metal compounds: a difficult case for time dependent density functional theory? RSC Adv 2015. [DOI: 10.1039/c5ra12962a] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
A series of 17 platinum(ii) and iridium(iii) complexes have been investigated theoretically and experimentally to elucidate the charge-transfer character in emission from the lowest triplet state. TDDFT is found to be surprisingly accurate.
Collapse
Affiliation(s)
- Thomas A. Niehaus
- Institut für Theoretische Physik
- Universität Regensburg
- 93040 Regensburg
- Germany
| | - Thomas Hofbeck
- Institut für Physikalische Chemie
- Universität Regensburg
- 93040 Regensburg
- Germany
| | - Hartmut Yersin
- Institut für Physikalische Chemie
- Universität Regensburg
- 93040 Regensburg
- Germany
| |
Collapse
|
43
|
Rüger R, van Lenthe E, Lu Y, Frenzel J, Heine T, Visscher L. Efficient Calculation of Electronic Absorption Spectra by Means of Intensity-Selected Time-Dependent Density Functional Tight Binding. J Chem Theory Comput 2014; 11:157-67. [DOI: 10.1021/ct500838h] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Robert Rüger
- Scientific Computing & Modelling NV, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
- Department
of Theoretical Chemistry, VU University Amsterdam, De Boelelaan
1083, 1081 HV Amsterdam, The Netherlands
| | - Erik van Lenthe
- Scientific Computing & Modelling NV, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - You Lu
- Scientific Computing & Modelling NV, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Johannes Frenzel
- Department
of Chemistry, University of Calgary, 2500 University Drive, N.W., T2N 1N4 Calgary, Canada
| | - Thomas Heine
- School
of Engineering and Science, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Lucas Visscher
- Department
of Theoretical Chemistry, VU University Amsterdam, De Boelelaan
1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
44
|
Negre CFA, Young KJ, Oviedo MB, Allen LJ, Sánchez CG, Jarzembska KN, Benedict JB, Crabtree RH, Coppens P, Brudvig GW, Batista VS. Photoelectrochemical Hole Injection Revealed in Polyoxotitanate Nanocrystals Functionalized with Organic Adsorbates. J Am Chem Soc 2014; 136:16420-9. [DOI: 10.1021/ja509270f] [Citation(s) in RCA: 63] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Christian F. A. Negre
- Department
of Chemistry and Energy Sciences Institute, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States
| | - Karin J. Young
- Department
of Chemistry and Energy Sciences Institute, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States
| | - Ma. Belén Oviedo
- Department
of Chemistry, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | - Laura J. Allen
- Department
of Chemistry and Energy Sciences Institute, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States
| | - Cristián G. Sánchez
- Departamento
de Matemática y Física, Facultad de Ciencias Químicas,
INFIQC, Universidad Nacional de Córdoba, Ciudad Universitaria, X5000HUA, Córdoba, Argentina
| | - Katarzyna N. Jarzembska
- Department
of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
- Department
of Chemistry, University of Warsaw, Pasteura 1, 02-093, Warszawa, Poland
| | - Jason B. Benedict
- Department
of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Robert H. Crabtree
- Department
of Chemistry and Energy Sciences Institute, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States
| | - Philip Coppens
- Department
of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Gary W. Brudvig
- Department
of Chemistry and Energy Sciences Institute, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States
| | - Victor S. Batista
- Department
of Chemistry and Energy Sciences Institute, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States
| |
Collapse
|
45
|
Plötz PA, Niehaus T, Kühn O. A new efficient method for calculation of Frenkel exciton parameters in molecular aggregates. J Chem Phys 2014; 140:174101. [DOI: 10.1063/1.4871658] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
|