キュリーの法則のソースを表示

[[常磁性]]物質においては、 その物質の[[磁化]]は、（ほぼ）かけられた[[磁場]]に正比例する。しかし、もし物質が熱せられていると、この線形性は消失する: 一定の磁場については、磁化は（ほぼ）[[温度]]に反比例する。この事実は'''キュリーの法則'''（{{lang-en|Curie's law}}）にまとめられる:

:<math>\boldsymbol{M} = C \frac{\boldsymbol{B}}{T}</math>

ここで
*<math>\boldsymbol{M}</math> は発生する磁化
*<math>\boldsymbol{B}</math> は磁場（単位は[[テスラ (単位)|テスラ]]）
*<math>T</math> は絶対温度（単位は[[ケルビン]]）
*<math>C</math> は物質固有の[[キュリー定数]]。

また、磁化率<math>\chi</math>を用いて以下のように書くこともできる。
{{Indent|<math>\chi = \frac{C}{T}</math>}}

この関係は[[1895年]]に[[ピエール・キュリー]]により（実験結果が想定されるモデルに適合するように調整されつつ）実験的に発見された。その後、[[ポール・ランジュバン]]が理論的に導出した（以下を参照）。そのため、'''キュリー・ランジュバンの法則'''とも呼ばれる。

この法則は高温または弱い磁場についてのみ成り立つ。以下で導く通り、低温または強磁場のような反対側の極限では磁化は飽和する。

なお、[[強磁性|強磁性体]]や[[反強磁性|反強磁性体]]では、キュリーの法則を拡張した[[キュリー・ワイスの法則]]が（ほぼ）成り立っている。

== 量子統計力学での導出 ==
[[ファイル:magnetization2.jpg|thumb|'''常磁性体の磁化'''<br>ただし横軸は温度の[[逆数]]|right|300px]]

単純な[[常磁性|常磁性体]]の[[数学モデル|モデル]]として、互いに相互作用をしない粒子で構成されている物質を考える。それぞれの粒子は<math>\boldsymbol{\mu}</math>で与えられる[[磁気モーメント]]をもつ。磁場中における磁気モーメントの[[エネルギー]]は以下で与えられる。

{{Indent|<math>E=-\boldsymbol{\mu}\cdot\boldsymbol{B}</math>}}

=== 2状態（スピン-1/2）の粒子 ===

計算を簡単にするために、'''2状態'''の粒子を考える。2状態とは、粒子の磁気モーメントが磁場に対して平行（平行で向きも同じ）と反平行（平行だが向きは逆）のどちらかをとることができるということを意味する。よって磁気モーメントに許される値は<math>\mu</math>または<math>-\mu</math>である。この模型を[[イジング模型]]とよぶ。この場合、各粒子は
{{Indent|<math>E_0 = - \mu B</math>}}
か
{{Indent|<math>E_1 = \mu B</math>}}
のどちらかのエネルギーをもつことができる。

次に、粒子が磁場の方向に向くときの向きやすさを考える。この向きやすさを磁化<math>\mu</math>の[[期待値]]で考える。
{{Indent|<math>\left\langle\mu\right\rangle =  \mu P\left(\mu\right) + (-\mu) P\left(-\mu\right) 
 = {1 \over Z} \left( \mu e^{ \mu B\beta} -  \mu e^{  - \mu B\beta} \right)
 = {2\mu \over Z} \sinh( \mu B\beta), </math>}}
ここで粒子の配向の確率はその[[ボルツマン因子]]によって与えられており、更に確率は[[分配関数]]<math>Z</math>で割ることによって[[正規化]]されている（よって全確率の和は1である）。 
一粒子の分配関数は以下で与えられる。
{{Indent|<math>Z = \sum_{n=0,1} e^{-E_n\beta} = e^{ \mu B\beta} + e^{-\mu B\beta} = 2 \cosh\left(\mu B\beta\right)</math>}}

以上より、この単純な場合には次式の一粒子当たりの磁化の[[期待値]]を得る。
{{Indent|<math>\left\langle\mu\right\rangle  = \mu \tanh\left(\mu B\beta\right)</math>}}

さらに、[[固体]]の磁化の総量は次の式で与えられる。
<blockquote style="border: 1px solid black; padding:10px;">
<math>M = N\left\langle\mu\right\rangle = N \mu \tanh\left(\frac{\mu B}{k T}\right)</math></blockquote>

上記の公式はランジュバンの常磁性方程式 (Langevin paramagnetic equation) などと呼ばれる。[[ピエール・キュリー]]は実験において、比較的高温の場合や低磁場の場合におけるこの法則の近似式を発見した。

このランジュバンの求めた磁化が、<math>T</math>が大きく<math>B</math>が小さいという特別な条件下でどのように表されるかを見てみよう。温度が上昇し、磁場が減少するに従い、[[双曲線関数|<math>\tanh</math>]]の引数が減少してゆく。すなわち<math>\mu B/k T \ll 1</math>となる。<!--this is sometimes called the '''Curie regime'''. (註:日本語では聞いたことがないので排除)-->ここで、<math>|x| \ll 1</math>の場合<math>\tanh x \approx x</math>という近似が成り立つため、磁化<math>\boldsymbol{M}</math>は
<blockquote style="border: 1px solid black; padding:10px 2em;">
<math>\boldsymbol{M}(T\rightarrow\infty)=\frac{N\mu^2}{k} \frac{\boldsymbol{B}}{T}</math></blockquote>
と表すことができる。

以上より、キュリーの法則を証明することができた。なお、キュリー定数<math>C</math>は<math>C= N\mu^2/k</math>である。また反対に、低温や高磁場の状況では磁化<math>M</math>は最大値<math>N\mu</math>に漸近する。これは全ての粒子が完全に磁場の方向へと整列していることを意味している。

=== 一般的な場合 ===

粒子が任意のスピン（たくさんのスピン状態）を持っている場合、磁化率の公式は少し複雑になる。このより一般的な公式とその導出には、[[ブリルアン関数]]を参照すること。スピンが無限に近づくにつれ、磁化の公式は以下の節で導出する古典的な値に近づいてゆく。

== 古典統計力学による導出 ==

常磁性磁子が古典的な自由に回転する磁気モーメントであると考えられる場合には、異なる扱いが適用される。この場合、磁気モーメントの状態は[[球座標]]における角度で表すことができる。またひとつ当たりのエネルギー<math>E</math>は以下で表される。

{{Indent|<math>E =  - \mu B\cos\theta, </math>}}

ここで<math>\theta</math>は磁気モーメントと磁場の間の角度である。なお、ここでは磁場は<math>z</math>軸方向を向いているとする。

対応する分配関数<math>Z</math>は以下で表される。
:<math>\begin{align}
Z &= \int_0^{2\pi} d\phi \int_0^{\pi}d\theta \sin\theta \exp( \mu B\beta \cos\theta) \\
  &= -2 \pi \frac{\exp(\mu B \beta \cos \theta)}{\mu B \beta} \Bigg|_{0}^{\pi} \\
&= \frac{4\pi\sinh( \mu B\beta )}{\mu B\beta }
\end{align}</math>
よって磁化の<math>z</math>成分の期待値は以下となる。（<math>z</math>以外は<math>\phi</math>についての積分より0とみなされるため。）{{Indent|<math>\left\langle\mu_z \right\rangle = \frac{1}{Z} \int_0^{2\pi} d\phi \int_0^{\pi}d\theta \sin\theta \, \mu\cos\theta \, \exp( \mu B\beta \cos\theta) </math>}}
計算を簡単にするために、<math>Z</math>の微分を用いて表すと以下となる。
{{Indent|<math>\left\langle\mu_z\right\rangle = {1 \over Z B} \frac{\partial Z}{\partial \beta}</math>}}
（この計算簡略化のアプローチは前述の量子統計力学のモデルの計算でも用いることができるが、もともと計算が複雑ではないため利用する利点は少ない。）

これを計算することにより、次の式を得る。
{{Indent|<math>\left\langle\mu_z\right\rangle = \mu L(\mu B\beta), </math>}}
ここで関数<math>L</math>は
{{Indent|<math> L(x)= \coth x -{1 \over x}</math>}}
で表される[[ランジュバン関数]]である。

ランジュバン関数は引数<math>x</math>が小さい場合には<math>L(x) \approx x/3</math>と近似され、反対に引数<math>x</math>が大きい場合には1に漸近する。これより、上記の<math>\left\langle\mu_z\right\rangle</math>は引数が小さいときにはキュリーの法則に従った振る舞いをするが、キュリー定数は1/3の大きさとなる。また、引数が大きいときには量子統計力学での導出と同様、最大値<math>N\mu</math>へと漸近する。

==応用==
キュリーの法則は、[[磁気温度計]]の原理として用いられている。磁気温度計を用いると、極低温を測定することができる。

==参照==
*[[キュリー・ワイスの法則]]
*[[常磁性]]
*[[ピエール・キュリー]]

<!--Categories-->
{{DEFAULTSORT:きゆりいのほうそく}}
[[Category:電磁気学]]
[[Category:磁気]]
[[Category:自然科学の法則]]
[[Category:ピエール・キュリー]]
[[Category:物理学のエポニム]]

{{sci-stub}}