C Nguyên tử được chọn làm gốc Đánh số 1, lần lượt H-2, N-
ĐO HHG DO TƯƠNG TÁC LASER-PHÂN TỬ
Ta biết cường độ của LASER thay đổi tuần hoàn theo thời gian, mỗi vòng tuần hoàn ta gọi là một xung. Trong một xung, tại một thời điểm nào đó LASER đạt cường độ cực đại gọi là cường độ đỉnh.
Khi dùng xung LASER tương tác ta phải xét đến các yếu tố sau của xung LASER:
Cường độđỉnh
Độ dài xung
Tần số
Pha
Ởđây, ta chọn xung Laser có: Cường độđỉnh: 2.1014 W/cm2, độ dài xung: 30 fs, tần số: có năng lượng 1,551 eV, bước sóng: =800 nm, pha: =0
Bước 1:Chạy get_wf_xie-modifyO2.f là một Source Code viết trên chương trình FORTRAN do nhóm nghiên cứu TSKH Lê Văn Hoàng viết.
B1 Chuẩn bị File Input:
Để chạy get_wf_xie-modifyO2.f ta cần hai file input: một file chứa thông tin về orbital nguyên tử (*.atom) và một file chứa thông tin về orbital phân tử
(*.molecule)
Hai file có được từ file *.out khi chạy chương trình Gaussian.
File 2: *.molecule. Chọn phần “Molecular Orbital Coefficients”
Chú ý: Khi lưu trong mục save as type ta phải chọn All Files, 19
19 3
B2 Quy định giá trị các tham số
Mở file: dipole.h file này quy định các giá trị của các hằng số. VD: Pi=3.141592654
B3 Chạyget_wf_CHN.f
Trong đó ta kiểm tra các thông tin NAB, NSHELL phù hợp với thông tin trong file *.molecule (các thông tin đó được đánh dấu trong phần trên)
Chú ý các nút trên màn hình.
Số thông tin về các phân tử (trong VD có 19 thông tin về phân tử 1, có 5 thông tin về phân tử 2, có 19 thông tin về phân tử 3)
Số lớp elctron trong phân tử(trong VD phân tử 1 có 5 lớp, phân tử 2 có 3 lớp, phân tử 3 có 5 lớp)
Ta lần lượt nhấn các nút sau:
Complie: để chuyển thông tin sang ngôn ngữ máy.
Nếu không có lỗi: phía dưới màn hình sẽ có thông tin như sau:
Nếu có lỗi: máy sẽ báo như sau và có chỉ dẫn đến vị trí có lỗi, nhờ đó ta có thể sửa lỗi dễ dàng (trong VD: máy báo lỗi ở dòng open1(1,file="coe.MOLECULE")
Built: kiểm tra cấu trúc thông tin. Execute Program: chạy HOMO File xuất có dạng input.wf
Bước 2: Chạy LewMol_2.2.f, đây cũng là một source code viết bằng FORTRAN sử dụng mô hình Lewenstein để tính toán HHG phát ra do phân tử, nguyên tử tương tác với laser cường độ mạnh.
B1 Chuẩn bị File input:
Khi chạy LewMol cần các file sau:
File omega-new.in chứa thông tin về laser chúng ta sử dụng: (có 5 số)
Số thứ nhất là bội số của chu kì, sau n chu kì, ta sẽ có xung là 30fs.
Số thứ hai(x 1014) là cường độ đỉnh. Số thứ ba: tần số omega (eV)
Số thứ tư là pha (pi) Số thứ năm là thế ion hóa.
File Rinfo.inp để kiểm tra xem thông tin có giống trong file *.gif của Guassian không?
NAtom: số nguyên tử của phân tử (trong ví dụ ta thấy có 3 nguyên tử)
Units: đơn vị dùng, nếu đơn vị Bohr ghi số 1, đơn vị Ams ghi số 2 (trong ví dụ ta chọn đơn vị Ams nên ghi số 2)
Dòng tiếp theo là tọa độ của các nguyên tử, được rút từ file *.gif
Dòng tiếp là các đại lượng mà LewMol_2.2.f có thể tính, nếu ta tính thì ghi 1, không tính thì ghi 0.
File Orient.in chứa thông tin về góc hợp bởi vectơ phân cực của laser và trục phân tử.
B2 Chạy LewMol_2.2.f
Sau khi đã chuẩn bị sẵn sàng các file input, làm các bước tương tự như khi chạy file: get_wf_CHN.f cho kết quả trong file hhg.dat, trong file kết quả có 5 cột: cột 1 là bậc của HHG, cột 2 là giá trị bình phương của dipole theo phương song song, cột 3 là giá trị log10 của dipole theo phương song song, cột 4 là giá trị bình phương của dipole theo phương song song, cột 5 là giá trị log10 của dipole theo phương song song.
Tuy nhiên, bậc HHG lẻ, khó sử dụng, ta sử dụng hhg-average-newH.f để tính gần đúng bậc HHG.
Bước 3: Chạy hhg-average-newH.f
Tương tự chạy: hhg-average-newH.f cho kết quả trong file hhg-power.dat.
Kết quả có 5 cột tương tự như trong file hhg.dat: cột 1 là bậc của HHG, cột 2 là giá trị bình phương của dipole theo phương song song, cột 3 là giá trị log10
của dipole theo phương song song, cột 4 là giá trị bình phương của dipole theo phương song song, cột 5 là giá trị log10 của dipole theo phương song song.