PROSTE A POZYTECZNE NIERÓWNOśCI MACIERZOWE
TOMASZ TKOCZ
Streszczenie. Tekst zawiera notatki do referatu z Seminarium Fizyki Teoretycznej. Jest to w zasadzie tłumaczenie (z uzupełnionymi pewnymi skrótami myślowymi i innym podejściem do triku Weyla) odpowiedniego rozdziału książki [3].
1. Wstęp
W fizyce statystycznej kluczowe znaczenie ma suma statystyczna. W me- chanice kwantowej ma ona postać Z = tr e
−β( ˆH−µ ˆN )(β =
kT1, ˆ H — ha- miltonian układu, µ — potencjał chemiczny, ˆ N — operator liczby cząstek
— jeśli w układzie jest ona zachowana). Zatem bardzo ważne znaczenie jest umieć szacować wyrażenia tr e
A, gdzie A jest macierzą i niech to będzie dla nas motywacją, choć nie wiem jaka ona była naprawdę, żeby sobie dzisiaj obejrzeć kilka klasycznych nierówności. Przyjrzymy się nierównościom:
• Peierlsa – Bogolubowa
• H¨oldera
• Goldena – Thompsona
O ich zastosowaniach nie będzie tu wspomniane. Aby jednak wszystkich (pragmatyków?) uspokoić: nierówności te stosują się! Hasłowo, niektóre za- stosowania, to
• dowody wypukłości różnych ważnych funkcji termodynamicznych, np. energii swobodnej
• dowody istnienia granicy termodynamicznej
• oszacowania i porównywanie sumy statystycznej różnych modeli (np.
Isinga i Heisenberga — przykład z [3])
• dowód braku uporządkowania dalekiego zasięgu w modelu Hubbar- da (w wymiarach 1 i 2 i temperaturach dodatnich) — twierdzenie Komy–Tasakiego (por. [2])
• ostro też te nierówności stosuje Pan Wojtkiewicz podając oszacowa- nia (z góry i z dołu) sumy statystycznej w modelu Hubbarda (por.
[4])
Chciałbym podziękować Panu J. Wojtkiewiczowi — opiekunowi referatu — za pokaza- nie mi tytułowych nierówności.
1
2. Nierówność Peierlsa–Bogolubowa
Rozważać będziemy operatory liniowe A : C
n−→ C
n(o macierzy A) oraz na C
nmamy zadany iloczyn skalarny <, >. Przypomnijmy sobie kilka faktów z algebry liniowej.
Twierdzenie 1 (spektralne). Jeśli A jest samosprzężony (A
†= A), to istnieje baza ortonormalna u
1, . . . , u
noraz liczby rzeczywiste λ
1, . . . λ
ntakie, że
Ax =
Xnk=1
λ
k< x, u
k> u
k.
Definicja 1. Niech A — operator samosprzężony. Dla f : R → R definiuje- my operator f (A) : C
n→ C
nf (A)x :=
Xn
k=1
f(λ
k) < x, u
k> u
k.
Uwaga 1. Jeśli f rozwija się w otoczeniu 0 w szereg (f(x) =
P∞k=0 f(k)k!(0)x
k,
|x| ¬ R), to f(A) =
P∞k=0f(k)(0)k!
A
k, dla A takiego, że kAk ¬ R. Zatem po- dana ad hoc definicja f (A) dla pożądnych f pokrywa się z tą znaną skądinąd (np. dla funkcji e, ln, √ ).
Dowód. Rozwijamy f (λ
k) w szereg mając na uwadze, że kAk = max
1¬i¬n|λ
i|
¬ R
f (A)x =
Xnk=1
f(λ
k) < x, u
k> u
k=
Xnk=1
X∞
l=0
f
(l)(0)
l! λ
lk< x, u
k> u
k=
X∞l=0
f
(l)(0) l!
Xn
k=1
λ
lk< x, u
k> u
k=
X∞l=0
f
(l)(0) l! A
l.
Koniec przypomnienia. Przechodzimy do sedna tego paragrafu.
Twierdzenie 2. Niech f : R → R będzie funkcją wypukłą, A operatorem samosprzężonym. Wówczas dla dowolnego x, kxk = 1 zachodzi
f (< Ax, x >) ¬<f(A)x, x> . Dowód. Mamy wobec rozkładu spektralnego dla A, że
f (< Ax, x > = f <
Xn
k=1
λ
k< x, u
k>
!
= f
Xnk=1
| <x,u
k> |
2λ
k!
¬
Xnk=1
| <x, u
k> |
2f (λ
k) =< f (A)x, x >,
gdzie skorzystaliśmy z nierówności Jensena (dla wag | < x,u
k> |
2sumują-
cych się do kxk
2= 1).
Z powyższego mamy natychmiast
Wniosek 1. Jeśli (v
k) jest bazą ortonormalną, to tr f (A)
Xn
k=1
f(< Av
k, v
k>).
Dowód. Wystarczy zauważyć, że ogólnie < Av
k, v
k> jest współrzędną (k,k) macierzy operatora A w bazie (v
k), więc
tr f (A) =
Xnk=1
< f(A)v
k, v
k>
Xnk=1
f (< Av
k, v
k>).
Twierdzenie 3 (Nierówność Peierlsa–Bogolubowa). Dla operatorów samo- sprzężonych A, B i funkcji wypukłej f : R → R zachodzi
tr f (A)e
Btr e
B f tr Ae
Btr e
B!
.
Dowód. Niech (v
k) będzie bazą orotonormalną z twierdzenia spektralnego dla B, zaś µ
1, . . . , µ
kodpowiadającymi im wartościami własnymi. Mamy wobec nierówności z twierdzenia 2 i nierówności Jensena, że
tr f (A)e
Btr e
B= 1
P
e
µk Xnk=1
< f(A)e
Bv
k, v
k>
=
Xnk=1
e
µkP
e
µk< f(A)v
k, v
k>
Xnk=1
e
µkP
e
µkf (< Av
k, v
k>)
f
Xnk=1
e
µkP
e
µk< Av
k, v
k>
!
= f tr Ae
Btr e
B!
.
Identycznie dowodzimy (ale nie wynika ono z powyższego twierdzenia, bo niekoniecznie e
A+B= e
Ae
B)
Twierdzenie 4 (Nierówność Peierlsa–Bogolubowa’). Dla operatorów samo- sprzężonych A, B zachodzi
tr e
A+Btr e
B exp tr (Ae
B) tr e
B!
.
Dowód. Z poprzedniego dowodu widać, że tr e
Bjest tylko czynnikiem nor- mującym, więc można bez utraty ogólności założyć, że tr e
B= 1 (w razie czego rozważając zamiast B operator B + cI, dla odpowiedniej stałej c).
3
Wówczas ((v
k), (µ
k) ma to samo znaczenie co powyżej) tr e
A+B=
Xn
k=1
< e
A+Bv
k, v
k>
Xnk=1
e
<(A+B)vk,vk>=
Xnk=1
e
µke
<Avk,vk> exp
Xnk=1
e
µk< Av
k, v
k>
!
= exp
tr (Ae
B)
.
Jako ciekawostkę zauważmy na zakończenie tego paragrafu takie przyjem- ne zastosowanie głównego narzędzia — twierdzenia 2
Wniosek 2. Niech C : C
n→ C
nbędzie operatorem dodatnio określonym i samosprzężonym. Wówczas
det C ¬
Ynk=1
c
kk,
gdzie [c
ij] jest macierzą C w bazie standardowej (e
k).
Dowód. Dodatnia określoność oznacza, że dla każdego x 6= 0 < Cx, x >> 0, więc w szczególności wartości własne C są liczbami rzeczywistymi dodatnimi oraz c
kk> 0. Biorąc f (x) = − ln x w twierdzeniu 2 mamy
− <ln Ce
k, e
k>=< f (C)e
k, e
k> f(<Ce
k, e
k> ) = − ln c
kk, więc ln c
kk<ln Ce
k, e
k> i powysumowaniu
ln
Ynk=1
c
kk tr ln C =
Xnk=1
ln λ
k= ln
Ynk=1
λ
k= ln det C,
czyli teza.
3. Nierówność H¨ oldera
Będziemy chcieli odejść na chwilę od operatorów samosprzężonych. Ale w ogólności jest też coś na kształt twierdzenia spektralnego (gdy dim V = ∞ i operator jest zwarty nazywa się to rozkładem Schaudera). Niech więc A będzie teraz dowolnym endomorfizmem przestrzeni C
n. Bierzemy T := A
†A.
Wtedy T jest oczywiście samosprzężony, ale też nieujemny, bo
< T x, x > =< Ax, Ax >= kAxk
2 0
(widać po co A
†A zamiast AA
†). Zatem z twierdzenia spektralnego istnieje baza ortonormalna (u
k) i wartości własne nieujemne (λ
k) dla T . Możemy więc zdefiniwać
|A| := √ T = √
A
†A.
Mamy
k |A| xk
2= < |A|x, |A|x>=<|A|
†|A|x, x>=<|A|
2x, x >
(1)
= < A
†Ax, x > =< Ax, Ax >= kAxk
2,
czyli definiując U : im|A| → imA wzorem U(|A|x) := Ax mamy poprawnie określoną (bo ker |A| = ker A — wobec (1)) izometrię na swój obraz, którą łatwo rozszerzamy do izometrii C
n. Zatem
Ax = U (|A|x) = U
Xnk=1
λ
k< x, u
k> u
k!
=
Xnk=1
λ
k< x, u
k> U u
k=
Xnk=1
λ
k< x, u
k> v
k,
gdzie (v
k) pewna baza ortonormalna (obraz bazy (u
k) przy U ). Zatem za- chodzi
Wniosek 3. Dla A: C
n→ C
nistnieją bazy ortonormalne (u
k), (v
k) i liczby µ
1(A) . . . µ
n(A) 0 — wartości własne |A| zwane wartościami osobliwymi A takie, że
Ax =
Xnk=1
µ
k(A) < x, u
k> v
kUwaga 2. µ
k(A
†) = µ
k(A)
Dowód. Wobec U |A| = A jest A
†= |A|U
†i mamy
AA
†v
k= A|A|U
†v
k= A|A|u
k= Aµ
k(A)u
k= µ
k(A)
2v
k,
więc |A
†|v
k= µ
k(A)v
k, co oznacza, że |A|, |A
†mają te same wartości własne.
Definicja 2. kAk
p:= (tr |A|
p)
1/p, p ∈ [1, ∞].
Uwaga 3.
(tr |A|
p)
1/p=
Xnk=1
µ
k(A)
p!1/p
= k(µ
1(A), . . . , µ
n(A))k
p−−−→
p→∞max
k
µ
k(A)
= µ
1(A) = k |Ak | = kAk, czyli kAk
∞= kAk.
Możemy już teraz sformułować główne twierdzenie tego rozdziału
Twierdzenie 5 (Nierówność H¨oldera). Dla r, p, q ∈ [1, ∞] takich, że
1r=
1
p
+
1qi dowolnych operatorów A, B : C
n→ C
nzachodzi kABk
r¬ kAk
pkBk
q5
Dowód można łatwo przeprowadzić w oparciu o zwykłą (dla p-tych norm ciągów) nierówność H¨ oldera jeśli tylko udowodnimy, że
(2)
Xn
k=1
µ
k(AB)
r¬
Xnk=1
µ
k(A)
rµ
k(B)
r, bowiem wówczas stąd i uwagi 3 mamy
kABk
rr=
Xnk=1
µ
k(AB)
r¬
Xnk=1
µ
k(A)
rµ
k(B)
r¬
Xnk=1
µ
k(A)
r)
p/r!r/p n
X
k=1
µ
k(B)
r)
p/r!r/p
= kAk
rpkBk
rq. Dowód (2) zaś wynika z dwóch rzeczy
Lemat 1. Dla liczb rzeczywistych 0 ¬ a
1¬ . . . ¬ a
n, 0 ¬ b
1¬ . . . ¬ b
n, takich, że dla każdego k = 1, . . . , n
Qkj=1b
j¬
Qkj=1a
jzachodzi
Xk
j=1
b
j¬
Xnj=1
a
j, dla każdego k = 1, . . . , n.
Dowód. Jest to w zasadzie zasada majoryzacyjna Hardy’ego – Littlewooda – Póly zastosowana do funkcji exp. Dowód można znaleźć w [1] lub (bardziej
geometryczny) w [3].
Lemat 2. Dla dowolnych operatorów A, B : C
n→ C
ni każdego k = 1, . . . , n zachodzi
Yk
j=1
µ
j(AB) ¬
Ykj=1
µ
j(A)µ
j(B).
Dowód. Użyjemy sztuczki Weyla z algebrą przekształceń antysymetrycz- nych, po to żeby zgrabnie wydobyć iloczyn
Qkj=1µ
j(A). W przestrzeni od- wzorowań k-liniowych antysymetrycznych
Vk(C
n)
∗definiujemy dla opera- tora A : C
n→ C
njego k-tą potęgę dziubkową A
∧k:
Vk(C
n)
∗→
Vk(C
n)
∗wzorem
(A
∧kω)(v
1, . . . , v
k) := ω(Av
1, . . . , Av
k),
dla ϕ ∈
Vk(C
n)
∗, v
1, . . . , v
k∈ C
n. (Można o tym myśleć np. jak o pool–backu stałej formy ω przy przekształceniu A). Zauważmy też, że dla k = 1 ta potęga dziubkowa operatora A staje się jego zwykłym sprzężeniem Banachowskim.
Po co to wszystko? Mamy
µ
1(A
∧k) = największa wartość własna|A
∧k| (3)
(5)
= największa wartość własna|A
†|
∧k(4)
= µ
1(A
†) · . . . · µ
k(A
†) = µ
1(A) · . . . · µ
k(A),
skąd
Ykj=1
µ
j(AB)
(3)= µ
1(AB)
∧k= µ
1(B
∧kA
∧k)
= kB
∧kA
∧kk ¬ kB
∧kkk
∧kk = µ
1(B
∧k)µ
1(A
∧k)
(3)
=
Ykj=1
µ
j(A)
Ykj=1
µ
j(B), czyli teza lematu. Pozostaje jeszcze udowdnić
(4) największa wartość własna|A|
∧k= µ
1(A) · . . . · µ
k(A).
(5) |A
∧k| = |A
†|
∧kDla dowodu (4) weźmy operator S samosprzężony nieujemny (takim jest zawsze |A|) o wartościach własnych λ
1 . . . λ
n 0 dla bazy ortonor- malnej (u
k) (tzn. Su
k= λ
ku
k). Zauważmy, że (pamiętamy, że jeśli (du
k) oznacza bazę dualną do (u
k), to du
i1∧ . . . ∧ du
ik, 1 ¬ i
1< . . . i
k¬ n jest bazą ortonormalną
Vk(C
n)
∗)
S
∧kdu
i1∧ . . . ∧ du
ik= (S
∗du
i1) ∧ . . . ∧ (S
∗du
ik), ale S
∗du
l= λ
ldu
l, więc
S
∧kdu
i1∧ . . . ∧ du
ik= λ
i1· . . . · λ
ikdu
i1∧ . . . ∧ du
ik,
czyli S
∧kma wektory własne du
i1∧ . . . ∧ du
iko wartościach własnych λ
i1· . . . · λ
ik. Ponieważ jest ich tyle co dim
Vk(C
n)
∗, więc są to wszystkie wektory własne, więc i wszystkie wartości własne, a największa, to λ
1· . . . · λ
ki to trzeba było tu pokazać.
(5) udowodnimy w dwóch małych krokach. Po pierwsze (A
∧k)
†= (A
†)
∧k, bo wystarczy sprawdzić, że
< A
∧kω, η > =< ω, A
∧kη >, ω, η ∈
^k(C
n)
∗,
ale tę równość wystarczy sprawdzić na bazie ortonormalnej ω = dx
I, η = dx
Jdualnej do bazy standardowej (e
k) (używamy dalej multiindeksów I = (i
1, . . . , i
k)); wówczas, skoro (det
I(v
1, . . . , v
k) oznacza wyznacznik macierzy powstałej przez wybranie wierszy i
1, . . . , i
kz kolum v
1, . . . v
k)
A
∧kdx
I=
XL
det
L(Ae
i1, . . . , Ae
ik)dx
L, to
<dx
I, (A
†)
∧kdx
J> = det
I(A
†e
j1, . . . , A
†e
jk) = det
I(A
Te
j1, . . . , A
Te
jk)
= det
J(Ae
i1, . . . , Ae
ik) =< A
∧kdx
I, dx
J> .
7
Po drugie
|A
∧k| =
q(A
∧k)
†(A
∧k) =
q(A
†)
∧kA
∧k=
q(AA
†)
∧k 4= ( √
AA
†)
∧k= |A
†|
∧k.
4. Nierówność Goldena–Thompsona
Celem tutaj jest dowód następującego
Twierdzenie 6 (Nierówność Goldena–Thompsona). Dla operatorów samo- sprzężonych A, B : C
n→ C
nzachodzi
tr e
A+B¬ tr (e
Ae
B).
Potrzebne będzie kilka lematów.
Lemat 3. Dla operatorów samosprzężonych A, B : C
n→ C
ndodatnich (tzn.
< Ax, x >> 0, dla każdego x 6= 0) zachodzi
kABk
psp¬ kA
pB
pk
s, p 1, s 2.
Dowód. Ustalmy p, s jak w założeniach i oznaczmy C := A
p, D := B
p, f (α) := kC
αD
αk
s/p— funkcja ciągła na (0, 1]. Sprawdźmy czy przedłuża się ona na [0, 1], tzn. chcemy obliczyć lim
α→0+f (α). Szacujemy
kC
αD
αk
s/αH¨older
¬ kC
αk
2s/αkD
αk
2s/α= (tr Ctr D)
α/2s−−−−→
α→0+
1,
dalej, biorąc za v dowolny jednostkowy wektor własny dla D o wartości własnej µ mamy z drugiej strony
kC
αD
αk
s/pporównanie p-tych norm ciągów
kC
αDα k
∞= kC
αD
αk
kC
αD
αv k =
Xn
k=1
µ
αλ
αk< v, u
k> u
k (|µ|min
1¬k¬n|λ
k|)
α−−−−→
α→0+
1, gdzie korzystaliśmy jeszcze z twierdzenia spektralnego dla C (oczywiście (u
k) oznacza bazę ortonormalna wektorów własnych odpowiadających war- tościom własnym (λ
k)). Zatem f (0) = 1. Zauważmy, że teza, to nierówność
f
1 p
¬ f(1) = f(1)f(0), co równoważnie można zapisać jako
f
1 p · 1 +
1 − 1
p
· 0
¬ f(1)
1/pf(0)
1−1/p,
czyli wystarczy udowodnić wypukłość funkcji ln f . Wobec jej ciągłości wy- niknie, to z nierówności
(6) f
α + β 2
¬ f(α)
1/2f (β)
1/2. Do jej dowodu przydadzą się trzy rzeczy.
(7) kXk
ppq= (tr |X|
pq)
1/q= k|X|
pk
q, p, q 1 X— dowolny
(8) kXY k
p¬ kY Xk
p, X, Y — dowolne takie, że XY samosprzężony, bo wobec lematu 1 (zasady majoryzacyjnej) i definicji p-tej normy macie- rzowej wystarczy pokazać, że dla każdego k = 1, . . . , n
Yk
j=1
µ
jXY ¬
Ykj=1
µ
jY X.
Jednak mamy
Ykj=1
µ
j(XY ) = µ
1(XY )
∧k) = max modułów wartości własnychY
∧kX
∧k= r(Y
∧kX
∧k) = r(X
∧kY
∧k) = r((Y X)
∧k)
¬ k(Y X)
∧kk = µ
1(Y X)
∧k=
Ykj=1
µ
jY X,
gdzie przez r(X) oznaczamy promień spektralny operatora X i korzystamy ze znanych faktów, że ogólnie r(XY ) = r(Y X), bo XY i Y X mają te same wartości własne oraz r(X) ¬ kXk, bo jeśli u jest jednostkowym wektorem własnym o wartości własnej λ, to |λ| = kXuk ¬ kXk.
(9) kX
†k
p= kXk
p,
oczywiste, wobec uwagi 2.
Wracamy do dowodu (6). Bez utraty ogólności załóżmy, że β α. Mamy f
α + β 2
= kC
α+β2D
α+β2k
2s/(α+β)(7)
= k|C
α+β2D
α+β2|
2k
s/(α+β)= kD
α+β2C
α+βC
α+β2k
1/2s/(α+β)= kD
αD β − α
2 C
α+βD
α+β2k
(8)
= kD
αC
α+βD
βk
1/2s/(α+β)H¨older¬ kD
αC
αk
1/2s/αkC
βD
βk
1/2s/β(9)
= f (α)
1/2f (β)
1/2.
Przy okazji można stąd udowodnić taki przyjemny
9
Wniosek 4. Przy założeniach lematu 3 zachodzi także nierówność tr (AB)
n¬ tr (A
nB
n).
Dowód. Korzystając z cykliczności śladu otrzymujemy
tr (AB)
n= tr
A
1/2(A
1/2BA
1/2) · . . . · (A
1/2BA
1/2)
| {z }
n−1
A
1/2B
= tr (A
1/2BA
1/2) = tr |B
1/2A
1/2|
2n= kB
1/2A
1/2k
2n2n lemat 3¬ kB
n/2A
n/2k
22= tr (A
n/2B
nA
n/2) = tr (A
nB
n).
Lemat 4. Dla dowolnych operatorów X, Y : C
n→ C
nzachodzi
e
X+Y= lim
n→∞
(e
n1Xe
1nY)
n= lim
n→∞
(e
2n1 Xe
n1Ye
2n1 X)
n, gdzie zbieżność jest w normie operatorowej.
Dowód. Niech A
n:= e
n1Xe
n1Y, B
n:= e
n1(X+Y ). Chemy udowodnić, że kA
nn− B
nnk −−−→
n→∞0, ale
kA
nk, kB
nk ¬ e
1n(kXk+kY k)kA
n− B
nk ¬ C
n
2,
bo pierwsze wyrazy odpowiednich szeregów się redukują i stała C zależy tylko od operatorów X, Y (nie zależy od n). Mamy stąd
kA
nn− B
nnk = kA
n−1n(A
n− B
n) + A
n−2n(A
n− B
n)B
n+ A
n−3n(A
n− B
n)B
n2+ . . . + A
n(A
n− B
n)B
nn−2+ (A
n− B
n)B
nn−1k
¬ ne
kXk+kY kkA
n− B
nk ¬ C
n e
kXk+kY k−−−→
n→∞0.
Dowód drugiej równości jest identyczny.
Lemat 5. Dla operatorów samosprzężonych A, B : C
n→ C
nzachodzi ke
A+Bk
p¬ ke
Ae
Bk
p, p 1.
Dowód. Niech najpierw p 2; wtedy ke
2n1 Ae
2n1 Bk
2n2nplemat 3
¬ ke
Ae
Bk
p, ale lewa strona wynosi
k|e
2n1 Ae
2n1 B|
2nk
p= k(e
2n1 Be
n1Ae
2n1B)
nk
p−−−−→
lemat 4n→∞ke
(A + B)k,
skąd teza w tym przypadku.
Jeśli 1 ¬ p < 2, to mamy z poprzedniego przypadku ke
A+Bk
p= ke
A+B2k
22p¬ ke
12Ae
12Bk
22p= k|e
12Ae
12B|
2k
p= ke
12Be
Ae
12Bk
p(8)
¬ ke
Ae
Bk
p.
Teraz łatwo zakończyć dowód wyjściowego twierdzenia.
Dowód twierdzenia 6. Pamiętając, że założenia A + B jest też samosprzężo- ny, mamy
tr e
A+B= ke
12(A+B)k
22 lemat 5¬ ke
12Ae
12Bk
22= tr
e
12Be
Ae
12B= tr (e
Ae
B).
Jako zastosowanie wszystkiego podamy jeszcze bardzo ważny
Wniosek 5. Funkcja f : {operatory samosprzężone C
n→ C
n} −→ R zada- na wzorem f (A) := ln
tr e
Ajest wypukła.
Dowód. Po prostu sprawdzamy definicję wypukłości
f (αA + (1 − α)B) = ln
tr e
αA+(1−α)BGolden–Thompson
¬ ln
tr (e
αAe
(1−α)B)
H¨older
¬ ln
(tr e
A)
α(tr e
B)
1−α= αf (A) + (1 − α)f(B).
Literatura
[1] M. Kuczma, An introduction to the theory of functional equations and inequalities.
Cauchy’s equation and Jensen’s inequality, Prace naukowe Uniwersytetu Śląskiego w Katowicach, nr 489, Państwowe Wydawnictwo Naukowe i Uniwersytet Śląski, Warsza- wa Kraków Katowice 1985
[2] T. Koma, H. Tasaki, Decay of Superconducting and Magnetic Correlations in One and Two–Dimensional Hubbard Models, Phys. Rev. Lett. 68, 2348(1992)
[3] B. Simon, The tatistical Mechanics of Lattice Gases, Princeton, New Jersey, Princeton University Press
[4] J. Wojtkiewicz, Estimations of free energy for the Hubbard model
11