1
|
Solovyova IV, Yang S, Starovoytov ON. Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water. J Chem Phys 2023; 158:084505. [PMID: 36859108 DOI: 10.1063/5.0134465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023] Open
Abstract
The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL-water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]+ cation and nitrate [NO3]- anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔVE, and liquid structures gr. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid-water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid-water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ∼0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. The ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.
Collapse
Affiliation(s)
- Iana V Solovyova
- Institute of Engineering Science, Ural Branch of the Russian Academy of Sciences, Yekaterinburg 620049, Russia
| | - Shizhong Yang
- Southern University and A&M College, Department of Computer Science, Baton Rouge, Louisiana 70807, USA
| | - Oleg N Starovoytov
- Southern University and A&M College, Department of Computer Science, Baton Rouge, Louisiana 70807, USA
| |
Collapse
|
2
|
Bernal DE, Ajagekar A, Harwood SM, Stober ST, Trenev D, You F. Perspectives of Quantum Computing for Chemical Engineering. AIChE J 2022. [DOI: 10.1002/aic.17651] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- David E. Bernal
- Research Institute for Advanced Computer Science Universities Space Research Association Mountain View California USA
- Quantum Artificial Intelligence Laboratory (QuAIL) NASA Ames Research Center Moffett Field California USA
- Department of Chemical Engineering Carnegie Mellon University Pittsburgh Pennsylvania USA
| | | | - Stuart M. Harwood
- Corporate Strategic Research ExxonMobil Research and Engineering Clinton New Jersey USA
| | - Spencer T. Stober
- Corporate Strategic Research ExxonMobil Research and Engineering Clinton New Jersey USA
| | - Dimitar Trenev
- Corporate Strategic Research ExxonMobil Research and Engineering Clinton New Jersey USA
| | - Fengqi You
- Systems Engineering Cornell University New York USA
- Robert Frederick Smith School of Chemical and Biomolecular Engineering Cornell University New York USA
| |
Collapse
|
3
|
Nieuwland C, Hamlin TA, Fonseca Guerra C, Barone G, Bickelhaupt FM. B-DNA Structure and Stability: The Role of Nucleotide Composition and Order. ChemistryOpen 2022; 11:e202100231. [PMID: 35083880 PMCID: PMC8805170 DOI: 10.1002/open.202100231] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2021] [Revised: 12/10/2021] [Indexed: 11/08/2022] Open
Abstract
We have quantum chemically analyzed the influence of nucleotide composition and sequence (that is, order) on the stability of double-stranded B-DNA triplets in aqueous solution. To this end, we have investigated the structure and bonding of all 32 possible DNA duplexes with Watson-Crick base pairing, using dispersion-corrected DFT at the BLYP-D3(BJ)/TZ2P level and COSMO for simulating aqueous solvation. We find enhanced stabilities for duplexes possessing a higher GC base pair content. Our activation strain analyses unexpectedly identify the loss of stacking interactions within individual strands as a destabilizing factor in the duplex formation, in addition to the better-known effects of partial desolvation. Furthermore, we show that the sequence-dependent differences in the interaction energy for duplexes of the same overall base pair composition result from the so-called "diagonal interactions" or "cross terms". Whether cross terms are stabilizing or destabilizing depends on the nature of the electrostatic interaction between polar functional groups in the pertinent nucleobases.
Collapse
Affiliation(s)
- Celine Nieuwland
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdam (TheNetherlands
| | - Trevor A. Hamlin
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdam (TheNetherlands
| | - Célia Fonseca Guerra
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdam (TheNetherlands
- Leiden Institute of ChemistryGorlaeus LaboratoriesLeiden UniversityEinsteinweg 552300 CCLeiden (TheNetherlands
| | - Giampaolo Barone
- Dipartimento di Scienze e Tecnologie Biologiche, Chimiche e FarmaceuticheUniversità degli Studi di PalermoViale delle Scienze, Edificio 1790128PalermoItaly
| | - F. Matthias Bickelhaupt
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdam (TheNetherlands
- Institute of Molecules and MaterialsRadboud University NijmegenHeyendaalseweg 1356525 AJNijmegen (TheNetherlands
| |
Collapse
|
4
|
Rostkowski M, Schürner HKV, Sowińska A, Vasquez L, Przydacz M, Elsner M, Dybala-Defratyka A. Isotope Effects on the Vaporization of Organic Compounds from an Aqueous Solution-Insight from Experiment and Computations. J Phys Chem B 2021; 125:13868-13885. [PMID: 34908428 PMCID: PMC8724799 DOI: 10.1021/acs.jpcb.1c05574] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
An isotope fractionation
analysis of organic groundwater pollutants
can assess the remediation at contaminated sites yet needs to consider
physical processes as potentially confounding factors. This study
explores the predictability of water–air partitioning isotope
effects from experiments and computational predictions for benzene
and trimethylamine (both H-bond acceptors) as well as chloroform (H-bond
donor). A small, but significant, isotope fractionation of different
direction and magnitude was measured with ε = −0.12‰
± 0.07‰ (benzene), εC = 0.49‰
± 0.23‰ (triethylamine), and εH = 1.79‰
± 0.54‰ (chloroform) demonstrating that effects do not
correlate with expected hydrogen-bond functionalities. Computations
revealed that the overall isotope effect arises from contributions
of different nature and extent: a weakening of intramolecular vibrations
in the condensed phase plus additional vibrational modes from a complexation
with surrounding water molecules. Subtle changes in benzene contrast
with a stronger coupling between intra- and intermolecular modes in
the chloroform–water system and a very local vibrational response
with few atoms involved in a specific mode of triethylamine. An energy
decomposition analysis revealed that each system was affected differently
by electrostatics and dispersion, where dispersion was dominant for
benzene and electrostatics dominated for chloroform and triethylamine.
Interestingly, overall stabilization patterns in all studied systems
originated from contributions of dispersion rather than other energy
terms.
Collapse
Affiliation(s)
- Michał Rostkowski
- Institute of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology, Zeromskiego 116, 90-924 Lodz, Poland
| | - Heide K V Schürner
- Chair of Analytical Chemistry and Water Chemistry, Technical University of Munich, Elisabeth-Winterhalter-Weg 6, 81377 Munich, Germany
| | - Agata Sowińska
- Institute of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology, Zeromskiego 116, 90-924 Lodz, Poland
| | - Luis Vasquez
- Institute of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology, Zeromskiego 116, 90-924 Lodz, Poland
| | - Martyna Przydacz
- Institute of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology, Zeromskiego 116, 90-924 Lodz, Poland
| | - Martin Elsner
- Chair of Analytical Chemistry and Water Chemistry, Technical University of Munich, Elisabeth-Winterhalter-Weg 6, 81377 Munich, Germany
| | - Agnieszka Dybala-Defratyka
- Institute of Applied Radiation Chemistry, Faculty of Chemistry, Lodz University of Technology, Zeromskiego 116, 90-924 Lodz, Poland
| |
Collapse
|
5
|
Alghanmi RM, Basha MT, Soliman SM, Alsaeedi RK. New Charge Transfer Complexes of K +-Channel-Blocker Drug (Amifampridine; AMFP) for Sensitive Detection; Solution Investigations and DFT Studies. Molecules 2021; 26:molecules26196037. [PMID: 34641581 PMCID: PMC8512129 DOI: 10.3390/molecules26196037] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Revised: 09/26/2021] [Accepted: 09/28/2021] [Indexed: 11/16/2022] Open
Abstract
UV-Vis spectroscopy was used to investigate two new charge transfer (CT) complexes formed between the K+-channel-blocker amifampridine (AMFP) drug and the two π-acceptors 2,3-dichloro-5,6-dicyano-p-benzoquinone (DDQ) and tetracyanoethylene (TCNE) in different solvents. The molecular composition of the new CT complexes was estimated using the continuous variations method and found to be 1:1 for both complexes. The formed CT complexes' electronic spectra data were further employed for calculating the formation constants (KCT), molar extinction coefficients (εCT), and physical parameters at various temperatures, and the results demonstrated the high stability of both complexes. In addition, sensitive spectrophotometric methods for quantifying AMFP in its pure form were proposed and statistically validated. Furthermore, DFT calculations were used to predict the molecular structures of AMFP-DDQ and AMFP-TCNE complexes in CHCl3. TD-DFT calculations were also used to predict the electronic spectra of both complexes. A CT-based transition band (exp. 399 and 417 nm) for the AMFP-TCNE complex was calculated at 411.5 nm (f = 0.105, HOMO-1 → LUMO). The two absorption bands at 459 nm (calc. 426.9 nm, f = 0.054) and 584 nm (calc. 628.1 nm, f = 0.111) of the AMFP-DDQ complex were theoretically assigned to HOMO-1 → LUMO and HOMO → LUMO excitations, respectively.
Collapse
Affiliation(s)
- Reem M. Alghanmi
- Department of Chemistry, College of Science, University of Jeddah, P.O. 80327, Jeddah 21589, Saudi Arabia; (M.T.B.); (R.K.A.)
- Correspondence:
| | - Maram T. Basha
- Department of Chemistry, College of Science, University of Jeddah, P.O. 80327, Jeddah 21589, Saudi Arabia; (M.T.B.); (R.K.A.)
| | - Saied M. Soliman
- Department of Chemistry, Faculty of Science, Alexandria University, P.O. 426, Ibrahimia, Alexandria 21525, Egypt;
| | - Razan K. Alsaeedi
- Department of Chemistry, College of Science, University of Jeddah, P.O. 80327, Jeddah 21589, Saudi Arabia; (M.T.B.); (R.K.A.)
| |
Collapse
|
6
|
Ahmad S, Bühl M. Computational modelling of Pd-catalysed alkoxycarbonylation of alkenes and alkynes. Phys Chem Chem Phys 2021; 23:15869-15880. [PMID: 34318843 DOI: 10.1039/d1cp02426d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
This perspective highlights the computational modelling of alkene and alkyne alkoxycarbonylation at palladium catalysts. We cover studies on Pd-catalysed alkoxycarbonylation of alkenes with bidentate diphosphine ligands, which reveal a hydride pathway is operating with an intermolecular alcoholysis step, where explicit solvation is mandatory to estimate the overall barriers correctly and model alcoholysis/copolymerisation selectivities. Subsequently, we discuss Pd-catalysed alkyne alkoxycarbonylation with P,N-chelating ligands, where an in situ base mechanism is operating involving ketene-type intermediates. We also discuss catalyst poisoning due to allene and designing a potential new catalyst tolerant towards allene poisoning.
Collapse
Affiliation(s)
- Shahbaz Ahmad
- School of Chemistry, University of St Andrews, North Haugh, St Andrews, Fife KY16 9ST, UK.
| | | |
Collapse
|
7
|
Berro Y, Badawi M, El Haj Hassan F, Kassir M, Tielens F. Water-silanol interactions on the amorphous silica surface: A dispersion-corrected DFT investigation. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.114496] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
8
|
Biswas S, Mallik BS. Probing the vibrational dynamics of amide bands of N-methylformamide, N, N-dimethylacetamide, and N-methylacetamide in water. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.113001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
9
|
Cabral Tenorio BN, Chaer Nascimento MA, Rocha AB. X-ray Photoionization Cross Section Spectra of Water and Ammonia Bonded on Polycyclic Aromatic Hydrocarbons: A Quantum Mechanical Interpretation to the Absorption Spectra on Graphene. J Phys Chem A 2020; 124:2591-2600. [PMID: 32187493 DOI: 10.1021/acs.jpca.9b11406] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A theoretical study of the K-shell total photoabsorption and photoionization cross section spectra of water and ammonia bonded to benzene (C6H6) and the polycyclic aromatic hydrocarbons (PAHs) naphthalene (C10H8), coronene (C24H12) and circumcoronene (C54H18) by van der Waals (vdW) forces is presented. The discretized electronic pseudospectra at the oxygen and nitrogen K-edges, covering the discrete and the continuum spectral regions, were obtained at the time-dependent density functional theory (TDDFT) level with dispersion correction. An analytic continuation procedure based on the Padé approximants was used in order to obtain the K-shell cross sections of the structures at the discrete and the continuum regions of the spectra. By examining the electronic spectra of water and ammonia bonded to coronene and circumcoronene, we observed that our results agree well with the experiments performed with graphene. This work provides a quantum mechanical interpretation to the NEXAFS experiments of water and ammonia adsorbed on graphene in terms of a physisorption model of these molecules by van der Waals forces.
Collapse
Affiliation(s)
- Bruno Nunes Cabral Tenorio
- Instituto de Química, Universidade Federal do Rio de Janeiro, UFRJ, Avenue Athos da Silveira Ramos, 149, Rio de Janeiro - RJ 21941-909, Brazil
| | - Marco Antonio Chaer Nascimento
- Instituto de Química, Universidade Federal do Rio de Janeiro, UFRJ, Avenue Athos da Silveira Ramos, 149, Rio de Janeiro - RJ 21941-909, Brazil
| | - Alexandre Braga Rocha
- Instituto de Química, Universidade Federal do Rio de Janeiro, UFRJ, Avenue Athos da Silveira Ramos, 149, Rio de Janeiro - RJ 21941-909, Brazil
| |
Collapse
|
10
|
Spiegelman F, Tarrat N, Cuny J, Dontot L, Posenitskiy E, Martí C, Simon A, Rapacioli M. Density-functional tight-binding: basic concepts and applications to molecules and clusters. ADVANCES IN PHYSICS: X 2020; 5:1710252. [PMID: 33154977 PMCID: PMC7116320 DOI: 10.1080/23746149.2019.1710252] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2019] [Accepted: 12/19/2019] [Indexed: 06/10/2023] Open
Abstract
The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.
Collapse
Affiliation(s)
- Fernand Spiegelman
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Nathalie Tarrat
- CEMES, Université de Toulouse (UPS), CNRS, UPR8011, Toulouse, Toulouse, France
| | - Jérôme Cuny
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Leo Dontot
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Evgeny Posenitskiy
- Laboratoire Collisions Agrégats et Réactivité LCAR/IRSAMC, UMR5589, Université de Toulouse (UPS) and CNRS, Toulouse, France
| | - Carles Martí
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
- Laboratoire de Chimie, UMR5182, Ecole Normale Supérieure de Lyon, Université de Lyon and CNRS, Lyon, France
| | - Aude Simon
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Mathias Rapacioli
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| |
Collapse
|
11
|
Shoji A, Schanzenbach D, Merrill R, Zhang J, Yang L, Sun R. Theoretical Study of the Potential Energy Profile of the HBr + + CO 2 → HOCO + + Br· Reaction. J Phys Chem A 2019; 123:9791-9799. [PMID: 31633930 DOI: 10.1021/acs.jpca.9b07651] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Recent guided ion beam experiments have revealed interesting reaction dynamics of the HBr+ + CO2 → HOCO+ + Br· reaction under different conditions. The hypothesis is that the predominant reaction mechanism depends on the collision energy between two reactants, the angular momentum of HBr+, and the spin-orbit coupling state of the system. The potential energy profile of the HBr+ + CO2 → HOCO+ + Br· reaction is studied in this research to lay the groundwork for an ab initio molecular dynamics simulation. First, a benchmark potential energy profile of this reaction is identified using coupled-cluster theory extrapolated to the complete basis set limit. A transition state connecting the previously reported intermediates is found, making the potential energy surface of the HBr+ + CO2 → HOCO+ + Br· reaction double-welled. Second, various single reference ab initio methods are compared with the benchmark potential energy profile to search for the most suitable ab initio method for the dynamics simulation. Two combinations of double-ζ basis sets (with effective core potentials) with MP2 and density functional theory have been identified to accurately represent the potential energy profile of this reaction.
Collapse
Affiliation(s)
- Alyson Shoji
- Department of Chemistry , The University of Hawai'i at Manoa , Honolulu , Hawaii 96822 , United States
| | - David Schanzenbach
- Information Technology Services , The University of Hawai'i , Honolulu , Hawaii 96822 , United States
| | - Ron Merrill
- Information Technology Services , The University of Hawai'i , Honolulu , Hawaii 96822 , United States
| | - Jiaxu Zhang
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering , Harbin Institute of Technology , Harbin 150001 , P. R. China
| | - Li Yang
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering , Harbin Institute of Technology , Harbin 150001 , P. R. China
| | - Rui Sun
- Department of Chemistry , The University of Hawai'i at Manoa , Honolulu , Hawaii 96822 , United States
| |
Collapse
|
12
|
Kim Y, Song S, Sim E, Burke K. Halogen and Chalcogen Binding Dominated by Density-Driven Errors. J Phys Chem Lett 2019; 10:295-301. [PMID: 30562033 DOI: 10.1021/acs.jpclett.8b03745] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
Dispersion corrections of various kinds usually improve DFT energetics of weak noncovalent interactions. However, in some cases involving molecules or halides, especially those with σ-hole interactions, the density-driven errors of uncorrected DFT are larger than the dispersion corrections. In these abnormal situations, HF-DFT (using Hartree-Fock densities instead of self-consistent densities) greatly improves bond energies, while dispersion corrections can even worsen the results. On the other hand, pnictogen bonds and the S22 data set are normal and are not improved by this procedure. Such effects should be accounted for when parametrizing dispersion interactions.
Collapse
Affiliation(s)
- Yeil Kim
- Department of Chemistry , Yonsei University , 50 Yonsei-ro , Seodaemun-gu, Seoul 03722 , Korea
| | - Suhwan Song
- Department of Chemistry , Yonsei University , 50 Yonsei-ro , Seodaemun-gu, Seoul 03722 , Korea
| | - Eunji Sim
- Department of Chemistry , Yonsei University , 50 Yonsei-ro , Seodaemun-gu, Seoul 03722 , Korea
| | - Kieron Burke
- Department of Chemistry , and Department of Physics , University of California , Irvine , California 92697 , United States
| |
Collapse
|
13
|
Caldeweyher E, Bannwarth C, Grimme S. Extension of the D3 dispersion coefficient model. J Chem Phys 2018; 147:034112. [PMID: 28734285 DOI: 10.1063/1.4993215] [Citation(s) in RCA: 460] [Impact Index Per Article: 76.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
A new model, termed D4, for the efficient computation of molecular dipole-dipole dispersion coefficients is presented. As in the related, well established D3 scheme, these are obtained as a sum of atom-in-molecule dispersion coefficients over atom pairs. Both models make use of dynamic polarizabilities obtained from first-principles time-dependent density functional theory calculations for atoms in different chemical environments employing fractional atomic coordination numbers for interpolation. Different from the D3 model, the coefficients are obtained on-the-fly by numerical Casimir-Polder integration of the dynamic, atomic polarizabilities α(iω). Most importantly, electronic density information is now incorporated via atomic partial charges computed at a semi-empirical quantum mechanical tight-binding level, which is used to scale the polarizabilities. Extended statistical measures show that errors for dispersion coefficients with the proposed D4 method are significantly lower than with D3 and other, computationally more involved schemes. Alongside, accurate isotropic charge and hybridization dependent, atom-in-molecule static polarizabilities are obtained with an unprecedented efficiency. Damping function parameters are provided for three standard density functionals, i.e., TPSS, PBE0, and B3LYP, allowing evaluation of the new DFT-D4 model for common non-covalent interaction energy benchmark sets.
Collapse
Affiliation(s)
- Eike Caldeweyher
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Christoph Bannwarth
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
14
|
Vitillo JG, Fjermestad T, D’Amore M, Milanesio M, Palin L, Ricchiardi G, Bordiga S. On the structure of superbasic (MgO)n sites solvated in a faujasite zeolite. Phys Chem Chem Phys 2018; 20:18503-18514. [DOI: 10.1039/c8cp01788c] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Theory and experiment reveal the structure of magnesium oxide nanoclusters in a superbasic faujasite zeolite.
Collapse
Affiliation(s)
- Jenny G. Vitillo
- Department of Chemistry
- Chemical Theory Center, and Minnesota Supercomputing Institute
- University of Minnesota
- Minneapolis
- USA
| | - Torstein Fjermestad
- Department of Chemistry and NIS Interdepartment Centre
- University of Turin
- 10125 Torino
- Italy
| | - Maddalena D’Amore
- Department of Chemistry and NIS Interdepartment Centre
- University of Turin
- 10125 Torino
- Italy
| | | | - Luca Palin
- Nova Res s.r.l
- Novara
- Italy
- Dipartimento di Scienze e Innovazione Tecnologica
- Università del Piemonte Orientale
| | - Gabriele Ricchiardi
- Department of Chemistry and NIS Interdepartment Centre
- University of Turin
- 10125 Torino
- Italy
| | - Silvia Bordiga
- Department of Chemistry and NIS Interdepartment Centre
- University of Turin
- 10125 Torino
- Italy
| |
Collapse
|
15
|
Hermann J, DiStasio RA, Tkatchenko A. First-Principles Models for van der Waals Interactions in Molecules and Materials: Concepts, Theory, and Applications. Chem Rev 2017; 117:4714-4758. [PMID: 28272886 DOI: 10.1021/acs.chemrev.6b00446] [Citation(s) in RCA: 255] [Impact Index Per Article: 36.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Noncovalent van der Waals (vdW) or dispersion forces are ubiquitous in nature and influence the structure, stability, dynamics, and function of molecules and materials throughout chemistry, biology, physics, and materials science. These forces are quantum mechanical in origin and arise from electrostatic interactions between fluctuations in the electronic charge density. Here, we explore the conceptual and mathematical ingredients required for an exact treatment of vdW interactions, and present a systematic and unified framework for classifying the current first-principles vdW methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem (namely the Rutgers-Chalmers vdW-DF, Vydrov-Van Voorhis (VV), exchange-hole dipole moment (XDM), Tkatchenko-Scheffler (TS), many-body dispersion (MBD), and random-phase approximation (RPA) approaches). Particular attention is paid to the intriguing nature of many-body vdW interactions, whose fundamental relevance has recently been highlighted in several landmark experiments. The performance of these models in predicting binding energetics as well as structural, electronic, and thermodynamic properties is connected with the theoretical concepts and provides a numerical summary of the state-of-the-art in the field. We conclude with a roadmap of the conceptual, methodological, practical, and numerical challenges that remain in obtaining a universally applicable and truly predictive vdW method for realistic molecular systems and materials.
Collapse
Affiliation(s)
- Jan Hermann
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University , Ithaca, New York 14853, United States
| | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany.,Physics and Materials Science Research Unit, University of Luxembourg , L-1511 Luxembourg, Luxembourg
| |
Collapse
|
16
|
Arabi AA. Evaluating dispersion forces for optimization of van der Waals complexes using a non-empirical functional. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2016; 374:rsta.2016.0145. [PMID: 27698041 PMCID: PMC5052729 DOI: 10.1098/rsta.2016.0145] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 07/05/2016] [Indexed: 05/27/2023]
Abstract
Modelling dispersion interactions with traditional density functional theory (DFT) is a challenge that has been extensively addressed in the past decade. The exchange-dipole moment (XDM), among others, is a non-empirical add-on dispersion correction model in DFT. The functional PW86+PBE+XDM for exchange, correlation and dispersion, respectively, compromises an accurate functional for thermochemistry and for van der Waals (vdW) complexes at equilibrium and non-equilibrium geometries. To use this functional in optimizing vdW complexes, rather than computing single point energies, it is necessary to evaluate accurate forces. The purpose of this paper is to validate that, along the potential energy surface, the distance at which the energy is minimum is commensurate with the distance at which the forces vanish to zero. This test was validated for 10 rare gas diatomic molecules using various integration grids and different convergence criteria. It was found that the use of either convergence criterion, 10-6 or 10-8, in Gaussian09, does not affect the accuracy of computed optimal distances and binding energies. An ultra-fine grid needs to be used when computing accurate energies using generalized gradient approximation functionals.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.
Collapse
Affiliation(s)
- Alya A Arabi
- College of Natural and Health Sciences, Zayed University, Abu Dhabi, PO Box 144534, United Arab Emirates
| |
Collapse
|
17
|
Ucak UV, Ji H, Singh Y, Jung Y. A soft damping function for dispersion corrections with less overfitting. J Chem Phys 2016; 145:174104. [DOI: 10.1063/1.4965818] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Umit V. Ucak
- Graduate School of EEWS, KAIST, Daejeon, South Korea
| | - Hyunjun Ji
- Graduate School of EEWS, KAIST, Daejeon, South Korea
| | - Yashpal Singh
- Graduate School of EEWS, KAIST, Daejeon, South Korea
| | - Yousung Jung
- Graduate School of EEWS, KAIST, Daejeon, South Korea
| |
Collapse
|
18
|
Grimme S, Hansen A, Brandenburg JG, Bannwarth C. Dispersion-Corrected Mean-Field Electronic Structure Methods. Chem Rev 2016; 116:5105-54. [DOI: 10.1021/acs.chemrev.5b00533] [Citation(s) in RCA: 799] [Impact Index Per Article: 99.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical
Chemistry, Universität Bonn, 53113 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical
Chemistry, Universität Bonn, 53113 Bonn, Germany
| | | | - Christoph Bannwarth
- Mulliken Center for Theoretical
Chemistry, Universität Bonn, 53113 Bonn, Germany
| |
Collapse
|
19
|
Daschakraborty S, Kiefer PM, Miller Y, Motro Y, Pines D, Pines E, Hynes JT. Reaction Mechanism for Direct Proton Transfer from Carbonic Acid to a Strong Base in Aqueous Solution I: Acid and Base Coordinate and Charge Dynamics. J Phys Chem B 2016; 120:2271-80. [PMID: 26879554 DOI: 10.1021/acs.jpcb.5b12742] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Protonation by carbonic acid H2CO3 of the strong base methylamine CH3NH2 in a neutral contact pair in aqueous solution is followed via Car-Parrinello molecular dynamics simulations. Proton transfer (PT) occurs to form an aqueous solvent-stabilized contact ion pair within 100 fs, a fast time scale associated with the compression of the acid-base hydrogen-bond (H-bond), a key reaction coordinate. This rapid barrierless PT is consistent with the carbonic acid-protonated base pKa difference that considerably favors the PT, and supports the view of intact carbonic acid as potentially important proton donor in assorted biological and environmental contexts. The charge redistribution within the H-bonded complex during PT supports a Mulliken picture of charge transfer from the nitrogen base to carbonic acid without altering the transferring hydrogen's charge from approximately midway between that of a hydrogen atom and that of a proton.
Collapse
Affiliation(s)
- Snehasis Daschakraborty
- Department of Chemistry and Biochemistry, University of Colorado , Boulder, Colorado 80309-0215, United States
| | - Philip M Kiefer
- Department of Chemistry and Biochemistry, University of Colorado , Boulder, Colorado 80309-0215, United States
| | - Yifat Miller
- Department of Chemistry, Ben-Gurion University of the Negev , P.O. Box 653, Beer-Sheva 84105, Israel
| | - Yair Motro
- Department of Chemistry, Ben-Gurion University of the Negev , P.O. Box 653, Beer-Sheva 84105, Israel
| | - Dina Pines
- Department of Chemistry, Ben-Gurion University of the Negev , P.O. Box 653, Beer-Sheva 84105, Israel
| | - Ehud Pines
- Department of Chemistry, Ben-Gurion University of the Negev , P.O. Box 653, Beer-Sheva 84105, Israel
| | - James T Hynes
- Department of Chemistry and Biochemistry, University of Colorado , Boulder, Colorado 80309-0215, United States.,Ecole Normale Supérieure-PSL Research University, Chemistry Department, Sorbonne Universités-UPMC University Paris 06, CNRS UMR 8640 Pasteur, 24 rue Lhomond, 75005 Paris, France
| |
Collapse
|
20
|
Interactions between metal cations with H2 in the M+- H2 complexes: Performance of DFT and DFT-D methods. J CHEM SCI 2016. [DOI: 10.1007/s12039-016-1054-8] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
21
|
Willow SY, Zeng XC, Xantheas SS, Kim KS, Hirata S. Why Is MP2-Water "Cooler" and "Denser" than DFT-Water? J Phys Chem Lett 2016; 7:680-684. [PMID: 26821830 DOI: 10.1021/acs.jpclett.5b02430] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Density functional theory (DFT) with a dispersionless generalized gradient approximation (GGA) needs much higher temperature and pressure than the ambient conditions to maintain water in the liquid phase at the correct (1 g/cm(3)) density during first-principles simulations. Conversely, ab initio second-order many-body perturbation (MP2) calculations of liquid water require lower temperature and pressure than DFT/GGA to keep water liquid. Here we present a unifying explanation of these trends derived from classical water simulations using a polarizable force field with different sets of parameters. We show that the different temperatures and pressures between DFT/GGA and MP2 at which the simulated water displays the experimentally observed liquid structure under the ambient conditions can be largely explained by their differences in polarizability and dispersion interaction, respectively. In DFT/GGA, the polarizability and thus the induced dipole moments and the hydrogen-bond strength are all overestimated. This hinders the rotational motion of molecules and requires a higher temperature for DFT-water to be liquid. MP2 gives a stronger dispersion interaction and thus shorter intermolecular distances than dispersionless DFT/GGA, which is why MP2-water is denser than DFT-water under the same external pressure.
Collapse
Affiliation(s)
- Soohaeng Yoo Willow
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Xiao Cheng Zeng
- Department of Chemistry, University of Nebraska at Lincoln , 536 Hamilton Hall, Lincoln, Nebraska 68588, United States
| | - Sotiris S Xantheas
- Physical Sciences Division, Pacific Northwest National Laboratory , 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States
| | - Kwang S Kim
- Center for Superfunctional Materials, Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST) , Ulsan 689-798, Korea
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| |
Collapse
|
22
|
Wojtaś M, Bil A, Ga̧gor A, Medycki W, Kholkin AL. Phase stability and dynamics of hybrid organic–inorganic crystals [(CH3)3PH][SbCl4] and [(CH3)3PH][SbBr4]: a computational and NMR approach. CrystEngComm 2016. [DOI: 10.1039/c6ce00160b] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
23
|
Du QS, Liu PJ, Deng J. Empirical Correction to Molecular Interaction Energies in Density Functional Theory (DFT) for Methane Hydrate Simulation. J Chem Theory Comput 2015; 3:1665-72. [PMID: 26627612 DOI: 10.1021/ct700026d] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
A general and empirical method is proposed for correction of London dispersion and other deficiencies in density functional theory (DFT). This method is based on the existing Lennard-Jones (L-J) equation and van der Waals parameters. The benchmark of energy correction is set as the energy difference between DFT and more accurate methods, for example CCSD(T). The energy correction includes all differences between CCSD(T) and DFT, dispersion energy, configuration interaction, induction interaction, residual correlation, and other effects. The energy correction is expressed as a linear combination of van der Waals potentials of nonbonded atomic pairs. The combination coefficients are determined using a least-squares approach in a training set of molecular pairs. The coefficients then can be used for the energy corrections in DFT calculations in a molecular family. Three correction equations of molecular pair interaction energy, water-water, water-methane, and methane-methane, are derived for methane hydrate simulation. The correction equation of the water-water pair is applied in the DFT calculation of water pentamer, yielding good intermolecular potential energy surfaces (PES), very close to the results of CCSD(T) over the active interaction range from 2.1 Å to 8.0 Å.
Collapse
Affiliation(s)
- Qi-Shi Du
- Department of Chemistry, Hainan Normal University, Haikou, Hainan 571158, China, Key Laboratory of Subtropical Bioresource Conservation and Utilization, Guangxi University, Nanning, Guangxi 530004, China, and Nanning New and High-Tech Incubator, 68 Keyuan Avenue, Nanning, Guangxi 530007, China
| | - Peng-Jun Liu
- Department of Chemistry, Hainan Normal University, Haikou, Hainan 571158, China, Key Laboratory of Subtropical Bioresource Conservation and Utilization, Guangxi University, Nanning, Guangxi 530004, China, and Nanning New and High-Tech Incubator, 68 Keyuan Avenue, Nanning, Guangxi 530007, China
| | - Jun Deng
- Department of Chemistry, Hainan Normal University, Haikou, Hainan 571158, China, Key Laboratory of Subtropical Bioresource Conservation and Utilization, Guangxi University, Nanning, Guangxi 530004, China, and Nanning New and High-Tech Incubator, 68 Keyuan Avenue, Nanning, Guangxi 530007, China
| |
Collapse
|
24
|
Biswas S, Mallik BS. Aqueous solvation of an amide molecule from first principles molecular simulations: Structure, hydrogen bond dynamics and spectral signature. J Mol Liq 2015. [DOI: 10.1016/j.molliq.2015.11.005] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
|
25
|
Howe EN, Ball GE, Thordarson P. Step-by-step DFT analysis of the cooperativity in the binding of cations and anions to a tetratopic ion-pairing host. Supramol Chem 2015. [DOI: 10.1080/10610278.2015.1088946] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- Ethan N.W. Howe
- School of Chemistry, The University of New South Wales, Sydney, Australia
- Australian Centre for Nanomedicine and the ARC Centre of Excellence in Convergent Bio-Nano Science and Technology, The University of New South Wales, Sydney, Australia
| | - Graham E. Ball
- School of Chemistry, The University of New South Wales, Sydney, Australia
| | - Pall Thordarson
- School of Chemistry, The University of New South Wales, Sydney, Australia
- Australian Centre for Nanomedicine and the ARC Centre of Excellence in Convergent Bio-Nano Science and Technology, The University of New South Wales, Sydney, Australia
| |
Collapse
|
26
|
Liakos DG, Neese F. Is It Possible To Obtain Coupled Cluster Quality Energies at near Density Functional Theory Cost? Domain-Based Local Pair Natural Orbital Coupled Cluster vs Modern Density Functional Theory. J Chem Theory Comput 2015; 11:4054-63. [DOI: 10.1021/acs.jctc.5b00359] [Citation(s) in RCA: 199] [Impact Index Per Article: 22.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Dimitrios G. Liakos
- Max Planck Institute for Chemical Energy Conversion, Stiftstr. 32-34, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max Planck Institute for Chemical Energy Conversion, Stiftstr. 32-34, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
27
|
Schröder H, Creon A, Schwabe T. Reformulation of the D3(Becke–Johnson) Dispersion Correction without Resorting to Higher than C6 Dispersion Coefficients. J Chem Theory Comput 2015; 11:3163-70. [DOI: 10.1021/acs.jctc.5b00400] [Citation(s) in RCA: 57] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Heiner Schröder
- Center for Bioinformatics
and Institute of Physical Chemistry, University of Hamburg, Bundesstraße
43, 20146 Hamburg, Germany
| | - Anne Creon
- Center for Bioinformatics
and Institute of Physical Chemistry, University of Hamburg, Bundesstraße
43, 20146 Hamburg, Germany
| | - Tobias Schwabe
- Center for Bioinformatics
and Institute of Physical Chemistry, University of Hamburg, Bundesstraße
43, 20146 Hamburg, Germany
| |
Collapse
|
28
|
|
29
|
Kocakaya SÖ, Turgut Y, Pirinççioglu N. Enantiomeric discrimination of chiral organic salts by chiral aza-15-crown-5 ether with C 1 symmetry: experimental and theoretical approaches. J Mol Model 2015; 21:55. [PMID: 25701087 DOI: 10.1007/s00894-015-2604-6] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2014] [Accepted: 02/01/2015] [Indexed: 11/25/2022]
Abstract
The work involves an experimental ((1)H NMR) and theoretical (MD, MM-PBSA and DFT) investigation of the molecular recognition and discrimination properties of a chiral aza-15-crown-5 against methyl esters of alanine, phenylalanine and valine hydrochloride salts. The results indicate that the receptor binds enantiomers with moderate binding constants (88-1,389 M(-1)), with phenylalanine being more discriminated. The difference in experimental binding free energies (ΔG(R) - ΔG(S)) for alanine, phenylalanine and valine enantiomers were calculated as -0.36, -1.58 and 0.80 kcal mol(-1), respectively. The differences in theoretical binding energies were calculated by MM-PBSA (ΔE(R)PB - ΔE(S)PB=) as -0.30, -1.45 and 0.88, by B3LYP/6-31+G(d) as -1.17, -0.84 and 0.74 and by M06-2X/6-31+G(d) as -1.40, -3.26 and 1.66 kcal mol(-1). The data obtained give valuable information regarding the molecular recognition mode of the organoammonium complexes of chiral aza-crown ether with C 1 symmetry, which may be relevant to biological systems.
Collapse
Affiliation(s)
- Safak Özhan Kocakaya
- Department of Chemistry, Faculty of Science, University of Dicle, Diyarbakir, 21280, Turkey
| | | | | |
Collapse
|
30
|
Ottonello G, Richet P, Vetuschi Zuccolini M. The wet solidus of silica: Predictions from the scaled particle theory and polarized continuum model. J Chem Phys 2015; 142:054503. [DOI: 10.1063/1.4906745] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- G. Ottonello
- DIPTERIS, Università di Genova, Corso Europa 26, 16132 Genoa, Italy
| | - P. Richet
- Institut de Physique du Globe, Rue Jussieu 2, 75005 Paris, France
| | | |
Collapse
|
31
|
Amelines-Sarria O, Basiuk VA, Duarte-Alaniz V, Rivera M. Properties of noncovalent tetraphenylporphine⋯C60 dyads as studied by different long-range and dispersion-corrected DFT functionals. Phys Chem Chem Phys 2015; 17:27399-408. [DOI: 10.1039/c5cp03921e] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
The exceptional ability of M05-2X/6-31G(d,p) calculations to describe binding energies and separation distances in porphyrin–fullerene complexes.
Collapse
Affiliation(s)
- Oscar Amelines-Sarria
- Instituto de Física
- Universidad Nacional Autónoma de México
- Sendero Bicipuma
- Cd. Universitaria
- Del. Coyoacán
| | - Vladimir A. Basiuk
- Instituto de Ciencias Nucleares
- Universidad Nacional Autónoma de México
- Circuito Exterior s/n
- Cd. Universitaria
- Del. Coyoacán
| | - Víctor Duarte-Alaniz
- Instituto de Química
- Universidad Nacional Autónoma de México
- Circuito Exterior
- Cd. Universitaria
- Del. Coyoacán
| | - Margarita Rivera
- Instituto de Física
- Universidad Nacional Autónoma de México
- Sendero Bicipuma
- Cd. Universitaria
- Del. Coyoacán
| |
Collapse
|
32
|
Chakrabarti D, Kusumaatmaja H, Rühle V, Wales DJ. Exploring energy landscapes: from molecular to mesoscopic systems. Phys Chem Chem Phys 2014; 16:5014-25. [PMID: 24067895 DOI: 10.1039/c3cp52603h] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
We review a comprehensive computational framework to survey the potential energy landscape for systems composed of rigid or partially rigid molecules. Illustrative case studies relevant to a wide range of molecular clusters and soft and condensed matter systems are discussed.
Collapse
Affiliation(s)
- Dwaipayan Chakrabarti
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| | | | | | | |
Collapse
|
33
|
Bergendahl LT, Paterson MJ. Excited states of porphyrin and porphycene aggregates: Computational insights. COMPUT THEOR CHEM 2014. [DOI: 10.1016/j.comptc.2014.03.014] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
34
|
Briggs EA, Besley NA. Modelling excited states of weakly bound complexes with density functional theory. Phys Chem Chem Phys 2014; 16:14455-62. [DOI: 10.1039/c3cp55361b] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Different dispersion correction parameters are required to describe the interaction when the molecule is in an excited Rydberg state.
Collapse
Affiliation(s)
- Edward A. Briggs
- School of Chemistry
- University of Nottingham
- University Park
- Nottingham, UK
| | - Nicholas A. Besley
- School of Chemistry
- University of Nottingham
- University Park
- Nottingham, UK
| |
Collapse
|
35
|
Navrátil V, Klusák V, Rulíšek L. Theoretical Aspects of Hydrolysis of Peptide Bonds by Zinc Metalloenzymes. Chemistry 2013; 19:16634-45. [DOI: 10.1002/chem.201302663] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2013] [Revised: 08/08/2013] [Indexed: 12/11/2022]
Affiliation(s)
- Václav Navrátil
- Institute of Organic Chemistry and Biochemistry, Gilead Sciences & IOCB Research Center, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 166 10 Praha 6 (Czech Republic)
- Charles University in Prague, Faculty of Science, Department of Biochemistry, Hlavova 2030, 128 43 Praha 2 (Czech Republic)
| | - Vojtěch Klusák
- Institute of Organic Chemistry and Biochemistry, Gilead Sciences & IOCB Research Center, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 166 10 Praha 6 (Czech Republic)
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry, Gilead Sciences & IOCB Research Center, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 166 10 Praha 6 (Czech Republic)
| |
Collapse
|
36
|
Barone G, Fonseca Guerra C, Bickelhaupt FM. B-DNA Structure and Stability as Function of Nucleic Acid Composition: Dispersion-Corrected DFT Study of Dinucleoside Monophosphate Single and Double Strands. ChemistryOpen 2013; 2:186-93. [PMID: 24551565 PMCID: PMC3892189 DOI: 10.1002/open.201300019] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2013] [Indexed: 11/16/2022] Open
Abstract
We have computationally investigated the structure and stability of all 16 combinations of two out of the four natural DNA bases A, T, G and C in a di-2′-deoxyribonucleoside-monophosphate model DNA strand as well as in 10 double-strand model complexes thereof, using dispersion-corrected density functional theory (DFT-D). Optimized geometries with B-DNA conformation were obtained through the inclusion of implicit water solvent and, in the DNA models, of sodium counterions, to neutralize the negative charge of the phosphate groups. The results obtained allowed us to compare the relative stability of isomeric single and double strands. Moreover, the energy of the Watson–Crick pairing of complementary single strands to form double-helical structures was calculated. The latter furnished the following increasing stability trend of the double-helix formation energy: d(TpA)2 <d(CpA)2 <d(ApT)2 <d(ApA)2 <d(GpT)2 <d(GpA)2 <d(ApG)2 <d(CpG)2 <d(GpG)2 <d(GpC)2, where the energy differences between the last four dimers, d(ApG)2, d(CpG)2, d(GpG)2 and d(GpC)2, is within 4.0 kcal mol−1, and the energy between the most and the least stable isomers is 13.4 kcal mol−1. This trend shows that the formation energy essentially increases with the number of hydrogen bonds per base pair, that is two between A and T and three between G and C. Superimposed on this main trend are more subtle effects that depend on the order in which bases occur within a strand from the 5’- to the 3’-end.
Collapse
Affiliation(s)
- Giampaolo Barone
- Dipartimento di Scienze e Tecnologie Biologiche, Chimiche e Farmaceutiche, Università di Palermo Viale delle Scienze, Edificio 17, 90128, Palermo (Italy) E-mail:
| | - Célia Fonseca Guerra
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, VU University De Boelelaan 1083, 1081 HV Amsterdam (The Netherlands) E-mail:
| | - F Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, VU University De Boelelaan 1083, 1081 HV Amsterdam (The Netherlands) E-mail: ; Institute for Molecules and Materials, Radboud University Nijmegen Heyendaalseweg 135, 6525 AJ Nijmegen (The Netherlands)
| |
Collapse
|
37
|
Pastore M, De Angelis F. Intermolecular Interactions in Dye-Sensitized Solar Cells: A Computational Modeling Perspective. J Phys Chem Lett 2013; 4:956-74. [PMID: 26291363 DOI: 10.1021/jz302147v] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
We present a unified overview of our recent activity on the modeling of relevant intermolecular interactions occurring in dye-sensitized solar cells (DSCs). The DSC is an inherent complex system, whose efficiency is essentially determined by the interrelated phenomena occurring at the multiple molecular-semiconductor-electrolyte heterointerfaces. In this Perspective, we illustrate the basic methodology and selected applications of computational modeling of dye-dye and dye-coadsorbent intermolecular interactions taking place at the dye-sensitized interface. We show that the proposed methodology offers a realistic picture of aggregation phenomena among surface-adsorbed dyes and nicely describes semiconductor surfaces cosensitized by different dyes. The information acquired from this type of studies might constitute the basis for an integrated multiscale computational description of the device functioning, including all of the possible interdependencies among the device constituents, which may further boost the DSCs efficiency. We believe that this direction should be the target of future computational research in the DSC field.
Collapse
Affiliation(s)
- Mariachiara Pastore
- Computational Laboratory for Hybrid Organic Photovoltaics (CLHYO), Istituto CNR di Scienze e Tecnologie Molecolari, via Elce di Sotto 8, I-06123 Perugia, Italy
| | - Filippo De Angelis
- Computational Laboratory for Hybrid Organic Photovoltaics (CLHYO), Istituto CNR di Scienze e Tecnologie Molecolari, via Elce di Sotto 8, I-06123 Perugia, Italy
| |
Collapse
|
38
|
Beyhan SM, Götz AW, Visscher L. Bond energy decomposition analysis for subsystem density functional theory. J Chem Phys 2013; 138:094113. [DOI: 10.1063/1.4793629] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
39
|
Abstract
Background Among the 20 natural amino acids histidine is the most active and versatile member that plays the multiple roles in protein interactions, often the key residue in enzyme catalytic reactions. A theoretical and comprehensive study on the structural features and interaction properties of histidine is certainly helpful. Results Four interaction types of histidine are quantitatively calculated, including: (1) Cation-π interactions, in which the histidine acts as the aromatic π-motif in neutral form (His), or plays the cation role in protonated form (His+); (2) π-π stacking interactions between histidine and other aromatic amino acids; (3) Hydrogen-π interactions between histidine and other aromatic amino acids; (4) Coordinate interactions between histidine and metallic cations. The energies of π-π stacking interactions and hydrogen-π interactions are calculated using CCSD/6-31+G(d,p). The energies of cation-π interactions and coordinate interactions are calculated using B3LYP/6-31+G(d,p) method and adjusted by empirical method for dispersion energy. Conclusions The coordinate interactions between histidine and metallic cations are the strongest one acting in broad range, followed by the cation-π, hydrogen-π, and π-π stacking interactions. When the histidine is in neutral form, the cation-π interactions are attractive; when it is protonated (His+), the interactions turn to repulsive. The two protonation forms (and pKa values) of histidine are reversibly switched by the attractive and repulsive cation-π interactions. In proteins the π-π stacking interaction between neutral histidine and aromatic amino acids (Phe, Tyr, Trp) are in the range from -3.0 to -4.0 kcal/mol, significantly larger than the van der Waals energies.
Collapse
|
40
|
Klimeš J, Michaelides A. Perspective: Advances and challenges in treating van der Waals dispersion forces in density functional theory. J Chem Phys 2013; 137:120901. [PMID: 23020317 DOI: 10.1063/1.4754130] [Citation(s) in RCA: 595] [Impact Index Per Article: 54.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Abstract
Electron dispersion forces play a crucial role in determining the structure and properties of biomolecules, molecular crystals, and many other systems. However, an accurate description of dispersion is highly challenging, with the most widely used electronic structure technique, density functional theory (DFT), failing to describe them with standard approximations. Therefore, applications of DFT to systems where dispersion is important have traditionally been of questionable accuracy. However, the last decade has seen a surge of enthusiasm in the DFT community to tackle this problem and in so-doing to extend the applicability of DFT-based methods. Here we discuss, classify, and evaluate some of the promising schemes to emerge in recent years. A brief perspective on the outstanding issues that remain to be resolved and some directions for future research are also provided.
Collapse
Affiliation(s)
- Jirí Klimeš
- Thomas Young Centre, London Centre for Nanotechnology and Department of Chemistry, University College London, London WC1E 6BT, United Kingdom
| | | |
Collapse
|
41
|
Smart P, Mason CA, Loader JR, Meijer AJHM, Florence AJ, Shankland K, Fletcher AJ, Thompson SP, Brunelli M, Hill AH, Brammer L. Zipping and Unzipping of a Paddlewheel Metal-Organic Framework to Enable Two-Step Synthetic and Structural Transformation. Chemistry 2013; 19:3552-7. [DOI: 10.1002/chem.201204492] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2012] [Indexed: 11/05/2022]
|
42
|
Isegawa M, Fiedler L, Leverentz HR, Wang Y, Nachimuthu S, Gao J, Truhlar DG. Polarized Molecular Orbital Model Chemistry 3. The PMO Method Extended to Organic Chemistry. J Chem Theory Comput 2013; 9:33-45. [PMID: 23704835 PMCID: PMC3658842 DOI: 10.1021/ct300509d] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The polarized molecular orbital (PMO) method, a neglect-of-diatomic-differential-overlap (NDDO) semiempirical molecular orbital method previously parameterized for systems composed of O and H, is here extended to carbon. We modified the formalism and optimized all the parameters in the PMO Hamiltonian by using a genetic algorithm and a database containing both electrostatic and energetic properties; the new parameter set is called PMO2. The quality of the resulting predictions is compared to results obtained by previous NDDO semiempirical molecular orbital methods, both including and excluding dispersion terms. We also compare the PMO2 properties to SCC-DFTB calculations. Within the class of semiempirical molecular orbital methods, the PMO2 method is found to be especially accurate for polarizabilities, atomization energies, proton transfer energies, noncovalent complexation energies, and chemical reaction barrier heights and to have good across-the-board accuracy for a range of other properties, including dipole moments, partial atomic charges, and molecular geometries.
Collapse
Affiliation(s)
- Miho Isegawa
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Luke Fiedler
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Hannah R. Leverentz
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Yingjie Wang
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Santhanamoorthi Nachimuthu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Jiali Gao
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
43
|
Zahn S, MacFarlane DR, Izgorodina EI. Assessment of Kohn–Sham density functional theory and Møller–Plesset perturbation theory for ionic liquids. Phys Chem Chem Phys 2013; 15:13664-75. [DOI: 10.1039/c3cp51682b] [Citation(s) in RCA: 90] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
44
|
Pastore M, De Angelis F. Modeling Materials and Processes in Dye-Sensitized Solar Cells: Understanding the Mechanism, Improving the Efficiency. Top Curr Chem (Cham) 2013; 352:151-236. [DOI: 10.1007/128_2013_468] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
|
45
|
Recent Progress in Density Functional Methodology for Biomolecular Modeling. STRUCTURE AND BONDING 2013. [DOI: 10.1007/978-3-642-32750-6_1] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
|
46
|
Milione S, Milano G, Cavallo L. Pentacoordinated Organoaluminum Complexes: A Computational Insight. Organometallics 2012. [DOI: 10.1021/om300689g] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Stefano Milione
- Dipartimento di Chimica e Biologia, Università di Salerno, via Ponte don Melillo.
I-84084 Fisciano (SA), Italy
| | - Giuseppe Milano
- Dipartimento di Chimica e Biologia, Università di Salerno, via Ponte don Melillo.
I-84084 Fisciano (SA), Italy
| | - Luigi Cavallo
- Kaust Catalysis Center, Physiscal
Sciences and Engineering Division, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900.
Saudi Arabia
| |
Collapse
|
47
|
Lin YS, Li GD, Mao SP, Chai JD. Long-Range Corrected Hybrid Density Functionals with Improved Dispersion Corrections. J Chem Theory Comput 2012; 9:263-72. [PMID: 26589028 DOI: 10.1021/ct300715s] [Citation(s) in RCA: 399] [Impact Index Per Article: 33.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
By incorporating the improved empirical atom-atom dispersion corrections from DFT-D3 [Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. J. Chem. Phys.2010, 132, 154104], two long-range corrected (LC) hybrid density functionals are proposed. Our resulting LC hybrid functionals, ωM06-D3 and ωB97X-D3, are shown to be accurate for a very wide range of applications, such as thermochemistry, kinetics, noncovalent interactions, frontier orbital energies, fundamental gaps, and long-range charge-transfer excitations, when compared with common global and LC hybrid functionals. Relative to ωB97X-D [Chai, J.-D.; Head-Gordon, M. Phys. Chem. Chem. Phys.2008, 10, 6615], ωB97X-D3 (reoptimization of ωB97X-D with improved dispersion corrections) is shown to be superior for nonbonded interactions, and similar in performance for bonded interactions, while ωM06-D3 is shown to be superior for general applications.
Collapse
Affiliation(s)
- You-Sheng Lin
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Guan-De Li
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Shan-Ping Mao
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Jeng-Da Chai
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan.,Center for Theoretical Sciences and Center for Quantum Science and Engineering, National Taiwan University, Taipei 10617, Taiwan
| |
Collapse
|
48
|
Factors that distort the heme structure in Heme-Nitric Oxide/OXygen-Binding (H-NOX) protein domains. A theoretical study. J Inorg Biochem 2012; 118:28-38. [PMID: 23123336 DOI: 10.1016/j.jinorgbio.2012.09.011] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2012] [Revised: 09/05/2012] [Accepted: 09/07/2012] [Indexed: 11/21/2022]
Abstract
DFT and dispersion-corrected DFT calculations were carried out to probe the factors that distort the heme structure in Heme-Nitric oxide/OXygen-binding (H-NOX) protein domains. Various model systems that include heme, heme+surrounding residues, and heme+surrounding residues+additional protein environment were examined; the latter system was calculated with a quantum mechanics/molecular mechanics (QM/MM) method. The computations were extended to a myoglobin (Mb) protein, in which the heme structure is quite planar, in contrast to that in H-NOX. The natural tendency of the heme is to be planar. The strong structural distortion in H-NOX is mainly brought about by the intermolecular interactions between the whole heme molecule (heme ring plus its peripheral substituents) and the surrounding residues, among which the polar residues (Tyr140, Pro115, Mse98) play major roles in distorting the heme structure. The two peripheral propionate substituents that are oriented on the same side of the heme plane can also make the molecule distort, but the distortion caused by this factor is not significant. In Mb, the surrounding residues considered are all nonpolar and do not cause a structural distortion. The different structural features of the heme macrocycle in the different proteins (H-NOX and Mb) are reproduced by the calculations. The dispersion correction is necessary, since it improves the calculated structures. The effects of the distortion on the binding affinity of the axial ligand to the heme were also examined.
Collapse
|
49
|
Austin A, Petersson GA, Frisch MJ, Dobek FJ, Scalmani G, Throssell K. A Density Functional with Spherical Atom Dispersion Terms. J Chem Theory Comput 2012; 8:4989-5007. [PMID: 26593191 DOI: 10.1021/ct300778e] [Citation(s) in RCA: 353] [Impact Index Per Article: 29.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A new hybrid density functional, APF, is introduced, which avoids the spurious long-range attractive or repulsive interactions that are found in most density functional theory (DFT) models. It therefore provides a sound baseline for the addition of an empirical dispersion correction term, which is developed from a spherical atom model (SAM). The APF-D empirical dispersion model contains nine adjustable parameters that were selected based on a very small training set (15 noble gas dimers and 4 small hydrocarbon dimers), along with two computed atomic properties (ionization potential and effective atomic polarizability) for each element. APF-D accurately describes a large portion of the potential energy surfaces of complexes of noble gas atoms with various diatomic molecules involving a wide range of elements and of dimers of small hydrocarbons, and it reproduces the relative conformational energies of organic molecules. The accuracy for these weak interactions is comparable to that of CCSD(T)/aug-cc-pVTZ calculations. The accuracy in predicting the geometry of hydrogen bond complexes is competitive with other models involving DFT and empirical dispersion.
Collapse
Affiliation(s)
- Amy Austin
- Hall-Atwater Laboratories of Chemistry, Wesleyan University , Middletown, Connecticut 06459, United States
| | - George A Petersson
- Hall-Atwater Laboratories of Chemistry, Wesleyan University , Middletown, Connecticut 06459, United States
| | - Michael J Frisch
- Hall-Atwater Laboratories of Chemistry, Wesleyan University , Middletown, Connecticut 06459, United States.,Gaussian, Inc. , 340 Quinnipiac Street Building 40, Wallingford, Connecticut 06492, United States
| | - Frank J Dobek
- Hall-Atwater Laboratories of Chemistry, Wesleyan University , Middletown, Connecticut 06459, United States
| | - Giovanni Scalmani
- Gaussian, Inc. , 340 Quinnipiac Street Building 40, Wallingford, Connecticut 06492, United States
| | - Kyle Throssell
- Hall-Atwater Laboratories of Chemistry, Wesleyan University , Middletown, Connecticut 06459, United States
| |
Collapse
|
50
|
Dappe YJ, Bolcatto PG, Ortega J, Flores F. Dynamical screening of the van der Waals interaction between graphene layers. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2012; 24:424208. [PMID: 23032606 DOI: 10.1088/0953-8984/24/42/424208] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
The interaction between graphene layers is analyzed combining local orbital DFT and second order perturbation theory. For this purpose we use the linear combination of atomic orbitals-orbital occupancy (LCAO-OO) formalism, that allows us to separate the interaction energy as the sum of a weak chemical interaction between graphene layers plus the van der Waals interaction (Dappe et al 2006 Phys. Rev. B 74 205434). In this work, the weak chemical interaction is calculated by means of corrected-LDA calculations using an atomic-like sp(3)d(5) basis set. The van der Waals interaction is calculated by means of second order perturbation theory using an atom-atom interaction approximation and the atomic-like-orbital occupancies. We also analyze the effect of dynamical screening in the van der Waals interaction using a simple model. We find that this dynamical screening reduces by 40% the van der Waals interaction. Taking this effect into account, we obtain a graphene-graphene interaction energy of 70 ± 5 meV/atom in reasonable agreement with the experimental evidence.
Collapse
Affiliation(s)
- Y J Dappe
- Service de Physique de l'Etat Condensé, CEA Saclay, DSM/IRAMIS/SPEC, URA CNRS 2464, F-91191 Gif-Sur-Yvette Cedex, France
| | | | | | | |
Collapse
|