Electronic Properties of Materials for Solar Cells: Which Ab Initio Approaches Can We Trust?

Total Page:16

File Type:pdf, Size:1020Kb

Electronic Properties of Materials for Solar Cells: Which Ab Initio Approaches Can We Trust? Electronic properties of materials for solar cells: Which ab initio approaches can we trust? Silvana Botti 1LSI, CNRS-CEA-Ecole´ Polytechnique, Palaiseau, France 2LPMCN, CNRS-Universite´ Lyon 1, France 3European Theoretical Spectroscopy Facility September 14, 2009 – CONARES, Helsinki Silvana Botti Electronic excitations in solar cells 1 / 57 Collaborators Ecole Polytechnique Universite´ Lyon 1 Julien Vidal, Lucia Reining Fabio Trani, Miguel Marques EDF Paris CEA Saclay Par¨ Olsson, J.-F. Guillemoles Fabien Bruneval Silvana Botti Electronic excitations in solar cells 2 / 57 European Theoretical Spectroscopy Facility (ETSF) The ETSF is a knowledge center for theoretical spectroscopy carrying out state-of-the-art research on theoretical and computational methods for studying electronic and optical properties of materials. The ETSF gathers the experience and know-how of more than 200 researchers in Europe and the United States. The ETSF offers its expertise to researchers, industry, and students in the form of collaborative projects, free scientific software and training. Proposals can be submitted at any moment! Further information: www.etsf.eu Silvana Botti Electronic excitations in solar cells 3 / 57 Outline 1 Thin-film photovoltaic materials 2 What can we calculate within standard DFT? 3 How to go beyond standard DFT? GW approaches! 4 How to compare with experiments? Silvana Botti Electronic excitations in solar cells 4 / 57 Thin-film photovoltaic materials Outline 1 Thin-film photovoltaic materials 2 What can we calculate within standard DFT? 3 How to go beyond standard DFT? GW approaches! 4 How to compare with experiments? Silvana Botti Electronic excitations in solar cells 5 / 57 Thin-film photovoltaic materials Present state of photovoltaic efficiency from National Renewable Energy Laboratory (USA) Silvana Botti Electronic excitations in solar cells 6 / 57 Thin-film photovoltaic materials CIGS solar cell Devices have to fulfill 2 functions: Photogeneration of electron-hole pairs Separation of charge carriers to generate a current Structure: Molybdenum back contact CIGS layer (p-type layer) CdS layer (n-type layer) Wurth¨ Elektronik GmbH & Co. ZnO:Al TCO contact Efficiency = 13 % Silvana Botti Electronic excitations in solar cells 7 / 57 Thin-film photovoltaic materials Modeling photovoltaic materials Objectives Predict accurate values for fundamental opto-electronical properties of materials Deal with complex materials (large unit cells, defects) Silvana Botti Electronic excitations in solar cells 8 / 57 What can we calculate within standard DFT? Outline 1 Thin-film photovoltaic materials 2 What can we calculate within standard DFT? 3 How to go beyond standard DFT? GW approaches! 4 How to compare with experiments? Silvana Botti Electronic excitations in solar cells 9 / 57 What can we calculate within standard DFT? Ground state densities vs potentials At the heart of density functional theory (DFT) Is there a 1-to-1 mapping between different external potentials v(r) and their corresponding ground state densities ρ(r)? Silvana Botti Electronic excitations in solar cells 10 / 57 What can we calculate within standard DFT? Density functional theory (DFT) If we can give a positive answer, then it can be proved that (i) all observable quantities of a quantum system are completely determined by the density. (ii) which means that the basic variable is no more the many-body wavefunction Ψ(fr)g but the electron density ρ(r). P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). You can find details in R. M. Dreizler and E.K.U. Gross, Density Functional Theory, Springer (Berlin, 1990). Silvana Botti Electronic excitations in solar cells 11 / 57 What can we calculate within standard DFT? Density functional theory (DFT) If we can give a positive answer, then it can be proved that (i) all observable quantities of a quantum system are completely determined by the density. (ii) which means that the basic variable is no more the many-body wavefunction Ψ(fr)g but the electron density ρ(r). P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). You can find details in R. M. Dreizler and E.K.U. Gross, Density Functional Theory, Springer (Berlin, 1990). Silvana Botti Electronic excitations in solar cells 11 / 57 What can we calculate within standard DFT? Density functional theory DFT in its standard form is a ground state theory: Structural parameters: lattice parameters, internal distortions are usually good in LDA or GGA Formation energies for defects calculated from total energies are often reliable ::: but ::: Kohn-Sham energies are not meant to reproduce quasiparticle band structures Kohn-Sham DOS is not meant to reproduce photoemission Silvana Botti Electronic excitations in solar cells 12 / 57 What can we calculate within standard DFT? Kohn-Sham band structure Kohn-Sham (KS) equations r2 − + v (r) 'KS (r) = "KS'KS (r) 2 KS i i i occ. X KS 2 ρ (r) = j'i (r) j i The KS states are not one-electron energy states for the quasi-electrons in the solid However, it is common to interpret the solutions of the Kohn-Sham equations as one-electron states Often one obtains good band dispersions but band gaps are systematically underestimated Silvana Botti Electronic excitations in solar cells 13 / 57 What can we calculate within standard DFT? Discontinuity in Vxc .6 .6 E = (E − E ) − (E − E ) ε ε g (N+1) (N) (N) (N−1) ε .6 1 1 KS KS ∆ = " (N + 1) − " (N) ;& N+1 N .6 ( ε JDS 1 .6 ( ε .6 JDS 1 DFT KS KS Eg = "N+1(N) − "N (N) N N 1HOHFWURQV 1HOHFWURQV DFT (N+1) (N) ∆xc = Eg − Eg = Vxc (r) − Vxc (r) Band gap error not due to LDA, but to the discontinuity in the exact Vxc L. J. Sham and M. Schluter,¨ PRL 51, 1888 (1983); PRB 32, 3883 (1985) J. P. Perdew and M. Levy, PRL 51, 1884 (1983) R. W. Godby, M. Schluter¨ and L. J. Sham, PRL 56, 2415 (1986) Silvana Botti Electronic excitations in solar cells 14 / 57 What can we calculate within standard DFT? An intuitive Picture for Absorption c unoccupied states occupied states v KS KS Independent particle KS picture using "i and 'i Silvana Botti Electronic excitations in solar cells 15 / 57 What can we calculate within standard DFT? An intuitive Picture for Absorption c unoccupied states occupied states v Photoemission process: hν − (Ekin + φ) = EN−1;v − EN;0 = −"v Silvana Botti Electronic excitations in solar cells 16 / 57 What can we calculate within standard DFT? An intuitive Picture for Absorption c unoccupied states occupied states v Optical absorption: electron-hole interaction (excitons) Silvana Botti Electronic excitations in solar cells 17 / 57 What can we calculate within standard DFT? Software supporting DFT Abinit EXCITING OpenMX ADF Fireball ORCA AIMPRO FSatom - list of codes ParaGauss Atomistix Toolkit GAMESS (UK) PLATO CADPAC GAMESS (US) PWscf (Quantum- CASTEP GAUSSIAN ESPRESSO) CPMD JAGUAR Q-Chem CRYSTAL06 MOLCAS SIESTA DACAPO MOLPRO Spartan DALTON MPQC S/PHI/nX deMon2K NRLMOL TURBOMOLE DFT++ NWChem VASP DMol3 OCTOPUS WIEN2k http://en.wikipedia.org/wiki/Density_functional_theory Silvana Botti Electronic excitations in solar cells 18 / 57 What can we calculate within standard DFT? The code ABINIT http://www.abinit.org ”First-principles computation of material properties : the ABINIT software project.” X. Gonze et al, Computational Materials Science 25, 478-492 (2002). ”A brief introduction to the ABINIT software package.” X. Gonze et al, Zeit. Kristallogr. 220, 558-562 (2005). Silvana Botti Electronic excitations in solar cells 19 / 57 What can we calculate within standard DFT? LDA Kohn-Sham energy gaps 8 AlN SrO diamond 6 ZnO,GaN,ZnS 4 ZnSe,CuBr AlAs,GaP,SiC,AlP,CdS Se,Cu2O MgO InP,GaAs,CdTe,AlSb 2 Si InN,Ge,GaSb,CdO CaO calculated gap (eV) InSb,P,InAs HgTe 0 :LDA experimental gap (eV) van Schilfgaarde, Kotani, and Faleev, PRL 96 (2006) Silvana Botti Electronic excitations in solar cells 20 / 57 What can we calculate within standard DFT? LDA Kohn-Sham energy gaps for CIS CuInS2 DFT-LDA exp. Eg -0.11 1.54 In-S 6.5 6.9 S s band 12.4 12.0 In 4 d band 14.6 18.2 CuInSe2 DFT-LDA exp. Eg -0.29 1.05 In-Se 5.8 6.5 Se s band 12.6 13.0 In 4 d band 14.7 18.0 www.abinit.org Silvana Botti Electronic excitations in solar cells 21 / 57 What can we calculate within standard DFT? Why do we need to go beyond standard DFT? For photovoltaic applications we are interested in evaluating quasiparticle band gap optical band gap defect energy levels optical absorption spectra All these quantities require going beyond standard DFT Silvana Botti Electronic excitations in solar cells 22 / 57 How to go beyond standard DFT? GW approaches! Outline 1 Thin-film photovoltaic materials 2 What can we calculate within standard DFT? 3 How to go beyond standard DFT? GW approaches! 4 How to compare with experiments? Silvana Botti Electronic excitations in solar cells 23 / 57 24 / 57 for the inclusion of ¡ ¡ ¢ £ ¢ £ £ ¡ ¡ ¡ ¡ ¢ £ ¢ ¢ ¢ ¡ £ £ £ £ ¢ ¡ ¡ ¡ ¡ ¢ ¢¤£ ¢ ¡ ¢ Electronic excitations in solar cells for quasi-particle properties ¢ ¢ ¢ ¡ ¡¢¡ ¡¢¡ ¡ ¢ GW Bethe-Salpeter equation electron-hole interaction ¢ ¢ ¢ ¢ ¡ ¢ ¢ ¢ ¡¥¡ ¡ ¨ Π In the many-body framework, wethese know problems: how to solve The first step can bethan substantially the more second, complicated so inGW the following we will focus on ¢ 0 ¢ ¢ ¨ ¡¢¡ ¡ ¡ ¢ ¢ ¨ ¨ ¢ ¢ ¢ ¢ ¡ Silvana Botti How to go beyond standard DFT? GW approaches! ¡ ¢ ¡ ¢ Solution How to go beyond standard DFT? GW approaches! Hedin’s equations Σ G=G Γ 0 GW +G 0 Σ Σ = G W G Γ W = v + vPW G)GG Γ Γ=1+(δΣ/δ P P = GGΓ L. Hedin, Phys. Rev. 139 (1965). Silvana Botti Electronic
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]
  • 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]
  • Robust Determination of the Chemical Potential in the Pole
    Robust determination of the chemical potential in the pole expansion and selected inversion method for solving Kohn-Sham density functional theory Weile Jia, and Lin Lin Citation: The Journal of Chemical Physics 147, 144107 (2017); doi: 10.1063/1.5000255 View online: http://dx.doi.org/10.1063/1.5000255 View Table of Contents: http://aip.scitation.org/toc/jcp/147/14 Published by the American Institute of Physics THE JOURNAL OF CHEMICAL PHYSICS 147, 144107 (2017) Robust determination of the chemical potential in the pole expansion and selected inversion method for solving Kohn-Sham density functional theory Weile Jia1,a) and Lin Lin1,2,b) 1Department of Mathematics, University of California, Berkeley, California 94720, USA 2Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA (Received 14 August 2017; accepted 24 September 2017; published online 11 October 2017) Fermi operator expansion (FOE) methods are powerful alternatives to diagonalization type methods for solving Kohn-Sham density functional theory (KSDFT). One example is the pole expansion and selected inversion (PEXSI) method, which approximates the Fermi operator by rational matrix functions and reduces the computational complexity to at most quadratic scaling for solving KSDFT. Unlike diagonalization type methods, the chemical potential often cannot be directly read off from the result of a single step of evaluation of the Fermi operator. Hence multiple evaluations are needed to be sequentially performed to compute the chemical potential to ensure the correct number of electrons within a given tolerance. This hinders the performance of FOE methods in practice. In this paper, we develop an efficient and robust strategy to determine the chemical potential in the context of the PEXSI method.
    [Show full text]
  • Camcasp 5.9 Alston J. Misquitta† and Anthony J
    CamCASP 5.9 Alston J. Misquittay and Anthony J. Stoneyy yDepartment of Physics and Astronomy, Queen Mary, University of London, 327 Mile End Road, London E1 4NS yy University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW September 1, 2015 Abstract CamCASP is a suite of programs designed to calculate molecular properties (multipoles and frequency- dependent polarizabilities) in single-site and distributed form, and interaction energies between pairs of molecules, and thence to construct atom–atom potentials. The CamCASP distribution also includes the programs Pfit, Casimir,Gdma 2.2, Cluster, and Process. Copyright c 2007–2014 Alston J. Misquitta and Anthony J. Stone Contents 1 Introduction 1 1.1 Authors . .1 1.2 Citations . .1 2 What’s new? 2 3 Outline of the capabilities of CamCASP and other programs 5 3.1 CamCASP limits . .7 4 Installation 7 4.1 Building CamCASP from source . .9 5 Using CamCASP 10 5.1 Workflows . 10 5.2 High-level scripts . 10 5.3 The runcamcasp.py script . 11 5.4 Low-level scripts . 13 6 Data conventions 13 7 CLUSTER: Detailed specification 14 7.1 Prologue . 15 7.2 Molecule definitions . 15 7.3 Geometry manipulations and other transformations . 16 7.4 Job specification . 18 7.5 Energy . 23 7.6 Crystal . 25 7.7 ORIENT ............................................... 26 7.8 Finally, . 29 8 Examples 29 8.1 A SAPT(DFT) calculation . 29 8.2 An example properties calculation . 30 8.3 Dispersion coefficients . 34 8.4 Using CLUSTER to obtain the dimer geometry . 36 9 CamCASP program specification 39 9.1 Global data .
    [Show full text]
  • The Impact of Density Functional Theory on Materials Research
    www.mrs.org/bulletin functional. This functional (i.e., a function whose argument is another function) de- scribes the complex kinetic and energetic interactions of an electron with other elec- Toward Computational trons. Although the form of this functional that would make the reformulation of the many-body Schrödinger equation exact is Materials Design: unknown, approximate functionals have proven highly successful in describing many material properties. Efficient algorithms devised for solving The Impact of the Kohn–Sham equations have been imple- mented in increasingly sophisticated codes, tremendously boosting the application of DFT methods. New doors are opening to in- Density Functional novative research on materials across phys- ics, chemistry, materials science, surface science, and nanotechnology, and extend- ing even to earth sciences and molecular Theory on Materials biology. The impact of this fascinating de- velopment has not been restricted to aca- demia, as DFT techniques also find application in many different areas of in- Research dustrial research. The development is so fast that many current applications could Jürgen Hafner, Christopher Wolverton, and not have been realized three years ago and were hardly dreamed of five years ago. Gerbrand Ceder, Guest Editors The articles collected in this issue of MRS Bulletin present a few of these suc- cess stories. However, even if the compu- Abstract tational tools necessary for performing The development of modern materials science has led to a growing need to complex quantum-mechanical calcula- tions relevant to real materials problems understand the phenomena determining the properties of materials and processes on are now readily available, designing a an atomistic level.
    [Show full text]
  • Practice: Quantum ESPRESSO I
    MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW FREE http://www.quantum-espresso.org PW-DFT, PP, PAW FREE http://www.abinit.org DFT PW, PP, Car-Parrinello FREE http://www.cpmd.org DFT PP, US-PP, PAW $3000 [moderate accuracy, fast] http://www.vasp.at DFT full-potential linearized augmented $500 plane-wave (FLAPW) [accurate, slow] http://www.wien2k.at Hartree-Fock, higher order correlated $3000 electron approaches http://www.gaussian.com 3 Quantum ESPRESSO 4 Quantum ESPRESSO Quantum ESPRESSO is an integrated suite of Open- Source computer codes for electronic-structure calculations and materials modeling at the nanoscale. It is based on density-functional theory, plane waves, and pseudopotentials. Core set of codes, plugins for more advanced tasks and third party packages Open initiative coordinated by the Quantum ESPRESSO Foundation, across Italy. Contributed to by developers across the world Regular hands-on workshops in Trieste, Italy Open-source code: FREE (unlike VASP...) 5 Performance Small jobs (a few atoms) can be run on single node Includes determining convergence parameters, lattice constants Can use OpenMP parallelization on multicore machines Large jobs (~10’s to ~100’s atoms) can run in parallel using MPI to 1000’s of cores Includes molecular dynamics, large geometry relaxation, phonons Parallel performance tied to BLAS/LAPACK (linear algebra routines) and 3D FFT (fast Fourier transform) New GPU-enabled version available 6 Usability Documented online:
    [Show full text]
  • The CECAM Electronic Structure Library and the Modular Software Development Paradigm
    The CECAM electronic structure library and the modular software development paradigm Cite as: J. Chem. Phys. 153, 024117 (2020); https://doi.org/10.1063/5.0012901 Submitted: 06 May 2020 . Accepted: 08 June 2020 . Published Online: 13 July 2020 Micael J. T. Oliveira , Nick Papior , Yann Pouillon , Volker Blum , Emilio Artacho , Damien Caliste , Fabiano Corsetti , Stefano de Gironcoli , Alin M. Elena , Alberto García , Víctor M. García-Suárez , Luigi Genovese , William P. Huhn , Georg Huhs , Sebastian Kokott , Emine Küçükbenli , Ask H. Larsen , Alfio Lazzaro , Irina V. Lebedeva , Yingzhou Li , David López- Durán , Pablo López-Tarifa , Martin Lüders , Miguel A. L. Marques , Jan Minar , Stephan Mohr , Arash A. Mostofi , Alan O’Cais , Mike C. Payne, Thomas Ruh, Daniel G. A. Smith , José M. Soler , David A. Strubbe , Nicolas Tancogne-Dejean , Dominic Tildesley, Marc Torrent , and Victor Wen-zhe Yu COLLECTIONS Paper published as part of the special topic on Electronic Structure Software Note: This article is part of the JCP Special Topic on Electronic Structure Software. This paper was selected as Featured ARTICLES YOU MAY BE INTERESTED IN Recent developments in the PySCF program package The Journal of Chemical Physics 153, 024109 (2020); https://doi.org/10.1063/5.0006074 An open-source coding paradigm for electronic structure calculations Scilight 2020, 291101 (2020); https://doi.org/10.1063/10.0001593 Siesta: Recent developments and applications The Journal of Chemical Physics 152, 204108 (2020); https://doi.org/10.1063/5.0005077 J. Chem. Phys. 153, 024117 (2020); https://doi.org/10.1063/5.0012901 153, 024117 © 2020 Author(s). The Journal ARTICLE of Chemical Physics scitation.org/journal/jcp The CECAM electronic structure library and the modular software development paradigm Cite as: J.
    [Show full text]
  • Kepler Gpus and NVIDIA's Life and Material Science
    LIFE AND MATERIAL SCIENCES Mark Berger; [email protected] Founded 1993 Invented GPU 1999 – Computer Graphics Visual Computing, Supercomputing, Cloud & Mobile Computing NVIDIA - Core Technologies and Brands GPU Mobile Cloud ® ® GeForce Tegra GRID Quadro® , Tesla® Accelerated Computing Multi-core plus Many-cores GPU Accelerator CPU Optimized for Many Optimized for Parallel Tasks Serial Tasks 3-10X+ Comp Thruput 7X Memory Bandwidth 5x Energy Efficiency How GPU Acceleration Works Application Code Compute-Intensive Functions Rest of Sequential 5% of Code CPU Code GPU CPU + GPUs : Two Year Heart Beat 32 Volta Stacked DRAM 16 Maxwell Unified Virtual Memory 8 Kepler Dynamic Parallelism 4 Fermi 2 FP64 DP GFLOPS GFLOPS per DP Watt 1 Tesla 0.5 CUDA 2008 2010 2012 2014 Kepler Features Make GPU Coding Easier Hyper-Q Dynamic Parallelism Speedup Legacy MPI Apps Less Back-Forth, Simpler Code FERMI 1 Work Queue CPU Fermi GPU CPU Kepler GPU KEPLER 32 Concurrent Work Queues Developer Momentum Continues to Grow 100M 430M CUDA –Capable GPUs CUDA-Capable GPUs 150K 1.6M CUDA Downloads CUDA Downloads 1 50 Supercomputer Supercomputers 60 640 University Courses University Courses 4,000 37,000 Academic Papers Academic Papers 2008 2013 Explosive Growth of GPU Accelerated Apps # of Apps Top Scientific Apps 200 61% Increase Molecular AMBER LAMMPS CHARMM NAMD Dynamics GROMACS DL_POLY 150 Quantum QMCPACK Gaussian 40% Increase Quantum Espresso NWChem Chemistry GAMESS-US VASP CAM-SE 100 Climate & COSMO NIM GEOS-5 Weather WRF Chroma GTS 50 Physics Denovo ENZO GTC MILC ANSYS Mechanical ANSYS Fluent 0 CAE MSC Nastran OpenFOAM 2010 2011 2012 SIMULIA Abaqus LS-DYNA Accelerated, In Development NVIDIA GPU Life Science Focus Molecular Dynamics: All codes are available AMBER, CHARMM, DESMOND, DL_POLY, GROMACS, LAMMPS, NAMD Great multi-GPU performance GPU codes: ACEMD, HOOMD-Blue Focus: scaling to large numbers of GPUs Quantum Chemistry: key codes ported or optimizing Active GPU acceleration projects: VASP, NWChem, Gaussian, GAMESS, ABINIT, Quantum Espresso, BigDFT, CP2K, GPAW, etc.
    [Show full text]
  • Quantum Chemistry (QC) on Gpus Feb
    Quantum Chemistry (QC) on GPUs Feb. 2, 2017 Overview of Life & Material Accelerated Apps MD: All key codes are GPU-accelerated QC: All key codes are ported or optimizing Great multi-GPU performance Focus on using GPU-accelerated math libraries, OpenACC directives Focus on dense (up to 16) GPU nodes &/or large # of GPU nodes GPU-accelerated and available today: ACEMD*, AMBER (PMEMD)*, BAND, CHARMM, DESMOND, ESPResso, ABINIT, ACES III, ADF, BigDFT, CP2K, GAMESS, GAMESS- Folding@Home, GPUgrid.net, GROMACS, HALMD, HOOMD-Blue*, UK, GPAW, LATTE, LSDalton, LSMS, MOLCAS, MOPAC2012, LAMMPS, Lattice Microbes*, mdcore, MELD, miniMD, NAMD, NWChem, OCTOPUS*, PEtot, QUICK, Q-Chem, QMCPack, OpenMM, PolyFTS, SOP-GPU* & more Quantum Espresso/PWscf, QUICK, TeraChem* Active GPU acceleration projects: CASTEP, GAMESS, Gaussian, ONETEP, Quantum Supercharger Library*, VASP & more green* = application where >90% of the workload is on GPU 2 MD vs. QC on GPUs “Classical” Molecular Dynamics Quantum Chemistry (MO, PW, DFT, Semi-Emp) Simulates positions of atoms over time; Calculates electronic properties; chemical-biological or ground state, excited states, spectral properties, chemical-material behaviors making/breaking bonds, physical properties Forces calculated from simple empirical formulas Forces derived from electron wave function (bond rearrangement generally forbidden) (bond rearrangement OK, e.g., bond energies) Up to millions of atoms Up to a few thousand atoms Solvent included without difficulty Generally in a vacuum but if needed, solvent treated classically
    [Show full text]
  • Jaguar 5.5 User Manual Copyright © 2003 Schrödinger, L.L.C
    Jaguar 5.5 User Manual Copyright © 2003 Schrödinger, L.L.C. All rights reserved. Schrödinger, FirstDiscovery, Glide, Impact, Jaguar, Liaison, LigPrep, Maestro, Prime, QSite, and QikProp are trademarks of Schrödinger, L.L.C. MacroModel is a registered trademark of Schrödinger, L.L.C. To the maximum extent permitted by applicable law, this publication is provided “as is” without warranty of any kind. This publication may contain trademarks of other companies. October 2003 Contents Chapter 1: Introduction.......................................................................................1 1.1 Conventions Used in This Manual.......................................................................2 1.2 Citing Jaguar in Publications ...............................................................................3 Chapter 2: The Maestro Graphical User Interface...........................................5 2.1 Starting Maestro...................................................................................................5 2.2 The Maestro Main Window .................................................................................7 2.3 Maestro Projects ..................................................................................................7 2.4 Building a Structure.............................................................................................9 2.5 Atom Selection ..................................................................................................10 2.6 Toolbar Controls ................................................................................................11
    [Show full text]
  • Thomas–Fermi–Dirac–Von Weizsäcker Models in Finite Systems Garnet Kin-Lic Chan, Aron J
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Caltech Authors - Main Thomas–Fermi–Dirac–von Weizsäcker models in finite systems Garnet Kin-Lic Chan, Aron J. Cohen, and Nicholas C. Handy Citation: 114, (2001); doi: 10.1063/1.1321308 View online: http://dx.doi.org/10.1063/1.1321308 View Table of Contents: http://aip.scitation.org/toc/jcp/114/2 Published by the American Institute of Physics JOURNAL OF CHEMICAL PHYSICS VOLUME 114, NUMBER 2 8 JANUARY 2001 Thomas–Fermi–Dirac–von Weizsa¨cker models in finite systems Garnet Kin-Lic Chan,a) Aron J. Cohen, and Nicholas C. Handy Department of Chemistry, Lensfield Road, Cambridge CB2 1EW, United Kingdom ͑Received 7 June 2000; accepted 8 September 2000͒ To gain an understanding of the variational behavior of kinetic energy functionals, we perform a numerical study of the Thomas–Fermi–Dirac–von Weizsa¨cker theory in finite systems. A general purpose Gaussian-based code is constructed to perform energy and geometry optimizations on polyatomic systems to high accuracy. We carry out benchmark studies on atomic and diatomic systems. Our results indicate that the Thomas–Fermi–Dirac–von Weizsa¨cker theory can give an approximate description of matter, with atomic energies, binding energies, and bond lengths of the correct order of magnitude, though not to the accuracy required of a qualitative chemical theory. We discuss the implications for the development of new kinetic functionals. © 2001 American Institute of Physics. ͓DOI: 10.1063/1.1321308͔ I. INTRODUCTION
    [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]