Tổng quan nghiên cứu

Quá trình phát xạ sóng điều hòa bậc cao (High-order Harmonic Generation - HHG) là một trong những hiệu ứng quang phi tuyến nổi bật khi cho chùm bức xạ laser cường độ cao trong khoảng $10^{13}$ đến $10^{16}$ W/cm2 tương tác với môi trường nguyên tử khí hiếm. Nhờ sự phát triển của công nghệ quang tử, hiệu suất chuyển đổi phát xạ HHG đã được cải thiện từ mức $10^{-7}$ lên khoảng $5 \times 10^{-5}$, giúp chùm bức xạ này trở thành công cụ đắc lực để tạo các chuỗi xung atto giây ($10^{-18}$ giây) và theo dõi động học electron với độ phân giải không gian dưới 1 Angstrom ($10^{-10}$ mét). Tuy nhiên, rào cản kỹ thuật lớn nhất hiện nay là việc tạo ra nguồn bức xạ HHG có độ phân cực elip hoặc phân cực tròn ở vùng năng lượng photon cao thuộc miền tử ngoại chân không và tia X mềm.

Vấn đề cốt lõi của nghiên cứu là khảo sát và kiểm soát trạng thái phân cực elip của chùm bức xạ HHG sinh ra khi chiếu chùm laser hai màu phân cực thẳng vuông góc vào nguyên tử Heli. Mục tiêu trọng tâm bao gồm việc xác định quy luật biến thiên của độ elip theo ba tham số điều khiển đầu vào: tỷ số cường độ chùm sáng thứ cấp trên sơ cấp ($I_2/I_1$), độ lệch pha tương đối ($\Delta\varphi$) và độ trễ thời gian bật xung ($\Delta t$). Toàn bộ mô phỏng số được thực hiện trên hệ thống máy tính hiệu năng cao tại Phòng thí nghiệm Vật lý tính toán với bước sóng laser cơ bản 800 nm trong năm 2018. Kết quả nghiên cứu có ý nghĩa quyết định trong việc tối ưu hóa độ elip cực đại lên mức 78% đến 80%, cung cấp cơ sở lý thuyết vững chắc cho các phép đo phổ hấp thụ lưỡng sắc từ tia X (XMCD) với dải năng lượng photon đạt trên 70 eV.

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

Khung lý thuyết áp dụng

Cơ sở phân tích tương tác quang phi tuyến trong nghiên cứu dựa trên mô hình ba bước bán cổ điển được phát triển bởi Corkum (1993) và Lewenstein (1994). Quá trình phát xạ HHG diễn ra tuần tự qua ba giai đoạn: ion hóa xuyên hầm electron qua rào thế nguyên tử dưới tác động của điện trường laser mạnh cỡ $5,1 \times 10^{11}$ V/m, electron tự do nhận gia tốc dọc theo quỹ đạo dao động trong trường điện từ ngoài, và cuối cùng electron quay trở lại va chạm tái kết hợp với ion mẹ để phát ra photon mang năng lượng cao. Hình thái phân cực của bức xạ phát ra phụ thuộc trực tiếp vào quỹ đạo chuyển động hai chiều của electron trước khi tái kết hợp.

Độ phân cực elip ($\epsilon$) của từng bậc sóng điều hòa được định lượng toán học thông qua tỷ số cường độ thành phần điện trường $I_y/I_z$ và độ lệch pha $\varphi_y - \varphi_z$ giữa hai phương phân cực vuông góc Oy và Oz. Giá trị độ elip biến thiên liên tục trong khoảng từ -1 đến +1, trong đó giá trị 0 đại diện cho sóng phân cực thẳng hoàn toàn và giá trị tuyệt đối bằng 1 tương ứng với sóng phân cực tròn lý tưởng. Để sóng điều hòa đạt độ elip cao, điều kiện tiên quyết là tỷ số cường độ giữa hai thành phần điện trường trên hai trục tọa độ phải tiến sát về giá trị 1 và hiệu số pha dao động phải xấp xỉ 90 độ hoặc 270 độ.

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

Nghiên cứu sử dụng phương pháp giải số trực tiếp phương trình Schrödinger phụ thuộc thời gian (TDSE) kết hợp điều kiện xấp xỉ một electron hoạt động (Single Active Electron - SAE). Không gian tọa độ tính toán được thiết lập với bán kính giới hạn cực đại RMAX = 140 đơn vị nguyên tử (a.u.), số điểm lưới rời rạc NDVRR = 550 điểm, số lượng mômen xung lượng góc tối đa LMAX = 120 và số bước tích phân thời gian NTMAX = 3000 bước nhằm đảm bảo độ hội tụ số học tuyệt đối.

Nguyên tử Heli ở trạng thái cơ bản được lựa chọn làm môi trường tương tác nhờ sở hữu thế ion hóa cao đạt 24,58 eV, cho phép áp dụng chùm tia laser có cường độ đỉnh lên đến $6,0 \times 10^{14}$ W/cm2 mà không làm môi trường bị ion hóa bão hòa thành plasma. Trường laser kích thích sử dụng xung bao phẳng (flat-topped envelope) với 6 chu kỳ quang học ở bước sóng 800 nm (tương đương năng lượng photon 1,553 eV), giúp cấu trúc các đỉnh phổ phân tán sắc nét hơn hẳn so với các dạng xung hàm Gauss truyền thống. Phương pháp TDSE kết hợp xấp xỉ SAE được ưu tiên lựa chọn vì khắc phục triệt để sai số của các phép gần đúng trường mạnh đơn giản, đồng thời tối ưu hóa thời gian tính toán hơn 50% so với phương pháp Hartree-Fock đa cấu hình.

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

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

Khảo sát trên cấu hình chuẩn $\omega + 2\omega$ cho thấy sự phân tách tính chất quang học rõ rệt giữa hai phương phân cực vuông góc. Phổ phát xạ theo phương Oz (trùng hướng phân cực laser cơ bản) tập trung chủ yếu các bậc sóng điều hòa lẻ với mật độ cường độ vượt trội, trong khi phổ phát xạ theo phương Oy (hướng phân cực chùm sóng điều hòa bậc hai) lại bị chi phối bởi các bậc chẵn với cường độ thấp hơn từ 2 đến 3 bậc độ lớn logarit. Vị trí ngưỡng cắt năng lượng cực đại (cut-off) ghi nhận tại bậc 53, tương ứng với mức năng lượng photon 82 eV.

Độ elip của các bậc chẵn đạt giá trị rất cao, dao động từ 50% đến 80% khi duy trì tỷ số cường độ laser thứ cấp trên sơ cấp $I_2/I_1 = 0,13$. Ngược lại, độ elip của tất cả các bậc lẻ luôn duy trì dưới ngưỡng 20%. Khi tăng tỷ số $I_2/I_1$ từ 0,13 lên 0,50, cường độ phát xạ tổng thể tăng lên nhưng số lượng các bậc chẵn có độ elip trên 50% bị sụt giảm hơn 65%. Về yếu tố thời gian, việc duy trì độ trễ xung hẹp trong khoảng $|\Delta t| \le 0,5$ chu kỳ quang học giúp số lượng bậc có độ elip cao đạt mức cực đại, trong khi việc kéo dài độ trễ lên mức 2,5 chu kỳ làm suy giảm độ elip xuống dưới 30%.

Đối với cấu hình đề xuất mới $\omega + 3\omega$, tính đối xứng không gian bị thay đổi căn bản giúp các bậc sóng điều hòa lẻ xuất hiện đồng thời trên cả hai trục Oy và Oz với cường độ tương đương, trong khi các bậc chẵn bị triệt tiêu hoàn toàn. Khi tăng tỷ số cường độ $I_2/I_1$ từ 0,75 lên 3,00 (tương ứng cường độ $I_2 = 6,0 \times 10^{14}$ W/cm2), độ elip cực đại của các bậc lẻ tăng đều đặn và đạt mức đỉnh 78% tại bậc 27 (năng lượng photon xấp xỉ 42 eV). Cấu hình $\omega + 3\omega$ đã giải quyết thành công bài toán đánh đổi giữa cường độ bức xạ và độ tinh khiết phân cực.

Thảo luận kết quả

Sự khác biệt căn bản giữa hai cấu hình bắt nguồn sâu xa từ quy tắc chọn lọc đối xứng thời gian nửa chu kỳ ($T/2$). Dưới cấu hình $\omega + 2\omega$, véctơ điện trường tại hai thời điểm cách nhau nửa chu kỳ đối xứng qua trục Oy, khiến quỹ đạo gia tốc của electron làm lệch pha tổng hợp giữa hai chùm tín hiệu theo trục Oy một lượng bằng $N\pi$. Điều này dẫn tới việc các sóng bậc lẻ bị triệt tiêu hoàn toàn trên trục Oy, đẩy tỷ số cường độ $I_y/I_z$ của bậc lẻ về gần mức 0 và làm giảm tính elip.

Khi trình bày dữ liệu qua đồ thị phổ công suất quang học và biểu đồ phân bố độ elip theo bậc điều hòa, cấu hình $\omega + 3\omega$ thể hiện ưu thế vượt trội khi các đỉnh phổ phát xạ trên cả hai trục tọa độ tiến lại gần nhau. Tính đối xứng tâm qua gốc tọa độ của trường điện $\omega + 3\omega$ bảo đảm rằng các electron giải phóng sau mỗi nửa chu kỳ đều đóng góp vào sự giao thoa tăng cường của các bậc lẻ trên cả hai phương. Mức độ phân cực elip cao 78% tại bậc 27 khẳng định tính khả thi của việc tạo ra nguồn tia X mềm phân cực elip có độ sáng cao, hoàn toàn tương thích với các quan sát thực nghiệm quốc tế gần đây.

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

Dựa trên các phân tích định lượng từ mô hình tính toán số, nghiên cứu đề xuất bốn giải pháp kỹ thuật cụ thể nhằm tối ưu hóa quá trình tạo chùm bức xạ HHG phân cực elip:

  1. Thiết lập tỷ số cường độ chùm sáng tối ưu: Đối với hệ thống sử dụng cấu hình $\omega + 2\omega$, các nhóm nghiên cứu thực nghiệm cần cố định tỷ số cường độ $I_2/I_1$ trong dải hẹp từ 0,10 đến 0,15 để giữ độ elip của các bậc chẵn trên mức 75%. Khi chuyển sang cấu hình $\omega + 3\omega$, các kỹ sư cần điều chỉnh tỷ số $I_2/I_1 \approx 3,0$ nhằm tối đa hóa đồng thời cả cường độ bức xạ lẫn độ elip cho các bậc sóng lẻ.
  2. Kiểm soát độ trễ thời gian xung siêu chính xác: Các phòng thí nghiệm quang học cần lắp đặt hệ thống dịch pha giao thoa kế nhằm duy trì độ trễ thời gian giữa hai chùm xung trong ngưỡng $|\Delta t| \le 0,5$ chu kỳ quang học (tương đương sai số thời gian dưới 1,33 femtogiây đối với chùm laser 800 nm), ngăn chặn hiện tượng suy thoái độ elip do sự bất đối xứng xung.
  3. Mở rộng dải bước sóng laser dẫn đường: Đơn vị nghiên cứu nên chuyển đổi nguồn phát laser kích thích từ vùng cận hồng ngoại 800 nm sang vùng hồng ngoại trung từ 1300 nm đến 1800 nm trong lộ trình nâng cấp thiết bị 2 năm tới, giúp mở rộng ngưỡng cắt quang phổ từ 82 eV lên trên 160 eV.
  4. Hoàn thiện mô hình tính toán tương tác đa electron: Đội ngũ phát triển phần mềm tính toán cần mở rộng mã nguồn FORTRAN hiện tại để tích hợp hiệu ứng tương tác tương quan electron đầy đủ cho các nguyên tử khí hiếm nặng như Argon và Krypton, hoàn thành trước quý 3 năm 2027.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị thực tiễn cho bốn nhóm đối tượng chuyên môn:

  1. Nhà nghiên cứu vật lý quang học và khoa học atto giây: Tiếp cận toàn diện thuật toán giải số TDSE và bộ thông số chuẩn hóa, ứng dụng trực tiếp vào việc thiết kế các nguồn phát xung siêu ngắn phục vụ việc chụp ảnh cấu trúc orbital phân tử.
  2. Kỹ sư phát triển vật liệu nano và thiết bị lưu trữ từ tính: Khai thác quy trình tối ưu hóa bức xạ elip để nâng cao độ nhạy của kỹ thuật đo phổ XMCD, hỗ trợ xác định chính xác sự phân bố spin và mômen động lượng quỹ đạo của các màng mỏng dị thể.
  3. Học viên cao học và nghiên cứu sinh ngành Sư phạm Vật lý và Vật lý tính toán: Sử dụng công trình như một tài liệu tham khảo mẫu mực về kỹ thuật mô hình hóa toán lý, phương pháp rời rạc hóa lưới không gian và xử lý dữ liệu phổ quang học phi tuyến.
  4. Kỹ thuật viên vận hành hệ thống laser công suất cao: Ứng dụng các khuyến nghị về việc điều chỉnh tinh thể phi tuyến BBO và bộ tách pha phản xạ để kiểm soát độ tinh khiết phân cực chùm sáng trong thực tế.

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

Tại sao cấu hình laser $\omega + 2\omega$ chỉ tạo ra độ elip cao ở các bậc chẵn?

Do tính chất đối xứng qua trục Oy của điện trường tổng hợp tại các thời điểm cách nhau nửa chu kỳ, các sóng điều hòa bậc lẻ bị triệt tiêu hoàn toàn theo phương Oy. Khi đó, chỉ có các bậc chẵn mới xuất hiện đồng thời trên cả hai phương phân cực Oz và Oy, tạo điều kiện cho tỷ số cường độ $I_y/I_z$ đạt giá trị gần bằng 1 và hình thành độ elip cao tới 80%.

Xấp xỉ một electron hoạt động (SAE) có bảo đảm độ chính xác cho nguyên tử Heli không?

Xấp xỉ SAE hoàn toàn chính xác đối với nguyên tử Heli vì thế ion hóa lần thứ nhất của Heli rất lớn (24,58 eV), cao hơn rất nhiều so với năng lượng một photon laser 800 nm (1,553 eV). Xác suất ion hóa đồng thời hai electron là cực kỳ nhỏ, do đó việc bỏ qua tương tác đa electron không làm sai lệch bức xạ phát xạ chính.

Độ trễ thời gian $\Delta t$ ảnh hưởng như thế nào đến chất lượng chùm sáng?

Khi độ trễ thời gian vượt quá 0,5 chu kỳ quang học, vùng chồng phủ giữa hai xung bị thu hẹp, khiến nguyên tử tương tác độc lập với từng thành phần laser phân cực thẳng. Quá trình này làm gia tăng cường độ phát xạ phân cực thẳng đơn thuần và kéo độ elip của toàn bộ phổ quang học sụt giảm hơn 40%.

Cấu hình đề xuất $\omega + 3\omega$ có ưu điểm gì vượt trội so với cấu hình cũ?

Cấu hình $\omega + 3\omega$ phá vỡ sự đối xứng gương và tạo ra tính đối xứng tâm, giúp các bậc lẻ xuất hiện đồng thời trên cả hai phương dao động vuông góc. Người nghiên cứu có thể đồng thời gia tăng cả cường độ bức xạ lẫn độ elip lên mức 78% thông qua việc tăng tỷ số cường độ chùm laser thứ cấp $I_2/I_1$.

Bức xạ HHG phân cực elip đóng vai trò gì trong phép đo phổ XMCD?

Phép đo XMCD đòi hỏi nguồn tia X phân cực tròn hoặc elip để đo lường sự chênh lệch hấp thụ quang học giữa các trạng thái từ hóa của mẫu vật liệu. Nguồn HHG phân cực elip năng lượng cao trên 70 eV cung cấp chùm bức xạ kết hợp pha có độ phân giải thời gian atto giây, giúp khảo sát động học từ tính cực nhanh ở cấp độ nano.

Kết luận

  • Luận văn đã xây dựng thành công quy trình giải số chính xác phương trình Schrödinger phụ thuộc thời gian (TDSE) kết hợp xấp xỉ một electron hoạt động (SAE) cho nguyên tử Heli với 3000 bước thời gian tích phân.
  • Nghiên cứu đã định lượng hóa một cách có hệ thống sự phụ thuộc của độ elip bức xạ HHG vào tỷ số cường độ chùm sáng, hiệu pha tương đối và độ trễ thời gian xung.
  • Công trình đã chỉ rõ hạn chế cố hữu của cấu hình $\omega + 2\omega$ khi các bậc sóng điều hòa có độ elip cao (bậc chẵn) lại sở hữu mật độ năng lượng phát xạ rất thấp.
  • Đề xuất cấu hình cải tiến $\omega + 3\omega$ đã giải quyết triệt để rào cản trên, cho phép đạt độ elip cực đại 78% tại bậc 27 (năng lượng 42 eV) đồng thời với mức cường độ phát xạ cao.
  • Các đơn vị nghiên cứu và viện chuyên ngành cần tiếp tục triển khai các thử nghiệm thực tế trên hệ thống laser hồng ngoại trung để mở rộng dải phổ bức xạ phân cực elip lên trên 150 eV trong thời gian tới.