1
|
Hillers-Bendtsen AE, Todarwal Y, Norman P, Mikkelsen KV. Dynamical Effects of Solvation on Norbornadiene/Quadricyclane Systems. J Phys Chem A 2024; 128:2602-2610. [PMID: 38511966 DOI: 10.1021/acs.jpca.4c00045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/22/2024]
Abstract
Molecules that can undergo reversible chemical transformations following the absorption of light, the so-called molecular photoswitches, have attracted increasing attention in technologies, such as solar energy storage. Here, the optical and thermochemical properties of the photoswitch are central to its applicability, and these properties are influenced significantly by solvation. We investigate the effects of solvation on two norbornadiene/quadricyclane photoswitches. Emphasis is put on the energy difference between the two isomers and the optical absorption as these are central to the application of the systems in solar energy storage. Using a combined classical molecular dynamics and quantum mechanical/molecular mechanical computational scheme, we showcase that the dynamic effects of solvation are important. In particular, it is found that standard implicit solvation models generally underestimate the energy difference between the two isomers and overestimate the strength of the absorption, while the explicit solvation spectra are also less red-shifted than those obtained using implicit solvation models. We also find that the absorption spectra of the two systems are strongly correlated with specific dihedral angles. Altogether, this highlights the importance of including the dynamic effects of solvation.
Collapse
Affiliation(s)
| | - Yogesh Todarwal
- Division of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden
| | - Patrick Norman
- Division of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden
| | - Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100 Copenhagen Ø, Denmark
| |
Collapse
|
2
|
Sülzner N, Hättig C. Role of Singles Amplitudes in ADC(2) and CC2 for Low-Lying Electronically Excited States. J Chem Theory Comput 2024; 20:2462-2474. [PMID: 38449383 DOI: 10.1021/acs.jctc.3c01355] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/08/2024]
Abstract
The closely related second-order methods CC2 and ADC(2) usually perform very similarly for single excitations of organic molecules. However, as rationalized in this work, significant deviations between these two methods can arise if the ground state and a low-lying singly excited state arise from a strong coupling between their leading configurations. Such a configuration mixing is partially accounted for in CC2 through the ground-state singles amplitudes but is omitted in ADC(2). This can cause unusual deviations between the results obtained with these methods. In this work, we study how severe this effect can become at the example of two solvatochromic dyes: the negatively solvatochromic betaine dye N1-tBu and the positively solvatochromic bithiophene P1. These two dyes allow one to study the limits of both small and somewhat larger excitation energies and configuration mixing by tuning the S0 → S1 transition energy through the polarity of the environment. Higher-level calculations at the CC3 level provide information on the accuracy of ADC(2) and CC2 in these cases. The most extreme deviation between ADC(2) and CC2 is found for N1-tBu in vacuum, where the ADC(2) result is 0.45 eV below that of CC2. In this case, the methodical error of CC2 with respect to CC3 is only 0.05 eV. With increasing excitation energy in polar solvents, the CC2-ADC(2) deviation decreases and reaches a value of only 0.15 eV. For P1, which has larger excitation energies, these effects are reversed due to the opposite solvatochromism but also smaller in magnitude: the deviation increases from 0.08 eV in vacuum to 0.16 eV in the so-called conductor limit of the continuum solvation model. Although for these two dyes larger deviations are observed for smaller excitation energies, the extent of configuration mixing does not generally correlate with only the size of excitation energy. For example, s-triazine (0.15 eV), formamide (0.19 eV), and formaldehyde (0.23 eV) also show large deviations between CC2 and ADC(2) despite their much higher excitation energies compared to those of N1-tBu and P1.
Collapse
Affiliation(s)
- Niklas Sülzner
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
3
|
Hillers-Bendtsen AE, Zhou Y, Mikkelsen KV. Investigation of Solvent Effects on the Molecular Energy Storage and Optical Properties of Bicyclooctadiene/Tetracyclooctane Photoswitches. J Phys Chem A 2024; 128:41-50. [PMID: 38152898 DOI: 10.1021/acs.jpca.3c04760] [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/2023]
Abstract
In this paper, we investigate the effects of solvation on the solar energy storage properties of bicyclooctadiene/tetracyclooctane (BOD/TCO) photoswitches. The solvent effects on the thermochemical and optical properties are studied in cyclohexane, toluene, dichloromethane, ethanol, acetonitrile, and a vacuum using density functional theory and coupled cluster theory. Our results show that the energy storage capacity of the BOD/TCO system increases as the solvent polarity increases, and the change is more significant with an unsubstituted system. The energy storage capacity of the substituted system is not dependent on the polarity of the solvent. As the solvent polarity increases, the absorption peaks shift away from each other and the absorption intensities increase. Overall, the solvents improve the performance of the optical properties and the energy storage capacities of the BOD/TCO molecular solar thermal systems.
Collapse
Affiliation(s)
| | - Yiwei Zhou
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100 Copenhagen Ø, Denmark
| | - Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100 Copenhagen Ø, Denmark
| |
Collapse
|
4
|
Pedersen J, Rasmussen MH, Mikkelsen KV. Redfield Propagation of Photoinduced Electron Transfer Reactions in Vacuum and Solution. J Chem Theory Comput 2022; 18:7052-7072. [PMID: 36413807 DOI: 10.1021/acs.jctc.2c00538] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Dynamical simulations of ultrafast electron transfer reactions are of utmost interest. To allow for energy dissipation directly into an external surrounding environment, a solvent coupling model has been deduced, implemented, and utilized to describe the photoinduced electron transfer dynamics within a model triad system herein. The model is based on Redfield theory, and the environment is represented by harmonic oscillators filled with bosonic quanta. To imitate real solvents, the oscillators have been equipped with frequencies and polarization lifetimes characteristic of the corresponding solvent. The population was found to transfer through the energetically lowest electron transfer route regardless of the medium. The condensed population transfer dynamics were observed to be highly dependent on the solvent parameters. In particular, an increase in the solvent coupling entailed a detainment in the population transfer from the initially prepared diabatic state and a promotion in the population transfer through the other electron transfer route. Two explanations based on the diagonal and off-diagonal matrix elements of the Kohn-Sham Fock matrix, respectively, have been provided. The lifetime of the populated partially charge-separated state was prolonged with increasing solvent polarity, and it was explained in terms of attractive interactions between the solvent's dipole moments and the fragments' charges. The high-frequency vibrational fine-structure in the correlation function was demonstrated to be important for the transfer dynamics, and the importance of dephasing effects in polar solvents was verified and precised to concern the optical polarization of the solvents.
Collapse
Affiliation(s)
- Jacob Pedersen
- Department of Chemistry, University of Copenhagen, Copenhagen, DK-2100, Denmark
| | - Maria H Rasmussen
- Department of Chemistry, University of Copenhagen, Copenhagen, DK-2100, Denmark
| | - Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, Copenhagen, DK-2100, Denmark
| |
Collapse
|
5
|
Høyer NM, Kjeldal FØ, Hillers-Bendtsen AE, Mikkelsen KV, Olsen J, Jorgensen P. Cluster perturbation theory. VI. Ground-state energy series using the Lagrangian. J Chem Phys 2022; 157:024106. [DOI: 10.1063/5.0082583] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we derive alternative cluster perturbation series to the energy cluster perturbation (ECP) series of Paper I. The ECP series were derived using the standard coupled cluster energy framework. Here, we use the coupled cluster Lagrangian framework to derive the Lagrangian cluster perturbation (LCP) series and show that a slightly modified order concept means that the ECP and the LCP series become identical. Using the Lagrangian, we also derive a perturbation series where total cluster amplitudes and multipliers are determined through the same orders as dictated by the 2n+1/2n+2 rule, the VCP series. The VCP energies have errors that are bilinear in the errors of the total cluster amplitudes and multipliers. Test calculations have been performed for S(D) and SD(T) orbital excitation spaces. The convergence of the test calculations can be divided into two groups. The first group contains molecules that exhibit a slow monotonic geometric convergence pattern in both orbital spaces. The second group, containing the rest of the molecules, exhibits a fast low order convergence. For the S(D) calculations, the deviations in fourth order from the coupled cluster singles and doubles energy are below 1.3 per cent for the LCP series and 0.3 per cent for the VCP series, respectively. For the SD(T) calculations, the second group hardly shows any difference between the low order convergence of the three series ECP, LCP, and LCP and at fourth order the deviations from the triples correlation energy are less than 1 per cent.
Collapse
Affiliation(s)
| | | | | | - Kurt V. Mikkelsen
- Department of Chemistry, University of Copenhagen Institute of Chemistry, Denmark
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University Department of Chemistry, Denmark
| | | |
Collapse
|
6
|
Garcia-Ratés M, Becker U, Neese F. Implicit solvation in domain based pair natural orbital coupled cluster (DLPNO-CCSD) theory. J Comput Chem 2021; 42:1959-1973. [PMID: 34347890 DOI: 10.1002/jcc.26726] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Revised: 07/08/2021] [Accepted: 07/14/2021] [Indexed: 11/11/2022]
Abstract
A nearly linear scaling implementation of coupled-cluster with singles and doubles excitations (CCSD) can be achieved by means of the domain-based local pair natural orbital (DLPNO) method. The combination of DLPNO-CCSD with implicit solvation methods allows the calculation of accurate energies and chemical properties of solvated systems at an affordable computational cost. We have efficiently implemented different schemes within the conductor-like polarizable continuum model (C-PCM) for DLPNO-CCSD in the ORCA quantum chemistry suite. In our implementation, the overhead due to the additional solvent terms amounts to less than 5% of the time the equivalent gas phase job takes. Our results for organic neutrals and open-shell ions in water show that for most systems, adding solvation terms to the coupled-cluster amplitudes equations and to the energy leads to small changes in the total energy compared to only considering solvated orbitals and corrections to the reference energy. However, when the solute contains certain functional groups, such as carbonyl or nitrile groups, the changes in the energy are larger and estimated to be around 0.04 and 0.02 kcal/mol for each carbonyl and nitrile group in the solute, respectively. For solutes containing metals, the use of accurate CC/C-PCM schemes is crucial to account for correlation solvation effects. Simultaneously, we have calculated the electrostatic component of the solvation energy for neutrals and ions in water for the different DLPNO-CCSD/C-PCM schemes. We observe negligible changes in the deviation between DLPNO-CCSD and canonical-CCSD data. Here, DLPNO-CCSD results outperform those for Hartree-Fock and density functional theory calculations.
Collapse
Affiliation(s)
| | - Ute Becker
- Max-Planck-Institut für Kohlenforschung, Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Mülheim an der Ruhr, Germany
| |
Collapse
|
7
|
Caricato M. Coupled cluster theory in the condensed phase within the singles‐T density scheme for the environment response. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1463] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Marco Caricato
- Department of Chemistry University of Kansas Lawrence Kansas
| |
Collapse
|
8
|
Hansen MB, Madsen NK, Zoccante A, Christiansen O. Time-dependent vibrational coupled cluster theory: Theory and implementation at the two-mode coupling level. J Chem Phys 2019; 151:154116. [DOI: 10.1063/1.5117207] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Mads Bøttger Hansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| | - Niels Kristian Madsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| | - Alberto Zoccante
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| |
Collapse
|
9
|
Izsák R. Single‐reference coupled cluster methods for computing excitation energies in large molecules: The efficiency and accuracy of approximations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1445] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Róbert Izsák
- Max‐Planck‐Institut für Kohlenforschung Mülheim an der Ruhr Germany
| |
Collapse
|
10
|
Ren S, Lipparini F, Mennucci B, Caricato M. Coupled Cluster Theory with Induced Dipole Polarizable Embedding for Ground and Excited States. J Chem Theory Comput 2019; 15:4485-4496. [PMID: 31265278 DOI: 10.1021/acs.jctc.9b00468] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
In this work, we present the theory and implementation of the coupled cluster single and double excitations (CCSD) method combined with a classical polarizable molecular mechanics force field (MMPol) based on the induced dipole model. The method is developed to compute electronic excitation energies within the state specific (SS) and linear response (LR) formalisms for the interaction of the quantum mechanical and classical regions. Furthermore, we consider an approximate expression of the correlation energy, originally developed for CCSD with implicit solvation models, where the interaction term is linear in the coupled cluster density. This approximation allows us to include the explicit contribution of the environment to the CC equations without increasing the computational effort. The test calculations on microsolvated systems, where the CCSD/MMPol method is compared to full CCSD calculations, demonstrates the reliability of this computational protocol for all interaction schemes (errors < 2%). We also show that it is important to include induced dipoles on all atom centers of the classical region and that too diffuse functions in the basis set may be problematic due to too strong interaction with the environment.
Collapse
Affiliation(s)
- Sijin Ren
- Department of Chemistry , University of Kansas , 1567 Irving Hill Road , Lawrence , Kansas 66044 , United States
| | - Filippo Lipparini
- Department of Chemistry , Università di Pisa , Via Giuseppe Moruzzi , 13 56124 Pisa , Italy
| | - Benedetta Mennucci
- Department of Chemistry , Università di Pisa , Via Giuseppe Moruzzi , 13 56124 Pisa , Italy
| | - Marco Caricato
- Department of Chemistry , University of Kansas , 1567 Irving Hill Road , Lawrence , Kansas 66044 , United States
| |
Collapse
|
11
|
Haqghu M, Irani M, Gholami MR. Theoretical Study of Kinetics and Mechanism of Reactions of Hydroxylamine and Amineoxide Anion with Methyl Iodide in Gas and Aqueous Phases. PROGRESS IN REACTION KINETICS AND MECHANISM 2019. [DOI: 10.3184/146867807x227462] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
The kinetics and mechanism of the reactions of hydroxylamine and aminoxide anion with methyl iodide were studied with ab initio calculations, Monte-Carlo and QM/MM simulations in gas and aqueous phases. Geometrical parameters and charge calculations show that these reactions proceed through the SN2 (bimolecular nucleophilic displacement) mechanism only. The solvent effects on these reactions were studied by inserting water molecules in reaction media, Onsager model, Monte-Carlo and QM/MM simulations. Activation parameters indicate the expected variation in activation energy and rate coefficient in aqueous phase in comparison to the gas phase. The shift of potential energy barrier through the reactants or products for the studied reactions in the gas phase is in the opposite direction in comparison to the aqueous phase.
Collapse
Affiliation(s)
- M. Haqghu
- Department of Chemistry, Sharif University of Technology, P.O. BOX 11365-9516,Tehran, Iran
| | - M. Irani
- Department of Chemistry, Sharif University of Technology, P.O. BOX 11365-9516,Tehran, Iran
| | - M. R. Gholami
- Department of Chemistry, Sharif University of Technology, P.O. BOX 11365-9516,Tehran, Iran
| |
Collapse
|
12
|
Heterogeneous nucleation of polymorphs on polymer surfaces: polymer-molecule interactions using a heterogeneous dielectric solvation model. J Mol Model 2018; 24:156. [PMID: 29886563 DOI: 10.1007/s00894-018-3657-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2017] [Accepted: 04/13/2018] [Indexed: 10/14/2022]
Abstract
We have investigated the mechanism of the nucleation of acetaminophen on poly(methyl-methacrylate) and poly(vinyl-acetate) utilizing a combination of quantum mechanical computations and electrostatic models. We have used a heterogeneous dielectric solvation model to determine the stability of different orientations of acetaminophen on polymer surfaces. We find that for the nucleation of acetaminophen on the polymer surfaces in vacuum, the most stable orientation is a flat orientation. For the nucleation process in solution where acetaminophen and the polymer surface are surrounded by a solvent, we find that the heterogeneous dielectric solvation model predicts that a sideways orientation is the most stable orientation.
Collapse
|
13
|
Caricato M. Linear response coupled cluster theory with the polarizable continuum model within the singles approximation for the solvent response. J Chem Phys 2018; 148:134113. [DOI: 10.1063/1.5021781] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Affiliation(s)
- Marco Caricato
- Department of Chemistry, University of Kansas, 1251 Wescoe Hall Dr., Lawrence, Kansas 66045, USA
| |
Collapse
|
14
|
Molecular Properties of Sandwiched Molecules Between Electrodes and Nanoparticles. ADVANCES IN QUANTUM CHEMISTRY 2017. [DOI: 10.1016/bs.aiq.2017.03.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
15
|
Ren S, Harms J, Caricato M. An EOM-CCSD-PCM Benchmark for Electronic Excitation Energies of Solvated Molecules. J Chem Theory Comput 2016; 13:117-124. [PMID: 27973775 DOI: 10.1021/acs.jctc.6b01053] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this work, we benchmark the equation of motion coupled cluster with single and double excitations (EOM-CCSD) method combined with the polarizable continuum model (PCM) for the calculation of electronic excitation energies of solvated molecules. EOM-CCSD is one of the most accurate methods for computing one-electron excitation energies, and accounting for the solvent effect on this property is a key challenge. PCM is one of the most widely employed solvation models due to its adaptability to virtually any solute and its efficient implementation with density functional theory methods (DFT). Our goal in this work is to evaluate the reliability of EOM-CCSD-PCM, especially compared to time-dependent DFT-PCM (TDDFT-PCM). Comparisons between calculated and experimental excitation energies show that EOM-CCSD-PCM consistently overestimates experimental results by 0.4-0.5 eV, which is larger than the expected EOM-CCSD error in vacuo. We attribute this decrease in accuracy to the approximated solvation model. Thus, we investigate a particularly important source of error: the lack of H-bonding interactions in PCM. We show that this issue can be addressed by computing an energy shift, ΔHB, from bare-PCM to microsolvation + PCM at DFT level. Our results show that such a shift is independent of the functional used, contrary to the absolute value of the excitation energy. Hence, we suggest an efficient protocol where the EOM-CCSD-PCM transition energy is corrected by ΔHB(DFT), which consistently improves the agreement with the experimental measurements.
Collapse
Affiliation(s)
- Sijin Ren
- Department of Chemistry, University of Kansas , 1251 Wescoe Hall Dr., Lawrence, Kansas 66045, United States
| | - Joseph Harms
- Lawrence High School , 1901 Louisiana St., Lawrence, Kansas 66046, United States
| | - Marco Caricato
- Department of Chemistry, University of Kansas , 1251 Wescoe Hall Dr., Lawrence, Kansas 66045, United States
| |
Collapse
|
16
|
Krause K, Klopper W. Communication: A simplified coupled-cluster Lagrangian for polarizable embedding. J Chem Phys 2016; 144:041101. [PMID: 26827193 DOI: 10.1063/1.4940895] [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/22/2022] Open
Abstract
A simplified coupled-cluster Lagrangian, which is linear in the Lagrangian multipliers, is proposed for the coupled-cluster treatment of a quantum mechanical system in a polarizable environment. In the simplified approach, the amplitude equations are decoupled from the Lagrangian multipliers and the energy obtained from the projected coupled-cluster equation corresponds to a stationary point of the Lagrangian.
Collapse
Affiliation(s)
- Katharina Krause
- Karlsruhe Institute of Technology (KIT), Institute of Physical Chemistry, Theoretical Chemistry Group, KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Wim Klopper
- Karlsruhe Institute of Technology (KIT), Institute of Physical Chemistry, Theoretical Chemistry Group, KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| |
Collapse
|
17
|
Lipparini F, Mennucci B. Perspective: Polarizable continuum models for quantum-mechanical descriptions. J Chem Phys 2016; 144:160901. [DOI: 10.1063/1.4947236] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Filippo Lipparini
- Institut für Physikalische Chemie, Universität Mainz, Duesbergweg 10-14, D55128 Mainz, Germany
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
18
|
Casanova D, Gusarov S, Kovalenko A, Ziegler T. Evaluation of the SCF Combination of KS-DFT and 3D-RISM-KH; Solvation Effect on Conformational Equilibria, Tautomerization Energies, and Activation Barriers. J Chem Theory Comput 2015; 3:458-76. [PMID: 26637028 DOI: 10.1021/ct6001785] [Citation(s) in RCA: 54] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The effect of solvation on conformational equilibria, tautomerization energies, and activation barriers in simple SN2 reactions is reproduced by using the self-consistent field coupling of the Kohn-Sham density functional theory (KS-DFT) for electronic structure and the three-dimensional reference interaction site model with the closure approximation of Kovalenko and Hirata (3D-RISM-KH) for molecular solvation structure. These examples are used in order to validate the implementation of the 3D-RISM-KH method in the Amsterdam Density Functional package. The computations of the free energy difference in the trans/gauche conformational equilibrium for 1,2-dichloroethane in different solvents; the relative tautomerization free energy for cytosine, isocytosine, and guanine; and the free energy activation barrier for a CH3X-type (X = F, Cl, Br) SN2 reaction exhibit agreement with the experimental data. The method is also applied to the electronic and hydration structure of carbon single-wall nanotubes of different diameters, including the effect of water located in the inner space of the nanotubes. A comparison with continuum models of solvation (including COSMO) as well as with other more precise and computationally more expensive calculations is made to demonstrate the accuracy and predictive capability of the new KS-DFT/3D-RISM-KH method.
Collapse
Affiliation(s)
- David Casanova
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| | - Sergey Gusarov
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| | - Andriy Kovalenko
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| | - Tom Ziegler
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| |
Collapse
|
19
|
List NH, Coriani S, Kongsted J, Christiansen O. Lanczos-driven coupled–cluster damped linear response theory for molecules in polarizable environments. J Chem Phys 2014; 141:244107. [DOI: 10.1063/1.4903981] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Nanna Holmgaard List
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, 5230 Odense M, Denmark
| | - Sonia Coriani
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Trieste, via Giorgieri 1, 34127 Trieste, Italy
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, 5230 Odense M, Denmark
| | - Ove Christiansen
- Center for Oxygen Microscopy and Imaging, Department of Chemistry, Aarhus University, Langelandsgade 140, Aarhus 8000, Denmark
| |
Collapse
|
20
|
Azizi A, Ebrahimi A, Habibi-Khorassani M, Rezazadeh S, Behazin R. The Effects of Interactions of Dicarboxylic Acids on the Stability of the Caffeine Molecule: A Theoretical Study. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2014. [DOI: 10.1246/bcsj.20130287] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Abolfazl Azizi
- Department of Chemistry, University of Sistan and Baluchestan
| | - Ali Ebrahimi
- Department of Chemistry, University of Sistan and Baluchestan
| | | | - Shiva Rezazadeh
- Department of Chemistry, University of Sistan and Baluchestan
| | - Roya Behazin
- Department of Chemistry, University of Sistan and Baluchestan
| |
Collapse
|
21
|
Caricato M. A corrected-linear response formalism for the calculation of electronic excitation energies of solvated molecules with the CCSD-PCM method. COMPUT THEOR CHEM 2014. [DOI: 10.1016/j.comptc.2014.02.001] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
22
|
Aidas K, Angeli C, Bak KL, Bakken V, Bast R, Boman L, Christiansen O, Cimiraglia R, Coriani S, Dahle P, Dalskov EK, Ekström U, Enevoldsen T, Eriksen JJ, Ettenhuber P, Fernández B, Ferrighi L, Fliegl H, Frediani L, Hald K, Halkier A, Hättig C, Heiberg H, Helgaker T, Hennum AC, Hettema H, Hjertenæs E, Høst S, Høyvik IM, Iozzi MF, Jansík B, Jensen HJA, Jonsson D, Jørgensen P, Kauczor J, Kirpekar S, Kjærgaard T, Klopper W, Knecht S, Kobayashi R, Koch H, Kongsted J, Krapp A, Kristensen K, Ligabue A, Lutnæs OB, Melo JI, Mikkelsen KV, Myhre RH, Neiss C, Nielsen CB, Norman P, Olsen J, Olsen JMH, Osted A, Packer MJ, Pawlowski F, Pedersen TB, Provasi PF, Reine S, Rinkevicius Z, Ruden TA, Ruud K, Rybkin VV, Sałek P, Samson CCM, de Merás AS, Saue T, Sauer SPA, Schimmelpfennig B, Sneskov K, Steindal AH, Sylvester-Hvid KO, Taylor PR, Teale AM, Tellgren EI, Tew DP, Thorvaldsen AJ, Thøgersen L, Vahtras O, Watson MA, Wilson DJD, Ziolkowski M, Agren H. The Dalton quantum chemistry program system. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013; 4:269-284. [PMID: 25309629 PMCID: PMC4171759 DOI: 10.1002/wcms.1172] [Citation(s) in RCA: 874] [Impact Index Per Article: 79.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, Møller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.
Collapse
Affiliation(s)
- Kestutis Aidas
- Department of General Physics and Spectroscopy, Faculty of Physics, Vilnius University Vilnius, Lithuania
| | | | - Keld L Bak
- Aarhus University School of Engineering Aarhus, Denmark
| | - Vebjørn Bakken
- Faculty of Mathematics and Natural Sciences, University of Oslo Oslo, Norway
| | - Radovan Bast
- Department of Theoretical Chemistry and Biology, School of Biotechnology, KTH Royal Institute of Technology Stockholm, Sweden
| | | | | | | | - Sonia Coriani
- Department of Chemical and Pharmaceutical Sciences, University of Trieste Trieste, Italy
| | - Pål Dahle
- Norwegian Computing Center Oslo, Norway
| | | | - Ulf Ekström
- CTCC, Department of Chemistry, University of Oslo Oslo, Norway
| | - Thomas Enevoldsen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark Odense, Denmark
| | | | | | - Berta Fernández
- Department of Physical Chemistry and Center for Research in Biological Chemistry and Molecular Materials (CIQUS), University of Santiago de Compostela Santiago de Compostela, Spain
| | - Lara Ferrighi
- CTCC, Department of Chemistry, UiT The Arctic University of Norway, Tromsø Norway
| | - Heike Fliegl
- CTCC, Department of Chemistry, University of Oslo Oslo, Norway
| | - Luca Frediani
- CTCC, Department of Chemistry, UiT The Arctic University of Norway, Tromsø Norway
| | | | | | - Christof Hättig
- Department of Theoretical Chemistry, Ruhr-University Bochum Bochum, Germany
| | | | - Trygve Helgaker
- CTCC, Department of Chemistry, University of Oslo Oslo, Norway
| | | | - Hinne Hettema
- Department of Philosophy, The University of Auckland Auckland, New Zealand
| | - Eirik Hjertenæs
- Department of Chemistry, Norwegian University of Science and Technology Trondheim, Norway
| | - Stinne Høst
- Department of Geoscience, Aarhus University Aarhus, Denmark
| | | | | | | | - Hans Jørgen Aa Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark Odense, Denmark
| | - Dan Jonsson
- High-Performance Computing Group, UiT The Arctic University of Norway, Tromsø Norway
| | - Poul Jørgensen
- Department of Chemistry, Aarhus University Aarhus, Denmark
| | - Joanna Kauczor
- Department of Physics, Chemistry and Biology, Linköping University Linköping, Sweden
| | | | | | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology Karlsruhe, Germany
| | - Stefan Knecht
- Laboratory of Physical Chemistry, ETH Zürich Zürich, Switzerland
| | - Rika Kobayashi
- Australian National University Supercomputer Facility Canberra, Australia
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology Trondheim, Norway
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark Odense, Denmark
| | | | | | - Andrea Ligabue
- Computer Services: Networks and Systems, University of Modena and Reggio Emilia Modena, Italy
| | | | - Juan I Melo
- Physics Department, FCEyN-UBA and IFIBA-CONICET, Universidad de Buenos Aires Buenos Aires, Argentina
| | - Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, Copenhagen Denmark
| | - Rolf H Myhre
- Department of Chemistry, Norwegian University of Science and Technology Trondheim, Norway
| | - Christian Neiss
- Department of Chemistry and Pharmacy, Friedrich-Alexander University Erlangen-Nürnberg Erlangen, Germany
| | | | - Patrick Norman
- Department of Physics, Chemistry and Biology, Linköping University Linköping, Sweden
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University Aarhus, Denmark
| | - Jógvan Magnus H Olsen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark Odense, Denmark
| | | | - Martin J Packer
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark Odense, Denmark
| | - Filip Pawlowski
- Institute of Physics, Kazimierz Wielki University Bydgoszcz, Poland
| | | | - Patricio F Provasi
- Department of Physics, University of Northeastern and IMIT-CONICET Corrientes, Argentina
| | - Simen Reine
- CTCC, Department of Chemistry, University of Oslo Oslo, Norway
| | - Zilvinas Rinkevicius
- Department of Theoretical Chemistry and Biology, School of Biotechnology and Swedish e-Science Research Center (SeRC), KTH Royal Institute of Technology Stockholm, Sweden
| | | | - Kenneth Ruud
- CTCC, Department of Chemistry, UiT The Arctic University of Norway, Tromsø Norway
| | - Vladimir V Rybkin
- Institute of Physical Chemistry, Karlsruhe Institute of Technology Karlsruhe, Germany
| | | | - Claire C M Samson
- Institute of Physical Chemistry, Karlsruhe Institute of Technology Karlsruhe, Germany
| | | | - Trond Saue
- Paul Sabatier University Toulouse, France
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, Copenhagen Denmark
| | - Bernd Schimmelpfennig
- Institute for Nuclear Waste Disposal, Karlsruhe Institute of Technology Karlsruhe, Germany
| | | | - Arnfinn H Steindal
- CTCC, Department of Chemistry, UiT The Arctic University of Norway, Tromsø Norway
| | | | - Peter R Taylor
- VLSCI and School of Chemistry, University of Melbourne Parkville, Australia
| | - Andrew M Teale
- School of Chemistry, University of Nottingham Nottingham, UK
| | - Erik I Tellgren
- CTCC, Department of Chemistry, University of Oslo Oslo, Norway
| | - David P Tew
- School of Chemistry, University of Bristol Bristol, UK
| | | | | | - Olav Vahtras
- Department of Theoretical Chemistry and Biology, School of Biotechnology, KTH Royal Institute of Technology Stockholm, Sweden
| | - Mark A Watson
- Department of Chemistry, Princeton University Princeton, New Jersey
| | - David J D Wilson
- Department of Chemistry and La Trobe Institute for Molecular Sciences, La Trobe University Melbourne, Australia
| | - Marcin Ziolkowski
- CoE for Next Generation Computing, Clemson University Clemson, South Carolina
| | - Hans Agren
- Department of Theoretical Chemistry and Biology, School of Biotechnology, KTH Royal Institute of Technology Stockholm, Sweden
| |
Collapse
|
23
|
Caricato M. Implementation of the CCSD-PCM linear response function for frequency dependent properties in solution: Application to polarizability and specific rotation. J Chem Phys 2013; 139:114103. [DOI: 10.1063/1.4821087] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
Caricato M. A comparison between state-specific and linear-response formalisms for the calculation of vertical electronic transition energy in solution with the CCSD-PCM method. J Chem Phys 2013; 139:044116. [DOI: 10.1063/1.4816482] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
Caricato M, Lipparini F, Scalmani G, Cappelli C, Barone V. Vertical Electronic Excitations in Solution with the EOM-CCSD Method Combined with a Polarizable Explicit/Implicit Solvent Model. J Chem Theory Comput 2013; 9:3035-3042. [PMID: 26504458 DOI: 10.1021/ct4003288] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The accurate calculation of electronic transition energies and properties of isolated chromophores is not sufficient to provide a realistic simulation of their excited states in solution. In fact, the solvent influences the solute geometry, electronic structure, and response to external fields. Therefore, a proper description of the solvent effect is fundamental. This can be achieved by combining polarizable explicit and implicit representations of the solvent. The former provides a realistic description of solvent molecules around the solute, while the latter introduces the electrostatic effect of the bulk solution and reduces the need of too large a number of explicit solvent molecules. This strategy is particularly appealing when an accurate method such as equation of motion coupled cluster singles and doubles (EOM-CCSD) is employed for the treatment of the chromophore. In this contribution, we present the coupling of EOM-CCSD with a fluctuating charges (FQ) model and polarizable continuum model (PCM) of solvation for vertical excitations in a state-specific framework. The theory, implementation, and prototypical applications of the method are presented. Numerical tests on small solute-water clusters show very good agreement between full EOM-CCSD and EOM-CCSD-FQ calculations, with and without PCM, with differences ≤ 0.1 eV. Additionally, approximated schemes that further reduce the computational cost of the method are introduced and showed to perform well compared to the full method (errors ≤ 0.1 eV).
Collapse
Affiliation(s)
- Marco Caricato
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| | | | - Giovanni Scalmani
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy. ; Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Risorgimento, 35 I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
26
|
Caricato M. Absorption and Emission Spectra of Solvated Molecules with the EOM-CCSD-PCM Method. J Chem Theory Comput 2012; 8:4494-502. [PMID: 26605609 DOI: 10.1021/ct3006997] [Citation(s) in RCA: 60] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The accurate calculation of transition energies and properties of isolated molecules is not enough for realistic simulations of their absorption and emission spectra in solution. In fact, the solvent influences the solute geometry, electronic structure, and response to external fields, and a proper description of the solvent effect is fundamental. However, the computational cost of including explicit solvent molecules around the solute becomes rather onerous when an accurate method such as the equation of motion coupled cluster singles and doubles (EOM-CCSD) is employed. The polarizable continuum model of solvation (PCM) may provide an efficient alternative to explicit models, since the sampling of solvent configurations is implicit and the solute-solvent mutual polarization is naturally accounted for. In this contribution, the absorption and emission spectra of molecules in solution are modeled through the EOM-CCSD-PCM method. The equilibrium solvation regime is employed for the geometry optimization of the solute molecule in the ground and excited states, while the nonequilibrium solvation regime is employed for vertical transitions. The theory, implementation, and prototypical applications of the method are presented. The numerical tests involve solvents that are particularly challenging for PCM: low-polar and protic polar solvents. Nonetheless, the experimental trends are well reproduced, and the overall agreement with the measured data is remarkable.
Collapse
Affiliation(s)
- Marco Caricato
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, United States
| |
Collapse
|
27
|
Dieterich JM, Oliveira JCA, Mata RA. Application of Local Second-Order Møller–Plesset Perturbation Theory to the Study of Structures in Solution. J Chem Theory Comput 2012; 8:3053-60. [DOI: 10.1021/ct300009n] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Johannes M. Dieterich
- Institut für
Physikalische
Chemie, Universität Göttingen, Tammannstrasse 6, D-37077 Göttingen, Germany
| | - João C. A. Oliveira
- Institut für
Physikalische
Chemie, Universität Göttingen, Tammannstrasse 6, D-37077 Göttingen, Germany
| | - Ricardo A. Mata
- Institut für
Physikalische
Chemie, Universität Göttingen, Tammannstrasse 6, D-37077 Göttingen, Germany
| |
Collapse
|
28
|
Schwabe T, Sneskov K, Haugaard Olsen JM, Kongsted J, Christiansen O, Hättig C. PERI–CC2: A Polarizable Embedded RI-CC2 Method. J Chem Theory Comput 2012; 8:3274-83. [DOI: 10.1021/ct3003749] [Citation(s) in RCA: 66] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Tobias Schwabe
- Center for Bioinformatics and
Institute of Physical Chemistry, University of Hamburg, Bundesstraße
43, D-20146 Hamburg, Germany
| | - Kristian Sneskov
- Center for Oxygen Microscopy
and Imaging, Department of Chemistry, Aarhus University, Langelandsgade
140, DK-8000 Aarhus C, Denmark
- The Lundbeck Foundation Center
for Theoretical Chemistry, Department of Chemistry, Aarhus University,
Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Jógvan Magnus Haugaard Olsen
- Department of Physics, Chemistry
and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230
Odense M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry
and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230
Odense M, Denmark
| | - Ove Christiansen
- Center for Oxygen Microscopy
and Imaging, Department of Chemistry, Aarhus University, Langelandsgade
140, DK-8000 Aarhus C, Denmark
- The Lundbeck Foundation Center
for Theoretical Chemistry, Department of Chemistry, Aarhus University,
Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Christof Hättig
- Lehrstuhl
für Theoretische
Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| |
Collapse
|
29
|
Caricato M. Exploring Potential Energy Surfaces of Electronic Excited States in Solution with the EOM-CCSD-PCM Method. J Chem Theory Comput 2012; 8:5081-91. [DOI: 10.1021/ct300382a] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Marco Caricato
- Gaussian, Inc., 340
Quinnipiac St., Bldg. 40, Wallingford,
Connecticut 06492, United States
| |
Collapse
|
30
|
Forti F, Cavasotto CN, Orozco M, Barril X, Luque FJ. A Multilevel Strategy for the Exploration of the Conformational Flexibility of Small Molecules. J Chem Theory Comput 2012; 8:1808-19. [DOI: 10.1021/ct300097s] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Flavio Forti
- Departament de Fisicoquímica
and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Av. Diagonal 643, E-08028,
Barcelona, Spain
| | - Claudio N. Cavasotto
- School
of Biomedical Informatics, University of Texas Health Science Center at Houston, 7000 Fannin Street Ste. 690,
Houston, Texas 77030, United States
| | - Modesto Orozco
- Joint IRB-BSC Program in Computational
Biology, Institute of Research in Biomedicine Barcelona, Barcelona Scientific Park, Josep Samitier 1-6, 08028 barcelona,
Spain
| | - Xavier Barril
- Departament de Fisicoquímica
and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Av. Diagonal 643, E-08028,
Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA),
Barcelona, Spain
| | - F. Javier Luque
- Departament de Fisicoquímica
and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Av. Diagonal 643, E-08028,
Barcelona, Spain
| |
Collapse
|
31
|
Helgaker T, Coriani S, Jørgensen P, Kristensen K, Olsen J, Ruud K. Recent Advances in Wave Function-Based Methods of Molecular-Property Calculations. Chem Rev 2012; 112:543-631. [DOI: 10.1021/cr2002239] [Citation(s) in RCA: 463] [Impact Index Per Article: 38.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Trygve Helgaker
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Sonia Coriani
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Trieste, Via Giorgieri 1, I-34127 Trieste, Italy
| | - Poul Jørgensen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Kasper Kristensen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Jeppe Olsen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Kenneth Ruud
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø, N-9037 Tromsø, Norway
| |
Collapse
|
32
|
Recent Advances in the Coupled-Cluster Analytical Derivatives Theory for Molecules in Solution Described With the Polarizable Continuum Model (PCM). ADVANCES IN QUANTUM CHEMISTRY 2012. [DOI: 10.1016/b978-0-12-396498-4.00001-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
33
|
Severo Pereira Gomes A, Jacob CR. Quantum-chemical embedding methods for treating local electronic excitations in complex chemical systems. ACTA ACUST UNITED AC 2012. [DOI: 10.1039/c2pc90007f] [Citation(s) in RCA: 165] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
34
|
Sneskov K, Christiansen O. Excited state coupled cluster methods. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.99] [Citation(s) in RCA: 123] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
35
|
Caricato M, Scalmani G. On the Importance of the Orbital Relaxation in Ground-State Coupled Cluster Calculations in Solution with the Polarizable Continuum Model of Solvation. J Chem Theory Comput 2011; 7:4012-8. [PMID: 26598347 DOI: 10.1021/ct2006677] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- Marco Caricato
- Gaussian, Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Giovanni Scalmani
- Gaussian, Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| |
Collapse
|
36
|
Caricato M. CCSD-PCM: Improving upon the reference reaction field approximation at no cost. J Chem Phys 2011; 135:074113. [DOI: 10.1063/1.3624373] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
|
37
|
Caricato M, Scalmani G, Frisch MJ. Brueckner doubles coupled cluster method with the polarizable continuum model of solvation. J Chem Phys 2011; 134:244113. [DOI: 10.1063/1.3604560] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
|
38
|
Sneskov K, Schwabe T, Kongsted J, Christiansen O. The polarizable embedding coupled cluster method. J Chem Phys 2011; 134:104108. [DOI: 10.1063/1.3560034] [Citation(s) in RCA: 107] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
39
|
Manzoni V, Lyra ML, Gester RM, Coutinho K, Canuto S. Study of the optical and magnetic properties of pyrimidine in water combining PCM and QM/MM methodologies. Phys Chem Chem Phys 2010; 12:14023-33. [PMID: 20856965 DOI: 10.1039/c0cp00122h] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The solvent effects on the low-lying absorption spectrum and on the (15)N chemical shielding of pyrimidine in water are calculated using the combined and sequential Monte Carlo simulation and quantum mechanical calculations. Special attention is devoted to the solute polarization. This is included by an iterative procedure previously developed where the solute is electrostatically equilibrated with the solvent. In addition, we verify the simple yet unexplored alternative of combining the polarizable continuum model (PCM) and the hybrid QM/MM method. We use PCM to obtain the average solute polarization and include this in the MM part of the sequential QM/MM methodology, PCM-MM/QM. These procedures are compared and further used in the discrete and the explicit solvent models. The use of the PCM polarization implemented in the MM part seems to generate a very good description of the average solute polarization leading to very good results for the n-π* excitation energy and the (15)N nuclear chemical shield of pyrimidine in aqueous environment. The best results obtained here using the solute pyrimidine surrounded by 28 explicit water molecules embedded in the electrostatic field of the remaining 472 molecules give the statistically converged values for the low lying n-π* absorption transition in water of 36 900 ± 100 (PCM polarization) and 36 950 ± 100 cm(-1) (iterative polarization), in excellent agreement among one another and with the experimental value observed with a band maximum at 36 900 cm(-1). For the nuclear shielding (15)N the corresponding gas-water chemical shift obtained using the solute pyrimidine surrounded by 9 explicit water molecules embedded in the electrostatic field of the remaining 491 molecules give the statistically converged values of 24.4 ± 0.8 and 28.5 ± 0.8 ppm, compared with the inferred experimental value of 19 ± 2 ppm. Considering the simplicity of the PCM over the iterative polarization this is an important aspect and the computational savings point to the possibility of dealing with larger solute molecules. This PCM-MM/QM approach reconciles the simplicity of the PCM model with the reliability of the combined QM/MM approaches.
Collapse
Affiliation(s)
- Vinícius Manzoni
- Instituto de Física, Universidade Federal de Alagoas, 57072-970, Maceió, AL, Brazil
| | | | | | | | | |
Collapse
|
40
|
Cammi R, Fukuda R, Ehara M, Nakatsuji H. Symmetry-adapted cluster and symmetry-adapted cluster-configuration interaction method in the polarizable continuum model: Theory of the solvent effect on the electronic excitation of molecules in solution. J Chem Phys 2010; 133:024104. [DOI: 10.1063/1.3456540] [Citation(s) in RCA: 71] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
|
41
|
Caricato M, Mennucci B, Scalmani G, Trucks GW, Frisch MJ. Electronic excitation energies in solution at equation of motion CCSD level within a state specific polarizable continuum model approach. J Chem Phys 2010; 132:084102. [PMID: 20192285 DOI: 10.1063/1.3314221] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
We present a study of excitation energies in solution at the equation of motion coupled cluster singles and doubles (EOM-CCSD) level of theory. The solvent effect is introduced with a state specific polarizable continuum model (PCM), where the solute-solvent interaction is specific for the state of interest. Three definitions of the excited state one-particle density matrix (1PDM) are tested in order to gain information for the development of an integrated EOM-CCSD/PCM method. The calculations show the accuracy of this approach for the computation of such property in solution. Solvent shifts between nonpolar and polar solvents are in good agreement with experiment for the test cases. The completely unrelaxed 1PDM is shown to be a balanced choice between computational effort and accuracy for vertical excitation energies, whereas the response of the ground state CCSD amplitudes and of the molecular orbitals is important for other properties, as for instance the dipole moment.
Collapse
Affiliation(s)
- Marco Caricato
- Department of Chemistry, Yale University, 225 Prospect St., New Haven, Connecticut 06511, USA.
| | | | | | | | | |
Collapse
|
42
|
Sneskov K, Matito E, Kongsted J, Christiansen O. Approximate Inclusion of Triple Excitations in Combined Coupled Cluster/Molecular Mechanics: Calculations of Electronic Excitation Energies in Solution for Acrolein, Water, Formamide, and N-Methylacetamide. J Chem Theory Comput 2010; 6:839-50. [DOI: 10.1021/ct900641w] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Kristian Sneskov
- The Lundbeck Foundation Center for Theoretical Chemistry, Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark, and Department of Physics and Chemistry, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Eduard Matito
- The Lundbeck Foundation Center for Theoretical Chemistry, Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark, and Department of Physics and Chemistry, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Jacob Kongsted
- The Lundbeck Foundation Center for Theoretical Chemistry, Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark, and Department of Physics and Chemistry, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Ove Christiansen
- The Lundbeck Foundation Center for Theoretical Chemistry, Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark, and Department of Physics and Chemistry, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| |
Collapse
|
43
|
Cammi R. Quantum cluster theory for the polarizable continuum model. I. The CCSD level with analytical first and second derivatives. J Chem Phys 2009; 131:164104. [DOI: 10.1063/1.3245400] [Citation(s) in RCA: 73] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Møgelhøj A, Aidas K, Mikkelsen KV, Sauer SPA, Kongsted J. Prediction of spin-spin coupling constants in solution based on combined density functional theory/molecular mechanics. J Chem Phys 2009; 130:134508. [PMID: 19355752 DOI: 10.1063/1.3098255] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present theory and implementation of calculation of spin-spin coupling constants within combined quantum mechanics/molecular mechanics methods. Special attention is given to the role of explicit solvent polarization as well as the molecular consequences due to hydrogen bonding. The model is generally applicable but is here implemented for the case of density functional theory. First applications to liquid water and acetylene in aqueous solution are presented. Good agreement between theory and experiment is obtained in both cases, thereby showing the strength of our approach. Finally, spin-spin coupling constants across hydrogen bonds are discussed considering for the first time the role of an explicit solvent on this class of spin-spin couplings.
Collapse
Affiliation(s)
- Andreas Møgelhøj
- Department of Chemistry, H. C. Ørsted Institute, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| | | | | | | | | |
Collapse
|
45
|
From Molecules to Droplets. ADVANCES IN QUANTUM CHEMISTRY 2008. [DOI: 10.1016/s0065-3276(07)00217-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register]
|
46
|
Georg HC, Coutinho K, Canuto S. Converged electronic polarization of acetone in liquid water and the role in the n–π∗ transition. Chem Phys Lett 2006. [DOI: 10.1016/j.cplett.2006.08.047] [Citation(s) in RCA: 77] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
47
|
Improta R, Barone V, Scalmani G, Frisch MJ. A state-specific polarizable continuum model time dependent density functional theory method for excited state calculations in solution. J Chem Phys 2006; 125:054103. [PMID: 16942199 DOI: 10.1063/1.2222364] [Citation(s) in RCA: 554] [Impact Index Per Article: 30.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An effective state specific (SS) model for the inclusion of solvent effects in time dependent density functional theory (TD-DFT) computations of excited electronic states has been developed and coded in the framework of the so-called polarizable continuum model (PCM). Different relaxation time regimes can be treated thus giving access to a number of different spectroscopic properties together with solvent relaxation energies of paramount relevance in electron transfer processes. SS and conventional linear response (LR) models have been compared for two benchmark systems (coumarin 153 and formaldehyde in different solvents) and in the limiting simple case of a dipolar solute embedded in a spherical cavity. The results point out the complementarity of LR and SS approaches and the advantages of the latter model especially for polar solvents. The favorable scaling properties of PCM-TD-DFT models in both SS and LR variants and their availability in effective quantum mechanical codes pave the route for the computation of reliable spectroscopic properties of large molecules of technological and/or biological interest in their natural environments.
Collapse
Affiliation(s)
- Roberto Improta
- Dipartimento di Chimica, Università Federico II, Complesso Monte S. Angelo, via Cintia, I-80126 Napoli, Italy.
| | | | | | | |
Collapse
|
48
|
Lee JY, Yoshida N, Hirata F. Conformational Equilibrium of 1,2-Dichloroethane in Water: Comparison of PCM and RISM-SCF Methods. J Phys Chem B 2006; 110:16018-25. [PMID: 16898759 DOI: 10.1021/jp0606762] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The RISM-SCF and polarizable continuum model (PCM) approaches have been applied to study the conformational equilibrium of 1,2-dichloroethane (DCE) in water. Both the electron correlation effect and basis sets play an important role in the relative energies of the gauche and trans conformers in gas and solution phases. Both PCM and RISM-MP2 methods resulted in a consistent trend with the previous experimental and theoretical studies that the population of the gauche conformer increases in going from the gas phase to the aqueous solution. However, the PCM treatment could not describe the solvent effect completely in that the sign of the relative free energy of the gauche and trans forms is opposite to the most recent experimental and theoretical data, while the RISM-MP2 gives the right sign in the free energy difference. We found that the larger excess chemical potential gain (by ca. -4.1 kcal/mol) for the gauche conformer is large enough to result in the gauche preference of DCE in water, though it has to compensate for more solute reorganization energy (approximately 1.6 kcal/mol) and overcome the energy difference (approximately 1.6 kcal/mol) in the gas phase. The radial distribution functions between DCE and the nearest water shows that the electrostatic repulsion between chlorine and oxygen atoms is higher in the trans conformer than in the gauche one, while the attractive interaction between chlorine and hydrogen of water is higher in the gauche conformer.
Collapse
Affiliation(s)
- Jin Yong Lee
- Department of Theoretical Molecular Science, Institute for Molecular Science, Okazaki 444-8585, Japan
| | | | | |
Collapse
|
49
|
Kongsted J, Osted A, Mikkelsen KV, Christiansen O. Second harmonic generation second hyperpolarizability of water calculated using the combined coupled cluster dielectric continuum or different molecular mechanics methods. J Chem Phys 2006; 120:3787-98. [PMID: 15268543 DOI: 10.1063/1.1642593] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this article we report the first calculations of second harmonic generation second hyperpolarizability of liquid water using coupled cluster/molecular mechanics (CC/MM) methods or coupled cluster/dielectric continuum (CC/DC) methods. The latter approach treats the solvent as an isotropic homogeneous fluid while the former accounts for the discrete nature of the solvent molecules. The CC/MM approach may include or exclude polarization effects explicitly. Alternatively, polarization effects may be included using perturbation theory. The CC descriptions implemented are the coupled cluster second-order approximate singles and doubles (CC2) and coupled cluster singles and doubles models. The second harmonic generation second hyperpolarizabilities are, depending on the model, obtained using either an analytical implementation of the cubic response function or using an analytical implementation of the quadratic response function combined with the finite field technique. The CC/MM results for the second harmonic generation second hyperpolarizability compare excellently with experimental data while a significant overestimation is found when using the CC/DC model. Particular, the cavity radius in the CC/DC calculations have an enormous effects on this fourth-order property.
Collapse
Affiliation(s)
- Jacob Kongsted
- Department of Chemistry, HC Ørsted Institute, University of Copenhagen, Universitetsparken 5, DK-2100 Ø, Denmark.
| | | | | | | |
Collapse
|
50
|
Mikkelsen KV. CORRELATED ELECTRONIC STRUCTURE NONLINEAR RESPONSE METHODS FOR STRUCTURED ENVIRONMENTS. Annu Rev Phys Chem 2006; 57:365-402. [PMID: 16599815 DOI: 10.1146/annurev.physchem.57.032905.104740] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
▪ Abstract This contribution concerns a brief outline of structural environment models where correlated electronic structure response methods are utilized for the determination of nonlinear optical properties of molecules. The presentation provides theory and applications of a heterogeneous dielectric media model and a quantum mechanical–classical mechanical model at the level of correlated electronic structure response methods. The correlated electronic structure response methods include the multiconfigurational self-consistent field method.
Collapse
Affiliation(s)
- Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen, Denmark.
| |
Collapse
|