Lecture on First-principles Computations (14): The Linear Combination of Atomic Orbitals (LCAO) Method

Size: px
Start display at page:

Download "Lecture on First-principles Computations (14): The Linear Combination of Atomic Orbitals (LCAO) Method"

Transcription

1 Lecture on First-principles Computations (14): The Linear Combination of Atomic Orbitals (LCAO) Method 任新国 (Xinguo Ren) 中国科学技术大学量子信息重点实验室 Key Laboratory of Quantum Information, USTC Hefei,

2 Recall: the plane wave basis set V r = V r + R I Periodic crystal potential 2 2m + V r ψ nk r = ε nk r Ritz variational principle: expanding ψ nk r in terms of basis functions: ψ nk r = c μ,n k χ μk r μ H μν k c υn k = ε nk c υn k v Basis-type I -- Plane waves : χ μk r = χ G k r = 1 Ω μ G ei k+g r

3 Linear combination of atomic orbitals (LCAO) V r = V r + R I Periodic crystal potential 2 2m + V r ψ nk r = ε nk r Expanding ψ nk r in terms of basis functions: ψ nk r = c μ,n μ k χ μk r Basis-type II localized atom-centered orbitals: The LCAO method χ μk r = 1 N R χ a,j,l,m r R r a e ik R μ = a, j, l, m Atomic position within the unit cell

4 General features of local orbitals χ a,j,l,m r r a = R a,j,l r r a Y lm r r a Centered at atomic positions. The angular part Y lm r are spherical harmonics, while the radial functions R a,j,l r are flexible -- they can be analytical functions (Gaussians, Slater-type orbitals), or completely numerical. j denotes different radial functions with the same angular momentum, and hence controls the size of the basis set (from tight-binding to full quantitative calculations). Single-ζ ( SZ ): 1s1p for O, C, N, etc. (minimal basis) Double-ζ ( DZ ): 2s2p for O, C, N, etc. Double- ζ plus polarization function ( DZP ): 2s2p1d for O, C, N, etc. Compact in size (Efficient), and suitable for local symmetry analysis Constructing systematically converging atomic basis set is a highly nontrivial task, and needs special attention.

5 The Hamiltonian and overlap matrices Hψ nk r = ε nk ψ nk r ψ nk r = c μ,n k χ μk r μ H μν k c ν,n k = ε nk S μν k c ν,n k ν ν H μν k = χ μk H χ νk S μν k = χ μk χ νk

6 The Hamiltonian and overlap matrices in the LCAO method H μν k c ν,n k = ε nk S μν k c ν,n k ν ν H μν k = χ μk H χ νk, S μν k = χ μk χ νk Now χ μk r = χ a,j,l,m r R r a e ik R R H μν k = 1 N R H μν R e ik R, S μν k = 1 N R S μν R e ik R H μν R = χ a,j,l,m r r a S μν R = χ a,j,l,m r r a H χ a,j,l,m r r a R χ a,j,l,m r r a R

7 Two-center integrals: Multi-center integrals S μν R = χ a,j,l,m r r a T μν R = χ a,j,l,m r r a χ a,j,l,m r r a R 2 2m χ a,j,l,m r r a R Can be computed efficiently in Fourier space. Three-center integrals: V μν R = χ a,j,l,m r r a V χ a,j,l,m r r a R V r = a,r V a r r a R Three centers For numerical orbitals, this is often evaluated by grid integration.

8 Slater-type orbitals (STO) χ STO abc x, y, z = Nx a y b z c e ζr Normalization factor a, b, c control the angular momentum: l = a + b + c. ζcontrols the spread of the orbital. These basis functions don't have radial nodes, and not pure spherical harmonics. For 1s, it is hydrogen-like orbital. Long-range and short-range behaviors are correct. Hydrogen 1s orbital: ψ 1s r = 1 π ( 1 a 0 ) 3/2 e r/a 0 a 0 : Bohr radius

9 Gaussian-type orbitals (GTO) χ GTO lmn x, y, z = Nx l y m z n e ζr2 Normalization factor N = 2l 1!! 2m 1!! 2n 1!! π3/2 (4α) l+m+n (2α) 3/2 l, m, n control the angular momentum. ζ controls the width of the orbital. They are not hydrogen-like orbitals for 1s. Long-range and short-range behaviors are incorrect. The Coulomb integrals involving GTOs can be easily computed

10 Comparison between STOs and GTOs φ STO r = e αr i a i e β ir 2

11 The computation integrals among GTOs Gaussian product theorem: χ GTO l A = N A x 1 m A y 1 n A z 1 A e ζ 1 r R A 2, χ GTO l B = N B x 2 m B y 2 n B z 2 B e ζ 2 r R B 2 r A = r R A = (x A, y A, z A ) e ζ 1 r R A 2 e ζ 2 r R B 2 = e ζ 1ζ 2 ζ 1 +ζ 2 R A R B 2 e ζ 1+ζ 2 r ζ 1R A +ζ 2 R B ζ 1 +ζ 2 χ A GTO χ B GTO = e ζ 1ζ 2 ζ 1 +ζ 2 R A R B 2 linear combination of Gaussians centered at ζ 1R A + ζ 2 R B ζ 1 + ζ 2 Multi-centered integrals can be reduced to linear combination of twocentered integrals, that can be computed analytically.

12 STOs, GTOs, and NAOs GTOs STOs NAOs Accuracy of the basis functions Number of basis functions at a given level of accuracy Number of integrals to be evaluated Difficulty of evaluating the integrals

13 Symmetry analysis for two-center integrals χ a,j,l,m χ a,j,l,m = d 3 rr a,j,l r r a Y lm r r a R a,j,l r r a R Y l m r r a R = S μ,ν δ m.m S ll m δ m.m l = s, p, d, f, m = σ, π, δ, D = R + r a r a z ssσ spσ ppσ ppπ pdπ

14 From tight-binding to full calculations J. M. Soler et al, J. Phys: Condens. Matter 14, 2745 (2002)

15 Spanish Initiative for Electronic Simulations with Thousands of Atoms

16 FHI-aims (

17 ABACUS (

18 Other computer codes based on NAOs OpenMX (Open source package Material explorer) Dmol (Bernard Delley, Swiss, all-electron) FPLO (Full-potential local orbital, all-electron, Dresden) ADF (Amsterdam, Slater-type orbitals ) Crystal (Gaussian orbitals for solids)

19 Quantum chemistry codes based on GTOs Gaussian (commercial code) ( NWChem (free code) ( Q-CHEM (commercial code) ( Turbomole (

20 -value project in materials science: reproducibility of density-functional theory calculations for solids. i a, b = 1.06V 0,i (Eb,i V E a,i (V)) 2 dv 0.12V 0,i 0.94V 0,i A set of 71 elemental solids. K. Lejaeghere, S. Cottenier, et al. Science 351, 1415 (2016). 71 Test accuracies of methods/codes (basis sets, pseudopotentials, relativisitic treatment etc. )

21 Comparing solid-state DFT codes, basis sets, and potentials Code Version Basis Electron treatment Δ-value WIEN2k 13.1 LAPW/APW+lo all-electron 0 mev/atom FHI-aims tier2 numerical orbitals all-electron (relativistic atomic_zora scalar) 0.2 mev/atom Exciting development version LAPW+xlo all-electron 0.2 mev/atom VASP plane waves PAW 2015 GW-ready (5.4) 0.3 mev/atom FHI-aims tier2 numerical orbitals all-electron (relativistic zora scalar 1e-12) 0.3 mev/atom Quantum ESPRESSO 5.1 plane waves SSSP Accuracy (link is external) (mixed NC/US/PAW potential library) 0.3 mev/atom Elk APW+lo all-electron 0.3 mev/atom ABINIT plane waves PAW JTH v1.0 (link is external) 0.4 mev/atom FLEUR 0.26 LAPW (+lo) all-electron 0.4 mev/atom Quantum ESPRESSO 5.1 plane waves SSSP Efficiency (link is external) (mixed NC/US/PAW potential library) 0.4 mev/atom CASTEP 9.0 plane waves OTFG CASTEP mev/atom

Lecture on First-principles Computations (11): The Linearized Augmented Plane-wave Method

Lecture on First-principles Computations (11): The Linearized Augmented Plane-wave Method Lecture on First-principles Computations (11): The Linearized Augmented Plane-wave Method 任新国 (Xinguo Ren) 中国科学技术大学量子信息重点实验室 Key Laboratory of Quantum Information, USTC Hefei, 2018.10.24 Recall: orthogonalized

More information

Lecture on First-Principles Computation (11): The LAPW Method

Lecture on First-Principles Computation (11): The LAPW Method Lecture on First-Principles Computation (11): The LAPW Method 任新国 (Xinguo Ren) 中国科学技术大学量子信息实验室 (Key Laboratory of Quantum Information, USTC) 2015-10-16 Recall: OPW and Pseudopotential Methods Problems:

More information

Lecture on First-principles Computation (2): The Hartree-Fock theory and the homogeneous electron gas

Lecture on First-principles Computation (2): The Hartree-Fock theory and the homogeneous electron gas Lecture on First-principles Computation (2): The Hartree-Fock theory and the homogeneous electron gas 任新国 (Xinguo Ren) 中国科学技术大学量子信息重点实验室 Key Laboratory of Quantum Information, USTC Hefei, 2016.9.18 Recall:

More information

Tight-Binding Model of Electronic Structures

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

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

Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals

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

Atomic orbitals of finite range as basis sets. Javier Junquera

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

Patterns, Correlations, and Causality in Big Data of Materials: Analytics for Novel Materials Discovery

Patterns, Correlations, and Causality in Big Data of Materials: Analytics for Novel Materials Discovery Patterns, Correlations, and Causality in Big Data of Materials: Analytics for Novel Materials Discovery From the periodic table of the elements to a chart (a map) of materials: Organize materials according

More information

Key concepts in Density Functional Theory (II)

Key concepts in Density Functional Theory (II) Kohn-Sham scheme and band structures European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Present Address: LPMCN Université

More information

Atomic orbitals of finite range as basis sets. Javier Junquera

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

Set the initial conditions r i. Update neighborlist. Get new forces F i

Set the initial conditions r i. Update neighborlist. Get new forces F i v Set the initial conditions r i ( t 0 ), v i ( t 0 ) Update neighborlist Quantum mechanical models Get new forces F i ( r i ) Solve the equations of motion numerically over time step Δt : r i ( t n )

More information

The Linearized Augmented Planewave (LAPW) Method

The Linearized Augmented Planewave (LAPW) Method The Linearized Augmented Planewave (LAPW) Method David J. Singh Oak Ridge National Laboratory E T [ ]=T s [ ]+E ei [ ]+E H [ ]+E xc [ ]+E ii {T s +V ks [,r]} I (r)= i i (r) Need tools that are reliable

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

Key concepts in Density Functional Theory (II) Silvana Botti

Key concepts in Density Functional Theory (II) Silvana Botti Kohn-Sham scheme, band structure and optical spectra European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address:

More information

Pseudopotentials and Basis Sets. How to generate and test them

Pseudopotentials and Basis Sets. How to generate and test them Pseudopotentials and Basis Sets How to generate and test them Pseudopotential idea Atomic Si Core electrons highly localized very depth energy are chemically inert 1s 2 2s 2 2p 6 3s 2 3p 2 Valence wave

More information

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

Comparison of various abinitio codes used in periodic calculations

Comparison of various abinitio codes used in periodic calculations Comparison of various abinitio codes used in periodic calculations 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India & Center for Materials Science and Nanotechnology,

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

Variational Methods for Electronic Structure

Variational Methods for Electronic Structure Variational Methods for Electronic Structure The hydrogen atom is a two-body system consisting of a proton and an electron. If spin and relativistic effects are ignored, then the Schrödinger equation for

More information

Introduction to DFT and its Application to Defects in Semiconductors

Introduction to DFT and its Application to Defects in Semiconductors Introduction to DFT and its Application to Defects in Semiconductors Noa Marom Physics and Engineering Physics Tulane University New Orleans The Future: Computer-Aided Materials Design Can access the space

More information

Pseudopotentials: design, testing, typical errors

Pseudopotentials: design, testing, typical errors Pseudopotentials: design, testing, typical errors Kevin F. Garrity Part 3 National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015 Testing Basis Sets

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

The Linearized Augmented Planewave (LAPW) Method (WIEN2k, ELK, FLEUR)

The Linearized Augmented Planewave (LAPW) Method (WIEN2k, ELK, FLEUR) The Linearized Augmented Planewave (LAPW) Method (WIEN2k, ELK, FLEUR) David J. Singh Oak Ridge National Laboratory E T [ρ]=t s [ρ]+e ei [ρ]+e H [ρ]+e xc [ρ]+e ii {T s +V ks [ρ,r]}ϕ I (r)=ε i ϕ i (r) Please

More information

Chemistry 120A 2nd Midterm. 1. (36 pts) For this question, recall the energy levels of the Hydrogenic Hamiltonian (1-electron):

Chemistry 120A 2nd Midterm. 1. (36 pts) For this question, recall the energy levels of the Hydrogenic Hamiltonian (1-electron): April 6th, 24 Chemistry 2A 2nd Midterm. (36 pts) For this question, recall the energy levels of the Hydrogenic Hamiltonian (-electron): E n = m e Z 2 e 4 /2 2 n 2 = E Z 2 /n 2, n =, 2, 3,... where Ze is

More information

A One-Slide Summary of Quantum Mechanics

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

Density Functional Theory Machinery

Density Functional Theory Machinery Solid State Theory Physics 545 Density Functional Theory- Density Functional Theory Machinery Calculating the Wave Function DFT (and other methods) iterate to self-consistency Guess the wave functions

More information

The electronic structure of materials 2 - DFT

The electronic structure of materials 2 - DFT Quantum mechanics 2 - Lecture 9 December 19, 2012 1 Density functional theory (DFT) 2 Literature Contents 1 Density functional theory (DFT) 2 Literature Historical background The beginnings: L. de Broglie

More information

Norm-conserving pseudopotentials and basis sets in electronic structure calculations. Javier Junquera. Universidad de Cantabria

Norm-conserving pseudopotentials and basis sets in electronic structure calculations. Javier Junquera. Universidad de Cantabria Norm-conserving pseudopotentials and basis sets in electronic structure calculations Javier Junquera Universidad de Cantabria Outline Pseudopotentials Why pseudopotential approach is useful Orthogonalized

More information

Potential energy, from Coulomb's law. Potential is spherically symmetric. Therefore, solutions must have form

Potential 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 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

Numerical quality control in computational materials databases

Numerical quality control in computational materials databases Numerical quality control in computational materials databases Björn Bieniek 1, Mansur Sahid 1, Christian Carbogno 1, Luca Ghiringhelli 1, and Matthias Scheffler 1 1 Fritz-Haber-Institut der MPG - Berlin

More information

Introduction to First-Principles Method

Introduction to First-Principles Method Joint ICTP/CAS/IAEA School & Workshop on Plasma-Materials Interaction in Fusion Devices, July 18-22, 2016, Hefei Introduction to First-Principles Method by Guang-Hong LU ( 吕广宏 ) Beihang University Computer

More information

The Projector Augmented Wave method

The Projector Augmented Wave method The Projector Augmented Wave method Advantages of PAW. The theory. Approximations. Convergence. 1 The PAW method is... What is PAW? A technique for doing DFT calculations efficiently and accurately. An

More information

ONE AND MANY ELECTRON ATOMS Chapter 15

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

Is Band better than ADF? The core of the issue

Is Band better than ADF? The core of the issue Is Band better than ADF? The core of the issue Introduction Introduction What is Band? Introduction What is Band? Slightly different basis set for ADF/Band Introduction What is Band? Slightly different

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

ψ s a ˆn a s b ˆn b ψ Hint: Because the state is spherically symmetric the answer can depend only on the angle between the two directions.

ψ s a ˆn a s b ˆn b ψ Hint: Because the state is spherically symmetric the answer can depend only on the angle between the two directions. 1. Quantum Mechanics (Fall 2004) Two spin-half particles are in a state with total spin zero. Let ˆn a and ˆn b be unit vectors in two arbitrary directions. Calculate the expectation value of the product

More information

2.4. Quantum Mechanical description of hydrogen atom

2.4. Quantum Mechanical description of hydrogen atom 2.4. Quantum Mechanical description of hydrogen atom Atomic units Quantity Atomic unit SI Conversion Ang. mom. h [J s] h = 1, 05459 10 34 Js Mass m e [kg] m e = 9, 1094 10 31 kg Charge e [C] e = 1, 6022

More information

Test Set for Materials Science and Engineering

Test Set for Materials Science and Engineering Test Set for Materials Science and Engineering Igor Ying Zhang Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195, Berlin, Germany 2015-07-20 FHI-aims workshop Definition of Test set

More information

Lecture 9. Hartree Fock Method and Koopman s Theorem

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

Atomic Structure and Atomic Spectra

Atomic Structure and Atomic Spectra Atomic Structure and Atomic Spectra Atomic Structure: Hydrogenic Atom Reading: Atkins, Ch. 10 (7 판 Ch. 13) The principles of quantum mechanics internal structure of atoms 1. Hydrogenic atom: one electron

More information

Addition of Angular Momenta

Addition of Angular Momenta Addition of Angular Momenta What we have so far considered to be an exact solution for the many electron problem, should really be called exact non-relativistic solution. A relativistic treatment is needed

More information

Dept of Mechanical Engineering MIT Nanoengineering group

Dept of Mechanical Engineering MIT Nanoengineering group 1 Dept of Mechanical Engineering MIT Nanoengineering group » Recap of HK theorems and KS equations» The physical meaning of the XC energy» Solution of a one-particle Schroedinger equation» Pseudo Potentials»

More information

Chemistry 532 Practice Final Exam Fall 2012 Solutions

Chemistry 532 Practice Final Exam Fall 2012 Solutions Chemistry 53 Practice Final Exam Fall Solutions x e ax dx π a 3/ ; π sin 3 xdx 4 3 π cos nx dx π; sin θ cos θ + K x n e ax dx n! a n+ ; r r r r ˆL h r ˆL z h i φ ˆL x i hsin φ + cot θ cos φ θ φ ) ˆLy i

More information

Projector augmented wave Implementation

Projector augmented wave Implementation Projector augmented wave Implementation Peter. E. Blöchl Institute for Theoretical Physics Clausthal University of Technology, Germany http://www.pt.tu-clausthal.de/atp/ 1 = Projector augmented wave +

More information

Quantum Theory of Angular Momentum and Atomic Structure

Quantum Theory of Angular Momentum and Atomic Structure Quantum Theory of Angular Momentum and Atomic Structure VBS/MRC Angular Momentum 0 Motivation...the questions Whence the periodic table? Concepts in Materials Science I VBS/MRC Angular Momentum 1 Motivation...the

More information

Solved radial equation: Last time For two simple cases: infinite and finite spherical wells Spherical analogs of 1D wells We introduced auxiliary func

Solved radial equation: Last time For two simple cases: infinite and finite spherical wells Spherical analogs of 1D wells We introduced auxiliary func Quantum Mechanics and Atomic Physics Lecture 16: The Coulomb Potential http://www.physics.rutgers.edu/ugrad/361 h / d/361 Prof. Sean Oh Solved radial equation: Last time For two simple cases: infinite

More information

Quantum Mechanical Simulations

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

Physics 228 Today: Ch 41: 1-3: 3D quantum mechanics, hydrogen atom

Physics 228 Today: Ch 41: 1-3: 3D quantum mechanics, hydrogen atom Physics 228 Today: Ch 41: 1-3: 3D quantum mechanics, hydrogen atom Website: Sakai 01:750:228 or www.physics.rutgers.edu/ugrad/228 Happy April Fools Day Example / Worked Problems What is the ratio of the

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

MODELING MATTER AT NANOSCALES

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

Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 27, December 5, 2014

Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 27, December 5, 2014 Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi Lecture 27, December 5, 2014 (Some material in this lecture has been adapted from Cramer, C. J.

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

Basics of APW-like methods. Lars Nordström Uppsala University

Basics of APW-like methods. Lars Nordström Uppsala University Basics of APW-like methods Lars Nordström Uppsala University APW-like methods All electron (no pseudo-potential) Full-potential Robust Accurate the benchmark to test against But, not too time consuming...

More information

5.111 Lecture Summary #6

5.111 Lecture Summary #6 5.111 Lecture Summary #6 Readings for today: Section 1.9 (1.8 in 3 rd ed) Atomic Orbitals. Read for Lecture #7: Section 1.10 (1.9 in 3 rd ed) Electron Spin, Section 1.11 (1.10 in 3 rd ed) The Electronic

More information

Pseudopotential generation and test by the ld1.x atomic code: an introduction

Pseudopotential generation and test by the ld1.x atomic code: an introduction and test by the ld1.x atomic code: an introduction SISSA and DEMOCRITOS Trieste (Italy) Outline 1 2 3 Spherical symmetry - I The Kohn and Sham (KS) equation is (in atomic units): [ 1 ] 2 2 + V ext (r)

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

Quantum Mechanics for Scientists and Engineers

Quantum Mechanics for Scientists and Engineers Quantum Mechanics for Scientists and Engineers Syllabus and Textbook references All the main lessons (e.g., 1.1) and units (e.g., 1.1.1) for this class are listed below. Mostly, there are three lessons

More information

Intro to ab initio methods

Intro to ab initio methods Lecture 2 Part A Intro to ab initio methods Recommended reading: Leach, Chapters 2 & 3 for QM methods For more QM methods: Essentials of Computational Chemistry by C.J. Cramer, Wiley (2002) 1 ab initio

More information

H!!!! = E! Lecture 7 - Atomic Structure. Chem 103, Section F0F Unit II - Quantum Theory and Atomic Structure Lecture 7. Lecture 7 - Introduction

H!!!! = E! Lecture 7 - Atomic Structure. Chem 103, Section F0F Unit II - Quantum Theory and Atomic Structure Lecture 7. Lecture 7 - Introduction Chem 103, Section F0F Unit II - Quantum Theory and Atomic Structure Lecture 7 Lecture 7 - Atomic Structure Reading in Silberberg - Chapter 7, Section 4 The Qunatum-Mechanical Model of the Atom The Quantum

More information

An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method

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

8.1 The hydrogen atom solutions

8.1 The hydrogen atom solutions 8.1 The hydrogen atom solutions Slides: Video 8.1.1 Separating for the radial equation Text reference: Quantum Mechanics for Scientists and Engineers Section 10.4 (up to Solution of the hydrogen radial

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

Approximation Methods in QM

Approximation Methods in QM Chapter 3 Approximation Methods in QM Contents 3.1 Time independent PT (nondegenerate)............... 5 3. Degenerate perturbation theory (PT)................. 59 3.3 Time dependent PT and Fermi s golden

More information

The Hydrogen Atom Chapter 20

The Hydrogen Atom Chapter 20 4/4/17 Quantum mechanical treatment of the H atom: Model; The Hydrogen Atom Chapter 1 r -1 Electron moving aroundpositively charged nucleus in a Coulombic field from the nucleus. Potential energy term

More information

Origin of the first Hund rule in He-like atoms and 2-electron quantum dots

Origin of the first Hund rule in He-like atoms and 2-electron quantum dots in He-like atoms and 2-electron quantum dots T Sako 1, A Ichimura 2, J Paldus 3 and GHF Diercksen 4 1 Nihon University, College of Science and Technology, Funabashi, JAPAN 2 Institute of Space and Astronautical

More information

Lecture 3: methods and terminology, part I

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

ASE: A programmable environment for calculations with many electronic structure codes

ASE: A programmable environment for calculations with many electronic structure codes ASE: A programmable environment for calculations with many electronic structure codes Ask Hjorth Larsen Center for Atomic-scale Materials Design, Technical University of Denmark 16th ETSF workshop September

More information

Introduction to numerical projects

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

Chemistry 881 Lecture Topics Fall 2001

Chemistry 881 Lecture Topics Fall 2001 Chemistry 881 Lecture Topics Fall 2001 Texts PHYSICAL CHEMISTRY A Molecular Approach McQuarrie and Simon MATHEMATICS for PHYSICAL CHEMISTRY, Mortimer i. Mathematics Review (M, Chapters 1,2,3 & 4; M&S,

More information

PHY492: Nuclear & Particle Physics. Lecture 4 Nature of the nuclear force. Reminder: Investigate

PHY492: Nuclear & Particle Physics. Lecture 4 Nature of the nuclear force. Reminder: Investigate PHY49: Nuclear & Particle Physics Lecture 4 Nature of the nuclear force Reminder: Investigate www.nndc.bnl.gov Topics to be covered size and shape mass and binding energy charge distribution angular momentum

More information

Electronic structure of solids: basic concepts and methods

Electronic structure of solids: basic concepts and methods Electronic structure of solids: basic concepts and methods Ondřej Šipr II. NEVF 514 Surface Physics Winter Term 2016-2017 Troja, 21st October 2016 Outline A bit of formal mathematics for the beginning

More information

Lecture #1. Review. Postulates of quantum mechanics (1-3) Postulate 1

Lecture #1. Review. Postulates of quantum mechanics (1-3) Postulate 1 L1.P1 Lecture #1 Review Postulates of quantum mechanics (1-3) Postulate 1 The state of a system at any instant of time may be represented by a wave function which is continuous and differentiable. Specifically,

More information

Comparison of Different Methods and Codes: (L)APW, LMTO, PAW, Pseudo Potentials, Gaussians, etc. Jörg Behler

Comparison of Different Methods and Codes: (L)APW, LMTO, PAW, Pseudo Potentials, Gaussians, etc. Jörg Behler Comparison of Different Methods and Codes: (L)APW, LMTO, PAW, Pseudo Potentials, Gaussians, etc. Jörg Behler Fritz-Haber-Institut der Max-Planck-Gesellschaft Faradayweg 4-6, D-14195 Berlin, Germany Density-Functional

More information

6: Plane waves, unit cells, k- points and all that

6: Plane waves, unit cells, k- points and all that The Nuts and Bolts of First-Principles Simulation 6: Plane waves, unit cells, k- points and all that Durham, 6th- 13th December 2001 CASTEP Developers Group with support from the ESF ψ k Network Overview

More information

Welcome back to PHY 3305

Welcome back to PHY 3305 Welcome back to PHY 3305 Today s Lecture: Hydrogen Atom Part I John von Neumann 1903-1957 One-Dimensional Atom To analyze the hydrogen atom, we must solve the Schrodinger equation for the Coulomb potential

More information

The 3 dimensional Schrödinger Equation

The 3 dimensional Schrödinger Equation Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum

More information

We now turn to our first quantum mechanical problems that represent real, as

We now turn to our first quantum mechanical problems that represent real, as 84 Lectures 16-17 We now turn to our first quantum mechanical problems that represent real, as opposed to idealized, systems. These problems are the structures of atoms. We will begin first with hydrogen-like

More information

7.1 Variational Principle

7.1 Variational Principle 7.1 Variational Principle Suppose that you want to determine the ground-state energy E g for a system described by H, but you are unable to solve the time-independent Schrödinger equation. It is possible

More information

ECE 535 Theory of Semiconductors and Semiconductor Devices Fall 2015 Homework # 5 Due Date: 11/17/2015

ECE 535 Theory of Semiconductors and Semiconductor Devices Fall 2015 Homework # 5 Due Date: 11/17/2015 ECE 535 Theory of Semiconductors and Semiconductor Devices Fall 2015 Homework # 5 Due Date: 11/17/2015 Problem # 1 Two carbon atoms and four hydrogen atoms form and ethane molecule with the chemical formula

More information

Vibrational motion. Harmonic oscillator ( 諧諧諧 ) - A particle undergoes harmonic motion. Parabolic ( 拋物線 ) (8.21) d 2 (8.23)

Vibrational motion. Harmonic oscillator ( 諧諧諧 ) - A particle undergoes harmonic motion. Parabolic ( 拋物線 ) (8.21) d 2 (8.23) Vibrational motion Harmonic oscillator ( 諧諧諧 ) - A particle undergoes harmonic motion F == dv where k Parabolic V = 1 f k / dx = is Schrodinge h m d dx ψ f k f x the force constant x r + ( 拋物線 ) 1 equation

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

eff (r) which contains the influence of angular momentum. On the left is

eff (r) which contains the influence of angular momentum. On the left is 1 Fig. 13.1. The radial eigenfunctions R nl (r) of bound states in a square-well potential for three angular-momentum values, l = 0, 1, 2, are shown as continuous lines in the left column. The form V (r)

More information

The Plane-Wave Pseudopotential Method

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

CHEM-UA 127: Advanced General Chemistry I

CHEM-UA 127: Advanced General Chemistry I 1 CHEM-UA 127: Advanced General Chemistry I Notes for Lecture 11 Nowthatwehaveintroducedthebasicconceptsofquantummechanics, wecanstarttoapplythese conceptsto build up matter, starting from its most elementary

More information

Key concepts in Density Functional Theory

Key concepts in Density Functional Theory From the many body problem to the Kohn-Sham scheme ILM (LPMCN) CNRS, Université Lyon 1 - France European Theoretical Spectroscopy Facility (ETSF) December 12, 2012 Lyon Outline 1 The many-body problem

More information

I. RADIAL PROBABILITY DISTRIBUTIONS (RPD) FOR S-ORBITALS

I. RADIAL PROBABILITY DISTRIBUTIONS (RPD) FOR S-ORBITALS 5. Lecture Summary #7 Readings for today: Section.0 (.9 in rd ed) Electron Spin, Section. (.0 in rd ed) The Electronic Structure of Hydrogen. Read for Lecture #8: Section. (. in rd ed) Orbital Energies

More information

Chemistry 483 Lecture Topics Fall 2009

Chemistry 483 Lecture Topics Fall 2009 Chemistry 483 Lecture Topics Fall 2009 Text PHYSICAL CHEMISTRY A Molecular Approach McQuarrie and Simon A. Background (M&S,Chapter 1) Blackbody Radiation Photoelectric effect DeBroglie Wavelength Atomic

More information

Computational Physics. J. M. Thijssen

Computational 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

Introduction to DFTB. Marcus Elstner. July 28, 2006

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

Lecture 10. Central potential

Lecture 10. Central potential Lecture 10 Central potential 89 90 LECTURE 10. CENTRAL POTENTIAL 10.1 Introduction We are now ready to study a generic class of three-dimensional physical systems. They are the systems that have a central

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

Lecture 3. Solving the Non-Relativistic Schroedinger Equation for a spherically symmetric potential

Lecture 3. Solving the Non-Relativistic Schroedinger Equation for a spherically symmetric potential Lecture 3 Last lecture we were in the middle of deriving the energies of the bound states of the Λ in the nucleus. We will continue with solving the non-relativistic Schroedinger equation for a spherically

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

The Electronic Structure of Atoms

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

Physics 606, Quantum Mechanics, Final Exam NAME ( ) ( ) + V ( x). ( ) and p( t) be the corresponding operators in ( ) and x( t) : ( ) / dt =...

Physics 606, Quantum Mechanics, Final Exam NAME ( ) ( ) + V ( x). ( ) and p( t) be the corresponding operators in ( ) and x( t) : ( ) / dt =... Physics 606, Quantum Mechanics, Final Exam NAME Please show all your work. (You are graded on your work, with partial credit where it is deserved.) All problems are, of course, nonrelativistic. 1. Consider

More information

Chem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.

Chem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components. Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative

More information

Workshop on: ATOMIC STRUCTURE AND TRANSITIONS: THEORY IN USING SUPERSTRUCTURE PROGRAM

Workshop on: ATOMIC STRUCTURE AND TRANSITIONS: THEORY IN USING SUPERSTRUCTURE PROGRAM Workshop on: ATOMIC STRUCTURE AND TRANSITIONS: THEORY IN USING SUPERSTRUCTURE PROGRAM PROF. SULTANA N. NAHAR Astronomy, Ohio State U, Columbus, Ohio, USA Email: nahar.1@osu.edu http://www.astronomy.ohio-state.edu/

More information

atoms and light. Chapter Goal: To understand the structure and properties of atoms.

atoms and light. Chapter Goal: To understand the structure and properties of atoms. Quantum mechanics provides us with an understanding of atomic structure and atomic properties. Lasers are one of the most important applications of the quantummechanical properties of atoms and light.

More information