DOCSLIB.ORG

How to Do Simple Calculations with Quantum ESPRESSO

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
How to Do Simple Calculations with Quantum ESPRESSO How To Do Simple Calculations With Quantum ESPRESSO Shobhana Narasimhan Theoretical Sciences Unit JNCASR, Bangalore [email protected] I. About The Quantum ESPRESSO Distribution 2 Shobhana Narasimhan, JNCASR Quantum ESPRESSO • www.quantum-espresso.org 3 Shobhana Narasimhan, JNCASR The Quantum ESPRESSO Software Distribution P. Giannozzi 4 Shobhana Narasimhan, JNCASR Why “Quantum ESPRESSO”?! 5 Shobhana Narasimhan, JNCASR Licence for Quantum ESPRESSO P. Giannozzi 6 Shobhana Narasimhan, JNCASR Quantum ESPRESSO: Organization P. Giannozzi 7 Shobhana Narasimhan, JNCASR Quantum ESPRESSO as a distribution P. Giannozzi 8 Shobhana Narasimhan, JNCASR Quantum ESPRESSO as a distribution OTHER PACKAGES WANNIER90: Maximally localized Wannier functions Pwcond: Ballistic conductance WanT: Coherent Transport from Maximally Localized Wannier Functions Xspectra: Calculation of x-ray near edge absorption spectra GIPAW: EPR and NMR Chemical Shifts Coming Soon: GWW: GW Band Structure with Ultralocalized Wannier Fns. TD-DFT: Time-Dependent Density Functional Pert. Theory 9 Shobhana Narasimhan, JNCASR What Can Quantum ESPRESSO Do? • Both point and k-point calculations. • Both insulators and metals, with smearing. • Any crystal structure or supercell form. • Norm conserving pseudopotentials, ultrasoft PPs, PAW. • LDA, GGA, DFT+U, hybrid functionals, exact exchange, meta GGA, van der Waals corrected functionals. • Spin polarized calculations, non-collinear magnetism, spin-orbit interactions. • Nudged elastic band to find saddle points. 10 Shobhana Narasimhan, JNCASR II. Doing a “Total Energy” Calculation with the PWscf Package of QE: The SCF Loop 11 Shobhana Narasimhan, JNCASR The Kohn-Sham problem Want to solve the Kohn-Sham equations: 1 2 V (r) V [n(r)] V [n(r)] (r) (r) 2 nuc H XC i i i H Note that self-consistent solution necessary, as H depends on solution: { i} n(r) H Convention: e me 1 12 Shobhana Narasimhan, JNCASR Self-consistent Iterative Solution Vnuc known/constructed How to solve the Initial guess n(r) Kohn-Sham eqns. for a set of fixed Generate Calculate VH[n] & VXC[n] new nuclear (ionic) n(r) Veff(r)= Vnuc(r) + VH (r) + VXC (r) positions. 2 H i(r) = [-1/2 + Veff(r)] i(r) = i i(r) 2 Calculate new n(r) = i| i(r)| No Self-consistent? Yes Problem solved! Can now calculate energy, forces, etc. 13 Shobhana Narasimhan, JNCASR Plane Waves & Periodic Systems • For a periodic system: 1 where G = reciprocal (r) c ei(k G) r k k,G lattice vector G • The plane waves that appear in this expansion can be represented as a grid in k-space: k y • Only true for periodic systems that grid is discrete. kx • In principle, still need infinite number of plane waves. 14 Shobhana Narasimhan, JNCASR Truncating the Plane Wave Expansion • In practice, the contribution from higher Fourier components (large |k+G|) is small. • So truncate the expansion at some value of |k+G|. • Traditional to express this cut-off in energy units: ky E=Ecut 2 | k G |2 E 2m cut kx Input parameter ecutwfc 15 Shobhana Narasimhan, JNCASR Checking Convergence wrt ecutwfc • Must always check. • Monotonic (variational). 16 Shobhana Narasimhan, JNCASR Step 0: Defining the (periodic) system Namelist ‘SYSTEM’ Shobhana Narasimhan, JNCASR How to Specify the System • All periodic systems can be specified by a Bravais Lattice and an atomic basis. + = 18 Shobhana Narasimhan, JNCASR How to Specify the Bravais Lattice / Unit Cell Input parameter ibrav - Gives the type of Bravais lattice (SC, BCC, Hex, etc.) Input parameters {celldm(i)} a - Give the lengths [& 2 directions, if necessary] of a1 the lattice vectors a1, a2, a3 • Note that one can choose a non-primitive unit cell (e.g., 4 atom SC cell for FCC structure). 19 Shobhana Narasimhan, JNCASR Atoms Within Unit Cell – How many, where? Input parameter nat - Number of atoms in the unit cell Input parameter ntyp - Number of types of atoms FIELD ATOMIC_POSITIONS - Initial positions of atoms (may vary when “relax” done). -Can choose to give in units of lattice vectors (“crystal”) or in Cartesian units (“alat” or “bohr” or “angstrom”) 20 Shobhana Narasimhan, JNCASR What if the system is not periodic? • Example 1: Want to study properties of a system with a surface. Surface atom • Presence of surface No periodicity along z. Bulk atom z x 21 Shobhana Narasimhan, JNCASR What if the system is not periodic? • Example 1: Want to study properties of a system with a surface. • Presence of surface No periodicity along z. • Use a supercell: artificial periodicity along z by repeating slabs separated by vacuum. • Have to check convergence z w.r.t. slab thickness & vacuum thickness. x 22 Shobhana Narasimhan, JNCASR What if the system is not periodic? • Example 2: Want to study properties of a nanowire. z y • Example 3: Want to study properties of a cluster 23 Shobhana Narasimhan, JNCASR What if the system is not periodic? • Example 2: Want to study properties of a nanowire introduce artificial periodicity along y & z. • Example 3: Want to study properties of a cluster introduce artificial periodicity along x, y & z. z 24 Shobhana Narasimhan, JNCASR y What if the system is not periodic? • Example 4: Want to study a system with a defect, e.g., a vacancy or impurity: 25 Shobhana Narasimhan, JNCASR What if the system is not periodic? • Example 4: Want to study a system with a defect, e.g., a vacancy or impurity: 26 Shobhana Narasimhan, JNCASR What if the system is not periodic? • Example 5: Want to study an amorphous or quasicrystalline system. 27 Shobhana Narasimhan, JNCASR What if the system is not periodic? • Example 5: Want to study an amorphous or quasicrystalline system: approximate by a periodic system (with large unit cell). 28 Shobhana Narasimhan, JNCASR Artificially Periodic Systems Large Unit Cells • Note: In all these cases, to minimize the effects of the artificially introduced periodicity, need a large unit cell. • Long a1, a2, a3 (primitive lattice vectors) • Short b1, b2, b3 (primitive reciprocal lattice vectors) • Many G‟s will fall within Ecut sphere! 29 Shobhana Narasimhan, JNCASR Step 1: Obtaining Vnuc Vnuc known/constructed Initial guess n(r) Generate Calculate VH[n] & VXC[n] new n(r) Veff(r)= Vnuc(r) + VH (r) + VXC (r) 2 H i(r) = [-1/2 + Veff(r)] i(r) = i i(r) 2 Calculate new n(r) = i| i(r)| No Self-consistent? Yes Problem solved! Can now calculate energy, forces, etc. 30 Shobhana Narasimhan, JNCASR Nuclear Potential • Electrons experience a Coulomb potential due to the nuclei. • This has a known and simple form: Z V nuc r • But this leads to computational problems! 31 Shobhana Narasimhan, JNCASR Problem for Plane-Wave Basis Core wavefunctions: Valence wavefunctions: sharply peaked near lots of wiggles near nucleus. nucleus. High Fourier components present i.e., need large Ecut 32 Shobhana Narasimhan, JNCASR Solutions for Plane-Wave Basis Core wavefunctions: Valence wavefunctions: sharply peaked near lots of wiggles near nucleus. nucleus. High Fourier components present i.e., need large Ecut Don‟t solve for the Remove wiggles from core electrons! valence electrons. 33 Shobhana Narasimhan, JNCASR Pseudopotentials • Replace nuclear potential by pseudopotential • This is a numerical trick that solves these problems • There are different kinds of pseudopotentials (Norm conserving pseudopotentials, ultrasoft pseudopotentials, etc.) • Which kind you use depends on the element. 34 Shobhana Narasimhan, JNCASR Pseudopotentials for Quantum Espresso - 1 • Go to http://www.quantum-espresso.org; Click on “PSEUDO” 35 Shobhana Narasimhan, JNCASR Pseudopotentials for Quantum Espresso - 2 • Click on element for which pseudopotential wanted. 36 Shobhana Narasimhan, JNCASR Pseudopotentials for Quantum-ESPRESSO Pseudopotential‟s name gives information about : • type of exchange- correlation functional • type of pseudopotential • e.g.: 37 Shobhana Narasimhan, JNCASR Element & Vion for Quantum-ESPRESSO e.g, for calculation on BaTiO3: ATOMIC_SPECIES Ba 137.327 Ba.pbe-nsp-van.UPF Ti 47.867 Ti.pbe-sp-van_ak.UPF O 15.999 O.pbe-van_ak.UPF • ecutwfc, ecutrho depend on type of pseudopotentials used (should test). • When using ultrasoft pseudopotentials, set ecutrho = 8-12 ecutwfc !! 38 Shobhana Narasimhan, JNCASR Element & Vion for Quantum-ESPRESSO • Should have same exchange-correlation functional for all pseudopotentials. oops! input output 39 Shobhana Narasimhan, JNCASR Step 2: Initial Guess for n(r) Vion known/constructed Initial guess n(r) Generate Calculate VH[n] & VXC[n] new n(r) Veff(r)= Vion(r) + VH (r) + VXC (r) 2 H i(r) = [-1/2 + Veff(r)] i(r) = i i(r) 2 Calculate new n(r) = i| i(r)| No Self-consistent? Yes Problem solved! Can now calculate energy, forces, etc. 40 Shobhana Narasimhan, JNCASR Starting Wavefunctions The closer your starting wavefunction is to the true wavefunction (which, of course, is something you don‟t necessarily know to start with!), the fewer the scf iterations needed. startingwfc ‘atomic’ Superposition of atomic orbitals ‘random’ ‘file’ “The beginning is the most important part of the work” - Plato 41 Shobhana Narasimhan, JNCASR Steps 3 & 4: Effective Potential Vion known/constructed Initial guess n(r) or i(r) Note that type of exchange- Generate Calculate VH[n] & VXC[n] correlation chosen new while specifying n(r) V (r)= V (r) + V (r) + V (r) eff nuc H XC pseudopotential 2 H i(r) = [-1/2 + Veff(r)] i(r) = i i(r) 2 Calculate new n(r) = i| i(r)| No Self-consistent? Yes Problem solved! Can now calculate energy, forces, etc. 42 Shobhana Narasimhan, JNCASR Exchange-Correlation Potential • VXC EXC/ n contains all the many-body information. • Known [numerically, from Quantum Monte Carlo ; various analytical approximations] for homogeneous electron gas. • Local Density Approximation: HOM Exc[n] = n(r) Vxc [n(r)] dr -surprisingly successful! Replace pz (in name of pseudopotential) • Generalized Gradient Approximation(s): Include terms involving gradients of n(r) pw91, pbe (in name of pseudopotential) 43 Shobhana Narasimhan, JNCASR Step 5: Diagonalization Vion known/constructed Initial guess n(r) Generate Calculate VH[n] & VXC[n] new n(r) Veff(r)= Vion(r) + VH (r) + VXC (r) Expensive! 2 H i(r) = [-1/2 + Veff(r)] i(r) = i i(r) 2 Calculate new n(r) = i| i(r)| No Self-consistent? Yes Problem solved! Can now calculate energy, forces, etc.
Load more

Recommended publications
  • GPAW, Gpus, and LUMI
    GPAW, GPUs, and LUMI Martti Louhivuori, CSC - IT Center for Science Jussi Enkovaara GPAW 2021: Users and Developers Meeting, 2021-06-01 Outline LUMI supercomputer Brief history of GPAW with GPUs GPUs and DFT Current status Roadmap LUMI - EuroHPC system of the North Pre-exascale system with AMD CPUs and GPUs ~ 550 Pflop/s performance Half of the resources dedicated to consortium members Programming for LUMI Finland, Belgium, Czechia, MPI between nodes / GPUs Denmark, Estonia, Iceland, HIP and OpenMP for GPUs Norway, Poland, Sweden, and how to use Python with AMD Switzerland GPUs? https://www.lumi-supercomputer.eu GPAW and GPUs: history (1/2) Early proof-of-concept implementation for NVIDIA GPUs in 2012 ground state DFT and real-time TD-DFT with finite-difference basis separate version for RPA with plane-waves Hakala et al. in "Electronic Structure Calculations on Graphics Processing Units", Wiley (2016), https://doi.org/10.1002/9781118670712 PyCUDA, cuBLAS, cuFFT, custom CUDA kernels Promising performance with factor of 4-8 speedup in best cases (CPU node vs. GPU node) GPAW and GPUs: history (2/2) Code base diverged from the main branch quite a bit proof-of-concept implementation had lots of quick and dirty hacks fixes and features were pulled from other branches and patches no proper unit tests for GPU functionality active development stopped soon after publications Before development re-started, code didn't even work anymore on modern GPUs without applying a few small patches Lesson learned: try to always get new functionality to the
    [Show full text]
  • Free and Open Source Software for Computational Chemistry Education
    Free and Open Source Software for Computational Chemistry Education Susi Lehtola∗,y and Antti J. Karttunenz yMolecular Sciences Software Institute, Blacksburg, Virginia 24061, United States zDepartment of Chemistry and Materials Science, Aalto University, Espoo, Finland E-mail: [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]
  • D6.1 Report on the Deployment of the Max Demonstrators and Feedback to WP1-5
    Ref. Ares(2020)2820381 - 31/05/2020 HORIZON2020 European Centre of Excellence ​ Deliverable D6.1 Report on the deployment of the MaX Demonstrators and feedback to WP1-5 D6.1 Report on the deployment of the MaX Demonstrators and feedback to WP1-5 Pablo Ordejón, Uliana Alekseeva, Stefano Baroni, Riccardo Bertossa, Miki Bonacci, Pietro Bonfà, Claudia Cardoso, Carlo Cavazzoni, Vladimir Dikan, Stefano de Gironcoli, Andrea Ferretti, Alberto García, Luigi Genovese, Federico Grasselli, Anton Kozhevnikov, Deborah Prezzi, Davide Sangalli, Joost VandeVondele, Daniele Varsano, Daniel Wortmann Due date of deliverable: 31/05/2020 Actual submission date: 31/05/2020 Final version: 31/05/2020 Lead beneficiary: ICN2 (participant number 3) Dissemination level: PU - Public www.max-centre.eu 1 HORIZON2020 European Centre of Excellence ​ Deliverable D6.1 Report on the deployment of the MaX Demonstrators and feedback to WP1-5 Document information Project acronym: MaX Project full title: Materials Design at the Exascale Research Action Project type: European Centre of Excellence in materials modelling, simulations and design EC Grant agreement no.: 824143 Project starting / end date: 01/12/2018 (month 1) / 30/11/2021 (month 36) Website: www.max-centre.eu Deliverable No.: D6.1 Authors: P. Ordejón, U. Alekseeva, S. Baroni, R. Bertossa, M. ​ Bonacci, P. Bonfà, C. Cardoso, C. Cavazzoni, V. Dikan, S. de Gironcoli, A. Ferretti, A. García, L. Genovese, F. Grasselli, A. Kozhevnikov, D. Prezzi, D. Sangalli, J. VandeVondele, D. Varsano, D. Wortmann To be cited as: Ordejón, et al., (2020): Report on the deployment of the MaX Demonstrators and feedback to WP1-5. Deliverable D6.1 of the H2020 project MaX (final version as of 31/05/2020).
    [Show full text]
  • Density Functional Theory and Nuclear Quantum Effects
    1 Grassmann manifold, gauge, and quantum chemistry Lin Lin Department of Mathematics, UC Berkeley; Lawrence Berkeley National Laboratory Flatiron Institute, April 2019 2 Stiefel manifold • × > ℂ • Stiefel manifold (1935): orthogonal vectors in × ℂ , = | = ∗ ∈ ℂ 3 Grassmann manifold • Grassmann manifold (1848): set of dimensional subspace in ℂ , = , / : dimensional unitary group (i.e. the set of unitary matrices in × ) ℂ • Any point in , can be viewed as a representation of a point in , 4 Example (in optimization) min ( ) . = ∗ • = , . invariant to the choice of basis Optimization on a Grassmann manifold. ∈ ⇒ • The representation is called the gauge in quantum physics / chemistry. Gauge invariant. 5 Example (in optimization) min ( ) . = ∗ • Otherwise, if ( ) is gauge dependent Optimization on a Stiefel manifold. ⇒ • [Edelman, Arias, Smith, SIMAX 1998] many works following 6 This talk: Quantum chemistry • Time-dependent density functional theory (TDDFT) • Kohn-Sham density functional theory (KSDFT) • Localization. • Gauge choice is the key • Used in community software packages: Quantum ESPRESSO, Wannier90, Octopus, PWDFT etc 7 Time-dependent density functional theory Parallel transport gauge 8 Time-dependent density functional theory [Runge-Gross, PRL 1984] (spin omitted, : number of electrons) , = , , , = 1, … , , , = ( ,) ( , ) ′ ∗ ′ � Hamiltonian =1 1 , = + ( ) + [ ] 2 − Δ 9 Density matrix • Key quantity throughout the talk × • (t) = , … , . : number of discretized grid points to represent each Ψ 1 ∈ ℂ = ∗ Ψ Ψ = 10 Gauge invariance • Unitary rotation = , = ∗ = Φ Ψ= ( ) = ∗ ∗ ∗ ∗ ΦΦ Ψ Ψ ΨΨ • Propagation of the density matrix, von Neumann equation (quantum Liouville equation) = , = , ( ) , − 11 Gauge choice affects the time step Extreme case (boring Hamiltonian) ≡ 0 Initial condition 0 = 0 0 Exact solution = (0) Small time step − = (0) (Formally) arbitrarily large time step 12 P v.s.
    [Show full text]
  • Arxiv:2005.05756V2
    “This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article appeared in Oliveira, M.J.T. [et al.]. The CECAM electronic structure library and the modular software development paradigm. "The Journal of Chemical Physics", 12020, vol. 153, núm. 2, and may be found at https://aip.scitation.org/doi/10.1063/5.0012901. The CECAM Electronic Structure Library and the modular software development paradigm Micael J. T. Oliveira,1, a) Nick Papior,2, b) Yann Pouillon,3, 4, c) Volker Blum,5, 6 Emilio Artacho,7, 8, 9 Damien Caliste,10 Fabiano Corsetti,11, 12 Stefano de Gironcoli,13 Alin M. Elena,14 Alberto Garc´ıa,15 V´ıctor M. Garc´ıa-Su´arez,16 Luigi Genovese,10 William P. Huhn,5 Georg Huhs,17 Sebastian Kokott,18 Emine K¨u¸c¨ukbenli,13, 19 Ask H. Larsen,20, 4 Alfio Lazzaro,21 Irina V. Lebedeva,22 Yingzhou Li,23 David L´opez-Dur´an,22 Pablo L´opez-Tarifa,24 Martin L¨uders,1, 14 Miguel A. L. Marques,25 Jan Minar,26 Stephan Mohr,17 Arash A. Mostofi,11 Alan O'Cais,27 Mike C. Payne,9 Thomas Ruh,28 Daniel G. A. Smith,29 Jos´eM. Soler,30 David A. Strubbe,31 Nicolas Tancogne-Dejean,1 Dominic Tildesley,32 Marc Torrent,33, 34 and Victor Wen-zhe Yu5 1)Max Planck Institute for the Structure and Dynamics of Matter, D-22761 Hamburg, Germany 2)DTU Computing Center, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark 3)Departamento CITIMAC, Universidad de Cantabria, Santander, Spain 4)Simune Atomistics, 20018 San Sebasti´an,Spain 5)Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC 27708, USA 6)Department of Chemistry, Duke University, Durham, NC 27708, USA 7)CIC Nanogune BRTA and DIPC, 20018 San Sebasti´an,Spain 8)Ikerbasque, Basque Foundation for Science, 48011 Bilbao, Spain 9)Theory of Condensed Matter, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, United Kingdom 10)Department of Physics, IRIG, Univ.
    [Show full text]
  • CHEM:4480 Introduction to Molecular Modeling
    CHEM:4480:001 Introduction to Molecular Modeling Fall 2018 Instructor Dr. Sara E. Mason Office: W339 CB, Phone: (319) 335-2761, Email: [email protected] Lecture 11:00 AM - 12:15 PM, Tuesday/Thursday. Primary Location: W228 CB. Please note that we will often meet in one of the Chemistry Computer Labs. Notice about class being held outside of the primary location will be posted to the course ICON site by 24 hours prior to the start of class. Office Tuesday 12:30-2 PM (EXCEPT for DoC faculty meeting days 9/4, 10/2, 10/9, 11/6, Hour and 12/4), Thursday 12:30-2, or by appointment Department Dr. James Gloer DEO Administrative Office: E331 CB Administrative Phone: (319) 335-1350 Administrative Email: [email protected] Text Instead of a course text, assigned readings will be given throughout the semester in forms such as journal articles or library reserves. Website http://icon.uiowa.edu Course To develop practical skills associated with chemical modeling based on quantum me- Objectives chanics and computational chemistry such as working with shell commands, math- ematical software, and modeling packages. We will also work to develop a basic understanding of quantum mechanics, approximate methods, and electronic struc- ture. Technical computing, pseudopotential generation, density functional theory software, (along with other optional packages) will be used hands-on to gain experi- ence in data manipulation and modeling. We will use both a traditional classroom setting and hands-on learning in computer labs to achieve course goals. Course Topics • Introduction: It’s a quantum world! • Dirac notation and matrix mechanics (a.k.a.
    [Show full text]
  • Trends in Atomistic Simulation Software Usage [1.3]
    A LiveCoMS Perpetual Review Trends in atomistic simulation software usage [1.3] Leopold Talirz1,2,3*, Luca M. Ghiringhelli4, Berend Smit1,3 1Laboratory of Molecular Simulation (LSMO), Institut des Sciences et Ingenierie Chimiques, Valais, École Polytechnique Fédérale de Lausanne, CH-1951 Sion, Switzerland; 2Theory and Simulation of Materials (THEOS), Faculté des Sciences et Techniques de l’Ingénieur, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland; 3National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland; 4The NOMAD Laboratory at the Fritz Haber Institute of the Max Planck Society and Humboldt University, Berlin, Germany This LiveCoMS document is Abstract Driven by the unprecedented computational power available to scientific research, the maintained online on GitHub at https: use of computers in solid-state physics, chemistry and materials science has been on a continuous //github.com/ltalirz/ rise. This review focuses on the software used for the simulation of matter at the atomic scale. We livecoms-atomistic-software; provide a comprehensive overview of major codes in the field, and analyze how citations to these to provide feedback, suggestions, or help codes in the academic literature have evolved since 2010. An interactive version of the underlying improve it, please visit the data set is available at https://atomistic.software. GitHub repository and participate via the issue tracker. This version dated August *For correspondence: 30, 2021 [email protected] (LT) 1 Introduction Gaussian [2], were already released in the 1970s, followed Scientists today have unprecedented access to computa- by force-field codes, such as GROMOS [3], and periodic tional power.
    [Show full text]
  • Porting the DFT Code CASTEP to Gpgpus
    Porting the DFT code CASTEP to GPGPUs Toni Collis [email protected] EPCC, University of Edinburgh CASTEP and GPGPUs Outline • Why are we interested in CASTEP and Density Functional Theory codes. • Brief introduction to CASTEP underlying computational problems. • The OpenACC implementation http://www.nu-fuse.com CASTEP: a DFT code • CASTEP is a commercial and academic software package • Capable of Density Functional Theory (DFT) and plane wave basis set calculations. • Calculates the structure and motions of materials by the use of electronic structure (atom positions are dictated by their electrons). • Modern CASTEP is a re-write of the original serial code, developed by Universities of York, Durham, St. Andrews, Cambridge and Rutherford Labs http://www.nu-fuse.com CASTEP: a DFT code • DFT/ab initio software packages are one of the largest users of HECToR (UK national supercomputing service, based at University of Edinburgh). • Codes such as CASTEP, VASP and CP2K. All involve solving a Hamiltonian to explain the electronic structure. • DFT codes are becoming more complex and with more functionality. http://www.nu-fuse.com HECToR • UK National HPC Service • Currently 30- cabinet Cray XE6 system – 90,112 cores • Each node has – 2×16-core AMD Opterons (2.3GHz Interlagos) – 32 GB memory • Peak of over 800 TF and 90 TB of memory http://www.nu-fuse.com HECToR usage statistics Phase 3 statistics (Nov 2011 - Apr 2013) Ab initio codes (VASP, CP2K, CASTEP, ONETEP, NWChem, Quantum Espresso, GAMESS-US, SIESTA, GAMESS-UK, MOLPRO) GS2NEMO ChemShell 2%2% SENGA2% 3% UM Others 4% 34% MITgcm 4% CASTEP 4% GROMACS 6% DL_POLY CP2K VASP 5% 8% 19% http://www.nu-fuse.com HECToR usage statistics Phase 3 statistics (Nov 2011 - Apr 2013) 35% of the Chemistry software on HECToR is using DFT methods.
    [Show full text]
  • Electronic and Structural Properties of Silicene and Graphene Layered Structures
    Wright State University CORE Scholar Browse all Theses and Dissertations Theses and Dissertations 2012 Electronic and Structural Properties of Silicene and Graphene Layered Structures Patrick B. Benasutti Wright State University Follow this and additional works at: https://corescholar.libraries.wright.edu/etd_all Part of the Physics Commons Repository Citation Benasutti, Patrick B., "Electronic and Structural Properties of Silicene and Graphene Layered Structures" (2012). Browse all Theses and Dissertations. 1339. https://corescholar.libraries.wright.edu/etd_all/1339 This Thesis is brought to you for free and open access by the Theses and Dissertations at CORE Scholar. It has been accepted for inclusion in Browse all Theses and Dissertations by an authorized administrator of CORE Scholar. For more information, please contact [email protected]. Electronic and Structural Properties of Silicene and Graphene Layered Structures A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science by Patrick B. Benasutti B.S.C.E. and B.S., Case Western Reserve and Wheeling Jesuit, 2009 2012 Wright State University Wright State University SCHOOL OF GRADUATE STUDIES August 31, 2012 I HEREBY RECOMMEND THAT THE THESIS PREPARED UNDER MY SUPERVI- SION BY Patrick B. Benasutti ENTITLED Electronic and Structural Properties of Silicene and Graphene Layered Structures BE ACCEPTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF Master of Science in Physics. Dr. Lok C. Lew Yan Voon Thesis Director Dr. Douglas Petkie Department Chair Committee on Final Examination Dr. Lok C. Lew Yan Voon Dr. Brent Foy Dr. Gregory Kozlowski Dr. Andrew Hsu, Ph.D., MD Dean, Graduate School ABSTRACT Benasutti, Patrick.
    [Show full text]
  • PLUMED User's Guide
    PLUMED User’s Guide A portable plugin for free-energy calculations with molecular dynamics Version 1.3.0 – Nov 2011 Contents 1 Introduction 5 1.1 What is PLUMED?.........................5 1.2 Supported codes . .6 1.3 Features . .7 1.4 New in version 1.3 . .8 1.5 Restrictions . .9 1.6 The PLUMED package . .9 1.7 Online resources . 10 1.8 Credits . 11 1.9 Citing PLUMED ........................... 11 1.10 License . 11 2 Installation 13 2.1 Compiling PLUMED ......................... 13 2.1.1 Compiling the ACEMD plugin with PLUMED ...... 17 2.2 Including reconnaissance metadynamics . 18 2.3 Testing the installation . 19 2.4 Back to the original code . 21 2.5 The Python interface to PLUMED ................. 21 3 Running free-energy simulations 23 3.1 How to activate PLUMED ..................... 23 3.2 The input file . 25 3.3 A note on units . 27 3.4 Metadynamics . 27 3.4.1 Typical output . 27 3.4.2 Bias potential . 28 3.4.3 Well-tempered metadynamics . 29 1 3.4.4 Restarting a metadynamics run . 30 3.4.5 Using GRID ........................ 30 3.4.6 Multiple walkers . 35 3.4.7 Monitoring a collective variable without biasing it . 35 3.4.8 Defining an interval . 36 3.4.9 Inversion condition . 38 3.5 Running in parallel . 40 3.6 Replica exchange methods . 40 3.6.1 Parallel tempering metadynamics . 41 3.6.2 Bias exchange simulations . 42 3.7 Umbrella sampling . 45 3.8 Steered MD . 46 3.8.1 Steerplan . 46 3.9 Adiabatic Bias MD .
    [Show full text]
  • "The Optoelectronic Properties of Nanoparticles
    THE UNIVERSITY OF CHICAGO THE OPTOELECTRONIC PROPERTIES OF NANOPARTICLES FROM FIRST PRINCIPLES CALCULATIONS A DISSERTATION SUBMITTED TO THE FACULTY OF THE INSTITUTE FOR MOLECULAR ENGINEERING IN CANDIDACY FOR THE DEGREE OF DOCTOR OF PHILOSOPHY BY NICHOLAS PETER BRAWAND CHICAGO, ILLINOIS JUNE 2017 Copyright c 2017 by Nicholas Peter Brawand All Rights Reserved To my wife, Tabatha TABLE OF CONTENTS LIST OF FIGURES . vi LIST OF TABLES . xii ABSTRACT . xiv 1 INTRODUCTION . 1 2 FIRST PRINCIPLES CALCULATIONS . 4 2.1 Methodological Foundations . .4 2.1.1 Density Functional Theory . .4 2.1.2 Many Body Perturbation Theory . .6 2.2 The Generalization of Dielectric Dependent Hybrid Functionals to Finite Sys- tems........................................8 2.2.1 Introduction . .8 2.2.2 Theory . 10 2.2.3 Results . 12 2.3 Self-Consistency of the Screened Exchange Constant Functional . 19 2.3.1 Introduction . 19 2.3.2 Results . 20 2.3.3 Effects of Self-Consistency . 30 2.4 Constrained Density Functional Theory and the Self-Consistent Image Charge Method . 35 2.4.1 Introduction . 36 2.4.2 Methodology . 39 2.4.3 Verification . 42 2.4.4 Application to silicon quantum dots . 52 3 OPTICAL PROPERTIES OF ISOLATED NANOPARTICLES . 57 3.1 Nanoparticle Blinking . 57 3.1.1 Introduction . 57 3.1.2 Structural Models . 60 3.1.3 Theory and Computational Methods . 61 3.1.4 Results . 64 3.1.5 Analysis of Blinking Processes . 68 3.2 Tuning the Optical Properties of Nanoparticles: Experiment and Theory . 70 3.2.1 Introduction . 70 3.2.2 Results . 72 4 CHARGE TRANSPORT PROPERTIES OF NANOPARTICLES .
    [Show full text]
  • Introduction to Quantum ESPRESSO for Would-Be Developers
    Introduction to quantum ESPRESSO for would-be developers P. Giannozzi Universit`adi Udine and IOM-Democritos, Trieste, Italy 25 March 2013 { Typeset by FoilTEX { Developers' Resources When everything else fails, read the manual! There isn't a real quantum ESPRESSO manual, but there are some useful resources for developers: • the pseudo-manual Doc/developers-man.pdf, plus other files in Doc/, install/, include/; • scattered documentation in source files; • www.quantum-espresso.org, Resources menu, under Developers' Resources • the q-e project on qe-forge.org, and in particular the bug tracker, forum, mailing list archives, web-svn interface. The quantum ESPRESSO ecosystem (zoo?) Quantum ESPRESSO package porpolio WanT Wannier90 PLUMED SaX YAMBO TDDFPT Xspectra GIPAW GWL NEB PHonon PWCOND Atomic PWscf CP CORE modules 7 (the above slide courtesy of Filippo Spiga) Core quantum ESPRESSO distribution In the following I will refer to the core distribution, including • scripts, installation tools, libraries, common source files; • basic packages { PWscf: self-consistent calculations, structural optimization, molecular dynamics on the ground state; { CP: Car-Parrinello molecular dynamics; { PostProc: data analysis and plotting (requires PWscf) All of the above is stored into a single ('main') SVN repository. Some libraries are however downloaded on demand from the web. Most of what will be said is valid also for additional packages. More quantum ESPRESSO distribution Additional packages that can be downloaded and installed on demand: • atomic: pseudopotential generation • PHonon: Density-Functional Perturbation Theory • NEB: reaction pathways and energy barriers • PWCOND: ballistic conductance • XSPECTRA: calculation of X-Ray spectra • TDDFPT: Time-dependent DFPT (requires PHonon) All these require that the core distribution is already installed.
    [Show full text]