![]() | проект "СиЗиФ" |
Е.E. Антонова, М.Ф. Бахарева, В.Н. Ломоносов, Б.А. ТверскойУскорительные механизмы в космосеУчебное пособие НИИЯФ и Физ.фак. МГУ, 1988 |
Глава 1 Бетатронное ускорение и формирование радиационных поясов Локальными космофизическими объектами, содержащих большое число энергичных частиц являются естественные магнитные ловушки, к которым относятся магнитосферы Земли и таких больших планет солнечной системы как Юпитер и Сатурн. Характерные времена изменения параметров таких ловушек связаны с медленными конвективными движениями проводящего вещества внутри планет. В атмосферах звезд и Солнца такие движения имеют гораздо меньшие временные масштабы, и времена существования локальных магнитных ловушек, в которых удерживаются частицы, соответственно меньше. В ряде случаев магнитные поля ловушечного типа могут образовываться в солнечном ветре. Как известно, плазма конечного давления является магнетиком и может существовать устойчиво только в областях минимумов магнитного поля. В то же время, отдельные энергичные частицы, составляющие хвост функции распределения, могут быть эффективно захвачены в ловушке с нарастающей к центру напряженностью магнитного поля, и время из жизни в ловушке будет определяться торможением при кулоновском вτзаимодействии с частицами плазмы, перезарядкой с нейтральными атомами и взаимодействием с волнами, генерируемыми как самими частицами (если их концентрация больше определенной величины), так и внешними источниками (менее энергичными частицами плазмы, молниевыми разрядами в атмосфере планеты, промышленными радиопомехами и др.). Для частиц боπльших энергий, для которых ларморовский радиус соизмерим с характерными масштабными неоднородности поля (для оценок можно принять rL~0.1r) магнитный момент перестает быть адиабатическим инвариантом и возникает диффузия по питч-углам, не зависящая от внешних возмущений. Удержание частиц имеет место,
если параметры ловушки стационарны или медленно меняются по сравнения с тремя
характерными временами движения частиц в ловушке: временем оборота по
ларморовской орбите 2π/ωB
(где ωB - циклотронная частота
частицы в магнитном поле B),
τθ
- времени осцилляций частицы между магнитными пробками и τφ -
времени обхода частицы вокруг центра ловушки в результате градиентного дрейфа
за счет кривизны магнитного поля. Как
правило, Последний механизм вносит основной вклад в формирование радиационных поясов Земли и планет – потоков энергичных заряженных частиц в областях захваченной радиации. Распределение потока частиц с определенной энергией в зависимости от расстояния от центра планеты имеет один или несколько максимумов, возникновение которых обусловлено конкуренцией процессов ускорения при сносе частиц в области больших магнитных полей и процессов потерь. На рис. 1.1 показаны зависимости потоков протонов, а на рис. 1.2 - электронов в магнитосфере Земли. Протоны и альфа-частицы попадают в геомагнитную ловушку из солнечного ветра. Наряду с солнечным ветром вторым источником протонов и более тяжелых ионов является земная ионосфера, из которой частицы вытягиваются и ускоряются электрическими полями вдоль магнитных силовых линий. Этот источник действует, в основном, на периферии ловушки и приводит к образованию ионных радиационных поясов Земли. В геомагнитной ловушке при энергиях >40 МэВ начинает играть роль нейтронный распад. Область повышенной интенсивности протонов с энергиями >100 КэВ имеет один максимум (на малых расстояниях от Земли при энергиях >40 МэВ наблюдается небольшой второй максимум). Характерной особенностью зависимости потоков протонов от геоцентрического расстояния является увеличение крутизны внешнего края пояса и приближение максимума к земле с ростом энергии. Дифференциальный спектр протонов вблизи плоскости экватора имеет вид ∂I0/∂ε=S0 exp(-ε/ε0). (1.1) При энергиях больших некоторого критического значения в энергетическом диапазоне 0,1-30 МэВ величина ε0 изменяется в соответствии с сохранением первых двух адиабатических инвариантов. Поэтому, если в качестве независимых переменных выбрать адиабатические инварианты, спектр частиц во внутренних областях геомагнитной ловушки не будет зависеть от геоцентрического расстояния, т.е. источник энергичных частиц радиационных поясов Земли в энергетическом интервале 0,1-30 МэВ находится на периферии ловушки. В отличии от протонов, область повышенной интенсивности электронов с энергией >1000 кэВ распадается на 2 максимума – внешний и внутренний пояса. Зазор между ними, в котором интенсивность на 2-3 порядка ниже чем в максимуме, расположен на геоцентрическом расстоянии 3-3,2 RE. Его положение слабо зависит от энергии. Во внутренней хоне и в зазоре временные вариации на больших высотах невелики. В то же время распределение интенсивности и спектр электронов во внешней зоне сильно изменяется, и амплитуды вариаций достигают трех порядков величины. Несмотря на значительные отличия ионных и электронных радиационных поясов, их объединяет одна важнейшая общая черта: возрастание интенсивности в экваториальной плоскости с уменьшением геоцентрического расстояния связано с увеличением жесткости спектра, а не с увеличением плотности частиц всех энергий, что является характерной особенностью адиабатического ускорения. Перенос частиц с периферии в
центр ловушки осуществляется электрическими полями, которые могут иметь
индукционный (при изменении магнитного поля) или электростатический характер.
Если характерное время изменения электрического поля удовлетворяет условиям τθ«τ<τφ, частиц движутся с сохранением первого и второго адиабатических инвариантов, а третий инвариант – сохраняется. Для того, чтобы происходило радиальное движение частицы, электрическое поле должно иметь долготную, зависящую от долготы, компоненту. Если эти поля имеют случайную природу (что, как правило, реализуется в эксперименте), движение частицы по радиусу будет иметь диффузионный характер (малые неупорядоченные смещения за характерное время изменения поля) и может быть описано в диффузионном приближении. При этом устанавливается некоторое стабильное равновесное распределение, соответствующее усредненному за большой период (характерное время переноса из области малых в область больших магнитных полей) распределению частиц, а резкие изменения параметров источника вызывает медленно движущиеся диффузионные волны. Среди радиационный поясов планет солнечной системы наиболее изучены радиационные пояса Юпитера. Дипольный момент Юпитера в 104 раз больше, чем дипольный момент Земли, а наклон юпитерианского диполя составляет 10о по отношению к оси вращения планеты. Самая внутренняя часть юпитерианской магнитосферы на расстоянии <10 RJ (где RJ - радиус Юпитера) имеет очень близкое сходство с внутренними областями магнитосферы Земли. В этой области зарегистрированы интенсивные потоки энергичных электронов и ионов, стабильно захваченных в дипольноподобном магнитном поле. Значительным отличием от магнитосферы Земли является присутствие 5 внутренних спутников Юпитера, орбиты которых расположены внутри радиационных поясов. Спутники влияют на интенсивность и угловые распределения потоков энергичных частиц как вследствие поглощения, так и инжекции частиц. Во внутренних областях магнитосферы Юпитера радиальная и питч-угловая диффузии являются преобладающими процессами, аналогично тому, как это происходит в радиационных поясах Земли. Ниже будут рассмотрены процессы формирования радиационных поясов на примере радиационных поясов Земли. 2. (L, B) координаты и функция распределения
частиц Напомним, что магнитное поле Земли на расстояниях ≤5RE близко к магнитному полю диполя с магнитным моментом ME=8.02∙1015Тл. В сферической системе координат (R,θ,φ) с центром в диполе и осью, антипараллельной магнитному моменту, напряженность поля на геомагнитной широте (π/2-θ) и расстоянии R от центра диполя B= ME(4-3sin2θ)1/2/R3,
(1,2) а уравнение силовой линии имеет вид R= RELsin2θ, (1,3) где параметр L определяет расстояние от вершины силовой линии до центра Земли, измеренное в земных радиусах RE. В поле диполя зависимость периода долготного дрейфа τφ=2πREL/ vφ (где vφ - скорость долготного дрейфа) от L и энергии частицы ε определяется соотношением τφ=cΨφ(θm)/Lε , (1,4) где константа c слабо зависит от энергии (изменяется в 2 раза при переходе от нерелятивистских к релятивистским энергиям), функция Ψφ(vm) определяет зависимость τφ от широты зеркальной точки (π/2-θ) и изменяется от Ψφ(π/2)=1 на экваторе до Ψφ(0)≈1.5. В случае
земного диполя для частиц с энергиями ~100 МэВ где интегрирование проводится
вдоль магнитной силовой линии, соединяющей магнитосопряженные точки. Т.к. в
стационарном магнитном поле энергия частицы сохраняется, эта величина
пропорциональна При произвольной инжекции частиц
на данную дрейфовую оболочку за время f(L,p,α0)=f0(L,p,arcsin[(sin2αBmin/B)1/2]) . (1,6) В важном частном случае, когда f0=ψ(p)sinvα0, функция f=ψ(p)sinvα(B/Bmin)-v/2, т.е. угловое и энергетическое распределение частиц будут одинаковыми по всей силовой линии, а их плотность изменяется как (B/Bmin)-v/2. При переходе с магнитной оболочки
L на магнитную
оболочку 3. Условия диффузии частиц поперек дрейфовых оболочек Далеко не всякое возмущение
приводит к диффузии частиц поперек дрейфовых оболочек. Если характерное время
возмущения достаточно велико (τ»τφ) или
возмущение имеет аксиально симметричный характер, движение частицы обратимо,
т.е. при исчезновении возмущения частица возвращается на ту же дрейфовую
оболочку. В первом случае сохраняется третий адиабатический инвариант –
магнитный поток через дрейфовую траектория, и при деформациях дрейфовой
оболочки дрейфовая траектория следует за изменениями оболочки. При аксиально
симметричных возмущениях сохраняется проекция обобщенного момента импульса M=[r(p+eA/c)] (где r, p, e- радиус вектор, импульс и
заряд частицы, Третьим условием диффузии является резонанс азимутального движения частиц и возмущений электрического поля, т.е. наличие в спектральных функциях долготных гармоник Фурье электрического поля достаточно больших возмущений на частотах Ωr=mΩD-mvd /LRE, (1,7) где m- номер долготной гармоники, vd- скорость азимутального дрейфа частицы. Получим это условие для случая движения частицы в экваториальной плоскости при наличии возмущений электрического поля. В рассматриваемом случае уравнения движения частицы имеют вид dr/dt=-cEφ/B, dφ/dt=cEr/B
, (1,8) где Er,Eφ - компоненты поля
Е. При малых смещениях за время одного оборота частицы с частотой Ω вокруг
Земли, т.е. при r=r0=const, φ=ΩD t+φ0
, (1,9) где
Разлагая Eφ в ряд Фурье, получаем Если среднее по времени поле <Eφ>=0 а его флуктуации не зависят от времени (стационарный эргодичный процесс), то где Несколько более сложное рассмотрение для флуктуаций геомагнитного поля со спектром мощности PM приводит к выражению 4. Радиальная диффузия под действием внезапных импульсов В магнитосфере постоянно
наблюдаются электрические поля достаточно больших амплитуд, но эти поля
сосредоточены на больших L(≥6-7),
где диффузия частиц по L носит, в основном, электростатический
характер. Во внутренней магнитосфере амплитуды флуктуаций электростатических
полей на частотах, кратных частоте вращения частицы вокруг оси диполя, резко
убывают в результате экранировки внутренней магнитосферы при движении частиц
плазменного слоя. Основной вклад в формирование радиационных поясов Земли дают
быстрые изменения магнитного поля, возникающие при взаимодействии магнитосферы
с межпланетными ударными волнами и разрывами – внезапные импульсы (рис. 1.3). В
результате сжатия или b=-b0(t)LL2sin2θcosφ . (1,15) Остальные квадрупольные члены
вызывают индукционные электрические поля, долготные компоненты которых
антисимметричны относительно плоскости экватора. При этом дрейфовые смещения
частицы в радиальном направлении в противоположных полушариях имеют различные
знаки и взаимно уничтожаются. Члены с более высокой мультипольностью быстро
убывают с уменьшением L.
Для однозначного определения поля Е из условия rotE=c-1∂b/∂t обычно налагается
дополнительное условие об идеальной проводимости вдоль магнитных силовых линий,
эквивалентное условие смещения частиц вместе с силовыми линиями. На рис.1,4
показано движение частиц при возникновении внезапного импульса. Во время
быстрой фазы DLL=D0L10, (1,16) де D0 зависит от энергии и сорта частиц, если их период дрейфа удовлетворяет неравенствам 1мин≤tD≤1час. Экспериментальные определения D0 приводят к D0≈(5-10)∙10-14RE/c. К соотношению (1,16)
приводит достаточно громоздкий рассчет, при котором сначала определяется
индукционное поле Eφ
от поля вида (1,15), а затем рассматривается движение частиц в электрическом и
искаженном магнитном полях. Полученные выражения усредняются по фазам
начального движения частиц от 0 до 2π и находятся среднее где D=0.5<ΔL2/T>, U=<ΔL/T>, Q - мощность источников частиц, G - скорость потерь.
Рассчет дает где D0=0.5α<(Δη)2/T>. Выражение (1,18) содержит регулярный поток частиц к Земле. Это число геометрический эффект, связанный с тем, что область, из которой частицы перейдут внутрь банной оболочки, больше, чем область, занятая уходящими из оболочки частицами. Порядок величины направленного потока ∂L∙∂(δL)/∂L, и, т.к. ∂L~L5, средний поток ~L9. Отсутствие зависимости D0 от энергии и сорта частиц обусловлено в случае внезапных импульсов тем, что на стадии смещения частицы движутся в электрическом поле, скорость дрейфа в котором одинакова для всех частиц, а на последующих стадиях движение обусловлено сохранением третьего инварианта, величина которого также не зависит от характера частицы. Расплывание оболочки складывается из смещения силовой линии в плоскости экватора на дневной стороне δL1 и искажения дрейфовой орбиты при возмущении поля δL2. Смещение δL1 одинаково для всех частиц на данной силовой линии, а δL2 зависит от питч-угла частицы. В случае внезапных импульсов, для которых все гармоники поля сфазированы, D(α=π/2)/D(α=0)≈1.7. В случае хаотических фаз отдельных временных гармоник, коэффициент диффузии для α=π/2на порядок больше, чем для α=0. Значение D0 определяется по прямым измерениям внезапных импульсов и скоростями распространения диффузионных волн. Последние возникают как результат скачкообразны изменений θ на границе диффузионной области. Замена переменной ξ=D0tL8 приводит уравнение (1,17) при Q=G=0 к автомодельному виду. Решение автомодельного уравнения дает одиночную диффузионную волну, распространяющуюся на Земле. Время распространения t~(D0L8)-1, и для L=1,5 время переноса составляет почти 200 лет. Получаемое по данным о диффузионных волнах значение D0 находится в хорошем согласии с данными по внезапным импульсам. 5. Источники и потери частиц поясов Наряду с солнечным ветром, наиболее мощным источником частиц в магнитосфере является ионосфера. Эти источники работают на больших L оболочках и входят в качестве граничного спектра при решении уравнения диффузии. Во внутренней магнитосфере действует источник от распада нейтронов альбедо космических лучей, который может быть задан выражением Q=Q0L-μJ-v. (1,19) Такая зависимость имеет место в случае степенного энергетического спектра и степенной зависимости полного числа рождающихся в единице объема частиц от L. Введение источника (1,19) в уравнение переноса позволяет с достаточной точностью описать самый внутренний протонный пояс, максимум которого находится на L≈1.5. Радиационные пояса в среднем стационарны, а максимумы интенсивности захваченных частиц, как правило, достаточно удалены от Земли. Такая ситуация означает, что имеет место динамическое равновесие между динамическим потоком и утечками частиц. Максимумы частиц определенной энергии лежат там, где время диффузии примерно равно времени жизни. Несмотря на большую угловую анизотропию протонные и ионные радиационные пояса за счет малого числа частиц в трубке, как правило, устойчивы относительно раскачки волн, и гибель частиц при больших энергиях происходит за счет ионизационого торможения, а при меньших (~200кэВ) начинают играть роль процессы перезарядки. Характерное время кулоновского торможения протона пропорционально кубу импульса. Поэтому при сносе вглубь магнитосферы протона с данным J время жизни при постоянной плотности холодной плазмы увеличивается как L-9/2. Так как время переноса растет как L-8, начиная с некоторого L, где время переноса сравнивается со временем торможения, рост интенсивности приостанавливается, а на меньших L интенсивность убывает – формируется максимум пояса. Положение максимумы пояса Lm связано с энергией при постоянной плотности холодной плазмы соотношением Lm~ε3/16. (1,20) При заданной энергии для
нерелятивистской частицы с атомным весом А и зарядом Z время торможения τс(ε)~A-1/2Z-2. В случае α частиц
(А=4, Z=2) время
торможения в 8 раз меньше, чем для протонов той же энергии, а время диффузии не
зависит от энергии. Поэтому максимум пояса Утечка электронов из внешнего радиационного пояса и из зазора на много порядков быстрее ионизационного торможения и связана с несохранением магнитного момента при взаимодействии с волнами. Полная теория такого взаимодействия еще далека от своего завершения. Она должна учетывать вклад черенковского циклотронного резонансов волн с частицами при выполнении резонансных условий ω-k||v|| -nωB , n=0,±1,±2… (1,21) (где ω,k|| - частота и продольная компонента волнового вектора, v|| - продольная компонента скорости частицы, ωB - циклотронная частота) для всех направлений k; генерацию волн за счет анизотропии распределения частиц и процессы их распространения, обратное влияние волн на частиц и наличие внешних источников волнового излучения. Многие особенности распределения электронов в радиальных поясах удалось объяснить в рамках квазилинейной теории циклотронной неустойчивости радиационных поясов с учетом зависимости отражающих свойств ионосферы от частоты волны. Для внутреннего электронного радиационного пояса из-за увеличения минимальной резонансной энергии частиц с уменьшением L не возникает циклотронных резонансных взаимодействий с постоянно наблюдаемыми внутри плазмосферы плазмосферными шипениями, обеспечивающими питч-угловую диффузию и потерю частиц в области зазора и во внешнем поясе, и основную роль играет рассеяние при кулоновских взаимодействиях. В районе максимума протонного пояса интенсивность частиц близка к критической, необходимой для раскачки альвеновских волн. Небольшие вариации потоков протонов или коэффициента отражения волн от ионосферы вызывают генерацию распространяющихся вдоль магнитных силовых линий волновых пакетов, регистрируемых на Земле в виде пульсаций типа Рс1 – жемчужин, взаимодействие с которым вызывает высыпание избытка протонов. |
вперед, гл.2 | ![]() | оглавление | ![]() | литература | ![]() | |||
другие обзоры | ![]() |