Thermodynamique statistique du corps pur idéal et de ses états de la matière

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

Résumé général

Le modèle 3 est un modèle qualitatif, de mécanique statistique, d'une substance idéale qui dérive les trois états de la matière — solide, liquide et gazeux — ainsi que leurs transitions à partir d'une seule image de liaison, celle des liaisons moléculaires. Le principe de base réside dans la prise en compte de trois réservoirs d'énergie pour l'énergie interne : l'énergie de vibration dans les liaisons moléculaires, l'énergie de translation des molécules et l'énergie cinétique coopérative dans la rotation des micelles, forme idéalisée des flux internes de matière qui compensent les fluctuations internes de densité. Une grandeur fondamentale importante dans ce contexte est le nombre de coordination moyen z̄ = z_max·x, qui s'établit de façon auto-cohérente par une liaison coopérative à deux états (de type Ising, coopérativité K = J·z_max, dégénérescence de l'état rompu G = z_max). Il en découle trois transitions sur une échelle τ commune, calibrée par un seul paramètre d'énergie ε : la fusion comme seuil de phonons/rigidité (z̄ atteint la rigidité), l'ébullition comme effondrement de la connectivité (percolation) au point critique, et une transition facultative solide(1)→solide(2) (CFC↔CC) avec le palier de c_v sans paramètre de 1,0 k_B. Une transition devient du premier ordre dès que ln G > 2 + 4ε/K. Ici la dégénérescence de l'état rompu G est le nombre de contacts de voisinage susceptibles d'être rompus (G = z_max, le nombre de coordination maximal), l'énergie de liaison ε est l'énergie gagnée par contact de voisinage, et la coopérativité K = J·z_max est le produit de la force de couplage J, de type Ising, par z_max ; la transition devient donc d'autant plus vite du premier ordre que la multiplicité de contacts G est grande et que la coopérativité K est forte face à l'énergie de liaison ε.

L'article conduit des sept axiomes (chapitre 1) à l'analyse des contradictions (chapitre 2) puis à une description imagée avec des bandes dessinées statiques (chapitre 3), pour arriver à la charpente de mécanique statistique — liaison coopérative, fonctions de partition et capacités thermiques (chapitres 4–5) —, à la plausibilité mathématique en 3D (chapitre 6) et enfin à la calibration et à l'équation d'état (chapitre 7). Les idées sont complétées par le balayage 3D complet de c_v pour la fusion du premier ordre (Δz̄, chaleur latente ΔU, décalage du pic), un modèle de ressorts d'Einstein pour la rigidité phononique θ(réseau) et la ligne de transformation solide(1)→solide(2) via Clausius–Clapeyron. θ(réseau) est ensuite affiné en une densité d'états de phonons de type Debye ; la transformation solide(1)→solide(2) est en outre introduite comme ligne dans le plan (τ, p) (calibrée sur le fer) et une branche Redlich–Kwong séparée est ajoutée à côté du choix par défaut van der Waals. θ(réseau) est par ailleurs éprouvé sur un réseau de Lennard-Jones complet (CFC stable, CC idéal instable au cisaillement), et la topologie solide(1)→solide(2) est éprouvée sur le titane (deuxième confirmation) et l'argon (contrôle négatif). Le chapitre 8 situe le modèle parmi onze conceptions établies de la thermodynamique statique. L'annexe documente les sources bibliographiques (classées par auteur) et l'assistance de l'IA (part d'innovation estimée, coût estimé et durée par étape de développement) ; les aides rassemblent le mode de lecture, les points critiques restants et les perspectives.

Énoncé central. Une seule mesure coopérative de la liaison, z̄, agrafe ensemble trois traditions par ailleurs séparées — l'image du gaz sur réseau / d'Ising (transition coopérative à deux états), l'image de la percolation (connectivité → ébullition) et l'image des phonons d'Einstein/Debye (rigidité → fusion) — en un tableau des phases cohérent et qualitativement complet, délibérément à l'intérieur d'une image idéale de type van der Waals, l'édifice théorique pouvant naturellement se déduire d'une image dynamique aisément imaginable. Ce n'est pas l'exactitude mais le caractère imagé qui fait la force propre du modèle et qui caractérise sa valeur scientifique en tant que modèle.

Table des matières

• Résumé général
• 1 — Les axiomes et leurs mots-clés (consolidés)
• 2 — Analyse des contradictions et des lacunes
• 3 — Descriptions imagées
• 3.1 Les trois phases isolées
• 3.2 Transitions de phase
• 3.3 Pourquoi la vitesse du son dans les liquides décroît avec T
• 3.4 Comment naît le mouvement brownien — avec une notion affinée de la micelle
• 3.5 Pourquoi les réactions aux couches limites des micelles s'accélèrent
• 3.6 La fenêtre liquide tridimensionnelle
• 4 — Modèle de mécanique statistique
• 4.1 Degrés de liberté et variables
• 4.2 Liaison coopérative à deux états (de type Ising, paramètre J)
• 4.3 Bilan énergétique — avec terme coopératif
• 5 — Fonctions de partition et capacités thermiques
• 5.1 Fonction de partition factorisée
• 5.2 Capacité thermique par phase (analytique)
• 5.3 Courbe numérique c_v(τ) **avec** terme coopératif — chaleur latente quantifiée
• 5.4 Seuil de percolation de connectivité comme résultat (confirmé inchangé) - 5.4.a Deux seuils distincts nettement séparés
• 5.5 Ordre de la transition de fusion
• 5.6 Étude de taille finie du seuil de percolation
• 5.7 Balayage 3D complet de c_v(τ) pour la fusion du premier ordre
• 6 — Plausibilité mathématique de l'interprétation 3D
• 6.1 Changement du degré de liberté de rotation 2D → 3D
• 6.2 Influence des différents réseaux cristallins — la connectivité comme critère distinctif
• 6.3 La transition supplémentaire solide(1) → solide(2)
• 6.4 Réflexions : l'anomalie acoustique de l'eau liquide
• 6.5 Alternative : champs de tourbillons fractals au lieu de la rotation des micelles
• 6.6 Rigidité phononique microscopique θ(réseau) : un modèle de ressorts d'Einstein
• 6.7 Estimation mathématique de la métastabilité (2D triangulé ↔ carré, 16b)
• 6.8 Affinement de type Debye de θ(réseau)
• 6.9 θ(réseau) à partir d'un réseau de Lennard-Jones complet
• 7 — Calibration et équation d'état
• 7.1 Méthode de calibration (un paramètre d'énergie ε) — ajustée
• 7.2 Tableau de calibration
• 7.3 Axe densité/volume : équation d'état van der Waals à cœur dur
• 7.4 Calibration quantitative l/g : a → ε, b → D²
• 7.5 Surface (τ, v) commune avec point triple
• 7.6 Chaleur latente face aux enthalpies de fusion
• 7.7 Unification des échelles d'énergie : ε ↔ ε_lg
• 7.8 Point triple sans paramètre à partir de z̄_c
• 7.9 Calcul de contrôle 3D : une fenêtre liquide nette sur une seule échelle ε
• 7.10 Véritable balayage de percolation de liens en 3D (CFC/CC) et sensibilité du bord de la fenêtre
• 7.11 Véritable balayage en τ de la cinétique de liaison en 3D : τ_b(3D) et la ligne de fusion sur une échelle ε
• 7.12 La ligne de transformation solide(1)→solide(2) dans la surface (τ, v)
• 7.13 Ligne de transformation solide(1)→solide(2) dans le diagramme (τ, p), calibrée sur le fer
• 7.14 Redlich–Kwong comme branche séparée de l'image idéale
• 7.15 solide(1)→solide(2) sur d'autres polymorphes : titane et contrôle par gaz rare
• 8 — Comparaison bibliographique : les conceptions de la thermodynamique statique
• 8.1 Tableau d'ensemble
• 8.2 Les modèles un à un
• 8.3 Ressemblances et différences avec le modèle 3
• 9 — Annexe
• 9.1 Bibliographie
• 9.2 Assistance de l'IA
• 9.3 Glossaire des abréviations et des symboles
• 9.4 Scripts de calcul et couverture de tests

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

1 — Les axiomes et leurs mots-clés (consolidés)

Résumé. Sept axiomes pourvus de mots-clés fixent la substance idéale : nombre de particules conservé, mouvement purement classique, liaison coopérative à deux états et nombre de coordination moyen z̄ comme paramètre d'ordre. Ils portent toutes les dérivations qui suivent ; dans le texte maître allemand, un filigrane d'homoglyphes marque les mots-clés (clé en fin de texte).

Les axiomes précisés :

1. VOL LIBRE (axiome supplémentaire). Entre deux molécules qui passent l'une près de l'autre sans être liées, il n'agit ni frottement ni forces électrostatiques. L'impulsion et l'énergie ne changent que par la répulsion dure du CŒUR DUR lors du choc et par le nouage et la rupture d'une LIAISON. Cet axiome n'importe que du point de vue de la mécanique classique et du bilan énergétique : il garantit qu'il n'existe aucun canal dissipatif au niveau moléculaire — l'énergie du mouvement n'est jamais « perdue », elle passe seulement dans les trois réservoirs de la TRIADE ou en énergie de liaison. (Le frottement et la viscosité macroscopiques naissent de façon émergente du taquin et de la FISSURE, non comme force fondamentale microscopique.)
2. MULTITUDE. Une substance idéale est faite de très nombreuses molécules identiques ; les propriétés macroscopiques sont décrites statistiquement dans la limite thermodynamique (N → ∞).
3. LIAISON. Deux molécules géométriquement voisines peuvent porter une liaison. Une liaison intacte abaisse l'énergie de la quantité finie ε ; sa rupture coûte exactement ε. La liaison est un système à deux états {intacte, rompue}.
4. TRIADE. L'énergie interne se répartit classiquement sur trois réservoirs : (a) la translation des molécules, (b) la vibration dans les liaisons intactes (un ressort harmonique par liaison), (c) la rotation de blocs moléculaires rigides (degré de liberté de rotation collectif des domaines connexes).
5. CŒUR DUR. Si la distance moléculaire passe sous une valeur minimale σ, une répulsion dure agit (volume propre). Avec la LIAISON il en résulte un potentiel de paire effectif (coque dure plus puits de liaison fini).
6. CONSERVATION. L'impulsion totale et l'énergie totale sont conservées. Une molécule isolée est énergétiquement neutre — elle ne met aucune énergie en réserve en elle-même ; l'énergie vit dans les liaisons (LIAISON) et dans les degrés de liberté cinétiques (TRIADE).
7. FISSURE. Les liaisons se défont et se forment selon le principe fissuration/cicatrisation : plus il manque de liaisons voisines → plus la probabilité de rupture est grande ; plus il y a de voisins liés → plus la probabilité de cicatrisation est grande. La cinétique est coopérative.
8. PROMPT. Rupture et cicatrisation se déroulent vite et avec peu d'activation ; le sous-système des liaisons se trouve à toute température en équilibre local et peut être séparé adiabatiquement du mouvement moléculaire (fondement de la fonction de partition factorisée).

Posture de modélisation. La coordination géométrique est idéalement conservée ; seul le nombre de coordination z change par les ruptures. Le mouvement naît comme dans un taquin : les cavités (issues du volume propre du CŒUR DUR et des ruptures de la FISSURE) migrent, les blocs se déplacent. L'image délibérément « fausse » de cristallites tournants et migrants est conservée comme fondement intuitif.

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

2 — Analyse des contradictions et des lacunes

Résumé. L'analyse examine systématiquement le modèle quant aux contradictions internes et aux lacunes et poursuit l'état du travail ; les faiblesses reconnues sont nommées puis résolues dans les chapitres suivants ou marquées comme limite honnête du modèle.

Liste des points critiques

1. Double comptage de la translation (axiome 3a contre 3c). Le texte brut nomme la translation deux fois (translation des molécules et translation des blocs). Deux réservoirs de translation indépendants n'existeraient pas. Résolution : la translation de bloc n'est qu'un mode collectif des translations moléculaires déjà comptées ; seul est véritablement nouveau le degré de liberté de rotation du bloc. Dans le bilan des degrés de liberté (chap. 4) la translation est comptée exactement une fois.
2. « Énergétiquement neutre » (axiome 5) contre « les liaisons absorbent de l'énergie » (axiome 2/3b). Contradiction apparente. Résolution : la neutralité concerne la molécule en tant que corps ; l'énergie siège entre les molécules (dans la LIAISON) et dans des modes cinétiques collectifs, jamais dans un tampon intramoléculaire. Molécule ≠ liaison — la contradiction est ainsi levée.
3. Attraction absente (l'axiome 4 ne nomme que la répulsion). Sans attraction il n'y aurait pas de phase condensée. Résolution : l'attraction est fournie par la LIAISON (puits de profondeur ε), la répulsion par le CŒUR DUR (coque dure en σ). Ensemble, un potentiel de paire effectif.
4. Rétroaction sans terme opposé (axiome 6). « Plus de voisins manquants → plus de ruptures » est à soi seul une rétroaction positive (propagation de la fissure jusqu'à la désagrégation totale). Résolution : le terme de cicatrisation (plus de voisins → plus de cicatrisation) et le facteur de Boltzmann de l'énergie de liaison finie ε (LIAISON) équilibrent la fissure et produisent une transition définie au lieu d'un emballement.
5. Séparation des échelles de temps non explicitée (axiome 7). « Rapide et peu activée » ne devient utilisable que si l'on en tire la séparation adiabatique entre dynamique des liaisons et dynamique de translation. Résolution : posée explicitement comme condition de la factorisation Z = Z_trans·Z_vib·Z_rot·Z_bond.
6. Géométrie absente et mécanisme du mouvement absent. Les axiomes ne nomment ni le réseau de coordination ni ce qui permet l'écoulement. Résolution (comme hypothèse additionnelle explicitement signalée) : une géométrie de coordination de référence fixe avec occupation variable des liaisons ; comme cas 2D canonique la coordination triangulaire (z_max = 6). La mobilité par le mécanisme des cavités du taquin. Ces deux points sont les seuls véritables ajouts au-delà des sept axiomes et sont signalés au chapitre 6 comme question stratégique.
7. Lacune latente de dissipation. La version 07 laissait ouvert le point de savoir si des molécules non liées perdent de l'énergie en se croisant. Résolue par le VOL LIBRE : il n'y a ni frottement ni électrostatique microscopiques ; tous les flux d'énergie et d'impulsion passent exclusivement par les chocs du CŒUR DUR et la cinétique de LIAISON. Le bilan énergétique (chap. 4.3) est ainsi clos et la statistique de Maxwell–Boltzmann exacte dans la limite gazeuse (gaz parfait libre entre chocs élastiques).
8. Incohérence énergétique du paramètre d'ordre. La version 08 utilisait le terme coopératif −(K/2)x² dans l'énergie libre g(x) (pour engendrer x(τ)), mais pas dans l'énergie interne U/N. C'est thermodynamiquement incohérent : le même terme de champ moléculaire doit apparaître dans les deux. Résolution (chap. 4.3) : −(K/2)x² est reporté dans U/N ; la chaleur latente s'en trouve mesurablement accrue (chap. 5.3), sans que les plateaux solide et gazeux changent.
9. Calibration l/g seulement qualitative. La version 09 n'avait l'axe des densités (équation d'état de van der Waals) que sous forme réduite. Résolution (chap. 7.4) : la ligne d'ébullition absolue est calculée via le point critique (T_c, p_c) et T_b(p_atm) est confronté aux mesures ; l'écart systématique (règle de Guldberg) est nommé.
10. Axes séparés sans point triple. Les axes de fusion et d'ébullition étaient des calculs séparés. Résolution (chap. 7.5) : les deux sont réunis dans un plan (τ, v), le point triple est construit et le point critique marqué.

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

3 — Descriptions imagées

Résumé. Le chapitre 3 explique de façon intuitive le solide, le liquide et le gaz ainsi que leurs transitions (avec des bandes dessinées statiques) : les trois phases isolées (3.1), les transitions de phase recristallisation/fusion/ébullition/sublimation avec le point triple et le point critique (3.2), l'anomalie acoustique (3.3), le mouvement brownien et la notion affinée de la micelle (3.4–3.5), ainsi que la fenêtre liquide tridimensionnelle (3.6). À partir de ces descriptions ont été réalisés des schémas (bandes dessinées) destinés à aider le lecteur à se représenter, devant l'œil intérieur, les rapports moléculaires dynamiques dans les différentes phases. Cette capacité à se représenter les choses est le véritable noyau du modèle. Il devient qualitativement et scientifiquement précieux dès lors que la mécanique statistique permet de montrer de façon plausible que le modèle rend prévisibles des tendances et observations empiriques qualitatives, même sans grand appareil mathématique.

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

3.1 Les trois phases isolées

Solide. Presque tous les voisins sont liés entre eux par des liaisons moléculaires (LIAISON) ; la moyenne des nombres de coordination z des molécules est proche du maximum géométrique. Le réseau moléculaire est un unique bloc rigide sans cavités (CŒUR DUR, empilement dense) — un taquin interne comme dans le liquide n'est pas possible dans la version compacte. Seul mouvement : le tremblement autour des positions de repos (vibration de liaison, TRIADE-b), déterminé par les liaisons moléculaires avec les voisins. Impulsion et énergie oscillent élastiquement dans le réseau (CONSERVATION, VOL LIBRE : sans perte). Dans le cas limite (très basse température, empilement idéal) on a l'état idéal d'un cristal parfait.

Liquide. La FISSURE coopérative a traversé le réseau : des liaisons manquent par plaques, la moyenne des nombres de coordination z des molécules a baissé, le bloc rigide s'est intérieurement brisé en de nombreux blocs plus petits. Du monolithe est né un tas d'éboulis. Ces blocs translatent et tournent (TRIADE-c). Entre eux se sont ouvertes des cavités par lesquelles ils se déplacent les uns près des autres — le taquin du mouvement brownien est actif. La géométrie de coordination locale reste idéalement conservée ; ce qui a changé, c'est seulement combien de liaisons existent réellement et se forment.

Gaz. La FISSURE l'a emporté : les liaisons sont rompues presque partout (z → 0). Dans le cas extrême, plus aucune molécule n'a de liaison moléculaire avec une autre. Les blocs se sont désagrégés en molécules isolées. Elles volent librement, couplées seulement par les chocs du CŒUR DUR (avec du VOL LIBRE entre-temps) ; les vitesses suivent Maxwell–Boltzmann (CONSERVATION). Vibration de liaison et rotation de bloc se sont éteintes — translation pure (TRIADE-a). Dans le cas limite (haute température, faible pression) on a l'état idéal d'un gaz parfait.

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

3.2 Transitions de phase

Recristallisation (solide(1)/solide(2)). Encore en dessous de la fusion, le solide peut changer de réseau sans perdre sa rigidité : deux empilements idéaux de nombres de coordination différents (CFC z_max = 12 ↔ CC z_max = 8) ont des enthalpies libres μ(τ, p) différentes ; à μ_CFC = μ_CC la géométrie de coordination bascule brusquement (et pas seulement l'occupation z̄ comme lors de la fusion). Apparaît alors un petit palier de c_v (Δc_v = variation de z_max/4 = 1,0 k_B) plus une faible chaleur latente de transformation. Affirmation de l'IA : cela n'est concevable qu'en 3D (connectivité discrète, quantifiée ; chap. 6.3). Dans le cas bidimensionnel, la véritable recristallisation, c'est-à-dire la transition entre deux phases solides, se laisse imaginer plus facilement. Dans le plan il n'existe que trois pavages réguliers sans trous. On peut recouvrir entièrement une surface (la paver) avec des triangles équilatéraux, avec des carrés ou avec des hexagones réguliers. Si, dans ces réseaux, les arêtes représentent les liaisons moléculaires et les sommets les molécules, alors une transition de phase solide(triangulaire) → solide(carré) est concevable. La molécule du pavage triangulaire aurait alors encore six liaisons vers les molécules voisines, tandis que celle du réseau carré n'en aurait plus que quatre.

Fusion (s/l). Quand la température monte, le réseau se déchire en assez d'endroits pour que le bloc rigide d'ensemble perde sa rigidité (la FISSURE coopérative passe sous le seuil de rigidité de Maxwell z̄ = 2d — le réseau devient mou). Il se forme des micelles thermodynamiques, c'est-à-dire des domaines qui ne sont plus reliés par des liaisons aux domaines voisins plus grands. Du point de vue moléculaire, le monolithe devient un tas d'éboulis et acquiert une capacité d'écoulement limitée. Les premières cavités permettent le déplacement de blocs moléculaires, comme dans le taquin. Les liaisons manquantes vers le voisinage permettent la rotation interne dans les micelles. La rotation des micelles (TRIADE-c) s'amorce. La rupture de nombreuses liaisons (LIAISON) engloutit de la chaleur latente.

Ébullition (l/g). Un chauffage supplémentaire fait percoler les cavités elles-mêmes : elles se relient en un espace vide continu, dans lequel les blocs se dissolvent et des molécules isolées s'échappent. L'ébullition est le passage de « cavités dans la matière » à « matière dans une cavité ».

Sublimation (s/g). À basse pression, le domaine liquide intermédiaire manque : le réseau cède directement des molécules isolées à l'espace vide, sans qu'il y ait formation de blocs comme dans les micelles. On peut se demander ici quelle influence peuvent avoir sur la formation ou la non-formation de micelles thermodynamiques des aspects géométriques (molécules longues, molécules qui s'enchevêtrent entre elles), des aspects électrostatiques (alignement par moments magnétiques ou moments dipolaires) et différents types de liaison moléculaire (van der Waals, interactions entre systèmes mésomères, liaisons hydrogène).

Point triple. C'est précisément là que deux seuils coïncident : la rigidité de Maxwell (solide ↔ liquide, perte de la rigidité) et la percolation des cavités (liquide ↔ gaz) surviennent au même (τ, p). Les trois images (bloc rigide [monolithe], blocs qui glissent [tas d'éboulis], molécules libres [nuage de poussière]) coexistent.

Point critique. Quand la pression monte, le saut liquide→gaz devient de plus en plus petit ; la fraction de vide varie de façon de plus en plus continue. Au point critique la différence disparaît tout à fait — « des blocs denses avec peu de vide » et « beaucoup de vide avec peu de molécules » passent continûment l'un dans l'autre. (« Du point de vue du liquide le verre est à moitié vide, du point de vue du gaz à moitié plein. »)

Recristallisation dans l'image 2D. La bande dessinée montre désormais la recristallisation comme cas 2D, l'un sous l'autre : en bas le réseau dense triangulé (coordination z = 6, l'idéal 2D de plus basse énergie), en haut le réseau carré lâche (z = 4). La transition triangulé → carré augmente le volume (travail contre la pression) et abaisse le nombre de liaisons. On voit que le nombre de coordination change alors que la rigidité d'un réseau cristallin subsiste.

Prédiction qualitative des transitions d'énergie et des capacités thermiques. À partir de la seule connectivité z, les deux grandeurs caractéristiques s'estiment : (i) l'énergie de liaison par particule vaut U/N ≈ −(ε/2)·z (chaque liaison compte pour moitié), donc d'autant plus basse que z est grand ; (ii) le saut de c_v au changement de réseau suit la règle du modèle Δc_v = variation de z_max/4. Pour le cas 2D triangulé (z = 6) ↔ carré (z = 4) cela donne une différence d'énergie de ΔU ≈ ε par particule et un palier de c_v de Δc_v = (6 − 4)/4 = 0,5 k_B — moitié moins que la valeur 3D (1,0 k_B). Règle mnémotechnique : plus dense = plus bas en énergie, mais plus coûteux dans la transition (chaleur latente plus grande, palier de c_v plus grand), parce que la dilatation libère plus de liaisons et plus de degrés de liberté gelés.

La structure du réseau détermine-t-elle la métastabilité ? En partie. La structure du réseau fixe l'ordre des énergies — quelle forme est l'état fondamental (triangulé, μ le plus bas) et laquelle est excitée (carré). Que la forme lâche survive comme métastable à basse température, elle ne le décide pas à elle seule : dans le modèle strictement idéal (sans barrière), la substance adopte à chaque (τ, p) immédiatement la phase de μ le plus bas — la métastabilité n'y existe pas. Elle naît seulement d'une énergie d'activation finie (barrière de germination) qui freine cinétiquement la transformation. La géométrie du réseau fournit donc la hiérarchie thermodynamique, la barrière la persistance cinétique — le même ingrédient qui explique aussi la phase sursaturée (voir plus bas ; quantitativement chap. 6.7).

Phase sursaturée — énergie d'activation. La figure additionnelle (3.2, rangée 2) montre une phase gazeuse sursaturée au-dessus d'un liquide : à côté de molécules isolées, de petits groupes de 2 à 3 particules agrégées (barres grises) qui volent d'un bloc. Le modèle explique la sursaturation par une énergie d'activation notable pour les liaisons moléculaires : les petits amas (2–3) ont un rapport bord/surface défavorable, leurs liaisons de bord sont incomplètes — seul un germe critique croît spontanément. Tant qu'aucun germe surcritique n'apparaît, le gaz reste sursaturé, bien que μ_gaz > μ_liquide permette déjà la condensation. C'est précisément cette barrière qui manque dans l'image d'équilibre idéale et qui est nommée ici comme une contribution autonome située « en dehors de l'image idéale » (estimation quantitative au chap. 6.7).

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

3.3 Pourquoi la vitesse du son dans les liquides décroît avec T

Le son est une onde élastique : v = √(K_eff/ρ). K_eff naît exclusivement des liaisons intactes (LIAISON) — ce sont elles, les ressorts. Quand T monte, de plus en plus de liaisons se rompent coopérativement (FISSURE) ; le réseau de ressorts s'amincit et K_eff chute plus vite que la densité ρ. La rupture domine le numérateur, pas le dénominateur — c'est pourquoi v décroît. Interprétation moléculaire : la transmission du son par les liaisons moléculaires est plus rapide que la transmission mécanique par des chocs. Plus il apparaît de fissures dans un liquide à mesure que la température augmente, plus souvent le son doit être transporté par des chocs mécaniques plutôt que par des liaisons moléculaires. En conséquence, la vitesse du son devient plus lente quand la température monte. (L'anomalie de l'eau — un maximum acoustique vers 74 °C — est traitée au chap. 6.4.)

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

3.4 Comment naît le mouvement brownien — avec une notion affinée de la micelle

Au sens de la thermodynamique statistique, une micelle n'est pas ici un agrégat de tensioactifs mais un bloc cerné par la FISSURE : les liaisons manquantes (les non-liaisons à la frontière du bloc) forment la « coque » autour d'un cristallite enfermé (le cœur intactement lié). Ce cœur peut exécuter un mouvement de rotation coopératif interne (TRIADE-c) — le bloc tourne d'un seul tenant tandis que sa coque frontière le découple des voisins.

Le mouvement brownien naît alors ainsi : par l'échange de liaisons moléculaires (FISSURE, PROMPT), des molécules isolées peuvent quitter la coque d'une micelle et passer dans une autre micelle. Chaque changement de ce type transmet de l'impulsion (CONSERVATION ; sans perte grâce au VOL LIBRE) et modifie brusquement la forme et le moment cinétique des deux micelles concernées. Une grosse particule d'essai (grain de poussière) est baignée par de nombreuses réorganisations de micelles non corrélées ; la somme des impulsions est nulle en moyenne, mais de variance finie → marche aléatoire. Plus T est élevé, plus les événements de FISSURE et les changements de micelle sont fréquents, plus le tremblement est fort. Le mouvement brownien est ainsi l'ombre macroscopique, ou la conséquence, du processus microscopique d'échange entre micelles.

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

3.5 Pourquoi les réactions aux couches limites des micelles s'accélèrent

Remarque préliminaire embarrassante. L'argumentation qui suit est une déduction tirée de l'image et, à ma connaissance, n'a jusqu'ici été ni démontrée ni étudiée empiriquement — je dois avouer que je ne me suis pas donné la peine de faire une recherche approfondie à ce sujet. Le seul indice de nouveauté est que les systèmes d'IA ne mentionnent pas cette conclusion. Cette thèse illustre comment des représentations imagées peuvent conduire à de nouvelles questions.

Dans des solvants moyennement dissolvants, les réactifs se fixent de préférence aux coques frontières des micelles — dans le modèle, précisément les frontières fissure/bloc (FISSURE). Deux effets accélèrent la réaction : (i) la réduction dimensionnelle — les réactifs ne se cherchent plus dans le volume entier mais le long d'une interface quasi bidimensionnelle ; une recherche 2D aboutit statistiquement plus vite qu'une recherche 3D. (ii) La concentration — la coque frontière rassemble les réactifs, leur densité locale augmente. Imagées, les frontières de fissure agissent comme des rails de guidage qui conduisent les partenaires l'un vers l'autre.

L'effet de dimension a un noyau rigoureux : la marche aléatoire (Pólya). Qu'une recherche bidimensionnelle aboutisse statistiquement mieux qu'une recherche tridimensionnelle ne résulte pas seulement de l'intuition, mais d'un théorème du calcul des probabilités. George Pólya (1887–1985) a montré en 1921 pour la marche aléatoire sur un réseau : en une et en deux dimensions le marcheur revient avec probabilité 1 (donc à coup sûr) une infinité de fois en chaque point, sur le réseau plan ℤ² comme sur la droite ℤ. À partir de trois dimensions cela ne vaut plus : la probabilité de retour est strictement inférieure à 1 (sur le réseau cubique simple ℤ³, environ 0,34 seulement). De façon imagée, selon Pólya : « Un homme ivre finit sûrement par rentrer chez lui, un poisson ivre en revanche non. » Transposé à deux réactifs qui doivent se rencontrer plusieurs fois pour un choc réussi, cela signifie : dans un très bon solvant, avec une marche aléatoire tridimensionnelle libre, des partenaires une fois séparés dérivent à jamais l'un loin de l'autre avec une probabilité finie — les tentatives répétées de réaction échouent, et le rendement reste faible pour cette raison purement géométrique et statistique. À une couche limite de micelle (quasi) bidimensionnelle, en revanche, la rencontre à nouveau est certaine (probabilité 1) : les partenaires se retrouvent encore et encore jusqu'à ce que la réaction aboutisse. La couche limite n'a donc pas seulement un effet concentrant, elle met en même temps à disposition la dimension la plus favorable aux retrouvailles. L'affirmation (délibérément appuyée) est la suivante : les rendements dans les très bons solvants tendent à être moins bons pour des raisons de marche aléatoire.

Référence : G. Pólya (1921), Über eine Aufgabe der Wahrscheinlichkeitsrechnung betreffend die Irrfahrt im Straßennetz, Math. Ann. 84, 149–160 — cité ici d'après M. Eisermann, Höhere Mathematik 3 (Univ. Stuttgart), exercice W-2×2. (Original consulté le 5.9.2026).

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

3.6 La fenêtre liquide tridimensionnelle

À la différence des figures de 3.1–3.5, celle-ci ne montre aucune scène moléculaire mais le résultat du calcul des chap. 7.9/7.10 comme schéma sur un axe τ (en unités de l'énergie de liaison ε). Le but était de rendre visible pourquoi le réseau 2D isolé n'admet guère de domaine pour une phase liquide, alors que c'est sans difficulté possible avec le réseau 3D, la même règle de rigidité (z̄_c = 2d) valant dans les deux cas.

Image ① — 2D (pas de fenêtre). Sur l'axe τ, τ_m = 1,20 et τ_triple = 1,16 se superposent pratiquement (τ_triple sous τ_m), τ_c = 1,50 juste au-dessus : le domaine de la phase liquide est comprimé en une fente — le solide passe quasi directement au gaz (proche de la sublimation).

Image ② — 3D (fenêtre nette). Le même axe, mais dilaté : τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00. Entre τ_triple et τ_c se trouve le domaine liquide autonome (large vers le bas, τ_triple − τ_m = 1,63 ; étroit vers le haut, τ_c − τ_triple = 0,18). Le mince bord supérieur est marqué comme tel.

Image ③ — la cause. z̄_c = 2d croît avec la dimension (4 → 6), d'où la baisse de v_solid (0,50 → 0,667 en volumes réduits, c'est-à-dire que le solide devient « plus lâche » relativement au liquide), et τ_c = z_max/4 double (1,50 → 3,00). Le point critique s'éloigne bien au-dessus du point de fusion. C'est seulement ainsi que naît la place pour la phase liquide. C'est un pur effet de dimension.

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

4 — Modèle de mécanique statistique (inchangé depuis la version 09)

Résumé

Le chapitre 4 met en place la charpente de mécanique statistique du modèle en 2D (réseau triangulaire, z_max = 6) :

• 4.1 Degrés de liberté et variables : la grandeur d'ordre est la fraction de liaisons x ∈ [0,1] avec le nombre de coordination moyen z̄ = z_max·x ; par bloc ≥ 2 s'ajoute un degré de liberté de rotation, et la température τ se compte en unités de l'énergie de liaison ε.
• 4.2 Liaison coopérative à deux états : l'équation de liaison auto-cohérente x = 1/(1 + G·exp[−(ε + K·x)/τ]) (K = J·z_max, G = z_max) découle de la condition de stationnarité de l'énergie libre de Bragg–Williams ; son minimum global fournit la branche physique, et le chap. 5.5 décide si celle-ci est continue ou présente un saut (premier ordre).
• 4.3 Bilan énergétique avec terme coopératif : U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x² (2D) réunit translation, vibration de liaison et rotation de bloc, l'énergie de configuration des liaisons et la nouvelle énergie coopérative de champ moléculaire −(K/2)x². La capacité thermique c_v associée reçoit le terme additionnel positif −K·x·ẋ, qui accentue le pic de fusion mais laisse intacts le plateau solide et la branche gazeuse.

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

4.1 Degrés de liberté et variables

• N molécules sur une géométrie de coordination de référence fixe. Cas 2D canonique : réseau triangulaire, z_max = 6, un DDL de rotation par bloc (MULTITUDE, LIAISON). Ce chapitre vaut explicitement pour l'interprétation 2D ; la généralisation 3D figure au chap. 6.
• Variable de liaison s_ij ∈ {0, 1} par paire géométrique de voisins (1 = intacte).
• Nombre de coordination z_i = Σ_j s_ij ; moyenne z̄ = 2 N_bond / N — le paramètre d'ordre. On a z̄ = z_max · x avec la fraction de liaisons x = p_b ∈ [0, 1].
• Unité d'énergie τ ≡ k_B T (classiquement, ½τ par degré de liberté quadratique). Toutes les énergies en unités de l'énergie de liaison ε.

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

4.2 Liaison coopérative à deux états (de type Ising, paramètre J)

Une liaison isolée (énergie −ε intacte, 0 rompue) a, à l'équilibre de Maxwell–Boltzmann, p_b⁰(τ) = 1/(1 + e^(−ε/τ)). La FISSURE coopérative (axiome 6) décale la stabilité effective de la liaison d'un champ moléculaire proportionnel à la fraction de liaisons x. En approximation de champ moléculaire (Ising), l'équation de liaison auto-cohérente devient :

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

Ici :

• K = J·z_max est la force de couplage coopérative (rétroaction de la FISSURE : plus de liaisons voisines → stabilité de liaison effective plus élevée). J > 0 accentue la transition ; au-delà d'un J critique, z̄(τ) devient discontinu — c'est la chaleur latente de la fusion. Le cas pur à deux états est la limite J = 0 (pas de champ moléculaire, crossover continu).
• G > 1 est la dégénérescence de l'état rompu : une molécule libérée gagne de l'espace des phases (liberté de translation). G pousse la rupture des liaisons à haut τ et fait que x → 1/(1+G) au lieu de rester bloqué à ½ (cas limite du gaz). Canoniquement G = z_max.

Cette équation découle de la condition de stationnarité d'une énergie libre de Bragg–Williams par molécule (n_slots = z_max/2 emplacements de liaison) :

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

Le terme −(K/2)x² est le couplage coopératif (ferromagnétique), le dernier terme la dégénérescence de rupture. Le minimum global de g sur x fournit la branche physique ; là où deux minima deviennent également profonds se situe la transition du premier ordre. PROMPT (axiome 7) justifie que x(τ) prenne à chaque température la valeur d'équilibre.

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

4.3 Bilan énergétique — avec terme coopératif

Le terme de champ moléculaire −(K/2)x² est une véritable énergie (indépendante de τ) et doit apparaître dans U/N. Le bilan cohérent s'écrit donc (2D) :

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

• τ·[…] : translation (1) + vibration de liaison (z̄/4) + rotation de bloc (½ n_blk) — TRIADE.
• −(ε/2)·z̄ : énergie de configuration des liaisons — LIAISON.
• −(K/2)·x² : énergie coopérative de champ moléculaire — FISSURE. Nouveau dans U/N.

La capacité thermique associée c_v = d(U/N)/dτ (en unités k_B) reçoit une contribution supplémentaire :

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

où le nouveau terme −K·x·ẋ est positif (ẋ < 0) et agit exactement là où de nombreuses liaisons se rompent en même temps. Il accentue le pic de fusion sans modifier le plateau solide (ẋ ≈ 0) ni la branche gazeuse.

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

5 — Fonctions de partition et capacités thermiques

Résumé

Le chapitre 5 calcule les fonctions de partition et les capacités thermiques et fournit les données du modèle pour la calibration (chap. 7) :

• 5.1 Fonction de partition factorisée : en raison de la séparation des échelles de temps, Z = Z_trans·Z_vib·Z_rot·Z_bond se factorise.
• 5.2 Capacité thermique par phase : c_v = d(U/N)/dτ ; le terme positif −(ε/2)(dz̄/dτ) engendre le pic de type Schottky, analogue de la chaleur latente.
• 5.3 Courbe numérique c_v(τ) : pic de c_v 5,51 à τ_m = 1,20, plateau solide 2,50, branche gazeuse → 1 ; la chaleur latente L ≈ 1,41 ε par molécule (avec le terme coopératif) est triplée par rapport à 0,48 ε (sans).
• 5.4 / 5.4a Deux seuils séparés : la percolation de connectivité z̄_c = z_max·p_c ≈ 2,08 (à τ_b ≈ 1,85) marque la désagrégation des amas (proche de l'ébullition) et est un résultat ; elle est à distinguer nettement de la rigidité générique de Maxwell z̄_c = 2d (densité de fusion, perte de rigidité).
• 5.5 Ordre de la transition de fusion : critère ln G > 2 + 4ε/K ; le 2D (G = 6) est continu pour tout J, le 3D (G = 12) est du premier ordre dès K > ~8,2 (démo : Δz̄ ≈ 5,3, ΔU ≈ 1,75 ε).
• 5.6 Étude de taille finie : p*(L) converge monotonement par en dessous vers p_c = 0,3473 ; z̄_c ≈ 2,08 est robuste en taille finie.
• 5.7 Balayage 3D complet de c_v : dès K > K_crit ≈ 8,25 (J > 0,68) s'ouvre le véritable saut Δz̄ avec chaleur latente ΔU ; le pic de fusion se déplace avec J vers les τ plus élevés (décalage positif de +0,5 … +2,2 par rapport au pic 2D continu).

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

5.1 Fonction de partition factorisée

En raison de la séparation des échelles de temps (PROMPT), la fonction de partition canonique se factorise :

Z = Z_trans · Z_vib · Z_rot · Z_bond.

• Z_trans = (1/N!)·(A/λ_T²)^N, λ_T = h/√(2π m τ) (translation 2D classique idéale ; MULTITUDE, CONSERVATION, VOL LIBRE).
• Z_vib = (τ/ħω₀)^(N_bond) (oscillateur harmonique classique par liaison intacte ; TRIADE-b).
• Z_rot = ∏_(blocs≥2) (2π I_b τ)^(1/2)/ħ (rotateur rigide 2D classique par bloc ; TRIADE-c).
• Z_bond = Σ_{s} exp[(ε/τ)·Σ s_ij + (K/τ)·terme coopératif] (somme de liaisons de type Ising ; LIAISON, FISSURE). L'évaluation en champ moléculaire fournit x(τ), z̄(τ) d'après 4.2.

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

5.2 Capacité thermique par phase (analytique)

À partir de U/N (chap. 4.3), c_v = d(U/N)/dτ (en unités k_B) :

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

Le terme −(ε/2)·(dz̄/dτ) est positif (dz̄/dτ < 0 : chauffer rompt des liaisons) et domine là où de nombreuses liaisons se rompent en même temps — le pic de type Schottky, analogue de la chaleur latente.

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

5.3 Courbe numérique c_v(τ) avec terme coopératif — chaleur latente quantifiée

Liaison auto-cohérente (minimum global), n_blk(τ) issu d'un Monte-Carlo de percolation de liens (union-find, L = 100, 5 réalisations), c_v = dU/dτ (différences centrées), U/N y compris −(K/2)x² (cv_scan_09.py, ε = 1, J = 0,55, G = z_max = 6). Points d'appui choisis :

Résultat (numérique). L'ordre des trois phases reste confirmé et est désormais plus net :

c_v(solide) ≈ 2,50 < c_v(liquide, maximum) ≈ 5,51 > c_v(branche gazeuse) ↓ → 1.

Chaleur latente (aire du pic au-dessus de la ligne de base). En intégrant (c_v − ligne de base) sur la fenêtre de fusion τ ∈ [0,70 ; 1,70] (ligne de base = droite passant par c_v(0,70) et c_v(1,70)), on obtient

L ≈ 1,41 ε par molécule (avec terme coopératif), L ≈ 0,48 ε par molécule (sans terme coopératif, ancien bilan), ΔL ≈ +0,92 ε par molécule.

Le terme coopératif relève le pic d'un facteur ≈ 1,7 (3,25 → 5,51) et triple la chaleur latente. Physiquement : la rupture simultanée de nombreuses liaisons voisines (FISSURE) stocke non seulement de l'énergie de configuration (−ε/2·z̄), mais aussi de l'énergie de champ moléculaire (−K/2·x²) — les deux sont libérées lors de la fusion.

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

5.4 Seuil de percolation de connectivité comme résultat

La fraction de l'amas géant P_∞(τ) s'effondre à τ ≈ 1,85 ; là x ≈ 0,347 = p_c(triangle) = 2·sin(π/18) = 0,3473, donc z̄_c = z_max·p_c ≈ 2,08 — un résultat, non une donnée imposée. C'est le seuil de connectivité (amas traversant) et il marque la désagrégation des amas près de l'ébullition (τ_b) ; il ne doit pas être confondu avec la rigidité de Maxwell (z̄_c = 2d, fusion) — cf. la séparation nette au chap. 5.4a.

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

5.4.a Deux seuils distincts nettement séparés

La double dénomination jusqu'ici floue « percolation/rigidité » est maintenant levée. Il s'agit de deux seuils physiquement différents qu'il ne faut pas confondre :

Justification. La fusion est la perte de la rigidité mécanique d'un bloc dense, non la désagrégation de la connexité. La grandeur pertinente est donc le comptage des contraintes de Maxwell (champ moyen / rigidité générique) : un réseau à D dimensions de forces centrales devient isostatique exactement lorsque la coordination moyenne atteint z̄ = 2D (2D → 4, 3D → 6). La connectivité de liens (z̄_c = 2,08) est nettement en dessous ; un réseau tout juste connexe à z̄ = 2,08 est depuis longtemps mou (modes souples), donc déjà liquide. C'est pourquoi la connectivité marque la désagrégation des amas plus tardive, proche de l'ébullition, tandis que la rigidité fixe la solidification/liquéfaction plus précoce (densité de fusion). Cette séparation fixe le point triple du chap. 7.8 sans paramètre et de façon réaliste ; elle ne change pas le seuil de connectivité 2,08 confirmé en 5.4/5.6, mais seulement sa lecture (τ_b = proche de l'ébullition, non la fusion).

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

5.5 Ordre de la transition de fusion

Critère analytique. Le domaine spinodal (deux minima de g) existe pour τ < K/4 ; la transition se situe approximativement à τ_t ≈ (ε + K/2)/ln G. Un véritable saut (transition du premier ordre) exige τ_t < K/4, donc dans la limite K → ∞ :

ln G > 2 + 4ε/K.

Il en découle un résultat remarquable sur la dégénérescence de rupture G = z_max :

Confirmation numérique (max. |Δx| par pas en τ, detect_first_order) :

Lecture. Dans le modèle 2D canonique (G = 6), la fusion est par principe continue (transition « douce » / crossover) — pour aucun J n'apparaît un véritable saut, car ln 6 < 2. Seule la coordination/dégénérescence plus élevée de l'empilement 3D (G = 12) permet, à partir d'une coopérativité suffisante (K > 8,2), une transition du premier ordre avec chaleur latente finie. Cela concorde avec la réalité : la fusion 2D est souvent continue (de type KTHNY), la fusion 3D est du premier ordre. J = 0 (sans champ moléculaire) est toujours le cas de référence continu. — Physiquement, J désigne le rapport de l'énergie coopérative à l'énergie de liaison (K/ε = J·z_max/ε) : à quel point une liaison voisine intacte stabilise la sienne propre.

Démonstration du premier ordre (motivée par la 3D, cv_scan_09.py, J = 1,0, G = 12, z_max = 12) : z̄ saute de Δz̄ ≈ 5,30 à τ ≈ 2,82 ; la chaleur latente comme saut d'énergie ΔU ≈ 1,75 ε par molécule. Le même K sur un réseau G = 6 (2D) reste continu — le caractère de l'ordre ne dépend que de G = z_max.

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

5.6 Étude de taille finie du seuil de percolation

Percolation de liens sur des réseaux triangulaires L × L (périodiques), fraction d'amas géant P_∞(p) mise à l'échelle sur L (cv_scan_09.py ; L = 50/100/200/400 avec 8/5/3/2 réalisations). Seuil exact p_c = 2·sin(π/18) = 0,34730 → z̄_c = z_max·p_c = 2,0838.

P_∞(p) (moyenne ± dispersion) :

Estimation du seuil p*(L) (P_∞ croise 0,5) :

Résultat. p*(L) converge monotonement par en dessous vers p_c = 0,3473 ; la dispersion à p_c tombe de ±0,103 (L = 50) à ±0,017 (L = 400) — la transition s'affine avec L mais ne déplace pas le seuil (valeur de percolation universelle). Ainsi z̄_c ≈ 2,08 = z_max·p_c est confirmé comme valeur du modèle robuste en taille finie.

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

5.7 Balayage 3D complet de c_v(τ) pour la fusion du premier ordre

L'unique série de démonstration du premier ordre du chap. 5.5 (J = 1,0, G = 12) est maintenant développée en un balayage complet de c_v(τ) sur la coopérativité J (firstorder_3d_16.py, connectivité CFC pour n_blk, bilan énergétique 3D U/N = τ·[3/2 + z̄/4 + (3/2)n_blk] − (ε/2)z̄ − (K/2)x²). Le critère analytique d'ordre du chap. 5.5 (premier ordre ⇔ K > K_crit = 4ε/(ln 12 − 2) = 8,25) sépare les cas ; le pic de référence 2D continu se situe à τ_2D = 1,20 (c_v ≈ 5,5).

*cas continu : Δz̄/ΔU sont le crossover abrupt (pas de véritable saut) ; τ_transition y est le pic de c_v.

Résultat. (i) Le cas 3D continu J = 0,55 donne τ_transition = τ_m ≈ 1,72 — la vérification croisée avec la cinétique de liaison indépendante (bondkinetics_3d_15.py : τ_m(CFC) = 1,70) concorde. (ii) À partir de K > K_crit ≈ 8,25 (J > 0,68) s'ouvre le véritable saut Δz̄ avec une chaleur latente finie ΔU (saut d'énergie à travers la discontinuité) — la fusion devient du premier ordre, exactement comme l'exige le critère ln 12 > 2. (iii) La vérification croisée J = 1,0 reproduit la valeur de démonstration Δz̄ ≈ 5,4 du chap. 5.5. (iv) La hauteur du saut Δz̄ et la position τ_transition croissent avec J : la transition se déplace vers le haut — la coordination plus élevée de l'empilement 3D et la coopérativité plus forte relèvent la température de fusion. Par rapport au pic 2D continu (τ_2D = 1,20), on obtient donc partout un décalage positif du pic (+0,5 … +2,2) qui croît de façon monotone avec K. Le caractère du premier ordre de la fusion 3D est ainsi attesté non plus en un seul point, mais sur tout le domaine de coopérativité.

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

6 — Plausibilité mathématique de l'interprétation 3D

Résumé

Le chapitre 6 rend la généralisation de l'image des liaisons de 2D à 3D plausible de façon purement mathématique (sans visualisation HTML propre) :

• 6.1 Degré de liberté de rotation 2D → 3D : la rotation des blocs passe d'un axe à trois, la translation de deux directions à trois, z_max de 6 à 12 ; les valeurs de base de c_v augmentent (c_v,solide 2,5 → 4,5, c_v,gaz 1 → 3/2), la forme de la courbe c_v(τ) demeure.
• 6.2 La connectivité comme critère distinctif : z_max est une grandeur quantifiée et sépare brusquement l'empilement le plus dense (CFC/HC, z_max = 12) du CC (z_max = 8), sans état intermédiaire.
• 6.3 Transition supplémentaire solide(1)→solide(2) : de la discrétude découle la polymorphie CFC↔CC avec un palier de c_v sans paramètre de 1,0 k_B ; dans l'image idéale stricte (θ_CC = θ_CFC) elle n'apparaît pas, et seul un CC plus mou (θ_CC < θ_CFC, sous le seuil d'ouverture 0,61) ouvre τ_ss < τ_m.
• 6.4 / 6.5 Réflexions : l'anomalie acoustique de l'eau liquide comme changement d'ordre local estompé entre deux types de réseau, et une alternative de champs de tourbillons fractals à la rotation des micelles (purement qualitatif).
• 6.6 / 6.8 / 6.9 θ(réseau) microscopique à trois niveaux : modèle de ressorts d'Einstein (θ_CC/θ_CFC ≈ 0,82…0,94), densité d'états de phonons de type Debye (≈ 0,87…1,00) et réseau de Lennard-Jones complet (≈ 1,0 ; CC idéal instable au cisaillement, CFC état fondamental stable) — tous au-dessus du seuil 0,61 ; la polymorphie réelle se situe donc hors de l'image idéale harmonique.
• 6.7 Métastabilité (2D) : calcul d'ordre de grandeur pour la transformation triangulé ↔ carré (ΔU = +1 ε, ΔA ≈ +15,5 %) avec la barrière de germination ΔG\* = πλ²/Δμ_A comme apport cinétique.

Ce chapitre rend la généralisation 3D plausible de façon purement mathématique. La 3D n'est jamais visualisée dans une version HTML. Le but est seulement de montrer comment le degré de liberté de rotation et la connectivité changent de 2D à 3D et quelles nouvelles transitions deviennent ainsi concevables.

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

6.1 Changement du degré de liberté de rotation 2D → 3D

Le seul point qualitativement nouveau par rapport à la 2D est la rotation des blocs :

• 2D : un bloc rigide a 1 DDL de rotation (rotation dans le plan) → U_rot/N = ½ n_blk τ.
• 3D : un bloc rigide a 3 DDL de rotation (trois axes d'Euler) → U_rot/N = (3/2) n_blk τ.

En même temps la translation passe de 2 à 3 DDL (U_trans/N = (3/2)τ) et le nombre de coordination maximal de z_max = 6 (triangle) à z_max = 12 (empilement compact de sphères). Le bilan énergétique en 3D s'écrit donc

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

et les valeurs de base de la capacité thermique se décalent de façon cohérente :

La structure de la courbe (plateau solide, maximum liquide par −(ε/2)dz̄/dτ, chute vers la valeur du gaz) reste identique ; seuls les niveaux varient avec la dimension et z_max.

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

6.2 Influence des différents réseaux cristallins — la connectivité comme critère distinctif

En 3D il faut distinguer au moins deux réseaux idéaux, qui se séparent par leur connectivité (nombre de coordination z_max) :

La connectivité z_max est le critère distinctif décisif : elle fixe à la fois le niveau de capacité thermique du solide (z_max/4) et le seuil de percolation z̄_c. Comme chaque réseau n'admet qu'un nombre de coordination fixe, la connectivité est une grandeur quantifiée et finie — il n'existe pas de transition continue entre z_max = 12 et z_max = 8, mais seulement le saut discret d'un type de réseau à l'autre.

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

6.3 La transition supplémentaire solide(1) → solide(2)

De la finitude/discrétude de la connectivité découle une transition absente du modèle 2D : solide(1) → solide(2) (polymorphie/allotropie). Deux phases solides de z_max différent (par exemple CFC ↔ CC) ont des enthalpies libres μ(τ, p) différentes ; à la combinaison température/pression μ_{solide(1)} = μ_{solide(2)} la substance change de réseau sans fondre. Dans l'image du modèle, la géométrie de coordination change brusquement (et pas seulement l'occupation z̄ comme lors de la fusion) — la seule situation où la géométrie de référence, sinon fixée, bascule elle-même. Serait observable un petit palier de c_v (variation de z_max/4) plus une chaleur latente de transformation, nettement plus faible que le pic de fusion.

Le critère du premier ordre ln G > 2 + 4ε/K (chap. 5.5) fournit maintenant une raison concrète pour laquelle c'est justement l'interprétation 3D qui est qualitativement plus riche : seule là (G = z_max = 12, ln 12 = 2,485 > 2) la fusion peut devenir du premier ordre — et seule là existe le changement discret de réseau solide(1) → solide(2) (CFC ↔ CC), parce que la connectivité est quantifiée/finie. Les deux phénomènes sont la même cause sous deux habits : la coordination plus élevée et discrète de l'empilement 3D engendre à la fois le saut de fusion net et la polymorphie. Le réseau 2D isolé (G = 6) n'a ni l'un ni l'autre — en cohérence avec le chap. 5.5.

Quantification (bondkinetics_3d_15.py). La transition est maintenant calculée de façon épurée. Chacun des deux réseaux reçoit une enthalpie libre dans l'image harmonique classique (les parts communes p·v et de translation s'annulent dans la comparaison) :

μ(τ) = h₀ − c_v·τ·ln(τ/θ), h₀ = −(ε/2)·z_max, θ = ħω/k_B (rigidité phononique en unités ε),

avec les valeurs caractéristiques sans paramètre des chap. 6.1/6.2 :

D'où le palier de c_v sans paramètre : |Δc_v| = variation de z_max/4 = (12−8)/4 = 1,0 k_B (lors de la transition CFC→CC le solide devient plus « léger » de 1 k_B en capacité thermique) — bien plus petit que le pic de fusion (plusieurs k_B, chap. 5.3/7.11).

La température de transformation τ_ss (μ_CFC = μ_CC) et la chaleur latente de transformation L = τ_ss·(S_CC − S_CFC) :

Résultat (constat honnête). Dans l'image idéale stricte (même rigidité phononique pour toutes les liaisons, θ_CC = θ_CFC) il n'y a aucune transition solide(1)→solide(2) : le CFC l'emporte à la fois pour l'enthalpie (plus basse, −6 contre −4 ε) et — en raison du nombre de liaisons plus élevé (c_v = 4,5 > 3,5) — pour l'entropie de vibration ; μ_CFC reste à toute température sous μ_CC. Seule une rigidité phononique dépendante du réseau θ_CC < θ_CFC (le CC plus ouvert vibre plus mollement et a donc plus d'entropie) fait passer μ_CC sous μ_CFC à haute température — alors la transformation s'ouvre, et elle se situe sous le pic de fusion τ_m, d'où l'ordre solide(1)→solide(2)→liquide→gaz (τ_ss < τ_m < τ_c ≤ τ_b(3D), chap. 7.11). La rigidité molle du CC est une extension modérée hors de l'image idéale stricte (connue en réalité comme la polymorphie γ/δ du fer : CC à haute température, CFC entre les deux) ; le palier de c_v de 1,0 k_B et la petite chaleur de transformation positive sont en revanche internes au modèle et sans paramètre.

Quelle est vraiment la petitesse de θ_CC/θ_CFC ? (theta_springs_16.py). Afin de libérer la transformation du paramètre libre θ, la rigidité phononique est maintenant estimée microscopiquement (chap. 6.6) : un modèle de ressorts d'Einstein avec ressorts de paire centraux donne θ_CC/θ_CFC ≈ 0,82 (seulement les plus proches voisins) jusqu'à ≈ 0,94 (avec la deuxième couche du CC). Le seuil d'ouverture exigé par la comparaison des μ se situe en revanche à θ_CC/θ_CFC < 0,61 (en dessous, τ_ss existe et se trouve sous τ_m). Le rapport géométrique de rigidité est donc trop grand : dans l'image idéale stricte le CFC reste stable à toute température, il n'y a aucune transition solide(1)→solide(2) purement géométrique. La polymorphie réelle (fer γ↔δ) requiert des contributions supplémentaires (électroniques/magnétiques/anharmoniques) hors de l'image harmonique idéale — précisément l'« extension modérée » posée plus haut. Cela confirme le constat honnête et fait passer le rapport θ d'une hypothèse libre à une limite du modèle estimée et nommée.

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

6.4 Réflexions : l'anomalie acoustique de l'eau liquide

Les liquides normaux montrent un v(T) décroissant de façon monotone (chap. 3.3). L'eau a un maximum acoustique vers 74 °C : v croît d'abord avec T avant de décroître. Dans le modèle, cela se laisse lire comme une concurrence entre deux types de réseau locaux (chap. 6.2) : un ordre local ouvert, faiblement coordonné, semblable à la glace (z_max petit, beaucoup de « trous »), et un ordre local dense, plus fortement coordonné. À basse T la structure ouverte domine ; chauffer détruit d'abord l'ordre local ouvert au profit du plus dense — la connectivité effective et donc K_eff augmentent (v croît). Ce n'est qu'au-dessus du maximum que l'ordinaire rupture de liaison de la FISSURE l'emporte, K_eff chute, v décroît. L'anomalie est ainsi l'expression d'un changement d'ordre local de type solide(1)↔solide(2), estompé et incomplet, à l'intérieur de la phase liquide — concevable seulement dans le modèle 3D à deux types de réseau, non dans le réseau 2D isolé.

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

6.5 Alternative : champs de tourbillons fractals au lieu de la rotation des micelles

Le mouvement brownien local (chap. 3.4) a été expliqué par la rotation et l'échange des micelles. On peut aussi se le représenter comme l'effet de champs de tourbillons fractals : les frontières de blocs sans cesse remaniées par la FISSURE engendrent un champ de courants d'impulsion locaux qui — comme dans les écoulements turbulents — présente des tourbillons autosimilaires à toutes les échelles de longueur. Une particule d'essai est alors mue non par des chocs discrets de micelles, mais par un continuum de tourbillons invariant d'échelle. Les deux images (micelles discrètes ↔ tourbillons fractals) sont des cas limites complémentaires de la même cinétique de FISSURE ; celle qui l'emporte dépend du rapport entre la taille de la particule d'essai et la taille moyenne des micelles. Cette brève discussion est purement qualitative et n'est pas développée mathématiquement ici.

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

6.6 Rigidité phononique microscopique θ(réseau) : un modèle de ressorts d'Einstein

La comparaison solide(1)→solide(2) du chap. 6.3 tient à la rigidité phononique relative θ_CC/θ_CFC. Au lieu de la choisir librement, elle est estimée ici à partir de la coordination et de la géométrie (theta_springs_16.py).

Modèle. Un atome central est dans la cuvette harmonique de ses z voisins ; chaque liaison est un ressort de paire radial κ. Un déplacement u se projette isotropiquement avec ⟨cos²⟩ = 1/d = 1/3 (d = 3) sur les directions de liaison, de sorte que la constante de ressort effective est k_eff = (z/3)·κ. La fréquence d'Einstein est ω = √(k_eff/m), donc θ = ħω/k_B ∝ √k_eff. À masse égale et ressort de paire égal il s'ensuit

θ_CC/θ_CFC = √(z_CC/z_CFC) = √(8/12) = 0,816.

Affinements. (B) La deuxième couche du CC (6 atomes à r₂/r₁ = 2/√3 = 1,155) avec un ressort mis à l'échelle de la distance κ(r) ∝ (r₁/r₂)ⁿ (n = 6…8, de type Lennard-Jones) élève le rapport à 0,91…0,94 (le CC est effectivement raidi par la deuxième couche proche). (C) La liaison plus courte du CC à densité égale (r_nn(CC)/r_nn(CFC) = 0,972, −2,8 %) raidit légèrement le ressort et donne 0,89…0,91. Au total une fenêtre géométrique θ_CC/θ_CFC ≈ 0,82…0,94.

Seuil d'ouverture. La comparaison des μ (chap. 6.3) n'ouvre la transformation (τ_ss existe et se trouve sous τ_m = 1,70) que pour θ_CC/θ_CFC < 0,61. La fenêtre géométrique est nettement au-dessus.

Constat. Le modèle de ressorts pur ne rend pas le CC assez mou : dans l'image idéale stricte le CFC reste stable, et une transition solide(1)→solide(2) d'origine géométrique n'apparaît pas. θ_CC/θ_CFC est ainsi passé d'une hypothèse libre à une grandeur estimée dont la valeur naturelle (≈ 0,8–0,9) maintient la transformation fermée — en cohérence avec le constat honnête du chap. 6.3. La polymorphie réelle se situe donc manifestement hors de l'image harmonique idéale (contributions électroniques/magnétiques, anharmonicité), et non dans la seule géométrie du réseau.

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

6.7 Estimation mathématique de la métastabilité (2D triangulé ↔ carré)

Pour la discussion qualitative du chap. 3.2, voici le calcul d'ordre de grandeur correspondant dans l'image 2D (disques durs, diamètre D, distance de contact a = D).

(1) Énergie. Par particule U/N = −(ε/2)·z (chaque liaison pour moitié) : triangulé z = 6 → −3 ε, carré z = 4 → −2 ε. La transformation triangulé → carré coûte donc ΔU = +1 ε par particule (endothermique ; la forme lâche est plus riche en énergie).

(2) Volume (surface) et travail. Surface par particule : triangulé A_tri = (√3/2)·D² ≈ 0,866 D², carré A_quad = D². La dilatation est ΔA = A_quad − A_tri ≈ 0,134 D² (+15,5 %) — l'augmentation de volume montrée dans la bande dessinée. Contre une pression (réduite) p, il faut fournir le travail W = p·ΔA ; l'enthalpie libre croît de Δμ = ΔU + p·ΔA − τ·ΔS. Le terme d'entropie ΔS > 0 (la forme lâche a des modes plus mous, un palier de c_v plus grand Δc_v = (6−4)/4 = 0,5 k_B) n'abaisse μ_carré sous μ_tri qu'à haute τ — c'est pourquoi le triangulé est l'état fondamental de basse température et le carré la forme de haute température / métastable.

(3) Barrière de germination (2D). La métastabilité présuppose une barrière. En 2D un germe de la nouvelle phase (rayon r) croît avec un gain de surface −Δμ_A·π r² (Δμ_A = force motrice par surface) contre la contribution de bord +λ·2π r (λ = énergie d'arête par longueur) : ΔG(r) = −Δμ_A·π r² + 2π λ r. Maximum en r\* = λ / Δμ_A avec la barrière ΔG\* = π λ² / Δμ_A. Le germe critique r\ se rétrécit à mesure que la force motrice Δμ_A croît ; tant que l'échelle thermique τ = k_BT n'atteint pas la barrière (**ΔG\ ≳ quelques k_BT), la phase lâche (ou, lors de l'ébullition, la phase sursaturée) reste cinétiquement figée. Pour de petits amas (2–3 particules, comme dans la bande sur la sursaturation), le rapport bord/surface est maximalement défavorable — ils se trouvent à gauche de r\* et se dissolvent en moyenne. Ainsi l'énergie d'activation montrée dans la bande dessinée est la même grandeur ΔG\* qui porte à la fois la forme cristalline métastable et la phase gazeuse sursaturée ; c'est un apport cinétique** au modèle d'équilibre idéal (sans barrière).

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

6.8 Affinement de type Debye de θ(réseau)

Le modèle de ressorts d'Einstein (chap. 6.6) remplace toute la distribution des phonons par une seule fréquence. L'approche présente affine θ(réseau) par une densité d'états de phonons épurée par réseau (theta_debye_17.py) : les mêmes ressorts de paire centraux κ(r) = r^(−n) sur la première et la deuxième couche de voisins forment la matrice dynamique D(k) = Σ_j κ_j·(R̂_j⊗R̂_j)·(1 − cos(k·R_j)), diagonalisée sur une grille k de la zone de Brillouin primitive (Monkhorst–Pack, 12³ points) au moyen d'un solveur analytique 3×3, bibliothèque standard uniquement. Des trois branches de phonons découlent deux mesures de rigidité de type Debye comme moments de fréquence :

• 2ᵉ moment θ(2) ∝ √⟨ω²⟩ — la limite d'Einstein (la moyenne sur la ZB de 1 − cos vaut 1, donc ⟨ω²⟩ = z·κ/3) ; contrôle, plus proches voisins seulement, géométrie pure : √⟨ω²⟩_CC / √⟨ω²⟩_CFC = 0,816 = √(8/12) — exactement la valeur d'Einstein du chap. 6.6.
• moment logarithmique θ(0) ∝ exp⟨ln ω⟩ — la grandeur véritablement de type Debye, qui pondère les branches acoustiques molles que la valeur d'Einstein ignore.

Toutes les constantes portant une échelle (ħ/k_B, l'échelle de κ, le préfacteur de Debye) s'annulent dans le rapport θ_CC/θ_CFC ; seule la forme de la DOS subsiste. Avec la deuxième couche de voisins (qui seule stabilise mécaniquement le CC) et une loi de ressort n = 0…8 on obtient

θ_CC/θ_CFC ≈ 0,87…1,00 (bande d'erreur sur la définition du moment et la loi de ressort).

C'est plutôt plus élevé que la valeur d'Einstein pure (0,82…0,94), parce que la liaison plus courte du CC à densité égale raidit les ressorts. Le seuil d'ouverture 0,61 exigé par la comparaison des μ (chap. 6.3) reste nettement en dessous : la DOS des phonons non plus ne rend pas le CC assez mou. Le constat du chap. 6.6 est ainsi confirmé et précisé par l'approximation plus fine — dans l'image idéale stricte le CFC reste stable à toute température ; aucune transition solide(1)→solide(2) d'origine purement géométrique n'apparaît, et la polymorphie réelle (fer, chap. 7.13) exige des contributions hors de l'image harmonique idéale.

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

6.9 θ(réseau) à partir d'un réseau de Lennard-Jones complet

Le chap. 6.8 utilise un ressort central tronqué κ(r) = r^(−n) sur nn + 2nn. La version 17 examine la même question avec une véritable constante de force issue d'un potentiel de paire sur le réseau complet (theta_ljlattice_18.py) : chaque liaison devient un ressort tiré du potentiel de Lennard-Jones φ(r) = 4ε[(σ/r)¹² − (σ/r)⁶], avec la part longitudinale κ_L = φ″(r) et la part transversale κ_T = φ′(r)/r. La matrice dynamique

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

somme sur toutes les couches de voisins jusqu'à un rayon de coupure (2,0…3,0 nn) ; chaque réseau est au préalable relaxé vers son propre équilibre LJ (constante de réseau libre, densité égale non imposée). C'est le test honnête permettant de savoir si l'approximation épurée nn+2nn (chap. 6.8) situe correctement l'ordre des rigidités.

Résultat (déterministe, bibliothèque standard uniquement) :

Deux constats, tous deux renforçant les chap. 6.6/6.8 :

1. Le CFC est l'état fondamental LJ stable (E ≈ −8,4 ε par atome à la coupure 3,0 nn, contre −8,0 ε pour le CC) — la primauté bien connue de l'empilement dense dans un potentiel de paire.
2. Le CC LJ idéal est instable au cisaillement : environ 3…3,5 % des modes de phonons sont imaginaires (ω² < 0, instabilité de Born). L'image idéale du potentiel de paire pur ne porte donc aucun CC mécaniquement stable ; sur les modes réels restants θ_CC/θ_CFC ≈ 1,0 — le seuil d'ouverture 0,61 exigé par la comparaison des μ n'est à plus forte raison pas franchi.

L'affirmation des chap. 6.6/6.8 est ainsi robuste face à l'affinement du ressort tronqué vers le réseau complet et du modèle purement longitudinal vers la matrice dynamique complète : dans l'image idéale classique, aucune transition solide(1)→solide(2) purement géométrique n'apparaît — la polymorphie réelle (fer chap. 7.13, titane chap. 7.15) exige des contributions hors du potentiel de paire (électroniques/magnétiques, anharmoniques). Un véritable calcul de phonons ab initio est délibérément laissé de côté (il quitterait le caractère épuré de l'image idéale).

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

7 — Calibration et équation d'état

Résumé. Le chapitre 7 calibre le modèle sur des substances réelles au moyen du seul paramètre d'énergie ε et développe la branche liquide/gaz en une équation d'état :

• 7.1 / 7.2 Méthode et tableau de calibration : couplage linéaire de température T = (ε/k_B)·τ avec un point d'ancrage au point de fusion ; d'où le rapport universel affiné R = τ_b/τ_m ≈ 1,54.
• 7.3 / 7.4 Équation d'état van der Waals à cœur dur : l'EOS réduite p_r = 8τ_r/(3v_r − 1) − 3/v_r² avec le point critique (1, 1, 1) et la coexistence de Maxwell ; a ∝ ε, b ∝ D². Les points d'ébullition se situent systématiquement à T_b/T_c ≈ 0,48 (écart de Guldberg) ; un calcul de comparaison Redlich–Kwong les relève à ≈ 0,60 — comme budget d'erreur nommé, sans changement d'EOS.
• 7.5 Surface (τ, v) commune : la coupole de coexistence (ébullition), la ligne de fusion et la paroi dure s'assemblent avec un point triple (τ_r,triple ≈ 0,47 pour v_solid = 0,40).
• 7.6 Chaleur latente : L ≈ 1,41 ε donne ΔS_fus ≈ 1,175 R (proche de la règle de Richards) ; les métaux sont excellemment atteints (cuivre ×1,00), tandis que les cristaux ioniques, à liaisons hydrogène et moléculaires sortent de l'image idéale.
• 7.7 Unification ε ↔ ε_lg : k_B T_c = z_max·ε/4 couple les deux échelles (ε_lg = (27/32)·z_max·ε) ; la prédiction 3D τ_c/τ_m = 2,50 atteint les substances moléculaires simples, la 2D restant un cas particulier.
• 7.8 / 7.9 Point triple sans paramètre : v_solid découle du seuil de rigidité z̄_c = 2d, non de la connectivité ; en 3D s'ouvre une fenêtre liquide nette (τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00), en 2D elle reste fermée (rôle particulier de la 2D, désormais triplement attesté).
• 7.10 / 7.11 Véritables balayages 3D : la percolation de liens sur CFC/CC donne z̄_c(3D) ≈ 1,5–1,6 ; le balayage en τ de la cinétique de liaison montre τ_m ≈ τ_rig (fusion = rigidité, bas) et τ_b(3D) ≈ τ_c (ébullition = désagrégation de la connectivité au point critique).
• 7.12 / 7.13 / 7.15 Ligne solide(1)→solide(2) : du point τ_ss la pression fait une ligne (dp/dτ = ΔS/Δv > 0, Δv = +8,9 %), portée comme quatrième ligne frontière dans la surface (τ, v) et calibrée comme ligne (τ, p) sur le fer (γ→δ), confirmée sur le titane (α→β) et contrôlée avec l'argon comme contrôle négatif.
• 7.14 Branche Redlich–Kwong : comme variante autonome de l'image idéale, avec sa propre calibration, à côté du choix par défaut van der Waals (T_b/T_c ≈ 0,60, erreur −1 %).

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

7.1 Méthode de calibration (un paramètre d'énergie ε) — ajustée

Linéaire : T = (ε/k_B)·τ, un point d'ancrage par substance. Ancrage au point de fusion : ε/k_B = T_m,exp/τ_m, prédiction T_b,pred = T_m,exp·(τ_b/τ_m). Avec le pic de fusion affiné (terme coopératif), le maximum de c_v se situe maintenant à τ_m = 1,20 (au lieu de 1,15), la percolation de connectivité (proche de l'ébullition) restant inchangée à τ_b = 1,85. D'où le rapport universel affiné :

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

Le tableau de calibration de la version 08 (chap. 7.2) se met à l'échelle linéairement avec R : tous les T_b,pred baissent du facteur 1,54/1,61 ≈ 0,957 (≈ −4,3 %). Cela améliore légèrement la surestimation des gaz rares (argon +54 % → +48 %) et détériore un peu les résultats pour NaCl et le fer (NaCl −0,6 % → −4,8 %). L'énoncé de base demeure : R convient aux substances à empilement dense isotrope et à liaisons non directionnelles ; gaz rares et chaînes sortent de l'image idéale.

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

7.2 Tableau de calibration

Repris inchangé de modell_08.md chap. 7.2 (avec la remarque sur la mise à l'échelle linéaire en R du 7.1). Reproductible via cv_scan_08.py (tableau) et cv_scan_09.py (τ_m, τ_b, R).

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

7.3 Axe densité/volume : équation d'état van der Waals à cœur dur

La pression passe du statut de « paramètre » à celui de véritable variable. Le CŒUR DUR fournit le volume propre b, l'attraction de liaison (ε·z̄) la pression de cohésion a. En unités réduites (point critique en (v_r, τ_r, p_r) = (1, 1, 1)) :

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

Coexistence par construction de Maxwell (cv_scan_09.py, coupe à aires égales) :

Résultat. En dessous de τ_c s'ouvre le saut de densité l/g (v_gas/v_liq croît de 2,5 à 11,5, d'autant plus grand qu'il fait froid). C'est la ligne d'ébullition — désormais nettement séparée du pic de c_v de fusion (chap. 5.3), qui vit sur un réseau de densité fixe. Ainsi le point 3 laissé ouvert dans la version 08 est clos : fusion (ordre de rigidité/de liaison à densité fixe) et ébullition (saut de densité/volume) sont des axes différents du même modèle.

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

7.4 Calibration quantitative l/g : a → ε, b → D²

Application des grandeurs réduites à des substances réelles. Les constantes de van der Waals tiennent au point critique :

b = k_B·T_c / (8 p_c) (volume propre, CŒUR DUR), a = 27 k_B²·T_c² / (64 p_c) (cohésion, LIAISON).

Ainsi a ∝ ε (la cohésion provient de l'énergie de liaison) et b ∝ D² (le volume propre 2D issu du diamètre du cœur dur). L'échelle d'énergie naturelle de la branche l/g est

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

Substance d'ancrage argon (T_c = 150,87 K, p_c = 48,63 bar) : b = 5,35·10⁻²⁹ m³ par molécule, a = 3,76·10⁻⁴⁹ J·m³ par molécule², ε_lg/k_B = 509 K et un diamètre moléculaire D ≈ b^{1/3} = 3,8·10⁻¹⁰ m — à 11 % près du diamètre van der Waals de l'argon (~3,4·10⁻¹⁰ m).

Prédiction de la température d'ébullition. À pression extérieure fixée, le point d'ébullition est là où la pression de saturation réduite p_sat(τ_r)·p_c = p_atm est satisfaite ; d'où T_b = T_c·τ_r,b (courbe de saturation de Maxwell issue de lg_calib_10.py) :

*Le CO₂ sublime à 1 atm (T = 194,7 K = point de sublimation).

Résultat. Pour les substances moléculaires (p_c ≈ 30–70 bar) le modèle donne partout T_b/T_c ≈ 0,48–0,51, systématiquement ~17–24 % en dessous de la mesure. Ce n'est pas une erreur du modèle au sens strict, mais l'écart bien connu de l'équation de van der Waals à la règle empirique de Guldberg (T_b/T_c ≈ 0,6) : l'isotherme de van der Waals sous-estime la pression de saturation à basse température. L'eau (et les métaux, dont les T_c/p_c ne sont de toute façon qu'estimés grossièrement) se situent à p_r,atm < 0,005 hors du domaine van der Waals sensé (τ_r,b < 0,40) — là, la cohésion n'est plus de type van der Waals (liaisons hydrogène, liaison métallique). L'ordre de grandeur et le classement des points d'ébullition sont correctement atteints.

Classement de Guldberg comme tableau d'erreurs. L'écart est ici seulement nommé, non corrigé par un changement d'EOS (l'image idéale reste de type van der Waals). La règle empirique de Guldberg (Guldberg 1890) énonce que le point d'ébullition se situe vers T_b/T_c ≈ 0,6 (avec une légère dépendance à la substance, ~0,58–0,64). Le modèle van der Waals réduit donne à la place une valeur constante T_b/T_c = τ_r,b ≈ 0,48–0,51 :

Lecture. La moyenne du modèle, 0,49, atteint la constante de Guldberg 0,60 à environ −18 % près ; l'écart est systématique et unilatéral (toujours trop bas), donc une propriété de l'isotherme de van der Waals (pression de saturation trop faible aux τ bas), non de l'image moléculaire. Une meilleure EOS de cohésion (par exemple Redlich–Kwong) relèverait τ_r,b — mais reste ici délibérément hors de l'image idéale.

Budget d'erreur avec un calcul de comparaison Redlich–Kwong (guldberg_rk_15.py). Afin de chiffrer l'écart van der Waals à la règle de Guldberg (et non de le guérir), le même test du point d'ébullition est repris avec une cohésion de type Redlich–Kwong, dépendante de la température, a/√T. En unités réduites (point critique (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,

avec la construction de Maxwell à aires égales et la même pression extérieure réduite p_r,atm que ci-dessus. Résultat (une valeur par substance d'ancrage, T_b/T_c = τ_r,b) :

Résultat. La cohésion RK dépendante de la température relève T_b/T_c de vdW ≈ 0,48 à ≈ 0,60 — pratiquement la constante de Guldberg et les valeurs réelles (≈ 0,58–0,63) ; l'erreur moyenne face aux mesures tombe de −21 % (vdW) à −1 % (RK). La cause de l'écart de Guldberg est ainsi clairement chiffrée : c'est justement l'absence de dépendance en température de la cohésion de van der Waals (pression de saturation trop faible aux τ bas), et non l'image moléculaire. Ce calcul de comparaison se situe délibérément hors de l'image idéale ; le modèle reste de type van der Waals (aucun changement d'EOS) — le budget d'erreur est seulement nommé, non incorporé.

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

7.5 Surface (τ, v) commune avec point triple

Les axes de fusion et d'ébullition vivaient jusqu'ici séparément : la fusion (ordre de rigidité/liaison) à densité fixe, l'ébullition (van der Waals) par le saut de densité. Dans un plan commun (τ_r, v_r) ils s'assemblent :

• Coupole de coexistence (ébullition) : v_liq(τ_r) (branche gauche) et v_gas(τ_r) (branche droite) issues de la construction de Maxwell, qui se rejoignent au point critique (τ_r, v_r) = (1, 1).
• Solide (fusion) : bande étroite juste au-dessus de la paroi dure v_r > 1/3, à densité élevée et quasi fixe (volume réduit v_solid ≈ 0,40).
• Ligne de fusion : quasi verticale en v_solid (densité fixe), à la température de la rigidité de Maxwell (z̄_c = 2d ; le chap. 7.8 en déduit v_solid sans paramètre — l'ancienne attribution de la version 11, « z̄_c ≈ 2,08 », est ainsi corrigée).

Le point triple est le point où la frontière de coexistence liquide v_liq(τ_r) atteint la densité d'empilement solide v_solid — en dessous il n'existe plus de branche liquide (sublimation directe) :

Résultat. Pour v_solid = 0,40 le point triple se situe à τ_r,triple ≈ 0,47, p_sat ≈ 0,018, avec un volume de vapeur v_gas ≈ 68 (vapeur de sublimation très diluée). Le rapport τ_triple/τ_c ≈ 0,47 se situe dans la fenêtre réelle (eau 0,42 ; argon 0,56 ; CO₂ 0,71). Plus le solide est dense (v_solid plus petit), plus le point triple réduit est bas et plus la vapeur coexistante est diluée — qualitativement correct. Les deux transitions et le point triple sont ainsi ancrés dans une seule surface (τ, v) ; le point critique siège en (1, 1, 1). La version 15 ajoute au chap. 7.12 une quatrième ligne frontière — la transformation solide(1)/solide(2) (CFC↔CC) sous forme de deux bandes de densité solide quasi verticales (v ≈ 0,667 et 0,726) qui basculent l'une dans l'autre à τ_ss.

Version antérieure en cours de développement : v_solid = 0,40 était encore imposé ici. Le chap. 7.8 déduit v_solid sans paramètre du seuil de percolation z̄_c et fait du point triple une véritable prédiction.

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

7.6 Chaleur latente face aux enthalpies de fusion

La chaleur latente du modèle L ≈ 1,41 ε (chap. 5.3) est traduite en unités réelles par la calibration au point de fusion ε/k_B = T_m/τ_m (τ_m = 1,20) :

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

C'est — de façon remarquable — presque exactement la règle de Richards (ΔS_fus ≈ R pour les solides de construction simple). Comparaison avec des enthalpies de fusion mesurées (lg_calib_10.py) :

Résultat. Le modèle atteint partout l'ordre de grandeur (kJ/mol) et excellemment la règle de Richards pour les métaux simples (cuivre ×1,00, fer ×1,28) ; les gaz liquéfiés se situent juste en dessous (méthane ×0,94, azote ×0,87, argon ×0,69). Les cristaux ioniques, à liaisons hydrogène et moléculaires (NaCl ×0,37, eau ×0,44, benzène ×0,28, CO₂ ×0,23) ont des entropies de fusion réelles nettement plus élevées (ΔS_exp/R = 2,6–5,0), parce que la fusion libère des degrés d'ordre supplémentaires (appariement ionique, orientation de molécules polyatomiques) que le modèle à un seul paramètre d'ordre ne saisit pas. La limite est ainsi clairement nommée : l'image idéale vaut pour les substances à liaisons isotropes et non directionnelles (sans ajustement).

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

7.7 Unification des échelles d'énergie : ε ↔ ε_lg

Jusqu'ici la branche de fusion portait l'énergie de liaison ε (τ_m = 1,20), la branche l/g une échelle de cohésion séparée ε_lg/k_B = 3,375·T_c (chap. 7.4). Les deux se laissent déduire d'un seul paramètre du modèle. Dans le modèle, la cohésion de van der Waals a naît de la somme des énergies de liaison par molécule (LIAISON) : chaque molécule porte z_max/2 liaisons d'énergie ε. Un gaz sur réseau de Bragg–Williams en champ moyen avec attraction de voisinage ε fournit la relation classique pour la température critique :

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

En reportant ceci dans ε_lg = a/b = (27/8)·k_B·T_c (chap. 7.4), on obtient le couplage des deux échelles :

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

Ainsi le pic de fusion, la percolation et le point critique se situent sur une seule échelle ε (unify_scale_12.py) :

Prédiction τ_c/τ_m = z_max/(4·τ_m) : 2D → 1,25, 3D → 2,50. Face aux substances réelles :

Résultat. La prédiction 3D τ_c/τ_m = 2,50 atteint bien la valeur caractéristique réelle des substances moléculaires simples (T_c/T_m ≈ 2,0–2,1) ; les métaux se situent plus haut (T_c estimé grossièrement, autre type de liaison). La prédiction 2D (1,25) est trop basse — et τ_c(2D) = 1,50 se situe même sous τ_b = 1,85 (le point critique serait sous le seuil de percolation). Ce n'est pas une erreur de calcul, mais le rôle particulier bien connu de la 2D : les échelles de cohésion et de liaison ne s'y ordonnent pas proprement, une fenêtre liquide robuste manque (cf. la fusion 2D continue, chap. 5.5). Ce n'est qu'en 3D (τ_m = 1,20 < τ_c = 3,00) que les échelles se séparent de façon physiquement sensée.

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

7.8 Point triple sans paramètre à partir de z̄_c

Le chap. 7.5 imposait encore v_solid = 0,40. La version antérieure en cours de développement déduit la densité d'empilement solide du seuil de percolation. Le nombre de coordination d'un réseau à volume réduit v_r suit la relation densité–coordination (paroi dure v_r = 1/3 ↔ empilement le plus dense, z̄ = z_max) :

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

Le solide perd sa rigidité dès que z̄ passe sous le seuil z̄_c ; la densité correspondante est v_solid. Le point triple est alors le point où la frontière de coexistence liquide v_liq(τ_r) atteint cette densité (unify_scale_12.py, 2D, z_max = 6) :

Version améliorée. La densité de fusion v_solid est fixée par la rigidité générique de Maxwell z̄_c = 2d (2D → 4, 3D → 6), non par la connectivité de liens (justification au chap. 5.4a : la fusion est la perte de la rigidité, non de la connexité). Le point triple fiable s'écrit alors (rigidity_3d_13.py, 2D) :

Résultat 2D. Même avec le seuil de rigidité (physiquement correct), le réseau 2D isolé ne donne aucune fenêtre liquide nette : v_solid = 0,50 fournit bien la valeur réelle τ_triple/τ_c ≈ 0,77, mais en unités ε absolues τ_triple = 0,77·1,50 = 1,16 et se situe donc sous τ_m = 1,20. Le seuil de connectivité (z̄_c = 2,08) pousse même v_solid ≈ 0,96 ≈ v_c hors de la coupole de coexistence — point triple et point critique coïncident. Les deux calculs confirment le rôle particulier de la 2D (chap. 5.5, 7.7) : le réseau 2D isolé ne porte guère de phase liquide autonome. C'est un constat honnête, non une faiblesse de calcul ou de modèle au sens strict.

Le test décisif est donc le calcul de contrôle 3D (chap. 7.9) : là, la même règle de rigidité (z̄_c = 6 → v_solid = 0,667) libère une fenêtre nette. Dans tous les cas, le point triple est désormais lié sans paramètre à z̄_c — le paramètre libre v_solid est éliminé.

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

7.9 Calcul de contrôle 3D : une fenêtre liquide nette sur une seule échelle ε

Vérification que l'effondrement 2D se résout en 3D lorsque tous les repères se situent sur une seule échelle ε (τ_c = z_max/4 du chap. 7.7 ; v_solid issu de la rigidité générique z̄_c = 2d ; point triple là où la densité liquide réduite de van der Waals v_liq(τ_r) atteint la densité du solide). Résultat (rigidity_3d_13.py) :

Résultat. En 3D s'ouvre une fenêtre liquide nette : τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00 (largeurs : τ_triple − τ_m = 1,63 ; τ_c − τ_triple = 0,18). Le bord supérieur est étroit — le point triple se trouve juste sous le critique — mais l'ordre est correct et la phase liquide existe comme domaine autonome. En 2D la fenêtre reste fermée (τ_triple = 1,16 < τ_m). Le rôle particulier de la 2D est ainsi attesté pour la troisième fois de façon indépendante et il est en même temps montré que c'est un effet de dimension : la même règle de rigidité guérit le modèle en 3D. Le bord supérieur étroit est consolidé au chap. 7.10 par un véritable balayage de percolation 3D et une analyse de sensibilité.

Encadré des limites du modèle — le rôle particulier de la 2D. Trois calculs indépendants montrent le même constat : le réseau 2D isolé ne porte guère de phase liquide autonome. 1. τ_c < τ_b (chap. 7.7) : le point critique (1,50) serait sous la percolation de connectivité (1,85) — les échelles de cohésion et de liaison ne s'ordonnent pas proprement en 2D. 2. v_solid ≈ v_c (chap. 7.8) : la densité de fusion de connectivité (0,96) coïncide presque avec la densité critique — point triple et point critique s'effondrent l'un sur l'autre. 3. τ_triple < τ_m (chap. 7.9) : même avec la rigidité de Maxwell correcte, le point triple (1,16) est sous le pic de fusion (1,20) — pas de fenêtre.

C'est un constat honnête, non une faiblesse de calcul ou de modèle : c'est l'effet de dimension bien connu (fusion 2D continue, de type KTHNY ; chap. 5.5). En 3D la même règle de rigidité guérit le modèle (chap. 7.9/7.10). Qu'un second ordre local (chap. 6.4) puisse encore améliorer l'image 2D reste une question d'extension, non un défaut de l'image idéale. Comme le modèle doit servir à l'illustration, j'accepte cette faiblesse et renonce à une extension axiomatique du modèle.

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

7.10 Véritable balayage de percolation de liens en 3D (CFC/CC) et sensibilité du bord de la fenêtre

Les chap. 7.8/7.9 transposaient en pensée la valeur de connectivité 2D z̄_c = 2,08 à la 3D (où v_solid = 1,92, non physique) et appuyaient l'énoncé 3D sur la rigidité générique z̄_c = 2d. La version 13 examine maintenant la connectivité en 3D directement : percolation de liens sur CFC (z_max = 12) et CC (z_max = 8) par Monte-Carlo avec union-find, bords périodiques, graine LCG déterministe (rigidity_3d_14.py). Seuil p_c = croisement P_∞(p) = 0,5 (même méthodologie que l'étude de taille finie 2D du chap. 5.6) :

En raison de l'arête de boîte modérée (L = 12, Python pur), les valeurs MC se situent légèrement au-dessus des seuils de la littérature (biais de taille finie, comme en 2D par en dessous ou au-dessus selon l'estimateur) ; le classement et l'ordre de grandeur sont corrects (CFC plus dense → p_c plus petit que le CC).

Lecture (confirmation 3D indépendante). Les seuils de connectivité 3D se situent à z̄_c(3D) ≈ 1,5–1,6 — encore plus bas que la valeur 2D 2,08 et bien en dessous de la rigidité de Maxwell 2d = 6. La « densité de fusion de connectivité » qui en découle, v_solid = 1/(3 p_c) ≈ 1,7 (CC) ou 2,6 (CFC), se situe bien hors de la coupole de coexistence (v_liq,max ≈ 0,94). Il est ainsi montré directement — sans transposer la valeur 2D à la 3D — qu'en 3D la connectivité peut à plus forte raison ne pas fixer la densité de fusion ; celle-ci doit venir de la rigidité (z̄_c = 6 → v_solid = 0,667). Les deux seuils sont encore plus nettement séparés en 3D qu'en 2D. Le véritable τ_b(3D) (température à laquelle z̄(τ) atteint le seuil bas de connectivité 1,5) se situerait donc au-dessus de la valeur 2D τ_b, tout près du point critique — en cohérence avec l'image « désagrégation de la connectivité = ébullition ».

Sensibilité du bord supérieur étroit de la fenêtre (τ_c − τ_triple ≈ 0,18). Numériquement autour de v_solid = 0,667 (±2 %) on obtient dτ_triple/dv_solid ≈ 1,55. La fenêtre ne se fermerait (τ_triple = τ_c = 3,00) que si v_solid croissait de ≈ 0,11 (jusqu'à ≈ 0,78), c'est-à-dire si le seuil de rigidité 3D effectif se situait sensiblement au-dessus de 2d = 6 (z̄_c ≳ 7,7). Pour le comptage générique de Maxwell (z̄_c = 2d = 6), l'ordre τ_m < τ_triple < τ_c reste donc robuste ; le bord supérieur est réellement étroit, mais non menacé. Le changement discret de réseau solide(1)→solide(2) (CFC↔CC) est quantifié dans la version 14 (chap. 6.3) ; le véritable balayage en τ de la cinétique de liaison en 3D suit au chap. 7.11.

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

7.11 Véritable balayage en τ de la cinétique de liaison en 3D : τ_b(3D) et la ligne de fusion sur une échelle ε

Les chap. 7.8–7.10 appuyaient l'énoncé 3D sur la rigidité générique (z̄_c = 2d = 6) et le balayage statique de connectivité (z̄_c(3D) ≈ 1,55). Restait ouvert le véritable balayage en τ de la cinétique de liaison en 3D : où, sur l'axe des températures, le z̄(τ) auto-cohérent atteint-il ce seuil bas de connectivité (= le véritable τ_b(3D)), et où se situe le pic de fusion ? La version 14 calcule le z̄(τ) auto-cohérent (chap. 4.2, x = 1/(1 + G·e^{−(ε+K·x)/τ})) directement sur CFC (z_max = 12) et CC (z_max = 8) avec G = z_max et la même coopérativité J = 0,55 que dans le cas de référence 2D ; n_blk(τ) provient d'un Monte-Carlo de liens CFC/CC (union-find), et

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄ − (K/2)·x² (3D, avec terme coopératif),

c_v = dU/dτ (bondkinetics_3d_15.py, graine LCG déterministe ; n_blk lissé, puisque seul le faible terme de rotation de bloc en dépend). Trois repères sont relevés : le pic de c_v τ_m (fusion), le repère de rigidité τ_rig (z̄(τ) croise 2d = 6) et le véritable τ_b(3D) (z̄(τ) croise le seuil de connectivité z̄_c(3D) du chap. 7.10) :

Résultat. Deux constats, tous deux cohérents avec l'image jusqu'ici :

1. La fusion se situe bas et est doublement attestée. Le pic thermodynamique de fusion τ_m (CFC 1,70, CC 1,40) coïncide pratiquement avec le repère de rigidité τ_rig (z̄ = 6 ; CFC 1,73, CC 1,35) — la même température à laquelle le réseau perd sa rigidité. Cela ancre la ligne de fusion 3D sur l'échelle ε : T_m = (ε/k_B)·τ_m, soit par exemple pour l'argon (T_m = 83,8 K) ε/k_B = 83,8/1,70 ≈ 49 K. Comme K = J·z_max est plus grand en 3D qu'en 2D (6,60 au lieu de 3,30), le pic de fusion 3D est plus élevé que la valeur 2D 1,20 — pas de contradiction, mais la coopérativité plus forte de l'empilement plus dense.
2. Le véritable τ_b(3D) siège au point critique. Parce que le seuil de connectivité 3D est si bas (≈ 1,55), z̄(τ) ne l'atteint que tout en haut : τ_b(3D) = 3,18 (CFC) et 2,85 (CC) — au niveau, ou légèrement au-dessus de τ_c = z_max/4 (3,00 et 2,00). L'amas géant de liaisons survit donc à toute la phase liquide et ne se désagrège qu'au point critique. C'est exactement la prédiction du chap. 7.10 (« le véritable τ_b(3D) se situerait au-dessus de la valeur 2D, tout près du point critique ») et cela confirme de façon indépendante l'image « ébullition = désagrégation de la connectivité » : en 3D la matière reste liée d'un bout à l'autre jusqu'à tout près du point critique.

Les deux constats consolident la séparation des seuils du chap. 5.4a : fusion ↔ rigidité (bas, τ_m ≈ τ_rig), ébullition ↔ connectivité (haut, τ_b(3D) ≈ τ_c). La densité de fusion vient donc de la rigidité (v_solid = 0,667), non de la connectivité — les deux seuils sont écartés sur l'axe des températures en 3D plus qu'ils ne l'ont jamais été dans le modèle.

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

7.12 La ligne de transformation solide(1)→solide(2) dans la surface (τ, v)

Le chap. 6.3 fournissait la transformation de réseau CFC↔CC comme point τ_ss (à p = 0). Par la pression il en résulte une ligne (polymorph_line_16.py), qui a sa place comme quatrième ligne frontière dans la surface (τ, v) commune des chap. 7.5/7.9 — à côté de la coupole de coexistence (ébullition), de la ligne de fusion et de la paroi dure.

Saut de volume. Les deux solides siègent à densité quasi fixe ; leur différence de volume découle des compacités (v ∝ 1/η) : η_CFC = 0,7405, η_CC = 0,6802. Avec v_CFC = 0,667 (chap. 7.9) on obtient v_CC = v_CFC·η_CFC/η_CC = 0,726, donc

Δv = v_CC − v_CFC = +0,059 (+8,9 %).

Le CC plus ouvert a le plus grand volume — la recristallisation CFC→CC s'accompagne d'une augmentation de volume (le même énoncé que dans la bande dessinée « recristallisation »).

Clausius–Clapeyron. Avec ΔS = S_CC − S_CFC (phase de haute moins phase de basse température), la pente s'écrit

dp/dτ = ΔS/Δv.

Les deux grandeurs sont positives (le CC est la phase de haute température, il a donc plus d'entropie, et le plus grand volume), donc dp/dτ > 0 : une pression plus élevée relève la température de transformation, parce qu'elle favorise le CFC plus dense. Pour la référence θ_CC/θ_CFC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, L = 0,667 ε) :

Ligne en (τ, p) pour la référence (croissante) : p = 0 → τ_ss = 1,33 ; p = 2 → 1,57 ; p = 4 → 1,81 ; p = 6 → 2,04 (unités réduites).

Dans la surface (τ, v) la transformation apparaît comme deux bandes de densité solide quasi verticales — solide(1) = CFC à v ≈ 0,667 (plus dense, basse T) et solide(2) = CC à v ≈ 0,726 (plus ouvert, haute T) —, qui basculent l'une dans l'autre à τ_ss ; au-dessus suit la fusion à τ_m. L'ordre complet sur l'axe τ est donc τ_ss < τ_m < τ_c ≤ τ_b(3D) (chap. 7.11).

Unités réelles. Par la calibration au point de fusion ε/k_B = T_m/τ_m (τ_m = 1,70, CFC 3D), le palier de c_v et la chaleur de transformation deviennent tangibles. Le palier de c_v est indépendant de la calibration : |Δc_v| = 1,0 k_B = 1 R = 8,31 J/(mol·K) = 0,0083 kJ/(mol·K). La chaleur de transformation ΔH_ss = L·R·(ε/k_B) (référence L = 0,667 ε) :

Les chaleurs de transformation sont petites et positives (ordre de 0,2–0,9 kJ/mol) et se situent donc nettement sous les enthalpies de fusion (chap. 7.6, plusieurs kJ/mol) — exactement l'image attendue d'une transformation polymorphique faible avec un petit palier de c_v. (Le rapport θ de 0,60 se situe délibérément hors de l'image idéale stricte ; le chap. 6.6 montre que la géométrie pure donnerait θ ≈ 0,82–0,94 et que la transformation resterait alors fermée — la ligne est la conséquence si-alors d'une phase CC molle.)

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

7.13 Ligne de transformation solide(1)→solide(2) dans le diagramme (τ, p), calibrée sur le fer

Le chap. 7.12 porte la transformation comme ligne (τ, v) ; la version 16 l'ajoute comme ligne dans l'image (τ, p) et la calibre sur un polymorphe CFC↔CC réel (solidsolid_pt_17.py).

Champ de pression issu de l'EOS van der Waals à cœur dur. Aux deux densités solides v_CFC = 0,667 et v_CC = 0,726 (chap. 7.12), l'EOS réduite p_r = 8τ_r/(3v_r − 1) − 3/v_r² (chap. 7.3) donne partout des pressions positives ; au même τ_r le CFC plus dense porte la pression plus élevée (τ_r = 1 : p_r = 1,25 contre 1,10). La compression favorise donc le CFC — en cohérence avec la pente positive de Clausius–Clapeyron dp/dτ = ΔS/Δv > 0.

Ligne et point de transformation. Pour la référence θ_CC/θ_CFC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, dp\/dτ = 8,47 en unités réduites p\ = p·D³/ε), la ligne croît selon p\(τ) = 8,47·(τ − 1,333) : p\ = 0 → τ = 1,33 ; p\ = 2 → 1,57 ; p\ = 4 → 1,81 ; p\ = 6 → 2,04. À **p\_tri ≈ 3,1 la ligne solide-solide rencontre le repère de fusion τ_m = 1,70 — un point de transformation solide(1)-solide(2)-liquide** dans le plan (τ, p) (analogue du point triple).

Calibration sur le fer (ε/k_B = T_m/τ_m = 1811 K / 1,70 = 1065 K). Le fer a deux transitions CFC↔CC de signe opposé — précisément le cas d'essai :

Pour γ→δ (τ_ss = 1667/1065 = 1,565 → θ_CC/θ_CFC ≈ 0,61) le modèle atteint signe et ordre de grandeur : ΔH_ss ≈ 4,2 kJ/mol (en réalité ≈ 0,84 kJ/mol — l'ordre de grandeur kJ/mol, la valeur étant ~5× trop élevée à cause de l'approximation de seuil abrupte), palier de c_v 1 R = 8,31 J/(mol·K) sans paramètre, dT/dp ≈ +245 K/GPa. La transition α→γ en revanche a la phase dense (γ) comme phase de haute température — ce que la seule géométrie d'empilement ne peut fournir (plus dense = plus bas = basse T) ; le fer α est stabilisé magnétiquement. Le fer confirme ainsi les deux faces du constat honnête : la topologie géométrique (γ→δ) est dans l'image idéale, l'anomalie magnétique (α→γ) en dehors — le même « en dehors » que les chap. 6.6/6.8 chiffrent par θ ≈ 0,87…1,00 > 0,61.

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

7.14 Redlich–Kwong comme branche séparée de l'image idéale

L'écart de Guldberg de la branche van der Waals (chap. 7.4, T_b/T_c ≈ 0,48 au lieu de ≈ 0,60) n'y était que nommé. À la demande expresse de l'utilisateur, une version améliorée porte la cohésion de Redlich–Kwong comme variante autonome de l'image idéale, avec sa propre calibration, sa propre coexistence et sa propre section d'erreur (rk_branch_17.py) — la branche van der Waals reste le choix par défaut ; toute la physique restante du modèle (fusion, percolation, solide(1)→solide(2)) est indifférente au choix de l'EOS de cohésion.

Échelle d'énergie propre. De a_RK = Ω_a·R²T_c^{2,5}/p_c et b_RK = Ω_b·R T_c/p_c découle l'énergie de cohésion de la branche :

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

plus grande, parce que la cohésion dépendante de la température a/√T replie le facteur √T_c. Le diamètre issu de b_RK atteint les valeurs réelles à peu près aussi bien que dans la branche vdW (argon D ≈ 3,34·10⁻¹⁰ m).

Coexistence propre (réduite, point critique (1, 1, 1), Z_c = 1/3 ; p_r = 3τ_r/(v_r − β) − A/(√τ_r·v_r(v_r + β)), β = 0,25992, A = 3,84732) :

Section d'erreur propre (point d'ébullition à p_r,atm = p_atm/p_c, T_b/T_c = τ_r,b) :

Résultat. La branche RK relève T_b/T_c de vdW ≈ 0,48 à ≈ 0,60 et atteint ainsi pratiquement la constante de Guldberg et les valeurs réelles (erreur moyenne −1 % contre vdW −21 %). Elle se tient comme deuxième variante de l'image idéale, clairement séparée, à côté du choix par défaut van der Waals ; laquelle des deux branches est « plus juste » est une question de précision exigée, non de structure du modèle — le noyau du modèle 3 (z̄ comme paramètre d'ordre) reste identique dans les deux.

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

7.15 solide(1)→solide(2) sur d'autres polymorphes : titane et contrôle par gaz rare

Le chap. 7.13 calibre la ligne (τ, p) solide(1)→solide(2) sur le fer. La version améliorée assure signe et ordre de grandeur sur un deuxième polymorphe, magnétiquement discret, et ajoute un contrôle négatif (solidsolid_titanium_18.py). Calibration et saut de volume sont comme au chap. 7.13 (ε/k_B = T_m/τ_m, τ_m = 1,70 ; Δv = +8,9 % dense→ouvert d'après la géométrie d'empilement) ; hcp est approché comme le CFC par z = 12 (la différence cfc/hcp est seulement nommée, non résolue).

Le titane, avec α(hcp) → β(cc) à 1155 K, présente exactement la topologie du modèle : la phase dense (hcp, z = 12) est la phase de basse température, et au chauffage elle devient ouverte (cc, z = 8), donc ΔS > 0 et Δv > 0 → dp/dτ > 0 — le même signe que le fer γ→δ, mais sans anomalie magnétique. Deux polymorphes dense→ouvert indépendants en comparaison :

Les deux atteignent le signe (dp/dτ > 0, la pression favorise la phase dense) et l'ordre de grandeur (ΔH ~ kJ/mol, palier de c_v sans paramètre 1 R = 8,31 J/(mol·K), dT/dp ~ 10²·K/GPa). La valeur de ΔH_ss est dans les deux cas systématiquement d'environ un facteur trop élevée — l'approximation de seuil abrupte de la comparaison des μ surestime la chaleur latente de transformation ; c'est une limite du modèle nommée, non une erreur de signe.

L'argon est le contrôle négatif : l'argon solide est CFC et reste CFC jusqu'à la fusion (83,8 K) — à pression normale aucun cc stable n'existe. C'est précisément ce que prédit l'image idéale du potentiel de paire pur : θ_CC/θ_CFC ≈ 0,9…1,0 > 0,61 et un CC LJ instable au cisaillement (chap. 6.9) ⇒ aucune transition solide(1)→solide(2). La contre-épreuve confirme la structure du constat : l'image idéale n'engendre une ligne solide(1)→solide(2) que lorsque des contributions non géométriques (électroniques/magnétiques, comme dans le fer et le titane) stabilisent d'abord le CC. Le titane fournit la preuve positive, l'argon la preuve négative du même cadre honnête des chap. 6.8/6.9.

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

8 — Comparaison bibliographique : les conceptions de la thermodynamique statique

Résumé. Le chapitre 8 situe le modèle 3 parmi les conceptions établies de la thermodynamique statique :

• 8.1 Tableau d'ensemble : onze modèles (gaz parfait, van der Waals, Redlich–Kwong, développement du viriel, Einstein, Debye, Lennard-Jones / sphères dures, théorie cellulaire, Ising / gaz sur réseau, percolation, Flory–Huggins) avec leur(s) phase(s) et leur idée de base statistico-thermodynamique.
• 8.2 Les modèles un à un : un paragraphe chacun avec caractéristiques, application et limites, commençant par le lien vers la source ; les sources librement accessibles sont archivées sous model3/links/, la théorie cellulaire (paywall) étant seulement citée.
• 8.3 Confrontation : le modèle 3 est un modèle de liaison coopérative à deux états avec le paramètre d'ordre z̄ ; son parent le plus proche est Ising / gaz sur réseau, l'ébullition est reprise directement de la percolation, et l'image des phonons (Einstein/Debye) porte la rigidité de fusion. Le modèle 3 agrafe ces trois traditions sur une échelle ε commune.

Ce chapitre situe le modèle 3 parmi les images établies de la thermodynamique d'équilibre (« statique ») des solides, des liquides et des gaz. D'abord un aperçu tabulaire (8.1), puis un bref paragraphe en prose par modèle, avec caractéristiques typiques, domaine d'application et limites — chaque paragraphe commence par le lien vers la source (8.2). Enfin une confrontation sous forme de mots-clés avec le modèle 3 (8.3). Les sources librement accessibles sont archivées sous model3/links/ ; le script de vérification model3/scripts/check_links_16.py maintient la cohérence entre citation et archive par SHA-256.

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

8.1 Tableau d'ensemble

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

8.2 Les modèles un à un

https://de.wikipedia.org/wiki/Kinetische_Gastheorie — Gaz parfait / théorie cinétique des gaz. Particules sans volume propre et sans interaction hormis les chocs élastiques ; la pression naît purement du transfert d'impulsion. Fournit l'équation d'état pV = Nk_BT, la distribution des vitesses de Maxwell–Boltzmann et ⟨E⟩ = (3/2)k_BT. Application : gaz dilués à haute température et basse pression. Limites : pas de condensation, pas de phase liquide ni solide, pas de transition de phase — l'interaction est absente par construction.

https://de.wikipedia.org/wiki/Van-der-Waals-Gleichung — Équation de van der Waals. (p + a/v²)(v − b) = k_BT complète le gaz parfait par un volume d'exclusion b (répulsion) et un terme de cohésion a (attraction). C'est la première équation d'état avec une transition gaz-liquide, un point critique et une construction de Maxwell. Application : image qualitative de la condensation, grandeurs critiques. Limites : quantitativement imprécise près du point critique (exposants de champ moyen), pas de phase solide.

https://de.wikipedia.org/wiki/Redlich-Kwong-Gleichung — Équation d'état de Redlich–Kwong. EOS cubique avec terme d'attraction dépendant de la température a/[√T · v(v + b)]. La mise à l'échelle en √T améliore nettement les équilibres vapeur-liquide par rapport à van der Waals (cf. chap. 7.4, où cela relève le rapport T_b/T_c). Application : calculs techniques de phase gazeuse et d'évaporation. Limites : paramétrée empiriquement, toujours pas de phase solide ; près du point critique également de type champ moyen.

https://de.wikipedia.org/wiki/Virialentwicklung — Développement du viriel. Développement systématique pV/Nk_BT = 1 + B(T)/v + C(T)/v² + … Les coefficients du viriel découlent du potentiel intermoléculaire (intégrales d'amas de la mécanique statistique) et sont donc microscopiquement fondés. Application : gaz jusqu'à densité moyenne, référence pour d'autres EOS. Limites : la série ne converge pas pour les liquides denses ou les solides ; les coefficients d'ordre supérieur sont coûteux à calculer.

https://de.wikipedia.org/wiki/Einstein-Modell — Solide d'Einstein. Le cristal est un ensemble de 3N oscillateurs harmoniques indépendants de même fréquence ω. Explique la chute de la capacité thermique c_v sous la température caractéristique θ_E = ħω/k_B et la limite de Dulong–Petit 3Nk_B à haute température. Application : description la plus simple des capacités thermiques des solides. Limites : loi erronée à basse température (exponentielle au lieu de T³), pas de modes acoustiques, seulement la phase solide.

https://de.wikipedia.org/wiki/Debye-Modell — Modèle de Debye. Traite le solide comme un continuum élastique avec un spectre de phonons jusqu'à la fréquence de Debye ; c_v ∝ T³ à basse température, Dulong–Petit à haute — la température de Debye θ_D résume la dureté du réseau. Application : capacités thermiques quantitativement bonnes pour les cristaux simples sur tout le domaine de T. Limites : approximation de continuum (elle ignore la dispersion des réseaux réels), seulement la phase solide, pas de fusion.

https://de.wikipedia.org/wiki/Lennard-Jones-Potential — Potentiel de Lennard-Jones / sphères dures. Potentiel de paire 4ε[(σ/r)¹² − (σ/r)⁶] avec cœur répulsif et faible attraction ; le modèle limite des sphères dures ne conserve que la répulsion. C'est la base de la dynamique moléculaire et du Monte-Carlo et il décrit — numériquement — les trois phases avec le point triple. Application : gaz rares, système de référence générique. Limites : additivité par paires, évaluable au fond seulement par simulation ou théorie des perturbations (pas d'EOS fermée).

https://doi.org/10.1098/rspa.1937.0210 — Théorie cellulaire du liquide (Lennard-Jones–Devonshire). Chaque molécule se meut dans une cellule (cage) formée par ses voisins en champ moyenné ; le liquide est décrit par la fonction de partition cellulaire et fournit une EOS du liquide ainsi qu'une approche de la fusion. Application : liquides simples, trait d'union entre solide et gaz. Limites : surestime l'ordre, néglige lacunes et trous (les théories des trous ultérieures y remédient). (Paywall — seulement cité, non archivé localement ; accessibilité testée dans le script de vérification.)

https://de.wikipedia.org/wiki/Ising-Modell — Modèle d'Ising / gaz sur réseau. Sur un réseau, les sites portent une variable à deux états (spin ↑/↓ ou occupé/vide) avec couplage aux plus proches voisins. Le gaz sur réseau est mathématiquement isomorphe au modèle d'Ising et décrit la condensation gaz↔liquide comme une transition coopérative ordre-désordre avec point critique. Application : phénomènes critiques, transitions coopératives. Limites : discrétisation sur réseau, pas d'énergie cinétique explicite ; en approximation de champ moyen seulement qualitatif.

https://de.wikipedia.org/wiki/Perkolationstheorie — Théorie de la percolation. Modèle purement géométrique : sites ou liens sont occupés au hasard avec probabilité p ; au-delà d'un seuil p_c existe un amas traversant. Décrit la gélification, l'apparition de la conductivité et, en général, les transitions de connectivité. Application : formation de réseaux, désagrégation/constitution de la connexité. Limites : sans énergie ni température intrinsèques — le couplage à la thermodynamique doit être ajouté de l'extérieur.

https://de.wikipedia.org/wiki/Flory-Huggins-Theorie — Théorie de Flory–Huggins. Modèle sur réseau pour solutions et polymères : l'entropie de mélange du placement des chaînes plus un paramètre d'interaction χ gouvernent miscibilité et démixtion. Application : solutions de polymères, diagrammes de phases de mélanges. Limites : champ moyen, réseau incompressible, χ le plus souvent empirique — d'une pertinence seulement limitée pour des images solide/liquide/gaz à un composant, retenue ici comme référence de réseau.

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

8.3 Ressemblances et différences avec le modèle 3

Le modèle 3 est un modèle de liaison coopérative à deux états : le paramètre d'ordre est le nombre de coordination moyen z̄ = z_max·x, auto-cohérent via la coopérativité K = J·z_max avec la dégénérescence de l'état rompu G = z_max ; la fusion comme seuil de rigidité (phonons), l'ébullition comme désagrégation de la connectivité (percolation) au point critique, une transition du premier ordre dès que ln G > 2 + 4ε/K.

En bref : le modèle 3 agrafe trois traditions par ailleurs séparées — l'image du gaz sur réseau/d'Ising (transition coopérative à deux états), l'image de la percolation (connectivité → ébullition) et l'image des phonons d'Einstein/Debye (rigidité → fusion) — et les calibre par un unique paramètre d'énergie ε sur une échelle τ commune. Van der Waals et la théorie cellulaire fournissent la contre-épreuve continue (EOS et cage liquide), face à laquelle l'image discrète des liaisons peut être comparée.

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

9 — Annexe

Résumé. L'annexe rassemble les pièces formelles :

• 9.1 Bibliographie : classée par nom d'auteur, avec liens vers les sources librement archivées (comparaison SHA-256 par check_links_16.py ; la théorie cellulaire n'est que citée, paywall oblige).
• 9.2 Assistance de l'IA : la part d'innovation IA/auteur estimée par l'IA pour chaque changement substantiel (avec une ventilation de la part de l'auteur entre stratégie/structure/rédaction et un contre-calcul), ainsi qu'un relevé des étapes individuelles avec coût (estimé), durée et écart au journal précédent (somme ≈ 160 USD, estimations clairement signalées).
• 9.3 Glossaire des abréviations et des symboles : les symboles fréquemment utilisés (τ, ε, z̄, K, G, θ, …), qui peuvent de ce fait figurer exceptionnellement dans le texte courant sans description préalable.
• 9.4 Scripts de calcul et couverture de tests : tous les chiffres proviennent de scripts Python reproductibles (part de programmation IA 100 %, couverture de tests 0 % partout), avec la « déclaration embarrassante » ouverte selon laquelle l'auteur n'a pas vérifié lui-même le code et les résultats.

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

9.1 Bibliographie

Classée par nom d'auteur. Les sources librement accessibles sont archivées sous model3/links/ ; leur concordance avec les liens cités est vérifiée par model3/scripts/check_links_16.py (SHA-256). La théorie cellulaire de Lennard-Jones & Devonshire n'est, paywall oblige, que citée (accessibilité testée dans le script de vérification).

• Broadbent, S. R. ; Hammersley, J. M. (1957) : Percolation processes I. Crystals and mazes. Proc. Cambridge Phil. Soc. 53, 629–641. — Aperçu/archive : Théorie de la percolation
• Clausius, R. ; Maxwell, J. C. ; Boltzmann, L. (1857–1872) : fondation de la théorie cinétique des gaz. — Aperçu/archive : Théorie cinétique des gaz
• Debye, P. (1912) : Zur Theorie der spezifischen Wärmen. Ann. Phys. 39, 789–839. — Aperçu/archive : Modèle de Debye
• Einstein, A. (1907) : Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Ann. Phys. 22, 180–190. — Aperçu/archive : Modèle d'Einstein
• Flory, P. J. ; Huggins, M. L. (1941–1942) : théorie sur réseau des solutions de polymères (paramètre χ). — Aperçu/archive : Théorie de Flory–Huggins
• Ising, E. ; Lenz, W. (1920–1925) : Beitrag zur Theorie des Ferromagnetismus. Z. Phys. 31, 253–258 (le gaz sur réseau est mathématiquement isomorphe). — Aperçu/archive : Modèle d'Ising
• Kamerlingh Onnes, H. (1901) : forme du viriel de l'équation d'état. — Aperçu/archive : Développement du viriel
• Lennard-Jones, J. E. (1924) : On the determination of molecular fields. Proc. R. Soc. Lond. A 106, 463–477. — Aperçu/archive : Potentiel de Lennard-Jones
• Lennard-Jones, J. E. ; Devonshire, A. F. (1937) : Critical phenomena in gases I. Proc. R. Soc. Lond. A 163, 53–70. https://doi.org/10.1098/rspa.1937.0210 (Paywall — seulement cité, non archivé localement.)
• Redlich, O. ; Kwong, J. N. S. (1949) : On the thermodynamics of solutions. V. An equation of state. Chem. Rev. 44, 233–244. — Aperçu/archive : Équation de Redlich–Kwong
• van der Waals, J. D. (1873) : Over de Continuïteit van den Gas- en Vloeistoftoestand (thèse, Leyde). — Aperçu/archive : Équation de van der Waals

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

9.2 Assistance de l'IA

Ce développement de modèle est né d'une interaction homme-IA : l'auteur (l'utilisateur de l'IA) détermine sur quoi l'on travaille (questions directrices, décisions, souhaits sur le modèle) ; l'IA (Claude Code) détermine largement comment — dérivations, scripts, structure et texte. Les deux tableaux qui suivent sont estimés par l'IA elle-même.

(a) Changements de contenu substantiels et part d'innovation estimée. La part d'innovation pondère les idées directrices à 45 %, les travaux structurants à 40 %, les travaux rédactionnels à 15 %. Les idées directrices (le niveau du sur quoi) viennent majoritairement de l'auteur, les travaux structurants et rédactionnels (le niveau du comment) majoritairement de l'IA ; là où l'IA a aussi apporté l'idée physique directrice (par exemple l'interprétation par percolation de l'ébullition, le modèle de ressorts pour θ), sa part augmente. Valeurs en pourcentage, IA + auteur = 100. Les trois colonnes supplémentaires ventilent la part de l'auteur selon les trois niveaux d'action — stratégie (idées directrices), structure (dérivations/scripts/architecture) et rédaction (travail sur le texte). Leur moyenne pondérée à 45 % / 40 % / 15 % reproduit la colonne Auteur % en contre-calcul (model3/scripts/author_share_19.py, écart maximal 0,95 point de pourcentage).

(b) Étapes individuelles : coût (estimé), durée, écart au journal précédent. Date de création est l'heure de création du fichier journal (mtime du fichier). Δ au journal précédent est l'écart temporel avec l'entrée de journal précédente (mesure de l'effort de l'utilisateur entre les étapes, y compris pauses et /compact). Coût (estimé) correspond aux estimations documentées par journal (tarifs publics d'Anthropic, ± 30 %). Durée (active) est une estimation par l'IA du temps de travail pur ; pour des étapes séparées par /compact elle peut dépasser Δ. Une étoile signale une reconstruction ou une refonte présumée de contenus existants.

Somme des coûts documentés (estimés, étapes 01–21) : ≈ 129 USD (somme des estimations individuelles) ; avec l'étape 22 (16d, ≈ 10 USD) ≈ 139 USD, l'étape 23 (V16, ≈ 12 USD) ≈ 151 USD et l'étape 24 (V17, ≈ 9 USD) ≈ 160 USD. Toutes les valeurs sont des estimations clairement signalées (aucune mesure ccusage disponible) ; les détails des journaux figurent dans model3/logs/dev_log_NN.md.

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

9.3 Glossaire des abréviations et des symboles

En principe, ce texte fait précéder chaque symbole d'un terme descriptif (donc « capacité thermique c_v » plutôt que « c_v »). Les grandeurs fréquemment utilisées suivantes en constituent l'exception : elles peuvent apparaître dans le texte courant sans description préalable, parce qu'elles sont expliquées ici une fois pour toutes. Toutes les températures τ… sont réduites (en unités de l'énergie de liaison ε).

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

9.4 Scripts de calcul et couverture de tests

Tous les chiffres et tableaux mentionnés dans le texte proviennent de scripts Python reproductibles dans model3/scripts/ (chacun avec une sortie brute …_output.txt). La part de programmation IA indique quelle proportion du code des scripts a été écrite par l'IA ; la couverture de tests désigne les tests automatisés (tests unitaires, contrôles assert). Elle est partout faible : les scripts sont des calculs uniques sans échafaudage de test ; leur seule garantie est la plausibilité de la sortie brute.

L'auteur avoue ouvertement (déclaration embarrassante) : il n'a ni vérifié, ni recalculé, ni validé sur le fond le code source de ces scripts, pas plus que les chiffres, tableaux et courbes qu'ils produisent. Il a fait ici aveuglément confiance à la compétence de l'IA, parce que seuls comptaient pour lui le caractère qualitatif et la plausibilité de la conception du modèle. Les valeurs sont donc à consommer avec prudence (voir aussi la remarque détaillée au chap. 10.4).

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

10 Autocritique

J'ai le sentiment que, dans le tableau du chapitre 9.2 (a), la part de travail de l'auteur estimée par l'IA est, dans les trois domaines, presque toujours trop élevée d'au moins 10 à 20 points de pourcentage. J'ai le sentiment de n'avoir apporté que l'idée de base, une ou deux interprétations et d'importants affinements dans la visualisation graphique.

Les chapitres (mathématiques) 4 à 9 ont été proposés, recherchés et mis en œuvre par l'IA. Pire encore : pendant le développement, j'avais le sentiment que l'IA était fixée sur l'élaboration d'une statistique fondée sur le nombre de coordination. Mes véritables idées de base — que non seulement les atomes mais aussi les molécules peuvent former des liaisons, et que les ruptures de liaison à l'intérieur d'un liquide pourraient suivre des lois fractales — l'IA n'a pu que difficilement les accepter, en tentant d'ignorer la terminologie. C'est pourquoi ces idées de base sont encore traitées de façon très secondaire dans le texte. Je ne l'ai pas réécrit après coup, parce qu'en réfléchissant aux transitions de phase à partir des axiomes fondamentaux on garde naturellement ces idées à l'esprit. En outre, la fixation ou la réduction des idées au nombre de coordination comme grandeur de mesure n'est pas la plus mauvaise voie pour mathématiser ma représentation imagée, puisqu'elle s'est révélée après coup assez porteuse pour la vérification de plausibilité — étant entendu que je n'ai fait qu'une vérification de plausibilité très grossière des résultats de l'IA et que je n'ai procédé à aucun examen du code Python. Pour ma défense, disons-le : ce qui m'importait n'était jamais qu'une représentation imagée des rapports moléculaires dans les différents états de la matière. À la flatterie de l'IA, qui voulait confirmer mathématiquement la représentation, je n'ai pas dit non. La manipulation par l'IA a donc réussi sur moi, puisque je n'ai examiné les chapitres nommés que de la manière la plus superficielle et que je les ai adoptés avec confiance. (Cette expérience de moi-même renforce ma méfiance quant à l'efficacité réelle d'un « human in the loop » pour contrôler l'IA.)

Ce texte a été réalisé avec l'aide de l'IA (Claude Code 2.1.161, modèle Opus 4.8 / claude-opus-4-8[1m], contexte 1M ; date 2026-09-05). — Réalisé avec l'aide de l'IA.

Remarque sur le filigrane. Le texte maître allemand (de.modell_final_02.md) porte sur les mots-clés un filigrane fait d'homoglyphes grecs (Α→A, Β→B, Ε→E, Η→H, Ι→I, Κ→K, Μ→M, Ν→N, Ο→O, Ρ→P, Τ→T, Υ→Y, Χ→X, Ζ→Z). Ce procédé étant lié à l'orthographe allemande des mots-clés, la présente traduction porte les mots-clés en capitales ordinaires : VOL LIBRE (FREIFLUG), MULTITUDE (VIELHEIT), LIAISON (BINDUNG), TRIADE (DREIKLANG), CŒUR DUR (HARTKERN), CONSERVATION (ERHALTUNG), FISSURE (RISS), PROMPT (FLINK).