Partition Function for Spin-1 One Dimensional Ising Model

pauladancer
Messages
26
Reaction score
0
Homework Statement
Hi everyone,
I'm writing a paper for my statistical mechanics course and require the partition function for the spin-1 Ising model. I've searched for a solution, but can't find one anywhere. I'm hoping to get some help!
Relevant Equations
See below
$$H=-J\sum_{i=1}^{N-1}\sigma_i\sigma_{i+1}$$ There is no external magnetic field, so the Hamiltonian is different than normal, and the spins $\sigma_i$ can be -1, 0, or 1. The boundary conditions are non-periodic (the chain just ends with the Nth spin)
$$Z=e^{-\beta H}$$
$$Z=\sum_{\sigma_1}...\sum_{\sigma_{N-1}}e^{\beta J\sum_{i=1}^{N-2}\sigma_i\sigma_{i+1}}\sum_{\sigma_N}e^{\beta J\sigma_{N-1}\sigma_N}$$
and here's where I get lost, I'm not sure how to evaluate this sum
 
Last edited:
Physics news on Phys.org
Update: It seems that the transfer matrix method is the best way to do this, but I cannot find the correct expression for the partition function in matrix form if the boundary conditions are not periodic.
 
I don't think I'll be of any real help, but let me try.
  • For your paper, are you explicitly not allowed to use periodic BCs?
  • Are you allowed to model it numerically?
  • Have you tried all the summation identities that seem reasonable? Maybe an index shift would work.
 
I am modelling DNA as a one dimensional Ising model, and so I don't think it would be wise to use periodic boundary conditions since DNA isn't circular (in eukaryotes anyway). The main part of my paper will be comparing this to a numerical/computational result :) Digging through the internet I've come to the conclusion that the first N-1 sums will contribute a 2, and the final sum contributes $$e^{-\beta J}+e^{\beta J}+1=2cosh(\beta J) +1$$ so the final result is $$Z=2^{N-1}(2cosh(\beta J)+1)^{N-1}$$ I'm not sure I completely understand where the $$2^{N-1}$$ comes from, so if someone could clarify that I'd really appreciate it!
 
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...
The value of H equals ## 10^{3}## in natural units, According to : https://en.wikipedia.org/wiki/Natural_units, ## t \sim 10^{-21} sec = 10^{21} Hz ##, and since ## \text{GeV} \sim 10^{24} \text{Hz } ##, ## GeV \sim 10^{24} \times 10^{-21} = 10^3 ## in natural units. So is this conversion correct? Also in the above formula, can I convert H to that natural units , since it’s a constant, while keeping k in Hz ?
Back
Top