Compression Algorithms for Electronic Structure Computations
|
|
- Jade Barker
- 5 years ago
- Views:
Transcription
1 Compression Algorithms for Electronic Structure Computations François Gygi University of California, Davis IPAM Workshop, Computation meets Big Data, Feb 2, 2017 Supported by DOE BES DE-SC and MICCoM 1
2 Context: First-Principles Molecular Dynamics Monte Carlo Molecular Dynamics Statistical Mechanics FPMD Electronic Structure Theory Quantum Chemistry Density Functional Theory 2
3 Context: First-Principles Molecular Dynamics Molecular Dynamics Density Functional Theory 2 d m R = dt F i 2 i i FPMD ( ) ρ Δ + V ϕ ( x) = εϕ ( x) ( x) eff n = ϕi ( x) i= 1 i i i 2 Newton equations Kohn-Sham equations R. Car and M. Parrinello (1985) 3
4 Electronic Structure and dynamical properties of complex structures Complex structures Nanoparticles Assemblies of nanoparticles Embedded nanoparticles Liquid/solid interfaces 4
5 Nanoparticles Multiplicity of locally stable structures requires extensive sampling Electronic structure requires accurate methods (beyond DFT) Finite temperature properties require first-principles molecular dynamics Cd 34 Se 34 5
6 Liquids and Liquid-Solid Interfaces Embedded nanoparticles, assemblies of nanoparticles Annealing of structures requires MD simulations Calculation of band gaps and band alignments requires accurate electronic structure (beyond DFT) Si/ZnS S. Wippermann, M. Vörös, A. Gali, F. Gygi, G. Zimanyi, G.Galli, Phys. Rev. Lett. 112, (2014).
7 Liquids and Liquid-Solid Interfaces H 2 O/Si(100)H Liquids require finite temperature simulations ab initio MD multiple samples/replica exchange simulations Electronic structure band alignment Spectroscopy requires calculation of IR and Raman spectra 7
8 Challenges Multiple length and time scales e.g. nanoparticle assembly process, non-equilibrium processes Finite temperature: MD/MC simulations ab initio MD for large systems (>1000 atoms) Need for accurate electronic structure hybrid DFT and/or GW/BSE level 8
9 Qbox code: DFT and hybrid DFT first-principles molecular dynamics massively parallel first-principles MD C++/MPI/OpenMP DFT and hybrid DFT MD GPL license All algorithms discussed here are available in the Qbox code 9
10 Computation meets Big Data: Why data compression is necessary Storage of restart files in simulations contain full information about wave functions size is O(N 2 ) for N atoms, reaches 100 GB-1 TB for large problems Acceleration of DFT and/or hybrid DFT calculations Reduce cost (e.g. from O(N 4 ) to O(N 3 )) Exploit locality of data on modern computer architectures 10
11 Desirable properties of compression schemes Accuracy control ideally a single parameter controlling the accuracy gradual reduction of the error to zero Efficiency tradeoff between space efficiency and cost of compression algorithm 11
12 Three compression approaches Reduced resolution Easy: in the Fourier basis: Reduce energy cutoff Efficient: requires only 1 FFT per orbital Problem: only moderate reduction is possible Compute Wannier functions Efficient (if using the right algorithm) Truncation procedure is ill-defined (no error control) Recursive Subspace Bisection (this work) Efficient (as fast as Wannier function calculation) Controllable error and systematic truncation scheme 12
13 Wannier functions Existence of exponentially localized Wannier functions Kohn (1959), des Cloiseaux (1964), Nenciu (1983), Helffer et al (1989), Brouder et al (2007), Panati (2007, 2013) Computation of Wannier functions: find an orthogonal transformation among orbitals that minimizes the spread σ 2 X = x x ( ) 2 Optimization problem (with local minima) Marzari, Vanderbilt (1997) Use conjugate gradients, etc. Approximate simultaneous diagonalization problem F.G, Fattebert, Schwegler (2003) F.G., J.L.Fattebert, E.Schwegler, Comput.. Phys. Comm. 155, 1 (2003) J.F.Cardoso and A. Souloumiac, SIAM J. Mat. Anal. Appl. 17, 161 (1996). 13
14 Spread Functionals Spread of an operator  (single orbital) σ 2 (  φ) = φ (  φ  φ ) 2 φ = φ Â2 φ φ  φ 2 Spread of a set of orbitals ({ }) ( ) ˆ σ φ = σ φi 2 2 Aˆ i A i 14
15 Spread Functionals The spread is not invariant under orthogonal transformations among orbitals ψ i = j x ij φ j X! n n orthogonal 2 σ Â ({ ψ }) 2 i σ Â ({ φ }) i There exists a matrix X that minimizes the spread 15
16 Spread Functionals Let A,B! n n a ij = i  j b ij = i Â2 j 2 σ  ({ ψ }) i = tr X T BX Minimize the spread = maximize = diagonalize A n i=1 ( ) ( X T AX ) 2 ii n i=1 ( X T AX ) 2 ii 16
17 Spread Functionals Case of multiple operators operators matrices 2 σ Â Â (k ) k =1,,m A (k ) k =1,,m ({ ψ }) 2 i = σ k ) Â( i k ( ψ ) i Minimize the spread = maximize n i=1 k = joint approximate diagonalization of the matrices A (k) ( X T A (k ) X ) 2 ii 17
18 Spread Functionals Example of multiple operators ( ϕ) ˆ (1) = ˆ ˆ (,, ) ϕ(,, ) A X X x y z x x y z ( ϕ) ˆ (2) = ˆ ˆ (,, ) ϕ(,, ) A Y Y x y z y x y z ( ϕ) ˆ (3) = ˆ ˆ (,, ) ϕ(,, ) A Z Z x y z z x y z The matrices A (k) do not necessarily commute, even if the operators  (k) do commute 18
19 Calculation of Wannier functions In periodic systems  (1) =  (3) = Ĉ y  (5) = Ĉx cos 2π L x ˆx  (2) = cos 2π L y ŷ  (4) = Ŝ y Ĉz cos 2π L z ẑ  (6) = Ŝx sin 2π L x sin 2π L y ˆx 2π Ŝz sin ẑ L z ŷ 19
20 Calculation of Wannier functions The spread is minimized by simultaneous diagonalization of the matrices C x,s x,c y,s y,c z,s z Positions of the center of mass of the localized solutions ( Wannier centers ) Spreads x θi Lx 2π y θi x τi = Ly θi arctan 2π = z θi Lz 2π 2 2 ( σ i ) = Lx 1 ( Cx) ( Sx) x ( 2 2 ) ii ii ( Sx ) ( C ) x ii ii F.G., J.L.Fattebert, E.Schwegler, Comput.. Phys. Comm. 155, 1 (2003) J.F.Cardoso and A. Souloumiac, SIAM J. Mat. Anal. Appl. 17, 161 (1996). 20
21 Maximally localized Wannier functions Extended orbital Wannier function 21
22 Truncation of Wannier functions Wannier functions (WFs) can be truncated in real space truncate orbital to zero below a given threshold truncate orbital to zero outside of a given radius WFs can have variable localization properties 22
23 Maximally localized Wannier functions Extended orbital Wannier function 23
24 Truncation of Wannier functions Wannier functions (WFs) can be truncated in real space truncate orbital to zero below a given threshold truncate orbital to zero outside of a given radius WFs can have variable localization properties Are there other ways to localize orbitals? 24
25 Recursive subspace bisection Localize some orbitals on domains of decreasing size divide the simulation domain into two subdomains Ω localize orbitals on 1) 2) or 3) keep extended on Ω L Ω L ΩR R apply recursively to smaller domains V ΩL ΩR ΩL ΩR F.G. Phys. Rev. Lett. 102, (2009) 25
26 Subspace Bisection Y = [φ 1 φ n ] ΩL ΩR Y YV 26
27 The CS decomposition A matrix Y having orthogonal columns, can be T decomposed as Y1 U1Σ1V Y = T Y = 2 U2Σ2V where U 1, U 2, V are orthogonal matrices,! Σ 1 = # C " 0 $! & Σ 2 = # % " C = diag(c 1,,c n ) S = diag(s 1,,s n ) c i 2 + s i 2 =1 S 0 $ & % Stewart (1982) 27
28 CS Decomposition ΩL ΩR Y YV 28
29 CS Decomposition norm = c 2 i ΩL ΩR Y norm = si = 1 c YV 2 2 i 29
30 CS Decomposition ΩL localized in Ω L!" # $# 2 ci < ε ΩR Y 2 si < ε!" # $#!" # $# extended states localized in Ω R 30
31 Subspace Bisection Algorithm 1. Choose the acceptable 2-norm error 2. Perform a CS decomposition of the matrix Y 3. For each vector of YV: if c i 2 < ε orbital localized in Ω R (truncate in Ω L ) else if s i 2 < ε orbital localized in Ω L (truncate in Ω R ) else orbital is extended Ideal limit: 2-fold data reduction Cost: The CS decomposition can be achieved by diagonalization of the matrix Y 1T Y 1 ε 31
32 Recursive Subspace Bisection Bisection is applied simultaneously in 3 directions Implementation: simultaneous (approximate) diagonalization of symmetric matrices F.G, J.-L.Fattebert, E.Schwegler, Comp. Phys. Comm. 155, 1 (2003). J.F.Cardoso and A. Souloumiac, SIAM J. Mat. Anal. Appl. 17, 161 (1996). 32
33 Recursive subspace bisection Use multiple bisecting planes in each direction W j (x) W 1 (x) W 3 (x) W 6 (x) W 12 (x) Optimal placement of bisecting planes is provided by Walsh functions W j (x) j=1,3,6,12, x 33
34 Recursive subspace bisection Phys. Rev. Lett. 102, (2009). 34
35 (H 2 O) 512 orbitals after bisection localized extended 35
36 Compression ratio (H 2 O) 512 ε = 10 Data size reduction: N % N 1/ % N 1/ % N 1/ % (19x0) Carbon nanotube (304 atoms) ε = 2048 orbitals 10 3 Data size reduction: 2.72 N % N 1/ % N 1/ % N 1/ % 608 orbitals 36
37 Compression ratio using recursive bisection recursion on l max levels, l max =1,2,3 1 4x compression ratio 0.1 l max = 1 l max = 2 r = x 36x l max = N mol r = truncated MLWFs 1/(N logn) slope F.G., Phys. Rev. Lett. 102, (2009) 37
38 Localization of orbitals in inhomogeneous systems W. Dawson and F.G, JCTC 11, 4655 (2015) 38
39 Free electron gas: Distribution of CS singular values 2 c i Y matrix = F n p (section of the Fourier matrix, i.e. eigenvectors of the Laplacian) 2 c i Y matrix = random orthogonal matrix A. Edelman et al. SIAM J. Sci. Comput. 20, 1094, (1999) 39
40 Hybrid density functionals Conventional density functionals are often insufficient to describe weak bonds (e.g. hydrogen bonds) or optical properties (band gap) Hybrid density functionals include a fraction of the Hartree-Fock exchange energy The Hartree-Fock exchange energy involves N(N-1)/2 exchange integrals (for all e - pairs) E x = 1 2 N i, j ϕ ( i r 1 )ϕ ( i r 2 )ϕ ( j r 1 )ϕ ( j r ) 2 dr 1 dr 2 r 1 r 2 Cost: O(N 3 log N) (with large prefactor) for plane waves For atom-centered basis sets: O(N) (Strain, Scuseria (1996), Burant, Scuseria, Frisch (1996), Schwegler, Challacombe, Head- Gordon (1997)) 40
41 Acceleration of Hartree-Fock and hybrid DFT calculations N(N-1)/2 exchange integrals (all e - pairs) E x = 1 2 N i, j ϕ ( i r 1 )ϕ ( i r 2 )ϕ ( j r 1 )ϕ ( j r ) 2 dr r 1 r 1 dr 2 2 using bisection: non-overlapping pairs can be skipped in the sum The error is positive F.G., I. Duchemin, JCTC 9, (2013). 41
42 Speedup of hybrid-dft calculations vs truncation threshold (H 2 O) 32, (H 2 O) 64 (H 2 O) 256 truncation threshold (%) ε time (s) speedup 0.0 (ref) % % % % % %
43 Truncation error due to bisection (H 2 O) 63 Cl Energy error E time( )/time(0) E (a.u.) time( )/time(0)
44 Energy error per orbital vs threshold 44
45 Force error vs threshold 45
46 Error in ionic forces for MD applications "fhist-0.01.dat" Error caused by bisection with 0.01 threshold "fhist-0-100ry.dat" Error caused by changing Ecut from 85Ry to 100Ry Force error (a.u.) Recursive bisection affects forces in a controlled way 46
47 Error in ionic forces for MD applications "fhist dat" Error caused by bisection with threshold "fhist-0-100ry.dat" Error caused by changing Ecut from 85Ry to 100Ry Force error (a.u.) Recursive bisection affects forces in a controlled way 47
48 HOMO-LUMO gap, band gaps Hybrid DFTs lead to large improvements in band gaps (Henderson, Paier, Scuseria, PhysStatSol 2011) Bisection algorithm: occupied and empty orbitals must not be mixed when localizing orbitals (H 2 O) 63 Cl - A (k ) = 0 threshold Egap (ev)
49 Localization of empty orbitals W. Dawson and F.G, JCTC 11, 4655 (2015) 49
50 Hybrid-DFT electronic structure of Si nanoparticles embedded in ZnS LDA PBE0 Si 66 Zn 228 S 218 DOS Si 66 Zn 228 S electrons 200 empty orbitals PBE0 DOS (arb. units) 1 scf step, 5 iterations BG/Q, 16k cores, 772 s E (ev) S. Wippermann, M. Vörös, A. Gali, F. Gygi, G. Zimanyi, G.Galli, Phys. Rev. Lett. 112, (2014).
51 Hybrid-DFT electronic structure of bulk SiC (4096 atoms) 4096 atoms electrons PBE0 electronic structure (hybrid) recursive subspace bisection ANL Mira (BG/Q) 64k cores 357s/self-consistent iteration hybrid DFT electronic structure for 4096 atoms
52 Where is the (Big) Data? Making data accessible is critical for verification and validation of simulation software Agreeing on data formats has proved difficult.. XML schemas for electronic structure data repository of reference simulations 52
53 Summary Simulation of complex materials Truncation of Maximally Localized Wannier functions Recursive subspace bisection Controlled error in inhomogeneous systems Acceleration of hybrid DFT simulations Supported by DOE BES DE-SC and the MICCoM DOE center 53
54 Acknowledgements Giulia Galli (UChicago) Marco Govoni (UChicago) Stefan Wippermann (MPI) Eric Schwegler (LLNL) Funding DOE BES computer time DOE INCITE/ANL-ALCF NSF XSEDE NERSC William Dawson (CS, UCDavis) Martin Schlipf (CS, UCDavis) Cui Zhang (UCDavis) Alex Gaiduk (UChicago) Marton Vörös (UCDavis, ANL) Quan Wan (UChicago) T.-Anh Pham (LLNL) Supported by DOE BES DE-SC and the DOE MICCoM center
First-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 informationSearch for materials to harvest light
Solar nanocomposites with complementary charge extraction pathways for electrons and holes: Si embedded in ZnS S. Wippermann, M. Vörös, F. Gygi, A. Gali, G. Zimanyi, G. Galli NanoMatFutur DPG-2014, 04/03/2014
More informationSi-nanoparticles embedded in solid matrices for solar energy conversion: electronic and optical properties from first principles
Si-nanoparticles embedded in solid matrices for solar energy conversion: electronic and optical properties from first principles S. Wippermann, M. Vörös, D. Rocca, T. Li, A. Gali, G. Zimanyi, F. Gygi,
More informationPerformance optimization of WEST and Qbox on Intel Knights Landing
Performance optimization of WEST and Qbox on Intel Knights Landing Huihuo Zheng 1, Christopher Knight 1, Giulia Galli 1,2, Marco Govoni 1,2, and Francois Gygi 3 1 Argonne National Laboratory 2 University
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 informationLarge scale GW calculations
M. Govoni, UoC & ANL T.A. Pham, LLNL Nguyen H.-Viet, Vietnam Acad. of Sci.&Technology D. Rocca, U.Lorraine, France H. Wilson RMIT Univ., Australia I.Hamada, P. Scherpelz, NIMS, UoC Japan & UoC Large scale
More informationMultiple Exciton Generation in Si and Ge Nanoparticles with High Pressure Core Structures
Multiple Exciton Generation in Si and Ge Nanoparticles with High Pressure Core Structures S. Wippermann, M. Vörös, D. Rocca, A. Gali, G. Zimanyi, G. Galli NanoMatFutur DPG-214, 4/3/214 Multiple Exciton
More informationHigh pressure core structures of Si nanoparticles for solar energy conversion
High pressure core structures of Si nanoparticles for solar energy conversion S. Wippermann, M. Vörös, D. Rocca, A. Gali, G. Zimanyi, G. Galli [Phys. Rev. Lett. 11, 4684 (213)] NSF/Solar DMR-135468 NISE-project
More informationSemiconducting nano-composites for solar energy conversion: insights from ab initio calculations. S. Wippermann, G. Galli
Semiconducting nano-composites for solar energy conversion: insights from ab initio calculations S. Wippermann, G. Galli ICAMP-12, 08/10/2012 Search for materials to harvest light: Desperately seeking
More informationHierarchical Coupling of First Principles. Molecular Dynamics with Advanced Sampling. Methods: Supplementary Information
Hierarchical Coupling of First Principles Molecular Dynamics with Advanced Sampling Methods: Supplementary Information Emre Sevgen,, Federico Giberti,, Hythem Sidky, Jonathan K. Whitmer, Giulia Galli,
More informationSupporting Information for Interfacial Effects on. the Band Edges of Functionalized Si Surfaces in. Liquid Water
Supporting Information for Interfacial Effects on the Band Edges of Functionalized Si Surfaces in Liquid Water Tuan Anh Pham,,, Donghwa Lee, Eric Schwegler, and Giulia Galli, Department of Chemistry, University
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 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 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 informationExchange-Correlation Functional
Exchange-Correlation Functional Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
More informationTDDFT in Chemistry and Biochemistry III
TDDFT in Chemistry and Biochemistry III Dmitrij Rappoport Department of Chemistry and Chemical Biology Harvard University TDDFT Winter School Benasque, January 2010 Dmitrij Rappoport (Harvard U.) TDDFT
More informationMD simulation: output
Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square
More informationQuantum Monte Carlo Benchmarks Density Functionals: Si Defects
Quantum Monte Carlo Benchmarks Density Functionals: Si Defects K P Driver, W D Parker, R G Hennig, J W Wilkins (OSU) C J Umrigar (Cornell), R Martin, E Batista, B Uberuaga (LANL), J Heyd, G Scuseria (Rice)
More informationLinking electronic and molecular structure: Insight into aqueous chloride solvation. Supplementary Information
Linking electronic and molecular structure: Insight into aqueous chloride solvation Ling Ge, Leonardo Bernasconi, and Patricia Hunt Department of Chemistry, Imperial College London, London SW7 2AZ, United
More informationSelf-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT
Self-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT Kiril Tsemekhman (a), Eric Bylaska (b), Hannes Jonsson (a,c) (a) Department of Chemistry,
More informationPath Integral Monte Carlo Simulations on the Blue Waters System. Burkhard Militzer University of California, Berkeley
Path Integral Monte Carlo Simulations on the Blue Waters System Burkhard Militzer University of California, Berkeley http://militzer.berkeley.edu Outline 1. Path integral Monte Carlo simulation method
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 informationUnderstanding hydrogen storage in metal organic frameworks from first principles
Understanding hydrogen storage in metal organic frameworks from first principles Understanding hydrogen storage in metal organic frameworks from first principles Sohrab Ismail-Beigi (Yale Applied Physics)
More informationChemical Dynamics of the First Proton Coupled Electron Transfer of Water Oxidation on TiO 2 Anatase
Supplementary Information Chemical Dynamics of the First Proton Coupled Electron Transfer of Water Oxidation on TiO 2 Anatase Jia Chen, Ye-Fei Li, Patrick Sit, and Annabella Selloni Department of Chemistry,
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 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 informationExact embedding of local defects in crystals
Exact embedding of local defects in crystals Mathieu LEWIN Mathieu.Lewin@math.cnrs.fr (CNRS & University of Cergy-Pontoise) joint work with É. Cancès & A. Deleurence (ENPC, Marne-La-Vallée) Minneapolis,
More informationOrbital Density Dependent Functionals
Orbital Density Dependent Functionals S. Kluepfel1, P. Kluepfel1, Hildur Guðmundsdóttir1 and Hannes Jónsson1,2 1. Univ. of Iceland; 2. Aalto University Outline: Problems with GGA approximation (PBE, RPBE,...)
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 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 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 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 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 informationMODULE 2: QUANTUM MECHANICS. Practice: Quantum ESPRESSO
MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW http://www.quantum-espresso.org FREE PW-DFT, PP, PAW http://www.abinit.org FREE
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 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 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 informationDFT: Exchange-Correlation
DFT: Local functionals, exact exchange and other post-dft methods Stewart Clark University of Outline Introduction What is exchange and correlation? Quick tour of XC functionals (Semi-)local: LDA, PBE,
More 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 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 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 informationPerformance ofhybrid density functional methods,screened exchange and EXX-OEP methodsin the PAW approach p.1/26
Performance of hybrid density functional methods, screened exchange and EXX-OEP methods in the PAW approach Georg Kresse, J Paier, R Hirschl, M Marsmann Institut für Materialphysik and Centre for Computational
More informationMulti-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005
Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson
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 informationMolecular and Electronic Dynamics Using the OpenAtom Software
Molecular and Electronic Dynamics Using the OpenAtom Software Sohrab Ismail-Beigi (Yale Applied Physics) Subhasish Mandal (Yale), Minjung Kim (Yale), Raghavendra Kanakagiri (UIUC), Kavitha Chandrasekar
More informationHECToR CSE technical meeting, Oxford Parallel Algorithms for the Materials Modelling code CRYSTAL
HECToR CSE technical meeting, Oxford 2009 Parallel Algorithms for the Materials Modelling code CRYSTAL Dr Stanko Tomi Computational Science & Engineering Department, STFC Daresbury Laboratory, UK Acknowledgements
More informationExcitation Dynamics in Quantum Dots. Oleg Prezhdo U. Washington, Seattle
Excitation Dynamics in Quantum Dots Oleg Prezhdo U. Washington, Seattle Warwick August 27, 2009 Outline Time-Domain Density Functional Theory & Nonadiabatic Molecular Dynamics Quantum backreaction, surface
More informationwien2wannier and woptic: From Wannier Functions to Optical Conductivity
wien2wannier and woptic: From Wannier Functions to Optical Conductivity Elias Assmann Institute of Solid State Physics, Vienna University of Technology AToMS-2014, Bariloche, Aug 4 Outline brief introduction
More informationNumerical construction of Wannier functions
July 12, 2017 Internship Tutor: A. Levitt School Tutor: É. Cancès 1/27 Introduction Context: Describe electrical properties of crystals (insulator, conductor, semi-conductor). Applications in electronics
More informationFirst-Principles Wannier Functions of Silicon and Gallium. Arsenide arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 22 Nov 1996.
First-Principles Wannier Functions of Silicon and Gallium Arsenide arxiv:cond-mat/9611176v1 [cond-mat.mtrl-sci] 22 Nov 1996 Pablo Fernández 1, Andrea Dal Corso 1, Francesco Mauri 2, and Alfonso Baldereschi
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 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 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 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 informationQbox User Guide version
Qbox User Guide version 1.62.2 2015-08-14 François Gygi fgygi@ucdavis.edu http://eslab.ucdavis.edu 1 Introduction Qbox is a First-Principles Molecular Dynamics code. It can be used to compute the electronic
More informationSynopsis of Numerical Linear Algebra
Synopsis of Numerical Linear Algebra Eric de Sturler Department of Mathematics, Virginia Tech sturler@vt.edu http://www.math.vt.edu/people/sturler Iterative Methods for Linear Systems: Basics to Research
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 informationThe Nature of the Interlayer Interaction in Bulk. and Few-Layer Phosphorus
Supporting Information for: The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus L. Shulenburger, A.D. Baczewski, Z. Zhu, J. Guan, and D. Tománek, Sandia National Laboratories, Albuquerque,
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 informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
More informationTime-dependent density functional theory (TDDFT)
04/05/16 Hands-on workshop and Humboldt-Kolleg: Density-Functional Theory and Beyond - Basic Principles and Modern Insights Isfahan University of Technology, Isfahan, Iran, May 2 to 13, 2016 Time-dependent
More informationElectronic structure theory: Fundamentals to frontiers. VI. Analysis and more.
Electronic structure theory: Fundamentals to frontiers. VI. Analysis and more. MARTIN HEAD-GORDON Department of Chemistry, University of California, Berkeley, and, Chemical Sciences Division, Lawrence
More informationarxiv: v1 [cond-mat.str-el] 7 Dec 2008
Order-N implementation of exact exchange in extended systems Xifan Wu, Annabella Selloni, and Roberto Car Chemistry Department, Princeton University, Princeton, NJ 08544-0001,USA (Dated: October 28, 2018)
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 informationBasic introduction of NWChem software
Basic introduction of NWChem software Background NWChem is part of the Molecular Science Software Suite Designed and developed to be a highly efficient and portable Massively Parallel computational chemistry
More informationTime-dependent density functional theory (TDDFT)
Advanced Workshop on High-Performance & High-Throughput Materials Simulations using Quantum ESPRESSO ICTP, Trieste, Italy, January 16 to 27, 2017 Time-dependent density functional theory (TDDFT) Ralph
More informationPeriodic Trends in Properties of Homonuclear
Chapter 8 Periodic Trends in Properties of Homonuclear Diatomic Molecules Up to now, we have discussed various physical properties of nanostructures, namely, two-dimensional - graphene-like structures:
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 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 informationJournal of Theoretical Physics
1 Journal of Theoretical Physics Founded and Edited by M. Apostol 53 (2000) ISSN 1453-4428 Ionization potential for metallic clusters L. C. Cune and M. Apostol Department of Theoretical Physics, Institute
More informationMinimization of the Kohn-Sham Energy with a Localized, Projected Search Direction
Minimization of the Kohn-Sham Energy with a Localized, Projected Search Direction Courant Institute of Mathematical Sciences, NYU Lawrence Berkeley National Lab 29 October 2009 Joint work with Michael
More informationTutorial on DFPT and TD-DFPT: calculations of phonons and absorption spectra
Tutorial on DFPT and TD-DFPT: calculations of phonons and absorption spectra Iurii Timrov SISSA Scuola Internazionale Superiore di Studi Avanzati, Trieste Italy itimrov@sissa.it Computer modelling of materials
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 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 informationRe-design of Higher level Matrix Algorithms for Multicore and Heterogeneous Architectures. Based on the presentation at UC Berkeley, October 7, 2009
III.1 Re-design of Higher level Matrix Algorithms for Multicore and Heterogeneous Architectures Based on the presentation at UC Berkeley, October 7, 2009 Background and motivation Running time of an algorithm
More informationKey concepts in Density Functional Theory (II) Silvana Botti
Kohn-Sham scheme, band structure and optical spectra European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address:
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 informationImporting ab-initio theory into DFT: Some applications of the Lieb variation principle
Importing ab-initio theory into DFT: Some applications of the Lieb variation principle Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department
More informationReactivity and Organocatalysis. (Adalgisa Sinicropi and Massimo Olivucci)
Reactivity and Organocatalysis (Adalgisa Sinicropi and Massimo Olivucci) The Aldol Reaction - O R 1 O R 1 + - O O OH * * H R 2 R 1 R 2 The (1957) Zimmerman-Traxler Transition State Counterion (e.g. Li
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 informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
More informationc 2017 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 39, No. 6, pp. B1178 B1198 c 17 Society for Industrial and Applied Mathematics COMPUTING LOCALIZED REPRESENTATIONS OF THE KOHN SHAM SUBSPACE VIA RANDOMIZATION AND REFINEMENT ANIL
More informationCITY UNIVERSITY OF HONG KONG. Theoretical Study of Electronic and Electrical Properties of Silicon Nanowires
CITY UNIVERSITY OF HONG KONG Ë Theoretical Study of Electronic and Electrical Properties of Silicon Nanowires u Ä öä ªqk u{ Submitted to Department of Physics and Materials Science gkö y in Partial Fulfillment
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
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 informationMethod development at SUNCAT in general
Bayesian error estimation functionals and further method development at SUNCAT GPAW 2016 University of Jyväskylä June 8th, 2016 Johannes Voss vossj@stanford.edu Method development at SUNCAT in general
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 informationIntroduction to spin and spin-dependent phenomenon
Introduction to spin and spin-dependent phenomenon Institut für Theoretische Physik Freie Universität Berlin, Germany and Fritz Haber Institute of the Max Planck Society, Berlin, Germany. May 16th, 2007
More informationNorm-conserving pseudopotentials and basis sets in electronic structure calculations. Javier Junquera. Universidad de Cantabria
Norm-conserving pseudopotentials and basis sets in electronic structure calculations Javier Junquera Universidad de Cantabria Outline Pseudopotentials Why pseudopotential approach is useful Orthogonalized
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 informationQuantum Phases in Bose-Hubbard Models with Spin-orbit Interactions
Quantum Phases in Bose-Hubbard Models with Spin-orbit Interactions Shizhong Zhang The University of Hong Kong Institute for Advanced Study, Tsinghua 24 October 2012 The plan 1. Introduction to Bose-Hubbard
More informationCalculations of band structures
Chemistry and Physics at Albany Planning for the Future Calculations of band structures using wave-function based correlation methods Elke Pahl Centre of Theoretical Chemistry and Physics Institute of
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 informationPredictive Computing for Solids and Liquids
Predictive Computing for Solids and Liquids So Hirata Department of Chemistry May 214 Blue Waters Symposium 1 Schrödinger equation for a water molecule 1-particle, 3-dimensional partial differential equation
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 informationGeometry Optimisation
Geometry Optimisation Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, UK http://www.cmt.york.ac.uk/cmd http://www.castep.org Motivation Overview of Talk Background
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 informationPreconditioning. Noisy, Ill-Conditioned Linear Systems
Preconditioning Noisy, Ill-Conditioned Linear Systems James G. Nagy Emory University Atlanta, GA Outline 1. The Basic Problem 2. Regularization / Iterative Methods 3. Preconditioning 4. Example: Image
More informationSupporting information for: Electronic chemical. potentials of porous metal-organic frameworks
Supporting information for: Electronic chemical potentials of porous metal-organic frameworks Keith T. Butler, Christopher H. Hendon, and Aron Walsh Department of Chemistry, University of Bath, Claverton
More information