Bốn thế nhiệt động cơ bản trong vật lý là những hàm nào?

Bốn thế nhiệt động cơ bản gồm nội năng U(S,V), năng lượng tự do Helmholtz F(T,V), enthalpy H(S,P) và năng lượng tự do Gibbs G(T,P).

Phép biến đổi Legendre có vai trò gì đối với các thế nhiệt động?

Phép biến đổi Legendre cho phép chuyển đổi giữa các biến số tự nhiên liên hợp mà vẫn bảo toàn nguyên vẹn toàn bộ nội dung thông tin nhiệt động lực học của hàm gốc.

Tại sao năng lượng tự do Gibbs được sử dụng phổ biến nhất trong hóa học?

Vì năng lượng tự do Gibbs có các biến tự nhiên là nhiệt độ và áp suất, phù hợp với điều kiện môi trường thực nghiệm thông thường của hầu hết các phản ứng hóa học.

Các hệ thức Maxwell được suy ra từ nguyên lý toán học nào?

Các hệ thức Maxwell được suy ra trực tiếp từ định lý Schwarz về sự bằng nhau của các đạo hàm riêng hỗn hợp bậc hai của các hàm thế nhiệt động khả vi liên tục.

Thế nhiệt là gì? Bản chất các hàm thế nhiệt động

Thế nhiệt (thế nhiệt động hay hàm thế nhiệt động) là đại lượng vô hướng trong nhiệt động lực học biểu diễn trạng thái năng lượng toàn phần hoặc phần năng lượng hữu ích có thể chuyển đổi của một hệ vật lý, có xu hướng đạt giá trị cực trị tại trạng thái cân bằng dưới các điều kiện ràng buộc xác định của môi trường. Các hàm thế nhiệt động đóng vai trò như các hàm thế trong cơ học cổ điển, cho phép xác định chiều hướng tự diễn biến của các quá trình vật lý, hóa học cũng như tính toán các thông số cân bằng pha mà không cần theo dõi chi tiết vi mô của từng phân tử. Bài viết này trình bày hệ thống về cơ sở toán học của phép biến đổi Legendre, bốn thế nhiệt động kinh điển, các hệ thức Maxwell, cấu trúc hình học trạng thái và những ứng dụng hiện đại trong khoa học vật liệu.

Bản chất vật lý và vai trò của các hàm thế nhiệt động

Trong cơ học cổ điển, một chất điểm hoặc hệ chất điểm chịu tác dụng của trường lực thế sẽ tự động dịch chuyển về vị trí có thế năng cực tiểu để đạt trạng thái cân bằng cơ học bền. Hoàn toàn tương tự, một hệ nhiệt động vĩ mô trao đổi năng lượng và vật chất với môi trường xung quanh sẽ tự phát biến đổi theo hướng làm tăng tổng entropy của vũ trụ theo nguyên lý thứ hai của nhiệt động lực học. Khi hệ được đặt dưới các điều kiện biên cố định từ môi trường bên ngoài như nhiệt độ, áp suất hoặc thể tích không đổi, xu hướng tăng entropy toàn phần tương đương với xu hướng tiến tới cực trị của một hàm trạng thái năng lượng tương ứng của hệ. Hàm trạng thái này được gọi là thế nhiệt động.

Giá trị và sự thay đổi của các thế nhiệt động cung cấp thông tin hoàn chỉnh về toàn bộ các tính chất nhiệt động học của hệ. Đạo hàm riêng bậc nhất của thế nhiệt động theo các biến tự nhiên xác định các thông số trạng thái liên hợp như áp suất, nhiệt độ hoặc thế hóa học. Đạo hàm riêng bậc hai xác định các hệ số phản ứng vĩ mô quan trọng như nhiệt dung, độ nén đẳng nhiệt và hệ số giãn nở nhiệt.

Bốn thế nhiệt động cơ bản và các biến số tự nhiên

Theo phân tích chuẩn hóa của Alberty (2001) trong báo cáo kỹ thuật của Hiệp hội Hóa học Thuần túy và Ứng dụng Quốc tế (IUPAC), có 4 thế nhiệt động cơ bản được sử dụng rộng rãi nhất trong vật lý và hóa học, mỗi thế tương ứng với một bộ biến số tự nhiên riêng biệt:

Thế nhiệt động	Ký hiệu	Định nghĩa toán học	Biến số tự nhiên	Điều kiện cân bằng cực tiểu
Nội năng	$U$	Thế năng và động năng vi mô toàn phần	Entropy $S$ , Thể tích $V$	Hệ cô lập ( $S$ , $V$ không đổi)
Năng lượng tự do Helmholtz	$F$ (hoặc $A$ )	$F = U - TS$	Nhiệt độ $T$ , Thể tích $V$	Hệ đẳng nhiệt, đẳng tích ( $T$ , $V$ không đổi)
Enthalpy	$H$	$H = U + PV$	Entropy $S$ , Áp suất $P$	Hệ đẳng áp, đoạn nhiệt ( $S$ , $P$ không đổi)
Năng lượng tự do Gibbs	$G$	$G = H - TS = U + PV - TS$	Nhiệt độ $T$ , Áp suất $P$	Hệ đẳng nhiệt, đẳng áp ( $T$ , $P$ không đổi)

Ý nghĩa vật lý của từng hàm thế

Nội năng ( $U$ ): Đại diện cho toàn bộ năng lượng dự trữ bên trong hệ gồm động năng chuyển động nhiệt hỗn loạn của các phân tử và thế năng tương tác giữa các hạt cấu thành. Dạng vi phân cơ bản đối với hệ một cấu tử là $dU = T dS - P dV + \mu dN$ , trong đó biến số tự nhiên là các đại lượng dung độ gồm entropy, thể tích và số hạt.
Năng lượng tự do Helmholtz ( $F$ ): Thể hiện phần năng lượng tối đa của hệ có thể chuyển hóa thành toàn bộ công cơ học cực đại (bao gồm công giãn nở thể tích) trong một quá trình biến đổi đẳng nhiệt. Phương trình vi phân là $dF = -S dT - P dV + \mu dN$ . Đây là hàm thế then chốt trong vật lý chất rắn và vật lý thống kê chính tắc.
Enthalpy ( $H$ ): Đo lường tổng nhiệt lượng trao đổi của hệ trong các quá trình biến đổi đẳng áp vốn phổ biến trong các lò phản ứng hở hoặc buồng đốt khí quyển. Phương trình vi phân là $dH = T dS + V dP + \mu dN$ .
Năng lượng tự do Gibbs ( $G$ ): Là thế nhiệt động quan trọng hàng đầu trong hóa học và địa chất học vì hầu hết các phản ứng hóa học và chuyển hóa sinh học diễn ra trong điều kiện nhiệt độ và áp suất môi trường không đổi. Phương trình vi phân là $dG = -S dT + V dP + \mu dN$ . Khi một phản ứng tự diễn biến ở nhiệt độ và áp suất cố định, độ biến thiên $\Delta G$ luôn mang giá trị âm, và hệ đạt trạng thái cân bằng bền khi $G$ đạt giá trị cực tiểu.

Cơ sở toán học của phép biến đổi Legendre

Một câu hỏi nền tảng đặt ra là tại sao nhiệt động lực học cần nhiều hàm thế năng lượng khác nhau thay vì chỉ sử dụng duy nhất một hàm nội năng $U$ . Công trình nghiên cứu của Zia và các cộng sự (2009) đã làm sáng tỏ bản chất toán học sâu sắc của vấn đề này thông qua phép biến đổi Legendre.

Nếu ta chỉ đơn thuần thay đổi biến số độc lập bằng phép thế đại số thông thường, ví dụ thay biến entropy $S$ bằng nhiệt độ $T = (\partial U / \partial S)_V$ trực tiếp vào hàm $U$ , ta sẽ làm mất mát thông tin nhiệt động lực học quan trọng và phương trình thu được không còn là một phương trình trạng thái cơ bản. Để thay thế một biến số dung độ mở rộng bằng một biến số cường độ liên hợp mà vẫn bảo toàn nguyên vẹn toàn bộ nội dung thông tin của hàm gốc, phép biến đổi Legendre là công cụ toán học duy nhất thỏa mãn yêu cầu này.

Đối với hàm nội năng $U(S, V)$ , phép biến đổi Legendre theo biến entropy $S$ dẫn đến hàm năng lượng tự do Helmholtz:

F(T, V) = U(S, V) - S \left(\frac{\partial U}{\partial S}\right)_V = U - TS

Tương tự, phép biến đổi Legendre của nội năng theo biến thể tích $V$ sinh ra enthalpy:

H(S, P) = U(S, V) - V \left(\frac{\partial U}{\partial V}\right)_S = U - V(-P) = U + PV

Và phép biến đổi Legendre đồng thời trên cả hai biến $S$ và $V$ tạo ra năng lượng tự do Gibbs:

G(T, P) = U - TS + PV

Các hệ thức Maxwell kinh điển

Một hệ quả toán học tuyệt đẹp xuất phát trực tiếp từ tính chất liên tục và khả vi của các hàm thế nhiệt động là định lý Schwarz về sự độc lập của đạo hàm riêng hỗn hợp bậc hai đối với thứ tự lấy vi phân. Từ bốn hàm thế cơ bản, ta suy ra bốn hệ thức Maxwell kinh điển:

Từ vi phân toàn phần của nội năng $dU$ :
$\left(\frac{\partial T}{\partial V}\right)_S = -\left(\frac{\partial P}{\partial S}\right)_V$
Từ vi phân toàn phần của năng lượng tự do Helmholtz $dF$ :
$\left(\frac{\partial S}{\partial V}\right)_T = \left(\frac{\partial P}{\partial T}\right)_V$
Từ vi phân toàn phần của enthalpy $dH$ :
$\left(\frac{\partial T}{\partial P}\right)_S = \left(\frac{\partial V}{\partial S}\right)_P$
Từ vi phân toàn phần của năng lượng tự do Gibbs $dG$ :
$-\left(\frac{\partial S}{\partial P}\right)_T = \left(\frac{\partial V}{\partial T}\right)_P$

Ý nghĩa thực tiễn to lớn của các hệ thức Maxwell nằm ở chỗ chúng cho phép liên kết các đại lượng rất khó hoặc không thể đo lường trực tiếp trong phòng thí nghiệm (như sự biến thiên của entropy theo áp suất hoặc thể tích) với các đại lượng vĩ mô hoàn toàn có thể đo đạc dễ dàng bằng áp kế và nhiệt kế (như độ dốc áp suất theo nhiệt độ hoặc hệ số nở khối nhiệt).

Cấu trúc hình học Riemann và độ ổn định nhiệt động

Không gian trạng thái của các hệ nhiệt động không chỉ là một tập hợp các biến số rời rạc mà sở hữu một cấu trúc hình học vi phân nội tại chặt chẽ. Trong công trình tiên phong của Weinhold (1975), tác giả đã đề xuất xây dựng một cấu trúc hình học Riemann trên đa tạp nhiệt động lực học bằng cách định nghĩa tenxơ mêtric xuất phát từ ma trận Hessian của hàm nội năng hoặc các thế nhiệt động tương đương.

Theo Weinhold (1975), ma trận Hessian của hàm nội năng đóng vai trò là tenxơ mêtric $g_{ij}$ xác lập cấu trúc hình học Riemann trên không gian trạng thái nhiệt động lực học. Để một hệ nhiệt động duy trì trạng thái cân bằng bền vững trước các nhiễu loạn vi mô ngẫu nhiên, mặt thế năng lượng phải thỏa mãn các điều kiện lồi toán học nghiêm ngặt:

Nội năng $U(S, V)$ phải là một hàm lồi theo cả hai biến entropy và thể tích, đòi hỏi nhiệt dung đẳng tích và độ nén đoạn nhiệt luôn mang giá trị dương.
Năng lượng tự do Gibbs $G(T, P)$ phải là một hàm lõm theo các biến cường độ nhiệt độ và áp suất, nghĩa là ma trận Hessian của $G$ theo $T$ và $P$ phải là nửa xác định âm.

Bất kỳ sự vi phạm nào đối với các điều kiện lồi lõm này đều báo hiệu sự mất ổn định nhiệt động lực học nội tại của hệ, kích hoạt hiện tượng phân tách pha tự phát hoặc chuyển pha bậc một.

Cấu trúc topo mặt thế và hiện tượng chuyển pha

Nghiên cứu hiện đại của Nitzke và cộng sự (2024) công bố trên tạp chí The Journal of Chemical Physics đã phân tích chi tiết cấu trúc topo của các mặt thế nhiệt động bao gồm thế Helmholtz, thế Gibbs và thế Grand Potential thông qua các mô hình vật lý chính xác. Các tác giả chỉ ra rằng hình thái học topo của các mặt thế năng lượng biến đổi đặc trưng khi hệ tiến gần đến các điểm kỳ dị chuyển pha.

Tại các ranh giới cùng tồn tại pha, mặt thế năng lượng hình thành các nếp gấp topo nơi các nhánh nghiệm khác nhau giao nhau. Việc biểu diễn trực quan và giải tích chính xác các mặt thế này phân tích topo của các thế nhiệt động toàn phần và phần đóng góp dư trên các chất lưu nguyên chất (như metan, etan, nitơ và metanol), giúp làm sáng tỏ mối liên hệ sâu sắc giữa tương tác liên phân tử và dạng hình học của các mặt thế nhiệt qua các tập hợp thống kê khác nhau.

Mở rộng trong hóa học lượng tử và lý thuyết hàm mật độ

Khái niệm hàm thế nhiệt động và phép biến đổi Legendre không chỉ dừng lại ở phạm vi nhiệt động học vĩ mô truyền thống mà còn được mở rộng thành công sang cơ học lượng tử vi mô. Công trình kinh điển của Nalewajski và Parr (1982) đã áp dụng thành công phép biến đổi Legendre và các hệ thức Maxwell vào khuôn khổ lý thuyết hàm mật độ điện tử.

Theo Nalewajski và Parr (1982), năng lượng trạng thái cơ bản của một hệ electron đóng vai trò tương tự như thế nhiệt động, trong đó thế hóa học điện tử và độ cứng hóa học lượng tử được định nghĩa một cách chặt chẽ thông qua các đạo hàm của năng lượng theo số lượng electron và thế năng hạt nhân bên ngoài. Mở rộng này thiết lập một cầu nối giải tích vững chắc giữa các định luật nhiệt động lực học vĩ mô và cấu trúc liên kết hóa học vi mô của phân tử.

Giới hạn áp dụng và các trường hợp đặc biệt

Dù là công cụ trung tâm của vật lý nhiệt và hóa lý, lý thuyết thế nhiệt động có những biên giới áp dụng rõ rệt:

Hệ kích thước nano và giới hạn nhiệt động học: Khung lý thuyết thế nhiệt động được xây dựng dựa trên giả định giới hạn nhiệt động học, khi số hạt và thể tích của hệ tiến tới vô cùng. Đối với các hệ cụm hạt nano hoặc phân tử đơn lẻ, mức độ thăng giáng năng lượng trở nên rất lớn so với giá trị trung bình, khiến các cực tiểu của thế nhiệt động không còn sắc nét và mất tính tất định vĩ mô.
Hệ phi cân bằng xa: Các hàm thế nhiệt động cổ điển chỉ xác định chặt chẽ cho các trạng thái cân bằng hoặc các quá trình biến đổi chuẩn tĩnh. Đối với các hệ phi cân bằng xa (như các cấu trúc tiêu tán tự tổ chức), các hàm thế này không thể áp dụng trực tiếp để dự đoán chiều hướng diễn tiến nếu không có sự mở rộng của nhiệt động học phi tuyến.
Kỳ dị toán học tại điểm chuyển pha: Tại các điểm chuyển pha liên tục hoặc bậc một, các đạo hàm riêng bậc cao của thế nhiệt động có thể phân kỳ hoặc gián đoạn, làm phá vỡ tính khả vi trơn của đa tạp năng lượng và khiến tenxơ metric Hessian mất tính xác định dương.
Hệ có tương tác tầm xa: Các hệ chịu lực tương tác tầm xa (như lực hấp dẫn thiên văn hoặc plasma không bị chắn điện) vi phạm tính chất cộng dồn dung độ, dẫn tới các hiện tượng nghịch lý như nhiệt dung âm và sự mất tương đương giữa các tập hợp thống kê vi chính tắc và chính tắc.

Ứng dụng thực tiễn trong khoa học và công nghệ hiện đại

Các hàm thế nhiệt động là công cụ tính toán cơ sở trong hầu hết các lĩnh vực kỹ thuật hiện đại:

Kỹ thuật hóa chất và lọc hóa dầu: Năng lượng tự do Gibbs là nền tảng để tính toán hằng số cân bằng phản ứng hóa học, hiệu suất chuyển hóa tối đa và thiết kế tháp chưng cất tách hỗn hợp hydrocacbon.
Luyện kim và khoa học vật liệu: Phương pháp tính toán giản đồ pha CALPHAD sử dụng mô hình hóa các hàm năng lượng tự do Gibbs của các pha rắn và lỏng để dự đoán cấu trúc tinh thể của các siêu hợp kim mới chịu nhiệt độ cao.
Năng lượng tái tạo và pin điện hóa: Suất điện động tiêu chuẩn của pin nhiên liệu hydro và pin lithium-ion được tính toán trực tiếp từ độ biến thiên năng lượng tự do Gibbs qua hệ thức $\Delta G^\circ = -nFE^\circ$ , và hiệu chỉnh theo điều kiện hoạt độ thực tế bằng phương trình Nernst.
Mô phỏng phân tử động lực học: Thế nhiệt Helmholtz và Grand Potential được sử dụng trực tiếp trong các tập hợp thống kê NVT và muVT để tính toán áp suất thẩm thấu, hấp phụ khí trên vật liệu xốp nano và cuộn gập protein sinh học.