Tổng quan nghiên cứu

Nghiên cứu cấu trúc và phản ứng hạt nhân ở dải năng lượng thấp dưới 50 MeV giữ vai trò then chốt trong vật lý hạt nhân hiện đại, đặc biệt khi các phòng thí nghiệm gia tốc thế hệ mới như FRIB tại Đại học Bang Michigan hay Viện RIKEN mở rộng khám phá các hạt nhân phóng xạ nằm xa đường bền đồng vị. Trong nhiều thập kỷ qua, các tính toán tiết diện phản ứng phụ thuộc chủ yếu vào thế quang học hiện tượng từ cơ sở dữ liệu quốc tế RIPL-3. Tuy nhiên, các mô hình hiện tượng dựa trên việc khớp tham số thực nghiệm thường bộc lộ độ bất định khó kiểm soát khi ngoại suy sang các vùng hạt nhân ngoại lai thiếu số liệu đo đạc.

Nghiên cứu này giải quyết bài toán cốt lõi: xây dựng mô hình thế quang học vi mô tự nhất quán hoàn toàn để mô tả tán xạ đàn hồi của nucleon trên hạt nhân ở dải năng lượng dưới 50 MeV mà không phụ thuộc vào các tham số điều chỉnh tùy tiện. Mục tiêu cụ thể của luận văn là phát triển thế quang học vi mô dựa trên tương tác hiệu dụng Skyrme SLy5 gồm 10 tham số chuẩn, kết hợp lý thuyết trường trung bình tự nhất quán, gần đúng pha ngẫu nhiên và cơ chế ghép hạt - dao động. Phạm vi khảo sát tập trung trên hệ thống 4 hạt nhân vỏ kín kép tiêu biểu từ nhẹ đến nặng gồm O-16, Ca-40, Ca-48 và Pb-208. Công trình mang ý nghĩa khoa học lớn khi tạo ra cầu nối trực tiếp giữa cấu trúc vi mô bên trong và các đại lượng động học phản ứng, góp phần chuẩn hóa dữ liệu hạt nhân với độ chính xác lý thuyết tiệm cận sai số thực nghiệm dưới 10 phần trăm ở vùng góc tán xạ trướ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 trên nền tảng của 3 lý thuyết vi mô chính thống: lý thuyết trường trung bình tự nhất quán Hartree-Fock, gần đúng pha ngẫu nhiên cải tiến và lý thuyết toán tử khối Green kết hợp cơ chế ghép hạt - dao động. Về mặt toán học, thế quang học vi mô được đồng nhất với toán tử khối của hàm Green một hạt theo hình thức luận phép chiếu Feshbach, bao gồm phần thực không cục bộ và phần ảo phức phụ thuộc năng lượng đại diện cho sự hấp thụ vào các kênh không đàn hồi.

Khái niệm trung tâm thứ nhất là tương tác hiệu dụng Skyrme không tầm với dạng mở rộng, tích hợp các số hạng phụ thuộc vận tốc, số hạng spin-quỹ đạo và các số hạng phụ thuộc mật độ - động lượng. Khái niệm thứ hai là gần đúng tựa boson trong không gian trạng thái kích thích hạt - lỗ, cho phép mô tả các mode dao động tập thể biên độ nhỏ như trạng thái kích thích thấp và cộng hưởng khổng lồ. Khái niệm thứ ba là hiệu chỉnh nguyên lý loại trừ Pauli thông qua tính toán thế động học bậc hai nhằm loại bỏ triệt để các trạng thái giả sinh ra từ chuyển động khối tâm.

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

Nguồn dữ liệu thực nghiệm phục vụ kiểm chứng được thu thập từ các bảng tiết diện tán xạ vi phân chuẩn quốc tế và đại lượng phân cực năng lượng phân tích. Cỡ mẫu nghiên cứu bao gồm 4 hệ hạt nhân vỏ kín kép đối xứng và bất đối xứng nơtron - proton đại diện cho 4 vùng khối lượng điển hình từ A = 16 đến A = 208. Phương pháp chọn mẫu tập trung vào các hạt nhân vỏ kín kép nhằm hạn chế tối đa ảnh hưởng phức tạp của tương tác ghép đôi hạt nhân, đảm bảo tính đối xứng cầu nghiêm ngặt trong không gian cấu hình.

Quy trình phân tích số trị được thiết lập chặt chẽ: giải hệ phương trình Hartree-Fock phi tuyến trong hộp cầu bán kính 15 fm với bước lưới không gian 0.1 fm để xác định các trạng thái đơn hạt và thế trường trung bình tĩnh. Tiếp theo, ma trận chuyển mức kích thích được thiết lập với đa cực phonon từ L = 0 đến L = 5, giới hạn năng lượng kích thích dưới 50 MeV. Phương pháp phân tích tích phân - vi phân này được lựa chọn vì tính tự nhất quán tuyệt đối, sử dụng chung một tương tác Skyrme SLy5 duy nhất cho cả trạng thái cơ bản và kích thích. Cuối cùng, phương trình Schrödinger vi phân chứa thế không cục bộ được giải thông qua việc tích hợp mã nguồn chuẩn DWBA98, hoàn thiện chu trình nghiên cứu trong vòng 36 tháng.

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

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

Thứ nhất, mô hình thế quang học vi mô tự nhất quán tái hiện thành công phân bố góc tán xạ đàn hồi của nơtron và proton trên cả 4 hạt nhân ở năng lượng dưới 50 MeV. Tại dải góc tán xạ trước dưới 60 độ, kết quả tính toán lý thuyết trùng khớp trên 95 phần trăm với số liệu thực nghiệm đo đạc được, phản ánh chính xác phần thế thực trường trung bình Hartree-Fock.

Thứ hai, việc đưa tương tác Coulomb vào cơ chế ghép hạt - dao động đã giải quyết thành công bài toán mô tả tán xạ đàn hồi proton, điều mà các mô hình vi mô trước đây thường bỏ qua. Độ hấp thụ phần ảo bề mặt được tái hiện rõ nét nhờ sự đóng góp của các mode kích thích tập thể 3- trừ và 4+ cộng, với độ cạn kiệt quy tắc tổng trọng số năng lượng đạt trên 99 phần trăm.

Thứ ba, phân tích đóng góp của các thành phần tương tác cho thấy các số hạng phụ thuộc vận tốc và tương tác spin-quỹ đạo đóng vai trò quyết định đến độ sâu của phần ảo thế quang học. Khi loại bỏ số hạng spin-quỹ đạo, năng lượng phân tích bị suy giảm hơn 60 phần trăm so với thực nghiệm.

Thứ tư, việc bổ sung các số hạng phụ thuộc mật độ - động lượng giúp điều chỉnh phân bố thế ảo ở vùng biên hạt nhân, cải thiện đáng kể độ khớp của tiết diện tán xạ tại vùng góc lùi trên 120 độ với mức cải thiện sai số định lượng đạt khoảng 25 phần trăm.

Thảo luận kết quả

Sự suy giảm tiết diện ở các góc tán xạ lớn được giải thích do cơ chế hấp thụ bề mặt mạnh mẽ khi hạt tới ghép với các trạng thái kích thích phonon tập thể. Dữ liệu nghiên cứu được biểu diễn trực quan thông qua các đồ thị phân bố góc tiết diện vi phân dạng logarit theo góc tán xạ và các bảng so sánh sai số tương đối giữa mô hình thế quang học vi mô với mô hình hiện tượng Koning-Delaroche.

Các biểu đồ hấp thụ thế ảo theo bán kính không gian chỉ ra rằng cơ chế ghép hạt - dao động tạo ra vùng hấp thụ tập trung chủ yếu tại bề mặt hạt nhân ở khoảng cách từ 3 fm đến 7 fm tùy theo khối lượng hạt nhân. So sánh với mô hình vật chất hạt nhân vô hạn sử dụng phép xấp xỉ mật độ cục bộ, mô hình vi mô cấu trúc hạt nhân thể hiện ưu thế vượt trội ở năng lượng thấp dưới 20 MeV vì đã bao hàm tường minh cấu trúc mức năng lượng rời rạc và các cộng hưởng khổng lồ của hạt nhân đích.

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

Thứ nhất, triển khai tích hợp tương tác ghép đôi Hartree-Fock-Bogoliubov vào mô hình thế quang học vi mô trong vòng 12 đến 18 tháng tới. Mục tiêu nhằm mở rộng phạm vi ứng dụng mô hình từ 4 hạt nhân vỏ kín sang hơn 200 hạt nhân mở vỏ và hạt nhân biến dạng, do các nhóm nghiên cứu vật lý lý thuyết hạt nhân chủ trì thực hiện.

Thứ hai, chuẩn hóa và xây dựng cơ sở dữ liệu thế quang học vi mô toàn cục phục vụ công nghệ tính toán phản ứng hạt nhân, hướng tới mục tiêu giảm độ bất định của tiết diện phản ứng trong thư viện RIPL xuống dưới 10 phần trăm đối với vùng hạt nhân giàu nơtron. Nhiệm vụ này cần sự phối hợp giữa các viện nghiên cứu nguyên tử quốc gia trong lộ trình 24 tháng.

Thứ three, nâng cấp thuật toán giải phương trình tích phân - vi phân không cục bộ trên nền tảng tính toán song song hiệu năng cao, tối ưu hóa thời gian xử lý dữ liệu ma trận từ vài giờ xuống dưới 15 phút cho mỗi phép tính tán xạ ở năng lượng 50 MeV.

Thứ tư, khuyến nghị các cơ quan quản lý và các phòng thí nghiệm vật lý hạt nhân ứng dụng trực tiếp bộ tham số vi mô SLy5 trong việc đánh giá an toàn bức xạ và thiết kế các bia bắn chùm ion phóng xạ trước năm 2026, thay thế dần việc dựa hoàn toàn vào các tham số hiện tượng kém tin cậy.

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

Nhóm học viên cao học và nghiên cứu sinh chuyên ngành Vật lý hạt nhân và Năng lượng nguyên tử: Luận văn cung cấp toàn bộ quy trình toán học chi tiết từ lý thuyết trường trung bình Hartree-Fock, gần đúng pha ngẫu nhiên đến phương pháp toán tử khối Green, phục vụ trực tiếp cho việc làm chủ các công cụ mô phỏng vi mô.

Nhóm kỹ sư thiết kế lò phản ứng và an toàn bức xạ: Cung cấp phương pháp đánh giá chính xác tiết diện tương tác nơtron trên các vật liệu kết cấu vỏ bọc như hợp kim chứa Canxi hay Chì-208 ở dải năng lượng dưới 50 MeV, giúp tối ưu hóa hệ số nhân nơtron và che chắn bức xạ.

Nhóm nhà thực nghiệm tại các trung tâm gia tốc chùm ion: Nắm bắt các dự đoán lý thuyết về phân bố góc và năng lượng phân tích để thiết lập cấu hình máy dò và góc đo tối ưu trong các thí nghiệm tán xạ hạt nhân ngoại lai.

Nhóm chuyên gia phát triển phần mềm tính toán khoa học: Tham khảo kiến trúc tích hợp giữa mô hình cấu trúc vi mô và mã nguồn động học DWBA98 để phát triển các gói phần mềm mã nguồn mở thế hệ mới phục vụ tính toán phản ứng hạt nhân tự nhất quán.

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

Thế quang học vi mô khác biệt như thế nào so với thế quang học hiện tượng truyền thống?

Thế hiện tượng như mô hình Koning-Delaroche sử dụng các hàm giải tích Woods-Saxon và khớp khoảng 10 đến 20 tham số tự do từ dữ liệu thực nghiệm sẵn có. Ngược lại, thế vi mô trong nghiên cứu này được xây dựng hoàn toàn từ tương tác hai hạt nucleon-nucleon cơ bản qua cấu trúc hạt nhân, có khả năng tiên đoán chính xác tại vùng chưa có số liệu đo đạc.

Tại sao nghiên cứu chỉ giới hạn năng lượng tán xạ dưới 50 MeV?

Dải năng lượng dưới 50 MeV là vùng mà các hiệu ứng cấu trúc hạt nhân, dao động tập thể và cộng hưởng khổng lồ chi phối mạnh nhất cơ chế phản ứng. Ở năng lượng này, phương pháp ghép hạt - dao động phát huy tối đa độ chính xác, trong khi trên 50 MeV các quá trình va chạm trực tiếp nhiều bước bắt đầu chiếm ưu thế.

Vai trò của việc chọn bộ tham số Skyrme SLy5 trong mô hình là gì?

Bộ tham số SLy5 gồm 10 tham số chuẩn được tối ưu hóa để mô tả đồng thời tính chất của vật chất nơtron tinh khiết và phổ đơn hạt của các hạt nhân giàu nơtron. Việc sử dụng SLy5 đảm bảo tính tự nhất quán 100 phần trăm xuyên suốt từ tính toán trạng thái cơ bản đến trạng thái kích thích RPA và tán xạ động học.

Làm thế nào để loại bỏ các trạng thái kích thích giả sinh ra trong gần đúng RPA?

Nghiên cứu áp dụng phép chiếu loại bỏ trực tiếp mode dao động lưỡng cực 1- trừ tương ứng với chuyển động khối tâm tịnh tiến của toàn bộ hệ hạt nhân. Đồng thời, các số hạng thế động học bậc hai được đưa vào để hiệu chỉnh vi phạm nguyên lý loại trừ Pauli trong gần đúng tựa boson.

Mô hình có thể mở rộng trực tiếp cho các hạt nhân biến dạng không đối xứng cầu được không?

Hoàn toàn có thể mở rộng bằng cách chuyển hệ phương trình Hartree-Fock và RPA sang hệ tọa độ trụ hoặc biến dạng trục. Tuy nhiên, không gian cấu hình ma trận sẽ tăng lên gấp nhiều lần, đòi hỏi kỹ thuật tính toán song song trên các cụm máy tính hiệu năng cao với bộ nhớ mở rộng.

Kết luận

Luận văn đã giải quyết trọn vẹn bài toán xây dựng mô hình thế quang học vi mô tự nhất quán cho tán xạ nơtron và proton năng lượng thấp dưới 50 MeV trên 4 hạt nhân vỏ kín kép tiêu biểu gồm O-16, Ca-40, Ca-48 và Pb-208. Năm đóng góp học thuật cốt lõi bao gồm:

  • Thiết lập thành công chu trình tính toán tự nhất quán hoàn chỉnh kết hợp giữa trường trung bình Hartree-Fock, gần đúng pha ngẫu nhiên và cơ chế ghép hạt - dao động với tương tác Skyrme SLy5.
  • Mô tả chính xác phân bố góc tán xạ vi phân và năng lượng phân tích của nơtron, đạt độ trùng khớp trên 95 phần trăm ở các góc tán xạ trước so với dữ liệu thực nghiệm.
  • Tích hợp thành công tương tác Coulomb vào phương trình ghép dao động để mô tả đồng thời tán xạ đàn hồi proton.
  • Định lượng rõ rệt vai trò của các số hạng spin-quỹ đạo và phụ thuộc mật độ - động lượng đối với phần ảo của thế quang học.
  • Khẳng định tính ưu việt của mô hình vi mô cấu trúc hạt nhân trong việc dự đoán dữ liệu tán xạ ở các vùng hạt nhân ngoại lai mà không cần tham số điều chỉnh tùy tiện.

Trong giai đoạn tiếp theo, nhóm tác giả định hướng mở rộng mô hình cho các hạt nhân mở vỏ có xét đến tương tác ghép đôi hạt nhân. Độc giả quan tâm và các đơn vị nghiên cứu có thể kết nối, khai thác mã nguồn tính toán và ứng dụng kết quả vào các dự án mô phỏng phản ứng hạt nhân chuyên sâu.