Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Лабораторный практикум по молекулярной физике

.pdf
Скачиваний:
43
Добавлен:
24.03.2015
Размер:
1.49 Mб
Скачать

12

D

ABC

1

4

16

4/16

13

C

ABD

 

 

 

 

14

B

ACD

 

 

 

 

15

A

BCD

 

 

 

 

16

-

ABCD

0

1

16

1/16

Σ

 

 

 

16

 

1

Легко заметить, что число способов размещения монет в данном состоянии (“ О”) равно числовому коэффициенту в разложении двучлена (бином Ньютона):

(p+q)4=1p4+4p3q+6p2q2+4pq3+1q4.

Как известно, биноминальные коэффициенты такого разложения равны числу сочетаний из 4-х элементов по k=4,3,2,1,0 (проверьте данные таблицы П1):

C4k = 4! k!(4 − k )!

Каждое из слагаемых бинома Ньютона определяет

вероятность конкретного состояния, т.е.: w(k 4 , p ) = C 4k p k q 4 − k .

Если эта формула справедлива для n=4, то она применима и для n =5,6, …, и т.д. (метод математической индукции). Тогда для n монет это выражение примет следующий вид: w(kn , p) = Cnk p k q n k .

Это соотношение называется биномиальным распределением. Оно определяет вероятность того, что при n испытаниях интересующее нас событие (выпадение монет "орлом" вверх) произойдет ровно k раз.

Применительно к физической системе, которая может находиться в двух состояниях, биномиальное распределение позволяет вычислять вероятность данного состояния.

ПРИЛОЖЕНИЕ 2

22

mk=np.

Характеристики биномиального распределения 1. Математическое ожидание

Найдем, чему равно математическое ожидание для биномиального распределения. Подставляя выражение (1) в формулу (2), получим:

m k =

n

 

 

n !

 

q n k =

0 +

 

n !

 

p q n − 1 +

k

 

 

 

p k

 

 

 

k !(n k )!

 

 

1)!

 

 

k = 0

 

 

 

 

 

1 !(n

 

+

 

n !

 

 

p 2 q n − 2

+

... +

n !

0 .

 

 

 

 

 

 

 

 

p n q

 

 

2 ! (n

2 )!

 

 

 

 

 

 

 

 

 

n !0 !

 

 

 

Это выражение можно упростить, используя свой-

ства факториала: n!=n(n-1)!, n!=n(n-1)(n-2)!, ...

Тогда будем иметь:

m k

=np

= n p q n − 1 + n ( n − 1 ) p 2 q n − 2 + . . .+ p n =

 

 

n −1

+(n −1) pq

n −2

+...+ p

n −1

=np(q+ p)

n −1

=np.

q

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь учтено, что (q+p)n-1=1 в силу того, что q+p=1. Таким образом, искомое выражение для математического ожидания биномиального распределения имеет следующий вид:

2. Дисперсия

Определим дисперсию биномиального распределения, воспользовавшись её математическим определением (4):

2

n

 

 

 

 

2

k

 

 

 

 

 

 

σ k =

(k m k

) w

 

, p ,

 

 

 

 

 

 

 

 

 

k

k = 0

 

(

 

 

 

n

 

 

 

)

 

k 2

2 k m

k

+

k

σ 2 =

 

 

m

2

W

k

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k = 0

 

 

 

 

k

 

 

 

 

 

n

 

 

n

 

 

 

 

 

2

n

− 2 m k kW

 

, p + m k

 

 

k = 0

 

 

n

 

 

 

 

 

k = 0

 

=

n

2

k

 

, p

k

 

W

 

, p

 

 

 

 

k = 0

 

n

 

 

 

k

 

 

 

 

 

 

W

 

, p .

 

 

n

 

 

 

 

 

 

 

 

 

Первое слагаемое в этом выражении представляет

23

собой математическое ожидание величины k2, сумма во втором слагаемом является математическим ожиданием величины k, а сумма в третьем слагаемом представляет собой полную вероятность того, что из n опытов "орел" выпадет 0, 1, 2, ..., n раз. Очевидно, что эта вероятность должна быть равна 1, так как выпадение из n опытов любого числа "орлов" является достоверным событием:

n

w (k n , p ) = 1.

k = 0

Это выражение называется условием нормировки. Учитывая сказанное, выражение для дисперсии можно записать следующим образом:

σ 2 = m

2 − 2m2

+ m2

= m

2

m2

= m

2 − (np)2 .

k

k

 

 

k

k

k

 

k

k

 

Для расчета m

2

воспользуемся тождеством:

 

 

 

k

k2 = k2 + k – k = k(k-1)+ k.

 

 

 

 

Тогда

m k 2

= m k ( k −1 )

+ m k

= m k ( k −1 )

+ np .

Подставляя это соотношение в выражение для дисперсии, получим:

σ k2 = m k ( k − 1 ) + np − ( np ) 2 .

Применяя определение математического ожидания, но уже для функции k(k-1), будем иметь:

m k ( k − 1 ) =

n

 

 

n !

 

p k q n k =

 

k ( k

1)

 

 

 

k !( n

 

 

 

 

k = 0

 

 

 

k ) !

 

 

n ! p 2 q n − 2

 

 

n ! p 3 q n − 3

 

 

n ! p n

 

= 0 + 0 + 2

 

 

+ 3 × 2

 

 

 

+L +n (n - 1)

 

=

 

 

 

3!(n - 3) !

n ! 0!

 

2 !(n - 2) !

 

 

 

 

 

=n(n −1) p2qn−2 + n(n −1)(n − 2) p3qn−3 +L+ n(n −1) pn =

=n(n − 1) p2 [qn−2 + (n − 2) pqn−2 + L+ pn−2 ]= n(n − 1) p2 ,

так как сумма в квадратных скобках равна единице:

q n − 2 + ( n − 2 ) p q n − 2 +L + p n − 2 = ( q + p ) n − 2 = 1 .

Таким образом, выражение для дисперсии биноми-

24

ального распределения получим в следующем виде:

σ 2k = n(n − 1) p2 + np − (np) 2 = (np)2 np2 + np − (np)2 = np(1 − p) .

 

Учитывая, что

1-p=q, окончательно

будем иметь:

σ k2

= npq .

 

 

 

 

 

 

ПРИЛОЖЕНИЕ 3

 

Вывод распределения Лапласа-Гаусса

 

Преобразуем биномиальное распределение, пользу-

ясь

формулой

Стирлинга

для

факториала:

n ! ≈

 

 

 

 

 

 

 

 

 

n n e n , которая справедлива при n >10.Тогда

 

n

 

 

 

 

 

 

 

 

k

 

 

 

n!

k

 

 

nk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

, p =

 

 

 

 

 

 

 

 

 

 

p

 

q

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

k!(n k )!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n n e n

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 π n

 

 

 

 

 

 

 

p k q n k

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 π k k k e k

 

 

2 π ( n k ) ( n k ) n k e − ( n k )

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n n k + k

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

p k q n k

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k ) k k

 

 

 

 

 

 

k ) n k

 

 

 

2

π

 

 

 

 

 

k ( n

 

( n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

np

 

k

 

 

 

 

n k

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

nq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

π

 

 

 

k ( n

k ) k

 

 

n

k

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем переменную u = k-np, и. заменив k через u:

k = u+np, будем иметь:

 

 

 

 

 

 

 

 

u+np

 

 

 

nunp

k

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

np

nq

w

 

 

, p

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(u + np)(n u np) u + np

 

n u np

 

Имея в виду, что p = 1- q, вынося за скобки np и nq и произведя далее очевидные математические действия, придем в виду:

k

 

=

w

 

, p

 

n

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n pq

 

+

u

u

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

np

 

nq

+1

u −(np+u)

np

−1

u −(nqu)

.

nq

Если n→∞, то u/np<<1 и u/nq<<1. Поэтому в вы-

25

ражении под радикалом можно оставить только единицу. Этого делать нельзя по отношению к двум последним сомножителям, так как они возводятся в более высокие степени.

Прологарифмируем последнее выражение:

k ln w

n

 

 

 

 

 

u

 

 

 

u

 

 

 

 

 

, p

= − ln

npq − (np + u) ln1

+

 

 

− (nq u) ln1

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

np

 

 

nq

Разложим логарифмические функции в ряд по

формуле:

ln (1 + x ) = x

x 2

+

x 3

−L

 

 

 

2

3

 

Ограничиваясь первыми двумя слагаемыми, полу-

чим:

 

 

u

 

u

 

1

u

2

,

 

u

+

=

 

ln 1 −

 

 

 

ln 1

 

 

n p

 

 

 

 

 

 

 

nq

 

 

n p

 

 

2

n p

 

 

 

 

= −

u

1

 

u

 

2

 

 

 

 

 

.

 

 

 

 

 

nq

 

2

 

 

 

 

 

 

 

nq

 

 

Подставив эти выражения в первоначальную формулу и произведя алгебраические преобразования, учтя при этом, что p+q =1, получим:

ln

w (k n , p

)

= −

u 2

.

 

 

 

 

 

 

npq

2 npq

Отсюда, возвращаясь к старой переменной и потенцируя, получим искомое выражение:

 

 

 

( k np ) 2

 

 

 

 

2npq

 

w(k n , p ) = 1 2πnpq × e

 

.

 

 

Это и есть распределение Лапласа-Гаусса. Поскольку при выводе этой зависимости использовалась приближенная формула Стирлинга для факториалов, которая справедлива только при n >10, то и распределение Лапла- са-Гаусса также справедливо только при n >10. При n→∞ распределение Лапласа-Гаусса совпадает с биномиальным распределением.

ПРИЛОЖЕНИЕ 4

26

Вывод распределения Максвелла

Термодинамическое равновесие устанавливается в результате громадного числа столкновений между молекулами. В результате каждого столкновения проекции скорости молекулы испытывают случайное изменение Δυx, Δυy, Δυz, причем изменения каждой проекции скорости независимы друг от друга.

Рассмотрим движение молекулы, скорость которой в начальный момент времени равна нулю. Изменение проекций ее скорости при i-том столкновении с другими молекулами обозначим Δυxi, Δυyi, Δυzi. По прошествии достаточно большого промежутка времени, в течение которого молекула испытала громадное число столкновений, ее проекции скорости станут равны:

υx = υxi ,

υy = υyi ,

υz = υzi .

i

i

i

Каждая из этих проекций является суммой большого числа случайных величин, для которых справедлив закон Гаусса (см., например, [1], §2). Следовательно, плотность распределения вероятности для проекций скорости можно представить в следующем виде:

ϕ(υ2 )= Ae−αυx2

, ϕ(υ2 )= Ae−αυy2

, ϕ(υ2 )= Ae−αυz2

, (П1)

x

y

z

 

где А и α - постоянные, одинаковые для всех трех проекций ввиду полной эквивалентности осей координат и независимости случайных величин υx, υy, υz .

Вероятность того, что проекция скорости на ось x заключена в интервале [υx,υx+dυx], равна:

dP(υ

x

) = ϕ (υ 2 )dυ

x

= Ae−αυ x2 dυ

x

.

(П2)

 

x

 

 

 

Аналогичные формулы справедливы и для других проекций скорости. Вероятность того, что скорость моле-

27

кулы заключена в интервале [υx, υy, υz, υx+dυx, υy+dυy, υz+dυz], по закону умножения вероятностей выражается в виде:

dP(υx y z )= A 3e[α(υx2 y2 z2 )]dυxdυydυz . (П3)

Постоянная А находится из условия нормировки:

∞ ∞ ∞

) = 1.

 

∫ ∫ ∫ dP(υx y z

(П4)

− ∞ − ∞ − ∞

Это условие означает следующее: вероятность того, что скорость молекулы может быть любой внутри интервала (-,), равна единице, то есть это событие является абсолютно достоверным.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−αυ x2

dυ x =

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая,

что

 

 

e

 

,

 

из условия нор-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− ∞

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

мировки (П4) находим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A =

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(П5)

 

 

 

 

 

Вычислим среднюю кинетическую энергию моле-

кулы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mυ

2

 

 

m

 

 

 

 

 

 

 

 

m

 

∞ ∞ ∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

υx2 y2 z2 =

 

∫ ∫ ∫ (υx2 y2 z2 )dP(υx y z ) =

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

−∞−∞−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

α

3 2 ∞ ∞ ∞

2

2

 

2

 

 

 

 

 

 

2

 

 

 

 

2

 

 

2

 

 

 

 

 

 

=

 

 

 

 

 

 

∫ ∫ (υ x

+ υ y

+ υ z

)exp [− α (υ x

+ υ y

+ υ z )]dυ x dυ y dυ z .

2

 

 

 

 

 

π

 

− ∞ − ∞ − ∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегралы такого типа находятся следующим обра-

зом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

π

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

υx2e−αυx dυx

=

 

 

 

e−αυx dυx = −

 

 

 

 

 

 

=

 

 

a 2 π .

 

 

 

∂α

 

∂α

a

2

 

 

−∞

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Применяя

этот

 

прием

к

 

каждому

интегралу, полу-

28

чим: mυ2 = 3m .

2 4a

Используя связь между средней кинетической энергией молекул и температурой: mυ22 = 32 kT , получим выраже-

ние для α: α = m 2 kT .

С учетом полученного выражения и формулы (П5) уравнение для вероятности (П3) примет вид:

dP (υ x y z ) = (m kT )

 

m (υ x2 y2 z2 )

dυ x dυ y dυ z . (П6)

e

2 kT

3 2

 

 

 

От распределения проекций скорости перейдем к распределению модуля скорости. Для этого, поскольку распределение скоростей изотропно, нужно перейти в сферическую систему координат в пространстве скоростей и

проинтегрировать уравнение (П6) по сферическому слою толщиной du, радиус которого равен u=(ux2+uy2+uz2)1/2.

Тогда имеем: duxduyduz=u2ddu,

где dW - телесный угол, под которым из начала координат виден элемент поверхности сферического слоя. Очевидно, интеграл по всей поверхности сферического

слоя равен:

υ 2dΩ = υ 2 dΩ = 4πυ 2 .

Ω = 4π Ω = 4π

Поэтому в результате интегрирования в (П6) по сферическому слою толщиной du получаем вероятность нахождения модуля скорости молекулы в интервале скоро-

стей [u, u+du]:

mυ 2

dP(u) = 4p (m2pkT )32 e 2 kT u 2 du .

Функция

29

mυ 2

 

f (υ ) = 4π (m kT )3 2 e

 

υ 2

 

2kT

(П7)

называется плотностью распределения вероятностей Максвелла. Поскольку молекулы движутся независимо и случайно, число молекул dn(u), скорости которых заключены в интервале [u, u+du], равно:

d n(u)=n×dP(u),

где n - полное число молекул в системе. Относительное число молекул в интервале

[u, u+du] равно: dn(u ) = dP(u ) = f (u )× du , или, с учетом n

(П7), получим закон Максвелла распределения молекул по скоростям:

 

n (υ )

 

(m 2 π kT )3 2

m υ 2

 

 

d

=

e 2 kT υ 2 d υ .

(П8)

n

 

 

 

 

 

 

 

30

ОПРЕДЕЛЕНИЕ ПОСТОЯННОЙ БОЛЬЦМАНА

1. Цель работы

1.1. Изучить экспериментальный метод определения постоянной Больцмана.

1.2. Определить величину постоянной Больцмана и оценить погрешность измерения.

2. Краткое теоретическое ведение

Согласно современным представлениям, все вещества состоят из молекул и атомов, которые находятся в непрерывном хаотическом движении. Между частицами вещества действуют силы взаимного притяжения и отталкивания, убывающие обратно пропорционально расстоянию в определенной степени. В зависимости от интенсивности хаотического движения и сил взаимодействия между молекулами, которые зависят от внешних условий (температуры, давления), различают три агрегатных состояния вещества: твердое, жидкое, газообразное. С изменением внешних условий наблюдается переход вещества из одного состояния в другое – фазовый переход первого рода. В газообразном состоянии тепловое движение наиболее интенсивно, а силы взаимодействия малы. Пренебрегая взаимодействием молекул и размерами, получаем модель идеального газа как совокупности материальных точек, не взаимодействующих друг с другом и находящихся в непрерывном хаотическом движении. Ясно, что такая модель газа выбрана для простоты расчета. Реальный газ, конечно, не удовлетворяет таким представлениям, однако, если газ находится под небольшим давлением, то по своим свойствам он приближается к идеальному. Несмотря на существенные упрощения, представления об идеальном газе позволяют находить зависимость макроскопических параметров

39