(PDF) DFT Software: Who, What, and How Much

Total Page:16

File Type:pdf, Size:1020Kb

(PDF) DFT Software: Who, What, and How Much DFT Software: Who, What, and How Much Brad Hubartt Introduction: was constructed [11]. KS-DFT takes the Doing pretty much anything in original energy functional: quantum mechanics is going to involve the E[ρ]=T [ρ]+ Vext(r)ρ(r)dr + VH [ρ]+Exc[ρ] wave equation of the system. Using the wave ! equation in a large system is going to be very where T is the kinetic energy, Vext is the difficult or impossible, if the wave equation is external potential, VH is the Hartree energy, even known. This essentially rules out using Exc is the exchange correlation energy, and ρ quantum mechanics in biophysics/chemistry. is the charge density. ρ is what Kohn and A new quantum model is needed for these Sham modified. They derived the following systems. Density Functional Theory (DFT) three equations: [4][6][9][12] is such a model. It involves N using the probability density functions of the ρ(r)= ψ (r, s) 2 electrons in a molecule to determine other | i | i s 2 ! ! properties of the molecule. Instead of solving ! 2ψ (r, s)+v (r)ψ (r, s)=# ψ (r, s) problems with wave equations for each − 2m∇ i eff i i i electron it uses the probabilities of electrons 2 ρ(r!) δExc[ρ] veff (r)=Vext(r)+e dr! + being at different locations. If the wave r r! δρ function is not known many implementations " | − | of DFT will allow for Gaussians to be used in The first one is exactly what anyone would their place. expect, the second is Schrödinger’s equation, DFT works by minimizing functions and the third is the actual potential at a given of functions, functionals, to determine the location based upon the externally applied ground state of a molecule. Instead of using potential, the charge density, and the functionals based upon each individual correlation energy. ρ can be solved for by electron, a pseudopotential is used to simplify integrating over these three equations with an the calculation. Additionally it can calculate RK4 [10] or something equivalent. Often the forces and has been expanded to include time- wave equation is replaced with a Gaussian to dependence [13]. Time-dependent DFT simplify things. makes simulations of molecular dynamics Another technique used to simplify possible. Molecular dynamics are of DFT is replacing all-electron systems with a particular interest as nothing in a biological pseudopotential. A pseudopotential is an system is static; proteins are floating around approximation of the Coulombic potential in a solution and constantly interacting with from the non-valence electrons. In chemistry different ions. the valence electrons are the ones that actually form bonds and cause reactions; the Methods: non-valence electrons are just there for the This is a seemingly wonderful theory, potential field (obviously a gross but how is it implemented? It is not oversimplification). Pseudopotentials are used particularly useful if every equation has to be when generating a self-consistent field. Pulay solved by hand at each iteration. First, where Mixing [3] (also called Direct Inversion of does the electron density come from? Kohn Iterative Subspace) is a non-linear mixing and Sham came up with an approach for this technique used for large SCF problems. which is how Kohn-Sham DFT (KS-DFT) Hubartt 1 Instead of having to store a giant hessian for unit cell, nband is the number of bands to be the system this algorithm works with an considered, and iscf determines which mixing amount n (user specified integer) of charge algorithm is used to find a self-consistent density vectors. Through a few linear field. The iscf variable by default is set to use operations a final charge density vector is Pulay Mixing. Other mixing algorithms are computed. The computational benefits of built into Abinit that are simpler than Pulay Pulay Mixing are increasingly noticeable the which could give mixed results based upon large the system becomes as the hessian the size and complexity of the system. Other grows as nelectrons squared versus needing only more advanced input variables are built into a handful of vectors for Pulay mixing. The Abinit including the Bravais lattice type, how actual algorithm for Pulay mixing is in to calculate the Berry phase, units for energy, reference [3] and [4]. various dielectric properties, etc... ABINIT No one wants to solve this by hand, has two different pseudopotentials built into it the natural place for this theory is in computer for each element. Additionally, it has multiple simulations. A quick internet search reveals a algorithms to calculate a new pseudopotential few dozen pieces of software for solving if one is desired. these systems. Each one comes with its own What does it output? Besides a quirks, learning curves, availability of confirmation of all the input parameters it has documentation, open/closed-source, and a hefty amount of items to output. The first is whether the originator of the software has the point symmetry group and the Bravais been banned for life from using it (amongst lattice type. The Bravais lattice, along with other people). The purpose of this paper is to many of the input parameters, are geared select a few of these programs to discuss how more for solid state physics than molecular they work. Featured programs were picked physics/chemistry. It has to do with the fairly randomly from a list of available ones crystal structure of the material being with little bias save one. analyzed. The point symmetry group has to do with how many rotational symmetries Abinit [2]: exist, again this is geared more for crystal This software was developed using a structures. After this comes the real output. It popular formulation of DFT known as KS- gives the atomic positions in cartesian DFT (Kohn-Sham). It is currently run by coordinates, energy gradients, forces in eV/ Xavier Gonze and written in FORTRAN. Angstrom, eigenvalues of each band for each Dozens of variables are available for k point, minimum and maximum charge customizing the input conditions of a system. densities, a list of all the different energy A few of the basic variables are: acell, components, and the stress (average force per angdeg, ecut, natom, nband, and iscf; area) tensor in GPa. essentially every property that can change from one system to the next is customizable if Gaussian: desired. Keep in mind that most of these have This company’s program has had quite fairly generic default settings. acell gives the the history. It was originally developed by size of the unit cell usually in Borh radii noble laureate John Pople in 1971 at Carnegie (other units are allowed if you specify them), Mellon in FORTRAN. In the beginning the angdeg sets the angles between the three source was open and freely distributed to any primitive vectors for the unit cell, ecut sets a university that wanted to use it. This is far maximum value for the kinetic energy of the from the present case. The program is now system, natom is the number of atoms in a somewhat closed-source and sold by Hubartt 2 Gaussian Inc. under a fairly restrictive forces calculated at a given time-step with license. Various stipulations of the license are: DFT and then integrated via a Verlet not being allowed to disclose how fast the algorithm. The timescale for these dynamics program is, work on any “competing” is small, on order of picoseconds, but it is a program, or allow anyone that is banned from start. using Gaussian to use your copy. Who exactly For visualizations the output files can would be banned from using this? Pople (the be imported to a variety of additional creator!) is banned from using it, as are programs. These allow for 3D overlays of the Princeton and CalTech along with many other electron densities with molecular structures. universities and individuals. The official list One such program, gOpenMol [8] produced of bannings is kept at Fig 1. www.bannedbygaussian.org. As it is illegal to discuss how fast this program is, any judgement as to whether Gaussian is a good program that is worth the $4500-36,000 (academic and industry prices for OS X version) for a site license is difficult to determine. It is especially hard to justify this price when there are open-source programs that are freely available. This company has established a huge barrier of entry for anyone wishing to get into the field. With the risk of being banned from using this software a clear Fig. 1 Electron density of the aspirin molecule danger, it is hard to imagine anyone using this software. Despite this, many groups do use Setting up the input for QE is not Gaussian for their DFT calculations[7]. particularly difficult. The user specifies the element(s) being used, which pseudopotential Quantum Espresso (QE) [15]: to use, the Bravais lattice, the k-point QE is another open-source, free locations of the atoms in the 1st Brillouin software suite for DFT. It is operated through zone (locations of the atoms in the inverse a collaboration between MIT, Princeton lattice), and other characteristics of their (probably why they are banned from system. Many common lattices can be found Gaussian), and multiple European on their website or on a few other databases Universities; it is coordinated by Paolo that they refer the user to. Giannozzi. It is written in a combination of C++ and FORTRAN. C++ is used to make Current Biophysics Uses: the program object-oriented, very modular, The purpose of this paper was to and easier to learn initially with bit of describe a few of the many choices in FORTRAN for speed boosts.
Recommended publications
  • GPAW, Gpus, and LUMI
    GPAW, GPUs, and LUMI Martti Louhivuori, CSC - IT Center for Science Jussi Enkovaara GPAW 2021: Users and Developers Meeting, 2021-06-01 Outline LUMI supercomputer Brief history of GPAW with GPUs GPUs and DFT Current status Roadmap LUMI - EuroHPC system of the North Pre-exascale system with AMD CPUs and GPUs ~ 550 Pflop/s performance Half of the resources dedicated to consortium members Programming for LUMI Finland, Belgium, Czechia, MPI between nodes / GPUs Denmark, Estonia, Iceland, HIP and OpenMP for GPUs Norway, Poland, Sweden, and how to use Python with AMD Switzerland GPUs? https://www.lumi-supercomputer.eu GPAW and GPUs: history (1/2) Early proof-of-concept implementation for NVIDIA GPUs in 2012 ground state DFT and real-time TD-DFT with finite-difference basis separate version for RPA with plane-waves Hakala et al. in "Electronic Structure Calculations on Graphics Processing Units", Wiley (2016), https://doi.org/10.1002/9781118670712 PyCUDA, cuBLAS, cuFFT, custom CUDA kernels Promising performance with factor of 4-8 speedup in best cases (CPU node vs. GPU node) GPAW and GPUs: history (2/2) Code base diverged from the main branch quite a bit proof-of-concept implementation had lots of quick and dirty hacks fixes and features were pulled from other branches and patches no proper unit tests for GPU functionality active development stopped soon after publications Before development re-started, code didn't even work anymore on modern GPUs without applying a few small patches Lesson learned: try to always get new functionality to the
    [Show full text]
  • Free and Open Source Software for Computational Chemistry Education
    Free and Open Source Software for Computational Chemistry Education Susi Lehtola∗,y and Antti J. Karttunenz yMolecular Sciences Software Institute, Blacksburg, Virginia 24061, United States zDepartment of Chemistry and Materials Science, Aalto University, Espoo, Finland E-mail: susi.lehtola@alumni.helsinki.fi Abstract Long in the making, computational chemistry for the masses [J. Chem. Educ. 1996, 73, 104] is finally here. We point out the existence of a variety of free and open source software (FOSS) packages for computational chemistry that offer a wide range of functionality all the way from approximate semiempirical calculations with tight- binding density functional theory to sophisticated ab initio wave function methods such as coupled-cluster theory, both for molecular and for solid-state systems. By their very definition, FOSS packages allow usage for whatever purpose by anyone, meaning they can also be used in industrial applications without limitation. Also, FOSS software has no limitations to redistribution in source or binary form, allowing their easy distribution and installation by third parties. Many FOSS scientific software packages are available as part of popular Linux distributions, and other package managers such as pip and conda. Combined with the remarkable increase in the power of personal devices—which rival that of the fastest supercomputers in the world of the 1990s—a decentralized model for teaching computational chemistry is now possible, enabling students to perform reasonable modeling on their own computing devices, in the bring your own device 1 (BYOD) scheme. In addition to the programs’ use for various applications, open access to the programs’ source code also enables comprehensive teaching strategies, as actual algorithms’ implementations can be used in teaching.
    [Show full text]
  • D6.1 Report on the Deployment of the Max Demonstrators and Feedback to WP1-5
    Ref. Ares(2020)2820381 - 31/05/2020 HORIZON2020 European Centre of Excellence ​ Deliverable D6.1 Report on the deployment of the MaX Demonstrators and feedback to WP1-5 D6.1 Report on the deployment of the MaX Demonstrators and feedback to WP1-5 Pablo Ordejón, Uliana Alekseeva, Stefano Baroni, Riccardo Bertossa, Miki Bonacci, Pietro Bonfà, Claudia Cardoso, Carlo Cavazzoni, Vladimir Dikan, Stefano de Gironcoli, Andrea Ferretti, Alberto García, Luigi Genovese, Federico Grasselli, Anton Kozhevnikov, Deborah Prezzi, Davide Sangalli, Joost VandeVondele, Daniele Varsano, Daniel Wortmann Due date of deliverable: 31/05/2020 Actual submission date: 31/05/2020 Final version: 31/05/2020 Lead beneficiary: ICN2 (participant number 3) Dissemination level: PU - Public www.max-centre.eu 1 HORIZON2020 European Centre of Excellence ​ Deliverable D6.1 Report on the deployment of the MaX Demonstrators and feedback to WP1-5 Document information Project acronym: MaX Project full title: Materials Design at the Exascale Research Action Project type: European Centre of Excellence in materials modelling, simulations and design EC Grant agreement no.: 824143 Project starting / end date: 01/12/2018 (month 1) / 30/11/2021 (month 36) Website: www.max-centre.eu Deliverable No.: D6.1 Authors: P. Ordejón, U. Alekseeva, S. Baroni, R. Bertossa, M. ​ Bonacci, P. Bonfà, C. Cardoso, C. Cavazzoni, V. Dikan, S. de Gironcoli, A. Ferretti, A. García, L. Genovese, F. Grasselli, A. Kozhevnikov, D. Prezzi, D. Sangalli, J. VandeVondele, D. Varsano, D. Wortmann To be cited as: Ordejón, et al., (2020): Report on the deployment of the MaX Demonstrators and feedback to WP1-5. Deliverable D6.1 of the H2020 project MaX (final version as of 31/05/2020).
    [Show full text]
  • The Impact of Density Functional Theory on Materials Research
    www.mrs.org/bulletin functional. This functional (i.e., a function whose argument is another function) de- scribes the complex kinetic and energetic interactions of an electron with other elec- Toward Computational trons. Although the form of this functional that would make the reformulation of the many-body Schrödinger equation exact is Materials Design: unknown, approximate functionals have proven highly successful in describing many material properties. Efficient algorithms devised for solving The Impact of the Kohn–Sham equations have been imple- mented in increasingly sophisticated codes, tremendously boosting the application of DFT methods. New doors are opening to in- Density Functional novative research on materials across phys- ics, chemistry, materials science, surface science, and nanotechnology, and extend- ing even to earth sciences and molecular Theory on Materials biology. The impact of this fascinating de- velopment has not been restricted to aca- demia, as DFT techniques also find application in many different areas of in- Research dustrial research. The development is so fast that many current applications could Jürgen Hafner, Christopher Wolverton, and not have been realized three years ago and were hardly dreamed of five years ago. Gerbrand Ceder, Guest Editors The articles collected in this issue of MRS Bulletin present a few of these suc- cess stories. However, even if the compu- Abstract tational tools necessary for performing The development of modern materials science has led to a growing need to complex quantum-mechanical calcula- tions relevant to real materials problems understand the phenomena determining the properties of materials and processes on are now readily available, designing a an atomistic level.
    [Show full text]
  • Practice: Quantum ESPRESSO I
    MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW FREE http://www.quantum-espresso.org PW-DFT, PP, PAW FREE http://www.abinit.org DFT PW, PP, Car-Parrinello FREE http://www.cpmd.org DFT PP, US-PP, PAW $3000 [moderate accuracy, fast] http://www.vasp.at DFT full-potential linearized augmented $500 plane-wave (FLAPW) [accurate, slow] http://www.wien2k.at Hartree-Fock, higher order correlated $3000 electron approaches http://www.gaussian.com 3 Quantum ESPRESSO 4 Quantum ESPRESSO Quantum ESPRESSO is an integrated suite of Open- Source computer codes for electronic-structure calculations and materials modeling at the nanoscale. It is based on density-functional theory, plane waves, and pseudopotentials. Core set of codes, plugins for more advanced tasks and third party packages Open initiative coordinated by the Quantum ESPRESSO Foundation, across Italy. Contributed to by developers across the world Regular hands-on workshops in Trieste, Italy Open-source code: FREE (unlike VASP...) 5 Performance Small jobs (a few atoms) can be run on single node Includes determining convergence parameters, lattice constants Can use OpenMP parallelization on multicore machines Large jobs (~10’s to ~100’s atoms) can run in parallel using MPI to 1000’s of cores Includes molecular dynamics, large geometry relaxation, phonons Parallel performance tied to BLAS/LAPACK (linear algebra routines) and 3D FFT (fast Fourier transform) New GPU-enabled version available 6 Usability Documented online:
    [Show full text]
  • The CECAM Electronic Structure Library and the Modular Software Development Paradigm
    The CECAM electronic structure library and the modular software development paradigm Cite as: J. Chem. Phys. 153, 024117 (2020); https://doi.org/10.1063/5.0012901 Submitted: 06 May 2020 . Accepted: 08 June 2020 . Published Online: 13 July 2020 Micael J. T. Oliveira , Nick Papior , Yann Pouillon , Volker Blum , Emilio Artacho , Damien Caliste , Fabiano Corsetti , Stefano de Gironcoli , Alin M. Elena , Alberto García , Víctor M. García-Suárez , Luigi Genovese , William P. Huhn , Georg Huhs , Sebastian Kokott , Emine Küçükbenli , Ask H. Larsen , Alfio Lazzaro , Irina V. Lebedeva , Yingzhou Li , David López- Durán , Pablo López-Tarifa , Martin Lüders , Miguel A. L. Marques , Jan Minar , Stephan Mohr , Arash A. Mostofi , Alan O’Cais , Mike C. Payne, Thomas Ruh, Daniel G. A. Smith , José M. Soler , David A. Strubbe , Nicolas Tancogne-Dejean , Dominic Tildesley, Marc Torrent , and Victor Wen-zhe Yu COLLECTIONS Paper published as part of the special topic on Electronic Structure Software Note: This article is part of the JCP Special Topic on Electronic Structure Software. This paper was selected as Featured ARTICLES YOU MAY BE INTERESTED IN Recent developments in the PySCF program package The Journal of Chemical Physics 153, 024109 (2020); https://doi.org/10.1063/5.0006074 An open-source coding paradigm for electronic structure calculations Scilight 2020, 291101 (2020); https://doi.org/10.1063/10.0001593 Siesta: Recent developments and applications The Journal of Chemical Physics 152, 204108 (2020); https://doi.org/10.1063/5.0005077 J. Chem. Phys. 153, 024117 (2020); https://doi.org/10.1063/5.0012901 153, 024117 © 2020 Author(s). The Journal ARTICLE of Chemical Physics scitation.org/journal/jcp The CECAM electronic structure library and the modular software development paradigm Cite as: J.
    [Show full text]
  • Kepler Gpus and NVIDIA's Life and Material Science
    LIFE AND MATERIAL SCIENCES Mark Berger; mberger@nvidia.com 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]
  • Density Functional Theory and Nuclear Quantum Effects
    1 Grassmann manifold, gauge, and quantum chemistry Lin Lin Department of Mathematics, UC Berkeley; Lawrence Berkeley National Laboratory Flatiron Institute, April 2019 2 Stiefel manifold • × > ℂ • Stiefel manifold (1935): orthogonal vectors in × ℂ , = | = ∗ ∈ ℂ 3 Grassmann manifold • Grassmann manifold (1848): set of dimensional subspace in ℂ , = , / : dimensional unitary group (i.e. the set of unitary matrices in × ) ℂ • Any point in , can be viewed as a representation of a point in , 4 Example (in optimization) min ( ) . = ∗ • = , . invariant to the choice of basis Optimization on a Grassmann manifold. ∈ ⇒ • The representation is called the gauge in quantum physics / chemistry. Gauge invariant. 5 Example (in optimization) min ( ) . = ∗ • Otherwise, if ( ) is gauge dependent Optimization on a Stiefel manifold. ⇒ • [Edelman, Arias, Smith, SIMAX 1998] many works following 6 This talk: Quantum chemistry • Time-dependent density functional theory (TDDFT) • Kohn-Sham density functional theory (KSDFT) • Localization. • Gauge choice is the key • Used in community software packages: Quantum ESPRESSO, Wannier90, Octopus, PWDFT etc 7 Time-dependent density functional theory Parallel transport gauge 8 Time-dependent density functional theory [Runge-Gross, PRL 1984] (spin omitted, : number of electrons) , = , , , = 1, … , , , = ( ,) ( , ) ′ ∗ ′ � Hamiltonian =1 1 , = + ( ) + [ ] 2 − Δ 9 Density matrix • Key quantity throughout the talk × • (t) = , … , . : number of discretized grid points to represent each Ψ 1 ∈ ℂ = ∗ Ψ Ψ = 10 Gauge invariance • Unitary rotation = , = ∗ = Φ Ψ= ( ) = ∗ ∗ ∗ ∗ ΦΦ Ψ Ψ ΨΨ • Propagation of the density matrix, von Neumann equation (quantum Liouville equation) = , = , ( ) , − 11 Gauge choice affects the time step Extreme case (boring Hamiltonian) ≡ 0 Initial condition 0 = 0 0 Exact solution = (0) Small time step − = (0) (Formally) arbitrarily large time step 12 P v.s.
    [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]
  • CPMD Tutorial Car–Parrinello Molecular Dynamics
    CPMD Tutorial Car–Parrinello Molecular Dynamics Ari P Seitsonen CNRS & Universit´ePierre et Marie Curie, Paris CSC, October 2006 Introduction CPMD — concept 1 • Car-Parrinello Molecular Dynamics – Roberto Car & Michele Parrinello, Physical Review Letters 55, 2471 (1985) — 20+1 years of Car-Parrinello method! Roberto Car Michele Parrinello CPMD — concept 2 • CPMD code — http://www.cpmd.org/; based on the original code of Roberto Car and Michele Parrinello CPMD What is it, what is it not? • Computer code for performing static, dynamic simulations and analysis of electronic structure • About 200’000 lines of code — FORTRAN77 with a flavour towards F90 • (Freely) available via http://www.cpmd.org/ with source code • Most suitable for dynamical simulations of condensed systems or large molecules — less for small molecules or bulk properties of small crystals • Computationally highly optimised for vector and scalar supercomputers, shared and distributed memory machines and combinations thereof (SMP) Capabilities • Solution of electronic and ionic structure • XC-functional: LDA, GGA, meta-GGA, hybrid • Molecular dynamics in NVE, NVT, NPT ensembles • General constraints (MD and geometry optimisation) • Metadynamics ‡ • Free energy functional ‡ • Path integrals ‡ • QM/MM ‡ • Wannier functions • Response properties: TDDFT, NMR, Raman, IR, . • Norm conserving and ultra-soft pseudo potentials ‡ : Will not be considered during this course CPMD: Examples Amorphous Si CPMD: Examples Amorphous Si CPMD: Examples Liquid water: Solvation and transport of
    [Show full text]
  • Pyscf: the Python-Based Simulations of Chemistry Framework
    PySCF: The Python-based Simulations of Chemistry Framework Qiming Sun∗1, Timothy C. Berkelbach2, Nick S. Blunt3,4, George H. Booth5, Sheng Guo1,6, Zhendong Li1, Junzi Liu7, James D. McClain1,6, Elvira R. Sayfutyarova1,6, Sandeep Sharma8, Sebastian Wouters9, and Garnet Kin-Lic Chany1 1Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena CA 91125, USA 2Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA 3Chemical Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 4Department of Chemistry, University of California, Berkeley, California 94720, USA 5Department of Physics, King's College London, Strand, London WC2R 2LS, United Kingdom 6Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA 7Institute of Chemistry Chinese Academy of Sciences, Beijing 100190, P. R. China 8Department of Chemistry and Biochemistry, University of Colorado Boulder, Boulder, CO 80302, USA 9Brantsandpatents, Pauline van Pottelsberghelaan 24, 9051 Sint-Denijs-Westrem, Belgium arXiv:1701.08223v2 [physics.chem-ph] 2 Mar 2017 Abstract PySCF is a general-purpose electronic structure platform designed from the ground up to emphasize code simplicity, so as to facilitate new method development and enable flexible ∗osirpt.sun@gmail.com ygkc1000@gmail.com 1 computational workflows. The package provides a wide range of tools to support simulations of finite-size systems, extended systems with periodic boundary conditions, low-dimensional periodic systems, and custom Hamiltonians, using mean-field and post-mean-field methods with standard Gaussian basis functions. To ensure ease of extensibility, PySCF uses the Python language to implement almost all of its features, while computationally critical paths are implemented with heavily optimized C routines.
    [Show full text]
  • Arxiv:2005.05756V2
    “This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article appeared in Oliveira, M.J.T. [et al.]. The CECAM electronic structure library and the modular software development paradigm. "The Journal of Chemical Physics", 12020, vol. 153, núm. 2, and may be found at https://aip.scitation.org/doi/10.1063/5.0012901. The CECAM Electronic Structure Library and the modular software development paradigm Micael J. T. Oliveira,1, a) Nick Papior,2, b) Yann Pouillon,3, 4, c) Volker Blum,5, 6 Emilio Artacho,7, 8, 9 Damien Caliste,10 Fabiano Corsetti,11, 12 Stefano de Gironcoli,13 Alin M. Elena,14 Alberto Garc´ıa,15 V´ıctor M. Garc´ıa-Su´arez,16 Luigi Genovese,10 William P. Huhn,5 Georg Huhs,17 Sebastian Kokott,18 Emine K¨u¸c¨ukbenli,13, 19 Ask H. Larsen,20, 4 Alfio Lazzaro,21 Irina V. Lebedeva,22 Yingzhou Li,23 David L´opez-Dur´an,22 Pablo L´opez-Tarifa,24 Martin L¨uders,1, 14 Miguel A. L. Marques,25 Jan Minar,26 Stephan Mohr,17 Arash A. Mostofi,11 Alan O'Cais,27 Mike C. Payne,9 Thomas Ruh,28 Daniel G. A. Smith,29 Jos´eM. Soler,30 David A. Strubbe,31 Nicolas Tancogne-Dejean,1 Dominic Tildesley,32 Marc Torrent,33, 34 and Victor Wen-zhe Yu5 1)Max Planck Institute for the Structure and Dynamics of Matter, D-22761 Hamburg, Germany 2)DTU Computing Center, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark 3)Departamento CITIMAC, Universidad de Cantabria, Santander, Spain 4)Simune Atomistics, 20018 San Sebasti´an,Spain 5)Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC 27708, USA 6)Department of Chemistry, Duke University, Durham, NC 27708, USA 7)CIC Nanogune BRTA and DIPC, 20018 San Sebasti´an,Spain 8)Ikerbasque, Basque Foundation for Science, 48011 Bilbao, Spain 9)Theory of Condensed Matter, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, United Kingdom 10)Department of Physics, IRIG, Univ.
    [Show full text]