Toward Quantum-Chemical Method Development for Arbitrary Basis Functions Michael F

Total Page:16

File Type:pdf, Size:1020Kb

Toward Quantum-Chemical Method Development for Arbitrary Basis Functions Michael F THE JOURNAL OF CHEMICAL PHYSICS 149, 084106 (2018) Toward quantum-chemical method development for arbitrary basis functions Michael F. Herbst,1,a) Andreas Dreuw,1,b) and James Emil Avery2,c) 1Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany 2Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 København, Denmark (Received 15 June 2018; accepted 2 August 2018; published online 27 August 2018) Wepresent the design of a flexible quantum-chemical method development framework, which supports employing any type of basis function. This design has been implemented in the light-weight program package molsturm, yielding a basis-function-independent self-consistent field scheme. Versatile interfaces, making use of open standards like python, mediate the integration of molsturm with existing third-party packages. In this way, both rapid extension of the present set of meth- ods for electronic structure calculations as well as adding new basis function types can be readily achieved. This makes molsturm well-suitable for testing novel approaches for discretising the electronic wave function and allows comparing them to existing methods using the same software stack. This is illustrated by two examples, an implementation of coupled-cluster doubles as well as a gradient-free geometry optimisation, where in both cases, arbitrary basis functions could be used. molsturm is open-sourced and can be obtained from http://molsturm.org. Published by AIP Publishing. https://doi.org/10.1063/1.5044765 I. INTRODUCTION computation per integral for having fewer more accurate basis functions. The central goal of electronic-structure theory is to find In many practical applications, the shortcomings of GTOs approximate solutions to the electronic wave equation numer- are not important, or one is able to compensate by employing ically. This requires a discretisation of the electronic wave specialised contracted basis sets.4,5 However, even contracted function. Typically, it is approximated as a linear combi- GTOs (cGTO) fail to describe both the nuclear cusp and the nation of Slater determinants: anti-symmetrised products of exponential decay of the electron density.6 In addition, no mat- single-particle functions. The latter are in turn constructed by ter the number of GTO basis functions used, the derivatives are expanding them in a basis set of a priori determined single- always wrong at the nucleus, which causes singularities and electron functions. Usually, such basis sets are not complete computational failure, for example, in quantum Monte Carlo and introduce basis set errors. Proper choice of the basis func- calculations.7–9 Furthermore the description of some proper- tion type and size of basis set is thus decisive for an accurate ties such as the nuclear-magnetic resonance (NMR) shield- description of the system under investigation. It is also clear ing tensors or a description of Rydberg-like or auto-ionising from the onset that different basis function types can be more states10–13 directly involves the nuclear cusp or the asymptotic or less suited for a specific problem, suggesting us to conduct tail, making physically accurate basis functions desirable.14,15 investigations across existing basis function types. The name of our implementation, “molsturm,” is a port- Gaussian-based methods are overwhelmingly predomi- manteau of molecular Sturmians: the project was born as a nant in computational electronic structure theory, which stems means to solve the problem of using state-of-the-art quantum from pragmatic reasons dating back to the founding years.1,2 chemistry methods together with generalised- and molecular It was well-known that bound state electronic wave func- Sturmian basis functions. The many promising results for gen- tions decay exponentially both at short and large distances eralised Sturmians were stranded due to the fact that only from the nuclei,3 but multi-centre electron-repulsion integrals electronic structure problems that were small enough to be (ERI) of products of exponential-type orbitals (ETO) were solved by direct configuration interaction methods could be impractically difficult to calculate. For Gaussian-type orbitals treated, preventing wider use. The existing mature quantum (GTO), on the other hand, ERI could be calculated efficiently chemistry software has been developed over hundreds of man- due to the Gaussian product theorem. However, the compu- years, and the methods are not easily reimplemented from tational challenges facing quantum chemists have changed scratch. since the 1970’s, and it may now be worth trading extra In theory, it should be a simple matter to include new basis function types in existing quantum chemistry software by swapping the integral calculator. In practice, it turned out to be exceedingly difficult due to the very large and com- a)Electronic mail: [email protected] b)Electronic mail: [email protected] plicated code bases of all the investigated quantum chem- c)Electronic mail: [email protected] istry software. Assumptions about the basis function type 0021-9606/2018/149(8)/084106/14/$30.00 149, 084106-1 Published by AIP Publishing. 084106-2 Herbst, Dreuw, and Avery J. Chem. Phys. 149, 084106 (2018) scattered around the source code only make this even task more back end in molsturm, which will then both provide simple difficult. stand-alone calculations and a common interface to hook into Our solution, which is presented in this paper, is to imple- existing quantum chemistry packages. ment a light-weight layer that makes it easy to experiment with many different basis function types and quantum-chemical B. Toward basis-type agnostic quantum chemistry methods. It is designed for researchers to both build simple stand-alone programs for prototyping and teaching purposes, In order to reach a basis-type agnostic design, there and make plug-in modules to be hosted in standard quan- are three fundamental components to consider: (i) an inte- tum chemistry software. To the best of our knowledge, such a gral interface accommodating a wide range of very dif- framework with the ability to explore quantum-chemical meth- ferent basis set types and discretisation, but providing a ods across multiple basis function types has been missing up uniform way of accessing them, (ii) simple discretisation- to today. agnostic implementations of the self-consistent field (SCF) algorithms, and (iii) a flexible interface to employ the result- A. Alternative basis function types ing SCF orbital basis further in existing third-party code. Once the SCF orbitals have been obtained, the remain- Many research groups have worked on alternative basis der of a calculation, e.g., a Post-Hartree-Fock (Post-HF) function types. This section will provide a brief overview with method, can usually be formulated entirely in the SCF orbital particular focus on exponential-type orbitals (ETO). For fur- basis, without reference to the underlying basis functions. ther details regarding the basis function mentioned, the reader Thus, a basis-function independent SCF scheme automati- is referred to the cited works. cally leads to basis-function independent Post-HF methods as Efforts on making various types of ETOs computation- well. 16–19 ally feasible were pioneered by Harris et al. A particular This structure has another advantageous side effect in the 20–24 form of complete ETO basis is Coulomb Sturmians (CS). context of developing new basis function types, as it allows us Their functional form is identical to the familiar hydrogen-like to perform comparisons between old and new methods using orbitals, just with all occurrences of the factors Z/r replaced exactly the same software stack. In other words, one can thus by the Sturmian exponent k—a parameter, which is the same be sure that, apart from the discretisation, all aspects of the cal- for all functions of the basis, culation, e.g., SCF algorithms or guess methods, are optimised CS 3=2 l −kr 2l+1 m at the same level leading to a fair apples-to-apples comparison ' r = k Nnl(2kr) e L (2kr)Y rˆ . (1) nlm n−l−1 l between old and new methods. m 2l+1 In this, Yl is a spherical harmonic, Ln−l−1 an associated Laguerre polynomial, and C. Paper outline s 2 (l + n)! The remainder of the paper is structured as follows: Sec.II N = (2) nl (2l + 1)! n(n − l − 1)! reviews existing projects with similar goals to molsturm. Section III provides a theoretical background for the program a normalisation constant. CS functions proved to be especially design choices, which are described in Sec.IV. SectionV easy to work with since their momentum-space representation provides example problems calculated using molsturm’s by hyperspherical harmonics allows efficient calculation of python interface, illustrating how to implement new methods multi-centre integrals, opening the way for efficient molecular on top of molsturm in a few lines of python. SectionVI calculations.25–28 The Coulomb Sturmian construction gen- outlines the current state of molsturm and what we hope to eralises well and generalised Sturmian basis sets preserving achieve in the future. many useful Sturmian properties can be constructed. These allow us, for example, to build N-particle basis functions that include important geometric properties of the physical II. RELATED QUANTUM-CHEMICAL SOFTWARE PACKAGES system under consideration at the level of the basis.29–37 Sim- ilarly, d-dimensional hyperspherical harmonic basis sets can This section reviews existing quantum chemistry software model collective motions of particles, for example, for treat- that shares some of the goals of molsturm. ing strongly interacting few-body systems or reactive scatter- The quantum Monte Carlo packages CASINO61 and ing.38–42 A particular type of one-particle Sturmians combines QMCPACK62 are among the few systems that support many a bound-state region and plane-wave asymptotics to model different basis function types.
Recommended publications
  • 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: [email protected].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]
  • User Manual for the Uppsala Quantum Chemistry Package UQUANTCHEM V.35
    User manual for the Uppsala Quantum Chemistry package UQUANTCHEM V.35 by Petros Souvatzis Uppsala 2016 Contents 1 Introduction 6 2 Compiling the code 7 3 What Can be done with UQUANTCHEM 9 3.1 Hartree Fock Calculations . 9 3.2 Configurational Interaction Calculations . 9 3.3 M¨ollerPlesset Calculations (MP2) . 9 3.4 Density Functional Theory Calculations (DFT)) . 9 3.5 Time Dependent Density Functional Theory Calculations (TDDFT)) . 10 3.6 Quantum Montecarlo Calculations . 10 3.7 Born Oppenheimer Molecular Dynamics . 10 4 Setting up a UQANTCHEM calculation 12 4.1 The input files . 12 4.1.1 The INPUTFILE-file . 12 4.1.2 The BASISFILE-file . 13 4.1.3 The BASISFILEAUX-file . 14 4.1.4 The DENSMATSTARTGUESS.dat-file . 14 4.1.5 The MOLDYNRESTART.dat-file . 14 4.1.6 The INITVELO.dat-file . 15 4.1.7 Running Uquantchem . 15 4.2 Input parameters . 15 4.2.1 CORRLEVEL ................................. 15 4.2.2 ADEF ..................................... 15 4.2.3 DOTDFT ................................... 16 4.2.4 NSCCORR ................................... 16 4.2.5 SCERR .................................... 16 4.2.6 MIXTDDFT .................................. 16 4.2.7 EPROFILE .................................. 16 4.2.8 DOABSSPECTRUM ............................... 17 4.2.9 NEPERIOD .................................. 17 4.2.10 EFIELDMAX ................................. 17 4.2.11 EDIR ..................................... 17 4.2.12 FIELDDIR .................................. 18 4.2.13 OMEGA .................................... 18 2 CONTENTS 3 4.2.14 NCHEBGAUSS ................................. 18 4.2.15 RIAPPROX .................................. 18 4.2.16 LIMPRECALC (Only openmp-version) . 19 4.2.17 DIAGDG ................................... 19 4.2.18 NLEBEDEV .................................. 19 4.2.19 MOLDYN ................................... 19 4.2.20 DAMPING ................................... 19 4.2.21 XLBOMD ..................................
    [Show full text]
  • Python Tools in Computational Chemistry (And Biology)
    Python Tools in Computational Chemistry (and Biology) Andrew Dalke Dalke Scientific, AB Göteborg, Sweden EuroSciPy, 26-27 July, 2008 “Why does ‘import numpy’ take 0.4 seconds? Does it need to import 228 libraries?” - My first Numpy-discussion post (paraphrased) Your use case isn't so typical and so suffers on the import time end of the balance. - Response from Robert Kern (Others did complain. Import time down to 0.28s.) 52,000 structures PDB doubles every 2½ years HEADER PHOTORECEPTOR 23-MAY-90 1BRD 1BRD 2 COMPND BACTERIORHODOPSIN 1BRD 3 SOURCE (HALOBACTERIUM $HALOBIUM) 1BRD 4 EXPDTA ELECTRON DIFFRACTION 1BRD 5 AUTHOR R.HENDERSON,J.M.BALDWIN,T.A.CESKA,F.ZEMLIN,E.BECKMANN, 1BRD 6 AUTHOR 2 K.H.DOWNING 1BRD 7 REVDAT 3 15-JAN-93 1BRDB 1 SEQRES 1BRDB 1 REVDAT 2 15-JUL-91 1BRDA 1 REMARK 1BRDA 1 .. ATOM 54 N PRO 8 20.397 -15.569 -13.739 1.00 20.00 1BRD 136 ATOM 55 CA PRO 8 21.592 -15.444 -12.900 1.00 20.00 1BRD 137 ATOM 56 C PRO 8 21.359 -15.206 -11.424 1.00 20.00 1BRD 138 ATOM 57 O PRO 8 21.904 -15.930 -10.563 1.00 20.00 1BRD 139 ATOM 58 CB PRO 8 22.367 -14.319 -13.591 1.00 20.00 1BRD 140 ATOM 59 CG PRO 8 22.089 -14.564 -15.053 1.00 20.00 1BRD 141 ATOM 60 CD PRO 8 20.647 -15.054 -15.103 1.00 20.00 1BRD 142 ATOM 61 N GLU 9 20.562 -14.211 -11.095 1.00 20.00 1BRD 143 ATOM 62 CA GLU 9 20.192 -13.808 -9.737 1.00 20.00 1BRD 144 ATOM 63 C GLU 9 19.567 -14.935 -8.932 1.00 20.00 1BRD 145 ATOM 64 O GLU 9 19.815 -15.104 -7.724 1.00 20.00 1BRD 146 ATOM 65 CB GLU 9 19.248 -12.591 -9.820 1.00 99.00 1 1BRD 147 ATOM 66 CG GLU 9 19.902 -11.351 -10.387 1.00 99.00 1 1BRD 148 ATOM 67 CD GLU 9 19.243 -10.169 -10.980 1.00 99.00 1 1BRD 149 ATOM 68 OE1 GLU 9 18.323 -10.191 -11.782 1.00 99.00 1 1BRD 150 ATOM 69 OE2 GLU 9 19.760 -9.089 -10.597 1.00 99.00 1 1BRD 151 ATOM 70 N TRP 10 18.764 -15.737 -9.597 1.00 20.00 1BRD 152 ATOM 71 CA TRP 10 18.034 -16.884 -9.090 1.00 20.00 1BRD 153 ATOM 72 C TRP 10 18.843 -17.908 -8.318 1.00 20.00 1BRD 154 ATOM 73 O TRP 10 18.376 -18.310 -7.230 1.00 20.00 1BRD 155 .
    [Show full text]
  • Core Software Blocks in Quantum Chemistry: Tensors and Integrals Workshop Program
    Core Software Blocks in Quantum Chemistry: Tensors and Integrals Workshop Program Start: Sunday, May 7, 2017 afternoon. Resort check-in at 4:00 pm. Finish: Wednesday, May 10, noon Lectures are in Scripps, posters are in Heather. Sunday 6:00-7:00: Dinner 7:30-9:00 Opening session 7:30–7:45 Anna Krylov (USC): “MolSSI and some lessons from previous workshop, goals of the workshop” 7:45-8:00 Theresa Windus (Iowa): “Mission of the Molecular Science Consortium” 8:00-9:00 Introduction of participants: 2 min presentation, can have one slide (send in advance) 9:00-10:30 Reception and posters Monday Breakfast: 7:30-9:00 9:00-11:45 Session I: Overview of tensors projects and current developments (moderator Daniel Smith) 9:00-9:10 Daniel Smith (MolSSI): Overview of tensor projects 9:10-9:30 Evgeny Epifanovsky (Q-Chem): Overview of Libtensor 9:30-9:50 Ed Solomonik: "An Overview of Cyclops Tensor Framework" 9:50-10:10 Ed Valeev (VT): "TiledArray: A composable massively parallel block-sparse tensor framework" 10:10-10:30 Coffee break 10:30-10:50 Devin Matthews (UT Austin): "Aquarius and TBLIS: Orthogonal Axes in Multilinear Algebra" 10:50-11:10 Karol Kowalskii (PNNL, NWChem): "NWChem, NWChemEX , and new tensor algebra systems for many-body methods" 11:10-11:30 Peng Chong (VT): "Many-body toolkit in re-designed Massively Parallel Quantum Chemistry package" 11:30-11:55 Moderated discussion: "What problems are we *still* solving?” Lunch: 12:00-1:00 Free time for unstructured discussions 5:00 Posters (coffee/tea) Dinner: 6-7:00 7:15-9:00 Session II: Computer
    [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 ∗[email protected] [email protected] 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]
  • Trends in Atomistic Simulation Software Usage [1.3]
    A LiveCoMS Perpetual Review Trends in atomistic simulation software usage [1.3] Leopold Talirz1,2,3*, Luca M. Ghiringhelli4, Berend Smit1,3 1Laboratory of Molecular Simulation (LSMO), Institut des Sciences et Ingenierie Chimiques, Valais, École Polytechnique Fédérale de Lausanne, CH-1951 Sion, Switzerland; 2Theory and Simulation of Materials (THEOS), Faculté des Sciences et Techniques de l’Ingénieur, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland; 3National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland; 4The NOMAD Laboratory at the Fritz Haber Institute of the Max Planck Society and Humboldt University, Berlin, Germany This LiveCoMS document is Abstract Driven by the unprecedented computational power available to scientific research, the maintained online on GitHub at https: use of computers in solid-state physics, chemistry and materials science has been on a continuous //github.com/ltalirz/ rise. This review focuses on the software used for the simulation of matter at the atomic scale. We livecoms-atomistic-software; provide a comprehensive overview of major codes in the field, and analyze how citations to these to provide feedback, suggestions, or help codes in the academic literature have evolved since 2010. An interactive version of the underlying improve it, please visit the data set is available at https://atomistic.software. GitHub repository and participate via the issue tracker. This version dated August *For correspondence: 30, 2021 [email protected] (LT) 1 Introduction Gaussian [2], were already released in the 1970s, followed Scientists today have unprecedented access to computa- by force-field codes, such as GROMOS [3], and periodic tional power.
    [Show full text]
  • Lawrence Berkeley National Laboratory Recent Work
    Lawrence Berkeley National Laboratory Recent Work Title From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape. Permalink https://escholarship.org/uc/item/4sm897jh Journal Chemical reviews, 121(8) ISSN 0009-2665 Authors Kowalski, Karol Bair, Raymond Bauman, Nicholas P et al. Publication Date 2021-04-01 DOI 10.1021/acs.chemrev.0c00998 Peer reviewed eScholarship.org Powered by the California Digital Library University of California From NWChem to NWChemEx: Evolving with the computational chemistry landscape Karol Kowalski,y Raymond Bair,z Nicholas P. Bauman,y Jeffery S. Boschen,{ Eric J. Bylaska,y Jeff Daily,y Wibe A. de Jong,x Thom Dunning, Jr,y Niranjan Govind,y Robert J. Harrison,k Murat Keçeli,z Kristopher Keipert,? Sriram Krishnamoorthy,y Suraj Kumar,y Erdal Mutlu,y Bruce Palmer,y Ajay Panyala,y Bo Peng,y Ryan M. Richard,{ T. P. Straatsma,# Peter Sushko,y Edward F. Valeev,@ Marat Valiev,y Hubertus J. J. van Dam,4 Jonathan M. Waldrop,{ David B. Williams-Young,x Chao Yang,x Marcin Zalewski,y and Theresa L. Windus*,r yPacific Northwest National Laboratory, Richland, WA 99352 zArgonne National Laboratory, Lemont, IL 60439 {Ames Laboratory, Ames, IA 50011 xLawrence Berkeley National Laboratory, Berkeley, 94720 kInstitute for Advanced Computational Science, Stony Brook University, Stony Brook, NY 11794 ?NVIDIA Inc, previously Argonne National Laboratory, Lemont, IL 60439 #National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37831-6373 @Department of Chemistry, Virginia Tech, Blacksburg, VA 24061 4Brookhaven National Laboratory, Upton, NY 11973 rDepartment of Chemistry, Iowa State University and Ames Laboratory, Ames, IA 50011 E-mail: [email protected] 1 Abstract Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems.
    [Show full text]
  • Openmolcas: from Source Code to Insight Ignacio Fdez
    Article Cite This: J. Chem. Theory Comput. XXXX, XXX, XXX−XXX pubs.acs.org/JCTC OpenMolcas: From Source Code to Insight Ignacio Fdez. Galvan,́ 1,2 Morgane Vacher,1 Ali Alavi,3 Celestino Angeli,4 Francesco Aquilante,5 Jochen Autschbach,6 Jie J. Bao,7 Sergey I. Bokarev,8 Nikolay A. Bogdanov,3 Rebecca K. Carlson,7,9 Liviu F. Chibotaru,10 Joel Creutzberg,11,12 Nike Dattani,13 Mickael̈ G. Delcey,1 Sijia S. Dong,7,14 Andreas Dreuw,15 Leon Freitag,16 Luis Manuel Frutos,17 Laura Gagliardi,7 Fredé rić Gendron,6 Angelo Giussani,18,19 Leticia Gonzalez,́ 20 Gilbert Grell,8 Meiyuan Guo,1,21 Chad E. Hoyer,7,22 Marcus Johansson,12 Sebastian Keller,16 Stefan Knecht,16 Goran Kovacevič ,́ 23 Erik Kallman,̈ 1 Giovanni Li Manni,3 Marcus Lundberg,1 Yingjin Ma,16 Sebastian Mai,20 Joaõ Pedro Malhado,24 Per Åke Malmqvist,12 Philipp Marquetand,20 Stefanie A. Mewes,15,25 Jesper Norell,11 Massimo Olivucci,26,27,28 Markus Oppel,20 Quan Manh Phung,29 Kristine Pierloot,29 Felix Plasser,30 Markus Reiher,16 Andrew M. Sand,7,31 Igor Schapiro,32 Prachi Sharma,7 Christopher J. Stein,16,33 Lasse Kragh Sørensen,1,34 Donald G. Truhlar,7 Mihkel Ugandi,1,35 Liviu Ungur,36 Alessio Valentini,37 Steven Vancoillie,12 Valera Veryazov,12 Oskar Weser,3 Tomasz A. Wesołowski,5 Per-Olof Widmark,12 Sebastian Wouters,38 Alexander Zech,5 J. Patrick Zobel,12 and Roland Lindh*,2,39 1Department of Chemistry − Ångström Laboratory, Uppsala University, P.O.
    [Show full text]
  • Virt&L-Comm.3.2012.1
    Virt&l-Comm.3.2012.1 A MODERN APPROACH TO AB INITIO COMPUTING IN CHEMISTRY, MOLECULAR AND MATERIALS SCIENCE AND TECHNOLOGIES ANTONIO LAGANA’, DEPARTMENT OF CHEMISTRY, UNIVERSITY OF PERUGIA, PERUGIA (IT)* ABSTRACT In this document we examine the present situation of Ab initio computing in Chemistry and Molecular and Materials Science and Technologies applications. To this end we give a short survey of the most popular quantum chemistry and quantum (as well as classical and semiclassical) molecular dynamics programs and packages. We then examine the need to move to higher complexity multiscale computational applications and the related need to adopt for them on the platform side cloud and grid computing. On this ground we examine also the need for reorganizing. The design of a possible roadmap to establishing a Chemistry Virtual Research Community is then sketched and some examples of Chemistry and Molecular and Materials Science and Technologies prototype applications exploiting the synergy between competences and distributed platforms are illustrated for these applications the middleware and work habits into cooperative schemes and virtual research communities (part of the first draft of this paper has been incorporated in the white paper issued by the Computational Chemistry Division of EUCHEMS in August 2012) INTRODUCTION The computational chemistry (CC) community is made of individuals (academics, affiliated to research institutions and operators of chemistry related companies) carrying out computational research in Chemistry, Molecular and Materials Science and Technology (CMMST). It is to a large extent registered into the CC division (DCC) of the European Chemistry and Molecular Science (EUCHEMS) Society and is connected to other chemistry related organizations operating in Chemical Engineering, Biochemistry, Chemometrics, Omics-sciences, Medicinal chemistry, Forensic chemistry, Food chemistry, etc.
    [Show full text]
  • Arxiv:2006.14956V1 [Physics.Comp-Ph]
    NECI: N-Electron Configuration Interaction with an emphasis on state-of-the-art stochastic methods Kai Guther,1, a) Robert J. Anderson,2 Nick S. Blunt,3 Nikolay A. Bogdanov,1 Deidre Cleland,4 Nike Dattani,5 Werner Dobrautz,1 Khaldoon Ghanem,1 Peter Jeszenski,6 Niklas Liebermann,1 Giovanni Li Manni,1 Alexander Y. Lozovoi,1 Hongjun Luo,1 Dongxia Ma,1 Florian Merz,7 Catherine Overy,3 Markus Rampp,8 Pradipta Kumar Samanta,1 Lauretta R. Schwarz,1, 3 James J. Shepherd,9 Simon D. Smart,3 Eugenio Vitale,1 Oskar Weser,1, 2 George H. Booth,2 and Ali Alavi1,3, b) 1)Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany 2)Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom 3)Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom 4)CSIRO Data61, Docklands VIC 3008, Australia 5)Department of Electrical and Computer Engineering, University of Waterloo, 200 University Avenue, Waterloo, Canada 6)Centre for Theoretical Chemistry and Physics, NZ Institute for Advanced Study, Massey University, New Zealand arXiv:2006.14956v1 [physics.comp-ph] 26 Jun 2020 7)Lenovo HPC&AI Innovation Center, Meitnerstr. 9, 70563 Stuttgart 8)Max Planck Computing and Data Facility (MPCDF), Gießenbachstr. 2, 85748 Garching, Germany 9)Department of Chemistry & Informatics Institute, University of Iowa (Dated: 29 June 2020) 1 We present NECI, a state-of-the-art implementation of the Full Configuration Inter- action Quantum Monte Carlo algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function.
    [Show full text]
  • Efficient Open-Source Implementations of Linear-Scaling
    Efficient Open-Source Implementations of Linear-Scaling Polarizable Embedding: Use Octrees to Save the Trees Maximilian Scheurer,∗,y,x Peter Reinholdt,∗,z,x Jógvan Magnus Haugaard Olsen,{ Andreas Dreuw,y and Jacob Kongstedz yInterdisciplinary Center for Scientific Computing, Heidelberg University, D-69120 Heidelberg, Germany zDepartment of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark {DTU Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark xContributed equally to this work E-mail: [email protected]; [email protected] Abstract We present open-source implementations of the linear-scaling Fast Multipole Method (FMM) within the Polarizable Embedding (PE) model for efficient treatment of large polarizable environments in computational spectroscopy simulations. The implementa- tions are tested for accuracy, efficiency, and usability on model systems as well as more realistic biomolecular systems. We explain how FMM parameters affect the calculation of molecular properties and show that PE calculations employing FMM can be carried out in a black-box manner. The efficiency of the linear-scaling approach is demonstrated by simulating the UV/Vis spectrum of a chromophore in an environment of more than 1 one million polarizable sites. Our implementations are interfaced to several open-source quantum chemistry programs, making computational spectroscopy simulations within the PE model and FMM available to a large variety of methods and a broad user base. 1 Introduction
    [Show full text]
  • Open Source Molecular Modeling
    Accepted Manuscript Title: Open Source Molecular Modeling Author: Somayeh Pirhadi Jocelyn Sunseri David Ryan Koes PII: S1093-3263(16)30118-8 DOI: http://dx.doi.org/doi:10.1016/j.jmgm.2016.07.008 Reference: JMG 6730 To appear in: Journal of Molecular Graphics and Modelling Received date: 4-5-2016 Accepted date: 25-7-2016 Please cite this article as: Somayeh Pirhadi, Jocelyn Sunseri, David Ryan Koes, Open Source Molecular Modeling, <![CDATA[Journal of Molecular Graphics and Modelling]]> (2016), http://dx.doi.org/10.1016/j.jmgm.2016.07.008 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. Open Source Molecular Modeling Somayeh Pirhadia, Jocelyn Sunseria, David Ryan Koesa,∗ aDepartment of Computational and Systems Biology, University of Pittsburgh Abstract The success of molecular modeling and computational chemistry efforts are, by definition, de- pendent on quality software applications. Open source software development provides many advantages to users of modeling applications, not the least of which is that the software is free and completely extendable. In this review we categorize, enumerate, and describe available open source software packages for molecular modeling and computational chemistry. 1. Introduction What is Open Source? Free and open source software (FOSS) is software that is both considered \free software," as defined by the Free Software Foundation (http://fsf.org) and \open source," as defined by the Open Source Initiative (http://opensource.org).
    [Show full text]