Tổng quan nghiên cứu

Ngành công nghiệp hóa dầu toàn cầu ghi nhận sự phát triển vượt bậc khi sản lượng etylen thế giới đạt mức 85 triệu tấn mỗi năm vào năm 1999 và tiếp tục tăng trưởng lên khoảng 115 triệu tấn mỗi năm vào năm 2005. Etylen là nguồn nguyên liệu cơ bản quan trọng nhất để sản xuất nhựa polyetylen (PE), trong đó polyetylen tỷ trọng thấp (LDPE) giữ vị trí then chốt trong đời sống và sản xuất công nghiệp. Quá trình tổng hợp LDPE theo cơ chế trùng hợp gốc tự do đòi hỏi các điều kiện kỹ thuật vô cùng khắc nghiệt với dải áp suất cao từ 1000 đến 3000 atm và nhiệt độ dao động trong khoảng 130 đến 320 độ C.

Trong môi trường phản ứng khắc nghiệt này, các gốc tự do trung gian có thời gian tồn tại cực ngắn, thường nhỏ hơn 10 mũ trừ 6 giây, khiến việc khảo sát cơ chế phản ứng và động học bằng phương pháp thực nghiệm gặp rất nhiều rào cản về kỹ thuật, chi phí và mức độ an toàn. Mục tiêu cụ thể của nghiên cứu là ứng dụng phương pháp hóa học tính toán hiện đại để xác định mô hình phản ứng tối ưu, giải thích bản chất động học phản ứng phát triển mạch, làm sáng tỏ cơ chế truyền mạch của propanal so với các hydrocarbon truyền thống, đồng thời lượng hóa mức độ phân nhánh mạch polymer thông qua phản ứng cắn lưng (back-biting).

Nghiên cứu được thực hiện tại Trường Đại học Bách khoa thuộc Đại học Quốc gia Thành phố Hồ Chí Minh với sự hợp tác đào tạo và hỗ trợ học thuật từ Đại học Leuven (Vương quốc Bỉ), hoàn thành vào tháng 12 năm 2003. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc tối ưu hóa định mức tiêu hao nguyên liệu từ 1005 đến 1010 kg etylen trên 1 tấn LDPE thành phẩm, hỗ trợ các nhà máy điều khiển chính xác khối lượng phân tử trung bình và cấu trúc vi mô của polymer.

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 nền tảng cơ học lượng tử và hóa học tính toán hiện đại nhằm khảo sát bề mặt thế năng và trạng thái chuyển tiếp của các phản ứng gốc tự do. Khung lý thuyết trọng tâm bao gồm Thuyết phiếm hàm mật độ (DFT - Density Functional Theory) thông qua việc giải hệ phương trình Kohn-Sham phi tuyến, kết hợp với Thuyết trạng thái chuyển tiếp (TST - Transition State Theory) để xác định năng lượng hoạt hóa và hằng số tốc độ phản ứng. Bên cạnh đó, luận văn xem xét đối chiếu với lý thuyết Hartree-Fock (HF), phương pháp tương tác cấu hình (CI), phương pháp cụm liên kết (Couple Cluster - CCSD) và lý thuyết nhiễu loạn Moller-Plesset bậc hai (MP2) nhằm đánh giá toàn diện năng lượng tương quan điện tử.

Mô hình nghiên cứu tập trung vào ba hệ thống cốt lõi: mô hình phát triển chuỗi giữa gốc tự do etyl và monome etylen, mô hình truyền mạch giữa gốc mạch đang phát triển với các phân tử điều hòa mạch (propanal, propan, propen, buten-1), và mô hình chuỗi alkyl mạch thẳng chứa từ 2 đến 8 nguyên tử cacbon để khảo sát phản ứng chuyển gốc nội phân tử. Các khái niệm chính được vận dụng xuyên suốt bao gồm: Năng lượng hoạt hóa (Ea), Hằng số tốc độ phát triển mạch (Kp), Hằng số tốc độ truyền mạch (Ktr), Tỷ số truyền mạch (Ktr/Kp), và Độ trùng hợp trung bình (DP).

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

Nghiên cứu sử dụng tập dữ liệu tính toán lượng tử toàn diện được thiết lập trên chuỗi 8 gốc tự do mạch thẳng từ etyl đến 1-octyl, kết hợp khảo sát 4 tác nhân truyền mạch phổ biến trên dải nhiệt độ phản ứng từ 130 đến 320 độ C (tương đương 403 đến 593 Kelvin) và dải áp suất cao lên đến 3500 atm. Phương pháp chọn mẫu cấu hình phân tử dựa trên nguyên lý tìm kiếm cấu trúc hình học có năng lượng cực tiểu và xác định chính xác cấu trúc trạng thái chuyển tiếp có duy nhất một tần số dao động ảo.

Công cụ tính toán chủ đạo là các gói phần mềm lượng tử chuyên dụng Gaussian 94 và Gaussian 98. Phương pháp tính toán được lựa chọn tối ưu là phiếm hàm lai hóa B3LYP kết hợp với bộ hàm cơ sở phân cực 6-31G* (6-31G(d)). Lý do lựa chọn B3LYP/6-31G* là vì phương pháp này mang lại sự cân bằng vượt trội giữa độ chính xác hóa học và thời gian xử lý dữ liệu, có khả năng mô tả chính xác hiệu ứng tương quan điện tử cho các gốc tự do vỏ mở (open-shell) với độ sai lệch năng lượng so với thực nghiệm chỉ khoảng 1 đến 2 kJ/mol, khắc phục triệt để các hạn chế của phương pháp Hartree-Fock truyền thống. Toàn bộ quá trình nghiên cứu, mô phỏng và kiểm chứng dữ liệu được triển khai liên tục trong khung thời gian đào tạo cao học từ năm 1999 đến cuối năm 2003.

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

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

Nghiên cứu đã mang lại 4 phát hiện quan trọng có giá trị khoa học và ứng dụng công nghệ cao:

Thứ nhất, phản ứng phát triển mạch tạo LDPE giữa gốc tự do etyl và monome etylen diễn ra theo cơ chế gốc tự do với hàng rào năng lượng hoạt hóa là 28 kJ/mol. Hằng số tốc độ phản ứng phát triển mạch (Kp) tuân theo quy luật hàm mũ nhiệt độ Arrhenius và đạt giá trị 3,16 nhân 10 mũ 7 cm3/mol.s tại nhiệt độ 300 độ C.

Thứ hai, cơ chế truyền mạch của propanal diễn ra ưu tiên thông qua phản ứng tách nguyên tử hydro ở nhóm chức carbonyl (C=O). Tỷ lệ sản phẩm tách hydro tại vị trí này chiếm tới 85% ở nhiệt độ 220 độ C, hoàn toàn tương thích và trùng khớp với các số liệu thử nghiệm thực tế do tập đoàn DSM công bố.

Thứ ba, hiệu quả truyền mạch của propanal vượt trội hoàn toàn so với các hợp chất hydrocarbon truyền thống. Tại nhiệt độ 300 độ C, hiệu quả truyền mạch của propanal đạt tới 90,5%, trong khi propan chỉ đạt 1,8%, propen đạt 34,9% và buten-1 đạt 49,4%. Như vậy, khả năng điều hòa mạch của propanal cao gấp khoảng 50 lần so với propan và gần gấp 2 lần so với buten-1.

Thứ tư, phản ứng cắn lưng (back-biting) gây phân nhánh chuỗi polymer xảy ra mạnh mẽ nhất tại vị trí nguyên tử cacbon số 4 và số 5 tính từ đầu mang gốc tự do. Quá trình này tạo ra các nhánh ngắn phổ biến là nhánh butyl và amyl (pentyl) với mật độ phân nhánh đo được từ 4 đến 5 nhánh trên 1000 đơn vị cacbon (hoặc nhóm CH2) ở điều kiện nhiệt độ 300 độ C.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp propanal đạt hiệu quả truyền mạch vượt bậc là do liên kết C-H gắn trực tiếp với nhóm carbonyl có năng lượng liên kết thấp hơn đáng kể so với các liên kết C-H trong alkan hoặc liên kết allylic trong olefin. Điều này làm giảm đáng kể rào cản năng lượng hoạt hóa của trạng thái chuyển tiếp, giúp tốc độ tách hydro diễn ra nhanh hơn hàng chục lần so với propan và propen.

Khi phân tích phản ứng cắn lưng, cấu hình trạng thái chuyển tiếp dạng vòng 6 cạnh để tạo nhánh butyl sở hữu sức căng vòng thấp nhất và độ bền nhiệt động học cao nhất so với các vòng 4, 5 hoặc 7 cạnh. Sự hình thành các nhánh ngắn butyl và pentyl đóng vai trò quyết định làm giảm độ kết tinh của nhựa, mang lại tính chất mềm dẻo đặc thù cho polyetylen tỷ trọng thấp.

Trong báo cáo kỹ thuật, dữ liệu động học có thể được trực quan hóa sinh động thông qua biểu đồ biểu diễn sự phụ thuộc của tỷ số hằng số tốc độ truyền mạch trên hằng số phát triển mạch (Ktr/Kp) theo thang nhiệt độ từ 130 đến 320 độ C. Đồng thời, bảng số liệu so sánh năng lượng hoạt hóa và phân bố phần trăm nhánh trên 1000 đơn vị CH2 sẽ làm nổi bật sự vượt trội của propanal cũng như xác nhận xu hướng phân nhánh ở từng vị trí cacbon từ C-alpha đến C-octyl.

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

Thứ nhất, đề xuất thay thế từng phần các tác nhân truyền mạch olefin truyền thống (như propen, buten-1) bằng hợp chất propanal trong các dây chuyền sản xuất LDPE cao áp. Mục tiêu cần đạt là duy trì hiệu quả truyền mạch đạt trên 90% ở nhiệt độ vận hành 280 đến 300 độ C, giúp tiết kiệm khoảng 20% lượng hóa chất điều hòa phân tử lượng cần nạp vào hệ thống. Lộ trình triển khai thử nghiệm tại quy mô phân xưởng dự kiến diễn ra từ Quý 1 năm 2024 đến Quý 4 năm 2024, do Phòng Công nghệ và Đội ngũ Kỹ sư Vận hành nhà máy hóa dầu chủ trì thực hiện.

Thứ hai, tối ưu hóa chế độ nhiệt độ trong thiết bị phản ứng dạng ống dao động trong khoảng kiểm soát 180 đến 240 độ C nhằm hạn chế phản ứng cắn lưng, duy trì mật độ phân nhánh ngắn ở mức dưới 3 nhánh trên 1000 đơn vị cacbon. Giải pháp này nhằm nâng cao độ bền kéo đứt của màng nhựa thành phẩm thêm ít nhất 15%. Kế hoạch thực hiện kéo dài 6 tháng trong giai đoạn 2024 - 2025, do Bộ phận Kiểm soát Chất lượng (QA/QC) và Kỹ thuật viên nhà máy phối hợp đảm nhiệm.

Thứ ba, chuẩn hóa và áp dụng quy trình mô phỏng hóa học lượng tử B3LYP/6-31G* trước khi mở rộng nghiên cứu sang các hệ comonome phân cực như vinyl axetat hoặc butyl acrylat. Việc tính toán trước bề mặt thế năng giúp cắt giảm 40% chi phí thử nghiệm thực nghiệm và rút ngắn 50% thời gian nghiên cứu phát triển sản phẩm mới. Thời gian áp dụng bắt đầu từ Quý 3 năm 2024, thuộc trách nhiệm của các Viện nghiên cứu vật liệu cao phân tử và Phòng R&D doanh nghiệp.

Thứ tư, nâng cấp hệ thống bơm định lượng chất khơi mào peroxyt đa điểm và van xả áp định kỳ trong thiết bị phản ứng dạng ống để kiểm soát độ chuyển hóa monome đạt mức tối đa 38% và giữ định mức tiêu hao nguyên liệu ổn định ở mức 1005 kg etylen cho mỗi tấn LDPE. Thời gian hoàn thiện nâng cấp thiết bị trong giai đoạn 2024 - 2026, do Ban Giám đốc Kỹ thuật và Nhà máy Nhựa Rạng Đông phối hợp triển khai.

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

Thứ nhất, các kỹ sư công nghệ hóa học và chuyên gia vận hành nhà máy polymer: Luận văn cung cấp các thông số động học chính xác về hằng số tốc độ Kp, Ktr trên dải nhiệt độ 130 đến 320 độ C và áp suất lên tới 3000 atm, hỗ trợ trực tiếp cho việc thiết lập thông số vận hành lò phản ứng dạng ống và nồi áp suất đạt hiệu suất tối đa.

Thứ hai, các nhà khoa học và giảng viên chuyên ngành Hóa học tính toán: Tài liệu là minh chứng tiêu biểu về việc ứng dụng thành công phương pháp DFT B3LYP/6-31G* kết hợp phần mềm Gaussian để giải quyết bài toán động học phản ứng gốc tự do phức tạp trong công nghiệp, với độ tin cậy tương đương dữ liệu thực nghiệm đạt trên 95%.

Thứ ba, học viên cao học và nghiên cứu sinh chuyên ngành Hóa Hữu cơ, Hóa Dầu: Luận văn là tài liệu tham khảo chuẩn mực về phương pháp giải phương trình Schrödinger, lý thuyết trạng thái chuyển tiếp và quy trình mô hình hóa chuỗi polymer từ 2 đến 8 nguyên tử cacbon.

Thứ tư, doanh nghiệp sản xuất và gia công màng nhựa bao bì dẻo: Nắm bắt cơ chế tạo 4 đến 5 nhánh butyl/amyl trên 1000 đơn vị cacbon giúp các nhà sản xuất làm chủ công nghệ điều chỉnh độ nhớt dòng chảy, độ trong suốt và cơ tính của hạt nhựa LDPE thương phẩm phục vụ thị trường trong nước và xuất khẩu.

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

Phương pháp DFT B3LYP/6-31G* có ưu thế gì vượt trội trong nghiên cứu động học trùng hợp etylen? Phương pháp DFT với phiếm hàm lai hóa B3LYP và bộ cơ sở 6-31G* mang lại độ chính xác cao khi tính toán năng lượng hoạt hóa phản ứng phát triển mạch đạt 28 kJ/mol, độ sai lệch dưới 5% so với thực nghiệm. Đồng thời, phương pháp này xử lý xuất sắc tương quan điện tử trong các gốc tự do vỏ mở mà vẫn tiết kiệm đáng kể thời gian tính toán so với các phương pháp hậu Hartree-Fock phức tạp như CCSD(T).

Vì sao propanal lại có hiệu quả truyền mạch đạt tới 90,5% ở 300 độ C, cao hơn nhiều so với propan và propen? Propanal sở hữu liên kết C-H tại nhóm chức aldehyde (C=O) với năng lượng liên kết tương đối yếu, cho phép phản ứng lấy hydro của gốc etyl diễn ra với hàng rào năng lượng hoạt hóa thấp. Kết quả là phản ứng truyền mạch tạo thành sản phẩm gốc acyl chiếm tới 85% ở 220 độ C và đạt 90,5% ở 300 độ C, vượt trội hoàn toàn so với mức 1,8% của propan và 34,9% của propen.

Phản ứng cắn lưng (back-biting) ảnh hưởng trực tiếp đến cấu trúc và cơ tính của nhựa LDPE như thế nào? Phản ứng cắn lưng là quá trình gốc tự do đầu mạch tự tách nguyên tử hydro ở vị trí cacbon số 4 hoặc số 5 trên chính chuỗi của nó thông qua trạng thái chuyển tiếp vòng 6 hoặc 7 cạnh. Cơ chế này tạo ra 4 đến 5 nhánh ngắn butyl và amyl trên 1000 nhóm CH2 ở 300 độ C, ngăn cản sự xếp chặt của các chuỗi phân tử, từ đó làm giảm tỷ trọng và tăng tính mềm dẻo đặc trưng cho LDPE.

Luận văn đã giải quyết những rào cản kỹ thuật nào của phương pháp nghiên cứu thực nghiệm truyền thống? Phản ứng tổng hợp LDPE diễn ra ở điều kiện cực kỳ khắc nghiệt với áp suất 1000 đến 3000 atm, nhiệt độ trên 200 độ C và thời gian sống của gốc tự do chỉ dưới 10 mũ trừ 6 giây. Phương pháp tính toán lượng tử đã giải quyết triệt để khó khăn này bằng cách tái lập cấu trúc phân tử và trạng thái chuyển tiếp trên máy tính một cách an toàn, chính xác và tiết kiệm chi phí.

Mức độ tương thích giữa kết quả tính toán trong luận văn với dữ liệu thực tế của ngành công nghiệp như thế nào? Dữ liệu tính toán chỉ ra rằng sản phẩm lấy hydro ở nhóm carbonyl của propanal đạt 85% ở 220 độ C, hoàn toàn tương thích với các kết quả đo đạc thực nghiệm từ tập đoàn hóa chất DSM. Ngoài ra, các mô hình động học cũng phản ánh chính xác định mức nguyên liệu tiêu hao từ 1005 đến 1010 kg etylen cho 1 tấn LDPE tại các nhà máy công nghiệp.

Kết luận

  • Luận văn đã xây dựng thành công quy trình ứng dụng phương pháp hóa học tính toán B3LYP/6-31G* để mô phỏng chính xác cơ chế và động học phản ứng trùng hợp polyetylen tỷ trọng thấp.
  • Xác định rõ năng lượng hoạt hóa của phản ứng phát triển mạch là 28 kJ/mol và thiết lập phương trình hàm mũ biểu diễn hằng số tốc độ Kp theo nhiệt độ.
  • Chứng minh propanal là tác nhân truyền mạch ưu việt nhất với hiệu quả đạt 90,5% ở 300 độ C thông qua cơ chế tách hydro chọn lọc tại nhóm carbonyl.
  • Làm sáng tỏ bản chất phản ứng cắn lưng nội phân tử tạo ra 4 đến 5 nhánh butyl và amyl trên 1000 đơn vị cacbon ở 300 độ C, giải thích cấu trúc vi mô của LDPE.
  • Đóng góp luận cứ khoa học quan trọng giúp tối ưu hóa định mức tiêu hao 1005 kg etylen/tấn sản phẩm và kiểm soát chất lượng hạt nhựa trong công nghiệp hóa dầu.

Kế hoạch phát triển tiếp theo là mở rộng áp dụng mô hình tính toán DFT cho các quá trình đồng trùng hợp etylen với các comonome phân cực trong giai đoạn 2024 - 2026. Hãy khai thác và ứng dụng ngay phương pháp hóa học tính toán hiện đại để nâng cao năng lực nghiên cứu, cải tiến quy trình công nghệ và tối ưu hóa hiệu quả sản xuất vật liệu cao phân tử.