What Is the Price of Open-Source Software? 8 9 Anna I

Total Page:16

File Type:pdf, Size:1020Kb

What Is the Price of Open-Source Software? 8 9 Anna I The Journal of Physical Chemistry Letters This document is confidential and is proprietary to the American Chemical Society and its authors. Do not copy or disclose without written permission. If you have received this item in error, notify the sender and delete all copies. What is the Price of Open -Source Software? Journal: The Journal of Physical Chemistry Letters Manuscript ID: jz-2015-01258h Manuscript Type: Viewpoint Date Submitted by the Author: 13-Jun-2015 Complete List of Authors: Krylov, Anna; University of Southern California, Dept. of Chemistry Herbert, John; The Ohio State University, Chemistry Furche, Filipp; University of California, Irvine, Chemistry Head-Gordon, Martin; University of California, Berkeley, Chemistry Knowles, Peter; Cardiff University, School of Chemistry Lindh, Roland; Chemistry - Ångström Laboratory, Theoretical Chemistry Manby, Frederick; University of Bristol, School of Chemistry Pulay, Peter; University of Arkansas, Chemistry and Biochemistry Skylaris, Chris-Kriton; University of Southampton, Werner, Hans-Joachim; Universitaet Stuttgart, Institut fuer Theoretische Chemie ACS Paragon Plus Environment Page 1 of 6 The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 What is the price of open-source software? 8 9 Anna I. Krylov,1∗ John M. Herbert,2† Filipp Furche,3 10 Martin Head-Gordon,4 Peter J. Knowles,5 Roland Lindh,6 Frederick R. Manby,7 11 8 9 10 12 Peter Pulay, Chris-Kriton Skylaris, and Hans-Joachim Werner 13 14 1Dept. of Chemistry, University of Southern California, Los Angeles, California, USA 15 2Dept. of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio, USA 16 3Dept. of Chemistry, University of California, Irvine, California, USA 17 4 18 Dept. of Chemistry, University of California, Berkeley, California, USA 5 19 School of Chemistry, Cardiff University, Cardiff, UK 20 6Dept. of Chemistry–Angstr¨omLaboratory,˚ Uppsala University, Uppsala, Sweden 21 7Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, UK 22 8Dept. of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas, USA 23 9 24 School of Chemistry, University of Southampton, Highfield, Southampton, UK 25 10Institut f¨urTheoretische Chemie, Universit¨atStuttgart, Stuttgart, Germany 26 27 June 12, 2015 28 29 30 The notion that all scientific software should be open-source and free has been actively promoted in 31 recent years, mostly from the top down via mandates from funding agencies [1] but occasionally from the 32 bottom up, as exemplified by a recent Viewpoint in this journal. [2] A commonly articulated rationale is that 33 the results of scientific research funded by government grants should be free for society and that the scientific 34 community benefits from free access. The purpose of this Viewpoint is to examine the consequences of these 35 opinions. 36 What is scientific software? Modern computational chemistry software is an extremely complex 37 product based on advanced scientific ideas (models and theories) and sophisticated algorithms that transform 38 these ideas from equations into useful tools. The development of practical software that can be used by 39 non-experts to solve contemporary research problems requires considerable technical effort to produce and 40 maintain robust, efficient, and validated code. Unlike the development of, e.g., a smart-phone app, where the 41 code base is small [3] and a relatively large community can easily write extensions and add-ons, production 42 of scientific software involves the curation of millions of lines of source code. The complexity of this code 43 demands long-term user and developer support to maintain its integrity and performance while keeping 44 45 up with new computer architectures, fixing bugs, and adding features. Recognizing the importance of these ideas, various funding agencies in the U.S. have made “sustainable software” a key priority in the distribution 46 [1] 47 of research support. Sustainability is a critical goal, but one that can be realized in various ways. 48 Good software is important to science. Computational chemistry software is an essential scientific 49 instrument that facilitates discovery and innovation far beyond the laboratories in which it is created, an [4] 50 achievement that was recognized by the 1998 and 2013 Nobel Prizes in Chemistry. Focusing on quantum 51 chemistry software in this Viewpoint, we note that today any chemist can (with very little training) use nu- 52 merous quantum chemistry programs as teaching and research tools that aid in the design and interpretation 53 of experiments. 54 ∗[email protected] 55 †[email protected] 56 57 58 1 59 60 ACS Paragon Plus Environment The Journal of Physical Chemistry Letters Page 2 of 6 1 2 3 A software package should be more than just a tool for end users, however; it should also be a platform 4 to develop and test new models and algorithms. Maintaining a code base requires extensive validation, and 5 given the complexity of modern computational methods, even testing of “pilot code” or a “proof-of-principle” 6 implementation requires access to basic software infrastructure, e.g., an integrals library, a self-consistent 7 field procedure, efficient I/O and memory management, tools for manipulating tensors, etc. Modularity 8 is a laudable goal, but in reality “interoperability” often comes at the expense of performance. In high- 9 performance codes, the aforementioned components are tightly interwoven, to the extent that expert help is 10 11 often required to modify key components or to develop non-standard interfaces to them. As such, the ability 12 to innovate along either applied or theoretical lines depends crucially on the quality of the software and the availability of documentation and expert support. 13 [5] [6] 14 As examples, consider two widely-used electronic structure programs, Q-Chem and Molpro. These ∼ . ∼ . 15 codes consist of 5 5 and 2 5 million lines of source code, respectively, written in multiple languages and 16 each in continuous development over several decades. Q-Chem incorporates scientific advances reported in 17 more than 300 peer-reviewed scientific publications, while methods implemented for the first time in Molpro 18 have led to 20 high-impact papers that have each been cited over 300 times. Neither code is static: more 19 than 70 scientists are actively contributing to Molpro, and the Q-Chem developer base numbers more than 20 100. Such agile innovation comes at a price, however. Significant effort is required to keep the code robust, 21 efficient, and sound, and to provide the documentation that ensures the usability of new methods and the 22 extensibility of older ones. 23 Software from academia is often developed with an emphasis on ideas rather than implementation, fed by 24 the need for timely peer-reviewed journal publications that provide ongoing grant support and future jobs 25 for graduate students. To bring new ideas to the production level, with software that is accessible to (and 26 useful for) the broader scientific community, contributions from expert programmers are required. These 27 technical tasks usually cannot—and generally should not—be conducted by graduate students or postdocs, 28 who should instead be focused on science and innovation. To this end, Q-Chem employs four scientific 29 programmers. Other quantum chemistry codes (e.g., Molpro, [6] Turbomole, [7] Jaguar, [8] Molcas, [9] 30 PQS, [10] and ONETEP [11]) face the same challenges and adopt similar models to ensure sustainability. 31 It is important to distinguish these academically-led software ventures from purely commercial endeavors. 32 The large majority of the code in a package like Q-Chem is funded by the government, either through grants 33 to academic groups or, in some cases, through technology grants to the company itself. The role of the 34 company programmers is to enable sustainability through bug fixes, user support, release management, and 35 the addition of features that academic developers either cannot or will not add themselves. Programmers 36 employed by the company place emphasis on functionality, robustness, and performance, more so than 37 scientific innovation. They are directly addressing the “reproducibility problem”. [2] Sales revenue cannot 38 support the entire development cost of an academic code, but it contributes critically to its sustainability. 39 The cost that the customer pays for a code like Q-Chem reflects this funding model: it is vastly lower than 40 the development cost, particularly for academic customers but also for industry. It primarily reflects the 41 sustainability cost. 42 Software is not data. In his Viewpoint, [2] Gezelter argues that both software and data should be 43 open, yet it is important not to conflate the two. Software is not data, and simply because it is feasible to 44 put software on the internet does not imply that it should be posted. Software is a product that contains an 45 intellectual component (models and algorithms) but owes its existence to additional technical efforts. Such 46 efforts include implementation of minor but useful (or requested) features that are not publishable in the 47 peer-reviewed literature. This is not to say that details should be withheld as proprietary information. Just 48 as models and algorithms should be described in full detail in scientific publications, so too should imple- 49 mentation details be specified, along with performance metrics (timings and scaling data) and benchmarks 50 (energies and other computed properties). Nevertheless, the software itself is a product, not a scientific find- 51 52 ing, more akin to, say, an NMR spectrometer—a sophisticated instrument—than to the spectra produced 53 by that instrument. 54 Consider an analogy from the field of photovoltaics. Scientific findings concerning the mechanistic details 55 of charge generation and exciton propagation in a given material are results that merit discussion in the 56 57 58 2 59 60 ACS Paragon Plus Environment Page 3 of 6 The Journal of Physical Chemistry Letters 1 2 3 peer-reviewed literature.
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]
  • Supporting Information
    Electronic Supplementary Material (ESI) for RSC Advances. This journal is © The Royal Society of Chemistry 2020 Supporting Information How to Select Ionic Liquids as Extracting Agent Systematically? Special Case Study for Extractive Denitrification Process Shurong Gaoa,b,c,*, Jiaxin Jina,b, Masroor Abroc, Ruozhen Songc, Miao Hed, Xiaochun Chenc,* a State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources, North China Electric Power University, Beijing, 102206, China b Research Center of Engineering Thermophysics, North China Electric Power University, Beijing, 102206, China c Beijing Key Laboratory of Membrane Science and Technology & College of Chemical Engineering, Beijing University of Chemical Technology, Beijing 100029, PR China d Office of Laboratory Safety Administration, Beijing University of Technology, Beijing 100124, China * Corresponding author, Tel./Fax: +86-10-6443-3570, E-mail: [email protected], [email protected] 1 COSMO-RS Computation COSMOtherm allows for simple and efficient processing of large numbers of compounds, i.e., a database of molecular COSMO files; e.g. the COSMObase database. COSMObase is a database of molecular COSMO files available from COSMOlogic GmbH & Co KG. Currently COSMObase consists of over 2000 compounds including a large number of industrial solvents plus a wide variety of common organic compounds. All compounds in COSMObase are indexed by their Chemical Abstracts / Registry Number (CAS/RN), by a trivial name and additionally by their sum formula and molecular weight, allowing a simple identification of the compounds. We obtained the anions and cations of different ILs and the molecular structure of typical N-compounds directly from the COSMObase database in this manuscript.
    [Show full text]
  • An Ab Initio Materials Simulation Code
    CP2K: An ab initio materials simulation code Lianheng Tong Physics on Boat Tutorial , Helsinki, Finland 2015-06-08 Faculty of Natural & Mathematical Sciences Department of Physics Brief Overview • Overview of CP2K - General information about the package - QUICKSTEP: DFT engine • Practical example of using CP2K to generate simulated STM images - Terminal states in AGNR segments with 1H and 2H termination groups What is CP2K? Swiss army knife of molecular simulation www.cp2k.org • Geometry and cell optimisation Energy and • Molecular dynamics (NVE, Force Engine NVT, NPT, Langevin) • STM Images • Sampling energy surfaces (metadynamics) • DFT (LDA, GGA, vdW, • Finding transition states Hybrid) (Nudged Elastic Band) • Quantum Chemistry (MP2, • Path integral molecular RPA) dynamics • Semi-Empirical (DFTB) • Monte Carlo • Classical Force Fields (FIST) • And many more… • Combinations (QM/MM) Development • Freely available, open source, GNU Public License • www.cp2k.org • FORTRAN 95, > 1,000,000 lines of code, very active development (daily commits) • Currently being developed and maintained by community of developers: • Switzerland: Paul Scherrer Institute Switzerland (PSI), Swiss Federal Institute of Technology in Zurich (ETHZ), Universität Zürich (UZH) • USA: IBM Research, Lawrence Livermore National Laboratory (LLNL), Pacific Northwest National Laboratory (PNL) • UK: Edinburgh Parallel Computing Centre (EPCC), King’s College London (KCL), University College London (UCL) • Germany: Ruhr-University Bochum • Others: We welcome contributions from
    [Show full text]
  • Automated Construction of Quantum–Classical Hybrid Models Arxiv:2102.09355V1 [Physics.Chem-Ph] 18 Feb 2021
    Automated construction of quantum{classical hybrid models Christoph Brunken and Markus Reiher∗ ETH Z¨urich, Laboratorium f¨urPhysikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Z¨urich, Switzerland February 18, 2021 Abstract We present a protocol for the fully automated construction of quantum mechanical-(QM){ classical hybrid models by extending our previously reported approach on self-parametri- zing system-focused atomistic models (SFAM) [J. Chem. Theory Comput. 2020, 16 (3), 1646{1665]. In this QM/SFAM approach, the size and composition of the QM region is evaluated in an automated manner based on first principles so that the hybrid model describes the atomic forces in the center of the QM region accurately. This entails the au- tomated construction and evaluation of differently sized QM regions with a bearable com- putational overhead that needs to be paid for automated validation procedures. Applying SFAM for the classical part of the model eliminates any dependence on pre-existing pa- rameters due to its system-focused quantum mechanically derived parametrization. Hence, QM/SFAM is capable of delivering a high fidelity and complete automation. Furthermore, since SFAM parameters are generated for the whole system, our ansatz allows for a con- venient re-definition of the QM region during a molecular exploration. For this purpose, a local re-parametrization scheme is introduced, which efficiently generates additional clas- sical parameters on the fly when new covalent bonds are formed (or broken) and moved to the classical region. arXiv:2102.09355v1 [physics.chem-ph] 18 Feb 2021 ∗Corresponding author; e-mail: [email protected] 1 1 Introduction In contrast to most protocols of computational quantum chemistry that consider isolated molecules, chemical processes can take place in a vast variety of complex environments.
    [Show full text]
  • Popular GPU-Accelerated Applications
    LIFE & MATERIALS SCIENCES GPU-ACCELERATED APPLICATIONS | CATALOG | AUG 12 LIFE & MATERIALS SCIENCES APPLICATIONS CATALOG Application Description Supported Features Expected Multi-GPU Release Status Speed Up* Support Bioinformatics BarraCUDA Sequence mapping software Alignment of short sequencing 6-10x Yes Available now reads Version 0.6.2 CUDASW++ Open source software for Smith-Waterman Parallel search of Smith- 10-50x Yes Available now protein database searches on GPUs Waterman database Version 2.0.8 CUSHAW Parallelized short read aligner Parallel, accurate long read 10x Yes Available now aligner - gapped alignments to Version 1.0.40 large genomes CATALOG GPU-BLAST Local search with fast k-tuple heuristic Protein alignment according to 3-4x Single Only Available now blastp, multi cpu threads Version 2.2.26 GPU-HMMER Parallelized local and global search with Parallel local and global search 60-100x Yes Available now profile Hidden Markov models of Hidden Markov Models Version 2.3.2 mCUDA-MEME Ultrafast scalable motif discovery algorithm Scalable motif discovery 4-10x Yes Available now based on MEME algorithm based on MEME Version 3.0.12 MUMmerGPU A high-throughput DNA sequence alignment Aligns multiple query sequences 3-10x Yes Available now LIFE & MATERIALS& LIFE SCIENCES APPLICATIONS program against reference sequence in Version 2 parallel SeqNFind A GPU Accelerated Sequence Analysis Toolset HW & SW for reference 400x Yes Available now assembly, blast, SW, HMM, de novo assembly UGENE Opensource Smith-Waterman for SSE/CUDA, Fast short
    [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]
  • Starting SCF Calculations by Superposition of Atomic Densities
    Starting SCF Calculations by Superposition of Atomic Densities J. H. VAN LENTHE,1 R. ZWAANS,1 H. J. J. VAN DAM,2 M. F. GUEST2 1Theoretical Chemistry Group (Associated with the Department of Organic Chemistry and Catalysis), Debye Institute, Utrecht University, Padualaan 8, 3584 CH Utrecht, The Netherlands 2CCLRC Daresbury Laboratory, Daresbury WA4 4AD, United Kingdom Received 5 July 2005; Accepted 20 December 2005 DOI 10.1002/jcc.20393 Published online in Wiley InterScience (www.interscience.wiley.com). Abstract: We describe the procedure to start an SCF calculation of the general type from a sum of atomic electron densities, as implemented in GAMESS-UK. Although the procedure is well known for closed-shell calculations and was already suggested when the Direct SCF procedure was proposed, the general procedure is less obvious. For instance, there is no need to converge the corresponding closed-shell Hartree–Fock calculation when dealing with an open-shell species. We describe the various choices and illustrate them with test calculations, showing that the procedure is easier, and on average better, than starting from a converged minimal basis calculation and much better than using a bare nucleus Hamiltonian. © 2006 Wiley Periodicals, Inc. J Comput Chem 27: 926–932, 2006 Key words: SCF calculations; atomic densities Introduction hrstuhl fur Theoretische Chemie, University of Kahrlsruhe, Tur- bomole; http://www.chem-bio.uni-karlsruhe.de/TheoChem/turbo- Any quantum chemical calculation requires properly defined one- mole/),12 GAMESS(US) (Gordon Research Group, GAMESS, electron orbitals. These orbitals are in general determined through http://www.msg.ameslab.gov/GAMESS/GAMESS.html, 2005),13 an iterative Hartree–Fock (HF) or Density Functional (DFT) pro- Spartan (Wavefunction Inc., SPARTAN: http://www.wavefun.
    [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]
  • Kepler Gpus and NVIDIA's Life and Material Science
    LIFE AND MATERIAL SCIENCES Mark Berger; [email protected] Founded 1993 Invented GPU 1999 – Computer Graphics Visual Computing, Supercomputing, Cloud & Mobile Computing NVIDIA - Core Technologies and Brands GPU Mobile Cloud ® ® GeForce Tegra GRID Quadro® , Tesla® Accelerated Computing Multi-core plus Many-cores GPU Accelerator CPU Optimized for Many Optimized for Parallel Tasks Serial Tasks 3-10X+ Comp Thruput 7X Memory Bandwidth 5x Energy Efficiency How GPU Acceleration Works Application Code Compute-Intensive Functions Rest of Sequential 5% of Code CPU Code GPU CPU + GPUs : Two Year Heart Beat 32 Volta Stacked DRAM 16 Maxwell Unified Virtual Memory 8 Kepler Dynamic Parallelism 4 Fermi 2 FP64 DP GFLOPS GFLOPS per DP Watt 1 Tesla 0.5 CUDA 2008 2010 2012 2014 Kepler Features Make GPU Coding Easier Hyper-Q Dynamic Parallelism Speedup Legacy MPI Apps Less Back-Forth, Simpler Code FERMI 1 Work Queue CPU Fermi GPU CPU Kepler GPU KEPLER 32 Concurrent Work Queues Developer Momentum Continues to Grow 100M 430M CUDA –Capable GPUs CUDA-Capable GPUs 150K 1.6M CUDA Downloads CUDA Downloads 1 50 Supercomputer Supercomputers 60 640 University Courses University Courses 4,000 37,000 Academic Papers Academic Papers 2008 2013 Explosive Growth of GPU Accelerated Apps # of Apps Top Scientific Apps 200 61% Increase Molecular AMBER LAMMPS CHARMM NAMD Dynamics GROMACS DL_POLY 150 Quantum QMCPACK Gaussian 40% Increase Quantum Espresso NWChem Chemistry GAMESS-US VASP CAM-SE 100 Climate & COSMO NIM GEOS-5 Weather WRF Chroma GTS 50 Physics Denovo ENZO GTC MILC ANSYS Mechanical ANSYS Fluent 0 CAE MSC Nastran OpenFOAM 2010 2011 2012 SIMULIA Abaqus LS-DYNA Accelerated, In Development NVIDIA GPU Life Science Focus Molecular Dynamics: All codes are available AMBER, CHARMM, DESMOND, DL_POLY, GROMACS, LAMMPS, NAMD Great multi-GPU performance GPU codes: ACEMD, HOOMD-Blue Focus: scaling to large numbers of GPUs Quantum Chemistry: key codes ported or optimizing Active GPU acceleration projects: VASP, NWChem, Gaussian, GAMESS, ABINIT, Quantum Espresso, BigDFT, CP2K, GPAW, etc.
    [Show full text]
  • Dr. David Danovich
    Dr. David Danovich ! 14 January 2021 (h-index: 34; 124 publications) Personal Information Name: Contact information Danovich David Private: Nissan Harpaz st., 4/9, Jerusalem, Date of Birth: 9371643 Israel June 29, 1959 Mobile: +972-54-4768669 E-mail: [email protected] Place of Birth: Irkutsk, USSR Work: Institute of Chemistry, The Hebrew Citizenship: University, Edmond Safra Campus, Israeli Givat Ram, 9190401 Jerusalem, Israel Marital status: Phone: +972-2-6586934 Married +2 FAX: +972-2-6584033 E-mail: [email protected] Current position 1992- present: Senior computational chemist and scientific programmer in the Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel Scientific skills Very experienced scientific researcher with demonstrated ability to implement ab initio molecular electronic structure and valence bond calculations for obtaining insights into physical and chemical properties of molecules and materials. Extensive expertise in implementation of Molecular Orbital Theory based on ab initio methods such as Hartree-Fock theory (HF), Density Functional Theory (DFT), Moller-Plesset Perturbation theory (MP2-MP5), Couple Cluster theories [CCST(T)], Complete Active Space (CAS, CASPT2) theories, Multireference Configuration Interaction (MRCI) theory and modern ab initio valence bond theory (VBT) methods such as Valence Bond Self Consistent Field theory (VBSCF), Breathing Orbital Valence Bond theory (BOVB) and its variants (like D-BOVB, S-BOVB and SD-BOVB), Valence Bond Configuration Interaction theory (VBCI), Valence Bond Perturbation theory (VBPT). Language skills Advanced Conversational English, Russian Hebrew Academic background 1. Undergraduate studies, 1978-1982 Master of Science in Physics, June 1982 Irkutsk State University, USSR Dr. David Danovich ! 2. Graduate studies, 1984-1989 Ph.
    [Show full text]
  • Computational Studies on Carbohydrates: I. Density Functional Ab Initio Geometry Optimization on Maltose Conformations
    Computational Studies on Carbohydrates: I. Density Functional Ab Initio Geometry Optimization on Maltose Conformations F. A. MOMANY, J. L. WILLETT Plant Polymer Research, National Center for Agricultural Utilization Research, USDA, Agricultural Research Service, 1815 N. University St., Peoria, Illinois 61604 Received 10 September 1999; accepted 10 May 2000 ABSTRACT: Ab initio geometry optimization was carried out on 10 selected conformations of maltose and two 2-methoxytetrahydropyran conformations using the density functional denoted B3LYP combined with two basis sets. The 6-31G∗ and 6-311CCG∗∗ basis sets make up the B3LYP/6-31G∗ and B3LYP/6-311CCG∗∗ procedures. Internal coordinates were fully relaxed, and structures were gradient optimized at both levels of theory. Ten conformations were studied at the B3LYP/6-31G∗ level, and five of these were continued with ∗∗ full gradient optimization at the B3LYP/6-311CCG level of theory. The details of the ab initio optimized geometries are presented here, with particular attention given to the positions of the atoms around the anomeric center and the effect of the particular anomer and hydrogen bonding pattern on the maltose ring structures and relative conformational energies. The size and complexity of the hydrogen-bonding network prevented a rigorous search of conformational space by ab initio calculations. However, using empirical force fields, low-energy conformers of maltose were found that were subsequently gradient optimized at the two ab initio levels of theory. Three classes of conformations were studied, as defined by the clockwise or counterclockwise direction of the hydroxyl groups, or a flipped conformer in which the -dihedral is rotated by ∼180◦.Different combinations of ! side-chain rotations gave energy differences of more than 6 kcal/mol above the lowest energy structure found.
    [Show full text]