브릴루앙 함수란?
브릴루앙 함수 \(B_J(x)\)는 총 각운동량 양자수가 J인 원자들로 이루어진 상자성체의 자화(磁化)를 기술하는 함수입니다. 통계역학에서 무차원 변수는 \(x = g\cdot\mu_B\cdot J\cdot B / (k_B\cdot T)\)로 정의되며, 이는 자기 에너지와 열 에너지의 비를 나타냅니다. 이 계산기는 물리 단위를 다루지 않는 순수 수학적 특수함수 계산기입니다. 즉, x 값을 직접 입력하면 \(B_J(x)\)의 표와 그래프를 돌려줍니다. J가 매우 커지면 이 함수는 고전적인 랑주뱅 함수 \(L(x)\)에 수렴하는데, J 값에 "inf"를 입력해 이 경우를 선택할 수 있습니다.
계산기 사용법
J는 정수, 1/2이나 3/2 같은 반정수 분수, 또는 0.5 같은 소수로 입력할 수 있습니다. 랑주뱅 극한을 원하면 "inf"라고 입력하세요. 그런 다음 첫 번째 x 값(x의 초깃값), 점 사이의 간격(증분), 그리고 생성할 행의 개수를 설정합니다. x 값은 \(x_i = \text{startX} + i\cdot\text{stepX}\) (i = 0부터 count\(-\)1까지) 규칙으로 만들어집니다. 계산기는 모든 \((x, B_J(x))\) 쌍을 출력하고 곡선을 그려 줍니다.
공식 설명
유한한 J에 대해 이 함수는 두 개의 쌍곡코탄젠트(coth)를 결합한 형태입니다. $$B_J(x) = \frac{2\text{J}+1}{2\text{J}}\coth\!\left(\frac{2\text{J}+1}{2\text{J}}\,x\right) - \frac{1}{2\text{J}}\coth\!\left(\frac{x}{2\text{J}}\right)$$ 이 함수는 기함수이므로 \(B_J(-x) = -B_J(x)\)가 성립하고, 원점을 지나며(\(B_J(0)=0\)), \(x \to \pm\infty\)일 때 \(\pm 1\)로 포화됩니다. coth는 0에서 발산하기 때문에, 계산기는 x가 원점이거나 원점에 극도로 가까울 때 0을 반환하며, 이것이 올바른 해석적 극한값입니다.
계산 예시
J = 1/2 (즉 2J = 1)이고 x = 1인 경우를 생각해 봅시다. 이때 \((2J+1)/2J = 2\), \(1/2J = 1\)이므로 $$B_{1/2}(1) = 2\cdot\coth(2) - \coth(1) = 2(1.037314) - 1.313035 = 0.761594$$가 됩니다. 검산해 보면, J = 1/2일 때 브릴루앙 함수는 \(\tanh(x)\)와 같고 \(\tanh(1) = 0.761594\)이므로 일치합니다. 랑주뱅의 경우 x = 2에서는 \(L(2) = \coth(2) - 1/2 = 1.037314 - 0.5 = 0.537314\)입니다.
브릴루앙 함수 결과 해석
계산기가 반환하는 값 \(B_J(x)\)는 무차원이며 0과 1 사이에 제한됩니다. 물리적으로 이는 상자성체의 부분 자화 — 실제 자화와 포화 자화의 비율을 나타냅니다:
$$B_J(x) = \frac{M}{M_\text{sat}}, \qquad 0 \le B_J(x) \le 1.$$결과가 0이면 자기 모멘트의 순정렬이 없음을 의미하고(영 자장 또는 무한 온도), 결과가 1에 접근하면 모든 모멘트가 인가된 자장과 완전히 정렬됨을 의미합니다(완전 포화).
인수 x: 자기 에너지 대 열 에너지
입력 \(x\)는 모멘트의 자기(제만) 에너지와 이용 가능한 열 에너지의 비율입니다:
$$x = \frac{g\,\mu_B\,J\,B}{k_B\,T}.$$\(x\)가 작을 때는 무작위 열 진동 \(k_B T\)가 정렬 자기 에너지를 압도하므로 모멘트가 거의 무작위화됩니다. \(x\)가 클 때는 자기 에너지가 우세하여 모멘트가 정렬로 고정됩니다.
저-x(큐리) 영역
\(x \ll 1\)일 때 브릴루앙 함수는 \(x\)에 대해 선형입니다:
$$B_J(x) \approx \frac{J+1}{3J}\,x.$$\(x = g\mu_B J B /(k_B T)\)를 대입하면 \(B/T\)에 비례하는 자화를 얻으며, 이는 정확히 큐리 법칙입니다: 감수율이 \(1/T\)로 감소합니다. 이는 실온의 실험실 자장에서 일반적인 상자성체에 적용되는 영역이며, 여기서 \(B_J(x)\)는 전형적으로 1보다 훨씬 작습니다.
고-x 포화 영역
\(x \gg 1\)일 때 두 쌍곡 코탄젠트가 모두 1로 향하고 함수가 포화됩니다:
$$B_J(x) \to 1.$$이는 강한 자장 및/또는 매우 낮은 온도에 해당하며, 본질적으로 모든 자기 모멘트가 자장을 따라 가리키고 자화가 더 이상 증가할 수 없습니다. 그래프에서 이는 수평선 \(B_J=1\)에 접근하는 평탄부로 나타납니다. \(J \to \infty\)일 때 곡선은 고전적 랑주뱅 함수 \(L(x)=\coth x - 1/x\)에 접근합니다.
주요 용어 및 변수
| 기호 / 용어 | 의미 |
|---|---|
| \(J\) | 자기 이온의 전체 각운동량 양자수(궤도 및 스핀 기여를 결합). 정수 또는 반정수일 수 있습니다(예: 1/2, 1, 3/2, 2). 곡선의 모양과 접근 가능한 \(m_J\) 상태의 개수 \(2J+1\)을 결정합니다. |
| \(x\) | 함수의 무차원 인수, \(x = g\mu_B J B/(k_B T)\) — 자기(제만) 에너지와 열 에너지의 비율입니다. 이는 표와 그래프의 수평축입니다. |
| \(g\) | 랑데 g-인수(분광학적 분열 인수), 자기 모멘트를 각운동량과 연결하는 무차원 수입니다. 순수 스핀의 경우 \(g \approx 2\); 궤도 및 스핀 각운동량이 결합된 경우 랑데 공식으로 주어집니다. |
| \(\mu_B\) | 보어 자성체, 원자 자기 모멘트의 자연 단위, \(\mu_B = e\hbar/(2m_e) \approx 9.274\times10^{-24}\ \text{J/T}\). |
| \(k_B\) | 볼츠만 상수, \(k_B \approx 1.381\times10^{-23}\ \text{J/K}\), 온도를 열 에너지 \(k_B T\)로 변환합니다. |
| \(B\) | 자기 플럭스 밀도(자장), 테슬라(T) 단위로 측정합니다. \(B\)가 클수록 \(x\)를 증가시키고 시스템을 포화로 유도합니다. |
| \(T\) | 켈빈(K) 단위의 절대 온도입니다. \(T\)가 높을수록 열 무작위화가 증가하여 \(x\)와 자화가 감소합니다. |
| \(\coth\) | 쌍곡 코탄젠트, \(\coth(u) = \cosh(u)/\sinh(u) = (e^{u}+e^{-u})/(e^{u}-e^{-u})\); 브릴루앙 함수에 두 번 나타나며 큰 \(u\)에 대해 1로 향합니다. |
| 랑주뱅 함수 \(L(x)\) | \(J \to \infty\)일 때 브릴루앙 함수의 고전적 극한: \(L(x) = \coth x - 1/x\). 자유롭게 회전하는 고전 자기 쌍극자(방향의 양자화 없음)를 설명합니다. |
자주 묻는 질문
\(B_J(0)\)이 왜 0으로 표시되나요? x = 0에서 두 coth 항은 모두 발산하지만 그 차이는 유한한 극한값 0을 갖습니다. 계산기는 바로 이 극한값을 표시합니다.
J에는 어떤 값을 넣을 수 있나요? 양의 정수와 반정수(1/2, 1, 3/2, 2, …)가 유효합니다. J = 0은 0으로 나누는 연산이 되므로 사용할 수 없으며, 반정수가 아닌 값을 입력하면 경고가 표시됩니다.
랑주뱅 함수는 어떻게 구하나요? J 값에 "inf"(또는 "infinity")를 입력하면 \(L(x) = \coth(x) - 1/x\)를 사용합니다.
