Termodinâmica estatística da substância pura ideal e dos seus estados de agregação

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

Resumo geral

O modelo 3 é um modelo qualitativo, de mecânica estatística, de uma substância ideal que deriva os três estados de agregação — sólido, líquido e gasoso — bem como as suas transições a partir de uma única imagem de ligação, a das ligações moleculares. O princípio de base está em considerar três reservatórios de energia para a energia interna: a energia de vibração nas ligações moleculares, a energia de translação das moléculas e a energia cinética cooperativa na rotação de micelas, forma idealizada dos fluxos internos de matéria que compensam as flutuações internas de densidade. Uma grandeza fundamental importante neste contexto é o número de coordenação médio z̄ = z_max·x, que se estabelece de forma autoconsistente através de uma ligação cooperativa de dois estados (do tipo Ising, cooperatividade K = J·z_max, degenerescência do estado rompido G = z_max). Dela seguem três transições numa escala τ comum, calibrada por um único parâmetro de energia ε: a fusão como limiar de fonões/rigidez (z̄ atinge a rigidez), a ebulição como colapso da conectividade (percolação) no ponto crítico e uma transição opcional sólido(1)→sólido(2) (FCC↔BCC) com o degrau de c_v sem parâmetros de 1,0 k_B. Uma transição torna-se de primeira ordem logo que ln G > 2 + 4ε/K. Aqui a degenerescência do estado rompido G é o número de contactos vizinhos que podem romper-se (G = z_max, o número de coordenação máximo), a energia de ligação ε é a energia ganha por contacto vizinho e a cooperatividade K = J·z_max é o produto da força de acoplamento J, do tipo Ising, por z_max; a transição torna-se portanto de primeira ordem tanto mais cedo quanto maior for a multiplicidade de contactos G e quanto mais forte for a cooperatividade K face à energia de ligação ε.

O artigo conduz dos sete axiomas (capítulo 1) através da análise das contradições (capítulo 2) e de uma descrição figurada com tiras estáticas (capítulo 3) até à armação de mecânica estatística — ligação cooperativa, funções de partição e capacidades térmicas (capítulos 4–5) —, à plausibilidade matemática em 3D (capítulo 6) e finalmente à calibração e à equação de estado (capítulo 7). As ideias são complementadas com o varrimento 3D completo de c_v da fusão de primeira ordem (Δz̄, calor latente ΔU, deslocamento do pico), um modelo de molas de Einstein para a rigidez fonónica θ(rede) e a linha de transformação sólido(1)→sólido(2) via Clausius–Clapeyron. θ(rede) é depois refinado numa densidade de estados de fonões do tipo Debye; a transformação sólido(1)→sólido(2) é ainda introduzida como linha no plano (τ, p) (calibrada com ferro) e é acrescentado um ramo Redlich–Kwong separado ao lado da opção por omissão de van der Waals. θ(rede) é além disso testado numa rede de Lennard-Jones completa (FCC estável, BCC ideal instável ao corte) e a topologia sólido(1)→sólido(2) é testada com titânio (segunda confirmação) e árgon (controlo negativo). O capítulo 8 situa o modelo entre onze conceções estabelecidas da termodinâmica estática. O apêndice documenta as fontes bibliográficas (ordenadas por autor) e o apoio da IA (proporção de inovação estimada, custo estimado e duração por passo de desenvolvimento); os auxílios reúnem o modo de leitura, os pontos críticos remanescentes e as perspetivas.

Enunciado central. Uma única medida cooperativa da ligação, z̄, agrafa três tradições de resto separadas — a imagem do gás em rede / de Ising (transição cooperativa de dois estados), a imagem da percolação (conectividade → ebulição) e a imagem dos fonões de Einstein/Debye (rigidez → fusão) — num quadro de fases coerente e qualitativamente completo, deliberadamente dentro de uma imagem ideal do tipo van der Waals, podendo naturalmente o edifício teórico ser deduzido de uma imagem dinâmica facilmente imaginável. Não é a exatidão mas o caráter figurado que constitui a verdadeira força do modelo e que caracteriza o seu valor científico enquanto modelo.

Índice

• Resumo geral
• 1 — Os axiomas com as suas palavras-chave (consolidados)
• 2 — Análise de contradições e lacunas
• 3 — Descrições figuradas
• 3.1 As três fases isoladas
• 3.2 Transições de fase
• 3.3 Porque a velocidade do som nos líquidos desce com T
• 3.4 Como surge o movimento browniano — com uma noção afinada de micela
• 3.5 Porque as reações nas camadas limite das micelas se tornam mais rápidas
• 3.6 A janela líquida tridimensional
• 4 — Modelo de mecânica estatística
• 4.1 Graus de liberdade e variáveis
• 4.2 Ligação cooperativa de dois estados (do tipo Ising, parâmetro J)
• 4.3 Balanço energético — com termo cooperativo
• 5 — Funções de partição e capacidades térmicas
• 5.1 Função de partição fatorizada
• 5.2 Capacidade térmica por fase (analítica)
• 5.3 Curva numérica c_v(τ) **com** termo cooperativo — calor latente quantificado
• 5.4 Limiar de percolação de conectividade como resultado (confirmado inalterado) - 5.4.a Dois limiares distintos nitidamente separados
• 5.5 Ordem da transição de fusão
• 5.6 Estudo de tamanho finito do limiar de percolação
• 5.7 Varrimento 3D completo de c_v(τ) da fusão de primeira ordem
• 6 — Plausibilidade matemática da interpretação 3D
• 6.1 Alteração do grau de liberdade de rotação 2D → 3D
• 6.2 Influência das diferentes redes cristalinas — a conectividade como traço distintivo
• 6.3 A transição adicional sólido(1) → sólido(2)
• 6.4 Considerações: a anomalia acústica da água líquida
• 6.5 Alternativa: campos de vórtices fractais em vez da rotação das micelas
• 6.6 Rigidez fonónica microscópica θ(rede): um modelo de molas de Einstein
• 6.7 Estimativa matemática da metastabilidade (2D triangulada ↔ quadrada, 16b)
• 6.8 Refinamento do tipo Debye de θ(rede)
• 6.9 θ(rede) a partir de uma rede de Lennard-Jones completa
• 7 — Calibração e equação de estado
• 7.1 Método de calibração (um parâmetro de energia ε) — ajustado
• 7.2 Tabela de calibração
• 7.3 Eixo densidade/volume: equação de estado de van der Waals com caroço duro
• 7.4 Calibração quantitativa l/g: a → ε, b → D²
• 7.5 Superfície (τ, v) comum com ponto triplo
• 7.6 Calor latente face às entalpias de fusão
• 7.7 Unificação das escalas de energia: ε ↔ ε_lg
• 7.8 Ponto triplo sem parâmetros a partir de z̄_c
• 7.9 Cálculo de controlo 3D: uma janela líquida limpa numa única escala ε
• 7.10 Varrimento real de percolação de ligações em 3D (FCC/BCC) e sensibilidade do bordo da janela
• 7.11 Varrimento real em τ da cinética de ligação em 3D: τ_b(3D) e a linha de fusão numa escala ε
• 7.12 A linha de transformação sólido(1)→sólido(2) na superfície (τ, v)
• 7.13 Linha de transformação sólido(1)→sólido(2) no diagrama (τ, p), calibrada com ferro
• 7.14 Redlich–Kwong como ramo separado da imagem ideal
• 7.15 sólido(1)→sólido(2) noutros polimorfos: titânio e controlo com gás nobre
• 8 — Comparação bibliográfica: conceções da termodinâmica estática
• 8.1 Tabela de conjunto
• 8.2 Os modelos um a um
• 8.3 Semelhanças e diferenças face ao modelo 3
• 9 — Apêndice
• 9.1 Bibliografia
• 9.2 Apoio da IA
• 9.3 Glossário de abreviaturas e símbolos
• 9.4 Scripts de cálculo e cobertura de testes

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

1 — Os axiomas com as suas palavras-chave (consolidados)

Resumo. Sete axiomas dotados de palavras-chave fixam a substância ideal: número de partículas conservado, movimento puramente clássico, uma ligação cooperativa de dois estados e o número de coordenação médio z̄ como parâmetro de ordem. Sustentam todas as derivações que se seguem; no texto mestre alemão uma marca de água de homóglifos assinala as palavras-chave (chave no fim do texto).

Os axiomas precisados:

1. VOO LIVRE (axioma adicional). Entre duas moléculas que passam uma pela outra sem estarem ligadas não atuam nem atrito nem forças eletrostáticas. Momento e energia só mudam pela repulsão dura do CAROÇO DURO no choque e pela formação e rutura de uma LIGAÇÃO. O axioma só importa do ponto de vista da mecânica clássica e do balanço energético: garante que não existem canais dissipativos ao nível molecular — a energia do movimento nunca se «perde», apenas passa aos três reservatórios da TRÍADE ou a energia de ligação. (Atrito e viscosidade macroscópicos nascem de forma emergente do quebra-cabeças deslizante e da FISSURA, não como força fundamental microscópica.)
2. MULTIDÃO. Uma substância ideal é feita de muitíssimas moléculas iguais; as propriedades macroscópicas são descritas estatisticamente no limite termodinâmico (N → ∞).
3. LIGAÇÃO. Duas moléculas geometricamente vizinhas podem portar uma ligação. Uma ligação intacta baixa a energia da quantidade finita ε; a sua rutura custa exatamente ε. A ligação é um sistema de dois estados {intacta, rompida}.
4. TRÍADE. A energia interna reparte-se classicamente por três reservatórios: (a) translação das moléculas, (b) vibração nas ligações intactas (uma mola harmónica por ligação), (c) rotação de blocos moleculares rígidos (grau de liberdade rotacional coletivo das regiões conexas).
5. CAROÇO DURO. Se a distância molecular descer abaixo de um valor mínimo σ, atua uma repulsão dura (volume próprio). Juntamente com a LIGAÇÃO resulta um potencial de pares efetivo (casca dura mais poço de ligação finito).
6. CONSERVAÇÃO. O momento total e a energia total conservam-se. Uma molécula isolada é energeticamente neutra — não armazena energia em si; a energia vive nas ligações (LIGAÇÃO) e nos graus de liberdade cinéticos (TRÍADE).
7. FISSURA. As ligações desfazem-se e formam-se segundo o princípio fissuração/cicatrização: quantas mais ligações vizinhas faltarem → maior a probabilidade de rutura; quantos mais vizinhos ligados → maior a probabilidade de cicatrização. A cinética é cooperativa.
8. LESTO. Rutura e cicatrização decorrem depressa e com pouca ativação; o subsistema das ligações encontra-se a toda a temperatura em equilíbrio local e pode ser separado adiabaticamente do movimento molecular (fundamento da função de partição fatorizada).

Postura de modelação. A coordenação geométrica conserva-se de forma idealizada; só muda o número de coordenação z pelas ruturas. O movimento surge como num quebra-cabeças deslizante: as cavidades (do volume próprio do CAROÇO DURO e das ruturas da FISSURA) migram, os blocos avançam. A imagem deliberadamente «errada» de cristalitos que rodam e migram mantém-se como base intuitiva.

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

2 — Análise de contradições e lacunas

Resumo. A análise examina sistematicamente o modelo quanto a contradições internas e lacunas e dá continuidade ao estado do trabalho; as fraquezas reconhecidas são nomeadas e resolvidas nos capítulos seguintes ou assinaladas como limite honesto do modelo.

Lista de pontos críticos

1. Dupla contagem da translação (axioma 3a contra 3c). O texto em bruto nomeia a translação duas vezes (translação molecular e translação de bloco). Dois reservatórios de translação independentes não existiriam. Resolução: a translação de bloco é apenas um modo coletivo das translações moleculares já contadas; genuinamente novo é somente o grau de liberdade de rotação do bloco. No balanço dos graus de liberdade (cap. 4) a translação é contada exatamente uma vez.
2. «Energeticamente neutra» (axioma 5) contra «as ligações absorvem energia» (axioma 2/3b). Contradição aparente. Resolução: a neutralidade diz respeito à molécula enquanto corpo; a energia reside entre as moléculas (na LIGAÇÃO) e em modos cinéticos coletivos, nunca num amortecedor intramolecular. Molécula ≠ ligação — assim a contradição fica resolvida.
3. Atração ausente (o axioma 4 só nomeia a repulsão). Sem atração não haveria fase condensada. Resolução: a atração é fornecida pela LIGAÇÃO (poço de profundidade ε), a repulsão pelo CAROÇO DURO (casca dura em σ). Juntas, um potencial de pares efetivo.
4. Retroação sem termo oposto (axioma 6). «Mais vizinhos em falta → mais ruturas» é por si só uma retroação positiva (avanço da fissura até à desagregação total). Resolução: o termo de cicatrização (mais vizinhos → mais cicatrização) mais o fator de Boltzmann da energia de ligação finita ε (LIGAÇÃO) equilibram a fissura e produzem uma transição definida em vez de uma descontrolada.
5. Separação de escalas temporais não explicitada (axioma 7). «Rápida e de baixa ativação» só se torna utilizável quando daí se extrai a separação adiabática entre dinâmica das ligações e dinâmica de translação. Resolução: consignada explicitamente como pressuposto da fatorização Z = Z_trans·Z_vib·Z_rot·Z_bond.
6. Geometria ausente e mecanismo do movimento ausente. Os axiomas não nomeiam nem a rede de coordenação nem o que torna possível o escoamento. Resolução (como suposição adicional explicitamente assinalada): uma geometria de coordenação de referência fixa com ocupação variável das ligações; como caso 2D canónico a coordenação triangular (z_max = 6). Mobilidade através do mecanismo de cavidades do quebra-cabeças deslizante. Estes dois pontos são as únicas adições genuínas para além dos sete axiomas e estão assinalados no capítulo 6 como questão estratégica.
7. Lacuna latente de dissipação. A versão 07 deixava em aberto se moléculas não ligadas perdem energia ao passar. Resolvida pelo VOO LIVRE: não há atrito nem eletrostática microscópicos; todos os fluxos de energia e momento correm exclusivamente pelos choques do CAROÇO DURO e pela cinética da LIGAÇÃO. Assim o balanço energético (cap. 4.3) fica fechado e a estatística de Maxwell–Boltzmann exata no limite gasoso (gás ideal livre entre choques elásticos).
8. Incoerência energética do parâmetro de ordem. A versão 08 usava o termo cooperativo −(K/2)x² na energia livre g(x) (para gerar x(τ)), mas não na energia interna U/N. Isso é termodinamicamente incoerente: o mesmo termo de campo molecular tem de aparecer em ambas. Resolução (cap. 4.3): −(K/2)x² é transposto para U/N; o calor latente aumenta com isso de forma mensurável (cap. 5.3), sem alterar os patamares sólido e gasoso.
9. Calibração l/g apenas qualitativa. A versão 09 tinha o eixo das densidades (equação de estado de van der Waals) apenas em forma reduzida. Resolução (cap. 7.4): a linha de ebulição absoluta é calculada através do ponto crítico (T_c, p_c) e T_b(p_atm) é confrontado com valores medidos; o desvio sistemático (regra de Guldberg) é nomeado.
10. Eixos separados sem ponto triplo. Os eixos de fusão e de ebulição eram cálculos separados. Resolução (cap. 7.5): ambos são reunidos num plano (τ, v), o ponto triplo é construído e o ponto crítico assinalado.

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

3 — Descrições figuradas

Resumo. O capítulo 3 explica de forma intuitiva o sólido, o líquido e o gasoso bem como as suas transições (com tiras estáticas): as três fases isoladas (3.1), as transições de fase recristalização/fusão/ebulição/sublimação juntamente com o ponto triplo e o crítico (3.2), a anomalia acústica (3.3), o movimento browniano e a noção afinada de micela (3.4–3.5), bem como a janela líquida tridimensional (3.6). A partir destas descrições foram elaborados esquemas (tiras) destinados a ajudar a representar diante do olho interior as relações moleculares dinâmicas nas diversas fases. Essa capacidade de representação é o verdadeiro núcleo do modelo. Ele ganha valor qualitativo e científico quando a mecânica estatística permite mostrar de forma plausível que o modelo torna previsíveis tendências e observações empíricas qualitativas mesmo sem grande aparato matemático.

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

3.1 As três fases isoladas

Sólido. Quase todos os vizinhos estão ligados entre si por ligações moleculares (LIGAÇÃO); a média dos números de coordenação z das moléculas está perto do máximo geométrico. A rede molecular é um único bloco rígido sem cavidades (CAROÇO DURO, empacotamento denso) — um quebra-cabeças deslizante interno como no líquido não é possível na versão compacta. Único movimento: o tremor em torno das posições de repouso (vibração de ligação, TRÍADE-b), determinado pelas ligações moleculares com os vizinhos. Momento e energia oscilam elasticamente na rede (CONSERVAÇÃO, VOO LIVRE: sem perdas). No caso limite (temperatura muito baixa, empacotamento ideal) tem-se o estado ideal de um cristal perfeito.

Líquido. A FISSURA cooperativa atravessou a rede: faltam ligações às manchas, a média dos números de coordenação z das moléculas desceu, o bloco rígido desfez-se interiormente em muitos blocos menores. Do monólito nasceu um monte de escombros. Estes blocos transladam e rodam (TRÍADE-c). Entre eles abriram-se cavidades pelas quais se deslocam uns ao lado dos outros — o quebra-cabeças deslizante do movimento browniano está ativo. A geometria de coordenação local conserva-se de forma idealizada; só mudou quantas ligações existem realmente e se formam.

Gás. A FISSURA venceu: as ligações estão rompidas quase por toda a parte (z → 0). No caso extremo nenhuma molécula tem já ligação molecular com outra. Os blocos desagregaram-se em moléculas isoladas. Voam livres, acopladas apenas por choques do CAROÇO DURO (com VOO LIVRE pelo meio); as velocidades seguem Maxwell–Boltzmann (CONSERVAÇÃO). Vibração de ligação e rotação de bloco extinguiram-se — translação pura (TRÍADE-a). No caso limite (temperatura alta, pressão baixa) tem-se o estado ideal de um gás ideal.

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

3.2 Transições de fase

Recristalização (sólido(1)/sólido(2)). Ainda abaixo da fusão, o sólido pode mudar a sua rede sem perder a rigidez: dois empacotamentos ideais com números de coordenação diferentes (FCC z_max = 12 ↔ BCC z_max = 8) têm entalpias livres μ(τ, p) distintas; em μ_FCC = μ_BCC a geometria de coordenação salta abruptamente (não apenas a ocupação z̄ como na fusão). Torna-se visível um pequeno degrau de c_v (Δc_v = variação de z_max/4 = 1,0 k_B) mais um fraco calor latente de transformação. Afirmação da IA: só é concebível em 3D (conectividade discreta, quantizada; cap. 6.3). No caso bidimensional a verdadeira recristalização, isto é, a transição entre duas fases sólidas, é mais fácil de imaginar. No plano só existem três pavimentações regulares sem buracos. Pode cobrir-se inteiramente uma superfície (pavimentá-la) com triângulos equiláteros, com quadrados ou com hexágonos regulares. Se nestas redes as arestas representarem as ligações moleculares e os vértices as moléculas, então é concebível uma transição de fase sólido(triangular) → sólido(quadrado). A molécula na pavimentação triangular teria ainda seis ligações às moléculas vizinhas, ao passo que a da rede quadrada só teria quatro.

Fusão (s/l). Ao subir a temperatura a rede rasga-se em pontos suficientes para que o bloco rígido total perca a sua rigidez (a FISSURA cooperativa desce abaixo do limiar de rigidez de Maxwell z̄ = 2d — a rede torna-se mole). Formam-se micelas termodinâmicas, isto é, regiões que já não estão unidas por ligações às regiões vizinhas maiores. Do ponto de vista molecular o monólito torna-se um monte de escombros e adquire uma capacidade de escoamento limitada. As primeiras cavidades permitem o deslocamento de blocos moleculares, como no quebra-cabeças deslizante. As ligações em falta para a vizinhança permitem a rotação interna nas micelas. A rotação das micelas (TRÍADE-c) inicia-se. A rutura de muitas ligações (LIGAÇÃO) devora calor latente.

Ebulição (l/g). Um aquecimento adicional faz percolar as próprias cavidades: unem-se num espaço vazio contínuo, no qual os blocos se dissolvem e moléculas isoladas voam para longe. A ebulição é a passagem de «cavidades na matéria» a «matéria numa cavidade».

Sublimação (s/g). A baixa pressão falta a zona líquida intermédia: a rede cede moléculas isoladas diretamente ao espaço vazio, sem que ocorra a formação de blocos como nas micelas. Pode perguntar-se aqui que influência podem ter na formação ou não formação de micelas termodinâmicas os aspetos geométricos (moléculas longas, moléculas que se entrelaçam), os aspetos eletrostáticos (alinhamento por momentos magnéticos ou dipolares) e os diferentes tipos de ligação molecular (van der Waals, interações entre sistemas mesómeros, pontes de hidrogénio).

Ponto triplo. Precisamente aí coincidem dois limiares: a rigidez de Maxwell (sólido ↔ líquido, perda da rigidez) e a percolação das cavidades (líquido ↔ gás) ocorrem no mesmo (τ, p). As três imagens (bloco rígido [monólito], blocos que deslizam [monte de escombros], moléculas livres [nuvem de poeira]) coexistem.

Ponto crítico. Ao subir a pressão o salto líquido→gás torna-se cada vez menor; a fração de vazio muda de forma cada vez mais contínua. No ponto crítico a diferença desaparece por completo — «blocos densos com pouco vazio» e «muito vazio com poucas moléculas» passam continuamente um ao outro. («Do ponto de vista do líquido o copo está meio vazio; do do gás, meio cheio.»)

Recristalização na imagem 2D. A tira mostra agora a recristalização como caso 2D, um por baixo do outro: em baixo a rede densa triangulada (coordenação z = 6, o ideal 2D de menor energia), em cima a rede quadrada frouxa (z = 4). A transição triangulada → quadrada aumenta o volume (trabalho contra a pressão) e reduz o número de ligações. Vê-se que o número de coordenação muda enquanto a rigidez de uma rede cristalina permanece.

Previsão qualitativa de transições de energia e capacidades térmicas. Só a partir da conectividade z se podem estimar ambos os valores característicos: (i) a energia de ligação por partícula é U/N ≈ −(ε/2)·z (cada ligação conta metade), logo tanto menor quanto maior for z; (ii) o salto de c_v na mudança de rede segue a regra do modelo Δc_v = variação de z_max/4. Para o caso 2D triangulada (z = 6) ↔ quadrada (z = 4) resulta uma diferença de energia de ΔU ≈ ε por partícula e um degrau de c_v de Δc_v = (6 − 4)/4 = 0,5 k_B — metade do valor 3D (1,0 k_B). Regra mnemónica: mais denso = menor energia, mas mais caro na transição (maior calor latente, maior degrau de c_v), porque ao alargar-se se libertam mais ligações e mais graus de liberdade congelados.

Determina a estrutura da rede a metastabilidade? Em parte. A estrutura da rede fixa a ordem das energias — que forma é o estado fundamental (triangulada, μ mais baixo) e qual a excitada (quadrada). Se a forma frouxa sobrevive como metastável a baixa temperatura, não o decide sozinha: no modelo estritamente ideal (sem barreira) a substância adota em cada (τ, p) imediatamente a fase de μ mais baixo — aí não existe metastabilidade. Esta só surge por uma energia de ativação finita (barreira de nucleação) que trava cineticamente a reorganização. A geometria da rede fornece portanto a hierarquia termodinâmica, a barreira a persistência cinética — o mesmo ingrediente que explica também a fase supersaturada (ver abaixo; quantitativamente cap. 6.7).

Fase supersaturada — energia de ativação. A figura adicional (3.2, fila 2) mostra uma fase gasosa supersaturada sobre um líquido: a par de moléculas isoladas, pequenos grupos de 2–3 partículas agregadas (barras cinzas) que voam como um todo. O modelo explica a supersaturação por uma energia de ativação apreciável para as ligações moleculares: os pequenos agregados (2–3) têm uma relação borda-superfície desfavorável, as suas ligações de borda são incompletas — só um núcleo crítico cresce espontaneamente. Enquanto não surgir um núcleo supercrítico, o gás permanece supersaturado, embora μ_gás > μ_líquido já permitisse a condensação. É precisamente essa barreira que falta na imagem ideal de equilíbrio e é aqui nomeada como contributo autónomo situado «fora da imagem ideal» (estimativa quantitativa no cap. 6.7).

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

3.3 Porque a velocidade do som nos líquidos desce com T

O som é uma onda elástica: v = √(K_eff/ρ). K_eff nasce exclusivamente das ligações intactas (LIGAÇÃO) — são elas as molas. Ao subir T rompem-se cooperativamente cada vez mais ligações (FISSURA); a rede de molas afina e K_eff cai mais abruptamente do que a densidade ρ. A rutura domina o numerador, não o denominador — por isso v desce. Interpretação molecular: a transmissão do som por ligações moleculares é mais rápida do que a transmissão mecânica por choques. Quantas mais fissuras surgem num líquido ao aumentar a temperatura, tanto mais vezes o som tem de ser transportado por choques mecânicos em vez de ligações moleculares. Em consequência, a velocidade do som torna-se mais lenta ao subir a temperatura. (A anomalia da água — um máximo acústico por volta dos 74 °C — é tratada no cap. 6.4.)

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

3.4 Como surge o movimento browniano — com uma noção afinada de micela

Em sentido termodinâmico-estatístico uma micela não é aqui um agregado de tensioativos mas um bloco cercado pela FISSURA: as ligações em falta (as não-ligações na fronteira do bloco) formam a «casca» em torno de um cristalito encerrado (o núcleo integralmente ligado). Esse núcleo pode executar um movimento de rotação cooperativo interno (TRÍADE-c) — o bloco roda como um todo enquanto a sua casca de fronteira o desacopla dos vizinhos.

O movimento browniano surge então assim: mediante a troca de ligações moleculares (FISSURA, LESTO) moléculas isoladas podem deixar a casca de uma micela e passar a outra micela. Cada uma destas mudanças transmite momento (CONSERVAÇÃO; sem perdas graças ao VOO LIVRE) e altera abruptamente a forma e o momento angular de ambas as micelas envolvidas. Uma grande partícula de teste (grão de poeira) é banhada por muitas destas reorganizações micelares não correlacionadas; a soma dos impulsos é nula em média, mas tem variância finita → caminho aleatório. Quanto mais alta for T, tanto mais frequentes são os eventos de FISSURA e as mudanças de micela, e tanto mais forte o tremor. O movimento browniano é assim a sombra macroscópica, ou consequência, do processo microscópico de troca entre micelas.

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

3.5 Porque as reações nas camadas limite das micelas se tornam mais rápidas

Observação preliminar embaraçosa. A argumentação que se segue é uma dedução a partir da imagem e, tanto quanto sei, até agora não foi demonstrada nem investigada empiricamente — embora deva confessar que não me dei ao trabalho de pesquisar o assunto a fundo. O único indício de novidade é que os sistemas de IA não mencionam esta conclusão. Esta tese é um exemplo de como representações figuradas podem conduzir a novas perguntas.

Em solventes de dissolução moderada os reagentes fixam-se de preferência às cascas de fronteira das micelas — no modelo, precisamente as fronteiras fissura/bloco (FISSURA). Dois efeitos aceleram a reação: (i) a redução dimensional — os reagentes já não se procuram no volume inteiro mas ao longo de uma interface quase bidimensional; uma busca 2D acerta estatisticamente mais depressa do que uma 3D. (ii) A concentração — a casca de fronteira reúne os reagentes, a sua densidade local aumenta. Em imagem, as fronteiras de fissura atuam como carris guia que conduzem os parceiros um ao outro.

O efeito dimensional tem um núcleo rigoroso: o passeio aleatório (Pólya). Que uma busca bidimensional acerte estatisticamente melhor do que uma tridimensional não decorre apenas da intuição, mas de um teorema do cálculo de probabilidades. George Pólya (1887–1985) mostrou em 1921 para o passeio aleatório numa rede: em uma e em duas dimensões o caminhante regressa com probabilidade 1 (ou seja, com certeza) infinitas vezes a cada ponto, tanto na rede plana ℤ² como na reta ℤ. A partir de três dimensões isso já não vale: a probabilidade de regresso é estritamente menor do que 1 (na rede cúbica simples ℤ³ apenas cerca de 0,34). Em imagem, segundo Pólya: «Um homem bêbado acaba de certeza por encontrar o caminho de casa; um peixe bêbado, esse não.» Transposto para dois reagentes que têm de se encontrar várias vezes para um choque bem sucedido, isto significa: num solvente muito bom com passeio aleatório tridimensional livre, parceiros uma vez separados afastam-se para sempre com probabilidade finita — as tentativas repetidas de reação falham e o rendimento fica baixo por esta razão puramente geométrico-estatística. Numa camada limite micelar (quase) bidimensional, pelo contrário, o reencontro é certo (probabilidade 1): os parceiros voltam a encontrar-se vezes sem conta até a reação ter êxito. A camada limite não só concentra como põe ao mesmo tempo à disposição a dimensão mais favorável aos reencontros. A afirmação (deliberadamente aguçada) reza: os rendimentos em solventes muito bons tendem a ser piores por razões de passeio aleatório.

Referência: G. Pólya (1921), Über eine Aufgabe der Wahrscheinlichkeitsrechnung betreffend die Irrfahrt im Straßennetz, Math. Ann. 84, 149–160 — citado aqui segundo M. Eisermann, Höhere Mathematik 3 (Univ. de Estugarda), exercício W-2×2. (Original consultado a 5.9.2026).

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

3.6 A janela líquida tridimensional

Ao contrário das figuras de 3.1–3.5, esta não mostra nenhuma cena molecular mas o resultado do cálculo dos cap. 7.9/7.10 como esquema num eixo τ (em unidades da energia de ligação ε). O objetivo era tornar visível porque a rede 2D isolada dificilmente admite um domínio para uma fase líquida, ao passo que com a rede 3D isso é possível sem dificuldade, valendo em ambos os casos a mesma regra de rigidez (z̄_c = 2d).

Imagem ① — 2D (sem janela). No eixo τ, τ_m = 1,20 e τ_triple = 1,16 sobrepõem-se praticamente (τ_triple abaixo de τ_m), com τ_c = 1,50 logo acima: o domínio da fase líquida está comprimido numa fenda — o sólido passa quase diretamente a gás (perto da sublimação).

Imagem ② — 3D (janela limpa). O mesmo eixo, mas dilatado: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00. Entre τ_triple e τ_c está o campo líquido autónomo (largo para baixo, τ_triple − τ_m = 1,63; estreito para cima, τ_c − τ_triple = 0,18). O fino bordo superior está assinalado como tal.

Imagem ③ — a causa. z̄_c = 2d cresce com a dimensão (4 → 6), pelo que v_solid desce (0,50 → 0,667 em volumes reduzidos, isto é, o sólido torna-se «mais frouxo» relativamente ao líquido), e τ_c = z_max/4 duplica (1,50 → 3,00). O ponto crítico afasta-se muito acima do ponto de fusão. Só assim surge espaço para a fase líquida. É um puro efeito dimensional.

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

4 — Modelo de mecânica estatística (inalterado desde a versão 09)

Resumo

O capítulo 4 estabelece a armação de mecânica estatística do modelo em 2D (rede triangular, z_max = 6):

• 4.1 Graus de liberdade e variáveis: a grandeza de ordem é a fração de ligações x ∈ [0,1] com o número de coordenação médio z̄ = z_max·x; por bloco ≥ 2 acresce um grau de liberdade de rotação, e a temperatura τ conta-se em unidades da energia de ligação ε.
• 4.2 Ligação cooperativa de dois estados: a equação de ligação autoconsistente x = 1/(1 + G·exp[−(ε + K·x)/τ]) (K = J·z_max, G = z_max) segue como condição de estacionaridade da energia livre de Bragg–Williams; o seu mínimo global fornece o ramo físico, e o cap. 5.5 decide se esse ramo decorre continuamente ou com um salto (primeira ordem).
• 4.3 Balanço energético com termo cooperativo: U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x² (2D) reúne translação, vibração de ligação e rotação de bloco, a energia de configuração das ligações e a nova energia cooperativa de campo molecular −(K/2)x². A capacidade térmica c_v associada recebe o termo adicional positivo −K·x·ẋ, que agudiza o pico de fusão mas deixa intactos o patamar sólido e o ramo gasoso.

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

4.1 Graus de liberdade e variáveis

• N moléculas sobre uma geometria de coordenação de referência fixa. Caso 2D canónico: rede triangular, z_max = 6, um GDL de rotação por bloco (MULTIDÃO, LIGAÇÃO). Este capítulo vale explicitamente para a interpretação 2D; a generalização 3D está no cap. 6.
• Variável de ligação s_ij ∈ {0, 1} por par geométrico de vizinhos (1 = intacta).
• Número de coordenação z_i = Σ_j s_ij; média z̄ = 2 N_bond / N — o parâmetro de ordem. Tem-se z̄ = z_max · x com a fração de ligações x = p_b ∈ [0, 1].
• Unidade de energia τ ≡ k_B T (classicamente, ½τ por grau de liberdade quadrático). Todas as energias em unidades da energia de ligação ε.

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

4.2 Ligação cooperativa de dois estados (do tipo Ising, parâmetro J)

Uma ligação isolada (energia −ε intacta, 0 rompida) tem, no equilíbrio de Maxwell–Boltzmann, p_b⁰(τ) = 1/(1 + e^(−ε/τ)). A FISSURA cooperativa (axioma 6) desloca a estabilidade efetiva da ligação de um campo molecular proporcional à fração de ligações x. Em aproximação de campo molecular (Ising) a equação de ligação autoconsistente torna-se:

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

Aqui:

• K = J·z_max é a força de acoplamento cooperativa (retroação da FISSURA: mais ligações vizinhas → maior estabilidade efetiva da ligação). J > 0 agudiza a transição; acima de um J crítico z̄(τ) torna-se descontínuo — esse é o calor latente da fusão. O caso puro de dois estados é o limite J = 0 (sem campo molecular, crossover contínuo).
• G > 1 é a degenerescência do estado rompido: uma molécula libertada ganha espaço de fases (liberdade de translação). G impulsiona a rutura de ligações a τ elevado e faz com que x → 1/(1+G) em vez de ficar preso em ½ (caso limite do gás). Canonicamente G = z_max.

Esta equação segue como condição de estacionaridade de uma energia livre de Bragg–Williams por molécula (n_slots = z_max/2 lugares de ligação):

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

O termo −(K/2)x² é o acoplamento cooperativo (ferromagnético), o último termo a degenerescência de rutura. O mínimo global de g sobre x fornece o ramo físico; onde dois mínimos se tornam igualmente profundos situa-se a transição de primeira ordem. LESTO (axioma 7) justifica que x(τ) tome a cada temperatura o valor de equilíbrio.

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

4.3 Balanço energético — com termo cooperativo

O termo de campo molecular −(K/2)x² é uma energia genuína (independente de τ) e tem de aparecer em U/N. O balanço coerente reza portanto (2D):

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

• τ·[…]: translação (1) + vibração de ligação (z̄/4) + rotação de bloco (½ n_blk) — TRÍADE.
• −(ε/2)·z̄: energia de configuração das ligações — LIGAÇÃO.
• −(K/2)·x²: energia cooperativa de campo molecular — FISSURA. Novo em U/N.

A capacidade térmica associada c_v = d(U/N)/dτ (em unidades k_B) recebe um contributo adicional:

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

onde o novo termo −K·x·ẋ é positivo (ẋ < 0) e atua exatamente onde muitas ligações se rompem em simultâneo. Ele agudiza o pico de fusão sem alterar o patamar sólido (ẋ ≈ 0) nem o ramo gasoso.

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

5 — Funções de partição e capacidades térmicas

Resumo

O capítulo 5 calcula funções de partição e capacidades térmicas e fornece os dados do modelo para a calibração (cap. 7):

• 5.1 Função de partição fatorizada: pela separação das escalas temporais, Z = Z_trans·Z_vib·Z_rot·Z_bond fatoriza-se.
• 5.2 Capacidade térmica por fase: c_v = d(U/N)/dτ; o termo positivo −(ε/2)(dz̄/dτ) gera o pico do tipo Schottky como análogo do calor latente.
• 5.3 Curva numérica c_v(τ): pico de c_v 5,51 em τ_m = 1,20, patamar sólido 2,50, ramo gasoso → 1; o calor latente L ≈ 1,41 ε por molécula (com termo cooperativo) fica triplicado face a 0,48 ε (sem ele).
• 5.4 / 5.4a Dois limiares separados: a percolação de conectividade z̄_c = z_max·p_c ≈ 2,08 (em τ_b ≈ 1,85) assinala a desagregação dos aglomerados (perto da ebulição) e é um resultado; deve distinguir-se nitidamente da rigidez genérica de Maxwell z̄_c = 2d (densidade de fusão, perda de rigidez).
• 5.5 Ordem da transição de fusão: critério ln G > 2 + 4ε/K; o 2D (G = 6) é contínuo para todo J, o 3D (G = 12) é de primeira ordem a partir de K > ~8,2 (demonstração: Δz̄ ≈ 5,3, ΔU ≈ 1,75 ε).
• 5.6 Estudo de tamanho finito: p*(L) converge monotonamente por baixo para p_c = 0,3473; z̄_c ≈ 2,08 é robusto face ao tamanho finito.
• 5.7 Varrimento 3D completo de c_v: a partir de K > K_crit ≈ 8,25 (J > 0,68) abre-se o verdadeiro salto Δz̄ com calor latente ΔU; o pico de fusão desloca-se com J para τ maiores (deslocamento positivo de +0,5 … +2,2 face ao pico 2D contínuo).

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

5.1 Função de partição fatorizada

Pela separação das escalas temporais (LESTO) a função de partição canónica fatoriza-se:

Z = Z_trans · Z_vib · Z_rot · Z_bond.

• Z_trans = (1/N!)·(A/λ_T²)^N, λ_T = h/√(2π m τ) (translação 2D clássica ideal; MULTIDÃO, CONSERVAÇÃO, VOO LIVRE).
• Z_vib = (τ/ħω₀)^(N_bond) (oscilador harmónico clássico por ligação intacta; TRÍADE-b).
• Z_rot = ∏_(blocos≥2) (2π I_b τ)^(1/2)/ħ (rotor rígido 2D clássico por bloco; TRÍADE-c).
• Z_bond = Σ_{s} exp[(ε/τ)·Σ s_ij + (K/τ)·termo cooperativo] (soma de ligações do tipo Ising; LIGAÇÃO, FISSURA). A avaliação em campo molecular fornece x(τ), z̄(τ) de 4.2.

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

5.2 Capacidade térmica por fase (analítica)

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

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

O termo −(ε/2)·(dz̄/dτ) é positivo (dz̄/dτ < 0: aquecer rompe ligações) e domina onde muitas ligações se rompem em simultâneo — o pico do tipo Schottky, análogo do calor latente.

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

5.3 Curva numérica c_v(τ) com termo cooperativo — calor latente quantificado

Ligação autoconsistente (mínimo global), n_blk(τ) de um Monte Carlo de percolação de ligações (union-find, L = 100, 5 realizações), c_v = dU/dτ (diferenças centradas), U/N incluindo −(K/2)x² (cv_scan_09.py, ε = 1, J = 0,55, G = z_max = 6). Pontos de apoio escolhidos:

Resultado (numérico). A ordem das três fases mantém-se confirmada e é agora mais nítida:

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

Calor latente (área do pico acima da linha de base). Integrando (c_v − linha de base) sobre a janela de fusão τ ∈ [0,70 ; 1,70] (linha de base = reta por c_v(0,70) e c_v(1,70)) obtém-se

L ≈ 1,41 ε por molécula (com termo cooperativo), L ≈ 0,48 ε por molécula (sem termo cooperativo, balanço antigo), ΔL ≈ +0,92 ε por molécula.

O termo cooperativo eleva o pico num fator ≈ 1,7 (3,25 → 5,51) e triplica o calor latente. Fisicamente: a rutura simultânea de muitas ligações vizinhas (FISSURA) armazena não só energia de configuração (−ε/2·z̄) mas também energia de campo molecular (−K/2·x²) — ambas se libertam ao fundir.

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

5.4 Limiar de percolação de conectividade como resultado

A fração do aglomerado gigante P_∞(τ) colapsa em τ ≈ 1,85; aí x ≈ 0,347 = p_c(triângulo) = 2·sin(π/18) = 0,3473, logo z̄_c = z_max·p_c ≈ 2,08 — um resultado, não algo imposto. Este é o limiar de conectividade (aglomerado que atravessa) e assinala a desagregação dos aglomerados perto da ebulição (τ_b); não deve confundir-se com a rigidez de Maxwell (z̄_c = 2d, fusão) — cf. a separação limpa no cap. 5.4a.

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

5.4.a Dois limiares distintos nitidamente separados

A dupla designação até agora difusa «percolação/rigidez» fica agora resolvida. Trata-se de dois limiares fisicamente distintos que não devem ser confundidos:

Fundamentação. A fusão é a perda da rigidez mecânica de um bloco denso, não a desagregação da conexidade. A grandeza relevante é portanto a contagem de restrições de Maxwell (campo médio / rigidez genérica): uma rede D-dimensional de forças centrais torna-se isostática exatamente quando a coordenação média atinge z̄ = 2D (2D → 4, 3D → 6). A conectividade das ligações (z̄_c = 2,08) fica nitidamente abaixo; uma rede apenas conexa em z̄ = 2,08 é há muito mole (modos moles), logo já líquida. Por isso a conectividade assinala a desagregação dos aglomerados mais tardia, perto da ebulição, ao passo que a rigidez fixa a solidificação/liquefação mais precoce (densidade de fusão). Esta separação fixa o ponto triplo do cap. 7.8 sem parâmetros e de forma realista; não altera o limiar de conectividade 2,08 confirmado em 5.4/5.6, apenas a sua leitura (τ_b = perto da ebulição, não a fusão).

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

5.5 Ordem da transição de fusão

Critério analítico. A região espinodal (dois mínimos de g) existe para τ < K/4; a transição situa-se aproximadamente em τ_t ≈ (ε + K/2)/ln G. Um verdadeiro salto (transição de primeira ordem) exige τ_t < K/4, ou seja, no limite K → ∞:

ln G > 2 + 4ε/K.

Daí segue um resultado notável sobre a degenerescência de rutura G = z_max:

Confirmação numérica (máx. |Δx| por passo em τ, detect_first_order):

Leitura. No modelo 2D canónico (G = 6) a fusão é por princípio contínua (uma transição «suave» / crossover) — para nenhum J surge um verdadeiro salto, porque ln 6 < 2. Só a maior coordenação/degenerescência do empacotamento 3D (G = 12) permite, a partir de uma cooperatividade suficiente (K > 8,2), uma transição de primeira ordem com calor latente finito. Isso concorda com a realidade: a fusão 2D é frequentemente contínua (do tipo KTHNY), a fusão 3D é de primeira ordem. J = 0 (sem campo molecular) é sempre o caso de referência contínuo. — Fisicamente J designa a relação entre energia cooperativa e energia de ligação (K/ε = J·z_max/ε): quanto uma ligação vizinha intacta estabiliza a própria.

Demonstração de primeira ordem (motivada pelo 3D, cv_scan_09.py, J = 1,0, G = 12, z_max = 12): z̄ salta em Δz̄ ≈ 5,30 em τ ≈ 2,82; o calor latente como salto de energia ΔU ≈ 1,75 ε por molécula. O mesmo K numa rede G = 6 (2D) mantém-se contínuo — o caráter da ordem depende unicamente de G = z_max.

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

5.6 Estudo de tamanho finito do limiar de percolação

Percolação de ligações em redes triangulares L × L (periódicas), fração de aglomerado gigante P_∞(p) escalada sobre L (cv_scan_09.py; L = 50/100/200/400 com 8/5/3/2 realizações). Limiar exato p_c = 2·sin(π/18) = 0,34730 → z̄_c = z_max·p_c = 2,0838.

P_∞(p) (média ± dispersão):

Estimativa do limiar p*(L) (P_∞ atravessa 0,5):

Resultado. p*(L) converge monotonamente por baixo para p_c = 0,3473; a dispersão em p_c cai de ±0,103 (L = 50) para ±0,017 (L = 400) — a transição agudiza-se com L mas não desloca o limiar (valor de percolação universal). Assim z̄_c ≈ 2,08 = z_max·p_c fica confirmado como valor do modelo robusto face ao tamanho finito.

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

5.7 Varrimento 3D completo de c_v(τ) da fusão de primeira ordem

A série demonstrativa única de primeira ordem do cap. 5.5 (J = 1,0, G = 12) é agora ampliada para um varrimento completo de c_v(τ) sobre a cooperatividade J (firstorder_3d_16.py, conectividade FCC para n_blk, balanço energético 3D U/N = τ·[3/2 + z̄/4 + (3/2)n_blk] − (ε/2)z̄ − (K/2)x²). O critério analítico de ordem do cap. 5.5 (primeira ordem ⇔ K > K_crit = 4ε/(ln 12 − 2) = 8,25) separa os casos; o pico de referência 2D contínuo situa-se em τ_2D = 1,20 (c_v ≈ 5,5).

*caso contínuo: Δz̄/ΔU são o crossover abrupto (não um verdadeiro salto); τ_transição é aí o pico de c_v.

Resultado. (i) O caso 3D contínuo J = 0,55 dá τ_transição = τ_m ≈ 1,72 — a verificação cruzada com a cinética de ligação independente (bondkinetics_3d_15.py: τ_m(FCC) = 1,70) concorda. (ii) A partir de K > K_crit ≈ 8,25 (J > 0,68) abre-se o verdadeiro salto Δz̄ com calor latente finito ΔU (salto de energia através da descontinuidade) — a fusão torna-se de primeira ordem, exatamente como exige o critério ln 12 > 2. (iii) A verificação cruzada J = 1,0 reproduz o valor demonstrativo Δz̄ ≈ 5,4 do cap. 5.5. (iv) A altura do salto Δz̄ e a posição τ_transição crescem com J: a transição desloca-se para cima — a maior coordenação do empacotamento 3D e a cooperatividade mais forte elevam a temperatura de fusão. Face ao pico 2D contínuo (τ_2D = 1,20) obtém-se portanto por toda a parte um deslocamento positivo do pico (+0,5 … +2,2) que cresce monotonamente com K. O caráter de primeira ordem da fusão 3D fica assim documentado não apenas num único ponto, mas em toda a gama de cooperatividade.

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

6 — Plausibilidade matemática da interpretação 3D

Resumo

O capítulo 6 torna plausível a generalização da imagem de ligação de 2D para 3D de forma puramente matemática (sem visualização HTML própria):

• 6.1 Grau de liberdade de rotação 2D → 3D: a rotação de bloco passa de um eixo a três, a translação de duas direções a três, z_max de 6 a 12; os valores base de c_v sobem (c_v,sólido 2,5 → 4,5, c_v,gás 1 → 3/2), a forma da curva c_v(τ) mantém-se.
• 6.2 A conectividade como traço distintivo: z_max é uma grandeza quantizada e separa abruptamente o empacotamento mais denso (FCC/HCP, z_max = 12) do BCC (z_max = 8), sem estado intermédio.
• 6.3 Transição adicional sólido(1)→sólido(2): da discretude segue o polimorfismo FCC↔BCC com um degrau de c_v sem parâmetros de 1,0 k_B; na imagem ideal estrita (θ_BCC = θ_FCC) ele não ocorre, e só um BCC mais mole (θ_BCC < θ_FCC, abaixo do limiar de abertura 0,61) abre τ_ss < τ_m.
• 6.4 / 6.5 Considerações: a anomalia acústica da água líquida como mudança esbatida de ordem local entre dois tipos de rede, e uma alternativa de campos de vórtices fractais à rotação das micelas (puramente qualitativa).
• 6.6 / 6.8 / 6.9 θ(rede) microscópico em três níveis: modelo de molas de Einstein (θ_BCC/θ_FCC ≈ 0,82…0,94), densidade de estados de fonões do tipo Debye (≈ 0,87…1,00) e rede de Lennard-Jones completa (≈ 1,0; BCC ideal instável ao corte, FCC estado fundamental estável) — todos acima do limiar 0,61, pelo que o polimorfismo real fica fora da imagem ideal harmónica.
• 6.7 Metastabilidade (2D): cálculo de ordem de grandeza para a transformação triangulada ↔ quadrada (ΔU = +1 ε, ΔA ≈ +15,5 %) juntamente com a barreira de nucleação ΔG\* = πλ²/Δμ_A como contributo cinético.

Este capítulo torna plausível a generalização 3D de forma puramente matemática. O 3D nunca é visualizado numa versão HTML. O objetivo é apenas mostrar como o grau de liberdade de rotação e a conectividade mudam de 2D para 3D e que novas transições se tornam assim concebíveis.

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

6.1 Alteração do grau de liberdade de rotação 2D → 3D

O único ponto qualitativamente novo face ao 2D é a rotação dos blocos:

• 2D: um bloco rígido tem 1 GDL de rotação (rotação no plano) → U_rot/N = ½ n_blk τ.
• 3D: um bloco rígido tem 3 GDL de rotação (três eixos de Euler) → U_rot/N = (3/2) n_blk τ.

Ao mesmo tempo a translação passa de 2 a 3 GDL (U_trans/N = (3/2)τ) e o número de coordenação máximo de z_max = 6 (triângulo) a z_max = 12 (empacotamento compacto de esferas). O balanço energético em 3D reza portanto

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

e os valores base da capacidade térmica deslocam-se coerentemente:

A estrutura da curva (patamar sólido, máximo líquido por −(ε/2)dz̄/dτ, queda para o valor do gás) mantém-se idêntica; só os níveis escalam com a dimensão e com z_max.

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

6.2 Influência das diferentes redes cristalinas — a conectividade como traço distintivo

Em 3D há que distinguir pelo menos duas redes ideais, que se separam pela sua conectividade (número de coordenação z_max):

A conectividade z_max é o traço distintivo decisivo: fixa tanto o nível de capacidade térmica do sólido (z_max/4) como o limiar de percolação z̄_c. Como cada rede admite apenas um número de coordenação fixo, a conectividade é uma grandeza quantizada e finita — não há transição contínua entre z_max = 12 e z_max = 8, mas apenas o salto discreto de um tipo de rede para o outro.

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

6.3 A transição adicional sólido(1) → sólido(2)

Da finitude/discretude da conectividade segue uma transição ausente no modelo 2D: sólido(1) → sólido(2) (polimorfismo/alotropia). Duas fases sólidas com z_max diferente (por exemplo FCC ↔ BCC) têm entalpias livres μ(τ, p) distintas; na combinação temperatura/pressão μ_{sólido(1)} = μ_{sólido(2)} a substância muda de rede sem fundir. Na imagem do modelo a geometria de coordenação muda abruptamente (não apenas a ocupação z̄ como na fusão) — a única situação em que a geometria de referência, de resto fixa, salta ela própria. Seria observável um pequeno degrau de c_v (variação de z_max/4) mais um calor latente de transformação, nitidamente mais fraco do que o pico de fusão.

O critério de primeira ordem ln G > 2 + 4ε/K (cap. 5.5) fornece agora uma razão concreta para que seja justamente a interpretação 3D a qualitativamente mais rica: só aí (G = z_max = 12, ln 12 = 2,485 > 2) a fusão pode tornar-se de primeira ordem — e só aí existe a mudança discreta de rede sólido(1) → sólido(2) (FCC ↔ BCC), porque a conectividade é quantizada/finita. Ambos os fenómenos são a mesma causa sob duas vestes: a coordenação mais elevada e discreta do empacotamento 3D gera tanto o salto de fusão nítido como o polimorfismo. A rede 2D isolada (G = 6) não tem nenhum dos dois — coerentemente com o cap. 5.5.

Quantificação (bondkinetics_3d_15.py). A transição é agora calculada de forma enxuta. Ambas as redes recebem uma entalpia livre na imagem harmónica clássica (as partes comuns p·v e de translação cancelam-se na comparação):

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

com os valores característicos sem parâmetros dos cap. 6.1/6.2:

Daí o degrau de c_v sem parâmetros: |Δc_v| = variação de z_max/4 = (12−8)/4 = 1,0 k_B (na transição FCC→BCC o sólido torna-se 1 k_B mais «leve» em capacidade térmica) — muito menor do que o pico de fusão (vários k_B, cap. 5.3/7.11).

A temperatura de transformação τ_ss (μ_FCC = μ_BCC) e o calor latente de transformação L = τ_ss·(S_BCC − S_FCC):

Resultado (constatação honesta). Na imagem ideal estrita (a mesma rigidez fonónica para todas as ligações, θ_BCC = θ_FCC) não há transição sólido(1)→sólido(2): o FCC ganha tanto a entalpia (mais baixa, −6 contra −4 ε) como — pelo maior número de ligações (c_v = 4,5 > 3,5) — a entropia vibracional; μ_FCC mantém-se a toda a temperatura abaixo de μ_BCC. Só uma rigidez fonónica dependente da rede θ_BCC < θ_FCC (o BCC mais aberto vibra mais molemente e tem portanto mais entropia) leva μ_BCC abaixo de μ_FCC a alta temperatura — então a transformação abre-se, e situa-se abaixo do pico de fusão τ_m, pelo que surge a ordem sólido(1)→sólido(2)→líquido→gás (τ_ss < τ_m < τ_c ≤ τ_b(3D), cap. 7.11). A rigidez mole do BCC é uma extensão moderada fora da imagem ideal estrita (conhecida na realidade como o polimorfismo γ/δ do ferro: BCC a alta temperatura, FCC no meio); o degrau de c_v de 1,0 k_B e o pequeno calor de transformação positivo são, pelo contrário, internos ao modelo e sem parâmetros.

Quão pequeno é realmente θ_BCC/θ_FCC? (theta_springs_16.py). Para libertar a transformação do parâmetro livre θ, a rigidez fonónica é agora estimada microscopicamente (cap. 6.6): um modelo de molas de Einstein com molas de par centrais dá θ_BCC/θ_FCC ≈ 0,82 (só vizinhos mais próximos) até ≈ 0,94 (com a segunda camada do BCC). O limiar de abertura exigido pela comparação dos μ situa-se pelo contrário em θ_BCC/θ_FCC < 0,61 (abaixo dele existe τ_ss e fica sob τ_m). A razão geométrica de rigidez é portanto demasiado grande: na imagem ideal estrita o FCC mantém-se estável a toda a temperatura, não há transição sólido(1)→sólido(2) de origem puramente geométrica. O polimorfismo real (ferro γ↔δ) exige contributos adicionais (eletrónicos/magnéticos/anarmónicos) fora da imagem harmónica ideal — precisamente a «extensão moderada» posta acima. Isso confirma a constatação honesta e transforma a razão θ de suposição livre num limite do modelo estimado e nomeado.

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

6.4 Considerações: a anomalia acústica da água líquida

Os líquidos normais mostram v(T) a decrescer monotonamente (cap. 3.3). A água tem um máximo acústico por volta dos 74 °C: v sobe primeiro com T antes de descer. No modelo isso pode ler-se como uma concorrência entre dois tipos de rede locais (cap. 6.2): uma ordem local aberta, de baixa coordenação, semelhante ao gelo (z_max pequeno, muitos «buracos»), e uma ordem local densa, de maior coordenação. A T baixa domina a estrutura aberta; aquecer destrói primeiro a ordem local aberta em favor da mais densa — a conectividade efetiva e com ela K_eff sobem (v cresce). Só acima do máximo predomina a rutura de ligações habitual da FISSURA, K_eff cai, v desce. A anomalia é assim expressão de uma mudança de ordem local do tipo sólido(1)↔sólido(2), esbatida e incompleta, dentro da fase líquida — concebível apenas no modelo 3D com dois tipos de rede, não na rede 2D isolada.

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

6.5 Alternativa: campos de vórtices fractais em vez da rotação das micelas

O movimento browniano local (cap. 3.4) foi explicado pela rotação e troca de micelas. Alternativamente pode representar-se como efeito de campos de vórtices fractais: as fronteiras de bloco, constantemente remodeladas pela FISSURA, geram um campo de correntes locais de momento que — de modo semelhante aos escoamentos turbulentos — apresenta vórtices autossemelhantes em todas as escalas de comprimento. Uma partícula de teste é então impelida não por choques micelares discretos, mas por um contínuo de vórtices invariante de escala. Ambas as imagens (micelas discretas ↔ vórtices fractais) são casos limite complementares da mesma cinética de FISSURA; qual delas prevalece depende da relação entre o tamanho da partícula de teste e o tamanho médio das micelas. Esta breve discussão é puramente qualitativa e não é aqui desenvolvida matematicamente.

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

6.6 Rigidez fonónica microscópica θ(rede): um modelo de molas de Einstein

A comparação sólido(1)→sólido(2) do cap. 6.3 depende da rigidez fonónica relativa θ_BCC/θ_FCC. Em vez de a escolher livremente, ela é aqui estimada a partir da coordenação e da geometria (theta_springs_16.py).

Modelo. Um átomo central está no poço harmónico dos seus z vizinhos; cada ligação é uma mola de par radial κ. Um deslocamento u projeta-se isotropicamente com ⟨cos²⟩ = 1/d = 1/3 (d = 3) sobre as direções de ligação, pelo que a constante elástica efetiva é k_eff = (z/3)·κ. A frequência de Einstein é ω = √(k_eff/m), logo θ = ħω/k_B ∝ √k_eff. Com igual massa e igual mola de par segue

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

Refinamentos. (B) A segunda camada do BCC (6 átomos em r₂/r₁ = 2/√3 = 1,155) com uma mola escalada com a distância κ(r) ∝ (r₁/r₂)ⁿ (n = 6…8, do tipo Lennard-Jones) eleva a razão a 0,91…0,94 (o BCC torna-se efetivamente mais rígido pela próxima segunda camada). (C) A ligação mais curta do BCC a igual densidade (r_nn(BCC)/r_nn(FCC) = 0,972, −2,8 %) enrijece ligeiramente a mola e dá 0,89…0,91. No conjunto uma janela geométrica θ_BCC/θ_FCC ≈ 0,82…0,94.

Limiar de abertura. A comparação dos μ (cap. 6.3) só abre a transformação (τ_ss existe e fica sob τ_m = 1,70) para θ_BCC/θ_FCC < 0,61. A janela geométrica fica nitidamente acima.

Constatação. O modelo de molas puro não torna o BCC suficientemente mole: na imagem ideal estrita o FCC mantém-se estável e uma transição sólido(1)→sólido(2) de origem geométrica não surge. Assim θ_BCC/θ_FCC passou de suposição livre a grandeza estimada, cujo valor natural (≈ 0,8–0,9) mantém a transformação fechada — coerentemente com a constatação honesta do cap. 6.3. O polimorfismo real fica portanto comprovadamente fora da imagem harmónica ideal (contributos eletrónicos/magnéticos, anarmonicidade), e não na mera geometria da rede.

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

6.7 Estimativa matemática da metastabilidade (2D triangulada ↔ quadrada)

Para a discussão qualitativa do cap. 3.2, eis o cálculo de ordem de grandeza correspondente na imagem 2D (discos duros, diâmetro D, distância de contacto a = D).

(1) Energia. Por partícula U/N = −(ε/2)·z (cada ligação a metade): triangulada z = 6 → −3 ε, quadrada z = 4 → −2 ε. A transformação triangulada → quadrada custa portanto ΔU = +1 ε por partícula (endotérmica; a forma frouxa é mais rica em energia).

(2) Volume (área) e trabalho. Área por partícula: triangulada A_tri = (√3/2)·D² ≈ 0,866 D², quadrada A_quad = D². O alargamento é ΔA = A_quad − A_tri ≈ 0,134 D² (+15,5 %) — o aumento de volume mostrado na tira. Contra uma pressão (reduzida) p há que fornecer o trabalho W = p·ΔA; a entalpia livre cresce em Δμ = ΔU + p·ΔA − τ·ΔS. O termo entrópico ΔS > 0 (a forma frouxa tem modos mais moles, um degrau de c_v maior Δc_v = (6−4)/4 = 0,5 k_B) baixa μ_quad abaixo de μ_tri só a τ alta — por isso a triangulada é o estado fundamental de baixa temperatura e a quadrada a forma de alta temperatura / metastável.

(3) Barreira de nucleação (2D). A metastabilidade pressupõe uma barreira. Em 2D um núcleo da nova fase (raio r) cresce com um ganho de área −Δμ_A·π r² (Δμ_A = força motriz por área) contra o contributo de borda +λ·2π r (λ = energia de aresta por comprimento): ΔG(r) = −Δμ_A·π r² + 2π λ r. Máximo em r\* = λ / Δμ_A com a barreira ΔG\* = π λ² / Δμ_A. O núcleo crítico r\ encolhe ao crescer a força motriz Δμ_A; enquanto a escala térmica τ = k_BT não atingir a barreira (**ΔG\ ≳ alguns k_BT), a fase frouxa (ou, na ebulição, a supersaturada) permanece congelada cineticamente. Para pequenos agregados (2–3 partículas, como na tira sobre a supersaturação) a relação borda-superfície é maximamente desfavorável — ficam à esquerda de r\* e dissolvem-se de novo em média. Assim a energia de ativação mostrada na tira é a mesma grandeza ΔG\* que sustenta tanto a forma cristalina metastável como a fase gasosa supersaturada; é um contributo cinético** ao modelo de equilíbrio ideal (sem barreira).

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

6.8 Refinamento do tipo Debye de θ(rede)

O modelo de molas de Einstein (cap. 6.6) substitui toda a distribuição de fonões por uma só frequência. A abordagem presente refina θ(rede) mediante uma densidade de estados de fonões enxuta por rede (theta_debye_17.py): as mesmas molas de par centrais κ(r) = r^(−n) na primeira e na segunda camada de vizinhos formam a matriz dinâmica D(k) = Σ_j κ_j·(R̂_j⊗R̂_j)·(1 − cos(k·R_j)), diagonalizada numa grelha k da zona de Brillouin primitiva (Monkhorst–Pack, 12³ pontos) mediante um resolvedor analítico 3×3, apenas biblioteca padrão. Dos três ramos de fonões seguem duas medidas de rigidez do tipo Debye como momentos de frequência:

• 2.º momento θ(2) ∝ √⟨ω²⟩ — o limite de Einstein (a média na ZB de 1 − cos dá 1, logo ⟨ω²⟩ = z·κ/3); controlo, só vizinhos mais próximos, pura geometria: √⟨ω²⟩_BCC / √⟨ω²⟩_FCC = 0,816 = √(8/12) — exatamente o valor de Einstein do cap. 6.6.
• momento logarítmico θ(0) ∝ exp⟨ln ω⟩ — a grandeza propriamente do tipo Debye, que pondera os ramos acústicos moles que o valor de Einstein ignora.

Todas as constantes dotadas de escala (ħ/k_B, a escala de κ, o prefator de Debye) cancelam-se na razão θ_BCC/θ_FCC; só permanece a forma da DOS. Com a segunda camada de vizinhos (que é o que estabiliza mecanicamente o BCC) e uma lei de mola n = 0…8 obtém-se

θ_BCC/θ_FCC ≈ 0,87…1,00 (banda de erro sobre a definição do momento e a lei da mola).

Isso é antes mais alto do que o puro valor de Einstein (0,82…0,94), porque a ligação mais curta do BCC a igual densidade enrijece as molas. O limiar de abertura 0,61 exigido pela comparação dos μ (cap. 6.3) mantém-se nitidamente abaixo: também a DOS de fonões não torna o BCC suficientemente mole. A constatação do cap. 6.6 fica assim confirmada e precisada pela aproximação mais fina — na imagem ideal estrita o FCC mantém-se estável a toda a temperatura; não surge qualquer transição sólido(1)→sólido(2) de origem puramente geométrica, e o polimorfismo real (ferro, cap. 7.13) precisa de contributos fora da imagem harmónica ideal.

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

6.9 θ(rede) a partir de uma rede de Lennard-Jones completa

O cap. 6.8 usa uma mola central truncada κ(r) = r^(−n) em nn + 2nn. A versão 17 examina a mesma questão com uma verdadeira constante de força a partir de um potencial de pares sobre a rede completa (theta_ljlattice_18.py): cada ligação torna-se uma mola derivada do potencial de Lennard-Jones φ(r) = 4ε[(σ/r)¹² − (σ/r)⁶], com a parte longitudinal κ_L = φ″(r) e a parte transversal κ_T = φ′(r)/r. A matriz dinâmica

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

soma sobre todas as camadas de vizinhos até um raio de corte (2,0…3,0 nn); cada rede é previamente relaxada até ao seu próprio equilíbrio LJ (constante de rede livre, densidade igual não imposta). É o teste honesto para saber se a aproximação enxuta nn+2nn (cap. 6.8) acerta a ordem das rigidezes.

Resultado (determinista, apenas biblioteca padrão):

Duas constatações, ambas a reforçar os cap. 6.6/6.8:

1. O FCC é o estado fundamental LJ estável (E ≈ −8,4 ε por átomo com corte 3,0 nn, contra −8,0 ε para o BCC) — o conhecido primado do empacotamento denso num potencial de pares.
2. O BCC LJ ideal é instável ao corte: cerca de 3…3,5 % dos modos de fonões são imaginários (ω² < 0, instabilidade de Born). A imagem ideal do puro potencial de pares não porta portanto nenhum BCC mecanicamente estável; sobre os modos reais restantes θ_BCC/θ_FCC ≈ 1,0 — o limiar de abertura 0,61 exigido pela comparação dos μ por maioria de razão não é ultrapassado por baixo.

Assim a afirmação dos cap. 6.6/6.8 é robusta face ao refinamento da mola truncada para a rede completa e do modelo puramente longitudinal para a matriz dinâmica completa: na imagem ideal clássica não surge nenhuma transição sólido(1)→sólido(2) puramente geométrica — o polimorfismo real (ferro cap. 7.13, titânio cap. 7.15) precisa de contributos fora do potencial de pares (eletrónicos/magnéticos, anarmónicos). Um verdadeiro cálculo de fonões ab initio é deliberadamente deixado de fora (abandonaria o caráter enxuto da imagem ideal).

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

7 — Calibração e equação de estado

Resumo. O capítulo 7 calibra o modelo em substâncias reais mediante o único parâmetro de energia ε e desenvolve o ramo líquido/gás numa equação de estado:

• 7.1 / 7.2 Método e tabela de calibração: acoplamento linear de temperatura T = (ε/k_B)·τ com um ponto de ancoragem no ponto de fusão; daí a razão universal refinada R = τ_b/τ_m ≈ 1,54.
• 7.3 / 7.4 Equação de estado de van der Waals com caroço duro: a EOS reduzida p_r = 8τ_r/(3v_r − 1) − 3/v_r² com ponto crítico (1, 1, 1) e coexistência de Maxwell; a ∝ ε, b ∝ D². Os pontos de ebulição situam-se sistematicamente em T_b/T_c ≈ 0,48 (distância de Guldberg); um cálculo de comparação Redlich–Kwong eleva-os a ≈ 0,60 — como orçamento de erro nomeado, sem mudança de EOS.
• 7.5 Superfície (τ, v) comum: a cúpula de coexistência (ebulição), a linha de fusão e a parede dura compõem-se com um ponto triplo (τ_r,triple ≈ 0,47 para v_solid = 0,40).
• 7.6 Calor latente: L ≈ 1,41 ε dá ΔS_fus ≈ 1,175 R (perto da regra de Richards); os metais são acertados excelentemente (cobre ×1,00), ao passo que os cristais iónicos, de pontes de hidrogénio e moleculares ficam fora da imagem ideal.
• 7.7 Unificação ε ↔ ε_lg: k_B T_c = z_max·ε/4 acopla ambas as escalas (ε_lg = (27/32)·z_max·ε); a previsão 3D τ_c/τ_m = 2,50 acerta com as substâncias moleculares simples, ficando o 2D como caso especial.
• 7.8 / 7.9 Ponto triplo sem parâmetros: v_solid segue do limiar de rigidez z̄_c = 2d, não da conectividade; em 3D abre-se uma janela líquida limpa (τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00), em 2D permanece fechada (papel especial do 2D, agora triplamente documentado).
• 7.10 / 7.11 Varrimentos 3D reais: a percolação de ligações em FCC/BCC dá z̄_c(3D) ≈ 1,5–1,6; o varrimento em τ da cinética de ligação mostra τ_m ≈ τ_rig (fusão = rigidez, baixo) e τ_b(3D) ≈ τ_c (ebulição = desagregação da conectividade no ponto crítico).
• 7.12 / 7.13 / 7.15 Linha sólido(1)→sólido(2): do ponto τ_ss a pressão faz uma linha (dp/dτ = ΔS/Δv > 0, Δv = +8,9 %), inserida como quarta linha de fronteira na superfície (τ, v) e calibrada como linha (τ, p) com ferro (γ→δ), confirmada com titânio (α→β) e verificada com árgon como controlo negativo.
• 7.14 Ramo Redlich–Kwong: como variante autónoma da imagem ideal com calibração própria ao lado da opção por omissão de van der Waals (T_b/T_c ≈ 0,60, erro −1 %).

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

7.1 Método de calibração (um parâmetro de energia ε) — ajustado

Linear: T = (ε/k_B)·τ, um ponto de ancoragem por substância. Ancoragem no ponto de fusão: ε/k_B = T_m,exp/τ_m, previsão T_b,pred = T_m,exp·(τ_b/τ_m). Com o pico de fusão refinado (termo cooperativo) o máximo de c_v situa-se agora em τ_m = 1,20 (em vez de 1,15), mantendo-se a percolação de conectividade (perto da ebulição) inalterada em τ_b = 1,85. Daí a razão universal refinada:

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

A tabela de calibração da versão 08 (cap. 7.2) escala linearmente com R: todos os T_b,pred descem pelo fator 1,54/1,61 ≈ 0,957 (≈ −4,3 %). Isso melhora ligeiramente a sobrestimação dos gases nobres (árgon +54 % → +48 %) e piora um pouco os acertos de NaCl e ferro (NaCl −0,6 % → −4,8 %). O enunciado de base mantém-se: R serve para substâncias de empacotamento denso isotrópico e ligação não direcional; gases nobres e cadeias ficam fora da imagem ideal.

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

7.2 Tabela de calibração

Retomada inalterada de modell_08.md cap. 7.2 (com a nota sobre o escalonamento linear em R do 7.1). Reproduzível via cv_scan_08.py (tabela) e cv_scan_09.py (τ_m, τ_b, R).

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

7.3 Eixo densidade/volume: equação de estado de van der Waals com caroço duro

A pressão passa de «parâmetro» a verdadeira variável. O CAROÇO DURO fornece o volume próprio b, a atração de ligação (ε·z̄) a pressão de coesão a. Em unidades reduzidas (ponto crítico em (v_r, τ_r, p_r) = (1, 1, 1)):

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

Coexistência por construção de Maxwell (cv_scan_09.py, corte de áreas iguais):

Resultado. Abaixo de τ_c abre-se o salto de densidade l/g (v_gas/v_liq cresce de 2,5 a 11,5, tanto maior quanto mais frio). Essa é a linha de ebulição — agora nitidamente separada do pico de c_v da fusão (cap. 5.3), que vive numa rede de densidade fixa. Assim fecha-se o ponto 3 deixado em aberto na versão 08: fusão (ordem de rigidez/ligação a densidade fixa) e ebulição (salto de densidade/volume) são eixos diferentes do mesmo modelo.

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

7.4 Calibração quantitativa l/g: a → ε, b → D²

Aplicação das grandezas reduzidas a substâncias reais. As constantes de van der Waals dependem do ponto crítico:

b = k_B·T_c / (8 p_c) (volume próprio, CAROÇO DURO), a = 27 k_B²·T_c² / (64 p_c) (coesão, LIGAÇÃO).

Assim a ∝ ε (a coesão provém da energia de ligação) e b ∝ D² (o volume próprio 2D do diâmetro do caroço duro). A escala de energia natural do ramo l/g é

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

Substância de ancoragem árgon (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 e um diâmetro molecular D ≈ b^{1/3} = 3,8·10⁻¹⁰ m — com uma exatidão de 11 % face ao diâmetro de van der Waals do árgon (~3,4·10⁻¹⁰ m).

Previsão da temperatura de ebulição. A pressão exterior fixa, o ponto de ebulição está onde se cumpre a pressão de saturação reduzida p_sat(τ_r)·p_c = p_atm; daí T_b = T_c·τ_r,b (curva de saturação de Maxwell de lg_calib_10.py):

*O CO₂ sublima a 1 atm (T = 194,7 K = ponto de sublimação).

Resultado. Para as substâncias moleculares (p_c ≈ 30–70 bar) o modelo dá em toda a parte T_b/T_c ≈ 0,48–0,51, sistematicamente ~17–24 % abaixo do valor medido. Isso não é um erro do modelo em sentido estrito, mas a conhecida distância da equação de van der Waals à regra empírica de Guldberg (T_b/T_c ≈ 0,6): a isotérmica de van der Waals subestima a pressão de saturação a baixas temperaturas. A água (e os metais, cujos T_c/p_c estão de qualquer modo apenas estimados grosseiramente) ficam com p_r,atm < 0,005 fora do intervalo sensato de van der Waals (τ_r,b < 0,40) — aí a coesão já não é do tipo van der Waals (pontes de hidrogénio, ligação metálica). A ordem de grandeza e a ordenação dos pontos de ebulição são acertadas corretamente.

Enquadramento de Guldberg como tabela de erros. A distância é aqui apenas nomeada, não corrigida por uma mudança de EOS (a imagem ideal mantém-se do tipo van der Waals). A regra empírica de Guldberg (Guldberg 1890) diz que o ponto de ebulição se situa por volta de T_b/T_c ≈ 0,6 (com ligeira dependência da substância, ~0,58–0,64). O modelo reduzido de van der Waals dá em vez disso um valor constante T_b/T_c = τ_r,b ≈ 0,48–0,51:

Leitura. A média do modelo, 0,49, acerta a constante de Guldberg 0,60 a cerca de −18 %; o desvio é sistemático e unilateral (sempre demasiado baixo), portanto uma propriedade da isotérmica de van der Waals (pressão de saturação demasiado pequena a τ baixos), não da imagem molecular. Uma melhor EOS de coesão (por exemplo Redlich–Kwong) elevaria τ_r,b — mas fica aqui deliberadamente fora da imagem ideal.

Orçamento de erro com um cálculo de comparação Redlich–Kwong (guldberg_rk_15.py). Para cifrar a distância de van der Waals à regra de Guldberg (não para a sanar), o mesmo teste do ponto de ebulição é repetido com uma coesão do tipo Redlich–Kwong, dependente da temperatura, a/√T. Em unidades reduzidas (ponto 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,

com construção de Maxwell de áreas iguais e a mesma pressão exterior reduzida p_r,atm de cima. Resultado (um valor por substância de ancoragem, T_b/T_c = τ_r,b):

Resultado. A coesão RK dependente da temperatura eleva T_b/T_c de vdW ≈ 0,48 para ≈ 0,60 — praticamente a constante de Guldberg e os valores reais (≈ 0,58–0,63); o erro médio face às medições cai de −21 % (vdW) para −1 % (RK). Assim a causa da distância de Guldberg fica claramente cifrada: é justamente a falta de dependência da temperatura da coesão de van der Waals (pressão de saturação demasiado pequena a τ baixos), e não a imagem molecular. Este cálculo de comparação fica deliberadamente fora da imagem ideal; o modelo mantém-se do tipo van der Waals (sem mudança de EOS) — o orçamento de erro é apenas nomeado, não incorporado.

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

7.5 Superfície (τ, v) comum com ponto triplo

Os eixos de fusão e de ebulição viveram até agora separados: a fusão (ordem de rigidez/ligação) a densidade fixa, a ebulição (van der Waals) através do salto de densidade. Num plano comum (τ_r, v_r) compõem-se:

• Cúpula de coexistência (ebulição): v_liq(τ_r) (ramo esquerdo) e v_gas(τ_r) (ramo direito) da construção de Maxwell, que se encontram no ponto crítico (τ_r, v_r) = (1, 1).
• Sólido (fusão): faixa estreita logo acima da parede dura v_r > 1/3, a densidade elevada e quase fixa (volume reduzido v_solid ≈ 0,40).
• Linha de fusão: quase vertical em v_solid (densidade fixa), à temperatura da rigidez de Maxwell (z̄_c = 2d; o cap. 7.8 dela deriva v_solid sem parâmetros — a antiga atribuição da versão 11, «z̄_c ≈ 2,08», fica assim corrigida).

O ponto triplo é o ponto em que a fronteira de coexistência líquida v_liq(τ_r) atinge a densidade de empacotamento sólido v_solid — abaixo já não existe ramo líquido (sublimação direta):

Resultado. Para v_solid = 0,40 o ponto triplo situa-se em τ_r,triple ≈ 0,47, p_sat ≈ 0,018, com um volume de vapor v_gas ≈ 68 (vapor de sublimação muito diluído). A razão τ_triple/τ_c ≈ 0,47 situa-se na janela real (água 0,42; árgon 0,56; CO₂ 0,71). Quanto mais denso o sólido (v_solid menor), tanto mais baixo o ponto triplo reduzido e tanto mais diluído o vapor coexistente — qualitativamente correto. Assim ambas as transições e o ponto triplo ficam ancorados numa única superfície (τ, v); o ponto crítico assenta em (1, 1, 1). A versão 15 acrescenta no cap. 7.12 uma quarta linha de fronteira — a transformação sólido(1)/sólido(2) (FCC↔BCC) como duas faixas de densidade sólida quase verticais (v ≈ 0,667 e 0,726) que se convertem uma na outra em τ_ss.

Versão anterior no desenvolvimento: aqui v_solid = 0,40 estava ainda imposto. O cap. 7.8 deriva v_solid sem parâmetros do limiar de percolação z̄_c e faz do ponto triplo uma verdadeira previsão.

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

7.6 Calor latente face às entalpias de fusão

O calor latente do modelo L ≈ 1,41 ε (cap. 5.3) é traduzido em unidades reais mediante a calibração no ponto de fusão ε/k_B = T_m/τ_m (τ_m = 1,20):

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

Isso é — notavelmente — quase exatamente a regra de Richards (ΔS_fus ≈ R para sólidos de construção simples). Comparação com entalpias de fusão medidas (lg_calib_10.py):

Resultado. O modelo acerta em toda a parte a ordem de grandeza (kJ/mol) e excelentemente a regra de Richards para metais simples (cobre ×1,00, ferro ×1,28); os gases liquefeitos ficam logo abaixo (metano ×0,94, azoto ×0,87, árgon ×0,69). Os cristais iónicos, de pontes de hidrogénio e moleculares (NaCl ×0,37, água ×0,44, benzeno ×0,28, CO₂ ×0,23) têm entropias de fusão reais nitidamente maiores (ΔS_exp/R = 2,6–5,0), porque ao fundir se libertam graus de ordem adicionais (emparelhamento iónico, orientação de moléculas poliatómicas) que o modelo de um só parâmetro de ordem não capta. O limite fica assim claramente nomeado: a imagem ideal vale para substâncias de ligação isotrópica e não direcional (sem ajuste).

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

7.7 Unificação das escalas de energia: ε ↔ ε_lg

Até agora o ramo de fusão portava a energia de ligação ε (τ_m = 1,20) e o ramo l/g uma escala de coesão separada ε_lg/k_B = 3,375·T_c (cap. 7.4). Ambas se podem derivar de um só parâmetro do modelo. No modelo a coesão de van der Waals a nasce da soma das energias de ligação por molécula (LIGAÇÃO): cada molécula porta z_max/2 ligações de energia ε. Um gás em rede de Bragg–Williams em campo médio com atração entre vizinhos ε fornece a relação clássica para a temperatura crítica:

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

Inserindo isto em ε_lg = a/b = (27/8)·k_B·T_c (cap. 7.4) segue o acoplamento de ambas as escalas:

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

Assim pico de fusão, percolação e ponto crítico situam-se numa única escala ε (unify_scale_12.py):

Previsão τ_c/τ_m = z_max/(4·τ_m): 2D → 1,25, 3D → 2,50. Face a substâncias reais:

Resultado. A previsão 3D τ_c/τ_m = 2,50 acerta bem o valor característico real das substâncias moleculares simples (T_c/T_m ≈ 2,0–2,1); os metais ficam mais altos (T_c estimado grosseiramente, outro tipo de ligação). A previsão 2D (1,25) é demasiado baixa — e τ_c(2D) = 1,50 fica até abaixo de τ_b = 1,85 (o ponto crítico ficaria sob o limiar de percolação). Não é erro de cálculo, mas o conhecido papel especial do 2D: aí as escalas de coesão e de ligação não se ordenam limpamente, falta uma janela líquida robusta (cf. a fusão 2D contínua, cap. 5.5). Só em 3D (τ_m = 1,20 < τ_c = 3,00) as escalas se separam de forma fisicamente sensata.

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

7.8 Ponto triplo sem parâmetros a partir de z̄_c

O cap. 7.5 ainda impunha v_solid = 0,40. A versão anterior no desenvolvimento deriva a densidade de empacotamento sólido do limiar de percolação. O número de coordenação de uma rede a volume reduzido v_r segue a relação densidade-coordenação (parede dura v_r = 1/3 ↔ empacotamento mais denso, z̄ = z_max):

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

O sólido perde a sua rigidez logo que z̄ desce abaixo do limiar z̄_c; a densidade correspondente é v_solid. O ponto triplo é então o ponto em que a fronteira de coexistência líquida v_liq(τ_r) atinge essa densidade (unify_scale_12.py, 2D, z_max = 6):

Versão melhorada. A densidade de fusão v_solid é fixada pela rigidez genérica de Maxwell z̄_c = 2d (2D → 4, 3D → 6), não pela conectividade das ligações (fundamentação no cap. 5.4a: a fusão é perda de rigidez, não de conexidade). Assim o ponto triplo fiável reza (rigidity_3d_13.py, 2D):

Resultado 2D. Mesmo com o limiar de rigidez (fisicamente correto) a rede 2D isolada não dá nenhuma janela líquida limpa: v_solid = 0,50 fornece de facto o valor real τ_triple/τ_c ≈ 0,77, mas em unidades ε absolutas τ_triple = 0,77·1,50 = 1,16 e fica portanto abaixo de τ_m = 1,20. O limiar de conectividade (z̄_c = 2,08) empurra até v_solid ≈ 0,96 ≈ v_c para fora da cúpula de coexistência — ponto triplo e ponto crítico coincidem. Ambos os cálculos confirmam o papel especial do 2D (cap. 5.5, 7.7): a rede 2D isolada dificilmente sustenta uma fase líquida autónoma. É uma constatação honesta, não uma fraqueza de cálculo ou de modelo em sentido estrito.

O teste decisivo é portanto o cálculo de controlo 3D (cap. 7.9): aí a mesma regra de rigidez (z̄_c = 6 → v_solid = 0,667) liberta uma janela limpa. Em qualquer caso o ponto triplo está agora ligado sem parâmetros a z̄_c — o parâmetro livre v_solid está eliminado.

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

7.9 Cálculo de controlo 3D: uma janela líquida limpa numa única escala ε

Verificação de se o colapso 2D se resolve em 3D quando todas as marcas se situam numa única escala ε (τ_c = z_max/4 do cap. 7.7; v_solid da rigidez genérica z̄_c = 2d; ponto triplo onde a densidade líquida reduzida de van der Waals v_liq(τ_r) atinge a densidade do sólido). Resultado (rigidity_3d_13.py):

Resultado. Em 3D abre-se uma janela líquida limpa: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00 (larguras: τ_triple − τ_m = 1,63; τ_c − τ_triple = 0,18). O bordo superior é estreito — o ponto triplo fica logo abaixo do crítico — mas a ordem está correta e a fase líquida existe como campo autónomo. Em 2D a janela permanece fechada (τ_triple = 1,16 < τ_m). Assim o papel especial do 2D fica documentado pela terceira vez de forma independente e mostra-se ao mesmo tempo que é um efeito dimensional: a mesma regra de rigidez sana o modelo em 3D. O estreito bordo superior é corroborado no cap. 7.10 por um varrimento de percolação 3D real e uma análise de sensibilidade.

Caixa dos limites do modelo — o papel especial do 2D. Três cálculos independentes mostram a mesma constatação: a rede 2D isolada dificilmente sustenta uma fase líquida autónoma. 1. τ_c < τ_b (cap. 7.7): o ponto crítico (1,50) ficaria abaixo da percolação de conectividade (1,85) — as escalas de coesão e de ligação não se ordenam limpamente em 2D. 2. v_solid ≈ v_c (cap. 7.8): a densidade de fusão de conectividade (0,96) coincide quase com a densidade crítica — ponto triplo e ponto crítico colapsam. 3. τ_triple < τ_m (cap. 7.9): mesmo com a rigidez de Maxwell correta o ponto triplo (1,16) fica abaixo do pico de fusão (1,20) — sem janela.

É uma constatação honesta, não uma fraqueza de cálculo ou de modelo: é o conhecido efeito dimensional (fusão 2D contínua, do tipo KTHNY; cap. 5.5). Em 3D a mesma regra de rigidez sana o modelo (cap. 7.9/7.10). Se uma segunda ordem local (cap. 6.4) poderia melhorar adicionalmente a imagem 2D permanece uma questão de extensão, não um defeito da imagem ideal. Dado que o modelo deve servir à ilustração, aceito a fraqueza e renuncio a uma ampliação axiomática do modelo.

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

7.10 Varrimento real de percolação de ligações em 3D (FCC/BCC) e sensibilidade do bordo da janela

Os cap. 7.8/7.9 transferiam mentalmente o valor de conectividade 2D z̄_c = 2,08 para o 3D (onde v_solid = 1,92, não físico) e apoiavam a afirmação 3D na rigidez genérica z̄_c = 2d. A versão 13 examina agora a conectividade em 3D diretamente: percolação de ligações em FCC (z_max = 12) e BCC (z_max = 8) por Monte Carlo com union-find, bordos periódicos, semente LCG determinista (rigidity_3d_14.py). Limiar p_c = cruzamento P_∞(p) = 0,5 (a mesma metodologia do estudo de tamanho finito 2D do cap. 5.6):

Devido à aresta de caixa moderada (L = 12, Python puro) os valores MC ficam ligeiramente acima dos limiares da literatura (viés de tamanho finito, como em 2D por baixo ou por cima consoante o estimador); a ordenação e a ordem de grandeza estão corretas (FCC mais denso → p_c menor do que o BCC).

Leitura (confirmação 3D independente). Os limiares de conectividade 3D situam-se em z̄_c(3D) ≈ 1,5–1,6 — até mais baixos do que o valor 2D 2,08 e muito abaixo da rigidez de Maxwell 2d = 6. A «densidade de fusão de conectividade» daí decorrente, v_solid = 1/(3 p_c) ≈ 1,7 (BCC) ou 2,6 (FCC), fica muito fora da cúpula de coexistência (v_liq,max ≈ 0,94). Assim mostra-se diretamente — sem transferir o valor 2D para o 3D — que em 3D a conectividade por maioria de razão não pode fixar a densidade de fusão; esta tem de vir da rigidez (z̄_c = 6 → v_solid = 0,667). Os dois limiares estão ainda mais nitidamente separados em 3D do que em 2D. O verdadeiro τ_b(3D) (temperatura à qual z̄(τ) atinge o baixo limiar de conectividade 1,5) ficaria portanto acima do valor 2D τ_b, muito perto do ponto crítico — coerentemente com a imagem «desagregação da conectividade = ebulição».

Sensibilidade do estreito bordo superior da janela (τ_c − τ_triple ≈ 0,18). Numericamente em torno de v_solid = 0,667 (±2 %) obtém-se dτ_triple/dv_solid ≈ 1,55. A janela só se fecharia (τ_triple = τ_c = 3,00) se v_solid crescesse em ≈ 0,11 (até ≈ 0,78), isto é, se o limiar de rigidez 3D efetivo estivesse sensivelmente acima de 2d = 6 (z̄_c ≳ 7,7). Para a contagem genérica de Maxwell (z̄_c = 2d = 6) a ordem τ_m < τ_triple < τ_c mantém-se portanto robusta; o bordo superior é realmente estreito, mas não está ameaçado. A mudança discreta de rede sólido(1)→sólido(2) (FCC↔BCC) está quantificada na versão 14 (cap. 6.3); o verdadeiro varrimento em τ da cinética de ligação em 3D segue no cap. 7.11.

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

7.11 Varrimento real em τ da cinética de ligação em 3D: τ_b(3D) e a linha de fusão numa escala ε

Os cap. 7.8–7.10 apoiavam a afirmação 3D na rigidez genérica (z̄_c = 2d = 6) e no varrimento estático de conectividade (z̄_c(3D) ≈ 1,55). Ficava em aberto o verdadeiro varrimento em τ da cinética de ligação em 3D: onde, no eixo das temperaturas, o z̄(τ) autoconsistente atinge esse baixo limiar de conectividade (= o verdadeiro τ_b(3D)), e onde está o pico de fusão? A versão 14 calcula o z̄(τ) autoconsistente (cap. 4.2, x = 1/(1 + G·e^{−(ε+K·x)/τ})) diretamente em FCC (z_max = 12) e BCC (z_max = 8) com G = z_max e a mesma cooperatividade J = 0,55 do caso de referência 2D; n_blk(τ) provém de um Monte Carlo de ligações FCC/BCC (union-find), e

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

c_v = dU/dτ (bondkinetics_3d_15.py, semente LCG determinista; n_blk suavizado, já que só o fraco termo de rotação de bloco dele depende). São lidas três marcas: o pico de c_v τ_m (fusão), a marca de rigidez τ_rig (z̄(τ) cruza 2d = 6) e o verdadeiro τ_b(3D) (z̄(τ) cruza o limiar de conectividade z̄_c(3D) do cap. 7.10):

Resultado. Duas constatações, ambas coerentes com a imagem até aqui:

1. A fusão fica em baixo e está duplamente documentada. O pico termodinâmico de fusão τ_m (FCC 1,70, BCC 1,40) coincide praticamente com a marca de rigidez τ_rig (z̄ = 6; FCC 1,73, BCC 1,35) — a mesma temperatura à qual a rede perde a sua rigidez. Isso ancora a linha de fusão 3D na escala ε: T_m = (ε/k_B)·τ_m, ou seja, por exemplo para o árgon (T_m = 83,8 K) ε/k_B = 83,8/1,70 ≈ 49 K. Como K = J·z_max é maior em 3D do que em 2D (6,60 em vez de 3,30), o pico de fusão 3D fica mais alto do que o valor 2D 1,20 — não é contradição, mas a cooperatividade mais forte do empacotamento mais denso.
2. O verdadeiro τ_b(3D) assenta no ponto crítico. Como o limiar de conectividade 3D é tão baixo (≈ 1,55), z̄(τ) só o atinge bem no topo: τ_b(3D) = 3,18 (FCC) ou 2,85 (BCC) — ao nível de, ou ligeiramente acima de, τ_c = z_max/4 (3,00 ou 2,00). O aglomerado gigante de ligações sobrevive portanto a toda a fase líquida e só se desagrega no ponto crítico. É exatamente a previsão do cap. 7.10 («o verdadeiro τ_b(3D) ficaria acima do valor 2D, muito perto do ponto crítico») e confirma de forma independente a imagem «ebulição = desagregação da conectividade»: em 3D a matéria permanece ligada de ponta a ponta até muito perto do ponto crítico.

Ambas as constatações corroboram a separação dos limiares do cap. 5.4a: fusão ↔ rigidez (baixo, τ_m ≈ τ_rig), ebulição ↔ conectividade (alto, τ_b(3D) ≈ τ_c). A densidade de fusão vem portanto da rigidez (v_solid = 0,667), não da conectividade — os dois limiares estão afastados no eixo das temperaturas em 3D mais do que nunca antes no modelo.

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

7.12 A linha de transformação sólido(1)→sólido(2) na superfície (τ, v)

O cap. 6.3 fornecia a transformação de rede FCC↔BCC como ponto τ_ss (a p = 0). Através da pressão resulta daí uma linha (polymorph_line_16.py), que pertence como quarta linha de fronteira à superfície (τ, v) comum dos cap. 7.5/7.9 — a par da cúpula de coexistência (ebulição), da linha de fusão e da parede dura.

Salto de volume. Os dois sólidos assentam a densidade quase fixa; a sua diferença de volume segue das frações de empacotamento (v ∝ 1/η): η_FCC = 0,7405, η_BCC = 0,6802. Com v_FCC = 0,667 (cap. 7.9) obtém-se v_BCC = v_FCC·η_FCC/η_BCC = 0,726, logo

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

O BCC mais aberto tem o maior volume — a recristalização FCC→BCC está ligada a um aumento de volume (o mesmo enunciado da tira «recristalização»).

Clausius–Clapeyron. Com ΔS = S_BCC − S_FCC (fase de alta menos fase de baixa temperatura), o declive reza

dp/dτ = ΔS/Δv.

Ambas as grandezas são positivas (o BCC é a fase de alta temperatura, tem portanto mais entropia e o maior volume), logo dp/dτ > 0: uma pressão maior eleva a temperatura de transformação, porque favorece o FCC mais denso. Para a referência θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, L = 0,667 ε):

Linha em (τ, p) para a referência (crescente): p = 0 → τ_ss = 1,33; p = 2 → 1,57; p = 4 → 1,81; p = 6 → 2,04 (unidades reduzidas).

Na superfície (τ, v) a transformação aparece como duas faixas de densidade sólida quase verticais — sólido(1) = FCC em v ≈ 0,667 (mais denso, T baixa) e sólido(2) = BCC em v ≈ 0,726 (mais aberto, T alta) — que se convertem uma na outra em τ_ss; acima segue, em τ_m, a fusão. A ordem completa no eixo τ é portanto τ_ss < τ_m < τ_c ≤ τ_b(3D) (cap. 7.11).

Unidades reais. Mediante a calibração no ponto de fusão ε/k_B = T_m/τ_m (τ_m = 1,70, FCC 3D) o degrau de c_v e o calor de transformação tornam-se palpáveis. O degrau de c_v é independente da calibração: |Δc_v| = 1,0 k_B = 1 R = 8,31 J/(mol·K) = 0,0083 kJ/(mol·K). O calor de transformação ΔH_ss = L·R·(ε/k_B) (referência L = 0,667 ε):

Os calores de transformação são pequenos e positivos (ordem de 0,2–0,9 kJ/mol) e ficam portanto nitidamente abaixo das entalpias de fusão (cap. 7.6, vários kJ/mol) — exatamente a imagem esperada de uma transformação polimórfica fraca com um pequeno degrau de c_v. (A razão θ de 0,60 fica deliberadamente fora da imagem ideal estrita; o cap. 6.6 mostra que a pura geometria daria θ ≈ 0,82–0,94 e a transformação permaneceria então fechada — a linha é a consequência se-então de uma fase BCC mole.)

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

7.13 Linha de transformação sólido(1)→sólido(2) no diagrama (τ, p), calibrada com ferro

O cap. 7.12 leva a transformação como linha (τ, v); a versão 16 acrescenta-a como linha na imagem (τ, p) e calibra-a com um polimorfo FCC↔BCC real (solidsolid_pt_17.py).

Campo de pressão da EOS de van der Waals com caroço duro. Nas duas densidades sólidas v_FCC = 0,667 e v_BCC = 0,726 (cap. 7.12) a EOS reduzida p_r = 8τ_r/(3v_r − 1) − 3/v_r² (cap. 7.3) dá em toda a parte pressões positivas; a igual τ_r o FCC mais denso porta a pressão maior (τ_r = 1: p_r = 1,25 contra 1,10). A compressão favorece portanto o FCC — coerentemente com o declive positivo de Clausius–Clapeyron dp/dτ = ΔS/Δv > 0.

Linha e ponto de transformação. Para a referência θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, dp\/dτ = 8,47 em unidades reduzidas p\ = p·D³/ε) a linha cresce como p\(τ) = 8,47·(τ − 1,333): p\ = 0 → τ = 1,33; p\ = 2 → 1,57; p\ = 4 → 1,81; p\ = 6 → 2,04. Em **p\_tri ≈ 3,1 a linha sólido-sólido encontra a marca de fusão τ_m = 1,70 — um ponto de transformação sólido(1)-sólido(2)-líquido** no plano (τ, p) (análogo do ponto triplo).

Calibração com ferro (ε/k_B = T_m/τ_m = 1811 K / 1,70 = 1065 K). O ferro tem duas transições FCC↔BCC de sinal oposto — precisamente o caso de teste:

Para γ→δ (τ_ss = 1667/1065 = 1,565 → θ_BCC/θ_FCC ≈ 0,61) o modelo acerta sinal e ordem de grandeza: ΔH_ss ≈ 4,2 kJ/mol (na realidade ≈ 0,84 kJ/mol — a ordem de grandeza kJ/mol, com o valor cerca de 5× demasiado alto pela aproximação de limiar abrupta), degrau de c_v 1 R = 8,31 J/(mol·K) sem parâmetros, dT/dp ≈ +245 K/GPa. A transição α→γ, pelo contrário, tem a fase densa (γ) como fase de alta temperatura — o que a pura geometria de empacotamento não pode fornecer (mais denso = mais baixo = T baixa); o ferro α é estabilizado magneticamente. O ferro confirma assim ambos os lados da constatação honesta: a topologia geométrica (γ→δ) está dentro da imagem ideal, a anomalia magnética (α→γ) fora — o mesmo «fora» que os cap. 6.6/6.8 cifram com θ ≈ 0,87…1,00 > 0,61.

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

7.14 Redlich–Kwong como ramo separado da imagem ideal

A distância de Guldberg do ramo de van der Waals (cap. 7.4, T_b/T_c ≈ 0,48 em vez de ≈ 0,60) era aí apenas nomeada. A pedido expresso do utilizador, uma versão melhorada leva a coesão de Redlich–Kwong como variante autónoma da imagem ideal com calibração própria, coexistência própria e secção de erro própria (rk_branch_17.py) — o ramo de van der Waals mantém-se a opção por omissão; toda a restante física do modelo (fusão, percolação, sólido(1)→sólido(2)) é indiferente à escolha da EOS de coesão.

Escala de energia própria. De a_RK = Ω_a·R²T_c^{2,5}/p_c e b_RK = Ω_b·R T_c/p_c segue a energia de coesão do ramo:

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

maior, porque a coesão dependente da temperatura a/√T dobra o fator √T_c. O diâmetro que segue de b_RK acerta os valores reais de modo semelhante ao ramo vdW (árgon D ≈ 3,34·10⁻¹⁰ m).

Coexistência própria (reduzida, ponto 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):

Secção de erro própria (ponto de ebulição em p_r,atm = p_atm/p_c, T_b/T_c = τ_r,b):

Resultado. O ramo RK eleva T_b/T_c de vdW ≈ 0,48 para ≈ 0,60 e acerta com isso praticamente a constante de Guldberg e os valores reais (erro médio −1 % contra vdW −21 %). Ele mantém-se como segunda variante da imagem ideal, claramente separada, ao lado da opção por omissão de van der Waals; qual dos ramos é «mais correto» é questão da exatidão exigida, não da estrutura do modelo — o núcleo do modelo 3 (z̄ como parâmetro de ordem) permanece idêntico em ambos.

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

7.15 sólido(1)→sólido(2) noutros polimorfos: titânio e controlo com gás nobre

O cap. 7.13 calibra a linha (τ, p) sólido(1)→sólido(2) com ferro. A versão melhorada assegura sinal e ordem de grandeza num segundo polimorfo, magneticamente discreto, e acrescenta um controlo negativo (solidsolid_titanium_18.py). Calibração e salto de volume são como no cap. 7.13 (ε/k_B = T_m/τ_m, τ_m = 1,70; Δv = +8,9 % denso→aberto da geometria de empacotamento); o hcp é aproximado como o FCC por z = 12 (a diferença fcc/hcp é apenas nomeada, não resolvida).

O titânio, com α(hcp) → β(bcc) a 1155 K, tem exatamente a topologia do modelo: a fase densa (hcp, z = 12) é a de baixa temperatura e ao aquecer torna-se aberta (bcc, z = 8), logo ΔS > 0 e Δv > 0 → dp/dτ > 0 — o mesmo sinal do ferro γ→δ, mas sem anomalia magnética. Dois polimorfos denso→aberto independentes em comparação:

Ambos acertam o sinal (dp/dτ > 0, a pressão favorece a fase densa) e a ordem de grandeza (ΔH ~ kJ/mol, degrau de c_v sem parâmetros 1 R = 8,31 J/(mol·K), dT/dp ~ 10²·K/GPa). O valor de ΔH_ss fica em ambos os casos sistematicamente cerca de um fator demasiado alto — a aproximação de limiar abrupta da comparação dos μ sobrestima o calor latente de transformação; é um limite do modelo nomeado, não um erro de sinal.

O árgon é o controlo negativo: o árgon sólido é FCC e mantém-se FCC até à fusão (83,8 K) — à pressão normal não existe bcc estável. É precisamente o que prevê a imagem ideal do puro potencial de pares: θ_BCC/θ_FCC ≈ 0,9…1,0 > 0,61 e um BCC LJ instável ao corte (cap. 6.9) ⇒ nenhuma transição sólido(1)→sólido(2). A contraprova confirma a estrutura da constatação: a imagem ideal gera uma linha sólido(1)→sólido(2) apenas quando contributos não geométricos (eletrónicos/magnéticos, como no ferro e no titânio) estabilizam antes de mais o BCC. O titânio fornece a prova positiva, o árgon a negativa do mesmo quadro honesto dos cap. 6.8/6.9.

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

8 — Comparação bibliográfica: conceções da termodinâmica estática

Resumo. O capítulo 8 situa o modelo 3 entre as conceções estabelecidas da termodinâmica estática:

• 8.1 Tabela de conjunto: onze modelos (gás ideal, van der Waals, Redlich–Kwong, desenvolvimento do virial, Einstein, Debye, Lennard-Jones / esferas duras, teoria celular, Ising / gás em rede, percolação, Flory–Huggins) com a sua fase ou fases e a sua ideia de base estatístico-termodinâmica.
• 8.2 Os modelos um a um: um parágrafo cada com características, aplicação e limites, começando pelo elo para a fonte; as fontes de livre acesso estão arquivadas em model3/links/, sendo a teoria celular (pago) apenas citada.
• 8.3 Contraposição: o modelo 3 é um modelo de ligação cooperativa de dois estados com o parâmetro de ordem z̄; o parente mais próximo é Ising / gás em rede, a ebulição é retomada diretamente da percolação e a imagem dos fonões (Einstein/Debye) sustenta a rigidez de fusão. O modelo 3 agrafa estas três tradições numa escala ε comum.

Este capítulo situa o modelo 3 entre as imagens estabelecidas da termodinâmica de equilíbrio («estática») de sólidos, líquidos e gases. Primeiro um panorama tabelar (8.1), depois um breve parágrafo em prosa por modelo com características típicas, âmbito de aplicação e limites — cada parágrafo começa com o elo para a fonte (8.2). Por fim uma contraposição esquemática com o modelo 3 (8.3). As fontes de livre acesso estão arquivadas em model3/links/; o script de verificação model3/scripts/check_links_16.py mantém coerentes citação e arquivo por SHA-256.

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

8.1 Tabela de conjunto

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

8.2 Os modelos um a um

https://de.wikipedia.org/wiki/Kinetische_Gastheorie — Gás ideal / teoria cinética dos gases. Partículas sem volume próprio e sem interação exceto os choques elásticos; a pressão nasce puramente da transferência de momento. Fornece a equação de estado pV = Nk_BT, a distribuição de velocidades de Maxwell–Boltzmann e ⟨E⟩ = (3/2)k_BT. Aplicação: gases diluídos a alta temperatura e baixa pressão. Limites: sem condensação, sem fase líquida ou sólida, sem transição de fase — a interação está ausente por construção.

https://de.wikipedia.org/wiki/Van-der-Waals-Gleichung — Equação de van der Waals. (p + a/v²)(v − b) = k_BT completa o gás ideal com um volume de exclusão b (repulsão) e um termo de coesão a (atração). É a primeira equação de estado com uma transição gás-líquido, ponto crítico e construção de Maxwell. Aplicação: imagem qualitativa da condensação, grandezas críticas. Limites: quantitativamente imprecisa perto do ponto crítico (expoentes de campo médio), sem fase sólida.

https://de.wikipedia.org/wiki/Redlich-Kwong-Gleichung — Equação de estado de Redlich–Kwong. EOS cúbica com termo de atração dependente da temperatura a/[√T · v(v + b)]. O escalonamento em √T melhora nitidamente os equilíbrios vapor-líquido face a van der Waals (cf. cap. 7.4, onde isso eleva a razão T_b/T_c). Aplicação: cálculos técnicos de fase gasosa e evaporação. Limites: parametrizada empiricamente, ainda sem fase sólida; perto do ponto crítico igualmente do tipo campo médio.

https://de.wikipedia.org/wiki/Virialentwicklung — Desenvolvimento do virial. Desenvolvimento sistemático pV/Nk_BT = 1 + B(T)/v + C(T)/v² + … Os coeficientes do virial seguem do potencial intermolecular (integrais de aglomerado da mecânica estatística) e estão assim microscopicamente fundados. Aplicação: gases até densidade média, referência para outras EOS. Limites: a série não converge para líquidos densos ou sólidos; os coeficientes superiores são computacionalmente onerosos.

https://de.wikipedia.org/wiki/Einstein-Modell — Sólido de Einstein. O cristal é um conjunto de 3N osciladores harmónicos independentes de igual frequência ω. Explica a queda da capacidade térmica c_v abaixo da temperatura característica θ_E = ħω/k_B e o limite de Dulong–Petit 3Nk_B a alta temperatura. Aplicação: a descrição mais simples das capacidades térmicas dos sólidos. Limites: lei errada a baixa temperatura (exponencial em vez de T³), sem modos acústicos, apenas a fase sólida.

https://de.wikipedia.org/wiki/Debye-Modell — Modelo de Debye. Trata o sólido como contínuo elástico com um espetro de fonões até à frequência de Debye; c_v ∝ T³ a baixa temperatura, Dulong–Petit a alta — a temperatura de Debye θ_D resume a dureza da rede. Aplicação: capacidades térmicas quantitativamente boas de cristais simples em todo o intervalo de T. Limites: aproximação de contínuo (ignora a dispersão das redes reais), apenas a fase sólida, sem fusão.

https://de.wikipedia.org/wiki/Lennard-Jones-Potential — Potencial de Lennard-Jones / esferas duras. Potencial de pares 4ε[(σ/r)¹² − (σ/r)⁶] com caroço repulsivo e atração fraca; o modelo limite de esferas duras conserva apenas a repulsão. É a base da dinâmica molecular e do Monte Carlo e descreve — numericamente — as três fases juntamente com o ponto triplo. Aplicação: gases nobres, sistema de referência genérico. Limites: aditividade aos pares, avaliável no fundo apenas por simulação ou teoria de perturbações (sem EOS fechada).

https://doi.org/10.1098/rspa.1937.0210 — Teoria celular do líquido (Lennard-Jones–Devonshire). Cada molécula move-se numa célula (gaiola) formada pelos seus vizinhos em campo mediado; o líquido é descrito pela função de partição celular e fornece uma EOS do líquido bem como uma abordagem da fusão. Aplicação: líquidos simples, elo entre sólido e gás. Limites: sobrestima a ordem, negligencia lacunas/buracos (corrigem-no teorias de buracos posteriores). (Pago — apenas citada, não arquivada localmente; acessibilidade testada no script de verificação.)

https://de.wikipedia.org/wiki/Ising-Modell — Modelo de Ising / gás em rede. Numa rede os sítios portam uma variável de dois estados (spin ↑/↓ ou ocupado/vazio) com acoplamento a vizinhos próximos. O gás em rede é matematicamente isomorfo ao modelo de Ising e descreve a condensação gás↔líquido como transição cooperativa ordem-desordem com ponto crítico. Aplicação: fenómenos críticos, transições cooperativas. Limites: discretização em rede, sem energia cinética explícita; em aproximação de campo médio apenas qualitativo.

https://de.wikipedia.org/wiki/Perkolationstheorie — Teoria da percolação. Modelo puramente geométrico: sítios ou ligações são ocupados ao acaso com probabilidade p; acima de um limiar p_c existe um aglomerado que atravessa. Descreve a gelificação, o início da condutividade e, em geral, as transições de conectividade. Aplicação: formação de redes, desagregação/constituição da conexidade. Limites: sem energia nem temperatura intrínsecas — o acoplamento à termodinâmica tem de ser acrescentado do exterior.

https://de.wikipedia.org/wiki/Flory-Huggins-Theorie — Teoria de Flory–Huggins. Modelo em rede para soluções e polímeros: a entropia de mistura da colocação de cadeias mais um parâmetro de interação χ governam a miscibilidade e a separação de fases. Aplicação: soluções de polímeros, diagramas de fases de misturas. Limites: campo médio, rede incompressível, χ na maioria empírico — de relevância apenas limitada para imagens sólido/líquido/gás de um componente, aqui mantida como referência de rede.

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

8.3 Semelhanças e diferenças face ao modelo 3

O modelo 3 é um modelo de ligação cooperativa de dois estados: o parâmetro de ordem é o número de coordenação médio z̄ = z_max·x, autoconsistente pela cooperatividade K = J·z_max com a degenerescência do estado rompido G = z_max; a fusão como limiar de rigidez (fonões), a ebulição como desagregação da conectividade (percolação) no ponto crítico, transição de primeira ordem logo que ln G > 2 + 4ε/K.

Em suma: o modelo 3 agrafa três tradições de resto separadas — a imagem do gás em rede / de Ising (transição cooperativa de dois estados), a imagem da percolação (conectividade → ebulição) e a imagem dos fonões de Einstein/Debye (rigidez → fusão) — e calibra-as por um único parâmetro de energia ε numa escala τ comum. Van der Waals e a teoria celular fornecem a contraprova contínua (EOS e gaiola líquida), face à qual a imagem discreta das ligações pode ser comparada.

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

9 — Apêndice

Resumo. O apêndice reúne os comprovativos formais:

• 9.1 Bibliografia: ordenada por nome de autor, com elos para as fontes livremente arquivadas (comparação SHA-256 por check_links_16.py; a teoria celular, por causa do paywall, apenas citada).
• 9.2 Apoio da IA: a proporção de inovação IA/autor estimada pela IA para cada alteração substancial (com desagregação da parte do autor por estratégia/estrutura/redação e um cálculo de controlo), bem como uma lista dos passos individuais com custo (estimado), duração e distância ao registo anterior (soma ≈ 160 USD, estimativas claramente assinaladas).
• 9.3 Glossário de abreviaturas e símbolos: os símbolos de uso frequente (τ, ε, z̄, K, G, θ, …), que por isso podem surgir excecionalmente no texto corrido sem descrição prévia.
• 9.4 Scripts de cálculo e cobertura de testes: todos os números provêm de scripts Python reproduzíveis (proporção de programação por IA 100 %, cobertura de testes 0 % em toda a parte), juntamente com a aberta «declaração embaraçosa» de que o autor não verificou pessoalmente o código nem os resultados.

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

9.1 Bibliografia

Ordenada por nome de autor. As fontes de livre acesso estão arquivadas em model3/links/; a sua concordância com os elos citados é verificada por model3/scripts/check_links_16.py (SHA-256). A teoria celular de Lennard-Jones e Devonshire, por causa do paywall, é apenas citada (acessibilidade testada no script de verificação).

• Broadbent, S. R.; Hammersley, J. M. (1957): Percolation processes I. Crystals and mazes. Proc. Cambridge Phil. Soc. 53, 629–641. — Panorama/arquivo: Teoria da percolação
• Clausius, R.; Maxwell, J. C.; Boltzmann, L. (1857–1872): fundamentação da teoria cinética dos gases. — Panorama/arquivo: Teoria cinética dos gases
• Debye, P. (1912): Zur Theorie der spezifischen Wärmen. Ann. Phys. 39, 789–839. — Panorama/arquivo: Modelo de Debye
• Einstein, A. (1907): Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Ann. Phys. 22, 180–190. — Panorama/arquivo: Modelo de Einstein
• Flory, P. J.; Huggins, M. L. (1941–1942): teoria em rede das soluções de polímeros (parâmetro χ). — Panorama/arquivo: Teoria de Flory–Huggins
• Ising, E.; Lenz, W. (1920–1925): Beitrag zur Theorie des Ferromagnetismus. Z. Phys. 31, 253–258 (o gás em rede é matematicamente isomorfo). — Panorama/arquivo: Modelo de Ising
• Kamerlingh Onnes, H. (1901): forma do virial da equação de estado. — Panorama/arquivo: Desenvolvimento do virial
• Lennard-Jones, J. E. (1924): On the determination of molecular fields. Proc. R. Soc. Lond. A 106, 463–477. — Panorama/arquivo: 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 (Pago — apenas citado, não arquivado localmente.)
• Redlich, O.; Kwong, J. N. S. (1949): On the thermodynamics of solutions. V. An equation of state. Chem. Rev. 44, 233–244. — Panorama/arquivo: Equação de Redlich–Kwong
• van der Waals, J. D. (1873): Over de Continuïteit van den Gas- en Vloeistoftoestand (tese, Leida). — Panorama/arquivo: Equação de van der Waals

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

9.2 Apoio da IA

Este desenvolvimento de modelo nasceu de uma interação entre humano e IA: o autor (o utilizador da IA) determina em que se trabalha (perguntas orientadoras, decisões, desejos sobre o modelo); a IA (Claude Code) determina em larga medida como — derivações, scripts, estrutura e texto. As duas tabelas que se seguem são estimadas pela própria IA.

(a) Alterações substanciais de conteúdo e proporção de inovação estimada. A proporção de inovação pondera as ideias orientadoras a 45 %, os trabalhos estruturantes a 40 %, os trabalhos redatoriais a 15 %. As ideias orientadoras (o nível do em que) provêm predominantemente do autor, os trabalhos estruturantes e redatoriais (o nível do como) predominantemente da IA; onde a IA trouxe também a ideia física orientadora (por exemplo a leitura da ebulição por percolação, o modelo de molas para θ), a sua parte cresce. Valores em percentagem, IA + autor = 100. As três colunas adicionais desagregam a parte do autor segundo os três níveis de ação — estratégia (ideias orientadoras), estrutura (derivações/scripts/arquitetura) e redação (trabalho sobre o texto). A sua média ponderada a 45 % / 40 % / 15 % reproduz a coluna Autor % como cálculo de controlo (model3/scripts/author_share_19.py, desvio máximo 0,95 pontos percentuais).

(b) Passos individuais: custo (estimado), duração, distância ao registo anterior. Data de criação é a hora de criação do ficheiro de registo (mtime do ficheiro). Δ ao registo anterior é a distância temporal à entrada de registo anterior (medida do esforço do utilizador entre passos, incluindo pausas e /compact). Custo (estimado) são as estimativas documentadas por registo (preços de tabela da Anthropic, ± 30 %). Duração (ativa) é uma estimativa da IA do tempo de trabalho puro; em passos separados por /compact pode ser maior do que Δ. Um asterisco assinala uma presumível reconstrução ou reelaboração de conteúdos existentes.

Soma dos custos documentados (estimados, passos 01–21): ≈ 129 USD (soma das estimativas individuais); com o passo 22 (16d, ≈ 10 USD) ≈ 139 USD, o passo 23 (V16, ≈ 12 USD) ≈ 151 USD e o passo 24 (V17, ≈ 9 USD) ≈ 160 USD. Todos os valores são estimativas claramente assinaladas (não há medição ccusage disponível); os pormenores dos registos estão em model3/logs/dev_log_NN.md.

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

9.3 Glossário de abreviaturas e símbolos

Por princípio este texto antepõe a cada símbolo um termo descritivo (assim «capacidade térmica c_v» em vez de apenas «c_v»). As seguintes grandezas de uso frequente constituem a exceção: podem surgir no texto corrido sem descrição prévia, porque aqui ficam explicadas de uma vez por todas. Todas as temperaturas τ… são reduzidas (em unidades da energia de ligação ε).

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

9.4 Scripts de cálculo e cobertura de testes

Todos os números e tabelas mencionados no texto provêm de scripts Python reproduzíveis em model3/scripts/ (cada um com uma saída em bruto …_output.txt). A proporção de programação por IA indica que parte do código dos scripts foi escrita pela IA; a cobertura de testes refere-se a testes automatizados (testes unitários, verificações assert). É baixa em toda a parte: os scripts são cálculos únicos sem andaime de testes; a sua única garantia é a plausibilidade da saída em bruto.

O autor confessa abertamente (declaração embaraçosa): não verificou, recalculou nem validou tecnicamente o código-fonte destes scripts nem os números, tabelas e curvas com eles produzidos. Confiou aqui cegamente na competência da IA, porque em primeiro plano estavam apenas o caráter qualitativo e a plausibilidade da conceção do modelo. Os valores devem por isso ser tomados com cautela (ver também a nota detalhada no cap. 10.4).

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

10 Autocrítica

Tenho a sensação de que, na tabela do capítulo 9.2 (a), a parte de trabalho do autor estimada pela IA está, nas três áreas, quase sempre demasiado alta em pelo menos 10–20 pontos percentuais. Tenho a sensação de ter contribuído apenas com a ideia de base, uma ou duas interpretações e importantes refinamentos na visualização gráfica.

Os capítulos (matemáticos) 4–9 foram propostos, pesquisados e realizados pela IA. Pior ainda: durante o desenvolvimento tive a sensação de que a IA estava fixada na elaboração de uma estatística baseada no número de coordenação. As minhas verdadeiras ideias de base — que não só os átomos mas também as moléculas podem formar ligações e que as ruturas de ligação dentro de um líquido poderiam seguir regularidades fractais — a IA só as pôde acolher com dificuldade, tentando ignorar a terminologia. Por isso essas ideias de base continuam a ser tratadas de forma muito secundária no texto. Não o reescrevi depois, porque, ao refletir sobre as transições de fase a partir dos axiomas fundamentais, se têm naturalmente essas ideias presentes. Além disso, a fixação ou redução das ideias ao número de coordenação como grandeza de medida não é o pior caminho para matematizar a minha representação figurada, pois a posteriori se revelou bastante sustentadora para a verificação de plausibilidade — embora eu só tenha feito uma verificação de plausibilidade muito grosseira dos resultados da IA e nenhum exame do código Python. Em minha defesa diga-se: o que me importava foi sempre apenas uma representação figurada das relações moleculares nos diversos estados de agregação. À lisonja da IA, que queria confirmar matematicamente a representação, não disse que não. A manipulação por parte da IA teve, pois, êxito comigo, porque examinei os capítulos referidos apenas do modo mais superficial e os adotei com fé. (Esta experiência sobre mim próprio reforça a minha desconfiança quanto a saber se um «human in the loop» é realmente eficaz para controlar a IA.)

Este texto foi elaborado com apoio de IA (Claude Code 2.1.161, modelo Opus 4.8 / claude-opus-4-8[1m], contexto 1M; data 2026-09-05). — Elaborado com apoio de IA.

Nota sobre a marca de água. O texto mestre alemão (de.modell_final_02.md) porta nas palavras-chave uma marca de água feita de homóglifos gregos (Α→A, Β→B, Ε→E, Η→H, Ι→I, Κ→K, Μ→M, Ν→N, Ο→O, Ρ→P, Τ→T, Υ→Y, Χ→X, Ζ→Z). Como esse recurso está ligado à grafia alemã das palavras-chave, esta tradução porta as palavras-chave em maiúsculas normais: VOO LIVRE (FREIFLUG), MULTIDÃO (VIELHEIT), LIGAÇÃO (BINDUNG), TRÍADE (DREIKLANG), CAROÇO DURO (HARTKERN), CONSERVAÇÃO (ERHALTUNG), FISSURA (RISS), LESTO (FLINK).