36 lines
2.2 KiB
TeX
36 lines
2.2 KiB
TeX
\documentclass[../ising_model.tex]{subfiles}
|
|
|
|
\begin{document}
|
|
\section{Conclusion}\label{sec:conclusion}
|
|
We studied the ferromagnetic behavior of the Ising model, and observed a phase
|
|
transition at temperatures close to the critical temperature. We used a Markov
|
|
chain Monte Carlo sampling method to generate spin configurations, while utilizing
|
|
parallel programming techniques. When increasing number of processes, and number
|
|
of threads we observed a speed-up in runtime.
|
|
|
|
We initialized the lattices using both ordered and unordered spin configuration,
|
|
and started sampling after the system reached an equilibrium. The number of Monte
|
|
Carlo cycles necessary to reach a system equilibrium, referred to as burn-in time,
|
|
was estimated to be 3000 cycles. We found that excluding the samples generated
|
|
during the burn-in time, improves the estimated expectation value of energy and
|
|
magnetization, in addition to the heat capacity and susceptibility, when samples
|
|
are sparse. However, when we increase number of samples, excluding the burn-in
|
|
samples does not affect the estimated values.
|
|
|
|
Continuing, we used the generated samples to compute energy per spin $\langle \epsilon \rangle$,
|
|
magnetization per spin $\langle |m| \rangle$, heat capacity $C_{V}$, and $\chi$.
|
|
In addition, we estimated the probability distribution for temperature $T_{1} = 1.0 J / k_{B}$,
|
|
and $T_{2} = 2.4 J / k_{B}$. We found that for $T_{1}$ the expected mean energy
|
|
per spin is $\langle \epsilon \rangle \approx -1.9969 J$, with a variation $\text{Var} (\epsilon) = 0.0001$.
|
|
And for $T_{2}$, the mean energy per spin is $\langle \epsilon \rangle \approx -1.2370 J$,
|
|
with a variation $\text{Var} (\epsilon) = 0.0203$.
|
|
|
|
We estimated the expected energy and magnetization per spin, in addition to the
|
|
heat capacity and susceptibility for lattices of size $L = {20, 40, 60, 80, 100}$.
|
|
We observed a phase transition in the temperature range $T \in [2.1, 2.4] J / k_{B}$.
|
|
Using the values from the finite lattices, we approximated the critical temperature
|
|
of a lattice of infinite size. Using linear regression, we numerically estimated
|
|
$T_{c}(L = \infty) \approx 2.2695 J/k_{B}$ which is close to the analytical solution
|
|
$T_{C}(L = \infty) \approx 2.269 J/k_{B}$ Lars Onsager found in 1944.
|
|
\end{document}
|