Large-scale real-space electronic structure calculations
|
|
- Philip Woods
- 5 years ago
- Views:
Transcription
1 Large-scale real-space electronic structure calculations YIP: Quasi-continuum reduction of field theories: A route to seamlessly bridge quantum and atomistic length-scales with continuum Grant no: FA PI: Vikram Gavini Slides: 32 Vikram Gavini and Phani Motamarri Department of Mechanical Engineering University of Michigan, Ann Arbor
2 Defective Crystals Defects play a crucial role in influencing a variety of materials properties mechanical, electronic, optical, chemical Dislocations Vacancies/Interstitials Interfaces/Surfaces TEM image of dislocation partial T.J. Balk, K.J. Hemker, Phil. Mag. A, 2001 Prismatic loops formed from vacancies, Giess et. al, Microsc Microanal, 2005 TEM image of Ni-Al interface, Mann et.al, J. Appl. Phys Metal Plasticity Renders the strength of materials to 1/1000 its theoretical strength Creep, Spall, Ageing, hardening due to radiation Phase stability, Energetics, Diffusion mediation, defect sources and sinks
3 Defective crystals: The challenge The energetics of defects: (i) core-energy; (ii) elastic energy The core of a defect is governed by electronic structure need electronic structure calculations! Defects result in a vast span of interacting length scales Electronic structure of the core (10-12 m) Complex rearrangements of atoms around the core (10-9 m) Long ranged elastic effects (10-6 m) But need single physics at all length scales! - No patching, seamless description Realistic defect concentration in materials is parts per million! Challenge : Need electronic structure calculations at macroscopic scales! (i) Development of computational techniques for large-scale electronic structure calculation that can explicitly treat systems up to 10,000 atoms (ii) Development of seamless coarse-graining schemes using adaptive numerical schemes
4 Kohn-Sham density-functional theory (KSDFT) The KSDFT energy functional is given by, Local density approximation (LDA) Non-local Classical electrostatic interaction energy : Computed in Fourier-space (reciprocalspace) in almost all DFT implementations
5 KSDFT Real-space formulation (Suryanarayana & Gavini et al. JMPS 2010, ) Non-local Green s function for Laplace operator Electrostatic interactions can be re-written locally as, Thus, (Regularized nuclear charges)
6 KSDFT Real-space formulation (Suryanarayana & Gavini et al. JMPS 2010, ) The saddle-point problem is given by, Define, Theorem : Proof : Sobolev embeddings; Poincaré inequality (Direct Method)
7 Kohn-Sham eigenvalue problem Consider the E-L equation corresponding to the variational problem: Self consistent iteration (Kohn-Sham map) To avoid charge-sloshing:
8 State of the art Solutions to Kohn-Sham Equations Fourier Space Formulations Key Features (plane-waves) Very efficient for periodic calculations Restrictive to periodic domains Provide only uniform spatial resolution Suitable only when the solution fields are smooth. Real Space Formulations (LCAO, FDM, FEM) Key Features (LCAO) Suitable for isolated systems Can handle both pseudopotential and all electron calculations Systematic convergence can not be ascertained Parallel scalability is a concern
9 KSDFT FE discretization Use finite-element basis for computing u 1 i=1 i=2 N 1 (r) N 2 (r) N 3 (r) i=1 i=2 r r Features of finite-element basis: 1. Unstructured coarse-graining 2. Complex geometries can be represented, and arbitrary boundary conditions can be imposed. 3. Systematic convergence 4. Ease of parallel implementation By changing the positioning of the nodes the spatial resolution of basis can be changed/adapted
10 KSDFT FE discretization Main Limitations: Previous attempts showed that the number of FE basis functions (linear) needed to obtain chemical accuracy is very large ~ 100,000-1,000,000 basis functions per atom. The finite-element discretization leads to a generalized eigenvalue problem, which is more challenging to solve than a standard eigenvalue problem Present Work: We demonstrate an efficient, scalable computational approach using adaptive higher-order finite-element discretization. We propose a linear scaling algorithm (in number of electrons) which treats both insulating and metallic systems on an equal footing.
11 KSDFT FE discretization Discrete eigenvalue problem: Transformation to a standard eigenvalue problem: Remark: denotes the projection of the Hamiltonian operator into a space spanned by Löwden orthonormalized finite-element basis
12 Can higher-order finite-elements do any better? Here, we investigate the viability and computational efficiency afforded by higher-order finite-element discretization in electronic structure calculations using density functional theory to answer the following questions: What is the numerical convergence rate for various orders of finite-element approximations in electronic structure calculations using DFT? What is the computational advantage derived by using higher-order finite element discretization in terms of the CPU time? We present some of the first studies which demonstrate the computational efficiency afforded by higher-order elements for Kohn-Sham DFT calculations.
13 Rate of convergence of the finite element approximation The study has been carried out by a suite of higher order elements: TET 10 (TETRAHEDRAL QUADRATIC ELEMENT) HEX 27 (TRI QUADRATIC HEXAHEDRAL ELEMENT) HEX 64 (TRI CUBIC HEXAHEDRAL ELEMENT) HEX 125 (TRI QUARTIC HEXAHEDRAL ELEMENT) HEX 64 SPECTRAL, HEX 125 SPECTRAL upto 10 th order (Lagrange Polynomials are constructed on Gauss-Lobatto Legendre Points for spectral elements) Elements have been tested against three types of problems: (a) CH 4 (b) Barium Cluster (35 atoms) (a) CH 4 : An all electron calculation (b) Barium Cluster: Pseudopotential calculation
14 Convergence rates (Motamarri et al. J. Comp. Phys (2013)) Barium Cluster Methane Molecule Optimal rate of convergence!
15 Computational efficiency of higher-order FE discretization Two Key questions: How do higher-order FE discretizations compare to lower order elements in computational efficiency? How do higher-order FE discretizations compare against plane-wave basis and Gaussian basis? Key ideas in improving computational efficiency: Developed a priori mesh adaption techniques Use of Gauss-Legendre-Lobatto quadrature rules for the overlap matrix in conjunction with Spectral FE discretization Developed a Chebyschev acceleration technique to directly compute the eigenspace
16 A priori mesh adaption (Motamarri et al. J. Comp. Phys (2013)) Error Estimate: Optimal mesh distribution:
17 Spectral FE and Gauss-Lobatto-Legendre quadrature Spectral-element basis functions: Constructed from Lagrange polynomials interpolated through nodes corresponding to the roots of the derivatives of the Legendre polynomials and boundary nodes (GLL points) Upon using a Gauss-Lobatto-Legendre quadrature rule, the quadrature points coincide with the FE nodes Remarks: Transformation to standard eigenvalue problem is trivial The reduced order quadrature rule is only employed for the computation of the overlap matrix, and the full Gauss quadrature is employed to compute the Hamiltonian matrix.
18 Eigen-space computation: Chebyschev acceleration (Motamarri et al. J. Comp. Phys (2013)) E Occupied eigen-space Unwanted eigen-space Occupied eigen-space -1 1
19 Numerical algorithm 1. Start with initial guess for electron density and the initial wavefunctions 2. Compute the discrete Hamiltonian using the input electron density 3. Compute the Chebyshev filtered basis : 4. Orthonormalize the basis and compute, the projected Hamiltonian into the subspace spanned by 5. Compute the Fermi-energy and the output electron density by diagonalizing and using the following equation 6. If, EXIT; else, compute new using a mixing scheme and go to (2).
20 Computational efficiency (Motamarri et al. J. Comp. Phys. 2013) Barium Cluster Methane Molecule
21 Aluminum clusters
22 All-electron calculations 100 atom Graphene sheet
23 Scalability
24 Computational complexity Complexity in each SCF iteration: M: Number of degrees of freedom N: Number of electrons (M N) Chebyshev filtering procedure : O(MN) Orthonormalization of Chebyschev filtered vectors : O(MN 2 ) Diagonalization of the projected Hamiltonian: O(N 3 ) Cubic Scaling in N!
25 Subspace projection technique (Motamarri & Gavini, Phys. Rev. B (2014)) Key features of the proposed method: Chebyshev filtering to generate the approximate occupied subspace Construct non-orthogonal localized basis functions that span the same space & truncate these localized functions beyond a prescribed tolerance Subsequently localized basis functions have a compact support Use an adaptive tolerance for the truncation: truncation tolerance tied to the error in the SCF iteration; ensures strict control on accuracy Project Hamiltonian into the occupied subspace expressed in the nonorthogonal basis Fermi-operator expansion of the projected Hamiltonian to estimate Fermienergy and compute the electron density Avoids diagonalization of the Hamiltonian to compute orbital occupancies Applicable for both metallic and insulating systems Applicable for both pseudopotenial and all-electron calculations
26 Subspace projection approach: Key ideas Eigen-space from Chebyschev filtering Localized basis spanning eigen-space (Garcia-Cervera et al.) Project Hamiltonian in the localized basis: Remarks: Locality of is sparse and can be computed in O(N) complexity can be computed using Newton-Schultz algorithm which has O(N) complexity Finally, can be computed in O(N), if is sparse.
27 Subspace projection approach: Key ideas Computation of electron density Recall: No diagonalization No knowledge of eigenvalues and eigenvectors Compute density matrix instead: The electron density is the diagonal of the density matrix Fermi-operator expansion techniques can be employed to compute the density matrix: Challenge: ; spectral width of the discrete Hamiltonian is about !
28 Subspace projection approach: Key ideas Fermi-operator expansion of the projected Hamiltonian: Compute the density matrix using the projected Hamiltonian in the nonorthogonal localized basis The spectral width of the projected Hamiltonian is ~ O(10) and thus can efficiently employ the Fermi-operator expansion This approach treats both insulating and metallic systems on equal footing This approach is applicable for both pseudopotential calculations and allelectron calculations
29 Case study: Al nano-clusters (3x3x3 9x9x9) Pseudopotential calculations Total computational time Subspace projection scaling: O(N 1.46 ) Electron density contours on the midplane of the 9x9x9 nano-cluster Accuracy of subspace projection method commensurate with chemical accuracy
30 Case study: Alkane chains (C 33 H 68 C 2350 H 4702 ) Pseudopotential calculations Numerical accuracy Total computational time Isocontours of alkane chain C 900 H 1802 Subspace projection scaling: O(N 1.18 )
31 Case study: Si nano-clusters (1x1x1 3x3x3) All-electron calculations Total computational time Subspace projection scaling: O(N 1.85 ) Electron density contours on the midplane of the 3x3x3 Si nano-cluster Accuracy of subspace projection method commensurate with chemical accuracy
32 Ongoing/future work Real-space DFT-FE: Develop robust geometry optimization techniques Incorporate more advanced exchange-correlation functionals (beyond LDA, GGA) Exploring tensor structured techniques and low rank approximations in conjunction with real-space formulation Coarse-graining KSDFT: Localization of the wavefunctions is key for extending the coarse-graining ideas O(N) formulations: Non-orthogonal localized orbitals QC-KSDFT: localization -> predictor-corrector approach -> QC Electronic structure calculations at macroscopic scales with Kohn-Sham DFT will enable a quantum-mechanically accurate study of defects in materials
33 Concluding remarks Developed real-space formulation for Kohn-Sham DFT Reformulation of electrostatics as a local variational problem Mathematical analysis Finite-element discretization of Kohn-Sham DFT & Numerical algorithms Optimal rates of convergence Spectral elements in conjunction with GLL quadratures (for overlap matrix) Chebyshev filtering to directly compute the eigenspace Large-scale calculations possible Algorithms exhibit good scalability Development of a linear-scaling algorithm Localized basis spanning the Chebyshev filtered subspace Project of Hamiltonian into subspace Use Fermi-operator expansion on the projected Hamiltonian
34 Publications from AFOSR Comp Math program support P. Motamarri, T. Blesgen, V. Gavini, Tucker-tensor algorithm for large-scale Kohn-Sham density functional theory calculations, submitted (2015). P. Motamarri, V. Gavini, Phys. Rev. B 90, (2014). M. Iyer, T.M. Pollock, V. Gavini, Phys. Rev. B 89, (2014). P. Motamarri, M.R. Nowak, K. Leiter, J. Knap, V. Gavini, J. Comput. Phys. 253, (2013). (Previous award) T. Blesgen, V. Gavini, V. Khoromskaia, J. Comput. Phys. 231, (2012). P. Mottamari, M. Iyer, J. Knap, V. Gavini, J. Comput. Phys. 231, (2012). V. Gavini, L. Liu, J. Mech. Phys. Solids 59, (2011). M. Iyer, V. Gavini, J. Mech. Phys. Solids 59, (2011). P. Suryanarayana, P., V. Gavini, T. Blesgen, K. Bhattacharya, M. Ortiz, J. Mech. Phys. Solids 58, (2010). B. Radhakrishnan, V. Gavini, Phys. Rev. B 82, (2010). V. Gavini, Proc. Roy. Soc. A. 465, (2009).
35 THANK YOU!
Electronic-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 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 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 informationChapter 1. Electronic Structure Calculations at Macroscopic Scales using Orbital-free DFT
Chapter 1 Electronic Structure Calculations at Macroscopic Scales using Orbital-free DFT Balachandran G. Radhakrishnan and Vikram Gavini Department of Mechanical Engineering, University of Michigan, Ann
More informationRole of the defect-core in energetics of vacancies arxiv: v1 [cond-mat.mtrl-sci] 4 Aug 2009
Role of the defect-core in energetics of vacancies arxiv:98.499v1 [cond-mat.mtrl-sci] 4 Aug 29 By Vikram Gavini Department of Mechanical Engineering, University of Michigan, Ann Arbor, USA Electronic structure
More informationIntroduction to Computational Materials Science
Introduction to Computational Materials Science Fundamentals to Applications RICHARD LESAR lowa State University.CAMBRIDGE ::: UNIVERSITY PRESS CONTENTS Preface 1 Introduction to materials modeling and
More informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
More informationELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers
ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers Victor Yu and the ELSI team Department of Mechanical Engineering & Materials Science Duke University Kohn-Sham Density-Functional
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 informationALMA: All-scale predictive design of heat management material structures
ALMA: All-scale predictive design of heat management material structures Version Date: 2015.11.13. Last updated 2015.12.02 Purpose of this document: Definition of a data organisation that is applicable
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationBlock Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems
Mitglied der Helmholtz-Gemeinschaft Block Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems Birkbeck University, London, June the 29th 2012 Edoardo Di Napoli Motivation and Goals
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 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 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 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 informationApplied Numerical Mathematics
Applied Numerical Mathematics 131 (2018) 1 15 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum The spectral-galerkin approximation of nonlinear eigenvalue
More informationProgram : M.A./M.Sc. (Mathematics) M.A./M.Sc. (Final) Paper Code:MT-08 Numerical Analysis Section A (Very Short Answers Questions)
Program : M../M.Sc. (Mathematics) M../M.Sc. (Final) Paper Code:MT-08 Numerical nalysis Section (Very Short nswers Questions) 1. Write two examples of transcendental equations. (i) x 3 + sin x = 0 (ii)
More informationFEAST eigenvalue algorithm and solver: review and perspectives
FEAST eigenvalue algorithm and solver: review and perspectives Eric Polizzi Department of Electrical and Computer Engineering University of Masachusetts, Amherst, USA Sparse Days, CERFACS, June 25, 2012
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationAccuracy 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 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 informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
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 informationBand calculations: Theory and Applications
Band calculations: Theory and Applications Lecture 2: Different approximations for the exchange-correlation correlation functional in DFT Local density approximation () Generalized gradient approximation
More informationSolid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationNumerical Mathematics
Alfio Quarteroni Riccardo Sacco Fausto Saleri Numerical Mathematics Second Edition With 135 Figures and 45 Tables 421 Springer Contents Part I Getting Started 1 Foundations of Matrix Analysis 3 1.1 Vector
More informationPart IB - Easter Term 2003 Numerical Analysis I
Part IB - Easter Term 2003 Numerical Analysis I 1. Course description Here is an approximative content of the course 1. LU factorization Introduction. Gaussian elimination. LU factorization. Pivoting.
More informationMany electrons: Density functional theory Part II. Bedřich Velický VI.
Many electrons: Density functional theory Part II. Bedřich Velický velicky@karlov.mff.cuni.cz VI. NEVF 514 Surface Physics Winter Term 013-014 Troja 1 st November 013 This class is the second devoted to
More informationQuantum wires, orthogonal polynomials and Diophantine approximation
Quantum wires, orthogonal polynomials and Diophantine approximation Introduction Quantum Mechanics (QM) is a linear theory Main idea behind Quantum Information (QI): use the superposition principle of
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationCheat Sheet for MATH461
Cheat Sheet for MATH46 Here is the stuff you really need to remember for the exams Linear systems Ax = b Problem: We consider a linear system of m equations for n unknowns x,,x n : For a given matrix A
More informationOn the Numerics of 3D Kohn-Sham System
On the of 3D Kohn-Sham System Weierstrass Institute for Applied Analysis and Stochastics Workshop: Mathematical Models for Transport in Macroscopic and Mesoscopic Systems, Feb. 2008 On the of 3D Kohn-Sham
More informationarxiv: v1 [math.na] 18 Mar 2012
DECAY PROPERTIES OF SPECTRAL PROJECTORS WITH APPLICATIONS TO ELECTRONIC STRUCTURE MICHELE BENZI, PAOLA BOITO, AND NADER RAZOUK arxiv:1203.3953v1 [math.na] 18 Mar 2012 Abstract. Motivated by applications
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More information7/29/2014. Electronic Structure. Electrons in Momentum Space. Electron Density Matrices FKF FKF. Ulrich Wedig
Electron Density Matrices Density matrices Γ, an alternative to the wavefunction Ψ, for the description of a quantum system Electronic Structure The N-particle density matrix Electrons in Momentum Space
More informationQuantum Chemistry. NC State University. Lecture 5. The electronic structure of molecules Absorption spectroscopy Fluorescence spectroscopy
Quantum Chemistry Lecture 5 The electronic structure of molecules Absorption spectroscopy Fluorescence spectroscopy NC State University 3.5 Selective absorption and emission by atmospheric gases (source:
More informationA Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions
Lin Lin A Posteriori DG using Non-Polynomial Basis 1 A Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions Lin Lin Department of Mathematics, UC Berkeley;
More informationQuasi-continuum orbital-free density-functional. theory: A route to multi-million atom non-periodic. DFT calculation
Quasi-continuum orbital-free density-functional theory: A route to multi-million atom non-periodic DFT calculation Vikram Gavini, Kaushik Bhattacharya, Michael Ortiz Division of Engineering and Applied
More informationhomogeneous 71 hyperplane 10 hyperplane 34 hyperplane 69 identity map 171 identity map 186 identity map 206 identity matrix 110 identity matrix 45
address 12 adjoint matrix 118 alternating 112 alternating 203 angle 159 angle 33 angle 60 area 120 associative 180 augmented matrix 11 axes 5 Axiom of Choice 153 basis 178 basis 210 basis 74 basis test
More informationCompressed representation of Kohn-Sham orbitals via selected columns of the density matrix
Lin Lin Compressed Kohn-Sham Orbitals 1 Compressed representation of Kohn-Sham orbitals via selected columns of the density matrix Lin Lin Department of Mathematics, UC Berkeley; Computational Research
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 informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
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 informationDFT in practice. Sergey V. Levchenko. Fritz-Haber-Institut der MPG, Berlin, DE
DFT in practice Sergey V. Levchenko Fritz-Haber-Institut der MPG, Berlin, DE Outline From fundamental theory to practical solutions General concepts: - basis sets - integrals and grids, electrostatics,
More informationContents. Preface to the Third Edition (2007) Preface to the Second Edition (1992) Preface to the First Edition (1985) License and Legal Information
Contents Preface to the Third Edition (2007) Preface to the Second Edition (1992) Preface to the First Edition (1985) License and Legal Information xi xiv xvii xix 1 Preliminaries 1 1.0 Introduction.............................
More informationPartial Differential Equations and the Finite Element Method
Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION
More informationSpectral Algorithms I. Slides based on Spectral Mesh Processing Siggraph 2010 course
Spectral Algorithms I Slides based on Spectral Mesh Processing Siggraph 2010 course Why Spectral? A different way to look at functions on a domain Why Spectral? Better representations lead to simpler solutions
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 informationLecture 16: DFT for Metallic Systems
The Nuts and Bolts of First-Principles Simulation Lecture 16: DFT for Metallic Systems Durham, 6th- 13th December 2001 CASTEP Developers Group with support from the ESF ψ k Network Overview of talk What
More informationNumerical Integra/on
Numerical Integra/on Applica/ons The Trapezoidal Rule is a technique to approximate the definite integral where For 1 st order: f(a) f(b) a b Error Es/mate of Trapezoidal Rule Truncation error: From Newton-Gregory
More informationDept of Mechanical Engineering MIT Nanoengineering group
1 Dept of Mechanical Engineering MIT Nanoengineering group » Recap of HK theorems and KS equations» The physical meaning of the XC energy» Solution of a one-particle Schroedinger equation» Pseudo Potentials»
More information2. Surface geometric and electronic structure: a primer
2. Surface geometric and electronic structure: a primer 2.1 Surface crystallography 2.1.1. Crystal structures - A crystal structure is made up of two basic elements: lattice + basis Basis: Lattice: simplest
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 informationTable of Contents. Foreword... Introduction...
Table of Contents Foreword.... Introduction.... xi xiii Chapter 1. Fundamentals of Heat Transfer... 1 1.1. Introduction... 1 1.2. A review of the principal modes of heat transfer... 1 1.2.1. Diffusion...
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 5 The Hartree-Fock method C.-K. Skylaris Learning outcomes Be able to use the variational principle in quantum calculations Be able to construct Fock operators
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 informationAb initio study of Mn doped BN nanosheets Tudor Luca Mitran
Ab initio study of Mn doped BN nanosheets Tudor Luca Mitran MDEO Research Center University of Bucharest, Faculty of Physics, Bucharest-Magurele, Romania Oldenburg 20.04.2012 Table of contents 1. Density
More informationTopics for the Qualifying Examination
Topics for the Qualifying Examination Quantum Mechanics I and II 1. Quantum kinematics and dynamics 1.1 Postulates of Quantum Mechanics. 1.2 Configuration space vs. Hilbert space, wave function vs. state
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationIntroduction. Chapter One
Chapter One Introduction The aim of this book is to describe and explain the beautiful mathematical relationships between matrices, moments, orthogonal polynomials, quadrature rules and the Lanczos and
More informationSpectral element agglomerate AMGe
Spectral element agglomerate AMGe T. Chartier 1, R. Falgout 2, V. E. Henson 2, J. E. Jones 4, T. A. Manteuffel 3, S. F. McCormick 3, J. W. Ruge 3, and P. S. Vassilevski 2 1 Department of Mathematics, Davidson
More informationBasics of density-functional theory and fast guide to actual calculations Matthias Scheffler
Basics of density-functional theory and fast guide to actual calculations Matthias Scheffler http://www.fhi-berlin.mpg.de/th/th.html I. From the many-particle problem to the Kohn-Sham functional II. From
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 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 informationCHEM6085: Density Functional Theory Lecture 10
CHEM6085: Density Functional Theory Lecture 10 1) Spin-polarised calculations 2) Geometry optimisation C.-K. Skylaris 1 Unpaired electrons So far we have developed Kohn-Sham DFT for the case of paired
More informationEigenvalues and Eigenvectors. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB
Eigenvalues and Eigenvectors Consider the equation A x = λ x, where A is an nxn matrix. We call x (must be non-zero) an eigenvector of A if this equation can be solved for some value of λ. We call λ an
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 informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationThe Gutzwiller Density Functional Theory
The Gutzwiller Density Functional Theory Jörg Bünemann, BTU Cottbus I) Introduction 1. Model for an H 2 -molecule 2. Transition metals and their compounds II) Gutzwiller variational theory 1. Gutzwiller
More informationPractical Guide to Density Functional Theory (DFT)
Practical Guide to Density Functional Theory (DFT) Brad Malone, Sadas Shankar Quick recap of where we left off last time BD Malone, S Shankar Therefore there is a direct one-to-one correspondence between
More informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
More informationNumerical integration and differentiation. Unit IV. Numerical Integration and Differentiation. Plan of attack. Numerical integration.
Unit IV Numerical Integration and Differentiation Numerical integration and differentiation quadrature classical formulas for equally spaced nodes improper integrals Gaussian quadrature and orthogonal
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationPreconditioned Eigenvalue Solvers for electronic structure calculations. Andrew V. Knyazev. Householder Symposium XVI May 26, 2005
1 Preconditioned Eigenvalue Solvers for electronic structure calculations Andrew V. Knyazev Department of Mathematics and Center for Computational Mathematics University of Colorado at Denver Householder
More informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
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 informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationRecent developments of Discontinuous Galerkin Density functional theory
Lin Lin Discontinuous Galerkin DFT 1 Recent developments of Discontinuous Galerkin Density functional theory Lin Lin Department of Mathematics, UC Berkeley; Computational Research Division, LBNL Numerical
More informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationMethod in 1D. Introduction of the Basic Concepts and Comparison to Optimal FD Operators
The Spectral Element Method in 1D Introduction of the Basic Concepts and Comparison to Optimal FD Operators Seismology and Seismics Seminar 17.06.2003 Ludwigs-Maximilians-Universität München by Bernhard
More informationVariational Methods for Electronic Structure
Variational Methods for Electronic Structure The hydrogen atom is a two-body system consisting of a proton and an electron. If spin and relativistic effects are ignored, then the Schrödinger equation for
More informationAtomic orbitals of finite range as basis sets. Javier Junquera
Atomic orbitals of finite range as basis sets Javier Junquera Most important reference followed in this lecture in previous chapters: the many body problem reduced to a problem of independent particles
More informationComputer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
More informationElectronic structure calculations with GPAW. Jussi Enkovaara CSC IT Center for Science, Finland
Electronic structure calculations with GPAW Jussi Enkovaara CSC IT Center for Science, Finland Basics of density-functional theory Density-functional theory Many-body Schrödinger equation Can be solved
More informationNon-linear least squares
Non-linear least squares Concept of non-linear least squares We have extensively studied linear least squares or linear regression. We see that there is a unique regression line that can be determined
More informationA Non-Linear Eigensolver-Based Alternative to Traditional Self-Consistent Electronic Structure Calculation Methods
University of Massachusetts Amherst ScholarWorks@UMass Amherst Masters Theses 1911 - February 2014 2013 A Non-Linear Eigensolver-Based Alternative to Traditional Self-Consistent Electronic Structure Calculation
More informationPRECONDITIONING ORBITAL MINIMIZATION METHOD FOR PLANEWAVE DISCRETIZATION 1. INTRODUCTION
PRECONDITIONING ORBITAL MINIMIZATION METHOD FOR PLANEWAVE DISCRETIZATION JIANFENG LU AND HAIZHAO YANG ABSTRACT. We present an efficient preconditioner for the orbital minimization method when the Hamiltonian
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 informationA MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC
More informationJournal of Computational Physics
Journal of Computational Physics 23 (211) 8573 862 Contents lists available at SciVerse ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp An eigen-based high-order
More informationEfficient hp-finite elements
Efficient hp-finite elements Ammon Washburn July 31, 2015 Abstract Ways to make an hp-finite element method efficient are presented. Standard FEMs and hp-fems with their various strengths and weaknesses
More information