Gaussian Basis Sets for Solid-State Calculations
|
|
- Alyson Hood
- 5 years ago
- Views:
Transcription
1 Gaussian Basis Sets for Solid-State Calculations K. Doll Molpro Quantum Chemistry Software Institute of Theoretical Chemistry, D Stuttgart, Germany MW-MSSC 2017, Minneapolis, July 10, 2017
2 Introduction about basis sets Outline How to make up your own basis set - all-electron example: carbon (for diamond) - ECP example: Sr (for SrO) - Miscellaneous hints: reoptimising exponents, exponents for free atom, brute force approach (even tempered), how to deal with error messages Note: Many, many numbers in this talk catch main ideas
3 Introduction
4 Basis set types Hartree-Fock/Kohn-Sham equations can in general not be solved analytically => expand molecular/crystalline orbitals in basis functions plane waves: ħ 2 ( k + K ) 2 E m cutoff el Φ k ( r )= c k, K exp( i( k+ K ) r) K with: ħ 2 2m el ( k+ K ) 2 E cutoff local basis sets (linear combination of atomic orbitals, LCAO): - Slater type orbitals (STO): ϕ STO ( r)=n x a y b z c exp( ζ r) solve H-atom exactly, core region properly described - Gaussian type orbitals (GTO): integrals easier than for STOs ϕ GTO ( r )= N x a y b z c exp( ζ r 2 )
5 Gaussian basis sets ϕ primitive ( r )=N (x A x ) a ( y A y ) b ( z A z ) c exp ( ζ ( r A ) 2 ) primitive basis function, atom-centered, atom at A Exponent: ζ: unit: 1/bohr 2 1 bohr = Å large ζ: tight important for core region Small ζ ( < 1 ) : diffuse, important for valence region, chemical bond
6 Contraction Idea: the core region does not change very much when a chemical bond is formed => make linear combination which describe core region well: contraction : n ϕ( r )= N i=1 c i x a y b z c exp( ζ i r 2 ) primitive Gaussian c i : contraction coefficients Total basis set: set of various ϕ μ ( r ) with different exponents, contractions coefficients
7 Orbitals Molecular orbital: Ψ i ( r )= μ a i μ ϕ μ ( r ) Crystalline orbital: Ψ i ( r, k )= μ a iμ Φ μ ( r, k) with: Φ μ ( r, k)= g ϕ μ ( r A μ g)exp(i k g) The orbitals are obtained when the Hartree-Fock/Kohn-Sham equations are solved. All basis functions contribute to each orbital.
8 Angular part, 1 Y lm (θ,φ) : spherical harmonics r l Y lm (θ,φ) : solid harmonics L z eigenstates only for atom or linear molecule => real combinations are preferably used (Condon-Shortley convention, normalisation ignored in the following!): l=0 : Y 00 => 1 l=1 : Y 11 - Y 1-1 => x Y 11 + Y 1-1 => y Y 10 => z
9 Angular part, 2 l=2 : Y 20 => 3z 2 -r 2 =2z 2 -x 2 -y 2 Y 21 Y 2-1 => xz Y 21 + Y 2-1 => yz Y 22 + Y 2-2 => x 2 -y 2 Y 22 Y 2-2 => xy The order is in the manual, page 22
10 Angular part, 3 Why does the order matter? Population analysis (keyword PPAN), projected density of states (keyword DOSS): this is how the orbitals are ordered Typical application of population analysis: Ion embedded in crystal field of other ions (octahedral, tetrahedral ) what is the d-occupancy? Use ROTREF if surrounding ions are not textbook-like arranged
11 Exponents: where from? Easiest: even tempered: exponents separated by common factor Example: factor 4, start from 0.25, stop at 4096: 0.25, 1, 4, 16, 64, 256, 1024, 4096 general: αβ k here: α=0.25 β=4 k=0...7 must be done for the different angular momenta: α s, β s, α p, β p, α d, β d Optimise α, β (lowest total energy of the respective atom)
12 Even tempered basis sets: optimise parameters α, β αβ k Hartree-Fock limit can be reached with more and more basis functions For practical purposes: use exponents separated by factors of 2-4
13 Basis set choice: databases Molecular case: databases available EMSL basis set exchange: MOLPRO code: Huge number of basis sets available Periodic systems: CRYSTAL website Basis sets for most atoms available, see also CRYSTAL manual chapter 12, page
14 Basis set choice: general recommendations Type of calculation is important: here Hartree-Fock or DFT, but no wavefunction based correlation (CCSD, CI, MP2 ) => relatively small basis set is sufficient, basis set does not have to be correlation-consistent Basis set depends on charge state (0, -1, +1?): negative ions have larger radius than positive ions => extra diffuse exponents required for negatively charged atoms (and metals in the periodic case), may be removed for positively charged atoms usually good to include polarisation function: e.g. if s and p orbitals are occupied, add d-exponents effective core potentials may be advantageous from 1st row transition metals on (to reduce calculational effort + include scalar-relativistic effects) but you can still do all-electron calculations even for heavy atoms (e.g. to look at core-levels)
15 Basis set: general recommendations adding more basis functions lowers the energy because of the variational principle! diffuse functions are important for the chemical bond very large exponents may improve the core region lead to lower energy, but only little improvement of the region where the chemical bond takes place => usually the part of the basis set which describes the region close to the nucleus can be kept as is diffuse functions may have to be altered (e.g. use educated guess such as even tempered basis functions, or reoptimise)
16 Basis set: personal favorites for light elements: a basis set similar to a molecular 6-31G basis set, usually already on CRYSTAL basis set website for heavy elements: Stuttgart/Cologne effective core potentials, diffuse exponents must be reoptimised
17 Basis set names huge confusion due to enormous number of basis sets, and enormous number of names most general: (number of s-, p-, d- exponents) / [number of s-, p-, d- contractions] e.g. (10s4p) / [3s2p]
18 Is my basis set good enough? Try it! If possible, compare with enhanced basis set (e.g. add diffuse function, polarisation function...): total energy, band structure, population analysis How large is the Mulliken population of the diffuse functions? Compute the band structure for a simple system where the element has a similar charge state - when it works for a simple system, it should work for a bigger system with similar charge state
19 Possible problems with Gaussian basis sets Very diffuse basis functions can cause numerical instabilities: => avoid exponents lower than for periodic systems (molecular codes use test on eigenvalues of overlap matrix and reduce basis set) But neither use a basis set which was calibrated for a positively charged ion for the neutral atom or negatively charged ion it needs extra diffuse functions e.g. a basis set calibrated for Li + in LiF may not work for Li metal without extra diffuse functions Exponents should be separated by factor > 2 remember even tempered basis sets: factor 2-4 between exponents
20 In the following: How to make up your own basis set
21 Example: Diamond
22 Example: Diamond, choose from EMSL web site
23 Example: Diamond Neutral => try standard basis set for molecules, e.g. 6-31G* * : polarisation function (here: d-shell) From EMSL website: 6-31G* basis set BASIS "ao basis" PRINT #BASIS SET: (10s,4p,1d) -> [3s,2p,1d] C S C SP C SP C D END s-shell sp-shell sp-shell d-shell exponent s-contraction p-contraction
24 Diamond quick analysis of the basis set From EMSL website: BASIS "ao basis" PRINT #BASIS SET: (10s,4p,1d) -> [3s,2p,1d] C S C SP C SP C D END Quick analysis: most diffuse exponent is => not too diffuse, try unchanged 6-31G* should be able to describe neutral carbon atom in molecular calculation => should be alright for diamond where carbon is neutral => try unchanged! exponent s-contraction p-contraction
25 Example: Diamond from EMSL to CRYSTAL format From EMSL website: BASIS "ao basis" PRINT #BASIS SET: (10s,4p,1d) -> [3s,2p,1d] C S C SP C SP C D END CRYSTAL format: This basis set had been previously used in CRYSTAL calculations and is on the web at
26 Diamond band structure: crystal input diamond.d12 part 1 part 2 Diamond 6-31G* basis set CRYSTAL END END DFT EXCHANGE LDA CORRELAT PZ END NODIRECT SHRINK 8 16 MAXCYCLE 60 FMIXING 70 NOSHIFT PPAN END Full input1 = put part1 + part2 together in 1 file: diamond.d12
27 Diamond band structure: properties input diamond.d3 BAND diamond band structure G-X X-W W-L L-G G-X END agrees well with H. Bross, R. Bader, phys. stat. sol. (b) 191, 369 (1995) Deviations possible for very high bands (note that we have only as many bands as basis functions!) => basis set reasonable for carbon in neutral, positive, or slightly negative state with a good band structure, the other properties are usually also reasonable, e.g. equilibrium geometry
28 ECP-Input and example: SrO
29 Effective core potentials Idea: reduce computational effort with a potential which describes nucleus and inner electrons + include a part of the relativistic effects Shape-consistent: fit shape of valence orbitals to all-electron shape Energy-consistent: fit to a set of energies of electronic states of an atom/ion In the following: Stuttgart/Cologne energy-consistent (ab initio) pseudopotentials
30 General recommendation Labeled as: ECPnXY ECP : effective core potential n: number of core electrons X: S or M S: single-valence-electron ion used in fit (use this only for one- or two valence electron atoms) M: multi-valence-electron ion used in fit (use M if available) Y: HF Hartree-Fock (only for light elements) WB Wood-Boring (simplification of DF) DF Dirac-Fock (fully relativistic, 4 components as in Dirac equation) use WB or DF if available e.g. ECP10MWB : nucleus + 1s 2 2s 2 2p 6 replaced by pseudopotential ( Ne core, 10 electrons) Multi-valence-electron atom used in fit, relativistic
31 Example: SrO O: use from CRYSTAL web site, here: O_8-411_towler_1994 Sr (38 electrons): select at Click on Sr, new window opens Options: ECP36SHF ECP36SDF ECP28MHF ECP28MWB ECP28MDF } } large core (36 electrons in core) small core (28 electrons in core) Let us choose ECP28MWB for the effective core potential and for the basis set also the one labeled ECP28MWB
32 Example: SrO ECP main web site
33 ECP for Sr Pseudopotential ECP28MWB for Sr # COMMENT # Q=10., MEFIT, WB, Ref 13. # COMMENT # Sr 238 <here, the number of valence-basis shells is needed> INPUT <please add the valence-basis input, without blank line> # COMMENT # References: # COMMENT # [13] M. Kaupp, P. v. R. Schleyer, H. Stoll, H. Preuss, J. Chem. Phys. 94, 1360 (1991). Available basis sets for Sr, ECP28MWB: ECP28MWB ECP28MDF_GUESS (ECP28MHF_GUESS ECP28MWB_GUESS) Core: 1s 2 2s 2 2p 6 3s 2 3p 6 3d 10 : 28 electrons in core Use ECP data as is no change when used with CRYSTAL
34 Basis set for Sr, compatible with ECP Basis set ECP28MWB for Sr, ECP28MWB # COMMENT # (6s6p5d)/[4s4p2d)-Basissatz fuer PP. von Ref X X X X X X X X X X X: here, the number of electrons per valence shell is needed CRYSTAL format: replace X with number X=2 electrons X=2 electrons X=0 electrons X=0 electrons X=6 electrons X=0 electrons X=0 electrons X=0 electrons X=0 electrons X=0 electrons 4s 2 4p 6 5s 2 (neutral)
35 Basis set for Sr, compatible with ECP Quick analysis: Many diffuse functions < 0.1 Will make problems => remove But then the outermost diffuse is for s for p for d not diffuse enough? => use educated guess, between Factor 2-4 smaller than smallest exponent e.g.: add extra diffuse s 0.1 p 0.1 d 0.1
36 Basis set: summary of the strategy for Sr Start from the initial basis set Remove diffuse exponents, e.g. exponents < 0.1 Fill up with new diffuse exponents, but only if by at least factor 2 smaller than smallest exponent Test!
37 Basis set for Sr, compatible with ECP Full basis set Remove diffuse Exponents, e.g. < 0.1 too diffuse too diffuse too diffuse too diffuse too diffuse Educated guess: add diffuse s,p,d with exponent 0.1 each (by at least factor 2 smaller than smallest exponent) Final basis set
38 Input for Sr and O INPUT :Sr, ECP 8 shells Effective core potential 10 valence electrons Sr: O: Basis set: s-shell, 2 electrons s-shell, 2 electrons s-shell, 0 electrons p-shell, 6 electrons p-shell, 0 electrons p-shell, 0 electrons d-shell, 0 electrons d-shell, 0 electrons 4s 2 4p 6 5s 2 : 10 electrons from Turin web page Starting configuration: 1s 2 2s 2 2p 4 : 8 electrons
39 SrO CRYSTAL END INPUT SrO: input for crystal sro.d12 Part 1 part 2 part 3 part 4 geometry, ECP Sr basis set O basis set SCF part END Put all 4 parts together in 1 file: sro.d12 DFT EXCHANGE LDA CORRELAT PZ END NODIRECT SHRINK 8 16 MAXCYCLE 60 FMIXING 70 NOSHIFT PPAN END
40 Test: SrO band structure, input sro.d3 BAND Input for fcc structure L-G G-X X-W W-K K-G END Good agreement with B. Pan, N.-P. Wang, M. Rohlfing, Appl. Phys. A (2015)
41 SrO: possible Sr basis set deficiencies full Sr basis set diffuse s omitted diffuse p omitted diffuse d omitted
42 SrO: possible Sr basis set deficiencies neutral Sr has configuration 5s 2 lowest unoccupied orbitals are 4d 5p => diffuse functions important for s, d ; less for p
43 diffuse p omitted diffuse d omitted SrO: possible Sr basis set deficiencies note: configuration of Sr 2+ is: 4s 2 4p 6 5s 0 4d 0 5p 0 full Sr basis set: ECP+[3s3p2d] diffuse s omitted
44 Miscellaneous hints
45 Improve an existing basis set: reoptimise exponents Assume we want to do calculation for Cu bulk Existing copper basis set, calibrated for KCuF 3, Cu 2+ state M.D. Towler, R. Dovesi and V.R. Saunders Phys. Rev. B 52, (1995) Quick analysis: sp is not diffuse enough for Cu metal typical values are sp d exponent? => educated guess: sp as initial guess then optimise: sp d criterion: lowest energy
46 Improve an existing basis set: reoptimise exponents Optimise sp d 0.430: Perform total energy calculation for Cu bulk, at experimental lattice constant, optimise iteratively: Initial energy (PWGGA) I 1) optimise sp 0.6 => 0.55 (keep other exponents fixed) E E+03 2) optimise sp 0.2 => E+03 3) optimise d 0.43 => E+03 II 1) optimise sp 0.55 => E+03 2) optimise sp 0.14 => 0.14 same energy 3) optimise d 0.40 => 0.40 same energy Done! Finally we have optimised: sp and d 0.40 Energy gain: about hartree band structure Thursday
47 Basis set for a single free atom Idea: linear dependence problem much less severe for a single free atom (keyword MOLECULE) => use basis set with many diffuse functions to describe the neutral atom (necessary for e.g. cohesive energy) but do not change the tight exponents, otherwise inconsistent with bulk!
48 Cu bulk Basis set for a free Cu atom Cu atom: extra diffuse functions
49 A basis set from scratch: even tempered for diamond If anything else fails, we can make up an even tempered basis set: Example for carbon, similar for other elements: factor ~3 between exponents LDA band structure with this basis set
50 Error message: number of electrons and nuclear charge do not match Necessary for periodic systems: number of electrons = nuclear charge Otherwise Coulomb energy diverges Error message is: ERROR **** INPBAS **** UNIT CELL NOT NEUTRAL remedy: properly count all the electrons in the basis set input simplest: make each atom neutral SrO example: Sr is neutral, O is neutral
51 Error message: linear dependence If this error message shows up: ERROR **** CHOLSK **** BASIS SET LINEARLY DEPENDENT Then most likely one or more exponents are too diffuse. Remedy: check all the diffuse exponents, remove small ones, possibly reoptimise remaining exponents
52 Summary Many good basis sets for periodic systems are already available on the web Charge state important, check diffuse exponents, use educated guess : - even tempered rule: factor 2-4 between exponents - avoid exponents below ~ compute band structure as test if basis set is large enough Optimise exponent: use total energy
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 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 informationThis 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 informationAb initio structure prediction for molecules and solids
Ab initio structure prediction for molecules and solids Klaus Doll Max-Planck-Institute for Solid State Research Stuttgart Chemnitz, June/July 2010 Contents structure prediction: 1) global search on potential
More informationOpen-shell and Magnetic Systems with CRYSTAL Tools
Open-shell and Magnetic Systems with CRYSTAL Tools Klaus Doll Molpro Quantum Chemistry Software Institute of Theoretical Chemistry, D-70569 Stuttgart, Germany MSSC 2018, Turin, September 2018 Contents
More informationGaussian: 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 informationComputational 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 informationMODELING MATTER AT NANOSCALES
MODELING MATTER AT NANOSCALES 6. The theory of molecular orbitals for the description of nanosystems (part II) 6.0. Ab initio methods. Basis functions. Luis A. Monte ro Firmado digitalmente por Luis A.
More informationBasis Set for Molecular Orbital Theory
Basis Set for Molecular Orbital Theory! Different Types of Basis Functions! Different Types of Atom Center Basis Functions! Classifications of Gaussian Basis Sets! Pseudopotentials! Molecular Properties
More informationModels and strategies for metallic systems
Models and strategies for metallic systems Klaus Doll Molpro Quantum Chemistry Software Institute of Theoretical Chemistry Pfaffenwaldring 55 D-70569 Stuttgart, Germany MW-MSSC 2017, Minneapolis, July
More informationDFT EXERCISES. FELIPE CERVANTES SODI January 2006
DFT EXERCISES FELIPE CERVANTES SODI January 2006 http://www.csanyi.net/wiki/space/dftexercises Dr. Gábor Csányi 1 Hydrogen atom Place a single H atom in the middle of a largish unit cell (start with a
More informationA One-Slide Summary of Quantum Mechanics
A One-Slide Summary of Quantum Mechanics Fundamental Postulate: O! = a! What is!?! is an oracle! operator wave function (scalar) observable Where does! come from?! is refined Variational Process H! = E!
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 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 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 informationAnswers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014
Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde
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 informationChem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 28, December 08, 2014
Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi Lecture 28, December 08, 2014 Solved Homework Water, H 2 O, involves 2 hydrogen atoms and an oxygen
More informationIntroduction 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 informationYingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2
Purdue University CHM 67300 Computational Quantum Chemistry REVIEW Yingwei Wang October 10, 2013 Review: Prof Slipchenko s class, Fall 2013 Contents 1 Hartree energy 2 2 Many-body system 2 3 Born-Oppenheimer
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More information( 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 informationElectronic Structure of Crystalline Solids
Electronic Structure of Crystalline Solids Computing the electronic structure of electrons in solid materials (insulators, conductors, semiconductors, superconductors) is in general a very difficult problem
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 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 informationCHEM6085: Density Functional Theory Lecture 10
CHEM6085: Density Functional Theory Lecture 10 1) Spin-polarised calculations 2) Geometry optimisation C.-K. Skylaris 1 Unpaired electrons So far we have developed Kohn-Sham DFT for the case of paired
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 informationAn Introduction to Basis Sets:
An Introduction to Basis Sets: The "Holy Grail" of computational chemistry is the calculation of the molecular orbitals (MOs) for a given molecule. IF we can calculate the MOs for a molecule, THEN we can
More informationSame 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 informationAn Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method
Fakultät für Mathematik und Naturwissenschaften - Lehrstuhl für Physikalische Chemie I / Theoretische Chemie An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method Jan-Ole Joswig
More information2 Electronic structure theory
Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are
More informationAN 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 informationarxiv:cond-mat/ v1 10 May 1996
Cohesive energies of cubic III-V semiconductors Beate Paulus, Peter Fulde Max-Planck-Institut für Physik komplexer Systeme, Bayreuther Str. 40, 01187 Dresden, Germany arxiv:cond-mat/9605064v1 10 May 1996
More informationLecture 3: methods and terminology, part I
Eritis sicut Deus, scientes bonum et malum. Lecture 3: methods and terminology, part I Molecular dynamics, internal coordinates, Z-matrices, force fields, periodic boundary conditions, Hartree-Fock theory,
More information5.61 Physical Chemistry Exam III 11/29/12. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry Physical Chemistry.
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry - 5.61 Physical Chemistry Exam III (1) PRINT your name on the cover page. (2) It is suggested that you READ THE ENTIRE EXAM before
More informationElectronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch
Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation There are 10 20 atoms in 1 mm 3
More informationTight-Binding Model of Electronic Structures
Tight-Binding Model of Electronic Structures Consider a collection of N atoms. The electronic structure of this system refers to its electronic wave function and the description of how it is related to
More informationHandbook 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 informationBasis 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 informationPotential energy, from Coulomb's law. Potential is spherically symmetric. Therefore, solutions must have form
Lecture 6 Page 1 Atoms L6.P1 Review of hydrogen atom Heavy proton (put at the origin), charge e and much lighter electron, charge -e. Potential energy, from Coulomb's law Potential is spherically symmetric.
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationPseudopotentials for hybrid density functionals and SCAN
Pseudopotentials for hybrid density functionals and SCAN Jing Yang, Liang Z. Tan, Julian Gebhardt, and Andrew M. Rappe Department of Chemistry University of Pennsylvania Why do we need pseudopotentials?
More informationThe Electronic Structure of Atoms
The Electronic Structure of Atoms Classical Hydrogen-like atoms: Atomic Scale: 10-10 m or 1 Å + - Proton mass : Electron mass 1836 : 1 Problems with classical interpretation: - Should not be stable (electron
More informationELECTRONIC STRUCTURE OF MAGNESIUM OXIDE
Int. J. Chem. Sci.: 8(3), 2010, 1749-1756 ELECTRONIC STRUCTURE OF MAGNESIUM OXIDE P. N. PIYUSH and KANCHAN LATA * Department of Chemistry, B. N. M. V. College, Sahugarh, MADHIPUR (Bihar) INDIA ABSTRACT
More informationInstitut Néel Institut Laue Langevin. Introduction to electronic structure calculations
Institut Néel Institut Laue Langevin Introduction to electronic structure calculations 1 Institut Néel - 25 rue des Martyrs - Grenoble - France 2 Institut Laue Langevin - 71 avenue des Martyrs - Grenoble
More informationJun Yang, Michael Dolg Institute of Theoretical Chemistry, University of Cologne Greinstr. 4, D-50939, Cologne, Germany
Ab initio Density Functional Theory Investigation of Monoclinic BiB 3 O 6 within a LCAO Scheme: the Bi lone pair and its role in electronic structure and vibrational modes Jun Yang, Michael Dolg Institute
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 informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for CrystEngComm. This journal is The Royal Society of Chemistry 2014 Electronic Supplementary Information Configurational and energetical study of the (100) and
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 informationIntermission: Let s review the essentials of the Helium Atom
PHYS3022 Applied Quantum Mechanics Problem Set 4 Due Date: 6 March 2018 (Tuesday) T+2 = 8 March 2018 All problem sets should be handed in not later than 5pm on the due date. Drop your assignments in the
More informationMolecular 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 informationI. CSFs Are Used to Express the Full N-Electron Wavefunction
Chapter 11 One Must be Able to Evaluate the Matrix Elements Among Properly Symmetry Adapted N- Electron Configuration Functions for Any Operator, the Electronic Hamiltonian in Particular. The Slater-Condon
More informationAtomic orbitals of finite range as basis sets. Javier Junquera
Atomic orbitals of finite range as basis sets Javier Junquera Most important reference followed in this lecture in previous chapters: the many body problem reduced to a problem of independent particles
More informationIntroduction to numerical projects
Introduction to numerical projects Here follows a brief recipe and recommendation on how to write a report for each project. Give a short description of the nature of the problem and the eventual numerical
More informationElectron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals
Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals Richard M. Martin UIUC Lecture at Summer School Hands-on introduction to Electronic Structure Materials Computation Center University
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 information4 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 informationChapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set
Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless
More informationECEN 5005 Crystals, Nanocrystals and Device Applications Class 20 Group Theory For Crystals
ECEN 5005 Crystals, Nanocrystals and Device Applications Class 20 Group Theory For Crystals Laporte Selection Rule Polarization Dependence Spin Selection Rule 1 Laporte Selection Rule We first apply this
More informationA 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 informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 5 The Hartree-Fock method C.-K. Skylaris Learning outcomes Be able to use the variational principle in quantum calculations Be able to construct Fock operators
More informationBasis sets for SIESTA. Emilio Artacho. Nanogune, Ikerbasque & DIPC, San Sebastian, Spain Cavendish Laboratory, University of Cambridge
Basis sets for SIESTA Emilio Artacho Nanogune, Ikerbasque & DIPC, San Sebastian, Spain Cavendish Laboratory, University of Cambridge Solving: Basis set Expand in terms of a finite set of basis functions
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationBeyond 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 informationElectron 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 informationTable of Contents. Table of Contents Spin-orbit splitting of semiconductor band structures
Table of Contents Table of Contents Spin-orbit splitting of semiconductor band structures Relavistic effects in Kohn-Sham DFT Silicon band splitting with ATK-DFT LSDA initial guess for the ground state
More informationPseudopotential. Meaning and role
Pseudopotential. Meaning and role Jean-Pierre Flament jean-pierre.flament@univ-lille.fr Laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM) Université de Lille-Sciences et technologies MSSC2018
More informationLecture 9. Hartree Fock Method and Koopman s Theorem
Lecture 9 Hartree Fock Method and Koopman s Theorem Ψ(N) is approximated as a single slater determinant Φ of N orthogonal One electron spin-orbitals. One electron orbital φ i = φ i (r) χ i (σ) χ i (σ)
More informationUsing 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 informationThe Plane-Wave Pseudopotential Method
Hands-on Workshop on Density Functional Theory and Beyond: Computational Materials Science for Real Materials Trieste, August 6-15, 2013 The Plane-Wave Pseudopotential Method Ralph Gebauer ICTP, Trieste
More informationPBS: FROM SOLIDS TO CLUSTERS
PBS: FROM SOLIDS TO CLUSTERS E. HOFFMANN AND P. ENTEL Theoretische Tieftemperaturphysik Gerhard-Mercator-Universität Duisburg, Lotharstraße 1 47048 Duisburg, Germany Semiconducting nanocrystallites like
More informationEnergy levels of group 10 transition metal atoms and ions
189 Appendix B Energy levels of group 10 transition metal atoms and ions B.1 Abstract The energies of the group 10 transition metals (Ni, Pd, and Pt) in different configurations (d 8 s, d 9 s 1, and d
More information7/29/2014. Electronic Structure. Electrons in Momentum Space. Electron Density Matrices FKF FKF. Ulrich Wedig
Electron Density Matrices Density matrices Γ, an alternative to the wavefunction Ψ, for the description of a quantum system Electronic Structure The N-particle density matrix Electrons in Momentum Space
More informationChapter 1: Chemical Bonding
Chapter 1: Chemical Bonding Linus Pauling (1901 1994) January 30, 2017 Contents 1 The development of Bands and their filling 3 2 Different Types of Bonds 7 2.1 Covalent Bonding.............................
More informationChem 344 Final Exam Tuesday, Dec. 11, 2007, 3-?? PM
Chem 344 Final Exam Tuesday, Dec. 11, 2007, 3-?? PM Closed book exam, only pencils and calculators permitted. You may bring and use one 8 1/2 x 11" paper with anything on it. No Computers. Put all of your
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 informationQuantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is,
Chapter : Molecules Quantum mechanics can be used to calculate any property of a molecule The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is, E = Ψ H Ψ Ψ Ψ 1) At first this seems like
More informationonly two orbitals, and therefore only two combinations to worry about, but things get
131 Lecture 1 It is fairly easy to write down an antisymmetric wavefunction for helium since there are only two orbitals, and therefore only two combinations to worry about, but things get complicated
More information1. Hydrogen atom in a box
1. Hydrogen atom in a box Recall H atom problem, V(r) = -1/r e r exact answer solved by expanding in Gaussian basis set, had to solve secular matrix involving matrix elements of basis functions place atom
More informationPractical Advice for Quantum Chemistry Computations. C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology
Practical Advice for Quantum Chemistry Computations C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Choice of Basis Set STO-3G is too small 6-31G* or 6-31G** 6 probably
More information(1/2) M α 2 α, ˆTe = i. 1 r i r j, ˆV NN = α>β
Chemistry 26 Spectroscopy Week # The Born-Oppenheimer Approximation, H + 2. Born-Oppenheimer approximation As for atoms, all information about a molecule is contained in the wave function Ψ, which is the
More informationBasis sets for electron correlation
Basis sets for electron correlation Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry
More informationIntroduction to first-principles modelling and CASTEP
to first-principles modelling and Phil Hasnip to + Atomistic Simulations If we know what the bonding in a material is beforehand, then we can often find good expressions for the forces between atoms, e.g.
More informationMany electrons: Density functional theory Part II. Bedřich Velický VI.
Many electrons: Density functional theory Part II. Bedřich Velický velicky@karlov.mff.cuni.cz VI. NEVF 514 Surface Physics Winter Term 013-014 Troja 1 st November 013 This class is the second devoted to
More informationComputational 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 informationONE AND MANY ELECTRON ATOMS Chapter 15
See Week 8 lecture notes. This is exactly the same as the Hamiltonian for nonrigid rotation. In Week 8 lecture notes it was shown that this is the operator for Lˆ 2, the square of the angular momentum.
More informationComputational 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 informationElectronic 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 informationATOMIC STRUCTURE. Atomic Structure. Atomic orbitals and their energies (a) Hydrogenic radial wavefunctions
ATOMIC STRUCTURE Atomic orbitals and their energies (a) Hydrogenic radial wavefunctions Bundet Boekfa Chem Div, Fac Lib Arts & Sci Kasetsart University Kamphaeng Saen Campus 1 2 Atomic orbitals and their
More informationInstructor background for the discussion points of Section 2
Supplementary Information for: Orbitals Some fiction and some facts Jochen Autschbach Department of Chemistry State University of New York at Buffalo Buffalo, NY 14260 3000, USA Instructor background for
More informationDFT with Hybrid Functionals
DFT with Hybrid Functionals Sanliang Ling University College London 4th CP2K Tutorial, 31st August 4th September 2015, Zurich What are hybrid functionals? Hybrid functionals: mixing non-local Hartree-Fock
More informationElectronic structure of correlated electron systems. Lecture 2
Electronic structure of correlated electron systems Lecture 2 Band Structure approach vs atomic Band structure Delocalized Bloch states Fill up states with electrons starting from the lowest energy No
More informationCondensed matter theory Lecture notes and problem sets 2012/2013
Condensed matter theory Lecture notes and problem sets 2012/2013 Dmitri Ivanov Recommended books and lecture notes: [AM] N. W. Ashcroft and N. D. Mermin, Solid State Physics. [Mar] M. P. Marder, Condensed
More informationGustavus 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 informationHECToR CSE technical meeting, Oxford Parallel Algorithms for the Materials Modelling code CRYSTAL
HECToR CSE technical meeting, Oxford 2009 Parallel Algorithms for the Materials Modelling code CRYSTAL Dr Stanko Tomi Computational Science & Engineering Department, STFC Daresbury Laboratory, UK Acknowledgements
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationCrystal structure prediction
Crystal structure prediction Klaus Doll Institute for Mathematical Physics, TU Braunschweig Max Planck Institute for Solid State Research, Stuttgart MSSC2011, Turin, September 2011 Motivation structure
More informationThe Gutzwiller Density Functional Theory
The Gutzwiller Density Functional Theory Jörg Bünemann, BTU Cottbus I) Introduction 1. Model for an H 2 -molecule 2. Transition metals and their compounds II) Gutzwiller variational theory 1. Gutzwiller
More informationComputational Physics. J. M. Thijssen
Computational Physics J. M. Thijssen Delft University of Technology CAMBRIDGE UNIVERSITY PRESS Contents Preface xi 1 Introduction 1 1.1 Physics and computational physics 1 1.2 Classical mechanics and statistical
More information