Wavelets for density functional calculations: Four families and three. applications
|
|
- Brian James
- 5 years ago
- Views:
Transcription
1 Wavelets for density functional calculations: Four families and three Haar wavelets Daubechies wavelets: BigDFT code applications Stefan Goedecker Interpolating wavelets: Poisson solver and direct exchange library Multi-wavelets: High accuracy atomization energies
2 Wavelets: A family of relatively new mathematical basis sets with astonishing properties All families share the properties: of being localized both in real and in Fourier space. This allows for the so-called multi-resolution analysis that gives information about the localization properties both in real and in Fourier space of a function (signal)
3 of being a systematic basis set, i.e the error is guaranteed to approach zero as the basis set tends to infinity an arbitrary high degree of adaptivity can be obtained.
4 Wavelet basis functions Each family is characterized by two functions. The mother scaling function ψ and the mother wavelet φ. A basis set is generated by translations and dilatation-s of these two functions j: Translation, Localization in real space k: Dilatation, Localization in Fourier space ψ k j(x) ψ(2 k x j) φ k j(x) φ(2 k x j) Both the wavelet and the scaling function at a certain resolution level can be written as a linear combination of scaling functions at a higher resolution level (refinement relations) φ(x)= m h j φ(2x j) j= m
5 Haar wavelet ψ and scaling function φ φ ψ Scaling function representation φ 4 0 x 1 f(x)= s 4 j φ 4 j(x) j
6 Haar wavelet basis set ψ φ ψ ψ ψ Wavelet representation: f(x)=s 0 1φ 0 1(x)+d 0 1ψ 0 1(x)+ 2 di 1 ψ 1 i(x)+ i=1 4 di 2 ψ 2 i(x)+ i=1 8 di 3 ψ 3 i(x) (1) i=1
7 Operators in a wavelet basis: The standard form c S 0 0 D D 1 A b S 0 0 D D 1 D 2 = * D 2 D 3 D 3 Coupling between different resolution levels Works fine for a limited number (two) of resolution levels. Complicated and numerically inefficient structure for a large number of resolution levels
8 Operators in a wavelet basis: The non-standard operator form S D S 0 D S D S D b * = S D S 0 D S D S D c A No coupling between different resolution levels For a small number of resolution levels larger prefactor than standard form Easy and numerically efficient structure if the transition region between the different resolution levels is large compared to the support length of the wavelet
9 Daubechies wavelets Ingrid Daubechies, 1988: ORTHONORMAL BASES OF COMPACTLY SUPPORTED WAVELETS Daubechies wavelet basis set combines advantages of Plane waves: Systematic, orthogonal basis set Localization in Fourier space allows for efficient preconditioning techniques. Gaussians: Localized in real space: well suited for molecules and other open structures, direct way to linear scaling Adaptivity
10 Daubechies wavelet of order DAUBECHIES-4 SCALING FUNCTION AND WAVELET scaling function wavelet Even though the function is not very smooth polynomials up to degree 4 can be represented exactly
11 High order Daubechies wavelets are used to represent wave functions in BigDFT LEAST ASYMMETRIC DAUBECHIES-16 scaling function wavelet Polynomial up to degree 16 can be represented exactly: convergence rate h 14
12 BigDFT solves the many-electron Schrödinger equation in the Kohn-Sham density functional approximation using a Daubechies wavelet basis All boundary conditions All XC functionals of the LibXC library (LDA, GGA, hybrid functionals) Local geometry optimizations (with constraints) Saddle point searches Global geometry optimization: Minima Hopping Implicit solvation models (dielectric, electrolyte) linear scaling version Vibrations Born Oppenheimer MD collinear and non-collinear magnetism Empirical van der Waals terms High degree of parallelization and GPU acceleration: fast time to solution Website: inac.cea.fr/l Sim/BigDFT/
13 Wavelet basis sets in three dimensions 1 scaling function 7 wavelets all are products of 1-dim scaling functions and wavelets φ i, j,k (x,y,z) = φ(x i)φ(y j)φ(z k) ψ 1 i, j,k (x,y,z) = φ(x i)φ(y j)ψ(z k) ψ 2 i, j,k (x,y,z) = φ(x i)ψ(y j)φ(z k) ψ 3 i, j,k (x,y,z) = φ(x i)ψ(y j)ψ(z k) ψ 4 i, j,k (x,y,z) = ψ(x i)φ(y j)φ(z k) ψ 5 i, j,k (x,y,z) = ψ(x i)φ(y j)ψ(z k) ψ 6 i, j,k (x,y,z) = ψ(x i)ψ(y j)φ(z k) ψ 7 i, j,k (x,y,z) = ψ(x i)ψ(y j)ψ(z k)
14 Pseudopotentials: two resolution levels are used in BigDFT
15 Kinetic energy matrix elements between scaling functions a i = Can be calculated exactly (Beylkin): a i = φ(x) 2 x 2 φ(x i)dx = 2h ν h µ ν,µ = 2h ν h µ ν,µ = h ν h µ 2 2 ν,µ = h ν h µ 2 2 a 2i ν+µ ν,µ φ(x) 2 x 2 φ(x i)dx φ(2x ν) 2 x 2 φ(2x 2i µ)dx φ(y ν) 2 y 2 φ(y 2i µ)dy φ(y) 2 y l φ(y 2i µ+ν)dy We thus have to find the eigenvector a associated with the eigenvalue of 2 2, ( ) 1 2 A i, j a j = a i j 2
16 where the matrix A i, j is given by A i, j = h ν h µ δ j,2i ν+µ ν,µ Operator approach: Applying the kinetic energy operator results in a convolution with the filter a (KΨ) i = Ψ j a i j j In the 3-dim case applying the kinetic energy operator requires a 3-dim convolution with a filter that is a product of 3 1-dim filters. The operation count is 3 N 1 N 2 N 3 L where L is the length of the filter.
17 Interpolating scaling functions: PSolver library of BigDFT for solution of Poisson s equation Compact support Many continuous derivatives for a given support length Trivial transformation from a real space data set to a scaling function expansion Not orthogonal
18 Construction of interpolating scaling functions Recursive interpolation from Kronecker data set: Example linear interpolation scf
19 High order interpolating scaling functions represent charge densities and potentials 1 Interpolating scf 14 Interpolating scf
20 Solution of Poisson s equation for free boundary conditions L. Genovese, T. Deutsch, A. Neelov, S. Goedecker, G. Beylkin, J. Chem. Phys. 125 (2006) Given the values of the charge density on a regular grid, ρ i, j,k, the continuous charge distribution is represented in terms of interpolating scaling functions ρ(r)= ρ i, j,k φ(x i)φ(y j)φ(z k) i, j,k The moments of the discrete and continuous charge distributions ρ i, j,k and ρ(r) are identical i l 1 j l 2 k l 3 ρ i, j,k = i, j,k dr x l 1 y l 2 z l 3 ρ(r) (2) if l 1,l 2,l 3 < m, where m is the order of the scaling functions. The potential at a grid point i 1,i 2,i 3 is given by V i1,i 2,i 3 = φ(x j 1 )φ(y j 2 )φ(z j 3 ) ρ j1, j 2, j 3 j 1, j 2, j 3 r j1, j 2, j 3 r dr = ρ j1, j 2, j 3 K i1 j 1,i 2 j 2,i 3 j 3 j 1, j 2, j 3 The above convolution can be calculated rapidly with Fourier methods
21 Boundary conditions The Poisson equation is solved exactly for all boundary conditions free wire surface periodic Highly accurate treatment of charged clusters dipolar surfaces clusters, surfaces in electric fields
22 Efficient and accurate hybrid functional calculations with the PSolver package of BigDFT The exact exchange energy E x is given by E X = N i=1 N j=i+1 dr dr ψ i (r) ψ j(r) ψ i (r ) ψ j (r ) r r (3) For a system of N electrons it requires the solution of N(N 1)/2 Poisson equations for all pairwise charge densities ψ i (r ) ψ j (r At greatly reduced cost hybrid functional calculations should become much more widespread in systematic basis sets since: High accuracy can nearly be reached in atomization energies Materials that are problematic with other functionals such as transition metal oxides can be well treated Gaps in solids are reasonably accurate
23 A direct hybrid functional calculation is about two orders of magnitude more expensive than a GGA calculation with all the state of the art plane wave codes BigDFT gives extremely high speed in the expensive exact exchange energy part: Since only a few basic operations have to be performed in our scaling function basis, all the computations are implemented in CUDA and are executed on the GPU Communication is overlapped with computation Direct GPU communication
24 Hybrid functional calculations possible up to 1000 atoms: Cray (Piz Daint) at CSCS
25 Exact ionic forces In contrast to other approaches that gain speed by evaluating the exact exchange based on localized orbitals, no cutoffs or other approximations are necessary and the forces are the exact derivative of the energy. This leads for instance to a perfect energy conservation in MD.
26 PBE0 pseudopotentials are available and should be used in PBE0 calculations XC used in calculation PBE PBE0 PBE0 XC of pseudopotential used PBE PBE PBE0 C 2 H s CH 3 Cl s CH 3 OH s CH s CO s H 2 CO s H 2 O s HOCl s OH d Atomization energy of molecules in kcal/mol for consistent (PBE/PBE and PBE0/PBE0) and inconsistent (PBE/PBE0) use of the exchange correlation functionals in the molecular calculation and for the generation of the pseudopotential. Dual space Gaussian pseudopotentials were used in the BigDFT code with free boundary conditions.
27 Overview: Poisson Solver package of BigDFT Solves the Poisson equation both for a constant and a spatially varying dielectric constant as well as the Poisson Boltzmann equation to describe electrolytes all standard boundary conditions Input and output are given on equally spaced Cartesian grids. In the periodic case both orthorombic and non-orthorombic cells can be treated The hybrid functional package is based on the Poisson solver package and the LibXC library
28 compact support Multi-wavelets: All electron calculations several basis functions per support interval representable as polynomials orthogonal symmetric Continuous derivatives within interval, possibly discontinuities among neighboring intervals (discontinuous Galerkin) Only integral equations can be solved
29 Solving Schrödingers equation in integral form Ψ i (r)+v(r)ψ i (r)=ε i Ψ i (r) ( 2 + 2ε i ) Ψi (r)=4π 1 2π V(r)Ψ i(r) Helmholtz equation: The inverse operator of 2 + 2ε i is exp( 2εi r r ) r r dr. Hence we obtain the following iteration scheme for the Kohn-Sham orbitals Ψ new i (r)= exp( 2εi r r ) r r Ψ old i (r )dr
30 MRChem Frediani et al.: Real-space numerical grid methods in quantum chemistry, Phys. Chem. Chem. Phys. 2015, 17, Developed in the group of Luca Frediani in Norway by Stig Jensen et al Performs non-relativistic all-electron density functional calculations for LDA, GGA and hybrid functionals by using multi-wavelets Allows for an arbitrary number of resolution levels and uses the non-standard operator form Any preset accuracy for the wave functions/energies can be obtained (if enough memory/cores are available) Same underlying method as in MADNESS, highly stable implementation Program under development to include more features
31 Accuracy of DFT calculations Nearly all codes use approximations that go beyond the XC functional: basis sets, pseudopotentials Kurt Lejaeghere et al.: Reproducibility in density functional theory calculations of solids, Science 351, 6280 (2016) A large number of electronic structure codes give more or less identical results for the energy versus volume curve For more difficult quantities significant disagreement between different codes can still be found Atomization energies with µha accuracy obtained with MRChem for a test set of nearly 300 molecules: S. Jensen et al.: The Elephant in the Room of Density Functional Theory Calculations. Phys. Chem. Lett., 2017, 8 (7), pp For codes with a systematic basis set the accuracy is only limited by the pseudopotential or PAW scheme.
32 Basis set errors
33 Pseudopotential errors
34 People involved in this work Basel: Bastian Schaefer, Augustin Degomme, Giuseppe Fisicaro, Santanu Saha, Huan Tran, Alireza Ghasemi Grenoble: Luigi Genovese, Damien Caliste, Thierry Deutsch Arctic University of Norway: Stig Jensen, Luca Frediani
35 Selected Applications of BigDFT Geometric ground state of metal decorated boron clusters Santanu Saha et al., accepted by Nature Scientific reports
36 Selected Applications of BigDFT Disconnectivity graph for a 26 atom gold cluster Au + 26 Bastian Schaefer et al. ACS Nano, (2014)
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 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 informationPoisson Solver, Pseudopotentials, Atomic Forces in the BigDFT code
CECAM Tutorial on Wavelets in DFT, CECAM - LYON,, in the BigDFT code Kernel Luigi Genovese L_Sim - CEA Grenoble 28 November 2007 Outline, Kernel 1 The with Interpolating Scaling Functions in DFT for Interpolating
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 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 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 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 informationA new future for finite-element methods in Quantum Chemistry?
A new future for finite-element methods in Quantum Chemistry? Kenneth Ruud Department of Chemistry University of Tromsø 9037 Tromsø Norway Geilo, 29/1 2007 A new future for finite-element methods in Quantum
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 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 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 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 informationIntroduction to Density Functional Theory
Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary
More informationWalter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.
Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm
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 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 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 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 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 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 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 informationOptimized energy landscape exploration for nanosciences using ab-initio based methods
Optimized energy landscape exploration for nanosciences using ab-initio based methods Eduardo Machado-Charry Nanosciences Foundation & Laboratoire de simulation atomistique (L Sim), SP2M, UMR-E CEA-Grenoble
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 informationModified Becke-Johnson (mbj) exchange potential
Modified Becke-Johnson (mbj) exchange potential Hideyuki Jippo Fujitsu Laboratories LTD. 2015.12.21-22 OpenMX developer s meeting @ Kobe Overview: mbj potential The semilocal exchange potential adding
More informationDFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory
DFT and beyond: Hands-on Tutorial Workshop 2011 Tutorial 1: Basics of Electronic Structure Theory V. Atalla, O. T. Hofmann, S. V. Levchenko Theory Department, Fritz-Haber-Institut der MPG Berlin July 13,
More informationPost Hartree-Fock: MP2 and RPA in CP2K
Post Hartree-Fock: MP2 and RPA in CP2K A tutorial Jan Wilhelm jan.wilhelm@chem.uzh.ch 4 September 2015 Further reading MP2 and RPA by Mauro Del Ben, Jürg Hutter, Joost VandeVondele Del Ben, M; Hutter,
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 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 informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
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 informationCP2K: Selected Developments
CP2K: Selected Developments Jürg Hutter Department of Chemistry University of Zurich Outline Introduction History and Performance Current and Future Developments Post-Hartree-Fock Methods GW Methods RPA
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 informationPractical calculations using first-principles QM Convergence, convergence, convergence
Practical calculations using first-principles QM Convergence, convergence, convergence Keith Refson STFC Rutherford Appleton Laboratory September 18, 2007 Results of First-Principles Simulations..........................................................
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 informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationInstitut Néel Institut Laue Langevin. Introduction to electronic structure calculations
Institut Néel Institut Laue Langevin Introduction to electronic structure calculations 1 Institut Néel - 25 rue des Martyrs - Grenoble - France 2 Institut Laue Langevin - 71 avenue des Martyrs - Grenoble
More informationAb initio molecular dynamics
Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field
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 informationIntroduction to Vibrational Spectroscopy
Introduction to Vibrational Spectroscopy Harmonic oscillators The classical harmonic oscillator The uantum mechanical harmonic oscillator Harmonic approximations in molecular vibrations Vibrational spectroscopy
More informationComputational Chemistry I
Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationPseudopotentials: design, testing, typical errors
Pseudopotentials: design, testing, typical errors Kevin F. Garrity Part 1 National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015 Parameter free calculations.
More informationDensity Functional Theory - II part
Density Functional Theory - II part antonino.polimeno@unipd.it Overview From theory to practice Implementation Functionals Local functionals Gradient Others From theory to practice From now on, if not
More informationWe have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma and (3.55) with j = 1. We can write any f W 1 as
88 CHAPTER 3. WAVELETS AND APPLICATIONS We have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma 3..7 and (3.55) with j =. We can write any f W as (3.58) f(ξ) = p(2ξ)ν(2ξ)
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 informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle
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 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 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 informationABSTRACT. Title: Wavelet-space solution of the Poisson equation: An algorithm for use in particle-incell simulations
ABSTRACT Name: Benjamin Sprague Department: Physics Title: Wavelet-space solution of the Poisson equation: An algorithm for use in particle-incell simulations Major: Physics Degree: Master of Science Approved
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 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 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 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 informationQ-Chem 4.0: Expanding the Frontiers. Jing Kong Q-Chem Inc. Pittsburgh, PA
Q-Chem 4.0: Expanding the Frontiers Jing Kong Q-Chem Inc. Pittsburgh, PA Q-Chem: Profile Q-Chem is a high performance quantum chemistry program; Contributed by best quantum chemists from 40 universities
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 informationAn Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method
Fakultät für Mathematik und Naturwissenschaften - Lehrstuhl für Physikalische Chemie I / Theoretische Chemie An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method Jan-Ole Joswig
More informationCLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122
CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 191 THANKS TO MANY COLLABORATORS, INCLUDING SY VOSKO DAVID LANGRETH ALEX
More informationIntroduction to first-principles modelling and CASTEP
to first-principles modelling and Phil Hasnip to + Atomistic Simulations If we know what the bonding in a material is beforehand, then we can often find good expressions for the forces between atoms, e.g.
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 informationDay 1 : Introduction to AIMD and PIMD
Day 1 : Introduction to AIMD and PIMD Aug 2018 In today s exercise, we perform ab-initio molecular dynamics (AIMD) and path integral molecular dynamics (PIMD) using CP2K[?]. We will use the Zundel s cation
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 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 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 information1 Introduction to Wavelet Analysis
Jim Lambers ENERGY 281 Spring Quarter 2007-08 Lecture 9 Notes 1 Introduction to Wavelet Analysis Wavelets were developed in the 80 s and 90 s as an alternative to Fourier analysis of signals. Some of the
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Chemistry 3502/4502 Final Exam Part I May 14, 2005 1. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle (e) The
More informationDENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY
DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES
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 informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
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 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 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 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 informationTransform methods. and its inverse can be used to analyze certain time-dependent PDEs. f(x) sin(sxπ/(n + 1))
AMSC/CMSC 661 Scientific Computing II Spring 2010 Transforms and Wavelets Dianne P. O Leary c 2005,2010 Some motivations: Transform methods The Fourier transform Fv(ξ) = ˆv(ξ) = v(x)e ix ξ dx, R d and
More informationAb-initio modeling of opto-electronic properties of molecules in solvents and in proximity to a semiconductor nanoparticle
Ab-initio modeling of opto-electronic properties of molecules in solvents and in proximity to a semiconductor nanoparticle Alain Delgado (a,b), Stefano Corni (b), Carlo Andrea Rozzi (b) Stefano Pittalis
More informationVASP: running on HPC resources. University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria
VASP: running on HPC resources University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria The Many-Body Schrödinger equation 0 @ 1 2 X i i + X i Ĥ (r 1,...,r
More informationBefore we start: Important setup of your Computer
Before we start: Important setup of your Computer change directory: cd /afs/ictp/public/shared/smr2475./setup-config.sh logout login again 1 st Tutorial: The Basics of DFT Lydia Nemec and Oliver T. Hofmann
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 informationLocalized Optimization: Exploiting non-orthogonality to efficiently minimize the Kohn-Sham Energy
Localized Optimization: Exploiting non-orthogonality to efficiently minimize the Kohn-Sham Energy Courant Institute of Mathematical Sciences, NYU Lawrence Berkeley National Lab 16 December 2009 Joint work
More informationOrbital dependent correlation potentials in ab initio density functional theory
Orbital dependent correlation potentials in ab initio density functional theory noniterative - one step - calculations Ireneusz Grabowski Institute of Physics Nicolaus Copernicus University Toruń, Poland
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 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 informationWavelets For Computer Graphics
{f g} := f(x) g(x) dx A collection of linearly independent functions Ψ j spanning W j are called wavelets. i J(x) := 6 x +2 x + x + x Ψ j (x) := Ψ j (2 j x i) i =,..., 2 j Res. Avge. Detail Coef 4 [9 7
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 informationComparison of various abinitio codes used in periodic calculations
Comparison of various abinitio codes used in periodic calculations 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India & Center for Materials Science and Nanotechnology,
More informationIntermolecular Forces in Density Functional Theory
Intermolecular Forces in Density Functional Theory Problems of DFT Peter Pulay at WATOC2005: There are 3 problems with DFT 1. Accuracy does not converge 2. Spin states of open shell systems often incorrect
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 informationSCF calculation on HeH +
SCF calculation on HeH + Markus Meuwly Department of Chemistry, University of Basel, Basel, Switzerland Abstract This document describes the main steps involved in carrying out a SCF calculation on the
More informationBoundary functions for wavelets and their properties
Graduate Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 009 Boundary functions for wavelets and their properties Ahmet Alturk Iowa State University Follow this and additional
More informationMGA Tutorial, September 08, 2004 Construction of Wavelets. Jan-Olov Strömberg
MGA Tutorial, September 08, 2004 Construction of Wavelets Jan-Olov Strömberg Department of Mathematics Royal Institute of Technology (KTH) Stockholm, Sweden Department of Numerical Analysis and Computer
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 informationLarge-scale real-space electronic structure calculations
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: FA9550-13-1-0113
More informationAl 2 O 3. Al Barrier Layer. Pores
Overview of DFT Calculations for Nanoscience and Nanotechnology Oğuz Gülseren Al 2 O 3 Al Barrier Layer Pores Jules Verne A Trip to the Moon, 1865 Nanostructures: a. Contain a countable number of atoms
More informationDensity Functional Theory
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function
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 informationCar-Parrinello Molecular Dynamics
Car-Parrinello Molecular Dynamics Eric J. Bylaska HPCC Group Moral: A man dreams of a miracle and wakes up with loaves of bread Erich Maria Remarque Molecular Dynamics Loop (1) Compute Forces on atoms,
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Yesterday you (should have) learned this: http://upload.wikimedia.org/wikipedia/commons/e/ea/ Simple_Harmonic_Motion_Orbit.gif 1. deterministic
More informationIs there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich
Is there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich Chemistry Chemistry is the science of atomic matter, especially its chemical reactions,
More informationIntroduction to DFT and its Application to Defects in Semiconductors
Introduction to DFT and its Application to Defects in Semiconductors Noa Marom Physics and Engineering Physics Tulane University New Orleans The Future: Computer-Aided Materials Design Can access the space
More information