DOCSLIB.ORG

Analytic MR-CISD and MR-AQCC Gradients and MR-AQCC-LRT for Excited States, GUGA Spin–Orbit CI and Parallel CI Density¤

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Analytic MR-CISD and MR-AQCC Gradients and MR-AQCC-LRT for Excited States, GUGA Spin–Orbit CI and Parallel CI Density¤ High-level multireference methods in the quantum-chemistry program system COLUMBUS: Analytic MR-CISD and MR-AQCC gradients and MR-AQCC-LRT for excited states, GUGA spin–orbit CI and parallel CI density¤ Hans Lischka,*a Ron Shepard,*b Russell M. Pitzer,*c Isaiah Shavitt,*cd Michal Dallos,a ThomasMu ller,a Pe ter G. Szalay,*e Michael Seth,f Gary S. Kedziora,g Satoshi Yabushitah and Zhiyong Zhangi a Institute for T heoretical Chemistry and Structural Biology, University of V ienna, Wa hringerstrasse17, A-1090 V ienna, Austria. E-mail: Hans.L ischka=univie.ac.at b T heoretical Chemistry Group, Chemistry Division, Argonne National L aboratory, Argonne, IL 60439, USA. E-mail: shepard.=tcg.anl.gov c Department of Chemistry, T he Ohio State University, 100 W est 18th Avenue, Columbus, OH 43210, USA. E-mail: pitzer=chemistry.ohio-state.edu, shavitt=chemistry.ohio-state.edu d Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA e Department of T heoretical Chemistry, Eo tvo s L ora nd University, P.O. Box 32, H-1518 Budapest, Hungary. E-mail: szalay=para.chem.elte.hu f Department of Chemistry, University of Calgary, University Drive 2500, Calgary, Alberta, Canada T 2N 1N4 g Department of Chemistry, Northwestern University, Evanston, IL 60208, USA h Department of Chemistry, Faculty of Science and T echnology, Keio University, 3-14-1 Hyoshi, Kohoku-ku, Y okohama 223-8522, Japan i PaciÐc Northwest National L aboratory, P.O. Box 999, Richland, W A 99352, USA Received 5th October 2000, Accepted 31st October 2000 First published as an Advance Article on the web 18th December 2000 Development of several new computational approaches within the framework of multi-reference ab initio molecular electronic structure methodology and their implementation in the COLUMBUS program system are reported. These new features are: calculation of the analytical MR-CI gradient for excited states based on state-averaged MCSCF orbitals, the extension of the MR-ACPF/AQCC methods to excited states in the framework of linear-response theory, spinÈorbit CI for molecules containing heavy atoms and the development of a massively-parallel code for the computation of the one- and two-particle density matrix elements. Illustrative examples are given for each of these cases. 1. Introduction of a calculation as compared to the SR case, especially in the selection of appropriate reference spaces. The popularity of quantum chemistry ab initio methods in An MR treatment usually starts with a multiconÐguration computational chemistry has increased rapidly in recent years SCF (MCSCF) calculation, which describes the nondynamical due to the development of new methods and of corresponding electron correlation by including the nearly degenerate elec- efficient computer programs that take advantage of modern tron conÐgurations (and often many others) in the MCSCF computer technology. The most popular methods by far have wavefunction. Dynamical electron correlation can then be been of the single-reference (SR) type, including approaches included in several ways. The classical method for that such as the HartreeÈFock self-consistent Ðeld (SCF) method, purpose is conÐguration interaction1,2 (CI), using the Ritz MÔllerÈPlesset perturbation theory, coupled-cluster theory variation principle. The most common form of this approach and density-functional theory (DFT). However, the SR is MR-CI singles and doubles (MR-CISD), in which the approaches often fail when nondynamical electron correlation expansion space for the wavefunction is constructed from (near-degeneracy e†ects) is important, as usually happens in single and double substitutions of occupied orbitals by virtual the description of bond stretching and dissociation, elec- orbitals in the individual conÐguration state functions (CSFs) tronically excited states and many open-shell cases. Multi- of the reference wavefunction. The truncation of the expansion reference (MR) methods can deal e†ectively with such space to single and double substitutions, while physically problems, but their successful application requires signiÐ- motivated, is usually mandated by the very steep increase in cantly more careful planning and understanding of the details the number of CSFs, and the consequent computational e†ort, with increases in the substitution level. ¤ Presented at the Third European Conference on Computational The CI method is very Ñexible and robust, but the greatest Chemistry, Budapest, Hungary, September 4È8, 2000. disadvantage of truncated CI is its lack of extensivity, i.e., its 664 Phys. Chem. Chem. Phys., 2001, 3, 664È673 DOI: 10.1039/b008063m This journal is( The Owner Societies 2001 incorrect scaling with the size of the system.3,4 A variety of this implementation has been included in the current version correction methods have been developed for the approximate of COLUMBUS,35 and has been applied to a number of treatment of this problem, but the use of extensivity correc- systems containing very heavy atoms. tions is by no means as straightforward in MR-CI as in the The COLUMBUS program project28,29,35 was started in SR case. We shall mention just a few of the many extensivity 1980 with the aim of developing an efficient and Ñexible tool correction approaches. The simplest and most commonly used for MR-CISD calculations, with emphasis on the multi- are the Davidson correction5h7 and its MR extension,8,9 reference character of the computational approaches. An early which provide a correction to the computed energy, but not to version included MR-LCCM, and the MR-ACPF and the wavefunction or to computed properties. These correc- MR-AQCC methods were added as soon as they became tions are obtained as relatively cheap byproducts of SR-CISD available. The implementation was constructed as a set of and MR-CISD calculations, respectively, but are of limited individual programs for the computation of AO integrals and reliability.10 In order to address the previously mentioned for the various stages of the SCF, MCSCF and CI procedures. single reference limitations, multireference extensions of These programs communicated with each other through Ðles. coupled cluster theory (see, e.g., ref. 11È15) and perturbation The modular and open structure facilitated the implementa- methods (see, e.g., ref. 16 and 17) have also been developed. tion of new methods and features and interaction with other However, there is still no widely applicable and generally program systems. It is the purpose of this work to present new available MR-coupled-cluster procedure. Various approx- methodological developments, their implementation in imations have been introduced, mostly as modiÐcations to the COLUMBUS, and illustrative applications. An important MR-CISD procedure, to deal with this problem. Of these aspect of the current version concerns excited state calcu- methods we shall mention MR-LCCM12 (linearized coupled lations, including implementation of analytic gradients for cluster), MR-ACPF18 (averaged coupled-pair functional) and MR-CISD and MR-AQCC calculations based on state- MC-CEPA19 (coupled electron-pair approximation). In many averaged MCSCF orbitals and the treatment of extensivity applications the MR-ACPF approach has been quite suc- corrections by means of MR-AQCC and MR-AQCC-LRT cessful, but it tends to overestimate the e†ect of the omitted (linear-response theory).36 These features signiÐcantly increase higher excitations, sometimes producing anomalous results.20 the capabilities for accurate geometry optimizations of excited The MR-AQCC (averaged quadratic coupled-cluster) states. The other major new development is the spinÈorbit CI method,20,21 which is closely related to MR-ACPF, was for- implementation, which opens the way for completely new mulated in an attempt to overcome these problems, and applications. In addition, the e†ort to develop a version of usually produces more reliable results. The usefulness and COLUMBUS for massively parallel computers has continued, e†ectiveness of MR-ACPF and MR-AQCC are due to a good adding an efficient parallel version of the one- and two- balance of simplicity and satisfactory treatment of the basic particle density-matrix program to the previously available size-extensivity e†ect. The availability of analytic gradients for parallel modules. This feature is very important for efficient these models (see below) is another important point in their analytic gradient calculations. favor. For a review of these and related methods see, e.g., ref. 22. The availability of analytic gradient methods is crucial for 2. Synopsis of COLUMBUS the geometry optimization of molecules. Analytic gradient methods are often unavailable for MR calculations because of The MCSCF and CI wavefunctions are speciÐed using the the much greater complexity of the formalisms23h27 compared graphical unitary group approach (GUGA).37h40 This allows to the SR case. The situation becomes particularly compli- the treatment of any pure spin state with no practical cated if the invariance properties of the wavefunction with restriction with respect to the number of open shells in the respect to orbital transformations di†er in the MCSCF and individual CSF expansion terms. Mixing of di†erent spin MR-CI wavefunctions. Such situations will arise, for example, states through the spinÈorbit interaction is included in the if orbitals are frozen at the CI stage or if the MR-CI reference spinÈorbit capability. The reference space and CI expansion wavefunction includes only a selected subset of the CSFs in space are deÐned by means of a distinct row table (DRT), the original MCSCF wavefunction. Additional requirements which is a tabular representation of a distinct row graph.39,40 arise in calculations of electronically excited states. For a bal- The graph comprises two parts: the internal part describes the anced description of a series of states, state-averaged MCSCF occupations and spin couplings of the internal orbitals (the calculations are very useful in order to provide a common set orbitals occupied in at least one of the reference of orbitals for all states of interest.
Load more

Recommended publications
  • The Molecular Sciences Software Institute
    The Molecular Sciences Software Institute T. Daniel Crawford, Cecilia Clementi, Robert Harrison, Teresa Head-Gordon, Shantenu Jha*, Anna Krylov, Vijay Pande, and Theresa Windus http://molssi.org S2I2 HEP/CS Workshop at NCSA/UIUC 07 Dec, 2016 1 Outline • Space and Scope of Computational Molecular Sciences. • “State of the art and practice” • Intellectual drivers • Conceptualization Phase: Identifying the community and needs • Bio-molecular Simulations (BMS) Conceptualization • Quantum Mechanics/Chemistry (QM) Conceptualization • Execution Phase. • Structure and Governance Model • Resource Distribution • Work Plan 2 The Molecular Sciences Software Institute (MolSSI) • New project (as of August 1st, 2016) funded by the National Science Foundation. • Collaborative effort by Virginia Tech, Rice U., Stony Brook U., U.C. Berkeley, Stanford U., Rutgers U., U. Southern California, and Iowa State U. • Total budget of $19.42M for five years, potentially renewable to ten years. • Joint support from numerous NSF divisions: Advanced Cyberinfrastructure (ACI), Chemistry (CHE), Division of Materials Research (DMR), Office of Multidisciplinary Activities (OMA) • Designed to serve and enhance the software development efforts of the broad field of computational molecular science. 3 Computational Molecular Sciences (CMS) • The history of CMS – the sub-fields of quantum chemistry, computational materials science, and biomolecular simulation – reaches back decades to the genesis of computational science. • CMS is now a “full partner with experiment”. • For an impressive array of chemical, biochemical, and materials challenges, our community has developed simulations and models that directly impact: • Development of new chiral drugs; • Elucidation of the functionalities of biological macromolecules; • Development of more advanced materials for solar-energy storage, technology for CO2 sequestration, etc.
    [Show full text]
  • SCF and CI Studies on Ground State of Water Molecule
    Ab initio SCF and CI studies on the ground state of the water molecule. II. Potential energy and property surfaces Bruce J. Rosenberg, Walter C. Ermler, and Isaiah Shavitt Citation: J. Chem. Phys. 65, 4072 (1976); doi: 10.1063/1.432861 View online: http://dx.doi.org/10.1063/1.432861 View Table of Contents: http://jcp.aip.org/resource/1/JCPSA6/v65/i10 Published by the American Institute of Physics. Additional information on J. Chem. Phys. Journal Homepage: http://jcp.aip.org/ Journal Information: http://jcp.aip.org/about/about_the_journal Top downloads: http://jcp.aip.org/features/most_downloaded Information for Authors: http://jcp.aip.org/authors Downloaded 29 Apr 2013 to 134.76.213.205. This article is copyrighted as indicated in the abstract. Reuse of AIP content is subject to the terms at: http://jcp.aip.org/about/rights_and_permissions Ab initio SCF and CI studies on the ground state of the water molecule. II. Potential energy and property surfaces Bruce J. Rosenberg* and Walter C. Ermlert Department of Chemistry. The Ohio State University. Columbus. Ohio 43210 Isaiah Shavitt Department of Chemistry. The Ohio State University. Columbus. Ohio 43210 and Battelle Memorial Institute. Columbus. Ohio 43201 (Received 29 June 1976) Self-consistent field and configuration interaction calculations for the energy and one-electron properties of the ground state of the water molecule were carried out with a (5s4p2d/3s1p) 39-function STO basis set. The CI treatment included all single and double excitation configurations (SD) relative to the SCF configuration. and a simple formula due to Davidson was used to estimate the energy contribution of quadruple excitations and thus produce a set of corrected (SDQ) energies.
    [Show full text]
  • Arxiv:1911.06836V1 [Physics.Chem-Ph] 8 Aug 2019 A
    Multireference electron correlation methods: Journeys along potential energy surfaces Jae Woo Park,1, ∗ Rachael Al-Saadon,2 Matthew K. MacLeod,3 Toru Shiozaki,2, 4 and Bess Vlaisavljevich5, y 1Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea. 2Department of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, IL 60208, USA. 3Workday, 4900 Pearl Circle East, Suite 100, Boulder, CO 80301, USA. 4Quantum Simulation Technologies, Inc., 625 Massachusetts Ave., Cambridge, MA 02139, USA. 5Department of Chemistry, University of South Dakota, 414 E. Clark Street, Vermillion, SD 57069, USA. (Dated: November 19, 2019) Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way, and therefore, can yield accurate and predictive results even when single-reference methods or multicon- figurational self-consistent field (MCSCF) theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces (PES). This includes the optimization of molecular ge- ometries, such as equilibrium geometries and conical intersections, and on-the-fly photodynamics simulations; both depend heavily on the ability of the method to properly explore the PES. Since such applications require the nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly improves the utility of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. To motivate the readers, we first summarize the notable applications of multireference elec- tron correlation methods to mainstream chemistry, including geometry optimizations and on-the-fly dynamics. Subsequently, we review the analytical nuclear gradient and derivative coupling theories for these methods, and the software infrastructure that allows one to make use of these quantities in applications.
    [Show full text]
  • Perturbational Treatment of Spin-Orbit Coupling for Generally Applicable High-Level Multi-Reference Methods
    Perturbational treatment of spin-orbit coupling for generally applicable high-level multi-reference methods Sebastian Mai1, Thomas Müller2,a), Felix Plasser3, Philipp Marquetand1, Hans Lischka and Leticia González1 1Institute of Theoretical Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria 2Institute for Advanced Simulation, Jülich Supercomputing Centre, Forschungszentrum Jülich, 52425 Jülich, Germany 3Interdisciplinary Center for Scientific Computing, University of Heidelberg, Im Neuenheimer Feld 368, 69120 Heidelberg, Germany 4Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA Abstract An efficient perturbational treatment of spin-orbit coupling within the framework of high-level multi-reference techniques has been implemented in the most recent version of the COLUMBUS quantum chemistry package, extending the existing fully variational two-component (2c) multi-reference configuration interaction singles and doubles (MRCISD) method. The proposed scheme follows related implementations of quasi-degenerate perturbation theory (QDPT) model space techniques. Our model space is built either from uncontracted, large-scale scalar relativistic MRCISD wavefunctions or based on the scalar-relativistic solutions of the linear-response-theory-based multi-configurational averaged quadratic coupled cluster method (LRT-MRAQCC). The latter approach allows for a consistent, approximatively size-consistent and size-extensive treatment of spin-orbit coupling. The approach is described in
    [Show full text]
  • Curriculum Vitae Rodney J
    Revised 4/29/2017 CURRICULUM VITAE RODNEY J. BARTLETT GRADUATE RESEARCH PROFESSOR OF CHEMISTRY AND PHYSICS QUANTUM THEORY PROJECT DEPARTMENTS OF CHEMISTRY AND PHYSICS UNIVERSITY OF FLORIDA PO BOX 118435 GAINESVILLE, FLORIDA USA 32611-8435 Phone: 352-392-6974 Email: [email protected] Web Page: http://www.clas.ufl.edu/users/rodbartl BIRTH DATE: March 31, 1944 CV INDEX AREA OF SPECIALIZATION Quantum chemistry, molecular electronic structure and spectra, ab initio many-electron methods Following the initial formulation of Jiri Cizek of coupled-cluster theory with double excitations, CCD, Rod Bartlett pioneered the development of coupled-cluster (CC) theory in quantum chemistry to offer highly accurate solutions of the Schroedinger equation for molecular structure and spectra. He introduced the term ‘size-extensivity’ for many-body methods like CC that scale properly with the number of electrons; now viewed as an essential element in quantum chemistry approximations. He and his co-workers were the first to formulate and implement CC theory with all single and double excitations (CCSD), to add triples non- iteratively (CCSD[T]), iteratively, CCSDT-1, and fully, CCSDT; followed by quadruple (CCSDTQ) and pentuple excitations (CCSDTQP). He also presented a density matrix formulation for the analytical gradients (forces) for the non- variational CC method, a necessity for any widely used method in quantum chemistry. This led to the CC functional, E=0|(1+)exp(-T)Hexp(T)|0 where the forces on an atom, , become E/X=0|(1+)exp(-T)(H/X)exp(T)|0. This functional also defines the CC response and relaxed density matrices for all properties, a non-Hermitian generalization of the conventional density matrices.
    [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]
  • Nwchem User Documentation Release 4.0.1
    NWChem User Documentation Release 4.0.1 High Performance Computational Chemistry Group W.R. Wiley Environmental Molecular Sciences Laboratory Pacific Northwest National Laboratory P.O. Box 999, Richland, WA 99352 January 2001 2 DISCLAIMER This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor Battelle, nor any of their employees, MAKES ANY WARRANTY, EXPRESS OR IMPLIED, OR ASSUMES ANY LEGAL LIABILITY OR RESPON- SIBILITY FOR THE ACCURACY, COMPLETENESS, OR USEFULNESS OF ANY INFORMATION, APPARA- TUS, PRODUCT, SOFTWARE, OR PROCESS DISCLOSED, OR REPRESENTS THAT ITS USE WOULD NOT INFRINGE PRIVATELY OWNED RIGHTS. LIMITED USE This software (including any documentation) is being made available to you for your internal use only, solely for use in performance of work directly for the U.S. Federal Government or work under contracts with the U.S. Department of Energy or other U.S. Federal Government agencies. This software is a version which has not yet been evaluated and cleared for commercialization. Adherence to this notice may be necessary for the author, Battelle Memorial Institute, to successfully assert copyright in and commercialize this software. This software is not intended for duplication or distribution to third parties without the permission of the Manager of Software Products at Pacific Northwest National Laboratory, Richland, Washington, 99352. ACKNOWLEDGMENT This software and its documentation were produced with Government support under Contract Number DE-AC06- 76RLO-1830 awarded by the United States Department of Energy. The Government retains a paid-up non-exclusive, irrevocable worldwide license to reproduce, prepare derivative works, perform publicly and display publicly by or for the Government, including the right to distribute to other Government contractors.
    [Show full text]
  • Mathematics and Computer Science
    ANL/MC S-TM-219 Workshop Report on Large-Scale Matrix Diagonalization Methods in Chemistry Theory Institute held at Argonne National Laboratory May 20-22, 1996 edited by Cliristian H. Bischof, Ron L. Shepard, and Steven Huss-Lederman MATHEMATICS AND COMPUTER SCIENCE DIVISION Argonne Xational Laboratory, with facilities in the states of Illinois and Idaho, is owned by the United States government, and operated by The University of Chicago under the provisions of a contract with the Department of Energy. DISCLAIMER This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of theiiemployees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or pro- cess disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. .. I. Reproduced from the best available copy. Available to and DOE contractors from the Office of ScientificDOE and Technical Information P.O.Box 62 Oak Ridge, TN 3783 1 Prices available from (423) 576-8401 Available to the public from the National Technical information Service U.S. Department of Commerce 5285 Port Royal Road Springfield.
    [Show full text]
  • Information to Users
    INFORMATION TO USERS This manuscript has been reproduced from the microfilm master. UMI films the text directly from the original or copy submitted. Thus, some thesis and dissertation copies are in typewriter face, while others may be from any type of computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion. Oversize materials (e.g., maps, drawings, charts) are reproduced by sectioning the original, beginning at the upper left-hand corner and continuing from left to right in equal sections with small overlaps. Each original is also photographed in one exposure and is included in reduced form at the back of the book. Photographs included in the original manuscript have been reproduced xerographically in this copy. Higher quality 6" x 9" black and white photographic prints are available for any photographs or illustrations appearing in this copy for an additional charge. Contact UMI directly to order. University Microfilms International A Bell & Howell Information Company 300 North Zeeb Road, Ann Arbor, Ml 48106-1346 USA 313/761-4700 800/521-0600 Order Number 9211227 Application of multireference-based correlation methods to the study of weak bonding interactions Stahlberg, Eric Alan, Ph.D. The Ohio State University, 1991 UMI 300 N.
    [Show full text]
  • Computational Chemistry: a Practical Guide for Applying Techniques to Real-World Problems
    Computational Chemistry: A Practical Guide for Applying Techniques to Real-World Problems. David C. Young Copyright ( 2001 John Wiley & Sons, Inc. ISBNs: 0-471-33368-9 (Hardback); 0-471-22065-5 (Electronic) COMPUTATIONAL CHEMISTRY COMPUTATIONAL CHEMISTRY A Practical Guide for Applying Techniques to Real-World Problems David C. Young Cytoclonal Pharmaceutics Inc. A JOHN WILEY & SONS, INC., PUBLICATION New York . Chichester . Weinheim . Brisbane . Singapore . Toronto Designations used by companies to distinguish their products are often claimed as trademarks. In all instances where John Wiley & Sons, Inc., is aware of a claim, the product names appear in initial capital or all capital letters. Readers, however, should contact the appropriate companies for more complete information regarding trademarks and registration. Copyright ( 2001 by John Wiley & Sons, Inc. All rights reserved. No part of this publication may be reproduced, stored in a retrieval system or transmitted in any form or by any means, electronic or mechanical, including uploading, downloading, printing, decompiling, recording or otherwise, except as permitted under Sections 107 or 108 of the 1976 United States Copyright Act, without the prior written permission of the Publisher. Requests to the Publisher for permission should be addressed to the Permissions Department, John Wiley & Sons, Inc., 605 Third Avenue, New York, NY 10158-0012, (212) 850-6011, fax (212) 850-6008, E-Mail: PERMREQ @ WILEY.COM. This publication is designed to provide accurate and authoritative information in regard to the subject matter covered. It is sold with the understanding that the publisher is not engaged in rendering professional services. If professional advice or other expert assistance is required, the services of a competent professional person should be sought.
    [Show full text]
  • A Case Study with Nwchem
    CONCURRENCY AND COMPUTATION: PRACTICE AND EXPERIENCE Concurrency Computat.: Pract. Exper. 2010; 00:1–17 Published online in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cpe Performance Characterization of Global Address Space Applications: A Case Study with NWChem Jeff R. Hammond∗, Sriram Krishnamoorthy+, Sameer Shende‡ Nichols A. Romero∗, Allen D. Malony‡ ∗ Argonne National Laboratory + Pacific Northwest National Laboratory ‡ University of Oregon SUMMARY The use of global address space languages and one-sided communication for complex applications is gaining attention in the parallel computing community. However, lack of good evaluative methods to observe multiple levels of performance makes it difficult to isolate the cause of performance deficiencies and to understand the fundamental limitations of system and application design for future improvement. NWChem is a popular computational chemistry package which depends on the Global Arrays / ARMCI suite for partitioned global address space functionality to deliver high-end molecular modeling capabilities. A workload characterization methodology was developed to support NWChem performance engineering on large-scale parallel platforms. The research involved both the integration of performance instrumentation and measurement in the NWChem software, as well as the analysis of one-sided communication performance in the context of NWChem workloads. Scaling studies were conducted for NWChem on Blue Gene/P and on two large-scale clusters using different generation Infiniband interconnects and x86 processors. The performance analysis and results show how subtle changes in the runtime parameters related to the communication subsystem could have significant impact on performance behavior. The tool has successfully identified several algorithmic bottlenecks which are already being tackled by computational chemists to improve NWChem performance.
    [Show full text]
  • Adiabatic and Nonadiabatic Molecular Dynamics with Multireference Ab
    Columbus in Rio: univie.ac.at/columbus/rio Rio de Janeiro, Nov. 27- Dec. 02, 2005 AdiabaticAdiabatic andand nonadiabaticnonadiabatic molecularmolecular dynamicsdynamics withwith multireferencemultireference abab initioinitio methodsmethods Mario Barbatti Institute for Theoretical Chemistry University of Vienna Outline First Lecture: An introduction to molecular dynamics 1. Dynamics, why? 2. Overview of the available approaches Second Lecture: Towards an implementation of surface hopping dynamics 1. Practical aspects to be addressed 2. The NEWTON-X program 3. Some applications: theory and experiment PartPart IIII TowardsTowards anan implementationimplementation ofof surfacesurface hoppinghopping dynamicsdynamics PracticalPractical aspectsaspects toto bebe addressedaddressed InitialInitial conditionsconditions 2 2 (0) 1 ⎡− 2α(R − Re ) ⎤ ⎡− P ⎤ PW (R, P) = exp⎢ ⎥exp⎢ ⎥ E πh ⎣ h ⎦ ⎣ 2αh ⎦ mω α = HO 2 E0 ± ∆E Sn Wigner distribution S0 accept R,P don’t accept ClassicalClassical dynamics:dynamics: integratorintegrator Any standard method can be used in the integration of the Newton equations. A good one is the Velocity Verlet (Swope et al. JCP 76, 637 (1982)): For each nucleus I 1 R (t + ∆t) = R (t) + v (t)∆t + a (t)∆t 2 I I I 2 I ⎛ ∆t ⎞ 1 vI ⎜t + ⎟ = vI (t) + aI (t)∆t ⎝ 2 ⎠ 2 1 aI (t) = − ∇R E[]R I (t + ∆t) Quantum chemistry calculation M I ⎛ ∆t ⎞ 1 vI (t + ∆t) = vI ⎜t + ⎟ + aI (t + ∆t)∆t ⎝ 2 ⎠ 2 TimeTime-step-step forfor thethe classicalclassicalequationsequations Schlick, Barth and Mandziuk, Annu. Rev. Biophys. Struct. 26, 181 (1997). TimeTime-step-step forfor thethe classicalclassicalequationsequations Time step should not be larger than 1 fs (1/10v). ∆t = 0.5 fs assures a good level of conservation of energy most of the time.
    [Show full text]