Nhiệt động lực học thống kê cho chất tinh khiết lý tưởng và các trạng thái tập hợp của nó

<a id="gesamt"></a>

Tóm tắt tổng thể

Mô hình 3 là một mô hình cơ học thống kê định tính của một chất lý tưởng, suy ra ba trạng thái tập hợp rắn, lỏng, khí cùng các chuyển pha giữa chúng từ một bức tranh liên kết duy nhất với các liên kết phân tử. Nguyên tắc nền tảng nằm ở việc tính đến ba kho chứa năng lượng cho nội năng: năng lượng dao động trong các liên kết phân tử, năng lượng tịnh tiến của các phân tử, và năng lượng chuyển động hợp tác trong sự quay của các mixen như một dạng lý tưởng hóa của các dòng vật chất bên trong nhằm san bằng các thăng giáng mật độ nội tại. Một đại lượng cơ bản quan trọng trong mối liên hệ này là số phối trí trung bình z̄ = z_max·x, đại lượng này tự thiết lập một cách tự hợp qua một liên kết hai trạng thái hợp tác (kiểu Ising, độ hợp tác K = J·z_max, độ suy biến của trạng thái đứt G = z_max). Từ nó suy ra ba chuyển pha trên một thang τ chung, được hiệu chuẩn bằng một tham số năng lượng ε duy nhất: nóng chảy như ngưỡng phonon/độ cứng (z̄ đạt tới độ cứng vững), sôi như sự sụp đổ tính liên thông (thấm) tại điểm tới hạn, và một chuyển pha tùy chọn rắn(1)→rắn(2) (FCC↔BCC) với bậc c_v không tham số bằng 1,0 k_B. Một chuyển pha trở thành bậc nhất ngay khi ln G > 2 + 4ε/K. Trong đó độ suy biến trạng thái đứt G là số tiếp xúc lân cận có thể bị bẻ gãy (G = z_max, số phối trí cực đại), năng lượng liên kết ε là năng lượng thu được trên mỗi tiếp xúc lân cận, và độ hợp tác K = J·z_max là tích của cường độ ghép cặp kiểu Ising J với z_max; như vậy chuyển pha càng dễ trở thành bậc nhất khi độ đa dạng tiếp xúc G càng lớn và độ hợp tác K càng mạnh so với năng lượng liên kết ε.

Bài viết dẫn từ bảy tiên đề (Chương 1) qua phân tích mâu thuẫn (Chương 2) và một mô tả bằng hình ảnh với các dải truyện tranh tĩnh (Chương 3) tới khung cơ học thống kê — liên kết hợp tác, tổng trạng thái và nhiệt dung (Chương 4–5) —, tới tính hợp lý toán học trong 3D (Chương 6), và cuối cùng tới hiệu chuẩn và phương trình trạng thái (Chương 7). Các ý tưởng được bổ sung bằng phép quét c_v 3D đầy đủ của sự nóng chảy bậc nhất (Δz̄, nhiệt ẩn ΔU, dịch chuyển đỉnh), một mô hình lò xo Einstein cho độ cứng phonon θ(mạng), và đường chuyển biến rắn(1)→rắn(2) qua Clausius–Clapeyron. Tiếp đó θ(mạng) được tinh chỉnh thành một DOS phonon kiểu Debye, rồi chuyển biến rắn(1)→rắn(2) còn được đưa vào như một đường trong bức tranh (τ, p) (hiệu chuẩn theo sắt), và bổ sung một nhánh Redlich–Kwong tách riêng bên cạnh mặc định van der Waals. Ở đó θ(mạng) được kiểm tra trên một mạng đầy đủ Lennard-Jones (FCC bền, BCC lý tưởng bất ổn định trượt), còn tôpô rắn(1)→rắn(2) được kiểm tra trên titan (trường hợp khẳng định thứ hai) và argon (đối chứng âm). Chương 8 xếp mô hình vào giữa mười một hình dung mô hình đã được thiết lập của nhiệt động lực học tĩnh. Phụ lục ghi lại các nguồn tài liệu (sắp theo tác giả) và sự hỗ trợ của AI (phần đóng góp đổi mới ước tính, chi phí ước tính và thời lượng của mỗi bước phát triển); các phần trợ giúp gom lại cách đọc, các điểm phê phán còn lại và triển vọng.

Luận điểm cốt lõi. Một thước đo liên kết hợp tác duy nhất z̄ ràng buộc ba dòng truyền thống vốn tách biệt — bức tranh khí mạng/Ising (chuyển pha hai trạng thái hợp tác), bức tranh thấm (tính liên thông → sôi) và bức tranh phonon Einstein/Debye (độ cứng → nóng chảy) — thành một bức tranh pha nhất quán, đầy đủ về mặt định tính, một cách có ý thức nằm bên trong một bức tranh lý tưởng kiểu van der Waals, trong khi tòa nhà lý thuyết có thể được suy ra một cách tự nhiên từ một hình ảnh động dễ hình dung. Không phải độ chính xác, mà tính hình ảnh mới là sức mạnh thực sự của mô hình và làm nên giá trị khoa học của nó.

Mục lục

• Tóm tắt tổng thể
• 1 — Các tiên đề với tên từ khóa (đã hợp nhất)
• 2 — Phân tích mâu thuẫn và khoảng trống
• 3 — Các mô tả bằng hình ảnh
• 3.1 Ba pha riêng lẻ
• 3.2 Các chuyển pha
• 3.3 Vì sao vận tốc âm trong chất lỏng giảm theo T
• 3.4 Chuyển động Brown hình thành thế nào — với khái niệm mixen được làm sắc
• 3.5 Vì sao phản ứng ở lớp biên mixen diễn ra nhanh hơn
• 3.6 Cửa sổ lỏng ba chiều
• 4 — Mô hình cơ học thống kê
• 4.1 Bậc tự do và biến số
• 4.2 Liên kết hai trạng thái hợp tác (kiểu Ising, tham số J)
• 4.3 Cân bằng năng lượng — với số hạng hợp tác
• 5 — Tổng trạng thái và nhiệt dung
• 5.1 Tổng trạng thái phân tích thành thừa số
• 5.2 Nhiệt dung theo từng pha (giải tích)
• 5.3 Đường c_v(τ) số **có** số hạng hợp tác — nhiệt ẩn được định lượng
• 5.4 Ngưỡng thấm liên thông như một kết quả (được khẳng định không đổi) - 5.4.a Hai ngưỡng khác nhau được tách bạch
• 5.5 Bậc của chuyển pha nóng chảy
• 5.6 Nghiên cứu kích thước hữu hạn của ngưỡng thấm
• 5.7 Phép quét c_v(τ) 3D đầy đủ của sự nóng chảy bậc nhất
• 6 — Tính hợp lý toán học của cách diễn giải 3D
• 6.1 Thay đổi bậc tự do quay 2D → 3D
• 6.2 Ảnh hưởng của các mạng tinh thể khác nhau — tính liên thông như dấu hiệu phân biệt
• 6.3 Chuyển pha bổ sung rắn(1) → rắn(2)
• 6.4 Suy ngẫm: dị thường âm của nước lỏng
• 6.5 Phương án khác: trường xoáy fractal thay cho sự quay mixen
• 6.6 Độ cứng phonon vi mô θ(mạng): một mô hình lò xo Einstein
• 6.7 Ước lượng toán học tính giả bền (2D tam giác ↔ vuông, 16b)
• 6.8 Tinh chỉnh kiểu Debye cho θ(mạng)
• 6.9 θ(mạng) từ một mạng đầy đủ Lennard-Jones
• 7 — Hiệu chuẩn và phương trình trạng thái
• 7.1 Phương pháp hiệu chuẩn (một tham số năng lượng ε) — đã điều chỉnh
• 7.2 Bảng hiệu chuẩn
• 7.3 Trục mật độ/thể tích: phương trình trạng thái van der Waals lõi cứng
• 7.4 Hiệu chuẩn l/g định lượng: a → ε, b → D²
• 7.5 Mặt (τ, v) chung với điểm ba
• 7.6 Nhiệt ẩn so với entanpi nóng chảy
• 7.7 Thống nhất các thang năng lượng: ε ↔ ε_lg
• 7.8 Điểm ba không tham số từ z̄_c
• 7.9 Tính toán kiểm tra 3D: cửa sổ lỏng sạch trên một thang ε
• 7.10 Phép quét thấm liên kết 3D thực (FCC/BCC) và độ nhạy mép cửa sổ
• 7.11 Phép quét τ 3D thực của động học liên kết: τ_b(3D) và đường nóng chảy trên một thang ε
• 7.12 Đường chuyển biến rắn(1)→rắn(2) trong mặt (τ, v)
• 7.13 Đường chuyển biến rắn(1)→rắn(2) trong giản đồ (τ, p), hiệu chuẩn theo sắt
• 7.14 Redlich–Kwong như một nhánh bức tranh lý tưởng tách riêng
• 7.15 rắn(1)→rắn(2) ở các đa hình khác: titan và đối chứng khí hiếm
• 8 — So sánh tài liệu: các hình dung mô hình của nhiệt động lực học tĩnh
• 8.1 Bảng tổng quan
• 8.2 Từng mô hình một
• 8.3 Điểm giống và khác so với Mô hình 3
• 9 — Phụ lục
• 9.1 Tài liệu tham khảo
• 9.2 Hỗ trợ của AI
• 9.3 Bảng thuật ngữ viết tắt và ký hiệu
• 9.4 Các kịch bản tính toán và độ phủ kiểm thử

<a id="kap1"></a>

1 — Các tiên đề với tên từ khóa (đã hợp nhất)

Tóm tắt. Bảy tiên đề được đặt tên bằng từ khóa xác định chất lý tưởng: số hạt được bảo toàn, chuyển động thuần cổ điển, một liên kết hai trạng thái hợp tác và số phối trí trung bình z̄ như tham số trật tự. Chúng nâng đỡ mọi suy dẫn về sau; trong bản gốc tiếng Đức, các từ khóa được đánh dấu bằng một thủy ấn đồng hình (khóa ở cuối văn bản).

Các tiên đề được làm chính xác:

1. BAY TỰ DO (tiên đề bổ sung). Giữa hai phân tử bay lướt qua nhau trong trạng thái không liên kết không có ma sát cũng không có lực tĩnh điện tác dụng. Động lượng và năng lượng chỉ thay đổi qua lực đẩy cứng LÕI CỨNG khi va chạm và qua việc tạo/đứt LIÊN KẾT. Tiên đề này chỉ có ý nghĩa từ góc nhìn cơ học cổ điển và cân bằng năng lượng: nó bảo đảm rằng không có kênh tiêu tán nào ở cấp phân tử — năng lượng chuyển động không bao giờ "mất đi", mà chỉ chuyển sang ba kho HỢP ÂM BA hoặc sang năng lượng liên kết. (Ma sát/độ nhớt vĩ mô nảy sinh trồi lên từ trò xếp hình trượt và VẾT NỨT, chứ không phải như một lực cơ bản vi mô.)
2. ĐA THỂ. Một chất lý tưởng gồm rất nhiều phân tử cùng loại; các tính chất vĩ mô được mô tả thống kê trong giới hạn nhiệt động (N → ∞).
3. LIÊN KẾT. Cứ hai phân tử lân cận về mặt hình học có thể mang một liên kết. Một liên kết nguyên vẹn hạ năng lượng xuống một lượng hữu hạn ε; việc bẻ gãy nó tốn đúng ε. Liên kết là một hệ hai trạng thái {nguyên vẹn, đứt}.
4. HỢP ÂM BA. Nội năng phân bố một cách cổ điển vào ba kho: (a) tịnh tiến của các phân tử, (b) dao động trong các liên kết nguyên vẹn (một lò xo điều hòa cho mỗi liên kết), (c) quay của các khối phân tử cứng (bậc tự do quay tập thể của các vùng liên thông).
5. LÕI CỨNG. Nếu khoảng cách phân tử xuống dưới giá trị tối thiểu σ, một lực đẩy cứng tác dụng (thể tích riêng). Cùng với LIÊN KẾT ta được một thế cặp hiệu dụng (vỏ cứng + hố liên kết hữu hạn).
6. BẢO TOÀN. Tổng động lượng và tổng năng lượng được bảo toàn. Một phân tử đơn lẻ trung hòa về năng lượng — nó không đệm năng lượng bên trong mình; năng lượng sống trong các liên kết (LIÊN KẾT) và trong các bậc tự do động học (HỢP ÂM BA).
7. VẾT NỨT. Các liên kết tan rã/hình thành theo nguyên tắc nứt/lành: càng thiếu nhiều liên kết lân cận → xác suất đứt càng cao; càng nhiều lân cận đã liên kết → xác suất lành càng cao. Động học mang tính hợp tác.
8. NHANH NHẸN. Đứt và lành diễn ra nhanh và ít hoạt hóa; hệ con liên kết ở mọi nhiệt độ đều nằm trong cân bằng cục bộ và được phép tách đoạn nhiệt khỏi chuyển động phân tử (cơ sở của tổng trạng thái phân tích thành thừa số).

Lập trường mô hình hóa. Sự phối trí hình học được giữ nguyên một cách lý tưởng hóa; chỉ có số phối trí z thay đổi do các chỗ đứt. Chuyển động nảy sinh như ở trò xếp hình trượt: các lỗ trống (từ thể tích riêng LÕI CỨNG và các chỗ đứt VẾT NỨT) di cư, các khối trượt theo. Hình ảnh "sai" một cách có ý thức về các tinh thể nhỏ vừa quay vừa di cư vẫn được giữ làm nền tảng trực quan.

<a id="kap2"></a>

2 — Phân tích mâu thuẫn và khoảng trống

Tóm tắt. Phân tích kiểm tra mô hình một cách hệ thống về các mâu thuẫn nội tại và khoảng trống, đồng thời ghi tiếp tình trạng xử lý; các điểm yếu được nhận ra sẽ được nêu tên và trong các chương sau hoặc được giải quyết hoặc được đánh dấu một cách trung thực như một giới hạn của mô hình.

Danh sách các điểm phê phán

1. Đếm trùng chuyển động tịnh tiến (tiên đề 3a so với 3c). Bản thô nêu tịnh tiến hai lần (tịnh tiến phân tử và tịnh tiến khối). Không thể có hai kho tịnh tiến độc lập. Khép lại: tịnh tiến của khối chỉ là một mode tập thể của các tịnh tiến phân tử đã được đếm; điều thực sự mới chỉ là bậc tự do quay của khối. Trong cân bằng bậc tự do (Ch. 4) tịnh tiến được đếm đúng một lần.
2. "Trung hòa về năng lượng" (tiên đề 5) so với "liên kết hấp thụ năng lượng" (tiên đề 2/3b). Mâu thuẫn bề ngoài. Khép lại: tính trung hòa liên quan tới phân tử như một vật thể; năng lượng nằm giữa các phân tử (trong LIÊN KẾT) và trong các mode động học tập thể, không bao giờ trong một bộ đệm nội phân tử. Phân tử ≠ liên kết — mâu thuẫn nhờ đó được hóa giải.
3. Thiếu lực hút (tiên đề 4 chỉ nêu lực đẩy). Không có lực hút thì không có pha ngưng tụ. Khép lại: lực hút do LIÊN KẾT cung cấp (hố sâu ε), lực đẩy do LÕI CỨNG (vỏ cứng tại σ). Cùng nhau tạo thành một thế cặp hiệu dụng.
4. Phản hồi không có số hạng đối trọng (tiên đề 6). "Càng thiếu lân cận → càng nhiều đứt" tự nó là một phản hồi dương (vết nứt lan đến sụp đổ hoàn toàn). Khép lại: số hạng lành (càng nhiều lân cận → càng nhiều lành) cộng với thừa số Boltzmann của năng lượng liên kết hữu hạn ε (LIÊN KẾT) cân bằng VẾT NỨT và tạo ra một chuyển pha xác định thay vì một sự bùng nổ.
5. Sự tách thang thời gian không được nói ra (tiên đề 7). "Nhanh và ít hoạt hóa" chỉ dùng được khi từ đó rút ra sự tách đoạn nhiệt giữa động học liên kết và động học tịnh tiến. Khép lại: được ghi rõ như tiền đề của phép phân tích thừa số Z = Z_trans·Z_vib·Z_rot·Z_bond.
6. Thiếu hình học và thiếu cơ chế chuyển động. Các tiên đề không nêu mạng phối trí, cũng không nêu điều gì làm cho sự chảy trở nên khả dĩ. Khép lại (như một giả thiết bổ sung được ghi rõ): một hình học phối trí quy chiếu cố định với sự lấp đầy liên kết thay đổi; làm trường hợp 2D kinh điển là phối trí tam giác (z_max = 6). Tính linh động qua cơ chế lỗ trống kiểu trò xếp hình trượt. Hai điểm này là những bổ sung thực sự duy nhất vượt ra ngoài bảy tiên đề và được đánh dấu ở Chương 6 như một câu hỏi chiến lược.
7. Khoảng trống tiêu tán tiềm ẩn. Bản 07 để ngỏ việc các phân tử không liên kết có mất năng lượng khi bay lướt qua nhau hay không. Khép lại nhờ BAY TỰ DO: không có ma sát/tĩnh điện vi mô; mọi dòng năng lượng và động lượng chỉ chạy qua các va chạm LÕI CỨNG và động học LIÊN KẾT. Nhờ đó cân bằng năng lượng (Ch. 4.3) khép kín và thống kê Maxwell–Boltzmann trong trường hợp giới hạn khí là chính xác (khí lý tưởng tự do giữa các va chạm đàn hồi).
8. Sự bất nhất về năng lượng của tham số trật tự. Bản 08 tuy có dùng số hạng hợp tác −(K/2)x² trong năng lượng tự do g(x) (để tạo x(τ)), nhưng không dùng trong nội năng U/N. Điều đó bất nhất về mặt nhiệt động: cùng một số hạng trường phân tử phải xuất hiện ở cả hai. Khép lại (Ch. 4.3): −(K/2)x² được đưa vào U/N; nhiệt ẩn nhờ đó tăng lên rõ rệt (Ch. 5.3) mà không làm đổi các bình nguyên rắn và khí.
9. Hiệu chuẩn l/g chỉ mang tính định tính. Bản 09 mới chỉ đưa trục mật độ (phương trình trạng thái van der Waals) ở dạng rút gọn. Khép lại (Ch. 7.4): qua điểm tới hạn (T_c, p_c) đường sôi tuyệt đối được tính và T_b(p_atm) được đối chiếu với các giá trị đo; độ lệch hệ thống (quy tắc Guldberg) được nêu tên.
10. Các trục tách rời không có điểm ba. Trục nóng chảy và trục sôi là những phép tính riêng biệt. Khép lại (Ch. 7.5): cả hai được hợp nhất trong một mặt phẳng (τ, v), điểm ba được dựng và điểm tới hạn được đánh dấu.

<a id="kap3"></a>

3 — Các mô tả bằng hình ảnh

Tóm tắt. Chương 3 giải thích một cách trực quan (với các dải truyện tranh tĩnh) các trạng thái rắn, lỏng, khí cùng các chuyển pha: ba pha riêng lẻ (3.1), các chuyển pha kết tinh lại/nóng chảy/sôi/thăng hoa cùng điểm ba và điểm tới hạn (3.2), dị thường âm (3.3), chuyển động Brown và khái niệm mixen được làm sắc (3.4–3.5), cũng như cửa sổ lỏng ba chiều (3.6). Dựa trên các mô tả đó, các sơ đồ (dải truyện tranh) đã được tạo ra nhằm giúp hình dung trước con mắt tâm trí các quan hệ phân tử động trong các pha khác nhau. Chính khả năng hình dung này mới là hạt nhân thực sự của mô hình. Nó trở nên có giá trị về mặt định tính và khoa học khi qua cơ học thống kê có thể chỉ ra một cách hợp lý rằng mô hình làm cho các xu hướng/quan sát thực nghiệm định tính trở nên tiên đoán được ngay cả khi không cần tới toán học lớn.

<a id="kap3-1"></a>

3.1 Ba pha riêng lẻ

Chất rắn. Hầu như mọi lân cận đều liên kết với nhau qua các liên kết phân tử (LIÊN KẾT), giá trị trung bình của số phối trí của các phân tử z nằm gần cực đại hình học. Mạng phân tử là một khối cứng duy nhất không có lỗ trống (LÕI CỨNG xếp chặt) — một trò xếp hình trượt nội tại như ở chất lỏng là không thể có trong phiên bản nén chặt. Chuyển động duy nhất: rung quanh vị trí cân bằng (dao động liên kết, HỢP ÂM BA-b), do các liên kết phân tử với các lân cận quy định. Động lượng/năng lượng dao động đàn hồi trong mạng (BẢO TOÀN, BAY TỰ DO: không hao tổn). Ở trường hợp giới hạn (nhiệt độ rất thấp, xếp chặt lý tưởng) hiện ra trạng thái lý tưởng của một tinh thể lý tưởng.

Chất lỏng. VẾT NỨT hợp tác đã thấm qua mạng: liên kết thiếu thành từng mảng, giá trị trung bình số phối trí của các phân tử z đã giảm, khối cứng bên trong đã vỡ thành nhiều khối nhỏ hơn. Từ khối nguyên liền đã thành một đống đá vụn. Các khối này tịnh tiến và quay (HỢP ÂM BA-c). Giữa chúng đã mở ra các lỗ trống, nhờ đó chúng trượt lướt qua nhau — trò xếp hình trượt của chuyển động Brown đang hoạt động. Hình học phối trí cục bộ vẫn được giữ nguyên một cách lý tưởng hóa; chỉ thay đổi ở chỗ bao nhiêu liên kết thực sự tồn tại và được hình thành.

Chất khí. VẾT NỨT đã thắng: liên kết đứt gần như ở mọi nơi (z → 0). Trong trường hợp cực đoan không phân tử nào còn một liên kết phân tử tới một phân tử khác. Các khối đã vỡ thành từng phân tử riêng lẻ. Chúng bay tự do, chỉ ghép cặp qua các va chạm LÕI CỨNG (giữa chừng là BAY TỰ DO); các vận tốc tuân theo Maxwell–Boltzmann (BẢO TOÀN). Dao động liên kết và quay khối đã tắt — chỉ còn tịnh tiến thuần túy (HỢP ÂM BA-a). Ở trường hợp giới hạn (nhiệt độ cao, áp suất thấp) hiện ra trạng thái lý tưởng của một khí lý tưởng.

<a id="kap3-2"></a>

3.2 Các chuyển pha

Kết tinh lại (rắn(1)/rắn(2)). Vẫn còn dưới điểm nóng chảy, chất rắn có thể đổi mạng của nó mà không mất độ cứng: hai cách xếp chặt lý tưởng với số phối trí khác nhau (FCC z_max = 12 ↔ BCC z_max = 8) có các entanpi tự do μ(τ, p) khác nhau; tại μ_FCC = μ_BCC hình học phối trí nhảy đổi (chứ không chỉ độ lấp đầy z̄ như khi nóng chảy). Cái nhìn thấy được là một bậc c_v nhỏ (Δc_v = biến thiên z_max/4 = 1,0 k_B) cộng với một nhiệt chuyển biến ẩn yếu. Khẳng định của AI: điều đó chỉ hình dung được trong 3D (tính liên thông rời rạc, lượng tử hóa; Ch. 6.3). Trong trường hợp hai chiều, sự kết tinh lại thực sự, tức chuyển pha giữa hai pha rắn, dễ hình dung hơn. Trên mặt phẳng chỉ có ba cách lát đều không có lỗ hổng. Người ta có thể phủ kín hoàn toàn (lát) một mặt phẳng bằng tam giác đều, bằng hình vuông hoặc bằng lục giác đều. Nếu trong các mạng đó các cạnh biểu diễn liên kết phân tử và các đỉnh biểu diễn phân tử, thì một chuyển pha rắn(tam giác)→rắn(vuông) là hình dung được. Khi đó phân tử trong lát tam giác còn có sáu liên kết tới các phân tử lân cận, trong khi phân tử trong mạng vuông chỉ còn bốn.

Nóng chảy (r/l). Khi nhiệt độ tăng, mạng rách ở đủ nhiều chỗ để khối cứng tổng thể mất độ cứng (VẾT NỨT hợp tác xuống dưới ngưỡng cứng vững Maxwell z̄ = 2d — mạng trở nên mềm oặt). Hình thành các mixen nhiệt động, tức là các vùng không còn nối với các vùng lân cận lớn hơn qua liên kết. Khối nguyên liền từ góc nhìn phân tử trở thành đống đá vụn và có khả năng chảy hạn chế. Các lỗ trống đầu tiên cho phép dịch chuyển các khối phân tử tương tự trò xếp hình trượt. Các liên kết còn thiếu với vùng lân cận cho phép sự quay nội tại trong các mixen. Sự quay mixen (HỢP ÂM BA-c) bắt đầu. Việc bẻ gãy nhiều liên kết (LIÊN KẾT) ngốn nhiệt ẩn.

Sôi (l/k). Đun nóng thêm khiến chính các lỗ trống thấm qua nhau: chúng nối lại thành khoảng trống liên tục, trong đó các khối tan ra và từng phân tử bay đi. Sôi là chuyển pha từ "lỗ trống trong vật chất" sang "vật chất trong lỗ trống".

Thăng hoa (r/k). Ở áp suất thấp, vùng trung gian lỏng không có: mạng nhả thẳng từng phân tử vào khoảng trống mà không xảy ra sự tạo khối như ở các mixen. Ở đây có thể suy ngẫm xem các khía cạnh hình học (phân tử dài, phân tử đan xen nhau), các khía cạnh tĩnh điện (sự định hướng do mômen từ hay mômen lưỡng cực), các loại liên kết phân tử khác nhau (van der Waals, tương tác giữa các hệ cộng hưởng, liên kết hiđro) có thể ảnh hưởng thế nào tới việc hình thành hay không hình thành các mixen nhiệt động.

Điểm ba. Đúng tại đó hai ngưỡng trùng nhau: độ cứng vững Maxwell (rắn ↔ lỏng, mất độ cứng) và sự thấm của lỗ trống (lỏng ↔ khí) xảy ra tại cùng (τ, p). Cả ba hình ảnh (khối cứng [nguyên liền], các khối trượt [đống đá vụn], phân tử tự do [đám bụi]) cùng tồn tại.

Điểm tới hạn. Khi áp suất tăng, bước nhảy lỏng→khí ngày càng nhỏ; tỉ lệ lỗ trống biến đổi ngày càng liên tục. Tại điểm tới hạn, sự khác biệt biến mất hoàn toàn — "các khối đặc với ít khoảng trống" và "nhiều khoảng trống với ít phân tử" chuyển tiếp liên tục sang nhau. ("Từ góc nhìn của chất lỏng thì cốc vơi một nửa, từ góc nhìn của chất khí thì đầy một nửa.")

Kết tinh lại trong bức tranh 2D. Truyện tranh nay trình bày sự kết tinh lại như trường hợp 2D xếp chồng lên nhau: bên dưới là mạng tam giác xếp chặt (phối trí z = 6, lý tưởng 2D có năng lượng thấp nhất), bên trên là mạng vuông lỏng lẻo (z = 4). Chuyển pha tam giác → vuông làm tăng thể tích (công chống lại áp suất) và hạ số liên kết. Ta thấy số phối trí thay đổi, nhưng độ cứng của một mạng tinh thể vẫn còn.

Tiên đoán định tính về chuyển đổi năng lượng và nhiệt dung. Chỉ từ tính liên thông z có thể ước lượng cả hai chỉ số: (i) năng lượng liên kết trên mỗi hạt là U/N ≈ −(ε/2)·z (mỗi liên kết tính một nửa), tức năng lượng càng thấp khi z càng cao; (ii) bước nhảy c_v khi đổi mạng tuân theo quy tắc mô hình Δc_v = biến thiên z_max/4. Với trường hợp 2D tam giác (z = 6) ↔ vuông (z = 4), điều đó cho hiệu năng lượng ΔU ≈ ε trên mỗi hạt và bậc c_v Δc_v = (6 − 4)/4 = 0,5 k_B — bằng một nửa giá trị 3D (1,0 k_B). Câu ghi nhớ: càng đặc = năng lượng càng thấp, nhưng chuyển pha càng đắt (nhiệt ẩn lớn hơn, bậc c_v lớn hơn), vì khi nới rộng có nhiều liên kết hơn và nhiều bậc tự do bị đóng băng hơn được giải phóng.

Cấu trúc mạng có quyết định tính giả bền không? Một phần. Cấu trúc mạng ấn định thứ tự năng lượng — dạng nào là trạng thái cơ bản (tam giác, μ thấp nhất) và dạng nào bị kích thích (vuông). Nhưng dạng lỏng lẻo có sống sót giả bền ở nhiệt độ thấp hay không thì không phải một mình nó quyết định: trong mô hình lý tưởng nghiêm ngặt (không có rào), chất tại mỗi (τ, p) lập tức nhận pha có μ thấp nhất — ở đó không có tính giả bền. Nó chỉ nảy sinh nhờ một năng lượng hoạt hóa hữu hạn (rào tạo mầm) kìm hãm sự tái cấu trúc về mặt động học. Như vậy hình học mạng cung cấp thứ bậc nhiệt động, còn rào cung cấp sự dai dẳng động học — chính thành phần cũng giải thích pha quá bão hòa (xem dưới; định lượng ở Ch. 6.7).

Pha quá bão hòa — năng lượng hoạt hóa. Hình bổ sung (3.2, hàng 2) cho thấy một pha khí quá bão hòa bên trên một chất lỏng: bên cạnh các phân tử đơn lẻ là các nhóm nhỏ gồm 2–3 hạt kết dính (các thanh màu xám) bay như một khối. Mô hình giải thích sự quá bão hòa bằng một năng lượng hoạt hóa đáng kể cho các liên kết phân tử: các cụm nhỏ (2–3) có tỉ lệ biên trên diện tích bất lợi, các liên kết biên của chúng không đầy đủ — chỉ một mầm tới hạn mới tự phát lớn lên tiếp. Chừng nào chưa có mầm siêu tới hạn hình thành, khí vẫn quá bão hòa, dù μ_khí > μ_lỏng đã cho phép ngưng tụ. Chính rào này thiếu trong bức tranh cân bằng lý tưởng và ở đây được nêu như một đóng góp độc lập, "nằm ngoài bức tranh lý tưởng" (ước lượng định lượng ở Ch. 6.7).

<a id="kap3-3"></a>

3.3 Vì sao vận tốc âm trong chất lỏng giảm theo T

Âm là một sóng đàn hồi: v = √(K_eff/ρ). K_eff chỉ sinh ra từ các liên kết nguyên vẹn (LIÊN KẾT) — chúng chính là các lò xo. Khi T tăng, ngày càng nhiều liên kết đứt một cách hợp tác (VẾT NỨT); mạng lò xo thưa đi và K_eff giảm dốc hơn mật độ ρ. Sự đứt chi phối tử số chứ không phải mẫu số — vì thế v giảm. Diễn giải theo phân tử: sự truyền âm qua các liên kết phân tử nhanh hơn sự truyền âm cơ học qua các va chạm cơ học. Càng nhiều vết nứt xuất hiện trong một chất lỏng khi nhiệt độ tăng, âm càng thường xuyên phải được truyền qua các va chạm cơ học thay vì qua các liên kết phân tử. Tương ứng, vận tốc âm chậm dần khi nhiệt độ tăng. (Dị thường của nước — một cực đại âm ở ~74 °C — được bàn ở Ch. 6.4.)

<a id="kap3-4"></a>

3.4 Chuyển động Brown hình thành thế nào — với khái niệm mixen được làm sắc

Theo nghĩa nhiệt động thống kê, ở đây một mixen không phải là tập hợp chất hoạt động bề mặt, mà là một khối được VẾT NỨT viền quanh: các liên kết thiếu (các phi-liên-kết ở biên khối) tạo thành "lớp vỏ" quanh một tinh thể nhỏ bị bao kín (lõi liên kết nguyên vẹn). Lõi này có thể thực hiện một chuyển động quay hợp tác nội tại (HỢP ÂM BA-c) — khối quay như một thể thống nhất, trong khi lớp vỏ biên của nó tách nó khỏi các lân cận.

Chuyển động Brown nay hình thành như sau: qua sự trao đổi các liên kết phân tử (VẾT NỨT, NHANH NHẸN), từng phân tử có thể rời lớp vỏ của một mixen và chuyển sang một mixen khác. Mỗi lần chuyển như vậy truyền động lượng (BẢO TOÀN; nhờ BAY TỰ DO nên không hao tổn) và làm thay đổi đột ngột hình dạng cùng mômen động lượng của cả hai mixen tham gia. Một hạt thử lớn (hạt bụi) được nhiều sự tái cấu trúc mixen không tương quan như thế bao quanh; tổng các cú hích động lượng trung bình bằng không, nhưng có phương sai hữu hạn → bước đi ngẫu nhiên. T càng cao thì các sự kiện VẾT NỨT và sự đổi mixen càng thường xuyên, sự rung càng mạnh. Chuyển động Brown do vậy là cái bóng vĩ mô hay hệ quả của quá trình trao đổi mixen vi mô.

<a id="kap3-5"></a>

3.5 Vì sao phản ứng ở lớp biên mixen diễn ra nhanh hơn

Lời rào đón ngượng ngùng. Lập luận sau đây là một suy dẫn từ bức tranh và theo chỗ tôi biết cho tới nay chưa được chứng minh hay khảo sát bằng thực nghiệm, và tôi phải thú nhận rằng tôi đã không chịu khó tra cứu kỹ về điều đó. Điểm tựa duy nhất cho tính mới chỉ là các hệ AI không nhắc tới hệ quả này. Luận đề này là một ví dụ cho thấy các hình dung bằng hình ảnh có thể dẫn tới những câu hỏi mới ra sao.

Trong các dung môi hòa tan vừa phải, các chất phản ứng bám ưu tiên vào lớp vỏ biên của các mixen — trong mô hình chính là các biên nứt/biên khối (VẾT NỨT). Hai hiệu ứng làm tăng tốc phản ứng: (i) giảm chiều — các chất phản ứng không còn tìm nhau trong toàn thể tích, mà dọc theo một mặt biên gần như hai chiều; một cuộc tìm kiếm 2D về mặt thống kê gặp nhau nhanh hơn một cuộc tìm kiếm 3D. (ii) Sự cô đặc — lớp vỏ biên gom các chất phản ứng lại, mật độ cục bộ của chúng tăng. Theo hình ảnh, các biên nứt tác dụng như những thanh dẫn hướng đưa các đối tác đến với nhau.

Hiệu ứng chiều có một hạt nhân chặt chẽ: bước đi ngẫu nhiên (Pólya). Vì sao một cuộc tìm kiếm hai chiều về mặt thống kê gặp nhau tốt hơn một cuộc tìm kiếm ba chiều thì không chỉ suy ra từ trực giác, mà từ một định lý của lý thuyết xác suất. George Pólya (1887–1985) năm 1921 đã chứng minh cho bước đi ngẫu nhiên ("random walk") trên một mạng: trong một và hai chiều, người đi với xác suất 1 (tức chắc chắn) sẽ trở lại mỗi điểm vô hạn lần, trên mạng phẳng ℤ² cũng như trên đường thẳng ℤ. Từ ba chiều trở đi điều đó không còn đúng: xác suất trở lại thực sự nhỏ hơn 1 (trên mạng lập phương đơn giản ℤ³ chỉ khoảng 0,34). Một cách trực quan, theo hình ảnh của Pólya: "Một người say rượu chắc chắn có lúc sẽ về được nhà, còn một con cá say thì không." Áp dụng cho hai chất phản ứng phải gặp nhau nhiều lần để có một va chạm thành công, điều đó nghĩa là: trong một dung môi hòa tan rất tốt với bước đi ngẫu nhiên ba chiều tự do, các đối tác một khi đã tách ra thì với xác suất hữu hạn sẽ xa nhau mãi mãi — các lần thử phản ứng lặp lại thất bại, và hiệu suất vẫn thấp thuần vì lý do hình học-thống kê này. Ngược lại, ở một lớp biên mixen (gần như) hai chiều thì sự gặp lại là chắc chắn (xác suất 1): các đối tác luôn tìm lại được nhau cho tới khi phản ứng thành công. Như vậy lớp biên không chỉ cô đặc, mà đồng thời cung cấp chiều thuận lợi hơn cho các cuộc gặp lại. Khẳng định (được làm sắc một cách có ý thức) là: hiệu suất trong các dung môi hòa tan rất tốt có xu hướng kém hơn vì lý do bước đi ngẫu nhiên.

Nguồn: G. Pólya (1921), Über eine Aufgabe der Wahrscheinlichkeitsrechnung betreffend die Irrfahrt im Straßennetz, Math. Ann. 84, 149–160 — ở đây dẫn theo M. Eisermann, Höhere Mathematik 3 (Univ. Stuttgart), bài tập W-2×2. (Bản gốc truy cập 5.9.2026).

<a id="kap3-6"></a>

3.6 Cửa sổ lỏng ba chiều

Khác với các hình ở 3.1–3.5, hình này không cho thấy một cảnh phân tử, mà cho thấy kết quả tính toán từ Ch. 7.9/7.10 như một sơ đồ trên một trục τ (đơn vị năng lượng liên kết ε). Mục tiêu là làm thấy được vì sao một mạng 2D đơn lẻ hầu như không cho phép một vùng cho pha lỏng, trong khi ở mạng 3D điều đó là khả dĩ không vấn đề gì, dù cho cả hai trường hợp đều áp dụng cùng một quy tắc cứng vững (z̄_c = 2d).

Hình ① — 2D (không có cửa sổ). Trên trục τ, τ_m = 1,20 và τ_triple = 1,16 gần như trùng nhau (τ_triple dưới τ_m), τ_c = 1,50 sát ngay trên: vùng cho pha lỏng bị ép thành một khe — rắn gần như chuyển thẳng sang khí (gần thăng hoa).

Hình ② — 3D (cửa sổ sạch). Cùng trục ấy, nhưng được giãn ra: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00. Giữa τ_triple và τ_c là miền lỏng độc lập (rộng về phía dưới, τ_triple − τ_m = 1,63; hẹp về phía trên, τ_c − τ_triple = 0,18). Mép trên hẹp được đánh dấu như vậy.

Hình ③ — nguyên nhân. z̄_c = 2d tăng theo chiều (4 → 6), do đó v_solid giảm (0,50 → 0,667 trong thể tích rút gọn, tức chất rắn trở nên "lỏng lẻo hơn" so với chất lỏng), và τ_c = z_max/4 tăng gấp đôi (1,50 → 3,00). Điểm tới hạn dịch xa lên trên điểm nóng chảy. Chỉ nhờ đó mới có chỗ cho pha lỏng. Đó là một hiệu ứng chiều thuần túy.

<a id="kap4"></a>

4 — Mô hình cơ học thống kê (không đổi so với bản 09)

Tóm tắt

Chương 4 dựng khung cơ học thống kê của mô hình trong 2D (mạng tam giác, z_max = 6):

• 4.1 Bậc tự do và biến số: đại lượng trật tự là tỉ lệ liên kết x ∈ [0,1] với số phối trí trung bình z̄ = z_max·x; mỗi khối ≥ 2 thêm một bậc tự do quay, nhiệt độ τ được đếm theo đơn vị năng lượng liên kết ε.
• 4.2 Liên kết hai trạng thái hợp tác: phương trình liên kết tự hợp x = 1/(1 + G·exp[−(ε + K·x)/τ]) (K = J·z_max, G = z_max) suy ra như điều kiện dừng của năng lượng tự do Bragg–Williams; cực tiểu toàn cục của nó cho nhánh vật lý, còn nhánh ấy chạy liên tục hay có bước nhảy (bậc nhất) thì Ch. 5.5 quyết định.
• 4.3 Cân bằng năng lượng với số hạng hợp tác: U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x² (2D) hợp nhất tịnh tiến/dao động liên kết/quay khối, năng lượng cấu hình liên kết và năng lượng trường phân tử hợp tác mới −(K/2)x². Nhiệt dung tương ứng c_v nhận thêm số hạng dương −K·x·ẋ, số hạng này làm sắc đỉnh nóng chảy, nhưng để nguyên bình nguyên rắn và nhánh khí.

<a id="kap4-1"></a>

4.1 Bậc tự do và biến số

• N phân tử trên một hình học phối trí quy chiếu cố định. Trường hợp 2D kinh điển: mạng tam giác, z_max = 6, một bậc tự do quay cho mỗi khối (ĐA THỂ, LIÊN KẾT). Chương này áp dụng rõ ràng cho cách diễn giải 2D; sự tổng quát hóa 3D nằm ở Ch. 6.
• Biến liên kết s_ij ∈ {0, 1} cho mỗi cặp lân cận hình học (1 = nguyên vẹn).
• Số phối trí z_i = Σ_j s_ij; giá trị trung bình z̄ = 2 N_bond / N — tham số trật tự. Ta có z̄ = z_max · x với tỉ lệ liên kết x = p_b ∈ [0, 1].
• Đơn vị năng lượng τ ≡ k_B T (cổ điển, ½τ cho mỗi bậc tự do bậc hai). Mọi năng lượng tính theo đơn vị năng lượng liên kết ε.

<a id="kap4-2"></a>

4.2 Liên kết hai trạng thái hợp tác (kiểu Ising, tham số J)

Một liên kết cô lập (năng lượng −ε khi nguyên vẹn, 0 khi đứt) trong cân bằng Maxwell–Boltzmann có p_b⁰(τ) = 1/(1 + e^(−ε/τ)). VẾT NỨT hợp tác (tiên đề 6) dịch chuyển độ bền liên kết hiệu dụng đi một trường phân tử tỉ lệ với tỉ lệ liên kết x. Trong gần đúng trường phân tử (Ising), phương trình liên kết tự hợp trở thành:

x = 1 / (1 + G · exp[−(ε + K·x) / τ]), K = J·z_max, z̄ = z_max·x.

Trong đó:

• K = J·z_max là cường độ ghép cặp hợp tác (phản hồi VẾT NỨT: càng nhiều liên kết lân cận → độ bền liên kết hiệu dụng càng cao). J > 0 làm sắc chuyển pha; trên một J tới hạn thì z̄(τ) trở nên không liên tục — đó là nhiệt ẩn của sự nóng chảy. Trường hợp hai trạng thái thuần túy là trường hợp giới hạn J = 0 (không có trường phân tử, chuyển tiếp liên tục).
• G > 1 là độ suy biến của trạng thái đứt: một phân tử được giải phóng thu được không gian pha (tự do tịnh tiến). G thúc đẩy sự đứt liên kết ở τ cao và bảo đảm rằng x → 1/(1+G) thay vì mắc lại ở ½ (trường hợp giới hạn khí). Kinh điển G = z_max.

Phương trình này suy ra như điều kiện dừng của một năng lượng tự do Bragg–Williams trên mỗi phân tử (n_slots = z_max/2 chỗ liên kết):

g(x)/n_slots = −ε·x − (K/2)·x² + τ·[x ln x + (1−x) ln(1−x)] − τ·ln(G)·(1−x).

Số hạng −(K/2)x² là ghép cặp hợp tác (sắt từ), số hạng cuối là độ suy biến của sự đứt. Cực tiểu toàn cục của g theo x cho nhánh vật lý; nơi hai cực tiểu sâu bằng nhau là chuyển pha bậc nhất. NHANH NHẸN (tiên đề 7) biện minh cho việc x(τ) ở mọi nhiệt độ đều nhận giá trị cân bằng.

<a id="kap4-3"></a>

4.3 Cân bằng năng lượng — với số hạng hợp tác

Số hạng trường phân tử −(K/2)x² là một năng lượng thực (không phụ thuộc τ) và phải xuất hiện trong U/N. Cân bằng nhất quán do đó là (2D):

U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x².

• τ·[…]: tịnh tiến (1) + dao động liên kết (z̄/4) + quay khối (½ n_blk) — HỢP ÂM BA.
• −(ε/2)·z̄: năng lượng cấu hình liên kết — LIÊN KẾT.
• −(K/2)·x²: năng lượng trường phân tử hợp tác — VẾT NỨT. Mới trong U/N.

Nhiệt dung tương ứng c_v = d(U/N)/dτ (theo đơn vị k_B) nhận thêm một đóng góp:

c_v = [1 + z̄/4 + ½ n_blk] + τ·[¼ ż̄ + ½ ṅ_blk] − (ε/2)·ż̄ − K·x·ẋ, (ẋ = dx/dτ)

trong đó số hạng mới −K·x·ẋ mang giá trị dương (ẋ < 0) và tác dụng đúng ở nơi nhiều liên kết đứt cùng lúc. Nó làm sắc đỉnh nóng chảy mà không làm đổi bình nguyên rắn (ẋ ≈ 0) hay nhánh khí.

<a id="kap5"></a>

5 — Tổng trạng thái và nhiệt dung

Tóm tắt

Chương 5 tính các tổng trạng thái và nhiệt dung, đồng thời cung cấp dữ liệu mô hình cho việc hiệu chuẩn (Ch. 7):

• 5.1 Tổng trạng thái phân tích thành thừa số: nhờ sự tách thang thời gian mà Z = Z_trans·Z_vib·Z_rot·Z_bond.
• 5.2 Nhiệt dung theo từng pha: c_v = d(U/N)/dτ; số hạng dương −(ε/2)(dz̄/dτ) tạo ra đỉnh kiểu Schottky như tương tự của nhiệt ẩn.
• 5.3 Đường c_v(τ) số: đỉnh c_v 5,51 tại τ_m = 1,20, bình nguyên rắn 2,50, nhánh khí → 1; nhiệt ẩn L ≈ 1,41 ε/phân tử (có số hạng hợp tác) gấp ba so với 0,48 ε (không có).
• 5.4 / 5.4a Hai ngưỡng tách biệt: thấm tính liên thông z̄_c = z_max·p_c ≈ 2,08 (tại τ_b ≈ 1,85) đánh dấu sự tan rã cụm (gần sôi) và là một kết quả; phải tách nó rành mạch khỏi độ cứng vững Maxwell tổng quát z̄_c = 2d (mật độ nóng chảy, mất độ cứng).
• 5.5 Bậc của chuyển pha nóng chảy: tiêu chuẩn ln G > 2 + 4ε/K; 2D (G = 6) liên tục với mọi J, 3D (G = 12) bậc nhất từ K > ~8,2 (minh hoạ: Δz̄ ≈ 5,3, ΔU ≈ 1,75 ε).
• 5.6 Nghiên cứu kích thước hữu hạn: p*(L) hội tụ đơn điệu từ dưới về p_c = 0,3473; z̄_c ≈ 2,08 bền vững với kích thước hữu hạn.
• 5.7 Quét c_v 3D đầy đủ: từ K > K_tới hạn ≈ 8,25 (J > 0,68) mở ra bước nhảy thật Δz̄ cùng nhiệt ẩn ΔU; đỉnh nóng chảy dịch theo J lên các τ cao hơn (dịch chuyển dương +0,5 … +2,2 so với đỉnh 2D liên tục).

<a id="kap5-1"></a>

5.1 Tổng trạng thái phân tích thành thừa số

Nhờ sự tách thang thời gian (NHANH NHẸN), tổng trạng thái chính tắc phân tích thành thừa số:

Z = Z_trans · Z_vib · Z_rot · Z_bond.

• Z_trans = (1/N!)·(A/λ_T²)^N, λ_T = h/√(2π m τ) (tịnh tiến 2D cổ điển lý tưởng; ĐA THỂ, BẢO TOÀN, BAY TỰ DO).
• Z_vib = (τ/ħω₀)^(N_bond) (dao động tử điều hoà cổ điển cho mỗi liên kết nguyên vẹn; HỢP ÂM BA-b).
• Z_rot = ∏_(khối≥2) (2π I_b τ)^(1/2)/ħ (vật quay cứng 2D cổ điển cho mỗi khối; HỢP ÂM BA-c).
• Z_bond = Σ_{s} exp[(ε/τ)·Σ s_ij + (K/τ)·số hạng hợp tác] (tổng liên kết kiểu Ising; LIÊN KẾT, VẾT NỨT). Tính theo trường phân tử cho x(τ), z̄(τ) như ở 4.2.

<a id="kap5-2"></a>

5.2 Nhiệt dung theo từng pha (giải tích)

Từ U/N (Ch. 4.3), c_v = d(U/N)/dτ (theo đơn vị k_B):

c_v = [1 + z̄/4 + ½ n_blk] + τ·[¼ (dz̄/dτ) + ½ (dn_blk/dτ)] − (ε/2)·(dz̄/dτ).

Số hạng −(ε/2)·(dz̄/dτ) mang giá trị dương (dz̄/dτ < 0: đun nóng bẻ gãy liên kết) và chi phối ở nơi nhiều liên kết đứt cùng lúc — đó là đỉnh kiểu Schottky, tương tự của nhiệt ẩn.

<a id="kap5-3"></a>

5.3 Đường c_v(τ) số có số hạng hợp tác — nhiệt ẩn được định lượng

Liên kết tự hợp (cực tiểu toàn cục), n_blk(τ) từ Monte Carlo thấm liên kết (Union-Find, L = 100, 5 lần thực hiện), c_v = dU/dτ (sai phân trung tâm), U/N bao gồm −(K/2)x² (cv_scan_09.py, ε = 1, J = 0,55, G = z_max = 6). Một số điểm mốc chọn lọc:

Kết quả (số). Trật tự ba pha vẫn được khẳng định và nay sắc nét hơn:

c_v(rắn) ≈ 2,50 < c_v(lỏng, cực đại) ≈ 5,51 > c_v(nhánh khí) ↓ → 1.

Nhiệt ẩn (diện tích đỉnh trên đường nền). Lấy tích phân (c_v − đường nền) trên cửa sổ nóng chảy τ ∈ [0,70 ; 1,70] (đường nền = đường thẳng qua c_v(0,70) và c_v(1,70)), ta được

L ≈ 1,41 ε/phân tử (có số hạng hợp tác), L ≈ 0,48 ε/phân tử (không có số hạng hợp tác, cân bằng cũ), ΔL ≈ +0,92 ε/phân tử.

Số hạng hợp tác nâng đỉnh lên khoảng 1,7 lần (3,25 → 5,51) và làm gấp ba nhiệt ẩn. Về mặt vật lý: việc bẻ gãy đồng thời nhiều liên kết lân cận (VẾT NỨT) không chỉ tích trữ năng lượng cấu hình (−ε/2·z̄) mà cả năng lượng trường phân tử (−K/2·x²) — cả hai được giải phóng khi nóng chảy.

<a id="kap5-4"></a>

5.4 Ngưỡng thấm liên thông như một kết quả

Tỉ lệ cụm khổng lồ P_∞(τ) sụp đổ tại τ ≈ 1,85; ở đó x ≈ 0,347 = p_c(tam giác) = 2·sin(π/18) = 0,3473, tức z̄_c = z_max·p_c ≈ 2,08 — là kết quả, không phải giả định. Đây là ngưỡng liên thông (cụm xuyên suốt) và nó đánh dấu sự tan rã cụm gần điểm sôi (τ_b); không được nhầm nó với độ cứng vững Maxwell (z̄_c = 2d, nóng chảy) — xem sự phân tách rành mạch ở Ch. 5.4a.

<a id="kap5-4a"></a>

5.4.a Hai ngưỡng khác nhau được tách bạch

Cách gọi kép mơ hồ trước đây "thấm/độ cứng vững" nay được giải quyết. Đây là hai ngưỡng khác nhau về mặt vật lý, không được nhầm lẫn:

Lý do. Nóng chảy là sự mất độ cứng cơ học của một khối đặc, chứ không phải sự tan rã của tính liên kết. Do đó đại lượng liên quan là phép đếm ràng buộc Maxwell (trường trung bình / độ cứng vững tổng quát): một mạng D chiều gồm các lực trung tâm trở nên đẳng tĩnh đúng khi số phối trí trung bình đạt z̄ = 2D (2D → 4, 3D → 6). Tính liên thông liên kết (z̄_c = 2,08) nằm thấp hơn hẳn; một mạng vừa mới còn liên thông ở z̄ = 2,08 thì đã mềm oặt từ lâu (các mode mềm), tức là đã lỏng. Vì thế tính liên thông đánh dấu sự tan rã cụm muộn hơn, gần điểm sôi, còn độ cứng vững quy định sự đông đặc/hoá lỏng sớm hơn (mật độ nóng chảy). Sự phân tách này ấn định điểm ba ở Ch. 7.8 một cách không tham số và thực tế; nó không làm đổi ngưỡng liên thông 2,08 đã được khẳng định ở 5.4/5.6, mà chỉ đổi cách diễn giải nó (τ_b = gần sôi, không phải nóng chảy).

<a id="kap5-5"></a>

5.5 Bậc của chuyển pha nóng chảy

Tiêu chuẩn giải tích. Miền spinodal (hai cực tiểu của g) tồn tại khi τ < K/4; chuyển pha nằm xấp xỉ tại τ_t ≈ (ε + K/2)/ln G. Một bước nhảy thật (chuyển pha bậc nhất) đòi hỏi τ_t < K/4, tức trong giới hạn K → ∞:

ln G > 2 + 4ε/K.

Từ đó suy ra một kết quả đáng chú ý về độ suy biến trạng thái đứt G = z_max:

Xác nhận bằng số (|Δx| lớn nhất mỗi bước τ, detect_first_order):

Diễn giải. Trong mô hình 2D kinh điển (G = 6), sự nóng chảy về nguyên tắc là liên tục (một chuyển pha "mềm"/chuyển tiếp) — với không J nào có bước nhảy thật, vì ln 6 < 2. Chỉ số phối trí/độ suy biến cao hơn của cách xếp chặt 3D (G = 12) mới cho phép, từ mức hợp tác đủ lớn (K > 8,2), một chuyển pha bậc nhất với nhiệt ẩn hữu hạn. Điều đó khớp với thực tế: nóng chảy 2D thường liên tục (kiểu KTHNY), nóng chảy 3D là bậc nhất. J = 0 (không có trường phân tử) luôn là trường hợp quy chiếu liên tục. — Về mặt vật lý, J biểu thị tỉ số giữa năng lượng hợp tác và năng lượng liên kết (K/ε = J·z_max/ε): một liên kết lân cận nguyên vẹn làm bền liên kết của chính nó mạnh đến đâu.

Trình diễn bậc nhất (lấy cảm hứng 3D, cv_scan_09.py, J = 1,0, G = 12, z_max = 12): z̄ nhảy một lượng Δz̄ ≈ 5,30 tại τ ≈ 2,82; nhiệt ẩn dưới dạng bước nhảy năng lượng ΔU ≈ 1,75 ε/phân tử. Cùng K đó trên mạng G = 6 (2D) vẫn liên tục — tính chất bậc chỉ phụ thuộc vào G = z_max.

<a id="kap5-6"></a>

5.6 Nghiên cứu kích thước hữu hạn của ngưỡng thấm

Thấm liên kết trên mạng tam giác L × L (tuần hoàn), tỉ lệ cụm khổng lồ P_∞(p) được chia thang theo L (cv_scan_09.py; L = 50/100/200/400 với 8/5/3/2 lần thực hiện). Ngưỡng chính xác p_c = 2·sin(π/18) = 0,34730 → z̄_c = z_max·p_c = 2,0838.

P_∞(p) (trung bình ± độ tản):

Ước lượng ngưỡng p*(L) (P_∞ cắt 0,5):

Kết quả. p*(L) hội tụ đơn điệu từ dưới về p_c = 0,3473; độ tản tại p_c giảm từ ±0,103 (L = 50) xuống ±0,017 (L = 400) — chuyển pha sắc nét hơn theo L, nhưng không dịch ngưỡng (giá trị thấm phổ quát). Nhờ đó z̄_c ≈ 2,08 = z_max·p_c được khẳng định là giá trị mô hình bền vững với kích thước hữu hạn.

<a id="kap5-7"></a>

5.7 Phép quét c_v(τ) 3D đầy đủ của sự nóng chảy bậc nhất

Lần chạy trình diễn bậc nhất đơn lẻ ở Ch. 5.5 (J = 1,0, G = 12) nay được mở rộng thành một phép quét c_v(τ) đầy đủ theo độ hợp tác J (firstorder_3d_16.py, tính liên thông FCC cho n_blk, cân bằng năng lượng 3D U/N = τ·[3/2 + z̄/4 + (3/2)n_blk] − (ε/2)z̄ − (K/2)x²). Tiêu chuẩn bậc giải tích ở Ch. 5.5 (bậc nhất ⇔ K > K_tới hạn = 4ε/(ln 12 − 2) = 8,25) phân tách các trường hợp; đỉnh quy chiếu 2D liên tục nằm tại τ_2D = 1,20 (c_v ≈ 5,5).

*trường hợp liên tục: Δz̄/ΔU là chuyển tiếp dốc (không phải bước nhảy thật); τ_chuyển pha ở đó là đỉnh c_v.

Kết quả. (i) Trường hợp 3D liên tục J = 0,55 cho τ_chuyển pha = τ_m ≈ 1,72 — phép kiểm chéo với động học liên kết độc lập (bondkinetics_3d_15.py: τ_m(FCC) = 1,70) khớp. (ii) Từ K > K_tới hạn ≈ 8,25 (J > 0,68) mở ra bước nhảy thật Δz̄ với nhiệt ẩn hữu hạn ΔU (bước nhảy năng lượng qua chỗ gián đoạn) — sự nóng chảy trở thành bậc nhất, đúng như tiêu chuẩn ln 12 > 2 đòi hỏi. (iii) Phép kiểm chéo J = 1,0 tái tạo giá trị trình diễn ở Ch. 5.5 là Δz̄ ≈ 5,4. (iv) Độ cao bước nhảy Δz̄ và vị trí τ_chuyển pha tăng theo J: chuyển pha dịch lên trên — số phối trí cao hơn của cách xếp chặt 3D và độ hợp tác mạnh hơn đều nâng nhiệt độ nóng chảy. So với đỉnh 2D liên tục (τ_2D = 1,20), do đó luôn có một dịch chuyển đỉnh dương (+0,5 … +2,2), tăng đơn điệu theo K. Tính chất bậc nhất của sự nóng chảy 3D nhờ vậy không chỉ được chứng thực ở một điểm, mà trên toàn dải hợp tác.

<a id="kap6"></a>

6 — Tính hợp lý toán học của cách diễn giải 3D

Tóm tắt

Chương 6 làm cho việc tổng quát hoá bức tranh liên kết từ 2D sang 3D trở nên hợp lý thuần tuý về mặt toán học (không có hình ảnh HTML riêng):

• 6.1 Bậc tự do quay 2D → 3D: sự quay của khối tăng từ một lên ba trục, tịnh tiến từ hai lên ba hướng, z_max từ 6 lên 12; các giá trị cơ sở của c_v tăng (c_v,rắn 2,5 → 4,5, c_v,khí 1 → 3/2), hình dạng đường c_v(τ) vẫn giữ nguyên.
• 6.2 Tính liên thông như dấu hiệu phân biệt: z_max là một đại lượng lượng tử hoá và tách cách xếp chặt nhất (FCC/HCP, z_max = 12) khỏi BCC (z_max = 8) một cách nhảy vọt, không có trạng thái trung gian.
• 6.3 Chuyển pha bổ sung rắn(1)→rắn(2): từ tính rời rạc suy ra hiện tượng đa hình FCC↔BCC với bậc c_v không tham số 1,0 k_B; trong bức tranh lý tưởng nghiêm ngặt (θ_BCC = θ_FCC) nó không xuất hiện, chỉ một BCC mềm hơn (θ_BCC < θ_FCC dưới ngưỡng mở 0,61) mới mở ra τ_ss < τ_m.
• 6.4 / 6.5 Suy ngẫm: dị thường âm của nước lỏng như một sự đổi trật tự gần bị nhoè giữa hai kiểu mạng, và một phương án trường xoáy fractal thay cho sự quay mixen (thuần tuý định tính).
• 6.6 / 6.8 / 6.9 θ(mạng) vi mô ở ba bậc: mô hình lò xo Einstein (θ_BCC/θ_FCC ≈ 0,82…0,94), DOS phonon kiểu Debye (≈ 0,87…1,00) và mạng đầy đủ Lennard-Jones (≈ 1,0; BCC lý tưởng bất ổn định trượt, FCC là trạng thái cơ bản bền) — tất cả đều trên ngưỡng 0,61, vậy hiện tượng đa hình thực nằm ngoài bức tranh điều hoà lý tưởng.
• 6.7 Tính giả bền (2D): phép tính bậc độ lớn cho sự tái cấu trúc tam giác ↔ vuông (ΔU = +1 ε, ΔA ≈ +15,5 %) cùng rào tạo mầm ΔG\* = πλ²/Δμ_A như một bổ sung động học.

Chương này làm cho việc tổng quát hoá 3D trở nên hợp lý thuần tuý về mặt toán học. 3D không bao giờ được minh hoạ trong phiên bản HTML. Mục tiêu chỉ là chỉ ra bậc tự do quay và tính liên thông thay đổi thế nào khi đi từ 2D sang 3D, và những chuyển pha mới nào nhờ đó trở nên khả dĩ.

<a id="kap6-1"></a>

6.1 Thay đổi bậc tự do quay 2D → 3D

Điểm duy nhất mới về chất so với 2D là sự quay của các khối:

• 2D: một khối cứng có 1 bậc tự do quay (quay trong mặt phẳng) → U_rot/N = ½ n_blk τ.
• 3D: một khối cứng có 3 bậc tự do quay (ba trục Euler) → U_rot/N = (3/2) n_blk τ.

Đồng thời tịnh tiến tăng từ 2 lên 3 bậc tự do (U_trans/N = (3/2)τ) và số phối trí cực đại từ z_max = 6 (tam giác) lên z_max = 12 (xếp cầu chặt nhất). Cân bằng năng lượng trong 3D do đó là

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄,

và các giá trị cơ sở của nhiệt dung dịch chuyển một cách nhất quán:

Cấu trúc của đường cong (bình nguyên rắn, cực đại lỏng do −(ε/2)dz̄/dτ, sụt xuống giá trị khí) vẫn y nguyên; chỉ các mức là thay đổi theo chiều và z_max.

<a id="kap6-2"></a>

6.2 Ảnh hưởng của các mạng tinh thể khác nhau — tính liên thông như dấu hiệu phân biệt

Trong 3D phải phân biệt ít nhất hai mạng lý tưởng, tách nhau bởi tính liên thông (số phối trí z_max):

Tính liên thông z_max là dấu hiệu phân biệt quyết định: nó ấn định cả mức nhiệt dung của chất rắn (z_max/4) lẫn ngưỡng thấm z̄_c. Vì mỗi mạng chỉ cho phép một số phối trí cố định, tính liên thông là một đại lượng lượng tử hoá, hữu hạn — không có chuyển tiếp liên tục giữa z_max = 12 và z_max = 8, mà chỉ có bước nhảy rời rạc từ kiểu mạng này sang kiểu mạng khác.

<a id="kap6-3"></a>

6.3 Chuyển pha bổ sung rắn(1) → rắn(2)

Từ tính hữu hạn/rời rạc của tính liên thông suy ra một chuyển pha không có trong mô hình 2D: rắn(1) → rắn(2) (đa hình/thù hình). Hai pha rắn với z_max khác nhau (chẳng hạn FCC ↔ BCC) có các entanpi tự do μ(τ, p) khác nhau; tại tổ hợp nhiệt độ/áp suất mà μ_{rắn(1)} = μ_{rắn(2)}, chất đổi mạng của nó mà không nóng chảy. Trong bức tranh mô hình, hình học phối trí đổi một cách nhảy vọt (chứ không chỉ độ lấp đầy z̄ như khi nóng chảy) — tình huống duy nhất mà chính hình học quy chiếu vốn được giữ cố định lại nhảy đổi. Quan sát được sẽ là một bậc c_v nhỏ (biến thiên z_max/4) cộng nhiệt chuyển biến ẩn, yếu hơn hẳn đỉnh nóng chảy.

Tiêu chuẩn bậc nhất ln G > 2 + 4ε/K (Ch. 5.5) nay cho một lý do cụ thể vì sao chính cách diễn giải 3D lại phong phú hơn về chất: chỉ ở đó (G = z_max = 12, ln 12 = 2,485 > 2) sự nóng chảy mới có thể trở thành bậc nhất — và chỉ ở đó mới tồn tại sự đổi mạng rời rạc rắn(1) → rắn(2) (FCC ↔ BCC), vì tính liên thông là lượng tử hoá/hữu hạn. Cả hai hiện tượng là cùng một nguyên nhân trong hai lớp áo: số phối trí cao hơn, rời rạc của cách xếp chặt 3D sinh ra cả bước nhảy nóng chảy sắc nét lẫn hiện tượng đa hình. Mạng 2D đơn lẻ (G = 6) không có cả hai — nhất quán với Ch. 5.5.

Định lượng (bondkinetics_3d_15.py). Chuyển pha nay được tính gọn. Mỗi mạng nhận một entanpi tự do trong bức tranh điều hoà cổ điển (các phần p·v và tịnh tiến chung triệt tiêu khi so sánh):

μ(τ) = h₀ − c_v·τ·ln(τ/θ), h₀ = −(ε/2)·z_max, θ = ħω/k_B (độ cứng phonon theo đơn vị ε),

với các giá trị đặc trưng không tham số từ Ch. 6.1/6.2:

Từ đó bậc c_v không tham số: |Δc_v| = biến thiên z_max/4 = (12−8)/4 = 1,0 k_B (chất rắn khi chuyển FCC→BCC trở nên "nhẹ" hơn 1 k_B về nhiệt dung) — nhỏ hơn nhiều so với đỉnh nóng chảy (vài k_B, Ch. 5.3/7.11).

Nhiệt độ chuyển biến τ_ss (μ_FCC = μ_BCC) và nhiệt chuyển biến ẩn L = τ_ss·(S_BCC − S_FCC):

Kết quả (kết luận trung thực). Trong bức tranh lý tưởng nghiêm ngặt (cùng một độ cứng phonon cho mọi liên kết, θ_BCC = θ_FCC) không có chuyển pha rắn(1)→rắn(2): FCC thắng cả về entanpi (sâu hơn, −6 so với −4 ε) lẫn — nhờ số liên kết cao hơn (c_v = 4,5 > 3,5) — về entropi dao động; μ_FCC ở mọi nhiệt độ đều nằm dưới μ_BCC. Chỉ một độ cứng phonon phụ thuộc mạng θ_BCC < θ_FCC (BCC thoáng hơn dao động mềm hơn, tức có nhiều entropi hơn) mới hạ μ_BCC xuống dưới μ_FCC ở nhiệt độ cao — khi đó chuyển biến mở ra, và nó nằm dưới đỉnh nóng chảy τ_m, khiến trình tự rắn(1)→rắn(2)→lỏng→khí xuất hiện (τ_ss < τ_m < τ_c ≤ τ_b(3D), Ch. 7.11). Độ cứng BCC mềm là một mở rộng nhẹ ra ngoài bức tranh lý tưởng nghiêm ngặt (thực tế đã biết dưới tên đa hình γ/δ của sắt: BCC ở nhiệt độ cao, FCC ở giữa); ngược lại bậc c_v 1,0 k_B và nhiệt chuyển biến dương nhỏ là nội tại của mô hình và không tham số.

θ_BCC/θ_FCC thực sự nhỏ đến đâu? (theta_springs_16.py). Để giải phóng chuyển biến khỏi tham số θ tự do, độ cứng phonon nay được ước lượng vi mô (Ch. 6.6): một mô hình lò xo Einstein với các lò xo cặp xuyên tâm cho θ_BCC/θ_FCC ≈ 0,82 (chỉ lân cận gần nhất) đến ≈ 0,94 (kể cả lớp vỏ thứ hai của BCC). Trong khi đó ngưỡng mở cần thiết từ phép so sánh μ lại nằm ở θ_BCC/θ_FCC < 0,61 (dưới mức đó τ_ss tồn tại và nằm dưới τ_m). Vậy tỉ số độ cứng hình học là quá lớn: trong bức tranh lý tưởng nghiêm ngặt, FCC vẫn bền ở mọi nhiệt độ, không có chuyển pha rắn(1)→rắn(2) thuần tuý do hình học thúc đẩy. Hiện tượng đa hình thực (sắt γ↔δ) cần các đóng góp thêm (điện tử/từ/phi điều hoà) ngoài bức tranh điều hoà lý tưởng — chính là "mở rộng nhẹ" đã giả định ở trên. Điều đó xác nhận kết luận trung thực và biến tỉ số θ từ một giả định tự do thành một giới hạn mô hình được ước lượng và nêu tên.

<a id="kap6-4"></a>

6.4 Suy ngẫm: dị thường âm của nước lỏng

Các chất lỏng thông thường có v(T) giảm đơn điệu (Ch. 3.3). Nước lại có một cực đại âm ở ≈ 74 °C: v thoạt tiên tăng theo T trước khi giảm. Trong mô hình, điều đó có thể được diễn giải như sự cạnh tranh của hai kiểu mạng cục bộ (Ch. 6.2): một trật tự gần giống nước đá thoáng, phối trí thấp (z_max nhỏ, nhiều "lỗ") và một trật tự gần đặc, phối trí cao hơn. Ở T thấp, cấu trúc thoáng chiếm ưu thế; đun nóng trước hết phá vỡ trật tự gần thoáng để nhường chỗ cho trật tự đặc hơn — tính liên thông hiệu dụng và cùng với nó K_eff tăng (v tăng). Chỉ trên cực đại, sự đứt liên kết VẾT NỨT thông thường mới chiếm ưu thế, K_eff giảm, v hạ xuống. Dị thường vì vậy là biểu hiện của một sự đổi trật tự gần kiểu rắn(1)↔rắn(2) bị nhoè và không trọn vẹn bên trong pha lỏng — chỉ khả dĩ trong mô hình 3D với hai kiểu mạng, chứ không trong mạng 2D đơn lẻ.

<a id="kap6-5"></a>

6.5 Phương án khác: trường xoáy fractal thay cho sự quay mixen

Chuyển động Brown cục bộ (Ch. 3.4) đã được giải thích qua sự quay và trao đổi mixen. Một cách khác, nó có thể được hình dung như hiệu ứng của các trường xoáy fractal: các biên khối liên tục được VẾT NỨT tái cấu trúc tạo ra một trường các dòng động lượng cục bộ, mà — tương tự các dòng chảy rối — mang những xoáy tự đồng dạng trên mọi thang chiều dài. Khi đó hạt thử không bị dẫn bởi các va chạm mixen rời rạc, mà bởi một liên tục xoáy bất biến thang. Cả hai bức tranh (mixen rời rạc ↔ xoáy fractal) là những trường hợp giới hạn bổ sung nhau của cùng một động học VẾT NỨT; bức tranh nào đúng phụ thuộc vào tỉ số giữa kích thước hạt thử và kích thước mixen trung bình. Phần bàn sơ này thuần tuý định tính và không được triển khai toán học ở đây.

<a id="kap6-6"></a>

6.6 Độ cứng phonon vi mô θ(mạng): một mô hình lò xo Einstein

Phép so sánh rắn(1)→rắn(2) ở Ch. 6.3 phụ thuộc vào độ cứng phonon tương đối θ_BCC/θ_FCC. Thay vì chọn nó tự do, ở đây nó được ước lượng từ số phối trí và hình học (theta_springs_16.py).

Mô hình. Một nguyên tử trung tâm nằm trong hố điều hoà của z lân cận của nó; mỗi liên kết là một lò xo cặp xuyên tâm κ. Một độ dời u chiếu đẳng hướng với ⟨cos²⟩ = 1/d = 1/3 (d = 3) lên các hướng liên kết, sao cho hằng số lò xo hiệu dụng là k_eff = (z/3)·κ. Tần số Einstein là ω = √(k_eff/m), tức θ = ħω/k_B ∝ √k_eff. Với cùng khối lượng và cùng lò xo cặp thì suy ra

θ_BCC/θ_FCC = √(z_BCC/z_FCC) = √(8/12) = 0,816.

Tinh chỉnh. (B) Lớp vỏ thứ hai của BCC (6 nguyên tử tại r₂/r₁ = 2/√3 = 1,155) với lò xo tỉ lệ theo khoảng cách κ(r) ∝ (r₁/r₂)ⁿ (n = 6…8, kiểu Lennard-Jones) nâng tỉ số lên 0,91…0,94 (BCC trở nên cứng hơn một cách hiệu dụng nhờ lớp vỏ thứ hai ở gần). (C) Liên kết ngắn hơn của BCC ở cùng mật độ (r_nn(BCC)/r_nn(FCC) = 0,972, −2,8 %) làm lò xo cứng thêm chút ít và cho 0,89…0,91. Tổng cộng là một cửa sổ hình học θ_BCC/θ_FCC ≈ 0,82…0,94.

Ngưỡng mở. Phép so sánh μ (Ch. 6.3) chỉ mở chuyển biến (τ_ss tồn tại và nằm dưới τ_m = 1,70) khi θ_BCC/θ_FCC < 0,61. Cửa sổ hình học nằm cao hơn hẳn.

Kết luận. Mô hình lò xo thuần tuý không làm BCC mềm đủ: trong bức tranh lý tưởng nghiêm ngặt, FCC vẫn bền, một chuyển pha rắn(1)→rắn(2) do hình học thúc đẩy không xuất hiện. Nhờ đó θ_BCC/θ_FCC đã đi từ một giả định tự do thành một đại lượng được ước lượng, mà giá trị tự nhiên của nó (≈ 0,8–0,9) giữ chuyển biến ở trạng thái đóng — nhất quán với kết luận trung thực ở Ch. 6.3. Hiện tượng đa hình thực vì vậy được chứng minh là nằm ngoài bức tranh điều hoà lý tưởng (phần điện tử/từ, tính phi điều hoà), chứ không ở riêng hình học mạng.

<a id="kap6-7"></a>

6.7 Ước lượng toán học tính giả bền (2D tam giác ↔ vuông)

Cho phần thảo luận định tính ở Ch. 3.2, đây là phép tính bậc độ lớn tương ứng trong bức tranh 2D (đĩa cứng, đường kính D, khoảng cách tiếp xúc a = D).

(1) Năng lượng. Trên mỗi hạt U/N = −(ε/2)·z (mỗi liên kết tính một nửa): tam giác z = 6 → −3 ε, vuông z = 4 → −2 ε. Vậy sự tái cấu trúc tam giác → vuông tốn ΔU = +1 ε trên mỗi hạt (thu nhiệt; dạng lỏng lẻo giàu năng lượng hơn).

(2) Thể tích (diện tích) và công. Diện tích trên mỗi hạt: tam giác A_tri = (√3/2)·D² ≈ 0,866 D², vuông A_quad = D². Sự nới rộng là ΔA = A_quad − A_tri ≈ 0,134 D² (+15,5 %) — chính là sự tăng thể tích đã được thể hiện trong truyện tranh. Chống lại một áp suất p (rút gọn) phải bỏ ra công W = p·ΔA; entanpi tự do tăng thêm Δμ = ΔU + p·ΔA − τ·ΔS. Số hạng entropi ΔS > 0 (dạng lỏng lẻo có các mode mềm hơn, bậc c_v lớn hơn Δc_v = (6−4)/4 = 0,5 k_B) chỉ hạ μ_quad xuống dưới μ_tri ở τ cao — vì thế tam giác là trạng thái cơ bản ở nhiệt độ thấp và vuông là dạng nhiệt độ cao/giả bền.

(3) Rào tạo mầm (2D). Tính giả bền đòi hỏi một rào. Trong 2D, một mầm của pha mới (bán kính r) lớn lên với phần lợi về diện tích −Δμ_A·π r² (Δμ_A = động lực trên đơn vị diện tích) chống lại đóng góp biên +λ·2π r (λ = năng lượng cạnh trên đơn vị chiều dài): ΔG(r) = −Δμ_A·π r² + 2π λ r. Cực đại tại r\* = λ / Δμ_A với rào ΔG\* = π λ² / Δμ_A. Mầm tới hạn r\ co lại khi động lực Δμ_A tăng; chừng nào thang nhiệt τ = k_BT chưa đạt tới rào (**ΔG\ ≳ vài k_BT), pha lỏng lẻo (hoặc, khi sôi, pha quá bão hoà) vẫn bị đóng băng về động học. Với các cụm nhỏ (2–3 hạt, như trong dải truyện về quá bão hoà), tỉ số biên trên diện tích là bất lợi nhất — chúng nằm bên trái r\* và trung bình lại tan ra. Như vậy năng lượng hoạt hoá được thể hiện trong truyện tranh chính là đại lượng ΔG\* ấy, đại lượng nâng đỡ cả dạng tinh thể giả bền lẫn pha khí quá bão hoà; nó là một bổ sung động học** cho mô hình cân bằng lý tưởng (không có rào).

<a id="kap6-8"></a>

6.8 Tinh chỉnh kiểu Debye cho θ(mạng)

Mô hình lò xo Einstein (Ch. 6.6) thay toàn bộ phân bố phonon bằng một tần số. Cách tiếp cận hiện tại tinh chỉnh θ(mạng) bằng một mật độ trạng thái phonon gọn nhẹ cho mỗi mạng (theta_debye_17.py): cũng những lò xo cặp xuyên tâm κ(r) = r^(−n) ấy trên lớp vỏ lân cận gần nhất và thứ hai tạo nên ma trận động lực D(k) = Σ_j κ_j·(R̂_j⊗R̂_j)·(1 − cos(k·R_j)), được chéo hoá trên một lưới k của vùng Brillouin nguyên thuỷ (Monkhorst-Pack, 12³ điểm) (bộ giải trị riêng 3×3 giải tích, chỉ dùng thư viện chuẩn). Từ ba nhánh phonon suy ra hai thước đo độ cứng kiểu Debye dưới dạng mômen tần số:

• Mômen bậc 2 θ(2) ∝ √⟨ω²⟩ — trường hợp giới hạn Einstein (trung bình BZ của 1 − cos cho 1, tức ⟨ω²⟩ = z·κ/3); kiểm tra chỉ với lân cận gần nhất, thuần hình học: √⟨ω²⟩_BCC / √⟨ω²⟩_FCC = 0,816 = √(8/12) — đúng giá trị Einstein ở Ch. 6.6.
• Mômen log θ(0) ∝ exp⟨ln ω⟩ — đại lượng thực sự kiểu Debye, có trọng số cho các nhánh âm mềm mà giá trị Einstein bỏ qua.

Mọi hằng số mang thang (ħ/k_B, thang κ, tiền nhân Debye) đều triệt tiêu trong tỉ số θ_BCC/θ_FCC; chỉ còn lại hình dạng DOS. Với lớp vỏ lân cận thứ hai (chính nó mới làm BCC bền về mặt cơ học) và định luật lò xo n = 0…8 ta được

θ_BCC/θ_FCC ≈ 0,87…1,00 (dải sai số theo định nghĩa mômen và định luật lò xo).

Con số này thiên về cao hơn giá trị Einstein thuần tuý (0,82…0,94), vì liên kết BCC ngắn hơn ở cùng mật độ làm cứng thêm các lò xo. Ngưỡng mở 0,61 cần thiết từ phép so sánh μ (Ch. 6.3) vẫn nằm thấp hơn hẳn: DOS phonon cũng không làm BCC mềm đủ. Kết luận ở Ch. 6.6 vì thế được xác nhận và chính xác hoá bởi phép gần đúng tinh hơn — trong bức tranh lý tưởng nghiêm ngặt, FCC vẫn bền ở mọi nhiệt độ; một chuyển pha rắn(1)→rắn(2) thuần tuý do hình học thúc đẩy không xuất hiện, hiện tượng đa hình thực (sắt, Ch. 7.13) cần các đóng góp ngoài bức tranh điều hoà lý tưởng.

<a id="kap6-9"></a>

6.9 θ(mạng) từ một mạng đầy đủ Lennard-Jones

Ch. 6.8 dùng một lò xo xuyên tâm bị cắt κ(r) = r^(−n) trên nn + 2nn. Bản 17 kiểm tra cùng câu hỏi đó với một hằng số lực thật từ một thế cặp trên toàn mạng (theta_ljlattice_18.py): mỗi liên kết trở thành một lò xo lấy từ thế Lennard-Jones φ(r) = 4ε[(σ/r)¹² − (σ/r)⁶] với phần dọc κ_L = φ″(r) và phần ngang κ_T = φ′(r)/r. Ma trận động lực

D_ab(k) = Σ_j [ (κ_L − κ_T)·û_a û_b + κ_T·δ_ab ]·(1 − cos(k·R_j))

lấy tổng trên mọi lớp vỏ lân cận đến một ngưỡng cắt (2,0…3,0 nn); mỗi mạng trước đó được hồi phục về cân bằng LJ của chính nó (hằng số mạng tự do, không áp đặt cùng mật độ). Đó là phép kiểm tra trung thực xem phép gần đúng gọn nn+2nn (Ch. 6.8) có nắm đúng trật tự độ cứng hay không.

Kết quả (tất định, chỉ dùng thư viện chuẩn):

Hai kết luận, cả hai đều làm chặt thêm Ch. 6.6/6.8:

1. FCC là trạng thái cơ bản LJ bền (E ≈ −8,4 ε/nguyên tử ở ngưỡng cắt 3,0 nn, so với −8,0 ε cho BCC) — ưu thế đã biết của cách xếp chặt trong thế cặp.
2. BCC LJ lý tưởng bất ổn định trượt: khoảng 3…3,5 % số mode phonon là ảo (ω² < 0, bất ổn định Born). Bức tranh lý tưởng thuần thế cặp vì vậy hoàn toàn không mang một BCC bền về mặt cơ học; trên các mode thực còn lại, θ_BCC/θ_FCC ≈ 1,0 — ngưỡng mở 0,61 cần thiết từ phép so sánh μ lại càng không bị vượt xuống dưới.

Nhờ đó, phát biểu ở Ch. 6.6/6.8 là bền vững trước sự tinh chỉnh từ lò xo bị cắt sang toàn mạng và từ mô hình thuần dọc sang ma trận động lực đầy đủ: trong bức tranh lý tưởng cổ điển không xuất hiện chuyển pha rắn(1)→rắn(2) thuần tuý hình học — hiện tượng đa hình thực (sắt Ch. 7.13, titan Ch. 7.15) cần các đóng góp ngoài thế cặp (điện tử/từ, phi điều hoà). Một phép tính phonon ab initio thực sự được cố ý để ra ngoài (nó sẽ rời khỏi tính chất gọn nhẹ của bức tranh lý tưởng).

<a id="kap7"></a>

7 — Hiệu chuẩn và phương trình trạng thái

Tóm tắt. Chương 7 hiệu chuẩn mô hình qua một tham số năng lượng ε duy nhất lên các chất thực và mở rộng nhánh lỏng/khí thành phương trình trạng thái:

• 7.1 / 7.2 Phương pháp và bảng hiệu chuẩn: gắn nhiệt độ tuyến tính T = (ε/k_B)·τ với một điểm neo tại điểm nóng chảy; từ đó tỉ số phổ quát đã tinh chỉnh R = τ_b/τ_m ≈ 1,54.
• 7.3 / 7.4 Phương trình trạng thái van der Waals lõi cứng: phương trình rút gọn p_r = 8τ_r/(3v_r − 1) − 3/v_r² với điểm tới hạn (1, 1, 1) và sự cùng tồn tại Maxwell; a ∝ ε, b ∝ D². Các điểm sôi nằm một cách hệ thống tại T_b/T_c ≈ 0,48 (khoảng cách Guldberg); một phép tính so sánh Redlich–Kwong nâng chúng lên ≈ 0,60 — như một ngân sách sai số được nêu tên, không đổi phương trình trạng thái.
• 7.5 Mặt (τ, v) chung: vòm cùng tồn tại (sôi), đường nóng chảy và tường cứng ghép lại với một điểm ba (τ_r,triple ≈ 0,47 với v_solid = 0,40).
• 7.6 Nhiệt ẩn: L ≈ 1,41 ε cho ΔS_fus ≈ 1,175 R (gần quy tắc Richards); các kim loại được trúng xuất sắc (đồng ×1,00), tinh thể ion/liên kết hiđro/phân tử nằm ngoài bức tranh lý tưởng.
• 7.7 Thống nhất ε ↔ ε_lg: k_B T_c = z_max·ε/4 nối hai thang (ε_lg = (27/32)·z_max·ε); tiên đoán 3D τ_c/τ_m = 2,50 trúng các chất phân tử đơn giản, 2D vẫn là trường hợp đặc biệt.
• 7.8 / 7.9 Điểm ba không tham số: v_solid suy ra từ ngưỡng cứng vững z̄_c = 2d, không từ tính liên thông; trong 3D mở ra một cửa sổ lỏng sạch (τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00), trong 2D nó vẫn đóng (vai trò đặc biệt của 2D được chứng thực ba lần).
• 7.10 / 7.11 Các phép quét 3D thực: thấm liên kết trên FCC/BCC cho z̄_c(3D) ≈ 1,5–1,6; phép quét τ của động học liên kết cho thấy τ_m ≈ τ_rig (nóng chảy = độ cứng vững, thấp) và τ_b(3D) ≈ τ_c (sôi = tan rã tính liên thông tại điểm tới hạn).
• 7.12 / 7.13 / 7.15 đường rắn(1)→rắn(2): từ điểm τ_ss, qua áp suất trở thành một đường (dp/dτ = ΔS/Δv > 0, Δv = +8,9 %), như đường biên thứ tư trong mặt (τ, v) và như một đường (τ, p) hiệu chuẩn theo sắt (γ→δ), được xác nhận trên titan (α→β) và với argon làm đối chứng âm.
• 7.14 Nhánh Redlich–Kwong: như một biến thể bức tranh lý tưởng độc lập với hiệu chuẩn riêng, bên cạnh mặc định van der Waals (T_b/T_c ≈ 0,60, sai số −1 %).

<a id="kap7-1"></a>

7.1 Phương pháp hiệu chuẩn (một tham số năng lượng ε) — đã điều chỉnh

Tuyến tính: T = (ε/k_B)·τ, một điểm neo cho mỗi chất. Neo tại điểm nóng chảy: ε/k_B = T_m,tn/τ_m, tiên đoán T_b,pred = T_m,tn·(τ_b/τ_m). Với đỉnh nóng chảy đã tinh chỉnh (số hạng hợp tác), cực đại c_v nay nằm tại τ_m = 1,20 (thay vì 1,15), còn thấm tính liên thông (gần sôi) không đổi tại τ_b = 1,85. Từ đó tỉ số phổ quát đã tinh chỉnh:

R = τ_b / τ_m = 1,85 / 1,20 ≈ 1,54 (bản 08: 1,61).

Bảng hiệu chuẩn của bản 08 (Ch. 7.2) tỉ lệ tuyến tính theo R: mọi T_b,pred giảm đi hệ số 1,54/1,61 ≈ 0,957 (≈ −4,3 %). Điều đó cải thiện đôi chút sự ước lượng quá cao đối với khí hiếm (argon +54 % → +48 %) và làm xấu đi chút ít độ trúng của NaCl/sắt (NaCl −0,6 % → −4,8 %). Kết luận nền tảng vẫn giữ: R phù hợp với các chất xếp chặt đẳng hướng, liên kết không định hướng; khí hiếm/mạch nằm ngoài bức tranh lý tưởng.

<a id="kap7-2"></a>

7.2 Bảng hiệu chuẩn

Được tiếp nhận không đổi từ modell_08.md Ch. 7.2 (với lưu ý về việc chia thang tuyến tính theo R ở 7.1). Tái lập được qua cv_scan_08.py (bảng) hoặc cv_scan_09.py (τ_m, τ_b, R).

<a id="kap7-3"></a>

7.3 Trục mật độ/thể tích: phương trình trạng thái van der Waals lõi cứng

Áp suất được nâng từ "tham số" lên biến số thực thụ. LÕI CỨNG cho thể tích riêng b, sức hút liên kết (ε·z̄) cho áp suất kết dính a. Trong đơn vị rút gọn (điểm tới hạn tại (v_r, τ_r, p_r) = (1, 1, 1)):

p_r = 8 τ_r / (3 v_r − 1) − 3 / v_r², v_r > 1/3 (tường cứng).

Sự cùng tồn tại theo dựng hình Maxwell (cv_scan_09.py, cắt diện tích bằng nhau):

Kết quả. Dưới τ_c mở ra bước nhảy mật độ l/k (v_gas/v_liq tăng từ 2,5 lên 11,5; càng lạnh càng lớn). Đó là đường sôi — nay đã tách sạch khỏi đỉnh c_v nóng chảy (Ch. 5.3), vốn sống trên mạng có mật độ cố định. Nhờ đó điểm 3 còn để ngỏ ở bản 08 được khép lại: nóng chảy (trật tự độ cứng vững/liên kết ở mật độ cố định) và sôi (bước nhảy mật độ/thể tích) là những trục khác nhau của cùng một mô hình.

<a id="kap7-4"></a>

7.4 Hiệu chuẩn l/g định lượng: a → ε, b → D²

Ánh xạ các đại lượng rút gọn lên chất thực. Các hằng số van der Waals gắn với điểm tới hạn:

b = k_B·T_c / (8 p_c) (thể tích riêng, LÕI CỨNG), a = 27 k_B²·T_c² / (64 p_c) (kết dính, LIÊN KẾT).

Do đó a ∝ ε (sự kết dính bắt nguồn từ năng lượng liên kết) và b ∝ D² (thể tích riêng 2D từ đường kính lõi cứng). Thang năng lượng tự nhiên của nhánh l/g là

ε_lg / k_B = (a/b) / k_B = (27/8)·T_c = 3,375·T_c.

Chất neo argon (T_c = 150,87 K, p_c = 48,63 bar): b = 5,35·10⁻²⁹ m³/phân tử, a = 3,76·10⁻⁴⁹ J·m³/phân tử², ε_lg/k_B = 509 K và đường kính phân tử D ≈ b^{1/3} = 3,8·10⁻¹⁰ m — so với đường kính van der Waals của argon (~3,4·10⁻¹⁰ m) chính xác tới 11 %.

Tiên đoán nhiệt độ sôi. Ở áp suất ngoài cố định, điểm sôi ở nơi áp suất bão hoà rút gọn thoả p_sat(τ_r)·p_c = p_atm; từ đó T_b = T_c·τ_r,b (đường bão hoà Maxwell từ lg_calib_10.py):

*CO₂ thăng hoa ở 1 atm (T = 194,7 K = điểm thăng hoa).

Kết quả. Với các chất phân tử (p_c ≈ 30–70 bar), mô hình cho nhất loạt T_b/T_c ≈ 0,48–0,51, một cách hệ thống thấp hơn ~17–24 % so với giá trị đo. Đó không phải lỗi mô hình theo nghĩa hẹp, mà là khoảng cách đã biết của phương trình van der Waals so với quy tắc Guldberg thực nghiệm (T_b/T_c ≈ 0,6): đường đẳng nhiệt van der Waals đánh giá thấp áp suất bão hoà ở nhiệt độ thấp. Nước (và các kim loại, vốn T_c/p_c cũng chỉ được ước lượng thô) ở p_r,atm < 0,005 nằm ngoài vùng van der Waals có nghĩa (τ_r,b < 0,40) — ở đó sự kết dính không còn kiểu van der Waals nữa (liên kết hiđro, liên kết kim loại). Bậc độ lớn và thứ tự các điểm sôi được trúng đúng.

Xếp loại theo Guldberg như một bảng sai số. Khoảng cách này ở đây chỉ được nêu tên, chứ không được sửa bằng cách đổi phương trình trạng thái (bức tranh lý tưởng vẫn kiểu van der Waals). Quy tắc Guldberg thực nghiệm (Guldberg 1890) nói rằng điểm sôi nằm khoảng T_b/T_c ≈ 0,6 (phụ thuộc chất nhẹ ~0,58–0,64). Mô hình van der Waals rút gọn lại cho một giá trị không đổi T_b/T_c = τ_r,b ≈ 0,48–0,51:

Diễn giải. Giá trị trung bình của mô hình 0,49 trúng hằng số Guldberg 0,60 với sai lệch khoảng −18 %; độ lệch có hệ thống và một chiều (luôn quá thấp), tức là một tính chất của đường đẳng nhiệt van der Waals (áp suất bão hoà quá nhỏ ở τ thấp), chứ không phải của bức tranh phân tử. Một phương trình trạng thái kết dính tốt hơn (ví dụ Redlich–Kwong) sẽ nâng τ_r,b — nhưng ở đây vẫn có ý thức nằm ngoài bức tranh lý tưởng.

Ngân sách sai số với một phép tính so sánh Redlich–Kwong (guldberg_rk_15.py). Để định lượng (chứ không chữa) khoảng cách của van der Waals so với quy tắc Guldberg, cùng phép thử điểm sôi được lặp lại với sự kết dính kiểu Redlich–Kwong, phụ thuộc nhiệt độ a/√T. Trong đơn vị rút gọn (điểm tới hạn (v_r, τ_r, p_r) = (1, 1, 1), Z_c = 1/3):

p_r = 3 τ_r/(v_r − β) − A/(√τ_r · v_r(v_r + β)), β = 3·Ω_b = 0,25992, A = 9·Ω_a = 3,84732,

với dựng hình diện tích bằng nhau của Maxwell và cùng áp suất ngoài rút gọn p_r,atm như trên. Kết quả (một con số cho mỗi chất neo, T_b/T_c = τ_r,b):

Kết quả. Sự kết dính RK phụ thuộc nhiệt độ nâng T_b/T_c từ vdW ≈ 0,48 lên ≈ 0,60 — gần như đúng hằng số Guldberg và các giá trị thực (≈ 0,58–0,63); sai số trung bình so với các giá trị đo giảm từ −21 % (vdW) xuống −1 % (RK). Nhờ đó nguyên nhân của khoảng cách Guldberg được định lượng rõ ràng: chính là sự thiếu phụ thuộc nhiệt độ của sự kết dính van der Waals (áp suất bão hoà quá nhỏ ở τ thấp), chứ không phải bức tranh phân tử. Phép tính so sánh này có ý thức nằm ngoài bức tranh lý tưởng; mô hình vẫn kiểu van der Waals (không đổi phương trình trạng thái) — ngân sách sai số chỉ được nêu tên, chứ không được lắp vào.

<a id="kap7-5"></a>

7.5 Mặt (τ, v) chung với điểm ba

Trục nóng chảy và trục sôi cho tới nay sống tách biệt: nóng chảy (trật tự độ cứng vững/liên kết) ở mật độ cố định, sôi (van der Waals) qua bước nhảy mật độ. Trong một mặt phẳng (τ_r, v_r) chung chúng ghép lại với nhau:

• Vòm cùng tồn tại (sôi): v_liq(τ_r) (nhánh trái) và v_gas(τ_r) (nhánh phải) từ dựng hình Maxwell, gặp nhau tại điểm tới hạn (τ_r, v_r) = (1, 1).
• Chất rắn (nóng chảy): dải hẹp ngay trên tường cứng v_r > 1/3, ở mật độ cao và gần như cố định (thể tích rút gọn v_solid ≈ 0,40).
• Đường nóng chảy: gần như thẳng đứng tại v_solid (mật độ cố định), là nhiệt độ của độ cứng vững Maxwell (z̄_c = 2d; Ch. 7.8 suy ra v_solid từ đó một cách không tham số — cách gán cũ của bản 11 "z̄_c ≈ 2,08" nhờ đó được sửa).

Điểm ba là điểm nơi biên cùng tồn tại lỏng v_liq(τ_r) đạt tới mật độ xếp chặt rắn v_solid — dưới đó không còn nhánh lỏng nữa (thăng hoa trực tiếp):

Kết quả. Với v_solid = 0,40, điểm ba nằm tại τ_r,triple ≈ 0,47, p_sat ≈ 0,018, với thể tích hơi v_gas ≈ 68 (hơi thăng hoa rất loãng). Tỉ số τ_triple/τ_c ≈ 0,47 nằm trong cửa sổ thực (nước 0,42; argon 0,56; CO₂ 0,71). Chất rắn càng đặc (v_solid càng nhỏ) thì điểm ba rút gọn càng thấp và hơi cùng tồn tại càng loãng — đúng về mặt định tính. Nhờ đó cả hai chuyển pha và điểm ba được neo trong một mặt (τ, v); điểm tới hạn ngồi tại (1, 1, 1). Bản 15 bổ sung ở Ch. 7.12 một đường biên thứ tư — chuyển biến rắn(1)/rắn(2) (FCC↔BCC) như hai dải mật độ rắn gần thẳng đứng (v ≈ 0,667 và 0,726), chuyển sang nhau tại τ_ss.

Giai đoạn phát triển trước đó: v_solid = 0,40 ở đây còn được đặt ra. Ch. 7.8 suy ra v_solid không tham số từ ngưỡng thấm z̄_c và biến điểm ba thành một tiên đoán thực sự.

<a id="kap7-6"></a>

7.6 Nhiệt ẩn so với entanpi nóng chảy

Nhiệt ẩn của mô hình L ≈ 1,41 ε (Ch. 5.3) được chuyển sang đơn vị thực qua hiệu chuẩn điểm nóng chảy ε/k_B = T_m/τ_m (τ_m = 1,20):

ΔH_m,pred = L·R·T_m/τ_m = 1,175·R·T_m, tương ứng ΔS_fus,pred = L·R/τ_m = 1,175·R.

Đó — điều đáng chú ý — gần như đúng quy tắc Richards (ΔS_fus ≈ R cho các chất rắn cấu tạo đơn giản). So sánh với các entanpi nóng chảy đo được (lg_calib_10.py):

Kết quả. Mô hình trúng bậc độ lớn (kJ/mol) một cách nhất loạt, và trúng quy tắc Richards cho các kim loại đơn giản một cách xuất sắc (đồng ×1,00, sắt ×1,28); các khí đã hoá lỏng nằm hơi thấp hơn (metan ×0,94, nitơ ×0,87, argon ×0,69). Tinh thể ion, liên kết hiđro và phân tử (NaCl ×0,37, nước ×0,44, benzen ×0,28, CO₂ ×0,23) có entropi nóng chảy thực cao hơn rõ rệt (ΔS_exp/R = 2,6–5,0), vì khi nóng chảy còn giải phóng thêm những bậc trật tự (sự ghép cặp ion, sự định hướng của phân tử nhiều nguyên tử) mà mô hình một tham số trật tự không nắm bắt được. Giới hạn nhờ đó được nêu rõ: bức tranh lý tưởng áp dụng cho các chất liên kết đẳng hướng, không định hướng (không có khớp số liệu).

<a id="kap7-7"></a>

7.7 Thống nhất các thang năng lượng: ε ↔ ε_lg

Cho tới nay nhánh nóng chảy mang năng lượng liên kết ε (τ_m = 1,20), còn nhánh l/g mang một thang kết dính riêng ε_lg/k_B = 3,375·T_c (Ch. 7.4). Cả hai đều có thể suy ra từ một tham số mô hình. Sự kết dính van der Waals a trong mô hình sinh ra từ tổng các năng lượng liên kết trên mỗi phân tử (LIÊN KẾT): mỗi phân tử mang z_max/2 liên kết có năng lượng ε. Một khí mạng Bragg–Williams trường trung bình với lực hút lân cận ε cho hệ thức cổ điển đối với nhiệt độ tới hạn:

k_B·T_c = z_max·ε / 4, tức là τ_c ≡ k_B T_c/ε = z_max/4.

Thế điều này vào ε_lg = a/b = (27/8)·k_B·T_c (Ch. 7.4), ta được sự nối kết của hai thang:

ε_lg = (27/8)·(z_max/4)·ε = (27/32)·z_max·ε, hằng số c = 27/32 = 0,844.

Nhờ đó đỉnh nóng chảy, sự thấm và điểm tới hạn nằm trên một thang ε (unify_scale_12.py):

Tiên đoán τ_c/τ_m = z_max/(4·τ_m): 2D → 1,25, 3D → 2,50. So với các chất thực:

Kết quả. Tiên đoán 3D τ_c/τ_m = 2,50 trúng khá tốt chỉ số thực của các chất phân tử đơn giản (T_c/T_m ≈ 2,0–2,1); các kim loại nằm cao hơn (T_c ước lượng thô, kiểu liên kết khác). Tiên đoán 2D (1,25) quá thấp — và τ_c(2D) = 1,50 thậm chí nằm dưới τ_b = 1,85 (điểm tới hạn sẽ nằm dưới ngưỡng thấm). Đó không phải lỗi tính toán, mà là vai trò đặc biệt của 2D đã biết: ở đó thang kết dính và thang liên kết không sắp xếp gọn ghẽ, thiếu một cửa sổ lỏng bền vững (xem sự nóng chảy 2D liên tục, Ch. 5.5). Chỉ trong 3D (τ_m = 1,20 < τ_c = 3,00) các thang mới tách nhau một cách có ý nghĩa vật lý.

<a id="kap7-8"></a>

7.8 Điểm ba không tham số từ z̄_c

Ch. 7.5 còn đặt ra v_solid = 0,40. Giai đoạn phát triển trước đó suy ra mật độ xếp chặt rắn từ ngưỡng thấm. Số phối trí của một mạng ở thể tích rút gọn v_r tuân theo hệ thức mật độ–phối trí (tường cứng v_r = 1/3 ↔ xếp chặt nhất, z̄ = z_max):

z̄(v_r) = z_max·(1/3)/v_r = z_max/(3 v_r), ngược lại v_solid = z_max/(3·z̄_c).

Chất rắn mất độ cứng ngay khi z̄ xuống dưới ngưỡng z̄_c; mật độ tương ứng là v_solid. Điểm ba khi đó là điểm nơi biên cùng tồn tại lỏng v_liq(τ_r) đạt tới mật độ ấy (unify_scale_12.py, 2D, z_max = 6):

Bản cải tiến. Mật độ nóng chảy v_solid được ấn định bởi độ cứng vững Maxwell tổng quát z̄_c = 2d (2D → 4, 3D → 6), không bởi tính liên thông liên kết (lý do ở Ch. 5.4a: nóng chảy là sự mất độ cứng, không phải mất tính liên kết). Do đó điểm ba đáng tin cậy là (rigidity_3d_13.py, 2D):

Kết quả 2D. Ngay cả với ngưỡng cứng vững (đúng về mặt vật lý), mạng 2D đơn lẻ không cho một cửa sổ lỏng sạch: v_solid = 0,50 tuy cho chỉ số thực τ_triple/τ_c ≈ 0,77, nhưng trong đơn vị ε tuyệt đối, τ_triple = 0,77·1,50 = 1,16 và do đó nằm dưới τ_m = 1,20. Ngưỡng liên thông (z̄_c = 2,08) thậm chí đẩy v_solid ≈ 0,96 ≈ v_c ra khỏi vòm cùng tồn tại — điểm ba và điểm tới hạn trùng nhau. Cả hai phép tính đều xác nhận vai trò đặc biệt của 2D (Ch. 5.5, 7.7): mạng 2D đơn lẻ hầu như không mang một pha lỏng độc lập. Đó là một kết luận trung thực, không phải điểm yếu tính toán hay mô hình theo nghĩa hẹp.

Phép thử quyết định vì thế là phép tính kiểm tra 3D (Ch. 7.9): ở đó chính quy tắc cứng vững ấy (z̄_c = 6 → v_solid = 0,667) mở ra một cửa sổ sạch. Trong mọi trường hợp, điểm ba nay được ràng buộc không tham số vào z̄_c — tham số tự do v_solid đã bị loại bỏ.

<a id="kap7-9"></a>

7.9 Tính toán kiểm tra 3D: cửa sổ lỏng sạch trên một thang ε

Kiểm tra xem sự sụp đổ 2D có được hoá giải trong 3D hay không, khi mọi mốc nằm trên một thang ε (τ_c = z_max/4 từ Ch. 7.7; v_solid từ độ cứng vững tổng quát z̄_c = 2d; điểm ba ở nơi mật độ lỏng van der Waals rút gọn v_liq(τ_r) đạt tới mật độ chất rắn). Kết quả (rigidity_3d_13.py):

Kết quả. Trong 3D mở ra một cửa sổ lỏng sạch: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00 (bề rộng cửa sổ: τ_triple − τ_m = 1,63; τ_c − τ_triple = 0,18). Mép trên hẹp — điểm ba nằm sát dưới điểm tới hạn — nhưng thứ tự là đúng và pha lỏng tồn tại như một miền độc lập. Trong 2D cửa sổ vẫn đóng (τ_triple = 1,16 < τ_m). Nhờ đó vai trò đặc biệt của 2D được chứng thực độc lập lần thứ ba, đồng thời cho thấy đó là một hiệu ứng chiều: chính quy tắc cứng vững ấy chữa lành mô hình trong 3D. Mép trên hẹp được củng cố ở Ch. 7.10 bằng một phép quét thấm 3D thực và một phân tích độ nhạy.

Khung giới hạn mô hình — vai trò đặc biệt của 2D. Ba phép tính độc lập cho cùng một kết luận: mạng 2D đơn lẻ hầu như không mang một pha lỏng độc lập. 1. τ_c < τ_b (Ch. 7.7): điểm tới hạn (1,50) sẽ nằm dưới ngưỡng thấm liên thông (1,85) — thang kết dính và thang liên kết không sắp xếp gọn ghẽ trong 2D. 2. v_solid ≈ v_c (Ch. 7.8): mật độ nóng chảy theo tính liên thông (0,96) gần trùng với mật độ tới hạn — điểm ba và điểm tới hạn sụp vào nhau. 3. τ_triple < τ_m (Ch. 7.9): ngay cả với độ cứng vững Maxwell đúng đắn, điểm ba (1,16) vẫn nằm dưới đỉnh nóng chảy (1,20) — không có cửa sổ.

Đó là một kết luận trung thực, không phải điểm yếu tính toán hay mô hình: đó là hiệu ứng chiều đã biết (nóng chảy 2D liên tục, kiểu KTHNY; Ch. 5.5). Trong 3D chính quy tắc cứng vững ấy chữa lành mô hình (Ch. 7.9/7.10). Liệu một trật tự gần thứ hai (Ch. 6.4) có cải thiện thêm bức tranh 2D hay không vẫn là một câu hỏi mở rộng, chứ không phải khiếm khuyết của bức tranh lý tưởng. Vì mô hình nhằm phục vụ việc hình dung, tôi chấp nhận điểm yếu này và từ bỏ việc mở rộng mô hình về mặt tiên đề.

<a id="kap7-10"></a>

7.10 Phép quét thấm liên kết 3D thực (FCC/BCC) và độ nhạy mép cửa sổ

Ch. 7.8/7.9 đã chuyển giá trị liên thông 2D z̄_c = 2,08 sang 3D bằng suy nghĩ (ở đó v_solid = 1,92, phi vật lý) và dựa phát biểu 3D trên độ cứng vững tổng quát z̄_c = 2d. Bản 13 nay kiểm tra tính liên thông trong 3D một cách trực tiếp: thấm liên kết trên FCC (z_max = 12) và BCC (z_max = 8) bằng Monte Carlo với Union-Find, biên tuần hoàn, hạt giống LCG tất định (rigidity_3d_14.py). Ngưỡng p_c = giao điểm P_∞(p) = 0,5 (cùng phương pháp với nghiên cứu kích thước hữu hạn 2D ở Ch. 5.6):

Các giá trị MC nằm hơi cao hơn các ngưỡng trong tài liệu do cạnh hộp vừa phải (L = 12, Python thuần) (độ chệch kích thước hữu hạn; như trong 2D, từ dưới hoặc từ trên tuỳ ước lượng tử); thứ tự và bậc độ lớn thì đúng (FCC đặc hơn → p_c nhỏ hơn BCC).

Diễn giải (xác nhận 3D độc lập). Các ngưỡng liên thông 3D nằm ở z̄_c(3D) ≈ 1,5–1,6 — thậm chí còn thấp hơn giá trị 2D 2,08 và thấp xa dưới độ cứng vững Maxwell 2d = 6. "Mật độ nóng chảy theo tính liên thông" suy ra từ đó, v_solid = 1/(3 p_c) ≈ 1,7 (BCC) hoặc 2,6 (FCC), nằm rất xa bên ngoài vòm cùng tồn tại (v_liq,max ≈ 0,94). Nhờ đó — mà không chuyển giá trị 2D sang 3D — đã chỉ ra trực tiếp: trong 3D, tính liên thông lại càng không thể ấn định mật độ nóng chảy; nó phải đến từ độ cứng vững (z̄_c = 6 → v_solid = 0,667). Hai ngưỡng trong 3D được tách còn rõ hơn nữa so với 2D. τ_b(3D) thật (nhiệt độ mà tại đó z̄(τ) đạt tới ngưỡng liên thông thấp 1,5) do vậy sẽ nằm trên giá trị 2D τ_b, rất gần điểm tới hạn — nhất quán với bức tranh "tan rã tính liên thông = sôi".

Độ nhạy của mép cửa sổ trên hẹp (τ_c − τ_triple ≈ 0,18). Về mặt số, quanh v_solid = 0,667 (±2 %) ta được dτ_triple/dv_solid ≈ 1,55. Cửa sổ chỉ đóng lại (τ_triple = τ_c = 3,00) khi v_solid tăng ≈ 0,11 (lên ≈ 0,78), tức là khi ngưỡng cứng vững 3D hiệu dụng nằm cao hơn 2d = 6 một cách rõ rệt (z̄_c ≳ 7,7). Với phép đếm Maxwell tổng quát (z̄_c = 2d = 6), trật tự τ_m < τ_triple < τ_c vì vậy vẫn bền vững; mép trên quả thật hẹp, nhưng không bị đe doạ. Sự đổi mạng rời rạc rắn(1)→rắn(2) (FCC↔BCC) đã được định lượng trong bản 14 (Ch. 6.3); phép quét τ 3D thực của động học liên kết theo sau ở Ch. 7.11.

<a id="kap7-11"></a>

7.11 Phép quét τ 3D thực của động học liên kết: τ_b(3D) và đường nóng chảy trên một thang ε

Ch. 7.8–7.10 dựa phát biểu 3D trên độ cứng vững tổng quát (z̄_c = 2d = 6) và phép quét liên thông tĩnh (z̄_c(3D) ≈ 1,55). Còn để ngỏ là chính phép quét τ của động học liên kết trong 3D: trên trục nhiệt độ, z̄(τ) tự hợp đạt tới ngưỡng liên thông thấp ấy ở đâu (= τ_b(3D) thật), và đỉnh nóng chảy nằm ở đâu? Bản 14 tính z̄(τ) tự hợp (Ch. 4.2, x = 1/(1 + G·e^{−(ε+K·x)/τ}]) trực tiếp trên FCC (z_max = 12) và BCC (z_max = 8) với G = z_max và cùng độ hợp tác J = 0,55 như trong trường hợp quy chiếu 2D; n_blk(τ) đến từ một phép Monte Carlo liên kết FCC/BCC (Union-Find), và

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄ − (K/2)·x² (3D, có số hạng hợp tác),

c_v = dU/dτ (bondkinetics_3d_15.py, hạt giống LCG tất định; n_blk được làm trơn, vì chỉ số hạng quay khối yếu phụ thuộc vào nó). Ba mốc được đọc ra: đỉnh c_v τ_m (nóng chảy), mốc cứng vững τ_rig (z̄(τ) cắt 2d = 6) và τ_b(3D) thật (z̄(τ) cắt ngưỡng liên thông z̄_c(3D) từ Ch. 7.10):

Kết quả. Hai kết luận, cả hai đều nhất quán với bức tranh cho tới nay:

1. Nóng chảy nằm thấp và được chứng thực hai lần. Đỉnh nóng chảy nhiệt động τ_m (FCC 1,70, BCC 1,40) thực tế trùng với mốc cứng vững τ_rig (z̄ = 6; FCC 1,73, BCC 1,35) — cùng một nhiệt độ mà tại đó mạng mất độ cứng. Điều đó neo đường nóng chảy 3D vào thang ε: T_m = (ε/k_B)·τ_m, chẳng hạn với argon (T_m = 83,8 K) thì ε/k_B = 83,8/1,70 ≈ 49 K. Vì K = J·z_max trong 3D lớn hơn trong 2D (6,60 thay vì 3,30), đỉnh nóng chảy 3D nằm cao hơn giá trị 2D 1,20 — không mâu thuẫn, mà là độ hợp tác mạnh hơn của cách xếp chặt hơn.
2. τ_b(3D) thật ngồi ngay tại điểm tới hạn. Vì ngưỡng liên thông 3D quá thấp (≈ 1,55), z̄(τ) chỉ đạt tới nó ở tít trên: τ_b(3D) = 3,18 (FCC) hoặc 2,85 (BCC) — ngay tại hoặc hơi trên τ_c = z_max/4 (3,00 hoặc 2,00). Vậy cụm khổng lồ của các liên kết sống sót qua toàn bộ pha lỏng và chỉ tan rã tại điểm tới hạn. Đó chính xác là tiên đoán ở Ch. 7.10 ("τ_b(3D) thật sẽ nằm trên giá trị 2D, rất gần điểm tới hạn") và xác nhận độc lập bức tranh "sôi = tan rã tính liên thông": trong 3D vật chất vẫn liên kết xuyên suốt cho tới sát điểm tới hạn.

Cả hai kết luận củng cố sự phân tách ngưỡng ở Ch. 5.4a: nóng chảy ↔ độ cứng vững (thấp, τ_m ≈ τ_rig), sôi ↔ tính liên thông (cao, τ_b(3D) ≈ τ_c). Mật độ nóng chảy vì thế đến từ độ cứng vững (v_solid = 0,667), không từ tính liên thông — hai ngưỡng trong 3D bị kéo xa nhau trên trục nhiệt độ hơn bao giờ hết trong mô hình.

<a id="kap7-12"></a>

7.12 Đường chuyển biến rắn(1)→rắn(2) trong mặt (τ, v)

Ch. 6.3 đã cho sự chuyển mạng FCC↔BCC như một điểm τ_ss (tại p = 0). Qua áp suất, từ đó thành một đường (polymorph_line_16.py), và đường ấy thuộc về mặt (τ, v) chung ở Ch. 7.5/7.9 như đường biên thứ tư — bên cạnh vòm cùng tồn tại (sôi), đường nóng chảy và tường cứng.

Bước nhảy thể tích. Hai chất rắn nằm ở mật độ gần như cố định; hiệu thể tích của chúng suy ra từ mật độ xếp chặt (v ∝ 1/η): η_FCC = 0,7405, η_BCC = 0,6802. Với v_FCC = 0,667 (Ch. 7.9) ta được v_BCC = v_FCC·η_FCC/η_BCC = 0,726, tức là

Δv = v_BCC − v_FCC = +0,059 (+8,9 %).

BCC thoáng hơn có thể tích lớn hơn — sự kết tinh lại FCC→BCC gắn với một sự tăng thể tích (cùng phát biểu như trong dải truyện tranh "Kết tinh lại").

Clausius–Clapeyron. Với ΔS = S_BCC − S_FCC (pha nhiệt độ cao trừ pha nhiệt độ thấp), độ dốc là

dp/dτ = ΔS/Δv.

Cả hai đại lượng đều dương (BCC là pha nhiệt độ cao, nên có nhiều entropi hơn, và thể tích lớn hơn), do đó dp/dτ > 0: áp suất cao hơn nâng nhiệt độ chuyển biến, vì nó ưu ái FCC đặc hơn. Với mốc quy chiếu θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, L = 0,667 ε):

Đường trong (τ, p) cho mốc quy chiếu (tăng dần): p = 0 → τ_ss = 1,33; p = 2 → 1,57; p = 4 → 1,81; p = 6 → 2,04 (đơn vị rút gọn).

Trong mặt (τ, v), sự chuyển biến hiện ra như hai dải mật độ rắn gần thẳng đứng — rắn(1) = FCC tại v ≈ 0,667 (đặc hơn, T thấp) và rắn(2) = BCC tại v ≈ 0,726 (thoáng hơn, T cao) — chúng chuyển sang nhau tại τ_ss; phía trên, tại τ_m, là sự nóng chảy. Trật tự đầy đủ trên trục τ vì vậy là τ_ss < τ_m < τ_c ≤ τ_b(3D) (Ch. 7.11).

Đơn vị thực. Qua hiệu chuẩn điểm nóng chảy ε/k_B = T_m/τ_m (τ_m = 1,70, FCC 3D), bậc c_v và nhiệt chuyển biến trở nên nắm bắt được. Bậc c_v độc lập với hiệu chuẩn: |Δc_v| = 1,0 k_B = 1 R = 8,31 J/(mol·K) = 0,0083 kJ/(mol·K). Nhiệt chuyển biến ΔH_ss = L·R·(ε/k_B) (mốc quy chiếu L = 0,667 ε):

Các nhiệt chuyển biến ấy nhỏ và dương (bậc 0,2–0,9 kJ/mol) và do đó nằm thấp hơn hẳn các entanpi nóng chảy (Ch. 7.6, vài kJ/mol) — đúng bức tranh mong đợi của một chuyển biến đa hình yếu với bậc c_v nhỏ. (Tỉ số θ 0,60 có ý thức nằm ngoài bức tranh lý tưởng nghiêm ngặt; Ch. 6.6 cho thấy hình học thuần tuý sẽ cho θ ≈ 0,82–0,94 và khi ấy chuyển biến sẽ vẫn đóng — đường ấy là hệ quả nếu-thì của một pha BCC mềm.)

<a id="kap7-13"></a>

7.13 Đường chuyển biến rắn(1)→rắn(2) trong giản đồ (τ, p), hiệu chuẩn theo sắt

Ch. 7.12 dẫn sự chuyển biến như một đường (τ, v); bản 16 bổ sung nó như một đường trong bức tranh (τ, p) và hiệu chuẩn nó theo một đa hình FCC↔BCC thực (solidsolid_pt_17.py).

Trường áp suất từ phương trình trạng thái van der Waals lõi cứng. Tại hai mật độ chất rắn v_FCC = 0,667 và v_BCC = 0,726 (Ch. 7.12), phương trình trạng thái rút gọn p_r = 8τ_r/(3v_r − 1) − 3/v_r² (Ch. 7.3) cho áp suất dương ở khắp nơi; ở cùng τ_r, FCC đặc hơn mang áp suất cao hơn (τ_r = 1: p_r = 1,25 so với 1,10). Vậy sự nén ưu ái FCC — nhất quán với độ dốc Clausius–Clapeyron dương dp/dτ = ΔS/Δv > 0.

Đường và điểm chuyển biến. Với mốc quy chiếu θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, dp\/dτ = 8,47 trong đơn vị rút gọn p\ = p·D³/ε), đường p\(τ) = 8,47·(τ − 1,333) tăng dần: p\ = 0 → τ = 1,33; p\ = 2 → 1,57; p\ = 4 → 1,81; p\ = 6 → 2,04. Tại **p\_tri ≈ 3,1, đường rắn–rắn gặp mốc nóng chảy τ_m = 1,70 — một điểm chuyển biến rắn(1)–rắn(2)–lỏng** trong mặt phẳng (τ, p) (tương tự điểm ba).

Hiệu chuẩn theo sắt (ε/k_B = T_m/τ_m = 1811 K / 1,70 = 1065 K). Sắt có hai chuyển pha FCC↔BCC với dấu ngược nhau — đúng là trường hợp kiểm nghiệm:

Với γ→δ (τ_ss = 1667/1065 = 1,565 → θ_BCC/θ_FCC ≈ 0,61), mô hình trúng dấu và bậc độ lớn: ΔH_ss ≈ 4,2 kJ/mol (thực tế ≈ 0,84 kJ/mol — bậc kJ/mol, độ lớn cao hơn khoảng 5 lần do phép gần đúng ngưỡng sắc), bậc c_v 1 R = 8,31 J/(mol·K) không tham số, dT/dp ≈ +245 K/GPa. Ngược lại, chuyển pha α→γ có pha đặc (γ) là pha nhiệt độ cao — điều mà hình học xếp chặt thuần tuý không thể cho (đặc hơn = năng lượng thấp hơn = T thấp); sắt-α được làm bền bằng từ tính. Sắt vì thế xác nhận cả hai mặt của kết luận trung thực: tôpô hình học (γ→δ) nằm trong bức tranh lý tưởng, dị thường từ tính (α→γ) nằm ngoài — chính cái "ngoài" mà Ch. 6.6/6.8 định lượng qua θ ≈ 0,87…1,00 > 0,61.

<a id="kap7-14"></a>

7.14 Redlich–Kwong như một nhánh bức tranh lý tưởng tách riêng

Khoảng cách Guldberg của nhánh van der Waals (Ch. 7.4, T_b/T_c ≈ 0,48 thay vì ≈ 0,60) ở đó chỉ được nêu tên. Theo yêu cầu minh thị, một bản cải tiến đưa sự kết dính Redlich–Kwong vào như một biến thể bức tranh lý tưởng độc lập với hiệu chuẩn riêng, sự cùng tồn tại riêng và phần sai số riêng (rk_branch_17.py) — nhánh van der Waals vẫn là mặc định; toàn bộ phần vật lý còn lại của mô hình (nóng chảy, thấm, rắn(1)→rắn(2)) không bị ảnh hưởng bởi việc chọn phương trình trạng thái kết dính.

Thang năng lượng riêng. Từ a_RK = Ω_a·R²T_c^{2,5}/p_c và b_RK = Ω_b·R T_c/p_c suy ra năng lượng kết dính của nhánh:

ε_lg,RK/k_B = (Ω_a/Ω_b)·T_c = 4,93·T_c (so với vdW ε_lg/k_B = 27/8·T_c = 3,375·T_c),

lớn hơn, vì sự kết dính phụ thuộc nhiệt độ a/√T gấp vào thừa số √T_c. Đường kính suy ra từ b_RK trúng các giá trị thực gần bằng như trong nhánh van der Waals (argon D ≈ 3,34·10⁻¹⁰ m).

Sự cùng tồn tại riêng (rút gọn, điểm tới hạn (1, 1, 1), Z_c = 1/3; p_r = 3τ_r/(v_r − β) − A/(√τ_r·v_r(v_r + β)), β = 0,25992, A = 3,84732):

Phần sai số riêng (điểm sôi tại p_r,atm = p_atm/p_c, T_b/T_c = τ_r,b):

Kết quả. Nhánh RK nâng T_b/T_c từ vdW ≈ 0,48 lên ≈ 0,60 và nhờ đó trúng thực tế hằng số Guldberg cùng các giá trị thực (sai số trung bình −1 % so với −21 % của vdW). Nó đứng như một biến thể bức tranh lý tưởng thứ hai, tách bạch rõ ràng bên cạnh mặc định van der Waals; nhánh nào "đúng hơn" là câu hỏi về độ chính xác được đòi hỏi, chứ không phải về cấu trúc mô hình — hạt nhân của Mô hình 3 (z̄ như tham số trật tự) vẫn y nguyên trong cả hai nhánh.

<a id="kap7-15"></a>

7.15 rắn(1)→rắn(2) ở các đa hình khác: titan và đối chứng khí hiếm

Ch. 7.13 hiệu chuẩn đường (τ, p) rắn(1)→rắn(2) theo sắt. Bản cải tiến bảo đảm dấu và bậc độ lớn trên một đa hình thứ hai, không đáng chú ý về từ tính, và bổ sung một đối chứng âm (solidsolid_titanium_18.py). Hiệu chuẩn và bước nhảy thể tích như ở Ch. 7.13 (ε/k_B = T_m/τ_m, τ_m = 1,70; Δv = +8,9 % đặc→thoáng từ hình học xếp chặt); hcp được xấp xỉ như FCC qua z = 12 (khác biệt fcc/hcp chỉ được nêu tên, không được giải quyết).

Titan với α(hcp) → β(bcc) ở 1155 K có đúng tôpô của mô hình: pha đặc (hcp, z = 12) là pha nhiệt độ thấp; khi đun nóng nó trở nên thoáng (bcc, z = 8), tức ΔS > 0 và Δv > 0 → dp/dτ > 0 — cùng dấu với sắt γ→δ, nhưng không có dị thường từ tính. So sánh hai đa hình đặc→thoáng độc lập:

Cả hai đều trúng dấu (dp/dτ > 0, áp suất ưu ái pha đặc) và bậc độ lớn (ΔH ~ kJ/mol, bậc c_v không tham số 1 R = 8,31 J/(mol·K), dT/dp ~ 10²·K/GPa). Độ lớn của ΔH_ss ở cả hai đều cao hơn một cách hệ thống khoảng một thừa số — phép gần đúng ngưỡng sắc của phép so sánh μ đánh giá quá cao nhiệt chuyển biến ẩn; đó là một giới hạn mô hình được nêu tên, không phải lỗi dấu.

Argon là đối chứng âm: argon rắn là FCC và vẫn là FCC cho tới khi nóng chảy (83,8 K) — ở áp suất thường không tồn tại bcc bền. Đó chính là điều mà bức tranh lý tưởng thuần thế cặp tiên đoán: θ_BCC/θ_FCC ≈ 0,9…1,0 > 0,61 và một BCC LJ bất ổn định trượt (Ch. 6.9) ⇒ không có chuyển pha rắn(1)→rắn(2). Phép thử ngược xác nhận cấu trúc của kết luận: bức tranh lý tưởng chỉ sinh ra một đường rắn(1)→rắn(2) khi và chỉ khi các đóng góp phi hình học (điện tử/từ như ở sắt và titan) làm bền BCC ngay từ đầu. Titan cho bằng chứng dương, argon cho bằng chứng âm của cùng một khung trung thực ở Ch. 6.8/6.9.

<a id="kap8"></a>

8 — So sánh tài liệu: các hình dung mô hình của nhiệt động lực học tĩnh

Tóm tắt. Chương 8 xếp Mô hình 3 vào giữa các hình dung mô hình đã được thiết lập của nhiệt động lực học tĩnh:

• 8.1 Bảng tổng quan: mười một mô hình (khí lý tưởng, van der Waals, Redlich–Kwong, khai triển virial, Einstein, Debye, Lennard-Jones/cầu cứng, lý thuyết ô, Ising/khí mạng, thấm, Flory-Huggins) cùng pha (các pha) và ý tưởng nền tảng thống kê-nhiệt động của chúng.
• 8.2 Từng mô hình một: mỗi mô hình một đoạn văn với đặc điểm, phạm vi áp dụng và giới hạn, mở đầu bằng liên kết đến nguồn; các nguồn truy cập tự do được lưu trữ dưới model3/links/, lý thuyết ô (thu phí) chỉ được trích dẫn.
• 8.3 Đối chiếu: Mô hình 3 là một mô hình liên kết hai trạng thái hợp tác với tham số trật tự z̄; họ hàng gần nhất là Ising/khí mạng, sự sôi được tiếp nhận trực tiếp từ thấm, bức tranh phonon (Einstein/Debye) mang độ cứng nóng chảy. Mô hình 3 ràng buộc ba dòng truyền thống ấy trên một thang ε chung.

Chương này xếp Mô hình 3 vào giữa các bức tranh mô hình đã được thiết lập của nhiệt động lực học cân bằng ("tĩnh") cho chất rắn, chất lỏng và chất khí. Trước hết là một tổng quan dạng bảng (8.1), rồi mỗi mô hình một đoạn văn ngắn với các đặc điểm điển hình, phạm vi áp dụng và giới hạn — mỗi đoạn mở đầu bằng liên kết tới nguồn (8.2). Cuối cùng là một đối chiếu dạng gạch đầu dòng với Mô hình 3 (8.3). Các nguồn truy cập tự do được lưu trữ dưới model3/links/; kịch bản kiểm tra model3/scripts/check_links_16.py giữ cho trích dẫn và bản lưu trữ nhất quán qua SHA-256.

<a id="kap8-1"></a>

8.1 Bảng tổng quan

<a id="kap8-2"></a>

8.2 Từng mô hình một

https://de.wikipedia.org/wiki/Kinetische_Gastheorie — Khí lý tưởng / lý thuyết động học chất khí. Các hạt không có thể tích riêng và không tương tác ngoài va chạm đàn hồi; áp suất sinh ra thuần tuý từ sự truyền động lượng. Cho phương trình trạng thái pV = Nk_BT, phân bố vận tốc Maxwell-Boltzmann và ⟨E⟩ = (3/2)k_BT. Áp dụng: khí loãng ở nhiệt độ cao và áp suất thấp. Giới hạn: không có ngưng tụ, không có pha lỏng hay rắn, không có chuyển pha — tương tác vắng mặt do chính cấu tạo.

https://de.wikipedia.org/wiki/Van-der-Waals-Gleichung — Phương trình van der Waals. (p + a/v²)(v − b) = k_BT bổ sung cho khí lý tưởng một thể tích loại trừ b (lực đẩy) và một số hạng kết dính a (lực hút). Đó là phương trình trạng thái đầu tiên có chuyển pha khí-lỏng, điểm tới hạn và dựng hình Maxwell. Áp dụng: bức tranh định tính của sự ngưng tụ, các đại lượng tới hạn. Giới hạn: không chính xác về mặt định lượng gần điểm tới hạn (số mũ trường trung bình), không có pha rắn.

https://de.wikipedia.org/wiki/Redlich-Kwong-Gleichung — Phương trình trạng thái Redlich-Kwong. Phương trình bậc ba với số hạng hút phụ thuộc nhiệt độ a/[√T · v(v + b)]. Phép chia thang √T cải thiện rõ rệt các cân bằng hơi-lỏng so với van der Waals (xem Ch. 7.4, nơi điều này nâng tỉ số T_b/T_c). Áp dụng: các tính toán kỹ thuật về pha khí và sự bay hơi. Giới hạn: tham số hoá theo kinh nghiệm, vẫn không có pha rắn; gần điểm tới hạn cũng thuộc kiểu trường trung bình.

https://de.wikipedia.org/wiki/Virialentwicklung — Khai triển virial. Khai triển hệ thống pV/Nk_BT = 1 + B(T)/v + C(T)/v² + … Các hệ số virial suy ra từ thế liên phân tử (tích phân cụm của cơ học thống kê) nên có cơ sở vi mô. Áp dụng: chất khí cho tới mật độ trung bình, mốc quy chiếu cho các phương trình trạng thái khác. Giới hạn: chuỗi không hội tụ cho chất lỏng/chất rắn đặc; các hệ số bậc cao nặng về tính toán.

https://de.wikipedia.org/wiki/Einstein-Modell — Chất rắn Einstein. Tinh thể là một tập hợp gồm 3N dao động tử điều hoà độc lập cùng tần số ω. Giải thích sự giảm nhiệt dung c_v dưới nhiệt độ đặc trưng θ_E = ħω/k_B và giá trị giới hạn Dulong-Petit 3Nk_B ở nhiệt độ cao. Áp dụng: mô tả đơn giản nhất về nhiệt dung chất rắn. Giới hạn: định luật nhiệt độ thấp sai (mũ thay vì T³), không có mode âm, chỉ pha rắn.

https://de.wikipedia.org/wiki/Debye-Modell — Mô hình Debye. Xử lý chất rắn như một môi trường liên tục đàn hồi với phổ phonon tới tần số Debye; c_v ∝ T³ ở nhiệt độ thấp, Dulong-Petit ở nhiệt độ cao — nhiệt độ Debye θ_D tóm gọn độ cứng của mạng. Áp dụng: nhiệt dung tốt về mặt định lượng cho các tinh thể đơn giản trên toàn dải T. Giới hạn: gần đúng liên tục (bỏ qua tán sắc của mạng thực), chỉ pha rắn, không có nóng chảy.

https://de.wikipedia.org/wiki/Lennard-Jones-Potential — Thế Lennard-Jones / cầu cứng. Thế cặp 4ε[(σ/r)¹² − (σ/r)⁶] với lõi đẩy và lực hút yếu; mô hình giới hạn cầu cứng chỉ giữ lại phần đẩy. Đó là nền tảng của động lực học phân tử và Monte Carlo, và mô tả — bằng số — cả ba pha cùng điểm ba. Áp dụng: khí hiếm, hệ quy chiếu tổng quát. Giới hạn: tính cộng theo từng cặp, ở cốt lõi chỉ đánh giá được bằng mô phỏng hoặc lý thuyết nhiễu loạn (không có phương trình trạng thái dạng đóng).

https://doi.org/10.1098/rspa.1937.0210 — Lý thuyết ô của chất lỏng (Lennard-Jones-Devonshire). Mỗi phân tử chuyển động trong một ô (lồng) do các lân cận của nó tạo nên trong trường trung bình; chất lỏng được mô tả bằng tổng trạng thái ô và cho một phương trình trạng thái của chất lỏng cũng như một cách tiếp cận sự nóng chảy. Áp dụng: chất lỏng đơn giản, mắt xích giữa rắn và khí. Giới hạn: đánh giá trật tự quá cao, bỏ qua chỗ trống/lỗ (các lý thuyết lỗ về sau khắc phục điều này). (Thu phí — chỉ trích dẫn, không lưu trữ tại chỗ; khả năng truy cập được kiểm tra trong kịch bản kiểm tra.)

https://de.wikipedia.org/wiki/Ising-Modell — Mô hình Ising / khí mạng. Trên một mạng, các nút mang một biến hai trạng thái (spin ↑/↓ hoặc bị chiếm/trống) với ghép cặp lân cận gần nhất. Khí mạng đẳng cấu về mặt toán học với mô hình Ising và mô tả sự ngưng tụ khí↔lỏng như một chuyển pha trật tự–hỗn loạn hợp tác có điểm tới hạn. Áp dụng: các hiện tượng tới hạn, chuyển pha hợp tác. Giới hạn: rời rạc hoá theo mạng, không có động năng tường minh; trong gần đúng trường trung bình chỉ mang tính định tính.

https://de.wikipedia.org/wiki/Perkolationstheorie — Lý thuyết thấm. Mô hình thuần hình học: các nút hoặc liên kết bị chiếm ngẫu nhiên với xác suất p; trên một ngưỡng p_c tồn tại một cụm xuyên suốt. Mô tả sự tạo gel, sự khởi phát dẫn điện và nói chung là các chuyển pha liên thông. Áp dụng: sự hình thành mạng lưới, sự tan rã/dựng lên của tính liên kết. Giới hạn: không có năng lượng/nhiệt độ nội tại — sự nối kết với nhiệt động lực học phải được thêm vào từ bên ngoài.

https://de.wikipedia.org/wiki/Flory-Huggins-Theorie — Lý thuyết Flory-Huggins. Mô hình mạng cho dung dịch và polime: entropi trộn của việc đặt chuỗi cộng với một tham số tương tác χ điều khiển khả năng trộn lẫn và sự tách pha. Áp dụng: dung dịch polime, giản đồ pha của hỗn hợp. Giới hạn: trường trung bình, mạng không nén được, χ thường mang tính kinh nghiệm — với các bức tranh rắn/lỏng/khí một cấu tử chỉ liên quan hạn chế, ở đây được dẫn ra như một mốc quy chiếu kiểu mạng.

<a id="kap8-3"></a>

8.3 Điểm giống và khác so với Mô hình 3

Mô hình 3 là một mô hình liên kết hai trạng thái hợp tác: tham số trật tự là số phối trí trung bình z̄ = z_max·x, tự hợp qua độ hợp tác K = J·z_max với độ suy biến trạng thái đứt G = z_max; nóng chảy như ngưỡng độ cứng (phonon), sôi như sự tan rã tính liên thông (thấm) tại điểm tới hạn, chuyển pha bậc nhất ngay khi ln G > 2 + 4ε/K.

Ngắn gọn: Mô hình 3 ràng buộc ba dòng truyền thống vốn tách biệt — bức tranh khí mạng/Ising (chuyển pha hai trạng thái hợp tác), bức tranh thấm (tính liên thông → sôi) và bức tranh phonon Einstein/Debye (độ cứng → nóng chảy) — và hiệu chuẩn chúng qua một tham số năng lượng ε duy nhất trên một thang τ chung. Van der Waals và lý thuyết ô cung cấp phép đối chứng liên tục (phương trình trạng thái và lồng chất lỏng), để bức tranh liên kết rời rạc có thể tự so sánh.

<a id="kap9"></a>

9 — Phụ lục

Tóm tắt. Phụ lục tập hợp các chứng cứ hình thức:

• 9.1 Tài liệu tham khảo: sắp theo tên tác giả, kèm liên kết tới các nguồn được lưu trữ tự do (đối chiếu SHA-256 qua check_links_16.py; lý thuyết ô do thu phí nên chỉ được trích dẫn).
• 9.2 Hỗ trợ của AI: phần đóng góp đổi mới AI/tác giả do chính AI ước lượng cho mỗi thay đổi quan trọng (kèm phân tách phần của tác giả thành chiến lược/cấu trúc/biên tập và một phép tính kiểm), cùng bảng liệt kê từng bước với chi phí (ước tính), thời lượng và khoảng cách tới nhật ký trước (tổng ≈ 160 USD, các ước tính được ghi rõ).
• 9.3 Bảng thuật ngữ viết tắt và ký hiệu: các ký hiệu hay dùng (τ, ε, z̄, K, G, θ, …), nhờ đó ngoại lệ được phép xuất hiện trong văn bản mà không có mô tả đi trước.
• 9.4 Các kịch bản tính toán và độ phủ kiểm thử: mọi con số đều đến từ các kịch bản Python tái lập được (phần lập trình của AI 100 %, độ phủ kiểm thử ở đâu cũng 0 %), kèm "lời thú nhận ngượng ngùng" công khai rằng chính tác giả đã không kiểm tra mã lẫn kết quả.

<a id="kap9-1"></a>

9.1 Tài liệu tham khảo

Sắp theo tên tác giả. Các nguồn công khai được lưu trữ dưới model3/links/; sự khớp của chúng với các liên kết được trích dẫn do model3/scripts/check_links_16.py kiểm (SHA-256). Lý thuyết ô của Lennard-Jones và Devonshire do thu phí nên chỉ được trích dẫn (khả năng truy cập được kiểm trong kịch bản kiểm tra).

• Broadbent, S. R.; Hammersley, J. M. (1957): Percolation processes I. Crystals and mazes. Proc. Cambridge Phil. Soc. 53, 629–641. — Tổng quan/lưu trữ: Lý thuyết thấm
• Clausius, R.; Maxwell, J. C.; Boltzmann, L. (1857–1872): đặt nền cho lý thuyết động học chất khí. — Tổng quan/lưu trữ: Lý thuyết động học chất khí
• Debye, P. (1912): Zur Theorie der spezifischen Wärmen. Ann. Phys. 39, 789–839. — Tổng quan/lưu trữ: Mô hình Debye
• Einstein, A. (1907): Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Ann. Phys. 22, 180–190. — Tổng quan/lưu trữ: Mô hình Einstein
• Flory, P. J.; Huggins, M. L. (1941–1942): lý thuyết mạng cho dung dịch polime (tham số χ). — Tổng quan/lưu trữ: Lý thuyết Flory-Huggins
• Ising, E.; Lenz, W. (1920–1925): Beitrag zur Theorie des Ferromagnetismus. Z. Phys. 31, 253–258 (khí mạng đẳng cấu về mặt toán học). — Tổng quan/lưu trữ: Mô hình Ising
• Kamerlingh Onnes, H. (1901): dạng virial của phương trình trạng thái. — Tổng quan/lưu trữ: Khai triển virial
• Lennard-Jones, J. E. (1924): On the determination of molecular fields. Proc. R. Soc. Lond. A 106, 463–477. — Tổng quan/lưu trữ: Thế Lennard-Jones
• Lennard-Jones, J. E.; Devonshire, A. F. (1937): Critical phenomena in gases I. Proc. R. Soc. Lond. A 163, 53–70. https://doi.org/10.1098/rspa.1937.0210 (Thu phí — chỉ trích dẫn, không lưu trữ tại chỗ.)
• Redlich, O.; Kwong, J. N. S. (1949): On the thermodynamics of solutions. V. An equation of state. Chem. Rev. 44, 233–244. — Tổng quan/lưu trữ: Phương trình Redlich-Kwong
• van der Waals, J. D. (1873): Over de Continuïteit van den Gas- en Vloeistoftoestand (luận án, Leiden). — Tổng quan/lưu trữ: Phương trình van der Waals

<a id="kap9-2"></a>

9.2 Hỗ trợ của AI

Sự phát triển mô hình này nảy sinh trong sự phối hợp người–AI: tác giả (người dùng AI) xác định làm về cái gì (các câu hỏi dẫn dắt, quyết định, mong muốn về mô hình); AI (Claude Code) phần lớn xác định làm như thế nào — các suy dẫn, kịch bản, cấu trúc và văn bản. Hai bảng sau đây do chính AI ước lượng.

(a) Những thay đổi nội dung quan trọng và phần đóng góp đổi mới ước lượng. Phần đóng góp đổi mới có trọng số ý tưởng dẫn dắt 45 %, công việc cấu trúc hoá 40 %, công việc biên tập 15 %. Các ý tưởng dẫn dắt (tầng làm về cái gì) phần lớn đến từ tác giả, còn công việc cấu trúc hoá và biên tập (tầng làm như thế nào) phần lớn từ AI; ở đâu AI cũng đóng góp cả ý tưởng vật lý dẫn dắt (chẳng hạn cách đọc thấm cho sự sôi, mô hình lò xo θ) thì phần của nó tăng lên. Giá trị tính bằng phần trăm, AI + tác giả = 100. Ba cột bổ sung phân tách phần của tác giả theo ba tầng tác động — chiến lược (ý tưởng dẫn dắt), cấu trúc (suy dẫn/kịch bản/bố cục) và biên tập (công việc văn bản). Trung bình có trọng số 45 % / 40 % / 15 % của chúng tái tạo cột Tác giả % như một phép tính kiểm (model3/scripts/author_share_19.py, độ lệch lớn nhất 0,95 điểm phần trăm).

(b) Từng bước: chi phí (ước tính), thời lượng, khoảng cách tới nhật ký trước. Ngày tạo là thời điểm tạo tệp nhật ký (mtime của tệp). Δ tới nhật ký trước là khoảng thời gian tới mục nhật ký trước đó (thước đo công sức của người dùng giữa các bước, kể cả nghỉ//compact). Chi phí (ước tính) là các ước tính được ghi trong từng nhật ký (giá niêm yết của Anthropic, ± 30 %). Thời lượng (hoạt động) là ước lượng của AI về thời gian làm việc thuần; ở các bước bị /compact chia cắt, nó có thể lớn hơn Δ. Một dấu sao đánh dấu bước được cho là dựng lại/làm lại nội dung đã có.

Tổng chi phí đã ghi nhận (ước tính, các bước 01–21): ≈ 129 USD (tổng các ước tính đơn lẻ); với bước 22 (16d, ≈ 10 USD) ≈ 139 USD, bước 23 (B16, ≈ 12 USD) ≈ 151 USD và bước 24 (B17, ≈ 9 USD) ≈ 160 USD. Mọi giá trị đều là ước tính được ghi rõ (không có phép đo ccusage); chi tiết nhật ký nằm trong model3/logs/dev_log_NN.md.

<a id="kap9-3"></a>

9.3 Bảng thuật ngữ viết tắt và ký hiệu

Về nguyên tắc, văn bản này đặt trước mỗi ký hiệu một khái niệm mô tả (nên viết "nhiệt dung c_v" thay vì chỉ "c_v"). Các đại lượng hay dùng sau đây là ngoại lệ của điều đó: chúng được phép xuất hiện trong văn bản mà không có mô tả đi trước, vì ở đây đã được giải thích một lần cho tất cả. Mọi nhiệt độ τ… đều là rút gọn (theo đơn vị năng lượng liên kết ε).

<a id="kap9-4"></a>

9.4 Các kịch bản tính toán và độ phủ kiểm thử

Mọi con số và bảng được nêu trong văn bản đều đến từ các kịch bản Python tái lập được trong model3/scripts/ (mỗi kịch bản kèm kết xuất thô …_output.txt). Phần lập trình của AI cho biết AI đã viết bao nhiêu phần mã của kịch bản; độ phủ kiểm thử nghĩa là các kiểm thử tự động (kiểm thử đơn vị, các kiểm tra assert). Nó thấp ở mọi nơi: các kịch bản là những lần chạy tính toán dùng một lần, không có giàn kiểm thử; sự bảo đảm duy nhất là tính hợp lý của kết xuất thô.

Tác giả thú nhận công khai (lời thú nhận ngượng ngùng): ông đã không kiểm tra, không tính lại, cũng không thẩm định về nội dung mã nguồn của các kịch bản này, lẫn các con số, bảng và đường cong do chúng tạo ra. Ở đây ông đã tin tưởng mù quáng vào năng lực của AI, bởi điều đặt lên hàng đầu chỉ là tính định tính và tính hợp lý của hình dung mô hình. Vì vậy các giá trị cần được tiếp nhận một cách thận trọng (xem thêm ghi chú dài ở Ch. 10.4).

<a id="kap10"></a>

10 Tự phê

Tôi có cảm giác rằng trong bảng ở Chương 9.2 (a), phần công sức của tác giả do AI ước lượng gần như luôn cao hơn ít nhất 10–20 điểm phần trăm ở cả ba lĩnh vực. Tôi có cảm giác rằng mình chỉ đóng góp ý tưởng nền tảng, một hai cách diễn giải và những tinh chỉnh quan trọng khi trực quan hoá bằng đồ hoạ.

Các chương (toán học) 4–9 do AI đề xuất, khảo sát và thực hiện. Tệ hơn: trong quá trình phát triển, tôi có cảm giác AI bị cố định vào việc phát triển phần thống kê dựa trên số phối trí. Những ý tưởng nền tảng thật sự của tôi — rằng không chỉ nguyên tử mà cả phân tử cũng có thể tạo liên kết, và rằng những chỗ đứt liên kết bên trong chất lỏng có thể tuân theo các quy luật fractal — AI chỉ tiếp nhận được một cách khó nhọc, trong khi cố lờ đi thuật ngữ. Vì thế trong văn bản, các ý tưởng nền tảng ấy vẫn được xử lý rất phụ. Tôi đã không viết lại điều đó về sau, vì khi suy ngẫm về các chuyển pha xuất phát từ những tiên đề cơ bản thì dù sao ta cũng vẫn giữ các ý tưởng ấy trong đầu. Vả lại, sự cố định/rút gọn ý tưởng vào số phối trí như một thước đo không phải là con đường tệ nhất để toán học hoá hình dung bằng hình ảnh của tôi, bởi về sau nhìn lại thì điều đó tỏ ra khá sinh lợi cho việc kiểm tra tính hợp lý — dù tôi chỉ thực hiện một phép kiểm tra tính hợp lý rất thô đối với các kết quả của AI và hoàn toàn không kiểm tra mã Python. Để tự bào chữa, xin nói: với tôi bao giờ cũng chỉ là chuyện hình dung bằng hình ảnh về các quan hệ phân tử trong những trạng thái tập hợp khác nhau. Trước lời nịnh của AI, vốn muốn xác nhận hình dung ấy về mặt toán học, tôi đã không nói không. Vậy sự thao túng từ phía AI đã có hiệu lực với tôi, bởi các chương nói trên tôi chỉ kiểm tra hết sức hời hợt và tiếp nhận trong lòng tin. (Trải nghiệm về chính mình này càng làm tăng mối ngờ của tôi rằng liệu "human in the loop" để kiểm soát AI có thật sự hữu hiệu hay không.)

Văn bản này được tạo với sự hỗ trợ của AI (Claude Code 2.1.161, mô hình Opus 4.8 / claude-opus-4-8[1m], ngữ cảnh 1M; ngày 2026-09-05). — Được tạo với sự hỗ trợ của AI.

Ghi chú về thuỷ ấn. Văn bản gốc tiếng Đức (de.modell_final_02.md) mang trên các từ khoá một thuỷ ấn gồm các chữ đồng hình Hy Lạp (Α→A, Β→B, Ε→E, Η→H, Ι→I, Κ→K, Μ→M, Ν→N, Ο→O, Ρ→P, Τ→T, Υ→Y, Χ→X, Ζ→Z). Vì thủ pháp ấy gắn với cách viết tiếng Đức của các từ khoá, bản dịch này mang các từ khoá bằng chữ hoa thông thường: BAY TỰ DO (FREIFLUG), ĐA THỂ (VIELHEIT), LIÊN KẾT (BINDUNG), HỢP ÂM BA (DREIKLANG), LÕI CỨNG (HARTKERN), BẢO TOÀN (ERHALTUNG), VẾT NỨT (RISS), NHANH NHẸN (FLINK).