Computational chemistry with GAMESS: a very brief overview with examples

Size: px
Start display at page:

Download "Computational chemistry with GAMESS: a very brief overview with examples"

Transcription

1 Computational chemistry with GAMESS: a very brief overview with examples PHY-6120 Molecular Physics (Spring 2015), UConn Phys. Dept. Feb 17 th 2015

2 H = ħ2 2μ i Intro: V(R) for diatomic molecules + k Z A Z B e2 R AB i 2 ħ2 2 I I,i 1 2 M I + I i> j k Z I e 2 Born-Oppenheimer approximation: separate electronic and nuclear motion k e 2 fix R = R AB, solve for the electronic motion, then solve for the slower nuclear motion higher-order electronic correlations (core-valence) neglected Energies and total wavefunctions r ij r A, i R AB r Ii +... A B M I m e i E = E elec + G (v)+ F ( J ) ψ i, v BO elec. = Φ i, Λ, S, Σ (r, R)χ nucl. v ( R) Non-trivial (14 terms total, different bases sets) plus non-born-oppenheimer terms, perturbations,...

3 LiNa molecule, 2 lowest electronic states V ex R = ±A R exp R V lr R = C 6 R C 8 6 R C 10 8 R 10

4 Quantum chemistry ab-initio codes Popular: GAMESS (free, opensource) NWCHEM (mostly commercial) GAUSSIAN (commercial) MOLPRO (commercial) Common methods: Hartree-Fock (HF) post-hartree-fock CC, CCSD, CCSD(T) MP2 CI and variants Density Functional Theory (DFT)

5 GAMESS on out machines our machines: hex.phys.uconn.edu (16 CPUs, 32 GB RAM) octane.phys.uconn.edu (8 CPUs, 8 GB RAM) pauli.phys.uconn.edu, fermi.phys.uconn.edu (Côté group research) on hex.phys.uconn.edu: /scratch/gacesa/gamess -> GAMESS installation /scratch/gacesa/gamess-class/_sample-input.zip -> sample inputs create folder corresponding to your username in /scratch/gacesa/gamess-class/ : e.g. mkdir /scratch/gacesa/gamess-class/gacesa

6 GAMESS input Main input groups: $CONTRL - (type of wavefunction, method, things to calculate) $SYSTEM - system parameters (memory, directories,...) $BASIS - the basis is specified here $DATA (required) - molecular composition and geometry selected optional groups $ZMAT (if NZVAR>0 in $CNTRL), $SCF (if $SCFTYP = RHF, UHF, or ROHF; required for GVB) $DFT, $TDDFT,..., $CIS (excited states), $MP2, $CP(MC)HF (solution of response HF equations), $MASS (isotopes),...

7 $CONTRL group SCFTYP = { RHF, UHF, ROHF, GVB, MCSCF, none } elec. correlation: MPLEVL, CITYP, CCTYP, DFTTYP, VBTYP,... RUNTYP = { ENERGY (def), GRADIENT, HESSIAN, OPTIMIZE, SADPOINT, OPTFMO, PROP, NACME, TRANSITN,... } ICHARG = molecular charge (default = 0) MULT = 1 (singlet = default), 2 (doublet), 3 (triplet),... Geometry: COORD (mol. coord.), UNITS, NZVAR Pseudopotentials: PP = { NONE, READ ($ECP), MCP,... } LOCAL = {NONE, BOYS, SVD,... } - orbital localization ($LOCAL)

8 $BASIS group Three ways to specify a basis (GBASIS, BASNAM, EXTFIL); or in $DATA GBASIS = MINI, DZV, TZV (good for start; ) ECP, MCP, empirical... check page of input.pdf custom: support.scinet.utoronto.ca/wiki/index.php/user_codes#tips_from_the_fekl_lab

9 $DATA group TITLE GROUP, NAXIS specified as a symmetry group (en.wikipedia.org/wiki/schoenflies_notation) e.g. Cnv 2 = C2v ; note: C1 = no symmetry; see e.g. Cotton "Group Theory" skip Lines 3 and 4 (for our purposes) and replace with a blank line for COORD=UNIQUE in ($CONTRL): (other formats possible) NAME (10 chars), ZNUC (nucl. charge), X, Y, Z (coordinates) Example for H-H molecule: $DATA hydrogen molecule, no symmetry C1 H H $END

10 A very simple input file for H 2 (X 1 Σ + ) H2-rhf.inp $CONTRL SCFTYP=RHF RUNTYP=ENERGY $END $BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 $END $DATA hydrogen molecule with 6-G31 C1 H H $END Pople notation for the basis: 6 Gaussians for the inner elec, 4 (3+1) for valence Polarized orbitals included in the basis set demonstration: run GAMESS for H2-rhf.inp (on hex):./rungms H2-rhf.inp >& H2-rhf.log

11 GAMESS output /scratch/gamess-class/$user/h2-rhf.log check for these: number of SCF cycles needed for convergence (section "RHF SCF CALCULATION") the HF energy (section "ENERGY COMPONENTS") the dipole moment (section "ELECTROSTATIC MOMENTS") (optional) visualization (using VMD or similar)

12 Example: LiH potential energy curve (PES) click on constants of diatomic molecules r e = Bohr (ground state) demonstration #1: run GAMESS for LiH-rhf.inp edit/create the file -> replace '???' with correct values./rungms lih-rhf.inp >& lih-rhf.log demonstration #2: lih-rhf-scan.inp edit/create the file lih-rhf-scan.inp./rungms lih-rhf-scan.inp >& lih-rhf-scan.log./scanlog2dat.sh./xmgrace -free lih-rhf-scan_pes.dat

13 Excited states simplest: Equation-of-motion (EOM) methods works with RHF, CCSD(T), and some higher level methods pay attention to the symmetry / electron multiplicity! parsing output becomes more complicated add CCTYP into $CONTRL, and $CCINP & $EOMINP groups $CONTRL SCFTYP=RHF CCTYP=EOM-CCSD RUNTYP=ENERGY $END $BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 $END $CCINP MULT=1 $END $DATA H2 with 6-G31 + first excited state C1 H H $END

14 Details about basis sets

15

AN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University

AN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,

More information

Chapter: 22. Visualization: Making INPUT File and Processing of Output Results

Chapter: 22. Visualization: Making INPUT File and Processing of Output Results Chapter: 22 Visualization: Making INPUT File and Processing of Output Results Keywords: visualization, input and output structure, molecular orbital, electron density. In the previous chapters, we have

More information

Chemistry 4560/5560 Molecular Modeling Fall 2014

Chemistry 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 information

Ab initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier

Ab 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 information

AN INTRODUCTION TO MCSCF: PART 2

AN INTRODUCTION TO MCSCF: PART 2 AN INTRODUCTION TO MCSCF: PART 2 ORBITAL APPROXIMATION Hartree product (hp) expressed as a product of spinorbitals ψ ι = φ i σ i φ i = space orbital, σ i = spin function (α,β) Pauli Principle requires

More information

DFT calculations of NMR indirect spin spin coupling constants

DFT 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 information

MO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )

MO 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 information

NWChem: Coupled Cluster Method (Tensor Contraction Engine)

NWChem: Coupled Cluster Method (Tensor Contraction Engine) NWChem: Coupled Cluster Method (Tensor Contraction Engine) Why CC is important?! Correlation effects are important!! CC is size-extensive theory: can be used to describe dissociation processes.! Higher-order

More information

Session 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education

Session 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 information

Chemistry 334 Part 2: Computational Quantum Chemistry

Chemistry 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 information

Electric properties of molecules

Electric 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 information

Using Web-Based Computations in Organic Chemistry

Using Web-Based Computations in Organic Chemistry 10/30/2017 1 Using Web-Based Computations in Organic Chemistry John Keller UAF Department of Chemistry & Biochemistry The UAF WebMO site Practical aspects of computational chemistry theory and nomenclature

More information

Electron Correlation

Electron Correlation Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation

More information

EXECUTION OF GAMESS BEGUN 16:18:

EXECUTION OF GAMESS BEGUN 16:18: 1 ****************************************************** * GAMESS VERSION = 6 MAY 1998 * * FROM IOWA STATE UNIVERSITY * * M.W.SCHMIDT, K.K.BALDRIDGE, J.A.BOATZ, S.T.ELBERT, * * M.S.GORDON, J.H.JENSEN,

More information

IFM Chemistry Computational Chemistry 2010, 7.5 hp LAB2. Computer laboratory exercise 1 (LAB2): Quantum chemical calculations

IFM 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 information

(23 Sep 99) Input to PC GAMESS may be in upper or lower case. There are three types of input groups in PC GAMESS:

(23 Sep 99) Input to PC GAMESS may be in upper or lower case. There are three types of input groups in PC GAMESS: ************************************ * * * CHAPTER I - Introduction * * * ************************************ (23 Sep 99) Input Philosophy Input to PC GAMESS may be in upper or lower case. There are three

More information

Input Description 2-1. (21 May 2009) ********************************* * * * Section 2 - Input Description * * * *********************************

Input Description 2-1. (21 May 2009) ********************************* * * * Section 2 - Input Description * * * ********************************* Input Description 2-1 ********************************* * * * Section 2 - Input Description * * * ********************************* (21 May 2009) This section of the manual describes the input to GAMESS.

More information

An Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine

An 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 information

Input Description 2-1. (21 April 2008) ********************************* * * * Section 2 - Input Description * * * *********************************

Input Description 2-1. (21 April 2008) ********************************* * * * Section 2 - Input Description * * * ********************************* Input Description 2-1 ********************************* * * * Section 2 - Input Description * * * ********************************* (21 April 2008) This section of the manual describes the input to GAMESS.

More information

NWChem: Coupled Cluster Method (Tensor Contraction Engine)

NWChem: Coupled Cluster Method (Tensor Contraction Engine) NWChem: Coupled Cluster Method (ensor Contraction Engine) What we want to solve H Ψ = E Ψ Many Particle Systems Molecular/Atomic Physics, Quantum Chemistry (electronic Schrödinger equations) Solid State

More information

Beyond the Hartree-Fock Approximation: Configuration Interaction

Beyond the Hartree-Fock Approximation: Configuration Interaction Beyond the Hartree-Fock Approximation: Configuration Interaction The Hartree-Fock (HF) method uses a single determinant (single electronic configuration) description of the electronic wavefunction. For

More information

This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry.

This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry. 1 Computational Chemistry (Quantum Chemistry) Primer This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry. TABLE OF CONTENTS Methods...1 Basis

More information

Section 2 - Input Description

Section 2 - Input Description 2-1 ( 7 Aug 02) Section 2 - Input Description This section of the manual describes the input to GAMESS. The section is written in a reference, rather than tutorial fashion. However, there are frequent

More information

Fragmentation methods

Fragmentation 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 information

The successful wavefunction can be written as a determinant: # 1 (2) # 2 (2) Electrons. This can be generalized to our 2N-electron wavefunction:

The successful wavefunction can be written as a determinant: # 1 (2) # 2 (2) Electrons. This can be generalized to our 2N-electron wavefunction: T2. CNDO to AM1: The Semiempirical Molecular Orbital Models The discussion in sections T2.1 T2.3 applies also to ab initio molecular orbital calculations. T2.1 Slater Determinants Consider the general

More information

Computational Material Science Part II. Ito Chao ( ) Institute of Chemistry Academia Sinica

Computational Material Science Part II. Ito Chao ( ) Institute of Chemistry Academia Sinica Computational Material Science Part II Ito Chao ( ) Institute of Chemistry Academia Sinica Ab Initio Implementations of Hartree-Fock Molecular Orbital Theory Fundamental assumption of HF theory: each electron

More information

( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r)

( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r) Born Oppenheimer Approximation: Ĥ el ( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r) For a molecule with N electrons and M nuclei: Ĥ el What is E el (R)? s* potential surface Reaction Barrier Unstable intermediate

More information

Handbook of Computational Quantum Chemistry

Handbook of Computational Quantum Chemistry Handbook of Computational Quantum Chemistry David B. Cook Dept. of Chemistry University of Sheffield DOVER PUBLICATIONS, INC. Mineola, New York F Contents 1 Mechanics and molecules 1 1.1 1.2 1.3 1.4 1.5

More information

Molecular Simulation I

Molecular Simulation I Molecular Simulation I Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences

More information

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Jack 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 information

Introduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule. Vesa Hänninen

Introduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule. Vesa Hänninen Introduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the electronic energy

More information

计算物理作业二. Excercise 1: Illustration of the convergence of the dissociation energy for H 2 toward HF limit.

计算物理作业二. Excercise 1: Illustration of the convergence of the dissociation energy for H 2 toward HF limit. 计算物理作业二 Excercise 1: Illustration of the convergence of the dissociation energy for H 2 toward HF limit. In this exercise, basis indicates one of the following basis sets: STO-3G, cc-pvdz, cc-pvtz, cc-pvqz

More information

Exercise 1: Structure and dipole moment of a small molecule

Exercise 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 information

Oslo node. Highly accurate calculations benchmarking and extrapolations

Oslo 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 information

MRCI calculations in MOLPRO

MRCI calculations in MOLPRO 1 MRCI calculations in MOLPRO Molpro is a software package written in Fortran and maintained by H.J. Werner and P.J. Knowles. It is often used for performing sophisticated electronic structure calculations,

More information

GEM4 Summer School OpenCourseWare

GEM4 Summer School OpenCourseWare GEM4 Summer School OpenCourseWare http://gem4.educommons.net/ http://www.gem4.org/ Lecture: Molecular Mechanics by Ju Li. Given August 9, 2006 during the GEM4 session at MIT in Cambridge, MA. Please use

More information

Multi-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005

Multi-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson

More information

Molecular Mechanics: The Ab Initio Foundation

Molecular Mechanics: The Ab Initio Foundation Molecular Mechanics: The Ab Initio Foundation Ju Li GEM4 Summer School 2006 Cell and Molecular Mechanics in BioMedicine August 7 18, 2006, MIT, Cambridge, MA, USA 2 Outline Why are electrons quantum? Born-Oppenheimer

More information

Jack Smith. Center for Environmental, Geotechnical and Applied Science. Marshall University

Jack Smith. Center for Environmental, Geotechnical and Applied Science. Marshall University Jack Smith Center for Environmental, Geotechnical and Applied Science Marshall University -- Division of Science and Research WV Higher Education Policy Commission WVU HPC Summer Institute June 20, 2014

More information

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby

Advanced 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 information

Handbook of Computational Quantum Chemistry. DAVID B. COOK The Department of Chemistry, University of Sheffield

Handbook of Computational Quantum Chemistry. DAVID B. COOK The Department of Chemistry, University of Sheffield Handbook of Computational Quantum Chemistry DAVID B. COOK The Department of Chemistry, University of Sheffield Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1998 CONTENTS 1 Mechanics and molecules 1 1.1

More information

NMR and IR spectra & vibrational analysis

NMR and IR spectra & vibrational analysis Lab 5: NMR and IR spectra & vibrational analysis A brief theoretical background 1 Some of the available chemical quantum methods for calculating NMR chemical shifts are based on the Hartree-Fock self-consistent

More information

Introduction to QM simulations

Introduction to QM simulations QM_tutorial Page 1 of 2 Next: Introduction University of Illinois at Urbana-Champaign Beckman Institute for Advanced Science and Technology Theoretical and Computational Biophysics Group Introduction to

More information

QUANTUM CHEMISTRY FOR TRANSITION METALS

QUANTUM 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 information

Computational Chemistry Using the University of Alaska WebMO Site

Computational Chemistry Using the University of Alaska WebMO Site 2/7/2017 1 Computational Chemistry Using the University of Alaska WebMO Site John Keller Department of Chemistry & Biochemistry University of Alaska Fairbanks Intro and operation of WebMO and MOPAC Basic

More information

Chemistry 3502/4502. Final Exam Part I. May 14, 2005

Chemistry 3502/4502. Final Exam Part I. May 14, 2005 Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle

More information

Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012

Joint 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 information

CHE3935. Lecture 4 Quantum Mechanical Simulation Methods Continued

CHE3935. Lecture 4 Quantum Mechanical Simulation Methods Continued CHE3935 Lecture 4 Quantum Mechanical Simulation Methods Continued 1 OUTLINE Review Introduction to CPMD MD and ensembles The functionals of density functional theory Return to ab initio methods Binding

More information

CHEM 545 Theory and Practice of Molecular Electronic Structure. Anna I. Krylov. DON T PANIC.

CHEM 545 Theory and Practice of Molecular Electronic Structure. Anna I. Krylov.   DON T PANIC. CHEM 545 Theory and Practice of Molecular Electronic Structure Anna I. Krylov http://iopenshell.usc.edu/chem545/ DON T PANIC USC Fall 2014 Things to do: 1. Install IQmol (by this Thursday). http://iqmol.org/.

More information

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Jack 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 information

This is called a singlet or spin singlet, because the so called multiplicity, or number of possible orientations of the total spin, which is

This is called a singlet or spin singlet, because the so called multiplicity, or number of possible orientations of the total spin, which is 9. Open shell systems The derivation of Hartree-Fock equations (Chapter 7) was done for a special case of a closed shell systems. Closed shell means that each MO is occupied by two electrons with the opposite

More information

Computational Chemistry I

Computational Chemistry I Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of

More information

Computational Chemistry. An Introduction to Molecular Dynamic Simulations

Computational 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 information

Computational Methods. Chem 561

Computational 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 information

Electronic structure calculations: fundamentals George C. Schatz Northwestern University

Electronic structure calculations: fundamentals George C. Schatz Northwestern University Electronic structure calculations: fundamentals George C. Schatz Northwestern University Electronic Structure (often called Quantum Chemistry) calculations use quantum mechanics to determine the wavefunctions

More information

QUANTUM CHEMISTRY PROJECT 3: PARTS B AND C

QUANTUM 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 information

NWChem: Hartree-Fock, Density Functional Theory, Time-Dependent Density Functional Theory

NWChem: Hartree-Fock, Density Functional Theory, Time-Dependent Density Functional Theory NWChem: Hartree-Fock, Density Functional Theory, Time-Depent Density Functional Theory Hartree-Fock! Functionality! Input! Wavefunctions! Initial MO vectors! Direct and semidirect algorithms! Convergence,

More information

Exchange 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 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 information

Basic introduction of NWChem software

Basic introduction of NWChem software Basic introduction of NWChem software Background NWChem is part of the Molecular Science Software Suite Designed and developed to be a highly efficient and portable Massively Parallel computational chemistry

More information

CASSCF and NEVPT2 calculations: Ground and excited states of multireference systems. A case study of Ni(CO)4 and the magnetic system NArO

CASSCF and NEVPT2 calculations: Ground and excited states of multireference systems. A case study of Ni(CO)4 and the magnetic system NArO CASSCF and NEVPT2 calculations: Ground and excited states of multireference systems. A case study of Ni(CO)4 and the magnetic system NArO The ground states of many molecules are often well described by

More information

Performance of Hartree-Fock and Correlated Methods

Performance of Hartree-Fock and Correlated Methods Chemistry 460 Fall 2017 Dr. Jean M. Standard December 4, 2017 Performance of Hartree-Fock and Correlated Methods Hartree-Fock Methods Hartree-Fock methods generally yield optimized geomtries and molecular

More information

Technical Note Calculations of Orbital Overlap Range Function EDR( r ; d) and Overlap Distance D(r )using Multiwfn

Technical Note Calculations of Orbital Overlap Range Function EDR( r ; d) and Overlap Distance D(r )using Multiwfn Technical Note Calculations of Orbital Overlap Range Function EDR( r ; d) and Overlap Distance D(r )using Multiwfn Abstract The orbital overlap range function EDR( r; d) (J. Chem. Phys. 2014, 141, 144104)

More information

Introduction to Computational Quantum Chemistry: Theory

Introduction to Computational Quantum Chemistry: Theory Introduction to Computational Quantum Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC 3108 Course Lectures 2007 Introduction Hartree Fock Theory Configuration Interaction Lectures 1 Introduction

More information

Quantum Chemical Simulations and Descriptors. Dr. Antonio Chana, Dr. Mosè Casalegno

Quantum 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 information

Density Func,onal Theory (Chapter 6, Jensen)

Density 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 information

Introduction to Hartree-Fock calculations in Spartan

Introduction to Hartree-Fock calculations in Spartan EE5 in 2008 Hannes Jónsson Introduction to Hartree-Fock calculations in Spartan In this exercise, you will get to use state of the art software for carrying out calculations of wavefunctions for molecues,

More information

Vol. 9 COMPUTATIONAL CHEMISTRY 319

Vol. 9 COMPUTATIONAL CHEMISTRY 319 Vol. 9 COMPUTATIONAL CHEMISTRY 319 COMPUTATIONAL QUANTUM CHEMISTRY FOR FREE-RADICAL POLYMERIZATION Introduction Chemistry is traditionally thought of as an experimental science, but recent rapid and continuing

More information

Gaussian Basis Sets for Solid-State Calculations

Gaussian Basis Sets for Solid-State Calculations Gaussian Basis Sets for Solid-State Calculations K. Doll Molpro Quantum Chemistry Software Institute of Theoretical Chemistry, D-70569 Stuttgart, Germany MW-MSSC 2017, Minneapolis, July 10, 2017 Introduction

More information

Lec20 Fri 3mar17

Lec20 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 information

On the Development of a New Computational Chemistry Software

On the Development of a New Computational Chemistry Software On the Development of a New Computational Chemistry Software Han Ung Lee, Hayan Lee and Wilfredo Credo Chung* Department of Chemistry, De La Salle University Manila, 2401 Taft Avenue, Manila, 1004 Philippines

More information

2~:J~ -ryej- r- 2 Jr. A - f3. sr(djk nv~tor rn~ +~ rvjs (::-CJ) ::;-1-.'--~ -. rhd. ('-.Ji.L.~ )- r'-d)c, -r/~ JJr - 2~d ~2-Jr fn'6.

2~:J~ -ryej- r- 2 Jr. A - f3. sr(djk nv~tor rn~ +~ rvjs (::-CJ) ::;-1-.'--~ -. rhd. ('-.Ji.L.~ )- r'-d)c, -r/~ JJr - 2~d ~2-Jr fn'6. .~, ~ I, sr(djk nv~tor rn~ +~ rvjs (::-CJ) ::;-1-.'--~ -. rhd. ('-.Ji.L.~ )- r'-d)c, -r/~ JJr - 2~d ~2-Jr fn'6.)1e'" 21t-ol Je C'...-------- lj-vi, J? Jr Jr \Ji 2~:J~ -ryej- r- 2 Jr A - f3 c _,~,= ~,.,w._..._.

More information

4 Post-Hartree Fock Methods: MPn and Configuration Interaction

4 Post-Hartree Fock Methods: MPn and Configuration Interaction 4 Post-Hartree Fock Methods: MPn and Configuration Interaction In the limit of a complete basis, the Hartree-Fock (HF) energy in the complete basis set limit (ECBS HF ) yields an upper boundary to the

More information

use the backs of pages as needed

use the backs of pages as needed CHEMISTRY 4021/8021 Q1) Propose a simple, united-atom molecular mechanics force-field needed to generate a potential energy surface for an isolated molecule of acetone (Me 2 CO). I.e., provide an energy

More information

Electron Correlation Methods

Electron Correlation Methods Electron Correlation Methods HF method: electron-electron interaction is replaced by an average interaction E HF c = E 0 E HF E 0 exact ground state energy E HF HF energy for a given basis set HF E c

More information

Introduction to Computational Chemistry: Theory

Introduction to Computational Chemistry: Theory Introduction to Computational Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC andrew.gilbert@anu.edu.au 3023 Course Lectures Introduction Hartree Fock Theory Basis Sets Lecture 1 1 Introduction

More information

Chemistry 3502/4502. Final Exam Part I. May 14, 2005

Chemistry 3502/4502. Final Exam Part I. May 14, 2005 Chemistry 3502/4502 Final Exam Part I May 14, 2005 1. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle (e) The

More information

Appendix C Calculating Excited States using Gaussian

Appendix C Calculating Excited States using Gaussian Appendix C Calculating Excited States using Gaussian 403 This appendix contains methods for using Gaussian 03 78 and Gaussian 09 121 to calculate excited states of molecules. Such methods are useful for

More information

Exam 4 Review. Exam Review: A exam review sheet for exam 4 will be posted on the course webpage. Additionally, a practice exam will also be posted.

Exam 4 Review. Exam Review: A exam review sheet for exam 4 will be posted on the course webpage. Additionally, a practice exam will also be posted. Chem 4502 Quantum Mechanics & Spectroscopy (Jason Goodpaster) Exam 4 Review Exam Review: A exam review sheet for exam 4 will be posted on the course webpage. Additionally, a practice exam will also be

More information

OVERVIEW OF QUANTUM CHEMISTRY METHODS

OVERVIEW 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 information

Molecular Orbitals for Ozone

Molecular Orbitals for Ozone Molecular Orbitals for Ozone Purpose: In this exercise you will do semi-empirical molecular orbital calculations on ozone with the goal of understanding the molecular orbital print out provided by Spartan

More information

Lecture 4: methods and terminology, part II

Lecture 4: methods and terminology, part II So theory guys have got it made in rooms free of pollution. Instead of problems with the reflux, they have only solutions... In other words, experimentalists will likely die of cancer From working hard,

More information

Introduction to Hartree-Fock calculations using ORCA and Chemcraft

Introduction to Hartree-Fock calculations using ORCA and Chemcraft Introduction to Hartree-Fock calculations using ORCA and Chemcraft In this exercise, you will get to use software for carrying out calculations of wavefunctions for molecules, the ORCA program. While ORCA

More information

Electron Correlation - Methods beyond Hartree-Fock

Electron 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 information

Study of Ozone in Tribhuvan University, Kathmandu, Nepal. Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal

Study of Ozone in Tribhuvan University, Kathmandu, Nepal. Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal Study of Ozone in Tribhuvan University, Kathmandu, Nepal Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal 1 Country of the Mt Everest 2 View of the Mt Everest 3 4 5

More information

INPUT DESCRIPTION FOR SQM version 1.0

INPUT DESCRIPTION FOR SQM version 1.0 INPUT DESCRIPTION FOR SQM version 1.0 INTRODUCTION SQM is an add-on module for the PQS program which scales force constants to produce a Scaled Quantum Mechanical (SQM) Force Field. This can correct for

More information

Density Functional Theory

Density Functional Theory Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function

More information

Chemistry, Physics and the Born- Oppenheimer Approximation

Chemistry, Physics and the Born- Oppenheimer Approximation Chemistry, Physics and the Born- Oppenheimer Approximation Hendrik J. Monkhorst Quantum Theory Project University of Florida Gainesville, FL 32611-8435 Outline 1. Structure of Matter 2. Shell Models 3.

More information

Walter 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 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 information

Q-Chem Workshop Examples Part 2

Q-Chem Workshop Examples Part 2 Marek Freindorf Q-Chem Workshop Examples Part 2 Louisville, KY March 2010 CO, MP2 Calculations Open Avogadro using the "Build" button of QUI Page 2 CO, MP2 Calculations Using Avogadro to create coordinates

More information

Input parameters for B3LYP/6-31G* stationary point searches in GAMESS. Coordinates (in bohr) and energies (in a.u.) of structure 1-E A

Input parameters for B3LYP/6-31G* stationary point searches in GAMESS. Coordinates (in bohr) and energies (in a.u.) of structure 1-E A Conformational Analysis via Calculations and NMR Spectroscopy for Isomers of the Mono(imino)pyridine Ligand, 2-{(2,6-Me 2 - C 6 H 3 )NC(i-Pr)}C 5 H 4 N Timothy J. Dudley*, Jennifer E. Beck, Earl E. P.

More information

QUANTUM CHEMISTRY PROJECT 3: ATOMIC AND MOLECULAR STRUCTURE

QUANTUM CHEMISTRY PROJECT 3: ATOMIC AND MOLECULAR STRUCTURE Chemistry 460 Fall 2017 Dr. Jean M. Standard November 1, 2017 QUANTUM CHEMISTRY PROJECT 3: ATOMIC AND MOLECULAR STRUCTURE OUTLINE In this project, you will carry out quantum mechanical calculations of

More information

EXAM INFORMATION. Radial Distribution Function: B is the normalization constant. d dx. p 2 Operator: Heisenberg Uncertainty Principle:

EXAM 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 information

INPUT DESCRIPTION FOR SQM version 2.0

INPUT DESCRIPTION FOR SQM version 2.0 INPUT DESCRIPTION FOR SQM version 2.0 INTRODUCTION SQM is an add-on module for the PQS program which scales force constants to produce a Scaled Quantum Mechanical (SQM) Force Field. This can correct for

More information

The TopMoD. user s manual. Bernard Silvi Laboratoire de Chime Théorique Université Pierre et Marie Curie

The TopMoD. user s manual. Bernard Silvi Laboratoire de Chime Théorique Université Pierre et Marie Curie The TopMoD user s manual Bernard Silvi Laboratoire de Chime Théorique Université Pierre et Marie Curie Contents 1 General Presentation 5 1.1 Useful reading and references................................

More information

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Jack 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 information

Pseudo-Hermitian eigenvalue equations in linear-response electronic-structure theory

Pseudo-Hermitian eigenvalue equations in linear-response electronic-structure theory 1/11 Pseudo-Hermitian eigenvalue equations in linear-response electronic-structure theory Julien Toulouse Université Pierre & Marie Curie and CNRS, 4 place Jussieu, Paris, France Web page: www.lct.jussieu.fr/pagesperso/toulouse/

More information

Introduction to Hartree-Fock Molecular Orbital Theory

Introduction to Hartree-Fock Molecular Orbital Theory Introduction to Hartree-Fock Molecular Orbital Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Origins of Mathematical Modeling in Chemistry Plato (ca. 428-347

More information

Coupled-Cluster Perturbative Triples for Bond Breaking

Coupled-Cluster Perturbative Triples for Bond Breaking Coupled-Cluster Perturbative Triples for Bond Breaking Andrew G. Taube and Rodney J. Bartlett Quantum Theory Project University of Florida INT CC Meeting Seattle July 8, 2008 Why does chemistry need triples?

More information

DFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory

DFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory DFT and beyond: Hands-on Tutorial Workshop 2011 Tutorial 1: Basics of Electronic Structure Theory V. Atalla, O. T. Hofmann, S. V. Levchenko Theory Department, Fritz-Haber-Institut der MPG Berlin July 13,

More information

Computational Chemistry and Quantum Tunneling

Computational Chemistry and Quantum Tunneling Computational Chemistry and Quantum Tunneling Key Words quantum chemistry, Hartree-Fock theory, density functional theory, variational principle, quantum tunneling, transition state Object It is expected

More information