Khảo sát ion hóa kép Argon bằng thuật toán Runge-Kutta

Khóa luận vật lý: Nghiên cứu ion hóa kép nguyên tử Argon bằng laser phân cực thẳng, ứng dụng thuật toán Runge Kutta. Chi tiết kết quả khảo sát.

• Tác giả: Trần Ngọc Liên Hương
• Người hướng dẫn: TS. Phạm Nguyễn Thành Vinh
• Trường học: Trường Đại học Sư phạm Thành phố Hồ Chí Minh
• Chuyên ngành: Vật lý
• Đề tài: Áp Dụng Thuật Toán Runge - Kutta Để Khảo Sát Sự Ion Hóa Kép Không Liên Tục Của Nguyên Tử Argon Dưới Tác Dụng Của Laser Phân Cực Thẳng
• Loại tài liệu: khóa luận tốt nghiệp
• Năm xuất bản: 2016
• Địa điểm: Thành phố Hồ Chí Minh

Mục lục chi tiết

• MỤC LỤC
• 1. CHƯƠNG 1: Quá trình tương tác giữa laser va nguyên tử, phân tử
• 1.1. Quá trình ion hóa
• 1.2. Quá trình ion hóa kép
• 1.3. Quá trình ion hóa kép liên tục
• 1.4. Quá trình ion hóa kép không liên tục
• 2. CHƯƠNG 2: CƠ SỞ LY THUYẾT VA PHƯƠNG PHÁP TÍNH
• 2.1. Mô hình tập hợp ba chiêu cô điển
• 2.2. Phương pháp tích phân liên tiếp
• 2.2.1. Phương pháp Runge — Kutta bậc 2
• 2.2.2. Phương pháp Runge — Kutta bậc 4
• 3. CHƯƠNG 3: KET QUA NGHIÊN CUU
• 3.1. Kiểm chứng tinh chính xác của thuật toán
• 3.1.1. Bài toán dao động
• 3.1.2. Tính chính xác của thuật ton
• 3.2. Khảo sát quá trình ion hóa kép không liên tục của nguyên tử Argon bằng mô hính tập Đợp ba chu cổ hans
• 3.3. Khảo sát sự phụ thuộc vào độ dài xung của quátrình NSDI của nguyên tử IẢÌDH
• 3.4. Khảo sát sự phụ thuộc vào cường độ laser của NSDI của nguyên từ Argon
• KET LUẬN VÀ HƯỚNG PHAT TRIÊN
• TATU THAM KẢ
• Phụ lục 1: Phương pháp Runge — Kutta bậc 2
• Phụ lục 2: Phương pháp Runge — Kutta bậc 4
• Danh mục ký hiệu và chữ viết tắt
• Danh mục hình vẽ, đồ thị
• Danh mục hình vẽ
• Danh mục đồ thị
• MỞ ĐÀU

Tóm tắt

I. Ion hóa kép Argon là gì Toàn cảnh cơ chế vật lý laser

Ion hóa kép Argon là một hiện tượng phức tạp trong lĩnh vực vật lý nguyên tử, nơi một nguyên tử Argon mất đi hai electron dưới tác dụng của một trường laser cường độ cao. Quá trình này không chỉ là một bài toán cơ bản của cơ học lượng tử mà còn là chìa khóa để hiểu sâu hơn về tương tác laser-vật chất. Khi một xung laser cực ngắn tương tác với nguyên tử, các hiệu ứng phi tuyến bắt đầu xảy ra, dẫn đến việc giải phóng electron. Tùy thuộc vào cường độ của trường laser so với trường Coulomb của nguyên tử mà cơ chế ion hóa sẽ khác nhau. Ở cường độ yếu, hiện tượng ion hóa đa photon chiếm ưu thế, nơi electron hấp thụ tuần tự nhiều photon để đạt đủ năng lượng thoát ra. Khi cường độ laser tăng lên, thế năng của nguyên tử bị biến dạng mạnh, tạo ra một rào thế. Electron có thể thoát ra thông qua ion hóa đường hầm (xuyên ngầm) hoặc vượt qua rào thế hoàn toàn. Hiện tượng ion hóa kép đặc biệt thu hút sự chú ý vì nó cung cấp thông tin quý giá về tương quan electron-electron, một trong những tương tác cơ bản nhưng khó mô tả nhất trong vật lý. Việc nghiên cứu quá trình này trên các khí hiếm như Argon giúp làm sáng tỏ các động lực học electron cực nhanh, diễn ra trong khoảng thời gian attosecond (10⁻¹⁸ giây).

1.1. Hiểu đúng về cơ chế ion hóa đa photon và xuyên hầm

Trong vật lý laser, có hai cơ chế ion hóa chính khi cường độ trường laser chưa đạt đến mức cực đại. Thứ nhất là ion hóa đa photon (Multiphoton Ionization), xảy ra khi cường độ laser yếu hơn trường Coulomb nội tại của nguyên tử. Trong cơ chế này, nguyên tử hấp thụ đồng thời hoặc liên tiếp nhiều photon để electron có đủ năng lượng vượt qua thế ion hóa. Quá trình này thuộc vùng nhiễu loạn của quang học phi tuyến. Thứ hai là ion hóa đường hầm (Tunneling Ionization), hay còn gọi là ion hóa xuyên ngầm. Khi cường độ laser đủ lớn, nó làm biến dạng mạnh thế năng của nguyên tử, hạ thấp và làm hẹp rào thế Coulomb. Theo nguyên lý cơ học lượng tử, electron có một xác suất nhất định để "xuyên" qua rào thế này và trở thành tự do, ngay cả khi năng lượng của nó thấp hơn đỉnh rào. Nếu cường độ tiếp tục tăng, rào thế có thể bị hạ thấp xuống dưới mức năng lượng liên kết, cho phép electron thoát ra một cách "cổ điển" - đây gọi là cơ chế ion hóa vượt rào.

1.2. Phân biệt ion hóa kép tuần tự và không tuần tự NSDI

Quá trình ion hóa kép được chia thành hai loại chính. Ion hóa kép tuần tự (Sequential Double Ionization - SDI) là một quá trình hai bước: electron thứ nhất bị ion hóa, tạo thành ion Ar⁺, sau đó trường laser tiếp tục tác động và ion hóa nốt electron thứ hai để tạo thành Ar²⁺. Quá trình này tương đối đơn giản. Ngược lại, ion hóa kép không liên tục (Non-sequential Double Ionization - NSDI) là một cơ chế phức tạp hơn và thu hút nhiều sự quan tâm. Theo "mô hình ba bước" kinh điển, electron thứ nhất sau khi được giải phóng qua ion hóa đường hầm, sẽ được gia tốc bởi trường laser. Khi trường laser đổi chiều, electron này có thể bị kéo ngược lại và va chạm với ion mẹ. Vụ tái va chạm này truyền đủ năng lượng để đánh bật electron thứ hai. Cơ chế NSDI cung cấp một cửa sổ độc đáo để nghiên cứu tương quan electron-electron và động lực học electron trong thời gian thực.

II. Các thách thức trong mô phỏng Ion hóa kép Argon hiện nay

Việc mô phỏng chính xác quá trình ion hóa kép không liên tục (NSDI) của Argon là một thách thức lớn đối với vật lý tính toán. Lý do chính nằm ở sự phức tạp của hệ nhiều hạt tương tác trong một trường điện từ mạnh và biến đổi nhanh. Phương pháp tiếp cận chính xác nhất về mặt lý thuyết là giải phương trình Schrödinger phụ thuộc thời gian (TDSE), tuy nhiên, việc giải phương trình này cho một hệ có nhiều hơn một electron là cực kỳ tốn kém về tài nguyên tính toán và thường chỉ khả thi cho các hệ đơn giản như Heli hoặc trong các mô hình giảm chiều. Một thách thức lớn khác là mô tả đúng đắn hiệu ứng tương quan electron-electron. Lực đẩy giữa hai electron đóng vai trò quyết định trong cơ chế NSDI, nhưng lại rất khó để đưa vào các mô hình lượng tử một cách chính xác mà không làm bài toán trở nên quá phức tạp. Các mô hình cổ điển, mặc dù hiệu quả hơn về mặt tính toán, lại gặp khó khăn trong việc mô tả các hiệu ứng lượng tử thuần túy như ion hóa đường hầm. Do đó, việc tìm kiếm một phương pháp dung hòa được độ chính xác và hiệu quả tính toán, có khả năng tái tạo các kết quả thực nghiệm về phổ năng lượng electron và sự phân bố động lượng, vẫn là một bài toán mở và hấp dẫn trong cộng đồng vật lý nguyên tử.

2.1. Độ phức tạp của phương trình Schrödinger phụ thuộc thời gian

Nền tảng của cơ học lượng tử để mô tả động lực học electron là phương trình Schrödinger phụ thuộc thời gian (TDSE). Về nguyên tắc, việc giải chính xác TDSE sẽ cung cấp mọi thông tin về hệ. Tuy nhiên, trong thực tế, đây là một phương trình vi phân đạo hàm riêng cực kỳ phức tạp. Đối với nguyên tử Argon có hai electron hóa trị tham gia vào quá trình, việc giải TDSE đòi hỏi phải rời rạc hóa không gian và thời gian trên một lưới tính toán khổng lồ. Điều này dẫn đến chi phí bộ nhớ và thời gian tính toán tăng theo cấp số nhân, khiến phương pháp này không thực tế cho các khảo sát quy mô lớn. Hơn nữa, kết quả từ TDSE thường là một hàm sóng phức tạp, việc trích xuất các thông tin vật lý trực quan như quỹ đạo của electron là không thể, gây khó khăn cho việc diễn giải cơ chế vật lý.

2.2. Vấn đề tương quan electron electron trong thực nghiệm

Hiệu ứng tương quan electron-electron là trái tim của hiện tượng NSDI. Chính sự tương tác giữa electron tái va chạm và electron còn lại đã dẫn đến việc giải phóng electron thứ hai. Tuy nhiên, việc đo lường và mô hình hóa tương tác này là một thách thức. Trong thực nghiệm, các nhà khoa học đo phổ năng lượng electron hoặc sự phân bố động lượng tương quan (CTEMD) để suy ra các cơ chế tương tác. Các kết quả thực nghiệm cho thấy sự phân bố này phụ thuộc rất nhạy vào các thông số của laser. Ví dụ, nghiên cứu của Yueming Zhou và cộng sự (2015) chỉ ra rằng CTEMD của Argon thay đổi mạnh mẽ từ dạng đối xứng sang bất đối xứng khi tăng cường độ laser. Việc tái tạo lại những đặc điểm tinh vi này đòi hỏi một mô hình lý thuyết có khả năng nắm bắt chính xác động lực học electron và lực đẩy Coulomb tức thời, một yêu cầu mà nhiều phương pháp đơn giản hóa không thể đáp ứng.

III. Phương pháp Runge Kutta Giải mã động lực học electron

Để vượt qua những thách thức của các phương pháp lượng tử, các mô hình bán cổ điển đã được phát triển, trong đó phương pháp giải số đóng vai trò trung tâm. Thuật toán Runge-Kutta, đặc biệt là Runge-Kutta bậc 4 (RK4), là một công cụ cực kỳ mạnh mẽ và phổ biến để giải phương trình vi phân mô tả chuyển động của các hạt. Trong bối cảnh nghiên cứu ion hóa kép Argon, thuật toán này được áp dụng để giải các phương trình chuyển động Newton cho từng electron, chịu tác dụng của lực từ hạt nhân, lực đẩy từ electron kia và lực từ trường laser cường độ cao. Ưu điểm lớn của phương pháp số này là hiệu quả tính toán cao và độ chính xác có thể kiểm soát được. Luận văn của tác giả Trần Ngọc Liên Hương đã kiểm chứng độ tin cậy của thuật toán bằng cách so sánh kết quả mô phỏng một hệ dao động tắt dần với nghiệm giải tích chính xác. Kết quả cho thấy khi sử dụng đủ số bước nhảy (N=1000), sai số của thuật toán RK4 chỉ khoảng 0,01%, chứng tỏ độ tin cậy rất cao. Điều này cho phép áp dụng thuật toán để thực hiện các mô phỏng số quy mô lớn với hàng triệu nguyên tử, từ đó thu được các kết quả thống kê có ý nghĩa về xác suất ion hóa và phân bố động lượng.

3.1. Nguyên lý cơ bản của thuật toán Runge Kutta bậc 4 RK4

Thuật toán Runge-Kutta bậc 4 (RK4) là một phương pháp số lặp để tìm nghiệm gần đúng của phương trình vi phân bậc nhất với một giá trị ban đầu cho trước (bài toán Cauchy). Thay vì chỉ sử dụng đạo hàm tại điểm đầu của một khoảng như phương pháp Euler, RK4 tính toán giá trị của đạo hàm tại bốn điểm khác nhau trong khoảng đó (điểm đầu, hai điểm ở giữa và điểm cuối) và lấy trung bình có trọng số của chúng để ước tính giá trị tiếp theo. Công thức tính toán của RK4 được thiết kế để khớp với chuỗi Taylor đến bậc 4, giúp giảm thiểu sai số cục bộ. Cụ thể, sai số trên mỗi bước tỉ lệ với h⁵, trong đó h là kích thước bước nhảy. Điều này làm cho RK4 trở thành một lựa chọn cân bằng tuyệt vời giữa độ chính xác và sự đơn giản trong lập trình, vượt trội hơn các phương pháp bậc thấp hơn như Euler.

3.2. Áp dụng RK4 để giải các phương trình vi phân chuyển động

Trong bài toán ion hóa kép Argon, hệ cần giải là một tập hợp các phương trình vi phân bậc hai của Newton (F=ma). Để áp dụng thuật toán RK4 (vốn được thiết kế cho phương trình bậc nhất), ta cần chuyển đổi hệ phương trình bậc hai thành một hệ lớn hơn gồm các phương trình bậc nhất. Điều này được thực hiện bằng cách đặt các biến mới, ví dụ, vận tốc v = dx/dt. Khi đó, phương trình a = F/m trở thành dv/dt = F/m. Bằng cách này, một phương trình bậc hai cho vị trí được chuyển thành hai phương trình bậc nhất cho vị trí và vận tốc. Thuật toán RK4 sau đó được áp dụng đồng thời cho tất cả các phương trình này để cập nhật vị trí và vận tốc của các electron sau mỗi bước thời gian nhỏ. Phương pháp này cho phép theo dõi quỹ đạo của từng electron trong suốt quá trình tương tác laser-vật chất.

IV. Hướng dẫn mô phỏng Ion hóa kép bằng mô hình cổ điển

Mô hình được sử dụng trong nghiên cứu này là mô hình tập hợp ba chiều cổ điển (classical three-dimensional ensemble model), một công cụ hiệu quả trong vật lý tính toán để nghiên cứu các quá trình trong trường mạnh. Mô hình này theo dõi quỹ đạo cổ điển của hai electron hóa trị của Argon bằng cách giải phương trình vi phân Newton. Các lực tác động bao gồm lực hút từ hạt nhân, lực đẩy Coulomb giữa hai electron, và lực từ trường laser ngoài. Một ưu điểm lớn của mô hình này là nó cung cấp một hình ảnh trực quan về động lực học electron, cho phép phân tích trạng thái của hệ tại bất kỳ thời điểm nào. Để bắt đầu mô phỏng số, một tập hợp lớn (hàng triệu) nguyên tử Argon được khởi tạo. Vị trí và động lượng ban đầu của các electron được phân bố ngẫu nhiên sao cho năng lượng của hệ tương ứng với năng lượng trạng thái cơ bản của Argon. Sau đó, hệ được để tiến hóa tự do trong một khoảng thời gian ngắn để đạt trạng thái ổn định. Cuối cùng, xung laser được "bật" lên và thuật toán Runge-Kutta bậc 4 được sử dụng để tính toán quỹ đạo của từng electron. Kết thúc mô phỏng, năng lượng của mỗi electron được kiểm tra để xác định xem sự kiện ion hóa kép có xảy ra hay không.

4.1. Thiết lập mô hình tập hợp ba chiều cổ điển cho Argon

Việc thiết lập mô hình đòi hỏi xác định các thế tương tác. Lực tương tác giữa electron và hạt nhân, cũng như giữa hai electron, về cơ bản là lực Coulomb. Tuy nhiên, thế Coulomb có một điểm kỳ dị tại r=0, gây ra sự bất ổn định trong mô phỏng số. Để giải quyết vấn đề này và tránh hiện tượng tự ion hóa phi vật lý, một thế Coulomb "làm mềm" (soft-core Coulomb potential) được sử dụng. Thế này có dạng -Z/sqrt(r² + a²) cho tương tác electron-hạt nhân và 1/sqrt(r₁₂² + b²) cho tương quan electron-electron, trong đó a và b là các tham số làm mềm. Các tham số này được lựa chọn cẩn thận để năng lượng trạng thái cơ bản của mô hình khớp với thế ion hóa thực nghiệm của nguyên tử Argon. Trường laser thường được mô tả dưới dạng hình thang, bao gồm giai đoạn bật, giai đoạn ổn định và giai đoạn tắt, để mô phỏng một xung laser thực tế.

4.2. Phân tích kết quả Xác suất ion hóa và phổ động lượng

Sau khi mô phỏng kết thúc, các dữ liệu quỹ đạo được phân tích. Tiêu chí chính để xác định một sự kiện ion hóa kép là cả hai electron đều có năng lượng dương vào cuối xung laser. Bằng cách đếm số lượng các sự kiện này trên toàn bộ tập hợp, ta có thể tính được xác suất ion hóa. Thông tin quan trọng hơn được trích xuất từ động lượng cuối cùng của các electron. Bằng cách vẽ biểu đồ động lượng của electron thứ nhất theo động lượng của electron thứ hai (dọc theo phương phân cực của laser), chúng ta thu được phổ động lượng tương quan hai electron (CTEMD). Cấu trúc của phổ này, chẳng hạn như sự tập trung ở các góc phần tư khác nhau hay sự xuất hiện của các cấu trúc hình chữ V, tiết lộ thông tin chi tiết về cơ chế chia sẻ năng lượng và động lực học electron trong quá trình tái va chạm.

V. Kết quả nghiên cứu Ion hóa kép Argon Bằng chứng thực tiễn

Nghiên cứu áp dụng thuật toán Runge-Kutta trên mô hình cổ điển đã tái tạo thành công nhiều đặc điểm quan trọng của quá trình ion hóa kép không liên tục của Argon, cho thấy sự phù hợp với các kết quả thực nghiệm. Một trong những phát hiện chính là sự phụ thuộc mạnh mẽ của quá trình ion hóa vào các thông số của xung laser. Khi thay đổi độ dài xung, kết quả mô phỏng số chỉ ra rằng hầu hết các sự kiện ion hóa kép xảy ra trong vài chu kỳ đầu tiên khi cường độ laser ở trạng thái ổn định. Điều này cho thấy có hiện tượng bão hòa: sau một vài chu kỳ, các nguyên tử dễ bị ion hóa nhất đã bị ion hóa. Quan trọng hơn, nghiên cứu đã khảo sát sự phụ thuộc của phổ năng lượng electron vào cường độ laser. Ở cường độ yếu (ví dụ 1x10¹⁴ W/cm²), phổ động lượng có tính đối xứng cao, cho thấy hai electron chia sẻ năng lượng một cách tương đối đồng đều. Nhưng khi cường độ tăng lên (ví dụ 2x10¹⁴ W/cm² trở lên), phổ trở nên bất đối xứng rõ rệt, tập trung ở góc phần tư thứ nhất và thứ ba. Sự bất đối xứng này là dấu hiệu của việc electron tái va chạm truyền rất ít năng lượng cho electron còn lại, phù hợp với lý thuyết về sự chia sẻ năng lượng bất đối xứng (AES).

5.1. Phân tích phổ năng lượng electron phụ thuộc độ dài xung

Kết quả khảo sát cho thấy khi giữ nguyên cường độ laser và tăng độ dài của giai đoạn ổn định trong xung (ví dụ từ 2 chu kỳ lên 4, rồi 6 chu kỳ), phổ động lượng tương quan hai electron ban đầu thay đổi nhưng sau đó nhanh chóng đạt trạng thái ổn định. Điều này chứng tỏ rằng quá trình ion hóa kép không liên tục của Argon chủ yếu xảy ra trong một khoảng thời gian ngắn ngay khi laser đạt cường độ đỉnh. Sau giai đoạn này, xác suất ion hóa thêm là rất thấp. Phát hiện này củng cố giả thuyết về sự bão hòa của quá trình ion hóa và cho thấy NSDI là một hiện tượng đặc biệt hiệu quả đối với các xung laser ngắn, phù hợp với các quan sát thực nghiệm.

5.2. Ảnh hưởng của cường độ trường laser lên quá trình NSDI

Đây là một trong những kết quả quan trọng nhất. Ở cường độ laser yếu, electron tái va chạm có vận tốc thấp, thời gian tương tác với ion mẹ đủ dài để xảy ra sự chia sẻ năng lượng đối xứng (Symmetric Energy Sharing - SES), dẫn đến phổ động lượng đối xứng. Ngược lại, trong trường laser cường độ cao, electron tái va chạm có vận tốc rất lớn. Thời gian tương tác trở nên cực ngắn, khiến nó gần như không truyền năng lượng cho electron liên kết. Kết quả là electron tái va chạm bay ra với động lượng lớn, trong khi electron thứ hai bị "lắc" ra với động lượng nhỏ. Đây chính là cơ chế chia sẻ năng lượng bất đối xứng (Asymmetric Energy Sharing - AES), tạo ra cấu trúc bất đối xứng trong phổ động lượng. Sự chuyển đổi từ SES sang AES khi tăng cường độ laser là một bằng chứng mạnh mẽ cho thấy mô hình cổ điển với thuật toán Runge-Kutta đã nắm bắt được bản chất vật lý của tương quan electron-electron.

VI. Tương lai nghiên cứu Ion hóa kép Argon và các ứng dụng

Nghiên cứu ion hóa kép Argon bằng thuật toán Runge-Kutta và mô hình cổ điển đã chứng tỏ là một hướng đi đầy hứa hẹn. Thành công trong việc mô phỏng các đặc điểm chính của phổ động lượng không chỉ xác nhận độ tin cậy của phương pháp số mà còn mở ra nhiều hướng phát triển mới. Việc hiểu rõ cơ chế tương quan electron-electron là nền tảng cho nhiều lĩnh vực tiên tiến. Ví dụ, nó giúp kiểm soát các phản ứng hóa học ở cấp độ phân tử bằng các xung laser định hình, hoặc tạo ra các nguồn phát tia X và tia cực tím (XUV) siêu ngắn thông qua quá trình phát sóng hài bậc cao (HHG). Trong tương lai, các mô phỏng số có thể được mở rộng để nghiên cứu các hệ phức tạp hơn, chẳng hạn như các phân tử (N₂, O₂) hoặc các nguyên tử nặng hơn. Việc kết hợp các yếu tố lượng tử vào mô hình cổ điển, như hiệu ứng ion hóa đường hầm, cũng là một hướng đi quan trọng để tăng độ chính xác. Cuối cùng, sự kết hợp giữa các mô phỏng tiên tiến và các thực nghiệm với công nghệ laser attosecond sẽ tiếp tục làm sáng tỏ những bí ẩn của thế giới vi mô, thúc đẩy sự phát triển của khoa học vật liệu, hóa học lượng tử và công nghệ nano.

6.1. Tổng kết những phát hiện chính từ mô phỏng số

Nghiên cứu đã thành công trong việc: (1) Kiểm chứng độ chính xác cao của thuật toán Runge-Kutta bậc 4 cho các bài toán động lực học. (2) Mô phỏng lại thành công sự phụ thuộc của quá trình NSDI vào độ dài xung laser, cho thấy hiện tượng bão hòa. (3) Tái tạo được sự chuyển đổi từ cơ chế chia sẻ năng lượng đối xứng sang bất đối xứng khi tăng cường độ trường laser cường độ cao. Những kết quả này khẳng định rằng mô hình tập hợp cổ điển, dù đơn giản hơn TDSE, vẫn là một công cụ mạnh mẽ và trực quan để khám phá vật lý nguyên tử trong trường mạnh.

6.2. Hướng phát triển Mở rộng sang các phân tử phức tạp

Một hướng phát triển tự nhiên là áp dụng phương pháp này để nghiên cứu ion hóa kép của các phân tử. Khác với nguyên tử có tính đối xứng cầu, phân tử có cấu trúc phức tạp hơn với nhiều tâm hạt nhân và các bậc tự do về dao động, quay. Điều này làm cho động lực học electron trở nên phong phú và đa dạng hơn nhiều. Việc mô phỏng sự tương tác laser-vật chất trên các phân tử như H₂ hay N₂ sẽ giúp làm rõ vai trò của cấu trúc phân tử đối với quá trình ion hóa, một vấn đề còn nhiều tranh cãi. Xa hơn nữa, việc kết hợp các thuật toán học máy (machine learning) để tối ưu hóa việc tìm kiếm các tham số mô hình hoặc phân tích các bộ dữ liệu lớn từ mô phỏng số cũng là một hướng đi đầy tiềm năng.

Trích đoạn nội dung tài liệu

Chương 1. Trình bày về tương tác giữa laser với vật chất, trong đó nhấn mạnh vào quá trình ion hóa của nguyên tử, phân tử thông qua ba cơ chế là ion hóa đa photon, ion hóa xuyên ngằm và ion hóa vượt rao. Ngoài ra trong chương này còn đi sâu tìm hiểu về quá trình ion hóa kép không liên tục của nguyên tử, phân tử. Cơ sở lý thuyết và phương pháp tính.

Trình bày về mô hình tập hợp ba chiều cé điển [3]. Đây là một công cụ toán học hiệu quả để khảo sat các quá trình vật lý xảy ra khi vật chất chịu tác dụng của trường laser. Ngoài ra trong chương này tác giả còn dé cập đến thuật toán Runge — Kutta dé giải quyết bài toán Cauchy cơ bản nhất, từ đó mở rộng khảo sát cho phương trình vi phân cắp cao. Ngoài ra trong chương này còn cho thấy mức độ tin cậy cao của thuật toán Runge — Kutta so với phương pháp giải tích phương trình vi phân cấp hai trong bài toán cổ điển Newton.

Kết quả nghiên cứu. Kiểm chứng tinh chính xác của chương trình tinh toán thông qua việc khảo sát cho hệ dao động tắt dần. Đồng thời trong chương này còn trình bày kết quả khảo sat sự ion hóa kép không liên tục của nguyên tử Argon dưới tác dụng của trường laser khi thay đổi độ dài xung và cường độ điện trường. Quá trình tương tác giữa laser và nguyên tử, phân tử Hiện nay, quá trình tương tác giữa trường laser với nguyên tử, phân tử là một trong những đề tài hap dẫn thu hút nhiều sự quan tâm của các nhà khoa học bởi nó cỏ nhiều ứng dụng trong khoa học công nghệ và trong cuộc sống.

Theo lý thuyết lượng tử, trường laser được xem như là những dòng hạt photon có năng lượng, động lượng và spin xác định. Vì vậy khi tương tac với nguyên tử phân tử, laser sẽ làm cho các trạng thái nguyên tử bị thay đổi. Khi đó nhiều hiệu img phi tuyến được xảy ra như quá trình ion hóa, phát xạ sóng điều hòa bậc cao, quá trình ion hóa kép không liên tục. Tùy thuộc vào cường độ của trường laser mà cơ ché tương tác giữa laser với nguyên tử, phân tử sẽ khác nhau.

Quá trình ion hóa Khi cường độ trường laser yếu hơn nhiều so với thế ion hóa của nguyên tử, trường laser chỉ có tác dụng gây ra sự nhiễu loạn các trạng thái electron của nguyên tử. Trong trường hợp nay, các mức năng lượng của nguyên tử bị thay đổi tỉ lệ với bình phương cường độ của trường laser, hiện tượng nay gọi là sự dich chuyển Stark. Do đó, vùng này được gọi là vùng nhiễu loạn của quang học phi tuyến. Trong vùng nay, sự ion hóa chủ yếu diễn ra theo cơ chế đa photon, nghĩa là nguyên tử hấp thụ liên tiếp nhiều photon làm cho năng lượng của nó tăng dần đến khi lớn hơn năng lượng liên kết thi electron được chuyển sang trạng thái tự do.

Như vậy, khi cường độ trường laser yếu hơn nhiều so với trường Coulomb của nguyên tử thì nguyên tử chỉ hấp thụ một cách tự phát N photon và xảy ra sự ion hóa đa photon. Trong trường hợp ion hóa đa photon thế của nguyên từ không bị biển dạng do điện trường không đủ mạnh. Mutiphoton tonization —_ Hình 1. Sự tan hóa da photon [13] Khi cường độ trường laser lớn hon so với trường Coulomb của nguyên tử.

Khi đỏ thé của nguyên tử sẽ bị biến dạng. Lúc đó các electron liên kết yếu với hạt nhân va sẽ thoát khỏi nguyên tử hoặc phân tử thông qua một rao thé xác định bởi cơ chế xuyên ngằm. Rao thé này trở nên thắp và móng hơn khi cường độ trường laser tăng lên.2a, đường mảnh ứng với thé của điện trường, đường cong day ứng với thé năng hiệu dụng và đường nằm ngang đặc trưng cho năng lượng liên kết của nguyên tử khi không có trường laser. Cơ chế ion hóa: xuyên ngắm (a), vượt rào (b) [4] Khi điện trường tiếp tục tăng rào thế sẽ bị hẹp lại và hạ thấp xuống, đến khi nó hạ xuống thấp hơn năng lượng liên kết của electron thì electron sẽ đi từ trạng thái liên kết sang trạng thái tự do.

Sự ion hóa lúc này điển ra theo cơ chế vượt rào (hình 1. Cường độ trường laser tới hạn F; xác định giao điểm giữa vùng ion hóa xuyên ngầm và vùng ion hóa vượt rao. Giá trị Fe này có thé được ước lượng bằng cách cho thé năng hiệu dụng cực đại bằng với thé năng ion hóa của electron liên kết. Quá trình ion hóa kép Như chúng ta đã biết, khi chiếu một chùm laser vào trong nguyên tử thì các electron trong nguyên tử có thé bị ion hóa theo các cơ chế như ion hóa da photon, xuyên ngầm hay vượt rào tùy thuộc vảo cường độ của chùm laser.

Trong đó khi các electron trong nguyên tử bị ion hóa xuyên ngâm qua rào thé tạo bởi thé ion hóa của nguyễn tử và cường độ chim laser sẽ có nhiều hiện tượng phi tuyến xảy ra. Khi đó electron tự do bứt ra khỏi ion mẹ được gia tốc trong trường laser, nó quay lại va chạm với ion mẹ khi trường laser đổi chiều. Quá trình ion hóa kép là quá trình hai electron bị bứt ra khỏi ion mẹ đưới tác dụng của trường laser khi bức xạ của trường laser tác dụng vào nguyên tử hay phân tử trung hòa. Quá trình ion hóa kép được chia thành hai cơ chế khác nhau: quá trình ion hóa kép liên tục và quá trình ion hóa kép không liên tục.

Quá trình ion hóa kép liên tục Trong đó quá trình ion hóa kép liên tục được hiểu là cả hai electron trong nguyên tử birt ra cùng một thời điểm dudi tác dụng của trường laser hoặc sau một thời gian electron thứ nhất bứt ra thi electron thứ hai cũng bứt ra khỏi ion mẹ. Trong quá trình ion hóa kép liên tục không xảy ra sự tái va chạm của electron và ion mẹ. Tu etectron (ián997x 9m Ww One 441260) (O4 w w ® lente ground state “@ w Newerat stom \ growed van Hình 1. Vi dụ về quá trình ion hóa kép liên tục [14] 1.

Quá trình ion hóa kép không liên tục Còn quá trình ion hóa kép không liên tục có cơ chế khác hoàn toàn, được hiểu đó là quả trinh electron thứ nhất bit khỏi ion mẹ, sau một thời gian khi trường laser đổi chiều, nó quay lại va chạm với ion mẹ, khi đó electron thứ hai mới được giải phóng. Khi 46 động năng cực đại dé electron quay lại va chạm với ion mẹ được tinh là E=3.17U, , do đó năng lượng của photon là: £SI,+3.1) trong đó /, lả thé ion hóa của nguyên tử, U, lả thé truyền động trong trường laser được xác định theo công thức U,, =//4ø` (trong đỏ !, @ là cường độ va tan số của trường laser). Hiện tượng ion hóa kép không liên tục được phát hiện bằng thực nghiệm bởi Suran và Zapesochny cho nguyên tử kiểm thé vào đâu năm 1975. Mặc đủ được nghiên cứu sâu 10 rong nhưng chi tiết về quá trình ion hóa kép của nguyên tử kiểm thé này vẫn chưa được tìm hiểu kỹ.

Quá trình ion hóa kép trong trường hợp này được thực hiện bởi sự chuyển mức của cả hai electron thông qua phố trạng thải nguyên tử, nằm giữa thé cua electron thứ nhất và thứ hai Đổi với các nguyên tử khí hiếm, quá trình ion hóa kép không liên tục đầu tiên được quan sát bởi L'Huillier. Hiện tượng này nhanh chóng thu hút sự quan tâm của các nhà nghiên cứu sau khi nó được tìm thấy trong trường hồng ngoại và đối với cường độ cao hơn. Cơ chế của quá trình ion hóa kép không liên tục của các nguyên tử khi hiểm khác với các nguyên tử kiểm thổ. Sau khi electron thứ nhất bị ion hóa, electron được giải phóng có thé quay lại va chạm với ion mẹ.

Electron nảy đóng vai trỏ như một “atomic antenna”, hap thu năng lượng tir trường laser giữa quá trình ion hóa và quá trình tải va chạm, cơ chế này được gọi là mô hình ba bước của quá trình ion hóa kép không liên tục [5].4, M6 hình ba bước cua NSDI [15] Trong tương tác giữa nguyên tử với laser, quá trình ion hóa kép không liên tục thu hút nhiều sự quan tâm của các nhà khoa học bởi những tín hiệu ghi nhận được từ phd động lượng tương quan của hai electron cung cấp rất nhiều thông tin về sự tương tác giữa laser với vật chất cũng như sự tương quan giữa electron — electron trong nguyên tử khi electron thứ nhất quay lại va chạm với electron thứ hai thông qua va chạm không đàn hỏi. CƠ SỞ LÝ THUYET VÀ PHƯƠNG PHÁP TÍNH 2. Mô hình tập hợp ba chiều cổ điển “Trong tương tác giữa nguyên tử, phân tử với trường laser, quá trình ion hóa kép không liên tục thu hút nhiều sự quan tâm của các nhà khoa học bởi những tín hiệu ghi nhận được từ phô động lượng tương quan của hai electron cung cấp rat nhiều thông tin về sự tương tác giữa laser với vật chất cũng như sự tương quan giữa electron — electron trong nguyên tử khi trải qua quá trình ion hóa kép. Hiện nay có hai cách phổ biến đẻ tiếp cận bai toán NSDI.

Cách thứ nhất là dựa theo nén tảng của cơ học lượng tử, trong đó những tính chất của các electron ion hóa được thu nhận từ việc sử dụng phương pháp giải chính xác nghiệm phương trình Schrödinger phụ thuộc thời gian (TDSE — Time Dependent Schrödinger Equation) với sự góp mặt của trường laser. Trong phương pháp giải số, đây là phương pháp cho kết quả chính xác. Tuy nhiên, với phương pháp này việc lập trình khá phức tạp và chi cho kết quả cuối cùng mà không cung cắp các thông tín vẻ quá trình động lực học vật lý của các electron dién ra trong sudt quá trình tương tác với laser. Vì vậy phương pháp này không được sử dụng phổ biến để khảo sát bai toán NSDI.

Cách tiếp cận thứ hai là sử dụng mô hình tập hợp ba chiều cé điển thông qua việc giải phương trình Newton cho từng electron chịu sự tac dụng của trường laser va trường tương tac hạt nhân của nguyên tử. Từ khi được giới thiệu vào năm 2001 [3], mô hình tập hợp ba chiều này được xem là một công cụ hiệu quả để nghiên cứu quá trình DI trong trường mạnh. Với phương pháp này, vai trò của lực đẩy electron — electron và lực hút electron ~ ion trong việc quan sát phổ electron của quá trình ion hóa kép không liên tục cũng được xác định một cách thành công và năng lượng bắt đối xứng trong quá trình tái va chạm cũng được phát hiện.