Tổng quan luận án

Bất biến tích phân điện tích - chẵn lẻ (CP) và bất biến tích điện tích - chẵn lẻ - đảo thời gian (CPT) là các nguyên lý đối xứng rời rạc trong lý thuyết trường lượng tử và vật lý hạt cơ bản. Việc kiểm chứng thực nghiệm các đối xứng này có ý nghĩa trong việc làm sáng tỏ bản chất của các tương tác cơ bản cũng như giải thích sự mất cân xứng vật chất - phản vật chất trong vũ trụ. Hiện tượng dao động neutrino, vốn là minh chứng thực nghiệm duy nhất về vật lý vượt ra ngoài Mô hình Chuẩn (Standard Model - SM) được xác nhận cho đến nay, cung cấp môi trường để đo đạc vi phạm đối xứng CP trong khu vực lepton và kiểm chứng độ chính xác của định lý CPT.

Luận án tiến sĩ vật lý của tác giả Trần Văn Ngọc tập trung giải quyết các khoảng trống nghiên cứu chính sau:

  1. Độ nhạy thống kê trong việc xác định pha vi phạm CP ($\delta_{CP}$) trong khu vực lepton ở các thí nghiệm hiện hành chưa đạt mức khám phá tiêu chuẩn ($5\sigma$), đồng thời bị ảnh hưởng bởi sự thoái biến tham số (parameter degeneracy) giữa pha $\delta_{CP}$, trật tự khối lượng neutrino (mass ordering/hierarchy) và góc phần tám của góc trộn $\theta_{23}$ (octant degeneracy).
  2. Kiểm chứng đối xứng CPT thông qua việc so sánh các tham số dao động của neutrino và phản neutrino ($\Delta m^2_{31}$ so với $\Delta \bar{m}^2_{31}$, $\sin^2\theta_{23}$ so với $\sin^2\bar{\theta}_{23}$) cần được khảo sát định lượng chi tiết với tập dữ liệu nâng tích lũy của thí nghiệm Tokai to Kamioka (T2K), đồng thời đánh giá tiềm năng kết hợp giữa các thí nghiệm thế hệ hiện tại và mở rộng.
  3. Kiểm soát các nguồn sai số hệ thống lớn từ thông lượng chùm neutrino và mô hình tương tác neutrino - hạt nhân thông qua hệ thống các detector gần (near detectors), bao gồm Interactive Neutrino GRID (INGRID) và WAter Grid SCIntillator Detector – prototype Magnetized Iron Neutrino Detector (WAGASCI-BabyMIND).

Mục tiêu nghiên cứu của luận án được xác định cụ thể:

  1. Thực hiện mô phỏng Monte Carlo và phân tích số liệu thực nghiệm nhằm giám sát tỷ lệ biến cố, hướng chùm và biên dạng chùm neutrino tại detector gần đặt thẳng trục INGRID dưới các cấu hình dòng điện nam thấu kính từ trường (horn) $250\text{ kA}$ và $320\text{ kA}$.
  2. Phân tích dữ liệu thực nghiệm của thí nghiệm T2K (tích lũy đến Run 9 và Run 10 với độ phơi chùm $3{,}13 \times 10^{21}\text{ POT}$) để đo đạc pha vi phạm CP Dirac $\delta_{CP}$, đồng thời kiểm tra giả thuyết bảo toàn đối xứng CPT trong các kênh biến mất của muon neutrino và phản neutrino.
  3. Đánh giá độ nhạy kết hợp (joint fit/synergy) giữa ba thí nghiệm T2K-II, phần mở rộng của NOvA (NOvA-II) và Jiangmen Underground Neutrino Observatory (JUNO) đến năm 2028 trong việc tìm kiếm vi phạm CP và thiết lập giới hạn vi phạm CPT.
  4. Xây dựng quy trình tái tạo vết ba chiều và phân tích phân bố động học của kênh tương tác neutrino tích điện tạo một pion ($CC1\pi$) tại hệ detector WAGASCI-BabyMIND.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng: Chùm muon (phản) neutrino tạo ra từ máy gia tốc proton $30\text{ GeV}$ tại Tổ hợp Nghiên cứu Máy gia tốc Proton Nhật Bản (J-PARC), được ghi nhận tại detector gần INGRID ($280\text{ m}$), hệ WAGASCI-BabyMIND ($1{,}5^\circ$ lệch trục) và detector xa Super-Kamiokande ($295\text{ km}$, $2{,}5^\circ$ lệch trục); dữ liệu từ thí nghiệm chùm hạt NOvA ($810\text{ km}$, $14\text{ mrad}$ lệch trục tại Hoa Kỳ) và thí nghiệm neutrino lò phản ứng JUNO ($52{,}5\text{ km}$ tại Trung Quốc).
  • Phạm vi dữ liệu: Dữ liệu thực nghiệm T2K thu thập từ Run 1 đến Run 10/11, kết hợp các ràng buộc từ các thí nghiệm lò phản ứng đường cơ sở ngắn (Daya Bay, RENO, Double Chooz), cùng các tập dữ liệu mô phỏng cho kế hoạch vận hành đến năm 2027–2028.

Tổng quan tài liệu và vị trí của luận án

Lịch sử nghiên cứu neutrino bắt đầu từ giả thuyết hạt trung hòa nhẹ của Wolfgang Pauli (1930) và tên gọi neutrino do Enrico Fermi đặt (1932). Clyde Cowan và Frederick Reines ghi nhận thực nghiệm đầu tiên về electron antineutrino vào năm 1956 qua phân rã beta ngược. Sự tồn tại của muon neutrino được Leon Lederman, Melvin Schwartz và Jack Steinberger xác nhận năm 1962 tại Alternating Gradient Synchrotron (AGS). Năm 2000, thí nghiệm DONUT tại Fermilab phát hiện trực tiếp tau neutrino. Hiện tượng dao động neutrino, khởi xướng từ ý tưởng chuyển pha neutrino - phản neutrino của Bruno Pontecorvo (1957) và mở rộng cho dao động mùi bởi Ziro Maki, Masami Nakagawa và Shoichi Sakata (1962), đã được xác nhận thực nghiệm tại Super-Kamiokande (1998) và SNO (2001), chứng minh neutrino có khối lượng khác không.

Về mặt đối xứng:

  • Vi phạm đối xứng chẵn lẻ (P) trong tương tác yếu được Tsung Dao Lee và Chen Ning Yang dự báo năm 1956, Chien Shiung Wu xác nhận thực nghiệm năm 1957.
  • Vi phạm đối xứng CP được Christenson, Cronin, Fitch và Turlay phát hiện trong hệ meson $K^0$ năm 1964, sau đó trong hệ meson $B$ tại BaBar và Belle (2001).
  • Đối xứng CPT là kết quả nền tảng của lý thuyết trường lượng tử dựa trên ba tiên đề: tính bất biến Lorentz, tính định xứ (locality) và tính Hermit của toán tử vi phân Hamiltonian (Schwinger 1951, Lüders 1954, Pauli 1955, Bell 1955).
  • Các thử nghiệm CPT trước đây tập trung vào hệ hạt - phản hạt, trong đó giới hạn đo đạc độ chênh lệch bình phương khối lượng của hệ kaon trung hòa là $|m^2(K^0) - m^2(\bar{K}^0)| < 0{,}25\text{ eV}^2$. Trong khi đó, các phép đo vi sai bình phương khối lượng neutrino $\Delta m^2_{21} \approx 7{,}45 \times 10^{-5}\text{ eV}^2$ và $|\Delta m^2_{31}| \approx 2{,}5 \times 10^{-3}\text{ eV}^2$ mang lại khả năng kiểm chứng CPT với độ nhạy cao hơn nhiều bậc độ lớn.

Vị trí của luận án: Luận án kết nối trực tiếp giữa công tác vận hành, mô phỏng detector gần tại J-PARC với việc phân tích dữ liệu dao động tầm xa của T2K và thiết lập mô hình tính toán độ nhạy liên hợp (T2K-II + NOvA-II + JUNO), giải quyết sự phụ thuộc chéo của các tham số dao động và triệt tiêu các nghiệm suy biến giả mạo do hiệu ứng vật chất gây ra.


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

Cơ sở lý thuyết

  1. Mô hình Chuẩn và tương tác yếu $V-A$: Tương tác dòng tích điện (CC) và dòng trung hòa (NC) có cấu trúc dạng véctơ trừ trục véctơ ($V-A$). Dưới phép chiếu chiral $P_L = \frac{1}{2}(1-\gamma^5)$ và $P_R = \frac{1}{2}(1+\gamma^5)$, chỉ các thành phần neutrino chiral trái ($\nu_L$) và phản neutrino chiral phải ($\nu_R^c$) tham gia vào tương tác dòng tích điện.
  2. Cơ chế sinh khối lượng Seesaw loại I: Bổ sung số hạng khối lượng Majorana cho neutrino phải ($\nu_R$) vào Lagrangian: $$\mathcal{L} = -\frac{1}{2}m_D(\bar{\nu}_L\nu_R + \bar{\nu}_R^c\nu_L^c) - \frac{1}{2}M\bar{\nu}R^c\nu_R + h.c.$$ Ma trận khối lượng cho ra hai trị riêng khối lượng: trạng thái neutrino nhẹ $|m\nu| \approx m_D^2/M$ và trạng thái neutrino nặng $m_N \approx M$. Với $m_D \sim 1\text{ GeV}$ và $M \sim 10^{10}\text{ GeV}$, khối lượng neutrino nhẹ thu được cỡ $0{,}1\text{ eV}$.
  3. Hình thức luận ma trận trộn PMNS trong chân không: Ba trạng thái mùi $(\nu_e, \nu_\mu, \nu_\tau)$ liên hệ với ba trạng thái khối lượng $(\nu_1, \nu_2, \nu_3)$ qua ma trận unita PMNS ($U_{PMNS}$), được tham số hóa bởi 3 góc trộn $(\theta_{12}, \theta_{13}, \theta_{23})$ và 1 pha Dirac $\delta_{CP}$: $$U = \begin{pmatrix} c_{12}c_{13} & s_{12}c_{13} & s_{13}e^{-i\delta_{CP}} \ -s_{12}c_{23}-c_{12}s_{23}s_{13}e^{i\delta_{CP}} & c_{12}c_{23}-s_{12}s_{23}s_{13}e^{i\delta_{CP}} & s_{23}c_{13} \ s_{12}s_{23}-c_{12}c_{23}s_{13}e^{i\delta_{CP}} & -c_{12}s_{23}-s_{12}c_{23}s_{13}e^{i\delta_{CP}} & c_{23}c_{13} \end{pmatrix}$$
  4. Hiệu ứng vật chất (MSW) và chuỗi nhiễu loạn Dyson: Khi truyền qua vật chất, tương tác tán xạ đàn hồi kết hợp của $\nu_e$ với electron tạo thế hiệu dụng $V_{CC} = \sqrt{2}G_F N_e$, còn tán xạ qua dòng trung hòa tạo thế $V_{NC} = -\frac{1}{2}\sqrt{2}G_F N_n$. Thế $V_{CC}$ làm biến đổi Hamiltonian và gây ra bất đối xứng CP biểu kiến (extrinsic/fake CP asymmetry) ngay cả khi $\delta_{CP} = 0$ hoặc $\pi$. Bất đối xứng CP tương đối gần cực đại dao động được tính bởi: $$A_{CP} \approx -\frac{\pi \sin 2\theta_{12}}{\tan \theta_{23} \sin 2\theta_{13}}\frac{\Delta m^2_{21}}{|\Delta m^2_{31}|}\sin \delta_{CP} \pm \frac{L}{2800\text{ km}}$$
  5. Kiểm chứng CPT: Theo định lý CPT, nếu bất biến CPT được bảo toàn thì khối lượng và góc trộn của hạt và phản hạt phải đồng nhất ($\Delta m^2_{31} = \Delta \bar{m}^2_{31}$ và $\theta_{23} = \bar{\theta}_{23}$). Bất kỳ sai lệch thực chất nào giữa hai bộ tham số này là dấu hiệu của vi phạm CPT.

Phương pháp nghiên cứu và công cụ phân tích

  • Mô phỏng chùm hạt và tương tác:
    • FLUKA 2011: Mô phỏng tương tác ban đầu của chùm proton $30\text{ GeV}$ với bia graphite, kết hợp chuẩn hóa dữ liệu thực nghiệm hadronic từ NA61/SHINE.
    • JNUBEAM (phiên bản 13a v1.1): Mô phỏng sự hội tụ qua hệ thấu kính từ trường horn và phân rã của các hạt mẹ ($\pi^\pm, K^\pm, K_L^0, \mu^\pm$) dọc ống phân rã.
    • NEUT (phiên bản 5.3.2 và 5.4.0): Mô phỏng tương tác neutrino với các bia hạt nhân ($H_2O, CH, Fe$) theo các mô hình tán xạ gần đàn hồi (CCQE - Llewellyn-Smith, Nieves 1p1h), kích thích cộng hưởng tạo meson đơn ($CC1\pi$ - Rein-Sehgal, Graczyk-Sobczyk có bổ sung hiệu chỉnh PCAC), tán xạ kết hợp tạo pion ($CC\text{coh}1\pi$ - Berger-Sehgal) và tán xạ không đàn hồi sâu (CCDIS/NCDIS).
    • GEANT4 (phiên bản 9.00, danh mục vật lý QGSP_BERT): Mô phỏng phản ứng của detector và tổn hao năng lượng qua vật chất.
  • Tái tạo biến cố: Thuật toán Cellular Automaton dùng để nối cụm tín hiệu (clusters) từ sợi chuyển bước sóng và cảm biến MPPC, tái tạo vết hai chiều (XZ, YZ) và phối hợp không gian ba chiều (3D matching).
  • Phân tích thống kê và độ nhạy dao động:
    • Bộ công cụ GLoBES (General Long Baseline Experiment Simulator) kết hợp ngôn ngữ định nghĩa thí nghiệm AEDL (Abstract Experiment Definition Language).
    • Xây dựng hàm $\Delta \chi^2$ kết hợp bao gồm các tham số hệ thống và ma trận tương quan để đánh giá khoảng tin cậy $1\sigma, 2\sigma, 3\sigma$ và $5\sigma$.

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

Chương 1: Neutrino oscillation phenomenon and experiments

Chương 1 thiết lập khung lý thuyết toán học về dao động neutrino trong chân không và môi trường vật chất, đồng thời trình bày cấu hình kỹ thuật của ba thí nghiệm chính:

  • T2K và T2K-II: Sử dụng chùm neutrino muon từ J-PARC truyền đến Super-Kamiokande ở khoảng cách $L = 295\text{ km}$, góc lệch trục $2{,}5^\circ$ để thu hẹp phổ năng lượng tại đỉnh $0{,}6\text{ GeV}$. Đề xuất T2K-II mở rộng thời gian chạy đến 2027 với tổng độ phơi chùm mục tiêu $20 \times 10^{21}\text{ POT}$.
  • NOvA và NOvA-II: Thí nghiệm đường cơ sở dài $L = 810\text{ km}$ tại Hoa Kỳ với góc lệch trục $14\text{ mrad}$. Do đường cơ sở dài, hiệu ứng vật chất tại NOvA đạt $\sim 28{,}9%$, cho phép phân giải trật tự khối lượng neutrino với độ nhạy cao.
  • JUNO: Thí nghiệm neutrino lò phản ứng đường cơ sở trung bình ($52{,}5\text{ km}$) tại Quảng Đông, Trung Quốc, sử dụng khối chất lỏng nhấp nháy $20\text{ kton}$ thu nhận electron antineutrino từ hai nhà máy điện hạt nhân Dương Giang ($6 \times 2{,}9\text{ GW}$) và Đài Sơn ($4 \times 4{,}59\text{ GW}$), tổng công suất nhiệt $36\text{ GW}$. Nhờ độ phân giải năng lượng cao, JUNO có khả năng đo đạc các tham số năng lượng mặt trời ($\Delta m^2_{21}, \sin^2\theta_{12}$) và $|\Delta m^2_{31}|$ với độ chính xác dưới $1%$.

Chương 2: Measurements at INGRID - the T2K on-axis near detector

Chương 2 trình bày kết quả đo đạc thực nghiệm và mô phỏng Monte Carlo tại detector gần đặt thẳng trục INGRID của T2K:

  • Cấu trúc detector INGRID: Gồm 14 module giống hệt nhau xếp thành hình chữ thập $10\text{ m} \times 10\text{ m}$ (7 module ngang, 7 module đứng). Mỗi module gồm 9 tấm sắt ($124 \times 124 \times 6{,}5\text{ cm}$) xen kẽ 11 lớp đo nhấp nháy (mỗi lớp 24 thanh), được bảo vệ chống tín hiệu ngoài bằng 6 mặt phẳng VETO (22 thanh nhấp nháy/mặt phẳng). Tín hiệu ánh sáng được thu gom qua sợi quang chuyển bước sóng và đọc bởi cảm biến mảng quang điện tử đa điểm (MPPC).
  • Mô phỏng thông lượng neutrino: Khảo sát 1100 tệp MC ($10^{21}\text{ POT}$/tệp) ở các chế độ dòng horn $0\text{ kA}$, $\pm 250\text{ kA}$ và $\pm 320\text{ kA}$. Kết quả chỉ ra khi áp dụng dòng horn, thông lượng neutrino tại đỉnh $1\text{ GeV}$ tăng từ 13–14 lần (với $250\text{ kA}$) và 14–15 lần (với $320\text{ kA}$) so với khi không có dòng horn ($0\text{ kA}$), giúp giảm thiểu phông tạp chất. Tăng dòng từ $250\text{ kA}$ lên $320\text{ kA}$ giúp tăng $10%$ thông lượng chùm tại đỉnh năng lượng cho cả chế độ neutrino (FHC) và phản neutrino (RHC).
  • Đo đạc tỷ lệ biến cố và biên dạng chùm (Run 10): Phân tích số liệu thực nghiệm Run 10 qua quy trình cắt 8 bước và thuật toán Cellular Automaton. Tỷ lệ biến cố ghi nhận hàng ngày và các thông số tâm chùm, độ rộng chùm theo phương ngang/phương đứng hoàn toàn ổn định và phù hợp với mô phỏng MC.
Cấu hình dòng Horn Tăng thông lượng so với 0 kA (tại đỉnh 1 GeV) Mức tăng khi nâng từ 250 kA lên 320 kA
Chế độ Neutrino (+250 kA / +320 kA) 13–14 lần (+250 kA) / 14–15 lần (+320 kA) $+10%$
Chế độ Phản Neutrino (-250 kA / -320 kA) 13–14 lần (-250 kA) / 14–15 lần (-320 kA) $+10%$

Chương 3: Testing CP and CPT invariances with neutrino oscillation measurements in T2K experiment

Chương 3 trình bày các kết quả phân tích số liệu thực nghiệm của T2K và tính toán độ nhạy liên hợp:

  • Kiểm chứng vi phạm CP:
    • Dữ liệu T2K thu thập từ Run 1 đến Run 9 với tổng độ phơi $3{,}13 \times 10^{21}\text{ POT}$ loại trừ các giá trị bảo toàn CP ($\delta_{CP} = 0, \pi$) ở mức độ tin cậy lớn hơn $95%$ ($>2\sigma$).
    • Giá trị ước lượng tốt nhất (best-fit) kèm sai số $\pm 1\sigma$ của pha Dirac $\delta_{CP}$ chỉ dùng số liệu T2K:
      • Trật tự khối lượng thuận (Normal Ordering - NO): $\delta_{CP} = -2{,}14^{+0{,}69}_{-0{,}69}$ radian.
      • Trật tự khối lượng nghịch (Inverted Ordering - IO): $\delta_{CP} = -1{,}26^{+0{,}61}_{-0{,}69}$ radian.
    • Khi kết hợp với các ràng buộc về góc $\theta_{13}$ từ các thí nghiệm lò phản ứng đường cơ sở ngắn:
      • Trật tự thuận (NO): $\delta_{CP} = -1{,}89^{+0{,}48}_{-0{,}58}$ radian.
      • Trật tự nghịch (IO): $\delta_{CP} = -1{,}41^{+0{,}55}_{-0{,}55}$ radian.
    • Phân tích kết hợp T2K-II + NOvA-II + JUNO: Đến năm 2028, phân tích đồng thời 3 thí nghiệm có khả năng khám phá vi phạm đối xứng CP ở mức ý nghĩa thống kê xấp xỉ $5\sigma$ C.L. nếu giá trị thực của $\delta_{CP}$ gần $-\pi/2$.
Tập số liệu phân tích Best-fit $\delta_{CP}$ (Trật tự thuận - NO) Best-fit $\delta_{CP}$ (Trật tự nghịch - IO) Mức loại trừ bảo toàn CP ($\delta_{CP}=0,\pi$)
Chỉ T2K ($3{,}13 \times 10^{21}\text{ POT}$) $-2{,}14^{+0{,}69}_{-0{,}69}$ $-1{,}26^{+0{,}61}_{-0{,}69}$ $> 95%$ C.L.
T2K + Ràng buộc lò phản ứng $-1{,}89^{+0{,}48}_{-0{,}58}$ $-1{,}41^{+0{,}55}_{-0{,}55}$ $> 95%$ C.L. (thu hẹp khoảng tin cậy)
T2K-II + NOvA-II + JUNO (Dự báo 2028) Nhạy cực đại tại $\sim -\pi/2$ Nhạy cực đại tại $\sim -\pi/2$ $\sim 5\sigma$ C.L.
  • Kiểm chứng vi phạm CPT:
    • Phân tích dữ liệu T2K với $3{,}13 \times 10^{21}\text{ POT}$ trong các kênh biến mất của $\nu_\mu$ và $\bar{\nu}_\mu$ cho thấy kết quả hoàn toàn phù hợp với giả thuyết bảo toàn đối xứng CPT.
    • Phân tích liên hợp T2K-II, NOvA-II và JUNO: Nếu các giá trị đo đạc best-fit hiện tại của T2K và NOvA đối với hiệu bình phương khối lượng $(\Delta m^2_{31}, \Delta \bar{m}^2_{31})$ và góc trộn $(\theta_{23}, \bar{\theta}_{23})$ là giá trị vật lý thực tế, sự kết hợp 3 thí nghiệm có khả năng loại trừ bảo toàn CPT ở mức ý nghĩa $1{,}6\sigma$ đến $2{,}2\sigma$.
    • Giới hạn trên đối với độ bất đối xứng CPT $|\delta(\Delta m^2_{31})| = |\Delta m^2_{31} - \Delta \bar{m}^2_{31}|$ được cải thiện đạt mức $5{,}0 \times 10^{-5}\text{ eV}^2$ ở mức tin cậy $3\sigma$ C.L. Đây là giới hạn chặt chẽ nhất thế giới trước khi các thí nghiệm DUNE hoặc Hyper-Kamiokande đi vào hoạt động (tốt hơn DUNE giai đoạn đầu và tốt hơn 1 bậc độ lớn so với các thí nghiệm riêng lẻ hiện hành).
    • Độ nhạy kiểm chứng CPT đối với tham số góc trộn $|\delta(\sin^2\theta_{23})|$ phụ thuộc lớn vào việc hai góc $\theta_{23}$ và $\bar{\theta}_{23}$ nằm cùng hay khác góc phần tám (octant). Khi bổ sung các mẫu biến cố xuất hiện electron neutrino ($\nu_e, \bar{\nu}_e$), độ nhạy thống kê loại trừ nghiệm bảo toàn CPT biểu kiến tăng lên rõ rệt.

Phụ lục A: Neutrino cross section measurements at WAGASCI BabyMIND

Phần phụ lục trình bày nghiên cứu thực nghiệm đo tiết diện tương tác neutrino trên bia nước ($H_2O$) và bia nhựa hydrocarbon ($CH$) nhằm giảm thiểu sai số hệ thống do sự khác biệt vật liệu giữa detector gần ND280 ($CH$) và detector xa Super-K ($H_2O$):

  • Cấu hình WAGASCI-BabyMIND đặt tại tầng B2 (góc lệch $1{,}5^\circ$):
    • Gồm 3 module trung tâm: Upstream WAGASCI (UWG), Proton Module (PM) và Downstream WAGASCI (DWG). Module WAGASCI chứa 1280 thanh nhấp nháy dạng lưới 3D, vận hành ở chế độ đổ đầy nước ($188\text{ kg}$ nước, chiếm $80%$ khối lượng vùng hoạt lực) hoặc không nước ($47\text{ kg}$ nhựa scintillator). Proton Module gồm 36 lớp theo dõi ($302\text{ kg}$ hydrocarbon).
    • Hệ thống đo quãng chạy muon: 2 module Wall-MRD đặt hai bên (11 tấm thép, 10 lớp nhấp nháy, nặng $8{,}6\text{ tấn}$) và detector từ hóa Baby MIND phía sau (33 module nam châm với tấm sắt dày $30\text{ mm}$, nặng 2 tấn/module, từ trường tối thiểu $1{,}5\text{ T}$).
  • Kết quả phân tích $CC1\pi$: Với $984$ tệp MC chuẩn hóa tương đương $5 \times 10^{20}\text{ POT}$ (khớp với độ phơi tích lũy Run 10 và 11 là $4{,}55 \times 10^{20}\text{ POT}$), các phân bố vết hạt ba chiều theo góc tán xạ và động lượng muon được tái tạo thành công, thể hiện sự phù hợp giữa mô phỏng và số liệu thực nghiệm.

Kết quả và những đóng góp mới

  1. Về mặt kỹ thuật và mô phỏng thực nghiệm:
    • Hoàn thành bộ tính toán mô phỏng thông lượng neutrino chi tiết cho detector INGRID dưới cấu hình nâng cấp dòng horn lên $320\text{ kA}$, chứng minh sự gia tăng $10%$ thông lượng chùm neutrino có ích.
    • Xác nhận sự ổn định của chùm hạt neutrino T2K trong suốt quá trình lấy dữ liệu Run 10 thông qua việc tái tạo biên dạng và tỷ lệ biến cố tại INGRID.
    • Xây dựng thành công thuật toán tái tạo vết 3D và phân tích động học kênh $CC1\pi$ cho hệ detector WAGASCI-BabyMIND, làm cơ sở để đo tỷ số tiết diện tương tác giữa $H_2O$ và $CH$ với độ chính xác mục tiêu $3%$.
  2. Về mặt kiểm chứng đối xứng CP:
    • Phân tích số liệu T2K tích lũy $3{,}13 \times 10^{21}\text{ POT}$ loại trừ miền giá trị bảo toàn CP ($\delta_{CP} = 0, \pi$) ở mức tin cậy trên $95%$.
    • Chứng minh về mặt định lượng rằng sự kết hợp giữa T2K-II, NOvA-II và JUNO vào năm 2028 có thể nâng độ nhạy phát hiện vi phạm CP lên mức tiêu chuẩn $5\sigma$ C.L., giải quyết triệt để sự suy biến trật tự khối lượng nhờ dữ liệu của JUNO và NOvA.
  3. Về mặt kiểm chứng đối xứng CPT:
    • Xác nhận dữ liệu thực nghiệm hiện tại của T2K hoàn toàn nhất quán với định lý CPT.
    • Xác lập rằng phân tích liên hợp ba thí nghiệm T2K-II, NOvA-II và JUNO sẽ thiết lập giới hạn vi phạm CPT đối với tham số hiệu bình phương khối lượng $|\Delta m^2_{31} - \Delta \bar{m}^2_{31}| \le 5{,}0 \times 10^{-5}\text{ eV}^2$ tại $3\sigma$ C.L., trở thành kết quả đo có độ chính xác cao nhất thế giới tại thời điểm năm 2028.

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

  • Hạn chế trong văn bản:
    • Độ nhạy kiểm chứng CPT đối với tham số góc trộn $\sin^2\theta_{23}$ còn bị giới hạn bởi sự thoái biến góc phần tám khi chỉ sử dụng các kênh biến mất muon neutrino/antineutrino.
    • Việc bổ sung các kênh xuất hiện ($\nu_e, \bar{\nu}e$) làm tăng độ nhạy phân tách CPT nhưng lại chịu sự phụ thuộc chéo vào pha $\delta{CP}$ và trật tự khối lượng chưa được xác định tuyệt đối.
  • Hướng nghiên cứu tiếp theo:
    • Thực hiện khớp dữ liệu thực (real data fit) đồng thời giữa các tập số liệu đo đạc thực tế của T2K, NOvA và MINOS/MINOS+.
    • Mở rộng phân tích độ nhạy kiểm chứng bất biến CPT cho các thí nghiệm thế hệ kế tiếp bao gồm Hyper-Kamiokande và Deep Underground Neutrino Experiment (DUNE).

Giá trị tham khảo

  • Đối tượng hưởng lợi:
    • Các nghiên cứu sinh, nhà nghiên cứu lý thuyết và thực nghiệm trong lĩnh vực vật lý hạt cơ bản, vật lý năng lượng cao và vật lý neutrino.
    • Các kỹ sư và chuyên gia phân tích dữ liệu tại các thí nghiệm máy gia tốc hạt và lò phản ứng hạt nhân.
  • Các phần nội dung đáng tham khảo:
    • Phương pháp xây dựng mô phỏng chuỗi chùm hạt kết hợp FLUKA, JNUBEAM, NEUT, GEANT4 trong Chương 2 và Phụ lục A.
    • Khung toán học và kỹ thuật triển khai GLoBES/AEDL để phân tích độ nhạy liên hợp đa thí nghiệm và khử suy biến tham số trong Chương 3.
    • Hình thức luận tường minh về kiểm chứng CPT thông qua các kênh dao động neutrino và phản neutrino.

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

1. Giá trị đo đạc pha vi phạm CP ($\delta_{CP}$) của thí nghiệm T2K đạt được là bao nhiêu khi có kết hợp với ràng buộc từ các thí nghiệm lò phản ứng?

Khi kết hợp dữ liệu T2K ($3{,}13 \times 10^{21}\text{ POT}$) với ràng buộc góc $\theta_{13}$ từ các thí nghiệm lò phản ứng đường cơ sở ngắn (Daya Bay, RENO, Double Chooz), giá trị best-fit kèm khoảng tin cậy $\pm 1\sigma$ của $\delta_{CP}$ đạt $-1{,}89^{+0{,}48}{-0{,}58}\text{ rad}$ đối với trật tự khối lượng thuận (NO) và $-1{,}41^{+0{,}55}{-0{,}55}\text{ rad}$ đối với trật tự khối lượng nghịch (IO).

2. Phân tích liên hợp T2K-II, NOvA-II và JUNO có thể đạt độ nhạy tìm kiếm vi phạm CP ở mức ý nghĩa thống kê nào vào năm 2028?

Vào năm 2028, sự kết hợp đồng thời số liệu từ T2K-II, NOvA-II và JUNO có khả năng loại trừ giả thuyết bảo toàn CP ($\delta_{CP} = 0, \pi$) ở mức ý nghĩa thống kê xấp xỉ $5\sigma$ C.L. nếu giá trị thực của pha $\delta_{CP}$ nằm gần $-\pi/2$.

3. Việc nâng dòng điện horn từ 250 kA lên 320 kA tại T2K mang lại hiệu quả gì đối với thông lượng neutrino tại INGRID?

Khi dòng điện horn tăng từ $250\text{ kA}$ lên $320\text{ kA}$, thông lượng neutrino ở vùng đỉnh năng lượng (khoảng $1\text{ GeV}$) tăng thêm $10%$ cho cả chế độ neutrino (FHC) và chế độ phản neutrino (RHC), đồng thời tăng cường độ hội tụ các hạt mẹ và giảm tỷ lệ phông tạp chất.

4. Giới hạn vi phạm đối xứng CPT đối với hiệu bình phương khối lượng $|\Delta m^2_{31} - \Delta \bar{m}^2_{31}|$ mà liên hợp ba thí nghiệm đạt được là bao nhiêu?

Phân tích liên hợp T2K-II, NOvA-II và JUNO có thể thiết lập giới hạn trên cho độ bất đối xứng $|\delta(\Delta m^2_{31})| = |\Delta m^2_{31} - \Delta \bar{m}^2_{31}|$ ở mức $5{,}0 \times 10^{-5}\text{ eV}^2$ tại $3\sigma$ C.L., mức giới hạn chặt chẽ nhất thế giới trước khi các thí nghiệm thế hệ mới vận hành.

5. Mục tiêu chính của hệ detector WAGASCI-BabyMIND trong tổ hợp T2K là gì?

WAGASCI-BabyMIND được thiết kế nhằm đo đạc chính xác tỷ số tiết diện tương tác dòng tích điện giữa bia nước ($H_2O$) và bia nhựa nhấp nháy hydrocarbon ($CH$) với độ chính xác mục tiêu $3%$, qua đó triệt tiêu nguồn sai số hệ thống lớn nhất bắt nguồn từ sự khác biệt vật liệu giữa detector gần ND280 ($CH$) và detector xa Super-Kamiokande ($H_2O$).


Kết luận

Luận án tiến sĩ của tác giả Trần Văn Ngọc đã thực hiện nghiên cứu thực nghiệm và tính toán mô phỏng về hiện tượng dao động neutrino nhằm kiểm chứng hai đối xứng rời rạc cơ bản CP và CPT. Tác giả đã đóng góp trực tiếp vào công tác giám sát chùm hạt tại detector gần INGRID, xử lý số liệu kênh $CC1\pi$ tại WAGASCI-BabyMIND, và phân tích tập dữ liệu $3{,}13 \times 10^{21}\text{ POT}$ của T2K để xác định pha Dirac $\delta_{CP}$ với mức loại trừ bảo toàn CP trên $95%$ C.L. Đồng thời, công trình chứng minh định lượng rằng phân tích liên hợp T2K-II, NOvA-II và JUNO đến năm 2028 có thể đạt độ nhạy khám phá vi phạm CP ở mức $\sim 5\sigma$ C.L. và thiết lập kỷ lục thế giới về giới hạn kiểm chứng bất biến CPT trong khu vực lepton ($|\Delta m^2_{31} - \Delta \bar{m}^2_{31}| \le 5{,}0 \times 10^{-5}\text{ eV}^2$).