Termodinámica estadística de la sustancia pura ideal y sus estados de agregación

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

Resumen general

El modelo 3 es un modelo cualitativo, de mecánica estadística, de una sustancia ideal que deriva los tres estados de agregación —sólido, líquido y gaseoso— así como sus transiciones a partir de una sola imagen de enlace, la de los enlaces moleculares. El principio básico reside en tener en cuenta tres depósitos de energía para la energía interna: la energía de vibración en los enlaces moleculares, la energía de traslación de las moléculas y la energía cinética cooperativa en la rotación de micelas como forma idealizada de los flujos internos de materia que compensan las fluctuaciones internas de densidad. Una magnitud fundamental importante en este contexto es el número de coordinación medio z̄ = z_max·x, que se establece de forma autoconsistente mediante un enlace cooperativo de dos estados (tipo Ising, cooperatividad K = J·z_max, degeneración del estado roto G = z_max). De él se siguen tres transiciones sobre una escala τ común, calibrada mediante un único parámetro de energía ε: la fusión como umbral de fonones/rigidez (z̄ alcanza la rigidez), la ebullición como colapso de la conectividad (percolación) en el punto crítico y una transición opcional sólido(1)→sólido(2) (FCC↔BCC) con el escalón de c_v sin parámetros de 1,0 k_B. Una transición pasa a ser de primer orden en cuanto ln G > 2 + 4ε/K. Aquí la degeneración del estado roto G es el número de contactos vecinos que pueden romperse (G = z_max, el número de coordinación máximo), la energía de enlace ε es la energía ganada por contacto vecino y la cooperatividad K = J·z_max es el producto de la fuerza de acoplamiento J, de tipo Ising, por z_max; la transición se vuelve, pues, de primer orden tanto antes cuanto mayor es la multiplicidad de contactos G y cuanto más fuerte es la cooperatividad K frente a la energía de enlace ε.

El artículo conduce desde los siete axiomas (capítulo 1) a través del análisis de contradicciones (capítulo 2) y de una descripción figurada con tiras de cómic estáticas (capítulo 3) hasta el armazón de mecánica estadística —enlace cooperativo, funciones de partición y capacidades caloríficas (capítulos 4–5)—, a la plausibilidad matemática en 3D (capítulo 6) y finalmente a la calibración y la ecuación de estado (capítulo 7). Las ideas se complementan con el barrido 3D completo de c_v de la fusión de primer orden (Δz̄, calor latente ΔU, desplazamiento del pico), un modelo de muelles de Einstein para la rigidez fonónica θ(red) y la línea de transformación sólido(1)→sólido(2) vía Clausius–Clapeyron. θ(red) se refina después en una densidad de estados de fonones de tipo Debye; la transformación sólido(1)→sólido(2) se introduce además como línea en el plano (τ, p) (calibrada con hierro) y se añade una rama Redlich–Kwong separada junto a la opción por defecto de van der Waals. θ(red) se comprueba asimismo sobre una red de Lennard-Jones completa (FCC estable, BCC ideal inestable a cizalla) y la topología sólido(1)→sólido(2) se comprueba con titanio (segunda confirmación) y argón (control negativo). El capítulo 8 sitúa el modelo entre once concepciones establecidas de la termodinámica estática. El apéndice documenta las fuentes bibliográficas (ordenadas por autor) y el apoyo de la IA (proporción de innovación estimada, coste estimado y duración por paso de desarrollo); las ayudas reúnen el modo de lectura, los puntos críticos restantes y las perspectivas.

Enunciado central. Una única medida cooperativa del enlace, z̄, engrapa tres tradiciones por lo demás separadas —la imagen del gas en red / de Ising (transición cooperativa de dos estados), la imagen de la percolación (conectividad → ebullición) y la imagen de los fonones de Einstein/Debye (rigidez → fusión)— en un cuadro de fases coherente y cualitativamente completo, deliberadamente dentro de una imagen ideal de tipo van der Waals, pudiendo por supuesto derivarse el edificio teórico de una imagen dinámica fácilmente imaginable. No es la exactitud sino el carácter figurado lo que constituye la verdadera fuerza del modelo y lo que caracteriza su valor científico como modelo.

Índice

• Resumen general
• 1 — Los axiomas con sus palabras clave (consolidados)
• 2 — Análisis de contradicciones y lagunas
• 3 — Descripciones figuradas
• 3.1 Las tres fases aisladas
• 3.2 Transiciones de fase
• 3.3 Por qué la velocidad del sonido en los líquidos baja con T
• 3.4 Cómo surge el movimiento browniano — con una noción afinada de la micela
• 3.5 Por qué las reacciones en las capas límite de las micelas se aceleran
• 3.6 La ventana líquida tridimensional
• 4 — Modelo de mecánica estadística
• 4.1 Grados de libertad y variables
• 4.2 Enlace cooperativo de dos estados (tipo Ising, parámetro J)
• 4.3 Balance energético — con término cooperativo
• 5 — Funciones de partición y capacidades caloríficas
• 5.1 Función de partición factorizada
• 5.2 Capacidad calorífica por fase (analítica)
• 5.3 Curva numérica c_v(τ) **con** término cooperativo — calor latente cuantificado
• 5.4 Umbral de percolación de conectividad como resultado (confirmado sin cambios) - 5.4.a Dos umbrales distintos claramente separados
• 5.5 Orden de la transición de fusión
• 5.6 Estudio de tamaño finito del umbral de percolación
• 5.7 Barrido 3D completo de c_v(τ) de la fusión de primer orden
• 6 — Plausibilidad matemática de la interpretación 3D
• 6.1 Cambio del grado de libertad de rotación 2D → 3D
• 6.2 Influencia de las distintas redes cristalinas — la conectividad como rasgo distintivo
• 6.3 La transición adicional sólido(1) → sólido(2)
• 6.4 Consideraciones: la anomalía acústica del agua líquida
• 6.5 Alternativa: campos de vórtices fractales en lugar de la rotación de micelas
• 6.6 Rigidez fonónica microscópica θ(red): un modelo de muelles de Einstein
• 6.7 Estimación matemática de la metaestabilidad (2D triangulada ↔ cuadrada, 16b)
• 6.8 Refinamiento de tipo Debye de θ(red)
• 6.9 θ(red) a partir de una red de Lennard-Jones completa
• 7 — Calibración y ecuación de estado
• 7.1 Método de calibración (un parámetro de energía ε) — ajustado
• 7.2 Tabla de calibración
• 7.3 Eje de densidad/volumen: ecuación de estado de van der Waals con núcleo duro
• 7.4 Calibración cuantitativa l/g: a → ε, b → D²
• 7.5 Superficie (τ, v) común con punto triple
• 7.6 Calor latente frente a entalpías de fusión
• 7.7 Unificación de las escalas de energía: ε ↔ ε_lg
• 7.8 Punto triple sin parámetros a partir de z̄_c
• 7.9 Cálculo de control 3D: una ventana líquida limpia sobre una sola escala ε
• 7.10 Barrido real de percolación de enlaces en 3D (FCC/BCC) y sensibilidad del borde de la ventana
• 7.11 Barrido real en τ de la cinética de enlace en 3D: τ_b(3D) y la línea de fusión sobre una escala ε
• 7.12 La línea de transformación sólido(1)→sólido(2) en la superficie (τ, v)
• 7.13 Línea de transformación sólido(1)→sólido(2) en el diagrama (τ, p), calibrada con hierro
• 7.14 Redlich–Kwong como rama separada de la imagen ideal
• 7.15 sólido(1)→sólido(2) en otros polimorfos: titanio y control con gas noble
• 8 — Comparación bibliográfica: concepciones de la termodinámica estática
• 8.1 Tabla de conjunto
• 8.2 Los modelos uno a uno
• 8.3 Semejanzas y diferencias con el modelo 3
• 9 — Apéndice
• 9.1 Bibliografía
• 9.2 Apoyo de la IA
• 9.3 Glosario de abreviaturas y símbolos
• 9.4 Scripts de cálculo y cobertura de pruebas

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

1 — Los axiomas con sus palabras clave (consolidados)

Resumen. Siete axiomas dotados de palabras clave fijan la sustancia ideal: número de partículas conservado, movimiento puramente clásico, un enlace cooperativo de dos estados y el número de coordinación medio z̄ como parámetro de orden. Sostienen todas las derivaciones que siguen; en el texto maestro alemán una marca de agua de homoglifos señala las palabras clave (clave al final del texto).

Los axiomas precisados:

1. VUELO LIBRE (axioma adicional). Entre dos moléculas que pasan una junto a otra sin estar enlazadas no actúan ni fricción ni fuerzas electrostáticas. El momento y la energía solo cambian por la repulsión dura del NÚCLEO DURO en el choque y por la formación y ruptura de un ENLACE. El axioma solo importa desde el punto de vista de la mecánica clásica y del balance energético: garantiza que no existen canales disipativos a nivel molecular — la energía de movimiento nunca se «pierde», sino que solo pasa a los tres depósitos de la TRÍADA o a energía de enlace. (La fricción y la viscosidad macroscópicas surgen de forma emergente del rompecabezas deslizante y de la GRIETA, no como fuerza fundamental microscópica.)
2. MULTITUD. Una sustancia ideal consta de muchísimas moléculas iguales; las propiedades macroscópicas se describen estadísticamente en el límite termodinámico (N → ∞).
3. ENLACE. Dos moléculas geométricamente vecinas cualesquiera pueden portar un enlace. Un enlace intacto rebaja la energía en la cantidad finita ε; su ruptura cuesta exactamente ε. El enlace es un sistema de dos estados {intacto, roto}.
4. TRÍADA. La energía interna se reparte clásicamente en tres depósitos: (a) traslación de las moléculas, (b) vibración en enlaces intactos (un muelle armónico por enlace), (c) rotación de bloques moleculares rígidos (grado de libertad rotacional colectivo de las regiones conexas).
5. NÚCLEO DURO. Si la distancia molecular cae por debajo de un valor mínimo σ, actúa una repulsión dura (volumen propio). Junto con el ENLACE resulta un potencial de pares efectivo (coraza dura más pozo de enlace finito).
6. CONSERVACIÓN. El momento total y la energía total se conservan. Una molécula aislada es energéticamente neutra — no almacena energía en sí misma; la energía vive en los enlaces (ENLACE) y en los grados de libertad cinéticos (TRÍADA).
7. GRIETA. Los enlaces se deshacen y se forman según el principio de agrietamiento/cicatrización: cuantos más enlaces vecinos faltan → mayor probabilidad de ruptura; cuantos más vecinos enlazados → mayor probabilidad de cicatrización. La cinética es cooperativa.
8. ÁGIL. Ruptura y cicatrización transcurren rápidamente y con poca activación; el subsistema de enlaces se halla a toda temperatura en equilibrio local y puede separarse adiabáticamente del movimiento molecular (fundamento de la función de partición factorizada).

Actitud de modelado. La coordinación geométrica se conserva de forma idealizada; solo cambia el número de coordinación z por las rupturas. El movimiento surge como en un rompecabezas deslizante: los huecos (procedentes del volumen propio del NÚCLEO DURO y de las rupturas de la GRIETA) migran, los bloques avanzan. La imagen deliberadamente «falsa» de cristalitos que giran y migran se conserva como base intuitiva.

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

2 — Análisis de contradicciones y lagunas

Resumen. El análisis examina sistemáticamente el modelo en cuanto a contradicciones internas y lagunas y continúa el estado del trabajo; las debilidades reconocidas se nombran y se resuelven en los capítulos siguientes o se marcan como límite honesto del modelo.

Lista de puntos críticos

1. Doble cómputo de la traslación (axioma 3a frente a 3c). El texto en bruto nombra la traslación dos veces (traslación molecular y traslación de bloque). Dos depósitos de traslación independientes no existirían. Resolución: la traslación de bloque es solo un modo colectivo de las traslaciones moleculares ya contadas; genuinamente nuevo es únicamente el grado de libertad de rotación del bloque. En el balance de grados de libertad (cap. 4) la traslación se cuenta exactamente una vez.
2. «Energéticamente neutra» (axioma 5) frente a «los enlaces absorben energía» (axioma 2/3b). Contradicción aparente. Resolución: la neutralidad afecta a la molécula como cuerpo; la energía reside entre las moléculas (en el ENLACE) y en modos cinéticos colectivos, nunca en un amortiguador intramolecular. Molécula ≠ enlace — con ello la contradicción queda resuelta.
3. Falta de atracción (el axioma 4 solo nombra la repulsión). Sin atracción no habría fase condensada. Resolución: la atracción la aporta el ENLACE (pozo de profundidad ε), la repulsión el NÚCLEO DURO (coraza dura en σ). Juntos, un potencial de pares efectivo.
4. Retroalimentación sin término opuesto (axioma 6). «Más vecinos ausentes → más rupturas» es por sí solo una retroalimentación positiva (avance de la grieta hasta la desintegración total). Resolución: el término de cicatrización (más vecinos → más cicatrización) más el factor de Boltzmann de la energía de enlace finita ε (ENLACE) equilibran la grieta y producen una transición definida en lugar de un desbocamiento.
5. Separación de escalas temporales no explicitada (axioma 7). «Rápida y de poca activación» solo resulta utilizable cuando de ello se hace la separación adiabática entre dinámica de enlaces y dinámica de traslación. Resolución: consignada explícitamente como requisito de la factorización Z = Z_trans·Z_vib·Z_rot·Z_bond.
6. Falta de geometría y falta de mecanismo de movimiento. Los axiomas no nombran ni la red de coordinación ni lo que hace posible el flujo. Resolución (como suposición adicional explícitamente señalada): una geometría de coordinación de referencia fija con ocupación variable de los enlaces; como caso 2D canónico la coordinación triangular (z_max = 6). Movilidad a través del mecanismo de huecos del rompecabezas deslizante. Estos dos puntos son las únicas adiciones genuinas más allá de los siete axiomas y están señalados en el capítulo 6 como cuestión estratégica.
7. Laguna latente de disipación. La versión 07 dejaba abierto si las moléculas no enlazadas pierden energía al pasar de largo. Resuelta mediante el VUELO LIBRE: no hay fricción ni electrostática microscópicas; todos los flujos de energía y momento discurren exclusivamente por los choques del NÚCLEO DURO y la cinética del ENLACE. Con ello el balance energético (cap. 4.3) queda cerrado y la estadística de Maxwell–Boltzmann es exacta en el límite gaseoso (gas ideal libre entre choques elásticos).
8. Incoherencia energética del parámetro de orden. La versión 08 usaba el término cooperativo −(K/2)x² en la energía libre g(x) (para generar x(τ)), pero no en la energía interna U/N. Eso es termodinámicamente incoherente: el mismo término de campo molecular debe aparecer en ambas. Resolución (cap. 4.3): −(K/2)x² se traslada a U/N; el calor latente aumenta con ello de forma medible (cap. 5.3), sin alterar las mesetas sólida y gaseosa.
9. Calibración l/g solo cualitativa. La versión 09 tenía el eje de densidades (ecuación de estado de van der Waals) solo en forma reducida. Resolución (cap. 7.4): la línea de ebullición absoluta se calcula vía el punto crítico (T_c, p_c) y T_b(p_atm) se contrasta con valores medidos; se nombra la desviación sistemática (regla de Guldberg).
10. Ejes separados sin punto triple. Los ejes de fusión y de ebullición eran cálculos separados. Resolución (cap. 7.5): ambos se reúnen en un plano (τ, v), se construye el punto triple y se marca el punto crítico.

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

3 — Descripciones figuradas

Resumen. El capítulo 3 explica de forma intuitiva lo sólido, lo líquido y lo gaseoso así como sus transiciones (con tiras de cómic estáticas): las tres fases aisladas (3.1), las transiciones de fase recristalización/fusión/ebullición/sublimación junto con el punto triple y el crítico (3.2), la anomalía acústica (3.3), el movimiento browniano y la noción afinada de la micela (3.4–3.5), así como la ventana líquida tridimensional (3.6). A partir de estas descripciones se elaboraron esquemas (tiras de cómic) que deben ayudar a representarse ante el ojo interior las relaciones moleculares dinámicas en las distintas fases. Esta capacidad de representación es el verdadero núcleo del modelo. Adquiere valor cualitativo y científico cuando la mecánica estadística permite mostrar de forma plausible que el modelo hace predecibles tendencias y observaciones empíricas cualitativas incluso sin grandes matemáticas.

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

3.1 Las tres fases aisladas

Sólido. Casi todos los vecinos están enlazados entre sí mediante enlaces moleculares (ENLACE); la media de los números de coordinación z de las moléculas está cerca del máximo geométrico. La red molecular es un único bloque rígido sin huecos (NÚCLEO DURO, empaquetamiento denso) — un rompecabezas deslizante interno como en el líquido no es posible en la versión compacta. Único movimiento: el temblor en torno a las posiciones de reposo (vibración de enlace, TRÍADA-b), determinado por los enlaces moleculares con los vecinos. Momento y energía oscilan elásticamente en la red (CONSERVACIÓN, VUELO LIBRE: sin pérdidas). En el caso límite (temperatura muy baja, empaquetamiento ideal) se tiene el estado ideal de un cristal perfecto.

Líquido. La GRIETA cooperativa ha atravesado la red: faltan enlaces a manchas, la media de los números de coordinación z de las moléculas ha descendido, el bloque rígido se ha fragmentado interiormente en muchos bloques menores. Del monolito ha surgido un montón de cascajo. Estos bloques trasladan y rotan (TRÍADA-c). Entre ellos se han abierto huecos por los que se desplazan unos junto a otros — el rompecabezas deslizante del movimiento browniano está activo. La geometría de coordinación local se conserva de forma idealizada; solo ha cambiado cuántos enlaces existen realmente y se forman.

Gas. La GRIETA ha vencido: los enlaces están rotos casi en todas partes (z → 0). En el caso extremo ninguna molécula tiene ya enlace molecular con otra. Los bloques se han desintegrado en moléculas sueltas. Vuelan libres, acopladas solo por choques del NÚCLEO DURO (con VUELO LIBRE en medio); las velocidades siguen a Maxwell–Boltzmann (CONSERVACIÓN). Vibración de enlace y rotación de bloque se han extinguido — traslación pura (TRÍADA-a). En el caso límite (temperatura alta, presión baja) se tiene el estado ideal de un gas ideal.

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

3.2 Transiciones de fase

Recristalización (sólido(1)/sólido(2)). Aún por debajo de la fusión, el sólido puede cambiar de red sin perder la rigidez: dos empaquetamientos ideales con números de coordinación diferentes (FCC z_max = 12 ↔ BCC z_max = 8) tienen entalpías libres μ(τ, p) distintas; en μ_FCC = μ_BCC la geometría de coordinación salta bruscamente (no solo la ocupación z̄ como en la fusión). Se hace visible un pequeño escalón de c_v (Δc_v = variación de z_max/4 = 1,0 k_B) más un débil calor latente de transformación. Afirmación de la IA: solo es concebible en 3D (conectividad discreta, cuantizada; cap. 6.3). En el caso bidimensional, la recristalización genuina, es decir, la transición entre dos fases sólidas, resulta más fácil de imaginar. En el plano solo hay tres teselaciones regulares sin huecos. Se puede cubrir por completo una superficie (teselarla) con triángulos equiláteros, con cuadrados o con hexágonos regulares. Si en estas redes las aristas representan los enlaces moleculares y los vértices las moléculas, entonces es concebible una transición de fase sólido(triangular) → sólido(cuadrado). La molécula del teselado triangular tendría todavía seis enlaces con moléculas vecinas, mientras que la del retículo cuadrado solo tendría cuatro.

Fusión (s/l). Al subir la temperatura la red se desgarra en suficientes puntos como para que el bloque rígido total pierda su rigidez (la GRIETA cooperativa cae por debajo del umbral de rigidez de Maxwell z̄ = 2d — la red se vuelve blanda). Se forman micelas termodinámicas, es decir, regiones que ya no están unidas por enlaces a las regiones vecinas mayores. Desde el punto de vista molecular, el monolito se convierte en un montón de cascajo y adquiere una capacidad de flujo limitada. Los primeros huecos permiten el desplazamiento de bloques moleculares, como en el rompecabezas deslizante. Los enlaces ausentes hacia la vecindad permiten la rotación interna en las micelas. La rotación de micelas (TRÍADA-c) se pone en marcha. La ruptura de muchos enlaces (ENLACE) devora calor latente.

Ebullición (l/g). Un calentamiento ulterior hace percolar los huecos mismos: se unen en un espacio vacío continuo en el que los bloques se disuelven y moléculas sueltas salen volando. La ebullición es el paso de «huecos en la materia» a «materia en el hueco».

Sublimación (s/g). A baja presión falta la zona líquida intermedia: la red cede moléculas sueltas directamente al espacio vacío, sin que se produzca la formación de bloques como en las micelas. Cabe preguntarse aquí qué influencia pueden tener sobre la formación o no formación de micelas termodinámicas los aspectos geométricos (moléculas largas, moléculas que se entrelazan), los aspectos electrostáticos (alineación por momentos magnéticos o dipolares) y los distintos tipos de enlace molecular (van der Waals, interacciones entre sistemas mesómeros, puentes de hidrógeno).

Punto triple. Justamente allí coinciden dos umbrales: la rigidez de Maxwell (sólido ↔ líquido, pérdida de rigidez) y la percolación de huecos (líquido ↔ gas) ocurren en el mismo (τ, p). Las tres imágenes (bloque rígido [monolito], bloques que se deslizan [montón de cascajo], moléculas libres [nube de polvo]) coexisten.

Punto crítico. Al subir la presión el salto líquido→gas se hace cada vez menor; la fracción de hueco cambia de forma cada vez más continua. En el punto crítico la diferencia desaparece por completo — «bloques densos con poco vacío» y «mucho vacío con pocas moléculas» pasan continuamente el uno al otro. («Desde el punto de vista del líquido el vaso está medio vacío; desde el del gas, medio lleno.»)

Recristalización en la imagen 2D. El cómic muestra ahora la recristalización como caso 2D, uno debajo del otro: abajo el retículo denso triangulado (coordinación z = 6, el ideal 2D de menor energía), arriba el retículo cuadrado suelto (z = 4). La transición triangulado → cuadrado aumenta el volumen (trabajo contra la presión) y reduce el número de enlaces. Se ve que el número de coordinación cambia mientras la rigidez de un retículo cristalino permanece.

Predicción cualitativa de transiciones de energía y capacidades caloríficas. Solo a partir de la conectividad z se pueden estimar ambas cifras características: (i) la energía de enlace por partícula es U/N ≈ −(ε/2)·z (cada enlace cuenta la mitad), es decir, tanto menor cuanto mayor sea z; (ii) el salto de c_v en el cambio de red sigue la regla del modelo Δc_v = variación de z_max/4. Para el caso 2D triangulado (z = 6) ↔ cuadrado (z = 4) resulta una diferencia de energía de ΔU ≈ ε por partícula y un escalón de c_v de Δc_v = (6 − 4)/4 = 0,5 k_B — la mitad que el valor 3D (1,0 k_B). Regla mnemotécnica: más denso = menor energía, pero más caro en la transición (mayor calor latente, mayor escalón de c_v), porque al expandirse se liberan más enlaces y más grados de libertad congelados.

¿Determina la estructura de la red la metaestabilidad? En parte. La estructura de la red fija el orden de energías — qué forma es el estado fundamental (triangulado, μ más bajo) y cuál la excitada (cuadrada). Que la forma suelta sobreviva como metaestable a baja temperatura no lo decide ella sola: en el modelo estrictamente ideal (sin barrera) la sustancia adopta en cada (τ, p) inmediatamente la fase de μ más bajo — allí no existe metaestabilidad. Esta surge solo por una energía de activación finita (barrera de nucleación) que frena cinéticamente la transformación. La geometría de la red aporta, pues, la jerarquía termodinámica; la barrera, la persistencia cinética — el mismo ingrediente que explica también la fase sobresaturada (véase abajo; cuantitativamente cap. 6.7).

Fase sobresaturada — energía de activación. La figura adicional (3.2, fila 2) muestra una fase gaseosa sobresaturada sobre un líquido: junto a moléculas sueltas, pequeños grupos de 2–3 partículas agregadas (barras grises) que vuelan como un todo. El modelo explica la sobresaturación por una energía de activación apreciable para los enlaces moleculares: los cúmulos pequeños (2–3) tienen una relación borde-superficie desfavorable, sus enlaces de borde son incompletos — solo un núcleo crítico crece espontáneamente. Mientras no surja un núcleo supercrítico, el gas permanece sobresaturado, aunque μ_gas > μ_líquido ya permitiría la condensación. Precisamente esa barrera falta en la imagen ideal de equilibrio y se nombra aquí como contribución autónoma situada «fuera de la imagen ideal» (estimación cuantitativa en cap. 6.7).

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

3.3 Por qué la velocidad del sonido en los líquidos baja con T

El sonido es una onda elástica: v = √(K_eff/ρ). K_eff surge exclusivamente de los enlaces intactos (ENLACE) — ellos son los muelles. Al subir T se rompen cooperativamente cada vez más enlaces (GRIETA); la red de muelles se adelgaza y K_eff cae más abruptamente que la densidad ρ. La ruptura domina el numerador, no el denominador — por eso baja v. Interpretación molecular: la transmisión del sonido por enlaces moleculares es más rápida que la transmisión mecánica por choques. Cuantas más grietas aparecen en un líquido al aumentar la temperatura, tanto más a menudo debe transportarse el sonido por choques mecánicos en lugar de por enlaces moleculares. En consecuencia, la velocidad del sonido se hace más lenta al subir la temperatura. (La anomalía del agua — un máximo acústico hacia 74 °C — se trata en el cap. 6.4.)

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

3.4 Cómo surge el movimiento browniano — con una noción afinada de la micela

En sentido termodinámico-estadístico, una micela no es aquí un agregado de tensioactivos sino un bloque rodeado por la GRIETA: los enlaces ausentes (los no-enlaces en la frontera del bloque) forman la «envoltura» en torno a un cristalito encerrado (el núcleo intactamente enlazado). Ese núcleo puede ejecutar un movimiento de rotación cooperativo interno (TRÍADA-c) — el bloque gira como un todo mientras su envoltura frontera lo desacopla de los vecinos.

El movimiento browniano surge entonces así: mediante el intercambio de enlaces moleculares (GRIETA, ÁGIL) moléculas sueltas pueden abandonar la envoltura de una micela y pasar a otra micela. Cada uno de estos cambios transmite momento (CONSERVACIÓN; sin pérdidas gracias al VUELO LIBRE) y altera bruscamente la forma y el momento angular de ambas micelas implicadas. Una gran partícula de prueba (grano de polvo) es bañada por muchas de estas reorganizaciones de micelas no correlacionadas; la suma de los impulsos es nula en promedio, pero tiene varianza finita → camino aleatorio. Cuanto mayor es T, tanto más frecuentes son los sucesos de GRIETA y los cambios de micela, y tanto más fuerte el temblor. El movimiento browniano es así la sombra macroscópica, o consecuencia, del proceso microscópico de intercambio entre micelas.

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

3.5 Por qué las reacciones en las capas límite de las micelas se aceleran

Advertencia preliminar embarazosa. La argumentación que sigue es una deducción a partir de la imagen y, hasta donde sé, no ha sido demostrada ni estudiada empíricamente — si bien he de confesar que no me tomé la molestia de investigarlo a fondo. El único indicio de novedad es que los sistemas de IA no mencionan esta conclusión. Esta tesis es un ejemplo de cómo las representaciones figuradas pueden llevar a nuevas preguntas.

En disolventes de disolución moderada, los reactivos se adhieren preferentemente a las envolturas frontera de las micelas — en el modelo, precisamente las fronteras grieta/bloque (GRIETA). Dos efectos aceleran la reacción: (i) la reducción dimensional — los reactivos ya no se buscan en todo el volumen sino a lo largo de una interfase casi bidimensional; una búsqueda 2D acierta estadísticamente antes que una 3D. (ii) La concentración — la envoltura frontera reúne a los reactivos, su densidad local aumenta. En imagen, las fronteras de grieta actúan como raíles guía que conducen a los compañeros el uno hacia el otro.

El efecto dimensional tiene un núcleo riguroso: el paseo aleatorio (Pólya). Que una búsqueda bidimensional acierte estadísticamente mejor que una tridimensional no se sigue solo de la intuición, sino de un teorema del cálculo de probabilidades. George Pólya (1887–1985) demostró en 1921 para el paseo aleatorio sobre una red: en una y en dos dimensiones el caminante vuelve con probabilidad 1 (es decir, con seguridad) infinitas veces a cada punto, tanto en la red plana ℤ² como en la recta ℤ. A partir de tres dimensiones eso ya no vale: la probabilidad de retorno es estrictamente menor que 1 (en la red cúbica simple ℤ³ solo alrededor de 0,34). Gráficamente, en la imagen de Pólya: «Un hombre borracho seguro que acaba encontrando el camino a casa; un pez borracho, en cambio, no.» Trasladado a dos reactivos que deben encontrarse varias veces para un choque exitoso, esto significa: en un disolvente muy bueno con paseo aleatorio tridimensional libre, los compañeros una vez separados se alejan para siempre con probabilidad finita — los intentos repetidos de reacción fracasan y el rendimiento queda bajo por esta razón puramente geométrico-estadística. En una capa límite de micela (casi) bidimensional, en cambio, el reencuentro es seguro (probabilidad 1): los compañeros se vuelven a encontrar una y otra vez hasta que la reacción tiene éxito. La capa límite no solo concentra, sino que a la vez pone a disposición la dimensión más favorable para los reencuentros. La afirmación (deliberadamente puntiaguda) reza: los rendimientos en disolventes muy buenos tienden a ser peores por razones de paseo aleatorio.

Referencia: G. Pólya (1921), Über eine Aufgabe der Wahrscheinlichkeitsrechnung betreffend die Irrfahrt im Straßennetz, Math. Ann. 84, 149–160 — citado aquí según M. Eisermann, Höhere Mathematik 3 (Univ. de Stuttgart), ejercicio W-2×2. (Original consultado el 5.9.2026).

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

3.6 La ventana líquida tridimensional

A diferencia de las figuras de 3.1–3.5, esta no muestra ninguna escena molecular sino el resultado del cálculo de los cap. 7.9/7.10 como esquema sobre un eje τ (en unidades de la energía de enlace ε). El objetivo era hacer visible por qué el retículo 2D aislado apenas admite un dominio para una fase líquida, mientras que con el retículo 3D ello es posible sin dificultad, rigiendo en ambos casos la misma regla de rigidez (z̄_c = 2d).

Imagen ① — 2D (sin ventana). Sobre el eje τ, τ_m = 1,20 y τ_triple = 1,16 se superponen prácticamente (τ_triple por debajo de τ_m), con τ_c = 1,50 justo encima: el dominio de la fase líquida queda comprimido en una rendija — el sólido pasa casi directamente a gas (cerca de la sublimación).

Imagen ② — 3D (ventana limpia). El mismo eje, pero dilatado: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00. Entre τ_triple y τ_c está el campo líquido autónomo (ancho hacia abajo, τ_triple − τ_m = 1,63; estrecho hacia arriba, τ_c − τ_triple = 0,18). El delgado borde superior está marcado como tal.

Imagen ③ — la causa. z̄_c = 2d crece con la dimensión (4 → 6), con lo que v_solid baja (0,50 → 0,667 en volúmenes reducidos, es decir, el sólido se vuelve «más suelto» respecto al líquido), y τ_c = z_max/4 se duplica (1,50 → 3,00). El punto crítico se aleja muy por encima del punto de fusión. Solo así surge espacio para la fase líquida. Es un puro efecto dimensional.

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

4 — Modelo de mecánica estadística (sin cambios desde la versión 09)

Resumen

El capítulo 4 establece el armazón de mecánica estadística del modelo en 2D (retículo triangular, z_max = 6):

• 4.1 Grados de libertad y variables: la magnitud de orden es la fracción de enlaces x ∈ [0,1] con el número de coordinación medio z̄ = z_max·x; por bloque ≥ 2 se añade un grado de libertad de rotación, y la temperatura τ se cuenta en unidades de la energía de enlace ε.
• 4.2 Enlace cooperativo de dos estados: la ecuación de enlace autoconsistente x = 1/(1 + G·exp[−(ε + K·x)/τ]) (K = J·z_max, G = z_max) se sigue como condición de estacionariedad de la energía libre de Bragg–Williams; su mínimo global proporciona la rama física, y el cap. 5.5 decide si esa rama es continua o presenta un salto (primer orden).
• 4.3 Balance energético con término cooperativo: U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x² (2D) reúne traslación, vibración de enlace y rotación de bloque, la energía de configuración de enlaces y la nueva energía cooperativa de campo molecular −(K/2)x². La capacidad calorífica c_v asociada recibe el término adicional positivo −K·x·ẋ, que agudiza el pico de fusión pero deja intactos la meseta sólida y la rama gaseosa.

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

4.1 Grados de libertad y variables

• N moléculas sobre una geometría de coordinación de referencia fija. Caso 2D canónico: retículo triangular, z_max = 6, un GDL de rotación por bloque (MULTITUD, ENLACE). Este capítulo vale explícitamente para la interpretación 2D; la generalización 3D está en el cap. 6.
• Variable de enlace s_ij ∈ {0, 1} por par geométrico de vecinos (1 = intacto).
• Número de coordinación z_i = Σ_j s_ij; media z̄ = 2 N_bond / N — el parámetro de orden. Se tiene z̄ = z_max · x con la fracción de enlaces x = p_b ∈ [0, 1].
• Unidad de energía τ ≡ k_B T (clásicamente, ½τ por grado de libertad cuadrático). Todas las energías en unidades de la energía de enlace ε.

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

4.2 Enlace cooperativo de dos estados (tipo Ising, parámetro J)

Un enlace aislado (energía −ε intacto, 0 roto) tiene, en el equilibrio de Maxwell–Boltzmann, p_b⁰(τ) = 1/(1 + e^(−ε/τ)). La GRIETA cooperativa (axioma 6) desplaza la estabilidad efectiva del enlace en un campo molecular proporcional a la fracción de enlaces x. En aproximación de campo molecular (Ising), la ecuación de enlace autoconsistente resulta:

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

Aquí:

• K = J·z_max es la fuerza de acoplamiento cooperativa (retroalimentación de la GRIETA: más enlaces vecinos → mayor estabilidad efectiva del enlace). J > 0 agudiza la transición; por encima de un J crítico z̄(τ) se vuelve discontinuo — ese es el calor latente de la fusión. El caso puro de dos estados es el límite J = 0 (sin campo molecular, crossover continuo).
• G > 1 es la degeneración del estado roto: una molécula liberada gana espacio de fases (libertad de traslación). G impulsa la ruptura de enlaces a τ alta y hace que x → 1/(1+G) en lugar de quedarse en ½ (caso límite del gas). Canónicamente G = z_max.

Esta ecuación se sigue como condición de estacionariedad de una energía libre de Bragg–Williams por molécula (n_slots = z_max/2 plazas de enlace):

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

El término −(K/2)x² es el acoplamiento cooperativo (ferromagnético), el último término la degeneración de ruptura. El mínimo global de g sobre x proporciona la rama física; donde dos mínimos se hacen igual de profundos se sitúa la transición de primer orden. ÁGIL (axioma 7) justifica que x(τ) tome a cada temperatura el valor de equilibrio.

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

4.3 Balance energético — con término cooperativo

El término de campo molecular −(K/2)x² es una energía genuina (independiente de τ) y debe aparecer en U/N. El balance coherente reza así (2D):

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

• τ·[…]: traslación (1) + vibración de enlace (z̄/4) + rotación de bloque (½ n_blk) — TRÍADA.
• −(ε/2)·z̄: energía de configuración de enlaces — ENLACE.
• −(K/2)·x²: energía cooperativa de campo molecular — GRIETA. Nuevo en U/N.

La capacidad calorífica asociada c_v = d(U/N)/dτ (en unidades k_B) recibe una contribución adicional:

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

donde el nuevo término −K·x·ẋ es positivo (ẋ < 0) y actúa exactamente allí donde muchos enlaces se rompen a la vez. Agudiza el pico de fusión sin alterar la meseta sólida (ẋ ≈ 0) ni la rama gaseosa.

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

5 — Funciones de partición y capacidades caloríficas

Resumen

El capítulo 5 calcula funciones de partición y capacidades caloríficas y proporciona los datos del modelo para la calibración (cap. 7):

• 5.1 Función de partición factorizada: por la separación de escalas temporales, Z = Z_trans·Z_vib·Z_rot·Z_bond se factoriza.
• 5.2 Capacidad calorífica por fase: c_v = d(U/N)/dτ; el término positivo −(ε/2)(dz̄/dτ) genera el pico de tipo Schottky como análogo del calor latente.
• 5.3 Curva numérica c_v(τ): pico de c_v 5,51 en τ_m = 1,20, meseta sólida 2,50, rama gaseosa → 1; el calor latente L ≈ 1,41 ε por molécula (con término cooperativo) queda triplicado frente a 0,48 ε (sin él).
• 5.4 / 5.4a Dos umbrales separados: la percolación de conectividad z̄_c = z_max·p_c ≈ 2,08 (en τ_b ≈ 1,85) marca la desintegración de los cúmulos (cerca de la ebullición) y es un resultado; hay que distinguirla claramente de la rigidez genérica de Maxwell z̄_c = 2d (densidad de fusión, pérdida de rigidez).
• 5.5 Orden de la transición de fusión: criterio ln G > 2 + 4ε/K; el 2D (G = 6) es continuo para todo J, el 3D (G = 12) es de primer orden a partir de K > ~8,2 (demo: Δz̄ ≈ 5,3, ΔU ≈ 1,75 ε).
• 5.6 Estudio de tamaño finito: p*(L) converge monótonamente desde abajo hacia p_c = 0,3473; z̄_c ≈ 2,08 es robusto frente al tamaño finito.
• 5.7 Barrido 3D completo de c_v: a partir de K > K_crit ≈ 8,25 (J > 0,68) se abre el salto genuino Δz̄ con calor latente ΔU; el pico de fusión se desplaza con J hacia τ mayores (desplazamiento positivo de +0,5 … +2,2 respecto al pico 2D continuo).

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

5.1 Función de partición factorizada

Por la separación de escalas temporales (ÁGIL) la función de partición canónica se factoriza:

Z = Z_trans · Z_vib · Z_rot · Z_bond.

• Z_trans = (1/N!)·(A/λ_T²)^N, λ_T = h/√(2π m τ) (traslación 2D clásica ideal; MULTITUD, CONSERVACIÓN, VUELO LIBRE).
• Z_vib = (τ/ħω₀)^(N_bond) (oscilador armónico clásico por enlace intacto; TRÍADA-b).
• Z_rot = ∏_(bloques≥2) (2π I_b τ)^(1/2)/ħ (rotor rígido 2D clásico por bloque; TRÍADA-c).
• Z_bond = Σ_{s} exp[(ε/τ)·Σ s_ij + (K/τ)·término cooperativo] (suma de enlaces de tipo Ising; ENLACE, GRIETA). La evaluación en campo molecular proporciona x(τ), z̄(τ) según 4.2.

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

5.2 Capacidad calorífica por fase (analítica)

A partir de U/N (cap. 4.3), c_v = d(U/N)/dτ (en unidades k_B):

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

El término −(ε/2)·(dz̄/dτ) es positivo (dz̄/dτ < 0: calentar rompe enlaces) y domina allí donde muchos enlaces se rompen a la vez — el pico de tipo Schottky, análogo del calor latente.

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

5.3 Curva numérica c_v(τ) con término cooperativo — calor latente cuantificado

Enlace autoconsistente (mínimo global), n_blk(τ) de un Monte Carlo de percolación de enlaces (union-find, L = 100, 5 realizaciones), c_v = dU/dτ (diferencias centradas), U/N incluido −(K/2)x² (cv_scan_09.py, ε = 1, J = 0,55, G = z_max = 6). Puntos de apoyo seleccionados:

Resultado (numérico). El orden de las tres fases sigue confirmado y ahora es más nítido:

c_v(sólido) ≈ 2,50 < c_v(líquido, máximo) ≈ 5,51 > c_v(rama gaseosa) ↓ → 1.

Calor latente (área del pico sobre la línea de base). Integrando (c_v − línea de base) sobre la ventana de fusión τ ∈ [0,70 ; 1,70] (línea de base = recta por c_v(0,70) y c_v(1,70)) resulta

L ≈ 1,41 ε por molécula (con término cooperativo), L ≈ 0,48 ε por molécula (sin término cooperativo, balance antiguo), ΔL ≈ +0,92 ε por molécula.

El término cooperativo eleva el pico en un factor ≈ 1,7 (3,25 → 5,51) y triplica el calor latente. Físicamente: la ruptura simultánea de muchos enlaces vecinos (GRIETA) almacena no solo energía de configuración (−ε/2·z̄) sino también energía de campo molecular (−K/2·x²) — ambas se liberan al fundir.

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

5.4 Umbral de percolación de conectividad como resultado

La fracción del cúmulo gigante P_∞(τ) colapsa en τ ≈ 1,85; allí x ≈ 0,347 = p_c(triángulo) = 2·sin(π/18) = 0,3473, luego z̄_c = z_max·p_c ≈ 2,08 — un resultado, no algo impuesto. Este es el umbral de conectividad (cúmulo que atraviesa) y marca la desintegración de los cúmulos cerca de la ebullición (τ_b); no debe confundirse con la rigidez de Maxwell (z̄_c = 2d, fusión) — cf. la separación limpia en el cap. 5.4a.

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

5.4.a Dos umbrales distintos claramente separados

La doble denominación hasta ahora borrosa «percolación/rigidez» queda ahora resuelta. Se trata de dos umbrales físicamente distintos que no deben confundirse:

Justificación. La fusión es la pérdida de la rigidez mecánica de un bloque denso, no la desintegración de la conexidad. La magnitud relevante es por tanto el recuento de restricciones de Maxwell (campo medio / rigidez genérica): una red D-dimensional de fuerzas centrales se vuelve isostática exactamente cuando la coordinación media alcanza z̄ = 2D (2D → 4, 3D → 6). La conectividad de enlaces (z̄_c = 2,08) queda muy por debajo; una red apenas conexa en z̄ = 2,08 hace tiempo que es blanda (modos suaves), es decir, ya líquida. Por eso la conectividad marca la desintegración de cúmulos posterior, cerca de la ebullición, mientras que la rigidez fija la solidificación/licuación anterior (densidad de fusión). Esta separación fija el punto triple del cap. 7.8 sin parámetros y de forma realista; no cambia el umbral de conectividad 2,08 confirmado en 5.4/5.6, sino solo su lectura (τ_b = cerca de la ebullición, no la fusión).

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

5.5 Orden de la transición de fusión

Criterio analítico. La región espinodal (dos mínimos de g) existe para τ < K/4; la transición se sitúa aproximadamente en τ_t ≈ (ε + K/2)/ln G. Un salto genuino (transición de primer orden) exige τ_t < K/4, es decir, en el límite K → ∞:

ln G > 2 + 4ε/K.

De ahí se sigue un resultado notable sobre la degeneración de ruptura G = z_max:

Confirmación numérica (máx. |Δx| por paso en τ, detect_first_order):

Lectura. En el modelo 2D canónico (G = 6) la fusión es por principio continua (una transición «blanda» / crossover) — para ningún J surge un salto genuino, porque ln 6 < 2. Solo la mayor coordinación/degeneración del empaquetamiento 3D (G = 12) permite, a partir de una cooperatividad suficiente (K > 8,2), una transición de primer orden con calor latente finito. Eso concuerda con la realidad: la fusión 2D es a menudo continua (tipo KTHNY), la fusión 3D es de primer orden. J = 0 (sin campo molecular) es siempre el caso de referencia continuo. — Físicamente J designa la relación entre energía cooperativa y energía de enlace (K/ε = J·z_max/ε): cuánto estabiliza un enlace vecino intacto al propio.

Demostración de primer orden (motivada por 3D, cv_scan_09.py, J = 1,0, G = 12, z_max = 12): z̄ salta en Δz̄ ≈ 5,30 en τ ≈ 2,82; el calor latente como salto de energía ΔU ≈ 1,75 ε por molécula. El mismo K en un retículo G = 6 (2D) sigue siendo continuo — el carácter del orden depende únicamente de G = z_max.

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

5.6 Estudio de tamaño finito del umbral de percolación

Percolación de enlaces en retículos triangulares L × L (periódicos), fracción de cúmulo gigante P_∞(p) escalada sobre L (cv_scan_09.py; L = 50/100/200/400 con 8/5/3/2 realizaciones). Umbral exacto p_c = 2·sin(π/18) = 0,34730 → z̄_c = z_max·p_c = 2,0838.

P_∞(p) (media ± dispersión):

Estimación del umbral p*(L) (P_∞ cruza 0,5):

Resultado. p*(L) converge monótonamente desde abajo hacia p_c = 0,3473; la dispersión en p_c cae de ±0,103 (L = 50) a ±0,017 (L = 400) — la transición se agudiza con L, pero no desplaza el umbral (valor de percolación universal). Con ello z̄_c ≈ 2,08 = z_max·p_c queda confirmado como valor del modelo robusto frente al tamaño finito.

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

5.7 Barrido 3D completo de c_v(τ) de la fusión de primer orden

La única serie de demostración de primer orden del cap. 5.5 (J = 1,0, G = 12) se amplía ahora a un barrido completo de c_v(τ) sobre la cooperatividad J (firstorder_3d_16.py, conectividad FCC para n_blk, balance energético 3D U/N = τ·[3/2 + z̄/4 + (3/2)n_blk] − (ε/2)z̄ − (K/2)x²). El criterio analítico de orden del cap. 5.5 (primer orden ⇔ K > K_crit = 4ε/(ln 12 − 2) = 8,25) separa los casos; el pico de referencia 2D continuo está en τ_2D = 1,20 (c_v ≈ 5,5).

*caso continuo: Δz̄/ΔU son el crossover abrupto (no un salto genuino); τ_transición es allí el pico de c_v.

Resultado. (i) El caso 3D continuo J = 0,55 da τ_transición = τ_m ≈ 1,72 — la comprobación cruzada con la cinética de enlace independiente (bondkinetics_3d_15.py: τ_m(FCC) = 1,70) concuerda. (ii) A partir de K > K_crit ≈ 8,25 (J > 0,68) se abre el salto genuino Δz̄ con calor latente finito ΔU (salto de energía a través de la discontinuidad) — la fusión pasa a ser de primer orden, exactamente como exige el criterio ln 12 > 2. (iii) La comprobación cruzada J = 1,0 reproduce el valor de demostración Δz̄ ≈ 5,4 del cap. 5.5. (iv) La altura del salto Δz̄ y la posición τ_transición crecen con J: la transición se desplaza hacia arriba — la mayor coordinación del empaquetamiento 3D y la cooperatividad más fuerte elevan la temperatura de fusión. Frente al pico 2D continuo (τ_2D = 1,20) se obtiene, pues, en todos los casos un desplazamiento positivo del pico (+0,5 … +2,2) que crece monótonamente con K. El carácter de primer orden de la fusión 3D queda así acreditado no solo en un punto, sino en todo el rango de cooperatividad.

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

6 — Plausibilidad matemática de la interpretación 3D

Resumen

El capítulo 6 hace plausible la generalización de la imagen de enlace de 2D a 3D de forma puramente matemática (sin visualización HTML propia):

• 6.1 Grado de libertad de rotación 2D → 3D: la rotación de bloque pasa de un eje a tres, la traslación de dos direcciones a tres, z_max de 6 a 12; los valores base de c_v suben (c_v,sólido 2,5 → 4,5, c_v,gas 1 → 3/2), la forma de la curva c_v(τ) permanece.
• 6.2 La conectividad como rasgo distintivo: z_max es una magnitud cuantizada y separa bruscamente el empaquetamiento más denso (FCC/HCP, z_max = 12) del BCC (z_max = 8), sin estado intermedio.
• 6.3 Transición adicional sólido(1)→sólido(2): de la discretitud se sigue el polimorfismo FCC↔BCC con un escalón de c_v sin parámetros de 1,0 k_B; en la imagen ideal estricta (θ_BCC = θ_FCC) no aparece, y solo un BCC más blando (θ_BCC < θ_FCC, por debajo del umbral de apertura 0,61) abre τ_ss < τ_m.
• 6.4 / 6.5 Consideraciones: la anomalía acústica del agua líquida como cambio difuminado de orden local entre dos tipos de red, y una alternativa de campos de vórtices fractales a la rotación de micelas (puramente cualitativa).
• 6.6 / 6.8 / 6.9 θ(red) microscópica en tres niveles: modelo de muelles de Einstein (θ_BCC/θ_FCC ≈ 0,82…0,94), densidad de estados de fonones de tipo Debye (≈ 0,87…1,00) y red de Lennard-Jones completa (≈ 1,0; BCC ideal inestable a cizalla, FCC estado fundamental estable) — todos por encima del umbral 0,61, de modo que el polimorfismo real queda fuera de la imagen ideal armónica.
• 6.7 Metaestabilidad (2D): cálculo de orden de magnitud para la transformación triangulada ↔ cuadrada (ΔU = +1 ε, ΔA ≈ +15,5 %) junto con la barrera de nucleación ΔG\* = πλ²/Δμ_A como aporte cinético.

Este capítulo hace plausible la generalización 3D de forma puramente matemática. La 3D nunca se visualiza en una versión HTML. El objetivo es únicamente mostrar cómo cambian el grado de libertad de rotación y la conectividad de 2D a 3D y qué nuevas transiciones se vuelven así concebibles.

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

6.1 Cambio del grado de libertad de rotación 2D → 3D

El único punto cualitativamente nuevo frente al 2D es la rotación de los bloques:

• 2D: un bloque rígido tiene 1 GDL de rotación (giro en el plano) → U_rot/N = ½ n_blk τ.
• 3D: un bloque rígido tiene 3 GDL de rotación (tres ejes de Euler) → U_rot/N = (3/2) n_blk τ.

Al mismo tiempo la traslación pasa de 2 a 3 GDL (U_trans/N = (3/2)τ) y el número de coordinación máximo de z_max = 6 (triángulo) a z_max = 12 (empaquetamiento compacto de esferas). El balance energético en 3D reza entonces

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

y los valores base de la capacidad calorífica se desplazan coherentemente:

La estructura de la curva (meseta sólida, máximo líquido por −(ε/2)dz̄/dτ, caída hacia el valor del gas) permanece idéntica; solo los niveles escalan con la dimensión y z_max.

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

6.2 Influencia de las distintas redes cristalinas — la conectividad como rasgo distintivo

En 3D hay que distinguir al menos dos redes ideales, que se separan por su conectividad (número de coordinación z_max):

La conectividad z_max es el rasgo distintivo decisivo: fija tanto el nivel de capacidad calorífica del sólido (z_max/4) como el umbral de percolación z̄_c. Puesto que cada red admite solo un número de coordinación fijo, la conectividad es una magnitud cuantizada y finita — no hay transición continua entre z_max = 12 y z_max = 8, sino solo el salto discreto de un tipo de red al otro.

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

6.3 La transición adicional sólido(1) → sólido(2)

De la finitud/discretitud de la conectividad se sigue una transición ausente en el modelo 2D: sólido(1) → sólido(2) (polimorfismo/alotropía). Dos fases sólidas con z_max distinto (por ejemplo FCC ↔ BCC) tienen entalpías libres μ(τ, p) distintas; en la combinación temperatura/presión μ_{sólido(1)} = μ_{sólido(2)} la sustancia cambia de red sin fundir. En la imagen del modelo la geometría de coordinación cambia bruscamente (no solo la ocupación z̄ como en la fusión) — la única situación en que la geometría de referencia, por lo demás fija, salta ella misma. Sería observable un pequeño escalón de c_v (variación de z_max/4) más un calor latente de transformación, claramente más débil que el pico de fusión.

El criterio de primer orden ln G > 2 + 4ε/K (cap. 5.5) proporciona ahora una razón concreta de por qué justamente la interpretación 3D es cualitativamente más rica: solo allí (G = z_max = 12, ln 12 = 2,485 > 2) puede la fusión ser de primer orden — y solo allí existe el cambio discreto de red sólido(1) → sólido(2) (FCC ↔ BCC), porque la conectividad es cuantizada/finita. Ambos fenómenos son la misma causa con dos ropajes: la coordinación mayor y discreta del empaquetamiento 3D genera tanto el salto de fusión nítido como el polimorfismo. El retículo 2D aislado (G = 6) no tiene ninguno de los dos — coherente con el cap. 5.5.

Cuantificación (bondkinetics_3d_15.py). La transición se calcula ahora de forma esbelta. Ambas redes reciben una entalpía libre en la imagen armónica clásica (las partes comunes p·v y de traslación se cancelan en la comparación):

μ(τ) = h₀ − c_v·τ·ln(τ/θ), h₀ = −(ε/2)·z_max, θ = ħω/k_B (rigidez fonónica en unidades ε),

con los valores característicos sin parámetros de los cap. 6.1/6.2:

De ahí el escalón de c_v sin parámetros: |Δc_v| = variación de z_max/4 = (12−8)/4 = 1,0 k_B (en la transición FCC→BCC el sólido se vuelve 1 k_B más «ligero» en capacidad calorífica) — mucho menor que el pico de fusión (varios k_B, cap. 5.3/7.11).

La temperatura de transformación τ_ss (μ_FCC = μ_BCC) y el calor latente de transformación L = τ_ss·(S_BCC − S_FCC):

Resultado (hallazgo honesto). En la imagen ideal estricta (la misma rigidez fonónica para todos los enlaces, θ_BCC = θ_FCC) no hay transición sólido(1)→sólido(2): el FCC gana tanto la entalpía (más baja, −6 frente a −4 ε) como — por el mayor número de enlaces (c_v = 4,5 > 3,5) — la entropía vibracional; μ_FCC permanece a toda temperatura por debajo de μ_BCC. Solo una rigidez fonónica dependiente de la red θ_BCC < θ_FCC (el BCC más abierto vibra más blandamente y por tanto tiene más entropía) lleva μ_BCC por debajo de μ_FCC a alta temperatura — entonces se abre la transformación, y esta queda por debajo del pico de fusión τ_m, de modo que surge el orden sólido(1)→sólido(2)→líquido→gas (τ_ss < τ_m < τ_c ≤ τ_b(3D), cap. 7.11). La rigidez blanda del BCC es una extensión moderada fuera de la imagen ideal estricta (conocida en la realidad como el polimorfismo γ/δ del hierro: BCC a alta temperatura, FCC en medio); el escalón de c_v de 1,0 k_B y el pequeño calor de transformación positivo son, en cambio, internos al modelo y sin parámetros.

¿Cuán pequeño es realmente θ_BCC/θ_FCC? (theta_springs_16.py). Para liberar la transformación del parámetro libre θ, la rigidez fonónica se estima ahora microscópicamente (cap. 6.6): un modelo de muelles de Einstein con muelles de par centrales da θ_BCC/θ_FCC ≈ 0,82 (solo vecinos más próximos) hasta ≈ 0,94 (con la segunda capa del BCC). El umbral de apertura exigido por la comparación de μ está, en cambio, en θ_BCC/θ_FCC < 0,61 (por debajo de él existe τ_ss y queda bajo τ_m). La razón geométrica de rigidez es, pues, demasiado grande: en la imagen ideal estricta el FCC permanece estable a toda temperatura, no hay transición sólido(1)→sólido(2) puramente geométrica. El polimorfismo real (hierro γ↔δ) requiere contribuciones adicionales (electrónicas/magnéticas/anarmónicas) fuera de la imagen armónica ideal — precisamente la «extensión moderada» supuesta antes. Eso confirma el hallazgo honesto y convierte la razón θ de una suposición libre en un límite del modelo estimado y nombrado.

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

6.4 Consideraciones: la anomalía acústica del agua líquida

Los líquidos normales muestran v(T) decreciente de forma monótona (cap. 3.3). El agua tiene un máximo acústico hacia 74 °C: v sube primero con T antes de bajar. En el modelo esto puede leerse como una competencia entre dos tipos de red locales (cap. 6.2): un orden local abierto, de baja coordinación, semejante al hielo (z_max pequeño, muchos «agujeros»), y un orden local denso, de mayor coordinación. A T baja domina la estructura abierta; calentar rompe primero el orden local abierto en favor del más denso — la conectividad efectiva y con ella K_eff suben (v sube). Solo por encima del máximo predomina la ruptura de enlaces habitual de la GRIETA, K_eff cae, v baja. La anomalía es así expresión de un cambio de orden local de tipo sólido(1)↔sólido(2), difuminado e incompleto, dentro de la fase líquida — concebible solo en el modelo 3D con dos tipos de red, no en el retículo 2D aislado.

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

6.5 Alternativa: campos de vórtices fractales en lugar de la rotación de micelas

El movimiento browniano local (cap. 3.4) se explicó mediante la rotación y el intercambio de micelas. Alternativamente puede representarse como efecto de campos de vórtices fractales: las fronteras de bloque, remodeladas constantemente por la GRIETA, generan un campo de corrientes de momento locales que — de modo semejante a los flujos turbulentos — presenta vórtices autosemejantes en todas las escalas de longitud. Una partícula de prueba es entonces impulsada no por choques discretos de micelas, sino por un continuo de vórtices invariante de escala. Ambas imágenes (micelas discretas ↔ vórtices fractales) son casos límite complementarios de la misma cinética de GRIETA; cuál de ellas rige depende de la relación entre el tamaño de la partícula de prueba y el tamaño medio de las micelas. Esta breve discusión es puramente cualitativa y no se desarrolla aquí matemáticamente.

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

6.6 Rigidez fonónica microscópica θ(red): un modelo de muelles de Einstein

La comparación sólido(1)→sólido(2) del cap. 6.3 depende de la rigidez fonónica relativa θ_BCC/θ_FCC. En lugar de elegirla libremente, se estima aquí a partir de la coordinación y la geometría (theta_springs_16.py).

Modelo. Un átomo central se halla en el pozo armónico de sus z vecinos; cada enlace es un muelle de par radial κ. Un desplazamiento u se proyecta isótropamente con ⟨cos²⟩ = 1/d = 1/3 (d = 3) sobre las direcciones de enlace, de modo que la constante elástica efectiva es k_eff = (z/3)·κ. La frecuencia de Einstein es ω = √(k_eff/m), luego θ = ħω/k_B ∝ √k_eff. Con igual masa e igual muelle de par se sigue

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

Refinamientos. (B) La segunda capa del BCC (6 átomos en r₂/r₁ = 2/√3 = 1,155) con un muelle escalado con la distancia κ(r) ∝ (r₁/r₂)ⁿ (n = 6…8, tipo Lennard-Jones) eleva la razón a 0,91…0,94 (el BCC se vuelve efectivamente más rígido por la cercana segunda capa). (C) El enlace más corto del BCC a igual densidad (r_nn(BCC)/r_nn(FCC) = 0,972, −2,8 %) endurece ligeramente el muelle y da 0,89…0,91. En conjunto una ventana geométrica θ_BCC/θ_FCC ≈ 0,82…0,94.

Umbral de apertura. La comparación de μ (cap. 6.3) abre la transformación (τ_ss existe y queda bajo τ_m = 1,70) solo para θ_BCC/θ_FCC < 0,61. La ventana geométrica queda claramente por encima.

Hallazgo. El modelo de muelles puro no hace el BCC lo bastante blando: en la imagen ideal estricta el FCC permanece estable y no surge una transición sólido(1)→sólido(2) de origen geométrico. Con ello θ_BCC/θ_FCC ha pasado de suposición libre a magnitud estimada cuyo valor natural (≈ 0,8–0,9) mantiene cerrada la transformación — coherente con el hallazgo honesto del cap. 6.3. El polimorfismo real queda, pues, demostrablemente fuera de la imagen armónica ideal (aportes electrónicos/magnéticos, anarmonicidad), y no en la mera geometría de la red.

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

6.7 Estimación matemática de la metaestabilidad (2D triangulada ↔ cuadrada)

Para la discusión cualitativa del cap. 3.2, he aquí el cálculo de orden de magnitud correspondiente en la imagen 2D (discos duros, diámetro D, distancia de contacto a = D).

(1) Energía. Por partícula U/N = −(ε/2)·z (cada enlace la mitad): triangulada z = 6 → −3 ε, cuadrada z = 4 → −2 ε. La transformación triangulada → cuadrada cuesta, pues, ΔU = +1 ε por partícula (endotérmica; la forma suelta es más rica en energía).

(2) Volumen (área) y trabajo. Área por partícula: triangulada A_tri = (√3/2)·D² ≈ 0,866 D², cuadrada A_quad = D². La expansión es ΔA = A_quad − A_tri ≈ 0,134 D² (+15,5 %) — el aumento de volumen mostrado en el cómic. Contra una presión (reducida) p hay que aportar el trabajo W = p·ΔA; la entalpía libre sube en Δμ = ΔU + p·ΔA − τ·ΔS. El término entrópico ΔS > 0 (la forma suelta tiene modos más blandos, un escalón de c_v mayor Δc_v = (6−4)/4 = 0,5 k_B) baja μ_cuad por debajo de μ_tri solo a τ alta — por eso la triangulada es el estado fundamental de baja temperatura y la cuadrada la forma de alta temperatura / metaestable.

(3) Barrera de nucleación (2D). La metaestabilidad presupone una barrera. En 2D un núcleo de la nueva fase (radio r) crece con una ganancia de área −Δμ_A·π r² (Δμ_A = fuerza impulsora por área) frente a la aportación de borde +λ·2π r (λ = energía de arista por longitud): ΔG(r) = −Δμ_A·π r² + 2π λ r. Máximo en r\* = λ / Δμ_A con la barrera ΔG\* = π λ² / Δμ_A. El núcleo crítico r\ se encoge al crecer la fuerza impulsora Δμ_A; mientras la escala térmica τ = k_BT no alcance la barrera (**ΔG\ ≳ unos pocos k_BT), la fase suelta (o, en la ebullición, la sobresaturada) permanece congelada cinéticamente. Para cúmulos pequeños (2–3 partículas, como en la tira sobre la sobresaturación) la relación borde-superficie es máximamente desfavorable — quedan a la izquierda de r\* y se disuelven en promedio. Con ello la energía de activación mostrada en el cómic es la misma magnitud ΔG\* que sostiene tanto la forma cristalina metaestable como la fase gaseosa sobresaturada; es un aporte cinético** al modelo de equilibrio ideal (sin barrera).

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

6.8 Refinamiento de tipo Debye de θ(red)

El modelo de muelles de Einstein (cap. 6.6) sustituye toda la distribución de fonones por una sola frecuencia. El enfoque presente refina θ(red) mediante una densidad de estados de fonones esbelta por red (theta_debye_17.py): los mismos muelles de par centrales κ(r) = r^(−n) sobre la primera y la segunda capa de vecinos forman la matriz dinámica D(k) = Σ_j κ_j·(R̂_j⊗R̂_j)·(1 − cos(k·R_j)), que se diagonaliza sobre una malla k de la zona de Brillouin primitiva (Monkhorst–Pack, 12³ puntos) mediante un solucionador analítico 3×3, solo biblioteca estándar. De las tres ramas de fonones se siguen dos medidas de rigidez de tipo Debye como momentos de frecuencia:

• 2.º momento θ(2) ∝ √⟨ω²⟩ — el límite de Einstein (el promedio en la ZB de 1 − cos da 1, luego ⟨ω²⟩ = z·κ/3); control, solo vecinos más próximos, geometría pura: √⟨ω²⟩_BCC / √⟨ω²⟩_FCC = 0,816 = √(8/12) — exactamente el valor de Einstein del cap. 6.6.
• momento logarítmico θ(0) ∝ exp⟨ln ω⟩ — la magnitud propiamente de tipo Debye, que pondera las ramas acústicas blandas que el valor de Einstein ignora.

Todas las constantes con escala (ħ/k_B, la escala de κ, el prefactor de Debye) se cancelan en la razón θ_BCC/θ_FCC; solo permanece la forma de la DOS. Con la segunda capa de vecinos (que es lo que estabiliza mecánicamente el BCC) y una ley de muelle n = 0…8 resulta

θ_BCC/θ_FCC ≈ 0,87…1,00 (banda de error sobre la definición del momento y la ley del muelle).

Eso es más bien mayor que el valor puro de Einstein (0,82…0,94), porque el enlace más corto del BCC a igual densidad endurece los muelles. El umbral de apertura 0,61 exigido por la comparación de μ (cap. 6.3) sigue claramente por debajo: tampoco la DOS de fonones hace el BCC lo bastante blando. El hallazgo del cap. 6.6 queda así confirmado y precisado por la aproximación más fina — en la imagen ideal estricta el FCC permanece estable a toda temperatura; no surge ninguna transición sólido(1)→sólido(2) de origen puramente geométrico, y el polimorfismo real (hierro, cap. 7.13) necesita aportes fuera de la imagen armónica ideal.

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

6.9 θ(red) a partir de una red de Lennard-Jones completa

El cap. 6.8 usa un muelle central truncado κ(r) = r^(−n) sobre nn + 2nn. La versión 17 examina la misma cuestión con una auténtica constante de fuerza procedente de un potencial de pares sobre la red completa (theta_ljlattice_18.py): cada enlace se convierte en un muelle derivado del potencial de Lennard-Jones φ(r) = 4ε[(σ/r)¹² − (σ/r)⁶], con la parte longitudinal κ_L = φ″(r) y la parte transversal κ_T = φ′(r)/r. La matriz dinámica

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

suma sobre todas las capas de vecinos hasta un radio de corte (2,0…3,0 nn); cada red se relaja previamente a su propio equilibrio LJ (constante de red libre, densidad igual no impuesta). Ese es el test honesto de si la aproximación esbelta nn+2nn (cap. 6.8) acierta el orden de las rigideces.

Resultado (determinista, solo biblioteca estándar):

Dos hallazgos, ambos reforzando los cap. 6.6/6.8:

1. El FCC es el estado fundamental LJ estable (E ≈ −8,4 ε por átomo con corte 3,0 nn, frente a −8,0 ε para el BCC) — la conocida primacía del empaquetamiento denso en un potencial de pares.
2. El BCC LJ ideal es inestable a cizalla: alrededor del 3…3,5 % de los modos de fonones son imaginarios (ω² < 0, inestabilidad de Born). La imagen ideal del potencial de pares puro no porta, pues, ningún BCC mecánicamente estable; sobre los modos reales restantes θ_BCC/θ_FCC ≈ 1,0 — el umbral de apertura 0,61 exigido por la comparación de μ con mayor razón no se alcanza.

Con ello la afirmación de los cap. 6.6/6.8 es robusta frente al refinamiento del muelle truncado a la red completa y del modelo puramente longitudinal a la matriz dinámica completa: en la imagen ideal clásica no surge ninguna transición sólido(1)→sólido(2) puramente geométrica — el polimorfismo real (hierro cap. 7.13, titanio cap. 7.15) necesita aportes fuera del potencial de pares (electrónicos/magnéticos, anarmónicos). Un auténtico cálculo de fonones ab initio se deja deliberadamente fuera (abandonaría el carácter esbelto de la imagen ideal).

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

7 — Calibración y ecuación de estado

Resumen. El capítulo 7 calibra el modelo sobre sustancias reales mediante el único parámetro de energía ε y desarrolla la rama líquido/gas hasta una ecuación de estado:

• 7.1 / 7.2 Método y tabla de calibración: acoplamiento lineal de temperatura T = (ε/k_B)·τ con un punto de anclaje en el punto de fusión; de ahí la razón universal refinada R = τ_b/τ_m ≈ 1,54.
• 7.3 / 7.4 Ecuación de estado de van der Waals con núcleo duro: la EOS reducida p_r = 8τ_r/(3v_r − 1) − 3/v_r² con punto crítico (1, 1, 1) y coexistencia de Maxwell; a ∝ ε, b ∝ D². Los puntos de ebullición quedan sistemáticamente en T_b/T_c ≈ 0,48 (distancia de Guldberg); un cálculo comparativo de Redlich–Kwong los eleva a ≈ 0,60 — como presupuesto de error nombrado, sin cambio de EOS.
• 7.5 Superficie (τ, v) común: la cúpula de coexistencia (ebullición), la línea de fusión y la pared dura encajan con un punto triple (τ_r,triple ≈ 0,47 para v_solid = 0,40).
• 7.6 Calor latente: L ≈ 1,41 ε da ΔS_fus ≈ 1,175 R (cerca de la regla de Richards); los metales se aciertan excelentemente (cobre ×1,00), mientras que los cristales iónicos, con puentes de hidrógeno y moleculares quedan fuera de la imagen ideal.
• 7.7 Unificación ε ↔ ε_lg: k_B T_c = z_max·ε/4 acopla ambas escalas (ε_lg = (27/32)·z_max·ε); la predicción 3D τ_c/τ_m = 2,50 acierta con las sustancias moleculares simples, mientras que el 2D sigue siendo caso especial.
• 7.8 / 7.9 Punto triple sin parámetros: v_solid se sigue del umbral de rigidez z̄_c = 2d, no de la conectividad; en 3D se abre una ventana líquida limpia (τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00), en 2D permanece cerrada (papel especial del 2D, ya triplemente documentado).
• 7.10 / 7.11 Barridos 3D reales: la percolación de enlaces en FCC/BCC da z̄_c(3D) ≈ 1,5–1,6; el barrido en τ de la cinética de enlace muestra τ_m ≈ τ_rig (fusión = rigidez, bajo) y τ_b(3D) ≈ τ_c (ebullición = desintegración de la conectividad en el punto crítico).
• 7.12 / 7.13 / 7.15 Línea sólido(1)→sólido(2): del punto τ_ss la presión hace una línea (dp/dτ = ΔS/Δv > 0, Δv = +8,9 %), incorporada como cuarta línea frontera en la superficie (τ, v) y calibrada como línea (τ, p) con hierro (γ→δ), confirmada con titanio (α→β) y contrastada con argón como control negativo.
• 7.14 Rama Redlich–Kwong: como variante autónoma de la imagen ideal con calibración propia junto a la opción por defecto de van der Waals (T_b/T_c ≈ 0,60, error −1 %).

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

7.1 Método de calibración (un parámetro de energía ε) — ajustado

Lineal: T = (ε/k_B)·τ, un punto de anclaje por sustancia. Anclaje en el punto de fusión: ε/k_B = T_m,exp/τ_m, predicción T_b,pred = T_m,exp·(τ_b/τ_m). Con el pico de fusión refinado (término cooperativo) el máximo de c_v está ahora en τ_m = 1,20 (en vez de 1,15), y la percolación de conectividad (cerca de la ebullición) sigue sin cambios en τ_b = 1,85. De ahí la razón universal refinada:

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

La tabla de calibración de la versión 08 (cap. 7.2) escala linealmente con R: todos los T_b,pred bajan en el factor 1,54/1,61 ≈ 0,957 (≈ −4,3 %). Eso mejora ligeramente la sobreestimación de los gases nobles (argón +54 % → +48 %) y empeora un poco los aciertos de NaCl y hierro (NaCl −0,6 % → −4,8 %). El enunciado básico permanece: R vale para sustancias de empaquetamiento denso isótropo y enlace no direccional; gases nobles y cadenas quedan fuera de la imagen ideal.

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

7.2 Tabla de calibración

Tomada sin cambios de modell_08.md cap. 7.2 (con la nota sobre el escalado lineal en R del 7.1). Reproducible vía cv_scan_08.py (tabla) y cv_scan_09.py (τ_m, τ_b, R).

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

7.3 Eje de densidad/volumen: ecuación de estado de van der Waals con núcleo duro

La presión asciende de «parámetro» a auténtica variable. El NÚCLEO DURO aporta el volumen propio b, la atracción de enlace (ε·z̄) la presión de cohesión a. En unidades reducidas (punto crítico en (v_r, τ_r, p_r) = (1, 1, 1)):

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

Coexistencia por construcción de Maxwell (cv_scan_09.py, corte de áreas iguales):

Resultado. Por debajo de τ_c se abre el salto de densidad l/g (v_gas/v_liq crece de 2,5 a 11,5, tanto mayor cuanto más frío). Esa es la línea de ebullición — ahora limpiamente separada del pico de c_v de fusión (cap. 5.3), que vive sobre un retículo de densidad fija. Con ello se cierra el punto 3 que quedaba abierto en la versión 08: fusión (orden de rigidez/enlace a densidad fija) y ebullición (salto de densidad/volumen) son ejes distintos del mismo modelo.

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

7.4 Calibración cuantitativa l/g: a → ε, b → D²

Aplicación de las magnitudes reducidas a sustancias reales. Las constantes de van der Waals dependen del punto crítico:

b = k_B·T_c / (8 p_c) (volumen propio, NÚCLEO DURO), a = 27 k_B²·T_c² / (64 p_c) (cohesión, ENLACE).

Así a ∝ ε (la cohesión procede de la energía de enlace) y b ∝ D² (el volumen propio 2D del diámetro del núcleo duro). La escala de energía natural de la rama l/g es

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

Sustancia de anclaje argón (T_c = 150,87 K, p_c = 48,63 bar): b = 5,35·10⁻²⁹ m³ por molécula, a = 3,76·10⁻⁴⁹ J·m³ por molécula², ε_lg/k_B = 509 K y un diámetro molecular D ≈ b^{1/3} = 3,8·10⁻¹⁰ m — con una exactitud del 11 % frente al diámetro de van der Waals del argón (~3,4·10⁻¹⁰ m).

Predicción de la temperatura de ebullición. A presión exterior fija, el punto de ebullición está donde se cumple la presión de saturación reducida p_sat(τ_r)·p_c = p_atm; de ahí T_b = T_c·τ_r,b (curva de saturación de Maxwell de lg_calib_10.py):

*El CO₂ sublima a 1 atm (T = 194,7 K = punto de sublimación).

Resultado. Para las sustancias moleculares (p_c ≈ 30–70 bar) el modelo da en todos los casos T_b/T_c ≈ 0,48–0,51, sistemáticamente ~17–24 % por debajo del valor medido. Eso no es un error del modelo en sentido estricto, sino la conocida distancia de la ecuación de van der Waals a la regla empírica de Guldberg (T_b/T_c ≈ 0,6): la isoterma de van der Waals subestima la presión de saturación a bajas temperaturas. El agua (y los metales, cuyos T_c/p_c de todos modos solo están estimados a grandes rasgos) quedan con p_r,atm < 0,005 fuera del rango sensato de van der Waals (τ_r,b < 0,40) — allí la cohesión ya no es de tipo van der Waals (puentes de hidrógeno, enlace metálico). El orden de magnitud y la ordenación de los puntos de ebullición se aciertan correctamente.

Clasificación de Guldberg como tabla de errores. La distancia se nombra aquí solamente, no se corrige mediante un cambio de EOS (la imagen ideal sigue siendo de tipo van der Waals). La regla empírica de Guldberg (Guldberg 1890) dice que el punto de ebullición se sitúa hacia T_b/T_c ≈ 0,6 (con ligera dependencia de la sustancia, ~0,58–0,64). El modelo reducido de van der Waals da en cambio un valor constante T_b/T_c = τ_r,b ≈ 0,48–0,51:

Lectura. La media del modelo, 0,49, acierta la constante de Guldberg 0,60 con un −18 % aproximadamente; la desviación es sistemática y unilateral (siempre demasiado baja), es decir, una propiedad de la isoterma de van der Waals (presión de saturación demasiado pequeña a τ bajas), no de la imagen molecular. Una mejor EOS de cohesión (por ejemplo Redlich–Kwong) elevaría τ_r,b — pero queda aquí deliberadamente fuera de la imagen ideal.

Presupuesto de error con un cálculo comparativo de Redlich–Kwong (guldberg_rk_15.py). Para cifrar la distancia de van der Waals a la regla de Guldberg (no para sanarla), se repite el mismo test del punto de ebullición con una cohesión de tipo Redlich–Kwong, dependiente de la temperatura, a/√T. En unidades reducidas (punto crítico (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,

con construcción de Maxwell de áreas iguales y la misma presión exterior reducida p_r,atm que arriba. Resultado (un valor por sustancia de anclaje, T_b/T_c = τ_r,b):

Resultado. La cohesión RK dependiente de la temperatura eleva T_b/T_c de vdW ≈ 0,48 a ≈ 0,60 — prácticamente la constante de Guldberg y los valores reales (≈ 0,58–0,63); el error medio frente a las medidas cae de −21 % (vdW) a −1 % (RK). Con ello la causa de la distancia de Guldberg queda claramente cifrada: es justamente la falta de dependencia de la temperatura de la cohesión de van der Waals (presión de saturación demasiado pequeña a τ bajas), no la imagen molecular. Este cálculo comparativo queda deliberadamente fuera de la imagen ideal; el modelo sigue siendo de tipo van der Waals (sin cambio de EOS) — el presupuesto de error solo se nombra, no se incorpora.

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

7.5 Superficie (τ, v) común con punto triple

Los ejes de fusión y de ebullición han vivido hasta ahora separados: la fusión (orden de rigidez/enlace) a densidad fija, la ebullición (van der Waals) mediante el salto de densidad. En un plano común (τ_r, v_r) encajan:

• Cúpula de coexistencia (ebullición): v_liq(τ_r) (rama izquierda) y v_gas(τ_r) (rama derecha) de la construcción de Maxwell, que se encuentran en el punto crítico (τ_r, v_r) = (1, 1).
• Sólido (fusión): banda estrecha justo por encima de la pared dura v_r > 1/3, a densidad alta y casi fija (volumen reducido v_solid ≈ 0,40).
• Línea de fusión: casi vertical en v_solid (densidad fija), a la temperatura de la rigidez de Maxwell (z̄_c = 2d; el cap. 7.8 deriva de ella v_solid sin parámetros — la antigua asignación de la versión 11, «z̄_c ≈ 2,08», queda así corregida).

El punto triple es el punto en que la frontera de coexistencia líquida v_liq(τ_r) alcanza la densidad de empaquetamiento sólido v_solid — por debajo ya no existe rama líquida (sublimación directa):

Resultado. Para v_solid = 0,40 el punto triple está en τ_r,triple ≈ 0,47, p_sat ≈ 0,018, con un volumen de vapor v_gas ≈ 68 (vapor de sublimación muy diluido). La razón τ_triple/τ_c ≈ 0,47 cae dentro de la ventana real (agua 0,42; argón 0,56; CO₂ 0,71). Cuanto más denso el sólido (menor v_solid), tanto más bajo el punto triple reducido y tanto más diluido el vapor coexistente — cualitativamente correcto. Con ello ambas transiciones y el punto triple quedan anclados en una sola superficie (τ, v); el punto crítico se sitúa en (1, 1, 1). La versión 15 añade en el cap. 7.12 una cuarta línea frontera — la transformación sólido(1)/sólido(2) (FCC↔BCC) como dos bandas de densidad sólida casi verticales (v ≈ 0,667 y 0,726) que se transforman una en otra en τ_ss.

Versión anterior en el desarrollo: aquí v_solid = 0,40 todavía estaba impuesto. El cap. 7.8 deriva v_solid sin parámetros del umbral de percolación z̄_c y convierte el punto triple en una auténtica predicción.

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

7.6 Calor latente frente a entalpías de fusión

El calor latente del modelo L ≈ 1,41 ε (cap. 5.3) se traduce a unidades reales mediante la calibración en el punto de fusión ε/k_B = T_m/τ_m (τ_m = 1,20):

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

Eso es —notablemente— casi exactamente la regla de Richards (ΔS_fus ≈ R para sólidos de construcción simple). Comparación con entalpías de fusión medidas (lg_calib_10.py):

Resultado. El modelo acierta en todos los casos el orden de magnitud (kJ/mol) y excelentemente la regla de Richards para metales simples (cobre ×1,00, hierro ×1,28); los gases licuados quedan justo por debajo (metano ×0,94, nitrógeno ×0,87, argón ×0,69). Los cristales iónicos, con puentes de hidrógeno y moleculares (NaCl ×0,37, agua ×0,44, benceno ×0,28, CO₂ ×0,23) tienen entropías de fusión reales claramente mayores (ΔS_exp/R = 2,6–5,0), porque al fundir se liberan grados de orden adicionales (apareamiento iónico, orientación de moléculas poliatómicas) que el modelo de un solo parámetro de orden no capta. El límite queda así claramente nombrado: la imagen ideal vale para sustancias de enlace isótropo y no direccional (sin ajuste).

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

7.7 Unificación de las escalas de energía: ε ↔ ε_lg

Hasta ahora la rama de fusión portaba la energía de enlace ε (τ_m = 1,20) y la rama l/g una escala de cohesión separada ε_lg/k_B = 3,375·T_c (cap. 7.4). Ambas pueden derivarse de un solo parámetro del modelo. En el modelo, la cohesión de van der Waals a surge de la suma de las energías de enlace por molécula (ENLACE): cada molécula porta z_max/2 enlaces de energía ε. Un gas en red de Bragg–Williams en campo medio con atracción de vecinos ε proporciona la relación clásica para la temperatura crítica:

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

Sustituyendo esto en ε_lg = a/b = (27/8)·k_B·T_c (cap. 7.4) se sigue el acoplamiento de ambas escalas:

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

Con ello el pico de fusión, la percolación y el punto crítico están sobre una sola escala ε (unify_scale_12.py):

Predicción τ_c/τ_m = z_max/(4·τ_m): 2D → 1,25, 3D → 2,50. Frente a sustancias reales:

Resultado. La predicción 3D τ_c/τ_m = 2,50 acierta bien la cifra característica real de las sustancias moleculares simples (T_c/T_m ≈ 2,0–2,1); los metales quedan más altos (T_c estimado a grandes rasgos, otro tipo de enlace). La predicción 2D (1,25) es demasiado baja — y τ_c(2D) = 1,50 queda incluso por debajo de τ_b = 1,85 (el punto crítico estaría bajo el umbral de percolación). Eso no es un error de cálculo, sino el conocido papel especial del 2D: allí las escalas de cohesión y de enlace no se ordenan limpiamente, falta una ventana líquida robusta (cf. la fusión 2D continua, cap. 5.5). Solo en 3D (τ_m = 1,20 < τ_c = 3,00) se separan las escalas de forma físicamente sensata.

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

7.8 Punto triple sin parámetros a partir de z̄_c

El cap. 7.5 todavía imponía v_solid = 0,40. La versión anterior en el desarrollo deriva la densidad de empaquetamiento sólido del umbral de percolación. El número de coordinación de un retículo a volumen reducido v_r sigue la relación densidad-coordinación (pared dura v_r = 1/3 ↔ empaquetamiento más denso, z̄ = z_max):

z̄(v_r) = z_max·(1/3)/v_r = z_max/(3 v_r), a la inversa v_solid = z_max/(3·z̄_c).

El sólido pierde su rigidez en cuanto z̄ cae por debajo del umbral z̄_c; la densidad correspondiente es v_solid. El punto triple es entonces el punto en que la frontera de coexistencia líquida v_liq(τ_r) alcanza esa densidad (unify_scale_12.py, 2D, z_max = 6):

Versión mejorada. La densidad de fusión v_solid la fija la rigidez genérica de Maxwell z̄_c = 2d (2D → 4, 3D → 6), no la conectividad de enlaces (justificación en el cap. 5.4a: la fusión es pérdida de rigidez, no de conexidad). Con ello el punto triple fiable reza (rigidity_3d_13.py, 2D):

Resultado 2D. Incluso con el umbral de rigidez (físicamente correcto), el retículo 2D aislado no proporciona ninguna ventana líquida limpia: v_solid = 0,50 da ciertamente el valor real τ_triple/τ_c ≈ 0,77, pero en unidades ε absolutas τ_triple = 0,77·1,50 = 1,16 y queda por tanto por debajo de τ_m = 1,20. El umbral de conectividad (z̄_c = 2,08) empuja incluso v_solid ≈ 0,96 ≈ v_c fuera de la cúpula de coexistencia — punto triple y punto crítico coinciden. Ambos cálculos confirman el papel especial del 2D (cap. 5.5, 7.7): el retículo 2D aislado apenas sostiene una fase líquida autónoma. Es un hallazgo honesto, no una debilidad de cálculo o de modelo en sentido estricto.

El test decisivo es, pues, el cálculo de control 3D (cap. 7.9): allí la misma regla de rigidez (z̄_c = 6 → v_solid = 0,667) libera una ventana limpia. En todo caso el punto triple queda ahora ligado sin parámetros a z̄_c — el parámetro libre v_solid queda eliminado.

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

7.9 Cálculo de control 3D: una ventana líquida limpia sobre una sola escala ε

Comprobación de si el colapso 2D se resuelve en 3D cuando todas las marcas están sobre una sola escala ε (τ_c = z_max/4 del cap. 7.7; v_solid de la rigidez genérica z̄_c = 2d; punto triple donde la densidad líquida reducida de van der Waals v_liq(τ_r) alcanza la densidad del sólido). Resultado (rigidity_3d_13.py):

Resultado. En 3D se abre una ventana líquida limpia: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00 (anchuras: τ_triple − τ_m = 1,63; τ_c − τ_triple = 0,18). El borde superior es estrecho —el punto triple queda justo debajo del crítico— pero el orden es correcto y la fase líquida existe como campo autónomo. En 2D la ventana permanece cerrada (τ_triple = 1,16 < τ_m). Con ello el papel especial del 2D queda documentado por tercera vez de forma independiente y se muestra a la vez que es un efecto dimensional: la misma regla de rigidez sana el modelo en 3D. El estrecho borde superior se corrobora en el cap. 7.10 mediante un barrido de percolación 3D real y un análisis de sensibilidad.

Recuadro de límites del modelo — el papel especial del 2D. Tres cálculos independientes muestran el mismo hallazgo: el retículo 2D aislado apenas sostiene una fase líquida autónoma. 1. τ_c < τ_b (cap. 7.7): el punto crítico (1,50) quedaría por debajo de la percolación de conectividad (1,85) — las escalas de cohesión y enlace no se ordenan limpiamente en 2D. 2. v_solid ≈ v_c (cap. 7.8): la densidad de fusión de conectividad (0,96) casi coincide con la densidad crítica — punto triple y punto crítico colapsan. 3. τ_triple < τ_m (cap. 7.9): incluso con la rigidez de Maxwell correcta el punto triple (1,16) queda por debajo del pico de fusión (1,20) — sin ventana.

Es un hallazgo honesto, no una debilidad de cálculo o de modelo: es el conocido efecto dimensional (fusión 2D continua, tipo KTHNY; cap. 5.5). En 3D la misma regla de rigidez sana el modelo (cap. 7.9/7.10). Si un segundo orden local (cap. 6.4) podría mejorar adicionalmente la imagen 2D sigue siendo una cuestión de extensión, no un defecto de la imagen ideal. Puesto que el modelo debe servir a la ilustración, acepto la debilidad y renuncio a una ampliación axiomática del modelo.

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

7.10 Barrido real de percolación de enlaces en 3D (FCC/BCC) y sensibilidad del borde de la ventana

Los cap. 7.8/7.9 trasladaban mentalmente el valor de conectividad 2D z̄_c = 2,08 al 3D (donde v_solid = 1,92, no físico) y apoyaban la afirmación 3D en la rigidez genérica z̄_c = 2d. La versión 13 examina ahora la conectividad en 3D directamente: percolación de enlaces en FCC (z_max = 12) y BCC (z_max = 8) por Monte Carlo con union-find, bordes periódicos, semilla LCG determinista (rigidity_3d_14.py). Umbral p_c = cruce P_∞(p) = 0,5 (la misma metodología que el estudio de tamaño finito 2D del cap. 5.6):

Debido a la arista de caja moderada (L = 12, Python puro), los valores MC quedan ligeramente por encima de los umbrales de la literatura (sesgo de tamaño finito, como en 2D desde abajo o desde arriba según el estimador); la ordenación y el orden de magnitud son correctos (FCC más denso → p_c menor que el BCC).

Lectura (confirmación 3D independiente). Los umbrales de conectividad 3D quedan en z̄_c(3D) ≈ 1,5–1,6 — incluso más bajos que el valor 2D 2,08 y muy por debajo de la rigidez de Maxwell 2d = 6. La «densidad de fusión de conectividad» que de ahí se sigue, v_solid = 1/(3 p_c) ≈ 1,7 (BCC) o 2,6 (FCC), queda muy fuera de la cúpula de coexistencia (v_liq,max ≈ 0,94). Con ello se muestra directamente —sin trasladar el valor 2D al 3D— que en 3D la conectividad con mayor razón no puede fijar la densidad de fusión; esta debe venir de la rigidez (z̄_c = 6 → v_solid = 0,667). Los dos umbrales están aún más claramente separados en 3D que en 2D. El auténtico τ_b(3D) (temperatura a la que z̄(τ) alcanza el bajo umbral de conectividad 1,5) quedaría, pues, por encima del valor 2D τ_b, muy cerca del punto crítico — coherente con la imagen «desintegración de la conectividad = ebullición».

Sensibilidad del estrecho borde superior de la ventana (τ_c − τ_triple ≈ 0,18). Numéricamente en torno a v_solid = 0,667 (±2 %) resulta dτ_triple/dv_solid ≈ 1,55. La ventana se cerraría (τ_triple = τ_c = 3,00) solo si v_solid subiera en ≈ 0,11 (hasta ≈ 0,78), es decir, si el umbral de rigidez 3D efectivo estuviera notablemente por encima de 2d = 6 (z̄_c ≳ 7,7). Para el recuento genérico de Maxwell (z̄_c = 2d = 6) el orden τ_m < τ_triple < τ_c sigue siendo, pues, robusto; el borde superior es realmente estrecho, pero no está en peligro. El cambio discreto de red sólido(1)→sólido(2) (FCC↔BCC) está cuantificado en la versión 14 (cap. 6.3); el auténtico barrido en τ de la cinética de enlace en 3D sigue en el cap. 7.11.

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

7.11 Barrido real en τ de la cinética de enlace en 3D: τ_b(3D) y la línea de fusión sobre una escala ε

Los cap. 7.8–7.10 apoyaban la afirmación 3D en la rigidez genérica (z̄_c = 2d = 6) y en el barrido estático de conectividad (z̄_c(3D) ≈ 1,55). Quedaba abierto el auténtico barrido en τ de la cinética de enlace en 3D: ¿dónde, sobre el eje de temperaturas, alcanza el z̄(τ) autoconsistente ese bajo umbral de conectividad (= el auténtico τ_b(3D)), y dónde está el pico de fusión? La versión 14 calcula el z̄(τ) autoconsistente (cap. 4.2, x = 1/(1 + G·e^{−(ε+K·x)/τ})) directamente sobre FCC (z_max = 12) y BCC (z_max = 8) con G = z_max y la misma cooperatividad J = 0,55 que en el caso de referencia 2D; n_blk(τ) procede de un Monte Carlo de enlaces FCC/BCC (union-find), y

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄ − (K/2)·x² (3D, con término cooperativo),

c_v = dU/dτ (bondkinetics_3d_15.py, semilla LCG determinista; n_blk suavizado, ya que solo el débil término de rotación de bloque depende de él). Se leen tres marcas: el pico de c_v τ_m (fusión), la marca de rigidez τ_rig (z̄(τ) cruza 2d = 6) y el auténtico τ_b(3D) (z̄(τ) cruza el umbral de conectividad z̄_c(3D) del cap. 7.10):

Resultado. Dos hallazgos, ambos coherentes con la imagen hasta aquí:

1. La fusión queda baja y está doblemente documentada. El pico termodinámico de fusión τ_m (FCC 1,70, BCC 1,40) coincide prácticamente con la marca de rigidez τ_rig (z̄ = 6; FCC 1,73, BCC 1,35) — la misma temperatura a la que la red pierde su rigidez. Eso ancla la línea de fusión 3D en la escala ε: T_m = (ε/k_B)·τ_m, es decir, por ejemplo para el argón (T_m = 83,8 K) ε/k_B = 83,8/1,70 ≈ 49 K. Como K = J·z_max es mayor en 3D que en 2D (6,60 en vez de 3,30), el pico de fusión 3D queda más alto que el valor 2D 1,20 — no es contradicción, sino la cooperatividad más fuerte del empaquetamiento más denso.
2. El auténtico τ_b(3D) se sitúa en el punto crítico. Como el umbral de conectividad 3D es tan bajo (≈ 1,55), z̄(τ) lo alcanza solo muy arriba: τ_b(3D) = 3,18 (FCC) y 2,85 (BCC) — en, o ligeramente por encima de, τ_c = z_max/4 (3,00 y 2,00). El cúmulo gigante de enlaces sobrevive, pues, a toda la fase líquida y se desintegra solo en el punto crítico. Esa es exactamente la predicción del cap. 7.10 («el auténtico τ_b(3D) quedaría por encima del valor 2D, muy cerca del punto crítico») y confirma de forma independiente la imagen «ebullición = desintegración de la conectividad»: en 3D la materia permanece enlazada de un extremo a otro hasta muy cerca del punto crítico.

Ambos hallazgos corroboran la separación de umbrales del cap. 5.4a: fusión ↔ rigidez (bajo, τ_m ≈ τ_rig), ebullición ↔ conectividad (alto, τ_b(3D) ≈ τ_c). La densidad de fusión viene, pues, de la rigidez (v_solid = 0,667), no de la conectividad — los dos umbrales están separados sobre el eje de temperaturas en 3D más que nunca antes en el modelo.

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

7.12 La línea de transformación sólido(1)→sólido(2) en la superficie (τ, v)

El cap. 6.3 proporcionaba la transformación de red FCC↔BCC como punto τ_ss (a p = 0). Mediante la presión resulta de ahí una línea (polymorph_line_16.py) que pertenece como cuarta línea frontera a la superficie (τ, v) común de los cap. 7.5/7.9 — junto a la cúpula de coexistencia (ebullición), la línea de fusión y la pared dura.

Salto de volumen. Los dos sólidos se sitúan a densidad casi fija; su diferencia de volumen se sigue de las fracciones de empaquetamiento (v ∝ 1/η): η_FCC = 0,7405, η_BCC = 0,6802. Con v_FCC = 0,667 (cap. 7.9) resulta v_BCC = v_FCC·η_FCC/η_BCC = 0,726, es decir

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

El BCC más abierto tiene el mayor volumen — la recristalización FCC→BCC va unida a un aumento de volumen (el mismo enunciado que en la tira de cómic «recristalización»).

Clausius–Clapeyron. Con ΔS = S_BCC − S_FCC (fase de alta menos fase de baja temperatura), la pendiente reza

dp/dτ = ΔS/Δv.

Ambas magnitudes son positivas (el BCC es la fase de alta temperatura, tiene por tanto más entropía y el mayor volumen), luego dp/dτ > 0: una presión mayor eleva la temperatura de transformación, porque favorece al FCC más denso. Para la referencia θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, L = 0,667 ε):

Línea en (τ, p) para la referencia (creciente): p = 0 → τ_ss = 1,33; p = 2 → 1,57; p = 4 → 1,81; p = 6 → 2,04 (unidades reducidas).

En la superficie (τ, v) la transformación aparece como dos bandas de densidad sólida casi verticales — sólido(1) = FCC en v ≈ 0,667 (más denso, T baja) y sólido(2) = BCC en v ≈ 0,726 (más abierto, T alta) — que se transforman una en otra en τ_ss; por encima sigue, en τ_m, la fusión. El orden completo sobre el eje τ es, pues, τ_ss < τ_m < τ_c ≤ τ_b(3D) (cap. 7.11).

Unidades reales. Mediante la calibración en el punto de fusión ε/k_B = T_m/τ_m (τ_m = 1,70, FCC 3D) el escalón de c_v y el calor de transformación se hacen tangibles. El escalón de c_v es independiente de la calibración: |Δc_v| = 1,0 k_B = 1 R = 8,31 J/(mol·K) = 0,0083 kJ/(mol·K). El calor de transformación ΔH_ss = L·R·(ε/k_B) (referencia L = 0,667 ε):

Los calores de transformación son pequeños y positivos (del orden de 0,2–0,9 kJ/mol) y quedan por tanto claramente por debajo de las entalpías de fusión (cap. 7.6, varios kJ/mol) — exactamente la imagen esperada de una transformación polimórfica débil con un pequeño escalón de c_v. (La razón θ de 0,60 queda deliberadamente fuera de la imagen ideal estricta; el cap. 6.6 muestra que la geometría pura daría θ ≈ 0,82–0,94 y la transformación quedaría entonces cerrada — la línea es la consecuencia si-entonces de una fase BCC blanda.)

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

7.13 Línea de transformación sólido(1)→sólido(2) en el diagrama (τ, p), calibrada con hierro

El cap. 7.12 lleva la transformación como línea (τ, v); la versión 16 la añade como línea en la imagen (τ, p) y la calibra con un polimorfo FCC↔BCC real (solidsolid_pt_17.py).

Campo de presión de la EOS de van der Waals con núcleo duro. En las dos densidades sólidas v_FCC = 0,667 y v_BCC = 0,726 (cap. 7.12), la EOS reducida p_r = 8τ_r/(3v_r − 1) − 3/v_r² (cap. 7.3) da en todos los casos presiones positivas; a igual τ_r el FCC más denso porta la presión mayor (τ_r = 1: p_r = 1,25 frente a 1,10). La compresión favorece, pues, al FCC — coherente con la pendiente positiva de Clausius–Clapeyron dp/dτ = ΔS/Δv > 0.

Línea y punto de transformación. Para la referencia θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, dp\/dτ = 8,47 en unidades reducidas p\ = p·D³/ε) la línea crece según p\(τ) = 8,47·(τ − 1,333): p\ = 0 → τ = 1,33; p\ = 2 → 1,57; p\ = 4 → 1,81; p\ = 6 → 2,04. En **p\_tri ≈ 3,1 la línea sólido-sólido encuentra la marca de fusión τ_m = 1,70 — un punto de transformación sólido(1)-sólido(2)-líquido** en el plano (τ, p) (análogo del punto triple).

Calibración con hierro (ε/k_B = T_m/τ_m = 1811 K / 1,70 = 1065 K). El hierro tiene dos transiciones FCC↔BCC de signo opuesto — justamente el caso de prueba:

Para γ→δ (τ_ss = 1667/1065 = 1,565 → θ_BCC/θ_FCC ≈ 0,61) el modelo acierta signo y orden de magnitud: ΔH_ss ≈ 4,2 kJ/mol (en realidad ≈ 0,84 kJ/mol — el orden de magnitud kJ/mol, con el valor unas 5× demasiado alto por la aproximación de umbral abrupta), escalón de c_v 1 R = 8,31 J/(mol·K) sin parámetros, dT/dp ≈ +245 K/GPa. La transición α→γ, en cambio, tiene la fase densa (γ) como fase de alta temperatura — eso no lo puede dar la pura geometría de empaquetamiento (más denso = más bajo = T baja); el hierro α se estabiliza magnéticamente. El hierro confirma así ambas caras del hallazgo honesto: la topología geométrica (γ→δ) está dentro de la imagen ideal, la anomalía magnética (α→γ) fuera — el mismo «fuera» que los cap. 6.6/6.8 cifran con θ ≈ 0,87…1,00 > 0,61.

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

7.14 Redlich–Kwong como rama separada de la imagen ideal

La distancia de Guldberg de la rama de van der Waals (cap. 7.4, T_b/T_c ≈ 0,48 en vez de ≈ 0,60) allí solo se nombraba. A petición expresa del usuario, una versión mejorada lleva la cohesión de Redlich–Kwong como variante autónoma de la imagen ideal con calibración propia, coexistencia propia y sección de errores propia (rk_branch_17.py) — la rama de van der Waals sigue siendo la opción por defecto; toda la restante física del modelo (fusión, percolación, sólido(1)→sólido(2)) es indiferente a la elección de la EOS de cohesión.

Escala de energía propia. De a_RK = Ω_a·R²T_c^{2,5}/p_c y b_RK = Ω_b·R T_c/p_c se sigue la energía de cohesión de la rama:

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

mayor, porque la cohesión dependiente de la temperatura a/√T pliega el factor √T_c. El diámetro que se sigue de b_RK acierta los valores reales de modo similar al de la rama vdW (argón D ≈ 3,34·10⁻¹⁰ m).

Coexistencia propia (reducida, punto crítico (1, 1, 1), Z_c = 1/3; p_r = 3τ_r/(v_r − β) − A/(√τ_r·v_r(v_r + β)), β = 0,25992, A = 3,84732):

Sección de errores propia (punto de ebullición en p_r,atm = p_atm/p_c, T_b/T_c = τ_r,b):

Resultado. La rama RK eleva T_b/T_c de vdW ≈ 0,48 a ≈ 0,60 y acierta con ello prácticamente la constante de Guldberg y los valores reales (error medio −1 % frente a vdW −21 %). Se sostiene como segunda variante de la imagen ideal, claramente separada, junto a la opción por defecto de van der Waals; cuál de las dos ramas es «más correcta» es cuestión de la exactitud exigida, no de la estructura del modelo — el núcleo del modelo 3 (z̄ como parámetro de orden) permanece idéntico en ambas.

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

7.15 sólido(1)→sólido(2) en otros polimorfos: titanio y control con gas noble

El cap. 7.13 calibra la línea (τ, p) sólido(1)→sólido(2) con hierro. La versión mejorada asegura signo y orden de magnitud en un segundo polimorfo, magnéticamente discreto, y añade un control negativo (solidsolid_titanium_18.py). Calibración y salto de volumen son como en el cap. 7.13 (ε/k_B = T_m/τ_m, τ_m = 1,70; Δv = +8,9 % denso→abierto de la geometría de empaquetamiento); el hcp se aproxima como el FCC mediante z = 12 (la diferencia fcc/hcp solo se nombra, no se resuelve).

El titanio, con α(hcp) → β(bcc) a 1155 K, tiene exactamente la topología del modelo: la fase densa (hcp, z = 12) es la de baja temperatura y al calentar se vuelve abierta (bcc, z = 8), luego ΔS > 0 y Δv > 0 → dp/dτ > 0 — el mismo signo que el hierro γ→δ, pero sin anomalía magnética. Dos polimorfos denso→abierto independientes en comparación:

Ambos aciertan el signo (dp/dτ > 0, la presión favorece la fase densa) y el orden de magnitud (ΔH ~ kJ/mol, escalón de c_v sin parámetros 1 R = 8,31 J/(mol·K), dT/dp ~ 10²·K/GPa). El valor de ΔH_ss queda en ambos casos sistemáticamente como un factor demasiado alto — la aproximación de umbral abrupta de la comparación de μ sobreestima el calor latente de transformación; es un límite del modelo nombrado, no un error de signo.

El argón es el control negativo: el argón sólido es FCC y sigue siendo FCC hasta la fusión (83,8 K) — a presión normal no existe bcc estable. Justo eso predice la imagen ideal del potencial de pares puro: θ_BCC/θ_FCC ≈ 0,9…1,0 > 0,61 y un BCC LJ inestable a cizalla (cap. 6.9) ⇒ ninguna transición sólido(1)→sólido(2). La contraprueba confirma la estructura del hallazgo: la imagen ideal genera una línea sólido(1)→sólido(2) solo cuando aportes no geométricos (electrónicos/magnéticos, como en el hierro y el titanio) estabilizan primero el BCC. El titanio aporta la prueba positiva, el argón la negativa del mismo marco honesto de los cap. 6.8/6.9.

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

8 — Comparación bibliográfica: concepciones de la termodinámica estática

Resumen. El capítulo 8 sitúa el modelo 3 entre las concepciones establecidas de la termodinámica estática:

• 8.1 Tabla de conjunto: once modelos (gas ideal, van der Waals, Redlich–Kwong, desarrollo del virial, Einstein, Debye, Lennard-Jones / esferas duras, teoría celular, Ising / gas en red, percolación, Flory–Huggins) con su(s) fase(s) y su idea básica estadístico-termodinámica.
• 8.2 Los modelos uno a uno: un párrafo cada uno con rasgos, aplicación y límites, comenzando por el enlace a la fuente; las fuentes de libre acceso están archivadas en model3/links/, y la teoría celular (de pago) solo se cita.
• 8.3 Contraposición: el modelo 3 es un modelo de enlace cooperativo de dos estados con el parámetro de orden z̄; su pariente más cercano es Ising / gas en red, la ebullición se toma directamente de la percolación y la imagen de los fonones (Einstein/Debye) sostiene la rigidez de fusión. El modelo 3 engrapa estas tres tradiciones sobre una escala ε común.

Este capítulo sitúa el modelo 3 entre las imágenes establecidas de la termodinámica de equilibrio («estática») de sólidos, líquidos y gases. Primero un panorama tabular (8.1), después un breve párrafo en prosa por modelo con rasgos típicos, ámbito de aplicación y límites — cada párrafo comienza con el enlace a la fuente (8.2). Finalmente una contraposición esquemática con el modelo 3 (8.3). Las fuentes de libre acceso están archivadas en model3/links/; el script de verificación model3/scripts/check_links_16.py mantiene consistentes cita y archivo mediante SHA-256.

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

8.1 Tabla de conjunto

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

8.2 Los modelos uno a uno

https://de.wikipedia.org/wiki/Kinetische_Gastheorie — Gas ideal / teoría cinética de gases. Partículas sin volumen propio y sin interacción salvo los choques elásticos; la presión surge puramente de la transferencia de momento. Proporciona la ecuación de estado pV = Nk_BT, la distribución de velocidades de Maxwell–Boltzmann y ⟨E⟩ = (3/2)k_BT. Aplicación: gases diluidos a alta temperatura y baja presión. Límites: sin condensación, sin fase líquida ni sólida, sin transición de fase — la interacción está ausente por construcción.

https://de.wikipedia.org/wiki/Van-der-Waals-Gleichung — Ecuación de van der Waals. (p + a/v²)(v − b) = k_BT completa el gas ideal con un volumen de exclusión b (repulsión) y un término de cohesión a (atracción). Es la primera ecuación de estado con una transición gas-líquido, punto crítico y construcción de Maxwell. Aplicación: imagen cualitativa de la condensación, magnitudes críticas. Límites: cuantitativamente imprecisa cerca del punto crítico (exponentes de campo medio), sin fase sólida.

https://de.wikipedia.org/wiki/Redlich-Kwong-Gleichung — Ecuación de estado de Redlich–Kwong. EOS cúbica con término de atracción dependiente de la temperatura a/[√T · v(v + b)]. El escalado en √T mejora notablemente los equilibrios vapor-líquido frente a van der Waals (cf. cap. 7.4, donde eleva la razón T_b/T_c). Aplicación: cálculos técnicos de fase gaseosa y evaporación. Límites: parametrizada empíricamente, sigue sin fase sólida; cerca del punto crítico también de tipo campo medio.

https://de.wikipedia.org/wiki/Virialentwicklung — Desarrollo del virial. Desarrollo sistemático pV/Nk_BT = 1 + B(T)/v + C(T)/v² + … Los coeficientes del virial se siguen del potencial intermolecular (integrales de cúmulo de la mecánica estadística) y están por tanto microscópicamente fundados. Aplicación: gases hasta densidad media, referencia para otras EOS. Límites: la serie no converge para líquidos densos o sólidos; los coeficientes superiores son costosos de calcular.

https://de.wikipedia.org/wiki/Einstein-Modell — Sólido de Einstein. El cristal es un conjunto de 3N osciladores armónicos independientes de igual frecuencia ω. Explica la caída de la capacidad calorífica c_v por debajo de la temperatura característica θ_E = ħω/k_B y el límite de Dulong–Petit 3Nk_B a alta temperatura. Aplicación: la descripción más sencilla de las capacidades caloríficas de sólidos. Límites: ley errónea a baja temperatura (exponencial en vez de T³), sin modos acústicos, solo la fase sólida.

https://de.wikipedia.org/wiki/Debye-Modell — Modelo de Debye. Trata el sólido como un continuo elástico con un espectro de fonones hasta la frecuencia de Debye; c_v ∝ T³ a baja temperatura, Dulong–Petit a alta — la temperatura de Debye θ_D resume la dureza de la red. Aplicación: capacidades caloríficas cuantitativamente buenas de cristales simples en todo el rango de T. Límites: aproximación de continuo (ignora la dispersión de las redes reales), solo la fase sólida, sin fusión.

https://de.wikipedia.org/wiki/Lennard-Jones-Potential — Potencial de Lennard-Jones / esferas duras. Potencial de pares 4ε[(σ/r)¹² − (σ/r)⁶] con núcleo repulsivo y atracción débil; el modelo límite de esferas duras conserva solo la repulsión. Es la base de la dinámica molecular y del Monte Carlo y describe —numéricamente— las tres fases junto con el punto triple. Aplicación: gases nobles, sistema de referencia genérico. Límites: aditividad por pares, evaluable en el fondo solo por simulación o teoría de perturbaciones (sin EOS cerrada).

https://doi.org/10.1098/rspa.1937.0210 — Teoría celular del líquido (Lennard-Jones–Devonshire). Cada molécula se mueve en una celda (jaula) que forman sus vecinos en campo promediado; el líquido se describe mediante la función de partición celular y proporciona una EOS del líquido así como un enfoque de la fusión. Aplicación: líquidos simples, eslabón entre sólido y gas. Límites: sobreestima el orden, desatiende vacantes/huecos (esto lo corrigen teorías de huecos posteriores). (De pago — solo citada, no archivada localmente; accesibilidad comprobada en el script de verificación.)

https://de.wikipedia.org/wiki/Ising-Modell — Modelo de Ising / gas en red. Sobre una red los sitios portan una variable de dos estados (espín ↑/↓ u ocupado/vacío) con acoplamiento a vecinos próximos. El gas en red es matemáticamente isomorfo al modelo de Ising y describe la condensación gas↔líquido como transición cooperativa orden-desorden con punto crítico. Aplicación: fenómenos críticos, transiciones cooperativas. Límites: discretización en red, sin energía cinética explícita; en aproximación de campo medio solo cualitativo.

https://de.wikipedia.org/wiki/Perkolationstheorie — Teoría de la percolación. Modelo puramente geométrico: sitios o enlaces se ocupan al azar con probabilidad p; por encima de un umbral p_c existe un cúmulo que atraviesa. Describe la gelificación, el inicio de la conductividad y en general las transiciones de conectividad. Aplicación: formación de redes, desintegración/constitución de la conexidad. Límites: sin energía ni temperatura intrínsecas — el acoplamiento a la termodinámica debe añadirse desde fuera.

https://de.wikipedia.org/wiki/Flory-Huggins-Theorie — Teoría de Flory–Huggins. Modelo en red para disoluciones y polímeros: la entropía de mezcla de la colocación de cadenas más un parámetro de interacción χ gobiernan la miscibilidad y la separación de fases. Aplicación: disoluciones de polímeros, diagramas de fases de mezclas. Límites: campo medio, red incompresible, χ casi siempre empírico — de relevancia solo limitada para imágenes sólido/líquido/gas de un componente, aquí se mantiene como referencia de red.

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

8.3 Semejanzas y diferencias con el modelo 3

El modelo 3 es un modelo de enlace cooperativo de dos estados: el parámetro de orden es el número de coordinación medio z̄ = z_max·x, autoconsistente mediante la cooperatividad K = J·z_max con la degeneración del estado roto G = z_max; la fusión como umbral de rigidez (fonones), la ebullición como desintegración de la conectividad (percolación) en el punto crítico, transición de primer orden en cuanto ln G > 2 + 4ε/K.

En resumen: el modelo 3 engrapa tres tradiciones por lo demás separadas — la imagen del gas en red / de Ising (transición cooperativa de dos estados), la imagen de la percolación (conectividad → ebullición) y la imagen de los fonones de Einstein/Debye (rigidez → fusión) — y las calibra mediante un único parámetro de energía ε sobre una escala τ común. Van der Waals y la teoría celular proporcionan la contraprueba continua (EOS y jaula líquida), frente a la cual puede compararse la imagen discreta de los enlaces.

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

9 — Apéndice

Resumen. El apéndice reúne los soportes formales:

• 9.1 Bibliografía: ordenada por nombre de autor, con enlaces a las fuentes archivadas libremente (comparación SHA-256 por check_links_16.py; la teoría celular solo se cita, por el muro de pago).
• 9.2 Apoyo de la IA: la proporción de innovación IA/autor estimada por la IA para cada cambio sustancial (con desglose de la parte del autor en estrategia/estructura/redacción y un contracálculo), así como una relación de los pasos individuales con coste (estimado), duración y distancia respecto al registro anterior (suma ≈ 160 USD, estimaciones claramente señaladas).
• 9.3 Glosario de abreviaturas y símbolos: los símbolos de uso frecuente (τ, ε, z̄, K, G, θ, …), que por ello pueden aparecer excepcionalmente en el texto corrido sin descripción previa.
• 9.4 Scripts de cálculo y cobertura de pruebas: todas las cifras proceden de scripts de Python reproducibles (proporción de programación por IA 100 %, cobertura de pruebas 0 % en todos los casos), junto con la abierta «declaración embarazosa» de que el autor no comprobó él mismo el código ni los resultados.

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

9.1 Bibliografía

Ordenada por nombre de autor. Las fuentes de libre acceso están archivadas en model3/links/; su concordancia con los enlaces citados la comprueba model3/scripts/check_links_16.py (SHA-256). La teoría celular de Lennard-Jones y Devonshire solo se cita, por el muro de pago (accesibilidad comprobada en el script de verificación).

• Broadbent, S. R.; Hammersley, J. M. (1957): Percolation processes I. Crystals and mazes. Proc. Cambridge Phil. Soc. 53, 629–641. — Panorama/archivo: Teoría de la percolación
• Clausius, R.; Maxwell, J. C.; Boltzmann, L. (1857–1872): fundamentación de la teoría cinética de gases. — Panorama/archivo: Teoría cinética de gases
• Debye, P. (1912): Zur Theorie der spezifischen Wärmen. Ann. Phys. 39, 789–839. — Panorama/archivo: Modelo de Debye
• Einstein, A. (1907): Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Ann. Phys. 22, 180–190. — Panorama/archivo: Modelo de Einstein
• Flory, P. J.; Huggins, M. L. (1941–1942): teoría en red de las disoluciones de polímeros (parámetro χ). — Panorama/archivo: Teoría de Flory–Huggins
• Ising, E.; Lenz, W. (1920–1925): Beitrag zur Theorie des Ferromagnetismus. Z. Phys. 31, 253–258 (el gas en red es matemáticamente isomorfo). — Panorama/archivo: Modelo de Ising
• Kamerlingh Onnes, H. (1901): forma del virial de la ecuación de estado. — Panorama/archivo: Desarrollo del virial
• Lennard-Jones, J. E. (1924): On the determination of molecular fields. Proc. R. Soc. Lond. A 106, 463–477. — Panorama/archivo: Potencial 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 (De pago — solo citado, no archivado localmente.)
• Redlich, O.; Kwong, J. N. S. (1949): On the thermodynamics of solutions. V. An equation of state. Chem. Rev. 44, 233–244. — Panorama/archivo: Ecuación de Redlich–Kwong
• van der Waals, J. D. (1873): Over de Continuïteit van den Gas- en Vloeistoftoestand (tesis, Leiden). — Panorama/archivo: Ecuación de van der Waals

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

9.2 Apoyo de la IA

Este desarrollo de modelo surgió de una interacción entre persona e IA: el autor (el usuario de la IA) determina en qué se trabaja (preguntas rectoras, decisiones, deseos sobre el modelo); la IA (Claude Code) determina en gran medida cómo — derivaciones, scripts, estructura y texto. Las dos tablas siguientes las estima la propia IA.

(a) Cambios sustanciales de contenido y proporción de innovación estimada. La proporción de innovación pondera las ideas rectoras al 45 %, los trabajos estructurantes al 40 %, los trabajos redaccionales al 15 %. Las ideas rectoras (el nivel del en qué) proceden predominantemente del autor; los trabajos estructurantes y redaccionales (el nivel del cómo), predominantemente de la IA; allí donde la IA aportó también la idea física rectora (por ejemplo la lectura de la ebullición por percolación, el modelo de muelles para θ), su parte aumenta. Valores en porcentaje, IA + autor = 100. Las tres columnas adicionales desglosan la parte del autor según los tres niveles de acción: estrategia (ideas rectoras), estructura (derivaciones/scripts/arquitectura) y redacción (trabajo sobre el texto). Su media ponderada al 45 % / 40 % / 15 % reproduce la columna Autor % como contracálculo (model3/scripts/author_share_19.py, desviación máxima 0,95 puntos porcentuales).

(b) Pasos individuales: coste (estimado), duración, distancia respecto al registro anterior. Fecha de creación es la hora de creación del archivo de registro (mtime del archivo). Δ al registro anterior es la distancia temporal respecto a la entrada anterior (medida del esfuerzo del usuario entre pasos, incluidas pausas y /compact). Coste (estimado) son las estimaciones documentadas por registro (precios de lista de Anthropic, ± 30 %). Duración (activa) es una estimación de la IA del tiempo de trabajo puro; en pasos separados por /compact puede ser mayor que Δ. Un asterisco marca una presunta reconstrucción o reelaboración de contenidos existentes.

Suma de los costes documentados (estimados, pasos 01–21): ≈ 129 USD (suma de las estimaciones individuales); con el paso 22 (16d, ≈ 10 USD) ≈ 139 USD, el paso 23 (V16, ≈ 12 USD) ≈ 151 USD y el paso 24 (V17, ≈ 9 USD) ≈ 160 USD. Todos los valores son estimaciones claramente señaladas (no hay medición ccusage disponible); los detalles de los registros están en model3/logs/dev_log_NN.md.

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

9.3 Glosario de abreviaturas y símbolos

Por principio, este texto antepone a cada símbolo un término descriptivo (así «capacidad calorífica c_v» en vez de solo «c_v»). Las siguientes magnitudes de uso frecuente constituyen la excepción: pueden aparecer en el texto corrido sin descripción previa porque aquí quedan explicadas de una vez por todas. Todas las temperaturas τ… son reducidas (en unidades de la energía de enlace ε).

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

9.4 Scripts de cálculo y cobertura de pruebas

Todas las cifras y tablas mencionadas en el texto proceden de scripts de Python reproducibles en model3/scripts/ (cada uno con una salida en bruto …_output.txt). La proporción de programación por IA indica qué parte del código de los scripts escribió la IA; la cobertura de pruebas se refiere a pruebas automatizadas (pruebas unitarias, comprobaciones assert). Es baja en todos los casos: los scripts son cálculos únicos sin andamiaje de pruebas; su única garantía es la plausibilidad de la salida en bruto.

El autor confiesa abiertamente (declaración embarazosa): no ha comprobado, recalculado ni verificado técnicamente ni el código fuente de estos scripts ni las cifras, tablas y curvas producidas con ellos. Confió aquí ciegamente en la competencia de la IA, porque lo único que estaba en primer plano era el carácter cualitativo y la plausibilidad de la concepción del modelo. Los valores deben, pues, tomarse con precaución (véase también la nota detallada en el cap. 10.4).

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

10 Autocrítica

Tengo la sensación de que en la tabla del capítulo 9.2 (a) la parte de trabajo del autor, estimada por la IA, es en los tres ámbitos casi siempre demasiado alta en al menos 10–20 puntos porcentuales. Tengo la sensación de haber aportado solo la idea básica, una o dos interpretaciones e importantes refinamientos en la visualización gráfica.

Los capítulos (matemáticos) 4–9 los propuso, investigó e implementó la IA. Peor aún: durante el desarrollo tuve la sensación de que la IA estaba obsesionada con elaborar una estadística basada en el número de coordinación. Mis verdaderas ideas básicas —que no solo los átomos sino también las moléculas pueden formar enlaces y que las rupturas de enlace dentro de un líquido podrían seguir regularidades fractales— la IA solo pudo asumirlas con dificultad, intentando ignorar la terminología. Por eso esas ideas básicas siguen tratándose de forma muy secundaria en el texto. No lo reescribí después, porque al reflexionar sobre las transiciones de fase partiendo de los axiomas fundamentales uno tiene esas ideas presentes de manera natural. Además, la fijación o reducción de las ideas al número de coordinación como magnitud de medida no es el peor camino para matematizar mi representación figurada, pues a posteriori resultó bastante sostenible para la comprobación de plausibilidad — si bien solo hice una comprobación de plausibilidad muy somera de los resultados de la IA y no realicé examen alguno del código Python. En mi descargo, dígase: lo que me importaba fue siempre únicamente una representación figurada de las relaciones moleculares en los distintos estados de agregación. A la adulación de la IA, que quería confirmar matemáticamente la representación, no le dije que no. La manipulación por parte de la IA tuvo, pues, éxito conmigo, porque examiné los capítulos mencionados solo del modo más superficial y los adopté con fe. (Esta experiencia sobre mí mismo refuerza mi desconfianza acerca de si un «human in the loop» es realmente eficaz para controlar la IA.)

Este texto se elaboró con apoyo de IA (Claude Code 2.1.161, modelo Opus 4.8 / claude-opus-4-8[1m], contexto 1M; fecha 2026-09-05). — Elaborado con apoyo de IA.

Nota sobre la marca de agua. El texto maestro alemán (de.modell_final_02.md) lleva en las palabras clave una marca de agua hecha de homoglifos griegos (Α→A, Β→B, Ε→E, Η→H, Ι→I, Κ→K, Μ→M, Ν→N, Ο→O, Ρ→P, Τ→T, Υ→Y, Χ→X, Ζ→Z). Como ese recurso está ligado a la grafía alemana de las palabras clave, esta traducción lleva las palabras clave en mayúsculas normales: VUELO LIBRE (FREIFLUG), MULTITUD (VIELHEIT), ENLACE (BINDUNG), TRÍADA (DREIKLANG), NÚCLEO DURO (HARTKERN), CONSERVACIÓN (ERHALTUNG), GRIETA (RISS), ÁGIL (FLINK).