Supercell method for the computation of energies of crystals
|
|
- Abner Curtis
- 5 years ago
- Views:
Transcription
1 Supercell method for the computation of energies of crystals David Gontier CEREMADE, Université Paris-Dauphine Warwick EPSRC Symposium: Density Functional Theory and Beyond: Analysis and Computation Warwick, July 5, 2017 Joint work with Salma Lahbabi David Gontier Simulation of crystals 1 / 33
2 Goals: Define and compute numerically the energy per cell of a crystal. Define and compute numerically the energy of a local defect inside a crystal. What is a crystal? Remarks A periodic arrangement of (fixed) nuclei on a lattice. Modeled by a periodic density charge µ per, or by a periodic potential V per. An arrangement of moving electrons around these nuclei. The number of electrons is infinite. There is no Schrödinger-like equations a priori. What is the energy (per cell) of the electrons? David Gontier Simulation of crystals 2 / 33
3 The thermodynamic limit David Gontier Simulation of crystals 3 / 33
4 Fact: There exist good models for finite systems. Thermodynamic limit: See the infinite periodic crystal as the limit of finite systems. L = 1 L = 2 L = 3 L = Roadmap: Choose a model for finite systems (Schrödinger equation, Hartree(-Fock), DFT, ). For each L N, calculate the ground state energy E L for the charge density µ L. Take the limit L. Questions: Does the sequence L 3 E L (energy per cell) converge as L? If yes, can we characterise the limit? David Gontier Simulation of crystals 4 / 33
5 The thermodynamic limit for different models Model Existence of the limit Characterisation of the limit Schrödinger 1,2 yes no non-interacting Schrödinger yes yes DFT no (yes?) reduced-hf (rhf) 3 yes yes Hartree-Fock (HF) 3 maybe (yes?) Remarks The main difficulty for the Existence/Characterisation is to prove that the final solution is periodic. This is the case for the rhf/non-interacting model: the model is convex. Otherwise, symmetry breaking may happen (periodic problem, non periodic minimiser). Speed of convergence? We expect a very slow speed of convergence: the edge effects only vanish at the limit. Idea If the final solution is periodic, perform a periodic thermodynamic limit. 1 C. Fefferman, Commun. Math. Phys. 98 (1985), no X. Blanc, C. Le Bris, and P.-L. Lions, Comm. Part. Diff. Eq. 28 (2003). 3 I. Catto, C. Le Bris, and P.-L. Lions, Ann. Inst. H. Poincaré (C) 18 (2001), no. 6. David Gontier Simulation of crystals 5 / 33
6 The non-interacting (linear) case David Gontier Simulation of crystals 6 / 33
7 Notations for the crystal Lattice: R = az 3. Unit cell: Γ = [ a/2, a/2) 3. Reciprocal lattice: R = (2π/a) Z 3. Reciprocal unit cell: Γ = [ π/a, π/a) 3. Supercell: Γ L := LΓ. The Hamiltonian we wish to study (mean field) potential: V per L 2 per (Γ). Non-interacting Hamiltonian: H per = V per, acting on L 2 (R 3 ). Number of electrons per cell: N N. Questions How to define the energy per cell? How to compute it? Idea: Perform a periodic thermodynamic limit. David Gontier Simulation of crystals 7 / 33
8 The supercell model Full model Supercell model L Laplacian with periodic boundary conditions: L. Supercell Hamiltonian H L := 1 2 L + V per acting on L 2 per(γ L ). Number of electrons = NL 3. Roadmap Compute the total energy E L of the supercell. Define the energy per cell as the limit of L 3 E L when L. David Gontier Simulation of crystals 8 / 33
9 1-body Hamiltonian h 1 := 1 2 L + V per self-adjoint operator acting on H 1 = L 2 (Γ L ). Fact: h 1 is compact resolvent (= discrete spectrum). ε 1... ε N the smallest eigenvalues. u 1,..., u N H 1 the corresponding normalised eigenvectors. N-body non-interacting Hamiltonian N H N := h 1 (x i ) i=1 Ground state energy of H N acting on the fermionic space N H1. E 0 N = ε ε N = u 1, h 1 u u N, h 1 u N = Tr H1 [h 1 ( u 1 u u N u N )] (cyclicity of the trace) = inf { Tr H1 [h 1 γ], γ is a projector of rank N } (min-max principle) = inf { Tr H1 [h 1 γ], γ P N } (min of linear function on convex set), where P N is the convex hull of the set of projectors of rank N: P N := { γ S(H 1 ), 0 γ 1, Tr H1 (γ) = N }. The minimiser is γ N := 1 (H N ε N ), and satisfies Tr H1 (γ N ) = N. David Gontier Simulation of crystals 9 / 33
10 Supercell minimiser γ L := 1 ( ) H L < ε L F Supercell energy per unit cell Exact minimiser In our case where ε L F is chosen so that 1 [ L 3 Tr L 2 per (Γ γ L] = N. L) I L := 1 L 3 EL = 1 L 3 Tr L 2 per (Γ L) [ H L γ L]. In the thermodynamic limit, we expect γ per of the form γ per := 1 (H per < ε F ) where ε F is chosen so that Tr (γ) = N. Exact energy per cell I per = Tr [H per γ per ]. Trace per cell 1 Tr := lim L L 3 Tr L 2 per (Γ L). Questions Does the sequence L 3 E L converge to I per as L? What is the speed of convergence? Lemma In the insulating case, there exists C R + and α > 0 such that, for all L N, I L I per Ce αl. David Gontier Simulation of crystals 10 / 33
11 The Bloch-Floquet transform David Gontier Simulation of crystals 11 / 33
12 L 2 per([0, L]) = Vect { e 2iπ0x, e 2iπ( L 1 )x, e 2iπ( ( ) L 2 )x,..., e 2iπ L 1 x L, e 2iπ( ( ) ( ) ( ) L L )x, e 2iπ L+1 x 2iπ L+2 x 2iπ 2L 1 x L, e L,..., e L, } = L 2 0 L 2 2π L 2 4π L 2 2π(L 1). L L L Fibers L 2 q = Vect { e ik x e iq x, k R } = { ψ L 2 loc (R3 ), R R, ψ( + R) = e iq R ψ( ) }. L 2 q does not depend on L, q Γ. 3D case L 2 per(γ L ) = Q Λ L L 2 Q. Λ L Γ David Gontier Simulation of crystals 12 / 33
13 Supercell set of translations : R L := R Γ L. (Supercell) Bloch-Floquet transform for functions f L 2 (Γ L ) = Q Λ L f Q with f Q L 2 Q. (Supercell) Bloch fibers f Q (x) := e iq R f(x R) = e iq R τ R f(x). R R L R R L Bloch-Floquet for operators Lemma Let A L : L 2 per(γ L ) L 2 per(γ L ) be such that τ R A L = A L τ R for all R R. Then, the operator A L is block-diagonal with respect to the L 2 Q decomposition. Proof [ A L f ]Q (x) = ( ) e iq R τ R A L f (x) = e iq R A L (τ R f)(x) R R L R R L = A L e iq R τ R f (x) = A L ( ) f Q (x). R R L David Gontier Simulation of crystals 13 / 33
14 In other words, writing ÃL Q : L2 Q L2 Q, we have A L Q A L 0 A L Q = (Bloch transform = block-decomposition) Ã L Q L 3 Covariance ) L 2 q = S q (L 2 per(γ) with S q [f](x) = e iq x f(x) and (S q ) 1 = S q. Twisted Bloch transform A L Q := S QÃL Q S Q, : L 2 per(γ) L 2 per(γ). Basic properties Spectrum Trace per cell σ(a L ) = Q Λ L σ Tr L (A L ) := 1 L d Tr L 2 per (Γ L) (AL ) = 1 L d (ÃL Q ) = Q Λ L σ ( ) A L Q. Tr L 2 per (Γ) (ÃL Q ) = 1 L d Q Λ L Q Λ L Tr L 2 per (Γ) (AL Q ). David Gontier Simulation of crystals 14 / 33
15 Example: The periodic Hamiltonian Bloch transform H q := H L q = 1 2 H L = 1 2 L + V per acting on L 2 per(γ L ). i 1 + q 2 + V per = q (i 1 ) + q2 2 + V per acting on L 2 per (Γ). Supercell thermodynamic limit (L ) Γ Γ Λ = Γ Λ L=4 Λ L=8 Trace per unit volume H L = Q Λ L H L Q Tr L (A) = 1 L 3 L ˆ H per = H qdq. Γ Q Λ L A Q L Tr (A) := Γ A q dq. Supercell model Regular sampling of Γ. David Gontier Simulation of crystals 15 / 33
16 Remark: The operator H q (acting on L 2 per (Γ)) is compact resolvent: H q = ε nq u nq u nq, ε 1,q ε 2,q, u nq, u mq L 2 per (Γ) = δnm. n=1 ε nq (ev) ε F N Band diagram of the silicon q Γ Fermi energy: ε F R s.t. n=1 Γ 1(ε nq ε F ) dq = N. Insulating system: The system is insulating if there exists a gap g > 0 such that q Γ, ε N,q + g 2 ε F ε N+1,q g 2. David Gontier Simulation of crystals 16 / 33
17 Full Hamiltonian One-body density ˆ γ per = 1(H per ε F ) = 1(H q ε F ) dq. Γ }{{} γ q Energy per unit cell (update) I per = Tr (H perγ per) = Supercell model Supercell one-body density Supercell energy per unit cell (update) I L = Tr (H L γ L) = 1 L 3 Tr L 2 Γ per (Γ) (Hqγq) dq = Γ γ L = 1(H L ε F ) = Q Λ L Tr L 2 per (Γ) Q Λ L γ Q. ( HQ γ Q ) = 1 L 3 ( N ) ε nq dq n=1 Q Λ L ( N ) ε nq. n=1 Error for the energy ( I L I per = 1 N ) L 3 ε nq Q Λ L n=1 ( N ) ε nq dq. Γ n=1 This is the difference between a Riemann sum and an integral. David Gontier Simulation of crystals 17 / 33
18 Tools in complex analysis David Gontier Simulation of crystals 18 / 33
19 Integrand Properties K(q) := N ε nq. n=1 The function K is R -periodic. The functions q ε nq are not smooth in general (conic singularities). However, the sum of the eigenvalues is smooth. Band diagram of graphene. Lemma ( 4,5,6 ) If the system is an insulator, then the (R -periodic) function K admits an analytic extension on some complex strip S A := R 3 + i[ A, A] 3 with A > 0. Link to the fact that γ per (x, y) Ce α x y. 4 W. Kohn, Phys. Rev. 115 (1959). 5 J. Des Cloizeaux, Phys. Rev. 135 (1964). 6 C. Brouder, G. Panati, M. Calandra, C. Mourougane, and N. Marzari, Phys. Rev. Lett. 98 (2007). David Gontier Simulation of crystals 19 / 33
20 Elements of the proof Cauchy residual formula γ q = 1(H q ε F ) = 1 2iπ Analytic continuation C dλ λ H q. C ε F σ(h) H q = q (i 1 ) + q2 2 + Vper = Hz = z (i 1 ) + zt z 2 + Vper Lemma There exists A > 0 such that z S A, γ z := 1 2iπ C is well-defined. Moreover, the map z γ z is analytic on S A. Integrand K(q) := dλ λ H z N ε nq = Tr L 2 per (Γ) (H qγ q ) = K(z) := Tr L 2 per (Γ) (H zγ z ). n=1 K(z) is R -periodic and analytic on S A. David Gontier Simulation of crystals 20 / 33
21 Lemma (classical) If f : S A C is analytic on S A for some A > 0 and satisfies f(z + k) = f(z) for all k R, then there exists C R + and α > 0 such that L N, 1 L 3 f(q) f(q)dq Q Λ Γ L Ce αl. Conclusion Lemma (DG, Salma Lahbabi) In the insulating case, there exists C R + and α > 0 such that, for all L N, L 3 E L I per Ce αl. (Convergence of the energy per cell) ρ γper ρ γ L Ce αl. (Convergence of the electronic density) L per David Gontier Simulation of crystals 21 / 33
22 Numerical illustration (Linear model for silicium taken in M.L. Cohen and T.K. Bergstresser, Phys. Rev. 141 (1966). Speed of convergence for the silicium (linear model). David Gontier Simulation of crystals 22 / 33
23 The non-linear case (reduced Hartree-Fock case) David Gontier Simulation of crystals 23 / 33
24 How to define an interacting model? Take an effective model (DFT, HF, ), and perform a (periodic) thermodynamic limit. Problem: symmetry breaking Lemma (current work with Mathieu Lewin and Faizan Nazar) Let { I α (L) := 1 ˆ } L 3 inf Tr L 2 (Γ L ) ˆΓ ( γ) + V per ρ γ + D L (ρ γ, ρ γ ) α ρ 4/3 γ, γ P(Γ L ), L Γ L then, for α > α c, it holds that I α (2) < I α (1). Idea Enforce the 1-periodicity (for instance using convexity) = reduced Hartree-Fock (rhf). Non-interacting energy as a minimisation problem Set of one-body density matrices Recall the non-interacting case I per = Tr (H per γ per ) = inf { Tr (H per γ), γ P per, Tr (γ) = N }. P per = { γ S(L 2 (R 3 )), 0 γ 1, R R, τ R γ = γτ R }. Set of supercell one-body density matrices } P per = {γ S(L 2 per(γ L )), 0 γ 1, R R, τ R γ = γτ R. David Gontier Simulation of crystals 24 / 33
25 Charge density: µ per reduced Hartree-Fock I rhf per { } = inf Iper rhf (γ), γ P per, Tr (γ) = N. with (µ per is the charge density) Periodic Green s function G 1 = 4π R R I rhf per (γ) := 1 2 Tr ( γ) D 1(ρ γ µ per, ρ γ µ per ). ( δr Γ 1) so that G 1 (x) = 4π Γ k R \{0} Periodic Coulomb potential f, g L 2 per(γ), D 1 (f, g) := f(x)g 1 (x y)g(y) dxdy. (Γ) 2 e ik x k 2. Lemma ( Catto, Le Bris, Lions and Cancès, Deleurence, Lewin ) The problem Iper rhf admits a unique minimizer γ per. This minimizer satisfies the Euler-Lagrange equations γ per = 1(H per ε F ) H per = Vper acting on L2 (R 3 ) V per = (ρ γper µ per) Γ G 1. David Gontier Simulation of crystals 25 / 33
26 Supercell energy (rhf) ( EL {E rhf = inf rhf (γ L ), γ L P L, Tr L 2 per (Γ L ) γ L) = NL 3}, with E rhf (γ L ) := 1 2 Tr L 2 per (Γ L) ( L γ L) D L(ρ γ L µ per, ρ γ L µ per). Supercell Green s function G L = 4π R LR ( δr Γ L 1) so that G L (x) = 4π Γ L k L 1 R \{0} Supercell Coulomb potential f, g L 2 per(γ L ), D L (f, g) := f(x)g L (x y)g(y) dxdy. (Γ L ) 2 e ik x k 2 = 1 ( x ) L G 1. L Lemma The problem EL rhf admits a unique minimizer γl. This minimizer satisfies the Euler-Lagrange equations γ L = 1(H L ε L F ) H L = 1 2 L + V L acting on L 2 per(γ L ) V L = (ρ γ L µ per ) ΓL G L. David Gontier Simulation of crystals 26 / 33
27 Question: What is the speed of convergence of L 3 EL rhf IrHF per? Remark. Non-linear setting: The potential V depends on L: H per = V per and H L = 1 2 L + V L. Idea: Construct good test functions for the minimization problems. Problem: γ per P per and γ L P L cannot be compared. By convexity of the problems w.r.t. the density ρ, we deduce that ρ γper and ρ γ L are R-periodic. V per and V L are R-periodic. = Construct test functions from ρ γper and ρ γ L, and use the results for the linear case. David Gontier Simulation of crystals 27 / 33
28 Lemma (DG, Salma Lahbabi) There exists C R + and α > 0 such that, for all L L gap, L 3 EL rhf IrHF per Ce αl. (Convergence of the energy per unit cell) ρ γper ρ γ L Ce αl. (Convergence of the electronic density) L per Rates of convergence for the reduced Hartree-Fock model. David Gontier Simulation of crystals 28 / 33
29 Crystals with local defects David Gontier Simulation of crystals 29 / 33
30 General picture Perfect crystal H per = Vper acting on L2 per(γ). Local defect (charge) ν L 2 (R 3 ) with compact support in Γ. Crystal with local defect H ν := 1 ˆ 2 + V ν(y) per + R 3 y dy. Local defects How to define the energy of the defect? The energy per cell is no longer a meaningful quantity. The energy of the crystal (with and without defect) is infinite. Idea: Perform a supercell thermodynamic limit. David Gontier Simulation of crystals 30 / 33
31 Supercell thermodynamic limit Supercell local defect ν ν L L 2 per(γ L ) such that ν L = ν on Γ L. Local defect Supercell defect The defect interacts with its periodic images! Roadmap For all L N, the supercell energy of the crystal with and without defect is finite. Define the supercell energy of the defect as Jν L = Eν L E L. Does Jν L have limit as L? Yes in rhf 7 : J ν. Can we characterize the limit? Yes, as a minimization problem. What is the speed of convergence of Jν L J ν? 7 E. Cancès, A. Deleurence, and M. Lewin, Commun. Math. Phys. 281 (2008) David Gontier Simulation of crystals 31 / 33
32 Lemma (DG, Salma Lahbabi) There exists η > 0, L N, C R +, and α > 0 such that, for all ν L 2 per (Γ) with compact support in Γ and satisfying ν L 2 < η, and for all L L, it holds that J ν Jν L 1 mq 2 ( L 2ϵ C ν 3 L 2 + ν 2 L 2 L 3 + ν L 2 e αl where q = R 3 ν is the net charge of the defect, m is the Madelung constant of the crystal, defined by m := lim x 0 G 1 (x) 1 x, and ϵ is the macroscopic dielectric constant of the crystal. Comments Recover the term predicted by Leslie and Gillan 8, and by Makov and Payne 9. Slow convergence in the defect case (interaction of the defect with its images). The term 1 mq 2 can be computed with low computational time. L 2ϵ = We can improve the numerical convergence by substracting by hand this term. ), 8 M. Leslie and M.J. Gillan. J. Phys. C 18 (1985). 9 G. Makov and M.C. Payne. Phys. Rev. B, 51 (1995). David Gontier Simulation of crystals 32 / 33
33 Conclusion Convergence of supercell models Convergence of Riemann sums. Exponential rate of convergence for perfect crystals. Slow rate of convergence for crystal with local defects. Identification of the L 1 term allows to speed the rate of convergence. Work in progress (see Antoine Levitt s talk) Speed of convergence in the metallic case. (with E. Cancès, V. Ehrlacher, A. Levitt and D. Lombardi) References DG, S. Lahbabi, Convergence rates of supercell calculations in the reduced Hartree-Fock model (M2AN, 50,5 (2016), arxiv ). DG, S. Lahbabi, Supercell calculations in the reduced Hartree-Fock model for crystals with local defects (Appl. Math. Res. Express (2016), arxiv ). Thank you for your attention. David Gontier Simulation of crystals 33 / 33
Exact embedding of local defects in crystals
Exact embedding of local defects in crystals Mathieu LEWIN Mathieu.Lewin@math.cnrs.fr (CNRS & University of Cergy-Pontoise) joint work with É. Cancès & A. Deleurence (ENPC, Marne-La-Vallée) Minneapolis,
More informationDerivation of Pekar s polaron from a microscopic model of quantum crystal
Derivation of Pekar s polaron from a microscopic model of quantum crystal Mathieu LEWIN Mathieu.Lewin@math.cnrs.fr (CNRS & University of Cergy-Pontoise) joint work with Nicolas Rougerie (Grenoble, France)
More informationEvolution en temps des défauts dans les cristaux
Evolution en temps des défauts dans les cristaux Gabriel STOLTZ stoltz@cermics.enpc.fr (CERMICS, Ecole des Ponts & MICMAC team, INRIA Rocquencourt) Travail en collaboration avec E. Cancès Séminaire de
More informationA perturbation-based reduced basis method for parametric eigenvalue problems
A perturbation-based reduced basis method for parametric eigenvalue problems Eric Cancès 1, Virginie Ehrlacher 1, David Gontier 2, Damiano Lombardi 3 and Antoine Levitt 1 1. CERMICS, Ecole des Ponts ParisTech
More informationNumerical construction of Wannier functions
July 12, 2017 Internship Tutor: A. Levitt School Tutor: É. Cancès 1/27 Introduction Context: Describe electrical properties of crystals (insulator, conductor, semi-conductor). Applications in electronics
More informationThe microscopic origin of the macroscopic dielectric permittivity of crystals
The microscopic origin of the macroscopic dielectric permittivity of crystals Gabriel STOLTZ stoltz@cermics.enpc.fr (CERMICS, Ecole des Ponts & MICMAC team, INRIA Rocquencourt) Work in collaboration with
More informationarxiv: v1 [math-ph] 30 Aug 2011
arxiv:8.593v [math-ph] 3 Aug 2 Derivation of Pekar s Polarons from a Microscopic Model of Quantum Crystals Mathieu LEWIN CNRS & Department of Mathematics (UMR 888), University of Cergy-Pontoise, 95 Cergy-Pontoise,
More informationDynamics of electrons in crystalline solids: Wannier functions,
Title: Dynamics of electrons in crystalline solids: Wannier functions, Berry curvature and related issues Name: Affil./Addr.: Gianluca Panati La Sapienza University of Rome Rome, 00185 Italy E-mail: panati@mat.uniroma1.it
More informationThe Localization Dichotomy for Periodic Schro dinger Operators
The Localization Dichotomy for Periodic Schro dinger Operators Gianluca Panati Dipartimento di Matematica G. Castelnuovo based on joint work with D. Monaco, A. Pisante and S. Teufel Tosio Kato Centennial
More informationarxiv: v1 [math-ph] 8 Oct 2018
SCREENING IN THE FINITE-TEMPERATURE REDUCED HARTREE-FOCK MODEL ANTOINE LEVITT Abstract. We prove the existence of solutions of the reduced Hartree-Fock equations at finite temperature for a periodic crystal
More informationarxiv: v1 [math-ph] 18 Oct 2010
The Microscopic Origin of the Macroscopic Dielectric Permittivity of Crystals: A Mathematical Viewpoint. Éric Cancès, Mathieu Lewin and Gabriel Stoltz arxiv:11.3494v1 [math-ph] 18 Oct 21 Abstract The purpose
More informationMagnetic symmetries:
Magnetic symmetries: applications and prospects Giuseppe De Nittis (AGM, Université de Cergy-Pontoise) Friedrich Schiller Universität, Mathematisches Institut, Jena. April 18, 2012 Joint work with: M.
More informationElectronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch
Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation There are 10 20 atoms in 1 mm 3
More informationMaximally localized Wannier functions for periodic Schrödinger operators: the MV functional and the geometry of the Bloch bundle
Maximally localized Wannier functions for periodic Schrödinger operators: the MV functional and the geometry of the Bloch bundle Gianluca Panati Università di Roma La Sapienza 4ème Rencontre du GDR Dynamique
More informationElectronic Structure of Crystalline Solids
Electronic Structure of Crystalline Solids Computing the electronic structure of electrons in solid materials (insulators, conductors, semiconductors, superconductors) is in general a very difficult problem
More informationSuperfluidity & Bogoliubov Theory: Rigorous Results
Superfluidity & Bogoliubov Theory: Rigorous Results Mathieu LEWIN mathieu.lewin@math.cnrs.fr (CNRS & Université Paris-Dauphine) collaborations with P.T. Nam (Vienna), N. Rougerie (Grenoble), B. Schlein
More informationThermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model
Thermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model Nikolaj Veniaminov (in collaboration with Frédéric Klopp) CEREMADE, University of Paris IX Dauphine
More informationApplications of Matrix Functions Part III: Quantum Chemistry
Applications of Matrix Functions Part III: Quantum Chemistry Emory University Department of Mathematics and Computer Science Atlanta, GA 30322, USA Prologue The main purpose of this lecture is to present
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 informationThe semiclassical. model for adiabatic slow-fast systems and the Hofstadter butterfly
The semiclassical. model for adiabatic slow-fast systems and the Hofstadter butterfly Stefan Teufel Mathematisches Institut, Universität Tübingen Archimedes Center Crete, 29.05.2012 1. Reminder: Semiclassics
More informationLocal Density Approximation for the Almost-bosonic Anyon Gas. Michele Correggi
Local Density Approximation for the Almost-bosonic Anyon Gas Michele Correggi Università degli Studi Roma Tre www.cond-math.it QMATH13 Many-body Systems and Statistical Mechanics joint work with D. Lundholm
More informationSELF-ADJOINTNESS OF DIRAC OPERATORS VIA HARDY-DIRAC INEQUALITIES
SELF-ADJOINTNESS OF DIRAC OPERATORS VIA HARDY-DIRAC INEQUALITIES MARIA J. ESTEBAN 1 AND MICHAEL LOSS Abstract. Distinguished selfadjoint extension of Dirac operators are constructed for a class of potentials
More informationGround-State Energies of Coulomb Systems and Reduced One-Particle Density Matrices
Ground-State Energies of Coulomb Systems and Reduced One-Particle Density Matrices Heinz Siedentop Chennai, August 16, 2010 I. Introduction Atomic Schrödinger operator H N,Z := N n=1 self-adjointly realized
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 informationEstimates for the resolvent and spectral gaps for non self-adjoint operators. University Bordeaux
for the resolvent and spectral gaps for non self-adjoint operators 1 / 29 Estimates for the resolvent and spectral gaps for non self-adjoint operators Vesselin Petkov University Bordeaux Mathematics Days
More informationMolecular propagation through crossings and avoided crossings of electron energy levels
Molecular propagation through crossings and avoided crossings of electron energy levels George A. Hagedorn Department of Mathematics and Center for Statistical Mechanics and Mathematical Physics Virginia
More informationPost-processing of the planewave approximation of Schrödinger equations. Part I: linear operators
Post-processing of the planewave approximation of Schrödinger equations. Part I: linear operators Eric Cancès, Geneviève Dusson, Yvon Maday, Benjamin Stamm, Martin Vohralík To cite this version: Eric Cancès,
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationThe spectral zeta function
The spectral zeta function Bernd Ammann June 4, 215 Abstract In this talk we introduce spectral zeta functions. The spectral zeta function of the Laplace-Beltrami operator was already introduced by Minakshisundaram
More informationFermions in the unitary regime at finite temperatures from path integral auxiliary field Monte Carlo simulations
Fermions in the unitary regime at finite temperatures from path integral auxiliary field Monte Carlo simulations Aurel Bulgac,, Joaquin E. Drut and Piotr Magierski University of Washington, Seattle, WA
More informationWeakly non-linear completely resonant hamiltonian PDEs and the problem of weak turbulence
Sergei Kuksin Weakly non-linear completely resonant hamiltonian PDEs and the problem of weak turbulence (Toronto, 10 January, 2014 ) 1 1 Introduction: weak turbulence (WT) (One of) origines: Rudolf Peierls,
More informationElectronic Structure Methodology 1
Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E
More informationElectronic-structure calculations at macroscopic scales
Electronic-structure calculations at macroscopic scales M. Ortiz California Institute of Technology In collaboration with: K. Bhattacharya, V. Gavini (Caltech), J. Knap (LLNL) BAMC, Bristol, March, 2007
More informationThe Borchers-Uhlmann Algebra
The Borchers-Uhlmann Algebra and its Descendants Jakob Yngvason Göttingen, July 29, 2009 Topics: 1. The BU algebra 2. Wightman functionals 3. Abundance of states satisfying parts of the W-axioms 4. Bounded
More informationSpectral Pollution and How to Avoid It (With Applications to Dirac and Periodic Schrödinger Operators)
Spectral Pollution and How to Avoid It With Applications to Dirac and Periodic Schrödinger Operators Mathieu Lewin, Eric Séré To cite this version: Mathieu Lewin, Eric Séré. Spectral Pollution and How
More informationNodal Sets of High-Energy Arithmetic Random Waves
Nodal Sets of High-Energy Arithmetic Random Waves Maurizia Rossi UR Mathématiques, Université du Luxembourg Probabilistic Methods in Spectral Geometry and PDE Montréal August 22-26, 2016 M. Rossi (Université
More informationWavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts.
Wavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts. Silvia Casassa Università degli Studi di Torino July 12, 2017 Minnesota Workshop on ab initio MSC Symmetry
More informationRecursion Systems and Recursion Operators for the Soliton Equations Related to Rational Linear Problem with Reductions
GMV The s Systems and for the Soliton Equations Related to Rational Linear Problem with Reductions Department of Mathematics & Applied Mathematics University of Cape Town XIV th International Conference
More informationSymmetry breaking in Caffarelli-Kohn-Nirenberg inequalities
Symmetry breaking in Caffarelli-Kohn-Nirenberg inequalities p. 1/47 Symmetry breaking in Caffarelli-Kohn-Nirenberg inequalities Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine
More informationTopological Physics in Band Insulators IV
Topological Physics in Band Insulators IV Gene Mele University of Pennsylvania Wannier representation and band projectors Modern view: Gapped electronic states are equivalent Kohn (1964): insulator is
More informationClusters and Percolation
Chapter 6 Clusters and Percolation c 2012 by W. Klein, Harvey Gould, and Jan Tobochnik 5 November 2012 6.1 Introduction In this chapter we continue our investigation of nucleation near the spinodal. We
More informationL 2 extension theorem for sections defined on non reduced analytic subvarieties
L 2 extension theorem for sections defined on non reduced analytic subvarieties Jean-Pierre Demailly Institut Fourier, Université de Grenoble Alpes & Académie des Sciences de Paris Conference on Geometry
More informationElectron levels in a periodic potential: general remarks. Daniele Toffoli January 11, / 46
Electron levels in a periodic potential: general remarks Daniele Toffoli January 11, 2017 1 / 46 Outline 1 Mathematical tools 2 The periodic potential 3 Bloch s theorem and Born-von Karman boundary conditions
More informationSpectral stability of periodic waves in dispersive models
in dispersive models Collaborators: Th. Gallay, E. Lombardi T. Kapitula, A. Scheel in dispersive models One-dimensional nonlinear waves Standing and travelling waves u(x ct) with c = 0 or c 0 periodic
More informationElectrons in Crystals. Chris J. Pickard
Electrons in Crystals Chris J. Pickard Electrons in Crystals The electrons in a crystal experience a potential with the periodicity of the Bravais lattice: U(r + R) = U(r) The scale of the periodicity
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 informationGROUND-STATES FOR THE LIQUID DROP AND TFDW MODELS WITH LONG-RANGE ATTRACTION
GROUND-STATES FOR THE LIQUID DROP AND TFDW MODELS WITH LONG-RANGE ATTRACTION STAN ALAMA, LIA BRONSARD, RUSTUM CHOKSI, AND IHSAN TOPALOGLU Abstract. We prove that both the liquid drop model in with an attractive
More informationPseudopotentials for hybrid density functionals and SCAN
Pseudopotentials for hybrid density functionals and SCAN Jing Yang, Liang Z. Tan, Julian Gebhardt, and Andrew M. Rappe Department of Chemistry University of Pennsylvania Why do we need pseudopotentials?
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 informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More informationCONTINUITY WITH RESPECT TO DISORDER OF THE INTEGRATED DENSITY OF STATES
Illinois Journal of Mathematics Volume 49, Number 3, Fall 2005, Pages 893 904 S 0019-2082 CONTINUITY WITH RESPECT TO DISORDER OF THE INTEGRATED DENSITY OF STATES PETER D. HISLOP, FRÉDÉRIC KLOPP, AND JEFFREY
More informationEffective Potentials Generated by Field Interaction in the Quasi-Classical Limit. Michele Correggi
Effective Potentials Generated by Field Interaction in the Quasi-Classical Limit Michele Correggi Dipartimento di Matematica Quantum Mean Field and Related Problems LAGA Université Paris 13 joint work
More information3: Density Functional Theory
The Nuts and Bolts of First-Principles Simulation 3: Density Functional Theory CASTEP Developers Group with support from the ESF ψ k Network Density functional theory Mike Gillan, University College London
More informationThe general reason for approach to thermal equilibrium of macroscopic quantu
The general reason for approach to thermal equilibrium of macroscopic quantum systems 10 May 2011 Joint work with S. Goldstein, J. L. Lebowitz, C. Mastrodonato, and N. Zanghì. Claim macroscopic quantum
More informationWaves in Honeycomb Structures
Waves in Honeycomb Structures Michael I. Weinstein Columbia University Nonlinear Schrödinger Equations: Theory and Applications Heraklion, Crete / Greece May 20-24, 2013 Joint work with C.L. Fefferman
More informationarxiv: v2 [math-ph] 13 Feb 2013
FACTORIZATIONS OF RATIONAL MATRIX FUNCTIONS WITH APPLICATION TO DISCRETE ISOMONODROMIC TRANSFORMATIONS AND DIFFERENCE PAINLEVÉ EQUATIONS ANTON DZHAMAY SCHOOL OF MATHEMATICAL SCIENCES UNIVERSITY OF NORTHERN
More informationMany-Body Fermion Density Matrix: Operator-Based Truncation Scheme
Many-Body Fermion Density Matrix: Operator-Based Truncation Scheme SIEW-ANN CHEONG and C. L. HENLEY, LASSP, Cornell U March 25, 2004 Support: NSF grants DMR-9981744, DMR-0079992 The Big Picture GOAL Ground
More informationThe Klein-Gordon Equation Meets the Cauchy Horizon
Enrico Fermi Institute and Department of Physics University of Chicago University of Mississippi May 10, 2005 Relativistic Wave Equations At the present time, our best theory for describing nature is Quantum
More informationNon-radial solutions to a bi-harmonic equation with negative exponent
Non-radial solutions to a bi-harmonic equation with negative exponent Ali Hyder Department of Mathematics, University of British Columbia, Vancouver BC V6TZ2, Canada ali.hyder@math.ubc.ca Juncheng Wei
More informationSSH Model. Alessandro David. November 3, 2016
SSH Model Alessandro David November 3, 2016 Adapted from Lecture Notes at: https://arxiv.org/abs/1509.02295 and from article: Nature Physics 9, 795 (2013) Motivations SSH = Su-Schrieffer-Heeger Polyacetylene
More informationFundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering
Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering Outline PART 1: Fundamentals of Density functional theory (DFT)
More informationDesingularization of an immersed self-shrinker.
Desingularization of an immersed self-shrinker. Xuan Hien Nguyen (joint work with S. Kleene and G. Drugan) Iowa State University May 29, 2014 1 / 28 Examples of gluing constructions Summary of a classical
More informationSTABLE STEADY STATES AND SELF-SIMILAR BLOW UP SOLUTIONS
STABLE STEADY STATES AND SELF-SIMILAR BLOW UP SOLUTIONS FOR THE RELATIVISTIC GRAVITATIONAL VLASOV- POISSON SYSTEM Mohammed Lemou CNRS and IRMAR, Rennes Florian Méhats University of Rennes 1 and IRMAR Pierre
More informationEnsembles and incomplete information
p. 1/32 Ensembles and incomplete information So far in this course, we have described quantum systems by states that are normalized vectors in a complex Hilbert space. This works so long as (a) the system
More informationPresenter: Noriyoshi Fukaya
Y. Martel, F. Merle, and T.-P. Tsai, Stability and Asymptotic Stability in the Energy Space of the Sum of N Solitons for Subcritical gkdv Equations, Comm. Math. Phys. 31 (00), 347-373. Presenter: Noriyoshi
More informationSurfaces in spheres: biharmonic and CMC
Surfaces in spheres: biharmonic and CMC E. Loubeau (joint with C. Oniciuc) Université de Bretagne Occidentale, France May 2014 Biharmonic maps Definition Bienergy functional at φ : (M, g) (N, h) E 2 (φ)
More informationGround State and Charge Renormalization in a Nonlinear Model of Relativistic Atoms
Ground State and Charge Renormalization in a Nonlinear Model of Relativistic Atoms Philippe GRAVEJAT a, Mathieu LEWIN b and Éric SÉRÉa a CEREMADE, UMR 7534, Université Paris-Dauphine, Place du Maréchal
More informationTHE ELECTRONIC STRUCTURE OF SMOOTHLY DEFORMED CRYSTALS: WANNIER FUNCTIONS AND THE CAUCHY-BORN RULE
THE ELECTRONIC STRUCTURE OF SMOOTHLY DEFORMED CRYSTALS: WANNIER FUNCTIONS AND THE CAUCHY-BORN RULE WEINAN E AND JIANFENG LU Abstract. The electronic structure of a smoothly deformed crystal is analyzed
More informationDFT: Exchange-Correlation
DFT: Local functionals, exact exchange and other post-dft methods Stewart Clark University of Outline Introduction What is exchange and correlation? Quick tour of XC functionals (Semi-)local: LDA, PBE,
More informationBRST and Dirac Cohomology
BRST and Dirac Cohomology Peter Woit Columbia University Dartmouth Math Dept., October 23, 2008 Peter Woit (Columbia University) BRST and Dirac Cohomology October 2008 1 / 23 Outline 1 Introduction 2 Representation
More informationModelling in photonic crystal structures
Modelling in photonic crystal structures Kersten Schmidt MATHEON Nachwuchsgruppe Multiscale Modelling and Scientific Computing with PDEs in collaboration with Dirk Klindworth (MATHEON, TU Berlin) Holger
More informationIntroduction to the Mathematics of the XY -Spin Chain
Introduction to the Mathematics of the XY -Spin Chain Günter Stolz June 9, 2014 Abstract In the following we present an introduction to the mathematical theory of the XY spin chain. The importance of this
More informationRepresentation of the polar cone of convex functions and applications
Representation of the polar cone of convex functions and applications G. Carlier, T. Lachand-Robert October 23, 2006 version 2.1 Abstract Using a result of Y. Brenier [1], we give a representation of the
More informationSloppy Nuclear Energy Density Functionals: effective model optimisation. T. Nikšić and D. Vretenar
Sloppy Nuclear Energy Density Functionals: effective model optimisation T. Nikšić and D. Vretenar Energy Density Functionals the nuclear many-body problem is effectively mapped onto a one-body problem
More informationON POISSON BRACKETS COMPATIBLE WITH ALGEBRAIC GEOMETRY AND KORTEWEG DE VRIES DYNAMICS ON THE SET OF FINITE-ZONE POTENTIALS
ON POISSON BRACKETS COMPATIBLE WITH ALGEBRAIC GEOMETRY AND KORTEWEG DE VRIES DYNAMICS ON THE SET OF FINITE-ZONE POTENTIALS A. P. VESELOV AND S. P. NOVIKOV I. Some information regarding finite-zone potentials.
More informationMicrolocal analysis and inverse problems Lecture 3 : Carleman estimates
Microlocal analysis and inverse problems ecture 3 : Carleman estimates David Dos Santos Ferreira AGA Université de Paris 13 Monday May 16 Instituto de Ciencias Matemáticas, Madrid David Dos Santos Ferreira
More informationQuantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University
PY502, Computational Physics, December 12, 2017 Quantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University Advancing Research in Basic Science and Mathematics Example:
More informationWelcome to the Solid State
Max Planck Institut für Mathematik Bonn 19 October 2015 The What 1700s 1900s Since 2005 Electrical forms of matter: conductors & insulators superconductors (& semimetals & semiconductors) topological insulators...
More informationA new proof of the analyticity of the electronic density of molecules.
A new proof of the analyticity of the electronic density of molecules. Thierry Jecko AGM, UMR 8088 du CNRS, Université de Cergy-Pontoise, Département de mathématiques, site de Saint Martin, 2 avenue Adolphe
More informationLecture 4: Band theory
Lecture 4: Band theory Very short introduction to modern computational solid state chemistry Band theory of solids Molecules vs. solids Band structures Analysis of chemical bonding in Reciprocal space
More informationTHE DIRAC EQUATION (A REVIEW) We will try to find the relativistic wave equation for a particle.
THE DIRAC EQUATION (A REVIEW) We will try to find the relativistic wave equation for a particle. First, we introduce four dimensional notation for a vector by writing x µ = (x, x 1, x 2, x 3 ) = (ct, x,
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 informationIntroduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić
Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824
More informationPekar s Ansatz and the Ground-State Symmetry of a Bound Polaron
Pekar s Ansatz and the Ground-State Symmetry of a Bound Polaron Rohan Ghanta Georgia Institute of Technology March 9, 018 The Polaron Model In 1937 H. Fröhlich proposed a model known today as the Fröhlich
More information1 Sectorial operators
1 1 Sectorial operators Definition 1.1 Let X and A : D(A) X X be a Banach space and a linear closed operator, respectively. If the relationships i) ρ(a) Σ φ = {λ C : arg λ < φ}, where φ (π/2, π); ii) R(λ,
More informationWave function methods for the electronic Schrödinger equation
Wave function methods for the electronic Schrödinger equation Zürich 2008 DFG Reseach Center Matheon: Mathematics in Key Technologies A7: Numerical Discretization Methods in Quantum Chemistry DFG Priority
More informationBosonic quadratic Hamiltonians. diagonalization
and their diagonalization P. T. Nam 1 M. Napiórkowski 2 J. P. Solovej 3 1 Masaryk University Brno 2 University of Warsaw 3 University of Copenhagen Macroscopic Limits of Quantum Systems 70th Birthday of
More informationSemicircle law on short scales and delocalization for Wigner random matrices
Semicircle law on short scales and delocalization for Wigner random matrices László Erdős University of Munich Weizmann Institute, December 2007 Joint work with H.T. Yau (Harvard), B. Schlein (Munich)
More informationNeighbor Tables Long-Range Potentials
Neighbor Tables Long-Range Potentials Today we learn how we can handle long range potentials. Neighbor tables Long-range potential Ewald sums MSE485/PHY466/CSE485 1 Periodic distances Minimum Image Convention:
More informationArithmetic properties of harmonic weak Maass forms for some small half integral weights
Arithmetic properties of harmonic weak Maass forms for some small half integral weights Soon-Yi Kang (Joint work with Jeon and Kim) Kangwon National University 11-08-2015 Pure and Applied Number Theory
More informationRégularité des équations de Hamilton-Jacobi du premier ordre et applications aux jeux à champ moyen
Régularité des équations de Hamilton-Jacobi du premier ordre et applications aux jeux à champ moyen Daniela Tonon en collaboration avec P. Cardaliaguet et A. Porretta CEREMADE, Université Paris-Dauphine,
More informationBayesian inverse problems with Laplacian noise
Bayesian inverse problems with Laplacian noise Remo Kretschmann Faculty of Mathematics, University of Duisburg-Essen Applied Inverse Problems 2017, M27 Hangzhou, 1 June 2017 1 / 33 Outline 1 Inverse heat
More informationElectrons in a weak periodic potential
Electrons in a weak periodic potential Assumptions: 1. Static defect-free lattice perfectly periodic potential. 2. Weak potential perturbative effect on the free electron states. Perfect periodicity of
More informationInequalities of Babuška-Aziz and Friedrichs-Velte for differential forms
Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz
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 information31st Jerusalem Winter School in Theoretical Physics: Problem Set 2
31st Jerusalem Winter School in Theoretical Physics: Problem Set Contents Frank Verstraete: Quantum Information and Quantum Matter : 3 : Solution to Problem 9 7 Daniel Harlow: Black Holes and Quantum Information
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationREMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID
REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID DRAGOŞ IFTIMIE AND JAMES P. KELLIHER Abstract. In [Math. Ann. 336 (2006), 449-489] the authors consider the two dimensional
More informationIntroduction to Quantum Spin Systems
1 Introduction to Quantum Spin Systems Lecture 2 Sven Bachmann (standing in for Bruno Nachtergaele) Mathematics, UC Davis MAT290-25, CRN 30216, Winter 2011, 01/10/11 2 Basic Setup For concreteness, consider
More information