High Temperature High Pressure Properties of Silica From Quantum Monte Carlo
|
|
- Felix Terry
- 5 years ago
- Views:
Transcription
1 High Temperature High Pressure Properties of Silica From Quantum Monte Carlo K.P. Driver, R.E. Cohen, Z. Wu, B. Militzer, P. Lopez Rios, M. Towler, R. Needs, and J.W. Wilkins Funding: NSF, DOE; Computation: NCAR, TeraGrid, NCSA, NERSC, OSC, CCNI Outline: 1) Physics of silica compression 2) Quantum Monte Carlo (QMC) as a benchmarking tool for first principles lattice dynamics. 3) Computing silica thermal EoS using QMC for static lattice and DFT for vibrations. 4) Compare QMC thermodynamic and elastic properties with DFT and experiment.
2 High-Pressure Silica Physics Simplest of Earth's silicates; free phases may exist at base of mantle. Complex series of phase transitions. Quartz to stishovite is a four fold to six fold coordination change. Stishovite undergoes a ferroelastic transition to a CaCl structured phase. 2 CaCl2 structure transforms to a PbO2 structure, which is stable to the core.
3 Motivation for QMC Computation Extreme materials research demands pushing the frontier of first-principle simulations. DFT is current workhorse, but exchange-correlation (XC) functionals leave room for doubt. QMC avoids XC approximation, but is too expensive replace DFT anytime soon. Use QMC as a benchmark tool to spot-check important results. DFT generally great, but can fail: LDA predicts stishovite as ambient ground state. GGA's may improve transition, but may worsen other properties. Goals of this work: Explore feasibility of using QMC with DFPT for high pressure/temperature properties of silica. Thermal equations of state Thermodynamic properties of silica Elastic constants
4 Density Functional Theory and Quantum Monte Carlo DFT WCGGA functional and potentials ABINIT, Opium PPs 9/6/12 atom cells 1)Single particle theory in effective potential. 2)Choose XC functional and pseudopotential. 3)Relax crystal structures. 4)Compute energy and wavefunction. DFT Vol/BulkMod at T=300 K, P= 0 GPa phase quartz stishovite apbo2 V0 QMC 254(2) 159.0(4) 154.8(1) V0 WCGGA V0 Experiment N/A K0 QMC 32(6) 305(20) 329(4) K0 WCGGA K0 Experiment N/A DFT(WCGGA) V too high, K too low QMC generally improves QMC 1)Explicit many body method. 2)Use DFT's relaxed crystal structures. 3)Optimize DFT wavefunction (fixed nodes). 4)Compute energy stochastically. Finite size errors < 30 mev MPC, 72/162/96 atom simulation cells CASINO, B spline basis, time step Jastrow with 2 & 3 body, PW (44 parameters)
5 Computing P V T Equations of State Compute quasi harmonic Helmholtz free energy QMC Static Energy vs Volume F QHA=E static V F vibration V, T Compute static lattice energy with QMC Dominant energy contribution Most accurate method available for solids CASINO code T=0 K Transition Pressure (GPa) Quartz Stishovite CaCl2 apbo2 Experiment QMC WCGGA 6 to 7 ~90 6.3(0.14) 88(8) Quartz Phonon Dispersion (P=0 GPa) Frequency (cm 1) Compute vibrational free energy with DFT Currently too costly for QMC Vibrational energy is small Typically well described in DFT ABINIT, Linear Response, Quasi harmonic Dispersion data from Burkel, et al. Physica B, , pp (1999).
6 Thermal Equations of State P= Quartz F QHA V T Thermal EoS determines fundamental thermodynamic parameters and phase relations. QMC improves agreement with experiment in each phase: quartz, stishovite, PbO2. Stishovite Only small number of measurements for PbO2, making QMC most accurate available. PbO2 Gray shading indicates one standard deviation of statistical error.
7 Enthalpy Difference (Ha/SiO2) Vquartz stish % Volume Difference VCaCl2 PbO2 H PbO2 CaCl2 (Ha/SiO2) Hstish quartz(ha/sio2) Enthalpy Difference and Volume Difference Pressure (GPa) Pressure (GPa) Enthalpy differences and errors determine equilibrium phase relations. CaCl PbO enthalpy change is not measurable; phases are not quenchable to zero pressure. 2 2 CaCl2 & PbO2 enthalpy curves lie very close together compared to quartz stishovite. Volume change in quartz stishovite is 20 times larger than in CaCl2 PbO2.
8 Silica Phase Boundaries Equilibrium phase boundaries computed from Gibbs free energies. Temperature (K) Metastable quartz stishovite transition measured with thermocal or shock. QMC agrees well. Temperature (K) QMC CaCl2 PbO2 transition lies between two measurements. DFT(WCGGA) boundary 4 GPa too low for quartz stishovite and within statistical error of QMC for CaCl2 PbO2 Pressure (GPa)
9 c11 c12 (GPa) Stishovite Shear Modulus c ij = 1 2 E V 2 Pressure (GPa) Stishovite to CaCl2 transition is driven by instability in the elastic shear modulus. VMC modulus computed at several pressures and DMC checks at endpoints. Shear modulus computed from strain energy relation (brute force 1000 CPU cost of DFT) All methods roughly agree, with the shear modulus going unstable around 50 GPa.
10 Bulk Modulus K = V P V T QMC and DFT(WCGGA) temperature and pressure dependence of K. Bulk Modulus (GPa) K decreases linearly with T K increases linearly with P QMC and DFT generally agree. (T dependence comes from DFT) Temperature (K)
11 Thermal Expansivity Thermal Expansivity (10 5K 1) = 1 V V T P QMC and DFT(WCGGA) temperature and pressure dependence of. QMC shows best agreement for quartz. Experimentally, quartz appears at 846 K we only consider quartz. QMC and DFT show good agreement with stishovite measurements. Temperature (K) PbO2 curves are the best available.
12 Heat Capacity C p= H T P QMC and DFT(WCGGA) temperature and pressure dependence of Cp. Cp/R QMC and DFT results are nearly identical. Good agreement with experiment for quartz and stishovite. Temperature (K) QMC PbO2 curves are best available.
13 Grüneisen Ratio =V P U V = KT V CV change in P due to increase in energy density Grüneisen ratio quantifies relation between thermal and elastic properties. Close agreement with experiment for quartz. QMC stishovite and PbO2 curves are the best available. Pressure (GPa)
14 Conclusions QMC with DFT phonons is a route to benchmark solid ab initio high P/T calculations. QMC gives excellent agreement with experiment for thermal equations of state and thermodynamic properties. QMC gives accurate phase transformation boundaries for quartz-stishovite and CaCl2 PbO2 transitions. QMC computation of shear modulus softening in stishovite confirms the stishovite CaCl2 transition is near 50 GPa at low temperatures. Publication is written, few details remain to be worked out.
15 Backup Slides 1)DMC time step convergence 2)DMC finite size convergence 3)Shear modulus strain energy technique 4)DMC and VMC strain energy curves 5)Statistical error propagation 6)Disordered Silica
16 DMC Time Step Convergence Time step of a.u. is converged to within 30 mev.
17 DMC Finite Size Convergence
18 Calculate elastic constants by straining the lattice Stishovite DFT(WC) Strain energy density relation: Strain the lattice: R ' =[ I ] R For a volume conserving strain: E 1 = c ijkl ij kl V processor hours c ijkl = 1 2 E V 2 Elastic constants obtained from curvature of energy strain curve Double well at 280 Bohr3 indicates elastic instability of stishovite CaCl becoming more stable that stishovite under pressure 2 Feasibility of Elastic Constants in QMC Elasticity is a tough problem for QMC: energy differences ~ ev Extremely expensive to get accurate error bars for large (100 atom) systems Through parallel computation on large supercomputers, it's possible to succeed.
19 QMC Energy vs. Strain Curves: The Brute Force Method 1) Take optimized input structures from DFT(WC) (we can't do forces in QMC yet) 2) Run QMC on thousands of processors for a few days until error bars are sufficiently small for each structure. VMC (500,000 hrs) DMC (additional 1.3 million hrs) QMC at this accuracy level is 1200 times more expensive than DFT. QMC error bars must be made much smaller than the strain energy differences. VMC error bars decrease twice as fast as DMC error bars. Highest pressure curves are most difficult to fit and require smallest error bars (work on high pressure curves is still in progress).
20 Statistical Error Propagation 1) Propagation of error equation from Taylor expanding a function about the mean values of its parmeters: F = i, j F F Ci C j Cov [i, j ] P = i 2) Propagation with Monte Carlo Pi Ei E i actual data set Allow random Gaussian fluctuations on QMC energies Ei with stdv Ei Fit Vinet Equation to new set of energies and compute property Allow random Gaussian fluctuations on QMC energies Ei with stdv Ei Allow random Gaussian fluctuations on QMC energies Ei with stdv Ei Fit Vinet Equation to new set of energies and compute property Fit Vinet Equation to new set of energies and compute property standard deviation of synthetic data sets gives uncertainty in property......
21 Does silica undergo order disorder transition at high T? Stishovite Stishovite and post stishovite phases can be described as distorted hcp arrays of oxygen ions with octahedral intersticies half filled by silicon ions. Suggests AB alloy in which silicons can be disordered. Goal: get a rough estimate of disorder temperature without doing any significant additional calculations. Approach: 1) Stishovite is unkinked, PbO2 is maximally kinked. Guess a transition structure and compute it's energy in DFT. PbO2 2)Assume we can estimate the enthalpy of a disordered silica phase by weighting three ordered phases (stish, PbO2, and the transition structure) according to their Boltzmann probabilities and including a configurational entropy term. 3) Compute entropy of disordered phase with lattice gas model?
Exploring deep Earth minerals with accurate theory
Exploring deep Earth minerals with accurate theory K.P. Driver, R.E. Cohen, Z. Wu, B. Militzer, P. Lopez Rios, M. Towler, R. Needs, and J.W. Wilkins Funding: NSF, DOE; Computation: NCAR, TeraGrid, NCSA,
More informationEstablishing Quantum Monte Carlo and Hybrid Density Functional Theory as Benchmarking Tools for Complex Solids
Establishing Quantum Monte Carlo and Hybrid Density Functional Theory as Benchmarking Tools for Complex Solids Kevin P. Driver, Ph.D. Defense, February 1, 011 Outline 1) Introduction: Quantum mechanics
More informationQuantum Monte Carlo Benchmarks Density Functionals: Si Defects
Quantum Monte Carlo Benchmarks Density Functionals: Si Defects K P Driver, W D Parker, R G Hennig, J W Wilkins (OSU) C J Umrigar (Cornell), R Martin, E Batista, B Uberuaga (LANL), J Heyd, G Scuseria (Rice)
More informationKevin Driver 1 Shuai Zhang 1 Burkhard Militzer 1 R. E. Cohen 2.
Quantum Monte Carlo Simulations of a Single Iron Impurity in MgO Kevin Driver 1 Shuai Zhang 1 Burkhard Militzer 1 R. E. Cohen 2 1 Department of Earth & Planetary Science University of California, Berkeley
More informationPhonon wavefunctions and electron phonon interactions in semiconductors
Phonon wavefunctions and electron phonon interactions in semiconductors Bartomeu Monserrat bm418@cam.ac.uk University of Cambridge Quantum Monte Carlo in the Apuan Alps VII QMC in the Apuan Alps VII Bartomeu
More informationAll-Electron Path Integral Monte Carlo (PIMC) Simulations of Warm Dense Matter: Application to Water and Carbon Plasmas
All-Electron Path Integral Monte Carlo (PIMC) Simulations of Warm Dense Matter: Application to Water and Carbon Plasmas Kevin Driver and Burkhard Militzer Department of Earth and Planetary Science University
More informationQuantum Monte Carlo backflow calculations of benzene dimers
Quantum Monte Carlo backflow calculations of benzene dimers Kathleen Schwarz*, Cyrus Umrigar**, Richard Hennig*** *Cornell University Department of Chemistry, **Cornell University Department of Physics,
More informationQMC dissociation energy of the water dimer: Time step errors and backflow calculations
QMC dissociation energy of the water dimer: Time step errors and backflow calculations Idoia G. de Gurtubay and Richard J. Needs TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay.
More informationQMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters
QMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters Idoia G. de Gurtubay TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay. Quantum
More informationThe QMC Petascale Project
The QMC Petascale Project Richard G. Hennig What will a petascale computer look like? What are the limitations of current QMC algorithms for petascale computers? How can Quantum Monte Carlo algorithms
More informationThermodynamics of Solids: Harmonic and Quasi-harmonic Approximations
Thermodynamics of Solids: Harmonic and Quasi-harmonic Approximations, USA, July 9-14, 2017 Alessandro Erba Dipartimento di Chimica, Università di Torino (Italy) alessandro.erba@unito.it 2017 Outline -
More informationBinding energy of 2D materials using Quantum Monte Carlo
Quantum Monte Carlo in the Apuan Alps IX International Workshop, 26th July to 2nd August 2014 The Apuan Alps Centre for Physics @ TTI, Vallico Sotto, Tuscany, Italy Binding energy of 2D materials using
More informationOn Dynamic and Elastic Stability of Lanthanum Carbide
Journal of Physics: Conference Series On Dynamic and Elastic Stability of Lanthanum Carbide To cite this article: B D Sahoo et al 212 J. Phys.: Conf. Ser. 377 1287 Recent citations - Theoretical prediction
More informationSpin crossovers in the Earth mantle. Spin crossovers in the Earth mantle
Spin crossovers in the Earth mantle Spin crossovers in the Earth mantle Renata M. Wentzcovitch Dept. of Chemical Engineering and Materials Science Minnesota Supercomputing Institute Collaborators Han Hsu
More informationStructural Calculations phase stability, surfaces, interfaces etc
Structural Calculations phase stability, surfaces, interfaces etc Keith Refson STFC Rutherford Appleton Laboratory September 19, 2007 Phase Equilibrium 2 Energy-Volume curves..................................................................
More informationFinite-temperature magnetism in bcc Fe under compression. Xianwei Sha* and R. E. Cohen
Finite-temperature magnetism in bcc Fe under compression Xianwei Sha* and R. E. Cohen Carnegie Institution of Washington, 5251 Broad Branch Road, NW, Washington, D. C. 20015, U. S. A. We investigate the
More informationPseudopotentials: design, testing, typical errors
Pseudopotentials: design, testing, typical errors Kevin F. Garrity Part 1 National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015 Parameter free calculations.
More informationComplementary approaches to high T- high p crystal structure stability and melting!
Complementary approaches to high T- high p crystal structure stability and melting! Dario ALFÈ Department of Earth Sciences & Department of Physics and Astronomy, Thomas Young Centre@UCL & London Centre
More informationQuantum Monte Carlo methods
Quantum Monte Carlo methods Lubos Mitas North Carolina State University Urbana, August 2006 Lubos_Mitas@ncsu.edu H= 1 2 i i 2 i, I Z I r ii i j 1 r ij E ion ion H r 1, r 2,... =E r 1, r 2,... - ground
More informationRustam Z. Khaliullin University of Zürich
Rustam Z. Khaliullin University of Zürich Molecular dynamics (MD) MD is a computational method for simulating time evolution of a collection of interacting atoms by numerically integrating Newton s equation
More informationMustafa Uludogan 1, Tahir Cagin, William A. Goddard, III Materials and Process Simulation Center, Caltech, Pasadena, CA 91125, U.S.A.
Ab Initio Studies On Phase Behavior of Barium Titanate Mustafa Uludogan 1, Tahir Cagin, William A. Goddard, III Materials and Process Simulation Center, Caltech, Pasadena, CA 91125, U.S.A. 1 Physics Department,
More informationRoger Johnson Structure and Dynamics: Displacive phase transition Lecture 9
9.1. Summary In this Lecture we will consider structural phase transitions characterised by atomic displacements, which result in a low temperature structure that is distorted compared to a higher temperature,
More informationElectronic structure of solid FeO at high pressures by quantum Monte Carlo methods
Physics Procedia 3 (2010) 1437 1441 www.elsevier.com/locate/procedia Electronic structure of solid FeO at high pressures by quantum Monte Carlo methods Jindřich Kolorenč a and Lubos Mitas a a Department
More informationComputational support for a pyrolitic lower mantle containing ferric iron
SUPPLEMENTARY INFORMATION DOI: 1.138/NGEO2458 Computational support for a pyrolitic lower mantle containing ferric iron Xianlong Wang, Taku Tsuchiya and Atsushi Hase NATURE GEOSCIENCE www.nature.com/naturegeoscience
More informationQuantum Monte Carlo simulations on silicate perovskite and other high pressure phases
Logo Quantum Monte Carlo simulations on silicate perovskite and other high pressure phases Carnegie Logo 2.1 RE 1,2 Cohen 1Geophysical Laboratory, Carnegie Institution 2Dept. Earth Science, University
More informationTHERMOPHYSICAL PROPERTIES OF THORIUM COMPOUNDS FROM FIRST PRINCIPLES
THERMOPHYSICAL PROPERTIES OF THORIUM COMPOUNDS FROM FIRST PRINCIPLES Vinayak Mishra a,* and Shashank Chaturvedi a a Computational Analysis Division, Bhabha Atomic Research Centre, Visakhapatnam 530012,
More informationFirst-principles calculations of bulk and interfacial thermodynamic properties for fcc-based Al-Sc alloys
PHYSICAL REVIEW B VOLUME 57, NUMBER 18 1 MAY 1998-II First-principles calculations of bulk and interfacial thermodynamic properties for fcc-based Al-Sc alloys Mark Asta, S. M. Foiles, and A. A. Quong*
More informationPhonon calculations with SCAN
Workshop on the SCAN density functional: Fundamentals, practices, and extensions Temple university, Philadelphia May 18th, 2017 Hands-on tutorial 3 Phonon calculations with SCAN Yubo Zhang and Jianwei
More informationKevin P. Driver University of California, Berkeley 307 McCone Hall Berkeley, CA Phone: (614)
Curriculum Vitae Kevin P. Driver University of California, Berkeley 307 McCone Hall Berkeley, CA 94720-4767 Phone: (614)214-8292 kevin.driver@gmail.com Career Summary Ph.D. in Physics with a research background
More informationStructural and thermal properties of Fe 2 (Zr,Nb) system in C15, C14 and C36 Laves phases: First-Principles study
Structural and thermal properties of Fe 2 (Zr,Nb) system in, and Laves phases: First-Principles study L. RABAHI 1, D. BRADAI 2 and A. KELLOU 3 1 Centre National de Recherche en Soudage et Contrôle, Route
More informationHigh-throughput Simulations of Shock-waves Using MD. Jacob Wilkins University of York
High-throughput Simulations of Using MD Jacob Wilkins University of York jacob.wilkins@york.ac.uk Mar 2018 Outline 1 2 3 4 5 What are shock-waves? are pressure-waves moving through material faster than
More informationarxiv: v1 [cond-mat.mtrl-sci] 21 Dec 2007
Quantum Monte Carlo calculations of structural properties of FeO under pressure Jindřich Kolorenč and Lubos Mitas Department of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina
More informationalloys at high pressures. The central peaks at 0 mev correspond to the elastic scattering.
1 2 3 4 5 Supplementary Figure 1 Representative NRIXS spectra of the basaltic glass and Fe-rich alloys at high pressures. The central peaks at 0 mev correspond to the elastic scattering. Open symbols:
More informationQuantum Monte Carlo Calculations of Point Defects in Alumina
Quantum Monte Carlo Calculations of Nicholas Hine Condensed Matter Theory Group Imperial College London Sunday 22nd July 2007 Dramatis Personae Collaborators Matthew Foulkes Mike Finnis Kilian Frensch
More informationThe Nature of the Interlayer Interaction in Bulk. and Few-Layer Phosphorus
Supporting Information for: The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus L. Shulenburger, A.D. Baczewski, Z. Zhu, J. Guan, and D. Tománek, Sandia National Laboratories, Albuquerque,
More informationarxiv:cond-mat/ v1 10 Jun 1994 K. M. Rabe
October 2, 2018 Phase transitions in BaTiO 3 from first principles W. Zhong and David Vanderbilt Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08855-0849 arxiv:cond-mat/9406049v1
More informationINTRODUCTION TO THE PHYSICS OF THE EARTH S INTERIOR
INTRODUCTION TO THE PHYSICS OF THE EARTH S INTERIOR SECOND EDITION JEAN-PAULPOIRIER Institut de Physique du Globe de Paris PUBLISHED BY THE PRESS SYNDICATE OF THE UNIVERSITY OF CAMBRIDGE The Pitt Building,
More informationFirst- principles studies of spin-crossover molecules
Vallico di Sotto, 30 th July 2012 First- principles studies of spin-crossover molecules Andrea Droghetti and Stefano Sanvito School of Physics and CRANN, Trinity College Dublin Dario Alfe' London Centre
More informationSupporting Online Material (1)
Supporting Online Material The density functional theory (DFT) calculations were carried out using the dacapo code (http://www.fysik.dtu.dk/campos), and the RPBE (1) generalized gradient correction (GGA)
More informationSMARTMET project: Towards breaking the inverse ductility-strength relation
SMARTMET project: Towards breaking the inverse ductility-strength relation B. Grabowski, C. Tasan SMARTMET ERC advanced grant 3.8 Mio Euro for 5 years (Raabe/Neugebauer) Adaptive Structural Materials group
More informationPressure Volume Temperature Equation of State
Pressure Volume Temperature Equation of State S.-H. Dan Shim ( ) Acknowledgement: NSF-CSEDI, NSF-FESD, NSF-EAR, NASA-NExSS, Keck Equations relating state variables (pressure, temperature, volume, or energy).
More informationA Thermal Gradient Approach for the Quasi-Harmonic Approximation and its Application to Improved Treatment of Anisotropic Expansion
A Thermal Gradient Approach for the Quasi-Harmonic Approximation and its Application to Improved Treatment of Anisotropic Expansion arxiv:1712.00936v3 [cond-mat.mtrl-sci] 7 Oct 2018 Nathan S. Abraham and
More informationAb Initio Study of the Mechanical, Electronic, Thermal and Optical Properties of Ge 2 Sb 2 Te 5
Ab Initio Study of the Mechanical, Electronic, Thermal and Optical Properties of Ge 2 Sb 2 Te 5 Odhiambo H. 1, Amolo G. 2, Makau N. 2, Othieno H. 1, and Oduor A. 1 1 Department of Physics, Maseno University,
More informationMaterials Chemistry A
Journal of Materials Chemistry A PAPER Cite this: J. Mater. Chem. A, 2015,3, 8002 Received 22nd December 2014 Accepted 4th March 2015 DOI: 10.1039/c4ta07062c www.rsc.org/materialsa View Journal View Issue
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationPredicting the Dislocation Nucleation Rate as a Function of Temperature and Stress. Abstract
Journal of Materials Research, Invited Feature Article, in press (2011) Predicting the Dislocation Nucleation Rate as a Function of Temperature and Stress Seunghwa Ryu 1, Keonwook Kang 2, and Wei Cai 3
More informationStructure, energetics, and mechanical stability of Fe-Cu bcc alloys from first-principles calculations
Structure, energetics, and mechanical stability of Fe-Cu bcc alloys from first-principles calculations Jefferson Z. Liu, A. van de Walle, G. Ghosh, and M. Asta* Department of Materials Science and Engineering,
More informationKinetic lattice Monte Carlo simulations of diffusion processes in Si and SiGe alloys
Kinetic lattice Monte Carlo simulations of diffusion processes in Si and SiGe alloys, Scott Dunham Department of Electrical Engineering Multiscale Modeling Hierarchy Configuration energies and transition
More informationMott metal-insulator transition on compressible lattices
Mott metal-insulator transition on compressible lattices Markus Garst Universität zu Köln T p in collaboration with : Mario Zacharias (Köln) Lorenz Bartosch (Frankfurt) T c Mott insulator p c T metal pressure
More informationQMC in the Apuan Alps IV Vallico Sotto, 28th July 2008
Van der Waals forces in graphitic nanostructures Norbert Nemec QMC in the Apuan Alps IV Vallico Sotto, 28th July 2008 1 History Johannes Diderik van der Waals (1837-1923): Postulation of intermolecular
More informationLecture 11 - Phonons II - Thermal Prop. Continued
Phonons II - hermal Properties - Continued (Kittel Ch. 5) Low High Outline Anharmonicity Crucial for hermal expansion other changes with pressure temperature Gruneisen Constant hermal Heat ransport Phonon
More informationKey concepts in Density Functional Theory (II)
Kohn-Sham scheme and band structures European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Present Address: LPMCN Université
More informationQuantum Monte Carlo calculations of the structural properties and the B1-B2 phase transition of MgO
Quantum Monte Carlo calculations of the structural properties and the B1-B2 phase transition of MgO D. Alfè, 1,2 M. Alfredsson, 1 J. Brodholt, 1 M. J. Gillan, 2 M. D. Towler, 3 and R. J. Needs 3 1 Department
More informationCrossing the barriers - simulations of activated processes
Crossing the barriers - simulations of activated processes Mgr. Ján Hreha for 6 th Student Colloquium and School on Mathematical Physics Faculty of Mathematics, Physics and Informatics Comenius University
More informationContents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21
Preface Reference tables Table A Counting and combinatorics formulae Table B Useful integrals, expansions, and approximations Table C Extensive thermodynamic potentials Table D Intensive per-particle thermodynamic
More informationPhonon Dispersion, Interatomic Force Constants Thermodynamic Quantities
Phonon Dispersion, Interatomic Force Constants Thermodynamic Quantities Umesh V. Waghmare Theoretical Sciences Unit J N C A S R Bangalore ICMR OUTLINE Vibrations and interatomic force constants (IFC) Extended
More informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
More informationFrom Atoms to Materials: Predictive Theory and Simulations
From Atoms to Materials: Predictive Theory and Simulations Week 3 Lecture 4 Potentials for metals and semiconductors Ale Strachan strachan@purdue.edu School of Materials Engineering & Birck anotechnology
More information1.050 Engineering Mechanics. Lecture 22: Isotropic elasticity
1.050 Engineering Mechanics Lecture 22: Isotropic elasticity 1.050 Content overview I. Dimensional analysis 1. On monsters, mice and mushrooms 2. Similarity relations: Important engineering tools II. Stresses
More informationAb initio calculations of the cohesive energy and the bulk modulus of Aluminium
Ab initio calculations of the cohesive energy and the bulk modulus of Aluminium R. Gaudoin and W. M. C. Foulkes CMTH Group, Blackett Laboratory, Imperial College of Science, Technology and Medicine, Prince
More informationUniversity of Eldoret May 13, 2013.
Materials for the hydrogen energy economy: Magnesium hydride (MgH 2 ) and Lithium hydride (LiH) Magero Denis Makau Nicholas Amolo George Lusweti Kituyi University of Eldoret May 13, 2013. Introduction
More informationVibrational modes of silicon metalattices from atomistic and finite-element calculations
Vibrational modes of silicon metalattices from atomistic and finite-element calculations Yihuang Xiong, Weinan hen, Ismaila Dabo Materials Science and Engineering, Penn State University October 6 th, 216
More informationThe high-pressure phase transitions of silicon and gallium nitride: a comparative study of Hartree Fock and density functional calculations
J. Phys.: Condens. Matter 8 (1996) 3993 4000. Printed in the UK The high-pressure phase transitions of silicon and gallium nitride: a comparative study of Hartree Fock and density functional calculations
More informationthe abdus salam international centre for theoretical physics
united nations educational, scientific and cultural organization the abdus salam international centre for theoretical physics international atomic energy agency SMR 1595-16 Joint DEMOCRITOS - ICTP School
More informationSupporting Information. Potential semiconducting and superconducting metastable Si 3 C. structures under pressure
Supporting Information Potential semiconducting and superconducting metastable Si 3 C structures under pressure Guoying Gao 1,3,* Xiaowei Liang, 1 Neil W. Ashcroft 2 and Roald Hoffmann 3,* 1 State Key
More informationFixed-Node quantum Monte Carlo for Chemistry
Fixed-Node quantum Monte Carlo for Chemistry Michel Caffarel Lab. Physique et Chimie Quantiques, CNRS-IRSAMC, Université de Toulouse e-mail : caffarel@irsamc.ups-tlse.fr. p.1/29 The N-body problem of Chemistry
More informationEstablishing Quantum Monte Carlo and Hybrid Density Functional Theory as benchmarking tools for complex solids
Establishing Quantum Monte Carlo and Hybrid Density Functional Theory as benchmarking tools for complex solids DISSERTATION Presented in Partial Fulfillment of the Requirements for the Degree Doctor of
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term 2013
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.044 Statistical Physics I Spring Term 2013 Problem 1: The Big Bang Problem Set #9 Due in hand-in box by 4;00 PM, Friday, April 19 Early in the
More informationComparison of various abinitio codes used in periodic calculations
Comparison of various abinitio codes used in periodic calculations 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India & Center for Materials Science and Nanotechnology,
More informationNeutron Matter: EOS, Spin and Density Response
Neutron Matter: EOS, Spin and Density Response LANL : A. Gezerlis, M. Dupuis, S. Reddy, J. Carlson ANL: S. Pieper, R.B. Wiringa How can microscopic theories constrain mean-field theories and properties
More informationPractical calculations using first-principles QM Convergence, convergence, convergence
Practical calculations using first-principles QM Convergence, convergence, convergence Keith Refson STFC Rutherford Appleton Laboratory September 18, 2007 Results of First-Principles Simulations..........................................................
More information* Theoretische Physik II, Universitat Dortmund, Dortmund, Germany
JOURNAL DE PHYSIQUE IV Colloque C8, suppl6ment au Journal de Physique III, Volume 5, dkembre 1995 Structural Phase Transformation and Phonon Softening in Iron-Based Alloys H.C. Herper, E. Hoffmann, P.
More informationSupplementary Information for. Universal elastic-hardening-driven mechanical instability in α-quartz and quartz. homeotypes under pressure
Supplementary Information for Universal elastic-hardening-driven mechanical instability in α-quartz and quartz homeotypes under pressure Juncai Dong, Hailiang Zhu, and Dongliang Chen * Beijing Synchrotron
More informationGeneration of Thermal Scattering Laws for YH 2 using Ab Initio Methods
Generation of Thermal Scattering Laws for YH 2 using Ab Initio Methods Michael L. Zerkle Bettis Atomic Power Laboratory WPEC SG42 Meeting May 18, 2015 May 18-19, 2015 WPEC SG42 Slide 1 Outline Motivation
More informationStrain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity
Strain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity Torino, Italy, September 4-9, 2016 Alessandro Erba Dipartimento di Chimica, Università di Torino (Italy) alessandro.erba@unito.it
More informationAb-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University
Ab-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University Hands-on Tutorial Workshop on Ab-Initio Molecular Simulations, Fritz- Haber-Institut, Berlin, July 12-21,
More informationCrystal Relaxation, Elasticity, and Lattice Dynamics
http://exciting-code.org Crystal Relaxation, Elasticity, and Lattice Dynamics Pasquale Pavone Humboldt-Universität zu Berlin http://exciting-code.org PART I: Structure Optimization Pasquale Pavone Humboldt-Universität
More informationarxiv: v1 [cond-mat.mtrl-sci] 18 Sep 2013
Accuracy of generalized gradient approximation functionals for density functional perturbation theory calculations arxiv:1309.4805v1 [cond-mat.mtrl-sci] 18 Sep 2013 Lianhua He, 1 Fang Liu, 2 Geoffroy Hautier,
More information(# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble
Recall from before: Internal energy (or Entropy): &, *, - (# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble & = /01Ω maximized Ω: fundamental statistical quantity
More informationAb-initio molecular dynamics for High pressure Hydrogen
Ab-initio molecular dynamics for High pressure Hydrogen Claudio Attaccalite Institut d'electronique, Microélectronique et Nanotechnologie (IEMN), Lille Outline A brief introduction to Quantum Monte Carlo
More informationT. Egami. Model System of Dense Random Packing (DRP)
Introduction to Metallic Glasses: How they are different/similar to other glasses T. Egami Model System of Dense Random Packing (DRP) Hard Sphere vs. Soft Sphere Glass transition Universal behavior History:
More informationequation of state and bulk modulus of AlP, AlAs and AlSb semiconductor compounds
PRAMANA c Indian Academy of Sciences Vol. 64, No. 1 journal of January 2005 physics pp. 153 158 Total energy, equation of state and bulk modulus of AlP, AlAs and AlSb semiconductors A R JIVANI, H J TRIVEDI,
More informationThermodynamics and Phase Transitions in Minerals
Studiengang Geowissenschaften M.Sc. Wintersemester 2004/05 Thermodynamics and Phase Transitions in Minerals Victor Vinograd & Andrew Putnis Basic thermodynamic concepts One of the central themes in Mineralogy
More informationab initio Lattice Vibrations: Calculating the Thermal Expansion Coeffcient Felix Hanke & Martin Fuchs June 30, 2009 This afternoon s plan
ab initio Lattice Vibrations: Calculating the Thermal Expansion Coeffcient Felix Hanke & Martin Fuchs June 3, 29 This afternoon s plan introductory talk Phonons: harmonic vibrations for solids Phonons:
More informationAb initio Rutile-Cristobalite Transitions in Silicon Dioxide and Titanium Dioxide
20 Ab initio Rutile-Cristobalite Transitions in Silicon Dioxide and Titanium Dioxide Moon, Timothy Y. ; Kroll, Peter Department of Chemistry and Biochemistry, The University of Texas at Arlington, Arlington,
More informationOutline. Introduction: graphene. Adsorption on graphene: - Chemisorption - Physisorption. Summary
Outline Introduction: graphene Adsorption on graphene: - Chemisorption - Physisorption Summary 1 Electronic band structure: Electronic properties K Γ M v F = 10 6 ms -1 = c/300 massless Dirac particles!
More informationDFT EXERCISES. FELIPE CERVANTES SODI January 2006
DFT EXERCISES FELIPE CERVANTES SODI January 2006 http://www.csanyi.net/wiki/space/dftexercises Dr. Gábor Csányi 1 Hydrogen atom Place a single H atom in the middle of a largish unit cell (start with a
More informationReconstruction and intermixing in thin Ge layers on Si 001
Reconstruction and intermixing in thin Ge layers on Si 001 L. Nurminen, 1 F. Tavazza, 2 D. P. Landau, 1,2 A. Kuronen, 1 and K. Kaski 1 1 Laboratory of Computational Engineering, Helsinki University of
More informationInternational Physics Course Entrance Examination Questions
International Physics Course Entrance Examination Questions (May 2010) Please answer the four questions from Problem 1 to Problem 4. You can use as many answer sheets you need. Your name, question numbers
More informationFIRST PRINCIPLES STUDY OF AlBi
Available at: http://publications.ictp.it IC/2008/025 United Nations Educational, Scientific and Cultural Organization and International Atomic Energy Agency THE ABDUS SALAM INTERNATIONAL CENTRE FOR THEORETICAL
More informationWave Function Optimization and VMC
Wave Function Optimization and VMC Jeremy McMinis The University of Illinois July 26, 2012 Outline Motivation History of Wave Function Optimization Optimization in QMCPACK Multideterminant Example Motivation:
More informationLangevin Dynamics in Constant Pressure Extended Systems
Langevin Dynamics in Constant Pressure Extended Systems D. Quigley and M.I.J. Probert CMMP 2004 1 Talk Outline Molecular dynamics and ensembles. Existing methods for sampling at NPT. Langevin dynamics
More informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationSecond harmonic generation in silicon waveguides strained by silicon nitride
Second harmonic generation in silicon waveguides strained by silicon nitride M. Cazzanelli, 1 F. Bianco, 1 E. Borga, 1 G. Pucker, 2 M. Ghulinyan, 2 E. Degoli, 3 E. Luppi, 4,7 V. Véniard, 4 S. Ossicini,
More informationNCTS One-Day Workshop. Yasutami Takada. cf. YT, Ryo Maezono, & Kanako Yoshizawa, Phys. Rev. B 92, (2015)
NCTS One-Day Workshop Emergence of a spin-singlet resonance state with Kondo temperature well beyond 1000K in the proton-embedded electron gas: New route to room-temperature superconductivity Yasutami
More informationTime-dependent linear-response variational Monte Carlo.
Time-dependent linear-response variational Monte Carlo. Bastien Mussard bastien.mussard@colorado.edu https://mussard.github.io/ Julien Toulouse julien.toulouse@upmc.fr Sorbonne University, Paris (web)
More informationPhysics of disordered materials. Gunnar A. Niklasson Solid State Physics Department of Engineering Sciences Uppsala University
Physics of disordered materials Gunnar A. Niklasson Solid State Physics Department of Engineering Sciences Uppsala University Course plan Familiarity with the basic description of disordered structures
More informationUNIVERSITY COLLEGE LONDON. University of London EXAMINATION FOR INTERNAL STUDENTS. For The Following Qualifications:-
UNIVERSITY COLLEGE LONDON University of London EXAMINATION FOR INTERNAL STUDENTS For The Following Qualifications:- B.Sc. M.Sci. Physics 1B28: Thermal Physics COURSE CODE : PHYSIB28 UNIT VALUE : 0.50 DATE
More informationPractical Guide to Density Functional Theory (DFT)
Practical Guide to Density Functional Theory (DFT) Brad Malone, Sadas Shankar Quick recap of where we left off last time BD Malone, S Shankar Therefore there is a direct one-to-one correspondence between
More informationMonte Carlo simulations of alloy segregation in PtAg octahedral nanoparticles
Monte Carlo simulations of alloy segregation in PtAg octahedral nanoparticles Louis C. Jones 6/8/12 Abstract Simulations were carried out to investigate phase segregation of insoluble alloy nanoparticles
More information