Catatan Kuliah Fisika Statistik 201008xx

Unduh sebagai pdf atau txt
Unduh sebagai pdf atau txt
Anda di halaman 1dari 194

Catatan Kuliah Fisika Statistik

Sparisoma Viridi, Siti Nurul Khotimah, dan Novitrian

Agustus 2010
ii
Isi

1 Pendahuluan 1

1.1 Ruang lingkup . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1

1.2 FI3202 Fisika Statistik (4 SKS) . . . . . . . . . . . . . . . . . . . 2

1.3 Versi catatan kuliah . . . . . . . . . . . . . . . . . . . . . . . . . 3

1.4 Buku rujukan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

2 Faktorial dan Fungsi Gamma 5

2.1 Faktorial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

2.2 Fungsi gamma untuk n bulat positif . . . . . . . . . . . . . . . . 5


1
2.3 Fungsi gamma untuk n kelipatan ganjil 2 . . . . . . . . . . . . . 7

2.4 Fungsi gamma yang lebih umum . . . . . . . . . . . . . . . . . . 9

2.5 Aproksimasi Striling . . . . . . . . . . . . . . . . . . . . . . . . . 9

2.6 Aproksimasi dengan grafik . . . . . . . . . . . . . . . . . . . . . . 9

2.7 Aproksimasi lain . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

2.8 Script menggambar grafik ln n . . . . . . . . . . . . . . . . . . . . 10

2.9 Referensi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11

3 Pengali Tak Tentu Lagrange 13

3.1 Maksimum dan minimum suatu fungsi . . . . . . . . . . . . . . . 13

3.2 Syarat tambahan . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

3.3 Referensi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

iii
iv ISI

4 Konfigurasi Paling Mungkin Suatu Statistik 19

4.1 Syarat batas suatu sistem . . . . . . . . . . . . . . . . . . . . . . 19

4.2 Memaksimumkan W . . . . . . . . . . . . . . . . . . . . . . . . . 21

4.3 Distribusi suatu statistik . . . . . . . . . . . . . . . . . . . . . . . 23

4.4 Referensi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24

5 Parameter 25

5.1 Dua buah sistem kontak secara termal . . . . . . . . . . . . . . . 25

5.2 Hukum pertama termodinamika . . . . . . . . . . . . . . . . . . 27

5.3 Teori kinetik gas . . . . . . . . . . . . . . . . . . . . . . . . . . . 29

5.4 Referensi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

6 Degenerasi dalam Ruang Fasa 31

6.1 Ruang fasa enam dimensi . . . . . . . . . . . . . . . . . . . . . . 31

6.2 Integral volume ruang momentum . . . . . . . . . . . . . . . . . 32

6.3 Integral volume ruang laju . . . . . . . . . . . . . . . . . . . . . . 33

6.4 Integral volume ruang energi . . . . . . . . . . . . . . . . . . . . 33

6.5 Integral volume ruang frekuensi . . . . . . . . . . . . . . . . . . . 34

6.6 Integral volume ruang panjang gelombang . . . . . . . . . . . . . 34

6.7 Referensi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35

7 Distribusi Suatu Statistik 37

7.1 Bentuk umum distribusi ketiga statistik . . . . . . . . . . . . . . 37

7.2 Statistik Maxwell-Boltzmann . . . . . . . . . . . . . . . . . . . . 38

7.3 Statistik Bose-Einstein . . . . . . . . . . . . . . . . . . . . . . . . 40

7.4 Statisti Fermi-Dirac . . . . . . . . . . . . . . . . . . . . . . . . . 41

7.5 Bentuk umum distribusi statistik lain . . . . . . . . . . . . . . . 42

7.6 Referensi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
ISI v

8 Termodinamika Gas Ideal Monoatomik 43

8.1 Peluang termodinamika Wmaks gas ideal klasik . . . . . . . . . . 43

8.2 Fungsi partisi Boltzmann . . . . . . . . . . . . . . . . . . . . . . 45

8.3 Tekanan dan kalor jenis . . . . . . . . . . . . . . . . . . . . . . . 47

8.4 Persamaan keadaan . . . . . . . . . . . . . . . . . . . . . . . . . 49

8.5 Referensi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50

9 Paradoks Gibb 51

9.1 Entropi gas klasik . . . . . . . . . . . . . . . . . . . . . . . . . . 51

9.2 Pencampuran dua gas berbeda jenis . . . . . . . . . . . . . . . . 52

9.3 Pencampuran gas sejenis: paradoks Gibb . . . . . . . . . . . . . 54

9.4 Gas ideal semi-klasik . . . . . . . . . . . . . . . . . . . . . . . . . 55

9.5 Referensi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56

10 Statistik Fermi-Dirac: N j dan S 57

10.1 Soal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58

10.2 Jawab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60

11 Tingkat dan Keadaan Energi 63

11.1 Tingkat Energi . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63

11.2 Keadaan Energi . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63

12 Keadaan Makro dan Mikro 65

13 Peluang Termodinamika 67

13.1 Postulat termodinamika statistik . . . . . . . . . . . . . . . . . . 67

13.2 Peluang termodinamika . . . . . . . . . . . . . . . . . . . . . . . 68

13.3 Observabel dan rata-rata bilangan okupasi . . . . . . . . . . . . . 68

14 Pengali dan 71
vi ISI

14.1 Peluang termodinamik suatu keadaan makro . . . . . . . . . . . 71

14.2 Keadaan makro yang paling mungkin . . . . . . . . . . . . . . . . 72

14.3 Fungsi distribusi dalam bentuk diferensial . . . . . . . . . . . . . 72

14.4 Pengali . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73

14.5 Ruang fasa enam dimensi . . . . . . . . . . . . . . . . . . . . . . 76

14.6 Degenerasi dalam volume ruang fasa . . . . . . . . . . . . . . . . 77

14.7 Teori kinetik gas dan . . . . . . . . . . . . . . . . . . . . . . . 78

14.8 Pengali . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79

15 Energi Bebas Helmholtz 83

15.1 Energy bebas Helmholtz . . . . . . . . . . . . . . . . . . . . . . . 83

15.2 Ekspansi reversibel . . . . . . . . . . . . . . . . . . . . . . . . . . 84

15.3 Energi sebagai fungsi dari energi bebas . . . . . . . . . . . . . . . 85

15.4 CV dari E . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85

16 Fungsi Partisi Boltzmann 87

16.1 Fungsi partisi Boltzmann . . . . . . . . . . . . . . . . . . . . . . 87

16.2 Fungsi partisi dan energi bebas Helmholtz . . . . . . . . . . . . . 88

16.3 Energi sistem dan fungsi partisi . . . . . . . . . . . . . . . . . . . 89

16.4 Entropi dan fungsi partisi . . . . . . . . . . . . . . . . . . . . . . 89

16.5 Energi bebas tiap partikel . . . . . . . . . . . . . . . . . . . . . . 90

16.6 Kapasitas panas spesifik pada volume tetap . . . . . . . . . . . . 90

16.7 Tekanan sistem . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90

17 Gas Ideal Monoatomik 91

17.1 Tingkat energi makro . . . . . . . . . . . . . . . . . . . . . . . . 91

17.2 Degenerasi tingkat energi makro . . . . . . . . . . . . . . . . . . 92

17.3 Fungsi partisi sistem . . . . . . . . . . . . . . . . . . . . . . . . . 93

17.4 Persamaan keadaan dan besaran-besaran termodinamika . . . . . 94


ISI vii

18 Distribusi Laju Molekular 95

18.1 Bilangan okupasi rata-rata tingkat energi makro . . . . . . . . . 95

18.2 Laju yang paling mungkin . . . . . . . . . . . . . . . . . . . . . . 96

18.3 Laju rata-rata . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97

18.4 Laju rms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97

18.5 Perbandingan ketiga laju . . . . . . . . . . . . . . . . . . . . . . 97

18.6 Ekipartisi energi . . . . . . . . . . . . . . . . . . . . . . . . . . . 98

19 Paradoks Gibb 99

19.1 Fungsi partisi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99

19.2 Beberapa besaran termodinamika . . . . . . . . . . . . . . . . . . 101

19.3 Paradoks Gibb . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103

19.4 Gas ideal semi-klasik . . . . . . . . . . . . . . . . . . . . . . . . . 104

20 Ekipartisi Energi 107

20.1 Bentuk-bentuk energi . . . . . . . . . . . . . . . . . . . . . . . . 107

20.2 Rata-rata energi kinetik . . . . . . . . . . . . . . . . . . . . . . . 108

20.3 Rata-rata energi potensial mirip pegas . . . . . . . . . . . . . . . 109

20.4 Rata-rata energi osilator harmonik . . . . . . . . . . . . . . . . . 110

20.5 Derajat kebebasan . . . . . . . . . . . . . . . . . . . . . . . . . . 111

20.6 Gas diatomik . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111

20.7 Bukan suku kuadrat koordinat . . . . . . . . . . . . . . . . . . . 111

21 Tambahan Informasi 1 113

21.1 Ilustrasi Cv bergantung T . . . . . . . . . . . . . . . . . . . . . . 113

21.2 Publikasi mengenai gas ideal dan ensembel mikrokanonik . . . . 114

22 Gas Ideal dalam Medan Gravitasi 117

22.1 Sistem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 117


viii ISI

22.2 Persamaan termodinamika . . . . . . . . . . . . . . . . . . . . . . 117

22.3 Energi total . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118

22.4 Fungsi partisi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118

22.5 Energi bebas Helmholtz . . . . . . . . . . . . . . . . . . . . . . . 120

22.6 Entropi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 121

22.7 Distribusi partikel sebagai fungsi ketinggian . . . . . . . . . . . . 121

22.8 Percobaan Perrin . . . . . . . . . . . . . . . . . . . . . . . . . . . 124

23 Gas diatomik 125

23.1 Suku-suku energi . . . . . . . . . . . . . . . . . . . . . . . . . . . 125

23.2 Fungsi-fungsi partisi . . . . . . . . . . . . . . . . . . . . . . . . . 126

23.3 Fungsi partisi gerak translasi . . . . . . . . . . . . . . . . . . . . 126

23.4 Fungsi partisi gerak rotasi . . . . . . . . . . . . . . . . . . . . . . 126

23.5 Fungsi partisi gerak vibrasi . . . . . . . . . . . . . . . . . . . . . 127

23.6 Fungsi partisi gerak elektron . . . . . . . . . . . . . . . . . . . . 127

23.7 Fungsi partisi spin nuklir . . . . . . . . . . . . . . . . . . . . . . 128

23.8 Fungsi partisi lengkap . . . . . . . . . . . . . . . . . . . . . . . . 128

23.9 Panas spesifik gas . . . . . . . . . . . . . . . . . . . . . . . . . . . 129

24 Gas Bose-Einstein 133

24.1 Distribusi molekul gas . . . . . . . . . . . . . . . . . . . . . . . . 133

24.2 Gas foton dan radiasi benda hitam . . . . . . . . . . . . . . . . . 135

24.3 Hukum radiasi Wien . . . . . . . . . . . . . . . . . . . . . . . . . 137

24.4 Formula Rayleigh-Jeans . . . . . . . . . . . . . . . . . . . . . . . 138

24.5 Hukum Stefan-Boltzmann . . . . . . . . . . . . . . . . . . . . . . 138

25 Gas Fermi-Dirac 141

25.1 Distribusi partikel . . . . . . . . . . . . . . . . . . . . . . . . . . 141

25.2 Fungsi Fermi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 142


ISI ix

26 Ensemble Kanonis 143

26.1 Ensemble . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 143

26.2 Ensemble yang bertemperatur konstan . . . . . . . . . . . . . . . 144

26.3 Sifat-sifat termodinamik ensemble kanonis . . . . . . . . . . . . . 145

26.4 Evaluasi Fungsi Partisi Total . . . . . . . . . . . . . . . . . . . . 147

26.5 Fungsi Partisi Klasik . . . . . . . . . . . . . . . . . . . . . . . . . 147

26.6 Fungsi partisi semi-klasik . . . . . . . . . . . . . . . . . . . . . . 148

26.7 Fungsi Partisi untuk Kasus Ada Interaksi . . . . . . . . . . . . . 149

26.8 Distribusi energi pada ensembel kanonik . . . . . . . . . . . . . . 149

26.9 Aplikasi ensemble kanonis untuk gas tak ideal . . . . . . . . . . . 150

27 Simulasi: Sistem Paramagnetik 153

28 Soal 1: Tingkat Energi dan Peluang Termodinamika 155

28.1 Soal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 155

28.2 Jawab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 156

29 Soal 2: Fungsi Distribusi dan Entropi 161

29.1 Soal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 161

29.2 Jawab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 162

30 Soal 3: Fungsi Partisi dan Tabulasi Keadaan Makro 165

30.1 Soal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 165

30.2 Jawab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 165

31 Soal 4: Distribusi Laju dan Persamaan Keadaaan 169

31.1 Soal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 169

31.2 Jawab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 170

32 Simulasi Keadaan Mikro dengan Kartu 175


x ISI

33 Berkas-berkas 181

33.1 Kuis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 181

33.2 Ujian . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 183


Catatan 1

Pendahuluan

Setiap cabang khusus fisika mula-mula dipelajari dengan memisahkan ruang


yang terbatas dari lingkungannnya. Bagian yang dipisahkan yang menjadi
pusat perhatian kita disebut sistem, dan segala sesuatu diluar sistem dise-
but lingkungan. Bila suatu sistem telah dipilih maka kelakuan sistem atau
antaraksinya dengan lingkungan atau keduanya dinyatakan dalam kuantitas-
kuantitas fisis. Pada umumnya terdapat dua pandangan yang dapat diambil,
pandangan makroskopik dan pandangan mikroskopik.

1.1 Ruang lingkup

Pemerian makroskopik suatu sistem meliputi perincian beberapa sifat pokok


sistem, atau sifat skala besar dari sistem, yang dapat diukur berdasarkan atas
penerimaan indera kita. Termodinamika adalah contoh cabang ilmu fisika yang
menerapkan pandangan makroskopik. Sedangkan, pemerian mikroskopik suatu
sistem meliputi beberapa ciri khas seperti adanya pengandaian bahwa sistem
terdiri atas sejumlah molekul, dan kuantitas-kuantitas yang diperinci tidak da-
pat diukur. Contoh penerapan pandangan mikroskopik untuk cabang ilmu fisika
yaitu dalam fisika statistik. Bila kedua pandangan itu diterapkan pada sistem
yang sama maka keduanya harus meghasilkan kesimpulan yang sama.

Ruang lingkup fisika statistik meliputi dua bagian besar, yaitu teori kinetik dan
mekanika statistik. Berdasarkan pada teori peluang dan hukum mekanika, teori
kinetik mampu menggambarkan sistem dalam keadaan tak seimbang, seperti:
proses efusi, viskositas, konduktivitas termal, dan difusi. Disini, molekul suatu
gas ideal tidak dianggap bebas sempurna tetapi ada antaraksi ketika bertum-
bukan dengan molekul lain atau dengan dinding. Bentuk antaraksi yang ter-
batas ini diacukan sebagai antaraksi lemah atau kuasi bebas. Ruang lingkup ini
tidak membahas partikel berantaraksi kuat

Tidak seperti pada teori kinetik, mekanika statistik tidak membahas perin-

1
2 CATATAN 1. PENDAHULUAN

cian mekanis gerak molekular, tetapi berurusan dengan segi energi molekul.
Mekanika statistik sangat mengandalkan teori peluang untuk menentukan
keadaan seimbang sistem. Dalam kuliah ini, bahasan ditekankan pada sistem
yang partikel-partikelnya berinteraksi sangat lemah baik untuk partikel-partikel
terbedakan maupun tak terbedakan. Selain memiliki sifat kuasi bebas, molekul-
molekul suatu gas ideal bersifat tak terbedakan karena molekul tidak berkecen-
derungan menempati tempat tertentu dalam ruang atau memiliki kecepatan
tertentu. Sedangkan, untuk partikel-partikel yang menempati kedudukan kisi
yang teratur dalam kristal, yakni partikel bergetar di sekitar titik tetap, dapat
dibedakan karena letaknya.

Materi kuliah mencakup probabilitas dan fungsi distribusi, teori kinetik, dan
mekanika statistik. Selain itu juga disentuh pengertian ensemble, terutama
ensemble kanonis untuk perluasan penerapan pada gas yang menyimpang dari
sifat ideal.

1.2 FI3202 Fisika Statistik (4 SKS)

Kuliah ini bertujuan untuk meletakkan dasar fisika statistik kepada mahasiswa
tingkat 3 jenjang stratum 1. Setelah mengikuti kuliah ini, mahasiswa diharap-
kan: (1) memahami peran dan kedudukan fisika statistik dalam bidang fisika,
(2) memahami dasar-dasar fisika statistik, (3) dapat menerapkannya dalam
masalah sederhana, dan (4) dapat memahami kuliah lanjut tentang sifat-sifat
zat maupun kuliah lain yang menggunakan fisika statistik.

Isi kuliah meliputi:

Probabilitas dan fungsi distribusi,

Teori kinetik gas: anggapan dasar, fluks molekul, tekanan, persamaan


keadaan, prinsip ekipartisi energi,

Fungsi distribusi laju menurut Maxwell,

Gejala transport: penampang tumbukan, jalan bebas rata-rata, viskositas


gas, konduktivitas termal gas, difusi gas.

Mekanika statistik: tingkat energi, keadaan energi, keadaan makro,


keadaan mikro,

Statistik Maxwell-Boltzmann: peluang termodinamik, penurunan dis-


tribusi partikel, fungsi partisi, entropi dan paradoks Gibbs,

Statistik semi-klasik: entropi, fungsi Helmholtz,

Statistik Bose-Einstein: peluang termodinamik, penurunan distribusi par-


tikel,

Statistik Fermi-Dirac: peluang termodinamik, penurunan distribusi par-


tikel,
1.3. VERSI CATATAN KULIAH 3

Keterbatasan ansambel mikrokanonis,


Ansambel kanonis: gas riel dengan interaksi lemah.

Prasyarat: FI2182 Fisika Moderen, FI2102 Fisika Matematika IA, FI2202


Fisika Matematika II, dan FI2202 Termodinamika.

Keempat prasayarat tersersebut sebaiknya telah dipenuhi agar peserta mataku-


liah ini tidak mengalami kesulitan dalam mengikuti materi-materi dalam perku-
liahan ini.

1.3 Versi catatan kuliah

Terdapat tiga versi catatan kuliah in sebelumnya, yaitu versi Mei 2010 yang
digunakan dalam kuliah pada Semester II Tahun 2009/2010, versi Juli 2010
yang digunakan dalam kuliah Semester III Tahun 2009/2010, dan versi draft
yang merupakan gabungan versi Mei 2010 ditambahkan dengan contoh simulasi
untuk diajukan pada hibah penulisan buku. Versi yang ada sekarang adalah
versi Agustus 2010 yang merupakan gabungan kesemua versi di atas. Oleh
karena itu versi ini terlihat agak tidak terintegrasi.

1.4 Buku rujukan

Buku rujukan utama kuliah ini adalah

1. A. J. Pointon, An Introduction to Statistical Physics for Students, Long-


mans, First Print, (1967)
2. Francis W. Sears and Gerhard L. Salinger, Thermodynamics, Kinetic
Theory, and Statistical Thermodynamics, Addison-Wesley, Third Edi-
tion, Fifth Print, (1980)
4 CATATAN 1. PENDAHULUAN
Catatan 2

Faktorial dan Fungsi


Gamma

Fungsi gamma atau biasa dituliskan sebagai (n) dan kaitannya dengan fak-
torial n! akan dibicarakan dalam tulisan ini. Detil mengenai relasi tersebut
dapat dilihat dalam literatur [1]. Faktorial untuk bilangan bulat dan setengah
bulat akan digunakan dalam distribusi Maxwell-Boltzmann untuk energi, mo-
mentum, dan laju dalam suatu asembli klasik [2] dan juga dalam penurunan
fungsi distribusi partikel yang memenuhi berbagai jenis statistik [3].

2.1 Faktorial

Faktorial dari suatu bilangan bulat, misalnya saja n, memiliki arti

n! = n (n 1) (n 2) (n 3) 3 2 1, (2.1)

di mana

0! = 1. (2.2)

Demikinlah nilai faktorial terdefinisi pada nilai bilangan bulat positif.

2.2 Fungsi gamma untuk n bulat positif

Fungsi gamma didefinisikan sebagai

5
6 CATATAN 2. FAKTORIAL DAN FUNGSI GAMMA

Z
(n + 1) = ex xn dx. (2.3)
0

Dengan melakukan integrasi parsial terhadap Persamaan (2.3) dapat diperoleh

(n + 1) = n(n), (2.4)

yang apabila dituliskan lebih jauh

(n + 1) = n (n 1) (n 2) (n 3) 2 1 (1) = n! (1). (2.5)

Dengan menggunakan Persamaan (2.3) dapat dihitung bahwa

Z
(1) = ex dx = 1, (2.6)
0

sehingga dapat diperoleh hubungan antara fungsi gamma dan faktorial, yaitu

(n + 1) = n!. (2.7)

Soal 1. Aturan LH
opital menyatakan bahwa

f (x) f (x)
lim = lim , (2.8)
xa g(x) xa g (x)

di mana f (a) dan g(a) keduanya bernilai nol. Gunakan Persamaan (2.8) untuk
menghitung


ex xn

x=0
. (2.9)

Jawab 1. Persamaan (2.9) dapat dituliskan dalam bentuk

 x n  ex ex ex
e x x=0 = lim n lim n = lim n 0. (2.10)
x x x0 x x x

Persamaan (2.10) dengan menggunakan aturan LH


opital akan menjadi

ex ex ex
lim = lim = lim
x x n x nx n+1 x n(n 1)xn+2

ex ex
= lim = = lim = 0.
x n(n 1)(n 2)xn+3 x n!
1
2.3. FUNGSI GAMMA UNTUK N KELIPATAN GANJIL 2 7

Penggunaan aturan LH opital dalam persamaan sebelumny dilakukan pada ben-


tuk xn /ex , baru kemudian pada setiap langkah dievaluasi untuk bentuk 0/0-nya
(ex /xn ).

Soal 2. Buktikan hubungan dalam Persamaan (2.4) dengan menggunakan in-


tegral parsial pada Persamaan (2.3) dan hasil dari Persamaan (2.9).

Jawab 2. Intergral pada ruas kanan Persamaan (2.3) dihitung melalui interasi
parsial

Z
(n + 1) = ex xn dx
0
Z Z
x n
 x n 
e x dx = e x x=0 + n ex xn1 dx
0 0
Z Z
ex xn dx = 0 + n ex xn1 dx
0
Z Z0
x n
e x dx = n ex xn1 dx
0 0
(n + 1) = n(n),

yang memberikan sifat rekusif dari fungsi gamma seperti dituliskan dalam Per-
samaan (2.4).

Soal 3. Hitunglah 0! dengan menggunakan fungsi gamma.

Jawab 3. Dari Persamaan (2.7) dapat diperoleh bahwa 0! = (1) dan dari
Persamaan (2.6) diperoleh bahwa (1) = 1. Dengan demikian dapat diperoleh
bahwa 0! = 1.

1
2.3 Fungsi gamma untuk n kelipatan ganjil 2

Dengan menggunakan Persamaan (2.3) dapat dituliskan bahwa


1
Z
1
( ) = ex x 2 dx. (2.11)
2 0

Soal 4. Turunkan Persamaan (2.11) dari Persamaan (2.3).

Jawab 4. Gunakan nilai n = 12 .

Bila n adalah setengah bilangan bulat maka fungsi gamma akan memberikan
hubungan

5 3 1 1
(n + 1) = n (n 1) (n 2) (n 3) ( ). (2.12)
2 2 2 2
8 CATATAN 2. FAKTORIAL DAN FUNGSI GAMMA

Persamaan (2.12) tetap memenuhi hubungan dalam Persamaan (2.4).

Nilai dari ( 21 ) sendiri dapat dihitung dengan menyelesaikan Persamaan (2.11)


sehingga diperoleh bahwa


Z
1
ex x 2 dx = . (2.13)
0

Soal 5. Buktikan Persamaan (2.13).

Jawab 5. Pertama-tama tuliskan Persamaan (2.13) dalam bentuk

Z
1
eu u 2 du.
0

Lalu misalkan u = x2 sehingga

Z Z Z
2 2 2
ex x1 (2xdx) = 2 ex dx = ex dx.
0 0

Misalkan bahwa
Z
2
I Ix = ex dx.

sehingga

Z Z Z Z
2 x2 y 2 2
+y 2 )
I Ix Iy = e dx e dy = e(x dxdy.
x= y=

Ubahlah elemen luas dalam sistem koordinat kartesian (dx)(dy) menjadi elemen
luas dalam sistem koordinat polar (dr)(rd) dan dengan hubunga r2 = x2 + y 2
sehingga

2 2
1 2 1
Z Z Z Z
2 r 2 r 2
I e rdrd = e d(r ) d = 1 2 = .
r=0 =0 r=0 2 =0 2

Dengan demikian dapat diperoleh bahwa

1
Z
1
I= eu u 2 du. = ( ) = .
0 2

Soal 6. Hitunglah ( 52 ).
3
Jawab 6. ( 52 ) = (1 + 23 ) = 3
2 1
2 ( 12 ) = 4 .
2.4. FUNGSI GAMMA YANG LEBIH UMUM 9

2.4 Fungsi gamma yang lebih umum

Secara umum dapat dituliskan bahwa

r
1 1
Z
1
ex x 2 dx = 1 ( ) = (2.14)
0 2 2

dan


1
Z Z
2
xn eax dx = y (n1)/2 ey dy
0 2a(n+1)/2 0
1
= [(n + 1)/2]. (2.15)
2a(n+1)/2

2.5 Aproksimasi Striling

Aproksimasi Striling yang berguna untuk menyederhanakan faktorial dan saat


menurunkannya adalah


n! nn en 2n (2.16)

atau

1 1
ln n! (n + ) ln n n + ln(2). (2.17)
2 2

2.6 Aproksimasi dengan grafik

Aproksimasi lain untuk ln n! dapat diperoleh lewat grafik seperti ditunjukkan


dalam Gambar 2.1. Dengan demikian dapat dituliskan bahwa aproksimasi untuk
ln n! adalah

n n
1 1
X Z
ln n! = ln i ln n + ln xdx = (n + ) ln n n + 1. (2.18)
i=1
2 1 2

Soal 7. Buktikan dari grafik aproksimasi dalam Persamaan (2.18) dengan meng-
gunakan Gambar 2.1.

Jawab 7. Bahas luas dari kurva di bawah ln x untuk kotak pertam, di mana
kelebihan kotak sebelah kanan titik tengah n adalah untuk bagian sebelah kiri di
10 CATATAN 2. FAKTORIAL DAN FUNGSI GAMMA

bawah kurva. Demikian seterusnyaR sehingga tersisa saat kotak ke n ada faktor
1
2 ln n yang belum dihitung dalam ln xdx. Lebar tiap kotak adalah 1.

2.5

1.5
ln n

0.5

0
0 1 2 3 4 5 6 7 8 9 10 11 12
n

Gambar 2.1: Histogram dari ln n dan kurva ln x yang saling digambarkan


bertumpang tindih.

2.7 Aproksimasi lain

Dengan melihat nilai yang besar dari n di mana umumnya merupakan daerah
kerja mekanika statistik dan umumnya yang dibahas adalah perubahan nilai
atau turunan dari ln n! maka aproksimasi lain digunakan, yaitu

ln n! n ln n n. (2.19)

Soal 8. Bandingkan nilai-nilai ln n! dengan menggunakan Persamaan (2.17),


(2.18), dan (2.19).

Jawab 8. Pembandingan nilai-nila ln n! yang diminta dapat dilihat dalam Tabel


2.1.

2.8 Script menggambar grafik ln n

set term post eps color enhanced 28 lw 1


set output "ln-n.eps"
set size 1.4, 1.2
set xrange [-0.2:12.2]
2.9. REFERENSI 11

Tabel 2.1: Nilai-nilai ln n! dengan menggunakan Persamaan (2.17), (2.18), dan


(2.19).
Persamaan
n ln n!
(2.17) (2.18) (2.19)
10 15.1 15.1 15.18 13.03
50 148.48 148.48 148.56 145.6
100 363.74 363.74 363.82 360.52
150 605.02 605.02 605.1 601.6
170 706.57 706.57 706.65 703.09

set yrange [-0.1:2.6]


set xtics 1
set ytics 0.5
set grid
#set label "{/Italics D}_{/Italics L}" \
#at 0, 0.12
set xlabel "{/Italics n}"
set ylabel "ln {/Italics n}"
plot \
"data.txt" u 1:(log($1)) t "" lw 4 w boxes, \
log(x) t "" lw 4

Script di atas dipanggil dengna menggunakan aplikas gnuplot sehingga berkas


ln-n.eps dihasilkan seperti dalam Gambar 2.1.

2.9 Referensi
1. Mary L. Boas, Mathematical Methods in the Physical Sciences, John
Wiley & Sons, New York, Second Edition, 457-462 (1983)
2. A. J. Pointon, An Introduction to Statistical Physics for Students, Long-
mans, First Print, 189-191 (1967)
3. Francis W. Sears and Gerhard L. Salinger, Thermodynamics, Kinetic
Theory, and Statistical Thermodynamics, Addison-Wesley, Third Edi-
tion, Fifth Print, 424-426 (1980)
12 CATATAN 2. FAKTORIAL DAN FUNGSI GAMMA
Catatan 3

Pengali Tak Tentu Lagrange

Pengali Tak Tentu Lagrange (the Lagrange method of undetermined multipli-


ers) adalah suatu metoda matematika untuk mencari maksimum atau minimum
suatu fungsi yang dibatasi oleh suatu syarat dari fungsi lain [1], di mana metode
ini termasuk dalam optimisasi matematika [2]. Bila jumlah variabel bebas dan
syarat yang membatasi sedikit, cukup dilakukan substitusi standar. Akan tetapi
bila jumlah variabel bebas dan syarat-syarat banyak, maka akan terdapat terlalu
banyak fungsi yang harus diselesaikan secara bersama-sama. Di sinilah metoda
ini ini berperan dengan memperkenalkan satu konstanta untuk setiap syarat
dari satu fungsi lain yang diperlukan [3]. Penjelasan yang cukup sederhana
dapat dilihat dalam literatur [4].

3.1 Maksimum dan minimum suatu fungsi

Bila terdapat suatu fungsi f (x, y, z) yang ingin dicari nilai maksimum atau
minimumnya, maka cukup dipenuhi

f f f
df = dx + dy + dz = 0. (3.1)
x y z

Soal 1. Suatu lingkaran yang terletak di pusat koordinat dengan jari-jari R


memiliki fungsi f (x, y) = x2 + y 2 R2 = 0. Tentukanlah nilai maksimum dan
minimum dari x.

Jawab 1. Gunakan Persamaan (3.1) sehingga diperoleh

df = 2xdx + 2ydy = 0.

Untuk mencari nilai maksimum dari x maka perlu dicari nilai y lewat

13
14 CATATAN 3. PENGALI TAK TENTU LAGRANGE

dx y
= = 0 y = 0.
dy x

Gunakan nilai ini ke persamaan lingkaran sehinga diperoleh bahwa

x2 R2 = 0 x = R.

Jadi nilai maksimum dan minimum x berturut-turut adalah R dan +R.

Soal 2. Tentukan nilai y minimum dari fungsi f (x, y) = y x2 = 0.

Jawab 2. Gunakan Persamaan (3.1) sehingga diperoleh

df = dy 2xdx = 0.

Untuk mencari nilai minimum y perlu dicari nilai x dari

dy
= 2x = 0 x = 0.
dx

Dengan demikian dapat diperoleh bahwa

y 02 = 0 y = 0,

yang merupakan nilai minimum y.

3.2 Syarat tambahan

Mencari minimum atau maksimum suatu fungsi f (x, y, z) tidak cukup dengan
menggunakan Persamaan (3.1) bila terdapat syarat tambahan berupa fungsi
lain, misalnya (x, y, z). Untuk itu diperkenalan dengan suatu metode yang
menggunakan pengali berupa konstanta yang belum diketahui nilainya, pen-
gali tak tentu Lagrange, sehingga perluasan dari Persamaan (3.1) yang meru-
pakan kondisi yang harus terpenuhi adalah

df + d = 0, (3.2)

dengan bentuk d mirip dengan bentuk df . Bila bentuk (x, y, z) dapat dit-
uliskan dalam bentuk
3.2. SYARAT TAMBAHAN 15

(x, y, z) = 0 (3.3)

maka


d = dx + dy + dz = 0. (3.4)
x y z

Selanjutnya dapat diperoleh dari Persamaan (3.4) bentuk

   

dx = dy + dz / , (3.5)
y z x

sebagaimana untuk dy maupun dz. Substitusi Persamaan (3.5) ke Persamaan


(3.1) akan memberikan

f f f
df = dx + dy + dz = 0,
x y z
    
f f f
df = dy + dz / + dy + dz = 0,
x y z x y z
     
f f /y f f /z
df = dy + dz = 0. (3.6)
y x /x z x /x

f
Dengan x / x pada titik stasioner diberi nilai maka

 
f
+ = 0. (3.7)
x x

Persamaan (3.6) dapat dituliskan menjadi

     
f f f
df + dx + + dy + + dz = 0. (3.8)
x x y y z z

Agar Persamaan (3.8) dapat menentukan suatu titik stasioner maka setiap suku
dalam tanda kurung harus bernilai nol, sebagaimana persamaan tersebut harus
dipenuhi untuk setiap nilai dari perubahan dy dan dz, maka kurung kedua dan
ketiga harus bernilai nol, sementara kurung pertama bernilai nol akibat definisi
dari dalam Persamaan (3.7).

Karena merupakan rasio dari f


x / x pada suatu titik stasioner, nilaniya
hanya dapat ditentukan dengan melakukan substitusi kembali solusi yang diper-
oleh ke persamaan yang merupakan syarat awal (x, y, z) = 0. Oleh karena itu
16 CATATAN 3. PENGALI TAK TENTU LAGRANGE

dikenal sebagai suatu pengali tak tentu atau lebih umum, pengali tak tenu
Lagrange (a Lagrange undetermined multiplier).

Soal 3. Suatu fungsi f (x, y, z) ingin diminimumkan dengan syarat-syarat


1 (x, y, z), 2 (x, y, z), dan 3 (x, y, z). Tentukanlah bentuk persamaan yang
harus dipecahkan dengan memperkenalkan tiga buah pengali tak tentu La-
grange.

Jawab 3. Dengan menggunakan Persamaan (3.2) dapat dituliskan

df = 1 + 2 + 3 = 0

adalah fungsi yang harus dipecahkan dengan pengali-pengali tak tentu La-
grangenya adalah , , dan , yang akan dicari kemudian.

Soal 4. Tentukan kuadrat jarak minimum dan maksimum dari suatu titik (2R,
0) terhadap lingkaran x2 + y 2 = R2 .

Jawab 4. Fungsi yang ingin dicari titik stasionernya adalah

f (x, y) = (x 2R)2 + y 2

dengan syarat batas

(x, y) = x2 + y 2 R2 = 0.

Dengan menggunakan Persamaan (3.2) dapat dituliskan bahwa

df + d = [2(x 2R)dx + 2ydy] + [2xdx + 2ydy] = 0.


[x(1 + ) 2R]dx + y(1 + )dy = 0.
dx y(1 + )
= . (3.9)
dy 2R x(1 + )

Nilai x minimum dan maksimum dapat dicari dengan membuat Persamaan (3.9)
menjadi nol, sehingga


y = 0, =?
0 = y(1 + )
= 1, y =?

Bila dipilih y = 0 maka diperoleh dari syarat batas bahwa x = R sehingga f


bernilai R2 dan (3R)2 . Dengan menggunakan x = 2R/(1 + ) dan = 1 tidak
memberikan solusi karena tidak memenuhi fungsi yang (x, y) yang membatasi.
3.3. REFERENSI 17

Soal 5. Tentukan kuadrat jarak minimum dan maksimum dari suatu titik (R,
R) terhadap lingkaran x2 + y 2 = R2 .

Jawab 5. Fungsi yang ingin dicari titik stasionernya adalah

f (x, y) = (x R)2 + (y R)2

dengan syarat batas

(x, y) = x2 + y 2 R2 = 0.

Dengan menggunakan Persamaan (3.2) dapat dituliskan bahwa

df + d = [2(x R)dx + 2(y R)dy] + [2xdx + 2ydy] = 0.


[x(1 + ) R]dx + [y(1 + ) R]dy = 0. (3.10)


Dapat dipilih = 2 1 agar Persamaan (3.10) dapat bernilai nol dan fungsi
yang membatasi tetap terpenuhi.
Dengan pilihan ini diperoleh bahwa kuadrat
2
jarak
minimum adalah (3 2 2)R dan kuadrat jarak maksimum adalah (3 +
2 2)R2 .

3.3 Referensi
1. Mary L. Boas, Mathematical Methods in the Physical Sciences, John
Wiley & Sons, New York, Second Edition, 174-181 (1983)

2. Wikipedia contributors, Lagrange multipliers, Wikipedia, The Free En-


cyclopedia, 26 May 2010, 20:32 UTC, http://en.wikipedia.org/w/
index.php?title=Lagrange multipliers&oldid=364362085 [accessed 6
July 2010]
3. Francis W. Sears and Gerhard L. Salinger, Thermodynamics, Kinetic
Theory, and Statistical Thermodynamics, Addison-Wesley, Third Edi-
tion, Fifth Print, 421-423 (1980)
4. A. J. Pointon, An Introduction to Statistical Physics for Students, Long-
mans, First Print, 189-191 (1967)
18 CATATAN 3. PENGALI TAK TENTU LAGRANGE
Catatan 4

Konfigurasi Paling Mungkin


Suatu Statistik

Dari ketiga statistik yang dipelajari, yaitu statistik Maxwell-Boltzmann, statis-


tik Bose-Einstein, dan statistik Fermi-Dirac, dapat diperoleh rumusan mengenai
peluang termodinamika suatu keadaan makro k, yaitu Wk . Dengan mengang-
gap bahwa suatu sistem tersusun atas banyak partikel maka terdapat suatu
puncak yang cukup tajam dari Wk terhadap nilai-nilai lain di sekelilingnya,
yang disebut sebagai Wk,maks dan dicari dengan memaksimumkan Wk [1].

4.1 Syarat batas suatu sistem

Sistem yang dibahas di sini dibatasi pada sistem tertutup dan terisolasi. Istilah
tertutup dan terisolasi terkait dengan jumlah total partikel dalam sistem N dan
energi total sistem U , di mana energi pada tingkat energi j adalah j dan jumlah
partikel yang menempati tingkat energi tersebut adalah Nj .

Soal 1. Bila terdapat M tingkat energi dengan masing-masing tingkat energi


ditempati oleh Nj partikel, tuliskan rumusan bagaimana menghitung jumlah
total partikel dalam sistem.

Jawab 1. Jumlah total partikel dalam sistem dihitung melalui

M
X
N= Nj . (4.1)
j=1

Soal 2. Energi tingkat energi j adalah j dan ditempati oleh Nj partikel.


Hitunglah energi total sistem U apabila terdapat M tingkat energi.

Jawab 2. Energi total sistem U dihitung melalui

19
20 CATATAN 4. KONFIGURASI PALING MUNGKIN SUATU STATISTIK

M
X
U= j Nj . (4.2)
j=1

Soal 3. Hitunglah energi rata-rata sistem bila energi tingkat energi j adalah j
dan ditempati oleh Nj partikel.

Jawab 3. Energi rata-rata sistem dihitung melalui

PM
U j=1 j Nj
= = PM . (4.3)
N j=1 Nj

Soal 4. Apa yang dimaksud dengan sistem tertutup dan terisolasi? Bagaimana
merumuskannya terkait dengan Persamaan (4.1) dan (4.2)?

Jawab 4. Sistem tertutup berarti bahwa jumlah partikel dalam sistem tetap.
Tidak terjadi perubahan jumlah partikel, jumlah partikel tidak berkurang
melalui keluarnya partikel dari sistem atau jumlah partikel tidak bertambah
melalui masuknya partikel ke dalam sistem. Syarat ini dirumuskan dengan


X X
dN = d Nj = dNj = 0. (4.4)
j j

Sedangkan sistem terisolasi berarti energi total sistem tetap yang dirumuskan
melalui


X X
dU = d j Nj = j dNj = 0. (4.5)
j j

Soal 5. Bagaimana cara mencari Wk,maks dari suatu sistem tertutup dan ter-
isolasi dengan memperkenalkan dua pengali tak tentu Lagrange a dan b?

Jawab 5. Fungsi yang harus dimaksimmukan adalah Wk dengan mencari tu-


runan parsialnya terhadap Nj dan syarat batas yang harus dipergunakan adalah
dN = 0 (sistem tertutup) dan dU = 0 (sistem terisolasi). Dengan demikian

dW + adN + bdU = 0, (4.6)

yang lebih eksplisitnya adalah

X Wk X (adN ) X (bdU )
dNj + dNj + dNj = 0,
j
Nj j
Nj j
Nj
4.2. MEMAKSIMUMKAN W 21

X Wk X (a P dNi ) X (b P i dNi )
i i
dNj + dNj + dNj = 0,
j
Nj j
N j j
N j

X Wk X Ni X Ni
dNj + aij dNj + bi ij dNj = 0,
j
Nj j
Nj j
N j
X Wk X X
dNj + a dNj + b j dNj = 0. (4.7)
j
Nj j j

Umumnya, yang dimaksimumkan bukanlah Wk akan tetapi ln Wk sehingga Per-


samaan (4.6) akan menjadi

dln W + dN + dU = 0, (4.8)

dengan memperkenalkan dan sebagai pengali tak tentu Lagrange. Dengan


menggunakan prosedur yang sama untuk menghasilkan Persamaan (4.7) dapat
diperoleh

X ln Wk X X
dNj + dNj + j dNj = 0. (4.9)
j
Nj j j

Selanjutnya Wk dalam Persamaan (4.9) akan dituliskan hanya sebagai W agar


tidak indeks k tidak membingungkan.

4.2 Memaksimumkan W

Telah diperoleh bahwa persamaan yang harus dipecahkan adalah Persamaan


(4.9), yang dapat dituliskan kembali menjadi

X  ln W 
+ + j dNj = 0. (4.10)
j
Nj

Untuk mencari nilai ln W maksimum (dapat juga minimum) maka harus pula
berlaku maksimum untuk setiap suku yang terkait dengan dNj , yang berarti
bahwa

ln W
+ + j = 0. (4.11)
Nj
Q
Soal 6. Bila W = j Nj !, selesaikan Persamaan (4.11) untuk setiap Nj dengan
menggunakan aproksimasi Stirling.

Jawab 6. Aproksimasi Strirling untuk ln n! dalam Persamaan (2.17) mem-


berikan ln N ! N ln N N . Dengan demikian dapat diperoleh bahwa
22 CATATAN 4. KONFIGURASI PALING MUNGKIN SUATU STATISTIK

X X
ln W = (ln Nj !) (Nj ln Nj Nj ).
j j

Soal 7. Tentukan ln W untuk statistik Maxwell-Boltzmann.

Jawab 7. Statistik Maxwell-Boltzmann

Y gjNj
WMB = N ! (4.12)
j
Nj !

sehingga

X
ln WMB = ln N ! + (Nj ln gj ln Nj !)
j
X
N ln N N + (Nj ln gj Nj ln Nj + Nj ) (4.13)
j

Soal 8. Tentukan ln W untuk statistik Bose-Einstein.

Jawab 8. Statistic Bose-Einstein

Y [(gj 1) + Nj ]!
WBE = (4.14)
j
(gj 1)!Nj !

sehingga

X X X
ln WBE = [(gj 1) + Nj ]! (gj 1)! Nj !
j j j
X
{[(gj 1) + Nj ] ln[(gj 1) + Nj ] [(gj 1) + Nj ]}
j
X X
{(gj 1) ln(gj 1) (gj 1)} {Nj ln Nj Nj }
j j
X
{[(gj 1) + Nj ] ln[(gj 1) + Nj ] (gj 1) ln(gj 1)
j
Nj ln Nj }. (4.15)

Soal 9 Tentukan ln W untuk statistik Fermi-Dirac.

Jawab 9. Statistik Fermi-Dirac


4.3. DISTRIBUSI SUATU STATISTIK 23

Y gj !
WFD = (4.16)
j
(gj Nj )!Nj !

sehingga

X X X
ln WFD = gj ! (gj Nj )! Nj !
j j j
X X
(gj ln gj gj ) [(gj Nj ) ln(gj Nj ) (gj Nj )]
j j
X
(Nj ln Nj Nj )
j
X
[gj ln gj (gj Nj ) ln(gj Nj ) Nj ln Nj ] (4.17)
j

4.3 Distribusi suatu statistik

Dengan menggunakan Persamaan (4.11) untuk masing-masing stastik seperti


dalam Persamaan (4.13), (4.15), dan (4.17), akan dapat diperoleh untuk kon-
figurasi yang paling mungkin distribusi partikel untuk masing-masing statistik,
yaitu

gj
Nj,MB = , (4.18)
e(+j )

gj
Nj,BE = , (4.19)
e(+j ) 1

gj
Nj,FD = . (4.20)
e (+ j) +1

Soal 10. Turunkan Persamaan (4.18).

Jawab 10. Dari Persamaan (4.13) dan Persamaan (4.11) dapat dituliskan dan
diperleh bahwa


X
ln N ln N + (Nj ln gj Nj ln Nj ) + + j = 0
Nj j

gj gj
ln + + j = 0 ln = ( + j )
Nj Nj
gj gj
= e(+j ) Nj = (+ ) ,
Nj e j
24 CATATAN 4. KONFIGURASI PALING MUNGKIN SUATU STATISTIK

seperti dalam Persamaan (4.18).

Soal 11. Turunkan Persamaan (4.19).

Jawab 11. Dari Persamaan (4.15) dan Persamaan (4.11) dapat dituliskan dan
diperleh bahwa


{[(gj 1) + Nj ] ln[(gj 1) + Nj ] (gj 1) ln(gj 1) Nj ln Nj }
Nj
+ + j = 0
ln[(gj 1) + Nj ] ln Nj + + j = 0
[(gj 1) + Nj ] gj 1
ln = ( + j ) = e(+j ) 1
Nj Nj
gj gj
gj >> 1 = e(+j ) 1 Nj = (+ ) ,
Nj e j 1

seperti dalam Persamaan (4.19).

Soal 12. Turunkan Persamaan (4.20).

Jawab 12. Dari Persamaan (4.17) dan Persamaan (4.11) dapat dituliskan dan
diperleh bahwa


{gj ln gj (gj Nj ) ln(gj Nj ) Nj ln Nj } + + j = 0
Nj
ln(gj Nj ) ln Nj + + j = 0
gj Nj gj
ln = ( + j ) = e(+j ) + 1
Nj Nj
gj
Nj = (+ ) ,
e j + 1

seperti dalam Persamaan (4.20).

4.4 Referensi
1. A. J. Pointon, An Introduction to Statistical Physics for Students, Long-
mans, First Print, 14-15 (1967)
Catatan 5

Parameter

Parameter yang digunakan sebagai salah satu pengali tak tentu Lagrange
untuk mencari nilai maksimum dari logaritma peluang termodinamika suatu
keadaan makro ln W, sebagaimana dituliskan dalam Persamaan (4.8), perlu
dicari artinya secara fisis. Distribusi partikel dari konfigurasi yang paling
mungkin untuk ketiga statistik, Maxwell-Boltzmann (MB), Bose-Einstein (BE),
dan Fermi-Dirac (FD), telah diperoleh dan masing-masing mengandung param-
eter sebagaimana dituangkan dalam Persamaan (4.18), (4.19), dan (4.20).
Bagaimana fungsi dari parameter dan bentuk eksplisitnya dapat dilihat pen-
jelasannya dalam [1] dan saduran bebasnya dalam [2].

5.1 Dua buah sistem kontak secara termal

Salah satu pendekatan yang digunakan untuk menunjukkan bagaimana intepre-


tasi secara fisis dari parameter adalah dengan memisalkan terdapatnya dua
buah sistem tertutup yang hanya dapat saling mempertukarkan energi, akan
tetapi tidak dapat mempertukarkan partikel. Kedua sistem yang dimaksud,
secara gabungan, dianggap sebagai sistem yang terisolasi.

Soal 1. Bagaimanakah rumusan dua buah sistem yang masing-masing tertutup


di mana keduanya dapat saling mempertukarkan energi, akan tetapi gabungan
keduanya merupakan suatu sistem terisolasi terhadap lingkungannya? Gunakan
untuk sistem pertama tanda dan untuk sistem kedua tanda .

Jawab 1. Kedua sistem merupakan sistem tertutup, sehingga dapat dituliskan


bahwa

dN = 0, (5.1)
dN = 0, (5.2)

25
26 CATATAN 5. PARAMETER

dan karena gabungan keduanya merupakan sistem yang terisolasi dengan


lingkungannya maka

dU = 0, (5.3)
dU = dU + dU . (5.4)

Saat dua buah sistem digabungkan maka ada parameter dalam sistem gabungan
yang merupakan hasil perkalian dari parameter masing-masing sistem. Salah
satu contoh parameter yang bersifat seperti ini adalah peluang termodinamika
suatu keadaan makro (yang mulai sekarang diambil tak lain adalah keadaan
makro yang paling mungkin muncul) W. Jadi bila peluang keadaan makro
yang paling mungkin muncul dari sistem pertama adalah W dan untuk sistem
kedua adalah W maka peluang keadaan makro sistem gabungan adalah

W = W W . (5.5)

Soal 2. Peluang termodinamika keadaan-keadaan makro suatu sistem adalah


20, 30, 4000, 35, 20, 5. Sedangan suatu sistem lain memiliki peluang termodi-
namika keadaan-keadaan makro 1, 5, 1500, 3, 1. Tentukanlah peluang keadaan
makro yang paling mungkin mumcul dari gabungan kedua sistem tersebut.

Jawab 2. Untuk sistem pertama W = 4000 dan untuk sistem kedua W =


1500, sehingga dengan menggunakan Persamaan (5.5) dapat diperoleh bahwa
W = 6000000.

Soal 3. Rumuskan dengan menggunakan pengali tak tentu Lagrange , ,


dan dua buah sistem tertutup yang dapat kontak secara termal dan meru-
pakan sistem gabungan yang terisolasi terhadap lingkungannya. Serta jelaskan
mengapa hanya perlu satu parameter .

Jawab 3. Dengan menggunakan W sistem gabungan dan syarat bahwa dN =


0, dN = 0, dan dU = 0 maka dapat dituliskan bahwa

d ln W + dN + dN + dU = 0, (5.6)

sehingga dapat diperoleh untuk tiap dNj dan dNj

W
+ + j = 0, (5.7)
Nj
W
+ + j = 0, (5.8)
Nj

karena W hanya bergantung dari Nj dan W hanya bergantung dari Nj .


5.2. HUKUM PERTAMA TERMODINAMIKA 27

Soal 4. Turunkan Persamaan (5.7) dan (5.8) dari Persamaan (5.6) dengan
menggunakan Persamaan (5.4).

Jawab 4. Dapat dituliskan dan diperoleh bahwa

d ln W + dN + dN + dU = 0,
d ln(W W ) + dN + dN + d(U + U ) = 0,
d ln W + d ln W + dN + dN + dU + dU = 0,
(d ln W + dN + dU ) + (d ln W + dN + dU ) = 0,
! !
X W X W

+ + j dNj + + + j dNj = 0
j
N j j
N j

seperti dalam Persamaan (5.7) dan (5.8) di mana masing-masing suku harus nol
untuk setiap perubahan dNj dan dNj .

Karena kedua sistem hanya dapat mempetukarkan kalor maka saat terjadinya
kesetimbangan hanya satu parameter yang akan berniali sama yaitu temperatur
T (hal ini sesuai dengan hukum ke-nol termodinamika). Dari (5.7) dan (5.8)
dapat dilihat bahwa hanya satu parameter yang sama untuk kedua sistem yaitu
. Dengan demikian dapat diperoleh bahwa seharusnya

= (T ) (5.9)

yang merupakan suatu intepretasi secara fisis kelakuan dari .

5.2 Hukum pertama termodinamika

Terdapat pula sudut pandang lain untuk melihat arti dari pengali yang me-
manfaatkan hubungan yang diungkapkan oleh hukum pertama termodinamika,
yaitu

dU = dQ pdV. (5.10)

Dengan menggunakan Persamaan (4.5) dalam bentuk yang lebih umum di mana
mungkin terdapat perubahan dj maka dapat dituliskan bahwa


X X X
dU = d j Nj = Nj dj + j dNj . (5.11)
j j j

Perubahan volume akan mengubah tingkat-tingkat energi sebagaimana kasus


partikel dalam kotak sedangkan perubahan kalor akan membuat terjadinya pe-
28 CATATAN 5. PARAMETER

rubahan susunan partikel dalam tingkat-tingkat energi. Dengan demikian dapat


dituliskan bahwa

X
Nj dj = pdV, (5.12)
j
X
j dNj = dQ. (5.13)
j

Soal 5. Pada saat tercapainya kesetimbangan sehingga tidak lagi terjadi pe-
rubahan volume, turunkan bentuk parameter secara eksplisit dengan meng-
gunakan rumusan pengali tak tentu Lagrange dalam mencari peluang termodi-
namika suatu keadaan makro W yang paling mungkin dan rumusan entropi dari
Boltzmann serta hubungan antara perubahan entropi dengan perubahan kalor.
Sistem merupakan sistem tertutup.

Jawab 5. Rumusan pengali tak tentu Lagrange dalam mencari peluang ter-
modinamika suatu keadaan makro W yang paling mungkin memberikan

dW + dN + dU = 0,

di mana bila tidak terjadi perubahan volume maka melalui hukum pertama
termodinamika

dU = dQ,

sehingga dapat diperoleh

dW + dN + dQ = 0.

Dengan menerapkan syarat sistem tertutup, yaitu dN = 0 maka dapat diperoleh


bahwa

d ln W = dQ.

Selanjutnya dengan menggunakan ungkapan Boltzmann untuk entropi, yaitu

S = k ln W, (5.14)

dan hubungan dS = dQ/T dapat dituliskan bahwa


5.3. TEORI KINETIK GAS 29

d ln W = dQ
 
S
d = dQ
k
dS
= dQ
k
dS
= k
dQ
1
= k
T
1
= . (5.15)
kT

Persamaan (5.15) menggambarkan hubungan eksplisit antara dang T .

5.3 Teori kinetik gas

Khusus untuk statistik Maxwell Boltzmann, terdapat hubungan

Nj = gj ej +j

seperti telah ditunjukkan oleh Persamaan (4.18). Sedangka teori kinetik gas
menyatakan bahwa energi rata-rata tiap partikel gas monoatomik adalah

3
= kT. (5.16)
2

Jumlah total partikel dapat diperoleh lewat

X X Z
j +j
N= Nj = gj e [2(2m)3/2 1/2 d/h3 ]e+ (5.17)
j j 0

dan energi sistem

X X Z
U= j Nj = j gj ej +j [2(2m)3/2 1/2 d/h3 ]e+ (5.18)
j j 0

di mana

2(2m)3/2 1/2 d
gj . (5.19)
h3
30 CATATAN 5. PARAMETER

Kemudian dengan menggunakan relasi yang diperoleh dari integral parsial


3
Z Z
3/2 e d = 1/2 e d (5.20)
0 2 0

dan bahwa = U/N maka

3 3 1
= = kT = ,
2 2 kT

seperti dalam Persamaan (5.15).

Soal 6. Buktikan Persamaan (5.20).

Jawab 6. Dengan melihat bentuk persamaan yang dimaksud maka dapat dit-
uliskan

1 3 1/2
Z Z
3/2 e d = 3/2 e e d
2
Z   Z
1 3/2 3
3/2 e d = e 1/2 e d
0 0 2 0
Z Z
3/2 3
e d = 0 1/2 e d
0 2
Z Z 0
3/2 3
e d = 1/2 e d,
0 2 0

di mana telah digunakan suatu asumsi mengenai nilai , yaitu bahwa < 0.

Soal 7.

5.4 Referensi
1. A. J. Pointon, An Introduction to Statistical Physics for Students, Long-
mans, First Print, 19-25 (1967)
2. Sparisoma Viridi dan Siti Nurul Khotimah, Catatan Kuliah Fisika Statis-
tik, Semester II Tahun 2009/2010, Mei, 21-27 (2010)
Catatan 6

Degenerasi dalam Ruang


Fasa

Degenerasi atau jumlah keadaan energi gj pada suatu tingkat energi j yang
memiliki energi antara j dan j + dj yang bersifat dikrit dapat dilihat menjadi
suatu besaran yang berharga kontinu [1]. Bagaimana hal itu dapat dilakukan,
akan diilustrasikan dalam catatan ini.

6.1 Ruang fasa enam dimensi

Saat sebuah partikel bergerak dalam ruang tiga dimensi (x, y, z) dan memiliki
momentum pada ketiga arah tersebut (px , py , pz ), keadaan partikel tersebut
setiap saat secara lengkap dispesifikasikan dengan enam koordinat yaitu (x, y,
z, px , py , pz ). Ruang di mana partikel dispesifikasikan dengan enam koordinat
tersebut disebut sebagai ruang enam dimensi atau ruang .

Soal 1. Bila elemen volume ruang koordinat tiga dimensi adalah dxdydz, ten-
tukanlah elemen volume ruang fasa enam dimensi .

Jawab 1. Ruang memiliki koordinat x, y, z, px , py , pz untuk tiap-tiap


partikel. Dengan demikian elemen volumenya adalah

d = (dV )(dVp ) = (dx, dy, dz)(dpx , dpy , dpz ) = dxdydzdpx dpy dpz . (6.1)

Kaitan antara gj dan d adalah

d
gj , (6.2)
h3

31
32 CATATAN 6. DEGENERASI DALAM RUANG FASA

di mana h adalah konstanta Planck, h = 6.626 1034 m2 kg s1 .

Bila fungsi yang akan diinteralkan, dalam hal ini adalah suku

1
, c = 1, 0, 1,
e+j +c

tidak bergantung pada koordinat spasial (x, y, z) maka d dapat dituliskan


menjadi

d = V dpx dpy dpz

yang artinya telah dilakukan integrasi terhadap elemen volume spasial.


Demikian pula bila suku tersebut tidak mengandung koordinat momentum (px ,
py , pz ) maka dapat dituliskan menjadi

d = Vp dxdydz

yang artinya telah dilakukan integrasi terhadap elemen volume momentum.

6.2 Integral volume ruang momentum

Elemen ruang momentum dpx dpy dpz dapat pula dituliskan sebagai

dVp = dpx dpy dpz = 4p2 dp

apabila sifat momentumnya dianggap isotropik, homogen ke semua arah.

Soal 2. Turunkan dVp = 4p2 dp.

Jawab 2. Dengan mengambil analogi seperti transformasi dari ruang spasial


dengan sistem koordinat kartesian ke sistem koordinat bola, maka dapat dit-
uliskan bahwa

dVp = dpx dpy dpz = (dp)(pd)(p sin )d.

Apabila momentum p bersifat isotropik, maka dapat dilakukan integral terhadap


variabel d dan d sehingga dapat diperoleh

Z Z 2
dVp = sin d d p2 dp = 4p2 dp.
0 0
6.3. INTEGRAL VOLUME RUANG LAJU 33

Dengan demikian dapat dituliskan bahwa

d = 4V p2 dp.

6.3 Integral volume ruang laju

Hubungan antara momentum dan laju adalah

p = mv dp = mdv

sehingga dapat diperoleh

d = 4V m3 v 2 dv.

Soal 3. Turunkan d = 4V m3 v 2 dv.

Jawab 3. Gunakan hubungan p = mv dan dp = mdv dalam d = 4V p2 dp.

6.4 Integral volume ruang energi

Energi setiap partikel dalam bentuk energi kinetik terkait dengan momentumnya
adalah melalui hubungan

p2
=
2m

sehingga dapat dituliskan bahwa

pdp
d = .
m

Soal 4. Rumuskan d dalam bentuk d.

Jawab 4. Dengan menggunakan = p2 /2m, d = pdp/m, dan d = 4V p2 dp,


dapat diperoleh

d = 4V (p2 )(dp)
 
m
d = 4V (2m) d
2m
d = 2V (2m)3/2 1/2 d
34 CATATAN 6. DEGENERASI DALAM RUANG FASA

6.5 Integral volume ruang frekuensi

Khusus untuk partikel yang merupakan foton, maka energinya dirumuskan se-
bagai

= h

sehingga

d = hd.

Perlu diingat bahwa foton tidak memiliki massa sehingga momentumnya adalah
p = h/c.

Soal 5. Rumuskan d dalam bentuk d.

Jawab 5. Dengan menggunakan d = 4V p2 dp dapat dituliskan bahwa

2 
h3 2
 
h hd
d = 4V (p2 )(dp) = 4V = 4V d.
c c c3

6.6 Integral volume ruang panjang gelombang

Selain dalam ruang frekuensi, untuk partikel yang merupakan foton, dapat pula
d dinyatakan dalam ruang panjang gelombang , dengan hubungan

c cd
= d = 2 .

Soal 6. Rumuskan d dalam bentuk d.

Jawab 6. Dengan menggunakan 4V h3 2 d/c3 dan = c/ (serta turunan-


nya) dapat diperoleh

h3
d = 4V 3 2 d
c
h3  c 2 4V h3
   
c 21
d = 4V 3 d = d.
c c 4

Bila diambil nilai positifnya dan sebuah foton memiliki dua arah polarisasi,
maka degenerasi gj tiap satuan volume akan menjadi
6.7. REFERENSI 35

gj 4h3
g()d = = d.
V 4

Umumnya hanya tanda negatif akibat penurunan tidak digunakan.

Soal 7. Gas foton memiliki statistik Bose-Einstein. Rumuskan bagaimana


bentuk g() dan n().

Jawab 7. Suatu foton dalam gas foton memiliki dua arah polarisasi sehingga
degenerasinya menjadi dua kali dari degenerasi yang diperoleh dari gj . Selain itu
umumnya jumlah denerasi atau keadaan yang diperbolehkan dinyatakan dalam
tiap satuan volume [2], sehingga

2gj 2d 8
g()d = = 3
= 4 d.
V Vh

Kemudian dengan menggunakan statistik Bose-Einstein dapat dituliskan bahwa

8 1
n()d = g()df () = d .
4 ehc/kT 1

6.7 Referensi
1. Sparisoma Viridi dan Siti Nurul Khotimah, Catatan Kuliah Fisika Statis-
tik, Semester II Tahun 2009/2010, Mei, 24-25 (2010)
2. A. J. Pointon, An Introduction to Statistical Physics for Students, Long-
mans, First Print, 51-55 (1967)
36 CATATAN 6. DEGENERASI DALAM RUANG FASA
Catatan 7

Distribusi Suatu Statistik

Telah diperkenalkan dalam suatu catatan sebelumnya yang berjudul Konfig-


urasi Paling Mungkin Suatu Statistik, bagaimana bentuk distribusi dari ketiga
statistik (Maxwell-Boltzmann, Bose-Einstein, dan Fermi-Dirac).

7.1 Bentuk umum distribusi ketiga statistik

Ketiga statistik memiliki bentuk umum distribusi partikel, yaitu

gj
Nj,MB = , (7.1)
e (+ j)

untuk statistik Maxwell-Boltzmann,

gj
Nj,BE = , (7.2)
e (+ j) 1

untuk statistik Bose-Einstein, dan

gj
Nj,FD = . (7.3)
e (+ j) +1

untuk statistik Fermi-Dirac. Ketiga bentuk dalam Persamaan (7.1), (7.2), dan
(7.3) dapa dituliskan dalam bentuk diferensialnya NX ()d, di mana X = MB,
BE, dan FD.

37
38 CATATAN 7. DISTRIBUSI SUATU STATISTIK

7.2 Statistik Maxwell-Boltzmann

Soal 1. Dengan menggunakan d = 2V (2m)3/2 1/2 d, gj d/h3 , =


1/kT , tentukanlah NMB ()d.

Jawab 1. Dapat dituliskan bahwa

gj 1
Nj,MB = = gj (+ )
e(+j ) e j

2V (2m)3/2 1/2 d 1
NMB ()d =
h3 e(+)
3/2
2V (2m)
NMB ()d = e e 1/2 d
h3
2V (2m)3/2 /kT 1/2
NMB ()d = e e d. (7.4)
h3

Soal 2. Bila diketahui bahwa bentuk distribusi Maxwell-Boltzmann dalam


bentuk diferensial secara lengkap adalah [1]

2N
NMB ()d = e/kT 1/2 d, (7.5)
(kT )3/2

tentukanlah nilai dari Persamaan (7.4).

Jawab 2. Dari Persamaan (7.4) dan (7.5) dapat dituliskan bahwa

2V (2m)3/2 2N
e =
h3 (kT )3/2
N h3
e =
V (2mkT )3/2
N h3
 
= ln . (7.6)
V (2mkT )3/2

Soal 3. Tunjukkan bahwa

Z
NMB ()d = N. (7.7)
0

Jawab 3. Dapat dituliskan bahwa


2N
Z Z
NMB ()d = e/kT 1/2 d
0 0 (kT )3/2
7.2. STATISTIK MAXWELL-BOLTZMANN 39


2N
Z
= e/kT 1/2 d.
(kT )3/2 0

Kemudian dengan menggunakan


  Z
1 1
= ex x 2 dx = ,
2 0

dan


1
1
Z Z
1 1
ex x 2 dx = ex x 2 dx = ,
0 2 0 2

dapat diperoleh

Z Z  1/2   1
e/kT 1/2 d = (kT )3/2 e/kT d = (kT )3/2 ,
0 0 kT kT 2

sehingga

2N
2N 1
Z
e/kT 1/2 d = (kT )3/2 = N.
(kT )3/2 0 (kT )3/2 2

Jadi, Persamaan (7.7) telah


R dapat dibuktikan. Sebenarnya nilai dapat di-
cari karena syarat bahwa 0 NMB ()d = N dari Persamaan (7.4) tanpa perlu
terlebih dahulu mengetahui bentuk lengkapnya seperti dalam Persamaan (7.5).

Soal 4. Gas ideal monoatomik memenuhi distribusi Maxwell-Boltzmann dalam


Persamaan (7.5). HitunglahR energi total sistem yang terdiri dari N partikel gas

dengan menggunakan U = 0 NMB ()d.

Jawab 4. Dapat dituliskan


2N
Z Z
U= NMB ()d = e/kT 3/2 d,
0 (kT )3/2 0

di mana

Z Z  3/2  
e/kT 3/2 d = (kT )5/2 e/kT d
0 0 kT kT
3 1 3
= (kT )5/2 = (kT )5/2 ,
2 2 4
40 CATATAN 7. DISTRIBUSI SUATU STATISTIK

sehingga

2N 3 3
U= 3/2
(kT )5/2 = N kT.
(kT ) 4 2

7.3 Statistik Bose-Einstein

Soal 5. Dengan menggunakan d = 2V (2m)3/2 1/2 d, gj d/h3 ,


= 1/kT , tentukanlah NBE ()d untuk foton dalam suatu ruang tertutup
berlubang, di mana foton memiliki dua arah polarisasi yang akan mempen-
garuhi jumlah keadaan energinya dan bahwa jumlah foton tidak tetap (ada
yang diserap dan dipancarkan kembali oleh dinding). Lubang pada ruang ter-
tutup tersebut akan berfungsi sebagai benda hitam. Apakah ada yang salah
dari hasil yang diperoleh?

Jawab 5. Dapat diperoleh bahwa

gj 1
Nj,BE = = gj
e (+ j)
1 e(+j )
1
2V (2m)3/2 1/2 d 1
NBE ()d = 2
h3 e(+) 1
3/2
4V (2m) 1
NBE ()d = 3
1/2 d
h e (+) 1
4V (2m)3/2 1
NBE ()d = 1/2 d. (7.8)
h3 e(/kT ) 1

Sebelum mencari nilai , sebaiknya hasil yang diperoleh dicermati terlebih


dahulu. Sekilas terlihat bahwa tidak ada yang salah dalam Persamaan (7.8)
akan tetapi perhatikan bahwa dalam bentuk NBE ()d yang dicari untuk foton
terdapat massa foton. Tidak ada arti fisis dari massa foton. Dengan demikian
ungkapan dalam Persamaan (7.8) adalah salah atau NBE tidak dapat dinyatakan
dalam bentuk seperti di atas. Ungkapan yang benar adalah apabila dinyatakan
dalam panjang gelombang atau frekuensi dari foton.

Soal 6. Perbaikilah Persamaan (7.8) dengan mencari NBE ()d menggunakan


representasi d dalam d.

Jawab 6. Dalam catatan sebelumnya dapat diperoleh bahwa

gj 4h3
g()d = = d
V 4

akan tetapi karena foton memiliki dua arah polarisasi yang menyebabkan jumlah
keadaan energi yang dimilikinya menjadi dua kalinya, maka ungkapan di atas
akan menjadi
7.4. STATISTI FERMI-DIRAC 41

8h3
d.
4

Ungkapan-ungkapan di atas diperoleh melalui hubungan = hc/ dan turunan-


nya. Dengan demikian dapat dituliskan

gj 1
Nj,BE = = gj
e 1
(+ j) 1e(+j )
8 1
NBE ()d = 4 d (+hc/)
e 1
8 1
NBE ()d = 4 (hc/kT ) d. (7.9)
e 1

Ungkapan dalam Persamaan (7.9) sudah dalam per satuan volume V . Selan-
jutnya adalah bagaimana mencari nilai . Dalam soal diinformasikan bahwa
jumlah foton dalam sistem tidak tetap karena ada foton yang diserap oleh
wadah tertutup dan ada foton yang dipancarkan kembali setelah diserap, dengan
demikian pada saat penurunan Nj,BE menggunakan pengali tak tentu Lagrange

d ln W + dN + dU = 0

tidak dapat dipenuhi bahwa dN = 0. Agar persamaan di atas dapat tetap


dipenuhi, dipilih = 0. Dengan demikian akan diperoleh untuk foton dalam
suatu ruang tertutup

8 1
NBE ()d = d. (7.10)
4 ehc/kT 1

7.4 Statisti Fermi-Dirac

Soal 7. Dengan menggunakan d = 2V (2m)3/2 1/2 d, gj d/h3 , =


1/kT , tentukanlah NFD ()d untuk gas elektron yang memiliki dua kemungk-
inan spin, yaitu + 21 dan 21 . Gunakan pula hubungan bahwa = F /kT .

Jawab 7. Dapat dituliskan bahwa

gj 1
Nj,FD = = gj (+ )
e(+j ) + 1 e j + 1

3/2 1/2
2V (2m) d 1
NFD ()d = 2
h3 e(+) + 1
42 CATATAN 7. DISTRIBUSI SUATU STATISTIK

4V (2m)3/2 1
NFD ()d = 3
1/2 d
h e (+) +1
"  3 #
2m 2 1/2 1
NFD ()d = V 4 d. (7.11)
h2 e(F )/kT + 1

Khusus untuk statistik Fermi-Dirac, distribusi partikel (dalam hal ini elektron)
dapat dituliskan dalam bentuk

N ()d = g()f ()d,

di mana

"   23 #
2m 1/2
g() = V 4
h2

dan

1
f () =
e(F )/kT + 1

yang dikenal sebagai fungsi Fermi. F disebut sebagai energi Fermi.

7.5 Bentuk umum distribusi statistik lain

Walaupun tidak lazim dituliskan, secara umum ketiga statistik seharusnya dapat
dituliskan dalam bentuk

N ()d = g()f ()d, (7.12)

dengan g() memiliki arti jumlah keadaan energi pada tiap tingkat energi atau
kerapatan keadaan energi (density of states).

7.6 Referensi
1. A. J. Pointon, An Introduction to Statistical Physics for Students, Long-
mans, First Print, 25-26 (1967)
Catatan 8

Termodinamika Gas Ideal


Monoatomik

Dalam gas ideal segala interaksi yang terjadi antara partikel-partikel gas, terma-
suk yang terjadi saat partikel-partikel gas saling bertumbukan, dianggap mem-
berikan pengaruh yang dapat diabaikan terhadap sifat-sifat termodinamika gas
[1].

8.1 Peluang termodinamika Wmaks gas ideal


klasik

Peluang termodinamika suatu keadaan makro dari gas ideal yang mengandung
N partikel gas tak berstruktur adalah

Y gjNj
W = N! , (8.1)
j
Nj !

dengan Nj adalah bilangan okupasi pada tingkat energi j, di mana tingkat


energi tersebut
P terdegenerasi sejumlah gj dan berenergi j . Terpenuhi pula
bahwa N = j Nj .

Soal 1. Gunakan aproksimasi Stirling ln x! x ln x x untuk mencari ln W.


Q P
Jawab 1. Bentuk ln j xj = j ln xj , sehingga


Y gjNj
ln W = ln N !
j
Nj !

43
44 CATATAN 8. TERMODINAMIKA GAS IDEAL MONOATOMIK

Y gjNj
= ln N ! + ln
j
Nj !
Nj
!
X gj
= ln N ! + ln
j
Nj !
X N

= ln N ! + ln gj j ln Nj !
j
X X
= ln N ! + Nj ln gj ln Nj !
j j
X X
(N ln N N ) + Nj ln gj (Nj ln Nj Nj )
j j
X X X
= N ln N N + Nj ln gj Nj ln Nj + Nj
j j j
X X
= N ln N N + Nj ln gj Nj ln Nj + N
j j
X X
= N ln N + Nj ln gj Nj ln Nj
j j
X gj
= N ln N + Nj ln . (8.2)
j
Nj

Soal 2. Dengan menggunakan rumusan untuk mencari Wmaks pada statis-


tik Maxwell-Boltzman yang memberikan Nj = gj /e(+j ) tentukan bentuk
ln Wmaks dalam N dan U .

Jawab 2. Dengan menggunakan Persamaan (8.2) untuk Wmaks sehingga


hubungan Nj = gj /e(+j ) dapat digunakan, diperoleh

X gj
ln W = N ln N + Nj ln
j
Nj
X
ln Wmaks = N ln N + Nj ln e(+j )
j
X
= N ln N + Nj [( + j )]
j
X X
= N ln N Nj Nj j
j j
= N ln N N U. (8.3)

Soal 3. Ubahlah Persamaan (8.3) dengan menggunakan A = e dan =


1/kT .

Jawab 3. Dapat dituliskan bahwa


8.2. FUNGSI PARTISI BOLTZMANN 45

ln Wmaks = N ln N N U
U
= N ln N (ln A)N +
kT
N U
= N ln + .
A kT

8.2 Fungsi partisi Boltzmann

Suatu fungsi partisi didefinisikan sebagai Z = N/A.


P
Soal 4. Rumuskan bentuk fungsi partisi Z dengan menggunakan j Nj dan
e .

Jawab 4. Dapat dituliskan bahwa

P
N j Nj
Z=
A e
(+j )
P
j gj /e
=
e
(+j )
P
j gj e
=
e
j
P
e j gj e
=
e
X X
j
= gj e = gj ej /kT . (8.4)
j j

Persamaan (8.4) ini disebut sebagai fungsi partisi Boltzmann (atau fungsi par-
tisi) sebuah partikel dalam suatu suatu sistem. Istilah ini digunakan karena
dalam ekspresi Z, setiap suku dalam somasi mementukan bagaimana partikel
dalam sistem didistribusikan atau dipartisikan di antara (pada) tingkat-tingkat
energi.

Soal 5. Dengan menggunakan perumusan Boltzmann untuk entropi, tentukan


bentuk dari S yang bergantung pada Z.

Jawab 5. Perumusan Boltzmann untuk entropi adalah S = k ln Wmaks sehingga

S = k ln Wmaks
 
N U
= k N ln +
A kT
 
U
= k N ln Z +
kT
46 CATATAN 8. TERMODINAMIKA GAS IDEAL MONOATOMIK

U
= N k ln Z + . (8.5)
T

Soal 6. Energi bebas Helmholtz didefinisikan sebagai F = U T S. Gunakan


ekspresi tersebut untuk membuat fungsi F = F (N, T, Z)/

Jawab 6. Dengan menggunakan Persamaan (8.5) dapat dituliskan

U
S = N k ln Z +
T
T S = N kT ln Z + U
N kT ln Z = U T S
N kT ln Z = F. (8.6)

Soal 7. Hitunglah energi dalam U dari energi bebas Helmholtz F dengan


menggunakan rumusan

 
(F/T )
U = T 2 (8.7)
T V

3
dan Z = V (2mkT ) 2 /h3 .

Jawab 7. Dengan menggunakan Persamaan (8.6) dan (8.7) dapat dituliskan

 
(F/T )
2
U = T
T
 V
(N kT ln Z/T )
= T 2
T
 V
(N k ln Z)
= T2
T
" 3
#V
3
ln[V (2mkT ) 2 /h ]
= T 2N k
T
" 3
#V
3
ln[V (2mkT ) ] ln h
2
= T 2N k
T T
V
 
1 3 1
= T 2N k 3 (2mkT ) 2 2mk 0

V (2mkT ) 2 2
 
3
= T 2N k
2T
3
= N kT.
2

Soal 8. Turunkan ekspresi


8.3. TEKANAN DAN KALOR JENIS 47

 
2 ln Z
U = N kT (8.8)
T V

Nj = gj Aej /kT , dan Z =


P P
dari U = N , N = j Nj , U = j j Nj ,
j /kT
P
j gj e .

Jawab 8. Dengan menggunakan persamaan-persamaan di atas dapat dituliskan


bahwa

P
U j j Nj
U = N = N =N P
N j Nj
j /kT j /kT
P P
j j gj Ae j j g j e
=N P /kT
=N P j /kT
j gj Ae j gj e
j

j /kT
 
NX j /kT N X 2 (e )
= j g j e = gj kT
Z j Z j T
 
N X N Z
= kT 2 gj ej /kT = kT 2
Z T j Z T V
 
ln Z
= N kT 2 .
T V

Soal 9. Turunkan Persamaan (8.8) dari Persamaan (8.7).

Jawab 9. Dengan menggunakan kedua persamaan yang disebutkan dalam soal,


dapat dituliskan bahwa

 
(F/T )
U = T 2
T
 V
2 (N kT ln Z/T )
= T
T
 V
(N k ln Z)
= T2
T
 V
ln Z
= N kT 2 .
T V

8.3 Tekanan dan kalor jenis

Melalui definisi energi bebas Helmholtz

F = U TS
48 CATATAN 8. TERMODINAMIKA GAS IDEAL MONOATOMIK

dapat diperoleh bahwa secara umum

dF = dU T dS SdT.

Dengan menggunakan hubungan

dQ = dU + pdV

dan

dQ = T dS

maka dapat diperoleh

dF = pdV SdT. (8.9)

Dari persamaan terakhir ini dapat diturunkan p dan S sebagai fungsi dari F .

Soal 10. Tentukanlah ungkapan p dan S dari F .

Jawab 10. Dengan menggunakan Persamaan (8.9) dapat dituliskan bahwa

 
F
p= (8.10)
V T

dan

 
F
S= . (8.11)
T V

Soal 11. Tentukanlah ungkapan U sebagai fungsi dari F . Bila perlu gunakan
pula hubungan = 1/kT .

Jawab 11. Dengan menggunakan Persamaan (8.11) dan definisi energi bebas
Helmholtz F = U T S dapat diperoleh

 
F
F = U TS = U + T
T V
     
F (F/T ) (F )
U =F T = T 2 = (8.12)
T V T V V

Soal 12. Dengan menggunakan definisi dari kapasitas panas pada volume tetap
8.4. PERSAMAAN KEADAAN 49

 
U
CV = (8.13)
T V

tentukanlah CV dari F . Bila perlu gunakan pula hubungan = 1/kT .

Jawab 12. Dengan segera dapat diperoleh bahwa

 
U
CV =
T
    V
(F/T )
= T 2
T T V
     V
(F/T ) (F/T )
= 2T T2
T V T T
        VV
F 1 F 2 F 1 F
= 2T +T
T2 T T V T T 2 T T V
       2 V
F F F F F F
=2 2 + 2 + T
T T V T V T T V T 2 V
 2   2 
F 2 (F )
= T = k .
T 2 V 2 V

8.4 Persamaan keadaan

Dengan menggunakan

3
V (2mkT ) 2
Z= ,
h3

 
F
p=
V T

dan

F = N kT ln Z

dapat diperoleh

N kT
p=
V

yang merupakan persamaan keadaan gas ideal monotomik.


50 CATATAN 8. TERMODINAMIKA GAS IDEAL MONOATOMIK

8.5 Referensi
1. A. J. Pointon, An Introduction to Statistical Physics for Students, Long-
mans, First Print, 86-93 (1967)
Catatan 9

Paradoks Gibb

Saat dua jenis gas berbeda dengan entropi masing-masing dicampur, maka en-
tropi campuran adalah penjumlahan kedua entropi semula. Lalu bagaimana
apabila kedua gas tersebut adalah jenis yang sama? Ternyata entropinya
bukanya hanya penjumlahan dari kedua entropi semula melainkan terdapat su-
atu suku tambahan. Untuk itu perumusan gas klasik perlu diperbaiki den-
gan menggunakan perumusan semi-klasik [1]. Dalam catatan ini gas yang
dibicarakan adalah gas ideal monoatomik tanpa adanya struktur di dalamnya.

9.1 Entropi gas klasik

Dengan menggunakan perumusan entropi S dari energi bebas Helmholtz F

 
F
S= , (9.1)
T V

kaitan antara energi bebas Helmholtz F dengan fungsi partisi Z

F = N kT ln Z, (9.2)

dan bentuk eksplisit fungsi partisi Boltzmann

3
V (2mkT ) 2
Z= , (9.3)
h3

dapat diperoleh bentuk eksplisit dari entropi yang bergantung dari jumlah par-
tikel N , volume gas V , dan temperatur gas T , yaitu

51
52 CATATAN 9. PARADOKS GIBB

" 3
#
V (2mkT ) 2 3
S = N k ln + N k. (9.4)
h3 2

Di sini m adalah massa satu partikel gas, k adalah konstanta Boltzmann, dan
h adalah konstanta Planck.

Soal 1. Turunkan Persamaan (9.4).

Jawab 1. Dapat dituliskan bahwa

 
F
S=
T V
   
(N kT ln Z) (N kT ln Z)
= =
T V T
 V
N kT Z
= N k ln Z +
Z T V
" 3
# ( " 3
#)
V (2mkT ) 2 N kT V (2mkT ) 2
= N k ln +
h3 3
V (2mkT ) 2 /h3 T h3
V
" 3
# 3 3
3
V (2mkT ) 2 N kT h V (2mk) d(T )
2 2
= N k ln +
h3 V (2mkT ) 2
3
h3 dT
" 3
#
V (2mkT ) 2 N kT 3 1
= N k ln 3
+ 3 T2
h T2 2
" 3
#
V (2mkT ) 2 3
= N k ln + N k.
h3 2

9.2 Pencampuran dua gas berbeda jenis

Sebuah sistem terdiri dari dua ruangan yang masing-masing terisi oleh satu
jenis gas. Gas 1 yang memiliki jumlah partikel N1 , dengan massa tiap par-
tikel m1 , menempati ruangan bervolume V , bertemperatur T , dan bertekanan
p. Sedangkan gas 2 yang menempati ruangan bervolume, bertemperatur, dan
bertekanan sama dengan gas 1, akan tetapi memiliki jumlah partikel N2 dan
massa tiap partikelnya adalah m2 . Terdapat sekat yang memisahkan ruangan
kedua jenis gas tersebut.

Soal 2. Hitunglah entropi total sistem sebelum kedua jenis gas bercampur.

Jawab 2. Dengan menggunakan Persamaan (9.4) dapat dihitung entropi


masing-masing gas, yaitu S1 dan S2 dan entropi total sistem S
9.2. PENCAMPURAN DUA GAS BERBEDA JENIS 53

" 3
#
V (2mkT ) 2 3
S = N k ln 3
+ N k,
h 2
" 3
#
V (2m1 kT ) 2 3
S1 = N1 k ln + N1 k,
h3 2
" 3
#
V (2m2 kT ) 2 3
S1 = N2 k ln 3
+ N2 k,
h 2
S = S1 + S2 .

Soal 3. Hitunglah entropi total sistem setelah kedua jenis gas bercampur.

Jawab 3. Setelah sekat pemisah ruangan kedua jenis gas dihilangkan maka
kedua jenis gas akan bercampur. Mengingat tekanan dan temperatur awal kedua
gas adalah sama, maka partikel-partikel kedua gas akan memiliki temperatur
dan tekanan campuran yang sama pula. Hanya saja setelah tercampur, masing-
masing partikel kedua gas akan melihat voume ruangan menjadi dua kali volume
semula. Dengan demikian

" 3
#
V (2mkT ) 2 3
S = N k ln + N k,
h3 2
" 3
#
2V (2m1 kT ) 2 3
S1 = N1 k ln 3
+ N1 k,
h 2
" 3
#
2V (2m2 kT ) 2 3
S1 = N2 k ln + N2 k,
h3 2
S = S1 + S2 .

Soal 4. Hitunglah perubahan entropi sistem S.

Jawab 4. Perubahan entropi sistem S = S S sehingga

S = S S = (S1 + S2 ) (S1 + S2 ) = (S1 S1 ) + (S2 S2 )


= S1 + S2 .
( " 3
# )
2V (2m1 kT ) 2 3
S1 = N1 k ln + N1 k
h3 2
( " 3
# )
V (2m1 kT ) 2 3
N1 k ln + N1 k
h3 2
= N1 k ln 2.
S2 = N2 k ln 2.
S = (N1 + N2 )k ln 2.
54 CATATAN 9. PARADOKS GIBB

9.3 Pencampuran gas sejenis: paradoks Gibb

Bagaimana bila gas yang dicampur memiliki jenis yang sama? Suatu fenomena
yang disebut sebagai paradoks Gibb muncul di sini. Sistem yang ditinjau sama
dengan sistem sebelumnya, hanya saja dalam hal ini kedua gas berjenis sama.
Dan karena dijaga agar tekanan p, temperatur T , dan volume V sama, maka
dengan m1 = m2 = m akan terpenuhi bahwa N1 = N2 = N .

Soal 5. Hitunglah entropi total sistem sebelum kedua gas berjenis sama bercam-
pur.

Jawab 5. Dengan menggunakan Persamaan (9.4) dapat dihitung entropi


masing-masing gas, yaitu S1 dan S2 dan entropi total sistem S

" 3
#
V (2mkT ) 2 3
S = N k ln + N k,
h3 2
" 3
#
V (2mkT ) 2 3
S1 = N k ln + N k,
h3 2
" 3
#
V (2mkT ) 2 3
S2 = N k ln + N k,
h3 2
S = S1 + S2 = 2S1 = 2S2 .

Soal 6. Hitunglah entropi total sistem setelah kedua gas berjenis sama bercam-
pur.

Jawab 6. Setelah sekat pemisah ruangan kedua jenis gas dihilangkan maka
kedua gas akan bercampur. Mengingat tekanan dan temperatur awal kedua
gas adalah sama, maka partikel-partikel kedua gas akan memiliki temperatur
dan tekanan campuran yang sama pula. Hanya saja setelah tercampur, masing-
masing partikel kedua gas akan melihat voume ruangan menjadi dua kali volume
semula. Dengan demikian

" 3
#
V (2mkT ) 2 3
S = N k ln + N k,
h3 2
" 3
#
(2V )(2mkT ) 2 3
S1 = N k ln + N k,
h3 2
" 3
#
(2V )(2mkT ) 2 3
S1 = N k ln + N k,
h3 2
S = S1 + S2 = 2S1 = 2S2 .

Soal 7. Hitunglah perubahan entropi sistem S.


9.4. GAS IDEAL SEMI-KLASIK 55

Jawab 7. Perubahan entropi sistem S

S = S1 + S2 .
( " 3
# )
2V (2mkT ) 2 3
S1 = N k ln + Nk
h3 2
( " 3
# )
V (2mkT ) 2 3
N k ln + Nk
h3 2
= N k ln 2.
S2 = N k ln 2.
S = 2N k ln 2.

Di sini diperoleh bahwa S = 2N k ln 2 dan bukan S = 0, padahal kedua gas


adalah jenis gas yang sama. Ketidakcocokan ini disebut sebagai paradok Gibb.

9.4 Gas ideal semi-klasik

Peluang suatu keadaan makro gas ideal klasik yang semula menggunakan statis-
tik Maxwell-Boltzmann

Y gjNj
WMB = N !
j
Nj !

dapat dikoreksi dengan menggunakan statistik kuantum, yang seharusnya tetap


memperhatikan sifat statistik dari partikel gas apakah bersifat sebagai boson
atau fermin, sehingga menjadi menjadi peluang termodinamika suatu keadaan
makro semi-klasik

Y gjNj
WSK = . (9.5)
j
Nj !

Dengan menggunakan dua pengali tak tentu Lagrange dan dapat diperoleh
bahwa

U
Wmaks = N + N,
kT

dan dengan S = k ln Wmaks dapat dituliskan

" 3
#
V (2mkT ) 2 5
S = N k ln + N k, (9.6)
N h3 2
56 CATATAN 9. PARADOKS GIBB

3
dengan menggunakan fungsi partisi Boltzmann yang sama Z = V (2mkT ) 2 /h3 .

Soal 8. Hitunglah entropi sistem yang terdiri dari dua gas berjenis sama seperti
dalam soal sebelumnya, saat sebelum dan sudah dicampur. Hitung pula peruba-
han entropinya.

Jawab 8. Saat sebelum dicampur, dengan menggunakan Persamaan (9.6) dapat


diperoleh

" 3
#
V (2mkT ) 2 5
S = N k ln + N k,
N h3 2
" 3
#
V (2mkT ) 2 5
S1 = N k ln + N k,
N h3 2
" 3
#
V (2mkT ) 2 5
S2 = N k ln + N k,
N h3 2
S = S1 + S2 = 2S1 = 2S2 ,

sedangkan saat setelah dicampur

" 3
#
V (2mkT ) 2 5
S = N k ln + N k,
N h3 2
" 3
#
(2V )(2mkT ) 2 5
S = S1 + S2 = (2N )k ln + (2N )k,
(2N )h3 2
= 2S2 = 2S1 ,

sehingga perubahan entropinya menjadi

S = S S,
= 2S1 2S1 = 2S2 2S2 = 0.

Dengan menggunakan statistik semi-klasik, telah ditunjukkan bahwa paradoks


Gibb tidak lagi muncul saat dua gas berjenis sama dicampurkan.

9.5 Referensi
1. A. J. Pointon, An Introduction to Statistical Physics for Students, Long-
mans, First Print, 93-99 (1967)
Catatan 10

Statistik Fermi-Dirac: N j
dan S

[20100629] Peluang termodinamika suatu keadaan makro-k dalam sistem yang


memenuhi statistik Fermi-Dirac diberikan oleh

Y gj !
Wk = , (10.1)
j
nj !(gj nj )!

dengan gj adalah degenerasi tingkat energi j yang memiliki energi j dalam


keadaan makro k dan nj adalah jumlah partikel yang menempati tingkat energi j
juga dalam keadaan makro k tersebut. Dalam statistik Fermi-Dirac hanya boleh
terdapat satu partikel atau tidak ada partikel yang menempati satu keadaan
energi. Jumlah keadaan energi dalam satu tingkat energi j ditunjukkan dengan
nilai degenerasi tingkat energi tersebut gj .

Bilangan okupasi rata-rata setiap tingkat energi j dapat diperoleh lewat

1X
Nj = Wk Njk . (10.2)

k

Terdapat suatu sistem yang terdiri dari 5 partikel mematuhi statistik Fermi-
Dirac. Terdapat empat tingkat energi yang diperhitungkan, yaitu 1 = 2, 2 =
3, 3 = 4, dan 4 = 5. Degenerasi masing-masing tingkat energi bergantung
dari volume sistem V dan energi total sistem tergantung dari temperatur sistem
T.

57
58 CATATAN 10. STATISTIK FERMI-DIRAC: N J DAN S

10.1 Soal

1. Tuliskan semua kemungkinan kelima partikel tersebut didistribusikan pada


keempat tingkat energi sehingga memberikan U = 19 dan U = 17 tanpa
memperhatikan statistik yang digunakan.

2. Pada titik a dalam ruang parameter V T , sistem memiliki temperatur


Ta , volume Va , dan energi total Ua = 19. Degenerasi tingkat-tingkat
energi sistem pada keadaan ini adalah g1 = 1, g2 = 3, g3 = 4, dan g4 = 6.
Lengkapilah tabel berikut ini.

k
j j / gj 1 2 3 4 5 6 7 Nj
Njk
4 5
3 4
2 3
1 2

Wk

Hitunglah entropi sistem Sa dengan menggunakan rumusan Planck.

3. Pada titik a dalam ruang parameter V T , sistem memiliki temperatur


Tb = Ta , volume Vb < Va , dan energi total Ua = 19. Degenerasi tingkat-
tingkat energi sistem pada keadaan ini adalah g1 = 1, g2 = 2, g3 = 3, dan
g4 = 5. Lengkapilah tabel berikut ini.

k
j j / gj 1 2 3 4 5 6 7 Nj
Njk
4 5
3 4
2 3
1 2

Wk

Hitunglah entropi sistem Sb dengan menggunakan rumusan Planck.

4. Pada titik a dalam ruang parameter V T , sistem memiliki temperatur


Tc < Tb , volume Vc < Vb , dan energi total Ua = 17. Degenerasi tingkat-
tingkat energi sistem pada keadaan ini adalah g1 = 1, g2 = 2, g3 = 3, dan
g4 = 5. Lengkapilah tabel berikut ini.
10.1. SOAL 59

k
j j / gj 1 2 3 4 5 6 7 Nj
Njk
4 5
3 4
2 3
1 2

Wk

Hitunglah entropi sistem Sc dengan menggunakan rumusan Planck.
5. Pada titik a dalam ruang parameter V T , sistem memiliki temperatur
Td = Tc , volume Vd = Va , dan energi total Ua = 17. Degenerasi tingkat-
tingkat energi sistem pada keadaan ini adalah g1 = 1, g2 = 3, g3 = 4, dan
g4 = 6. Lengkapilah tabel berikut ini.
k
j j / gj 1 2 3 4 5 6 7 Nj
Njk
4 5
3 4
2 3
1 2

Wk

Hitunglah entropi sistem Sd dengan menggunakan rumusan Planck.
6. Gambarkan keempat titik a, b, c, dan d dalam ruang parameter V T
dan tentukanlah proses dari titik mana ke titik mana yang mungkin terjadi
apabila hanya entropi sistem yang ditinjau. Apa syaratnya?
60 CATATAN 10. STATISTIK FERMI-DIRAC: N J DAN S

10.2 Jawab
1. Agar diperoleh U = 19 kelima partikel dapat disusun seperti tampak
dalam Tabel 10.1 berikut. Sedangkan untuk U = 17 dapat dilihat dalam
Tabel 10.2.

Tabel 10.1: Susunan yang mungkin kelima partikel pada empat tingkat energi
dengan U = 19.
k
j j / 1 2 3 4 5 6
Njk
4 5 2 1 0 1 2 3
3 4 1 3 4 2 0 0
2 3 1 0 1 2 3 0
1 2 1 1 0 0 0 2
Uk / 19 19 19 19 19 19

Tabel 10.2: Susunan yang mungkin kelima partikel pada empat tingkat energi
dengan U = 17.
k
j j / 1 2 3 4 5 6
Njk
4 5 0 1 2 1 0 1
3 4 3 1 0 0 2 2
2 3 1 2 1 4 3 0
1 2 1 1 2 0 0 2
Uk / 17 17 17 17 17 17

2. Dengan degenerasi tingkat-tingkat energi sistem pada Ta dan Va adalah


g1 = 1, g2 = 3, g3 = 4, dan g4 = 6 maka dapat diperoleh penempatan
yang mungkin bagi kelima partikel adalah seperti dalam Tabel 10.3.

6! 4! 3! 1!
W1 = = 15 4 3 1 = 180
2!(6 2)! 1!(4 1)! 1!(3 1)! 1!(1 1)!
6! 4! 3! 1!
W2 = = 6 4 1 1 = 24
1!(6 1)! 3!(4 3)! 0!(3 0)! 1!(1 1)!
6! 4! 3! 1!
W3 = =1131= 3
0!(6 0)! 4!(4 4)! 1!(3 1)! 0!(1 0)!
6! 4! 3! 1!
W4 = = 6 6 3 1 = 108
1!(6 1)! 2!(4 2)! 2!(3 2)! 0!(1 0)!
6! 4! 3! 1!
W5 = = 15 1 1 1 = 15
2!(6 2)! 0!(4 0)! 3!(3 3)! 0!(1 0)!

= 180 + 24 + 3 + 108 + 15 = 330


10.2. JAWAB 61

1 204
N1 = (180 1 + 24 1 + 3 0 + 108 0 + 15 0) = = 0.618
270 330
1 444
N2 = (180 1 + 24 0 + 3 1 + 108 2 + 15 3) = = 1.345
270 330
1 480
N3 = (180 1 + 24 3 + 3 4 + 108 2 + 15 0) = = 1.455
270 330
1 522
N4 = (180 2 + 24 1 + 3 0 + 108 1 + 15 2) = = 1.582
270 330

N = 0.618 + 1.345 + 1.455 + 1.582 = 5

Tabel 10.3: Susunan yang mungkin kelima partikel pada empat tingkat energi
dengan U = 19 dan g1 = 1, g2 = 3, g3 = 4, dan g4 = 6.
k
j j / gj 1 2 3 4 5 6 Nj
Njk
4 5 6 2 1 0 1 2 - 0.618
3 4 4 1 3 4 2 0 - 1.345
2 3 3 1 0 1 2 3 - 1.455
1 2 1 1 1 0 0 0 - 1.582

Wk 180 24 3 108 15 - 330


3. Dengan degenerasi tingkat-tingkat energi sistem pada Tb = Ta dan Vb < Va


adalah g1 = 1, g2 = 2, g3 = 3, dan g4 = 5 maka dapat diperoleh pen-
empatan yang mungkin bagi kelima partikel adalah seperti dalam Tabel
10.4.

5! 3! 2! 1!
W1 = = 10 3 2 1 = 60
2!(5 2)! 1!(3 1)! 1!(2 1)! 1!(1 1)!
5! 3! 2! 1!
W2 = = 5111 =5
1!(5 1)! 3!(3 3)! 0!(2 0)! 1!(1 1)!
5! 3! 2! 1!
W4 = = 5 3 1 1 = 15
1!(5 1)! 2!(3 2)! 2!(2 2)! 0!(1 0)!

= 60 + 5 + 15 = 80

1 65
N1 = (60 1 + 5 1 + 15 0) = = 0.8125
80 80
1 90
N2 = (60 1 + 5 0 + 15 2) = = 1.1250
80 80
62 CATATAN 10. STATISTIK FERMI-DIRAC: N J DAN S

1 105
N3 = (60 1 + 5 3 + 15 2) = = 1.3125
80 80
1 140
N4 = (60 2 + 5 1 + 15 1) = = 1.7500
80 80

N = 0.8125 + 1.1250 + 1.3125 + 1.7500 = 5

Tabel 10.4: Susunan yang mungkin kelima partikel pada empat tingkat energi
dengan U = 19 dan g1 = 1, g2 = 2, g3 = 3, dan g4 = 5.
k
j j / gj 1 2 3 4 5 6 Nj
Njk
4 5 5 2 1 - 1 - - 1.7500s
3 4 3 1 3 - 2 - - 1.3125
2 3 2 1 0 - 2 - - 1.1250
1 2 1 1 1 - 0 - - 0.8125

Wk 60 5 - 15 - - 80

Catatan 11

Tingkat dan Keadaan


Energi

Suatu sistem kuantum memiliki diskritisasi energi. Dapat dibedakan antara


tingkat energi (energy levels) dan keadaan energi (energy sates). Sebagai ilus-
trasi beberapa sistem dengan konfigurasi yang berbeda akan ditunjukkan. Isti-
lah degenerasi pun akan digunakan dalam bab ini.

11.1 Tingkat Energi

Tingkat energi atau level energi (energy level) adalah susunan tingkat-tingkat
di mana energi pada tingkat-tingkat tersebut berbeda. Dalam buku ini suatu
tingkat energi diberi label j dan besar energi pada suatu tingkat adalah j .

11.2 Keadaan Energi

Dalam satu tingkat energi terdapat semacam ruang-ruang yang memiliki energi
hampir sama dan dinamakan sebagai keadaan-keadaan energi.

63
64 CATATAN 11. TINGKAT DAN KEADAAN ENERGI
Catatan 12

Keadaan Makro dan Mikro

65
66 CATATAN 12. KEADAAN MAKRO DAN MIKRO
Catatan 13

Peluang Termodinamika

Dalam suatu sistem yang terisolasi dan tertutup jumlah energi sistem E tetap
dan jumlah partikel dalam sistem N tetap. Dengan berevolusinya waktu, in-
teraksi antar partikel dalam suatu sistem yang terisolasi dan tertutup men-
gakibatkan perubahan jumlah partikel yang menempati suatu tingkat energi
dan dapat juga terjadi perubahan keadaan energi dari setiap partikel. Untuk
sistem berupa gas interaksi yang dimaksud dapat berupa tumbukan antar par-
tikel gas atau dengan wadahnya sedangkan untuk molekul-molekul kristas dapat
berupa pertukaran energi. Berbagai bentuk interaksi ini menghasilkan peruba-
han keadaan mikro dari sistem yang tetap harus memenuhi syarat tetapnya E
dan N .

13.1 Postulat termodinamika statistik

Postulat fundamental dair termodinamika statistik menyatakan bahwa semua


keadaan mikro yang mungkin muncul dari suatu sistem terisolai adalah sama
peluangnya. Terdapat dua cara untuk melakukan intepretasi dari postulat ini.

Cara pertama adalah dengan membayangkan sistem telah diamati dalam suatu
rentang waktu t yang cukup lama sehingga setiap keadan mikro dari suatu sis-
tem yang terisolasi telah muncul amat sering. Bila t adalah total waktu sistem
berada pada suatu keadaan mikro yang mungkin, maka postulat ini menyatakan
bahwa rentang waktu t adalah sama untuk semua keadaan mikro.

Sebagai alternatif, cara kedua dapat dipergunakan di mana dibayangkan ter-


dapat sejumlah salinan atau replika dari sistem (sebuah ensemble) yang jum-
lahnya adalah N . Pada suatu saat pengamatan, terdapat sejumlah N replika
yang berada dalam keadaan mikro yang sama. Postulat termodinamika statistik
menyatakan bahwa jumlah N adalah sama untuk semua keadaan mikro.

Postulat ini terlihat tidak diturunkan suatu prinsip fundamental apapun se-
hingga tidak dapat diverifikasi menggunakan eksperimen. Justifikasi kebenaran

67
68 CATATAN 13. PELUANG TERMODINAMIKA

postulat ini terletak pada ketepatan kesimpulan yang dapat ditarik.

13.2 Peluang termodinamika

Sejumlah keadaan mikro akan membentuk satu keadaan makro. Jumlah dari
semua keadaan mikro yang mungkin bagi suatu keadaan makro k disebut sebagai
peluang termodinamika Wk dari keadaan makro tersebut. Suatu asembli dengan
banyak partikel, peluang termodinamika akan bernilai besar.

Jumlah total keadaan mikro yang mungkin untuk suatu asembli, atau dapat
dikatakan sebagai peluang termodinamika asembli tersebut, adalah jumlah pelu-
ang termodinammika keadaan makro dari semua keadaan makro dalam asembli
tersebut

X
= Wk . (13.1)
k

Persamaan (13.1) dapat dijelaskan dengan ilustrasi sebagai berikut. Misalkan


saja dalam suatu sistem terdapat keadaan mikro. Jumlah keadaan mikro
yang dapat membentuk keadaan makro pertama adalah W1 (peluang termodi-
namika keadaan makro pertama), jumlah keadaan mikro yang dapat memben-
tuk keadaan makro kedua adalah W2 (peluang termodinamika keadaan makro
kedua), dan seterusnya. Dengan demikian jumlah seluruh keadaan mikro dalam
sistem tersebut tak lain adalah jumlah peluang termodinammika keadaan makro
dari semua keadaan makro dalam asembli tersebut.

Untuk sistem dengan aturan yang berbeda, peluang termodinamika suatu


keadaan makro Wk , akan berbeda pula cara perhitungannya. Pada bagian
statistik Fermi-Dirac, Bose-Einstein, dan Maxwell-Boltzmann, akan diperli-
hatkan bagaimana menghitung Wk untuk ketiga kasus tersebut.

13.3 Observabel dan rata-rata bilangan okupasi

Sifat atau properti suatu observabel suatu sistem makroskopik bergantung pada
nilai rata-rata terhadap waktu dari properti atau sifat mikroskopik sistem terse-
but. Sebagai contoh, tekanan suatu gas bergantung pada nilai rata-rata ter-
hadap waktu dari transpor momentum pada suatu luasan. Melalui postulat
fundmental yang telah dibahas pada bagian sebelumnya, properti observabel
suatu sistem makroskopik akan pula bergantung pada nilai rata-rata properti
mikroskopik dari banyak replika suatu asembli yang diamati hanya pada suatu
waktu.

Kemudian tujuan dari teori statistik adalah mencari bagaimana menurunkan


ekspresi jumlah rata-rata dari partikel N j yang menempati tingkat energi j
13.3. OBSERVABEL DAN RATA-RATA BILANGAN OKUPASI 69

yang diperbolehkan dalam suatu asembli. Ekspresi yang akan diturunkan ini
disebut sebagai rata-rata bilangan okupasi pada tingkat (energi) j.

Misalkan Njk adalah bilangan okupasi tingkat j dalam keadaan makro k. Nila
g
rata-rata kelompok (grup) bilangan okupasi pada tingkat j, N j , diperoleh den-
gan mengalikan Njk dengan jumlah replika pada keadaan makro k, Wk N
dan dijumlahkan untuk seluruh keadaan makro dalam asembli, dibagi dengan
jumlah replika N , yaitu

g 1 X
Nj = Njk Wk N . (13.2)
N
k

Akan tetapi

X
N = Wk N , (13.3)
k

di mana N sama untuk semua keadaan makro sehingga

P
g k Njk Wk 1X
Nj = P = Njk Wk , (13.4)
k Wk
k

di mana rumusan untuk menghitung diperoleh dari Persamaan (13.1).

Dengan cara yang serupa dapat dicari rata-rata waktu dari bilangan okupasi
pada tingkat (energi) j. Sebagaimana telah dijelaskan dalam postulat funda-
mental termodinamika statistik bahwa semua keadaan mikro memiliki peluang
yang sama untuk muncul, yang artinya bahwa apabila sistem diamati untuk su-
atu rentang waktu yang lama t maka setiap keadaan mikro akan muncul dalam
rentang waktu total t yang sama. Total durasi waktu suatu asembli berada
pada keadaan makro k tak lain adalah perkalian dari rentang waktu t dengan
jumlah keadaan mikro Wk dalam keadaan makro tersebut. Jumlah dari semua
hasi perkalian ini untuk seluruh keadaan makro adalah sama dengan total waktu
pengamatan t,

X
t= Wk t. (13.5)
k

t
Kemudian nilai rata-rata waktu bilangan okupasi pada tingkat j, N j , diperoleh
dengan mengalikan bilakang okupasi pada tingkat j pada keadaan makro k, Njk
dengan waktu asembli tersebut pada keadaan makro k, Wk t, dijumlahkan un-
tuk seluruh keadaan makro dalam asembli tersebut, dan hasilnya dibagi dengan
total waktu pengamatan t, yaitu

t 1X
Nj = Njk Wk t. (13.6)
t
k
70 CATATAN 13. PELUANG TERMODINAMIKA

Dengan mempergunakan Persamaan (13.5) dan postulat bahwa t sama untuk


semua keadaan mikro, maka Persamaan (13.6) dapat dituliskan kembali menjadi

P
t k Njk Wk 1 X
Nj = = Njk Wk . (13.7)
Wk
k

Jadi apabila semua keadaan mikro memiliki peluang yang sama untuk muncul
maka rata-rata kelompok bilangan okupasi pada tingkat j sama dengan rata-
rata waktu bilangan okupasi pada tingkat j,

t g
Nj = Nj. (13.8)

seperti telah ditunjukkan dalam Persamaan (13.4) dan (13.7). Selanjutnya ke-
dua besaran yang sama ini akan dirujuk sebagai rata-rata bilangan okupasi pada
tingkat j, yaitu N j .
Catatan 14

Pengali dan

Pada bab sebelumnya telah ditunjukkan bagaimana bentuk peluang termodi-


namik Wk setiap keadaan makro k untuk ketiga jenis statistik, yaitu statistik
Maxwell-Boltzmann (MB), statistik Bose-Einstein (BE), dan statistik Fermi-
Dirac (FD).

14.1 Peluang termodinamik suatu keadaan


makro

Bila terdapat sejumlah tingkat energi j yang memiliki energi j dengan jumlah
keadaan energi atau degenerasi pada masing-masing tingkat energi adalah gj ,
maka untuk statistik MB bentuk peluang termodinamik suatu keadaan makro-
nya adalah

Y gjNj
WMB = N ! , (14.1)
j
Nj !

untuk statistik BE adalah

Y (gj + Nj 1)!
WBE = , (14.2)
j
(gj 1)!Nj !

dan untuk statistik FD adalah

Y gj !
WFD = . (14.3)
j
(gj Nj )!Nj !

71
72 CATATAN 14. PENGALI DAN

14.2 Keadaan makro yang paling mungkin

Dengan menggunakan pengali tak tentu Lagrange dan untuk mencari


keadaan makro yang memiliki keadaan mikro yang paling besar, digunakan
hubungan

d ln W + dN + dE = 0, (14.4)

dengan syarat

X X
N= Nj dN = dNj = 0, (14.5)
j j
X X
E= j Nj E = j dNj = 0. (14.6)
j j

Selanjutnya dapat diperoleh keadaan makro yang paling mungkin dari ketiga
statistik, atau disebut distribusi dari statistik tersebut. Distribusi MB memiliki
bentuk

gj
NjMB = , (14.7)
e(+j )

distribusi BE memiliki bentuk

gj
NjBE = , (14.8)
e(+j ) 1

dan distribusi FD memiliki bentuk

gj
NjFD = . (14.9)
e(+j ) + 1

14.3 Fungsi distribusi dalam bentuk diferensial

Setelah konstanta pengali dan diintepretasikan secara fisis dan diterap-


kan pada gas, masing-masing distribusi dari masing-masing statistik dapat dit-
uliskan dalam bentuk diferensial, yaitu untuk distribusi MB menjadi,

  1/2
2N d
NMB ()d = 3/2
(14.10)
(kT ) e/kT

distribusi BE menjadi
14.4. PENGALI 73

2(2m)3/2 V 1/2 d N h3
 
NBE ()d = 1 /kT , A = , (14.11)
h3 Ae 1 V (2mkT )3/2

dan distribusi FD menjadi

2/3
4(2m)3/2 V 1/2 d h2
  
3N
NFD ()d = , F (0) = . (14.12)
h3 e(F )/kT 1 2m 8V

14.4 Pengali

Terdapat berbagai kriteria untuk menentukan bagaimana arti sebenarnya dari


pengali . Dikarenakan jumlah partikel yang memiliki energi tak hingga harus-
lah nol maka Persamaan (14.7), (14.8), dan (14.9) memperkirakan bahwa nilai
haruslah lebih kecil dari nol dengan syarat dalam bagian kanan Persamaan
(14.5) dan (14.6).

Pendekatan dengan salah satu sudut pandang termodinamika dapat men-


gungkapkan bagaimana sifat dari pengali . Untuk itu dimisalkan terdapat
dua buah sistem, yang masing-masing tersusun atas N dan N partikel, yang
saling kontak sehingga dapat bertukar energi tetapi tidak bertukar partikel atau
dua buah sistem yang memenuhi kondisi

dN = 0, dN = 0, dE = 0. (14.13)

Dengan demikian energi total kedua sistem tak lain adalah

X X
E= j Nj + j Nj . (14.14)
j j

Selanjutnya kondisi dalam Persamaan (14.5) dan (14.6) akan menjadi

X X
dN = dNj = 0, dN = dNj = 0, (14.15)
j j

dan

X X
dE = j dNj + j dNj = 0. (14.16)
j j

Dalam bab sebelumnya telah ditunjukkan bahwa peluang termodinamik suatu


keadaan makro sistem gabungan tak lain adalah perkalian peluang termodi-
namik suatu keadaan makro dari masing-masing sistem, yaitu
74 CATATAN 14. PENGALI DAN

WT = W W . (14.17)

Dengan kembali menggunakan pengali tak tentu Lagrange, yang dalam hal ini
menjadi , , dan , maka diperoleh

d ln WT + dN + dN + dE = 0. (14.18)

Dikarenakan W hanya bergantung dari nj dan juga W hanya bergantung dari


nj maka dapat diperoleh bahwa

ln W
+ + j = 0 (14.19)
Nj

dan

ln W
+ + j = 0. (14.20)
Nj

Persamaan (14.19) dan (14.20) mendefinisikan suatu keadaan makro yang paling
mungkin muncul bagi kedua sistem penyusun sistem gabungan dan terlihat
bahwa kedua keduanya bergantung dari pengali . Dari kedua sistem hanya
terdapat satu parameter fisis yang perlu bernilai sama, karena keduanya kontak
secara termal, yaitu temperatur sesuai dengan hukum kenol termodinamika.
Dengan demikian dapat disimpulkan bahwa hanya fungsi dari temperatur

= (T ). (14.21)

Pengali dapat pula dilihat dari sudut pandang lain apabila dikaitkan dengan
perubahan energi dE. Untuk itu misalkan dalam suatu sistem diasupkan panas
sebesar dQ sehingga sebagian energi tersebut digunakan untuk melakukan kerja
dalam bentuk ekspansi dV . Menurut hukum pertama termodinamika

dE = dQ pdV, (14.22)

di mana dalam hal ini

X X X
dE = d j Nj = j dNj + Nj dj . (14.23)
j j j

Suku pertama pada ruas paling kanan Persamaan (14.23) menyatakan kerja
yang dilakukan sistem, di mana perubahan volume akan mengubah tingkat-
tingkat energi sistem. Dengan sendirinya dj pada tingkat energi j akan
14.4. PENGALI 75

berubah. Sedangkan suku kedua terkait dari perubahan susunan partikel pada
masing-masing tingkat energi dan hal ini dapat terjadi karena adanya asupan
panas. Perbandingan Persamaan (14.23) dengan Persamaan (14.22) akan mem-
berikan

X
Nj dj = pdV (14.24)
j

dan

X
j dNj = dQ. (14.25)
j

Saat kondisi kesetimbangan tercapai di mana tidak lagi terdapat perubahan


volume, substitusi Persamaan (14.25) ke dalam Persamaan (14.4) melalui Per-
samaan (14.23) akan memberikan

d ln W + dN + dQ = 0. (14.26)

Dengan menerapkan syarat bahwa jumlah partikel dalam sistem adalah tetap
akan diperoleh bahwa

d ln W = dQ. (14.27)

Selanjutnya dengan menggunakan hubungan bahwa

S = k ln , (14.28)

dan untuk sistem dengan jumlah partikel banyak sehingga

W, (14.29)

serta

dQ
dS = , (14.30)
T

maka dapat diperoleh bahwa

1
= . (14.31)
kT
76 CATATAN 14. PENGALI DAN

14.5 Ruang fasa enam dimensi

Suatu elemen ruang fasa enam dimensi d didefinisikan melalui relasi (Pointon,
1967) dalam bentuk

d = dxdydzdpx dpy dpz . (14.32)

Dalam Persamaan (14.32) terdapat elemen volume dalam ruang momentum


dan elemen volume dalam ruang koordinat. Pertama-tama, misalkan bahwa
volume dalam ruang momentum terletak antara dua nilai momentum, yaitu p
dan p + dp. Bila momentum total dinyatakan dalam koordinat polar (p, , )
maka elemen dari ruang momentum dengan rentang koordinat antara p sampai
p + dp, antara sampai + d, dan antara sampai + d adalah

dVp = (dp)(pd)(p sin d) = p2 sin dddp. (14.33)

Dengan cara yang sama apabila posisi terletak rentang koordinat antara x sam-
pai x + dx, antara y sampai y + dy, dan antara z sampai z + dz, maka elemen
ruang koordinat tak lain adalah

dV = (dx)(dx)(dz) = dxdydz. (14.34)

Dengan demikian volume dalam ruang fasa yang berkorespondensi dengan mo-
mentum dalam rentang koordinat antara p sampai p + dp, antara sampai
+ d, dan antara sampai + d dan posisi dalam rentang koordinat antara
x sampai x + dx, antara y sampai y + dy, dan antara z sampai z + dz adalah

d = dxdydz p2 sin dddp. (14.35)

Volume ruang momentum Vp yang terletak antara p dan p + dp dan tidak


lagi bergantung arah diperoleh dengan melakukan inegrasi Persamaan (14.33)
terhadap seluruh nilai dan , yaitu

Z Z 2
Vp = p2 dp sin d d = 4p2 dp, (14.36)
0 0

yang tak lain adalah volume dari kulit bola antara p dan p + dp.

Dengan demikian volume dalam ruang fasa yang terkait dengan Vp dalam
ruang momentum dan volume V dalam ruang koordinat diberikan oleh

Z
= 4p2 dp dxdydz = 4p2 dp V. (14.37)
V
14.6. DEGENERASI DALAM VOLUME RUANG FASA 77

Dengan menggunakan hubungan antara momentum dan kecepatan (pi =


mvi , i = x, y, z), elemen volume dalam ruang fasa untuk rentang kecepatan
antara vx sampai vx + dvx , antara vy sampai vy + dvy , dan antara vz sampai
vz + dvz dapat dituliskan dalam bentuk

d = dxdydz m3 dvx dvy dvz = m3 dvx dvy dvz V. (14.38)

Baik dengan menggunakan hubungan p = mv dan dp = mdv dalam Persamaan


(14.37) atau dvx dvy dvz = 4v 2 dv dalam Persamaan (14.38) dapat diperoleh
hubungan

= m3 4v 2 dv V. (14.39)

Selanjutnya adalah bagaimana mendifinisikan elemen ruang fasa dalam rentang


energi kinetik antara sampai + d. Dengan
menggunakan hubungan
p antara
momentum dan energi kinetik melalui p = 2m sehingga dp = m/(2)d,
diperoleh

p
= 4 2m m/(2) V = 2(2m)3/2 1/2 d V. (14.40)

14.6 Degenerasi dalam volume ruang fasa

Degenerasi atau jumlah keadaan energi pada suatu tingkat energi j , yaitu gj
dapat diungkapkan sebagai fungsi dari j , di mana umumnya suatu tingkat
energi memiliki energi antara j sampai j + dj . Dengan menggunakan asumsi
bahwa volume ruang fasa yang sama akan memberikan jumlah keadaan energi,
yang diperbolehkan, yang sama pula. Asumsi ini dapat dijustifikasi dalam kasus
mekanika kuantum, misalknya pada contoh partikel dalam kotak. Bila terdapat
B keadaan energi tiap satuan volume ruang fasa sehingga sebuah elemen ruang
fasa d akan mengandung Bd keadaan energi. Degenerasi dari tingkat energi
j tak lain adalah

gj = B()j , (14.41)

dengan ()j adalah volume dari ruang fasa enam dimensi yang terletak dalam
rentan energi antara j sampai j + dj dan dalam volume koordinat V dalam
sistem.

Dengan menggunakan Persamaan (14.40) untuk tingkat energi j Persamaan


(14.41) dapat dituliskan menjadi

1/2
gj BV 2(2m)3/2 j dj (14.42)
78 CATATAN 14. PENGALI DAN

14.7 Teori kinetik gas dan

Teori kinetik gas menyatakan bahwa energi rata-rata tiap partikel gas adalah

3
= kT, (14.43)
2

dengan k = R/NA adalah konstanta Boltzmann, R adalah konstanta gas uni-


versal, dan NP
A adalah bilangan Avogadro. Dengan menggunakan
P energi sistem
adalah E = j Nj j dan jumlah partikel adalah N = j Nj , maka dengan
menggunakan Persamaan (14.7) dapat dituliskan bahwa

X
E= (gj e+j )(j ) (14.44)
j

dan

X
N= (gj e+j ). (14.45)
j

Dengan demikian rata-rata energi tiap partikel adalah

+j j
P P
E j j g j e j j g j e
= = P +j
= P j
. (14.46)
N j gj e j gj e

Kemudian dengan menggunakan rumusan untuk gj sebagai fungsi dari j dalam


Persamaan (14.41) maka dapat diperoleh bahwa

R 1/2 R 3/2
0
j (BV 2(2m)3/2 j dj )ej 0
j ej dj
= R 1/2
= . R 1/2 (14.47)
0 (BV 2(2m)3/2 j dj )ej 0 j e
j d
j

dengan mengganti penjumlahan menjadi integral terhadap semua nilai energi


yang mungkin. Kemudian dengan mengingat bahwa kuantitas adalah lebih
kecil dari nol, maka integral parsial akan memberikan


3
Z Z
3/2 1/2
j ej dj = j ej dj (14.48)
0 2 0

sehingga Persamaan (14.47) akan menjadi

3
= . (14.49)
2
14.8. PENGALI 79

Dengan membandingkan Persamaan (14.49) dengan Persamaan (14.43) dapat


diperoleh bahwa

1
=
kT

seperti dalam Persamaan(14.31).

14.8 Pengali

Secara umum, pengali tidak bernilai sama untuk ketiga statistik, melainkan
bergantung dari pada kasus yang ditinjau. Berikut ini akan dibahas dengan
menggunakan statistik MB. Bila dituliskan bahwa

A = e (14.50)

maka Persamaan (14.7) akan menjadi

Nj = Agj ej (14.51)

dan jumlah total partikel adalah

X
N =A gj ej (14.52)
j

sehingga

N
A= P j
. (14.53)
j je
g

Dengan menggunakan hasil dari Persamaan (14.41) maka

N
A = R 1/2
, (14.54)
0 (BV 2(2m)3/2 j dj )ej

yang kembali diperoleh dengan mengganti penjumlahan dengan integrasi untuk


semua nilai energi yang mungkin. Kemudian dengan menggunakan integral
fungsi , yaitu
80 CATATAN 14. PENGALI DAN

Z
(n + 1) = ex xn dx, (14.55)
0
(n + 1) = n(n), (14.56)
(n + 1) = n!, (14.57)

(1/2) = , (14.58)

dan dengan x , maka

Z Z
1/2
j ej dj = () x1/2 ex dx
1/2
0 0

1/2 1/2
= () (3/2) = () . (14.59)
2

Dengan demikian

N N
A= = (14.60)
BV (2m/)3/2 BV (2mkT )3/2

dan

 
N
= ln A = ln . (14.61)
BV (2mkT )3/2

Untuk gas BE akan diperoleh bahwa

N h3
A= (14.62)
V (2mkT )3/2

atau

N h3
 
= ln A = ln , (14.63)
V (2mkT )3/2

di mana dalam hal ini degenerasi, atau lebih tepatnya jumlah keadaan energi
yang menempati elemen ruang fasa enam dimensi didefinisikan sebagai

d
g= . (14.64)
h3

Sedangkan dalam statistik FD maka akan terkait dengan apa yang dikenal
sebagai energi Fermi F , yaitu lewat hubungan
14.8. PENGALI 81

F
= . (14.65)
kT

Gambar 14.1: Ilustrasi band gap dan energi Fermi pada semikonduktor.
82 CATATAN 14. PENGALI DAN
Catatan 15

Energi Bebas Helmholtz

Dengan mengetahui temperatur dan entropi suatu sistem dalam deskripsi statis-
tiknya ada gunanya pula untuk mengaitkan fungsi termodinamika lain dengan
sifat-sifat statistik. Dalam hal ini, misalnya saja adala energi bebas Helmholtz
F sistem yang didefinisikan sebagai

15.1 Energy bebas Helmholtz

F = E T S. (15.1)

Bila suatu sistem mengalami perubahan kecil keadaannya pada temperatur


tetap sehingga energinya berubah dari E menjadi E + dE dan entropinya
berubah dari S menjadi S + dS, maka perubahan energi bebas Helmholtznya
menjadi

dF = dE T dS. (15.2)

Dari hukum pertama dan kedua termodinamika dapat dituliskan bahwa

T ds dE + dW, (15.3)

di mana dW adala kerja yang dilakukan oleh sistem terhadap lingkungannya


dan tanda sama dengan berlaku hanya untuk proses termodinamika reversibel.
Dengan mengabungkan Persamaan (15.2) dan (15.3) dapat diperoleh ketidak-
samaan

dF dW. (15.4)

83
84 CATATAN 15. ENERGI BEBAS HELMHOLTZ

Selama terjadinya perubahan, energi bebas F akan berkurang sejumlah sama


atau lebih besar dari kerja yang dilakukan oleh sistem. Bila tidak ada kerja yang
dilakukan (dW = 0) setiap perubahan isotermal dalam energi bebas akan kurang
dari atau sama dengan nol. Keadaan setimbang sistem dalam kondisi ini berada
pada keadaan di mana energi bebas telah berkurang sampai nilai minimumnya
karena perubahan parameter dari sistem akan memberikan perubahan energi
bebas sebesar

dF = 0 (15.5)

dan tidak ada kerja yang dilakukan sistem.

Saat temperatur sistem bernilai tetap, penerapan Persamaan (15.5) akan mem-
perbolehkan keadaan kesetimbangan sistem ditentukan bila energi bebas dike-
tahui bentuknya dalam berbagai parameter termodinamika.

15.2 Ekspansi reversibel

Kegunaan dari diketahuinya energi bebas terkait pula dengan hubungannya den-
gan fungsi termodinamika lainnya dari sistem. Sebagai contoh, misalnya sebuah
berubahan reversibel sistem berlangsung dengan perubahan dT dan kerja yang
dilakukan hanya merupakan ekspansi sehingga volumenya bertambah sebesar
dV . Kerja yang dilakukan sistem adalah pdV dengan p adalah tekanan sis-
tem. Pertidaksamaan (15.3) akan menjadi persamaan untuk proses reversible
sehingga

T ds = dE + dW. (15.6)

Dengan nilai T dS ini dan bahwa temperatu boleh berubah maka perubahan
energi bebas F dapat diperoleh dari Persamaan (15.1), yaitu

dF = dE T dS SdT = dE (dE + pdV ) SdT = pdV SdT. (15.7)

Dari Persamaan (15.7) dapat diperoleh bahwa

 
F
p= (15.8)
V T

dan

 
F
S= . (15.9)
T V
15.3. ENERGI SEBAGAI FUNGSI DARI ENERGI BEBAS 85

Persamaan (15.8) akan berguna saat mencari persamaan keadaan sistem yang
akan memberikan tekanan sistem dalam fungsi volume dan temperatur sistem.
Salah satu contoh persamaan keadaan misalnya untuk gas ideal adalah pV =
N RT .

15.3 Energi sebagai fungsi dari energi bebas

Substitusi Persamaan (15.9) ke dalam Persamaan (15.1) akan memberikan

 
F
F =E+T . (15.10)
T V

Selanjutnya dapat diperoleh bahwa

  

F (F/T )
2
E =F T = T (15.11)
T V T
 V
(F )
, (15.12)
V

dengan = 1/kT (di mana = T /kT 2 ).

15.4 CV dari E

Dengan menggunakan definisi dari CV

 
E
CV = (15.13)
T V

dapat diperoleh

2F
 
CV = T (15.14)
T 2 V

atau

2 (F )
 
2
CV = k . (15.15)
2 V

Dengan menggunakan pengetahuan mengenai distribusi statistik sistem, energi


bebas dapat ditentukan dalam representasi statis tik.
86 CATATAN 15. ENERGI BEBAS HELMHOLTZ
Catatan 16

Fungsi Partisi Boltzmann

Peluang termodinamika untuk keadaan makro yang paling sering muncul Wmax
pada sistem gas sempurna memiliki bentuk

 
N E
ln Wmax = N ln + . (16.1)
A kT

Suatu kuantitas Z, yang nantinya akan disebut sebagai fungsi partisi, diperke-
nalkan sebagai pengganti suku N/A.

16.1 Fungsi partisi Boltzmann

Kuantitas Z didefinisikan sebagai

N
Z= . (16.2)
A

Dengan menggunakan Persamaan (14.7) dan (14.50) yang disubtsitusikan ke


Persamaan (16.2) akan diperoleh

gj e(j /kT )
P
j
X
Z= = gj ej /kT (16.3)
e j

yang disebut sebagai fungsi partisi Boltzmann atau fungsi partisi, untuk suatu
partikel dalam suatu sistem. Istilah ini digunakan karena dalam ekspresi Z,
suku-suku dalam penjumlahan di atas menentukan bagaimana partikel-partikel
dalam sistem terdistribusi atau terpartisi di antara berbagai tingkat-tingkat
energi.

87
88 CATATAN 16. FUNGSI PARTISI BOLTZMANN

Selain dalam bentuk Persamaan (16.3) fungsi partisi dapat pula diungkapkan
untuk keadaan-keadaan energi yang ada secara individual. Bila energi dari
keadaaan energi i adalah i , dan karena degenerasi dari suatu keadaan energi
adalah 1, maka

X
Z= ei /kT . (16.4)
i

(Dikarenakan terdapat kemungkinan untuk menuliskan Z dalam bentuk ini, Z


kadang dirujuk sebagai jumlah meliputi seluruh keadaan energi atau jumlah
keadaan energi bagi suatu sistem. Lambang Z diambil dari ekspresi ekivalen
dalam bahasa Jerman Zustandsumme.)

Fungsi partisi yang diperoleh baik dari Persamaan (16.3) maupun Persamaan
(16.4) adalah bukan besaran termodinamika yang terukur atau dapat diukur
secara umum, ataupun muncul dalam persamaan termodinamika yang wajar.
Tetapi, hasil tersebut merupakan suatu jembatan yang penting antara ekspresi
statistik untuk suatu keadaan energi suatu sistem dengan fungsi-fungsi termod-
inamika terkait, yang akan ditunjukkan kemudian.

16.2 Fungsi partisi dan energi bebas Helmholtz

Dengan menggunakan definisi dari fungsi partisi dalam Persamaan (16.2) ke


dalam Persamaan (16.1), asumsi Wmax , dan kaitan antara entropi S dan
maka

E
S = N k ln Z + . (16.5)
T

(Boltzmann mendefinisikan bahwa S = k ln Wmax sedangkan Planck menggu-


nakan S = k ln . Untuk sistem dengan jumlah keadaan mikro yang banyak
kedua ungkapan tersebut hampir tidak memiliki perbedaan, akan tetapi apa-
bila terdapat pembedaan maka ekspresi yang diungkapkan oleh Planck adalah
yang lebih tepat.)

Dengan menggunakan Persamaan (15.1) maka dapat diperoleh bahwa

F = N KT ln Z. (16.6)

Dengan demikian besaran-besaran lain seperti energi E, tekanan p, dan capa-


sitas panas CV dapat diperoleh.
16.3. ENERGI SISTEM DAN FUNGSI PARTISI 89

16.3 Energi sistem dan fungsi partisi

Dengan menggunakan Persamaan (15.11) atau melalui energi rata-rata dapat


diperoleh kaitan antara E dan Z. Energi rata-rata setiap partikel diperoleh
lewat

P
E j nj j
= = P
N j nj
j /kT j /kT
P P
j g j j e j g j j e
= P /kT
= P j /kT
j gj e j gj e
j

j /kT
P
j g j j e
= .
Z

Dengan menggunakan Persamaan (16.3) dapat diturunkan bahwa

 
Z 1 X
= gj j ej /kT (16.7)
T V kT 2 j

sehingga

Z
kT 2
  
T V ln Z
= = kT 2 . (16.8)
Z T V

Kemudian dengan menggunakan E = N dapat diperoleh bahwa

   
2 ln Z ln Z
E = N kT =N (16.9)
T V V

16.4 Entropi dan fungsi partisi

Dengan melakukan substitusi Persamaan (16.9) ke dalam Persamaan (16.5) da-


pat diperoleh

   
ln Z
S = N k ln Z + (16.10)
ln V

yang menunjukkan bahwa entropi S dapat dihitung hanya dengan meman-


faatkan informasi mengenai fungsi partisi Z.
90 CATATAN 16. FUNGSI PARTISI BOLTZMANN

16.5 Energi bebas tiap partikel

Persamaan (16.6) dapat diubah bentuknya menjadi

Z = eF/N kT (16.11)

di mana dengan f = F/N yang menyatakan energi bebas tiap partikel maka
fungsi partisi untuk tiap keadaan energi bagi partikel dalam sistem dapat dit-
uliskan menjadi

X
Z= ei /kT = ef /kT . (16.12)
i

16.6 Kapasitas panas spesifik pada volume tetap

Dengan menggunakan definisi dari Persamaan (15.13) dapat diperoleh bahwa

2 ln Z
   
E ln Z
CV = = N K 2T + T2 (16.13)
T V T T 2

16.7 Tekanan sistem

Tekanan dapat diperoleh dari Persamaan (15.8) dan Persamaan (16.6)

   
F ln Z
p= = N kT (16.14)
V T V T
Catatan 17

Gas Ideal Monoatomik

Suatu gas ideal yang terdiri dari N molekul identik yang masing-masing
bermassa m. Molekul-molekul gas tak-terbedakan dan jumlah molekul dalam
tiap keadaan energi yang mungkin, kecuali pada temperatur amat rendah se-
hingga semua gas mencair, adalah amat kecil. Statistik yang cocok adalah
statistik klasik.

Langkah pertama adalah menghitung fungsi partisi dari sistem ini, sebagaimana
telah diungkapkan dalam Persamaan (16.3), yaitu

X
Z= gj ej /kT
j

17.1 Tingkat energi makro

Untuk memperoleh fungsi partisi diperlukan informasi mengenai energi j


dan degenerasi gj pada tiap tingkat energi j. Di sini diasumsikan bahwan
molekul-molekul tidak berinteraksi kecuali saat saling bertumbukan secara in-
stan, sehingga tiap molekul dapat dianggap sebagai partikel bebas dan memiliki
susunan tingkat energi yang sama sebagaimana sebuah partikel dalam kotak.
Dengan menggunakan mekanika kuantum dapat diperoleh bahwa tingkat energi
partikel dalam kotak adalah

n2j h2 V 2/3
j = , (17.1)
8m

di mana n2j = n2x + n2y + n2z , dan nx , ny , nz adalah bilangan bulat yang masing-
masing dapat bernilai 1, 2, 3, ....

Degenerasi masing-masing tingkat energi atau jumlah keadaan energi dalam


tiap tingkat energi dapat dengan mudah dihitung apabila bilangan kuantum

91
92 CATATAN 17. GAS IDEAL MONOATOMIK

bernilai kecil. Dalam banyak contoh, sayangnya, tingkat-tingkat energi suatu


asembli berjarak amat rapat dibandingkan dengan nilai energi tingkat energi itu
sendiri. Dengan demikian tingkat-tingkat energi dapat dikelompokkan menjadi
suatu grup dengan lebar energi j untuk tingkat-tingkat energi yang memiliki
yang energi antara j dan j + j . Setiap grup ini dirujuk sebagai tingkat
energi makro.

Nyatakan bahwa Gj merepresentasikan jumlah total keadaan energi yang


mungkin dalam tingkat energi sampai dan termasuk yang memiliki energi j .
Jumlah keadaaan energi yang mungkin Gj dalam tingkat energi makro terse-
but adalah sama dengan jumlah keadaan-keadaan energi dalam semua tingkat
energi yang termasuk ke dalam tingkat energi makro tersebut. Jadi, Gj tak
lain adalah degenerasi dari tingkat energi makro, di mana degenerasi ini muncul
karena terjadinya pengelompokan sejumlah besar tingkat-tingkat energi, akan
tetapi gj ditetapkan dari sifat alami asembli.

17.2 Degenerasi tingkat energi makro

Bayangkan bilangan kuantum nx , ny , nz merupakan label tiga buah sumbu


yang saling tegak lurus. Setiap kombinasi tiga angka nx , ny , nz akan mendefin-
isikan sebuah titik dalam ruang-n. Setiap titik seperti itu terkait dengan suatu
keadaan energi yang mungkin, memberikan nilai bilangan kuantum positif. Da-
pat diambil bahwa setiap titik terdapat dalam kubus dengan panjang rusuk
dalam satuan panjang sehingga volumenya adalah satu.

Bilangan kuantum nj berkaitan dengan sebuah vektor dalam ruang-n yang di-
hitung dari pusat koordinat ke setiap titik, karena n2j = n2x + n2y + n2z . Dalam
sistem, dengan volume yang diberikan, energi hanya bergantun dari nj , sehingga
semua kedaan energi dengan energi yang saya berada pada permukaan bola den-
gan jari-jari nj dari pusat koordinat. Karena nx , ny , nz semua bernilai positif,
dan terdapat satu titik setiap satuan volume dalam ruang-n, maka jumlah to-
tal Gj dari semua keadaan energi yang mungkin, dalam semua tingkat energi
sampai dan termasuk yang memiliki energi j , adalah sama dengan volume satu
oktan dari bola dengan jari-jari nj , yaitu

1 4 3
Gj = n = n3j . (17.2)
8 3 j 6

Kulit bola sudah tentu akan memotong beberapa sel satuan sehingga tidak jelas
apakah titik yang merepresentasikan suatu keadaan energi terdapat di dalam
atau di luar permukaan bola. Akan tetapi, untuk nilai nj besar, yang meru-
pakan kasus kebanyakan molekul-molekul suatu gas pada temperatur umum,
ketidakpastian ini menjadi dapat diabaikan karena kecil.

Selanjutnya, jumlah keadaan energi dalam tingkat keadaan makro yang memiliki
energi antara j dan j + j , atau degenerasi tingkat energi makro Gj adalah
17.3. FUNGSI PARTISI SISTEM 93


Gj = 3n2j nj = n2j nj . (17.3)
6 2

Secara geometri, nilai ini terkait dengan jumlah titik-titik dalam suatu kulit bola
dengan jari-jari nj dan tebal nj . Dengan melihat bentuk dalam Persamaan
(17.3), degenerasi akan meningkat dengan pertambahan nj secara kuadratik
untuk harga nj yang sama.

17.3 Fungsi partisi sistem

Dengan demikian fungsi partisi sistem dapat dituliskan sebagai

X
Z= Gj ej /kT , (17.4)
j

di mana dengan menggunakan Persamaan (17.3) dan (17.1) akan menjadi

  2 2/3  
X 2 h V
Z= n exp n2j nj . (17.5)
2 j j 8mkT

Dengan

  2 2/3  
h V
f (nj ) = n2j exp n2j (17.6)
8mkT

maka f (nj )nj tak lain adalah luas di bawah kurva dari grafik f (nj ) terhadap
nj . Nilai Z terkait dengan jumlah seluruh luas seperti itu untuk seluruh nilai
nj dari j = 1 sampai j = karena tidak terdapat batas atas untuk nilai yang
diperbolehkan bagi nj .

Untuk suatu aproksimasi yang baik, jumlah dalam Persamaan (17.5) sama den-
gan luas di bawah kurva kontinu yang melewati puncak nilai-nilai f (nj ), antara
batas 0 dan adalah

  2 2/3  
h V
Z
Z= n2j exp n2j dnj . (17.7)
2 0 8mkT

Dengan menggunakan tabel integral tentu dapat diperoleh bahwa

 3/2
2mkT
Z=V . (17.8)
h2
94 CATATAN 17. GAS IDEAL MONOATOMIK

Hasil dalam Persamaan (17.8) ini akan dapat digunakan untuk mencari besaran-
besaran termodinamika. Atau lebih umum apabila digunakan

 
3 2mkT
ln Z = ln V + ln . (17.9)
2 h2

17.4 Persamaan keadaan dan besaran-besaran


termodinamika

Dengan menggunakan Persamaan (15.8) dan (16.6) dapat diperoleh bahwa

   
F
p= = (N kT ln Z)
V V
T  T
ln Z N kT nRT
= N kT = = , (17.10)
V T V V

yang tak lain adalah persamaan keadaan untuk gas ideal. Kemudian dengan
menggunakan Persamaan (15.11)

    
(F/T ) N kT ln Z
U = E = T 2 = T 2
T T T
  V  2
V

2 ln Z 2 3 h 2mk
= N kT = N kT
T V 2 2mkT h2
3 3
= N kT = nRT, (17.11)
2 2

yang juga cocok dengan hasil teori kinetik untuk gas monoatomik dengan tiga
derajat kebebasan.

Kapasitas panas pada volume tetap pun dapat diperoleh melalui

 
U 3 3
CV = = N k = nR (17.12)
T V 2 2

dan juga kapasitas spesifik molal

CV 3
cV = = R. (17.13)
n 2
Catatan 18

Distribusi Laju Molekular

Sebagaimana penurunan dalam Bab Gas Ideal Monoatomik di sini, di mana


digunakan pula statistik klasik, tetap digunakan tingkat energi makro.

18.1 Bilangan okupasi rata-rata tingkat energi


makro

Distribusi akan dinyatakan dalam bilangan okupasi rata-rata suatu tingkat en-
ergi makro, yang meliputi suatu rentang energi antara j sampai j + j .

Variabel N menyatakan jumlah molekul yang memiliki energi sampai dan ter-
masuk berenergi j . Rata-rata jumlah molekul yang termasuk ke dalam tingkat
energi makro atau rata-rata bilangan okupasi tingkat energi makro tersebut
adalah Nj . Kuantitas Nj dan Gj terkait dengan bilangan okupasi rata-
rata suatu tingkat energi N j dan degenerasi suatu tingkat energi tunggal gj .
Untuk kedua fungsi distribusi, baik Maxwell-Boltzmann maupun klasik, memi-
liki bentuk

N
Nj = Gj e/kT . (18.1)
Z

Dengan menggunakan

n2j h2 V 2/3 1
j = = mvj2
8m 2
2 2
4m vj 2mv j
n2j = 2 2/3 nj = .
h V hV 1/3
2m
nj = vj .
hV 1/3

95
96 CATATAN 18. DISTRIBUSI LAJU MOLEKULAR
!
4m2 vj2 4m3 V 2

2m
Gj = n2j nj = vj = vj vj . (18.2)
2 2 h2 V 2/3 hV 1/3 h3

Dengan menghilangkan indeks j pada Persamaan (18.2) di atas dan memberikan


indeks v dapat dituliskan

4m3 V 2
Gv = v v. (18.3)
h3

Gunakan hasil dalam Persamaan (18.3) ke dalam Persamaan (18.1) dan hasil
untuk Z dari Persamaan (17.8) sehingga

N
Nj = Gj e/kT
Z
4m3 V 2

N 2
Nv = v v emv /2kT
V 2mkT 3/2
 h3
h2
4N  m 3/2 2 mv2 /2kT
= v e v. (18.4)
2kT

Kuantitas Nv merepresentasikan rata-rata total jumlah molekul yang memiliki


laju sampai dan termasuk v, dan Nv adalah rata-rata jumlah dengan laju
antara v dan v + v.

Persamaan (18.4) disebut sebagai fungsi distribusi laju Maxwell-Boltzmann.

18.2 Laju yang paling mungkin

Laju yang paling mungkin vm untuk molekul-molekul yang memiliki laju antara
v dan v+v dapat diperoleh dengan mencari nilai maksimum dari (18.4) dengan
hanya memperhatikan

d  2 mv2 /2kT 
v e =0
dv

yang akan memberikan hasil

r
2kT
vm = . (18.5)
m

Dengan demikian Persamaan (18.4) dapat dituliskan dalam bentuk yang lebih
kompak, dengan menggunakan vm , yaitu
18.3. LAJU RATA-RATA 97

Nv 4N 2 2
= 3 v 2 ev /vm . (18.6)
v vm

Fungsi distribusi laju ini bergantung terhadap temperatur T melalui kuantitas


vm .

18.3 Laju rata-rata

Laju rata-rata diperoleh melalui

1 X
v= vNv .
N

Dengan menggunakan Persamaan (18.6) dan melakukan aproksimasi jumlah


dengan integral dapat diperoleh bahwa


r
4 2 8kT
Z
3 v 2 /vm
2
v= 3 v e dv = vm = (18.7)
vm 0 m

18.4 Laju rms

Dengan menggunakan cara yang sama untuk mencari v dapat dicari vrms , yaitu

1/2


1 X 2
vrms = v= v Nv
N
 Z 1/2 r
4 4 v 2 /vm
2 3 3kT
= 3 v e dv = vm = (18.8)
vm 0 2 m

18.5 Perbandingan ketiga laju

Dengan melihat Persamaan (18.5), (18.7), dan (18.8) dapat dituliskan bahwa

r r r
2kT 8kT 3kT
vm : v : vrms = : : = 1 : 1.128 : 1.224. (18.9)
m m m
98 CATATAN 18. DISTRIBUSI LAJU MOLEKULAR

18.6 Ekipartisi energi

Bila energi partikel-partikel dalam suatu sistem terbuat dari suku-suku


kuadratik dalam koordinat posisi dan momentum partikel, maka untuk setiap
suku kuadratik, kontribusinya kepada energi rata-rata akan berbentuk 21 kT di
mana T adalah temperatur dari sistem.
Catatan 19

Paradoks Gibb

Bab ini akan menceritakan bagaimana entropi bertambah setelah dua jenis gas
yang semula terpisahkan oleh dinding diatermal dicampurkan. Pertambahan
entropi ini sayangnya juga terjadi saat kedua jenis gas merupakan jenis gas
yang sama, yang menjadikannya suatu paradoks yang dikenal sebagai paradoks
Gibb.

19.1 Fungsi partisi

Secara umum, degenerasi pada setiap tingkat energi j , yaitu gj , dapat di-
ungkapkan dalam bentuk kontinu di ruang- menjadi

d
gj (19.1)
h3

sehingga fungsi partisi dapat dituliskan menjadi

d
X Z
j /kT
Z= gj e e/kT . (19.2)
j
h3

Elemen ruang fasa enam-dimensi d terdiri dari tiga elemen ruang spasial dan
tiga elemen ruang momentum, yaitu

d = (dV )(dVp ) = (dxdydz)(dpx dpy dpz ). (19.3)

Suku pertama dalam ruas paling kanan pada Persamaan (19.2), yaitu e/kT
merupakan fungsi dari energi, dengan demikian elemen ruang fasa enam-dimensi
d dalam Persamaan (19.3) harus pula dirumuskan dalam variabel energi .

99
100 CATATAN 19. PARADOKS GIBB

Untuk partikel-partikel gas yang bebas bergerak tanpa pengaruh suatu medan
maka energinya hanya merupkana energi kinetik sehingga

p2
= . (19.4)
2m

Oleh karena itu elemen ruang momentum dVp = dpx dpy dpz dalam Persamaan
(19.3) harus dirumuskan dalam momentum total p dan bukan dalam px , py , dan
pz . Dengan mengingat bentuk bahwa dxdydz dalam sistem koordinat kartesian
memiliki representasi r2 dr sin dd dalam sistem koordinat bola maka dengan
cara yang sama dapat diperoleh

dVp = dpx dpy dpz = p2 dp sin dd. (19.5)

Dan karena energi dalam Persamaan (19.4) hanya bergantung dari p dan tidak
dan maka Persamaan (19.5) dapat diungkapkan dalam bentuk

Z Z 2
dVp = p2 dp sin d d = 4p2 dp. (19.6)
0 0

Ungkapan untuk dp dapat diperoleh dari d dengan menurunkan Persamaan


(19.4) terhadap sehingga diperoleh

r r
d p 2m 2 m
= = = dp = d. (19.7)
dp m m m 2

Dengan demikian Persamaan (19.6) akan menjadi

r 
2 m
dVp = 4p dp = 4(2m) d = 2(2m)3/2 1/2 d. (19.8)
2

Selanjutnya dengan menuliskan elemen volume dalam ruang spasial dapat di-
integralkan menjadi V karena fungsi energi yang dibahas tidak tergantung dari
koordinat spasial. Akhirnya elemen ruang fasa enam-dimensi d dapat dit-
uliskan menjadi

d = 2V (2m)3/2 1/2 d (19.9)

sehingga bentuk kontinu gj dalam Persamaan (19.1) dalam ruang energi akan
menjadi

gj g()d = 2V (2m)3/2 1/2 d g() = 2V (2m)3/2 1/2 . (19.10)


19.2. BEBERAPA BESARAN TERMODINAMIKA 101

Dengan menggunakan Persamaan (19.9) dan V maka Persamaan (19.2) dapat


dituliskan menjadi

d 2V (2m)3/2
Z Z
Z e/kT 3
= 1/2 e/kT d (19.11)
h h3 0

Persamaan (19.11) dapat diselesaikan dengan menggunakan integral fungsi (x)


sehingga



Z
1/2 e/kT d = (kT )3/2 (3/2) = (kT )3/2 (19.12)
0 2

Dengan menggunakan hasil dalam Persamaan (19.12) maka Persamaan (19.11)


akan menjadi


2V (2m)3/2 3/2 V (2mkT )3/2
Z= (kT ) = (19.13)
h3 2 h3

19.2 Beberapa besaran termodinamika

Energi bebas Helmholtz F dapat dihitung dengan menggunakan hubungan

F = N kT ln Z (19.14)

sehingga diperoleh

V (2mkT )3/2
 
F = N kT ln . (19.15)
h3

Entropi S dapat diperoleh lewat

 
F
S= (19.16)
T V

sehingga diperoleh

V (2mkT )3/2

  

S= N kT ln
T h3 V
V (2mkT )3/2
   

= N kT ln
T h3 V
102 CATATAN 19. PARADOKS GIBB

V (2mkT )3/2 h3 V (2mk)3/2


 
= N k ln 3
+ N kT
h V (2mkT ) 3/2 h3
3
T
2
V (2mkT )3/2
 
3
= N k ln + N k. (19.17)
h3 2

Energi total sistem U adalah

U = F + TS (19.18)

sehingga dengan menggunakan Persamaan (19.15) dan (19.17) dapat dituliskan


kembali menjadi

V (2mkT )3/2 V (2mkT )3/2


    
U = N kT ln + T N k ln
h3 h3

3
+ Nk
2
3
= N kT (19.19)
2

yang cocok dengan energi dalam suatu gas ideal monoatomik. Sedangkan
tekanan dapat diperoleh dari

 
F
p= (19.20)
V T

yang akan memberikan

V (2mkT )3/2
   
N kT
p= N kT ln 3
= . (19.21)
V h T V

Persamaan (19.21) akan memberikan persamaan keadaan untuk gas ideal klasik
monoatomik. Selanjutnya, kapasitas panas pada volume konstan CV dapat
diperoleh lewat

   
Q U
CV = = (19.22)
T V T V

sehingga dengan menggunakan Persamaan (19.19) akan diperoleh

3
CV = N k. (19.23)
2
19.3. PARADOKS GIBB 103

Sedangkan untuk kapasitas panas pada tekanan konstan Cp

     
Q U V
Cp = = +p . (19.24)
T p T p T p

Dengan menggunakan Persamaan (19.21) akan diperoleh bahwa

Cp = CV + N k. (19.25)

19.3 Paradoks Gibb

Misalkan terdapat gas berjenis sama yang menempati dua ruang berbeda dengan
volume V yang sama, temperatur T yang sama, tekanan p yang sama, jumlah
partikel N yang sama, sehingga akan memiliki entropi S yang sama dan energi
total U yang sama.

Apabila kedua ruang tersebut digabungkan maka besaran-besaran yang dise-


butkan sebelumnya akan menjadi

T = T, (19.26)

U = 2U, (19.27)
V = 2V, (19.28)
N = 2N, (19.29)
3/2
 
(2V )(2mkT ) 3
S = (2N )k ln 3
+ (2N )k
h 2
3/2
   
V (2mkT ) 3
= 2 N k ln + N k + 2N k ln 2
h3 2
= 2S + 2N k ln 2 > 2S. (19.30)

Persamaan (19.30) menceritakan bahwa dengan menghilangkan pembatas an-


tara kedua ruang yang berisi gas yang sama, akan terjadi peningkatan entropi
sebesar 2N k ln 2 dari seharusnya yang diperkirakan, yaitu 2S. Sebaliknya, bila
pembatas dipasang kembali maka entropi kedua sistem yang terpisah akan kem-
bali menjadi 2S. Seharusnya dinding pemisah tidak memiliki pengaruh pada
entropi gas apabila kedua jenis gas adalah sama. Kontradiksi kelakuan entropi
ini disebut sebagai paradoks Gibb.

Bila kedua jenis gas berbeda, tidak terdapat kontradiksi ini karena setelah dind-
ing pembatas dihilangkan dan entropi bertambah, pemasangan kembali dinding
pembatas tidak akan mengambalikan dinding pembatas tidak akan mengamba-
likan gas kepada kondisi masing-masing jenis gas terpisah. Masing-masing gas
104 CATATAN 19. PARADOKS GIBB

memiliki entropi awal yang berbeda karena memiliki massa partikel gas yang
bebeda.

Suku campuran apabila kedua gas berjenis sama ini dapat pula dipandang dari
sisi statistik klasik bahwa ruang yang baru akan membuat seakan-akan gj men-
jadi 2gj sehingga peluang termodinamika keadaan makro paling mungkin akan
menjadi


Wmax = (2N Wmax )(2N Wmax ). (19.31)

Dengan menggunakan hubungan entropi menurut Boltzmann

= k ln (2N Wmax )2
 
S = k ln Wmax

= 2k ln Wmax + 2kN ln 2
= 2S + 2kN ln 2. (19.32)

19.4 Gas ideal semi-klasik

Gibb kemudian mempostulatkan bahwa terdapat kesalahan dalam menghi-


tung jumlah keadaan sehingga peluang termodinamika untuk statistik Maxwell-
Boltzmann perlu dikoreksi sehingga menjadi N ! kali lebih kecil dari seharusnya,
yaitu dari

Y gjNj
WMB = N ! (19.33)
j
Nj !

menjadi

Y gjNj
Wsemiklasik = . (19.34)
j
Nj !

Yang disebut sebagai statistik semi-klasik. Koreksi ini akan memberikan kon-
sekuensi dalam perumusan entropi sehingga menjadi

V (2mkT )3/2
 
5
S = N k ln + N k. (19.35)
N h3 2

yang dikenal sebagai persamaan Sackur-Tetrode. Banding Persamaan (19.35)


dengan Persamaan (19.17) semula. Dapat pula dilihat bahwa volume yang
berperan di sini sekarang adalah volume untuk tiap partikel atau N/V sehingga
19.4. GAS IDEAL SEMI-KLASIK 105

apabila dua jenis gas dicampur dengan kondisi seperti telah disebabkan sebelum-
nya, volume untuk tiap partikel tidak akan berubah. Dengan Persamaan (19.35)
akan diperoleh bahwa

(2V )(2mkT )3/2


 
5
S = (2N )k ln 3
+ (2N )k
(2N )h 2
V (2mkT )3/2
   
5
= 2 N k ln + N k = 2S.
N h3 2

Persamaan (19.35) juga akan membuat perubahan hubungan antara besaran-


besaran termodinamika terhadap fungsi partisi. Akan tetapi hasi akhir besaran-
besaran termodinamika yang diturunkan tidak akan berubah. Hal ini akan
dibahas pada bab khusus mengenai gas semi-klasik.
106 CATATAN 19. PARADOKS GIBB
Catatan 20

Ekipartisi Energi

Bila energi partikel-partikel dalam suatu sistem berbentuk kuadrat dari koordi-
nat posisi dan momentum sistem maka setiap suku yang mengandung kuadrat
tersebut tersebut akan berkontribusi terhadap energi rata-rata sebesar 12 kT di
mana T adalah temperatur sistem. Hal ini akan dibahas sebagai suatu aplikasi
dari statistik Maxwell-Boltzmann.

20.1 Bentuk-bentuk energi

Energi suatu partikel dapat berbentuk murni energi kinetik, misalnya dalam
arah-x

p2x
x = , (20.1)
2m

yang juga berlaku sama dengan dalam arah-y dan arah-z. Dapat pula berbentuk
energi kinetik dan energi potensial, misalnya pada osilator harmonik, yang untuk
arah-x-nya adalah

p2x 1
x = + x2 , (20.2)
2m 2

yang akan memiliki bentuk yang sama pula dalam kedua arah lainnya.

Dengan demikian, secara lengkap suku energi merupakan fungsi dari x, y, z,


px , py , dan pz atau koordinat dari ruang fasa enam-dimensi . Suatu bentuk
lengkap yang bergantung kuadrat dari koordinat-koordinat ruang adalah

107
108 CATATAN 20. EKIPARTISI ENERGI

!
p2x p2y p2z
   
1 1 1
= + x2 + + y 2 + + z 2 . (20.3)
2m 2 2m 2 2m 2

20.2 Rata-rata energi kinetik

Dengan demikian, apabila ingin dihitung x untuk suatu partikel gas


monoatomik yang bebas dari pengaruh medan apapun, dapat dilakukan lewat

2
R /kT
(pRx /2m)e d
x = , (20.4)

e /kT d

dengan d = dxdydxdpx dpy dpz . Dengan melihat bentuk umum dalam Per-
samaan (20.3) maka apabila dituliskan

(p2x /2m) (20.5)

merupakan suatu suku yang tidak lagi bergantung dari px . Dengan menggu-
nakan cara ini maka Persamaan (20.7) dapat dituliskan menjadi

R [(p2 /2m)]/kT R 2


e x dV dpy dpz (p2x /2m)epx /2mkT dpx
x = R [(p2 /2m)]/kT R 2

e x dV dpy dpz
epx /2mkT dpx
R 2 2
(px /2m)epx /2mkT dpx
= R p2 /2mkT . (20.6)
e dpx
x

Dengan melakukan subtsitusi u2 = p2x /2mkT maka Persamaan (20.6) akan men-
jadi

R 2
kT
u2 eu du 1
x = R = kT, (20.7)

eu2 du 2

karena dengan menggunakan integral parsial dapat diperoleh bahwa


1 u2
Z Z
2
u2 eu du = ue d(u2 )
2
1h 1 u2 1 u2
Z Z
2
i
= ueu + e du = e du. (20.8)
2 u= 2 2

Dengan cara yang sama akan dapat diperoleh bahwa


20.3. RATA-RATA ENERGI POTENSIAL MIRIP PEGAS 109

1
y = kT, (20.9)
2
1
z = kT. (20.10)
2

20.3 Rata-rata energi potensial mirip pegas

Bila partikel memiliki energi potensial yang bergantun posisi seperti dalam Per-
samaan (20.2), maka rata-rata energi potensial dalam arah-x dapat dicari, mis-
alnya saja

1 2
ux = x . (20.11)
2

Kemudian dengan menggunakan prosedur yang sama seperti sebelumnya, yaitu


membuat suatu suku yang bebas terhadap x, yaitu

1
x2 , (20.12)
2

maka dapat diperoleh bahwa

R [(x2 /2)]/kT R 2

e dVp dydz (x2 /2)ex /2kT dx


ux = R [(x2 /2)]/kT R
x2 /2kT dx
e dVp dydz e
R 2
(x2 /2)ex /2kT dx
= R x2 /2kT . (20.13)

e dx

Kali ini dengan menggunakan u2 = x2 /2kT maka Persamaan (20.13) akan


menjadi

R 2
kT
u2 eu du 1
ux = R = kT, (20.14)

eu2 du 2

sehingga untuk potensial pada arah-y dan arah-z akan diperoleh pula

1
uy = kT, (20.15)
2
1
uz = kT. (20.16)
2
110 CATATAN 20. EKIPARTISI ENERGI

20.4 Rata-rata energi osilator harmonik

Suatu osilator harmonik yang memiliki energi pada arah-x seperti dalam Per-
samaan (20.2) dapat pula dihitung energi rata-ratanya pada arah-x, yaitu

R R 2 2 2 2
(px /2m + x /2) exp[(px /2m + x /2)/kT ]dxdpx
x = R R (20.17)
x
(p2 /2m + x2 /2)dxdpx

dengan kembali menggunakan prosedur yang sama dalam membahas integral


dalam d. Untuk menyelesaikan Persamaan (20.2), nyatakan koordinat x dan
px dalam bentuk polar sehingga

p2x
= r2 sin2 , (20.18)
2m
1 2
x = r2 cos2 , (20.19)
2
dxdpx = 2(m/)1/2 rdrd. (20.20)

Persamaan (20.20) diperoleh dengan memisalkan dx (/2)1/2 dx, dy =


(2m)1/2 dpx , dan dx dy = (dr)(rd). Dengan demikian Persamaan (20.17)
akan menjadi

R 2 R 2
/kT 3
d er r dr
x = R02 0
R = kT. (20.21)
d e r 2 /kT rdr
0 0

2
eau u3 du dapat dipecahkan lewat
R
Integral

Z
1
Z
2 2
eau u3 du = u2 eau d(au2 )
0 2a
  0 Z
1 2 au2 1 2
= u e + ueau du
2a u=0 a 0
1 au2 1 au2
Z Z
=0+ ue du = ue du.
a 0 a 0
Z Z
1 2 2
a= er /kT r3 dr = kT er /kT rdr. (20.22)
kT 0 0

Hasil dari Persamaan (20.21) cocok dengan hasil sebelumnya di mana rata-rata
dari suku kuadrat dari x dan px akan memberikan kontribusi energi 21 kT .
20.5. DERAJAT KEBEBASAN 111

20.5 Derajat kebebasan

Umumnya, dan lebih berguna, apabila setiap kontribusi saling bebas dalam
energi berupa suku kuadrat yang bergantung pada koordinat ruang dirujuk
sebagai suatu derajat kebebasan sebuah partikel gas. Energi rata-rata 12 kT dim-
iliki oleh setiap derajat kebebasan atau mode saling bebas dalam menyimpang
energi. Sebagai contoh, misalnya terdapat N partikel yang merupakan osila-
tor harmonik, maka derajat kebebasannya adalah 6, biasa dinyatakan dengan
f , sehingga energi sistem tak lain adalah U = 6 N 12 kT = 3N kT . Dengan
demikian energi sistem adalah

 
1
U = fN kT . (20.23)
2

Kapasitas panas pada volume konstan CV adalah

 
U 1 1
CV = = N kf = nRf. (20.24)
T V 2 2

20.6 Gas diatomik

Untuk gas diatomik selain bertranslasi, molekul gas dapat pula berotasi sehingga
ada tambahan energi dari

1 1
rot = Ix x2 + Iy y2 . (20.25)
2 2

Dengan demikian dalam hal ini derajat kebebasannya f = 5. Untuk molekul


gas diatomik atau yang lebih kompleks umumnya terdapat, secara umum, enam
derajat kebebasan.

20.7 Bukan suku kuadrat koordinat

Bila suku energi bukan berbentuk kuadrat dari koordinat dari ruang , seperti

uz = mgz, (20.26)

maka prinsip ekipertisi energi tidak akan berlaku di sini, yang berarti

1
uz 6= kT. (20.27)
2
112 CATATAN 20. EKIPARTISI ENERGI

Hal ini akan dibahas pada gas monoatomik dalam pengaruh energi potensial
gravitasi.
Catatan 21

Tambahan Informasi 1

Bab ini memberikan tambahan informasi yang mendukung kuliah dan diharap-
kan dapat memberikan motivasi bagi peserta kuliah.

21.1 Ilustrasi Cv bergantung T

Pada bab sebelumnya telah diceritakan bahwa terkait dengan prinsip ekipartisi
energi bahwa setiap suku energi yang mengandung kuadrat dari koordinat ruang
fasa akan memberikan kontribusi 12 kT terhadap energi rata-rata, yang dirujuk
sebagai suatu derajat kebebasan. Akan tetapi kapan suatu derajat kebebasan
muncul pada molekul gas, dalam hal ini harus diatomi, triatomik, atau lebih
kompleks, masih tidak jelas.

Gambar 21.1: Ilustrasi secara kualitaf perubahan CV terhadap temperatur.

113
114 CATATAN 21. TAMBAHAN INFORMASI 1

Gambar 21.1 (AbsoluteAstronomy.com, http://www.absoluteastronomy.com/


topics/Equipartition theorem [accessed 2010.04.16.06.22]) menggambarkan
bahwa CV berubah dengan temperatur. Skala temperatur sebenarnya dalam
skala logaritmik. Ilustrasi yang lebih baik karena merupakan hasil eksperimen
untuk gas hidrogen dapa dilihat dalam literatur (Francis W. Sears and Gerhard
L. Salinger, Thermodynamics, Kinetic Theory, and Statistical Thermodynam-
ics, Addison-Wesley Pub. Co., Reading, Fifth Printing, 1980, pp. 379, Fig.
12-16) dengan temperatur rotasi kira-kira pada 90 K dan vibrasi pada 7500 K.

Gambar 21.2: Sketsa mode-mode energi molekul gas triatomik, H2 O.

Gambar 21.2 memperlihatkan model-mode energi yang memberikan kontribusi


1
2 kT pada energi rata-rata tiap molekul gas (Chemistry Department Occidental
College http://departments.oxy.edu/chemistry/enthropy.htm [accessed 2010.04.
16.06.23]).

21.2 Publikasi mengenai gas ideal dan ensembel


mikrokanonik

Dengan menggunakan sedikit partikel gas misalnya terdapat publikasi oleh S.


Velasco, J. A. White, dan J. G uemez dalam European Journal of Physics 14
(4), 166-170 (1993) dengan judul Single-particle energy and velocity distribu-
tions for finite simple systems in the microcanonical ensemble yang melakukan
simulasi dengan sedikit partikel.
21.2. PUBLIKASI MENGENAI GAS IDEAL DAN ENSEMBEL MIKROKANONIK115

Gambar 21.3: Hasil yang diperoleh oleh Velasco et. al (1993).

Selain itu juga terdapat publikasi oleh F. L. Roman, J. A. White, dan S. Velasco
dalam European Journal of Physics 16 (2), 83-90 (1995) dengan judul Micro-
canonical single-particle distributions for an ideal gas in a gravitational field
yang terkait dengan topik yang akan dibahas pada bab berikutnya.

Selain itu terdapat pula publikas-publikasi dengan topik-topik yang lebih ad-
vanced, misalnya saja:

Kourosh Nozari and S. Hamid Mehdipour, Implications of Minimal


Length Scale on the Statistical Mechanics of Ideal Gas, Arxiv 0601096
(2006)

F. Becattini and L. Ferroni, The microcanonical ensemble of the ideal


relativistic quantum gas, Arxiv 0704.1967 (2007)

F. Becattini and L. Ferroni, The microcanonical ensemble of the ideal


relativistic quantum gas with angular momentum conservation, Arxiv
0707.0793 (2007)

Ying-Qiu Gu, Thermodynamics of Ideal Gas in Cosmology, Arxiv


0708.2962 (2009)
116 CATATAN 21. TAMBAHAN INFORMASI 1

Felipe Asenjo, Cristian A. Faras and Pablo S. Moya, Statistical rela-


tivistic temperature transformation for ideal gas of bradyons, luxons and
tachyons, Arxiv 0712.4368 (2009)
Catatan 22

Gas Ideal dalam Medan


Gravitasi

Dalam bab-bab sebelumnya energi gas umumnya dianggap seluruhnya dalam


bentuk energi kinetik, yang berarti walaupun wadah gas memiliki ketinggian,
energi potensial gravitasi dari molekul-molekul gas diabaikan. Dalam bab ini
energi potensial jenis ini akan diperhitungkan dan akan ditunjukkan bahwa gas
berfungsi sebagai suatu sistem multi variabel.

22.1 Sistem

Sistem yang dibahas di sini adalah terdapat suatu wadah dengan tinggi L yang
dapat diatur menggunakan piston dan luas penampang A sehingga volume ruang
yang ditempati gas pada suatu saat adalah AL. Sumbu-y diambil ke arah atas
dengan percepatan gravitasi diambil ke arah bawah. Temperatur T dalam gas
diasumsikan seragam.

Dengan demikian gas menjadi suatu sistem multivariabel yang dideskripsikan


oleh tiga variabel saling bebas T , L, dang g.

22.2 Persamaan termodinamika

Untuk sistem ini dapat dituliskan bahwa

T dS = dU + Y1 dX1 X2 dY2 (22.1)

dengan variabel ekstensif X1 adalah L dan variabel intensif Y2 adalah medan

117
118 CATATAN 22. GAS IDEAL DALAM MEDAN GRAVITASI

gravitasi dengan intensitas g. Dengan menggunakan sebagai Y1 dan sebagai


X2 maka Persamaan (22.1) akan menjadi

T dS = dU + dL dg. (22.2)

22.3 Energi total

Sistem ini memiliki energi potensiap gravitasi UG sebagai mana energi internal
gas yang hanya merupakan fungsi dari temperatur UT . Energi tiap partikel
akan menjadi

p2
= + mgy (22.3)
2m

dengan titik acuan nol untuk y diambil pada dasar ruang.

22.4 Fungsi partisi

Seperti telah dijelaskan dalam bab-bab sebelumnya, fungsi partisi dapat diny-
atakan dalam bentuk diskrit

X
Z= gj ej /kT (22.4)
j

ataupun kontinu


d
Z
Z= e/kT (22.5)
0 h3

dengan

d
gj . (22.6)
h3

Saat energi terdiri dari dua suku, yaitu suku energi kinetik dan suku energi
potensial, maka dapat dituliskan bahwa masing-masing suku akan memberikan
satu fungsi partisi. Untuk suku energi kinetik digunakan Zp yang hanya
bergantung dari koordinat momentum

d
Z
2
/2mkT
Zp = ep (22.7)
h3
22.4. FUNGSI PARTISI 119

dan untuk suku energi potensial gravitasi digunakan Zy yang hanya bergantung
dari koordinat spasial y

d
Z
Zy = emgy/kT . (22.8)
h3

Akan tetapi karena partisi sistem secara keseluruhan, tak lain adalah

Z = Zp Zy , (22.9)

maka akan lebih tepat apabila

d
Z
2
/2mkT mgy/kT
Z= ep e (22.10)
h3

sehingga Persamaan (22.7) dapat tetap digunakan akan tetapi Persamaan (22.8)
harus dimodifikasi menjadi

1
Z
Zy = emgy/kT dy. (22.11)
L

Dengan menggunakan

d = 4V p2 dp (22.12)

maka Persamaan (22.7) akan menjadi


p2
 
 
4V p
Z
p2 /2mkT
Zp = 3 (2mkT )3/2 e d
h 0 2mKT
2mkT

 
4V 1
= 3 (2mkT )3/2 ( ) =
h 4
 3/2  3/2
2mkT 2mkT
V = AL . (22.13)
h2 h2

R 2
Pemecahan ini dilakukan dengan menggunakan ex dx = dan Per-
samaan (20.8)

Sedangkan Persamaan (22.8) akan menjadi

L
1 kT  mgy  kT h mgy/kT iL
Z
Zy = emgy/kT d = e
L mg 0 kT mgL y=0
kT  
= 1 emgL/kT . (22.14)
mgL
120 CATATAN 22. GAS IDEAL DALAM MEDAN GRAVITASI

Faktor 1/L muncul karena dalam Persamaan (22.13) dihitung volume spasial
padahal untuk komponen dy telah dihitung pada Zy . Dengan demikian fungsi
partisi total adalah

"  3/2 # 
2mkT kT  
Z = AL 1 emgL/kT . (22.15)
h2 mgL

22.5 Energi bebas Helmholtz

Dengan menggunakan rumusan energi bebas Helmholtz untuk gas semi-klasik,


yaitu

F = N kT (ln Z ln N + 1), (22.16)

(Ingat bahwa rumusan F = N kT ln Z adalah untuk gas klasik, yang masih


dapat menyebabkan paradok Gibb) dapat diperoleh dua besaran baru, yaitu

   
F ln Z
= = N kT (22.17)
L T,g L T,g

dan

   
F ln Z
= = N kT , (22.18)
g T,L g T,L

yang akan memberikan

N mg
= , (22.19)
exp(mgL/kT ) 1

N kT N mL
= . (22.20)
g exp(mgL/kT ) 1

Jadi sistem ini memiliki dua persamaan keadaan yaitu (T, L, g) dan (T, L, g).
= UG /g sehingga

N mgL
UG = N kT (22.21)
exp(mgL/kT ) 1

sehingga dengan menggunakan


22.6. ENTROPI 121

 
2 ln Z 5 N mgL
U = N kT = N kT (22.22)
T L,g 2 exp(mgL/kT ) 1

dan U = UG + UT dapat diperoleh bahwa

3
UT = N kT (22.23)
2

yang merupakan energi dalam gas saat tidak terdapat medan gravitasi, yang
hanya bergantung dari temperatur.

22.6 Entropi

Dengan menggunakan hubungan antara entropi S, energi dalam sistem U , dan


energi bebas Helmholtz F untuk gas semi-klasik dari Persamaan (22.16)

F = U T S, (22.24)

maka dapat diperoleh entropi sistem adalah

U F
S=
  T T
1 5 N mgL 1
= N kT [N kT (ln Z ln N + 1)]
T 2 exp(mgL/kT ) 1 T
 
5 N mgL/T
= Nk + N k(ln Z ln N + 1). (22.25)
2 exp(mgL/kT ) 1

22.7 Distribusi partikel sebagai fungsi keting-


gian

Jumlah total partikel dalam sistem secara umum adalah

X X
N= nj = gj ej /kT (22.26)
j j

atau dalam bentuk integralnya adalah


d
Z
N= e/kT . (22.27)
0 h3
122 CATATAN 22. GAS IDEAL DALAM MEDAN GRAVITASI

Dalam hal ini, apabila ingin dibahas bagaimana distribusinya terhadap keting-
gian y maka, suku energi akan diperhatikan hanya bagian yang bergantung y
dari Persamaan (22.3). Dengan demikian Persamaan (22.27) akan menjadi

d
Z
2
N= ep /2mkT mgy/kT
h3
Z Z L
2 d dy
= e ep /2mkT 3 emgy/kT
0 h 0 L
Z L
dy
= e Zp emgy/kT . (22.28)
0 L

Jumlah partikel yang berada antara ketinggian y dan y + y adalah

dy
dNy = e Zp emgy/kT (22.29)
L

Persammaan (22.28) dapat diubah menjadi

L
dy
Z

N = e Zp emgy/kT = e Zp Zy , (22.30)
0 L

selanjutnya, dengan menggunakan Persamaan (22.30), Persamaan (22.29) dapat


dituliskan kembali menjadi

N mgy/kT dy
dNy = e . (22.31)
Zy L

Volume sebuah daerah tipis adalah Ady sehingga jumlah partikel per satuan
volume pada ketinggian y adalah

dNy 1 dNy
ny = = . (22.32)
V A dy

Dengan menggunakan persamaan keadaaan untuk gas idel (pV = N kT ) maka


tekanan pada ketinggial y adalah

py = ny KT. (22.33)

Bentuk eksplisit dari ny dapat dihitung menggunakan Persamaan (22.31), yaitu

N mgy/kT 1
ny = e . (22.34)
AZy L
22.7. DISTRIBUSI PARTIKEL SEBAGAI FUNGSI KETINGGIAN 123

Nilai dari Zy telah dihitung dalam Persamaan (22.14) sehingga dapat dituliskan
bahwa

N 1 mgy/kT 1
ny = e
A Zy L
 
N mgL 1 1
= emgy/kT
A kT 1 emgL/kT L
N mg 1 emgy/kT
= . (22.35)
A kT 1 emgL/kT

Substitusikan Persamaan (22.35) ke dalam Persamaan (22.33) sehingga diper-


oleh

N mg emgy/kT
py = (22.36)
A 1 emgL/kT

Pada y = 0 tekanan dinyatakan sebagai p0 sehingga

N mg 1
p0 = , (22.37)
A 1e mgL/kT

sehingga Persamaan (22.36) dapat dituliskan dalam bentuk yang lebih kompak,
yaitu menjadi

py = p0 emgy/kT . (22.38)

Persamaan (22.38) ini menceritakan bahwa tekanan berkuran secara eksponen-


sial dengan bertambahnya ketinggian. Persamaan ini dikenal sebagai persamaan
barometrik atau hukum atmosfer (the law of atmospheres), yang dapat pula di-
turunkan dari prinsip hidrostatik dan persamaan keadaan gas ideal, yaitu

p = p0 gy dp = gdy
M Nm
= =
V V
1 p
pV = N kT =
V N kT
p m
= Nm = p
N kT kT
m dp
dp = gdy = pgdy = (mg/kT )dy
kT p
ln p ln p0 = mgy/kT
p = p0 emgy/kT . (22.39)
124 CATATAN 22. GAS IDEAL DALAM MEDAN GRAVITASI

Pada bagian atas wadah gas di mana y = L akan diperoleh bahwa

N mg emgL/kT N mg 1
pL = = = , (22.40)
A 1 emgL/kT A emgL/kT 1 A

sehingga

= pL A, (22.41)

dengan kuantitas adalah gaya yang melawan piston pada bagian atas wadah
gas. Bila terjadi perubahan tinggi piston sebesar dL ke arah atas maka kerja
adalah

dW = dL = pL AdL = pL dV, (22.42)

adalah kerja yang dihasilkan saat gas mengembang.

22.8 Percobaan Perrin

Jean Perrin (1870-1942), seorang fisikawan Perancis, menggunakan Persamaan


(22.38) untuk menentukan bilangan Avogadro NA . Caranya ini termasuk salah
satu dari cara-cara terawal yang presisi.

Perrin tidak menggunakan molekul-molekul gas melainkan partikel-partikel


mikroskopik yang tersuspensi dalam fluida dengan densitas sedikit lebih kecil
dari densitas partikel-partikel mikroskopik tersebut, yang akan mereduksi nilai
efektif dari percepatan gravitasi g. Jumlah partikel pada ketinggian-ketinggian
yang berbeda dihitung menggunakan sebuah mikroskop.

Jika N1 dan N2 adalah jumlah rata-rata partikel pada ketinggian y1 dan y2 ,


maka dengan menggunakan Persamaan (22.31) akan dapat diperoleh bahwa

N1
= emg(y1 y2 )/kT . (22.43)
N2

Semua kuantitas dalam Persamaan (22.43) dapat diukur secara eksperimen ke-
cuali konstanta Boltzmann k. Jadi persamaan tersebut dapat diselesaikan untuk
mendapatkan k. Kemudian NA dapat dicari melalui R/k, di mana konstanta
universal gas R diperoleh dari eksperimen lain pada saat itu.

Perrin menyimpulkan bahwa nilai NA terletak antara 6.5 1026 dan 7.2 1026 ,
yang terbandingkan dengan hasil eksperimen terbaik saat ini 6.022 1026
molekul per kmol.
Catatan 23

Gas diatomik

Dalam bab ini akan dibahas mengenai gas diatomik, di mana partikel gas bukan
lagi merupakan satu atom (gas monoatomik) melainkan merupakan molekul
yang terdiri dari dua atom (gas diatomik).

23.1 Suku-suku energi

Dalam sebuah molekul diatomik energi dapat dianggap terbangun dari lima
buah kontribusi yang saling bebas satu sama lain. Kelima kontribusi tersebut
muncul akibat

translasi molekul sebagai suatu kesatuan,

rotasi molekul sebagai suatu kesatuan,

gerak vibrasi dua atom sepanjang sumbu molekul,

gerak elektron-elekron di sekeliling inti, dan

spin nuklir.

Tingkat-tingkat energi yang tersedia bagi kelima jenis gerak perlu terkuanti-
sasi dan hanya pada kasus gerak translasi dapat diperlakukan tingkat-tingkat
energinya sebagai suatu tingkat kontinu klasik. Pada kasus yang lain tingkat-
tingkat energi harus diperlakukan secara benar-benar diskrit kecuali, mungkin,
pada temperatur amat tinggi.

Dengan melihat ini, dapat disimpulkan untuk sementara, bahwa untuk sistem
yang lebih rumit, dapat dilakukan prosedur yang sama asalkan deskripsi sistem
telah lengkap untuk seluruh faktor yang dapat menyumbangkan energi.

125
126 CATATAN 23. GAS DIATOMIK

23.2 Fungsi-fungsi partisi

Untuk mempelajari properti termodinamika suatu gas ideal diatomik perlu diru-
muskan fungsi partisi sebuah molekul gas dalam suku-suku fungsi-fungsi partisi
yang terpisah untuk setiap bentuk gerak yang berkontribusi. Dengan mengikuti
cara dalam bab gas ideal dalam pengaruh medan gravitasi, adalah mungkin un-
tuk melaukan faktorisasi fungsi partisi total Z seperti

Z = Zt Zr Zv Ze Zn , (23.1)

di mana Zt adalah fungsi partisi untuk gerak translasi, Zr adalah fungsi partisi
untuk gerak rotasi, Zv adalah fungsi partisi untuk gerak vibrasi, Ze adalah
fungsi partisi untuk gerak elektron-elektron, dan Zn adalah fungsi partisi untuk
gerak spin inti.

23.3 Fungsi partisi gerak translasi

Gerak translasi molekul diatomik cukup mirip sehingga dapat dianalogikan


dengan gerak translasi molekul tak-berstruktur (structureless) yang telah ditu-
runkan untuk gas ideal monoatomik, sehingga Zt tak lain adalah seperti dalam
Persamaan (22.13), yaitu

 1/2
2mkT
Zt = V . (23.2)
h2

23.4 Fungsi partisi gerak rotasi

Untuk mementukan fungsi partisi gerak rotasi molekul diatomik perlu dituliskan
tingkat energi rotasi yang diperbolehkan j dalam bentuk mekanika kuantum,
yaitu

h2
j = j(j + 1) , (23.3)
8 2 I

di mana I adalah momen inersia molekul terhada suatu sumbu yang melewati
titik pusat massa molekul dan tegak lurus terhadap garis yang menghubungkan
kedua atom dalam molekul diatomik. Sedangkan j adalah bilangan kuantum
momentum angular total gerakan rotasi. Salah satu penurunan Persamaan
(23.3) dapat dilihat dalam S. Glasstone, Theoretical Chemistry, D. Van Nos-
trand Co. Inc., (1947).

Untuk setiap nilai bilangan kuantum j bilangan kuantum magnetik mj hanya


dapat memiliki nilai bilangan bulat antara j dan j, yang berarti terdapat
23.5. FUNGSI PARTISI GERAK VIBRASI 127

(2j + 1) nilai. Dengan demikian setiap tingkat energi merepresentasikan (2j + 1)


keadaan energi dan, dengan menggunakan informasi ini sebagai degenerasi, akan
memberikan fungsi partisi untuk gerak rotasi menjadi


X
X
Zr = (2j + 1)ej /kT = (2j + 1)ej(j+1)K/kT , (23.4)
j j

dengan

h2
K= . (23.5)
8 2 I

23.5 Fungsi partisi gerak vibrasi

Gerak vibrasi molekul dapat diasumsikan, pada suatu aproksimasi yang baik,
sebagai suatu bentuk osilator harmonik dan bebas dari segala distorsi tak-
harmonik. Oleh karena itu memungkinkan untuk menggunakan hasil fungsi
partisi untuk osilator harmonik satu-dimensi, yaitu

1
e 2 h/kT
Zv = , (23.6)
1 eh/kT

dengan adalah frekuensi karakteristik vibrasi molekular, yang ditentukan den-


gan menggunakan massa atom yang membentuk molekul dan kopling alami
antara keduanya.

23.6 Fungsi partisi gerak elektron

Fungsi partisi elektronik biasanya dapat direpresentasikan dengan akurasi yang


cukup dalam dua suku pertama energi penjumlahan normal. Bila energi diper-
lukan untuk mengeksitasi sebuah elektron dari tingkat energi dasar ke tingkat
energi tereksitasi pertama, kedua, ketiga, dan seterusnya, dengan energi e1 , e2 ,
e3 , dan seterusnya, maka fungsi partisi dapat secara efektif dituliskan sebagai

Ze = g0 + g1 ee1 /kT + g2 ee2 /kT + g3 ee3 /kT + .., (23.7)

di mana g0 adalah degenerasi untuk tingkat energi dasar dan g1 , g2 , g3 , dan


seterusnya adalah denerasi tingkat energi tereksitasi pertama, kedua, ketiga,
dan seterusnya.
128 CATATAN 23. GAS DIATOMIK

Energi e1 dan e2 secara umum amat sangat besar apabila dibandingkan dengan
energi termal kT kecuali dalam kasus pada temperatur amat tinggi. Oleh karena
itu umumnya mungkin untuk mendekati Persamaan (23.7) dengan

Ze g0 + g1 ee1 /kT . (23.8)

23.7 Fungsi partisi spin nuklir

Umumnya fungsi partisi spin nuklir hanya melibatkan perkalian dengan suatu
faktor konstan, misalnya

Zn . (23.9)

23.8 Fungsi partisi lengkap

Dengan demikian, fungsi partisi lengkap molekul diatomik dapat dituliskan den-
gan melakukan subtsitusi dari Persamaan (23.2), (23.4), (23.6), (23.8), dan
(23.9) ke dalam Persamaan (23.1), yaitu


 1/2 X
2mkT (2j + 1)ej(j+1)K/kT
Z=V
h2 j
1
!
e 2 h/kT 
e1 /kT

g 0 + g 1 e Zn (23.10)
1 eh/kT

Fungsi partisi total suatu gas yang tersusun atas N molekul diatomik identik
adalah

ZN
Ztotal = . (23.11)
N!

Dengan menggunakan hubungan antara fungsi partisi dan energi dalam sistem,
yaitu

 
ln Z
U = N kT 2 , (23.12)
T V

maka dapat diperoleh bahwa


23.9. PANAS SPESIFIK GAS 129

     
2 ln Zt ln Zr ln Zv
U = N kT + +
T
V T T
  V   V
ln Ze ln Zn
+ +
T V T V

3
ln X
=N kT + kT 2 (2j + 1)ej(j+1)K/kT
2 T j

g1 e1 ee1 /kT
  
1 1
+h + h/kT + + 0 . (23.13)
2 e 1 g0 + g1 e1 ee1 /kT

Persamaan (23.13) menggambarkan bahwa tidak ada kontribusi dari suku Zn


yang dianggap tidak bergantung dari temperatur.

Besarnya berbagai kontribusi terhadap U dalam Persamaan (23.13) dapat lebih


mudah diparesiasikan dan umum untuk mendefisikan temperatur karakteristik
untuk rotasi rot , vibrasi vib , dan energi-energi elektronik e .

K h2
rot = = , (23.14)
k 8 2 Ik
h
vib = , (23.15)
k
e1
e = , (23.16)
k

yang nilai-nilainya dapat diperoleh dari parameter-parameter fisis molekul dan


secara umum ditemukan bahwa terdapat pertidaksamaan

rot vib e . (23.17)

Beberapa nilai rot , vib , dan e dapat dilihat dalam Tabel 23.1 (gabungan dari
Sears and Salinger, 1975 dan Pointon, 1967; data yang kosong tidak terdapat
pada kedua sumber tersebut).

23.9 Panas spesifik gas

Pada temperatur rendah di mana T rot jelas dari bentuk fungsi partisi
lengkap dalam Persamaan (23.10) bahwa hanya sejumlah kecil molekul, yang
merupakan fraksi tak-signifikan, yang berotasi, bervibrasi, atau energi elek-
troniknya berada di atas tingkat energi dasarnya. Dalam kasus energi total,
seperti dalam rumusan Persamaan (23.13), hanya merupakan kontribusi dari
gerak translasi ditambah dengan energi titik nol gerak rotasi ( 12 N h), dan en-
ergi tingkat dasar dari elektron-elektron.
130 CATATAN 23. GAS DIATOMIK

Tabel 23.1: Temperatur karakteristik bagi beberapa molekul diatomik.


Bahan rot (K) vib (K) e (K)
H2 85.5 6140 118000
Cl2 0.35 810 25500
O2 2.09 2260 11000
OH 27.5 5360
HCl 15.3 4300
CH 20.7 4100
CO 2.77 3120
NO 2.47 2740
Br2 0.117 470
Na2 0.224 230
K2 0.081 140

Dalam hal ini dapat dituliskan bahwa

 3/2
2mKT
Zt = V , (23.18)
h2
Zr 1, (23.19)
12 h/kT
Zv e , (23.20)
Ze g 0 , (23.21)
Zn , (23.22)
 3/2
2mKT 1
Z=V 1 e 2 h/kT g0 Zn , (23.23)
h2

sehingga

3 1 h 1 h/kT
U (T rot vib ) = N kT + 0 e 2 +0+0
2 2 k
3 1 1 3
= N kT vib e 2 vib /T N kT (23.24)
2 2   2
U 3
CV = = N k. (23.25)
T V 2

Bila termperatur dinaikkan sampai orde rot maka sebagian molekul akan terek-
sitasi sehingga menempati tingkat energi rotasi di atas tingkat energi dasar. En-
ergi rotasi akan mulai berkontribusi kepada energi total sistem dan panas spesifik
gas, sementara gerak vobrasi masih berada pada tingkat energi dasarnya.

Pada keadaan T rot ini dapat pula dilakukan aproksimasi untuk fungsi
partisi gerak rotasi Zr dengan mengubah somasi menjadi integrasi
23.9. PANAS SPESIFIK GAS 131


X
Zr = (2j + 1)ej(j+1)rot /T ,
j
Z
Zr = (2j + 1)ej(j+1)rot /T dj
0
Z    
(j+ 12 )2 rot /T 1 1
= erot /4T e 2 j+ d j+
0 2 2
T
= erot /4T . (23.26)
rot

Kemudian dengan

rot T
1 erot /4T 1 Zr . (23.27)
T rot

akan membuat fungsi partisi lengkap menjadi

 3/2
2mKT T 1
Z =V e 2 h/kT g0 Zn , (23.28)
h2 rot

sehingga

3 1 h 1 h/kT
U (rot T vib ) = N kT + N kT e 2 +0+0
2 2 k
5 1 1 5
= N kT vib e 2 vib /T N kT (23.29)
2 2   2
U 5
CV = = N k. (23.30)
T V 2

Bertambahnya nilai CV sebesar 2( 21 kT ) ini sesuai dengan bertambahnya dua


derajat kebebasan pada rotasi molekul gas diatomik. Pada temperatur ini energi
vibrasi belum berkontribusi terhadap energi total sistem sehingga juga gerak
vibrasi tidak berkontribusi pada panas spesifik gas sampai temperatur mencapai
temperatur karakteristik vibrasi vib .

Bila temperatur dinaikkan sehingga T vib maka molekul gas akan mulai
bervibrasi dengan fungsi partisi vibrasi Zv seperti dalam Persamaan (23.6) mulai
berperan yang akan memberikan sumbangan ke energi total sebesar 2( 21 kT )
karena terdapat dua mode vibrasi yaitu simetri dan asimetri. Dalam hal ini
tetap berada dalam kondisi T e . Dengan demikian dapat dituliskan bahwa

7
U= N kT, (23.31)
2
132 CATATAN 23. GAS DIATOMIK

7
CV = N k. (23.32)
2

Peningkatan CV ini sebagai fungsi temperatur akan seperti pernah diilustrasikan


dalam Gambar 21.1.

Penurunan panas spesifik sampai sejauh ini telah mengabaikan setiap kontribusi
yang mungkin dari tingkat-tingkat energi elektronik. Pengabaian ini mengimp-
likasikan bahwa untuk semua molekul diatomik yang memenuhi kondisi T e
akan terpenuhi kecuali pada temperatur amat tinggi.

Untuk oksigen misalnya, e kira-kira 11000 K, yang termasuk terendah untuk


molekul diatomik. Teramati bahwa gas yang tersusun atas molekul oksigen
tidak akan terdapat kontribusi yang signifikan dari tingkat energi elektronik
sampai temperatur gas mencapai kira-kira 2000 K. Dalam kasus ini kontribusi
tingkat-tingkat energ elektronik cukup dengan dalam bentuk faktor mutiplikasi
g0 , degenerasi tingkat energi dasar.

Pada kasus lain, molekul NO memiliki nilai e yang jauh lebih kecil dari molekul
gas oksigen dan tingkat-tingkat energi elektronik memberikan kontribusi ke
panas spesifik saat temperatur lingkungan 80 K. Kasus ini, akan tetapi, meru-
pakan kasus yang tidak umum.

Kontribusi dari spin nuklir akan muncul dari degenerasi keadaan energi spin
nuklir. Bila jumlah total spin yang disumbangkan oleh dua inti suatu molekul
memiliki bilangan kuantum I maka keadaan spin akan berjumlah (2I + 1),
atau degenerasinya adalah (2I + 1). Nilai (2I + 1) jumlah keadaan ini memiliki
energi yang hampir sama berkaitan dengan nilai-nilai yang diperbolehkan untuk
dimiliki bilangan kuantum magnetik mI , I, (I 1), .., 0, (I 1), I. Dengan
demikian fungsi partisi Zn dapat digantikan dengan suku tingkat dasar 2I + 1),
sehingga Z = Zt Zr Zv g0 (2I + 1) yang mengatakan bahwa suku g0 (2I + 1) ini
akan berkontribusi pada entropi sistem S, tetapi tidak pada energi sistem U
atau panas spesifik CV .
Catatan 24

Gas Bose-Einstein

Dalam bab ini akan dibahas mengenai gas Bose-Einstein yang memenuhi statis-
tik Bose-Einstein. Photon dan phonon termasuk di dalamnya. Untuk saat ini
hanya photon yang akan dibahas

Bila molekul-molekul dalam suatu gas biasa memiliki momentum angular inte-
gral dalam satuan h/2 maka, dapat dikatakan dengan sangat, bahwa molekul-
molekul tersebut adalah boson dan akan mematuhi statistik Bose-Einstein.

24.1 Distribusi molekul gas

Distribusi molekul-molekul gas meliputi berbagai tingkat energi diberikan oleh


oleh

gj
Nj = , (24.1)
e ej /kT 1

di mana

A = e . (24.2)

Karena setiap keadaan energi yang diperbolehkan membutuhkan suatu volume


h3 dalam ruang fasa, bobot suatu tingkat energi, atau keadaan-keadaan yang
dapat dipertimbangkan memenuhi suatu volume d dalam ruang fasa akan men-
jadi

d
gj . (24.3)
h3

Dengan menggunakan bahwa elemen ruang fasa enam-dimensi bahwa

133
134 CATATAN 24. GAS BOSE-EINSTEIN

d = (dV )(dVp ) = (dxdydz)(dpx dpy dpz), (24.4)

maka dapat dituliskan, apabila momentum dilihat dalam koordinat polar, men-
jadi

d = 4V p2 dp. (24.5)

Dengan hanya memperhitungkan energi kinetik

p2
= p = 2m, (24.6)
2m

dapat diperoleh bahwa

2md = 2pdp. (24.7)

Substitusi Persamaan (24.7) dan (24.6) ke dalam Persamaan (24.5) akan mem-
berikan

3 1
d = 2V (2m) 2 2 d. (24.8)

Dengan demikian Persamaan (24.3), dengan substitusi dari Persamaan (24.8),


akan menjadi

3 1
d 2V (2m) 2 2 d
g()d = = , (24.9)
h3 h3

yang menyatakan jumlah keadaan energi yang tersedia dalam rentang energi
antara dan + d untuk suatu volume V . Di sini g() adalah kerapatan
keadaan energi (density of states).

Jumlah molekul-molekul yang memiliki energi dalam rentang dan + d


diberikan oleh Persamaan (24.1) dan (24.10), yaitu

3 1
2V (2m) 2 2 d 1
N ()d = 3 /kT
. (24.10)
h e e
1

Nilai parameter A atau untuk distribusi ini dapat ditentukan untuk kondisi
bahwa

Z
N ()d = N, (24.11)
0
24.2. GAS FOTON DAN RADIASI BENDA HITAM 135

dengan N adalah jumlah total molekul dalam volume V .

Secara umum bentuk integral dalam Persamaan (24.11) sulit untuk dipecahkan
secara eksak, akan tetapi dapat dilihat, bahwa dalam beberapa kasus prak-
tis, nilai A untuk gas cukup kecil sehingga menyebabkan suku bernilai 1 pada
penyebut dalam Persamaan (24.10) dapat diabaikan. Bila kondisi ini dipenuhi
distribusi akan mendekati distribusi Maxwell-Boltzmann, dan karena molekul-
molekul gas akan tersebar di antara keadaan-keadaan energi, gas dikatakan tidak
terdegenerasi. Dengan demikian integrasi Persamaan (24.10) akan menghasilkan
seperti integrasi dalam distribusi Maxwell-Boltzmann yang memberikan

N h3
A= 3 (24.12)
V (2mkT ) 2

dan

N h3
 
= ln A = ln 3 . (24.13)
V (2mkT ) 2

Dikarenakan nilai exponen e/kT selalu lebih besar (atau setidaknya sama den-
gan) satu untuk semua nilai energi kondisi yang akan didekati oleh Persamaan
(24.10) untuk menjadi distribusi klasik adalah membuat A 1.

Bila digunakan nilai-nilai N , V , dan m untuk helium maka akan diperoleh nilai
A untuk tekanan atmosfer, yaitu

A 3 106 untuk T = 300 K, dan

A 1.5 101 untuk T = 4 K.

Jadi bahkan untuk temperaturn 4 K pun, kondisi (1/A)e/kT 1 tetap ter-


penuhi dan gas helium akan berlaku, untuk suatu aproksimasi yang baik, seba-
gai suatu gas klasik.

24.2 Gas foton dan radiasi benda hitam

Radiasi elektromagnetik yang berada dalam suatu ruang tertutup bertemper-


atur tetap dapat dipertimbangkan sebagai suatu sistem foton-foton dengan
berbagai nilai energi. Dan karena foton-foton memiliki momentum angular in-
tegral dalam satuan h/2 maka mereka akan secara alami berkelakuan sebagai
boson dan dapat diasumsikan bahwa suatu gas foton akan memiliki distribusi
energi yang diberikan oleh statistik Bose-Einstein. Akan tetapi, terdapat dua
hal yang harus diperhatikan.

Pertama, foton dapat diserap dan dipancarkan kembali oleh dinding lingkun-
gan tertutup yang bertemperatur tetap, dengan demikian jumlah foton dalam
136 CATATAN 24. GAS BOSE-EINSTEIN

P
lingkungan tersebut tidaklah tetap. Dengan demikian kondisi
P j Nj = N atau
dNj = 0 dalam

X X
d ln W + dNj + j dNj = 0 (24.14)
j j

tidak dapat terpenuhi. Agar Persamaan (24.14) masih dapat berlaku maka
perlu dipilih bahwa = 0 sehingga A = 1.

Kedua, energi foton berbentuk h, di mana adalah frekuensi radiasi. Oleh


karena itu lebih memudahkan apabila distribusi energi diungkapkan dalam
freku-ensi atau panjang gelombang foton.

Dengan menggunakan rumusan panjang gelombang de Broglie

h
= (24.15)
p

maka

h
dp = d (24.16)
2

Persamaan (24.5) akan menjadi

 2 
h3

2 h h
d = 4V p dp = 4V 2 d = 4V 4 d (24.17)

Selanjutnya karena setiap foton memiliki kemungkinan polarisasi pada dua arah
maka jumlah keadaan energi yang diperbolehkan atau mode, dalam rentang
antara dan + d, untuk setiap satuan volume adalah (dengan mengambil
nilai positifnya)

h3
    
d 1 1 d
g()d = 3 = (2) 4V 4 d = 8 4 . (24.18)
h V h3

Kemudian dengan menggunakan = hc dalam Persamaan (24.1) akan diper-


oleh bahwa

8 d
N ()d = 4 hc/kT
. (24.19)
e 1

Rumusan untuk E() dapat diperoleh dengan mengalikan Persamaan (24.19)


dengan hc/ sehingga diperoleh
24.3. HUKUM RADIASI WIEN 137

8hc d
E()d = . (24.20)
5 (ehc/kT 1)

Ekspresi dalam Persamaan (24.20) dikenal sebagai hukum radias Planck untuk
distribusi energi spektral dari energi radiasi dalam suatu lingkungan tertutup
bertemperatur konstan. Ilustrasi distribusi energi spektral dapat dilihat dalam
Gambar 24.1.

6e+31 T3
T3 > T2 > T1

4e+31
E()

T2

2e+31
T1

0
2e-07 4e-07 6e-07 8e-07 1e-06

Gambar 24.1: Distribusi energi terhadap panjang gelombang E() untuk fungsi
c1 5 (ec2 /T 1)1 dengan c1 = 25 dan c2 = 4 104 , untuk berbagai tem-
peratur: T1 = 250K, T2 = 300K, dan T3 = 330K.

Beberapa pengamatan dapat dibuat terkait dengan Persamaan (24.20).

24.3 Hukum radiasi Wien

Ekspresi untuk E() dalam bentuk

 
1
E() = f (T ) (24.21)
5

diprediksikan oleh hukum radiasi Wien yang hanya berdasarkan argumen-


argumen termodinamika. Aproksimasi Persamaan (24.20) pada daerah panjang
gelombang dengan nilai kecil, di mana ehc/kT 1, yang akan menghasilkan
formula distribusi Wien

8hc hc/kT
E()d e d, (24.22)
5

yang awalnya diusulkan sebagai suatu pencocokan empiris terhadap dapat


eksperimen pada pengukuran di daerah panjang gelombang kecil.
138 CATATAN 24. GAS BOSE-EINSTEIN

24.4 Formula Rayleigh-Jeans

Pada daerah dengan panjang gelombang bernilai besar di mana dapat dilakukan
aproksimasi ehc/kT 1 + hc/kT sehingga Persamaan (24.20) akan menjadi

8kT
E()d d, (24.23)
4

yang diturunkan pada asumsi bahwa setiap (8/4 )d foton memiliki energi
osilator harmonik klasik kT . Rumusan ini tidak baik untuk menjelaskan radiasi
pada panjang gelombang pendek karena akan menghasilkan E()

6e+31

Hukum Wien
4e+31 Hukum Rayleigh-Jeans
E()

2e+31
Hukum Plank

0
2e-07 4e-07 6e-07 8e-07 1e-06

Gambar 24.2: Ketiga hukum radiasi: Wien (cocok sekali pada daerah panjang
gelombang kecil), Rayleigh-Jeans (berkelakukan baik pada panjang gelombang
besar), dan Plank (yang cocok untuk semua daerah panjang gelombang).

Dengan menggunakan parameter seperti dalam Gambar 24.1 ketiga hukum ra-
diasi diilustrasikan dalam Gambar 24.2.

24.5 Hukum Stefan-Boltzmann

Bila lingkungan tertutup dengan temperatur tetap, yang digunakan dalam


memi-salkan gelombang elektromagnetik sebagai gas boson, dibuat berlubang
maka gelombang elektromagnetik akan teradiasi ke luar melewati lubang terse-
but.

Dari teori kinetik diketahui bahwa gas dengan jumlah molekul setiap satuan vol-
ume n jumlah molekul yang keluar pada setiap satuan luas dalam satuan waktu
adalah 41 nv, di mana v adalah laju rata-rata molekul. Kemudian, bila tidak ada
seleksi terhadap panjang gelombang tertentu untuk peristiwa absorbsi ataupun
emisi dari radiasi oleh lubang (yang berperan sebagai suatu peradiasi benda
24.5. HUKUM STEFAN-BOLTZMANN 139

hitam) maka jumlah foto yang diemisikan dalam rentang panjang gelombang
dan + d tiap satuan luas lubang tiap satuan waktu diungkapkan dalam
Nrad ()d.

Untuk foton yang laju rata-ratanya adalah c maka dapat dituliskan bahwa

c
Nrad ()d = N ()d. (24.24)
4

Dengan menggunakan Persamaan (24.19) maka diperoleh bahwa

2c d
Nrad ()d = (24.25)
4 ehc/kT 1

dan

2hc2 d
Erad ()d = 5 hc/kT
. (24.26)
e 1

Total energi tiap satuan volume suatu lingkungan tertutup dengan temperatur
tetap diperoleh dengan melakukan integrasi Persamaan (24.20) terhadap seluruh
rentang panjang gelombang


8hc d
Z Z
E= E()d =5 (ehc/kT 1)
0 0
 4 Z 3  5 4
8h kT t dt 8 k
= 3 t
= T 4, (24.27)
c h 0 e 1 15h3 c3
Z 3
t dt X 1 4
= 6 = , (24.28)
0 et 1 n=1
n4 15
c
Erad = E = T 4 , (24.29)
4
2 5 k 4
= . (24.30)
15h3 c2

Hukum Stefan-Boltzmann dalam Persamaan (24.29) menggambarkan


bagaimana kalor dirambatkan secara radiatif (tanpa perantaraan medium,
sebagaimana medium diperlukan pada peristiwa kondusi atau konveksi).
Konstant Stefan diungkapkan dalam Persamaan (24.30).
140 CATATAN 24. GAS BOSE-EINSTEIN
Catatan 25

Gas Fermi-Dirac

Partikel-partikel yang termasuk dalam keluarga fermion, yang memiliki dua


jenis spin (spin up dan down), dapat pula dimisalkan sebagai gas fermin.

25.1 Distribusi partikel

Distribusi partikel dalam statistik Fermi-Dirac adalah

gj
Nj = . (25.1)
e(+j ) + 1

Setiap fermion memiliki kemungkinan dua keadaan spin (up dengan spin 12 dan
down dengan spin 21 ) sehingga dalam representasi ruang enam-dimensi d
Persamaan (25.1) dengan menggunakan Persamaan (24.8) akan menjadi

 3 1
 1  1

N ()d = (2) 2V (2m) d
2 2 . (25.2)
h3 e(+) + 1

Persamaan (25.2) dapat dituliskan menjadi

N ()d = g()f ()d, (25.3)

dengan

  23
2m 1
g() = 4V 2 , (25.4)
h2

141
142 CATATAN 25. GAS FERMI-DIRAC

1
f () = , (25.5)
e(F )/kT ) + 1
1
= , (25.6)
kT
F = kT . (25.7)

Persamaan (25.5) merupakan fungsi Fermi dengan energi Fermi tak lain didefin-
isikan untuk menggantikan parameter seperti dalam Persamaan (25.7).

25.2 Fungsi Fermi

Sebelum membahas gas fermin perlu terlebih dahulu melihat fungsi Fermi dan
bagaimana kelakukannya.
Catatan 26

Ensemble Kanonis

Semua pembahasan yang telah dilakukan adalah terbatas pada sistem-sistem


yang mempunyai energi total tetap dan banyaknya partikel dalam sistem juga
tetap, dan tidak ada interaksi antar partikel. Meskipun demikian, studi tentang
sistem-sistem tersebut dapat memberikan hasil-hasil yang berguna pada daerah
permasalahan fisika yang cukup luas. Pembatasan yang diterapkan pada sistem
juga akan membatasi daerah berlakunya hasil yang diperoleh.

26.1 Ensemble

Agar metoda statistik dapat diperlebar ke daerah permasalahan fisika yang


lebih luas, akan diperkenalkan konsep ensemble. Untuk maksud ini, tinjau
sekumpulan sistem yang semuanya mempunyai volum sama dan partikelnya
sejenis. Lalu, bergantung pada kondisi yang diterapkan pada sistem-sistem
tersebut, ada banyak keadaan mikro dalam sistem-sistem itu. Jika untuk se-
tiap keadaan mikro, ada paling tidak satu sistem lain dalam kumpulan sis-
tem tersebut yang mempunyai keadaan mikro sama, maka kumpulan sistem itu
dikatakan membentuk sebuah ensemble. Sekarang, karena setiap kemungkinan
keadaan dari sebuah sistem akan direpresentasikan dalam ensamble, mempela-
jari sebuah ensemble akan ekivalen dengan mempelajari sebuah sistem yang
susunan partikel-partikelnya bervariasi terhadap waktu untuk semua kemungk-
inan keadaan mikronya. Maka, perata-rataan sifat sistem untuk seluruh sistem
yang membentuk ensemble akan memberikan hasil yang sama sebagaimana jika
perata-rataan ini dilakukan untuk sebuah sistem yang keadaannya bervariasi
dengan waktu.

Dalam bab ini, akan dibahas sistem-sistem yang membentuk ensemble kano-
nis (canonical ensemble). Dalam ensemble ini, banyaknya partikel pada tiap
sistem adalah sama dan merupakan bilangan konstan, dan temperatur tiap sis-
tem (bukan energinya) adalah sama dan merupakan bilangan konstan,. Kondisi
seperti ini memperbolehkan kemungkinan adanya pertukaran energi antar sis-

143
144 CATATAN 26. ENSEMBLE KANONIS

tem dalam sebuah ensemble dan juga interaksi antar partikel dalam sebuah
sistem. Untuk sistem-sistem yang boleh mengalami perubahan energi dan jum-
lah partikelnya tidak akan dibahas dalam bab ini. Sebagai catatan , ensemble
yang sistem-sistemnya dibatasi dengan N (banyaknya partikel), V (volum) dan
U (energi) yang tetap dinamakan ensemble mikrokanonis.

26.2 Ensemble yang bertemperatur konstan

Karena sistem-sistem dalamsebuah ensemble kanonis didefinisikan bertemper-


atur sama, maka sistem-sistem ini kontak termal dengan sistem tetangganya
dalam sebuah ensemble. Bila sistem-sistem ini dalam keadaan keseimbangan
termodinamik, keseluruhan ensemble dapat dipandang membentuk sebuah en-
closure bertemperatur konstan.

Ilustrasikan ensemble kanonis atau ensemble bertemperatur konstan ditun-


jukkan pada Gambar 26.1. Sistem mempunyai nilai N (banyaknya partikel),
V (volum) dan T (temperatur) yang tetap. Batas antar sistem adalah dinding
diatermik yang permeabel terhadap transfer energi (kalor). Secara teori, energi
sebuah sistem dalam ensemble kanonis berubah dengan waktu dari nilai energi
sama dengan nol sampai nilai total energi ensemble tersebut.

Gambar 26.1: Ilustrasi suatu ensembel kanonis.

Sekarang, tinjaulah sebuah sistem dalam ensemble kanonis yang ada dalam
keadaan i dan berenergi Ui . Keadaan ini akan didefinisikan oleh nilai dari
6N koordinat momentum dan posisi dari N buah partikel sistem yang ada
dalam keadaan ini. Peluang bahwa sebuah sistem ada dalam keadaan i den-
gan energi Ui mungkin diperoleh dengan memperlakukan masing-masing sistem
dalam ensemble seolah-olah mereka merupakan partikel dari sebuah sistem
yang besar. Dalam hal ini, sistem yang besar tersebut adalah ensemble kanonis
itu sendiri yang dipandang mempunyai nilai yang tetap untuk energi dan juga
untuk temperatur. Sekarang, tinjau bahwa sistem-sistem dalam ensemble kano-
nis mempunyai ukuran yang cukup sedemikian sehingga energi interaksi antara
dua sistem dapat diabaikan terhadap energi total sistem-sistem yang berinter-
aksi tersebut. Maka, meskipun interaksi antar partikel dalam sebuah sistem
26.3. SIFAT-SIFAT TERMODINAMIK ENSEMBLE KANONIS 145

mungkin tidak dapat diabaikan, sistem-sistem ini dapat diperlakukan sebagai


partiklel-partikel yang tidak saling berinteraksi dalam sistem yang besar yang
dibentuk oleh ensemble kanonis.

Dengan mengingat kembali hasil yang diperoleh untuk distribusi energi dalam
sistem klasik, dapat dituliskan bahwa peluang sebuah sistem berada dalam
keadaan adalah

pi = p(0)eUi /kT , (26.1)

dengan p(0)
Padalah fungsi dari temperatur ensemble T dan total seluruh peluang
adalah 1, i pi = 1 dengan penjumlahan dilaukan untuk seluruh keadaan i.
Jadi, dapat diperoleh

1
p(0) = P Ui /kT
. (26.2)
i e

Definisikan kembali fungsi partisi untuk sistem-sistem dalam ensemble kanonis


sebagai berikut (sering disebut sebagai constant temperature partition function
of the ensemble)

X
Z= eUi /kT . (26.3)
i

P i /kT
Bandingkan Persamaan (26.3) dengan bentuk Z = ie , penjumlahan
dilakukan untuk seluruh keadaan energi. Dapat dilihat bahwa fungsi partisi ini
mempunyai sifat-sifat yang serupa dengan fungsi partisi total meskipun harus
diingat bahwa fungsi partisi total diberikan untuk kasus sebuah sistem dimana
tidak ada interaksi antar partikelnya.
eUi /kT
Jadi peluang bahwa sebuah sistem berada dalam keadaan i adalah pi = Z .

26.3 Sifat-sifat termodinamik ensemble kanonis

Jika keadaan-keadaan energi dari sistem-sistem dalam ensemble kanonis adalah


Ui maka energi rata-rata sebuah sistem dalam ensemble pada temperatur T
adalah

X
U= pi Ui , (26.4)
i

1 X Ui /kT
U= e Ui , (26.5)
Z
i
146 CATATAN 26. ENSEMBLE KANONIS

kT 2 Z ln Z
U= = kT 2 . (26.6)
Z T T

Hubungan termodinamik antara energi F dan energi U diberikan oleh per-


samaan

 
2 (F/T )
U = T . (26.7)
T

Dari Persamaan (26.6) dan (26.7) dapat diperoleh

(F/T ) ln Z
= k (26.8)
T T

yang mempunyai solusi dalam bentuk

F = kT ln Z + C, (26.9)

dengan C adalah sebarang fungsi yang tidak bergantung pada temperatur.

Sekarang jika hubungan F = U T S dan S = F



T V digunakan maka jelas
bahwa nilai C adalah nol. Jadi energi bebas Helmholtz sebuah sistem dalam
ensembel kanonis sama seperti dalam statistik klasik dan semi-klasis dalam bab-
bab sebelumnya.

Ungkapan untuk fungsi partisi Z dan peluang pi adalah

Z = eF/kT , (26.10)

pi = e(F Ui )/kT (26.11)

Persamaan (26.11) kadang-kadang dipakai sebagai definisi suatu ensembel kano-


nis. Jika peluang sebuah sistem ada dalam keadaan i dengan energi Ui diny-
atakan oleh peramaan pi = e(F Ui )/kT maka sistem tersebut dikatakan sebagai
anggota sebuah ensembel kanonis.

Entropi sebuah sistem dalam ensemble kanonis dapat dievaluasi dari ungkapan
S = F

T V . Dengan menggunakan ungkapan energi Helmholtz F = kT ln Z
maka ungkapan untuk entropi dapat dituliskan sebagai berikut

 
ln Z
S = k ln Z + T . (26.12)
T
26.4. EVALUASI FUNGSI PARTISI TOTAL 147

Tetapi, ungkapan lain untuk entropi dapat diperoleh dengan menggunakan


hubungan lain dalam termodinamik, yaitu

U F
S= . (26.13)
T
P P P
Tulislah energi U sebagai U = i pi Ui dan F = F i pi (ingat bahwa i pi =
1) maka:

X Ui F X F Ui X
S= pi = k pi = k pi ln pi . (26.14)
i
T i
kT i

Ungkapan dalam Persamaan (26.14) untuk entropi ini digunakan oleh beberapa
pengarang sebagai titik awal untuk menyelesaikan semua permasalah mekanika
statistik. Hai itu dilakukan karena sembarang formula mekanika statistik harus
dijustifikasi oleh a posteriori, sering pula dinyatakan bahwa titik awal ini sebagai
postulat yang dapat langsung P diterapkan pada sejumlah besar sistem. Logika
dasar menerapkan S = k i pi ln pi sebagai persamaan fundamental dalam
fisika statistik terletak pada studi teori informasi yang memungkinkan untuk
menghubungkan entropi dengan informasi.

26.4 Evaluasi Fungsi Partisi Total

Dalam bagian ini metoda penjumlahan dan integrasi dalam penentuanfungsi


partisi total akan dibahas untuk kasus klasik dan semi-klasik saja.

26.5 Fungsi Partisi Klasik

Tinjau sebuah sistem pada temperatur T yang terdiri dari N partikel yang
terbedakan dan tidak saling berinteraksi. Jika s , gs , dan Ns masing-masing
menyatakan energi, degenerasi, dan banyaknya partikel pada tingkat ke-s maka
:

P
banyaknya partikel pada tiap tingkat energi dibatasi oleh kondisi s Ns =
N,
P
total energi sistem ketika berada dalam keadaan i adalah Ui = s Ns s

bobot wi untuk keadaan i menyatakan banyaknya susunan yang berbeda


dari partikel-partikel atau banyaknya keadaan mikro sistem dalam
Q  gsNs 
keadaan makro i, yaitu: wi = N s Ns ! .
148 CATATAN 26. ENSEMBLE KANONIS

Bila bobot ini dimasukkan dalam fungsi partisi yang didefinisikan dalam Per-
samaan 26.3 dengan penjumlahan dilakukan untuk semua keadaan i yang
mungkin dari sistem tersebut, maka diperoleh hasil:

!
X X Y  g Ns  X Ns s
s
Z= wi eUi /kT = N exp , (26.15)
i s
Ns ! s
kT
Ns

sekarang, penjumlahan dilakukan


P untuk semua himpunan nilai Ns yang
mungkin yang memenuhi syarat s NS = N . Penjumlahan ini akan mema-
sukkan semua konfigurasi yang mungkin, yaitu semua keadaan makro yang
mungkin.
Q  xNs 
Dengan menggunakan formula ( s xs )N = Ns N s Nss ! dan menuliskan
P P

xs = gs es /kT maka

!N
X X Y  g Ns eNs s /kT  X Y  g Ns 
s /kT s s
gs e = N! = N! eNs s /kT .
s s
Ns ! s
Ns !
Ns Ns

Dan diperoleh hasil bahwa:

!N
X
Z= gs es /kT = ZN
s

s /kT
P
dengan fungsi partisi sebuah partikel dalam sistem s gs e = Z.

26.6 Fungsi partisi semi-klasik

Bila partikel-partikel sistem tak terbedakan dan tak saling berinteraksi, maka
bobot wi untuk keadaan i harus diganti dengan bobot semi-klasik yaitu wi =
Q  gsNs 
s Ns ! . Ungkapan untuk fungsi partisi total menjadi:

!
X X g Ns
Uj /kT s Ns s /kT
P
Z= wj e = e s (26.16)
j s
Ns !

!N
1 X
s /kT ZN
Z= gs e = . (26.17)
N! s
N!
26.7. FUNGSI PARTISI UNTUK KASUS ADA INTERAKSI 149

26.7 Fungsi Partisi untuk Kasus Ada Interaksi

Bila interaksi antar partikel tidak dapat diabaikan, energi partikel akan bergan-
tung pada koordinat posisi dan koordinat momentum. Fungsi partisi diperoleh
dengan mengintegrasikan koordinat -koordinat partikel untuk semua nilai yang
mungkin. Bila d6N adalah elemen volum dalam ruang fasa berdimensi 6N
dari N partikel sistem, maka banyaknya keadaan dalam elemen volum tersebut
adalah d6N /h3N . Besaran ini ekivalen dengan bobot w untuk sebuah sistem
yang partikel-partikelnya ada dalam keadaan tersebut sedemikian sehingga jika
energi sistem dalam keadaan ini adalah U maka:

untuk sistem yang partikel-partikelnya terbedakan


Fungsi partisi klasik:

d6N
Z
Z= eU/kT
6N h3N

untuk sistem yang partikel-partikelnya tak terbedakan


Fungsi partisi semi-klasik:

d6N
Z
Z= eU/kT
6N h3N N !

Perhatikan bahwa untuk partikel-partikel tak terbedakan, ada N ! cara


penyusunan untuk tiga koordinat posisi dan tiga koordinat momentum
dari N partikel sistem yang menyatakan keadaan sistem yang sama. Den-
gan demikian, ada faktor N ! pada penyebut fungsi partisi total pada kasus
ini.

untuk sistem yang partikel-partikelnya tak terbedakan dan tidak ada in-
teraksi antar partikel, maka

1 X 2
pxj + p2yj + p2zj

U=
2m j=1

dan dengan menggunakan d6N = N


Q
j=1 dxj dyj dzj dpxj dpyj dpzj , maka
dengan mudah dapat ditunjukkan bahwa fungsi partisi total adalah
N
1 N (2mkT )3/2 ZN

U/kT d6N
Z=e = = .
h3N N ! N! h3 N!

26.8 Distribusi energi pada ensembel kanonik

Peluang bahwa sebuah sistem pada temperatur T ada dalam keadaan energi Ui
telah ditunjukkan oleh persamaan:
150 CATATAN 26. ENSEMBLE KANONIS

pi = e(F Ui )/kT

Bila (U )dU menyatakan banyaknya keadaan sistem-sistem dalam ensemble


yang mempunyai energi antara U dan U + dU , maka banyaknya sistem yang
mempunyai energi antara U dan U + dU adalah:

N (U )dU = e(F U)/kT (U )dU.

Kurva N (U ) diperlihatkan pada Gambar 26.2. Energi pada nilai maksimum


kurva akan sangat dekat dengan nilai energi rata-rata sistem karena ketajaman
distribusi kurva tersebut sangat menyolok.

Gambar 26.2: Distribusi energi sistem-sistem.

26.9 Aplikasi ensemble kanonis untuk gas tak


ideal

Bagian ini akan membahas satu contoh aplikasi ensemble kanonis untuk gas tak
ideal, yaitu untuk kasus interaksi antar molekul tak dapat diabaikan. Molekul-
molekul yang ditinjau adalah molekul semi-klasik. Energi sistem gas ini diten-
tukan oleh komponen yang bergantung pada momentum dan posisi molekul-
molekulnya. Jika dianggap bahwa energi interaksi antara dua buah molekul
tidak bergantung pada momentum kedua buah molekul yang berinteraksi dan
juga tidak bergantung pada posisi molekul-molekul lain, maka energi total sis-
tem dapat dituliskan sebagai:

N N X
1 X 2  X
U= pxj + p2yj + p2zj + + Uij
2m j=1 j=1 l>j
26.9. APLIKASI ENSEMBLE KANONIS UNTUK GAS TAK IDEAL 151

engan N adalah jumlah molekul, pxj , pyj , pzj adalah komponen dalam arah x,
y dan z dari momentum molekul ke-j. Sedangkan, Ujl adalah energi interaksi
antara molekul ke-j dan molekul ke-l4 dan kondisi l > j dalam salah satu
somasi dari dobel somasi tersebut berlaku untuk semua nilai l untuk mencegah
penghitungan energi interaksi dua kali, misal sekali untuk Ujl dan sekali untuk
Ulj . Fungsi partisi untuk gas semi-klasik ini adalah:

d6N
Z
Z= eU/kT
6N h3N N !


N
1 1
Z X X
p2xj + p2yj + p2zj +

Z= exp Ujl /kT d6N
N !h3N 6N 2m
j=1

l>j

QN
dengan d6N = j=1 dxj dyj dzj dpxj dpyj dpzj .

Integrasi yang melibatkan

(26.18)
152 CATATAN 26. ENSEMBLE KANONIS
Catatan 27

Simulasi: Sistem
Paramagnetik

Sebuah sistem paramegnetik dengan memperhitungkan dua jenis spin atom-


atom, yaitu up (+ 21 ) dan down ( 12 ), dapat disimulasi dengan menggunakan
piranti lunak yang telah akrab dengan peserta kuliah. Piranti lunak yang di-
maksud adalah spreadsheet, misalnya saja Microsoft Excel atau OpenOffice
Calc.

153
154 CATATAN 27. SIMULASI: SISTEM PARAMAGNETIK

Gambar 27.1: Ilustrasi simulasi sistem paramagnetik dengan menggunakan


spreadsheet.
Catatan 28

Soal 1: Tingkat Energi dan


Peluang Termodinamika

28.1 Soal
1. (a) Tabulasikan nilai-nilai bilangan kuantum nx , ny , dan nz untuk dua be-
las tingkat energi terendah dari suatu partikel bebas yang berada dalam se-
buah wadah dengan volume V ! Nilai nx , ny , dan nz = 0, 1, 2, 3, .. (tetapi
tidak boleh semuanya bernilai nol bersama-sama). (b) Bagaimana degen-
erasi dari tiap tingkat (energi)? (c) Hitunglah energi pada tiap tingkat
dalam satuan h2 /(8mV 2/3 ! (d) Apakah tingkat-tingkat energi tersebut
memiliki perbedaan energi yang sama satu sama lain?

2. Hitunglah nilai nj sehingga sebuah atom oksigen yang berada di dalam


suatu kubus dengan rusuk 1 cm akan memiliki energi yang sama seperti
energi terendah yang diperbolehkan bagi sebuah atom helium yang berada
di dalam sebuah kubus dengan rusuk 2 1010 m! Perhatikan bahwa nx ,
ny , dan nz = 0, 1, 2, 3, .. (tetapi tidak boleh semuanya bernilai nol
bersama-sama).

3. Terdapat 30 partikel terbedakan yang terdistribusi pada tiga tingkat energi


yang tidak terdegenerasi, yang dilabelkan dengan 1, 2, 3, sehingga N1 =
N2 = N3 = 10. Energi pada masing-masing tingkat energi adalah 1 = 2
eV, 2 = 4 eV, dan 3 = 6 eV. (a) Bila terjadi perubahan bilangan okupasi
pada tingkat 2, di mana dN2 = 2, tentukan dN1 dan dN3 sehingga
dU = 0! (b) Hitunglah peluang termodinamika keadaan makro sebelum
dan sesudah perubahan terjadi!

4. Lima partikel terdistribusi pada keadaan-keadaan energi dari empat


tingkat energi yang perbedaan jarak energi satu dengan lainnya sama,
1 = 1 , g1 = 1, 2 = 21 , g2 = 3, 3 = 31 , g3 = 4, dan 4 = 41 ,
g4 = 5 sehingga total energi sistem U = 121 . Hitunglah peluang ter-
modinamika untuk setiap keadaan makro dan rata-rata bilangan okupasi

155
156CATATAN 28. SOAL 1: TINGKAT ENERGI DAN PELUANG TERMODINAMIKA

untuk setiap tingkat (energi) apabila partike-partikel tersebut memenuhi


(a) statitstik Bose-Einstein, (b) statitstik Fermi-Dirac, dan (c) statistik
Maxwell-Boltzmann!

28.2 Jawab
1. (a), (b), dan (c) dapat dilihat pada Tabel 28.1, (d) Tingkat-tingkat energi
tidak selalu berjarak sama satu sama lain untuk dua belas tingkat pertama
h2
energi. Tingkat energi ke-6 dan ke-7 berjarak 2 8mV 3/2 sementara yang
 2

h
lain berjarak 8mV 3/2 .

2. Energi sebuah atom oksigen O dalam kubus berusuk 1 cm sama dengan


energi terendah sebuah atom helium He dalam kubus berusuk 2 1010
m, sehingga

O = He
h2 h2
n2j 2/3
= 12 2/3
8mO VO 8mHe VHe
2/3
mO VO mO L2O
n2j = =
mHe V 2/3 mHe L2He d
He
16 (102 )2
= = 1016
4 (2 1010 )2
nj = 108 .

Tingkat energi terendah helium adalah saat n2j = n2x + n2y + n2z = 1.

3. N = 30
3 = 6 eV, g3 = 1, N3 = 10, dN3 =?
2 = 4 eV, g3 = 1, N2 = 10, dN2 = 2
1 = 2 eV, g3 = 1, N3 = 10, dN1 =?
P
(a) dU = j dNj = 1 dN1 + 2 dN2 + 3 dN3
0 = 2dN1 + 4(2) + 6dN3
dN = dN1 + dN2 + dN3
0 = dN1 2 + dN3
2dN1 + 6dN3 = 8
2dN1 + 2dN3 = 4
4dN3 = 4
dN3 1
dN1 + dN3 = 2
dN1 = 2 dN3 = 2 1 = 1
dN1 = 1, dN3 = 1
28.2. JAWAB 157

(b) Sebelum:

Y gjNj g1N1 g2N2 g3N3 110 110 110 30!


Wk = N ! = N! = 30! =
j
Nj ! N1 !N2 !N3 ! 10!10!10! 10!10!10!

Sesudah:
Y gjNj g1N1 g2N2 g3N3 111 18 111 30!
Wk = N ! = N! = 30! =
j
Nj ! N1 !N2 !N3 ! 11!8!11! 11!8!11!

4. N = 5 dan U = 121
Jawab (a), (b), dan (c) dapat dilihat dalam Tabel 28.2, 28.3, dan 28.4.
Catatan:
P
N= Nj = N1 + N2 + N3 + N4 = 0 + 2 + 3 + 0 = 5.
Pj
U = j Nj j = N1 1 + N2 2 + N3 3 + N4 4
= 0 1 + 3 21 + 2 31 + 0 41 = 121
Statistik BE:
Y (Nj + gj 1)!
WBE =
j
Nj !(gj 1)!
e.g.
0! 5! 5! 4!
W1 = = 1 10 10 1 = 100
0!0! 3!2! 2!3! 0!4!
Statistik FD:
Y gj !
WF D = ; gj Nj
j
Nj !(gj Nj )!

e.g.
1! 3! 4! 5!
W1 = = 1161 = 6
0!1! 3!0! 2!2! 0!5!
Statistik MB:
Y gjNj
WMB = N !
j
Nj !
e.g.
10 33 42 50 9
W1 = 5! = 120 1 8 1 = 4320
1! 3! 2! 0! 2

1 X
Nj = Njk Wk

k

e.g.
1
N 1,BE = (0 100 + 0 75 + 1 60 + 1 120 + 2 50 + +2 45)
450
0 + 0 + 60 + 120 + 100 + 90 370
= = = 0.82
450 450
158CATATAN 28. SOAL 1: TINGKAT ENERGI DAN PELUANG TERMODINAMIKA

Tabel 28.1: Tingkat energi, degenerasi, dan keadaan mikro yang mungkin untuk
dua belas tingkat energi pertama partikel dalam kotak tiga dimensi.

Nomor  
h2
Tingkat nx ny nz j 8mV 2/3
gj
Energi
1 0 0
1 0 1 0 1 3
0 0 1
1 1 0
2 1 0 1 2 3
0 1 1
3 1 1 1 3 1
2 0 0
4 0 2 0 4 3
0 0 2
2 1 0
2 0 1
1 2 0
5 5 6
1 0 2
0 1 2
0 2 1
2 1 1
6 1 2 1 6 3
1 1 2
2 2 0
7 2 0 2 8 3
0 2 2
2 2 1
2 1 2
1 2 2
8 9 6
3 0 0
0 3 0
0 0 3
3 1 0
3 0 1
1 0 3
9 10 6
1 3 0
0 1 3
0 3 1
3 1 1
10 1 3 1 11 3
1 1 3
11 2 2 2 12 1
3 2 0
3 0 2
2 3 0
12 13 6
2 0 3
0 3 2
0 2 3
28.2. JAWAB 159

Tabel 28.2: Keadaan makro, keadaan mikro, bilangan okupasi rata-rata dengan
N = 5 dan U = 121 untuk statistik B-E.

Nomor
Tingkat Njk Keadaan Makro k
j /1 gj Nj
Energi j 1 2 3 4 5 6
4 4 5 0 1 0 1 1 2 0.74
3 3 4 2 0 3 1 2 0 1.33
2 2 3 3 4 1 2 0 1 2.10
1 1 1 0 0 1 1 2 2 0.82

Wk 100 75 60 120 50 45 450


Tabel 28.3: Keadaan makro, keadaan mikro, bilangan okupasi rata-rata dengan
N = 5 dan U = 121 untuk statistik F-D.

Nomor
Tingkat Njk Keadaan Makro k
j /1 gj Nj
Energi j 1 2 3 4 5 6
4 4 5 0 1 0 1 1 2 0.77
3 3 4 2 0 3 1 2 0 1.38
2 2 3 3 4 1 2 0 1 1.92
1 1 1 0 0 1 1 2 2 0.92

Wk 6 0 12 60 0 0 78

Tabel 28.4: Keadaan makro, keadaan mikro, bilangan okupasi rata-rata dengan
N = 5 dan U = 121 untuk statistik M-B.

Nomor
Tingkat Njk Keadaan Makro k
j /1 gj Nj
Energi j 1 2 3 4 5 6
4 4 5 0 1 0 1 1 2 0.77
3 3 4 2 0 3 1 2 0 1.39
2 2 3 3 4 1 2 0 1 1.90
1 1 1 0 0 1 1 2 2 0.93

Wk 4320 2025 3840 10800 2400 2250 25635



160CATATAN 28. SOAL 1: TINGKAT ENERGI DAN PELUANG TERMODINAMIKA
Catatan 29

Soal 2: Fungsi Distribusi


dan Entropi

29.1 Soal
1. Tunjukkan bahwa rumusan entropi statistik Bose-Einstein dalam batasan
klasik (gj >> Nj >> 1) akan tereduksi menjadi

X 
gj
 
Sk Nj ln + Nj .
j
Nj

Gunakan aproksimasi Strirling.


2. Fungsi distribusi untuk partikel tak terbedakan dapat representasikan oleh
satu persamaan, yaitu

Nj 1
=   .
gj j
exp kT + a

(a) Apa arti dari ? (b) Apakah nilai dari a untuk (i) statistik Bose-
Einstein, (ii) statistik Fermi-Dirac,
  dan (iii)
 statistik
 klasik? (c) Gam-
Nj j
barkan dalam diagram gj terhadap kT fungsi distribusi Bose-
Einstein, Fermi-Dirac, dan klasik! (d) Dengan merujuk pada pertanyaan
sebelumnya, dalam kondisi apa statistik klasik dapat diterapkan?
3. Misalkan terdapat sebuah sistem dengan N partikel terbedakan. Partike-
partikel tersebut terdistribusi dalam dua tingkat energi tak terdegenerasi.
Bilangan okupasi pada tingkat energi 1 adalah N1 dan partikel-partikel
lain berada pada tingkat energi 2. (a) Tuliskan peluang termodinamika
untuk sistem ini. (b) Tuliskan entropi sistem ini (dalam variabel N , N1 ,
dan konstanta Boltzmann k).

161
162 CATATAN 29. SOAL 2: FUNGSI DISTRIBUSI DAN ENTROPI

29.2 Jawab
1. Statistik Bose-Einstein (BE) memiliki peluang termodinamika untuk su-
atu keadaan makro k
Y (Nj + gj 1)!
WBE = ,
j
Nj !(gj 1)!

di mana dalam limit klasik (gj >> Nj >> 1) suku-suku dalam persamaan
di atas aka menjadi

gj + Nj 1 gj + Nj ,

gj 1 gj .

sehingga

Y (Nj + gj )!
WBE = .
j
Nj !gj !

Kemudian


Y (Nj + gj )!
ln WBE = ln
j
Nj !gj !
X
= [(Nj + gj ) ln(Nj + gj ) Nj ln Nj gj ln gj ]
j
X 
Nj + gj
 
Nj + gj

= Nj ln + gj ln .
j
Nj gj

Kemudian dengan gj >> Nj akan diperoleh

Nj + gj gj

Nj Nj

dan

Nj + gj Nj
=1+ .
gj gj

Dari syarat limit klasik dapat dilihat bahwa Nj /gj << 1 sehingga
 
Nj Nj
ln 1 + .
gj gj

Dengan demikian peluang termodinamika suatu keadaan makro statistik


BE dengan limit klasik akan menjadi
29.2. JAWAB 163

X  
gj
  X 
Nj
 
gj

ln WBE Nj ln + gj = Nj ln + Nj ,
j
Nj gj j
Nj

sehingga entropinya (dengan W) akan mejadi

X  
gj

S = k ln WBE = k Nj ln + Nj .
j
Nj

2. Fungsi distribusi partikel tak terbedakan

Nj 1
=   .
gj j
exp kT + a

(a) adalah potensial kimia setiap partikel.


(b) Nilai a tidak sama untuk ketiga statistik, yaitu
a = 1 untuk statistik BE,
a = +1 untuk statistik FD,
a = 0 untuk statistik MB.
Lihat Sears dan Salinger (Cetakan ke-5, 1980) Gambar 11.11 Ha-
laman 334.
Statistik klasik dapat diterapkan hanya dengan syarat Nj << gj .

3. Sistem partikel terbedakan mengikuti statistik BM.

Tabel 29.1: Sistem partikel terbedakan dengan dua tingkat energi yang masing-
masing memiliki bilangan okupasi N1 dan N N1 dengan degenerasi yang sama
g1 = g2 = 1.

j j gj Nj
1 1 1 N1
2 2 1 N N1

(a) Peluang termodinamika sistem ini adalah

Y gjNj 1N1 1N N1 N!
WMB = N ! = N! =
j
Nj ! N1 !(N N1 )! N1 !(N N1 )!

(b) Entropi sistem adalah


 
N!
S = k ln k ln W = k ln .
N1 !(N N1 )!
164 CATATAN 29. SOAL 2: FUNGSI DISTRIBUSI DAN ENTROPI
Catatan 30

Soal 3: Fungsi Partisi dan


Tabulasi Keadaan Makro

30.1 Soal

1. Sebuah sistem yang terdiri dari N partikel terbedakan terdistribusi dalam


dua tingkat energi yang tak terdegenerasi. Tingkat pertama memiliki en-
ergi nol sedangkan tingkat energi kedua memiliki energi . Sistem ini
berada dalam kesetimbangan termal dengan sebuah reservoir pada tem-
peratur T . Tentukan: (dalam variabel , T , N , dan konstanta Boltzmann
k) (a) fungsi partisi, (b) fraksi dari N1 /N dan N2 /N dari partikel dalam
setiap keadaan, (c) energi internal sistem U (atau E), dan (d) rata-rata
energi sebuah partikel.

2. Enam partikel tak-terbedakan yang mematuhi statistik Fermi-Dirac ter-


distribusi dalam lima tingkat energi, j = (j 1); j = 1, 2, .., 5. Setiap
tingkat energi memiliki degenerasi dengan gj = 3. Energi total sistem
adalah U = 6. (a) Tabulasikan nilai-nilai bilangan okupasi pada setiap
tingkat energi untuk keadaan-keadaan makro yang mungkin bagi sistem
ini, (b) hitunglah peluang termodinamika untuk tiap-tiap keadaan makro,
dan (c) tentukan bilangan okupasi rata-rata dari tiap tingkat energi.

30.2 Jawab

1. Statistik untuk partikel terbedakan menggunakan statistik Maxwell-


Boltzmann (MB)

Jumlah partikel N = N1 + N2 , dengan N1 partikel berada dalam tingkat


energi pertama dan N2 partikel berada dalam tingkat energi kedua.

165
166CATATAN 30. SOAL 3: FUNGSI PARTISI DAN TABULASI KEADAAN MAKRO

Tabel 30.1: Sistem partikel terbedakan dengan dua tingkat energi yang masing-
masing memiliki bilangan okupasi N1 dan N2 dengan degenerasi yang sama
g1 = g2 = 1.

j j gj Nj
1 0 1 N1
2 1 N2

(a) Fungsi partisi setiap partikel


X
Z= gj ej /kT = g1 e1 /kT + g2 e2 /kT
j

= 1e0/kT + 1e/kT = 1 + e/kT .

(b) Fraksi N1 /N dan N2 /N adalah

Nj gj ej /kT gj ej /kT gj ej /kT


=P j /kT
= P j /kT
= .
N j gj e j gj e Z
N1 1 1 h
/kT
i1
= = = 1 + e
N Z 1 + e/kT
/kT /kT
N2 e e 1 h
/kT
i1
= = = = 1 + e ,
N Z 1 + e/kT e/kT + 1

(c) Dari jawab sebelumnya dapat diperoleh


 
N1 h i1
N1 = N = N 1 + e/kT
N

dan  
N2 h i1
N2 = N = N 1 + e/kT
N
sehingga
X
U =E= Nj j = N1 1 + N2 2
j
h i1 h i1 h i1
= N 1 + e/kT 0 + N 1 + e/kT = N 1 + e/kT .

(d) Energi rata-rata tiap partikel

U h i1
= = 1 + e/kT .
N

2. (a) Tabulasi bilangan okupasi pada tiap-tiap keadaan makro yang


mungkin bagi sistem ini, (b) peluang termodinamika tiap-tiap keadaan
makro, dan (c) bilangan okupasi rata-rata tiap tingkat energi dapat dili-
hat dalam Tabel 30.2 berikut ini.
30.2. JAWAB 167

Y gj !
Wk =
j
Nj !(gj Nj )!
3! 3! 3! 3! 3!
W1 =
1!(3 1)! 0!(3 0)! 0!(3 0)! 2!(3 2)! 3!(3 3)!
= 3 1 1 3 1 = 9.
3! 3! 3! 3! 3!
W2 =
0!(3 0)! 1!(3 1)! 1!(3 1)! 1!(3 1)! 3!(3 3)!
= 1 3 3 3 1 = 27.
3! 3! 3! 3! 3!
W3 =
0!(3 0)! 1!(3 1)! 0!(3 0)! 3!(3 3)! 2!(3 2)!
= 1 3 1 1 3 = 9.
3! 3! 3! 3! 3!
W4 =
0!(3 0)! 0!(3 0)! 2!(3 2)! 2!(3 2)! 2!(3 2)!
= 1 1 3 3 3 = 27.
3! 3! 3! 3! 3!
W5 =
0!(3 0)! 0!(3 0)! 3!(3 3)! 0!(3 0)! 3!(3 3)!
= 1 1 1 1 1 = 1.
X
= Wk = 9 + 27 + 9 + 27 + 1 = 73.
k

Tabel 30.2: Keadaan makro, keadaan mikro, bilangan okupasi rata-rata dengan
N = 6 dan U = 6 untuk statistik F-D.
Nomor
Tingkat Njk untuk tiap k
j / gj Nj
Energi j 1 2 3 4 5
5 4 3 1 0 0 0 0 0.123
4 3 3 0 1 1 0 0 0.493
3 2 3 0 1 0 2 3 1.151
2 1 3 2 1 3 2 0 1.726
1 0 3 3 3 2 2 3 2.507

Wk 9 27 9 27 1 73

P P
k Njk Wk Njk Wk
Nj = P = k
k Wk
N11 W1 + N12 W2 + N13 W3 + N14 W4 + N15 W5
N1 =

3 9 + 3 27 + 2 9 + 2 27 + 3 1 183
= = = 2.507.
73 73
N21 W1 + N22 W2 + N23 W3 + N24 W4 + N25 W5
N2 =

2 9 + 1 27 + 3 9 + 2 27 + 0 1 126
= = = 1.726.
73 73
168CATATAN 30. SOAL 3: FUNGSI PARTISI DAN TABULASI KEADAAN MAKRO

N31 W1 + N32 W2 + N33 W3 + N34 W4 + N35 W5


N3 =

0 9 + 1 27 + 0 9 + 2 27 + 3 1 84
= = = 1.151.
73 73
N41 W1 + N42 W2 + N43 W3 + N44 W4 + N45 W5
N4 =

0 9 + 1 27 + 1 9 + 0 27 + 0 1 36
= = = 0.493.
73 73
N51 W1 + N52 W2 + N53 W3 + N54 W4 + N55 W5
N5 =

1 9 + 0 27 + 0 9 + 0 27 + 0 1 9
= = = 0.123.
73 73

X
N= Nj = N1 + N2 + N3 + N4 + N5
j
= 2.507 + 1.726 + 1.151 + 0.493 + 0.123 = 6.
Catatan 31

Soal 4: Distribusi Laju dan


Persamaan Keadaaan

Persamaan-persamaan yang diberikan:

Y gjNj
WMB = N !
j
Nj !

X
Z= gj ej /kT
j

N
Nj = gj ej /kT
Z

r
1
Z
ax2
e dx =
0 4 a3

Y gj !
WFD =
j
Nj !(gj Nj )!

31.1 Soal
1. Lima partikel terdistribusi dalam keadaan-keadaan dengan empat tingkat
energi yang berjarak energi sama satu sama lain: 1 = , g1 = 2, 2 = 2,

169
170CATATAN 31. SOAL 4: DISTRIBUSI LAJU DAN PERSAMAAN KEADAAAN

g2 = 3, 3 = 3, g3 = 4, dan 4 = 4, g4 = 5 dengan energi total U adalah


12. Partikel-partikel tersebut mematuhi statistik Fermi-Dirac statistics.
Hitunglah: (a) peluang termodinamika untuk setiap keadaan makro (dari
keseluruhan 5 keadaan makro), (b) bilangan okupasi rata-rata untuk tiap
tingkat energi, dan (c) entropi sistem.
2. Dalam suatu gas dua-dimensi, molekul-molekul gas dapat bergerak bebas
pada sebuah bidang, akan tetapi dibatasi dalam suatu luas A. Terdapat
suatu sistem dengan N partikel untuk gas monoatomik dua-dimensi pada
temperatur T , massa tiap partikel adalah m. Konsep yang berkaitan
dengan tekanan P menjadi gaya per satuan panjang dan volume V
menjadi luas A. Tentukan: (a) energi tiap partikel, (b) fungsi partisi, dan
(c) persamaan keadaaan gas dari fungsi Helmholtz gas tersebut.
3. Untuk gas ideal monoatomik, jumlah partikel dengan laju antara v dan
v + dv adalah
mv 2
 m 3/2  
2
dNv = N 4 v exp dv
2kT 2kT
(a) Dari persamaan di atas, tentukan ekspresi untuk laju rata-rata kuadrat
v, (b) Telah diturunkan ekspresi untuk tekanan P dari teori kinetik, yaitu
1 2
3 m(N/V )v . Temukan persamaan keadaan dari persamaan ini.

4. Sebuah sistem N partikel mematuhi statistik Maxwell-Boltzmann atau


klasik. (a) Temukan ekspresi untuk rata-rata jumlah partikel pada tingkat
energi j (Nj ) (dalam variable N , T , ln Z, j , dan konstanta Boltzmann
k). (b) Temukan ekspresi untuk energi sistem (U ) (dalam variable N , T ,
ln Z, dan konstanta Boltzmann k).

31.2 Jawab
1. Terdapat lima keadaan makro yang mungkin seperti ditunjukkan dalam
Tabel 31.1. Statistik FD:
(a) Peluang termodinamika masing-masing keadaan makro:
Y gj !
WF D = ; gj Nj
j
N j !(g j Nj )!

2! 3! 4! 5!
W1 = = 1161 = 6
0!(2 0)! 3!(3 3)! 2!(4 2)! 0!(5 0)!
2! 3! 4! 5!
W2 = = 2 3 4 1 = 24
1!(2 1)! 1!(3 1)! 3!(4 3)! 0!(5 0)!
2! 3! 4! 5!
W3 = = 2 3 4 5 = 120
1!(2 1)! 2!(3 2)! 1!(4 1)! 1!(5 1)!
2! 3! 4! 5!
W4 = = 1 1 6 5 = 30
2!(2 2)! 0!(3 0)! 2!(4 2)! 1!(5 1)!
2! 3! 4! 5!
W5 = = 1 3 1 10 = 30
2!(2 2)! 1!(3 1)! 0!(4 0)! 2!(5 2)!
31.2. JAWAB 171

(b) Bilangan okupasi rata-rata tiap tingkat energi:

1 X
Nj = Wk Njk

k

1 264
N1 = (0 6 + 1 24 + 1 120 + 2 30 + 2 30) = = 1.257
210 210
1 312
N2 = (3 6 + 1 24 + 2 120 + 0 30 + 1 30) = = 1.486
210 210
1 264
N3 = (2 6 + 3 24 + 1 120 + 2 30 + 0 30) = = 1.257
210 210
1 210
N4 = (0 6 + 0 24 + 1 120 + 1 30 + 2 30) = = 1.000
210 210
(31.1)

Tabel 31.1: Keadaan makro, keadaan mikro, bilangan okupasi rata-rata dengan
N = 5 dan U = 121 untuk statistik F-D.
Nomor
Tingkat Njk Kead. Mak. k
j / gj Nj
Energi j 1 2 3 4 5
4 4 5 0 0 1 1 2 1.257
3 3 4 2 3 1 2 0 1.486
2 2 3 3 1 2 0 1 1.257
1 1 2 0 1 1 2 2 1.000

Wk 6 24 120 30 30 210

(c) Entropi sistem:


menurut Boltzmann S = k ln Wmax = k ln 120 = 4.787k
menurut Planck S = k ln = k ln 210 = 5.347k

2. (a) Partikel hanya bergerak dalam bidang sehingga hanya terdapat gerak
pada arah-x dan arah-y. Dengan demikian energi tiap partikel adalah

p2x + p2y
=
2m

(b) Fungsi partisi dihitung dengan mengambil bentuk kontinunya, yaitu


Z 2
d 1
Z Z Z
/kT p2x /2mkT
Z= e = e dpx dxdy
h3 h2
A A
= 2 [(2mkT )]2 = 2 (2mkT )
h h
172CATATAN 31. SOAL 4: DISTRIBUSI LAJU DAN PERSAMAAN KEADAAAN

(c) Dengan menggunakan analogi untuk p dapat dicari , yaitu


   
F F
p= =
V T A T

di mana 
A
F = N kT ln Z = N kT ln 2 (2mkT )
h
sehingga
N kT
=
A
Lalu A = N kT merupakan persamaan keadaan yang dicari (mirip
seperti pV = N kT ).

3. Jumlah partikel yang memiliki laju antara v dan v + dv

dNv = N f (v)dV

mv 2
 m 3/2  
f (v) = 4 v 2 exp
2kT 2kT

(a) Rata-rata laju adalah



mv 2
Z Z  m 3/2  
2 4
v2 = v f (v)dv = 4 v exp dv
0 0 2kT 2kT
 m 3/2 mv 2
Z  
= 4 v 2 exp dv
2kT 0 2kT
 m 3/2 3  2kT 5/2
= 4
2kT 8 m
3 2kT 3kT
= =
2 m m

(b) Persamaan keadaaan menjadi

1 1 3kT N kT
p= m(N/V )v 2 = m(N/V ) = pV = N kT
3 3 m V

4. (a) Rata-rata jumlah partikel pada tingkat energi j, yaitu Nj telah


diberikan akan tetapi belum dinyatakan dalam N , T , ln Z, j , dan
konstanta Boltzmann k (tanpa lagi adanya rumusan gj

N X
Nj = gj ej /kT dan Z = gj ej /kT
Z i
Z 1 Z
= gj ej /kT gj ej /kT = kT
j kT j
 
N Z ln Z
Nj = kT Nj = N kT
Z j j
31.2. JAWAB 173

(b) Energi sistem U dalam variable N , T , ln Z, dan konstanta Boltzmann


k
   
F
U = F + TS = F + T
T V
   
F (F/T )
=F T = T 2
T V T V
  
2 N kT ln Z 2
= T =T (N k ln Z)
T T T

= N kT 2 ln Z
T
174CATATAN 31. SOAL 4: DISTRIBUSI LAJU DAN PERSAMAAN KEADAAAN
Catatan 32

Simulasi Keadaan Mikro


dengan Kartu

Keadaan-keadaan mikro dari suatu sistem yang memiliki 28 keadaan mikro dibu-
atkan dalam bentuk kartu dan para peserta kuliah diberikan kesempatan untuk
memilih kartu yang akan keluar dan seorang dari mereka diminta mencatat dan
melakukan tabulasi dari hasil-hasil yang diperoleh.

Gambar 32.1: Ilustrasi 28 keadaan mikro dari suatu sistem yang akan disimu-
lasikan.

175
176 CATATAN 32. SIMULASI KEADAAN MIKRO DENGAN KARTU
177
178 CATATAN 32. SIMULASI KEADAAN MIKRO DENGAN KARTU
179
180 CATATAN 32. SIMULASI KEADAAN MIKRO DENGAN KARTU

Dari hasil yang diperoleh, walaupun hanya dibuat sekitar 4 set, hasil yang diper-
oleh untuk kemungkinan suatu keadaan makro muncul, sudah mendekati. Suatu
hasil yang mencengangkan. Perlu diakan telaah lebih lanjut mengenai hal ini.
Catatan 33

Berkas-berkas

Beberapa kuis dan ujian pada Semester III Tahun 2009/2010 diarsipkan dalam
catatann ini.

33.1 Kuis

Suatu sistem yang terdiri dari 5 partikel mematuhi statistik Fermi-Dirac. Ter-
dapat empat tingkat energi yang diperhitungkan, yaitu 1 = 2, 2 = 3, 3 = 4, dan
4 = 5. Degenerasi masing-masing tingkat energi bergantung dari volume sistem
V dan energi total sistem tergantung dari temperatur sistem T.

Soal 1. Pada titik a dalam ruang parameter V-T, temperatur memiliki nilai
Ta dan volume memiliki nilai Va sistem memiliki energi total Ua = 19 dan
degenerasi tingkat-tingkat energi adalah g1 = 1, g2 = 3, g3 = 4, dan g4 = 6.
Lengkapilah tabel di bawah ini.

Hitunglah entropi sistem Sa dengan menggunakan rumusan Planck.

Soal 2. Pada titik b dalam ruang parameter V-T, temperatur memiliki nilai Tb
= Ta dan volume memiliki nilai Vb Va sistem memiliki energi total Ub = 19
dan degenerasi tingkat-tingkat energi adalah g1 = 1, g2 = 2, g3 = 3, dan g4 =
5. Lengkapilah tabel di bawah ini.

181
182 CATATAN 33. BERKAS-BERKAS

Hitunglah entropi sistem Sb dengan menggunakan rumusan Planck.

Soal 3. Pada titik b dalam ruang parameter V-T, temperatur memiliki nilai Tc
Tb dan volume memiliki nilai Vc = Vb sistem memiliki energi total Uc = 17
dan degenerasi tingkat-tingkat energi adalah g1 = 1, g2 = 2, g3 = 3, dan g4 =
5. Lengkapilah tabel di bawah ini.

Hitunglah entropi sistem Sc dengan menggunakan rumusan Planck.

Soal 4. Pada titik b dalam ruang parameter V-T, temperatur memiliki nilai Td
= Tc dan volume memiliki nilai Vd = Va sistem memiliki energi total Ud = 17
dan degenerasi tingkat-tingkat energi adalah g1 = 1, g2 = 3, g3 = 4, dan g4 =
6. Lengkapilah tabel di bawah ini.

Hitunglah entropi sistem Sd dengan menggunakan rumusan Planck.

Soal 5. Gambarkan keempat titik a, b, c, dan d dalam ruang parameter V-T dan
tentukanlah proses dari titik mana ke titik mana yang mungkin terjadi apabila
hanya entropi sistem yang ditinjau. Apa syaratnya?
33.2. UJIAN 1 183

33.2 Ujian 1

Soal 1. Dengan mengacu pada tabel di atas, a. jelaskan apa yang dimaksud
dengan variabel j, j, gj, k, Njk, N, , Wk dan , b. tentukanlah jumlah partikel
dan jenis statistik dari sistem partikel di atas disertai dengan alasannya.

Soal 2. Terkait dengan Jawab 1.b, a. apakah keadaan makro sistem partikel di
atas telah lengkap? bila belum, lengkapilah, b. hitung pula nilai-nilai Njk, N, ,
Wk dan d. lengkapi pada tabel di atas (tuliskan jalannya pada lembar terpisah,
soal kembali dikumpulkan).

Soal 3. Hitunglah entropi dari sistem partikel di atas dengan menggunakan


a. rumusan Boltzman, nyatakan dengan SBoltzmann, b. rumusan Planck, ny-
atakan dengan SPlanck.

Soal 4. Bila diketahui bahwa S = S1 + S2 dan = 1 2 untuk dua sistem


yang dicampur dalam kesetimbangan termal, a. buktikan hubungan ini dengan
menggunakan rumusan yang diketahui, b. buktikan hubungan ini dengan meng-
gunakan contoh di bawah ini (konstruksi keadaan makro, mikro, dan gabugan
sistem), c. bagaimanakah dengan U, V, dan T? Lengkapi pula tabel berikut ini.
184 CATATAN 33. BERKAS-BERKAS

33.3 Ujian 2

Ketiga statistik diujikan dengan contohnya masing-masing. Diminta untuk men-


jelaskan mengapai terdapat faktor pengali seperti arah polarisasi, arah spin, dan
lainnya.

Soal terakhir membahas mengenai paradok Gibb.

Anda mungkin juga menyukai