O (N) Methods in Electronic Structure Calculations

Total Page:16

File Type:pdf, Size:1020Kb

O (N) Methods in Electronic Structure Calculations O(N) Methods in electronic structure calculations D R Bowler1;2;3 and T Miyazaki4 1London Centre for Nanotechnology, UCL, 17-19 Gordon St, London WC1H 0AH, UK 2Department of Physics & Astronomy, UCL, Gower St, London WC1E 6BT, UK 3Thomas Young Centre, UCL, Gower St, London WC1E 6BT, UK 4National Institute for Materials Science, 1-2-1 Sengen, Tsukuba, Ibaraki 305-0047, JAPAN E-mail: [email protected] E-mail: [email protected] Abstract. Linear scaling methods, or O(N) methods, have computational and memory requirements which scale linearly with the number of atoms in the system, N, in contrast to standard approaches which scale with the cube of the number of atoms. These methods, which rely on the short-ranged nature of electronic structure, will allow accurate, ab initio simulations of systems of unprecedented size. The theory behind the locality of electronic structure is described and related to physical properties of systems to be modelled, along with a survey of recent developments in real-space methods which are important for efficient use of high performance computers. The linear scaling methods proposed to date can be divided into seven different areas, and the applicability, efficiency and advantages of the methods proposed in these areas is then discussed. The applications of linear scaling methods, as well as the implementations available as computer programs, are considered. Finally, the prospects for and the challenges facing linear scaling methods are discussed. Submitted to: Rep. Prog. Phys. arXiv:1108.5976v5 [cond-mat.mtrl-sci] 3 Nov 2011 O(N) Methods 2 1. Introduction Electronic structure calculation methods based on the density functional theory (DFT) have been playing important roles in condensed matter physics for more than forty years. In the early stages, DFT calculations were employed mainly for the study of the electronic structure of simple solids, using a few atoms in a unit cell, with the use of periodic boundary conditions. Since then, there has been a huge effort to improve the accuracy and efficiency of the calculation techniques. In terms of efficiency, after the pioneering work by Car and Parrinello [1], the size of the target systems has increased dramatically and more and more examples of the DFT studies, especially on aperiodic systems like surface structures, have emerged. DFT calculations on systems containing hundreds of atoms are currently ubiquitous. As the system size for DFT calculations has become larger, the variety of materials and phenomena investigated by the method has increased. The information of the total energy and atomic forces calculated by DFT methods can provide reliable data independently from experiments, and the methods are nowadays considered as one of the established research tools in many fields, like physics, chemistry, materials science, and many others. Recently, there have been DFT studies in the complex fields of nano-structured materials and biological systems. In the study of these classes of materials, we need to treat systems containing at least thousands of atoms. However, as is well known, once the number of atoms N in a system reaches around one thousand, the cost of standard DFT calculations increases very rapidly as a cube of N. To overcome this problem, the methods known as linear-scaling or O(N) DFT methods have been developed [2]. The progress of these methods in the last ten to fifteen years is remarkable and the purpose of this review paper is to overview the recent progress of O(N) DFT methods. We will start with an overview of the conventional DFT method and its advantages. In the normal DFT approach, we solve for the Kohn-Sham (KS) orbitals Ψνk(r), which are the eigenstates of the KS equation [3]. " h¯ # H^ KSΨ (r) = − r2 + V (r) + V (r) + V (r) Ψ (r) = Ψ (r)(1) νk 2m ext H XC νk νk νk Here H^ KS is the Kohn-Sham Hamiltonian, and ν and k are the band index and k points in the first Brillouin zone, respectively. Hereafter, we omit k for clarity because we consider large systems and the number of k points is small. Vext(r) is the potential from nuclei, VH (r) is the Hartree potential, and VXC(r) is the exchange-correlation potential in the Kohn-Sham formalism. The most accurate DFT calculations often use a plane-wave basis set to express the KS orbitals: X Ψν(r) = cν(G) exp(iG · r) (2) jGj<Gmax A plane-wave basis set has two main advantages. First, the accuracy of the basis set can be systematically improved. In Eq. (2), Gmax is obtained from the cutoff energy Ecut 2 2 ¯h Gmax as 2m = Ecut. The number of plane-waves, NG, is controlled only by the number O(N) Methods 3 Ecut. The accuracy of the basis set can be improved simply by increasing Ecut, and a variational principle with respect to Ecut is satisfied. The other advantage is that forces can be calculated easily without the Pulay correction term because the basis set is independent on atomic positions (though such basis-set-dependent corrections become necessary when changing the unit cell size or shape). These two advantages make it possible to calculate both energy and forces accurately with plane-wave basis sets. In order to realise accurate plane-wave calculations, we need to introduce several theoretical techniques. First of all, plane-wave calculations rely on the idea of pseudopotentials [4]. With this method, it is possible to work only with valence electrons and their pseudo-wavefunctions, which are much smoother than the real wave functions which oscillate strongly in the core region, and to replace the nuclear potential and the core electrons with a pseudopotential. There have been several kinds of techniques proposed to make pseudo-wavefunctions smoother [5{7]. Using the method of ultra-soft pseudopotentials [8], even the cutoff energy for the localised 3d orbitals of transition metals can be reduced dramatically. With these improvements in theoretical techniques, the total energy converges quickly with respect to the cutoff energy and this is essential to make the accurate DFT calculations feasible. In addition, the major part of the error in the total energy usually comes from the expression of KS orbitals in the core region. Hence, the relative energetic stability of two states (e.g. two different atomic structures) can be reproduced without the absolute convergence because most of the errors are cancelled in the energy difference. Note that it is also possible to reduce the number of plane-waves by using augmentation for the wavefunctions in the core region as in the linearised augmented plane-wave (LAPW) or the projector augmented wave (PAW) method [9]. It is essential that we reduce the number of plane-waves by introducing the pseudopotential or other similar techniques. However, even with very smooth pseudo wavefunctions, NG is typically one hundred times larger than the number of electrons. KS 0 2 When we want to diagonalise < GjH jG >, the required memory scales as O(NG) and 3 CPU time as O(NG). Hence it is impossible to employ direct (exact) diagonalisation except for very small systems. Instead of using exact diagonalisation, we can obtain the Kohn-Sham orbitals by minimising the DFT total energy with respect to the coefficients fcν(G)g, as shown in the work by Car and Parrinello [1]. Since we only need the occupied Kohn-Sham orbitals in such iterative methods, the memory requirement to 2 store fcν(G)g is proportional to NBNG, which is roughly 100 times smaller than NG. Then, we update the coefficients fcν(G)g by calculating the gradient of the total energy with a constraint to keep the KS orbitals orthogonal to each other. This is done by KS KS calculating (H −Λν,ν0 ) or H −ν with Gram-Schmidt orthogonalisation of fcν(G)g. In the calculation of the Kohn-Sham Hamiltonian, we need to calculate the density n(r). For this, we first calculate n(r) as X ∗ n(r) = fνΨν(r)Ψν(r) (3) ν 2 If we perform this in a straightforward way, we need the operations of O(NG) for each O(N) Methods 4 band ν, and the total number of operations needed for the transformation from fcν(G)g 2 to Ψν(r) is of the order of NBNG, which is quite expensive. However, we can dramatically reduce the cost of the calculation using the fast Fourier transform (FFT) method, and the number of operations in Eq. (3) becomes NBNG ln(NG). Although Eq. (3) is still the most expensive part in the calculations of small systems for many plane-wave DFT codes, the reduction of the computational cost by the FFT method is essential for the success of plane-wave DFT calculations. The orthogonalisation of the KS orbitals is also an expensive operation, which R 0 includes the calculations drΨν(r)Ψν0 (r) for all pairs of band indices fν; ν g. The total 2 cost of the operations is O(NBNG), but we can see that it is only proportional to NG. As we have seen, in the iterative method with the plane-wave basis set, there are no 2 operations where the cost increases as fast as NG. This is the reason why we can do efficient calculations even with large NG. The iterative diagonalisation technique, FFTs, and ab initio pseudopotentials used in the plane-wave calculations are the key factors which make it possible to employ accurate but efficient DFT calculations. Using these techniques, with the increase of the computer power, the time for solving KS equations has become smaller and smaller, and the system size for the target of DFT studies has become larger and larger.
Recommended publications
  • Introducing ONETEP: Linear-Scaling Density Functional Simulations on Parallel Computers Chris-Kriton Skylaris,A) Peter D
    THE JOURNAL OF CHEMICAL PHYSICS 122, 084119 ͑2005͒ Introducing ONETEP: Linear-scaling density functional simulations on parallel computers Chris-Kriton Skylaris,a) Peter D. Haynes, Arash A. Mostofi, and Mike C. Payne Theory of Condensed Matter, Cavendish Laboratory, Madingley Road, Cambridge CB3 0HE, United Kingdom ͑Received 29 September 2004; accepted 4 November 2004; published online 23 February 2005͒ We present ONETEP ͑order-N electronic total energy package͒, a density functional program for parallel computers whose computational cost scales linearly with the number of atoms and the number of processors. ONETEP is based on our reformulation of the plane wave pseudopotential method which exploits the electronic localization that is inherent in systems with a nonvanishing band gap. We summarize the theoretical developments that enable the direct optimization of strictly localized quantities expressed in terms of a delocalized plane wave basis. These same localized quantities lead us to a physical way of dividing the computational effort among many processors to allow calculations to be performed efficiently on parallel supercomputers. We show with examples that ONETEP achieves excellent speedups with increasing numbers of processors and confirm that the time taken by ONETEP as a function of increasing number of atoms for a given number of processors is indeed linear. What distinguishes our approach is that the localization is achieved in a controlled and mathematically consistent manner so that ONETEP obtains the same accuracy as conventional cubic-scaling plane wave approaches and offers fast and stable convergence. We expect that calculations with ONETEP have the potential to provide quantitative theoretical predictions for problems involving thousands of atoms such as those often encountered in nanoscience and biophysics.
    [Show full text]
  • Natural Bond Orbital Analysis in the ONETEP Code: Applications to Large Protein Systems Louis P
    WWW.C-CHEM.ORG FULL PAPER Natural Bond Orbital Analysis in the ONETEP Code: Applications to Large Protein Systems Louis P. Lee,*[a] Daniel J. Cole,[a] Mike C. Payne,[a] and Chris-Kriton Skylaris[b] First principles electronic structure calculations are typically Generalized Wannier Functions of ONETEP to natural atomic performed in terms of molecular orbitals (or bands), providing a orbitals, NBO analysis can be performed within a localized straightforward theoretical avenue for approximations of region in such a way that ensures the results are identical to an increasing sophistication, but do not usually provide any analysis on the full system. We demonstrate the capabilities of qualitative chemical information about the system. We can this approach by performing illustrative studies of large derive such information via post-processing using natural bond proteins—namely, investigating changes in charge transfer orbital (NBO) analysis, which produces a chemical picture of between the heme group of myoglobin and its ligands with bonding in terms of localized Lewis-type bond and lone pair increasing system size and between a protein and its explicit orbitals that we can use to understand molecular structure and solvent, estimating the contribution of electronic delocalization interactions. We present NBO analysis of large-scale calculations to the stabilization of hydrogen bonds in the binding pocket of with the ONETEP linear-scaling density functional theory package, a drug-receptor complex, and observing, in situ, the n ! p* which we have interfaced with the NBO 5 analysis program. In hyperconjugative interactions between carbonyl groups that ONETEP calculations involving thousands of atoms, one is typically stabilize protein backbones.
    [Show full text]
  • Metalloboranes from First-Principles Calculations: a Candidate for High-Density Hydrogen Storage
    Metalloboranes from first-principles calculations: A candidate for high-density hydrogen storage A. R. Akbarzadeh, D. Vrinceanu, and C.J. Tymczak Department of Physics, Texas Southern University, Houston, Texas 77004, USA Abstract Using first principles calculations, we show the high hydrogen storage capacity of a new class of compounds, metalloboranes. Metalloboranes are transition metal (TM) and borane compounds that obey a novel-bonding scheme. We have found that the transition metal atoms can bind up to 10 H2-molecules with an average binding energy of 30 kJ/mole of H2, which lies favorably within the reversible adsorption range. Among the first row TM atoms, Sc and Ti are found to be the optimum in maximizing the H2 storage on the metalloborane cluster. Additionally, being ionically bonded to the borane molecule, the TMs do not suffer from the aggregation problem, which has been the biggest hurdle for the success of TM- decorated graphitic materials for hydrogen storage. Furthermore, since the borane 6-atom ring has identical bonding properties as carbon rings it is possible to link the metalloboranes into metal organic frameworks (MOF’s), which are thus able to adsorb hydrogen via Kubas interaction as well as the well-known van der Waals interaction. Finally, we construct a simple Monte-Catlo algorithm for Hydrogen uptake and show that Titanium metalloboranes in a MOF5 structure can absorb up to 11.5% hydrogen per weight at 100 bar of pressure. I. Introduction Due to the geographical distribution and supply limitation of fossil fuel resources and the increasing perceived negative impact of carbon dioxide emission in the environment, it is advantageous for the societies to take action in the search and implementation of alternative energy systems.
    [Show full text]
  • Density Functional Theory (DFT)
    Herramientas mecano-cuánticas basadas en DFT para el estudio de moléculas y materiales en Materials Studio 7.0 Javier Ramos Biophysics of Macromolecular Systems group (BIOPHYM) Departamento de Física Macromolecular Instituto de Estructura de la Materia – CSIC [email protected] Webinar, 26 de Junio 2014 Anteriores webinars Como conseguir los videos y las presentaciones de anteriores webminars: Linkedin: Grupo de Química Computacional http://www.linkedin.com/groups/Química-computacional-7487634 Índice Density Functional Theory (DFT) The Jacob’s ladder DFT modules in Maretials Studio DMOL3, CASTEP and ONETEP XC functionals Basis functions Interfaces in Materials Studio Tasks Properties Example: n-butane conformations Density Functional Theory (DFT) DFT is built around the premise that the energy of an electronic system can be defined in terms of its electron probability density (ρ). (Hohenberg-Kohn Theorem) E 0 [ 0 ] Te [ 0 ] E ne [ 0 ] E ee [ 0 ] (easy) Kinetic Energy for ????? noninteracting (r )v (r ) dr electrons(easy) 1 E[]()()[]1 r r d r d r E e e2 1 2 1 2 X C r12 Classic Term(Coulomb) Non-classic Kohn-Sham orbitals Exchange & By minimizing the total energy functional applying the variational principle it is Correlation possible to get the SCF equations (Kohn-Sham) The Jacob’s Ladder Accurate form of XC potential Meta GGA Empirical (Fitting to Non-Empirical Generalized Gradient Approx. atomic properties) (physics rules) Local Density Approximation DFT modules in Materials Studio DMol3: Combine computational speed with the accuracy of quantum mechanical methods to predict materials properties reliably and quickly CASTEP: CASTEP offers simulation capabilities not found elsewhere, such as accurate prediction of phonon spectra, dielectric constants, and optical properties.
    [Show full text]
  • Compact Orbitals Enable Low-Cost Linear-Scaling Ab Initio Molecular Dynamics for Weakly-Interacting Systems Hayden Scheiber,1, A) Yifei Shi,1 and Rustam Z
    Compact orbitals enable low-cost linear-scaling ab initio molecular dynamics for weakly-interacting systems Hayden Scheiber,1, a) Yifei Shi,1 and Rustam Z. Khaliullin1, b) Department of Chemistry, McGill University, 801 Sherbrooke St. West, Montreal, QC H3A 0B8, Canada Today, ab initio molecular dynamics (AIMD) relies on the locality of one-electron density matrices to achieve linear growth of computation time with systems size, crucial in large-scale simulations. While Kohn-Sham orbitals strictly localized within predefined radii can offer substantial computational advantages over density matrices, such compact orbitals are not used in AIMD because a compact representation of the electronic ground state is difficult to find. Here, a robust method for maintaining compact orbitals close to the ground state is coupled with a modified Langevin integrator to produce stable nuclear dynamics for molecular and ionic systems. This eliminates a density matrix optimization and enables first orbital-only linear-scaling AIMD. An application to liquid water demonstrates that low computational overhead of the new method makes it ideal for routine medium-scale simulations while its linear-scaling complexity allows to extend first- principle studies of molecular systems to completely new physical phenomena on previously inaccessible length scales. Since the unification of molecular dynamics and den- LS methods restrict their use in dynamical simulations sity functional theory (DFT)1, ab initio molecular dy- to very short time scales, systems of low dimensions, namics (AIMD) has become an important tool to study and low-quality minimal basis sets6,18–20. On typical processes in molecules and materials. Unfortunately, the length and time scales required in practical and accurate computational cost of the conventional Kohn-Sham (KS) AIMD simulations, LS DFT still cannot compete with DFT grows cubically with the number of atoms, which the straightforward low-cost cubically-scaling KS DFT.
    [Show full text]
  • Quantum Chemistry (QC) on Gpus Feb
    Quantum Chemistry (QC) on GPUs Feb. 2, 2017 Overview of Life & Material Accelerated Apps MD: All key codes are GPU-accelerated QC: All key codes are ported or optimizing Great multi-GPU performance Focus on using GPU-accelerated math libraries, OpenACC directives Focus on dense (up to 16) GPU nodes &/or large # of GPU nodes GPU-accelerated and available today: ACEMD*, AMBER (PMEMD)*, BAND, CHARMM, DESMOND, ESPResso, ABINIT, ACES III, ADF, BigDFT, CP2K, GAMESS, GAMESS- Folding@Home, GPUgrid.net, GROMACS, HALMD, HOOMD-Blue*, UK, GPAW, LATTE, LSDalton, LSMS, MOLCAS, MOPAC2012, LAMMPS, Lattice Microbes*, mdcore, MELD, miniMD, NAMD, NWChem, OCTOPUS*, PEtot, QUICK, Q-Chem, QMCPack, OpenMM, PolyFTS, SOP-GPU* & more Quantum Espresso/PWscf, QUICK, TeraChem* Active GPU acceleration projects: CASTEP, GAMESS, Gaussian, ONETEP, Quantum Supercharger Library*, VASP & more green* = application where >90% of the workload is on GPU 2 MD vs. QC on GPUs “Classical” Molecular Dynamics Quantum Chemistry (MO, PW, DFT, Semi-Emp) Simulates positions of atoms over time; Calculates electronic properties; chemical-biological or ground state, excited states, spectral properties, chemical-material behaviors making/breaking bonds, physical properties Forces calculated from simple empirical formulas Forces derived from electron wave function (bond rearrangement generally forbidden) (bond rearrangement OK, e.g., bond energies) Up to millions of atoms Up to a few thousand atoms Solvent included without difficulty Generally in a vacuum but if needed, solvent treated classically
    [Show full text]
  • Introduction to DFT and the Plane-Wave Pseudopotential Method
    Introduction to DFT and the plane-wave pseudopotential method Keith Refson STFC Rutherford Appleton Laboratory Chilton, Didcot, OXON OX11 0QX 23 Apr 2014 Parallel Materials Modelling Packages @ EPCC 1 / 55 Introduction Synopsis Motivation Some ab initio codes Quantum-mechanical approaches Density Functional Theory Electronic Structure of Condensed Phases Total-energy calculations Introduction Basis sets Plane-waves and Pseudopotentials How to solve the equations Parallel Materials Modelling Packages @ EPCC 2 / 55 Synopsis Introduction A guided tour inside the “black box” of ab-initio simulation. Synopsis • Motivation • The rise of quantum-mechanical simulations. Some ab initio codes Wavefunction-based theory • Density-functional theory (DFT) Quantum-mechanical • approaches Quantum theory in periodic boundaries • Plane-wave and other basis sets Density Functional • Theory SCF solvers • Molecular Dynamics Electronic Structure of Condensed Phases Recommended Reading and Further Study Total-energy calculations • Basis sets Jorge Kohanoff Electronic Structure Calculations for Solids and Molecules, Plane-waves and Theory and Computational Methods, Cambridge, ISBN-13: 9780521815918 Pseudopotentials • Dominik Marx, J¨urg Hutter Ab Initio Molecular Dynamics: Basic Theory and How to solve the Advanced Methods Cambridge University Press, ISBN: 0521898633 equations • Richard M. Martin Electronic Structure: Basic Theory and Practical Methods: Basic Theory and Practical Density Functional Approaches Vol 1 Cambridge University Press, ISBN: 0521782856
    [Show full text]
  • Presidential Search Committee Meets, Reviews Upcoming Work
    ETHEMedford, MA 02155 TUFTSTuesday, November 26,1991 DAILYJVol XXIII, Number 58 Faculty tables vote on Presidential search committee World Civ proposal meets, reviews upcoming work courses iUC to be intcrdiscipli by SCHAEF1cR nary. intcgr;ltblg disciplines fr()n by MAUREEN LENIHAN Senior Staff Writer Daily Editorial Board and PATRICK HEALY at leiist three of the five distribu The 19-member search com- Dally Edil(~r1alhtUd tion areas. mittee fonned to begin the pro- Although the proposed World Since 1986. threeteamsoffac cess of selecting anew university CivilizationsFoundation rcquire- ult y. including one professor frorr president met for the first time inent was briefly presented at the the School ofEngineering,devel yesterday, familiariz,ing them- Liberal Arts and Jackson facully oped three two-semester course: selves with their task and the meeting yesterday. the LAM fac- which were ultiinately opened tc committee‘s initial duties includ- ulty postponed further discussion students to take on a voluntce ing formali~inga Job description and voting on the prognun until basis. These courses, titled “P for the post. Dec. 11. Sense of Place.” “Time and Cal The committee. representing Dean of College of Liberal endars,“and “Memory and Iden trustees. administrators, faculty Arts and Jackson Mary Ella tity.“wereestablishedwith fincan members. staff members. stu- Fcinleibmoved tocrcate the Dec. cia1 support from ihe Davis Foun dents. and alumni, was fonncd in 1 I meeting since other meeting dation. the National Endowmen response toai miouiiccmcnt last topics preceding the World Civ for the Humimities, and the Uni May that University president presentation took up most of the versit y.
    [Show full text]
  • What's New in Biovia Materials Studio 2020
    WHAT’S NEW IN BIOVIA MATERIALS STUDIO 2020 Datasheet BIOVIA Materials Studio 2020 is the latest release of BIOVIA’s predictive science tools for chemistry and materials science research. Materials Studio empowers researchers to understand the relationships between a material’s mo- lecular or crystal structure and its properties in order to make more informed decisions about materials research and development. More often than not materials performance is influenced by phenomena at multiple scales. Scientists using Materials Studio 2020 have an extensive suite of world class solvers and parameter sets operating from atoms to microscale for simulating more materials and more properties than ever before. Furthermore, the predicted properties can now be used in multi-physics modeling and in systems modeling software such as SIMULIA Abaqus and CATIA Dymola to predict macroscopic behaviors. In this way multiscale simulations can be used to solve some of the toughest pro- blems in materials design and product optimization. BETTER MATERIALS - BETTER BATTERIES Safe, fast charging batteries with high energy density and long life are urgently needed for a host of applications - not least for the electrification of all modes of transportation as an alternative to fossil fuel energy sources. Battery design relies on a complex interplay between thermal, mechanical and chemical processes from the smallest scales of the material (electronic structure) through to the geometry of the battery cell and pack design. Improvements to the component materials used in batteries and capacitors are fundamental to providing the advances in performance needed. Materials Studio provides new functionality to enable the simula- tion of key materials parameters for both liquid electrolytes and electrode components.
    [Show full text]
  • Porting the DFT Code CASTEP to Gpgpus
    Porting the DFT code CASTEP to GPGPUs Toni Collis [email protected] EPCC, University of Edinburgh CASTEP and GPGPUs Outline • Why are we interested in CASTEP and Density Functional Theory codes. • Brief introduction to CASTEP underlying computational problems. • The OpenACC implementation http://www.nu-fuse.com CASTEP: a DFT code • CASTEP is a commercial and academic software package • Capable of Density Functional Theory (DFT) and plane wave basis set calculations. • Calculates the structure and motions of materials by the use of electronic structure (atom positions are dictated by their electrons). • Modern CASTEP is a re-write of the original serial code, developed by Universities of York, Durham, St. Andrews, Cambridge and Rutherford Labs http://www.nu-fuse.com CASTEP: a DFT code • DFT/ab initio software packages are one of the largest users of HECToR (UK national supercomputing service, based at University of Edinburgh). • Codes such as CASTEP, VASP and CP2K. All involve solving a Hamiltonian to explain the electronic structure. • DFT codes are becoming more complex and with more functionality. http://www.nu-fuse.com HECToR • UK National HPC Service • Currently 30- cabinet Cray XE6 system – 90,112 cores • Each node has – 2×16-core AMD Opterons (2.3GHz Interlagos) – 32 GB memory • Peak of over 800 TF and 90 TB of memory http://www.nu-fuse.com HECToR usage statistics Phase 3 statistics (Nov 2011 - Apr 2013) Ab initio codes (VASP, CP2K, CASTEP, ONETEP, NWChem, Quantum Espresso, GAMESS-US, SIESTA, GAMESS-UK, MOLPRO) GS2NEMO ChemShell 2%2% SENGA2% 3% UM Others 4% 34% MITgcm 4% CASTEP 4% GROMACS 6% DL_POLY CP2K VASP 5% 8% 19% http://www.nu-fuse.com HECToR usage statistics Phase 3 statistics (Nov 2011 - Apr 2013) 35% of the Chemistry software on HECToR is using DFT methods.
    [Show full text]
  • HPC Issues for DFT Calculations
    HPC Issues for DFT Calculations Adrian Jackson EPCC Scientific Simulation • Simulation fast becoming 4 th pillar of science – Observation, Theory, Experimentation, Simulation • Explore universe through simulation rather than experimentation – Test theories – Predict or validate experiments – Simulate “untestable” science • Reproduce “real world” in computers – Generally simplified – Dimensions and timescales restricted – Simulation of scientific problem or environment – Input of real data – Output of simulated data – Parameter space studies – Wide range of approaches http://www.nu-fuse.com Reduce runtime • Serial code optimisations – Reduce runtime through efficiencies – Unlikely to produce required savings • Upgrade hardware – 1965: Moore’s law predicts growth in complexity of processors – Doubling of CPU performance – Performance often improved through on chip parallelism http://www.nu-fuse.com Parallel Background • Why not just make a faster chip? – Theoretical • Physical limitations to size and speed of a single chip • Capacitance increases with complexity • Speed of light, size of atoms, dissipation of heat • The power used by a CPU core is proportional to Clock Frequency x Voltage 2 • Voltage reduction vs Clock speed for power requirements – Voltages become too small for “digital” differences – Practical • Developing new chips is incredibly expensive • Must make maximum use of existing technology http://www.nu-fuse.com Parallel Systems • Different types of parallel systems P M – Shared memory P M P M – Distributed memory P M Interconnect
    [Show full text]
  • Curriculum Vitae
    CURRICULUM VITAE Christopher J. Tymczak Tenured Associate Professor Director Department of Physics Texas Southern University High Texas Southern University Performance Computing Center Houston, Texas 77004 Texas Southern University Office: (713) 313-1849 Houston, Texas 77004 Email: [email protected] (http://hpcc.tsu.edu/) (http://physics.tsu.edu/) Co-Director Visiting Research Scientist NSF CREST Center for Research on Group T-1 Complex Networks Los Alamos National Laboratory Texas Southern University Los Alamos, New Mexico, 87507 Houston, Texas 77004 I specialize in the identification, development, and implementation of new scientific codes for exploiting advanced computing resources impacting large scale computation in diverse areas in Many-body physics, Quantum Chemistry and ab initio Molecular Dynamics. I am a tenured associate professor of physics at Texas Southern University, where I am spearheading the integration of supercomputing resources into various STEM programs. I am also the founder and director of the Texas Southern High Performance Computing Center (TSU-HPCC). At Los Alamos National Laboratory, I am a permanent scientific associate within the FreeON initiative, involving the development of massively parallel linear scaling quantum chemistry methods, currently under development in collaboration with Dr. Matt Challacombe (T-12) and Dr. Anders Niklasson (T-1). I was one of the first to exploit wavelet-based methods in large scale computing for understanding the electronic structure of materials. For the last six years, I have advanced the FreeON development through the exploitation of advanced data structures and advanced machine architectures. FreeON is now recognized as one of the first quantum chemical codes with demonstrable scalable parallelism within large- scale parallel clusters.
    [Show full text]