Calculating Magnetization for 2x2 Ising Model Lattice

LagrangeEuler
Messages
711
Reaction score
22

Homework Statement


Calculate magnetisation for partition function
##Z=12+4\cosh (8\beta J)## for Ising model 2x2 lattice.

Homework Equations


F=-k_BTln Z
M(H,T)=-\frac{\partial}{\partial H}(\frac{F}{k_BT})

The Attempt at a Solution


For me it looks that magnetisation is zero. By just doing the derivatives. But solution of problem is ##\frac{16+8e^{\beta J}}{Z}##. If I combine relation from relevant equation it looks like for me that solution is zero. Where I making the mistake.
 
Physics news on Phys.org
Where does the partition function come from?

There is no H in the partition function, so the system does not react to any external field, just as you calculated.

What are the units of J?
 
##J## is constant. Well you have 2x2 lattice. You calculate
Z=\sum e^{-\beta H}
 
You are confusing notation here. The derivative to find the magnetization should be with respect to h, the external magnetic field. H is the Hamiltonian used to find the partition function. If there is no term dependent on h in your partition function expression, then you have used an incorrect Hamiltonian. Assuming a constant interaction J and constant external field h, the 2D Ising Model has a Hamiltonian of the form

H=-J \sum_{i,j} \sigma_i \sigma_j - h\sum_j \sigma_j

The first sum is over all nearest neighbors.As written your partition function trivially has no magnetization.
 
Last edited:
I think there are a few sources of confusion here...

It seems that the problem presents a certain partition function as a given. Presumably it's some kind of approximation for the 2-d Ising model Z.

Also it looks like J is some parameter that's proportional to the applied magnetic field. So dZ/dH won't be zero.
 
Oxvillian said:
I think there are a few sources of confusion here...

It seems that the problem presents a certain partition function as a given. Presumably it's some kind of approximation for the 2-d Ising model Z.

Also it looks like J is some parameter that's proportional to the applied magnetic field. So dZ/dH won't be zero.

J is typically the spin-spin coupling in Ising models, and I suspect that is the case here. The partition function given simply doesn't depend on the magnetic field.

To the OP: why should one expect the magnetization to be non-zero in the absence of a magnetic field? You have a 2x2 Ising model at finite temperature. Furthermore, I assume that the spin-spin interactions are nearest neighbor only, is that correct?

If so, consider this question: what is the difference between your 2x2 Ising model and a 4-spin 1d Ising model with periodic boundary conditions (and no external field)?
 
##|\Psi|^2=\frac{1}{\sqrt{\pi b^2}}\exp(\frac{-(x-x_0)^2}{b^2}).## ##\braket{x}=\frac{1}{\sqrt{\pi b^2}}\int_{-\infty}^{\infty}dx\,x\exp(-\frac{(x-x_0)^2}{b^2}).## ##y=x-x_0 \quad x=y+x_0 \quad dy=dx.## The boundaries remain infinite, I believe. ##\frac{1}{\sqrt{\pi b^2}}\int_{-\infty}^{\infty}dy(y+x_0)\exp(\frac{-y^2}{b^2}).## ##\frac{2}{\sqrt{\pi b^2}}\int_0^{\infty}dy\,y\exp(\frac{-y^2}{b^2})+\frac{2x_0}{\sqrt{\pi b^2}}\int_0^{\infty}dy\,\exp(-\frac{y^2}{b^2}).## I then resolved the two...
Hello everyone, I’m considering a point charge q that oscillates harmonically about the origin along the z-axis, e.g. $$z_{q}(t)= A\sin(wt)$$ In a strongly simplified / quasi-instantaneous approximation I ignore retardation and take the electric field at the position ##r=(x,y,z)## simply to be the “Coulomb field at the charge’s instantaneous position”: $$E(r,t)=\frac{q}{4\pi\varepsilon_{0}}\frac{r-r_{q}(t)}{||r-r_{q}(t)||^{3}}$$ with $$r_{q}(t)=(0,0,z_{q}(t))$$ (I’m aware this isn’t...
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
Back
Top