Non-covalent force fields computed ab initio
|
|
- Evelyn Small
- 5 years ago
- Views:
Transcription
1 Non-covalent force fields computed ab initio Supermolecule calculations Symmetry-adapted perturbation theory (SAPT)
2 Supermolecule calculations Requirements: E = E AB E A E B. Include electron correlation, intra- and inter-molecular (dispersion energy = intermolecular correlation). Choose good basis, with diffuse orbitals (and bond functions ) especially to converge the dispersion energy 3. Size consistency. Currently best method: CCSD(T) 4. Correct for basis set superposition error (BSSE) by computing E A and E B in dimer basis
3 Symmetry-adapted perturbation theory (SAPT) Combine perturbation theory with antisymmetrization A (Pauli) to include short-range exchange effects. Advantages:. E calculated directly.. Contributions (electrostatic, induction, dispersion, exchange) computed individually. Useful in analytic fits of potential surface. Advantage of supermolecule method: Easy, use any black-box molecular electronic structure program
4 Problems in SAPT:. Pauli: AH = HA. Antisymmetrizer commutes with total Hamiltonian H = H (0) + H (), but not with H (0) and H () separately. Has led to different definitions of second (and higher) order energies.. Free monomer wavefunctions Φ A k and Φ B k not exactly known. Use Hartree-Fock wave functions and apply double perturbation theory to include intra-molecular correlation, or use CCSD wave functions of monomers Many-body SAPT. Program packages: - SAPT for pair potentials - SAPT3 for 3-body interactions
5 Most difficult: dispersion interactions First ab initio calculation of He He binding curve: Phys. Rev. Letters, 5 (970) H.F. Schaefer, D.R. McLaughlin, F.E. Harris, and B.J. Alder page 988: D e =.0 K P. Bertoncini and A. C. Wahl page 99: D e =.4 K Present value: D e =.0 K = 7.65 cm 0. kj/mol 0 3 ev Hartree
6 Can one use DFT methods? Reviews by Johnson & DiLabio, Zhao & Truhlar Many different functionals tested with different basis sets type dimer mean error in D e Van der Waals Rg, (CH 4 ), (C H ), (benzene) % dipole-induced dipole CH 4 HF, H O benzene, etc % dipole-dipole (H CO), etc % hydrogen bonded (NH 3 ), (H O), (HCOOH), etc. 3-0% Some VdW dimers, (benzene) for example, not bound B97 best, B3LYP worst Often best results without BSSE correction, or smaller basis sets Right for wrong reason
7 Basic problems with DFT. Exchange repulsion Incorrect asymptotic behavior of one-electron potential: v(r) exp( αr) instead of /r In intermolecular Coulomb energy no self-term present self-exchange spurious attraction. Dispersion Intrinsically non-local: cannot be described by local LDA or semi-local GGA methods
8 DFT with dispersion explicitly included vdw-df: M. Dion, H. Rydberg, E. Schroder, D.C. Langreth, B.I. Lundqvist, Phys. Rev. Lett. 9 (004) Ar Ar interaction E int (cm ) SAPT(DFT)/PBE0 SAPT(DFT)/B97 CCSD(T)/CBS Benchmark vdw DF R (Å) From: R. Podeszwa and K. Szalewicz, Chem. Phys. Lett. 4 (005) 488
9 DFT with dispersion included semi-empirically Becke, Johnson, and others: - Include atom-atom C 6 R 6 C 8 R 8 C 0 R 0 term in addition to the semi-local correlation energy - Approximate C 6, C 8, C 0 from atomic squared multipole moment expectation values and (empirical) static polarizabilities - Empirically optimized damping parameters
10 Alternative: SAPT-DFT Implemented by G. Jansen et al. (Delaware) (Essen) and K. Szalewicz et al. First order SAPT energy (electrostatic + exchange) computed with monomer densities and density matrices from Kohn-Sham DFT Second-order SAPT energy (induction, dispersion + exchange) from (time-dependent) coupled perturbed Kohn-Sham response functions Only Hartree-Fock like expressions from Many-Body SAPT needed better scaling
11 Caution! SAPT-DFT requires XC potential that is good in inner region and has correct /r behavior for r Coupled time-dependent DFT must be used for (frequencydependent) density-density polarizabilities α(r, r, ω) Both groups, K. Szalewicz (Delaware) and G. Jansen (Essen), further improved efficiency by implementation of density fitting.
12 Ar Ar interaction From: R. Podeszwa and K. Szalewicz, Chem. Phys. Lett. 4 (005) E int (cm ) SAPT(DFT)/PBE0 SAPT(DFT)/B97 CCSD(T)/CBS Benchmark vdw DF R (Å)
13 A. Heßelmann, G. Jansen, M. Schütz, J. Chem. Phys. (005) 0403 (benzene) 5 GTOs (aug-cc-pvqz), extrapolation to basis set limit MP kj/mol CCSD(T) DF-SAPT-DFT standard DFT unbound metastable unbound
14 Adenine-Thymine (G. Jansen et al.) DF-SAPT-DFT up to aug-cc-pvqz level E [kj/mol] E [kj/mol] E () el E () el E () exch E () exch E () ind E () ind E () exch-ind E () exch-ind E () disp E () disp E () exch-disp E () exch-disp δ(hf) DFT-SAPT MP δ(hf) DFT-SAPT MP E int CCSD(T) E int CCSD(T)
15 Pair interactions in water trimer exch kcal/mole 0 0 disper elec induct total
16 Many-body interactions (per hydrogen bond) Trimer Tetramer pair 3-body total pair 3-body 4-body total Pentamer pair 3-body 4-body 5-body total kcal/mole
17 Illustration: New water potential Tested by spectroscopy on dimer and trimer Used in MD simulations for liquid water R. Bukowski, K. Szalewicz, G.C. Groenenboom, and A. van der Avoird, Science, 35, 49 (007)
18 New polarizable water pair potential: CC-pol From CCSD(T) calculations in aug-cc-pvtz + bond function basis Extrapolated to complete basis set (CBS) limit at MP level 50 carefully selected water dimer geometries Estimated uncertainty < 0.07 kcal/mol (same as best single-point calculations published)
19 CC-pol: Analytic representation Site-site model with 8 sites (5 symmetry distinct) per molecule Coulomb interaction, dispersion interaction, exponential overlap terms: first-order exchange repulsion, second-order exchange induction + dispersion Extra polarizability site for induction interaction Long range R n contributions computed by perturbation theory, subtracted before fit of short range terms Good compromise between accuracy of reproducing computed points (rmsd of 0.09 kcal/mol for E < 0) and simplicity needed for molecular simulations
20 How to test non-covalent force fields?
21 Molecular beam spectroscopy of Van der Waals molecules
22 W.L. Meerts, Molecular and Laser Physics, Nijmegen
23 Intermolecular potential Cluster (Van der Waals molecule) quantum levels, i.e., eigenstates of nuclear motion Hamiltonian Van der Waals spectra
24 Nuclear motion Hamiltonian H = T + V for normal (= semi-rigid) molecules single equilibrium structure small amplitude vibrations Use rigid rotor/harmonic oscillator model For (harmonic) vibrations Wilson GF -matrix method frequencies, normal coordinates Rigid rotor model fine structure (high resolution spectra)
25 Nuclear motion Hamiltonian H = T + V for weakly bound complexes (Van der Waals or hydrogen bonded) multiple equivalent equilibrium structures (= global minima in the potential surface V ) small barriers tunneling between minima large amplitude (VRT) motions: vibrations, internal rotations, tunneling (more or less rigid monomers) curvilinear coordinates complicated kinetic energy operator T
26 Method for molecule-molecule dimers H O H O, NH 3 NH 3, etc. Hamiltonian H = T A + T B h µr R R + (J j A j B ) µr + V (R, ω A, ω B ) Monomer Hamiltonians (X = A, B): T X = A X (j X ) a + B X (j X ) b + C X (j X ) c Basis for bound level calculations χ n (R) D (J) MK (α, β, 0) m A,m B D (j A) m A k A (ω A ) D (j B) m B k B (ω B ) j A, m A ; j B, m B j AB, K
27 The permutation-inversion (PI) symmetry group For semi-rigid molecules Use Point Group of Equilibrium Geometry to describe the (normal coordinate) vibrations N.B. This point group is isomorphic to the P I group, which contains all feasible permutations of identical nuclei, combined with inversion E. Molecule Point group P I group C v {E, E, (), () } C 3v {E, (3), (3), (), (3), (3) }
28 Example: H O P I operation frame rotation point group operation () R z (π) = C z C z E R y (π) = C y σ xz reflection () R x (π) = C x σ yz reflection permutation frame rotation + point group rotation permutation-inversion frame rotation + reflection Hence: P I-group point group
29 For floppy molecules/complexes multiple equivalent minima in V low barriers: tunneling between these minima is feasible. additional feasible P I-operations Example NH 3 semi-rigid NH 3 P I(C 3v ) = {E, (3), (3), (), (3), (3) } + inversion tunneling (umbrella up down) P I(D 3h ) = P I(C 3v ) {E, E } Also (), (3), (3) and E are feasible Additional feasible P I-operations observable tunneling splittings in spectrum
30 For H O H O PI group G 6 = {E, P } {E, P 34 } {E, P AB } {E, E } Equilibrium geometry has C s symmetry 8-fold tunneling splitting of rovib levels
31 Water cluster spectra (far-infrared, high-resolution) from Saykally group (UC Berkeley) Used for test of potential: Pair potential Dimer VRT levels Dimer spectrum Pair + 3-body potential Trimer VRT levels Trimer spectrum
32 Water dimer tunneling pathways Acceptor Tunneling Donor-Acceptor Interchange Tunneling Donor (Bifurcation) Tunneling
33 B E Water dimer A tunneling levels (J = K = 0) B + E + A + Acceptor tunneling Interchange tunneling Bifurcation tunneling
34 0 5 (H O) levels computed from various potentials (J=K=0) B R. S. Fellers et al., J. Chem. Phys. 0, 6306 (999) E A B E A cm 0 B E A B E A B E 5 0 B + E + + A A + B E + + A B + E + + A B + E + + A B + E + + A Experiment ASP W ASP S NEMO3 MCY
35 Acceptor tunneling a(k=0) + a(k=) = 6.58 (ab initio) 3.9 (experiment) Rotational constants B,A E +,E + A,B H O dimer A = (B+C)/ = B+C = B +,A E,E + A +,B tunneling levels from CC-pol (red) experiment (black) B E A B +,A A +,B B,A A,B a a E +,E E,E B + E + + A J=K=0 J=K= J=K=
36 H O dimer intermolecular vibrations from CC-pol (red) cm calculated experiment calculated experiment K = 0 K = 0 K = K = DT AT AW DT DT AT AW DT AT AW DT AW DT experiment (black) GS GS GS GS
37 Acceptor tunneling D O dimer a(k=0) a(k=) a(k=) (ab initio) (experiment) Rotational constants A = (B+C)/ = B+C = a(k=) B +,A E +,E + A,B B,A + E,E + + A,B tunneling levels from CC-pol (red) experiment (black) B E A a(k=) B +,A + A,B B +,A A +,B E +,E E,E + a(k=0) B + E + + A J=K=0 J=K= J=K=
38 40 calculated experiment calculated experiment K = 0 K = 0 K = K = D O dimer 0 DT DT DT intermolecular vibrations from CC-pol (red) experiment (black) cm AT AW DT AT AW DT AT AW DT AW DT GS GS GS GS
39 0 00 Second virial coefficient of water vapor B(T) [cm 3 /mol] Experiment CC pol SAPT 5s VRT(ASP W)III T [K]
40 Water trimer tunneling pathways 4 Torsional tunneling (flipping) 3 B 6 C 5 dud A uud Bifurcation (donor) tunneling
41 Water trimer tunneling levels (J = 0)
42 H O trimer torsional levels cm 00 k ± ± Experiment CC pol TTM. VRT(ASP W)III
43 D O trimer torsional levels cm k 00 ± (k=0) ± ± Experiment CC pol TTM. VRT(ASP W)III
44 MD simulations of liquid water, T = 98 K
45 3.5 CC-pol, B CC-pol+NB CC-pol+NB(ind) Experiment.6. g OH g OO.5 0 r [Angstroms] g HH r [Angstroms] r [Angstroms] Atom-atom radial distribution functions
46 Conclusions CC-pol pure ab initio Predicts dimer spectra better than semi-empirical potentials fitted to these spectra Second virial coefficients in excellent agreement with experiment CC-pol + 3-body potential gives good trimer spectrum Simulations of liquid water with CC-pol + N-body forces predict the neutron and X-ray diffraction data equally well as the best empirical potentials fitted to these data Important role of many-body forces in liquid water Nearly tetrahedral coordination: 3.8 hydrogen bonds, only.8 with pure pair potential
47 General conclusion CC-pol, with the accompanying many-body interaction model, provides the first water force field that recovers the dimer, trimer, and liquid properties well
48 Latest development Water pair potential with flexible monomers: -dimensional X. Huang, B. J. Braams, and J. M. Bowman, J. Phys. Chem. A 0 (006) 445 Improved version (J. Chem. Phys., accepted) From CCSD(T) calculations in aug-cc-pvtz basis water dimer geometries. Accurate fit in terms of symmetry-adapted polynomials of scaled atom-atom distances. Produces even more accurate VRT levels of (H O) and (D O). Also acceptor splitting a 0 + a accurate in d calculations. Good VRT levels of (H O) 3 and (D O) 3 when 3-body interactions added.
Noncovalent interactions between molecules. spectra of weakly bound molecular clusters
Noncovalent interactions between molecules and spectra of weakly bound molecular clusters Ad van der Avoird Theoretical Chemistry Outline: Preamble: Born-Oppenheimer approximation What are non-covalent
More informationInteraction potential for water dimer from symmetry-adapted perturbation theory based on density functional description of monomers
THE JOURNAL OF CHEMICAL PHYSICS 125, 044301 2006 Interaction potential for water dimer from symmetry-adapted perturbation theory based on density functional description of monomers Robert Bukowski and
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 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 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 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 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 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 informationEXAM INFORMATION. Radial Distribution Function: B is the normalization constant. d dx. p 2 Operator: Heisenberg Uncertainty Principle:
EXAM INFORMATION Radial Distribution Function: P() r RDF() r Br R() r B is the normalization constant., p Operator: p ^ d dx Heisenberg Uncertainty Principle: n ax n! Integrals: xe dx n1 a x p Particle
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 informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationDispersion 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 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 informationClose coupling results for inelastic collisions of NH3 and Ar. A stringent test of a spectroscopic potential
12 August 1994 ELSEVIER Chemical Physics Letters 226 ( 1994) 22-26 CHEMICAL PHYSICS LETTERS Close coupling results for inelastic collisions of NH3 and Ar. A stringent test of a spectroscopic potential
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
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 informationElectronic quantum effect on hydrogen bond geometry in. water dimer
Electronic quantum effect on hydrogen bond geometry in water dimer Danhui Li 1,2, Zhiyuan Zhang 1,2 Wanrun Jiang 1,2 Depeng Zhang 1,2 Yu Zhu 1,2 and Zhigang Wang 1,2* 1 Institute of Atomic and Molecular
More informationStructural Bioinformatics (C3210) Molecular Mechanics
Structural Bioinformatics (C3210) Molecular Mechanics How to Calculate Energies Calculation of molecular energies is of key importance in protein folding, molecular modelling etc. There are two main computational
More informationQMC dissociation energy of the water dimer: Time step errors and backflow calculations
QMC dissociation energy of the water dimer: Time step errors and backflow calculations Idoia G. de Gurtubay and Richard J. Needs TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay.
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 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 informationTowards a force field based on density fitting
THE JOURNAL OF CHEMICAL PHYSICS 124, 104101 2006 Towards a force field based on density fitting Jean-Philip Piquemal a and G. Andrés Cisneros Laboratory of Structural Biology, National Institute of Environmental
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationAb initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier
Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential
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 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 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 informationUptake of OH radical to aqueous aerosol: a computational study
Uptake of OH radical to aqueous aerosol: a computational study Grigory Andreev Karpov Institute of Physical Chemistry 10 Vorontsovo pole, Moscow, 105064, Russia Institute of Physical Chemistry and Electrochemistry
More informationCharacteristics of the interaction in azulene (H 2 X) n=1,2 (X=O,S) clusters.
Characteristics of the interaction in azulene (H 2 X) n=1,2 (X=O,S) clusters. Enrique M. Cabaleiro-Lago (a), Ángeles Peña-Gallego (b), Jesús Rodríguez-Otero (b), M. Merced Montero-Campillo (b) (a) Departamento
More informationReport on Atomistic Modeling of Bonding in Carbon-Based Nanostructures
Report on Atomistic Modeling of Bonding in Carbon-Based Nanostructures Timothy Stillings Department of Physics, Astronomy and Materials Science Missouri State University Advisor: Ridwan Sakidja Abstract
More informationMETHODS FOR TREATING SOLVENT EFFECTS AND INTERMOLECULAR FORCES. Mark S. Gordon Iowa State University Ames Laboratory
METHODS FOR TREATING SOLVENT EFFECTS AND INTERMOLECULAR FORCES Mark S. Gordon Iowa State University Ames Laboratory OUTLINE Solvation Methods Explicit vs. implicit methods Explicit Methods TIP3P, TIP4P
More informationXANES spectra of Bulk Water (Ph. Wernet, A. Nilsson, et al. Science 304, 995 (2004)) Pre-edge peak Post-edge
XANES spectra of Bulk Water (Ph. Wernet, A. Nilsson, et al. Science 304, 995 (2004)) Pre-edge peak Post-edge Transition State Potential Method Consider an electronic transition: Initial state Final state
More informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More informationElectron Correlation - Methods beyond Hartree-Fock
Electron Correlation - Methods beyond Hartree-Fock how to approach chemical accuracy Alexander A. Auer Max-Planck-Institute for Chemical Energy Conversion, Mülheim September 4, 2014 MMER Summerschool 2014
More informationA CSOV Study of the Difference between HF and DFT Intermolecular Interaction Energy Values: The Importance of the Charge Transfer Contribution
A CSOV Study of the Difference between HF and DFT Intermolecular Interaction Energy Values: The Importance of the Charge Transfer Contribution JEAN-PHILIP PIQUEMAL, 1 * ANTONIO MARQUEZ, 2 OLIVIER PARISEL,
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 informationMolecular Aggregation
Molecular Aggregation Structure Analysis and Molecular Simulation of Crystals and Liquids ANGELO GAVEZZOTTI University of Milano OXFORD UNIVERSITY PRESS Contents PART I FUNDAMENTALS 1 The molecule: structure,
More informationQuantum Chemical Simulations and Descriptors. Dr. Antonio Chana, Dr. Mosè Casalegno
Quantum Chemical Simulations and Descriptors Dr. Antonio Chana, Dr. Mosè Casalegno Classical Mechanics: basics It models real-world objects as point particles, objects with negligible size. The motion
More informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
More informationFeet on the potential energy surface, head in the π clouds
Graduate Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 2011 Feet on the potential energy surface, head in the π clouds Quentin Anthony Smith Iowa State University Follow
More informationClose agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations
Close agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations Alexandre V. Morozov, Tanja Kortemme, Kiril Tsemekhman, David Baker
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 informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationDispersion Interactions in Density-Functional Theory
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
More informationThe Schrödinger equation for many-electron systems
The Schrödinger equation for many-electron systems Ĥ!( x,, x ) = E!( x,, x ) 1 N 1 1 Z 1 Ĥ = " $ # " $ + $ 2 r 2 A j j A, j RAj i, j < i a linear differential equation in 4N variables (atomic units) (3
More informationChem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.
Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative
More informationOrbital dependent correlation potentials in ab initio density functional theory
Orbital dependent correlation potentials in ab initio density functional theory noniterative - one step - calculations Ireneusz Grabowski Institute of Physics Nicolaus Copernicus University Toruń, Poland
More informationChapter 2 Experimental sources of intermolecular potentials
Chapter 2 Experimental sources of intermolecular potentials 2.1 Overview thermodynamical properties: heat of vaporization (Trouton s rule) crystal structures ionic crystals rare gas solids physico-chemical
More informationThe far-infrared vibration rotation tunneling spectrum of the water tetramer-d8
The far-infrared vibration rotation tunneling spectrum of the water tetramer-d8 J. D. Cruzan, M. G. Brown, K. Liu, L. B. Braly, and R. J. Saykally Department of Chemistry, University of California at Berkeley,
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 informationHigh-level Quantum Chemistry Methods and Benchmark Datasets for Molecules
High-level Quantum Chemistry Methods and Benchmark Datasets for Molecules Markus Schneider Fritz Haber Institute of the MPS, Berlin, Germany École Polytechnique Fédérale de Lausanne, Switzerland دانشگاه
More informationUsing contracted basis functions and the Lanczos algorithm to compute spectra of flexible Van der Waals dimers
Using contracted basis functions and the Lanczos algorithm to compute spectra of flexible Van der Waals dimers Xiao-Gang Wang and Tucker Carrington Chemistry Department Queen s University 26 janvier 2016
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 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 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 informationComputational Chemistry. An Introduction to Molecular Dynamic Simulations
Computational Chemistry An Introduction to Molecular Dynamic Simulations Computational chemistry simulates chemical structures and reactions numerically, based in full or in part on the fundamental laws
More informationCCSD(T) benchmarks of non-equilibrium water clusters: the importance of monomer deformation
CCSD(T) benchmarks of non-equilibrium water clusters: the importance of monomer deformation Biswajit Santra 1, Angelos Michaelides 1,2, and Matthias Scheffler 1 1 Fritz-Haber-Institut der MPG, Berlin,
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 informationJoint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012
2327-3 Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data 23-27 January 2012 Qunatum Methods for Plasma-Facing Materials Alain ALLOUCHE Univ.de Provence, Lab.de la Phys.
More informationPrinciples of Quantum Mechanics
Principles of Quantum Mechanics - indistinguishability of particles: bosons & fermions bosons: total wavefunction is symmetric upon interchange of particle coordinates (space,spin) fermions: total wavefuncftion
More informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
More informationNew Perspective on structure and bonding in water using XAS and XRS
New Perspective on structure and bonding in water using XAS and XRS Anders Nilsson Stanford Synchrotron Radiation Laboratory (SSRL) and Stockholm University, Sweden R. Ludwig Angew. Chem. 40, 1808 (2001)
More informationAccurate Description of Intermolecular Interactions Involving Ions Using Symmetry-Adapted Perturbation Theory
pubs.acs.org/jctc Accurate Description of Intermolecular Interactions Involving Ions Using Symmetry-Adapted Perturbation Theory Ka Un Lao, Rainer Scha ffer, Georg Jansen,*, and John M. Herbert*, Department
More informationQUANTUM CHEMISTRY FOR TRANSITION METALS
QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline
More informationThe Monomer Electron Density Force Field. (MEDFF): A Physically Inspired Model for. Non-Covalent Interactions
The Monomer Electron Density Force Field (MEDFF): A Physically Inspired Model for Non-Covalent Interactions Steven Vandenbrande, Michel Waroquier, Veronique Van Speybroeck, and Toon Verstraelen Center
More informationTheoretical study of the water tetramer
Theoretical study of the water tetramer David J. Wales and Tiffany R. Walsh University Chemical Laboratories, Lensfield Road, Cambridge CB2 1EW, United Kingdom Received 18 December 1996; accepted 22 January
More informationQUANTUM CHEMISTRY PROJECT 3: PARTS B AND C
Chemistry 460 Fall 2017 Dr. Jean M. Standard November 6, 2017 QUANTUM CHEMISTRY PROJECT 3: PARTS B AND C PART B: POTENTIAL CURVE, SPECTROSCOPIC CONSTANTS, AND DISSOCIATION ENERGY OF DIATOMIC HYDROGEN (20
More informationWater dimer hydrogen bond stretch, donor torsion overtone, and in-plane bend vibrations
JOURNAL OF CHEMICAL PHYSICS VOLUM19, NUMBER 17 1 NOVEMBER 2003 Water dimer hydrogen bond stretch, donor torsion overtone, and in-plane bend vibrations Frank N. Keutsch Department of Chemistry and Chemical
More informationFragmentation methods
Fragmentation methods Scaling of QM Methods HF, DFT scale as N 4 MP2 scales as N 5 CC methods scale as N 7 What if we could freeze the value of N regardless of the size of the system? Then each method
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 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 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 informationDielectric polarization of 2-pyrrolidinone molecules in benzene solution - a quantum-chemical study
Dielectric polarization of 2-pyrrolidinone molecules in benzene solution - a quantum-chemical study L. Gorb* ), J. Jadżyn $) and K. W. Wojciechowski #) Institute of Molecular Physics, Polish Academy of
More informationCharacterization of methanol as a magnetic field tracer in star-forming regions
SUPPLEMENTARY INFORMATION Letters https://doi.org/10.1038/s41550-017-0341-8 In the format provided by the authors and unedited. Characterization of methanol as a magnetic field tracer in star-forming regions
More informationInvestigating Popular Water Models
Investigating Popular Water odels Albert DeFusco, Ken Jordan Center for olecular and aterials Simulations What are we modeling? OOOOOLLLLLDDDDDEDEEEENENNNNNN OOLL D E OOOOOLLLLLDDDDDEDEEEENENNNNNN OOLL
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 information1 Density functional theory (DFT)
1 Density functional theory (DFT) 1.1 Introduction Density functional theory is an alternative to ab initio methods for solving the nonrelativistic, time-independent Schrödinger equation H Φ = E Φ. The
More informationSolution of the Electronic Schrödinger Equation. Using Basis Sets to Solve the Electronic Schrödinger Equation with Electron Correlation
Solution of the Electronic Schrödinger Equation Using Basis Sets to Solve the Electronic Schrödinger Equation with Electron Correlation Errors in HF Predictions: Binding Energies D e (kcal/mol) HF Expt
More informationExercise 1: Structure and dipole moment of a small molecule
Introduction to computational chemistry Exercise 1: Structure and dipole moment of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the dipole moment of a small
More informationNoncovalent interactions in biochemistry
Noncovalent interactions in biochemistry Kevin E. Riley 1 and Pavel Hobza 2,3 Noncovalent interactions are known to play a key role in biochemistry. The knowledge of stabilization (relative) energies and
More informationSolutions and Non-Covalent Binding Forces
Chapter 3 Solutions and Non-Covalent Binding Forces 3.1 Solvent and solution properties Molecules stick together using the following forces: dipole-dipole, dipole-induced dipole, hydrogen bond, van der
More informationTerahertz vibration-rotation-tunneling spectroscopy of the water tetramer-d8: Combined analysis of vibrational bands at 4.1 and 2.
THE JOURNAL OF CHEMICAL PHYSICS 128, 094302 2008 Terahertz vibration-rotation-tunneling spectroscopy of the water tetramer-d8: Combined analysis of vibrational bands at 4.1 and 2.0 THz Wei Lin, a Jia-Xiang
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 informationChemistry 4560/5560 Molecular Modeling Fall 2014
Final Exam Name:. User s guide: 1. Read questions carefully and make sure you understand them before answering (if not, ask). 2. Answer only the question that is asked, not a different question. 3. Unless
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 informationIntroduction to DFTB. Marcus Elstner. July 28, 2006
Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD
More information3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D.
3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D. Thierry Langer, Ph.D. Jana Vrbková, Ph.D. UP Olomouc, 23.1.-26.1. 2018
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationSupplementary information
Matthias Heger, Tina Scharge, and Martin A. Suhm Institute of Physical Chemistry, Georg-August-Universität, Tammannstraße 6, 37077 Göttingen, Germany. E-mail: msuhm@gwdg.de Current address: Gesellschaft
More informationStructure of diatomic molecules
Structure of diatomic molecules January 8, 00 1 Nature of molecules; energies of molecular motions Molecules are of course atoms that are held together by shared valence electrons. That is, most of each
More informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Hylleraas Centre, Department of Chemistry, University of Oslo, Norway and Centre for Advanced Study at the Norwegian Academy of Science and Letters, Oslo,
More informationIFM Chemistry Computational Chemistry 2010, 7.5 hp LAB2. Computer laboratory exercise 1 (LAB2): Quantum chemical calculations
Computer laboratory exercise 1 (LAB2): Quantum chemical calculations Introduction: The objective of the second computer laboratory exercise is to get acquainted with a program for performing quantum chemical
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 informationThe Self Interaction Correction revisited
The Self Interaction Correction revisited Explicit dynamics of clusters and molecules under irradiation Spectroscopic accuracy at low energy SIC problem : one electron interacts with its own mean-field!
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 informationValence electronic structure of isopropyl iodide investigated by electron momentum spectroscopy. --- Influence of intramolecular interactions
Valence electronic structure of isopropyl iodide investigated by electron momentum spectroscopy --- Influence of intramolecular interactions Minfu Zhao, Xu Shan, Shanshan Niu, Yaguo Tang, Zhaohui Liu,
More informationElectronic structure theory: Fundamentals to frontiers. VI. Analysis and more.
Electronic structure theory: Fundamentals to frontiers. VI. Analysis and more. MARTIN HEAD-GORDON Department of Chemistry, University of California, Berkeley, and, Chemical Sciences Division, Lawrence
More informationCharge-Transfer and Dispersion Energies in Water Clusters
II.26 Charge-Transfer and Dispersion Energies in Water Clusters Suehiro Iwata 1,2, Pradipta Bandyopadhyay 3, Sotiris S. Xantheas 4 1)Toyota Physical and Chemical Research Institute (2008-2012, fellow)
More informationLow-temperature isomers of the water hexamer
Low-temperature isomers of the water hexamer Volodymyr Babin and Francesco Paesani Department of Chemistry & Biochemistry, University of California, San Diego, La Jolla, CA 92093-0314 Hundreds of years
More informationMolecular Physics. Attraction between the ions causes the chemical bond.
Molecular Physics A molecule is a stable configuration of electron(s) and more than one nucleus. Two types of bonds: covalent and ionic (two extremes of same process) Covalent Bond Electron is in a molecular
More information