14
.doc
21.1Дебаевское экранирование Рассмотрим потенциал неподвижного точечного заряда в изотропной плазме и запишем Фурье-образы уравнений , (1.43) и E(r) = – grad (r). (4.4) : ikD(k, ) = 4e(k, ), (6.7) E(k, ) = –ik(k, ). (6.8) Уравнение (6.7) с учетом соотношения (3.32) можно преобразовать к виду , который совместно с уравнением (6.8) дает . (6.9) Если поле создано одним точечным зарядом, находящимся в точке r0, то e(r, t) = e(r0), соответственно . Подставляя это соотношение в уравнение (6.9), получим . Проведем обратное преобразование Фурье вида (3.27), вернувшись к переменным r и t: |
21.2равной плотности вероятности обнаружить в момент времени t частицу с импульсом p в окрестности точки r. Функция распределения замкнутой системы удовлетворяет уравнению Лиувилля df/dt = 0. Если электроны испытывают рассеяние, например, при столкновениях с ионами или нейтральными атомами, уравнение Лиувилля принимает вид кинетического уравнения Больцмана (6.11) Здесь = 1/ – время релаксации системы, – частота столкновений, f0 – равновесная, невозмущенная электромагнитным полем, функция распределения системы, к которой она релаксирует с постоянной времени при выключении поля. Обозначим в уравнении (6.11) f = f – f0 – вызванное электрическим полем малое возмущение функции распределения. Пусть магнитное поле отсутствует. Тогда производная импульса электрона, то есть сила, действующая на него, определяется только электрическим полем . В рамках линейной электродинамики в уравнении (6.11) следует оставить только линейные по полю слагаемые, а слагаемые, содержащие произведение Еf, квадратичные по полю, должны быть отброшены. Тогда Фурье-образ уравнения (6.11) принимает вид . Решение этого уравнения есть . (6.12) |
21.3 что . (6.16) В случае бесстолкновительной плазмы, то есть для = 0, что соответствует бесконечной проводимости по постоянному току, неопределенность, возникающую в формуле (6.16) при = 0, раскроем по правилу Лопиталя: . Таким образом, , (6.17) где обозначено – (6.18) |
|
|
22. 1Обобщенная проницаемость плазмы Формулу (6.16) для продольной диэлектрической проницаемости изотропной плазмы удобно записывать в инвариантном виде . (6.20) Поперечная проницаемость плазмы, как следует из формулы (3.32) является коэффициентом пропорциональности между электрическим полем E, перпендикулярным волновому вектору k, и обобщенной индукцией D. Запишем в инвариантном виде связь (3.29) обобщенной индукции с электрическим полем, учитывая формулу (6.15): . Преобразование подынтегрального выражения в этом соотношении можно обосновать, записав интегралы в координатной форме и направив ось х вдоль вектора Е, а ось у вдоль вектора k. Сравнивая правую часть формулы с левой, непосредственно получаем |
22.2 – (6.25) плазменная (ленгмюровская) частота электронов. Проводя аналогичное интегрирование в формуле (6.18), получим: . (6.26) В металлах дебаевский радиус экранирования порядка 10–7 см, в полупроводниках порядка 10–5 см, в термоядерной плазме 10–3 ... 10–4 см, в ионосфере 10–1 см, а в межпланетной плазме 103 см. В предельном случае длинноволновых колебаний знаменатель подынтегрального выражения в формулах (6.23) и (6.24) можно разложить по степеням малого параметра . Для бесстолкновительной плазмы при 0 получим . (6.27) Мнимая часть обобщенной продольной проницаемости плазмы (6.27) растет с ростом волнового вектора. Так как в изотропной среде именно мнимая часть обобщенной проницаемости обусловливает диссипацию энергии поля, можно сделать вывод, что с уменьшением длины волны будет увеличиваться затухание продольных колебаний поля. |
|
|
|
23.1Магнитогидродинамика плазмы Помимо кинетического рассмотрения, использованного в предыдущих двух разделах, часто используется другой, более грубый, гидродинамический метод описания. В сложных задачах гидродинамический метод часто используется для оценочного рассмотрения и основан на том, что в уравнения Максвелла функция распределения непосредственно не входит. В эти уравнения входят плотность индуцированного заряда и плотность индуцированного тока, которые выражаются через интеграл вида от функции распределения , (6.13) . Иначе говоря, в уравнения Максвелла входят только среднее значение скорости частиц среды. В уравнения же механики, определяющие характер реакции среды на поле, скорость входит непосредственно. Если детали фактического распределения частиц по скоростям мало сказываются на характере реакции среды на внешнее поле, возможно гидродинамическое описание движения среды, в которое с самого начала входит только средняя скорость. В гидродинамическом приближении свойства среды описываются пятью скалярными величинами: тремя компонентами средней скорости частиц среды и двумя термодинамическими величинами, описывающими внутреннее состояние среды в каждой точке. |
23.2. Продольная и поперечная компоненты плотности тока при этом имеют вид . (6.32) Формулы (6.28) – (6.32) написаны для одной из компонент плазмы. Поскольку ток обусловлен всеми ее компонентами, в формулах (6.32) следует подразумевать суммирование по всем компонентам. Из соотношений (6.32) непосредственно получаются выражения для продольной и поперечной проводимостей плазмы , (6.33) где производится суммирование по всем компонентам плазмы. В соответствии с формулами – (3.16) и (6.25) для продольной и поперечной обобщенных проницаемостей изотропной немагнитной плазмы получаем |
23.3 Тогда уравнение движения (6.28) принимает вид уравнения Навье – Стокса [v/t + (v ) v] = –p + [j B]/c + Fc. (6.35) Уравнение движения (6.35) помимо уравнения непрерывности (6.29) и уравнений Максвелла c rot E = –B/t, c rot B = 4j, div B = 0, (6.36) в которых для низкочастотных колебаний опущен ток смещения и положено для немагнитной плазмы В = Н, необходимо дополнить материальным уравнением j = (E + [v B]/c) (6.37) Отметим, что в силу уравнений (2.7) комбинация E + [v B]/c есть электрическое поле в системе, движущейся вместе с плазмой. Именно в такой системе определен закон Ома ji = ikEk, (1.40) . Выразим из уравнения (6.37) электрическое поле E = j/ – [v B]/c и подставим это значение в первое из уравнений Максвелла (6.36), учтя при этом второе и третье уравнения: B/t = rot [v B] + с2В/(4). (6.38) Для неподвижной среды уравнение (6.38) принимает вид уравнения диффузии: с2В – 4B/t = 0. (6.39) Из структуры уравнения (6.39) можно определить время диффузии магнитного поля на характерное расстояние L, на котором существенно изменяется магнитное поле: td = 4L2/с2. (6.40) Другому предельному случаю отвечает движущаяся среда с очень высокой проводимостью. При уравнение (6.38) переходит в B/t = rot [v B] = (B)v – (v)B – B div v + v div B Последнее слагаемое в правой части этого уравнения равно нулю в силу третьего уравнения Максвелла (6.36), а третье слагаемое может быть преобразовано с учетом уравнения непрерывности (6.29): . |
23.4 В этом случае в правой части уравнения (6.35) отсутствует сила Fс, а в правой части уравнения (6.38) отсутствует второе слагаемое. К этим уравнениям необходимо еще добавить уравнение состояния. При колебательных процессах в жидкости давление р обычно связано с плотностью соотношением p = s2, (6.44) где s – скорость звука в среде. Линеаризуем уравнения (6.35), (6.36) и (6.44) совместно с уравнением непрерывности (6.29) для малых отклонений от положения равновесия = 0 + 1, р = р0 + р1, В = В0 + В1. Полагая, что v0 = 0, получим . Эту систему можно разрешить относительно скорости 2v/t2 – s2 grad div v + [va rot rot [v va]] = 0. (6.45) Здесь – альвеновская скорость. Будем искать решение уравнения (6. 45) в виде плоской магнитогидродинамической волны v = va exp(ikr – it). Подстановка в уравнение (6.45) дает . (6.46) Если волна распространяется перпендикулярно магнитному полю, то (kva) = 0 и уравнение (6.46) упрощается: . Поскольку направления векторов k и v должны совпадать, |
|
24. 1Поверхностные волны Возьмем ротор от уравнения (1.28) и подставим в него уравнение . (1.44) с учетом материального уравнения (1.42). Для среды без сторонних зарядов и токов получим . Пусть среда является неоднородной, но изотропной и немагнитной ( = 1). Для гармонической волны в такой среде без пространственной дисперсии получаем: . (7.1) Можно также получить уравнение относительно магнитного поля Н = В: . (7.2) Важной особенностью неоднородных сред являются новые решения уравнений (7.1) и (7.2), которые не существуют в пространственно-однородных средах. Наиболее простой вид такие новые моды электромагнитных волн имеют в предельном случае резкой границы раздела двух сред. |
24.2 B1n = B2n. (1.34) D2n – D1n = 4. (1.35) HmS = HvS, EmS = EvS, DmnS = DvnS. (7.5) Решения уравнений (7.3) и (7.4), описывающие поверхностные волны, должны стремиться к нулю при z , поэтому будем их искать в виде Hm = Hm0exp(iqx – kmz), Hv = Hv0exp(iqx + kvz). (7.6) Подставляя соотношения (7.6) в уравнения (7.3) и (7.4), получим: , , . (7.7) Из равенства тангенциальных компонент магнитного поля (7.5) следует, что Hm0 = Hv0 = H0. Равенство тангенциальных компонент электрического поля означает, что . Из этого условия получается , или . (7.8) Уравнение (7.8) имеет решение лишь при . закон дисперсии: . (7.9) |
|
|
|
|
|
|
При стремлении R к нулю значение |AS|R также стремится к нулю. Таким образом, все компоненты векторного потенциала оказываются непрерывной функцией координат и в том случае, когда имеется резкая граница раздела сред. Рассмотрим систему линейных токов, сосредоточенных в ограниченной области, причем все пространство между проводниками заполнено однородным магнетиком с проницаемостью . при такой постановке задачи все линейные токи представляют собой замкнутые витки. Наряду с током Ii, важной характеристикой такого витка является магнитный поток через контур витка . (5.12) Подставляя формулу (5.8) в уравнение (5.12), видим, что потоки и токи связаны между собой линейно , (5.13) где – (5.14) взаимная индуктивность i-го и j-го проводников (при i j). Величина Lii называется индуктивностью i-го проводника. Коэффициенты Lij зависят от геометрии системы: формы, размеров и взаимного расположения проводников. По аналогии с электростатической теоремой взаимности (4.13) и (4.14) можно доказать, что матрица взаимных индуктивностей симметричная: Lij = Lji. (5.15). Вводя матрицу (симметричную) , обратную к матрице Lij, можно переписать формулу (5.13) в виде . (5.16) |
Соотношение (5.6), таким образом, позволяет найти векторный потенциал А по заданному распределению сторонних токов. Взяв ротор от выражения (5.6) с учетом соотношения (5.3) и материального уравнения (1.37), получим . (5.7) Из формулы (5.7) видно, что напряженность поля Н в пространственно-неоднородном магнетике при заданном распределении сторонних токов je такая же, какой она была бы в отсутствии всякой среды. Будем рассматривать распределение поля только в области вне проводников, достаточно тонких, так что направление тока в данной точке проводника совпадает с направлением dl – элемента длины проводника. Такую систему принято называть системой линейных токов. Тогда, выполняя в соотношениях (5.6) и (5.7) интегрирование по поперечному сечению проводников, получим , (5.8) , (5.9) где Ii – полный ток, протекающий по i-му проводнику. Выражение (5.9) представляет собой закон Био – Савара. |
|
|
|
сти свободной энергии F = W – TS равно . (5.20) Соответственно, изменение плотности потенциала G = F – BH/(4) составляет . (5.21) Кинетические явления в магнитном поле Магнитное поле оказывает влияние на поток электрических зарядов и тепла, возникающих в среде под действием температуры и электрического поля. Как показано в п. 4.6, эти потоки могут быть записаны с учетом соотношения (4.70) в виде уравнений (4.71) и (4.75) E = j + grad T, (5.22) s = j – kgrad T, (5.23) где = 1/, = –/, = /. При включении магнитного поля с напряженностью Н кинетические коэффициенты , , и k в уравнениях (5.22) и (5.23) начинают зависеть от Н. Если магнитное поле слабое, эти коэффициенты можно разложить в ряд по степеням Н. Ограничиваясь линейными членами, получим E = 0j + 0grad T + RH j + NH grad T, (5.24) s = 0j – k0grad T + NTH j + LH grad T. (5.25) Нулевой индекс у коэффициентов , , и k означает, что их значения берутся в отсутствии магнитного поля. В коэффициентах при H grad T в уравнении (5.24) и H j в уравнении (5.25) учтен принцип Онсагера вида (4.70). Слагаемое RH j в уравнении (5.24) описывает влияние магнитного поля на ток. В соответствии с этим выражением, магнитное поле приводит к появлению электрическому поля, перпендикулярного току и пропорционального напряженности магнитного поля и току. Это явление называется эффектом Холла. |
(5.17) где интегрирование выполняется только по объему проводников. Если магнитное поле создается системой линейных токов, то в соотношении (5.17) можно выполнить интегрирование по поперечному сечению проводников. С учетом соотношения (5.12) получаем . (5.18) Выражение (5.18) можно преобразовать так, чтобы энергия выражалась только через токи Ii или потоки Фi. Воспользовавшись соотношениями (5.13) и (5.15), получаем . (5.19) Рассмотрим вопрос о замкнутости системы, магнитная энергия которой определена выражением (5.18). Энергия замкнутой системы должна сохраняться. С другой стороны, проинтегрировав по всему пространству уравнение w/t + divS = –(j + je)E, (1.46) и учитывая, что на бесконечности поля, созданные ограниченной системой проводников, затухают, получим, что энергия (5.18) системы сохраняется, если обращается в нуль интеграл . В общем случае этот интеграл определяет работу электромагнитного поля над средой, то есть уменьшение энергии электромагнитного поля. При протекании постоянного тока по обычным проводникам уменьшение энергии поля связано с выделением джоулева тепла и компенсируется за счет источников: электрических батарей, аккумуляторов и т. п |
|
|
Поэтому даже в переменных полях уравнения (5.31) и (5.32) принимают вид: , (5.36) . (5.37) Уравнения (5.36) и (5.37) называются, соответственно, первым и вторым уравнениями Лондонов. Обычно в правой части уравнения (5.36) пренебрегают вторым слагаемым, считая плотности тока достаточно малыми. Уравнение (5.37), записанное в виде , (5.38) является материальным уравнением сверхпроводника. Для сверхпроводящего полупространства x > 0 его можно записать в виде 2d 2B/dx2 = B. Затухающее на бесконечности в силу эффекта Мейсснера – Оксенфельда решение этого уравнения имеет вид B(x > 0) = B0exp(–x/), то есть величина , определенная соотношением (5.33), описывает глубину проникновения постоянного магнитного поля в сверхпроводник и называется лондоновской глубиной проникновения |
. (1.8) принимает вид . Если следить не за фиксированными электронами, а полем их скоростей в фиксированной точке пространства, то, учитывая, что dv/dt = v/t + (v grad) v = v/t + grad (v2/2) – [v rot v], уравнению движения электронов можно придать вид уравнения Эйлера: . Учитывая, что j = env, где n – концентрация электронов, перепишем уравнение движения электронной жидкости в сверхпроводнике в виде , (5.31) где обозначено , (5.32) . (5.33) В соотношении (5.32) учтено, что сверхпроводник является парамагнитным металлом и можно принять В = Н. Возьмем ротор от обеих частей уравнения (5.31) и с учетом уравнения Максвелла (1.28) получим, что: enw/t = rot [j w]. (5.34) Найдем теперь полную внутреннюю энергию сверхпроводника как сумму энергии магнитного поля и кинетической энергии сверхпроводящих электронов: |
Такие оболочки (3d) имеются у элементов подгруппы железа (железо, кобальт, никель). По аналогии с соотношением (4.58), полученным для сегнетоэлектриков, свободную энергию однородного ферромагнетика во внешнем магнитном поле, опуская энергию самого магнитного поля, можно записать в виде: G(M, H, T) = G0(M2, T) – MH, (5.28) где первое слагаемое учитывает свободную энергию ферромагнетика в отсутствии поля, а второе – энергию его взаимодействия с магнитным полем. Величина G0 в формуле (5.28) зависит от квадрата модуля магнитного момента, так как появление спонтанной намагниченности связано с обменными силами, которые определяются только взаимной ориентацией соседних спинов и не зависят от их ориентации относительно кристаллографических осей. Разлагая функцию G0(M2, T) в ряд по степеням M2 и повторяя рассуждения, приведенные в п. 4.5 для сегнетоэлектриков, нетрудно получить формулы (5.27). В постоянном и однородном поле на магнетик действует сила, анализ которой аналогичен проведенному в п. 4.3 и 4.5. Заменяя в формулах . (4.45) и , (4.46) на –1 и Е на В, получаем: . (5.29) |
|
|
При таких столкновениях энергия электрона сохраняется, а импульс меняется на противоположный, то есть происходит рассеяние назад. Скорость электрона в переменном высокочастотном поле меняется по закону (6.2), а для его средней за период поля энергии получаем с учетом соотношения (6.3): , где p0 = mv0. Пусть в момент t = t0 происходит мгновенное упругое столкновение с изменением импульса на обратный, то есть p(t0 + 0) = p(t0 – 0). Тогда . Средняя за период энергия электрона после столкновения равна , а изменение энергии при столкновении составит . . . |
. (6.4) Здесь учтено, что поток энергии в бегущей волне I = cE2/(8), r0 – классический радиус электрона. В первом порядке по v/c следует учесть силу Лоренца в уравнении движения по оси z: . (6.5) Это уравнение показывает, что в первом приближении по v/c электрон колеблется вдоль оси z с удвоенной частотой, описывая фигуру Лиссажу – восьмерку в плоскости xz. Эта восьмерка сильно вытянута вдоль оси х, так как отношение амплитуд колебаний по осям z и х порядка v/c. Сечение рассеяния для второй гармоники составляет (2) = ()Ek/(mc2) << (). В следующих приближениях по v/c появляются и следующие гармоники, но вычислять их надо на основе релятивистских уравнений движения, так как они порядка v2/c2 и выше. Для круговой поляризации поля все нелинейные эффекты исчезают, так как скорость электрона все время параллельна направлению магнитного поля и сила Лоренца равна нулю. Для свободного электрона, взаимодействующего с электромагнитным полем, классическая теория, включая и формулу Томпсона (6.4), применима при условии ħ << mc2, которое выполняется вплоть до гамма-излучения. Для коротковолнового излучения с длиной волны порядка ħ/mc 10–10 см (комптоновская длина волны) рассеяние происходит с изменением частоты рассеянного света (эффект Комптона). |
, (6.2) где v0 – начальная скорость электрона. Электрон колеблется в поле, и его средняя кинетическая энергия колебательного движения равна . (6.3) Смещение электрона за период поля r, отнесенное к длине волны = с/, при этом порядка v/c, поэтому в нерелятивистских задачах поле, действующее на электрон, мало отличается от гармонического. Для того чтобы учесть члены первого порядка по v/c, разложим электрическое поле по степеням r и ограничимся линейным слагаемым: . Подставляя это выражение в уравнение (6.1) и учитывая, что в силу уравнения Максвелла (1.28) в гармонической волне rot E = –B/c, получим после усреднения по периоду поля: |