Computational Chemistry Using the University of Alaska WebMO Site

Size: px
Start display at page:

Download "Computational Chemistry Using the University of Alaska WebMO Site"

Transcription

1 2/7/ 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 Gaussian

2 2/7/ Click here to log in (guest/webmo or yourusername/password)

3 2/7/ How does a WebMO job work? webserver world wide web User Draw molecule & choose options WebMO creates input file for MOPAC or Gaussian When the job is done, the output file is parsed and shown on Results page. The updated, intermediate results can be visualized in the Job Manager. compute servers chemlinux1 chemlinux2 (& others). Running Gaussian, MOPAC, and NWChem

4 2/7/ User Accounts MOPAC NWChem Gaussian CPU time limit (h) guest 4 Non-UAF unlimited UAF unlimited Machines Server Cores Procs Mem (GB) CPU Gaussian NWChem MOPAC Chemlinux Xeon (2) (6.3) Chemlinux Xeon (2) Chemlinux Xeon (2) (6.3) Antec Core i7 5820K (6.6) Corsair Core i7 6950X (6.6) Corsair Core i7 930 (WebMO)

5 2/7/ Levels of theory Correlated ab initio MO (MP2 & others) Molecular mechanics (MM+, OPLS, etc) Ball-andspring model no orbitals no electronics Semi-empirical MO (PM3, PM6) Uses 27 orbital parameters for each element Parameters derived from stable molecules Electron correlation implicit, since the parameters are fr real molecules, where actual electron correlation occurs. Ab initio MO (Hartree-Fock HF) MOs = linear combination of basis orbitals. Parameters derived from isolated atoms Some electron correlation Density Functional Theory (DFT) (B3LYP ) MOs = linear combination of basis orbitals Parameters derived from atoms and molecular databases More electron correlation MOs = linear combination of basis orbitals. Parameters derived from isolated atoms Includes excited state configurations Most electron correlation

6 2/7/ Basic WebMO and MOPAC 1) Build Example : ClOOF or similar, starting with trans conformation. 2) Optimize Use MOPAC PM3 or PM6. 3) Symmetry 4) Vibrations 5) Adjust and re-optimize. 6) MO surfaces

7 2/7/ Log in

8 2/7/ Begin with New Job, Create New Job

9 2/7/ Choose atom type ( F key for fluorine, O key for oxygens.. Use the Build tool to add atoms and bonds. Then do Cleanup, Comprehensive- Idealized Cleanup adjusts angles and bond lengths to their nominal values

10 2/7/ Optimize your molecule 1) Right Arrow 2) Engine: MOPAC 3) Select Server: First Available

11 2/7/ Enter an informative Job Name Type of calculation Geometry Optimization Theory: PM6 (a recent version of PM3 semi-empirical) Charge: (0 for molecule) Multiplicity: generally singlet (unless odd # of valence e - )

12 2/7/

13 2/7/ Generate on the Preview tab shows the input file automatically created by WebMO. It can be modified manually. Submit job by clicking the right arrow

14 2/7/ When the job is Complete, click the job name or icon to see the Results View Job: 1) The symmetry elements in your optimized molecule (Show Symmetry Elements button 4 th from top left) 2) Heat of Formation 3) Geometry Sequence Energies plot Now click New Job Using This Geometry 1) Job Name: JK-FOOF-PM6-Freq 2) Calculation: Vibrational Frequencies 3) Theory: PM6 4) Charge: 0 5) Multiplicity: Singlet Submit job

15 2/7/ ! Click here to animate the vibrational mode. The negative frequency is just a symbol for an imaginary frequency.

16 Energy 2/7/ Surprisingly often, WebMO and HyperChem build molecules with imaginary vibrations. These are transition states, not stable states. Common in organic molecules with rotatable bonds. When these very symmetric structures are optimized, they sometimes remain symmetric! Why? The geometry optimization algorithm tests the E vs coordinate slope, which = 0 at both maxima and minima. How do you know if this has happened? A vibration frequency calculation will show a negative (imaginary) frequency. The cure? Manually nudge the system off the maximum, and restart the optimization. Geometry coordinate

17 2/7/ Basic WebMO and MOPAC 1) Build 2) Optimize 3) Symmetry 4) Vibrations For stable conformation, all frequencies must be > 0. If one frequency is < 0, you have arrived at a transition state, not a stable molecule. 5) Adjust - 6) MO surfaces A frequency calculation is valid only if it uses the same theoretical method and basis set as the geometry optimization.

18 2/7/ # * * = 90 FOOCl has two different transition states: one anti (*) and one eclipsed (#).

19 2/7/ Basic WebMO and MOPAC 1) Build 2) Optimize 3) Symmetry 4) Vibrations 5) Adjust and re-optimize 6) MO surfaces

20 2/7/ Open X-O-O-Y job 1) Click New Job Using This Geometry 2) Choose Molecular Orbitals 3) Change Job name; keep other settings 4) Right Arrow yes 5) Open Completed job 6) Find HOMO, LUMO etc. in MO Table

21 2/7/ Use File, Preferences to adjust the appearance of MOs

22 2/7/ Basic Gaussian

23 2/7/ Gaussian software Review MOs Basis sets Performance comparison: system size and basis set choice Typical Gaussian jobs: Opt, Freq, MO, ESP, NBO, Transition metals

24 2/7/ On-line Gaussian 09 User s Reference

25 2/7/ Excellent website and book!

26 Increased accuracy (and calculation time) 2/7/ Levels of theory- Gaussian Correlated ab initio MO (MP2 & others) Density Functional Theory (DFT) (B3LYP ) Ab initio MO (Hartree-Fock HF) Semi-empirical MO (PM3, PM6) Molecular mechanics (MM+, OPLS, UFF, etc)

27 H 2 MOs: formed by linear combination of 1s atomic orbitals. 27 Electron density (0.008 e/å 3 ) node node anti-bonding MO H atom H atom bonding MO 1s AO 1s AO d = 0.76 Å d = 2.8 Å d =

28 28 In semi-empirical MO theory (PM6 and others), molecular orbitals are formed by linear combinations of s, (and on heavy atoms) p x, p y, and p z AOs. s p x p y p z The orbital functions on the constituent atoms are the basis set.

29 2/7/ Gaussian and other QM programs use larger basis sets. The most common are split-valence basis sets. Basis orbital for core electrons of 2 nd row atom contains 6 Gaussian terms. 6-31G(d) 2 nd row atoms have an additional set of six d-type orbitals. Inner valence shell orbital contains 3 Gaussian terms. Outer valence shell orbital contains 1 Gaussian. H 2 O has 19 orbitals in the 6-31G(d) basis. O: 1s + 2s,2p x,2p y,2p z + (2s,2p x,2p y,2p z ) + six d H: - 1s + (1s) H: - 1s + (1s)

30 2/7/ r e - n

31 Wall clock time 2/7/ Assuming optimized in ~20 iterations, these would require ~10 min using 12 proc s on Antec12 Tetracycline C 22 H 24 O 8 N 2 HF/3-21G: 336 basis orbs Optimization required 26 iterations in 10 min Hartree-Fock; Antec12; 12 processors

32 2/7/ Gaussian input (.com) file format - energy (single point) - optimization - vibrations - molecular orbital name the checkpoint file theory/basisset jobtype charge multiplicity molecular description

33 2/7/ Gaussian input (.com) file format - energy (single point) - optimization - opt + freq - molecular orbital You may wish to add terms to modify the default optimization.

34 2/7/ Gaussian input (.com) file format - energy (single point) - optimization - opt + freq - molecular orbital Optimization + Vibrations jobs run sequentially In Gaussian, these are often run together because it is important to prove that the final structure is an energy minimum.

35 2/7/ What to do in case of a Failed Gaussian job. ***Check the.log output file**** The most common reason is a typo when you manually edit the input file. opt-vtight instead of opt=vtight Another fairly common error is adding an extra H when building. Optimization (Opt) jobs fail most often because of slow convergence of the energy, where the default maximum number of iterations is reached. - in this case use opt=(maxcycle=50) Other tactics for failed optimizations: - optimize at lower level of theory, then re-start. - for DFT methods, use ultrafine grid. - if you encounter small oscillations, test vibrations of one structure it usually is a minimum w/o imaginary frequency

36 2/7/ Gaussian input (.com) file format - energy (single point) - optimization - opt + freq - molecular orbital Pop is population analysis, which gives the electron populations of all orbitals. MOs and electron density graphics require orbital densities.

37 2/7/

38 38 The electron density at a molecule s outer surface is about electrons/å 3. The outer surface is sometimes defined as the density at intermolecular contacts within a crystal. 0.2 electrons/å electrons/å electrons/å 3 Electron density ~ electrons/å 3 Electron density increases closer to nuclei, and gradually goes to zero at greater distances.

39 Electrostatic potential mapped onto the electron density surface 39 Nuc: Color at each point on the surface reflects work required to bring +1 unit of electric charge from infinity to that point. = sum of potential due to the nuclei (+) and potential due to the electron density (-). H +

40 2/7/ Natural Bond Orbitals (NBO) NBO re-mixes the calculated MOs into linear combinations, making localized orbitals similar to, *,, * orbitals. Do CH 3 NH 2 example. %NProcShared=2 %CHK=CH3NH2-631Gss-B3LYP-NBO.chk #N B3LYP/6-31G(d,p) SP GFINPUT POP=(FULL,NBORead)formcheck CH3NH2-631Gss-B3LYP-NBO 0 1 C N H H H H H $NBO FILE=NBODATA AOPNAO=W32 AOPNHO=W34 AOPNBO=W36 DMNAO=W82 DMNHO=W84 DMNBO=W86 FNAO=W92 FNHO=W94 FNBO=W96 $END

41 2/7/ Lone pair on N. This is the donor in the intramolecular donor-acceptor interaction. Or, in other words, hyperconjugation between N and CH 3 group. C-H* anti-bond. This is the acceptor orbital.

42 2/7/ Natural Bond Orbitals (NBOs) in GaussView

43 2/7/ Transition metals (and other larger atoms): Effective Core Potentials The issue: The high nuclear charge of a larger atoms pulls in inner core electrons and increases their speed to near speed of light. Hence this is a relativistic effect on the core electrons. Especially important for 2 nd and 3 rd row transition metals. This relativistic effect changes the energy and shape of inner orbitals, which in turn affects the valence shell electrons. Inner electrons are represented by pseudo potential functions, or effective core potentials. Mo(CO) 4 (CS) 2 /Sally/7235

44 2/7/ %NProcShared=8 %CHK=MoC6O4S2trans-lan2dz-B3-Opt.chk #N B3LYP/gen pseudo=read OPT MoC6O4S2trans-lan2dz-B3-Opt Gen keyword says describe basis sets down here. 0 1 Mo C O C O C S C O C S C O Pseudo=read keyword says check at end for the type of ECP basis to used for heavy atom Mo. Blank lines required here. C O S G(d,p) **** Mo 0 lanl2dz **** Mo 0 lanl2dz zero

45 2/7/ Your questions

John Keller Department of Chemistry & Biochemistry University of Alaska Fairbanks

John Keller Department of Chemistry & Biochemistry University of Alaska Fairbanks 10/15/2016 1 WebMO & Gaussian John Keller Department of Chemistry & Biochemistry University of Alaska Fairbanks Corrections and updates 9-5-2017 SCHEDULE 9-10 Intro and basic operation of WebMO and MOPAC

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

General Chemistry Lab Molecular Modeling

General Chemistry Lab Molecular Modeling PURPOSE The objectives of this experiment are PROCEDURE General Chemistry Lab Molecular Modeling To learn how to use molecular modeling software, a commonly used tool in chemical research and industry.

More information

Introductory WebMO Exercises

Introductory WebMO Exercises Introductory WebMO Exercises These directions assume no prior knowledge of e WebMO interface and provide detailed, click by click instructions on building molecules and setting up calculations. Use your

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

Lab #3: Choice of Theoretical Method

Lab #3: Choice of Theoretical Method Lab #3: Choice of Theoretical Method These directions assume the user is familiar with the WebMO interface and can build molecules, set up calculations, etc. Exercise 1 - Determine the Proton Affinity

More information

1 Introduction to Computational Chemistry (Spartan)

1 Introduction to Computational Chemistry (Spartan) 1 Introduction to Computational Chemistry (Spartan) Start Spartan by clicking Start / Programs / Spartan Then click File / New Exercise 1 Study of H-X-H Bond Angles (Suitable for general chemistry) Structure

More information

Calculating Bond Enthalpies of the Hydrides

Calculating Bond Enthalpies of the Hydrides Proposed Exercise for the General Chemistry Section of the Teaching with Cache Workbook: Calculating Bond Enthalpies of the Hydrides Contributed by James Foresman, Rachel Fogle, and Jeremy Beck, York College

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

Transition states and reaction paths

Transition states and reaction paths Transition states and reaction paths Lab 4 Theoretical background Transition state A transition structure is the molecular configuration that separates reactants and products. In a system with a single

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

CHEM 344 Molecular Modeling

CHEM 344 Molecular Modeling CHEM 344 Molecular Modeling The Use of Computational Chemistry to Support Experimental Organic Chemistry Part 1: Molecular Orbital Theory, Hybridization, & Formal Charge * all calculation data obtained

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

QUANTUM CHEMISTRY WITH GAUSSIAN : A VERY BRIEF INTRODUCTION (PART 2)

QUANTUM CHEMISTRY WITH GAUSSIAN : A VERY BRIEF INTRODUCTION (PART 2) QUANTUM CHEMISTRY WITH GAUSSIAN : A VERY BRIEF INTRODUCTION (PART 2) TARAS V. POGORELOV AND MIKE HALLOCK SCHOOL OF CHEMICAL SCIENCES, UIUC This tutorial continues introduction to Gaussian [2]. Here we

More information

Literature values: ΔH f, gas = % error Source: ΔH f, solid = % error. For comparison, your experimental value was ΔH f = phase:

Literature values: ΔH f, gas = % error Source: ΔH f, solid = % error. For comparison, your experimental value was ΔH f = phase: 1 Molecular Calculations Lab: Some guideline given at the bottom of page 3. 1. Use the semi-empirical AM1 method to calculate ΔH f for the compound you used in the heat of combustion experiment. Be sure

More information

Project 1 Report File: Chem4PB3_project_1_2017-solutions last changed: 02-Feb-2017

Project 1 Report File: Chem4PB3_project_1_2017-solutions last changed: 02-Feb-2017 Project 1 Report File: Chem4PB3_project_1_2017-solutions last changed: 02-Feb-2017 1. Formaldehyde comparison of results from different methods Yellow shaded boxes are closest to experimental Method #e-

More information

Figure 1: Transition State, Saddle Point, Reaction Pathway

Figure 1: Transition State, Saddle Point, Reaction Pathway Computational Chemistry Workshops West Ridge Research Building-UAF Campus 9:00am-4:00pm, Room 009 Electronic Structure - July 19-21, 2016 Molecular Dynamics - July 26-28, 2016 Potential Energy Surfaces

More information

Instructions for Using Spartan 14

Instructions for Using Spartan 14 Instructions for Using Spartan 14 Log in to the computer with your Colby ID and password. Click on the Spartan 14 icon in the dock at the bottom of your screen. I. Building Molecules Spartan has one main

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

Chem 105 Friday Bond length qualitative Infrared spectroscopy Bond and molecular polarity Polarity and geometry Bond energy

Chem 105 Friday Bond length qualitative Infrared spectroscopy Bond and molecular polarity Polarity and geometry Bond energy 11/7/2011 1 Chem 105 Friday 11-7-11 Bond length qualitative Infrared spectroscopy Bond and molecular polarity Polarity and geometry Bond energy Hour Exam 3 on Chap 6, 7, and 8 is this Friday, 11/11/11.

More information

Modeling the UV-Vis Absorption of a Series of Dyes CH342L: Spectroscopy February 15, 2016

Modeling the UV-Vis Absorption of a Series of Dyes CH342L: Spectroscopy February 15, 2016 Modeling the UV-Vis Absorption of a Series of Dyes CH342L: Spectroscopy February 15, 2016 We ll correlate the absorbance maximum of a series of dyes with structural changes between them 1. Chemicals absorb

More information

Name. Chem Organic Chemistry II Laboratory Exercise Molecular Modeling Part 2

Name. Chem Organic Chemistry II Laboratory Exercise Molecular Modeling Part 2 Name Chem 322 - Organic Chemistry II Laboratory Exercise Molecular Modeling Part 2 Click on Titan in the Start menu. When it boots, click on the right corner to make the window full-screen. icon in the

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

Learning to Use Scigress Wagner, Eugene P. (revised May 15, 2018)

Learning to Use Scigress Wagner, Eugene P. (revised May 15, 2018) Learning to Use Scigress Wagner, Eugene P. (revised May 15, 2018) Abstract Students are introduced to basic features of Scigress by building molecules and performing calculations on them using semi-empirical

More information

Gaussian: Basic Tutorial

Gaussian: Basic Tutorial Input file: # hf sto-g pop=full Water - Single Point Energy 0 H.0 H.0 H 04.5 Route Section Start with # Contains the keywords Gaussian: Basic Tutorial Route Section Title Section Charge-Multiplicity Molecule

More information

Calculating NMR Chemical Shifts for beta-ionone O

Calculating NMR Chemical Shifts for beta-ionone O Calculating NMR Chemical Shifts for beta-ionone O Molecular orbital calculations can be used to get good estimates for chemical shifts. In this exercise we will calculate the chemical shifts for beta-ionone.

More information

Chemistry 14CL. Worksheet for the Molecular Modeling Workshop. (Revised FULL Version 2012 J.W. Pang) (Modified A. A. Russell)

Chemistry 14CL. Worksheet for the Molecular Modeling Workshop. (Revised FULL Version 2012 J.W. Pang) (Modified A. A. Russell) Chemistry 14CL Worksheet for the Molecular Modeling Workshop (Revised FULL Version 2012 J.W. Pang) (Modified A. A. Russell) Structure of the Molecular Modeling Assignment The molecular modeling assignment

More information

Conformational energy analysis

Conformational energy analysis Lab 3 Conformational energy analysis Objective This computational project deals with molecular conformations the spatial arrangement of atoms of molecules. Conformations are determined by energy, so the

More information

Tutorial I: IQ MOL and Basic DFT and MP2 Calculations 1 / 30

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

Gustavus Adolphus College. Lab #5: Computational Chemistry

Gustavus Adolphus College. Lab #5: Computational Chemistry CHE 372 Gustavus Adolphus College Lab #5: Computational Chemistry Introduction In this investigation we will apply the techniques of computational chemistry to several of the molecular systems that we

More information

Exercises for Windows

Exercises for Windows Exercises for Windows CAChe User Interface for Windows Select tool Application window Document window (workspace) Style bar Tool palette Select entire molecule Select Similar Group Select Atom tool Rotate

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

Introduc)on to IQmol: Part I.!!! Shirin Faraji, Ilya Kaliman, and Anna Krylov

Introduc)on to IQmol: Part I.!!! Shirin Faraji, Ilya Kaliman, and Anna Krylov Introduc)on to IQmol: Part I!!! Shirin Faraji, Ilya Kaliman, and Anna Krylov! 1 Resources! Written by Dr. Andrew Gilbert Keep yourself up to date with IQmol website: http://iqmol.org! IQmol Youtube channel:

More information

Excited States Calculations for Protonated PAHs

Excited States Calculations for Protonated PAHs 52 Chapter 3 Excited States Calculations for Protonated PAHs 3.1 Introduction Protonated PAHs are closed shell ions. Their electronic structure should therefore be similar to that of neutral PAHs, but

More information

Jaguar DFT Optimizations and Transition State Searches

Jaguar DFT Optimizations and Transition State Searches Jaguar DFT Optimizations and Transition State Searches Density Functional Theory (DFT) is a quantum mechanical (QM) method that gives results superior to Hartree Fock (HF) in less computational time. A

More information

NH 3 inversion: Potential energy surfaces and transition states CH342L March 28, 2016

NH 3 inversion: Potential energy surfaces and transition states CH342L March 28, 2016 N 3 inversion: Potential energy surfaces and transition states C342L March 28, 2016 Last week, we used the IR spectrum of ammonia to determine the splitting of energy levels due to inversion of the umbrella

More information

The Hückel Approximation Consider a conjugated molecule i.e. a molecule with alternating double and single bonds, as shown in Figure 1.

The Hückel Approximation Consider a conjugated molecule i.e. a molecule with alternating double and single bonds, as shown in Figure 1. The Hückel Approximation In this exercise you will use a program called Hückel to look at the p molecular orbitals in conjugated molecules. The program calculates the energies and shapes of p (pi) molecular

More information

Ethene. Introduction. The ethene molecule is planar (i.e. all the six atoms lie in the same plane) and has a high degree of symmetry:

Ethene. Introduction. The ethene molecule is planar (i.e. all the six atoms lie in the same plane) and has a high degree of symmetry: FY1006 Innføring i kvantefysikk og TFY4215 Kjemisk fysikk og kvantemekanikk Spring 2012 Chemical Physics Exercise 1 To be delivered by Friday 27.04.12 Introduction Ethene. Ethylene, C 2 H 4, or ethene,

More information

Basis Set for Molecular Orbital Theory

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

Appendix D Simulating Spectroscopic Bands Using Gaussian and PGopher

Appendix D Simulating Spectroscopic Bands Using Gaussian and PGopher 429 Appendix D Simulating Spectroscopic Bands Using Gaussian and PGopher This appendix contains methods for using Gaussian 09 121 and PGopher 120 to simulate vibrational and electronic bands of molecules.

More information

Be H. Delocalized Bonding. Localized Bonding. σ 2. σ 1. Two (sp-1s) Be-H σ bonds. The two σ bonding MO s in BeH 2. MO diagram for BeH 2

Be H. Delocalized Bonding. Localized Bonding. σ 2. σ 1. Two (sp-1s) Be-H σ bonds. The two σ bonding MO s in BeH 2. MO diagram for BeH 2 The Delocalized Approach to Bonding: The localized models for bonding we have examined (Lewis and VBT) assume that all electrons are restricted to specific bonds between atoms or in lone pairs. In contrast,

More information

one ν im: transition state saddle point

one ν im: transition state saddle point Hypothetical Potential Energy Surface Ethane conformations Hartree-Fock theory, basis set stationary points all ν s >0: minimum eclipsed one ν im: transition state saddle point multiple ν im: hilltop 1

More information

Lewis Structures and Molecular Shapes

Lewis Structures and Molecular Shapes Lewis Structures and Molecular Shapes Rules for Writing Lewis Structures 1. Determine the correct skeleton structure (connectivity of atoms). Usually, the most electronegative atoms go around the edges

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

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

Theoretical UV/VIS Spectroscopy

Theoretical UV/VIS Spectroscopy Theoretical UV/VIS Spectroscopy Why is a Ruby Red When Chromium Oxide is Green? How Does a Ruby Laser Work? Goals of this Exercise: - Calculation of the energy of electronically excited states - Understanding

More information

Introduction to Ab Initio Quantum Chemical Computation

Introduction to Ab Initio Quantum Chemical Computation c:\374-17\computation\computation17.doc for 9mar17 Prof. Patrik Callis 8mar17 Introduction to Ab Initio Quantum Chemical Computation Purpose: 1. To become acquainted with basic concepts of ab initio quantum

More information

Chapter 5. Molecular Orbitals

Chapter 5. Molecular Orbitals Chapter 5. Molecular Orbitals MO from s, p, d, orbitals: - Fig.5.1, 5.2, 5.3 Homonuclear diatomic molecules: - Fig. 5.7 - Para- vs. Diamagnetic Heteronuclear diatomic molecules: - Fig. 5.14 - ex. CO Hybrid

More information

LUMO + 1 LUMO. Tómas Arnar Guðmundsson Report 2 Reikniefnafræði G

LUMO + 1 LUMO. Tómas Arnar Guðmundsson Report 2 Reikniefnafræði G Q1: Display all the MOs for N2 in your report and classify each one of them as bonding, antibonding or non-bonding, and say whether the symmetry of the orbital is σ or π. Sketch a molecular orbital diagram

More information

Department of Chemistry. The Intersection of Computational Chemistry and Experiment

Department of Chemistry. The Intersection of Computational Chemistry and Experiment Department of Chemistry The Intersection of Computational Chemistry and Experiment Structure, Vibrational and Electronic Spectra of Organic Molecules Angelo R. Rossi Department of Chemistry The University

More information

XYZ file format Protein Data Bank (pdb) file format Z - matrix

XYZ file format Protein Data Bank (pdb) file format Z - matrix Chemistry block (exercise 1) In this exercise, students will be introduced how to preform simple quantum chemical calculations. Input files for Gaussian09. Output file structure. Geometry optimization,

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

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

Molecular Orbital Theory This means that the coefficients in the MO will not be the same!

Molecular Orbital Theory This means that the coefficients in the MO will not be the same! Diatomic molecules: Heteronuclear molecules In heteronuclear diatomic molecules, the relative contribution of atomic orbitals to each MO is not equal. Some MO s will have more contribution from AO s on

More information

The Potential Energy Surface

The Potential Energy Surface The Potential Energy Surface In this section we will explore the information that can be obtained by solving the Schrödinger equation for a molecule, or series of molecules. Of course, the accuracy of

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

Be H. Delocalized Bonding. Localized Bonding. σ 2. σ 1. Two (sp-1s) Be-H σ bonds. The two σ bonding MO s in BeH 2. MO diagram for BeH 2

Be H. Delocalized Bonding. Localized Bonding. σ 2. σ 1. Two (sp-1s) Be-H σ bonds. The two σ bonding MO s in BeH 2. MO diagram for BeH 2 The Delocalized Approach to Bonding: The localized models for bonding we have examined (Lewis and VBT) assume that all electrons are restricted to specific bonds between atoms or in lone pairs. In contrast,

More information

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

Same idea for polyatomics, keep track of identical atom e.g. NH 3 consider only valence electrons F(2s,2p) H(1s)

Same idea for polyatomics, keep track of identical atom e.g. NH 3 consider only valence electrons F(2s,2p) H(1s) XIII 63 Polyatomic bonding -09 -mod, Notes (13) Engel 16-17 Balance: nuclear repulsion, positive e-n attraction, neg. united atom AO ε i applies to all bonding, just more nuclei repulsion biggest at low

More information

Project 2. Chemistry of Transient Species in Planetary Atmospheres: Exploring the Potential Energy Surfaces of CH 2 S

Project 2. Chemistry of Transient Species in Planetary Atmospheres: Exploring the Potential Energy Surfaces of CH 2 S Chemistry 362 Spring 2018 Dr. Jean M. Standard March 21, 2018 Project 2. Chemistry of Transient Species in Planetary Atmospheres: Exploring the Potential Energy Surfaces of CH 2 S In this project, you

More information

Patrick: An Introduction to Medicinal Chemistry 5e MOLECULAR MODELLING EXERCISES CHAPTER 17

Patrick: An Introduction to Medicinal Chemistry 5e MOLECULAR MODELLING EXERCISES CHAPTER 17 MOLECULAR MODELLING EXERCISES CHAPTER 17 Exercise 17.6 Conformational analysis of n-butane Introduction Figure 1 Butane Me Me In this exercise, we will consider the possible stable conformations of butane

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

Practical 1: Structure and electronic properties of organic molecules. B/ Structure, electronic and vibrational properties of the water molecule

Practical 1: Structure and electronic properties of organic molecules. B/ Structure, electronic and vibrational properties of the water molecule D1CH9116 - MMCO Molecular Modelling applied to organic chemistry Practical 1: tructure and electronic properties of organic molecules B/ tructure, electronic and vibrational properties of the water molecule

More information

with the larger dimerization energy also exhibits the larger structural changes.

with the larger dimerization energy also exhibits the larger structural changes. A7. Looking at the image and table provided below, it is apparent that the monomer and dimer are structurally almost identical. Although angular and dihedral data were not included, these data are also

More information

A Computer Study of Molecular Electronic Structure

A Computer Study of Molecular Electronic Structure A Computer Study of Molecular Electronic Structure The following exercises are designed to give you a brief introduction to some of the types of information that are now readily accessible from electronic

More information

Example: H 2 O (the car file)

Example: H 2 O (the car file) Example: H 2 O (the car file) As a practical example of DFT methods we calculate the energy and electronic properties of the water molecule. In order to carry out the DFT calculation you will need a set

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

The symmetry properties & relative energies of atomic orbitals determine how they react to form molecular orbitals. These molecular orbitals are then

The symmetry properties & relative energies of atomic orbitals determine how they react to form molecular orbitals. These molecular orbitals are then 1 The symmetry properties & relative energies of atomic orbitals determine how they react to form molecular orbitals. These molecular orbitals are then filled with the available electrons according to

More information

Contents. 1. Building a molecular model 2. Additional Procedures and Options: 3. Selecting a calculation method

Contents. 1. Building a molecular model 2. Additional Procedures and Options: 3. Selecting a calculation method Getting started with HyperChem. Program description. In this project you will create simple organic molecules and learn how to manipulate them. However, before starting you should create your home directory

More information

Lab #5: Electron Densities, Electrostatic Potentials, and Reactivity Indices

Lab #5: Electron Densities, Electrostatic Potentials, and Reactivity Indices Lab #5: Electron Densities, Electrostatic Potentials, and Reactivity Indices Exercise 1 - Visualizing Different Bond Types Build H 2 and perform a geometry optimization (Mopac) using the choices shown

More information

The Hückel Approximation

The Hückel Approximation The ückel Approximation 1 In this exercise you will use a program called ückel to look at the π molecular orbitals in conjugated molecules. The program calculates the energies and shapes of π (pi) 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

Theoretical Chemistry - Level II - Practical Class Molecular Orbitals in Diatomics

Theoretical Chemistry - Level II - Practical Class Molecular Orbitals in Diatomics Theoretical Chemistry - Level II - Practical Class Molecular Orbitals in Diatomics Problem 1 Draw molecular orbital diagrams for O 2 and O 2 +. E / ev dioxygen molecule, O 2 dioxygenyl cation, O 2 + 25

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

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

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

Computational chemistry with GAMESS: a very brief overview with examples

Computational chemistry with GAMESS: a very brief overview with examples Computational chemistry with GAMESS: a very brief overview with examples PHY-6120 Molecular Physics (Spring 2015), UConn Phys. Dept. Feb 17 th 2015 H = ħ2 2μ i Intro: V(R) for diatomic molecules + k Z

More information

Project 3: Molecular Orbital Calculations of Diatomic Molecules. This project is worth 30 points and is due on Wednesday, May 2, 2018.

Project 3: Molecular Orbital Calculations of Diatomic Molecules. This project is worth 30 points and is due on Wednesday, May 2, 2018. Chemistry 362 Spring 2018 Dr. Jean M. Standard April 20, 2018 Project 3: Molecular Orbital Calculations of Diatomic Molecules In this project, you will investigate the molecular orbitals and molecular

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

Computational Chemistry Problem Set 4 Due Monday, February 22, 2011 (at the start of class) Total Number of Points = 85

Computational Chemistry Problem Set 4 Due Monday, February 22, 2011 (at the start of class) Total Number of Points = 85 Computational Chemistry Problem Set 4 Due Monday, February 22, 2011 (at the start of class) Total Number of Points = 85 Basic Technical Notes: (1) For security reasons, you are allowed to log into the

More information

CHEM J-5 June 2014

CHEM J-5 June 2014 CHEM1101 2014-J-5 June 2014 The molecular orbital energy level diagrams for H 2, H 2 +, H 2 and O 2 are shown below. Fill in the valence electrons for each species in its ground state and label the types

More information

Molecular Modeling and Conformational Analysis with PC Spartan

Molecular Modeling and Conformational Analysis with PC Spartan Molecular Modeling and Conformational Analysis with PC Spartan Introduction Molecular modeling can be done in a variety of ways, from using simple hand-held models to doing sophisticated calculations on

More information

Introduction to Computational Chemistry Exercise 2

Introduction to Computational Chemistry Exercise 2 Introduction to Computational Chemistry Exercise 2 Intermolecular interactions and vibrational motion Lecturer: Antti Lignell Name Introduction In this computer exercise, we model intermolecular interactions

More information

Basis Sets and Basis Set Notation

Basis Sets and Basis Set Notation Chemistry 46 Fall 215 Dr. Jean M. Standard November 29, 217 Basis Sets and Basis Set Notation Using the LCAO-MO approximation, molecular orbitals can be represented as linear combinations of atomic orbitals,

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

Graphical User Interface Simplified Graphics Requests. The most commonly requested graphics (HOMO,

Graphical User Interface Simplified Graphics Requests. The most commonly requested graphics (HOMO, Spartan 10 ( Spartan ) has been designed to address the ever increasing role that calculations play in chemistry and related fields. It represents a continued collaboration between Wavefunction, Inc.,

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

Lewis Structures and Molecular Shapes

Lewis Structures and Molecular Shapes Lewis Structures and Molecular Shapes Rules for Writing Lewis Structures 1. Determine the correct skeleton structure (connectivity of atoms). Usually, the most electronegative atoms go around the edges

More information

Molecular Modeling 2: Computational Chemistry

Molecular Modeling 2: Computational Chemistry Molecular Modeling 2: Computational Chemistry Authors: J. M. McCormick, E. V. Patterson and M. C. Nagan* Update: December 29, 2011 and July 3, 2012 by V. M. Pultz Introduction Computational chemists use

More information

Quick Tutorial on Natural Bond Order 3 Calculations Within Gaussian 09

Quick Tutorial on Natural Bond Order 3 Calculations Within Gaussian 09 Quick Tutorial on Natural Bond Order 3 Calculations Within Gaussian 09 Benjamin Rudshteyn, and Victor S. Batista* victor.batista@yale.edu Yale University, Department of Chemistry 225 Prospect Street, New

More information

The use of solvation models and the ONIOM layered approach in Gaussian.

The use of solvation models and the ONIOM layered approach in Gaussian. The use of solvation models and the NIM layered approach in Gaussian. In this lab we will consider two techniques that are very useful to model larger systems: the use of solvation models to mimic systems

More information

Homework Problem Set 5 Solutions. E e + H corr (a.u.) a.) Determine the bond dissociation enthalpy of ethane in kcal/mol (1 a.u. = kcal/mol).

Homework Problem Set 5 Solutions. E e + H corr (a.u.) a.) Determine the bond dissociation enthalpy of ethane in kcal/mol (1 a.u. = kcal/mol). Chemistry 380.37 Dr. Jean M. Standard Homework Problem Set 5 Solutions 1. Given below are the sum of electronic and thermal enthalpies, E e + H corr, from Hartree-Fock calculations using a 6-31G(d) basis

More information

Physical Chemistry II Recommended Problems Chapter 12( 23)

Physical Chemistry II Recommended Problems Chapter 12( 23) Physical Chemistry II Recommended Problems Chapter 1( 3) Chapter 1(3) Problem. Overlap of two 1s orbitals recprobchap1.odt 1 Physical Chemistry II Recommended Problems, Chapter 1(3), continued Chapter

More information

Experiment 15: Atomic Orbitals, Bond Length, and Molecular Orbitals

Experiment 15: Atomic Orbitals, Bond Length, and Molecular Orbitals Experiment 15: Atomic Orbitals, Bond Length, and Molecular Orbitals Introduction Molecular orbitals result from the mixing of atomic orbitals that overlap during the bonding process allowing the delocalization

More information

In this lecture we will understand how the molecular orbitals are formed from the interaction of atomic orbitals.

In this lecture we will understand how the molecular orbitals are formed from the interaction of atomic orbitals. Lecture 7 Title: Understanding of Molecular Orbital Page-1 In this lecture we will understand how the molecular orbitals are formed from the interaction of atomic orbitals. We will see how the electrons

More information

First Semester Organic Chemistry Five Sets of Practice Exams

First Semester Organic Chemistry Five Sets of Practice Exams First Semester Organic Chemistry Five Sets of Practice Exams 2000-2006 N Cl Prof. J. W. Keller i Cover: Molecular orbital (MO) description of the transition state of the nucleophilic substitution reaction

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

Lecture 9 Electronic Spectroscopy

Lecture 9 Electronic Spectroscopy Lecture 9 Electronic Spectroscopy Molecular Orbital Theory: A Review - LCAO approximaton & AO overlap - Variation Principle & Secular Determinant - Homonuclear Diatomic MOs - Energy Levels, Bond Order

More information

Assignment A02: Geometry Definition: File Formats, Redundant Coordinates, PES Scans

Assignment A02: Geometry Definition: File Formats, Redundant Coordinates, PES Scans Assignment A02: Geometry Definition: File Formats, Redundant Coordinates, PES Scans In Assignments A00 and A01, you familiarized yourself with GaussView and G09W, you learned the basics about input (GJF)

More information

Lecture 4 Model Answers to Problems

Lecture 4 Model Answers to Problems Self-Study Problems / Exam Preparation Construct and annotate a valence MO diagram for H 2 CN -. Use your diagram to explain why the neutral radical is more stable than the anion. (this is an old exam

More information