Notes on Ewald summation techniques

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "Notes on Ewald summation techniques"

Transcription

1 February 3, 011 Notes on Ewald summation techniques Adapted from similar presentation in PHY 71 he total electrostatic potential energy of interaction between point charges {q i } at the positions {r i } is given by W = 1 4πε 0 i<j q i q j r i r j = 1 q i q j r i r j. (1) In this expression, the first form explicitly counts all pairs, while the second form counts all interactions and divides by to compensate for double counting. For a finite system of charges, this expression can be evaluated directly, however, for a large or infinite system, the expression 1 does not converge and numerical tricks must be used to evaluate the energy. In fact, for the infinite system, one has an infinite amount of charge and the energy of interaction is undefined. If the system is neutral, it is possible to define a meaningful interaction energy by use of an Ewald transformation. he basic idea of the Ewald approach is as follows. he error function erf(x) and its complement erfc(x) are defined as: erf(x) = x e t dt and erfc(x) = 1 erf(x) = e t dt. () π π 0 x Ewald noted that 1 r = erf( 1 ηr) + erfc( 1 ηr). (3) r r In this expression, the first term goes to a constant ( η/π) as r 0, but has a long range tail as r. he second term has a singular behavior as r 0, but vanishes exponentially as r. his is illustrated in the plot below where η was chosen to be η = 4.

2 10 9 erf(x)/x (1-erf(x))/x 1/x hus, Ewald s idea is to replace a single divergent summation with two convergent summations. he first summation has a convergent summation in the form of its Fourier transform and the second has a convergent direct summation. hus the calculation of the electrostatic energy would be evaluated using: W = 1 q i q j r i r j = 1 q i q j erf( 1 η ri r j ) + r i r j q i q j erfc( 1 η ri r j ). r i r j (4) For an appropriate choice of the parameter η, the second summation in Eq. 4 converges quickly and can be evaluated directly. he first term in the summation of Eq. 4 must be transformed into Fourier space. In order to described these summations explicitly, we assume that we have a periodic lattice so that every ion can be located by r i = τ α + a location τ α within a unit cell and a periodic translation vector, In this way, the summation becomes: 1 = N (5) ij αβ, where N denotes the number of unit cells in the system. Since we have a periodic system, N is infinite, but the energy per unit cell W/N is well defined. he other identity that we 1 Note that for any two lattice translations i and j, i j = k, where k is also a lattice translation.

3 must use is that a sum over lattice translations may be transformed into an equivalent sum over reciprocal lattice translations according to the identity: δ 3 (r ) = 1 e i r, (6) Ω where Ω denotes the unit cell volume. he proof of this relation is given in the appendix. he first term of Eq. 4 thus becomes q i q j erf( 1 η ri r j )) = N q α q β r i r j αβ erf( 1 η τα τ β + ) τ α τ β + α q α η, (7) π where the last term in Eq. 7 comes from subtracting out the self-interaction (i = j) term from the complete lattice sum. Using the short hand notation, τ αβ τ α τ β, the lattice sum can be evaluated: his becomes, 1 Ω erf( 1 η ταβ + ) τ αβ + = d 3 re i r erf( 1 η ταβ + r ) τ αβ + r d 3 r = 4π Ω δ 3 (r ) erf( 1 η ταβ + r ). (8) τ αβ + r 0 e i τ αβ e /η + 1 where the last term, which is infinite, comes from the = 0 contribution. 1 η 0 du, (9) If the last term of Eq. 9 cannot be eliminated, it is clear that the electrostatic energy is infinite. he term can be eliminated if and only if the system is neutral. hus, it is only meaningful to calculate the electrostatic energy of a neutral periodic system. If the actual system has a net charge of Q α q α, we could calculate a meaningful energy if we add to the actual charge density, a compensating uniform density charge density of Q/Ω. aking all of the terms into account, we find the final Ewald expression to be: u 3 W N = αβ q α q β 4π Ω e i τ αβ e /η 0 η π δ αβ + erfc( 1 η ταβ + ) 4πQ τ αβ + Ωη, (10) where the in the summation over lattice translations indicates that all self-interaction terms should be omitted. Appendix I comments on lattice vectors and reciprocal lattice vectors In this discussion, will assume we have a 3-dimensional periodic system. It can be easily generalized to 1- or - dimensional systems. In general, a translation vector can be described a linear combination of the three primitive translation vectors 1,, and 3 : = n n + n 3 3, (11)

4 where {n 1, n, n 3 } are integers. Note that the unit cell volume Ω can be expressed in terms of the primitive translation vectors according to: Ω = 1 ( 3 ). (1) he reciprocal lattice vectors can generally be written as a linear combination of the three primitive reciprocal lattice vectors 1,, and 3 : = m m + m 3 3, (13) where {m 1, m, m 3 } are integers. he primitive reciprocal lattice vectors are determined from the primitive translation vectors according to the identities: i j = πδ ij. (14) Note that the volume of the primitive reciprocal lattice is given by 1 ( 3 ) = (π)3 Ω. (15) Some examples of this are given below. Proof of Eq. 6 Consider the geometric series +M k= M e ik( 1 r) = sin ( (M + 1 ) 1 r ) sin( 1 r/). (16) he behavior of the right hand side of Eq. (16) is that it is small in magnitude very except when the denominator vanishes. his occurs whenever 1 r/ = n 1 π, where n 1 represents any integer. If we take the limit M, we find that the function represents the behavior of a sum of delta functions: sin ( (M + 1 lim ) 1 r ) = π δ( 1 (r n 1 1 ). (17) M sin( 1 r/) n 1 he summation over all lattice translations n 1 1 is due to the fact that sin( 1 r/) = 0 whenever r = n 1 1. Carrying out the geometric summations in the right hand side of Eq. 6 for all three reciprocal lattice vectors and taking the limit as in Eq. 17, e i r = (π) 3 δ( 1 (r n 1 1 ))δ( (r n ))δ( 3 (r n 3 3 )). (18) n 1,n,n 3 Finally, the right hand side of Eq. 18 can be simplified by eliminating the reciprocal lattice vectors from the δ functions: e i r = which is consistent with Eq. 6. (π) 3 1 ( 3 ) δ 3 (r ) = Ωδ 3 (r ), (19)

5 Appendix II examples In these examples, denote the length of the unit cell by a. CsCl structure here are two kinds of sites τ Cs = 0 and τ Cl = a (ˆx + ŷ + ẑ). 1 = aˆx = aŷ 3 = aẑ. (0) 1 = π a ˆx = π a ŷ 3 = π a ẑ. (1) In these terms, Eq. 10 can be broken into two summations: W N = q ( 1 + ), () where and 1 4π Ω 0 erfc( η 0 (1 ei τ Cl )e /η η π (3) erfc( η τ Cl + ). (4) τ Cl + In this expression, the sum over α and β has a total of 4 contributions pairs of identical contributions. For Cs-Cs or Cl-Cl interactions, τ αβ = 0 and q α = q β resulting in repulsive contributions. For Cs-Cl or Cl-Cs interations, τ αβ = τ Cl and q α = q β resulting in attractive contributions. his expression can be evaluated using Maple, which gives the result W N = q ( 4.071). (5) a NaCl structure here are two kinds of sites τ Na = 0 and τ Cl = a ˆx. he primitive lattice has the translation vectors: and the reciprocal vectors: 1 = a (ˆx + ŷ) = a (ŷ + ẑ) 3 = a (ˆx + ẑ), (6) 1 = π a (ˆx + ŷ ẑ) = π a ( ˆx + ŷ + ẑ) 3 = π (ˆx ŷ + ẑ). (7) a

6 Note that there are many materials that have this structure, including MgO in which Mg has a charge of +e and O has a charge of -e. Zinc-blende structure his structure has the same bravais lattice as NaCl, but the basis vectors are τ Zn = 0 and τ S = a (ˆx + ŷ + ẑ). In this case Zn has charge of +e and S has a charge of -e. 4

7 Figure 1: Diagram of CsCl structure, indicating translation vectors and atomic site vectors τ. τ

E12 UNDERSTANDING CRYSTAL STRUCTURES

E12 UNDERSTANDING CRYSTAL STRUCTURES E1 UNDERSTANDING CRYSTAL STRUCTURES 1 Introduction In this experiment, the structures of many elements and compounds are rationalized using simple packing models. The pre-work revises and extends the material

More information

Reciprocal Lattice. Let A(r) be some function featuring discrete translation symmetry:

Reciprocal Lattice. Let A(r) be some function featuring discrete translation symmetry: Reciprocal Lattice From Quantum Mechanics we know that symmetries have kinematic implications. Each symmetry predetermines its specific quantum numbers and leads to certain constraints (conservation laws/selection

More information

Lattice energy of ionic solids

Lattice energy of ionic solids 1 Lattice energy of ionic solids Interatomic Forces Solids are aggregates of atoms, ions or molecules. The bonding between these particles may be understood in terms of forces that play between them. Attractive

More information

Keble College - Hilary 2012 Section VI: Condensed matter physics Tutorial 2 - Lattices and scattering

Keble College - Hilary 2012 Section VI: Condensed matter physics Tutorial 2 - Lattices and scattering Tomi Johnson Keble College - Hilary 2012 Section VI: Condensed matter physics Tutorial 2 - Lattices and scattering Please leave your work in the Clarendon laboratory s J pigeon hole by 5pm on Monday of

More information

Atomic Arrangement. Primer in Materials Spring

Atomic Arrangement. Primer in Materials Spring Atomic Arrangement Primer in Materials Spring 2017 30.4.2017 1 Levels of atomic arrangements No order In gases, for example the atoms have no order, they are randomly distributed filling the volume to

More information

Taylor Series and Series Convergence (Online)

Taylor Series and Series Convergence (Online) 7in 0in Felder c02_online.te V3 - February 9, 205 9:5 A.M. Page CHAPTER 2 Taylor Series and Series Convergence (Online) 2.8 Asymptotic Epansions In introductory calculus classes the statement this series

More information

Ceramic Bonding. CaF 2 : large SiC: small

Ceramic Bonding. CaF 2 : large SiC: small Recall ceramic bonding: - Mixed ionic and covalent. - % ionic character ( f ) increases with difference in electronegativity Large vs small ionic bond character: Ceramic Bonding CaF : large SiC: small

More information

Introduction to Series and Sequences Math 121 Calculus II Spring 2015

Introduction to Series and Sequences Math 121 Calculus II Spring 2015 Introduction to Series and Sequences Math Calculus II Spring 05 The goal. The main purpose of our study of series and sequences is to understand power series. A power series is like a polynomial of infinite

More information

Bravais Lattice + Basis = Crystal Structure

Bravais Lattice + Basis = Crystal Structure Bravais Lattice + Basis = Crystal Structure A crystal structure is obtained when identical copies of a basis are located at all of the points of a Bravais lattice. Consider the structure of Cr, a I-cubic

More information

Geometry of Crystal Lattice

Geometry of Crystal Lattice 0 Geometry of Crystal Lattice 0.1 Translational Symmetry The crystalline state of substances is different from other states (gaseous, liquid, amorphous) in that the atoms are in an ordered and symmetrical

More information

1 Solutions in cylindrical coordinates: Bessel functions

1 Solutions in cylindrical coordinates: Bessel functions 1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates

More information

... 3, , = a (1) 3 3 a 2 = a (2) The reciprocal lattice vectors are defined by the condition a b = 2πδ ij, which gives

... 3, , = a (1) 3 3 a 2 = a (2) The reciprocal lattice vectors are defined by the condition a b = 2πδ ij, which gives PHZ646: Fall 013 Problem set # 4: Crystal Structure due Monday, 10/14 at the time of the class Instructor: D. L. Maslov maslov@phys.ufl.edu 39-0513 Rm. 114 Office hours: TR 3 pm-4 pm Please help your instructor

More information

CHEM6085: Density Functional Theory

CHEM6085: Density Functional Theory Lecture 11 CHEM6085: Density Functional Theory DFT for periodic crystalline solids C.-K. Skylaris 1 Electron in a one-dimensional periodic box (in atomic units) Schrödinger equation Energy eigenvalues

More information

In particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with

In particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient

More information

Nearly Free Electron Gas model - I

Nearly Free Electron Gas model - I Nearly Free Electron Gas model - I Contents 1 Free electron gas model summary 1 2 Electron effective mass 3 2.1 FEG model for sodium...................... 4 3 Nearly free electron model 5 3.1 Primitive

More information

FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017

FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II November 5, 207 Prof. Alan Guth FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 207 A few items below are marked

More information

This ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0

This ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0 Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x

More information

Phys 412 Solid State Physics. Lecturer: Réka Albert

Phys 412 Solid State Physics. Lecturer: Réka Albert Phys 412 Solid State Physics Lecturer: Réka Albert What is a solid? A material that keeps its shape Can be deformed by stress Returns to original shape if it is not strained too much Solid structure

More information

Atomic Arrangement. Primer Materials For Science Teaching Spring

Atomic Arrangement. Primer Materials For Science Teaching Spring Atomic Arrangement Primer Materials For Science Teaching Spring 2016 31.3.2015 Levels of atomic arrangements No order In gases, for example the atoms have no order, they are randomly distributed filling

More information

Failures and successes of the free electron model

Failures and successes of the free electron model Failures and successes of the free electron model Temperature dependence on the electric conduc1vity j e = E = ne2 m electron gas electron- phonon coupling Drude s model predicts that σ is independent

More information

UNIT I SOLID STATE PHYSICS

UNIT I SOLID STATE PHYSICS UNIT I SOLID STATE PHYSICS CHAPTER 1 CRYSTAL STRUCTURE 1.1 INTRODUCTION When two atoms are brought together, two kinds of forces: attraction and repulsion come into play. The force of attraction increases

More information

Computer simulation methods (2) Dr. Vania Calandrini

Computer simulation methods (2) Dr. Vania Calandrini Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing

More information

Generalized sources (Sect. 6.5). The Dirac delta generalized function. Definition Consider the sequence of functions for n 1, Remarks:

Generalized sources (Sect. 6.5). The Dirac delta generalized function. Definition Consider the sequence of functions for n 1, Remarks: Generalized sources (Sect. 6.5). The Dirac delta generalized function. Definition Consider the sequence of functions for n, d n, t < δ n (t) = n, t 3 d3 d n, t > n. d t The Dirac delta generalized function

More information

Math review. Math review

Math review. Math review Math review 1 Math review 3 1 series approximations 3 Taylor s Theorem 3 Binomial approximation 3 sin(x), for x in radians and x close to zero 4 cos(x), for x in radians and x close to zero 5 2 some geometry

More information

Convergence of sequences and series

Convergence of sequences and series Convergence of sequences and series A sequence f is a map from N the positive integers to a set. We often write the map outputs as f n rather than f(n). Often we just list the outputs in order and leave

More information

Metallic and Ionic Structures and Bonding

Metallic and Ionic Structures and Bonding Metallic and Ionic Structures and Bonding Ionic compounds are formed between elements having an electronegativity difference of about 2.0 or greater. Simple ionic compounds are characterized by high melting

More information

Notes for Expansions/Series and Differential Equations

Notes for Expansions/Series and Differential Equations Notes for Expansions/Series and Differential Equations In the last discussion, we considered perturbation methods for constructing solutions/roots of algebraic equations. Three types of problems were illustrated

More information

Limits and Continuity

Limits and Continuity Chapter Limits and Continuity. Limits of Sequences.. The Concept of Limit and Its Properties A sequence { } is an ordered infinite list x,x,...,,... The n-th term of the sequence is, and n is the index

More information

Assignment 16 Assigned Weds Oct 11

Assignment 16 Assigned Weds Oct 11 Assignment 6 Assigned Weds Oct Section 8, Problem 3 a, a 3, a 3 5, a 4 7 Section 8, Problem 4 a, a 3, a 3, a 4 3 Section 8, Problem 9 a, a, a 3, a 4 4, a 5 8, a 6 6, a 7 3, a 8 64, a 9 8, a 0 56 Section

More information

Crystal Structure Determination II

Crystal Structure Determination II Crystal Structure Determination II Dr. Falak Sher Pakistan Institute of Engineering and Applied Sciences 09/10/2010 Diffraction Intensities The integrated intensity, I (hkl) (peak area) of each powder

More information

Review of Linear Time-Invariant Network Analysis

Review of Linear Time-Invariant Network Analysis D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x

More information

Physics 116C The Distribution of the Sum of Random Variables

Physics 116C The Distribution of the Sum of Random Variables Physics 116C The Distribution of the Sum of Random Variables Peter Young (Dated: December 2, 2013) Consider a random variable with a probability distribution P(x). The distribution is normalized, i.e.

More information

Differential of the Exponential Map

Differential of the Exponential Map Differential of the Exponential Map Ethan Eade May 20, 207 Introduction This document computes ɛ0 log x + ɛ x ɛ where and log are the onential mapping and its inverse in a Lie group, and x and ɛ are elements

More information

Homework #1 Solutions version 3 6 problems

Homework #1 Solutions version 3 6 problems Homework #1 Solutions version 6 problems Exercises 1. Charges A and lie 5 cm apart, on a line. Charge A is q A = + µc, and charge is q = mc. (e careful of the metrix prefixes!) (a) efore doing any calculation,

More information

8.1 Sequences. Example: A sequence is a function f(n) whose domain is a subset of the integers. Notation: *Note: n = 0 vs. n = 1.

8.1 Sequences. Example: A sequence is a function f(n) whose domain is a subset of the integers. Notation: *Note: n = 0 vs. n = 1. 8. Sequences Example: A sequence is a function f(n) whose domain is a subset of the integers. Notation: *Note: n = 0 vs. n = Examples: 6. Find a formula for the general term a n of the sequence, assuming

More information

Electron levels in a periodic potential: general remarks. Daniele Toffoli January 11, / 46

Electron levels in a periodic potential: general remarks. Daniele Toffoli January 11, / 46 Electron levels in a periodic potential: general remarks Daniele Toffoli January 11, 2017 1 / 46 Outline 1 Mathematical tools 2 The periodic potential 3 Bloch s theorem and Born-von Karman boundary conditions

More information

Taylor and Laurent Series

Taylor and Laurent Series Chapter 4 Taylor and Laurent Series 4.. Taylor Series 4... Taylor Series for Holomorphic Functions. In Real Analysis, the Taylor series of a given function f : R R is given by: f (x + f (x (x x + f (x

More information

Chemical Bonding -- Lewis Theory (Chapter 9)

Chemical Bonding -- Lewis Theory (Chapter 9) Chemical Bonding -- Lewis Theory (Chapter 9) Ionic Bonding 1. Ionic Bond Electrostatic attraction of positive (cation) and negative (anion) ions Neutral Atoms e - transfer (IE and EA) cation + anion Ionic

More information

MTS 105 LECTURE 5: SEQUENCE AND SERIES 1.0 SEQUENCE

MTS 105 LECTURE 5: SEQUENCE AND SERIES 1.0 SEQUENCE MTS 105 LECTURE 5: SEQUENCE AND SERIES 1.0 SEQUENCE A sequence is an endless succession of numbers placed in a certain order so that there is a first number, a second and so on. Consider, for example,

More information

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )

More information

Landau s Fermi Liquid Theory

Landau s Fermi Liquid Theory Thors Hans Hansson Stockholm University Outline 1 Fermi Liquids Why, What, and How? Why Fermi liquids? What is a Fermi liquids? Fermi Liquids How? 2 Landau s Phenomenological Approach The free Fermi gas

More information

Sample Spaces, Random Variables

Sample Spaces, Random Variables Sample Spaces, Random Variables Moulinath Banerjee University of Michigan August 3, 22 Probabilities In talking about probabilities, the fundamental object is Ω, the sample space. (elements) in Ω are denoted

More information

Chapter 6 INORGANIC THERMODYNAMICS. Exercises

Chapter 6 INORGANIC THERMODYNAMICS. Exercises Chapter 6 INORGANIC THERMODYNAMICS Exercises 6. (a) A reaction that occurs without external help, or a reaction for which G is negative. (b) A measure of disorder. (c) The enthalpy change when a mole of

More information

Lecture Note on Crystal structures Masatsugu Sei Suzuki and Itsuko S. Suzuki Department of Physics, SUNY at Binghamton (Date: February 03, 2012)

Lecture Note on Crystal structures Masatsugu Sei Suzuki and Itsuko S. Suzuki Department of Physics, SUNY at Binghamton (Date: February 03, 2012) Lecture Note on Crystal structures Masatsugu Sei Suzuki and Itsuko S. Suzuki Department of Physics, SUNY at Binghamton (Date: February 03, 2012) This is a part of lecture note on solid state physics (Phys.472/572)

More information

Problem Set # 1 Solutions CHAPTERS 2 & 3 ATOMIC STRUCTURE AND INTERATOMIC BONDING and THE STRUCTURE OF CRYSTALLINE SOLIDS

Problem Set # 1 Solutions CHAPTERS 2 & 3 ATOMIC STRUCTURE AND INTERATOMIC BONDING and THE STRUCTURE OF CRYSTALLINE SOLIDS Problem Set # Solutions CHAPTERS & ATOMIC STRUCTURE AND INTERATOMIC BONDING and THE STRUCTURE OF CRYSTALLINE SOLIDS Assigned:.7(a),.9,.,.6,.8,.7(a),.7,.9,. (graded problems indicated in bold).7 (a) The

More information

12.1 Arithmetic Progression Geometric Progression General things about sequences

12.1 Arithmetic Progression Geometric Progression General things about sequences ENGR11 Engineering Mathematics Lecture Notes SMS, Victoria University of Wellington Week Five. 1.1 Arithmetic Progression An arithmetic progression is a sequence where each term is found by adding a fixed

More information

Multipole moments. Dipole moment. The second moment µ is more commonly called the dipole moment, of the charge. distribution and is a vector

Multipole moments. Dipole moment. The second moment µ is more commonly called the dipole moment, of the charge. distribution and is a vector Dipole moment Multipole moments The second moment µ is more commonly called the dipole moment, of the charge distribution and is a vector µ = µ x ˆx + µ y ŷ + µ z ẑ where the α component is given by µ

More information

Crystal planes. Neutrons: magnetic moment - interacts with magnetic materials or nuclei of non-magnetic materials. (in Å)

Crystal planes. Neutrons: magnetic moment - interacts with magnetic materials or nuclei of non-magnetic materials. (in Å) Crystallography: neutron, electron, and X-ray scattering from periodic lattice, scattering of waves by periodic structures, Miller indices, reciprocal space, Ewald construction. Diffraction: Specular,

More information

PHY Fall HW6 Solutions

PHY Fall HW6 Solutions PHY249 - Fall 216 - HW6 Solutions Allen Majewski Department Of Physics, University of Florida 21 Museum Rd. Gainesville, FL 32611 October 11, 216 These are solutions to Halliday, Resnick, Walker Chapter

More information

Green s functions for planarly layered media

Green s functions for planarly layered media Green s functions for planarly layered media Massachusetts Institute of Technology 6.635 lecture notes Introduction: Green s functions The Green s functions is the solution of the wave equation for a point

More information

SECTION 9.2: ARITHMETIC SEQUENCES and PARTIAL SUMS

SECTION 9.2: ARITHMETIC SEQUENCES and PARTIAL SUMS (Chapter 9: Discrete Math) 9.11 SECTION 9.2: ARITHMETIC SEQUENCES and PARTIAL SUMS PART A: WHAT IS AN ARITHMETIC SEQUENCE? The following appears to be an example of an arithmetic (stress on the me ) sequence:

More information

Name: Date: Blk: Examine your periodic table to answer these questions and fill-in-the-blanks. Use drawings to support your answers where needed:

Name: Date: Blk: Examine your periodic table to answer these questions and fill-in-the-blanks. Use drawings to support your answers where needed: Name: Date: Blk: NOTES: BONDING Examine your periodic table to answer these questions and fill-in-the-blanks. Use drawings to support your answers where needed: I. IONIC BONDING Ionic bond: formed by the

More information

Physics 221B Spring 2018 Notes 34 The Photoelectric Effect

Physics 221B Spring 2018 Notes 34 The Photoelectric Effect Copyright c 2018 by Robert G. Littlejohn Physics 221B Spring 2018 Notes 34 The Photoelectric Effect 1. Introduction In these notes we consider the ejection of an atomic electron by an incident photon,

More information

Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0

Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0 Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can

More information

5. Orthogonal matrices

5. Orthogonal matrices L Vandenberghe EE133A (Spring 2017) 5 Orthogonal matrices matrices with orthonormal columns orthogonal matrices tall matrices with orthonormal columns complex matrices with orthonormal columns 5-1 Orthonormal

More information

Ionic Bonding - Electrostatic Interactions and Polarization

Ionic Bonding - Electrostatic Interactions and Polarization Ionic Bonding - Electrostatic Interactions and Polarization Chemistry 754 Solid State Chemistry Dr. Patrick Woodward Lecture #13 Born-Haber Cycle for NaCl It is energetically unfavorable for Na metal and

More information

Two-Body Problem. Central Potential. 1D Motion

Two-Body Problem. Central Potential. 1D Motion Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of

More information

Handout 7 Reciprocal Space

Handout 7 Reciprocal Space Handout 7 Reciprocal Space Useful concepts for the analysis of diffraction data http://homepages.utoledo.edu/clind/ Concepts versus reality Reflection from lattice planes is just a concept that helps us

More information

where a is the lattice constant of the triangular Bravais lattice. reciprocal space is spanned by

where a is the lattice constant of the triangular Bravais lattice. reciprocal space is spanned by Contents 5 Topological States of Matter 1 5.1 Intro.......................................... 1 5.2 Integer Quantum Hall Effect..................... 1 5.3 Graphene......................................

More information

Continuous Time Signal Analysis: the Fourier Transform. Lathi Chapter 4

Continuous Time Signal Analysis: the Fourier Transform. Lathi Chapter 4 Continuous Time Signal Analysis: the Fourier Transform Lathi Chapter 4 Topics Aperiodic signal representation by the Fourier integral (CTFT) Continuous-time Fourier transform Transforms of some useful

More information

Experimental Determination of Crystal Structure

Experimental Determination of Crystal Structure Experimental Determination of Crystal Structure Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, U.S.A. PHYS 624: Introduction to Solid State Physics http://www.physics.udel.edu/~bnikolic/teaching/phys624/phys624.html

More information

Calculating Fractal Dimension of Attracting Sets of the Lorenz System

Calculating Fractal Dimension of Attracting Sets of the Lorenz System Dynamics at the Horsetooth Volume 6, 2014. Calculating Fractal Dimension of Attracting Sets of the Lorenz System Jamie Department of Mathematics Colorado State University Report submitted to Prof. P. Shipman

More information

1 Multiplicity of the ideal gas

1 Multiplicity of the ideal gas Reading assignment. Schroeder, section.6. 1 Multiplicity of the ideal gas Our evaluation of the numbers of microstates corresponding to each macrostate of the two-state paramagnet and the Einstein model

More information

PHYS4210 Electromagnetic Theory Quiz 1 Feb 2010

PHYS4210 Electromagnetic Theory Quiz 1 Feb 2010 PHYS4210 Electromagnetic Theory Quiz 1 Feb 2010 1. An electric dipole is formed from two charges ±q separated by a distance b. For large distances r b from the dipole, the electric potential falls like

More information

Wave diffraction and the reciprocal lattice

Wave diffraction and the reciprocal lattice Wave diffraction and the reciprocal lattice Dept of Phys M.C. Chang Braggs theory of diffraction Reciprocal lattice von Laue s theory of diffraction Braggs view of the diffraction (1912, father and son)

More information

Graphene and Planar Dirac Equation

Graphene and Planar Dirac Equation Graphene and Planar Dirac Equation Marina de la Torre Mayado 2016 Marina de la Torre Mayado Graphene and Planar Dirac Equation June 2016 1 / 48 Outline 1 Introduction 2 The Dirac Model Tight-binding model

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises

More information

8.323 Relativistic Quantum Field Theory I

8.323 Relativistic Quantum Field Theory I MIT OpenCourseWare http://ocw.mit.edu 8.33 Relativistic Quantum Field Theory I Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Prof.Alan Guth

More information

Algebra 2 (2006) Correlation of the ALEKS Course Algebra 2 to the California Content Standards for Algebra 2

Algebra 2 (2006) Correlation of the ALEKS Course Algebra 2 to the California Content Standards for Algebra 2 Algebra 2 (2006) Correlation of the ALEKS Course Algebra 2 to the California Content Standards for Algebra 2 Algebra II - This discipline complements and expands the mathematical content and concepts of

More information

A number that can be written as, where p and q are integers and q Number.

A number that can be written as, where p and q are integers and q Number. RATIONAL NUMBERS 1.1 Definition of Rational Numbers: What are rational numbers? A number that can be written as, where p and q are integers and q Number. 0, is known as Rational Example:, 12, -18 etc.

More information

TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1

TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1 COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2

More information

Infinite series, improper integrals, and Taylor series

Infinite series, improper integrals, and Taylor series Chapter 2 Infinite series, improper integrals, and Taylor series 2. Introduction to series In studying calculus, we have explored a variety of functions. Among the most basic are polynomials, i.e. functions

More information

Ergodic Theory. Constantine Caramanis. May 6, 1999

Ergodic Theory. Constantine Caramanis. May 6, 1999 Ergodic Theory Constantine Caramanis ay 6, 1999 1 Introduction Ergodic theory involves the study of transformations on measure spaces. Interchanging the words measurable function and probability density

More information

Chapter 1. Introduction to Electrostatics

Chapter 1. Introduction to Electrostatics Chapter. Introduction to Electrostatics. Electric charge, Coulomb s Law, and Electric field Electric charge Fundamental and characteristic property of the elementary particles There are two and only two

More information

The Feynman-Kac formula

The Feynman-Kac formula The Feynman-Kac formula William G. Faris February, 24 The Wiener process (Brownian motion) Consider the Hilbert space L 2 (R d ) and the self-adjoint operator H = σ2, () 2 where is the Laplace operator.

More information

Absence of Luttinger s Theorem

Absence of Luttinger s Theorem Absence of Luttinger s Theorem Kiaran dave Charlie Kane n = µ N(ω)dω µ IR UV n = µ N(ω)dω µ IR UV can n be deduced entirely from the IR(lowenergy) scale? Luttinger s theorem n =2 k Θ(G(k,ω = 0)) Luttinger

More information

Translation Symmetry, Space Groups, Bloch functions, Fermi energy

Translation Symmetry, Space Groups, Bloch functions, Fermi energy Translation Symmetry, Space Groups, Bloch functions, Fermi energy Roberto Orlando and Silvia Casassa Università degli Studi di Torino July 20, 2015 School Ab initio Modelling of Solids (UniTo) Symmetry

More information

Complex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,

Complex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )

More information

Section 10 Metals: Electron Dynamics and Fermi Surfaces

Section 10 Metals: Electron Dynamics and Fermi Surfaces Electron dynamics Section 10 Metals: Electron Dynamics and Fermi Surfaces The next important subject we address is electron dynamics in metals. Our consideration will be based on a semiclassical model.

More information

Complex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,

Complex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )

More information

Physics 2212 GH Quiz #4 Solutions Spring 2016

Physics 2212 GH Quiz #4 Solutions Spring 2016 Physics 2212 GH Quiz #4 Solutions Spring 2016 I. (18 points) A bar (mass m, length L) is connected to two frictionless vertical conducting rails with loops of wire, in the presence of a uniform magnetic

More information

Solutions to Complex Analysis Prelims Ben Strasser

Solutions to Complex Analysis Prelims Ben Strasser Solutions to Complex Analysis Prelims Ben Strasser In preparation for the complex analysis prelim, I typed up solutions to some old exams. This document includes complete solutions to both exams in 23,

More information

Fermi surfaces which produce large transverse magnetoresistance. Abstract

Fermi surfaces which produce large transverse magnetoresistance. Abstract Fermi surfaces which produce large transverse magnetoresistance Stephen Hicks University of Florida, Department of Physics (Dated: August 1, ) Abstract The Boltzmann equation is used with elastic s-wave

More information

The Matrix Representation of a Three-Dimensional Rotation Revisited

The Matrix Representation of a Three-Dimensional Rotation Revisited Physics 116A Winter 2010 The Matrix Representation of a Three-Dimensional Rotation Revisited In a handout entitled The Matrix Representation of a Three-Dimensional Rotation, I provided a derivation of

More information

Class 8: Numbers Exercise 3B

Class 8: Numbers Exercise 3B Class : Numbers Exercise B 1. Compare the following pairs of rational numbers: 1 1 i First take the LCM of. LCM = 96 Therefore: 1 = 96 Hence we see that < 6 96 96 1 1 1 1 = 6 96 1 or we can say that

More information

1 Infinite-Dimensional Vector Spaces

1 Infinite-Dimensional Vector Spaces Theoretical Physics Notes 4: Linear Operators In this installment of the notes, we move from linear operators in a finitedimensional vector space (which can be represented as matrices) to linear operators

More information

Chapter 7: Chemical Bonding and Molecular Structure

Chapter 7: Chemical Bonding and Molecular Structure Chapter 7: Chemical Bonding and Molecular Structure Ionic Bond Covalent Bond Electronegativity and Bond Polarity Lewis Structures Orbital Overlap Hybrid Orbitals The Shapes of Molecules (VSEPR Model) Molecular

More information

Fractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials

Fractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials Applied Mathematical Sciences, Vol. 5, 211, no. 45, 227-2216 Fractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials Z. Avazzadeh, B. Shafiee and G. B. Loghmani Department

More information

Appendix A. Sequences and series. A.1 Sequences. Definition A.1 A sequence is a function N R.

Appendix A. Sequences and series. A.1 Sequences. Definition A.1 A sequence is a function N R. Appendix A Sequences and series This course has for prerequisite a course (or two) of calculus. The purpose of this appendix is to review basic definitions and facts concerning sequences and series, which

More information

Stochastic Processes

Stochastic Processes Stochastic Processes A very simple introduction Péter Medvegyev 2009, January Medvegyev (CEU) Stochastic Processes 2009, January 1 / 54 Summary from measure theory De nition (X, A) is a measurable space

More information

Algebraic aspects of some Riordan arrays related to binary words avoiding a pattern

Algebraic aspects of some Riordan arrays related to binary words avoiding a pattern Algebraic aspects of some Riordan arrays related to binary words avoiding a pattern D. Merlini and R. Sprugnoli Dipartimento di Sistemi e Informatica viale Morgagni 65, 5034, Firenze, Italia merlini@dsi.unifi.it,

More information

A Simply Regularized Derivation of the Casimir Force

A Simply Regularized Derivation of the Casimir Force EJTP 3, No. 13 (2006) 121 126 Electronic Journal of Theoretical Physics A Simply Regularized Derivation of the Casimir Force H. Razmi Habibollah Razmi, 37185-359, Qom, I. R. Iran Department of Physics,

More information

arxiv:physics/ v1 [physics.comp-ph] 19 Dec 2005

arxiv:physics/ v1 [physics.comp-ph] 19 Dec 2005 Analytical Hartree-Fock gradients with respect to the cell parameter: systems periodic in one and two dimensions K. Doll Institut für Mathematische Physik, TU Braunschweig, Mendelssohnstraße 3, D-38106

More information

We are going to discuss what it means for a sequence to converge in three stages: First, we define what it means for a sequence to converge to zero

We are going to discuss what it means for a sequence to converge in three stages: First, we define what it means for a sequence to converge to zero Chapter Limits of Sequences Calculus Student: lim s n = 0 means the s n are getting closer and closer to zero but never gets there. Instructor: ARGHHHHH! Exercise. Think of a better response for the instructor.

More information

Review of Fourier Transform

Review of Fourier Transform Review of Fourier Transform Fourier series works for periodic signals only. What s about aperiodic signals? This is very large & important class of signals Aperiodic signal can be considered as periodic

More information

10.1 Sequences. A sequence is an ordered list of numbers: a 1, a 2, a 3,..., a n, a n+1,... Each of the numbers is called a term of the sequence.

10.1 Sequences. A sequence is an ordered list of numbers: a 1, a 2, a 3,..., a n, a n+1,... Each of the numbers is called a term of the sequence. 10.1 Sequences A sequence is an ordered list of numbers: a 1, a 2, a 3,..., a n, a n+1,... Each of the numbers is called a term of the sequence. Notation: A sequence {a 1, a 2, a 3,...} can be denoted

More information

The infinite series is written using sigma notation as: lim u k. lim. better yet, we can say if the

The infinite series is written using sigma notation as: lim u k. lim. better yet, we can say if the Divergence and Integral Test With the previous content, we used the idea of forming a closed form for the n th partial sum and taking its limit to determine the SUM of the series (if it exists). *** It

More information

Saxon Calculus Scope and Sequence

Saxon Calculus Scope and Sequence hmhco.com Saxon Calculus Scope and Sequence Foundations Real Numbers Identify the subsets of the real numbers Identify the order properties of the real numbers Identify the properties of the real number

More information

Experiment 7: Understanding Crystal Structures

Experiment 7: Understanding Crystal Structures Experiment 7: Understanding Crystal Structures To do well in this laboratory experiment you need to be familiar with the concepts of lattice, crystal structure, unit cell, coordination number, the different

More information

Surface Structure and Morphology 2D Crystallography

Surface Structure and Morphology 2D Crystallography Surface Structure and Morphology 2D Crystallography Selvage (or selvedge (it. cimosa)): Region in the solid in the vicinity of the mathematical surface Surface = Substrate (3D periodicity) + Selvage (few

More information