Path Integral Monte Carlo Simulations on the Blue Waters System. Burkhard Militzer University of California, Berkeley
|
|
- Oliver Nickolas Willis
- 6 years ago
- Views:
Transcription
1 Path Integral Monte Carlo Simulations on the Blue Waters System Burkhard Militzer University of California, Berkeley
2 Outline 1. Path integral Monte Carlo simulation method 2. Application of CH plastic ablator materials 3. Application to warm, dense silicon Kevin Driver, Francois Soubiran, Shuai Zhang Lorin Benedict (LLNL) and Suxing Hu (LLE) Supported by NSF-DOE partnership for plasma science. PRAC computer time on Bluewaters machine at Urbana-Champaign, IL.
3 Phase Diagram of Hydrogen Solar GP: Jupiter, Saturn n=10 18 cm -3 n=10 26 cm -3 Density
4 Juno Mission Now in orbit around Jupiter Interior model to match gravity data Equation of state calculations for hydrogen-helium mixtures Thermodynamics of heavier elements (Z) Mission Timeline: Launch - August 2011 Earth flyby gravity assist - October 2013 Jupiter arrived in July 2016 End of mission much past October 2017
5 1) Path integral Monte Carlo for T>5000K n=10 18 cm -3 n=10 26 cm -3
6 1) Path integral Monte Carlo for T>5000K 2) Density functional molecular dynamics below Born-Oppenheimer approx. MD with classical nuclei: F = m a Forces derived DFT with electrons in the instantaneous ground state.
7 I. Path Integral Monte Carlo
8 Canonical Ensembles: Classical Quantum Boltzmann factor e E / k BT Density matrix ˆρ = e β Ĥ ρ(r, R #,β) = R e β H ˆ R' ρ(r, R #,β) = Thermodynamic averages: S e β E S Ψ S * (R) Ψ S (R') Z Cl = e β E S S Z Q = Tr[ ˆρ] = Ô = Tr[Ô ˆρ] Tr[ ˆρ] dr R e β Ĥ R
9 Step 1 towards the path integral Matrix squaring property of the density matrix Matrix squaring in operator notation: ˆ ρ = e β H ˆ = ( (β / 2) H ˆ e ) ( / 2) H ˆ ) e (β, β = 1 k B T Matrix squaring in real-space notation: R ˆ ρ R $ = dr 1 (β / 2) H ˆ (β / 2) H ˆ R e R 1 R 1 e R $ Matrix squaring in matrix notation: " R' % $ ' $ R ' # $ &' = " R 1 % $ ' $ R ' # $ &' * " R ( % $ ' $ R 1 ' # $ &'
10 Repeat the matrix squaring step Matrix squaring in operator notation: ˆ ρ = e β H ˆ = ( (β / 4) H ˆ e ) 4, β = 1 k B T Matrix squaring in real-space notation: R ˆ ρ R $ = dr 1 dr 2 dr 3 (β / 4) H ˆ R e R 1 R 1 e (β / 4 ) H ˆ (β / 4) H ˆ R 2 R 2 e R 3 R 3 e (β / 4 ) H ˆ R $
11 Path Integrals in Imaginary Time Every particle is represented by a path, a ring polymer. Density matrix: ˆ ρ =e β H ˆ & = e τh ˆ ) M ( + ( + ˆ ρ = e β H ˆ = ( ( ', β=1/k H 1 T,τ =β/m B O ˆ = Tr[ O ˆ ˆ ρ ] k B T, τ = β M Tr[ ˆ ρ ] ( τ ˆ e ) M, β = + * + Trotter break-up: R ˆ ρ R $ = R (e τ H ˆ ) M R $ = dr 1... dr M 1 R e τ ˆ H R 1 R 1 e τ ˆ H R 2... R M 1 e τ ˆ H $ R
12 Path Integrals in Imaginary Time Simplest form for the paths action: primitive approx. Density matrix: ˆ ρ =e β H ˆ & = e τh ˆ ) M ( + ( + ˆ ρ = e β H ˆ = ( ( ', β=1/k H 1 T,τ =β/m B O ˆ = Tr[ O ˆ ˆ ρ ] k B T, τ = β M Tr[ ˆ ρ ] ( τ ˆ e ) M, β = + * + Trotter break-up: R ˆ ρ R $ = R (e τ H ˆ ) M R $ = dr 1... dr M 1 R e τ ˆ H R 1 R 1 e τ ˆ H R 2... R M 1 e τ ˆ H $ R Path integral and primitive action S : R ˆ ρ $ R = dr t e S[R t ] R R' M (R S[R t ] = i+1 R i ) 2 + τ i=1 4λτ 2 V (R i ) + V (R i +1 ) [ ] Pair action: Militzer, Comp. Phys. Comm. (2016)
13 Douglas Adams: Infinite Improbability Drive of spaceship Heart of Gold
14 Restricted PIMC for fermions: How is the restriction applied? Construct a fermionic trial density matrix in form of a Slater determinant of single-particle density matrices: Imaginary time x ρ T (R, R #,β) = ρ(r 1, r 1 #,β) ρ(r 1, r N #,β) ρ(r N, r 1 #,β) ρ(r N, r N #,β) Enforce the following nodal condition for all time slices along the paths: ρ T [ R(t),R(0),t] > 0 This 3N-dimensional conditions eliminates all negative and some positive contribution to the path Solves the fermion sign problem approx. Free-particle nodes:
15 Starting from Restricted PIMC Simulations of Hydrogen
16 Water
17 DFT-MD Simulations of a Water Plasma at Temperatures of K and Megabar Pressures (P-P 0 )/P 0 (P-P 0 )/P Water Density=3.18 g/cm 3 24-atom cell Water Density=11.18 g/cm 3 24-atom cell Water DFT-MD Temperature (K)
18 Debye Plasma model added Big Gap in Temperature remains (P vs T for a Water Plasma) (P-P 0 )/P 0 (P-P 0 )/P Water Density=3.18 g/cm 3 24-atom cell Water Density=11.18 g/cm 3 24-atom cell Water DFT-MD Debye model Temperature (K)
19 First Path Integral Monte Carlo Simulations for Heavier Elements Fill this Gap in Temperature (P-P 0 )/P 0 (P-P 0 )/P Water Density=3.18 g/cm 3 24-atom cell Water Density=11.18 g/cm 3 24-atom cell Water PIMC DFT-MD Debye model Temperature (K)
20 Study Structural Properties: Pair Correlation Functions for Water and Carbon Plasmas Carbon Density =12.64 g/cm 3 Water Density=11.18 g/cm 3 T=750,000 K PIMC DFT-MD g(r) Water (H-H) Water (O-H) Carbon Water (O-O) r(å)
21 Path Integral Monte Carlo and DFT-MD are in very good agreement Carbon Density =12.64 g/cm 3 Water Density=11.18 g/cm 3 T=750,000 K PIMC DFT-MD g(r) Water (H-H) Water (O-H) Carbon Water (O-O) Driver, Militzer, Phys. Rev. Lett. (2012) r(å)
22 New Path Integral Monte Carlo Simulations of Heavier Elements Aid Fusion Capsule Design ICF Hohlraum ICF Capsule
23 Study planetary interiors in the laboratory: shock wave experiments Two-stage gas gun (Livermore) 0.2 Mbar Z capacitor bank (Sandia) 2 Mbar Nova laser (Livermore) 3.4 Mbar National Ignition Facility 700 Mbar
24 Shock wave measurements determine the Equation of State on the Hugoniot curve Conservation of mass, momentum and energy yields: ρ = ρ 0 u s u s u p P = P 0 + ρ 0 u s u p E = E (V 0 V )(P + P 0 )
25 CH plastics
26 Inertial confinement fusion experiments with plastic coated spheres of liquid H 2 CH H 2 (liquid) (Graphics: Bachmann et al. LLNL)
27 PIMC and DFT-MD simulations performed for C 2 H, CH, C 2 H 3, CH 3 and CH 4.
28 CH Shock Hugoniot Curves: Comparison of Theory and Experiments
29 Silicon
30 Path Integral Monte Carlo with localized nodal surfaces and application to silicon plasmas How the nodes are enforced: Nodes are a Slater determinant: Before we used only free-particle orbitals (plane waves): New idea: Add Hartree-Fock orbitals Militzer, Driver, Phys. Rev. Lett. (2015)
31 Silicon: Energy and Pressure Comparison of Path Integral Monte Carlo and Kohn-Sham DFT Militzer, Driver, Phys. Rev. Lett. (2015)
32 Silicon Hugoniot Curve: Experiments and Semi-analytical EOS models Hu, Militzer, et al. Rev. B 94 (2016)
33 Silicon Hugoniot Curve: Path Integral Monte Carlo, Orbital-Free DFT Hu, Militzer, et al. Rev. B 94 (2016)
34 Silicon Hugoniot Curve: Average-Atom models Hu, Militzer, et al. Rev. B 94 (2016)
35 Comparison of Hugoniot Curves for Different Materials Driver, Militzer, PRB 93, (2016)
36 Comparison of Hugoniot Curves for Different Materials Driver, Militzer, PRB 93, (2016)
37 Comparison of Hugoniot Curves for Different Materials Driver, Militzer, PRB 93, (2016)
38 Comparison of Hugoniot Curves for Different Materials Driver, Militzer, PRB 93, (2016)
39 Comparison of Hugoniot Curves for Different Materials Driver, Militzer, PRB 93, (2016)
40 Conclusions For elements from hydrogen to through silicon, the energy, pressure, and g(r) functions computed with PIMC and KS-DFT agree well at K. Internal energy agrees to better than 5 Ha/atom. Pressure to 2%. So far, we could not detect any problem with the zerotemperature exchange-correlation functionals. We will provide more EOS data to improve orbit-free methods. Second row elements: More comparison between simulations and experiments Extending PIMC to third row elements in particular iron.
41 Conclusions Blue Waters provided a 1 order of magnitude of increase in productivity. Allowed us to solve more complex problems (heavier elements and more accurate nodal surfaces) Enabled us to study a wider parameter space (vary material s composition in addition to density and temperature)
42 The End
All-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 informationDensity Functional Theory Methods for Transport and Optical Properties: Application to Warm Dense Silicon
Density Functional Theory Methods for Transport and Optical Properties: Application to Warm Dense Silicon 2200 Si, T = 62.5 kk K-edge position (ev) 2100 2000 1900 DFT (shifted by 50 ev) AOT Significant
More informationPath integral Monte Carlo simulations of dense carbon-hydrogen plasmas
Path integral Monte Carlo simulations of dense carbon-hydrogen plasmas Shuai Zhang, 1, 2, a) Burkhard Militzer, 1, 3, b) Lorin X. Benedict, 2 François Soubiran, 1 Kevin P. Driver, 1, 2 and Philip A. Sterne
More informationComputational methods: Coupled Electron Ion Monte Carlo Path Integral Monte Carlo Examples: Electron gas Hydrogen
! Computational methods: Coupled Electron Ion Monte Carlo Path Integral Monte Carlo Examples: Electron gas Hydrogen WHO DID THE WORK? Miguel Morales, Livermore Carlo Pierleoni: L Aquila, Italy Jeff McMahon
More informationPath Integral Monte Carlo and Density Functional Molecular Dynamics Simulations of Hot, Dense Helium
Path Integral Monte Carlo and Density Functional Molecular Dynamics Simulations of Hot, Dense Helium B. Militzer Departments of Earth and Planetary Science and Astronomy, University of California, Berkeley,
More informationFPEOS: A First-Principles Equation of State Table of Deuterium for Inertial Confinement Fusion Applications
FPEOS: A First-Principles Equation of State Table of Deuterium for Inertial Confinement Fusion Applications S. X. Hu 1,, B. Militzer 2, V. N. Goncharov 1, S. Skupsky 1 1. Laboratory for Laser Energetics,
More informationPath Integral Monte Carlo Calculation of the Momentum Distribution of the Homogeneous Electron Gas at Finite Temperature
Path Integral Monte Carlo Calculation of the Momentum Distribution of the Homogeneous Electron Gas at Finite Temperature B. Militzer a, E. L. Pollock b, D. M. Ceperley c a Department of Earth and Planetary
More informationCorrelations in Hot Dense Helium
Correlations in Hot Dense Helium Burkhard Militzer University of California, Berkeley, Departments of Earth and Planetary Science and Astronomy, Berkeley, CA 947, USA (Dated: February 4, 9) Hot dense helium
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 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 informationHexatic and microemulsion phases in a 2D quantum Coulomb gas
Hexatic and microemulsion phases in a 2D quantum Coulomb gas Bryan K Clark (University of Illinois at Urbana Champaign) Michele Casula (Ecole Polytechnique, Paris) David M Ceperley (University of Illinois
More informationFilippo Tramonto. Miniworkshop talk: Quantum Monte Carlo simula9ons of low temperature many- body systems
Miniworkshop talk: Quantum Monte Carlo simulations of low temperature many-body systems Physics, Astrophysics and Applied Physics Phd school Supervisor: Dott. Davide E. Galli Outline Interests in quantum
More informationIntroduction to Path Integral Monte Carlo. Part I.
Introduction to Path Integral Monte Carlo. Part I. Alexey Filinov, Jens Böning, Michael Bonitz Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
More informationAb Initio EOS for Planetary Matter and Implications for Giant Planets
Introduction Ab initio EOS H-He phase separation Ab Initio EOS for Planetary Matter and Implications for Giant Planets Winfried Lorenzen, Bastian Holst, Nadine Nettelmann, Ronald Redmer Planet Formation
More informationPath Integral Monte Carlo for Quantum Boltzmannons P R O GRAM I N T B I N ST I T UTE FOR N U C L EAR T H EO RY AU G U ST 2 9,
Path Integral Monte Carlo for Quantum Boltzmannons WILLIAM DAWKINS ADVISOR: ALEXANDROS GEZERLIS UNIVERSITY OF GUELPH P R O GRAM I N T - 18-2B I N ST I T UTE FOR N U C L EAR T H EO RY AU G U ST 2 9, 2 0
More informationComputation of the High Temperature Coulomb Density Matrix in Periodic Boundary Conditions
Computation of the High Temperature Coulomb Density Matrix in Periodic Boundary Conditions B. Militzer Department of Earth and Planetary Science and Department of Astronomy, University of California, Berkeley,
More informationSupplementary Material: Path Integral Monte Carlo Simulation of the Warm-Dense Homogeneous Electron Gas
Supplementary Material: Path Integral Monte Carlo Simulation of the Warm-Dense Homogeneous Electron Gas (Dated: February 26, 2013) In the main text, we describe our calculations of the 3D finite-temperature
More informationChallenges in Path Integral Monte Carlo. David Ceperley Physics UIUC
Challenges in Path Integral Monte Carlo David Ceperley Physics UIUC Imaginary time path integrals Exchange in quantum crystals The fermion sign problem Restricted paths The Penalty method Quantum dynamics?
More informationIntroduction to Density Functional Theory
Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary
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 informationWhere are we with laser fusion?
Where are we with laser fusion? R. Betti Laboratory for Laser Energetics Fusion Science Center Dept. Mechanical Engineering and Physics & Astronomy University of Rochester HEDSA HEDP Summer School August
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Yesterday you (should have) learned this: http://upload.wikimedia.org/wikipedia/commons/e/ea/ Simple_Harmonic_Motion_Orbit.gif 1. deterministic
More informationSimulations of Dense Atomic Hydrogen in the Wigner Crystal Phase
Published in Journal of Physics and Chemistry of Solids, 67 (2006) 2136 Simulations of Dense Atomic Hydrogen in the Wigner Crystal Phase Burkhard Militzer Geophysical Laboratory, Carnegie Institution of
More informationD-D FUSION NEUTRONS FROM A STRONG SPHERICAL SHOCK WAVE FOCUSED ON A DEUTERIUM BUBBLE IN WATER. Dr. Michel Laberge General Fusion Inc.
D-D FUSION NEUTRONS FROM A STRONG SPHERICAL SHOCK WAVE FOCUSED ON A DEUTERIUM BUBBLE IN WATER Dr. Michel Laberge General Fusion Inc. SONOFUSION Sonofusion is making some noise A bit short in energy, ~mj
More informationAn Overview of Quantum Monte Carlo Methods. David M. Ceperley
An Overview of Quantum Monte Carlo Methods David M. Ceperley Department of Physics and National Center for Supercomputing Applications University of Illinois Urbana-Champaign Urbana, Illinois 61801 In
More informationIs the homogeneous electron gas homogeneous?
Is the homogeneous electron gas homogeneous? Electron gas (jellium): simplest way to view a metal homogeneous and normal Hartree-Fock: simplest method for many-electron systems a single Slater determinant
More informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 17 Nov 2004
arxiv:cond-mat/411458v1 [cond-mat.mtrl-sci] 17 Nov 24 Proceedings article of the 5th Conference on Cryocrystals and Quantum Crystals in Wroclaw, Poland, submitted to J. Low. Temp. Phys. (24). Hydrogen-Helium
More informationIon Acceleration from the Interaction of Ultra-Intense Laser Pulse with a Thin Foil
Ion Acceleration from the Interaction of Ultra-Intense Laser Pulse with a Thin Foil Matthew Allen Department of Nuclear Engineering UC Berkeley mallen@nuc.berkeley.edu March 15, 2004 8th Nuclear Energy
More informationDensity Functional Theory (DFT)
Density Functional Theory (DFT) An Introduction by A.I. Al-Sharif Irbid, Aug, 2 nd, 2009 Density Functional Theory Revolutionized our approach to the electronic structure of atoms, molecules and solid
More informationSimulating Quantum High Energy Density Plasmas: Approximations According to the Time-Dependent Variational Principle
Simulating Quantum High Energy Density Plasmas: Approximations According to the Time-Dependent Variational Principle (LANL) with contributions from Andreas Markmann (Yale), Mike Surh (LLNL), Michael S.
More informationKey concepts in Density Functional Theory (I) Silvana Botti
From the many body problem to the Kohn-Sham scheme European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre
More informationOuter Solar System. Jupiter. PHY outer-solar-system - J. Hedberg
Outer Solar System 1. Jupiter 1. Pressure & Density & size 2. Jupiter's Magnetosphere 3. Juno Mission 4. Jovian Satellites 2. Saturn 1. The Rings! 2. Saturn's Moons 3. Titan 3. Uranus 4. Neptune 5. Dwarf
More informationTwo-Dimensional Simulations of Electron Shock Ignition at the Megajoule Scale
Two-Dimensional Simulations of Electron Shock Ignition at the Megajoule Scale Laser intensity ( 1 15 W/cm 2 ) 5 4 3 2 1 Laser spike is replaced with hot-electron spike 2 4 6 8 1 Gain 2 15 1 5 1. 1.2 1.4
More informationHigh-order Chin actions in path integral Monte Carlo
High-order Chin actions in path integral Monte Carlo High-order actions and their applications, Barcelona 2009 Jordi Boronat Departament de Física i Enginyeria Nuclear Universitat Politècnica de Catalunya
More informationComputational Physics. J. M. Thijssen
Computational Physics J. M. Thijssen Delft University of Technology CAMBRIDGE UNIVERSITY PRESS Contents Preface xi 1 Introduction 1 1.1 Physics and computational physics 1 1.2 Classical mechanics and statistical
More informationDensity matrix functional theory vis-á-vis density functional theory
Density matrix functional theory vis-á-vis density functional theory 16.4.007 Ryan Requist Oleg Pankratov 1 Introduction Recently, there has been renewed interest in density matrix functional theory (DMFT)
More informationThe Interior of Giant Planets. Cyrill Milkau
The Interior of Giant Planets Cyrill Milkau 01.12.15 Outline 1. What is a planet? 2. Nuclear fusion 3. Properties of Jupiter 4. Summary 5. Sources Cyrill Milkau 2 1. What is a Planet? Definition by International
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 informationTaxonomy: Computational Methods and Challenges in High Energy Density Physics
LLNL-PRES-537592 Taxonomy: Computational Methods and Challenges in High Energy Density Physics IPAM Computational Methods in High Energy Density Plasmas Tutorials March 13-16, 2012 Frank R. Graziani, P.
More informationExploration of the Feasibility of Polar Drive on the LMJ. Lindsay M. Mitchel. Spencerport High School. Spencerport, New York
Exploration of the Feasibility of Polar Drive on the LMJ Lindsay M. Mitchel Spencerport High School Spencerport, New York Advisor: Dr. R. S. Craxton Laboratory for Laser Energetics University of Rochester
More informationThe National Ignition Campaign: Status and Progress
1 The National Ignition Campaign: Status and Progress E. I. Moses Lawrence Livermore National Laboratory, Livermore, CA 94450 Abstract. The National Ignition Facility (NIF) at Lawrence Livermore National
More informationEquation of state of warm-dense boron nitride combining computation, modeling, and experiment
Equation of state of warm-dense boron nitride combining computation, modeling, and experiment At ambient conditions, BN exists in two stable, nearly degenerate phases: hexagonal BN (h-bn) and cubic BN
More informationWorm Algorithm PIMC Application to Liquid and Solid 4 He
Worm Algorithm PIMC Application to Liquid and Solid 4 He KITP - Santa Barbara, 2/8/06 Massimo Boninsegni University of Alberta Nikolay Prokof ev and Boris Svistunov University of Massachusetts Outline
More informationInertial Confinement Fusion
Inertial Confinement Fusion Prof. Dr. Mathias Groth Aalto University School of Science, Department of Applied Physics Outline Principles of inertial confinement fusion Implosion/compression physics Direct
More informationTutorial on Path Integral Monte Carlo Simulations of Hydrogen and Warm Dense Carbon
Tutorial on Path Integral Monte Carlo Simulations of Hydrogen and Warm Dense Carbon Kevin Driver 1 Burkhard Militzer 1,2 University of California, Berkeley 1 Department of Earth and Planetary Science 2
More informationReferences and Figures from: - Basdevant, Fundamentals in Nuclear Physics
Lecture 22 Fusion Experimental Nuclear Physics PHYS 741 heeger@wisc.edu References and Figures from: - Basdevant, Fundamentals in Nuclear Physics 1 Reading for Next Week Phys. Rev. D 57, 3873-3889 (1998)
More informationThe Interior of Giant Planets
YETI Workshop in Jena, 15-17 November 2010 The Interior of Giant Planets Ronald Redmer Universität Rostock, Institut für Physik D-18051 Rostock, Germany ronald.redmer@uni-rostock.de CO-WORKERS AIMD simulations:
More informationFinite-Temperature Hartree-Fock Exchange and Exchange- Correlation Free Energy Functionals. Travis Sjostrom. IPAM 2012 Workshop IV
1 of 45 Finite-Temperature Hartree-Fock Exchange and Exchange- Correlation Free Energy Functionals Travis Sjostrom Quantum Theory Project Depts. of Physics and Chemistry IPAM 2012 Workshop IV 2012 2 of
More informationFirst-Principles Prediction of the Softening of the Silicon Shock. Hugoniot Curve
First-Principles Prediction of the Softening of the Silicon Shock Hugoniot Curve S. X. Hu ( 胡素兴 ), 1,* B. Militzer, 2,3 L. A. Collins, 4 K. P. Driver 2, and J. D. Kress 4 1 Laboratory for Laser Energetics,
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationCoupled Electron-Ion Monte Carlo of high pressure hydrogen
RPMBT 14; July 16-20, 2007, Barcelona Spain p. 1/6 Coupled Electron-Ion Monte Carlo of high pressure hydrogen Carlo Pierleoni Physics Department, University of L Aquila, Italy RPMBT 14; July 16-20, 2007,
More informationQuantum Monte Carlo. QMC methods in the continuum
Quantum Monte Carlo Premise: need to use simulation techniques to solve manybody quantum problems just as you need them classically. Both the wavefunction and expectation values are determined by the simulations.
More informationDept of Mechanical Engineering MIT Nanoengineering group
1 Dept of Mechanical Engineering MIT Nanoengineering group » To calculate all the properties of a molecule or crystalline system knowing its atomic information: Atomic species Their coordinates The Symmetry
More informationMany electrons: Density functional theory Part II. Bedřich Velický VI.
Many electrons: Density functional theory Part II. Bedřich Velický velicky@karlov.mff.cuni.cz VI. NEVF 514 Surface Physics Winter Term 013-014 Troja 1 st November 013 This class is the second devoted to
More informationHydrodynamic instability measurements in DTlayered ICF capsules using the layered-hgr platform
Journal of Physics: Conference Series PAPER OPEN ACCESS Hydrodynamic instability measurements in DTlayered ICF capsules using the layered-hgr platform Related content - Mix and hydrodynamic instabilities
More informationAnalysis of Experimental Asymmetries using Uncertainty Quantification: Inertial Confinement Fusion (ICF) & its Applications
Analysis of Experimental Asymmetries using Uncertainty Quantification: Inertial Confinement Fusion (ICF) & its Applications Joshua Levin January 9, 2009 (Edited: June 15, 2009) 1 Contents 1. Uncertainty
More informationChapter IX: Nuclear fusion
Chapter IX: Nuclear fusion 1 Summary 1. General remarks 2. Basic processes 3. Characteristics of fusion 4. Solar fusion 5. Controlled fusion 2 General remarks (1) Maximum of binding energy per nucleon
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationHNRS 227 Fall 2006 Chapter 13. What is Pluto? What is a Planet? There are two broad categories of planets: Terrestrial and Jovian
Key Points of Chapter 13 HNRS 227 Fall 2006 Chapter 13 The Solar System presented by Prof. Geller 24 October 2006 Planets Mercury, Venus, Earth, Mars, Jupiter, Saturn, Uranus, Neptune Dwarf Planets Pluto,
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More informationThe Uniform Electron Gas at Warm Dense Matter Conditions
The Uniform Electron Gas at Warm Dense Matter Conditions Tobias Dornheim, Simon Groth, and Michael Bonitz, Physics Reports 744, 1-86 (2018) Institute of Theoretical Physics and Astrophysics Kiel University
More informationSqueezing and superposing many-body states of Bose gases in confining potentials
Squeezing and superposing many-body states of Bose gases in confining potentials K. B. Whaley Department of Chemistry, Kenneth S. Pitzer Center for Theoretical Chemistry, Berkeley Quantum Information and
More informationMethods of Fermion Monte Carlo
Methods of Fermion Monte Carlo Malvin H. Kalos kalos@llnl.gov Institute for Nuclear Theory University of Washington Seattle, USA July, 2013 LLNL-PRES-XXXXXX This work was performed under the auspices of
More informationICF ignition and the Lawson criterion
ICF ignition and the Lawson criterion Riccardo Betti Fusion Science Center Laboratory for Laser Energetics, University of Rochester Seminar Massachusetts Institute of Technology, January 0, 010, Cambridge
More informationarxiv:cond-mat/ v1 5 Aug 2002
Superfluidity in a Doped Helium Droplet E. W. Draeger and D. M. Ceperley Department of Physics and National Center for Supercomputing Applications, University of Illinois Urbana-Champaign, 68 Path Integral
More informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationPath Integral Methods for Continuum Quantum Systems
Path Integral Methods for Continuum Quantum Systems D. Ceperley University of Illinois Urbana-Champaign Imaginary time path integrals Path integrals for bosons Restricted path integral method for fermions
More informationModule 6 1. Density functional theory
Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft
More informationDensity Functional Theory
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function
More informationThe Coupled Electron-Ion Monte Carlo Method and its application to high pressure hydrogen
The Coupled Electron-Ion Monte Carlo Method and its application to high pressure hydrogen Carlo Pierleoni Maison de la Simulation, CEA-Saclay Paris, France and Dept. of Physical and Chemical Sciences University
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 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 informationExchange-Correlation Functional
Exchange-Correlation Functional Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
More informationAyan Chattopadhyay Mainak Mustafi 3 rd yr Undergraduates Integrated MSc Chemistry IIT Kharagpur
Ayan Chattopadhyay Mainak Mustafi 3 rd yr Undergraduates Integrated MSc Chemistry IIT Kharagpur Under the supervision of: Dr. Marcel Nooijen Associate Professor Department of Chemistry University of Waterloo
More informationInertial Confinement Fusion Experiments & Modeling
Inertial Confinement Fusion Experiments & Modeling Using X-ray Absorption Spectroscopy of Thin Tracer Layers to Diagnose the Time-Dependent Properties of ICF Ablator Materials David Cohen (Swarthmore College,
More informationProgress in Warm Dense Matter and Planetary Physics
Progress in Warm Dense Matter and Planetary Physics Winfried Lorenzen, Andreas Becker, and Ronald Redmer Abstract We give an introduction into the method of quantum molecular dynamics simulations which
More informationDiagnosing OMEGA and NIF Implosions Using the D 3 He Spectrum Line Width
Introduction Diagnosing OMEGA and NIF Implosions Using the D 3 He Spectrum Line Width A. B. Zylstra, M. Rosenberg, N. Sinenian, C. Li, F. Seguin, J. Frenje, R. Petrasso (MIT) R. Rygg, D. Hicks, S. Friedrich,
More informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
More informationElectronic Structure Calculations and Density Functional Theory
Electronic Structure Calculations and Density Functional Theory Rodolphe Vuilleumier Pôle de chimie théorique Département de chimie de l ENS CNRS Ecole normale supérieure UPMC Formation ModPhyChem Lyon,
More informationWhat is. Inertial Confinement Fusion?
What is Inertial Confinement Fusion? Inertial Confinement Fusion: dense & short-lived plasma Fusing D and T requires temperature to overcome Coulomb repulsion density & confinement time to maximize number
More informationSolid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
More informationAuxiliary-field quantum Monte Carlo calculations of excited states and strongly correlated systems
Auxiliary-field quantum Monte Carlo calculations of excited states and strongly correlated systems Formally simple -- a framework for going beyond DFT? Random walks in non-orthogonal Slater determinant
More informationBlock Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems
Mitglied der Helmholtz-Gemeinschaft Block Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems Birkbeck University, London, June the 29th 2012 Edoardo Di Napoli Motivation and Goals
More informationNew Frontiers in Nuclear Structure Theory
New Frontiers in Nuclear Structure Theory From Realistic Interactions to the Nuclear Chart Robert Roth Institut für Kernphysik Technical University Darmstadt Overview Motivation Nucleon-Nucleon Interactions
More informationAdvanced Electronic Structure Theory Density functional theory. Dr Fred Manby
Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)
More informationDensity Functional Theory for Electrons in Materials
Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory for
More informationInertial Confinement Fusion DR KATE LANCASTER YORK PLASMA INSTITUTE
Inertial Confinement Fusion DR KATE LANCASTER YORK PLASMA INSTITUTE In the beginning In the late fifties, alternative applications of nuclear explosions were being considered the number one suggestion
More informationarxiv: v1 [astro-ph.ep] 9 Aug 2016
JOURNAL OF GEOPHYSICAL RESEARCH, VOL.???, XXXX, DOI:1.12/, Understanding Jupiter s Interior Burkhard Militzer 1,2, François Soubiran 1, Sean M. Wahl 1, William Hubbard 3 1 Department of Earth and Planetary
More informationIntro to ab initio methods
Lecture 2 Part A Intro to ab initio methods Recommended reading: Leach, Chapters 2 & 3 for QM methods For more QM methods: Essentials of Computational Chemistry by C.J. Cramer, Wiley (2002) 1 ab initio
More informationThe Jovian Planets. Why do we expect planets like this in the outer reaches of the solar system?(lc)
The Jovian Planets Beyond Mars and the Asteroid belt are the Jovian or Gas Giant Planets that are totally different than the terrestrial planets: They are composed almost entirely of gas They do not have
More informationQuantum Monte Carlo Simulations of Exciton Condensates
Quantum Monte Carlo Simulations of Exciton Condensates J. Shumway a and D. M. Ceperley b a Dept. of Physics and Astronomy, Arizona State University, Tempe, AZ 8583 b Dept. of Physics, University of Illinois,
More informationPulsed-power science and engineering, and accelerator design
Department of Energy National Nuclear Security Administration Laboratory Residency Graduate Fellowship (DOE NNSA LRGF) Residency Opportunities at Sandia National Laboratories, New Mexico For more than
More informationThermodynamics of finite quantum systems
Thermodynamics of finite quantum systems Shell structure in finite quantum systems Erice summer school, July 25-30, 2010 Klavs Hansen Department of Physics, University of Gothenburg Topics Metal clusters
More informationTHE NATIONAL IGNITION FACILITY AND THE GOLDEN AGE OF HIGH ENERGY DENSITY SCIENCE
THE NATIONAL IGNITION FACILITY AND THE GOLDEN AGE OF HIGH ENERGY DENSITY SCIENCE Edward I. Moses Lawrence Livermore National Laboratory Livermore, CA 94550 ABSTRACT The National Ignition Facility (NIF)
More informationarxiv: v1 [cond-mat.str-el] 10 Sep 2014
Overcoming the fermion sign problem in homogeneous systems Jonathan L DuBois and Berni J. Alder Lawrence Livermore National Laboratory, Livermore, CA 94550, USA Ethan W. Brown Department of Physics, University
More informationX-ray optics for Laboratory Astrophysics
X-ray optics for Laboratory Astrophysics Abstract: One of the goals of the inertial confinement fusion program at the Lawrence Livermore National Laboratory is to achieve ignition through the indirect
More informationImproving Equation of State Tables
Improving Equation of State Tables Priya Rajasethupathy Advisor: Dr. Jacques Delettrez University of Rochester Laboratory for Laser Energetics Summer High School Research Program 2000 Inertial Confinement
More informationFukuoka, Japan. 23 August National Ignition Facility (NIF) Laboratory for Laser Energetics (OPERA)
Fukuoka, Japan 23 August 2012 National Ignition Facility (NIF) LLNL-PRES-562760 This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under
More informationT T Fusion Neutron Spectrum Measured in Inertial Confinement Fusion Experiment
T T Fusion Neutron Spectrum Measured in Inertial Confinement Fusion Experiment V. Yu. Glebov University of Rochester Laboratory for Laser Energetics 48th Annual Meeting of the American Physical Society
More information