Hanoi 7/11/2018. Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam.
|
|
- Marlene Webb
- 5 years ago
- Views:
Transcription
1 Hanoi 7/11/2018 Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam.
2 Finite size effects and Reweighting methods 1. Finite size effects 2. Single histogram method 3. Multiple histogram method 4. Wang-Landau method 5. The applications
3 A Guide to Monte Carlo Simulations in Statistical Physics D. Landau and K. Binder, (Cambridge University Press, 2009). Monte Carlo Simulation in Statistical Physics: An Introduction K. Binder and D. W. Heermann (Springer-Verlag Berlin Heidelberg, 2010). Understanding Molecular Simulation : From Algorithms to Applications D. Frenkel, (Academic Press, 2002). Frustrated Spin Systems H. T. Diep, 2 nd Ed. (World Scientific, 2013). Lecture notes Example code :
4 Phase transitions finite size effects Second order Critical exponents First order Double peak histogram Single histogram Multiple histogram Flat histogram
5 1. Finite size effects Using for determining the properties of the corresponding infinite system distinguish the order of the transition (first or second order transition) test the results of simulation of finite size systems Finite size scaling and critical exponents Consider a spin model on the lattice of finite size L Scaling forms of thermodynamic quantities (3.1) : are scaling functions Scaling relations (3.2)
6 The thermodynamic properties at the transition reduce to proportionality constants Scaling relations (3.3) If, then The cumulants fourth order cumulant of the order parameter (3.4) As the system size For For
7 fourth order cumulant of the energy (3.5) The logarithmic derivative of the n th power of the magnetization (3.6) For a continuous transition (second order transition) (3.7) where
8 Calculate the critical exponents exponent : (3.8) exponent : (3.9) exponent : exponent : (3.10) (3.11) For a discontinuous transition (first order transition) finite-size scaling laws (3.12) with
9 Ferromagnetic Ising spin model on simple cubic lattice
10 Ferromagnetic Ising spin model on simple cubic lattice
11 The behavior of V L for the q = 10 Potts model in two dimensions A Guide to Monte Carlo Simulations in Statistical Physics D. Landau and K. Binder
12 2. Single histogram method Introduced by Ferrenberg and Swendsen (1988), Ferrenberg (1991) using histograms to extract information from Monte Carlo simulations Applying to the calculation of critical exponents The theory Consider a Monte Carlo simulation performed at T = T 0 generates system configurations with a frequency proportional to the Boltzmann weight, exp(-e/k B T). Histogram of energy and magnetization H(E, M) The probability of simultaneously observing the system : is the number of configurations (density of states) with energy E and magnetization M (3.13) Partition function: (3.14)
13 : provides an estimate for over the range of E and M values generated during the simulation N is the number of measurements made We have (3.15) : is an estimate for the true density of states From the histogram H(E, M) we can invert Eq. (3.15) to determine (3.16) Replace in Eq. (3.13) with the expression for, and normalize the distribution. the relationship between the histogram measured at K = K 0 and the estimated probability distribution for arbitrary K (3.17)
14 The average value of any function of E and M : (3.18) Energy histogram vs energy for ferromagnetic Ising spin model on square lattice of size L = 64, at temperature K =
15 Case study - one dimension single histogram From the histogram of energy H(E) The estimated probability distribution for arbitrary K (3.19) The average value of any function of E (3.20) For example
16 The simulations Find the temperature T 0 which is close to transition temperature Simulate the model by using the standard Monte Carlo simulation (Metropolis algorithm). Plot specific heat or susceptibility versus temperature, peaking the value of temperature at maxima of or. Choose the range of energy T 0 = 2.26 E min = -1.8, E max = -1.0
17 Do the simulation by using the single histogram method at T = T 0 Program structure - Gererating an array of N_BIN for energy range of interest - initialize the lattice SUBROUTINE equilibrating! Equilibrating process DO ieq = 1 to N_eq CALL monte_carlo_step() ENDDO END SUBROUTINE SUBROUTINE averaging! Averaging process DO iav = 1 to N_av CALL monte_carlo_step() do-analysis ENDDO END SUBROUTINE SUBROUTINE monte_carlo_step generate-one-sweep END SUBROUTINE
18 Source : single_histogram_ising.f90 PROGRAM SINGLE_HISTOGRAM_ISING IMPLICIT NONE INTEGER, PARAMETER :: n = 40,nms = n*n REAL,DIMENSION(n,n) :: s REAL,ALLOCATABLE,DIMENSION(:) :: e1,e2,m1,m2,he INTEGER :: i,j,ip,im,jp,jm,ie,ia,ms,neq,nav,ib,nbin REAL :: r,emin,emax,t,dt,temp,et1,et2! CALL RANDOM_SEED() neq = nav = t = 2.3; dt = 4.0 emax = -1.0*float(n*n); emin = -2.0*float(n*n) nbin = INT((emax-emin)/dt) + 1 ALLOCATE(he(nbin),e1(nbin),e2(nbin),m1(nbin),m2(nbin))
19 CALL spin_conf()! initial configuration CALL equilibrating() he =0.0; e1 = 0.0; e2 = 0.0; m1 = 0.0; m2 = 0.0 CALL averaging() OPEN(UNIT=12,FILE='sh_ising.dat') i = 0 DO ib = 1,nbin IF (he(ib) > 0) THEN write(12,1) he(ib),e1(ib)/he(ib),e2(ib)/he(ib)&,m1(ib)/he(ib),m2(ib)/he(ib) i = i + 1 ENDIF ENDDO 1 format(5(f16.6,1x)) CLOSE(12) OPEN(UNIT=12,FILE='sh_ising.info') write(12,*) t,i,n CLOSE(12) DEALLOCATE(he,e1,e2,m1,m2)! CONTAINS! list of subroutine/function
20 SUBROUTINE spin_conf! CALL RANDOM_NUMBER(s) DO j = 1,n DO i = 1,n IF (s(i,j) > 0.5) THEN s(i,j) = 1.0 ELSE s(i,j) = -1.0 ENDIF ENDDO ENDDO END SUBROUTINE spin_conf! SUBROUTINE equilibrating DO ie = 1, neq CALL monte_carlo_step() ENDDO END SUBROUTINE equilibrating
21 SUBROUTINE averaging DO ia = 1, nav CALL monte_carlo_step() temp = 0.0 DO i = 1,n ip = i + 1; im = i - 1 IF (i == n) ip = 1 IF (i == 1) im = n DO j = 1,n jp = j + 1; jm = j - 1 IF (j == n) jp = 1 IF (j == 1) jm = n temp = temp & - s(i,j)*(s(ip,j)+s(im,j)+s(i,jp)+s(i,jm)) ENDDO ENDDO temp = 0.5*temp ib = NINT((temp-emin)/dt) + 1
22 IF ((ib > 0).AND. (ib <= nbin)) THEN he(ib) = he(ib) temp = temp/float(n*n) e1(ib) = e1(ib) + temp e2(ib) = e2(ib) + temp*temp! temp = 0.0 DO j = 1,n DO i = 1,n temp = temp + S(i,j) ENDDO ENDDO temp = temp/float(n*n) m1(ib) = m1(ib)+ abs(temp) m2(ib) = m2(ib)+ temp*temp ENDIF ENDDO END SUBROUTINE averaging
23 SUBROUTINE monte_carlo_step DO ms = 1,nms CALL RANDOM_NUMBER(r) i = INT(r*float(n))+1 CALL RANDOM_NUMBER(r) j = INT(r*float(n))+1 ip = i + 1; im = i - 1 IF (i == n) ip = 1 IF (i == 1) im = n jp = j + 1; jm = j - 1 IF (j == n) jp = 1 IF (j == 1) jm = n et1 = -s(i,j)*(s(ip,j)+s(im,j)+s(i,jp)+s(i,jm)) et2 = -et1 CALL RANDOM_NUMBER(r) IF (r < exp(-(et2-et1)/t)) s(i,j) = -s(i,j) ENDDO END SUBROUTINE monte_carlo_step! END PROGRAM SINGLE_HISTOGRAM_ISING
24 Do analysis Read the data from histogram MC simulation Calculate the partition function Calculate the average value of physical quantities
25 Source : single_histogram_analysis.f90 PROGRAM SINGLE_HISTOGRAM_ANALYSIS IMPLICIT NONE REAL,ALLOCATABLE,DIMENSION(:) :: e,e1,e2,m1,m2,he INTEGER :: n,ib,nbin,it,nt REAL :: r,tmin,tmax,t,t0,dt,z REAL :: et1,et2,mt1,mt2,cv,chi! nt = 80 tmin = 2.1 tmax = 2.4 dt = (tmax-tmin)/float(nt-1)!
26 OPEN(UNIT=12,FILE='sh_ising.info') read(12,*) t0,nbin,n CLOSE(12) ALLOCATE(he(nbin),e(nbin),e1(nbin)&,e2(nbin),m1(nbin),m2(nbin))! OPEN(UNIT=12,FILE='sh_ising.dat') DO ib = 1,nbin read(12,1) he(ib),e1(ib),e2(ib),m1(ib),m2(ib) ENDDO 1 format(5(f16.6,1x)) CLOSE(12) e = e1*float(n*n)! OPEN(UNIT=12,FILE='ising_analysis.dat') t = tmin
27 DO it = 1,nt z = 0.0 DO ib = 1,nbin z = z + he(ib)*exp(e(ib)*(1/t0-1/t)) ENDDO et1 = 0.0; et2 = 0.0; mt1 = 0.0; mt2 = 0.0 DO ib = 1, nbin et1 = et1 + he(ib)*e1(ib)*exp(e(ib)*(1/t0-1/t)) et2 = et2 + he(ib)*e2(ib)*exp(e(ib)*(1/t0-1/t)) mt1 = mt1 + he(ib)*m1(ib)*exp(e(ib)*(1/t0-1/t)) mt2 = mt2 + he(ib)*m2(ib)*exp(e(ib)*(1/t0-1/t)) ENDDO et1 = et1/z; et2 = et2/z; mt1 = mt1/z; mt2 = mt2/z cv = float(n*n)*(et2-et1*et1)/t/t chi = float(n*n)*(mt2-mt1*mt1)/t write(12,1) t,et1,mt1,cv,chi t = t + dt ENDDO CLOSE(12) DEALLOCATE(he,e,e1,e2,m1,m2) END PROGRAM SINGLE_HISTOGRAM_ANALYSIS
28 Ferromagnetic Ising spin model on the square lattice Specific heat and susceptibility vs temperature with N = 40 red points : Standard MC result green points : Single histogram result
29 Ferromagnetic Ising spin model on the square lattice Specific heat and susceptibility vs temperature with N = 40, 50, 60
30 3. Multiple histogram method Introduced by Ferrenberg and Swendsen (Phys. Rev. Lett. 63,1195, 1989) Developed from single histogram method The technique is known to reproduce with very high accuracy the critical exponents of second order phase transitions The theory Consider n independent MC simulation Each simulation at temperature T j, the system has N j configurations The overall probability distribution at temperature T (3.21) where (3.22) f i is chosen self-consistently using Eq. 3.21
31 The thermal average of a physical quantity A is then calculated by (3.23) Monte Carlo simulations Choose the range of temperature [T min, T max ] Divide the teperature into n points: T i (i = 1,, n) Choose the range of energy [E min, E max ] Divide the energy into nbin Perform n independent MC simulation Calculate the energy histogram, average energy, magnetization Perform the analysis Read the simulation data Calculate self-consistently f i by iterating Eqs. (3.21) and (3.22) Do the calculations for the physical quantities by using (3.23)
32 4. Wang-Landau method Introduced by Wang and Landau, Phys. Rev. Lett. 86, 2050 (2001) It permits one to detect with efficiency weak first-order transitions Use to study classical statistical models with difficultly accessed microscopic states. The theory Wang Landau sampling The algorithm uses a random walk in energy space to obtain an accurate estimate for the density of states g(e) g(e) is defined as the number of spin configurations for any given E The classical partition function can either be written as a sum over all states or over all energies (3.24) g(e) independent of temperature, it can be used to find all properties of the system at all temperatures.
33 We begin with some simple guess for the density of states: g(e) = 1 For improve g(e), spins are overturned according to the probability (3.25) Following each spin-flip trial the density of states is updateds f i is a modification factor that is initially greater than 1 (3.26) At the beginning of the random walk, the modification factor f can be as large as A histogram H(E) records the number of times a state of energy E is visited Each time the energy histogram satisfies a certain flatness criterion, f is reduced according to and reset histogram H(E) = 0 for all E
34 The reduction process of the modification factor f is repeated several times until the final value f final is close enough to 1.0 The histogram is considered as flat if x% is flatness criterion, it can be chosen between 70% and 95% is the average (mean value) histogram. The thermodynamic quantities can be evaluated by (3.27) (3.28) (3.29) where Z is the partition function defined by (3.29) The canonical distribution at a temperature T (3.30)
35 Wang-Landau Monte Carlo scheme (1) Set g(e) = 1; choose a modification factor (e.g. f 0 = e 1 ) (2) Choose an initial state (3) Choose a site i (4) Calculate the ratio of the density of states (5) Generate a random number r such that 0 < r < 1 (6) If r < : then flip the spin (7) Set (8) If the histogram is not flat, go to the next site and go to (4) (9) If the histogram is flat, decrease f, e.g. (10) Repeat steps (3) (9) until (11) Calculate properties using final density of states g(e)
36 Critical behavior of magnetic thin films Study the critical behavior of magnetic thin films as a function of the film thickness Lattice: simple cubic Spin: Ising Methods: Standard Monte Carlo method Multiple histogram method Model We consider a thin film made made from a ferromagnetic simple cubic lattice The size of the film is L L N z Hamiltonian The periodic boundary conditions (PBC) is applied in the xy planes X.T. Pham Phu, V. T. Ngo and H. T. Diep, Surface Science 603, (2009).
37 Monte Carlo results Magnetization and susceptibility with
38 Susceptibility and V 1 with
39 Critical exponents Two-dimensional Ising model = 1
40 Critical exponents Two-dimensional Ising model = 1.75
41 Critical exponent 3D = Three-dimensional Ising model = 0.63
42 Critical exponent 3D = 1.25 Three-dimension Ising model = 1.24
43 Critical exponent vs film thickness
44 Critical exponent vs film thickness
45 From, we calculate the effective dimension of thin film
46 Stacked triangular antiferromagnets Study the phase transition in the frustrated XY and Heisenberg spin model Lattice: Stacked triangular lattice (3D) Spin: XY and Heisenberg Methods: Standard Monte Carlo method Wang-Landau method Model We consider the stacking of triangular lattices in the z direction The system size is N N N Lattice size : N = 12, 18, 24,, 120, 150 Hamiltonian V. T. Ngo and H. T. Diep, J. Appl. Phys. 103, 07C712 (2008). V. T. Ngo and H. T. Diep, Phys. Rev. E 78, (2008).
47 XY antiferromagnets Energy vs temperature for N = 84
48 XY antiferromagnets magnetization vs temperature for N = 84
49 XY antiferromagnets Energy histogram
50 Heisenberg antiferromagnets Energy and magnetization vs temperature for N = 120
51 Heisenberg antiferromagnets Energy histogram for N = 96 and 120
52 Heisenberg antiferromagnets Energy histogram for N = 150 Susceptibility vs temperature
53 Heisenberg antiferromagnets Susceptibility vs temperature Maximum of susceptibility versus N = 96, 108, 120, and 150 on a ln-ln scale
54 Crossover from first- to second-order transition in frustrated Ising antiferromagnetic films Study the nature of this phase transition in the case of a thin film as a function of the film thickness Lattice: Face center cubic thin films Spin: Ising antiferromagnetic Methods: Standard Monte Carlo method Wang-Landau method Green s function technique Model Consider a film of FCC lattice structure Hamiltonian X. T. Pham Phu, V. T. Ngo and H. T. Diep, Phys. Rev. E 79, (2009).
55 Ground state: The spin configuration depends on the interface interaction J s (a) ordering of type I for (b) ordering of type II for The energy of a surface spin for two configurations
56 Monte Carlo results for J s = J = -1 For the bulk case, L = N z = 12 Energy vs temperature Energy histogram with PBC (a) and without PBC (b) in z direction
57 For the N z = 4 Energy histogram for L = 20, 30 and 40
58 For the case N z = 2 and L = 120 Energy vs temperature Energy histogram
59 For the case N z = 2 and L =120 Specific heat vs temperature Susceptibilities of sublattices (a) 1 and (b) 3
60 The latent heat E as a function of film thickness
61 Critical exponents for N z = 2 Maximum sublattice susceptibility The maximum value of V 1
62 Fully frustrated simple cubic lattice Study the nature of the phase transition in the fully frustrated Lattice: simple cubic lattice Spin: Ising, XY and Heisenberg Methods: Standard Monte Carlo method Wang-Landau method Model Hamiltonian : for antiferromagnetic bond : for antiferromagnetic bond V. T. Ngo, D. Tien Hoang and H. T. Diep, Phys. Rev. E 82, (2010) V. T. Ngo, D. Tien Hoang and H. T. Diep, Mod. Phys. Lett. 25, 929 (2011) V. T. Ngo, D. Tien Hoang and H. T. Diep, J. Phys.: Condens. Matter 23, (2011)
63 Monte Carlo results for the case XY spin Energy and specific heat with N =24 Magnetization and susceptibility
64 Monte Carlo results for the case XY spin Energy histogram for N =24 Energy histogram The transition is clearly of first order at L = 36
65 Monte Carlo results for the case Heisenberg spin Energy and specific heat Magnetization and susceptibility
66 Monte Carlo results for the case Heisenberg spin Energy histogram The transition is clearly of first order at L = 70
67 Monte Carlo results for the case Ising spin Energy and specific heat
68 Monte Carlo results for the case Ising spin Energy histogram Maximum of the specific heat The transition is clearly of first order at L = 120 We conclude that a fully frustrated simple cubic lattice undergoes a first-order transition for Ising, XY and Heisenberg spin models
69
VSOP19, Quy Nhon 3-18/08/2013. Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam.
VSOP19, Quy Nhon 3-18/08/2013 Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam. Part III. Finite size effects and Reweighting methods III.1. Finite size effects III.2. Single histogram method III.3.
More informationVSOP19, Quy Nhon 3-18/08/2013. Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam.
VSOP19, Quy Nhon 3-18/08/2013 Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam. Part II. Monte Carlo Simulation Methods II.1. The spin models II.2. Boundary conditions II.3. Simple sampling Monte Carlo
More informationSurface 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 informationStacked Triangular XY Antiferromagnets: End of a Controversial Issue on the Phase Transition
Stacked Triangular XY Antiferromagnets: End of a Controversial Issue on the Phase Transition V. Thanh Ngo, H. T. Diep To cite this version: V. Thanh Ngo, H. T. Diep. Stacked Triangular XY Antiferromagnets:
More informationMonte 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 informationBlume-Emery-Griffiths Model In Thin Films of Stacked Triangular Lattices
Blume-Emery-Griffiths Model In Thin Films of Stacked Triangular Lattices Sahbi El Hog, H. T. Diep To cite this version: Sahbi El Hog, H. T. Diep. Blume-Emery-Griffiths Model In Thin Films of Stacked Triangular
More informationLecture V: Multicanonical Simulations.
Lecture V: Multicanonical Simulations. 1. Multicanonical Ensemble 2. How to get the Weights? 3. Example Runs (2d Ising and Potts models) 4. Re-Weighting to the Canonical Ensemble 5. Energy and Specific
More informationHexagonal-Close-Packed Lattice: Phase Transition and Spin Transport
Hexagonal-Close-Packed Lattice: Phase ransition and Spin ransport Danh-ai Hoang, Hung he Diep o cite this version: Danh-ai Hoang, Hung he Diep. Hexagonal-Close-Packed Lattice: Phase ransition and Spin
More informationWang-Landau Sampling of an Asymmetric Ising Model: A Study of the Critical Endpoint Behavior
Brazilian Journal of Physics, vol. 36, no. 3A, September, 26 635 Wang-andau Sampling of an Asymmetric Ising Model: A Study of the Critical Endpoint Behavior Shan-o sai a,b, Fugao Wang a,, and D.P. andau
More informationEvaluation of Wang-Landau Monte Carlo Simulations
2012 4th International Conference on Computer Modeling and Simulation (ICCMS 2012) IPCSIT vol.22 (2012) (2012) IACSIT Press, Singapore Evaluation of Wang-Landau Monte Carlo Simulations Seung-Yeon Kim School
More informationMonte Carlo Study of Planar Rotator Model with Weak Dzyaloshinsky Moriya Interaction
Commun. Theor. Phys. (Beijing, China) 46 (2006) pp. 663 667 c International Academic Publishers Vol. 46, No. 4, October 15, 2006 Monte Carlo Study of Planar Rotator Model with Weak Dzyaloshinsky Moriya
More informationMONTE 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 informationPhase transitions and finite-size scaling
Phase transitions and finite-size scaling Critical slowing down and cluster methods. Theory of phase transitions/ RNG Finite-size scaling Detailed treatment: Lectures on Phase Transitions and the Renormalization
More informationNumerical 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 informationA New Method to Determine First-Order Transition Points from Finite-Size Data
A New Method to Determine First-Order Transition Points from Finite-Size Data Christian Borgs and Wolfhard Janke Institut für Theoretische Physik Freie Universität Berlin Arnimallee 14, 1000 Berlin 33,
More informationWORLD SCIENTIFIC (2014)
WORLD SCIENTIFIC (2014) LIST OF PROBLEMS Chapter 1: Magnetism of Free Electrons and Atoms 1. Orbital and spin moments of an electron: Using the theory of angular momentum, calculate the orbital
More informationCluster Algorithms to Reduce Critical Slowing Down
Cluster Algorithms to Reduce Critical Slowing Down Monte Carlo simulations close to a phase transition are affected by critical slowing down. In the 2-D Ising system, the correlation length ξ becomes very
More informationGPU-based computation of the Monte Carlo simulation of classical spin systems
Perspectives of GPU Computing in Physics and Astrophysics, Sapienza University of Rome, Rome, Italy, September 15-17, 2014 GPU-based computation of the Monte Carlo simulation of classical spin systems
More informationarxiv: v1 [cond-mat.dis-nn] 12 Nov 2014
Representation for the Pyrochlore Lattice arxiv:1411.3050v1 [cond-mat.dis-nn] 12 Nov 2014 André Luis Passos a, Douglas F. de Albuquerque b, João Batista Santos Filho c Abstract a DFI, CCET, Universidade
More informationPopulation 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 informationGraphical 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 informationMonte 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 informationPotts And XY, Together At Last
Potts And XY, Together At Last Daniel Kolodrubetz Massachusetts Institute of Technology, Center for Theoretical Physics (Dated: May 16, 212) We investigate the behavior of an XY model coupled multiplicatively
More information3.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 informationMonte Carlo Simulation of Spins
Monte Carlo Simulation of Spins Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Department of Computer Science Department of Physics & Astronomy Department of Chemical Engineering &
More informationJ 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 informationInvaded 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 informationPotts percolation Gauss model of a solid
Cleveland State University From the SelectedWorks of Miron Kaufman January, 2008 Potts percolation Gauss model of a solid Miron Kaufman H T Diep Available at: https://works.bepress.com/miron_kaufman/21/
More informationStatistical description of magnetic domains in the Ising model
arxiv:0804.3522v1 [cond-mat.stat-mech] 22 Apr 2008 Statistical description of magnetic domains in the Ising model K. Lukierska-Walasek Institute of Physics University of Zielona Góra ul. Z. Szafrana 4a,
More informationLecture 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 informationWang-Landau sampling for Quantum Monte Carlo. Stefan Wessel Institut für Theoretische Physik III Universität Stuttgart
Wang-Landau sampling for Quantum Monte Carlo Stefan Wessel Institut für Theoretische Physik III Universität Stuttgart Overview Classical Monte Carlo First order phase transitions Classical Wang-Landau
More informationClassical Monte Carlo Simulations
Classical Monte Carlo Simulations Hyejin Ju April 17, 2012 1 Introduction Why do we need numerics? One of the main goals of condensed matter is to compute expectation values O = 1 Z Tr{O e βĥ} (1) and
More informationStochastic series expansion (SSE) and ground-state projection
Institute of Physics, Chinese Academy of Sciences, Beijing, October 31, 2014 Stochastic series expansion (SSE) and ground-state projection Anders W Sandvik, Boston University Review article on quantum
More informationarxiv: v1 [cond-mat.stat-mech] 22 Sep 2009
Phase diagram and critical behavior of the square-lattice Ising model with competing nearest- and next-nearest-neighbor interactions Junqi Yin and D. P. Landau Center for Simulational Physics, University
More informationInterface tension of the 3d 4-state Potts model using the Wang-Landau algorithm
Interface tension of the 3d 4-state Potts model using the Wang-Landau algorithm CP3-Origins and the Danish Institute for Advanced Study DIAS, University of Southern Denmark, Campusvej 55, DK-5230 Odense
More informationCONTINUOUS- AND FIRST-ORDER PHASE TRANSITIONS IN ISING ANTIFERROMAGNETS WITH NEXT-NEAREST- NEIGHBOUR INTERACTIONS
Continuous- Rev.Adv.Mater.Sci. and first-order 14(2007) phase 1-10 transitions in ising antiferromagnets with next-nearest-... 1 CONTINUOUS- AND FIRST-ORDER PHASE TRANSITIONS IN ISING ANTIFERROMAGNETS
More informationarxiv:cond-mat/ v1 22 Sep 1998
Scaling properties of the cluster distribution of a critical nonequilibrium model Marta Chaves and Maria Augusta Santos Departamento de Física and Centro de Física do Porto, Faculdade de Ciências, Universidade
More informationREVIEW: 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 informationBrazilian Journal of Physics, vol. 26, no. 4, december, P. M.C.deOliveira,T.J.P.Penna. Instituto de Fsica, Universidade Federal Fluminense
Brazilian Journal of Physics, vol. 26, no. 4, december, 1996 677 Broad Histogram Method P. M.C.deOliveira,T.J.P.Penna Instituto de Fsica, Universidade Federal Fluminense Av. Litor^anea s/n, Boa Viagem,
More informationUncovering the Secrets of Unusual Phase Diagrams: Applications of Two-Dimensional Wang-Landau Sampling
6 Brazilian Journal of Physics, vol. 38, no. 1, March, 28 Uncovering the Secrets of Unusual Phase Diagrams: Applications of wo-dimensional Wang-Landau Sampling Shan-Ho sai a,b, Fugao Wang a,, and D. P.
More informationCritical 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 informationQuantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University
PY502, Computational Physics, December 12, 2017 Quantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University Advancing Research in Basic Science and Mathematics Example:
More informationPhase transitions in the Potts spin-glass model
PHYSICAL REVIEW E VOLUME 58, NUMBER 3 SEPTEMBER 1998 Phase transitions in the Potts spin-glass model Giancarlo Franzese 1 and Antonio Coniglio 1,2 1 Dipartimento di Scienze Fisiche, Università di Napoli,
More informationPhase Transitions of Random Binary Magnetic Square Lattice Ising Systems
I. Q. Sikakana Department of Physics and Non-Destructive Testing, Vaal University of Technology, Vanderbijlpark, 1900, South Africa e-mail: ike@vut.ac.za Abstract Binary magnetic square lattice Ising system
More informationThe 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 informationAdvanced Monte Carlo Methods Problems
Advanced Monte Carlo Methods Problems September-November, 2012 Contents 1 Integration with the Monte Carlo method 2 1.1 Non-uniform random numbers.......................... 2 1.2 Gaussian RNG..................................
More informationMetropolis 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 informationOvercoming the slowing down of flat-histogram Monte Carlo simulations: Cluster updates and optimized broad-histogram ensembles
PHYSICAL REVIEW E 72, 046704 2005 Overcoming the slowing down of flat-histogram Monte Carlo simulations: Cluster updates and optimized broad-histogram ensembles Yong Wu, 1 Mathias Körner, 2 Louis Colonna-Romano,
More informationGeneralized Ensembles: Multicanonical Simulations
Generalized Ensembles: Multicanonical Simulations 1. Multicanonical Ensemble 2. How to get the Weights? 3. Example Runs and Re-Weighting to the Canonical Ensemble 4. Energy and Specific Heat Calculation
More informationTopological defects and its role in the phase transition of a dense defect system
Topological defects and its role in the phase transition of a dense defect system Suman Sinha * and Soumen Kumar Roy Depatrment of Physics, Jadavpur University Kolkata- 70003, India Abstract Monte Carlo
More informationarxiv: 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 informationLECTURE 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 informationPartition function zeros and finite size scaling for polymer adsorption
CompPhys13: New Developments in Computational Physics, Leipzig, Nov. 28-30, 2013 Partition function zeros and finite size scaling for polymer adsorption Mark P. Taylor 1 and Jutta Luettmer-Strathmann 2
More informationWang-Landau algorithm: A theoretical analysis of the saturation of the error
THE JOURNAL OF CHEMICAL PHYSICS 127, 184105 2007 Wang-Landau algorithm: A theoretical analysis of the saturation of the error R. E. Belardinelli a and V. D. Pereyra b Departamento de Física, Laboratorio
More informationProgress 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 informationNumerical Study of the Antiferromagnetic Ising Model in Hyperdimensions
Adv. Studies Theor. Phys., Vol. 3, 2009, no. 12, 481-488 Numerical Study of the Antiferromagnetic Ising Model in Hyperdimensions N. Olivi-Tran GES, UMR-CNRS 5650, Universite Montpellier II cc 074, place
More informationThe 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 informationA Monte Carlo Study of the Specific Heat of Diluted Antiferromagnets
A Monte Carlo Study of the Specific Heat of Diluted Antiferromagnets M. Staats and K. D. Usadel email: michael@thp.uni-duisburg.de Theoretische Physik and SFB 166 Gerhard-Mercator-Universität Gesamthochschule
More informationPhysics 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 informationMonte 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 informationThe critical behaviour of the long-range Potts chain from the largest cluster probability distribution
Physica A 314 (2002) 448 453 www.elsevier.com/locate/physa The critical behaviour of the long-range Potts chain from the largest cluster probability distribution Katarina Uzelac a;, Zvonko Glumac b a Institute
More informationFinite-size analysis via the critical energy -subspace method in the Ising models
Materials Science-Poland, Vol. 23, No. 4, 25 Finite-size analysis via the critical energy -subspace method in the Ising models A. C. MAAKIS *, I. A. HADJIAGAPIOU, S. S. MARTINOS, N. G. FYTAS Faculty of
More informationPHYSICAL REVIEW LETTERS
PHYSICAL REVIEW LETTERS VOLUME 76 4 MARCH 1996 NUMBER 10 Finite-Size Scaling and Universality above the Upper Critical Dimensionality Erik Luijten* and Henk W. J. Blöte Faculty of Applied Physics, Delft
More informationMassively parallel Monte Carlo simulation of a possible topological phase transition in two-dimensional frustrated spin systems
Massively parallel Monte Carlo simulation of a possible topological phase transition in two-dimensional frustrated spin systems Tsuyoshi OKUBO Institute for Solid State Physics, University of Tokyo Kashiwa-no-ha,
More informationICCP 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 informationMulticanonical 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 informationMagnetism at finite temperature: molecular field, phase transitions
Magnetism at finite temperature: molecular field, phase transitions -The Heisenberg model in molecular field approximation: ferro, antiferromagnetism. Ordering temperature; thermodynamics - Mean field
More informationQuantum and classical annealing in spin glasses and quantum computing. Anders W Sandvik, Boston University
NATIONAL TAIWAN UNIVERSITY, COLLOQUIUM, MARCH 10, 2015 Quantum and classical annealing in spin glasses and quantum computing Anders W Sandvik, Boston University Cheng-Wei Liu (BU) Anatoli Polkovnikov (BU)
More information3.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 informationARTICLES. Monte Carlo study of a compressible Ising antiferromagnet on a triangular lattice
PHYSICAL REVIEW B VOLUME 53, NUMBER 18 ARTICLES 1 MAY 1996-II Monte Carlo study of a compressible Ising antiferromagnet on a triangular lattice Lei Gu and Bulbul Chakraborty The Martin Fisher School of
More informationParallel Tempering Algorithm in Monte Carlo Simulation
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) Metropolis Algorithm on Ising Model
More informationDecimation Technique on Sierpinski Gasket in External Magnetic Field
Egypt.. Solids, Vol. (), No. (), (009 ) 5 Decimation Technique on Sierpinski Gasket in External Magnetic Field Khalid Bannora, G. Ismail and M. Abu Zeid ) Mathematics Department, Faculty of Science, Zagazig
More informationMonte Carlo and spin dynamics simulation of the fully frustrated anisotropic two-dimensional Heisenberg model
Journal of Magnetism and Magnetic Materials 6 () 4 Monte Carlo and spin dynamics simulation of the fully frustrated anisotropic two-dimensional Heisenberg model A.B. Lima*, B.V. Costa Departamento de F!ısica,
More informationAGuideto Monte Carlo Simulations in Statistical Physics
AGuideto Monte Carlo Simulations in Statistical Physics Second Edition David P. Landau Center for Simulational Physics, The University of Georgia Kurt Binder Institut für Physik, Johannes-Gutenberg-Universität
More informationOptimized 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 information8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization
8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization This problem set is partly intended to introduce the transfer matrix method, which is used
More informationPY 502, Computational Physics, Fall 2017 Monte Carlo simulations in classical statistical physics
PY 502, Computational Physics, Fall 2017 Monte Carlo simulations in classical statistical physics Anders W. Sandvik, Department of Physics, Boston University 1 Introduction Monte Carlo simulation is a
More informationMonte Carlo and cold gases. Lode Pollet.
Monte Carlo and cold gases Lode Pollet lpollet@physics.harvard.edu 1 Outline Classical Monte Carlo The Monte Carlo trick Markov chains Metropolis algorithm Ising model critical slowing down Quantum Monte
More information(1) Consider the ferromagnetic XY model, with
PHYSICS 10A : STATISTICAL PHYSICS HW ASSIGNMENT #7 (1 Consider the ferromagnetic XY model, with Ĥ = i
More informationMarkov Chain Monte Carlo Simulations and Their Statistical Analysis An Overview
Markov Chain Monte Carlo Simulations and Their Statistical Analysis An Overview Bernd Berg FSU, August 30, 2005 Content 1. Statistics as needed 2. Markov Chain Monte Carlo (MC) 3. Statistical Analysis
More informationUniversality class of triad dynamics on a triangular lattice
Universality class of triad dynamics on a triangular lattice Filippo Radicchi,* Daniele Vilone, and Hildegard Meyer-Ortmanns School of Engineering and Science, International University Bremen, P. O. Box
More informationRensselaer Polytechnic Institute, 110 8th Street, Troy, NY , USA. Florida State University, Tallahassee, Florida , USA
5 Dynamic Phase Diagram for a Periodically Driven Kinetic Square-lattice Ising Ferromagnet: Finite-size Scaling Evidence for the Absence of a Tri-critical Point G. Korniss 1, P.A. Rikvold 2,3, and M.A.
More informationarxiv: v2 [cond-mat.stat-mech] 10 May 2017
Macroscopically constrained Wang-Landau method for systems with multiple order parameters and its application to drawing arxiv:1703.03019v2 [cond-mat.stat-mech] 10 May 2017 complex phase diagrams C. H.
More informationComplex Systems Methods 9. Critical Phenomena: The Renormalization Group
Complex Systems Methods 9. Critical Phenomena: The Renormalization Group Eckehard Olbrich e.olbrich@gmx.de http://personal-homepages.mis.mpg.de/olbrich/complex systems.html Potsdam WS 2007/08 Olbrich (Leipzig)
More informationPhase diagram of the Ising square lattice with competing interactions
Eur. Phys. J. B6 (8) 7 arxiv:8.98 Phase diagram of the Ising square lattice with competing interactions A. Kalz, A. Honecker, S. Fuchs and T. Pruschke Institut für Theoretische Physik, Georg-August-Universität
More informationPhase Transitions in Condensed Matter Spontaneous Symmetry Breaking and Universality. Hans-Henning Klauss. Institut für Festkörperphysik TU Dresden
Phase Transitions in Condensed Matter Spontaneous Symmetry Breaking and Universality Hans-Henning Klauss Institut für Festkörperphysik TU Dresden 1 References [1] Stephen Blundell, Magnetism in Condensed
More informationMulticanonical 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 informationPhysics 127b: Statistical Mechanics. Landau Theory of Second Order Phase Transitions. Order Parameter
Physics 127b: Statistical Mechanics Landau Theory of Second Order Phase Transitions Order Parameter Second order phase transitions occur when a new state of reduced symmetry develops continuously from
More informationAvailable online at ScienceDirect. Physics Procedia 57 (2014 ) 82 86
Available online at www.sciencedirect.com ScienceDirect Physics Procedia 57 (2014 ) 82 86 27th Annual CSP Workshops on Recent Developments in Computer Simulation Studies in Condensed Matter Physics, CSP
More informationarxiv: v1 [cond-mat.stat-mech] 9 Feb 2012
Magnetic properties and critical behavior of disordered Fe 1 x Ru x alloys: a Monte Carlo approach I. J. L. Diaz 1, and N. S. Branco 2, arxiv:1202.2104v1 [cond-mat.stat-mech] 9 Feb 2012 1 Universidade
More informationComputer simulations as concrete models for student reasoning
Computer simulations as concrete models for student reasoning Jan Tobochnik Department of Physics Kalamazoo College Kalamazoo MI 49006 In many thermal physics courses, students become preoccupied with
More information4. 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 informationCover Page. The handle holds various files of this Leiden University dissertation.
Cover Page The handle http://hdl.handle.net/1887/49403 holds various files of this Leiden University dissertation. Author: Keesman, R. Title: Topological phases and phase transitions in magnets and ice
More informationQuantum Monte Carlo Simulations in the Valence Bond Basis. Anders Sandvik, Boston University
Quantum Monte Carlo Simulations in the Valence Bond Basis Anders Sandvik, Boston University Outline The valence bond basis for S=1/2 spins Projector QMC in the valence bond basis Heisenberg model with
More informationSpin liquids in frustrated magnets
May 20, 2010 Contents 1 Frustration 2 3 4 Exotic excitations 5 Frustration The presence of competing forces that cannot be simultaneously satisfied. Heisenberg-Hamiltonian H = 1 J ij S i S j 2 ij The ground
More informationABSTRACT. We measure the efficiency of the Metropolis, Swendsen-Wang, and Wolff algorithms in
ABSTRACT KYIMBA, ELOI-ALAIN KYALONDAWA. Comparison of Monte Carlo Metropolis, Swendsen-Wang, and Wolff Algorithms in the Critical Region for the 2-Dimensional Ising Model. (Under the direction of Dean
More informationAvailable online at ScienceDirect. Physics Procedia 53 (2014 ) Andressa A. Bertolazzo and Marcia C.
Available online at www.sciencedirect.com ScienceDirect Physics Procedia 53 (2014 ) 7 15 Density and Diffusion Anomalies in a repulsive lattice gas Andressa A. Bertolazzo and Marcia C. Barbosa Instituto
More informationSolitonic elliptical solutions in the classical XY model
Solitonic elliptical solutions in the classical XY model Rodrigo Ferrer, José Rogan, Sergio Davis, Gonzalo Gutiérrez Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago,
More informationEquilibrium, 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 information8.3.2 The finite size scaling method
232 Chapter 8: Analysing Monte Carlo data In general we don t know this value, which makes it difficult to perform the fit. It is possible to guess T c and then vary the guess to make the line in Figure
More information