Tổng quan nghiên cứu

Nước chiếm khoảng 70% bề mặt Trái Đất và đóng vai trò dung môi thiết yếu trong mọi cơ thể sống. Trong điều kiện khối thông thường, nước sở hữu các tính chất nhiệt động quen thuộc như nhiệt độ sôi ở 100 độ C và nhiệt độ đóng băng ở 0 độ C tại áp suất khí quyển tiêu chuẩn 1 bar. Tuy nhiên, khi bị giam hãm trong không gian hẹp cấp độ nanomet dưới 2 nm, các liên kết hydro và cấu trúc mạng lưới phân tử bị biến đổi sâu sắc, dẫn đến những hành vi chuyển pha dị thường chưa từng được ghi nhận ở vật liệu dạng khối.

Vấn đề cốt lõi của nghiên cứu là làm sáng tỏ cơ chế chuyển pha của nước bên trong ống nano carbon đơn lớp khi nhiệt độ hạ từ 300 K xuống 200 K. Mục tiêu cụ thể là phân tích sự biến đổi của độ lệch bình phương trung bình và nhiệt dung mol dưới tác động của 3 chế độ làm lạnh khác nhau cùng 3 mức áp suất môi trường gồm 0 bar, 1 bar và 2 bar.

Đề tài được thực hiện trong khung thời gian từ tháng 9 năm 2021 đến tháng 12 năm 2022 tại Phòng thí nghiệm Vật lý Tính toán thuộc Trường Đại học Bách Khoa, Đại học Quốc gia Thành phố Hồ Chí Minh. Tác giả Lâm Toan Vĩ dưới sự hướng dẫn khoa học của Phó Giáo sư Trần Thị Thu Hạnh đã hoàn thành và bảo vệ thành công công trình vào tháng 2 năm 2023.

Nghiên cứu mang ý nghĩa then chốt trong việc cung cấp bức tranh nhiệt động học chi tiết về chất lỏng bị giam hãm, tạo tiền đề định lượng cho việc thiết kế các kênh dẫn vi lưu nano, màng lọc khử muối với đường kính ống quanh mức 14.5 Angstrom và các thiết bị nano phục vụ công nghệ sinh học trong tương lai.

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 nhiệt động lực học thống kê và các mô hình chuyển pha cổ điển kết hợp hiện đại:

  • Tiêu chuẩn chuyển pha Ehrenfest: Phân loại chuyển pha dựa trên đạo hàm của thế nhiệt động. Chuyển pha loại một thể hiện sự nhảy bậc gián đoạn của đạo hàm bậc một như entropy và thể tích, đi kèm nhiệt lượng chuyển pha xác định. Chuyển pha loại hai đặc trưng bởi sự biến thiên liên tục của entropy nhưng lại xuất hiện bước nhảy đột ngột ở các đạo hàm bậc hai như nhiệt dung đẳng áp, hệ số nở nhiệt và hệ số nén đẳng nhiệt.
  • Tiêu chuẩn hiện đại và lý thuyết Landau: Bổ sung khái niệm tham số trật tự và giải thích hiện tượng mất trật tự ngẫu nhiên trong cấu trúc. Khi tham số trật tự biến đổi liên tục về 0, hệ trải qua quá trình chuyển pha loại hai với vùng chuyển tiếp kéo dài.
  • Tương tác thế vi mô: Mô hình hóa tương tác giữa các nguyên tử qua hàm thế Lennard-Jones 12-6 cho lực phân tán Van der Waals tầm ngắn và định luật Coulomb cho tương tác tĩnh điện tầm xa giữa các điện tích điểm hiệu dụng.

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

Nghiên cứu sử dụng phương pháp mô phỏng động lực học phân tử thông qua bộ công cụ mã nguồn mở LAMMPS, tích hợp các thuật toán tích phân phương trình chuyển động Newton với bước thời gian 1 femtogiây:

  • Cỡ mẫu và cách thức chọn mẫu: Hệ mô phỏng gồm các ống nano carbon đơn lớp có chỉ số mạng lần lượt là (14,14), (15,15), (16,16) và (17,17) với đường kính dao động từ 1.4 nm đến 2.3 nm, chứa từ 300 đến hơn 1200 phân tử nước. Mẫu cấu hình ban đầu được phân bố ngẫu nhiên và tối thiểu hóa năng lượng tĩnh bằng thuật toán suy giảm độ dốc Steepest Descent trước khi đưa vào trạng thái cân bằng nhiệt động.
  • Mô hình nước: Lựa chọn mô hình nước 3 điểm TIP3P với độ dài liên kết O-H cố định 0.9572 Angstrom và góc liên kết H-O-H là 104.52 độ, điện tích hiệu dụng nguyên tử oxy là -0.834e và hydro là +0.417e.
  • Điều kiện biên và tập hợp thống kê: Áp dụng điều kiện biên tuần hoàn dọc theo trục đối xứng của ống nano để tái tạo hệ vô hạn, sử dụng tập hợp chính tắc NVT và tập hợp đẳng áp đẳng nhiệt NPT. Toàn bộ quy trình mô phỏng hạ nhiệt từ 300 K xuống 200 K diễn ra trong vòng 15 tháng nghiên cứu. Phương pháp mô phỏng động lực học phân tử được lựa chọn vì đây là công cụ tối ưu cho phép theo dõi trực tiếp quỹ đạo vi mô của từng nguyên tử, điều mà các kỹ thuật thực nghiệm đo đạc ở cấp độ nanomet hiện nay rất khó can thiệp.

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

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

  1. Tốc độ làm lạnh quyết định tính chất đột ngột của quá trình đông đặc: Khi áp dụng tốc độ làm lạnh nhanh, sự chuyển pha của nước diễn ra tức thời và tập trung trong một dải nhiệt độ hẹp từ 220 K đến 240 K. Độ lệch bình phương trung bình sụt giảm đột ngột hơn 80%, phản ánh sự suy giảm tức thì của khả năng khuếch tán phân tử.
  2. Tốc độ làm lạnh chậm tạo ra các cấu trúc băng giả bền: Khi hạ nhiệt từ từ, hệ nước trải qua quá trình chuyển pha liên tục. Năng lượng tự do và enthalpy của hệ xuất hiện nhiều điểm cực tiểu cục bộ, hình thành các pha tinh thể giả bền trước khi đạt tới trạng thái năng lượng thấp nhất.
  3. Áp suất môi trường thấp thúc đẩy chuyển pha dứt khoát: Tại mức áp suất 0 bar và 1 bar, sự nhảy bậc nhiệt động diễn ra rõ nét hơn khoảng 35% so với mức áp suất 2 bar. Khi tăng áp suất lên 2 bar, quá trình chuyển pha có xu hướng trở nên mượt mà và dải nhiệt độ chuyển tiếp bị nới rộng.
  4. Kích thước ống định hình cấu trúc tinh thể băng: Đường kính ống nano carbon quanh ngưỡng 14.5 Angstrom quyết định trực tiếp việc hình thành mạng lưới nước đá dạng một lớp, hai lớp hay ba lớp đồng trục tại nhiệt độ 0 K.

Thảo luận kết quả

Cơ chế vật lý của hiện tượng bắt nguồn từ sự cạnh tranh giữa tương tác nước-nước và tương tác nước-carbon. Dưới tốc độ làm lạnh nhanh, các phân tử nước không có đủ thời gian hồi phục cấu trúc để sắp xếp vào các vị trí giả bền, buộc toàn bộ hệ thống phải chuyển trạng thái tập thể đồng loạt.

Dữ liệu mô phỏng được biểu diễn trực quan qua hệ thống đồ thị độ lệch bình phương trung bình theo nhiệt độ và đồ thị nhiệt dung riêng theo nhiệt độ. Trên đồ thị khuếch tán, các đường cong chuyển từ trạng thái dốc đứng sang phương ngang phẳng lì tại điểm hóa rắn, đồng thời đồ thị nhiệt dung riêng ghi nhận các đỉnh nhọn cực đại tại lân cận nhiệt độ tới hạn, minh chứng cho hiện tượng chuyển pha loại hai theo tiêu chuẩn Landau. Bảng tổng hợp các điểm chuyển pha xác nhận sự dịch chuyển rõ rệt của nhiệt độ đông đặc khi thay đổi tốc độ hạ nhiệt từ nhanh sang chậm.

So sánh với các công bố của nhóm tác giả Koga trong điều kiện áp suất siêu cao trên 50 MPa và nghiên cứu của tác giả Walther về dòng chảy mao dẫn ở 120 bar đến 1000 bar, công trình này chứng minh rằng ngay cả ở áp suất thấp từ 0 đến 2 bar, hiệu ứng giam hãm nano vẫn tạo ra các cấu trúc băng hình học dị thường mà dạng khối không thể có được.

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

  • Nâng cấp mô hình tương tác phân tử: Nhóm nghiên cứu thuộc Phòng thí nghiệm Vật lý Tính toán cần tích hợp các mô hình nước 4 điểm hoặc 5 điểm tiên tiến như TIP4P-Ew và OPC trong vòng 6 đến 12 tháng tới nhằm giảm sai số tính toán phân cực điện môi xuống dưới 5%.
  • Mở rộng phạm vi khảo sát điều kiện biên: Các nghiên cứu sinh chuyên ngành Vật lý kỹ thuật cần mở rộng dải áp suất khảo sát lên 50 MPa đến 500 MPa và kiểm tra dải tốc độ làm lạnh từ 10 mũ 9 K/s đến 10 mũ 11 K/s trong giai đoạn 2024 đến 2025 để hoàn thiện giản đồ pha toàn diện.
  • Ứng dụng vào chế tạo màng lọc nano: Các kỹ sư công nghệ vật liệu và môi trường nên ứng dụng dữ liệu chuyển pha để tối ưu hóa thiết kế màng thẩm thấu ngược bằng ống nano carbon, hướng tới mục tiêu tăng thông lượng dòng chảy nước thêm 25% đến 30% trước năm 2026.
  • Phát triển thuật toán dự đoán cấu trúc băng: Nhóm chuyên gia khoa học tính toán cần xây dựng mô hình học máy kết hợp động lực học phân tử nhằm dự đoán chính xác trên 95% trạng thái tinh thể giả bền của nước trong các kênh nano phức tạp trong vòng 18 tháng.

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

  • Học viên cao học và nghiên cứu sinh chuyên ngành Vật lý Kỹ thuật, Khoa học Vật liệu: Luận văn cung cấp quy trình chuẩn mực về thiết lập mô phỏng động lực học phân tử bằng LAMMPS, phương pháp tối thiểu hóa năng lượng và kỹ thuật phân tích chuyển pha nhiệt động.
  • Kỹ sư phát triển màng lọc công nghệ cao và vật liệu nano: Tài liệu mang lại dữ liệu nền tảng về động học dòng chảy và hành vi đóng băng của nước trong ống carbon, hỗ trợ tối ưu hóa thiết kế màng lọc khử mặn và xử lý nước tinh khiết.
  • Chuyên gia mô phỏng tính toán vi mô và hóa lý: Cung cấp góc nhìn so sánh chuyên sâu giữa các mô hình nước TIP3P, SPC/E, TIP4P, TIP5P và ảnh hưởng của các tập hợp thống kê NVT, NPT trong hệ phi cân bằng.
  • Giảng viên và nhà nghiên cứu vật lý chất rắn, nhiệt động học nano: Nguồn tư liệu học thuật giá trị để xây dựng bài giảng chuyên đề về hiện tượng chuyển pha trong không gian hẹp và vật lý màng mỏng.

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

  • Chuyển pha của nước trong ống nano carbon khác biệt thế nào so với nước khối thông thường? Trong không gian hẹp dưới 2 nm, nước không đóng băng cố định ở 273.15 K mà trải qua quá trình kết tinh trong khoảng 200 K đến 300 K. Quá trình này có thể diễn ra liên tục qua nhiều trạng thái giả bền tùy thuộc vào đường kính ống 14.5 Angstrom thay vì biến đổi tức thời như nước dạng khối.

  • Tại sao nghiên cứu lại lựa chọn mô hình nước 3 điểm TIP3P trong mô phỏng LAMMPS? Mô hình TIP3P với góc liên kết 104.52 độ giúp giảm thiểu chi phí và thời gian tính toán khoảng 40% so với mô hình 5 điểm phức tạp. Đây là bước tiếp cận nền tảng tối ưu để đánh giá tổng thể hành vi chuyển pha trước khi mở rộng sang các mô hình phân cực cao hơn.

  • Tốc độ làm lạnh nhanh ảnh hưởng như thế nào đến trạng thái đông đặc của nước? Làm lạnh nhanh khiến thời gian hồi phục cấu trúc của phân tử nước bị rút ngắn, dẫn đến hiện tượng chuyển pha đột ngột trong dải nhiệt độ hẹp từ 220 K đến 240 K. Độ lệch bình phương trung bình sụt giảm trên 80% chỉ trong vài pico-giây, phản ánh sự đóng băng tức thời.

  • Áp suất môi trường từ 0 bar đến 2 bar đóng vai trò gì trong quá trình hóa rắn nano? Áp suất thấp từ 0 bar đến 1 bar tạo điều kiện cho các phân tử nước sắp xếp trật tự dứt khoát hơn, làm bước nhảy nhiệt động sắc nét hơn 35%. Khi tăng áp suất lên 2 bar, tương tác đẩy giữa thành ống và nước tăng lên, làm mượt đường cong chuyển pha.

  • Đồ thị MSD và nhiệt dung riêng thể hiện điểm chuyển pha qua dấu hiệu nào? Trên đồ thị MSD, điểm chuyển pha là vị trí đường cong gãy khúc đột ngột và chuyển sang trạng thái nằm ngang phẳng. Trên đồ thị nhiệt dung riêng theo nhiệt độ, điểm chuyển pha được xác định bởi đỉnh cực đại cục bộ tương ứng với sự thay đổi đạo hàm nhiệt động theo tiêu chuẩn Ehrenfest loại hai.

Kết luận

  • Luận văn đã mô phỏng thành công quá trình làm lạnh nước từ 300 K xuống 200 K trong ống nano carbon đơn lớp bằng phần mềm LAMMPS.
  • Tốc độ làm lạnh nhanh tạo ra hiện tượng chuyển pha đột ngột và tập trung trong khoảng nhiệt độ hẹp từ 220 K đến 240 K.
  • Tốc độ làm lạnh chậm dẫn đến quá trình chuyển pha liên tục và hình thành các cấu trúc tinh thể băng giả bền.
  • Áp suất môi trường thấp từ 0 bar đến 1 bar thúc đẩy sự chuyển pha diễn ra dứt khoát hơn so với mức áp suất 2 bar.
  • Xác lập mối liên hệ định lượng chặt chẽ giữa độ lệch bình phương trung bình và đỉnh nhiệt dung mol trong việc xác định điểm chuyển pha.

Đóng góp lớn nhất của công trình là làm sáng tỏ bản chất nhiệt động học của nước bị giam hãm trong không gian nano dưới các tốc độ làm lạnh khác nhau. Kế hoạch tiếp theo là mở rộng khảo sát ở dải áp suất cao 50 MPa đến 500 MPa và ứng dụng mô hình nước TIP4P-Ew trong giai đoạn 2024-2025. Các nhà khoa học và kỹ sư công nghệ nano có thể khai thác bộ số liệu này để phát triển thế hệ màng lọc khử muối và linh kiện vi lưu nano tiên tiến.