1
|
Lee IS, Filatov M, Min SK. Formulation of transition dipole gradients for non-adiabatic dynamics with polaritonic states. J Chem Phys 2024; 160:154103. [PMID: 38624116 DOI: 10.1063/5.0202095] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Accepted: 03/31/2024] [Indexed: 04/17/2024] Open
Abstract
A general formulation of the strong coupling between photons confined in a cavity and molecular electronic states is developed for the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham method. The light-matter interaction is included in the Jaynes-Cummings model, which requires the derivation and implementation of the analytical derivatives of the transition dipole moments between the molecular electronic states. The developed formalism is tested in the simulations of the nonadiabatic dynamics in the polaritonic states resulting from the strong coupling between the cavity photon mode and the ground and excited states of the penta-2,4-dieniminium cation, also known as PSB3. Comparison with the field-free simulations of the excited-state decay dynamics in PSB3 reveals that the light-matter coupling can considerably alter the decay dynamics by increasing the excited state lifetime and hindering photochemically induced torsion about the C=C double bonds of PSB3. The necessity of obtaining analytical transition dipole gradients for the accurate propagation of the dynamics is underlined.
Collapse
Affiliation(s)
- In Seong Lee
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
| | - Michael Filatov
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
| | - Seung Kyu Min
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, Republic of Korea
| |
Collapse
|
2
|
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
|
3
|
Effect of Orientational Isomerism in Neutral Water Hexamers on Their Thermodynamic Properties and Concentrations in the Gas Phase. J CLUST SCI 2022. [DOI: 10.1007/s10876-022-02365-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
4
|
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
|
5
|
Theoretical study of the formation of C18H and C18H2 molecules by low energy irradiation with atomic and molecular hydrogen. Radiat Phys Chem Oxf Engl 1993 2021. [DOI: 10.1016/j.radphyschem.2020.109166] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
6
|
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
|
7
|
Carlson RK, Cawkwell MJ, Batista ER, Yang P. Tight-Binding Modeling of Uranium in an Aqueous Environment. J Chem Theory Comput 2020; 16:3073-3083. [DOI: 10.1021/acs.jctc.0c00089] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rebecca K. Carlson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - M. J. Cawkwell
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Enrique R. Batista
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Ping Yang
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
8
|
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: 387] [Impact Index Per Article: 96.8] [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
|
9
|
Lodola A, Callegari D, Scalvini L, Rivara S, Mor M. Design and SAR Analysis of Covalent Inhibitors Driven by Hybrid QM/MM Simulations. Methods Mol Biol 2020; 2114:307-337. [PMID: 32016901 DOI: 10.1007/978-1-0716-0282-9_19] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Quantum mechanics/molecular mechanics (QM/MM) hybrid technique is emerging as a reliable computational method to investigate and characterize chemical reactions occurring in enzymes. From a drug discovery perspective, a thorough understanding of enzyme catalysis appears pivotal to assist the design of inhibitors able to covalently bind one of the residues belonging to the enzyme catalytic machinery. Thanks to the current advances in computer power, and the availability of more efficient algorithms for QM-based simulations, the use of QM/MM methodology is becoming a viable option in the field of covalent inhibitor design. In the present review, we summarized our experience in the field of QM/MM simulations applied to drug design problems which involved the optimization of agents working on two well-known drug targets, namely fatty acid amide hydrolase (FAAH) and epidermal growth factor receptor (EGFR). In this context, QM/MM simulations gave valuable information in terms of geometry (i.e., of transition states and metastable intermediates) and reaction energetics that allowed to correctly predict inhibitor binding orientation and substituent effect on enzyme inhibition. What is more, enzyme reaction modelling with QM/MM provided insights that were translated into the synthesis of new covalent inhibitor featured by a unique combination of intrinsic reactivity, on-target activity, and selectivity.
Collapse
Affiliation(s)
- Alessio Lodola
- Drug Design and Discovery Group, Department of Food and Drug, University of Parma, Parma, Italy.
| | - Donatella Callegari
- Drug Design and Discovery Group, Department of Food and Drug, University of Parma, Parma, Italy
| | - Laura Scalvini
- Drug Design and Discovery Group, Department of Food and Drug, University of Parma, Parma, Italy
| | - Silvia Rivara
- Drug Design and Discovery Group, Department of Food and Drug, University of Parma, Parma, Italy
| | - Marco Mor
- Drug Design and Discovery Group, Department of Food and Drug, University of Parma, Parma, Italy
| |
Collapse
|
10
|
Sakti AW, Nishimura Y, Nakai H. Recent advances in quantum‐mechanical molecular dynamics simulations of proton transfer mechanism in various water‐based environments. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1419] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Aditya W. Sakti
- Element Strategy Initiative for Catalysts and Batteries (ESICB) Kyoto University Kyoto Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE) Waseda University Tokyo Japan
| | - Hiromi Nakai
- Element Strategy Initiative for Catalysts and Batteries (ESICB) Kyoto University Kyoto Japan
- Waseda Research Institute for Science and Engineering (WISE) Waseda University Tokyo Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering Waseda University Tokyo Japan
| |
Collapse
|
11
|
Lee KH, Schnupf U, Sumpter BG, Irle S. Performance of Density-Functional Tight-Binding in Comparison to Ab Initio and First-Principles Methods for Isomer Geometries and Energies of Glucose Epimers in Vacuo and Solution. ACS OMEGA 2018; 3:16899-16915. [PMID: 31458314 PMCID: PMC6643604 DOI: 10.1021/acsomega.8b02213] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/29/2018] [Accepted: 11/07/2018] [Indexed: 05/12/2023]
Abstract
Density functional theory (DFT) is a widely used methodology for the computation of molecular and electronic structure, and we confirm that B3LYP and the high-level ab initio G3B3 method are in excellent agreement for the lowest-energy isomers of the 16 glucose epimers. Density-functional tight-binding (DFTB) is an approximate version of DFT with typically comparable accuracy that is 2 to 3 orders of magnitude faster, therefore generally very suitable for processing large numbers of complex structures. Conformational isomerism in sugars is well known to give rise to a large number of isomer structures. On the basis of a comprehensive study of glucose epimers in vacuo and aqueous solution, we found that the performance of DFTB is on par to B3LYP in terms of geometrical parameters excluding hydrogen bonds and isomer energies. However, DFTB underestimates both hydrogen bonding interactions as well as torsional barriers associated with rotations of the hydroxy groups, resulting in a counterintuitive overemphasis of hydrogen bonding in both gas phase as well as in water. Although the associated root mean squared deviation from B3LYP within epimer isomer groups is only on the order of 1 kcal/mol, this deviation affects the correct assignment of major isomer ordering, which span less than 10 kcal/mol. Both second- as well as third-order DFTB methods are exhibiting similar deviations from B3LYP. Even after the inclusion of empirical dispersion corrections in vacuum, these deviations remain for a large majority of isomer energies and geometries when compared to dispersion-corrected B3LYP.
Collapse
Affiliation(s)
- Ka Hung Lee
- Department
of Chemistry, Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
- Bredesen
Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Udo Schnupf
- Department
of Chemistry and Biochemistry, Bradley University, Peoria, Illinois 61625, United States
- E-mail: (U.S.)
| | - Bobby G. Sumpter
- Computational Sciences and Engineering Division &
Chemical Sciences
Division and Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Stephan Irle
- Department
of Chemistry, Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
- Computational Sciences and Engineering Division &
Chemical Sciences
Division and Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
- E-mail: (S.I.)
| |
Collapse
|
12
|
Simon A, Rapacioli M, Michoulier E, Zheng L, Korchagina K, Cuny J. Contribution of the density-functional-based tight-binding scheme to the description of water clusters: methods, applications and extension to bulk systems. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1554903] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- A. Simon
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| | - M. Rapacioli
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| | - E. Michoulier
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
- Laboratoire Collisions Agrégats et Réactivité LCAR/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| | - L. Zheng
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| | - K. Korchagina
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| | - J. Cuny
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| |
Collapse
|
13
|
Sakti AW, Nishimura Y, Chou CP, Nakai H. Density-Functional Tight-Binding Molecular Dynamics Simulations of Excess Proton Diffusion in Ice Ih, Ice Ic, Ice III, and Melted Ice VI Phases. J Phys Chem A 2017; 122:33-40. [DOI: 10.1021/acs.jpca.7b10664] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
| | | | | | - Hiromi Nakai
- CREST, Japan Science and Technology Agency, Tokyo 102-0075, Japan
- ESICB, Kyoto University, Kyoto 615-8520, Japan
| |
Collapse
|
14
|
Giese TJ, York DM. Quantum mechanical force fields for condensed phase molecular simulations. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:383002. [PMID: 28817382 PMCID: PMC5821073 DOI: 10.1088/1361-648x/aa7c5c] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Molecular simulations are powerful tools for providing atomic-level details into complex chemical and physical processes that occur in the condensed phase. For strongly interacting systems where quantum many-body effects are known to play an important role, density-functional methods are often used to provide the model with the potential energy used to drive dynamics. These methods, however, suffer from two major drawbacks. First, they are often too computationally intensive to practically apply to large systems over long time scales, limiting their scope of application. Second, there remain challenges for these models to obtain the necessary level of accuracy for weak non-bonded interactions to obtain quantitative accuracy for a wide range of condensed phase properties. Quantum mechanical force fields (QMFFs) provide a potential solution to both of these limitations. In this review, we address recent advances in the development of QMFFs for condensed phase simulations. In particular, we examine the development of QMFF models using both approximate and ab initio density-functional models, the treatment of short-ranged non-bonded and long-ranged electrostatic interactions, and stability issues in molecular dynamics calculations. Example calculations are provided for crystalline systems, liquid water, and ionic liquids. We conclude with a perspective for emerging challenges and future research directions.
Collapse
|
15
|
Řezáč J. Empirical Self-Consistent Correction for the Description of Hydrogen Bonds in DFTB3. J Chem Theory Comput 2017; 13:4804-4817. [DOI: 10.1021/acs.jctc.7b00629] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Jan Řezáč
- Institute of Organic Chemistry
and Biochemistry, Czech Academy of Sciences, 166 10 Prague, Czech Republic
| |
Collapse
|
16
|
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
|
17
|
Christensen A, Kubař T, Cui Q, Elstner M. Semiempirical Quantum Mechanical Methods for Noncovalent Interactions for Chemical and Biochemical Applications. Chem Rev 2016; 116:5301-37. [PMID: 27074247 PMCID: PMC4867870 DOI: 10.1021/acs.chemrev.5b00584] [Citation(s) in RCA: 243] [Impact Index Per Article: 30.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2015] [Indexed: 12/28/2022]
Abstract
Semiempirical (SE) methods can be derived from either Hartree-Fock or density functional theory by applying systematic approximations, leading to efficient computational schemes that are several orders of magnitude faster than ab initio calculations. Such numerical efficiency, in combination with modern computational facilities and linear scaling algorithms, allows application of SE methods to very large molecular systems with extensive conformational sampling. To reliably model the structure, dynamics, and reactivity of biological and other soft matter systems, however, good accuracy for the description of noncovalent interactions is required. In this review, we analyze popular SE approaches in terms of their ability to model noncovalent interactions, especially in the context of describing biomolecules, water solution, and organic materials. We discuss the most significant errors and proposed correction schemes, and we review their performance using standard test sets of molecular systems for quantum chemical methods and several recent applications. The general goal is to highlight both the value and limitations of SE methods and stimulate further developments that allow them to effectively complement ab initio methods in the analysis of complex molecular systems.
Collapse
Affiliation(s)
- Anders
S. Christensen
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Tomáš Kubař
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Qiang Cui
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Marcus Elstner
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| |
Collapse
|
18
|
Christensen AS, Elstner M, Cui Q. Improving intermolecular interactions in DFTB3 using extended polarization from chemical-potential equalization. J Chem Phys 2015; 143:084123. [PMID: 26328834 PMCID: PMC4552706 DOI: 10.1063/1.4929335] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2015] [Accepted: 08/10/2015] [Indexed: 11/14/2022] Open
Abstract
Semi-empirical quantum mechanical methods traditionally expand the electron density in a minimal, valence-only electron basis set. The minimal-basis approximation causes molecular polarization to be underestimated, and hence intermolecular interaction energies are also underestimated, especially for intermolecular interactions involving charged species. In this work, the third-order self-consistent charge density functional tight-binding method (DFTB3) is augmented with an auxiliary response density using the chemical-potential equalization (CPE) method and an empirical dispersion correction (D3). The parameters in the CPE and D3 models are fitted to high-level CCSD(T) reference interaction energies for a broad range of chemical species, as well as dipole moments calculated at the DFT level; the impact of including polarizabilities of molecules in the parameterization is also considered. Parameters for the elements H, C, N, O, and S are presented. The Root Mean Square Deviation (RMSD) interaction energy is improved from 6.07 kcal/mol to 1.49 kcal/mol for interactions with one charged species, whereas the RMSD is improved from 5.60 kcal/mol to 1.73 for a set of 9 salt bridges, compared to uncorrected DFTB3. For large water clusters and complexes that are dominated by dispersion interactions, the already satisfactory performance of the DFTB3-D3 model is retained; polarizabilities of neutral molecules are also notably improved. Overall, the CPE extension of DFTB3-D3 provides a more balanced description of different types of non-covalent interactions than Neglect of Diatomic Differential Overlap type of semi-empirical methods (e.g., PM6-D3H4) and PBE-D3 with modest basis sets.
Collapse
Affiliation(s)
- Anders S Christensen
- Department of Chemistry, University of Wisconsin-Madison, 1101 University Ave., Madison, Wisconsin 53706, USA
| | - Marcus Elstner
- Theoretische Chemische Biologie, Universität Karlsruhe, Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - Qiang Cui
- Department of Chemistry, University of Wisconsin-Madison, 1101 University Ave., Madison, Wisconsin 53706, USA
| |
Collapse
|