Tổng quan nghiên cứu

Thí nghiệm vật lý hạt nhân Belle II tại Tổ chức Nghiên cứu Máy gia tốc Năng lượng cao (KEK), Tsukuba, Nhật Bản đánh dấu bước tiến vượt bậc trong việc tìm kiếm các hiện tượng vật lý mới ngoài Mô hình Chuẩn. Với việc nâng cấp máy gia tốc KEKB thành SuperKEKB, độ sáng thiết kế (Luminosity) đạt mức $8 \times 10^{35} \text{ cm}^{-2}\text{s}^{-1}$, cao gấp 40 lần so với phiên bản tiền nhiệm. Cùng với sự gia tăng độ sáng, lượng meson B tạo thành từ quá trình va chạm giữa chùm hạt electron 7 GeV và positron 4 GeV ở vùng năng lượng cộng hưởng $1.21 \times 10^{-20}$ giây của hạt Upsilon tăng lên đáng kể, cho phép ghi nhận các phân rã hiếm với độ chính xác vị trí phân rã đạt khoảng 200 micromet xung quanh ống va chạm beryllium đường kính 1 cm.

Tuy nhiên, cường độ chùm tia cao đồng nghĩa với việc mức độ bức xạ phông do chùm tia tạo ra tăng vọt từ 10 đến 20 lần. Bức xạ neutron sinh ra từ các quá trình tán xạ phông này đe dọa trực tiếp đến tuổi thọ và độ ổn định của detector buồng trôi trung tâm (Central Drift Chamber - CDC), đặc biệt là hệ thống bảng mạch điện tử đọc tín hiệu (Front-End Electronics) gồm 14.336 kênh truyền dẫn. Mục tiêu trọng tâm của nghiên cứu là mô phỏng, định lượng thông lượng neutron tác động lên detector CDC, xác định nguồn gốc không gian sinh ra hạt neutron và xây dựng cấu trúc hình học che chắn tối ưu. Nghiên cứu được thực hiện trong giai đoạn chuẩn bị vận hành từ năm 2013 đến năm 2014, mang ý nghĩa quyết định nhằm bảo vệ hệ thống thiết bị đo đạc, kiểm soát sai số đo lường góc của tam giác unita ở mức dưới 10%, đồng thời đánh dấu sự đóng góp học thuật của nhóm Vật lý năng lượng cao Việt Nam trong dự án hạt nhân quy mô quốc tế.

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

Khung lý thuyết áp dụng

Cơ sở lý thuyết của luận văn dựa trên Mô hình Chuẩn của vật lý hạt cơ bản và cơ chế vi phạm đối xứng CP (Charge Conjugation - Parity). Sự vi phạm đối xứng CP trong phân rã của meson B được mô tả thông qua ma trận CKM (Cabibbo-Kobayashi-Maskawa) gồm 3 thế hệ quark. Điều kiện trực giao của ma trận tạo nên tam giác unita db với phương trình $V_{ud}V_{ub}^* + V_{cd}V_{cb}^* + V_{td}V_{tb}^* = 0$, trong đó các góc của tam giác phản ánh mức độ bất đối xứng vật chất - phản vật chất trong vũ trụ.

Quá trình hình thành phông bức xạ trong máy gia tốc SuperKEKB được giải thích qua ba hiện tượng vật lý chủ đạo:

  • Hiệu ứng Touschek: Tán xạ đàn hồi giữa các hạt mang điện trong cùng một bó chùm tia chuyển động, làm biến đổi xung lượng ngang thành xung lượng dọc, khiến hạt bị lệch quỹ đạo và va chạm vào thành ống chân không. Tốc độ tán xạ Touschek tỷ lệ thuận với bình phương mật độ hạt trong bó và tỷ lệ nghịch với năng lượng chùm tia.
  • Tán xạ Bhabha phát xạ photon: Tương tác tán xạ giữa electron và positron kèm theo phát bức xạ hãm năng lượng cao ($e^+e^- \to e^+e^-\gamma$). Tia gamma góc hẹp này va chạm vào thành thiết bị và kích hoạt các phản ứng quang hạt nhân $(\gamma, n)$ tạo ra chùm neutron thứ cấp.
  • Tán xạ chùm tia với khí dư (Beam-gas): Tán xạ Coulomb và bức xạ hãm của chùm hạt với các phân tử khí dư như hydro và carbon monoxide tồn tại trong đường ống ở áp suất khoảng $10^{-7}\text{ Pa}$.

Các neutron tự do sinh ra từ các quá trình trên tương tác với cấu trúc bán dẫn CMOS và vi mạch FPGA trên hệ bảng mạch của detector CDC thông qua hai cơ chế chính: sự phá hủy dịch chuyển mạng tinh thể (displacement damage) và hiệu ứng liều ion hóa tích lũy (Total Ionizing Dose - TID).

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

Nghiên cứu sử dụng nguồn dữ liệu mô phỏng từ cấu hình máy gia tốc SuperKEKB với chùm tia electron 7 GeV tại vòng năng lượng cao (HER) và chùm positron 4 GeV tại vòng năng lượng thấp (LER), va chạm chéo nhau ở góc 41 mrad. Cỡ mẫu nghiên cứu bao gồm hàng triệu sự kiện tán xạ hạt được lấy mẫu ngẫu nhiên theo phương pháp Monte Carlo phân bố không gian và động lượng, mô phỏng tương đương thời gian vận hành liên tục trong 1 năm của tổ hợp máy gia tốc.

Phương pháp phân tích dựa trên khung phần mềm chuyên dụng Basf2 (Belle II Analysis Framework) tích hợp công cụ mô phỏng tương tác hạt Geant4. Các mô-đun chuyên biệt như BBBrem được sử dụng để phát sinh sự kiện bức xạ Bhabha, trong khi phần mềm SAD (Strategic Accelerator Design) tính toán quỹ đạo các hạt lệch chùm do hiệu ứng Touschek và tán xạ khí dư. Lý do lựa chọn tổ hợp phần mềm Basf2 và Geant4 là khả năng tái hiện chi tiết hình học 3D phức tạp của toàn bộ 51.456 sợi dây điện cực và các khối vật liệu che chắn trong buồng CDC, cho phép theo dõi vết (tracking) của từng hạt neutron từ vị trí sinh ra cho đến khi truyền năng lượng vào bảng mạch. Toàn bộ chuỗi tính toán và phân tích dữ liệu được hoàn thành theo mốc thời gian từ tháng 6 năm 2013 đến đầu năm 2014.

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

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

Quá trình mô phỏng và phân tích dữ liệu trên phần mềm Basf2 đã mang lại bốn phát hiện quan trọng:

  • Phân bố thông lượng neutron trên hệ bảng mạch: Mật độ neutron tập trung cao nhất tại các khu vực rìa ngoài và nắp đầu cuối (End-Cap) của detector CDC. Tổng thông lượng neutron chiếu xạ lên hệ thống bảng mạch Front-End trong 1 năm vận hành vượt mức $10^{11} \text{ neutron/cm}^2$, gây nguy cơ cao dẫn đến các sự cố lật bit dữ liệu (Single Event Upset) trong vi mạch FPGA.
  • Nguồn gốc phát sinh hạt neutron: Bức xạ Bhabha và hiệu ứng Touschek là hai tác nhân chủ yếu, đóng góp tới hơn 85% tổng lượng neutron sinh ra trong toàn bộ hệ thống. Trong đó, các photon năng lượng cao từ bức xạ Bhabha va đập vào bề mặt các nam châm siêu dẫn hội tụ cuối (Final Focusing Magnets) là điểm phát xạ quang neutron lớn nhất.
  • Tốc độ đếm phông (Hit rate) trong buồng trôi: Mức độ chiếm dụng tín hiệu (occupancy) của detector CDC tăng khoảng 15 đến 20 lần so với thí nghiệm Belle trước đây. Với 14.336 kênh đọc tín hiệu, sự gia tăng này có thể làm giảm hiệu suất tái tạo quỹ đạo hạt tích điện nếu không có biện pháp triệt tiêu phông nền.
  • Hiệu quả của vật liệu suy giảm bức xạ: Việc thử nghiệm các cấu hình vật liệu khác nhau cho thấy các hợp chất giàu hydro kết hợp nguyên tố có tiết diện bắt neutron lớn có khả năng giảm từ 60% đến hơn 70% mật độ dòng neutron truyền tới khu vực bảng mạch điện tử 48 kênh.

Thảo luận kết quả

Sự gia tăng đột biến của thông lượng neutron bắt nguồn từ việc kích thước bó hạt theo phương đứng tại điểm tương tác bị thu hẹp 20 lần kết hợp cùng cường độ dòng điện chùm tia tăng gấp 2 lần (đạt 3,6 A ở vòng LER và 2,6 A ở vòng HER). Mật độ hạt dày đặc làm tăng xác suất va chạm nội bó Touschek và bức xạ Bhabha theo hàm mũ.

Dữ liệu mô phỏng được trực quan hóa sinh động thông qua các bản đồ nhiệt 3D thể hiện mật độ phân bố hạt trong không gian buồng trôi và các biểu đồ phân bố năng lượng neutron theo trục dọc z. Các bảng số liệu so sánh thông lượng neutron trước và sau khi bổ sung các lớp vật liệu che chắn đã định lượng chính xác độ suy giảm liều chiếu, cho thấy dạng suy giảm hàm mũ theo bề dày lớp chắn.

Khi so sánh với các nghiên cứu phông bức xạ tại máy gia tốc LHCb thuộc CERN, kết quả nghiên cứu hoàn toàn tương thích về mặt quy luật phân bố quang neutron trong các tương tác va chạm lepton năng lượng cao. Ý nghĩa thực tiễn của phát hiện này là cung cấp cơ sở dữ liệu vật lý chính xác để các kỹ sư thiết kế lại vỏ bọc cơ khí và tối ưu hóa hệ thống làm mát cho detector trước thời điểm lắp đặt thực tế.

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

Nhằm bảo vệ hệ thống thiết bị và đảm bảo tính toàn vẹn của dữ liệu thực nghiệm, luận văn đưa ra bốn nhóm giải pháp hành động cụ thể:

  • Lắp đặt hệ thống tấm chắn polyetylen borat hóa: Tiến hành chế tạo và lắp ráp các tấm chắn polyetylen chứa 5% nguyên tố Bo với độ dày từ 5 cm đến 10 cm bao bọc trực tiếp quanh cụm bảng mạch Front-End của detector CDC. Mục tiêu nhằm làm chậm và hấp thụ neutron nhiệt, giảm 65% đến 75% thông lượng neutron tại các vi mạch nhạy cảm trong vòng 6 tháng trước khi đóng nắp buồng đo. Đơn vị thực hiện: Nhóm kỹ thuật cơ khí KEK phối hợp cùng nhóm nghiên cứu Belle II Việt Nam.
  • Bổ sung lớp chắn chì tại vùng tương tác và nam châm hội tụ: Bố trí các khối chì tiêu chuẩn có bề dày 2 cm đến 4 cm tại các vị trí chùm tia lệch va đập nhiều nhất gần ống va chạm beryllium nhằm hấp thụ tia gamma cứng, cắt đứt phản ứng $(\gamma, n)$ ngay tại nguồn phát trong quý I năm 2015. Chủ thể thực hiện: Bộ phận An toàn Bức xạ SuperKEKB.
  • Nâng cấp firmware và thiết kế dự phòng trên chip FPGA: Triển khai kiến trúc dự phòng bộ ba mô-đun (Triple Modular Redundancy - TMR) trên toàn bộ 14.336 kênh xử lý số liệu của bảng mạch để tự động sửa lỗi lật bit nhớ do neutron gây ra, hoàn thành quy trình nạp và kiểm tra code mô phỏng trong vòng 90 ngày. Đơn vị phụ trách: Nhóm Kỹ thuật Điện tử CDC.
  • Tăng cường kiểm soát độ chân không đường ống dẫn chùm: Duy trì áp suất buồng chân không ở mức dưới $10^{-7}\text{ Pa}$ thông qua hệ thống bơm hút ion định kỳ, giảm thiểu 30% lượng hạt tán xạ với khí dư CO và H2 trước mỗi chu kỳ vận hành máy gia tốc. Trách nhiệm thực hiện: Ban Quản lý Vận hành Máy gia tốc KEK.

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

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị học thuật và ứng dụng sâu rộng cho 4 nhóm đối tượng chính:

  • Nhà nghiên cứu và chuyên gia vật lý năng lượng cao: Cung cấp bức tranh toàn diện về phương pháp mô phỏng phông bức xạ trên các máy gia tốc thế hệ mới, ứng dụng hiệu quả khung phần mềm Basf2 kết hợp Geant4 trong các bài toán phân tích va chạm hạt.
  • Kỹ sư thiết kế điện tử hạt nhân và hệ thống đo lường: Tiếp cận dữ liệu thực nghiệm về cơ chế phá hủy bức xạ trên linh kiện bán dẫn CMOS và FPGA, làm cơ sở xây dựng tiêu chuẩn kháng bức xạ cho các khối xử lý tín hiệu Front-End trong môi trường khắc nghiệt.
  • Giảng viên và học viên cao học chuyên ngành Vật lý nguyên tử: Sử dụng làm tài liệu tham khảo chuyên khảo về lý thuyết vi phạm đối xứng CP, ma trận pha trộn quark CKM và nguyên lý vận hành của buồng trôi trung tâm CDC 56 lớp dây.
  • Chuyên gia an toàn bức xạ và công nghệ vật liệu che chắn: Khai thác các mô hình tính toán tán xạ neutron và giải pháp kết hợp vật liệu composite borat hóa để ứng dụng trong che chắn an toàn cho các lò phản ứng nghiên cứu và thiết bị xạ trị y tế hạt nhân.

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

Tại sao bức xạ neutron lại gây nguy hiểm nghiêm trọng cho hệ điện tử của detector CDC?
Neutron là hạt trung hòa điện tích nên có khả năng đâm xuyên sâu vào cấu trúc bán dẫn. Khi va chạm với mạng tinh thể silicon, neutron gây ra hiện tượng dịch chuyển nguyên tử (displacement damage) và tạo ra các xung điện tích tức thời dẫn đến lỗi lật bit bộ nhớ (Single Event Upset) trong vi mạch FPGA, đe dọa làm tê liệt 14.336 kênh đọc tín hiệu.

Nguồn gốc chính sinh ra neutron trong thí nghiệm Belle II là gì?
Neutron chủ yếu sinh ra từ các phản ứng quang hạt nhân $(\gamma, n)$ khi tia gamma năng lượng cao phát ra từ quá trình tán xạ Bhabha và bức xạ hãm va chạm với thành ống chùm tia hoặc khối nam châm kim loại nặng, cùng với sự phân rã của các hạt thứ cấp sau tán xạ Touschek.

Đặc điểm cấu tạo nổi bật của detector CDC trong Belle II là gì?
Detector CDC sử dụng hỗn hợp khí gồm 50% Heli và 50% Ethane ($C_2H_6$), bao gồm 51.456 sợi dây điện cực chia thành 56 lớp và 9 siêu lớp sắp xếp đan xen giữa phương dọc và phương nghiêng (AUAVAUAVA), giúp đo đạc chính xác quỹ đạo hạt tích điện và phân biệt loại hạt qua độ suy giảm năng lượng.

Làm thế nào để mô phỏng chính xác phông chùm tia trong phần mềm Basf2?
Hệ thống Basf2 kết hợp các mô-đun sinh sự kiện như BBBrem cho bức xạ Bhabha, SAD cho hiệu ứng Touschek và tán xạ khí dư, sau đó truyền thông tin vị trí và động lượng hạt vào nhân Geant4 để mô phỏng tương tác vật lý chi tiết của từng hạt xuyên qua các lớp vật liệu của máy dò.

Giải pháp che chắn nào đem lại hiệu quả giảm phông neutron tối ưu nhất?
Sự kết hợp giữa vật liệu polyetylen borat hóa chứa 5% Bo có bề dày 5 đến 10 cm để làm chậm và hấp thụ neutron nhiệt, cùng với các tấm chì dày 2 đến 4 cm để hấp thụ bức xạ gamma sơ cấp, giúp giảm hơn 70% thông lượng hạt bức xạ đến bảng mạch điện tử.

Kết luận

  • Luận văn đã đánh giá toàn diện cơ chế phát sinh và mức độ gia tăng từ 10 đến 20 lần của phông chùm tia trong thí nghiệm Belle II khi nâng cấp độ sáng lên mức $8 \times 10^{35} \text{ cm}^{-2}\text{s}^{-1}$.
  • Xác định định lượng thông lượng neutron và tác động phá hủy bức xạ lên hệ thống 14.336 kênh điện tử Front-End của buồng trôi trung tâm CDC.
  • Ứng dụng thành công hệ công cụ mô phỏng Basf2 kết hợp Geant4 để tái hiện không gian tương tác bức xạ phức tạp trong cấu trúc máy gia tốc SuperKEKB.
  • Đề xuất phương án thiết kế hình học che chắn đa tầng bằng polyetylen borat hóa và chì, giúp suy giảm hơn 70% mật độ neutron tại vị trí linh kiện nhạy cảm.
  • Khẳng định năng lực hội nhập khoa học quốc tế của đội ngũ vật lý năng lượng cao Việt Nam trong việc đóng góp trực tiếp vào các dự án nghiên cứu hạt nhân hàng đầu thế giới.

Bước tiếp theo của hướng nghiên cứu là tiến hành đo đạc đối chứng thực tế khi chùm tia SuperKEKB chính thức hoạt động trong giai đoạn 2015-2016 nhằm tinh chỉnh mô hình mô phỏng. Các nhà nghiên cứu và kỹ sư quan tâm có thể ứng dụng khung phương pháp luận của đề tài để phát triển các cấu trúc cảm biến bức xạ tiên tiến trong tương lai.