Project 4: Difference between revisions
(40 intermediate revisions by the same user not shown) | |||
Line 1: | Line 1: | ||
'''First and second order phase | '''First and second order phase transitions in the Ising model of ferromagnetism''' | ||
==Introduction== | ==Introduction== | ||
Write a program which simulates via Metropolis Monte Carlo algorithm the 2D Ising model on a <math> | Write a program which simulates via Metropolis Monte Carlo algorithm the 2D Ising model on a <math> 128 \times 128 </math> square lattice with ''periodic'' boundary conditions. The model is assumed to describe a system placed in an applied external magnetic field <math> H </math>. The energy of a particular state of this system is then given by: | ||
Line 13: | Line 13: | ||
==Part I for both PHYS460 and PHYS660 students: First order phase transition in the ferromagnetic Ising model== | ==Part I for both PHYS460 and PHYS660 students: First order phase transition in the ferromagnetic Ising model== | ||
Consider the ferromagnetic Ising model with <math> J | Consider the ferromagnetic Ising model with <math> J =1 </math>. Start from a case in which temperature is below the critical one <math> T_c =2.269 </math>, such as <math> T=1.0 </math>, and the field is large and pointing down, e.g., <math> H = -5 </math>. The selected initial state of the Ising model should be the one in which all spins <math> S_i = -1 </math> are parallel to H, so that magnetization per spin is <math> M = -1 </math>. | ||
Sweep the magnetic field up to a large positive value <math> H=5 </math> and then back down to your starting value <math> H = -5 </math>. Do this in steps of <math> \delta H = 0.5 </math>. After each change in the magnetic field <math> H </math>, allow sufficient Monte Carlo trials | Sweep the magnetic field up to a large positive value <math> H=5 </math> and then back down to your starting value <math> H = -5 </math>. Do this in steps of <math> \delta H = 0.5 </math>. After each change in the magnetic field <math> H </math>, allow sufficient Monte Carlo trials for the system to reach the equilibrium and thermalize. | ||
for the system to reach the equilibrium and thermalize. | |||
For each value of <math> H </math>, | For each value of <math> H </math>, the mean magnetization per spin | ||
<math> \langle M \rangle = \frac{1}{N_{MC}} \sum_{\alpha=1}^{N_{MC}} M_\alpha </math> | |||
The sum in this expression run are over <math> N_{MC} </math> microstates generated by Monte Carlo simulation. | |||
'''(a)''' Plot <math> M </math> vs. <math> H </math>. You should see a discontinuous change in as H is increased, and a second discontinuous change in <math> M </math> as <math> H </math> is decreased back to its initial value. These discontinuous changes indicate a ''first order phase transition''. Also note that the values of <math> H </math> at which discontinuities occur are not identical. Such behavior is called hysteresis. Perform at least one thousand Monte Carlo time steps, where '''one time step is one complete pass through the lattice''', at each value of the magnetic field. | |||
'''(b)''' Repeat the calculations for a temperature <math> T=4.0 </math> above the critical point <math> T_c </math>. Is there a first order phase transition? Is there hysteresis? Provide a physical explanation of any difference in behavior for the two temperatures. | |||
==Part II for both PHYS460 and PHYS660 students: Second order phase transition in the ferromagnetic Ising model== | ==Part II for both PHYS460 and PHYS660 students: Second order phase transition in the ferromagnetic Ising model== | ||
Investigate the behavior of the specific heat as the temperature | '''(a)''' Plot magnetization per spin and energy per spin for each sample obtained after one Monte Carlo sweep through the lattice as a function of Monte Carlo time step at three temperatures <math> T<T_c</math>, <math> T=T_c</math>, and <math> T>T_c</math> in the absence of magnetic field <math> H=0 </math>. These plots should reveal how many Monte Carlo steps is necessary to thermalize (i.e., bring it into the thermodynamic equilibrium with external heat bath which is assumed to set the temperature <math> T </math>) the Ising ferromagnet. Thermodynamic average values should be computed ''only'' for the samples generated after the thermalization has been achieved. | ||
'''(b)''' Investigate the behavior of the specific heat | |||
<math> C_V = \frac{1}{k_B T^2} ( \langle E^2 \rangle - \langle E \rangle^2) </math> | |||
and magnetic susceptibility | |||
<math> \chi = \frac{1}{k_B T} ( \langle M^2 \rangle - \langle M \rangle^2) </math> | |||
as the function of increasing temperature in the absence of external magnetic field <math> H=0 </math>. In this case we expect a second order or continuous paramagnet-ferromanget phase transition at the critical temperature <math> T_c </math> around which these two quantities diverge in the thermodynamic limit. Since your calculations is for finite size system, repeat computation of these two quantities while increasing the lattice size (e.g., <math> 16 \times 16 </math>, <math> 32 \times 32 </math>, <math> 64 \times 64 </math>, ...) in order to observe how the peak in <math> C_V </math> and <math> \chi </math> around <math> T_c </math> increases while its center shifts because toward <math> T_c </math>. The two quantities should be divided by the lattice size to obtain their values "per spin". | |||
The connection between these two response functions and fluctuations | |||
<math> (\Delta X)^2 = \langle X^2 \rangle - \langle X \rangle^2 </math> | |||
where | |||
<math> \langle X^2 \rangle = \frac{1}{N_{MC}} \sum_{\alpha=1}^{N_{MC}} X_\alpha^2 </math> | |||
<math> \langle X \rangle^2 = \left(\frac{1}{N_{MC}} \sum_{\alpha=1}^{N_{MC}} X_\alpha \right)^2 </math> | |||
utilizes the so-called [http://iopscience.iop.org/0034-4885/29/1/306 fluctuation-dissipation theorem] of statistical physics. | |||
==Part III for PHYS660 students only: Correlation length in the ferromagnetic Ising model== | ==Part III for PHYS660 students only: Correlation length in the ferromagnetic Ising model== | ||
The correlation lenght can be obtained from the r-dependence of the spin-spin correlation function | |||
<math> G(r) = \langle S_i S_j \rangle - \langle S_i \rangle \langle S_j \rangle </math> | |||
where <math> \langle S_i \rangle \langle S_j \rangle = M^2</math> since the system is translationally invariant. Here <math> r=|i-j| </math> is the distance between the sites | |||
<math> i </math> and <math> j </math>. The average here <math> \langle \ldots \rangle </math> is | |||
defined over all sites for a given configuration and over many configurations. Because the spins are not correlated for large <math> r </math>, <math> G(r) \rightarrow 0 </math> in this limit. Assume that <math> G(r) \sim \exp(-r/\xi) </math> for sufficiently large <math> r </math> and estimate as a function of temperature <math> T </math> in zero external field <math> H </math>. | |||
How does your estimate for <math> \xi </math> compare with the size of the domains of spins | |||
pointing in the same direction? We expect that <math> \xi </math> should be of the order of a lattice spacing at <math> T \gg T_c </math> and approaching the system size (or diverging | |||
for infinite systems in the thermodynamic limit) as <math> T \rightarrow T_c </math>. | |||
==Part IV for PHYS660 students only: Second order phase transition in the antiferromagnetic Ising model== | ==Part IV for PHYS660 students only: Second order phase transition in the antiferromagnetic Ising model== | ||
Compute temperature dependence of <math> C_V </math>, <math> \chi </math>, and <math> M </math> in the absence of magnetic field <math> H =0 </math> for the antiferromagnetic Ising model with <math> J=-1 </math>. The initial configuration can be chosen as all spins up <math> S_i = +1 </math> on a <math> 128 \times 128 </math> lattice. What configuration of spins corresponds to the thermodynamic state below the Neel temperature <math> T_N = 2.269 </math>? At the so-called Neel temperature this type of materials exhibit continous antiferromagnet-paramagnet phase transition. | |||
==References== | |||
*I. Vattulainen, T. Ala-Nissila, and K. Kankaala, ''Physical tests for random numbers in simulations'', Phys. Rev. Lett. '''73''', 2513 (1994). [http://prl.aps.org/abstract/PRL/v73/i19/p2513_1 [PDF]] | |||
*D. P. Landau, S.-H. Tsai, and M. Exler, ''A new approach to Monte Carlo simulations in statistical physics: Wang-Landau sampling'', Am. J. Phys. '''72''', 1294 (2004). [http://ajp.aapt.org/resource/1/ajpias/v72/i10/p1294_s1 [PDF]] |
Latest revision as of 23:23, 22 April 2014
First and second order phase transitions in the Ising model of ferromagnetism
Introduction
Write a program which simulates via Metropolis Monte Carlo algorithm the 2D Ising model on a square lattice with periodic boundary conditions. The model is assumed to describe a system placed in an applied external magnetic field . The energy of a particular state of this system is then given by:
where means that the sum runs over all pairs of nearest neighbor spins (attached to the lattice sites), and is the exchange coupling. Positive interaction between the spins leads to ferromagnetism, will negative gives rise to antiferromagnetism in the Ising model. Choose units so that and the Boltzmann constant is . The free parameters of the model are then just the temperature and the applied field .
Part I for both PHYS460 and PHYS660 students: First order phase transition in the ferromagnetic Ising model
Consider the ferromagnetic Ising model with . Start from a case in which temperature is below the critical one , such as , and the field is large and pointing down, e.g., . The selected initial state of the Ising model should be the one in which all spins are parallel to H, so that magnetization per spin is .
Sweep the magnetic field up to a large positive value and then back down to your starting value . Do this in steps of . After each change in the magnetic field , allow sufficient Monte Carlo trials for the system to reach the equilibrium and thermalize.
For each value of , the mean magnetization per spin
The sum in this expression run are over microstates generated by Monte Carlo simulation.
(a) Plot vs. . You should see a discontinuous change in as H is increased, and a second discontinuous change in as is decreased back to its initial value. These discontinuous changes indicate a first order phase transition. Also note that the values of at which discontinuities occur are not identical. Such behavior is called hysteresis. Perform at least one thousand Monte Carlo time steps, where one time step is one complete pass through the lattice, at each value of the magnetic field.
(b) Repeat the calculations for a temperature above the critical point . Is there a first order phase transition? Is there hysteresis? Provide a physical explanation of any difference in behavior for the two temperatures.
Part II for both PHYS460 and PHYS660 students: Second order phase transition in the ferromagnetic Ising model
(a) Plot magnetization per spin and energy per spin for each sample obtained after one Monte Carlo sweep through the lattice as a function of Monte Carlo time step at three temperatures , , and in the absence of magnetic field . These plots should reveal how many Monte Carlo steps is necessary to thermalize (i.e., bring it into the thermodynamic equilibrium with external heat bath which is assumed to set the temperature ) the Ising ferromagnet. Thermodynamic average values should be computed only for the samples generated after the thermalization has been achieved.
(b) Investigate the behavior of the specific heat
and magnetic susceptibility
as the function of increasing temperature in the absence of external magnetic field . In this case we expect a second order or continuous paramagnet-ferromanget phase transition at the critical temperature around which these two quantities diverge in the thermodynamic limit. Since your calculations is for finite size system, repeat computation of these two quantities while increasing the lattice size (e.g., , , , ...) in order to observe how the peak in and around increases while its center shifts because toward . The two quantities should be divided by the lattice size to obtain their values "per spin".
The connection between these two response functions and fluctuations
where
utilizes the so-called fluctuation-dissipation theorem of statistical physics.
Part III for PHYS660 students only: Correlation length in the ferromagnetic Ising model
The correlation lenght can be obtained from the r-dependence of the spin-spin correlation function
where since the system is translationally invariant. Here is the distance between the sites and . The average here is defined over all sites for a given configuration and over many configurations. Because the spins are not correlated for large , in this limit. Assume that for sufficiently large and estimate as a function of temperature in zero external field .
How does your estimate for compare with the size of the domains of spins pointing in the same direction? We expect that should be of the order of a lattice spacing at and approaching the system size (or diverging for infinite systems in the thermodynamic limit) as .
Part IV for PHYS660 students only: Second order phase transition in the antiferromagnetic Ising model
Compute temperature dependence of , , and in the absence of magnetic field for the antiferromagnetic Ising model with . The initial configuration can be chosen as all spins up on a lattice. What configuration of spins corresponds to the thermodynamic state below the Neel temperature ? At the so-called Neel temperature this type of materials exhibit continous antiferromagnet-paramagnet phase transition.
References
- I. Vattulainen, T. Ala-Nissila, and K. Kankaala, Physical tests for random numbers in simulations, Phys. Rev. Lett. 73, 2513 (1994). [PDF]
- D. P. Landau, S.-H. Tsai, and M. Exler, A new approach to Monte Carlo simulations in statistical physics: Wang-Landau sampling, Am. J. Phys. 72, 1294 (2004). [PDF]