How can I calculate the eigenvalues of a Hamiltonian with spin 1/2 objects?

JohanL
Messages
154
Reaction score
0
Find the eigenvalues of the hamiltonian

<br /> H=a(S_A \cdot S_B+S_B \cdot S_C+S_C \cdot S_D+S_D \cdot S_A)<br />

where S_A, S_B, S_C, S_D are spin 1/2 objects
_________________________

I rewrite it as

<br /> H=(1/2)*a*[(S_A+S_B+S_C+S_D)^2-(S_A+S_C)^2-(S_B+S_D)^2]<br />

then i define

<br /> J_1=S_A+S_B+S_C+S_D<br />

<br /> J_2=S_A+S_C<br />

<br /> J_3=S_B+S_D<br />

and uses

<br /> J^2_i |j_1j_2j_3;m_1m_2m_3&gt; = (h^2) j_i(j_i+1)|j_1j_2j_3;m_1m_2m_3&gt;<br />

which gives the energies

<br /> E(j_1,j_2,j_3)=(h^2/2)*a*[j_1(j_1+1)-j_2(j_2+1)-j_3(j_3+1)]<br />

Where j_1 is addition of four angular momentum of 1/2 which gives it values of 0 1, 2 and in the same way j_2 and j_3 have values of 0 1.

Am i doing this the right way? It doesn't feel so :smile:
 
Last edited:
Physics news on Phys.org
From the structure of your hamiltonian it almost looks like you could adapt transfer matrix methods, unless your spin things are vectors (I'm not clear on that). I would also say that there are restrictions on j_2 and j_3 based on j_1, but the thought process seems right.
 
Your method is completely correct. Just include the a.
Did it just seem too easy?
 
Meir Achuz said:
Your method is completely correct. Just include the a.
Did it just seem too easy?

Thank you. Yes it seemed too easy :smile:
 
How about the degeneracy of the energy levels.
For example E(010)=E(001)=E(111) and then m_1 can take on 9 different values , m_2 and m_3 5 different values. So the degeneracy of this level is 3*9*5*5 ? Is it correct so far?

But then the j_i in turn are addtions of angular momentums.
Does this add even more to the degeneracy?
 
I would count the degeneracy of each estate as the product of 2j+1 for each sub j. So I think it is 3+3+27 for your example. You shouldn't just count the m's because they are correlated to give the j estates.
 
To solve this, I first used the units to work out that a= m* a/m, i.e. t=z/λ. This would allow you to determine the time duration within an interval section by section and then add this to the previous ones to obtain the age of the respective layer. However, this would require a constant thickness per year for each interval. However, since this is most likely not the case, my next consideration was that the age must be the integral of a 1/λ(z) function, which I cannot model.
Back
Top