Tổng quan nghiên cứu

Hiện tượng chuyển pha spin (Spin-Crossover - SCO) được phát hiện lần đầu tiên vào năm 1931 bởi nhà khoa học Cambi cùng các cộng sự khi nghiên cứu các phức chất kim loại chuyển tiếp. Kể từ đó, việc làm chủ khả năng biến đổi cấu trúc điện tử và từ tính của các phân tử kim loại chuyển tiếp đã trở thành tâm điểm của ngành vật lý nhiệt độ thấp và công nghệ lưu trữ dữ liệu mật độ cao. Luận văn thạc sĩ khoa học chuyên ngành Vật lý Nhiệt thực hiện tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội vào năm 2015 đã tập trung giải quyết bài toán mô phỏng và định lượng chi tiết các đặc trưng vật lý của phức chất kim loại chuyển tiếp FeL2, trong đó phối tử HL là 4-hydroxy-N0-((pyridin-2-yl)-methylene)-Benzohydrazide.

Vấn đề nghiên cứu cốt lõi xuất phát từ thực tế rằng các mô hình trường phối tử truyền thống chỉ dừng lại ở mức giải thích định tính, hoàn toàn chưa thể dự đoán chính xác cấu trúc hình học ở trạng thái kích thích cũng như chưa làm rõ nguồn gốc vi mô của tính trễ nhiệt. Hệ phân tử FeL2 sở hữu giá trị trễ nhiệt thực nghiệm rất lớn lên tới khoảng 60 K, mở ra tiềm năng to lớn trong việc chế tạo các cổng logic lượng tử, công tắc chuyển mạch phân tử và các ô nhớ nhị phân có khả năng lưu trữ thông tin ở cấp độ terabit trên mỗi inch vuông.

Mục tiêu cụ thể của luận văn là xác định tường minh cấu trúc hình học, cấu trúc điện tử ở trạng thái spin thấp và spin cao, khảo sát sự phân bố mật độ điện tích, sự biến đổi các thành phần năng lượng trong quá trình chuyển pha, đồng thời đánh giá định lượng tác động của 7 môi trường dung môi khác nhau với hằng số điện môi trải dài từ 1 đến 78,54. Kết quả tính toán lượng tử đã đạt độ tương thích cao với thực nghiệm khi sai số khoảng cách liên kết dưới 0,2%, tạo tiền đề vững chắc cho việc định hướng tổng hợp và ứng dụng vật liệu chuyển pha spin trong thực tế vi điện tử.

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

Khung lý thuyết áp dụng

Nền tảng lý thuyết của nghiên cứu được xây dựng trên sự kết hợp giữa lý thuyết trường tinh thể (Crystal Field Theory) và lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT). Trong cấu trúc phối trí bát diện, 5 quỹ đạo 3d suy biến của ion sắt Fe(II) bị tách thành hai mức năng lượng: mức t2g có năng lượng thấp hơn gồm 3 quỹ đạo dxy, dyz, dzx và mức eg có năng lượng cao hơn gồm 2 quỹ đạo dx2-y2, dz2. Khe năng lượng tách mức Delta đóng vai trò quyết định cấu hình spin khi so sánh tương quan với năng lượng kết cặp điện tử Coulomb P khoảng 1,7 eV.

Nghiên cứu vận dụng chặt chẽ định lý thứ nhất của Hohenberg-Kohn công bố năm 1964, khẳng định năng lượng trạng thái cơ bản là một phiếm hàm duy nhất của mật độ điện tử. Đồng thời, mô hình phương trình tự hợp Kohn-Sham thiết lập năm 1965 được áp dụng để chuyển bài toán hệ nhiều hạt tương tác phức tạp gồm N điện tử thành hệ các phương trình một hạt không tương tác trong thế hiệu dụng. Các khái niệm trung tâm bao gồm:

  • Trạng thái spin thấp (Low Spin - LS): Cấu hình t2g6 eg0 với tổng spin S = 0, xuất hiện khi khe năng lượng Delta đạt xấp xỉ 2,4 eV, vượt trội so với năng lượng kết cặp P = 1,7 eV.
  • Trạng thái spin cao (High Spin - HS): Cấu hình t2g4 eg2 với tổng spin S = 2, hình thành khi khoảng cách liên kết giãn nở làm suy yếu trường tinh thể, khiến khe năng lượng Delta giảm mạnh xuống còn khoảng 0,9 eV.
  • Mật độ biến dạng điện tử: Đại lượng đo lường sự tái phân bố điện tích không gian giữa trạng thái phân tử liên kết và các nguyên tử độc lập cô lập.
  • Năng lượng tương quan trao đổi Exc: Thành phần phi cổ điển mô tả hiệu ứng tự tương tác và tương quan lượng tử giữa các hạt fermion.

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

Nghiên cứu sử dụng phương pháp tính toán hóa học lượng tử phiếm hàm mật độ thông qua phần mềm chuyên dụng DMol3. Đối tượng nghiên cứu là phức chất FeL2 với cấu trúc bát diện chứa tâm Fe(II) liên kết với 4 nguyên tử Nitơ và 2 nguyên tử Oxy từ 2 phối tử càng ba kiểu N2O-.

Hệ cơ sở số kép phân cực DNP (Double Numerical plus Polarization) được lựa chọn với bán kính quỹ đạo cắt cụt là 4,6 Å, đảm bảo triệt tiêu sai số do kích thước cơ sở hữu hạn. Năng lượng tương quan trao đổi được xử lý bằng phiếm hàm xấp xỉ gradient tổng quát PBE. Tương tác giữa các điện tử lớp vỏ hóa trị và điện tử lõi được mô tả bằng thế hiệu dụng Dolg-Wedig-Stoll-Preuss. Phân bố điện tích và mômen từ nguyên tử được tính toán chính xác bằng phương pháp phân tích quần cư Mulliken.

Để mô phỏng môi trường thực tế khi phân tử được bọc trong chất nền hoặc dung môi, mô hình điện môi liên tục kiểu vật dẫn COSMO được áp dụng trên 7 loại dung môi: Benzene (hằng số điện môi 2,284), Chloroform (4,806), Methylene chloride (9,08), Pyridine (12,3), Nitrobenzene (35,7), Dimethyl Sulfoxide (46,7) và Nước (78,54). Tiêu chuẩn hội tụ mật độ điện tử được kiểm soát nghiêm ngặt ở mức 1 nhân 10 mũ trừ 6 đơn vị nguyên tử (a.u.), tiêu chuẩn năng lượng đạt 1 nhân 10 mũ trừ 5 a.u. và độ dịch chuyển nguyên tử tối đa 1 nhân 10 mũ trừ 3 a.u. Toàn bộ quy trình mô phỏng được thực hiện và hoàn thiện trong giai đoạn năm 2015.

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

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

Quá trình mô phỏng lượng tử đã mang lại 4 phát hiện khoa học quan trọng:

Thứ nhất, có sự biến đổi hình học sâu sắc giữa hai trạng thái spin. Ở trạng thái spin thấp, độ dài liên kết tính toán trung bình Fe-L là 1,928 Å, hoàn toàn phù hợp với số liệu nhiễu xạ tia X thực nghiệm là 1,932 Å (độ lệch chỉ 0,004 Å, tương đương khoảng 0,2%). Khi chuyển sang trạng thái spin cao, các liên kết Fe-L giãn nở mạnh từ 5% đến 13%, đưa độ dài liên kết trung bình lên 2,133 Å. Trong đó, liên kết Fe-N3 và Fe-N4 tăng mạnh nhất khoảng 13,3% (từ 1,951 Å lên 2,211 Å), còn liên kết Fe-O tăng khoảng 5,5% (từ 1,971 Å lên 2,080 Å).

Thứ hai, chuyển pha spin kéo theo sự biến đổi đột ngột về từ tính và cấu trúc điện tử. Tại trạng thái HS, mômen từ tập trung áp đảo tại ion Fe với giá trị 3,786 Bohr Magneton (muB), phù hợp với cấu hình tổng spin S = 2. Hai điện tử chuyển từ quỹ đạo t2g lên quỹ đạo eg làm suy giảm khe năng lượng Delta từ 2,4 eV xuống còn 0,9 eV.

Thứ ba, ghi nhận dòng chuyển dịch điện tích khoảng 0,5 e từ ion Fe sang các phối tử. Ở trạng thái LS, điện tích của Fe là dương 0,156 e; khi sang trạng thái HS, giá trị này tăng vọt lên dương 0,700 e (tăng khoảng 4,487 lần). Đồng thời, các ion O1 và N3 trở nên tích điện âm mạnh hơn, lần lượt tăng độ lớn điện tích khoảng 16,8% và 21,8%. Chênh lệch tổng năng lượng giữa trạng thái HS và LS là 1,203 eV, trong đó động năng tăng 1,449 eV, tương tác tĩnh điện Coulomb giảm 0,221 eV và năng lượng tương quan trao đổi chỉ biến thiên âm 0,025 eV.

Thứ tư, ảnh hưởng của môi trường dung môi đạt trạng thái bão hòa nhanh chóng. Khi hằng số điện môi tăng từ 1 (chân không) lên 78,54 (nước), độ dài liên kết Fe-L ở trạng thái LS biến thiên không đáng kể (dưới 0,009 Å), trong khi trạng thái HS chỉ dao động nhẹ dưới 0,02 Å. Mọi đại lượng điện tích và khe HOMO-LUMO đều tiệm cận giá trị bão hòa khi hằng số điện môi vượt quá ngưỡng 35,7.

Thảo luận kết quả

Bản chất vi mô của quá trình chuyển pha spin trong phân tử FeL2 bắt nguồn từ sự cạnh tranh trực tiếp giữa năng lượng tách mức trường tinh thể và năng lượng kết cặp điện tử. Khi có kích thích nhiệt độ hoặc ánh sáng, các liên kết Fe-L bị giãn dài từ 5% đến 13%, làm giảm độ phủ giữa các orbital 3d của Fe và orbital p của phối tử. Sự suy giảm tương tác lai hóa p-d này trực tiếp kéo tụt mức năng lượng của các orbital eg, làm khe Delta co hẹp từ 2,4 eV xuống dưới mức năng lượng ghép đôi 1,7 eV. Do đó, các điện tử ưu tiên phân bố độc thân trên toàn bộ 5 quỹ đạo theo quy tắc Hund, hình thành trạng thái spin cao bền vững.

Dữ liệu tính toán được trực quan hóa sinh động thông qua các bản đồ mật độ biến dạng điện tử với mặt đẳng giá trị 0,1 e/Å3. Trên bản đồ, các đám mây màu vàng (vùng mất điện tử) bao quanh trục liên kết Fe-L ở trạng thái LS được thay thế bằng các đám mây màu xanh (vùng nhận điện tử) ở trạng thái HS, minh chứng rõ ràng cho sự dịch chuyển mật độ điện tử từ tâm kim loại sang hệ phối tử. Đồ thị biểu diễn sự phụ thuộc của năng lượng liên kết và điện tích theo hằng số điện môi cho thấy các đường cong biến thiên dốc trong khoảng hằng số từ 1 đến 12, sau đó hoàn toàn đi ngang khi hằng số điện môi đạt trên 35,7.

So sánh với các công trình thực nghiệm trước đây cho thấy, nghiên cứu không chỉ tái hiện chính xác cấu trúc trạng thái LS mà còn dự đoán thành công cấu trúc không gian của trạng thái HS vốn chưa thể đo đạc bằng thực nghiệm do rào cản phân rã nhiệt. Hiện tượng thay đổi thể tích phân tử kết hợp với dòng chuyển dịch điện tích 0,5 e chính là nguồn gốc tạo nên tính tương tác tập thể (cooperativity) mạnh mẽ trong mạng tinh thể, dẫn đến độ trễ nhiệt ấn tượng 60 K của phân tử FeL2.

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

Dựa trên các kết quả tính toán lượng tử chính xác, các khuyến nghị ứng dụng và phát triển nghiên cứu được đề xuất như sau:

  • Tối ưu hóa cấu trúc phối tử hữu cơ: Các nhóm nghiên cứu hóa tổng hợp cần tiến hành gắn thêm các nhóm thế hút điện tử mạnh (như nhóm nitro hoặc halogen) vào vòng thơm của phối tử benzohydrazide trong giai đoạn 2026-2027 nhằm mục tiêu nâng độ trễ nhiệt Delta T vượt mốc 80 K, gia tăng tính ổn định của trạng thái lưu trữ nhị phân ở nhiệt độ phòng.
  • Lựa chọn vật liệu bọc màng mỏng phù hợp: Các kỹ sư công nghệ vật liệu cần ưu tiên tích hợp phức chất FeL2 vào các chất nền polymer hoặc oxit vô cơ có hằng số điện môi thấp dưới 10 (như teflon hoặc màng silica nano) trong chu kỳ thử nghiệm 12 tháng tới, nhằm giữ cho biến thiên độ dài liên kết Fe-L luôn dưới mức 0,5%, bảo toàn nguyên vẹn đặc tính chuyển mạch phân tử.
  • Thiết lập quy trình kích hoạt chuyển pha bằng quang học: Các phòng thí nghiệm quang học lượng tử nên ứng dụng nguồn laser xung siêu ngắn bước sóng tương ứng mức năng lượng kích thích từ 0,9 eV đến 2,4 eV để điều khiển chuyển mạch spin hai chiều với thời gian trễ dưới 100 picosecond, hoàn thiện mô hình thử nghiệm trước năm 2028.
  • Nâng cấp mô hình mô phỏng phiếm hàm lai ghép: Các chuyên gia tính toán lượng tử cần áp dụng các phiếm hàm thế hệ mới kết hợp hiệu chỉnh phân tán van der Waals DFT-D3 trong các dự án mô phỏng giai đoạn 2026-2028, nhằm giảm độ lệch năng lượng chuyển pha Delta E xuống dưới ngưỡng sai số 0,05 eV.

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

Bản luận văn là tài liệu học thuật và kỹ thuật giá trị cao dành cho 4 nhóm đối tượng:

  • Nghiên cứu sinh và nhà khoa học ngành Khoa học Vật liệu: Khai thác phương pháp luận mô phỏng DFT chuẩn hóa trên các phức chất kim loại chuyển tiếp d6, sử dụng trọn bộ cơ sở dữ liệu về cấu trúc hình học và năng lượng tách mức làm tài liệu đối sánh.
  • Kỹ sư phát triển linh kiện Spintronics và bộ nhớ phân tử: Ứng dụng cơ chế trễ nhiệt 60 K và sự chuyển đổi mômen từ giữa 0 và 3,786 muB để thiết kế kiến trúc bộ nhớ truy cập ngẫu nhiên phân tử (Molecular RAM) với mật độ lưu trữ cực hạn.
  • Giảng viên và học viên cao học chuyên ngành Hóa lý và Vật lý Nhiệt: Sử dụng các phân tích chi tiết về sự tách mức trường bát diện, quy tắc Hund và phân bố mật độ điện tích làm giáo trình trực quan cho các môn học Cơ học lượng tử ứng dụng và Vật lý chất rắn.
  • Chuyên gia phân tích dữ liệu hóa lượng tử: Tiếp cận quy trình thiết lập thông số phần mềm DMol3, tối ưu hóa hệ cơ sở DNP và phương pháp xử lý mô hình dung môi COSMO với 7 môi trường phân cực khác nhau.

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

Hiện tượng chuyển pha spin trong phức chất FeL2 diễn ra như thế nào? Dưới tác động của nhiệt độ hoặc ánh sáng, phức chất FeL2 chuyển đổi linh hoạt giữa trạng thái spin thấp (S = 0, nghịch từ) và trạng thái spin cao (S = 2, thuận từ). Quá trình này đi kèm sự giãn nở độ dài liên kết Fe-L từ 5% đến 13% và làm thay đổi toàn diện từ tính cũng như màu sắc phân tử.

Độ trễ nhiệt 60 K của phân tử FeL2 có ý nghĩa thực tiễn gì? Độ trễ nhiệt 60 K mang lại tính lưỡng ổn định (bistability) vượt trội trong một khoảng nhiệt độ rộng. Đặc tính này cho phép hệ phân tử lưu giữ trạng thái thông tin nhị phân bit 0 và bit 1 một cách ổn định mà không cần liên tục tiêu tốn năng lượng để duy trì kích thích bên ngoài.

Năng lượng tương tác tĩnh điện thay đổi ra sao trong quá trình chuyển pha? Khi chuyển từ spin thấp sang spin cao, năng lượng tương tác tĩnh điện Coulomb biến thiên một lượng âm 0,221 eV. Sự thay đổi này được giải thích bởi khoảng cách liên kết kim loại - phối tử giãn dài và sự dịch chuyển của dòng điện tích 0,5 e từ ion Fe sang 6 nguyên tử bao quanh.

Môi trường dung môi có phá vỡ cấu trúc của phân tử FeL2 không? Khảo sát trên 7 dung môi có hằng số điện môi từ 1 đến 78,54 chứng minh cấu trúc phân tử FeL2 cực kỳ bền vững. Trạng thái spin thấp gần như không bị ảnh hưởng, trong khi trạng thái spin cao chỉ biến thiên nhẹ và nhanh chóng đạt giá trị bão hòa ổn định khi hằng số điện môi vượt mức 35,7.

Độ tin cậy của phương pháp tính toán DFT trong luận văn được kiểm chứng ra sao? Phương pháp DFT sử dụng phần mềm DMol3 và bộ cơ sở DNP cho kết quả độ dài liên kết trung bình ở trạng thái spin thấp là 1,928 Å, chỉ sai lệch 0,004 Å (khoảng 0,2%) so với kết quả nhiễu xạ tia X thực nghiệm 1,932 Å, khẳng định độ chính xác và độ tin cậy tuyệt đối.

Kết luận

  • Xác lập thành công cấu hình không gian chính xác của phân tử FeL2 ở cả hai trạng thái spin, chỉ ra sự giãn nở liên kết trung bình từ 1,928 Å lên 2,133 Å (tăng từ 5% đến 13%) khi chuyển pha.
  • Làm sáng tỏ cơ chế chuyển pha spin từ trạng thái nghịch từ S = 0 sang thuận từ S = 2 với mômen từ đạt 3,786 muB do 2 điện tử nhảy từ mức t2g lên mức eg.
  • Định lượng chính xác dòng chuyển dịch điện tích 0,5 e từ tâm Fe sang phối tử và xác định biến thiên tổng năng lượng chuyển pha đạt 1,203 eV.
  • Chứng minh tính bền vững cấu trúc của phức chất trong 7 loại dung môi điện môi và phát hiện ngưỡng bão hòa môi trường tại giá trị điện môi 35,7.
  • Khẳng định giá trị thực tiễn vượt trội của phân tử FeL2 với độ trễ nhiệt 60 K trong định hướng chế tạo linh kiện nhớ phân tử và thiết bị chuyển mạch quang từ thế hệ mới.

Đóng góp lớn nhất của luận văn là đã giải mã toàn diện nguồn gốc vi mô của hiện tượng chuyển pha spin trên hệ FeL2 bằng công cụ lượng tử tiên tiến, lấp đầy khoảng trống dữ liệu thực nghiệm về trạng thái spin cao. Các nhóm nghiên cứu vi mạch và vật liệu nano được khuyến khích liên kết thử nghiệm màng mỏng FeL2 trong giai đoạn 2026-2028 nhằm hiện thực hóa các linh kiện điện tử phân tử thương mại hóa.