Dispersion Interactions in Density-Functional Theory
|
|
- Dulcie Hampton
- 5 years ago
- Views:
Transcription
1 Dispersion Interactions in Density-Functional Theory Erin R. Johnson and Alberto Otero-de-la-Roza Chemistry and Chemical Biology, University of California, Merced E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 1 / 37
2 Dispersion interactions Biological molecules. Surface adsorption. Molecular crystal packing. Crystal structure prediction. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 2 / 37
3 Dispersion interactions Long-range non-local correlation effect, not captured by semi-local functionals. Johnson et al. CPL 394 (2004) 334. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 3 / 37
4 Review of dispersion methods Fluctuation-dissipation theorem/rpa. Non-local correlation functionals (vdw-dfx). Dispersion correcting potentials (DCP). Meta-GGAs fit to binding energies. Post-SCF pairwise energy, fixed C 6 : DFT-D2. Variable C 6 : Tkatchenko-Scheffler, DFT-D3, XDM. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 4 / 37
5 The XDM method Dispersion arises from interaction of induced dipoles. The source of the instantaneous dipole moments is taken to be the dipole moment of the exchange hole. Becke and Johnson JCP 127 (2007) E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 5 / 37
6 The exchange hole The exchange hole measures the depletion in probability of finding another same-spin electron in the vicinity of a reference electron. dipole reference electron nucleus hole center An electron plus its exchange hole has zero total charge, but a non-zero dipole moment in general. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 6 / 37
7 The exchange-hole model The magnitude d X of the exchange-hole dipole moment is obtained using the Becke-Roussel exchange-hole model; PRA 39 (1989) Ae -ar reference point hole center b Parameters (A,a,b) obtained from normalization, density, and curvature at reference point. Advantages: semi-local (meta-gga) model of the dipole: d x = b. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 7 / 37
8 The XDM method The dispersion energy comes from second-order perturbation theory E (2) = ˆV 2 int E V int (r A, r B ) = multipole moments of electron + hole at r A interacting with multipole moments of electron + hole at r B E is the average excitation energy, obtained from second-order pertubation theory applied to polarizability. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 8 / 37
9 The XDM equations The XDM dispersion energy is: E disp = 1 2 ij C 6 f 6 (R ij ) R 6 ij + C 8f 8 (R ij ) R 8 ij The dispersion coefficients are non-empirical. + C 10f 10 (R ij ) R 10 ij +... C 6,ij = α iα j M 2 1 M2 1 j M 2 1 α j + M 2 1 jα i Atomic multipole moment integrals use Hirshfeld atomic partitioning. Ml 2 i = ω i (r)ρ σ (r)[ri l (r i d Xσ ) l ] 2 dr σ E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 9 / 37
10 Damping function Corrects for the multipolar-expansion error and avoids discontinuities. f n (R) = R n R n + R n vdw R vdw = a 1 R c,ij + a 2 a 1 and a 2 are parameters fit for use with a particular XC functional. R c,ij are proportional to atomic volumes. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 10 / 37
11 Higher-order dispersion coefficients XDM is easily extended to n-body dispersion coefficients to any order. The pairwise dispersion coefficients are: C2n+2 AB = 1 n/2 λ ε A + ε k(n k) Mk 2 A Mn k 2 B B k=1 and the average excitation energies are: ε A = 2 M1 2 A 3 αa 0 where αa 0 is the atom-in-molecule static polarizability. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 11 / 37
12 Generalization to XDM to higher orders The 3-body dispersion coefficients are: Z la l B l C = ξ la l B l C M 2 l A A M 2 l B B M 2 l C C (ε A + ε B + ε C ) (ε A + ε B )(ε A + ε C )(ε B + ε C ) The triple-dipole Axilrod-Teller-Muto term has coefficient C 9 = Z 111 : E (3) 9 = C 9 3 cos θ A cos θ B cos θ C + 1 R 3 AB R3 AC R3 BC and involves a geometrical factor depending on the atomic positions. Otero-de-la-Roza and Johnson, JCP 138 (2013) E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 12 / 37
13 3-body damping function The choice of 3-body damping function is not unique. Alternatives are: f (3) 9 (R AB, R AC, R BC ) = f 3 (R AB )f 3 (R AC )f 3 (R BC ) f (3) 9 (R AB, R AC, R BC ) = f 6 (R AB )f 6 (R AC )f 6 (R BC ) The latter damping is chosen to recover the united-atom limit. It gives improved agreement with SAPT for benzene and noble gas crystals. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 13 / 37
14 Effects on molecular dimer binding energies MAPE (%) E disp = 1 2 n ij C n f n (R ij ) R n ij no C 9 with C n C 9 and pairwise terms beyond C 10 do not improve performance. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 14 / 37
15 Implementation for molecules Pair XDM with GGA, hybrid, and range-separated functionals. XDM calculations used Gaussian 09 and the postg program. From the wfn file, postg gives: XDM dispersion energy forces for geometry optimization (fixed coefficients) second derivatives for frequencies Download postg from E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 15 / 37
16 Parametrization set 49 gas-phase dimers from Kannemann and Becke; JCTC 6 (2010) dispersion π-stacking dipole - induced dipole mixed dipole - dipole hydrogen-bonding E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 16 / 37
17 Binding energies XDM with aug-cc-pvtz; mean absolute errors in kcal/mol. Pure functionals: Quantity BLYP PW86 PBE MAE MA%E Hybrid and range-separated functionals: Quantity B3LYP BH&HLYP PBE0 CAM-B3LYP LC-ωPBE MAE MA%E E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 17 / 37
18 Role of exchange The exact exchange potential decays as 1/r far from a molecule. In terms of the exchange hole, h X remains on the molecule as the reference point moves away from it. The 1/r asymptotic dependence was used to design the B88 exchange functional. Functionals based on B88 or range-separated hybrids with the full exact-exchange limit (LC-ωPBE) give more accurate intermolecular exchange contributions. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 18 / 37
19 Thermochemistry Mean absolute errors in XDM thermochemical benchmarks. (kcal/mol) BLYP B3LYP LC-ωPBE G3/99 atomization energies Base functional With XDM Main group bond energies Base functional With XDM Charge-transfer complexes Base functional With XDM JCP 138 (2013) E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 19 / 37
20 Implementation for solids Quantum ESPRESSO, plane waves/psuedopotentials d xσ, valence τ, ρ. ω i, all-electron ρ, ρ at. E disp fast compared to E DFT. Errors in calculated XDM binding energies (kcal/mol) PBE PW86 B86b BLYP MAE MA%E JCP 136 (2012) E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 20 / 37
21 Graphite E ex (mev/atom) c (Å) B86b PW86 PBE revpbe Expt. Expt. ACFDT RPA LDA GGA+vdw QMC VDW DF E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 21 / 37
22 Sublimation enthalpies E crys el Eel mol H sub (V, T) = Eel mol + E trans + E rot + Evib mol + pv ( E crys el + E crys ) vib DFT + dispersion DFT + dispersion, supercell E trans + E rot + pv 4RT (7/2RT) Rigid molecules: E mol vib = Ecrys vib for intramolecular modes Dulong-Petit limit: E crys vib Zero-point vibrational contributions neglected 6RT (5RT) for intramolecular modes E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 22 / 37
23 Sublimation enthalpies Mean absolute errors in sublimation enthalpies, in kj/mol Quantity XDM DFT-D2 TS vdw-df B86b PW86 PBE PBE PBE v1 v2 MAE MA%E molecular crystals, low polymorphism, varied interaction types. Well known sublimation enthalpies at or below room temperature. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 23 / 37
24 Crystal structures Vibrational Helmholtz free energy: Thermal pressure: 3n [ ωj F vib (V, T) = 2 + k BT ln (1 )] e ω j/k B T p th = F vib V Equilibrium condition: j= E V = p th = p sta 0 Relax the crystal under negative pressure, p th DOS (1/cm -1 ) Frequency (cm -1 ) E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 24 / 37
25 Crystal structures Mean absolute errors in cell lengths, in atomic units. XDM DFT-D2 TS vdw-df Quantity B86b PW86 PBE PBE PBE v1 v2 MAE MA%E JCP 137 (2012) E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 25 / 37
26 Chiral crystals E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 26 / 37
27 Enantiomeric excess Enantiomers have the same solvation energies. Assume same temperature effects. Predicted ee: E = E dl E l ee = β2 1 β β = e E/RT Enantiomeric excess values (%): Amino acid DFT Expt. Serine Histidine Leucine Alanine Cysteine Tyrosine Valine Proline Aspartic acid Glutamic acid E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 27 / 37
28 Enantiomeric excess 1 Enantiomeric Excess E (kcal/mol) E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 28 / 37
29 Electrides An electride is an ionic substance in which a localized electron acts as an anion. Existing electrides require a cage like structure to stabilise the cation: crown ethers and cryptands. High magnetic susceptibilities, variable conductivities, very strong reducing agents. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 29 / 37
30 NCI analysis Reduced density gradient: s = 1 ρ 2(3π 2 ) 1/3 ρ 4/3 Maps intermolecular contacts (i.e. Pauli repulsion zones). Simpler to calculate than Bader s topology Color-mapping and region extent: interaction strength. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 30 / 37
31 NCI analysis Implemented in nciplot for molecules and critic2 for solids. JACS 132 (2010) 6498 and PCCP 14 (2012) E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 31 / 37
32 Electrides Use the NCI index to visualize the electrons - non-nuclear attractors will have low density and reduced gradient. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 32 / 37
33 Electrides J. L. Dye used van der Waals radii to generate approximate channels and vacancies of electrons - JACS 1996, 118 (1996) E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 33 / 37
34 Introduction XDM Molecules Solids Summary Graphite surfaces 2 2 f p p s h e d. o e - Applications FIG. 2. The two-layer finite DFT model for the step edge is shown, together with the methane molecule representing the tip. graphite two-dimensional cells (containing two carbon atoms) were enough to represent a step. However, in order to spurious self-interactions between adjacent FIG. 1. avoid Illustration of the molecular dynamics simulation of molecules, at tip least 4 cells required the amethane 5.6 nm diameter AFM sliding overwere graphite with in a step edge. Details are given in the exceeding text. complementary direction, our computational FIG. 2. The two-layer finite DFT model for the step edge is resources at present. Because of this limitation, we deshown, together with the methane molecule representing the cided to use a finite model. Several finite graphite steps R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 34 / 37 tip. weree.tested and we decided to use a model containing two
35 Introduction XDM Molecules Solids Applications Summary Graphite surfaces Appl. Phys. Lett. In press. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 35 / 37
36 XDM summary XDM dispersion coefficients are obtained from first-principles. Computational cost is comparable to semi-local DFT. Excellent binding energies. Improved thermochemistry. Applications for molecules and the solid state. E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 36 / 37
37 Acknolwdgements Group members: Alberto Otero-de-la-Roza Matthew Christian Stephen Dale Joseph Dizon Joel Mallory E. R. Johnson (UC Merced) Dispersion from XDM GRC: TD-DFT (Aug 2013) 37 / 37
Dispersion Interactions from the Exchange-Hole Dipole Moment
Dispersion Interactions from the Exchange-Hole Dipole Moment Erin R. Johnson and Alberto Otero-de-la-Roza Chemistry and Chemical Biology, University of California, Merced E. R. Johnson (UC Merced) Dispersion
More informationIntermolecular Forces in Density Functional Theory
Intermolecular Forces in Density Functional Theory Problems of DFT Peter Pulay at WATOC2005: There are 3 problems with DFT 1. Accuracy does not converge 2. Spin states of open shell systems often incorrect
More informationMethods for van der Waals Interactions
Methods for van der Waals Interactions Alexandre Tkatchenko Theory Department, Fritz Haber Institut der MPG Berlin, Germany tkatchen@fhi berlin.mpg.de Haber Institute FHI DFT and Beyond Workshop, Jul.
More informationNon-covalent force fields computed ab initio
Non-covalent force fields computed ab initio Supermolecule calculations Symmetry-adapted perturbation theory (SAPT) Supermolecule calculations Requirements: E = E AB E A E B. Include electron correlation,
More informationDensity driven delocalisation error
Density driven delocalisation error Stephen G. Dale, Alberto Otero-de-la-Roza, Erin R. Johnson Department of Chemistry and Chemical Biology, University of California, Merced, 5200 North Lake Road, Merced,
More informationIncluding dispersive interactions in DFT calculations
04/09/18 1 MSSC2018 Ab initio Modelling in Solid State Chemistry Torino, 2 7/09/2016 Including dispersive interactions in DFT calculations Bartolomeo Civalleri Department of Chemistry NIS Centre University
More informationAccurate van der Waals interactions from ground state electron density
Accurate van der Waals interactions from ground state electron density Alexandre Tkatchenko Theory Department, Fritz Haber Institut der MPG Berlin, Germany tkatchen@fhi berlin.mpg.de Haber Institute EXCITCM09,
More informationAtomic and molecular interaction forces in biology
Atomic and molecular interaction forces in biology 1 Outline Types of interactions relevant to biology Van der Waals interactions H-bond interactions Some properties of water Hydrophobic effect 2 Types
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationCombining density-functional theory and many-body methods
Combining density-functional theory and many-body methods Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France Vrije Universiteit Amsterdam, Netherlands November 2017 Outline 2/23 1
More informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
More informationApplications: Molecular crystals Graphite MgO(001)/CO MIL-53(Al) 2
Bartolomeo Civalleri Voice: Loredana Valenzano B3LYP augmented with an empirical dispersion term (B3LYP-D*) as applied to solids Università di Torino Dipartimento di Chimica IFM & NIS Torino - MSSC2009-10/09/09
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationDispersion Correcting Atom Centered Potentials (DCACP) and Many-Body Dispersion (MBD) contributions to interatomic vdw forces
Dispersion Correcting Atom Centered Potentials (DCACP) and Many-Body Dispersion (MBD) contributions to interatomic vdw forces O. Anatole von Lilienfeld Current: Argonne Leadership Computing Facility From
More informationMany-body dispersion interactions between semiconducting wires
Many-body dispersion interactions between semiconducting wires Alston J. Misquitta TCM Group, Cavendish Laboratory University of Cambridge 26 th July, 2010 Collaborators James Spencer (Imperial College,
More informationAdvanced Electronic Structure Theory Density functional theory. Dr Fred Manby
Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)
More informationJudging Density-Functional Approximations: Some Pitfalls of Statistics
Top Curr Chem (2015) 365: 81 96 DOI: 10.1007/128_2014_600 # Springer International Publishing Switzerland 2014 Published online: 24 February 2015 Judging Density-Functional Approximations: Some Pitfalls
More informationThe Electronic Structure of Dye- Sensitized TiO 2 Clusters from Many- Body Perturbation Theory
The Electronic Structure of Dye- Sensitized TiO 2 Clusters from Many- Body Perturbation Theory Noa Marom Center for Computational Materials Institute for Computational Engineering and Sciences The University
More informationDFT augmented with an empirical dispersion term as applied to solids Bartolomeo Civalleri
MSSC2011 - Ab initio Modelling in Solid State Chemistry Torino, 05/09/2011 1 DFT augmented with an empirical dispersion term as applied to solids Bartolomeo Civalleri Theoretical Chemistry Group Department
More informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationExtending Kohn-Sham density-functional theory
Extending Kohn-Sham density-functional theory Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France Email: julien.toulouse@upmc.fr Web page: www.lct.jussieu.fr/pagesperso/toulouse/ January
More informationReport on TS-vdW Method, and Code Development, and. Results
Michael Pawley July, 21, 2010 Report on TS-vdW Method, and Code Development, and Results Introduction: Non-covalent Interactions have long been difficult to account for using Density Functional Theory(DFT),
More informationDensity Functional Theory - II part
Density Functional Theory - II part antonino.polimeno@unipd.it Overview From theory to practice Implementation Functionals Local functionals Gradient Others From theory to practice From now on, if not
More informationElectric properties of molecules
Electric properties of molecules For a molecule in a uniform electric fielde the Hamiltonian has the form: Ĥ(E) = Ĥ + E ˆµ x where we assume that the field is directed along the x axis and ˆµ x is the
More informationSupporting information for Polymer interactions with Reduced Graphene Oxide: Van der Waals binding energies of Benzene on defected Graphene
Supporting information for Polymer interactions with Reduced Graphene Oxide: Van der Waals binding energies of Benzene on defected Graphene Mohamed Hassan, Michael Walter *,,, and Michael Moseler, Freiburg
More informationRole of van der Waals Interactions in Physics, Chemistry, and Biology
Role of van der Waals Interactions in Physics, Chemistry, and Biology How can we describe vdw forces in materials accurately? Failure of DFT Approximations for (Long-Range) Van der Waals Interactions 1
More informationWhen (2 + 2) 4? Alexandre Tkatchenko Fritz Haber Institut der Max Planck Gesellschaft, Berlin, Germany
When (2 + 2) 4? Alexandre Tkatchenko Fritz Haber Institut der Max Planck Gesellschaft, Berlin, Germany QMC@TTI, Apuan Alps, Jul 29, 2013 When (2 + 2) 4 or Van der Waals Interactions in Complex (and Simple)
More informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
More informationMolecular mechanics. classical description of molecules. Marcus Elstner and Tomáš Kubař. April 29, 2016
classical description of molecules April 29, 2016 Chemical bond Conceptual and chemical basis quantum effect solution of the SR numerically expensive (only small molecules can be treated) approximations
More informationThe broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals.
Physical Metallurgy The broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals. Crystal Binding In our discussions
More informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
More informationvan der Waals corrections to approximative DFT functionals
van der Waals corrections to approximative DFT functionals Ari P Seitsonen Physikalisch Chemisches Institut, Universität Zürich February 8, 2011 van der Waals interactions van der Waals interactions important
More informationEvaluating van der Waals energies from dielectric response functions
Evaluating van der Waals energies from dielectric response functions Deyu Lu Department of chemistry, UC Davis Electronic structure workshop, Davis, CA 06/23/2009 Motivation vdw forces play important roles
More informationSupplementary Figure 1(a) The trajectory of the levitated pyrolytic graphite test sample (blue curve) and
Supplementary Figure 1(a) The trajectory of the levitated pyrolytic graphite test sample (blue curve) and the envelope from free vibration theory (red curve). (b) The FFT of the displacement-time curve
More informationCORRECTING DENSITY FUNCTIONAL THEORY METHODS FOR DISPERSION INTERACTIONS USING PSEUDOPOTENTIALS
CORRECTING DENSITY FUNCTIONAL THEORY METHODS FOR DISPERSION INTERACTIONS USING PSEUDOPOTENTIALS by Ozan Karalti B.S. Chemistry, Bilkent University, 2002 Submitted to the Graduate Faculty of the Kenneth
More informationExchange-Correlation Functional
Exchange-Correlation Functional Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
More informationSUPPLEMENTARY INFORMATION
Calculations predict a stable molecular crystal of N 8 : Barak Hirshberg a, R. Benny Gerber a,b, and Anna I. Krylov c a Institute of Chemistry and The Fritz Haber Center for Molecular Dynamics, The Hebrew
More informationvan der Waals dispersion interactions in molecular materials: Beyond pairwise additivity
van der Waals dispersion interactions in molecular materials: Beyond pairwise additivity Anthony M. Reilly a, b and Alexandre Tkatchenko b, Received Xth XXXXXXXXXX 20XX, Accepted Xth XXXXXXXXX 20XX First
More informationAtomic structure & interatomic bonding. Chapter two
Atomic structure & interatomic bonding Chapter two 1 Atomic Structure Mass Charge Proton 1.67 х 10-27 kg + 1.60 х 10-19 C Neutron 1.67 х 10-27 kg Neutral Electron 9.11 х 10-31 kg - 1.60 х 10-19 C Electron
More informationCLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122
CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 191 THANKS TO MANY COLLABORATORS, INCLUDING SY VOSKO DAVID LANGRETH ALEX
More informationChapter 3. Crystal Binding
Chapter 3. Crystal Binding Energy of a crystal and crystal binding Cohesive energy of Molecular crystals Ionic crystals Metallic crystals Elasticity What causes matter to exist in three different forms?
More informationSubject of the Lecture:
Subject of the Lecture: Conceptual basis for the development of force fields. Implementation/validation Water - a worked example Extensions - combining molecular mechanics and quantum mechanics (QM/MM)
More informationAdvanced Quantum Chemistry III: Part 3. Haruyuki Nakano. Kyushu University
Advanced Quantum Chemistry III: Part 3 Haruyuki Nakano Kyushu University 2013 Winter Term 1. Hartree-Fock theory Density Functional Theory 2. Hohenberg-Kohn theorem 3. Kohn-Sham method 4. Exchange-correlation
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationQ-Chem 4.0: Expanding the Frontiers. Jing Kong Q-Chem Inc. Pittsburgh, PA
Q-Chem 4.0: Expanding the Frontiers Jing Kong Q-Chem Inc. Pittsburgh, PA Q-Chem: Profile Q-Chem is a high performance quantum chemistry program; Contributed by best quantum chemists from 40 universities
More informationQMC in the Apuan Alps IV Vallico Sotto, 28th July 2008
Van der Waals forces in graphitic nanostructures Norbert Nemec QMC in the Apuan Alps IV Vallico Sotto, 28th July 2008 1 History Johannes Diderik van der Waals (1837-1923): Postulation of intermolecular
More information5/27/2012. Role of van der Waals Interactions in Physics, Chemistry, and Biology
Role of van der Waals Interactions in Physics, Chemistry, and Biology 1 Role of van der Waals Interactions in Physics, Chemistry, and Biology Karin Jacobs: Take van der Waals forces into account in theory,
More informationDesign of Efficient Catalysts with Double Transition Metal. Atoms on C 2 N Layer
Supporting Information Design of Efficient Catalysts with Double Transition Metal Atoms on C 2 N Layer Xiyu Li, 1, Wenhui Zhong, 2, Peng Cui, 1 Jun Li, 1 Jun Jiang 1, * 1 Hefei National Laboratory for
More informationCHAPTER 2. Atomic Structure And Bonding 2-1
CHAPTER 2 Atomic Structure And Bonding 2-1 Structure of Atoms ATOM Basic Unit of an Element Diameter : 10 10 m. Neutrally Charged Nucleus Diameter : 10 14 m Accounts for almost all mass Positive Charge
More informationCanadian Journal of Chemistry. Van der Waals Potential Energy Surfaces from the Exchange-Hole Dipole Moment Dispersion Model
Van der Waals Potential Energy Surfaces from the Exchange-Hole Dipole Moment Dispersion Model Journal: Manuscript ID cjc-2016-021.r1 Manuscript Type: Article Date Submitted by the Author: 0-Jul-2016 Complete
More informationModule 6 1. Density functional theory
Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft
More informationDENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY
DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES
More informationVan der Waals interactions in DFT
Van der Waals interactions in DFT Maxime Dion*, Aaron Puzder*, T. Thonhauser,* Valentino R. Cooper*, Shen Li*, Eamonn Murray, Lingzhu Kong, Kyuho Lee, and David C. Langreth Department of Physics and Astronomy,
More information- intermolecular forces forces that exist between molecules
Chapter 11: Intermolecular Forces, Liquids, and Solids - intermolecular forces forces that exist between molecules 11.1 A Molecular Comparison of Liquids and Solids - gases - average kinetic energy of
More informationAtoms can form stable units called molecules by sharing electrons.
Atoms can form stable units called molecules by sharing electrons. The formation of molecules is the result of intramolecular bonding (within the molecule) e.g. ionic, covalent. Forces that cause the aggregation
More informationSolids, Liquids and Gases
WHY? Why is water usually a liquid and not a gas? Why does liquid water boil at such a high temperature for such a small molecule? Why does ice float on water? Why do snowflakes have 6 sides? Why is I
More informationThe change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:
The Born Energy of an Ion The free energy density of an electric field E arising from a charge is ½(ε 0 ε E 2 ) per unit volume Integrating the energy density of an ion over all of space = Born energy:
More informationWhy Is Molecular Interaction Important in Our Life
Why Is Molecular Interaction Important in ur Life QuLiS and Graduate School of Science iroshima University http://www.nabit.hiroshima-u.ac.jp/iwatasue/indexe.htm Suehiro Iwata Sept. 29, 2007 Department
More informationPublicat ion I. Publication I
Publicat ion I I Publication I K. J. Jalkanen, V. W. Jurgensen and I. M. Degtyarenko, Linear response properties required to simulate vibrational spectra of biomolecules in various media: (R)-Phenyloxirane
More informationIntermolecular and Intramolecular Forces. Introduction
Intermolecular and Intramolecular Forces Introduction Atoms can form stable units called molecules by sharing electrons. The formation of molecules is the result of intramolecular bonding (within the molecule)
More informationFundamental Interactions: 6 Forces
Fundamental Interactions: 6 Forces In nuclear and high-energy physics 6 fundamental forces are recognized, which describe the structure of matter. - the strong interaction - the weak interaction act inside
More informationSupplementary Information
Supplementary Information Supplementary Figure 1: After structural optimization of the CH 3 NH 3 PbI 3 unit cell, the resulting relaxed volumes for three distinct orientation of the MA molecules are shown.
More informationLecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions. Let's get together.
Lecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions Let's get together. Most gases are NOT ideal except at very low pressures: Z=1 for ideal gases Intermolecular interactions come
More informationBiomolecular modeling I
2015, December 15 Biomolecular simulation Elementary body atom Each atom x, y, z coordinates A protein is a set of coordinates. (Gromacs, A. P. Heiner) Usually one molecule/complex of interest (e.g. protein,
More informationChapter 11 Intermolecular Forces, Liquids, and Solids. Intermolecular Forces
Chapter 11, Liquids, and Solids States of Matter The fundamental difference between states of matter is the distance between particles. States of Matter Because in the solid and liquid states particles
More informationTutorial I: IQ MOL and Basic DFT and MP2 Calculations 1 / 30
Tutorial I: IQ MOL and Basic DFT and MP2 Calculations Q-Chem User Workshop, Denver March 21, 2015 1 / 30 2 / 30 Introduction to IQMOL DFT and MP2 Calculations 3 / 30 IQMOL and Q-CHEM IQMOL is an open-source
More informationMCB100A/Chem130 MidTerm Exam 2 April 4, 2013
MCBA/Chem Miderm Exam 2 April 4, 2 Name Student ID rue/false (2 points each).. he Boltzmann constant, k b sets the energy scale for observing energy microstates 2. Atoms with favorable electronic configurations
More informationLONG-RANGE VAN DER WAALS INTERACTION
International Journal of Modern Physics B Vol. 27, No. 18 (2013) 1330011 (31 pages) c World Scientific Publishing Company DOI: 10.1142/S0217979213300119 LONG-RANGE VAN DER WAALS INTERACTION JIANMIN TAO,
More informationMD simulation: output
Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square
More information2. As gas P increases and/or T is lowered, intermolecular forces become significant, and deviations from ideal gas laws occur (van der Waal equation).
A. Introduction. (Section 11.1) CHAPTER 11: STATES OF MATTER, LIQUIDS AND SOLIDS 1. Gases are easily treated mathematically because molecules behave independently. 2. As gas P increases and/or T is lowered,
More informationAdiabatic-connection fluctuation-dissipation density-functional theory based on range separation
Adiabatic-connection fluctuation-dissipation density-functional theory based on range separation Julien Toulouse 1 I. Gerber 2, G. Jansen 3, A. Savin 1, W. Zhu 1, J. Ángyán 4 1 Laboratoire de Chimie Théorique,
More informationModified Becke-Johnson (mbj) exchange potential
Modified Becke-Johnson (mbj) exchange potential Hideyuki Jippo Fujitsu Laboratories LTD. 2015.12.21-22 OpenMX developer s meeting @ Kobe Overview: mbj potential The semilocal exchange potential adding
More informationCombining density-functional theory and many-body methods
Combining density-functional theory and many-body methods Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France Albuquerque New Mexico, USA June 2016 Outline 2/23 1 A brief overview of
More information9/2/10 TYPES OF INTERMOLECULAR INTERACTIONS
Tro Chpt. 11 Liquids, solids and intermolecular forces Solids, liquids and gases - A Molecular Comparison Intermolecular forces Intermolecular forces in action: surface tension, viscosity and capillary
More informationRandom-phase approximation and beyond for materials: concepts, practice, and future perspectives. Xinguo Ren
Random-phase approximation and beyond for materials: concepts, practice, and future perspectives Xinguo Ren University of Science and Technology of China, Hefei USTC-FHI workshop on frontiers of Advanced
More informationCHAPTER ELEVEN KINETIC MOLECULAR THEORY OF LIQUIDS AND SOLIDS KINETIC MOLECULAR THEORY OF LIQUIDS AND SOLIDS
CHAPTER ELEVEN AND LIQUIDS AND SOLIDS KINETIC MOLECULAR THEORY OF LIQUIDS AND SOLIDS Differences between condensed states and gases? KINETIC MOLECULAR THEORY OF LIQUIDS AND SOLIDS Phase Homogeneous part
More informationSpring College on Computational Nanoscience
2145-29 Spring College on Computational Nanoscience 17-28 May 2010 At the Fifth Rung of Jacob's Ladder: A Discussion of Exact Exchange plus Local- and Nonlocal-density Approximations to the Correlation
More informationLondon Dispersion Forces (LDFs) Intermolecular Forces Attractions BETWEEN molecules. London Dispersion Forces (LDFs) London Dispersion Forces (LDFs)
LIQUIDS / SOLIDS / IMFs Intermolecular Forces (IMFs) Attractions BETWEEN molecules NOT within molecules NOT true bonds weaker attractions Represented by dashed lines Physical properties (melting points,
More informationNew applications of Diffusion Quantum Monte Carlo
New applications of Diffusion Quantum Monte Carlo Paul R. C. Kent (ORNL) Graphite: P. Ganesh, J. Kim, C. Park, M. Yoon, F. A. Reboredo (ORNL) Platinum: W. Parker, A. Benali, N. Romero (ANL), J. Greeley
More informationCHAPTER 11: INTERMOLECULAR FORCES AND LIQUIDS AND SOLIDS. Chemistry 1411 Joanna Sabey
CHAPTER 11: INTERMOLECULAR FORCES AND LIQUIDS AND SOLIDS Chemistry 1411 Joanna Sabey Forces Phase: homogeneous part of the system in contact with other parts of the system but separated from them by a
More informationLecture 2: Atom and Bonding Semester /2013
EMT 110 Engineering Materials Lecture 2: Atom and Bonding Semester 1 2012/2013 Atomic Structure Fundamental Concept Atoms are the structural unit of all engineering materials! Each atoms consist of nucleus
More informationAtoms, Molecules and Solids (selected topics)
Atoms, Molecules and Solids (selected topics) Part I: Electronic configurations and transitions Transitions between atomic states (Hydrogen atom) Transition probabilities are different depending on the
More informationSupporting Materials
Supporting Materials Figure S1 Experimental Setup Page Figure S (a) (b) (c) Feynman Diagrams Page 3-6 Figure S3 D IR Spectra Page 7 Figure S4 Kinetic Model Page 8 Figure S5 Van t Hoff Plots Page 9 1 k
More informationXYZ of ground-state DFT
XYZ of ground-state DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA January 5-9th, 2014 Kieron (UC Irvine) XYZ of ground-state
More informationIntermediate DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA
Intermediate DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA October 10-19th, 2012 Kieron (UC Irvine) Intermediate DFT Lausanne12
More information2. As gas P increases and/or T is lowered, intermolecular forces become significant, and deviations from ideal gas laws occur (van der Waal equation).
A. Introduction. (Section 11.1) CHAPTER 11: STATES OF MATTER, LIQUIDS AND SOLIDS 1. Gases are easily treated mathematically because molecules behave independently. 2. As gas P increases and/or T is lowered,
More informationSoot - Developing anisotropic potentials from first principles for PAH molecules. Tim Totton, Alston Misquitta and Markus Kraft 12/11/2009
Soot - Developing anisotropic potentials from first principles for PAH molecules. Tim Totton, Alston Misquitta and 12/11/2009 HRTEM images of soot Some evidence for different soot structures based on different
More informationImperfect Gases. NC State University
Chemistry 431 Lecture 3 Imperfect Gases NC State University The Compression Factor One way to represent the relationship between ideal and real gases is to plot the deviation from ideality as the gas is
More informationBasis Set for Molecular Orbital Theory
Basis Set for Molecular Orbital Theory! Different Types of Basis Functions! Different Types of Atom Center Basis Functions! Classifications of Gaussian Basis Sets! Pseudopotentials! Molecular Properties
More informationWalter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.
Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm
More informationOne-Electron Hamiltonians
One-Electron Hamiltonians Hartree-Fock and Density Func7onal Theory Christopher J. Cramer @ChemProfCramer 2017 MSSC, July 10, 2017 REVIEW A One-Slide Summary of Quantum Mechanics Fundamental Postulate:
More informationTowards gas-phase accuracy for condensed phase problems
Towards gas-phase accuracy for condensed phase problems Fred Manby Centre for Computational Chemistry, School of Chemistry University of Bristol STC 2006: Quantum Chemistry Methods and Applications Erkner,
More information510 Subject Index. Hamiltonian 33, 86, 88, 89 Hamilton operator 34, 164, 166
Subject Index Ab-initio calculation 24, 122, 161. 165 Acentric factor 279, 338 Activity absolute 258, 295 coefficient 7 definition 7 Atom 23 Atomic units 93 Avogadro number 5, 92 Axilrod-Teller-forces
More informationDensity Func,onal Theory (Chapter 6, Jensen)
Chem 580: DFT Density Func,onal Theory (Chapter 6, Jensen) Hohenberg- Kohn Theorem (Phys. Rev., 136,B864 (1964)): For molecules with a non degenerate ground state, the ground state molecular energy and
More informationLec20 Fri 3mar17
564-17 Lec20 Fri 3mar17 [PDF]GAUSSIAN 09W TUTORIAL www.molcalx.com.cn/wp-content/uploads/2015/01/gaussian09w_tutorial.pdf by A Tomberg - Cited by 8 - Related articles GAUSSIAN 09W TUTORIAL. AN INTRODUCTION
More informationOutline. Introduction: graphene. Adsorption on graphene: - Chemisorption - Physisorption. Summary
Outline Introduction: graphene Adsorption on graphene: - Chemisorption - Physisorption Summary 1 Electronic band structure: Electronic properties K Γ M v F = 10 6 ms -1 = c/300 massless Dirac particles!
More informationSolutions to Assignment #4 Getting Started with HyperChem
Solutions to Assignment #4 Getting Started with HyperChem 1. This first exercise is meant to familiarize you with the different methods for visualizing molecules available in HyperChem. (a) Create a molecule
More information