Bối cảnh và vấn đề nghiên cứu

Khóa luận tốt nghiệp "Phương pháp R-matrix cho nghiên cứu hiệu ứng phi định xứ của thế quang học nucleon – hạt nhân" do sinh viên Võ Thị Hải Nhật thực hiện dưới sự hướng dẫn khoa học của CN. Nguyễn Lê Anh tại Khoa Vật lí, Trường Đại học Sư phạm Thành phố Hồ Chí Minh (tháng 4 năm 2018), thuộc chuyên ngành Sư phạm Vật lí. Đề tài tập trung vào lĩnh vực vật lí hạt nhân lý thuyết và tính toán mô phỏng phản ứng hạt nhân.

Trong lịch sử phát triển của vật lí hạt nhân kể từ thí nghiệm tán xạ Rutherford, tán xạ nucleon lên hạt nhân là phương pháp thực nghiệm then chốt nhằm thăm dò cấu trúc nội tại của hạt nhân. Để mô tả trạng thái của hệ nucleon – hạt nhân trong quá trình tán xạ, việc thiết lập và giải phương trình Schrödinger đối với thế tương tác phi định xứ (nonlocal potential) đóng vai trò trung tâm. Các mô hình truyền thống như trường tự hợp Hartree-Fock hay gần đúng khối lượng hiệu dụng thường đối mặt với sự phức tạp tính toán hoặc hạn chế về độ chính xác.

Năm 1962, Perey và Buck đề xuất mẫu thế quang học phi định xứ (NLOMP) với nhân phi định xứ dạng Gauss, giải thích thành công tán xạ đàn hồi nucleon trên bia $^{208}\text{Pb}$ ở các năng lượng thấp 7 MeV và 14 MeV. Tuy nhiên, khi đối chiếu với dữ liệu thực nghiệm mới hơn trên dải năng lượng rộng, bộ thông số nguyên bản của Perey–Buck bộc lộ nhiều sai lệch. Đồng thời, việc giải số phương trình Schrödinger vi tích phân bằng các phương pháp lặp thông thường gặp trở ngại lớn do tốc độ và khả năng hội tụ của hàm sóng phụ thuộc chặt chẽ vào từng dạng thế cụ thể. Nhằm khắc phục hạn chế này, các nghiên cứu kế thừa (như hướng tiếp cận của Udagawa và các hiệu chỉnh gần đây kết hợp thông số hình học từ mẫu định xứ KD03 của Koning – Delaroche) đã xây dựng bộ thông số NLOMP mới độc lập với năng lượng tới trong khoảng 10 – 30 MeV cho các hạt nhân từ $^{27}\text{Al}$ đến $^{208}\text{Pb}$.

Song song đó, lý thuyết R-matrix – khởi nguồn từ công trình của Wigner và Eisenbud (1947) – ban đầu dùng để mô tả hiện tượng cộng hưởng hạt nhân, đã được mở rộng thành công cụ giải phương trình vi phân và vi tích phân trong cơ học lượng tử. Lý thuyết này phát triển theo hai nhánh: R-matrix hiện tượng luận (phenomenological R-matrix) và R-matrix tính toán (calculable R-matrix). Nhược điểm cố hữu của R-matrix hiện tượng luận là kết quả phụ thuộc vào việc lựa chọn bán kính kênh, trong khi R-matrix tính toán (được đề xuất từ năm 1965 bởi Haglund và Robson, sau đó kết hợp với kỹ thuật lưới Lagrange) mang lại lời giải số trực tiếp, nhanh chóng và chính xác cho các thế phi định xứ mà không cần địa phương hóa thế năng.

Từ bối cảnh học thuật trên, tác giả xác định các nhiệm vụ nghiên cứu cụ thể được đánh số như sau:

  1. Trình bày tổng quan cơ sở lý thuyết tán xạ cổ điển và tán xạ lượng tử, thiết lập mối liên hệ giải tích giữa ma trận tán xạ cơ sở ($S$-matrix) và ma trận phản ứng ($R$-matrix).
  2. Xây dựng thế tương tác quang học phi định xứ (NLOMP) cho hệ nucleon – hạt nhân, phân tích dạng hàm giải tích của thế Perey–Buck gồm các thành phần thế khối (thực và ảo), thế ảo bề mặt, thế spin – quỹ đạo và thế Coulomb.
  3. Ứng dụng hình thức luận R-matrix tính toán kết hợp phương pháp lưới Lagrange (Lagrange-mesh method) để giải phương trình Schrödinger vi tích phân có chứa thế phi định xứ.
  4. Triển khai thuật toán tính toán số, xác định độ lệch pha và tiết diện tán xạ đàn hồi phân bố theo góc tán xạ vi phân và theo xung lượng chuyển đối với chùm hạt proton và neutron.
  5. So sánh, đánh giá kết quả tính toán số với cơ sở dữ liệu thực nghiệm EXFOR của Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA).

Đối tượng nghiên cứu của đề tài là quá trình tán xạ đàn hồi của nucleon (proton và neutron) lên các hạt nhân bia chẵn – chẵn: $^{40}\text{Ca}$, $^{56}\text{Fe}$, $^{120}\text{Sn}$, $^{208}\text{Pb}$. Phạm vi nghiên cứu giới hạn ở bài toán đơn kênh (single-channel), bỏ qua các hiệu ứng cộng hưởng, trạng thái kích thích do phản ứng trực tiếp và các phản ứng chuyển kênh hạt. Dải năng lượng khảo sát tập trung chủ yếu vào vùng năng lượng trung bình từ 10 MeV đến 40 MeV, đồng thời mở rộng kiểm chứng ở một số mức năng lượng thấp (9,943 MeV đối với xung lượng chuyển) và năng lượng cao (lên đến 65 MeV đối với proton và 96 MeV đối với neutron).

Cơ sở lý thuyết và phương pháp

Khung lý thuyết của khóa luận được xây dựng dựa trên sự kết hợp giữa cơ học lượng tử tán xạ, mẫu thế quang học phi định xứ và phương pháp R-matrix tính toán.

Trong lý thuyết tán xạ lượng tử, hàm sóng tiệm cận ở khoảng cách lớn $r \to \infty$ được phân tích thành sóng phẳng tới và sóng cầu tán xạ: $$\psi^{(+)}(\mathbf{r}) \xrightarrow{r \to \infty} e^{ikz} + f(\theta)\frac{e^{ikr}}{r}$$ với $k = \sqrt{2\mu E}/\hbar$ là số sóng, $f(\theta)$ là biên độ tán xạ. Bằng phép khai triển sóng từng phần theo đa thức Legendre $P_l(\cos\theta)$, bài toán quy về giải phương trình vi phân hướng tâm cho hàm sóng thành phần $u_l(r)$. Mối liên hệ giữa các phần tử của ma trận tán xạ $S_l$ và độ lệch pha $\delta_l$ được xác định qua $S_l = e^{2i\delta_l}$, từ đó tính được tiết diện tán xạ đàn hồi vi phân: $$\frac{d\sigma}{d\Omega} = |f(\theta)|^2 = \left| \frac{1}{2ik} \sum_{l=0}^\infty (2l+1) P_l(\cos\theta) (S_l - 1) \right|^2$$

Để mô tả tương tác giữa nucleon và hạt nhân bia, tác giả áp dụng mẫu quang học phi định xứ (NLOMP) Perey–Buck. Thế tương tác toàn phần có dạng: $$V_{\text{opt}} = V_C(r) + \int V(\mathbf{r}, \mathbf{r}') \psi(\mathbf{r}') d\mathbf{r}' + V_{SO}(r) (\mathbf{L} \cdot \mathbf{s})$$ Trong đó, thế Coulomb $V_C(r)$ áp dụng cho proton với bán kính điện tích $R_C = 1{,}2 A^{1/3}$ fm; thế spin – quỹ đạo $V_{SO}(r)$ có dạng đạo hàm của hàm Woods–Saxon; thế phi định xứ $V(\mathbf{r}, \mathbf{r}')$ được tách thành thừa số định xứ $U(p)$ (với $p = |\mathbf{r} + \mathbf{r}'|/2$) nhân với hàm phi định xứ chuẩn hóa dạng Gauss: $$H(|\mathbf{r} - \mathbf{r}'|) = \frac{1}{\pi^{3/2} \beta^3} \exp\left( - \frac{|\mathbf{r} - \mathbf{r}'|^2}{\beta^2} \right)$$ với $\beta$ là thông số đặc trưng cho tầm phi định xứ.

Phương pháp R-matrix tính toán chia không gian tương tác thành hai miền phân cách bởi bán kính kênh $a$:

  • Miền ngoài ($r > a$): Lực hạt nhân tầm ngắn triệt tiêu, hạt chỉ chịu tương tác Coulomb tầm xa. Hàm sóng ngoài $u_l^{\text{ext}}(r)$ là tổ hợp tuyến tính của hàm Coulomb tới $I_l(kr)$ và hàm Coulomb lui $O_l(kr)$ thông qua ma trận tán xạ $S_l$: $$u_l^{\text{ext}}(r) = I_l(kr) - S_l O_l(kr)$$
  • Miền trong ($r \le a$): Tồn tại đầy đủ các tương tác hạt nhân phi định xứ. Để khắc phục tính không Hermite của toán tử động năng do biên hữu hạn, toán tử Bloch $\mathcal{L}(B)$ được đưa vào: $$\mathcal{L}(B) = \frac{\hbar^2}{2\mu} \delta(r - a) \left( \frac{d}{dr} - \frac{B}{r} \right)$$ với $B$ là thông số biên. Khi đó toán tử $H_l + \mathcal{L}(B)$ trở thành toán tử Hermite trong miền trong. Hàm sóng miền trong được khai triển qua hệ cơ sở Lagrange gồm $N$ hàm $\varphi_j(r)$ xây dựng từ các đa thức Legendre trực giao trên các điểm cầu phương Gauss–Legendre $x_n \in [0, 1]$.

Nhờ tính chất trực giao của hàm cơ sở Lagrange, tích phân của thế phi định xứ được rút gọn thành phép nhân ma trận đại số mà không cần tính tích phân số phức tạp: $$\langle \varphi_n | V | \varphi_m \rangle \approx \lambda_n^{1/2} \lambda_m^{1/2} a V(a x_n, a x_m)$$ Từ ma trận $C_{ij}(E, B) = \langle \varphi_i | T_l + V(r) + \mathcal{L}(B) - E | \varphi_j \rangle$, nghịch đảo ma trận $C^{-1}$ cho phép tính trực tiếp phần tử ma trận $R_l(E, B)$ tại biên $r = a$: $$R_l(E, B) = \frac{\hbar^2}{2\mu a} \sum_{i,j=1}^N \varphi_i(a) C_{ij}^{-1}(E, B) \varphi_j(a)$$ Qua đó, ma trận tán xạ $S_l$ được xác định giải tích từ $R_l$ và các hàm sóng tiệm cận tại biên.

Phương pháp nghiên cứu thực tế là phương pháp giải số trên máy tính. Tác giả sử dụng chương trình nguồn viết bằng ngôn ngữ lập trình Fortran được phát triển bởi nhóm nghiên cứu tại Viện Khoa học và Kỹ thuật Hạt nhân (Hà Nội) kết hợp các tài liệu chuyên khảo. Số liệu đầu ra được xử lý, nội suy và biểu diễn đồ thị thông qua phần mềm Origin.

Nguồn dữ liệu thực nghiệm dùng để kiểm chứng mô hình được trích xuất từ hệ thống cơ sở dữ liệu phản ứng hạt nhân thực nghiệm EXFOR thuộc Trung tâm Dữ liệu Hạt nhân của Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA).

Nội dung chính theo từng chương

Chương 1: Lí thuyết tán xạ

Chương 1 thiết lập toàn bộ cơ sở lý thuyết vi mô, các công thức toán học và thuật toán số của đề tài qua 3 khối nội dung chính:

  1. Lí thuyết tán xạ cổ điển và lượng tử: Thiết lập định nghĩa tiết diện tán xạ vi phân $d\sigma/d\Omega$ và toàn phần $\sigma$. Trong cơ học lượng tử, phương trình Schrödinger ba chiều được tách biến trong hệ tọa độ cầu thành phần góc (biểu diễn qua hàm cầu $Y_l^m(\theta, \phi)$ và đa thức Legendre liên kết $P_l^m(\cos\theta)$) và phương trình vi phân hướng tâm cho $u_l(r)$. Nghiệm của phương trình tự do được biểu diễn qua hàm cầu Bessel $j_l(kr)$ và hàm cầu Neumann $n_l(kr)$. Chương này định nghĩa ma trận tán xạ $S_l$, góc lệch pha $\delta_l$ và thiết lập biểu thức toán học liên hệ giữa ma trận $S_l$ và ma trận $R_l$.
  2. Thế tương tác phi định xứ (NLOMP): Phân tích sự hình thành trường thế quang học từ tương tác đa hạt hiệu dụng. Trình bày dạng giải tích của hàm xuyên tâm Woods–Saxon $f(r, R, a) = [1 + \exp((r-R)/a)]^{-1}$. Trình bày chi tiết mẫu thế phi định xứ Perey–Buck với thành phần thực xuyên tâm $V_R$, thế ảo xuyên tâm $W_I$, thế ảo bề mặt $W_D$ (dạng đạo hàm Woods–Saxon), thế spin – quỹ đạo $V_{SO}$ và thế Coulomb $V_C(r)$.
  3. Hình thức luận R-matrix và phương pháp Lagrange-mesh: Trình bày nguyên lý chia không gian bởi bán kính kênh $a$; định nghĩa toán tử Bloch $\mathcal{L}(B)$; thiết lập hàm Green $G_l(r, r')$ và chứng minh giải tích biểu thức xác định ma trận $C(E, B)$. Trình bày kỹ thuật tính ma trận động năng và thế năng trên lưới tọa độ Gauss–Legendre. Mô tả quy trình thuật toán R-matrix, định dạng dữ liệu đầu vào (file input gồm năng lượng $E$, số khối $A$, điện tích $Z$, bán kính Coulomb $R_C$, các thông số lưới $N_R$, số miền $N_S$, bán kính kênh $a$, mô men quán tính cực đại $l_{\max}$ và 12 thông số thế NLOMP) cùng cấu trúc dữ liệu đầu ra (file output chứa độ lệch pha và tiết diện tán xạ vi phân).

Chương 2: Kết quả

Chương 2 trình bày toàn bộ số liệu tính toán số và kết quả so sánh với dữ liệu thực nghiệm IAEA qua 4 tiểu mục:

  1. Tán xạ proton lên hạt nhân: Khảo sát tiết diện tán xạ đàn hồi phân bố theo góc của proton trên 4 hạt nhân bia chẵn – chẵn:
    • $^{40}\text{Ca}$ ở các mức năng lượng trung bình.
    • $^{56}\text{Fe}$ ở 16 MeV và 65 MeV.
    • $^{208}\text{Pb}$ ở 30,3 MeV, 39,6 MeV và 45 MeV.
    • $^{120}\text{Sn}$ ở 30,6 MeV và 36,4 MeV. Kết quả cho thấy đường cong lý thuyết R-matrix tái hiện rất chính xác các đỉnh và đáy dao động nhiễu xạ của thực nghiệm trong dải năng lượng 10 – 40 MeV. Ở năng lượng cao trên 40 MeV (như 45 MeV trên $^{208}\text{Pb}$ hay 65 MeV trên $^{56}\text{Fe}$), xuất hiện sự lệch pha và chênh lệch biên độ ở các góc tán xạ lớn do hiệu ứng truyền qua rào Coulomb gia tăng và sự mở của các kênh phản ứng phi đàn hồi.
  2. Tán xạ neutron lên hạt nhân: Do neutron không mang điện tích ($Z=0$), thế tương tác không chứa thành phần Coulomb ($R_C = 0$), giúp giảm bớt số lượng thông số đầu vào. Khảo sát được thực hiện trên:
    • $^{120}\text{Sn}$ ở năng lượng thấp và trung bình.
    • $^{56}\text{Fe}$ ở 11 MeV, 20 MeV, 26 MeV và 96 MeV.
    • $^{208}\text{Pb}$ ở 30,3 MeV và 40 MeV.
    • $^{40}\text{Ca}$ ở 16,8 MeV. Độ khớp giữa tính toán R-matrix và thực nghiệm của neutron đạt mức độ chính xác cao hơn so với proton. Riêng đối với hạt nhân $^{40}\text{Ca}$, kết quả tính toán có sự sai khác đáng kể so với thực nghiệm do $^{40}\text{Ca}$ là hạt nhân tương đối nhẹ, cấu trúc "mềm", khả năng đâm xuyên của neutron lớn dễ kích hoạt các phản ứng liên kênh làm thay đổi trạng thái bia.
  3. Tiết diện tán xạ phân bố theo xung lượng chuyển: Tiến hành kiểm chứng mở rộng cho tán xạ đàn hồi neutron trên hạt nhân bia $^{120}\text{Sn}$ ở mức năng lượng $E_n = 9{,}943$ MeV theo biến số xung lượng chuyển $q = 2k \sin(\theta/2)$ trong dải $0{,}0 - 3{,}5\text{ fm}^{-1}$. Đường cong tính toán hoàn toàn trùng khớp với phân bố thực nghiệm EXFOR.
  4. Tổng kết: Khẳng định tính hiệu quả của phương pháp R-matrix tính toán kết hợp lưới Lagrange trong việc giải trực tiếp bài toán tán xạ đơn kênh với thế phi định xứ; nêu bật triển vọng mở rộng mô hình cho các bài toán liên kênh phức tạp.

Kết quả và đóng góp

Khóa luận đã chứng minh rằng phương pháp R-matrix tính toán giải quyết triệt để bài toán tán xạ đàn hồi nucleon trên hạt nhân có tính đến hiệu ứng phi định xứ của thế tương tác. Các tính toán số không gặp phải hiện tượng phân kỳ hay phụ thuộc vào dạng thế như các phương pháp lặp giải tích trước đây.

Bảng dưới đây tổng hợp các trường hợp hạt nhân bia, loại nucleon tới và mức năng lượng đã được tính toán và kiểm chứng trong khóa luận:

Hạt nhân bia Loại hạt tới Năng lượng hạt tới khảo sát ($E$) Mức độ phù hợp với dữ liệu thực nghiệm IAEA
$^{40}\text{Ca}$ Proton ($p$) Dải năng lượng trung bình Phù hợp tốt ở góc nhỏ và vừa
$^{40}\text{Ca}$ Neutron ($n$) 16,8 MeV Xuất hiện sai lệch do hiệu ứng cấu trúc mềm và liên kênh
$^{56}\text{Fe}$ Proton ($p$) 16 MeV; 65 MeV Khớp rất tốt ở 16 MeV; lệch ở góc lớn tại 65 MeV
$^{56}\text{Fe}$ Neutron ($n$) 11 MeV; 20 MeV; 26 MeV; 96 MeV Phù hợp tốt từ 11 đến 26 MeV; giảm chính xác ở 96 MeV
$^{120}\text{Sn}$ Proton ($p$) 30,6 MeV; 36,4 MeV Phù hợp chính xác cao toàn bộ dải góc đo
$^{120}\text{Sn}$ Neutron ($n$) 9,943 MeV (theo xung lượng chuyển) Trùng khớp hoàn hảo với các điểm thực nghiệm
$^{208}\text{Pb}$ Proton ($p$) 30,3 MeV; 39,6 MeV; 45 MeV Khớp tốt ở 30,3 và 39,6 MeV; sai lệch nhẹ ở 45 MeV
$^{208}\text{Pb}$ Neutron ($n$) 30,3 MeV; 40 MeV Mô tả chính xác các đỉnh nhiễu xạ thực nghiệm

Bảng thông số thế quang học phi định xứ (NLOMP) cho bài toán tán xạ neutron được tác giả sử dụng (trích xuất từ Phụ chú 4 của văn bản):

Loại thông số Độ sâu thế (MeV) Bán kính hình học (fm) Độ nhòe bề mặt (fm) Thông số phi định xứ $\beta$ (fm)
Thế thực xuyên tâm ($V_R$) 70,74 $r_R$ xác định từ KD03 $a_R$ xác định từ KD03 0,85
Thế ảo xuyên tâm ($W_I$) Xác định theo bảng $r_I$ xác định từ KD03 $a_I$ xác định từ KD03 0,85
Thế ảo bề mặt ($W_D$) Xác định theo bảng $r_D$ xác định từ KD03 $a_D$ xác định từ KD03 0,85
Thế spin – quỹ đạo ($V_{SO}$) Xác định theo bảng $r_{SO}$ xác định từ KD03 $a_{SO}$ xác định từ KD03 -

Đóng góp khoa học chính mà tác giả tự nêu so với các nghiên cứu trước bao gồm:

  • Khắc phục triệt để nhược điểm phụ thuộc vào bán kính kênh $a$ của lý thuyết R-matrix hiện tượng luận, bảo toàn đầy đủ ý nghĩa vật lý của các đại lượng tương tác.
  • Giải quyết bài toán hội tụ hàm sóng khi giải phương trình vi tích phân thế Perey–Buck bằng cách đưa vào toán tử Hermite Bloch $\mathcal{L}(B)$ và kỹ thuật tính toán lưới Lagrange, cho phép giải trực tiếp thế phi định xứ mà không cần dùng gần đúng địa phương hóa tương đương.
  • Khẳng định tính độc lập với năng lượng của bộ thông số NLOMP hiệu chỉnh trong vùng năng lượng 10 – 40 MeV đối với các hạt nhân có số khối $A \ge 27$.

Hạn chế và hướng nghiên cứu tiếp

Văn bản khóa luận đã thẳng thắn chỉ ra các giới hạn nghiên cứu cụ thể:

  • Phạm vi đối tượng bia: Đề tài chỉ áp dụng cho các hạt nhân bia chẵn – chẵn có spin trạng thái cơ bản $J^\pi = 0^+$ ($^{40}\text{Ca}$, $^{56}\text{Fe}$, $^{120}\text{Sn}$, $^{208}\text{Pb}$) nhằm triệt tiêu thành phần tương tác spin giữa nucleon tới và hạt nhân bia. Do đó, mô hình chưa thể khẳng định tính đúng đắn khi mở rộng cho các hạt nhân có spin khác không.
  • Giới hạn gần đúng đơn kênh: Khóa luận chỉ khảo sát bài toán đơn kênh thuần túy cho tán xạ đàn hồi. Khi năng lượng nucleon vượt quá 40 MeV hoặc khi tương tác với các hạt nhân có mật độ nucleon "mềm" như $^{40}\text{Ca}$, xác suất đâm xuyên qua rào Coulomb và xác suất xảy ra phản ứng không đàn hồi, phản ứng chuyển kênh tăng cao, khiến kết quả tính toán đơn kênh bị lệch so với thực nghiệm.
  • Dữ liệu góc tán xạ lớn: Ở các góc tán xạ lớn, một số bài toán tính toán bị lệch so với dữ liệu thực nghiệm, hoặc gặp tình trạng thiếu hụt số liệu đo đạc thực nghiệm từ EXFOR để đối sánh đầy đủ.

Các hướng phát triển tiếp theo được đề xuất trong khóa luận:

  • Tiếp tục hiệu chỉnh tinh vi bộ thông số NLOMP để nâng cao độ chính xác và mở rộng khả năng mô tả cho các nucleon ở vùng năng lượng cao hơn 40 MeV.
  • Thiết lập tương tác spin giữa nucleon và hạt nhân bia để mở rộng bài toán tán xạ đàn hồi trên các hạt nhân có số khối lẻ hoặc hạt nhân có spin trạng thái cơ bản khác 0.
  • Phát triển phương pháp R-matrix tính toán để giải hệ phương trình Schrödinger liên kênh (coupled-channel Schrödinger equations), mô tả đồng thời tán xạ đàn hồi và các kênh kích thích trạng thái liên tục.

Giá trị tham khảo

Bản khóa luận là tài liệu học thuật có giá trị tham khảo chuyên sâu đối với:

  • Sinh viên và học viên cao học ngành Sư phạm Vật lí, Vật lí hạt nhân và Vật lí tính toán: Cung cấp tài liệu tham khảo chi tiết về kỹ thuật giải số phương trình Schrödinger vi tích phân, phương pháp xây dựng ma trận R tính toán và phương pháp lưới Lagrange.
  • Người làm nghiên cứu về phản ứng hạt nhân: Các phần trình bày về toán tử Bloch, phép biến đổi ma trận $C_{ij}$, cách thiết lập các thông số thế Perey–Buck và bộ thông số quang học kết hợp KD03 là tài liệu hữu ích cho việc xây dựng mã nguồn mô phỏng phản ứng hạt nhân.
  • Phần đáng tham khảo nhất:
    • Mục 1.3 và 1.4 của Chương 1 cùng Phụ chú 5 trình bày tường minh các bước biến đổi toán học từ phương trình Schrödinger qua hàm Green đến công thức giải tích xác định ma trận $C(E, B)$ và phần tử $R_l(E, B)$.
    • Cấu trúc mẫu các tệp dữ liệu đầu vào (input) và đầu ra (output) của chương trình Fortran trong Chương 1 (Hình 3.2, 3.3, 3.4), cung cấp quy chuẩn thực hành tính toán số rõ ràng.

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

1. Vì sao khóa luận chỉ lựa chọn các hạt nhân bia chẵn – chẵn ($^{40}\text{Ca}$, $^{56}\text{Fe}$, $^{120}\text{Sn}$, $^{208}\text{Pb}$) để khảo sát?
Các hạt nhân chẵn – chẵn có spin và độ chẵn lẻ của trạng thái cơ bản bằng $0^+$. Việc chọn các hạt nhân này cho phép bài toán bỏ qua tương tác spin – spin giữa nucleon bắn tới và hạt nhân bia, giúp đơn giản hóa hệ phương trình vi phân hướng tâm và tập trung kiểm chứng hiệu ứng phi định xứ của thế xuyên tâm và thế spin – quỹ đạo.

2. Toán tử Bloch $\mathcal{L}(B)$ đóng vai trò gì trong phương pháp R-matrix tính toán?
Trong miền trong hữu hạn ($r \le a$), toán tử động năng thông thường $T_l$ không thỏa mãn tính chất Hermite do sự xuất hiện của các số hạng tích phân mặt tại biên $r = a$. Toán tử Bloch $\mathcal{L}(B) = \frac{\hbar^2}{2\mu} \delta(r - a) \left( \frac{d}{dr} - \frac{B}{r} \right)$ được cộng vào Hamiltonian để bù trừ chính xác các số hạng tích phân mặt không Hermite đó, đảm bảo toán tử toàn phần $H_l + \mathcal{L}(B)$ là Hermite, từ đó xác định nghiệm hàm sóng duy nhất thông qua hàm Green.

3. Phương pháp lưới Lagrange (Lagrange-mesh) mang lại ưu thế kỹ thuật gì khi xử lý thế phi định xứ?
Thế phi định xứ Perey–Buck có dạng tích phân nhân $\int V(\mathbf{r}, \mathbf{r}') \psi(\mathbf{r}') d\mathbf{r}'$, thông thường đòi hỏi giải thuật tích phân số hai lớp rất tốn kém tài nguyên tính toán. Khi khai triển trên hệ cơ sở Lagrange liên kết với các điểm cầu phương Gauss–Legendre, ma trận thế năng phi định xứ được tính trực tiếp từ giá trị thế tại các điểm nút lưới theo công thức đại số $V_{nm} \approx \lambda_n^{1/2} \lambda_m^{1/2} a V(a x_n, a x_m)$, giúp tăng tốc độ tính toán và đảm bảo độ hội tụ cao.

4. Vì sao kết quả tính toán tán xạ neutron phù hợp với thực nghiệm tốt hơn so với proton?
Hạt neutron không mang điện tích nên trong tương tác neutron – hạt nhân hoàn toàn không có thế Coulomb ($V_C(r) = 0$). Điều này giúp loại bỏ thông số bán kính Coulomb $R_C$ ra khỏi tệp thông số đầu vào. Việc giảm bớt số lượng thông số tự do giúp hệ phương trình ít bị nhiễu sai số và nâng cao độ chính xác của đường khớp lý thuyết so với thực nghiệm.

5. Yếu tố vật lý nào dẫn đến sự sai lệch giữa mô hình tính toán và thực nghiệm ở năng lượng trên 40 MeV và tại bia $^{40}\text{Ca}$?
Khi năng lượng nucleon vượt quá 40 MeV, xác suất truyền qua rào Coulomb tăng cao, làm kích hoạt các phản ứng phi đàn hồi và phản ứng chuyển hạt khiến bài toán chuyển thành bài toán liên kênh mà mô hình đơn kênh không mô tả được. Đối với $^{40}\text{Ca}$, đây là hạt nhân tương đối nhẹ có cấu trúc "mềm", khả năng đâm xuyên sâu của neutron dễ làm biến dạng cấu trúc hạt nhân bia, vi phạm giả thiết hạt nhân bia giữ nguyên trạng thái cơ bản của bài toán tán xạ đàn hồi thuần túy.

Kết luận

Khóa luận tốt nghiệp của tác giả Võ Thị Hải Nhật đã xây dựng thành công quy trình tính toán số giải phương trình Schrödinger vi tích phân có thế quang học phi định xứ Perey–Buck bằng phương pháp R-matrix tính toán kết hợp lưới Lagrange. Kết quả tính toán tiết diện tán xạ đàn hồi của proton và neutron trên các hạt nhân bia chẵn – chẵn ($^{40}\text{Ca}$, $^{56}\text{Fe}$, $^{120}\text{Sn}$, $^{208}\text{Pb}$) tái hiện chính xác các dữ liệu thực nghiệm từ cơ sở dữ liệu EXFOR của IAEA trong khoảng năng lượng từ 10 đến 40 MeV. Công trình khẳng định tính ưu việt của phương pháp R-matrix tính toán trong việc xử lý hiệu ứng phi định xứ mà không cần thông qua gần đúng địa phương hóa, đồng thời tạo tiền đề phương pháp luận vững chắc cho việc phát triển các mô hình phản ứng hạt nhân liên kênh phức tạp hơn.