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

книги из ГПНТБ / Смирнова Н.А. Методы статистической термодинамики в физической химии учеб. пособие

.pdf
Скачиваний:
151
Добавлен:
27.10.2023
Размер:
20.12 Mб
Скачать

Так как

k T 2

то

дЕ

(Д £ ) 2 = feT2 _ = кГ-Сѵ ; (VI .89)

\дТ ]v,N

дисперсия энергии в системе канонического ансамбля определяется теплоемкостью системы и зависит от температуры. Среднее квадра­ тичное отклонение энергии равно

ѴіЩ* = Т Vk с . (VI.90)

Найдем относительную флуктуацию энергии идеального одноатом­

ного газа, для которого

t — 3/2NkT

и Сѵ = 3lzNk.

Согласно (VI.89)

 

( О Т 2 = - у

Nk*T*\

 

 

У ( Д £ ) 2

. / 2

 

ъ

* = і -

= Ѵ ш -

( V L 9 1 )

В макроскопической системе относительные флуктуации энергии, как правило, очень малы. В некоторых случаях они, однако, явля­ ются значительными, в частности, для гетерогенных систем. Так, если в равновесии находятся две фазы одинакового состава (случай однокомпонентной или азеотропной системы), то подвод тепла к сис­ теме при р = const будет иметь следствием только изменение их масс, но не состава. Состояния фаз (и их температура) будут оставаться не­ изменными, пока присутствуют обе фазы. Поэтому теплоемкость такой гетерогенной системы имеет бесконечную величину. Для системы при заданной общей массе компонентов, заданных р и Т возможны любые значения энергии в интервале между значениями энергии каждой из фаз, флуктуации энергии очень велики, что согласуется с соотноше­ нием (VI.89). Велики также относительные флуктуации энергии вбли­ зи абсолютного нуля. Причина этого становится ясной при рассмот­ рении распределения по уровням энергии в согласии с квантовомеханическими закономерностями.

Формулы, рассмотренные в настоящем параграфе, иллюстрируют

соотношение (1.46) для аддитивных величин [8,.

-=.) . Аналогич-

V х

Ун)

ная зависимость от числа частиц найдена также для относительных флуктуации температуры и плотности, — эти параметры являются нормальными в термодинамическом смысле [см. условие (111.55)1.

VII. НВЛНТОВОМЕХАНИЧЕСКИЕ ПРЕДСТАВЛЕНИЯ

ВСТАТИСТИЧЕСКОЙ ТЕРМОДИНАМИКЕ

§1. Квантовомеханическое описание состояния системы *

Вклассической механике, на которой до сих пор было основано наше рассмотрение, мгновенное состояние системы с F степенями сво­ боды определяется совокупностью значений F обобщенных координат

qv qF (сокращенно q) и сопряженных им импульсов ри pF (сокращенно р). Заданному состоянию системы отвечает точка в фа­ зовом пространстве. Изменение состояния системы во времени опи­ сывается фазовой траекторией, которая, согласно классическим за­ конам движения, однозначно определена при заданных начальных условиях.

Квантовая механика использует понятия классической механики (энергия, координаты, импульсы и др.), но квантовомеханическое опи­ сание состояния системы и ее движения основано на принципах, существенно отличающихся от классических. Классическая механика вытекает из квантовой как предельный случай.

При квантовомеханическом рассмотрении частица наделяется по­ мимо корпускулярных волновыми свойствами (де Бройль, 1924 г.). Согласно принципу де Бройля движение свободной материальной час­ тицы, обладающей импульсом р, связано с распространением моно­

хроматического колебания с длиной

волны

х = ~

'

(VI I .1)

где h — 6,62-Ю- 2 7 эрг-сек-— постоянная

Планка.

Следствием волновых свойств частиц является соотношение неопре­ деленностей Гейзенберга:

A P i A q i > h ;

(VII. 2)

здесь Aqt — неопределенность задания і-й составляющей

координа­

ты; Арі — неопределенность задания сопряженной (і-й) составляющей импульса; h = /г/2тс. Иначе говоря, произведение неопределенностей

задания сопряженных координат и импульса не может

быть меньше

h = hl2iz. Чем точнее фиксируется в опыте положение

частицы, тем

больше становится неопределенность в ее скорости (точная фиксация положения частицы связана с бесконечно большой неопределенностью в задании ее скорости); если же точно задана скорость, то совершенно неопределенным становится положение частицы (случай монохрома­ тической плоской волны). Одновременное точное задание координат и импульсов тела, согласно соотношению Гейзенберга, принципиально невозможно. Для больших и тяжелых тел это соотношение дает лишь пренебрежимо малые поправки к классическому описанию, посколь­

ку величина % очень мала. Однако

в описание систем атомного раз-

* См., например, [32].

 

6—119

161

мера соотношение неопределенностей вносит существенно новые черты

по сравнению с

классическим случаем.

В системе с F степенями свободы соотношение неопределенностей

выполняется для

каждой из F пар канонически сопряженных коорди­

нат и импульсов. Приходится отказаться от классического представ­ ления о том, что все динамические переменные системы (координаты и импульсы, а также функции этих величин) могут быть одновремен­ но определены с любой желаемой степенью точности. Нельзя состоя­ нию системы в данный момент времени сопоставить точку в фазовом пространстве, как это делается в классической механике; можно ука­ зать лишь конечный интервал значений импульсов и координат.

Квантовая теория, в отличие от классической, дает в основном ве­ роятностные предсказания относительно параметров системы в дан­ ный момент времени. Состояние системы с заданным числом частиц

определяется волновой

функцией г|) (q, t), где q

набор обобщенных

координат qi,

qF

Волновая функция, в общем случае

комплекс­

ная, интерпретируется следующим образом: величина г|)* (q,

t)ty(q,

t)dq

пропорциональна

вероятности

того,

что значения

координат для

дан­

ной системы в момент времени

t заключены в интервале от

q до q -f-

- f dq. Если движение

системы финитно (происходит в ограниченном

объеме), то интеграл от

по всем возможным значениям

координат

сходится. Допустим, для волновой функции г|/

dq — К- Можем

ввести функцию

\|) = yp'K~t/2,

для

которой

 

 

 

 

 

J

Ф*4-^

= 1.

 

(ѵп.з)

Выражение (VII.3) носит название условия нормировки волновой фун­ кции, а функции, удовлетворяющие условию (VII.3), называют нор­

мированными.

Квадрат модуля

нормированной

волновой функции

І^І2 =

есть плотность вероятности (величина |ij>|2 dq — вероятность

заданного

значения координат)*.

 

величине М(р, q, t)

В квантовой

механике каждой

динамической

сопоставляется оператор М, который находят, заменяя в выражении М{р, q, t) обобщенные импульсы и координаты на соответствующие операторы: импульсу рк сопоставляют дифференциальный оператор

— — ,

координате qK

— оператор

умножения на координату

qK**.

i dqK

функции

Гамильтона H (р,

q, t)

сопоставляют

оператор

Га­

Так,

мильтона н(—~,

q,

t] —квантовомеханический

гамильтониан.

 

V і

dq

J

 

 

 

 

* Возможно также квантовомеханическое описание на языке функции <р(р, t),

определяющей вероятность

заданного значения

импульсов в момент времени г:

dw(p,t) =tp* (р,<)<р(р, 0 dp-

Функции ty(q, t) и <р(р, t) взаимно связаны.

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

162

Для одной частицы

РІ +

РІ+РІ

Я =

+ И ( г ) .

где г — радиус-вектор частицы;

•v.

h2 I d 2

ô 2

2

\

 

Я

= - 2 ^ Ы +

Ѵ

+ ^ )

+ " ( Г ) -

( Ѵ І Ы )

Первое слагаемое в правой части выражения (VI 1.4)-—оператор ки­ нетической энергии, второе слагаемое—оператор потенциальной энер­ гии*.

В состоянии, описываемом волновой функцией \J), среднее значение

величины М, которой соответствует оператор М, дается

формулой

M = J ф*Мф dq.

(VII.5)

Формула (VII.5) является исходной для расчета измеряемых на опыте физических величин.

Так как оператор импульса есть — — , среднее значение импуль- i dqK

са рк в заданном ^-состоянии определяется как

Среднее значение энергии определяется в виде

В квантовой механике используются только такие операторы физических ве­ личин, которые являются линейными и самосопряженными. Линейными на­ зывают операторы, удовлетворяющие условию

A4 (cj 4-1 + с 2 ф 2 ) = d M 4ft + с 2 M 4>2,

где о и ft постоянные, фі и фг — волновые функции. Требование линейности

оператора связано с одним из основных принципов квантовой механики — прин­ ципом суперпозиции состояний, состоящем в следующем: если система может находиться в состояниях, описываемых волновыми функциями фі и фг, то она может находиться и в состоянии, которое описывается волновой функцией ф = = сіфі + сафг, где ci и с% — постоянные. Оператор является самосопряженным (эрмитовым), если

j ^M^dq = j фі (А< Фі)* Л/-

Требование самосопряженности связано с тем, что средние значения физических величин должны быть вещественными (не мнимыми).

Описание с помощью волновой функции ty(q, t) представляет наи­ более полное описание, возможное в рамках квантовой механики.

* При нахождении оператора кинетической

энергии делается замена рг

h д I h д \

4 д *

2

2

Т а Д Т а Г = ~ *

а н а л о г и ч н о д л я

P* и

V

б*

163

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

д ф

 

ih~~ = H^.

(VI 1.6)

Уравнение Шредингера (VII.6) является одним из основных уравнений квантовой механики и занимает в квантовой механике такое же мес­ то, как уравнения движения в классической механике*.

Если внешнее поле отсутствует или постоянно, то функция Га­ мильтона системы и оператор H явно от времени не зависят; энергия системы постоянна. Состояния, в которых энергия имеет определен­ ное, постоянное значение, называют в квантовой механике стационар­ ными состояниями системы. Стационарные состояния описываются волновой функцией вида

ф(<М) = ф ( д ) Л ( 0 ,

где A (t) = ехр ^ — i E - ^ - j , а функция ty(q) удовлетворяет стационар­

ному уравнению Шредингера

//ф(<7) = £ ф ( < ? ) .

(VII.7)

Так как \А (t)\ = 1, то плотность

вероятности

 

I ФОь ОI 2

= IФ (<?) I2

 

от времени не зависит. Не зависят от времени и средние значения ве­ личин М, операторы которых от времени явно не зависят [см. выра­ жение (VII.5)1. В дальнейшем речь будет идти о стационарных состоя­

ниях и будет рассматриваться

волновая

функция г]з (q).

Интерпретация величины

как

плотности вероятности и, в

частности, вытекающее отсюда условие (ѴІІ.З) налагают определенные ограничения на функцию vj). Требования состоят в следующем: фун­ кция ф должна быть однозначна, конечна и непрерывна. Математи­ ческая задача нахождения возможных состояний системы с заданной энергией Е [возможных зависимостей ty{q)] сводится к решению диф-

* Описание системы с помоіцью волновой функции предусматривает нали­ чие полных сведений о системе и взаимодействиях в ней; если же система на­ ходится во внешнем поле, необходимо знание параметров этого поля в их зави­ симости от времени (все это отражается в операторе Гамильтона системы). Име­ ется аналогия с постановкой задачи в классической механике, когда требуется однозначно описать изменение состояния системы во времени. Разница состоит в том, что в квантовой механике состояние системы в данный момент времени задается волновой функцией <\>(q) и описывается статистически, тогда как в клас­ сической механике состояние определяется совокупностью значений импульсов и координат. Изменения состояния системы во времени однозначно описываются уравнением (VII.6) в квантовом случае и уравнениями движения (11.28) в клас­ сическом. Состояния, описываемые волновой функцией (так называемые «чис­ тые» состояния), представляют, однако, теоретическую абстракцию, о чем под­ робнее см. § 5 этой главы.

164

ференциального уравнения (VII.7) при заданном значении Е, причем требуется найти такие решения, которые удовлетворяют условиям однозначности, конечности и непрерывности. Из математической тео­ рии уравнений типа (VI 1.7) известно, что функция ty, удовлетворяю­ щая названным условиям, в случае если система заключена в конечном объеме, может быть найдена только для определенных дискретных значений Е. Эти значения энергии носят название собственных зна­ чений оператора Гамильтона. Совокупность всех возможных значений ЕиЕ2, ... называют энергетическим спектром системы. Функции tyi(q), удовлетворяющие уравнению (VI 1.7), называют собственными функ­ циями. Задание волновой функции tyi(q) есть определение квантовомеханического состояния системы; энергия системы в заданном кван­

товом

состоянии

фиксирована.

 

Et отвечает

Во

многих случаях

одному собственному

значению

не одна, а несколько линейно независимых

волновых

функций

- ф а ,

\|)j2 , ...

(в этом

случае

любая линейная комбинация функций

tyix,

tyi2, ... также

описывает возможное состояние системы). Число g(

ли­

нейно независимых волновых функций (число

квантовых состояний),

отвечающих

заданному

значению энергии,

называют

кратностью

вырождения

данного уровня.

 

 

 

Одним из основных принципов квантовой механики является прин­

цип неразличимости тождественных частиц*. При описании состояния

системы

частицы обычно условно нумеруются; допустим,

волновая

функция N частиц записывается как г|з(/*і,

rN), где rt

— радиус-

вектор

і-й частицы. Однако перестановка

пронумерованных частиц

не дает нового физического состояния и, следовательно, не должна изменять величины Щ*. Это налагает следующее требование на вол­

новую

функцию

при перестановке пары тождественных частиц

функция г[з либо

остается неизменной (волновую функцию в таком

случае

называют

симметричной), либо изменяет только знак (анти­

симметричная волновая функция).

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

как предельные закономерности. Во многих случаях

квантовые члены

дают

лишь небольшие поправки к классическим

уравнениям

(для

 

 

h

 

« 1 ,

частицы, движущейся одномерно, это имеет место, если —

d p ]

dx

 

 

 

Р2

 

 

*

Возможность различить частицы предполагает, что мы можем

проследить

за движением каждой из частиц в отдельности,т. е. наблюдать траекторию час­ тицы. Таким является классическое описание движения. Однако в квантовой механике в силу принципа неопределенности понятие траектории частицы ли­ шено смысла. Чтобы «пометить» частицы, мы должны определить их положения в какой-то момент времени. Но если при измерении точно определены координаты частиц, то совершенно неопределенными оказываются скорости. Поэтому не­ определенными будут положения частиц после измерения. Повторив измерение координат частиц через какой-то промежуток времени, мы не сможем сказать, какая именно частица находится в заданной точке пространства. Частицы как бы теряют свою индивидуальность.

165

т. е. если постоянная h очень мала по сравнению с величиной р г

dx

Надо сказать, что принципиальную роль играет наличие постоянной h в квантовых уравнениях; фундаментальное значение этой величины проявляется, в частности, через соотношение неопределенностей. Если положить h = 0, квантовые уравнения движения переходят в строго классические. Широкая область применимости классической механики обусловлена тем, что обычно рассматриваются объекты и процессы, по сравнению с масштабами которых постоянная h — очень малая величина.

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

§ 2. Квантовые состояния некоторых простых систем

Частица в прямоугольном потенциальном ящике. Рассмотрим сво­ бодное движение частицы внутри ящика кубической формы с идеальноотражающими стенками. Внешнее поле внутри ящика отсутствует и потенциальная энергия частицы постоянна. Примем, что внутри ящика и(х, у, z) = 0. Стенки ящика представляют потенциальный барь­ ер бгсконечной высоты, так что на стенках происходит скачок потен­

циала

от и = 0 до и — со. Поэтому вероятность нахождения частицы

вне ящика равна

нулю; вне ящика = 0.

Найдем допустимые зна­

чения

энергии и

собственные функции

частицы, движущейся

внутри куба, длина ребра которого равна l(V — I 3 ) . Масса частицы т.

Уравнение Шредингера (VII.7) для

частицы в отсутствие

внешнего

поля [в выражении (ѴІІ.4) положим и = 0] запишется в виде

 

Ъ? / О Н

ô H

Ô H \

^ Н

 

 

-^(^ + ^ +

( V I L 8

(энергию частицы

обозначим через г ) .

Начало

координат

совместим

с одной из вершин

куба, оси х,

у, я z

направим

по ребрам.

В такой

системе координат граничными условиями задачи являются следую­ щие:

ф ( 0 , у,

г)

= 0 ;

ф(1,

у,г)

= 0 ;

ф(лг, 0,

2) = 0 ;

і<(х,

/ , 2 ) = 0 ;

(VII.9)

</.

0)

= 0 ;

 

у,

/ ) = 0

;

[поскольку функция ty(x,y, z) непрерывна и равна нулю вне ящика, она должна быть равна нулю и на стенках ящика]. Задача, таким образом, сформулирована математически и состоит в отыскании ре­ шения дифференциального уравнения (VI 1.8), удовлетворяющего гра­ ничным условиям (VI 1.9).

Решение будем искать в форме

«К*, у, z) = Ы * ) Ы < / ) < Ы г ) . (VII.10)

166

Подставив

( V I I . 10) в уравнение (VII.8) и поделив обе части

уравнения

на ty(x, у ,

г), найдем

 

 

 

 

 

 

 

± * ^

+

± J ^ +

±

=

 

( V I I Ц)

 

Фі

àx*

^

ф2 ду*

^ ф3

<?22

h2

У • >

Три слагаемых

в левой

части равенства

( V I I . 11)

зависят

соответст­

венно только от X, у и z. Сумма трех независимых функций равна постоянной лишь в том случае, если каждая из функций постоянна. Следовательно,

фх дх* ~~ *' фа ду* - *«* ф3 dz2 г > l V U>

где еж, Ау, kz — некоторые постоянные, с которыми энергия

частицы

связана

следующим соотношением:

 

 

 

 

 

 

 

 

 

 

 

—Ê(4

 

 

+ *y + Q-

 

 

 

<ѴІІЛЗ>

Согласно ( V I I . 12)

 

аа Фі (*)

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

fc2

 

(ж) - 0.

 

 

 

 

 

 

 

 

дх2

+

k\ фх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение этого уравнения можно записать в форме

 

 

 

 

 

 

 

фі (X) = Аг sin (ед. x)

+

ВІ cos

д;).

 

 

(VI1.14)

Из граничных условий (VII.9)

и

выражения ( V I I . 10)

вытекает,

что

 

 

 

 

Фх(0) =

0

 

и

ф г ( / ) = 0 .

 

 

 

 

Чтобы удовлетворить первому из написанных условий,

требуется

в

выражении

(VI 1.14) положить Вг

 

=

0.

Второе условие

будет выпол­

нено, если

 

kxl=nxn,

где

rtjf

= 1, 2, 3, ...

 

 

(VII.15)

 

 

 

 

 

Постоянная

kx =

nxrjl,

следовательно, принимает дискретные зна­

чения. Волновая

функция tyi(x)

определена

выражением

 

 

 

 

 

 

ФІ (X) = Ai

sin i

ç

^ - X

j.

 

 

( V I I . 16)

Аналогичным образом

можно

показать, что

 

 

 

 

 

 

 

 

А у = - 2 і 1 .

 

и

 

 

 

 

 

(VII.17)

где пу

и nz

X целые

положительные

числа. Для

волновой

функции

ij)(x, у , г) после подстановки в

(VI 1.10)

выражений

для

сомножителей

<Ы*). Ыі/)>

Ыг)

получим

 

 

 

 

 

 

 

 

 

 

 

ф (X,

у, z)

= A sin (2Ç-x}

 

sin

(JbÇLy ) sin ( - ^ J L г j ,

(vii.18)

где nx, ny, пг — целые положительные числа. Таким образом, квантовомеханическое состояние частицы в потенциальном ящике опре-

167

деляется тремя числами соответственно трем степеням свободы час­

тицы. Квантовые числа

пх,

пу,

nz

можно

рассматривать как

состав­

ляющие вектора я = (пх,

 

nv, nz).

дискретный.

Учитывая

 

условия

Энергетический

спектр

частицы

 

( V I I . 15) и (VII.17)

из выражения

(VII.13), находим

 

 

£

= —

(

ni

+ ni

+ ni) =

- Д г

г

"2 -

 

(VII.19)

 

8m;2

v *

y

г /

8mV / s

 

 

v

'

Большинство уровней энергии вырождено. Кратность вырождения уровня равна числу наборов чисел пх, пу, пг, которым отвечает один и тот же модуль вектора п. Так, одну и ту же энергию имеют состоя­ ния со следующими значениями квантовых чисел:

пх

пу

nz

1

1

2

1

2

1

2

1

1

(для всех трех состояний п 2 = 6).

Одномерный гармонический осциллятор. Энергетический спектр одномерного гармонического осциллятора может быть найден при решении уравнения Шредингера, являющегося аналогом классичес­ кого уравнения

Н(р, х) =

Уравнение Шредингера

для одномерного гармонического осциллятора

запишется в виде

 

 

 

 

 

 

/

д2

1

Л

,

 

 

+

 

к х 2

ф =

еф

V

дх*

2

/

т

т

(выражение в скобках есть оператор Гамильтона для одномерного

гармонического осциллятора)

или

 

 

дЧ

I

1

\

— - +

е —

кх2

ф = 0.

дх*

\

2

у

Квантовомеханическое состояние одномерного гармонического осциллятора определяется заданием одного квантового числа (ѵ). Уровни энергии определены формулой

 

Е

„ = А ѵ ( 0 + -і-|,

(V1I.20)

где v— частота

колебания,

ѵ — любое целое

неотрицательное

число

(ѵ = 0, 1, 2...).

Как видно из выражения (VII.20), расстояние

между

любыми соседними уровнями равно h (уровни энергии эквидистантны); все уровни энергии не вырождены. Наинизшему энергетическому состоянию осциллятора = 0) отвечает энергия Аѵ/2, которую мож­ но назвать нулевой энергией осциллятора (заметим, что согласно клас-

168

сической теории наинизшее возможное значение энергии осциллятора соответствует покою и равно нулю).

Ротатор. Ротатор, по определению, представляет систему, прост­ ранственное положение которой в любой момент времени полностью определяется двумя углами Ѳ и <р (например, материальная точка, связанная с неподвижным центром невесомым жестким стерж­ нем, так что движение ее происходит по сфере). Классическое опи­ сание движения ротатора см. гл. IV, § 5. Движение классического ротатора, на который не действуют какие-либо внешние силы, есть вращение в плоскости с неизменной угловой скоростью, так что вектор момента количества движения M постоянен.

Для изолированного квантового ротатора (ротатора в стационарном

состоянии) определены величина M

момента количества

движения

(а следовательно, и энергия s = МУ2І)

и проекция Мг

вектора M

на фиксированную ось г. Квантовомеханическое состояние ротатора (волновая функция) характеризуется двумя целыми числами / и т, где / может принимать все целые неотрицательные значения от 0 до

со, a m принимает значения—/, —/ +

1, ...,0

/ — 1 ,

/ (каждому

зна­

чению /

отвечает

2/ + 1 значений числа т).

Число /

определяет

ве­

личину

момента

количества движения:

 

 

 

 

 

М- = Vi U +

1)

ft-

(VII.21)

Число m характеризует величину проекции момента количества дви­ жения на фиксированную ось г:

Мг = hm.

Энергия ротатора зависит только от / и определяется выражением

е, = M2 = / ( / + 1 ) h2 . (VII.22)

Все уровни энергии, кроме наинизшего (/ = 0; е0 = 0), вырождены;

степень вырождения

определяется

числом значений т, возможных

при заданном /, и

равна gj = 2/

+ 1.

§ 3. Число квантовых состояний для заданного интервала значений

энергии. Квазиклассическое приближение

Число собственных состояний, которым отвечает значение энер­ гии частицы в интервале от е до г + Д г , обозначим

A Q = с (e) A e,

(VII.23)

где с(г) = ДО/Де — энергетическая плотность

состояний.

Определим энергетическую плотность состояний для частицы, дви­ жущейся в потенциальном ящике объема V, свободном от действия внешних сил. Мы показали в предыдущем параграфе, что состояние

169

Соседние файлы в папке книги из ГПНТБ