Register to reply 
Partition function as a description of the system 
Share this thread: 
#1
Dec1313, 06:19 PM

P: 96

Hi,
Let, $$\hat{H} = a\hat{S_x} + b\hat{S_z}$$where Sx,Sz are the spin operators, a,b constants. Assume the system is coupled to a reservoir. For clarity, Let $$\hbar=\beta=1$$ The density matrix is $$ρ=\frac{e^{β\hat{H}}}{Z}= \frac{1}{Z} \left(\begin{array}{cc}0&e^{a/2}\\e^{a/2}&0\end{array}\right) + \frac{1}{Z} \left(\begin{array}{cc}e^{b/2}&0\\0&e^{b/2}\end{array}\right) = \frac{1}{Z} \left(\begin{array}{cc}e^{b/2}&e^{a/2}\\e^{a/2}&e^{b/2}\end{array}\right)$$ And $$Z=tr(e^{β\hat{H}})=2\cosh{\frac{b}{2}}$$ The partition function is a description of the statistical properties of the system. But here it is ignorant to the fact that H also contains Sx part. The constant 'a' does not come into play at all. How, then, can we depend on Z to derive our statistical functions for the system? 


#2
Dec1413, 07:39 AM

HW Helper
P: 3,443

I don't understand how ##\exp(\beta \hat{H} )## is equal the thing you have written out on the right hand side. When there is an exponential of a matrix, you need to write it out as the Taylor series of an exponential, where each term will contain the matrix, to some power.
This will give an infinite series, which is not something that is usually 'nice'. But you can use the properties of the Pauli matrices to get a 'nice' answer. This is not obvious if you have not been shown how to do this problem before. 


#3
Dec1413, 11:37 AM

Thanks
P: 1,948

Also, the calculation of the trace doesn't look right. The trace is supposed to be taken over the space of eigenfunctions of your solutions. That has nothing to do with the trace of spin matrix. For instance, a harmonic oscillator doesn't have spin matrix at all but it has an infinite number of eigenstates. the trace is taken over that infinite space.



#4
Dec1413, 02:25 PM

Sci Advisor
P: 1,947

Partition function as a description of the system
Your computation of [itex]e^{\beta H}[/itex] isn't correct, I don't think.
An easy way to compute [itex]e^{\beta H}[/itex] is to first figure out how to diagonalize [itex]H[/itex]. That is, you find some matrix [itex]U[/itex] such that [itex]U H U^{1} = D[/itex] where [itex]D[/itex] is a diagonal matrix, with entries [itex]\lambda_i[/itex]. (Those [itex]\lambda_i[/itex] are just the eigenvalues of [itex]H[/itex]) Then: [itex]e^{\beta H} = U^{1} (e^{\beta D}) U[/itex] Now, one of the magical properties of the trace is that for any matrices [itex]A, B, C[/itex], [itex]tr(A B C) = tr(B C A)[/itex] So [itex]tr(e^{\beta H}) = tr(U^{1} (e^{\beta D}) U) = tr(e^{\beta D} U U^{1}) = tr(e^{\beta D})[/itex] Computing [itex] e^{\beta D}[/itex] is easy: it's just a diagonal matrix with entries [itex]e^{ \beta \lambda_i}[/itex]. The trace of a matrix is just the sum of the diagonal entries, so So [itex]tr(e^{\beta H}) = tr(e^{\beta D}) = \sum_i e^{\beta \lambda_i}[/itex] For the particular [itex]H[/itex] that you're talking about, I believe it's true that the two eigenvalues are: [itex]\lambda_1 = \sqrt{a^2 + b^2}[/itex] and [itex]\lambda_2 =  \sqrt{a^2 + b^2}[/itex] So [itex]tr(e^{\beta H}) = e^{\beta \sqrt{a^2 + b^2}} + e^{+\beta \sqrt{a^2 + b^2}} = 2 cosh(\beta \sqrt{a^2 + b^2})[/itex] 


#5
Dec1413, 02:46 PM

P: 96

Yes I think I missed up things while trying to use Latex.



#6
Dec1413, 03:57 PM

HW Helper
P: 3,443




Register to reply 
Related Discussions  
Construct a partition function for the system  Advanced Physics Homework  1  
Partition Function for Thermodynamic System  Advanced Physics Homework  1  
Magnetic system, partition function  Advanced Physics Homework  2  
Partition Function of 2 State System  Introductory Physics Homework  1  
Partition Function of Fermion System  Advanced Physics Homework  0 