Hàng rào thế Coulomb ảnh hưởng như thế nào đến phản ứng kích thích proton?

Do proton mang điện tích dương cùng dấu với hạt nhân bia, lực đẩy tĩnh điện Coulomb tạo nên một hàng rào thế năng cản trở sự tiếp cận của hạt đạn. Proton phải sở hữu động năng đủ lớn để vượt qua hoặc xuyên hầm lượng tử qua hàng rào này trước khi lực hạt nhân mạnh tầm ngắn có thể khởi phát các phản ứng hạt nhân.

Phản ứng kích thích proton khác phản ứng do neutron kích thích ở điểm cơ bản nào?

Khác biệt cơ bản nhất nằm ở tương tác tĩnh điện: neutron là hạt trung hòa điện tích nên không bị hàng rào thế Coulomb ngăn cản và có thể gây ra phản ứng hạt nhân ở mức năng lượng nhiệt rất thấp, trong khi proton luôn đòi hỏi một ngưỡng động năng tối thiểu để khắc phục lực đẩy tĩnh điện của hạt nhân bia.

Tại sao phản ứng kích thích proton được ứng dụng rộng rãi trong y học hạt nhân?

Các máy gia tốc cyclotron y tế gia tốc chùm proton đến năng lượng tối ưu vài chục MeV để kích hoạt các phản ứng hạt nhân có tiết diện cao trên các tấm bia chuyên dụng, tạo ra các đồng vị phóng xạ có độ tinh khiết cao và thời gian sống phù hợp cho kỹ thuật chụp cắt lớp phát xạ positron (PET) và xạ trị điều trị ung thư.

Phản ứng kích thích proton (proton-induced reaction) là gì?

Định nghĩa và nguyên lý cơ bản của phản ứng kích thích proton

Phản ứng kích thích proton (tiếng Anh: proton-induced reaction) là phản ứng hạt nhân cảm ứng xảy ra khi hạt nhân nguyên tử bia bị bắn phá bởi chùm hạt proton mang động năng xác định từ các máy gia tốc hạt, dẫn đến quá trình tái cấu trúc hạt nhân, phát xạ các hạt thứ cấp hoặc giải phóng bức xạ điện từ. Trong vật lý hạt nhân thực nghiệm và vật lý ứng dụng, quá trình này đóng vai trò then chốt để nghiên cứu cấu trúc mức năng lượng hạt nhân, đo đạc tiết diện phản ứng và sản xuất các hạt nhân phóng xạ phục vụ khoa học và đời sống.

Khác với các phản ứng cảm ứng bằng neutron vốn là hạt trung hòa điện tích dễ dàng thâm nhập vào hạt nhân bia mà không chịu lực cản tĩnh điện, phản ứng kích thích proton chịu sự chi phối mạnh mẽ của tương tác đẩy Coulomb giữa điện tích dương của proton và điện tích hạt nhân bia. Để phản ứng hạt nhân có thể diễn ra, chùm hạt proton tới phải mang động năng đủ lớn để vượt qua hoặc xuyên hầm lượng tử qua hàng rào thế Coulomb. Độ cao của hàng rào thế Coulomb được mô tả thông qua biểu thức vật lý:

V_C = \frac{Z_1 Z_2 e^2}{4 \pi \varepsilon_0 R}

Trong biểu thức trên, đại lượng $V_C$ biểu diễn thế năng hàng rào Coulomb, đại lượng $Z_1$ là số proton của hạt tới mang giá trị đơn vị, đại lượng $Z_2$ là số hiệu nguyên tử của hạt nhân bia, đại lượng $e$ là điện tích nguyên tố, đại lượng $\varepsilon_0$ là hằng số điện môi chân không, và đại lượng $R$ đại diện cho bán kính tiếp xúc hiệu dụng giữa hạt proton và hạt nhân bia. Khi động năng của hạt proton thấp hơn đáng kể so với hàng rào thế, xác suất xảy ra phản ứng hạt nhân suy giảm rõ rệt theo quy luật xuyên hầm hiệu ứng Gamow; ngược lại, khi năng lượng proton vượt qua ngưỡng rào cản thế, lực hạt nhân mạnh tầm ngắn sẽ chiếm ưu thế và khởi động các kênh phản ứng hạt nhân đa dạng.

Phân loại các cơ chế phản ứng hạt nhân do proton kích thích

Tùy thuộc vào dải động năng của hạt proton tới và thời gian tương tác đặc trưng giữa hạt đạn và hạt nhân bia, phản ứng kích thích proton diễn ra theo ba cơ chế vật lý cơ bản: cơ chế hình thành nhân hợp phần, cơ chế phát xạ tiền cân bằng, và cơ chế phản ứng trực tiếp.

1. Cơ chế tạo nhân hợp phần (Compound Nucleus Mechanism)

Ở dải động năng thấp và trung bình thấp ngay phía trên hàng rào thế Coulomb, hạt proton khi xâm nhập vào vùng tác dụng của lực hạt nhân sẽ va chạm liên tiếp với các nucleon bên trong hạt nhân bia, truyền toàn bộ năng lượng của mình cho toàn thể hệ nucleon. Quá trình này hình thành nên một trạng thái kích thích cao độ gọi là nhân hợp phần (compound nucleus). Trạng thái nhân hợp phần có thời gian sống tương đối dài so với thời gian bay xuyên qua hạt nhân của một nucleon tự do. Trong suốt thời gian tồn tại này, năng lượng kích thích được phân phối ngẫu nhiên giữa các nucleon thông qua vô số va chạm nội hạt nhân cho đến khi đạt trạng thái cân bằng nhiệt động thống kê.

Sự phân rã của nhân hợp phần hoàn toàn độc lập với phương thức hình thành ban đầu của nó, tuân theo định đề độc lập Bohr và được tính toán lượng hóa thông qua lý thuyết thống kê Hauser-Feshbach. Khi một nucleon hoặc cụm nucleon ngẫu nhiên tập trung đủ năng lượng vượt quá năng lượng liên kết, nhân hợp phần sẽ khử kích thích thông qua việc bay hơi các hạt nhẹ như neutron, proton hoặc hạt alpha, đồng thời giải phóng phần năng lượng kích thích còn lại dưới dạng các bước chuyển bức xạ gamma về trạng thái cơ bản.

2. Cơ chế phát xạ tiền cân bằng (Pre-equilibrium Emission)

Khi năng lượng chùm proton tăng lên dải trung bình, quá trình tương tác nội hạt nhân diễn ra quá nhanh khiến hệ nucleon chưa kịp đạt tới trạng thái cân bằng thống kê hoàn toàn trước khi phát xạ hạt. Cơ chế này được gọi là sự phát xạ tiền cân bằng. Về mặt lý thuyết, quá trình tiến tới cân bằng được xem như một chuỗi các trạng thái trung gian kế tiếp nhau, trong đó số lượng hạt và lỗ trống (gọi chung là exciton) tăng dần qua từng bước va chạm nhị phân.

Nền tảng vi mô hoàn chỉnh để mô tả hiện tượng này được xác lập bởi Herman Feshbach, Arthur Kerman và Steven Koonin vào năm 1980 thông qua lý thuyết cơ học lượng tử đa bước FKK. Lý thuyết FKK phân định rõ ràng quá trình phản ứng thành hai thành phần cơ bản: tán xạ hợp phần đa bước (multistep compound - MSC) tương ứng với các trạng thái trung gian bị giam giữ trong giếng thế hạt nhân, và tán xạ trực tiếp đa bước (multistep direct - MSD) tương ứng với các trạng thái liên tục nơi ít nhất một nucleon bảo lưu động năng lớn để thoát ra ngoài trước khi cân bằng nhiệt được thiết lập. Sự xuất hiện của các hạt phát xạ tiền cân bằng giải thích chính xác phần phổ năng lượng cứng và tính bất đối xứng về hướng phát xạ phía trước quan sát được trong thực nghiệm phản ứng kích thích proton.

3. Cơ chế phản ứng trực tiếp và phản ứng vỡ vụn hạt nhân (Direct Reactions and Spallation)

Ở dải năng lượng cao, bước sóng de Broglie của hạt proton trở nên rất ngắn so với kích thước hạt nhân bia, khiến proton không còn tương tác với toàn bộ hạt nhân như một khối thống nhất mà va chạm trực tiếp với từng nucleon riêng lẻ trên bề mặt hoặc bên trong thể tích hạt nhân. Quá trình này diễn ra trong khoảng thời gian cực ngắn tương đương thời gian bay của proton qua đường kính hạt nhân. Các kênh phản ứng trực tiếp điển hình bao gồm tán xạ không đàn hồi kích thích các trạng thái tập thể, phản ứng bứt nucleon (stripping), phản ứng bắt nucleon (pick-up), hoặc phản ứng đánh bật nucleon (knock-out).

Khi năng lượng proton vượt lên trên dải hàng trăm MeV, phản ứng bước vào giai đoạn vỡ vụn hạt nhân (spallation). Hạt proton sơ cấp năng lượng cao khởi phát một chuỗi va chạm thác nội hạt nhân (intranuclear cascade), đẩy bật liên tiếp các nucleon thứ cấp năng lượng cao ra khỏi bia, để lại một hạt nhân tàn dư kích thích cao độ tiếp tục phân rã thông qua quá trình bay hơi hạt hoặc vỡ vụn thành các mảnh hạt nhân kích thước trung bình.

Mô hình quang học và hệ thống đánh giá dữ liệu hạt nhân

Để tính toán chính xác tiết diện phản ứng tổng cộng, tiết diện hấp thụ và sự phân bố góc của các hạt phát xạ trong phản ứng kích thích proton, các nhà vật lý lý thuyết sử dụng mô hình quang học hạt nhân (optical model). Trong mô hình này, tương tác phức tạp giữa proton tới và hàng chục nucleon bia được quy đổi tương đương về bài toán tán xạ một hạt trên một thế năng phức:

U(r) = -V_R f_R(r) - i W_V f_V(r) + 4 i a_D W_D \frac{d f_D(r)}{dr} + V_{so} \frac{1}{r} \frac{d f_{so}(r)}{dr} (\mathbf{L} \cdot \mathbf{\sigma}) + V_C(r)

Trong phương trình thế quang học, thành phần thực $V_R$ thể hiện lực hút hạt nhân trung bình gây ra hiện tượng tán xạ đàn hồi, các thành phần ảo thể tích $W_V$ và ảo bề mặt $W_D$ mô tả sự suy giảm dòng proton do bị hấp thụ vào các kênh phản ứng không đàn hồi, thành phần spin - quỹ đạo $V_{so}$ liên kết mômen quỹ đạo $\mathbf{L}$ với spin hạt $\mathbf{\sigma}$ , và thành phần $V_C(r)$ là thế tĩnh điện Coulomb của khối điện tích phân bố.

Bước đột phá trong việc chuẩn hóa các thông số thế quang học được thiết lập bởi Arjan Koning và Jean-Paul Delaroche vào năm 2003. Công trình của Koning và Delaroche (2003) đã xây dựng bộ thông số thế quang học hiện tượng luận toàn cục áp dụng cho cả proton và neutron trên dải năng lượng rộng từ 1 keV đến 200 MeV, bao quát hệ thống hạt nhân bia dạng cầu và gần cầu có số khối từ 24 đến 209. Bộ thông số toàn cục này loại bỏ sự phụ thuộc cục bộ của từng phép đo riêng lẻ, mang lại độ chính xác cao trong việc dự báo tiết diện phản ứng kích thích proton trên các đồng vị chưa có số liệu thực nghiệm.

Trên cơ sở mô hình quang học hiện đại và các lý thuyết phân rã thống kê, hệ thống phần mềm tính toán phản ứng hạt nhân TALYS do Arjan Koning và Dimitri Rochman công bố năm 2012 đã trở thành công cụ tiêu chuẩn toàn cầu trong việc đánh giá và hoàn thiện dữ liệu hạt nhân. Hệ thống TALYS tích hợp đồng bộ các mô hình tán xạ trực tiếp, mô hình tiền cân bằng exciton và mô hình nhân hợp phần Hauser-Feshbach, cho phép mô phỏng toàn diện tất cả các kênh phản ứng kích thích bởi proton từ mức năng lượng 1 keV đến 200 MeV.

Ứng dụng thực tiễn của phản ứng kích thích proton

Nhờ khả năng biến đổi hạt nhân và kích thích đặc trưng, phản ứng do proton kích thích được ứng dụng rộng rãi trong nhiều lĩnh vực công nghệ cao, y học bức xạ và phân tích vật liệu.

1. Sản xuất đồng vị phóng xạ y tế

Ứng dụng quan trọng và phổ biến nhất của phản ứng kích thích proton trong đời sống là sản xuất các hạt nhân phóng xạ dùng trong y học hạt nhân tại các máy gia tốc hạt cyclotron bệnh viện và viện nghiên cứu. Nhờ tiết diện phản ứng lớn ở dải năng lượng vài chục MeV, chùm proton bắn phá các tấm bia chuyên dụng để tạo ra các đồng vị có thời gian sống phù hợp cho chẩn đoán hình ảnh và xạ trị đích.

Để đảm bảo hiệu suất tạo đồng vị tối ưu và kiểm soát tạp chất phóng xạ đi kèm, Cơ quan Năng lượng Nguyên tử Quốc tế IAEA đã liên tục thẩm định và công bố các tập dữ liệu hạt nhân chuẩn hóa. Năm 2022, Ferenc Tárkányi cùng các chuyên gia quốc tế trong khuôn khổ dự án điều phối của IAEA đã hoàn thành bản nâng cấp toàn diện cơ sở dữ liệu tiết diện phản ứng khuyến nghị cho việc sản xuất các hạt nhân phóng xạ điều trị bằng chùm hạt mang điện, bao gồm các phản ứng kích thích proton quan trọng. Các dữ liệu hàm kích thích chuẩn hóa này giúp các cơ sở y tế xác định chính xác dải năng lượng chùm proton chiếu xạ để tối đa hóa sản lượng hạt nhân đích và hạn chế tối đa các đồng vị phụ không mong muốn.

2. Kỹ thuật phân tích phát xạ tia X kích thích bằng proton (PIXE)

Bên cạnh các tương tác hạt nhân sâu, chùm proton năng lượng khoảng vài MeV khi bắn phá vật chất còn gây ra tương tác ion hóa mạnh mẽ với các electron thuộc các lớp vỏ trong của nguyên tử bia. Khi các electron lớp vỏ trong bị đánh bật ra ngoài, các electron ở phân lớp ngoài sẽ nhảy xuống lấp chỗ trống, giải phóng năng lượng dưới dạng bức xạ tia X đặc trưng. Phương pháp phân tích phát xạ tia X kích thích bằng proton (PIXE) là một kỹ thuật phân tích hạt nhân không phá hủy có độ nhạy cực cao, cho phép xác định đồng thời hàm lượng vết của nhiều nguyên tố hóa học trong mẫu khảo sát môi trường, địa chất, sinh học và cổ vật văn hóa.

Nghiên cứu thực nghiệm tại Việt Nam

Tại Việt Nam, các hướng nghiên cứu liên quan đến phản ứng kích thích proton đã và đang được thúc đẩy thông qua các chương trình hợp tác quốc tế về vật lý hạt nhân và kỹ thuật máy gia tốc. Các nhà khoa học trong nước tập trung vào các hướng nghiên cứu trọng điểm: ứng dụng kỹ thuật phân tích chùm hạt ion và nghiên cứu thực nghiệm đo đạc hàm kích thích phản ứng hạt nhân phục vụ dữ liệu hạt nhân ứng dụng.

Về nghiên cứu thực nghiệm tiết diện phản ứng, nhóm nghiên cứu do Nguyễn Văn Đỗ dẫn đầu cùng các cộng sự năm 2021 đã công bố công trình đo đạc thực nghiệm hàm kích thích tạo các đồng vị phóng xạ niobium (bao gồm niobium 90, niobium 91m và niobium 92m) thông qua phản ứng kích thích proton trên bia zirconium tự nhiên. Nghiên cứu thực nghiệm được tiến hành trong dải năng lượng chùm hạt proton từ 10,58 MeV đến 43,61 MeV bằng phương pháp kích hoạt lá bia xếp lớp kết hợp phổ kế gamma bán dẫn siêu tinh khiết. Kết quả đo đạc thực nghiệm của Nguyễn Văn Đỗ cùng các cộng sự (2021) đã cung cấp dữ liệu tiết diện chính xác, đồng thời được đối sánh chi tiết với các phép tính toán lý thuyết từ hệ chương trình mô phỏng TALYS, đóng góp thiết thực vào cơ sở dữ liệu hạt nhân quốc tế phục vụ các thiết kế bia bắn phá cyclotron.

Bảng so sánh các cơ chế phản ứng kích thích proton

Bảng dưới đây tổng hợp các đặc trưng vật lý cơ bản phân biệt giữa ba cơ chế tương tác chính của hạt proton với hạt nhân nguyên tử:

Cơ chế phản ứng	Dải năng lượng đặc trưng	Thời gian tương tác	Mô hình lý thuyết chủ đạo	Đặc điểm hạt phát xạ thứ cấp
Nhân hợp phần (Compound Nucleus)	Năng lượng thấp, gần hàng rào Coulomb	Kéo dài, khoảng mười lũy thừa âm mười sáu đến âm mười tám giây	Lý thuyết thống kê Hauser-Feshbach	Phát xạ đẳng hướng, phân bố năng lượng theo quy luật bay hơi nhiệt Maxwell
Tiền cân bằng (Pre-equilibrium)	Năng lượng trung bình	Trung gian, khoảng mười lũy thừa âm hai mươi mốt giây	Mô hình exciton, lý thuyết đa bước FKK (1980)	Phát xạ ưu tiên góc phía trước, phổ năng lượng cứng hơn phân bố nhiệt
Phản ứng trực tiếp (Direct Reactions)	Năng lượng cao	Cực ngắn, khoảng mười lũy thừa âm hai mươi hai giây	Mô hình sóng méo DWBA, mô hình quang học Koning - Delaroche (2003)	Bảo lưu mạnh tính chất trạng thái hạt đạn, đỉnh tán xạ góc hẹp phía trước

Hạn chế thực nghiệm và triển vọng nghiên cứu

Mặc dù lý thuyết và thực nghiệm về phản ứng kích thích proton đã đạt được nhiều thành tựu vững chắc, các nhà nghiên cứu vẫn đối mặt với một số rào cản kỹ thuật và bất định lý thuyết đáng kể:

Hiệu ứng tỏa nhiệt trên bia bắn phá: Khi vận hành máy gia tốc với chùm proton cường độ cao để thu nhận sản lượng đồng vị lớn, sự mất mát năng lượng do ion hóa Coulomb sinh ra mật độ nhiệt khổng lồ trên bia kim loại, đòi hỏi các hệ thống làm mát tản nhiệt phức tạp để ngăn bia bị nóng chảy hoặc biến dạng cơ học.
Bất định của mật độ mức hạt nhân: Trong các mô hình thống kê Hauser-Feshbach và TALYS, độ chính xác của tiết diện tính toán phụ thuộc rất lớn vào các tham số mật độ mức hạt nhân và hàm lực bức xạ gamma ở trạng thái kích thích cao, vốn vẫn còn độ lệch nhất định đối với các đồng vị nằm xa dải ổn định.
Nhu cầu dữ liệu tiết diện năng lượng cao: Với sự phát triển của các nguồn neutron phân rã vỡ vụn và các lò phản ứng hạt nhân thế hệ mới điều khiển bằng máy gia tốc (Accelerator Driven Systems - ADS), việc mở rộng đo đạc thực nghiệm phản ứng kích thích proton lên dải năng lượng trên hàng trăm MeV đối với các vật liệu kết cấu và nhiên liệu hạt nhân đang là trọng tâm nghiên cứu sôi động trên thế giới.