Ab initio statistical mechanics

Size: px
Start display at page:

Download "Ab initio statistical mechanics"

Transcription

1 Ab initio statistical mechanics Luca M. Ghiringhelli Fritz-Haber-Institut der MPG, Berlin Hands-on Workshop Density-Functional Theory and Beyond: Accuracy, Efficiency and Reproducibility in Computational Materials Science Humboldt University, Berlin, Germany, July 31 to August 11, 2017

2 Extending the scale

3 Extending the scale

4 Chemical energy conversion: catalysis

5 Entropy?

6 Free energy: one quantity, many definitions

7 Free energy: one quantity, many definitions

8 Free energy: one quantity, many definitions

9 Statistical mechanics: free energy as a probabilistic concept

10 Statistical mechanics, quantities derived from Z

11 Statistical mechanics, quantities derived from Z Evaluation of pressure

12 Ensemble averages on discrete machines

13 The problem of free energy sampling

14 Computational free-energy evaluation: the zoo - Analytic: ab initio atomistic thermodynamics - Canonical sampling: thermodynamic integration - Canonical sampling: thermodynamic perturbation - Generalized sampling: biased sampling / biased dynamics - Unbiased (canonical) sampling re-weighting techniques - Evaluation: Parallel or >>> Serial <<<

15 Computational free-energy evaluation: the zoo - Analytic: ab initio atomistic thermodynamics - Canonical sampling: thermodynamic integration Lecture Lectureof of Sergey SergeyLevchenko, Levchenko,th Tuesday, Tuesday,August August88th,, at at10:00 10:00 - Canonical sampling: thermodynamic perturbation - Generalized sampling: biased sampling / biased dynamics - Unbiased (canonical) sampling re-weighting techniques - Evaluation: Parallel or >>> Serial <<<

16 Computational free-energy evaluation: the zoo - Analytic: ab initio atomistic thermodynamics - Canonical sampling: thermodynamic integration - Canonical sampling: thermodynamic perturbation - Generalized sampling: biased sampling / biased dynamics - Unbiased (canonical) sampling re-weighting techniques - Evaluation: Parallel or >>> Serial <<<

17 Free energy: physical path thermodynamic integration

18 Free energy: unphysical path thermodynamic integration

19 Case study: phase diagram of pure carbon

20 Case study: phase diagram of pure carbon

21 Case study: phase diagram of pure carbon

22 Case study: phase diagram of pure carbon

23 Case study: λ ensemble sampling and integration

24 Case study: integration of P(ρ) equations of state

25 Case study: equating Gibbs free energies

26 Carbon phase diagram

27 Alternative method for finding phase coexistence via F(V)

28 Notable cases (at 0 K): Silicon (1980)

29 Notable cases (at 0 K): Cerium (2013) Casadei et al. PRL (2013)

30 Computational free-energy evaluation: the zoo - Analytic: ab initio atomistic thermodynamics - Canonical sampling: thermodynamic integration - Canonical sampling: thermodynamic perturbation - Generalized sampling: biased sampling / biased dynamics - Unbiased (canonical) sampling re-weighting techniques - Evaluation: Parallel or >>> Serial <<<

31 Thermodynamic perturbation

32 Thermodynamic perturbation

33 Thermodynamic perturbation

34 Thermodynamic perturbation

35 Thermodynamic perturbation

36 Thermodynamic perturbation: recycling data

37 Thermodynamic perturbation: recycling data

38 Computational free-energy evaluation: the zoo - Analytic: ab initio atomistic thermodynamics - Canonical sampling: thermodynamic integration - Canonical sampling: thermodynamic perturbation - Generalized sampling: biased sampling / biased dynamics - Unbiased (canonical) sampling re-weighting techniques - Evaluation: Parallel or >>> Serial <<<

39 Umbrella sampling

40 Umbrella sampling (multiplying by 1 few more times...)

41 Umbrella sampling (multiplying by 1 few more times...) Parallel (over biasing potentials)

42 Umbrella sampling

43 Computational free-energy evaluation: the zoo - Analytic: ab initio atomistic thermodynamics - Canonical sampling: thermodynamic integration - Canonical sampling: thermodynamic perturbation - Generalized sampling: biased sampling / biased dynamics - Unbiased (canonical) sampling re-weighting techniques - Evaluation: Parallel or >>> Serial <<<

44 Parallel tempering: the concept

45 Parallel tempering: the implementation Parallel (over temperatures, and more) To be tuned for efficient sampling: number of temperatures, list of temperatures, attempted swap frequency

46 Parallel tempering: free energy?

47 Au4: coexistence of several isomers

48 Au4: coexistence of several isomers

49 Relative population of larger clusters: Au n, 2D vs 3D

50 The Nested Sampling Obtaining the partition function Consider cumulative density E

51 The Nested Sampling: the main trick Instead of χ(e), compute E(χ) χ(e) E(χ) χ E At E =, we have an ideal gas, χ0 = VN Constrained uniform sampling E(χ) E1 E χ1 χ0 J. Skilling (2006) G. Csányi (2011)

52 The Nested Sampling: the main trick E(χ) E1 χ1 χ0 1. obtain K uniform samples such that E(q) < Elimit 2. compute median: E(χ1) = E1, χ1 χ0/2, Elimit E1 3. repeat... E(χ) E 2 E 1 χ2 χ1 χ0 (linked with LAMMPS), Csányi et al.

53 The Nested Sampling: application to (EAM) Aluminum Temperature [K] Parallel over walkers Pressure [GPa] R. Baldock,, G. Csányi, PRB 93, (2016)

54 Computational free-energy evaluation: the zoo - Analytic: ab initio atomistic thermodynamics - Canonical sampling: thermodynamic integration - Canonical sampling: thermodynamic perturbation - Generalized sampling: biased sampling / biased dynamics - Unbiased (canonical) sampling (replica exchange, nested sampling) - Evaluation: Parallel or >>> Serial <<<

Computing free energy: Replica exchange

Computing free energy: Replica exchange Computing free energy: Replica exchange Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls Macroscopic i a t e regime

More information

Free metal and metal oxide clusters: beyond the static, monostructure description

Free metal and metal oxide clusters: beyond the static, monostructure description Free metal and metal oxide clusters: beyond the static, monostructure description Luca M. Ghiringhelli Fritz Haber Institute of the Max Planck Society, Berlin IPAM Workshop: Materials for a Sustainable

More information

Computing free energy: Thermodynamic perturbation and beyond

Computing free energy: Thermodynamic perturbation and beyond Computing free energy: Thermodynamic perturbation and beyond Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls

More information

Monte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky

Monte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)

More information

Lecture 14: Advanced Conformational Sampling

Lecture 14: Advanced Conformational Sampling Lecture 14: Advanced Conformational Sampling Dr. Ronald M. Levy ronlevy@temple.edu Multidimensional Rough Energy Landscapes MD ~ ns, conformational motion in macromolecules ~µs to sec Interconversions

More information

(# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble

(# = %(& )(* +,(- 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 information

Free energy calculations and the potential of mean force

Free energy calculations and the potential of mean force Free energy calculations and the potential of mean force IMA Workshop on Classical and Quantum Approaches in Molecular Modeling Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science

More information

Ab initio molecular dynamics and nuclear quantum effects

Ab initio molecular dynamics and nuclear quantum effects Ab initio molecular dynamics and nuclear quantum effects Luca M. Ghiringhelli Fritz Haber Institute Hands on workshop density functional theory and beyond: First principles simulations of molecules and

More information

Monte Carlo Methods. Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects

Monte Carlo Methods. Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects Monte Carlo Methods Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects Lecture 1 2 Lecture 1&2 3 Lecture 1&3 4 Different Ensembles Ensemble ame Constant (Imposed) VT Canonical,V,T P PT

More information

Advanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)

Advanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds) Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is

More information

Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry

Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore

More information

Combining quasiparticle energy calculations with exact-exchange density-functional theory

Combining quasiparticle energy calculations with exact-exchange density-functional theory Combining quasiparticle energy calculations with exact-exchange density-functional theory Patrick Rinke 1, Abdallah Qteish 1,2, Jörg Neugebauer 1,3,4, Christoph Freysoldt 1 and Matthias Scheffler 1 1 Fritz-Haber-Institut

More information

Road map for tutorial 5: Big-data analytics

Road map for tutorial 5: Big-data analytics Road map for tutorial 5: Big-data analytics Compressed sensing and screening for properties Emre Ahmetcik, Angelo Ziletti, and Luca M. Ghiringhelli Fritz-Haber-Institut der MPG, Berlin Hands-on Workshop

More information

Multi-Ensemble Markov Models and TRAM. Fabian Paul 21-Feb-2018

Multi-Ensemble Markov Models and TRAM. Fabian Paul 21-Feb-2018 Multi-Ensemble Markov Models and TRAM Fabian Paul 21-Feb-2018 Outline Free energies Simulation types Boltzmann reweighting Umbrella sampling multi-temperature simulation accelerated MD Analysis methods

More information

Copyright 2001 University of Cambridge. Not to be quoted or copied without permission.

Copyright 2001 University of Cambridge. Not to be quoted or copied without permission. Course MP3 Lecture 4 13/11/2006 Monte Carlo method I An introduction to the use of the Monte Carlo method in materials modelling Dr James Elliott 4.1 Why Monte Carlo? The name derives from the association

More information

1. Thermodynamics 1.1. A macroscopic view of matter

1. Thermodynamics 1.1. A macroscopic view of matter 1. Thermodynamics 1.1. A macroscopic view of matter Intensive: independent of the amount of substance, e.g. temperature,pressure. Extensive: depends on the amount of substance, e.g. internal energy, enthalpy.

More information

(i) T, p, N Gibbs free energy G (ii) T, p, µ no thermodynamic potential, since T, p, µ are not independent of each other (iii) S, p, N Enthalpy H

(i) T, p, N Gibbs free energy G (ii) T, p, µ no thermodynamic potential, since T, p, µ are not independent of each other (iii) S, p, N Enthalpy H Solutions exam 2 roblem 1 a Which of those quantities defines a thermodynamic potential Why? 2 points i T, p, N Gibbs free energy G ii T, p, µ no thermodynamic potential, since T, p, µ are not independent

More information

Nanoscale simulation lectures Statistical Mechanics

Nanoscale simulation lectures Statistical Mechanics Nanoscale simulation lectures 2008 Lectures: Thursdays 4 to 6 PM Course contents: - Thermodynamics and statistical mechanics - Structure and scattering - Mean-field approaches - Inhomogeneous systems -

More information

Ab Initio Atomistic Thermodynamics

Ab Initio Atomistic Thermodynamics Ab Initio Atomistic Thermodynamics (Elizabeth C. Beret) Luca M. Ghiringhelli Fritz-Haber-Institut Berlin Extending the scale Thermodynamics: p, T, V, N Length (m) continuum 0 e or m -3 Potential Energy

More information

Stability, Composition and Function of Palladium Surfaces in Oxidizing Environments: A First-Principles Statistical Mechanics Approach

Stability, Composition and Function of Palladium Surfaces in Oxidizing Environments: A First-Principles Statistical Mechanics Approach Stability, Composition and Function of Palladium Surfaces in Oxidizing Environments: A First-Principles Statistical Mechanics Approach von Diplom-Chemikerin Jutta Rogal im Fachbereich Physik der Freien

More information

Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany

Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics

More information

Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces

Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces Minh Triet Dang, 1 Ana Vila Verde, 2 Van Duc Nguyen, 1 Peter G. Bolhuis, 3 and Peter Schall 1

More information

Monte Carlo strategies for first-principles simulations of elemental systems

Monte Carlo strategies for first-principles simulations of elemental systems Monte Carlo strategies for first-principles simulations of elemental systems Lev Gelb Department of Materials Science and Engineering, University of Texas at Dallas XSEDE12 Lev Gelb (UT Dallas) Monte Carlo

More information

A Brief Introduction to Statistical Mechanics

A Brief Introduction to Statistical Mechanics A Brief Introduction to Statistical Mechanics E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop Universidade

More information

Au-C Au-Au. g(r) r/a. Supplementary Figures

Au-C Au-Au. g(r) r/a. Supplementary Figures g(r) Supplementary Figures 60 50 40 30 20 10 0 Au-C Au-Au 2 4 r/a 6 8 Supplementary Figure 1 Radial bond distributions for Au-C and Au-Au bond. The zero density regime between the first two peaks in g

More information

Free energy calculations using molecular dynamics simulations. Anna Johansson

Free energy calculations using molecular dynamics simulations. Anna Johansson Free energy calculations using molecular dynamics simulations Anna Johansson 2007-03-13 Outline Introduction to concepts Why is free energy important? Calculating free energy using MD Thermodynamical Integration

More information

Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany

Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical

More information

Collective behavior, from particles to fields

Collective behavior, from particles to fields 978-0-51-87341-3 - Statistical Physics of Fields 1 Collective behavior, from particles to fields 1.1 Introduction One of the most successful aspects of physics in the twentieth century was revealing the

More information

UNIVERSITY OF OSLO FACULTY OF MATHEMATICS AND NATURAL SCIENCES

UNIVERSITY OF OSLO FACULTY OF MATHEMATICS AND NATURAL SCIENCES UNIVERSITY OF OSLO FCULTY OF MTHEMTICS ND NTURL SCIENCES Exam in: FYS430, Statistical Mechanics Day of exam: Jun.6. 203 Problem :. The relative fluctuations in an extensive quantity, like the energy, depends

More information

Lecture 2+3: Simulations of Soft Matter. 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications

Lecture 2+3: Simulations of Soft Matter. 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications Lecture 2+3: Simulations of Soft Matter 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications D. Frenkel, Boulder, July 6, 2006 What distinguishes Colloids from atoms or

More information

(De)-localization in mean-field quantum glasses

(De)-localization in mean-field quantum glasses QMATH13, Georgia Tech, Atlanta October 10, 2016 (De)-localization in mean-field quantum glasses Chris R. Laumann (Boston U) Chris Baldwin (UW/BU) Arijeet Pal (Oxford) Antonello Scardicchio (ICTP) CRL,

More information

Multicanonical parallel tempering

Multicanonical parallel tempering JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 13 1 APRIL 2002 Multicanonical parallel tempering Roland Faller, Qiliang Yan, and Juan J. de Pablo Department of Chemical Engineering, University of Wisconsin,

More information

Grand Canonical Formalism

Grand Canonical Formalism Grand Canonical Formalism Grand Canonical Ensebmle For the gases of ideal Bosons and Fermions each single-particle mode behaves almost like an independent subsystem, with the only reservation that the

More information

exciting in a nutshell

exciting in a nutshell http://exciting-code.org exciting in a nutshell Pasquale Pavone Humboldt-Universität zu Berlin http://exciting-code.org exciting in a (coco)nutshell Pasquale Pavone Humboldt-Universität zu Berlin Outline

More information

Introduction to Restricted Boltzmann Machines

Introduction to Restricted Boltzmann Machines Introduction to Restricted Boltzmann Machines Ilija Bogunovic and Edo Collins EPFL {ilija.bogunovic,edo.collins}@epfl.ch October 13, 2014 Introduction Ingredients: 1. Probabilistic graphical models (undirected,

More information

Statistical Mechanics

Statistical Mechanics 42 My God, He Plays Dice! Statistical Mechanics Statistical Mechanics 43 Statistical Mechanics Statistical mechanics and thermodynamics are nineteenthcentury classical physics, but they contain the seeds

More information

PHYS 352 Homework 2 Solutions

PHYS 352 Homework 2 Solutions PHYS 352 Homework 2 Solutions Aaron Mowitz (, 2, and 3) and Nachi Stern (4 and 5) Problem The purpose of doing a Legendre transform is to change a function of one or more variables into a function of variables

More information

Electron-phonon scattering from green s function transport combined with molecular dynamics: Applications to mobility predictions

Electron-phonon scattering from green s function transport combined with molecular dynamics: Applications to mobility predictions Electron-phonon scattering from green s function transport combined with molecular dynamics: Applications to mobility predictions Daniele Stradi www.quantumwise.com daniele.stradi@quantumwise.com Introduction

More information

Statistical Mechanics. Atomistic view of Materials

Statistical Mechanics. Atomistic view of Materials Statistical Mechanics Atomistic view of Materials What is statistical mechanics? Microscopic (atoms, electrons, etc.) Statistical mechanics Macroscopic (Thermodynamics) Sample with constrains Fixed thermodynamics

More information

Thermal pure quantum state

Thermal pure quantum state Thermal pure quantum state Sho Sugiura ( 杉浦祥 ) Institute for Solid State Physics, Univ. Tokyo Collaborator: Akira Shimizu (Univ. Tokyo) SS and A.Shimizu, PRL 108, 240401 (2012) SS and A.Shimizu, PRL 111,

More information

MONTE CARLO METHODS IN SEQUENTIAL AND PARALLEL COMPUTING OF 2D AND 3D ISING MODEL

MONTE CARLO METHODS IN SEQUENTIAL AND PARALLEL COMPUTING OF 2D AND 3D ISING MODEL Journal of Optoelectronics and Advanced Materials Vol. 5, No. 4, December 003, p. 971-976 MONTE CARLO METHODS IN SEQUENTIAL AND PARALLEL COMPUTING OF D AND 3D ISING MODEL M. Diaconu *, R. Puscasu, A. Stancu

More information

The Second Virial Coefficient & van der Waals Equation

The Second Virial Coefficient & van der Waals Equation V.C The Second Virial Coefficient & van der Waals Equation Let us study the second virial coefficient B, for a typical gas using eq.v.33). As discussed before, the two-body potential is characterized by

More information

Critical like behavior in the lattice gas model

Critical like behavior in the lattice gas model Critical like behavior in the lattice gas model A. Wielocha, J. Brzychczyka, J. Łukasikb, P. Pawłowskib, T. Pietrzaka, W. Trautmannc a M. Smoluchowski Institute of Physics, Jagiellonian University, Reymonta

More information

Contents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21

Contents. 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 information

The Electronic Structure of Dye- Sensitized TiO 2 Clusters from Many- Body Perturbation Theory

The Electronic Structure of Dye- Sensitized TiO 2 Clusters from Many- Body Perturbation Theory The Electronic Structure of Dye- Sensitized TiO 2 Clusters from Many- Body Perturbation Theory Noa Marom Center for Computational Materials Institute for Computational Engineering and Sciences The University

More information

Phys Midterm. March 17

Phys Midterm. March 17 Phys 7230 Midterm March 17 Consider a spin 1/2 particle fixed in space in the presence of magnetic field H he energy E of such a system can take one of the two values given by E s = µhs, where µ is the

More information

Why thermodynamics for materials?

Why thermodynamics for materials? Why thermodynamics for materials? Example p 2mkT T For = 300 K, = 1 atm ~ 10 8 site -1 s -1 p p Requires 10-12 atm to keep a clean surface clean; surface can also lose atoms Example Thermodynamic potentials

More information

6.730 Physics for Solid State Applications

6.730 Physics for Solid State Applications 6.730 Physics for Solid State Applications Lecture 25: Chemical Potential and Equilibrium Outline Microstates and Counting System and Reservoir Microstates Constants in Equilibrium Temperature & Chemical

More information

Rustam Z. Khaliullin University of Zürich

Rustam 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 information

Free energy simulations

Free energy simulations Free energy simulations Marcus Elstner and Tomáš Kubař January 14, 2013 Motivation a physical quantity that is of most interest in chemistry? free energies Helmholtz F or Gibbs G holy grail of computational

More information

Computer Simulation of Shock Waves in Condensed Matter. Matthew R. Farrow 2 November 2007

Computer Simulation of Shock Waves in Condensed Matter. Matthew R. Farrow 2 November 2007 Computer Simulation of Shock Waves in Condensed Matter Matthew R. Farrow 2 November 2007 Outline of talk Shock wave theory Results Conclusion Computer simulation of shock waves Shock Wave Theory Shock

More information

An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation

An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation J. Comput. Chem. Jpn., Vol. 8, o. 3, pp. 97 14 (9) c 9 Society of Computer Chemistry, Japan An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation Yosuke KATAOKA* and Yuri YAMADA

More information

Bimodality and latent heat of gold nuclei

Bimodality and latent heat of gold nuclei Bimodality and latent heat of gold nuclei Eric Bonnet (GANIL) for the INDRA and ALADIN collaboration IWM 09, Catania Outline First order phase transition in finite system Definition and related signals

More information

Table of Contents [ttc]

Table of Contents [ttc] Table of Contents [ttc] 1. Equilibrium Thermodynamics I: Introduction Thermodynamics overview. [tln2] Preliminary list of state variables. [tln1] Physical constants. [tsl47] Equations of state. [tln78]

More information

Pressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster

Pressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster University of Rhode Island DigitalCommons@URI Chemistry Faculty Publications Chemistry 2005 Pressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster Dubravko Sabo University

More information

PH4211 Statistical Mechanics Brian Cowan

PH4211 Statistical Mechanics Brian Cowan PH4211 Statistical Mechanics Brian Cowan Contents 1 The Methodology of Statistical Mechanics 1.1 Terminology and Methodology 1.1.1 Approaches to the subject 1.1.2 Description of states 1.1.3 Extensivity

More information

Accurate van der Waals interactions from ground state electron density

Accurate van der Waals interactions from ground state electron density Accurate van der Waals interactions from ground state electron density Alexandre Tkatchenko Theory Department, Fritz Haber Institut der MPG Berlin, Germany tkatchen@fhi berlin.mpg.de Haber Institute EXCITCM09,

More information

The PLUMED plugin and free energy methods in electronic-structure-based molecular dynamics

The PLUMED plugin and free energy methods in electronic-structure-based molecular dynamics The PLUMED plugin and free energy methods in electronic-structure-based molecular dynamics Davide Branduardi, Theoretical Molecular Biophysics Group Max Planck for Biophysics, Frankfurt am Main, Germany

More information

Lecture 8: Computer Simulations of Generalized Ensembles

Lecture 8: Computer Simulations of Generalized Ensembles Lecture 8: Computer Simulations of Generalized Ensembles Bernd A. Berg Florida State University November 6, 2008 Bernd A. Berg (FSU) Generalized Ensembles November 6, 2008 1 / 33 Overview 1. Reweighting

More information

V.C The Second Virial Coefficient & van der Waals Equation

V.C The Second Virial Coefficient & van der Waals Equation V.C The Second Virial Coefficient & van der Waals Equation Let us study the second virial coefficient B, for a typical gas using eq.(v.33). As discussed before, the two-body potential is characterized

More information

Skyrmions à la carte

Skyrmions à la carte This project has received funding from the European Union's Horizon 2020 research and innovation programme FET under grant agreement No 665095 Bertrand Dupé Institute of Theoretical Physics and Astrophysics,

More information

Lecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods

Lecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods Statistical Thermodynamics Lecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods Dr. Ronald M. Levy ronlevy@temple.edu Free energy calculations Free energy

More information

Lecture 20. Chemical Potential

Lecture 20. Chemical Potential Lecture 20 Chemical Potential Reading: Lecture 20, today: Chapter 10, sections A and B Lecture 21, Wednesday: Chapter 10: 10 17 end 3/21/16 1 Pop Question 7 Boltzmann Distribution Two systems with lowest

More information

Random processes and probability distributions. Phys 420/580 Lecture 20

Random processes and probability distributions. Phys 420/580 Lecture 20 Random processes and probability distributions Phys 420/580 Lecture 20 Random processes Many physical processes are random in character: e.g., nuclear decay (Poisson distributed event count) P (k, τ) =

More information

Road map for tutorial 8: Big-data analytics Compressed sensing and screening for properties

Road map for tutorial 8: Big-data analytics Compressed sensing and screening for properties Road map for tutorial 8: Big-data analytics Compressed sensing and screening for properties Emre Ahmetcik, Angelo Ziletti, and Luca M. Ghiringhelli Fritz-Haber-Institut der MPG, Berlin Hands-On DFT and

More information

Lecture Topic 4: Chapter 7 Sampling and Sampling Distributions

Lecture Topic 4: Chapter 7 Sampling and Sampling Distributions Lecture Topic 4: Chapter 7 Sampling and Sampling Distributions Statistical Inference: The aim is to obtain information about a population from information contained in a sample. A population is the set

More information

Full counting statistics in quantum thermodynamics Absorption refrigerator at strong coupling

Full counting statistics in quantum thermodynamics Absorption refrigerator at strong coupling Full counting statistics in quantum thermodynamics Absorption refrigerator at strong coupling Dvira Segal Department of Chemistry University of Toronto Publications Method Full counting statistics for

More information

Equilibrium Study of Spin-Glass Models: Monte Carlo Simulation and Multi-variate Analysis. Koji Hukushima

Equilibrium Study of Spin-Glass Models: Monte Carlo Simulation and Multi-variate Analysis. Koji Hukushima Equilibrium Study of Spin-Glass Models: Monte Carlo Simulation and Multi-variate Analysis Koji Hukushima Department of Basic Science, University of Tokyo mailto:hukusima@phys.c.u-tokyo.ac.jp July 11 15,

More information

Methods of Data Analysis Learning probability distributions

Methods of Data Analysis Learning probability distributions Methods of Data Analysis Learning probability distributions Week 5 1 Motivation One of the key problems in non-parametric data analysis is to infer good models of probability distributions, assuming we

More information

Thermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning

Thermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning Thermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning Ben Leimkuhler University of Edinburgh Joint work with C. Matthews (Chicago), G. Stoltz (ENPC-Paris), M. Tretyakov

More information

Thermodynamics (Lecture Notes) Heat and Thermodynamics (7 th Edition) by Mark W. Zemansky & Richard H. Dittman

Thermodynamics (Lecture Notes) Heat and Thermodynamics (7 th Edition) by Mark W. Zemansky & Richard H. Dittman Thermodynamics (Lecture Notes Heat and Thermodynamics (7 th Edition by Mark W. Zemansky & Richard H. Dittman 2 Chapter 1 Temperature and the Zeroth Law of Thermodynamics 1.1 Macroscopic Point of View If

More information

Long-Term Atomistic Simulation of Heat and Mass Transport

Long-Term Atomistic Simulation of Heat and Mass Transport Long-Term Atomistic Simulation of Heat and Mass Transport Kevin G. Wang 1, Mauricio Ponga 2, Michael Ortiz 2 1 Department of Aerospace and Ocean Engineering, Virginia Polytechnic Institute and State University

More information

Solving the sign problem for a class of frustrated antiferromagnets

Solving the sign problem for a class of frustrated antiferromagnets Solving the sign problem for a class of frustrated antiferromagnets Fabien Alet Laboratoire de Physique Théorique Toulouse with : Kedar Damle (TIFR Mumbai), Sumiran Pujari (Toulouse Kentucky TIFR Mumbai)

More information

of Long-Term Transport

of Long-Term Transport Atomistic Modeling and Simulation of Long-Term Transport Phenomena in Nanomaterials K.G. Wang, M. and M. Ortiz Virginia Polytechnic Institute and State University Universidad de Sevilla California Institute

More information

Supporting Information. Dissociative Water Adsorption by Al 3 O 4 + in the Gas Phase. Linnéstrasse 2, D Leipzig, Germany.

Supporting Information. Dissociative Water Adsorption by Al 3 O 4 + in the Gas Phase. Linnéstrasse 2, D Leipzig, Germany. Supporting Information Dissociative Water Adsorption by Al 3 O 4 + in the Gas Phase Matias R. Fagiani, 1,2, Xiaowei Song, 1,2 Sreekanta Debnath, 1,2 Sandy Gewinner, 2 Wieland Schöllkopf, 2 Knut R. Asmis,

More information

Phase Transitions. µ a (P c (T ), T ) µ b (P c (T ), T ), (3) µ a (P, T c (P )) µ b (P, T c (P )). (4)

Phase Transitions. µ a (P c (T ), T ) µ b (P c (T ), T ), (3) µ a (P, T c (P )) µ b (P, T c (P )). (4) Phase Transitions A homogeneous equilibrium state of matter is the most natural one, given the fact that the interparticle interactions are translationally invariant. Nevertheless there is no contradiction

More information

EQUATION OF STATE AND FLUCTUATIONS FROM THE LATTICE Claudia Ratti University of Houston (USA)

EQUATION OF STATE AND FLUCTUATIONS FROM THE LATTICE Claudia Ratti University of Houston (USA) EQUATION OF STATE AND FLUCTUATIONS FROM THE LATTICE Claudia Ratti University of Houston (USA) Collaborators: Paolo Alba, Rene Bellwied, Szabolcs Borsanyi, Zoltan Fodor, Jana Guenther, Sandor Katz, Stefan

More information

PHY 5524: Statistical Mechanics, Spring February 11 th, 2013 Midterm Exam # 1

PHY 5524: Statistical Mechanics, Spring February 11 th, 2013 Midterm Exam # 1 PHY 554: Statistical Mechanics, Spring 013 February 11 th, 013 Midterm Exam # 1 Always remember to write full work for what you do. This will help your grade in case of incomplete or wrong answers. Also,

More information

The nature of superfluidity in the cold atomic unitary Fermi gas

The nature of superfluidity in the cold atomic unitary Fermi gas The nature of superfluidity in the cold atomic unitary Fermi gas Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo (AFMC) method The trapped unitary Fermi gas

More information

U = γkt N. is indeed an exact one form. Thus there exists a function U(S, V, N) so that. du = T ds pdv. U S = T, U V = p

U = γkt N. is indeed an exact one form. Thus there exists a function U(S, V, N) so that. du = T ds pdv. U S = T, U V = p A short introduction to statistical mechanics Thermodynamics consists of the application of two principles, the first law which says that heat is energy hence the law of conservation of energy must apply

More information

Sampling the free energy surfaces of collective variables

Sampling the free energy surfaces of collective variables Sampling the free energy surfaces of collective variables Jérôme Hénin Enhanced Sampling and Free-Energy Calculations Urbana, 12 September 2018 Please interrupt! struct bioinform phys chem theoretical

More information

fiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES

fiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES Content-Thermodynamics & Statistical Mechanics 1. Kinetic theory of gases..(1-13) 1.1 Basic assumption of kinetic theory 1.1.1 Pressure exerted by a gas 1.2 Gas Law for Ideal gases: 1.2.1 Boyle s Law 1.2.2

More information

Ab-initio molecular dynamics for High pressure Hydrogen

Ab-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 information

Information Theory and Predictability Lecture 6: Maximum Entropy Techniques

Information Theory and Predictability Lecture 6: Maximum Entropy Techniques Information Theory and Predictability Lecture 6: Maximum Entropy Techniques 1 Philosophy Often with random variables of high dimensional systems it is difficult to deduce the appropriate probability distribution

More information

1 The fundamental equation of equilibrium statistical mechanics. 3 General overview on the method of ensembles 10

1 The fundamental equation of equilibrium statistical mechanics. 3 General overview on the method of ensembles 10 Contents EQUILIBRIUM STATISTICAL MECHANICS 1 The fundamental equation of equilibrium statistical mechanics 2 2 Conjugate representations 6 3 General overview on the method of ensembles 10 4 The relation

More information

2.4 Quantum confined electrons

2.4 Quantum confined electrons 2.4. Quantum confined electrons 5 2.4 Quantum confined electrons We will now focus our attention on the electron charge densities in case of one, two and three-dimensional confinement. All the relations

More information

V.E Mean Field Theory of Condensation

V.E Mean Field Theory of Condensation V.E Mean Field heory of Condensation In principle, all properties of the interacting system, including phase separation, are contained within the thermodynamic potentials that can be obtained by evaluating

More information

Design of Oscillator Networks for Generating Signal with Prescribed Statistical Property

Design of Oscillator Networks for Generating Signal with Prescribed Statistical Property Journal of Physics: Conference Series PAPER OPEN ACCESS Design of Oscillator Networks for Generating Signal with Prescribed Statistical Property To cite this article: Tatsuo Yanagita 2017 J. Phys.: Conf.

More information

( ) ( )( k B ( ) ( ) ( ) ( ) ( ) ( k T B ) 2 = ε F. ( ) π 2. ( ) 1+ π 2. ( k T B ) 2 = 2 3 Nε 1+ π 2. ( ) = ε /( ε 0 ).

( ) ( )( k B ( ) ( ) ( ) ( ) ( ) ( k T B ) 2 = ε F. ( ) π 2. ( ) 1+ π 2. ( k T B ) 2 = 2 3 Nε 1+ π 2. ( ) = ε /( ε 0 ). PHYS47-Statistical Mechanics and hermal Physics all 7 Assignment #5 Due on November, 7 Problem ) Problem 4 Chapter 9 points) his problem consider a system with density of state D / ) A) o find the ermi

More information

Protection of excited spin states by a superconducting energy gap

Protection of excited spin states by a superconducting energy gap Protection of excited spin states by a superconducting energy gap B. W. Heinrich, 1 L. Braun, 1, J. I. Pascual, 1, 2, 3 and K. J. Franke 1 1 Institut für Experimentalphysik, Freie Universität Berlin, Arnimallee

More information

Metadynamics-biased ab initio Molecular Dynamics Study of

Metadynamics-biased ab initio Molecular Dynamics Study of Metadynamics-biased ab initio Molecular Dynamics Study of Heterogeneous CO 2 Reduction via Surface Frustrated Lewis Pairs Mireille Ghoussoub 1, Shwetank Yadav 1, Kulbir Kaur Ghuman 1, Geoffrey A. Ozin

More information

Long-Term Atomistic Simulation of Hydrogen Diffusion in Metals

Long-Term Atomistic Simulation of Hydrogen Diffusion in Metals Long-Term Atomistic Simulation of Hydrogen Diffusion in Metals Kevin Wang (1), Pilar Ariza (2), and Michael Ortiz (3) (1) Department of Aerospace and Ocean Engineering Virginia Polytechnic Institute and

More information

Molecular Modeling of Matter

Molecular Modeling of Matter Molecular Modeling of Matter Keith E. Gubbins Lecture 1: Introduction to Statistical Mechanics and Molecular Simulation Common Assumptions Can treat kinetic energy of molecular motion and potential energy

More information

2m + U( q i), (IV.26) i=1

2m + U( q i), (IV.26) i=1 I.D The Ideal Gas As discussed in chapter II, micro-states of a gas of N particles correspond to points { p i, q i }, in the 6N-dimensional phase space. Ignoring the potential energy of interactions, the

More information

Biomolecular modeling. Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010)

Biomolecular modeling. Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010) Biomolecular modeling Marcus Elstner and Tomáš Kubař Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010) IX. FREE ENERGY SIMULATIONS When searching for a physical quantity that is of most

More information

Information, Physics, and Computation

Information, Physics, and Computation Information, Physics, and Computation Marc Mezard Laboratoire de Physique Thdorique et Moales Statistiques, CNRS, and Universit y Paris Sud Andrea Montanari Department of Electrical Engineering and Department

More information

Optimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics

Optimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics Optimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics Sangeeta Sharma 1,2, J. K. Dewhurst 3, C. Ambrosch-Draxl 4, S. Pittalis 2, S. Kurth 2, N. Helbig 2, S. Shallcross

More information

On the number of circuits in random graphs. Guilhem Semerjian. [ joint work with Enzo Marinari and Rémi Monasson ]

On the number of circuits in random graphs. Guilhem Semerjian. [ joint work with Enzo Marinari and Rémi Monasson ] On the number of circuits in random graphs Guilhem Semerjian [ joint work with Enzo Marinari and Rémi Monasson ] [ Europhys. Lett. 73, 8 (2006) ] and [ cond-mat/0603657 ] Orsay 13-04-2006 Outline of the

More information

Elementary Lectures in Statistical Mechanics

Elementary Lectures in Statistical Mechanics George DJ. Phillies Elementary Lectures in Statistical Mechanics With 51 Illustrations Springer Contents Preface References v vii I Fundamentals: Separable Classical Systems 1 Lecture 1. Introduction 3

More information

CHEM-UA 652: Thermodynamics and Kinetics

CHEM-UA 652: Thermodynamics and Kinetics 1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 4 I. THE ISOTHERMAL-ISOBARIC ENSEMBLE The isothermal-isobaric ensemble is the closest mimic to the conditions under which most experiments are

More information