A Primer to Electronic Structure Computation
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
Density Functional Theory
Density Functional Theory Fundamentals Video V.i Density Functional Theory: New Developments Donald G. Truhlar Department of Chemistry, University of Minnesota Support: AFOSR, NSF, EMSL Why is electronic structure theory important? Most of the information we want to know about chemistry is in the electron density and electronic energy. dipole moment, Born-Oppenheimer charge distribution, approximation 1927 ... potential energy surface molecular geometry barrier heights bond energies spectra How do we calculate the electronic structure? Example: electronic structure of benzene (42 electrons) Erwin Schrödinger 1925 — wave function theory All the information is contained in the wave function, an antisymmetric function of 126 coordinates and 42 electronic spin components. Theoretical Musings! ● Ψ is complicated. ● Difficult to interpret. ● Can we simplify things? 1/2 ● Ψ has strange units: (prob. density) , ● Can we not use a physical observable? ● What particular physical observable is useful? ● Physical observable that allows us to construct the Hamiltonian a priori. How do we calculate the electronic structure? Example: electronic structure of benzene (42 electrons) Erwin Schrödinger 1925 — wave function theory All the information is contained in the wave function, an antisymmetric function of 126 coordinates and 42 electronic spin components. Pierre Hohenberg and Walter Kohn 1964 — density functional theory All the information is contained in the density, a simple function of 3 coordinates. Electronic structure (continued) Erwin Schrödinger -
Chapter 18 the Single Slater Determinant Wavefunction (Properly Spin and Symmetry Adapted) Is the Starting Point of the Most Common Mean Field Potential
Chapter 18 The single Slater determinant wavefunction (properly spin and symmetry adapted) is the starting point of the most common mean field potential. It is also the origin of the molecular orbital concept. I. Optimization of the Energy for a Multiconfiguration Wavefunction A. The Energy Expression The most straightforward way to introduce the concept of optimal molecular orbitals is to consider a trial wavefunction of the form which was introduced earlier in Chapter 9.II. The expectation value of the Hamiltonian for a wavefunction of the multiconfigurational form Y = S I CIF I , where F I is a space- and spin-adapted CSF which consists of determinental wavefunctions |f I1f I2f I3...fIN| , can be written as: E =S I,J = 1, M CICJ < F I | H | F J > . The spin- and space-symmetry of the F I determine the symmetry of the state Y whose energy is to be optimized. In this form, it is clear that E is a quadratic function of the CI amplitudes CJ ; it is a quartic functional of the spin-orbitals because the Slater-Condon rules express each < F I | H | F J > CI matrix element in terms of one- and two-electron integrals < f i | f | f j > and < f if j | g | f kf l > over these spin-orbitals. B. Application of the Variational Method The variational method can be used to optimize the above expectation value expression for the electronic energy (i.e., to make the functional stationary) as a function of the CI coefficients CJ and the LCAO-MO coefficients {Cn,i} that characterize the spin- orbitals. -
Arxiv:1902.07690V2 [Cond-Mat.Str-El] 12 Apr 2019
Symmetry Projected Jastrow Mean-field Wavefunction in Variational Monte Carlo Ankit Mahajan1, ∗ and Sandeep Sharma1, y 1Department of Chemistry and Biochemistry, University of Colorado Boulder, Boulder, CO 80302, USA We extend our low-scaling variational Monte Carlo (VMC) algorithm to optimize the symmetry projected Jastrow mean field (SJMF) wavefunctions. These wavefunctions consist of a symmetry- projected product of a Jastrow and a general broken-symmetry mean field reference. Examples include Jastrow antisymmetrized geminal power (JAGP), Jastrow-Pfaffians, and resonating valence bond (RVB) states among others, all of which can be treated with our algorithm. We will demon- strate using benchmark systems including the nitrogen molecule, a chain of hydrogen atoms, and 2-D Hubbard model that a significant amount of correlation can be obtained by optimizing the energy of the SJMF wavefunction. This can be achieved at a relatively small cost when multiple symmetries including spin, particle number, and complex conjugation are simultaneously broken and projected. We also show that reduced density matrices can be calculated using the optimized wavefunctions, which allows us to calculate other observables such as correlation functions and will enable us to embed the VMC algorithm in a complete active space self-consistent field (CASSCF) calculation. 1. INTRODUCTION Recently, we have demonstrated that this cost discrep- ancy can be removed by introducing three algorithmic 16 Variational Monte Carlo (VMC) is one of the most ver- improvements . The most significant of these allowed us satile methods available for obtaining the wavefunction to efficiently screen the two-electron integrals which are and energy of a system.1{7 Compared to deterministic obtained by projecting the ab-initio Hamiltonian onto variational methods, VMC allows much greater flexibil- a finite orbital basis. -
“Band Structure” for Electronic Structure Description
www.fhi-berlin.mpg.de Fundamentals of heterogeneous catalysis The very basics of “band structure” for electronic structure description R. Schlögl www.fhi-berlin.mpg.de 1 Disclaimer www.fhi-berlin.mpg.de • This lecture is a qualitative introduction into the concept of electronic structure descriptions different from “chemical bonds”. • No mathematics, unfortunately not rigorous. • For real understanding read texts on solid state physics. • Hellwege, Marfunin, Ibach+Lüth • The intention is to give an impression about concepts of electronic structure theory of solids. 2 Electronic structure www.fhi-berlin.mpg.de • For chemists well described by Lewis formalism. • “chemical bonds”. • Theorists derive chemical bonding from interaction of atoms with all their electrons. • Quantum chemistry with fundamentally different approaches: – DFT (density functional theory), see introduction in this series – Wavefunction-based (many variants). • All solve the Schrödinger equation. • Arrive at description of energetics and geometry of atomic arrangements; “molecules” or unit cells. 3 Catalysis and energy structure www.fhi-berlin.mpg.de • The surface of a solid is traditionally not described by electronic structure theory: • No periodicity with the bulk: cluster or slab models. • Bulk structure infinite periodically and defect-free. • Surface electronic structure theory independent field of research with similar basics. • Never assume to explain reactivity with bulk electronic structure: accuracy, model assumptions, termination issue. • But: electronic -
Simple Molecular Orbital Theory Chapter 5
Simple Molecular Orbital Theory Chapter 5 Wednesday, October 7, 2015 Using Symmetry: Molecular Orbitals One approach to understanding the electronic structure of molecules is called Molecular Orbital Theory. • MO theory assumes that the valence electrons of the atoms within a molecule become the valence electrons of the entire molecule. • Molecular orbitals are constructed by taking linear combinations of the valence orbitals of atoms within the molecule. For example, consider H2: 1s + 1s + • Symmetry will allow us to treat more complex molecules by helping us to determine which AOs combine to make MOs LCAO MO Theory MO Math for Diatomic Molecules 1 2 A ------ B Each MO may be written as an LCAO: cc11 2 2 Since the probability density is given by the square of the wavefunction: probability of finding the ditto atom B overlap term, important electron close to atom A between the atoms MO Math for Diatomic Molecules 1 1 S The individual AOs are normalized: 100% probability of finding electron somewhere for each free atom MO Math for Homonuclear Diatomic Molecules For two identical AOs on identical atoms, the electrons are equally shared, so: 22 cc11 2 2 cc12 In other words: cc12 So we have two MOs from the two AOs: c ,1() 1 2 c ,1() 1 2 2 2 After normalization (setting d 1 and _ d 1 ): 1 1 () () [2(1 S )]1/2 12 [2(1 S )]1/2 12 where S is the overlap integral: 01 S LCAO MO Energy Diagram for H2 H2 molecule: two 1s atomic orbitals combine to make one bonding and one antibonding molecular orbital. -
Iia. Systems of Many Electrons
EE5 in 2008 Lecture notes Univ. of Iceland Hannes J´onsson IIa. Systems of many electrons Exchange of particle labels Consider a system of n identical particles. By identical particles we mean that all intrinsic properties of the particles are the same, such as mass, spin, charge, etc. Figure II.1 A pairwise permutation of the labels on identical particles. In classical mechanics we can in principle follow the trajectory of each individual particle. They are, therefore, distinguishable even though they are identical. In quantum mechanics the particles are not distinguishable if they are close enough or if they interact strongly enough. The quantum mechanical wave function must reflect this fact. When the particles are indistinguishable, the act of labeling the particles is an arbi- trary operation without physical significance. Therefore, all observables must be unaffected by interchange of particle labels. The operators are said to be symmetric under interchange of labels. The Hamiltonian, for example, is symmetric, since the intrinsic properties of all the particles are the same. Let ψ(1, 2,...,n) be a solution to the Schr¨odinger equation, (Here n represents all the coordinates, both spatial and spin, of the particle labeled with n). Let Pij be an operator that permutes (or ‘exchanges’) the labels i and j: Pijψ(1, 2,...,i,...,j,...,n) = ψ(1, 2,...,j,...,i,...,n). This means that the function Pij ψ depends on the coordinates of particle j in the same way that ψ depends on the coordinates of particle i. Since H is symmetric under interchange of labels: H(Pijψ) = Pij Hψ, i.e., [Pij,H]=0. -
Electronic Structure Methods
Electronic Structure Methods • One-electron models – e.g., Huckel theory • Semiempirical methods – e.g., AM1, PM3, MNDO • Single-reference based ab initio methods •Hartree-Fock • Perturbation theory – MP2, MP3, MP4 • Coupled cluster theory – e.g., CCSD(T) • Multi-configurational based ab initio methods • MCSCF and CASSCF (also GVB) • MR-MP2 and CASPT2 • MR-CI •Ab Initio Methods for IPs, EAs, excitation energies •CI singles (CIS) •TDHF •EOM and Greens function methods •CC2 • Density functional theory (DFT)- Combine functionals for exchange and for correlation • Local density approximation (LDA) •Perdew-Wang 91 (PW91) • Becke-Perdew (BP) •BeckeLYP(BLYP) • Becke3LYP (B3LYP) • Time dependent DFT (TDDFT) (for excited states) • Hybrid methods • QM/MM • Solvation models Why so many methods to solve Hψ = Eψ? Electronic Hamiltonian, BO approximation 1 2 Z A 1 Z AZ B H = − ∑∇i − ∑∑∑+ + (in a.u.) 2 riA iB〈〈j rij A RAB 1/rij is what makes it tough (nonseparable)!! Hartree-Fock method: • Wavefunction antisymmetrized product of orbitals Ψ =|φ11 (τφτ ) ⋅⋅⋅NN ( ) | ← Slater determinant accounts for indistinguishability of electrons In general, τ refers to both spatial and spin coordinates For the 2-electron case 1 Ψ = ϕ (τ )ϕ (τ ) = []ϕ (τ )ϕ (τ ) −ϕ (τ )ϕ (τ ) 1 1 2 2 2 1 1 2 2 2 1 1 2 Energy minimized – variational principle 〈Ψ | H | Ψ〉 Vary orbitals to minimize E, E = 〈Ψ | Ψ〉 keeping orbitals orthogonal hi ϕi = εi ϕi → orbitals and orbital energies 2 φ j (r2 ) 1 2 Z J = dr h = − ∇ − A + (J − K ) i ∑ ∫ 2 i i ∑ i i j ≠ i r12 2 riA ϕϕ()rr () Kϕϕ()rdrr= ji22 () ii121∑∫ j ji≠ r12 Characteristics • Each e- “sees” average charge distribution due to other e-. -
Slater Decomposition of Laughlin States
Slater Decomposition of Laughlin States∗ Gerald V. Dunne Department of Physics University of Connecticut 2152 Hillside Road Storrs, CT 06269 USA [email protected] 14 May, 1993 Abstract The second-quantized form of the Laughlin states for the fractional quantum Hall effect is discussed by decomposing the Laughlin wavefunctions into the N-particle Slater basis. A general formula is given for the expansion coefficients in terms of the characters of the symmetric group, and the expansion coefficients are shown to possess numerous interesting symmetries. For expectation values of the density operator it is possible to identify individual dominant Slater states of the correct uniform bulk density and filling fraction in the physically relevant N limit. →∞ arXiv:cond-mat/9306022v1 8 Jun 1993 1 Introduction The variational trial wavefunctions introduced by Laughlin [1] form the basis for the theo- retical understanding of the quantum Hall effect [2, 3]. The Laughlin wavefunctions describe especially stable strongly correlated states, known as incompressible quantum fluids, of a two dimensional electron gas in a strong magnetic field. They also form the foundation of the hierarchy structure of the fractional quantum Hall effect [4, 5, 6]. In the extreme low temperature and strong magnetic field limit, the single particle elec- tron states are restricted to the lowest Landau level. In the absence of interactions, this Landau level has a high degeneracy determined by the magnetic flux through the sample [7]. ∗UCONN-93-4 1 Much (but not all) of the physics of the quantum Hall effect may be understood in terms of the restriction of the dynamics to fixed Landau levels [8, 3, 9]. -
The Electronic Structure of Molecules an Experiment Using Molecular Orbital Theory
The Electronic Structure of Molecules An Experiment Using Molecular Orbital Theory Adapted from S.F. Sontum, S. Walstrum, and K. Jewett Reading: Olmstead and Williams Sections 10.4-10.5 Purpose: Obtain molecular orbital results for total electronic energies, dipole moments, and bond orders for HCl, H 2, NaH, O2, NO, and O 3. The computationally derived bond orders are correlated with the qualitative molecular orbital diagram and the bond order obtained using the NO bond force constant. Introduction to Spartan/Q-Chem The molecule, above, is an anti-HIV protease inhibitor that was developed at Abbot Laboratories, Inc. Modeling of molecular structures has been instrumental in the design of new drug molecules like these inhibitors. The molecular model kits you used last week are very useful for visualizing chemical structures, but they do not give the quantitative information needed to design and build new molecules. Q-Chem is another "modeling kit" that we use to augment our information about molecules, visualize the molecules, and give us quantitative information about the structure of molecules. Q-Chem calculates the molecular orbitals and energies for molecules. If the calculation is set to "Equilibrium Geometry," Q-Chem will also vary the bond lengths and angles to find the lowest possible energy for your molecule. Q- Chem is a text-only program. A “front-end” program is necessary to set-up the input files for Q-Chem and to visualize the results. We will use the Spartan program as a graphical “front-end” for Q-Chem. In other words Spartan doesn’t do the calculations, but rather acts as a go-between you and the computational package. -
PART 1 Identical Particles, Fermions and Bosons. Pauli Exclusion
O. P. Sushkov School of Physics, The University of New South Wales, Sydney, NSW 2052, Australia (Dated: March 1, 2016) PART 1 Identical particles, fermions and bosons. Pauli exclusion principle. Slater determinant. Variational method. He atom. Multielectron atoms, effective potential. Exchange interaction. 2 Identical particles and quantum statistics. Consider two identical particles. 2 electrons 2 protons 2 12C nuclei 2 protons .... The wave function of the pair reads ψ = ψ(~r1, ~s1,~t1...; ~r2, ~s2,~t2...) where r1,s1, t1, ... are variables of the 1st particle. r2,s2, t2, ... are variables of the 2nd particle. r - spatial coordinate s - spin t - isospin . - other internal quantum numbers, if any. Omit isospin and other internal quantum numbers for now. The particles are identical hence the quantum state is not changed under the permutation : ψ (r2,s2; r1,s1) = Rψ(r1,s1; r2,s2) , where R is a coefficient. Double permutation returns the wave function back, R2 = 1, hence R = 1. ± The spin-statistics theorem claims: * Particles with integer spin have R =1, they are called bosons (Bose - Einstein statistics). * Particles with half-integer spins have R = 1, they are called fermions (Fermi − - Dirac statistics). 3 The spin-statistics theorem can be proven in relativistic quantum mechanics. Technically the theorem is based on the fact that due to the structure of the Lorentz transformation the wave equation for a particle with spin 1/2 is of the first order in time derivative (see discussion of the Dirac equation later in the course). At the same time the wave equation for a particle with integer spin is of the second order in time derivative. -
Electronic Structure Calculations and Density Functional Theory
Electronic Structure Calculations and Density Functional Theory Rodolphe Vuilleumier P^olede chimie th´eorique D´epartement de chimie de l'ENS CNRS { Ecole normale sup´erieure{ UPMC Formation ModPhyChem { Lyon, 09/10/2017 Rodolphe Vuilleumier (ENS { UPMC) Electronic Structure and DFT ModPhyChem 1 / 56 Many applications of electronic structure calculations Geometries and energies, equilibrium constants Reaction mechanisms, reaction rates... Vibrational spectroscopies Properties, NMR spectra,... Excited states Force fields From material science to biology through astrochemistry, organic chemistry... ... Rodolphe Vuilleumier (ENS { UPMC) Electronic Structure and DFT ModPhyChem 2 / 56 Length/,("-%0,11,("2,"3,4)(",3"2,"1'56+,+&" and time scales 0,1 nm 10 nm 1µm 1ps 1ns 1µs 1s Variables noyaux et atomes et molécules macroscopiques Modèles Gros-Grains électrons (champs) ! #$%&'(%')$*+," #-('(%')$*+," #.%&'(%')$*+," !" Rodolphe Vuilleumier (ENS { UPMC) Electronic Structure and DFT ModPhyChem 3 / 56 Potential energy surface for the nuclei c The LibreTexts libraries Rodolphe Vuilleumier (ENS { UPMC) Electronic Structure and DFT ModPhyChem 4 / 56 Born-Oppenheimer approximation Total Hamiltonian of the system atoms+electrons: H^T = T^N + V^NN + H^ MI me approximate system wavefunction: ! ~ ~ ~ ΨS (RI ; ~r1;:::; ~rN ) = Ψ(RI )Ψ0(~r1;:::; ~rN ; RI ) ~ Ψ0(~r1;:::; ~rN ; RI ): the ground state electronic wavefunction of the ~ electronic Hamiltonian H^ at fixed ionic configuration RI n o Energy of the electrons at this ionic configuration: ~ E0(RI ) = Ψ0 -
The Atomic and Electronic Structure of Dislocations in Ga Based Nitride Semiconductors Imad Belabbas, Pierre Ruterana, Jun Chen, Gérard Nouet
The atomic and electronic structure of dislocations in Ga based nitride semiconductors Imad Belabbas, Pierre Ruterana, Jun Chen, Gérard Nouet To cite this version: Imad Belabbas, Pierre Ruterana, Jun Chen, Gérard Nouet. The atomic and electronic structure of dislocations in Ga based nitride semiconductors. Philosophical Magazine, Taylor & Francis, 2006, 86 (15), pp.2245-2273. 10.1080/14786430600651996. hal-00513682 HAL Id: hal-00513682 https://hal.archives-ouvertes.fr/hal-00513682 Submitted on 1 Sep 2010 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Philosophical Magazine & Philosophical Magazine Letters For Peer Review Only The atomic and electronic structure of dislocations in Ga based nitride semiconductors Journal: Philosophical Magazine & Philosophical Magazine Letters Manuscript ID: TPHM-05-Sep-0420.R2 Journal Selection: Philosophical Magazine Date Submitted by the 08-Dec-2005 Author: Complete List of Authors: BELABBAS, Imad; ENSICAEN, SIFCOM Ruterana, Pierre; ENSICAEN, SIFCOM CHEN, Jun; IUT, LRPMN NOUET, Gérard; ENSICAEN, SIFCOM Keywords: GaN, atomic structure, dislocations, electronic structure Keywords (user supplied): http://mc.manuscriptcentral.com/pm-pml Page 1 of 87 Philosophical Magazine & Philosophical Magazine Letters 1 2 3 4 The atomic and electronic structure of dislocations in Ga 5 6 based nitride semiconductors 7 8 9 10 11 12 1,3 1 2 1 13 I.