Tổng quan nghiên cứu

Trong khoa học vật liệu nano hiện đại, các cụm nguyên tử (cluster) kích thước siêu nhỏ gồm từ 2 đến khoảng 100 nguyên tử kim loại chuyển tiếp kết hợp bán dẫn đang là tâm điểm chú ý nhờ những đặc tính quang điện từ vượt trội. Đặc biệt, germanium (Ge) với cấu hình electron hóa trị 4s2 4p2 và vùng cấm hẹp đang được định hướng ứng dụng nhằm thay thế silicon trong các linh kiện bán dẫn thế hệ mới. Tuy nhiên, các cluster germanium tinh khiết kém bền nhiệt và có hoạt tính hóa học cao, đòi hỏi phải pha tạp các kim loại chuyển tiếp như sắt (Fe) để cải thiện độ bền và tạo ra các tính chất từ đặc biệt. Vấn đề nghiên cứu cốt lõi nảy sinh khi các trạng thái spin của cluster Fe-Ge có mật độ năng lượng gần như suy biến do sự đóng góp đồng thời của các orbital 3d, 4s của Fe và 4p của Ge, khiến các phương pháp tính đơn cấu hình truyền thống như lý thuyết phiếm hàm mật độ (DFT) gặp sai số đáng kể và không thể mô tả các trạng thái kích thích.

Nhằm giải quyết thách thức này, luận văn thạc sĩ Hóa học của tác giả Nguyễn Thị Lan Phương, thực hiện tại Trường Đại học Đồng Tháp vào năm 2019, tập trung mô tả chi tiết hình học, cấu hình electron, tần số dao động và xác định chính xác trạng thái spin cơ bản cũng như kích thích của cluster FeGen−/0 với n = 1 và n = 2. Nghiên cứu áp dụng các phương pháp hóa học lượng tử đa cấu hình tiên tiến với không gian hoạt động mở rộng lên tới 17 orbital. Kết quả nghiên cứu đã cung cấp năng lượng tách electron chuẩn xác với sai số dưới 0,05 eV, làm cơ sở dữ liệu lý thuyết vững chắc định hướng cho việc tổng hợp vật liệu nano từ tính và linh kiện spintronics trong thực tế.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng cơ học lượng tử phân tử, khởi đầu từ phương trình sóng dừng Schrödinger kết hợp với phép gần đúng Born-Oppenheimer nhằm tách biệt chuyển động của hạt nhân và electron. Để giải quyết tương tác đa điện tử, luận văn kết hợp phương pháp trường tự hợp Hartree-Fock với lý thuyết phiếm hàm mật độ (DFT) dựa trên hai định lý nền tảng của Hohenberg-Kohn và hệ phương trình Kohn-Sham. Các phiếm hàm trao đổi tương quan gradient suy rộng (GGA) như BP86 cùng các phiếm hàm lai ghép hiện đại như B3LYP, TPSS và TPSSh được áp dụng để tối ưu hóa cấu trúc không gian ba chiều.

Đặc biệt, để khắc phục nhược điểm đơn cấu hình của DFT khi khảo sát kim loại chuyển tiếp, khung lý thuyết đa cấu hình CASSCF (Complete Active Space Self-Consistent Field) và CASPT2 (Second-Order Perturbation Theory) được thiết lập nhằm tính toán đồng thời cả năng lượng tương quan tĩnh và tương quan động. Hơn nữa, phương pháp không gian hoạt động giới hạn RASSCF/RASPT2 với sự phân chia 3 phân vùng không gian RAS1, RAS2 và RAS3 được tích hợp, cho phép mở rộng không gian tính toán từ 14 lên đến 40 orbital. Các khái niệm chính bao gồm: trạng thái spin (doublet, triplet, quartet, quintet, sextet), năng lượng tách electron đoạn nhiệt ADE (Adiabatic Detachment Energy), năng lượng tách electron thẳng đứng VDE, tính chất đa cấu hình và hiệu ứng vỏ kép 4d của nguyên tử sắt.

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

Cỡ mẫu nghiên cứu bao gồm 4 hệ phân tử chính: FeGe-, FeGe, FeGe2- và FeGe2. Phương pháp chọn mẫu tập trung khảo sát có chủ đích toàn diện hơn 20 trạng thái spin khác nhau trên 2 dạng đồng phân cấu trúc không gian đặc trưng là dạng hình đường thẳng (nhóm đối xứng C∞v) và dạng hình vòng tam giác (nhóm đối xứng C2v). Cỡ mẫu này mang tính đại diện cốt lõi cho các cluster kích thước nhỏ n = 1-2, giúp làm sáng tỏ bản chất liên kết Fe-Ge trước khi phát triển lên các hệ cụm nguyên tử lớn hơn.

Nguồn dữ liệu tính toán lượng tử được thu thập trực tiếp thông qua việc thực thi các thuật toán hóa toán trên gói phần mềm OpenMolcas và ORCA phiên bản 4.0, kết hợp phần mềm Gabedit để trực quan hóa hình ảnh orbital và kiểu dao động. Tối ưu hóa cấu hình hình học và tính tần số dao động điều hòa được thực hiện bằng phiếm hàm BP86 với bộ hàm cơ sở def2-TZVPP cho cả Fe và Ge. Năng lượng tương đối và bề mặt thế năng được tính toán chuẩn xác bằng CASPT2 và RASPT2 với bộ hàm tương đối tính aug-cc-pwCVTZ-DK cho Fe và aug-cc-pVTZ-DK cho Ge, có tích hợp toán tử hiệu chỉnh tương đối tính Douglas-Kroll bậc hai. Lý do lựa chọn kết hợp RASPT2 và DFT là vì phương pháp này vừa đảm bảo độ chính xác nhiệt động tuyệt đối nhờ xử lý đa cấu hình, vừa tối ưu hóa được chi phí thời gian tính toán trên hệ thống máy chủ hiệu năng cao. Toàn bộ quy trình tính toán và thẩm định được hoàn thành trong tháng 07 năm 2019.

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

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

Thứ nhất, đối với hệ cluster FeGe-/0, nghiên cứu xác định trạng thái cơ bản của anion FeGe- là trạng thái doublet 2A2 với độ dài liên kết Fe-Ge là 2,232 Å và tần số dao động điều hòa đạt 280 cm−1. Trạng thái kích thích đầu tiên là quartet 4B1 có năng lượng tương đối cao hơn trạng thái cơ bản chỉ 0,13 eV. Với cụm FeGe trung hòa, trạng thái cơ bản chuyển thành triplet 3A2 với độ dài liên kết co ngắn lại còn 2,180 Å (giảm 2,33% so với dạng anion) và tần số dao động tăng lên 341 cm−1 (tăng 21,79%). Năng lượng tách electron đoạn nhiệt ADE của bước chuyển từ 2A2 sang 3A2 đạt 0,96 eV theo phương pháp chuẩn CASPT2.

Thứ hai, phân tích hàm sóng lượng tử phát hiện trạng thái cơ bản 2A2 của FeGe- mang bản chất đa cấu hình cực mạnh, trong đó cấu hình electron chính chiếm tỷ trọng chỉ 35% trong tổng hàm sóng. Các trạng thái kích thích khác như 4A1, 2B1 và 6A2 có tỷ trọng cấu hình chính dao động từ 22% đến 35%. Sự xen phủ giữa orbital 3dxz, 3dyz của Fe và 4px, 4py của Ge đóng vai trò chủ đạo tạo nên liên kết bền vững 6π.

Thứ ba, đối với cluster FeGe2-/0, đồng phân hình vòng (đối xứng C2v) thể hiện độ bền vượt trội so với đồng phân hình đường thẳng (đối xứng C∞v), với mức chênh lệch năng lượng nhiệt động lên tới 1,63 eV ở dạng anion và 1,73 eV ở dạng trung hòa điện. Trạng thái cơ bản của FeGe2- là quartet 4B1 (độ dài Fe-Ge là 2,285 Å, Ge-Ge là 2,463 Å), gần như suy biến với trạng thái kích thích 4B2 khi mức chênh lệch năng lượng chỉ là 0,04 eV. Đối với FeGe2 trung hòa, trạng thái cơ bản là quintet 5B2 (độ dài Fe-Ge là 2,235 Å, Ge-Ge là 2,284 Å).

Thứ tư, giá trị ADE của bước chuyển electron ở cluster FeGe2- hình vòng được tính toán bằng RASPT2 đạt 1,68 eV. Khi so sánh với các phiếm hàm DFT, phiếm hàm BP86 cho kết quả sát nhất với 1,56 eV (sai lệch 7,14%), trong khi phiếm hàm TPSS đạt 1,45 eV (sai lệch 13,69%), TPSSh đạt 1,40 eV (sai lệch 16,67%) và B3LYP đạt 1,27 eV (sai lệch cao nhất lên tới 24,40%).

Thảo luận kết quả

Nguyên nhân giúp đồng phân hình vòng FeGe2 đạt độ bền nhiệt động vượt trội hơn dạng đường thẳng bắt nguồn từ hiệu ứng tạo vòng ba cạnh kín, tối ưu hóa việc phân tán mật độ điện tích từ nguyên tử sắt sang hai nguyên tử germanium xung quanh. Việc đưa 5 orbital 4d của Fe vào vùng không gian hoạt động RAS3 trong phương pháp RASPT2 đã mô tả trọn vẹn hiệu ứng vỏ kép, giúp khắc phục triệt để sai số năng lượng 0,22 eV từng xuất hiện trong phương pháp CASPT2 không gian hẹp.

Toàn bộ dữ liệu nghiên cứu có thể được hệ thống hóa trực quan thông qua biểu đồ đường cong thế năng hai chiều mô tả sự biến thiên năng lượng theo khoảng cách liên kết Fe-Ge từ 1,8 Å đến 3,2 Å, kết hợp với bảng ma trận so sánh 8 trạng thái spin của FeGe- và 9 trạng thái spin của FeGe2. So sánh với các công bố quốc tế về cluster kim loại chuyển tiếp cùng phân nhóm như VGen- và CoGen-, quy luật biến đổi cấu hình spin và khoảng cách liên kết của FeGen- hoàn toàn tương thích và làm rõ cơ chế chuyển mức quang electron ở bước sóng 226 nm trong các thí nghiệm thực tế.

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

Thứ nhất, mở rộng quy mô tính toán lượng tử đa cấu hình RASPT2 cho các cụm cluster FeGen kích thước lớn hơn từ n = 3 đến n = 12. Mục tiêu kỹ thuật cần hướng đến là tối ưu hóa không gian hoạt động với hơn 20 orbital nhằm khống chế sai số năng lượng tách electron ADE dưới 0,03 eV. Lộ trình thực hiện khuyến nghị trong vòng 12 đến 18 tháng, do các nhóm nghiên cứu Hóa lý thuyết và Hóa tính toán tại các trường đại học chủ trì thực hiện.

Thứ hai, tích hợp hiệu ứng tương tác spin-quỹ đạo (Spin-Orbit Coupling) và hiệu chỉnh tương đối tính mở rộng vào việc mô phỏng phổ quang electron (PES). Mục tiêu đạt độ trùng khớp đường phổ lý thuyết với thực nghiệm trên 95% đối với các trạng thái spin cao như quintet và sextet. Thời gian triển khai dự kiến từ 6 đến 12 tháng, do các chuyên gia vật lý lượng tử và hóa tin học đảm nhận.

Thứ ba, ứng dụng các tham số cấu trúc hình học và mật độ spin của cụm FeGe2 hình vòng vào quy trình thiết kế vật liệu màng mỏng nano từ tính và các linh kiện bộ nhớ skyrmion. Mục tiêu nhằm gia tăng mật độ lưu trữ dữ liệu lên gấp 2 đến 3 lần và cải thiện độ ổn định từ tính ở dải nhiệt độ phòng 300 K. Thời gian nghiên cứu ứng dụng từ 24 đến 36 tháng, phối hợp giữa các viện nghiên cứu vật liệu và các doanh nghiệp công nghệ bán dẫn.

Thứ tư, thiết lập cơ sở dữ liệu số hóa mở về bề mặt thế năng và tần số dao động của các hệ cluster kim loại chuyển tiếp pha tạp nhóm IVA (Si, Ge, Sn). Mục tiêu số hóa hơn 100 cấu trúc cluster với hơn 500 trạng thái electron kích thích trong thời gian 12 tháng, phục vụ việc tra cứu và đối chuẩn cho cộng đồng khoa học quốc tế.

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

Nhóm thứ nhất là học viên cao học và nghiên cứu sinh chuyên ngành Hóa lý thuyết, Hóa lý và Vật lý nguyên tử. Tài liệu mang lại quy trình chuẩn mực về cách thiết lập không gian hoạt động cho phương pháp CASSCF và RASSCF, giúp người học áp dụng trực tiếp vào việc giải quyết các bài toán hệ đa điện tử phức tạp chứa kim loại chuyển tiếp có trên 14 electron hóa trị.

Nhóm thứ hai là các giảng viên và nhà nghiên cứu trong lĩnh vực hóa học tính toán. Luận văn cung cấp bộ dữ liệu đối chuẩn chính xác cao về năng lượng ADE, độ dài liên kết và tần số dao động, hỗ trợ đánh giá và tinh chỉnh độ tin cậy của các phiếm hàm DFT mới như BP86, TPSS hay TPSSh.

Nhóm thứ ba là các kỹ sư và chuyên gia phát triển vật liệu nano, vi điện tử và spintronics. Việc thấu hiểu cơ chế tương tác electron và tính chất từ tính của cụm Fe-Ge giúp định hình thiết kế các cấu trúc chấm lượng tử và cảm biến nano thế hệ mới thay thế silicon.

Nhóm thứ tư là sinh viên đại học ngành Hóa học, Vật lý kỹ thuật và Khoa học vật liệu. Luận văn là tài liệu tham khảo chi tiết về việc sử dụng các phần mềm ORCA phiên bản 4.0, OpenMolcas và Gabedit, hỗ trợ đắc lực trong quá trình thực hiện khóa luận tốt nghiệp đạt thang điểm xuất sắc từ 9,0 trở lên.

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

Tại sao cluster FeGe- lại có trạng thái cơ bản là doublet 2A2 thay vì quartet 4B1?

Phương pháp tính toán CASPT2 với 14 orbital chứng minh trạng thái doublet 2A2 bền hơn trạng thái quartet 4B1 khoảng 0,13 eV. Nguyên nhân là do cấu hình electron đóng kín 14σ2 và 6π4 kết hợp với 3 electron trên orbital 2δ tạo nên sự xen phủ cực đại giữa orbital 3d của Fe và 4p của Ge, mang lại độ bền nhiệt động lực học cao nhất cho phân tử.

Vì sao phương pháp CASSCF truyền thống bị hạn chế khi tính toán cluster FeGe2-/0?

Phương pháp CASSCF truyền thống bị giới hạn ở không gian hoạt động khoảng 14 electron trên 14 orbital do sự bùng nổ tổ hợp cấu hình vượt quá bộ nhớ máy tính tính toán. Với FeGe2 cần tới 17 orbital hóa trị, việc dùng CASPT2 không gian hẹp sẽ bỏ qua orbital 4d của Fe, dẫn đến sai số năng lượng lên tới 0,22 eV. Phương pháp RASSCF/RASPT2 giải quyết triệt để vấn đề này.

Cấu trúc hình học ảnh hưởng thế nào đến độ bền của cluster FeGe2?

Đồng phân hình vòng đối xứng C2v có độ dài liên kết Fe-Ge là 2,235 Å và Ge-Ge là 2,284 Å, bền hơn đồng phân thẳng đối xứng C∞v từ 1,63 eV đến 1,73 eV. Cấu trúc tam giác ba cạnh này cho phép nguyên tử sắt trung tâm tạo liên kết đồng thời với cả hai nguyên tử germanium, tối ưu hóa việc phân bố điện tích và giải tỏa sức căng hình học.

Phiếm hàm DFT nào cho kết quả tính toán năng lượng tách electron FeGe2- tốt nhất?

Phiếm hàm GGA gradient suy rộng BP86 cho kết quả năng lượng tách electron ADE đạt 1,56 eV, chỉ sai lệch 7,14% so với giá trị chuẩn 1,68 eV của phương pháp RASPT2. Trái lại, các phiếm hàm lai ghép như B3LYP có sai lệch lên tới 24,40% (đạt 1,27 eV), cho thấy phiếm hàm BP86 là lựa chọn tối ưu nhất trong tính toán DFT cho hệ cluster Fe-Ge.

Giá trị năng lượng tách electron ADE có ý nghĩa thực tiễn như thế nào trong phân tích quang phổ?

Năng lượng tách electron đoạn nhiệt ADE (đạt 0,96 eV ở FeGe- và 1,68 eV ở FeGe2-) xác định ngưỡng năng lượng tối thiểu để ion hóa anion thành phân tử trung hòa ở trạng thái bền nhất. Giá trị này là cơ sở trực tiếp để các nhà thực nghiệm đối chiếu, giải thích vạch phổ quang electron quang phổ tử ngoại đo tại bước sóng 226 nm và 266 nm.

Kết luận

  • Xác định chính xác trạng thái cơ bản của FeGe- là doublet 2A2 (độ dài liên kết 2,232 Å) và của FeGe là triplet 3A2 (độ dài liên kết 2,180 Å, tần số dao động 341 cm−1).
  • Khẳng định đồng phân hình vòng của FeGe2- (trạng thái 4B1) và FeGe2 (trạng thái 5B2) bền hơn đồng phân hình đường thẳng từ 1,63 eV đến 1,73 eV.
  • Chứng minh tính ưu việt vượt trội của phương pháp đa cấu hình RASPT2 với 17 orbital so với lý thuyết phiếm hàm mật độ đơn cấu hình trong việc mô tả các trạng thái kích thích gần suy biến.
  • Cung cấp bộ dữ liệu năng lượng tách electron đoạn nhiệt ADE chuẩn mực đạt 0,96 eV cho FeGe- và 1,68 eV cho FeGe2-.
  • Thiết lập kế hoạch chi tiết trong 12 đến 24 tháng tới nhằm phát triển mô phỏng lượng tử cho các cluster FeGen kích thước lớn từ n = 3 đến n = 12.

Luận văn thạc sĩ của tác giả Nguyễn Thị Lan Phương là một công trình nghiên cứu công phu, giải quyết trọn vẹn bài toán cấu trúc electron và trạng thái spin phức tạp của cluster sắt - germanium, đóng góp nguồn dữ liệu lý thuyết nền tảng cho ngành hóa học lượng tử Việt Nam. Quý độc giả, giảng viên và các nhà nghiên cứu hãy tham khảo toàn văn tài liệu để ứng dụng hiệu quả các bộ tham số chuẩn xác này vào các dự án phát triển vật liệu nano tiên tiến.