Giới thiệu dự án

Nghiên cứu tương tác giữa trường laser cường độ cao ($10^{14} - 10^{15}\text{ W/cm}^2$) với nguyên tử và phân tử là một trong những trụ cột của vật lý quang học phi tuyến và khoa học atto-giây (Attosecond Science) – lĩnh vực đã được vinh danh tại giải Nobel Vật lý 2023. Khi đặt hệ nguyên tử trong điện trường ngoài cực mạnh, rào thế Coulomb bị biến dạng sâu sắc, kích hoạt hàng loạt hiệu ứng phi tuyến phi nhiễu loạn như phát xạ sóng điều hòa bậc cao (High-order Harmonic Generation - HHG), ion hóa trên ngưỡng (Above Threshold Ionization - ATI) và ion hóa hai lần (Double Ionization - DI).

   [Laser Pulse: E(t)]

Trong quá trình ion hóa hai lần, cơ chế ion hóa hai lần không liên tiếp (Nonsequential Double Ionization - NSDI) đóng vai trò chìa khóa để giải mã bản chất tương quan tương tác giữa electron - electron ($e-e$) và electron - ion ($e-i$). Mô hình 3 bước kinh điển (Corkum, 1993) chỉ ra rằng: electron thứ nhất sau khi bứt ra sẽ được gia tốc bởi trường laser, khi điện trường đổi chiều sẽ quay trở lại va chạm với ion mẹ. Tuy nhiên, một bài toán vật lý then chốt vẫn tồn tại nhiều tranh luận: Liệu electron thứ nhất có ion hóa electron thứ hai ngay trong lần va chạm đầu tiên, hay phải trải qua quá trình tái va chạm nhiều lần (Multiple Recollisions) thông qua trạng thái kích thích kép?

Đề tài khóa luận "Ảnh hưởng của sự tái va chạm nhiều lần lên quá trình ion hóa kép không liên tục" (thực hiện bởi Trần Ngọc Liên Hương, hướng dẫn bởi TS. Phạm Nguyễn Thành Vinh tại Trường ĐH Sư phạm TP.HCM) tập trung giải quyết triệt để bài toán động lực học này trên nguyên tử Argon ($\text{Ar}$).

Mục tiêu nghiên cứu cụ thể

  1. Xây dựng chương trình mô phỏng số 3 chiều cổ điển (3D Classical Ensemble Model) tích hợp thuật toán Runge-Kutta bậc 4 (RK4) để giải hệ phương trình vi phân chuyển động Newton cho $2 \times 10^6$ nguyên tử Argon.
  2. Khảo sát chi tiết phổ phân bố động lượng tương quan hai electron (Correlated Two-Electron Momentum Distribution - CTEMD) dọc theo trục phân cực laser.
  3. Định lượng tỷ lệ đóng góp và phân tích vi mô quỹ đạo của các sự kiện tái va chạm 1 lần (Single Recollision) so với tái va chạm nhiều lần (Multiple Recollisions).
  4. Khảo sát sự phụ thuộc phi tuyến của cơ chế tái va chạm vào cường độ laser ($I = 0.8 \times 10^{14}\text{ W/cm}^2$, $1.5 \times 10^{14}\text{ W/cm}^2$, $2.5 \times 10^{14}\text{ W/cm}^2$) và bước sóng trường bức xạ ($\lambda = 780\text{ nm}$ và $1200\text{ nm}$).

Phạm vi và giới hạn

  • Đối tượng: $2.000.000$ nguyên tử Argon trong hệ tọa độ Descartes 3D.
  • Điện trường laser: Phân cực thẳng dọc trục $x$, dạng hình bao xung thang gồm 10 chu kỳ quang học ($N_1=2$ chu kỳ bật, $N_2=6$ chu kỳ đỉnh, $N_3=2$ chu kỳ tắt).
  • Giới hạn mô hình: Tiếp cận hoàn toàn theo cơ chế cơ học cổ điển và ion hóa vượt rào (Over-the-barrier ionization), sử dụng thế Coulomb làm mềm (soft-core potential) để loại bỏ điểm kỳ dị tính toán.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí so sánh Phương trình Schrödinger phụ thuộc thời gian (TDSE) Gần đúng trường mạnh (SFA - Lewenstein) Mô hình tập hợp 3D cổ điển (Đề tài áp dụng)
Độ chính xác lượng tử Tuyệt đối (Exact solution) Trung bình (Bỏ qua tương tác Coulomb trạng thái liên tục) Cao ở vùng trường mạnh ($I > I_{th}$)
Khả năng truy vết quỹ đạo Kém (Chỉ cho thông tin hàm sóng trạng thái cuối) Kém (Tích phân pha lượng tử) Rất cao (Truy vết liên tục vị trí, vận tốc, năng lượng từng thời điểm)
Chi phí tính toán (RAM/CPU) Rất lớn (Lưới 6D đòi hỏi cụm siêu máy tính HPC) Thấp đến trung bình Tối ưu, song song hóa cao trên tập hợp triệu hạt
Phân tích va chạm nhiều lần Rất khó tách biệt cơ chế Không hỗ trợ đầy đủ tái va chạm bậc cao Tách biệt chính xác số lần quay về ($1^{st}, 2^{nd}, 3^{rd}, ...$)

Phân tích yêu cầu kỹ thuật (MoSCoW)

  • Must have: Hệ phương trình vi phân Newton tương tác Coulomb 3 thân ($e_1 - e_2 - \text{Ion}$); Thuật toán tích phân số RK4 bảo toàn năng lượng; Thuật toán phân loại chính xác sự kiện ion hóa kép (DI) khi $E_1 > 0$ và $E_2 > 0$.
  • Should have: Bộ lọc nhận diện tái va chạm khi khoảng cách giữa 2 electron $r_{12} < 2.0\text{ a.u.}$; Biểu đồ mật độ 2D cho phổ CTEMD $(p_{1x}, p_{2x})$.
  • Could have: Tự động thống kê phân bố thời gian trễ $\Delta t = t_{DI} - t_r$ và năng lượng hồi quy $E_r$.
  • Won't have: Hiệu ứng giao thoa pha lượng tử thuần túy (Quantum Interference fringes) của cơ học sóng.

Thiết kế hệ thống mô phỏng

Công thức toán học và công nghệ áp dụng

Hệ phương trình vi phân vector chi phối chuyển động của hai electron trong trường laser phân cực thẳng $\mathbf{E}(t) = E_x(t)\mathbf{e}_x$:

$$\frac{d^2\mathbf{r}i}{dt^2} = -\nabla{\mathbf{r}i} \left[ V{ne}(r_i) + V_{ee}(r_1, r_2) \right] - \mathbf{E}(t)$$

Trong đó, đơn vị nguyên tử ($\text{a.u.}: e = \hbar = m_e = 1$) được áp dụng. Thế tương tác hút hạt nhân - electron ($V_{ne}$) và tương tác đẩy giữa hai electron ($V_{ee}$) sử dụng tham số làm mềm $a = 1.0\text{ a.u.}$ và $b = 1.0\text{ a.u.}$ để chống tự ion hóa:

$$V_{ne}(r_i) = -\frac{2}{\sqrt{r_i^2 + a^2}}, \quad V_{ee}(r_1, r_2) = \frac{1}{\sqrt{(\mathbf{r}_1 - \mathbf{r}_2)^2 + b^2}}$$

Năng lượng tức thời của mỗi electron $E_i(t)$ để xác định điều kiện ion hóa:

$$E_i(t) = \frac{v_{xi}^2 + v_{yi}^2 + v_{zi}^2}{2} - \frac{2}{\sqrt{x_i^2 + y_i^2 + z_i^2 + a^2}} + \frac{1}{2\sqrt{(x_1-x_2)^2 + (y_1-y_2)^2 + (z_1-z_2)^2 + b^2}}$$

  • Technology Stack:
    • Ngôn ngữ lập trình tính toán: C++17 (tối ưu hóa chỉ thị SIMD AVX-512) kết hợp Python 3.9 (phân tích dữ liệu).
    • Thư viện: OpenMP 4.5 (song song hóa đa lõi), NumPy 1.21, SciPy 1.7, Matplotlib 3.5 (kết xuất đồ thị động lượng).
    • Môi trường tính toán: Linux Ubuntu 20.04 LTS, bộ nhớ RAM yêu cầu 16 GB, xử lý $2 \times 10^6$ hạt trong ~4.2 giờ CPU.

Implementation và kết quả

Chi tiết thuật toán và mã nguồn mô phỏng

Thuật toán Runge-Kutta bậc 4 (RK4) giải hệ phương trình vi phân chuyển động vi mô cho mỗi tọa độ không gian:

$$y_{n+1} = y_n + \frac{h}{6}(K_1 + 2K_2 + 2K_3 + K_4)$$

// Core Classical Ensemble RK4 Engine for NSDI Argon Simulation
#include <iostream>
#include <vector>
#include <cmath>

struct ElectronState {
    double x, y, z;
    double vx, vy, vz;
};

struct SystemState {
    ElectronState e1;
    ElectronState e2;
    double time;
};

// Compute acceleration considering Laser Field & Soft-core Coulomb potentials
void computeAcceleration(const SystemState& state, double laserField, 
                         double a2, double b2,
                         double& ax1, double& ay1, double& az1,
                         double& ax2, double& ay2, double& az2) {
    // Distance from nucleus
    double r1_sq = state.e1.x * state.e1.x + state.e1.y * state.e1.y + state.e1.z * state.e1.z;
    double r2_sq = state.e2.x * state.e2.x + state.e2.y * state.e2.y + state.e2.z * state.e2.z;
    double term_ne1 = 2.0 / std::pow(r1_sq + a2, 1.5);
    double term_ne2 = 2.0 / std::pow(r2_sq + a2, 1.5);

    // Distance between two electrons
    double dx = state.e1.x - state.e2.x;
    double dy = state.e1.y - state.e2.y;
    double dz = state.e1.z - state.e2.z;
    double r12_sq = dx * dx + dy * dy + dz * dz;
    double term_ee = 1.0 / std::pow(r12_sq + b2, 1.5);

    // Acceleration on Electron 1
    ax1 = -state.e1.x * term_ne1 + dx * term_ee - laserField;
    ay1 = -state.e1.y * term_ne1 + dy * term_ee;
    az1 = -state.e1.z * term_ne1 + dz * term_ee;

    // Acceleration on Electron 2
    ax2 = -state.e2.x * term_ne2 - dx * term_ee - laserField;
    ay2 = -state.e2.y * term_ne2 - dy * term_ee;
    az2 = -state.e2.z * term_ne2 - dz * term_ee;
}

Dữ liệu thực nghiệm và kết quả phân tích quỹ đạo

1. Vùng cường độ thấp ($I = 0.8 \times 10^{14}\text{ W/cm}^2, \lambda = 780\text{ nm}$)

  • Đặc trưng phân bố CTEMD: Phổ động lượng tập trung chủ yếu trên đường chéo phụ tại góc phần tư thứ hai ($p_{1x} < 0, p_{2x} > 0$) và thứ tư ($p_{1x} > 0, p_{2x} < 0$). Hai electron thoát ra với động lượng bằng nhau về độ lớn nhưng ngược chiều ($p_{1x} \approx -p_{2x}$).
  • Cơ chế chiếm ưu thế: Ion hóa thông qua trạng thái kích thích do tái va chạm (RESI - Recollision Excitation with Subsequent Ionization). Electron thứ nhất không đủ động năng ($< 3.17 U_p$) để bứt electron thứ hai ngay lập tức, hình thành trạng thái kích thích kép.
  • Tái va chạm 1 lần vs 2 lần:
    • Tái va chạm 1 lần chiếm 88% tổng số sự kiện DI, phân bố hoàn toàn ở góc phần tư II và IV.
    • Tái va chạm 2 lần chiếm 12% tổng số sự kiện DI, phân bố đều đối xứng ở cả 4 góc phần tư.
                  p2x (Động lượng e2)

2. Phân bố thời gian động học và độ trễ

  • Khoảng thời gian từ lúc ion hóa đơn đến khi tái va chạm ($t_r - t_{SI}$) ghi nhận 3 đỉnh cực đại rõ rệt tại $0.65 T_0$, $1.75 T_0$, và $2.65 T_0$ (tương ứng với lần quay về thứ nhất, thứ ba và thứ năm). Ở các lần quay về chẵn ($2^{nd}, 4^{th}$), động năng electron hồi quy triệt tiêu nên không thể gây kích thích ion mẹ.
  • Thời gian trễ giữa hai lần tái va chạm liên tiếp $\Delta t_r = t_{r2} - t_{r1}$ hội tụ chính xác tại $0.5 T_0$ (nửa chu kỳ quang học).
  • Thời gian trễ ion hóa $\Delta t_{delay} = t_{DI} - t_r > 0.3 T_0$ ở cường độ thấp, chứng minh electron phải chờ điện trường đạt đỉnh chu kỳ tiếp theo mới bứt khỏi rào thế.

3. Sự phụ thuộc vào cường độ trường laser ($I = 1.5 \times 10^{14}$ và $2.5 \times 10^{14}\text{ W/cm}^2$)

  • Khi tăng cường độ lên $2.5 \times 10^{14}\text{ W/cm}^2$, động lượng tương quan dịch chuyển đột ngột sang góc phần tư thứ nhất và thứ ba ($p_{1x} \cdot p_{2x} > 0$). Hai electron thoát ra cùng hướng (side-by-side emission).
  • Cơ chế va chạm trực tiếp $(e, 2e)$ áp đảo hoàn toàn. Thời gian trễ $t_{DI} - t_r$ co ngắn xuống $< 0.3 T_0$. Hiện tượng chia sẻ năng lượng bất đối xứng (AES - Asymmetric Energy Sharing) xuất hiện rõ nét. Đóng góp của tái va chạm nhiều lần giảm mạnh do electron mang năng lượng vượt trội ngay trong lần quay về đầu tiên.

4. Sự phụ thuộc vào bước sóng ($\lambda = 1200\text{ nm}$ vs $780\text{ nm}$)

  • Năng lượng truyền động Ponderomotive $U_p = I / (4\omega^2) \propto I \lambda^2$ tăng gấp $2.37$ lần khi chuyển từ $780\text{ nm}$ sang $1200\text{ nm}$ tại cùng cường độ $0.8 \times 10^{14}\text{ W/cm}^2$.
  • Tỷ lệ sự kiện tái va chạm nhiều lần giảm từ $12%$ xuống chỉ còn $5%$, đồng thời phổ CTEMD chuyển dịch từ phân bố góc phần tư II/IV sang phân bố dọc các trục tọa độ ở góc phần tư I/III.

Đổi mới và đóng góp

  1. Làm sáng tỏ vai trò vi mô của tái va chạm nhiều lần: Lần đầu tiên chỉ ra rằng cấu trúc phân bố động lượng tương quan đồng đều ở cả 4 góc phần tư tại cường độ thấp bắt nguồn từ các electron trải qua từ 2 lần tái va chạm trở lên ($\Delta t \approx 0.5 T_0$).
  2. Khẳng định tính chất lượng tử qua mô hình tập hợp cổ điển: Chứng minh mô hình tập hợp 3D cổ điển hoàn toàn có thể tái lập chính xác các cấu trúc phổ phân bố động lượng thực nghiệm phức tạp mà trước đây bắt buộc phải giải hệ TDSE 6D tốn kém.
  3. Tiết kiệm tài nguyên tính toán vượt bậc: Giảm thiểu hơn $98.5%$ thời gian tính toán so với việc giải phương trình sóng TDSE trên mạng lưới không gian 3 chiều, cho phép thực hiện khảo sát thống kê tin cậy trên quy mô $2.000.000$ quỹ đạo hạt.
  4. Đóng góp học thuật: Cung cấp cơ sở lý thuyết định lượng hỗ trợ các nhóm thực nghiệm đo đạc phổ COLTRIMS (Cold Target Recoil Ion Momentum Spectroscopy) trong việc hiệu chỉnh thông số xung laser atto-giây.

Ứng dụng thực tế và triển khai

  • Tối ưu hóa nguồn phát xung cực ngắn Atto-giây: Điều khiển thời điểm tái va chạm và kiểm soát quỹ đạo electron giúp cực đại hóa tần số cắt (Cut-off frequency) trong quá trình phát xạ sóng điều hòa bậc cao (HHG), tạo tiền đề cho các nguồn laser tia X mềm công suất cao.
  • Kỹ thuật chụp ảnh động học phân tử (Molecular Ultrafast Imaging): Định lượng chính xác tương quan va chạm $e-e$ cho phép tái tạo cấu trúc không gian và theo dõi sự dịch chuyển mật độ điện tích trong liên kết hóa học với độ phân giải thời gian cỡ sub-femtosecond.
  • Tích hợp vào phần mềm mô phỏng vật lý chuyên dụng: Thuật toán có thể đóng gói thành thư viện mô phỏng tương tác trường mạnh cho các viện nghiên cứu vật lý plasma và công nghệ bán dẫn nano.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật:
    • Mô hình cổ điển thuần túy chưa tính toán được hiệu ứng xuyên hầm lượng tử dưới rào thế (Tunneling ionization) ở giai đoạn khởi tạo ban đầu, mà sử dụng cơ chế vượt rào.
    • Chưa tính đến hiệu ứng giao thoa hàm sóng nhiều đường đi (Quantum trajectory interference) tạo nên các vân giao thoa vi mô trên phổ động lượng.
  • Hướng nâng cấp tương lai:
    • Kết hợp mô hình bán cổ điển tích hợp pha lượng tử (Semiclassical SCTS - Semiclassical Two-Step Model).
    • Mở rộng thuật toán cho các phân tử đa nguyên tử bất đối xứng ($CO_2, N_2, O_2$) và xung laser phân cực elip hoặc phân cực tròn hai màu (Two-color orthogonal laser fields).

Đối tượng hưởng lợi

+-------------------------------------------------------------------------------+
|                             ĐỐI TƯỢNG HƯỞNG LỢI                              |
+-------------------------------------------------------------------------------+
|  [Sinh viên / Học viên cao học]                                               |
|  - Nắm vững cơ chế vật lý tương tác laser trường mạnh (RESI vs (e,2e)).       |
|  - Tài liệu mẫu chuẩn mực về giải số phương trình vi phân và mô hình hóa 3D.  |
+-------------------------------------------------------------------------------+
|  [Kỹ sư phần mềm / Nhà tính toán khoa học]                                    |
|  - Tiếp cận kiến trúc mã nguồn C++/Python tích hợp thuật toán RK4 tối ưu.     |
|  - Kỹ thuật xử lý dữ liệu lớn (Big Data) trên 2 triệu hạt vật lý.             |
+-------------------------------------------------------------------------------+
|  [Nhà nghiên cứu / Phòng thí nghiệm Laser]                                    |
|  - Tiết kiệm 98% chi phí tính toán lý thuyết trước khi thực nghiệm COLTRIMS.  |
|  - Bộ thông số chuẩn về thời gian trễ và phân bố năng lượng electron.        |
+-------------------------------------------------------------------------------+

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

1. Yêu cầu phần cứng để triển khai mô phỏng tập hợp 2 triệu nguyên tử là gì?

Hệ thống yêu cầu tối thiểu CPU 8 nhân (hỗ trợ tập lệnh AVX2/AVX-512), 16 GB RAM DDR4 và ổ cứng SSD 50 GB để lưu trữ dữ liệu quỹ đạo chi tiết. Với OpenMP, thời gian hoàn thành một bộ thông số ($I, \lambda$) dao động từ 3 đến 5 giờ tính toán.

2. Vì sao hiện tượng tái va chạm 2 lần lại giảm mạnh khi tăng bước sóng từ 780 nm lên 1200 nm?

Năng lượng truyền động $U_p \propto I\lambda^2$. Khi bước sóng tăng lên $1200\text{ nm}$, động năng electron quay về ở chu kỳ đầu tiên tăng gấp 2.37 lần, vượt xa thế ion hóa của electron thứ hai. Do đó, quá trình ion hóa kép diễn ra ngay lập tức qua phản ứng $(e, 2e)$ mà không cần tái va chạm lần hai, khiến tỷ lệ giảm từ $12%$ xuống $5%$.

3. Tham số làm mềm Coulomb ($a = 1.0\text{ a.u.}$, $b = 1.0\text{ a.u.}$) có ảnh hưởng như thế nào đến tính chính xác của mô hình?

Tham số làm mềm triệt tiêu điểm kỳ dị khi khoảng cách $r \to 0$, ngăn chặn hiện tượng lực hút/đẩy tiến tới vô cùng gây sụp đổ tích phân số trong cơ học cổ điển. Giá trị $a = b = 1.0\text{ a.u.}$ được chuẩn hóa theo thực nghiệm để tái lập chính xác thế ion hóa trạng thái cơ bản của Argon ($-1.59\text{ a.u.}$).

4. Phổ CTEMD phân bố ở góc phần tư II & IV khác biệt như thế nào so với góc phần tư I & III?

  • Góc phần tư II & IV ($p_{1x} \cdot p_{2x} < 0$): Hai electron bay ra theo hai hướng đối nghịch nhau (Back-to-back emission), đặc trưng của cơ chế RESI qua trạng thái kích thích kép ở vùng cường độ thấp.
  • Góc phần tư I & III ($p_{1x} \cdot p_{2x} > 0$): Hai electron bay ra cùng một phía (Side-by-side emission), đặc trưng của cơ chế va chạm trực tiếp $(e, 2e)$ ở vùng cường độ cao.

5. Làm thế nào để xác định chính xác thời điểm xảy ra sự kiện tái va chạm trong thuật toán?

Trong quá trình tích phân từng bước thời gian $dt$, thuật toán liên tục kiểm tra khoảng cách tương đối giữa 2 electron. Một sự kiện được định nghĩa là tái va chạm khi khoảng cách $r_{12} = |\mathbf{r}_1 - \mathbf{r}_2| < 2.0\text{ a.u.}$ sau khi electron thứ nhất đã bị ion hóa ($E_1(t) > 0$).

Kết luận

Đề tài khóa luận tốt nghiệp của tác giả Trần Ngọc Liên Hương đã giải quyết xuất sắc bài toán động lực học tương tác phi tuyến trường mạnh, chứng minh tầm ảnh hưởng quyết định của sự tái va chạm nhiều lần lên cấu trúc phổ động lượng tương quan hai electron trong quá trình ion hóa kép không liên tiếp (NSDI) của nguyên tử Argon. Bằng việc kết hợp mô hình tập hợp 3D cổ điển và thuật toán tích phân Runge-Kutta bậc 4 trên tập dữ liệu lớn ($2 \times 10^6$ hạt), công trình đã cung cấp bức tranh vật lý toàn diện về sự chuyển dịch cơ chế từ RESI (kích thích trễ) sang $(e, 2e)$ (va chạm trực tiếp) dưới sự điều biến của cường độ và bước sóng laser. Đây là tài liệu khoa học mẫu mực, kết nối chặt chẽ giữa lý thuyết toán lý, kỹ thuật mô phỏng số hiệu năng cao và ứng dụng thực tiễn trong ngành vật lý quang học atto-giây hiện đại.