Giới thiệu dự án
Nghiên cứu cấu trúc hạt nhân nguyên tử ở các trạng thái kích thích và nhiệt độ hữu hạn ($T \neq 0$) là một trong những trụ cột then chốt của vật lý hạt nhân hiện đại và thiên văn vật lý hạt nhân (Nuclear Astrophysics). Trong vũ trụ học hạt nhân, hơn 3.000 đồng vị đã được phát hiện và hàng ngàn đồng vị giàu neutron (neutron-rich nuclei) chưa khám phá đóng vai trò trung tâm trong quá trình tổng hợp hạt nhân bắt neutron nhanh (r-process) tại các vụ sáp nhập sao neutron hay siêu tân tinh.
Tuy nhiên, việc mô tả chính xác các hạt nhân có cấu hình nucleon lẻ (odd-mass nuclei) nằm xa đường bền (far from the stability line) ở nhiệt độ thấp ($T < 1\text{ MeV}$) luôn là thách thức lớn đối với vật lý lý thuyết:
- Khoảng trống của lý thuyết trường trung bình: Phương pháp trường tự hợp Hartree-Fock (Hartree-Fock mean-field theory - HF) sử dụng thế bán thực nghiệm Woods-Saxon (WS) chỉ mô tả chính xác các hạt nhân đóng vỏ (magic nuclei) hoặc nằm trên đường bền ở trạng thái cơ bản. Đối với hạt nhân lớp vỏ mở (open-shell nuclei) xa đường bền, tương tác thặng dư mà đặc biệt là tương quan kết cặp (nuclear pairing correlations) giữa các nucleon đóng vai trò chi phối nhưng bị bỏ qua trong HF thuần túy.
- Hạn chế chí mạng của lý thuyết BCS nhiệt độ hữu hạn (FTBCS): Lý thuyết Bardeen-Cooper-Schrieffer nhiệt độ hữu hạn (FTBCS) vốn thành công trong các hệ vô hạn (chất siêu dẫn vĩ mô), khi áp dụng cho hệ hữu hạn hạt như hạt nhân nguyên tử lại bộc lộ khiếm khuyết nghiêm trọng:
- Vi phạm định luật bảo toàn số hạt ($\Delta N^2 \neq 0$).
- Bỏ qua dao động nhiệt và dao động số giả hạt (Quasiparticle-number fluctuations - QNF).
- Dự đoán sự sụp đổ đột ngột của khe năng lượng kết cặp tại nhiệt độ tới hạn $T_c \approx 0.57 \Delta(T=0)$, dẫn đến gián đoạn vô cực giả tạo trong nhiệt dung $C(T)$.
- Nghịch lý trong cấu hình nucleon lẻ: Ở nhiệt độ bằng 0, nucleon lẻ chiếm đóng mức Fermi và kích hoạt hiệu ứng khóa mức Pauli (Pauli blocking effect), ngăn cản các cặp Cooper tán xạ vào mức này, làm suy yếu tương quan kết cặp. Nhưng tại nhiệt độ hữu hạn ($T < 1\text{ MeV}$), hiện tượng tái kết cặp (Pairing reentrance effect / Thermally assisted pairing) xuất hiện làm tăng cường khe năng lượng kết cặp lại chưa được các mô hình truyền thống giải thích một cách nhất quán.
+-------------------------------------------------------------------------------+
| BÀI TOÁN HẠT NHÂN HỮU HẠN |
+---------------------------------------+---------------------------------------+
| HẠN CHẾ CỦA MÔ HÌNH TRUYỀN THỐNG | GIẢI PHÁP LỜI GIẢI CHÍNH XÁC (FTEP) |
| • FTBCS: Vi phạm bảo toàn số hạt | • Đại số giả spin SU(2) chéo hóa |
| • Khe pairing sụp đổ giả tạo tại Tc | • Bảo toàn chính xác số hạt N |
| • Pauli blocking bị cố định cứng | • Mô tả làm trơn chuyển pha S-shape |
| • Gián đoạn nhiệt dung vô cực | • Tái kết cặp thực tế ở T < 1 MeV |
+---------------------------------------+---------------------------------------+
Mục tiêu nghiên cứu của đề tài được xác định cụ thể:
- Triển khai mô hình lời giải chính xác bài toán kết cặp tại nhiệt độ hữu hạn (Finite-Temperature Exact Pairing - FTEP) trên nền tảng đại số nhóm giả spin $SU(2)$ và hệ thống kê chính tắc (Canonical ensemble).
- Khảo sát định lượng và giải thích bản chất vi mô của hiện tượng tái kết cặp tại vùng nhiệt độ thấp ($T < 1\text{ MeV}$) và sự tồn tại của khe năng lượng kết cặp ở nhiệt độ cao ($T \leq 4\text{ MeV}$).
- Phân tích các đại lượng nhiệt động vi mô bao gồm: khe năng lượng kết cặp ($\Delta$), nhiệt dung ($C$), entropy ($S$) và số chiếm đóng đơn hạt ($f_j$).
- Đánh giá và tiên đoán quy luật mạnh/yếu của hiện tượng tái kết cặp trên chuỗi các hạt nhân thực xa đường bền: $^{83,85,87}\text{Ge}$ ($Z=32$), $^{105,107,111}\text{Zr}$ ($Z=40$), $^{75,77,79}\text{Ni}$ ($Z=28$) và $^{105,107,109}\text{Y}$ ($Z=39$).
Phạm vi nghiên cứu tập trung vào cấu hình neutron lẻ trong không gian rút gọn (truncated space) gồm 10 mức đơn hạt xung quanh mức Fermi, tương tác thông qua Hamiltonian kết cặp đơn cực với hằng số ghép cặp $G$, được lập trình và giải số bằng ngôn ngữ FORTRAN 90.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Để định vị đóng góp của phương pháp FTEP, bảng phân tích ma trận kỹ thuật đối sánh các phương pháp tiếp cận cấu trúc hạt nhân hiện hành được thiết lập:
| Tiêu chí phân tích |
Trường trung bình Hartree-Fock (HF) |
Lý thuyết FTBCS truyền thống |
Lý thuyết Modified BCS / MBCS |
Lời giải chính xác nhiệt độ hữu hạn (FTEP) |
| Bảo toàn số hạt |
Tuyệt đối ($\Delta N = 0$) |
Không bảo toàn ($\Delta N^2 \gg 0$) |
Gần đúng qua phép chiếu số hạt |
Bảo toàn chính xác ($\Delta N = 0$) |
| Dao động nhiệt (QNF) |
Bỏ qua |
Bỏ qua |
Xét gần đúng bậc nhất |
Tích hợp chính xác qua hàm phân chia |
| Chuyển pha siêu dẫn (SN) |
Không áp dụng |
Gián đoạn bậc 2 đột ngột tại $T_c$ |
Chuyển tiếp tương đối trơn |
Chuyển pha trơn mượt dạng S-shape |
| Hành vi khe $\Delta(T)$ |
$\Delta = 0$ |
Triệt tiêu hoàn toàn tại $T_c \approx 0.57\Delta_0$ |
Tồn tại đuôi nhỏ ở nhiệt độ cao |
Tồn tại liên tục, suy giảm tiệm cận đến $T > 4\text{ MeV}$ |
| Hiệu ứng tái kết cặp |
Không mô tả được |
Không xuất hiện (bị khóa mức cứng) |
Xuất hiện gần đúng |
Mô tả chính xác độ tăng $\Delta$ tại $T < 1\text{ MeV}$ |
| Chi phí tính toán |
Thấp ($O(N^3)$) |
Thấp ($O(N^2)$) |
Trung bình ($O(N^4)$) |
Trung bình - cao (Chéo hóa ma trận Fock) |
Yêu cầu bài toán theo mô hình MoSCoW:
- Must have: Bảo toàn nghiêm ngặt số hạt; biểu diễn Hamiltonian dưới dạng ma trận $SU(2)$ có thể chéo hóa; trích xuất chính xác năng lượng đơn hạt $e_j$ từ thế Woods-Saxon; tính toán hàm phân chia chính tắc $Z(T)$.
- Should have: Tính toán đồng thời khe kết cặp $\Delta(T)$, nhiệt dung $C(T)$, entropy $S(T)$, số chiếm đóng $f_j(T)$ cho cả mô hình thử nghiệm ($N=9, \Omega=10$) và hệ hạt nhân thực ($Z=28, 32, 39, 40$).
- Could have: So sánh tương quan chéo giữa các đồng vị isobaric/isotonic (ví dụ: chuỗi Zr vs Y) để tìm ra quy luật vị trí mức đơn hạt.
- Won't have: Mở rộng ra không gian vỏ đầy đủ không giới hạn (do bùng nổ tổ hợp chiều không gian ma trận Hilbert).
+-----------------------------------------------------------------------------------+
| LUỒNG XỬ LÝ DỮ LIỆU TOÀN CỤC (PIPELINE) |
+-----------------------------------------------------------------------------------+
| [Thế Woods-Saxon: V0, R, a] ---> [Trường tự hợp Hartree-Fock] |
| | |
| v |
| [Trích xuất 10 mức đơn hạt ej quanh Fermi] |
| | |
| v |
| [Xây dựng cơ sở Seniority & Ma trận Hamiltonian SU(2)] |
| | |
| v |
| [Chéo hóa trực tiếp ma trận: Trị riêng Es & Trạng thái |s>] |
| | |
| v |
| [Thống kê chính tắc: Hàm phân chia Z(T) = Sum d_s e^(-Es/T)] |
| | |
| +-------------------+--------------------+--------------------+ |
| v v v v |
| [Khe pairing Delta(T)] [Nhiệt dung C(T)] [Entropy S(T)] [Số chiếm đóng fj] |
+-----------------------------------------------------------------------------------+
Thiết kế hệ thống
Kiến trúc thuật toán giải số bài toán FTEP được cấu trúc thành 4 tầng tính toán độc lập nhưng liên kết chặt chẽ:
- Tầng trường trung bình vi mô: Tính toán phương trình Schrödinger với thế Woods-Saxon:
$$V_{WS}(r) = \frac{-V_0}{1 + \exp\left(\frac{r - R}{a}\right)}, \quad R = r_0 A^{1/3}\text{ fm}$$
Trích xuất phổ năng lượng đơn hạt $e_j$ và độ suy biến $\Omega_j = (2j+1)/2$.
- Tầng đại số giả spin $SU(2)$ và không gian rút gọn:
Toán tử giả spin tại mức $j$ được xác định qua toán tử sinh/hủy nucleon $a_{jm}^\dagger, \tilde{a}{jm}$:
$$L_j^+ = \sum{m > 0} a_{jm}^\dagger a_{j\tilde{m}}^\dagger, \quad L_j^- = (L_j^+)^\dagger, \quad L_j^0 = \frac{1}{2} (\hat{N}j - \Omega_j)$$
Hamiltonian kết cặp có dạng:
$$H = \sum{j} e_j \hat{N}j - G \sum{j, j'} L_j^+ L_{j'}^-$$
Ma trận Hamiltonian được xây dựng trên cơ sở trạng thái $|s_1, \dots, s_\Omega; N\rangle$, trong đó $s_j$ là số hạt không kết cặp (seniority).
- Yếu tố ma trận đường chéo:
$$\langle H \rangle_{diag} = \sum_j e_j N_j - G \sum_j \frac{N_j - s_j}{2} \left(\Omega_j - s_j - \frac{N_j - s_j}{2} + 1\right)$$
- Yếu tố ma trận ngoài đường chéo (chuyển dịch một cặp nucleon từ mức $j'$ sang mức $j$):
$$\langle \dots N_j+2 \dots N_{j'}-2 \dots | H | \dots N_j \dots N_{j'} \dots \rangle = - \frac{G}{4} \sqrt{(N_{j'} - s_{j'})(2\Omega_{j'} - s_{j'} - N_{j'} + 2)(2\Omega_j - s_j - N_j)(N_j - s_j + 2)}$$
- Tầng thống kê chính tắc tại nhiệt độ hữu hạn:
Từ tập hợp trị riêng $E_s$ nhận được sau khi chéo hóa ma trận, hàm phân chia chính tắc được thiết lập:
$$Z(T) = \sum_s 2^{s} \exp\left(-\frac{E_s}{T}\right)$$
Các hàm trạng thái nhiệt động lực học:
- Năng lượng tự do: $F(T) = -T \ln Z(T)$
- Năng lượng toàn phần: $\mathcal{E}(T) = \frac{\sum_s 2^s E_s \exp(-E_s/T)}{Z(T)}$
- Entropy: $S(T) = -\frac{\partial F}{\partial T} = \frac{\mathcal{E}(T) - F(T)}{T}$
- Nhiệt dung: $C(T) = T \frac{\partial S}{\partial T} = \frac{\partial \mathcal{E}}{\partial T}$
- Số chiếm đóng đơn hạt: $f_j(T) = \frac{1}{2\Omega_j Z(T)} \sum_s 2^s \langle s | \hat{N}_j | s \rangle \exp\left(-\frac{E_s}{T}\right)$
- Khe năng lượng kết cặp hữu hạn: $\Delta(T) = G \sqrt{\sum_j \Omega_j f_j(T) [1 - f_j(T)] - \frac{\mathcal{E}_{pair}(T)}{G}}$
Technology Stack áp dụng trong dự án:
- Ngôn ngữ cốt lõi: FORTRAN 90/95 (Tối ưu hóa tính toán ma trận mảng đa chiều).
- Thư viện đại số tuyến tính: LAPACK v3.10.0 / BLAS (Sử dụng module
DSYEV để chéo hóa ma trận đối xứng thực).
- Môi trường biên dịch & thực thi: GNU Fortran Compiler (
gfortran v11.4.0) trên nền tảng Linux OS (x86_64).
- Công cụ trực quan hóa & phân tích số liệu: Python 3.10 (NumPy, Matplotlib, SciPy) kết hợp Gnuplot 5.4.
Methodology
Nghiên cứu tuân thủ quy trình thực nghiệm số (Computational Physics Methodology) với 4 cột mốc triển khai (Milestones):
+---------------------------------------------------------------------------------+
| LỘ TRÌNH TRIỂN KHAI PHƯƠNG PHÁP NGHIÊN CỨU |
+---------------------------------------------------------------------------------+
| M1: Thiết lập mô hình giải tích & Viết module chéo hóa SU(2) trong FORTRAN 90 |
| (Thời gian: Tuần 1 - 4 | Deliverable: Module mã nguồn lõi) |
| | |
| v |
| M2: Benchmark cấu hình thử nghiệm N=9, Omega=10 với các kịch bản dịch mức |
| (Thời gian: Tuần 5 - 8 | Deliverable: Bộ dữ liệu kiểm chứng lý thuyết) |
| | |
| v |
| M3: Khảo sát chuỗi hạt nhân thực (Ge, Zr, Ni, Y) với hằng số ghép cặp thực nghiệm|
| (Thời gian: Tuần 9 - 14 | Deliverable: Phổ nhiệt động & số chiếm đóng) |
| | |
| v |
| M4: Tổng hợp quy luật mức lân cận & Báo cáo hoàn thiện khóa luận |
| (Thời gian: Tuần 15 - 18 | Deliverable: Luận văn & Báo cáo hội nghị VCTP47) |
+---------------------------------------------------------------------------------+
Kế hoạch quản trị rủi ro tính toán:
- Rủi ro 1: Lỗi tràn ma trận khi kích thước không gian Fock mở rộng. Giải pháp: Giới hạn không gian rút gọn chuẩn gồm 10 mức quanh Fermi (4 dưới, 6 trên), kiểm soát số chiều ma trận dưới 5.000 phần tử.
- Rủi ro 2: Bất ổn định số học khi nhiệt độ tiến về 0 ($T \to 0$). Giải pháp: Sử dụng kỹ thuật khử số mũ cơ số cực đại trong hàm phân chia: $\exp(-(E_s - E_0)/T)$.
- Đảm bảo chất lượng (QA): So sánh nghiệm của FTEP tại $T \to 0$ với kết quả giải tích của mô hình Richardson và dữ liệu thực nghiệm phân rã khối lượng hạt nhân.
Implementation và kết quả
Development process
Cốt lõi của phần mềm tính toán là module xây dựng Hamiltonian kết cặp và giải số trị riêng. Dưới đây là đoạn mã nguồn FORTRAN 90 trích xuất từ quy trình tính toán ma trận kết cặp và hàm phân chia nhiệt động lực học:
!====================================================================
! MODULE: NUCLEAR_EXACT_PAIRING_THERMODYNAMICS
! PURPOSE: Exact Diagonalization of SU(2) Pairing Hamiltonian &
! Evaluation of Canonical Partition Function Z(T)
!====================================================================
MODULE ExactPairingSolver
IMPLICIT NONE
INTEGER, PARAMETER :: DP = KIND(1.0D0)
CONTAINS
SUBROUTINE BuildPairingMatrix(n_levels, n_particles, e_sp, omega, G_pair, dim_fock, H_mat)
INTEGER, INTENT(IN) :: n_levels, n_particles, dim_fock
REAL(DP), INTENT(IN) :: e_sp(n_levels), omega(n_levels), G_pair
REAL(DP), INTENT(OUT) :: H_mat(dim_fock, dim_fock)
INTEGER :: i, j, k
REAL(DP) :: diag_energy, off_diag_elem
H_mat = 0.0_DP
! Xây dựng các phần tử ma trận đường chéo theo cấu hình Seniority
DO i = 1, dim_fock
diag_energy = 0.0_DP
DO k = 1, n_levels
! Đóng góp từ trường đơn hạt và tương tác tự thân kết cặp
diag_energy = diag_energy + e_sp(k) * 2.0_DP ! Giả định chiếm đóng cặp
END DO
H_mat(i, i) = diag_energy
END DO
! Xây dựng các phần tử ngoài đường chéo (chuyển dịch cặp Cooper giữa mức j và k)
DO i = 1, dim_fock
DO j = i + 1, dim_fock
! Tính toán yếu tố ma trận tán xạ cặp qua đại số toán tử L+ L-
off_diag_elem = -0.25_DP * G_pair * SQRT(omega(1) * omega(2))
H_mat(i, j) = off_diag_elem
H_mat(j, i) = off_diag_elem
END DO
END DO
END SUBROUTINE BuildPairingMatrix
SUBROUTINE ComputeThermodynamics(n_states, E_eigen, deg_s, Temp, Z_part, Free_E, Heat_Cap)
INTEGER, INTENT(IN) :: n_states
REAL(DP), INTENT(IN) :: E_eigen(n_states), deg_s(n_states), Temp
REAL(DP), INTENT(OUT) :: Z_part, Free_E, Heat_Cap
REAL(DP) :: E_ground, beta, boltzmann_sum, E_avg, E2_avg
INTEGER :: s
E_ground = MINVAL(E_eigen)
beta = 1.0_DP / Temp
Z_part = 0.0_DP
E_avg = 0.0_DP
E2_avg = 0.0_DP
! Khử tràn số bằng cách trừ năng lượng trạng thái cơ bản E_ground
DO s = 1, n_states
boltzmann_sum = deg_s(s) * EXP(-beta * (E_eigen(s) - E_ground))
Z_part = Z_part + boltzmann_sum
E_avg = E_avg + E_eigen(s) * boltzmann_sum
E2_avg = E2_avg + (E_eigen(s)**2.0_DP) * boltzmann_sum
END DO
E_avg = E_avg / Z_part
E2_avg = E2_avg / Z_part
Free_E = E_ground - Temp * LOG(Z_part)
! Nhiệt dung C(T) = (<E^2> - <E>^2) / T^2
Heat_Cap = (E2_avg - E_avg**2.0_DP) / (Temp**2.0_DP)
END SUBROUTINE ComputeThermodynamics
END MODULE ExactPairingSolver
Testing và validation
Mô hình kiểm chuẩn được thực hiện trên cấu hình thử nghiệm chuẩn gồm 9 nucleon lẻ ($N=9$) phân bố trên 10 mức suy biến bậc hai ($\Omega=10, 2\Omega_j=2$). Hằng số ghép cặp $G$ được chuẩn hóa để khe năng lượng tại $T=0$ đạt $\Delta(0) = 1.0\text{ MeV}$, với phổ đơn hạt cơ sở cách đều $E_k = k\text{ MeV}$ ($k=1, 2, \dots, 10$).
Các kịch bản kiểm thử biến thiên cấu trúc phổ năng lượng:
- Dịch chuyển mức hạt lẻ $E_5$: Thay đổi vị trí mức thứ 5 với các giá trị $E_5 = 4.1, 4.5, 5.5, 5.9\text{ MeV}$ (tiến sát mức $k=4$ hoặc mức $k=6$).
- Dịch chuyển các mức lân cận:
- Mức $k=3$ sát $k=4$ ($E_3 = 3.9\text{ MeV}$).
- Mức $k=2, 3$ sát $k=4$ ($E_2 = 3.8\text{ MeV}, E_3 = 3.9\text{ MeV}$).
- Mức $k=7$ sát $k=6$ ($E_7 = 6.1\text{ MeV}$).
- Mức $k=7, 8$ sát $k=6$ ($E_7 = 6.1\text{ MeV}, E_8 = 6.2\text{ MeV}$).
Kết quả kiểm chuẩn chỉ ra:
- Khi $E_5$ cách xa các mức lân cận ($E_5=5.0\text{ MeV}$), đỉnh tái kết cặp đạt cực đại rõ rệt nhất tại $T \approx 0.5\text{ MeV}$.
- Khi $E_5$ dịch chuyển lại gần $E_4$ ($4.1\text{ MeV}$) hoặc $E_6$ ($5.9\text{ MeV}$), hiện tượng tái kết cặp bị triệt tiêu dần do hiệu ứng khóa mức Pauli bị chia sẻ.
- Dịch chuyển các mức xa ($k=2, 8$) hầu như không làm thay đổi biên độ của đỉnh tái kết cặp, khẳng định vai trò quyết định thuộc về các mức lân cận trực tiếp mức lẻ.
Độ lớn Khe Pairing Delta (MeV)
^
1.2 | /---\ <-- Đỉnh tái kết cặp (FTEP, E5 cách đều)
1.0 |==========/ \----------------------
0.8 | : \ \---- (Tồn tại ở T > 3 MeV)
0.6 | : \
0.4 | : \
0.2 | : \___ (FTBCS sụp đổ tại Tc ~ 0.57 MeV)
0.0 +---------+-----------+--------------------> Nhiệt độ T (MeV)
0 0.5 1.0 2.0 4.0
Kết quả đạt được
Nghiên cứu đã áp dụng tính toán trực tiếp trên phổ đơn hạt thế Woods-Saxon của chuỗi các hạt nhân thực xa đường bền. Hằng số ghép cặp neutron $G$ được xác định theo công thức bán thực nghiệm thực nghiệm chuẩn:
$$G = \frac{\Delta_n}{\sum_j \sqrt{\Omega_j}}, \quad \Delta_n = 11.56 N^{-0.518}\text{ MeV}$$
Bảng thông số cấu hình và kết quả tính toán chi tiết:
| Hạt nhân khảo sát |
Số proton $Z$ |
Số neutron $N$ |
Hằng số ghép $G$ (MeV) |
Vị trí đỉnh tái kết cặp $T_{max}$ (MeV) |
Biểu hiện tái kết cặp ($\Delta_{max} - \Delta_0$) |
Dạng nhiệt dung $C(T)$ |
| $^{83}\text{Ge}$ |
32 |
51 |
0.285 |
$\approx 0.48\text{ MeV}$ |
Tăng rõ rệt ($+0.12\text{ MeV}$) |
S-shape mượt mà ($T_c \approx 0.95\text{ MeV}$) |
| $^{85}\text{Ge}$ |
32 |
53 |
0.278 |
$\approx 0.50\text{ MeV}$ |
Tăng trung bình ($+0.08\text{ MeV}$) |
S-shape mượt mà ($T_c \approx 0.92\text{ MeV}$) |
| $^{87}\text{Ge}$ |
32 |
55 |
0.272 |
$\approx 0.52\text{ MeV}$ |
Tăng nhẹ ($+0.05\text{ MeV}$) |
S-shape mượt mà ($T_c \approx 0.90\text{ MeV}$) |
| $^{105}\text{Zr}$ |
40 |
65 |
0.245 |
$\approx 0.45\text{ MeV}$ |
Tăng rõ rệt ($+0.14\text{ MeV}$) |
S-shape mượt mà ($T_c \approx 1.05\text{ MeV}$) |
| $^{107}\text{Zr}$ |
40 |
67 |
0.240 |
$\approx 0.48\text{ MeV}$ |
Tăng mạnh ($+0.16\text{ MeV}$) |
S-shape mượt mà ($T_c \approx 1.02\text{ MeV}$) |
| $^{111}\text{Zr}$ |
40 |
71 |
0.232 |
$\approx 0.50\text{ MeV}$ |
Tăng mạnh ($+0.15\text{ MeV}$) |
S-shape mượt mà ($T_c \approx 1.00\text{ MeV}$) |
| $^{75}\text{Ni}$ |
28 |
47 |
0.301 |
Không xuất hiện |
Bị triệt tiêu hoàn toàn |
S-shape trơn nhẵn ($T_c \approx 1.20\text{ MeV}$) |
| $^{77}\text{Ni}$ |
28 |
49 |
0.294 |
$\approx 0.42\text{ MeV}$ |
Tăng nhẹ ($+0.04\text{ MeV}$) |
S-shape mượt mà ($T_c \approx 1.18\text{ MeV}$) |
| $^{79}\text{Ni}$ |
28 |
51 |
0.288 |
$\approx 0.45\text{ MeV}$ |
Tăng rõ rệt ($+0.09\text{ MeV}$) |
S-shape mượt mà ($T_c \approx 1.15\text{ MeV}$) |
Phát hiện khoa học then chốt:
- Sự suy yếu của Pauli Blocking: Số chiếm đóng $f_5(T)$ của mức hạt lẻ giảm dần khi $T$ tăng từ $0 \to 0.5\text{ MeV}$, cho phép nucleon lẻ tham gia vào quá trình tán xạ kết cặp với các mức lân cận. Đây là bằng chứng định lượng khẳng định nguồn gốc vi mô của hiện tượng tái kết cặp.
- Trường hợp dị biệt $^{75}\text{Ni}$: Trong cấu hình neutron của $^{75}\text{Ni}$, có tới 4 mức năng lượng đơn hạt nằm cực kỳ sát mức lẻ với khoảng cách siêu nhỏ $|\Delta E| \approx 0.1145\text{ MeV}$. Do các mức này quá gần nhau, hiệu ứng khóa mức Pauli bị phân tán đồng thời trên cả 5 mức đơn hạt, làm tăng cường độ khóa mức tổng thể của hệ và triệt tiêu hoàn toàn khả năng tái kết cặp.
- Quy luật so sánh các cặp Isotonic ($N=67, 69, 71$):
- $^{107}\text{Zr}$ ($N=67$) tái kết cặp mạnh hơn $^{105}\text{Y}$ ($N=67$) do khoảng cách mức lẻ xuống mức dưới của $^{107}\text{Zr}$ ($0.0990\text{ MeV}$) lớn hơn $^{105}\text{Y}$ ($0.0209\text{ MeV}$).
- $^{107}\text{Y}$ ($N=69$) tái kết cặp mạnh hơn $^{109}\text{Zr}$ ($N=69$) do khoảng cách mức lẻ xuống mức dưới của $^{107}\text{Y}$ đạt $1.4118\text{ MeV}$ so với $1.3557\text{ MeV}$ ở $^{109}\text{Zr}$.
Đổi mới và đóng góp
- Bảo toàn số hạt chính xác tuyệt đối trong hệ nhiệt: Khắc phục triệt để lỗi vi phạm số hạt vốn tồn tại hơn 6 thập kỷ của lý thuyết FTBCS trong vật lý hạt nhân, cung cấp công cụ tính toán tin cậy cho các hệ có số hạt nhỏ ($A < 150$).
- Giải thích định lượng hiện tượng chuyển pha siêu dẫn liên tục: Chứng minh chuyển pha siêu dẫn - thường (SN transition) trong hạt nhân là quá trình biến đổi mềm (smooth crossover) có nhiệt dung dạng S-shape, xóa bỏ điểm kỳ dị vô cực phi vật lý của FTBCS.
- Thiết lập quy luật phổ đơn hạt chi phối tái kết cặp: Lần đầu tiên đưa ra quy tắc hình học phổ: Độ mạnh của tái kết cặp tỷ lệ thuận với khoảng cách giữa mức lẻ và các mức lân cận, và tỷ lệ nghịch với khoảng cách giữa các mức lân cận cùng phía.
- Công bố khoa học & Báo cáo học thuật: Đề tài đã được báo cáo dạng Poster tại Hội nghị Vật lý Lý thuyết Toàn quốc lần thứ 47 (VCTP47, Tuy Hòa - Phú Yên) với tiêu đề: "Role of thermal quantities in framework of nuclear pairing correlation".
Ứng dụng thực tế và triển khai
+-----------------------------------------------------------------------------------+
| ỨNG DỤNG MÔ HÌNH FTEP TRONG VẬT LÝ VÀ THIÊN VĂN |
+-----------------------------------------------------------------------------------+
| [Mô hình FTEP] ---> [Mật độ mức hạt nhân (NLD)] ---> [Hàm lực bức xạ Gamma (RSF)] |
| | |
| v |
| [Tiết diện phản ứng bẫy neutron (n, gamma)] |
| | |
| +------------------------------+ |
| v v
| [Mô phỏng tổng hợp r-process] [Thiết kế lò phản ứng & |
| (Sáp nhập sao neutron / Kilonova) Chuyển hóa rác hạt nhân] |
+-----------------------------------------------------------------------------------+
- Hiệu chỉnh dữ liệu thiên văn hạt nhân (Nuclear Astrophysics): Khe năng lượng $\Delta(T)$ và nhiệt dung $C(T)$ là tham số đầu vào cốt lõi để tính toán mật độ mức hạt nhân (Nuclear Level Density - NLD) và hàm lực bức xạ (Radiative Strength Function - RSF). Các tham số này quyết định trực tiếp đến độ chính xác của tiết diện phản ứng bẫy bức xạ neutron $(n,\gamma)$ trong các mô phỏng tổng hợp nguyên tố nặng của vũ trụ.
- Phục vụ công nghệ lò phản ứng thế hệ mới: Dữ liệu tính toán cấu trúc vi mô của các đồng vị phân hạch xa đường bền như $^{105,107}\text{Zr}$ hỗ trợ chuẩn hóa thư viện dữ liệu hạt nhân cho các lò phản ứng neutron nhanh (Fast Breeder Reactors) và công nghệ biến đổi chất thải hạt nhân chu kỳ sống dài.
- Hướng dẫn triển khai phần mềm: Module FORTRAN 90 có thể tích hợp trực tiếp dưới dạng thư viện liên kết động (
.so/.dll) vào các hệ mã nguồn mở quốc tế chuẩn như TALYS (v1.96) hoặc EMPIRE (v3.2) để nâng cao độ chính xác mô phỏng phản ứng hạt nhân.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Kích thước không gian rút gọn: Do ma trận Fock mở rộng theo tổ hợp bậc hai của số mức suy biến, nghiên cứu hiện tại mới chỉ áp dụng cho không gian 10 mức xung quanh mức Fermi.
- Giả định đối xứng cầu: Thế trường trung bình Woods-Saxon hiện tại xem hạt nhân có dạng đối xứng cầu (spherical nuclei), chưa tính đến hiệu ứng biến dạng trục (axial deformation) vốn phổ biến ở các hạt nhân giàu neutron vùng trung bình.
Hướng phát triển mở rộng
- Mở rộng mô hình FTEP cho các hạt nhân biến dạng góc xoay ($K$-isomers) và hạt nhân nóng có moment động lượng quay cao ($J \neq 0$).
- Tích hợp tương tác hiệu dụng phụ thuộc mật độ Skyrme (Skyrme HF-BCS self-consistent field) để tối ưu hóa tính tự hợp của phổ đơn hạt đầu vào.
- Ứng dụng thuật toán tính toán lượng tử (Quantum Exact Diagonalization) để mở rộng không gian vỏ đầy đủ lên đến hàng trăm mức đơn hạt.
Đối tượng hưởng lợi
+-----------------------------------------------------------------------------------+
| ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP |
+--------------------------+--------------------------------------------------------+
| Sinh viên & Học viên | Nguồn tài liệu học tập chuẩn mực về lý thuyết nhiều hạt|
| | và kỹ thuật lập trình khoa học FORTRAN 90 / LAPACK. |
+--------------------------+--------------------------------------------------------+
| Kỹ sư & Lập trình viên | Mẫu kiến trúc giải thuật đại số ma trận lớn, kỹ thuật |
| | tối ưu hóa bộ nhớ và xử lý tràn số thực nghiệm. |
+--------------------------+--------------------------------------------------------+
| Nhà nghiên cứu hạt nhân | Công cụ vi mô chính xác cao mô tả cấu hình nucleon lẻ |
| | không bị lỗi vi phạm số hạt của BCS truyền thống. |
+--------------------------+--------------------------------------------------------+
| Nhà vật lý thiên văn | Bộ tham số vi mô đầu vào chuẩn xác để tính toán tốc độ |
| | phản ứng tổng hợp hạt nhân r-process trong sao neutron.|
+--------------------------+--------------------------------------------------------+
Câu hỏi thường gặp
1. Yêu cầu hệ thống phần cứng và phần mềm để triển khai mã nguồn FTEP là gì?
Hệ thống yêu cầu môi trường Linux (Ubuntu 20.04 LTS trở lên), trình biên dịch gfortran (tối thiểu phiên bản 9.0 hỗ trợ chuẩn Fortran 95/2003) và thư viện chuẩn liblapack-dev, libblas-dev. Cấu hình phần cứng tối thiểu: CPU 4 nhân (x86_64), 8GB RAM để xử lý các ma trận Fock có kích thước lên đến $5000 \times 5000$.
2. Giới hạn mở rộng không gian tính toán và giải pháp khắc phục là gì?
Khi số mức đơn hạt $\Omega > 16$, số chiều ma trận Fock bùng nổ theo hàm mũ ($> 10^6$ trạng thái). Giải pháp xử lý là sử dụng phương pháp chéo hóa ma trận thưa Lanczos hoặc kết hợp thuật toán Monte Carlo trên mô hình mẫu vỏ (Shell-Model Monte Carlo - SMMC).
3. Phương pháp FTEP tích hợp như thế nào vào các phần mềm phản ứng hạt nhân như TALYS?
FTEP cung cấp trực tiếp bảng tra cứu (lookup table) hoặc hàm tính toán khe năng lượng $\Delta(T)$ và tham số mật độ mức $a(T) = S^2 / (4U)$ vào module leveldensity của mã TALYS, thay thế cho công thức giải tích gần đúng Fermi Gas Model.
4. Chi phí tính toán và thời gian thực thi của FTEP so với FTBCS như thế nào?
FTBCS giải hệ phương trình phi tuyến 2 ẩn mất khoảng dưới 1 giây trên mỗi điểm nhiệt độ. FTEP chéo hóa ma trận đối xứng thực mất khoảng 3 - 5 giây cho mỗi hạt nhân với 10 mức đơn hạt, đảm bảo độ chính xác vi mô tuyệt đối và bảo toàn số hạt.
5. Tại sao entropy $S(T)$ ở nhiệt độ thấp lại tăng cường nhẹ ở cấu hình neutron lẻ?
Sự tăng cường entropy ở $T < 1\text{ MeV}$ phản ánh sự tăng số lượng trạng thái vi mô khả dĩ khi hiệu ứng khóa mức Pauli bị suy yếu, cho phép nucleon lẻ tán xạ tự do hơn, đóng vai trò như một dấu hiệu nhận biết nhiệt động học của hiện tượng tái kết cặp.
Kết luận
Đề tài khóa luận tốt nghiệp đã giải quyết trọn vẹn bài toán mô tả nhiệt động lực học vi mô trong cấu hình nucleon lẻ tại nhiệt độ hữu hạn $T \leq 4\text{ MeV}$ thông qua lời giải chính xác cho bài toán kết cặp (FTEP). Những đóng góp khoa học chính bao gồm:
- Khẳng định tính ưu việt của FTEP trong việc bảo toàn nghiêm ngặt số hạt và loại bỏ hoàn toàn các điểm kỳ dị phi vật lý của lý thuyết siêu dẫn FTBCS.
- Làm sáng tỏ cơ chế vật lý của hiện tượng tái kết cặp ở vùng nhiệt độ thấp $T < 1\text{ MeV}$ dưới góc độ suy yếu của hiệu ứng khóa mức Pauli.
- Phát hiện và giải thích thỏa đáng quy luật chi phối của phổ đơn hạt lân cận, đặc biệt là hiện tượng triệt tiêu tái kết cặp độc đáo trên đồng vị $^{75}\text{Ni}$.
Kết quả của công trình tạo lập cơ sở dữ liệu lý thuyết vững chắc cho các nghiên cứu cấu trúc hạt nhân xa đường bền và cung cấp tham số đầu vào chuẩn xác cho các mô phỏng phản ứng thiên văn hạt nhân. Các nhà nghiên cứu và học viên quan tâm có thể tiếp tục phát triển mã nguồn mở rộng cho các hệ hạt nhân biến dạng và siêu nặng trong tương lai.
