Tổng quan nghiên cứu

Nhu cầu chuyển dịch năng lượng toàn cầu đang thúc đẩy mạnh mẽ các công nghệ phát điện sạch, trong đó pin nhiên liệu màng trao đổi proton (PEMFC) nổi lên như một giải pháp đột phá. Trong khi động cơ đốt trong truyền thống chỉ đạt hiệu suất nhiệt khoảng 25%, hệ thống PEMFC có thể đạt hiệu suất chuyển đổi điện năng trực tiếp từ 40% đến 60%, thậm chí vượt mức 85% khi tận dụng nhiệt lượng đồng phát. Tuy nhiên, rào cản lớn nhất đối với việc thương mại hóa đại trà nằm ở giá thành và độ bền của màng điện phân. Vật liệu tiêu chuẩn hiện nay là Nafion có chi phí sản xuất đắt đỏ và bị suy giảm độ dẫn proton nghiêm trọng khi nhiệt độ vận hành vượt quá 80°C.

Để giải quyết thách thức này, nghiên cứu tập trung vào phát triển và đánh giá màng polymer lai hóa trên nền poly(ethylene-co-tetrafluoroethylene) (ETFE) biến tính qua hai giai đoạn: ghép mạch monomer styrene bằng phương pháp chiếu xạ bức xạ và sunfo hóa để gắn nhóm chức dẫn proton -SO3H. Mục tiêu cốt lõi của luận văn là làm sáng tỏ bản chất cấu trúc vi mô, phát hiện các liên kết hóa học đặc trưng, các mode dao động phân tử và quy luật dịch chuyển số sóng khi thay đổi tỷ lệ ghép mạch cũng như tỷ lệ sunfo hóa.

Phạm vi nghiên cứu được thực hiện tại Hà Nội và Thành phố Hồ Chí Minh trong giai đoạn năm 2018, kết hợp giữa đo đạc quang phổ thực nghiệm và mô phỏng hóa tính toán lượng tử. Công trình mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp hệ cơ sở dữ liệu dao động quang học chính xác, góp phần định hướng chế tạo màng dẫn proton thay thế Nafion với độ bền cơ học cao hơn 1,5 lần và hạ giá thành sản xuất pin nhiên liệu xuống dưới ngưỡng 45 USD/kW.

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 dựa trên hai trụ cột lý thuyết vật lý vững chắc: lý thuyết tán xạ quang phổ Raman và lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT).

Hiệu ứng tán xạ Raman dựa trên sự tương tác không đàn hồi giữa photon của chùm bức xạ kích thích đơn sắc và các mode dao động nội tại của phân tử. Dưới góc độ cơ học lượng tử, khi phân tử hấp thụ năng lượng photon kích thích $h\nu_0$, nó chuyển lên trạng thái năng lượng ảo không ổn định rồi nhanh chóng hồi phục về mức dao động $v=1$ hoặc $v=0$, phát ra photon tán xạ Stokes ($h\nu_0 - h\nu_m$) hoặc Anti-Stokes ($h\nu_0 + h\nu_m$). Điều kiện cần để một mode dao động thể hiện hoạt tính Raman là đạo hàm của tenxơ độ phân cực theo tọa độ chuẩn phải khác không ($\partial\alpha/\partial q \neq 0$).

Bên cạnh đó, lý thuyết phiếm hàm mật độ DFT với gần đúng Born-Oppenheimer và hệ phương trình tự hợp Kohn-Sham đóng vai trò then chốt trong việc mô phỏng cấu trúc điện tử. Phương pháp này thay thế việc giải hàm sóng đa hạt phức tạp bằng việc biểu diễn năng lượng trạng thái cơ bản thông qua mật độ electron $\rho(r)$, cho phép tối ưu hóa cấu trúc hình học và tính toán chính xác tần số dao động phân tử cũng như hoạt tính Raman của từng nhóm chức.

Các khái niệm trọng tâm bao gồm: màng dẫn proton PEMFC, polymer nền ETFE, phản ứng ghép mạch bức xạ (radiation grafting), quá trình sunfo hóa (sulfonation) tạo nhóm axit sulfonic (-SO3H) và tenxơ phân cực Raman.

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

Nguồn dữ liệu của luận văn được thu thập từ quy trình thực nghiệm đồng bộ kết hợp tính toán số học trên máy tính:

  1. Cỡ mẫu và quy cách chọn mẫu: Nghiên cứu khảo sát 7 trạng thái mẫu đại diện, bao gồm màng ETFE nguyên thủy, 3 mẫu ETFE ghép mạch styrene với các mức độ ghép (grafting degree - GD) lần lượt là 15,8%, 21,3% và 35,9%, cùng 3 mẫu ETFE sunfo hóa với các mức độ sunfo hóa (degree of sulfonation - DS) là 15,9%, 22,1% và 38,0%. Phương pháp chọn mẫu định lượng hệ thống theo thang tỷ lệ tăng dần giúp đánh giá chính xác mối tương quan tuyến tính giữa thành phần hóa học và cường độ tín hiệu quang học.
  2. Kỹ thuật phân tích phổ thực nghiệm: Sử dụng hệ phổ kế FT-Raman hiện đại trang bị nguồn laser kích thích bước sóng 632,81 nm, dải quét số sóng trải rộng từ 0 đến 4000 cm⁻¹. Việc sử dụng kỹ thuật biến đổi Fourier (FT-Raman) với giao thoa kế Michelson và đầu dò CCD độ nhạy cao giúp triệt tiêu triệt để phông huỳnh quang ký sinh phát sinh từ tạp chất polymer, đem lại độ phân giải không gian quang học đạt mức 1 µm. Lý do lựa chọn phổ Raman là nhờ ưu thế phân tích không phá hủy mẫu, độ nhạy cao với các liên kết đối xứng và không bị ảnh hưởng bởi hơi nước $H_2O$ hay $CO_2$ trong môi trường.
  3. Mô hình hóa và tính toán mô phỏng: Thực hiện trên phần mềm tính toán lượng tử DMol3, khảo sát các chuỗi oligomer ETFE nguyên thủy với số mắt xích lặp lại $n = 2, 3, 4$, ETFE ghép mạch styrene và ETFE sunfo hóa gắn axit benzensulfonic. Tần số dao động và cường độ tán xạ Raman Stokes được tính toán thông qua dẫn xuất sai phân hữu hạn của tenxơ phân cực dưới tác dụng của gradient điện trường ngoài.
  4. Thời gian thực hiện: Toàn bộ quy trình tổng hợp, đo đạc quang phổ và xử lý dữ liệu mô phỏng được hoàn thành trong chu kỳ 12 tháng của năm 2018.

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

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

Quá trình phân tích thực nghiệm và mô phỏng lý thuyết đã làm sáng tỏ cấu trúc vi mô của màng ETFE qua từng giai đoạn biến tính:

  1. Đặc trưng dao động của màng ETFE nguyên thủy: Trên phổ Raman thực nghiệm trong vùng 0 - 4000 cm⁻¹, màng ETFE chưa biến tính thể hiện rõ các dải dao động đặc trưng của khung aliphatic fluoropolymer, nổi bật là dao động biến dạng của nhóm $-CF_2-$ tại vùng số sóng 835 cm⁻¹, dao động kéo dãn $C-C$ tại 1060 cm⁻¹ và dao động kéo dãn đối xứng của nhóm $-CH_2-$ tại vùng 2975 - 3000 cm⁻¹.
  2. Sự hình thành liên kết sau ghép mạch styrene: Khi chiếu xạ ghép mạch styrene vào ETFE với các tỷ lệ 15,8%, 21,3% và 35,9%, phổ Raman xuất hiện các đỉnh phổ hoàn toàn mới minh chứng cho sự hiện diện của vòng thơm benzene. Cụ thể, đỉnh dao động kéo dãn liên kết $C=C$ vòng thơm tại 1602 cm⁻¹ và dao động "ring breathing" tại 1001 cm⁻¹ tăng mạnh về cường độ. Khi tỷ lệ ghép mạch tăng từ 15,8% lên 35,9%, cường độ tích phân của vạch 1001 cm⁻¹ tăng hơn 220%, phản ánh mật độ chuỗi nhánh polystyrene phát triển tỷ lệ thuận với mức độ ghép.
  3. Dấu vết hóa học của quá trình sunfo hóa: Trên các mẫu ETFE sunfo hóa với tỷ lệ 15,9%, 22,1% và 38,0%, xuất hiện rõ nét mode dao động kéo dãn đối xứng của nhóm $-SO_3^-$ tại vị trí số sóng 1128 cm⁻¹ và dao động biến dạng $S-O$ tại vùng 620 - 640 cm⁻¹. Tỷ số cường độ giữa đỉnh 1128 cm⁻¹ và đỉnh nền ETFE tăng đều đặn theo tỷ lệ sunfo hóa, chứng minh các nhóm axit sulfonic đã được gắn thành công vào vị trí para của vòng thơm styrene.
  4. Mức độ tương thích giữa mô phỏng DFT và thực nghiệm: Kết quả tính toán mô phỏng từ phần mềm DMol3 cho các chuỗi cấu trúc $n = 2, 3, 4$ thể hiện sự trùng khớp rất cao với phổ đo thực tế. Độ lệch số sóng giữa giá trị tính toán lý thuyết và giá trị đo bằng laser 632,81 nm kiểm soát trong biên độ dưới 15 cm⁻¹, khẳng định độ tin cậy tuyệt đối của mô hình cấu trúc phân tử được đề xuất.

Thảo luận kết quả

Cơ chế dịch chuyển số sóng và thay đổi cường độ phổ phản ánh sâu sắc tương tác liên phân tử bên trong màng polymer. Khi mạch styrene được ghép vào khung ETFE chính, hiện tượng cộng hưởng electron pi ($\pi$) trong vòng thơm làm biến đổi mật độ điện tích cục bộ, dẫn đến sự dịch chuyển nhẹ khoảng 3 - 5 cm⁻¹ ở các mode dao động của khung $C-C$. Khi tiếp tục sunfo hóa, tính hút electron mạnh của nhóm $-SO_3H$ gây ra hiệu ứng phân cực hóa mạnh mẽ, làm biến dạng trường thế cục bộ và tăng cường độ tán xạ Raman của các liên kết lân cận.

Dữ liệu thực nghiệm và mô phỏng có thể được trực quan hóa tối ưu thông qua biểu đồ xếp chồng đa phổ (stacked Raman spectra) từ 0 đến 4000 cm⁻¹ nhằm so sánh trực tiếp ba trạng thái màng, kết hợp cùng bảng thống kê gán mode dao động chi tiết đối chiếu giữa tần số thực nghiệm và tần số DFT. Khi so sánh với các nghiên cứu phổ FT-IR trước đây của tác giả Lê Y Khoa hay nghiên cứu phổ hủy positron của tác giả Trần Thị Thùy Trang, quang phổ Raman cho thấy tính ưu việt vượt trội trong việc phân tách các mode dao động đối xứng của vòng thơm và liên kết lưu huỳnh mà phổ hồng ngoại thường bị nhòe do hấp thụ nước. Sự kết hợp giữa thực nghiệm và mô phỏng DFT đã lấp đầy khoảng trống nghiên cứu lý thuyết, giúp xác nhận chính xác cấu trúc vi mô ở cấp độ phân tử của vật liệu màng lai hóa.

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

Dựa trên các phát hiện quang phổ và mô hình cấu trúc phân tử, luận văn đề xuất 4 nhóm giải pháp công nghệ nhằm thúc đẩy ứng dụng màng PEM ETFE:

  1. Chuẩn hóa tỷ lệ ghép mạch và sunfo hóa tối ưu: Các viện nghiên cứu vật liệu cần điều chỉnh quy trình chế tạo để kiểm soát tỷ lệ ghép mạch styrene trong dải 20% đến 25% và tỷ lệ sunfo hóa đạt mức 22% đến 25%. Mục tiêu nhằm đạt độ dẫn proton trên 0,08 S/cm tại nhiệt độ 80°C trong khi vẫn duy trì độ bền kéo cơ học trên 25 MPa, triển khai hoàn thiện trong vòng 6 tháng tới.
  2. Ứng dụng giám sát chất lượng trực tuyến bằng quang phổ FT-Raman: Doanh nghiệp sản xuất pin nhiên liệu nên tích hợp module cảm biến quang phổ Raman laser bước sóng 785 nm hoặc 632,81 nm vào dây chuyền tổng hợp màng liên tục. Giải pháp này cho phép kiểm tra không phá hủy độ đồng đều của nhóm $-SO_3H$ trên diện tích màng với thời gian phản hồi dưới 30 giây/mẫu, áp dụng trong giai đoạn 2024 - 2025.
  3. Mở rộng mô phỏng động lực học phân tử và DFT cấu trúc mở rộng: Nhóm nghiên cứu tính toán lý thuyết tại các trường đại học cần mở rộng mô hình DFT từ oligomer ngắn ($n=4$) lên các hệ siêu phân tử có mặt các cụm phân tử nước ($H_2O$) từ 4 đến 12 phân tử để khảo sát chi tiết cơ chế truyền proton dạng Grotthuss và Vehicular, hoàn thành trong lộ trình 12 tháng.
  4. Thử nghiệm độ bền vận hành thực tế trong cell pin đơn (Single Cell): Trung tâm nghiên cứu năng lượng cần tiến hành lắp ráp hệ màng - điện cực (MEA) sử dụng màng ETFE đã sunfo hóa vào cell pin đo kiểm với công suất 1 kW, kiểm tra tuổi thọ liên tục qua 5000 giờ làm việc để đánh giá tốc độ thoái hóa hóa học của liên kết $C-S$, thực hiện trước năm 2026.

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

Tài liệu học thuật này mang lại giá trị thực tiễn và chiều sâu lý thuyết cho 4 nhóm đối tượng chuyên môn:

  1. Học viên cao học, nghiên cứu sinh ngành Quang học và Vật lý chất rắn: Tiếp cận phương pháp luận mẫu mực trong việc kết hợp giữa kỹ thuật thực nghiệm quang phổ FT-Raman và công cụ mô phỏng lượng tử phiếm hàm mật độ DFT/Kohn-Sham để giải quyết bài toán cấu trúc vật liệu.
  2. Kỹ sư và chuyên gia R&D pin nhiên liệu (Fuel Cell Engineers): Nắm vững quy trình tổng hợp, đặc tính hóa học và cơ chế suy thoái của màng trao đổi proton gốc hydrocacbon-fluoropolymer nhằm nghiên cứu vật liệu mới thay thế màng thương mại đắt tiền Nafion.
  3. Giảng viên và nhà nghiên cứu tại các trường đại học khoa học tự nhiên: Sử dụng tài liệu làm bài giảng chuyên đề, tài liệu tham khảo giảng dạy trong các môn học Quang phổ ứng dụng, Hóa lý polyme và Khoa học vật liệu tiên tiến.
  4. Chuyên viên phân tích quang phổ tại các phòng kiểm định chất lượng: Khai thác thư viện số sóng và các mode dao động đặc trưng của hệ polymer ghép mạch để xây dựng quy trình phân tích định tính, định lượng vật liệu polyme biến tính.

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

1. Tại sao kỹ thuật quang phổ Raman lại được ưu tiên lựa chọn để khảo sát màng ETFE thay vì phổ hồng ngoại FT-IR? Phổ Raman có độ nhạy vượt trội đối với các dao động làm thay đổi độ phân cực phân tử như liên kết đối xứng $C=C$ vòng thơm và khung $C-C$, đồng thời tín hiệu tán xạ của nước tự do rất yếu. Do đó, kỹ thuật này cho phép phân tích chính xác cấu trúc màng trong điều kiện ẩm mà không bị che khuất tín hiệu bởi dải hấp thụ rộng của nước như trong phổ FT-IR.

2. Màng ETFE biến tính có những ưu điểm nào vượt trội hơn so với màng Nafion thương mại? Màng ETFE ghép mạch sunfo hóa sở hữu bộ khung copolymer fluorinated có độ bền cơ học cao hơn khoảng 1,5 lần so với màng Nafion 212, duy trì độ ổn định cấu trúc ở nhiệt độ trên 80°C. Đặc biệt, chi phí nguyên liệu thô của ETFE thấp hơn đáng kể, giúp giảm tổng chi phí chế tạo hệ màng - điện cực MEA tới 40%.

3. Vai trò của việc tính toán mô phỏng DFT bằng phần mềm DMol3 trong luận văn là gì? Mô phỏng DFT giúp xác lập cấu hình hình học tối ưu năng lượng của các mắt xích oligomer ETFE từ $n=2$ đến $n=4$. Phần mềm tính toán chính xác tần số dao động lý thuyết và tenxơ phân cực, giúp gán chính xác từng đỉnh phổ thực nghiệm với mode dao động cơ bản tương ứng của các nhóm chức styrene và axit sulfonic.

4. Sự biến đổi của phổ Raman khi tỷ lệ ghép mạch styrene tăng từ 15,8% lên 35,9% diễn ra như thế nào? Khi mức độ ghép mạch tăng từ 15,8% lên 35,9%, mật độ nhóm phenyl trong cấu trúc màng tăng mạnh, thể hiện qua sự gia tăng tuyến tính cường độ vạch 1001 cm⁻¹ (dao động hít thở vòng) và 1602 cm⁻¹ (dao động kéo dãn $C=C$). Các đỉnh dao động của khung ETFE ban đầu giảm tỷ trọng cường độ tương đối do bị che phủ bởi các nhánh polystyrene.

5. Bước sóng laser 632,81 nm mang lại lợi thế gì trong phép đo phổ Raman của vật liệu polymer? Laser He-Ne bước sóng 632,81 nm cung cấp năng lượng kích thích vừa đủ trong vùng ánh sáng đỏ khả kiến, giảm thiểu tối đa hiện tượng phát huỳnh quang nền thường gặp ở bước sóng ngắn (532 nm hay 488 nm). Đồng thời, bước sóng này đảm bảo tín hiệu tán xạ Raman có cường độ đủ lớn để detector CCD ghi nhận với tỷ số tín hiệu trên nhiễu (SNR) cao.

Kết luận

  • Luận văn đã xác định toàn diện cấu trúc dao động quang học của màng ETFE qua 3 trạng thái: nguyên thủy, ghép mạch styrene (15,8% - 35,9%) và sunfo hóa (15,9% - 38,0%).
  • Thiết lập thành công mối tương quan định lượng giữa tỷ lệ biến tính hóa học với cường độ các vạch phổ Raman đặc trưng tại 1001 cm⁻¹ (vòng thơm) và 1128 cm⁻¹ (nhóm $-SO_3^-$).
  • Ứng dụng xuất sắc lý thuyết phiếm hàm mật độ DFT với mô hình chuỗi $n=2, 3, 4$ trên phần mềm DMol3, đạt độ sai lệch tần số dao động dưới 15 cm⁻¹ so với thực nghiệm đo bằng laser 632,81 nm.
  • Lộ trình tiếp theo giai đoạn 2024 - 2026 cần tập trung thử nghiệm độ bền MEA trên hệ pin đơn 1 kW qua 5000 giờ và tối ưu hóa quy trình ghép mạch bức xạ quy mô công nghiệp.
  • Mời quý độc giả, nhà nghiên cứu và doanh nghiệp quan tâm liên hệ, tải toàn văn luận văn và hợp tác phát triển công nghệ vật liệu màng pin nhiên liệu thế hệ mới.