Comp. Mat. Science School 2001 1 Linear Scaling ‘Order-N’ Methods in Electronic Structure Theory Richard M. Martin University of Illinois Acknowledgements: Pablo Ordejon David Drabold Matthew Grumbach Uwe Stephan Daniel Sanchez-Portal Satoshi Itoh Thanks to: Jose Soler, Emilio Artacho, Giulia Galli, ...
30
Embed
Comp. Mat. Science School 20011 Linear Scaling ‘Order-N’ Methods in Electronic Structure Theory Richard M. Martin University of Illinois Acknowledgements:
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Comp. Mat. Science School 2001 1
Linear Scaling ‘Order-N’ Methodsin Electronic Structure Theory
Richard M. Martin
University of Illinois
Acknowledgements:Pablo Ordejon David Drabold
Matthew Grumbach Uwe StephanDaniel Sanchez-PortalSatoshi Itoh
Thanks to: Jose Soler, Emilio Artacho, Giulia Galli, ...
Comp. Mat. Science School 2001 2
Linear Scaling ‘Order-N’ Methodsand Car-Parrinello Simulations
• Fundamental Issues of locality in quantum mechanics
• Paradigm for view of electronic properties• Practical Algorithms • Results
Comp. Mat. Science School 2001 3
Locality in Quantum Mechanics
• V. Heine (Sol. St. Phys. Vol. 35, 1980)“Throwing out k-space”Based on ideas of Friedel (1954) , . . .
• Many properties of electrons in any region are independent of distant regions
• Walter Kohn “Nearsightness”
Comp. Mat. Science School 2001 4
Locality in Quantum Mechanics
• Which properties of electrons are independent of distant regions?
• Total integrated quantitiesDensity, Forceson atoms, . . .
• Coulomb Forces are long range but they can be handled in O(N) fashion just as in classical systems
Comp. Mat. Science School 2001 5
Non-Locality in Quantum Mechanics
• Which properties of electrons are non-local?• Individual Eigenstates in crystals• Sharp features of the Fermi surface at low T • Electrical Conductivity at T=0
Metals vs insulators: distinguished by delocalization of eigenstates at the Fermi energy (metals) vs localization of the entire many-electron system (insulators)
• Approach in the Order-N methods: Identify localized and delocalized aspects
Comp. Mat. Science School 2001 6
Density Matrix I
• Key property that describes the range of the non-locality is the density matrix (r,r’)
• In an insulator (r,r’) is exponentially localized• In a metal (r,r’) decays as a power law at T = 0, exponentially
for T > 0. (Goedecker, Ismail-Beigi)
• For non-interacting Bosons or Fermions, Landau and Lifshitz show that the correlation function g(r,r’) is uniquely related to the square of (r,r’)
• Thus correlation lengths and the density matrix generally become shorter range at high T
Comp. Mat. Science School 2001 7
Density Matrix II
• Key property that describes the range of the non-locality is the density matrix (r,r’)
• Definition: (r,r’) = i i *(r) i (r’)
• Can be localized even if each i *(r) is not!
r fixed at r =0
r’
Atom positions
Comp. Mat. Science School 2001 8
Toward Working Algorithms I
(My own personal view)
Heine and Haydock laid the groundwork - but it was applied only to limited Hamiltonians, ….
1985 - Car-Parrinello Methods changed the picture
Key quantity is the total energy E[{i}] which does not require eigenstates - only traces over the occupied states - the {i} can be linear combinations of eigenstates
Comp. Mat. Science School 2001 9
Toward Working Algorithms II
How can we use the advantages of the Car-Parrinello and the local approaches?
1992 - Galli and Parrinello pointed out the key idea -
to make a Car-Parrinello algorithm that takes advantage of the locality
Require that the states in localized.
Note this does not require a localized basis - it may be very convenient, but a localized basis is not essential to construct localized states (example: sum of plane waves can be localized)
Comp. Mat. Science School 2001 10
Toward Working Algorithms III
What are localized combinations of the eigenfunctions?
Wannier Functions (generalized)!
Wannier Functions span the same space as the eigenstates - all traces are the same
Wannier FunctionsOne localized Wannier Ftn centered on each site
Extended Bloch Eigenfunctions
Comp. Mat. Science School 2001 11
Toward Working Algorithms IV
Can work with either localized Wannier functions wi (r)
or
localized density matrix (r,r’) = i i
*(r) i (r’) = i wi *(r) wi (r’)
Functions of one variableBut not unique
Functions of two variables - more complexBut unique
Comp. Mat. Science School 2001 12
Linear Scaling ‘Order-N’ Methods
• Computational complexity ~ N = number of atoms (Current methods scale as N2 or N3)
• Intrinsically Parallel• “Divide and Conquer” • Green’s Functions• Fermi Operator Expansion• Density matrix “purification”• Generalized Wannier Functions• Spectral “Telescoping”
(Review by S. Goedecker in Rev Mod Phys)
Comp. Mat. Science School 2001 13
Divide and Conquer (Yang, 1991)
• Divide System into (Overlapping) Spatial Regions. Solve each region in terms only of its neighbors.(Terminate regions suitably)
• Use standard methods for each region
• Sum charge densities to get total density, Coulomb terms
Comp. Mat. Science School 2001 14
Expansion of the Fermi function
• Sankey, et al (1994); Goedecker, Colombo (1994); Wang et al (1995)
• Explicit T nonzero• Projection into the occupied Subspace• Multiply trial function by “Fermi operator”:
F = [(H - EF)/KBT +1]-1 • Localized leads to localized projection since the Fermi
operator (density matrix) is localized• Accomplish by expanding F in power series in H operator -
• Relaxation - ~15-60 min/step (~ 1 day with diagonalization)
Iso-density surfaces
Comp. Mat. Science School 2001 27
HOMO and LUMO in DNA (SIESTA code)
• Eigenstates found by N3 method after relaxation
• Could be O(N) for each state
Comp. Mat. Science School 2001 28
O(N) Simulation of Magnets at T > 0
• Collaboration at ORNL, Ames, Brookhaven• Snapshot of magnetic order in a finite temperature simulation of paramagnetic Fe. These
calculations represent significant progress towards the goal of full implementation of a first principles theory of the finite temperature and non-equilibrium properties of magnetic materials.
• Record setting performance for large unit cell models (up to 1024-atoms) led to the award of the 1998 Gordon Bell prize.
• The calculations that were the basis for the award were performed using the locally self-consistent multiple scattering method, which is an O(N) Density Functional method
• Web Site: http://oldpc.ms.ornl.gov/~gms/MShome.html
Comp. Mat. Science School 2001 29
FUTURE! ---- Biological Systems• Examples of oriented pigment molecules that today are being simulated by empirical potentials
Comp. Mat. Science School 2001 30
Conclusions• It is possible to treat many thousands of atoms in a
full simulation - on a workstation with approximate methods - intrinsically parallel for a supercomputer
• Why treat many thousands of atoms?
• Large scale structures in materials - defects, boundaries, ….
• Biological molecules
• The ideas are also relevant to understanding even small systems