GIỚI THIỆU
Vật lý tính toán ngày càng phát triển mạnh mẽ, cống hiến cho nhân loại một cái nhìn toàn diện hơn về vật liệu ở những điều kiện, trạng thái mà thực nghiệm vẫn chưa thực hiện được Tính toán mô phỏng có thể đặt vật liệu vào các điều kiện khắc nghiệt như ở áp suất cực lớn, nhiệt độ cực thấp, với các chất độc hoặc chất phóng xạ Như vậy, tính toán mô phỏng là bước tiên phong trong việc xác định và dự đoán các tính chất của vật liệu, đồng thời giúp định hướng cho các thí nghiệm thực tiễn, tiết kiệm thời gian và chi phí cho thực nghiệm Việc tiên đoán trước các tính chất của vật liệu cũng như, tìm ra vật liệu bền vững có tính chất như mong muốn nằm trong khả năng của tính toán mô phỏng vật liệu Các tính toán mô phỏng cũng giúp giải thích tốt hơn các hiện tượng xảy ra trong thực nghiệm
Trong thực tế, nhiều kỹ thuật mô phỏng khác nhau có thể được sử dụng Liệt kê đơn giản theo thứ tự từ phương pháp cho mô hình lớn đến mô hình ít nguyên tử gồm: phương pháp Phần tử hữu hạn, phương pháp Monte Carlo, phương pháp Động lực học phân tử (MD), phương pháp Lý thuyết phiếm hàm mật độ (DFT) Trong các kỹ thuật mô phỏng, phương pháp DFT là mạnh mẽ nhất vì có thể giải thích các tính chất lượng tử của vật liệu Phương pháp này thích hợp để nghiên cứu các vật liệu có ít nguyên tử, nằm trong kích thước nanomét, như graphene, germanene và trong bài luận văn này là SiC hai chiều (2D-SiC) Mục tiêu cốt lõi của phương pháp DFT là giải phương trỡnh Schrửdinger của cỏc điện tử bằng lý thuyết cơ học lượng tử Với sự phát triển của khoa học máy tính và vật lý lý thuyết, DFT trở thành phương pháp được áp dụng nhiều nhất để đơn giản hóa bài toán nhiều hạt thành phương trình Kohn–Sham theo mật độ điện tử Việc giải phương trình Schordinger dừng của điện tử trở nên dễ dàng hơn và có thể thực hiện với mô hình cỡ dưới 1000 nguyên tử Từ đó, thu được cấu trúc điện tử và các tính chất liên quan như tổng năng lượng, mật độ trạng thái, tính dao động mạng Tất cả kết quả này đều có ý nghĩa quan trọng để xác định và dự đoán các tính chất của vật liệu, ứng dụng cho xã hội
Một trong những ứng dụng sáng giá và giàu tiềm năng là pin nhiên liệu, để thay thế cho nguồn nhiên liệu hóa thạch Pin nhiên liệu sử dụng nhiên liệu hydro, với các đặc tính nổi bật như thân thiện với môi trường, có thể tái tạo, ứng dụng được trên phương tiện di chuyển và trong công nghiệp năng lượng [1-3] Và cũng chính là câu trả lời cho việc thiếu hụt nhiên liệu và ô nhiễm môi trường do nguồn nhiên liệu hóa thạch gây ra Thách thức nêu trên thúc đẩy các nghiên cứu về vật liệu mới, ứng dụng được trong pin nhiên liệu Để chế tạo ra vật liệu hai chiều có tính chất phù hợp với từng loại ứng dụng trong pin nhiên liệu, việc tính toán mô phỏng là vô cùng cần thiết Các tính chất tổng năng lượng, cấu trúc vùng điện tử, dao động mạng quyết định khả năng tồn tại của loại vật liệu 2D-SiC khi bị hydro hóa, có thể được thu được bằng phương pháp Lý thuyết phiếm hàm mật độ
Luận văn này nghiên cứu tương tác của vật liệu SiC hai chiều với hydro bằng phương pháp Lý thuyết phiếm hàm mật độ Nhằm đưa ra các kết luận về khả năng ứng dụng của 2D-SiC trong pin nhiên liệu thông qua các kết quả bao gồm: Vị trí và cấu hình hấp phụ hydro bền vững nhất, sự thay đổi cấu trúc của 2D-SiC khi hấp phụ hydro, năng lượng hấp phụ hydro, năng lượng tương tác giữa nhiều hydro, tính bền nhiệt động qua phổ phonon Ngoài ra, còn có khảo sát cấu trúc điện tử của 2D-SiC trước và sau khi hấp phụ hydro, năng lượng điểm không và các dạng sai hỏng cấu trúc của 2D-SiC có hấp phụ hydro.
TỔNG QUAN
Vật liệu hai chiều
Sự phát triển không ngừng của khoa học và công nghệ trong thời gian gần đây đã giúp con người liên tục có những phát minh và khám phá khoa học quan trọng trong nhiều lĩnh vực khác nhau Một lĩnh vực được nghiên cứu rộng rãi là vật liệu nano, loại vật liệu có kích thước đặc trưng khoảng từ 1nm đến 100nm Mục đích để phát triển các vật liệu mới có hiệu suất làm việc tốt hơn cho nhiều ứng dụng khác nhau Vật liệu nano đã mang lại nhiều ứng dụng trong lĩnh vực y học, vật liệu, điện tử Trong lĩnh vực y sinh, vật liệu nano được sử dụng trong điều trị ung thư, băng bó vết thương, dẫn truyền thuốc…[4] Một cách phổ biến để phân loại vật liệu nano là dựa trên kích thước các chiều cấu trúc, gồm các nhóm: vật liệu không chiều, có ba chiều cấu trúc đều nằm trong kích thước nano (như chấm lượng tử, fullerenes…); vật liệu một chiều có một chiều không gian nằm ngoài kích thước nano (như ống nano carbon
…); vật liệu hai chiều (2D) có hai chiều nằm ngoài kích thước nano (như graphene…); tương tự, vật liệu ba chiều (3D) hay vật liệu khối thông thường, hoàn toàn không nằm trong kích thước nano (xem Hình 2.1) Tùy thuộc vào kích thước và cấu trúc mà vật liệu thể hiện các tính chất hoàn toàn khác nhau
Vật liệu hai chiều, gọi tắt là vật liệu 2D, là loại vật liệu cực mỏng, thường bao gồm một lớp nguyên tử duy nhất hoặc một vài lớp nguyên tử được sắp xếp trong một mặt phẳng hai chiều Những vật liệu này có các tính chất đặc biệt do có diện tích bề mặt lớn dễ dàng tăng tốc độ các phản ứng Việc có kích thước nano, mỏng trên một chiều không gian, cũng giúp loại vật liệu này mang những tính chất điện, quang đặc biệt Ví dụ như là đặc tính có cấu trúc điện tử dễ bị thay đổi khi có chất hấp phụ lên, thích hợp ứng dụng vào các loại cảm biến khí [5] Việc khám phá ra graphene đã mở ra con đường nghiên cứu và phát triển vật liệu hai chiều mới Các tính chất vật lý và hóa học độc đáo của các vật liệu 2D đã thúc đẩy mạnh mẽ các nghiên cứu thí nghiệm, giúp vật liệu 2D ngày càng ưu việt hơn so với các vật liệu khối trong nhiều loại ứng dụng khác nhau Tương tác giữa vật liệu hai chiều và hydro được các nghiên cứu đặc biệt quan tâm chú ý bởi vì các tương tác này ảnh hưởng trực tiếp đến hiệu suất của các ứng dụng, như ứng dụng pin nhiên liệu màng trao đổi proton, trong dự trữ hydro, và trong linh kiện điện tử, mang đến nhiều giá trị thực tiễn
Hình 2.1 Carbon có đủ các thù hình từ (a) không chiều tăng dần đến (d) ba chiều.
Giới thiệu về SiC hai chiều
Trước đây có các ý kiến cho rằng vật liệu hai chiều không tồn tại nếu thiếu một cơ sở 3D Vào năm 2004, graphene dưới dạng lục giác 2D đã được khám phá trong thí nghiệm bóc tách lớp của Geim và Novoselov [6] Kể từ khám phá đó, rất nhiều nỗ lực nghiên cứu đã được dành để tìm hiểu các tính chất của vật liệu hai chiều Trong đó có SiC hai chiều (2D-SiC), là loại vật liệu hai chiều có một lớp nguyên tử duy nhất, gồm các nguyên tử silicon (Si) và carbon (C) được nối chặt trong một mạng lục giác (Hình 2.2) Khó có thể tổng hợp được SiC hai chiều một lớp thuần túy, trong thực tế 2D-SiC được tổng hợp có từ hai đến ba lớp Cấu trúc nguyên tử của SiC dạng khối thường được thấy là dạng lập phương 3C-SiC hoặc dạng lục giác như 4H-SiC,
6H-SiC Đây đều là các chất bán dẫn sử dụng phổ biến trong các ứng dụng điện và quang ở nhiệt độ cao [7] Các mô phỏng động lực học trước đó của nhóm tác giả Tranh, cho SiC hai chiều nung chảy ở nhiệt độ cao,
SiC hai chiều cũng thừa hưởng các đặc tính bền cơ học của SiC và có thêm có thêm các tính chất điện độc đáo do có một chiều không gian cấu trúc nằm trong kích thước nanomét, nơi các hiệu ứng lượng tử ảnh hưởng mạnh mẽ đến tính chất vật liệu Nhiều tính chất của vật liệu SiC cũng được bắt nguồn từ hai nguyên tố Si và C
Hình 2.2 Mô hình tinh thể 2D-SiC, với Si màu xanh và C màu nâu [8] Đầu tiên, silicon là chất bán dẫn kí hiệu là Si và số nguyên tử là 14, nằm trong nhóm 14 trong bảng tuần hoàn hóa học Do đặc tính dễ phản ứng hóa học với ôxi, nờn đến tận năm 1823 nguyờn tố này mới được Jửns Jakob Berzelius tổng hợp dưới dạng tinh khiết Silicon là chất bán dẫn được sử dụng nhiều nhất trong các linh kiện điện tử hiện đại, là loại vật liệu không thể thiếu cho việc chế tạo chip xử lý, transistor, điót, vi mạch Các phức chất sủa Si cũng có tính bền cơ học cao, sử dụng các ứng dụng gốm sứ, thủy tinh (silica - SiO2), chế tạo cao su (silanes - SiH4) và polymer sinh học (silicones) Tiếp đến là carbon, kí hiệu hóa học C, là chất nằm thứ 6 trong trong bảng tuần hoàn hóa học, chung nhóm 14 với nguyên tố Si Carbon là nguyên tố phi kim, mang nhiều thù hình khác nhau như than chì, kim cương, carbon vô định hình Đây là nguyên tố quan trọng nhất cho mọi sự sống, có mặt trong mọi hợp chất hữu cơ quan trọng như protein, đường, chất béo, axit nucleic… Khả năng đặc biệt của carbon là có thể dễ dàng liên kết bền vững với nhiều nguyên tố khác, tạo ra những hợp chất vô cùng phong phú, ứng dụng được trong mọi mặt của đời sống như trong nhiên liệu hóa thạch, sợi carbon, luyện kim Có thể thấy điểm chung lớn của hai loại nguyên tố này, là có tính bền cơ học cao, dễ dàng tạo liên kết với các nguyên tố khác Đều này có thể được giải thích thông qua cấu trúc điện tử vùng hóa trị của hai nguyên tố
Với cấu trúc lớp điện tử ngoài cùng là 2s 2 2p 2 , carbon có bốn điện tử trong vùng hóa trị Khi các nguyên tử carbon liên kết với nhau, năng lượng toàn phần giảm đi do một số vùng điện tử chồng lên nhau Năng lượng này đủ để đưa điện tử ở phân lớp 2s lên 2p Vì vậy các thù hình khối của carbon thường liên kết với nhau bằng liên kết lai hóa sp 3 có hình dạng tứ diện Còn các thù hình hai chiều của carbon (như graphene, ống nanocarbon) thì chỉ gồm các điện tử nằm ở các phân lớp 2s, 2px và 2py, không nằm trên 2pz nên có liên kết lai hóa là sp 2 [1], lí giải cho cấu trúc lục giác tổ ong của graphene Tương tự trong trường hợp của Si cùng nhóm nguyên tố với carbon, có bốn điện tử ở lớp ngoài cùng là 3s 2 3p 2 , thường tạo lai hóa sp 3 trong cấu trúc khối, và cũng được cho là mang lai hóa gần giống sp 2 trong cấu trúc hai chiều là silicene [9] Tuy mang cùng mang hình dạng lục giác tổ ong nhưng silicene lại có độ vênh giữa các nguyên tử, do lai hóa sp 2 không hoàn toàn, vẫn có xác suất điện tử xuất hiện tại phân lớp 3pz
Như vậy, với khả năng tạo ra tối đa 4 liên kết hóa học, với 4 điện tử ở các lớp hóa trị, nguyên tố C và Si có thể dễ dàng liên kết với nhau bằng liên kết sp 3 , tạo ra hợp chất bền vững SiC Các cấu trúc khối của SiC như 3C-SiC, 4H-SiC, 6H-SiC đều rất bền về cơ học lẫn hóa học và thường được thấy trong các ứng dụng điện tử hoặc phát quang ở nhiệt độ cao Sự tồn tại của thù hình hai chiều của C và Si là graphene và silicene, cũng đặt ra khả năng tồn tại loại thù hình hai chiều của SiC với liên kết lai hóa sp 2
Sự tồn tại của 2D-SiC đã được dự đoán bằng các tính toán nguyên lý ban đầu kết hợp với vật lý thống kê ở nhiệt độ cao của nhóm nghiên cứu của Yakobson và Shi vào năm 2015 [10] Cùng năm đó, Shisheng Lin với các cộng sự đã tổng hợp ra những tấm SiC mỏng đầu tiên với độ dày từ 0.3–5 μm, bằng phản ứng hóa học nhiệt độ cao giữa graphene và Si [11] Đến năm 2021, thì 2D-SiC mới thực sự ra đời, nhóm nghiên cứu của Sakineh Chabi đã bóc tách lớp ướt 2D-SiC thành công từ thù hình lục giác 6H-SiC, thu được với tấm 2D-SiC mỏng với độ dày chỉ 0.25 nm Kết quả này gây nhiều ngạc nhiên do liên kết sp 3 trong 6H-SiC là bền hơn liên kết Van der Waals của graphene rất nhiều, gây khó khăn trong việc bóc tách lớp Kết quả của họ cho thấy, 2D-SiC có thể tồn tại bền vững ở điều kiện chuẩn, không bị thay đổi cấu trúc theo thời gian [12] Các nghiên cứu trên cùng chỉ ra rằng, cấu trúc của 2D-SiC là hoàn toàn phẳng, với cả Si và C đều nằm trên cùng một bề mặt, không có độ vênh như silicene và germanene [9, 13] Đều này chứng tỏ giữa Si và C có liên kết lai hóa sp 2 hoàn toàn, tương tự như trong graphene Tuy nhiên, độ dài giữa liên kết Si-C dao động từ 1.77 đến 1.85 Å, lớn hơn liên kết của carbon trong graphene là 1.42 Å, giúp cho mạng lục giác của 2D-SiC cũng lớn hơn cỡ 20 – 30% Diện tích bề mặt lớn của 2D-SiC có thể giúp cải thiện khả năng hấp phụ hydro so với graphene Vì ở graphene, liên kết C-C quá ngắn, không có đủ vị trí cho hydro hấp phụ, khiến tỉ lệ hấp phụ hydro trên graphene chỉ mật độ đạt dưới 1 wt% ở điều kiện chuẩn [15,16] Với liên kết dài hơn, có cơ sở để tin rằng SiC có thể hấp phụ hydro dễ dàng với tỉ lệ hấp phụ là 1:1, một hydro trên mỗi nguyên tử bề mặt Si hoặc C Độ dài liên kết Si-C cũng quyết định nhiều tính chất khác của vật liệu 2D-SiC Bằng mô phỏng nguyên lý ban đầu và động lực học phân tử, nhiều nhóm nghiên cứu chỉ ra các độ dài liên kết Si-C khác nhau như của nhóm tác giả Sun với độ dài liên kết là 1.77 Å [17], nhóm của tác giả Susi thì cho ra độ dài liên kết là 1.89 Å [18], dẫn đến độ rộng vùng cấm cũng khác nhau, thay đổi từ 2.3 đến 2.6 eV [17-19] Điều này, chứng tỏ 2D-SiC là loại vật liệu bán dẫn có vòng cấm lớn Để làm sáng tỏ về độ dài liên kết của 2D-SiC, nhóm nghiên cứu của Sharma đã tiến hành tính toán ứng suất giãn nén, với kết quả cho thấy mô hình 2D-SiC có thể bị giãn nén đến 16% mà vẫn giữ được sự ổn định cấu trúc và hình dáng của vùng điện tử, riêng độ rộng vùng cấm là tăng khi bị nén, và giảm khi bị kéo [19] Chứng minh cho khả năng điều chỉnh độ rộng vùng cấm, tính bền cơ học của 2D-SiC, giúp cho vật liệu này thích hợp cho các ứng dụng điện tử bán dẫn công suất cao
Trong luận văn này, chiều dài liên kết giữa Si-C được thu từ kết quả mô phỏng động lực học trước đó của nhóm, tại nhiệt độ phòng (300K) [8] Điều này giúp DFT thể hiện đúng hơn tính chất của vật liệu ở nhiệt độ phòng, vì các phương pháp giải phương trình sóng đều tính toán ở nhiệt độ 0K Dựa trên cơ sở các nghiên cứu trước đây của nhóm về phản ứng với hydro của vật liệu hai chiều germanene [13], những bước tối ưu hóa mô hình tương tự cũng được kế thừa và mở rộng Sau khi thu được kết quả mô hình 2D-SiC chính xác bằng các kỹ thuật hội tụ năng lượng, các tính toán DFT tiếp tục được thực hiện để tính năng lượng hấp phụ của hydro trên bề mặt 2D-SiC Các xấp xỉ dao động điểm không cũng được áp dụng làm tăng độ chính xác của năng lượng hấp phụ Ngoài ra, nhiều vị trí hấp phụ của hydro được tính toán từ đó chỉ ra vị trí hấp phụ hydro tốt nhất trên 2D-SiC Những biến đổi trong cấu trúc của 2D-SiC khi hấp phụ hydro cũng được ghi nhận Các cấu trúc sai hỏng mạng của 2D-SiC cũng được cho hấp phụ hydro, gồm sai hỏng dãy (8-4), sai hỏng thêm nút, sai hỏng khuyết đôi (4-10-4), sai hỏng khuyết đôi (8-4-4-8), sai hỏng khuyết đôi (4-8-8-4) Cần xác định rõ tương tác giữa hydro và 2D-SiC vì năng lượng hấp phụ hydro quyết định được nhiều tính chất trong pin nhiên liệu như tốc độ phản ứng, chiều phản ứng, mật độ dự trữ hydro.
Ứng dụng của hydro
Hydro được phát hiện rất sớm, vào năm 1766 bởi Henry Cavendish, đây là nguyên tố đầu tiên của bảng tuần hoàn, có số nguyên tử là 1 Đây là nguyên tố nhẹ nhất và phong phú nhất trong Vũ Trụ, vì nó là sản phẩm của các phản ứng phân rã nhiệt hạch Trên bề mặt Trái Đất có khá hiếm hydro ở dạng phân tử tinh khiết, do phân tử này quá nhẹ và dễ dàng thoát khỏi trọng lực của Trái Đất Tuy nhiên, hydro có khả năng dễ tạo ra liên kết hóa học giúp cho chúng vẫn tồn tại nhiều trong các hợp chất như nước (H2O) và hợp chất hữu cơ (CH4 ) Hydro là nguyên liệu lí tưởng, vì có thể tái tạo được, thân thiện với môi trường, cho hiệu suất năng lượng cao Trước khi pin nhiên liệu ra đời, hydro được xếp vào nguồn nhiên liệu không an toàn, do ngọn lửa hydro có nhiệt độ cao, dễ cháy, và dễ lan rộng, không thích hợp cho động cơ đốt trong Khi ứng dụng của pin nhiên liệu được đưa vào hoạt động như xe hydro, nhiên liệu hydro được nhìn nhận lại Khác với động cơ đốt trong có khả năng gây cháy nổ khi rò rỉ nhiên liệu, hydro dùng trong pin nhiên liệu chỉ tỏa nhiệt nhẹ, và không có phản ứng cháy nổ trong lúc hoạt động [3] Các nước trên thế giới đã nhanh chóng mở ra các hướng nghiên cứu về dự trữ hydro, thiết bị sử dụng hydro, vận chuyển hydro, tái tạo hydro… để sử dụng nguồn nhiên liệu này hiệu quả Các cách để tái tạo hydro có thể kể đến như, điện phân nước, pin nhiên liệu methanol, phản ứng CO và nước, phản ứng phân hủy ammonia… Trong đó, dòng điện kích hoạt cho cách phản ứng có thể lấy từ năng lượng tự nhiên, như năng lượng mặt trời hoặc thủy điện, giúp cho hydro là nguồn nhiên liệu bền vững
Tuy nhiên, nhiên liệu hydro cũng có những hạn chế, liên quan đến các tính chất vật lí của phân tử H2 Ở điều kiện tiêu chuẩn, hydro là chất khí, các nguyên tử hydro kết hợp với nhau tạo thành phân tử H2 Trọng lượng phân tử hydro là 2.016 g/mol và có khoảng cách giữa hai nguyên tử hydro thay đổi trong khoảng từ 3.1 đến 2.4 Å [20- 22] Trong thực tế có thể bắt gặp hydro ở dạng khí hoặc dạng lỏng, với điểm ba trạng thái ở -259.1°C/0.07 bar và điểm tới hạn ở -239.8°C/13 bar (Hình 2.3) Ở áp suất tiêu chuẩn 1 bar, nhiệt độ đông đặc của hydro là -259°C, còn nhiệt độ bay hơi là -253°C Nhiệt độ tới hạn của hydro ở cỡ -239.8°C, điều đó có nghĩa ở mọi nhiệt độ cao hơn - 239.8°C, dù áp suất tăng bao nhiêu hydro vẫn là chất khí và cũng sẽ không bao giờ trở thành chất lỏng Có thể nhận xét là ở điều kiện chuẩn, hydro chỉ tồn tại ở thể khí Ngoài ra, nhiệt độ bay hơi và nhiệt độ đông đặc của hyrdo cũng cách nhau không qua nhiều, dẫn đến khả năng hydro ở trạng thái đông đặc hoặc trạng thái lỏng dễ bay hơi rò rỉ khi nhiệt độ biến động Để giữ cho hydro ở trạng thái lỏng, nhiệt độ cần được giữ trong khoảng -240°C đến -259°C với áp suất tương ứng Ngoài ra, dưới áp suất khí quyển, khối lượng riêng của hydro ở nhiệt độ phòng là 0.09 kg/m 3 ở dạng khí, còn ở dạng lỏng là 70.8 kg/m 3 (-253°C), dạng đông đặc là 70.6 kg/m 3 (-262°C) [22]
Như vậy, khối lượng riêng của hydro tăng gần một ngàn lần khi hóa lỏng, giúp hydro có thể được trữ nhiều hơn trong không gian nhỏ hơn Tuy nhiên do nhiệt độ đông đặc và hóa lỏng quá thấp, hydro vẫn thường được dự trữ dưới dạng khí bằng các bình nén áp suất cao
Ngoài ra, một lưu ý khác là hydro rất dễ cháy, chỉ cần một lượng hydro nhỏ là đã có thể xuất hiện ngọn lửa Ngọn lửa hydro màu xanh nhạt, gần như vô hình, khiến cho đám cháy hydro khó nhìn thấy, tiềm ẩn nhiều nguy hiểm Ở áp suất khí quyển, nhiệt độ tự bốc cháy tối thiểu của hydro trong không khí là 570°C Hydro có thể cháy khi chiếm từ 4% đến 75% thể tích không khí Nằm ngoài khoảng này nghĩa là hydro đã đạt giới hạn phát nổ hoặc không đủ nhiên liệu đốt cháy [20] Khả năng dễ cháy của hydro đặt ra những tranh cãi về độ an toàn của nguồn nhiên liệu này, đặc biệt khi hydro thường được dự trữ, vận chuyển trong các bình nén áp suất cao
Như vậy những trở ngại chính của hydro là rất khó dự trữ và vận chuyển Các kỹ thuật dự trữ hydro hiện tại đòi hỏi các yêu cầu an toàn đặc biệt, sử dụng công nghệ, kỹ thuật cao, như bình nén áp suất cao, bình siêu cách nhiệt Cần phải chi trả cho các thiết bị dự trữ, vận chuyển hydro đắt tiền, và thiết kế lại toàn bộ cơ sở hạ tầng, từ trạm phát hydro đến xe chạy bằng nhiên liệu hydro, mới có thể sử dụng nguồn nhiên liệu hydro an toàn hiệu quả Các nghiên cứu đã chỉ ra, vật liệu hai chiều có thể giúp giải quyết các vấn đề liên quan tới dự trữ hydro, và cải tiến hiệu quả của pin nhiên liệu [23] Với khả năng dễ tạo ra liên kết hóa học, hydro có thể dễ dàng được các vật liệu hai chiều như germanene, graphene hấp phụ và cố định [13] Quá trình này cũng làm thay đổi cấu trúc và tính chất điện của vật liệu, giúp chúng thích hợp hơn các ứng dụng quang, và điện
Hình 2.3 Giản đồ pha cho hydro [21].
Giới thiệu về pin nhiên liệu
Pin nhiên liệu là thiết bị biến đổi nhiên liệu hóa học trực tiếp thành năng lượng điện thông qua các phản ứng điện hóa Tuy nhiên không giống như các bình ác quy thông thường, phải sạc điện hoặc thay mới thành phần phản ứng để có thể tiếp tục sử dụng, pin nhiên liệu có thể liên tục tạo ra điện khi nào nguồn nhiên liệu vẫn còn Dựa trên lý thuyết này, nhiều loại pin nhiên liệu khác nhau được phát triển như: pin màng trao đổi proton (Proton exchange membrane - PEM), pin methanol (DMFC), pin nhiên liệu oxít rắn (SOFC), pin nhiên liệu carbon nóng chảy (MCFC), pin nhiên liệu kiềm (AFC), và pin acid phosphoric (PAFC) Mỗi loại đều có ưu điểm riêng, sử dụng nguồn nhiên liệu, lớp xúc tác khác nhau và có những ứng dụng nhất định Dựa vào ứng dụng, các loại pin nhiên liệu được chia ra các loại là: cầm tay, trạm cố định và di động (Bảng 2.1)
Bảng 2.1 Ứng dụng của pin nhiên liệu trong thực tế [2,3]
Trong đó, nổi bật là pin nhiên liệu màng trao đổi proton, được tin tưởng là đã sẵn sàng ứng dụng vào phương tiện di chuyển Với nguồn nhiên liệu sử dụng trong phản ứng điện hóa là hydro và oxi, và sản phẩm thải ra sau cùng là nước Pin PEM không xảy ra phản ứng đốt cháy nhiên liệu với hiệu suất chuyển hóa năng lượng cao hơn phản ứng đốt cháy Quá trình trình hoạt động của pin PEM cũng không tạo ra tiếng ồn, hay thải ra khí thải nhà kính Với hydro là nguyên tố phong phú bậc nhất trong vũ trụ, và hiện tại đã có nhiều phương pháp tái tạo hydro sạch như: điện phân nước bằng năng lượng mặt trời, phân hủy vật chất hữu cơ, chuyển hóa sinh học… Nguồn nhiên liệu hydro sử dụng trong pin nhiên liệu PEM là một nguồn nhiên liệu bền vững, an toàn cho con người, giúp cho môi trường sống tự nhiên không còn bị ô nhiễm bởi khí thải từ việc đốt cháy nhiên liệu hóa thạch
Thiết kế cơ bản của pin nhiên liệu PEM gồm có ba thành phần chính gồm: anode, cathode, và quan trọng nhất là màng điện phân Ở anode, xảy ra phản ứng ôxi hóa hydro khiến, hydro bị mất electron, tạo thành ion dương H + Ion này chỉ gồm proton dẫn đến tên gọi màng dẫn proton
Loại ứng dụng Cầm tay Trạm cố định Di động Định nghĩa
Pin nhiên liệu tích hợp sạc nhiều lần, được thiết kế để có thể cầm tay mang đi được, là thiết bị phát điện hỗ trợ
Pin nhiên liệu cung cấp điện (và có thể cung cấp nhiệt) được đặt cố định, không thể di chuyển được
Pin nhiên liệu cung cấp điện cho các loại phương tiện di chuyển
Công suất 1 W – 20 kW 0.5 kW – 2 MW 1 – 300 kW
Loại công nghệ PEM, DMFC,
Thiết bị GPS, điện thoại, laptop, thiết bị liên lạc quân sự
Trạm phát điện, kết hợp giữa điện và nhiệt, nguồn cung cấp điện liên tục, xe tải, tàu
Xe chạy bằng pin nhiên liệu, xe tải và xe bus, xe nâng, phương tiện đường sắt
Chất điện phân nằm giữa anode và cathode, có tác dụng bẻ gãy liên kết hóa học của hydro, và dẫn truyền proton từ anode đến cathode Ion H + lúc này đi theo chất điện phân đến cathode, còn dòng điện tử được dẫn ra mạch điện bên ngoài để sử dụng Ở cathode thì xảy ra phản ứng khử ôxi thành nước
Phản ứng ở cathode: ẵO2 + 2H + + 2e - → H2O (2.2) Phản ứng tổng: ẵO2 + H2 → H2O (2.3)
Như vậy với sản phẩm là nước và dòng điện, pin nhiên liệu PEM có thể tạo ra dòng điện hoàn toàn thân thiện với môi trường Các ưu điểm của pin nhiên liệu PEM so với các pin nhiên liệu khác bao gồm: Công suất đa dạng, khối lượng nhẹ, tốc độ phản ứng nhanh, dãy nhiệt độ hoạt động rộng (-40 đến 90 o C) Thích hợp cho các ứng dụng công suất trung bình và thấp (như điện thoại, máy tính), đến các hoạt động đòi hỏi nhiều năng lượng (như xe ô tô, trạm phát điện)
Thành phần quan trọng nhất của pin nhiên liệu là màng điện phân, gồm nhiều tế bào pin nhiên liệu ép lại với nhau Mỗi tế bào pin nhiên liệu lại gồm 3 lớp vật liệu: Lớp đế màng, có tác dụng giữ cho tế bào pin nhiên liệu cố định và vững chắc; Lớp dẫn khí, gồm nhiều lỗ thông nhỏ, chủ yếu làm từ polytetrafluoroethylene hoặc các chất kị nước khác, có tác dụng tăng bề mặt tiếp xúc thúc đẩy tốc độ phản ứng, dẫn khí nhiên liệu và đào thải nước; cuối cùng là lớp xúc tác, nơi xảy ra phản ứng ôxi hóa khử, biến hydro và ôxi thành nước và điện, độ dày rất mỏng từ 50-100 μm, với lỗ rỗng cỡ 40–70%, cho phép các hạt bán kính 1–10 nm đi qua dễ dàng [3] Một pin nhiên liệu PEM tốt phải có khả năng dẫn truyền proton nhanh, có thể dẫn khí và đào thải nước tốt, trơ với các chất hóa học, và có cấu trúc bền vững, hoạt động được ở nhiệt độ cao Điểm hạn chế hiện tại của pin nhiên liệu PEM đó là lớp xúc tác của pin nhiên liệu chủ yếu làm từ kim loại hiếm đắt tiền, bạch kim (Pt) Kim loại Pt được chọn bởi vì khả năng bẻ gãy liên kết hydro mà không sản xuất ra quá nhiều nhiệt, cộng với khả năng cho dòng ion H + đi qua dễ dàng, giúp tốc độ phản ứng tăng nhanh đáng kể Vì thuộc nhóm kim loại hiếm, Pt cũng trơ về mặt hóa học và có sức bền cơ học tốt Trong các ứng dụng của pin nhiên liệu hiện nay thì Pt là một nguyên tố không thể thiếu Tuy nhiên nhược điểm hết sức quan trọng của kim loại này là có giá thành cao và quý hiếm Giá thành của vật liệu Pt thường chiếm hơn một nửa chi phí sản xuất một tế bào pin nhiên liệu Khi kim loại Pt được phủ mỏng lên lớp dẫn khí, việc điều khiển cho độ dày Pt là thích hợp cũng gặp nhiều khó khăn, chỉ tận dụng được 20% lớp xúc tác Các nghiên cứu đã cố gắng phát triển các kĩ thuật khác nhau để tìm ra loại vật liệu thay thế, tăng cường mật độ dẫn proton cho Pt, để cải thiện hiệu quả làm việc của lớp màng xúc tác Vật liệu hai chiều cũng được chú ý đến trong lĩnh vực này do tính chất dẫn điện dị hướng, diện tích bề mặt lớn, mỏng nhẹ, và bền Có tiềm năng làm lớp xúc tác độc lập hoặc dùng làm chất nền cho lớp xúc tác là Pt hoặc các kim loại trong nhóm chuyển tiếp như Mo, Fe, Co và Ni [24], giúp tăng khả năng hóa điện, tăng tốc độ tách hydro và dẫn H + , làm ổn định bề mặt Pt
Vật liệu hai chiều như MoS2 đặt trên cấu cúc cầu nano carbon đã được thí nghiệm làm lớp xúc tác và đạt tốc độ phản ứng cao [25] Khả năng hoạt động của một lớp xúc tác được đánh giá qua các chỉ số như sau: thế kích hoạt, gồm năng lượng kích hoạt để xảy ra phản ứng; mật độ dòng, chỉ khả năng dẫn proton của vật liệu; độ ổn định với môi trường và thời gian hoạt động; phổ trở điện hoá, điện trở dẫn càng nhỏ thì phản ứng càng dễ xảy ra; và năng lượng tự do của hydro hấp phụ Phương pháp phiếm hàm mật độ sẽ giúp khảo sát được hai chỉ số là độ ổn định và năng lượng tự do của hydro hấp phụ hệ trên vật liệu, giúp một phần sáng tỏ khả năng làm lớp xúc tác của 2D-SiC trong ứng dụng của pin nhiên liệu.
Các kỹ thuật dự trữ hydro
Hydro ứng dụng trong pin nhiên liệu có thể được dự trữ bằng phương pháp vật lý (nhiệt độ, áp suất) hoặc phương pháp hóa học (phản ứng hóa học, vật liệu hấp phụ) Phương pháp vật lý phù hợp cho việc dự trữ hydro ở trạng thái khí hoặc trạng thái lỏng Phương pháp hóa học thì hydro được dự trữ bằng các phức chức như hợp chất hữu cơ (CH4) hoặc vật liệu hấp thụ/hấp phụ hydro Các phương pháp dự trữ hydro bằng vật liệu hấp phụ được mong đợi là có thể dự trữ nhiều hydro hơn và an toàn hơn để vận chuyển so với cách dự trữ truyền thống Vì có thể hấp phụ hydro dưới dạng nguyên tử, bỏ qua các tính chất nhiệt động lực học khó chịu của phân tử khí H2 Khả năng dự trữ hydro thường được đánh giá bằng phần trăm hydro theo khối lượng (wt%)
Trong đó, 𝑀 𝐻 là khối lượng hydro được dự trữ, chia cho tổng khối lượng của hệ dự trữ 𝑀 𝑆 và hydro Năm 2022, Bộ Năng Lượng Mỹ (DOE) đã đặt ra các mục tiêu cho công nghệ dự trữ hydro, để nhiên liệu này có khả năng được ứng dụng hiệu quả [26] Trong đó, mục tiêu phần trăm hydro theo khối lượng là 6.5 wt% với nhiệt độ hoạt động trong khoảng -40 đến 60 o C, tốc độ nạp xả hydro từ 3-5 phút Trong mô phỏng, DFT cần phải hơn 10 wt% mới thì có khả năng đạt được mục tiêu này ở nhiệt độ phòng Chưa công nghệ dự trữ nào đáp ứng được tất cả các tiêu chí này, tất cả đều có ưu và khuyết điểm
Các phương pháp vật lý có lợi thế là giải phóng hydro nhanh, hoạt động được trong khoảng nhiệt độ cho phép, nhưng lại có mật độ hydro dự trữ thấp Hydro có thể được dự trữ dưới dạng khí nén trong các bình chứa áp suất cao, làm từ thép, nhôm, hoặc các vật liệu composite mới như sợi carbon, polyethylene… giúp các loại bình chứa hydro hiện đại chịu được áp suất cực kì cao lên đến 300 đến 750 bar mà vẫn đảm bảo an toàn [27] Khi tăng áp suất, khối lượng riêng của hydro trong cũng được cải thiện từ 0.30 đến 1.16 kg/L Tuy nhiên, do khối lượng của các bình nén khí rất nặng, được thiết kế chắc chắn, mới đảm bảo được độ bền và ổn định Nên mật độ hydro dự trữ theo khối lượng của phương pháp này không được cao, chỉ đạt mật độ theo khối lượng từ 2 wt% cho bình nén sắt hoặc nhôm, và tối đa 10 wt% cho các bình từ vật liệu composite mới [27]
Ngoài ra, hydro còn có thể được dự trữ dưới trạng thái lỏng So với trạng thái khí, khối lượng riêng của hydro lỏng tăng lên gần một ngàn lần, giúp dự trữ được nhiều hydro hơn trong cùng một thể tích Tuy nhiên, nhiệt độ để hóa lỏng hydro là cực kì thấp, phải ít nhất -239.8 °C, hydro mới chuyển pha sang thể lỏng [22] Một hạn chế khác của phương pháp này là hydro rất dễ hóa hơi, do các biến động nhiệt độ Khiến cho bình chứa hydro lỏng lúc nào cũng có một phần khí hydro, hình thành do sự trao đổi nhiệt với môi trường Lượng hydro bay hơi tùy thuộc vào khả năng cách nhiệt của bình chứa Để dự trữ hydro lỏng trong thời gian dài, cần phải thường xuyên rút bớt khí hydro, để duy trì giảm áp suất trong bình Như vậy phương pháp dự trữ hydro dưới trạng thái lỏng tốn nhiểu chi phí cho năng lượng hóa lỏng hydro, đòi hỏi bình chứa siêu cách nhiệt kỹ thuật cao và phải có sự giám sát vận hành của con người Vì các lý do đó, phương pháp dự trữ vật lý thông dụng nhất vẫn là các bình nén khí áp suất cao
Một kỹ thuật dự trữ hydro khác là bằng phương pháp hóa học Hydro sẽ phản ứng hóa học tạo liên kết, hấp thụ trong mạng tinh thể hoặc hấp phụ lên trên bề mặt vật liệu Phương pháp này đầy hứa hẹn do có mật độ hydro dự trữ cao và dễ dàng vận chuyển, được tin là phù hợp với các ứng dụng phương tiện di chuyển sử dụng nhiên liệu hydro Tuy nhiên chúng có hạn chế là tốc độ nạp xả chậm, đòi hỏi các điều kiện hoạt động khắc nghiệt, nằm ngoài vùng nhiệt độ cho phép Các vật liệu dự trữ hydro thường được nghiên cứu và thí nghiệm là phức chất kim loại nhẹ bị hydride hóa như M(BH4)n, M(AlH4)n, M(NH2)n,… với M là các kim loại nhẹ như Li, Be, Na, Mg, B và Al Mật độ hydro theo khối lượng của các vật liệu này có thể đạt đến tối đa là 18.5 wt% của LiBH4 [28], các vật liệu khác trong nhóm này cũng đạt khả năng dự trữ hydro tốt như NaBH4 (10.7 wt%) [29], LiAlH4 (10.5 wt %) [30] Tuy nhiên, nhiệt độ để giải hấp hydro quá cao, lên đến 268 o C cho LiBH4 hoặc nhỏ nhất là 150 o C cho
NaBH4 và LiAlH4 Khiến các vật liệu này kém bền khi giải hấp nhiều lần, đồng thời không ứng dụng được trên phương tiện di chuyển
Vật liệu nano như graphene, ống nano carbon, cũng được xem xét làm vật liệu dự trữ hydro, vì các vật liệu này có diện tích bề mặt lớn cho hydro dễ dàng hấp phụ Diện tích bề mặt là một yếu tố quan trọng, ảnh hưởng đến khả năng hấp phụ hydro của vật liệu Diện tích lớn giúp hydro dễ tiếp cận bề mặt vật liệu, nhiều phân tử hydro được hấp phụ hơn, tăng tốc độ phản ứng tạo liên kết Ở graphene, đáng ngạc nhiên là vật liệu này có khả năng dự trữ hydro kém Ở nhiệt độ phòng mật độ dự trữ hydro của graphene chỉ cỡ 1 wt% [15,16], nếu giảm nhiệt độ về -196 o C lượng hydro hấp phụ tối đa vẫn chỉ đạt hơn 1 wt% chút đỉnh [16] Tuy nhiên, các nhà nghiên cứu không bỏ cuộc vì vật liệu graphene có tính chất cơ học bền bỉ, có liên kết hóa học mạnh mẽ, có thể dễ dàng uốn cong, tạo ra nhiều dạng thù hình khác nhau Các dạng thù hình này có khả năng dự trữ hydro cao hơn, như thù hình hai lớp graphene bọc lên nhiều nguyên tử hydro được nhóm Slepchenkov [31] mô phỏng động lực học, cho ra mật độ dự trữ đạt 6.86 wt%, với phần lớn hydro vẫn ở trạng thái phân tử, thuận lợi cho phản ứng giải hấp Tính chất cơ học bền bỉ của graphene cũng dẫn tới phương pháp giải hấp bằng cách uốn cong bề mặt graphene bằng điện trường ngoài hoặc tác dụng lực, giúp có thể giải hấp nhanh chóng hydro trên bề mặt graphene bị hydride hóa [32] Khi bị hydride hóa, graphene có thể đạt mật độ dự trữ tối đa đến 8.3 wt% [33] Ngoài ra, hydro còn có thể được dự trữ giữa các lớp graphene, các nghiên cứu đã chỉ ra khoảng cách tối ưu để dự trữ hydro là 0.7 nm, khoảng cách này được duy trì bằng cách kẹp giữa các loại phân tử hữu cơ, cho ra mật độ dự trữ đạt 7.04 wt % [34] Như vậy, bằng sự nỗ lực nghiên cứu, các nhà khoa học đã cải thiện được khả năng dự trữ của hydro, nhờ vào tính bền và dễ tạo liên kết của graphene Hiện tại, vẫn chưa có kết luận nào về khả năng dự trữ hydro của 2D-SiC là loại vật liệu mang nhiều tính chất cơ học thuận lợi giống như graphene, gồm tính bền, nhẹ, diện tích bề mặt lớn, chịu được áp suất và nhiệt độ cao Có thể tin rằng, vật liệu 2D-SiC sẽ còn phát triển hơn nữa trong tương lai và được sử dụng trong nhiều ứng dụng khác nhau, trong đó có khả năng cải tiến công nghệ dự trữ hydro.
Phương pháp lý thuyết phiếm hàm mật độ
Lý thuyết phiếm hàm mật độ (density functional theory - DFT) là một trong những lý thuyết quan trọng nhất hiện nay, được sử dụng phổ biến như một công cụ tính toán nguyên lí ban đầu mạnh mẽ, tìm ra trạng thái lượng tử của nguyên tử, phân tử, vật liệu Lý thuyết phiếm hàm mật độ được Thomas và Fermi cho ra đời vào năm 1927 Lúc này Thomas và Fermi biểu diễn các hệ điện tử bằng hàm phân phối xác suất Tiếp theo đó, dựa trên nền tảng phát triển của cơ học lượng tử, Hohenberg, Kohn và Sham đã phát triển lý thuyết phiếm hàm mật độ vượt trội hơn cả lý thuyết của Thomas-Fermi, mở ra cánh cửa rộng lớn để tính toán các tính chất lượng tử cho các hệ thực tế Trong đó, phương pháp lý thuyết phiếm hàm mật độ cho phép giải gần đỳng phương trỡnh Schrửdinger cho trạng thỏi cơ bản cho hệ nhiều điện tử bằng cỏch biểu diễn các tính chất của hệ N điện tử thông qua hàm mật độ điện tử Vì vậy, phương pháp DFT có thể mô tả chính xác tính chất của hệ dưới mức độ mô tả chính xác của điện tử và hạt nhân, giới hạn trong trạng thái cơ bản Ưu điểm lớn này cũng đủ để DFT được sử dụng trong nhiều nghiên cứu các tính chất hóa học, vật lý của các hệ vật liệu, từ nguyên tử, phân tử cho tới chất rắn Ngày nay, lý thuyết phiếm hàm mật độ đã phát triển vượt bậc và được tích hợp sẵn sàng trong nhiều chương trình tính toán mô phỏng Siesta là một trong những chương trình máy tính, ứng dụng DFT để có thể thực hiện các tính toán nguyên lí ban đầu của hệ vật lý, số lượng hàng trăm nguyên tử
2.6.1 Phương trỡnh Schrửdinger - Bài toỏn hệ nhiều hạt
Mọi vật liệu đều cấu thành nên từ các nguyên tử, phân tử, đi sâu hơn là hạt nhân và điện tử Trọng tõm chớnh của DFT là giải gần đỳng phương trỡnh Schrửdinger dừng với hệ nhiều điện tử và hạt nhân Trong cơ học lượng tử, các điện tử và hạt nhân không chỉ đơn thuần là các hệ riêng lẻ, mà còn có thể tương tác ảnh hưởng đến nhau qua tương tác tĩnh điện, hay lực hấp dẫn, khiến cho việc giải phương trình Schrửdinger của hệ nhiều điện tử trở nờn bất khả thi Hiện tại chưa cú phương phỏp nào giải chớnh xỏc được phương trỡnh Schrửdinger của hệ nhiều điện tử
Phương trỡnh Schrửdinger dừng của một hệ gồm điện tử và hạt nhõn là phương trình đầu tiên nhất của hệ vật liệu có dạng:
Trong đó 𝐻̂ là toán tử Hamilton của hệ, 𝜓 𝑛 là hàm sóng của hệ ứng với trị riêng là mức năng lượng 𝐸 𝑛 Lần lượt 𝑟⃗ 𝑖 và 𝑅⃗⃗ α là véc tơ xác định vị trí của các điện tử và hạt nhân tương ứng Toán tử Hamilton đại diện cho mọi tương tác của các hạt nhân và điện tử với nhau Trong cơ học lượng tử, toán tử Hamilton cấu thành từ các tương tác động năng, và tĩnh điện Toán tử Hamilton được thể hiện như sau
𝐻̂ = 𝑇̂ 𝑒 + 𝑇̂ 𝑛 + 𝑉̂ 𝑒𝑒 + 𝑉̂ 𝑒𝑛 + 𝑉̂ 𝑛𝑛 (2.6) Trong đó, 𝑇̂ 𝑒 𝑇̂ 𝑛 lần lượt là toán tử động năng của eletron và hạt nhân Còn 𝑉̂ 𝑒𝑒 𝑉̂ 𝑒𝑛
𝑉̂ 𝑛𝑛 lần lượt là tương tác điện giữa điện tử với điện tử, điện tử với hạt nhân, và hạt nhân với hạt nhân Các toán tử trên được thể hiện thông qua: vị trí điện tử và hạt nhân (𝑟⃗, 𝑅⃗⃗), cùng với khối lượng điện tử và hạt nhân (m, M) và điện tích của điện tử e như sau
Như vậy trên mặt lý thuyết, để giải được phương trình này cần đầy đủ các biến tọa độ của điện tử và hạt nhân 𝑟⃗ và 𝑅⃗⃗ Khiến cho phương trình này có bộ biến lên đến ba lần số lượng điện tử với hạt nhân cộng lại, và không thể giải chính xác được với các công cụ tính toàn hiện đại
Cỏch đầu tiờn để giảm số biến của phương trỡnh Schrửdinger hệ nhiều hạt, bao gồm hạt nhân và điện tử, là xấp xỉ Born – Oppenheimer Có thể thấy, do proton (+1e) và điện tử (-1e) có độ lớn điện tích giống nhau, dẫn tới tương tác điện lên hạt nhân và điện tử cũng có độ lớn gần giống nhau Tuy nhiên, khối lượng hạt nhân lại lớn hơn điện tử hàng ngàn lần (proton và nơtron đều nặng cỡ 2000 lần điện tử), nên lực tác dụng sẽ làm hạt nhân thay đổi vận tốc chậm hơn điện tử theo tỉ lệ chênh lệch giữa khối lượng điện tử và hạt nhân Như vậy trong một khoảng thời gian rất ngắn, có thể xấp xỉ là hạt nhõn đứng yờn Thớch hợp trong phương trỡnh Schrửdinger dừng, khụng xét hệ số thời gian, bỏ qua động năng của hạt nhân, chỉ tập trung vào điện tử Ý tưởng chính của phương pháp này là tách phương trình sóng điện tử và hạt nhân Dẫn đến việc năng lượng của trạng thái cơ bản không phụ thuộc hạt nhân mà chỉ phụ thuộc vào hàm sóng của điện tử Hàm sóng của điện tử được mô tả như sau
𝑒𝑙(𝑟, 𝑅)𝜓 𝑒𝑙 (𝑟, 𝑅) = 𝐸 𝑒𝑙 (𝑅)𝜓 𝑒𝑙 (𝑟, 𝑅) (2.11) với 𝐻̂ 𝑒𝑙 là toán tử Hamilton của điện tử, 𝜓 𝑒𝑙 (𝑟, 𝑅) là hàm sóng điện tử và 𝐸 𝑒𝑙 là năng lượng điện tử Toán tử Hamilton của điện tử có dạng
Trong đó 𝑇̂ là động năng của điện tử, 𝑉̂ 𝑛𝑒 là thế tương tác tĩnh điện giữa hạt nhân và điện tử, và 𝑉̂ 𝑒𝑒 là thế tương tác tĩnh điện giữa điện tử − điện tử Toán tử Hamilton của điện tử chỉ bao gồm các tương tác của điện tử, giảm được thành phần động năng và tương tỏc tĩnh điện của hạt nhõn với nhau khi so sỏnh với phương trỡnh Schrửdinger
Vì vậy khi tính tổng năng lượng cuối cùng của hệ, cần bổ sung:
Việc giải phương trỡnh Schrửdinger cho ở trạng thỏi cơ bản vật liệu nằm chủ yếu ở việc tính Etot, có nghĩa là tính 𝐸 𝑒𝑙 , với phần 𝐸 𝑛𝑢𝑐 là năng lượng của hạt nhân được bỏ qua và được tớnh riờng như là một hệ số Từ đú việc giải phương trỡnh Schrửdingercho hệ nhiều hạt xoay quanh việc giải phương trỡnh Schrửdinger cho điện tử Việc giảm bớt ảnh hưởng của hạt nhân phần nào giúp bài toán đơn giản hơn, nhưng phức tạp lớn nhất là bộ biến vẫn còn phụ thuộc vào 3N điện tử vẫn chưa giải quyết được
2.6.3 Các định lý Hohenberg-Kohn
Dựa trên cơ sở của mô hình Thomas–Fermi, Hohenberg và Kohn đã thành công trong việc thể hiện toàn bộ năng lượng của hệ thành hàm phụ thuộc mật độ và tổng quát lên cho thấy rằng nguyên tắc này có thể được áp dụng cho mọi hệ điện tử bất kỳ Hàm mật độ điện tử 𝜌 tại vị trí 𝑟⃗ được mô tả như sau
Trong đó N là số lượng điện tử, nhân với tích phân tất cả xác xuất tìm được của điện tử tại vị trí 𝑟⃗ Trong toàn bộ vùng không gian, tích phân của hàm mật độ điện là số lượng điện tử
Hai định lý của Hohenberg và Kohn phát triển vào năm 1964 [35] liên hệ giữa tổng năng lượng và mật độ điện tử, giúp đặt nền tảng cho phương pháp DFT tiếp tục phát triển Định lý Hohenberg-Kohn thứ nhất phát biểu như sau: Trong hệ nhiều hạt tương tác, chỉ có một mật độ hạt duy nhất ứng mỗi trường thế ngoài 𝑉 𝑒𝑥𝑡 Thế ngoài là gộp chung của các tương tác với hạt nhân và các trường thế ngoài tác động lên hệ Ý nghĩa của định luật này là trạng thái cơ bản của hạt xác định hoàn toàn bởi thế ngoài và phiếm hàm mật độ Nếu hai hệ có cùng mật độ điện tử, thì cũng sẽ có cùng một mức năng lượng cơ bản Trong hệ không có trường thế phụ, năng lượng tổng được viết lại như sau
Với: 𝐹 𝐻𝐾 [𝜌] = 𝑇 𝑒 [𝜌] + 𝑉 𝑒𝑒 [𝜌] (2.17) Trong đó, 𝐹 𝐻𝐾 là tổng của động năng điện tử và thế tương tác giữa các điện tử, phụ thuộc hoàn toàn vào số lượng điện tử trong hệ Có nghĩa là các hệ có cùng số lượng điện tử sẽ có giá trị 𝐹 𝐻𝐾 giống nhau Năng lượng tổng của hệ lúc này chỉ còn phụ thuộc vào của mật độ điện tử và số lượng điện tử trong hệ Định lý Hohenberg-Kohn thứ hai phát biểu: Tất cả các thành phần năng lượng của hệ nhiều hạt đều thể hiện được dưới dạng phiếm hàm của mật độ Mức năng lượng cơ bản thấp nhất chỉ khi mật độ hạt ở trạng thái cơ bản
Hai định lớ này giỳp đơn giản đỏng kể bộ biến của phương trỡnh Schrửdinger dừng cho điện tử, từ 3N điện tử thành bộ biến theo mật độ trong không gian ba chiều Đồng thời đề ra cỏch giải đơn giản hơn cho phương trỡnh Schrửdinger, thay vỡ tớnh trực tiếp mức năng lượng cơ bản của hệ, ta có thể lấy cực tiểu của hàm 𝑉 𝑒𝑥𝑡 theo mật độ điện tử, trong trường hợp đã tính được 𝐹 𝐻𝐾 Việc tính ra được 𝐹 𝐻𝐾 là rất khó, khi phiếm hàm này phụ thuộc vào số lượng điện tử trong hệ, các điện tử này có thể va chạm, ảnh hưởng, tương tác với nhau, khiến cho việc tính toán vẫn còn điểm phức tạp chưa thể giải được
2.6.4 Phương trình Kohn-Sham Để áp dụng các phát triển mới trong phương pháp DFT của Hohenberg-Kohn, mục tiêu lúc này là cần phải xác định được năng lượng tương tác của N điện tử, 𝐹 𝐻𝐾 Trước khó khăn này, Kohn-Sham đã đưa ra cách tiếp cận mới đối với một hệ N điện tử, có thể thay thế một hệ tương tác, thành một hệ không tương tác, đồng thời thêm một thế ngoài để tương đương với tổng năng lượng của hệ [35] Thế tương tác giữa các điện tử của Kohn-Sham gồm có ba thành phần:
Trong đó 𝑇 𝑠 , là toàn bộ động năng của hệ điện tử không tương tác , 𝐸 𝑥𝑐 là năng lượng trao đổi tương quan, bổ sung cho tất cả các hiệu ứng lượng tử hoặc tương tác giữa các điện từ còn thiếu trong phương trình Tương tác tĩnh điện cổ điển trong công thức, còn có thể viết dưới dạng mật độ điện tử tính toán theo phương pháp Hartree-Fock như sau:
Khi áp dụng cả lý thuyết của Hohenberg-Kohn và xấp xỉ Kohn-Sham, tổng năng lượng có thể viết lại trở thành như sau:
Chương trình Siesta
Ngày nay, việc giải phương trình Kohn – Sham có thể thực hiện bằng nhiều phần mềm khác nhau, như Quantum Espresso và Vasp Trong đó, có Siesta (Spanish Initiative for Electronic Simulations with Thousands of Atoms) là một chương trình tính toán mô phỏng nguyên lí ban đầu cho vật liệu uy tín và hiệu quả Dựa trên lý thuyết phiếm hàm mật độ Siesta sử dụng phương pháp trường tự hợp để giải phương trình Kohn – Sham theo phiếm hàm mật độ, đồng thời tích hợp tối ưu hóa cấu trúc bằng phương pháp động lực học phân tử, thích hợp cho đa dạng các hệ nguyên tử, phân tử, vật liệu Chương trình cho phép sử dụng giả thế duy trình định mức (norm- coserving), xấp xỉ mật độ định xứ (LDA) hoặc xấp xỉ gradient tổng quát (GGA), sử dụng đa dạng bộ cơ sở obital nguyên tử như DZ, SZP, DZP để đạt được khả năng tính toán hiệu quả Siesta còn là chương trình mã nguồn mở, sử dụng miễn phí trên nền tảng Linux, xuất hiện trong hàng ngàn các nghiên cứu khoa học trong lĩnh vực vật liệu, hóa học, sinh học nên được chọn làm công cụ mô phỏng chính trong luận văn này Ngoài ra, luận văn còn sử dụng các chương trình hỗ trợ gồm: trực quan hóa dữ liệu (Origin, Gnuplot), trực quan hóa vật liệu (Avogadro, Vesta), lập trình (Python, Bash) và các chương trình đi kèm với Siesta (Vibra, Gnubands, Fcbuild).