Hàm Brillouin là gì?
Hàm Brillouin \(B_J(x)\) mô tả độ từ hóa của một chất thuận từ cấu tạo từ các nguyên tử có số lượng tử mômen động lượng toàn phần J. Trong cơ học thống kê, biến không thứ nguyên được tính bằng \(x = g\cdot\mu_B\cdot J\cdot B / (k_B\cdot T)\), tức tỉ số giữa năng lượng từ và năng lượng nhiệt. Công cụ này là một bộ tính hàm đặc biệt thuần toán học: bạn nhập trực tiếp giá trị x (không cần đơn vị vật lý), và máy sẽ trả về bảng giá trị cùng đồ thị của \(B_J(x)\). Khi J trở nên rất lớn, hàm này tiến dần tới hàm Langevin cổ điển \(L(x)\) — bạn có thể chọn chế độ này bằng cách nhập "inf".
Cách sử dụng máy tính
Nhập J dưới dạng số nguyên, phân số bán nguyên như 1/2 hoặc 3/2, hoặc số thập phân như 0.5. Gõ "inf" để chuyển sang giới hạn Langevin. Sau đó đặt giá trị x đầu tiên (Giá trị ban đầu của x), khoảng cách giữa các điểm (Bước nhảy) và số dòng cần tạo. Các giá trị x được sinh ra theo công thức $$x_i = \text{startX} + i\cdot\text{stepX}$$ với i chạy từ 0 đến \(\text{count}-1\). Công cụ sẽ in ra từng cặp \((x, B_J(x))\) và vẽ đường cong tương ứng.
Giải thích công thức
Với J hữu hạn, hàm này kết hợp hai hàm cotang hyperbolic: $$B_J(x) = \frac{2J+1}{2J}\coth\!\left(\frac{2J+1}{2J}\,x\right) - \frac{1}{2J}\coth\!\left(\frac{x}{2J}\right)$$ Đây là một hàm lẻ, nên \(B_J(-x) = -B_J(x)\), đi qua gốc tọa độ (\(B_J(0)=0\)) và bão hòa về \(\pm 1\) khi \(x \to \pm\infty\). Vì coth có điểm kỳ dị tại 0 nên máy tính trả về 0 khi x bằng (hoặc cực kỳ gần) gốc tọa độ — đây chính là giới hạn giải tích đúng.
Ví dụ minh họa
Lấy \(J = 1/2\) (tức \(2J = 1\)) tại \(x = 1\). Khi đó \(\frac{2J+1}{2J} = 2\) và \(\frac{1}{2J} = 1\), ta được $$B_{1/2}(1) = 2\cdot\coth(2) - \coth(1) = 2(1.037314) - 1.313035 = 0.761594$$ Để kiểm chứng: với \(J = 1/2\), hàm Brillouin chính bằng \(\tanh(x)\), và \(\tanh(1) = 0.761594\). Còn với trường hợp Langevin tại \(x = 2\): $$L(2) = \coth(2) - \frac{1}{2} = 1.037314 - 0.5 = 0.537314$$
Giải thích Kết quả Hàm Brillouin
Giá trị trả về bởi máy tính, \(B_J(x)\), không có đơn vị và bị giới hạn giữa 0 và 1. Về mặt vật lý, nó bằng từ hóa phân số của một chất thuận từ — tỷ lệ giữa từ hóa thực tế và từ hóa bão hòa:
$$B_J(x) = \frac{M}{M_\text{bão hòa}}, \qquad 0 \le B_J(x) \le 1.$$Kết quả bằng 0 có nghĩa là không có căn chỉnh ròng của các momen từ (trường không hoặc nhiệt độ vô hạn), trong khi kết quả tiến gần 1 có nghĩa là mọi momen được căn chỉnh hoàn toàn theo trường tác dụng (bão hòa hoàn toàn).
Đối số x: năng lượng từ và năng lượng nhiệt
Đầu vào \(x\) là tỷ lệ giữa năng lượng từ (Zeeman) của một momen và năng lượng nhiệt khả dụng:
$$x = \frac{g\,\mu_B\,J\,B}{k_B\,T}.$$Khi \(x\) nhỏ, sự rung động nhiệt ngẫu nhiên \(k_B T\) chiếm ưu thế so với năng lượng từ căn chỉnh, do đó các momen gần như bị ngẫu nhiên hóa; khi \(x\) lớn, năng lượng từ thắng và các momen bị khóa vào vị trí căn chỉnh.
Chế độ x nhỏ (Curie)
Đối với \(x \ll 1\) hàm Brillouin tuyến tính trong \(x\):
$$B_J(x) \approx \frac{J+1}{3J}\,x.$$Thay \(x = g\mu_B J B /(k_B T)\) cho một từ hóa tỷ lệ với \(B/T\), đây chính là định luật Curie: độ cảm ứng từ giảm theo \(1/T\). Đây là chế độ áp dụng cho các chất thuận từ thông thường trong trường từ phòng thí nghiệm ở nhiệt độ phòng, trong đó \(B_J(x)\) thường nằm xa dưới 1.
Chế độ bão hòa x lớn
Đối với \(x \gg 1\) cả hai hàm côtang hyperbolic đều có xu hướng bằng 1 và hàm bão hòa:
$$B_J(x) \to 1.$$Điều này tương ứng với các trường mạnh và/hoặc nhiệt độ rất thấp, trong đó về cơ bản tất cả các momen từ chỉ theo trường và từ hóa không thể tăng lên được nữa. Trên biểu đồ, điều này xuất hiện dưới dạng một cao nguyên tiến gần đến đường nằm ngang \(B_J=1\). Khi \(J \to \infty\) đường cong tiến gần đến hàm Langevin cổ điển \(L(x)=\coth x - 1/x\).
Các Thuật ngữ và Biến Chính
| Ký hiệu / Thuật ngữ | Ý nghĩa |
|---|---|
| \(J\) | Số lượng tử momen góc toàn phần của ion từ (kết hợp các đóng góp quỹ đạo và spin). Có thể là số nguyên hoặc nửa nguyên (ví dụ: 1/2, 1, 3/2, 2). Nó xác định hình dạng của đường cong và số lượng trạng thái \(m_J\) có thể truy cập được, \(2J+1\). |
| \(x\) | Đối số không có đơn vị của hàm, \(x = g\mu_B J B/(k_B T)\) — tỷ lệ giữa năng lượng từ (Zeeman) và năng lượng nhiệt. Đây là trục ngang của bảng và biểu đồ. |
| \(g\) | Hệ số g của Landé (hệ số tách rã phổ), một số không có đơn vị liên hệ momen từ với momen góc. Đối với spin thuần túy \(g \approx 2\); đối với momen góc quỹ đạo và spin kết hợp, nó được cho bởi công thức Landé. |
| \(\mu_B\) | Magneton Bohr, đơn vị tự nhiên của momen từ nguyên tử, \(\mu_B = e\hbar/(2m_e) \approx 9,274\times10^{-24}\ \text{J/T}\). |
| \(k_B\) | Hằng số Boltzmann, \(k_B \approx 1,381\times10^{-23}\ \text{J/K}\), chuyển đổi nhiệt độ thành năng lượng nhiệt \(k_B T\). |
| \(B\) | Mật độ từ thông (trường từ), được đo bằng tesla (T). \(B\) lớn hơn làm tăng \(x\) và đẩy hệ thống về phía bão hòa. |
| \(T\) | Nhiệt độ tuyệt đối tính bằng kelvin (K). \(T\) cao hơn làm tăng ngẫu nhiên hóa nhiệt, làm giảm \(x\) và từ hóa. |
| \(\coth\) | Côtang hyperbolic, \(\coth(u) = \cosh(u)/\sinh(u) = (e^{u}+e^{-u})/(e^{u}-e^{-u})\); nó xuất hiện hai lần trong hàm Brillouin và có xu hướng bằng 1 đối với \(u\) lớn. |
| Hàm Langevin \(L(x)\) | Giới hạn cổ điển của hàm Brillouin khi \(J \to \infty\): \(L(x) = \coth x - 1/x\). Nó mô tả các dipôl từ cổ điển quay tự do (không có lượng tử hóa của hướng). |
Câu hỏi thường gặp
Tại sao \(B_J(0)\) lại bằng 0? Cả hai số hạng coth đều phân kỳ tại \(x = 0\), nhưng hiệu của chúng có giới hạn hữu hạn bằng 0; công cụ sẽ báo giá trị giới hạn đó.
Những giá trị J nào hợp lệ? Các số nguyên dương và bán nguyên (1/2, 1, 3/2, 2, ...). \(J = 0\) không hợp lệ vì sẽ gây chia cho 0, và máy tính sẽ cảnh báo khi bạn nhập giá trị không phải bán nguyên.
Làm thế nào để dùng hàm Langevin? Nhập "inf" (hoặc "infinity") cho J để sử dụng \(L(x) = \coth(x) - \frac{1}{x}\).
