Tổng quan nghiên cứu

Trong vật lý thiên văn hiện đại, quá trình tổng hợp hạt nhân trong các vụ nổ sao như siêu tân tinh loại II và bùng nổ tia X diễn ra dưới điều kiện nhiệt độ cực hạn vượt trên 1,0 tỷ Kelvin (tương ứng vùng nhiệt độ sao $T_9 > 1,0\text{ GK}$). Tại môi trường khắc nghiệt này, chu trình phản ứng hydro bùng nổ và quá trình bắt proton nhanh (rp-process) đóng vai trò quyết định đến sự hình thành các nguyên tố hóa học trong vũ trụ. Hạt nhân $^{22}\text{Mg}$ được xem là một điểm nghẽn tiềm năng trong chuỗi phản ứng tổng hợp từ $^{14}\text{O}$ đến $^{26}\text{Si}$. Phản ứng nhiệt hạch $^{22}\text{Mg}(\alpha,p)^{25}\text{Al}$ chính là mắt xích quan trọng giúp dòng phản ứng phá vỡ điểm nghẽn này để tiếp tục tổng hợp các nguyên tố nặng hơn, đồng thời chi phối trực tiếp đến hàm lượng đồng vị $^{22}\text{Na}$ và bức xạ tia gamma đặc trưng 1,275 MeV quan sát được từ không gian.

Tuy nhiên, các nghiên cứu trước đây về tốc độ phản ứng $^{22}\text{Mg}(\alpha,p)^{25}\text{Al}$ chủ yếu dựa trên các phép đo gián tiếp hoặc ngoại suy từ hạt nhân đối xứng gương $^{26}\text{Mg}$, để lại độ bất định rất lớn vượt quá 50%. Nguyên nhân là do thiếu hụt dữ liệu thực nghiệm trực tiếp về cấu trúc mức kích thích của hạt nhân $^{26}\text{Si}$ nằm trên ngưỡng tách hạt alpha tại mức năng lượng 9,164 MeV. Luận văn này tập trung giải quyết bài toán then chốt: thực hiện phép đo thực nghiệm trực tiếp phản ứng tán xạ $^{22}\text{Mg}(\alpha,\alpha)^{22}\text{Mg}$ bằng chùm ion phóng xạ tại hệ phổ kế CRIB thuộc Viện Nghiên cứu Vật lý và Hóa học RIKEN (Nhật Bản). Mục tiêu cụ thể là xác định các trạng thái cộng hưởng mới của $^{26}\text{Si}$ trong dải năng lượng khối tâm từ 0,5 MeV đến 3,0 MeV, từ đó tính toán chính xác tốc độ phản ứng $^{22}\text{Mg}(\alpha,p)^{25}\text{Al}$ ở nhiệt độ sao $T_9 = 0,5 - 3,0\text{ GK}$. Kết quả nghiên cứu giúp giảm độ bất định của tốc độ phản ứng xuống dưới 20%, cung cấp cơ sở dữ liệu thực nghiệm chuẩn xác để giải thích hiện tượng dị thường tỷ lệ đồng vị $^{22}\text{Ne}/^{20}\text{Ne}$ trong các mẫu thiên thạch cổ.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên nền tảng lý thuyết tổng hợp hạt nhân sao và cơ chế phản ứng hạt nhân qua trạng thái hợp phần. Trong môi trường plasma sao, các hạt mang điện chuyển động theo phân bố vận tốc Maxwell-Boltzmann và tương tác với nhau qua hiệu ứng xuyên ngầm hàng rào Coulomb lượng tử. Vùng năng lượng có xác suất phản ứng cao nhất được xác định bởi cửa sổ Gamow, được đặc trưng bởi hệ số thiên văn $S(E)$ và thông số Sommerfeld. Đối với các phản ứng hạt nhân ở vùng nhiệt độ cao, cơ chế cộng hưởng qua hạt nhân hợp phần đóng vai trò áp đảo, được mô tả thông qua công thức cộng hưởng Breit-Wigner với các tham số năng lượng cộng hưởng $E_R$, bề rộng phân rã hạt $\Gamma_\alpha$, $\Gamma_p$ và spin-tính chẵn lẻ $J^\pi$.

Để trích xuất các thông số cấu trúc vi mô từ phổ tán xạ thực nghiệm, khung lý thuyết ma trận R do Wigner và Eisenbud phát triển được áp dụng toàn diện. Lý thuyết ma trận R chia không gian cấu hình hạt nhân thành hai miền: miền bên trong (bán kính hạt nhân khoảng 5,5 fm) phản ánh cấu trúc mức năng lượng nội tại của hạt nhân hợp phần $^{26}\text{Si}$, và miền bên ngoài mô tả chuyển động tương đối của các kênh phản ứng phân rã. Việc phân tích ma trận R cho phép xác định chính xác vị trí các mức kích thích, độ rộng phân rã riêng phần và cấu trúc cụm alpha ($\alpha$-cluster) của $^{26}\text{Si}$ trên ngưỡng tách hạt 9,164 MeV mà không phụ thuộc vào các giả định mô hình thống kê hạt nhân thuần túy.

Phương pháp nghiên cứu

Nghiên cứu sử dụng phương pháp bia khí dày động học nghịch đảo (Thick Target in Inverse Kinematics - TTIK), một kỹ thuật tiên tiến giải quyết triệt để rào cản không thể chế tạo bia cố định từ đồng vị phóng xạ không bền $^{22}\text{Mg}$ (thời gian bán rã chỉ khoảng 3,86 giây). Trong cơ chế này, chùm ion phóng xạ $^{22}\text{Mg}^{12+}$ năng lượng cao được gia tốc và bắn phá vào bia khí Heli tĩnh có áp suất kiểm soát. Bia khí dày vừa đóng vai trò môi trường phản ứng, vừa làm suy giảm năng lượng chùm tia liên tục dọc theo quỹ đạo bay, cho phép quét toàn bộ dải năng lượng quan tâm từ 0,5 MeV đến 3,0 MeV trong hệ khối tâm chỉ qua một lần chạy chùm.

Dữ liệu thực nghiệm được thu thập từ hơn 1.500.000 sự kiện tương tác hạt tại hệ phổ kế CRIB của Đại học Tokyo đặt tại RIKEN. Chùm sơ cấp $^{20}\text{Ne}$ từ nguồn ion HyperECR được gia tốc qua máy gia tốc AVF Cyclotron đạt thông số $K = 70$, cường độ dòng đạt từ 100 enA đến 2 e$\mu$A để tạo chùm thứ cấp $^{22}\text{Mg}$ qua phản ứng hạt nhân tại bia sản xuất $F_0$. Hệ thống phân định hạt (Particle Identification - PID) kết hợp máy đo thời gian bay (ToF), buồng ion hóa tỷ lệ tích cực GEM-MSTPC và hệ kính thiên văn bán dẫn silicon nhiều dải (SSD telescopes) với độ phân giải năng lượng đạt dưới 50 keV. Việc lựa chọn phương pháp phân tích ma trận R thông qua mã chuyên dụng AZURE và tính toán tổn hao năng lượng bằng SRIM-2010 bảo đảm xử lý chính xác sự mất mát năng lượng của chùm hạt và trích xuất phổ tán xạ vi phân với độ tin cậy thống kê cao. Toàn bộ lộ trình thực nghiệm và xử lý dữ liệu được tiến hành liên tục trong giai đoạn 2011-2013.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Phép đo thực nghiệm trực tiếp đã mang lại những kết quả có tính đột phá về cấu trúc hạt nhân và tốc độ phản ứng thiên văn:

  • Xác lập phổ hàm kích thích tán xạ đàn hồi: Nghiên cứu đã đo đạc và xây dựng thành công hàm kích thích tiết diện vi phân của phản ứng $^{22}\text{Mg}(\alpha,\alpha)^{22}\text{Mg}$ trong hệ tọa độ khối tâm tại các góc tán xạ phòng thí nghiệm $\theta_{lab} = 0^\circ - 5^\circ$ và $5^\circ - 10^\circ$, bao phủ toàn diện dải năng lượng từ $E_{cm} = 0,8\text{ MeV}$ đến $2,8\text{ MeV}$.
  • Phát hiện các trạng thái cộng hưởng mới: Đã xác định được nhiều mức năng lượng kích thích mới của hạt nhân $^{26}\text{Si}$ nằm trên ngưỡng tách hạt alpha 9,164 MeV. Phân tích ma trận R chỉ ra các mức cộng hưởng quan trọng với gán trị spin-tính chẵn lẻ rõ ràng, tiêu biểu là các trạng thái $J^\pi = 2^+$ và $0^+$ ở vùng năng lượng tương ứng với cửa sổ Gamow của các siêu tân tinh.
  • Chứng minh sự tồn tại của cấu trúc cụm alpha: Tỷ số bề rộng phân rã thực nghiệm so với giới hạn Wigner đạt giá trị cao ở một số mức kích thích, khẳng định sự hiện diện mạnh mẽ của cấu trúc cụm $\alpha$ bên trong hạt nhân $^{26}\text{Si}$. Cấu trúc này làm tăng cường tiết diện phản ứng cộng hưởng lên gấp nhiều lần so với dự đoán lý thuyết thông thường.
  • Định lượng chính xác tốc độ phản ứng thiên văn: Tốc độ phản ứng $^{22}\text{Mg}(\alpha,p)^{25}\text{Al}$ tính toán từ các tham số thực nghiệm tại vùng nhiệt độ $T_9 = 1,0 - 3,0\text{ GK}$ cao hơn từ 200% đến 500% (tức gấp 2 đến 5 lần) so với các giá trị tính từ mô hình thống kê Hauser-Feshbach trong cơ sở dữ liệu NON-SMOKER.

Thảo luận kết quả

Các kết quả thực nghiệm có thể được biểu diễn trực quan thông qua đồ thị hàm kích thích tán xạ (đường cong biểu diễn tiết diện vi phân $d\sigma/d\Omega$ theo năng lượng khối tâm $E_{cm}$) với các đỉnh cộng hưởng sắc nét được khớp tối ưu bằng đường cong ma trận R. Đồng thời, bảng số liệu so sánh năng lượng mức, độ rộng $\Gamma_\alpha$ và hệ số thiên văn $S(E)$ giữa $^{26}\text{Si}$ và hạt nhân gương $^{26}\text{Mg}$ thể hiện rõ sự tương thích về cấu trúc lượng tử nhưng có sự dịch chuyển năng lượng Coulomb đặc trưng.

Nguyên nhân dẫn đến sự chênh lệch lớn giữa tốc độ phản ứng thực nghiệm và mô hình thống kê bắt nguồn từ mật độ mức ở vùng năng lượng thấp trên ngưỡng alpha chưa đủ dày đặc để áp dụng mô hình thống kê liên tục. Sự hiện diện của các cộng hưởng đơn lẻ có cấu trúc cụm alpha chiếm ưu thế tuyệt đối trong việc đóng góp vào tiết diện phản ứng tổng. So với các công trình trước đây chỉ đo đạc phân rã beta của $^{26}\text{P}$ với dải năng lượng hẹp, dữ liệu thực nghiệm trực tiếp này đã loại bỏ hoàn toàn các giả định ngoại suy không chắc chắn. Về mặt vật lý thiên văn, tốc độ phản ứng $^{22}\text{Mg}(\alpha,p)^{25}\text{Al}$ tăng cao khẳng định rằng ở nhiệt độ $T_9 > 1,2\text{ GK}$, nhánh phản ứng bắt hạt alpha chiếm ưu thế hoàn toàn so với phân rã beta của $^{22}\text{Mg}$. Điều này đồng nghĩa với việc dòng tổng hợp hạt nhân đã vượt qua điểm nghẽn $^{22}\text{Mg}$ một cách nhanh chóng, làm suy giảm lượng $^{22}\text{Na}$ tạo thành và giải thích thỏa đáng vì sao các vệ tinh thiên văn học không quan sát thấy tín hiệu gamma 1,275 MeV mạnh như các mô hình cũ dự báo.

Đề xuất và khuyến nghị

Dựa trên những phát hiện mới về tốc độ phản ứng hạt nhân thiên văn, luận văn đưa ra bốn đề xuất hành động cụ thể cho cộng đồng nghiên cứu:

  1. Mở rộng đo đạc trực tiếp các kênh phản ứng lân cận: Các nhóm nghiên cứu thực nghiệm tại các trung tâm gia tốc quốc tế cần xúc tiến đo đạc trực tiếp kênh phát proton $^{22}\text{Mg}(\alpha,p)^{25}\text{Al}$ và phản ứng bắt proton $^{25}\text{Al}(p,\gamma)^{26}\text{Si}$. Mục tiêu là kiểm chứng chéo các thông số bề rộng phân rã proton $\Gamma_p$ và đưa độ bất định của tốc độ phản ứng xuống dưới mức 10% trong vòng 24 tháng tới.
  2. Cập nhật mạng lưới phản ứng trong mô hình thiên văn: Các nhà vật lý thiên văn lý thuyết cần tích hợp ngay bộ tham số tốc độ phản ứng mới vào các phần mềm mô phỏng thủy động lực học 3D (như mã MESA hoặc FLASH). Việc này nhằm nâng cao độ chính xác dự báo tỷ lệ quang thông và phân bố đồng vị trong tàn dư siêu tân tinh lên trên 95% trong giai đoạn 2024-2026.
  3. Phát triển công nghệ phát hiện hạt tích cực trong nước: Viện Năng lượng Nguyên tử Việt Nam và các trường đại học trọng điểm cần đầu tư làm chủ công nghệ buồng ion hóa tích cực GEM-MSTPC kết hợp hệ bia khí dày. Mục tiêu là chế tạo và thử nghiệm thành công hệ đo hạt nhân độ nhạy cao tại máy gia tốc Tandem trong nước trước năm 2027 nhằm phục vụ đào tạo và nghiên cứu chuyên sâu.
  4. Xây dựng cơ sở dữ liệu số về phản ứng hạt nhân thiên văn: Bộ Khoa học và Công nghệ cùng Viện Nghiên cứu Hạt nhân cần chủ trì thiết lập cổng dữ liệu chuyên ngành kết nối với mạng lưới quốc tế JINA-CEE. Mục tiêu là chuẩn hóa và lưu trữ hơn 50 bộ số liệu phản ứng hạt nhân then chốt của các chu trình CNO và NeNa-MgAl định kỳ hàng năm.

Đối tượng nên tham khảo luận văn

Công trình nghiên cứu mang giá trị học thuật và thực tiễn cao, đặc biệt hữu ích cho bốn nhóm đối tượng sau:

  • Nghiên cứu sinh và học viên cao học ngành Vật lý hạt nhân: Luận văn là tài liệu tham khảo mẫu mực về quy trình thiết kế thí nghiệm chùm ion phóng xạ, phương pháp bia dày động học nghịch đảo và kỹ thuật phân tích ma trận R trên mã AZURE, giúp tiết kiệm khoảng 40% thời gian tiếp cận kỹ thuật phân tích phổ.
  • Kỹ sư và chuyên gia vận hành hệ thống máy gia tốc: Tài liệu cung cấp đầy đủ thông số kỹ thuật chi tiết của hệ phổ kế CRIB, từ trường quét 1,28 Tm, cấu trúc bộ lọc Wien Filter và hệ thống chân không, phục vụ trực tiếp cho việc thiết kế, bảo dưỡng và tối ưu hóa hiệu suất truyền chùm ion phụ.
  • Các nhà mô hình hóa vật lý thiên văn và vũ trụ học: Cung cấp bộ dữ liệu thực nghiệm chuẩn xác về tốc độ phản ứng $^{22}\text{Mg}(\alpha,p)^{25}\text{Al}$ ở dải nhiệt độ $T_9 = 0,5 - 3,0\text{ GK}$, hỗ trợ hiệu chỉnh chính xác các mô hình tiến hóa sao và tính toán độ phổ biến của các nguyên tố nhẹ trong vũ trụ.
  • Giảng viên đại học chuyên ngành Vật lý nguyên tử và hạt nhân: Nguồn tư liệu chuyên sâu gắn liền với 11 công trình công bố quốc tế, là giáo trình tham khảo giá trị cho các học phần Vật lý hạt nhân thực nghiệm nâng cao và Vật lý thiên văn hạt nhân tại bậc sau đại học.

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

Tại sao phương pháp bia dày động học nghịch đảo (TTIK) lại bắt buộc trong nghiên cứu này?
Hạt nhân $^{22}\text{Mg}$ là đồng vị phóng xạ có chu kỳ bán rã rất ngắn (3,86 giây), không thể dùng làm bia cố định trong phòng thí nghiệm. Kỹ thuật TTIK sử dụng chùm ion $^{22}\text{Mg}$ bắn vào bia khí Heli tĩnh, cho phép quét liên tục năng lượng khối tâm từ 0,5 MeV đến 3,0 MeV chỉ trong một lượt đo, bảo đảm độ chính xác cao và tối ưu hóa thời gian chùm máy gia tốc.

Ý nghĩa của ngưỡng năng lượng tách hạt alpha 9,164 MeV của $^{26}\text{Si}$ là gì?
Mức 9,164 MeV là ngưỡng năng lượng tối thiểu để hạt nhân hợp phần $^{26}\text{Si}$ phân rã phát hạt alpha hoặc được tạo thành từ tương tác $^{22}\text{Mg} + \alpha$. Việc xác định chính xác các trạng thái kích thích nằm trên ngưỡng 9,164 MeV là chìa khóa then chốt để tính toán các tham số cộng hưởng trong cửa sổ Gamow của các vụ nổ sao nhiệt độ cao.

Phản ứng $^{22}\text{Mg}(\alpha,p)^{25}\text{Al}$ giải thích vấn đề dị thường đồng vị Ne-E như thế nào?
Phản ứng $(\alpha,p)$ cạnh tranh trực tiếp với kênh phân rã beta của $^{22}\text{Mg}$ thành $^{22}\text{Na}$ (nguồn sinh ra $^{22}\text{Ne}$). Khi tốc độ phản ứng $(\alpha,p)$ đo được cao gấp 2 đến 5 lần so với lý thuyết ở nhiệt độ trên 1,0 GK, lượng $^{22}\text{Na}$ bị tiêu hao nhanh chóng, giải thích nguyên nhân tỷ lệ $^{22}\text{Ne}/^{20}\text{Ne}$ bị biến đổi mạnh trong các hạt bụi thiên thạch.

Phương pháp ma trận R hỗ trợ xử lý dữ liệu tán xạ như thế nào?
Ma trận R phân chia không gian phản ứng thành miền cấu hình trong và ngoài dựa trên bán kính tương tác hạt nhân (khoảng 5,5 fm). Phương pháp này giúp khớp đường cong tiết diện vi phân thực nghiệm một cách chặt chẽ, trích xuất chính xác các thông số lượng tử cơ bản gồm vị trí mức năng lượng, spin-tính chẵn lẻ $J^\pi$ và bề rộng phân rã $\Gamma_\alpha, \Gamma_p$.

Kết quả thực nghiệm này khác biệt gì so với các tính toán theo mô hình thống kê NON-SMOKER?
Mô hình thống kê NON-SMOKER giả định mật độ mức biến thiên liên tục và trơn nhẵn, trong khi thực nghiệm trực tiếp phát hiện các mức cộng hưởng cụm alpha riêng lẻ sắc nét. Do đó, tốc độ phản ứng thực tế ở dải $T_9 = 1,5 - 2,5\text{ GK}$ cao hơn ước tính thống kê từ 200% đến 500%, làm thay đổi hoàn toàn bức tranh động lực học trong vụ nổ siêu tân tinh.

Kết luận

  • Thực hiện thành công phép đo thực nghiệm trực tiếp đầu tiên trên thế giới đối với tương tác hạt nhân $^{22}\text{Mg} + \alpha$ tại hệ phổ kế CRIB (RIKEN, Nhật Bản).
  • Phát hiện và xác lập hệ thống các mức kích thích cộng hưởng mới của hạt nhân $^{26}\text{Si}$ nằm trên ngưỡng tách hạt alpha 9,164 MeV.
  • Xác định chính xác tốc độ phản ứng nhiệt hạch $^{22}\text{Mg}(\alpha,p)^{25}\text{Al}$ trong dải nhiệt độ sao $T_9 = 0,5 - 3,0\text{ GK}$, chỉ ra giá trị thực tế cao hơn từ 2 đến 5 lần so với mô hình lý thuyết thống kê.
  • Làm sáng tỏ cơ chế phá vỡ điểm nghẽn $^{22}\text{Mg}$ trong chu trình NeNa-MgAl và cung cấp lời giải thực nghiệm cho hiện tượng suy giảm bức xạ gamma 1,275 MeV trong vũ trụ.
  • Khẳng định tính ưu việt của phương pháp bia khí dày động học nghịch đảo kết hợp buồng đo GEM-MSTPC trong nghiên cứu vật lý hạt nhân năng lượng thấp.

Đóng góp cốt lõi của công trình là thiết lập cầu nối thực nghiệm vững chắc giữa vật lý cấu trúc hạt nhân vi mô và các mô hình tiến hóa thiên văn học vĩ mô. Hướng nghiên cứu tiếp theo sẽ tập trung triển khai đo đạc trực tiếp kênh phản ứng bắt proton của các đồng vị lân cận trong giai đoạn 2024-2026. Quý độc giả, các nhà nghiên cứu và học viên hãy cùng khai thác dữ liệu và phương pháp luận của luận văn để mở rộng các đề tài nghiên cứu tiên phong trong lĩnh vực vật lý thiên văn hạt nhân!