Poisson Solver, Pseudopotentials, Atomic Forces in the BigDFT code
|
|
- Jonathan Thompson
- 5 years ago
- Views:
Transcription
1 CECAM Tutorial on Wavelets in DFT, CECAM - LYON,, in the BigDFT code Kernel Luigi Genovese L_Sim - CEA Grenoble 28 November 2007
2 Outline, Kernel 1 The with Interpolating Scaling Functions in DFT for Interpolating scaling functions for electrostatic problems Caclulation of the Kernel Benchmarks and performances A solver for surfaces BC 2 GTH-HGH pseudopoentials in BigDFT Convenience of wavelet basis representation 3 Calculation of Treatment of the different terms
3 The Hartree potential In the DFT calculation in the Kohn-Sham formalism we have:, Find a set of orthonormal orbitals Ψ i (r) that minimizes: E = 1 2 N/2 i=1 Ψ i (r) 2 Ψ i (r)dr + 1 ρ(r)v H (r)dr 2 + E xc [ρ(r)] + V ext (r)ρ(r)dr Kernel where N/2 ρ(r) = 2 i=1 2 V H (r) = 4πρ(r) Ψ i (r)ψ i (r)
4 The in electronic structure calculation During the minimization procedure we need to perform, Kernel s equation Calculation of the self-consistent potential: 2 V H (r) = 4πρ(r) Such equation should be solved at each minimisation iteration. Need of having an efficient and accurate formalism. Plane waves approach The most common approach. Uses the fourier components f(x,y,z) = p x,p y,p z e 2πi ( px Lx ) py pz x+ y+ Ly Lz z f px,p y,p z The equation is algebraic in the Fourier coefficients
5 with plane waves treatment The Laplacian is diagonal in Plane waves representation, Kernel Immediate solution V px,p y,p z = 1 1 ( ) π 2 ( ) 2 ( ) 2 ρ px,p y,p z, p x py L x + L y + p z L z Characteristics Simple and fast, easy to parallelize (FFT) Automatically implement Periodic BC on a finite volume Do not fix the value of V px,p y,p z May result in problems for systems with other BC How to solve this equation for other BC?
6 Isolated BC: the Green function treatment Consider the equation for isolated BC. In this case the solution is given by, Kernel Green s function for (kernel) 2 1 r = 4πδ(r) = V(r) = An unambiguous solution dr ρ(r ) r r Such prescription is unique and is compatible with. How to implement it? Plane-wave based approaches: Truncated kernel, numerically or analytically Screening Functions Approximate treatment, large box required. Need of an accurate and efficient algorithm
7 with Interpolating Scaling Functions A convenient basis for an electrostatic problem, Kernel Interpolating Scaling Functions A set of localised functions centered on the nodes of a uniform mesh 1.5 ϕ j (x) = ϕ 0 (x j) Undergo multiscale relation ϕ(x) = m h j j= m}{{} filters ϕ(2x j) LEAST ASYMMETRIC DAUBECHIES-16 scaling function wavelet The expansion coefficients are the real space values ρ jx,j y,j z = ρ(h x j x,h y j y,h z j z ) Represents exactly an order m polynomial The first m discrete and continuous moments coincide
8 A finite three-dimensional convolution, Kernel The expression of the potential in this basis is thus intuitive: V(i) = K ij ρ j j Where the central object is the Kernel in the K ij = K i j, K i = K ( r )ϕ i (r)dr, K (r) = 1 r Values of the potential are obtained via a convolution V(i) = K i j ρ j j It can be treated via a zero-padded FFT algorithm Exact, easy to parallelize (FFT) For a box of N 3 points, it reduces the scaling from O(N 6 ) to O(N logn)
9 Characteristics of the approach This approach is: Explicit, guarantees the good BC, Kernel Real space-based, immediate interpretation of the expansion coefficients Can be combined with other real-space treatments of the density (e.g. XC) combined with ABINIT XC routines Can be used independently from the DFT code Preserves the first m (multipole) moments of the electrostatic potential Requires only the evaluation of the kernel We need to evaluate N 3 integrals K i = K ( r )ϕ i (r)dr
10 Gaussian tensor product decomposition It can be shown that (Beylkin et al.), Kernel Approximation with gaussians 1 r ω k e p k r 2 k with k = 1,,89, p k, ω k suitably chosen Accuracy of 10 8 for r (10 9,1) We can rescale for R (0,L). The computational cost is reduced N 3 89 N K j = 89 k=1 ω k K jx (p k )K jy (p k )K jz (p k ) K j (p) = ϕ 0 (x)e p(x j)2 dx
11 Other properties of the scaling functions The computational cost is reduced N 3 89 N., Moreover The scaling property of the interpolets ϕ 0 (x/2) = h j ϕ j (x), j Implies similar condition of the one-dimensional function Kernel K j (4p) = 1 2 j h j K 2i j (p). Thus we can evaluate the integrals for low p (not too sharp gaussians), then rescaling. ISF properties allows us to gain in accuracy
12 Characteristics Very fast with moderate memory occupation: Elapsed Time for a grid on a Cray XT3, proc s More precise than other existing s e-04 Kernel 1e-05 1e-06 Max Error 1e-07 1e-08 Hockney Tuckerman 1e-09 8th 14th 20th 30th 1e-10 40th 50th 60th 100th 1e Grid step
13 Characteristics of the, Kernel In summary, we have developed a technique Free boundary conditions Very high accuracy Good computational performance, easy to parallelize Can be used also in other contexts and/or combined with other treatment (e.g. XC) Coupled with ABINIT XC routines In BigDFT code, for big systems it represents a small amount (typical values are less than 2-3% for big systems) of the overall computation L. Genovese, T. Deutsch, A. Neelov, S. Goedecker, G. Beylkin Efficient solution of s equation with free boundary conditions, [arxiv: cond-mat/ ], J. Chem. Phys. 125, (2006)
14 The Surfaces boundary conditions The same formalism can be applied to other BC, Kernel A domain isolated in one direction (say y) and periodic in x and z, with periods L x and L z. A function f which lives in such a domain can be expanded as px pz f(x,y,z) = e 2πi( x+ z) Lx Lz f px,p z (y) p x,p z without any loss of generality. Mixed representation For such functions the s equation become ( 2 y µ 2 p x,p z ) Vpx,p z (y) = ρ px,p z (y), where µ 2 p x,p z = 4π 2 (p x /L x ) 2 + (p z /L z ) 2.
15 A Green s function formalism, Kernel The Green s function for the one-dimensional Helmoltz equation can be used ( 2 y µ 2) G(µ;y) = δ(y) ; { 1 2µ G(µ;y) = e µ y µ > 0 1 y µ =, 2 0 the components of the potential can be carried out: A Green s function for each Fourier component V px,p z (y) = dy G(µ px p z ;y y )ρ px,p z (y ). We can use Interpolating Scaling Functions for the Isolated direction.
16 with Mixed representation, Kernel By defining the one-dimensional kernels K (µ px,p z ;j) = G(µ px p z ;y)ϕ j (y)dy We can have a treatment similar to the case 1-dim convolution for each reciprocal space component V px,p z (i) = K (µ px,p z ;i j)ρ px,p z (j) j From mixed representation to full real space The calculation can be performed with a (semi-)zero-padded FFT algorithm. The I/O are function of the real domain.
17 Speed-up the calculation The calculation of the kernel can be improved, Kernel Speed The analytic form of the Green functions allows for recursion relations improvement in speed of a factor of (N z + 1)/3 Accuracy The scaling relation will again help in the kernel K (2µ;j) = 1 2 j h j K (2µ,2i j) A Free-Lunch case We have the same advantages of the treatment with an explicit formalism conceived for surfaces BC
18 for surface boundary conditions Elapsed Time on a Cray XT3, grid, #proc sec More precise than other treatments e-04 Kernel Absolute relative error 1e-06 1e-08 1e-10 1e-12 Mortensen 8th 16th 1e-14 24th 40th h 8 curve 1e grid spacing L. Genovese, T. Deutsch, S. Goedecker Efficient and accurate three dimensinal solver for surface problems, [arxiv: cond-mat/ ], J. Chem. Phys. 127, (2007)
19 A solver for surface problems, Kernel Like the case We developed a technique Accurate and fast, easy to parallelize Can be applied both in real or reciprocal space codes Explicit treatment No supercell or screening functions More precise than other existing approaches Allows comparisons betweeen different backgrounds (charged systems) Example of the plane capacitor: Periodic Hockney Our approach V V V y y y x x x
20 Pseudopotential treatment Smooth the wavefunction, Kernel Eliminate the rapid variations of the valence wavefunction in the core region. Thanks to them, only two levels of resolutions are enough for good accuracy in BigDFT The pseudopotential used are of GTH-HGH type V loc (r) = ( ) [ ( ) ] eff erf(r) + e 1 r 2 2 r r l P r r r l V nonloc = l ij l m= l h (l) ij p l,m i r p l,m i = Y l,m (θ,ϕ)f (l) i (r l )e 1 2 p l,m j ( ) r r l
21 of the pseudopotential operators These terms act differently on the energy calculation, Kernel Different treatment E loc = R V loc (r)ρ(r)dr V loc (r) can be expressed in interpolating scaling functions basis (real space values) E nonloc = l ij l m= l h(l) ij ψ p l,m i p l,m j ψ The important quantities are scalar products The general term which appear in the application of the nonlocal operator is of the form drψ(r)p(r) where P(r) is a gaussian times a separable polynomial The projector is a tensor product P(r) = P x (x)p y (y)p z (z)
22 Agrees with the wavefunction representation, Kernel For a wavefunction expressed in Daubechies wavelets ψ(r) = i x,i y,i z c ix,i y,i z φ ix (x)φ iy (y)φ iz (z) the calculation of the scalar product can be decomposed From N 3 to 3N calculations dr ψ(r)p(r) = i x,i y,i z c ix,i y,i z φ ix (x)p x (x)dx φ iy (y)p y (y)dy φ iz (z)p z (z)dz The numbers R φ ix (x)p x (x) are the expansion coefficients in Daubechies basis of the projectors. Treatment simple and accurate The application of the nonlocal projectors is easy and exact in the basis representation (no intrinsic approx.)
23 Suppose we find the ground state energy of an electronic system. The forces on an atom in the position R are, Kernel Gradient of the energy F i = de dr i The DFT Hamiltonian is given in the Bohr-Oppenheimer approximation. For this approach it is valid the Feynman-Hellmann theorem de = ψ H dr i R i ψ for the normalised ground state eigenfunction Equivalently, for this case we can write F i = E ext R i
24 Which terms depend explicitly of the atomic positions? The external ionic energy:, E ion = 1 2 a b a b R a R b Kernel The local pseudopotential energy: E loc = drv loc (r,{r a })ρ(r) The nonlocal (separable) one: E nonloc = h i,j ψ p i (R) p j (R) ψ i,j Let us see the procedure used for each term separately
25 The calculation of each term Immediate ionic force calculation The contribution coming from the ionic energy is analytical, Kernel Local pseudopotentials The local pseudopotential term can be separated in two parts: V loc (r) = ( ) [ ( ) ] eff erf(r) + e 1 r 2 2 r r l P = V e + V g r r r l where V g is a gaussian times a polynomial V e is the solution of the isolated s equation where ρ e e V e = 4πρ e ( ) r r l is a gaussian
26 Local pseudopotential energy contribution The contribution from the local term is thus, Kernel Local pseudopotentials F a = a V g ρ + a = a V g ρ + a = a V g ρ + drv e (r)ρ(r) drdr ρ e (r ) r r ρ(r) dr a ρ e (r)v H (r) This approach Preserves locality Easy to parallelize Is of O(n) scaling
27 Contribution from the Nonlocal pseudopotentials, Kernel Same treatment of the energy calculation The projectors of the separable part are gaussian times polynomials (spherical harmonics). Their derivatives are still gaussian times polynomials. For each projector its derivative p(r) p R is a function of the same kind and can thus be expressed with the same technique. We calculate the forces contribution by applying the Leibnitz rule to the energy expression
28 In Summary We implemented the calculation of, where E ext = E I + E L + E S F i = E ext R i Ionic E I = 1 2 i j i j R ij Local E L = R dr V loc (r)ρ el (r) Separable E S = o n o ψ o H S ψ o Kernel Each term is calculated with a different method. The complete implementation Preserves the locality (easy to parallelize) Compatible with the adaptivity The calculation can be performed with arbitrary accuracy (no intrinsic approximations)
Accuracy benchmarking of DFT results, domain libraries for electrostatics, hybrid functional and solvation
Accuracy benchmarking of DFT results, domain libraries for electrostatics, hybrid functional and solvation Stefan Goedecker Stefan.Goedecker@unibas.ch http://comphys.unibas.ch/ Multi-wavelets: High accuracy
More informationBigDFT tutorial. BigDFT tutorial FIRST YARMOUK SCHOOL. Thierry Deutsch. November, L_Sim - CEA Grenoble. Introduction DFT.
Big tutorial FIRST YARMOUK SCHOOL Big tutorial Thierry Deutsch L_Sim - CEA Grenoble November, 2. 2010 Outline 1 2 Density Functional Theory (quick view) 3 4 A basis for nanosciences: the Big project STREP
More informationWavelets for density functional calculations: Four families and three. applications
Wavelets for density functional calculations: Four families and three Haar wavelets Daubechies wavelets: BigDFT code applications Stefan Goedecker Stefan.Goedecker@unibas.ch http://comphys.unibas.ch/ Interpolating
More informationPseudopotential generation and test by the ld1.x atomic code: an introduction
and test by the ld1.x atomic code: an introduction SISSA and DEMOCRITOS Trieste (Italy) Outline 1 2 3 Spherical symmetry - I The Kohn and Sham (KS) equation is (in atomic units): [ 1 ] 2 2 + V ext (r)
More informationDaubechies wavelets as a basis set for density functional pseudopotential calculations
Daubechies wavelets as a basis set for density functional pseudopotential calculations Luigi Genovese, Alexey Neelov, 2 Stefan Goedecker, 2 Thierry Deutsch, Seyed Alireza Ghasemi, 2 Alexander Willand,
More informationElectronic Structure Calculations, Density Functional Theory and its Modern Implementations
Tutoriel Big RENOBLE Electronic Structure Calculations, Density Functional Theory and its Modern Implementations Thierry Deutsch L_Sim - CEA renoble 19 October 2011 Outline 1 of Atomistic calculations
More informationCharge Analysis: Atoms in Molecules
Daubechies Wavelets in Electronic Structure Calculation: BigDFT Code Tutorial CECAM - GRENOBLE : Atoms in Molecules Ali Sadeghi Basel University 21 November 2011 An output of electronic structure calculations
More informationIntroduction to density functional perturbation theory for lattice dynamics
Introduction to density functional perturbation theory for lattice dynamics SISSA and DEMOCRITOS Trieste (Italy) Outline 1 Lattice dynamic of a solid: phonons Description of a solid Equations of motion
More informationAlgorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part I Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
More informationIntroduction to Electronic Structure Calculations with BigDFT
Multiscale Modelling Methods for Applications in Materials Science CECAM JÜLICH, GERMANY run to Electronic Structure Calculations with Thierry Deutsch, Damien Caliste, Luigi Genovese L_Sim - CEA Grenoble
More informationCP2K: the gaussian plane wave (GPW) method
CP2K: the gaussian plane wave (GPW) method Basis sets and Kohn-Sham energy calculation R. Vuilleumier Département de chimie Ecole normale supérieure Paris Tutorial CPMD-CP2K CPMD and CP2K CPMD CP2K http://www.cpmd.org
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 informationAdaptive and Localized Basis Functions for Linear Scaling, Large Systems and Complex QM Simulations
Max Conference on the Materials Design Ecosystem at the Exascale: High-Performance and High-Throughput Computing TRIESTE Adaptive and Localized Basis Functions for Linear Scaling, Large Systems and Complex
More informationReinhold Schneider, MATHEON TU Berlin. Direct Minimization for effective single particle models (DFT)
Direct Minimization for effective single particle models (DFT) Partners joint work with: J. Blauert, T. Rohwedder (TU Berlin), A. Neelov (U Basel) joint EU NEST project BigDFT together with Dr. Thierry
More informationOpportunities from Accurate and Efficient Density Functional Theory Calculations for Large Systems
Seminar CENTRE FOR PREDICTIVE MODELLING, WARWICK Opportunities from Accurate and Efficient Density Functional Theory Calculations for Large Systems Luigi Genovese L_Sim CEA Grenoble October 30, 2017 http://bigdft.org
More information1. Hydrogen atom in a box
1. Hydrogen atom in a box Recall H atom problem, V(r) = -1/r e r exact answer solved by expanding in Gaussian basis set, had to solve secular matrix involving matrix elements of basis functions place atom
More informationRepresentation of operators in MADNESS (Multiresolution ADaptive Numerical Environment for Scientific Simulation)
Representation of operators in MADNESS (Multiresolution ADaptive Numerical Environment for Scientific Simulation) Gregory Beylkin University of Colorado at Boulder IPAM: Big Data Meets Computation workshop
More informationIntroduction to DFTB. Marcus Elstner. July 28, 2006
Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD
More informationThe Plane-Wave Pseudopotential Method
Hands-on Workshop on Density Functional Theory and Beyond: Computational Materials Science for Real Materials Trieste, August 6-15, 2013 The Plane-Wave Pseudopotential Method Ralph Gebauer ICTP, Trieste
More informationA Tutorial on Wavelets and their Applications. Martin J. Mohlenkamp
A Tutorial on Wavelets and their Applications Martin J. Mohlenkamp University of Colorado at Boulder Department of Applied Mathematics mjm@colorado.edu This tutorial is designed for people with little
More informationClassical Field Theory: Electrostatics-Magnetostatics
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 1-5 Electrostatics The behavior of an electrostatic field can be described
More informationChapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set
Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless
More informationThe Projector Augmented Wave method
The Projector Augmented Wave method Advantages of PAW. The theory. Approximations. Convergence. 1 The PAW method is... What is PAW? A technique for doing DFT calculations efficiently and accurately. An
More informationarxiv: v1 [physics.comp-ph] 22 Nov 2012
A Customized 3D GPU Poisson Solver for Free BCs Nazim Dugan a, Luigi Genovese b, Stefan Goedecker a, a Department of Physics, University of Basel, Klingelbergstr. 82, 4056 Basel, Switzerland b Laboratoire
More informationDFT / SIESTA algorithms
DFT / SIESTA algorithms Javier Junquera José M. Soler References http://siesta.icmab.es Documentation Tutorials Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang
More informationWhy use pseudo potentials?
Pseudo potentials Why use pseudo potentials? Reduction of basis set size effective speedup of calculation Reduction of number of electrons reduces the number of degrees of freedom For example in Pt: 10
More informationOn the adaptive finite element analysis of the Kohn-Sham equations
On the adaptive finite element analysis of the Kohn-Sham equations Denis Davydov, Toby Young, Paul Steinmann Denis Davydov, LTM, Erlangen, Germany August 2015 Denis Davydov, LTM, Erlangen, Germany College
More informationarxiv: v1 [physics.comp-ph] 22 Oct 2015
arxiv:1510.06489v1 [physics.comp-ph] 22 Oct 2015 Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework II: Force, vibration, and molecular dynamics calculations
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 informationPlane waves, pseudopotentials and PAW. X. Gonze Université catholique de Louvain, Louvain-la-neuve, Belgium
Plane waves, pseudopotentials and PAW X. Gonze Université catholique de Louvain, Louvain-la-neuve, Belgium 1 Basic equations in DFT Solve self-consistently the Kohn-Sham equation H ψ n = ε n ψ n!!! ρ(r
More informationDensity Functional Theory
Density Functional Theory Iain Bethune EPCC ibethune@epcc.ed.ac.uk Overview Background Classical Atomistic Simulation Essential Quantum Mechanics DFT: Approximations and Theory DFT: Implementation using
More informationA Fast N-Body Solver for the Poisson(-Boltzmann) Equation
A Fast N-Body Solver for the Poisson(-Boltzmann) Equation Robert D. Skeel Departments of Computer Science (and Mathematics) Purdue University http://bionum.cs.purdue.edu/2008december.pdf 1 Thesis Poisson(-Boltzmann)
More informationDFT in practice : Part II. Ersen Mete
pseudopotentials Department of Physics Balıkesir University, Balıkesir - Turkey August 13, 2009 - NanoDFT 09, İzmir Institute of Technology, İzmir Outline Pseudopotentials Basic Ideas Norm-conserving pseudopotentials
More informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationFirst-principles Molecular Dynamics Simulations
First-principles Molecular Dynamics Simulations François Gygi University of California, Davis fgygi@ucdavis.edu http://eslab.ucdavis.edu http://www.quantum-simulation.org MICCoM Computational School, Jul
More informationDGDFT: A Massively Parallel Method for Large Scale Density Functional Theory Calculations
DGDFT: A Massively Parallel Method for Large Scale Density Functional Theory Calculations The recently developed discontinuous Galerkin density functional theory (DGDFT)[21] aims at reducing the number
More informationThe Linearized Augmented Planewave (LAPW) Method
The Linearized Augmented Planewave (LAPW) Method David J. Singh Oak Ridge National Laboratory E T [ ]=T s [ ]+E ei [ ]+E H [ ]+E xc [ ]+E ii {T s +V ks [,r]} I (r)= i i (r) Need tools that are reliable
More informationSpring College on Computational Nanoscience May Variational Principles, the Hellmann-Feynman Theorem, Density Functional Theor
2145-25 Spring College on Computational Nanoscience 17-28 May 2010 Variational Principles, the Hellmann-Feynman Theorem, Density Functional Theor Stefano BARONI SISSA & CNR-IOM DEMOCRITOS Simulation Center
More informationLinear Diffusion and Image Processing. Outline
Outline Linear Diffusion and Image Processing Fourier Transform Convolution Image Restoration: Linear Filtering Diffusion Processes for Noise Filtering linear scale space theory Gauss-Laplace pyramid for
More informationElectron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals
Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals Richard M. Martin UIUC Lecture at Summer School Hands-on introduction to Electronic Structure Materials Computation Center University
More informationDirect Minimization in Density Functional Theory
Direct Minimization in Density Functional Theory FOCM Hongkong 2008 Partners joint paper with: J. Blauert, T. Rohwedder (TU Berlin), A. Neelov (U Basel) joint EU NEST project BigDFT together with Dr. Thierry
More informationLarge Scale Electronic Structure Calculations
Large Scale Electronic Structure Calculations Jürg Hutter University of Zurich 8. September, 2008 / Speedup08 CP2K Program System GNU General Public License Community Developers Platform on "Berlios" (cp2k.berlios.de)
More informationRepresentation theory and quantum mechanics tutorial Spin and the hydrogen atom
Representation theory and quantum mechanics tutorial Spin and the hydrogen atom Justin Campbell August 3, 2017 1 Representations of SU 2 and SO 3 (R) 1.1 The following observation is long overdue. Proposition
More informationTight-Binding Model of Electronic Structures
Tight-Binding Model of Electronic Structures Consider a collection of N atoms. The electronic structure of this system refers to its electronic wave function and the description of how it is related to
More informationElectronic band structure, sx-lda, Hybrid DFT, LDA+U and all that. Keith Refson STFC Rutherford Appleton Laboratory
Electronic band structure, sx-lda, Hybrid DFT, LDA+U and all that Keith Refson STFC Rutherford Appleton Laboratory LDA/GGA DFT is good but... Naive LDA/GGA calculation severely underestimates band-gaps.
More information1. Thomas-Fermi method
1. Thomas-Fermi method We consider a system of N electrons in a stationary state, that would obey the stationary Schrödinger equation: h i m + 1 v(r i,r j ) Ψ(r 1,...,r N ) = E i Ψ(r 1,...,r N ). (1.1)
More informationSelf Consistent Cycle
Self Consistent Cycle Step 0 : defining your system namelist SYSTEM How to specify the System All periodic systems can be specified by a Bravais Lattice and and atomic basis How to specify the Bravais
More informationSolutions to selected problems from Giuliani, Vignale : Quantum Theory of the Electron Liquid
Solutions to selected problems from Giuliani, Vignale : Quantum Theory of the Electron Liquid These problems have been solved as a part of the Independent study class with prof. Neepa Maitra. Compiled
More informationDunkl operators and Clifford algebras II
Translation operator for the Clifford Research Group Department of Mathematical Analysis Ghent University Hong Kong, March, 2011 Translation operator for the Hermite polynomials Translation operator for
More informationExpansion of 1/r potential in Legendre polynomials
Expansion of 1/r potential in Legendre polynomials In electrostatics and gravitation, we see scalar potentials of the form V = K d Take d = R r = R 2 2Rr cos θ + r 2 = R 1 2 r R cos θ + r R )2 Use h =
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
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 informationModule 6 1. Density functional theory
Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft
More information1 r 2 sin 2 θ. This must be the case as we can see by the following argument + L2
PHYS 4 3. The momentum operator in three dimensions is p = i Therefore the momentum-squared operator is [ p 2 = 2 2 = 2 r 2 ) + r 2 r r r 2 sin θ We notice that this can be written as sin θ ) + θ θ r 2
More informationKey concepts in Density Functional Theory (I) Silvana Botti
From the many body problem to the Kohn-Sham scheme European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre
More informationReferences. Documentation Manuals Tutorials Publications
References http://siesta.icmab.es Documentation Manuals Tutorials Publications Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang atomic energy unit = 1 Hartree =
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 information( nonlinear constraints)
Wavelet Design & Applications Basic requirements: Admissibility (single constraint) Orthogonality ( nonlinear constraints) Sparse Representation Smooth functions well approx. by Fourier High-frequency
More informationarxiv: v1 [cond-mat.mtrl-sci] 29 Jan 2015
Interatomic potentials for ionic systems with density functional accuracy based on charge densities obtained by a neural network S. Alireza Ghasemi 1,, Albert Hofstetter 2, Santanu Saha 2 and Stefan Goedecker
More informationPhani Motamarri and Vikram Gavini Department of Mechanical Engineering, University of Michigan, Ann Arbor, MI 48109, USA
A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization Phani Motamarri and Vikram Gavini Department
More informationMore On Carbon Monoxide
More On Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results Jerry Gilfoyle The Configurations of CO 1 / 26 More On Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results Jerry Gilfoyle The Configurations
More informationDirect Learning: Linear Classification. Donglin Zeng, Department of Biostatistics, University of North Carolina
Direct Learning: Linear Classification Logistic regression models for classification problem We consider two class problem: Y {0, 1}. The Bayes rule for the classification is I(P(Y = 1 X = x) > 1/2) so
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationarxiv: v1 [cond-mat.str-el] 18 Jul 2007
arxiv:0707.2704v1 [cond-mat.str-el] 18 Jul 2007 Determination of the Mott insulating transition by the multi-reference density functional theory 1. Introduction K. Kusakabe Graduate School of Engineering
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationDensity matrix functional theory vis-á-vis density functional theory
Density matrix functional theory vis-á-vis density functional theory 16.4.007 Ryan Requist Oleg Pankratov 1 Introduction Recently, there has been renewed interest in density matrix functional theory (DMFT)
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationPseudopotential methods for DFT calculations
Pseudopotential methods for DFT calculations Lorenzo Paulatto Scuola Internazionale Superiore di Studi Avanzati and CNR-INFM DEMOCRITOS National Simulation Center Tieste Italy July 9, 2008 Outline pseudopotential
More informationPolynomial Approximation: The Fourier System
Polynomial Approximation: The Fourier System Charles B. I. Chilaka CASA Seminar 17th October, 2007 Outline 1 Introduction and problem formulation 2 The continuous Fourier expansion 3 The discrete Fourier
More informationCHEM6085: Density Functional Theory
Lecture 11 CHEM6085: Density Functional Theory DFT for periodic crystalline solids C.-K. Skylaris 1 Electron in a one-dimensional periodic box (in atomic units) Schrödinger equation Energy eigenvalues
More information2. TranSIESTA 1. SIESTA. DFT In a Nutshell. Introduction to SIESTA. Boundary Conditions: Open systems. Greens functions and charge density
1. SIESTA DFT In a Nutshell Introduction to SIESTA Atomic Orbitals Capabilities Resources 2. TranSIESTA Transport in the Nanoscale - motivation Boundary Conditions: Open systems Greens functions and charge
More informationElectronic structure, plane waves and pseudopotentials
Electronic structure, plane waves and pseudopotentials P.J. Hasnip Spectroscopy Workshop 2009 We want to be able to predict what electrons and nuclei will do from first principles, without needing to know
More informationLectures notes. Rheology and Fluid Dynamics
ÉC O L E P O L Y T E C H N IQ U E FÉ DÉR A L E D E L A U S A N N E Christophe Ancey Laboratoire hydraulique environnementale (LHE) École Polytechnique Fédérale de Lausanne Écublens CH-05 Lausanne Lectures
More informationAlgorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part II Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
More informationProblem 1: A 3-D Spherical Well(10 Points)
Problem : A 3-D Spherical Well( Points) For this problem, consider a particle of mass m in a three-dimensional spherical potential well, V (r), given as, V = r a/2 V = W r > a/2. with W >. All of the following
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationSpin in Density-Functional Theory
Spin in Density-Functional Theory Christoph R. Jacob 1 and Markus Reiher 2 1 Karlsruhe Institute of Technology (KIT), Center for Functional Nanostructures and Institute of Physical Chemistry, Wolfgang-Gaede-Straße
More informationNotes on Density Functional Theory
Notes on Density Functional Theory 1 Basic Theorems The energy, E, of a system with a given Hamiltonian H is a functional of the (normalized, many-particle) wave function Ψ. We write this functional as
More informationELECTRONIC STRUCTURE CALCULATIONS FOR THE SOLID STATE PHYSICS
FROM RESEARCH TO INDUSTRY 32 ème forum ORAP 10 octobre 2013 Maison de la Simulation, Saclay, France ELECTRONIC STRUCTURE CALCULATIONS FOR THE SOLID STATE PHYSICS APPLICATION ON HPC, BLOCKING POINTS, Marc
More information20 The Hydrogen Atom. Ze2 r R (20.1) H( r, R) = h2 2m 2 r h2 2M 2 R
20 The Hydrogen Atom 1. We want to solve the time independent Schrödinger Equation for the hydrogen atom. 2. There are two particles in the system, an electron and a nucleus, and so we can write the Hamiltonian
More informationShort-time Fourier transform for quaternionic signals
Short-time Fourier transform for quaternionic signals Joint work with Y. Fu and U. Kähler P. Cerejeiras Departamento de Matemática Universidade de Aveiro pceres@ua.pt New Trends and Directions in Harmonic
More informationHilbert Space Methods for Reduced-Rank Gaussian Process Regression
Hilbert Space Methods for Reduced-Rank Gaussian Process Regression Arno Solin and Simo Särkkä Aalto University, Finland Workshop on Gaussian Process Approximation Copenhagen, Denmark, May 2015 Solin &
More informationThe Helmholtz theorem at last!
Problem. The Helmholtz theorem at last! Recall in class the Helmholtz theorem that says that if if E =0 then E can be written as E = φ () if B =0 then B can be written as B = A (2) (a) Let n be a unit
More informationAtoms 09 update-- start with single electron: H-atom
Atoms 09 update-- start with single electron: H-atom VII 33 x z φ θ e -1 y 3-D problem - free move in x, y, z - handy to change systems: Cartesian Spherical Coordinate (x, y, z) (r, θ, φ) Reason: V(r)
More informationLecture 11 Spin, orbital, and total angular momentum Mechanics. 1 Very brief background. 2 General properties of angular momentum operators
Lecture Spin, orbital, and total angular momentum 70.00 Mechanics Very brief background MATH-GA In 9, a famous experiment conducted by Otto Stern and Walther Gerlach, involving particles subject to a nonuniform
More informationTime-dependent density functional theory : direct computation of excitation energies
CECAM Tutorial Electronic excitations and spectroscopies : Theory and Codes Lyon 007 Time-dependent density functional theory : direct computation of excitation energies X. Gonze Université Catholique
More informationElectronic energy optimisation in ONETEP
Electronic energy optimisation in ONETEP Chris-Kriton Skylaris cks@soton.ac.uk 1 Outline 1. Kohn-Sham calculations Direct energy minimisation versus density mixing 2. ONETEP scheme: optimise both the density
More informationSpectral Processing. Misha Kazhdan
Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry
More informationElectronic State Calculation for Hydrogen Atom under Inhomogeneous Field
26 Special Issue Multiscale Simulations for Materials Research Report Electronic State Calculation for Hydrogen Atom under Inhomogeneous Field Shunsuke Yamakawa and Shi-aki Hyodo We discuss herein the
More information1.2 Rare earth atoms H Ψ=E Ψ, (1.2.1) where the non-relativistic Hamiltonian operator is. H = h2 2m. v ext (r i ) (1.2.
8 1. ELEMENTS OF RARE EARTH MAGNETISM 1.2 Rare earth atoms The starting point for the understanding of the magnetism of the rare earths is the description of the electronic states, particularly of the
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 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 informationIntegrated Computational Materials Engineering Education
Integrated Computational Materials Engineering Education Lecture on Density Functional Theory An Introduction Mark Asta Dept. of Materials Science and Engineering, University of California, Berkeley &
More informationElectrochemistry project, Chemistry Department, November Ab-initio Molecular Dynamics Simulation
Electrochemistry project, Chemistry Department, November 2006 Ab-initio Molecular Dynamics Simulation Outline Introduction Ab-initio concepts Total energy concepts Adsorption energy calculation Project
More informationWavelets in Scattering Calculations
Wavelets in Scattering Calculations W. P., Brian M. Kessler, Gerald L. Payne polyzou@uiowa.edu The University of Iowa Wavelets in Scattering Calculations p.1/43 What are Wavelets? Orthonormal basis functions.
More informationab initio Electronic Structure Calculations
ab initio Electronic Structure Calculations New scalability frontiers using the BG/L Supercomputer C. Bekas, A. Curioni and W. Andreoni IBM, Zurich Research Laboratory Rueschlikon 8803, Switzerland ab
More informationTime-Dependent Density-Functional Theory
Summer School on First Principles Calculations for Condensed Matter and Nanoscience August 21 September 3, 2005 Santa Barbara, California Time-Dependent Density-Functional Theory X. Gonze, Université Catholique
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More information