물리2010. 2. 10. 16:33

이 글은 Dung-Hai Lee 교수의 critical phenomena라는 책(?), 강의노트(?)를 기반으로 스터디한 글입니다.

고수님들 틀린 부분있으면 확실하게 지적해주시기 바랍니다. ('-')/

2.2절 입니다. 전 절에서 썼던 식 

<S>=\tanh\beta(h+Jz<S>) 로 시작해 보도록 하죠. 

우선 우리는 온도 T\to T_c^- 일때 m_0\to 0  즉, <S>\to 0 이므로 \beta Jz<S>도 0으로 수렴합니다. 그러면 tanh함수를 테일러 급수 전개(talyor expansion)를 할수 있습니다. 

<S>\approx\beta Jz<S>-\frac{1}{3}(\beta Jz<S>)^3+\ldots

여기서 <S>를 m_0라고 생각하고 바꿈니다. 스핀의 기대값이나 자기밀도나 같은 값이라고 이야기해도 괜찮겠죠. 여기서 m_0는 자발적 자기밀도(spontaneous magnetization)입니다. 그리고

T_c \equiv \frac{Jz}{k_B}

라고 정의 했었죠. 

\left(\beta=\frac{1}{k_B T}\right) 이구요.

Jz=k_B T_c 로 바꾸면

m_0=\frac{T_c}{T}m_0-\frac{1}{3}\left(\frac{T_c}{T}\right)^3 m_0^3 + \ldots

1=\frac{T_c}{T}-\frac{1}{3}\left(\frac{T_c}{T}\right)^3 m_0^2 + \mathcal{O}(m_0^4) + \ldots

의 형태가 됩니다. 여기서 m_0\approx0 이므로 4승항 이하로는 제거하도록 하죠. 그리고 나서 m_0에 대해서 정리하게 되면 

m_0=\sqrt{3\frac{T^3}{T_c^3}\left( \frac{T_c}{T} -1 \right)}

이렇게 됨니다. 우리는 T\to T_c^-일때를 이야기 했기 때문에 앞에 \frac{T^3}{T_c^3} 항은 1로 가게 되고 뒤에 괄호 안은 \frac{T_c-T}{T_c}이 됨니다. 온도가 임계온도에 가까워지는 경우 값이 0으로 수렴하게 되죠. 

t\equiv\frac{T-T_c}{T_c}  라는 값으로 1장에서 정의 했습니다. 이 양은 0이 될수록 임계온도에 온도가 가까워지게 되는 것이죠.

결론은 T\to T_c^-일때 m_0값은 0으로 수렴하게 됨니다. 온도가 임계온도에 가까워지면 자성적 성질이 완전히 사라 지게 되는것이죠. 

m_0 \sim |t|^{\frac{1}{2}}

라는 식으로 1장에서 이야기한 m_0 \sim |t|^{\beta} 에서 \beta=\frac{1}{2} 값을 가지게 됨니다. 실제 실험 값이 0.32 정도 인데요. 조금 멀어보이긴 하지만 그래도 대략적으로 조금은 가까워진거 같네요.

다음으로는 임계온도 때에서 외부 자기장 값이 존재할때의 경우를 생각해 보도록 하죠. 여기선 스핀 기대값과 자기장의 관계를 보기위해 이런 접근을 합니다.  그리고 온도가 임계온도이기 때문에 \beta를  \beta_c 라고 합니다. 

<S>=\tanh\beta_c(h+Jz<S>)

=\tanh(<S>+\beta_c h)

여기서는 h\neq0 이지만 후에 h\to0 이기 때문에 외부 자기장 h를 충분히 작다고 생각하고 문제에 접근하면 같은 방식으로 테일러 전개를 할수 있게 됨니다. 

<S>=<S>+\beta_c h-\frac{1}{3}\left(<S>+\beta_c h\right)^3+\ldots

뒤에 5승 이하를 제거 하고 정리하면

\beta_c h =\frac{1}{3}\left(<S>+\beta_c h\right)^3

3(\beta_c h)^\frac{1}{3} -\beta_c h =<S>

여기서  h\to0라고 가정 하기 때문에 1/3승인 h가 그냥 h보다 훨신 큰 영향을 끼치게 됨니다. (dominant)

<S>\sim h^\frac{1}{3}

이라는 결론이 나옴니다. 이런 경우 임계온도에서 외부 자기장의 세기에 따른 스핀의 기대값이 외부 자기장의 1/3승에 비례한다는 것을 알수 있는데요. 실제 실험값은 1/4.80이니까 대략 비슷한긴 하지만 조금 멀어 보이긴 합니다. 조금씩 평균장 이론(mean-field theory)의 한계가 보이네요.

이제 온도가 임계온도 보다 조금더 높은 부분에서의 상황을 보겠습니다. 여기서 \beta h\ll1 인 상황에서만 생각합니다.  이는 열적 현상에 비해 자기장을 지나치게 세게 가하지 않는 다는 것인데요.. 상자성도 강력한 자기장을 가하게 되면 어느정도 자성을 띄기 때문입니다. 

<S>=\tanh\beta(h+Jz<S>)

<S>=0이기 때문에 위의 식은 다시 쉽게 근사 됨니다. 

<S>\approx\beta(h+ Jz<S>)-\frac{1}{3}\left[\beta(h+ Jz<S>)\right]^3+\ldots

\beta(h+Jz<S>)\ll1 이라고 위에서 가정 했기때문에 우선 3승항 이하로는 모두 제거해 보도록 하죠. 

<S>\approx\beta(h+ Jz<S>)

자 이녀석을 잘 요리해 봅니다. 

k_B T<S>\approx(h+ Jz<S>)

k_B T<S>\approx h+ k_B T_c <S>

k_B (T-T_c)<S>\approx h

<S>\approx \frac{h}{k_B (T-T_c)}

온도과 스핀의 기대값은 임계온도보다 조금 높은 부분에서는 이러한 관계를 가집니다. 여기서 스핀의 기대값은 자화율(magnetic susceptibility)와 관계가 있습니다. 상자성을 띄는 상황에서 온도를 얼마나 높여 주느냐에 따라(외부 자기장이 있는 상황에서) 자화가 급격히 낮아지니까요. 강자성에는 스핀이 임계온도에 가까워지지 않으면 모두 한방향을 가르키기 때문에 자화율을 따지는게 의미가 없습니다. 

이 식에서는 온도와 자화율은 \chi(T)\sim T^{-1}를 띄는군요. 하지만 실제 실험값은 1.24의 지수를 가진 멱급수(pow-law)를 따름니다. 그래도 조금은 비슷하네요.

자.. 책에선 이 부분이 가장 논란이 되는 부분인데요. 다음이 1절에서 언급한 헤밀토니안입니다.  

\mathcal{H}=-J\sum_{ij}C_{i j}S_iS_j-h\sum_{i}S_i

처음에 썼던 에너지 연산자(Hamiltonian)을 수정해 줍니다. 전에는 두번 계산(double counting)이 됬덨다고 하네요.  이쪽 공부 해보신 분들은 알겠지만 실제로 각 격자자리(lattices site)들을 지나면서 계산을 하다보면 1번 위치에서 1번 스핀과 2번 스핀 사이에 관계를 계산하고 나서 다시 2번의 위치에서 1번스핀과 2번 스핀의 관계를 또 계산하게 되기때문에 그런 반복된 계산을 제거하기 위해 반으로 나누어 줍니다.

문제는 지금까지 이전의 헤밀토니안으로 이것도 하고 저것도 하면서 잘 써오다가 갑자기 이런식으로 변덕을 주는게 좀 많이 기이하지 않을수 없습니다. 

\mathcal{H}_{mf}=-\sum_i h_{eff,i}S_i + \frac{J}{2}\sum_{i,j}C_{ij}<S_i><S_j>

이전 헤밀토니안이 음의 부호가 있기 때문에 저 항을 더해주면 두번 더해준 부분을 제거 할수 있습니다. 책에서 보면은 저 부분의 항은 자화(magnetization)에는 영향을 주지 않는다고 하면서 이르 프리에너지(free energy)에 넣습니다. 이 말 부분도 석연치가 않습니다. 확실히 저 부분이 임계온도에 영향을 주기때문에.. 자화에도 어느정도 영향을 줄텐데 말이죠. 물런 임계온도와의 상대적인 값(t)만을 고려한다면 상관은 없습니다.  이제 슬슬 고통스러운 계산의 시작입니다. 

f=-\frac{k_B T}{N}\ln Z

기본적으로 자유에너지(free energy)는 위와 같이 정의 됨니다. 분배함수(partition function)을 적으면,

f_{mf}=-\frac{k_B T}{N}\ln\lbrace Tr[e^{-\beta\mathcal{H}_{mf}}] \rbrace

이와 같이 됨니다. 그럼 헤밀토니안을 넣고 계산을 해보죠. 

f_{mf}=-\frac{k_B T}{N}\ln\lbrace \sum_{S_i}\left[e^{\beta \sum_i h_{eff,i}S_i - \frac{J}{2}\beta\sum_{i,j}C_{ij}<S_i><S_j> }\right] \rbrace

글씨가 작아져서 잘 안보이기 시작하네요. 

f_{mf}=-\frac{k_B T}{N}\ln\lbrace \sum_{S_i}\left[e^{\beta \sum_i h_{eff,i}S_i} e^{- \frac{J}{2}\beta\sum_{i,j}C_{ij}<S_i><S_j> }\right] \rbrace

뒤에 두번째 지수(exp)항은 앞에 S_i랑은 관계가 없으므로 로그를 이용해서 제외 시켜 줌니다.

f_{mf}=-\frac{k_B T}{N}\ln\lbrace \sum_{S_i}\left[e^{\beta \sum_i h_{eff,i}S_i} \right] \rbrace -\frac{k_B T}{N}\ln e^{-\frac{J}{2}\beta\sum_{i,j}C_{ij}<S_i><S_j> }

f_{mf}=-\frac{k_B T}{N}\ln\lbrace \sum_{S_i}\left[e^{\beta \sum_i h_{eff,i}S_i} \right] \rbrace + \frac{J}{2N} \sum_{i,j}C_{ij}<S_i><S_j> }

f_{mf}=-\frac{k_B T}{N}\ln\lbrace \sum_{S_1}\sum_{S_2}\sum_{S_3}\ldots\sum_{S_N}\left[e^{\beta h_{eff,1}S_1}e^{\beta h_{eff,2}S_2}e^{\beta h_{eff,3}S_3}\ldots e^{\beta h_{eff,N}S_N} \right] \rbrace \\+ \frac{J}{2N} \sum_{i,j}C_{ij}<S_i><S_j> }

f_{mf}=-\frac{k_B T}{N}\ln\lbrace \sum_{S_1}e^{\beta h_{eff,1}S_1}\sum_{S_2}e^{\beta h_{eff,2}S_2}\sum_{S_3}e^{\beta h_{eff,3}S_3}\ldots \sum_{S_N}e^{\beta h_{eff,N}S_N} \rbrace\\ + \frac{J}{2N} \sum_{i,j}C_{ij}<S_i><S_j> }

f_{mf}=-\frac{k_B T}{N}\left[\ln \sum_{S_1}e^{\beta h_{eff,1}S_1}+\ln\sum_{S_2}e^{\beta h_{eff,2}S_2}+\ln\sum_{S_3}e^{\beta h_{eff,3}S_3}+\ldots+\ln \sum_{S_N}e^{\beta h_{eff,N}S_N}\right] \\+ \frac{J}{2N} \sum_{i,j}C_{ij}<S_i><S_j> }

f_{mf}=-\frac{k_B T}{N}\left[\ln(e^{\beta h_{eff,1}}-e^{-\beta h_{eff,1}})+\ln(e^{\beta h_{eff,2}}-e^{-\beta h_{eff,2}} ) + \ldots +\ln(e^{\beta h_{eff,N}}-e^{-\beta h_{eff,N}})\right] \\ + \frac{J}{2N} \sum_{i,j}C_{ij}<S_i><S_j> }f_{mf}=-\frac{k_B T}{N}\left[\ln(2\cosh{\beta h_{eff,1}})+\ln(2\cosh{\beta h_{eff,2}})+ \ldots +\ln(2\cosh{\beta h_{eff,N}})\right] \\+ \frac{J}{2N} \sum_{i,j}C_{ij}<S_i><S_j> }

f_{mf}=\frac{J}{2N} \sum_{i,j}C_{ij}<S_i><S_j> }-\frac{k_B T}{N}\sum_i\ln\left[2\cosh{\beta h_{eff,i}}\right]

으로 간단하게 정리가 되네요. 휴... 솔직히 중간에 곱기호를 써서 하면 더 편하지만 이해를 돕기위해서 펼쳤는데 앞으론 간단하게 하겠습니다. 

외부 자기장이 균일한 경우에는 전에 봤던대로 간단한 형태로 바꿀수 있습니다. 좀더 복잡한 것을 간단한 것에서 부터 잘 다룰수 있게 됨니다. 

f_{mf}=\frac{Jz}{2}<S>^2 -k_B T\ln\left[2\cosh{\beta \left(h+Jz<S>\right) }\right]

i 값과 관계 없이 모든 자리(site)들이 같기 때문에 같은것을 N번 더해주는 것이 되서 N은 소거 됨니다. 

T<T_c 이고 h=0일때는 전부터 언급했던것처럼 물질이 <S>=\pm m_0를 가지게 됨니다. 

f_{mf}=\frac{Jz}{2}<S>^2 -k_B T\ln\left[2\cosh{\beta \left(Jz<S>\right) }\right]

우선 자유 에너지의 대략적인 모습을 살펴볼 필요가 있는데요. 이를 위해 <S>이 충분히 작다고 가정하고 그때의 자유 에너지를 근사 할수 있습니다. Jz=k_B T로 바꾸고 두번째 항을 ln와 cosh항으로 테일러 급수 전개(talyor expansion)를 합니다. 그러면

\cosh x=1+\frac{x^2}{2!}+\frac{x^4}{4!}+\ldots

\ln(1+x)=x-\frac{x^2}{2}+\frac{x^3}{3}+\ldots

이런 전개를 적용해서

f_{mf}=\frac{Jz}{2}<S>^2 -k_B T\ln\left[2 + 2\frac{\beta^2 \left(Jz<S>\right)^2}{2!} }+\ldots\right]

 

\approx \frac{k_B T_c}{2}<S>^2 -k_B T\ln2 -k_B T\ln\left[1 + \frac{\beta^2 \left(Jz<S>\right)^2}{2} }\right]

\approx \frac{k_B T_c}{2}<S>^2 -k_B T\ln2 -k_B T\left[\frac{\beta^2 \left(Jz<S>\right)^2}{2} -\frac{\beta^4 \left(Jz<S>\right)^4}{8} \right]

=\frac{k_B T_c}{2}<S>^2 -k_B T\ln2 -k_B T\left[\frac{\beta^2 \left(k_B T_c<S>\right)^2}{2} -\frac{\beta^4 \left(k_B T_c<S>\right)^4}{8} \right]

=\frac{k_B T_c}{2}<S>^2 -k_B T\ln2 -k_B T\frac{\beta^2 k_B^2 T_c^2 <S>^2}{2} -k_B T\frac{\beta^4 k_B^4 T_c^4 <S>^4}{8}

=\frac{k_B T_c}{2}<S>^2 -k_B T\ln2 -\frac{ k_B T_c^2 <S>^2}{2T} -\frac{k_B T_c^4 <S>^4}{8T^3}

=\frac{k_B}{2}\left(\frac{T_c}{T}\right)T<S>^2 -k_B T\ln2 -\frac{k_B}{2}\left(\frac{T_c}{T}\right)T_c<S>^2 -\frac{k_B T}{8} \left(\frac{T_c^4}{T^4}\right)<S>^4}

=-k_B T\ln2 -\frac{k_B}{2}\left(\frac{T_c}{T}\right)(T_c-T)<S>^2 -\frac{k_B T}{8} \left(\frac{T_c^4}{T^4}\right)<S>^4}

=-k_B \lbrace T\ln2 +\frac{T_c}{2T}(T_c-T)<S>^2 +\frac{T}{8} \left(\frac{T_c^4}{T^4}\right)<S>^4}\rbrace

f_{mf}=-k_B \lbrace T\ln2 +\frac{T_c}{2T}(T_c-T)<S>^2 +\frac{T}{8} \left(\frac{T_c}{T}\right)^4<S>^4}\rbrace

이것이 계산한 식이고.

f_{mf}=-k_B \lbrace T\ln2+\frac{T_c}{2T}(T-T_c) <S>^2 +\frac{T}{12}\left(\frac{T_c}{T}\right)^4<S>^4 \rbrace

이건 책에 있는 식인데요. 확실히 계산한 식이 맞는거 같습니다. 책이 틀렸습니다. 

여기서 알수 있는것은 free energy가 우함수(even function)이라는 것입니다. 그리고 온도가 임계온도보다 낮아 지는 경우, 두번째 항이 음수가 되서 <S>\neq0 이어도 free engergy가 최소가 될수 있습니다. 임계온도가 높을때는 <S>=0 인 경우에만 에너지가 최소가 되죠. 그럼 임계온도보다 낮은 부분에서 free energy를 온도에 변화에 따라 생각해 보도록 하죠. 여기서 중요한건 온도는 임계온도 보다 낮지만 <S>\neq0 이기 때문에 온도는 임계온도에 충분히 가까워야 하겠죠?

\left. \frac{\partial{f_{mf}}}{\partial{<S>}} \right|_{<S>=m}=0

이런 조건을 만족하는 m을 찻습니다. 이것도 한번 미분한 다음에 <S>의 제곱에 대해서 정리하면 돼는데요. 이건 인내심을 가지고 풀어보도록 하죠.

\frac{\partial{f_{mf}}}{\partial{<S>}}=-k_B \lbrace 2\frac{T_c}{2T}(T_c-T) <S> +3\frac{T}{8}\left(\frac{T_c}{T}\right)^4<S>^3 \rbrace =0

\lbrace 2\frac{T_c}{2T}(T_c-T) <S> +3\frac{T}{8}\left(\frac{T_c}{T}\right)^4<S>^3 \rbrace =0

\frac{T_c}{T}(T_c-T) <S> +\frac{3T}{8}\left(\frac{T_c}{T}\right)^4<S>^3 =0

\frac{T_c}{T}(T-T_c) <S> =\frac{3T}{8}\left(\frac{T_c}{T}\right)^4<S>^3

(T-T_c) =\frac{3T}{8}\left(\frac{T_c}{T}\right)^3<S>^2

<S>^2=\frac{8}{3T}\left(\frac{T}{T_c}\right)^3(T-T_c)

제가 계산한 식인데요. 책이랑 달라지기 시작합니다.

이것을 다시 처음 free energy식에 넣어주면은

 

f_{mf}=-k_B \lbrace T\ln2 +\frac{T_c}{2T}(T_c-T)<S>^2 +\frac{T}{8} \left(\frac{T_c}{T}\right)^4<S>^4}\rbrace

f_{mf}=-k_B \lbrace T\ln2 +\frac{T_c}{2T}(T_c-T)\frac{8}{3T}\left(\frac{T}{T_c}\right)^3(T-T_c) +\frac{T}{8} \left(\frac{T_c}{T}\right)^4\left(\frac{8}{3T}\left(\frac{T}{T_c}\right)^3(T-T_c)\right)^2}\rbrace

f_{mf}=-k_B \lbrace T\ln2 -(T_c-T)^2\frac{4}{3T}\left(\frac{T}{T_c}\right)^2 +\frac{T}{8} \left(\frac{T_c}{T}\right)^4\frac{8^2}{9T^2}\left(\frac{T}{T_c}\right)^6(T_c-T)^2}\rbrace

f_{mf}=-k_B \lbrace T\ln2 -(T_c-T)^2\frac{4}{3T}\left(\frac{T}{T_c}\right)^2 + \frac{8}{9T}\left(\frac{T}{T_c}\right)^2(T_c-T)^2}\rbrace

f_{mf}=-k_B \lbrace T\ln2 - \frac{8}{9T}\left(\frac{T}{T_c}\right)^2(T_c-T)^2}\rbrace

f_{mf}=-k_B \lbrace T\ln2 - \frac{8}{9}\left(\frac{T}{T_c^2}\right)(T_c-T)^2}\rbrace

다음과 같은 결과를 얻게 됨니다.

f_{mf}=-k_b T\ln 2 \ \ \ ; T>T_c

f_{mf}=-k_B\left[T\ln 2-\frac{8}{9}\frac{T}{T_c^2}(T-T_c)^2\right] \ \ ;\ T<T_c

책이랑 다르지만 이게 확실히 맞습니다. 

임계온도보다 온도가 높을때는 확연하게 스핀항이 없는 것이 확실하게 free energy를 낮추는 것이구요. 임계온도보다 온도가 낮을 때는 스핀의 기대값이 (자화 magnetization)이 확실하게 존재 하므로 두번째 세번째 항이 존재하게 되서 온도가 임계온도보다 낮을때 저런 관계가 존재하게 됨니다. 이 결과에서 우리는 비열(specific heat)를 알수 있게 되는데요. 

C_{h=0}=-T\frac{\partial^2f}{\partial T^2} =0 \ \ ; T>T_c

C_{h=0}=-T\frac{\partial^2f}{\partial T^2} =k_B\frac{T}{T_c^2}\left[\frac{9}{2}T-3T_c\right] \ \ ;T<T_c

이 부분이 책에 나와있는 식입니다. 네.. 저 식은 당연히 틀린 것입니다. 실제 계산은 다르겠죠?

C_{h=0}=-T\frac{\partial^2 f}{\partial T^2} =-T\frac{8}{9}\frac{\partial f}{\partial T} \left[k_B\frac{1}{T_c^2}(T-T_c)^2+k_B\frac{T}{T^2}2(T-T_c)\right]

=-T\frac{8}{9} k_B\left[ \frac{2}{T_c^2}(T-T_c)+\frac{2}{T_c^2}(T-T_c)+\frac{2T}{T_c^2} \right]

=-T\frac{8}{9} k_B\left[ \frac{4}{T_c^2}(T-T_c)+\frac{2T}{T_c^2} \right]

=-T\frac{8}{9} k_B\left[ \frac{(6T-4T_c)}{T_c^2} \right]

=-k_B\frac{T}{T_c^2} \left[ \frac{16}{3}T-\frac{32}{9}T_c \right]

C_{h=0}=k_B\frac{T}{T_c^2} \left[ \frac{32}{9}T_c-\frac{16}{3}T \right] \ \ \ \ ;T<T_c

완성입니다. 

두번 미분하는게 좀 까다롭긴 하지만 생각 보다 어려운 것은 아닙니다. 결과값을 실험값과 비교해 보도록 하죠.  실험 값에는 

C_h \sim A_+\left(\frac{T-T_c}{T_c}\right)^{-\alpha} \ \ \ ; T>T_c

C_h \sim A_- \left(\frac{T_c-T}{T_c}\right)^{-\alpha} \ \ \ ; T<T_c

여기서 A는 그냥 실험적인 상수구요. 비교를 해보니까 비슷한거 같으면서 어떤부분은 많이 다르네요. 평균장 이론의 식만 본다면 비열은 상전이 점에서 연속적이지 않습니다. 우리가 대부분의 온도인 임계온도 위에서만이라도 맞춰본다면 \alpha=0 이라고 이야기 할수 있습니다. 실제값인 \alpha=0.11 과는 조금.. 머네요. 확실히 0과 그보다 조금이라도 큰값은 많이 다른 의미를 가지고 있는 것이니까요.

 

2.2절은 여기 까지네요. 생각 보다 책이 틀린 부분이 많네요. 이제 좀더 어려운 2.3절을 공부해보도록 합니다. !! 힘내시길~!

이 글은 스프링노트에서 작성되었습니다.

Posted by blindfish