Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Fizika_tyazhelykh_ionov

.pdf
Скачиваний:
47
Добавлен:
05.06.2015
Размер:
7.02 Mб
Скачать

тема распадается на невзаимодействующие вторичные адроны. Легко видеть, что скейлинговые решения уравнений релятивистской гидродинамики приводят к

 

равномерному

распределению

 

вторичных частиц по быстротам

 

— плато Фейнмана. Из рис. 9.15

 

очевидны

трудности экспери-

 

ментального

детектирования

 

эффектов

образования

кварк-

 

глюонной плазмы. Действи-

Рис. 9.15. Эволюция ядерной материи

тельно,

время

жизни

КГП

в скейлинговых переменных

τ ≈ τQ , эта величина оказывает-

ся гораздо меньше времени τf .

Поэтому сигналы КГП в значительной степени «маскируются» испусканием частиц из адронной и смешанной фаз.

9.2.3. Начальная плотность энергии термализованной ядерной материи

Рассмотрим центральное (с прицельным параметром, равным нулю) столкновение двух одинаковых ядер с массовым числом A, имеющих в системе центра масс лоренц-фактор γ ( γ = Em >>1, где

m — масса нуклона, E — энергия, приходящаяся на нуклон). В продольном направлении ядра испытывают лоренц-сжатие до раз-

меров 2RA γ ≈ 2A13mπγ , где RA — радиус ядра A и mπ — масса π- мезона. В начальный момент t = 0 нуклоны сталкивающихся ядер взаимодействуют между собой, образуя адронный файербол. В последующие моменты времени файербол расширяется со скоростью порядка скорости света. Оценим начальную плотность энергии файербола в столкновениях релятивистских тяжелых ионов.

Для простоты заменим ядро «нуклоном», движущимся с лоренцфактором γ. При энергиях SPS высота плато в распределении вто-

ричных частиц по быстротам составляет dndych 3 , где nch — число заряженных вторичных частиц.

190

Предполагая, что средняя энергия, приходящаяся на одну час-

тицу в центральной

области быстрот, E 400 МэВ и отношение

числа нейтральных

частиц

к заряженным n0 nch 0.5 ,

находим

плотность энергии в интервале быстрот от y до y + dy :

 

 

d E

3 dnch E 1.8 ГэВ.

(9.46)

 

dy

2 dy

 

Если налетающие частицы не нуклоны, а ядра, то в приближении невзаимодействующих нуклонов предыдущее выражение следует умножить на число нуклонов A. Оценим плотность энергии файербола, толщина которого 2 x , между двумя плоскими слоями,

движущимися со скоростями v = ± xt . Согласно скейлинговому решению, обсуждавшемуся выше, вторичные адроны будут огра-

ничены интервалом быстрот

y = 2arctanh

 

x

2

x

. Тогда энер-

 

 

t

 

 

 

t

 

 

гия частиц

 

 

 

 

 

 

 

 

 

E = A d E

 

 

y .

(9.47)

 

 

 

dy

 

 

 

 

Чтобы найти

начальную плотность энергии εi ,

надо разделить

энергию E на начальный объем:

εi

 

= E Vi . В скейлинговых пере-

менных элемент 4-мерного объема

 

 

 

 

 

 

d 4 x = d 2 x dχτdτ .

(9.48)

 

 

 

t

 

 

 

 

 

Интегрируя d 4 x по χ в пределах ±Δy 2 , получим:

 

 

V = πR2

τ

i

y ,

(9.49)

 

 

i

A

 

 

 

 

где πRA2 — площадь поперечного сечения сталкивающихся ядер

( RA 1.2 A1 3

фм). Тогда начальная плотность энергии файербола

 

εi

A1 3

d E

1 .

(9.50)

 

 

4.5фм

 

dy

τi

 

При столкновении ядер свинца начальная плотность энергии для значений времени τi 1 фм составляет

191

ε

i

5

ГэВ .

(9.51)

 

 

фм3

 

Эта плотность энергии примерно на порядок величины больше нуклонной плотности энергии в ядрах, и как было отмечено в первой главе, при таких плотностях энергии следует ожидать перехода ядерной материи в фазу кварк-глюонной плазмы.

Важной особенностью формулы (9.50) для плотности энергии

является ее пропорциональность A13 . Причина такой зависимости состоит в том, что в начальный момент гидродинамического рас-

ширения вклад в плотность энергии дают A13 пар нуклонов. Важными являются вопросы:

1)в каком объеме образуется ядерный файербол?

2)достаточны ли энергии имеющихся ускорителей для формирования плотностей энергии согласно выражения (9.51)?

Например, для SPS коллайдера ЦЕРН характерные энергии в системе центра масс ~500 ГэВ. Если вся имеющаяся энергия трансформируется в тепловую энергию файербола и плотность энергии достигает значений согласно формуле (9.51), то объем

файербола должен составлять Vi 500 фм3. Эта величина сравнима

с объемом тяжелых ядер (типа свинца). Однако имеющаяся энергия столкновения не полностью трансформируется в тепловую энергию файербола, поскольку заметная доля энергии переносится лидирующими частицами, не входящими в состав файербола. Это обстоятельство понижает оценку плотности энергии примерно в два раза, если вторичные лидирующие частицы несут ~50% энергии сталкивающихся частиц. С другой стороны, для повышения начальной плотности энергии следует отбирать центральные столкновения ядер, т.е. иметь специальный экспериментальный триггер на такие события.

9.2.4.Начальная плотность энтропии

Вэтом разделе попытаемся связать начальную температуру Ti

ядерного файербола с характеристиками вторичных адронов и оценить времена τQ , τH и τf в предположении сохранения энтропии

в процессе расширения.

192

Условие сохранения энтропии означает отсутствие диссипации в процессе расширения. Заметим, что это предположение не означает, что взаимодействие между кварками и глюонами мало. Как известно, условие идеальности ядерной жидкости применимо тогда, когда характеристические времена и длины свободного пробега частиц (кварков, антикварков и глюонов) удовлетворяют условиям:

τ<< t,

(9.52)

λ << L,

где t и L — временные и масштабные факторы, на которых такие гидродинамические величины, как плотность энергии, давление и т. д. испытывают значительные изменения. В какой степени эти условия выполнимы в столкновениях релятивистских ионов? Так как характерная длина свободного пробега кварков λ ≤1 фм, а раз-

мер системы в поперечном направлении RA ~ A13 фм, условие

(9.52) выполняется с хорошей точностью. Характерные времена взаимодействия кварков и глюонов τ ~ 1фм, поэтому первое из условий (9.52) также выполняется с разумной точностью, если учесть, что время t ~ 10 фм — порядка времени взаимодействия ядер. Используя уравнение состояния идеальной кварк-глюонной плазмы, найдем для плотности энтропии

s =

(ε+ p) = dp = 4π2 g T 3 .

(9.53)

 

T

dT

Q

 

 

 

 

Согласно скейлинг-решению, плотность энтропии зависит от собственного времени s(x,t )= s(τ), тогда в момент времени τ пол-

ная энтропия S = s(τ)V (τ), где V (τ)= d

4

Интегрируя

x δ(τ−τ ).

по d 4 x, получим V (τ)= πRA2τ 2Ym . Поэтому

 

 

S = πRA2 2Yms(τ)τ = πRA2 4π2 gQT 3τ 2Ym .

(9.54)

Для того чтобы связать энтропию S с наблюдаемыми величинами, предположим, что кварк-глюонная плазма расширяется в смешанной и адронной фазах адиабатически, т. е. с сохранением энтропии. Так как в адронной фазе пионный газ имеет уравнения состояния (9.1)–(9.3), то по аналогии с (9.54) имеем в этой фазе

S = 4π2 g

H

πR2

2Y T 3τ,

(9.55)

 

A

m

 

193

где τ > τH . В частности, в момент распада файербола на наблюдаемые адроны (τ = τf ):

S 2Y

dS

c dnπ ,

(9.56)

m

dy

dy

 

 

 

где c 3.6 ; dndyπ — распределение по быстротам π-мезонов, рож-

денных в ион-ионных взаимодействиях.

Сравнивая выражения (9.54) и (9.56), получаем соотношение, связывающее плотность энтропии с наблюдаемыми величинами

[9.24]:

s

=

c

 

1

dnπ =

1.2 fm2

dnπ .

(9.57)

 

 

τi A2 3

i

 

τi πRA2 dy

dy

 

Заметим, что плотность энтропии можно определить только в состоянии локального термодинамического равновесия, в то время как плотность энергии можно определить и в неравновесном состоянии. Из выражения (9.57) находим связь начальной температуры с τi :

 

 

 

1.2 fm

2

 

1 3

 

 

T

=

 

dnπ .

(9.58)

 

4π2 g

 

τ

A2 3

 

i

 

Q

dy

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

Время термализации КГП τi

может зависеть, в частности, от

массового числа ядра.

Например,

в параметризации

[9.25]

τiAA

= τipp Aδ , 0 < δ <1 3

и τipp = 0.5 фм. Точность этих оценок вре-

мен

τi составляет порядок величины. Тогда из уравнения (9.58)

следует, что Ti определяется с точностью до фактора два.

Например, при столкновениях ядер 16 O с ядрами Pb при энергиях ЦЕРН [9.24] плотность рожденных π -мезонов составляет

dnπ

42 . Тогда получаем для времени термализации τi = 0.1фм,

 

dy

 

 

 

T 320

МэВ и начальную плотность энергии ε

i

17 ГэВ/фм3.

i

 

 

Время жизни в кварковой фазе τQ = 0.8 фм, Tc =160 МэВ. Для адронной фазы τf 14 фм, Tf mπ 140 МэВ. Как видно из этих

194

оценок, время жизни системы в кварковой фазе гораздо меньше, чем в адронной фазе, поэтому весьма непросто выделить сигналы КГП на фоне смешанной и адронной фаз.

9.2.5.Статистическая модель рождения частиц

встолкновениях тяжелых ионов

Как известно, равновесное поведение термодинамических величин можно получить путем усреднения по статистическому ансамблю. При этом равновесное термодинамическое состояние ансамбля — это такое состояние, в котором плотность в фазовом пространстве однородна. В этом смысле однородное заполнение фазового пространства является необходимым и достаточным условием равновесия.

Основная величина, необходимая для вычисления выходов частиц в столкновениях тяжелых ионов при статистическом подходе

— производящий функционал z (T ,V ) . В каноническом ансамбле

для частицы i, несущей

странность si

,

барионное число Bi ,

электрический

заряд

 

 

Qi

и

спин-изоспиновый

фактор

gi = (2Ji +1)(2Ii +1) производящий функционал

 

 

 

 

 

 

 

Vg

i

 

 

 

 

 

E −μ

i

 

 

 

 

ln z = ±

 

dpp2 ln 1

±exp

 

i

 

 

,

(9.59)

 

2

 

 

 

 

i

2π

0

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где знак “+” соответствует фермионам, а “–” — бозонам.

 

Тогда плотность частиц

 

 

 

 

 

 

 

 

 

 

 

 

 

μi

= μb Bi

s Si 3i I3i ,

 

 

 

 

(9.60)

где T — температура,

E =

p2 + m2

 

полная

энергия,

 

 

 

 

 

i

 

 

i

 

 

 

 

 

 

 

μi = μb Bi s Si

3i I3i

— химический потенциал, связанный с хи-

мическими потенциалами барионного заряда, странности и электрического заряда. При этом предполагается сохранение барионно-

го заряда V ni Bi = Z + N , странности V ni Si = 0 и электриче-

i

i

ского заряда V ni I3i =

Z N .

i

2

195

Эти ограничения оставляют только два свободных параметра модели — температуру T и барионный химический потенциал μb .

На самом деле, конечно, не определен объем системы, и поэтому следует сравнивать отношения выходов различных частиц, которые от объема не зависят.

Очевидно, что большой канонический ансамбль — это простейшая реализация статистического подхода, приемлемая для систем с большим числом рожденных частиц. Однако для небольших систем (например, рожденных в периферических ядро-ядерных столкновениях) или при столкновениях ядер при низких энергиях описание системы каноническим ансамблем является более адекватным. Переход от большого канонического ансамбля к каноническому ансамблю приводит к уменьшению числа частиц («каноническое подавление»). Например, в каноническом подходе хорошо

описываются выходы частиц в e+e-аннигиляции, хотя необходимо ввести феноменологический фактор подавления странности.

На рис. 9.16 представлены отношения выходов различных частиц в Pb+Pb столкновениях при энергиях 158 ГэВ/нуклон (ЦЕРН). Как оказалось, экспериментальные данные хорошо описываются в

статистическом

подходе, если

использовать

параметры

T = (170 ±5) МэВ,

μb = (225 ±10) МэВ.

На рис. 9.17

изображены

отношения выходов частиц в условиях RHIC (при энергиях 130 и

200ГэВ/нуклон).

Вэтом случае данные описываются статистической моделью,

если выбрать параметры T =176 МэВ и μb = 30 ÷40МэВ.

Рис. 9.16. Отношения выходов частиц в Pb+Pb столкновениях при энергиях 158 ГэВ/нуклон (ЦЕРН)

196

Рис. 9.17. Отношения выходов частиц в условиях RHIC

Если нанести извлеченные из экспериментальных данных значения T и μb на фазовую диаграмму ядерной материи, то она бу-

дет выглядеть так, как представлено на рис. 9.18.

На этом рисунке извлеченные из экспериментальных данных значения (T , μb ) изображены точками, приведены также траектории, на которых происходит вымораживание адронного газа, соответствующие постоянной плотности энергии ( ε = 500 МэВ/фм3 ) и

постоянной барионной плотности ( nb = 0.12 фм3 ). Граница между

адронной материей и кварк-глюонной плазмой, рассчитанная в решеточной КХД, обозначена

(LQCD). Квадратики определяют положения критических точек для различных способов вычислений в решеточной теории [9.26].

Как видно из рис. 9.18 при энергиях SPS извлеченные из экспериментальных данных параметры

(T ,μb ) близки к фазовой границе,

полученной в решеточных вычислениях. Означает ли это, что уже при энергиях SPS достигается равновесие? Как отмечено в работе [9.27] скорости столкновений частиц и временные масштабы рас-

Рис. 9.18. Фазовая диаграмма по данным термальной модели

197

ширения адронной материи не позволяют сделать вывод об установлении равновесия при энергиях SPS, хотя наличие фазового перехода при этом не исключается.

9.3. Корреляции между поперечной энергией и прицельным параметром в столкновениях релятивистских тяжелых ионов

9.3.1. «Жесткая» и «мягкая» компоненты поперечной энергии

Как видно из (9.50), начальная плотность энергии пропорцио-

нальна A13 для центральных столкновений ядер (с прицельным параметром, равным нулю). Для нецентральных соударений зависимость плотности энергии от массового числа гораздо слабее, чем

A13 . Таким образом, наибольших εi следует ожидать именно в

центральных соударениях. Но как экспериментально отбирать эти соударения, ведь прицельный параметр не является величиной, непосредственно измеряемой в столкновениях ионов?

В этом случае нужно найти измеряемую величину, которая бы коррелировала с прицельным параметром. В этой главе обсуждается метод косвенного измерения прицельного параметра по попе-

речной энергии в столкновении E = E i ,

где

суммирование

производится по всем частицам в событии E

i

=

m2

+ p2

,

m и

 

 

i

i

 

i

p i — масса и поперечный импульс i -й частицы.

Отметим, что кроме измерения поперечной энергии в событии, существуют и другие методы определения прицельного параметра. В работе [9.28] предлагается извлекать прицельный параметр путем регистрации ядерных осколков после развала ядра в калориметре, установленном под углом M R к оси сталкивающихся пуч-

ков (ZDC калориметр).

В формирование поперечной энергии вносит вклад как «жесткая», так и «мягкая» составляющие в столкновении. «Жесткая» компонента обусловлена столкновениями кварков и глюонов с характерными переданными импульсами, превышающими p0 12 ГэВ/с — это масштаб, при котором применима теория

возмущений КХД. «Мягкие» процессы с характерными p < p0 ,

198

также дающие вклад в поперечную энергию, не описываются в рамках пертурбативной КХД, и их учет должен проводиться феноменологически. Будем предполагать, что «мягкая» часть протон-

протонного сечения σspp = σнеупрpp , где Ι +Ι — сечение неупругого

pp-взаимодействия. Если жесткая компонента сечения в столкновениях ядер A и B формируется за счет независимых партон-

партонных взаимодействий, то среднее число жестких столкновений как функция прицельного параметра b представляется в виде суммы по всем типам партонов f (кваркам, антикваркам и глюонам):

 

 

 

 

 

 

ABH (b)= N

ABf

(b, sNN , p0 ),

 

 

(9.61)

 

 

 

 

N

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dσf

s

, p

 

 

 

 

 

 

 

ABf (b, sNN , p0 )= 2

AB (

NN

0

)

,

(9.62)

 

 

N

 

 

d

2

 

 

 

dσf

 

 

 

 

 

 

b

 

 

 

 

и сечение

 

 

рождения партона сорта f вычисляется в теории

AB

 

 

d 2b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

возмущений КХД [9.29]. Тогда распределение по поперечной энергии E можно записать в виде:

 

H

 

 

 

 

 

 

H

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dσ

 

 

2

 

 

 

NAB (b)

 

 

 

H

 

 

 

= d b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dE

 

 

 

 

N !

exp NAB

 

 

N =1

 

 

 

 

 

(9.63)

N

 

1

d

σpp

 

 

N

 

 

 

 

∫∏dE i

 

 

 

 

H δ E E i .

 

 

pp

 

 

 

i=1

 

σH

dE i

 

 

i=1

 

 

Как показано в работе [9.30], при числе жестких столкновений NABH 5 ÷7 , распределение (9.63) при фиксированном значении

прицельного параметра хорошо аппроксимируется распределением Гаусса:

 

 

 

 

 

 

 

(E

 

ABH (b))2

 

 

dσH

= d

2

b

1

 

E

 

 

 

 

exp

2σ2H (b)

,

(9.64)

dE

 

2πσ2H (b)

 

 

 

 

 

 

 

 

 

 

 

 

199

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]