SF - Zbirka, Vladimir Miljković
SF - Zbirka, Vladimir Miljković
SF - Zbirka, Vladimir Miljković
UNIVERZITET U BEOGRADU
FIZIČKI FAKULTET
Zbirka zadataka iz
Statističke fizike
verzija 3.48
§1 Matematički uvod . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
§2 Statistička termodinamika . . . . . . . . . . . . . . . . . . . . . . . . . . 5
I. Prvi i drugi princip termodinamike . . . . . . . . . . . . . . . . . . 5
a) Termomehanički sistemi . . . . . . . . . . . . . . . . . . . . . 6
b) Elastične trake . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
c) Termodinamika kondenzovanog stanja materije . . . . . . . . 39
d) Termodinamika magnetika . . . . . . . . . . . . . . . . . . . . 44
e) Termodinamika sistema sa promenljivim brojem čestica . . . 54
f) Termodinamičke nejednakosti . . . . . . . . . . . . . . . . . . 59
II. Treći princip termodinamike . . . . . . . . . . . . . . . . . . . . . . 63
§3 Fazni prelazi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
I. Clausius–Clapeironova jednačina . . . . . . . . . . . . . . . . . . . . 67
II. Osnove Landau-ove teorije faznih prelaza . . . . . . . . . . . . . . . 80
III. Hipoteza skaliranja za termodinamičke funkcije . . . . . . . . . . . 93
IV. Primeri faznih prelaza . . . . . . . . . . . . . . . . . . . . . . . . . . 99
§8 Dodatak . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 296
I. Gama funkcija . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 296
II. O Riemann-ovoj ζ funkciji . . . . . . . . . . . . . . . . . . . . . . . . 297
III. Bessel-ove funkcije . . . . . . . . . . . . . . . . . . . . . . . . . . . . 298
§1 Matematički uvod
( ) ( ∂z )2
-
2. Pretprostaviti da je z = z(x, y), odrediti funkcijsku vezu izmedu ∂x
i +
∂y ∂x y
z
2y.
(x, b) → (x, y), i (b) integraliti duž puta (a, b) → (a, y) → (x, y). Diskutovati dobi-
jene rezultate.]
̸dF
može se naći funkcija w(x, y), tako da je w(x,y)
totalni diferencijal. Odrediti takve
klase funkcija w(x, y) koje odgovaraju funkcijama
u(x, y) = x + y; v(x, y) = 1.
§2 Statistička termodinamika
p
T1 =const
S1=const S2=const
T2 =const
Slika 2.1:
1. Idealan gas je podvgnut kvazistatičkim procesima koji čine Carnot-ov ciklus. Kao
što je na grafiku predstavljeno, ciklus se sastoji iz dve adijabate i dve izoterme. Na
grafiku (vidi sliku 2.1) treba izdvojiti ”1-2”, ”2-3”, ”3-4” i ”4-1” procese koji redom
predstavljaju izotermski, adijabatski, a zatim izotermski i na kraju adijabatski
proces. Dokazati da važi
Q1 Q2
+ = 0, (2.1)
T1 T2
6 § 2. STATISTIČKA TERMODINAMIKA
(p1 , V1 , T1 ) →
|{z} (p2 , V2 , T2 ) → (p2 , V1 , T3 )
|{z} (2.2)
adijabatsko izobarno
sirenje sabijanje
a) Termomehanički sistemi
( )
3. 2 - adijabatske kompresibilnosti Ks = − 1
Pokazati da izmedu ∂V
i izotermske
( ) V ∂p
S
kompresibilnosti KT = − V ∂p
1 ∂V
postoji veza
T
Cp
KT = KS , (2.3)
Cv
jakobijana prelaza
( ) ( ) ( ∂S )
∂V ∂(V, p) ∂(V, T ) ∂(V, S) ∂(p, S) ∂V ∂T
= = = ( ∂S ) p . (2.4)
∂T p ∂(T, p) ∂(V, S) ∂(p, S) ∂(p, T ) ∂p S ∂T V
I.. PRVI I DRUGI PRINCIP TERMODINAMIKE 7
( ∂V ) ( ∂p )
(a) Cp − Cv = T ∂T p ∂T V
,
(( ∂p ) )2 ( ∂V )
(b) Cp − Cv = −T ∂T V ∂p
.
T
Rešenje:
∂(S, p)
Cp = T , (2.6)
∂(T, p)
∂(S, p) ∂(T, V )
Cp = T =,
∂(T, V ) ∂(T, p)
( ) ( ) ( ) ( ) ( )
∂S ∂p ∂S ∂p ∂V
= T( − ) , (2.7)
∂T V ∂V T ∂V T ∂T V ∂p T
( ) ( ) ( )
∂S ∂p ∂V
= Cv − T ,
∂V T ∂T V ∂p T
( ∂S ) ( ∂V ) ( )
∂S - da važi i
gde možemo iskoristiti da je ∂V = , a takode
( ) T
( ∂V )
∂p
T
∂p
T
Maxwell-ova relacija ∂S∂p
= − ∂T p
, i time možemo potvrditi traženu
T
relaciju
( ) ( )
∂V ∂p
Cp − Cv = T (2.8)
∂T p ∂T V
8 § 2. STATISTIČKA TERMODINAMIKA
( ) ( ∂V )
(b) Polazeći od izraza (2.7) i upotrebom Maxwell-ove relacije ∂S
∂p
=− ∂T p
,
T
dobija se
(( ) )2 ( )
∂p ∂V
Cp − Cv = −T (2.9)
∂T V ∂p T
αp = pαV KT ,
( )
gde je αp = V1 ∂V – koeficijent termalnog širenja pri konstantnom pritisku,
( )
∂T p
( )
αV = p ∂T V koeficijent termalne promene pritiska, dok je KT = − V ∂V
1 ∂p 1
∂p
T
iyotermske kompresibilnosti.
Rešenje:
Ova relacija se može dokazati, ako se zapiše koeficijent termalnog širenja u obliku
jakobijana ( )
1 ∂V 1 ∂(V, p)
αp = = , (2.10)
V ∂T p V ∂(T, p)
što se zatim može izraziti
1 ∂(V, p) ∂(V, T )
αp = , (2.11)
V ∂(V, T ) ∂(T, p)
(a) Izračunati Cp .
Rešenje:
(a) Cp = CV + kB .
(c) U (S, V ) = Uo + Cv T .
1 CV 1
(d) kT = p
i kS = Cp p
.
Rešenje:
Promenu unutrašnje energije sa temperaturom i pritiskom možemo napisati
( ) ( )
∂U ∂U
dU = dp + dT, (2.14)
∂p T ∂T p
što možemo zameniti u izraz za prvi princip termodinamike
δQ = dU + δA =
( ) ( ) (( ) ( ) )
∂U ∂U ∂V ∂V
= dp + dT + p dp + dT , (2.15)
∂p T ∂T p ∂p T ∂T p
gde smo i promenu zapremine napisali u zavisnosti od preomene p i T . Pošto je
po definiciji toplotni kapacitet pri konstantnom pritisku
( )
∂Q
Cp = , (2.16)
∂T p
onda nalazimo ( ) ( )
∂U ∂V
Cp = +p . (2.17)
∂T p ∂T p
10 § 2. STATISTIČKA TERMODINAMIKA
Rešenje:
5 R2 T v 3
Cp = R + . (2.19)
2 RT v 3 − 2(v − b)2 a
du = T ds − pdv. (2.21)
I.. PRVI I DRUGI PRINCIP TERMODINAMIKE 11
Cv ln T = −R ln(v − b) + C, (2.33)
9. Termička jednačina stanja i toplotni kapacitet gasa bozona koji pripadaju Bose
kondenzatu se mogu opisati izrazima
Rešenje:
( ∂p ) 3 1 3
a. dU (T, V ) = Cv dT + (T ∂T V
− p)dV = [(gT 2 + eT 2 )V + f T 2 ]dT + [ 32 aT 2 +
2bT 3 − c
V2
]dV .
15
b. g = 4
a i e = 6b.
5 5
c. U (T, V ) = U (To , Vo ) − ( 32 aTo2 + 2bTo3 )Vo − C
Vo
+ ( 23 aT 2 + 2bT 3 )V + VC .
10. Kalorička jednačina proizvoljnog idealnog kvantnog gasa može biti data izrazom
gU = pV , gde je g konstanta koja zavisi od osobina tog gasa. Pokazati da važi
jednačina U = V −g f (T V g ), gde je f (x) neka analitička funkcija argumenta x.
- odrediti jednačinu adijabastkog procesa ovog idealnog kvantnog gasa.
Takode,
-
11. Jedan od mogućih izbora medusobno nezavisnih termodinamičkih parametara
koji karakterišu stanje gasova i homogenih materijala su pritisak p, empirijska tem-
peratura t i zapremina v. U tom slučaju, razmenjena toplota nekog materijala sa
okolinom se može prikazati pomoću relacija
dQ = Cv dt + lv dv = Cp dt + lp dp = mv dv + mp dp, (2.36)
a. mv = Clpv−C
Cp
v
, mp = − Clpp−C
Cv
v
, mlvv + mp
lp
= 1;
( ) Cp −Cv ( ∂v )
b. ∂p
∂t v
= − lp
, ∂t p
= Cpl−C
v
v
.
14 § 2. STATISTIČKA TERMODINAMIKA
12. Pokazati da mala količina toplote δQ nije totalni izvod, tj da nije funkcija stanja
sistema.
dU = T dS − pdV, (2.37)
= 0. (2.42)
-
Time je potvrdeno da unutrašnja energija ne zavisi od zapremine pri konstantnoj
temperaturi.
14. Temperatura idealnog gasa diatomskih molekula molarne mase M , koji se nalazi
sudu velike visine, se menja sa visinom u skladu s izrazom T = To (1 − az), gde
je a poznata pozitivna konstanta, a To temperatura gasa na visini z = 0. Odrediti
kako se menja pritisak sa visinom, ako je p(z = 0) = po .
I.. PRVI I DRUGI PRINCIP TERMODINAMIKE 15
p+dp
dm
d mg
Slika 2.2:
Rešenje:
Usled mehaničke ravnoteže slojeva vazduha (vidi sliku (2.2)), može se zaključiti
da je
(p + dp)S + dmg = pS, (2.43)
Cp
15. Odrediti kolika je γ = Cv
, koristeći eksperimentalna merenja brzine zvuka u
nekom gasu na temperaturi T . Iskoristiti da je brzina zvuka se može shvatiti
kao √ -
brzina prenosa poremećaja kroz sredinu, i da se odreduje pomoću izraza
( )
∂p
c= ∂ρ
.
S
u = σT 4 . (2.46)
Pored toga uzeti u obzir da je pritisak za taj sistem povezan sa gustinom un-
utrašnje energije p = 13 u.
Rešenje:
S jedne strane znamo da je
( ) ( )
1 p 1 ∂U 1 ∂U p
dS = dU + dV = dT + ( + )dV, (2.47)
T T T ∂T V T ∂V T T
V T 3 = const. (2.54)
19. Gustina unutrašnje energije zračenja za apsolutno crno telo je dato jednačinom
U
u = V
= σT 4 . Uzimajući u obzir da za apsolutno crno telo važi µ = 0, odred-
iti termičku jednačinu stanja. Odrediti razmenjenu toplotu sistema u slučaju da
To
temperatura poraste sa 2
na To .
Rešenje:
Koristeći da je ( )
∂F
U (T, V ) = F − T (2.55)
∂T V,,
σV T 4
F =− + T · f (V ), (2.57)
3
f (V ) = V f ′ (V ), (2.59)
To
uzimajući u obzir da se temperatura menja od 2
do To , dobijamo da je
15 αTo
Q = σTo4 ( Vo − ). (2.61)
16 30
u pv
s = + =
T T
= aT 3 + R, (2.63)
u = T s − pv, (2.64)
daje izraz
sdT = vdp. (2.67)
T a (T 3 −To3 )
p = po e 3R . (2.69)
To
b. Izvesti izraz za rad koji gas vrši u adijabatskom procesu nad okolinom iz
termodinamičkog stanja (p1 , V1 , T1 ) u (p2 , V2 , T2 ). Petpostaviti da se molarna
specifična toplota pri konstantnoj zapremini ne menja.
Rešenje:
22. 3 Izračunati (a) brzinu promene temperature gasa sa pritiskom pri konstant-
( ) ( )
noj entropiji ∂T
∂p
, (b) odnosno pri konstantnoj entalpiji 1 ∂T
∂p
, u zavisnosti od
S E
koeficijenta toplotnog širenja αp i toplotnog kapaciteta pri kontantnom pritisku
Cp .
1
Odnosno tzv. Joule-Thomon-ov koeficijent.
I.. PRVI I DRUGI PRINCIP TERMODINAMIKE 21
Rešenje:
( )
∂S ∂(S, p)
Cp = T =T , (2.74)
∂T p ∂(T, p)
Metodom jakobijana, koeficijent toplotnog pritiska se tranformiše u
( )
∂T ∂(T, S) ∂(T, S) ∂(T, p)
= = =
∂p S ∂(p, S) ∂(T, p) ∂(p, S)
∂(T, S) ∂(T, p)
= =
∂(T, p) ∂(p, S)
( ) [ ( )]
∂S ∂T
= − =
∂p T ∂S p
[ ( ) ][ ]
∂V T
= − − =
∂T p Cp
αp
= VT , (2.75)
Cp
- ovih funkcija odgovora sistema.
čime smo odredili traženu vezu izmedu
sa temperturom i pritiskom
( ) ( )
∂S ∂S
dE = (V + T )dp + T dT =
∂p T ∂T p
( )
∂V
= (V − T )dp + Cp dT. (2.78)
∂T p
24. Kod izotermnog zračenja gustina unutrašnje energije u je monotona rastuća funkcija
temperature, a pritisak zavisi od gusine unutrašnje energije p = u3 . Kakav oblik
- na osnovu ovih podataka fenomenološka
funkcionalne zavisnosti u(T ) predvida
termodinamika? Da li je rezultat saglasan sa zakonima zračenja? Primeniti isto
- monotono rastuća
rezonovanje na idealan gas, kod koga je p = 2u i u je takode 3
RešenjeN S :
3
U = N kB T (2.81)
2
24 § 2. STATISTIČKA TERMODINAMIKA
pV = N kB T (2.82)
TS 5
µ=− + kB T. (2.83)
N 2
Helmholtz-ova slobodna energija je:
F (T, V, N ) = U − T S =
( [( ) 32 ( )− 52 ])
3 N 3N kB T V N
= N kB T − T S0 + N kB ln · · .
2 N0 2U0 V0 N0
G(T, p, N ) = U − T S + pV =
( [( )− 23 ( )− 52 ])
5 N 3N kB T N kB T N
= N kB T − T S0 + +N kB ln .
2 N0 2U0 pV0 N0
(( ) 3 ( )− 25 )
( ∂S ) 2
S0
Pošto je ∂N U,V = N0 + kB ln U
U0
· V0 · NN0
V
− 52 kB = − Tµ , dobijamo
korišćenjem (2.81):
( ) 32 ( )
3kB T V 25 µ S0 5
N= N exp + − (2.84)
2U0 V0 0 T kB N0 kB 2
Koristeći redom (2.83) i izraz za U (T, N ), a potom (2.84) dobijamo:
( )
TS 5
Ω(T, V, µ) =U − T S − − + kB T N =
N 2
| {z }
µ
= − kB T N (2.85)
( )3
3kB T 2 V 52
= − kB T N =
2U0 V0 0
( )
µ S0 5
· exp + − .
T kB N0 kB 2
I.. PRVI I DRUGI PRINCIP TERMODINAMIKE 25
X(T, p, µ) =Ω + pV =
= − N kB T + N k B T =
=0.
-
26. Virialni razvoj za slučaj razredenog gasa se može zadržati na drugom članu u
razvoju
N kB T N
p= [1 + B2 (T )], (2.86)
V V
gde je B2 (T ) drugi virialni koeficijent. Toplotni kapacitet će imati korekciju na
vrednost koju ima za slučaj idealnog gasa
3 N 2 kB
CV,N = N kB − F (T ). (2.87)
2 V
(a) Odrediti kakav oblik funkcija F (T ) ima da bi ove dve date relacije bile ter-
modinamički odgovarajuće.
I.. PRVI I DRUGI PRINCIP TERMODINAMIKE 27
kT = T vf (P ) (2.90)
N S
dp = dµ + dT.
V V
( ∂p ) ( )
S ∂p
Pošto je pritisak potpuni diferencijal, dobijamo tražene izraze ∂T µ
= V
i ∂µ
=
T
N
V
. Zamenom (2.84) u (2.82) dobijamo:
( ) 23 5 ( )
3kB T N02 µ S0 5
p(µ, T ) = exp + − · kB T
2U0 V0 T kB N0 kB 2
Slika 2.3:
-
(b) Odredujući dU (T, V ), pokazati da je unutrašnja energija U jedino funkcija
temperature T i broja čestica N .
Cp
(c) Pokazati da je γ = Cv
= 1 + N kCBv , odnosno da je nezavisna od T i V .
32. U ovom zadatku razmatramo infitezimalni Carnot-ov ciklus fotonskog gasa koja
je prikazana na Slici 2.3.
(b) Izraziti apsorbovanu toplotu Q, pri širenju gasa duž izoterme u zavisnosti p,
dV i odgovarajućih izvoda unutrašnje energije U (T, V ).
33. Za klasični približno idealni gas termička jednačina stanja ima sledeći vid viri-
jalnog razvoja ( )
RT 1
p= 1 − β1 (T ) + . . . (2.94)
v 2v
gde je β1 (T ) poznata funkcija temperature T . Dok je u slučaju jonizaovanog gasa
termička jednačina stanja:
( 1 ( ) 12 )
RT e3 π 2 1
p= 1− + ... (2.95)
v 2 vT
Odrediti toplotni kapacitet pri konstantnoj zapremini Cv (T, v), ako je poznato da
je u limesu kada v → ∞ toplotni kapacitet iznosi Cv (T ).
RešenjeN S :
Pošto je u tekstu zadatka jedino naveden izraz za pritisak, potrebno je povezati
promenu entropije sa promenom pritiska. S tim ciljem, koristeći da je Helmholtz-
ova slobodna energija funkcija stanja dF = −S dT − p dv, dobijamo Maxwell-
( ) ( ∂p )
ovu relaciju: ∂S
∂v T
= ∂T v
. Parcijalnim deferenciranjem leve i desne strane po
temperaturi nalazimo izraz
( ( ) ) ( )
∂ ∂S ∂ 2p
= . (2.96)
∂T ∂v T v
∂T 2 v
R dβ1 RT 2 d2 β1
Cv (T, v) = Cv (T ) + T +
v dT 2v dT 2
gde je nm broj molova gasa. Slično, uvršćavanjem izraza (2.95) u jednačinu (2.98),
nakon diferenciranja i integracije dobijamo:
1
Re3 π 2
Cv (T, v) = Cv (T ) − √
4 Tv
b) Elastične trake
(a) Količinu toplote Q koju traka dobija od okoline pri sporom istezanju od lo
do l na temperaturi To .
ATo 2 2 1 1
Q=− (l − lo ) + AT lo (αAT lo2 (To − 2T ) − 1)( − ) + Cl (T − To ), (2.101)
2lo lo l
ATo 2 1 1
Q=− (l − lo2 ) − ATo lo (αATo2 lo2 − 1)( − )(1 + αTo ). (2.102)
2lo lo l
(b) Pošto je istezanje trake brzo proces možemo smatrati adijabatskim. Tada
imamo ( ) ( )
∂S ∂S
dl + dT = 0. (2.103)
∂l T ∂T l
ATo 2 Al2 To 1 1
∆T = (l − lo2 ) + o ( − ). (2.104)
2lo Cl Cl lo l
35. Razmotriti elastičnu oprugu čija sila data jednačinom F (x, T ) = k(T )x, gde je x
istezanje opruge. Odrediti unutrašnju energiju U i entropiju S.
Rešenje:
Na osnovu date zavisnosti slobodne energije, nalazimo da je unutrašnja energija
data oblikom
x2
U (x, T ) = U (0, T ) + (k(T ) − T k ′ (T )), (2.105)
2
dok je entropija
S(x, T ) = S(0, T ) − xk(T ). (2.106)
I.. PRVI I DRUGI PRINCIP TERMODINAMIKE 33
36. Kod elastičnih gumenih traka, sila zatezanja F zavisi od dužine l i temperature T
kao
f = αx − AT + BT x, (2.107)
da(x)
a. Pokazati da važi dx
= 0, dnosno da a(x) ne zavisiod x.
Rešenje:
dU = T dS + f dx, (2.108)
dU = Cx dT + αxdx. (2.111)
( x) ( ∂αx )
Da bi dU bio totalni diferencijal, potrebno je da važi da je ∂C ∂x T
= ∂T x
,
( ∂Cx )
odakle se dobija da je ∂x T = 0, što implicira da je da(x)
dx
= 0.
34 § 2. STATISTIČKA TERMODINAMIKA
x2
S = So + aT + Ax − B . (2.114)
2
c. Iskoristićemo identitete
( )
∂S ∂(S, x) ∂(f, T )
= · (2.115)
∂T x ∂(f, T ) ∂(T, x)
( ) ( ) ( )
∂S ∂S ∂x
= − , (2.116)
∂T f ∂x T ∂T f
- toplotnih kapaciteta Cx pri konstantnoj dužini, odnosno Cf
za vezu izmedu
pri konstnoj sili. Posle nekoliko jednostavnih računskih koraka, dobija se
(A − Bx)2
Cf = T (a + ). (2.117)
α + BT
Rešenje:
Entropija elastiňe trake je
lo x2 lo
S(T, x) = C ln T − ( 2 + ). (2.119)
Tl 2lo x
Slika 2.4:
( ) [( ) ]
1 ∂U 1 ∂U
dS = dT + − axT dx
T ∂T x T ∂x T
X
Nakon integracije, korišćenjem izraza T = ax
dobijamo
Cx X0
Q23 = , (2.121)
2al0
izvršeni rad
η=
primljena energija
Q12 + Q23 + Q31
=
Q23 + Q31
al02
=
Cx + 3al02
i relacija
( ) ( ( x )3 )
∂X ax 0
= 1− , (2.124)
∂T x xo x
gde je X – sila koja deluje na razmatranu traku, dok je sa x obeležena dužina
trake. Pored toga, pretpostaviti da su a = 1 N
K
i x0 = 0.5m poznate i pozitivne
- i da se masa trake ne menja.
konstante, a takode
1
( ∂x )
(a) Odrediti termalni koeficijent širenja trake αX = x ∂T X
. Prodiskutovati
promenu αX .
40. Razmotriti štap koji je dužine lo u neistegnutom stanju, i koji se istegne za dužinu
X
x= , (2.125)
aT 2
- uzeti
gde je X sila koja deluje na štap, T temperatura štapa i a konstanta. Takode,
da se toplotni kapacitet pri konstantnoj dužini ponaša u skladu sa jednačinom
Cx (x = 0, T ) = bT, (2.126)
gde je b konstanta.
X
x= , (2.127)
aT 2
- uzeti
gde je X sila koja deluje na štap, T temperatura štapa i a konstanta. Takode,
da se toplotni kapacitet pri konstantnoj dužini ponaša u skladu sa jednačinom
Cx (x = 0, T ) = bT, (2.128)
gde je b konstanta.
41. Pokazati da je slobodna energija po jedinici površine neke tečnosti jednaka površinskom
naponu γ te tečnosti. Pokazati, takodje, da je unutrašnja energija po jedinici
površine jednaka γ − T dT
dγ
.
Rešenje:
Krenimo od izraza za unutrašnju energiju površine
s F T
U + - G
S E g
Slika 2.5:
-
koja se može uz Legendre-ovu relaciju iskoristiti Fσ = Uσ − ST u odredivanju
relacije
dFσ = −SdT + γdσ. (2.130)
40 § 2. STATISTIČKA TERMODINAMIKA
Fσ
γ= , (2.135)
σ
42. U slučaju da je izvršeni rad nad slobodnom površinom tečnosti dA = γdσ, gde je
∂γ
γ(T, σ) površinski napon i ∂T
< 0.
-
43. Jedna od potvrdenih i često korišćenih fenomenoloških izraza za unutrašnju en-
ergiju čvrstog tela, je data izrazom
AT 3
44. Termička jednačina stanja neke nove materije je p = 3V
, dok je kalorička jednačina
stanja U = BT n ln( VVo ) + f (T ), gde su B, n i Vo konstante a f (T ) zavisi od temera-
ture. Odrediti B i n.
Rešenje:
Polazeći od izraza za entropiju sistema
( ) ( )
∂S ∂S
dS = dV + dT, (2.138)
∂V T ∂T V
može se iskoristiti uslov da je sa leve strane izraza totalni diferencijal, odnosno
da je [( ) ] [( ) ]
∂ ∂S ∂ ∂S
= . (2.139)
∂T ∂V T ∂V ∂T V
Partcijalni izvodi entropije po temperaturi i zapremini se dobijaju iz izraza za
neki od termodinamičkih potencijala, dobijamo da je
( ) ( )
∂S 1 ∂U AT 2
= +
∂V T T ∂V T V
( ) ( )
∂S 1 ∂U
= ,
∂T V T ∂T V
42 § 2. STATISTIČKA TERMODINAMIKA
-
45. U odredenim sistemima kvazičestica, unutrašnja energija je povezana sa entropi-
jom izrazom
( )d
N dS
U (V, S) = AN e N kB , (2.140)
V
- odrediti
gde su A i d konstante. Napisati temičku jednačinu stanja. Takode,
toplotne kapacitete Cp i Cv sistema.
51. Tanki sloj sapunice je pridržavan od strane četiri srane pravougaonika ABCD
(vidi sliku), čije su ivice tanki štapovi. Štap AD, dužine a, se može pokretati
paralelno štapu BC tako da se sloj sapunice time može rastezati. Intenzitet sile
|f⃗| = σl kojom se deluje u pozitivnom smeru x–ose, omogućuje da se štap ostane
nepokretan. Sa σ je obeležen površinski napon.
(b) Ako je poznat toplotni kapacitet pri konstantnoj dužini Cx (T ) napisati kako
( ) ( ∂U )
zavise ∂U
∂T x
i ∂x T preko parcijalnih izvoda površinskog napona (po istim
paramettrima).
σ = σo (1 + a(T − To )) (2.145)
d) Termodinamika magnetika
52. Odrediti koliki je rad izvršen nad paramagnetnim dijelektrikom pri njegovoj po-
larizaciji i namagnetisavanju.
Rešenje:
Za diksontinualno naelektrisanje dijamagnetne elektrčne i magnetne osobine su
-
odredene dielektričnim i magnetnim momentom
∑ ∑
P
⃗ = qi⃗ri , ⃗ =1
M qi (⃗ri × ⃗vi ), (2.146)
i
2 i
dok se za kontinualno naelektrisanje, sume prelaze u integralima
∫ ∫
1
P = ⃗rρ(⃗r)d ⃗r,
⃗ 3
M=
⃗ (⃗r × ⃗v )ρ(⃗r)d3⃗r. (2.147)
2
Električna polarizacija za kontinualno naelektrisanje se može uvesti izrazom
dP⃗
⃗ = dM .
⃗
P⃗ = 3 , M (2.148)
d ⃗r d3⃗r
Ukupni rad koji izvrši električno i magnetno polje nad nekim magnetikom se sas-
toji od ”elementarnog” rada koje izvrši magnetno polje na svaku pojedinačnom
česticom putem Lorentz-ove sile
⃗ + ⃗v × B),
F⃗L = q(E ⃗ (2.149)
I.. PRVI I DRUGI PRINCIP TERMODINAMIKE 45
što nam omogućuje da uvedemo polje Lorentz-ove sile u svakoj tački magnetika
dF⃗L
f⃗L = 3 = ρ(⃗r)(E
⃗ + ⃗v × B).
⃗ (2.150)
d ⃗r
Time je rad koji izvrši fluks naelektrisanja neke male, odnosno delića zapremine
magnetika d3⃗r, u nekoj proizvoljnoj tački ⃗r je
δw = f⃗L d⃗r,
Pošto je ⃗v · (⃗v × B)
⃗ = 0, dobijamo da je
∫
δW = ⃗ v dtd3⃗r =
ρ(⃗r)E⃗
∫
= ⃗
⃗j(⃗r)Edtd 3
⃗r, (2.151)
M F T
U + - G
S E H
Slika 2.6:
dU = T dS + HdM, (2.156)
odkle se dolazi do
( ) ( ) ( )
∂s ∂u ∂m
T = −H . (2.160)
∂T ∂T ∂T
| {z }
H H H
CH =
( ∂u )
Parcijalni izvod ∂T H
se može dobiti prelaskom u zapis sa jakobijanima prelaska
( )
∂u ∂(u, H)
= =
∂T H ∂(T, H)
∂(u, H) ∂(T, m)
= =
∂(T, m) ∂(T, H)
[( ) ( ) ( ) ( ) ] ( )
∂u ∂H ∂u ∂h ∂m
= − · =
∂T m ∂T m ∂m T ∂T m ∂H T
( ) ( ) ( ) ( )
∂u ∂u ∂H ∂m
= − =
∂T m ∂m T ∂T m ∂H T
| {z }
∂(H,m) ∂(M,T ) ∂(H,m)
=− ∂(T,m) =− ∂(H,T )
)
( ( ) ( )
∂(H,T )
∂u ∂u ∂m
= + . (2.161)
∂T m ∂m T ∂T H
| {z }
=Cm
tako nalazimo ( ) ( )
∂u ∂H
=T + H. (2.163)
∂m T ∂T m
48 § 2. STATISTIČKA TERMODINAMIKA
dU = T dS + HdM, (2.167)
dobijamo
( ) ( ) ( ) ( )
∂U ∂S ∂U ∂S
=T , =T + H. (2.168)
∂T M ∂T M ∂M T ∂M T
i nakon integraljenja
∫ ∫ H/T
dU (T )
S = − xf ′ (x)dx + const =
T 0
∫ H/T
H H
= g(T ) − f ( ) + f (x)dx + const. (2.171)
T T 0
Uzimajući da je
( )
d F
U (T, M = 0) = −T 2
, (2.175)
dT T
dobijamo da je
( )
M 2T 2 d 1
U (T, M ) = U (T, M = 0) + . (2.176)
2 dT T χT (T )
Rešenje:
δQ = T dS = dU − HdM, (2.179)
χ
QT =const = − H 2 . (2.181)
2
χH 2
U (T2 ) − U (T1 ) = − , (2.182)
2
C
57. Razmatrati paramagnetni materijal koji ima izotermsku susceptibilnost χT = T
,
K
i ima toplotni kapacitet C(T, M = 0) = T2
.
(a) Odrediti toplotni kapacitet pri konstantnom magnetnom polju CH , i pri kon-
stantnoj nenultoj magnetizaciji CM .
C
59. Razmatrati paramagnetni materijal koji ima izotermsku susceptibilnost χT = T
, i
K
ima toplotni kapacitet C(T, M = 0) = (T +θ) 2 . Odrediti adijabatsku susceptibilnost
( ∂M )
pri χS = ∂H S u zavisnosti od magnetnog polja H i temperature T .
C
60. Magnetna susceptibilnost paramgnetika je data Kirijevim zakonom χ = T
, a
topolotni kapacitet po jedinici zapremine pri konstantnom polju CH = (a+bH 2 ) TV2
(gde su a i b konstante). Odrediti kolika je temperatura T ′ posle adijabatskog raz-
magnetisavanja sa temperature T , pri samnjivanju polja sa H na 0.
A
61. Magnetna susceptibilnost antiferomagnetika se može predstaviti kao χ = Tα
, a
gde je A > 0, i gde je 0 < α < 1. Unutrašnja energija van spoljašnjeg polja je
U (T, M = 0) = BT β , gde su B, β > 0. Ako je sistem imao temperaturu T u
′
polju jačine H, posle adijabatskog isključenja ohladiće se do temperature T < T .
I.. PRVI I DRUGI PRINCIP TERMODINAMIKE 53
′
Odrediti temperaturu T . Uzeti u obzir da je A i α ne zavise od H (odnosno, da
su slaba polja), što znači da su efekti magnetostrikcije mali.
( )
∂µ
= −s, (2.186)
∂T H,p
( )
∂µ µo H
= −m = − . (2.187)
∂H p,T 3kB T
Koristeći prethodne relacije možemo dobiti hemijski potencijal u funkciji svojih
prirodnih parametara
µo 2
µ(p, T, H) = kB T ln p − H − cp T ln T + cp T. (2.188)
kB T
Gas koji se nalazi u magnetnom polju i izvan magnetnog polja je u termodi-
namičkoj ravnoteži, tako da je T1 = T2 , i µ1 = µ2 . Koristeći izraz za hemijski
potencijal, u ovom slučaju termodinamičke ravnoteže možemo napisati relaciju
µ(p1 , T1 , H1 ) = µ(p2 , T2 = T1 , H2 = 0), gde su veličine indeksirane sa 1 su un-
utar solenoida, dok veličine indeksirane sa 2 su izvan solenoida. Ako iskoristimo
izraz za hemijski potencijal dobijamo relaciju
µo H
p1 = p2 e 3(kB T )3 , (2.189)
( ) ( ) ( ) ( )
∂S ∂S ∂S ∂µ
= − , (2.193)
∂T V,µ ∂T V,N ∂µ T,V ∂T V,N
gde je x = Tµ .
Rešenje:
Rešenje:
dU = T dS − pdV, (2.203)
odakle se dobija da je
( ) ( )
∂µ ∂p
= −V . (2.207)
∂N T,V ∂N T,V
a
dH = Cv dT + dV + d(pV ). (2.212)
V2
f) Termodinamičke nejednakosti
( ∂S )
(a) Cv = T ≥ 0,
∂T p
) (
(b) KS = − V ∂p
1 ∂V
≥ 0,
S
( )
(c) KT = − V1 ∂V
∂p
≥ 0.
T
Rešenje:
(b) Polazeći od toga da kvadratna forma (2.213) mora biti pozitivna, dobijamo
( ) ( )
∂ 2U ∂ ∂E ∂p
= =− ≥ 0, (2.215)
∂V 2 ∂V ∂V S ∂V S
odakle zaključujemo da je KS ≥ 0.
- funkcija odgovora sistema (2.3), i na
(c) Pošto znamo da važi odnos izmedu
osnovu toga da je KS ≥ 0, i na osnovu Cp > CV ≥ 0, dobijamo da je i
KT ≥ KS ≥ 0.
odakle se koristeći
∑ ∂2U ∑ ∂Yi
δyk = δyk = δYi ,
k
∂y i ∂yk
k
∂y k
I.. PRVI I DRUGI PRINCIP TERMODINAMIKE 61
dolazi do
1∑
δ2 U = δYi δyi . (2.218)
2 i
δ2 U > 0, (2.219)
odnosno ako se izraz (2.221) razmatra kao kvadratna forma, onda determnanta
koeficijenata u toj kvadratnoj formi mora biti pozitivna
∂2U ∂2U ∂2U
···
∂yk1 ∂yk1 ∂yk1 ∂yk2 ∂yk1 ∂ykn
∂2U ∂2U ∂2U
···
δ2 U =
∂yk2 ∂yk1 ∂yk2 ∂yk2 ∂yk2 ∂ykn ,
(2.220)
··· ··· ··· ···
∂2U ∂2U ∂2U
∂ykn ∂yk1 ∂ykn ∂yk2
··· ∂ykn ∂ykn
( ) ( )
∂yi ∂yk
za veličine za koje važi ∂Yk
= ∂Yi
.
Yi Yk
Rešenje:
62 § 2. STATISTIČKA TERMODINAMIKA
( ∂X )
(a) ∂x S,p
> 0,
( ∂X )
(b) ∂x T,p
> 0.
II.. TREĆI PRINCIP TERMODINAMIKE 63
Slika 2.7: Entropija paramagnetnih soli kao funkcija temperature za nekoliko ra-
- B = 0
zličitih jačina magnetnih polja (magnetnih indukcija u intervalu izmedu
-
i maksimalne vrednosti Bb ). Magnetno hladenje paramagnetnih soli od tempera-
ture T = Ti do temperature T = Tf se sprovelo u dva koraka: prvo, izotermskom
magnetizacijom koja se vrši povećanjem magnetnog polja od B = 0 do B = Bb
na konstantnoj temperaturi Ti ; a zatim adijabatskom demagnetizacijom pri kojoj se
smanjuje temperatura od T = Ti do T = Tf .
75. 5 Pokazati da toplotni kapacitet teži nuli, u limesu kada T → 0. Koristiti III
princip termodinamike.
Rešenje:
Iz trećeg principa termodinamike, odnosno pri nekom konstantnom parametru
64 § 2. STATISTIČKA TERMODINAMIKA
lim CX = 0. (2.227)
T →0
( )
∂V
= 0. (2.232)
∂T p
T →0
( ∂p )
S ciljem da se dokaže ∂T V
−−−→ 0, treba krenuti od izraza za toplotni kapacitet
T →0
pri konstantnoj zapremini CV , a zatim sličnim putem dokazati traženu graničnu
vrednost.
dχ(T ) H 2
S(T, H) = S(T, 0) + . (2.237)
dT 2
I. Clausius–Clapeironova jednačina
2. 5
(a) Kako zavisi pritisak na liniji koegzistencije gasne i čvrste faze, uzimajući da
se jednačinom idealnog gasa može predstaviti stanje gasne faze razmatra-
nog materijala.
(b) Ovaj problem se može primeniti na H2 u blizini trojne tačke (p0 , T0 ) = (90M P a, 14K).
kJ
Latentna toplota prelaza čvrsto-gas stanje u tom slučaju q = 1.01 mol .
Rešenje:
dp sg − ss q
= = , (3.5)
dT vg − vs T (vg − vs )
dp q p
= · 2, (3.6)
dT R T
(b) Uzimajući rezultat koji je dobijen u prethodnog delu ovog zadatka možemo
napisati kakvog je oblika jednačina koegzistencije faza
Rešenje:
(a) Na trojnoj tački, pritisak i temperatura pare, tečnosti i čvrstog stanja su iste.
Zbog toga, temperatura trojne tačke T0 se dobija iz preseka krivih subli-
macije i isparavanja, odnosno rešavajući jednačinu
3726K 3005K
27.79 − = 24.10 , (3.14)
T0 T0
(b) Iz nagiba krive koegzistencije faza čvrste i gasne faze i korišćenjem relacije
(3.6) dobijene pod pretpostavkom da se gas može smatrati idealnim, dobija
se latentna tolota sublimacije qsg ≈ 31 mol
kJ
. I na isti način, iz krive koegzis-
tencije tečne i gasne faze se može odrediti latentna toplota isparavanja qlg ≈
kJ
25 mol .
s l
Slika 3.1:
dq ∂q ∂q dp
= ( )p + ( )T ( )p(T ) .
dT ∂T ∂p dT
dq q
= sB − sA + CpB − CpA + (αp,B vB − αp,A vA ) ,
dT vB − vA
odakle možemo dobiti traženi izraz. Ukoliko je jedna od faza ima za istu količinu
materije ima mnogo manju zapreminu vA << vB , a takode - i entropiju sA << sB i
termalni koeficijenti širenja αp,A << αp,B može se doći do aproksimativnog izraza
za promenu latentne toplote
dq q
= CpB − CpA + − αp,B q. (3.19)
dT T
U slučaju da se faza B može razmatrati kao idealni gas, onda se termalni koefici-
jent širenja može odrediti αp,B = T1 , odakle dobijamo
dq
= CpB − CpA . (3.20)
dT
Toplotni kapacitet duž linije koegzistencije faza se može izračunati pomoću izraza
( )
∂s
C=T , (3.22)
∂T p(T )
možemo dobiti desnu stranu izraza (3.22) kada prethodni izraz podelimo sa dT ,
i zahtevamo da se ta jednačina može dobiti duž linije koegzistencije faza p(T ),
odnosno ( ) ( ) ( ) ( )
∂s ∂s ∂s dp
= + . (3.24)
∂T p(T ) ∂T p ∂p T dT p(T )
( )
∂s
Uz korišćenje Clausius-Clapeiron-ove jednačine, i Maxwell-ove relacije ∂p
=
( ∂v ) T
− ∂T p
dolazimo do izraza
q
C = Cpg − v g αp , (3.25)
vg − vl
1
( ∂v )
gde je αp = v ∂T p
termalni koeficijent širenja gasa. Za slučaj idealnog gasa,
u ovoj relaciji možemo zameniti konkretni oblik zavisnosti koeficijenta širenja
1
αp = T
i da u potrebimo izraz 3.20, čime nalazimo
dq q d q
C = Cpl + − = Cpl + T . (3.26)
dT T dT T
(a) Odrediti kako pritsak p duž linije faznih prelaza zavisi od T , zapremine
tečnosti vL i vg .
(c) Proceniti vrednost latentne toplote q za T << Tc , kada možemo uzeti vg >>
b i vl ∼ b.
(d) Na osnovu prethodnih delova ovog problema, naći jednačinu linije koegzist-
ncija faza.
U = T S − pV + µA NA + µb NB .
Na krivi koegzistenciji dve faze p(T ) ova jednakost se održava. Drugim rečima,
pri promeni temperature važi:
Sg V g dp NAg g Sl V l dp NAl l
− + − dµ = − + − dµ (3.31)
NBg NBg dT NBg A NBl NBl dT NBl A
Pošto je entropija aditivna veličina važi S g = sgA NAg + sgb NBg . Slično nalazimo
entropiju tečne faze. Pretpostavimo da zapreminu odgovarajućih faza možemo
da nadjemo kao V l = vAl NAl + vBl NBl , i slično za gasnu fazu. Zamenićemo dobijene
izraze u (3.31) i napisati izraz u obliku simetričan deo po A i B + nesimetričan
2
Sa indeksima 1 i 2 označavam redom veličine koje se odnose na neku prvu i drugu fazu. Ovo se
sve može uopštiti i na prozvoljan broj faza.
76 § 3. FAZNI PRELAZI
∂µgA
deo. Pošto se po uslovu zadatka broj čestica A u gasnoj fazi ne menja, važi ∂T
=
∂µlA
∂T
. Tako dobijamo:
dp [ g ]
(vA − vAl )NAg + (vBg − vBl )NBg = (sgA − slA )NAg
dT
+(sgB − slB )NBg
[ ]
∂µgA 1
+ ∂T NBl − dT vA (NAg NBl − NAl NBg )
dp l
(3.32)
11. Termomehanički sistem ima liniju faznog prelaza drugog reda (tzv. lambda lin-
iju) koja razdvaja dve faze A i B tog sistema. Molarni toplotni kapacitet Cp i
termički koeficijent širenja αp je drugačiji u te dve faze. Izračunati nagib te krive
dp
( dT )koegz te λ − − linije u zavisnosti od temperature T , molarne zapremine v,
p − Cp i ∆αp = αp − αp .
∆cp = cA B A B
Rešenje:
-
Nagib krive se odreduje relacijom
dp 1 CpB − CpA
= .
dT vT αpB − αpA
-
gde oznaka koegz ukazuje da se izvod odreduje duž krive koegzistencije, duž
koje je gas u ravnoteži sa tečnom fazom. Odrediti aproksimativnu i eksplicitnu
vrednost ovog koeficijenta, koristeći jednačinu gasa idealnog gasa.
15. Helmholtz-ova slobodna energija (po jednom molu) čvrste faze je data izrazom
B B
fs = T v3
, dok je slobodna energija tečne faze (po molu) izrazom fl = T v2
. U ovim
izrazima, su uvedne konstante A i B, dok je sa v zadata molarna zapremina.
a
(p + )(v − b) = RT. (3.34)
T v2
prelaza.
19. Za van der Waals-ov gas, skicirati izotermu u p̄ − V̄ ravni (gde su p̄ i V̄ reduko-
vani pritisak i zapremina) za redukovane temperature T̄ = 0.5, T̄ = 1.0, i T̄ = 1.5.
Uzeti da je za T̄ = 0.5, vrednost ravnotežnog pritiska p̄ = 0.1 u oblasti koegzis-
tencije tečnost–gas.
a
(p + 1 )(v − b) = RT. (3.35)
T v(v + b)
2
P (v − b) = kT e−a/vkT , (3.36)
PV
a. Odrediti odnos N kB T
na kritičnoj tački za gas koja opisuje tu jednačinu. Daj
numeričku vrednost tih triju figura.
N kb T N
P = [1 + B(T )] (3.37)
V V
kB T ∑ m2i+1
∞
Tc
h= ( − m). (3.41)
µo i=0 2i + 1 T
kb Tc 1 2 1 kb T 4 1 kB T 6
g(m, T ) = f (T, m = 0) + ( − )m − m − m + ··· . (3.44)
µo 2 4 µo 6 µo
( )
∂G
h(T, m) = = 2b(T )m + 4c(T )m3 , (3.46)
∂m T
1
χ(T, m) = . (3.47)
2b(T ) + 12c(T )m2
Razvijajući u red zavisnost koeficijenata b(T ) i c(T ) u okolini faznih prelaza dobija
se da je
′
• Uvodeći pogodnije parametre b̃ = 2b′ Tc i c̃ = b2c
′ T , za h = 0 i τ < 0, dobija se
c
√ τ
m = − c̃ , što ukazuje da je kritični eksponent β = 12 .
82 § 3. FAZNI PRELAZI
Slika 3.2:
-
24. Rešavanjem jednačine (3.38) se odreduje magnetizacija u aproksimaciji teoriji sred-
zJ
njeg polja. Kritična temperatura je povezana sa parametrom interakcije Tc = kB
,
gde je z koordinacioni broj rešetke, što je zapravo broj čvorova sa koliko je svaki
čvor rešetke spojen3 . Grafički rešiti jednačinu
b(t) 2 c(t) 4
GM F (h, t; m) = a(t) + m + m − hm, (3.52)
2 4
3
Kod regularnih rešetki koordinacioni broj z = 2d, gde je d–dimenzionalnost sistema koji se raz-
matra. Kritične temperature se mogu oddrediti i bez primene aproksimacije srednjeg polja.
Tako, za d = 1 postoji imamo Tc = 0, što je manje nego u aproksimaciji srednjeg polja
Tc (d = 1) < Tc,M F (d = 1) = 2J. U slučaju d = 2 sistema Lars Onsager je 1944 dobio da
postoji nenulta Tc ≈ 2.269J, Tc (d = 2) < Tc,M F (d = 2) = 4J.
II.. OSNOVE LANDAU-OVE TEORIJE FAZNIH PRELAZA 83
Slika 3.3:
-
(a) Na Slici (3.2) je prikazano kako se odreduje magnetizacija za T > Tc . Dobija
se jedna nenulta magnetizacija, dok je Landau-ova energija predstavljena na
insetu grafika. Jedan globalni minimum se može dobiti i za dovoljno veliko
h.
-
(b) Na Slici (3.3) je prikazano kako se odreduje magnetizacija za T < Tc . Dobi-
jaju se dve ili tri nenulte magnetizacije u zavinosti koliko je magnetno polje
h. Globalni minimum se dobija za magnetizaciju m > 0, dok ostali ekt-
stremumi su lokalni minimumi Landau-ove slobodne energije.
-
(c) Slika (3.4) pokazuje tipično ponašanje neuredene faze, odnosno za tempera-
ture T > Tc . Dobija se samo jedan globalni minimum Landau-ove slobodne
energije za m = 0.
-
(d) Slika (3.5) pokazuje ponašanje u uredenoj fazi, odnosno za temperature T <
84 § 3. FAZNI PRELAZI
Slika 3.4:
t 1
G(m, t) = m2 + m6 − hm, (3.53)
2 6
T −Tc
gde je t = Tc
.
26. Pod uslovima prethodnog zadatka, odrditi zavisnost razlike toplotnih kapaciteta
T −Tc
Cm − Ch od τ = Tc
.
Rešenje:
Koristeći identitet
( ) ( )
∂S ∂m
Cm − Ch = −T , (3.54)
∂m T ∂T h
II.. OSNOVE LANDAU-OVE TEORIJE FAZNIH PRELAZA 85
Slika 3.5:
b̃T τ
Cm − Ch = − . (3.55)
2Tc2 c̃
a b c
G(m, t) = m2 − m3 + m4 , (3.56)
2 3 4
V
28. Tečnost ima ravnotežnu gustinu koja odgovara zapremini vo (T ) = N
. Njena slo-
bodna energija je oblika
(c) Nematički fazni prelazi se dešavaju kada se nematički faza (S ̸= 0) ima nižu
slobodnu energiju nego izotropna faza (S = 0). Odrediti vrednosti rc temper-
aturskog parametra na kojoj se prelaz dešava. Odrediti vrednost Sc parame-
-
tra uredenja na tom prelazu.
II.. OSNOVE LANDAU-OVE TEORIJE FAZNIH PRELAZA 87
Rešenje:N S
0.4
0.0
-0.4
f
(r -r ) / r =
critical critical
-0.8
0.350
0.069
0.00
-1.2
-0.438
-1.11
-1.6
-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5
Slika 3.6: Landauova slobodna energija pri različitim vrednostima parametra r. Cr-
2w2
vena (treća puna odozgo) linja odgovara vrednosti rkriticno = 9u
pri kojoj dolazi do
w2
faznog prelaza. Iznad te linije, za r < 4u
(isprekidana linija odgovara jednakosti), uz
2
ravnotežno stanje S = 0, postoji metastabilno stanje za koje je S < 0. Za r ∈ ( 2w
9u
, 0)
(ljubičasta linija) metastabilno stanje je odredjeno sa S = 0 a stabilno sa S < 0.
Pošto su oba lokalna minimuma stabilna na male promene S, da li će sistem ostati
u metastabilnom stanju ili će preći u stabilno zavisi od veličine fluktruacija koje ter-
modinamička teorija bez statističke analize ne može da odredi. Za r < 0 S = 0 više
nije stabilan minimum (svetlo plava linija), te će i najmanja fluktruacija da sistem
dovede u stanje S ̸= 0. Tada postoji još jedno metastabilno stanje za S > 0
(a =)f = −T ln 2
zJ
(b =)r = (T − zJ) ≈ T − Tc
T
zJ
(c =)n = zT ( )4 ≈ 2Tc
T
b c ∂u
G = a + M 2 + M 4 −→ H = = bM + cM 3
2 4 ∂M
∂H 1
χ−1 = = b + 3cM 2 ⇒ χ =
∂M b + 3cM 2
Sada ćemo odrediti kritične eksponente:
√
b
H = 0 ⇒ b = −cM ⇒ M =2
∼ (Tc − T )−1 ⇒ β = 1
−c
t > 0 ⇒ χ ∼ t−1 ⇒ γ = 1
t = 0 ⇒ H ∼ M3 ⇒ δ = 3
⇒ LM F ∼ Tc (T − Tc )0 ⇒ α = 0
1 1 1
L = am2 + bm4 + cm6 − hm
2 4 6
Rešenje:
∂L
0 = γL = dm
∂m
am + bm3 + cm5 − h = 0
Za h = 0 imamo da je:
a + bm2s + cm4s = 0
1 √
m2s = (−b ± b2 − 4ac)
2c
Uzimamo samo pozitivno rešenje:
√
−b + b2 − 4ac
m2s =
2c
II.. OSNOVE LANDAU-OVE TEORIJE FAZNIH PRELAZA 91
(d) SLIKE
ams + cm5s − h = 0
Za h = 0 i a < 0 imamo:
a 1 a1 t + a2 p 1
ms = (− ) 4 = (− )4
c c
Odavde dobijamo da je β = 41 .
Za p = t = 0, odnosno a = 0 imamo:
h 1
ms = ( ) 5
c
Odakle zaključujemo da je δ = 5.
Susceptibilnos
∂ms
(a + 5cm4s ) −1=0
∂h
92 § 3. FAZNI PRELAZI
Odavde sledi:
∂m 1
=
∂h a + 5cm4s
∂m 1 1
Za a > 0 i ms = 0 imamo da je ∂h
= a
, a za a < 0 i ms = (− ac ) 4 imamo
∂m
∂h
= − 4a
1
. Pošto je a = a1 t + a2 p, dobijamo γ = γ ′ = 1.
Tada je Landau-ova slobodna energija:
1 1 1 c
L = am2 + cm6 = m2 (a + m4 )
2 6 2 3
1
Za a > 0 i ms = 0 dobijamo da je L = 0, dok za a < 0 i ms = (− ac ) 4 imamo:
1 a 1 1 a
L = (− ) 2 (a + c(− )) =
2 c 3 c
1 a 1 1 a 3 1 a1 t + a2 p 3
= a(− ) 2 = − c(− ) 2 = − c(− )2
3 c 3 c 3 c
Toplotni kapacitet se nalazi kao drugi izvod po L:
a1 t + a2 p − 1
C ∼ (− ) 2
c
Za a > 0 imamo da je C = 0 pa odatle sledi da je α = 0, a za a < 0 je
C ∼ (−t)− 2 pa je onda α′ = 12 .
1
1
gde je Φ generalna homogena funkcija Φ(λτ, λM 1/β ) = λγ Φ(τ, M β ). Pošto je rel-
T −Tc
ativno rastojanje temperature od faznog prelaza je obeleženo sa τ = Tc
, onda
pokazati da se skalirajuće relacije mogu zapisati sa
(a) M ≈ τ β , za H = 0 i τ < 0,
Rešenje:
koja implicira da je
λM 1/β = φ(λτ ). (3.66)
III.. HIPOTEZA SKALIRANJA ZA TERMODINAMIČKE FUNKCIJE 95
H = M Φ(t, M α ),
gde je Φ funkcija
Φ(λt, λM α ) = λb Φ(t, M ).
1 1
Za T < Tc i za H = 0 imamo da je 0 = Φ(t, M β ), odnosno M β = φ(t).
1 1
M = t a iβ = .
a
-
Odredujući prvi izvod po magnetizaciji
( ) 1
∂H 1 ∂Φ(t, M β ) 1 β1 −1
= Φ(t, M ) + M
β
1 M .
∂M t ∂M β β
Za T > Tc imamo da je i za H ̸= 0 imamo magnetizaciju M = 0, odnosno da je
( )
∂H
= Φ(t, 0) = λ−β Φ(λt, 0),
∂M t
III.. HIPOTEZA SKALIRANJA ZA TERMODINAMIČKE FUNKCIJE 97
1
odaberemo da je λt = 1 tako da imamo χ
= tb , pa imamo:
χ ∼ t−b ∧ χ ∼ t−γ ⇒ b = γ.
1
Za slučaj T < Tc i H = 0 znamo da je M ∼ tβ , odnosno M β = ct i tada je
1 1 1
Φ(t, M β ) = 0, pa zbog Φ(t, M β ) = λ−γ Φ(λt, λM β ) dobijamo
1 1 1
1 ∂Φ(t, M β ) 1 β1 −1 M β −γ ∂Φ(λt, λM β )
=M 1 M = λ 1 .
χT ∂M β β β ∂M β
Uzmemo da je λ|t| = 1, odnosno −λt = 1, uvodimo parametar u = cλt = −c što
1
se još moze zapisati izrazom u = λM β = λct = −c, pa se tako dobija
1 ( )
1 Mβ 1 −γ ∂Φ(−1, u)
= (− ) λ .
χT β t ∂u u=−c
odakle nalazimo
χ ∼ (−t)γ .
γ γ
H = M M β Φ(0, 1) = M 1+ β Φ(0, 1),
γ = β(δ − 1).
98 § 3. FAZNI PRELAZI
gde su parametri t = T
Tc
− 1 temperaturno rastojanje i H primenjeno magnetno
polje.
(a) Pokazati da magnetni kritični eksponent β i δ mogu biti napisani kao funkcije
at i aH .
-
(b) Pokazati da kritični eksponenti γ i γ ′ , odredeni na osnovu ponašanja suscep-
-
tibilnosti ispod i iznad kritične temperature, medusobno jednaki.
Slika 3.7:
39. U Tabeli br. 3.7 su data merenja magnetiizacije nikla. S kolikom greškom ova
-
merenja potvrduju hipotezu skaliranja? Za nikl kritična temperatura je Tc =
627.2K, β = 0.368 i γ = 4.22.
40. U blizini kritične tačke, jednačina stanja magnetika se može modelirati jednačinom
41. (a) f (r) je homogena funkcija u slučaju da je ispunjeno f (λr) = g(λ)f (r). Pret-
postavljajući da je g(λ) diferencijabilna funkcija, pokazati da je g(λ) = λp .
i. Pokazati da su
(A) f (x, y) = x3 + y 2 ,
(B) f (x, y) = x3 y 2 + x2 y 3
generalisane homogene funkcije. Odrediti odgovrajuće parametre a i b.
1 −b 1 −a
f (x, y) = x a f1 (x a y) = y b f2 (y b x). (3.75)
42. Odrediti latentnu toplotu Van der Waals-ovog gasa, uzimajući da je CV poznato
i konstantno na temperaturi T < Tc . Smatrati da se može uzeti ista jednačina
stanja iznad i ispod faznog prelaza.
Rešenje:
100 § 3. FAZNI PRELAZI
q = T (sg − sl ), (3.80)
43. Odrediti kritičnu tačku (pc , Tc , vc ) za van der Waals-ov gas. Rešenje:
Kritičnu tačku van der Waals-ovog gasa odredujemo - na osnovu jednačina
( )
∂p kB T 2a
=− + 3 = 0, (3.82)
∂v T (v − b)2 v
2
∂ p 2kB T 6a
( 2 )T = − 4 =0 (3.83)
∂v (v − b) 3 v
2a kB T
= ,
v 3 (v − b)2
2 2a 6a
− 4 =0
v − b v3 v
Odavde se dobija
4 6
= .
v−b v
Tako da na kraju za kritičnu zapreminu nalazimo
vc = 3b
102 § 3. FAZNI PRELAZI
-
Vraćajući kritičnu zapreminu vc u jednačinu (3.82) odredujemo i kritičnu temper-
aturu
kB Tc 2a 1 2a
= 3 ⇒ Tc = (2b)2
(vc − b)2 vc kB (3b)3
što nam na kraju daje da je kritična temperatura:
8a
Tc =
27kB b
44. Odrediti zavisnost razmene toplotne energije koja prodje kroz jediničnu površinu
neke tečnosti od temperature, razmatrajući problem slobodne površine tečnosti
kao Karnoov ciklus.
tj. da je promena Gibbsove slobodne energije u jednoj i drugoj fazi jednaka. Tako
se dobija da je
−Ss dT − Ms dH = −Sn dT − Mn dH , (3.85)
odakle dobijamo da je
dH Ss − Sn
=− . (3.86)
dT Ms − Mn
U opštem slučaju važi sledeća relacija5
⃗
B ⃗ + M),
⃗ = µ0 (H (3.87)
V
gde je V zapremina sistema. U slučaju superprovodne faze
⃗s = 0 ,
B (3.88)
pa iz (3.87) sledi
Ms = −V H . (3.89)
RT
47. Posmatrajmo jedan mol realnog gasa, čija je jednačina stanja oblika p = (v−b)eaRT v
.
-
48. Posmatrati gas čija su stanja odredena jednačinom
nRT [− na ]
p= e RT v , (3.94)
v − nb
Slika 3.8:
T −Tλ (p)
gde je relativno rastojanje temperature τ = Tλ (p)
od λ-faznog prelaza Tλ (p).
Hemijski potencijal se menja u skladu sa relacijom
1
µ(T, p) = µo (T, p) + α(p)Tλ (p)τ 2 ln τ + .... (3.96)
2
52. N
a. Pokazati da se nad der Waals-ova jednačina može biti napisana u obliku re-
dukovanog izraza:
8T̃ 3
P̃ = − 2, (3.97)
3ṽ − 1 ṽ
p V T
gde je P̃ = pc
, ṽ = Vc
i T̃ = Tc
.
53. N
(a) Pokazati da se nad der Waals-ova jednačina može biti napisana u obliku
redukovanog izraza:
8T̃ 3
P̃ = − 2, (3.98)
3ṽ − 1 ṽ
p V T
gde je P̃ = pc
, ṽ = Vc
i T̃ = Tc
.
(b) Za T < Tc ,van der Waals-ova jednačina prikazuje tečnost-gas faznih prelaza.
Termodinamičke granice za ovaj fazni prelaz je odreden- izrazom µl (P, T ) =
µg (P, T ) (ovde, i u daljem delu teksta nećemo koristiti onaku tilde). Koristeći
sledeće uslove
∫ vg
P dv = pf.p. (vg − vl ). (3.99)
vl
PV
a. Odrediti odnos N kB T
na kritičnoj tački za gas koja opisuje tu jednačinu. Daj
numeričku vrednost tih triju figura.
N kb T N
P = [1 + B(T )] (3.101)
V V
- u zapremini v.
a. Odrediti distribuciju verovatnoće fn da se n molekula nade
b. Izračunati srednju vrednost < n > i disperziju < (n− < n >)2 >.
d. Pokazati da u limesu p = v
V
→ 0, V → ∞, sa V → ∞ sa N
V
= const, f (n) je
oblika Poisson-ove raspodele.
Rešenje:
∑
N ∑
N
N!
f (n) = pn (1 − p)N −n = (p + 1 − p)N = 1. (4.1)
n=0 n=0
n!(N − n)!
∑
N
F (x) = f (n)xn =
n=0
= [(1 − p) + px]N ).
∑N
Polazeći od izraza za srednju vrednost < n >= n=1 nf (n), što se može
odrediti iz prvog izvoda fukcije generatrise po parametru x
dF (x) ∑
N
= nf (n)xn−1 , (4.2)
dx n=1
dF (x) ∑ N
|x=1 = nf (n) =< n > . (4.3)
dx n=1
dF (x)
< n >= |x=1 = pN [(1 − p) − px]N −1 |x=1 = pN. (4.4)
dx
< (n− < n >)2 >=< n2 − 2n < n > + < n >2 >=< n2 > − < n >2 . (4.5)
I.. OSNOVE NERAVNOTEŽNE STATISTIČKE FIZIKE 111
d2 F (x) ∑
N
= n(n − 1)xn−2 f (n) =< n2 > − < n >, (4.6)
dx2 n=2
d2 F (x) d
2
= pN [(1 − p) − px]N −1 =
dx dx
= p2 N (N − 1)[(1 − p) − px]N −2 , (4.7)
N! NN
f (n) = pn (1 − p)N −n ≈ n pn (1 − p)N −n .
n!(N − n)! n (N − n)N −n
Ako za neku vrednost n∗ funkcija f (n∗ ) ima maksimum, tada će i njen loga-
ritam ln f (n∗ ) imati maksimum. Sledi
(N − n)p
(ln f (n))′ = − ln n − 1 + ln(N − n) + 1 + ln p − ln(1 − p) = ln .
n(1 − p)
112 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
1
ln f (n) ≈ ln f (n∗ ) + (ln f (n∗ ))′′ · (n − n∗ )2 ,
2!
∗ ))′′ (n−n∗ )2
f (n) ≈ f (n∗ ) e 2 (ln f (n
1
.
Može se izraziti da je
1 1
(ln f (n∗ ))′′ = − · ,
p(1 − p) N
pa konačno imamo da
(n−n∗ )2
∗ − 2N p(1−p)
f (n) ≈ f (n ) e
pa je f (k) = nk
k!
e−n Poasonova raspodela.
II.. SLUČAJNE ŠETNJE 113
a. Uniformne P (x) = 1
2a
, za −a < x < a,
1 − |x|
b. Laplace-ova raspodela P (x) = 2a
e a,
√ 2
2 x2 − 2a
x
c. Maxwell-ova raspodela P (x) = π a2
e 2.
3. Kretanje čestica u 3D je niz nezavisnih koraka dužine l, gde svaki korak zaklapa
ugao θ za z-osom u skladu sa gustinom verovatnoće P (θ) = π2 cos2 2θ , dok je ugao
φ uniformno distribuiran izmedu- 0 i 2π. Uzeti da se čestice kreću od centra i
(a) Odrediti očekivane vrednosti < x >, < y >, < z >, < x2 >, < y 2 >, < z 2 > i
kovarijanse < xy >, < xz > i < yz >.
∂n ∑ τ i−1 ∂ i n ∑
∞ ∞
∂
+ i
= (− )i [Di (x, t, τ )n(x, t)], (4.10)
∂t i=2
i! ∂t i=1
∂x
Mi (x,t,τ )
gde je Di (x, t, τ ) = i!τ
.
114 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
Rešenje:
∂n(x, t) ∂ 2 n(x, t)
=D , (4.12)
∂t ∂x2
∂n(x, t)
n(x, 0) = n0 i − D = τ,
∂x
i granični uslov:
n(∞, t) = n0 ,
n − n0 D x
θ= ,η = t 2,ϕ = ,
n0 l l
∂θ(ϕ, η) ∂ 2 θ(ϕ, η)
= ,
∂η ∂ϕ2
116 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
θ(ϕ, 0) = 0,
θ(∞, η) = 0,
∂θ(0, η)
= −α,
∂ϕ
τl
gde je α = Ako funciju θ(ϕ, η) prepišemo kao inverznu Laplace-ovu transfor-
Dn0
.
∫
maciju po η: θ(ϕ, η) = θ̄(ϕ, s)e−sη ds i vratimo u prethodnu jednačinu dobijamo
sledeći jednačinu za θ̄(ϕ, s):
∂ 2 θ(ϕ, s)
= sθ̄(ϕ, s).
∂ϕ2
θ̄(ϕ, 0) = 0,
θ̄(∞, s) = 0,
∂ θ̄(0, s) α
=− .
∂ϕ s
3
c1 (s) = αs 2 , c2 (s) = 0.
Kada je fluks negativan, tada se Sands-ovo vreme može posmatrati kao bezdi-
menziona koncentracija:
π n0 2
ts = D( ) .
4 τ
Posle Sands-ovog vremena koncentracija je nula.
p 1-p
x-a x x+a x
Slika 4.1:
6. Čovek hoda duž prave linije, praveći korak dužine a svakih τ sekundi. Svaki
korak je napravljen unapred, sa verovatnoćom p, ili unazad sa verovatnoćom
q = 1 − p.
a. Za veliki broj koraka, odrediti Gausovu raspodelu n(x, t), tako da je n(x, t)dx
- x i x + dx u trenutku t, uzimajući u
je vervatnoća nalaženja u tački izmedu
obzir da je žetač krenuo iz x = 0 u trenutku t = 0.
Rešenje:
118 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
N!
PN (N+ ) = pN+ (1 − p)N −N+ . (4.13)
N+ !(N − N+ )!
Np N (1 − p) 1 1 N
lnPN (N+ ) = N+ ln + (N − N+ )ln + ln( ).
N+ N − N+ 2 2π N+ (N − N+ )
∂lnPN (N+ )
= 0 ⇒ Nmax = N+ = N p
∂N+
∂ 2 lnPN (N+ ) 1
=−
∂N+2
N p(1 − p)
1
n(x, t) = PN (N+ ) dx.
2a
2a 2
4π 2p(1−p)a
τ
t
t (2p−1)a 2p(1−p)a2
Koristimo sledeće relacije: N = τ
, v = τ
iD = τ
, pa dobijamo
sledeću Gauss-ovu raspodelu:
1 1 (x−vt)
2
n(x, t) = √ e− 4 Dt ,
4πDt
gde je D difuziona konstanta a v brzina pomeranja.
∫ ∞ ∫ ∞ √
1 (x−vt)2
− 41 Dt 1 π
n(x, t)dx = √ e dx = √ 1 = 1
−∞ 4πDt −∞ 4πDt 4Dt
∫ ∞
n(x, t)dx = 1
−∞
1 x−vt
√
uvedemo smenu z = 2 Dt
pa imamo:
√ ∫
8Dt ∞ 2 −z2 dz
ydm = z e .
4π −∞ 2
∫∞ √
z 2 e−z
2 dz π
Integral I = −∞ 2
rešavamo smenom ω = z 2 i dobijamo I = Γ( 23 ) = 2
.
Pa imamo da je:
√
ydm = 2Dt.
1 a2 a2
D = 2( − λ2 a2 ) ≈ ,
4 τ 2τ
2( 12 + λa) − 1 2λa2
v= a= .
τ τ
8. Razmotrimo jedan jednostavan model po kome se cene neke akcije može menjati
svakih τ sekundi. Ta cena se može povećavati i da uzima vrednost zn+1 = (1 +
α) ∗ zn u n + 1-vom trenutku vremena. A cena te akcije se može i smanjivati
zn+1 = (1 − β) ∗ zn , gde su α i β < 1 pozitivni unapred parametri. Početna
vrednost akcije je zo u trenutku t = 0.
c. Pokazati da se srednja vrednost cene z, menja po zakonu < log z >=< log zo >
√
+vt, i gde je ∆(log z) = σ t i odrediti v i σ u modela koji je dat gore.
1 (x − a)2
P (x)dx = √ exp(− ). (4.16)
2πσ 2 2σ 2
122 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
10. Verovatnoća pojedinačnog koraka jedno dimenzione slučajne šetnje za korake čije
- d − a < x < d + a je uniformna (a < d). Posle N koraka
su dužine izmedu
pomeranješetača je dato izrazom S = X1 + X2 + · · · + XN , gde je Xi pomeraj
-
šetača tokom i-tog koraka. Pretpostavljajući da su koraci medusobno nezavisni,
kolika je gustina verovatnoće da se posle N koraka šetač pomerio za s. Pored
toga odrediti kolika je srednja vrednost pomeraja < s >, a kolika je standardna
devijacija šetača.
Prvi član opisuje brzinu pomeraja gustine v, dok drugi član −D ∂n(x,t)
∂x
opusije
difuziju čestica u suprotnom smeru od gradijenta koncentracije.
12. Čovek hoda slučajno duž ulice, koračajući koracima dužine a unapred i unazad,
svakih τ sekundi. Ako je x(t = 0) = x0 > 0, odrediti verovatnoću P (t) da se
- na mestu x = 0 u trenutku t.
čovek nade
Rešenje: Za a << x0 potrebno je mnogo više koraka. U jednom od prethodnih
no(x+xo,t)
no(x-xo ,t)
n (x ,t)
-xo 0 xo x
Slika 4.2:
1 ∂n ∂2n a2
zadataka za p = 2
je dobijeno ∂t
= ∂x2
D, gde je D = 2τ
. Moramo uzeti u obzir
da je početni uslov
n(x, 0) = δ(x − x0 ),
1 x2
n0 (x, t) = √ e− 4Dt .
4πDt
∫ ∞ ∫ ∞
P (t) = 1 − n0 (x − x0 , t)dx +
n0 (x + x0 , t)dx =
∫ ∞ 0
∫ ∞ 0
1 (x−x )2
− 4Dt0 1 (x+x0 )2
=1− √ e dx + √ e− 4Dt dx.
4πDt 0 4πDt 0
-
Uvodenjem smene y = x − x0 dobijamo:
∫ ∞ ∫ ∞
1 y2
− 4Dt 1 y2
P (t) = 1 − √ e dy + √ e− 4Dt dy.
4πDt −x0 4πDt x0
y2
Ako jedinicu prepišemo kao normu Gausijana i izvršimo smenu z 2 = 4Dt
dobi-
jamo: ∫ ∫
−z0 ∞
1 −z 2 1
e−z dz.
2
P (t) = √ e dz + √
π −∞ π z0
Na kraju dobijamo
2 2 x0
P (t) = √ erf c(z0 ) = √ erf c( √ ),
π π 4Dt
gde je erfc funkcija komplementarna error funkciji. Verovatnoća ”preživljavanja”
je prikazana na Slici ??.
x = l1 + l2 + . . . + lt .
II.. SLUČAJNE ŠETNJE 125
Tako imamo:
∑
N ∑
N ∑
N
Pt (x) = ... δx,l1 +l2 +...+lt p(l1 )p(l2 ) . . . p(lt ). (4.17)
l1 =−N l2 =−N lt =−N
1 ∑N
2πinx
Pt (x) = ft (n)exp(− ).
2N + 1 n=−N 2N + 1
1 ∑N
2πn
p(x) = λ(n)exp(−i x)
2N + 1 n=−N 2N + 1
2πx
Zamenjujući (4.18) u (4.17) i množenjem sa exp(i 2N +1
n) dobijamo:
∑
N
2πx ∑
N ∑
N
2πx
Pt (x) exp i n = δx,l1 +l2 +...+lt exp i np(l1 )
x=−N
2N + 1 x=−N x=−N
2N + 1
∑
N
ft (n) = δx,l1 +l2 +...+lt λ(n)t → λ(n)t
x=−N
U termodinamičkom limesu je
∑N ∫ N ∫
1 1 1
= dn ≈ πdφ,
2N + 1 n=−N 2N + 1 −N 2π −π
126 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
imamo
∫ π
1
Pt (x) = ft (φ)e−ixφ dφ
2π −π
∑∞
ft (φ) = Pt (x)eixφ
x=−∞
da n sijalica budu upaljena u trenutku t) objasniti zašto ova evolucija pokazuje os-
obine Markov-ljevog lanca. Odrediti matricu prelaza, i pokazati da postoji jedin-
stveno stacionarno stanje, koje je dato binomalnom distribucijom.
Rešenje:
Neka je n(t) broj svetla koji je uključen u trenutku t. Pošto je to slučajna promenljiva,
n
onda je verovatnoća da se isključi jedno svetlo jednako N
, dok je verovatnoća da
N −n
se uključi jedno svetlo N
. Onda je uslovna verovatnoća data izrazom
n N −n
P (n, t|n′ , t + 1) = δn,n−1 + δn′ ,n+1 , (4.19)
N N
gde je n′ - broj uključenih svetijkiu trenutku t+1. Ova uslovna verovatnoća karak-
teriše Markovljeve slučajne procese, pošto je vrednost slučajne promenljive za-
vise samo od vrednosti slučajne promenljive u pretohdnom trenutku, odnosno
P (n(0), 0; n(1), 1; . . . ; n(t), t|n(t + 1), t + 1) = P (n(t), t|n(t + 1), t + 1). Ako sa
W (n, t) obeležimo verovatnoću da se trenutku t bude uključeno n sijalica, onda
je verovatnoća da u trenutku t + 1 bude uključeno n′ sijalica data izrazom
∑
N
′
W (n , t + 1) = W (n, t)P (n, t|n′ , t + 1). (4.20)
n=0
w(t + 1) = w · P, (4.21)
gde je
n N −n
Pnn′ = δn′ ,n−1 + δn′ ,n+1 , (4.23)
N N
128 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
∑
N
N!
1 = w0 , (4.27)
n=0
n!(N − n)!
-
odakle se vidi da je w0 = 1, medutim prethodna jednačina važi ako se uvede neki
pomoćni parametar z
∑
N
N!
1 = w0 z n (1 − z)N −n , (4.28)
n=0
n!(N − n)!
III.. MARKOVLJEVI PROCESI 129
-
pošto je ovako uvedena suma ima istu konačnu vrednost kao i pre uvodenja
paramtra z:
∑
N
N!
z n (1 − z)N −n = (z + (1 − z))N = 1. (4.29)
n=0
n!(N − n)!
15. Razmotriti dve kutije A i B u kojima se nalaze tri svetle i dve tamne kuglice.
U kutiji A mogu da stanu dve kuglice, dok u kuglice mogu se naći tri kuglice.
-
Prelaz iz jednog stanja u drugo se sastoji u zameni bilo koje dve kuglice izmedu
njih, tako što se jedna kuglica premesti iz kutije A u kutiju B.
(a) Odrediti matricu prelaza Q i matricu uslovne verovatnoće P (m, to |n, t).
(a) Odrediti matricu prelaza, njene svojstvene vrednosti, i leve i desne svojstvene
vektore.
130 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
17. Razmotriti asimetričnu slučajnu šetnju na rešetki sa četri čvorišta. Matrične vred-
nosti prelaza su w12 = w43 = 1, w23 = w34 = 34 , w23 = w34 = 1, w21 = w32 = 1 i
wij = 0 za sve ostale prelaze.
(c) Napisati P (n, t) za slučaj P (n, 0) = δn,1 . Kolika je vrednost P1 (2, t)?
(b) Ako šetač startuje na čvorištu 3 u trenutku 0 odrediti srednje vreme prvog
prolaska.
IV.. LANGEVIN-OVA JEDNAČINA 131
- matricom
20. Razmotriti Markov-ljev proces koji je odreden
a 0 0
T̂ = 12 (1 − a) 0 1 . (4.30)
1
2
(1 − a) 1 0
gde je η(t) Beli šum sa ⟨η(t)⟩ = 0 i ⟨η(t)η(t′ )⟩ = δ(t − t′ ). Odrediti ⟨x2 (t)⟩ i ⟨x(t)⟩
ako je ẋ(t = 0) = 0 i x(t = 0) = 0.
Rešenje:
Uradimo Furier-ov transfrom ẋ :
∫ ∞ ∫ ∞
1 1 1
√ dtẋ e
|{z}
iωt
= √ xeiωt |∞
−∞ − √ dtx(iω)eiωt ,
2π −∞
dx
2π 2π −∞
xe−st |t=∞
t=0 − s x(t)e−st dt =
0
−x(t = 0) + sX(s).
Tako dobijamo:
4 g g
2 3 3 +w
2
g
- w 4
2
-w w
g
- 2 +w 1
2 -g 1 g
2 2 - -w
2
Slika 4.3:
∫ ∞ ∫ ∞
′ dω dω ′ ′ ′
⟨x(t)x(t )⟩ = √ √ ⟨x(ω)x(ω ′ )⟩e−iωt e−iω t ,
−∞ 2π −∞ 2π
Pa kada zamenimo pretposlednu relaciju u poslednju i prointegralimo po ω ′ do-
bijamo:
∫ ∞ ′
′1 e−iω(t−t )
⟨x(t)x(t )⟩ = dω .
2π −∞ P (ω)P (−ω)
134 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
√
Nule P (ω) su ω1,2 = −i γ2 ± ω02 − ( γ2 )2 , a za nule P (−ω)imamo ω3,4 = i γ2 ±
√
ω02 − ( γ2 )2 .
γ2
U prvom slučaju imamo da je ω02 > 4
, tada su rezidijumi rasporedjeni u kom-
pleksnoj ravni postavljeni na način kao što je prikazano na Slici 4.3a, dok su u
γ2
drugom slučaju, kada je ω02 < 4
, rezidijumi su rasporedjeni duž imaginarne ose
(slika 4.3b).
1 1 11
⟨x2 (t)⟩ = −i( + ) = .
2ω ′ (−iγ)(2ω ′ − iγ) (−2ω ′ )(−iγ)(−2ω ′ − iγ) γ ω0
22. Na materijalnu česticu deluje elastična sila duž ose x, tako da se Langevin-ova
jednačina može zapisati relacijom
d2 x dx 1
2
+γ + ωo2 x = R(t), (4.33)
dt dt m
⃗
R(t)
⃗u˙ = −γ⃗u + , (4.34)
m
⃗
gde je R(t) sila kojom deluje fluid na tu Brown-ovu česticu. Srednja vrednost te
⃗
sile < R(t) ⃗ R(t
>= 0, dok je autokorelaciona funkcija data izrazom < R(t) ⃗ ′ ) >=
IV.. LANGEVIN-OVA JEDNAČINA 135
gde beli šum η ima srednju vrednost < η(t) >= 0 i korelacionu funkciju <
η(t)η(t′ ) >= Γδ(t − t′ ). Svi ostali viši momenti šuma uzeti da su jednaki nuli.
(b) Odrediti kako četvrti < v(t)4 > i šesti moment < v(t)6 > zavisi od vremena.
Odrediti asimptotske vrednosti ovih momenata za t → ∞.
25. Iz otvora koji se nalazi na visini h iznad podloge, u homogenom polju sile Zemljine
teže ⃗g , ispuštaju se male čestice mase m bez početne brzine. Pretpostavljajući da
se kretanje tih četsica može opisati Langevin-ovom jednačinama
d⃗v
m = −mγ⃗v + f⃗(t) + m⃗g , (4.36)
dt
za koje važi da je < f⃗ >= 0, dok je < f⃗(t)f⃗(t′ ) >= 6mγkB T δt,t′ . Oceniti radijus
kruga na podlozi u koji upadaju spomenute čestice dodu - do podloge, znatno
gde je ζ(t) Gauss-ova slučajna varijabla, sa < ζ >= 0 i < ζ 2 >= 1. Zapaziti
da greška u diskrtizaciji je O(∆t), što je preciznije od standardne diferencijalne
kB T
jednačine O(∆T 2 ). U prethodnom izrazu je uvedena konstanta B = α
. Radi
jednostavnosti uzeti da je kB = 1 i α = 1.
x2
a. Napisati program koji uključuje jednačinu (4.38), uzeti da je V (x) = 2
)).
Voditi računa da je slučajna varijabla Gausova varijabla.
c. Nacrtati numeričke izračunati < x(t)2 >, zatim nacrtati ananlitičko rešenje
koji si dobio pod (a), i uporedi te dve funkcije. To uradi za nekoliko različitih
vrdnosti temperature T . Ovaj način omogućuje da podesis vremenski korak
∆t.
Langevin-ova jednačina:
dv γ 1
= − v − ωo2 x + ξ(t), (4.39)
dt m m
√
k
gde je ωo = m
. Ako su xo i vo početni položaj i početna brzina Brown-ove
čestice, i ako pretpostavimo da je u ravnoteži sa fluidom pokazati da se stanje
procesa može uzeti stacionarnim ako je g = 4γkB T . Odrediti korelacionu funkciju
m
brzine << v(t2 )v(t1 ) >ξ >T . Uzeti da važi ekviparticiona teorema < vo2 >T =
2
1 -
k T i da je vibraciona potencijalna energija < 12 ωo2 < x2o >T >= 21 kB T. Takode
2 B
gde je η(t) sila kojom deluje fluid na tu Brown-ovu česticu. Srednja vrednost te
sile < η(t) >= 0, dok je autokorelaciona funkcija data izrazom < η(t)η(t′ ) >=
Γδt,t′ . Odrediti srednju vrednost < v(t)2 > i < x(t)2 >. Pretpostaviti x(0) = 0 i
v(0) = 0.
Rešenje:N S
dv
∫t
Pošto je dt
= η(t), i koristeći početne uslove v(0) = 0, dobijamo: v(t) = 0
dt η(t).
Odatle, koristeći linearnost očekivane vrednosti dobijamo:
∫ t ∫ t
⟨v (t)⟩ =
2
dt dt′ ⟨η(t)η(t′ )⟩ = Γt
0 0
-
29. Razmatrati polimer, čiji monomeri se mogu smatrati bezmasenim, i koji su medusobno
povezani jedan za drugi.
gde je ⟨ηn (t)ηn (t′ )⟩ = Γδnn′ δ(t − t′ ). Odrediti ⟨(xn (t) − xn (t = 0))2 ⟩.
Rešenje:
Radi jednostavnosti razmotriémo slučaj kada je xn (t = 0) = 0. Uvesti Fourier-ov
transform:
∫
1
xn (t) = √ x(h, t)eihn dh
2π
sama jednačina je onda:
∫ ∫
1 1
√ ikn
ẋ(k, t)e dk = √ x(h, t)eik(n+1) dk +
2π 2π ∫
1
+ √ x(h, t)eik(n−1) dk
2π ∫ ∫
1 1
− 2√ x(k, t)e dk + √
ikn
η(k, t)eikn dk.
2π 2π
Odavde se dobija da je
dx
− a(k)x(k, t) + η(k, t),
dt
⟨η(⃗r, t)⟩ = 0
Izračunati korelacionu funkciju visine deponovanih čestica ⟨h(⃗r, t)h(⃗r′′ , t′′ )⟩.
Rešenje:
Izvršimo Fourier-ovu transformaciju funkcije visine deponovanih čestica h(⃗r, t) =
∫∫ 3 ⃗
1
(2π)2 d ⃗kdteik⃗r eiωt h(⃗k, ω), tako da Langevin-ova jednačina može da se prepiše
kao:
iωh(⃗k, ω) = −k 2 h(⃗k, ω) + η(⃗k, ω).
Koristeći izraz koji sledi iz jednačine kretanja i gore dobijenog izraza za korelaciju
imamo:
(2π)2 δ(⃗k + k⃗′ )δ(ω + ω ′ )
⟨h(⃗k, ω)h(k⃗′ , ω ′ )⟩ = ,
(k 2 + iω)(k ′2 + iω ′ )
pa kada ovaj izraz vratimo u integral imamo:
∫ ∫
′ 1 3⃗ r −r⃗′ ) iω(t−t′ )
i⃗k(⃗ 1
⟨h(⃗r, t)h(r⃗ , t )⟩ =
′ d kdωe e .
(2π)2 (k 2 + iω)(k ′2 + iω ′ )
Singulariteti se nazlaze na vrednostima ω1 = −ik 2 i ω2 = +ik 2 .
(2π)2 δ(⃗k + k⃗′ )δ(ω + ω ′ )
⟨h(⃗k, ω)h(k⃗′ , ω ′ )⟩ = ,
(k 2 + iω)(k ′2 + iω ′ )
pa kada ovaj izraz vratimo u integral imamo:
∫ ′
′ 1 r −r⃗′ )
3⃗ i⃗k(⃗ eiω(t−t )
⟨h(⃗r, t)h(r⃗′ , t )⟩ = d ke 2πRes ω2 =ik 2 ( )
(2π)2 (k 2 + iω)(k ′2 + iω ′ )
Kada sračunamo rezidum i uvedemo smenu R ⃗ = ⃗r − r⃗′ i τ = t − t′ dobijamo:
∫ −k2 τ
′ 1 ⃗⃗ e
⟨h(⃗r, t)h(r⃗′ , t )⟩ = d3⃗keikR 2
4π k
Sada razdvajamo integral na normalnu i paralelnu komponentu po k:
∫ ∫ ∞
′ 1 ⃗ −k 2 |τ | eik∥ R −k∥2 τ
⟨h(⃗r, t)h(r , t )⟩ =
⃗′ 2
d k⊥ e ⊥ dk∥ 2 2
e .
4π −∞ k∥ + k⊥
Posle integracije dobijamo:
∫ ⃗ ∫
′ 1 2 τ cosh(|k⊥ |R)
2 ⃗ −2k⊥ 1
dk⊥ e−2k⊥ τ cosh(|k⃗⊥ |R),
2
⟨h(⃗r, t)h(r , t )⟩ =
⃗′ d k⊥ e =
4π k⊥ 4
gde smo izvršili integraciju po polarnom delu cilindričnih koordinata. Uvodimo
2|τ |
smenu |k⃗⊥ |R i α = R2
,
pa dobijamo izraz:
∫
⟨h(⃗r, t)h(r⃗ , t )⟩ = R dzcosh(z)e−αz .
′ 2
′
∫
U slučaju kada R → 0, tada dobijamo da α → ∞ i cosh z ∼ 1, pa je onda dreαr ∼
2
√1 odnosno imamo:
α
R2
⟨h(⃗r, t)h(r⃗′ , t′ )⟩ ∼ .
|τ |
142 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
gde beli šum η ima srednju vrednost < η(t) >= 0 i korelacionu funkciju <
η(t)η(t′ ) >= Γδ(t − t′ ). Svi ostali viši momenti šuma uzeti da su jednaki nuli.
(b) Odrediti kako četvrti < v(t)4 > i šesti moment < v(t)6 > zavisi od vremena.
Odrediti asimptotske vrednosti ovih momenata za t → ∞.
gde beli šum η ima srednju vrednost < η(t) >= 0 i korelacionu funkciju <
η(t)η(t′ ) >= Γδ(t − t′ ). Svi ostali viši momenti šuma uzeti da su jednaki nuli.
(b) Odrediti kako četvrti < v(t)4 > i šesti moment < v(t)6 > zavisi od vremena.
Odrediti asimptotske vrednosti ovih momenata za t → ∞.
φ̈ + γφ = M (t). (4.43)
Odrediti korelacionu funkciju x-komponenti dipolnog momenta molekula < cos φ(t1 ) cos φ(t2 ) >.
V.. FOKKER-PLANCK-OVA JEDNAČINA 143
V. Fokker-Planck-ova jednačina
1 (x−a(t))2
P (x, t) = √ e− 2σ(t) , (4.44)
2πσ(t)
gde je a(t) =< x(t) > srednja vrednost slučajnog položaja x i varijansa σ(t) =<
x2 (t) > − < x(t) >2 računata za istu slučajnu promenljivu. Raspodela verovatnoća
se menja u skladu sa jednačinom Smoluchovski-og
∂P (x, t) 1 ∂ ∂2
+ (F (x)P (x, t)) − D 2 P (x, t) = 0. (4.45)
∂t mγ ∂x ∂x
36. Rešiti Foker-Plankovu jednačinu braunovske čestice koja se kreće u polju, poten-
cijala U (x) = ax2 , sa početnim položajem u tački xo .
1 (x −x1 )2
− 4D(t
2
P (x1 , t1 |x2 , t2 ) = √ e 2 −t1 ) , (4.46)
4πD(t2 − t1 )
1 x2
W (x, t) = √ e− 4Dt , (4.47)
4πDt
144 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
⃗ o e−iωo t+iφ(t) ,
⃗ =E
38. Jednačina promene elektromagnetskog polja je data jednačinom E
gde je φ(t) Wiener-ov stohastični proces, koji je dat difuzionom jednačinom sa
konstantom difuzije D. Odrediti korelacionu funckiju G(t1 −t2 ) =< E(t1 )E ∗ (t2 ) >
∫∞
i spektar snage Iω = −∞ dteiωt G(t).
Rešenje:
Korelaciona funkcija je data jednačinom
G(t1 − t2 ) =< E(t1 )E ∗ (t2 ) >= |E0 |2 e−iω(t1 −t2 ) < ei(φ(t1 )−φ(t2 )) > . (4.48)
∫∞
Koristeći da se može preći u sledeći integralni zapis oblika φ(t2 )−φ(t1 ) = −∞
φ(t)(δ(t−
t2 ))−δ(t−t1 ))dt funkcija generatrise za Wiener-ove procese data jednačinom, gde
radi jednostavnosti se može uzeti da je g(t) = (δ(t − t2 )) − δ(t − t1 )), odakle se
dobija da je
∫∞ ∫∞
< ei −∞ g(t)φ(t)dt
>= e− −∞ dt1 dt2 g(t1 )g(t2 )min(t1 ,t2 )
= e−D|t2 −t1 | . (4.49)
39. Izračunati spektar snage Iv,v (ω), za harmonijski vezanu Braunovsku česticu. Nacr-
γ
tati korelacionu funkciju Cv,v (t) i Iv,v (ω) za slučaj ωo > m
.
V.. FOKKER-PLANCK-OVA JEDNAČINA 145
40. Razmotrimo stohastičku varijablu Y sa tri realizacije y(1), y(2) i y(3). Razmot-
- dva stanja Q1,1 (s) = Q2,2 (s) = Q3,3 (s) = 0,
rimo tranzicionu verovatnoću izmedu
Q1,1 (s) = Q2,2 (s) = Q3,3 (s) = cos2 ( 2πs
3
) i Q1,1 (s) = Q2,2 (s) = Q3,3 (s) = sin2 ( 2πs
3
).
Ako je inicijalno sistem u stanju y(1), koja je verovatnoća da se nade - stanje y(2)
Pa imamo:
√ √
Pa imamo a± = 1
2D
(v ± v 2 + 4sD i b± = 1
2D
(v ± v 2 + 4sD, gde zbog konver-
gencije odbacujemo a− i b+ , i uzimamo nove oznake gde su a = a+ i b = b− .Dalje
zbog neprekidnosti funkcije p(x, s) imamo sledeći uslov:
p(x> ) = p(x< ) ⇒ A = B.
dp(x, s) d2 p(x, s)
sp(x, s) − δ(x) = −v +D ,
dx dx2
∫ϵ ∫ϵ
sa −ϵ
, i primetimo da je −ϵ
p(x, s)dx = 0 kao integral neparne funkcije na simetričnim
granicama, pa dobijamo:
dp(x, s) dp(x, s)
−1 = −v(p(ϵ, s) − p(−ϵ, s)) + D( |ϵ − |−ϵ ).
dx dx
1
0− = Bb − Aa,
D
VI.. MASTER JEDNAČINE 147
v 2 + 2sD
43. Pokazati da master jednačina uslovljava stalni rast entropije sistema pri bilo kojoj
početnoj neravnotežnoj raspodeli.
Rešenje:
Entropija sistema je uvedena
∑
S = − < ln ωi >= − ωi log ωi . (4.52)
i
44. Dokazati da pri bilo kojoj početnoj raspodeli ωi (0) vodi do relaksacije stanja.
Rešenje:
ω1
ω2
Uvedimo vektor stanja W =
..
, i matricu prelaza sa njenim komponen-
.
ωn
∑
tama Lij = Pij − δij k Pik . Projektovanjem jednačinu po bazisu funkcija čistih
L† = L. (4.57)
-
Ako se uvedu bazis, sa medusobno ortonormiranim vektorima φi , tako što je
Lφn = Ln φn , (4.61)
onda dobijamo da je
∑ ∑
exp[tL] Cn φn = Cn exp[tLn ]φn . (4.64)
n n
ψ1
ψ2
U slučaju da odaberemo proizvoljni operator ψ =
..
. Ispitajmo znak izraza
.
ψn
(ψ, Lψ), on je jednak
∑
(ψ, Lψ) = ψi∗ (Pij ψj − Pij ψi )
i,j
1∑
= − Pij |ψi − ψj |2 ≤ 0,
2 ij
φo prirpada skupu svojstvenih vektora. Ako izdvojimo iz sume (4.59) sabirak koji
odgovara tom svijstvenom vektoru
∑
W (t) = C0 φ0 + Cn exp[tLn ]φn . (4.65)
n̸=0
∂Pn (t)
= gn−1 Pn−1 (t) + rn+1 Pn+1 (t) − (gn + rn )Pn (t) (4.66)
∂t
46. Veliki broj uzoraka radioaktivnih uzoraka sa različitim brojem atoma. Radioat-
ivni materijal koja se raspada sa koeficientom raspada α.
- susednih čvorišta sa
47. Čestice se pomeraju po jednodimenzionoj rešetki izmedu
verovatnoćom α u jednom smeru, a sa verovatnoćom β u drugom smeru (u je-
-
dinici vremena). Napisati master jednačinu koja odreduje evoluciju skupa verovatnoća
- na n-tom čvorištu (gde n ∈ {−∞, ..., ∞} u trenutku t).
pn (t) da se čestica nade
Uzeti da se čestica u trenutku t = 0 nalazi na čvorištu n = 0. Odrediti < n >t u t
√
i standardnu devijaciju ∆n = < n2 >t − < n >2t .
Rešenje:
Čestica sa verovatnoćom αpn−1 (t)δt prelazi n − 1-vog čvorišta rešetke na susedno
čvorište u vremenskom intervalu (t, t + δt). U suprotnom smeru, sa čvorišta n + 1
na čvorište n, čestica preskoči sa verovatnoćom βpn+1 (t)δt u intervalu (t, t + δt).
Te verovatoće odreduju - koliki će biti doprinos čestica sa susednih čvorišta na n-
to čvorište, dok verovatno ca (α + β)pn δt predstavlja kolika će biti verovatnoća
”odliva” čestica. Master jednačina je onda oblika
∑
∞
F (x, t) = pn (t)xn , (4.68)
n=−∞
152 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
∂F (x, t) ∑∞
dpn (t) n
= x , (4.69)
∂t n=−∞
dt
gde ćemo verovatnoć pn (t) prepoznati kao koeficijente pored članova sa xn , zbog
čega uvodimo parametar n = k − l, tako da funkcija generatrise postaje
∑
∞ ∑
∞
αl+n β l tn+2l
F (x, t) = exp[−(α + β)t] xn , (4.73)
n=0 l=0
(n + l)!l!
odakle dobijamo da je
∑∞
αj β j t2j
−n
pn (t) = exp[−(α + β)t](βt) , (4.74)
j=0
(j − n)!j!
- izrazom
dok je drugi moment raspodele odreden
( )
∂ ∂F
2
< n (t) >= x |x=1 = (α − β)2 t2 + (α + β)t. (4.76)
∂x ∂x
∂Pn (t)
= (n + 1)Pn+1 (t) + (n − 1)Pn−1 (t) − 2nPn (t). (4.77)
∂t
n
Uzimajući da je početna raspodela Poisson-ovog oblika Pn (0) = e−µ µn! , odrediti
verovatnoću preživljavanja populacije do vremena t, odnosno funkciju S(t) =
- srednju vrednost kvadrata broja jedinki populacije
1 − Pn=0 (t). Odrediti takode
< n2 (t) >.
Rešenje:
Prvo ćemo odrediti funkciju generatrise
∑
G(z, t) = Pn (t)z n . (4.78)
n
∂G(z, t) ∑ ∑ ∑
= nPn (t)z n−1 + nPn (t)z n+1 − 2nPn (t)z n , (4.80)
∂t n n n
∂G(z, t) ∑ ∂z n ∑ ∂z n ∑ ∂z n
= Pn (t) + z2 Pn (t) − 2z Pn (t) =
∂t n
∂z n
∂z n
∂z
∂G
= (1 − z)2 .
∂z
154 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
-
Rešenje prethodne parcijalne diferencijalne jednačine se može svesti na odredivanje
funckije generatrise koja je funkcija samo jedne promenljive g(t̃). Onda se difer-
encijalna jednačina po G, može napisati u obliku
dg(t̃) ∂ t̃ dg(t̃) ∂ t̃
= (1 − z)2 . (4.81)
dt̃ ∂t dt̃ ∂z
Uzimajući u obzir da je zavisnost g(t̃) nije trivijalna zavisnost1 , dobija se
∂ t̃ ∂ t̃
= (1 − z)2 , (4.82)
∂t ∂z
1
Time smo dobili i funkciju generatrise g( 1−z ) = e−µ(1−z) , koja se može trans-
formisati u
µ
g(x) = e− x , (4.84)
-
a što može koristiti za odredivanje funkcije generatrise
1
) = e−µ 1+t(1−z) .
1−z
G(z, t) = g(t + (4.85)
1−z
µ
S(t) = 1 − Pn=0 (t) = 1 − e− 1+t . (4.87)
-
Srednja vrednost kvadrata broja jedinki se odreduje izrazom
∑
< n2 (t) >= Pn (t)n2 (t). (4.88)
n
∂G(z, t) ∑
|z=1 = nPn (t) =< n(t) >, (4.89)
∂z n
odnosno
∂ 2 G(z, t) ∑
|z=1 = n(n − 1)Pn (t) =< n2 (t) > − < n(t) > . (4.90)
∂z 2 n
∂ 2 G(z, t) ∂G(z, t)
< n2 (t) >= |z=1 + |z=1 . (4.91)
∂z 2 ∂z
3 4
1 2
Slika 4.4:
156 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
- adsorbovan u čvorištima 1, 4 i
(d) Odrediti koju su verovatnoće da se šetač nade
5.
Rešenje:N S
dp1
= αp2 ,
dt
dp2
= −3αp2 + αp3 ,
dt
dp3
= −2αp3 + αp2 ,
dt
dp4
= αp3 ,
dt
dp5
= αp2 .
dt
sp1 = αp2
sp4 = αp3
sp4 = αp2
α
p1 = p 5 = p2
s
s + 2α
=α √ √
5+ 5 5− 5
(s + 2
α)(s + 2
α)s
α
p4 = p3
s
α2
= √ √
5+ 5 5− 5
(s + 2
α)(s + 2
α)s
Inverznim Laplasovim transformom dobijamo:
√ [ √ ] √ [ √ ]
5−1 5− 5 1+ 5 5+ 5
p2 = √ exp − αt + √ exp − αt
2 5 2 2 5 2
( [√ ] [ √ ])
1 5− 5 5+ 5
p3 = √ exp − αt − exp − αt
5 2 2
2
Radi jednostavnosti zapisa, Laplasov transform funkcije označavamo isto kao i samu funkciju.
Iz konteksta i korišćenjih promenljivih (t ili s) vidi se šta je u pitanju.
158 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
( [ √ ] [ √ ])
1 5+ 5 5− 5
p1 = p5 = 2 − exp − αt − exp − αt
5 2 2
( [ √ ] [ √ ])
1 2 5+ 5 2 5− 5
p4 = 1+ √ exp − αt − √ exp − αt
5 5( 5 + 1) 2 5( 5 − 1) 2
(d) Verovatnoće da će se šetač na kraju (posle dugo vremena) naći u tačkama 1,
5. i 4, dobijamo iz formula za p1 (t), p5 (t) i p4 (t) dobijenih u delu zadatka pod
(b), posmatrajući šta se dešava kada posle dugo vremena (tj. kada t → ∞).
One su redom 25 , 2
5
i 15 .
4
a a
b b 4
2 3 2 3
b b b
a a a a
a a
1 1
ṗ1 = pN + p2 − 2p1
ṗ2 = p1 + p3 − 2p2
..
.
ṗN = pN −1 + p1 − 2pN
(s + 2)p1 = pN + p2 − 1
(s + 2)p2 = p1 + p3
..
.
(s + 2)pN = pN −1 + p1
160 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
(j)
gde je ei i-ta komponenta j-tog svojstvenog vektora.
†
Posto je M̂ p⃗ = ⃗i ubacivanjem jedinica dobijamo T̂ † T̂ ⃗
{zT̂ } T̂ p⃗ = i, pa odatle lako
| M̂
M̂diag
VI.. MASTER JEDNAČINE 161
dobijemo izraz:
λ−1
(1) (1) (1) (1)
e1 ··· e1 0 e1 ··· e1 1
1
. .. .. . ..
−1 ⃗
p⃗ = T̂ † Mdiag T̂ i = .. . . .. . 0 ,
(N ) (N ) −1 (N ) (N ) ..
eN · · · eN 0 λN eN · · · eN .
kada se gornji izraz raspiše po komponentima dobija sledeći oblik vektora verovatnoće:
p1 = λ−1 −1 −1
(1)
2 2 (2) 2 (N )
1 (e1 ) + λ2 (e1 ) + . . . + λN (e1 )
p2 = λ−1 −1 −1
(1) (1) (2) (2) (N ) (N )
1 e1 e2 + λ2 e1 e2 + . . . + λN e1 e2
..
.
(s + 2 − λ)x1 = x2 + xN
(s + 2 − λ)x2 = x1 + x3
..
.
(s + 2 − λ)xN = xN −1 + x1
2π
λn = s + 2 − 2cos( n),
N
1 2πn
enm = √ ei N m .
N
1 1 1 1 s+1
p1 = ( + + )= .
3 s+3 s+3 s s(s + 3)
w-
1 2 3
w+ w+
Slika 4.7:
sp1 − 1 = −ω+ p1 + ω− p2
sp2 = ω+ p1 − ω− p2 + ω+ p3
sp3 = ω+ p2 .
∂f ∂f ∂f f − fo
+ ⃗v + F⃗ =− .
∂t ∂⃗r ∂⃗p τ
p2
energija elektrona ϵ(p) = 2m∗
(gde je m∗ efektivna masa elektrona), i da je relak-
saciono vreme funkcija brzine elektrona data izrazom τ = A|v|s gde su A > 0
i s > −7 poznate konstante. Pri tome uzeti da je ∂f
∂t
= 0, i da je ∂f
∂⃗
r
= 0, dok se
f −fo
kolizoni integral svodi na τ
, gde je fo ravnotežna raspodela, a τ je relaksaciono
vreme.
Rešenje:
Iz uslova stacionarnosti i homogenosti funkcije raspodele primenjene za slučaj
idealnih gasova, dobijamo da je Boltzmann-ova jednačina oblika:
∂f f − fo
−eE
⃗ =− . (4.92)
∂⃗p τ
Pošto nije potrebno razmatrati neravnotežnu raspodelu elektrona f od impulsa
p⃗, prelazimo na zavisnost od ϵ, što je moguće pošto radimo u klasičnom slučaju.
Tako se dobija izraz:
⃗ v ∂f .
f = fo + τ eE⃗ (4.93)
∂ϵ
Da bi odredili provodnost, potrebno je pomoću raspodele iz prethodnog izraza
odrediti gustinu struje, koja se uvodi kao
∑ p⃗
⃗j = −e f (⃗p) = (4.94)
m
p
⃗
∫
2
=− 3 e⃗v f d3 p⃗, (4.95)
h
zamenom izraza za neravnotežnu raspodelu f , uz aproksimaciju da se ta raspodela
malo razlikuje od ravnotežne raspodele fo , dobijamo da je
∫
⃗j = − 2 e2 ( ∂f )τ⃗v (⃗v E)d
⃗ 3 p⃗. (4.96)
h3 ∂ϵ
Tenzor provodnosti ćemo dobiti ako prethodni izraz napišemo u obliku
∑
3
⃗j = jβ ⃗eβ , (4.97)
β=1
166 § 4. NERAVNOTEŽNA STATISTIČKA FIZIKA
tako dobijamo da je
∫
2 ∂fo
σα,β = − 3 e2 ( )τ < vα vβ > D(ϵ)dϵ, (4.98)
h ∂ϵ
e2 s+5 s+5
σα,β = 8A (k B T ) 2 Γ( ). (4.99)
3kB T m∗ π 2 h3 2
§5 Ravnotežna statistička fizika
I. Mikrokanonski ansambl
On se svodi na
∫ ∞ ∫ ∞ (∫ ∞ )n ( π ) n2
−ax21 −ax2n −ax2
I(a) = e dx1 · · · e dxn = e dx = .
−∞ −∞ −∞ a
-
Rešavamo uvodenjem smene t = ar2 , te je
∫ ( ) n−1 ∫ ∞
∞
−t t 2 dt 1 ncn ( n )
t 2 −1 e−t dt = n/2 Γ
n
I(a) = ncn e √ = ncn n/2 .
0 a 2 at 2a 0 2a 2
π n/2 ncn ( n )
= Γ ,
an/2 2an/2 2
odakle se dobija da je
π n/2
cn = ( ).
Γ n2 + 1
I.. MIKROKANONSKI ANSAMBL 169
T dS = dE + pdV,
1
Koja se često naziva i de Broglie-eva talasna dužina. De Broglie je prvi izneo pretpostavku da
materija može imati i drugu prirodu, da pored čestične može imati i talasnu. Centralna tačka
njegove pretpotavke √ je zapravo njegov izraz za koji povezuje impuls čestica√p i talasnu dužinu
λB = hp . Pošto je p = 2mϵ, dobijamo da je λB = √2mϵh
, odnosno da je λB = 2πλT .
170 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
3. Koji oblik ima površina konstante energije u faznom prostoru jednog oscilatora
frekvencije ν? Odrediti fazni integral dela tog faznog prostora Γ0 (E) u faznom
prostoru sa energijom manjom od E. Zatim nadi - broj kvantnih stanja Ω0 (E) sa
E − 21 hν
n≤ .
hν
E − 21 hν
Ω0 (E) = + 1,
hν
što kada E → ∞ =⇒ Ω0 (E) ≈ E
hν
, te je zaista
Γ0 (E)
∼ Ω0 (E).
h
(c) Pokazati da Γo (E) je se održava konstantnom kada se zid kreće sporo prema
koordinatnom početku.
(d) Primenjujući formalizam kvantne mehanike odrediti Ωo (E) koji zapravo pred-
stavlja broj kvantnih stanja čija energija je manja od E. Dobijeni broj stanja
uporediti sa γo (E).
172 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
n+ ϵ0 − n− ϵ0 = M ϵ 0 =⇒ n+ − n− = M.
Sledi da je
1
n+ = (N + M ) ,
2
1
n− = (N − M ) .
2
Možemo primeniti rezon sličan prošlom zadatku. Zamislimo da su čestice na ϵ0
bele, a čestice na −ϵ0 crne. Tada je broj načina na koje ih možemo razmestiti
N! N!
WM = = ( N +M ) ( N −M ) ,
n+ !n− ! 2
! 2
!
NN
WM ≈ ( ) 2 ( N −M ) N −M
N +M .
N +M 2
2 2
1 ∂S ∂S ∂M 1 ∂S
= = · = · .
T ∂E ∂M ∂E ϵ0 ∂M
I.. MIKROKANONSKI ANSAMBL 173
1 kB 1 N − M
= · ln ,
T ϵ0 2 N + M
iz čega sledi
2ϵ0
1 − e kB T
M =N 2ϵ0 .
1+e kB T
Pošto je E = M ϵ0 biće
2ϵ0 ϵ0 ϵ0
−k
1 − e kB T e BT − e kB T ϵ0
E = ϵ0 N 2ϵ0 = ϵ0 N ϵ0 2ϵ = −ϵ0 N tanh .
− k 0T kB T
1+e kB T
e kB T
+e B
∂E N ϵ20
C= = .
∂T kB T 2 cosh2 ϵ0
kB T
Rešenje:
∑ ( p2i )
(a) Hamiltonijan ovakvog sistema je H = i 2m
+ 12 m(2πν)2 qi2 , a broj mogućih
stanja
∫∫ ∏ (∫ )N (∫ )N
1 1
Ω0 (E) = N dqi dpi = dq dp .
h N! H≤E i
N !hN
tada sledi da je
√
2 1 √
dq = dx dp = 2mdy,
m 2πν
pa
(∫ ∫ )N
1 1
Ω0 (E) = · dx dy .
N !hN (πν)N x2 +y 2 ≤E
-
Nakon sredivanja imamo da je broj mogućih stanja
( )N
1 E
Ω0 (E) = .
N !hN ν
1 ∂S
Kalorička jednačina se nalazi preko T
= ∂E
, pa je
E = N kB T.
1 3 1
ϵ = hν, hν, ..., (n + )hν. (5.2)
2 2 2
1
E = N hν + M hν. (5.3)
2
što iznosi
1
E = N hν + M hν,
2
∑N
gde je M = i=1 ni .
Potrebno je naći termodinamičku težinu WM . Modelirajmo oscilatore na sledeći
način. Imamo N kutija tako da i-ta kutija odgovara i-tom oscilatoru. U kutije
- odgovarajući oscilator.
stavljamo bele kuglice i to onoliko njih koliko je pobuden
Oscilatorima koji su u osnovnom stanju odgovaraju prazne kutije. Ukupno ima
N −1 pregradnih zidova i njih modeliramo crnim kuglicama. Dakle, ukupno M +
N − 1 kuglica. Na koliko načina možemo namestiti ove kuglice? Ukoliko je svaka
kuglica različita tada bi to bila permutacija, no kako su bele/crne nerazlučive
- sobom, broj permutacija delimo sa njihovim permutacijama
medu
(M + N − 1)!
WM = .
M !(N − 1)!
(M + N − 1)M +N −1
WM ≈ ,
M M (N − 1)N −1
te je entropija S = kB ln WM
1 ∂S
Kako je T dS = dE + pdV pa T
= ∂E
sledi
( )
1 ∂S ∂S ∂M kB M +N
= = · ≈ ln .
T ∂E ∂M ∂E hν M
hν 1 E 1 1
e kB T = 1 + =⇒ = + hν .
E
hνN
− 1
2
N hν 2 e kB T − 1
176 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
što daje
E kB T
≈ =⇒ E ≈ N kB T.
hνN hν
Rešenje:
(1 )
(a) Ako za jedan oscilator važi En = hν 2
+ n , tada je za sistem E = 21 N hν +
M hν. Pokazano je da je termodinamička težina ovakvog sistema
(M + N − 1)!
WM = .
M !(N − 1)!
jamo
M (M − 1) · · · (M − n + 1)(N − 1)
W = ,
(M + N − 1) · · · (M + N − n − 1)
što ukoliko pretpostavimo da N >> 1 i M >> n prelazi u
M nN
W ≈ .
(M + N )n+1
I.. MIKROKANONSKI ANSAMBL 177
− k ϵT
W ≈ e− 2E ϵ = e
3N
B .
9. Naći broj kvantnih stanja za gas čestica koje su smešteni u kocku stranice l i
-
uporediti dobijeni rezultat sa zapreminom u klasičnom faznom prostoru. Takode,
odrediti gustinu stanja.
Rešenje:
Šredingerova jednačina za ovu česticu je
( 2 )
~2 ∂ ∂2 ∂2
− + + ψ(x, y, z) = Eψ(x, y, z).
2m ∂x2 ∂y 2 ∂z 2
Rešenje pretpostavljamo u obliku ψ(x, y, z) = X(x)Y (y)Z(z), pa je
~2 1 d2 X ~2 1 d2 Y ~2 1 d2 Z
− · · − · · − · · = E,
2m X(x) dx2 2m Y (y) dy 2 2m Z(z) dz 2
što može da se napiše kao
~2 d2 X
− · = E (1) X(x),
2m dx2
178 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
~2 d2 Y
− · = E (2) Y (y),
2m dy 2
~2 d2 Z
− · = E (3) Z(z),
2m dz 2
pri čemu važi E = E (1) + E (2) + E (3) . Rešenje za
d2 X 2mE (1)
+ X(x) = 0
dx2 ~2
je
X(x) = A cos k1 x + B sin k1 x,
√
2mE (1)
gde je k1 = ~2
. Kako talasna funkcija u tačkama x = 0 i x = l mora biti
jednaka nuli, sledi
pa je k1 l = n1 π, te
~2 2 2 h2 2
E (1) = n π = n.
2ml2 1 8ml2 1
Analogno za E (2) i E (3) . Ukupna energija čestice je
primetite da
Γ0 (E)
Ω0 (E) ∼ .
h3
Gustina stanja je
dΩ0 (E) V 3√
Ω(E) = = 2π 3 (2mE) 2 E.
dE h
n+ + n− = n
(n+ − n− )a = x.
Odatle sledi da je
1( x)
n+ = n+ ,
2 a
1( x)
n− = n− .
2 a
Termodinamička težina će biti
n! n!
W = = (1 ( )) ( 1 ( )) ,
n+ ! n− ! 2
n + x
a
! 2
n − x
a
!
180 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
nn n!
W ≈ n+ n− = ( ( )) 1 x ( ( )) 12 (n− xa ) .
2(
n+ n− n+ a ) 1
1
n+
2
x
a
n−
2
x
a
Entropija sistema je
( ( ) ( ( )
kB ( x) 1 x) kB ( x) 1 x)
S = kB ln W = kB n ln n− n+ ln n+ − n− ln n− .
2 a 2 a 2 a 2 a
F = U − T S,
kB T n + xa
X= ln .
2a n − xa
kB T 2x kB T
X≈− · = − 2 · x.
2a na na
11. Razmatrati idealni gas N jednoatomskih molekula masa m koji se nalazi u sudu
zapremine V . Dinamika gasa je opisana hamiltonijanom
∑
3N
p2i
H= , (5.4)
i
2m
Rešenje:
(a) Statistička težina koja odgovara sloju stanja oko hipersfere sa konstantnom
energijom je jednaka
∫ ∫ ∏
3N ∫ ∫ ∏
3N ∫ ∫ ∏
3N ∫ ∫ ∏
3N
∆Λ(E) = ··· dqi ··· dpi − ··· dqi ··· dpi .
| {z } i | {z } i | {z } i | {z } i
3N 3N 3N 3N
H<E+∆ H<E
(5.5)
Nakon integracije po energijama dobijamo da je zapremina faznog prostora
u tom sloju faznog prostora
3N [ ]− 3N
V Nπ 2 3N ∆ 2
∆Λ(E) = N (2m(E + ∆)) 2 1 − (1 + ) . (5.6)
Γ( 2 + 1) E
Kako su V
N
i E
N
intenzivne veličine, samim tim izraz u [· · · ] ne zavisi od broja
čestica to je onda odavde očigledno da entropija S predstavlja ekstenzivnu
veličinu pošto ostaje zavisnost po broju čestica.
13. Strune ”kvantne violine” mogu da vibriraju na frekvencama ω, 2ω, 3ω i tako dalje.
Svaki mod vibracije može se tretirati kao nezavisni harmonijski oscilator. Razma-
traćemo mogućnost da sistem može imati nultu energiju, a da struna koja vibrira
frekvencom pω ima energiju n~pω, gde je n ∈ Z.
I.. MIKROKANONSKI ANSAMBL 183
π 2 kB
2 2
T
F =− (5.13)
6 ~ω
a. Koliko je broj mikrostanja, koji imaju najmanju vrednost energije? Koja je na-
jveća energija mikrostanja?
Rešenje:
(b) Ako je energija sistema E, onda je broj pobudjenih stanja E/ϵ. Tada je broj
realizacija takvog makrostanja jednak
N
Ω= · 2N −E/ϵ . (5.17)
E/ϵ
∑
Entropija sistema je data sa S = −kB r pr ln pr == kB ln Ω. Ako zamenimo
predhodni izraz i iskoristimo najgrublju formu Stirlingove aproksimacije
I.. MIKROKANONSKI ANSAMBL 185
N ! ∼ N N dobijamo
S ≃ kB N ln N − kB (N − E/ϵ) ln (N − E/ϵ)
E
−kB ln E/ϵ + kB (N − E/ϵ) ln 2
ϵ
(c) Nadjimo temperaturu sistema
( )
1 ∂S
=
T ∂E
( )
kB N ϵ/E − 1
= ln .
ϵ 2
Primetimo da je temperatura negativna za E > 31 N ϵ. Odatle možemo izraziti
energiju E, čijim diferenciranjem po temperaturi dobijamo toplotni kapacitet
( )
∂E 2N ϵ2 eϵ/kB T
C= =
∂T kB T 2 (2eϵ/kB T + 1)2
Isti rezultat bi dobili i u kanonskom ansamblu. Grafik izgleda tako što se za
T ≪ ϵ/kB C odozgo lepi za nulu. Za T ≫ ϵ/kB C se opet odozgo lepi za
nulu ∼ 1
T2
. Izmedju je ispupčen.
-
16. Jedan složen sistem se sastoji od podsistema A, B, C, ..., koji medusobno slabo in-
teraguju. Uzeti da sistem, zajedno sa podsistemima nalazi se u toplotnoj ravnoteži
sa okolinom. Pokazati da particiona funkcija sistema i slobodna energija sistema
zadovoljavaju sledeće relacije:
17. Izračunati toplotni kapacitet, termičku i kaloričku jednačinu stanja idealnog klasičnog
∑ p2i
gasa (H = N i=1 2m ) na visokim temperaturama.
te je
( )N
(4πV )N 2N 1 8πV
Z(V, T, N ) = · = .
N !h3N (βc)3N N! (βch)3
Slobodna energija je
8πV
F (V, T, N ) = −kB T ln Z = N kB T ln N − N kB T − N kB T ln .
(βch)3
Unutrašnja energija ( )
∂ F
U = −T 2
= 3N kB T.
∂T T
188 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
Dalje G = F + pV , a pritisak je
( )
∂F N kB T
p=− = ,
∂V T,N V
pa
8πV
G = N kB T ln N − N kB T ln .
(βch)3
Konačno, entalpija iznosi
E = U + pV = 4N kB T.
Rešenje:
∑
∞
− k~ωT (n+ 12 )
Zqu = e B ,
n=0
− 2k~ωT
∑
∞
− k~ωT n
Zqu = e B e B ,
n=0
− 2k~ωT
e B
Zqu = ,
− k~ωT
1−e B
1
Zqu = .
2 sinh 2k~ω
BT
-
(c) Kada imamo sistem koji se sastoji od N medusobno nezavisnih oscilatora,
particiona funkcija je
( )−N
N ~ω
Z= Zqu = 2 sinh ,
2kB T
pa unutrašnja energija može da se dobije iz F = −kB T ln Z, jer F = U − T S,
pa ( ) ( )
∂F ∂ F
U =F −T = −T 2
,
∂T
V ∂T T V
∂
U = kB T 2 ln Z,
∂T
1
a kako je β = kB T
,
∂ ∂β
(ln Z) ·
U = kB T 2 ,
∂β ∂T
∂
U =− ln Z.
∂β
-
Uvrštavanjem izraza za Z i sredivanjem, imamo da je unutrašnja energija
konačno
N ~ω
U= ~ω coth .
2 2kB T
190 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
Dok je entropija
∂F ∂
S=− =− (−kB T ln Z) ,
∂T ∂T
-
što nakon sredivanja iznosi
( ( ))
~ω ~ω ~ω
S = N kB coth − ln 2 sinh .
2kB T 2kB T 2kB T
H = Ho + λH1 , (5.22)
21. Idealni gas se nalazi u cilindru poluprečnika R i dužine L, rotira oko svoje ose
simetrije ugaonom brzinom ω. Zanemarujući uticaj homogenog polja Zemljine
teže, odrediti slobodnu energiju sistema na temperaturi T .
23. Izračunati toplotni kapacitet, termičku i kaloričku jednačinu stanja idealnog ul-
∑
trarelativističkog gasa (H = N i=1 cpi ) na visokim temperaturama.
II.. KANONSKI ANSAMBL 191
Rešenje:
(a) Broj domenskih zidova je n << N , jer domenski zid ne može stajati na kraju.
Problem postavljanja domenskih zidova ekvivalentan je problemu postavl-
janja n identičnih kuglica u N različitih kutija, pri čemu u svakoj kutiji sme
stajati najviše jedna kuglica jer izmedu - dva spina sme stajati samo jedan
domenski zid. Ukupan broj načina na koji se kutije mogu popuniti dobi-
jamo tako što za svaku kutiju treba da odlučimo hoće li kuglica stajati u njoj
ili ne. Pošto imamo n kuglica i N različitih kutija, u pitanju su kombinacije
n-te klase od N elemenata
( )
N N!
W = = .
n n!(N − n)!
(b) Pošto je
N!
W = (E ) ( ),
ϵ
! N − Eϵ !
entropija sistema će biti
( ( ) ( ))
E E
S(E) = kB T ln W = kB ln N ! − ln ! − ln N − ! .
ϵ ϵ
192 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
Nϵ
E(T ) = ϵ .
1 + e kB T
(c) Toplotni kapacitet je
∂E N ϵ2 1
c= = · ( )2 .
∂T kB T 2 ϵ
− 2k ϵ T
e 2kB T
+e B
E C
T T
lim E(T ) = 0,
T →0
a za T → ∞
Nϵ
lim E(T ) ≈ ϵ ,
T →∞ 1+1+ kB T
Nϵ
lim E(T ) = .
T →∞ 2
II.. KANONSKI ANSAMBL 193
lim c(T ) = 0,
T →0
lim c(T ) = 0,
T →∞
-
rešavanjem za prostorne veličine i prevodenjem u sferne koordinate imamo
[∫ ∫ ∫ (
p2
) ]N
∞ p2
x + y + pz
2
1 −β
Z= VN dpx dpy dpz e 2m 2m 2m
· (2π)N ·
N !h5N −∞
( ) N
∫ π ∫∫ ∞ p2
−β 2I +
θ
p2ϕ
−µE cos θ
· ,
2I sin2 θ
dθ dpθ dpϕ e
0 −∞
194 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
što daje
(∫ ∞ )3N [∫ π ∫ ∞ ∫ ∞ p2
]N
(2πV )N p2
−β 2m
p2
−β 2Iθ −β ϕ
Z= · dpe · dθ dpθ e dpϕ e 2I sin2 θ e βµE cos θ
.
N !h5N −∞ 0 −∞ −∞
-
Sredivanjem prethodnog izraza dobijamo
( )N (∫ π )N
(2πV )N 3N 2πI
Z= · (2πmkB T ) 2 · · dθ sin θe βµE cos θ
,
N !h5N β 0
D = εE = εo E + P,
P
odakle se dobija da je εr = 1 + εo E
. Da bi odredili relativinu dielektričnu pro-
pustljivost, potrebno je potražiti izraz (??) u limesu malih polja. Koristeći da je
x3
coth x − 1
x
≈ x
3
− 45
+ · · · , nalazimo
N µ2
ϵr = 1 + . (5.26)
V 3ϵo kB T
II.. KANONSKI ANSAMBL 195
(b) entropiju;
Rešenje:
Ukupna particiona funckija sistema ZN jednaka je N -tom stepenu particione funkcije
koja odgovara jednoj čestici Z1
µH
Z1 = eβµH + e−βµH = 2 cosh ,
kB T
µH
ZN = 2N coshN .
kB T
Sve potrebne veličine lako nalazimo ako odredimo Helmholtz-ovu slobodnu en-
ergiju F = −kB T ln ZN ,
( )
µH
F = −N kB T ln 2 cosh .
kB T
(b) Entropija,
∂F
S=− ,
∂T
( )
N µH µH µH
S= tanh − N kB ln 2 cosh .
T kB T kB T
(c) Toplotni kapacitet, ( )
∂U
c= ,
∂T H
( )2
µH 1
c = N kB · .
kB T cosh2 kµH
BT
∂F
dF = −SdT − M dH ⇒ M =− ,
∂H
µH
M = N µ tanh .
kB T
-
27. Dat je gas koji se sastoji od medusobno neinteragujućih čestica kod kojih je veza
izmedu- energije čestice ε i intenziteta njenog impulsa p oblika ε = ε0 +αp3 , gde su
ε0 i α konstante. Naći statističku sumu ovog gasa i na osnovu toga izračunati ko-
liki rad vrši jedan mol ovakvog gasa pri adijabatskoj ekspanziji iz stanja u kome
je zapremina sistema V1 , temperatura T1 u stanje u kome je zapremina V2 . Kolika
je temperatura gasa nakon ekspanzije?
28. Idealan gas N jednoatomskih molekula (pojedninačne mase m). Gas se nalazi
u beskona čno visokom cilindru (površine bazisa A) pod uticajem homogenog
gravitacionog polja. Pretpostavljajući da se gas nalazi u stanju termodinamičke
ravnoteže, odrediti specifičnu toplotu gasa.
1 N
Z(N, T, V ) = V Z1 (T )N . (5.27)
N!
II.. KANONSKI ANSAMBL 197
Cv = nN kB . (5.28)
Rešenje:
Particiona funkcija N -čestičnog sistema je
1 N N 1 N N nN
Z= V Z1 = V C T . (5.29)
N! N!
-
Nakon toga se odreduje toplotni kapacitet
( )
∂S
CV = T =
∂T V
( )
kB T ∂Z
= = nN kB . (5.31)
Z ∂T V
-
30. Jednodimenzioni polimer sastoji se od N medusobono povezanih monomera.
Svaki monomer je u obliku elipse, čija je veća poluosa dužine 2a, dok je kraća
poluosa dužine a. Monomeri mogu da se postave tako da pravac polimera može
da se poklopi sa dužom poluosom (longitudinalno postavljni), ili sa kraćom (transverzalno
postavljeni). Energijski doprinos ϵ se uvršćuje u ukupnu energiju ukoliko je
monomer trasverzalno postavljen, i zbog toga ukupna energija polimera je jed-
naka E = nϵ, gde je n broj transverrzalnih monomera.
(b) Odrediti kolika je verovatnoća da se monomeri postave duž kraće ili duže
poluose.
Rešenje:
∑
N
N!
Z(T, N ) = e−βnϵ (5.32)
n=0
n!(N − n)!
= (1 + e−βε )N . (5.33)
e−βε
pT = , (5.34)
1 + e−βε
- u longitudinalnom položaju
je verovatnoća da se monomer nade
1
pT = . (5.35)
1 + e−βε
31. Klasični gas molekula (koji se mogu aproksimovati tvrdim sferama zapremine
v) se nalazi u sudu zapremine V . Odrediti entropiju u funkciji energije sistema.
- odrediti termičku i kaloričku jednačinu stanja. Pokazati da je izotermska
Takode,
kompresibilnost uvek pozitivna.
200 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
a) Maxwell-ova raspodela
Slika 5.3:
( ) 23
m − 2km T (vx2 +vy2 +vz2 )
f (⃗v ) = e B .
2πkB T
Važi da je δV = Svx δt, pa
∑∑
−δN = n(⃗v )Svx δt,
vx ,vy vx >0
gde smo integralili po x–osi samo u pozitivnom smeru, pošto jedino ti molekuli
- kroz taj otvor (videti Sliku 5.3). Rešavanjem integrala dobijamo
mogu da produ
√
δN kB T
− =S .
δt 2πm
Za idealan gas jednačina stanja je p = nkB T , zbog čega prethodna relacija oblika
√
δN np
− =S .
δt 2πm
34. Za Maksvelov razredjeni idealni gas odrediti srednje vreme kretanja čestica u
gasu, uzimajući da su čestice oblika sfere poluprečnika a.
-
35. Razredene gas koji se nalazi u zapremini V na pritisku p. Pretpostavljajući da
je rapodela molekula gasa Maxwell-ova, izračunajte meru (brzinu) isticanja gasa
kroz otvor male površine S.
(a) koji imaju jednu komponentu veću od 2vm , bez obzira na vrednosti druge
dve komponente.
37. Koliko procenata od ukupnog broja udara molekula o jedinicu površine zida
-
suda u jedinicu vremena otpada na molekule koji su odredeni uslovima koji su
dati u prethodnom zadatku pod (a) i pod (b)?
- energija ϵ i
38. Odrediti broj molekula po jedinici zapremine sa energijama izmedu
ϵ + dϵ u gasu koji se nalazi u stanju termodinamičke ravnoteže na temperaturi T .
Na osnovu dobijenog rezultata izračunati (nerelativistički):
39. 8 Klasičan idealni gas se nalazi u sudu na temperaturi T . Čestice gasa (svaka
mase m) mogu da izlaze kroz mali otvor na zidu suda. Smatrajući da isticanje
gasa ne remeti termodinamičku ravnotežu gasa u sudu, odrediti oblik funkcije
raspodele po intenzitetima brzina čestica gasa van suda. Izračunati zatim nji-
hovu brzinu van suda. izračunati zatim najverovatniju brzinu i srednju vrednost
kinetičke energije.
Rešenje:
- iz kutije u nekom intervalu δt su medusobno
Broj molekula gasa koji izadu - povezani
relacijama
∑ ∑
δN = − n(⃗v )δV = − nvz f (⃗v )dSdt. (5.47)
⃗v ⃗v
- izvan suda je jednak ovom broju,
Kao što se očekuje, broj molekula koji izade
2
− 2k
mv
i raspodela je samim proporcionalna Cvz e T
B . Konstanta C se može odrediti
normiranjem
∫ ∞ ∫ π ∫ π 2
− 2k
mv
Cv cos θe T
B = 1. (5.48)
0 0 0
II.. KANONSKI ANSAMBL 203
( )2
2
Tako dobijamo da je C = 2 2kB T
, što znači d je raspodela po brzinama oblika
( )2
2 − 2k
mv 2
mvi
pi = √
1 − vc2
2
2
Koristeći Maxwell-ovu raspodelu odrediti < √ mv >.
2
1−( vc )
-
kojom se odreduje verovatnoća dp(v1 , v2 ) da se brzina jedne čestice nalazi u inter-
valu (v1 , v1 + dv1 ) i brzina druge u intervalu (v2 , v2 + dv2 ), i koja se može prikazati
izrazom
dp(v1 , v2 ) = f (v1 )f (v2 )dv1 dv2 . (5.53)
Tako dobijamo da je
( ) 12 √
1 m mv 2
− 2k m mv 2
F (v) = e B T
[2 v( − 3) +
2π 2kB T 2kB T 2kB T
( )2 (( ) )
mv 2 √ mv 2
mv 2
mv 2
1
e 4kB T 2π(3 − + )]Erf .
kB T 2kB T 2kB T 2
Srednja brzina je √
∫ ∞
8 kB T
< v >= v · F (v)dv = 2 · , (5.55)
0 π m
dok je srednja kvadratna brzina jednaka
∫ ∞
6kB T 8 kB T
2
< v >= v 2 · F (v)dv = + . (5.56)
0 m π m
Time se dobije da je srednjs vrednost zbira brzine čestica < v >= 2 < v1 >, dok je
srednja vrednost kvadrata brzine jednak < v 2 >= 2(< v12 > + < v1 >2 ).
(a) Odrediti koliko je udara molekula po jedinici površine zida suda u jedinici
√
vremena, čiji je intenzitet brzine veći od 3vm , gde je vm = 2kmB T najverovat-
nija brzina.
(b) Odrediti koliko je udara molekula po jedinici površine zida suda u jedinici
√
vremena, čija je jedna od komponenti brzine veći od 4vm , gde je vm = 2kmB T
najverovatnija brzina.
II.. KANONSKI ANSAMBL 205
∫x
e−t dt.
2
Prikaži rezultate koristeći error funkciju Υ(x) = √2
π 0
Rešenje:
dS
dq
Slika 5.4:
(a) Samo atomi sa brzinom v, koji dolase iz koordinata odredjenih sa (θ, θ + dθ)
i (ϕ, ϕ + dϕ), koji imaju od svih mogućih brzina usmerenu tako da pro-
lazi kroz dS će se sudariti sa tom površinom. Od svih atoma u zapremini
dV = vdt · rdθ · r sin θdϕ sa intenzitetom brzine u intervalu (v, v + dv), kojih
ukupno ima dV · n P (v)dv, samo njih dS cos θ
4πr 2
ima brzinu u dobrom smeru
(videti sliku 5.4), tako da će za vreme dt proči kroz deo sfere na rastojanju r i
posle odredjenog vremena sudariti se sa površinom. U predhodnoj formuli,
P (v)dv je verovatnoća da atom ima intenzitet brzine v (bez obzira na smer
brzine!), a n je koncentracija gasa u sudu. Da bismo pronašli sve atome koji
imaju brzinu 3vm i veću, a koji se sudaraju sa površinom, treba samo proin-
tegraliti predhodno dobijenu jednačinu po svim brzinama većim od 3vm :
∫ ∞
dS cos θ
dN = r dθ r sin θ dϕ · vdt · n P (v)dv · . (5.57)
3vm 4πr2
206 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
(b) Postavimo isti integral kao i predhodnom slučaju, samo sada, ako želimo da
je komponenta brzine duž z ose veća od nekog zadatog v ′ , od svih atoma koji
dolaze pod uglom θ na površinu dS treba uračunati samo one sa brzinama
v ′ cos θ i većim:
∫ π/2 ∫ ∞
Cn
v 3 e−mv
2 /2k
bT
dN = dtdS dθ sin θ cos θ dv
2 0 4vm sin θ
∫
Cn π/2
= dtdS dθ sin θ cos θ ·
2 0
( )2
2kb T 1 v2 cos θ 2
· e (v cos θ + 1),
m 2
mv ′2
gde smo u poslednjem redu uveli oznaki v ≡ kb T
= 42 = 16. Ako uvedemo
oznaku x ≡ v 2 cos2 θ dobijamo rešivi integral tipa
∫ v
1 2 1
e−x (x2 + 1)dx2 = e−v (v 2 + 1) + (1 − e−v ).
2 2
(5.60)
0 2 2
Zamenom u predhodnu jednačinu, nakon zamene konstante C, zamene v =
4 i sredjivanja, dobijamo:
√ ( )
dN 2kB T 1 1
=n − e−16 (5.61)
dSdt πm 32 2
II.. KANONSKI ANSAMBL 207
√
43. Pokazati da je relativna brzina dva molekula u Maxwell-ovoj raspodeli 2 puta
veća nego brzina tih molekula u odnosu na zid suda.
-
45. Ravna daščica, površine S se kreće brzinom ⃗u u jako razredenom gasu kon-
centracije n na temperaturi T . Masa molekula gasa je m. Uzeti da je brzina
daščice mnogo manja od srednje brzine molekula gasa v, i da je brzina daščice
u proizvoljnom pravcu u odnosu na površinu daščice. Odrediti kolika je vred-
nost F⃗ koja deluje na tu daščicu.
⃗ U
46. Sistem koji se sastoji od N čestica, spina 12 , nalazi se u magnetnom polju H.
slučaju da čestice poseduju magnetni moment µ, koji može biti orijentisan samo u
dva smera ↑ ili ↓ (koji je u istom ili suprotnom smeru od polje koje je npr. ↑), onda
su odgovarajući energijski nivoi svake čestice +µH, ili −µH. Odrediti particionu
funkciju, slobodnu energiju, entropiju i srednju vrednost magnetnog momenta i
njegovu disperziju.
47. Atomi mogu posedovati magnetne momente, koji mogu uzeti vrednosti iz diskretnog
skupa vrednosti skupa gµB m (magnetni moment uzima brojeve J, J − 1, J −
2, ..., −J) kao njihove komponente duž magnetnog polja H. Izračunajte mag-
netizaciju M tela koje sadrži n takvih atoma po jedinici zapremine. Izračunati
susceptibilnost χ za slabo polje interakcije na visokim temperaturama (kB T >>
gµB JH).
208 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
(d) Nadji specifičnu toplotu uzorka kao funkciju temperature na konstantnoj za-
premini. Diskutovati specifičnu toplotu na niskim i visokim temperaturskim
granicama. Skicirati Cv (T ).
Rešenje:
ϵ ϵ
− Kϵo T − K 1T − K 2T
F − KB T N ln Z1 = −N kB T ln(e B + 3e B + 5e B ) (5.63)
tako dobijamo da je za J = 2:
5N e−βϵ2
n2 = − Kϵo T − K 1T
ϵ ϵ
− K 2T
, (5.65)
e B + 3e B + 5e B
n1 3
= ϵ1 −ϵo ϵ1 −ϵ2 , (5.66)
N e KB T
+ 3 + 5e KB T
i za J = 0
no 1
= ϵo −ϵ1 ϵo −ϵ2 . (5.67)
N 1 + 3e KB T + 5e KB T
(d) Unutrašnja energija se dobija pomoću sledećih identiteta
∂ ∂
UN = − ln ZN = −N ln Z1 . (5.68)
∂β ∂β
Istovremeno, važi:
∑ ∂Z
(n + 1/2)~ω e−(n+1/2)~ω/kB T = −kB T ω
n
∂ω
~ω eβ + e−β
= , (5.73)
2 (eβ − e−β )2
~ω
gde je β ≡ 2kB T
. Zamenom (5.73) i (5.72) u (5.71) dobijamo:
( )
~ω ~ω
E= coth
2 2kB T
c) T − p raspodela
53. Polimerni lanac se sastoji od N monomera, koji svi imaju dužinu a. Monomeri
-
su medusobno povezani svojim krajevima, tako da slobodno mogu rotirati oko
-
ovih krajeva. Odrediti zavisnost sile zatezanja polimera f od rastojanja z izmedu
krajeva polimera.
Rešenje:
Monomeri, konstannte dužine a, su postavljeni u proizvoljnom pravcu u trodi-
menzionom prostoru. z-osu postavimo tao da prolazi kroz dva krajnja monomera
tog polimera. Ugao θi koji zaklapa prava na kojoj se nalazi taj monomer sa z-
osom. Unutrašnja energija U se menja u skladu sa zakonom
dU = T dS + f dl, (5.74)
gde je prostorni ugao dΩi = sin θi dθi dφi . Takod dobijamo, posle integracije po θi i
φi , dobija se da je particiona funkcija
[ ]
4λ sinh(βf a)
G(T, f, N ) = −kb T N ln , (5.77)
βf a
odakle se dobija da je
( )
∂G kb T N
z=− = (coth(βf a)βf a − 1). (5.78)
∂f T f
212 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
a a b b a
a a b b a
Slika 5.5:
54. 8 N monomernih jedinica se nalaze duž prave linije, i formiraju lanac molekula.
Svaki molekul može da bude u α i β stanju. Ako se molekuli nalaze u stanju α,
onda je energija tih molekula Eα i njihova dužina a (vidi sliku (5.5)), dok en-
-
ergija molekula koji se nalazi u stanju β je Eβ i dužina b. Odrediti relaciju izmedu
dužine L i sile zatezanja X koja deluje na jednom kraju lanca molekula.
Rešenje:
Neka je Nα monomera u stanju α, dok je Nβ monomera u stanju β. Ukupan broj
monomera je N, tako da relacija N = Nα +Nβ mora biti ispunjen. Takodje i dužina
lanca je povezana sa brojem monomera relacijom L = Nα a + Nβ b, a energija lanca
relacijom E = Nα Eα + Nβ Eβ . Pošto nas interesuje veza izmedju dužine lanca L i
sile X, koristićemo T − X raspodelu. Gustina verovatnoće pomenute raspodele
je
e−β(E−XL)
P r(Nα , Nβ ) = , (5.80)
Y
gde je Y particiona funkcija raspodele. Particiconu funkciju Y odredićemo iz
II.. KANONSKI ANSAMBL 213
uslova
∑
P r(Nα , Nβ ) = 1, (5.81)
po svim konfiguracijama
odakle dobijamo da je
∑
N
N!
Y = e−β[E(Nα ,Nβ )−XL(Nα ,Nβ )] . (5.82)
Nα =0
Nα !Nβ !
Y = [e kB T
+e kB T
]N , (5.83)
G = −kB T ln Y = −N kB T ln(e kB T
+e kB T
). (5.84)
55. Odrediti srednju vrednost kvadrata i kuba fluktuacije < ∆U 2 > i < ∆U 3 > ide-
alnog klasičnog gasa N dvoatomskih molekula, mase m, momenta inercije I i
dielektričnog momenta µ0 . Gas se nalazi u sudu zapremine V i jačine električnog
momenta E.
Rešenje:
A B C
4
1 2 3 4 1 2 3 1 2 3
4
D E F
4
3 4 3 4 3
1 2 1 2 1 2
G H I
1 2 1 2 1 2
3 4 3 4 3
4
Slika 5.6:
Na slici (5.6) su prikazane moguće konfiguracije linearnog polimera koji ima os-
obinu da ne može da se samo preseca. Broj konfiguracija je prikazan u Tabeli I.
II.. KANONSKI ANSAMBL 217
9 + 20e−βε + 12e−2βϵ
< r2 >= a2 . (5.87)
4 + 16e−βε + 16e−2βε
A 4 0 3a
√
B 4 ε a 5
√
C 4 ε a 5
√
D 4 2ε a 5
√
E 4 ε a 5
F 4 2ε a
√
G 4 2ε a 5
√
H 4 ε a 5
I 4 2ε a
gde je < E1 > srednja kinetička energija jednog molekula tog gasa < E1 >=<
p2
2m
>. Pored toga koristeći vezu srednje energije i pritiska P , dobija se da je
pritisak
1 p2 p2 Γ[7/4]
P = N< >= N o . (5.92)
V 2m 2m Γ[5/4]
gde se posmatrana veza može naći u n-tom stanju (n ∈ {0, 1, 2, ...}). Polimer
se nalazi u rastvoru, temperature T i pod dejstvom sile zatezanja F . Odred-
iti srednju energiju, entropiju i dužinu polimernog lanca. Odrediti asimptotiku
II.. KANONSKI ANSAMBL 219
YN = Y1N , (5.94)
-
a zatim posle odredivanja prvog izvoda dobijamo
( )
β~ωa β~ωb
eβF a ~ωa cosh( 2 ) eβF b ~ωb cosh ( 2 )
( 2 sinhaFβ~ωa − 4 sinh2 ( β~ωa )
) + ( bF
(
β~ωb
) − 4 sinh2 β~ωb
)
2 ( 2 ) 2
2 2 sinh
< E >= −N
2 2
eβF a eβF b .
+ ( )
2 sinh( β~ω
2 )
a
2 sinh
β~ωb
2
(5.103)
Entropija se dobija izrazom
( )
∂G
S = − = (5.104)
∂T F
( )
eβF a eβF b
= kB N ln ( a) + ( β~ωb ) + (5.105)
2 sinh β~ω 2
2 sinh 2
( )
β~ωa β~ωb
~ωa cosh( 2 ) ~ωb cosh ( 2 )
eβF a ( 2 sinhaFβ~ωa − 4 sinh ( β~ωa )
) + eβF b ( bF
(
β~ωb
) − 4 sinh2 β~ωb
)
N ( 2 ) 2
2 2 sinh 2 2
+ eβF a e(βF b )
(5.106)
,
T +
sinh( β~ω
2 )
a
sinh
β~ωb
2
√ ( ) 12
2 p2
E= m2 c4 + p2 c2 = mc q+ 2 2 .
mc
Razvojem u red imamo da je približno
( ( 4 ))
p2 1 p4 p
E ≈ mc q +
2
2 2
− · 4 4 +o ,
2m c 8 mc m4 c4
II.. KANONSKI ANSAMBL 221
odnosno
p2 p4
E ≈ mc2 + − .
2m 8m3 c2
Pošto zanemarujemo unutrašnju energiju, onda od ukupne energije možemo oduzeti
- gas je idealan, pa je hamiltonijan sistema N
energiju mirovanja mc2 . Takode,
čestica ( → →
)
∑
N
| pi |2 | pi |4
H= − .
i=1
2m 8m3 c2
Particiona funkcija je
N
∫ (∫ ∫ ∫ )N ∫ ∫ ∫ ∞ p4
1 p2
− 2mk 3
Z(V, T, N ) = e−βH dΓ = · dx dy dz · e BT
|e
8m3 c2 kB T
{z } dp .
N !h3N V −∞
ovo je popravka
p4
Kako je 8m3 c2 kB T
malo, možemo razviti u red
( ∫ ∞ ( ))N
VN 2
2
p
− 2mk p4
Z(V, T, N ) = 4π p dp e BT 1+ ,
N !h3N 0 8m3 c2 kB T
(∫ ∞ ∫ ∞ )N
(4πV )N p
2 − 2mkB T
2
1 p 2
6 − 2mkB T
Z(V, T, N ) = p e dp + p e .
N !h3N 0 8m c2 kB T
3
0
te je
( ( ) ( ))N
(4πV )N 1 3 3 1 7 7
Z(V, T, N ) = (2mkB T ) 2 Γ + (2mk B T ) 2Γ ,
N !h3N 2 2 8m3 c2 kB T 2
( ) 3N ( )N
VN 2πmkB T 2 15 kB T
Z(V, T, N ) = 1+ · .
N! h2 8 mc2
Tada je
( ( )3 ( ))
2πmkB T 2 N 15 kB T
F (V, T, N ) = −N kB ln Z = −N kB T ln − ln + ln V + ln 1 + · ,
h2 e 8 mc2
222 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
61. Odrediti particionu funkciju, slobodnu energiju i termičku jednačinu stanja, rela-
tivističkog gasa čija je disperziona relacija
√
ϵ= p 2 c2 + m2 c4 . (5.109)
U
Pokazati da je pritisak p u zavisnosti od unutrašnje energije po čestici u = N
prikazano izrazima:
2
a. nisko-temperaturskom slučaju kB T << mc2 : p = − mc2 )(1 − 65 u−mc
2
3v
(u mc2
+ ...)
( 2 )2
u
b. visoko-temperaturskom slučaju kB T >> mc2 : p = 3v (1 − 23 mcu + ...)
62. Sistem se sastoji od N neintereagujućih atoma koji se mogu naći u dva stanja, čije
su energije E1 = 0 i E2 = ϵ.
(c) Izračunati toplotni kapacitet sistema kada je broj ukupan atoma fiksiran.
Rešenje:
Radimo u kanonskom ansamblu (možemo raditi i u mikrokanonskom ansam-
blu, ali u termodinamičkom limitu moramo dobiti iste rezultate). Pošto se sistem
II.. KANONSKI ANSAMBL 223
F = −kB T ln Z
( )
= −kB T N ln 1 + e−ϵ/kB T .
( ∂F )
Entropiju nalazimo kao S = − ∂T v,N
( ) e−ϵ/kB T ϵN
S = kB N ln 1 + e−ϵ/kB T + −ϵ/k T
· . (5.110)
1+e B T
Energiju sistema nalazimo kao srednju energiju sistema:
1 ∑
E = E e−E{states}/kB T
Z
{states}
( )
2 ∂ ln Z
= kB T
∂T
E
⇒ e−ϵ/kB T = . (5.111)
Nϵ − E
Zamenom (5.111) u (5.110) dobijamo traženi rezultat:
( )
Nϵ E
S = kB N ln + (5.112)
Nϵ − E T
Temperaturu nalazimo direktno iz (5.111)
ϵ 1
T =− ( E ),
kB ln N ϵ−E
5 kB T
a. nisko-temperaturskom slučaju kB T << mc2 : C = kB N (1 + 2 mc2
+ ...)
( 2
)2
b. visoko-temperaturskom slučaju kB T >> mc2 : C = 3kB N (1 − 1
6
mc
kB T
+ ...)
-
65. Razmotrimo sistem N medusobno neinteragujućih dimera, čiji je spin s = 12 . Sis-
tem dimera je opisan hamiltonijanom
∑
Ĥ = ⃗i,1 · S
(S ⃗i,2 ), (5.116)
i
gde je i broj dimera, a broj 1 i 2 predstavljaju brojeve spina koji pripadaju i-tom
dimeru.
-
66. Razmotrimo sistem N medusobno neinteragujućih dimera, čiji je spin s = 12 . Sis-
tem dimera je opisan hamiltonijanom
∑
Ĥ = ⃗i,1 · S
(S ⃗i,2 ), (5.118)
i
gde je i broj dimera, a broj 1 i 2 predstavljaju brojeve spina koji pripadaju i-tom
dimeru.
226 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
-
67. Hamiltonijan dugog lanca koji je sastavljen od n medusobno različitih atoma je
dat jednačinom
∑n
p2i ∑
n−1
mω 2
H= + (qi+1 − qi − a)2 , (5.120)
i=1
2m i=1
2
∑
n−1
′
V = −c (qi+1 − qi − a)3 , c > 0, (5.122)
i=1
c. Izračunati srednju dužinu lanca l =< qn − q >, i na osnovu nje odrediti ko-
1 ∂l
eficijent termmalnog širenja α = u prvom redu razvoja po parmetru c.
l ∂T
-
UPUTSTVO: Koristiti prilaz kako je uradeno pod b.
-
68. Razmotrimo atome rasporedene -
po pravoj tako da je njihov položaj odreden
- u jednom od tri stanja, obeleženih in-
parametrom x. Čestice mogu da se nadu
-
deksom n = 1, 2, 3. Energije čestice su odredene izrazima
1
ϵ1 = ϵ3 = ax2 − bx
2
2
ϵ2 = ax + bx,
F (x,T )
a. Odrediti Helmholtz-ovu slobodnu energiju f (x, T ) = N
po čestici uzi-
-
majući u obzir da ne postoji medusobna interakcija.
Za idealan gas je
∫ ∫ ∑
1 p⃗i 2
ZN (T, V ) = dx1 dy1 dz1 · · · dxN dyN dzN dpx1 dpy1 dpz1 · · · dpxN dpyN dpzN e−β i 2m ,
N !h3N
(∫ ∞ )3N
VN p
−β 2m
2
ZN = dp e ,
N !h3N −∞
VN 3N
ZN (T, V ) = 3N
(2πmkB T ) 2 .
N !h
Sledi da je
∑∞ ( )N
1 V βµ 3
Ξ(µ, T, V ) = · 3
e (2πmkB T ) 2 ,
N =0
N ! h
( )3
2πmkB T 2
eβµ V
Ξ(µ, T, V ) = e h2 = eλf .
Rešenje:
Kada n atoma postavimo na površinu imaćemo n rupa u kristalu i ukupno N + n
čvorova rešetke. Entropija ovog sistema je
(N + n)
S(n) = kB ln ,
N !n!
što je približno S(n) = kB [(N + n) ln(N + n) − n ln n − N ln N ]. Slobodna energija
F (n) = E(n) − T S(n) je
∂F
Stanju ravnoteže odgovara minimum slobodne energije 0 = ∂n
, pa se dobija da
je
n − ω
= e kB T .
n+N
73. U ravnotežnom sistemu velikog broja čestica postoji adsorpcioni centar u kojoj
može da bude jedna ili nijedna čestica energije ϵ i energije 0. Izračunajte un-
utrašnju enrgiju i entropiju tog podsistema kao funckiju temperature T i očekivanog
broja čestica N (< 1). Izračunati relativnu fluktuaciju broja čestica.
Rešenje:
Particiona funkcija da je Na molekula adsorbovanih na N čvorova je jednaka
∑
′ −β(ϵi1 +ϵi2 +·+ϵiN )
ZN a = e , (5.124)
i1 ,i2 ,...,iN
ili da se jedan se molekul postavi duž noramle na površinu (duž z-ose). Energija
molekula u z-smeru je ε > 0, dok je energija molekula koji su usmereni duž ili y-
ose je jednaka 0. Odrediti broj mikrostanja sistema energije E, odredi broj stanja
Ω(E, N ) i entropiju S(E, N ). Odredi takode- toplotni kapacitet C i skicirati tu
78. Čvrsta faza i gasna faza (para) su u ravnoteži u zatvorenom sudu zapremine V ,
i na temperaturi T . Pretpostavimo da se statistička suma čvrste faze može pred-
staviti u obliku Zs (T, Ns ) = zs (T )Ns , gde je Ns broj molekula (koji se uzima da
su jednoatomni zbog jednostavnosti) čvrste faze, i da gasnu fazu možemo sma-
trati idealnim gasom koji se sastoji Ng molekula. Pokaži da se uslov ravnoteže,
za Ns >> 1 i Ng >> 1, može aproksimativno prikazati u obliku izraza:
zg (T, V )
Ng = (5.128)
zs (T )
79. Uvodeći novu veličinu ω u slučaju idealnog gasa opisanog velikim kanonskim
ansamblom
∑
ω = ln exp(βµN − βEf,N ). (5.137)
N,f
∑ ( )
( ) − N
N,f ζ β
∂Ef,N
exp(−βEf,N )
∂ω ∂V
= ∑ ζ,β
. (5.138)
∂V z,β N,f ζ N exp(−βEf,N )
( )
∂Ef,N
Ovde izvod − ∂V
predstavlja pritisak sistema koji se nalazi u mikrostanju
ζ,β
f sa brojemčestica N . Srednja vrednost ove veličine predstavlja pritisak sistema,
tako da važi
( )
∂ω
= βp. (5.139)
∂V ζ,β
∑
ω = ln ζ N zN = (5.140)
N
∑ 1 ( V ζ )N
= ln = (5.141)
N
N ! λ3T
Vζ
= . (5.142)
λ3T
ζ
βp = . (5.144)
λ3T
pV
ω= . (5.145)
kB T
81. (a) Pokazati da za sistem u velikom kanonskom ansamblu važi sledeća relacija
( )
∂U
< N E > − < N > < E >= < (∆N )2 > . (5.146)
∂N T,V
III.. VELIKI KANONSKI ANSAMBL 235
Rešenje:
1 ∑ βµN −βEk,N
< N >= e e , (5.147)
Ξ k,N
Tako dobijamo da je
( )
∂<N >
= β(< N 2 > − < N >2 ) = β < (∆N )2 > . (5.150)
∂µ T,V
da je
( ) (
1 ∂<E> 1 1∑
= N βEk,N eβµN −βEk,N
β ∂µ T,V β Ξ k,N
( ) ∑ )
1 ∂Ξ
− Ek,N eβµN −βEk,N
Ξ2 ∂µ T,V k,N
1
= [β < N · E > −β < N >< E >]
β
= < N · E > − < N > · < E >, (5.153)
(b) Koristićemo prethodno dobijeni izraz, tako da je relacija koju ćemo iskoristiti
( )
∂U
< N E >=< N > < E > + < (∆N )2 > . (5.154)
∂N T,V
Velika particiona funkcija Ξ se može izračunati
∑
+∞
Ξ= eβµN Z(N, T, V ), (5.155)
N =0
VN 3N
gde je pariciona funkcija za klaični idelani gas Z(N, T, V ) = N !h3N
(2πmkB T ) 2 .
1 N
Ova particiona funkcija se može zapisati kao Z(N, T, V ) = z ,
N! 1
zbog čega
velika statisitčka suma za idealana gas
3
< E >= < N > kB T, (5.159)
2
odakle zaključujemo da je
( )
∂<E> 3
= kB T. (5.160)
∂<N > T,V 2
3 V 3 V 3
< N · E >= kB T · E βµ 3 (2πmkB T ) 2 (1 + 3 eβµ (2πmkB T ) 2 ). (5.162)
2 h h
82. Odrediti srednju vrednost kvadrata i kuba fluktuacije < ∆U 2 > i < ∆U 3 > ide-
alnog klasičnog gasa N dvoatomskih molekula, mase m, momenta inercije I i
dielektričnog momenta µ0 . Gas se nalazi u sudu zapremine V i jačine električnog
momenta E.
- un-
b. Posmatrati podsistem, koji se zapravo predstavlja kružni prsten odreden
utrašnjim poluprečnikom r, i spoljašnjim r+dr, i koji obuhvata sloj vazduha na
visini (z, z + dz). Svaki od tih prstenova predstavlja veliki kanonski ansambl.
Pokazati da velika particiona funkcija zadovoljava Ξ(T, Vp , µ, r, z) = ΞO (T, Vp , µef f (r, z)),
gde ΞO velika particiona funkcija za slučaj da cilindar ne rotira, i da se jačina
gravitacionog polja se može zanemraiti (ili da je jednaka 0).
p1 p2
S = −kB [p1 ln( ) + p2 ln( )], (5.163)
g1 g2
gde je uzeto da je p1 + p2 = 1.
b. Pokazati da je entropija
g2 −x x
S = kB [ln g1 + ln(1 + e )+ ], (5.164)
g1 1 + gg12 ex
u slučaju da sistem ispoljava osobine kanonskog ansambla. U prethodnom
ϵ2 −ϵ1
izrazu parametar je korišten parametar x = kB T
, gde je T temperatura sis-
tema.
J = U − N µ = T S − pV. (5.166)
III.. VELIKI KANONSKI ANSAMBL 239
C deo izraza je već u obliku koji odgovara jednom delu rešenja (??). Ostalo je da
transformišemo A i B u oblik koji odgovara rešenju. Kako je
1 ∑ −β(E−µN )
U= Ee ,
Z
vidimo da je
( )
∂U 1 ∑ ... ∑ 1∑
=− 2 Ee (−E)e... + (−E 2 )e... ,
∂β ξ | Z {z } | Z {z }
⟨E⟩2 −⟨E 2 ⟩
( )
2 ∂U
odakle dobijamo kB T ∂T ξ
= (∆U )2 . Sada hoćemo da raspišemo ovaj parcijalni
izvod po T u obliku bližem rešenju. Kako je U = U (T, N, V ) = U (T, N (T, ξ, V ), V ),
dobijamo ( ) ( ) ( ) ( )
∂U ∂U ∂U ∂N
= + . (5.168)
∂T ξ ∂T N,V ∂N T,V ∂T ξ
Prvi član u gornjem zbiru je zapravo CV odakle dobijamo prvi član rešenja kB T 2 CV
(??). Drugi član u zbiru sadrži jedan član koji se javlja i u konačnom rešenju, ali
drugi je potrebno transformisati. Kako je
1∑
N = ⟨N ⟩ = N e−βE ξ N ,
Z
240 § 5. RAVNOTEŽNA STATISTIČKA FIZIKA
dobijamo:
( ) ∑( )
∂N 1 ∑ −βE N E
= − 2 Ne ξ · − 2
e−βE ξ N
∂T ξ Z kB T
∑ ( )
1 NE
+ − e−βE ξ N
Z kB T 2
1 1 ∑ N −βE N ∑ −βE N
= e ξ · Ee ξ
T Z2 kB T
1 1 ∑ N E −βE N
− e ξ
TZ k T
( ) B
1 ∂U
=
T ∂µ T,V
( ) ( )
1 ∂U ∂N
= . (5.169)
T ∂N T,V ∂µ T,V
( )
∂N 1 ∑ ∑ 1∑
= − 2 N e... βN e... + N · βN e...
∂µ T,V | Z {z } |Z {z }
β⟨N ⟩2 β⟨N 2 ⟩
1
= ⟨(∆N )2 ⟩.
kB T
( )
∂U
⟨(∆U ) ⟩ = kB T
2 2
∂T ξ
( )2
∂U
2
= kB T CV + ⟨(∆N )2 ⟩.
∂N T,V
u(r)
0
ε
-ε
Slika 6.1:
1. Odrediti drugi virijalni koeficijent za realni (neidalni) gas, čiji molekuli interaguju
potencijalom
∞, za r < r0 ;
u(r) = ( )
−ϵ r0 6 , za r > r ,
(6.1)
r 0
-
gde je ϵ minimalna vrednost potencijalne enrgije (vidi sliku 6.1). Koristeći odredeni
drugi virijalni koeficijent, naći prvu popravku termičke jednačine stanja u odnosu
243
i tako dolazimo do
ro3 π
b2 = −4π + 4πβεro3 + β 2 ε2 ro3 . (6.7)
3 3
Za interagujuće sistema, kada patricionu funkciju raspišemo koristeći Majerov
razvoj, dobijamo [ ]
∑
∞
ξn
Z = exp V bl 3
l=1
λ
( )
∂ ln Z
N = ξ
∂ξ
V ∑
∞
= 3
l bl ξ l
λ l=1
Odatle dobijamo:
n λ 3 = ξ + 2 b2 ξ 2 + . . .
|{z}
≡x
244 § 6. RAVNOTEŽNA STATISTIČKA FIZIKA INTERAGUJUĆIH SISTEMA
pV ∑V ∞
ln Z = = b ξl.
3 l
kB T l=1
λ
pλ3
= b1 ξ + b2 ξ 2
kB T
= b1 (x + C2 x2 + . . .) + b2 (x + C2 x2 . . .)2
Odakle dobijamo:
p ∑ ∞
=n+ Bl nl ,
kB T l=2
1 ∂V 1
kT = − = ( 3
( ) )
V ∂p r0
kB T n + 2π 3 1 + ϵ
n2
kB T
α
u(r) = , san > 0. (6.8)
rn
- gasa
3. U skladu sa Lenard-Jones-ovom jednačinom interakcije izmedu
A B
u(r) = 12
− 6, (6.13)
r r
Odrediti drugi virijalni koeficijent Lennard-Jones-ovog gasa i uporediti dobijeni
izraz sa termičkom jednačinom van der Waals-ovog gasa.
ro r1
-u o
Slika 6.2:
-
gde se drugi virijalni koeficijent odreduje izrazom
∫ ∫
1 1
B2 (T ) = − d⃗rf (r) = − d⃗r(e−βU − 1), (6.16)
2 2
odnosno da je ∫ ∞
B2 (T ) = −2π r2 dr(e−βU (r) − 1). (6.17)
0
2π 3 uo
B2′ (T ) = − 2
(ro − r13 )eβuo . (6.21)
3T kB
-
5. Čestice gasa u d = 3 dimenzija medusobno interaguju parnom interakcijom u(|⃗r2 −
⃗r1 |) = u(r), koja je data funkcijom
∞, za
0 < r < r0 ;
u(r) = −u0 , za r0 < r < r1 ; (6.22)
0, za r1 < r < ∞.
Rešenje:
1 4πro3
B2 (T ) ≈ , (6.23)
2 3
dok se za T → 0 dobija
2π ( 2 )
B2 (T ) ≈ − r1 − ro2 exp(βuo ). (6.24)
3
248 § 6. RAVNOTEŽNA STATISTIČKA FIZIKA INTERAGUJUĆIH SISTEMA
odakle se dobija
( )2
1 uo 4π 3 2πro3
(p + N2 (r1 − ro3 ))(V − N ) = N kB T. (6.30)
V | 2 3{z } 3}
| {z
a b
x
ro
-ro ro x
-ro
∫+∞
∫
1
B3 = − dxdx′ f (x)f (x′ )f (x − x′ ), (6.32)
3
−∞
∫+r
∫o
1
B3 = − dxdx′ (−1), (6.35)
3
−ro
−ro <x−x′ <ro
B3 = ro3 . (6.36)
9. Odrediti popravku najnižeg reda koji pokazuje razliku u odnosu na idealni gas za
Helmholtz-ovu slobodnu energiju F , Gibbs-ovu slobodnu energiju G, entropiju
S, unutrašnju energiju U specifičnu toplotu pri konstantnoj zapremini CV i specifičnu
toplotu pri konstantnom pritisku Cp .
11. Graf g poseduje n = nl +n2 −1 čvorova koji se mogu razdvojiti na dva podgrafa g1
i g2 sa n1 i n2 čvora, respektivno. Odrediti izraz kojim se povezuje klaster integral
a(g) sa a(g1 ) i a(g2 ).
12. (a) Izračunati drugi virijalni koeficijent za d-dimenzioni sistem, kada je poten-
- dve čestice dat relacijom
cijal izmedu
( a )j
u(r) = ε , (6.44)
r
sa a > 0 i j > 3.
∑
∞
Ξ= eβµN ZN
N =0
∑
∞ ∫ ∑
∞ ∑ −βE{nf }
Ξ= e βµN
e−βHN dΓN = eβµN e .
N =0 N =0 {nf }
∑
Kako je E{nf } = f nf ϵf , biće
∑
∞ ∑
∑
Ξ= e−β f (ϵ−µ)nf ,
N =0 {nf }
∑
∞ ∏
Ξ= e−β(ϵf −µ)nf .
∑
nf =0 f
∑
∞ ∏ n
∑
∞
n
af f = an1 1 an2 2 · · · af f · · ·
∑ ∑
nf =0 f nf =0
Grupisanjem dobijamo
∑
∞ ∏ n
af f = (1 + a1 + a21 + · · · )(1 + a2 + a22 + · · · ) · · · (1 + af + a2f + · · · ) · · ·
∑
nf =0 f
∑
∞ ∏ n
∑ ∑ ∑ n
∏∑
∞
n
af f = an1 1 an2 2 ··· af f ··· = af f .
∑
nf =0 f n1 n2 nf f nf =0
∑
1 ∏( )
e−β(ϵf −µ)nf = 1 + e−β(ϵf −µ) ⇒ Ξ= 1 + e−β(ϵf −µ) .
nf =0 f
I. Bose-Einstein-ova statistika
m
Tc
T
Slika 7.1:
∑ ( )
Za Bose idealan gas imamo da je Ω = kB T f ln 1 − e−β(ϵf −µ) . Kako je ϵf ≡ ϵp⃗ =
⃗2
p
2m
, pa je
∑ ( )
Ω = kB T ln 1 − e−β(ϵp⃗ −µ) .
p
⃗
∫∞ ( )
Integral I = ± √4π x2 ln 1 ∓ ze−x dx ne može da se izračuna u opštem slučaju,
2
0
Za naš slučaj to je
∫ ∞ ∫
4 ∑ (±1)n+1 z n ∞ 2 −nx2
∞
4 (±1)n+1 n −nx2
I = −√ z e dx = − √ xe dx.
π 0 n π n=1 n 0
∫∞ √
x2 e−nx dx = · ·
2 1 1 π
Budući da je 0 2 3 2
, imaćemo
n2
∑
∞
zn
I=− (±1)n+1 5 .
n=1 n2
pa se vidi da je I = −G±
5/2 (z). Ovde uvedena klasa funkcija je u vezi sa Riemann-
ovom ζ-funkcijom
1 1
ζ(s) = 1 + s
+ s + ···
2 3
Funkcije G±
s (z) i Riemann-ova funkcija ζ(s) su vezane
1 ±
p = kB T G (z).
λ3T 5/2
Unutrašnja energija je
( )
∂
U (z, T, V ) = − ln Ξ ,
∂β z,V
I.. BOSE-EINSTEIN-OVA STATISTIKA 257
3
U = pV.
2
Da bi smo našli termičku jednačinu stanja potreban nam je srednji broj čestica,
( )
∂Ω ∂Ω ∂z
⟨N ⟩ = − =− · .
∂µ T,V ∂z ∂µ
Može se pokazati da je
∂ ±
G (z) = −G±
3/2 (z),
∂z 5/2
-
te se sredivanjem dobija
V ±
⟨N ⟩ = G (z).
λ3T 3/2
Pritisak je ranije dobijen da je
( )
kB T ∑
∞ n
1 ± n+1 z kB T 1 2 1 3
p = kB T 3 G5/2 (z) = 3 (±1) 5 = · z ± 5 z + 5 z + ··· .
λT λT n=1 n2 λ3T 22 32
λ3T ⟨N ⟩ ∑ ∞ n
n+1 z 1 2 1 3
= G±
3/2 (z) = (±1) 3 = z ± 3 z + 3 z + ···
V n=1 n 2 2 2 3 2
λ3T ⟨N ⟩
Označimo ξ = V
i izrazimo pritisak preko te veličine.
Može se pokazati da ako imamo dva polinoma oblika
ϕ(y) = x = Ay + By 2 + Cy 3 + · · ·
258 § 7. OSNOVE KVANTNE STATISTIČKE FIZIKE
možemo pisati
1 b 2b2 − ac
A= ; B=− ; C= .
a a3 a5
1 1 1 1
ξ=z± 3 z2 + 3 + ··· a = 1; b=± 3 ; c= 3 ;
2 2 3 2 2 2 32
te stoga
1 1 1 1
A = 1; B=∓ 3 ; C= − 3; z=ξ∓ 3 ξ2 + · · ·
2 2 2 32 2 2
Konačno, pritisak je
( ( )2 )
⟨N ⟩ kB T 1 1 1
p= · 1∓ 3 ξ± 5 ξ 1∓ 3 ξ + ···
V 2 2 2 2 22
3
pV = ⟨N ⟩ kB T ⇒ U= ⟨N ⟩ kB T
2
Na temperaturi T = Tc je µ = 0, pa imamo
∑ ∫ ∫ ∞
1 V 1 3⃗ 4πV 1
N= βϵf − 1
= 3 βϵ − 1
d k = 3 βϵ − 1
k 2 dk.
f
e (2π) e (2π) 0 e
~2 k2
Kako je ϵ = 2m
, biće
√ 3 ∫ ∞ √
2 V m2 ϵ
N= · dϵ.
2 2π 2 ~3 0 eβϵ −1
N= N (ϵ = 0) + N (ϵ ̸= 0) ,
| {z } | {z }
broj bozona broj bozona koji
u osnovnom stanju nisu u osnovnom stanju
1
no = , (7.4)
eβ(µo H+µ) − 1
1 1
µ = −µo H − ln(1 + ), (7.5)
β no
∑ 1 ∑ 1 ∑ 1
N= p2
+ p2
+ p2
.
−µo H+µo H)
p
⃗
β( 2m
e −1 p
⃗ e β( 2m +0+µo H)
−1 p
⃗
β( 2m
e +µo H+µo H)
−1
(7.6)
µo H
Prelazeći na termodinamički limes, i koristeći uslov kB T
>> 1, dobijamo
2π 3 3 3 − µo H − 2µ oH
n= (2mk B T ) 2 Γ( )[ζ( ) + e kB T + e kB T
]. (7.7)
h3 2 2
µo H
Gde, takodje, možemo iskoristiti da je kB T
>> 1, što dovodi do
N h3 2 1
Tc = [ 3 ]
3 . (7.8)
V ζ( 2 ) 2πmkB
Drugo rešenje:
Hamiltonijan jednog bozona je
H = cp −σµ0 H.
|{z}
ϵ
262 § 7. OSNOVE KVANTNE STATISTIČKE FIZIKE
1
N0 = ,
eβ(−µ0 H−µ) −1
odakle sledi ( )
1
µ = −µ0 H − kB T ln 1 + .
N0
Ako je T ≈ 0 ⇒ N0 ≈ N , tada
µ ≈ −µ0 H.
(a) H = 0 ∫
4πV ∑ ∞ p2 dp ∑
N = N0 + 3 β(ϵ−µ) − 1
, 1 = 3,
h σ 0 e σ
pa je ∫ ∞
12πV p2 dp
N = N0 + .
h3 0 eβ(ϵ−µ) − 1
Za T = Tc ⇒ µ=0 ⇒ N0 = 0,
∫
N 12π ∞ p2 dp
n= = 3 .
V h 0 eβ(ϵ−µ) − 1
(b) H ̸= 0
∑ 4πV ∫ ∞
p2 dp
N= ,
σ
h3 0 eβ(ϵ−σµ0 H+µ0 H) − 1
a kako je σ = {−1, 0, 1}
∫ ∞ 2 ∫ ∞ ∫ ∞
4πV
p dp p2 dp p2 dp
.
N= 3 + +
h eβϵ − 1 eβ(ϵ+µ0 H) − 1 eβ(ϵ+2µ0 H) − 1
| 0
{z } | 0
{z } | 0
{z }
I1 I2 I3
(a) 2D gas
∫ ∫ ∞
gS d2⃗k gS k dk
N = N0 + = N0 + 2π .
(2π)2 e −1
βϵ (2π)2 0 eβϵ−1
~2 k 2
Kako je ϵ = 2m
dolazimo do
∫ ∞
gS m dϵ
N = N0 + 2π 2 = ∞,
(2π)2 ~ 0 eβϵ−1
odnosno
N
n= → ∞.
S
(b) 1D gas
∫ ∞ ∫ ∞
gL dk gL dk
N = N0 + = N0 + 2 = ∞,
2π −∞ e −1
βϵ 2π 0 e
2 k2
β ~2m
−1
tj.
N
n= → ∞.
L
Rešenje:
Bose kondenzacija se javlja kada je temperatura Bose gasa manja od kritične tem-
perature Tc . Tada je linijski potencijal jednak nuli. Ako pak postoji magnetno
polje, hemijski potencijal mora da bude manji od nule. Ukupan broj čestica je
∑ 1 ∑V ∫ d3 p⃗
N= ( ) = N0 + .
2
β p −σµ0 H−µ h3 eβ(ϵ−µ0 σH−µ) − 1
⃗,σ e 2m
p −1 σ
I.. BOSE-EINSTEIN-OVA STATISTIKA 265
µ ≈ −µ0 H.
(a) H = 0
∑ 4πV ∫ ∞
p2 dp ∑
N = N0 + , 1=3
σ
h3 0 eβ(ϵ−µ) − 1 σ
∫ ∞ 2 ∫ ∞
12πV p dp 12πV p2 dp
N = N0 + ( ) = N0 + ( ) .
p2 p2
h3 β −µ h3 β +µ0 H
0 e 2m
−1 0 e 2m
−1
Kako je H = 0, ∫ ∞
12πV p2 dp
N = N0 + p2
,
h3 0 eβ 2m − 1
na T → Tc je µ = 0 ⇒ N0 = 0, imamo
∫
12π ∞ p2 dp
n= 3 p2
,
h 0 eβ 2m − 1
rešavajući integral biće
√ ( )3 ( ) ( )
12 2π m 2 3 3
n= 3
Γ ζ ,
h β 2 2
odakle sledi ( ) 23
1 nh3
Tc = ( ) .
2πmkB 3 ζ 32
(b) H ̸= 0; µ0 H ≫ kB T
∑ 4πV ∫ ∞
p2 dp
N= ( ) ,
p2
h3 β +µ0 H(σ+1)
σ 0 e 2m
−1
∫ ∞ ∫ ∞ ∫ ∞
4πV p 2
dp p 2
dp p 2
dp
N= 3 + ( 2 ) + ( 2 ) .
h p 2 p p
| 0 e{z − 1}
β 2m 0 e
β 2m
+µ 0 H
− 1 0 e
β 2m
+2µ 0 H
− 1
| {z } | {z }
I1 I2 I3
266 § 7. OSNOVE KVANTNE STATISTIČKE FIZIKE
Rešavanjem imamo
√ ( ) 32 ( )
2π m 3
I1 = ζ ;
2 β 2
√ ( ) 32
2π m
I2 = e−βµ0 H ;
2 β
√ ( ) 32
2π m
I3 = e−2βµ0 H .
2 β
što daje
( ) 32
1 nh3
Tc = ( ) .
2πmkB ζ 23
~2 k2
8. Razmotriti trodimenzionalni gas bozona sa disperzionom relacijom ϵ = 2m
.
Odrediti kritičnu temperaturu Tc Boze kondezacije, i diskontinuitet toplotnog ka-
( v)
paciteta ∂C
∂T V
na toj temperaturi.
Rešenje:
Odredimo kritičnu temperaturu. Da bi makroskopski deo čestica prešao u kon-
Ncond
denzat potrebno je da Ne
bude nenulto, gde je Ncond broj bozona u stanju sa
minimalnom energijom (kondezovani), a Ne broj bozona u preostalim, pobud-
jenim energetskim nivoima. Pošto je minimalna energija za datu disperzionu
relaciju ϵ = 0 za k = 0, treba da važi ξ
(1−ξ)N
̸= 0 za makroskopske uzorke N → ∞,
što je jedino moguće ako je ξ → 1 kada T → Tc . Broj atoma Ne u pobudjenim stan-
jima, na kritičnoj temperaturi (granici kada dolazi do Boze kondezacije) jednak je
I.. BOSE-EINSTEIN-OVA STATISTIKA 267
~2 k 2 1
Nakon uvedjenja smene x ≡ 2m kb T
dobijamo:
∫ ∞
V x3/2−1 dx
N = 3 2π (2 kB T )3/2
h ex − 1
|0 {z }
ζ(3/2) Γ(3/2)
(2π m kB T )3/2
= V ζ(3/2),
h3
(a)
−
f5/2 (z)
pV = N kB T − (7.11)
f3/2 (z)
(b)
( ) −
1 ∂z 5 f5/2 (z)
=− − . (7.12)
z ∂T P 2T f3/2 (z)
∫∞
Koristiti da je Γ(s)fs− (z) = 0
xs−1 dx
z −1 ex −1
, gde je Γ(s)- gama funkcija, dok je fs− (z)
Riemann-ova zeta funkcija.
Rešenje:
∑ -
pov́ sto je suma po spinovima σ 1 = g degeneracija ⃗k-tog nivoa i ako predemo
u termodinamički limes, dobijamo
∫
V d3 k
N = g 3 (7.13)
2π eβ(ϵ−µ) − 1
∫ ∞
V k2 dk
= g 3 4π , (7.14)
2π 0 z eβϵ − 1
−1
~2 k2
gde je z = eβµ fugacitivnost i ϵ = 2m
. Prelaskom na integraciju po ϵ dobijamo
N g −
= 3 f3/2 (z), (7.15)
V λ
h
gde je λ = 1 . Koristeći da u trodimenzionom slučaju idelanog gasa
(2πkB T ) 2
imamo
p g −
= 3 f5/2 (z). (7.16)
kB T λ
270 § 7. OSNOVE KVANTNE STATISTIČKE FIZIKE
2 2
( )
5 gkB − gkB T − ∂z
0= 3
f5/2 (z) + 3
f3/2 (z) , (7.18)
2 λ λ ∂T p
10. Razmatrati gas neinteragujućih identičnih Bose kvantnih čestica, čija je energijska
disperziona relacija istog oblika kao u prethodnom zadatku. Odrediti minimalnu
dimenziju sistema d(s) na kojoj se mo ze zapaziti Bose-Einsein-ova kondenzacija.
12. Razmatrati idealni Bose gas u trodimenzionom sudu u blizini kritične tempera-
( v)
ture T = Tc . Odrediti vrednost diskontinuiteta veličine ∂C
∂T V
na kritičnoj tem-
peraturi.
-
gde je g degeneracija. Uvodenjem sfernih koordinata imamo
∫ ∞
2 V ϵk 2 dk
pV = g 4π ,
3 (2π)3 0 eβ(ϵ−µ) + 1
~2 k 2
a kako je ϵ = 2m
, biće
√ 3 ∫ ∞ 3
2 gV 2m 2 ϵ 2 dϵ
pV = 4π .
3 (2π)3 ~3 0 eβ(ϵ−µ) + 1
Integral oblika
∫ ∞ 3
ϵ 2 dϵ
I= β(ϵ−µ)
0 e +1
272 § 7. OSNOVE KVANTNE STATISTIČKE FIZIKE
pa je
(∫ ∫ ∫ ∫ )
3 3 ∞ 3
1 βµ
3
βµ
(βµ − x) 2 βµ
(x + βµ) 2 (x + βµ) 2
I= 5 (βµ − x) dx −2 dx + dx + dx .
β2 0 0 ex + 1 0 ex + 1 βµ ex + 1
biće ( )
∫ ∞
1 2 5 3 3 x dx
I= 5 (βµ) 2 + (βµ) 2 + ··· .
β2 5 βµ 0 ex + 1
Kako je ( )
∫ ∞
x dx 1 1 π2
= Γ(2) 1 − ζ(2) = · ,
0 ex + 1 2 2 6
II.. FERMI-DIRAC-OVA RASPODELA 273
imamo konačno
( )
1 2 5 3 3 π2
I= 5 · (βµ) 2 + (βµ) 2 · · + ··· ,
β2 5 βµ 12
pa je
( ) 32 ( ( ) 52 ( ) 12 )
2 V 2m 25 µ π2 µ
pV = g 2 (kB T ) · 2 + + ··· .
3 4π ~2 5 kB T 4 kB T
( )
Kako važi Ω = −pV , a N = − ∂Ω
∂µ
, biće
T,V
( ) 32 ( 3
)
gV 2m 5 µ2 π2 1
N= 2 (kB T ) 2
5 + · 1 + ··· .
6π ~2 (kB T ) 2 8 µ 2 (kB T ) 12
Sa druge strane je
∑ ∫ ∫ kF
1 V 3⃗ V V 3
NT =0 = = g· d k θ(kF −k) = g ·4π k 2 dk = g ·k .
eβ(ϵ−µ)+1 (2π)3 (2π)3 0 6π 2 F
p
⃗,σ
Ako je ϵF = µ0 , tad
~2 kF2 1√
= µ0 ⇒ kF = 2mµ0 , µ0 = µ(T = 0)
2m ~
Pretpostavimo da je
( ( )B )
kB T
µ = µ0 1+A + ··· ,
µ0
274 § 7. OSNOVE KVANTNE STATISTIČKE FIZIKE
π2 π2
0 = A(kB T )B µ1−B
0 + (kB T )2 µ−1
0 − A · (kB T )B+2 µ−1−B
0 +···
12 | 12 {z }
≈0
mora biti 1 − B = −1 ⇒ B = 2, te
1 π2 1 π2
0 = A(kB T )2 + (kB T )2 ⇒ A=− ,
µ0 12 µ0 12
pa je konačno ( )
( )2
π2 kB T
µ(T ) = µ0 · 1− + ··· .
12 µ0
17. Izračunati razliku izmedju toplotnih kapaciteta cp i cv za idealni Fermi gas čestica
spina 1/2 u trodimenzionom slučaju na temperaturama bliskim 0.
↑ T →0⇒ µ
→∞
∫ ∞
kB T
f (µ + kB T z)dz
+ kB T =
ez + 1
∫ µ ∫ 0
∞
f (µ + kB T z) − f (µ − kB T z)
= dεf (ε) + kB T dz ≈
0 0 ez + 1
∫ µ ∫ ∞
2 ′ z
≈ dεf (ε) + 2(kB T ) f (µ) dz z +
e +1
0
|0 {z }
π2
∫
2·6
∞
(kB T )4 ′′′ z3
+ f (µ) dz z +o(T 5 ) . (7.22)
3 e +1
| 0
{z }
3·7π 4
360
∫∞ s−1 ( ) -
Gde je iskorišćena jednakost 0
dz ezaz +1 = Γ(s)
as
1− 2s−1
1
ζ(s) .1 Neznatnim sredivanjem
jednakosti (7.22) dobijamo
∫ ∞ ∫ µ
dεf (ε) π2 2 ′ 7π 4
1 ≈ dεf (ε) + (kB T ) f (µ) + (kB T )4 f ′′′ (µ) + o(T 4 ) .
0 e kB T
(ε−µ)
+1 0 3 360
(7.23)
1
20. Odrediti hemijski potencijal jednodimenzionog fermionskog gasa, spina 2
i mag-
netnog momenta u zavisnosti od spoljašnjeg magnetnog polja H na temperaturi
T = 0.
Rešenje:
Energija fermiona je
p2
− σµo H.
ϵp,σ = (7.24)
2m
Pošto je spin σ ∈ {1/2, −1/2}, broj fermiona u sistemu je
1
Jednakost poznata iz Statističke fizike - Uvod u kvantnu statističku fiziku. Dokazuje se jednos-
tavno, odgovarajućim razvojem u red i integracijom ”član po član”.
276 § 7. OSNOVE KVANTNE STATISTIČKE FIZIKE
ε ε
εk
εk +µB B
εF εF
εk−µB B g (ε)
∑ ∑ 1 ∑ 1
N= np⃗,σ = ϵ+ −µ + ϵ− µ , (7.25)
p,σ p
⃗ e kB T
+1 p
⃗ e kB T
+1
gde je µ hemijski potencijal sistema, T temperatura sistema, dok su ϵ+ i ϵ− energije
fermiona koji ima spin + 21 , odnosno − 12 . Ako predjemo u termodinamički limes
i zamenimo za ϵ+ , tj ϵ− dobijamo
∫ ∞ ∫ ∞
2L dp dp
N= { p2
+ p2
}. (7.26)
h 0 µo
2m − 2 H−µ 0
µo
2m + 2 H−µ
e kB T +1 e kB T +1
Prethodni izraz ćemo izračunati, tako što najpre predjemo na integraciju po ϵ,
pa zatim integral izračunamo za T = 0. U tom graničnomslučaju dobićemo
raspodelu broja v cestica po energijama u obliku Heviside-ove Θ funkcije, tako
da se prethodna relacija svodi na
∫ ∞ ∫ ∞
2L dϵ p2 µo dϵ p2 µo
N= { √ Θ(− + H + µ) + √ Θ(− − H + µ)} = (7.27)
h 0 ϵ 2m 2 0 ϵ 2m 2
II.. FERMI-DIRAC-OVA RASPODELA 277
2L √ µo 1 µo 1
= 2m{(µ(T = 0, H) + H) 2 + (µ(T − 0, H) − H) 2 }. (7.28)
h 2 2
U slučaju slabog spoljašnjeg magnetnog polja, gornje izraze možemo da razvi-
µo H N
jemo po 2µ(0,H)
. Takodje ako uvedemo koncentraciju fermiona n = V
dobijamo
izraz
√ 1 4 1 µo H 2
n= 2mµ(0, H) 2 {1 − ( ) }. (7.29)
h 16 2µ(0, H)
Koncentraciju broja fermiona možemo prikazati i preko hemijskog potencijala za
√ √
T = 0 i H = 0. Tako dobijamo da je n = h4 2m µ(0, 0), koji možemo da koris-
timo prilikom zamene koncentracije fermiona u okviru sistema, tako da nizom
jednstavnih matematičkih transformacija dobijamo da hemijski potencijal fermi
sistema za H ̸= 0 se menja sa uključenjem magnetnog polja na sledeći način:
1 µo H 2
µ(0, H) = µ(0, 0){1 + ( ) }. (7.30)
32 2µ(0, H)
µ0
− 02 −µ(0,0)
∑ π2 ∑
β 2m
1 kB T βe
M= ( ) − µ(0, 0) ( ( 2 ) )2 −
2L p2 µ H
− 02 −µ(0,0) µ(0, 0) 12 β p
−
µ0 H
−µ(0,0)
p
⃗ e 2m
+1 p
⃗ e 2m 2 +1
( )
∑ ( )2 p2
µ H
+ 02 −µ(0,0)
π2 ∑
β 2m
1 kB T βe
− ( ) − µ(0, 0) ( ( 2 ) )2 .
p2 µ H
+ 02 −µ(0,0) µ(0, 0) 12 µ H
β p + 0 −µ(0,0)
p
⃗ e 2m
+1 p
⃗ e 2m 2 +1
te je
( ) 21 ∫ βµ ( )
2L m ex dx − 12 − 12
I(µ) = · (βµ − x) + (βµ + x) ,
h 2(kB T )5 0 (ex + 1)2
√ ∫ ( ( )2 ( )2 )
βµ
2Lβ 2 m ex dx 1 x 3 x 1 x 3 x
I(µ) = · 1+ + +1− + ,
h 2µ 0 (ex + 1)2 2 βµ 8 βµ 2 βµ 8 βµ
√ ∫ ( )
4Lβ 2 m βµ ex dx 3 x2
I(µ) = · 1+
h 2µ 0 (ex + 1)2 8 (βµ)2
što na kraju daje √ ( )
4Lβ 2 m 1 π2 3 1
I(µ) = + .
h 2µ 2 6 8 (βµ)2
( )
Za slučaj I µ(0, 0) ± µ02H , imamo
( ) √ ( )− 12 ( ( )2 )
µ0 H 2Lβ 2 m µ0 H π2 µ0 H
I µ(0, 0) ± = µ(0, 0) ± · 1 + 2 µ(0, 0) ± ,
2 h 2 2 8β 2
a razvijanjem u red
( ) √ ( )( ( ))
µ0 H 2Lβ 2 m 1 µ0 H π2 µ0 H
I µ(0, 0) ± = √ 1∓ · 1+ 1∓ ,
2 h 2 µ(0, 0) 4µ(0, 0) 8(βµ(0, 0))2 µ(0, 0)
dolazimo do
( ) √
µ0 H 2Lβ 2 m 1 µ0 H
I µ(0, 0) ± = √ (1 ∓ +
2 h 2 µ(0, 0) 2µ(0, 0)
)
π2 π 2 µ0 H π 2 µ0 H
+ ∓ ∓ .
8(βµ(0, 0))2 32β 2 µ(0, 0)3 8β 2 µ(0, 0)3
∂N N
22. Pokazati da je u 2D slučaju za fermione ispunjeno ∂T
= T
.
23. Razmatrati sud sa pregradom, koji je prikazan na slici. Zapremina na levoj strani
3
suda je 4
ukupne zapremine suda. Deo suda s leve strane sadrži fermione, spina
1
2
, na temperaturi T = 0 sa Fermi-jevom energijom od 1 eV. Desna strana suda je
prazna. Sistem je potpuno izolovan. U trenutku t = 0 pregradni zid je pomeren
280 § 7. OSNOVE KVANTNE STATISTIČKE FIZIKE
kB T mc2 1
Cv = N kπ 2 (1 + )2 . (7.31)
~kF ~kF
tada
∂N ∂N ∂β 1 ∂N
= · =− 2
· ,
∂T ∂β ∂T kB T ∂β
a drugi činioc je
∫ ( )
∂N gS ∂ ∞
k dk gS ∂ m ( )
= · = · ln 1 + e βµ
,
∂β 2π ∂β 0 eβ(ϵ−µ) + 1 2π ∂β ~2 β
( )
∂N gS m ( ) gS m µeβµ
= − 2 2 ln 1 + e βµ
+ · 2 · ,
∂β 2π ~β 2π ~ β 1 + eβµ
II.. FERMI-DIRAC-OVA RASPODELA 281
dok je
gS m ( )
N= · 2 ln 1 + eβµ .
2π ~ β
Za visoke temperature će važiti N ≈ gS
2π
· m βµ
~2 β
e ,
∂N ∂N
= −kB β 2 ,
∂T ∂β
gde je
∂N gS m gS m ( )
≈− · 2 2 eβµ + · 2 · µeβµ 1 − eβµ ,
∂β 2π ~ β 2π ~ β
∂N gSm gSm gSm gSmµ
≈− 2 2
(1 + βµ) + 2
µ(1 + βµ)(−βµ) = − 2 2
− .
∂β 2π~ β 2π~ β 2π~ β 2π~2 β
Rešenje:
Velika particiona funckija je data izrazom
∏
Ξ= (1 − e−β(ϵg −µ) )−1 . (7.33)
g
∫ ∞ d
−β(ks −µ) π2
Ω = kB T d ln(1 − e ) d k d−1 dk. (7.36)
0 Γ( 2 + 1)
V Γ( ds + 1) eβµ 1
Ω=− d d d f d +1 (e
βµ
). (7.37)
(4π) 2 Γ( 2 + 1) β s s
Γ( ds + 1) d eβµ 1
V βµ
E= d d d f d +1 (e ) (7.40)
s
(4π) 2 Γ( 2 + 1) β s s
dobijamo da je
E d
= . (7.41)
pV s
Rešenje:
Polazeći od formule za srednji broj fermiona
∫ ∞ 1
N 2gπ(2kB T )3/2 x2
= , (7.42)
V h3 0 z −1 ex + 1
∑∞ ∫ ∞
− 2s
e−t dt,
s−1
r
= (−1) (r + 1) t 2
r=0 0
∑
∞
(−1)r (r + 1)− 2 .
s
I(s) = Γ(s) (7.44)
r=0
∑
∞
(−1)r (r + 1)− 2 = (1− 2 + 3− 2 + 5− 2 + 7− 2 + · · · + (2r − 1)− 2 ) +
s s s s s s
r=0
+ (2− 2 + 4− 2 + 6− 2 + 8− 2 + · · · + (2r)− 2 ) =
s s s s s
∑∞
−s
∑
∞
− 2s
(r)− 2 =
s
= (r) − 2 2
+1
r=1 r=1
−s s
= (1 − 2 2
+1
)ζ( ),
2
284 § 7. OSNOVE KVANTNE STATISTIČKE FIZIKE
2gπ(2kB T0 )3/2 3 −1 3
n= 3
Γ( )(1 − 2 2 )ζ( ). (7.45)
h 2 2
∑ ∫ kF
V
N = g n̄k = g 2 k 2 dk =
k
2π 0
( )3
V kF3 g 2mεF 2
= g 2 = 2 . (7.46)
2π 3 6π ~2
4πϵF 1
To = √ . (7.47)
kB 6π 2 (1 − 2 )ζ( 3 )
2 2
29. Odrediti razliku toplotnih kapaciteta Cp − Cv za idealni Fermi gas čestica spina
1
2
u trodimenzionom slučaju na niskim temperaturama. Uzeti da je disperziona
relacija ε = Ak s , gde su A i s konstante.
N
a. Odrediti izraz za Fermi-jev talasni broj kF za taj gas u funkciji gustine n = V
.
2
Videti na kraju zbirke Dodatak C deo.
II.. FERMI-DIRAC-OVA RASPODELA 285
~kF
gde je xF = mc
.
31. Pokazati da u slučaju idealnog Fermi gasa važe za izotermsku i adijabatsku kom-
presibilnost sledeće relacije:
+
1 ζ1/2 (z)
KT = + ,
nkT ζ3/2 (z)
+
3 ζ3/2 (z)
KS = + ,
5nkT ζ5/2 (z)
∫∞ xj−1 dx
gde je n koncentracija broja čestica. Sa ζj+ (z) = 1
Γ(j) 0 z −1 ex +1
je dat tzv. general-
isani logaritamski red. Pokazati da na niskim temperaturama važe relacije:
[ ( )2 ]
3 π 2 kB T
KT ≈ 1− ,
2nϵF 12 ϵF
[ ( )2 ]
2
3 5π kB T
KS ≈ 1− .
5nϵF 12 ϵF
286 § 7. OSNOVE KVANTNE STATISTIČKE FIZIKE
M = µ0 n+ − µ0 n−
[ d ∫ ∞
L 1
= µ0 d
d Cd k d−1 dk · β(ϵ+ −µ(H,T ))
(2π) 0 e +1
d ∫ ∞
]
L 1
− d Cd k dk · β(ϵ− −µ(H,T ))
d−1
,
(2π)d 0 e +1
( 2π )d
gde je L
zapremina jednog stanja sistema u d-dimenzionalnom k prostoru, a
Cd zapremina sfere jediničnog radijusa u d dimenzionalnom prostoru. Hoćemo
da u uračunamo prvu korekciju koju unose magnetno polje i temperatura u broj
popunjenosti stanja. Već vidimo da magnetno polje utiče na energiju sistema, ali
želimo da vidimo i kako utiče na hemijski potencijal. Podjimo od
( )d [∫ ∞
L k d−1 dk
N = d Cd
2π eβ(ϵ+ −µ(H,T )) + 1
∫ ∞
0
]
k d−1 dk
+ . (7.50)
0 eβ(ϵ− −µ(H,T )) + 1
II.. FERMI-DIRAC-OVA RASPODELA 287
( √ )d ∫ ∞
L 2m 1 dϵ
N = d Cd ϵd/2−1 β(ϵ−µ)+1 .
2π ~ 2 e
|0 {z }
≡I
∫ ∞ ( )d/2−1
x 1 1
I = +µ dx
−βµ β ex + 1 β
∫ βµ ( )d/2−1 x ∫ ∞
x e +1−11
= µ− dx + ...
0 β 1 + ex β 0
∫ βµ ( )d/2−1
x 1
= µ− dx
0 β β
( )d/2−1 ( )d/2−1
∫ βµ µ + x − µ − x ∫ ∞
β β ...
+ x
dx + x
dx
0 e +1 βµ e + 1
∫ ∞
d/2 2 µd/2−2 x dx
≃ µ + (d − 2) 2
d β ex + 1
| 0 {z }
1
Γ(2)ζ(2)(1− )
22−1
( )
d/2 d π 2 (d − 2) kB
2 2
T
= µ + ,
2 12 µ2
Sada se vraćamo na početan izraz za magnetizaciju sistema. Kao što smo ranije
zaključili, nije potrebno u prvoj korekciji srednjeg broja popunjenosti uračunavati
promenu hemijskog potencijala usled postojanja magnetnog polja, ali treba uračunati
promenu hemijskog potencijala usled nenulte temperature. Razvijanjem funkcije
srednjeg broja popunjenosti u red oko µ(0, 0) dobijamo:
( √ )d [∫ (
∞
L 2m 1 d/2 1
M = µ0 K dK +
2π ~ 2 0 eβ(K−µ0 H/2−µ(0,0)) +1
( )2 )
βeβ(ϵ−µ0 H/2−µ(0,0)) π 2 (2 − d) kB T
+ β(ϵ−µ0 H/2−µ(0,0)) · µ(0, 0)
e +1 12 µ(0, 0)
+slično sa drugim integralom]
34. Stanje termodinamičke ravnoteže kristalne rešetke od N atoma, koja osciluje, pri-
bližno je ekvivalentno stanju ravnoteže skupa od 3N neinteragujućih kvantnome-
haničkih oscilatora. Prema tzv. Debajevoj aproksimaciji postoji: dnω = ω9N3 ω 2 dω
D
-
za ω < ωD i dnω = 0 za ω > ωD , oscilatora sa učestanošću izmedu ω i ω + dω.
- ∫ω
Konstanta ωD , tzv. Debajeva učestanost odredena je uslovom: 0 D dnω = 3N .
Prihvatajući ovu aproksimaciju odrediti toplotni kapacitet kristalne rešetke pri
veoma niskim i veoma visokim temperaturama.
35. Pri faznom prelazu gasa elektrona iz normalnog u superprovodno stanje, e− us-
trojavaju se po energijskim nivoima
~2 (k − kF )2
ϵk = ∆ + , (7.53)
2m∗
III.. PRIMENA KVANTNIH RASPODELA 291
36. Uzeti da je fononski deo spektra čvrstog tela dat relacijom ϵ = Ak s . Koristeći
Debay-evu aproksimaciju odrediti toplotni kapacitet trodimenzionog čvrstog tela.
- odrediti gustinu stanja, i pokazati da se u limesu niskih temperatura
Takode
3
toplotni kapacitet ponaša u skladu sa relacijom Cv ≈ T s .
Rešenje:
Uzimajući u obzir da je hemijski potencijal fonona jednak nuli i da su bozoni, za
f(ω)
ωD ω
Slika 7.3: Raspodela normalnih moda f (ω) za aluminijum. Puna linija prikazuje
rezultate dobijene rasejanjem x-zraka, dok isprekidana linija prikazuje raspodelu
normalnih moda u Debay-evoj aproksimaciji.
V 2
g(k) = g k , (7.55)
2π 2
V − 3 3 −1 1
g(ϵ) = g A s εs . (7.56)
2π 2 s
gde smo uveli smenu x = βε. U cilju da pojednostavimo zapis, uvešćemo oznaku
∫ xm
TD x3/s dx
D(xm ) = D( ) = g , (7.61)
T 0 ex − 1
εD
gde je definisana tzv. Debay-eva temeratura TD = kB
. Za ukupni toplotni ka-
pacitet dobijamo
( )
∂E
Cv = = (7.62)
∂T V
[ ( )]
gV 1 −3/s 1+ 3s 3 3 TD 3
′ TD TD
= A kB ( + 1)T s D( ) + T s D ( ) − 2
+1
. (7.63)
2π 2 s s T T T
III.. PRIMENA KVANTNIH RASPODELA 293
38. Disperziona relacija talasnog kretanja unutar čvrstog stanja je data relacijom ω =
αk n . Odrediti specifičnu toplotu na visokim temperaturama.
-
39. Grafen je pojedinačan list (sloj) ugljenikovih atoma, koji su medusobono spo-
jeni u dvodimenzionoj ravni. On se može dobiti skidanjem (ljuštenjem) grafita.
Struktura provodnih zona grafena su takve da ekscitacije pojedinačnih čestica se
mogu ponašati kao relativistički Dirac-ovi fermioni, čija je disperziona jednačina
na nisikim energijama data jednačinom
a. Pokazati da je u ravnoteži
µ2 = 2µ1 , (7.67)
A + A → A2 . (7.68)
ature T = 0. I na kraju, odrediti odnos izmedu- tih dvaju vrsta čestica i za slučaj
42. Elektroni u nekom komadu metala bakra se mogu razmatrati kao Fermi gas.
gr
Metal bakra u čvrstom stanju ima masenu gustinu 9 cm3 . Pretpostavljajući da
je T = 0, izračunati Fermi-jevu energiju εF i Fermi-jevu temperaturu TF = kεFB .
43. He3 na niskim temperaturama može preći iz tečnog u čvrsto stanje promenom
( dp )
pritiska. Posebna osobina te granice faza je u tome da dT topljenja
je negativno na
( dp )
temperaturama od 0.3K (dok je na temeraturama oko 0.1K dostiže čak dT topljenja
≈
−3 · 106 PKa .
a. U čvrstoj fazi, He3 formira kristalnu rešetku. Spin nukleusa svakg atma je
1 - spinova, odrediti entropiju ss po atomu
. Zanemarujući interakciju izmedu
2
III.. PRIMENA KVANTNIH RASPODELA 295
(koristeći jedino ”spinske” stepene slobode). Svaki atom ima nuklearni spin 21 .
b. Tečni He3 je modelovan kao idealni gas fermiona sa zapreminom od 46·10−30 m−3
po atomu. Odrediti kolika je Fermijeva temperatura TF , u stepenima Kelvina?
I. Gama funkcija
∫∞
Γ(n + 1) = 0
xn e−x dx, (8.1)
π
Γ(n)Γ(1 − n) = zan ∈ (0, 1); (8.5)
sin(πx)
√
1 π
Γ(n)Γ(n + ) = 2n−1 Γ(2n). (8.6)
2 2
√
Γ(n) ∼ nn e−n 2πn. (8.7)
II.. O RIEMANN-OVOJ ζ FUNKCIJI 297
i predstavlja specijalan slučaj tzv. Hurwitz-ove ζ funkcije, koja se uvodi kao ζH (s, q) =
∑∞ 1
n=1 (n+q)s . Vrednosti Riemann-ove funkcije se dobijaju kao specijalan slučaj, odnosno
-
predemo sa integracije po x na integraciju po sx, dobijamo
∫ ∞
Γ(n) = s n
e−sx xn−1 dx. (8.11)
0
3 5 7
s 2
2 2
3 2
4 5 6 7
π2 π4 π6
ζ(s) 2.612 6
1.341 1.202 1.127 90
1.037 945
1.008
d
Zahtevajući da je dx
(xJν (x)) = xJν−1 (x), dobijamo da je ν = 1, što dovodi do izraza
d
xJo (x) = (xJ1 (x)). (8.17)
dx
Ovome je ekvalitetno ∫ x
zJo (z)dz = xJ1 (x). (8.18)
0
1 ( x )1 x
J1 (x) ≈ = . (8.20)
Γ(1 + 1) 2 2
[1] Dalvit D.A.R, Frastai J. i Lawrie I.D. Problems on Statistical Mechanics, Institute
of Physics Publishing, Bristol and Philadelphia, 1999.
[3] Milić B, Milošević S. i Dobrosavljević Lj. Zbirka Zadataka iz Teorijske Fizike III deo:
Statistička Fizika, Naučna Knjiga, Beograd 1979.
[7] Gibbs J.W. Elementary Principles in Statistical Mechanics, Dover Publication, New
York 1960 (Reprint of original 1902. edition).
[8] Hill T.L. Statistical Mechanics: Principles and Selected Applications, McGraw-Hill,
New York 1956.
BIBLIOGRAPHY 301
[9] Huang K. Statistical mechanics, 2nd ed., John Wiley, New York 1987.
[10] Huang K. Introduction to Statistical Physics, Taylor & Francis, London 2001.
[11] Isihara A. Statistical Physics, Academic Press, New York – London 1971.
[13] Landau L.D. and Lifshitz E.M. Statistical Physics, Addison-Wesley 1969.
[18] Reichl L. A Modern Course in Statistical Physics, 2nd ed., J. Wiley and Sons, New
York 1998.
[19] Reif F. Statistical Physics, Berkeley Physics Course, V5, McGraw-Hill, New York
1965.
[22] Stanley H.E. Introduction to Phase Transitions and Critical Phenomena, Oxford Uni-
versity Press, Oxford 1971.
302 BIBLIOGRAPHY