CASSCF and NEVPT2 calculations: Ground and excited states of multireference systems. A case study of Ni(CO)4 and the magnetic system NArO
|
|
- Emily Black
- 5 years ago
- Views:
Transcription
1 CASSCF and NEVPT2 calculations: Ground and excited states of multireference systems. A case study of Ni(CO)4 and the magnetic system NArO The ground states of many molecules are often well described by a single determinant. In such cases it is often most appropriate to use so-called single reference methods. Density functional methods (DFT) are widely used, and excited states can be described by the corresponding time-depent DFT (TDDFT approaches). Alternatively one can start from Hartree-Fock (HF), and apply coupled cluster methods. For ground states the CCSD(T) aproach is most suitable, and in ORCA one has very efficient local CC approaches that go by the name DLPNO-CCSD(T). For excited states one can use EOM-CC or STEOM-CC. In ORCA, one can run efficient DLPNO-STEOM or bt-steom calculations, that employ the efficient DLPNO technology. These techniques are described in other labs. Here I will focus on systems or states that require a different (more complicated) treatment. They are characterized by a (small) set of active orbitals, and wave functions that are qualitatively well described by determinants that vary in the occupation of these active orbitals. They can have additional (deeper lying) orbitals that are all doubly occupied. Such cases often arise if we are intersted in both ground and excited states that can both be described by active space configurations. Examples of such systems are biradicals (e.g benzene with two H-atoms removed), open-shell atoms, magnetic systems (essentiall open-shell atoms connected by bridging ligands), and many systems involving transition metal atoms. The latter systems often have a number of close lying orbitals comprised of the d orbitals of the metal. The same can happen for f-shell atoms (actanides or lathanides). The calculation proceeds in a number of steps. The first step is the selection of a space of active orbitals, and optimization of the orbitals. This is done using a complete active space CASSCF calculation. This type of calculation can provide qualitatively accurate wavefunctions, and geometries, similar to HF for single reference molecules. However, to get more accurate energetics and excitation energies one should include additional electron correlation effects. In ORCA the NEVPT2 approach is an efficient and accurate way to do this. In the Nooijen lab we have developed the MREOM Coupled cluster approach. This is more expensive, but is expected to be more accurate than the perturbative NEVPT2 approach. This will be discussed in other labs. Other MRCI and DDCI like approaches are available in ORCA. However to pursue MRCI calculation the MOLPRO package may be more suitable. Here I will focus on the CASSCF and NEVPT2 approach and some of their variations. Detailed Case study Ni(CO)4 I will treat one example in detail, the tetrahedral Ni(CO)4 molecule. This molecule has a high symmetry, and we would like to preserve the degeneracy patterns of the true energy levels. The first step in the calculation is the creation of a suitable CASSCF calculation to calculate a number of low-lying electronic states using one type of so-called state-averaged calculation. All of the calculations can be found in the NEVPT2 directory on chem400a/b.
2 The calculations discussed in this section are all in NEVPT2/NiCO4/10e_8o. I will also briefly discuss a different type of CAS, using 10 electrons in 10 (spatial) orbitals. A good starting point to get an idea of the orbitals is often a proceeding DFT calculation. I started with a fixed geometry (that I obtained from elsewhere), and run the following DFT bp86 calculation, asking to print the orbitals (normalprint): Input file bp86_orb.inp! bp86 def2-svp NormalPrint * xyz 0 1 Ni C C C C O O O O * This yields the following orbital energies around the Fermi level: We notice a nice degeneray pattern showing the E and T representations which are 2- and 3-fold degenerate respectively. The 5 occupied orbitals (2.0000) are d-orbitals on Ni, bonded to the carbonyl COs. Likewise the first five orbitals above the fermi level (0.0000) have significant d-type character. This can be seen by inspecting the orbital character in this section of the output: LOEWDIN REDUCED ORBITAL POPULATIONS PER MO
3 THRESHOLD FOR PRINTING IS 0.1% If we scroll down towards the fermi level we see: Ni pz Ni px Ni py Ni dz Ni dxz Ni dyz Ni dx2y Ni dxy C s C pz C px Showing the predominant d-character of the last 5 orbitals (37-41). Likewise if we continue and examine the virtual orbitals: we see some d-character in orbitals 42-46, but not in other low-lying virtuals, e.g Orbital 50 is the Ni 4s orbital: Ni pz Ni px Ni py Ni dz Ni dxz Ni dyz Ni dx2y Ni dxy This suggests we can use an active space with 10 orbitals (5 occ + 5 vrt), and 10 electrons. Let us now do an initial CASSCF calculation, optimizing just one singlet state Input file CAS_a.inp!CASSCF def2-svp norb 10 nroots 1
4 We read in the orbitals from the previous pb86 calculation ( section), and calculate the ground state CAS. I am doing this all with a small basis set, as we do this to learn the procedure. I would prefer to use the def2-tzvp basis set in general for such calculations. Even then it makes sense to first explore things with a smaller basis set first. We can examine the relevant parts of the output. Look for the statement Final CASSCF energies, Final CASSCF energy : Eh ev ORBITAL ENERGIES NO OCC E(Eh) E(eV) Scrolling down, to were the occupation numbers are partial (no longer 2.000): We will be interested in a number of low-lying singlet states, and we need to figure out how many, and what the size of the final active space will be. To this we use the orbitals from the previous CAS_a calculation to start from, and ask for many states in the singlet manifold. However we do not optimize the orbitals (NoIter). Here is the input file Input File CASCI_a.inp!CASSCF def2-svp NoIter
5 norb 10 nroots 20 We asked for 20 roots here. To examine the relevant part of the output search for Final CASSCF energy, and then look for (or scroll down to) CAS-SCF STATES FOR BLOCK 1 MULT= 1 NROOTS= ROOT 0: E= Eh [ 0]: [ 2916]: [ 1060]: [ 841]: [ 168]: [ 27]: [ 217]: [ 2965]: [ 3184]: [ 792]: [ 76]: [ 436]: [ 6]: [ 3936]: [ 1812]: [ 6060]: ROOT 1: E= Eh ev cm** [ 2907]: [ 785]: [ 2]: [ 21]: This lists a number of excitation energies. We see nice degeneracy patterns, 3 states at ev, 3 at , 2 at , and so on. These energies are very high. Let us make a second attempt, by only asking for 4 roots in a CASSCF calculation (optimize orbitals).!casscf def2-svp norb 10
6 nroots 4 We can run the CASCI_a2 calculation afterwards, asking to calculate 20 roots, but now with the optimized orbtals from the above 4-state CAS(10e,10o). Looking for the converged energies we now find the excitation energy pattern: 1* ev, 3*6.028 ev, 2*6.169, 3*6.475, 3* Then there is a jump to ev This indicates we might look at 12 low-lying states. Moreover looking at the occupation numbers before final CASSCF energy, we see N(occ)= This indicates that the last two orbitals are not strongly occupied, and we might reduce the CAS to 12 states, 10e, 8o This leads to the following CAS_b input file:!casscf def2-svp norb 8 nroots 12 Which has final occupations: N(occ)= You can also find the final state energies, 0.000, 4.938, and so on. We can observe that the degeneracy pattern is not perfect. In particular the last ones ( etc.). The reson is that the CASSCF is not sufficiently converged. As long as we are stll exploring this makes perfect sense. However, once we are happy with the CAS it may make sense to converge energies more sharply (ExtermeSCF), as in CAS_c.inp!CASSCF def2-svp ExtremeSCF
7 norb 8 nroots 12 You can check occupation numbers and excitation energies to see that degeneracy patterns are nearly perfect now which shows perfect degeneracies. You can also examine the energies of the 12 states (just above Final CASSCF energy): CI-ITERATION 16: ( 0.03) CI-PROBLEM SOLVED DENSITIES MADE E(CAS)= Eh DE= Energy gap subspaces: Ext-Act = Act-Int = N(occ)= g = Max(G)= Rot=46,43 We have found the CASSCF solution we want. Now it is strightforward to run NEVPT2 calculations: Input file NEVPT2.inp!CASSCF def2-svp NoIter
8 norb 8 nroots 12 NEVPT2 SC We read in the orbitals from the previous CAS_c calculation, and can use NoIter therefore. The NEVPT2 keyword (lots of documentation in the ORCA manual) can be included in the CASSCF input section. NEVPT2 SC is the fastest version of the NEVPT2 variants. This yields the following output in the NEVPT2 section: NEVPT2 TRANSITION ENERGIES LOWEST ROOT (ROOT 0, MULT 1) = Eh ev STATE ROOT MULT DE/a.u. DE/eV DE/cm**-1 1: : : : : : : : : : : You can also look at the absorbtion (or CD) spectrum. At thi spoint we can examine some other variations of NEVPT2 calculations. The first is a potentially more accurate variant, quasi-degenerate NEVPT2, using the input file qd_nevpt2.inp!casscf def2-svp NoIter
9 norb 8 nroots 12 NEVPT2 SC NEVPT qdtype 1 which yields both the original NEVPT2 and the corrected QD-NEVPT2 energies: NEVPT2 TRANSITION ENERGIES LOWEST ROOT (ROOT 0, MULT 1) = Eh ev STATE ROOT MULT DE/a.u. DE/eV DE/cm**-1 1: : : : : : : : : : : QD-NEVPT2 TRANSITION ENERGIES LOWEST ROOT (ROOT 0, MULT 1) = Eh ev STATE ROOT MULT DE/a.u. DE/eV DE/cm**-1 1: : : : : : : : : : :
10 A final possibly very useful option is to use dlpno_nevpt2. This can be a bit expensive for small molecules, but it is much more efficient when you get to really large molecules. Here is the input file:!dlpno-nevpt2 def2-svp def2-svp/c Nofrozencore norb 8 nroots 12 You need to specify a fitting basis set (def2_svp/c here), and you can set the method on the first input line. The calculation takes more than an hour. The efficieny kicks in for larger molecules. Nofrozencore is recommed for all NEVPT2 calculations. More Examples. For the Ni(CO)4 molecule the choice of CAS was not so clear. We decided to do 12 states in a CAS(10e, 8o), but a clear alternative was to use 12 states with (10e,10o). These calculations are also done in the NEVPT2/NiCO4/10e_10o directory. You can examine the differences. This shows such calculations can be delicate. Finally the calculations may be suitable for magnetic systems. The NEVPT2 calculations will run the fastest among available multireference methods. I provided results for a simple NArO system. This is considered in more detail in the MREOM write-up and directory. The NEVPT2 results compare fairly well with MREOM energies and they are often an enormous improvement over CASSCF energies themselves. The DDCI approach is also a well-established approach for magnetic systems, and inclusion of spin orbit effects are straigtforward. It is possible to use Spin-Orbit Coupling with NEVPT2, but this needs a detailed investigation. I did provide in an input file avg_nevpt2_soc, which uses the avg_nevpt2 energies in a CASCI_soc calculation, as discussed in the SOC section of the ORCA manual. I did not make a detailed write-up for this system. You can follow the input and output files in the NEVPT2/NArO directory, as well as the joball file (which has comments), which reveals the full history.
11
Let me go through the basic steps of an MREOM calculation in ORCA. On chem440a/b the calculations can be found under ~nooijen/orca_examples_2017/mreom
MREOM calculations in ORCA. The procedure to run MREOM calculations in ORCA is similar to running MREOM calculations in ACES. There are essentially three steps. 1. Determine a good set of CASSCF reference
More informationMRCI calculations in MOLPRO
1 MRCI calculations in MOLPRO Molpro is a software package written in Fortran and maintained by H.J. Werner and P.J. Knowles. It is often used for performing sophisticated electronic structure calculations,
More informationHints on Using the Orca Program
Computational Chemistry Workshops West Ridge Research Building-UAF Campus 9:00am-4:00pm, Room 009 Electronic Structure - July 19-21, 2016 Molecular Dynamics - July 26-28, 2016 Hints on Using the Orca Program
More informationPractical Issues on the Use of the CASPT2/CASSCF Method in Modeling Photochemistry: the Selection and Protection of an Active Space
Practical Issues on the Use of the CASPT2/CASSCF Method in Modeling Photochemistry: the Selection and Protection of an Active Space Roland Lindh Dept. of Chemistry Ångström The Theoretical Chemistry Programme
More informationTheoretical UV/VIS Spectroscopy
Theoretical UV/VIS Spectroscopy Why is a Ruby Red When Chromium Oxide is Green? How Does a Ruby Laser Work? Goals of this Exercise: - Calculation of the energy of electronically excited states - Understanding
More informationPerformance of Hartree-Fock and Correlated Methods
Chemistry 460 Fall 2017 Dr. Jean M. Standard December 4, 2017 Performance of Hartree-Fock and Correlated Methods Hartree-Fock Methods Hartree-Fock methods generally yield optimized geomtries and molecular
More informationBuilding a wavefunction within the Complete-Active. Cluster with Singles and Doubles formalism: straightforward description of quasidegeneracy
Building a wavefunction within the Complete-Active Active-Space Coupled-Cluster Cluster with Singles and Doubles formalism: straightforward description of quasidegeneracy Dmitry I. Lyakh (Karazin Kharkiv
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 informationRethinking Hybridization
Rethinking Hybridization For more than 60 years, one of the most used concepts to come out of the valence bond model developed by Pauling was that of hybrid orbitals. The ideas of hybridization seemed
More informationQUANTUM CHEMISTRY FOR TRANSITION METALS
QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline
More informationPractical Advice for Quantum Chemistry Computations. C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology
Practical Advice for Quantum Chemistry Computations C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Choice of Basis Set STO-3G is too small 6-31G* or 6-31G** 6 probably
More informationMolecular Orbital Theory (MOT)
Molecular Orbital Theory (MOT) In this section, There are another approach to the bonding in metal complexes: the use of molecular orbital theory (MOT). In contrast to crystal field theory, the molecular
More informationComputational Chemistry I
Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of
More informationElectron Correlation Methods
Electron Correlation Methods HF method: electron-electron interaction is replaced by an average interaction E HF c = E 0 E HF E 0 exact ground state energy E HF HF energy for a given basis set HF E c
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 informationMethods for Treating Electron Correlation CHEM 430
Methods for Treating Electron Correlation CHEM 430 Electron Correlation Energy in the Hartree-Fock approximation, each electron sees the average density of all of the other electrons two electrons cannot
More informationLECTURE 2 DEGENERACY AND DESCENT IN SYMMETRY: LIGAND FIELD SPLITTINGS AND RELATED MATTERS
SYMMETRY II. J. M. GOICOECHEA. LECTURE 2. 1 LECTURE 2 DEGENERACY AND DESCENT IN SYMMETRY: LIGAND FIELD SPLITTINGS AND RELATED MATTERS 2.1 Degeneracy When dealing with non-degenerate symmetry adapted wavefunctions
More informationAN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University
AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,
More informationSupporting Information. Nonclassical Single-State Reactivity of an Oxo- Iron(IV) Complex Confined to Triplet Pathways
Supporting Information for Nonclassical Single-State Reactivity of an Oxo- Iron(IV) Complex Confined to Triplet Pathways Claudia Kupper, ǁ Bhaskar Mondal, ǁ Joan Serrano-Plana, Iris Klawitter, Frank Neese,
More information5.4. Electronic structure of water
5.4. Electronic structure of water Water belongs to C 2v point group, we have discussed the corresponding character table. Here it is again: C 2v E C 2 σ v (yz) σ v (xz) A 1 1 1 1 1 A 2 1 1-1 -1 B 1 1-1
More informationThe MCSCF Method *, Molecular Orbitals, Reference Spaces and COLUMBUS Input
The MCSCF Method *, Molecular Orbitals, Reference Spaces and COLUMBUS Input Hans Lischka University of Vienna *Excerpt of a course presented by R. Shepard, Argonne National Laboratory, at the Workshop
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationNWChem: Coupled Cluster Method (Tensor Contraction Engine)
NWChem: Coupled Cluster Method (Tensor Contraction Engine) Why CC is important?! Correlation effects are important!! CC is size-extensive theory: can be used to describe dissociation processes.! Higher-order
More informationOther Crystal Fields
Other Crystal Fields! We can deduce the CFT splitting of d orbitals in virtually any ligand field by " Noting the direct product listings in the appropriate character table to determine the ways in which
More informationABC of DFT: Hands-on session 1 Introduction into calculations on molecules
ABC of DFT: Hands-on session 1 Introduction into calculations on molecules Tutor: Alexej Bagrets Wann? 09.11.2012, 11:30-13:00 Wo? KIT Campus Nord, Flachbau Physik, Geb. 30.22, Computerpool, Raum FE-6
More informationIntroduction to Electronic Structure Theory
CSC/PRACE Spring School in Computational Chemistry 2017 Introduction to Electronic Structure Theory Mikael Johansson http://www.iki.fi/~mpjohans Objective: To get familiarised with the, subjectively chosen,
More informationXMVB 3.0 A Right Way To Do Valence Bond Calculations. Zhenhua Chen Xiamen University
XMVB 3.0 A Right Way To Do Valence Bond Calculations Zhenhua Chen Xiamen University XMVB 3.0 Xiamen Valence Bond An Ab Initio Non-orthogonal Valence Bond Program Version 3.0 Lingchun Song, Zhenhua Chen,
More informationσ u * 1s g - gerade u - ungerade * - antibonding σ g 1s
One of these two states is a repulsive (dissociative) state. Other excited states can be constructed using linear combinations of other orbitals. Some will be binding and others will be repulsive. Thus
More informationElectron Correlation - Methods beyond Hartree-Fock
Electron Correlation - Methods beyond Hartree-Fock how to approach chemical accuracy Alexander A. Auer Max-Planck-Institute for Chemical Energy Conversion, Mülheim September 4, 2014 MMER Summerschool 2014
More informationExchange coupling can frequently be understood using a simple molecular orbital approach.
6.4 Exchange Coupling, a different perspective So far, we ve only been looking at the effects of J on the magnetic susceptibility but haven t said anything about how one might predict the sign and magnitude
More informationElectronic structure of correlated electron systems. Lecture 2
Electronic structure of correlated electron systems Lecture 2 Band Structure approach vs atomic Band structure Delocalized Bloch states Fill up states with electrons starting from the lowest energy No
More 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 informationB. Electron Deficient (less than an octet) H-Be-H. Be does not need an octet Total of 4 valence electrons
B. Electron Deficient (less than an octet) e.g. BeH 2 H-Be-H Be does not need an octet Total of 4 valence electrons Not the same as unsaturated systems that achieve the 8e - (octet) through the formation
More informationExtended Wavefunction Analysis for Multireference Methods
Extended Wavefunction Analysis for Multireference Methods Felix Plasser González Research Group Institute for Theoretical Chemistry, University of Vienna, Austria Vienna, 1 st April 2016 Introduction Analysis
More informationComputational chemistry with GAMESS: a very brief overview with examples
Computational chemistry with GAMESS: a very brief overview with examples PHY-6120 Molecular Physics (Spring 2015), UConn Phys. Dept. Feb 17 th 2015 H = ħ2 2μ i Intro: V(R) for diatomic molecules + k Z
More 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 informationThis is a simple input file for the calculation of NMR chemical shieldings for a given molecule using the B3LYP functional and def2-tzvpp basis set:
Computing NMR parameters using ORCA This practical comes with a short lecture on the basics of the computation of NMR parameters using standard electronic structure theory methods. By now you should have
More informationNWChem: Coupled Cluster Method (Tensor Contraction Engine)
NWChem: Coupled Cluster Method (ensor Contraction Engine) What we want to solve H Ψ = E Ψ Many Particle Systems Molecular/Atomic Physics, Quantum Chemistry (electronic Schrödinger equations) Solid State
More informationAb initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier
Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential
More informationAdvanced usage of MOLCAS. Ex. 1. Usage of Symmetry in Molcas
Advanced usage of MOLCAS Ex 1. Symmetry in Molcas (20 min). Ex 2. Transition state optimization (30 min). Ex 3. Ways to run RASSCF program (30 min). Ex 4. Optional. Using EXPBAS module (20 min) Ex 5. Optional.
More informationMolecular Orbitals for Ozone
Molecular Orbitals for Ozone Purpose: In this exercise you will do semi-empirical molecular orbital calculations on ozone with the goal of understanding the molecular orbital print out provided by Spartan
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 informationRelativistic and correlated calculations on the ground, excited, and ionized states of iodine
Relativistic and correlated calculations on the ground, excited, and ionized states of iodine W. A. de Jong, L. Visscher, a) and W. C. Nieuwpoort Laboratory for Chemical Physics and Materials Science Centre,
More informationSupporting Information: Predicting the Ionic Product of Water
Supporting Information: Predicting the Ionic Product of Water Eva Perlt 1,+, Michael von Domaros 1,+, Barbara Kirchner 1, Ralf Ludwig 2, and Frank Weinhold 3,* 1 Mulliken Center for Theoretical Chemistry,
More informationSupplemental Material: Experimental and Theoretical Investigations of the Electronic Band Structure of Metal-Organic Framework of HKUST-1 Type
Supplemental Material: Experimental and Theoretical Investigations of the Electronic Band Structure of Metal-Organic Framework of HKUST-1 Type Zhigang Gu, a Lars Heinke, a,* Christof Wöll a, Tobias Neumann,
More informationDensity matrix renormalization group calculations on relative energies of transition metal complexes and clusters
THE JOURNAL OF CHEMICAL PHYSICS 128, 014104 2008 Density matrix renormalization group calculations on relative energies of transition metal complexes and clusters Konrad H. Marti, Irina Malkin Ondík, Gerrit
More informationChemistry 4560/5560 Molecular Modeling Fall 2014
Final Exam Name:. User s guide: 1. Read questions carefully and make sure you understand them before answering (if not, ask). 2. Answer only the question that is asked, not a different question. 3. Unless
More informationMolecular orbitals for σbonding in T d complexes
Molecular orbitals for σbonding in T d complexes The set of n A B σ bonds in AB n (T d n = 4) molecules are often thought of as independent entities. The concept of MO s allows us to begin with a very
More informationMulticonfigurational Quantum Chemistry. Björn O. Roos as told by RL Department of Theoretical Chemistry Chemical Center Lund University Sweden
Multiconfigurational Quantum Chemistry Björn O. Roos as told by RL Department of Theoretical Chemistry Chemical Center Lund University Sweden April 20, 2009 1 The Slater determinant Using the spin-orbitals,
More informationElectron transport through Shiba states induced by magnetic adsorbates on a superconductor
Electron transport through Shiba states induced by magnetic adsorbates on a superconductor Michael Ruby, Nino Hatter, Benjamin Heinrich Falko Pientka, Yang Peng, Felix von Oppen, Nacho Pascual, Katharina
More informationFigure 1: Transition State, Saddle Point, Reaction Pathway
Computational Chemistry Workshops West Ridge Research Building-UAF Campus 9:00am-4:00pm, Room 009 Electronic Structure - July 19-21, 2016 Molecular Dynamics - July 26-28, 2016 Potential Energy Surfaces
More informationQUANTUM CHEMISTRY PROJECT 3: ATOMIC AND MOLECULAR STRUCTURE
Chemistry 460 Fall 2017 Dr. Jean M. Standard November 1, 2017 QUANTUM CHEMISTRY PROJECT 3: ATOMIC AND MOLECULAR STRUCTURE OUTLINE In this project, you will carry out quantum mechanical calculations of
More informationElectronic Supplementary Material (ESI) for Chemical Science This journal is The Royal Society of Chemistry 2012
Fig. S1 CASSCF (13,10) active space orbitals with Ru-Ru distance of 2.4 Å. Occupation numbers are on the left and energies in Hartrees are on the right of each orbital. The δ orbital is also included here,
More informationHow to compute the projected density of states (PDOS) Javier Junquera
How to compute the projected density of states (PDOS) Javier Junquera Density Of States (DOS) the number of one-electron levels between E and E + de SrTiO 3 bulk Units: (Energy) -1 Projected Density Of
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 5 The Hartree-Fock method C.-K. Skylaris Learning outcomes Be able to use the variational principle in quantum calculations Be able to construct Fock operators
More information2~:J~ -ryej- r- 2 Jr. A - f3. sr(djk nv~tor rn~ +~ rvjs (::-CJ) ::;-1-.'--~ -. rhd. ('-.Ji.L.~ )- r'-d)c, -r/~ JJr - 2~d ~2-Jr fn'6.
.~, ~ I, sr(djk nv~tor rn~ +~ rvjs (::-CJ) ::;-1-.'--~ -. rhd. ('-.Ji.L.~ )- r'-d)c, -r/~ JJr - 2~d ~2-Jr fn'6.)1e'" 21t-ol Je C'...-------- lj-vi, J? Jr Jr \Ji 2~:J~ -ryej- r- 2 Jr A - f3 c _,~,= ~,.,w._..._.
More informationThis is called a singlet or spin singlet, because the so called multiplicity, or number of possible orientations of the total spin, which is
9. Open shell systems The derivation of Hartree-Fock equations (Chapter 7) was done for a special case of a closed shell systems. Closed shell means that each MO is occupied by two electrons with the opposite
More informationChem 673, Problem Set 5 Due Thursday, November 29, 2007
Chem 673, Problem Set 5 Due Thursday, November 29, 2007 (1) Trigonal prismatic coordination is fairly common in solid-state inorganic chemistry. In most cases the geometry of the trigonal prism is such
More informationACES2 Labs. Part B: Calculating electronically excited, ionized, and electro-attached states using equation of motion coupled cluster calculations.
ACES2 Labs. Part B: Calculating electronically excited, ionized, and electro-attached states using equation of motion coupled cluster calculations. I. Electronically excited states. Purpose: - Calculation
More informationStudy of Iron Dimers Reveals Angular Dependence of Valence- to- Core X- ray Emission Spectra
Supporting Information for: Study of Iron Dimers Reveals Angular Dependence of Valence- to- Core X- ray Emission Spectra Christopher J. Pollock, a Kyle M. Lancaster, b Kenneth D. Finkelstein, c Serena
More informationABC of DFT: Hands-on session 2 Molecules: structure optimization, visualization of orbitals, charge & spin densities
ABC of DFT: Hands-on session 2 Molecules: structure optimization, visualization of orbitals, charge & spin densities Tutor: Alexej Bagrets Wann? 16.11.2012, 11:30-13:00 Wo? KIT Campus Nord, Flachbau Physik,
More informationElectronic Spectroscopy of Polyatomics
Electronic Spectroscopy of Polyatomics We shall discuss the electronic spectroscopy of the following types of polyatomic molecules: 1. general AH 2 molecules, A = first-row element 2. formaldehyde 3. benzene
More informationChem 673, Problem Set 5 Due Thursday, December 1, 2005
otton, Problem 9.3 (assume D 4h symmetry) Additional Problems: hem 673, Problem Set 5 Due Thursday, December 1, 2005 (1) Infrared and Raman spectra of Benzene (a) Determine the symmetries (irreducible
More informationLecture 4: methods and terminology, part II
So theory guys have got it made in rooms free of pollution. Instead of problems with the reflux, they have only solutions... In other words, experimentalists will likely die of cancer From working hard,
More informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More informationFree-Ion Terms to Ligand-field Terms
Free-Ion Terms to Ligand-field Terms! Orbital term symbols for free atoms and ions are identical to symbols for irreducible representations in R 3. " The irreducible representations of R 3 include all
More informationThe Gutzwiller Density Functional Theory
The Gutzwiller Density Functional Theory Jörg Bünemann, BTU Cottbus I) Introduction 1. Model for an H 2 -molecule 2. Transition metals and their compounds II) Gutzwiller variational theory 1. Gutzwiller
More informationElectronic structure / bonding in d-block complexes
LN05-1 Electronic structure / bonding in d-block complexes Many, many properties of transition metal complexes (coordination number, structure, colour, magnetism, reactivity) are very sensitive to the
More informationLecture 9: Molecular Orbital theory for hydrogen molecule ion
Lecture 9: Molecular Orbital theory for hydrogen molecule ion Molecular Orbital Theory for Hydrogen Molecule Ion We have seen that the Schrödinger equation cannot be solved for many electron systems. The
More informationMulti-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005
Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson
More information5.61 Physical Chemistry Exam III 11/29/12. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry Physical Chemistry.
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry - 5.61 Physical Chemistry Exam III (1) PRINT your name on the cover page. (2) It is suggested that you READ THE ENTIRE EXAM before
More informationIFM Chemistry Computational Chemistry 2010, 7.5 hp LAB2. Computer laboratory exercise 1 (LAB2): Quantum chemical calculations
Computer laboratory exercise 1 (LAB2): Quantum chemical calculations Introduction: The objective of the second computer laboratory exercise is to get acquainted with a program for performing quantum chemical
More informationIntroduction to multiconfigurational quantum chemistry. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Notations
More informationChemistry 3502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2003 Christopher J. Cramer. Lecture 25, November 5, 2003
Chemistry 3502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2003 Christopher J. Cramer Lecture 25, November 5, 2003 (Some material in this lecture has been adapted from Cramer, C.
More informationNWChem: Hartree-Fock, Density Functional Theory, Time-Dependent Density Functional Theory
NWChem: Hartree-Fock, Density Functional Theory, Time-Depent Density Functional Theory Hartree-Fock! Functionality! Input! Wavefunctions! Initial MO vectors! Direct and semidirect algorithms! Convergence,
More information( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r)
Born Oppenheimer Approximation: Ĥ el ( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r) For a molecule with N electrons and M nuclei: Ĥ el What is E el (R)? s* potential surface Reaction Barrier Unstable intermediate
More informationA Rigorous Introduction to Molecular Orbital Theory and its Applications in Chemistry. Zachary Chin, Alex Li, Alex Liu
A Rigorous Introduction to Molecular Orbital Theory and its Applications in Chemistry Zachary Chin, Alex Li, Alex Liu Quantum Mechanics Atomic Orbitals and Early Bonding Theory Quantum Numbers n: principal
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 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 informationActivity Molecular Orbital Theory
Activity 201 9 Molecular Orbital Theory Directions: This Guided Learning Activity (GLA) discusses the Molecular Orbital Theory and its application to homonuclear diatomic molecules. Part A describes the
More informationBridging Scales Through Wavefunction Analysis
Bridging Scales Through Wavefunction Analysis Felix Plasser Institute for Theoretical Chemistry, University of Vienna Excited States Bridging Scales Marseille, November 7 10, 2016 F. Plasser Wavefunction
More informationUsing Symmetry to Generate Molecular Orbital Diagrams
Using Symmetry to Generate Molecular Orbital Diagrams review a few MO concepts generate MO for XH 2, H 2 O, SF 6 Formation of a bond occurs when electron density collects between the two bonded nuclei
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 informationActivity Molecular Orbital Theory
Activity 201 9 Molecular Orbital Theory Directions: This Guided Learning Activity (GLA) discusses the Molecular Orbital Theory and its application to homonuclear diatomic molecules. Part A describes the
More informationCalculations of band structures
Chemistry and Physics at Albany Planning for the Future Calculations of band structures using wave-function based correlation methods Elke Pahl Centre of Theoretical Chemistry and Physics Institute of
More informationLecture B6 Molecular Orbital Theory. Sometimes it's good to be alone.
Lecture B6 Molecular Orbital Theory Sometimes it's good to be alone. Covalent Bond Theories 1. VSEPR (valence shell electron pair repulsion model). A set of empirical rules for predicting a molecular geometry
More informationSUPPLEMENTARY INFORMATION
DOI: 10.1038/NCHEM.1677 Entangled quantum electronic wavefunctions of the Mn 4 CaO 5 cluster in photosystem II Yuki Kurashige 1 *, Garnet Kin-Lic Chan 2, Takeshi Yanai 1 1 Department of Theoretical and
More information0 belonging to the unperturbed Hamiltonian H 0 are known
Time Independent Perturbation Theory D Perturbation theory is used in two qualitatively different contexts in quantum chemistry. It allows one to estimate (because perturbation theory is usually employed
More informationReikniefnafræði - Verkefni 2 Haustmisseri 2013 Kennari - Hannes Jónsson
Háskóli Íslands, raunvísindasvið Reikniefnafræði - Verkefni 2 Haustmisseri 2013 Kennari - Hannes Jónsson Guðjón Henning 18. september 2013 1 A. Molecular orbitals of N 2 Q1: Display all the MOs for N 2
More informationIntroduction to Hartree-Fock calculations in Spartan
EE5 in 2008 Hannes Jónsson Introduction to Hartree-Fock calculations in Spartan In this exercise, you will get to use state of the art software for carrying out calculations of wavefunctions for molecues,
More information4 Post-Hartree Fock Methods: MPn and Configuration Interaction
4 Post-Hartree Fock Methods: MPn and Configuration Interaction In the limit of a complete basis, the Hartree-Fock (HF) energy in the complete basis set limit (ECBS HF ) yields an upper boundary to the
More informationElectronic communication through molecular bridges Supporting Information
Electronic communication through molecular bridges Supporting Information Carmen Herrmann and Jan Elmisz Institute of Inorganic and Applied Chemistry, University of Hamburg, Martin-Luther-King-Platz 6,
More informationAN INTRODUCTION TO MCSCF: PART 2
AN INTRODUCTION TO MCSCF: PART 2 ORBITAL APPROXIMATION Hartree product (hp) expressed as a product of spinorbitals ψ ι = φ i σ i φ i = space orbital, σ i = spin function (α,β) Pauli Principle requires
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 informationChem120a : Exam 3 (Chem Bio) Solutions
Chem10a : Exam 3 (Chem Bio) Solutions November 7, 006 Problem 1 This problem will basically involve us doing two Hückel calculations: one for the linear geometry, and one for the triangular geometry. We
More informationChapter 9: Multi- Electron Atoms Ground States and X- ray Excitation
Chapter 9: Multi- Electron Atoms Ground States and X- ray Excitation Up to now we have considered one-electron atoms. Almost all atoms are multiple-electron atoms and their description is more complicated
More informationBe H. Delocalized Bonding. Localized Bonding. σ 2. σ 1. Two (sp-1s) Be-H σ bonds. The two σ bonding MO s in BeH 2. MO diagram for BeH 2
The Delocalized Approach to Bonding: The localized models for bonding we have examined (Lewis and VBT) assume that all electrons are restricted to specific bonds between atoms or in lone pairs. In contrast,
More informationPeriodic Trends in Properties of Homonuclear
Chapter 8 Periodic Trends in Properties of Homonuclear Diatomic Molecules Up to now, we have discussed various physical properties of nanostructures, namely, two-dimensional - graphene-like structures:
More informationTheoretical study of the low-lying excited singlet states of furan
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 2 8 JULY 2003 Theoretical study of the low-lying excited singlet states of furan E. V. Gromov, A. B. Trofimov, and N. M. Vitkovskaya Laboratory of Quantum
More informationMolecular Orbital Theory. WX AP Chemistry Chapter 9 Adapted from: Luis Bonilla Abel Perez University of Texas at El Paso
Molecular Orbital Theory WX AP Chemistry Chapter 9 Adapted from: Luis Bonilla Abel Perez University of Texas at El Paso Molecular Orbital Theory The goal of molecular orbital theory is to describe molecules
More information