Introduction to numerical projects
|
|
- Julianna Ryan
- 5 years ago
- Views:
Transcription
1 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 methods you have used. Describe the algorithm you have used and/or developed. Here you may find it convenient to use pseudocoding. In many cases you can describe the algorithm in the program itself. Include the source code of your program. Comment your program properly. If possible, try to find analytic solutions, or known limits in order to test your program when developing the code. Include your results either in figure form or in a table. Remember to label your results. All tables and figures should have relevant captions and labels on the axes. Try to evaluate the reliabilty and numerical stability/precision of your results. If possible, include a qualitative and/or quantitative discussion of the numerical stability, eventual loss of precision etc. Try to give an interpretation of you results in your answers to the problems. Critique: if possible include your comments and reflections about the exercise, whether you felt you learnt something, ideas for improvements and other thoughts you ve made when solving the exercise. We wish to keep this course at the interactive level and your comments can help us improve it. We do appreciate your comments. Try to establish a practice where you log your work at the computerlab. You may find such a logbook very handy at later stages in your work, especially when you don t properly remember what a previous test version of your program did. Here you could also record the time spent on solving the exercise, various algorithms you may have tested or other topics which you feel worthy of mentioning. Format for electronic delivery of report and programs The preferred format for the report is a PDF file. As programming language we prefer that you choose between C/C++ and Fortran. You could also use Java or Python as programming languages. Matlab/Maple/Mathematica/IDL are not accepted, but you can use them to check your results where possible. Finally, we do prefer that you work together. Optimal working groups consist of 2-3 students, but more people can collaborate. You can then hand in a common report. Below you will also find a guide on how to write up the project as a scientific article. We recommend this option when writing the final report. 1
2 Project for FYS4411/9411 Spring 2014, Advanced Hartree- Fock Part 1: Hartree-Fock calculations of closed-shell atoms The Hartree-Fock functional is written as E[Φ] = + µ=1 µ h µ 1 2 [ µ=1 ν=1 µν 1 r ij µν µν 1 r ij νµ With the given functional, we can perform at least two types of variational strategies. Vary the Slater determinant by changing the spatial part of the single-particle wave functions themselves. Expand the single-particle functions in a known basis and vary the coefficients, that is, the new function single-particle wave function a is written as a linear expansion in terms of a fixed basis φ (harmonic oscillator, Laguerre polynomials etc) ]. ψ a = λ C aλ φ λ, Both cases lead to a new Slater determinant which is related to the previous via a unitary transformation. The second one is the one we will use in this project. 1a) Consider a Slater determinant built up of single-particle orbitals ψ λ, with λ = 1, 2,..., N. The unitary transformation ψ a = C aλ φ λ, λ brings us into the new basis. Show that the new basis is orthonormal. Show that the new Slater determinant constructed from the new single-particle wave functions can be written as the determinant based on the previous basis and the determinant of the matrix C. Show that the old and the new Slater determinants are equal up to a complex constant with absolute value unity. (Hint, C is a unitary matrix). 1b) Minimizing with respect to Ckα, remembering that Ckα and C kα are independent and defining h HF αγ = α h γ + CaβC aδ αβ V γδ AS, a=1 show that you can write the Hartree-Fock equations as h HF αγ C kγ = ɛ k C kα. Explain the meaning of the different terms. γ 2 βδ
3 1c) We will now set up the Hartree-Fock program by using hydrogen-like single-particle functions for the 1s, 2s and 3s orbitals. We will apply these basis functions only to the helium and beryllium atoms. For all other systems, we will use so-called Gaussian type of orbitals (GTO). These will be introduced after we have set up our first Hartree-Fock programs. For Hydrogen like-orbitals, the matrix elements involving the Coulomb interaction, can all be evaluated in a closed form, see the table below for all integrals involving these orbitals. Your first step is to write a function which reads in these integrals and sets up the antisymmetrized matrix elements. (Hint: the table lists only the integrals, there is no spin). 1d) With the above ingredients we are now ready to solve the Hartree-Fock equations for the helium and beryllium atoms. Set up the Hartree-Fock equations for the ground states of helium and beryllium with the electrons occupying the respective hydrogen-like orbitals 1s, 2s and 3s. There is no spin-orbit part in the two-body Hamiltonian. Find the single-particle energies iteratively as explained in the slides and lecture notes. Find also the total binding ground state energy and compare this to the energy obtaining without performing the Hartree-Fock calculations. The experimental values are atomic units (a.u.) and a.u. for the helium and beryllium atoms, respectively. In setting up the equations make sure your Hartree-Fock matrix is block diagonal in the spin quantum numbers. Since there is no spin-dependent part in the Hamiltonian, states with the same l and n will be degenerate in energy with respect to the spin projections. Part 2: Hartree-Fock calculations of atoms and molecules using GTOs 2a) Our next steps involves an object orientation of the code and the inclusion of Gaussiantype of orbitals (GTOs) as detailed in the slides. Your task here is to repeat the calculations of the helium and beryllium atoms using GTOs using three and five primitive functions, a so-called STO-3G and STO-5L basis. Compare your results with those obtained with the hydrogen-like basis discussed above. Comment your results. 2b) Use thereafter a chosen GTO basis and compute the ground state energy of the neon and argon atoms and compare with existing theoretical calculations in the literature. Comment your results. 2c) We switch now to the computation of molecules and will use the hydrogen and beryllium molecules as our first test case with GTOs. The H 2 molecule consists of two protons and two electrons with a ground state energy E = ev or
4 a.u. and the equilibrium distance between the two hydrogen atoms of r 0 = 1.40 Bohr radii (recall that a Bohr radius is m. The Be 2 molecule contains eight electrons and eight protons (in addition we have eight neutrons but their charge is zero). Useful benchmarks for these molecules are the articles of Moskowitz and Kalos, Int. Journal of Quantum Chemistry XX, 1107 (1981) for He and H 2 and Filippi, Singh and Umrigar, J. Chemical Physics 105, 123 (1996) for the Be 2 molecule. You may also find Jørgen Høgberget s Master of Science thesis at the University of Oslo, 2013 a useful reference. 2d) Our final aim is to compute the ground state energy and equilibrium configuration for the water molecule H 2 O and Si 2 O using GTOs. Your task here is to figure out which is the optimal basis and perform calculations of the ground state energy and compare with the literature. With the optimal orbitals we can in turn also compute the charge distribution of electrons, which is given by the the so-called density distribution. Multiplying with the electrical charge gives the charge distribution. The density distribution is defined as ρ(r) = ρ(r 1 ) = dr 2... dr N Ψ(r 1, r 2... r N ) 2, which in a Hartree-Fock based approach becomes ρ(r) = ψ i (r) 2, i=1 with ψ the single-particle wave functions (our best basis). A useful check that the numerics works properly is to integrate the density distributions since it has to give us the total number of electrons in the system, that is N = drρ(r). Compute the density distribution for the H 2, Be 2, H 2 O and Si 2 O molecules at their equilbrium positions (the value of R where the ground state energy is at its minimum). Make plots of the density distributions. Give an interpretation of your results and discuss the physics behind the results. How to write a scientific report What should it contain? A possible structure If you wish to write up your findings as a scientific article, here are some simple guidelines An introduction where you explain the rational for the physics case and what you have done. At the end of the introduction you should give a brief summary of the structure of the report 4
5 Theoretical models and technicalities. This is the methods section. Results and discussion Conclusions and perspectives Appendix with extra material Bibliography You don t need to answer all questions in a chronological order. When you write the introduction you could focus on the following aspects A central aim is to study the role of correlations due to the repulsion between the electrons. To do this we have singled out several systems from atomic and molecular physics. We test also the wave functions by computing onebody densities. The methods section could contain the following elements: Describe the methods (quantum mechanical and algorithms) You need to explain Hartree-Fock theory. The trial wave functions. Why do you choose the functions you do? You can also plug in some calculations to demonstrate your code, such as selected runs from exercises 1a-1d for helium and beryllium. The results section should include results which support your conclusions. You don t need to present all results, focus on those which you find the most relevant. Provide also benchmarks which you have used in order to validate your code. Present your discussions of results and compare with other calculations and theories and give critical assessments of what you have done. In the concluding section you could focus on State your main findings and interpretations Try as far as possible to present perspectives for future work Try to discuss the pros and cons of the method and possible improvements Additional material, code listings etc can be placed in various appendices. Finally, you should provide a proper list of relevant references. Give always references to material you base your work on, either scientific articles/reports or books. When citing articles, refer as: name(s) of author(s), journal, volume (boldfaced), page and year. For books, refer as: name(s) of author(s), title of book, publisher, place and year, eventual page numbers. 5
6 Appendix Table 1 contains the matrix elements for the radial integrals to be used for the direct part and the exchange part. Note again that these integrals do not include spin. You will need these integrals in solving the Hartree-Fock equations for the helium and beryllium atoms when using hydrogen-like single-particle states. 6
7 Table 1: Closed form expressions for the Coulomb matrix elements. The nomenclature is 1 = 1s, 2 = 2s and 3 = 3s, with no spin degrees of freedom. 11 V 11 = (5Z)/8 11 V 12 = (4096 2Z)/ V 13 = (1269 3Z)/ V 21 = (4096 2Z)/ V 22 = (16Z)/ V 23 = ( Z)/ V 31 = (1269 3Z)/ V 32 = ( Z)/ V 33 = (189Z)/ V 11 = (4096 2Z)/ V 12 = (17Z)/81 12 V 13 = ( Z)/ V 21 = (16Z)/ V 22 = (512 2Z)/ V 23 = (2160 3Z)/ V 31 = ( Z)/ V 32 = ( Z)/ V 33 = ( Z)/ V 11 = (1269 3Z)/ V 12 = ( Z)/ V 13 = (815Z)/ V 21 = ( Z)/ V 22 = (2160 3Z)/ V 23 = ( Z)/ V 31 = (189Z)/ V 32 = ( Z)/ V 33 = (617Z)/( ) 21 V 11 = (4096 2Z)/ V 12 = (16Z)/ V 13 = ( Z)/ V 21 = (17Z)/81 21 V 22 = (512 2Z)/ V 23 = ( Z)/ V 31 = ( Z)/ V 32 = (2160 3Z)/ V 33 = ( Z)/ V 11 = (16Z)/ V 12 = (512 2Z)/ V 13 = (2160 3Z)/ V 21 = (512 2Z)/ V 22 = (77Z)/ V 23 = ( Z)/ V 31 = (2160 3Z)/ V 32 = ( Z)/ V 33 = (73008Z)/ V 11 = ( Z)/ V 12 = (2160 3Z)/ V 13 = ( Z)/ V 21 = ( Z)/ V 22 = ( Z)/ V 23 = (32857Z)/ V 31 = ( Z)/ V 32 = (73008Z)/ V 33 = ( /3Z)/ V 11 = (1269 3Z)/ V 12 = ( Z)/ V 13 = (189Z)/ V 21 = ( Z)/ V 22 = (2160 3Z)/ V 23 = ( Z)/ V 31 = (815Z)/ V 32 = ( Z)/ V 33 = (617Z)/( ) 32 V 11 = ( Z)/ V 12 = ( Z)/ V 13 = ( Z)/ V 21 = (2160 3Z)/ V 22 = ( Z)/ V 23 = (73008Z)/ V 31 = ( Z)/ V 32 = (32857Z)/ V 33 = ( /3Z)/ V 11 = (189Z)/ V 12 = ( Z)/ V 13 = (617Z)/( ) 33 V 21 = ( Z)/ V 22 = (73008Z)/ V 23 = ( /3Z)/ V 31 = (617Z)/( ) 33 V 7 32 = ( /3Z)/ V 33 = (17Z)/256
CHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 5 The Hartree-Fock method C.-K. Skylaris Learning outcomes Be able to use the variational principle in quantum calculations Be able to construct Fock operators
More information2 Electronic structure theory
Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are
More informationHartree-Fock-Roothan Self-Consistent Field Method
Hartree-Fock-Roothan Self-Consistent Field Method 1. Helium Here is a summary of the derivation of the Hartree-Fock equations presented in class. First consider the ground state of He and start with with
More informationNext topic: Quantum Field Theories for Quantum Many-Particle Systems; or "Second Quantization"
Next topic: Quantum Field Theories for Quantum Many-Particle Systems; or "Second Quantization" Outline 1 Bosons and Fermions 2 N-particle wave functions ("first quantization" 3 The method of quantized
More information7.1 Creation and annihilation operators
Chapter 7 Second Quantization Creation and annihilation operators. Occupation number. Anticommutation relations. Normal product. Wick s theorem. One-body operator in second quantization. Hartree- Fock
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 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 informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry Chemicum 4th floor vesa.hanninen@helsinki.fi September 10, 2013 Lecture 3. Electron correlation methods September
More informationSelf-consistent Field
Chapter 6 Self-consistent Field A way to solve a system of many electrons is to consider each electron under the electrostatic field generated by all other electrons. The many-body problem is thus reduced
More informationIntermission: Let s review the essentials of the Helium Atom
PHYS3022 Applied Quantum Mechanics Problem Set 4 Due Date: 6 March 2018 (Tuesday) T+2 = 8 March 2018 All problem sets should be handed in not later than 5pm on the due date. Drop your assignments in the
More informationExercises FYS-KJM4480
Exercises FYS-KJM4480 Exercise 1 Consider the Slater determinant Φ AS λ (x 1x 2...x N ;α 1 α 2...α N ) = 1 N ( ) p P ψ αi (x i ). N! where P is an operator which permutes the coordinates of two particles.
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 informationChem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 21, November 12, 2014
Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 204 Laura Gagliardi Lecture 2, November 2, 204 (Some material in this lecture has been adapted from Cramer, C. J. Essentials
More informationonly two orbitals, and therefore only two combinations to worry about, but things get
131 Lecture 1 It is fairly easy to write down an antisymmetric wavefunction for helium since there are only two orbitals, and therefore only two combinations to worry about, but things get complicated
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 4 Molecular orbitals C.-K. Skylaris Learning outcomes Be able to manipulate expressions involving spin orbitals and molecular orbitals Be able to write down
More informationHartree, Hartree-Fock and post-hf methods
Hartree, Hartree-Fock and post-hf methods MSE697 fall 2015 Nicolas Onofrio School of Materials Engineering DLR 428 Purdue University nonofrio@purdue.edu 1 The curse of dimensionality Let s consider a multi
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More 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 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 informationr R A 1 r R B + 1 ψ(r) = αψ A (r)+βψ B (r) (5) where we assume that ψ A and ψ B are ground states: ψ A (r) = π 1/2 e r R A ψ B (r) = π 1/2 e r R B.
Molecules Initial questions: What are the new aspects of molecules compared to atoms? What part of the electromagnetic spectrum can we probe? What can we learn from molecular spectra? How large a molecule
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 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 informationPreliminary Quantum Questions
Preliminary Quantum Questions Thomas Ouldridge October 01 1. Certain quantities that appear in the theory of hydrogen have wider application in atomic physics: the Bohr radius a 0, the Rydberg constant
More informationBrief review of Quantum Mechanics (QM)
Brief review of Quantum Mechanics (QM) Note: This is a collection of several formulae and facts that we will use throughout the course. It is by no means a complete discussion of QM, nor will I attempt
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 informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 20, March 8, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 20, March 8, 2006 Solved Homework We determined that the two coefficients in our two-gaussian
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 informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry room B430, Chemicum 4th floor vesa.hanninen@helsinki.fi September 3, 2013 Introduction and theoretical backround September
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 information( ( ; R H = 109,677 cm -1
CHAPTER 9 Atomic Structure and Spectra I. The Hydrogenic Atoms (one electron species). H, He +1, Li 2+, A. Clues from Line Spectra. Reminder: fundamental equations of spectroscopy: ε Photon = hν relation
More informationA Computer Study of Molecular Electronic Structure
A Computer Study of Molecular Electronic Structure The following exercises are designed to give you a brief introduction to some of the types of information that are now readily accessible from electronic
More informationHartree-Fock. It can also be viewed as the lowest order term in perturbative expansion (linear in U) with respect to interaction between electrons.
Hartree-Fock It is probably the simplest method to treat the many-particle system. The dynamic many particle problem is replaced by an effective one-electron problem: electron is moving in an effective
More informationChemistry 883 Computational Quantum Chemistry
Chemistry 883 Computational Quantum Chemistry Instructor Contact Information Professor Benjamin G. Levine levine@chemistry.msu.edu 215 Chemistry Building 517-353-1113 Office Hours Tuesday 9:00-11:00 am
More informationChem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 28, December 08, 2014
Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi Lecture 28, December 08, 2014 Solved Homework Water, H 2 O, involves 2 hydrogen atoms and an oxygen
More informationExercises FYS-KJM4480
Exercises FYS-KJM4480 Exercise 1 Consider the Slater determinant Φ AS λ (x 1x 2...x N ;α 1 α 2...α N ) = 1 N ( ) p P ψ αi (x i ). N! where P is an operator which permutes the coordinates of two particles.
More informationAnswers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014
Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde
More information3: Many electrons. Orbital symmetries. l =2 1. m l
3: Many electrons Orbital symmetries Atomic orbitals are labelled according to the principal quantum number, n, and the orbital angular momentum quantum number, l. Electrons in a diatomic molecule experience
More informationThe Hartree-Fock approximation
Contents The Born-Oppenheimer approximation Literature Quantum mechanics 2 - Lecture 7 November 21, 2012 Contents The Born-Oppenheimer approximation Literature 1 The Born-Oppenheimer approximation 2 3
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 informationAtomic Systems (PART I)
Atomic Systems (PART I) Lecturer: Location: Recommended Text: Dr. D.J. Miller Room 535, Kelvin Building d.miller@physics.gla.ac.uk Joseph Black C407 (except 15/1/10 which is in Kelvin 312) Physics of Atoms
More informationCHEMISTRY 4021/8021 MIDTERM EXAM 1 SPRING 2014
CHEMISTRY 4021/8021 Q1) Propose a simple, united-atom molecular mechanics force-field needed to generate a potential energy surface for an isolated molecule of acetone (Me 2 CO). I.e., provide an energy
More informationQuantum Physics II (8.05) Fall 2002 Assignment 12 and Study Aid
Quantum Physics II (8.05) Fall 2002 Assignment 12 and Study Aid Announcement This handout includes 9 problems. The first 5 are the problem set due. The last 4 cover material from the final few lectures
More informationSCF calculation on HeH +
SCF calculation on HeH + Markus Meuwly Department of Chemistry, University of Basel, Basel, Switzerland Abstract This document describes the main steps involved in carrying out a SCF calculation on the
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 informationOBJECTIVES: By the end of class, students will be able to. SWBAT find the number of shells and number of valence electrons for the different elements.
7 th Grade Science Name: Unit: Matter and Periodic Table Lesson: Atoms 4 Valance Elections I OBJECTIVES: By the end of class, students will be able to SWBAT find the number of shells and number of valence
More informationAtomic Structure Atoms are very small ~ metres All atoms are made up of three sub-atomic particles: protons, neutrons and electrons
IB Chemistry (unit ) ATOMIC THEORY Atomic Structure - Recap Questions Define the following words: Atom Element Molecule Compound Atomic Structure Atoms are very small ~ - metres All atoms are made up of
More informationBorn-Oppenheimer Approximation
Born-Oppenheimer Approximation Adiabatic Assumption: Nuclei move so much more slowly than electron that the electrons that the electrons are assumed to be obtained if the nuclear kinetic energy is ignored,
More informationP3317 HW from Lecture 7+8 and Recitation 4
P3317 HW from Lecture 7+8 and Recitation 4 Due Friday Tuesday September 25 Problem 1. In class we argued that an ammonia atom in an electric field can be modeled by a two-level system, described by a Schrodinger
More informationChm 331 Fall 2015, Exercise Set 4 NMR Review Problems
Chm 331 Fall 015, Exercise Set 4 NMR Review Problems Mr. Linck Version.0. Compiled December 1, 015 at 11:04:44 4.1 Diagonal Matrix Elements for the nmr H 0 Find the diagonal matrix elements for H 0 (the
More informationLUMO + 1 LUMO. Tómas Arnar Guðmundsson Report 2 Reikniefnafræði G
Q1: Display all the MOs for N2 in your report and classify each one of them as bonding, antibonding or non-bonding, and say whether the symmetry of the orbital is σ or π. Sketch a molecular orbital diagram
More informationQuantum Theory of Many-Particle Systems, Phys. 540
Quantum Theory of Many-Particle Systems, Phys. 540 Questions about organization Second quantization Questions about last class? Comments? Similar strategy N-particles Consider Two-body operators in Fock
More informationECE440 Nanoelectronics. Lecture 07 Atomic Orbitals
ECE44 Nanoelectronics Lecture 7 Atomic Orbitals Atoms and atomic orbitals It is instructive to compare the simple model of a spherically symmetrical potential for r R V ( r) for r R and the simplest hydrogen
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 informationIntroduction to Path Integral Monte Carlo. Part I.
Introduction to Path Integral Monte Carlo. Part I. Alexey Filinov, Jens Böning, Michael Bonitz Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
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 informationChem 467 Supplement to Lecture 19 Hydrogen Atom, Atomic Orbitals
Chem 467 Supplement to Lecture 19 Hydrogen Atom, Atomic Orbitals Pre-Quantum Atomic Structure The existence of atoms and molecules had long been theorized, but never rigorously proven until the late 19
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 information1 Introduction. 2 The hadronic many body problem
Models Lecture 18 1 Introduction In the next series of lectures we discuss various models, in particluar models that are used to describe strong interaction problems. We introduce this by discussing the
More informationHartree-Fock Theory. ˆf(r)χ i (x) = ɛ i χ i (x) (1) Z A K=1 I. DERIVATION OF THE HARTREE-FOCK EQUATIONS. A. The energy of a Slater determinant
Hartree-Fock Theory The HF approximation plays a crucial role in chemistry and constitutes the starting point for more elaborate treatments of electron correlation. Furthermore, many semi-empirical methods
More informationNew Frontiers in Nuclear Structure Theory
New Frontiers in Nuclear Structure Theory From Realistic Interactions to the Nuclear Chart Robert Roth Institut für Kernphysik Technical University Darmstadt Overview Motivation Nucleon-Nucleon Interactions
More informationMean-field concept. (Ref: Isotope Science Facility at Michigan State University, MSUCL-1345, p. 41, Nov. 2006) 1/5/16 Volker Oberacker, Vanderbilt 1
Mean-field concept (Ref: Isotope Science Facility at Michigan State University, MSUCL-1345, p. 41, Nov. 2006) 1/5/16 Volker Oberacker, Vanderbilt 1 Static Hartree-Fock (HF) theory Fundamental puzzle: The
More informationQuantum Mechanics - I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras
Quantum Mechanics - I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 14 Exercises on Quantum Expectation Values (Refer Slide Time: 00:07) In the last couple
More informationYingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2
Purdue University CHM 67300 Computational Quantum Chemistry REVIEW Yingwei Wang October 10, 2013 Review: Prof Slipchenko s class, Fall 2013 Contents 1 Hartree energy 2 2 Many-body system 2 3 Born-Oppenheimer
More informationQuantum Mechanics-I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras. Lecture - 21 Square-Integrable Functions
Quantum Mechanics-I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 21 Square-Integrable Functions (Refer Slide Time: 00:06) (Refer Slide Time: 00:14) We
More informationStudent Exploration: Electron Configuration
Name: Date: Student Exploration: Electron Configuration Vocabulary: atomic number, atomic radius, Aufbau principle, chemical family, diagonal rule, electron configuration, Hund s rule, orbital, Pauli exclusion
More informationLS coupling. 2 2 n + H s o + H h f + H B. (1) 2m
LS coupling 1 The big picture We start from the Hamiltonian of an atomic system: H = [ ] 2 2 n Ze2 1 + 1 e 2 1 + H s o + H h f + H B. (1) 2m n e 4πɛ 0 r n 2 4πɛ 0 r nm n,m Here n runs pver the electrons,
More informationLecture 4: Hartree-Fock Theory
Lecture 4: Hartree-Fock Theory One determinant to rule them all, One determinant to find them, One determinant to bring them all and in the darkness bind them Second quantization rehearsal The formalism
More informationMolecular Physics. Attraction between the ions causes the chemical bond.
Molecular Physics A molecule is a stable configuration of electron(s) and more than one nucleus. Two types of bonds: covalent and ionic (two extremes of same process) Covalent Bond Electron is in a molecular
More informationConsequently, the exact eigenfunctions of the Hamiltonian are also eigenfunctions of the two spin operators
VI. SPIN-ADAPTED CONFIGURATIONS A. Preliminary Considerations We have described the spin of a single electron by the two spin functions α(ω) α and β(ω) β. In this Sect. we will discuss spin in more detail
More informationLecture #13 1. Incorporating a vector potential into the Hamiltonian 2. Spin postulates 3. Description of spin states 4. Identical particles in
Lecture #3. Incorporating a vector potential into the Hamiltonian. Spin postulates 3. Description of spin states 4. Identical particles in classical and QM 5. Exchange degeneracy - the fundamental problem
More informatione L 2m e the Bohr magneton
e L μl = L = μb 2m with : μ B e e 2m e the Bohr magneton Classical interation of magnetic moment and B field: (Young and Freedman, Ch. 27) E = potential energy = μ i B = μbcosθ τ = torque = μ B, perpendicular
More informationOptimization of quantum Monte Carlo (QMC) wave functions by energy minimization
Optimization of quantum Monte Carlo (QMC) wave functions by energy minimization Julien Toulouse, Cyrus Umrigar, Roland Assaraf Cornell Theory Center, Cornell University, Ithaca, New York, USA. Laboratoire
More informationLecture 4. Magic Numbers
Lecture 4 In last lecture we discussed how one can solve the Schroedinger equation to determine the allowed energies of a particle trapped in a spherically symmetric potential. We will apply these methods
More informationP3317 HW from Lecture and Recitation 7
P3317 HW from Lecture 1+13 and Recitation 7 Due Oct 16, 018 Problem 1. Separation of variables Suppose we have two masses that can move in 1D. They are attached by a spring, yielding a Hamiltonian where
More informationwith a proper choice of the potential U(r). Clearly, we should include the potential of the ions, U ion (r):.
The Hartree Equations So far we have ignored the effects of electron-electron (e-e) interactions by working in the independent electron approximation. In this lecture, we shall discuss how this effect
More information4πε. me 1,2,3,... 1 n. H atom 4. in a.u. atomic units. energy: 1 a.u. = ev distance 1 a.u. = Å
H atom 4 E a me =, n=,,3,... 8ε 0 0 π me e e 0 hn ε h = = 0.59Å E = me (4 πε ) 4 e 0 n n in a.u. atomic units E = r = Z n nao Z = e = me = 4πε = 0 energy: a.u. = 7. ev distance a.u. = 0.59 Å General results
More informationLecture 19: Building Atoms and Molecules
Lecture 19: Building Atoms and Molecules +e r n = 3 n = 2 n = 1 +e +e r y even Lecture 19, p 1 Today Nuclear Magnetic Resonance Using RF photons to drive transitions between nuclear spin orientations in
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 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 informationQuantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is,
Chapter : Molecules Quantum mechanics can be used to calculate any property of a molecule The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is, E = Ψ H Ψ Ψ Ψ 1) At first this seems like
More informationLecture 13: The Actual Shell Model Recalling schematic SM Single-particle wavefunctions Basis states Truncation Sketch view of a calculation
Lecture 13: The Actual Shell Model Recalling schematic SM Single-particle wavefunctions Basis states Truncation Sketch view of a calculation Lecture 13: Ohio University PHYS7501, Fall 2017, Z. Meisel (meisel@ohio.edu)
More informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationBeyond the Hartree-Fock Approximation: Configuration Interaction
Beyond the Hartree-Fock Approximation: Configuration Interaction The Hartree-Fock (HF) method uses a single determinant (single electronic configuration) description of the electronic wavefunction. For
More information5.61 Physical Chemistry Final Exam 12/16/09. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry Physical Chemistry
5.6 Physical Chemistry Final Exam 2/6/09 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry - 5.6 Physical Chemistry Final Examination () PRINT your name on the cover page. (2) It
More informationAtoms, Molecules and Solids. From Last Time Superposition of quantum states Philosophy of quantum mechanics Interpretation of the wave function:
Essay outline and Ref to main article due next Wed. HW 9: M Chap 5: Exercise 4 M Chap 7: Question A M Chap 8: Question A From Last Time Superposition of quantum states Philosophy of quantum mechanics Interpretation
More informationLecture 5. Hartree-Fock Theory. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 5 Hartree-Fock Theory WS2010/11: Introduction to Nuclear and Particle Physics Particle-number representation: General formalism The simplest starting point for a many-body state is a system of
More informationComparing and Improving Quark Models for the Triply Bottom Baryon Spectrum
Comparing and Improving Quark Models for the Triply Bottom Baryon Spectrum A thesis submitted in partial fulfillment of the requirements for the degree of Bachelor of Science degree in Physics from the
More informationLecture 4: Nuclear Energy Generation
Lecture 4: Nuclear Energy Generation Literature: Prialnik chapter 4.1 & 4.2!" 1 a) Some properties of atomic nuclei Let: Z = atomic number = # of protons in nucleus A = atomic mass number = # of nucleons
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 informationReikniefnafræði - Verkefni 1 Haustmisseri 2013 Kennari - Hannes Jónsson
Háskóli Íslands, raunvísindasvið Reikniefnafræði - Verkefni 1 Haustmisseri 2013 Kennari - Hannes Jónsson Guðjón Henning 11. september 2013 A. Molecular orbitals and electron density of H 2 Q1: Present
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More informationIdentical Particles in Quantum Mechanics
Identical Particles in Quantum Mechanics Chapter 20 P. J. Grandinetti Chem. 4300 Nov 17, 2017 P. J. Grandinetti (Chem. 4300) Identical Particles in Quantum Mechanics Nov 17, 2017 1 / 20 Wolfgang Pauli
More informationQuantum Mechanics Solutions. λ i λ j v j v j v i v i.
Quantum Mechanics Solutions 1. (a) If H has an orthonormal basis consisting of the eigenvectors { v i } of A with eigenvalues λ i C, then A can be written in terms of its spectral decomposition as A =
More informationPauli Deformation APPENDIX Y
APPENDIX Y Two molecules, when isolated say at infinite distance, are independent and the wave function of the total system might be taken as a product of the wave functions for the individual molecules.
More informationMany-Body Approaches to Quantum Dots
Many-Body Approaches to Patrick Merlot Department of Computational Physics November 4, 2009 Many-Body Approaches to Table of contents 1 2 3 Many-Body Approaches to What is a quantum dot? Physics of QDs
More informationQUANTUM MECHANICS I PHYS 516. Solutions to Problem Set # 1
QUANTUM MECHANICS I PHYS 516 Solutions to Problem Set # 1 1. Scaling: Bohr computed the energy level spectrum of the hydrogen atom using the Old Quantum Theory, Heisenberg did the same using Matrix Mechanics,
More informationUniversity of Michigan Physics Department Graduate Qualifying Examination
Name: University of Michigan Physics Department Graduate Qualifying Examination Part II: Modern Physics Saturday 17 May 2014 9:30 am 2:30 pm Exam Number: This is a closed book exam, but a number of useful
More informationLECTURE 4: HOW TO GENERATE ELECTRONIC CONFIGURATIONS FOR ATOMS AND IONS
LECTURE 4: HOW TO GENERATE ELECTRONIC CONFIGURATIONS FOR ATOMS AND IONS We are about to do something incredibly useful create the electronic configurations for every kind of neutral or charged atom we
More informationAlkali metals show splitting of spectral lines in absence of magnetic field. s lines not split p, d lines split
Electron Spin Electron spin hypothesis Solution to H atom problem gave three quantum numbers, n,, m. These apply to all atoms. Experiments show not complete description. Something missing. Alkali metals
More information