Electronic Hamiltonians - part I
|
|
- Janice Thomas
- 5 years ago
- Views:
Transcription
1 Electronic Hmiltonins - prt I In this section we discuss how to del with electronic Hmiltonins, i.e. the first step in our sequence of pproimtions. Becuse we consider the ions to be frozen in their equilibrium positions, the electrons re moving in their resulting periodic potentil, nd lso eperience electron-electron interction. We will first ignore the electron-electron interctions, i.e. set them to zero. Although this my seem very strnge things to do, it is ctully firly good pproimtion in some cses. We will revisit this issue towrd the end of the section, but roughly speing, wht hppens is tht the repulsion between two electrons is nerly compensted (screened) by the ttrction with the positive chrges present in the sme region. This compenstion wors very well for electrons tht re fr prt, less so if they re sptilly close. evertheless, we will begin by ssuming tht we cn completely ignore electron-electron interctions. This llows us to study ectly some models of non-intercting electrons moving in periodic potentils, to understnd the typicl spects of electronic bnd-structures. Then, we will consider the role of (we) electron-electron repulsions using the Hrtree-Foc pproimtion, nd we will revisit screening using the Rndom Phse Approimtion. on-intercting electrons in periodic potentil For simplicity, let s ssume tht we put together mny Li toms ( s 2 2s electronic structure) in crystl, nd see how we would model this if we cn ignore electron-electron interctions. From wht we hve discussed, we epect the 2s electrons to be the vlence electrons, so we re considering vlence electron per Li + ion. For simplicity, let me consider D model of chin of such ions, plced equidistntly from ech other (of course, in relity Li would crystllize s 3D crystl, but you should be ble to del with 3D lttices if you understnd how the D cse wors). We cn lbel these with n inde n =,...,, where we will let, in the end. Let be the lttice constnt, so the position of the n th ion is R n = n. To describe the vlence electrons, we introduce cretion nd nnihiltion opertors c nσ,c nσ for n electron t site n with spin σ. To be more precise, c nσ 0 is the one-electron stte for which,σ c nσ 0 = δ σ,σ φ 2s ( R ) = δ σ,σ n, i.e. in fct the electron is creted into the 2s orbitl of the n th ion. But becuse this is the only ind of orbitl tht we include in our model, we cn sy tht the electron is t site n without creting confusion. This 2s orbitl hs some chrcteristic lenghtscle B (its Bohr rdius) over which it spreds. Our pproimtion to eep only these orbitls is resonble pproimtion if the distnce between neighbour ions is much bigger thn the size of the 2s orbitl, so tht there is no confusion to which orbitl n electron belongs. This is setched below to the left, nd is nown s tight-binding limit. One could lso consider the nerly-free electron limit, which would be vlid if the tomic orbitls spred over mny lttice sites, so tht you couldn t be sure to which orbitl n electron which is close to ion n, let s sy, belongs in relity. This is setched to the right. We will discuss this second cse lter on. tight binding limit nerly free electron limit Figure : Contrst between tightly-bound model, with orbitls smll compred to the lttice constnt (left) nd nerly-free model, with orbitls spred over mny lttice constnts.. Tight-binding models For our tight-binding model, we ssume n m = dφ 2s ( R n )φ 2s ( R m ) = δ nm. This orthogonlity is necessry in order to be ble to ssocite cretion nd nnihiltion opertors tht stisfy the proper lgebr to these sttes
2 (remember tht we cn do so only for bsis of orthonorml one-prticle sttes. If the orbitls were not orthogonl, we would first hve to define new bsis of orthogonl sttes). Deriving Hmiltonin for this pproimtion proceeds just lie we did in the Hubbrd cse, nd we will get the sme result. If we me the dditionl pproimtion tht we cn ignore even the on-site repulsion U 0, we obtin the D tight-binding Hmiltonin: Ĥ = ǫ n,σ c nσc nσ t n,σ ( c nσ c n+,σ +h.c. ) () The first term is just the energy to hve the electron t ny given site, which up to some smll corrections cn be thought of s the tomic energy of this orbitl. If ech tom eeps its electron 00% of the time, for instnce when nd we hve collection of isolted toms, the energy would be ǫ. The second term describes hopping: if the electron is t site n, ttrction from the neighbouring ions cn me it move either to n+ or to n. This is nown s nerest-neighbour hopping. As we discussed, the electron will feel ttrction from ions n+2 nd n 2 nd ions further wy, s well, nd could hop directly two or more lttice sites, to one of those, so in principle longer-rnge hopping should be present in the Hmiltonin. However, remember tht the hopping mtri decreses eponentilly with the distnce between sites, so nerest-neighbour hopping is by fr the lrgest term. In mny cses it is good pproimtion to ignore longer-rnge hopping, nd we will do so. It cn be esily included nd you should now how to do tht nd be ble to figure out precisely wht its effects re. One issue we hve to be creful with is wht hppens t the ends of the chin. I will use the so-clled periodic boundry conditions (PBC), which stte tht if the electron is t site, it cn hop either to 2 or to, nd if the electron is t site, it cn hop either to - or to. In other words, it s s if the chin is ctully closed into circle, nd nd re nerest neighbours s well. The reson we prefer to use PBC in ll our problems is becuse it mes the clcultions simpler, s we ll see soon. The logic is tht for mcroscopic crystl, with 0 23, it should not mtter how we end the crystl, becuse the vst mjority of the contributions comes from the bul sites nd the smll contributions from the surfce must be negligible by comprison. So we cn choose whichever sitution mes clcultions esiest, nd PBC do. (As good physicist you should be septic of everything I sy nd you should chec this sttement, by repeting the clcultion we do net for chin with open ends, i.e. where n electron t site cn only hop to 2, nd one t cn only hop to -. Let me now if you need some help with tht). The nice thing bout Hmiltonins lie Eq. (), which re nown s qudrtic Hmiltonins becuse ll terms re products of two opertors, is tht they cn be digonlized ectly (in other words, we cn find ectly ll their eigensttes nd eigenenergies). Here we don t even hve to wor hrd to figure out the solution, becuse it is directly pointed t by the symmetry of the problem: if I plce one electron in the system, it is eqully liely to be found t ll sites becuse they re identicl, so we should wor in momentum spce nd therefore we define opertors ssocited with sttes with well defined momentum: c σ = n e irn c nσ c σ = n e irn c nσ which cn be inverted to (chec!): c nσ = e irn c σ c nσ = n e irn c σ The normliztion is defined such tht the new opertors lso stisfy the proper lgebr {c,σ,c,σ } = δ, δ σσ etc, in other words their corresponding wvefunctions re properly orthogonl. If we plug these into Eq. () nd do bit of clening up, we find tht Ĥ =,σ E()c,σ c,σ where E() = ǫ 2tcos(). In other words, to the one-electron eigenstte c σ 0 = n= ein c nσ 0 corresponds the eigenenergy E() = ǫ 2t cos(). Becuse the electrons re non-intercting, mny-electron solutions re trivil s well: the eigenstte c σ...c pσ p 0 hs the eigenenergy E( )+...+E( p ), i.e. the sum of individul contributions. 2
3 But we hve one more question to nswer: cn the prmeter be nything, or do we hve some restrictions? We sed for PBC, which cn be written more generlly s c n+,σ = c n,σ,c n+,σ = c n,σ, for ny n (we re removing/dding the electron from/to the sme site). If we put this into their dependence of c,σ, we find tht we must hve e i = = 2m,m = 0,±,±2,... So it seems tht we hve n infinite (but countble) number of eigensttes, one for ech integer. More creful thining shows tht this is not true, in fct we hve distinct vlues only for consecutive vlues of m, nd then we repet the sme eigenfunctions nd eigenvlues. Becuse we don t find something new, we cn t pretend those to be distinct sttes/opertors, so we throw them out. To see this, consider m = 0 nd m =. The first corresponds to = 0, the second to = 2/. It is esy to see tht c =0,σ = n c n,σ = c =2/,σ = c =4/,σ =..., nd in generl sttes such tht = +mg where m is n integer nd G = 2 re in fct identicl (sme cretion nd sme nnihiltion opertors). Sowectullyonlyget distinctllowed vlues. Wecouldchoosethesesbeingm = 0,,2,..., ornyother choice of consecutive integers. Becuse we lie symmetry, the trditionl choice is m = 2 +,...,,0,,... 2, i.e. vlues symmetriclly plced bout the origin. In other words, we only consider vlues such tht <. This intervl is the Brillouin zone. The llowed vlues for nd the corresponding eigenenergies in chin with = 2 sites re setched below. E() ε llowed vlue, if =2 corresponding E() Figure 2: Allowed moment nd corresponding eigenenergies E() (white nd blue circles, respectively) inside the first Brillouin zone <. For chin with sites there re precisely llowed vlues, equidistntly plced. Snity chec: is it resonble tht we hve 2 eigensttes ( possible vlues for, times 2 for spin-up nd spindown)? Mybe we didn t find ll of them we lwys hve to worry bout this when we guess the solution, lie we did here. The nswer is: ll is good. Our Hilbert spce is 2-dimensionl, becuse it ws spnned by the 2s sttes n =,..., nd σ = ±. The eigensttes form bsis, so there must be precisely 2 of them, s well. We found 2, so we hve ll of them (by the wy, I forgot to mention tht these eigensttes must be orthonorml. You cn chec by direct clcultion tht this is indeed the cse). So wht is the solution for our chin of Li, i.e. when we hve electrons in the system so tht precisely hlf of the sttes re filled? We hve to plce the electrons in these one-prticle sttes, obeying Puli s principle. In the ground-stte they occupy the lowest vilble levels, to minimize the totl energy. Tht mens tht electrons occupy the hlf of the sttes with energies E() < ǫ (these re the sttes with moment 2 < 2 ), nd we leve empty the other hlf of the sttes, with energy E() > ǫ. So clerly the totl energy is lower thn the ǫ we would hve if ech Li ept its electron. So the ground-stte eigenstte is: GS = c,σ 0 2 2,σ 3
4 nd the ground-stte energy is: E GS = 2 2,σ E() = ǫ 4t 2 <ǫ< 2 cos() In the limit, the sum cn be turned into n integrl. This should not be surprising, becuse the llowed vlues of come closer nd closer together. As you should be ble to show, in this limit: m f() = 2 = min m min df() where min nd m re the limiting vlues for the llowed, in our cse ±/2. If you do the integrl, you ll find tht: E GS = ǫ 4 t < ǫ So there is gin of 4t/ per tom in cohesion energy, becuse the vlence electrons become deloclized over the whole system (ech one cn be found nywhere in the system with equl probbility). This gin increses if we bring the toms closer nd t becomes lrger, but then we strt to py the price due to repulsion between the clouds of core electrons (this is included in ǫ) so there is some equilibrium distnce between ions where the totl energy is minimized. For other types of lttices (nd in higher dimension), the epression of E() will be different, nd s result the prefctor multiplying t chnges, so different lttices led to different cohesion energies. The one with lowest energy will be the preferred structure tht will be observed in relity..2 Free electrons (jellium model without interctions) Remember tht the jellium model ssumes tht the ions re spred out in uniform density of positive chrge. If we ignore the electron-electron interctions, the resulting Hmiltonin (see previous chpter, but now for D system) becomes: Ĥ =,σ ǫ,σ,σ where ǫ = h2 2 2m is the free-electron inetic energy, nd the opertors re ssocited with plne-wves with momentum :,σ,σ 0 = δ σ,σ ei / L, where L is the length of the D crystl. ote tht these opertors re very different from the c,σ we found in the previous section: those creted/nnihilted electrons in very different sttes, which were liner combintions of 2s orbitls, not these plne-wves! This is why I ll cll them insted of c. This difference is lso cler if we loo t the llowed vlues in this cse. If we gin impose PBC, this mens tht nothing should chnge if we trnslte + L. From the definition bove, this obviously hppens only if e il = = 2 L m,m = 0,±,±2,... This loos similr to wht we hd before (there L = ) ecept now there is no restriction on the vlues of m! For the tight-binding chin, the restriction ppered becuse the vlues of the plne-wves only mttered t R n = n. However, in this jellium model they mtter t ech [0,L], so ech llowed vlue gives different, cceptble solution. Thus ny integer m nd so ny = 2 L m is llowed, in other words there is no restriction to Brillouin zone lie before. If you thin bout it, this jellium model ctully mimics well the nerly-free electrons limit shown on the right-hnd side of Fig.. The limit B is obtined when we let 0 indeed, in this cse insted of distinct ions, we see continuous distribution of positive chrge. Hving 0 lso mens tht the Brillouin zone etends to cover the whole is since ± ±, which is wht we just found; nd finlly, in this limit we find E() = ǫ 2tcos() ǫ 2t + (t 2 ) 2. If we define n effective mss m such tht h2 2m = t2, then up to n unimportnt constnt, E() ǫ s well (it is resonble tht t 2 might sty finite s 0, becuse remember tht hopping integrls increse with decresing distnce). So formlly the results mp onto one nother, lthough the wvefunctions re somewht more problemtic. If we hve electrons in this system, in the GS they will gin occupy ll sttes up to F, so: GS = < F,σ,σ 0 4
5 nd where the Fermi momentum is such tht: E GS = ǫ = 2 ǫ < F,σ < F = = 2 < F,σ < F Agin, in the limit L, the llowed vlues of become dense nd we cn replce the sums with integrls: 2 f() = L 2 2 df(). We then find: = 2 F L/ F = n/2, where n = /L is the electron density, nd we cn then go nd clculte the GS energy EGS = h2 2 n 2 24m nd ny other quntities tht we re interested in. ote: I epect you to be ble to do similr clcultions in higher dimensions. I epect you sw them before in undergrdute-level solid stte nd lso sttisticl mechnics courses. If they give you trouble, come nd tl to me..3 Free electrons in periodic potentil: bnd structures At this point you re probbly feeling bit confused becuse these two cses loo very different from one nother, so it s not ll tht cler wht they hve in common nd how one my go from one to the other by vrying the rtio B /. I thin the best wy to understnd this is to ctully solve problem we ll do this net for the simplest non-trivil emple tht I now (still in D). Let s gin ssume tht we hve simple crystl with single tom in the bsis, lthough gin, things generlize strightforwrdly to more complicted cses. Suppose tht we hve ions plced t their lttice loctions X n = n, where is the lttice constnt nd n =,..,, nd vlence electrons whose behvior we re trying to understnd. We will gin ssume PBC nd t the end we will let, so tht the boundry condition becomes irrelevnt. Then, the totl potentil felt by ny of the vlence electrons is: V tot () = V ( n) n= where V () is the potentil creted by n ion locted t the origin, nd the sum is simply the contribution of ll the ions. Becuse there re no electron-electron interctions we cn solve the problem for single electrons, nd then me Slter determinnt from the one-electron eigensttes. So we need to solve Schrödinger s eqution: [ h2 d 2 ] 2md 2 +V tot() φ() = Eφ() with the periodic boundry condition φ(+) = φ(), for ny. At this point I epect you remember tht such periodic potentils lwys led to so-clled Bloch eigensttes. I could jump directly to using them, but becuse this is ey point, let me briefly review one wy to derive them (if this is not enough to jog your memory, plese do some reding!). Let s try to guess wht the solution of this Schrödinger eqution must loo lie: (i) the solution must hve plne-wve lie chrcter, becuse the symmetry of the problem tells us tht the probbility to be in ny unit cell must be the sme, irrespective of which unit cell we re tling bout. However, (ii) we epect tht the wvefunction vries from plce to plce inside the unit cell; for emple, the probbility to find the electron should be higher ner n ion thn in between two ions. Putting these together, let s try the guess: φ () = e i u () where u () = u ( + ) is periodic function. This hs the right chrcter, becuse it leds to probbility φ () 2 = u () 2 which is periodic function, just s we wnt it to be. Another wy to see this is tht, with our guess: φ (+) = e i(+) u (+) = e i e i u () = e i φ (), 5
6 i.e. moving one unit cell over leves the mgnitude of the wvefunction unchnged, but it dds phse, just s it should (remember, we wnt mny eigenfunctions, not just the one with φ(+) = φ(), nd we need the plne-wve for tht). By the wy, in 3D, the guess would be φ ( r) = e i r u ( r), where now u ( r) is periodic in the 3D unit cell. It turns out tht this guess wors perfectly. In fct, one cn prove tht ny single-prticle Schrödinger s eqution with periodic potentil hs eigenfunctions of this type, which re nown s Bloch wvefunctions. Wht is left to do is to figure out wht is u () nd its corresponding eigenenergy E, so tht the eqution is stisfied. This is mtter of mth. Before discussing how to do tht, let s loo more crefully t the momentum. For our D crystl, the wvefunction must lso stisfy the periodic boundry condition φ(+) = φ(). For Bloch wvefunction, though, φ(+) = e i φ(), so we must hve = 2 m, where m is n integer. One cn gin show tht only of these solutions re distinct (just lie the discussion we hd for the tight-binding model), so we must restrict (, ], i.e. to the first Brillouin zone. This is lwys the cse when we hve model with lttice. Before looing t how this wors in simple but non-trivil emple, let s first me sure tht we regin the correct solution if U () 0. This, of course, leds us bc to free electrons, but zero is lso periodic potentil so we should be ble to put the free electron solution in the Bloch form, nd with the restriction on the vlues of if we dopt the periodic boundry condition. The wy this wors is illustrted in the figure below. On the left, I show the dispersion in free spce: the usul prbol with ny momentum llowed. If I choose to view the free spce s periodic lttice (without ny potentil), then the correct description is shown on the right. I plced prbol t ech reciprocl lttice point, so now the eigensttes re periodic but I only count the sttes inside the first Brillouin zone, since the other unit cells re just repeting these nd don t give nything new. So now we ve restricted the momentum, but for ech vlue of we hve mny solutions E n (0) (), n = 0,,2,... Overll we hve the sme spectrum, it s just different wy to inde the sttes. By the wy, this re-rrngement of the sttes to fll ll inside the Brillouin zone is now s folding. E() E() 4 2 Brillouin zone 2 4 Figure 3: Free spce eigensttes. Left: the usul prbol for free spce, where ny is llowed nd we hve single eigenstte for given. Right: the corresponding bnd-structure if we te the vcuum to be periodic crystl with lttice constnt (nd no potentil). ow the moment re restricted to the BZ, but for ech we hve n infinite number of bnds E n (0) (). The thic coloured lines show the n = 0,,2 bnds. For ny given in the BZ, we cn see to wht K of the originl prbol tht prticulr energy corresponds. I showed 2 such correspondences (horizontl lines with rrows). To be bit more precise, let s consider the first 3 bnds, n = 0,,2, for given (see right figure). I illustrted the corresponding sttes in both representtions in the usul free spce lnguge, they correspond to different moment, K 0 =,K = 2,K 2 = + 2. The energies, then, re E(0) 0 () = h2 2 2m = E 0(K 0 ); E (0) () = h2 ( 2 )2 2m = E 0 (K ); E (0) 2 () = h2 (+ 2 )2 2m = E 0 (K 2 ), nd so on nd so forth, I hope you gree tht we cn mp ll the points in one spectrum to points in the other. How bout the wvefunctions, cn we put them in the Bloch form? For n = 0 this is trivil, since here φ (0) K 0 () = e ik0 = e i = φ 0, (), so it hs the epected Bloch form φ 0, () = e i u 0, () with u 0, () =, which is certinly periodic. How bout for K? In this cse, the eigenstte is: φ (0) K () = e ik = e i( 2 ) = e i u, () = φ, () 6
7 where for the n = bnd, it follows tht we must hve u, () = e i2. Is this periodic function? The nswer is yes: u, ( + ) = e i2 (+) = e i2 e i2 = u, (). And so on nd so forth, you cn chec tht for ech bnd n, the phse difference between the corresponding K n nd is multiple of 2/ nd gives different periodic prt u n, () for the Bloch wvefunction φ n, () = e i u n, () of tht bnd. By the wy, it is good tht different bnds hve different u n, (), becuse wvefunctions with different n nd/or different hve to be orthogonl. So this scheme wors for free spce, if we choose to view it s periodic potentil. As the rule sys, we find tht is indeed restricted to the Brillouin zone. The price is tht we hve mny bnds (different eigensttes for the sme ), ech of which hs Bloch wvefunction with different periodic prt. Turns out tht this is wht hppens in the generl cse, too we will hve (infinitely) mny bnds in the Brillouin zone, however their energies E n () will, of course, depend on the prticulr potentil we choose, s will the periodic prts of their Bloch wvefunctions, φ n, () = e i u n, (). To the find the specific E n () nd u n, (), we need to solve Schrödinger s eqution. Let me show you how it wors in simple cse when we hve the simplest, δ-function ttrction between ions nd electrons, i.e. V () = U 0 δ() (By the wy, why does the constnt in front hve units of energy distnce if V () is n energy? Figure tht out. Once you do tht, you ll see tht we need to use the lttice constnt, which is our nturl lengthscle in this problem, s the length. Tht leves U 0 s n energy which we cn djust t will: the lrger it is, the stronger the ttrction). First, let me quicly review the solution for single ion. I epect you hve solved this problem in n introductory QM course, nd even if not tht you cn figure it out but s for pointers if you need help. Of course, there re eigensttes for ll positive energies E > 0 (when the electron is not bound to the ion), but it turns out tht there is one bound stte with wvefunction (up to normliztion constnt): nd energy φ B () e κ = e B E B = h2 κ 2 2m = U h 2 m 2 where the Bohr rdius is B = κ = h2 mu 0 defines the spred of the wvefunction. This tells us tht the lrger U 0 is (the stronger the ttrctive potentil), the more negtive the energy of the bound stte, nd the smller B is. This mens tht by decresing U 0, we cn go from the limit where B, which is the tight-binding limit where our tight-binding model should wor well, to the limit with B, which is the nerly free electron limit. These re setched below. B B n n n+ n n n+ tight binding limit >> B nerly free electron limit << B Figure 4: Left: tight-binding limit, where B. Right: nerly free electron limit, where B. For single ion, then, the spectrum for negtive energies hs only this one level t energy E B. Let s see wht hppens for lttice of such ions, t negtive energies E = h2 λ 2 2m ; for wht vlues of λ do we now find eigensttes? We need to solve the SE for the potentil U tot () = U 0 n=0 δ( n). The Bloch theorem gurntees tht the wvefunctions must be such tht φ ( + ) = e i φ () this is very nice, becuse it mens tht we only need solve SE inside one unit cell, nd mtch the boundries ccordingly. So consider the unit cell 0 < <. The potentil t ll these vlues is zero, so the wvefunction in this intervl must be φ () = αe λ +βe λ. According to Bloch s theorem, the wvefunction in the net unit cell, < + < 2, 7
8 must be φ (+) = e i φ () = e i( αe λ +βe λ). All we need to do now, is to sew together the wvefunctions t =, where we hve delt function in the potentil. Continuity t = implies: αe λ +βe λ = e i (α+β) while the jump in the derivtive becuse of the δ-function results in: h2 [ e i (αλ βλ) ( αλe λ βλe λ)] = U 0 e i (α+β) 2m So we hve two homogeneous equtions for α,β, nd we now tht non-trivil solution is possible only if the determinnt is zero. After some boring lgebr, this condition becomes: cosh(λ) sinh(λ) λ B = cos() where cosh() = 2 (e +e ), sinh() = 2 (e e ). So wht hppens here is tht for ech vlue of (inside the Brillouin zone), this eqution will give us one or more solutions λ n (),n =,2,.. (however mny solutions we get) nd therefore we find the llowed energies E n () = h2 λ n() 2 2m. This eqution is too complicted to dmit simple nlyticl solutions, but we cn still figure out lot bout its solutions in vrious symptotic cses. Let me define the function: f() = cosh() sinh B Then λ = / is the solution when f() = cos(). Since cos(), we need to figure out wht s the shpe of f(), nd in prticulr when it hs vlues between,. At the origin, f(0) = B. This is definitely less tht. If > 2 B, then f(0) <. As increses, you should be ble to convince yourself tht this function increses monotoniclly. In prticulr, t lrge we cn discrd the e eponentils s being smll, nd we find f() e 2 ( ), which diverges s. B f() =0 f() =0 B = 2 λ( 2 ) λ(0) λ( ) λ( ) 2 λ(0) = 2 B = Figure 5: f() setched when (left) < 2 B, nd (right) > 2 B. The solutions for f() = cos() with (, ] lie in between between -,. This function is setched bove, for both < 2 B nd > 2 B cses. Let me strt with the second cse which contins the tight-binding limit B. Here, there is single solution for f() = cos() [,], so we epect single solution for λ = /, nd therefore single bnd t negtive energies E = h2 λ 2 2m. Let s find it in the limit B, i.e. in the tight-binding limit. ow f(0) strts t very negtive vlue, so by the time it crosses the [,] intervl the vlues of re very lrge nd we cn use the simpler symptotic epression, i.e. here: f() e 2 ( ) B = cos() = B [ 2e cos()] This is still trnscendentl eqution so we don t hve n ect solution, but we cn get n ccurte one by itertions. If the eponentil wsn t there, the solution would be = B in this limit (which is lrge vlue, consistent with 8
9 our symptotic pproimtion). As result, e is smll nd leds to tiny correction, so the ctul solution won t be too fr from B. So I cn replce the in the eponentil with / B, to find: = λ() [ ] [ ] +2e B 2e B cos() B cos() B (this is just Tylor epnsion, since the eponentil continues to be smll). Remembering tht / B = κ for the bound stte, we find: ] λ() = κ [+2e B cos() leding to the energies: E B () = h2 λ() 2 2m E B ) (+4e B cos = E B 2tcos() if we define the hopping t = 2 E B e B. Before continuing, note tht I dropped the term tht goes lie e 2 B from the squre this is becuse I now tht it is much smller thn the ones I ept, nd moreover I lredy neglected such terms when I did the Tylor series bove, so it s simply not meningful to eep it here. So for B we regin the tight-binding eigenenergy we epected to be vlid here. Moreover, we lso hve the correct epression for the hopping t (if you remember how this is defined s the overlp between neighboring orbitls, you cn chec this sttement). The left pnel below shows the llowed energy bnd, i.e. the energy intervl where we find solutions in the limit when, so tht ll inside the BZ re llowed. The right picture shows the dispersion for one prticulr vlue of / B, in the tight-binding limit. The lrger / B is, the nrrower the bnd becomes, becuse the hopping between nerest sites becomes eponentilly smll. In the limit / B we bsiclly hve isolted toms, ech with their llowed level t E B. E 2 llowed energies B E 0() E 2tcos() B E B Figure 6: Left: llowed bnd of energies with E < 0, s function of / B. Right: For prticulr vlue of / B, indicted by the dshed line, I showed the dispersion E B () vs, in the Brillouin zone. How bout the wvefunctions, do they come out correctly? In the tight-binding limit we obtined the solutions c,σ 0 = n ein c n,σ 0, whose wvefunctions re (ignoring the spin) φ () = c 0 = n=0 ein φ B ( n). ote tht these re indeed Bloch sttes! (s they should be, the tight-binding model corresponds to periodic potentil.) Remember tht the bound level eigenfunctions re simple eponentils. For 0 < <, the ones tht contribute most re the e /B prt from the ion t n = 0, nd the e B prt from the ion t n = ll other ions contributions re very much smller thn these terms, since they re further wy from this intervl. Keeping only these two lrgest contribution, the tight-binding pproimtion sys tht for 0 < < we should find: φ () e B +e i e B 9
10 Let sseewhttheectsolutionpredictsinthislimit. Thisrequiresustofindα,β sothtwehvethewvefunction. We cn use either eqution, let s use the simpler one: αe λ +βe λ = e i (α+β) α β = ei e λ e λ e i But in this limit λ B e λ = e B, thus we cn simplify the frction: α β ei e λ α = βei e B. Putting this into our solution on the 0 < < intervl, φ () = αe λ + βe λ αe B + βe B (since λ B ) gives us, up to normliztion constnt, the epected solution. So our simple tight-binding model gives the correct eigenenergies nd eigenfunctions in the correct limit. Wht hppens when we decrese / B, i.e. we bring the ions closer nd closer together? For one, the hopping increses so the bndwidth 4t becomes lrger. From the discussion of the ect f(), you cn see tht something hppens t / B = 2 only bove it re we gurnteed to find solutions for ll in the BZ. Wht hppens for smller / B, s the figure bove hopefully suggests, is tht we need to worry bout positive energies (so fr we only discussed E < 0 eigensttes). So let s loo t the spectrum there, too. The E = h2 K 2 2m > 0 cse is studied just lie before. Everything mirrors wht we did for E = h2 λ 2 2m < 0 if we replce λ ik. In prticulr, we now find tht for ech, the llowed vlues of K (which determine the energy) re given by the eqution: cos(k) sin(k) = cos() K B Let s define gin function g() = cos() sin B, nd solve g() = cos() K = /. Wht does g() loo lie? Its two components re setched below; the cos function oscilltes, while the sin oscilltes but lso decys slowly with incresing. cos sin 2 Figure 7: Left: cos() vs. Right: sin()/ vs. Let s see first wht hppens for B. In this cse, the sin term is lrge ner the origin nd domintes there, so g() loos lie in the net plot. I lso drew the lines for g() = (this gives solutions for = 0) nd g() =, which corresponds to solutions for = ±. As we increse we sweep between the two. Unlie for E < 0, where we found single solution for ech, here we hve n infinity of them, so they predict n infinity of bnds t positive energies. The full spectrum loos lie setched on the right, with the tight-binding bnd t negtive energy, nd ll these other llowed bnds seprted by gps t positive energies. ote how the curvture of the bnds lterntes. O, so now finlly let s see wht hppens s / B becomes smll. In this limit, the function g() loos lie below. The contribution of the sin term is very smll, so g() brely crosses bove or below -. We still get llowed bnds seprted by forbidden energy rnges, but they re wider. The bottom of the lowest bnd nicely mtches up with wht we found there for the negtive spectrum, so the whole spectrum here loos lie shown in the right pnel: This is strting to loo lie the free electron spectrum, which it should evolve into when / B 0 (for B when U 0 0, i.e. when there is no potentil). Indeed, s we me / B smller, the bottom of the lowest bnd comes closer to zero while the gps bove strt to close, nd we regin the folded bnd structure from Fig.. If you wish to chec this sttement, you cn clculte, for emple, the energies of the first two bnds t = /. I did this, nd found E 0 ( = ) = ( [ ] h2 2 2m ) 4 2 B, while for the second bnd, I find E ( = ) = ( [ h2 2 2m ) B ]. So the 0
11 E() E g() =0 = /2 = / bnd of llowed energies forbidden gp B E B Figure 8: Left: Solution of g() = cos() in the limit / B. Right: Full spectrum in the limit B. ( first bnds ends just below the h2 ) 2 2m epected in vcuum for =, while the second bnd strts just bit bove it. The gp between them is proportionl to / B, nd indeed closes s U 0 0. This clcultion tells us tht in the presence of periodic potentil, we genericlly epect to see spectrum consisting of llowed energy bnds which re seprted by gps. Of course, rel toms hve mny bound sttes, not just one lie our pretend toms. However, this is not going to chnge things significntly. Ech level hs chrcteristic distnce over which it decys eponentilly, nd deeper levels decy fster becuse they re locted closer to the nucleus. For given lttice constnt, the low-lying levels decy much fster thn so they re well-described by tight-binding models. For higher-energy orbitls, which re more etended, the tight-binding pproimtion is not going to be ccurte nymore, nd we ll need to do better. In ny event, if we consider wht hppens s we decresed the lttice constnt from infinity (corresponding to isolted toms, ech with its own tomic spectrum), we should get something lie setched below: The precise loctions of the llowed energy bnds will, of course, depend on wht type of crystl we hve, nd wht is the equilibrium vlue for, etc etc. But there will certinly be energy bnds chrcterized by momentum restricted to the first Brillouin zone, seprted by forbidden gps. Higher bnds hve lrger hopping constnts becuse their orbitls re more spred-out, so their width increses fster; they re lso closer together, so we epect tht they strt overlpping t some point or other (when this hppens, the simple tight-binding pproimtion definitely fils nd we need to wor hrder to find n ccurte solution). So this is the generic picture for the single-electron spectrum of crystl. ow we hve to fill the lowest levels with ll the electrons contributed by ll the toms. There re two possible outcomes: (i) the lst bnd tht contins electrons is only prtilly full. In this cse, the mteril is metl if we pply smll electric field, there re free levels just bove the occupied ones, so we cn esily ecite electrons nd me them move preferentilly in the direction of the field, to get n electric current (we ll discuss this more crefully very soon). (ii) the lst bnd tht contins electrons is completely full, in the ground-stte. In this cse, the mteril will not E() E bnd of llowed energies g() =0 B = /2 = / forbidden gp Figure 9: Left: Solution of g() = cos() in the limit / B. Right: Full spectrum in the limit B.
12 E gps E E 2p 2s bnds of llowed energies E s tomic lie core sttes Figure 0: Generl evolution of the energy bnds in solid, s the lttice constnt is decresed. be metl; to ecite electrons to empty levels, we need to give them n energy t lest equl to the gp so tht we cn move them into the net bnd. This is why pplying smll electric field will not strt current flowing, so this is not conductor. If the gp is lrger thn 2eV, we cll this n insultor; if it is less thn 2eV, we cll such mteril semiconductor. Bsed on the solution for the tight-binding model, you might epect tht ny time we hve one (or n odd number) of vlence electrons per ion we should get metl. This is not true! To obtin the ground-stte with the lower-hlf of the tight-binding bnd full, we mde the crucil pproimtion tht we cn ignore electron-electron interctions. If those interctions cnnot be ignored, then the mny-electron wvefunction is not Slter determinnt of one-electron sttes nd these rguments fil, nd we hve to clculte the mny-body wvefunction somehow else. It turns out tht for Hubbrd model in 3D, if U is sufficiently lrge, gp opens up nd splits the occupied sttes from the empty sttes, so the system is n insultor (in D, this hppens for ny U > 0). We distinguish such insultors from the ones discussed bove, which occur even in the bsence of interctions. Those re clled bnd insultors (nd re bit boring lthough etremely importnt for current technologies), wheres the ones tht rise becuse of electron-electron interctions re clled Mott insultors. We ll discuss bit more bout these lter on. Before finishing this section, let me mention tht we cn mesure directly the energy bnds using photo-emission spectroscopy. This is when bem of high-energy photons is shined onto smple, nd some of them re bsorbed by electrons which therefore receive enough energy to escpe the crystl. Detectors mesure the energy nd momentum of the escping electrons, nd from conservtion of energy nd momentum, one cn etrct informtion bout wht ws the energy nd momentum, i.e. E n (), of the electrons while in the crystl. This technique is nown s Angle-Resolved Photoemission Spectroscopy (ARPES), nd the Dmscelli lb here t UBC does such mesurements. You cn loo t their webpge for more informtion nd some nice pictures of eperimentlly mesured electronic bnds. In terms of clculting bnd structures going beyond simple models lie tight-binding, one uses density functionl theory, lso nown s b-initio methods. These re bsed on the Hohenberg-Kohn theorem, which is n mzing result plese te bit of time to red bout this in stndrd tetboo (for emple, Tylor nd Heinonen pp ). These dys there re free (nd commercil) softwre pcges to del with this pproch, nd some of you my end up using some of these. Even if you don t, you should hve bsic understnding of the underlying ides, which re firly simple to grsp. Let me just sy tht lot of cre is needed with the interprettion of these results, though, becuse this theory only gurntees to predict the (totl) ground-stte energy nd density of electrons in the ground-stte, nd this only if certin functionl is nown. We do not now tht functionl but we cn compute n pproimtion for the jellium model, nd tht is used for ny model (the so-clled LDA pproimtion). This wors pretty well for systems with we electron-electron interctions, but not for systems with strong interctions. Vrious improvements nd wys to fi the problem re proposed for the ltter problems, but this is still very much wor in progress. 2
Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon
Phys463.nb 49 7 Energy Bnds Ref: textbook, Chpter 7 Q: Why re there insultors nd conductors? Q: Wht will hppen when n electron moves in crystl? In the previous chpter, we discussed free electron gses,
More informationName Solutions to Test 3 November 8, 2017
Nme Solutions to Test 3 November 8, 07 This test consists of three prts. Plese note tht in prts II nd III, you cn skip one question of those offered. Some possibly useful formuls cn be found below. Brrier
More informationand that at t = 0 the object is at position 5. Find the position of the object at t = 2.
7.2 The Fundmentl Theorem of Clculus 49 re mny, mny problems tht pper much different on the surfce but tht turn out to be the sme s these problems, in the sense tht when we try to pproimte solutions we
More informationThe First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).
The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples
More informationChapter 0. What is the Lebesgue integral about?
Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous
More information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
More informationMath Lecture 23
Mth 8 - Lecture 3 Dyln Zwick Fll 3 In our lst lecture we delt with solutions to the system: x = Ax where A is n n n mtrix with n distinct eigenvlues. As promised, tody we will del with the question of
More informationProblem 3: Band Structure of YBa 2 Cu 3 O 7
HW 5 SSP 601-2017. here is very relistic clcultion which uses the concepts of lttice, reciprocl spce, Brillouin zone nd tight-binding pproximtion. Go over the solution nd fill up every step nd every detil
More informationChapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1
Tor Kjellsson Stockholm University Chpter 5 5. Strting with the following informtion: R = m r + m r m + m, r = r r we wnt to derive: µ = m m m + m r = R + µ m r, r = R µ m r 3 = µ m R + r, = µ m R r. 4
More informationChapter 3 The Schrödinger Equation and a Particle in a Box
Chpter 3 The Schrödinger Eqution nd Prticle in Bo Bckground: We re finlly ble to introduce the Schrödinger eqution nd the first quntum mechnicl model prticle in bo. This eqution is the bsis of quntum mechnics
More information5.04 Principles of Inorganic Chemistry II
MIT OpenCourseWre http://ocw.mit.edu 5.04 Principles of Inorgnic Chemistry II Fll 2008 For informtion bout citing these mterils or our Terms of Use, visit: http://ocw.mit.edu/terms. 5.04, Principles of
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk bout solving systems of liner equtions. These re problems tht give couple of equtions with couple of unknowns, like: 6 2 3 7 4
More informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More informationLecture 13 - Linking E, ϕ, and ρ
Lecture 13 - Linking E, ϕ, nd ρ A Puzzle... Inner-Surfce Chrge Density A positive point chrge q is locted off-center inside neutrl conducting sphericl shell. We know from Guss s lw tht the totl chrge on
More informationThe solutions of the single electron Hamiltonian were shown to be Bloch wave of the form: ( ) ( ) ikr
Lecture #1 Progrm 1. Bloch solutions. Reciprocl spce 3. Alternte derivtion of Bloch s theorem 4. Trnsforming the serch for egenfunctions nd eigenvlues from solving PDE to finding the e-vectors nd e-vlues
More information7.2 The Definite Integral
7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where
More information13: Diffusion in 2 Energy Groups
3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups
More informationSolutions to Problems in Merzbacher, Quantum Mechanics, Third Edition. Chapter 7
Solutions to Problems in Merzbcher, Quntum Mechnics, Third Edition Homer Reid April 5, 200 Chpter 7 Before strting on these problems I found it useful to review how the WKB pproimtion works in the first
More informationRecitation 3: More Applications of the Derivative
Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech
More informationUniversity of Alabama Department of Physics and Astronomy. PH126: Exam 1
University of Albm Deprtment of Physics nd Astronomy PH 16 LeClir Fll 011 Instructions: PH16: Exm 1 1. Answer four of the five questions below. All problems hve equl weight.. You must show your work for
More informationMassachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 6
Msschusetts Institute of Technology Quntum Mechnics I (8.) Spring 5 Solutions to Problem Set 6 By Kit Mtn. Prctice with delt functions ( points) The Dirc delt function my be defined s such tht () (b) 3
More informationNOTES ON HILBERT SPACE
NOTES ON HILBERT SPACE 1 DEFINITION: by Prof C-I Tn Deprtment of Physics Brown University A Hilbert spce is n inner product spce which, s metric spce, is complete We will not present n exhustive mthemticl
More informationImproper Integrals. Type I Improper Integrals How do we evaluate an integral such as
Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationf(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral
Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one
More informationDuality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.
Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More information1 Part II: Numerical Integration
Mth 4 Lb 1 Prt II: Numericl Integrtion This section includes severl techniques for getting pproimte numericl vlues for definite integrls without using ntiderivtives. Mthemticll, ect nswers re preferble
More informationQuantum Physics I (8.04) Spring 2016 Assignment 8
Quntum Physics I (8.04) Spring 206 Assignment 8 MIT Physics Deprtment Due Fridy, April 22, 206 April 3, 206 2:00 noon Problem Set 8 Reding: Griffiths, pges 73-76, 8-82 (on scttering sttes). Ohnin, Chpter
More informationARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac
REVIEW OF ALGEBRA Here we review the bsic rules nd procedures of lgebr tht you need to know in order to be successful in clculus. ARITHMETIC OPERATIONS The rel numbers hve the following properties: b b
More informationMath 113 Exam 2 Practice
Mth Em Prctice Februry, 8 Em will cover sections 6.5, 7.-7.5 nd 7.8. This sheet hs three sections. The first section will remind you bout techniques nd formuls tht you should know. The second gives number
More information5.5 The Substitution Rule
5.5 The Substitution Rule Given the usefulness of the Fundmentl Theorem, we wnt some helpful methods for finding ntiderivtives. At the moment, if n nti-derivtive is not esily recognizble, then we re in
More informationInfinite Geometric Series
Infinite Geometric Series Finite Geometric Series ( finite SUM) Let 0 < r < 1, nd let n be positive integer. Consider the finite sum It turns out there is simple lgebric expression tht is equivlent to
More informationProblem Set 3 Solutions
Chemistry 36 Dr Jen M Stndrd Problem Set 3 Solutions 1 Verify for the prticle in one-dimensionl box by explicit integrtion tht the wvefunction ψ ( x) π x is normlized To verify tht ψ ( x) is normlized,
More informationUNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3
UNIFORM CONVERGENCE Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 Suppose f n : Ω R or f n : Ω C is sequence of rel or complex functions, nd f n f s n in some sense. Furthermore,
More informationQuadratic Forms. Quadratic Forms
Qudrtic Forms Recll the Simon & Blume excerpt from n erlier lecture which sid tht the min tsk of clculus is to pproximte nonliner functions with liner functions. It s ctully more ccurte to sy tht we pproximte
More informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
More informationPolynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230
Polynomil Approimtions for the Nturl Logrithm nd Arctngent Functions Mth 23 You recll from first semester clculus how one cn use the derivtive to find n eqution for the tngent line to function t given
More informationPhysics 202H - Introductory Quantum Physics I Homework #08 - Solutions Fall 2004 Due 5:01 PM, Monday 2004/11/15
Physics H - Introductory Quntum Physics I Homework #8 - Solutions Fll 4 Due 5:1 PM, Mondy 4/11/15 [55 points totl] Journl questions. Briefly shre your thoughts on the following questions: Of the mteril
More informationIntegration. 148 Chapter 7 Integration
48 Chpter 7 Integrtion 7 Integrtion t ech, by supposing tht during ech tenth of second the object is going t constnt speed Since the object initilly hs speed, we gin suppose it mintins this speed, but
More informationLecture 1. Functional series. Pointwise and uniform convergence.
1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is
More informationClassical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011
Clssicl Mechnics From Moleculr to Con/nuum Physics I WS 11/12 Emilino Ippoli/ October, 2011 Wednesdy, October 12, 2011 Review Mthemtics... Physics Bsic thermodynmics Temperture, idel gs, kinetic gs theory,
More information2.4 Linear Inequalities and Interval Notation
.4 Liner Inequlities nd Intervl Nottion We wnt to solve equtions tht hve n inequlity symol insted of n equl sign. There re four inequlity symols tht we will look t: Less thn , Less thn or
More informationMath 360: A primitive integral and elementary functions
Mth 360: A primitive integrl nd elementry functions D. DeTurck University of Pennsylvni October 16, 2017 D. DeTurck Mth 360 001 2017C: Integrl/functions 1 / 32 Setup for the integrl prtitions Definition:
More informationModule 6: LINEAR TRANSFORMATIONS
Module 6: LINEAR TRANSFORMATIONS. Trnsformtions nd mtrices Trnsformtions re generliztions of functions. A vector x in some set S n is mpped into m nother vector y T( x). A trnsformtion is liner if, for
More information( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that
Arc Length of Curves in Three Dimensionl Spce If the vector function r(t) f(t) i + g(t) j + h(t) k trces out the curve C s t vries, we cn mesure distnces long C using formul nerly identicl to one tht we
More informationp-adic Egyptian Fractions
p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction
More informationP 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)
1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationAbstract inner product spaces
WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the
More information221B Lecture Notes WKB Method
Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using
More informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
More informationUniversity of Washington Department of Chemistry Chemistry 453 Winter Quarter 2010 Homework Assignment 4; Due at 5p.m. on 2/01/10
University of Wshington Deprtment of Chemistry Chemistry 45 Winter Qurter Homework Assignment 4; Due t 5p.m. on // We lerned tht the Hmiltonin for the quntized hrmonic oscilltor is ˆ d κ H. You cn obtin
More informationContinuous Quantum Systems
Chpter 8 Continuous Quntum Systems 8.1 The wvefunction So fr, we hve been tlking bout finite dimensionl Hilbert spces: if our system hs k qubits, then our Hilbert spce hs n dimensions, nd is equivlent
More informationapproaches as n becomes larger and larger. Since e > 1, the graph of the natural exponential function is as below
. Eponentil nd rithmic functions.1 Eponentil Functions A function of the form f() =, > 0, 1 is clled n eponentil function. Its domin is the set of ll rel f ( 1) numbers. For n eponentil function f we hve.
More informationLecture 8. Band theory con.nued
Lecture 8 Bnd theory con.nued Recp: Solved Schrodinger qu.on for free electrons, for electrons bound in poten.l box, nd bound by proton. Discrete energy levels rouse. The Schrodinger qu.on pplied to periodic
More informationSummary: Method of Separation of Variables
Physics 246 Electricity nd Mgnetism I, Fll 26, Lecture 22 1 Summry: Method of Seprtion of Vribles 1. Seprtion of Vribles in Crtesin Coordintes 2. Fourier Series Suggested Reding: Griffiths: Chpter 3, Section
More informationMath 131. Numerical Integration Larson Section 4.6
Mth. Numericl Integrtion Lrson Section. This section looks t couple of methods for pproimting definite integrls numericlly. The gol is to get good pproimtion of the definite integrl in problems where n
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk out solving systems of liner equtions. These re prolems tht give couple of equtions with couple of unknowns, like: 6= x + x 7=
More information221A Lecture Notes WKB Method
A Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using ψ x, t = e
More informationExam 1 Solutions (1) C, D, A, B (2) C, A, D, B (3) C, B, D, A (4) A, C, D, B (5) D, C, A, B
PHY 249, Fll 216 Exm 1 Solutions nswer 1 is correct for ll problems. 1. Two uniformly chrged spheres, nd B, re plced t lrge distnce from ech other, with their centers on the x xis. The chrge on sphere
More information2 b. , a. area is S= 2π xds. Again, understand where these formulas came from (pages ).
AP Clculus BC Review Chpter 8 Prt nd Chpter 9 Things to Know nd Be Ale to Do Know everything from the first prt of Chpter 8 Given n integrnd figure out how to ntidifferentite it using ny of the following
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term Solutions to Problem Set #1
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Deprtment 8.044 Sttisticl Physics I Spring Term 03 Problem : Doping Semiconductor Solutions to Problem Set # ) Mentlly integrte the function p(x) given in
More informationMATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1
MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 7.. Introduction to Integrtion. 7. Integrl Clculus As ws the cse with the chpter on differentil
More informationQuantum Mechanics Qualifying Exam - August 2016 Notes and Instructions
Quntum Mechnics Qulifying Exm - August 016 Notes nd Instructions There re 6 problems. Attempt them ll s prtil credit will be given. Write on only one side of the pper for your solutions. Write your lis
More informationf(a+h) f(a) x a h 0. This is the rate at which
M408S Concept Inventory smple nswers These questions re open-ended, nd re intended to cover the min topics tht we lerned in M408S. These re not crnk-out-n-nswer problems! (There re plenty of those in the
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct
More informationIntermediate Math Circles Wednesday, November 14, 2018 Finite Automata II. Nickolas Rollick a b b. a b 4
Intermedite Mth Circles Wednesdy, Novemer 14, 2018 Finite Automt II Nickols Rollick nrollick@uwterloo.c Regulr Lnguges Lst time, we were introduced to the ide of DFA (deterministic finite utomton), one
More informationQuantum Analogs Chapter 4 Student Manual
Quntum Anlogs Chpter 4 Student Mnul Modeling One Dimensionl Solid Professor Rene Mtzdorf Universitet Kssel Stud. Mn. Rev 2.0 12/09 4. Modeling one-dimensionl solid There re two different wys to explin
More information1 Probability Density Functions
Lis Yn CS 9 Continuous Distributions Lecture Notes #9 July 6, 28 Bsed on chpter by Chris Piech So fr, ll rndom vribles we hve seen hve been discrete. In ll the cses we hve seen in CS 9, this ment tht our
More informationChapter 6 Notes, Larson/Hostetler 3e
Contents 6. Antiderivtives nd the Rules of Integrtion.......................... 6. Are nd the Definite Integrl.................................. 6.. Are............................................ 6. Reimnn
More informationMath 113 Exam 2 Practice
Mth 3 Exm Prctice Februry 8, 03 Exm will cover 7.4, 7.5, 7.7, 7.8, 8.-3 nd 8.5. Plese note tht integrtion skills lerned in erlier sections will still be needed for the mteril in 7.5, 7.8 nd chpter 8. This
More informationIn Section 5.3 we considered initial value problems for the linear second order equation. y.a/ C ˇy 0.a/ D k 1 (13.1.4)
678 Chpter 13 Boundry Vlue Problems for Second Order Ordinry Differentil Equtions 13.1 TWO-POINT BOUNDARY VALUE PROBLEMS In Section 5.3 we considered initil vlue problems for the liner second order eqution
More informationAppendix 3, Rises and runs, slopes and sums: tools from calculus
Appendi 3, Rises nd runs, slopes nd sums: tools from clculus Sometimes we will wnt to eplore how quntity chnges s condition is vried. Clculus ws invented to do just this. We certinly do not need the full
More informationAP Calculus Multiple Choice: BC Edition Solutions
AP Clculus Multiple Choice: BC Edition Solutions J. Slon Mrch 8, 04 ) 0 dx ( x) is A) B) C) D) E) Divergent This function inside the integrl hs verticl symptotes t x =, nd the integrl bounds contin this
More informationBernoulli Numbers Jeff Morton
Bernoulli Numbers Jeff Morton. We re interested in the opertor e t k d k t k, which is to sy k tk. Applying this to some function f E to get e t f d k k tk d k f f + d k k tk dk f, we note tht since f
More informationdx x x = 1 and + dx α x x α x = + dx α ˆx x x α = α ˆx α as required, in the last equality we used completeness relation +
Physics 5 Assignment #5 Solutions Due My 5, 009. -Dim Wvefunctions Wvefunctions ψ α nd φp p α re the wvefunctions of some stte α in position-spce nd momentum-spce, or position representtion nd momentum
More informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
More informationCHM Physical Chemistry I Chapter 1 - Supplementary Material
CHM 3410 - Physicl Chemistry I Chpter 1 - Supplementry Mteril For review of some bsic concepts in mth, see Atkins "Mthemticl Bckground 1 (pp 59-6), nd "Mthemticl Bckground " (pp 109-111). 1. Derivtion
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationQuantum Physics III (8.06) Spring 2005 Solution Set 5
Quntum Physics III (8.06 Spring 005 Solution Set 5 Mrch 8, 004. The frctionl quntum Hll effect (5 points As we increse the flux going through the solenoid, we increse the mgnetic field, nd thus the vector
More informationCS667 Lecture 6: Monte Carlo Integration 02/10/05
CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number of
More informationWeek 10: Line Integrals
Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.
More informationW. We shall do so one by one, starting with I 1, and we shall do it greedily, trying
Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)
More informationChapter 3 Polynomials
Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling
More informationSUMMER KNOWHOW STUDY AND LEARNING CENTRE
SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationContinuous Random Variables Class 5, Jeremy Orloff and Jonathan Bloom
Lerning Gols Continuous Rndom Vriles Clss 5, 8.05 Jeremy Orloff nd Jonthn Bloom. Know the definition of continuous rndom vrile. 2. Know the definition of the proility density function (pdf) nd cumultive
More informationDo the one-dimensional kinetic energy and momentum operators commute? If not, what operator does their commutator represent?
1 Problem 1 Do the one-dimensionl kinetic energy nd momentum opertors commute? If not, wht opertor does their commuttor represent? KE ˆ h m d ˆP i h d 1.1 Solution This question requires clculting the
More informationHomework 11. Andrew Ma November 30, sin x (1+x) (1+x)
Homewor Andrew M November 3, 4 Problem 9 Clim: Pf: + + d = d = sin b +b + sin (+) d sin (+) d using integrtion by prts. By pplying + d = lim b sin b +b + sin (+) d. Since limits to both sides, lim b sin
More informationThe practical version
Roerto s Notes on Integrl Clculus Chpter 4: Definite integrls nd the FTC Section 7 The Fundmentl Theorem of Clculus: The prcticl version Wht you need to know lredy: The theoreticl version of the FTC. Wht
More informationLecture 1: Introduction to integration theory and bounded variation
Lecture 1: Introduction to integrtion theory nd bounded vrition Wht is this course bout? Integrtion theory. The first question you might hve is why there is nything you need to lern bout integrtion. You
More informationQuantum Physics II (8.05) Fall 2013 Assignment 2
Quntum Physics II (8.05) Fll 2013 Assignment 2 Msschusetts Institute of Technology Physics Deprtment Due Fridy September 20, 2013 September 13, 2013 3:00 pm Suggested Reding Continued from lst week: 1.
More informationPhysics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2011
Physics 74 Grdute Quntum Mechnics Solutions to Finl Exm, Fll 0 You my use () clss notes, () former homeworks nd solutions (vilble online), (3) online routines, such s Clebsch, provided by me, or (4) ny
More informationState space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies
Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More information