Parallel Tempering Algorithm in Monte Carlo Simulation

Size: px
Start display at page:

Download "Parallel Tempering Algorithm in Monte Carlo Simulation"

Transcription

1 Parallel Tempering Algorithm in Monte Carlo Simulation Tony Cheung (CUHK) Kevin Zhao (CUHK) Mentors: Ying Wai Li (ORNL) Markus Eisenbach (ORNL) Kwai Wong (UTK/ORNL)

2 Metropolis Algorithm on Ising Model Reason: difficulty of direct sampling Objective: compute average physical quantities of interest Idea: generate microstates according to Boltzmann distribution (canonical ensemble) after sufficient number of steps Boltzmann distribution: P s; T = exp ( βe s) Z Underlying principle: detailed balance, β = 1 KT

3 Metropolis Algorithm on Ising Model Simulation process 1. Randomly initialize the model 2. Choose a spin at random & make a trial flip 3. Accept the flip with probability P flip = min*1, exp β E +, β = 1 KT 4. If the flip is accepted, determine the desired physical quantities 5. Repeat steps 2-4 to obtain a sufficient number of microstates 6. Calculate the ensemble average of quantities

4 Parallel Tempering Recall: P flip = min*1, exp β E +, β = 1 Drawback: Low temperature Unlikely to accept flips with positive energy difference Trapped in energy local minimum Motivation: Run Metropolis algorithm on different temperatures & allow exchange of microstates High-temperature configuration at low-temperature system The probability of accepting an exchange is given by P excange = min*1, exp β E + KT

5 Mean Magnetization Per Spin dependence of mean magnetization per spin Total # of MC steps fixed to be 1 9 ; equilibration time set to N=

6 Mean Magnetization Per Spin dependence of mean magnetization per spin Total # of MC steps fixed to be 1 9 ; equilibration time set to N= N=

7 Mean Magnetization Per Spin dependence of mean magnetization per spin Total # of MC steps fixed to be 1 9 ; equilibration time set to N=9 N=

8 Mean Magnetization Per Spin dependence of mean magnetization per spin Total # of MC steps fixed to be 1 9 ; equilibration time set to N=249 N=1^

9 Mean Magnetization Per Spin dependence of mean magnetization per spin Total # of MC steps fixed to be 1 9 ; equilibration time set to N=1^4 N=1^

10 Mean Magnetization Per Spin dependence of mean magnetization per spin Total # of MC steps fixed to be 1 9 ; equilibration time set to N= N=1^

11 Energy Standard Deviation dependence of energy fluctuation (SD) Total # of MC steps fixed to be 1 9 ; equilibration time set to N=

12 Energy Standard Deviation dependence of energy fluctuation (SD) Total # of MC steps fixed to be 1 9 ; equilibration time set to N= N=1^

13 Energy Standard Deviation dependence of energy fluctuation (SD) Total # of MC steps fixed to be 1 9 ; equilibration time set to N=1^4 N=1^

14 Total running time Running Time Dependence on # of Replica Exchanges Monte Carlo Simulation takes up most of the time Parallel tempering becoming the dominant factor ^4 1^5 1^6 1^7 Number of replica exchange (total MC step fixed)

15 Coming soon Parallel tempering Metropolis running with: Various temperature spacing (# of processors) Different exchange patterns Geometric temperature sequence Implementation on other models Goal: optimize the algorithm Better convergence with less time Self adjusting algorithms

Parallel Tempering Algorithm in Monte Carlo Simula5on

Parallel Tempering Algorithm in Monte Carlo Simula5on Parallel Tempering Algorithm in Monte Carlo Simula5on Tony Cheung (CUHK) Kevin Zhao (CUHK) Mentors: Ying Wai Li (ORNL) Markus Eisenbach (ORNL) Kwai Wong (UTK/ORNL) Monte Carlo Algorithms Mo5va5on: Idea:

More information

Parallel Tempering And Adaptive Spacing. For Monte Carlo Simulation

Parallel Tempering And Adaptive Spacing. For Monte Carlo Simulation Parallel Tempering And Adaptive pacing For Monte Carlo imulation Authors: iu Wun Cheung Department of Mathematics, CUHK Zhao, Yiwei Department of Mathematics, CUHK upervisors: Dr. Markus isenbach ORNL

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

Intro. Each particle has energy that we assume to be an integer E i. Any single-particle energy is equally probable for exchange, except zero, assume

Intro. Each particle has energy that we assume to be an integer E i. Any single-particle energy is equally probable for exchange, except zero, assume Intro Take N particles 5 5 5 5 5 5 Each particle has energy that we assume to be an integer E i (above, all have 5) Particle pairs can exchange energy E i! E i +1andE j! E j 1 5 4 5 6 5 5 Total number

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

Graphical Representations and Cluster Algorithms

Graphical Representations and Cluster Algorithms Graphical Representations and Cluster Algorithms Jon Machta University of Massachusetts Amherst Newton Institute, March 27, 2008 Outline Introduction to graphical representations and cluster algorithms

More information

3.320 Lecture 18 (4/12/05)

3.320 Lecture 18 (4/12/05) 3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,

More information

Monte Carlo (MC) Simulation Methods. Elisa Fadda

Monte Carlo (MC) Simulation Methods. Elisa Fadda Monte Carlo (MC) Simulation Methods Elisa Fadda 1011-CH328, Molecular Modelling & Drug Design 2011 Experimental Observables A system observable is a property of the system state. The system state i is

More information

Computational Physics (6810): Session 13

Computational Physics (6810): Session 13 Computational Physics (6810): Session 13 Dick Furnstahl Nuclear Theory Group OSU Physics Department April 14, 2017 6810 Endgame Various recaps and followups Random stuff (like RNGs :) Session 13 stuff

More information

Their Statistical Analvsis. With Web-Based Fortran Code. Berg

Their Statistical Analvsis. With Web-Based Fortran Code. Berg Markov Chain Monter rlo Simulations and Their Statistical Analvsis With Web-Based Fortran Code Bernd A. Berg Florida State Univeisitfi USA World Scientific NEW JERSEY + LONDON " SINGAPORE " BEIJING " SHANGHAI

More information

A Monte Carlo Implementation of the Ising Model in Python

A Monte Carlo Implementation of the Ising Model in Python A Monte Carlo Implementation of the Ising Model in Python Alexey Khorev alexey.s.khorev@gmail.com 2017.08.29 Contents 1 Theory 1 1.1 Introduction...................................... 1 1.2 Model.........................................

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

3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ

3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ 3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple

More information

Understanding temperature and chemical potential using computer simulations

Understanding temperature and chemical potential using computer simulations University of Massachusetts Amherst ScholarWorks@UMass Amherst Physics Department Faculty Publication Series Physics 2005 Understanding temperature and chemical potential using computer simulations J Tobochnik

More information

LECTURE 10: Monte Carlo Methods II

LECTURE 10: Monte Carlo Methods II 1 LECTURE 10: Monte Carlo Methods II In this chapter, we discuss more advanced Monte Carlo techniques, starting with the topics of percolation and random walks, and then continuing to equilibrium statistical

More information

arxiv: v1 [cond-mat.dis-nn] 25 Apr 2018

arxiv: v1 [cond-mat.dis-nn] 25 Apr 2018 arxiv:1804.09453v1 [cond-mat.dis-nn] 25 Apr 2018 Critical properties of the antiferromagnetic Ising model on rewired square lattices Tasrief Surungan 1, Bansawang BJ 1 and Muhammad Yusuf 2 1 Department

More information

Metropolis Monte Carlo simulation of the Ising Model

Metropolis Monte Carlo simulation of the Ising Model Metropolis Monte Carlo simulation of the Ising Model Krishna Shrinivas (CH10B026) Swaroop Ramaswamy (CH10B068) May 10, 2013 Modelling and Simulation of Particulate Processes (CH5012) Introduction The Ising

More information

Monte Caro simulations

Monte Caro simulations Monte Caro simulations Monte Carlo methods - based on random numbers Stanislav Ulam s terminology - his uncle frequented the Casino in Monte Carlo Random (pseudo random) number generator on the computer

More information

arxiv: v1 [cond-mat.dis-nn] 13 Jul 2015

arxiv: v1 [cond-mat.dis-nn] 13 Jul 2015 Spin glass behavior of the antiferromagnetic Heisenberg model on scale free network arxiv:1507.03305v1 [cond-mat.dis-nn] 13 Jul 2015 Tasrief Surungan 1,3, Freddy P. Zen 2,3, and Anthony G. Williams 4 1

More information

Introduction to the Renormalization Group

Introduction to the Renormalization Group Introduction to the Renormalization Group Gregory Petropoulos University of Colorado Boulder March 4, 2015 1 / 17 Summary Flavor of Statistical Physics Universality / Critical Exponents Ising Model Renormalization

More information

Monte-Carlo simulation of small 2D Ising lattice with Metropolis dynamics

Monte-Carlo simulation of small 2D Ising lattice with Metropolis dynamics Monte-Carlo simulation of small 2D Ising lattice with Metropolis dynamics Paul Secular Imperial College London (Dated: 6th February 2015) Results of a Monte-Carlo simulation of the nearest-neighbour Ising

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

Efficiency of Parallel Tempering for Ising Systems

Efficiency of Parallel Tempering for Ising Systems University of Massachusetts Amherst ScholarWorks@UMass Amherst Masters Theses 1911 - February 2014 2010 Efficiency of Parallel Tempering for Ising Systems Stephan Burkhardt University of Massachusetts

More information

Monte Carlo study of the Baxter-Wu model

Monte Carlo study of the Baxter-Wu model Monte Carlo study of the Baxter-Wu model Nir Schreiber and Dr. Joan Adler Monte Carlo study of the Baxter-Wu model p.1/40 Outline Theory of phase transitions, Monte Carlo simulations and finite size scaling

More information

Monte Carlo Simulation of the Ising Model. Abstract

Monte Carlo Simulation of the Ising Model. Abstract Monte Carlo Simulation of the Ising Model Saryu Jindal 1 1 Department of Chemical Engineering and Material Sciences, University of California, Davis, CA 95616 (Dated: June 9, 2007) Abstract This paper

More information

Population Annealing: An Effective Monte Carlo Method for Rough Free Energy Landscapes

Population Annealing: An Effective Monte Carlo Method for Rough Free Energy Landscapes Population Annealing: An Effective Monte Carlo Method for Rough Free Energy Landscapes Jon Machta SFI WORKING PAPER: 21-7-12 SFI Working Papers contain accounts of scientific work of the author(s) and

More information

Simulation of the two-dimensional square-lattice Lenz-Ising model in Python

Simulation of the two-dimensional square-lattice Lenz-Ising model in Python Senior Thesis Simulation of the two-dimensional square-lattice Lenz-Ising model in Python Author: Richard Munroe Advisor: Dr. Edward Brash 2012-10-17 Abstract A program to simulate the the two-dimensional

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

Monte Carlo importance sampling and Markov chain

Monte Carlo importance sampling and Markov chain Monte Carlo importance sampling and Markov chain If a configuration in phase space is denoted by X, the probability for configuration according to Boltzman is ρ(x) e βe(x) β = 1 T (1) How to sample over

More information

Critical Dynamics of Two-Replica Cluster Algorithms

Critical Dynamics of Two-Replica Cluster Algorithms University of Massachusetts Amherst From the SelectedWorks of Jonathan Machta 2001 Critical Dynamics of Two-Replica Cluster Algorithms X. N. Li Jonathan Machta, University of Massachusetts Amherst Available

More information

Electronic Supplementary Information

Electronic Supplementary Information Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2018 Electronic Supplementary Information Self-assembly of conformationaly flexible molecules under

More information

Janus: FPGA Based System for Scientific Computing Filippo Mantovani

Janus: FPGA Based System for Scientific Computing Filippo Mantovani Janus: FPGA Based System for Scientific Computing Filippo Mantovani Physics Department Università degli Studi di Ferrara Ferrara, 28/09/2009 Overview: 1. The physical problem: - Ising model and Spin Glass

More information

4/18/2011. Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania

4/18/2011. Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania 1. Introduction Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania Bibliography Computer experiments Ensemble averages and time averages Molecular

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

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

arxiv: v1 [cond-mat.dis-nn] 13 Oct 2015

arxiv: v1 [cond-mat.dis-nn] 13 Oct 2015 Search for the Heisenberg spin glass on rewired square lattices with antiferromagnetic interaction arxiv:1510.03825v1 [cond-mat.dis-nn] 13 Oct 2015 asrief Surungan, Bansawang BJ, and Dahlang ahir Department

More information

Multicanonical methods

Multicanonical methods Multicanonical methods Normal Monte Carlo algorithms sample configurations with the Boltzmann weight p exp( βe). Sometimes this is not desirable. Example: if a system has a first order phase transitions

More information

Choosing weights for simulated tempering

Choosing weights for simulated tempering PHYSICAL REVIEW E 7, 173 7 Choosing weights for simulated tempering Sanghyun Park 1 and Vijay S. Pande 1, 1 Department of Chemistry, Stanford University, Stanford, California 935, USA Department of Structural

More information

Markov Chain Monte Carlo Method

Markov Chain Monte Carlo Method Markov Chain Monte Carlo Method Macoto Kikuchi Cybermedia Center, Osaka University 6th July 2017 Thermal Simulations 1 Why temperature 2 Statistical mechanics in a nutshell 3 Temperature in computers 4

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

Principles of Equilibrium Statistical Mechanics

Principles of Equilibrium Statistical Mechanics Debashish Chowdhury, Dietrich Stauffer Principles of Equilibrium Statistical Mechanics WILEY-VCH Weinheim New York Chichester Brisbane Singapore Toronto Table of Contents Part I: THERMOSTATICS 1 1 BASIC

More information

Markov Processes. Stochastic process. Markov process

Markov Processes. Stochastic process. Markov process Markov Processes Stochastic process movement through a series of well-defined states in a way that involves some element of randomness for our purposes, states are microstates in the governing ensemble

More information

REVIEW: Derivation of the Mean Field Result

REVIEW: Derivation of the Mean Field Result Lecture 18: Mean Field and Metropolis Ising Model Solutions 1 REVIEW: Derivation of the Mean Field Result The Critical Temperature Dependence The Mean Field Approximation assumes that there will be an

More information

The Role of Limited Valency and Bonding Orientation on Phase Behavior

The Role of Limited Valency and Bonding Orientation on Phase Behavior Wesleyan University The Honors College The Role of Limited Valency and Bonding Orientation on Phase Behavior by Tan Zhi Da Henry Class of 2008 A thesis submitted to the faculty of Wesleyan University in

More information

Monte Carlo Methods in Statistical Mechanics

Monte Carlo Methods in Statistical Mechanics Monte Carlo Methods in Statistical Mechanics Mario G. Del Pópolo Atomistic Simulation Centre School of Mathematics and Physics Queen s University Belfast Belfast Mario G. Del Pópolo Statistical Mechanics

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

Assignment 6: MCMC Simulation and Review Problems

Assignment 6: MCMC Simulation and Review Problems Massachusetts Institute of Technology MITES 2018 Physics III Assignment 6: MCMC Simulation and Review Problems Due Monday July 30 at 11:59PM in Instructor s Inbox Preface: In this assignment, you will

More information

Guiding Monte Carlo Simulations with Machine Learning

Guiding Monte Carlo Simulations with Machine Learning Guiding Monte Carlo Simulations with Machine Learning Yang Qi Department of Physics, Massachusetts Institute of Technology Joining Fudan University in 2017 KITS, June 29 2017. 1/ 33 References Works led

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

TSTC Lectures: Theoretical & Computational Chemistry

TSTC Lectures: Theoretical & Computational Chemistry TSTC Lectures: Theoretical & Computational Chemistry Rigoberto Hernandez May 2009, Lecture #2 Statistical Mechanics: Structure S = k B log W! Boltzmann's headstone at the Vienna Zentralfriedhof.!! Photo

More information

Optimized statistical ensembles for slowly equilibrating classical and quantum systems

Optimized statistical ensembles for slowly equilibrating classical and quantum systems Optimized statistical ensembles for slowly equilibrating classical and quantum systems IPAM, January 2009 Simon Trebst Microsoft Station Q University of California, Santa Barbara Collaborators: David Huse,

More information

Physics 115/242 Monte Carlo simulations in Statistical Physics

Physics 115/242 Monte Carlo simulations in Statistical Physics Physics 115/242 Monte Carlo simulations in Statistical Physics Peter Young (Dated: May 12, 2007) For additional information on the statistical Physics part of this handout, the first two sections, I strongly

More information

The Monte Carlo Method of Random Sampling in Statistical Physics

The Monte Carlo Method of Random Sampling in Statistical Physics IOSR Journal of Applied Physics (IOSR-JAP) e-issn: 2278-4861. Volume 4, Issue 4 (Sep. - Oct. 213), PP 15-23 The Monte Carlo Method of Random Sampling in Statistical Physics 1 Akabuogu E. U., 2 Chiemeka

More information

Two case studies of Monte Carlo simulation on GPU

Two case studies of Monte Carlo simulation on GPU Two case studies of Monte Carlo simulation on GPU National Institute for Computational Sciences University of Tennessee Seminar series on HPC, Feb. 27, 2014 Outline 1 Introduction 2 Discrete energy lattice

More information

Population annealing study of the frustrated Ising antiferromagnet on the stacked triangular lattice

Population annealing study of the frustrated Ising antiferromagnet on the stacked triangular lattice Population annealing study of the frustrated Ising antiferromagnet on the stacked triangular lattice Michal Borovský Department of Theoretical Physics and Astrophysics, University of P. J. Šafárik in Košice,

More information

The Phase Transition of the 2D-Ising Model

The Phase Transition of the 2D-Ising Model The Phase Transition of the 2D-Ising Model Lilian Witthauer and Manuel Dieterle Summer Term 2007 Contents 1 2D-Ising Model 2 1.1 Calculation of the Physical Quantities............... 2 2 Location of the

More information

Hanoi 7/11/2018. Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam.

Hanoi 7/11/2018. Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam. Hanoi 7/11/2018 Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam. Finite size effects and Reweighting methods 1. Finite size effects 2. Single histogram method 3. Multiple histogram method 4. Wang-Landau

More information

Trivially parallel computing

Trivially parallel computing Parallel Computing After briefly discussing the often neglected, but in praxis frequently encountered, issue of trivially parallel computing, we turn to parallel computing with information exchange. Our

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

J ij S i S j B i S i (1)

J ij S i S j B i S i (1) LECTURE 18 The Ising Model (References: Kerson Huang, Statistical Mechanics, Wiley and Sons (1963) and Colin Thompson, Mathematical Statistical Mechanics, Princeton Univ. Press (1972)). One of the simplest

More information

Equilibrium, out of equilibrium and consequences

Equilibrium, out of equilibrium and consequences Equilibrium, out of equilibrium and consequences Katarzyna Sznajd-Weron Institute of Physics Wroc law University of Technology, Poland SF-MTPT Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium

More information

Progress toward a Monte Carlo Simulation of the Ice VI-VII Phase Transition

Progress toward a Monte Carlo Simulation of the Ice VI-VII Phase Transition Progress toward a Monte Carlo Simulation of the Ice VI-VII Phase Transition Christina Gower 2010 NSF/REU PROJECT Physics Department University of Notre Dame Advisor: Dr. Kathie E. Newman August 6, 2010

More information

Invaded cluster dynamics for frustrated models

Invaded cluster dynamics for frustrated models PHYSICAL REVIEW E VOLUME 57, NUMBER 1 JANUARY 1998 Invaded cluster dynamics for frustrated models Giancarlo Franzese, 1, * Vittorio Cataudella, 1, * and Antonio Coniglio 1,2, * 1 INFM, Unità di Napoli,

More information

Classical Monte-Carlo simulations

Classical Monte-Carlo simulations Classical Monte-Carlo simulations Graduate Summer Institute Complex Plasmas at the Stevens Insitute of Technology Henning Baumgartner, A. Filinov, H. Kählert, P. Ludwig and M. Bonitz Christian-Albrechts-University

More information

Numerical Analysis of 2-D Ising Model. Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011

Numerical Analysis of 2-D Ising Model. Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011 Numerical Analysis of 2-D Ising Model By Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011 Contents Abstract Acknowledgment Introduction Computational techniques Numerical Analysis

More information

NATIONAL UNIVERSITY OF SINGAPORE PC5215 NUMERICAL RECIPES WITH APPLICATIONS. (Semester I: AY ) Time Allowed: 2 Hours

NATIONAL UNIVERSITY OF SINGAPORE PC5215 NUMERICAL RECIPES WITH APPLICATIONS. (Semester I: AY ) Time Allowed: 2 Hours NATIONAL UNIVERSITY OF SINGAPORE PC55 NUMERICAL RECIPES WITH APPLICATIONS (Semester I: AY 0-) Time Allowed: Hours INSTRUCTIONS TO CANDIDATES. This examination paper contains FIVE questions and comprises

More information

Markov Chain Monte Carlo Lecture 6

Markov Chain Monte Carlo Lecture 6 Sequential parallel tempering With the development of science and technology, we more and more need to deal with high dimensional systems. For example, we need to align a group of protein or DNA sequences

More information

Lab 70 in TFFM08. Curie & Ising

Lab 70 in TFFM08. Curie & Ising IFM The Department of Physics, Chemistry and Biology Lab 70 in TFFM08 Curie & Ising NAME PERS. -NUMBER DATE APPROVED Rev Aug 09 Agne 1 Introduction Magnetic materials are all around us, and understanding

More information

3. General properties of phase transitions and the Landau theory

3. General properties of phase transitions and the Landau theory 3. General properties of phase transitions and the Landau theory In this Section we review the general properties and the terminology used to characterise phase transitions, which you will have already

More information

c 2007 by Harvey Gould and Jan Tobochnik 28 May 2007

c 2007 by Harvey Gould and Jan Tobochnik 28 May 2007 Chapter 5 Magnetic Systems c 2007 by Harvey Gould and Jan Tobochnik 28 May 2007 We apply the general formalism of statistical mechanics developed in Chapter 4 to the Ising model, a model magnetic system

More information

The XY-Model. David-Alexander Robinson Sch th January 2012

The XY-Model. David-Alexander Robinson Sch th January 2012 The XY-Model David-Alexander Robinson Sch. 08332461 17th January 2012 Contents 1 Introduction & Theory 2 1.1 The XY-Model............................... 2 1.2 Markov Chains...............................

More information

arxiv:cond-mat/ v4 [cond-mat.dis-nn] 23 May 2001

arxiv:cond-mat/ v4 [cond-mat.dis-nn] 23 May 2001 Phase Diagram of the three-dimensional Gaussian andom Field Ising Model: A Monte Carlo enormalization Group Study arxiv:cond-mat/488v4 [cond-mat.dis-nn] 3 May M. Itakura JS Domestic esearch Fellow, Center

More information

Teaching Statistical and Thermal Physics Using Computer Simulations

Teaching Statistical and Thermal Physics Using Computer Simulations Teaching Statistical and Thermal Physics Using Computer Simulations Tutorial T2, 4 March 2007 Harvey Gould, Clark University Collaborators: Wolfgang Christian, Davidson College Jan Tobochnik,

More information

A New Algorithm for Monte Carlo Simulation of king Spin Systems* A. B. BORTZ M. H. KALOS J. L. LEBOWITZ

A New Algorithm for Monte Carlo Simulation of king Spin Systems* A. B. BORTZ M. H. KALOS J. L. LEBOWITZ JOURNAL OF COMPUTATIONAL PHYSICS 17, lo-18 (1975) A New Algorithm for Monte Carlo Simulation of king Spin Systems* A. B. BORTZ Berfer Graduate School of Science, Yeshiva University, New York, New York

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

Random Walks A&T and F&S 3.1.2

Random Walks A&T and F&S 3.1.2 Random Walks A&T 110-123 and F&S 3.1.2 As we explained last time, it is very difficult to sample directly a general probability distribution. - If we sample from another distribution, the overlap will

More information

Paramagnetic phases of Kagome lattice quantum Ising models p.1/16

Paramagnetic phases of Kagome lattice quantum Ising models p.1/16 Paramagnetic phases of Kagome lattice quantum Ising models Predrag Nikolić In collaboration with T. Senthil Massachusetts Institute of Technology Paramagnetic phases of Kagome lattice quantum Ising models

More information

Statistical Mechanics Victor Naden Robinson vlnr500 3 rd Year MPhys 17/2/12 Lectured by Rex Godby

Statistical Mechanics Victor Naden Robinson vlnr500 3 rd Year MPhys 17/2/12 Lectured by Rex Godby Statistical Mechanics Victor Naden Robinson vlnr500 3 rd Year MPhys 17/2/12 Lectured by Rex Godby Lecture 1: Probabilities Lecture 2: Microstates for system of N harmonic oscillators Lecture 3: More Thermodynamics,

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

2. Verification of the Boltzmann distribution

2. Verification of the Boltzmann distribution Exercises Lecture VIII Macrostates and microstates: equilibrium and entropy. Metropolis algorithm in the canonical ensemble: verifying Boltzmann s distribution 1. MC simulation of a simple N-particles

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

arxiv: v1 [cond-mat.stat-mech] 11 Oct 2014

arxiv: v1 [cond-mat.stat-mech] 11 Oct 2014 Heuristic Monte Carlo Method Applied to Cooperative Motion Algorithm for Binary Lattice Fluid P. Knycha la 1, and M. Banaszak 2 1 Chair of Computer Science, The President Stanislaw arxiv:1410.3015v1 [cond-mat.stat-mech]

More information

ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below

ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940

More information

Optimized multicanonical simulations: A proposal based on classical fluctuation theory

Optimized multicanonical simulations: A proposal based on classical fluctuation theory Optimized multicanonical simulations: A proposal based on classical fluctuation theory J. Viana Lopes Centro de Física do Porto and Departamento de Física, Faculdade de Ciências, Universidade do Porto,

More information

arxiv: v3 [cond-mat.stat-mech] 20 Nov 2014

arxiv: v3 [cond-mat.stat-mech] 20 Nov 2014 Effective ergodicity in single-spin-flip dynamics Mehmet Süzen Applied Mathematical Physiology Lab Bonn University, Sigmund-Freud-Str.25, 53127 Bonn, Germany (Dated: March 9, 2018) arxiv:1405.4497v3 [cond-mat.stat-mech]

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

Statistics and Quantum Computing

Statistics and Quantum Computing Statistics and Quantum Computing Yazhen Wang Department of Statistics University of Wisconsin-Madison http://www.stat.wisc.edu/ yzwang Workshop on Quantum Computing and Its Application George Washington

More information

Monte Carlo Methods in High Energy Physics I

Monte Carlo Methods in High Energy Physics I Helmholtz International Workshop -- CALC 2009, July 10--20, Dubna Monte Carlo Methods in High Energy Physics CALC2009 - July 20 10, Dubna 2 Contents 3 Introduction Simple definition: A Monte Carlo technique

More information

Distributed Processing of the Lattice in Monte Carlo Simulations of the Ising Type Spin Model

Distributed Processing of the Lattice in Monte Carlo Simulations of the Ising Type Spin Model CMST 21(3) 117-121 (2015) DOI:10.12921/cmst.2015.21.03.002 Distributed Processing of the Lattice in Monte Carlo Simulations of the Ising Type Spin Model Szymon Murawski *, Grzegorz Musiał, Grzegorz Pawłowski

More information

Surface effects in frustrated magnetic materials: phase transition and spin resistivity

Surface effects in frustrated magnetic materials: phase transition and spin resistivity Surface effects in frustrated magnetic materials: phase transition and spin resistivity H T Diep (lptm, ucp) in collaboration with Yann Magnin, V. T. Ngo, K. Akabli Plan: I. Introduction II. Surface spin-waves,

More information

Markov Chain Monte Carlo The Metropolis-Hastings Algorithm

Markov Chain Monte Carlo The Metropolis-Hastings Algorithm Markov Chain Monte Carlo The Metropolis-Hastings Algorithm Anthony Trubiano April 11th, 2018 1 Introduction Markov Chain Monte Carlo (MCMC) methods are a class of algorithms for sampling from a probability

More information

4. Cluster update algorithms

4. Cluster update algorithms 4. Cluster update algorithms Cluster update algorithms are the most succesful global update methods in use. These methods update the variables globally, in one step, whereas the standard local methods

More information

Parallel Tempering Algorithm for. Conformational Studies of Biological Molecules

Parallel Tempering Algorithm for. Conformational Studies of Biological Molecules Parallel Tempering Algorithm for Conformational Studies of Biological Molecules Ulrich H.E. Hansmann 1 Department of Theoretical Studies, Institute for Molecular Science Okazaki, Aichi 444, Japan ABSTRACT

More information

Metropolis/Variational Monte Carlo. Microscopic Theories of Nuclear Structure, Dynamics and Electroweak Currents June 12-30, 2017, ECT*, Trento, Italy

Metropolis/Variational Monte Carlo. Microscopic Theories of Nuclear Structure, Dynamics and Electroweak Currents June 12-30, 2017, ECT*, Trento, Italy Metropolis/Variational Monte Carlo Stefano Gandolfi Los Alamos National Laboratory (LANL) Microscopic Theories of Nuclear Structure, Dynamics and Electroweak Currents June 12-30, 2017, ECT*, Trento, Italy

More information

Equilibrium Molecular Thermodynamics from Kirkwood Sampling

Equilibrium Molecular Thermodynamics from Kirkwood Sampling This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source

More information

Title Theory of solutions in the energy r of the molecular flexibility Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2003), 9702 Issue Date 2003-11-08 URL http://hdl.handle.net/2433/50354

More information

Cluster Monte Carlo study of multicomponent fluids of the Stillinger-Helfand and Widom- Rowlinson type

Cluster Monte Carlo study of multicomponent fluids of the Stillinger-Helfand and Widom- Rowlinson type University of Massachusetts Amherst From the SelectedWorks of Jonathan Machta 2000 Cluster Monte Carlo study of multicomponent fluids of the Stillinger-Helfand and Widom- Rowlinson type Rongfeng Sun Harvey

More information

arxiv: v2 [cond-mat.stat-mech] 23 Mar 2009

arxiv: v2 [cond-mat.stat-mech] 23 Mar 2009 arxiv:0806.1054v2 [cond-mat.stat-mech] 23 Mar 2009 The ground state energy of the Edwards-Anderson spin glass model with a parallel tempering Monte Carlo algorithm F. Romá 1,2, S. Risau-Gusman 1, A. J.

More information

Equilibrium Molecular Thermodynamics from Kirkwood Sampling

Equilibrium Molecular Thermodynamics from Kirkwood Sampling This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source

More information

The Monte Carlo Method An Introduction to Markov Chain MC

The Monte Carlo Method An Introduction to Markov Chain MC The Monte Carlo Method An Introduction to Markov Chain MC Sachin Shanbhag Department of Scientific Computing, Florida State University Sources. Introduction to Monte Carlo Algorithms : Werner Krauth. Available

More information