Ensemble Generalized Kohn-Sham Theory: the Good, the Bad, and the Ugly

Total Page:16

File Type:pdf, Size:1020Kb

Ensemble Generalized Kohn-Sham Theory: the Good, the Bad, and the Ugly Ensemble generalized Kohn-Sham theory: the good, the bad, and the ugly Tim Gould1 and Leeor Kronik2 1)Qld Micro- and Nanotechnology Centre, Griffith University, Nathan, Qld 4111, Australia 2)Department of Materials and Interfaces, Weizmann Institute of Science, Rehovoth 76100, Israel Two important extensions of Kohn-Sham (KS) theory are generalized KS theory and ensemble KS theory. The former allows for non-multiplicative potential operators and greatly facilitates practical calculations with advanced, orbital-dependent functionals. The latter allows for quantum ensembles and enables the treatment of, e.g., open systems and excited states. Here, we combine the two extensions, both formally and practically, first via an exact yet complicated formalism, then via a computationally tractable variant that involves a controlled approximation of ensemble \ghost interactions" by means of an iterative algorithm. The resulting formalism is illustrated using selected examples. This opens the door to the application of generalized KS theory in more challenging quantum scenarios and to the improvement of ensemble theories for the purpose of practical and accurate calculations. I. INTRODUCTION tures of GKS and EKS theories, though some ad hoc solutions39,41 and approximations.54,59 have been re- Kohn-Sham1,2 (KS) density functional theory (DFT) ported. This hampers the application of ensemble theo- has proven to be an indispensable tool for first principles ries, which must be solved non-self-consistently, or with calculations across a wide range of disciplines.3,4 Many less-sophisticated approximations, or using OEPs. widely used density functional approximations (DFAs), Here, we demonstrate how to combine EKS and GKS e.g. Refs. 5{10, are \hybrid functionals", i.e., they mix theories into an ensemble GKS (EGKS) theory. We show, exact (Fock) exchange with KS exchange. These are al- however, that the resulting approach is ill-suited to exist- most always applied using the non-multiplicative Fock- ing GKS implementations in its direct form. We then in- exchange operator, which places them outside the realm troduce simple, formally-motivated approximations that of KS theory, but well within generalized Kohn-Sham make EGKS tractable without losing its good features. (GKS) theory.11{14 The popularity of hybrid function- als reflcts three useful properties of hybrids: first, the inclusion of exact Fermionic exchange allows for a sys- tematically higher predictive accuracy;15 second, this is II. BACKGROUND THEORY done in a numerically straightforward fashion, whereas use of Fock exchange within \pure" KS theory must A. Ensemble Kohn-Sham theory rely on difficult to calculate optimized effective poten- tials (OEPs);16{21 third, the non-multplicative potential allows for a reduction and in some circumstances even the We first briefly introduce EKS theory before turning elimination of the KS derivative discontinuity,22,allowing to GKS and EGKS theory. For simplicity, we use a spin- one to overcome notoriously difficult problems for KS unpolarized formalism throughout. To begin, we intro- theory such as the bandgap problem23,24 or the charge duce the universal \Levy" ensemble density functional, transfer excitation problem.25{27 λ ^ ^ ^ Recent years have seen renewed interest in ensemble F [n; w] = min Tr[(T + λW )Γw] ; (1) Γ^ !n KS (EKS) theory, which broadens the scope of traditional w KS theory to include mapping to an ensemble of refer- ence non-interacting electron states, and thereby extends where λ 2 [0; 1] is an adjustable interacting-strength pa- the reach of DFT to a more diverse range of quantum rameter. Here, n is the density, w is a set of statis- tical weights that define the ensemble, T^ is the kinetic scenarios, e.g., systems with non-integer electron num- ^ ber22,28{34 and mixtures of ground- and low-lying excited energy operator and W is the electron-electron inter- 35{37 action energy operator. The ensemble density matrix, states. It can also improve the treatment of quantum ^ P systems that are technically within the reach of KS the- Γw = κ wκjκihκj, involves a weighted sum over or- ^ ory, such as those with significant spin-contamination.38 thonormal wave functions jκi. Γw ! n is short-hand for 28,33,35{59 Excited-state EDFT, In particular, it has at- Tr[^n(r)Γ^w] = n(r), wheren ^(r) is the electron-density tracted significant recent attention as a route to low cost operator. By varying the weights, w, we can obtain use- prediction of charge transfer40,56 and double40,50,53 exci- ful properties (e.g. excitation energies) of the system that tations, which are hard to predict using traditional time- are inaccessible to ground-state DFT. dependent KS theory. The Levy functional F λ provides a versatile framework To date, we are not aware of any formal theory that for defining key functionals in EKS52,55,57, specifically has been able to combine rigorously the attractive fea- the KS kinetic energy (Ts), Hartree-exchange (Hx), and 2 correlation (c), functionals, respectively: where v is the external potential and vR is a multiplica- tive effective potential, can then be used to find orbitals 0 Ts[n; w] =F [n; w] ; (2) f'g and density n. We obtain different operatorsv ^S and λ F [n; w] − Ts[n; w] potentials vR for different orbital functionals S[f'g] that EHx[n; w] = lim ; (3) are invariant to unitary transformations of the orbitals. λ!0+ λ 1 We use f'g to indicate orbital solutions of Eq. (8), as Ec[n; w] =F [n; w] − Ts[n; w] − EHx[n; w] : (4) opposed to fφg of Eq. (7). Each is a functional of the electron density n and weights The existence of the GKS equation is deduced by defin- 1 w, denoted by [n; w]. ing: S[n] = minf'g!n S[f'g], and RS[n] = F [n] − S[n], where we use non-calligraphic F 1 to represent the special It has recently been shown that Ts[n; w] ≡ Tr[Γ^s;wT^] 55 case of Eq. (1) for pure ground states (w0 = 1, wκ6=0 = 0). and EHx[n; w] ≡ Tr[Γ^s;wW^ ], for some non-interacting Then, the ground-state energy, E0 = minnfS[n]+RS[n]+ ensemble density matrix Γ^s;w, which is formed on con- V [n]g = minf'gfS[f'g] + Rs[n'] + V [n']g, is minimized figuration state functions. Ts and EHx are thereby func- for the density n ≡ n' ≡ n0. Uniqueness (in non- tionals of a set of orbitals, fφg: degenerate systems) follows from the external potential Z v being uniquely defined by the ground state density n0. X w ∗ ^ X w R 1 Ts[fφg] = fi drφi (r)tφi(r) ≡ fi ti,φ ; (5) Here V [n] = drn(r)v(r) and we used S + RS = F . i i Importantly, almost all practical hybrid approxima- 1 X J;w W K;w W tions can be re-expressed as functionals of the first- EHx[fφg] = F (ijjij)φ + F (ijjji)φ ; (6) 2 ij ij 0 ij order reduced density matrix (1-RDM), γ1(r; r ) = P ∗ 0 iσ2occ 'iσ(r)'iσ(r ). This includes all Hartree-Fock- W R 0 involving electron-repulsion integrals (ijjkl)φ = drdr like expressions, S[f'g] := hΦjT^ + W^ SjΦi, where Φ is 0 ∗ 0 ∗ 0 1 60 ^ W (r; r )φi (r)φj(r )φk(r)φl (r ) with W = jr−r0j . The a Slater determinant formed from f'g, and WS is a w coefficients fi 2 [0; 2] are average occupancies of orbitals (modified) interaction term involving interactions of form J;w 0 0 0 < WS(jr − r j) ≤ W (jr − r j), which can be a simple i across the ensemble. The pair-coefficients, Fij and fraction W = αW of the full interaction term, or a more F K;w, encode all information about the non-interacting S ij complex range-separated 8,61 expression.12,62 Then, ensemble and are ensemble-specific. They are more com- plicated to obtain { some examples are provided in ap- Z 1 2 0 pendix A. Importantly, these pair-coefficients are not 0 S[γ1] ≡ − 2 r γ1(r; r )jr=r dr necessarily (scaled) products of f w, i.e., generally F J 6= i ij Z K 0 0 1 0 2 drdr0 −2Fij 6= fjfj. Recent work has used the fluctuation- + WS(r; r ) n(r)n(r ) − 2 γ1(r; r ) 2 : (9) dissipation theorem to show that separated EH and Ex also yield functionals of similar form to Eq. (6).58 For pedagogical simplicity we shall not consider the separated We restrict ourselves to this broad and popular class of terms further. hybrids, though many results derived below are general. Just as in KS theory, the EKS orbitals, φi[n; w](r), are Our first step toward EGKS theory is to recognise that eigen-solutions of the KS equation, Eq. (9) can also accept more general 1-RDMs, ^ ft + v + vHxc[n; w]gφi[n; w] =i[n; w]φi[n; w] : (7) 0 X ∗ 0 γ1(r; r ) = fi'i (r)'i(r ) ; (10) δEHx[n;w] δEc[n;w] i Here, vHxc[n; w] = δn + δn is the multiplica- tive Hxc effective potential of EKS theory, and v is the ex- in which f 2 [0; 2] can be non-integer [note, n(r) = ternal potential.36 The effective potential, v := v +v , i s Hxc γ (r; r)]. Such 1-RDMs appear naturally in ensemble applied to the electrons thus depends both on the overall 1 theories, e.g., the 1-RDM for a lithium cation with a ensemble density, n, and the set of weights, w.
Recommended publications
  • Singlet/Triplet State Anti/Aromaticity of Cyclopentadienylcation: Sensitivity to Substituent Effect
    Article Singlet/Triplet State Anti/Aromaticity of CyclopentadienylCation: Sensitivity to Substituent Effect Milovan Stojanovi´c 1, Jovana Aleksi´c 1 and Marija Baranac-Stojanovi´c 2,* 1 Institute of Chemistry, Technology and Metallurgy, Center for Chemistry, University of Belgrade, Njegoševa 12, P.O. Box 173, 11000 Belgrade, Serbia; [email protected] (M.S.); [email protected] (J.A.) 2 Faculty of Chemistry, University of Belgrade, Studentski trg 12-16, P.O. Box 158, 11000 Belgrade, Serbia * Correspondence: [email protected]; Tel.: +381-11-3336741 Abstract: It is well known that singlet state aromaticity is quite insensitive to substituent effects, in the case of monosubstitution. In this work, we use density functional theory (DFT) calculations to examine the sensitivity of triplet state aromaticity to substituent effects. For this purpose, we chose the singlet state antiaromatic cyclopentadienyl cation, antiaromaticity of which reverses to triplet state aromaticity, conforming to Baird’s rule. The extent of (anti)aromaticity was evaluated by using structural (HOMA), magnetic (NICS), energetic (ISE), and electronic (EDDBp) criteria. We find that the extent of triplet state aromaticity of monosubstituted cyclopentadienyl cations is weaker than the singlet state aromaticity of benzene and is, thus, slightly more sensitive to substituent effects. As an addition to the existing literature data, we also discuss substituent effects on singlet state antiaromaticity of cyclopentadienyl cation. Citation: Stojanovi´c,M.; Aleksi´c,J.; Baranac-Stojanovi´c,M. Keywords: antiaromaticity; aromaticity; singlet state; triplet state; cyclopentadienyl cation; substituent Singlet/Triplet State effect Anti/Aromaticity of CyclopentadienylCation: Sensitivity to Substituent Effect.
    [Show full text]
  • Colour Structures in Supersymmetric QCD
    Colour structures in Supersymmetric QCD Jorrit van den Boogaard August 23, 2010 In this thesis we will explain a novel method for calculating the colour struc- tures of Quantum ChromoDynamical processes, or QCD prosesses for short, with the eventual goal to apply such knowledge to Supersymmetric theory. QCD is the theory of the strong interactions, among quarks and gluons. From a Supersymmetric viewpoint, something we will explain in the first chapter, it is interesting to know as much as possible about QCD theory, for it will also tell us something about this Supersymmetric theory. Colour structures in particular, which give the strength of the strong force, can tell us something about the interactions with the Supersymmetric sector. There are a multitude of different ways to calculate these colour structures. The one described here will make use of diagrammatic methods and is therefore somewhat more user friendly. We therefore hope it will see more use in future calculations regarding these processes. We do warn the reader for the level of knowledge required to fully grasp the subject. Although this thesis is written with students with a bachelor degree in mind, one cannot escape touching somewhat dificult matterial. We have tried to explain things as much as possible in a way which should be clear to the reader. Sometimes however we saw it neccesary to skip some of the finer details of the subject. If one wants to learn more about these hiatuses we suggest looking at the material in the bibliography. For convenience we will use some common conventions in this thesis.
    [Show full text]
  • Deuterium Isotope Effect in the Radiative Triplet Decay of Heavy Atom Substituted Aromatic Molecules
    Deuterium Isotope Effect in the Radiative Triplet Decay of Heavy Atom Substituted Aromatic Molecules J. Friedrich, J. Vogel, W. Windhager, and F. Dörr Institut für Physikalische und Theoretische Chemie der Technischen Universität München, D-8000 München 2, Germany (Z. Naturforsch. 31a, 61-70 [1976] ; received December 6, 1975) We studied the effect of deuteration on the radiative decay of the triplet sublevels of naph- thalene and some halogenated derivatives. We found that the influence of deuteration is much more pronounced in the heavy atom substituted than in the parent hydrocarbons. The strongest change upon deuteration is in the radiative decay of the out-of-plane polarized spin state Tx. These findings are consistently related to a second order Herzberg-Teller (HT) spin-orbit coupling. Though we found only a small influence of deuteration on the total radiative rate in naphthalene, a signi- ficantly larger effect is observed in the rate of the 00-transition of the phosphorescence. This result is discussed in terms of a change of the overlap integral of the vibrational groundstates of Tx and S0 upon deuteration. I. Introduction presence of a halogene tends to destroy the selective spin-orbit coupling of the individual triplet sub- Deuterium (d) substitution has proved to be a levels, which is originally present in the parent powerful tool in the investigation of the decay hydrocarbon. If the interpretation via higher order mechanisms of excited states The change in life- HT-coupling is correct, then a heavy atom must time on deuteration provides information on the have a great influence on the radiative d-isotope electronic relaxation processes in large molecules.
    [Show full text]
  • Introduction to Unconventional Superconductivity Manfred Sigrist
    Introduction to Unconventional Superconductivity Manfred Sigrist Theoretische Physik, ETH-Hönggerberg, 8093 Zürich, Switzerland Abstract. This lecture gives a basic introduction into some aspects of the unconventionalsupercon- ductivity. First we analyze the conditions to realized unconventional superconductivity in strongly correlated electron systems. Then an introduction of the generalized BCS theory is given and sev- eral key properties of unconventional pairing states are discussed. The phenomenological treatment based on the Ginzburg-Landau formulations provides a view on unconventional superconductivity based on the conceptof symmetry breaking.Finally some aspects of two examples will be discussed: high-temperature superconductivity and spin-triplet superconductivity in Sr2RuO4. Keywords: Unconventional superconductivity, high-temperature superconductivity, Sr2RuO4 INTRODUCTION Superconductivity remains to be one of the most fascinating and intriguing phases of matter even nearly hundred years after its first observation. Owing to the breakthrough in 1957 by Bardeen, Cooper and Schrieffer we understand superconductivity as a conden- sate of electron pairs, so-called Cooper pairs, which form due to an attractive interaction among electrons. In the superconducting materials known until the mid-seventies this interaction is mediated by electron-phonon coupling which gises rise to Cooper pairs in the most symmetric form, i.e. vanishing relative orbital angular momentum and spin sin- glet configuration (nowadays called s-wave pairing). After the introduction of the BCS concept, also studies of alternative pairing forms started. Early on Anderson and Morel [1] as well as Balian and Werthamer [2] investigated superconducting phases which later would be identified as the A- and the B-phase of superfluid 3He [3]. In contrast to the s-wave superconductors the A- and B-phase are characterized by Cooper pairs with an- gular momentum 1 and spin-triplet configuration.
    [Show full text]
  • Family Symmetry Grand Unified Theories
    Family symmetry Grand Unified Theories January 13, 2019 Abstract The mixing amongst the fermions of the Standard Model is a feature that has up to yet not been understood. It has been suggested that an A4 symmetry can be responsible for the specific mixing within the lepton sector. This thesis investigates the embedding of this A4 symmetry in a continuous family symmetry, SU(3)F . This continuous embedding allows for the construction of GUTs which include a family symmetry. Constructing and outlining these GUTs is the ultimate objective of this thesis. Author: Jelle Thole (s2183110) University of Groningen [email protected] Supervisor: Prof. Dr. D. Boer 2nd examinator: Prof. Dr. E. Pallante Contents 1 Introduction 3 2 The Standard Model of particle physics 4 2.1 Structure of the Standard model . .4 2.2 Standard Model Higgs . .6 2.2.1 Breaking U(1) . .7 2.2.2 Breaking SU(2) . .8 2.2.3 Electroweak symmetry breaking . .8 2.2.4 The Yukawa mechanism . 10 2.3 Quark masses . 11 2.3.1 The CKM Matrix . 11 2.4 Lepton masses . 13 2.4.1 PMNS matrix . 15 3 Family Symmetries 16 3.1 An A4 family symmetry . 17 3.2 An A4 invariant Lagrangian . 18 3.3 Explaining neutrino mixing . 19 4 Group theoretical tools 22 4.1 Lie groups and algebras . 22 4.2 Roots . 23 4.3 Dynkin diagrams . 24 4.4 Symmetry breaking using Dynkin diagrams . 28 4.5 Method of highest weights . 30 4.6 Symmetry breaking in grand unified theories .
    [Show full text]
  • Probing the C60 Triplet State Coupling to Nuclear Spins Inside and Out
    Downloaded from rsta.royalsocietypublishing.org on November 29, 2013 Probing the C60 triplet state coupling to nuclear spins inside and out Vasileia Filidou, Salvatore Mamone, Stephanie Simmons, Steven D. Karlen, Harry L. Anderson, Christopher W. M. Kay, Alessandro Bagno, Federico Rastrelli, Yasujiro Murata, Koichi Komatsu, Xuegong Lei, Yongjun Li, Nicholas J. Turro, Malcolm H. Levitt and John J. L. Morton Phil. Trans. R. Soc. A 2013 371, 20120475, published 5 August 2013 References This article cites 48 articles, 4 of which can be accessed free http://rsta.royalsocietypublishing.org/content/371/1998/201204 75.full.html#ref-list-1 Subject collections Articles on similar topics can be found in the following collections quantum physics (83 articles) Receive free email alerts when new articles cite this article - sign up Email alerting service in the box at the top right-hand corner of the article or click here To subscribe to Phil. Trans. R. Soc. A go to: http://rsta.royalsocietypublishing.org/subscriptions Downloaded from rsta.royalsocietypublishing.org on November 29, 2013 Probing the C60 triplet state coupling to nuclear spins inside and out Vasileia Filidou1,2,SalvatoreMamone3, Stephanie rsta.royalsocietypublishing.org Simmons1,2,StevenD.Karlen4, Harry L. Anderson4, ChristopherW.M.Kay1,5, Alessandro Bagno6, Federico Rastrelli6, Yasujiro Murata7, Koichi Research Komatsu7,XuegongLei8,YongjunLi8, Nicholas J. Cite this article: Filidou V, Mamone S, 8,† 3 1,9 Simmons S, Karlen SD, Anderson HL, Kay Turro , Malcolm H. Levitt andJohnJ.L.Morton CWM, Bagno A, Rastrelli F, Murata Y, Komatsu 1London Centre for Nanotechnology, University College London, K,LeiX,LiY,TurroNJ,LevittMH,MortonJJL.
    [Show full text]
  • Pauli Exclusion Principle and Multimulti--Electronelectron Atoms
    Quantum Mechanics and Atomic Physics Lecture 21: Pauli Exclusion Principle and MultiMulti--electronelectron atoms http://www. physics. rutgers. edu/ugrad/361 Prof. Sean Oh Last time : Fine structure constant Electron mass : mc2 : ~0.5 MeV Bohifdhr energies: of order α2mc2 : ~10 eV Fine structure: of order α4mc2 : ~10--44 eV Lamb shift: of order α5mc2 : ~10--66 eV 4 2 Hyperfine splitting: of order (m/mp)α mc ::: ~10--66 eV Multi-Electron Atoms Atoms with 2 or more electrons have a new feature: Electrons are indistinguishable! ΨA(1) There is no way to tell them apart! Any measurable quantity (probability , Ψ (2) expectation value, etc.) must not B depend on which electron is labeled 1, 2, etc. S.E. for Multi-electron atoms Let’s consider two electrons in Helium with coordinates: The total Hamiltonian operator for this system is 1 2 So the Schrodinger equation is: S.E. for Multi-electron atoms The total potential Vtot has 3 contributions: 1. V between electron 1 and the nucleus 2. V between electron 2 and the nucleus 3. V between electron 1 and electron 2 For now, let’s consider only #1 and #2 So, Note that the potential function is the same for both electrons S.E. for Multi-electron atoms We get the usual separation of variables Each Ψ will depend on quantum numbers n, ll , mll , ms So, A and B stand for the particular sets of quantum numbers So, let’s call ΨA(1) eigenfunction for electron #1 and has the quantum numbers symbolized by A .
    [Show full text]
  • Entanglement
    Intermediate Quantum Mechanics Lecture 22 Notes (4/15/15) Entanglement The importance of entanglement • Forty years ago (when I was learning quantum mechanics) no \real" physicist (which doesn't necessarily mean the best) cared about issues of entanglement and questions of what quantum mechanics \means". Textbooks were, by and large, devoid of the subject. Nowadays almost every textbook devotes at least a chapter to the subject. What has changed? • I can think of the following reasons why entanglement is now considered impor- tant. There are likely others − Entanglement is closely related to quantum computing that is currently a rather hot topic. − There is a possibility (though very unlikely) that quantum mechanics is not a fundament theory but an “effective" theory reflecting something that is deeper. Thinking about what quantum mechanics \means" might then give a hint of what the deeper theory is. − A rather radical view taken by some people working in the area of quantum computing is that quantum mechanics is not actually a theory of physics but rather it is a theory of probability involving complex numbers. They would say that physicists just stumbled upon this theory in trying to find a solution to their physics problems. As I mentioned, this is a rather radical point of view. − We all enjoy talking about philosophical issues such as what is objective reality or does it even exist. Although most physicists think these discussion are best had over coffee or even better over a few beers. − Nevertheless, there are several important and interesting issues involving entanglement that we will discuss in the next four lectures.
    [Show full text]
  • Critical Comparison of Triplet State Mechanisms, Singlet State Mechanisms, and Photoionization Mechanisms Will Be Deferred Until a Later Time
    VOL. 49, 1963 BIOCHEMISTRY: G. W. ROBINSON 521 some way, but we have no reason to believe that the enzymes initially present in the egg have been subjected to a significantly reduced possibility of error. Various solutions seem possible; for example, selection in the growing embryo may be strong enough or the process of embryogenesis may demand such a high accuracy that its successful completion guarantees the necessary accuracy of protein synthesis. One can also conceive of special mechanisms of quality control; for example, special proteins might be synthesized which are converted by a certain class of errors into lethal polypeptides. This would guarantee that the frequency of this class of errors in viable cells is kept low. At present there is no evidence available which enables one to select among these possibilities. I wish to make it quite clear that I am not proposing here that the accumulation of protein transcription errors is "the mechanism of ageing." My object is the more modest one of pointing out one source of progressive deterioration of cells and cell lines. Since I am unable to estimate the time scale of this process, I can only suggest experiments which should show where, if anywhere, it contributes to the ageing process in higher organisms. I am indebted to Professor H. C. Longuet-Higgins and Dr. F. H. C. Crick for valuable criticisms of my original manuscript. 1 Szilard, L., these PROCEEDINGS, 45, 30 (1959); Maynard Smith, J., Proc. Roy. Soc., B157, 115 (1962). 2Perutz, M., Proteins and Nucleic Acids: Structure and Function (Amsterdam: Elsevier, 1962).
    [Show full text]
  • Chapter 6: Nucelon-Nucleon Interactions and Deuteron
    Chapter 6 Nucleon-Nucleon Interaction, Deuteron Protons and neutrons are the lowest-energy bound states of quarks and gluons. When we put two or more of these particles together, they interact, scatter and sometimes form bound states due to the strong interactions. If one is interested in the low-energy region where the nucleons hardly get excited internally, we can treat the nucleons as inert, structureless elementary particles, and we can understand many of the properties of the multi-nucleon systems by the nucleon-nucleon interactions. If the nucleons are non-relativistic, the interaction can be described by a potential. Since the fundamental theory governing the nucleon-nucleon interactions is QCD, the interactions shall be calculable from the physics of quarks and gluons. Nonetheless, the problem is quite complicated and only limited progress has been made from the first principles so far. Therefore, the approach we are going take is a phenomenological one: One first tries to extract the nucleon- nucleon interaction from the nucleon-nucleon scattering data or few nucleon properties, and then one tries to use these interactions to make predictions for the nuclear many-body system. In this and the following few chapters, we are going to discuss some of results of this approach. 6.1 Yukawa Interaction and One-Pion Exchange at Long Distance If we consider the nucleons as elementary particles, their strong interactions are known have a short range from α particle scattering experiments of Rutherford. In fact, the range of the interaction is 15 roughly the size of the atomic nuclei, namely on the order of a few fermis (fm) (1 fm = 10− m).
    [Show full text]
  • Two Electrons: Excited States
    5.61 Physical Chemistry Lecture #26 page 1 TWO ELECTRONS: EXCITED STATES In the last lecture, we learned that the independent particle model gives a reasonable description of the ground state energy of the Helium atom. Before moving on to talk about many­electron atoms, it is important to point out that we can describe many more properties of the system using the same type of approximation. By using the same independent particle prescription we can come up with wavefunctions for excited states and determine their energies, their dependence on electron spin, etc. by examining the wavefunctions themselves. That is to say, there is much we can determine from simply looking at Ψ without doing any significant computation. We will use the excited state 1s2s configuration of Helium as an example. For the ground state we had: r r Ψspace ( 1, 2) × Ψspin (σ1,σ2) 1 ⎛ ⎞ ⇒ ψ r ψ r ⎜α σ β σ − β σ α σ ⎟ 1s 1s ( 1) 1s ( 2 ) 2 ⎝ ( 1) ( 2 ) ( 1) ( 2 ) ⎠ In constructing excited states it is useful to extend the stick diagrams we have used before to describe electronic configurations. Then there are four different configurations we can come up with: r r Ψspace ( 1, 2) × Ψspin (σ1,σ2) 2s ? ? 1s 2s ? ? 1s 2s ? ? 1s 2s ? ? 1s 5.61 Physical Chemistry Lecture #26 page 2 Where the question marks indicate that we need to determine the space and spin wavefunctions that correspond to these stick diagrams. It is fairly easy to come up with a reasonable guess for each configuration. For example, in the first case we might write down a wavefunction like: r r r r ψ 1sα ( 1σ1)ψ 2sα ( 2σ 2 ) =ψ 1s ( 1)ψ 2s ( 2 )α (σ 1)α (σ 2 ) However, we immediately note that this wavefunction is not antisymmetric.
    [Show full text]
  • Introduction to Heisenberg Model
    Introduction to Heisenberg model Javier Junquera Most important reference followed in this lecture Magnetism in Condensed Matter Physics Stephen Blundell Oxford Master Series in Condensed Matter Physics Exchange interactions If relativistic effects are not considered, then the electric interaction between particles does not depend on their spins In the absence of a magnetic field, the Hamiltonian of a system of particles in electric interaction does not contain spin operators A look to the hamiltonian: a difficult interacting many-body system. Kinetic energy operator for the electrons Potential acting on the electrons due to the nuclei Electron-electron interaction Kinetic energy operator for the nuclei Nucleus-nucleus interaction Exchange interactions If relativistic effects are not considered, then the electric interaction between particles does not depend on their spins In the absence of a magnetic field, the Hamiltonian of a system of particles in electric interaction does not contain spin operators When it is applied to a wavefunction, it does not affect the spinoidal variables The wave function of a system of particles can be written as a product It depends only on the spin It depends only on the variables spatial coordinates of the particles The Schrödinger equation determines only , leaving arbitrary The indistinguisibility of the particles leads to symmetry conditions on the wave function Despite the fact that the electric interaction is spin-independent, there is a dependendency of the energy of the system with respect the total
    [Show full text]