Lecture on First-principles Computations (14): The Linear Combination of Atomic Orbitals (LCAO) Method
|
|
- Merryl Young
- 5 years ago
- Views:
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 任新国 (Xinguo Ren) 中国科学技术大学量子信息重点实验室 Key Laboratory of Quantum Information, USTC Hefei, 2018.10.24 Recall: orthogonalized
More informationLecture 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 informationLecture 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 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 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 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 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 informationPatterns, 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 informationKey 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 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 informationSet 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 informationThe 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 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 informationKey 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 informationPseudopotentials 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 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 informationComparison 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 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 informationVariational 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 informationIntroduction 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 informationPseudopotentials: 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 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 informationThe 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 informationChemistry 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 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 informationDensity 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 informationThe 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 informationNorm-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 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 informationChemistry 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 informationNumerical 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 informationIntroduction 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 informationThe 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 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 informationIs 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 informationIntroduction 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.
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 information2.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 informationTest 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 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 informationAtomic 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 informationAddition 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 informationDept 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 informationChemistry 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 informationProjector 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 informationQuantum 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 informationSolved 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 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 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 informationPhysics 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 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 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 informationIntroduction 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 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 informationBasics 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 information5.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 informationPseudopotential 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 informationChemistry 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 informationQuantum 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 informationIntro 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 informationH!!!! = 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 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 information8.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 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 informationApproximation 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 informationThe 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 informationOrigin 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 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 informationASE: 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 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 informationChemistry 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 informationPHY492: 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 informationElectronic 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 informationLecture #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 informationComparison 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 information6: 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 informationWelcome 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 informationThe 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 informationWe 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 information7.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 informationECE 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 informationVibrational 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 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 informationeff (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 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 informationCHEM-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 informationKey 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 informationI. 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 informationChemistry 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 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 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 informationLecture 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 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 informationLecture 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 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 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 informationPhysics 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 informationChem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.
Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative
More informationWorkshop 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 informationatoms 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