Коллективные Е1-возбуждения в ядре отвечают колебанию
электрического дипольного момента ядра, т.е. колебанию всех протонов
относительно всех нейтронов. Не случайно, поэтому, первые попытки интерпретации
ГДР были связаны с коллективными моделями. Причем, без детальных расчетов было
ясно, что энергия Е1-возбуждений, обусловленных движением всех нуклонов ядра,
должна быть во много раз больше, чем энергия поверхностных колебаний ядра (< 1
МэВ), в которых участвует лишь несколько нуклонов, находящихся на поверхности.
Первой коллективной (и вообще первой) моделью ГДР явилась
модель
Мигдала
[6],
давшая правильное объяснение этому явлению. В этой модели ядро рассматривается
как совокупность взаимопроникающих сжимаемых протонной и нейтронной жидкостей.
Внешнее электрическое поле вызывает колебания протонной жидкости относительно
нейтронной и одновременно изменение их плотностей. Ниже мы более подробно
рассмотрим гидродинамический подход Мигдала в том варианте, который позже
развили Штейнведел и Йенсен [7]
и Данос [8],
а перед этим коснёмся модели
Гольдхабера
и
Теллера
[9], в которой в процессе электрических дипольных
колебаний несжимаемые протонная и нейтронная сферы двигаются как целые друг
относительно друга (так как это изображено на рис. 2).
Резонансная частота такого осциллятора может быть оценена из
классического соотношения
(13) |
где k - коэффициент упругости, а m - масса осциллятора. Упругую силу создает притяжение сдвинутых (неперекрывающихся) частей протонной и нейтронной сфер с остальной частью ядра. При малых смещениях протонной и нейтронной сфер число неперекрывающихся (оголенных) нуклонов, а, следовательно, и коэффициент упругости пропорциональны площади поверхности ядра, т.е. R2. Масса ядра пропорциональна R3. Отсюда получаем
(14) |
Гольдхабер и Теллер нашли для константы в этом выражении значение 45 МэВ и, таким образом, в рассматриваемой модели для энергии ГДР получается следующая формула
Em = 45.A-1/6 МэВ, | (15) |
которая, во-первых, дает слишком медленное падение энергии ГДР с ростом
массового числа, а, во-вторых, сильно (в 1.5-2 раза для средних и тяжелых ядер)
занижает эту энергию. Отметим, что максимальное смещение протонной сферы
относительно нейтронной (амплитуда колебаний) оказывается небольшим, всего лишь
≈ 0.1 Фм.
Значительно более успешной в объяснении ГДР оказалась
гидродинамическая модель двухкомпонентной ядерной жидкости Мигдала, Штейнведела
и Йенсена, к которой мы и переходим.
В этой модели нейтроны и протоны рассматриваются как две
взаимопроникающие сжимаемые жидкости, находящиеся внутри фиксированной ядерной
поверхности. В основном состоянии обе жидкости равномерно распределены по объему
ядра, так как ядерные силы благоприятствуют однородному перемешиванию нейтронов
и протонов. Электромагнитное поле возмущает это равновесное распределение,
поэтому при поглощении γ-кванта возникают колебания плотностей нейтронной и
протонной жидкостей (поляризационные колебания), которые и проявляются как
широкие пики (гигантские резонансы) в сечениях реакций под действием фотонов.
Нейтрон-протонные поляризационные (изовекторных) колебания
можно трактовать как колебания сжатия-расширения нейтронной и протонной
жидкостей, которые колеблются в противофазе, так что общая ядерная плотность
ρ(,t) = ρp(,t) + ρn(,t) = ρ0 | (16) |
остается постоянной в любой точке ядра (здесь ρp(,t)
– плотность протонной жидкости, ρn(,t)
– плотность нейтронной жидкости и ρ0 – константа,
характеризующая полную плотность). Разделение протонной и нейтронной жидкостей
под действием внешнего поля описывается флуктуацией плотности η(,t)
. При этом
ρp(,t)
= ρp(0)[1 + η(,t)], ρn(,t) = ρn(0)[1 - (Z/N)η(,t)] |
(17) |
где ρp(0) и ρn(0) = (N/Z)ρp(0) - равновесные протонные и нейтронные плотности, а η( ,t) играет роль отклонения плотности от равновесного значения. Видно, что для любого момента времени выполняется условие (16), отвечающее пренебрежению компрессионными колебаниями (сжатия и разряжения ядра). η( ,t) можно представить в виде
η(,t) = η()e-iωt, | (18) |
где ω - частота колебаний плотности (предполагается, что колебания гармонические). Тогда колебания протонной и нейтронной плотностей будут удовлетворять обычному волновому уравнению гидродинамики. При этом η( ) подчиняется уравнению Гельмгольца
Δη() + k2η() = 0 | (19) |
где Δ - лапласиан, с граничным условием
.grad[η()]r = R = 0, | (20) |
означающим отсутствие потока нуклонов через поверхность ядра (
- вектор, нормальный поверхности ядра).
Итак, колебания плотности описываются волновым уравнением.
Волны плотности - своего рода звуковые волны в ядерной материи. Волновое число k
связано с константой энергии симметрии β = 23.6 МэВ в формуле Вайцзеккера
(21) |
В свою очередь
(22) |
где
(23) |
- скорость распространения в ядре колебаний плотности (аналог
скорости звука в ядерной среде).
Нормальные решения волнового уравнения (19), удовлетворяющие
граничному условию (20), имеют вид
(24) |
где J - мультипольность колебаний, - сферические функции, - сферические функции Бесселя, а - нормировочные константы. Волновые числа получаются из граничного условия (20), которое теперь принимает вид
(25) |
Решения этого трансцендентного уравнения получаются численно. Они даны в
табл. 2 до четвертого (n = 4) обертона мультипольности J = 4 (в таблице
приведены
R).
Энергии колебаний даются соотношением
(26) |
Таблица 2. Решения R уравнения (25) |
n | J = 0 | J = 1 | J = 2 | J = 3 | J = 4 |
1 | 4.493 | 2.081 | 3.342 | 4.514 | 5.646 |
2 | 7.725 | 5.940 | 7.289 | 8.583 | 9.840 |
3 | 10.904 | 9.205 | 10.613 | 11.972 | 13.295 |
4 | 14.066 | 12.404 | 13.846 | 15.244 | 16.609 |
Из (26) для энергии основных (n = 1) поляризационных колебаний, характеризующихся самыми низкими угловыми моментами J = 0, 1 и 2, получаем (сравни с данными табл. 1):
(27) |
Оценка сделана для N/Z ≈ 1.3, что соответствует
типичному ядру с А ≈ 100. Наименьшую собственную энергию имеет
дипольная поляризационная мода ( J = 1). Все остальные поляризационные моды
– квадрупольная ( J = 2), монопольная (J = 0), а также октупольная ( J = 3)
и др. располагаются значительно выше по шкале энергий. Если оценить энергию
дипольных колебаний по формуле (27) для ядра с массовым числом А =
100, то получится 16.4 МэВ. Более высокие (n > 1) обертоны для всех J
располагаются при существенно более высоких энергиях и возбуждаются с
существенно меньшей вероятностью. Так, основная Е1-мода ( n = 1) вбирает 86%
вероятности электрических дипольных возбуждений. На Е1-обертоны с n
= 2 и 3, располагающиеся при энергиях ≈ 47 МэВ
и 72 МэВ, приходится соответственно 6% и 2% вероятности Е1-переходов.
Отметим, что указанные Е1-обертоны на языке оболочечных переходов (см. рис.
3) отвечают возбуждениям с энергиями соответственно 3 ћω и 5 ћω .
На рис. 6 показано, как примерно выглядят поляризационные
колебания нижайшей мультипольности для ядра с А ≈ 60 в крайней фазе
колебаний, отвечающей максимальному разделению протонной и нейтронной
жидкостей (протоны и нейтроны по разному окрашены). Для наглядности эффект
разделения протонной и нейтронной жидкостей сильно утрирован. На самом деле
максимальное их разделение (амплитуда поляризационных колебаний) составляет,
как мы отмечали в связи с рассмотрением модели Гольдхабера-Теллера, всего
лишь около 0.1 Фм.
Отметим также, что рассмотренные коллективные модели
применимы к достаточно массивным ядрам, содержащим не менее 50 нуклонов.
Рис. 6. Поляризационные колебания нижайшей мультипольности. Показана крайняя фаза колебаний, отвечающая максимальному разделению протонной и нейтронной жидкостей (протоны и нейтроны по разному окрашены)
Как видно из рис. 4 а, гидродинамическая модель, давая для энергии ГДР величину Еm ≈ 75А-1/3 МэВ, в целом довольно успешно объясняет зависимость Еm от массового числа (по крайней мере в области А ≈ 50-150). Более того, эта модель обеспечивает естественное и безупречное объяснение формы гигантского резонанса несферических (деформированных) ядер. Для таких ядер, имеющих, как правило, форму аксиально-симметричного эллипсоида, сечение поглощения Е1-фотонов должно иметь два гигантских резонанса, а не один, как в случае сферических ядер. Это непосредственно следует из зависимости Еm от радиуса ядра
Еm ~ A-1/3 ~ 1/R. | (28) |
Для сферического ядра наблюдается одна резонансная частота электрических дипольных колебаний. Если же ядро несферическое (например, аксиально-симметричное), то оно имеет два характерных размера b и а - длины полуосей ядерного эллипсоида соответственно вдоль и перпендикулярно оси симметрии ядра. В соответствии с этим у такого ядра существует две резонансные частоты дипольных колебаний (см. также рис. 7)
Ea ≈ 75r0/a МэВ, Eb ≈ 75r0/b МэВ, | (29) |
где r0 = 1.2 Фм. Таким образом,
Ea/Eb ≈ b/a. | (30) |
Более детальное рассмотрение Даноса [8] даёт
Ea/Eb ≈ 0.911b/a + 0.089. | (31) |
Если, например, ядро имеет положительный электрический квадрупольный момент, т. е. вытянуто вдоль оси симметрии, то для такого ядра b > а и дипольные колебания вдоль короткой оси совершаются с большей частотой, чем вдоль длинной. Различие в энергиях ΔЕm (частотах) этих колебаний можно получить из следующего соотношения
ΔЕm = Еm(b) - Еm(a) = 75r0(1/a - 1/b) МэВ ≈ 75r0β/<R> МэВ ≈ 75A-1/3β МэВ, | (32) |
где β = (b - a)/<R> - параметр деформации ядра, а величину <R> = (ab)1/2 можно оценить по формуле радиуса ядра <R> ≈ r0 A-1/3 с r0 = 1.2 Фм.
Рис. 7. Электрические дипольные колебания вдоль короткой и длинной осей ядерного эллипсоида