CC BY
19
УДК 533.9
ДЕБАЕВСКАЯ МОЛЕКУЛА В ПЛАЗМЕ
В. А. Гундиенков, С. И. Яковленко
С использованием координат Кассини численно исследовано решение уравнения Пуассона-Больцмана для облака зарядов, окружающих две заряженных пылинки, рассматриваемые как дебаевские атомы, образующие деба-евскую молекулу. Сила воздействия электрического поля на пылинку со стороны другой пылинки определялась интегрированием электростатического давления на поверхность пылинки. Показано, что при малых расстояниях между пылинками имеет место отталкивание. Силы притяжения возникают при выполнении следующих двух условий. Во-первых, дебаевский радиус (соответствующий плотности электронов на половине среднего расстояния между пылинками) должен быть примерно равен половине среднего расстояния между пылинками. При этом притяжение пылинок возникает на расстоянии, примерно равном среднему расстоянию между пылинками. Во-вторых, притяжение имеет место, если заряды одного из знаков сосредоточены преимущественно на пылинках. Если пылинки несут малую долю заряда одного из знаков, на всех расстояниях имеет место отталкивание пылинок.
Эта работа является продолжением нашего рассмотрения [1] сил взаимодействия заряженных пылинок в облаках термодинамически равновесных зарядов. В [1] основное внимание уделено термоэмиссионной плазме, когда положительно заряженные пылинки окружены электронными облаками. Здесь мы больше внимания уделяем ситуации, когда в облаках, окружающих пылинки, имеются заряды обоих знаков.
Уравнение Пуассона-Болъцмана. Как и в работах [1 - 9] мы будем исходить из уравнения Пуассона-Болыдмана в безразмерном виде:
Аv? = ехр(</?) — <5ехр( —(р), VE = —(ехр(ур) — <$ехр(—(¿>)), Е = — Vy?. (1)
Здесь длина измеряется в единицах дебаевского радиуса гд = (T/47re2iVeo)1/2, соответствующего плотности электронов Neо в тех точках, где потенциал равен нулю; Т температура газа. Связь безразмерных величин (потенциала ip, напряженности поля Е и плотности электронов пе) с соответствующими размерными величинами (ф, — Ne) определяется соотношениями: <р = фе/Т] Е = — S/ф ■ erp/T] ne = r^,Ne = п^ехр(у), где Ti£) = rzDNeQ. В уравнении (1) 8 = Nl0/Nc0 - параметр, характеризующий дополнительную ионизацию газа; iV,o - плотность положительно заряженных ионов в тех точках, где потенциал равен нулю. В силу квазинейтральности плазмы: 0 < 8 < 1.
Ниже при оценках мы, как правило, будем ориентироваться на условия экспериментов [10], в которых Ne0 — 2.5 ■ Ю10 см~3, Т = 0.146 эВ = 1700 для характерных величин имеем: гд = 18 мкм, Т/е = 0.146 ß, Т/егр = 80 В/см. При среднем радиусе пылинки гр = 0.4 мкм (г0 = Гр/гр = 2.23 • Ю-2) и ее заряде Zpe = 500е имеем напряженность поля на поверхности частицы Zpe/rl = 4.5 • 104 В/см (Eq — E(r0) = 550).
Граничные условия. Одиночную заряженную пылинку, окруженную облаком более легких зарядов, находящихся в термодинамическом равновесии, следуя [2], будем называть дебаевским атомом, а две или несколько пылинок - дебаевской молекулой [3, 4]. Рассмотрение дебаевского атома и дебаевской молекулы с формальной точки зрения отличается лишь геометрией задачи.
В реальной физической задаче задан заряд пылинок Zpe и их радиус гр (о формировании заряда пылинки см. [4, 6, 7]). Следовательно, одним из граничных условий является условие на напряженность поля на поверхности пылинок S: Е0 = — V</>|,s. При этом заряд частицы определяется выражением:
5 5
Здесь Zp - безразмерный заряд частицы, он связан с зарядом частицы в единицах электронного заряда Zp выражением Zv — iw • zp ■ ud', площадь поверхности измеряется в квадратах дебаевского радиуса.
Вторым граничным условием должно быть задание поверхности S', на которой поле равно нулю:
V^ls' = 0. (3)
Нулевое значение электрического поля на границе следует из квазинейтральности рассматриваемой системы зарядов. Поверхность S' задает границу рассматриваемой деба-евской системы.
Основная цель рассмотрения дебаевской молекулы состоит в том, чтобы найти зависимость результирующей электростатической силы, воздействующей на пылинки, от расстояния d между ними. В этом случае удобнее исходить из других граничных условий. На поверхности пылинок задается не поле, а постоянный потенциал
V?|s = <¿>0 = const. (4)
Из решения уравнения Пуассона-Больцмана находится напряженность поля Е0 на поверхности пылинки. Результирующая сила, определяется с помощью интеграла от электростатического давления по поверхности пылинки. Для того, чтобы получить нужное значение заряда zp (2), надо соответствующим образом изменить значение потенциала ip0.
В рассматриваемом случае сила взаимодействия пылинок направлена вдоль оси z и определяется выражением:
F = 11f = I E°dS*■ (5)
S S
Здесь dsz - проекция элемента поверхности ds на ось z\ сила F связана с безразмерной силой / выражением F = (Т2/87ге2) • /.
Энергия взаимодействия пылинок вычислялась по формуле
оо
U(d) = J f(x)dx + const. (6)
d
Радиус дебаевского атома определяется как расстояние г = а0, на котором заряд пылинки полностью компенсируется свободными зарядами плазмы (Е(а0) = 0). Считая, что каждая пылинка компенсируется "своими" свободными зарядами, величину ао мы выберем равной половине среднего расстояния между пылинками Go = öp/гд = (-/V~1/,3/2r£>), где Np - плотность пылинок (рис. 1). Величина ар = Лг~1^3/2 на 24% меньше радиуса Вигнера-Зейца rwz = (47r7Vp/3)-1/3.
Рассмотрим наиболее интересную ситуацию, когда радиус пылинки гр много меньше расстояния между пылинками г0 = гр/гд << а0 (подробнее см. [9]). В экспериментах [10] Np = 5 • 107 см~3, соответственно, ар = 0.755, при этом условие малости радиуса пылинки выполняется: ao/rQ = 34.
Зарядовая оболочка дебаевского атома
Пылинка Кор дебаевского атома
Дебаевская молекула
d ~ 2а<)
Рис. 1. Схема совокупности дебаевских атомов и дебаевская молекула.
При 8 —» 1 радиус дебаевского атома стремится к бесконечности: а0 оо. Это соответствует уединенной заряженной частице в квазинейтральной плазме. Дело в том, что конечный заряд частицы г0 может полностью компенсироваться квазинейтральной плазмой только при ее бесконечных размерах. Если 8 < 1, радиус дебаевского атома конечен.
О структуре дебаевского атома. Результаты рассмотрения уравнения Пуассона Больцмана для сферически симметричного случал показывают [9], что при большом заряде пылинки дебаевский атом имеет некоторый кор из зарядовой оболочки вблизи поверхности пылинки. Заряд пылинки вместе с кором равен г^ = а'о/З. Экранировка этого "оставшегося" заряда имеет место на большом расстоянии г, близком к а0. В связи с тем, что дебаевский атом имеет кор, экранирующий заряд пылинки, рассматривая взаимодействие дебаевских атомов, нельзя придавать пылинке неэкранированное значение заряда.
Условие большого заряда частицы гр = Zpe2¡г^Т > гсот может быть переписано для заряда пылинки, измеренного в единицах электронного заряда: Zv > Zcoт = (тгЛ^ео)/(6А^). Согласно измерениям [10] заряд пылинок велик: Zv = 500 > Zcoт = 262,
zp = 0.273 > zcor = 0.143. Однако из расчетов следует, что при измеренных в [10] значениях плотности и температуры электронов для данного радиуса пылинок их заряд в состоянии теплового равновесия должен иметь значение Zv = 286 (zp — 0.156), меньшее измеренного Zp « 500. Следовательно, либо измерения параметров плазмы существенно неточны, либо заряд пылинок в условиях экспериментов [10] неравновесен (см. также [6]).
Зарядовый состав дебаевского атома. Электронный и ионный безразмерный заряды, содержащиеся в зарядовой оболочке, определяются выражениями:
= /ехрМ-ИЛ, = i]* - - (.0)
г0 го
Величина 8\ = z0,/zoe дает отношение свободного заряда ионов в дебаевском атоме к заряду электронов. Вообще говоря, величина ¿i должна быть сложной функцией параметров 8, а0 и (fo- Однако, в тех случаях, когда основной вклад в интегрирование (10) вносит область малых значений потенциала <,£,(г) << можно приближенно положить 6i & 6.
Зависимости zoe, zqи Si от 8 иллюстрирует рис. 2. В результатах, представленных на рис. 2, величина а0 для разных значений 6 выбиралась максимально большой для радиуса пылинки, соответствующего экспериментам [10]: r0 = rp/rD = 2.23 • Ю-2. Это осуществлялось пристрелкой: при выборе значения ао больше того, которое представлено на рис. 2, заряд частицы становится бесконечно большим: z(r0) —► оо. Полученные зависимости z(r), ср(г) использовались для определения z0 = z(r0), <р0 = <^(го) ПРИ г0 = 0.1.
Видно, что с ростом 6 за счет увеличения объема дебаевского атома растет число как положительных, так и отрицательных зарядов в его оболочке. В то же время, число некомпенсированных зарядов zq = z0e — z0, практически не меняется с изменением 8. В рассмотренной области параметров ¿i ~ 8.
Результаты расчетов. Для точного нахождения величины силы, действующей на частицу, метод решения уравнения Пуассона-Больцмана должен обеспечивать максимальную точность в области вблизи поверхности пылинок. При этом основной интерес представляют расстояния между пылинками намного превышающие их диаметр. В обычных системах координат этого трудно добиться. Мы использовали координаты, построенные на основе известного овала Кассини [11, 12] для некоторого его частного случая (подробнее см. [1]).
ZÖe'ZQi'ZO 51>а0
Рис. 2. Связь параметров, характеризующих пылевую плазму с параметрами дебаевского атома: а) зависимость безразмерного заряда электронной, ионной оболочки z0e (сплошная кривая), z0i (пунктирная кривая) от параметра 6; штриховая прямая - их разность, б) зависимость параметра 6характеризующего отношение полного числа ионов к полному числу электронов в оболочке дебаевского атома (сплошная кривая) и радиуса дебаевского атома а0 (пунктирная кривая) от параметра 6, характеризующего отношение числа ионов к числу электронов на границе дебаевского атома. Параметры пылинки г0 = 0.1, <р(г0) = 2.4, z(r) = 0.28.
Расчеты проводились для таких параметров ip0, г0, а0, которые при d » а0 соответствуют уединенному дебаевскому атому. Для этого сначала решалась сферически симметричная задача, в которой значения поля и потенциала в точке г = а0 полагались равными нулю. Из решения этой задачи определялся потенциал <ро на частице заданного радиуса г0. Затем с этими значениями Пъ ао решалась двуцентровая задача для d = Юао. Результаты решения для сферически симметричной и для двуцентровой задачи совпадали с высокой точностью. В дальнейшей серии расчетов переходили к меньшим значениям d.
Для определения зависимости силы взаимодействия пылинок от расстояния между ними d проведены серии расчетов с заданными значениями у>0, го, При этом заряд частицы z0 также оказывается функцией d. Поэтому проводились дополнительные расчеты со значениями ip0 или а0, измененными таким образом, чтобы заряд частицы zq не зависел от d.
Серия расчетов, ориентированная на параметры плазмы работы [10], представлена
ß [13-
Расчеты показали, что при малых расстояниях между частицами d ~ г0 имеет место отталкивание. Это не согласуется с результатами численных расчетов [4, 5], в которых при d ~ г о, имело место притяжение. По-видимому, в расчетах [4, 5] была высокой погрешность вычисления производной потенциала вблизи поверхности пылинки [1].
В расчетах, ориентированных на параметры плазмы работы [10] (а0 — 7.55), смена притяжения отталкиванием происходит на расстоянии do ~ 1.3, что несколько меньше среднего расстояния между пылинками 2а0 = 1.51. Положение точки г = d0, в которой происходит смена знака силы, слабо зависит от того, какие величины ао или z0, а0) сохранялись в расчетах при изменении d. Изменение а0 (при постоянных значениях Zo, tfo) влияет на значение do более существенно. По-видимому, удержание постоянного заряда zo = const за счет изменения потенциала частицы = <fo(d) больше соответствует физике взаимодействия заряженных пылинок.
Зная силу притяжения пылинок F(2a0) на среднем межчастичном расстоянии 2ао, можно оценить электростатическое давление сжимающее газ пылинок
РЕ « F(2a0) ■ N2J3 = (Г2/8тге2) • N2p'z ■ f(2a0) (11)
и поверхностное натяжение "пылевой жидкости" сге ~ F • N^3 = (N^3T2/Sire2) ■ /(2а0).
Сравнивая электростатическое давление на пылинки с газокинетическим давлением пылинок и газокинетическим давлением свободных электронов имеем: Pe/NpT =
(Г/Втге2^/3) - /(2ао).
В условиях экспериментов [10]: |/(2а0)| « 0.2; РЕ = 9.7 • Ю-7 • |/(2а0)| • Topp к 2 • 10~7 • Topp; (тЕ = 3.5 • Ю-9 ■ |/(2а0)| • Н/м « 7 • Ю"10 • Н/м; PE/NpT к 20.
Из приведенных оценок следует, что газ дебаевских атомов в смеси с инертным газом должен проявлять в условиях экспериментов [10] тенденцию к сжатию. Такого рода ситуация рассмотрена в [13 - 15]. Рассмотрение вопросов влияния взаимодействия дебаевских атомов на газодинамические свойства пылевой плазмы выходит за рамки данной работы.
Было проведено несколько серий расчетов для различных значений ао (рис. 3). Расчеты показали, что притяжение имеет место лишь при а0 < 1. Уже при ао > 1.12 точка смены знака удаляется на большое расстояние do > 4ао-
Условие а0 = ар/г£) < 1 можно переписать для размерных величин в виде:
Ne0 > Necr = —Ж'2. (12)
TT f>L г
do(ao), 4(а0)
Рис. 3. Зависимость координаты с10 точки смены знака силы (жирные кривые) и крутизны силы £ в точке (10 (тонкие кривые с квадратиками) от размера дебаевского атома а0. Сплошные кривые соответствуют не зависящему от (1 потенциалу на поверхностях пылинок <р0; пунктирные кривые соответствуют постоянному заряду пылинки, обеспеченному подбором (р0{(1). Потенциал ¡р0 определялся при г0 = 1 для заданного значения а0.
Электростатические силы сжатия обращаются в нуль, когда do = 2ао, т.е. при ао = 1. Соответственно условие ао = 1 или Ne = Necr является условием равновесия сил отталкивания и притяжения в газе дебаевских атомов.
При этом условие большого заряда частицы zp > 1/3 может быть переписано для заряда пылинки в единицах электронного заряда следующим образом:
7Г NP,r
zp > ZecT = -
В условиях эксперимента имеем Necт = 4.4 • Ю10 см~3, £есг = 460. Эти величины порядка измеренных в экспериментах [10]: Л^о ~ 2.5 • Ю10 сл«-3, Zv ~ 500.
Ввиду того, что на основании рассмотрения двух дебаевских атомов нельзя определить глубину потенциальной ямы без учета воздействия других частиц, силу взаимодействия дебаевских атомов будем характеризовать крутизной в точке пересечения с осью абсцисс: £ = /'(¿)|¿=в.0 = ~и"
Через величину £ выражается частота колебаний пылинок вокруг положения равновесия: ш = • где шо — ит/ар, Vт — (2Т/тпРУ^2 - тепловая скорость пылинок, тр - их масса. В условиях экспериментов [10] имеем тпр ~ 2 • Ю-12 г, 1>т ~ 0.5 см/с, ар ~ 1.4 • 10_3 см. Отсюда следует и>о = 357с-1 частота колебаний, 27г/о;о = 18 мс - период колебаний. Из рис. 3 видно, что наиболее сильная связь имеет место при 0.5 < а0 < 1. В этих условиях газ дебаевских атомов должен стремиться к сжатию (ср. [15]).
Для пылинок малого радиуса меньше и размер дебаевского атома. Например, максимальное значение радиуса дебаевского атома а0тах = а0(г0 —> оо), как функция г0 может быть аппроксимирована выражением: а0тах = З-Гд-3, или г0 = (аотах/3)10/3 при г0 < 0.02 [1]. Поэтому радиус пылинки должен быть не слишком большим и не слишком малень ким. При 0.5 < аотах < 1 имеем условие 2.5 • 10~3 < г0 < 2.6 • Ю-2. В экспериментах [10] Го = 2.23 • Ю-2, и это условие выполняется.
т
Рис. 4. Иллюстрация зависимости положения координаты ¿0 точки смены силы от радиуса пылинки г0: а) зависимость силы от расстояния между частицами (I вблизи точки «¿0 для различных значений г0 (сплошная кривая - г0 = 0.2, пунктирная кривая - г0 — 0.25, сплошная кривая - г0 = 0.3, штрих-пунктирная кривая - г0 = 0.4, сплошная тонкая кривая г0 = 0.5^; б) зависимость координаты (10 точки смены знака силы от г0/а0. Здесь г0 можно рассматривать, как радиус области, в которой пренебрегаете я поляризацией зарядового облака. Потенциал (¿>0 для г — г0 определялся при а0 = 0.755.
В расчетах небольшой заряженный проводящий шарик заменяется проводящим шариком большего размера с зарядом, частично компенсированным свободными зарядами оболочки дебаевского атома. Возникает естественный вопрос, насколько такая замена
адекватна. Для этого было проведено несколько серий расчетов с разными значениями Го и, соответственно, </?о- При радиусах пылинки, малых по сравнению с радиусом дебаевского атома ао, различие результатов расчетов невелико. Например, в случае ао = 0.755 (рис. 4) при изменении радиуса пылинки в диапазоне г0 = 0.1 — 0.2 (и выборе значения <р0, соответствующего данному значению г0) отличия в положении точки смены притяжения отталкиванием ¿0 = 1.28 имеют разброс менее 2%, что соответствует имеющейся точности расчета.
Существенное влияние размеров пылинки имеет место при г0 > 0.3ао. При г0 > 0.4 поляризационное притяжение настолько уменьшается, что точка смены знака силы становится больше среднего межчастичного расстояния <10 > 2ао- Из этого можно заключить, что существенный вклад в поляризационные силы вносят заряды, находящиеся на расстоянии т % О.Зао от центра пылинки, а не только периферия дебаевского атома г яа а0. Иначе говоря, сила притяжения формируется за счет поляризации большинства электронов зарядовой оболочки. В связи с этим сложно надеяться на вычисление сил притяжения с помощью каких-либо приближенных методов.
В случае двух дебаевских атомов, находящихся в почти квазинейтральной плазме (1 — ¿<<1), расстояние между ними можно брать сколь угодно большим. Однако для того, чтобы взаимодействие можно было считать бинарным, необходимо выполнить условие » <1 • гд.
Как и в случае 8 = 0, для определения зависимости силы взаимодействия пылинок от расстояния (1 проводились серии расчетов с заданными значениями <р0, г0, а0 для дебаевского атома. При этом проводились дополнительные расчеты со значениями </?о или а0, измененными таким образом, чтобы заряд частиц г0 не зависел от ¿. Как и в случае 8 = 0 мы выбирали Го больше радиуса кора атома, моделируя пылинку проводящим шариком большего размера с зарядом, частично компенсированным свободными зарядами оболочки дебаевского атома. Тем самым поляризацией кора пренебрегалось.
В результатах представленных на рис. 5, величина а0 для разных значений 8 (см. табл.) выбиралась соответствующей предельно большому заряду для радиуса пылинки, соответствующего экспериментам [10]: Го = 2.23-Ю-2. Это осуществлялось пристрелкой (ср. выше). Полученные зависимости г (г), <р(г) использовались для определения го = г(г0), (р0 = <р(го), при г0 = 0.1.
При 1 — 6 << 1 в рассмотренном диапазоне параметров притяжение пылинок на расстоянии, меньшем среднего межчастичного состояния, обнаружить не удалось. Притяжение при с1 < 2ао возникает лишь в тех случаях, когда заметную долю положитель-
ного заряда плазмы несут пылинки (при 8 < 0.5, см. рис. 5). При этом, чем меньше доля зарядов плазмы, тем больше максимальная сила притяжения и глубина потенциальной ямы.
и(ё)
0.1 1 10 100 й
Щс!)
10 100 1 с1
10 100 с!
Рис. 5. Зависимость проекции силы на ось х (а, в) и потенциальной энергии (б, г) взаимодействия пылинок от расстояния между ними в, для разных значений 6 ф 0. Нормировка потенциальной энергии выбрана так, чтобы выполнялось условие: 11(75) = 0 для (б), 11(25) — 0 для (г). Во всех случаях г0 = 0.1. На рис. а, б сплошные кривые соответствуют аналитическим выражениям работы [17] при с = 12. Пунктирные кривые на рис. а, б соответствуют 6 - 0.999, штриховые кривые соответствует 6 = 0.9. На рис. в, г сплошные кривые соответствуют 6 = 0.1; пунктирные кривые - 8 — 0.3; штриховые кривые - 6 = 0.5; штрих-пунктирные кривые - 6 = 0.7.
Ослабление сил притяжения с ростом 8 имеет простое объяснение. Силы притяжения возникают за счет того, что электроны скапливаются вблизи оси х между центрами пылинок, и обеспечивают притяжение к центру дебаевской молекулы. Это притяжение превышает силу отталкивания зарядов пылинок, экранированную внутренними слоями
электронных оболочек дебаевских атомов. При 1 — 8 << 1 эффект экранировки заряда пылинки внутренними слоями электронных оболочек остается прежним. Однако притяжение к центру молекулы существенно ослабевает из-за того, что в этой области концентрируются не только электроны, но и положительные заряды, отталкивающие пылинки.
Таблица Параметры расчетов для различных значений 8
8 а0 ¥>о zo
0.999 4.1 2.455 0.282
0.9 1.71 2.426 0.283
0.7 1.288 2.413 0.272
0.5 1.12 2.378 0.286
0.3 1.018 2.364 0.286
0.1 0.94 2.292 0.279
В случае малых значений заряда плазмы 8 << 1 глубина потенциальной ямы довольно велика, порядка нескольких температур. Однако следует помнить, что бинарное рассмотрение ограничено величиной порядка диаметра дебаевского атома 2ао- При с1 » 2ао становится существенным отталкивание другими дебаевскими атомами (см. рис. 1). Поэтому глубину потенциальной ямы оценить нельзя. Расчеты, проведенные для больших расстояний >> 2а0, могут лишь качественно характеризовать энергию удаления дебаевского атома с поверхности дебаевской жидкости.
Изложенный выше вывод об отсутствии притяжения при 1 — 8 << 1 не согласуется с данными недавних приближенных аналитических рассмотрений [16, 17, рис. 11]. Из результатов этих работ следует, что при 8 — 1 в области г > (31//2 + = 1.93 при
рассмотрении линеаризованного уравнения Пуассона-Больцмана имеет место притяжение пылинок.
Этот результат вызывает удивление. Дело в том, что в линейном приближении [16, 17] воздействие точечных зарядов (находящихся в точках гг и г2) на плазму считается независимым, так что значение в точке г определяется суммой экранированных потенциалов точечных зарядов:
¥>(г) = Ф(|г - п|) + Ф(|г - г2|), Ф(х) = (г^/х) ■ ех^-я/г1'2), |г, - гг| = д.
Согласно простым соображениям, изложенным в [1], в отсутствие перестройки зарядовой оболочки одной пылинки под воздействием другой пылинки для одноименно заряженных пылинок не могут иметь место силы притяжения. В линейном приближении [16, 17] должно быть:
№ = -гр ■ дФ(х)дх\х=а = г2р ■ ((21'2 + <1)/<12) ■ ехрМ/21/2) > 0, (13)
что соответствует отталкиванию.
Неточность результатов [16, 17], по-видимому, связана со следующим обстоятельством. К силе (13), действующей непосредственно на пылинку, в работах [16, 17] прибавлена сила притяжения, действующая на электронную оболочку пылинки со стороны второй пылинки. Такое приближение было бы оправданным, если бы зарядовые оболочки пылинок были бы жестко связаны с пылинками какими-то другими силами. Однако таких сторонних жестких сил в рассматриваемой задаче нет. Наличие вычисленной в работах [16, 17] силы притяжения электронной оболочки одного заряда к другому заряду говорит лишь о том, что данная конфигурация зарядовой оболочки не является равновесной, и сила притяжения к другому заряду должна приводить к поляризации зарядовой оболочки, которая в работах [16, 17] не учтена. Нет никаких оснований складывать эту поляризующую силу с силой, действующей непосредственно на пылинку.
Здесь уместна аналогия с поляризационными силами притяжения атомов, существенными на больших расстояниях между ядрами по сравнению с размерами электронных оболочек. Как известно, для сферически симметричных атомов поляризационное взаимодействие в первом порядке теории возмущений не имеет места. Оно возникает лишь во втором порядке теории возмущений, когда учитывается поляризация электронной оболочки одного атома зарядами другого атома. Обычный атом отличается от дебаевского атома лишь тем, что в нем электроны оболочки двигаются по квантовым, а не классическим законам. Природа же возникновения поляризационных сил у обычного атома и дебаевского атома одинакова.
Подведем итоги проведенного рассмотрения.
1. Дебаевский атом имеет структуру: кор и электронную оболочку. При больших зарядах пылинки область электронного облака высокой плотности (кор) существенно экранирует большой заряд пылинки вблизи ее поверхности. В связи с этим, рассматривая взаимодействие дебаевских атомов, нельзя приписывать пылинке неэкранированное значение заряда. Заряд пылинки, экранированной кором, имеет универсальное значение, определяемое расстоянием между пылинками. Он экранируется электронной оболочкой дебаевского атома.
2. Силы притяжения обусловлены поляризацией зарядовых оболочек дебаевских атомов. В отсутствие поляризации отсутствует и притяжение. Сила притяжения формируется за счет поляризации большей части электронов зарядовой оболочки. В связи с этим сложно надеяться на вычисление сил притяжения с помощью каких-либо приближенных методов. Поляризация кора несущественна.
3. Силы притяжения пылинок возникают на сравнительно большом расстоянии, примерно равном среднему расстоянию между пылинками. При этом дебаевский радиус должен быть примерно равен половине среднего расстояния между пылинками.
4. Притяжение имеет место, если заряды одного из знаков сосредоточены преимущественно на пылинках. Если пылинки несут малую долю заряда одного из знаков, на всех расстояниях имеет место отталкивание пылинок.
5. При некотором соотношении между плотностью электронов и плотностью пылинок имеет место равновесие дебаевской "жидкости": электростатические силы взаимодействия пылинок обращаются в нуль.
Ввиду того, что силы притяжения возникают на расстояниях порядка среднего расстояния между пылинками, корректно решить задачу о формировании пылевых жидкостей и кристаллов можно лишь, учитывал многочастичное взаимодействие пылинок. На самом деле поляризационные силы, действующие на пылинку, флуктуируют и зависят от конфигурации пылинок, ее окружающих. Однако из проведенного рассмотрения ясно, что равновесная конфигурация пылинок должна иметь место.
Авторы признательны А. Н. Ткачеву за обсуждение результатов данной работы и работ [12, 13], а также Ю. И. Сыцько за обсуждения вычислительных аспектов задачи.
ЛИТЕРАТУРА
[1] Гундиенков В. А., Яковленко С. И. Краткие сообщения по физике ФИАН, N 12, 3 (2001).
[2] Т к а ч е в А. Н., Яковленко С. И. ЖТФ, 69, N 1, 53 (1999).
[3] Т к а ч е в А. Н., Яковленко С. И. Письма в ЖТФ, 25, N 1, 52 (1999).
[4] Я к о в л е н к о с. И. Письма в ЖТФ, 25, N 16, 83 (1999).
[5] Я к о в л е н к о с. И. Краткие сообщения по физике ФИАН, N 9, 3 (1999)
[6] Я к о в л е н к о с. И. Письма в ЖТФ, 26, N 16, 47 (1999).
[7] я к о в л е н к о с. И. Письма в ЖТФ, 26, N 26, 38 (2000).
[8] я к о в л е н к о С. И. Письма в ЖТФ, 27, N 9, 83 (2001).
[9] Я к о в л е н к о С. И. Краткие сообщения по физике ФИАН, N 1,3 (2002).
10] Фортов В. Е., Нефедов А. П., Петров О. Ф. и др. ЖЭТФ, 111, N 2, 467 (1997).
11] Математический энциклопедический словарь. Гл. ред. Ю. В. Прохорова. М., Научн. Изд. "Большая российская энциклопедия", 1995.
12] Б р о н ш т е й н И. Н., Семендяев К. А. Справочник по математике. М., Наука, 1964.
13] М а й о р о в С. А., Ткачев А. Н., Яковленко С. И. Усп. физич. наук, 164, N 3, 298 (1994).
14] М а у о г о V S. A., Tkachev A. N., Yakovlenko S. I. Physica Scripta, 51, 498 (1995).
15] Я к о в л е н к о С. И. Изв. ВУЗов, Физика. 38, N 4, 3 (1995) [in Russian, for English translation see: Russian Physics J., 38(4), 329 (1995)].
16] Герасимов Д. H., Синкевич О. А. Теплофизика высоких температур, 37, N 6, 853 (1999).
17] I V а п о V A. S. Phys. Letters, А290, 304 (2001).
Поступила в редакцию 14 марта 2002 г.
