arxiv: v1 [cond-mat.str-el] 4 Jul 2008
|
|
- Aron Phillips
- 5 years ago
- Views:
Transcription
1 Eistence of a metallic phase and upper bounds of the Hartree-Fock energy in the homogeneous electron gas F. Delyon, M. Duneau, B. Bernu, and M. Holzmann, 3 CPHT, UMR 7644 of CNRS, Ecole Polytechnique, Palaiseau, France LPTMC, UMR 76 of CNRS, Université P. et M. Curie, Paris, France 3 LPMMC, UMR 5493 of CNRS, Université J. Fourier, Grenoble, France Dated: July 8, 8 We consider the ground state of the homogeneous electron gas and we prove that a Hartree-Fock solution, motivated by previous simulations, has lower energy than the Fermi gas in the large density limit. This solution is a metallic phase : the density modulation corresponds to a partially occupied crystal the number of sites is larger than the number of electrons. I. INTRODUCTION arxiv:87.77v [cond-mat.str-el] 4 Jul 8 The homogeneous electron gas is one of the fundamental models in condensed matter physics. Despite its simplicity - the system consists of electrons interacting with each other through a /r potential to which a uniform positive background is added for charge neutrality - the phase diagram at zero temperature is nontrivial,,3. In general, it is given in terms of the dimensionless parameter r s = α D a B n /D, where D is the space dimension, n is the electronic density, a B the Bohr radius and αd D is the volume of the sphere unity. At r s =, the ground state is the Fermi gas. At large r s, the Hartree-Fock ground state of the electron gas is a Wigner crystal, that is a state where the charge density forms a triangular crystal with eactly one electron per lattice site. In the intermediate region the Hartree-Fock approimation is not relevant and more sophisticated methods show that the ground state of the Hamiltonian is quite different. At small r s, the ground state of the electron gas within Hartree-Fock is still not known 8. Although the ground state of the ideal Fermi gas r s = remains an eigenstate of the Hartree-Fock Hamiltonian of the electron gas at any density, already Wigner 5 argued that the unpolarized Fermi gas is unstable even in the limit r s. Later, Overhauser showed the instability of the unpolarized Fermi gas with respect to spin-density waves within the Hartree- Fock approimation 6. Only recently, a Hartree-Fock study of the unpolarized three-dimensional electron gas was performed which proposes a more complicated structure of a ground state with spin-density waves 9. In Ref., we describe the results of our numerical simulations of the two-dimensional D, polarized electron gas at small r s, and show that the ground state is neither a Fermi gas nor a Wigner crystal: the charge density modulation we find corresponds to a partially occupied crystal the number of sites is larger than the number N of electrons. We refer to this solution as a metallic phase. The charge density of this metallic phase represents a triangular lattice with reciprocal generators Q i of modulus k F. This modulation is mainly carried by the wave vectors close to the Fermi surface. In order to observe the metallic phase, the number of electrons N has to eceed a threshold ranging from N > at r s.6 up to N > at r s. As we shall see below, the threshold increases eponentially as r s goes to which may eplain why this metallic phase has not been observed in previous simulations. In this study, we obtain rigorous upper bounds on the energy of the metallic phase in the two-dimensional polarized case with a class of states mimicking the real states obtained numerically. These bounds are obtained in the limit r s where the calculation is simplified by the long range behavior of the interaction potential. However, our numerical simulations show that this metallic phase also eists for screened potential, but strictly speaking we cannot prove in this case the eistence of such a metallic phase in the limit r s. This bounds are easily etended to the unpolarized D case. Then we show how this proof etends to the 3D case, whereas we have no numerical simulations indicating that this phase may correspond to the actual ground state of the 3D electron gas in Hartree-Fock. Finally, let us point out that our bounds for instance Eq. 4 must be considered as mathematical bounds; they only give an understanding of the behavior of the energy of the ground state, but the actual constants involved in the epressions have to be evaluated by other means. II. NUMERICAL RESULTS In a previous paper, we have computed D Hartree-Fock states of lower energy than the Fermi gas for values of r s about. Such a state Ψ is obtained as a Slater determinant N i= ψ i where N is the number of electrons and {ψ i } is a set of orthonormal vectors. Only the space generated by the single particle wavefunctions ψ i s is relevant, and in order to
2 b k N=499 r s = k y / k k / k FIG. : Numerical values of b k for 499 electrons in two dimensions at r s =. understand our numerical results we need to choose a canonical representation of the ψ i s. The natural choice is to use the natural base of the Fermi gas {φ i } i=...n corresponding to some indeation of the plane waves of the Fermi sphere. Let M be the square matri defined by M ij = φ i ψ j. The Singular Value Decomposition SVD of M is M = UσV where U and V are unitary matrices and σ is a diagonal positive matri. Let us define φ i = j U ij φ j ψ i = j V ji ψ j. Then {φ i } and {ψ i } are also two sets of N orthonormal vectors satisfying φ i ψ j = δij σ i. They are the nearest bases of Span{φ i } and Span{ψ i } and they define a canonical unitary operator W from Span{φ i } onto Span{ψ i } by: W i λ i φ i = i λ iψ i 3 In particular, we have: W φ i = W j U ij φ j 4 = j = j = j U ij ψ j 5 U ij k V kj ψ k 6 UV ik ψ k 7 Thus the natural basis of Span{ψ i } associated to the basis {φ i } i=...n is {W φ i } i=...n. Numerically we have chosen the ψ i s in this way and it appears that, if r s is not too large, ψ i is close to φ i at least for i associated to a wave vector not too close to the Fermi surface. Thus the largest amplitude of ψ i, in the k space, is for the k i corresponding to φ i. Fig. represents the net largest amplitude of ψ i, that we denote b ki, for 499 electrons at r s = in D. As r s decreases, the other amplitudes of the ψ i s become very small. Furthermore, the wave vector k i corresponding to b ki satisfies k i k i = Q α where the generators {Q α } α=...6 belong to the si-fold star of a triangular lattice.
3 Notice that the modulus of b k is maimal for k on a si-fold star and close to the Fermi surface. In the direct space, the Q α s induce a modulation corresponding to a triangular lattice for the charge density. But while this lattice forms a Wigner crystal with one electron per site for r s >.6, for smaller r s it is a denser lattice corresponding to a crystal with an incomplete band filling. Indeed, at any r s the electrostatic interaction favors a periodic distribution of charges and at small r s only k-vectors close to Fermi surface can be modulated. The optimal solution is to choose Q α = k F which is larger than the reciprocal lattice vector of the Wigner case. In the following, we consider analytic solutions analogous to those of our numerical results and we prove that their energies are lower than the energies of the Fermi gas as r s goes to. 3 III. ENERGY OF THE FERMI GAS: POLARIZED CASE We consider the Hamiltonian of N electrons in a D or 3D square bo of volume Ω with periodic boundary conditions. H = h m + e V 8 where V is the -body Coulomb potential i j / r i r j, the electron mass is m, and e is its charge. It is convenient to choose Hartree as the unit of energy, Ha = h /ma B, where a B = h /me is the Bohr radius. We get: H = a B + a B V 9 Let ψ n be an orthonormalized set of N vectors of L Ω. They define the N-particle Slater determinant Ψ = n ψ n. And the energy of Ψ is: E = Ψ H Ψ = a B ψ n ψ n + ψ n ψ n v ψ n ψ n a n B n,n where v is defined as: ϕ ϕ v ψ ψ = d dy ϕ ϕ y y ψ ψ y. In order to avoid problems due to the Coulomb singularity, we introduce the jellium model and define the potential acting on the plane waves φ k as: φ k φ k v φ k+q φ k q = π Ω D q for q and otherwise, so that the total charge of the electrons is compensated by a positive background charge. The Fermi gas is defined by Φ = k <k F φ k where α D kf D = πd N/Ω and αd D is the volume of the unit sphere. E F G = Φ H Φ = a B k D π a B Ω k <k F k, k <k F k k D 3 As Ω goes to with Ω/N fied, the thermodynamic limit for the energy per particle is obtained by the substitution k Ω π D dk: E F G N = a B = a B Ω Nπ D Ω Nπ D kd+ F dk k k <k F dk k k < a B π D a B k F π D dkdk k, k <k F k k D dkdk k k D k, k < 4
4 4 k y S S k+q k k' S k S 3 k S 4 S 5 FIG. : The circle is the Fermi surface. The shaded surfaces are the regions where bk is nonzero. The new state ψ k miing φ k and φ k+qk is now resonant with ψ k. From the definition of r s = α D a B n /D and k F, it follows that k F αd r sa B = π. Thus, we have: E F G N = π α D+4 dk k r sα D k < 4π D dkdk k, k < k k D which gives for D = α = π: D r s E F G N = πr s dk k r s k < 4π k, k < dkdk k k 5 6 IV. HARTREE-FOCK UPPER BOUNDS: POLARIZED D CASE We restrict ourself to the D polarized case and we want to estimate the energy for a class of states inspired by our numerical results. Let us consider a state Ψ = k <k F ψ k where: ψ k = a k φ k + b k φ k+qk with Q k in { k F cos pπ/3, sin pπ/3} p=...5. For k = k cos θ, sin θ we choose Q k such that k + Q k is minimal; that is, we choose p as the integer part of 3θ/π + / and we must assume b k is zero if k is zero or θ = π/6 + nπ/3. Furthermore, we assume that a k and b k are real positive number and invariant thru the rotation of nπ/6 and the symmetry θ θ i.e. the dihedral group D 6. The ψ k s are normalized, so that a k + b k = and b k = if k Q k < kf ɛ i.e. b k is not zero only in the vicinity of {k F cos pπ/3, sin pπ/3} p=...5, see Fig.. Thus the limit energy per particle is given by: E N = πrs dk ψ k ψ k + r s k < 4π k, k < dkdk where, as in 6, the k s have been renormalized by k F and thus Q k =. We define E by: Then E = k < E E F G N Ω π ψ k ψ k v ψ k ψ k 7 = πrs E 8 dk [ ψ k ψ k k ] + r s 4π E V 9
5 5 where E V = k, k < dkdk Ω π ψ k ψ k v ψ k ψ k + k k A. Potential energy contribution: E V Setting v q = / q : Eq., may be divided into 4 parts: Ω π ψ k ψ k v ψ k ψ k + v k k = v k k v k k Q k b k a k {b k =, b k = }: the contribution is zero. + v k k v k+qk k b ka k {b k =, b k }, {b k, b k = }: both cases are equivalent. For {b k =, b k }, the integrant of Eq. is: + v k k v k+qk Q k k b kb k + v Qk a k b k a k b k δ Qk +Q k + δ Qk Q k v k k a k b k a k b k δ Qk Q k v k+qk k + v k k Q k a k a k b k b k δ Qk +Q k Ω π ψ k φ k, v ψ k φ k + k k = v k k v k+q k k b k Let S be the sector of unit disk between π/6 and π/6 see Fig. ; then in this sector Q k =, and by symmetry: dkdk v k k v k+qk k b k = 6 dkdk v k k v k+qk k b k 3 b k = = 6 k S,b k = k S,b k = dkdk v k k v k k b k 4 Cɛ dkdk v k k v k k b k S, k < ɛ k 5 where k = k, k y. In S, k = k, k y where k is close to and setting k =, we assume from now that b k = b/ɛ. In Appendi A we prove that: {b k, b k }: k S, k < ɛ dkdk v k k v k+qk k b k 8ɛ ɛ [ ln ɛ + O ] db 6 By symmetry we can assume that k belongs to S. If k S S 3 all the v appearing in are uniformly bounded. And since the k-volume for each sector goes like ɛ ɛ, the contribution of these terms is bounded by Cɛ 3. In the same way v k k is bounded when k S 3 and v k+qk k is bounded when k S. Thus setting: f := a kb k + b ka k a ka k b k b k = a k b k b k a k 7 g := a kb k + b ka k + a ka k b k b k = a k b k + b k a k 8
6 one can check that: b k,b k In Appendi B we prove that Ω π ψ k ψ k v ψ k ψ k + v k k Cɛ dkdk v k k f v k+qk +k g 9 k,k S k,k S dkdk v k k f v k+qk +k g 4ɛ ɛ [ ln ɛ + O ] d d fɛ, ɛ gɛ, ɛ 3 6 Thus, summing the four contribution gives: E V Cɛ 3 + 6ɛ ɛ [ ln ɛ + O ] d 6b + 4 d f, g, 3 B. Kinetic energy contribution: The variation of the kinetic energy is given by: dk [ ψ k ψ k k ] = 6 dk ψ k ψ k k 3 k < k S = 6 ɛ 6 8ɛ ɛ d y m 4b /ɛ 33 d b 34 C. Total energy: Inserting Eqs.3,34 in Eq. 9, the variation of the total energy from the Fermi gas energy becomes: E 6ɛ ɛ d 8b + r s [ ln ɛ + O ] 6b + 4 d f, g, 4π = 6 8ɛ ɛ d b + r s [ ln ɛ + O ] b ab d a b π Let us set 35 Then δ = ɛ ɛ I = I = 5π 36 d b 37 d b + ab d a b 38 E 6 8 δ [ I r s I ln δ + O ] 39 If I >, as r s goes to, E is minimal in Eq. 39 for δ defined by: δ min = e ep I I r s 4 and finally inserting δ min in Eq. 39 gives: E < 6 8 e ep I r s I 4 I r s
7 We now have to find a solution b such that I is positive. Choosing b = b or b = b leads to negative I. In the Appendi C, as r s goes to we find a family of b leading to : E < r s ep 5π + O 4 3r s rs Though such a bound is correct in the thermodynamic limit, this behavior in not so relevant for finite systems. Indeed, numerical systems consider about 3 electrons and, as we shall see later, the Fermi gas becomes then the ground state for r s <. Thus the asymptotic bound 4 is not very helpful in the real word. Nevertheless, for realistic r s, on can choose a suitable function b and evaluate numerically I and I. For instance, with b = b η as in 88 of Appendi C and η =. we get E <.6 4 r s ep r s 7 V. HF UPPER BOUNDS: UNPOLARIZED D CASE We consider now the spin of the electrons, and we restrict to the case where half the electrons have a spin up and the others have a spin down. We can choose a solution as the product of Slaters Ψ + and Ψ for the spins up and down. a B ψ + n ψ n + + ψ + a n ψ + n v ψ+ n ψ + n B E NP = n n,n ψ n ψ n + ψ a n ψ n v ψ n ψ n n B n,n + ψ + a n ψ n v ψ+ n ψ n + B n,n n,n ψ n ψ + n v ψ n ψ + n 46 The first terms 44 and 45 can be evaluated as in 8 noticing that the definition of k F becomes k F = 4πN ±/Ω = πn/ω. So the r s in the evaluation of E is now r s and: E NP E NP F G N = πrs Er s The etra terms 46 coming from the direct potential provide contributions like a + n b + n a n b n times a positive factor. These contributions are regular as r s goes to, and does not modify the asymptotic energy. Nevertheless, these contributions depends on the signs of b ± n in the two Slaters Ψ + and Ψ. Thus, since the energy of one Slater does not depend on the global sign of the b ± n s, the minimal state is obtained by choosing b n = b + n. For such a state, the main difference with the polarized case is that the charge density becomes flat, though the spin density is not. Finally, we can compare our solutions with previous solutions proposed by Fedders and Martin 7. They consider a more complicated modulation of the states near the center of Fermi sphere. While their results may be correct, they missed the dominant contribution coming from the states near the Fermi surface. Indeed, our computations give nonzero contribution for the states near the center, but for r s =.5 b k is about 5 and the energy benefit, following their formula, must be of order of the energy benefit of the eternal states. VI. THE METALLIC PHASE IN FINITE SYSTEMS For finite systems of N electrons, the minimum requirement is that the surface k Q k > kf ɛ, so that it contains at least one plane wave of the finite, discrete system. This gives the condition Nɛ ɛ > and from Eqs. 36,4, this leads to: 3I N > ep 48 5I r s
8 8 With our approimate behavior for I /I in two dimensions, we find: 3π N > ep r s 49 i.e. N > 5 for r s =.8. This bound is compatible with our numerical simulations where the metallic phase disappears at rs = for N = 5. Furthermore, this may eplain why this metallic phase has not been observed in previous simulations. VII. HARTREE-FOCK UPPER BOUNDS: POLARIZED 3D CASE Mutatis mutandis, the 3D case works in the same way. The potential term v k is now / k but the dominant terms are provide now from integrals like: Ik, k = dk y dk z dk ydk z k k 5 where k is close to k. Let us suppose that k > k, then k, k k k Ik, k = µdu k k + u 5 where u stands for k y k y + k z k z / and µdu stands for the distribution of u. Here we are interested in small values of u. One can check that µu, u + ɛ π k ɛu provided that u is sufficiently small k y, k z may be any point in disk of radius k and k y, k z may belong to a thin ring around k y, k z Ik, k π k ɛ udu k k + u 5 π k lnk k 53 This must be compared with the D case Ik, k k lnk k. Mutatis mutandis, Eq. 9 becomes: E = dk ψ k ψ k k + κr s E V 54 k < where κ = α 3/4π 3 and α 3 3 = 4π/3. Now, in E V E V Cɛ 4 + 6zπ ɛ 3 [ ln ɛ + O ] the ɛ factors have to be replaced by π ɛ: d b ab d a b 55 where z is the number of sectors of the Fermi sphere. The variation of the kinetic energy is now Thus setting I = k < I = κπ 3 dk ψ k ψ k k = z8πɛ 3 d b 56 d b 57 d b + ab E z8πɛ 3 I 3r s I [ ln ɛ + O ] d a b 58
9 Once more, if we find a b such that I is positive, we have: E < z8π e ep I r s I 59 I r s The operator A of Appendi C is now modified. But the same analysis show that the spectrum of A is still, 4, so that the asymptotic bound for small r s is: E < z8π e ep + o/r s r s 6 6κπr s 9 VIII. CONCLUSION We have proven upper bounds for the energy of a metallic state below the Fermi gas energy. To our knowledge, it is the first time that rigorous upper bounds for the ground state energy of the polarized electron gas are obtained going beyond the simple Fermi gas states. The main interest of this proof is that it focuses on the small r s region where the Hartree-Fock approimation is relevant and that the upper bounds were obtained using metallic states which correspond to solutions obtained by numerical simulations. In order to judge their relevance for the true ground state of the electron gas in the high density region, we have to consider correlation effects beyond the Hartree-Fock approimation. Due to the long-range, singular behavior of the Coulomb-potential, the perturbation epansion has to be rearranged which amounts to an effective screening of the electron interaction. Unfortunately a rigorous etension of our proof to the case of a screened Coulomb potential is not straightforward. However, numerical calculations provide evidence for the stability of the metallic states in the D polarized system within Hartree-Fock, so that these states are relevant candidates for the true ground state of the two-dimensional electron gas. IX. APPENDIX A We have to estimate If = dkdk v k, k k k v k k f k 6 k <ɛ, k < ɛ where k = k, k y and f is a positive function. where y m = dk v k k v k k = dk asinh y m k y k k + asinh y m + k y k k k. And since asinh asinh y ln /y for > y > : dk v k k v k k ɛ +ɛ asinh y m k y k k asinh y m + k y k k dk ln k k k k 6 63 We set k = and y m =, If = 4 4 ɛ ɛ ɛ df y m ɛ +ɛ df y m ɛ ɛ df y m du ɛ ɛ dk ln + k k du ln u + u u = 8ɛ ɛ [ ln ɛ + O ] dfɛ 67
10 X. APPENDIX B We have to estimate If, g = dkdk k, k vk k f v k+k g k, k <ɛ 68 where k = k, k y and f and g are positive functions of k and k. Setting k =, k =, k y = y, k y = y, and r ± = ± + y y, Eq.68 can be rewritten: ɛ ɛ If, g = d d dydy f g r r + ɛ ɛ = d d dydy f g r r + 69 where y and y must satisfy + y and + y. Since asinh ln +, the first term in 69 is bounded by: ɛ ɛ d d dydy ɛ ɛ ym f = d d f dy asinh y m + y r y m + asinh y m y 7 4 ɛ ɛ d d ɛ ɛ d f ym y m dy asinh y m 7 d fy m ln + 4y m 7 4ɛ ɛ [ lnɛ + O ] d d fɛ, ɛ 73 On the other hand, using asinh ln, the last term of 69 is: ɛ ɛ d d g dydy ɛ ɛ ym = d d g dy asinh y m y r + y m + + asinh y m + y + 74 ɛ ɛ ym d d g dy ln 4 y m y y m + 75 And we have: 4 ɛ d ɛ d gy m [ ln ɛ + O ] 76 4ɛ ɛ [ ln ɛ + O ] d If, g 4ɛ ɛ [ ln ɛ + O ] d d gɛ, ɛ 77 d fɛ, ɛ gɛ, ɛ 78 XI. APPENDIX C Here we provide eact bounds on I and I given by 37, 57. In order to estimate I we introduce the linear operator A : Af = fydy + defined on the Hilbert space of the functions on [, ] with the scalar product: f g = fy ydy 79 fgd 8
11 Then A is a bounded symmetric operator and: I /I = 5π ab Aab b 8 The unitary operator f gy = fe y e 5/4y from L [, ], d onto L [, + ], d maps the operator A onto the operator Ã: Then Ãg = e /4 Ãe ik = e y/4 gydy + e/4 + e y/4 gydy 8 4/6 + k eik / + ik e /4 83 Thus setting: g k = {g k } k> is a full set of pseudo-eigenvectors satisfying: [ + i4ke ik i4ke ik] 84 + i4k Ãg k = 4/6 + k g k Thus the spectrum of à is, 4 and the spectral measure is purely absolutely continuous; the largest spectral value is 4 with a pseudo-eigenvector g 4 = + 4 corresponding to f 4 = 5/4 4 ln. But f 4 is infinite and f 4 diverges at. The net step is to choose a family of functions b η such that a η = b η is defined and a η b η Aa η b η / b η is close to 4. Thus setting f η = minf 4, f 4 η for < η, we have: f η = [ ln 3 η 66 ] 3 5 ln η + Oln η 85 f η Af η = 4 [ ln 3 η 4 ] 3 5 ln η + Oln η 86 Then: f η Af η f η = 4 ln η + Oln η 87 Thus f η is a good candidate for the linear part of the problem. Now, by the simple scaling: b η = f η fη η 88 we get the nonlinear candidate satisfying b η /, a η = b η is well defined, a η / and b η satisfies 87. We must now estimate the simultaneous convergence of I/I 8 and I as η decreases. b η Ab η b η a η Ab η a η = b η b η a η A b η b η a η + b η b η a η A b η b η b η a η A b η = 8 b η b η a η b η + b η b η a η A 4b η 8 b η b η a η b η + b η b η a η A 4b η 8 b η b η a η b η + 8 b η b η a η b η A 4b η
12 where: b η b η a η b η = b η [ a η ] d b η [ a η ] d sup a η b η b η a η and since b η b η a η > b η b η a η = b η b η a η b η b η a η b η a η b η b η a η = b η By direct computation: for η small enough, and thus: b η b η 6 5 ln η b η Ab η b η b ηa η Ab η a η b η 8 6 ln η And finally, from 87 for b η and 89 and I = b η, 8 gives: I /I 3 C + Oln η 9 5π ln η I [ ] 5π η5/ ln 3 η + O 9 ln η where C = Choosing η sufficiently small, this proves that E is strictly negative for any r s > ; furthermore choosing η to minimize E Eq. 4, i.e. 9r s ln η = Cπ, we obtain as r s goes to : E < r s ep 5π + O 93 3r s rs B. Tanatar and D. M. Ceperley, Phys. Rev. B 39, ; C. Attaccalite, S. Moroni, P. Gori-Giorgi, and G. B. Bachelet, Phys. Rev. Lett. 88, 566. B. Bernu, L. Cândido, and D. M. Ceperley, Phys. Rev. Lett. 86, H. Falakshahi and X. Waintal, Phys. Rev. Lett. 94, 468 5, X. Waintal, Phys. Rev. B 73, J. R. Trail, M. D. Towler, and R. J. Needs, Phys. Rev. B 68, E. P. Wigner, Trans. Faraday Soc. 34, ; Phys. Rev. 46, A. W. Overhauser, Phys. Rev. Lett. 4, 46 96; Phys. Rev. 8, A. Fedders and P. C. Martin, Phys. Rev. 43, G. F. Giuliani and G. Vignale, Quantum Theory of the Electron Liquid, Cambridge University Press, Cambridge 5. 9 S. Zhang and D. M. Ceperley, Phys. Rev. Lett., , arxiv: B. Bernu, F. Delyon, M. Duneau, and M. Holzmann, cond-mat. arxiv:84.5
arxiv: v2 [cond-mat.str-el] 16 Sep 2013
Hartree-Fock Ground State Phase Diagram of Jellium arxiv:1307.3081v2 [cond-mat.str-el] 16 Sep 2013 L. Baguet, 1 F. Delyon, 2 B. Bernu, 1 and. Holzmann 1, 3 1 LPTC, UR 7600 of CNRS, Université P. et. Curie,
More informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
More informationIs the homogeneous electron gas homogeneous?
Is the homogeneous electron gas homogeneous? Electron gas (jellium): simplest way to view a metal homogeneous and normal Hartree-Fock: simplest method for many-electron systems a single Slater determinant
More informationThe Overhauser Instability
The Overhauser Instability Zoltán Radnai and Richard Needs TCM Group ESDG Talk 14th February 2007 Typeset by FoilTEX Introduction Hartree-Fock theory and Homogeneous Electron Gas Noncollinear spins and
More informationarxiv: v1 [cond-mat.str-el] 15 May 2011
arxiv:1105.2964v1 cond-mat.str-el] 15 May 2011 Finite-size analysis of the Fermi liquid properties of the homogeneous electron gas Markus Holzmann a,b, Bernard Bernu a, and David M. Ceperley c a LPTMC,
More informationCondensed matter theory Lecture notes and problem sets 2012/2013
Condensed matter theory Lecture notes and problem sets 2012/2013 Dmitri Ivanov Recommended books and lecture notes: [AM] N. W. Ashcroft and N. D. Mermin, Solid State Physics. [Mar] M. P. Marder, Condensed
More informationQuantum Dynamics. March 10, 2017
Quantum Dynamics March 0, 07 As in classical mechanics, time is a parameter in quantum mechanics. It is distinct from space in the sense that, while we have Hermitian operators, X, for position and therefore
More informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationV. ELECTRON GAS, RANDOM PHASE APPROXIMATION
V. ELECTRON GAS, RANDOM PHASE APPROXIMATION Marus Holzmann LPMMC, Maison de Magistère, Grenoble, and LPTMC, Jussieu, Paris marus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/marus/cours.html (Dated:
More informationdoes not change the dynamics of the system, i.e. that it leaves the Schrödinger equation invariant,
FYST5 Quantum Mechanics II 9..212 1. intermediate eam (1. välikoe): 4 problems, 4 hours 1. As you remember, the Hamilton operator for a charged particle interacting with an electromagentic field can be
More informationarxiv:quant-ph/ v5 10 Feb 2003
Quantum entanglement of identical particles Yu Shi Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, United Kingdom and Theory of
More informationarxiv: v1 [cond-mat.other] 20 Apr 2010
Characterization of 3d topological insulators by 2d invariants Rahul Roy Rudolf Peierls Centre for Theoretical Physics, 1 Keble Road, Oxford, OX1 3NP, UK arxiv:1004.3507v1 [cond-mat.other] 20 Apr 2010
More informationThe Quantum Heisenberg Ferromagnet
The Quantum Heisenberg Ferromagnet Soon after Schrödinger discovered the wave equation of quantum mechanics, Heisenberg and Dirac developed the first successful quantum theory of ferromagnetism W. Heisenberg,
More informationPhysics 127c: Statistical Mechanics. Weakly Interacting Fermi Gas. The Electron Gas
Physics 7c: Statistical Mechanics Wealy Interacting Fermi Gas Unlie the Boson case, there is usually no ualitative change in behavior going from the noninteracting to the wealy interacting Fermi gas for
More informationPhysics 443, Solutions to PS 4
Physics, Solutions to PS. Neutrino Oscillations a Energy eigenvalues and eigenvectors The eigenvalues of H are E E + A, and E E A and the eigenvectors are ν, ν And ν ν b Similarity transformation S S ν
More informationInteracting Fermi Gases
Interacting Fermi Gases Mike Hermele (Dated: February 11, 010) Notes on Interacting Fermi Gas for Physics 7450, Spring 010 I. FERMI GAS WITH DELTA-FUNCTION INTERACTION Since it is easier to illustrate
More information1 Fluctuations of the number of particles in a Bose-Einstein condensate
Exam of Quantum Fluids M1 ICFP 217-218 Alice Sinatra and Alexander Evrard The exam consists of two independant exercises. The duration is 3 hours. 1 Fluctuations of the number of particles in a Bose-Einstein
More informationMean-Field Theory: HF and BCS
Mean-Field Theory: HF and BCS Erik Koch Institute for Advanced Simulation, Jülich N (x) = p N! ' (x ) ' 2 (x ) ' N (x ) ' (x 2 ) ' 2 (x 2 ) ' N (x 2 )........ ' (x N ) ' 2 (x N ) ' N (x N ) Slater determinant
More information18.303: Introduction to Green s functions and operator inverses
8.33: Introduction to Green s functions and operator inverses S. G. Johnson October 9, 2 Abstract In analogy with the inverse A of a matri A, we try to construct an analogous inverse  of differential
More informationFloquet Topological Insulator:
Floquet Topological Insulator: Understanding Floquet topological insulator in semiconductor quantum wells by Lindner et al. Condensed Matter Journal Club Caltech February 12 2014 Motivation Motivation
More informationPhase Diagram of the Two-Dimensional Homogeneous Electron Gas
Phase Diagram of the Two-Dimensional Homogeneous Electron Gas Neil D. Drummond and Richard J. Needs TCM Group, Cavendish Laboratory, University of Cambridge ESDG Meeting Wednesday 13th February, 2008 Two-Dimensional
More informationLight-Cone Quantization of Electrodynamics
Light-Cone Quantization of Electrodynamics David G. Robertson Department of Physics, The Ohio State University Columbus, OH 43210 Abstract Light-cone quantization of (3+1)-dimensional electrodynamics is
More informationSelf-consistent Field
Chapter 6 Self-consistent Field A way to solve a system of many electrons is to consider each electron under the electrostatic field generated by all other electrons. The many-body problem is thus reduced
More informationGeneralization of the matrix product ansatz for integrable chains
arxiv:cond-mat/0608177v1 [cond-mat.str-el] 7 Aug 006 Generalization of the matrix product ansatz for integrable chains F. C. Alcaraz, M. J. Lazo Instituto de Física de São Carlos, Universidade de São Paulo,
More informationMultidimensional partitions of unity and Gaussian terrains
and Gaussian terrains Richard A. Bale, Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureu ABSTRACT Partitions of unity play an important rôle as amplitude-preserving windows for nonstationary
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationJournal of Theoretical Physics
1 Journal of Theoretical Physics Founded and Edited by M. Apostol 53 (2000) ISSN 1453-4428 Ionization potential for metallic clusters L. C. Cune and M. Apostol Department of Theoretical Physics, Institute
More informationSpontaneous Spin Polarization in Quantum Wires
Spontaneous Spin Polarization in Quantum Wires Julia S. Meyer The Ohio State University with A.D. Klironomos K.A. Matveev 1 Why ask this question at all GaAs/AlGaAs heterostucture 2D electron gas Quantum
More informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More informationQuantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid
Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid N. D. Drummond and R. J. Needs TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge
More informationLecture 5 (Sep. 20, 2017)
Lecture 5 8.321 Quantum Theory I, Fall 2017 22 Lecture 5 (Sep. 20, 2017) 5.1 The Position Operator In the last class, we talked about operators with a continuous spectrum. A prime eample is the position
More informationarxiv:quant-ph/ v2 24 Dec 2003
Quantum Entanglement in Heisenberg Antiferromagnets V. Subrahmanyam Department of Physics, Indian Institute of Technology, Kanpur, India. arxiv:quant-ph/0309004 v2 24 Dec 2003 Entanglement sharing among
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationElectrons in periodic potential
Chapter 9 Electrons in periodic potential The computation of electronic states in a solid is a nontrivial problem. A great simplification can be achieved if the solid is a crystal, i.e. if it can be described
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationChapter 4 (Lecture 6-7) Schrodinger equation for some simple systems Table: List of various one dimensional potentials System Physical correspondence
V, E, Chapter (Lecture 6-7) Schrodinger equation for some simple systems Table: List of various one dimensional potentials System Physical correspondence Potential Total Energies and Probability density
More informationWe already came across a form of indistinguishably in the canonical partition function: V N Q =
Bosons en fermions Indistinguishability We already came across a form of indistinguishably in the canonical partition function: for distinguishable particles Q = Λ 3N βe p r, r 2,..., r N ))dτ dτ 2...
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz November 29, 2010 The self interaction error and its correction Perdew-Zunger SIC Average-density approximation Weighted density approximation
More informationFoundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
More informationGiven the vectors u, v, w and real numbers α, β, γ. Calculate vector a, which is equal to the linear combination α u + β v + γ w.
Selected problems from the tetbook J. Neustupa, S. Kračmar: Sbírka příkladů z Matematiky I Problems in Mathematics I I. LINEAR ALGEBRA I.. Vectors, vector spaces Given the vectors u, v, w and real numbers
More informationProperties of Matrices and Operations on Matrices
Properties of Matrices and Operations on Matrices A common data structure for statistical analysis is a rectangular array or matris. Rows represent individual observational units, or just observations,
More informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
More informationPhysics 221B Spring 2018 Notes 30 The Thomas-Fermi Model
Copyright c 217 by Robert G. Littlejohn Physics 221B Spring 218 Notes 3 The Thomas-Fermi Model 1. Introduction The Thomas-Fermi model is a relatively crude model of multi-electron atoms that is useful
More informationLOCAL MOMENTS NEAR THE METAL-INSULATOR TRANSITION
LOCAL MOMENTS NEAR THE METAL-INSULATOR TRANSITION Subir Sachdev Center for Theoretical Physics, P.O. Box 6666 Yale University, New Haven, CT 06511 This paper reviews recent progress in understanding the
More informationFermi liquid theory Can we apply the free fermion approximation to a real metal? Phys540.nb Strong interaction vs.
Phys540.nb 7 Fermi liquid theory.. Can we apply the free fermion approximation to a real metal? Can we ignore interactions in a real metal? Experiments says yes (free fermion models work very well), but
More informationA matrix over a field F is a rectangular array of elements from F. The symbol
Chapter MATRICES Matrix arithmetic A matrix over a field F is a rectangular array of elements from F The symbol M m n (F ) denotes the collection of all m n matrices over F Matrices will usually be denoted
More informationPhysics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension
Physics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension In these notes we examine Bloch s theorem and band structure in problems with periodic potentials, as a part of our survey
More informationPhonons and lattice dynamics
Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed
More informationQuantum Cluster Methods (CPT/CDMFT)
Quantum Cluster Methods (CPT/CDMFT) David Sénéchal Département de physique Université de Sherbrooke Sherbrooke (Québec) Canada Autumn School on Correlated Electrons Forschungszentrum Jülich, Sept. 24,
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationIn particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with
Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient
More information7.4. Why we have two different types of materials: conductors and insulators?
Phys463.nb 55 7.3.5. Folding, Reduced Brillouin zone and extended Brillouin zone for free particles without lattices In the presence of a lattice, we can also unfold the extended Brillouin zone to get
More informationQuantum Chaos and Nonunitary Dynamics
Quantum Chaos and Nonunitary Dynamics Karol Życzkowski in collaboration with W. Bruzda, V. Cappellini, H.-J. Sommers, M. Smaczyński Phys. Lett. A 373, 320 (2009) Institute of Physics, Jagiellonian University,
More informationFourier Analysis Fourier Series C H A P T E R 1 1
C H A P T E R Fourier Analysis 474 This chapter on Fourier analysis covers three broad areas: Fourier series in Secs...4, more general orthonormal series called Sturm iouville epansions in Secs..5 and.6
More informationMomentum Distribution and Effective Mass of Jellium and Simple Metals. Markus Holzmann. LPTMC, CNRS, Universite Pierre et Marie Curie, Paris
Momentum Distribution and Effective Mass of Jellium and Simple Metals Markus Holzmann LPTMC, CNRS, Universite Pierre et Marie Curie, Paris LPMMC, CNRS, Université Joseph Fourier, Grenoble, France QMC:
More informationLinear Algebra and Dirac Notation, Pt. 3
Linear Algebra and Dirac Notation, Pt. 3 PHYS 500 - Southern Illinois University February 1, 2017 PHYS 500 - Southern Illinois University Linear Algebra and Dirac Notation, Pt. 3 February 1, 2017 1 / 16
More informationMany-body wavefunctions for normal liquid 3 He
Many-body wavefunctions for normal liquid 3 He Markus Holzmann, 1 Bernard Bernu, 1 and D. M. Ceperley 2 1 LPTMC, UMR 7600 of CNRS, Université Pierre et Marie Curie, Paris, France 2 Department of Physics
More informationThe Large-Spin Asymptotics of the Ferromagnetic XXZ Chain
Markov Processes Relat. Fields, 237 266 (2005) Markov Processes and M &RF P Related Fields c Polymat, Moscow 2005 The Large-Spin Asymptotics of the Ferromagnetic XXZ Chain T. Michoel and B. Nachtergaele
More informationarxiv: v2 [cond-mat.str-el] 20 Apr 2015
Gauging time reversal symmetry in tensor network states ie Chen, 2 and Ashvin Vishwanath 2 Department of Physics and Institute for Quantum Information and Matter, California Institute of echnology, Pasadena,
More informationAlgebraic Theory of Entanglement
Algebraic Theory of (arxiv: 1205.2882) 1 (in collaboration with T.R. Govindarajan, A. Queiroz and A.F. Reyes-Lega) 1 Physics Department, Syracuse University, Syracuse, N.Y. and The Institute of Mathematical
More informationNext topic: Quantum Field Theories for Quantum Many-Particle Systems; or "Second Quantization"
Next topic: Quantum Field Theories for Quantum Many-Particle Systems; or "Second Quantization" Outline 1 Bosons and Fermions 2 N-particle wave functions ("first quantization" 3 The method of quantized
More informationCanonical lossless state-space systems: staircase forms and the Schur algorithm
Canonical lossless state-space systems: staircase forms and the Schur algorithm Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics School of Mathematical Sciences Projet APICS Universiteit
More informationLinear Algebra Review. Vectors
Linear Algebra Review 9/4/7 Linear Algebra Review By Tim K. Marks UCSD Borrows heavily from: Jana Kosecka http://cs.gmu.edu/~kosecka/cs682.html Virginia de Sa (UCSD) Cogsci 8F Linear Algebra review Vectors
More informationV. FERMI LIQUID THEORY
V. FERMI LIQUID THEORY Markus Holzmann LPMMC, Maison de Magistère, Grenoble, and LPTMC, Jussieu, Paris markus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/markus/cours.html (Dated: March 30, 2010)
More informationPseudo-Hermitian eigenvalue equations in linear-response electronic-structure theory
1/11 Pseudo-Hermitian eigenvalue equations in linear-response electronic-structure theory Julien Toulouse Université Pierre & Marie Curie and CNRS, 4 place Jussieu, Paris, France Web page: www.lct.jussieu.fr/pagesperso/toulouse/
More informationLines and points. Lines and points
omogeneous coordinates in the plane Homogeneous coordinates in the plane A line in the plane a + by + c is represented as (a, b, c). A line is a subset of points in the plane. All vectors (ka, kb, kc)
More informationEntanglement in Topological Phases
Entanglement in Topological Phases Dylan Liu August 31, 2012 Abstract In this report, the research conducted on entanglement in topological phases is detailed and summarized. This includes background developed
More informationEntanglement in Many-Body Fermion Systems
Entanglement in Many-Body Fermion Systems Michelle Storms 1, 2 1 Department of Physics, University of California Davis, CA 95616, USA 2 Department of Physics and Astronomy, Ohio Wesleyan University, Delaware,
More informationElectron-electron interactions and Dirac liquid behavior in graphene bilayers
Electron-electron interactions and Dirac liquid behavior in graphene bilayers arxiv:85.35 S. Viola Kusminskiy, D. K. Campbell, A. H. Castro Neto Boston University Workshop on Correlations and Coherence
More informationThe 4th Windsor Summer School on Condensed Matter Theory Quantum Transport and Dynamics in Nanostructures Great Park, Windsor, UK, August 6-18, 2007
The 4th Windsor Summer School on Condensed Matter Theory Quantum Transport and Dynamics in Nanostructures Great Park, Windsor, UK, August 6-18, 2007 Kondo Effect in Metals and Quantum Dots Jan von Delft
More informationSession 1 - XY model in two dimensions
Ecole de Physique des Houches August, 007 Predoctoral school in statistical physics Quantum Statistical Mechanics Session 1 - XY model in two dimensions The XY model is defined as H = J r ˆµ=ˆ ŷ S r+ˆµ
More informationChapter 29. Quantum Chaos
Chapter 29 Quantum Chaos What happens to a Hamiltonian system that for classical mechanics is chaotic when we include a nonzero h? There is no problem in principle to answering this question: given a classical
More informationLecture 4 Recap: normal metals and the clectron-phonon interaction *
Phys. 598SC Fall 2015 Prof. A. J. Leggett Lecture 4 Recap: normal metals and the clectron-phonon interaction * 1. Normal metals: Sommerfeld-Bloch picture 2. Screening 3. Fermi liquid theory 4. Electron-phonon
More informationAnswers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014
Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde
More informationν=0 Quantum Hall state in Bilayer graphene: collective modes
ν= Quantum Hall state in Bilayer graphene: collective modes Bilayer graphene: Band structure Quantum Hall effect ν= state: Phase diagram Time-dependent Hartree-Fock approximation Neutral collective excitations
More informationFerromagnetism of electron gas. Abstract
Ferromagnetism of electron gas H.Gholizade Engineering Science Department, University of Tehran, Enghelab Sq., Tehran, Iran. arxiv:090.533v4 [cond-mat.stat-mech] 6 Nov 00 D.Momeni Department of physics,faculty
More informationSupersymmetric Gauge Theories, Matrix Models and Geometric Transitions
Supersymmetric Gauge Theories, Matrix Models and Geometric Transitions Frank FERRARI Université Libre de Bruxelles and International Solvay Institutes XVth Oporto meeting on Geometry, Topology and Physics:
More informationThe University of Texas at Austin Department of Electrical and Computer Engineering. EE381V: Large Scale Learning Spring 2013.
The University of Texas at Austin Department of Electrical and Computer Engineering EE381V: Large Scale Learning Spring 2013 Assignment Two Caramanis/Sanghavi Due: Tuesday, Feb. 19, 2013. Computational
More informationThe fate of the Wigner crystal in solids part II: low dimensional materials. S. Fratini LEPES-CNRS, Grenoble. Outline
The fate of the Wigner crystal in solids part II: low dimensional materials S. Fratini LEPES-CNRS, Grenoble G. Rastelli (Università dell Aquila, Italy) P. Quémerais (LEPES-CNRS Grenoble) Outline competing
More informationThe Hartree-Fock Method
1/41 18. December, 2015 Problem Statement First Quantization Second Quantization References Overview 1. Problem Statement 2. First Quantization 3. Second Quantization 4. 5. References 2/41 Problem Statement
More informationI. QUANTUM MONTE CARLO METHODS: INTRODUCTION AND BASICS
I. QUANTUM MONTE CARLO METHODS: INTRODUCTION AND BASICS Markus Holzmann LPMMC, UJF, Grenoble, and LPTMC, UPMC, Paris markus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/markus (Dated: January 24, 2012)
More informationLecture on First-principles Computation (2): The Hartree-Fock theory and the homogeneous electron gas
Lecture on First-principles Computation (2): The Hartree-Fock theory and the homogeneous electron gas 任新国 (Xinguo Ren) 中国科学技术大学量子信息重点实验室 Key Laboratory of Quantum Information, USTC Hefei, 2016.9.18 Recall:
More informationClassical field theory 2012 (NS-364B) Feynman propagator
Classical field theory 212 (NS-364B Feynman propagator 1. Introduction States in quantum mechanics in Schrödinger picture evolve as ( Ψt = Û(t,t Ψt, Û(t,t = T exp ı t dt Ĥ(t, (1 t where Û(t,t denotes the
More informationTranslation Symmetry, Space Groups, Bloch functions, Fermi energy
Translation Symmetry, Space Groups, Bloch functions, Fermi energy Roberto Orlando and Silvia Casassa Università degli Studi di Torino July 20, 2015 School Ab initio Modelling of Solids (UniTo) Symmetry
More informationExchange energy and generalized polarization in the presence of spin-orbit coupling in two dimensions
Physics Physics Research Publications Purdue University Year 007 Exchange energy and generalized polarization in the presence of spin-orbit coupling in two dimensions S. Chesi G. F. Giuliani This paper
More information3 Symmetry Protected Topological Phase
Physics 3b Lecture 16 Caltech, 05/30/18 3 Symmetry Protected Topological Phase 3.1 Breakdown of noninteracting SPT phases with interaction Building on our previous discussion of the Majorana chain and
More informationDisordered Hyperuniformity: Liquid-like Behaviour in Structural Solids, A New Phase of Matter?
Disordered Hyperuniformity: Liquid-like Behaviour in Structural Solids, A New Phase of Matter? Kabir Ramola Martin Fisher School of Physics, Brandeis University August 19, 2016 Kabir Ramola Disordered
More informationSYLLABUS FOR ENTRANCE EXAMINATION NANYANG TECHNOLOGICAL UNIVERSITY FOR INTERNATIONAL STUDENTS A-LEVEL MATHEMATICS
SYLLABUS FOR ENTRANCE EXAMINATION NANYANG TECHNOLOGICAL UNIVERSITY FOR INTERNATIONAL STUDENTS A-LEVEL MATHEMATICS STRUCTURE OF EXAMINATION PAPER. There will be one -hour paper consisting of 4 questions..
More informationPhysics 221A Fall 2017 Notes 27 The Variational Method
Copyright c 2018 by Robert G. Littlejohn Physics 221A Fall 2017 Notes 27 The Variational Method 1. Introduction Very few realistic problems in quantum mechanics are exactly solvable, so approximation methods
More informationLandau s Fermi Liquid Theory
Thors Hans Hansson Stockholm University Outline 1 Fermi Liquids Why, What, and How? Why Fermi liquids? What is a Fermi liquids? Fermi Liquids How? 2 Landau s Phenomenological Approach The free Fermi gas
More informationFirst-Principles Calculation of Exchange Interactions
Chapter 2 First-Principles Calculation of Exchange Interactions Before introducing the first-principles methods for the calculation of exchange interactions in magnetic systems we will briefly review two
More information1 Equal-time and Time-ordered Green Functions
1 Equal-time and Time-ordered Green Functions Predictions for observables in quantum field theories are made by computing expectation values of products of field operators, which are called Green functions
More informationThe Simple Harmonic Oscillator
The Simple Harmonic Oscillator Michael Fowler, University of Virginia Einstein s Solution of the Specific Heat Puzzle The simple harmonic oscillator, a nonrelativistic particle in a potential ½C, is a
More informationRenormalization Group for the Two-Dimensional Ising Model
Chapter 8 Renormalization Group for the Two-Dimensional Ising Model The two-dimensional (2D) Ising model is arguably the most important in statistical physics. This special status is due to Lars Onsager
More informationPart E1. Transient Fields: Leapfrog Integration. Prof. Dr.-Ing. Rolf Schuhmann
Part E1 Transient Fields: Leapfrog Integration Prof. Dr.-Ing. Rolf Schuhmann MAXWELL Grid Equations in time domain d 1 h() t MC e( t) dt d 1 e() t M Ch() t j( t) dt Transient Fields system of 1 st order
More informationLecture notes for QFT I (662)
Preprint typeset in JHEP style - PAPER VERSION Lecture notes for QFT I (66) Martin Kruczenski Department of Physics, Purdue University, 55 Northwestern Avenue, W. Lafayette, IN 47907-036. E-mail: markru@purdue.edu
More informationThree Most Important Topics (MIT) Today
Three Most Important Topics (MIT) Today Electrons in periodic potential Energy gap nearly free electron Bloch Theorem Energy gap tight binding Chapter 1 1 Electrons in Periodic Potential We now know the
More informationTime-Independent Perturbation Theory
4 Phys46.nb Time-Independent Perturbation Theory.. Overview... General question Assuming that we have a Hamiltonian, H = H + λ H (.) where λ is a very small real number. The eigenstates of the Hamiltonian
More informationLecture 4: Hartree-Fock Theory
Lecture 4: Hartree-Fock Theory One determinant to rule them all, One determinant to find them, One determinant to bring them all and in the darkness bind them Second quantization rehearsal The formalism
More information