Processing math: 93%
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

7.4: Variational Density Matrix Method

( \newcommand{\kernel}{\mathrm{null}\,}\)




























































































































































































































































































































\( \newcommand\Dalpha

ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[1], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Dbeta
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[2], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Dgamma
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[3], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Ddelta
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[4], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Depsilon
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[5], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Dvarepsilon
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[6], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Dzeta
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[7], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Deta
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[8], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Dtheta
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[9], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Dvartheta
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[10], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Diota
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[11], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Dkappa
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[12], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Dlambda
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[13], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)





\( \newcommand\Dvarpi
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[14], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)











\( \newcommand\DGamma
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[15], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\DDelta
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[16], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\DTheta
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[17], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)









































































\( \newcommand\Vmu
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[18], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vnu
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[19], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vxi
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[20], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vom
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[21], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vpi
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[22], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vvarpi
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[23], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vrho
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[24], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vvarrho
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[25], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vsigma
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[26], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vvarsigma
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[27], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vtau
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[28], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vupsilon
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[29], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vphi
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[30], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vvarphi
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[31], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vchi
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[32], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vpsi
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[33], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\Vomega
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[34], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\VGamma
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[35], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)
\( \newcommand\VDelta
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/p[1]/span[36], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)

















\newcommand\BI{\mib I}}










































\)










































































\newcommand { M}

























\newcommand { m}














































}


















\( \newcommand\tcb{\textcolor{blue}\)
\( \newcommand\tcr{\textcolor{red}\)



































1$#1_$






















































































\newcommand\SZ{\textsf Z}} \( \newcommand\kFd{k\ns_{\RF\dar}\)

\newcommand\mutB{\tilde\mu}\ns_\ssr{B}



\( \newcommand\xhihOZ
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/span[1], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)



\( \newcommand\labar
ParseError: invalid DekiScript (click for details)
Callstack:
    at (Template:MathJaxArovas), /content/body/div/span[2], line 1, column 1
    at template()
    at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/07:_Mean_Field_Theory_of_Phase_Transitions/7.04:_Variational_Density_Matrix_Method), /content/body/p/span, line 1, column 23
\)





















The variational principle

Suppose we are given a Hamiltonian ˆH. From this we construct the free energy, F:

F=ETS=Tr(ϱˆH)+kBTTr(ϱlnϱ) .

Here, ϱ is the density matrix13. A physical density matrix must be (i) normalized ( Trϱ=1), (ii) Hermitian, and (iii) non-negative definite ( all the eigenvalues of ϱ must be non-negative).

Our goal is to extremize the free energy subject to the various constraints on ϱ. Let us assume that ϱ is diagonal in the basis of eigenstates of ˆH,

ϱ=γPγ|γγ| ,

where Pγ is the probability that the system is in state |γ. Then

F=γEγPγ+kBTγPγlnPγ .

Thus, the free energy is a function of the set {Pγ}. We now extremize F subject to the normalization constraint. This means we form the extended function

F({Pγ},λ)=F({Pγ})+λ(γPγ1) ,

and then freely extremize over both the probabilities {Pγ} as well as the Lagrange multiplier λ. This yields the Boltzmann distribution,

Peqγ=1Zexp(Eγ/kBT) ,

where Z=γeEγ/kBT=TreˆH/kBT is the canonical partition function, which is related to λ through

λ=kBT(lnZ1) .

Note that the Boltzmann weights are, appropriately, all positive.

If the spectrum of ˆH is bounded from below, our extremum should in fact yield a minimum for the free energy F. Furthermore, since we have freely minimized over all the probabilities, subject to the single normalization constraint, any distribution {Pγ} other than the equilibrium one must yield a greater value of F.

Alas, the Boltzmann distribution, while exact, is often intractable to evaluate. For one-dimensional systems, there are general methods such as the transfer matrix approach which do permit an exact evaluation of the free energy. However, beyond one dimension the situation is in general hopeless. A family of solvable (“integrable") models exists in two dimensions, but their solutions require specialized techniques and are extremely difficult. The idea behind the variational density matrix approximation is to construct a tractable trial density matrix ϱ which depends on a set of variational parameters {xα}, and to minimize with respect to this (finite) set.

Variational density matrix for the Ising model

Consider once again the Ising model Hamiltonian,

ˆH=i<jJijσiσjHiσi .

The states of the system |γ may be labeled by the values of the spin variables: |γ|σ1,σ2,. We assume the density matrix is diagonal in this basis,

ϱN(γ|γ)ϱ(γ)δγ,γ ,

where

δγ,γ=iδσi,σi .

Indeed, this is the case for the exact density matrix, which is to say the Boltzmann weight,

ϱN(σ1,σ2,)=1ZeβˆH(σ1,,σN) .

We now write a trial density matrix which is a product over contributions from independent single sites:

ϱN(σ1,σ2,)=iϱ(σi) ,

where

ϱ(σ)=(1+m2)δσ,1+(1m2)δσ,1 .

Note that we’ve changed our notation slightly. We are denoting by ϱ(σ) the corresponding diagonal element of the matrix

ϱ=(1+m2001m2) ,

and the full density matrix is a tensor product over the single site matrices:

ϱN=ϱϱϱ .

Note that ϱ and hence ϱN are appropriately normalized. The variational parameter here is m, which, if ρ is to be non-negative definite, must satisfy 1m1. The quantity m has the physical interpretation of the average spin on any given site, since

σi=σϱ(σ)σ=m.

We may now evaluate the average energy:

E=Tr(ϱNˆH)=i<jJijm2Him=12NˆJ(0)m2NHm ,

where once again ˆJ(0) is the discrete Fourier transform of J(R) at wavevector q=0. The entropy is given by

S=kBTr(ϱNlnϱN)=NkBTr(ϱlnϱ)i=NkB{(1+m2)ln(1+m2)+(1m2)ln(1m2)} .

We now define the dimensionless free energy per site: fF/NˆJ(0). We have

f(m,h,θ)=12m2hm+θ{(1+m2)ln(1+m2)+(1m2)ln(1m2)} ,

where θkBT/ˆJ(0) is the dimensionless temperature, and hH/ˆJ(0) the dimensionless magnetic field, as before. We extremize f(m) by setting

fm=0=mh+θ2ln(1+m1m) .

Solving for m, we obtain

m=tanh(m+hθ) ,

which is precisely what we found in Equation [isingmft].

[ferg] Variational field free energy Δf=f(m,h,θ)+θln2 versus magnetization m at six equally spaced temperatures interpolating between ‘high’ (θ=1.25, red) and ‘low’ (θ=0.75, blue) values. Top panel: h=0. Bottom panel: h=0.06.

Note that the optimal value of m indeed satisfies the requirement |m|1 of non-negative probability. This nonlinear equation may be solved graphically. For h=0, the unmagnetized solution m=0 always applies. However, for θ<1 there are two additional solutions at \boldsymbol{m=\pm m\nd_\ssr{A}(\theta)}, with \boldsymbol{m\nd_\ssr{A}(\theta) =\sqrt{3(1-\theta)} + \CO\big( (1-\theta)^{3/2}\big)} for t close to (but less than) one. These solutions, which are related by the Z2 symmetry of the h=0 model, are in fact the low energy solutions. This is shown clearly in figure [ferg], where the variational free energy f(m,t) is plotted as a function of m for a range of temperatures interpolating between ‘high’ and ‘low’ values. At the critical temperature θc=1, the lowest energy state changes from being unmagnetized (high temperature) to magnetized (low temperature).

For h>0, there is no longer a Z2 symmetry ( σiσi i). The high temperature solution now has m>0 (or m<0 if h<0), and this smoothly varies as t is lowered, approaching the completely polarized limit m=1 as θ0. At very high temperatures, the argument of the tanh function is small, and we may approximate tanh(x)x, in which case

m(h,θ)=hθθc .

This is called the Curie-Weiss law. One can infer θc from the high temperature susceptibility χ(θ)=(m/h)h=0 by plotting χ1 versus θ and extrapolating to obtain the θ-intercept. In our case, χ(θ)=(θθc)1. For low θ and weak h, there are two inequivalent minima in the free energy.

When m is small, it is appropriate to expand f(m,h,θ), obtaining

f(m,h,θ)=θln2hm+12(θ1)m2+θ12m4+θ30m6+θ56m8+ .

This is known as the Landau expansion of the free energy in terms of the order parameter m. An order parameter is a thermodynamic variable ϕ which distinguishes ordered and disordered phases. Typically ϕ=0 in the disordered (high temperature) phase, and ϕ0 in the ordered (low temperature) phase. When the order sets in continuously, when ϕ is continuous across θc, the phase transition is said to be second order. When ϕ changes abruptly, the transition is first order. It is also quite commonplace to observe phase transitions between two ordered states. For example, a crystal, which is an ordered state, may change its lattice structure, say from a high temperature tetragonal phase to a low temperature orthorhombic phase. When the high T phase possesses the same symmetries as the low T phase, as in the tetragonal-to-orthorhombic example, the transition may be second order. When the two symmetries are completely unrelated, for example in a hexagonal-to-tetragonal transition, or in a transition between a ferromagnet and an antiferromagnet, the transition is in general first order.

Throughout this discussion, we have assumed that the interactions Jij are predominantly ferromagnetic, Jij>0, so that all the spins prefer to align. When Jij<0, the interaction is said to be antiferromagnetic and prefers anti-alignment of the spins ( σiσj=1). Clearly not every pair of spins can be anti-aligned – there are two possible spin states and a thermodynamically extensive number of spins. But on the square lattice, for example, if the only interactions Jij are between nearest neighbors and the interactions are antiferromagnetic, then the lowest energy configuration (T=0 ground state) will be one in which spins on opposite sublattices are anti-aligned. The square lattice is bipartite – it breaks up into two interpenetrating sublattices A and B (which are themselves square lattices, rotated by 45 with respect to the original, and with a larger lattice constant by a factor of 2), such that any site in A has nearest neighbors in B, and vice versa. The honeycomb lattice is another example of a bipartite lattice. So is the simple cubic lattice. The triangular lattice, however, is not bipartite (it is tripartite). Consequently, with nearest neighbor antiferromagnetic interactions, the triangular lattice Ising model is highly frustrated. The moral of the story is this: antiferromagnetic interactions can give rise to complicated magnetic ordering, and, when frustrated by the lattice geometry, may have finite specific entropy even at T=0.

Mean Field Theory of the Potts Model

The Hamiltonian for the Potts model is

ˆH=i<jJijδσi,σjHiδσi,1 .

Here, σi{1,,q}, with integer q. This is the so-called ‘q-state Potts model’. The quantity H is analogous to an external magnetic field, and preferentially aligns (for H>0) the local spins in the σ=1 direction. We will assume H0.

The q-component set is conveniently taken to be the integers from 1 to q, but it could be anything, such as

σi{tomato, penny, ostrich, Grateful Dead ticket from 1987, } .

The interaction energy is Jij if sites i and j contain the same object (q possibilities), and 0 if i and j contain different objects (q2q possibilities).

The two-state Potts model is equivalent to the Ising model. Let the allowed values of σ be ±1. Then the quantity

\[\delta_{\sigma,\sigma'}=\half + \half\,\sigma\sigma'\\]

equals 1 if σ=σ, and is zero otherwise. The three-state Potts model cannot be written as a simple three-state Ising model, one with a bilinear interaction σσ where σ{1,0,+1}. However, it is straightforward to verify the identity

δσ,σ=1+12σσ+32σ2σ2(σ2+σ2) .

Thus, the q=3-state Potts model is equivalent to a S=1 (three-state) Ising model which includes both bilinear (σσ) and biquadratic (σ2σ2) interactions, as well as a local field term which couples to the square of the spin, σ2. In general one can find such correspondences for higher q Potts models, but, as should be expected, the interactions become increasingly complex, with bi-cubic, bi-quartic, bi-quintic, terms. Such a formulation, however, obscures the beautiful Sq symmetry inherent in the model, where Sq is the permutation group on q symbols, which has q! elements.

Getting back to the mean field theory, we write the single site variational density matrix ϱ as a diagonal matrix with entries

ϱ(σ)=xδσ,1+(1xq1)(1δσ,1) ,

with ϱN(σ1,,σN)=ϱ(σ1)ϱ(σN). Note that Tr(ϱ)=1. The variational parameter is x. When x=q1, all states are equally probable. But for x>q1, the state σ=1 is preferred, and the other (q1) states have identical but smaller probabilities. It is a simple matter to compute the energy and entropy:

E=Tr(ϱNˆH)=12NˆJ(0){x2+(1x)2q1}NHxS=kBTr(ϱNlnϱN)=NkB{xlnx+(1x)ln(1xq1)} .

The dimensionless free energy per site is then

f(x,θ,h)=12{x2+(1x)2q1}+θ{xlnx+(1x)ln(1xq1)}hx ,

where h=H/ˆJ(0). We now extremize with respect to x to obtain the mean field equation,

fx=0=x+1xq1+θlnxθln(1xq1)h .

Note that for h=0, x=q1 is a solution, corresponding to a disordered state in which all states are equally probable. At high temperatures, for small h, we expect xq1h. Indeed, using Mathematica  one can set

xq1+s ,

and expand the mean field equation in powers of s. One obtains

h=q(qθ1)q1s+q3(q2)θ2(q1)2s2+O(s3) .

For weak fields, |h|1, and we have

s(θ)=(q1)hq(qθ1)+O(h2) ,

which again is of the Curie-Weiss form. The difference s=xq1 is the order parameter for the transition.

Finally, one can expand the free energy in powers of s, obtaining the Landau expansion,

f(s,θ,h)=2h+12qθlnqhs+q(qθ1)2(q1)s2(q2)q3θ6(q1)2s3+q3θ12[1+(q1)3]s4q4θ20[1(q1)4]s5+q5θ30[1+(q1)5]s6+.

Note that, for q=2, the coefficients of s3, s5, and higher order odd powers of s vanish in the Landau expansion. This is consistent with what we found for the Ising model, and is related to the Z2 symmetry of that model. For q>3, there is a cubic term in the mean field free energy, and thus we generically expect a first order transition, as we shall see below when we discuss Landau theory.

Mean Field Theory of the XY Model

Consider the so-called XY model, in which each site contains a continuous planar spin, represented by an angular variable \phi_i\in [-\pi,\pi] :

\HH=-{1\over 2}\sum_{i\ne j} J_{ij}\cos\!\big(\phi_i-\phi_j\big)-H\sum_i\cos\phi_i\ . \label{XYmodel}

We write the (diagonal elements of the) full density matrix once again as a product:

\vrh\ns_N(\phi_1,\phi_2,\ldots)=\prod_i\vrh(\phi_i)\ .

Our goal will be to extremize the free energy with respect to the function \vrh(\phi). To this end, we compute

E=\Tra (\vrh\ns_N\,\HH) = - \half N\jhz\> \Big|\! \Tra\!\big(\vrh\,e^{i\phi}\big) \Big|^2 - NH\Tra\!\big(\vrh\,\cos\phi\big)\ .

The entropy is

S=-N\kB\Tra (\vrh\,\ln\vrh)\ .

Note that for any function A(\phi), we have14

\Tra\big(\vrh \, A)\equiv\int\limits_{-\pi}^\pi\!{d\phi\over 2\pi}\>\vrh(\phi)\,A(\phi)\ .

We now extremize the functional F\big[\vrh(\phi)\big]=E-TS with respect to \vrh(\phi), under the condition that \Tra\vrh=1. We therefore use Lagrange’s method of undetermined multipliers, writing

F^*=F-N\kT\,\lambda\>\Big(\Tra\!\vrh -1\Big)\ .

Note that F^* is a function of the Lagrange multiplier \lambda and a functional of the density matrix \vrh(\phi). The prefactor N\kT which multiplies \lambda is of no mathematical consequence – we could always redefine the multiplier to be \lambda'\equiv N\kT\lambda. It is present only to maintain homogeneity and proper dimensionality of F^* with \lambda itself dimensionless and of order N^0. We now have

\begin{split} {\delta F^*\over\delta\vrh(\phi)}&={\delta\over\delta\vrh(\phi)} \Bigg\{\!- \half N\jhz\> \Big|\! \Tra\!\! \big(\vrh\,e^{i\phi}\big) \Big|^2 - NH\Tra\!\!\big(\vrh\,\cos\phi\big)\\ &\hskip 1.25in+ N\kT\,\Tra\!\big(\vrh\,\ln\vrh\big)-N\kT\>\lambda\, \Big(\Tra\!\vrh -1\Big)\Bigg\} \ . \end{split}

To this end, we note that

{\delta\over\delta\vrh(\phi)} \Tra(\vrh\,A)={\delta\over\delta\vrh(\phi)} \int\limits_{-\pi}^\pi\!{d\phi\over 2\pi}\>\vrh(\phi)\,A(\phi)={1\over 2\pi}\>A(\phi)\ .

Thus, we have

\begin{split} {\delta {\tilde F}\over\delta\vrh(\phi)}&=-\half N\jhz\cdot{1\over 2\pi}\Bigg[ \mathop{\textsf{Tr}}_{\phi'}\!\big(\vrh \>e^{i\phi'}\big) \>e^{-i\phi} +\mathop{\textsf{Tr}}_{\phi'}\!\big(\vrh \>e^{-i\phi'}\big)\>e^{i\phi} \Bigg]-NH\cdot{\cos\phi\over 2\pi}\\ &\hskip 1.5in + N\kT\cdot {1\over 2\pi}\,\Big[\ln\vrh(\phi)+1\Big] -N\kT\cdot{\lambda\over 2\pi}\ . \end{split}

Now let us define

\mathop{\textsf{Tr}}_{\phi}\!\big(\vrh\,e^{i\phi}\big)=\int\limits_{-\pi}^\pi\!{d\phi\over 2\pi}\> \vrh(\phi)\>e^{i\phi}\equiv m\,e^{i\phi\ns_0}\ .

We then have

\ln\vrh(\phi)={\jhz\over\kT}\>m\,\cos(\phi-\phi\ns_0) + {H\over\kT}\,\cos\phi+ \lambda-1 .

Clearly the free energy will be reduced if \phi\ns_0=0 so that the mean field is maximal and aligns with the external field, which prefers \phi=0. Thus, we conclude

\vrh(\phi)=\CC\>\exp\!\bigg({H\ns_{eff}\over\kT}\>\cos\phi\bigg)\ ,

where

H\ns_{eff}=\jhz\, m + H

and \CC=e^{\lambda-1}. The value of \lambda is then determined by invoking the constraint,

\Tra\!\vrh=1=\CC\!\!\int\limits_{-\pi}^\pi\!{d\phi\over 2\pi}\>\exp\!\bigg({H\ns_{eff}\over\kT}\>\cos\phi\bigg) =\CC\,I\ns_0(H\ns_{eff}/\kT)\ ,

where I_0(z) is the Bessel function. We are free to define \ve\equiv H\ns_{eff}/\kT, and treat \ve as our single variational parameter. We then have the normalized single site density matrix

\vrh(\phi)={\exp(\ve\cos\phi)\over\int\limits_{-\pi}^\pi \!\!{d\phi'\over 2\pi}\,\exp(\ve\cos\phi')} = {\exp(\ve\cos\phi)\over I_0(\ve)}\ .

We next compute the following averages:

\begin{aligned} \blangle\, e^{\pm i\phi}\,\brangle & = \int\limits_{-\pi}^\pi\!{d\phi\over 2\pi}\,\vrh(\phi)\,e^{\pm i\phi} = {I_1(\ve)\over I_0(\ve)}\\ \big\langle\!\cos(\phi-\phi')\big\rangle&= \Rep\,\big\langle e^{i\phi}\,e^{-i\phi'}\big\rangle =\bigg({I_1(\ve)\over I_0(\ve)}\bigg)^{\!2}\ ,\end{aligned}

as well as

\Tra(\vrh\,\ln\vrh)=\int\limits_{-\pi}^\pi\!{d\phi\over 2\pi}\,{e^{\ve\cos\phi}\over I_0(\ve)}\, \Big\{\ve\cos\phi - \ln I_0(\ve)\Big\} = \ve\,{I_1(\ve)\over I_0(\ve)} - \ln I_0(\ve)\ .

The dimensionless free energy per site is therefore

f(\ve,h,\theta)=-{1\over 2}\bigg({I_1(\ve)\over I_0(\ve)}\bigg)^{\!2} +(\theta\ve-h)\,{I_1(\ve)\over I_0(\ve)} - \theta\,\ln I_0(\ve)\ ,

with \theta=\kT/\jhz and h=H/\jhz and f=F/N\jhz as before. Note that the mean field equation is m=\theta\ve-h=\blangle\, e^{i\phi}\,\brangle ,

\theta\ve - h = {I_1(\ve)\over I_0(\ve)}\quad. \label{XYvdm}

For small \ve, we may expand the Bessel functions, using

I_\nu(z)=(\half z)^\nu\sum_{k=0}^\infty {(\frac{1}{4} z^2)^k\over k!\,\RGamma(k+\nu+1)}\ ,

to obtain

f(\ve,h,\theta)=\frac{1}{4}\,\big(\theta-\half\big)\,\ve^2 + \frac{1}{64}\,\big(2-3\theta\big)\,\ve^4 -\half\,h\ve + \frac{1}{16}\,h\ve^3 + \ldots\ .

This predicts a second order phase transition at \theta_\Rc=\half.15 Note also the Curie-Weiss form of the susceptibility at high \theta:

{\pz f\over\pz \ve}=0 \quad \Longrightarrow\quad \ve={h\over \theta-\theta_\Rc} + \ldots\ .

XY model via neglect of fluctuations method

Consider again the Hamiltonian of Equation [XYmodel]. Define z\ns_i\equiv \exp(i\phi\ns_i) and write

z\ns_i=w + \delta z\ns_i\quad,

where w\equiv\langle z\ns_i\rangle and \delta z\ns_i\equiv z\ns_i-w. Of course we also have the complex conjugate relations z^*_i=w^*+\delta z^*_i and w^*=\langle z^*_i\rangle. Writing \cos(\phi\ns_i-\phi\ns_j)=\Rep(z^*_i z\ns_j) , by neglecting the terms proportional to \delta z^*_i\,\delta z\ns_j in \HH we arrive at the mean field Hamiltonian,

\HH^\ssr{MF}=\half N \HJ(0)\,|w|^2 - \half\HJ(0)\,|w|\sum_i\big(w^* z\ns_i + w z^*_i\big) - \half H\sum_i \big(z^*_i+z\ns_i\big)

It is clear that the free energy will be minimized if the mean field w breaks the \SO(2) symmetry in the same direction as the external field H, which means w\in\MR and

\HH^\ssr{MF}=\half N \HJ(0)\,|w|^2 - \big(H+\HJ(0)\,|w|\big)\sum_i\cos\phi\ns_i\quad.

The dimensionless free energy per site is then

f=\half |w|^2 -\theta\ln I\ns_0\bigg({h+|w|\over\theta}\bigg)\quad.

Differentiating with respect to |w| , one obtains

|w|\equiv m = {I\ns_1\!\left({h+m\over\theta}\right)\over I\ns_0\!\left({h+m\over\theta}\right)}\quad,

which is the same equation as Equation [XYvdm]. The two mean field theories yield the same results in every detail (see §10).


This page titled 7.4: Variational Density Matrix Method is shared under a CC BY-NC-SA license and was authored, remixed, and/or curated by Daniel Arovas.

Support Center

How can we help?