İdeal saf madde ve maddenin hâlleri için istatistiksel termodinamik

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

Genel özet

Model 3, katı, sıvı ve gaz olmak üzere üç maddenin hâlini ve bunlar arasındaki geçişleri moleküler bağlardan oluşan tek bir bağ resminden türeten niteliksel, istatistiksel-mekanik bir ideal madde modelidir. Temel ilke, iç enerji için üç enerji deposunun hesaba katılmasında yatar: moleküler bağlardaki titreşim enerjisi, moleküllerin öteleme enerjisi ve iç yoğunluk dalgalanmalarını dengeleyen iç madde akışlarının idealleştirilmiş biçimi olarak misellerin dönmesindeki işbirlikçi hareket enerjisi. Bu bağlamdaki önemli bir temel büyüklük, işbirlikçi bir iki-durumlu bağ üzerinden (Ising benzeri, işbirliği K = J·z_max, kırık durum dejenereliği G = z_max) kendiyle tutarlı biçimde yerleşen ortalama koordinasyon sayısı z̄ = z_max·x'tir. Ondan, tek bir enerji parametresi ε ile kalibre edilmiş ortak bir τ ölçeğinde üç geçiş çıkar: erime bir fonon/sertlik eşiği olarak (z̄ rijitliğe ulaşır), kaynama kritik noktadaki bağlantılılık (perkolasyon) çöküşü olarak ve isteğe bağlı bir katı(1)→katı(2) geçişi (YMK↔HMK), 1,0 k_B'lik parametresiz c_v basamağıyla. Bir geçiş, ln G > 2 + 4ε/K olur olmaz birinci mertebeden olur. Burada kırık durum dejenereliği G, kopartılabilir komşu temaslarının sayısıdır (G = z_max, en büyük koordinasyon sayısı), bağ enerjisi ε her komşu temas başına kazanılan enerjidir ve işbirliği K = J·z_max, Ising benzeri eşleşme gücü J ile z_max'ın çarpımıdır; geçiş, temas çeşitliliği G ne kadar büyük ve işbirliği K bağ enerjisi ε karşısında ne kadar güçlü olursa o kadar kolay birinci mertebeden olur.

Yazı, yedi aksiyomdan (Bölüm 1) çelişki çözümlemesine (Bölüm 2) ve durağan çizgi-roman şeritleriyle resimli bir betimlemeye (Bölüm 3), oradan istatistiksel-mekanik iskelete — işbirlikçi bağ, durum toplamları ve ısı sığaları (Bölüm 4–5) —, üç boyutlu matematiksel akla yatkınlığa (Bölüm 6) ve son olarak kalibrasyon ile hâl denklemine (Bölüm 7) götürür. Düşünceler, birinci mertebeden erimenin tam 3B c_v taramasıyla (Δz̄, gizli ısı ΔU, tepe kayması), fonon sertliği θ(örgü) için bir Einstein yay modeliyle ve Clausius–Clapeyron üzerinden katı(1)→katı(2) dönüşüm çizgisiyle tamamlanır. Ardından θ(örgü), Debye benzeri bir fonon DOS'una inceltilir; sonra katı(1)→katı(2) dönüşümü ek olarak (demire kalibre edilmiş) (τ, p) resminde bir çizgi olarak sokulur ve van der Waals varsayılanının yanına ayrı bir Redlich–Kwong dalı eklenir. Bu arada θ(örgü) bir Lennard-Jones tam örgüsünde (YMK kararlı, ideal HMK kaymaya karşı kararsız), katı(1)→katı(2) topolojisi ise titanyumda (ikinci doğrulayıcı) ve argonda (negatif kontrol) sınanır. Bölüm 8, modeli statik termodinamiğin on bir yerleşik model tasavvuru arasına yerleştirir. Ek, kaynakları (yazarlara göre sıralı) ve YZ desteğini (tahmini yenilik payı, tahmini maliyet ve her geliştirme adımının süresi) belgeler; yardımcı bölümler okuma biçimini, kalan eleştirel noktaları ve ileriye bakışı toplar.

Çekirdek sav. Tek bir işbirlikçi bağ ölçüsü z̄, başka türlü ayrı duran üç gelenek çizgisini — örgü gazı/Ising resmi (işbirlikçi iki-durumlu geçiş), perkolasyon resmi (bağlantılılık → kaynama) ve Einstein/Debye'nin fonon resmi (sertlik → erime) — bilinçli olarak van der Waals benzeri bir ideal resim içinde tutarlı, niteliksel olarak eksiksiz bir faz resmine kenetler; kuramsal yapı doğal olarak kolayca canlandırılabilir, dinamik bir resimden türetilebilir. Modelin asıl gücü kesinlik değil, resimselliktir ve bilimsel model değerini bu belirler.

İçindekiler

• Genel özet
• 1 — Anahtar sözcük adlı aksiyomlar (birleştirilmiş)
• 2 — Çelişki ve boşluk çözümlemesi
• 3 — Resimli betimlemeler
• 3.1 Üç tekil faz
• 3.2 Faz geçişleri
• 3.3 Sıvılarda ses hızı neden T ile düşer
• 3.4 Brown hareketi nasıl doğar — keskinleştirilmiş misel kavramıyla
• 3.5 Misel sınır katmanlarında tepkimeler neden hızlanır
• 3.6 Üç boyutlu sıvı penceresi
• 4 — İstatistiksel-mekanik model
• 4.1 Serbestlik dereceleri ve değişkenler
• 4.2 İşbirlikçi iki-durumlu bağ (Ising benzeri, parametre J)
• 4.3 Enerji dengesi — işbirlikçi terimle
• 5 — Durum toplamları ve ısı sığaları
• 5.1 Çarpanlara ayrılmış durum toplamı
• 5.2 Faz başına ısı sığası (analitik)
• 5.3 İşbirlikçi terim **içeren** sayısal c_v(τ) eğrisi — gizli ısı niceliksel
• 5.4 Sonuç olarak bağlantılılık perkolasyon eşiği (değişmeden doğrulandı) - 5.4.a İki farklı eşik temiz biçimde ayrıldı
• 5.5 Erime geçişinin mertebesi
• 5.6 Perkolasyon eşiğinin sonlu boyut incelemesi
• 5.7 Birinci mertebeden erimenin tam 3B c_v(τ) taraması
• 6 — 3B yorumun matematiksel akla yatkınlığı
• 6.1 Dönme serbestlik derecesinin 2B → 3B değişimi
• 6.2 Farklı kristal örgülerin etkisi — ayırt edici özellik olarak bağlantılılık
• 6.3 Ek geçiş katı(1) → katı(2)
• 6.4 Düşünceler: sıvı suyun ses anomalisi
• 6.5 Seçenek: misel dönmesi yerine fraktal girdap alanları
• 6.6 Mikroskobik fonon sertliği θ(örgü): bir Einstein yay modeli
• 6.7 Yarı kararlılığın matematiksel kestirimi (2B üçgen ↔ kare, 16b)
• 6.8 θ(örgü)'nün Debye benzeri inceltilmesi
• 6.9 Lennard-Jones tam örgüsünden θ(örgü)
• 7 — Kalibrasyon ve hâl denklemi
• 7.1 Kalibrasyon yöntemi (tek enerji parametresi ε) — uyarlandı
• 7.2 Kalibrasyon tablosu
• 7.3 Yoğunluk/hacim ekseni: sert çekirdekli van der Waals hâl denklemi
• 7.4 Niceliksel s/g kalibrasyonu: a → ε, b → D²
• 7.5 Üçlü noktalı ortak (τ, v) yüzeyi
• 7.6 Gizli ısıya karşı erime entalpileri
• 7.7 Enerji ölçeklerinin birleştirilmesi: ε ↔ ε_lg
• 7.8 z̄_c'den parametresiz üçlü nokta
• 7.9 3B kontrol hesabı: tek bir ε ölçeğinde temiz sıvı penceresi
• 7.10 Gerçek 3B bağ perkolasyon taraması (YMK/HMK) ve pencere kenarı duyarlılığı
• 7.11 Bağ kinetiğinin gerçek 3B τ taraması: τ_b(3B) ve tek ε ölçeğinde erime çizgisi
• 7.12 (τ, v) yüzeyinde katı(1)→katı(2) dönüşüm çizgisi
• 7.13 Demire kalibre edilmiş (τ, p) diyagramında katı(1)→katı(2) dönüşüm çizgisi
• 7.14 Ayrı bir ideal resim dalı olarak Redlich–Kwong
• 7.15 Başka polimorflarda katı(1)→katı(2): titanyum ve soy gaz kontrolü
• 8 — Yazın karşılaştırması: statik termodinamiğin model tasavvurları
• 8.1 Genel bakış tablosu
• 8.2 Modeller tek tek
• 8.3 Model 3 ile benzerlikler ve farklar
• 9 — Ek
• 9.1 Kaynakça
• 9.2 YZ desteği
• 9.3 Kısaltma ve simge sözlüğü
• 9.4 Hesaplama betikleri ve test kapsamı

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

1 — Anahtar sözcük adlı aksiyomlar (birleştirilmiş)

Özet. Anahtar sözcüklerle adlandırılmış yedi aksiyom ideal maddeyi belirler: korunan parçacık sayısı, saf klasik hareket, işbirlikçi bir iki-durumlu bağ ve düzen parametresi olarak ortalama koordinasyon sayısı z̄. Bunlar sonraki bütün türetimleri taşır; Almanca özgün metinde anahtar sözcükleri bir homoglif filigranı işaretler (anahtar metnin sonunda).

Kesinleştirilmiş aksiyomlar:

1. SERBEST UÇUŞ (ek aksiyom). Birbirinin yanından bağsız geçen iki molekül arasında ne sürtünme ne de elektrostatik kuvvetler etkir. Momentum ve enerji yalnız çarpışmadaki sert SERT ÇEKİRDEK itmesiyle ve BAĞ'ın kurulması/kopmasıyla değişir. Aksiyom yalnızca klasik mekanik ve enerji dengesi açısından önemlidir: moleküler düzeyde hiçbir yitimli kanal bulunmadığını güvenceye alır — hareket enerjisi asla "yitmez", yalnızca üç ÜÇLÜ AKOR deposuna ya da bağ enerjisine geçer. (Makroskobik sürtünme/ağdalık, mikroskobik bir temel kuvvet olarak değil, kaydırmalı bulmacadan ve ÇATLAK'tan belirerek doğar.)
2. ÇOKLUK. İdeal bir madde çok sayıda benzer molekülden oluşur; makroskobik özellikler termodinamik limitte (N → ∞) istatistiksel olarak betimlenir.
3. BAĞ. Geometrik olarak komşu her iki molekül bir bağ taşıyabilir. Sağlam bir bağ enerjiyi sonlu ε kadar düşürür; kopması tam ε'ya mal olur. Bağ, {sağlam, kopuk} iki-durumlu bir sistemdir.
4. ÜÇLÜ AKOR. İç enerji klasik olarak üç depoya dağılır: (a) moleküllerin ötelenmesi, (b) sağlam bağlardaki titreşim (bağ başına harmonik yay), (c) katı molekül bloklarının dönmesi (bağlantılı bölgelerin kolektif dönme serbestlik derecesi).
5. SERT ÇEKİRDEK. Molekül uzaklığı en küçük σ değerinin altına inerse sert bir itme etkir (öz hacim). BAĞ ile birlikte etkin bir çift potansiyeli çıkar (sert kabuk + sonlu bağ çukuru).
6. KORUNUM. Toplam momentum ve toplam enerji korunur. Tek bir molekül enerjik olarak yansızdır — içinde enerji tamponlamaz; enerji bağlarda (BAĞ) ve kinetik serbestlik derecelerinde (ÜÇLÜ AKOR) yaşar.
7. ÇATLAK. Bağlar çatlak/onarım ilkesine göre bozulur/oluşur: eksik komşu bağ ne kadar çoksa → kopma olasılığı o kadar yüksek; bağlı komşu ne kadar çoksa → onarım olasılığı o kadar yüksek. Kinetik işbirlikçidir.
8. ÇEVİK. Kopma ve onarım hızlı ve düşük aktivasyonlu geçer; bağ alt sistemi her sıcaklıkta yerel dengededir ve molekül hareketine karşı adyabatik olarak ayrılabilir (çarpanlara ayrılmış durum toplamının temeli).

Modelleme tutumu. Geometrik koordinasyon idealleştirilmiş biçimde korunur; yalnız kopmalarla koordinasyon sayısı z değişir. Hareket kaydırmalı bulmacadaki gibi doğar: boşluklar (SERT ÇEKİRDEK öz hacminden ve ÇATLAK kopmalarından) göç eder, bloklar ardından kayar. Dönen, göç eden küçük kristaller biçimindeki bilinçli olarak "yanlış" resim, sezgisel temel olarak korunur.

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

2 — Çelişki ve boşluk çözümlemesi

Özet. Çözümleme modeli iç çelişkiler ve boşluklar bakımından dizgesel olarak sınar ve işlenme durumunu sürdürür; saptanan zayıf noktalar adlandırılır ve sonraki bölümlerde ya çözülür ya da dürüstçe bir model sınırı olarak işaretlenir.

Eleştiri noktaları listesi

1. Ötelenmenin iki kez sayılması (aksiyom 3a'ya karşı 3c). Ham metin ötelenmeyi iki kez anar (molekül ötelenmesi ve blok ötelenmesi). Birbirinden bağımsız iki öteleme deposu olmazdı. Kapanış: blok ötelenmesi yalnızca zaten sayılmış molekül ötelenmelerinin kolektif bir kipidir; gerçekten yeni olan yalnız bloğun dönme serbestlik derecesidir. Serbestlik derecesi dengesinde (Böl. 4) ötelenme tam bir kez sayılır.
2. "Enerjik olarak yansız" (aksiyom 5) ile "bağlar enerji alır" (aksiyom 2/3b). Görünürde bir çelişki. Kapanış: yansızlık bir cisim olarak molekülü ilgilendirir; enerji moleküller arasında (BAĞ'da) ve kolektif kinetik kiplerde oturur, asla molekül içi bir tamponda değil. Molekül ≠ bağ — böylece çelişki çözülür.
3. Eksik çekim (aksiyom 4 yalnız itmeyi anar). Çekim olmasa yoğun faz olmazdı. Kapanış: çekimi BAĞ (ε derinlikli çukur), itmeyi SERT ÇEKİRDEK (σ'da sert kabuk) sağlar. Birlikte etkin bir çift potansiyeli.
4. Karşı terimsiz geri besleme (aksiyom 6). "Eksik komşu arttıkça daha çok kopma" tek başına olumlu bir geri beslemedir (tam çözülmeye dek çatlak ilerlemesi). Kapanış: onarım terimi (daha çok komşu → daha çok onarım) artı sonlu bağ enerjisi ε'nın (BAĞ) Boltzmann çarpanı ÇATLAK'ı dengeler ve bir kaçıp gitme yerine tanımlı bir geçiş üretir.
5. Dile getirilmemiş zaman ölçeği ayrımı (aksiyom 7). "Hızlı ve düşük aktivasyonlu", ancak ondan bağ ile öteleme dinamiğinin adyabatik ayrımı çıkarıldığında kullanışlı olur. Kapanış: Z = Z_trans·Z_vib·Z_rot·Z_bond çarpanlarına ayrılışının önkoşulu olarak açıkça saptanmıştır.
6. Eksik geometri ve eksik hareket düzeneği. Aksiyomlar ne koordinasyon örgüsünü ne de akmayı neyin olanaklı kıldığını anar. Kapanış (belirtilmiş ek varsayım olarak): değişken bağ doluluklu sabit bir referans koordinasyon geometrisi; kanonik 2B durum olarak üçgen koordinasyon (z_max = 6). Devingenlik kaydırmalı bulmaca boşluk düzeneği üzerinden. Bu iki nokta yedi aksiyomun ötesindeki tek gerçek eklemedir ve Bölüm 6'da stratejik soru olarak işaretlenmiştir.
7. Gizli yitim boşluğu. Sürüm 07, bağsız moleküllerin yanından geçerken enerji yitirip yitirmediğini açık bırakıyordu. SERBEST UÇUŞ ile kapanış: mikroskobik sürtünme/elektrostatik yoktur; bütün enerji ve momentum akışları yalnız SERT ÇEKİRDEK çarpışmaları ve BAĞ kinetiği üzerinden gider. Böylece enerji dengesi (Böl. 4.3) kapanır ve Maxwell–Boltzmann istatistiği gaz sınır durumunda tam olur (esnek çarpışmalar arası serbest ideal gaz).
8. Düzen parametresinin enerjik tutarsızlığı. Sürüm 08, işbirlikçi terim −(K/2)x²'yi serbest enerji g(x)'te (x(τ)'yu üretmek için) kullanıyor, ama iç enerji U/N'de kullanmıyordu. Bu termodinamik olarak tutarsızdır: aynı moleküler alan terimi ikisinde de görünmelidir. Kapanış (Böl. 4.3): −(K/2)x² U/N'ye eklenir; gizli ısı bununla ölçülebilir biçimde artar (Böl. 5.3), katı ve gaz platolarını değiştirmeden.
9. s/g kalibrasyonu yalnız niteliksel. Sürüm 09 yoğunluk eksenini (van der Waals hâl denklemi) yalnızca indirgenmiş biçimde vermişti. Kapanış (Böl. 7.4): kritik nokta (T_c, p_c) üzerinden mutlak kaynama çizgisi hesaplanır ve T_b(p_atm) ölçüm değerleriyle karşılaştırılır; dizgesel sapma (Guldberg kuralı) adlandırılır.
10. Üçlü noktasız ayrı eksenler. Erime ve kaynama ekseni ayrı hesaplardı. Kapanış (Böl. 7.5): her ikisi bir (τ, v) düzleminde birleştirilir, üçlü nokta kurulur ve kritik nokta işaretlenir.

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

3 — Resimli betimlemeler

Özet. Bölüm 3, katı, sıvı ve gaz hâlleri ile geçişlerini sezgisel olarak (durağan çizgi-roman şeritleriyle) açıklar: üç tekil faz (3.1), yeniden kristalleşme/erime/kaynama/süblimleşme faz geçişleri ile üçlü ve kritik nokta (3.2), ses anomalisi (3.3), Brown hareketi ve keskinleştirilmiş misel kavramı (3.4–3.5) ve üç boyutlu sıvı penceresi (3.6). Betimlemelere dayanarak, farklı fazlardaki devingen moleküler ilişkileri zihinde canlandırmaya yardım etmesi amaçlanan şemalar (çizgi-roman şeritleri) hazırlanmıştır. Bu canlandırılabilirlik modelin asıl çekirdeğidir. Nitelikli ve bilimsel olarak değerli olması, istatistiksel mekanik üzerinden modelin niteliksel görgül eğilimleri/gözlemleri büyük matematik olmadan da öngörülebilir kıldığının akla yatkın biçimde gösterilebilmesiyle olur.

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

3.1 Üç tekil faz

Katı madde. Neredeyse bütün komşular moleküler bağlarla birbirine bağlıdır (BAĞ), moleküllerin koordinasyon sayılarının ortalaması z geometrik en büyük değere yakındır. Molekül ağı boşluksuz tek bir katı bloktur (SERT ÇEKİRDEK sıkı paketlenmiş) — sıvıdaki gibi bir iç kaydırmalı bulmaca sıkışık sürümde olanaksızdır. Tek hareket: denge konumları çevresinde titreşim (bağ titreşimi, ÜÇLÜ AKOR-b); bunu komşularla olan moleküler bağlar belirler. Momentum/enerji ağda esnek biçimde salınır (KORUNUM, SERBEST UÇUŞ: yitimsiz). Sınır durumda (çok düşük sıcaklık, ideal paketleme) ideal bir kristalin ideal hâli ortaya çıkar.

Sıvı. İşbirlikçi ÇATLAK ağı delik deşik etmiştir: bağlar yer yer eksiktir, moleküllerin koordinasyon sayılarının ortalaması z düşmüştür, katı blok içeriden birçok küçüğe ayrılmıştır. Monolitten bir moloz yığını olmuştur. Bu bloklar ötelenir ve döner (ÜÇLÜ AKOR-c). Aralarında boşluklar açılmıştır ve onlardan geçerek birbirinin yanından kayarlar — Brown hareketinin kaydırmalı bulmacası etkindir. Yerel koordinasyon geometrisi idealleştirilmiş biçimde korunur; değişen yalnızca bağların kaçının gerçekten var olduğu ve oluştuğudur.

Gaz. ÇATLAK utku kazanmıştır: bağlar hemen her yerde kopuktur (z → 0). Uç durumda hiçbir molekülün başka bir molekülle moleküler bağı kalmamıştır. Bloklar tek tek moleküllere ayrılmıştır. Serbestçe uçarlar, yalnız SERT ÇEKİRDEK çarpışmalarıyla eşleşerek (aralarında SERBEST UÇUŞ); hızlar Maxwell–Boltzmann'ı izler (KORUNUM). Bağ titreşimi ve blok dönmesi sönmüştür — saf öteleme (ÜÇLÜ AKOR-a). Sınır durumda (yüksek sıcaklık, düşük basınç) ideal bir gazın ideal hâli ortaya çıkar.

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

3.2 Faz geçişleri

Yeniden kristalleşme (katı(1)/katı(2)). Daha erimenin altında katı madde, sertliğini yitirmeden örgüsünü değiştirebilir: farklı koordinasyon sayılı iki ideal paketlemenin (YMK z_max = 12 ↔ HMK z_max = 8) serbest entalpileri μ(τ, p) farklıdır; μ_YMK = μ_HMK olduğunda koordinasyon geometrisi sıçramalı olarak değişir (erimedeki gibi yalnız doluluk z̄ değil). Görünür olan küçük bir c_v basamağıdır (Δc_v = z_max/4 değişimi = 1,0 k_B) artı zayıf bir gizli dönüşüm ısısı. YZ savı: yalnız 3B'de düşünülebilir (kesikli, nicemlenmiş bağlantılılık; Böl. 6.3). İki boyutlu durumda gerçek yeniden kristalleşmeyi, yani iki katı faz arasındaki geçişi canlandırmak daha kolaydır. Düzlemde deliksiz yalnız üç düzenli döşeme vardır. Bir düzlemi eşkenar üçgenlerle, karelerle ya da düzgün altıgenlerle deliksiz olarak tümüyle kaplamak (döşemek) olanaklıdır. Bu örgülerde kenarlar moleküler bağları, köşeler molekülleri temsil ediyorsa, katı(üçgen)→katı(kare) bir faz geçişi düşünülebilir. Bu durumda üçgen döşemedeki molekülün komşu moleküllere altı bağı olurken, kare örgüdeki molekülün yalnız dört bağı olur.

Erime (k/s). Sıcaklık yükselirken ağ yeterince çok yerde yırtılır ve katı bütün blok sertliğini yitirir (işbirlikçi ÇATLAK Maxwell rijitlik eşiğinin z̄ = 2d altına iner — ağ gevşek olur). Termodinamik miseller, yani daha büyük komşu bölgelerle artık bağlarla bağlı olmayan bölgeler oluşur. Monolit moleküler bakışta moloz yığınına dönüşür ve sınırlı ölçüde akabilir. İlk boşluklar molekül bloklarının kaydırmalı bulmacadaki gibi yer değiştirmesine izin verir. Komşuluğa eksik bağlar misellerde iç dönmeye izin verir. Misel dönmesi (ÜÇLÜ AKOR-c) devreye girer. Birçok bağın kopması (BAĞ) gizli ısı yutar.

Kaynama (s/g). Daha fazla ısıtma boşlukların kendilerini perkole ettirir: sürekli bir boş uzaya birleşirler, bu uzayda bloklar çözülür ve tek tek moleküller uçup gider. Kaynama, "maddedeki boşluklar"dan "boşluktaki madde"ye geçiştir.

Süblimleşme (k/g). Düşük basınçta sıvı ara bölge yoktur: ağ, misellerdeki gibi bir blok oluşumu olmadan tek tek molekülleri doğrudan boş uzaya verir. Burada geometrik yönlerin (uzun moleküller, birbirine dolanan moleküller), elektrostatik yönlerin (manyetik ya da dipol momentleriyle yönelim), farklı moleküler bağ türlerinin (van der Waals, mezomeri dizgeleri arası etkileşimler, hidrojen köprüleri) termodinamik misellerin oluşup oluşmaması üzerinde nasıl bir etkisi olabileceği düşünülebilir.

Üçlü nokta. Tam orada iki eşik çakışır: Maxwell rijitliği (katı ↔ sıvı, sertliğin yitimi) ve boşluk perkolasyonu (sıvı ↔ gaz) aynı (τ, p)'de gerçekleşir. Üç resim de (katı blok [monolit], kayan bloklar [moloz yığını], serbest moleküller [toz bulutu]) bir arada bulunur.

Kritik nokta. Basınç yükseldikçe sıvı→gaz sıçraması gitgide küçülür; boşluk oranı gitgide daha sürekli değişir. Kritik noktada fark tümüyle yiter — "az boşluklu yoğun bloklar" ile "az moleküllü çok boşluk" sürekli olarak birbirine geçer. ("Sıvının bakışından bardak yarı boş, gazın bakışından yarı doludur.")

Yeniden kristalleşme 2B resimde. Çizgi-roman yeniden kristalleşmeyi artık 2B durum olarak alt alta gösterir: altta sıkı üçgenlenmiş örgü (koordinasyon z = 6, en düşük enerjili 2B ideal), üstte gevşek kare örgü (z = 4). Üçgen → kare geçişi hacmi büyütür (basınca karşı iş) ve bağ sayısını düşürür. Koordinasyon sayısının değiştiği, ama bir kristal örgünün sertliğinin sürdüğü görülür.

Enerji geçişlerinin ve ısı sığalarının niteliksel öngörüsü. Yalnız bağlantılılık z'den her iki büyüklük de kestirilebilir: (i) parçacık başına bağ enerjisi U/N ≈ −(ε/2)·z'dir (her bağ yarı sayılır), yani z ne kadar yüksekse enerji o kadar düşüktür; (ii) örgü değişimindeki c_v sıçraması model kuralı Δc_v = z_max/4 değişimini izler. Üçgen (z = 6) ↔ kare (z = 4) 2B durumu için bu, parçacık başına ΔU ≈ ε enerji farkı ve Δc_v = (6 − 4)/4 = 0,5 k_B c_v basamağı verir — 3B değerin (1,0 k_B) yarısı kadar. Belletici: daha sıkı = daha düşük enerjili, ama geçişte daha pahalı (daha büyük gizli ısı, daha büyük c_v basamağı), çünkü genişlerken daha çok bağ ve daha çok donmuş serbestlik derecesi serbest kalır.

Örgü yapısı yarı kararlılığı belirler mi? Kısmen. Örgü yapısı enerji sıralamasını saptar — hangi biçimin taban durum (üçgen, en düşük μ) ve hangisinin uyarılmış (kare) olduğunu. Gevşek biçimin düşük sıcaklıkta yarı kararlı kalıp kalmayacağına ise tek başına o karar vermez: sıkı ideal modelde (engelsiz) madde her (τ, p)'de hemen en düşük μ'lü fazı alır — orada yarı kararlılık yoktur. O ancak yeniden yapılanmayı kinetik olarak engelleyen sonlu bir etkinleşme enerjisiyle (çekirdeklenme engeli) doğar. Örgü geometrisi böylece termodinamik sıra düzenini, engel ise kinetik kalıcılığı sağlar — aşırı doymuş fazı da açıklayan aynı bileşen (bkz. aşağıda; niceliksel Böl. 6.7).

Aşırı doymuş faz — etkinleşme enerjisi. Ek çizim (3.2, sıra 2) bir sıvının üstünde aşırı doymuş gaz fazını gösterir: tek moleküllerin yanında bütün olarak uçan, 2–3 birleşmiş parçacıktan oluşan küçük öbekler (gri çubuklar). Model aşırı doymayı moleküler bağlar için gözle görülür bir etkinleşme enerjisiyle açıklar: küçük öbeklerin (2–3) kenar/yüzey oranı elverişsizdir, kenar bağları eksiktir — ancak kritik bir çekirdek kendiliğinden büyümeyi sürdürür. Kritik üstü bir çekirdek oluşmadığı sürece, μ_gaz > μ_sıvı yoğuşmaya izin verse bile gaz aşırı doymuş kalır. Tam bu engel ideal denge resminde yoktur ve burada bağımsız, "ideal resmin dışında duran" bir katkı olarak adlandırılır (niceliksel kestirim Böl. 6.7'de).

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

3.3 Sıvılarda ses hızı neden T ile düşer

Ses esnek bir dalgadır: v = √(K_eff/ρ). K_eff yalnız sağlam bağlardan doğar (BAĞ) — yaylar onlardır. T yükselirse işbirlikçi biçimde gitgide daha çok bağ kopar (ÇATLAK); yay ağı seyrelir ve K_eff yoğunluk ρ'dan daha dik düşer. Kopma paydayı değil, payı egemenler — bu yüzden v düşer. Moleküler yorum: sesin moleküler bağlar üzerinden iletimi, mekanik çarpışmalarla mekanik iletiminden hızlıdır. Sıcaklık artarken bir sıvıda ne kadar çok çatlak ortaya çıkarsa, ses o kadar sık moleküler bağlar yerine mekanik çarpışmalarla taşınmak zorunda kalır. Buna uygun olarak sesin hızı sıcaklık yükseldikçe yavaşlar. (Suyun anomalisini — ~74 °C'de bir ses en büyük değeri — Böl. 6.4 ele alır.)

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

3.4 Brown hareketi nasıl doğar — keskinleştirilmiş misel kavramıyla

İstatistiksel-termodinamik anlamda misel burada bir yüzey etkin madde kümesi değil, ÇATLAK ile çevrelenmiş bir bloktur: eksik bağlar (blok sınırındaki bağ-olmayanlar), kapatılmış küçük kristalin (sağlam bağlı çekirdeğin) çevresinde "kabuğu" oluşturur. Bu çekirdek bir iç işbirlikçi dönme hareketi yapabilir (ÜÇLÜ AKOR-c) — blok bütün olarak döner, sınır kabuğu ise onu komşulardan ayırır.

Brown hareketi şöyle doğar: moleküler bağların değiş tokuşuyla (ÇATLAK, ÇEVİK) tek tek moleküller bir miselin kabuğunu terk edip başka bir misele geçebilir. Böyle her geçiş momentum aktarır (KORUNUM; SERBEST UÇUŞ nedeniyle yitimsiz) ve ilgili iki miselin biçimini ve açısal momentumunu sıçramalı olarak değiştirir. Büyük bir sınama parçacığını (toz tanesi) bu tür birçok ilintisiz misel yeniden yapılanması yalar; momentum darbelerinin toplamı ortalamada sıfırdır, ama sonlu değişkelidir → rastgele yürüyüş. T ne kadar yüksekse ÇATLAK olayları ve misel değişimleri o kadar sık, titreşim o kadar güçlüdür. Brown hareketi böylece mikroskobik misel değiş tokuş sürecinin makroskobik gölgesi ya da sonucudur.

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

3.5 Misel sınır katmanlarında tepkimeler neden hızlanır

Utandırıcı ön not. Aşağıdaki uslamlama resimden bir çıkarımdır ve bildiğim kadarıyla şimdiye dek görgül olarak kanıtlanmamış ya da incelenmemiştir; itiraf etmeliyim ki bunu iyice araştırma zahmetine girmedim. Yeniliğe ilişkin tek dayanak, YZ dizgelerinin bu sonucu anmamasıdır. Bu sav, resimli tasavvurların nasıl yeni soru koyuşlarına götürebileceğinin bir örneğidir.

Orta derecede çözen çözücülerde tepkenler yeğlemeli olarak misellerin sınır kabuklarına yerleşir — modelde tam da çatlak/blok sınırlarına (ÇATLAK). İki etki tepkimeyi hızlandırır: (i) boyut indirgemesi — tepkenler birbirini artık tüm hacimde değil, yarı iki boyutlu bir sınır yüzeyi boyunca arar; iki boyutlu bir arama istatistiksel olarak üç boyutludan daha hızlı bulur. (ii) Derişikleşme — sınır kabuğu tepkenleri toplar, yerel yoğunlukları artar. Resimli olarak çatlak sınırları, eşleri birbirine getiren kılavuz raylar gibi etkir.

Boyut etkisinin sıkı bir çekirdeği vardır: rastgele yürüyüş (Pólya). İki boyutlu bir aramanın neden istatistiksel olarak üç boyutludan daha iyi bulduğu yalnız sezgiden değil, olasılık kuramının bir teoreminden çıkar. George Pólya (1887–1985) 1921'de bir örgü üzerindeki rastgele yürüyüş için şunu gösterdi: bir ve iki boyutta yürüyen, 1 olasılıkla (yani kesinlikle) her noktaya sonsuz kez döner; düzlem örgü ℤ²'de de doğru ℤ'de olduğu gibi. Üç boyuttan başlayarak bu artık geçerli değildir: dönüş olasılığı gerçekten 1'den küçüktür (yalın kübik örgü ℤ³'te ancak yaklaşık 0,34). Sezgisel olarak, Pólya'nın resminde: "Sarhoş bir insan er geç kesinlikle evini bulur, sarhoş bir balık ise bulmaz." Başarılı bir çarpışma için birçok kez buluşması gereken iki tepkene uygulandığında bu şu demektir: serbest üç boyutlu yürüyüşlü, çok iyi çözen bir çözücüde bir kez ayrılmış eşler sonlu olasılıkla sonsuza dek uzaklaşır — yinelenen tepkime denemeleri başarısız olur ve verim, salt bu geometrik-istatistiksel nedenden düşük kalır. Buna karşılık (yarı) iki boyutlu bir misel sınır katmanında yeniden karşılaşma kesindir (olasılık 1): eşler tepkime başarılana dek hep yeniden bir araya gelir. Sınır katmanı böylece yalnız derişikleştirmez, aynı zamanda yeniden karşılaşmalar için daha elverişli boyutu sağlar. (Bilinçli olarak keskinleştirilmiş) sav şudur: çok iyi çözen çözücülerdeki verimler rastgele yürüyüş nedenleriyle eğilim olarak daha kötüdür.

Kaynak: G. Pólya (1921), Über eine Aufgabe der Wahrscheinlichkeitsrechnung betreffend die Irrfahrt im Straßennetz, Math. Ann. 84, 149–160 — burada M. Eisermann, Höhere Mathematik 3 (Univ. Stuttgart), alıştırma W-2×2'ye göre anılmıştır. (Özgün kaynak 5.9.2026'da açıldı).

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

3.6 Üç boyutlu sıvı penceresi

3.1–3.5'teki çizimlerden farklı olarak bu çizim bir molekül sahnesi değil, Böl. 7.9/7.10'daki hesabın sonucunu bir τ ekseni üzerinde şema olarak gösterir (bağ enerjisi ε birimlerinde). Amaç, tek bir 2B örgünün sıvı faz için neden neredeyse hiç alan bırakmadığını, 3B örgüde ise bunun sorunsuz olanaklı olduğunu görünür kılmaktı; her iki durumda da aynı rijitlik kuralı (z̄_c = 2d) geçerlidir.

Resim ① — 2B (pencere yok). τ ekseninde τ_m = 1,20 ile τ_triple = 1,16 hemen hemen üst üste düşer (τ_triple, τ_m'nin altında), τ_c = 1,50 hemen üstünde: sıvı faz alanı bir yarığa sıkışmıştır — katı neredeyse doğrudan gaza geçer (süblimleşmeye yakın).

Resim ② — 3B (temiz pencere). Aynı eksen, ama gerilmiş: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00. τ_triple ile τ_c arasında bağımsız sıvı alan yatar (aşağıya geniş, τ_triple − τ_m = 1,63; yukarıya dar, τ_c − τ_triple = 0,18). Dar üst kenar böyle olarak işaretlenmiştir.

Resim ③ — nedeni. z̄_c = 2d boyutla büyür (4 → 6), böylece v_solid düşer (indirgenmiş hacimlerde 0,50 → 0,667, yani katı madde sıvıya göre "daha gevşek" olur) ve τ_c = z_max/4 ikiye katlanır (1,50 → 3,00). Kritik nokta erime noktasının çok üstüne çıkar. Sıvı faz için yer ancak böyle doğar. Bu saf bir boyut etkisidir.

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

4 — İstatistiksel-mekanik model (sürüm 09'dan değişmeden)

Özet

Bölüm 4, modelin istatistiksel-mekanik iskeletini 2B'de kurar (üçgen örgü, z_max = 6):

• 4.1 Serbestlik dereceleri ve değişkenler: düzen büyüklüğü, ortalama koordinasyon sayısı z̄ = z_max·x ile birlikte bağ oranı x ∈ [0,1]'dir; ≥ 2 blok başına bir dönme serbestlik derecesi eklenir, sıcaklık τ bağ enerjisi ε birimlerinde sayılır.
• 4.2 İşbirlikçi iki-durumlu bağ: kendiyle tutarlı bağ denklemi x = 1/(1 + G·exp[−(ε + K·x)/τ]) (K = J·z_max, G = z_max), Bragg–Williams serbest enerjisinin durağanlık koşulu olarak çıkar; onun genel en küçüğü fiziksel dalı verir, bunun sürekli mi yoksa sıçramalı mı (birinci mertebeden) gittiğine Böl. 5.5 karar verir.
• 4.3 İşbirlikçi terimli enerji dengesi: U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x² (2B), ötelenme/bağ titreşimi/blok dönmesini, bağ konfigürasyon enerjisini ve yeni işbirlikçi moleküler alan enerjisi −(K/2)x²'yi birleştirir. İlgili ısı sığası c_v, erime tepesini keskinleştiren, ama katı platoyu ve gaz dalını dokunulmadan bırakan pozitif ek terim −K·x·ẋ'i alır.

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

4.1 Serbestlik dereceleri ve değişkenler

• Sabit bir referans koordinasyon geometrisi üzerinde N molekül. Kanonik 2B durum: üçgen örgü, z_max = 6, blok başına bir dönme SD (ÇOKLUK, BAĞ). Bu bölüm açıkça 2B yorum için geçerlidir; 3B genelleme Böl. 6'dadır.
• Her geometrik komşu çifti için bağ değişkeni s_ij ∈ {0, 1} (1 = sağlam).
• Koordinasyon sayısı z_i = Σ_j s_ij; ortalama z̄ = 2 N_bond / N — düzen parametresi. x = p_b ∈ [0, 1] bağ oranıyla z̄ = z_max · x'tir.
• Enerji birimi τ ≡ k_B T (klasik, karesel serbestlik derecesi başına ½τ). Bütün enerjiler bağ enerjisi ε birimlerinde.

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

4.2 İşbirlikçi iki-durumlu bağ (Ising benzeri, parametre J)

Yalıtılmış bir bağın (enerji sağlamken −ε, kopukken 0) Maxwell–Boltzmann dengesinde p_b⁰(τ) = 1/(1 + e^(−ε/τ)) olur. İşbirlikçi ÇATLAK (aksiyom 6), etkin bağ kararlılığını bağ oranı x ile orantılı bir moleküler alan kadar kaydırır. Moleküler alan (Ising) yaklaşımında kendiyle tutarlı bağ denklemi şu olur:

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

Burada:

• K = J·z_max işbirlikçi eşleşme gücüdür (ÇATLAK geri beslemesi: daha çok komşu bağ → daha yüksek etkin bağ kararlılığı). J > 0 geçişi keskinleştirir; kritik bir J'nin üstünde z̄(τ) süreksiz olur — bu, erimenin gizli ısısıdır. Saf iki-durumlu durum, J = 0 sınır durumudur (moleküler alan yok, sürekli geçiş).
• G > 1 kopuk durumun dejenereliğidir: serbest kalan bir molekül faz uzayı (öteleme özgürlüğü) kazanır. G yüksek τ'da bağ kopmasını sürer ve x'in ½'de takılıp kalmak yerine 1/(1+G)'ye gitmesini sağlar (gaz sınır durumu). Kanonik olarak G = z_max.

Bu denklem, molekül başına bir Bragg–Williams serbest enerjisinin durağanlık koşulu olarak çıkar (n_slots = z_max/2 bağ yeri):

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

−(K/2)x² terimi işbirlikçi (ferromıknatıssal) eşleşme, son terim kopma dejenereliğidir. g'nin x üzerindeki genel en küçüğü fiziksel dalı verir; iki en küçüğün eşit derinlikte olduğu yerde birinci mertebeden geçiş yatar. ÇEVİK (aksiyom 7), x(τ)'nun her sıcaklıkta denge değerini almasını haklı çıkarır.

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

4.3 Enerji dengesi — işbirlikçi terimle

Moleküler alan terimi −(K/2)x² gerçek bir enerjidir (τ'dan bağımsız) ve U/N'de görünmelidir. Tutarlı denge böylece şudur (2B):

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

• τ·[…]: ötelenme (1) + bağ titreşimi (z̄/4) + blok dönmesi (½ n_blk) — ÜÇLÜ AKOR.
• −(ε/2)·z̄: bağ konfigürasyon enerjisi — BAĞ.
• −(K/2)·x²: işbirlikçi moleküler alan enerjisi — ÇATLAK. U/N'de yeni.

İlgili ısı sığası c_v = d(U/N)/dτ (k_B birimlerinde) ek bir katkı alır:

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

burada yeni terim −K·x·ẋ pozitiftir (ẋ < 0) ve tam olarak birçok bağın aynı anda koptuğu yerde etkir. Erime tepesini keskinleştirir, katı platoyu (ẋ ≈ 0) ya da gaz dalını değiştirmeden.

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

5 — Durum toplamları ve ısı sığaları

Özet

Bölüm 5, durum toplamlarını ve ısı sığalarını hesaplar ve kalibrasyon için model verilerini sağlar (Böl. 7):

• 5.1 Çarpanlara ayrılmış durum toplamı: zaman ölçeği ayrımı sayesinde Z = Z_trans·Z_vib·Z_rot·Z_bond biçiminde çarpanlarına ayrılır.
• 5.2 Faz başına ısı sığası: c_v = d(U/N)/dτ; pozitif terim −(ε/2)(dz̄/dτ) gizli ısının karşılığı olan Schottky benzeri tepeyi üretir.
• 5.3 Sayısal c_v(τ) eğrisi: c_v tepesi 5,51, τ_m = 1,20'de, katı plato 2,50, gaz dalı → 1; gizli ısı L ≈ 1,41 ε/molekül (işbirlikçi terimle) 0,48 ε'ya (terimsiz) göre üç katına çıkmıştır.
• 5.4 / 5.4a İki ayrı eşik: bağlantılılık perkolasyonu z̄_c = z_max·p_c ≈ 2,08 (τ_b ≈ 1,85'te) öbek çözülmesini (kaynamaya yakın) işaretler ve bir sonuçtur; bunu genel Maxwell rijitliğinden z̄_c = 2d (erime yoğunluğu, sertlik yitimi) temiz biçimde ayırmak gerekir.
• 5.5 Erime geçişinin mertebesi: ölçüt ln G > 2 + 4ε/K; 2B (G = 6) her J için süreklidir, 3B (G = 12) K > ~8,2 olunca birinci mertebedendir (gösterim: Δz̄ ≈ 5,3, ΔU ≈ 1,75 ε).
• 5.6 Sonlu boyut incelemesi: p*(L) aşağıdan tekdüze biçimde p_c = 0,3473'e yakınsar; z̄_c ≈ 2,08 sonlu boyuta karşı gürbüzdür.
• 5.7 Tam 3B c_v taraması: K > K_kritik ≈ 8,25 (J > 0,68) olunca gerçek sıçrama Δz̄ ile gizli ısı ΔU açılır; erime tepesi J ile daha yüksek τ'lara göçer (sürekli 2B tepeye göre +0,5 … +2,2'lik pozitif kayma).

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

5.1 Çarpanlara ayrılmış durum toplamı

Zaman ölçeği ayrımı (ÇEVİK) nedeniyle kanonik durum toplamı çarpanlarına ayrılır:

Z = Z_trans · Z_vib · Z_rot · Z_bond.

• Z_trans = (1/N!)·(A/λ_T²)^N, λ_T = h/√(2π m τ) (klasik ideal 2B öteleme; ÇOKLUK, KORUNUM, SERBEST UÇUŞ).
• Z_vib = (τ/ħω₀)^(N_bond) (her sağlam bağ için klasik harmonik salınıcı; ÜÇLÜ AKOR-b).
• Z_rot = ∏_(bloklar≥2) (2π I_b τ)^(1/2)/ħ (her blok için klasik 2B katı döner; ÜÇLÜ AKOR-c).
• Z_bond = Σ_{s} exp[(ε/τ)·Σ s_ij + (K/τ)·işbirlikçi terim] (Ising benzeri bağ toplamı; BAĞ, ÇATLAK). Moleküler alan değerlendirmesi 4.2'deki x(τ), z̄(τ)'yu verir.

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

5.2 Faz başına ısı sığası (analitik)

U/N'den (Böl. 4.3), c_v = d(U/N)/dτ (k_B birimlerinde):

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

−(ε/2)·(dz̄/dτ) terimi pozitiftir (dz̄/dτ < 0: ısıtmak bağları kırar) ve birçok bağın aynı anda kırıldığı yerde egemendir — gizli ısının karşılığı olan Schottky benzeri tepe.

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

5.3 İşbirlikçi terim içeren sayısal c_v(τ) eğrisi — gizli ısı niceliksel

Kendiyle tutarlı bağ (genel en küçük), n_blk(τ) bağ perkolasyonu Monte Carlo'dan (Union-Find, L = 100, 5 gerçekleme), c_v = dU/dτ (merkezî farklar), U/N −(K/2)x² dâhil (cv_scan_09.py, ε = 1, J = 0,55, G = z_max = 6). Seçilmiş dayanak noktaları:

Sonuç (sayısal). Üç faz sıralaması doğrulanmış kalır ve artık daha keskindir:

c_v(katı) ≈ 2,50 < c_v(sıvı, en büyük) ≈ 5,51 > c_v(gaz dalı) ↓ → 1.

Gizli ısı (taban çizgisi üstündeki tepe alanı). (c_v − taban çizgisi) erime penceresi τ ∈ [0,70 ; 1,70] üzerinde integre edilirse (taban çizgisi = c_v(0,70) ile c_v(1,70) üzerinden geçen doğru), şu çıkar:

L ≈ 1,41 ε/molekül (işbirlikçi terimle), L ≈ 0,48 ε/molekül (işbirlikçi terimsiz, eski denge), ΔL ≈ +0,92 ε/molekül.

İşbirlikçi terim tepeyi yaklaşık 1,7 kat yükseltir (3,25 → 5,51) ve gizli ısıyı üçe katlar. Fiziksel olarak: birçok komşu bağın aynı anda kırılması (ÇATLAK) yalnız konfigürasyon enerjisini (−ε/2·z̄) değil, moleküler alan enerjisini de (−K/2·x²) depolar — ikisi de erimede serbest kalır.

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

5.4 Sonuç olarak bağlantılılık perkolasyon eşiği

Dev öbek oranı P_∞(τ) τ ≈ 1,85'te çöker; orada x ≈ 0,347 = p_c(üçgen) = 2·sin(π/18) = 0,3473, yani z̄_c = z_max·p_c ≈ 2,08 — konulmuş değil, sonuç. Bu, bağlantılılık eşiğidir (sürekli öbek) ve kaynamaya yakın (τ_b) öbek çözülmesini işaretler; Maxwell rijitliğiyle (z̄_c = 2d, erime) karıştırılmamalıdır — bkz. Böl. 5.4a'daki temiz ayrım.

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

5.4.a İki farklı eşik temiz biçimde ayrıldı

Şimdiye dek bulanık olan "perkolasyon/rijitlik" ikili adlandırması artık çözülmüştür. Söz konusu olan, karıştırılmaması gereken fiziksel olarak farklı iki eşiktir:

Gerekçe. Erime, yoğun bir bloğun mekanik sertliğinin yitimidir, bağlılığın çözülmesi değil. Bu yüzden ilgili büyüklük Maxwell kısıt sayımıdır (ortalama alan / genel rijitlik): merkezî kuvvetlerden oluşan D boyutlu bir ağ tam olarak ortalama koordinasyon z̄ = 2D'ye ulaştığında izostatik olur (2B → 4, 3B → 6). Bağ bağlantılılığı (z̄_c = 2,08) bunun belirgin biçimde altındadır; z̄ = 2,08'de daha yeni bağlı olan bir ağ çoktan gevşektir (yumuşak kipler), yani zaten sıvıdır. Bu nedenle bağlantılılık kaynamaya yakın daha geç gelen öbek çözülmesini işaretlerken, rijitlik daha erken katılaşmayı/sıvılaşmayı (erime yoğunluğunu) belirler. Bu ayrım Böl. 7.8'deki üçlü noktayı parametresiz ve gerçekçi biçimde sabitler; 5.4/5.6'da doğrulanan 2,08'lik bağlantılılık eşiğini değiştirmez, yalnız onun yorumunu (τ_b = kaynamaya yakın, erime değil).

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

5.5 Erime geçişinin mertebesi

Analitik ölçüt. Spinodal bölge (g'nin iki en küçüğü) τ < K/4 için vardır; geçiş yaklaşık τ_t ≈ (ε + K/2)/ln G'de yatar. Gerçek bir sıçrama (birinci mertebeden geçiş) τ_t < K/4 gerektirir, yani K → ∞ sınırında:

ln G > 2 + 4ε/K.

Bundan kırılma dejenereliği G = z_max hakkında dikkate değer bir sonuç çıkar:

Sayısal doğrulama (τ adımı başına en büyük |Δx|, detect_first_order):

Yorum. Kanonik 2B modelde (G = 6) erime ilkece süreklidir ("yumuşak" bir geçiş/kroskover) — hiçbir J için gerçek bir sıçrama doğmaz, çünkü ln 6 < 2. Ancak 3B paketlemenin daha yüksek koordinasyonu/dejenereliği (G = 12), yeterli işbirliğinden (K > 8,2) itibaren sonlu gizli ısılı birinci mertebeden bir geçişe izin verir. Bu gerçeklikle örtüşür: 2B erime çoğu kez süreklidir (KTHNY benzeri), 3B erime birinci mertebedendir. J = 0 (moleküler alan yok) her zaman sürekli referans durumdur. — Fiziksel olarak J, işbirlikçi enerjinin bağ enerjisine oranını (K/ε = J·z_max/ε) gösterir: sağlam bir komşu bağın kendi bağını ne kadar kararlı kıldığını.

Birinci mertebe gösterimi (3B esinli, cv_scan_09.py, J = 1,0, G = 12, z_max = 12): z̄, τ ≈ 2,82'de Δz̄ ≈ 5,30 kadar sıçrar; gizli ısı enerji sıçraması olarak ΔU ≈ 1,75 ε/molekül. Aynı K bir G = 6 örgüsünde (2B) sürekli kalır — mertebe niteliği yalnız G = z_max'a bağlıdır.

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

5.6 Perkolasyon eşiğinin sonlu boyut incelemesi

L × L üçgen örgülerde (periyodik) bağ perkolasyonu, dev öbek oranı P_∞(p) L üzerinden ölçeklendi (cv_scan_09.py; L = 50/100/200/400, 8/5/3/2 gerçeklemeyle). Kesin eşik p_c = 2·sin(π/18) = 0,34730 → z̄_c = z_max·p_c = 2,0838.

P_∞(p) (ortalama ± yayılım):

Eşik kestirimi p*(L) (P_∞ 0,5'i keser):

Sonuç. p*(L) tekdüze biçimde aşağıdan p_c = 0,3473'e yakınsar; p_c'deki yayılım ±0,103'ten (L = 50) ±0,017'ye (L = 400) düşer — geçiş L ile keskinleşir, ama eşiği kaydırmaz (evrensel perkolasyon değeri). Böylece z̄_c ≈ 2,08 = z_max·p_c sonlu boyuta gürbüz bir model değeri olarak doğrulanmıştır.

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

5.7 Birinci mertebeden erimenin tam 3B c_v(τ) taraması

Böl. 5.5'teki tek birinci mertebe gösterim koşusu (J = 1,0, G = 12) şimdi işbirlikçilik J üzerinde tam bir c_v(τ) taramasına genişletilmiştir (firstorder_3d_16.py, n_blk için YMK bağlantılılığı, 3B enerji dengesi U/N = τ·[3/2 + z̄/4 + (3/2)n_blk] − (ε/2)z̄ − (K/2)x²). Böl. 5.5'teki analitik mertebe ölçütü (birinci mertebeden ⇔ K > K_kritik = 4ε/(ln 12 − 2) = 8,25) durumları ayırır; sürekli 2B referans tepesi τ_2D = 1,20'de yatar (c_v ≈ 5,5).

*sürekli durum: Δz̄/ΔU dik kroskoverdir (gerçek sıçrama değil); τ_geçiş orada c_v tepesidir.

Sonuç. (i) Sürekli 3B durum J = 0,55, τ_geçiş = τ_m ≈ 1,72 verir — bağımsız bağ kinetiğiyle (bondkinetics_3d_15.py: τ_m(YMK) = 1,70) yapılan çapraz sınama tutar. (ii) K > K_kritik ≈ 8,25 (J > 0,68) olunca sonlu gizli ısı ΔU ile (süreksizlik üzerinden enerji sıçraması) gerçek sıçrama Δz̄ açılır — erime birinci mertebeden olur, tam da ln 12 > 2 ölçütünün istediği gibi. (iii) J = 1,0 çapraz sınaması Böl. 5.5'teki gösterim değerini Δz̄ ≈ 5,4 yeniden üretir. (iv) Sıçrama yüksekliği Δz̄ ve τ_geçiş konumu J ile büyür: geçiş yukarı göçer — 3B paketlemenin daha yüksek koordinasyonu ve daha güçlü işbirlikçilik erime sıcaklığını yükseltir. Sürekli 2B tepeye göre (τ_2D = 1,20) böylece baştan sona pozitif bir tepe kayması (+0,5 … +2,2) çıkar ve bu K ile tekdüze artar. 3B erimenin birinci mertebe niteliği böylece artık yalnız bir noktada değil, tüm işbirlikçilik aralığında belgelenmiştir.

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

6 — 3B yorumun matematiksel akla yatkınlığı

Özet

Bölüm 6, bağ resminin 2B'den 3B'ye genelleştirilmesini salt matematiksel olarak akla yatkın kılar (kendi HTML görselleştirmesi olmaksızın):

• 6.1 Dönme serbestlik derecesi 2B → 3B: blok dönmesi bir eksenden üç eksene, öteleme iki yönden üç yöne, z_max 6'dan 12'ye büyür; c_v temel değerleri yükselir (c_v,katı 2,5 → 4,5, c_v,gaz 1 → 3/2), c_v(τ) eğrisinin biçimi kalır.
• 6.2 Ayırt edici özellik olarak bağlantılılık: z_max nicemlenmiş bir büyüklüktür ve en sıkı paketlemeyi (YMK/HSP, z_max = 12) HMK'den (z_max = 8) ara durum olmaksızın sıçramalı olarak ayırır.
• 6.3 Ek geçiş katı(1)→katı(2): kesiklilikten YMK↔HMK polimorfizmi ve parametresiz 1,0 k_B'lik c_v basamağı çıkar; sıkı ideal resimde (θ_HMK = θ_YMK) bu oluşmaz, ancak daha yumuşak bir HMK (0,61 açılma eşiğinin altında θ_HMK < θ_YMK) τ_ss < τ_m'yi açar.
• 6.4 / 6.5 Düşünceler: sıvı suyun ses anomalisi iki örgü tipinin bulanıklaşmış yakın düzen değişimi olarak, ve misel dönmesine fraktal girdap alanı seçeneği (salt niteliksel).
• 6.6 / 6.8 / 6.9 Mikroskobik θ(örgü) üç basamakta: Einstein yay modeli (θ_HMK/θ_YMK ≈ 0,82…0,94), Debye benzeri fonon DOS (≈ 0,87…1,00) ve Lennard-Jones tam örgüsü (≈ 1,0; ideal HMK kaymaya kararsız, YMK kararlı taban durum) — hepsi 0,61 eşiğinin üstünde, demek ki gerçek polimorfizm harmonik ideal resmin dışındadır.
• 6.7 Yarı kararlılık (2B): üçgen ↔ kare yeniden yapılanması için büyüklük mertebesi hesabı (ΔU = +1 ε, ΔA ≈ +15,5 %) ile kinetik ek olarak çekirdeklenme engeli ΔG\* = πλ²/Δμ_A.

Bu bölüm 3B genelleştirmeyi salt matematiksel olarak akla yatkın kılar. 3B hiçbir zaman bir HTML sürümünde görselleştirilmez. Amaç yalnızca dönme serbestlik derecesinin ve bağlantılılığın 2B'den 3B'ye nasıl değiştiğini ve bununla hangi yeni geçişlerin düşünülebilir olduğunu göstermektir.

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

6.1 Dönme serbestlik derecesinin 2B → 3B değişimi

2B'ye göre niteliksel olarak tek yeni nokta blokların dönmesidir:

• 2B: katı bir bloğun 1 dönme SD'si vardır (düzlemde dönme) → U_rot/N = ½ n_blk τ.
• 3B: katı bir bloğun 3 dönme SD'si vardır (üç Euler ekseni) → U_rot/N = (3/2) n_blk τ.

Aynı zamanda öteleme 2'den 3 SD'ye (U_trans/N = (3/2)τ) ve en büyük koordinasyon sayısı z_max = 6'dan (üçgen) z_max = 12'ye (en sıkı küre paketlemesi) büyür. Böylece 3B'deki enerji dengesi şudur:

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

ve ısı sığasının temel değerleri tutarlı biçimde kayar:

Eğrinin yapısı (katı plato, −(ε/2)dz̄/dτ nedeniyle sıvı en büyüğü, gaz değerine düşüş) aynı kalır; yalnız düzeyler boyutla ve z_max ile ölçeklenir.

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

6.2 Farklı kristal örgülerin etkisi — ayırt edici özellik olarak bağlantılılık

3B'de en az iki ideal örgü ayırt edilmelidir ve bunlar bağlantılılıklarıyla (koordinasyon sayısı z_max) ayrılır:

Belirleyici ayırt edici özellik bağlantılılık z_max'tır: hem katının ısı sığası düzeyini (z_max/4) hem de perkolasyon eşiği z̄_c'yi belirler. Her örgü yalnız bir sabit koordinasyon sayısına izin verdiği için bağlantılılık nicemlenmiş, sonlu bir büyüklüktür — z_max = 12 ile z_max = 8 arasında sürekli bir geçiş yoktur, yalnız bir örgü tipinden diğerine kesikli sıçrama vardır.

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

6.3 Ek geçiş katı(1) → katı(2)

Bağlantılılığın sonluluğundan/kesikliliğinden, 2B modelde bulunmayan bir geçiş çıkar: katı(1) → katı(2) (polimorfizm/allotropi). Farklı z_max'lı iki katı faz (örneğin YMK ↔ HMK) farklı serbest entalpilere μ(τ, p) sahiptir; μ_{katı(1)} = μ_{katı(2)} olan sıcaklık/basınç birleşiminde madde erimeden örgüsünü değiştirir. Model resminde bu sırada koordinasyon geometrisi sıçramalı olarak değişir (erimedeki gibi yalnız doluluk z̄ değil) — normalde sabit tutulan referans geometrinin kendisinin sıçradığı tek durum. Gözlenebilir olan küçük bir c_v basamağı (z_max/4 değişimi) artı gizli dönüşüm ısısıdır; erime tepesinden belirgin biçimde zayıftır.

Birinci mertebe ölçütü ln G > 2 + 4ε/K (Böl. 5.5) şimdi 3B yorumun niteliksel olarak neden daha zengin olduğuna somut bir gerekçe verir: yalnız orada (G = z_max = 12, ln 12 = 2,485 > 2) erime birinci mertebeden olabilir — ve yalnız orada, bağlantılılık nicemlenmiş/sonlu olduğu için, kesikli örgü değişimi katı(1) → katı(2) (YMK ↔ HMK) vardır. Her iki olgu aynı nedenin iki kılığıdır: 3B paketlemenin daha yüksek, kesikli koordinasyonu hem keskin erime sıçramasını hem polimorfizmi üretir. Tek 2B örgüde (G = 6) bunların hiçbiri yoktur — Böl. 5.5 ile tutarlı.

Niceleme (bondkinetics_3d_15.py). Geçiş şimdi yalın biçimde hesaplanır. Her iki örgü klasik harmonik resimde birer serbest entalpi alır (ortak p·v ve öteleme payları karşılaştırmada düşer):

μ(τ) = h₀ − c_v·τ·ln(τ/θ), h₀ = −(ε/2)·z_max, θ = ħω/k_B (ε birimlerinde fonon sertliği),

Böl. 6.1/6.2'deki parametresiz karakteristik değerlerle:

Buradan c_v basamağı parametresiz olarak: |Δc_v| = z_max/4 değişimi = (12−8)/4 = 1,0 k_B (katı, YMK→HMK geçişinde ısı sığasında 1 k_B "hafifler") — erime tepesinden (birkaç k_B, Böl. 5.3/7.11) çok daha küçük.

Dönüşüm sıcaklığı τ_ss (μ_YMK = μ_HMK) ve gizli dönüşüm ısısı L = τ_ss·(S_HMK − S_YMK):

Sonuç (dürüst bulgu). Sıkı ideal resimde (bütün bağlar için aynı fonon sertliği, θ_HMK = θ_YMK) katı(1)→katı(2) geçişi yoktur: YMK hem entalpiyi (daha derin, −6'ya karşı −4 ε) hem de — daha yüksek bağ sayısı nedeniyle (c_v = 4,5 > 3,5) — titreşim entropisini kazanır; μ_YMK her sıcaklıkta μ_HMK'nin altında kalır. Ancak örgüye bağlı bir fonon sertliği θ_HMK < θ_YMK (daha açık HMK daha yumuşak titreşir, yani daha çok entropiye sahiptir) yüksek sıcaklıkta μ_HMK'yi μ_YMK'nin altına indirir — o zaman dönüşüm açılır ve erime tepesi τ_m'nin altında yatar, böylece katı(1)→katı(2)→sıvı→gaz sırası doğar (τ_ss < τ_m < τ_c ≤ τ_b(3B), Böl. 7.11). Yumuşak HMK sertliği sıkı ideal resmin dışında ılımlı bir genişletmedir (gerçekte demirin γ/δ polimorfizmi olarak bilinir: yüksek sıcaklıkta HMK, arada YMK); buna karşılık 1,0 k_B'lik c_v basamağı ve küçük pozitif dönüşüm ısısı model içi ve parametresizdir.

θ_HMK/θ_YMK gerçekte ne kadar küçük? (theta_springs_16.py). Dönüşümü serbest θ parametresinden kurtarmak için fonon sertliği artık mikroskobik olarak kestirilir (Böl. 6.6): merkezî çift yaylı bir Einstein yay modeli θ_HMK/θ_YMK ≈ 0,82 (yalnız en yakın komşular) ile ≈ 0,94 (HMK ikinci kabuğuyla) arasını verir. Buna karşılık μ karşılaştırmasından gereken açılma eşiği θ_HMK/θ_YMK < 0,61'dedir (bunun altında τ_ss vardır ve τ_m'nin altında yatar). Demek ki geometrik sertlik oranı çok büyüktür: sıkı ideal resimde YMK her sıcaklıkta kararlı kalır, salt geometrik olarak sürülen bir katı(1)→katı(2) geçişi yoktur. Gerçek polimorfizm (demir γ↔δ) ideal harmonik resmin dışında ek (elektronik/manyetik/anharmonik) katkılar gerektirir — tam da yukarıda varsayılan "ılımlı genişletme". Bu, dürüst bulguyu doğrular ve θ oranını serbest bir varsayımdan kestirilmiş, adlandırılmış bir model sınırına dönüştürür.

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

6.4 Düşünceler: sıvı suyun ses anomalisi

Normal sıvılarda v(T) tekdüze azalır (Böl. 3.3). Suyun ise ≈ 74 °C'de bir ses en büyüğü vardır: v önce T ile artar, sonra düşer. Modelde bu, iki yerel örgü tipinin yarışı olarak yorumlanabilir (Böl. 6.2): açık, düşük koordinasyonlu buz benzeri bir yakın düzen (küçük z_max, çok "delik") ve sıkı, daha yüksek koordinasyonlu bir yakın düzen. Düşük T'de açık yapı egemendir; ısıtmak önce açık yakın düzeni daha sıkı olan lehine kırar — etkin bağlantılılık ve onunla K_eff artar (v artar). Ancak en büyüğün üstünde olağan ÇATLAK bağ kırılması ağır basar, K_eff düşer, v azalır. Anomali böylece sıvı faz içinde bulanıklaşmış, eksik bir katı(1)↔katı(2) benzeri yakın düzen değişiminin ifadesidir — yalnız iki örgü tipli 3B modelde düşünülebilir, tek 2B örgüde değil.

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

6.5 Seçenek: misel dönmesi yerine fraktal girdap alanları

Yerel Brown hareketi (Böl. 3.4) misel dönmesi ve değiş tokuşu üzerinden açıklanmıştı. Seçenek olarak onu fraktal girdap alanlarının etkisi olarak da düşünmek olanaklıdır: ÇATLAK tarafından sürekli yeniden yapılan blok sınırları, yerel momentum akımlarından bir alan üretir; bu alan — türbülanslı akışlara benzer biçimde — bütün uzunluk ölçeklerinde özbenzer girdaplar taşır. O zaman bir sınama parçacığını kesikli misel çarpışmaları değil, ölçek değişmez bir girdap sürekliliği sürer. İki resim de (kesikli miseller ↔ fraktal girdaplar) aynı ÇATLAK kinetiğinin tümleyici sınır durumlarıdır; hangisinin geçerli olduğu sınama parçacığının boyu ile ortalama misel boyunun oranına bağlıdır. Bu ön tartışma salt nitelikseldir ve burada matematiksel olarak işlenmemiştir.

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

6.6 Mikroskobik fonon sertliği θ(örgü): bir Einstein yay modeli

Böl. 6.3'teki katı(1)→katı(2) karşılaştırması göreli fonon sertliği θ_HMK/θ_YMK'ye bağlıdır. Onu serbestçe seçmek yerine burada koordinasyon ve geometriden kestirilir (theta_springs_16.py).

Model. Bir merkez atom z komşusunun harmonik çanağında oturur; her bağ radyal bir çift yay κ'dır. Bir u yer değiştirmesi bağ doğrultularına ⟨cos²⟩ = 1/d = 1/3 (d = 3) ile izotropik olarak izdüşer, öyle ki etkin yay sabiti k_eff = (z/3)·κ olur. Einstein frekansı ω = √(k_eff/m), yani θ = ħω/k_B ∝ √k_eff. Aynı kütle ve aynı çift yay için şu çıkar:

θ_HMK/θ_YMK = √(z_HMK/z_YMK) = √(8/12) = 0,816.

İncelikler. (B) HMK'nin ikinci kabuğu (r₂/r₁ = 2/√3 = 1,155'te 6 atom) uzaklıkla ölçeklenen κ(r) ∝ (r₁/r₂)ⁿ yayıyla (n = 6…8, Lennard-Jones benzeri) oranı 0,91…0,94'e yükseltir (yakın ikinci kabuk HMK'yi etkin biçimde sertleştirir). (C) HMK'nin aynı yoğunlukta daha kısa bağı (r_nn(HMK)/r_nn(YMK) = 0,972, −2,8 %) yayı hafifçe sertleştirir ve 0,89…0,91 verir. Toplamda geometrik bir pencere θ_HMK/θ_YMK ≈ 0,82…0,94.

Açılma eşiği. μ karşılaştırması (Böl. 6.3) dönüşümü (τ_ss vardır ve τ_m = 1,70'in altındadır) ancak θ_HMK/θ_YMK < 0,61 için açar. Geometrik pencere bunun belirgin biçimde üstündedir.

Bulgu. Salt yay modeli HMK'yi yeterince yumuşak yapmaz: sıkı ideal resimde YMK kararlı kalır, geometrik olarak sürülen bir katı(1)→katı(2) geçişi doğmaz. Böylece θ_HMK/θ_YMK serbest bir varsayımdan, doğal değeri (≈ 0,8–0,9) dönüşümü kapalı tutan kestirilmiş bir büyüklüğe dönüşmüştür — Böl. 6.3'teki dürüst bulguyla tutarlı olarak. Gerçek polimorfizm böylece kanıtlanabilir biçimde ideal harmonik resmin dışındadır (elektronik/manyetik paylar, anharmoniklik), salt örgü geometrisinde değil.

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

6.7 Yarı kararlılığın matematiksel kestirimi (2B üçgen ↔ kare)

Böl. 3.2'deki niteliksel tartışma için burada 2B resimdeki ilgili büyüklük mertebesi hesabı verilir (sert diskler, çap D, temas uzaklığı a = D).

(1) Enerji. Parçacık başına U/N = −(ε/2)·z (her bağ yarı): üçgen z = 6 → −3 ε, kare z = 4 → −2 ε. Üçgen → kare yeniden yapılanması demek ki parçacık başına ΔU = +1 ε tutar (endotermik; gevşek biçim enerjice daha zengindir).

(2) Hacim (alan) ve iş. Parçacık başına alan: üçgen A_tri = (√3/2)·D² ≈ 0,866 D², kare A_quad = D². Genişleme ΔA = A_quad − A_tri ≈ 0,134 D² (+15,5 %)'dir — çizgi romanda gösterilen hacim büyümesi. (İndirgenmiş) bir p basıncına karşı bunun için W = p·ΔA işi harcanmalıdır; serbest entalpi Δμ = ΔU + p·ΔA − τ·ΔS kadar artar. Entropi terimi ΔS > 0 (gevşek biçimin daha yumuşak kipleri, daha büyük c_v basamağı Δc_v = (6−4)/4 = 0,5 k_B vardır) μ_quad'ı ancak yüksek τ'da μ_tri'nin altına indirir — bu yüzden üçgen düşük sıcaklık taban durumu, kare ise yüksek sıcaklık/yarı kararlı biçimdir.

(3) Çekirdeklenme engeli (2B). Yarı kararlılık bir engel gerektirir. 2B'de yeni fazın bir çekirdeği (yarıçap r) alan kazancı −Δμ_A·π r² ile (Δμ_A = alan başına itici güç) kenar katkısına +λ·2π r karşı büyür (λ = uzunluk başına kenar enerjisi): ΔG(r) = −Δμ_A·π r² + 2π λ r. En büyük değer r\* = λ / Δμ_A'da, engel ΔG\* = π λ² / Δμ_A. Kritik çekirdek r\, itici güç Δμ_A büyüdükçe küçülür; ısıl ölçek τ = k_BT engele ulaşmadığı sürece (**ΔG\ ≳ birkaç k_BT), gevşek (ya da kaynamada aşırı doymuş) faz kinetik olarak donmuş kalır. Küçük öbekler için (2–3 parçacık, aşırı doyma şeridindeki gibi) kenar/alan oranı en elverişsizdir — bunlar r\*'ın solunda yatar ve ortalamada yeniden çözülür. Böylece çizgi romanda gösterilen etkinleşme enerjisi, hem yarı kararlı kristal biçimini hem de aşırı doymuş gaz fazını taşıyan aynı ΔG\* büyüklüğüdür; o, ideal (engelsiz) denge modeline kinetik** bir ektir.

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

6.8 θ(örgü)'nün Debye benzeri inceltilmesi

Einstein yay modeli (Böl. 6.6) bütün fonon dağılımını tek bir frekansla değiştirir. Buradaki yaklaşım θ(örgü)'yü her örgü için yalın bir fonon durum yoğunluğuyla inceltir (theta_debye_17.py): aynı merkezî çift yaylar κ(r) = r^(−n), en yakın ve ikinci komşu kabukta dinamik matrisi D(k) = Σ_j κ_j·(R̂_j⊗R̂_j)·(1 − cos(k·R_j)) oluşturur; bu matris ilkel Brillouin bölgesinin bir k örgüsünde (Monkhorst-Pack, 12³ nokta) köşegenleştirilir (analitik 3×3 özdeğer çözücü, yalnız standart kütüphane). Üç fonon dalından iki Debye benzeri sertlik ölçüsü frekans momentleri olarak çıkar:

• 2. moment θ(2) ∝ √⟨ω²⟩ — Einstein sınır durumu (1 − cos'un BZ ortalaması 1 verir, yani ⟨ω²⟩ = z·κ/3); yalnız en yakın komşu denetimi, salt geometri: √⟨ω²⟩_HMK / √⟨ω²⟩_YMK = 0,816 = √(8/12) — tam da Böl. 6.6'daki Einstein değeri.
• log momenti θ(0) ∝ exp⟨ln ω⟩ — asıl Debye benzeri büyüklük; Einstein değerinin yok saydığı yumuşak akustik dalları ağırlıklandırır.

Ölçek taşıyan bütün sabitler (ħ/k_B, κ ölçeği, Debye önçarpanı) θ_HMK/θ_YMK oranında düşer; yalnız DOS biçimi kalır. HMK'yi mekanik olarak ancak stabilize eden ikinci komşu kabukla ve n = 0…8 yay yasasıyla şu çıkar:

θ_HMK/θ_YMK ≈ 0,87…1,00 (moment tanımı ve yay yasası üzerinden hata bandı).

Bu, salt Einstein değerinden (0,82…0,94) daha yüksektir, çünkü aynı yoğunlukta daha kısa olan HMK bağı yayları sertleştirir. μ karşılaştırmasından (Böl. 6.3) gereken 0,61 açılma eşiği hâlâ belirgin biçimde altındadır: fonon DOS'u da HMK'yi yeterince yumuşak yapmaz. Böl. 6.6'daki bulgu böylece daha ince yaklaşımla doğrulanır ve kesinleşir — sıkı ideal resimde YMK her sıcaklıkta kararlı kalır; salt geometrik olarak sürülen bir katı(1)→katı(2) geçişi doğmaz, gerçek polimorfizm (demir, Böl. 7.13) ideal harmonik resmin dışında katkılar gerektirir.

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

6.9 Lennard-Jones tam örgüsünden θ(örgü)

Böl. 6.8 nn + 2nn üzerinde kesilmiş bir merkezî yay κ(r) = r^(−n) kullanır. Sürüm 17 aynı soruyu tam örgü üzerinde bir çift potansiyelden gelen gerçek kuvvet sabitiyle sınar (theta_ljlattice_18.py): her bağ, Lennard-Jones potansiyeli φ(r) = 4ε[(σ/r)¹² − (σ/r)⁶] içinden boyuna pay κ_L = φ″(r) ve enine pay κ_T = φ′(r)/r ile bir yaya dönüşür. Dinamik matris

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

bir kesime kadar (2,0…3,0 nn) bütün komşu kabuklar üzerinden toplanır; her örgü önceden kendi LJ dengesine gevşetilir (serbest örgü sabiti, aynı yoğunluk dayatılmaz). Bu, yalın nn+2nn yaklaşımının (Böl. 6.8) sertlik sıralamasını doğru verip vermediğinin dürüst sınamasıdır.

Sonuç (belirlenimci, yalnız standart kütüphane):

İki bulgu, ikisi de Böl. 6.6/6.8'i keskinleştirir:

1. YMK, kararlı LJ taban durumudur (3,0 nn kesiminde atom başına E ≈ −8,4 ε, HMK için −8,0 ε'ya karşı) — çift potansiyelde sıkı paketlemenin bilinen önceliği.
2. İdeal LJ HMK'si kaymaya kararsızdır: fonon kiplerinin yaklaşık %3…3,5'i sanaldır (ω² < 0, Born kararsızlığı). Salt çift potansiyelli ideal resim demek ki mekanik olarak kararlı bir HMK'yi hiç taşımaz; geriye kalan gerçel kipler üzerinden θ_HMK/θ_YMK ≈ 1,0 kalır — μ karşılaştırmasından gereken 0,61 açılma eşiğinin altına hiç inilmez.

Böylece Böl. 6.6/6.8'in savı, kesilmiş yaydan tam örgüye ve salt boyuna modelden tam dinamik matrise geçen inceltmeye karşı gürbüzdür: klasik ideal resimde salt geometrik bir katı(1)→katı(2) geçişi doğmaz — gerçek polimorfizm (demir Böl. 7.13, titanyum Böl. 7.15) çift potansiyelin dışında katkılar gerektirir (elektronik/manyetik, anharmonik). Gerçek bir ab initio fonon hesabı bilinçli olarak dışarıda bırakılmıştır (yalın ideal resim niteliğinden çıkardı).

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

7 — Kalibrasyon ve hâl denklemi

Özet. Bölüm 7, modeli tek enerji parametresi ε üzerinden gerçek maddelere kalibre eder ve sıvı/gaz dalını hâl denklemine dek genişletir:

• 7.1 / 7.2 Kalibrasyon yöntemi ve tablosu: erime noktasında tek bir çapa noktasıyla doğrusal sıcaklık bağlantısı T = (ε/k_B)·τ; buradan inceltilmiş evrensel oran R = τ_b/τ_m ≈ 1,54.
• 7.3 / 7.4 Sert çekirdekli van der Waals hâl denklemi: indirgenmiş denklem p_r = 8τ_r/(3v_r − 1) − 3/v_r², kritik nokta (1, 1, 1) ve Maxwell birlikte varoluşu; a ∝ ε, b ∝ D². Kaynama noktaları dizgesel olarak T_b/T_c ≈ 0,48'de yatar (Guldberg açıklığı); bir Redlich–Kwong karşılaştırma hesabı bunları ≈ 0,60'a yükseltir — adlandırılmış bir hata bütçesi olarak, hâl denklemi değiştirilmeden.
• 7.5 Ortak (τ, v) yüzeyi: birlikte varoluş kubbesi (kaynama), erime çizgisi ve sert duvar bir üçlü noktayla birleşir (v_solid = 0,40 için τ_r,triple ≈ 0,47).
• 7.6 Gizli ısı: L ≈ 1,41 ε, ΔS_fus ≈ 1,175 R verir (Richards kuralına yakın); metaller mükemmel tutturulur (bakır ×1,00), iyonik/H-köprülü/moleküler kristaller ideal resmin dışındadır.
• 7.7 ε ↔ ε_lg birleştirmesi: k_B T_c = z_max·ε/4 iki ölçeği eşler (ε_lg = (27/32)·z_max·ε); 3B öngörüsü τ_c/τ_m = 2,50 basit moleküler maddeleri tutturur, 2B özel durum olarak kalır.
• 7.8 / 7.9 Parametresiz üçlü nokta: v_solid bağlantılılıktan değil rijitlik eşiğinden z̄_c = 2d çıkar; 3B'de temiz bir sıvı penceresi açılır (τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00), 2B'de kapalı kalır (üç kez doğrulanmış 2B özel rolü).
• 7.10 / 7.11 Gerçek 3B taramalar: YMK/HMK üzerinde bağ perkolasyonu z̄_c(3B) ≈ 1,5–1,6 verir; bağ kinetiğinin τ taraması τ_m ≈ τ_rig (erime = rijitlik, düşük) ve τ_b(3B) ≈ τ_c (kaynama = kritik noktada bağlantılılık çözülmesi) gösterir.
• 7.12 / 7.13 / 7.15 katı(1)→katı(2) çizgisi: τ_ss noktasından basınç üzerinden bir çizgi olur (dp/dτ = ΔS/Δv > 0, Δv = +8,9 %); (τ, v) yüzeyine dördüncü sınır çizgisi olarak ve demire (γ→δ) kalibre edilmiş bir (τ, p) çizgisi olarak, titanyumda (α→β) doğrulanmış ve argonla negatif kontrollü.
• 7.14 Redlich–Kwong dalı: van der Waals varsayılanının yanında kendi kalibrasyonlu bağımsız bir ideal resim türü olarak (T_b/T_c ≈ 0,60, hata −1 %).

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

7.1 Kalibrasyon yöntemi (tek enerji parametresi ε) — uyarlandı

Doğrusal: T = (ε/k_B)·τ, madde başına bir çapa noktası. Erime noktasına çapalama: ε/k_B = T_m,deney/τ_m, öngörü T_b,pred = T_m,deney·(τ_b/τ_m). İnceltilmiş erime tepesiyle (işbirlikçi terim) c_v en büyüğü artık τ_m = 1,20'de (1,15 yerine), bağlantılılık perkolasyonu (kaynamaya yakın) değişmeden τ_b = 1,85'te. Buradan inceltilmiş evrensel oran:

R = τ_b / τ_m = 1,85 / 1,20 ≈ 1,54 (sürüm 08: 1,61).

Sürüm 08'in kalibrasyon tablosu (Böl. 7.2) R ile doğrusal ölçeklenir: bütün T_b,pred değerleri 1,54/1,61 ≈ 0,957 çarpanıyla düşer (≈ −4,3 %). Bu, soy gazların fazla tahminini hafifçe iyileştirir (argon +54 % → +48 %) ve NaCl/demir isabetlerini az miktarda kötüleştirir (NaCl −0,6 % → −4,8 %). Temel sav aynı kalır: R, eşyönlü sıkı paketlenmiş, yönsüz bağlanan maddeler için uyar; soy gazlar/zincirler ideal resmin dışındadır.

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

7.2 Kalibrasyon tablosu

modell_08.md Böl. 7.2'den değiştirilmeden alındı (7.1'deki doğrusal R ölçeklemesi uyarısıyla). cv_scan_08.py (tablo) ya da cv_scan_09.py (τ_m, τ_b, R) ile yeniden üretilebilir.

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

7.3 Yoğunluk/hacim ekseni: sert çekirdekli van der Waals hâl denklemi

Basınç "parametre" olmaktan çıkıp gerçek değişkene yükseltilir. SERT ÇEKİRDEK öz hacim b'yi, bağ çekimi (ε·z̄) kohezyon basıncı a'yı verir. İndirgenmiş birimlerde (kritik nokta (v_r, τ_r, p_r) = (1, 1, 1)'de):

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

Maxwell kurgusuyla birlikte varoluş (cv_scan_09.py, eşit alan kesiti):

Sonuç. τ_c'nin altında s/g yoğunluk sıçraması açılır (v_gas/v_liq 2,5'ten 11,5'e büyür; ne kadar soğuksa o kadar büyük). Bu, kaynama çizgisidir — artık sabit yoğunluklu örgü üzerinde yaşayan erime c_v tepesinden (Böl. 5.3) temiz biçimde ayrılmıştır. Böylece sürüm 08'de açık kalan 3. nokta kapanır: erime (sabit yoğunlukta rijitlik/bağ düzeni) ile kaynama (yoğunluk/hacim sıçraması) aynı modelin farklı eksenleridir.

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

7.4 Niceliksel s/g kalibrasyonu: a → ε, b → D²

İndirgenmiş büyüklüklerin gerçek maddelere eşlenmesi. van der Waals sabitleri kritik noktaya bağlıdır:

b = k_B·T_c / (8 p_c) (öz hacim, SERT ÇEKİRDEK), a = 27 k_B²·T_c² / (64 p_c) (kohezyon, BAĞ).

Böylece a ∝ ε (kohezyon bağ enerjisinden gelir) ve b ∝ D² (2B öz hacim sert çekirdek çapından). s/g dalının doğal enerji ölçeği

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

Çapa maddesi argon (T_c = 150,87 K, p_c = 48,63 bar): b = 5,35·10⁻²⁹ m³/molekül, a = 3,76·10⁻⁴⁹ J·m³/molekül², ε_lg/k_B = 509 K ve bir molekül çapı D ≈ b^{1/3} = 3,8·10⁻¹⁰ m — argonun van der Waals çapına (~3,4·10⁻¹⁰ m) göre %11 duyarlıkta.

Kaynama sıcaklığının öngörüsü. Sabit dış basınçta kaynama noktası, indirgenmiş doyma basıncının p_sat(τ_r)·p_c = p_atm koşulunu sağladığı yerdedir; buradan T_b = T_c·τ_r,b (lg_calib_10.py'den Maxwell doyma eğrisi):

*CO₂ 1 atm'de süblimleşir (T = 194,7 K = süblimleşme noktası).

Sonuç. Moleküler maddeler için (p_c ≈ 30–70 bar) model baştan sona T_b/T_c ≈ 0,48–0,51 verir; bu, ölçüm değerinin dizgesel olarak ~%17–24 altındadır. Bu dar anlamda bir model hatası değil, van der Waals denkleminin görgül Guldberg kuralına (T_b/T_c ≈ 0,6) bilinen uzaklığıdır: van der Waals izotermi düşük sıcaklıklarda doyma basıncını olduğundan az verir. Su (ve T_c/p_c değerleri zaten yalnız kabaca kestirilmiş metaller) p_r,atm < 0,005'te anlamlı van der Waals alanının dışında kalır (τ_r,b < 0,40) — orada kohezyon artık van der Waals türü değildir (hidrojen köprüleri, metalik bağ). Kaynama noktalarının büyüklük mertebesi ve sıralaması doğru tutturulur.

Hata tablosu olarak Guldberg yerleştirmesi. Uzaklık burada yalnız adlandırılır, bir hâl denklemi değişimiyle düzeltilmez (ideal resim van der Waals türü kalır). Görgül Guldberg kuralı (Guldberg 1890), kaynama noktasının yaklaşık T_b/T_c ≈ 0,6'da yattığını söyler (hafif madde bağımlılığı ~0,58–0,64). İndirgenmiş van der Waals modeli bunun yerine sabit bir değer T_b/T_c = τ_r,b ≈ 0,48–0,51 verir:

Yorum. Model ortalaması 0,49, Guldberg sabiti 0,60'ı yaklaşık −%18 ile tutturur; sapma dizgesel ve tek yönlüdür (hep çok düşük), yani moleküler resmin değil van der Waals izoterminin bir özelliğidir (düşük τ'da çok küçük doyma basıncı). Daha iyi bir kohezyon hâl denklemi (örn. Redlich–Kwong) τ_r,b'yi yükseltirdi — ama burada bilinçli olarak ideal resmin dışında kalır.

Bir Redlich–Kwong karşılaştırma hesabıyla hata bütçesi (guldberg_rk_15.py). van der Waals'ın Guldberg kuralına uzaklığını sayısallaştırmak için (iyileştirmek için değil), aynı kaynama noktası sınaması Redlich–Kwong türü, sıcaklığa bağlı bir a/√T kohezyonuyla yinelenir. İndirgenmiş birimlerde (kritik nokta (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,

Maxwell eşit alan kurgusu ve yukarıdakiyle aynı indirgenmiş dış basınç p_r,atm ile. Sonuç (çapa madde başına bir sayı, T_b/T_c = τ_r,b):

Sonuç. Sıcaklığa bağlı RK kohezyonu T_b/T_c'yi vdW ≈ 0,48'den ≈ 0,60'a yükseltir — neredeyse tam Guldberg sabiti ve gerçek değerler (≈ 0,58–0,63); ölçüm değerlerine karşı ortalama hata −%21'den (vdW) −%1'e (RK) düşer. Böylece Guldberg açıklığının nedeni açıkça sayısallaşmıştır: van der Waals kohezyonunun eksik sıcaklık bağımlılığıdır (düşük τ'da çok küçük doyma basıncı), moleküler resim değil. Bu karşılaştırma hesabı bilinçli olarak ideal resmin dışındadır; model van der Waals türü kalır (hâl denklemi değişikliği yok) — hata bütçesi yalnız adlandırılmış, içeri alınmamıştır.

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

7.5 Üçlü noktalı ortak (τ, v) yüzeyi

Erime ve kaynama ekseni şimdiye dek ayrı yaşadı: erime (rijitlik/bağ düzeni) sabit yoğunlukta, kaynama (van der Waals) yoğunluk sıçraması üzerinden. Ortak bir (τ_r, v_r) düzleminde birleşirler:

• Birlikte varoluş kubbesi (kaynama): Maxwell kurgusundan v_liq(τ_r) (sol dal) ve v_gas(τ_r) (sağ dal); bunlar kritik noktada (τ_r, v_r) = (1, 1) buluşur.
• Katı (erime): sert duvarın hemen üstünde v_r > 1/3'te, yüksek ve neredeyse sabit yoğunlukta dar bir şerit (indirgenmiş hacim v_solid ≈ 0,40).
• Erime çizgisi: v_solid'de neredeyse dikey (sabit yoğunluk), Maxwell rijitliğinin sıcaklığı (z̄_c = 2d; Böl. 7.8 v_solid'i bundan parametresiz türetir — eski sürüm-11 eşlemesi "z̄_c ≈ 2,08" böylece düzeltilmiştir).

Üçlü nokta, sıvı birlikte varoluş sınırı v_liq(τ_r)'nin katı paketleme yoğunluğu v_solid'e ulaştığı noktadır — bunun altında artık sıvı dal yoktur (doğrudan süblimleşme):

Sonuç. v_solid = 0,40 için üçlü nokta τ_r,triple ≈ 0,47'de, p_sat ≈ 0,018 ile ve buhar hacmi v_gas ≈ 68 ile yatar (çok seyreltik süblimleşme buharı). τ_triple/τ_c ≈ 0,47 oranı gerçek pencerede yatar (su 0,42; argon 0,56; CO₂ 0,71). Katı ne kadar yoğunsa (v_solid ne kadar küçükse), indirgenmiş üçlü nokta o kadar düşük ve birlikte var olan buhar o kadar seyreltiktir — niteliksel olarak doğru. Böylece her iki geçiş ve üçlü nokta tek bir (τ, v) yüzeyinde çapalanmıştır; kritik nokta (1, 1, 1)'de oturur. Sürüm 15, Böl. 7.12'de bir dördüncü sınır çizgisi ekler — katı(1)/katı(2) dönüşümü (YMK↔HMK), τ_ss'de birbirine dönüşen iki neredeyse dikey katı yoğunluk şeridi olarak (v ≈ 0,667 ve 0,726).

Daha önceki geliştirme aşaması: v_solid = 0,40 burada henüz konulmuştu. Böl. 7.8 v_solid'i perkolasyon eşiği z̄_c'den parametresiz türetir ve üçlü noktayı gerçek bir öngörüye dönüştürür.

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

7.6 Gizli ısıya karşı erime entalpileri

Model gizli ısısı L ≈ 1,41 ε (Böl. 5.3), erime noktası kalibrasyonu ε/k_B = T_m/τ_m (τ_m = 1,20) üzerinden gerçek birimlere çevrilir:

ΔH_m,pred = L·R·T_m/τ_m = 1,175·R·T_m, buna karşılık ΔS_fus,pred = L·R/τ_m = 1,175·R.

Bu — kayda değer biçimde — neredeyse tam olarak Richards kuralıdır (basit yapılı katılar için ΔS_fus ≈ R). Ölçülen erime entalpileriyle karşılaştırma (lg_calib_10.py):

Sonuç. Model büyüklük mertebesini (kJ/mol) baştan sona ve basit metaller için Richards kuralını mükemmel tutturur (bakır ×1,00, demir ×1,28); sıvılaştırılmış gazlar hemen altında kalır (metan ×0,94, azot ×0,87, argon ×0,69). İyonik, H-köprülü ve moleküler kristallerin (NaCl ×0,37, su ×0,44, benzen ×0,28, CO₂ ×0,23) gerçek erime entropileri belirgin biçimde daha yüksektir (ΔS_exp/R = 2,6–5,0), çünkü erimede tek düzen parametreli modelin kavramadığı ek düzen dereceleri serbest kalır (iyon çiftlenimi, çok atomlu moleküllerin yönelimi). Sınır böylece açıkça adlandırılmıştır: ideal resim eşyönlü, yönsüz bağlanan maddeler için geçerlidir (uydurma yok).

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

7.7 Enerji ölçeklerinin birleştirilmesi: ε ↔ ε_lg

Şimdiye dek erime dalı bağ enerjisi ε'yi (τ_m = 1,20), s/g dalı ayrı bir kohezyon ölçeğini ε_lg/k_B = 3,375·T_c taşıdı (Böl. 7.4). Her ikisi tek bir model parametresinden türetilebilir. van der Waals kohezyonu a, modelde molekül başına bağ enerjilerinin toplamından doğar (BAĞ): her molekül ε enerjili z_max/2 bağ taşır. Komşu çekimi ε olan bir ortalama alan Bragg–Williams örgü gazı kritik sıcaklık için klasik bağıntıyı verir:

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

Bu, ε_lg = a/b = (27/8)·k_B·T_c'de (Böl. 7.4) yerine konursa iki ölçeğin eşleşmesi çıkar:

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

Böylece erime tepesi, perkolasyon ve kritik nokta tek bir ε ölçeğinde yatar (unify_scale_12.py):

Öngörü τ_c/τ_m = z_max/(4·τ_m): 2B → 1,25, 3B → 2,50. Gerçek maddelere karşı:

Sonuç. 3B öngörüsü τ_c/τ_m = 2,50, basit moleküler maddelerin gerçek göstergesini (T_c/T_m ≈ 2,0–2,1) iyi tutturur; metaller daha yüksektedir (T_c kabaca kestirilmiş, farklı bağ türü). 2B öngörüsü (1,25) çok düşüktür — ve τ_c(2B) = 1,50 hatta τ_b = 1,85'in altındadır (kritik nokta perkolasyon eşiğinin altında kalırdı). Bu bir hesap hatası değil, bilinen 2B özel rolüdür: kohezyon ve bağ ölçekleri orada temiz sıralanmaz, gürbüz bir sıvı penceresi yoktur (bkz. sürekli 2B erime, Böl. 5.5). Ölçekler ancak 3B'de (τ_m = 1,20 < τ_c = 3,00) fiziksel olarak anlamlı biçimde ayrılır.

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

7.8 z̄_c'den parametresiz üçlü nokta

Böl. 7.5 v_solid = 0,40'ı henüz koymuştu. Daha önceki geliştirme aşaması katı paketleme yoğunluğunu perkolasyon eşiğinden türetir. Bir örgünün indirgenmiş hacim v_r'deki koordinasyon sayısı yoğunluk–koordinasyon bağıntısını izler (sert duvar v_r = 1/3 ↔ en sıkı paketleme, z̄ = z_max):

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

Katı, z̄ eşik z̄_c'nin altına iner inmez sertliğini yitirir; buna karşılık gelen yoğunluk v_solid'dir. Üçlü nokta o zaman sıvı birlikte varoluş sınırı v_liq(τ_r)'nin bu yoğunluğa ulaştığı noktadır (unify_scale_12.py, 2B, z_max = 6):

İyileştirilmiş sürüm. Erime yoğunluğu v_solid'i genel Maxwell rijitliği z̄_c = 2d belirler (2B → 4, 3B → 6), bağ bağlantılılığı değil (gerekçe Böl. 5.4a: erime sertliğin yitimidir, bağlılığın değil). Böylece güvenilir üçlü nokta şudur (rigidity_3d_13.py, 2B):

Sonuç 2B. (Fiziksel olarak doğru) rijitlik eşiğiyle bile tek 2B örgü temiz bir sıvı penceresi vermez: v_solid = 0,50 gerçi gerçek gösterge τ_triple/τ_c ≈ 0,77'yi verir, ama mutlak ε birimlerinde τ_triple = 0,77·1,50 = 1,16'dır ve böylece τ_m = 1,20'nin altında kalır. Bağlantılılık eşiği (z̄_c = 2,08) v_solid ≈ 0,96 ≈ v_c'yi birlikte varoluş kubbesinin dışına bile iter — üçlü ve kritik nokta çakışır. İki hesap da 2B özel rolünü doğrular (Böl. 5.5, 7.7): tek 2B örgü bağımsız bir sıvı fazı pek taşımaz. Bu dürüst bir bulgudur, dar anlamda bir hesap ya da model zayıflığı değil.

Belirleyici sınama bu yüzden 3B kontrol hesabıdır (Böl. 7.9): orada aynı rijitlik kuralı (z̄_c = 6 → v_solid = 0,667) temiz bir pencere açar. Her durumda üçlü nokta artık parametresiz olarak z̄_c'ye bağlıdır — serbest parametre v_solid ortadan kalkmıştır.

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

7.9 3B kontrol hesabı: tek bir ε ölçeğinde temiz sıvı penceresi

Bütün işaretler tek bir ε ölçeğinde yattığında 2B çöküşünün 3B'de çözülüp çözülmediğinin sınanması (Böl. 7.7'den τ_c = z_max/4; genel rijitlikten v_solid, z̄_c = 2d; üçlü nokta, indirgenmiş van der Waals sıvı yoğunluğu v_liq(τ_r)'nin katı yoğunluğa ulaştığı yerde). Sonuç (rigidity_3d_13.py):

Sonuç. 3B'de temiz bir sıvı penceresi açılır: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00 (pencere genişlikleri: τ_triple − τ_m = 1,63; τ_c − τ_triple = 0,18). Üst kenar dardır — üçlü nokta kritik noktanın hemen altındadır — ama sıra doğrudur ve sıvı faz bağımsız bir alan olarak vardır. 2B'de pencere kapalı kalır (τ_triple = 1,16 < τ_m). Böylece 2B özel rolü üçüncü kez bağımsız olarak belgelenmiş ve aynı zamanda bunun bir boyut etkisi olduğu gösterilmiştir: aynı rijitlik kuralı modeli 3B'de iyileştirir. Dar üst kenar Böl. 7.10'da gerçek bir 3B perkolasyon taramasıyla ve bir duyarlık çözümlemesiyle pekiştirilir.

Model sınırları kutusu — 2B özel rolü. Üç bağımsız hesap aynı bulguyu gösterir: tek 2B örgü bağımsız bir sıvı fazı pek taşımaz. 1. τ_c < τ_b (Böl. 7.7): kritik nokta (1,50) bağlantılılık perkolasyonunun (1,85) altında kalırdı — kohezyon ve bağ ölçekleri 2B'de temiz sıralanmaz. 2. v_solid ≈ v_c (Böl. 7.8): bağlantılılık erime yoğunluğu (0,96) neredeyse kritik yoğunlukla çakışır — üçlü ve kritik nokta çöker. 3. τ_triple < τ_m (Böl. 7.9): doğru Maxwell rijitliğiyle bile üçlü nokta (1,16) erime tepesinin (1,20) altındadır — pencere yok.

Bu dürüst bir bulgudur, hesap ya da model zayıflığı değil: bilinen boyut etkisidir (sürekli 2B erime, KTHNY benzeri; Böl. 5.5). 3B'de aynı rijitlik kuralı modeli iyileştirir (Böl. 7.9/7.10). İkinci bir yakın düzenin (Böl. 6.4) 2B resmi ayrıca iyileştirip iyileştiremeyeceği bir genişletme sorusudur, ideal resmin bir eksikliği değil. Model canlandırmaya hizmet edeceği için bu zayıflığı kabul ediyor ve modelin aksiyomatik bir genişletmesinden vazgeçiyorum.

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

7.10 Gerçek 3B bağ perkolasyon taraması (YMK/HMK) ve pencere kenarı duyarlılığı

Böl. 7.8/7.9, 2B bağlantılılık değeri z̄_c = 2,08'i düşüncede 3B'ye taşıdı (orada v_solid = 1,92, fiziksel değil) ve 3B savını genel rijitliğe z̄_c = 2d dayadı. Sürüm 13 bağlantılılığı 3B'de artık doğrudan sınar: YMK (z_max = 12) ve HMK (z_max = 8) üzerinde Union-Find ile Monte Carlo bağ perkolasyonu, periyodik sınırlar, belirlenimci LCG tohumu (rigidity_3d_14.py). Eşik p_c = P_∞(p) = 0,5 kesişimi (Böl. 5.6'daki 2B sonlu boyut incelemesiyle aynı yöntem):

MC değerleri ılımlı kutu kenarı nedeniyle (L = 12, saf Python) yazın eşiklerinin biraz üstündedir (sonlu boyut yanlılığı; 2B'de olduğu gibi kestiriciye göre aşağıdan ya da yukarıdan); sıralama ve büyüklük mertebesi tutar (YMK daha yoğun → HMK'den küçük p_c).

Yorum (bağımsız 3B doğrulama). 3B bağlantılılık eşikleri z̄_c(3B) ≈ 1,5–1,6'dadır — 2B değeri 2,08'den bile daha düşük ve Maxwell rijitliği 2d = 6'nın çok altında. Bundan çıkan "bağlantılılık erime yoğunluğu" v_solid = 1/(3 p_c) ≈ 1,7 (HMK) ya da 2,6 (YMK), birlikte varoluş kubbesinin (v_liq,max ≈ 0,94) çok dışındadır. Böylece — 2B değerini 3B'ye taşımadan — doğrudan gösterilmiştir: 3B'de bağlantılılık erime yoğunluğunu hele hiç belirleyemez; o, rijitlikten gelmelidir (z̄_c = 6 → v_solid = 0,667). İki eşik 3B'de 2B'dekinden daha da net ayrılmıştır. Gerçek τ_b(3B) (z̄(τ)'nun düşük bağlantılılık eşiği 1,5'e ulaştığı sıcaklık) böylece 2B değeri τ_b'nin üstünde, kritik noktaya çok yakın olurdu — "bağlantılılık çözülmesi = kaynama" resmiyle tutarlı.

Dar üst pencere kenarının duyarlılığı (τ_c − τ_triple ≈ 0,18). v_solid = 0,667 çevresinde sayısal olarak (±%2) dτ_triple/dv_solid ≈ 1,55 çıkar. Pencere ancak v_solid ≈ 0,11 artarsa (≈ 0,78'e) kapanırdı (τ_triple = τ_c = 3,00), yani etkin 3B rijitlik eşiği 2d = 6'nın belirgin biçimde üstünde olsaydı (z̄_c ≳ 7,7). Genel Maxwell sayımı için (z̄_c = 2d = 6) τ_m < τ_triple < τ_c sırası demek ki gürbüz kalır; üst kenar gerçekten dardır ama tehlikede değildir. Kesikli örgü değişimi katı(1)→katı(2) (YMK↔HMK) sürüm 14'te nicelenmiştir (Böl. 6.3); bağ kinetiğinin gerçek 3B τ taraması Böl. 7.11'de gelir.

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

7.11 Bağ kinetiğinin gerçek 3B τ taraması: τ_b(3B) ve tek ε ölçeğinde erime çizgisi

Böl. 7.8–7.10, 3B savını genel rijitliğe (z̄_c = 2d = 6) ve durağan bağlantılılık taramasına (z̄_c(3B) ≈ 1,55) dayadı. Açık kalan, bağ kinetiğinin asıl τ taramasıydı: kendiyle tutarlı z̄(τ) bu düşük bağlantılılık eşiğine sıcaklık ekseninin neresinde ulaşır (= gerçek τ_b(3B)) ve erime tepesi nerede yatar? Sürüm 14, kendiyle tutarlı z̄(τ)'yu (Böl. 4.2, x = 1/(1 + G·e^{−(ε+K·x)/τ}]) doğrudan YMK (z_max = 12) ve HMK (z_max = 8) üzerinde, G = z_max ve 2B referans durumdaki aynı işbirliklikle J = 0,55 hesaplar; n_blk(τ) bir YMK/HMK bağ Monte Carlo'sundan (Union-Find) gelir ve

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄ − (K/2)·x² (3B, işbirlikçi terimle),

c_v = dU/dτ (bondkinetics_3d_15.py, belirlenimci LCG tohumu; n_blk yumuşatılmıştır, çünkü ona yalnız zayıf blok dönme terimi bağlıdır). Üç işaret okunur: c_v tepesi τ_m (erime), rijitlik işareti τ_rig (z̄(τ) 2d = 6'yı keser) ve gerçek τ_b(3B) (z̄(τ) Böl. 7.10'daki bağlantılılık eşiği z̄_c(3B)'yi keser):

Sonuç. İki bulgu, ikisi de şimdiye dekki resimle tutarlı:

1. Erime aşağıda yatar ve iki kez belgelenmiştir. Termodinamik erime tepesi τ_m (YMK 1,70, HMK 1,40) pratikte rijitlik işareti τ_rig ile (z̄ = 6; YMK 1,73, HMK 1,35) çakışır — ağın sertliğini yitirdiği aynı sıcaklık. Bu, 3B erime çizgisini ε ölçeğine çapalar: T_m = (ε/k_B)·τ_m, örneğin argon için (T_m = 83,8 K) ε/k_B = 83,8/1,70 ≈ 49 K. K = J·z_max 3B'de 2B'dekinden büyük olduğundan (3,30 yerine 6,60), 3B erime tepesi 2B değeri 1,20'nin üstünde yatar — çelişki değil, daha yoğun paketlemenin daha güçlü işbirlikçiliğidir.
2. Gerçek τ_b(3B) kritik noktada oturur. 3B bağlantılılık eşiği o kadar düşük olduğundan (≈ 1,55), z̄(τ) ona ancak en yukarıda ulaşır: τ_b(3B) = 3,18 (YMK) ya da 2,85 (HMK) — τ_c = z_max/4'ün (3,00 ya da 2,00) üstünde ya da hemen üstünde. Demek ki bağ dev öbeği bütün sıvı fazı atlatır ve ancak kritik noktada çözülür. Bu, tam olarak Böl. 7.10'daki öngörüdür ("gerçek τ_b(3B) 2B değerinin üstünde, kritik noktaya çok yakın olurdu") ve "kaynama = bağlantılılık çözülmesi" resmini bağımsız biçimde doğrular: 3B'de madde kritik noktaya çok yaklaşana dek sürekli bağlı kalır.

Her iki bulgu Böl. 5.4a'daki eşik ayrımını pekiştirir: erime ↔ rijitlik (düşük, τ_m ≈ τ_rig), kaynama ↔ bağlantılılık (yüksek, τ_b(3B) ≈ τ_c). Erime yoğunluğu bu yüzden rijitlikten gelir (v_solid = 0,667), bağlantılılıktan değil — iki eşik 3B'de sıcaklık ekseninde modelde hiç olmadığı kadar birbirinden uzaklaşmıştır.

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

7.12 (τ, v) yüzeyinde katı(1)→katı(2) dönüşüm çizgisi

Böl. 6.3, YMK↔HMK örgü dönüşümünü bir nokta τ_ss olarak vermişti (p = 0'da). Basınç üzerinden bundan bir çizgi olur (polymorph_line_16.py) ve bu çizgi Böl. 7.5/7.9'daki ortak (τ, v) yüzeyine dördüncü sınır çizgisi olarak aittir — birlikte varoluş kubbesi (kaynama), erime çizgisi ve sert duvarın yanına.

Hacim sıçraması. İki katı da neredeyse sabit yoğunlukta oturur; hacim farkları paketleme yoğunluklarından çıkar (v ∝ 1/η): η_YMK = 0,7405, η_HMK = 0,6802. v_YMK = 0,667 (Böl. 7.9) ile v_HMK = v_YMK·η_YMK/η_HMK = 0,726 çıkar, yani

Δv = v_HMK − v_YMK = +0,059 (+8,9 %).

Daha açık olan HMK'nin hacmi daha büyüktür — YMK→HMK yeniden kristalleşmesi bir hacim artışıyla bağlıdır (çizgi romandaki "Yeniden kristalleşme" ile aynı sav).

Clausius–Clapeyron. ΔS = S_HMK − S_YMK ile (yüksek eksi düşük sıcaklık fazı) eğim şudur:

dp/dτ = ΔS/Δv.

Her iki büyüklük de pozitiftir (HMK yüksek sıcaklık fazıdır, yani daha çok entropisi ve daha büyük hacmi vardır), dolayısıyla dp/dτ > 0: daha yüksek basınç dönüşüm sıcaklığını yükseltir, çünkü daha yoğun YMK'yi kayırır. θ_HMK/θ_YMK = 0,60 referansı için (τ_ss = 1,333, ΔS = 0,50 k_B, L = 0,667 ε):

Referans için (τ, p)'de çizgi (yükselen): p = 0 → τ_ss = 1,33; p = 2 → 1,57; p = 4 → 1,81; p = 6 → 2,04 (indirgenmiş birimler).

(τ, v) yüzeyinde dönüşüm neredeyse dikey iki katı yoğunluk şeridi olarak görünür — katı(1) = YMK v ≈ 0,667'de (daha yoğun, düşük T) ve katı(2) = HMK v ≈ 0,726'da (daha açık, yüksek T) —, bunlar τ_ss'de birbirine döner; üstünde τ_m'de erime gelir. τ ekseni üzerindeki tam sıra böylece τ_ss < τ_m < τ_c ≤ τ_b(3B)'dir (Böl. 7.11).

Gerçek birimler. Erime noktası kalibrasyonu ε/k_B = T_m/τ_m (τ_m = 1,70, 3B YMK) üzerinden c_v basamağı ve dönüşüm ısısı elle tutulur olur. c_v basamağı kalibrasyondan bağımsızdır: |Δc_v| = 1,0 k_B = 1 R = 8,31 J/(mol·K) = 0,0083 kJ/(mol·K). Dönüşüm ısısı ΔH_ss = L·R·(ε/k_B) (referans L = 0,667 ε):

Dönüşüm ısıları küçük ve pozitiftir (0,2–0,9 kJ/mol mertebesi) ve böylece erime entalpilerinin (Böl. 7.6, birkaç kJ/mol) belirgin biçimde altında kalır — zayıf bir polimorfizm dönüşümünün küçük c_v basamağıyla beklenen resmi tam da budur. (θ oranı 0,60 bilinçli olarak sıkı ideal resmin dışındadır; Böl. 6.6, salt geometrinin θ ≈ 0,82–0,94 vereceğini ve dönüşümün o zaman kapalı kalacağını gösterir — çizgi, yumuşak bir HMK fazının eğer-ise sonucudur.)

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

7.13 Demire kalibre edilmiş (τ, p) diyagramında katı(1)→katı(2) dönüşüm çizgisi

Böl. 7.12 dönüşümü bir (τ, v) çizgisi olarak verir; sürüm 16 onu (τ, p) resminde bir çizgi olarak tamamlar ve gerçek bir YMK↔HMK polimorfuna kalibre eder (solidsolid_pt_17.py).

Sert çekirdekli van der Waals hâl denkleminden basınç alanı. İki katı yoğunluğunda v_YMK = 0,667 ve v_HMK = 0,726 (Böl. 7.12) indirgenmiş hâl denklemi p_r = 8τ_r/(3v_r − 1) − 3/v_r² (Böl. 7.3) baştan sona pozitif basınçlar verir; aynı τ_r'de daha yoğun YMK daha yüksek basıncı taşır (τ_r = 1: p_r = 1,25'e karşı 1,10). Sıkıştırma demek ki YMK'yi kayırır — pozitif Clausius–Clapeyron eğimi dp/dτ = ΔS/Δv > 0 ile tutarlı.

Çizgi ve dönüşüm noktası. θ_HMK/θ_YMK = 0,60 referansı için (τ_ss = 1,333, ΔS = 0,50 k_B, dp\/dτ = 8,47; indirgenmiş birimlerde p\ = p·D³/ε) çizgi p\(τ) = 8,47·(τ − 1,333) yükselir: p\ = 0 → τ = 1,33; p\ = 2 → 1,57; p\ = 4 → 1,81; p\ = 6 → 2,04. **p\_tri ≈ 3,1'de katı-katı çizgisi erime işareti τ_m = 1,70'e çarpar — (τ, p) düzleminde bir katı(1)-katı(2)-sıvı dönüşüm noktası** (üçlü noktanın karşılığı).

Demire kalibrasyon (ε/k_B = T_m/τ_m = 1811 K / 1,70 = 1065 K). Demirin ters işaretli iki YMK↔HMK geçişi vardır — tam da sınama durumu:

γ→δ için (τ_ss = 1667/1065 = 1,565 → θ_HMK/θ_YMK ≈ 0,61) model işareti ve büyüklük mertebesini tutturur: ΔH_ss ≈ 4,2 kJ/mol (gerçek ≈ 0,84 kJ/mol — mertebe kJ/mol, keskin eşik yaklaşımının değeri ~5 kat yüksek), c_v basamağı parametresiz 1 R = 8,31 J/(mol·K), dT/dp ≈ +245 K/GPa. Buna karşılık α→γ geçişinde yoğun faz (γ) yüksek sıcaklık fazıdır — bunu salt paketleme geometrisi veremez (daha yoğun = daha düşük enerjili = düşük T); α-demir manyetik olarak kararlı kılınır. Demir böylece dürüst bulgunun her iki yanını da doğrular: geometrik topoloji (γ→δ) ideal resmin içindedir, manyetik anomali (α→γ) dışındadır — Böl. 6.6/6.8'in θ ≈ 0,87…1,00 > 0,61 ile nicelediği aynı "dışarısı".

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

7.14 Ayrı bir ideal resim dalı olarak Redlich–Kwong

van der Waals dalının Guldberg açıklığı (Böl. 7.4, ≈ 0,60 yerine T_b/T_c ≈ 0,48) orada yalnız adlandırılmıştı. Açık istek üzerine geliştirilmiş bir sürüm, Redlich–Kwong kohezyonunu kendi kalibrasyonu, kendi birlikte varoluşu ve kendi hata bölümüyle bağımsız bir ideal resim türü olarak sokar (rk_branch_17.py) — van der Waals dalı varsayılan olarak kalır; modelin geri kalan bütün fiziği (erime, perkolasyon, katı(1)→katı(2)) kohezyon hâl denklemi seçiminden etkilenmez.

Kendi enerji ölçeği. a_RK = Ω_a·R²T_c^{2,5}/p_c ve b_RK = Ω_b·R T_c/p_c'den dalın kohezyon enerjisi şöyle çıkar:

ε_lg,RK/k_B = (Ω_a/Ω_b)·T_c = 4,93·T_c (vdW ε_lg/k_B = 27/8·T_c = 3,375·T_c'ye karşı),

daha büyüktür, çünkü sıcaklığa bağlı kohezyon a/√T √T_c çarpanını içeri katlar. b_RK'den çıkan çap gerçek değerleri van der Waals dalındaki kadar iyi tutturur (argon D ≈ 3,34·10⁻¹⁰ m).

Kendi birlikte varoluşu (indirgenmiş, kritik nokta (1, 1, 1), Z_c = 1/3; p_r = 3τ_r/(v_r − β) − A/(√τ_r·v_r(v_r + β)), β = 0,25992, A = 3,84732):

Kendi hata bölümü (p_r,atm = p_atm/p_c'de kaynama noktası, T_b/T_c = τ_r,b):

Sonuç. RK dalı T_b/T_c'yi vdW ≈ 0,48'den ≈ 0,60'a yükseltir ve böylece Guldberg sabitini ve gerçek değerleri pratikte tutturur (ortalama hata vdW'nin −%21'ine karşı −%1). O, van der Waals varsayılanının yanında ikinci, açıkça ayrılmış bir ideal resim türü olarak durur; hangi dalın "daha doğru" olduğu, model yapısının değil istenen duyarlığın sorusudur — Model 3'ün çekirdeği (düzen parametresi olarak z̄) her iki dalda da aynıdır.

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

7.15 Başka polimorflarda katı(1)→katı(2): titanyum ve soy gaz kontrolü

Böl. 7.13 katı(1)→katı(2) (τ, p) çizgisini demire kalibre eder. Geliştirilmiş sürüm işaret ve büyüklük mertebesini ikinci, manyetik olarak dikkat çekmeyen bir polimorfta güvenceye alır ve bir negatif kontrol ekler (solidsolid_titanium_18.py). Kalibrasyon ve hacim sıçraması Böl. 7.13'teki gibidir (ε/k_B = T_m/τ_m, τ_m = 1,70; paketleme geometrisinden Δv = +8,9 % yoğun→açık); hsp, YMK gibi z = 12 üzerinden yaklaştırılır (ymk/hsp farkı yalnız adlandırılır, çözülmez).

Titanyum 1155 K'de α(hsp) → β(hmk) ile tam model topolojisine sahiptir: yoğun faz (hsp, z = 12) düşük sıcaklık fazıdır, ısıtınca açık olur (hmk, z = 8), yani ΔS > 0 ve Δv > 0 → dp/dτ > 0 — demir γ→δ ile aynı işaret, ama manyetik anomali olmadan. Karşılaştırmada iki bağımsız yoğun→açık polimorf:

Her ikisi de işareti (dp/dτ > 0, basınç yoğun fazı kayırır) ve büyüklük mertebesini tutturur (ΔH ~ kJ/mol, c_v basamağı parametresiz 1 R = 8,31 J/(mol·K), dT/dp ~ 10²·K/GPa). ΔH_ss'nin değeri her ikisinde de dizgesel olarak ~bir çarpan kadar yüksektir — μ karşılaştırmasının keskin eşik yaklaşımı gizli dönüşüm ısısını abartır; bu, adlandırılmış bir model sınırıdır, işaret hatası değil.

Argon negatif kontroldür: katı argon YMK'dir ve erimeye dek (83,8 K) YMK kalır — normal basınçta kararlı bir hmk yoktur. Salt çift potansiyelli ideal resmin öngördüğü tam da budur: θ_HMK/θ_YMK ≈ 0,9…1,0 > 0,61 ve kaymaya kararsız bir LJ HMK'si (Böl. 6.9) ⇒ katı(1)→katı(2) geçişi yok. Karşı sınama bulgunun yapısını doğrular: ideal resim bir katı(1)→katı(2) çizgisini ancak geometrik olmayan katkılar (demir ve titanyumdaki gibi elektronik/manyetik) HMK'yi genel olarak kararlı kıldığında üretir. Titanyum, Böl. 6.8/6.9'daki aynı dürüst çerçevenin pozitif, argon ise negatif kanıtını verir.

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

8 — Yazın karşılaştırması: statik termodinamiğin model tasavvurları

Özet. Bölüm 8, Model 3'ü statik termodinamiğin yerleşik model tasavvurları arasına yerleştirir:

• 8.1 Genel bakış tablosu: on bir model (ideal gaz, van der Waals, Redlich–Kwong, virial açılım, Einstein, Debye, Lennard-Jones/sert küreler, hücre kuramı, Ising/örgü gazı, perkolasyon, Flory-Huggins) faz(lar)ı ve istatistiksel-termodinamik temel fikriyle.
• 8.2 Modeller tek tek: her biri için özellikler, uygulama ve sınırlar içeren birer paragraf, kaynak bağlantısıyla başlayarak; serbest erişilebilir kaynaklar model3/links/ altında arşivlenmiştir, hücre kuramı (ödeme duvarı) yalnız anılmıştır.
• 8.3 Karşılaştırma: Model 3, düzen parametresi z̄ olan bir işbirlikçi iki-durumlu bağ modelidir; en yakın akrabası Ising/örgü gazıdır, kaynama doğrudan perkolasyondan alınmıştır, fonon resmi (Einstein/Debye) erime sertliğini taşır. Model 3 bu üç gelenek çizgisini ortak bir ε ölçeğinde kenetler.

Bu bölüm Model 3'ü katıların, sıvıların ve gazların denge ("statik") termodinamiğinin yerleşik model resimleri arasına yerleştirir. Önce tablo biçiminde bir genel bakış (8.1), sonra her model için tipik özellikler, uygulama alanı ve sınırlar içeren kısa bir düz metin paragrafı — her paragraf kaynağın bağlantısıyla başlar (8.2). Sonunda Model 3 ile madde madde bir karşılaştırma (8.3). Serbest erişilebilir kaynaklar model3/links/ altında arşivlenmiştir; sınama betiği model3/scripts/check_links_16.py alıntı ile arşivi SHA-256 üzerinden tutarlı tutar.

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

8.1 Genel bakış tablosu

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

8.2 Modeller tek tek

https://de.wikipedia.org/wiki/Kinetische_Gastheorie — İdeal gaz / kinetik gaz kuramı. Öz hacimsiz ve esnek çarpışmalar dışında etkileşimsiz parçacıklar; basınç salt momentum aktarımından doğar. pV = Nk_BT hâl denklemini, Maxwell-Boltzmann hız dağılımını ve ⟨E⟩ = (3/2)k_BT'yi verir. Uygulama: yüksek sıcaklık ve düşük basınçta seyreltik gazlar. Sınırlar: yoğuşma yok, sıvı ya da katı faz yok, faz geçişi yok — etkileşim yapı gereği yoktur.

https://de.wikipedia.org/wiki/Van-der-Waals-Gleichung — van der Waals denklemi. (p + a/v²)(v − b) = k_BT, ideal gaza bir dışlama hacmi b (itme) ve bir kohezyon terimi a (çekim) ekler. Bu, gaz-sıvı geçişi, kritik nokta ve Maxwell kurgusu içeren ilk hâl denklemidir. Uygulama: yoğuşmanın niteliksel resmi, kritik büyüklükler. Sınırlar: kritik noktaya yakın niceliksel olarak kesin değil (ortalama alan üsleri), katı faz yok.

https://de.wikipedia.org/wiki/Redlich-Kwong-Gleichung — Redlich-Kwong hâl denklemi. Sıcaklığa bağlı çekim terimi a/[√T · v(v + b)] içeren kübik hâl denklemi. √T ölçeklemesi buhar-sıvı dengelerini van der Waals'a göre belirgin biçimde iyileştirir (bkz. Böl. 7.4; orada bu, T_b/T_c oranını yükseltir). Uygulama: teknik gaz fazı ve buharlaşma hesapları. Sınırlar: görgül olarak parametrelendirilmiş, hâlâ katı faz yok; kritik nokta yakınında yine ortalama alan türü.

https://de.wikipedia.org/wiki/Virialentwicklung — Virial açılım. Dizgesel açılım pV/Nk_BT = 1 + B(T)/v + C(T)/v² + … Virial katsayılar moleküller arası potansiyelden çıkar (istatistiksel mekaniğin öbek integralleri) ve böylece mikroskobik olarak temellendirilmiştir. Uygulama: orta yoğunluğa dek gazlar, başka hâl denklemleri için referans. Sınırlar: seri yoğun sıvılar/katılar için yakınsamaz; yüksek katsayılar hesapça ağırdır.

https://de.wikipedia.org/wiki/Einstein-Modell — Einstein katısı. Kristal, aynı ω frekanslı 3N bağımsız harmonik salınıcıdan oluşan bir topluluktur. Isı sığası c_v'nin θ_E = ħω/k_B karakteristik sıcaklığının altındaki düşüşünü ve yüksek sıcaklıkta Dulong-Petit sınır değeri 3Nk_B'yi açıklar. Uygulama: katı ısı sığalarının en yalın betimi. Sınırlar: yanlış düşük sıcaklık yasası (T³ yerine üstel), akustik kipler yok, yalnız katı faz.

https://de.wikipedia.org/wiki/Debye-Modell — Debye modeli. Katıyı, Debye frekansına dek fonon spektrumu olan esnek bir süreklilik olarak ele alır; düşük sıcaklıkta c_v ∝ T³, yüksekte Dulong-Petit — Debye sıcaklığı θ_D örgü sertliğini özetler. Uygulama: basit kristallerin bütün T aralığında niceliksel olarak iyi ısı sığaları. Sınırlar: süreklilik yaklaşımı (gerçek örgülerin dağınımını yok sayar), yalnız katı faz, erime yok.

https://de.wikipedia.org/wiki/Lennard-Jones-Potential — Lennard-Jones potansiyeli / sert küreler. İtme çekirdeği ve zayıf çekim içeren çift potansiyeli 4ε[(σ/r)¹² − (σ/r)⁶]; sert küreler sınır modeli yalnız itmeyi tutar. Moleküler dinamik ve Monte Carlo'nun temelidir ve — sayısal olarak — üç fazı da üçlü noktayla birlikte betimler. Uygulama: soy gazlar, genel referans dizgesi. Sınırlar: ikili toplanabilirlik, çekirdekte yalnız benzetim ya da tedirginlik kuramıyla değerlendirilebilir (kapalı bir hâl denklemi yok).

https://doi.org/10.1098/rspa.1937.0210 — Sıvının hücre kuramı (Lennard-Jones-Devonshire). Her molekül, komşularının ortalanmış alanda oluşturduğu bir hücrede (kafeste) devinir; sıvı, hücre durum toplamıyla betimlenir ve bir sıvı hâl denklemi ile bir erime yaklaşımı verir. Uygulama: basit sıvılar, katı ile gaz arasında bağlantı halkası. Sınırlar: düzeni fazla tahmin eder, boşlukları/delikleri ihmal eder (bunu sonraki delik kuramları giderir). (Ödeme duvarı — yalnız anılmıştır, yerel olarak arşivlenmemiştir; erişilebilirliği sınama betiğinde sınanır.)

https://de.wikipedia.org/wiki/Ising-Modell — Ising modeli / örgü gazı. Bir örgüde yerler, en yakın komşu bağlaşımıyla iki-durumlu bir değişken (spin ↑/↓ ya da dolu/boş) taşır. Örgü gazı matematiksel olarak Ising modeline eşyapılıdır ve gaz↔sıvı yoğuşmasını kritik noktalı işbirlikçi bir düzen-düzensizlik geçişi olarak betimler. Uygulama: kritik olgular, işbirlikçi geçişler. Sınırlar: örgü kesiklileştirmesi, açık kinetik enerji yok; ortalama alan yaklaşımında yalnız niteliksel.

https://de.wikipedia.org/wiki/Perkolationstheorie — Perkolasyon kuramı. Salt geometrik model: yerler ya da bağlar p olasılığıyla rastgele doludur; bir p_c eşiğinin üstünde boydan boya bir öbek vardır. Jelleşmeyi, iletkenliğin başlamasını ve genel olarak bağlantılılık geçişlerini betimler. Uygulama: ağ oluşumu, bağlılığın çözülmesi/kurulması. Sınırlar: kendinden enerji/sıcaklık yok — termodinamiğe bağlaşım dışarıdan eklenmelidir.

https://de.wikipedia.org/wiki/Flory-Huggins-Theorie — Flory-Huggins kuramı. Çözeltiler ve polimerler için örgü modeli: zincir yerleşiminin karışma entropisi artı bir etkileşim parametresi χ, karışabilirliği ve ayrışmayı yönetir. Uygulama: polimer çözeltileri, karışımların faz diyagramları. Sınırlar: ortalama alan, sıkıştırılamaz örgü, χ çoğunlukla görgül — tek bileşenli katı/sıvı/gaz resimleri için ancak sınırlı ölçüde ilgili, burada örgü referansı olarak anılmıştır.

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

8.3 Model 3 ile benzerlikler ve farklar

Model 3 bir işbirlikçi iki-durumlu bağ modelidir: düzen parametresi ortalama koordinasyon sayısı z̄ = z_max·x'tir; kırık durum dejenereliği G = z_max ile işbirlikçilik K = J·z_max üzerinden kendiyle tutarlıdır; erime bir sertlik (fonon) eşiği olarak, kaynama kritik noktada bağlantılılık (perkolasyon) çözülmesi olarak, ln G > 2 + 4ε/K olur olmaz birinci mertebeden geçiş.

Kısaca: Model 3, başka türlü ayrı duran üç gelenek çizgisini — örgü gazı/Ising resmi (işbirlikçi iki-durumlu geçiş), perkolasyon resmi (bağlantılılık → kaynama) ve Einstein/Debye'nin fonon resmi (sertlik → erime) — kenetler ve bunları tek bir enerji parametresi ε üzerinden ortak bir τ ölçeğinde kalibre eder. van der Waals ve hücre kuramı, kesikli bağ resminin kendini karşılaştırabileceği sürekli karşı sınamayı sağlar (sırasıyla hâl denklemi ve sıvı kafesi).

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

9 — Ek

Özet. Ek, biçimsel kanıtları toplar:

• 9.1 Kaynakça: yazar adlarına göre sıralı, serbestçe arşivlenmiş kaynaklara bağlantılarla (check_links_16.py ile SHA-256 karşılaştırması; hücre kuramı ödeme duvarı nedeniyle yalnız anılmıştır).
• 9.2 YZ desteği: her önemli değişiklik için YZ'nin kendi kestirdiği YZ/yazar yenilik payı (yazar payının strateji/yapı/redaksiyon olarak ayrıştırılması ve bir kontrol hesabıyla) ve tek tek adımların maliyet (kestirim), süre ve önceki günlüğe uzaklıkla dökümü (toplam ≈ 160 USD, açıkça işaretlenmiş kestirimler).
• 9.3 Kısaltma ve simge sözlüğü: sık kullanılan simgeler (τ, ε, z̄, K, G, θ, …); bunlar sayesinde istisnai olarak metinde önlerinde açıklama olmaksızın geçebilirler.
• 9.4 Hesaplama betikleri ve test kapsamı: bütün sayılar yeniden üretilebilir Python betiklerinden gelir (YZ programlama payı %100, test kapsamı her yerde %0), ayrıca yazarın kodu ve sonuçları kendisi denetlemediğine dair açık bir "utandırıcı bildirim".

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

9.1 Kaynakça

Yazar adlarına göre sıralı. Serbest erişilebilir kaynaklar model3/links/ altında arşivlenmiştir; bunların anılan bağlantılarla uyumunu model3/scripts/check_links_16.py (SHA-256) sınar. Lennard-Jones ve Devonshire'ın hücre kuramı ödeme duvarı nedeniyle yalnız anılmıştır (erişilebilirlik sınama betiğinde sınanır).

• Broadbent, S. R.; Hammersley, J. M. (1957): Percolation processes I. Crystals and mazes. Proc. Cambridge Phil. Soc. 53, 629–641. — Genel bakış/arşiv: Perkolasyon kuramı
• Clausius, R.; Maxwell, J. C.; Boltzmann, L. (1857–1872): kinetik gaz kuramının temellendirilmesi. — Genel bakış/arşiv: Kinetik gaz kuramı
• Debye, P. (1912): Zur Theorie der spezifischen Wärmen. Ann. Phys. 39, 789–839. — Genel bakış/arşiv: Debye modeli
• Einstein, A. (1907): Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Ann. Phys. 22, 180–190. — Genel bakış/arşiv: Einstein modeli
• Flory, P. J.; Huggins, M. L. (1941–1942): polimer çözeltilerinin örgü kuramı (χ parametresi). — Genel bakış/arşiv: Flory-Huggins kuramı
• Ising, E.; Lenz, W. (1920–1925): Beitrag zur Theorie des Ferromagnetismus. Z. Phys. 31, 253–258 (örgü gazı matematiksel olarak eşyapılı). — Genel bakış/arşiv: Ising modeli
• Kamerlingh Onnes, H. (1901): hâl denkleminin virial biçimi. — Genel bakış/arşiv: Virial açılım
• Lennard-Jones, J. E. (1924): On the determination of molecular fields. Proc. R. Soc. Lond. A 106, 463–477. — Genel bakış/arşiv: Lennard-Jones potansiyeli
• 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 (Ödeme duvarı — yalnız anılmıştır, yerel olarak arşivlenmemiştir.)
• Redlich, O.; Kwong, J. N. S. (1949): On the thermodynamics of solutions. V. An equation of state. Chem. Rev. 44, 233–244. — Genel bakış/arşiv: Redlich-Kwong denklemi
• van der Waals, J. D. (1873): Over de Continuïteit van den Gas- en Vloeistoftoestand (doktora tezi, Leiden). — Genel bakış/arşiv: van der Waals denklemi

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

9.2 YZ desteği

Bu model geliştirmesi bir insan-YZ birlikteliğinde doğdu: yazar (YZ kullanıcısı) ne üzerinde çalışılacağını belirler (yol gösterici sorular, kararlar, model istekleri); YZ (Claude Code) büyük ölçüde nasıl olacağını belirler — türetimler, betikler, yapı ve metin. Aşağıdaki iki tabloyu YZ kendisi kestirmektedir.

(a) Önemli içerik değişiklikleri ve kestirilen yenilik payı. Yenilik payı yol gösterici fikirleri %45, yapılandırıcı işleri %40, redaksiyon işlerini %15 ağırlıklandırır. Yol gösterici fikirler (ne üzerinde düzlemi) ağırlıkla yazardan, yapılandırıcı ve redaksiyon işleri (nasıl düzlemi) ağırlıkla YZ'den gelir; YZ'nin fiziksel yol gösterici fikri de kattığı yerde (örn. kaynamanın perkolasyon okuması, θ yay modeli) onun payı artar. Değerler yüzde, YZ + yazar = 100. Üç ek sütun yazar payını üç etki düzlemine ayırır — strateji (yol gösterici fikirler), yapı (türetimler/betikler/kuruluş) ve redaksiyon (metin işi). Bunların %45 / %40 / %15 ağırlıklı ortalaması Yazar % sütununu bir kontrol hesabı olarak yeniden üretir (model3/scripts/author_share_19.py, en büyük sapma 0,95 yüzde puanı).

(b) Tek tek adımlar: maliyet (kestirim), süre, önceki günlüğe uzaklık. Oluşturma tarihi, günlük dosyasının oluşturulma zamanıdır (dosya mtime'ı). Önceki günlüğe Δ, bir önceki günlük girdisine olan zaman aralığıdır (adımlar arasındaki kullanıcı emeğinin ölçüsü, aralar//compact dâhil). Maliyet (kestirim), her günlükte belgelenen kestirimlerdir (Anthropic liste fiyatları, ± %30). Süre (etkin), salt işlem süresine dair bir YZ kestirimidir; /compact ile ayrılmış adımlarda Δ'dan büyük olabilir. Bir yıldız, var olan içeriğin yeniden kurulduğu/elden geçirildiği tahmin edilen adımı işaretler.

Belgelenen maliyetlerin toplamı (kestirim, adımlar 01–21): ≈ 129 USD (tek tek kestirimlerin toplamı);

1. adımla (16d, ≈ 10 USD) ≈ 139 USD, 23. adımla (S16, ≈ 12 USD) ≈ 151 USD ve 24. adımla (S17, ≈ 9 USD)

≈ 160 USD. Bütün değerler açıkça işaretlenmiş kestirimlerdir (ccusage ölçümü yoktur); günlük ayrıntıları model3/logs/dev_log_NN.md içindedir.

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

9.3 Kısaltma ve simge sözlüğü

Bu metin ilkesel olarak her simgenin önüne betimleyici bir kavram koyar ("c_v" yerine "ısı sığası c_v"). Aşağıdaki sık kullanılan büyüklükler bunun istisnasıdır: burada bir kez ve hep açıklandıkları için düz metinde önlerinde betimleme olmaksızın geçebilirler. Bütün τ… sıcaklıkları indirgenmiştir (bağ enerjisi ε birimlerinde).

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

9.4 Hesaplama betikleri ve test kapsamı

Metinde anılan bütün sayılar ve tablolar model3/scripts/ içindeki yeniden üretilebilir Python betiklerinden gelir (her biri kendi ham çıktısı …_output.txt ile). YZ programlama payı betik kodunun ne kadarını YZ'nin yazdığını gösterir; test kapsamı otomatik testleri (birim testler, assert denetimleri) anlatır. Bu her yerde düşüktür: betikler test iskelesi olmayan tek seferlik hesap koşularıdır; tek güvence ham çıktının akla yatkınlığıdır.

Yazar açıkça itiraf ediyor (utandırıcı bildirim): bu betiklerin kaynak kodunu da, onlarla üretilen sayıları, tabloları ve eğrileri de ne denetledi, ne yeniden hesapladı, ne de içerikçe doğruladı. Burada YZ'nin yetkinliğine körü körüne güvendi, çünkü ön planda yalnızca model tasavvurunun niteliksel karakteri ve akla yatkınlığı vardı. Bu yüzden değerler ihtiyatla alınmalıdır (ayrıca bkz. Böl. 10.4'teki geniş uyarı).

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

10 Özeleştiri

Bölüm 9.2 (a)'daki tabloda YZ'nin kestirdiği yazar emeği payının üç alanda da neredeyse hep en az 10–20 yüzde puanı fazla olduğu izlenimindeyim. Yalnızca temel fikri, bir iki yorumu ve grafik görselleştirmedeki önemli incelikleri kattığım izlenimindeyim.

(Matematiksel) 4–9. bölümleri YZ önerdi, araştırdı ve gerçekleştirdi. Daha kötüsü: geliştirme sırasında YZ'nin koordinasyon sayısına dayalı istatistiği geliştirmeye saplandığı izlenimindeydim. Asıl temel fikirlerimi — yalnız atomların değil moleküllerin de bağ kurabileceğini ve sıvı içindeki bağ kopmalarının fraktal düzenliliklere uyabileceğini — YZ ancak güçlükle alabildi ve terminolojiyi yok saymaya çalıştı. Bu yüzden bu temel fikirler metinde hâlâ çok ikincil ele alınıyor. Bunu sonradan yeniden yazmadım, çünkü temel aksiyomlardan yola çıkarak faz geçişlerini düşünürken bu fikirler zaten akılda tutuluyor. Ayrıca fikrin bir ölçü olarak koordinasyon sayısına saplanması/indirgenmesi, resimli tasavvurumu matematikleştirmek için en kötü yol değil; çünkü bu, sonradan bakıldığında akla yatkınlık sınaması için epeyce verimli çıktı — gerçi YZ'nin sonuçlarının akla yatkınlığını yalnızca çok kabaca sınadım ve Python kodunun hiçbir denetimini yapmadım. Mazeretim olarak şu söylensin: benim için hep yalnızca farklı maddenin hâllerindeki moleküler ilişkilerin resimli tasavvuru söz konusuydu. Bu tasavvuru matematiksel olarak doğrulamak isteyen YZ'nin pohpohlamasına hayır demedim. YZ'nin yönlendirmesi bende demek ki etkili oldu, çünkü anılan bölümleri en üstünkörü biçimde sınadım ve inanarak kabul ettim. (Kendime dair bu deneyim, YZ denetimi için "human in the loop"un gerçekten etkili olup olmadığına dair güvensizliğimi güçlendiriyor.)

Bu metin YZ desteğiyle oluşmuştur (Claude Code 2.1.161, model Opus 4.8 / claude-opus-4-8[1m], 1M bağlam; tarih 2026-09-05). — YZ desteğiyle oluşmuştur.

Filigran hakkında not. Almanca kaynak metin (de.modell_final_02.md) anahtar sözcüklerde Yunan homogliflerinden bir filigran taşır (Α→A, Β→B, Ε→E, Η→H, Ι→I, Κ→K, Μ→M, Ν→N, Ο→O, Ρ→P, Τ→T, Υ→Y, Χ→X, Ζ→Z). Bu hüner Almanca yazımlara bağlı olduğundan, bu çeviri anahtar sözcükleri olağan büyük harflerle taşır: SERBEST UÇUŞ (FREIFLUG), ÇOKLUK (VIELHEIT), BAĞ (BINDUNG), ÜÇLÜ AKOR (DREIKLANG), SERT ÇEKİRDEK (HARTKERN), KORUNUM (ERHALTUNG), ÇATLAK (RISS), ÇEVİK (FLINK).