In 1935, Pauling proposed an estimate for the number of Eulerian orientations of a graph in the context of the theoretical behaviour of water ice. The logarithm of the number of Eulerian orientations, normalised by the number of vertices, is called the residual entropy. In an earlier paper, we conjectured that the residual entropy of a sequence of regular graphs of increasing degree was asymptotically equal to Pauling's estimate. Here we prove the conjecture under constraints on the number of short circuits. These constraints hold under weak eigenvalue conditions and apply to sequences of increasing girth and repeated Cartesian products such as hypercubes.
- ID articolo: 2509.20671
- Titolo: On Pauling's residual entropy estimate for regular graphs with growing degree
- Autori: Mahdieh Hasheminezhad (Yazd University), Mikhail Isaev (UNSW Sydney), Brendan D. McKay (Australian National University), Rui-Ray Zhang
- Classificazione: math.CO (Matematica combinatoria)
- Data di pubblicazione: 5 novembre 2025 (arXiv v2)
- Link articolo: https://arxiv.org/abs/2509.20671
Nel 1935, Pauling propose una stima del numero di orientamenti euleriani di grafi durante lo studio del comportamento teorico del ghiaccio d'acqua. Il logaritmo del numero di orientamenti euleriani diviso per il numero di vertici è denominato entropia residua. In articoli precedenti, gli autori hanno congetturato che l'entropia residua di sequenze di grafi regolari con grado crescente sia asintoticamente uguale alla stima di Pauling. Questo articolo dimostra tale congettura sotto vincoli sul numero di cicli brevi. Questi vincoli valgono sotto deboli condizioni di autovalori e si applicano a sequenze di grafi con circonferenza crescente e prodotti cartesiani ripetuti (come gli ipercubi).
Questo articolo studia il problema del conteggio degli orientamenti euleriani (Eulerian orientations) di grafi. Per un grafo d-regolare G (d pari), un orientamento euleriano è un orientamento degli spigoli tale che ogni vertice abbia grado entrante uguale al grado uscente. L'entropia residua è definita come:
ρ(G):=n1logEO(G)
dove EO(G) è il numero di orientamenti euleriani e n è il numero di vertici.
- Sfondo fisico: L'entropia residua ha importanza significativa nei modelli di ghiaccio (ice-type models) della fisica statistica, correlata al numero di microstati dei cristalli di ghiaccio d'acqua
- Significato matematico: Per reticoli specifici (reticolo quadrato, reticolo triangolare, monostrato di ghiaccio esagonale), il valore esatto dell'entropia residua è noto, ma il comportamento asintotico per grafi generali rimane un problema aperto
- Valore teorico: Connette la matematica combinatoria, la fisica statistica e la teoria della probabilità
Pauling nel 1935 propose una stima euristica: assumendo che ogni spigolo sia orientato casualmente in modo indipendente, la probabilità che il vertice i abbia grado entrante uguale al grado uscente è 2−d(d/2d). Se gli eventi per n vertici fossero indipendenti, si otterrebbe:
ρ^(G)=log(d/2d)−2dlog2
Lieb e Wu hanno dimostrato che per qualsiasi grafo: ρ(G)≥ρ^(G)
- I risultati precedenti (Isaev, McKay, Zhang 5) valgono solo per grafi con buone proprietà di espansione e d≥log8n
- I risultati su grafi casuali (6) dipendono da proprietà ad alta probabilità
- Il risultato di Las Vergnas 7 richiede che la circonferenza cresca più velocemente di logd
Congettura 1.1: Per una sequenza di grafi d-regolari G(n), se l'intero pari d=d(n)→∞, allora ρ(G)=ρ^(G)+o(1)
L'obiettivo di questo articolo è provare tale congettura sotto condizioni più deboli.
- Teorema principale (Theorem 2.1): Dimostra la Congettura 1.1 sotto vincoli sul numero di tracce chiuse brevi, con condizioni:
- Esiste ℓmax=ω(logd) e una costante C > 0
- Per tutti 3≤ℓ≤ℓmax: cℓ(G)≤Ce−(l+1)dℓ−1n
- Corollario su condizioni di autovalori (Corollary 2.2): Se al massimo nd−ω(1) autovalori si trovano al di fuori di [−d1−δ,d1−δ] (per una costante δ > 0), allora la congettura vale
- Corollario su condizioni di circonferenza (Corollary 2.3): Per sequenze di grafi d-regolari con circonferenza g→∞, la congettura vale (non richiede d→∞)
- Corollario su prodotti cartesiani (Corollary 2.4): Per prodotti cartesiani Gt=H1(t)□⋯□Ht(t), la congettura vale quando la somma dei quadrati dei gradi soddisfa ∑i=1t(hi(t))2=O(dt2−δ), in particolare per ipercubi Qd
- Innovazione tecnica: Trasforma il problema in un'analisi della funzione generatrice dei momenti del numero di tracce chiuse indotte da partizioni euleriane casuali
Dato un grafo d-regolare G (d pari), calcolare l'entropia residua ρ(G)=n1logEO(G) e provare che è asintoticamente uguale alla stima di Pauling ρ^(G).
Definizione: Una partizione euleriana P è una partizione degli spigoli del grafo in diverse tracce chiuse. L'accoppiamento di spigoli di ogni vertice determina univocamente una partizione euleriana.
Formula chiave:
ρ(G)=ρ^(G)+n1logE2∣T(P)∣
dove P è una partizione euleriana uniformemente casuale e T(P) è l'insieme di tracce chiuse indotte da P.
Equivalenza: La congettura è equivalente a provare E2∣T(P)∣=eo(n)
Definire k=⌊min{ℓmax/2,log2d}⌋, dividere le tracce chiuse in:
- Tracce lunghe Lk(P): numero di tracce chiuse contenenti più di k vertici distinti
- Tracce brevi Sk(P): numero di tracce chiuse contenenti al massimo k vertici distinti
Utilizzare la disuguaglianza di Hölder:
E2∣T(P)∣=E2Lk(P)+Sk(P)≤(E227Lk(P))2/7(E257Sk(P))5/7
Lemma chiave (Lemma 3.1): Per m pari e λ ≥ 0,
EeλX(m)≤(Bm)(eλ−1)/2
dove X(m) è la somma di variabili casuali di Bernoulli indipendenti con parametri 1/(2j−1) (j=1,...,m/2).
Indipendenza degli accoppiamenti (Lemma 3.2): Negli accoppiamenti casuali, il numero di tracce chiuse distinte attraverso il vertice v segue la distribuzione X(m) ed è indipendente dagli accoppiamenti di altri vertici.
Limite della funzione generatrice dei momenti (Lemma 3.3): Per un sottoinsieme di vertici U,
EeλY(U)=dO(∣U∣)
Limite finale: Scegliere un sottoinsieme casuale di ∣U∣=⌊n/k⌋, utilizzare il fatto che Lk≤2Y(U) con probabilità almeno 1/2, ottenere:
EeλLk≤(edk)O(n/k)=eo(n)
Utilizzare il metodo di switching (Theorem 4.1) — uno strumento potente di conteggio combinatorio.
Definizione di switching: Data una partizione euleriana P e una traccia chiusa T, uno T-switching modifica l'accoppiamento in ogni vertice attraversato da T:
- Il vertice v è visitato t volte da T, corrispondente alle coppie di spigoli (e1,e1′),…,(et,et′)
- Selezionare t coppie di spigoli aggiuntive (et+1,et+1′),…,(e2t,e2t′)
- Sostituire con il nuovo accoppiamento {e1,et+1},…,{et,e2t} e {e1′,et+1′},…,{et′,e2t′}
Grafo di switching: Costruire un multigrafo orientato Γ:
- Vertici: m=(m3,…,mL), rappresentando classi di partizioni euleriane con esattamente mℓ tracce chiuse di lunghezza ℓ
- Spigoli: uno spigolo orientato colorato ℓ da (m,m′) rappresenta l'esistenza di uno ℓ-switching da C(m) a C(m')
Stima chiave (Lemma 4.3):
- Numero di ℓ-switching da P ∈ C(m): ≥∥m∥1(d−2L)ℓ/L
- Numero di ℓ-switching verso P' ∈ C(m'): ≤ck,ℓ
Funzione di peso:
α^(e)≤dℓm2ck,ℓL2
Applicazione di Theorem 4.1: Definire
M0:=2CnL2/d,Y=⋃m>MCm,Z=⋃m≤M0Cm
Poiché α^(e)<1/2 per tutti gli spigoli con ∥m∥1>M0, ottenere:
∣Y∣≤2e−M+M0∣Z∣
Limite di probabilità della coda:
P(Sk(P)>M)≤2e−M+M0
Limite dei momenti: Attraverso la rappresentazione integrale,
EeλSk(P)≤eλM0+1−λ2eλM0=eo(n)
dove λ=57log2≈0.97<1.
- Eleganza della strategia di decomposizione: Attraverso la scelta di k si bilancia il contributo delle tracce lunghe e brevi, la scelta degli esponenti della disuguaglianza di Hölder (7/2 e 7/5) è attentamente progettata
- Applicazione innovativa del metodo di switching:
- Trasforma il problema di conteggio delle partizioni euleriane in un problema di flusso su grafi orientati
- Utilizza il vincolo sul numero di tracce brevi per controllare i pesi di switching
- Attraverso Theorem 4.1 ottiene il decadimento esponenziale della probabilità della coda
- Indebolimento delle condizioni:
- Da d≥log8n a d→∞
- Da forte espansività a vincoli sul numero di tracce brevi
- La condizione di autovalori richiede solo nd−ω(1) outlier
- Struttura unificata: Tratta simultaneamente circonferenza crescente, vincoli di autovalori, prodotti cartesiani e altre situazioni
Questo articolo è un articolo puramente teorico e non contiene esperimenti numerici o computazionali. Tutti i risultati sono dimostrazioni matematiche rigorose.
L'articolo verifica la congettura attraverso corollari su le seguenti classi di grafi:
- Grafi con circonferenza crescente (Corollary 2.3)
- Grafi con vincoli spettrali (Corollary 2.2)
- Ipercubi Qd (d pari)
- Prodotti cartesiani iterati: convergono localmente a reticoli di dimensione superiore
- Grafi regolari casuali (risultati citati da 6)
Verifica delle condizioni di Theorem 2.1:
- Grafi con circonferenza g → ∞:
- Il numero di tracce chiuse di lunghezza ℓ < g è zero, soddisfa automaticamente le condizioni
- Vale anche per d = O(1)
- Grafi con vincoli di autovalori:
- Il numero di cammini chiusi ≤∑iλiℓ
- Se al massimo nd−f(n) autovalori si trovano al di fuori di [−d1−δ,d1−δ]
- Scegliendo ℓmax=21f(n)logd, le condizioni sono soddisfatte
- Prodotti cartesiani:
- Utilizzare la disuguaglianza di Hoeffding per controllare la distribuzione degli autovalori
- Per ∑i(hi(t))2=O(dt2−δ), le condizioni sono soddisfatte
- Ipercubi: hi=1, soddisfa ∑hi2=t=o(dt2)
Limite delle tracce lunghe (Lemma 3.4):
EeλLk(P)=eo(n),∀λ≥0,k≫logd
Limite delle tracce brevi (Lemma 4.5):
E257Sk(P)=eo(n)
\mathbb{E} 2^{|T(P)|} &\leq \left(\mathbb{E} 2^{\frac{7}{2}L_k(P)}\right)^{2/7}\left(\mathbb{E} 2^{\frac{7}{5}S_k(P)}\right)^{5/7}\\
&= (e^{o(n)})^{2/7} \cdot (e^{o(n)})^{5/7}\\
&= e^{o(n)}
\end{align}$$
Pertanto $\rho(G) = \hat{\rho}(G) + o(1)$.
## Lavori correlati
### 1. Contesto storico
- **Pauling (1935)**: Propone la stima euristica dell'entropia residua, utilizzata per spiegare l'entropia zero del ghiaccio d'acqua
- **Lieb & Wu (1972)**: Dimostrano il limite inferiore $\rho(G) \geq \hat{\rho}(G)$
### 2. Risultati esatti
- **Lieb (1967)**: Valore esatto per il reticolo quadrato
- **Baxter (1969)**: Valore esatto per il reticolo triangolare
- **Li et al. (2023)**: Valore esatto per il monostrato di ghiaccio esagonale
### 3. Risultati asintotici
- **Las Vergnas (1990)**: La congettura vale quando la circonferenza cresce più velocemente di $\log d$
- **Isaev, McKay, Zhang (2025, [5])**: La congettura vale per grafi espandenti con $d \geq \log^8 n$
- **Isaev, McKay, Zhang (2025, [6])**: Vale ad alta probabilità per grafi casuali
### 4. Metodologia
- **Metodo di switching**: Greenhill & McKay (2013), Hasheminezhad & McKay (2010)
- **Espansione dei cumulanti**: Metodo in [5]
- **Accoppiamenti probabilistici**: Lugo (2009) sulla struttura ciclica di involuzione casuale
### 5. Vantaggi relativi di questo articolo
- **Condizioni più deboli**: Richiede solo vincoli sul numero di tracce brevi
- **Applicazioni più ampie**: Copre circonferenza, autovalori, prodotti cartesiani e altre situazioni
- **Struttura tecnica unificata**: Un unico framework per trattare diverse classi di grafi
## Conclusioni e discussione
### Conclusioni principali
1. **Livello teorematico**: Sotto il vincolo sul numero di tracce chiuse brevi $c_\ell(G) \leq Ce^{-(l+1)}d^{\ell-1}n$ ($3 \leq \ell \leq \ell_{max} = \omega(\log d)$), la congettura di Pauling vale
2. **Livello dei corollari**:
- Sequenze di grafi con circonferenza crescente soddisfano la congettura
- Grafi con vincoli spettrali soddisfano la congettura
- Prodotti cartesiani (inclusi gli ipercubi) soddisfano la congettura
3. **Livello metodologico**: Il metodo di switching combinato con l'analisi della funzione generatrice dei momenti è uno strumento efficace per affrontare questo tipo di problemi
### Limitazioni
1. **Dipendenza dalle condizioni**:
- Richiede ancora $d \to \infty$ (eccetto nel caso di circonferenza)
- Il vincolo sul numero di tracce brevi, sebbene debole, non è banale
- Il requisito $\ell_{max} = \omega(\log d)$
2. **Limitazioni tecniche**:
- La scelta degli esponenti della disuguaglianza di Hölder (7/2, 7/5) sembra essere tecnica e potrebbe non essere ottimale
- La scelta di k dipende dal bilanciamento tra $\ell_{max}$ e $\log^2 d$
3. **Ambito di applicazione**:
- Non si applica a sequenze di grafi con grado limitato ($d = O(1)$ vale solo nel caso di circonferenza)
- L'estensione a grafi irregolari non è ovvia
4. **Complessità computazionale**:
- I risultati sono asintotici e non forniscono limiti di errore per n finito
- Non ci sono contributi a livello algoritmico
### Direzioni future
1. **Indebolimento delle condizioni**:
- È possibile eliminare completamente il vincolo sulle tracce brevi?
- È possibile trattare il caso generale di $d = O(1)$?
2. **Termini di errore**:
- Qual è l'ordine esatto del termine o(1)?
- È possibile ottenere una stima più fine $\rho(G) = \hat{\rho}(G) + O(1/d)$?
3. **Grafi irregolari**:
- Estensione a grafi con sequenze di gradi non completamente uniformi
- Trattamento del caso di grafi bipartiti
4. **Applicazioni algoritmiche**:
- È possibile progettare algoritmi di conteggio approssimato?
- Analisi del tempo di mescolanza di MCMC
5. **Connessioni fisiche**:
- Connessioni più profonde con la teoria delle transizioni di fase
- Applicazioni ad altri modelli di reticolo
## Valutazione approfondita
### Punti di forza
1. **Profondità teorica**:
- Tecniche di dimostrazione sofisticate, combinazione abile di strumenti da diversi campi
- L'applicazione del metodo di switching dimostra la potenza dell'ottimizzazione combinatoria
- L'analisi della funzione generatrice dei momenti è rigorosa e dettagliata
2. **Ampiezza dei risultati**:
- Trattamento unificato di diverse classi di grafi (circonferenza, autovalori, prodotti cartesiani)
- I corollari coprono applicazioni importanti (ipercubi, reticoli di dimensione superiore)
- Connessione tra matematica combinatoria e fisica statistica
3. **Innovazione tecnica**:
- La strategia di decomposizione delle tracce lunghe e brevi è innovativa
- La costruzione del grafo di switching è ingegnosa
- L'uso della disuguaglianza di Hölder è appropriato
4. **Qualità della scrittura**:
- Struttura chiara, logica rigorosa
- Dettagli tecnici completi, facili da verificare
- Introduzione al contesto sufficiente, motivazione chiara
### Insufficienze
1. **Condizioni non ottimali**:
- Il requisito $\ell_{max} = \omega(\log d)$ potrebbe essere indebolito
- La scelta degli esponenti di Hölder (7/2, 7/5) sembra avere spazio di ottimizzazione
2. **Limitazioni di applicazione**:
- Copertura insufficiente del caso di grafi con grado limitato (come d fisso)
- L'estensione a grafi irregolari non è ovvia
3. **Assenza di calcoli**:
- Nessuna verifica numerica o calcolo di esempi specifici
- I limiti di errore per grafi finiti sono sconosciuti
4. **Interpretazione fisica**:
- Sebbene abbia sfondo fisico, la discussione della connessione con transizioni di fase e fenomeni critici è insufficiente
- Il significato fisico dell'entropia residua potrebbe essere esplorato più profondamente
### Impatto
1. **Contributo accademico**:
- Risolve una congettura importante della matematica combinatoria (sotto condizioni di vincolo)
- Fornisce strumenti e framework potenti per ricerche successive
- Potrebbe stimolare ulteriore sviluppo del metodo di switching
2. **Valore pratico**:
- Fornisce supporto teorico ai modelli di ghiaccio della fisica statistica
- Ispira il problema di orientamento nei network science
- Potrebbe applicarsi alla progettazione di algoritmi di conteggio approssimato
3. **Riproducibilità**:
- La dimostrazione è completa e verificabile passo dopo passo
- Gli strumenti tecnici (Theorem 4.1) sono stabiliti in letteratura
- La verifica dei corollari è relativamente diretta
4. **Ricerca successiva**:
- Stimola l'ulteriore indebolimento delle condizioni
- Promuove la ricerca su grafi irregolari
- Potrebbe stimolare lavori su algoritmi e complessità computazionale
### Scenari applicabili
1. **Ricerca teorica**:
- Problemi di conteggio in matematica combinatoria
- Studio delle proprietà euleriane in teoria dei grafi
- Combinatoria probabilistica
2. **Applicazioni fisiche**:
- Meccanica statistica del modello di ghiaccio
- Funzione di partizione di sistemi reticolari
- Teoria delle transizioni di fase
3. **Network science**:
- Problemi di orientamento in reti su larga scala
- Analisi della proprietà di bilanciamento in reti di flusso
- Ricerca sulla reciprocità in reti sociali
4. **Progettazione di algoritmi**:
- Base teorica per algoritmi di conteggio approssimato
- Analisi del tempo di mescolanza di metodi MCMC
- Garanzie di prestazione di algoritmi casuali
## Bibliografia
### Citazioni chiave
1. **Pauling (1935)**: Proposizione della congettura originale, entropia residua dei cristalli di ghiaccio
2. **Lieb & Wu (1972)**: Dimostrazione del limite inferiore e studio sistematico del modello ferroelettrico
3. **Greenhill & McKay (2013)**: Sistematizzazione del metodo di switching
4. **Isaev, McKay, Zhang (2025, [5])**: Metodo di espansione dei cumulanti, risultato per $d \geq \log^8 n$
5. **Isaev, McKay, Zhang (2025, [6])**: Risultato per grafi casuali, connessione con l'entropia dell'albero generatore
6. **Las Vergnas (1990)**: Limite superiore sotto condizioni di circonferenza
7. **Lugo (2009)**: Struttura ciclica di involuzione casuale, indipendenza dell'accoppiamento
### Letteratura metodologica
- **Hoeffding (1963)**: Disuguaglianze probabilistiche
- **McKay (1981)**: Distribuzione attesa degli autovalori di grafi regolari
- **Merris (1998)**: Autovettori laplaciani di grafi, proprietà spettrali dei prodotti cartesiani
---
**Valutazione complessiva**: Questo è un articolo teorico di alta qualità che risolve parzialmente una congettura importante della matematica combinatoria e della fisica statistica sotto condizioni indebolite. Le tecniche di dimostrazione sono sofisticate, i risultati hanno ampia copertura e forniscono contributi importanti al campo. Sebbene le condizioni non siano ancora ottimali, l'articolo prepara il terreno per una risoluzione completa della congettura. L'articolo dimostra la potenza del metodo di switching e merita uno studio approfondito da parte dei ricercatori di matematica combinatoria.