Giới thiệu dự án

Nghiên cứu cấu trúc điện tử và phổ năng lượng của nguyên tử hydro ($H$) đặt trong các điều kiện môi trường cực hạn là nền tảng của vật lý nguyên tử, vật lý plasma và vật lý thiên văn hiện đại. Trong các môi trường vật lý khắc nghiệt như từ quyển của sao lùn trắng ($B \sim 10^2 - 10^4\text{ T}$), sao neutron ($B \sim 10^7 - 10^9\text{ T}$) hay buồng giam giữ từ tính của lò phản ứng nhiệt hạch Tokamak (như KSTAR đạt mốc $100\text{ triệu }^\circ\text{C}$ duy trì $48\text{ giây}$), các tương tác tĩnh điện Coulomb đơn thuần bị biến tính mạnh mẽ bởi hai yếu tố: hiệu ứng màn chắn điện tích Debye-Yukawa và trường thế nhiệt - từ (thermal-magnetic field).

                      +------------------------------------------+
                      |       Nguyên tử Hydro trong Plasma       |
                      +------------------------------------------+
                                           |
                  +------------------------+------------------------+
                  |                                                 |
                  v                                                 v
   +------------------------------+                  +------------------------------+
   |   Thế Màn Chắn Debye-Yukawa  |                  |     Từ Trường & Nhiệt Độ     |
   |   V(r) = -(e^2/r) exp(-λr)   |                  |  B = 0.0 - 1.0 a.u. (10^5 T) |
   |   λ = 1/D (0.0 - 0.1 a.u.)   |                  |  T = 0 - 10^6 K (K = √(3MkT))|
   +------------------------------+                  +------------------------------+
                  |                                                 |
                  +------------------------+------------------------+
                                           |
                                           v
                      +------------------------------------------+
                      |   Toán tử Giả Động Lượng Ps & Parabolic  |
                      |   H_rel Ψ(ξ,η,φ) = E Ψ(ξ,η,φ)            |
                      +------------------------------------------+
                                           |
                                           v
                      +------------------------------------------+
                      |   Phổ Năng Lượng 1s, 2s, 2p0, 2p±1       |
                      +------------------------------------------+

Vấn đề cốt lõi xuất phát từ việc khi đặt nguyên tử hydro trong từ trường ngoài, động lượng của khối tâm $\mathbf{P}$ không còn là một tích phân chuyển động bảo toàn. Sự ghép cặp phi tuyến giữa chuyển động khối tâm và chuyển động tương đối của cặp electron - proton dẫn đến việc phá vỡ tính đối xứng cầu, làm mất hiệu lực của các phép giải giải tích thông thường trong hệ tọa độ cầu. Hơn nữa, sự hiện diện của plasma làm xuất hiện thế màn chắn Yukawa $V(r) = -\frac{e^2}{r}e^{-\lambda r}$ (với $\lambda = 1/D$ là tham số màn chắn tỉ lệ nghịch với bán kính Debye $D$), làm biến đổi sâu sắc các mức năng lượng liên kết và độ dịch chuyển vạch quang phổ.

Mục tiêu cụ thể của đề tài bao gồm:

  1. Xây dựng toán tử giả động lượng $\hat{\mathbf{P}}_s = \mathbf{P} - \frac{e}{2}\mathbf{B}\times\mathbf{r}$ nhằm tách biến chuyển động khối tâm và chuyển động tương đối trong không gian tọa độ parabolic $(\xi, \eta, \phi)$.
  2. Thiết lập biểu thức giải tích tường minh cho các yếu tố ma trận của động năng $\hat{T}$, thế năng màn chắn $\hat{V}\lambda$ và tương tác từ trường đối với các trạng thái kích thích bậc nhất ($n = 2: 2s, 2p_0, 2p{\pm 1}$) và trạng thái cơ bản ($1s$).
  3. Ứng dụng lý thuyết nhiễu loạn dừng không suy biến (Non-degenerate stationary perturbation theory) kết hợp phương pháp biến phân tham số để tính toán phổ năng lượng chính xác đến sai số dưới $0.05%$.
  4. Định lượng hóa ảnh hưởng của hiệu ứng nhiệt độ $T$ thông qua động lượng khối tâm $K = \sqrt{3 M k_B T}$ lên độ dịch chuyển phổ năng lượng trong dải từ trường $\gamma = 0 - 1.0\text{ a.u.}$ và tham số màn chắn $\lambda = 0.0 - 0.1\text{ a.u.}$

Phạm vi nghiên cứu tập trung vào nguyên tử hydro đơn lẻ trong plasma cân bằng tĩnh nhiệt, xét trong hệ quy chiếu phi tương đối tính với các mức kích thích $n=1, 2$.


Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Nghiên cứu nguyên tử hydro trong trường ngoài đã trải qua nhiều giai đoạn phát triển với các mô hình tính toán khác nhau.

Tiêu chí so sánh Phương pháp Tọa độ cầu giải tích Phương pháp Hartree-Fock 2D số trị Giải pháp đề xuất: Tọa độ Parabolic & Toán tử Giả động lượng
Tính đối xứng Giả định đối xứng cầu đối xứng hoàn toàn Lưới số 2D rời rạc Khai thác trực tiếp đối xứng trục của từ trường dọc trục $Oz$
Xử lý khối tâm Bỏ qua ghép cặp chuyển động khối tâm ($\infty$ khối lượng hạt nhân) Xử lý số học gần đúng, chi phí tính toán cao $O(N^3)$ Tách biến chính xác qua véc-tơ giả động lượng $\hat{\mathbf{P}}_s$, bảo toàn đại số giao hoán
Thế màn chắn Yukawa Khó tích phân giải tích hàm sóng cầu Xấp xỉ tích phân số dễ bất ổn định tại $r \to 0$ Khai triển giải tích qua hàm siêu bội hợp lưu $F(-n,
Hiệu ứng nhiệt độ Không tích hợp tham số nhiệt độ $T$ Phải mô phỏng động lực học phân tử bổ sung Tích hợp trực tiếp vector nhiệt $K = \sqrt{3Mk_BT}$ vào phần tử nhiễu loạn

Yêu cầu tính toán được lượng hóa theo chuẩn MoSCoW:

  • Must-have: Tách biến Hamiltonian bằng biến đổi uniter $U = \exp\left[\frac{i}{\hbar}\left(\mathbf{K} - \frac{e}{2}\mathbf{B}\times\mathbf{r}\right)\cdot\mathbf{R}\right]$; công thức giải tích ma trận thế Yukawa bậc $n=1, 2$.
  • Should-have: Khảo sát liên tục dải thông số $\gamma \in [0, 1.0]$ và $\lambda \in [0, 0.1]$; đối chuẩn dữ liệu với công trình chuẩn quốc tế.
  • Could-have: Mở rộng cho các trạng thái kích thích bậc cao hơn ($n=3, 4$).
  • Won't-have: Tính toán hiệu ứng điện động lực học lượng tử (QED) và cấu trúc siêu tinh tế trong pha này.

Thiết kế hệ thống

Kiến trúc giải thuật tính toán bao gồm 4 khối chức năng chính:

+-----------------------------------------------------------------------------------+
|                            LUỒNG XỬ LÝ TOÁN TỬ VÀ MA TRẬN                         |
+-----------------------------------------------------------------------------------+
| 1. Khối Chuyển đổi Tọa độ:                                                        |
|    (x, y, z) ---> (ξ = r + z, η = r - z, φ = arctan(y/x))                         |
|    dV = 1/4 (ξ + η) dξ dη dφ                                                      |
+-----------------------------------------------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
| 2. Khối Chuẩn hóa Hàm sóng:                                                       |
|    Ψ_n1,n2,m(ξ,η,φ) = C * exp(-α(ξ+η)/2) * (ξη)^(|m|/2) * L_n1^|m|(αξ) *          |
|                       L_n2^|m|(αη) * exp(imφ)                                     |
+-----------------------------------------------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
| 3. Khối Tích phân Yếu tố Ma trận:                                                 |
|    <n1',n2',m'| T + V_Yukawa + H_B + H_T | n1,n2,m>                               |
+-----------------------------------------------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
| 4. Khối Tối ưu hóa Biến phân & Nhiễu loạn:                                        |
|    ∂E^(0)(α)/∂α = 0  ===> α_opt  ===> E = E^(0)(α_opt) + ΔE^(1) + ΔE^(2)          |
+-----------------------------------------------------------------------------------+

Hệ thống tính toán khoa học được triển khai trên nền tảng:

  • Core Computational Kernel: GNU Fortran 13.2 / Python 3.11.8
  • Algebraic Engine: SymPy 1.13.1 (tích phân giải tích tham số)
  • Numerical Library: NumPy 1.26.4, SciPy 1.12.0 (tối ưu hóa phi tuyến bậc cao)
  • Visualization: Matplotlib 3.8.3, Seaborn 0.13.2

Methodology

Phương pháp luận nghiên cứu kết hợp giữa phương pháp giải tích giải toán tử và phương pháp mô phỏng số chuẩn hóa:

  • Pha 1 (Phân tích lý thuyết): Xây dựng đại số giao hoán giữa toán tử giả động lượng $\hat{\mathbf{P}}_s$ và Hamiltonian tương tác $\hat{H}$.
  • Pha 2 (Giải tích hóa yếu tố ma trận): Áp dụng hệ tọa độ parabolic biến đổi vi phân Laplace $\nabla^2 = \frac{4}{\xi+\eta}\left[\frac{\partial}{\partial\xi}\left(\xi\frac{\partial}{\partial\xi}\right) + \frac{\partial}{\partial\eta}\left(\eta\frac{\partial}{\partial\eta}\right)\right] + \frac{1}{\xi\eta}\frac{\partial^2}{\partial\phi^2}$.
  • Pha 3 (Thẩm định số): Kiểm soát rủi ro phân kỳ bằng thuật toán tìm nghiệm biến phân $\alpha$ đảm bảo đạo hàm năng lượng triệt tiêu $\frac{\partial E}{\partial \alpha} = 0$.

Implementation và kết quả

Development process

Quá trình cài đặt thuật toán được phân chia thành các module toán tử độc lập. Dưới đây là đoạn mã nguồn thực thi thuật toán giải tích ma trận Yukawa và tối ưu hóa tham số biến phân $\alpha$ viết bằng Python:

import numpy as np
from scipy.optimize import root_scalar
import sympy as sp

class ScreenedHydrogenSolver:
    """
    Module tinh toan pho nang luong nguyen tu Hydro trong tu truong 
    co the man chan Yukawa su dung he toa do Parabolic va ly thuyet nhieu loan.
    """
    def __init__(self, gamma: float, lam: float, temp_k: float = 0.0):
        self.gamma = gamma          # Tu truong quy chuan (a.u.)
        self.lam = lam              # Tham so man chan Debye (a.u.)
        self.temp_k = temp_k        # Nhiet do plasma (Kelvin)
        self.k_b = 3.166811563e-6   # Hang so Boltzmann (Hartree/K)
        self.M = 1836.152673 + 1.0  # Khoi luong khoi tam (a.u.)
        
    def energy_zero_order_1s(self, alpha: float) -> float:
        """
        Nang luong gan dung bac 0 muc 1s: E^(0)(alpha)
        E^(0) = alpha^2 / 2 - 4 * alpha^3 / (2*alpha + lam)^2 + (gamma^2 / (2 * alpha^2))
        """
        v_screened = 4.0 * (alpha**3) / ((2.0 * alpha + self.lam)**2)
        h_kinetic = (alpha**2) / 2.0
        h_magnetic = (self.gamma**2) / (2.0 * (alpha**2))
        return h_kinetic - v_screened + h_magnetic

    def optimize_alpha_1s(self) -> float:
        """Tim tham so bien phan alpha thoa man dE^(0)/dalpha = 0"""
        def d_energy(a: float) -> float:
            d_t = a
            d_v = (12.0 * (a**2) * (2.0 * a + self.lam) - 16.0 * (a**3)) / ((2.0 * a + self.lam)**3)
            d_mag = - (self.gamma**2) / (a**3)
            return d_t - d_v + d_mag
        
        res = root_scalar(d_energy, bracket=[0.1, 3.0], method='brentq')
        return res.root

    def compute_ground_state(self) -> dict:
        alpha_opt = self.optimize_alpha_1s()
        e_total = self.energy_zero_order_1s(alpha_opt)
        
        # Bo chinh nhiet - tu: K = sqrt(3 * M * k_b * T)
        k_thermal = np.sqrt(3.0 * self.M * self.k_b * self.temp_k)
        delta_e_thermal = (self.gamma * k_thermal) / (self.M * alpha_opt) if self.temp_k > 0 else 0.0
        
        return {
            "state": "1s",
            "alpha_opt": alpha_opt,
            "energy_hartree": e_total + delta_e_thermal,
            "binding_energy": (self.gamma / 2.0) - (e_total + delta_e_thermal)
        }

# Khoi chay mo phong cho trang thai co ban va kich thich
solver = ScreenedHydrogenSolver(gamma=0.01, lam=0.02, temp_k=10000.0)
result = solver.compute_ground_state()
print(f"Trang thai 1s: E = {result['energy_hartree']:.6f} a.u., alpha = {result['alpha_opt']:.4f}")

Testing và validation

Kết quả tính toán của luận văn được thẩm định thông qua việc đối chuẩn trực tiếp với các công trình chuẩn mực quốc tế:

  1. Nghiên cứu của Montgomery et al. [1] (khảo sát theo tham số màn chắn $\lambda$ tại $\gamma = 0$).
  2. Nghiên cứu của Schmelcher & Cederbaum [48] (khảo sát theo từ trường $\gamma$ tại $\lambda = 0$).
       Mức Năng Lượng 1s (Hartree) theo Tham Số Màn Chắn λ (tại γ = 0)
  -0.30 |                                                  
        |                                       * Đề tài
  -0.35 |                                       o Montgomery [1]
        |                             * o
  -0.40 |                   * o
        |         * o
  -0.45 |   * o
        |
  -0.50 +--------------------------------------------------> λ (a.u.)
       0.00      0.02      0.04      0.06      0.08      0.10
Mức trạng thái Tham số $\lambda$ Từ trường $\gamma$ Kết quả Đề tài (a.u.) Đối chuẩn Văn liệu [1], [48] Độ lệch tương đối (%)
$1s$ $0.00$ $0.000$ $-0.500000$ $-0.500000$ [1] $0.000%$
$1s$ $0.02$ $0.000$ $-0.480198$ $-0.480200$ [1] $0.0004%$
$1s$ $0.10$ $0.000$ $-0.409091$ $-0.409100$ [1] $0.0022%$
$2s$ $0.00$ $0.000$ $-0.125000$ $-0.125000$ [1] $0.000%$
$2s$ $0.05$ $0.000$ $-0.089415$ $-0.089410$ [1] $0.0056%$
$2p_0$ $0.05$ $0.000$ $-0.089415$ $-0.089410$ [1] $0.0056%$
$1s$ $0.00$ $0.005$ $-0.499975$ $-0.499974$ [48] $0.0002%$
$1s$ $0.00$ $0.010$ $-0.499900$ $-0.499899$ [48] $0.0002%$
$2s$ $0.00$ $0.010$ $-0.125398$ $-0.125395$ [48] $0.0024%$

Kết quả đạt được

  • Hoàn thành $100%$ các mục tiêu tính toán giải tích cho 4 mức trạng thái: $1s, 2s, 2p_0, 2p_{\pm 1}$.
  • Xác lập được quy luật suy biến bậc nhất: khi $\gamma = 0$, mức $2s$ và $2p_0$ hoàn toàn trùng khớp năng lượng ở mọi giá trị $\lambda$ màn chắn.
  • Khi có sự xuất hiện đồng thời của từ trường ($\gamma > 0$) và nhiệt độ ($T > 0$), mức năng lượng kích thích $2s$ dịch chuyển rõ nét (lên đến $5.8%$ tại $T = 10^5\text{ K}$, $\gamma = 0.1\text{ a.u.}$), khẳng định vai trò chi phối của hiệu ứng nhiệt - từ.

Đổi mới và đóng góp

  1. Ứng dụng toán tử giả động lượng toàn phần trong tọa độ Parabolic: Giải quyết dứt điểm sự bất đối xứng của chuyển động khối tâm trong từ trường đều mà không cần giả định hạt nhân đứng yên vô hạn.
  2. Khai triển giải tích ma trận thế Yukawa đa chiều: Thiết lập công thức tường minh dạng tổng hữu hạn chứa đa thức tích phân suy rộng: $$D = \sum_{k=0}^{n_1} \sum_{s=0}^{n_2} \frac{[(n_1)_k (n_2)_s]^2}{(|m|+k)! k! (|m|+s)! s!} \frac{2\alpha (|m|+k+s)!}{(2\alpha + \lambda)^{|m|+k+s+1}}$$ cho phép triệt tiêu hoàn toàn sai số tích phân số tại vùng kỳ dị lõi ($r \to 0$).
  3. Mô hình hóa định lượng hiệu ứng ghép cặp nhiệt - từ: Đưa số hạng kích thích nhiệt $K = \sqrt{3 M k_B T}$ trực tiếp vào phương trình trị riêng, chứng minh sự phá vỡ tính bảo toàn phổ quang học tĩnh truyền thống trong môi trường plasma thiên văn.

Ứng dụng thực tế và triển khai

+-----------------------------------------------------------------------------+
|                      KHÔNG GIAN ỨNG DỤNG THỰC TIỄN                          |
+-----------------------------------------------------------------------------+
|                                                                             |
|   [1] VẬT LÝ THIÊN VĂN & VŨ TRỤ                                             |
|   - Giải mã quang phổ phát xạ sao lùn trắng (B ~ 10^2 - 10^4 T).            |
|   - Định chuẩn mô hình khí quyển sao neutron mật độ cao.                    |
|                                                                             |
|   [2] PHẢN ỨNG NHIỆT HẠCH TOKAMAK (KSTAR / ITER)                            |
|   - Đo đạc mật độ electron plasma qua độ dịch vạch Zeeman-Yukawa.           |
|   - Đánh giá thông lượng tổn hao năng lượng do bức xạ nguyên tử tạp chất.   |
|                                                                             |
|   [3] VẬT LIỆU BÁN DẪN & CHẤM LƯỢNG TỬ (QUANTUM DOTS)                      |
|   - Mô phỏng trạng thái Exciton bị giam giữ trong điện từ trường.           |
|   - Tối ưu hóa linh kiện quang điện tử thế hệ mới.                          |
|                                                                             |
+-----------------------------------------------------------------------------+

Lộ trình ứng dụng được triển khai qua 3 giai đoạn:

  • Giai đoạn 1 (Hiện tại): Tích hợp module vào thư viện phân tích quang phổ Plasma Diagnostic Suite.
  • Giai đoạn 2 (6-12 tháng): Mở rộng tính toán cho các ion dạng hydro ($He^+, Li^{2+}, C^{5+}$) trong lò Tokamak.
  • Giai đoạn 3 (12-24 tháng): Xây dựng cơ sở dữ liệu mở (Open Atomic Database) phục vụ các đài thiên văn không gian.

Hạn chế và hướng phát triển

  • Hạn chế: Mô hình dừng lại ở lý thuyết nhiễu loạn bậc 1 và biến phân bậc 0; chưa xét hiệu ứng tự quay spin-quỹ đạo (Spin-Orbit Coupling) và hiệu ứng tương đối tính Dirac khi từ trường đạt ngưỡng tới hạn $B \gg 10^9\text{ T}$.
  • Hướng phát triển: Mở rộng cho hệ nguyên tử nhiều electron ($He, Li$) trong thế chắn nhiều tâm; áp dụng mạng nơ-ron vật lý (PINNs - Physics-Informed Neural Networks) để giải phương trình Schrödinger đa chiều dưới trường bức xạ biến thiên thời gian.

Đối tượng hưởng lợi

  • Sinh viên & Học viên cao học: Tiếp cận tài liệu chuẩn mực về phương pháp tách biến giả động lượng và cơ chế lượng tử trong hệ tọa độ parabolic.
  • Kỹ sư vận hành Plasma / Fusion: Có công cụ giải tích chuẩn xác để hiệu chuẩn phổ đo đạc nhiệt độ và mật độ Debye trong buồng phản ứng.
  • Nhà nghiên cứu Vật lý Thiên văn: Nhận bộ dữ liệu năng lượng trạng thái kích thích $n=2$ độ chính xác cao để phân tích quang phổ sao đặc.

Câu hỏi thường gặp

1. Yêu cầu cấu hình tính toán và triển khai hệ thống là gì?

Hệ thống sử dụng các công thức giải tích đóng nên yêu cầu tài nguyên cực kỳ tối ưu: CPU lõi đơn $\ge 2.0\text{ GHz}$, RAM $2\text{ GB}$, hỗ trợ nền tảng Python 3.10+ hoặc trình biên dịch Fortran 90 trở lên. Thời gian giải một trạng thái mất dưới $5\text{ ms}$.

2. Giới hạn tham số màn chắn $\lambda$ và từ trường $\gamma$ là bao nhiêu?

Mô hình nhiễu loạn đạt độ chính xác cao nhất trong vùng $\lambda \in [0.0, 0.1\text{ a.u.}]$ (tương ứng bán kính Debye $D \ge 10 a_0$) và từ trường yếu - trung bình $\gamma \in [0.0, 0.2\text{ a.u.}]$ ($B \le 4.7 \times 10^4\text{ T}$). Vượt quá ngưỡng này, cần chuyển đổi sang phương pháp biến phân lưới số cơ sở B-spline.

3. Làm thế nào để tích hợp thuật toán vào phần mềm phân tích quang phổ hiện có?

Mã nguồn được đóng gói dạng Python Module chuẩn PEP-8 hoặc C-Shared Library (.so / .dll), dễ dàng liên kết với các hệ thống phân tích dữ liệu như ROOT CERN hay OriginLab qua API wrapper.

4. Chi phí tính toán so với các phương pháp giải số lưới không gian 3D?

Phương pháp giải tích biến phân trong đề tài giúp giảm hơn $99.4%$ thời gian CPU so với phương pháp sai phân hữu hạn 3D (FDM) hoặc phần tử hữu hạn (FEM), từ hàng chục giây xuống dưới vài mili-giây cho mỗi điểm phổ.

5. Khả năng bảo trì và mở rộng cho các hệ nguyên tử khác như thế nào?

Cấu trúc module hóa cao cho phép thay đổi điện tích hạt nhân $Z$ để mô phỏng tức thì các ion tương tự hydro ($He^+, Li^{2+}$) chỉ bằng cách tái định nghĩa hàm thế $V(r) = -Z \frac{e^2}{r} e^{-\lambda r}$.


Kết luận

Khóa luận tốt nghiệp của tác giả Huỳnh Gia Phú (hướng dẫn bởi TS. Lý Duy Nhất, Trường ĐH Sư phạm TP.HCM) đã giải quyết thành công bài toán phức tạp về phổ năng lượng nguyên tử hydro dưới tác động đồng thời của từ trường mạnh, thế màn chắn plasma Yukawa và hiệu ứng nhiệt độ khối tâm. Bằng việc kết hợp tinh tế toán tử giả động lượng trong hệ tọa độ parabolic với lý thuyết nhiễu loạn biến phân, đề tài không chỉ đạt độ chuẩn xác cao tương đương các công bố quốc tế mà còn mở ra phương pháp tiếp cận giải tích hiệu quả, phục vụ trực tiếp cho chẩn đoán plasma nhiệt hạch và nghiên cứu vật lý thiên văn hiện đại.