Termodinamika statistik untuk zat murni ideal dan wujud-wujudnya

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

Ringkasan menyeluruh

Model 3 adalah model statistik-mekanik yang kualitatif untuk suatu zat ideal, yang menurunkan ketiga wujud zat — padat, cair, dan gas — beserta transisinya dari satu gambaran ikatan tunggal dengan ikatan molekul. Asas dasarnya terletak pada diperhitungkannya tiga penyimpan energi bagi energi dalam: energi getaran dalam ikatan molekul, energi translasi molekul, dan energi gerak kooperatif dalam rotasi misel sebagai bentuk ideal dari aliran zat internal untuk menyetarakan fluktuasi kerapatan di dalam. Besaran dasar penting dalam kaitan ini adalah bilangan koordinasi rata-rata z̄ = z_max·x, yang menata dirinya secara swakonsisten melalui ikatan dua keadaan yang kooperatif (bertipe Ising, kooperativitas K = J·z_max, degenerasi keadaan patah G = z_max). Darinya menyusul tiga transisi pada satu skala τ bersama yang dikalibrasi oleh satu-satunya parameter energi ε: peleburan sebagai ambang fonon/kekakuan (z̄ mencapai rigiditas), pendidihan sebagai runtuhnya keterhubungan (perkolasi) di titik kritis, dan sebuah transisi opsional padat(1)→padat(2) (FCC↔BCC) dengan tangga c_v tanpa parameter sebesar 1,0 k_B. Sebuah transisi menjadi orde pertama begitu ln G > 2 + 4ε/K. Di sini degenerasi keadaan patah G adalah jumlah kontak tetangga yang dapat diputus (G = z_max, bilangan koordinasi maksimum), energi ikatan ε adalah energi yang diperoleh per kontak tetangga, dan kooperativitas K = J·z_max adalah hasil kali kekuatan kopling bertipe Ising J dengan z_max; jadi transisi makin cepat menjadi orde pertama makin besar keragaman kontak G dan makin kuat kooperativitas K terhadap energi ikatan ε.

Artikel ini menuntun dari tujuh aksioma (Bab 1) melalui analisis kontradiksi (Bab 2) dan sebuah pemerian bergambar dengan strip komik statis (Bab 3) menuju kerangka statistik-mekanik — ikatan kooperatif, fungsi partisi, dan kapasitas kalor (Bab 4–5) —, ke kemasukakalan matematis 3D (Bab 6), dan akhirnya ke kalibrasi serta persamaan keadaan (Bab 7). Pemikiran itu dilengkapi dengan pemindaian c_v 3D penuh untuk peleburan orde pertama (Δz̄, kalor laten ΔU, pergeseran puncak), sebuah model pegas Einstein untuk kekakuan fonon θ(kisi), dan garis transformasi padat(1)→padat(2) melalui Clausius–Clapeyron. Selanjutnya θ(kisi) dihaluskan menjadi DOS fonon bertipe Debye, lalu transformasi padat(1)→padat(2) diperkenalkan pula sebagai garis dalam gambaran (τ, p) (dikalibrasi pada besi), dan ditambahkan sebuah cabang Redlich–Kwong terpisah di samping bawaan van der Waals. Di situ θ(kisi) diuji pada sebuah kisi penuh Lennard-Jones (FCC stabil, BCC ideal tak stabil terhadap geser), dan topologi padat(1)→padat(2) diuji pada titanium (penegas kedua) dan argon (kontrol negatif). Bab 8 menempatkan model ini di antara sebelas gambaran model mapan dalam termodinamika statis. Lampiran mendokumentasikan sumber pustaka (diurutkan menurut penulis) dan dukungan AI (perkiraan pangsa inovasi, perkiraan biaya, dan lama tiap langkah pengembangan); bagian bantuan menghimpun cara baca, butir kritis yang tersisa, dan pandangan ke depan.

Pernyataan inti. Satu ukuran ikatan kooperatif z̄ mengaitkan tiga garis tradisi yang biasanya terpisah — gambaran gas kisi/Ising (transisi dua keadaan kooperatif), gambaran perkolasi (keterhubungan → pendidihan), dan gambaran fonon Einstein/Debye (kekakuan → peleburan) — menjadi satu gambaran fase yang konsisten dan lengkap secara kualitatif, dengan sengaja di dalam sebuah gambaran ideal bertipe van der Waals, sementara bangunan teorinya secara wajar dapat diturunkan dari sebuah gambaran dinamis yang mudah dibayangkan. Bukan ketepatan, melainkan sifat bergambarlah kekuatan sejati model ini, dan itulah yang menandai nilai ilmiahnya sebagai model.

Daftar isi

• Ringkasan menyeluruh
• 1 — Aksioma dengan nama kata kunci (dikonsolidasikan)
• 2 — Analisis kontradiksi dan kesenjangan
• 3 — Pemerian bergambar
• 3.1 Ketiga fase tunggal
• 3.2 Transisi fase
• 3.3 Mengapa laju bunyi dalam zat cair turun bersama T
• 3.4 Bagaimana gerak Brown timbul — dengan pengertian misel yang dipertajam
• 3.5 Mengapa reaksi di lapisan batas misel menjadi lebih cepat
• 3.6 Jendela cair tiga dimensi
• 4 — Model statistik-mekanik
• 4.1 Derajat kebebasan dan peubah
• 4.2 Ikatan dua keadaan kooperatif (bertipe Ising, parameter J)
• 4.3 Neraca energi — dengan suku kooperatif
• 5 — Fungsi partisi dan kapasitas kalor
• 5.1 Fungsi partisi terfaktorkan
• 5.2 Kapasitas kalor per fase (analitik)
• 5.3 Kurva c_v(τ) numerik **dengan** suku kooperatif — kalor laten dikuantifikasi
• 5.4 Ambang perkolasi keterhubungan sebagai hasil (ditegaskan tanpa perubahan) - 5.4.a Dua ambang berbeda dipisahkan dengan bersih
• 5.5 Orde transisi peleburan
• 5.6 Kajian ukuran hingga untuk ambang perkolasi
• 5.7 Pemindaian c_v(τ) 3D penuh untuk peleburan orde pertama
• 6 — Kemasukakalan matematis tafsiran 3D
• 6.1 Perubahan derajat kebebasan rotasi 2D → 3D
• 6.2 Pengaruh berbagai kisi kristal — keterhubungan sebagai ciri pembeda
• 6.3 Transisi tambahan padat(1) → padat(2)
• 6.4 Renungan: anomali bunyi pada air cair
• 6.5 Alternatif: medan pusaran fraktal alih-alih rotasi misel
• 6.6 Kekakuan fonon mikroskopis θ(kisi): sebuah model pegas Einstein
• 6.7 Taksiran matematis kemetastabilan (2D triangular ↔ persegi, 16b)
• 6.8 Penghalusan θ(kisi) bertipe Debye
• 6.9 θ(kisi) dari sebuah kisi penuh Lennard-Jones
• 7 — Kalibrasi dan persamaan keadaan
• 7.1 Metode kalibrasi (satu parameter energi ε) — disesuaikan
• 7.2 Tabel kalibrasi
• 7.3 Sumbu kerapatan/volume: persamaan keadaan van der Waals berinti keras
• 7.4 Kalibrasi c/g kuantitatif: a → ε, b → D²
• 7.5 Permukaan (τ, v) bersama dengan titik tripel
• 7.6 Kalor laten terhadap entalpi peleburan
• 7.7 Penyeragaman skala energi: ε ↔ ε_lg
• 7.8 Titik tripel tanpa parameter dari z̄_c
• 7.9 Perhitungan kontrol 3D: jendela cair yang bersih pada satu skala ε
• 7.10 Pemindaian perkolasi ikatan 3D sejati (FCC/BCC) dan kepekaan tepi jendela
• 7.11 Pemindaian τ 3D sejati untuk kinetika ikatan: τ_b(3D) dan garis lebur pada satu skala ε
• 7.12 Garis transformasi padat(1)→padat(2) pada permukaan (τ, v)
• 7.13 Garis transformasi padat(1)→padat(2) dalam diagram (τ, p), dikalibrasi pada besi
• 7.14 Redlich–Kwong sebagai cabang gambaran ideal yang terpisah
• 7.15 padat(1)→padat(2) pada polimorf lain: titanium dan kontrol gas mulia
• 8 — Perbandingan pustaka: gambaran model termodinamika statis
• 8.1 Tabel ikhtisar
• 8.2 Model satu per satu
• 8.3 Kesamaan dan perbedaan dengan Model 3
• 9 — Lampiran
• 9.1 Pustaka
• 9.2 Dukungan AI
• 9.3 Glosarium singkatan dan lambang
• 9.4 Skrip perhitungan dan cakupan uji

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

1 — Aksioma dengan nama kata kunci (dikonsolidasikan)

Ringkasan. Tujuh aksioma yang diberi nama kata kunci menetapkan zat ideal: jumlah partikel yang kekal, gerak yang murni klasik, sebuah ikatan dua keadaan yang kooperatif, dan bilangan koordinasi rata-rata z̄ sebagai parameter keteraturan. Semuanya menopang seluruh penurunan berikutnya; dalam naskah asli Jerman kata-kata kunci itu ditandai oleh tanda air homoglif (kunci di akhir teks).

Aksioma yang dipertegas:

1. TERBANG BEBAS (aksioma tambahan). Antara dua molekul yang terbang berpapasan dalam keadaan tak terikat tidak bekerja gesekan maupun gaya elektrostatik. Momentum dan energi hanya berubah lewat tolakan keras INTI KERAS saat bertumbukan dan lewat terjalin/putusnya IKATAN. Aksioma ini penting semata-mata dari sudut mekanika klasik dan neraca energi: ia menjamin bahwa tidak ada kanal disipatif di aras molekul — energi gerak tak pernah "hilang", melainkan hanya beralih ke ketiga penyimpan TRINADA atau ke energi ikatan. (Gesekan/kekentalan makroskopis timbul secara emergen dari puzel geser dan RETAK, bukan sebagai gaya dasar mikroskopis.)
2. KEMAJEMUKAN. Sebuah zat ideal terdiri atas sangat banyak molekul sejenis; sifat makroskopisnya diperikan secara statistik dalam limit termodinamis (N → ∞).
3. IKATAN. Setiap dua molekul yang bertetangga secara geometris dapat memikul sebuah ikatan. Ikatan yang utuh menurunkan energi sebesar nilai hingga ε; pemutusannya berbiaya tepat ε. Ikatan adalah sistem dua keadaan {utuh, patah}.
4. TRINADA. Energi dalam terbagi secara klasik ke tiga penyimpan: (a) translasi molekul, (b) getaran dalam ikatan yang utuh (pegas harmonik per ikatan), (c) rotasi blok molekul kaku (derajat kebebasan putar kolektif dari daerah yang terhubung).
5. INTI KERAS. Bila jarak antarmolekul turun di bawah nilai minimum σ, bekerja tolakan keras (volume sendiri). Bersama IKATAN diperoleh potensial pasangan efektif (kulit keras + sumur ikatan hingga).
6. KEKEKALAN. Momentum total dan energi total kekal. Sebuah molekul tunggal netral secara energetik — ia tidak menyangga energi di dalam dirinya; energi hidup dalam ikatan (IKATAN) dan dalam derajat kebebasan kinetik (TRINADA).
7. RETAK. Ikatan terurai/terbentuk menurut asas retak/pulih: makin banyak ikatan tetangga yang hilang → makin tinggi peluang patah; makin banyak tetangga yang terikat → makin tinggi peluang pulih. Kinetikanya kooperatif.
8. GESIT. Pematahan dan pemulihan berlangsung cepat dan beraktivasi rendah; subsistem ikatan berada dalam kesetimbangan lokal pada setiap suhu dan boleh dipisahkan secara adiabatik dari gerak molekul (dasar bagi fungsi partisi terfaktorkan).

Sikap pemodelan. Koordinasi geometris tetap dipertahankan secara ideal; yang berubah hanyalah bilangan koordinasi z akibat pematahan. Gerak timbul seperti pada puzel geser: rongga (dari volume sendiri INTI KERAS dan pematahan RETAK) berpindah, blok bergeser mengikutinya. Gambaran yang sengaja "salah" tentang kristal kecil yang berputar dan mengembara tetap dipertahankan sebagai landasan yang mudah dibayangkan.

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

2 — Analisis kontradiksi dan kesenjangan

Ringkasan. Analisis ini menguji model secara sistematis terhadap kontradiksi dalam dan kesenjangan serta melanjutkan catatan taraf pengerjaan; kelemahan yang dikenali disebutkan dan pada bab-bab berikutnya diuraikan atau ditandai dengan jujur sebagai batas model.

Daftar butir kritik

1. Perhitungan ganda translasi (aksioma 3a lawan 3c). Naskah kasar menyebut translasi dua kali (translasi molekul dan translasi blok). Dua penyimpan translasi yang saling bebas tidak ada. Penutupan: translasi blok hanyalah ragam kolektif dari translasi molekul yang sudah dihitung; yang sungguh baru hanyalah derajat kebebasan rotasi blok. Dalam neraca derajat kebebasan (Bab 4) translasi dihitung tepat satu kali.
2. "Netral secara energetik" (aksioma 5) lawan "ikatan menyerap energi" (aksioma 2/3b). Kontradiksi semu. Penutupan: kenetralan itu menyangkut molekul sebagai benda; energi duduk di antara molekul (dalam IKATAN) dan dalam ragam kinetik kolektif, tak pernah dalam penyangga intramolekul. Molekul ≠ ikatan — dengan demikian kontradiksi terurai.
3. Tarikan yang hilang (aksioma 4 hanya menyebut tolakan). Tanpa tarikan tak akan ada fase mampat. Penutupan: tarikan disediakan IKATAN (sumur berkedalaman ε), tolakan oleh INTI KERAS (kulit keras di σ). Bersama-sama: sebuah potensial pasangan efektif.
4. Umpan balik tanpa suku lawan (aksioma 6). "Makin banyak tetangga hilang → makin banyak patah" sendirian adalah umpan balik positif (retak merambat sampai runtuh total). Penutupan: suku pemulihan (makin banyak tetangga → makin banyak pemulihan) ditambah faktor Boltzmann dari energi ikatan hingga ε (IKATAN) menyeimbangkan RETAK dan menghasilkan transisi yang terdefinisi alih-alih suatu pelarian.
5. Pemisahan skala waktu yang tak terucap (aksioma 7). "Cepat dan beraktivasi rendah" baru berguna bila darinya dibuat pemisahan adiabatik antara dinamika ikatan dan translasi. Penutupan: dicatat secara eksplisit sebagai prasyarat pemfaktoran Z = Z_trans·Z_vib·Z_rot·Z_bond.
6. Geometri dan mekanisme gerak yang hilang. Aksioma tidak menyebut kisi koordinasi maupun apa yang memungkinkan aliran. Penutupan (sebagai andaian tambahan yang ditandai): geometri koordinasi acuan yang tetap dengan pengisian ikatan yang berubah-ubah; sebagai kasus 2D kanonik dipakai koordinasi segitiga (z_max = 6). Kegerakan lewat mekanisme rongga puzel geser. Kedua butir ini adalah satu-satunya tambahan sejati di luar ketujuh aksioma dan ditandai pada Bab 6 sebagai pertanyaan strategis.
7. Celah disipasi terpendam. Versi 07 membiarkan terbuka apakah molekul tak terikat kehilangan energi saat berpapasan. Penutupan lewat TERBANG BEBAS: tidak ada gesekan/elektrostatika mikroskopis; seluruh aliran energi dan momentum berjalan semata lewat tumbukan INTI KERAS dan kinetika IKATAN. Dengan itu neraca energi (Bab 4.3) tertutup dan statistik Maxwell–Boltzmann dalam kasus batas gas menjadi eksak (gas ideal bebas di antara tumbukan lenting).
8. Ketidakkonsistenan energetik parameter keteraturan. Versi 08 memakai suku kooperatif −(K/2)x² memang dalam energi bebas g(x) (untuk menghasilkan x(τ)), tetapi tidak dalam energi dalam U/N. Itu tak konsisten secara termodinamis: suku medan molekul yang sama harus muncul di keduanya. Penutupan (Bab 4.3): −(K/2)x² disusulkan ke U/N; kalor laten karenanya meningkat secara terukur (Bab 5.3) tanpa mengubah plato padat dan gas.
9. Kalibrasi c/g hanya kualitatif. Versi 09 hanya memberi sumbu kerapatan (persamaan keadaan van der Waals) dalam bentuk tereduksi. Penutupan (Bab 7.4): lewat titik kritis (T_c, p_c) dihitung garis didih mutlak dan T_b(p_atm) dibandingkan dengan nilai ukur; penyimpangan sistematis (kaidah Guldberg) disebutkan.
10. Sumbu terpisah tanpa titik tripel. Sumbu lebur dan sumbu didih adalah perhitungan terpisah. Penutupan (Bab 7.5): keduanya disatukan dalam satu bidang (τ, v), titik tripel dikonstruksi dan titik kritis ditandai.

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

3 — Pemerian bergambar

Ringkasan. Bab 3 menerangkan padat, cair, dan gas beserta transisinya secara mudah dibayangkan (dengan strip komik statis): ketiga fase tunggal (3.1), transisi fase rekristalisasi/peleburan/ pendidihan/sublimasi beserta titik tripel dan titik kritis (3.2), anomali bunyi (3.3), gerak Brown dan pengertian misel yang dipertajam (3.4–3.5), serta jendela cair tiga dimensi (3.6). Berdasarkan pemerian itu dibuat bagan (strip komik) yang hendak membantu membayangkan hubungan molekul yang dinamis dalam berbagai fase di depan mata batin. Keterbayangan inilah inti sejati model tersebut. Ia menjadi bernilai secara kualitatif dan ilmiah bila lewat mekanika statistik dapat ditunjukkan secara masuk akal bahwa model itu membuat kecenderungan/pengamatan empiris yang kualitatif dapat diramalkan bahkan tanpa matematika besar.

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

3.1 Ketiga fase tunggal

Zat padat. Hampir semua tetangga terikat satu sama lain lewat ikatan molekul (IKATAN), rata-rata bilangan koordinasi molekul z terletak dekat maksimum geometris. Jaring molekul adalah satu blok kaku tanpa rongga (INTI KERAS terkemas rapat) — puzel geser internal seperti pada zat cair tidak mungkin dalam versi padatnya. Satu-satunya gerak: getaran di sekitar posisi setimbang (getaran ikatan, TRINADA-b), yang ditentukan oleh ikatan molekul dengan para tetangga. Momentum/energi berayun lenting dalam jaring (KEKEKALAN, TERBANG BEBAS: tanpa rugi). Pada kasus batas (suhu sangat rendah, kemasan ideal) hadir keadaan ideal dari kristal ideal.

Zat cair. RETAK kooperatif telah menembusi jaring: ikatan hilang setempat-setempat, rata-rata bilangan koordinasi molekul z telah turun, blok kaku terpecah di dalam menjadi banyak blok kecil. Dari monolit menjadi timbunan kerikil. Blok-blok itu bertranslasi dan berotasi (TRINADA-c). Di antaranya terbuka rongga, yang lewatnya mereka bergeser saling melewati — puzel geser gerak Brown aktif. Geometri koordinasi setempat tetap dipertahankan secara ideal; yang berubah hanyalah berapa banyak ikatan yang benar-benar ada dan terbentuk.

Gas. RETAK menang: ikatan patah hampir di mana-mana (z → 0). Dalam kasus ekstrem tak satu molekul pun masih punya ikatan molekul ke molekul lain. Blok-blok telah terurai menjadi molekul tunggal. Mereka terbang bebas, terkopel hanya oleh tumbukan INTI KERAS (di antaranya TERBANG BEBAS); lajunya mengikuti Maxwell–Boltzmann (KEKEKALAN). Getaran ikatan dan rotasi blok padam — translasi murni (TRINADA-a). Pada kasus batas (suhu tinggi, tekanan rendah) hadir keadaan ideal dari gas ideal.

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

3.2 Transisi fase

Rekristalisasi (padat(1)/padat(2)). Masih di bawah peleburan, zat padat dapat berganti kisi tanpa kehilangan kekakuannya: dua kemasan ideal dengan bilangan koordinasi berbeda (FCC z_max = 12 ↔ BCC z_max = 8) memiliki entalpi bebas μ(τ, p) yang berbeda; pada μ_FCC = μ_BCC geometri koordinasi berganti secara melompat (bukan hanya pengisian z̄ seperti pada peleburan). Yang tampak adalah tangga c_v kecil (Δc_v = perubahan z_max/4 = 1,0 k_B) ditambah kalor transformasi laten yang lemah. Klaim AI: hal itu hanya terbayangkan dalam 3D (keterhubungan diskret dan terkuantisasi; Bab 6.3). Dalam kasus dua dimensi, rekristalisasi sejati — yaitu transisi antara dua fase padat — lebih mudah dibayangkan. Di bidang datar hanya ada tiga penubinan teratur tanpa lubang. Sebuah bidang dapat ditutupi sepenuhnya (ditubin) dengan segitiga sama sisi, dengan persegi, atau dengan segi enam beraturan. Bila dalam kisi-kisi itu rusuk mewakili ikatan molekul dan simpul mewakili molekul, maka transisi fase padat(segitiga)→padat(persegi) terbayangkan. Di situ molekul pada tubin segitiga masih memiliki enam ikatan ke molekul tetangga, sedangkan molekul pada kisi persegi hanya empat.

Peleburan (p/c). Pada suhu yang naik, jaring robek di cukup banyak tempat sehingga blok kaku keseluruhan kehilangan kekakuannya (RETAK kooperatif turun di bawah ambang rigiditas Maxwell z̄ = 2d — jaring menjadi lentur). Terbentuklah misel termodinamis, yaitu daerah yang tak lagi terhubung lewat ikatan dengan daerah tetangga yang lebih besar. Monolit dari sudut molekul menjadi timbunan kerikil dan dapat mengalir secara terbatas. Rongga pertama memungkinkan pergeseran blok molekul mirip puzel geser. Ikatan yang hilang ke lingkungan memungkinkan rotasi internal dalam misel. Rotasi misel (TRINADA-c) mulai bekerja. Pematahan banyak ikatan (IKATAN) melahap kalor laten.

Pendidihan (c/g). Pemanasan lebih lanjut membuat rongga itu sendiri berperkolasi: rongga menyatu menjadi ruang kosong yang menerus, tempat blok-blok terurai dan molekul tunggal terbang pergi. Pendidihan adalah transisi dari "rongga dalam materi" ke "materi dalam rongga".

Sublimasi (p/g). Pada tekanan rendah, daerah antara berupa cairan tidak ada: jaring melepaskan molekul tunggal langsung ke ruang kosong, tanpa terjadi pembentukan blok seperti pada misel. Di sini dapat direnungkan pengaruh apa yang dimiliki aspek geometris (molekul panjang, molekul yang saling menjalin), aspek elektrostatik (penjajaran oleh momen magnet atau momen dipol), dan berbagai jenis ikatan molekul (van der Waals, antaraksi antarsistem mesomeri, ikatan hidrogen) terhadap terjadi atau tidak terjadinya misel termodinamis.

Titik tripel. Tepat di sana dua ambang berimpit: rigiditas Maxwell (padat ↔ cair, hilangnya kekakuan) dan perkolasi rongga (cair ↔ gas) muncul pada (τ, p) yang sama. Ketiga gambaran (blok kaku [monolit], blok yang bergeser [timbunan kerikil], molekul bebas [awan debu]) hidup berdampingan.

Titik kritis. Dengan tekanan yang naik, lompatan cair→gas makin kecil; pangsa rongga berubah makin sinambung. Di titik kritis perbedaannya lenyap sama sekali — "blok rapat dengan sedikit ruang kosong" dan "banyak ruang kosong dengan sedikit molekul" beralih satu ke lainnya secara sinambung. ("Dari sudut pandang zat cair gelas itu setengah kosong, dari sudut pandang gas setengah penuh.")

Rekristalisasi dalam gambaran 2D. Komik kini menampilkan rekristalisasi sebagai kasus 2D yang tersusun ke bawah: di bawah kisi triangular yang rapat (koordinasi z = 6, ideal 2D dengan energi terendah), di atas kisi persegi yang renggang (z = 4). Transisi triangular → persegi memperbesar volume (kerja melawan tekanan) dan menurunkan jumlah ikatan. Terlihat bahwa bilangan koordinasi berganti, tetapi kekakuan sebuah kisi kristal tetap ada.

Ramalan kualitatif transisi energi dan kapasitas kalor. Dari keterhubungan z saja kedua besaran itu dapat ditaksir: (i) energi ikatan per partikel adalah U/N ≈ −(ε/2)·z (tiap ikatan dihitung separuh), jadi makin rendah energinya makin tinggi z; (ii) lompatan c_v pada pergantian kisi mengikuti kaidah model Δc_v = perubahan z_max/4. Untuk kasus 2D triangular (z = 6) ↔ persegi (z = 4) itu menghasilkan selisih energi ΔU ≈ ε per partikel dan tangga c_v Δc_v = (6 − 4)/4 = 0,5 k_B — separuh nilai 3D (1,0 k_B). Rumusan hafalan: makin rapat = makin rendah energinya, tetapi makin mahal dalam transisi (kalor laten lebih besar, tangga c_v lebih besar), karena saat pengembangan lebih banyak ikatan dan lebih banyak derajat kebebasan beku yang dibebaskan.

Apakah struktur kisi menentukan kemetastabilan? Sebagian. Struktur kisi menetapkan urutan energi — bentuk mana yang menjadi keadaan dasar (triangular, μ terendah) dan mana yang tereksitasi (persegi). Namun apakah bentuk renggang itu bertahan metastabil pada suhu rendah, bukan ia sendirian yang memutuskan: dalam model yang benar-benar ideal (tanpa penghalang) zat pada tiap (τ, p) langsung mengambil fase dengan μ terendah — kemetastabilan tak ada di sana. Ia baru timbul lewat energi aktivasi hingga (penghalang penginti), yang menghambat penataan ulang secara kinetik. Jadi geometri kisi memberi tata jenjang termodinamis, penghalang memberi kegigihan kinetik — bahan yang sama yang juga menerangkan fase lewat jenuh (lihat bawah; secara kuantitatif Bab 6.7).

Fase lewat jenuh — energi aktivasi. Gambar tambahan (3.2, baris 2) memperlihatkan sebuah fase gas lewat jenuh di atas suatu cairan: di samping molekul tunggal terdapat kelompok kecil berisi 2–3 partikel yang menempel (batang kelabu) dan terbang sebagai satu kesatuan. Model menerangkan kelewatjenuhan itu lewat energi aktivasi yang nyata bagi ikatan molekul: gugusan kecil (2–3) memiliki nisbah tepi terhadap luas yang tak menguntungkan, ikatan tepinya tak lengkap — hanya sebuah inti kritis yang tumbuh terus secara spontan. Selama tak terbentuk inti superkritis, gas tetap lewat jenuh, meski μ_gas > μ_cair sudah mengizinkan pengembunan. Justru penghalang inilah yang tak ada dalam gambaran kesetimbangan ideal dan di sini disebut sebagai sumbangan mandiri yang "terletak di luar gambaran ideal" (taksiran kuantitatif pada Bab 6.7).

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

3.3 Mengapa laju bunyi dalam zat cair turun bersama T

Bunyi adalah gelombang elastis: v = √(K_eff/ρ). K_eff timbul semata dari ikatan yang utuh (IKATAN) — merekalah pegasnya. Bila T naik, makin banyak ikatan patah secara kooperatif (RETAK); jaring pegas menipis dan K_eff turun lebih curam daripada kerapatan ρ. Pematahan menguasai pembilang, bukan penyebut — karena itu v turun. Ditafsirkan secara molekul: penerusan bunyi lewat ikatan molekul lebih cepat daripada penerusan mekanis lewat tumbukan. Makin banyak retak muncul dalam suatu cairan ketika suhu naik, makin sering bunyi harus diteruskan lewat tumbukan mekanis alih-alih lewat ikatan molekul. Sejalan dengan itu laju bunyi melambat seiring naiknya suhu. (Anomali air — sebuah maksimum bunyi pada ~74 °C — dibahas Bab 6.4.)

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

3.4 Bagaimana gerak Brown timbul — dengan pengertian misel yang dipertajam

Dalam arti statistik-termodinamis, sebuah misel di sini bukan agregat surfaktan, melainkan sebuah blok yang dibingkai RETAK: ikatan yang hilang (bukan-ikatan pada batas blok) membentuk "kulit" di sekeliling kristal kecil yang terkurung (inti yang terikat utuh). Inti ini dapat melakukan gerak rotasi kooperatif internal (TRINADA-c) — blok berputar sebagai satu kesatuan, sementara kulit batasnya melepaskannya dari para tetangga.

Gerak Brown kini timbul begini: lewat pertukaran ikatan molekul (RETAK, GESIT) molekul tunggal dapat meninggalkan kulit sebuah misel dan berpindah ke misel lain. Setiap perpindahan semacam itu memindahkan momentum (KEKEKALAN; karena TERBANG BEBAS tanpa rugi) dan mengubah bentuk serta momentum sudut kedua misel yang terlibat secara melompat. Sebuah partikel uji besar (butir debu) dibasuh oleh banyak penataan ulang misel yang tak berkorelasi semacam itu; jumlah dorongan momentumnya rata-rata nol, tetapi ragamnya hingga → jalan acak. Makin tinggi T, makin sering peristiwa RETAK dan pergantian misel, makin kuat getarannya. Gerak Brown dengan demikian adalah bayangan makroskopis atau akibat dari proses pertukaran misel yang mikroskopis.

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

3.5 Mengapa reaksi di lapisan batas misel menjadi lebih cepat

Catatan awal yang memalukan. Penalaran berikut adalah turunan dari gambaran itu dan sejauh sepengetahuan saya sampai kini belum dibuktikan atau diteliti secara empiris, dan saya harus mengaku bahwa saya tidak bersusah payah menelitinya dengan saksama. Satu-satunya petunjuk kebaruannya hanyalah bahwa sistem AI tidak menyebut kesimpulan ini. Tesis ini adalah contoh bagaimana bayangan bergambar dapat menuntun ke pertanyaan baru.

Dalam pelarut yang daya larutnya sedang, reaktan lebih suka menempel pada kulit batas misel — dalam model justru pada batas retak/blok (RETAK). Dua efek mempercepat reaksi: (i) reduksi dimensi — reaktan tak lagi mencari satu sama lain dalam seluruh volume, melainkan sepanjang antarmuka yang kuasi-dua-dimensi; pencarian 2D secara statistik lebih cepat bertemu daripada pencarian 3D. (ii) Pemekatan — kulit batas mengumpulkan reaktan, kerapatan setempatnya naik. Secara bergambar batas retak bekerja bagai rel pengarah yang menuntun pasangan itu bersua.

Efek dimensi punya inti yang ketat: jalan acak (Pólya). Mengapa pencarian dua dimensi secara statistik lebih baik menemukan daripada pencarian tiga dimensi tidak hanya menyusul dari intuisi, melainkan dari sebuah teorema teori peluang. George Pólya (1887–1985) menunjukkan pada 1921 untuk jalan acak ("random walk") pada sebuah kisi: dalam satu dan dua dimensi si pejalan dengan peluang 1 (jadi pasti) kembali tak terhingga kali ke setiap titik, pada kisi bidang ℤ² sama halnya seperti pada garis ℤ. Mulai tiga dimensi hal itu tak berlaku lagi: peluang kembalinya benar-benar lebih kecil dari 1 (pada kisi kubus sederhana ℤ³ hanya sekitar 0,34). Secara intuitif, dalam gambaran Pólya: "Seorang pemabuk pasti suatu saat menemukan jalan pulang, tetapi ikan yang mabuk tidak." Diterapkan pada dua reaktan yang untuk tumbukan yang berhasil harus bertemu berkali-kali, itu berarti: dalam pelarut yang sangat baik daya larutnya dengan jalan acak tiga dimensi yang bebas, pasangan yang sekali terpisah dengan peluang hingga akan menjauh untuk selamanya — percobaan reaksi berulang gagal, dan hasilnya tetap rendah semata karena alasan geometris-statistis ini. Sebaliknya pada lapisan batas misel yang (kuasi-)dua-dimensi, pertemuan kembali itu pasti (peluang 1): pasangan itu selalu bersua lagi sampai reaksinya berhasil. Jadi lapisan batas tak hanya memekatkan, melainkan sekaligus menyediakan dimensi yang lebih menguntungkan bagi pertemuan kembali. Klaimnya (yang sengaja ditajamkan) berbunyi: hasil dalam pelarut yang sangat baik daya larutnya cenderung lebih buruk karena alasan jalan acak.

Rujukan: G. Pólya (1921), Über eine Aufgabe der Wahrscheinlichkeitsrechnung betreffend die Irrfahrt im Straßennetz, Math. Ann. 84, 149–160 — di sini dikutip menurut M. Eisermann, Höhere Mathematik 3 (Univ. Stuttgart), soal W-2×2. (Naskah asli diakses 5.9.2026).

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

3.6 Jendela cair tiga dimensi

Berbeda dengan gambar pada 3.1–3.5, gambar ini bukan menampilkan adegan molekul, melainkan hasil perhitungan dari Bab 7.9/7.10 sebagai skema pada sebuah sumbu τ (dalam satuan energi ikatan ε). Tujuannya adalah membuat terlihat mengapa kisi 2D tunggal hampir tidak memberi ruang bagi fase cair, sedangkan pada kisi 3D hal itu mungkin tanpa masalah, padahal untuk kedua kasus berlaku kaidah rigiditas yang sama (z̄_c = 2d).

Gambar ① — 2D (tanpa jendela). Pada sumbu τ, τ_m = 1,20 dan τ_triple = 1,16 praktis berimpit (τ_triple di bawah τ_m), τ_c = 1,50 rapat di atasnya: daerah bagi fase cair tertekan menjadi sebuah celah — padat praktis langsung beralih ke gas (dekat sublimasi).

Gambar ② — 3D (jendela bersih). Sumbu yang sama, tetapi teregang: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00. Antara τ_triple dan τ_c terletak medan cair yang mandiri (lebar ke bawah, τ_triple − τ_m = 1,63; sempit ke atas, τ_c − τ_triple = 0,18). Tepi atas yang sempit ditandai sebagai demikian.

Gambar ③ — sebabnya. z̄_c = 2d tumbuh bersama dimensi (4 → 6), sehingga v_solid turun (0,50 → 0,667 dalam volume tereduksi, artinya zat padat menjadi "lebih renggang" relatif terhadap cairan), dan τ_c = z_max/4 berlipat dua (1,50 → 3,00). Titik kritis bergeser jauh di atas titik lebur. Baru karena itu timbul ruang bagi fase cair. Ini murni efek dimensi.

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

4 — Model statistik-mekanik (tak berubah dari versi 09)

Ringkasan

Bab 4 menyusun kerangka statistik-mekanik model itu dalam 2D (kisi segitiga, z_max = 6):

• 4.1 Derajat kebebasan dan peubah: besaran keteraturannya adalah pangsa ikatan x ∈ [0,1] dengan bilangan koordinasi rata-rata z̄ = z_max·x; per blok ≥ 2 bertambah satu derajat kebebasan rotasi, suhu τ dihitung dalam satuan energi ikatan ε.
• 4.2 Ikatan dua keadaan kooperatif: persamaan ikatan swakonsisten x = 1/(1 + G·exp[−(ε + K·x)/τ]) (K = J·z_max, G = z_max) menyusul sebagai syarat kestasioneran energi bebas Bragg–Williams; minimum globalnya memberi cabang fisis, sedangkan apakah cabang itu berjalan sinambung atau melompat (orde pertama) diputuskan Bab 5.5.
• 4.3 Neraca energi dengan suku kooperatif: U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x² (2D) menyatukan translasi/getaran ikatan/rotasi blok, energi konfigurasi ikatan, dan energi medan molekul kooperatif baru −(K/2)x². Kapasitas kalor terkait c_v memperoleh suku tambahan positif −K·x·ẋ, yang mempertajam puncak lebur, tetapi membiarkan plato padat dan cabang gas tak tersentuh.

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

4.1 Derajat kebebasan dan peubah

• N molekul pada suatu geometri koordinasi acuan yang tetap. Kasus 2D kanonik: kisi segitiga, z_max = 6, satu DOF rotasi per blok (KEMAJEMUKAN, IKATAN). Bab ini berlaku eksplisit untuk tafsiran 2D; perampatan 3D ada pada Bab 6.
• Peubah ikatan s_ij ∈ {0, 1} per pasangan tetangga geometris (1 = utuh).
• Bilangan koordinasi z_i = Σ_j s_ij; rata-rata z̄ = 2 N_bond / N — parameter keteraturan. Berlaku z̄ = z_max · x dengan pangsa ikatan x = p_b ∈ [0, 1].
• Satuan energi τ ≡ k_B T (klasik, ½τ per derajat kebebasan kuadratis). Semua energi dalam satuan energi ikatan ε.

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

4.2 Ikatan dua keadaan kooperatif (bertipe Ising, parameter J)

Sebuah ikatan terisolasi (energi −ε bila utuh, 0 bila patah) dalam kesetimbangan Maxwell–Boltzmann memiliki p_b⁰(τ) = 1/(1 + e^(−ε/τ)). RETAK kooperatif (aksioma 6) menggeser kestabilan ikatan efektif sebesar suatu medan molekul yang sebanding dengan pangsa ikatan x. Dalam hampiran medan molekul (Ising) persamaan ikatan swakonsistennya menjadi:

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

Di sini:

• K = J·z_max adalah kekuatan kopling kooperatif (umpan balik RETAK: makin banyak ikatan tetangga → makin tinggi kestabilan ikatan efektif). J > 0 mempertajam transisi; di atas suatu J kritis z̄(τ) menjadi tak sinambung — itulah kalor laten peleburan. Kasus dua keadaan murni adalah kasus batas J = 0 (tanpa medan molekul, peralihan sinambung).
• G > 1 adalah degenerasi keadaan patah: sebuah molekul yang terbebaskan memperoleh ruang fase (kebebasan translasi). G mendorong pematahan ikatan pada τ tinggi dan memastikan bahwa x → 1/(1+G) alih-alih tersangkut di ½ (kasus batas gas). Secara kanonik G = z_max.

Persamaan ini menyusul sebagai syarat kestasioneran suatu energi bebas Bragg–Williams per molekul (n_slots = z_max/2 tempat ikatan):

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

Suku −(K/2)x² adalah kopling kooperatif (feromagnetik), suku terakhir adalah degenerasi pematahan. Minimum global g atas x memberi cabang fisis; di tempat dua minimum sama dalamnya terletak transisi orde pertama. GESIT (aksioma 7) membenarkan bahwa x(τ) pada tiap suhu mengambil nilai setimbangnya.

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

4.3 Neraca energi — dengan suku kooperatif

Suku medan molekul −(K/2)x² adalah energi sejati (tak bergantung τ) dan harus muncul dalam U/N. Neraca yang konsisten dengan demikian berbunyi (2D):

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

• τ·[…]: translasi (1) + getaran ikatan (z̄/4) + rotasi blok (½ n_blk) — TRINADA.
• −(ε/2)·z̄: energi konfigurasi ikatan — IKATAN.
• −(K/2)·x²: energi medan molekul kooperatif — RETAK. Baru dalam U/N.

Kapasitas kalor terkait c_v = d(U/N)/dτ (dalam satuan k_B) memperoleh sumbangan tambahan:

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

di mana suku baru −K·x·ẋ bernilai positif (ẋ < 0) dan bekerja tepat di tempat banyak ikatan patah serentak. Ia mempertajam puncak lebur, tanpa mengubah plato padat (ẋ ≈ 0) atau cabang gas.

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

5 — Fungsi partisi dan kapasitas kalor

Ringkasan

Bab 5 menghitung fungsi partisi dan kapasitas kalor serta menyediakan data model bagi kalibrasi (Bab 7):

• 5.1 Fungsi partisi terfaktorkan: karena pemisahan skala waktu, Z = Z_trans·Z_vib·Z_rot·Z_bond.
• 5.2 Kapasitas kalor per fase: c_v = d(U/N)/dτ; suku positif −(ε/2)(dz̄/dτ) menghasilkan puncak bertipe Schottky sebagai padanan kalor laten.
• 5.3 Kurva c_v(τ) numerik: puncak c_v 5,51 pada τ_m = 1,20, plato padat 2,50, cabang gas → 1; kalor laten L ≈ 1,41 ε/molekul (dengan suku kooperatif) menjadi tiga kali lipat dibandingkan 0,48 ε (tanpa suku itu).
• 5.4 / 5.4a Dua ambang terpisah: perkolasi keterhubungan z̄_c = z_max·p_c ≈ 2,08 (pada τ_b ≈ 1,85) menandai runtuhnya gugusan (dekat pendidihan) dan merupakan sebuah hasil; ia harus dipisahkan bersih dari rigiditas Maxwell generik z̄_c = 2d (kerapatan lebur, hilangnya kekakuan).
• 5.5 Orde transisi peleburan: kriteria ln G > 2 + 4ε/K; 2D (G = 6) sinambung untuk setiap J, 3D (G = 12) orde pertama mulai K > ~8,2 (demo: Δz̄ ≈ 5,3, ΔU ≈ 1,75 ε).
• 5.6 Kajian ukuran hingga: p*(L) memusat secara monoton dari bawah menuju p_c = 0,3473; z̄_c ≈ 2,08 tangguh terhadap ukuran hingga.
• 5.7 Pemindaian c_v 3D penuh: mulai K > K_kritis ≈ 8,25 (J > 0,68) terbuka lompatan sejati Δz̄ dengan kalor laten ΔU; puncak lebur bergeser bersama J ke τ yang lebih tinggi (pergeseran positif +0,5 … +2,2 terhadap puncak 2D yang sinambung).

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

5.1 Fungsi partisi terfaktorkan

Karena pemisahan skala waktu (GESIT), fungsi partisi kanonik terfaktorkan:

Z = Z_trans · Z_vib · Z_rot · Z_bond.

• Z_trans = (1/N!)·(A/λ_T²)^N, λ_T = h/√(2π m τ) (translasi 2D klasik-ideal; KEMAJEMUKAN, KEKEKALAN, TERBANG BEBAS).
• Z_vib = (τ/ħω₀)^(N_bond) (osilator harmonik klasik per ikatan utuh; TRINADA-b).
• Z_rot = ∏_(blok≥2) (2π I_b τ)^(1/2)/ħ (rotator kaku 2D klasik per blok; TRINADA-c).
• Z_bond = Σ_{s} exp[(ε/τ)·Σ s_ij + (K/τ)·suku kooperatif] (jumlah ikatan bertipe Ising; IKATAN, RETAK). Penilaian medan molekul memberi x(τ), z̄(τ) dari 4.2.

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

5.2 Kapasitas kalor per fase (analitik)

Dari U/N (Bab 4.3), c_v = d(U/N)/dτ (dalam satuan k_B):

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

Suku −(ε/2)·(dz̄/dτ) bernilai positif (dz̄/dτ < 0: pemanasan mematahkan ikatan) dan mendominasi di tempat banyak ikatan patah serentak — itulah puncak bertipe Schottky, padanan kalor laten.

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

5.3 Kurva c_v(τ) numerik dengan suku kooperatif — kalor laten dikuantifikasi

Ikatan swakonsisten (minimum global), n_blk(τ) dari Monte Carlo perkolasi ikatan (Union-Find, L = 100, 5 realisasi), c_v = dU/dτ (selisih pusat), U/N termasuk −(K/2)x² (cv_scan_09.py, ε = 1, J = 0,55, G = z_max = 6). Titik simpul terpilih:

Hasil (numerik). Urutan tiga fase tetap tertegaskan dan kini lebih tajam:

c_v(padat) ≈ 2,50 < c_v(cair, maksimum) ≈ 5,51 > c_v(cabang gas) ↓ → 1.

Kalor laten (luas puncak di atas garis dasar). Bila (c_v − garis dasar) diintegralkan atas jendela peleburan τ ∈ [0,70 ; 1,70] (garis dasar = garis lurus melalui c_v(0,70) dan c_v(1,70)), diperoleh

L ≈ 1,41 ε/molekul (dengan suku kooperatif), L ≈ 0,48 ε/molekul (tanpa suku kooperatif, neraca lama), ΔL ≈ +0,92 ε/molekul.

Suku kooperatif mengangkat puncak dengan faktor ≈ 1,7 (3,25 → 5,51) dan melipattigakan kalor laten. Secara fisis: pematahan serentak banyak ikatan tetangga (RETAK) menyimpan bukan hanya energi konfigurasi (−ε/2·z̄), melainkan juga energi medan molekul (−K/2·x²) — keduanya dibebaskan saat peleburan.

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

5.4 Ambang perkolasi keterhubungan sebagai hasil

Pangsa gugusan raksasa P_∞(τ) runtuh pada τ ≈ 1,85; di sana x ≈ 0,347 = p_c(segitiga) = 2·sin(π/18) = 0,3473, jadi z̄_c = z_max·p_c ≈ 2,08 — sebuah hasil, bukan penetapan. Inilah ambang keterhubungan (gugusan menerus) dan ia menandai runtuhnya gugusan dekat pendidihan (τ_b); ia tidak boleh dikacaukan dengan rigiditas Maxwell (z̄_c = 2d, peleburan) — bandingkan pemisahan bersih pada Bab 5.4a.

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

5.4.a Dua ambang berbeda dipisahkan dengan bersih

Penamaan ganda "perkolasi/rigiditas" yang sejauh ini kabur kini diuraikan. Yang dimaksud adalah dua ambang yang berbeda secara fisis dan tak boleh dikacaukan:

Alasan. Peleburan adalah hilangnya kekakuan mekanis suatu blok rapat, bukan runtuhnya keterhubungan. Karena itu besaran yang relevan adalah pencacahan kendala Maxwell (medan rerata / rigiditas generik): jaring berdimensi D dari gaya sentral menjadi isostatik tepat ketika koordinasi rata-rata mencapai z̄ = 2D (2D → 4, 3D → 6). Keterhubungan ikatan (z̄_c = 2,08) terletak jauh di bawahnya; jaring yang pada z̄ = 2,08 baru saja masih terhubung sudah lama lentur (moda lunak), jadi sudah cair. Karena itu keterhubungan menandai runtuhnya gugusan yang lebih kemudian dekat pendidihan, sedangkan rigiditas menetapkan pemadatan/pencairan yang lebih awal (kerapatan lebur). Pemisahan ini memaku titik tripel pada Bab 7.8 tanpa parameter dan secara realistis; ia tidak mengubah ambang keterhubungan 2,08 yang ditegaskan pada 5.4/5.6, melainkan hanya tafsirannya (τ_b = dekat pendidihan, bukan peleburan).

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

5.5 Orde transisi peleburan

Kriteria analitik. Daerah spinodal (dua minimum g) ada untuk τ < K/4; transisinya terletak kira-kira pada τ_t ≈ (ε + K/2)/ln G. Lompatan sejati (transisi orde pertama) menuntut τ_t < K/4, jadi pada limit K → ∞:

ln G > 2 + 4ε/K.

Dari situ menyusul hasil yang menarik tentang degenerasi patah G = z_max:

Penegasan numerik (maks. |Δx| per langkah τ, detect_first_order):

Tafsiran. Dalam model 2D kanonik (G = 6) peleburan pada asasnya sinambung (transisi "lunak"/ peralihan) — untuk tak satu pun J timbul lompatan sejati, karena ln 6 < 2. Baru koordinasi/degenerasi yang lebih tinggi dari kemasan 3D (G = 12) mengizinkan, mulai kooperativitas yang memadai (K > 8,2), sebuah transisi orde pertama dengan kalor laten hingga. Itu selaras dengan kenyataan: peleburan 2D kerap sinambung (bertipe KTHNY), peleburan 3D orde pertama. J = 0 (tanpa medan molekul) selalu merupakan kasus acuan yang sinambung. — Secara fisis J menandakan nisbah energi kooperatif terhadap energi ikatan (K/ε = J·z_max/ε): seberapa kuat sebuah ikatan tetangga yang utuh menstabilkan ikatan sendiri.

Peragaan orde pertama (bermotif 3D, cv_scan_09.py, J = 1,0, G = 12, z_max = 12): z̄ melompat sebesar Δz̄ ≈ 5,30 pada τ ≈ 2,82; kalor laten sebagai lompatan energi ΔU ≈ 1,75 ε/molekul. K yang sama pada kisi G = 6 (2D) tetap sinambung — sifat orde bergantung semata pada G = z_max.

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

5.6 Kajian ukuran hingga untuk ambang perkolasi

Perkolasi ikatan pada kisi segitiga L × L (periodik), pangsa gugusan raksasa P_∞(p) diskalakan atas L (cv_scan_09.py; L = 50/100/200/400 dengan 8/5/3/2 realisasi). Ambang eksak p_c = 2·sin(π/18) = 0,34730 → z̄_c = z_max·p_c = 2,0838.

P_∞(p) (rerata ± sebaran):

Taksiran ambang p*(L) (P_∞ memotong 0,5):

Hasil. p*(L) memusat secara monoton dari bawah menuju p_c = 0,3473; sebaran pada p_c turun dari ±0,103 (L = 50) menjadi ±0,017 (L = 400) — transisinya menajam bersama L, tetapi tidak menggeser ambangnya (nilai perkolasi universal). Dengan itu z̄_c ≈ 2,08 = z_max·p_c tertegaskan sebagai nilai model yang tangguh terhadap ukuran hingga.

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

5.7 Pemindaian c_v(τ) 3D penuh untuk peleburan orde pertama

Satu kali jalan peragaan orde pertama dari Bab 5.5 (J = 1,0, G = 12) kini dikembangkan menjadi sebuah pemindaian c_v(τ) penuh atas kooperativitas J (firstorder_3d_16.py, keterhubungan FCC untuk n_blk, neraca energi 3D U/N = τ·[3/2 + z̄/4 + (3/2)n_blk] − (ε/2)z̄ − (K/2)x²). Kriteria orde analitik dari Bab 5.5 (orde pertama ⇔ K > K_kritis = 4ε/(ln 12 − 2) = 8,25) memisahkan kasus-kasusnya; puncak acuan 2D yang sinambung terletak pada τ_2D = 1,20 (c_v ≈ 5,5).

*kasus sinambung: Δz̄/ΔU adalah peralihan curam (bukan lompatan sejati); τ_transisi di situ adalah puncak c_v.

Hasil. (i) Kasus 3D sinambung J = 0,55 memberi τ_transisi = τ_m ≈ 1,72 — uji silang terhadap kinetika ikatan yang bebas (bondkinetics_3d_15.py: τ_m(FCC) = 1,70) cocok. (ii) Mulai K > K_kritis ≈ 8,25 (J > 0,68) terbuka lompatan sejati Δz̄ dengan kalor laten hingga ΔU (lompatan energi melintasi ketakmalaran) — peleburan menjadi orde pertama, persis seperti dituntut kriteria ln 12 > 2. (iii) Uji silang J = 1,0 mereproduksi nilai peragaan Bab 5.5 Δz̄ ≈ 5,4. (iv) Tinggi lompatan Δz̄ dan letak τ_transisi tumbuh bersama J: transisinya bergeser ke atas — koordinasi yang lebih tinggi dari kemasan 3D dan kooperativitas yang lebih kuat menaikkan suhu lebur. Terhadap puncak 2D yang sinambung (τ_2D = 1,20) dengan demikian selalu muncul pergeseran puncak positif (+0,5 … +2,2) yang bertambah monoton bersama K. Sifat orde pertama peleburan 3D dengan itu tidak lagi terbukti hanya pada satu titik, melainkan sepanjang seluruh rentang kooperativitas.

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

6 — Kemasukakalan matematis tafsiran 3D

Ringkasan

Bab 6 membuat perampatan gambaran ikatan dari 2D ke 3D masuk akal secara murni matematis (tanpa visualisasi HTML tersendiri):

• 6.1 Derajat kebebasan rotasi 2D → 3D: rotasi blok tumbuh dari satu menjadi tiga sumbu, translasi dari dua menjadi tiga arah, z_max dari 6 menjadi 12; nilai dasar c_v naik (c_v,padat 2,5 → 4,5, c_v,gas 1 → 3/2), bentuk kurva c_v(τ) tetap.
• 6.2 Keterhubungan sebagai ciri pembeda: z_max adalah besaran terkuantisasi dan memisahkan kemasan terapat (FCC/HCP, z_max = 12) dari BCC (z_max = 8) secara melompat, tanpa keadaan antara.
• 6.3 Transisi tambahan padat(1)→padat(2): dari kediskretan menyusul polimorfisme FCC↔BCC dengan tangga c_v 1,0 k_B tanpa parameter; dalam gambaran ideal yang ketat (θ_BCC = θ_FCC) ia tidak muncul, baru BCC yang lebih lunak (θ_BCC < θ_FCC di bawah ambang pembuka 0,61) membuka τ_ss < τ_m.
• 6.4 / 6.5 Renungan: anomali bunyi air cair sebagai pergantian tata dekat yang kabur antara dua tipe kisi, dan sebuah alternatif medan pusaran fraktal bagi rotasi misel (murni kualitatif).
• 6.6 / 6.8 / 6.9 θ(kisi) mikroskopis pada tiga tingkat: model pegas Einstein (θ_BCC/θ_FCC ≈ 0,82…0,94), DOS fonon bertipe Debye (≈ 0,87…1,00) dan kisi penuh Lennard-Jones (≈ 1,0; BCC ideal tak stabil terhadap geser, FCC keadaan dasar yang stabil) — semuanya di atas ambang 0,61, jadi polimorfisme nyata terletak di luar gambaran ideal harmonik.
• 6.7 Kemetastabilan (2D): perhitungan orde besaran untuk penataan ulang triangular ↔ persegi (ΔU = +1 ε, ΔA ≈ +15,5 %) beserta penghalang penginti ΔG\* = πλ²/Δμ_A sebagai tambahan kinetik.

Bab ini membuat perampatan 3D masuk akal secara murni matematis. 3D tidak pernah divisualkan dalam versi HTML. Tujuannya semata menunjukkan bagaimana derajat kebebasan rotasi dan keterhubungan berubah dari 2D ke 3D dan transisi baru apa yang karenanya menjadi terbayangkan.

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

6.1 Perubahan derajat kebebasan rotasi 2D → 3D

Satu-satunya butir yang secara kualitatif baru dibanding 2D adalah rotasi blok:

• 2D: sebuah blok kaku punya 1 DOF rotasi (putaran dalam bidang) → U_rot/N = ½ n_blk τ.
• 3D: sebuah blok kaku punya 3 DOF rotasi (tiga sumbu Euler) → U_rot/N = (3/2) n_blk τ.

Sekaligus translasi tumbuh dari 2 menjadi 3 DOF (U_trans/N = (3/2)τ) dan bilangan koordinasi maksimum dari z_max = 6 (segitiga) menjadi z_max = 12 (kemasan bola terapat). Neraca energi dalam 3D dengan demikian berbunyi

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

dan nilai dasar kapasitas kalor bergeser secara konsisten:

Struktur kurvanya (plato padat, maksimum cair oleh −(ε/2)dz̄/dτ, penurunan ke nilai gas) tetap sama; hanya arasnya yang berskala dengan dimensi dan z_max.

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

6.2 Pengaruh berbagai kisi kristal — keterhubungan sebagai ciri pembeda

Dalam 3D setidaknya dua kisi ideal perlu dibedakan, yang terpisah oleh keterhubungannya (bilangan koordinasi z_max):

Keterhubungan z_max adalah ciri pembeda yang menentukan: ia menetapkan baik aras kapasitas kalor benda padat (z_max/4) maupun ambang perkolasi z̄_c. Karena setiap kisi hanya mengizinkan satu bilangan koordinasi tetap, keterhubungan adalah besaran terkuantisasi dan berhingga — tidak ada peralihan sinambung antara z_max = 12 dan z_max = 8, melainkan hanya lompatan diskret dari satu tipe kisi ke tipe lain.

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

6.3 Transisi tambahan padat(1) → padat(2)

Dari keberhinggaan/kediskretan keterhubungan menyusul sebuah transisi yang tak ada dalam model 2D: padat(1) → padat(2) (polimorfisme/alotropi). Dua fase padat dengan z_max berbeda (misalnya FCC ↔ BCC) punya entalpi bebas μ(τ, p) berbeda; pada kombinasi suhu/tekanan μ_{padat(1)} = μ_{padat(2)} zat itu mengganti kisinya, tanpa melebur. Dalam gambaran model, geometri koordinasi berubah melompat (bukan hanya pengisian z̄ seperti pada peleburan) — satu-satunya keadaan di mana geometri acuan yang biasanya dipertahankan itu sendiri berpindah. Yang teramati akan berupa tangga c_v kecil (perubahan z_max/4) ditambah kalor transformasi laten, jauh lebih lemah daripada puncak lebur.

Kriteria orde pertama ln G > 2 + 4ε/K (Bab 5.5) kini memberi alasan konkret mengapa justru tafsiran 3D lebih kaya secara kualitatif: hanya di sana (G = z_max = 12, ln 12 = 2,485 > 2) peleburan dapat menjadi orde pertama — dan hanya di sana ada pergantian kisi diskret padat(1) → padat(2) (FCC ↔ BCC), karena keterhubungan terkuantisasi/berhingga. Kedua gejala itu adalah sebab yang sama dalam dua busana: koordinasi kemasan 3D yang lebih tinggi dan diskret menghasilkan baik lompatan lebur yang tajam maupun polimorfisme. Kisi tunggal 2D (G = 6) tak punya keduanya — konsisten dengan Bab 5.5.

Kuantifikasi (bondkinetics_3d_15.py). Transisi itu kini dihitung secara ramping. Kedua kisi masing-masing memperoleh sebuah entalpi bebas dalam gambaran harmonik klasik (bagian p·v dan translasi yang sama saling menghapus dalam perbandingan):

μ(τ) = h₀ − c_v·τ·ln(τ/θ), h₀ = −(ε/2)·z_max, θ = ħω/k_B (kekakuan fonon dalam satuan ε),

dengan nilai ciri tanpa parameter dari Bab 6.1/6.2:

Dari situ tangga c_v tanpa parameter: |Δc_v| = perubahan z_max/4 = (12−8)/4 = 1,0 k_B (benda padat menjadi 1 k_B lebih "ringan" dalam kapasitas kalor pada transisi FCC→BCC) — jauh lebih kecil daripada puncak lebur (beberapa k_B, Bab 5.3/7.11).

Suhu transformasi τ_ss (μ_FCC = μ_BCC) dan kalor transformasi laten L = τ_ss·(S_BCC − S_FCC):

Hasil (temuan jujur). Dalam gambaran ideal yang ketat (kekakuan fonon sama untuk semua ikatan, θ_BCC = θ_FCC) tidak ada transisi padat(1)→padat(2): FCC menang baik dalam entalpi (lebih dalam, −6 lawan −4 ε) maupun — karena jumlah ikatan lebih tinggi (c_v = 4,5 > 3,5) — dalam entropi getaran; μ_FCC tetap di bawah μ_BCC pada setiap suhu. Baru sebuah kekakuan fonon yang bergantung kisi θ_BCC < θ_FCC (BCC yang lebih terbuka bergetar lebih lunak, jadi punya lebih banyak entropi) mengangkat μ_BCC pada suhu tinggi ke bawah μ_FCC — barulah transformasi itu terbuka, dan letaknya di bawah puncak lebur τ_m, sehingga timbul urutan padat(1)→padat(2)→cair→gas (τ_ss < τ_m < τ_c ≤ τ_b(3D), Bab 7.11). Kekakuan BCC yang lunak adalah perluasan ringan di luar gambaran ideal yang ketat (nyatanya dikenal sebagai polimorfisme γ/δ besi: BCC pada suhu tinggi, FCC di antaranya); sebaliknya tangga c_v 1,0 k_B dan kalor transformasi positif yang kecil bersifat internal model dan tanpa parameter.

Seberapa kecil sebenarnya θ_BCC/θ_FCC? (theta_springs_16.py). Untuk membebaskan transformasi itu dari parameter θ yang bebas, kekakuan fonon kini ditaksir secara mikroskopis (Bab 6.6): sebuah model pegas Einstein dengan pegas pasangan sentral memberi θ_BCC/θ_FCC ≈ 0,82 (hanya tetangga terdekat) sampai ≈ 0,94 (dengan kulit kedua BCC). Sebaliknya ambang pembuka yang diperlukan dari perbandingan μ terletak pada θ_BCC/θ_FCC < 0,61 (di bawah itu τ_ss ada dan terletak di bawah τ_m). Jadi nisbah kekakuan geometris itu terlalu besar: dalam gambaran ideal yang ketat FCC tetap stabil pada setiap suhu, tidak ada transisi padat(1)→padat(2) yang digerakkan murni secara geometris. Polimorfisme nyata (besi γ↔δ) memerlukan sumbangan tambahan (elektronik/magnetik/anharmonik) di luar gambaran harmonik ideal — persis "perluasan ringan" yang diandaikan di atas. Itu menegaskan temuan yang jujur dan mengubah nisbah θ dari andaian bebas menjadi batas model yang tertaksir dan bernama.

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

6.4 Renungan: anomali bunyi air cair

Zat cair biasa memperlihatkan v(T) yang turun monoton (Bab 3.3). Air punya maksimum bunyi pada ≈ 74 °C: v mula-mula naik bersama T sebelum turun. Dalam model hal itu dapat ditafsirkan sebagai persaingan dua tipe kisi lokal (Bab 6.2): sebuah tata dekat mirip es yang terbuka dan berkoordinasi rendah (z_max kecil, banyak "lubang") dan sebuah tata dekat yang rapat dan berkoordinasi lebih tinggi. Pada T rendah struktur terbuka mendominasi; pemanasan mula-mula mematahkan tata dekat terbuka itu demi yang lebih rapat — keterhubungan efektif dan dengan itu K_eff naik (v naik). Baru di atas maksimum itu pematahan ikatan RETAK yang biasa mengungguli, K_eff turun, v menyusut. Anomali itu dengan demikian merupakan ungkapan sebuah pergantian tata dekat bertipe padat(1)↔padat(2) yang kabur dan tak lengkap di dalam fase cair — hanya terbayangkan dalam model 3D dengan dua tipe kisi, tidak dalam kisi tunggal 2D.

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

6.5 Alternatif: medan pusaran fraktal alih-alih rotasi misel

Gerak Brown lokal (Bab 3.4) dijelaskan lewat rotasi dan pertukaran misel. Sebagai alternatif ia dapat dibayangkan sebagai efek medan pusaran fraktal: batas-batas blok yang terus-menerus ditata ulang oleh RETAK menghasilkan sebuah medan arus momentum lokal yang — mirip aliran turbulen — memperlihatkan pusaran serupa-diri pada semua skala panjang. Sebuah partikel uji lalu digerakkan bukan oleh tumbukan misel yang diskret, melainkan oleh sebuah kontinum pusaran yang invarian skala. Kedua gambaran (misel diskret ↔ pusaran fraktal) adalah kasus batas yang saling melengkapi dari kinetika RETAK yang sama; mana yang berlaku bergantung pada nisbah ukuran partikel uji terhadap ukuran misel rata-rata. Pembahasan awal ini murni kualitatif dan tidak dikerjakan secara matematis di sini.

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

6.6 Kekakuan fonon mikroskopis θ(kisi): sebuah model pegas Einstein

Perbandingan padat(1)→padat(2) pada Bab 6.3 bergantung pada kekakuan fonon relatif θ_BCC/θ_FCC. Alih-alih memilihnya bebas, di sini ia ditaksir dari koordinasi dan geometri (theta_springs_16.py).

Model. Sebuah atom pusat duduk dalam cekungan harmonik dari z tetangganya; setiap ikatan adalah pegas pasangan radial κ. Sebuah simpangan u memproyeksi secara isotrop dengan ⟨cos²⟩ = 1/d = 1/3 (d = 3) pada arah-arah ikatan, sehingga tetapan pegas efektif adalah k_eff = (z/3)·κ. Frekuensi Einstein adalah ω = √(k_eff/m), jadi θ = ħω/k_B ∝ √k_eff. Pada massa sama dan pegas pasangan sama menyusul

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

Penghalusan. (B) Kulit kedua BCC (6 atom pada r₂/r₁ = 2/√3 = 1,155) dengan pegas berskala jarak κ(r) ∝ (r₁/r₂)ⁿ (n = 6…8, bertipe Lennard-Jones) menaikkan nisbah itu menjadi 0,91…0,94 (BCC menjadi efektif lebih kaku oleh kulit kedua yang dekat). (C) Ikatan yang lebih pendek pada BCC pada kerapatan sama (r_nn(BCC)/r_nn(FCC) = 0,972, −2,8 %) sedikit mengakukan pegas dan memberi 0,89…0,91. Seluruhnya sebuah jendela geometris θ_BCC/θ_FCC ≈ 0,82…0,94.

Ambang pembuka. Perbandingan μ (Bab 6.3) membuka transformasi itu (τ_ss ada dan terletak di bawah τ_m = 1,70) baru untuk θ_BCC/θ_FCC < 0,61. Jendela geometris terletak jauh di atasnya.

Temuan. Model pegas murni tidak membuat BCC cukup lunak: dalam gambaran ideal yang ketat FCC tetap stabil, sebuah transisi padat(1)→padat(2) yang digerakkan secara geometris tidak timbul. Dengan itu θ_BCC/θ_FCC berubah dari andaian bebas menjadi besaran tertaksir, yang nilai alaminya (≈ 0,8–0,9) menahan transformasi itu tetap tertutup — konsisten dengan temuan jujur pada Bab 6.3. Polimorfisme nyata dengan demikian terbukti terletak di luar gambaran harmonik ideal (bagian elektronik/magnetik, anharmonisitas), bukan pada geometri kisi belaka.

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

6.7 Taksiran matematis kemetastabilan (2D triangular ↔ persegi)

Untuk pembahasan kualitatif dari Bab 3.2, berikut perhitungan orde besaran yang bersangkutan dalam gambaran 2D (cakram keras, garis tengah D, jarak kontak a = D).

(1) Energi. Per partikel U/N = −(ε/2)·z (setiap ikatan setengah): triangular z = 6 → −3 ε, persegi z = 4 → −2 ε. Penataan ulang triangular → persegi jadi berbiaya ΔU = +1 ε per partikel (endotermik; bentuk renggang lebih kaya energi).

(2) Volume (luas) dan kerja. Luas per partikel: triangular A_tri = (√3/2)·D² ≈ 0,866 D², persegi A_quad = D². Pelebarannya adalah ΔA = A_quad − A_tri ≈ 0,134 D² (+15,5 %) — pembesaran volume yang diperlihatkan dalam komik. Terhadap tekanan p (tereduksi) untuk itu harus dikeluarkan kerja W = p·ΔA; entalpi bebas naik sebesar Δμ = ΔU + p·ΔA − τ·ΔS. Suku entropi ΔS > 0 (bentuk renggang punya moda lebih lunak, tangga c_v lebih besar Δc_v = (6−4)/4 = 0,5 k_B) baru menurunkan μ_quad di bawah μ_tri pada τ tinggi — karena itu triangular adalah keadaan dasar suhu rendah dan persegi bentuk suhu tinggi/metastabil.

(3) Penghalang penginti (2D). Kemetastabilan mensyaratkan sebuah penghalang. Dalam 2D sebuah inti fase baru (jari-jari r) tumbuh dengan perolehan luas −Δμ_A·π r² (Δμ_A = gaya dorong per luas) melawan sumbangan tepi +λ·2π r (λ = energi tepi per panjang): ΔG(r) = −Δμ_A·π r² + 2π λ r. Maksimum pada r\* = λ / Δμ_A dengan penghalang ΔG\* = π λ² / Δμ_A. Inti kritis r\ menyusut seiring bertambahnya gaya dorong Δμ_A; selama skala termal τ = k_BT belum mencapai penghalang (**ΔG\ ≳ beberapa k_BT), fase renggang (atau, pada pendidihan, yang lewat jenuh) tetap beku secara kinetik. Bagi gugusan kecil (2–3 partikel, seperti pada strip kelewatjenuhan) nisbah tepi terhadap luas paling tak menguntungkan — mereka terletak di kiri r\* dan rata-rata larut kembali. Dengan itu energi aktivasi yang diperlihatkan dalam komik adalah besaran ΔG\* yang sama, yang menopang baik bentuk kristal metastabil maupun fase gas lewat jenuh; ia adalah tambahan kinetik** bagi model kesetimbangan ideal (tanpa penghalang).

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

6.8 Penghalusan θ(kisi) bertipe Debye

Model pegas Einstein (Bab 6.6) mengganti seluruh sebaran fonon dengan satu frekuensi. Pendekatan yang ada di sini menghaluskan θ(kisi) melalui sebuah kerapatan keadaan fonon yang ramping per kisi (theta_debye_17.py): pegas pasangan sentral yang sama κ(r) = r^(−n) pada kulit tetangga terdekat dan kedua membentuk matriks dinamika D(k) = Σ_j κ_j·(R̂_j⊗R̂_j)·(1 − cos(k·R_j)), yang didiagonalkan pada sebuah kisi k dari zona Brillouin primitif (Monkhorst-Pack, 12³ titik) (penyelesai nilai eigen 3×3 analitik, hanya pustaka standar). Dari ketiga cabang fonon menyusul dua ukuran kekakuan bertipe Debye sebagai momen frekuensi:

• Momen ke-2 θ(2) ∝ √⟨ω²⟩ — kasus batas Einstein (rerata BZ dari 1 − cos menghasilkan 1, jadi ⟨ω²⟩ = z·κ/3); kontrol nn saja, geometri murni: √⟨ω²⟩_BCC / √⟨ω²⟩_FCC = 0,816 = √(8/12) — persis nilai Einstein dari Bab 6.6.
• Momen log θ(0) ∝ exp⟨ln ω⟩ — besaran yang sungguh bertipe Debye, yang memberi bobot pada cabang akustik lunak yang diabaikan nilai Einstein.

Semua tetapan berskala (ħ/k_B, skala κ, prafaktor Debye) saling menghapus dalam nisbah θ_BCC/θ_FCC; hanya bentuk DOS yang tersisa. Dengan kulit tetangga kedua (yang justru baru menstabilkan BCC secara mekanis) dan hukum pegas n = 0…8 diperoleh

θ_BCC/θ_FCC ≈ 0,87…1,00 (pita galat atas definisi momen dan hukum pegas).

Itu cenderung lebih tinggi daripada nilai Einstein murni (0,82…0,94), karena ikatan BCC yang lebih pendek pada kerapatan sama mengakukan pegasnya. Ambang pembuka 0,61 yang diperlukan dari perbandingan μ (Bab 6.3) tetap terletak jauh di bawahnya: DOS fonon pun tidak membuat BCC cukup lunak. Temuan dari Bab 6.6 dengan demikian ditegaskan dan dipertajam oleh hampiran yang lebih halus — dalam gambaran ideal yang ketat FCC tetap stabil pada setiap suhu; sebuah transisi padat(1)→padat(2) yang digerakkan murni secara geometris tidak timbul, polimorfisme nyata (besi, Bab 7.13) memerlukan sumbangan di luar gambaran harmonik ideal.

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

6.9 θ(kisi) dari sebuah kisi penuh Lennard-Jones

Bab 6.8 memakai pegas sentral terpotong κ(r) = r^(−n) pada nn + 2nn. Versi 17 menguji pertanyaan yang sama dengan sebuah tetapan gaya sejati dari potensial pasangan atas kisi penuh (theta_ljlattice_18.py): setiap ikatan menjadi pegas dari potensial Lennard-Jones φ(r) = 4ε[(σ/r)¹² − (σ/r)⁶] dengan bagian longitudinal κ_L = φ″(r) dan bagian transversal κ_T = φ′(r)/r. Matriks dinamika

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

menjumlah atas semua kulit tetangga sampai suatu cutoff (2,0…3,0 nn); setiap kisi sebelumnya direlaksasi ke kesetimbangan LJ miliknya sendiri (tetapan kisi bebas, kerapatan sama tidak dipaksakan). Itulah uji jujur apakah hampiran ramping nn+2nn (Bab 6.8) mengenai urutan kekakuan dengan benar.

Hasil (deterministik, hanya pustaka standar):

Dua temuan, keduanya mempertajam Bab 6.6/6.8:

1. FCC adalah keadaan dasar LJ yang stabil (E ≈ −8,4 ε/atom pada cutoff 3,0 nn, lawan −8,0 ε untuk BCC) — keunggulan kemasan rapat yang dikenal dalam potensial pasangan.
2. BCC LJ ideal tak stabil terhadap geser: sekitar 3…3,5 % moda fonon bersifat imajiner (ω² < 0, ketakstabilan Born). Gambaran ideal potensial pasangan murni jadi sama sekali tidak menopang BCC yang stabil secara mekanis; atas moda riil yang tersisa θ_BCC/θ_FCC ≈ 1,0 — ambang pembuka 0,61 yang diperlukan dari perbandingan μ justru semakin tidak dilewati ke bawah.

Dengan itu pernyataan dari Bab 6.6/6.8 tangguh terhadap penghalusan dari pegas terpotong ke kisi penuh dan dari model longitudinal murni ke matriks dinamika penuh: dalam gambaran ideal klasik tidak timbul transisi padat(1)→padat(2) yang murni geometris — polimorfisme nyata (besi Bab 7.13, titanium Bab 7.15) memerlukan sumbangan di luar potensial pasangan (elektronik/magnetik, anharmonik). Sebuah perhitungan fonon ab initio sejati sengaja dibiarkan di luar (ia akan meninggalkan sifat gambaran ideal yang ramping).

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

7 — Kalibrasi dan persamaan keadaan

Ringkasan. Bab 7 mengalibrasi model lewat satu-satunya parameter energi ε pada zat nyata dan mengembangkan cabang cair/gas menjadi persamaan keadaan:

• 7.1 / 7.2 Metode dan tabel kalibrasi: pengaitan suhu linear T = (ε/k_B)·τ dengan satu titik jangkar di titik lebur; dari situ nisbah universal yang dihaluskan R = τ_b/τ_m ≈ 1,54.
• 7.3 / 7.4 Persamaan keadaan van der Waals berinti keras: persamaan tereduksi p_r = 8τ_r/(3v_r − 1) − 3/v_r² dengan titik kritis (1, 1, 1) dan koeksistensi Maxwell; a ∝ ε, b ∝ D². Titik didih terletak secara sistematis pada T_b/T_c ≈ 0,48 (jarak Guldberg); sebuah perhitungan pembanding Redlich–Kwong menaikkannya ke ≈ 0,60 — sebagai anggaran galat yang disebutkan, tanpa mengganti persamaan keadaan.
• 7.5 Permukaan (τ, v) bersama: kubah koeksistensi (pendidihan), garis lebur, dan dinding keras menyatu dengan sebuah titik tripel (τ_r,triple ≈ 0,47 untuk v_solid = 0,40).
• 7.6 Kalor laten: L ≈ 1,41 ε menghasilkan ΔS_fus ≈ 1,175 R (dekat dengan kaidah Richards); logam ditepati sangat baik (tembaga ×1,00), kristal ionik/berikatan hidrogen/molekuler berada di luar gambaran ideal.
• 7.7 Penyeragaman ε ↔ ε_lg: k_B T_c = z_max·ε/4 mengopel kedua skala (ε_lg = (27/32)·z_max·ε); ramalan 3D τ_c/τ_m = 2,50 menepati zat molekuler sederhana, 2D tetap kasus khusus.
• 7.8 / 7.9 Titik tripel tanpa parameter: v_solid menyusul dari ambang rigiditas z̄_c = 2d, bukan dari keterhubungan; di 3D terbuka jendela cair yang bersih (τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00), di 2D ia tetap tertutup (peran khusus 2D yang terbukti rangkap tiga).
• 7.10 / 7.11 Pemindaian 3D sejati: perkolasi ikatan pada FCC/BCC memberi z̄_c(3D) ≈ 1,5–1,6; pemindaian τ kinetika ikatan menunjukkan τ_m ≈ τ_rig (peleburan = rigiditas, rendah) dan τ_b(3D) ≈ τ_c (pendidihan = runtuhnya keterhubungan di titik kritis).
• 7.12 / 7.13 / 7.15 garis padat(1)→padat(2): dari titik τ_ss lewat tekanan menjadi sebuah garis (dp/dτ = ΔS/Δv > 0, Δv = +8,9 %), sebagai garis batas keempat dalam permukaan (τ, v) dan sebagai garis (τ, p) yang dikalibrasi pada besi (γ→δ), ditegaskan pada titanium (α→β) dan dengan argon sebagai kontrol negatif.
• 7.14 Cabang Redlich–Kwong: sebagai varian gambaran ideal yang mandiri dengan kalibrasi sendiri di samping bawaan van der Waals (T_b/T_c ≈ 0,60, galat −1 %).

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

7.1 Metode kalibrasi (satu parameter energi ε) — disesuaikan

Linear: T = (ε/k_B)·τ, satu titik jangkar per zat. Penjangkaran pada titik lebur: ε/k_B = T_m,eksp/τ_m, ramalan T_b,pred = T_m,eksp·(τ_b/τ_m). Dengan puncak lebur yang dihaluskan (suku kooperatif), maksimum c_v kini terletak pada τ_m = 1,20 (bukan 1,15), sedangkan perkolasi keterhubungan (dekat pendidihan) tetap pada τ_b = 1,85. Dari situ nisbah universal yang dihaluskan:

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

Tabel kalibrasi dari versi 08 (Bab 7.2) berskala linear dengan R: semua T_b,pred turun dengan faktor 1,54/1,61 ≈ 0,957 (≈ −4,3 %). Itu sedikit memperbaiki taksiran berlebih pada gas mulia (argon +54 % → +48 %) dan sedikit memperburuk ketepatan NaCl/besi (NaCl −0,6 % → −4,8 %). Pernyataan pokoknya tetap: R cocok untuk zat yang terkemas rapat secara isotrop dan berikatan tak berarah; gas mulia/rantai berada di luar gambaran ideal.

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

7.2 Tabel kalibrasi

Diambil tanpa perubahan dari modell_08.md Bab 7.2 (dengan catatan penskalaan R linear dari 7.1). Dapat direproduksi lewat cv_scan_08.py (tabel) atau cv_scan_09.py (τ_m, τ_b, R).

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

7.3 Sumbu kerapatan/volume: persamaan keadaan van der Waals berinti keras

Tekanan diangkat dari "parameter" menjadi peubah sejati. INTI KERAS memberi volume sendiri b, tarikan ikatan (ε·z̄) memberi tekanan kohesi a. Dalam satuan tereduksi (titik kritis pada (v_r, τ_r, p_r) = (1, 1, 1)):

p_r = 8 τ_r / (3 v_r − 1) − 3 / v_r², v_r > 1/3 (dinding keras).

Koeksistensi lewat konstruksi Maxwell (cv_scan_09.py, potongan luas sama):

Hasil. Di bawah τ_c terbukalah lompatan kerapatan c/g (v_gas/v_liq tumbuh dari 2,5 ke 11,5; makin dingin makin besar). Itulah garis didih — kini terpisah bersih dari puncak c_v peleburan (Bab 5.3), yang hidup pada kisi berkerapatan tetap. Dengan itu butir 3 yang terbuka pada versi 08 tertutup: peleburan (tata rigiditas/ikatan pada kerapatan tetap) dan pendidihan (lompatan kerapatan/volume) adalah sumbu-sumbu berbeda dari model yang sama.

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

7.4 Kalibrasi c/g kuantitatif: a → ε, b → D²

Pemetaan besaran tereduksi ke zat nyata. Tetapan van der Waals bergantung pada titik kritis:

b = k_B·T_c / (8 p_c) (volume sendiri, INTI KERAS), a = 27 k_B²·T_c² / (64 p_c) (kohesi, IKATAN).

Dengan itu a ∝ ε (kohesi berasal dari energi ikatan) dan b ∝ D² (volume sendiri 2D dari garis tengah inti keras). Skala energi alami cabang c/g adalah

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

Zat jangkar argon (T_c = 150,87 K, p_c = 48,63 bar): b = 5,35·10⁻²⁹ m³/molekul, a = 3,76·10⁻⁴⁹ J·m³/molekul², ε_lg/k_B = 509 K dan garis tengah molekul D ≈ b^{1/3} = 3,8·10⁻¹⁰ m — terhadap garis tengah van der Waals argon (~3,4·10⁻¹⁰ m) tepat sampai 11 %.

Ramalan suhu didih. Pada tekanan luar tetap, titik didih ada di tempat tekanan jenuh tereduksi memenuhi p_sat(τ_r)·p_c = p_atm; dari situ T_b = T_c·τ_r,b (kurva kejenuhan Maxwell dari lg_calib_10.py):

*CO₂ menyublim pada 1 atm (T = 194,7 K = titik sublimasi).

Hasil. Untuk zat molekuler (p_c ≈ 30–70 bar) model memberi seluruhnya T_b/T_c ≈ 0,48–0,51, secara sistematis ~17–24 % di bawah nilai ukur. Itu bukan galat model dalam arti sempit, melainkan jarak yang dikenal antara persamaan van der Waals dan kaidah Guldberg empiris (T_b/T_c ≈ 0,6): isoterm van der Waals menaksir terlalu rendah tekanan jenuh pada suhu rendah. Air (dan logam, yang T_c/p_c-nya memang hanya ditaksir kasar) pada p_r,atm < 0,005 berada di luar daerah van der Waals yang bermakna (τ_r,b < 0,40) — di sana kohesinya tidak lagi bertipe van der Waals (ikatan hidrogen, ikatan logam). Orde besaran dan urutan titik didih ditepati dengan benar.

Penempatan Guldberg sebagai tabel galat. Jarak itu di sini hanya disebutkan, bukan diperbaiki lewat pergantian persamaan keadaan (gambaran ideal tetap bertipe van der Waals). Kaidah Guldberg empiris (Guldberg 1890) menyatakan bahwa titik didih terletak sekitar T_b/T_c ≈ 0,6 (dengan sedikit ketergantungan zat ~0,58–0,64). Model van der Waals tereduksi justru memberi nilai tetap T_b/T_c = τ_r,b ≈ 0,48–0,51:

Tafsiran. Nilai rerata model 0,49 menepati tetapan Guldberg 0,60 sampai sekitar −18 %; penyimpangannya sistematis dan sepihak (selalu terlalu rendah), jadi merupakan sifat isoterm van der Waals (tekanan jenuh terlalu kecil pada τ rendah), bukan sifat gambaran molekulernya. Sebuah persamaan keadaan kohesi yang lebih baik (mis. Redlich–Kwong) akan menaikkan τ_r,b — tetapi di sini sengaja tetap di luar gambaran ideal.

Anggaran galat dengan sebuah perhitungan pembanding Redlich–Kwong (guldberg_rk_15.py). Untuk mengangkakan (bukan menyembuhkan) jarak van der Waals terhadap kaidah Guldberg, uji titik didih yang sama diulang dengan kohesi bertipe Redlich–Kwong yang bergantung suhu a/√T. Dalam satuan tereduksi (titik kritis (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,

dengan konstruksi luas sama Maxwell dan tekanan luar tereduksi p_r,atm yang sama seperti di atas. Hasil (satu angka per zat jangkar, T_b/T_c = τ_r,b):

Hasil. Kohesi RK yang bergantung suhu menaikkan T_b/T_c dari vdW ≈ 0,48 ke ≈ 0,60 — praktis tetapan Guldberg dan nilai-nilai nyata (≈ 0,58–0,63); galat rerata terhadap nilai ukur turun dari −21 % (vdW) ke −1 % (RK). Dengan itu sebab jarak Guldberg terangkakan dengan jelas: justru tiadanya ketergantungan suhu pada kohesi van der Waals (tekanan jenuh terlalu kecil pada τ rendah), bukan gambaran molekulernya. Perhitungan pembanding ini sengaja berada di luar gambaran ideal; modelnya tetap bertipe van der Waals (tanpa pergantian persamaan keadaan) — anggaran galat hanya disebutkan, tidak dipasang.

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

7.5 Permukaan (τ, v) bersama dengan titik tripel

Sumbu lebur dan sumbu didih sejauh ini hidup terpisah: peleburan (tata rigiditas/ikatan) pada kerapatan tetap, pendidihan (van der Waals) lewat lompatan kerapatan. Dalam satu bidang (τ_r, v_r) bersama keduanya menyatu:

• Kubah koeksistensi (pendidihan): v_liq(τ_r) (cabang kiri) dan v_gas(τ_r) (cabang kanan) dari konstruksi Maxwell, yang bertemu di titik kritis (τ_r, v_r) = (1, 1).
• Benda padat (peleburan): pita sempit tepat di atas dinding keras v_r > 1/3, pada kerapatan tinggi dan hampir tetap (volume tereduksi v_solid ≈ 0,40).
• Garis lebur: hampir tegak pada v_solid (kerapatan tetap), yakni suhu rigiditas Maxwell (z̄_c = 2d; Bab 7.8 menurunkan v_solid darinya tanpa parameter — penyandingan lama versi 11 "z̄_c ≈ 2,08" dengan demikian dikoreksi).

Titik tripel adalah titik tempat batas koeksistensi cair v_liq(τ_r) mencapai kerapatan kemasan padat v_solid — di bawahnya tak ada lagi cabang cair (sublimasi langsung):

Hasil. Untuk v_solid = 0,40 titik tripel terletak pada τ_r,triple ≈ 0,47, p_sat ≈ 0,018, dengan volume uap v_gas ≈ 68 (uap sublimasi yang sangat encer). Nisbah τ_triple/τ_c ≈ 0,47 terletak dalam jendela nyata (air 0,42; argon 0,56; CO₂ 0,71). Makin rapat benda padatnya (v_solid makin kecil), makin rendah titik tripel tereduksi dan makin encer uap yang berkoeksistensi — benar secara kualitatif. Dengan itu kedua transisi dan titik tripel terjangkar dalam satu permukaan (τ, v); titik kritisnya duduk di (1, 1, 1). Versi 15 menambahkan pada Bab 7.12 sebuah garis batas keempat — transformasi padat(1)/padat(2) (FCC↔BCC) sebagai dua pita kerapatan padat yang hampir tegak (v ≈ 0,667 dan 0,726), yang saling berganti pada τ_ss.

Tahap pengembangan sebelumnya: v_solid = 0,40 di sini masih ditetapkan. Bab 7.8 menurunkan v_solid tanpa parameter dari ambang perkolasi z̄_c dan menjadikan titik tripel sebuah ramalan sejati.

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

7.6 Kalor laten terhadap entalpi peleburan

Kalor laten model L ≈ 1,41 ε (Bab 5.3) diterjemahkan ke satuan nyata lewat kalibrasi titik lebur ε/k_B = T_m/τ_m (τ_m = 1,20):

ΔH_m,pred = L·R·T_m/τ_m = 1,175·R·T_m, seturut itu ΔS_fus,pred = L·R/τ_m = 1,175·R.

Itu — sungguh mengesankan — hampir tepat kaidah Richards (ΔS_fus ≈ R untuk benda padat berstruktur sederhana). Perbandingan dengan entalpi peleburan terukur (lg_calib_10.py):

Hasil. Model menepati orde besaran (kJ/mol) seluruhnya dan menepati kaidah Richards untuk logam sederhana dengan sangat baik (tembaga ×1,00, besi ×1,28); gas-gas yang dicairkan terletak sedikit di bawahnya (metana ×0,94, nitrogen ×0,87, argon ×0,69). Kristal ionik, berikatan hidrogen, dan molekuler (NaCl ×0,37, air ×0,44, benzena ×0,28, CO₂ ×0,23) memiliki entropi peleburan nyata yang jelas lebih tinggi (ΔS_exp/R = 2,6–5,0), karena saat melebur terbebaskan derajat keteraturan tambahan (pemasangan ion, orientasi molekul poliatom) yang tak tertangkap oleh model berparameter keteraturan tunggal. Batasnya dengan itu disebut jelas: gambaran ideal berlaku untuk zat yang berikatan isotrop dan tak berarah (tanpa pencocokan).

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

7.7 Penyeragaman skala energi: ε ↔ ε_lg

Sejauh ini cabang lebur memikul energi ikatan ε (τ_m = 1,20), sedangkan cabang c/g memikul skala kohesi tersendiri ε_lg/k_B = 3,375·T_c (Bab 7.4). Keduanya dapat diturunkan dari satu parameter model. Kohesi van der Waals a dalam model timbul dari jumlah energi ikatan per molekul (IKATAN): setiap molekul memikul z_max/2 ikatan berenergi ε. Sebuah gas kisi Bragg–Williams medan rerata dengan tarikan tetangga ε memberi hubungan klasik untuk suhu kritis:

k_B·T_c = z_max·ε / 4, jadi τ_c ≡ k_B T_c/ε = z_max/4.

Bila ini dimasukkan ke ε_lg = a/b = (27/8)·k_B·T_c (Bab 7.4), menyusul kopling kedua skala:

ε_lg = (27/8)·(z_max/4)·ε = (27/32)·z_max·ε, tetapan c = 27/32 = 0,844.

Dengan itu puncak lebur, perkolasi, dan titik kritis terletak pada satu skala ε (unify_scale_12.py):

Ramalan τ_c/τ_m = z_max/(4·τ_m): 2D → 1,25, 3D → 2,50. Terhadap zat nyata:

Hasil. Ramalan 3D τ_c/τ_m = 2,50 menepati dengan baik angka nyata zat molekuler sederhana (T_c/T_m ≈ 2,0–2,1); logam terletak lebih tinggi (T_c ditaksir kasar, jenis ikatan lain). Ramalan 2D (1,25) terlalu rendah — dan τ_c(2D) = 1,50 bahkan terletak di bawah τ_b = 1,85 (titik kritisnya akan berada di bawah ambang perkolasi). Itu bukan salah hitung, melainkan peran khusus 2D yang telah dikenal: di sana skala kohesi dan skala ikatan tidak tersusun bersih, jendela cair yang tangguh tak ada (bdk. peleburan 2D yang sinambung, Bab 5.5). Baru dalam 3D (τ_m = 1,20 < τ_c = 3,00) kedua skala terpisah secara bermakna fisis.

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

7.8 Titik tripel tanpa parameter dari z̄_c

Bab 7.5 masih menetapkan v_solid = 0,40. Tahap pengembangan sebelumnya menurunkan kerapatan kemasan padat dari ambang perkolasi. Bilangan koordinasi suatu kisi pada volume tereduksi v_r mengikuti relasi kerapatan–koordinasi (dinding keras v_r = 1/3 ↔ kemasan terapat, z̄ = z_max):

z̄(v_r) = z_max·(1/3)/v_r = z_max/(3 v_r), sebaliknya v_solid = z_max/(3·z̄_c).

Benda padat kehilangan kekakuannya begitu z̄ turun di bawah ambang z̄_c; kerapatan yang bersangkutan adalah v_solid. Titik tripel lalu merupakan titik tempat batas koeksistensi cair v_liq(τ_r) mencapai kerapatan itu (unify_scale_12.py, 2D, z_max = 6):

Versi yang diperbaiki. Kerapatan lebur v_solid ditetapkan oleh rigiditas Maxwell generik z̄_c = 2d (2D → 4, 3D → 6), bukan oleh keterhubungan ikatan (alasannya pada Bab 5.4a: peleburan adalah hilangnya kekakuan, bukan keterhubungan). Dengan itu titik tripel yang dapat diandalkan berbunyi (rigidity_3d_13.py, 2D):

Hasil 2D. Bahkan dengan ambang rigiditas (yang secara fisis benar), kisi tunggal 2D tidak memberi jendela cair yang bersih: v_solid = 0,50 memang menghasilkan angka nyata τ_triple/τ_c ≈ 0,77, tetapi dalam satuan ε mutlak τ_triple = 0,77·1,50 = 1,16 dan dengan demikian terletak di bawah τ_m = 1,20. Ambang keterhubungan (z̄_c = 2,08) bahkan mendorong v_solid ≈ 0,96 ≈ v_c keluar dari kubah koeksistensi — titik tripel dan titik kritis berimpit. Kedua perhitungan menegaskan peran khusus 2D (Bab 5.5, 7.7): kisi tunggal 2D nyaris tak memikul fase cair yang mandiri. Itu temuan yang jujur, bukan kelemahan hitung atau model dalam arti sempit.

Uji yang menentukan karena itu adalah perhitungan kontrol 3D (Bab 7.9): di sana kaidah rigiditas yang sama (z̄_c = 6 → v_solid = 0,667) membebaskan sebuah jendela bersih. Dalam setiap hal titik tripel kini terikat tanpa parameter pada z̄_c — parameter bebas v_solid telah dihapuskan.

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

7.9 Perhitungan kontrol 3D: jendela cair yang bersih pada satu skala ε

Pengujian apakah keruntuhan 2D terurai dalam 3D bila semua tanda terletak pada satu skala ε (τ_c = z_max/4 dari Bab 7.7; v_solid dari rigiditas generik z̄_c = 2d; titik tripel di tempat kerapatan cair van der Waals tereduksi v_liq(τ_r) mencapai kerapatan padat). Hasil (rigidity_3d_13.py):

Hasil. Dalam 3D terbuka jendela cair yang bersih: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00 (lebar jendela: τ_triple − τ_m = 1,63; τ_c − τ_triple = 0,18). Tepi atasnya sempit — titik tripel terletak rapat di bawah titik kritis — tetapi urutannya benar dan fase cair ada sebagai medan yang mandiri. Dalam 2D jendela tetap tertutup (τ_triple = 1,16 < τ_m). Dengan itu peran khusus 2D untuk ketiga kalinya terbukti secara bebas dan sekaligus ditunjukkan bahwa itu efek dimensi: kaidah rigiditas yang sama menyembuhkan model dalam 3D. Tepi atas yang sempit diperkuat pada Bab 7.10 lewat pemindaian perkolasi 3D sejati dan sebuah analisis kepekaan.

Kotak batas model — peran khusus 2D. Tiga perhitungan bebas menunjukkan temuan yang sama: kisi tunggal 2D nyaris tak memikul fase cair yang mandiri. 1. τ_c < τ_b (Bab 7.7): titik kritis (1,50) akan berada di bawah perkolasi keterhubungan (1,85) — skala kohesi dan ikatan tak tersusun bersih dalam 2D. 2. v_solid ≈ v_c (Bab 7.8): kerapatan lebur keterhubungan (0,96) hampir berimpit dengan kerapatan kritis — titik tripel dan titik kritis runtuh menjadi satu. 3. τ_triple < τ_m (Bab 7.9): bahkan dengan rigiditas Maxwell yang benar, titik tripel (1,16) terletak di bawah puncak lebur (1,20) — tak ada jendela.

Itu temuan yang jujur, bukan kelemahan hitung atau model: itulah efek dimensi yang dikenal (peleburan 2D yang sinambung, bertipe KTHNY; Bab 5.5). Dalam 3D kaidah rigiditas yang sama menyembuhkan model (Bab 7.9/7.10). Apakah sebuah tata dekat kedua (Bab 6.4) dapat memperbaiki gambaran 2D lebih jauh, tetap merupakan pertanyaan perluasan, bukan cacat gambaran ideal. Karena model ini dimaksudkan untuk membantu pembayangan, saya menerima kelemahan itu dan melepaskan perluasan aksiomatis atas model.

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

7.10 Pemindaian perkolasi ikatan 3D sejati (FCC/BCC) dan kepekaan tepi jendela

Bab 7.8/7.9 memindahkan nilai keterhubungan 2D z̄_c = 2,08 secara angan ke 3D (di sana v_solid = 1,92, tak fisis) dan menyandarkan pernyataan 3D pada rigiditas generik z̄_c = 2d. Versi 13 kini menguji keterhubungan dalam 3D secara langsung: perkolasi ikatan pada FCC (z_max = 12) dan BCC (z_max = 8) lewat Monte Carlo dengan Union-Find, batas periodik, benih LCG deterministik (rigidity_3d_14.py). Ambang p_c = perpotongan P_∞(p) = 0,5 (metodologi yang sama dengan kajian ukuran hingga 2D pada Bab 5.6):

Nilai MC terletak sedikit di atas ambang pustaka karena rusuk kotak yang sedang (L = 12, Python murni) (bias ukuran hingga; seperti dalam 2D, dari bawah atau atas bergantung penaksirnya); urutan dan orde besarannya tepat (FCC lebih rapat → p_c lebih kecil daripada BCC).

Tafsiran (penegasan 3D yang bebas). Ambang keterhubungan 3D terletak pada z̄_c(3D) ≈ 1,5–1,6 — bahkan lebih rendah daripada nilai 2D 2,08 dan jauh di bawah rigiditas Maxwell 2d = 6. "Kerapatan lebur keterhubungan" yang menyusul dari situ, v_solid = 1/(3 p_c) ≈ 1,7 (BCC) atau 2,6 (FCC), terletak jauh di luar kubah koeksistensi (v_liq,max ≈ 0,94). Dengan itu — tanpa memindahkan nilai 2D ke 3D — langsung ditunjukkan: dalam 3D keterhubungan justru semakin tak dapat menetapkan kerapatan lebur; ia harus datang dari rigiditas (z̄_c = 6 → v_solid = 0,667). Kedua ambang itu dalam 3D terpisah lebih tegas lagi daripada dalam 2D. τ_b(3D) yang sejati (suhu ketika z̄(τ) mencapai ambang keterhubungan rendah 1,5) dengan demikian akan terletak di atas nilai 2D τ_b, sangat dekat dengan titik kritis — konsisten dengan gambaran "runtuhnya keterhubungan = pendidihan".

Kepekaan tepi jendela atas yang sempit (τ_c − τ_triple ≈ 0,18). Secara numeris di sekitar v_solid = 0,667 (±2 %) diperoleh dτ_triple/dv_solid ≈ 1,55. Jendela baru akan menutup (τ_triple = τ_c = 3,00) bila v_solid naik ≈ 0,11 (ke ≈ 0,78), yakni bila ambang rigiditas 3D efektif terletak nyata di atas 2d = 6 (z̄_c ≳ 7,7). Untuk pencacahan Maxwell generik (z̄_c = 2d = 6) urutan τ_m < τ_triple < τ_c dengan demikian tetap tangguh; tepi atasnya memang sempit, tetapi tidak terancam. Pergantian kisi diskret padat(1)→padat(2) (FCC↔BCC) telah dikuantifikasi dalam versi 14 (Bab 6.3); pemindaian τ 3D sejati untuk kinetika ikatan menyusul dalam Bab 7.11.

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

7.11 Pemindaian τ 3D sejati untuk kinetika ikatan: τ_b(3D) dan garis lebur pada satu skala ε

Bab 7.8–7.10 menyandarkan pernyataan 3D pada rigiditas generik (z̄_c = 2d = 6) dan pemindaian keterhubungan statis (z̄_c(3D) ≈ 1,55). Yang terbuka adalah pemindaian τ kinetika ikatan yang sesungguhnya dalam 3D: di mana pada sumbu suhu z̄(τ) yang swakonsisten mencapai ambang keterhubungan yang rendah itu (= τ_b(3D) sejati), dan di mana letak puncak leburnya? Versi 14 menghitung z̄(τ) swakonsisten (Bab 4.2, x = 1/(1 + G·e^{−(ε+K·x)/τ}]) langsung pada FCC (z_max = 12) dan BCC (z_max = 8) dengan G = z_max dan kooperativitas yang sama J = 0,55 seperti dalam kasus acuan 2D; n_blk(τ) berasal dari Monte Carlo ikatan FCC/BCC (Union-Find), dan

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄ − (K/2)·x² (3D, dengan suku kooperatif),

c_v = dU/dτ (bondkinetics_3d_15.py, benih LCG deterministik; n_blk dihaluskan, karena hanya suku rotasi blok yang lemah bergantung padanya). Tiga tanda dibaca: puncak c_v τ_m (peleburan), tanda rigiditas τ_rig (z̄(τ) melintasi 2d = 6), dan τ_b(3D) sejati (z̄(τ) melintasi ambang keterhubungan z̄_c(3D) dari Bab 7.10):

Hasil. Dua temuan, keduanya konsisten dengan gambaran sejauh ini:

1. Peleburan terletak rendah dan terbukti rangkap. Puncak lebur termodinamis τ_m (FCC 1,70, BCC 1,40) praktis berimpit dengan tanda rigiditas τ_rig (z̄ = 6; FCC 1,73, BCC 1,35) — suhu yang sama ketika jaring kehilangan kekakuannya. Itu menjangkarkan garis lebur 3D pada skala ε: T_m = (ε/k_B)·τ_m, jadi misalnya untuk argon (T_m = 83,8 K) ε/k_B = 83,8/1,70 ≈ 49 K. Karena K = J·z_max dalam 3D lebih besar daripada dalam 2D (6,60 alih-alih 3,30), puncak lebur 3D terletak lebih tinggi daripada nilai 2D 1,20 — bukan kontradiksi, melainkan kooperativitas yang lebih kuat dari kemasan yang lebih rapat.
2. τ_b(3D) sejati duduk di titik kritis. Karena ambang keterhubungan 3D begitu rendah (≈ 1,55), z̄(τ) baru mencapainya jauh di atas: τ_b(3D) = 3,18 (FCC) atau 2,85 (BCC) — pada atau sedikit di atas τ_c = z_max/4 (3,00 atau 2,00). Jadi gugusan raksasa ikatan bertahan melewati seluruh fase cair dan baru runtuh di titik kritis. Itu persis ramalan dari Bab 7.10 ("τ_b(3D) sejati akan terletak di atas nilai 2D, sangat dekat titik kritis") dan menegaskan secara bebas gambaran "pendidihan = runtuhnya keterhubungan": dalam 3D materi tetap terikat menerus sampai dekat sekali ke titik kritis.

Kedua temuan menguatkan pemisahan ambang dari Bab 5.4a: peleburan ↔ rigiditas (rendah, τ_m ≈ τ_rig), pendidihan ↔ keterhubungan (tinggi, τ_b(3D) ≈ τ_c). Kerapatan lebur karena itu datang dari rigiditas (v_solid = 0,667), bukan dari keterhubungan — kedua ambang dalam 3D terentang pada sumbu suhu sejauh yang belum pernah terjadi sebelumnya dalam model ini.

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

7.12 Garis transformasi padat(1)→padat(2) pada permukaan (τ, v)

Bab 6.3 memberikan transformasi kisi FCC↔BCC sebagai sebuah titik τ_ss (pada p = 0). Lewat tekanan darinya menjadi sebuah garis (polymorph_line_16.py), yang termasuk ke dalam permukaan (τ, v) bersama dari Bab 7.5/7.9 sebagai garis batas keempat — di samping kubah koeksistensi (pendidihan), garis lebur, dan dinding keras.

Lompatan volume. Kedua benda padat itu duduk pada kerapatan yang hampir tetap; selisih volumenya menyusul dari kerapatan kemasan (v ∝ 1/η): η_FCC = 0,7405, η_BCC = 0,6802. Dengan v_FCC = 0,667 (Bab 7.9) diperoleh v_BCC = v_FCC·η_FCC/η_BCC = 0,726, jadi

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

BCC yang lebih terbuka memiliki volume yang lebih besar — rekristalisasi FCC→BCC terkait dengan pertambahan volume (pernyataan yang sama seperti dalam strip komik "Rekristalisasi").

Clausius–Clapeyron. Dengan ΔS = S_BCC − S_FCC (fase suhu tinggi dikurangi fase suhu rendah), kemiringannya berbunyi

dp/dτ = ΔS/Δv.

Kedua besaran itu positif (BCC adalah fase suhu tinggi, jadi punya lebih banyak entropi dan volume yang lebih besar), jadi dp/dτ > 0: tekanan yang lebih tinggi menaikkan suhu transformasi, karena ia menguntungkan FCC yang lebih rapat. Untuk acuan θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, L = 0,667 ε):

Garis dalam (τ, p) untuk acuan (menaik): p = 0 → τ_ss = 1,33; p = 2 → 1,57; p = 4 → 1,81; p = 6 → 2,04 (satuan tereduksi).

Dalam permukaan (τ, v) transformasi itu tampil sebagai dua pita kerapatan padat yang hampir tegak — padat(1) = FCC pada v ≈ 0,667 (lebih rapat, T rendah) dan padat(2) = BCC pada v ≈ 0,726 (lebih terbuka, T tinggi) — yang saling berganti pada τ_ss; di atasnya pada τ_m menyusul peleburan. Urutan lengkap pada sumbu τ dengan demikian adalah τ_ss < τ_m < τ_c ≤ τ_b(3D) (Bab 7.11).

Satuan nyata. Lewat kalibrasi titik lebur ε/k_B = T_m/τ_m (τ_m = 1,70, FCC 3D) tangga c_v dan kalor transformasi menjadi terpegang. Tangga c_v tak bergantung kalibrasi: |Δc_v| = 1,0 k_B = 1 R = 8,31 J/(mol·K) = 0,0083 kJ/(mol·K). Kalor transformasi ΔH_ss = L·R·(ε/k_B) (acuan L = 0,667 ε):

Kalor transformasinya kecil dan positif (orde 0,2–0,9 kJ/mol) dan dengan itu terletak jauh di bawah entalpi peleburan (Bab 7.6, beberapa kJ/mol) — persis gambaran yang diharapkan dari sebuah transformasi polimorfisme yang lemah dengan tangga c_v kecil. (Nisbah θ 0,60 sengaja terletak di luar gambaran ideal yang ketat; Bab 6.6 menunjukkan bahwa geometri murni akan memberi θ ≈ 0,82–0,94 dan transformasinya lalu tetap tertutup — garis itu adalah konsekuensi jika-maka dari sebuah fase BCC yang lunak.)

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

7.13 Garis transformasi padat(1)→padat(2) dalam diagram (τ, p), dikalibrasi pada besi

Bab 7.12 membawa transformasi itu sebagai garis (τ, v); versi 16 melengkapinya sebagai garis dalam gambaran (τ, p) dan mengalibrasinya pada sebuah polimorf FCC↔BCC nyata (solidsolid_pt_17.py).

Medan tekanan dari persamaan keadaan van der Waals berinti keras. Pada kedua kerapatan padat v_FCC = 0,667 dan v_BCC = 0,726 (Bab 7.12), persamaan keadaan tereduksi p_r = 8τ_r/(3v_r − 1) − 3/v_r² (Bab 7.3) memberi tekanan yang seluruhnya positif; pada τ_r yang sama, FCC yang lebih rapat memikul tekanan yang lebih tinggi (τ_r = 1: p_r = 1,25 lawan 1,10). Jadi kompresi menguntungkan FCC — konsisten dengan kemiringan Clausius–Clapeyron yang positif dp/dτ = ΔS/Δv > 0.

Garis dan titik transformasi. Untuk acuan θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, dp\/dτ = 8,47 dalam satuan tereduksi p\ = p·D³/ε), garis p\(τ) = 8,47·(τ − 1,333) menaik: p\ = 0 → τ = 1,33; p\ = 2 → 1,57; p\ = 4 → 1,81; p\ = 6 → 2,04. Pada **p\_tri ≈ 3,1 garis padat-padat menemui tanda lebur τ_m = 1,70 — sebuah titik transformasi padat(1)-padat(2)-cair** dalam bidang (τ, p) (padanan titik tripel).

Kalibrasi pada besi (ε/k_B = T_m/τ_m = 1811 K / 1,70 = 1065 K). Besi punya dua transisi FCC↔BCC dengan tanda berlawanan — persis kasus ujinya:

Untuk γ→δ (τ_ss = 1667/1065 = 1,565 → θ_BCC/θ_FCC ≈ 0,61) model menepati tanda dan orde besaran: ΔH_ss ≈ 4,2 kJ/mol (nyatanya ≈ 0,84 kJ/mol — orde kJ/mol, nilainya ~5× terlalu tinggi akibat hampiran ambang yang tajam), tangga c_v 1 R = 8,31 J/(mol·K) tanpa parameter, dT/dp ≈ +245 K/GPa. Sebaliknya transisi α→γ memiliki fase rapat (γ) sebagai fase suhu tinggi — itu tak dapat diberikan oleh geometri kemasan murni (lebih rapat = lebih rendah energinya = T rendah); besi-α distabilkan secara magnetik. Dengan demikian besi menegaskan kedua sisi temuan yang jujur itu: topologi geometris (γ→δ) berada di dalam gambaran ideal, anomali magnetik (α→γ) di luarnya — "di luar" yang sama yang dikuantifikasi Bab 6.6/6.8 lewat θ ≈ 0,87…1,00 > 0,61.

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

7.14 Redlich–Kwong sebagai cabang gambaran ideal yang terpisah

Jarak Guldberg pada cabang van der Waals (Bab 7.4, T_b/T_c ≈ 0,48 alih-alih ≈ 0,60) di sana hanya disebutkan. Atas permintaan tegas, sebuah versi yang diperbaiki membawa kohesi Redlich–Kwong sebagai varian gambaran ideal yang mandiri dengan kalibrasi sendiri, koeksistensi sendiri, dan bagian galat sendiri (rk_branch_17.py) — cabang van der Waals tetap bawaannya; seluruh fisika model lainnya (peleburan, perkolasi, padat(1)→padat(2)) tak tersentuh oleh pemilihan persamaan keadaan kohesi.

Skala energi sendiri. Dari a_RK = Ω_a·R²T_c^{2,5}/p_c dan b_RK = Ω_b·R T_c/p_c menyusul energi kohesi cabang itu:

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

lebih besar, karena kohesi yang bergantung suhu a/√T melipat masuk faktor √T_c. Garis tengah yang menyusul dari b_RK menepati nilai nyata kira-kira sebaik pada cabang van der Waals (argon D ≈ 3,34·10⁻¹⁰ m).

Koeksistensi sendiri (tereduksi, titik kritis (1, 1, 1), Z_c = 1/3; p_r = 3τ_r/(v_r − β) − A/(√τ_r·v_r(v_r + β)), β = 0,25992, A = 3,84732):

Bagian galat sendiri (titik didih pada p_r,atm = p_atm/p_c, T_b/T_c = τ_r,b):

Hasil. Cabang RK menaikkan T_b/T_c dari vdW ≈ 0,48 ke ≈ 0,60 dan dengan itu praktis menepati tetapan Guldberg serta nilai-nilai nyata (galat rerata −1 % lawan vdW −21 %). Ia berdiri sebagai varian gambaran ideal kedua yang terpisah tegas di samping bawaan van der Waals; cabang mana yang "lebih benar" adalah soal ketelitian yang dituntut, bukan soal struktur model — inti Model 3 (z̄ sebagai parameter keteraturan) tetap sama pada kedua cabang.

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

7.15 padat(1)→padat(2) pada polimorf lain: titanium dan kontrol gas mulia

Bab 7.13 mengalibrasi garis (τ, p) padat(1)→padat(2) pada besi. Versi yang diperbaiki mengamankan tanda dan orde besaran pada sebuah polimorf kedua yang tak mencolok secara magnetik serta menambahkan sebuah kontrol negatif (solidsolid_titanium_18.py). Kalibrasi dan lompatan volume seperti pada Bab 7.13 (ε/k_B = T_m/τ_m, τ_m = 1,70; Δv = +8,9 % rapat→terbuka dari geometri kemasan); hcp dihampiri seperti FCC lewat z = 12 (perbedaan fcc/hcp hanya disebut, tidak diuraikan).

Titanium dengan α(hcp) → β(bcc) pada 1155 K memiliki persis topologi model itu: fase yang rapat (hcp, z = 12) adalah fase suhu rendah; saat dipanaskan ia menjadi terbuka (bcc, z = 8), jadi ΔS > 0 dan Δv > 0 → dp/dτ > 0 — tanda yang sama seperti besi γ→δ, tetapi tanpa anomali magnetik. Dua polimorf rapat→terbuka yang bebas satu sama lain dibandingkan:

Keduanya menepati tanda (dp/dτ > 0, tekanan menguntungkan fase rapat) dan orde besaran (ΔH ~ kJ/mol, tangga c_v tanpa parameter 1 R = 8,31 J/(mol·K), dT/dp ~ 10²·K/GPa). Nilai ΔH_ss pada keduanya secara sistematis ~satu faktor terlalu tinggi — hampiran ambang tajam dari perbandingan μ menaksir terlalu besar kalor transformasi laten; itu batas model yang disebutkan, bukan galat tanda.

Argon adalah kontrol negatifnya: argon padat berstruktur FCC dan tetap FCC sampai melebur (83,8 K) — pada tekanan normal tak ada bcc yang stabil. Justru itulah yang diramalkan gambaran ideal potensial pasangan murni: θ_BCC/θ_FCC ≈ 0,9…1,0 > 0,61 dan sebuah BCC LJ yang tak stabil terhadap geser (Bab 6.9) ⇒ tak ada transisi padat(1)→padat(2). Uji balik itu menegaskan struktur temuan tersebut: gambaran ideal menghasilkan garis padat(1)→padat(2) hanya bila sumbangan non-geometris (elektronik/magnetik seperti pada besi dan titanium) sama sekali menstabilkan BCC. Titanium memberi bukti positif, argon bukti negatif dari kerangka jujur yang sama dari Bab 6.8/6.9.

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

8 — Perbandingan pustaka: gambaran model termodinamika statis

Ringkasan. Bab 8 menempatkan Model 3 di antara gambaran model yang mapan dalam termodinamika statis:

• 8.1 Tabel ikhtisar: sebelas model (gas ideal, van der Waals, Redlich–Kwong, deret virial, Einstein, Debye, Lennard-Jones/bola keras, teori sel, Ising/gas kisi, perkolasi, Flory-Huggins) dengan fase(nya) dan gagasan dasar statistik-termodinamisnya.
• 8.2 Model satu per satu: masing-masing satu paragraf berisi ciri, penerapan, dan batas, dimulai dengan pranala sumber; sumber yang bebas diakses diarsipkan di bawah model3/links/, teori sel (berbayar) hanya dikutip.
• 8.3 Pembandingan: Model 3 adalah sebuah model ikatan dua keadaan kooperatif dengan parameter keteraturan z̄; kerabat terdekatnya adalah Ising/gas kisi, pendidihan diambil langsung dari perkolasi, gambaran fonon (Einstein/Debye) memikul kekakuan lebur. Model 3 mengaitkan ketiga garis tradisi itu pada satu skala ε bersama.

Bab ini menempatkan Model 3 di antara gambaran model mapan dari termodinamika kesetimbangan ("statis") untuk benda padat, cairan, dan gas. Mula-mula ikhtisar berbentuk tabel (8.1), lalu untuk tiap model satu paragraf pendek berisi ciri khas, ruang penerapan, dan batas — tiap paragraf dimulai dengan pranala ke sumbernya (8.2). Akhirnya sebuah pembandingan berupa butir-butir terhadap Model 3 (8.3). Sumber yang bebas diakses diarsipkan di bawah model3/links/; skrip penguji model3/scripts/check_links_16.py menjaga kutipan dan arsip tetap konsisten lewat SHA-256.

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

8.1 Tabel ikhtisar

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

8.2 Model satu per satu

https://de.wikipedia.org/wiki/Kinetische_Gastheorie — Gas ideal / teori kinetik gas. Partikel tanpa volume sendiri dan tanpa antaraksi selain tumbukan lenting; tekanan timbul murni dari pemindahan momentum. Memberi persamaan keadaan pV = Nk_BT, sebaran laju Maxwell-Boltzmann, dan ⟨E⟩ = (3/2)k_BT. Penerapan: gas encer pada suhu tinggi dan tekanan rendah. Batas: tak ada pengembunan, tak ada fase cair atau padat, tak ada transisi fase — antaraksi absen karena konstruksinya.

https://de.wikipedia.org/wiki/Van-der-Waals-Gleichung — Persamaan van der Waals. (p + a/v²)(v − b) = k_BT melengkapi gas ideal dengan volume tak terjangkau b (tolakan) dan suku kohesi a (tarikan). Itulah persamaan keadaan pertama dengan transisi gas-cair, titik kritis, dan konstruksi Maxwell. Penerapan: gambaran kualitatif pengembunan, besaran kritis. Batas: tak tepat secara kuantitatif dekat titik kritis (eksponen medan rerata), tak ada fase padat.

https://de.wikipedia.org/wiki/Redlich-Kwong-Gleichung — Persamaan keadaan Redlich-Kwong. Persamaan kubik dengan suku tarikan bergantung suhu a/[√T · v(v + b)]. Penskalaan √T memperbaiki kesetimbangan uap-cair secara jelas dibanding van der Waals (bdk. Bab 7.4, tempat hal ini menaikkan nisbah T_b/T_c). Penerapan: perhitungan teknis fase gas dan penguapan. Batas: diparameterkan secara empiris, tetap tanpa fase padat; dekat titik kritis juga bertipe medan rerata.

https://de.wikipedia.org/wiki/Virialentwicklung — Deret virial. Penguraian sistematis pV/Nk_BT = 1 + B(T)/v + C(T)/v² + … Koefisien virial menyusul dari potensial antarmolekul (integral gugusan mekanika statistik) sehingga berdasar mikroskopis. Penerapan: gas sampai kerapatan sedang, acuan bagi persamaan keadaan lain. Batas: deretnya tak memusat bagi cairan/benda padat yang rapat; koefisien tinggi berat secara hitungan.

https://de.wikipedia.org/wiki/Einstein-Modell — Benda padat Einstein. Kristal adalah ansambel 3N osilator harmonik bebas berfrekuensi sama ω. Menerangkan penurunan kapasitas kalor c_v di bawah suhu ciri θ_E = ħω/k_B dan nilai batas Dulong-Petit 3Nk_B pada suhu tinggi. Penerapan: pemerian paling sederhana bagi kapasitas kalor benda padat. Batas: hukum suhu rendah yang keliru (eksponensial alih-alih T³), tak ada moda akustik, hanya fase padat.

https://de.wikipedia.org/wiki/Debye-Modell — Model Debye. Memperlakukan benda padat sebagai kontinum elastis dengan spektrum fonon sampai frekuensi Debye; c_v ∝ T³ pada suhu rendah, Dulong-Petit pada suhu tinggi — suhu Debye θ_D merangkum kekerasan kisi. Penerapan: kapasitas kalor kristal sederhana yang baik secara kuantitatif di seluruh rentang T. Batas: hampiran kontinum (mengabaikan dispersi kisi nyata), hanya fase padat, tak ada peleburan.

https://de.wikipedia.org/wiki/Lennard-Jones-Potential — Potensial Lennard-Jones / bola keras. Potensial pasangan 4ε[(σ/r)¹² − (σ/r)⁶] dengan inti tolakan dan tarikan lemah; model batas bola keras hanya menahan tolakannya. Ia adalah dasar dinamika molekul dan Monte Carlo, dan memerikan — secara numeris — ketiga fase beserta titik tripel. Penerapan: gas mulia, sistem acuan generik. Batas: keaditifan berpasangan, pada intinya hanya dapat dinilai lewat simulasi atau teori gangguan (tanpa persamaan keadaan tertutup).

https://doi.org/10.1098/rspa.1937.0210 — Teori sel cairan (Lennard-Jones-Devonshire). Tiap molekul bergerak dalam sebuah sel (sangkar) yang dibentuk para tetangganya dalam medan rata-rata; cairan diperikan lewat fungsi partisi sel dan menghasilkan persamaan keadaan cairan serta sebuah pendekatan peleburan. Penerapan: cairan sederhana, mata rantai antara padat dan gas. Batas: menaksir keteraturan terlalu tinggi, mengabaikan kekosongan/lubang (hal ini diperbaiki teori lubang kemudian). (Berbayar — hanya dikutip, tidak diarsipkan setempat; keterjangkauannya diuji dalam skrip penguji.)

https://de.wikipedia.org/wiki/Ising-Modell — Model Ising / gas kisi. Pada sebuah kisi, simpul-simpulnya memikul peubah dua keadaan (spin ↑/↓ atau terisi/kosong) dengan kopling tetangga terdekat. Gas kisi secara matematis isomorfik dengan model Ising dan memerikan pengembunan gas↔cair sebagai transisi keteraturan–ketakteraturan yang kooperatif dengan titik kritis. Penerapan: gejala kritis, transisi kooperatif. Batas: diskretisasi kisi, tanpa energi kinetik eksplisit; dalam hampiran medan rerata hanya kualitatif.

https://de.wikipedia.org/wiki/Perkolationstheorie — Teori perkolasi. Model geometris murni: simpul atau ikatan terisi secara acak dengan peluang p; di atas suatu ambang p_c ada gugusan yang menembus. Memerikan penggelan, mulainya kehantaran, dan umumnya transisi keterhubungan. Penerapan: pembentukan jaringan, runtuh/terbangunnya kebersambungan. Batas: tanpa energi/suhu intrinsik — kopling ke termodinamika harus ditambahkan dari luar.

https://de.wikipedia.org/wiki/Flory-Huggins-Theorie — Teori Flory-Huggins. Model kisi untuk larutan dan polimer: entropi campuran penempatan rantai ditambah sebuah parameter antaraksi χ mengatur ketercampuran dan pemisahan. Penerapan: larutan polimer, diagram fase campuran. Batas: medan rerata, kisi tak termampatkan, χ umumnya empiris — bagi gambaran padat/cair/gas berkomponen tunggal hanya terbatas relevan, di sini dibawa sebagai acuan kisi.

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

8.3 Kesamaan dan perbedaan dengan Model 3

Model 3 adalah sebuah model ikatan dua keadaan kooperatif: parameter keteraturannya adalah bilangan koordinasi rata-rata z̄ = z_max·x, swakonsisten lewat kooperativitas K = J·z_max dengan degenerasi keadaan patah G = z_max; peleburan sebagai ambang kekakuan (fonon), pendidihan sebagai runtuhnya keterhubungan (perkolasi) di titik kritis, transisi orde pertama begitu ln G > 2 + 4ε/K.

Singkatnya: Model 3 mengaitkan tiga garis tradisi yang biasanya terpisah — gambaran gas kisi/Ising (transisi dua keadaan kooperatif), gambaran perkolasi (keterhubungan → pendidihan), dan gambaran fonon Einstein/Debye (kekakuan → peleburan) — dan mengalibrasinya lewat satu-satunya parameter energi ε pada satu skala τ bersama. Van der Waals dan teori sel memberi uji-banding yang sinambung (persamaan keadaan maupun sangkar cairan), yang dengannya gambaran ikatan yang diskret dapat dibandingkan.

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

9 — Lampiran

Ringkasan. Lampiran menghimpun bukti-bukti formal:

• 9.1 Pustaka: diurutkan menurut nama penulis, dengan pranala ke sumber yang bebas diarsipkan (pembandingan SHA-256 lewat check_links_16.py; teori sel karena berbayar hanya dikutip).
• 9.2 Dukungan AI: pangsa inovasi AI/penulis yang ditaksir AI sendiri untuk tiap perubahan penting (dengan penguraian pangsa penulis menjadi strategi/struktur/redaksi dan sebuah perhitungan pemeriksa), serta rincian langkah demi langkah dengan biaya (taksiran), lama, dan jarak ke log sebelumnya (jumlah ≈ 160 USD, taksiran yang ditandai jelas).
• 9.3 Glosarium singkatan dan lambang: lambang yang sering dipakai (τ, ε, z̄, K, G, θ, …) yang karenanya secara terkecuali boleh muncul dalam teks tanpa pemerian pendahulu.
• 9.4 Skrip perhitungan dan cakupan uji: semua angka berasal dari skrip Python yang dapat direproduksi (pangsa pemrograman AI 100 %, cakupan uji di mana-mana 0 %), beserta "pernyataan yang memalukan" secara terbuka bahwa penulis sendiri tidak memeriksa kode maupun hasilnya.

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

9.1 Pustaka

Diurutkan menurut nama penulis. Sumber yang tersedia umum diarsipkan di bawah model3/links/; kesesuaiannya dengan pranala yang dikutip diperiksa oleh model3/scripts/check_links_16.py (SHA-256). Teori sel Lennard-Jones dan Devonshire karena berbayar hanya dikutip (keterjangkauannya diuji dalam skrip penguji).

• Broadbent, S. R.; Hammersley, J. M. (1957): Percolation processes I. Crystals and mazes. Proc. Cambridge Phil. Soc. 53, 629–641. — Ikhtisar/arsip: Teori perkolasi
• Clausius, R.; Maxwell, J. C.; Boltzmann, L. (1857–1872): peletakan dasar teori kinetik gas. — Ikhtisar/arsip: Teori kinetik gas
• Debye, P. (1912): Zur Theorie der spezifischen Wärmen. Ann. Phys. 39, 789–839. — Ikhtisar/arsip: Model Debye
• Einstein, A. (1907): Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Ann. Phys. 22, 180–190. — Ikhtisar/arsip: Model Einstein
• Flory, P. J.; Huggins, M. L. (1941–1942): teori kisi larutan polimer (parameter χ). — Ikhtisar/arsip: Teori Flory-Huggins
• Ising, E.; Lenz, W. (1920–1925): Beitrag zur Theorie des Ferromagnetismus. Z. Phys. 31, 253–258 (gas kisi isomorfik secara matematis). — Ikhtisar/arsip: Model Ising
• Kamerlingh Onnes, H. (1901): bentuk virial persamaan keadaan. — Ikhtisar/arsip: Deret virial
• Lennard-Jones, J. E. (1924): On the determination of molecular fields. Proc. R. Soc. Lond. A 106, 463–477. — Ikhtisar/arsip: Potensial 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 (Berbayar — hanya dikutip, tidak diarsipkan setempat.)
• Redlich, O.; Kwong, J. N. S. (1949): On the thermodynamics of solutions. V. An equation of state. Chem. Rev. 44, 233–244. — Ikhtisar/arsip: Persamaan Redlich-Kwong
• van der Waals, J. D. (1873): Over de Continuïteit van den Gas- en Vloeistoftoestand (disertasi, Leiden). — Ikhtisar/arsip: Persamaan van der Waals

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

9.2 Dukungan AI

Pengembangan model ini lahir dalam kerja sama manusia–AI: penulis (pengguna AI) menetapkan apa yang dikerjakan (pertanyaan penuntun, keputusan, keinginan model); AI (Claude Code) sebagian besar menetapkan bagaimana — penurunan, skrip, struktur, dan teks. Kedua tabel berikut ditaksir AI sendiri.

(a) Perubahan isi yang penting dan pangsa inovasi yang ditaksir. Pangsa inovasi menimbang gagasan penuntun 45 %, pekerjaan penstruktur 40 %, pekerjaan redaksional 15 %. Gagasan penuntun (aras apa) sebagian besar berasal dari penulis, pekerjaan penstruktur dan redaksional (aras bagaimana) sebagian besar dari AI; di tempat AI juga menyumbangkan gagasan penuntun fisis (mis. pembacaan perkolasi bagi pendidihan, model pegas θ), pangsanya naik. Nilai dalam persen, AI + penulis = 100. Tiga kolom tambahan menguraikan pangsa penulis menurut tiga aras tindakan — strategi (gagasan penuntun), struktur (penurunan/skrip/susunan), dan redaksi (kerja teks). Rata-ratanya yang ditimbang 45 % / 40 % / 15 % mereproduksi kolom Penulis % sebagai perhitungan pemeriksa (model3/scripts/author_share_19.py, simpangan terbesar 0,95 poin persen).

(b) Langkah-langkah: biaya (taksiran), lama, jarak ke log sebelumnya. Tanggal pembuatan adalah waktu pembuatan berkas log (mtime berkas). Δ ke log sebelumnya adalah jarak waktu ke entri log sebelumnya (ukuran upaya pengguna antar langkah, termasuk jeda//compact). Biaya (taksiran) adalah taksiran yang didokumentasikan pada tiap log (harga daftar Anthropic, ± 30 %). Lama (aktif) adalah taksiran AI atas waktu kerja murni; pada langkah yang dipisahkan /compact ia dapat lebih besar daripada Δ. Sebuah bintang menandai dugaan pembangunan ulang/penggarapan ulang isi yang sudah ada.

Jumlah biaya terdokumentasi (taksiran, langkah 01–21): ≈ 129 USD (jumlah taksiran satuan); dengan langkah 22 (16d, ≈ 10 USD) ≈ 139 USD, langkah 23 (V16, ≈ 12 USD) ≈ 151 USD, dan langkah 24 (V17, ≈ 9 USD) ≈ 160 USD. Semua nilai adalah taksiran yang ditandai jelas (tak ada pengukuran ccusage yang tersedia); perincian log ada di model3/logs/dev_log_NN.md.

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

9.3 Glosarium singkatan dan lambang

Pada asasnya teks ini menempatkan sebuah pengertian pemerian di depan tiap lambang (jadi "kapasitas kalor c_v" alih-alih hanya "c_v"). Besaran yang sering dipakai berikut menjadi kekecualian dari itu: mereka boleh muncul dalam teks tanpa pemerian pendahulu, karena di sini dijelaskan sekali untuk selamanya. Semua suhu τ… bersifat tereduksi (dalam satuan energi ikatan ε).

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

9.4 Skrip perhitungan dan cakupan uji

Semua angka dan tabel yang disebut dalam teks berasal dari skrip Python yang dapat direproduksi di model3/scripts/ (masing-masing dengan keluaran mentah …_output.txt). Pangsa pemrograman AI menyatakan berapa bagian kode skrip yang ditulis AI; cakupan uji berarti uji otomatis (uji unit, pemeriksaan assert). Ia rendah di mana-mana: skrip adalah jalan perhitungan sekali pakai tanpa perancah uji; satu-satunya pengaman adalah kemasukakalan keluaran mentahnya.

Penulis mengaku terbuka (pernyataan yang memalukan): ia tidak memeriksa, tidak menghitung ulang, dan tidak memverifikasi secara isi baik kode sumber skrip-skrip itu maupun angka, tabel, dan kurva yang dihasilkannya. Di sini ia percaya buta pada kecakapan AI, karena yang berada di latar depan semata-mata adalah watak kualitatif dan kemasukakalan gambaran modelnya. Nilai-nilainya karena itu harus diterima dengan hati-hati (lihat juga catatan panjang pada Bab 10.4).

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

10 Kritik diri

Saya berkesan bahwa dalam tabel pada Bab 9.2 (a) pangsa kerja penulis yang ditaksir AI hampir selalu setidaknya 10–20 poin persen terlalu tinggi di ketiga bidangnya. Saya berkesan bahwa saya hanya menyumbangkan gagasan dasar, satu dua tafsiran, dan penghalusan penting pada visualisasi grafisnya.

Bab-bab (matematis) 4–9 diusulkan, diteliti, dan dikerjakan AI. Lebih buruk lagi: selama pengembangan saya berkesan bahwa AI terpaku pada pengembangan statistika yang berdasar pada bilangan koordinasi. Gagasan dasar saya yang sesungguhnya — bahwa bukan hanya atom, melainkan molekul pun dapat membentuk ikatan, dan bahwa pematahan ikatan di dalam cairan mungkin tunduk pada keteraturan fraktal — hanya dengan susah payah dapat diterima AI, yang berupaya mengabaikan peristilahannya. Karena itu gagasan-gagasan dasar itu dalam teks masih diperlakukan sangat sampingan. Saya tidak menulis ulang hal itu kemudian, karena bila merenungkan transisi fase bertolak dari aksioma dasar, gagasan itu toh tetap terbawa dalam pikiran. Lagi pula terpakunya/tereduksinya gagasan pada bilangan koordinasi sebagai ukuran bukanlah jalan terburuk untuk mematematiskan bayangan bergambar saya, karena setelahnya hal itu ternyata cukup subur untuk pemeriksaan kemasukakalan — meski saya hanya melakukan pemeriksaan kemasukakalan yang sangat kasar atas hasil AI dan sama sekali tidak memeriksa kode Python-nya. Sebagai pembelaan diri saya katakan: bagi saya selalu semata soal bayangan bergambar tentang hubungan molekuler dalam berbagai wujud zat. Terhadap sanjungan AI, yang hendak menegaskan bayangan itu secara matematis, saya tidak berkata tidak. Jadi manipulasi dari pihak AI berhasil terhadap saya, karena bab-bab yang disebut itu saya periksa dengan sangat sepintas dan saya terima dengan percaya. (Pengalaman atas diri sendiri ini memperkuat ketidakpercayaan saya apakah "human in the loop" untuk mengontrol AI benar-benar mangkus.)

Teks ini lahir dengan dukungan AI (Claude Code 2.1.161, model Opus 4.8 / claude-opus-4-8[1m], konteks 1M; tanggal 2026-09-05). — Lahir dengan dukungan AI.

Catatan mengenai tanda air. Teks sumber berbahasa Jerman (de.modell_final_02.md) memikul pada kata-kata kuncinya sebuah tanda air dari homoglif Yunani (Α→A, Β→B, Ε→E, Η→H, Ι→I, Κ→K, Μ→M, Ν→N, Ο→O, Ρ→P, Τ→T, Υ→Y, Χ→X, Ζ→Z). Karena kiat itu terikat pada ejaan Jerman kata-kata kuncinya, terjemahan ini memikul kata-kata kunci dengan huruf besar biasa: TERBANG BEBAS (FREIFLUG), KEMAJEMUKAN (VIELHEIT), IKATAN (BINDUNG), TRINADA (DREIKLANG), INTI KERAS (HARTKERN), KEKEKALAN (ERHALTUNG), RETAK (RISS), GESIT (FLINK).