где С0о и SC,» определяются так, чтобы |
для |
бесконечного |
ато |
ма энергия |
и поправка к ней |
обращались |
в нуль. |
При |
этом- |
потенциал |
атома и поправка |
к нему |
вычисляются |
по |
фор |
мулам: |
|
|
|
|
|
|
|
eU (г) = р [Ф(х) х - 1— 1]; |
еШ (г) = |
[6Ф (*) х~1— 6Ф (1)]; |
|
6Ф (х) = |
{[*Ф (Х)]'7- -ь |
|
Ф (*)}• |
|
|
Разлагая выражения (41.8) в ряды по степеням параметра ос, получаем асимптотические выражения для случая сильного* сжатия:
ФМ |
|
|
|
■4“* + |
4 4 + 0 (<*!); |
|
ф* = |
т а + > |
, + 0 ( а '): |
|
ф (х) — х + О(а2); |
а = г0 |
(3Z)‘/j |
Г |
1 |
— гО |
(3Z)'!' |
а |
ха0 |
|
20ха0 |
' |
|
L |
|
Подставляя эти выражения в формулы (41.8), получаем асимп тотические выражения и для термодинамических величин си стемы. Не будем их здесь выписывать ввиду громоздкости по лучающихся формул [4].
Ранее уже отмечалось, что уравнение Т—Ф универсально в том смысле, что его нужно решить один раз, получив решения
|
|
|
|
|
|
для любых Z изменением масштаба. |
Такие |
уравнения назы |
ваются а в т о м о д е л ь н ы м и . Если |
за неизвестные |
перемен |
ные принять ZV и |
р-1, а за |
функции от |
них |
величины |
2- 10/а р у Z -4/» р, Z-7/a Е, |
то легко |
видеть, что и более сложные |
уравнения, рассматриваемые здесь, не содержат Z в комбина циях, отличных от приведенных. Если же дополнительно вве сти функции Z-"/. бР, Z~b/s бЕ, Z_2/a р, то автомодельность распространяется и на уравнения для поправок первого по рядка. Это обстоятельство сильно облегчает численные расче ты, так как в этом случае любое решение уравнений годится для всех элементов. Напомним, что уравнение Т—Ф—Д необхо димо решать отдельно для каждого Z.
Численное интегрирование уравнений (41.8) выполнено в работе [4]. Интересно сравнить полученные результаты по сжи-
маемости железа с экспериментальными данными [1]. На рис. 56 и 57 представлены теоретические и экспериментальные
кривые давления и энергии в функции от |
сжатия |
п/п0, где |
п0— нормальная плотность вещества. Как |
видно |
из |
рисунков, |
вклад квантовых эффектов в статистической теории |
системы |
сильно сжатых атомов является существенным. |
Пунктиром |
показана экстраполяция экспериментальных кривых в область, где поправка к давлению, вычисленному в приближении Т—Ф, составляет 30%. Эти экстраполированные кривые, переходя-
Рис. 56. Теоретическая и эксперимеи- |
Рис. 57. Теоретическая и эксперимен |
тальная записимостьдавления железа |
тальная зависимость энергии системы |
от степени сжатия. |
(железо) от степени сжатия. |
Т _ ф _ П — приближение Т —Ф с квантовы
м и п о п р а в к а м и .
щие в кривые, полученные на основе рассмотренного прибли жения, по-видимому, с точностью до 10—15% описывают сжи маемость железа при давлениях, еще не достигнутых в экспе рименте, если картина не осложняется при этом фазовым пере ходом.
К настоящему времени в лабораторных условиях достигну ты давления порядка 107 атм. В природе сильно сжатое веще ство встречается в недрах многих небесных тел. Особенно больших давлений (порядка 1017—1020 атм) можно ожидать в звездах (так называемых белых карликах).
В области нормальных плотностей, где нельзя пренебре гать оболочечной структурой атома, статистическая модель должна приводить к неточным результатам. Однако в отличие от простейшей модели Т—Ф рассмотренная модель, хотя и ка чественно, но описывает некоторые свойства твердого вещест ва. Поэтому интересно рассмотреть экстраполяцию в область нормальных плотностей. Поправки к вычисленным термодина мическим величинам отрицательны и при У -*-оо стремятся к
нулю медленнее, чем неисправленные величины. Поэтому ис правленное давление обращается в нуль при некоторой конеч
ной плотности «о(2 ). |
минимум |
при |
Поскольку P=(dE/dV)T=о, энергия Е имеет |
плотности п0(Z), значение которого— энергия |
связи твердого |
тела в рассматриваемой модели. На рис. 58 представлены |
экс- |
Рис. 58. Теоретическая и экспериментальная зависимости плотности несжатого вещества от Z.
периментальная и теоретическая зависимости плотности не сжатого вещества от Z. Статистическая модель не может, ко нечно, объяснить ход экспериментальной кривой, однако и здесь видно, что квантовые поправки приводят к лучшему опи санию усредненной зависимости плотности от заряда ядра.
Обсудим вкратце физический смысл статистической теории и пределы ее применимости. В наиболее простом виде стати стическая теория описывается приближением Т—Ф, которое является квазиклассическим. Условие квазиклассичности дви жения частицы требует, чтобы длина волны была малой по сравнению с характерными неоднородностями задачи:
1 = |
« |
1. |
(41.10) |
Поскольку |
dx |
|
|
|
|
|
гд е |
|
|
|
PF = [2m(eF — u)y/‘ , |
|
то граничный импульс pF и эффективный потенциал |
должны |
быть слабомепяющимися функциями |
координат. Разумеется, |
квазиклассика будет хорошо описывать свойства |
системы, |
когда условие (41.10) выполняется если не во всей |
области, |
где система определена, то |
хотя бы |
в преобладающей ее ча |
сти. В этом случае можно ожидать, что интегральные харак теристики системы будут описываться статистической теорией достаточно корректно.
Если параметр g достаточно мал, то в нулевом приближе нии можно вообще пренебречь всеми градиентами граничного импульса и потенциала, что и приводит к теории Т—Ф. Однако при описании взаимодействия сложных атомов или при изу чении детального хода плотности электронов вблизи ядра по лучаются грубые результаты, если пользоваться статистиче
ской теорией атома в ее наипростейшем виде. |
Вблизи |
ядра |
существенна |
неоднородность системы, а вблизи |
границы |
ато |
ма или иона |
существенны как неоднородность, |
так и то, что |
плотность электронов становится малой и применимость само го статистического описания можно поставить под сомнение.
Оцепим относительную роль полученных выше обменных и квантовых поправок. Вклад этих поправок определяется пара
метром | 2, поскольку |
при разложении физических |
величин |
в |
ряд по | |
члены нечетного порядка по £ исчезают. |
Обменные |
и |
квантовые |
поправки |
имеют одинаковый порядок |
величины |
в |
том важном случае, когда область, непосредственно прилегаю
щая к точечному источнику внешнего поля (ядру), |
не играет |
определяющей роли. |
Сравнивая |
порядки |
величин |
отдельных |
членов обобщенного уравнения Т—Ф, получаем |
|
|
|
Др2 ~ р2/х2 — р*1а0, |
|
|
|
откуда |
|
|
|
|
(41.11) |
|
( * о Р р ) ~ 2 ~ |
( a o P F ) ~ 2 . |
|
где х0— характерная |
длина, определяющая расстояние, |
на |
котором не происходит заметного |
изменения граничного |
им |
пульса pF. Поскольку |
(х0P f ) ~ 2 есть |
не что |
иное, как |
| 2, то |
по |
порядку величины квантовая и обменная поправки совпадают. Из выражения (41.11) следует, что с увеличением плотности системы вклад квантовых и обменных поправок уменьшается. В самом деле, поскольку плотность
п(х) = р3(х)1 Зл2,
то
М -1
( а 0п и )
Если область вблизи ядра выделена по каким-либо причи нам, то приведенные оценки теряют смысл. В этой области
и (г)------ |
Ze2/r; р% (г) ~ 2Z/a0r, |
где г — расстояние от ядра. Отсюда
д(\/рР) 12 |
Qq |
(аоРр) 1 ( a 0Z ) |
дг |
тг |
Следовательно, при ro^a0/Z квантовые поправки велики, а об менные малы.
Модель Т—Ф предполагает независимое заполнение ячеек в фазовом пространстве. При этом в значительной мере игно рируется кваптовомеханический принцип неопределенности. Формальным источником квантовых поправок является, таким образом, некоммутативность операторов координаты и им пульса.
Статистическая теория существенно лучше описывает свой ства сильно сжатого вещества в отличие от несжатого веще ства. Действительно, вещество, находящееся под действием невысоких давлений, обладает большим разнообразием свойств. Величины, характеризующие эти свойства, обнаружи вают сильную и немонотонную зависимость от химического со става вещества. Последний же, в свою очередь, определяется свойствами периферийной области атомов, строением элек тронных оболочек и т. д. При сжатии вещества проявляется тенденция сглаживания его свойств, а при достаточно больших давлениях эта зависимость становится плавной и монотонной. При Р~^>ео наружные электронные оболочки атомов перестают существовать, поскольку входящие в их состав электроны от рываются от атомов и коллективизируются. Это обстоятельст во существенно облегчает теоретическое рассмотрение стати стической модели.
Внутреннее давление ео нетрудно оценить из следующих соображений. Обозначая энергию электрона е, число частиц N и характерную длину L, можно положить eq~ N e/L3. Пусть
для простоты |
вещество построено из |
однотипных атомов с |
зарядом ядра |
Z > 1. При этом следует |
различать три |
харак |
терных масштаба ео. В периферийной |
области атома |
е ~ е 2/ао, |
JV~ 1, L ~ a 0, поэтому |
|
|
ео1)