DOCSLIB.ORG

Intro to Real-Space, Linear-Response, and TD Methods July 19, 2011 1 / 62 Outline

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Intro to Real-Space, Linear-Response, and TD Methods July 19, 2011 1 / 62 Outline Introduction to real-space, linear-response, and time-dependent methods Heiko Appel Fritz-Haber-Institut der Max-Planck-Gesellschaft, Berlin Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 1 / 62 Outline Linear Response in DFT I Response functions I Casida equation I Sternheimer equation Real-space representation and real-time propagation I Real-space representation for wavefunctions and Hamiltonians I Time-propagation schemes I Optimal control of electronic motion Time-evolution of open quantum systems I Stochastic Schr¨odingerequations, master equations I Stochastic current DFT I Stochastic quantum molecular dynamics Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 2 / 62 Motivation Where is electron dynamics important? I Electron-hole pair creation in solar cells I Photosynthesis and energy transfer in light-harvesting antenna complexes I Quantum computing (e.g. electronic transitions in ultracold atoms) I Molecular electronics, quantum transport Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 3 / 62 Time-dependent density-functional theory I One-to-one correspondence of time-dependent densities and potentials v(r; t) 1−!1 ρ(r; t) For fixed inital states, the time-dependent density determines uniquely the time-dependent external potential and hence all physical observables. I Time-dependent Kohn-Sham system The time-dependent density of an interacting many-electron system can be calculated as density N X 2 ρ(r; t) = j'j (r; r)j j=1 of an auxiliary non-interacting Kohn-Sham system 2r2 i @ ' (r; t) = − ~ + v [ρ](r; t) ' (r; t) ~ t j 2m S j with a local multiplicative potential Z ρ(r0; t) v [ρ(r0; t0)](r; t)= v(r; t) + d3r0 + v [ρ(r0; t0)](r; t) S jr − r0j xc E. Runge, and E.K.U. Gross, Phys. Rev. Lett. 52, 997 (1984). Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 4 / 62 Linear Response Theory I Hamiltonian H^ (t) = H^0 + Θ(t − t0)v1(r; t) I Initial condition: for times t < t0 the system is in the ground-state of the unperturbed Hamiltonian H^0 with potential v0 and density ρ0(r) I For times t > t0, switch on perturbation v1(r; t). Leads to time-dependent density ρ(r; t) = ρ0(r) + δρ(rt) I Functional Taylor expansion of ρ[v](r; t) around v0: ρ[v](r; t) = ρ[v0 + v1](r; t) = ρ[v0](r; t) Z δρ[v](rt) + v (r0t0)d3r0dt0 0 0 1 δv(r t ) v0 ZZ δ2ρ[v](rt) + v (r0t0)v (r00t00)d3r0dt0d3r00dt00 0 0 00 00 1 1 δv(r t )δv(r t ) v0 + ::: Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 5 / 62 Computing Linear Response Different ways to compute first order response in DFT I Response functions, Casida equation I (frequency-dependent) perturbation theory, Sternheimer equation I real-time propagation with weak external perturbation Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 6 / 62 Response functions I Functional Taylor expansion of ρ[v](r; t) around external potential v0: Z δρ[v](rt) ρ[v + v ](r; t) = ρ[v ](r) + v (r0t0)d3r0dt0 + ::: 0 1 0 0 0 1 δv(r t ) v0 I Density-density response function of interacting system δρ[v](rt) χ(rt; r0t0) := 0 0 δv(r t ) v0 I Response of non-interacting Kohn-Sham system: Z δρ[v ](rt) ρ[v + v ](r; t) = ρ[v ](r) + S v (r0t0)d3r0dt0 + ::: S;0 S;1 S;0 0 0 S δvS (r t ) v0 I Density-density response function of time-dependent Kohn-Sham system δρ [v ](rt) χ (rt; r0t0) := S S S 0 0 δvS (r t ) vS;0 Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 7 / 62 Derivation of response equation I Definition of time-dependent xc potential vxc(rt) = vKS (rt) − vext(rt) − vH (rt) I Take functional derivative δv (rt) δv (rt) δv (rt) δ(t − t0) xc = KS − ext − δρ(r0t0) δρ(r0t0) δρ(r0t0) jr − r0j 0 0 −1 0 0 −1 0 0 0 0 fxc(rt; r t )= χS (rt; r t ) − χ (rt; r t ) − Wc(rt; r t ) I Act with reponse functions from left and right −1 −1 χS ·j Wc + fxc = χS − χ j· χ χS (Wc + fxc)χ = χ − χS I Dyson-type equation for response functions χ = χS + χS (Wc + fxc)χ Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 8 / 62 First order density response I Exact density response to first order ρ1 = χv1 = χS v1 + χS (Wc + fxc)ρ1 I In integral notation Z 3 0 0 0 0 h 0 0 ρ1(rt) = d r dt χS (rt; r t ) v1(r t ) Z 3 00 00 0 0 00 00 0 0 00 00 00 00 i + d r dt (Wc(r t ; r t ) + fxc(r t ; r t ))ρ1(r t ) I For practical application: iterative solution with approximate kernel fxc δv [ρ](r0t0) f (r0t0; r00t00) = xc xc 00 00 δρ(r t ) ρ0 Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 9 / 62 Lehmann representation of linear response function I Exact many-body eigenstates H^ (t = t0)jmi = Emjmi I Lehmann representation of linear response function: 0 0 0 X h0jρ^(r)jmihmjρ^(r)j0i h0jρ^(r )jmihmjρ^(r )j0i χ(r; r ; !) = lim − η!0+ ! − (Em − E0) + iη ! + (Em − E0) + iη m Neutral excitation energies are poles of the linear response function! I Exact linear density response to perturbation v1(!) ρ1(!) =χ ^(!)v1(!) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 10 / 62 Excitation energies I Dyson-type equation for response functions in frequency space [1^ − χ^S (!)(W^ c + f^xc(!))]ρ1(!) = χS v1(!) I ρ1(!) has poles for exact excitation energies Ωj ρ1(!) ! 1 for ! ! Ωj I On the other hand, rhs χS v1(!) stays finite for ! ! Ωj hence the eigenvalues of the integral operator [1^ − χ^S (!)(W^ c + f^xc(!))]ξ(!) = λ(!)ξ(!) vanish, λ(!) ! 0 for ! ! Ωj . I Determines rigorously the exact excitation energies [1^ − χ^S (Ωj )(W^ c + f^xc(Ωj ))]ξ(Ωj ) = 0 Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 11 / 62 Casida equation I (Non-linear) eigenvalue equation for excitation energies 2 ΩFj = !j Fj with 2 p p Ωiaσ;jbτ = δσ,τ δi;j δa;b(a − i) + 2 (a − i)Kiaσ;jbτ (b − j ) and Z Z h 1 i K (!)= d3r d3r0φ (r)φ (r) + f (r; r0;!) φ (r)φ (r) iaσ;jbτ iσ jσ jr − r0j xc kτ lτ I Eigenvalues !j are exact vertical excitation energies I Eigenvectors can be used to compute oscillator strength I Drawback: need occupied and unoccupied orbitals Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 12 / 62 Adiabatic approximation I Adiabatic approximation: evaluate static Kohn-Sham potential at time-dependent density adiab static DFT vxc [ρ](rt) := vxc [ρ(t)](rt) I Example: adiabatic LDA ALDA LDA 1=3 vxc [ρ](rt) := vxc (ρ(t)) = −αρ(r; t) + ::: I Exchange-correlation kernel δvALDA[ρ](rt) @vALDA f ALDA(rt; r0t0) = xc = δ(t − t0)δ(r − r0) xc xc 0 0 δρ(r t ) @ρ(r) ρ0(r) @2ehom = δ(t − t0)δ(r − r0) xc 2 @n ρ0(r) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 13 / 62 Failures of the adiabatic approximation in linear response I H2 dissociation is incorrect 1 + 1 + R!1 E( Σu ) − E( Σg ) −! 0 (in ALDA) Gritsenko, van Gisbergen, Grling, Baerends, JCP 113, 8478 (2000). I sometimes problematic close to conical intersections I response of long chains strongly overestimated Champagne et al., JCP 109, 10489 (1998) and 110, 11664 (1999). I in periodic solids fxc(q; !; ρ) = c(ρ), whereas for insulators, exact q!0 2 fxc −! 1=q divergent I charge transfer excitations not properly described Dreuw et al., JCP 119, 2943 (2003). Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 14 / 62 Sternheimer equation Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 15 / 62 Sternheimer equation I Perturbed Hamiltonian and states (zero frequency) (H^0 + λH1 + :::)( 0 + λ 1 + :::) = (E0 + λE1 + :::)( 0 + λ 1 + :::) I Expand and keep terms to first order in λ 2 H^0 0 + λH1 0 + λH0 1 = E0 0 + λE0 1 + λE1 0 + O(λ ) I Use H^0 0 = E0 0 (H^0 − E0) 1 = −(H^1 − E1) 0; Sternheimer equation Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 16 / 62 Sternheimer equation in TDDFT I (Weak) monochromatic perturbation v1(r; t) = λri cos(!t) I Expand time-dependent Kohn-Sham wavefunctions in powers of λ (0) (1) m(r; t) = exp(−i(m + λm )t)× n 1 o (0)(r) + λ[exp(i!t) (1)(r;!) + exp(−i!t) (1)(r; −!)] m 2 m m I Insert in time-dependent Kohn-Sham equation and keep terms up to first order in λ Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods July 19, 2011 17 / 62 Sternheimer equation in DFT I Frequency-dependent response (self-consistent solution!) h (0) i (1) (1) (0) H^ − j ± ! + iη (r; ±!) = H^ (±!) (r); with first-oder frequency-dependent perturbation Z ρ (r;!) Z H^ (1)(!) = v(r) + 1 d3r0 + f (r; r0;!)ρ (r0;!)d3r0 jr − r0j xc 1 and first-order density response occ.
Load more

Recommended publications
  • Accessing the Accuracy of Density Functional Theory Through Structure
    This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source are cited. Letter Cite This: J. Phys. Chem. Lett. 2019, 10, 4914−4919 pubs.acs.org/JPCL Accessing the Accuracy of Density Functional Theory through Structure and Dynamics of the Water−Air Interface † # ‡ # § ‡ § ∥ Tatsuhiko Ohto, , Mayank Dodia, , Jianhang Xu, Sho Imoto, Fujie Tang, Frederik Zysk, ∥ ⊥ ∇ ‡ § ‡ Thomas D. Kühne, Yasuteru Shigeta, , Mischa Bonn, Xifan Wu, and Yuki Nagata*, † Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan ‡ Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany § Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States ∥ Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Warburger Strasse 100, 33098 Paderborn, Germany ⊥ Graduate School of Pure and Applied Sciences, University of Tsukuba, Tennodai 1-1-1, Tsukuba, Ibaraki 305-8571, Japan ∇ Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan *S Supporting Information ABSTRACT: Density functional theory-based molecular dynamics simulations are increasingly being used for simulating aqueous interfaces. Nonetheless, the choice of the appropriate density functional, critically affecting the outcome of the simulation, has remained arbitrary. Here, we assess the performance of various exchange−correlation (XC) functionals, based on the metrics relevant to sum-frequency generation spectroscopy. The structure and dynamics of water at the water−air interface are governed by heterogeneous intermolecular interactions, thereby providing a critical benchmark for XC functionals.
    [Show full text]
  • Density Functional Theory
    Density Functional Approach Francesco Sottile Ecole Polytechnique, Palaiseau - France European Theoretical Spectroscopy Facility (ETSF) 22 October 2010 Density Functional Theory 1. Any observable of a quantum system can be obtained from the density of the system alone. < O >= O[n] Hohenberg, P. and W. Kohn, 1964, Phys. Rev. 136, B864 Density Functional Theory 1. Any observable of a quantum system can be obtained from the density of the system alone. < O >= O[n] 2. The density of an interacting-particles system can be calculated as the density of an auxiliary system of non-interacting particles. Hohenberg, P. and W. Kohn, 1964, Phys. Rev. 136, B864 Kohn, W. and L. Sham, 1965, Phys. Rev. 140, A1133 Density Functional ... Why ? Basic ideas of DFT Importance of the density Example: atom of Nitrogen (7 electron) 1. Any observable of a quantum Ψ(r1; ::; r7) 21 coordinates system can be obtained from 10 entries/coordinate ) 1021 entries the density of the system alone. 8 bytes/entry ) 8 · 1021 bytes 4:7 × 109 bytes/DVD ) 2 × 1012 DVDs 2. The density of an interacting-particles system can be calculated as the density of an auxiliary system of non-interacting particles. Density Functional ... Why ? Density Functional ... Why ? Density Functional ... Why ? Basic ideas of DFT Importance of the density Example: atom of Oxygen (8 electron) 1. Any (ground-state) observable Ψ(r1; ::; r8) 24 coordinates of a quantum system can be 24 obtained from the density of the 10 entries/coordinate ) 10 entries 8 bytes/entry ) 8 · 1024 bytes system alone. 5 · 109 bytes/DVD ) 1015 DVDs 2.
    [Show full text]
  • Octopus: a First-Principles Tool for Excited Electron-Ion Dynamics
    octopus: a first-principles tool for excited electron-ion dynamics. ¡ Miguel A. L. Marques a, Alberto Castro b ¡ a c, George F. Bertsch c and Angel Rubio a aDepartamento de F´ısica de Materiales, Facultad de Qu´ımicas, Universidad del Pa´ıs Vasco, Centro Mixto CSIC-UPV/EHU and Donostia International Physics Center (DIPC), 20080 San Sebastian,´ Spain bDepartamento de F´ısica Teorica,´ Universidad de Valladolid, E-47011 Valladolid, Spain cPhysics Department and Institute for Nuclear Theory, University of Washington, Seattle WA 98195 USA Abstract We present a computer package aimed at the simulation of the electron-ion dynamics of finite systems, both in one and three dimensions, under the influence of time-dependent electromagnetic fields. The electronic degrees of freedom are treated quantum mechani- cally within the time-dependent Kohn-Sham formalism, while the ions are handled classi- cally. All quantities are expanded in a regular mesh in real space, and the simulations are performed in real time. Although not optimized for that purpose, the program is also able to obtain static properties like ground-state geometries, or static polarizabilities. The method employed proved quite reliable and general, and has been successfully used to calculate linear and non-linear absorption spectra, harmonic spectra, laser induced fragmentation, etc. of a variety of systems, from small clusters to medium sized quantum dots. Key words: Electronic structure, time-dependent, density-functional theory, non-linear optics, response functions PACS: 33.20.-t, 78.67.-n, 82.53.-k PROGRAM SUMMARY Title of program: octopus Catalogue identifier: Program obtainable from: CPC Program Library, Queen’s University of Belfast, N.
    [Show full text]
  • 1-DFT Introduction
    MBPT and TDDFT Theory and Tools for Electronic-Optical Properties Calculations in Material Science Dott.ssa Letizia Chiodo Nano-bio Spectroscopy Group & ETSF - European Theoretical Spectroscopy Facility, Dipartemento de Física de Materiales, Facultad de Químicas, Universidad del País Vasco UPV/EHU, San Sebastián-Donostia, Spain Outline of the Lectures • Many Body Problem • DFT elements; examples • DFT drawbacks • excited properties: electronic and optical spectroscopies. elements of theory • Many Body Perturbation Theory: GW • codes, examples of GW calculations • Many Body Perturbation Theory: BSE • codes, examples of BSE calculations • Time Dependent DFT • codes, examples of TDDFT calculations • state of the art, open problems Main References theory • P. Hohenberg & W. Kohn, Phys. Rev. 136 (1964) B864; W. Kohn & L. J. Sham, Phys. Rev. 140 , A1133 (1965); (Nobel Prize in Chemistry1998) • Richard M. Martin, Electronic Structure: Basic Theory and Practical Methods, Cambridge University Press, 2004 • M. C. Payne, Rev. Mod. Phys.64 , 1045 (1992) • E. Runge and E.K.U. Gross, Phys. Rev. Lett. 52 (1984) 997 • M. A. L.Marques, C. A. Ullrich, F. Nogueira, A. Rubio, K. Burke, E. K. U. Gross, Time-Dependent Density Functional Theory. (Springer-Verlag, 2006). • L. Hedin, Phys. Rev. 139 , A796 (1965) • R.W. Godby, M. Schluter, L. J. Sham. Phys. Rev. B 37 , 10159 (1988) • G. Onida, L. Reining, A. Rubio, Rev. Mod. Phys. 74 , 601 (2002) codes & tutorials • Q-Espresso, http://www.pwscf.org/ • Abinit, http://www.abinit.org/ • Yambo, http://www.yambo-code.org • Octopus, http://www.tddft.org/programs/octopus more info at http://www.etsf.eu, www.nanobio.ehu.es Outline of the Lectures • Many Body Problem • DFT elements; examples • DFT drawbacks • excited properties: electronic and optical spectroscopies.
    [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]
  • Accelerating Performance and Scalability with NVIDIA Gpus on HPC Applications
    Accelerating Performance and Scalability with NVIDIA GPUs on HPC Applications Pak Lui The HPC Advisory Council Update • World-wide HPC non-profit organization • ~425 member companies / universities / organizations • Bridges the gap between HPC usage and its potential • Provides best practices and a support/development center • Explores future technologies and future developments • Leading edge solutions and technology demonstrations 2 HPC Advisory Council Members 3 HPC Advisory Council Centers HPC ADVISORY COUNCIL CENTERS HPCAC HQ SWISS (CSCS) CHINA AUSTIN 4 HPC Advisory Council HPC Center Dell™ PowerEdge™ Dell PowerVault MD3420 HPE Apollo 6000 HPE ProLiant SL230s HPE Cluster Platform R730 GPU Dell PowerVault MD3460 10-node cluster Gen8 3000SL 36-node cluster 4-node cluster 16-node cluster InfiniBand Storage (Lustre) Dell™ PowerEdge™ C6145 Dell™ PowerEdge™ R815 Dell™ PowerEdge™ Dell™ PowerEdge™ M610 InfiniBand-based 6-node cluster 11-node cluster R720xd/R720 32-node GPU 38-node cluster Storage (Lustre) cluster Dell™ PowerEdge™ C6100 4-node cluster 4-node GPU cluster 4-node GPU cluster 5 Exploring All Platforms / Technologies X86, Power, GPU, FPGA and ARM based Platforms x86 Power GPU FPGA ARM 6 HPC Training • HPC Training Center – CPUs – GPUs – Interconnects – Clustering – Storage – Cables – Programming – Applications • Network of Experts – Ask the experts 7 University Award Program • University award program – Universities / individuals are encouraged to submit proposals for advanced research • Selected proposal will be provided with: – Exclusive computation time on the HPC Advisory Council’s Compute Center – Invitation to present in one of the HPC Advisory Council’s worldwide workshops – Publication of the research results on the HPC Advisory Council website • 2010 award winner is Dr.
    [Show full text]
  • Application Profiling at the HPCAC High Performance Center Pak Lui 157 Applications Best Practices Published
    Best Practices: Application Profiling at the HPCAC High Performance Center Pak Lui 157 Applications Best Practices Published • Abaqus • COSMO • HPCC • Nekbone • RFD tNavigator • ABySS • CP2K • HPCG • NEMO • SNAP • AcuSolve • CPMD • HYCOM • NWChem • SPECFEM3D • Amber • Dacapo • ICON • Octopus • STAR-CCM+ • AMG • Desmond • Lattice QCD • OpenAtom • STAR-CD • AMR • DL-POLY • LAMMPS • OpenFOAM • VASP • ANSYS CFX • Eclipse • LS-DYNA • OpenMX • WRF • ANSYS Fluent • FLOW-3D • miniFE • OptiStruct • ANSYS Mechanical• GADGET-2 • MILC • PAM-CRASH / VPS • BQCD • Graph500 • MSC Nastran • PARATEC • BSMBench • GROMACS • MR Bayes • Pretty Fast Analysis • CAM-SE • Himeno • MM5 • PFLOTRAN • CCSM 4.0 • HIT3D • MPQC • Quantum ESPRESSO • CESM • HOOMD-blue • NAMD • RADIOSS For more information, visit: http://www.hpcadvisorycouncil.com/best_practices.php 2 35 Applications Installation Best Practices Published • Adaptive Mesh Refinement (AMR) • ESI PAM-CRASH / VPS 2013.1 • NEMO • Amber (for GPU/CUDA) • GADGET-2 • NWChem • Amber (for CPU) • GROMACS 5.1.2 • Octopus • ANSYS Fluent 15.0.7 • GROMACS 4.5.4 • OpenFOAM • ANSYS Fluent 17.1 • GROMACS 5.0.4 (GPU/CUDA) • OpenMX • BQCD • Himeno • PyFR • CASTEP 16.1 • HOOMD Blue • Quantum ESPRESSO 4.1.2 • CESM • LAMMPS • Quantum ESPRESSO 5.1.1 • CP2K • LAMMPS-KOKKOS • Quantum ESPRESSO 5.3.0 • CPMD • LS-DYNA • WRF 3.2.1 • DL-POLY 4 • MrBayes • WRF 3.8 • ESI PAM-CRASH 2015.1 • NAMD For more information, visit: http://www.hpcadvisorycouncil.com/subgroups_hpc_works.php 3 HPC Advisory Council HPC Center HPE Apollo 6000 HPE ProLiant
    [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 SoftwareESS2020 This paper was selected as Featured This paper was selected as Scilight 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 Electronic structure software The Journal of Chemical Physics 153, 070401 (2020); https://doi.org/10.1063/5.0023185 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]
  • 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]
  • Simulation Studies of the Mechanisms of Interaction Between Carbon Nanotubes and Amino Acids
    Simulation Studies of the Mechanisms of Interaction between Carbon Nanotubes and Amino Acids by George Bassem Botros Abadir B.Sc., Ain Shams University, Cairo, Egypt, 2000 M.Sc., Ain Shams University, Cairo, Egypt, 2005 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in The Faculty of Graduate Studies (Electrical and Computer Engineering) THE UNIVERSITY OF BRITISH COLUMBIA (Vancouver) August 2010 c George Bassem Botros Abadir 2010 Abstract In this thesis, molecular dynamics and ab initio density functional theory/non- equilibrium Green’s function simulations are used to study the interaction between carbon nanotubes and amino acids. Firstly, rules for the proper choice of the parameters used in these simulations are established. It is demonstrated how the improper choice of these parameters (particularly the basis set used in ab initio simulations) can lead to quantitatively and qualitatively erroneous conclusions regarding the bandgap of the nanotubes. It is then shown that the major forces responsible for amino-acid adsorp- tion on carbon nanotubes are van der Waals forces, and that hydrophobic interactions may accelerate the adsorption process, but are not necessary for it to occur. The mechanisms of interaction between carbon nanotubes and amino acids are elucidated. It is found that geometrical deformations do not play a major role in the sensing process, and that electrostatic inter- actions represent the major interaction mechanism between the tubes and amino acids. Fully metallic armchair tubes are found to be insensitive to various amino acids, while small-radius nanotubes are shown to be inade- quate for sensing in aqueous media, as their response to the motion of the atoms resulting from the immersion in water is comparable to that of an- alyte adsorption.
    [Show full text]
  • Newsletter 85
    Ψk Newsletter AB INITIO (FROM ELECTRONIC STRUCTURE) CALCULATION OF COMPLEX PROCESSES IN MATERIALS Number 85 February 2008 Editor: Z. (Dzidka) Szotek Sponsored by: UK’s CCP9 E-mail: [email protected] and ESF Psi-k Programme 1 Contents 1 Editorial 4 2 General News 5 2.1 Psi-k Portal and Psi-k Mailing List . ........ 5 3 Psi-k School/Workshop/Conference Announcements 6 3.1 THIRD Psi-k/NANOQUANTA SCHOOL& WORKSHOP . 6 3.2 Psi-kTraining-MarieCurieSummerSchool . ........ 11 4 General School/Workshop/Conference Announcements 12 4.1 2nd International Workshop on “Ab initio Description of Iron and Steel (ADIS2008): MagnetismandPhaseDiagrams” . 12 4.2 TheCAMDSummerSchool............................. 13 4.3 International Center for Materials Research (ICMR) Summer School on Multi- ferroicMaterialsandBeyond . 16 4.4 15th WIEN2khands-onworkshop . 18 4.5 Meeting on Molecular Dynamics for Non-Adiabatic Processes ........... 19 4.6 12th European Conference ”Physics of Magnetism” . ........... 20 4.7 SymposiumG,EMRSFallMeeting. 21 4.8 Summer School on Ab Initio Modelling in Solid State Chemistry - MSSC2008 . 23 4.9 Computational Molecular Science 2008 . ......... 25 5 General Job Announcements 26 6 Abstracts 55 7 New Book Announcements 73 8 SCIENTIFIC HIGHLIGHT OF THE MONTH: Local Self-Interaction Cor- rection of a Scattering Resonance: The Atom in Jellium Model 74 1 Introduction 74 2 Local Self-interaction Correction (LSIC) 76 2 3 LSDA solution of an atom in jellium 77 4 LSIC solution for an atom in jellium 78 5 Self-interaction correction of scattering resonances for Ce in jellium 79 6 Effective Medium Theory of Ce metal 82 7 Conclusions 84 3 1 Editorial As we are currently trying to change the format of the future Psi-k newsletters, the present newsletter in a way already reflects this transition period.
    [Show full text]
  • Quantum Chemical Calculations of NMR Parameters
    Quantum Chemical Calculations of NMR Parameters Tatyana Polenova University of Delaware Newark, DE Winter School on Biomolecular NMR January 20-25, 2008 Stowe, Vermont OUTLINE INTRODUCTION Relating NMR parameters to geometric and electronic structure Classical calculations of EFG tensors Molecular properties from quantum chemical calculations Quantum chemistry methods DENSITY FUNCTIONAL THEORY FOR CALCULATIONS OF NMR PARAMETERS Introduction to DFT Software Practical examples Tutorial RELATING NMR OBSERVABLES TO MOLECULAR STRUCTURE NMR Spectrum NMR Parameters Local geometry Chemical structure (reactivity) I. Calculation of experimental NMR parameters Find unique solution to CQ, Q, , , , , II. Theoretical prediction of fine structure constants from molecular geometry Classical electrostatic model (EFG)- only in simple ionic compounds Quantum mechanical calculations (Density Functional Theory) (EFG, CSA) ELECTRIC FIELD GRADIENT (EFG) TENSOR: POINT CHARGE MODEL EFG TENSOR IS DETERMINED BY THE COMBINED ELECTRONIC AND NUCLEAR WAVEFUNCTION, NO ANALYTICAL EXPRESSION IN THE GENERAL CASE THE SIMPLEST APPROXIMATION: CLASSICAL POINT CHARGE MODEL n Zie 4 V2,k = 3 Y2,k ()i,i i=1 di 5 ATOMS CONTRIBUTING TO THE EFG TENSOR ARE TREATED AS POINT CHARGES, THE RESULTING EFG TENSOR IS THE SUM WITH RESPECT TO ALL ATOMS VERY CRUDE MODEL, WORKS QUANTITATIVELY ONLY IN SIMPLEST IONIC SYSTEMS, BUT YIELDS QUALITATIVE TRENDS AND GENERAL UNDERSTANDING OF THE SYMMETRY AND MAGNITUDE OF THE EXPECTED TENSOR ELECTRIC FIELD GRADIENT (EFG) TENSOR: POINT CHARGE MODEL n Zie 4 V2,k = 3 Y2,k ()i,i i=1 di 5 Ze V = ; V = 0; V = 0 2,0 d 3 2,±1 2,±2 2Ze V = ; V = 0; V = 0 2,0 d 3 2,±1 2,±2 3 Ze V = ; V = 0; V = 0 2,0 2 d 3 2,±1 2,±2 V2,0 = 0; V2,±1 = 0; V2,±2 = 0 MOLECULAR PROPERTIES FROM QUANTUM CHEMICAL CALCULATIONS H = E See for example M.
    [Show full text]