1
|
Kumawat RL, Sherrill CD. High-Order Quantum-Mechanical Analysis of Hydrogen Bonding in Hachimoji and Natural DNA Base Pairs. J Chem Inf Model 2023; 63:3150-3157. [PMID: 37125692 DOI: 10.1021/acs.jcim.3c00428] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
High-order quantum chemistry is applied to hydrogen-bonded natural DNA nucleobase pairs [adenine:thymine (A:T) and guanine:cytosine (G:C)] and non-natural Hachimoji nucleobase pairs [isoguanine:1-methylcytosine (B:S) and 2-aminoimidazo[1,2a][1,3,5]triazin-4(1H)-one:6-amino-5-nitropyridin-2-one (P:Z)] to see how the intermolecular interaction energies and their energetic components (electrostatics, exchange-repulsion, induction/polarization, and London dispersion interactions) vary among the base pairs. We examined the Hoogsteen (HG) geometries in addition to the traditional Watson-Crick (WC) geometries. Coupled-cluster theory through perturbative triples [CCSD(T)] extrapolated to the complete basis set (CBS) limit and high-order symmetry-adapted perturbation theory (SAPT) at the SAPT2+(3)(CCD)δMP2/aug-cc-pVTZ level are used to estimate highly accurate noncovalent interaction energies. Electrostatic interactions are the most attractive component of the interaction energies, but the sum of induction/polarization and London dispersion is nearly as large, for all base pairs and geometries considered. Interestingly, the non-natural Hachimoji base pairs interact more strongly than the corresponding natural base pairs, by -21.8 (B:S) and -0.3 (P:Z) kcal mol-1 in the WC geometries, according to CCSD(T)/CBS. This is consistent with the H-bond distances being generally shorter in the non-natural base pairs. The natural base pairs are energetically more stabilized in their Hoogsteen geometries than in their WC geometries. The Hoogsteen geometry makes the A:T base pair slightly more stable, by -0.8 kcal mol-1, and it greatly stabilizes the G:C+ base pair, by -15.3 kcal mol-1. The G:C+ stabilization is mainly due to the fact that C has typically added a proton when found in Hoogsteen geometries. By contrast, Hoogsteen geometries are substantially less favorable than WC geometries for non-natural Hachimoji base pairs, by 17.3 (B:S) and 13.8 (P:Z) kcal mol-1.
Collapse
Affiliation(s)
- Rameshwar L Kumawat
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
2
|
Volk KR, Casabianca LB. Quantum mechanical study of interactions between sunscreen ingredients and nucleotide bases. J Mol Model 2022; 28:243. [PMID: 35925497 DOI: 10.1007/s00894-022-05253-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2022] [Accepted: 07/27/2022] [Indexed: 11/26/2022]
Abstract
Interactions between the popular sunscreen ingredients oxybenzone and homosalate and DNA bases have been studied using density functional theory and ab initio methods. Low-energy structures for each sunscreen ingredient interacting with each nucleotide base in either a pi-stacked or hydrogen-bonded fashion were found. The binding energies are comparable to those for the Watson-Crick-Franklin Ade-Thy and Cyt-Gua pairs. Pi-stacked and hydrogen-bonded structures are comparable in energy, with hydrogen-bonded structures having a more negative counterpoise-corrected binding energy, while the final pi-stacked structures are lower in energy. This is due to a geometrical rearrangement required to form the hydrogen bonds that raise the total energy of the complex. It was also found that when using the M06-2X density functional, the STO-3G basis set favors hydrogen bonding, but 6-31G(d) and 6-31 + G(s) basis sets predict similar binding geometries.
Collapse
Affiliation(s)
- Kyle R Volk
- Department of Chemistry, Clemson University, Clemson, SC, USA
| | | |
Collapse
|
3
|
|
4
|
Mammadova F, Hamarat B, Ahmadli D, Şahin O, Bozkaya U, Türkmen YE. Polarization‐Enhanced Hydrogen Bonding in 1,8‐Dihydroxynaphthalene: Conformational Analysis, Binding Studies and Hydrogen Bonding Catalysis. ChemistrySelect 2020. [DOI: 10.1002/slct.202002960] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Flora Mammadova
- Department of Chemistry, Faculty of Science Bilkent University Ankara 06800 Turkey
| | - Büşra Hamarat
- Department of Chemistry Hacettepe University Ankara 06800 Turkey
| | - Dilgam Ahmadli
- Department of Chemistry, Faculty of Science Bilkent University Ankara 06800 Turkey
| | - Onur Şahin
- Scientific and Technological Research Application and Research Center Sinop University Sinop 57000 Turkey
| | - Uğur Bozkaya
- Department of Chemistry Hacettepe University Ankara 06800 Turkey
| | - Yunus E. Türkmen
- Department of Chemistry, Faculty of Science Bilkent University Ankara 06800 Turkey
- UNAM-National Nanotechnology Research Center, Institute of Materials Science and Nanotechnology Bilkent University Ankara 06800 Turkey
| |
Collapse
|
5
|
Fedorov IA. Elastic properties of the molecular crystals of hydrocarbons from first principles calculations. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:085704. [PMID: 31698348 DOI: 10.1088/1361-648x/ab554e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/19/2023]
Abstract
I studied the elastic properties of oligoacenes and polycyclic aromatic hydrocarbons crystals within the framework of density functional theory with van der Waals interactions. The full sets of elastic constants were computed. The computed parameters have good agreement with experimental data. The study of four forms (α, β, γ and dimeric) demonstrate that the stacking of molecules leads to different character of interaction between molecules. Crystals with dimeric form have lower anisotropy. If symmetry of crystals is the same, the increase in number of aromatic rings results in anisotropy increase.
Collapse
|
6
|
Giovannini T, Lafiosca P, Cappelli C. A General Route to Include Pauli Repulsion and Quantum Dispersion Effects in QM/MM Approaches. J Chem Theory Comput 2017; 13:4854-4870. [PMID: 28898079 DOI: 10.1021/acs.jctc.7b00776] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
A methodology to account for nonelectrostatic interactions in Quantum Mechanical (QM)/Molecular Mechanics (MM) approaches is developed. Formulations for Pauli repulsion and dispersion energy, explicitly depending on the QM density, are derived. Such expressions are based on the definition of an auxiliary density on the MM portion and the Tkatchenko-Scheffler (TS) approach, respectively. The developed method is general enough to be applied to any QM/MM method and partition, provided an accurate tuning of a small number of parameters is obtained. The coupling of the method with both nonpolarizable and fully polarizable QM/fluctuating charge (FQ) approaches is reported and applied. A suitable parametrization for the aqueous solution, so that its most representative features are well reproduced, is outlined. Then, the obtained parametrization and method are applied to calculate the nonelectrostatic (repulsion and dispersion) interaction energy of nicotine in aqueous solution.
Collapse
Affiliation(s)
| | - Piero Lafiosca
- Scuola Normale Superiore , Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore , Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
7
|
Fortenberry RC, Thackston R, Francisco JS, Lee TJ. Toward the laboratory identification of the not-so-simple NS2neutral and anion isomers. J Chem Phys 2017; 147:074303. [DOI: 10.1063/1.4985901] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ryan C. Fortenberry
- Department of Chemistry and Biochemistry, Georgia Southern University, Statesboro, Georgia 30460-8064, USA
| | - Russell Thackston
- Department of Information Technology, Georgia Southern University, Statesboro, Georgia 30460-8150, USA
| | - Joseph S. Francisco
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| | - Timothy J. Lee
- MS 245-1, NASA Ames Research Center, Moffett Field, California 94035-1000, USA
| |
Collapse
|
8
|
Shen LQ, Kundu S, Collins TJ, Bominaar EL. Analysis of Hydrogen Atom Abstraction from Ethylbenzene by an FeVO(TAML) Complex. Inorg Chem 2017; 56:4347-4356. [DOI: 10.1021/acs.inorgchem.6b02796] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Longzhu Q. Shen
- Department of Chemistry, Carnegie Mellon University, 4400 Fifth Avenue, Pittsburgh, Pennsylvania 15213, United States
| | - Soumen Kundu
- Department of Chemistry, Carnegie Mellon University, 4400 Fifth Avenue, Pittsburgh, Pennsylvania 15213, United States
| | - Terrence J. Collins
- Department of Chemistry, Carnegie Mellon University, 4400 Fifth Avenue, Pittsburgh, Pennsylvania 15213, United States
| | - Emile L. Bominaar
- Department of Chemistry, Carnegie Mellon University, 4400 Fifth Avenue, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
9
|
Filipek G, Fortenberry RC. Formation of Potential Interstellar Noble Gas Molecules in Gas and Adsorbed Phases. ACS OMEGA 2016; 1:765-772. [PMID: 31457160 PMCID: PMC6640802 DOI: 10.1021/acsomega.6b00249] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/16/2016] [Accepted: 10/20/2016] [Indexed: 06/10/2023]
Abstract
The discovery of naturally occurring ArH+ in various regions of the interstellar medium has shown the need for more understanding of the reactions that lead to covalently bonded noble gas molecules. The test comes with trying to predict the formation of other small noble gas molecules. Many molecules have been observed in various interstellar environments, which possess the possibility of bonding with noble gases. This work explores how both argon and neon can form bonds to ligands made of these species through quantum chemical computations. Argon and neon are chosen as they are among the most abundant atoms in the universe but are more polarizable than the more common but smaller helium atom. Reactions leading to noble gas molecules are modeled in the gas phase as well as through the adsorbed phase by catalysis with a polycyclic aromatic hydrocarbon (PAH) surface. The adsorption energy of the neutral noble gas atoms to the surface increases as the size of the PAH also increases but this is still less than 10 kcal/mol. It is proposed and supported herein that an incoming molecule can bond with the noble gas atom adsorbed onto the PAH, form a stable structure, and allow the PAH to function as the leaving group. This work shows that the noble gas molecules ArCCH+, ArOH+, ArNH+, and NeCCH+ are not only stable minima on their respective potential energy surfaces but also can be formed in either the gas phase or through PAH adsorption with known or hypothesized interstellar molecules. Most notably, NeCCH+ does not appear to form in the gas phase but could be catalyzed on PAH surfaces. Hence, the interstellar detection of such molecules could also serve as a probe for the observation of interstellar PAHs.
Collapse
|
10
|
Forni A, Pieraccini S, Franchini D, Sironi M. Assessment of DFT Functionals for QTAIM Topological Analysis of Halogen Bonds with Benzene. J Phys Chem A 2016; 120:9071-9080. [PMID: 27718571 DOI: 10.1021/acs.jpca.6b07578] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Halogen bonding, a noncovalent interaction between a halogen atom and a nucleophilic site, is receiving a growing attention in the chemical community stimulating a large number of theoretical investigations. The density functional theory (DFT) approach revealed to be one of the most suitable methods owing to its accuracy and low computational cost. We report here a detailed analysis of the performance of an extensive set of DFT functionals in reproducing accurate binding energies and topological properties for the halogen-bonding interaction of either NCX or PhX molecules (X = F, Cl, Br, I) with the aromatic system of benzene in the T-shaped configuration. It was found that the better performance for both sets of properties is provided by a small subset of functionals able to take into account, implicitly or explicitly (by inclusion of an additive pairwise potential), the dispersion contribution, that is, ωB97X, M06-2X, M11, mPW2PLYP-D, and B2PLYP-D3.
Collapse
Affiliation(s)
- Alessandra Forni
- Istituto di Scienze e Tecnologie Molecolari, CNR (CNR-ISTM) and INSTM UdR , Via Golgi 19, 20133 Milano, Italy
| | - Stefano Pieraccini
- Istituto di Scienze e Tecnologie Molecolari, CNR (CNR-ISTM) and INSTM UdR , Via Golgi 19, 20133 Milano, Italy.,Dipartimento di Chimica and INSTM UdR, Università degli Studi di Milano , Via Golgi 19, 20133 Milano, Italy
| | - Davide Franchini
- Dipartimento di Chimica and INSTM UdR, Università degli Studi di Milano , Via Golgi 19, 20133 Milano, Italy
| | - Maurizio Sironi
- Istituto di Scienze e Tecnologie Molecolari, CNR (CNR-ISTM) and INSTM UdR , Via Golgi 19, 20133 Milano, Italy.,Dipartimento di Chimica and INSTM UdR, Università degli Studi di Milano , Via Golgi 19, 20133 Milano, Italy
| |
Collapse
|
11
|
Ghosh D. Hybrid Equation-of-Motion Coupled-Cluster/Effective Fragment Potential Method: A Route toward Understanding Photoprocesses in the Condensed Phase. J Phys Chem A 2016; 121:741-752. [DOI: 10.1021/acs.jpca.6b08263] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Debashree Ghosh
- Physical
and Materials Chemistry
Division, CSIR-National Chemical Laboratory, Pune, India 411008
| |
Collapse
|
12
|
Fortenberry RC, Moore MM, Lee TJ. Excited State Trends in Bidirectionally Expanded Closed-Shell PAH and PANH Anions. J Phys Chem A 2016; 120:7327-34. [PMID: 27585793 DOI: 10.1021/acs.jpca.6b06654] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Some anions are known to exhibit excited states independent of external forces such as dipole moments and induced polarizabilities. Such states exist simply as a result of the stabilization of valence accepting orbitals whereby the binding energy of the extra electron is greater than the valence excitation energy. Closed-shell anions are interesting candidates for such transitions since their ground-state, spin-paired nature makes the anions more stable from the beginning. Consequently, this work shows the point beyond which deprotonated, closed-shell polycyclic aromatic hydrocarbons (PAHs) and those PAHs containing nitrogen heteroatoms (PANHs) will exhibit valence excited states. This behavior has already been demonstrated in some PANHs and for anistropically extended PAHs. This work establishes a general trend for PAHs/PANHs of arbitrary size and directional extension, whether in one dimension or two. Once seven six-membered rings make up a PAH/PANH, valence excited states are present. For most classes of PAHs/PANHs, this number is closer to four. Even though most of these excited states are weak absorbers, the sheer number of PAHs present in various astronomical environments should make them significant contributors to astronomical spectra.
Collapse
Affiliation(s)
- Ryan C Fortenberry
- Department of Chemistry, Georgia Southern University , Statesboro, Georgia 30460, United States
| | - Megan M Moore
- Department of Chemistry, Georgia Southern University , Statesboro, Georgia 30460, United States
| | - Timothy J Lee
- NASA Ames Research Center , MS 245-1, Moffett Field, California 94035-1000, United States
| |
Collapse
|
13
|
Stout HD, Kleespies ST, Chiang CW, Lee WZ, Que L, Münck E, Bominaar EL. Spectroscopic and Theoretical Study of Spin-Dependent Electron Transfer in an Iron(III) Superoxo Complex. Inorg Chem 2016; 55:5215-26. [PMID: 27159412 DOI: 10.1021/acs.inorgchem.6b00134] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
It was shown previously (J. Am. Chem. Soc. 2014, 136, 10846) that bubbling of O2 into a solution of Fe(II)(BDPP) (H2BDPP = 2,6-bis[[(S)-2-(diphenylhydroxymethyl)-1-pyrrolidinyl]methyl]pyridine) in tetrahydrofuran at -80 °C generates a high-spin (SFe = (5)/2) iron(III) superoxo adduct, 1. Mössbauer studies revealed that 1 is an exchange-coupled system, [Formula: see text], where SR = (1)/2 is the spin of the superoxo radical, of which the spectra were not well enough resolved to determine whether the coupling was ferromagnetic (S = 3 ground state) or antiferromagnetic (S = 2). The glass-forming 2-methyltetrahydrofuran solvent yields highly resolved Mössbauer spectra from which the following data have been extracted: (i) the ground state of 1 has S = 3 (J < 0); (ii) |J| > 15 cm(-1); (iii) the zero-field-splitting parameters are D = -1.1 cm(-1) and E/D = 0.02; (iv) the major component of the electric-field-gradient tensor is tilted ≈7° relative to the easy axis of magnetization determined by the MS = ±3 and ±2 doublets. The excited-state MS = ±2 doublet yields a narrow parallel-mode electron paramagnetic resonance signal at g = 8.03, which was used to probe the magnetic hyperfine splitting of (17)O-enriched O2. A theoretical model that considers spin-dependent electron transfer for the cases where the doubly occupied π* orbital of the superoxo ligand is either "in" or "out" of the plane defined by the bent Fe-OO moiety correctly predicts that 1 has an S = 3 ground state, in contrast to the density functional theory calculations for 1, which give a ground state with both the wrong spin and orbital configuration. This failure has been traced to a basis set superposition error in the interactions between the superoxo moiety and the adjacent five-membered rings of the BDPP ligand and signals a fundamental problem in the quantum chemistry of O2 activation.
Collapse
Affiliation(s)
- Heather D Stout
- Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| | - Scott T Kleespies
- Department of Chemistry and Center for Metals in Biocatalysis, University of Minnesota , Minneapolis, Minnesota 55455, United States
| | - Chien-Wei Chiang
- Department of Chemistry, National Taiwan Normal University , 88, Section 4, Ting-Chow Road, Taipei 11677, Taiwan (R.O.C.)
| | - Way-Zen Lee
- Department of Chemistry, National Taiwan Normal University , 88, Section 4, Ting-Chow Road, Taipei 11677, Taiwan (R.O.C.)
| | - Lawrence Que
- Department of Chemistry and Center for Metals in Biocatalysis, University of Minnesota , Minneapolis, Minnesota 55455, United States
| | - Eckard Münck
- Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| | - Emile L Bominaar
- Department of Chemistry, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
14
|
Trouillas P, Sancho-García JC, De Freitas V, Gierschner J, Otyepka M, Dangles O. Stabilizing and Modulating Color by Copigmentation: Insights from Theory and Experiment. Chem Rev 2016; 116:4937-82. [PMID: 26959943 DOI: 10.1021/acs.chemrev.5b00507] [Citation(s) in RCA: 286] [Impact Index Per Article: 35.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Natural anthocyanin pigments/dyes and phenolic copigments/co-dyes form noncovalent complexes, which stabilize and modulate (in particular blue, violet, and red) colors in flowers, berries, and food products derived from them (including wines, jams, purees, and syrups). This noncovalent association and their electronic and optical implications constitute the copigmentation phenomenon. Over the past decade, experimental and theoretical studies have enabled a molecular understanding of copigmentation. This review revisits this phenomenon to provide a comprehensive description of the nature of binding (the dispersion and electrostatic components of π-π stacking, the hydrophobic effect, and possible hydrogen-bonding between pigment and copigment) and of spectral modifications occurring in copigmentation complexes, in which charge transfer plays an important role. Particular attention is paid to applications of copigmentation in food chemistry.
Collapse
Affiliation(s)
- Patrick Trouillas
- INSERM UMR 850, Univ. Limoges , Faculty of Pharmacy, 2 rue du Dr. Marcland, F-87025 Limoges, France.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc , tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Juan C Sancho-García
- Departamento de Química Física, Universidad de Alicante , Apartado de Correos 99, E-03080 Alicante, Spain
| | - Victor De Freitas
- REQUIMTE/LAQV - Research Unit, Faculty of Science, Porto University , Rua do Campo Alegre, 4169-007 Porto, Portugal
| | - Johannes Gierschner
- Madrid Institute for Advanced Studies - IMDEA Nanoscience , C/Faraday 9, Ciudad Universitaria de Cantoblanco, E-28049 Madrid, Spain
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc , tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Olivier Dangles
- University of Avignon, INRA, UMR408 SQPOV , F-84000 Avignon, France
| |
Collapse
|
15
|
Simmonett AC, Pickard FC, Schaefer HF, Brooks BR. An efficient algorithm for multipole energies and derivatives based on spherical harmonics and extensions to particle mesh Ewald. J Chem Phys 2015; 140:184101. [PMID: 24832247 DOI: 10.1063/1.4873920] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
Next-generation molecular force fields deliver accurate descriptions of non-covalent interactions by employing more elaborate functional forms than their predecessors. Much work has been dedicated to improving the description of the electrostatic potential (ESP) generated by these force fields. A common approach to improving the ESP is by augmenting the point charges on each center with higher-order multipole moments. The resulting anisotropy greatly improves the directionality of the non-covalent bonding, with a concomitant increase in computational cost. In this work, we develop an efficient strategy for enumerating multipole interactions, by casting an efficient spherical harmonic based approach within a particle mesh Ewald (PME) framework. Although the derivation involves lengthy algebra, the final expressions are relatively compact, yielding an approach that can efficiently handle both finite and periodic systems without imposing any approximations beyond PME. Forces and torques are readily obtained, making our method well suited to modern molecular dynamics simulations.
Collapse
Affiliation(s)
- Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| |
Collapse
|
16
|
Calbo J, Ortí E, Sancho-García JC, Aragó J. The Nonlocal Correlation Density Functional VV10. ACTA ACUST UNITED AC 2015. [DOI: 10.1016/bs.arcc.2015.09.002] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
|
17
|
Bozkaya U, Sherrill CD. Orbital-optimized MP2.5 and its analytic gradients: Approaching CCSD(T) quality for noncovalent interactions. J Chem Phys 2014; 141:204105. [DOI: 10.1063/1.4902226] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Atatürk University, Erzurum 25240, Turkey
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
18
|
Esrafili MD, Yourdkhani S, Bahrami A. Characteristics and nature of the halogen-bonding interactions between CCl3F and ozone: a supermolecular and SAPT study. Mol Phys 2013. [DOI: 10.1080/00268976.2013.788740] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
19
|
Howard JC, Tschumper GS. Wavefunction methods for the accurate characterization of water clusters. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013. [DOI: 10.1002/wcms.1168] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
| | - Gregory S. Tschumper
- Department of Chemistry and Biochemistry University of Mississippi, University Mississippi USA
| |
Collapse
|
20
|
Raju RK, Bloom JWG, Wheeler SE. Broad Transferability of Substituent Effects in π-Stacking Interactions Provides New Insights into Their Origin. J Chem Theory Comput 2013; 9:3479-90. [DOI: 10.1021/ct400481r] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Rajesh K. Raju
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
| | - Jacob W. G. Bloom
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
| | - Steven E. Wheeler
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
| |
Collapse
|
21
|
Wheeler SE. Understanding substituent effects in noncovalent interactions involving aromatic rings. Acc Chem Res 2013; 46:1029-38. [PMID: 22725832 DOI: 10.1021/ar300109n] [Citation(s) in RCA: 359] [Impact Index Per Article: 32.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Noncovalent interactions involving aromatic rings such as π-stacking, cation/π, and anion/π interactions are central to many areas of modern chemistry. Decades of experimental studies have provided key insights into the impact of substituents on these interactions, leading to the development of simple intuitive models. However, gas-phase computational studies have raised some doubts about the physical underpinnings of these widespread models. In this Account we review our recent efforts to unravel the origin of substituent effects in π-stacking and ion/π interactions through computational studies of model noncovalent dimers. First, however, we dispel the notion that so-called aromatic interactions depend on the aromaticity of the interacting rings by studying model π-stacked dimers in which the aromaticity of one of the monomers can be "switched off". Somewhat surprisingly, the results show that not only is aromaticity unnecessary for π-stacking interactions, but it actually hinders these interactions to some extent. Consequently, when thinking about π-stacking interactions, researchers should consider broader classes of planar molecules, not just aromatic systems. Conventional models maintain that substituent effects in π-stacking interactions result from changes in the aryl π-system. This view suggests that π-stacking interactions are maximized when one ring is substituted with electron-withdrawing groups and the other with electron donors. In contrast to these prevailing models, we have shown that substituent effects in π-stacking interactions can be described in terms of direct, local interactions between the substituents and the nearby vertex of the other arene. As a result, in polysubstituted π-stacked dimers the substituents operate independently unless they are in each other's local environment. This means that in π-stacked dimers in which one arene is substituted with electron donors and the other with electron acceptors the interactions will be enhanced only to the extent provided by each substituent on its own, unless the substituents on opposing rings are in close proximity. Overall, this local, direct interaction model predicts that substituent effects in π-stacking interactions will be additive and transferable and will also depend on the relative position of substituents on opposing rings. For cation/π and anion/π interactions, similar π-resonance-based models pervade the literature. Again, computational results indicate that substituent effects in model ion/π complexes can be described primarily in terms of direct interactions between the ion and the substituent. Changes in the aryl π-system do not significantly affect these interactions. We also present a simple electrostatic model that further demonstrates this effect and suggests that the dominant interaction for simple substituents is the interaction of the charged ion with the local dipole associated with the substituents. Finally, we discuss substituent effects in electrostatic potentials (ESPs), which are widely used in discussions of noncovalent interactions. In the past, widespread misconceptions have confused the relationship between changes in ESPs and local changes in the electron density. We have shown that computed ESP plots of diverse substituted arenes can be reproduced without altering the aryl π-density. This is because substituent-induced changes in the ESP above the center of aryl rings result primarily from through-space effects of substituents rather than through changes in the distribution of the π-electron density.
Collapse
Affiliation(s)
- Steven E. Wheeler
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
| |
Collapse
|
22
|
Fedorov IA, Zhuravlev YN, Berveno VP. Structural and electronic properties of perylene from first principles calculations. J Chem Phys 2013; 138:094509. [DOI: 10.1063/1.4794046] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
|
23
|
Bozkaya U, Schaefer HF. Symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles method: improving upon CCSD(T) and CCSD(T)(Λ): preliminary application. J Chem Phys 2012; 136:204114. [PMID: 22667547 DOI: 10.1063/1.4720382] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles (OO-CCD or simply "OD" for short) method are investigated. The conventional symmetric and asymmetric perturbative triples corrections [(T) and (T)(Λ)] are implemented, the latter one for the first time. Additionally, two new triples corrections, denoted as OD(Λ) and OD(Λ)(T), are introduced. We applied the new methods to potential energy surfaces of the BH, HF, C(2), N(2), and CH(4) molecules, and compare the errors in total energies, with respect to full configuration interaction, with those from the standard coupled-cluster singles and doubles (CCSD), with perturbative triples [CCSD(T)], and asymmetric triples correction (CCSD(T)(Λ)) methods. The CCSD(T) method fails badly at stretched geometries, the corresponding nonparallelity error is 7-281 kcal mol(-1), although it gives reliable results near equilibrium geometries. The new symmetric triples correction, CCSD(Λ), noticeably improves upon CCSD(T) (by 4-14 kcal mol(-1)) for BH, HF, and CH(4); however, its performance is worse than CCSD(T) (by 1.6-4.2 kcal mol(-1)) for C(2) and N(2). The asymmetric triples corrections, CCSD(T)(Λ) and CCSD(Λ)(T), perform remarkably better than CCSD(T) (by 5-18 kcal mol(-1)) for the BH, HF, and CH(4) molecules, while for C(2) and N(2) their results are similar to those of CCSD(T). Although the performance of CCSD and OD is similar, the situation is significantly different in the case of triples corrections, especially at stretched geometries. The OD(T) method improves upon CCSD(T) by 1-279 kcal mol(-1). The new symmetric triples correction, OD(Λ), enhances the OD(T) results (by 0.01-2.0 kcal mol(-1)) for BH, HF, and CH(4); however, its performance is worse than OD(T) (by 1.9-2.3 kcal mol(-1)) for C(2) and N(2). The asymmetric triples corrections, OD(T)(Λ) and OD(Λ)(T), perform better than OD(T) (by 2.0-6.2 kcal mol(-1)). The latter method is slightly better for the BH, HF, and CH(4) molecules. However, for C(2) and N(2) the new results are similar to those of OD(T). For the BH, HF, and CH(4) molecules, OD(Λ)(T) provides the best potential energy curves among the considered methods, while for C(2) and N(2) the OD(T) method prevails. Hence, for single-bond breaking the OD(Λ)(T) method appears to be superior, whereas for multiple-bond breaking the OD(T) method is better.
Collapse
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Atatürk University, Erzurum 25240, Turkey.
| | | |
Collapse
|
24
|
Effect of stepwise microhydration on the methylammonium···phenol and ammonium···phenol interaction. J Mol Model 2012; 19:1985-94. [DOI: 10.1007/s00894-012-1579-9] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2012] [Accepted: 08/21/2012] [Indexed: 10/27/2022]
|
25
|
Grimme S. Supramolecular binding thermodynamics by dispersion-corrected density functional theory. Chemistry 2012; 18:9955-64. [PMID: 22782805 DOI: 10.1002/chem.201200497] [Citation(s) in RCA: 1183] [Impact Index Per Article: 98.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2012] [Indexed: 11/07/2022]
Abstract
The equilibrium association free enthalpies ΔG(a) for typical supramolecular complexes in solution are calculated by ab initio quantum chemical methods. Ten neutral and three positively charged complexes with experimental ΔG(a) values in the range 0 to -21 kcal mol(-1) (on average -6 kcal mol(-1)) are investigated. The theoretical approach employs a (nondynamic) single-structure model, but computes the various energy terms accurately without any special empirical adjustments. Dispersion corrected density functional theory (DFT-D3) with extended basis sets (triple-ζ and quadruple-ζ quality) is used to determine structures and gas-phase interaction energies (ΔE), the COSMO-RS continuum solvation model (based on DFT data) provides solvation free enthalpies and the remaining ro-vibrational enthalpic/entropic contributions are obtained from harmonic frequency calculations. Low-lying vibrational modes are treated by a free-rotor approximation. The accurate account of London dispersion interactions is mandatory with contributions in the range -5 to -60 kcal mol(-1) (up to 200% of ΔE). Inclusion of three-body dispersion effects improves the results considerably. A semilocal (TPSS) and a hybrid density functional (PW6B95) have been tested. Although the ΔG(a) values result as a sum of individually large terms with opposite sign (ΔE vs. solvation and entropy change), the approach provides unprecedented accuracy for ΔG(a) values with errors of only 2 kcal mol(-1) on average. Relative affinities for different guests inside the same host are always obtained correctly. The procedure is suggested as a predictive tool in supramolecular chemistry and can be applied routinely to semirigid systems with 300-400 atoms. The various contributions to binding and enthalpy-entropy compensations are discussed.
Collapse
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| |
Collapse
|
26
|
Gao W, Feng H, Xuan X, Chen L. The assessment and application of an approach to noncovalent interactions: the energy decomposition analysis (EDA) in combination with DFT of revised dispersion correction (DFT-D3) with Slater-type orbital (STO) basis set. J Mol Model 2012; 18:4577-89. [PMID: 22643965 DOI: 10.1007/s00894-012-1425-0] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2012] [Accepted: 04/03/2012] [Indexed: 11/28/2022]
Abstract
An assessment study is presented about energy decomposition analysis (EDA) in combination with DFT including revised dispersion correction (DFT-D3) with Slater-type orbital (STO) basis set. There has been little knowledge about the performance of the EDA + DFT-D3 concerning STOs. In this assessment such an approach was applied to calculate noncovalent interaction energies and their corresponding components. Complexes in S22 set were used to evaluate the performance of EDA in conjunction with four representative types of GGA-functionals of DFT-D3 (BP86-D3, BLYP-D3, PBE-D3 and SSB-D3) with three STO basis sets ranging in complexity from DZP, TZ2P to QZ4P. The results showed that the approach of EDA + BLYP-D3/TZ2P has a better performance not only in terms of calculating noncovalent interaction energy quantitatively but also in analyzing corresponding energy components qualitatively. This approach (EDA + BLYP-D3/TZ2P) was thus applied further to two representative large-system complexes including porphine dimers and fullerene aggregates to gain a better insight into binding characteristics.
Collapse
Affiliation(s)
- Wei Gao
- KLGHEI of Environment and Energy Chemistry, School of Chemistry and Chemical Engineering, Sun Yat-sen University, Guangzhou 510275, China
| | | | | | | |
Collapse
|
27
|
Bozkaya U. Orbital-optimized third-order Møller-Plesset perturbation theory and its spin-component and spin-opposite scaled variants: Application to symmetry breaking problems. J Chem Phys 2011; 135:224103. [DOI: 10.1063/1.3665134] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Cabaleiro-Lago EM, Rodríguez-Otero J, Peña-Gallego Á. Effect of microhydration on the guanidinium⋯benzene interaction. J Chem Phys 2011; 135:214301. [DOI: 10.1063/1.3663277] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
|
29
|
Cabaleiro-Lago EM, Rodríguez-Otero J, Peña-Gallego Á. Computational study of the interaction of indole-like molecules with water and hydrogen sulfide. J Chem Phys 2011; 135:134310. [DOI: 10.1063/1.3643840] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
Abstract
Popular explanations of substituent effects in π-stacking interactions hinge upon substituent-induced changes in the aryl π-system. This entrenched view has been used to explain substituent effects in countless stacking interactions over the past 2 decades. However, for a broad range of stacked dimers, it is shown that substituent effects are better described as arising from local, direct interactions of the substituent with the proximal vertex of the other ring. Consequently, substituent effects in stacking interactions are additive, regardless of whether the substituents are on the same or opposite rings. Substituent effects are also insensitive to the introduction of heteroatoms on distant parts of either stacked ring. This local, direct interaction viewpoint provides clear, unambiguous explanations of substituent effects for myriad stacking interactions that are in accord with robust computational data, including DFT-D and new benchmark CCSD(T) results. Many of these computational results cannot be readily explained using traditional π-polarization-based models. Analyses of stacking interactions based solely on the sign of the electrostatic potential above the face of an aromatic ring or the molecular quadrupole moment face a similar fate. The local, direct interaction model provides a simple means of analyzing substituent effects in complex aromatic systems and also offers simple explanations of the crystal packing of fluorinated benzenes and the recently published dependence of the stability of protein-RNA complexes on the regiochemistry of fluorinated base analogues [J. Am. Chem. Soc.2011, 133, 3687-3689].
Collapse
Affiliation(s)
- Steven E Wheeler
- Department of Chemistry, Texas A&M University, College Station, Texas 77843, USA.
| |
Collapse
|
31
|
Burns LA, Mayagoitia ÁV, Sumpter BG, Sherrill CD. Density-functional approaches to noncovalent interactions: A comparison of dispersion corrections (DFT-D), exchange-hole dipole moment (XDM) theory, and specialized functionals. J Chem Phys 2011; 134:084107. [DOI: 10.1063/1.3545971] [Citation(s) in RCA: 543] [Impact Index Per Article: 41.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
|
32
|
Cabaleiro-Lago EM, Rodríguez-Otero J, Peña-Gallego Á. Study of the interaction between aniline and CH3CN, CH3Cl and CH3F. Theor Chem Acc 2010. [DOI: 10.1007/s00214-010-0789-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
33
|
Grimme S, Antony J, Ehrlich S, Krieg H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J Chem Phys 2010; 132:154104. [PMID: 20423165 DOI: 10.1063/1.3382344] [Citation(s) in RCA: 23959] [Impact Index Per Article: 1711.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
The method of dispersion correction as an add-on to standard Kohn-Sham density functional theory (DFT-D) has been refined regarding higher accuracy, broader range of applicability, and less empiricism. The main new ingredients are atom-pairwise specific dispersion coefficients and cutoff radii that are both computed from first principles. The coefficients for new eighth-order dispersion terms are computed using established recursion relations. System (geometry) dependent information is used for the first time in a DFT-D type approach by employing the new concept of fractional coordination numbers (CN). They are used to interpolate between dispersion coefficients of atoms in different chemical environments. The method only requires adjustment of two global parameters for each density functional, is asymptotically exact for a gas of weakly interacting neutral atoms, and easily allows the computation of atomic forces. Three-body nonadditivity terms are considered. The method has been assessed on standard benchmark sets for inter- and intramolecular noncovalent interactions with a particular emphasis on a consistent description of light and heavy element systems. The mean absolute deviations for the S22 benchmark set of noncovalent interactions for 11 standard density functionals decrease by 15%-40% compared to the previous (already accurate) DFT-D version. Spectacular improvements are found for a tripeptide-folding model and all tested metallic systems. The rectification of the long-range behavior and the use of more accurate C(6) coefficients also lead to a much better description of large (infinite) systems as shown for graphene sheets and the adsorption of benzene on an Ag(111) surface. For graphene it is found that the inclusion of three-body terms substantially (by about 10%) weakens the interlayer binding. We propose the revised DFT-D method as a general tool for the computation of the dispersion energy in molecules and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems.
Collapse
Affiliation(s)
- Stefan Grimme
- Theoretische Organische Chemie, Organisch-Chemisches Institut, Universität Münster, Germany.
| | | | | | | |
Collapse
|
34
|
Wheeler SE, McNeil AJ, Müller P, Swager TM, Houk KN. Probing substituent effects in aryl-aryl interactions using stereoselective Diels-Alder cycloadditions. J Am Chem Soc 2010; 132:3304-11. [PMID: 20158182 DOI: 10.1021/ja903653j] [Citation(s) in RCA: 147] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Stereoselective Diels-Alder cycloadditions that probe substituent effects in aryl-aryl sandwich complexes were studied experimentally and theoretically. Computations on model systems demonstrate that the stereoselectivity in these reactions is mediated by differential pi-stacking interactions in competing transition states. This allows relative stacking free energies of substituted and unsubstituted sandwich complexes to be derived from measured product distributions. In contrast to gas-phase computations, dispersion effects do not appear to play a significant role in the substituent effects, in accord with previous experiments. The experimental pi-stacking free energies are shown to correlate well with Hammett sigma(m) constants (r = 0.96). These substituent constants primarily provide a measure of the inductive electron-donating and -withdrawing character of the substituents, not donation into or out of the benzene pi-system. The present experimental results are most readily explained using a recently proposed model of substituent effects in the benzene sandwich dimer in which the pi-system of the substituted benzene is relatively unimportant and substituent effects arise from direct through-space interactions. Specifically, these results are the first experiments to clearly show that OMe enhances these pi-stacking interactions, despite being a pi-electron donor. This is in conflict with popular models in which substituent effects in aryl-aryl interactions are modulated by polarization of the aryl pi-system.
Collapse
Affiliation(s)
- Steven E Wheeler
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, USA
| | | | | | | | | |
Collapse
|
35
|
Cabaleiro-Lago EM, Carrazana-García JA, Rodríguez-Otero J. Study of the interaction between water and hydrogen sulfide with polycyclic aromatic hydrocarbons. J Chem Phys 2009; 130:234307. [PMID: 19548727 DOI: 10.1063/1.3152577] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A computational study has been carried out for determining the characteristics of the interaction between one water and hydrogen sulfide molecule with a series of polycyclic aromatic hydrocarbons of increasing size, namely, benzene, anthracene, triphenylene, coronene, circumcoronene, and dicircumcoronene. Potential energy curves were calculated for structures where H(2)X (X=O,S) molecule is located over the central six-membered ring with its hydrogen atoms pointing toward to (mode A) or away from (mode B) the hydrocarbon. The accuracy of different methods has been tested against the results of coupled cluster calculations extrapolated to basis set limit for the smaller hydrocarbons. The spin component scaled MP2 (SCS-MP2) method and a density functional theory method empirically corrected for dispersion (DFT-D) reproduce fairly well the results of high level calculations and therefore were employed for studying the larger systems, though DFT-D seems to underestimate the interaction in hydrogen sulfide clusters. Water complexes in mode A have interaction energies that hardly change with the size of the hydrocarbon due to compensation between the increase in the correlation contribution to the interaction energy and the increase in the repulsive character of the Hartree-Fock energy. For all the other clusters studied, there is a continuous increase in the intensity of the interaction as the size of the hydrocarbon increases, suggesting already converged values for circumcoronene. The interaction energy for water clusters extrapolated to an infinite number of carbon atoms amounts to -13.0 and -15.8 kJ/mol with SCS-MP2 and DFT-D, respectively. Hydrogen sulfide interacts more strongly than water with the hydrocarbons studied, leading to a limiting value of -21.7 kJ/mol with the SCS-MP2 method. Also, complexes in mode B are less stable than the corresponding A structures, with interaction energies amounting to -8.2 and -18.2 kJ/mol for water and hydrogen sulfide, respectively. The DFT-D calculations give values of -16.2 and -9.3 kJ/mol for hydrogen sulfide complexes in modes A and B, less negative than those predicted by the SCS-MP2 method, probably indicating problems with sulfur dispersion parameters.
Collapse
Affiliation(s)
- Enrique M Cabaleiro-Lago
- Departamento de Química Física, Facultade de Ciencias, Universidade de Santiago de Compostela, Campus de Lugo, Avda. Alfonso X El Sabio s/n, 27002 Lugo, Galicia, Spain.
| | | | | |
Collapse
|
36
|
Janesko BG, Henderson TM, Scuseria GE. Long-range-corrected hybrid density functionals including random phase approximation correlation: Application to noncovalent interactions. J Chem Phys 2009; 131:034110. [DOI: 10.1063/1.3176514] [Citation(s) in RCA: 81] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
|