Lecture VIII : The pseudopotential
|
|
- Warren Riley
- 5 years ago
- Views:
Transcription
1 Lecture VIII : The pseudopotentia I. KOHN-SHAM PROBLEM FOR AN ISOLATED ATOM For a one-eectron atom, the Couombic potentia, V ( r) = V (r) = Z/r is sphericay symmetric. The soutions may then be spit into a radia and an anguar part ψ nm ( r) = ψ n (r)y m (θ, φ) = r 1 φ n (r)y m (θ, φ) (1) The ast equaity in Eq. 1 is given in order to simpify the sphericay-symmetric Schrödinger equation and Y m (θ, φ) are normaized spherica harmonics. The sphericay-symmetric Schrödinger equation for principe quantum number n can be reduced to the radia equation 1 d 2 [ ] ( + 1) 2 dr 2 φ n + 2r 2 + V ext (r) ǫ n, φ n, = 0. (2) In the Kohn-Sham approach to the many-partice system, the form of the singe-partice equations are identica to the above radia Schrödinger equation, with an effective potentia, V eff repacing the pure Couomb potentia. The substitution can be done directy provided V eff is aso sphericay symmetric. Let s see if that is the case : Density : n( r) = occ n,,m Externa potentia : Ionic Couomb potentia occ ψ n,,m (r) 2 = (2 + 1) ψ n, (r) 2 = n(r) (3) n, V ext ( r) = Z/r = V ext (r) (4) The Hartree potentia V H ( r) = = n( r ) r r d r n(r ) [(r r cosθ sin φ) 2 + (r cosθ sinθ) 2 + (r sin θ) 2 ] 1/2 r 2 dr d(cos θ)dφ = V H (r) (5) The exchange-correation potentia V xc ( r) = ǫ xc [n(r)] + n(r) dǫ xc dn [n(r)] = V xc(r) (6) Tota effective potentia V eff (r) = V ext (r) + V H (r) + V xc (r) (7) Thus the independent-partice Kohn-Sham states may be cassified by the anguar quantum numbers L =, m. The one-partice equations are then anaogous to the Schrödinger equation for the one-eectron atom, 1 d 2 [ ] ( + 1) 2 dr 2 φ n + 2r 2 + V eff (r) ǫ n, φ n, = 0. (8)
2 2 II. THE PSEUDOPOTENTIAL IDEA The concept of a pseudopotentia is reated to repacing the effects of the core eectrons with an effective potentia. The pseudopotentia generation procedure starts with the soution of the atomic probem using the Kohn-Sham approach. Once the KS orbitas are obtained, we make an arbitrary distinction between vaence and core states. The core states are assumed to change very itte due to changes in the environment so their effect is repaced by a mode potentia derived in the atomic configuration and it is assumed to be transferabe. The vaence states are seen to osciate rapidy cose to the core regions. With the introduction of the new potentia, the vaence states are made smoother. Let s now work out the pseudopotentia transformation in its most genera framework : Assume that we have a Hamitonian Ĥ, core states { χ n } and core eigenvaues {E n }. We are ooking at a singe vaence state ψ Let s repace the vaence state by the smoother φ and expand the remaining portion in terms of core states : core ψ = φ + a n χ n (9) Next, we take the inner product of Eq. 9 with one of the core states. Because the vaence state has to be orthogona to the core states, we have n core χ n ψ = χ n φ + a n χ m χ n = 0 (10) n } {{ } a m We now write the right-hand side of Eq. 10 in terms of the pseudofunction, φ ψ = φ n χ n φ χ n (11) Appying the Hamitonian onto the expression in Eq. 11, we get core Ĥ φ + (E E n ) χ n χ n φ = E φ (12) n As a resut, the smooth, pseudowavefunction satisfies an effective equation with the same eigenenergy of the rea vaence wavefunction. In the case of isoated atoms, the indices n corresponds to the combined index nm incuding the principa quantum number n, the anguar momentum quantum number and the magnetic quantum number m. However, since the potentia does not incude any terms (such as L S couping) that ift the m degenaracy, we incude in what foows ony n and. The reason why the -degenaracy is broken instead wi become cear ater on. Eq. 12 may then be written as an eigenvaue equation for the smooth pseudo wavefunction as (Ĥ + ˆVn ) φ = E φ (13) where the extra potentia ˆV n depends on due to spherica symmetry and its effect is ocaized to the core. Since E > E n by definition, it is a repusive potentia and it partiay cances the effect of the attractive Couombic potentia. The resuting potentia is then a much weaker one than the origina potentia. This shows that the eigenstates of this new potentia are smoother. III. NORM-CONSERVING PSEUDOPOTENTIALS In order to ensure the optimum smoothness and transferabiity Hamann, Schüter and Chiang (Phys Rev Lett (1979)) came up with four criteria that the pseudopotentias are expected to obey in the reference configuration (i.e. the atomic configuration for which the pseudopotentia is generated.) : 1. A-eectron and pseudo eigenvaues agree for the reference configuration. Ĥ ψn AE = ǫ n ψn AE (14) ) (Ĥ + ˆVn ψn PS = ǫ n ψn PS (15) (16)
3 3 2. AE and PS wavefunctions agree beyond a certain cutoff, r c. 3. Rea and pseudo norm squares integrated from 0 to R for a R < r c agree. ψ AE n (r) = ψ PS n (r) for r r c (17) R 0 φ AE n 2 r 2 dr = R 0 φ PS n 2 r 2 dr (18) 4. The ogarithmic derivatives and the energy derivative of the ogarithmic derivative agree for a R < r c. [ (rφ AE (r)) 2 d ] [ d de dr nφae (r) = (rφ PS (r)) 2 d ] d de dr nφps (r) R R (19) Items 1 and 2 are obvious. Item 3 is necessary so that by Gauss theorem, outside the core region where r < r c, the eectrostatic potentia is identica for the rea and pseudo wavefunctions. Finay, 4 is necessary in order to improve transferabiity, reproducing phase shifts of the rea potentia in the scattering sense. When we go from an atom to a different environment such as a soid or a moecue, the eigenvaues change. A pseudopotentia satisfying this item reproduces these changes to inear order. Using the Friede sum rue, one can show that item 4 is impied by item 3. In order to generate a norm-conserving pseudopotentia, we must foow the steps beow : 1. Sove the a-eectron atomic system. 2. Decide which of the states are to be taken as core and which are to be taken as vaence. 3. Generate the pseudopotentias, V n (r), for the vaence states. There are severa methods for doing that. We describe here the method due to Kerker. In this method, we assume a parametrized anaytica function for the core region. The anaytic form chosen by Kerker is φ PS n (r) = ep(r), where p(r) is a fourth-order poynomia matching the true soution and its various derivatives at r c. 4. Once we have a form for the pseudowavefunction, invert the Schro dinger equation, yieding [ d ( + 1) 2 dr φ PS V n,tota (r) = ǫ 2r 2 2 n (r) ] φ PS n (r) (20) Keep in mind that we are deaing with φ(r), which is defined by ψ(r) = rφ(r). 5. The potentia V,tota (r) found by inversion of the Schrödinger equation is not ony the ionic potentia that we want, but incudes aso contributions from the Hartree and exchange-correation terms. We must therefore subtract these two contributions. V n (r) = V,tota (r) VHxc PS (r) (21) where the Hartree and xc contributions are computed for the pseudowavefunctions. This operation is caed unscreening. There is no unique way to determine the V n s. There are two opposing considerations : 1. Good transferabiity sma r c. 2. Large r c smoother pseudopotentias. A good pseudopotentia is one that strikes a baance between these two contraints. Athough the generation procedure is simiar, the form of the pseudopotentia may be assumed differenty depending on considerations such as ease of computation, accuracy etc. In the foowing subsections, we sha see two different ways of forming norm-conserving pseudopotentias. Keep in mind that these pseudopotentias are formed within the atomic configuration and then exported to the soid with the hope of decent transferabiity. An actua exampe of the construction of the pseudopotentia can be seen in Fig. 1 for the radia parts of 2S and 2P orbitas of the eement B. The agreement of the wavefunctions after the cutoff radius r c is evident.
4 rea 2S for B pseudo 2S for B rc= rea 2P for B pseudo 2P for B rc= FIG. 1: Semiog pot of the rea and pseudo wavefunctions for 2S (a) and 2P (b) orbitas of B. A. Semioca pseudopotentias Since a mode pseudopotentia repaces the potentia of a nuceus, it is sphericay symmetric and each anguar momentum may be treated separatey. This eads to the nonoca, i.e. -dependent mode pseudopotentias V n (r). Thus, the genera form of a pseudopotentia may be written in the semioca form as ˆV SL = m Y m V n (r) Y m. (22) Here the anguar and the radia parts are separated due to spherica symmetry and the Y m s are there because any function that depends upon anges ony may be written as an expansion. This form is caed semioca because it is nonoca ony in the anguar coordinates and it s oca in the radia coordinate. Athough conceptuay simpe, semioca pseudopotentias create a severe probem in the panewave basis representation due to the form of the nonoca part. This can be seen by considering a particuar matrix eement of the pseudopotentia and using a ot of straighforward but tedious agebra q ˆV SL q = 1 e i q r Y Ω m(θ, φ)v n (r)y m (θ, φ)e i q r r 2 dω dω dr (23) m where r = (r, θ, φ ). Notice the fact that we cacuate a singe integra over the radia component r but two integras over the anges because the nonoca part of the pseudopotentia is in the form of projectors instead of a reguar function as for the radia coordinate. In order to make use of their orthonormaity, we expand the panewaves in the integra in Eq. 23 in terms of spherica harmonics e i q r = 4π m i j (qr)y m (ˆq)Y m(ˆr) (24) where j (qr) are spherica Besse functions and ˆq and ˆr denote the anges associated with the vectors q and r. Substituting Eq. 24 into Eq. 23 we have q ˆV NL q = (4π)2 ( i) i j (qr)j (q r)y m Ω m (ˆq )V n (r)r 2 dr m m m Y m (ˆr )Y m (ˆr )dω Ym m (ˆr)dΩ } {{ }} {{ } δ,mm δ,mm (25) where we have used the orthonormaity of spherica harmonics as mentioned before. When the deta functions are imposed the two sums over m and m are gone and the ony indices that remain are and m. Next in order to dea with the remaining spherica harmonics, we utiize the summation theorem for spherica harmonics. m Y m (ˆq)Y m (ˆq ) = π P (cos θ q q ) (26)
5 where P is the generaized Legendre poynomia and θ q q is the ange between the two wavevectors q and q. Putting everything together the expectation vaue in Eq. 25 becomes q ˆV NL q = 4π (2 + 1) j (qr)j (q r)p (cos θ q q )V (r)r 2 dr (27) Ω As seen in Eq. 27, the semioca pseudopotentia mixes panewaves with different q vectors. This means that if we have N pw panewaves in the basis that we choose we have to evauate Npw 2 such integras as in Eq. 27, which brings a arge computationa oad. In their semina work Keinman and Byander (Physica Review Letters, 48, 1425) proposed a different scheme for composing the anguar momentum-dependent part of the pseudopotentia such that the number of integras to be cacuated coud be reduced from Npw 2 to just N pw. This form is caed the fuy nonoca form and we sha see how it is formed in the next section. 5 B. The Keinman-Byander Transformation The difficuty with the semioca pseudopotentias above is that due to the ocaity in the radia coordinate there is a singe radia variabe for both panewaves whie the anguar part is nonoca and the integras reated to each anguar coordinate may be evauated separatey. The fuy nonoca formuation incudes the radia part in the projectors thereby separating the integra reated to each projector. The starting point for the transformation is the usua semioca form expained above where the oca radia part is given by V. We then find the pseudowavefunctions in the atomic configuration, φ PS corresponding to this pseudopotentia. Using these wavefunctions, we construct a new pseudopotentia who is composed of a oca part, V L, which is a competey arbitrary function that has the ong-range behavior of the rea potentia and a nonoca part V NL, which is the remaining part of the pseudopotentia. After the subtraction of the oca part, the nonoca part is eft behind ocaized in the core region. The nonoca part is then given by ˆV NL = ˆV NL = m ˆV SL φ PS m φps SL m ˆV φ PS m ˆV φ PS m (28) This particuar form is chosen such that its action on a pseudowavefunction yieds the same resut as the action of the semioca form ˆV NL φ PS m In order to see this, et s appy this operator on the pseudowavefunction ˆV NL φ PS m = SL ˆV φ PS SL = ˆV φ PS m (29) m φps m φ PS m ˆV SL SL ˆV φ PS φ PS m m = SL ˆV φ PS m (30) In the panewave representation, an arbitrary matrix eement of the Keinman-Byanderized pseudopotentia is q ˆV NL q = q ˆV φ PS m ˆV φ PS m q φ PS m ˆV φ PS m. (31) The two integras in Eq. 31 do not share any common integra variabes and may thus be performed separatey. One coud then compute N PW such integras, i.e. q ˆV φ PS m in advance and cacuate the appropriate product q ˆV φ PS m ˆV φ PS m q for each q, q pair. IV. ULTRASOFT PSEUDOPOTENTIALS In 1990, Vanderbit (see Physica Review B, 41, 7892) proposed a new and radica method for generating pseudopotentias by reaxing the norm-conservation constraint. Pseudopotentias generated in this way (due to their softness) require a much smaer panewave cutoff and thus a much smaer number of panewaves. For this reason they are usuay caed utrasoft pseudopotentias. First et s see how utrasoft pseudopotentia is integrated into the Kohn-Sham equations and then see how it is generated. The tota energy in the presence of utrasoft pseudopotentias is written as
6 6 E e = i φ i V NL φ i i d rd r n( r)n( r ) r r + E xc [n] + d rv oc ( r)n( r) (32) where the Hartree, exchange-correation, oca and kinetic terms have the usua form. The nonoca part of the pseudopotentia instead is given in terms of a new set of projection operators as V NL = D (0) nm βi n βi m (33) where the superscript I indicates that the projector is centered on atom I. The norm-conserving condition is reaxed by introducing a generaized orthonormaity condition through the overap matrix S where φ i S φ j = δ ij (34) S = 1 + nm,i q nm β I n βi m (35) with q nm = d rq nm ( r). Because the norm-conserving condition is now reaxed, the eectronic density needs to be augmented in order to preserve the correct number of eectrons. The new density is then n( r) = i φ i ( r) 2 + Q I nm ( r) φ i βn I βi m φ i (36) nm,i In order to verify that the density as given in Eq. 36 integrates to the correct number of eectrons, consider n( r)d r = φ i ( r) 2 + d rq I nm( r) φ i βn β I m φ I i i nm,i = ij ij δ ij φ i φ j + nm,i,j δ ij q nm φ i βn β I m φ I j δ ij φ i S φ j = δ ij δ ij = δ ii = 1 = N v (37) ij i i where N v stands for the number of vaence eectrons in the system. The new form of the density, athough it gives us the fexibiity to create softer pseudopotentias, inevitaby makes the Kohn-Sham equations more compex. The functiona derivative of the density with respect to the orbitas is δn( r ) δφ i ( r) = φ i( r)δ( r r ) + Q I nm ( r )βn I ( r) βi m φ i. (38) nm,i Utiizing Eq. 38 and Eq. 32 we arrive at the modified Kohn-Sham equations δe e δφ = d r δe e δn( r [ ) i δn( r ) δφ i ( r) = 1 ( ) ] V eff + D nm (0) + d r V eff ( r )Q nm ( r ) βn β I m I φ i. (39) where V eff is the reguar potentia V eff = V H + V oc + V xc. (40) As easiy seen in Eq. 39, unike norm-conserving pseudopotentias, the coefficients of the projectors in the nonoca part of the pseudopotentia gets updated in a sef-consistent manner at every iteration. This introduces a sight overhead to the cacuation but the gain in the panewave expansion surpasses the oss thus introduced. Now, we take a step back and see the justification and generation of such a radica impementation. The generation of utrasoft potentias start out in the same way as norm-conserving pseudopotentias. For a given species, the atomic
7 Kohn-Sham system is soved sef-consistenty, resuting in the screened a-eectron potentia, V AE. Then for each anguar momentum channe, a few energy vaues, ǫ τ are chosen (τ=number of such energies, usuay max 3). For the V AE aready obtained, the Schrödinger equation is found [T + V AE ] ψ n (r) = ǫ n ψ n (r) (41) where n is a composite index hoding for the anguar momentum quantum number, m for the projection and τ. In the genera case, ψ n is not going to be normaizabe but we use the bra-ket notation anyway. Next, we choose three cutoff radii : 1. A smooth oca (-independent) potentia is created which matches V AE (r) after the radius r oc c. 2. A second cutoff radius, r c, is chosen for each anguar momentum channe and a pseudowavefunction, φ n is constructed which matches the rea wavefunction ψ n at r c. 3. Finay a diagnostic radius R is chosen which is sighty arger than the maximum of the maximum of a the r c and r oc C. Then new orbitas are formed that satisfy which are zero beyond R. Next, an auxiiary matrix of inner products and the projectors necessary for the definition of the nonoca part of the potentia χ n = (ǫ n T V oc ) φ n (42) B nm = φ n χ m (43) β n = m (B 1 ) mn χ m (44) 7 are defined. Whie determining the pseudowavefunctions mentioned above, norm-conservation has not been imposed. We thus define a set of generaized augmentation functions and the reated augmentation charges Q nm ( r) = ψ n( r)ψ m ( r) φ n( r)φ m ( r) (45) q nm = ψ n ψ m R φ n φ m R. (46) where the subscript R denotes an integration within the cutoff radius, R. With a these descriptions, we can verify that the pseudowavefunctions φ n obey the equation ( ˆT + V oc + ) ( D nm β n β m φ i = ǫ i 1 + ) q nm β n β m φ i (47) nm nm where D nm = B nm + ǫ m q nm. In order to prove the centra equation Eq. 47, et s ook at ( ) B nm β n β m φ i = B nm (B 1 ) n χ (B 1 km ) χ k phi i }{{} nm nm k Bik = m χ n (B 1 ) n B nm (Bik (B 1 ) km ) = m δ m δ im χ = χ i = (ǫ i ˆT V oc ) φ i (48) Sustituting Eq. 48 into Eq. 47 competes the proof.
8 Thus the fundamenta equation of the pseudopotentia theory has been satisfied, that is pseudowavefunctions that satisfy the same equation as the rea wavefunctions have been found. These pseudopotentias, however, may be proven to satisfy, instead of being orthonorma in the usua sense, satisfy a generaized orthonormaity equation through an overap matrix S = 1 + nm q nm β n β m (49) 8 Once the D nm are obtained, the oca part and the nonoca coefficients of the bare pseudopotentia are obtained through a descreening procedure V ion oc = V oc V H V xc (50) D nm (0) = D nm d r V oc ( r )n( r ) (51) Whie generating pseudopotentias, we don t enforce norm conservation. This aows us to increase the cutoff radii to arger vaues. This convenient fact however brings about the necessity of defining a separate panewave cutoff for the density because as seen in the expression in Eq. 36, the density now has a smooth part and a sharp part. This sharp part needs to be expressed using panewaves that have high frequencies, i.e. a pane wave cutoff that is typicay severa times the panewave cutoff used to express the wavefunctions.
DFT in practice : Part II. Ersen Mete
pseudopotentials Department of Physics Balıkesir University, Balıkesir - Turkey August 13, 2009 - NanoDFT 09, İzmir Institute of Technology, İzmir Outline Pseudopotentials Basic Ideas Norm-conserving pseudopotentials
More informationMore Scattering: the Partial Wave Expansion
More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction
More informationLegendre Polynomials - Lecture 8
Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha
More information221B Lecture Notes Notes on Spherical Bessel Functions
Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationJoel Broida UCSD Fall 2009 Phys 130B QM II. Homework Set 2 DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS.
Joe Broida UCSD Fa 009 Phys 30B QM II Homework Set DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS. You may need to use one or more of these: Y 0 0 = 4π Y 0 = 3 4π cos Y
More informationPhysics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions
Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p
More informationMidterm 2 Review. Drew Rollins
Midterm 2 Review Drew Roins 1 Centra Potentias and Spherica Coordinates 1.1 separation of variabes Soving centra force probems in physics (physica systems described by two objects with a force between
More informationIntroduction to LMTO method
1 Introduction to MTO method 24 February 2011; V172 P.Ravindran, FME-course on Ab initio Modeing of soar ce Materias 24 February 2011 Introduction to MTO method Ab initio Eectronic Structure Cacuations
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More information$, (2.1) n="# #. (2.2)
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationIntroduction to DFT and Density Functionals. by Michel Côté Université de Montréal Département de physique
Introduction to DFT and Density Functionas by Miche Côté Université de Montréa Département de physique Eampes Carbazoe moecue Inside of diamant Réf: Jean-François Brière http://www.phys.umontrea.ca/ ~miche_cote/images_scientifiques/
More informationPhysics 235 Chapter 8. Chapter 8 Central-Force Motion
Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion
More informationSEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l
Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed
More informationLecture 8 February 18, 2010
Sources of Eectromagnetic Fieds Lecture 8 February 18, 2010 We now start to discuss radiation in free space. We wi reorder the materia of Chapter 9, bringing sections 6 7 up front. We wi aso cover some
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationPhysics 116C Helmholtz s and Laplace s Equations in Spherical Polar Coordinates: Spherical Harmonics and Spherical Bessel Functions
Physics 116C Hemhotz s an Lapace s Equations in Spherica Poar Coorinates: Spherica Harmonics an Spherica Besse Functions Peter Young Date: October 28, 2013) I. HELMHOLTZ S EQUATION As iscusse in cass,
More informationThis shoud be true for any, so the equations are equivaent to the Schrodingerequation (H E) : () Another important property of the (2) functiona is, t
I. Variationa method Introduction Usefu for the determination of energies and wavefunctions of different atomic states. H the time-independent Hamitonian of the system E n eigenenergies Ψ n - eigenfunctions
More informationhole h vs. e configurations: l l for N > 2 l + 1 J = H as example of localization, delocalization, tunneling ikx k
Infinite 1-D Lattice CTDL, pages 1156-1168 37-1 LAST TIME: ( ) ( ) + N + 1 N hoe h vs. e configurations: for N > + 1 e rij unchanged ζ( NLS) ζ( NLS) [ ζn unchanged ] Hund s 3rd Rue (Lowest L - S term of
More informationJackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.
More informationSection 6: Magnetostatics
agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The
More informationAn approximate method for solving the inverse scattering problem with fixed-energy data
J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999
More informationPHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased
PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization
More informationPhysics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27
Physics 55 Fa 7 Homework Assignment #5 Soutions Textook proems: Ch. 3: 3.3, 3.7, 3.6, 3.7 3.3 Sove for the potentia in Proem 3., using the appropriate Green function otained in the text, and verify that
More informationAgenda Administrative Matters Atomic Physics Molecules
Fromm Institute for Lifeong Learning University of San Francisco Modern Physics for Frommies IV The Universe - Sma to Large Lecture 3 Agenda Administrative Matters Atomic Physics Moecues Administrative
More informationHYDROGEN ATOM SELECTION RULES TRANSITION RATES
DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS
More informationThe Accuracy of the Pseudopotential Approximation within Density-Functional Theory
D. Porezag et a.: The Accuracy of the Pseudopotentia Approximation 219 phys. stat. so. (b) 217, 219 (2000) Subject cassification: 71.15.Hx; 71.15.Mb The Accuracy of the Pseudopotentia Approximation within
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationBohr s atomic model. 1 Ze 2 = mv2. n 2 Z
Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated
More informationQuantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18
Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms
More informationPhysics 506 Winter 2006 Homework Assignment #6 Solutions
Physics 506 Winter 006 Homework Assignment #6 Soutions Textbook probems: Ch. 10: 10., 10.3, 10.7, 10.10 10. Eectromagnetic radiation with eiptic poarization, described (in the notation of Section 7. by
More informationPhysics 505 Fall Homework Assignment #4 Solutions
Physics 505 Fa 2005 Homework Assignment #4 Soutions Textbook probems: Ch. 3: 3.4, 3.6, 3.9, 3.0 3.4 The surface of a hoow conducting sphere of inner radius a is divided into an even number of equa segments
More informationElectron-impact ionization of diatomic molecules using a configuration-average distorted-wave method
PHYSICAL REVIEW A 76, 12714 27 Eectron-impact ionization of diatomic moecues using a configuration-average distorted-wave method M. S. Pindzoa and F. Robicheaux Department of Physics, Auburn University,
More informationEffects of energy loss on interaction dynamics of energetic electrons with plasmas. C. K. Li and R. D. Petrasso. 1 November 2008
PSFC/JA-8-3 ffects of energy oss on interaction dynamics of energetic ctrons with pasmas C. K. Li and R. D. Petrasso November 8 Pasma Science and Fusion Center Massachusetts Institute of Technoogy Cambridge,
More informationLECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October, 202 Prof. Aan Guth LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS
More informationFFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection
FFTs in Graphics and Vision Spherica Convoution and Axia Symmetry Detection Outine Math Review Symmetry Genera Convoution Spherica Convoution Axia Symmetry Detection Math Review Symmetry: Given a unitary
More informationGaussian Curvature in a p-orbital, Hydrogen-like Atoms
Advanced Studies in Theoretica Physics Vo. 9, 015, no. 6, 81-85 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/astp.015.5115 Gaussian Curvature in a p-orbita, Hydrogen-ike Atoms Sandro-Jose Berrio-Guzman
More informationLECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October 7, 202 Prof. Aan Guth LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationApplied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation
22.101 Appied Nucear Physics (Fa 2006) Lecture 7 (10/2/06) Overview of Cross Section Cacuation References P. Roman, Advanced Quantum Theory (Addison-Wesey, Reading, 1965), Chap 3. A. Foderaro, The Eements
More informationDavid Eigen. MA112 Final Paper. May 10, 2002
David Eigen MA112 Fina Paper May 1, 22 The Schrodinger equation describes the position of an eectron as a wave. The wave function Ψ(t, x is interpreted as a probabiity density for the position of the eectron.
More informationNuclear Size and Density
Nucear Size and Density How does the imited range of the nucear force affect the size and density of the nucei? Assume a Vecro ba mode, each having radius r, voume V = 4/3π r 3. Then the voume of the entire
More informationPhysics 566: Quantum Optics Quantization of the Electromagnetic Field
Physics 566: Quantum Optics Quantization of the Eectromagnetic Fied Maxwe's Equations and Gauge invariance In ecture we earned how to quantize a one dimensiona scaar fied corresponding to vibrations on
More informationarxiv: v1 [hep-lat] 23 Nov 2017
arxiv:1711.08830v1 [hep-at] 23 Nov 2017 Tetraquark resonances computed with static attice QCD potentias and scattering theory Pedro Bicudo 1,, Marco Cardoso 1, Antje Peters 2, Martin Pfaumer 2, and Marc
More informationCS229 Lecture notes. Andrew Ng
CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view
More informationCluster modelling. Collisions. Stellar Dynamics & Structure of Galaxies handout #2. Just as a self-gravitating collection of objects.
Stear Dynamics & Structure of Gaaxies handout # Custer modeing Just as a sef-gravitating coection of objects. Coisions Do we have to worry about coisions? Gobuar custers ook densest, so obtain a rough
More informationParallel-Axis Theorem
Parae-Axis Theorem In the previous exampes, the axis of rotation coincided with the axis of symmetry of the object For an arbitrary axis, the paraeaxis theorem often simpifies cacuations The theorem states
More informationLecture 17 - The Secrets we have Swept Under the Rug
Lecture 17 - The Secrets we have Swept Under the Rug Today s ectures examines some of the uirky features of eectrostatics that we have negected up unti this point A Puzze... Let s go back to the basics
More informationRadiation Fields. Lecture 12
Radiation Fieds Lecture 12 1 Mutipoe expansion Separate Maxwe s equations into two sets of equations, each set separatey invoving either the eectric or the magnetic fied. After remova of the time dependence
More informationTRANSFORMATION OF REAL SPHERICAL HARMONICS UNDER ROTATIONS
Vo. 39 (008) ACTA PHYSICA POLONICA B No 8 TRANSFORMATION OF REAL SPHERICAL HARMONICS UNDER ROTATIONS Zbigniew Romanowski Interdiscipinary Centre for Materias Modeing Pawińskiego 5a, 0-106 Warsaw, Poand
More informationModule 22: Simple Harmonic Oscillation and Torque
Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque
More information1.2 Partial Wave Analysis
February, 205 Lecture X.2 Partia Wave Anaysis We have described scattering in terms of an incoming pane wave, a momentum eigenet, and and outgoing spherica wave, aso with definite momentum. We now consider
More informationLecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential
Lecture 6 Povh Krane Enge Wiiams Properties of -nuceon potentia 16.1 4.4 3.6 9.9 Meson Theory of Nucear potentia 4.5 3.11 9.10 I recommend Eisberg and Resnik notes as distributed Probems, Lecture 6 1 Consider
More informationSupporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers
Supporting Information for Suppressing Kein tunneing in graphene using a one-dimensiona array of ocaized scatterers Jamie D Was, and Danie Hadad Department of Chemistry, University of Miami, Cora Gabes,
More informationScattering of Particles by Potentials
Scattering of Partices by Potentias 2 Introduction The three prototypes of spectra of sef-adjoint operators, the discrete spectrum, with or without degeneracy, the continuous spectrum, andthe mixed spectrum,
More informationDECAY THEORY BEYOND THE GAMOW PICTURE
Dedicated to Academician Aureiu Sanduescu s 8 th Anniversary DECAY THEORY BEYOND THE GAMOW PICTURE D. S. DELION Horia Huubei Nationa Institute for Physics and Nucear Engineering, P.O. Box MG-6, Bucharest,
More informationASummaryofGaussianProcesses Coryn A.L. Bailer-Jones
ASummaryofGaussianProcesses Coryn A.L. Baier-Jones Cavendish Laboratory University of Cambridge caj@mrao.cam.ac.uk Introduction A genera prediction probem can be posed as foows. We consider that the variabe
More informationWhy use pseudo potentials?
Pseudo potentials Why use pseudo potentials? Reduction of basis set size effective speedup of calculation Reduction of number of electrons reduces the number of degrees of freedom For example in Pt: 10
More informationThree Dimensional systems. In three dimensions the time-independent Schrödinger Equation becomes
Three Dimensiona systems In three dimensions the time-independent Schrödinger Equation becomes ] [ m r + V r) ψr) = Eψr), r = x + x + x. The easiest case one can imagine in 3D is when the potentia energy
More informationXSAT of linear CNF formulas
XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open
More informationarxiv: v1 [math.ca] 6 Mar 2017
Indefinite Integras of Spherica Besse Functions MIT-CTP/487 arxiv:703.0648v [math.ca] 6 Mar 07 Joyon K. Boomfied,, Stephen H. P. Face,, and Zander Moss, Center for Theoretica Physics, Laboratory for Nucear
More information4 1-D Boundary Value Problems Heat Equation
4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x
More informationarxiv: v1 [cond-mat.mtrl-sci] 10 Mar 2015
Constructing optima oca pseudopotentias from first principes arxiv:153.2992v1 [cond-mat.mtr-sci] 1 Mar 215 Wenhui Mi, 1,2 Shoutao Zhang, 1,2 Yanming Ma, 1,2, and Maosheng Miao 3,2, 1 State Key Laboratory
More informationarxiv:nlin/ v2 [nlin.cd] 30 Jan 2006
expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische
More information14 Separation of Variables Method
14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt
More information11.1 One-dimensional Helmholtz Equation
Chapter Green s Functions. One-dimensiona Hemhotz Equation Suppose we have a string driven by an externa force, periodic with frequency ω. The differentia equation here f is some prescribed function) 2
More informationarxiv: v1 [nucl-th] 25 Nov 2011
Deveopment of a Cox-Thompson inverse scattering method to charged partices Tamás Pámai 1, Barnabás Apagyi 1 and Werner Scheid 1 Department of Theoretica Physics Budapest University of Technoogy and Economics,
More informationWeek 6 Lectures, Math 6451, Tanveer
Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n
More informationSupplemental Information
Suppementa Information A Singe-eve Tunne Mode to Account for Eectrica Transport through Singe Moecue- Sef-Assembed Monoayer-based Junctions by A. R. Garrigues,. Yuan,. Wang, E. R. Muccioo, D. Thompson,
More informationCopyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU
Copyright information to be inserted by the Pubishers Unspitting BGK-type Schemes for the Shaow Water Equations KUN XU Mathematics Department, Hong Kong University of Science and Technoogy, Cear Water
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),
More informationHigher dimensional PDEs and multidimensional eigenvalue problems
Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u
More informationMONTE CARLO SIMULATIONS
MONTE CARLO SIMULATIONS Current physics research 1) Theoretica 2) Experimenta 3) Computationa Monte Caro (MC) Method (1953) used to study 1) Discrete spin systems 2) Fuids 3) Poymers, membranes, soft matter
More informationVIII. Addition of Angular Momenta
VIII Addition of Anguar Momenta a Couped and Uncouped Bae When deaing with two different ource of anguar momentum, Ĵ and Ĵ, there are two obviou bae that one might chooe to work in The firt i caed the
More informationVolume 13, MAIN ARTICLES
Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More information17 Lecture 17: Recombination and Dark Matter Production
PYS 652: Astrophysics 88 17 Lecture 17: Recombination and Dark Matter Production New ideas pass through three periods: It can t be done. It probaby can be done, but it s not worth doing. I knew it was
More informationAtom-Atom Interactions in Ultracold Quantum Gases
Atom-Atom Interactions in Utracod Quantum Gases Caude Cohen-Tannoudji Lectures on Quantum Gases Institut Henri Poincaré, Paris, 5 Apri 7 Coège de France 1 Lecture 1 (5 Apri 7) Quantum description of eastic
More informationProblem set 6 The Perron Frobenius theorem.
Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator
More informationDiscrete Techniques. Chapter Introduction
Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, we as various
More informationA nodal collocation approximation for the multidimensional P L equations. 3D applications.
XXI Congreso de Ecuaciones Diferenciaes y Apicaciones XI Congreso de Matemática Apicada Ciudad Rea, 1-5 septiembre 9 (pp. 1 8) A noda coocation approximation for the mutidimensiona P L equations. 3D appications.
More informationDiscrete Techniques. Chapter Introduction
Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, as we as various
More informationVI.G Exact free energy of the Square Lattice Ising model
VI.G Exact free energy of the Square Lattice Ising mode As indicated in eq.(vi.35), the Ising partition function is reated to a sum S, over coections of paths on the attice. The aowed graphs for a square
More informationThe state vectors j, m transform in rotations like D(R) j, m = m j, m j, m D(R) j, m. m m (R) = j, m exp. where. d (j) m m (β) j, m exp ij )
Anguar momentum agebra It is easy to see that the operat J J x J x + J y J y + J z J z commutes with the operats J x, J y and J z, [J, J i ] 0 We choose the component J z and denote the common eigenstate
More informationIn Coulomb gauge, the vector potential is then given by
Physics 505 Fa 007 Homework Assignment #8 Soutions Textbook probems: Ch. 5: 5.13, 5.14, 5.15, 5.16 5.13 A sphere of raius a carries a uniform surface-charge istribution σ. The sphere is rotate about a
More information18. Atmospheric scattering details
8. Atmospheric scattering detais See Chandrasekhar for copious detais and aso Goody & Yung Chapters 7 (Mie scattering) and 8. Legendre poynomias are often convenient in scattering probems to expand the
More information1D Heat Propagation Problems
Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2
More informationV.B The Cluster Expansion
V.B The Custer Expansion For short range interactions, speciay with a hard core, it is much better to repace the expansion parameter V( q ) by f( q ) = exp ( βv( q )), which is obtained by summing over
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More informationis the scattering phase shift associated with the th continuum molecular orbital, and M
Moecuar-orbita decomposition of the ionization continuum for a diatomic moecue by ange- and energy-resoved photoeectron spectroscopy. I. Formaism Hongkun Park and Richard. Zare Department of Chemistry,
More informationStochastic Variational Inference with Gradient Linearization
Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,
More informationChapter 5. Wave equation. 5.1 Physical derivation
Chapter 5 Wave equation In this chapter, we discuss the wave equation u tt a 2 u = f, (5.1) where a > is a constant. We wi discover that soutions of the wave equation behave in a different way comparing
More informationMotion in Spherically Symmetric Potentials
Chapter 7 Motion in Sphericay Symmetric Potentias We escribe in this section the stationary boun states of quantum mechanica partices in sphericay symmetric potentias V (r), i.e., in potentias which are
More informationFoundations of Complex Mechanics
Foundations of Compex Mechanics Chapter Outine ------------------------------------------------------------------------------------------------------------------------ 3. Quantum Hamiton Mechanics 56 3.
More informationAssignment 7 Due Tuessday, March 29, 2016
Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges
More informationAPPENDIX C FLEXING OF LENGTH BARS
Fexing of ength bars 83 APPENDIX C FLEXING OF LENGTH BARS C.1 FLEXING OF A LENGTH BAR DUE TO ITS OWN WEIGHT Any object ying in a horizonta pane wi sag under its own weight uness it is infinitey stiff or
More informationTHE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS
ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP
More informationJost Function for Singular Potentials
Jost Function for Singuar Potentias S. A. Sofianos, S. A. Rakityansky, and S. E. Massen Physics Department, University of South Africa, P.O.Box 392, Pretoria 0003, South Africa (January 2, 999) An exact
More informationSTABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION
Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,
More information