Аннотация:
В рамках модели, учитывающей слияние кварк-глюонных струн и образование струнных кластеров, для случая pp-рассеяния при высоких энергиях рассчитана сильноинтенсивная переменная Σ, характеризующая корреляции между числом частиц, образующихся в двух разнесенных по быстроте интервалах наблюдения. Проведено сравнение результатов расчетов с экспериментальными данными коллаборации ALICE на Большом адронном коллайдере. Показано, что экспериментально наблюдаемый рост этой переменной с увеличением начальной энергии можно объяснить только с учетом образования струнных кластеров, состоящих из возрастающего числа слившихся струн.
Ключевые слова:
взаимодействия адронов, высокие энергии, множественное рождение частиц, кварк-глюонные струны, слияние струн, кластеры, флуктуации, корреляции, cильноинтенсивные переменные.
Известно, что изучение флуктуаций и корреляций различных наблюдаемых в процессах множественного рождения частиц при высоких энергиях дает информацию о начальном этапе взаимодействия адронов, отвечающем наибольшей плотности образующейся кварк-глюонной материи [1]. Это позволяет получить уникальные данные о новых объектах, которые образуются на этом этапе, в том числе о кварк-глюонных струнах и струнных кластерах.
Поскольку данные физические явления относятся в основном к так называемой мягкой области сильного взаимодействия с малыми передачами импульса, в которой невозможно провести расчеты по теории возмущения КХД, то обоснованным является использование модели кварк-глюонных (цветных) струн [2], [3]. В этом подходе используется имеющее качественное обоснование в рамках КХД представление о струнах как о трубках цветового потока (см. более подробно [4]–[8]).
Дальнейшие исследования показали необходимость учета процессов взаимодействия и слияния струн [9]–[11]. Эти процессы были включены в генератор событий DIPSY в виде формирования так называемых “цветных веревок”, приводящих к увеличенному выходу странных частиц [12]. В статье [13], посвященной экспериментальным исследованиям выхода мультистранных адронов в pp-взаимодействиях как функции центральности столкновения, при сравнении этих результатов с предсказаниями различных теоретических моделей коллаборация ALICE делает вывод, что модель DIPSY с учетом процессов образования “цветных веревок” описывает данные лучше, чем другие генераторы событий. Аналогичные результаты [14] были получены и в рамках мультипомеронной модели, в которой эффективным образом также учтен вклад процессов слияния струн [15], [16].
Чтобы извлечь информацию о процессах слияния струн и образовании струнных кластеров из анализа флуктуаций и корреляций между различными наблюдаемыми, необходимо уметь отделить в экспериментальных данных тривиальный затеняющий вклад так называемых “объемных” флуктуаций, возникающих из-за флуктуаций объема области взаимодействия (например, числа образующихся струн) от события к событию, в частности по причине принципиально неустранимых флуктуаций прицельного параметра в реальном эксперименте. Поэтому особое значение имеет изучение флуктуаций и корреляций с использованием интенсивных и сильноинтенсивных наблюдаемых, которые позволяют минимизировать вклад “объемных” флуктуаций и дают возможность получить информацию о свойствах объектов, образующихся на начальной стадии сильного взаимодействия. В настоящей работе мы используем для этой цели так называемую сильноинтенсивную переменную Σ, введенную в работах [17], [18]. Мы рассчитываем поведение этой переменной в рамках модели pp-взаимодействия со слиянием струн и образованием струнных кластеров для различных начальных энергий.
Сравнение полученных теоретических зависимостей для переменной Σ с предварительными экспериментальными данными [19] коллаборации ALICE на Большом адронном коллайдере (БАК) при трех начальных энергиях рр-столкновений позволяет сделать однозначный вывод о наличии эффектов слияния струн и образования струнных кластеров при энергиях БАК уже в рр-взаимодействиях. При этом, как и ожидалось, вклад струнных кластеров в переменную Σ растет с ростом начальной энергии. Сравнение с экспериментальными данными позволяет также фиксировать параметры, характеризующие кластеры с различным числом слившихся струн.
2. Сильноинтенсивные переменные
В статистической физике используемые величины обычно делят на интенсивные и экстенсивные. К первым относятся давление, температура и тому подобные величины, значение которых остается неизменным при выделении из рассматриваемой системы какой-то ее части. Ко вторым относят объем, энтропию и другие величины, значение которых равно сумме их значений для отдельных частей системы.
При перенесении этих понятий на физику высоких энергий, например на случай ядро-ядерных столкновений, объем образующейся системы считается пропорциональным суммарному числу нуклонов в сталкивающихся ядрах, принявших участие в процессе взаимодействия. Ясно, что в этом случае объем зависит от прицельного параметра, при котором происходит столкновение ядер. В более общем случае взаимодействия адронов высоких энергий считается, что объем образующейся системы пропорционален числу N кварк-глюонных струн, формирующихся на начальном этапе взаимодействия.
С этой точки зрения число частиц n с импульсом, принадлежащим некоторой заданной области (в заданном аксептансе), представляет собой пример экстенсивной величины, так как ее значение пропорционально объему образующейся системы. В каждом столкновении адронов значение n равно сумме вкладов, n=n1+⋯+nN, от каждой из N образующихся струн. Очевидно, что величина n флуктуирует от события к событию вместе с величиной объема образующейся системы (числом образующихся первичных струн), в том числе из-за неизбежных флуктуаций прицельного параметра.
Однако, как отмечалось в разделе 1, для нахождения вклада физически интересных флуктуаций и корреляций, несущих информацию о фундаментальных свойствах образующейся в результате взаимодействия адронов материи, необходимо уметь выделять их вклад на фоне тривиального вклада, возникающего из-за флуктуаций объема образующейся системы. Для этого в физике высоких энергий, наряду с экстенсивными и интенсивными величинами, вводится понятие сильноинтенсивных переменных, под которыми понимаются такие наблюдаемые величины, значение которых не зависит не только от объема образующейся системы, но и от флуктуаций этого объема от события к событию.
Общие методы построения таких наблюдаемых исследовались в работе [17], где, в частности, было показано, что в рамках определенного класса статистических моделей, рассмотренных в работе, величина
Σ(A,B)≡⟨A⟩ωB+⟨B⟩ωA−2cov(A,B)⟨A⟩+⟨B⟩,
составленная из любых двух экстенсивных величин A и B, является сильноинтенсивной переменной. В этой формуле ⟨A⟩ и ⟨B⟩ – средние значения величин A и B, а ωA и ωB – их приведенные дисперсии:
ωA≡DA⟨A⟩=⟨A2⟩−⟨A⟩2⟨A⟩,ωB≡DB⟨B⟩=⟨B2⟩−⟨B⟩2⟨B⟩.
В формулу (1) входит также коррелятор этих переменных:
cov(A,B)=⟨AB⟩−⟨A⟩⟨B⟩.
Позднее в работе [18] было предложено при построении переменной Σ по формуле (1) использовать в качестве этих экстенсивных переменных A и B значения числа частиц nF и nB, регистрируемых в данном событии в двух интервалах импульсов этих частиц, в так называемых окнах наблюдения. Часто эти окна наблюдения обозначаются как переднее (forward) и заднее (backward), так как в первых экспериментах им соответствовали частицы, вылетающие соответственно в переднюю и заднюю полусферы реакции в системе центра масс.
При рассмотрении процессов множественного рождения частиц при высоких энергиях удобно перейти от продольных, направленных вдоль оси столкновения, компонент pz импульсов частиц к быстротам этих частиц:
y=12lnp+p−=12lnp0+pzp0−pz,
где p0=√p2z+p2⊥. В нашем рассмотрении мы считаем, что частица попадает в передний интервал наблюдения, если ее быстрота принадлежит заданному интервалу: y∈(yF,yF+δyF), поперечный импульс p⊥ при этом может быть произвольным, и аналогично для заднего окна наблюдения. Соответственно δyF и δyB – ширины этих окон наблюдения по быстроте.
В дальнейшем мы ограничимся случаем симметричных реакций, к которым, в частности, относится изучаемое pp-взаимодействие, и одинаковых, δyF=δyB≡δy, симметрично расположенных относительно y=0 быстротных интервалов наблюдения. В этом случае
⟨nF⟩=⟨nB⟩≡⟨n⟩,DnF=DnB≡Dn,ωnF=ωnB≡ωn,
и выражение (1) для этого частного случая принимает вид
Σ(nF,nB)=Dn−cov(nF,nB)⟨n⟩=⟨n2⟩−⟨nFnB⟩⟨n⟩.
3. Связь наблюдаемой Σ(nF,nB) со свойствами струнных кластеров
В работе [26] мы получили выражение для переменной Σ(nF,nB) в рамках модели, учитывающей процессы слияния струн и образование струнных кластеров. В этой работе для учета эффектов от слияния струн мы использовали решетку в плоскости прицельного параметра (см. более подробное обсуждение этого метода ниже в разделе 5). Было показано, что в этом случае
Σ(nF,nB)=∞∑k=1α(k)Σk(μF,μB),
где
α(k)≡⟨n⟩k⟨n⟩.
В этих формулах переменная Σk(μF,μB) определена согласно общим выражениям (1) и (6), но для какого-то одного кластера, образованного слиянием k струн. Соответственно, величины μF и μB – это числа частиц, которые попадают в переднее и заднее окна наблюдения в данном событии от распада этого кластера. Разумеется, экспериментально их невозможно выделить из общего числа частиц nF и nB, регистрируемых в этих окнах в данном событии.
Весовые множители α(k) имеют смысл средней доли частиц, рожденных от всех кластеров, образованных слиянием ровно k струн; ⟨n⟩k – среднее число частиц, рождающихся от всех таких кластеров. Поэтому выполняется нормировочное условие
∞∑k=1α(k)=1⟨n⟩∞∑k=1⟨n⟩k=1.
Значению k=1 отвечает вклад одиночных струн.
Из формул (7) и (8) мы видим, что без учета процессов слияния струн вклад дают только одиночные струны, и мы имеем Σ(nF,nB)=Σ1(μF,μB). В этом случае подтверждается сильноинтенсивный характер переменной Σ(nF,nB) – ее значение полностью определяется только свойствами одиночной струны и не зависит от числа N образующихся в данном столкновении струн и флуктуаций этого числа от события к событию.
В то же время в общем случае с учетом слияния струн и формирования новых типов источников частиц – струнных кластеров – переменная Σ(nF,nB), как мы увидим ниже, оказывается зависящей не только от свойств струнных кластеров, определяющих величины Σk(μF,μB), но и от деталей столкновения – его начальной энергии и центральности, через весовые коэффициенты α(k).
В работе [20] было найдено следующее выражение для переменных Σk(μF,μB), входящих в формулу (7), через величины, характеризующие свойства струнных кластеров, образованных слиянием k струн:
Σk(μF,μB)=1+μ(k)0δy[J(k)−J(k)FB],
где
J(k)=1(δy)2∫δydy1∫δydy2Λ(k)(y1−y2)
и
J(k)FB=1(δy)2∫δyFdy1∫δyBdy2Λ(k)(y1−y2).
В последней формуле при интегрировании переменная y1 принадлежит переднему, а переменная y2 – заднему окну.
Здесь мы приводим эти формулы для случая симметричных по быстроте окон наблюдения одинаковой ширины δyF=δyB=δy, выбранных в центральной области быстрот. Мы также предполагаем наличие трансляционной инвариантности по быстроте в этой области. Отметим, что такое приближение хорошо работает при энергиях БАК, когда формирующиеся струны дают вклад в достаточно широкий интервал быстроты.
В этом приближении одночастичные распределения λ(k)(y) частиц от распада данного кластера, сформированного слиянием k струн, постоянны, а двухчастичные распределения λ(k)2(y1,y2) зависят только от разности быстрот образующихся частиц:
λ(k)(y)=μ(k)0,λ(k)2(y1,y2)=λ(k)2(y1−y2),
где μ(k)0 – среднее число частиц, рождающихся на единицу быстроты от распада этого струнного кластера.
Парная корреляционная функция кластера Λ(k) в формуле (10) определена стандартным образом:
Λ(k)(y1,y2)≡λ(k)2(y1,y2)λ(k)(y1)λ(k)(y2)−1.
При наличии трансляционной инвариантности по быстроте это определение принимает вид
Λ(k)(y1−y2)=λ(k)2(y1−y2)(μ(k)0)2−1.
Особенно простой вид формула (10) принимает в случае быстротных окон наблюдения малой ширины, когда δyF=δyB≡δy≪ycorr, где ycorr≃1÷2 – характерная длина корреляции в пространстве быстрот. В этом случае формула (10) принимает следующий вид:
Σk(μF,μB)=1+μ(k)0δy[Λ(k)(0)−Λ(k)(Δy)],
где Δy – расстояние по быстроте между центрами этих малых окон наблюдения δyF и δyB.
Эта упрощенная формула позволяет качественно понять основные черты зависимости переменной Σk(μF,μB) для данного кластера из k струн от ширины окон наблюдения δy и расстояния по быстроте между ними Δy. При малом расстоянии между окнами наблюдения, когда Δy≪ycorr, мы видим, что Σk(μF,μB)→1. С ростом Δy значение Σk(μF,μB) возрастает, выходя на уровень 1+μ(k)0δyΛ(k)(0), так как Λ(k)(Δy)→0 при Δy≫ycorr. При этом скорость выхода на этот уровень и его величина пропорциональны ширине окон наблюдения δy.
Это значение 1+μ(k)0δyΛ(k)(0), на которое выходит переменная Σk(16) для удаленных друг от друга по быстроте окон наблюдения, есть просто приведенная дисперсия ω(k)μ≡D(k)μ/⟨μ⟩(k), ⟨μ⟩(k)=μ(k)0δy, числа частиц, попадающих при распаде струнного кластера в быстротный интервал ширины δy. Действительно, в работах [21], [22] в общем случае источника любого типа было показано, что ω(k)μ дается первым слагаемым в формулах (10) и (16):
ω(k)μ=1+μ(k)0δyJ(k)≈1+μ(k)0δyΛ(k)(0),
где J(k) определяется по формуле (11), а последний переход справедлив, когда ширина окна много меньше характерной корреляционной длины по быстроте, δy≪ycorr≃1÷2.
Важным физическим следствием формулы (17) является неизбежное отклонение распределения по числу частиц в данном фиксированном интервале быстроты δy от распределения Пуассона, для которого ω(k)μ должно быть равно единице, при наличии любых, даже ближних, корреляций между частицами, образующимися в процессе распада данного кластера. Из формулы (17) видно также, что распределение Пуассона по числу частиц от распада кластера в данном интервале быстроты δy возможно только для очень узких окон, когда δy⋅μ(k)0Λ(k)(0)≪1.
4. Моделирование распределения струн в поперечной плоскости pp-столкновения
Расчеты величины Σ(nF,nB) для pp-столкновений при энергиях БАК проводились на основе формулы (7). Остановимся сначала на расчете весовых коэффициентов α(k)(8), входящих в эту формулу и определяющих соотношение вкладов кластеров с различным числом слившихся струн. Ясно, что это соотношение зависит от конкретных деталей pp-столкновения, в частности от начальной энергии и от степени его центральности.
Расчет весовых коэффициентов α(k) в настоящей работе проводился в два этапа. Сначала проводилось моделирование распределения первичных струн в поперечной плоскости с учетом реальных условий pp-столкновения согласно методике, предложенной в работе [23]. На втором этапе мы моделировали процесс слияния первичных струн и образования струнных кластеров, вводя конечную решетку (грид) в плоскости прицельного параметра (см. раздел 5).
Чтобы иметь возможность учесть эти процессы наложения струн и образования струнных кластеров, мы должны уметь на первом этапе моделировать не только общее количество первичных струн, образующихся при данной начальной энергии и заданном прицельном параметре pp-столкновения, но и их распределение в поперечной плоскости. Mы предполагаем, что и в случае столкновения двух протонов при прицельном параметре b плотность распределения струн в точке s поперечной плоскости пропорциональна произведению профильных функций сталкивающихся протонов:
wstr(s,b)∼1σpp(b)T(s−b2)T(s+b2),
где T(s) – партонная профильная функция нуклона, а σpp(b) – вероятность недифракционного pp-взаимодействия (по крайней мере, с одним рассеченным помероном) при заданном значении прицельного параметра b. Интеграл
σpp=∫σpp(b)d2b
дает полное сечение недифракционного pp-взаимодействия.
По аналогии с легкими ядрами для партонной профильной функции протона мы будем использовать простейшее гауссово распределение с параметром r0:
Видно, что в этом приближении зависимость от переменных s и b факторизуется. Интегрируя (21) по s, мы находим, что для среднего числа струн, образующихся в pp-столкновении на прицельном параметре b, справедливо следующее соотношение:
⟨Nstr(b)⟩∼e−b2/2r20σpp(b).
Поскольку один рассеченный померон приводит к формированию двух кварк-глюонных струн, то Nstr=2N, где N – число рассеченных померонов в данном событии. Поэтому в силу (22) мы также имеем
⟨N(b)⟩∼e−b2/2r20σpp(b).
Эта формула дает зависимость среднего числа рассеченных померонов от прицельного параметра для недифракционных pp-столкновений.
Для расчета величины флуктуаций и корреляций различных наблюдаемых требуется знание не только среднего значения числа рассеченных померонов в pp-столкновениях при заданном значении прицельного параметра b, но и распределения числа померонов от события к событию вокруг этого среднего значения. Мы будем предполагать, что при заданном прицельном параметре b это распределение ˜P(N,b), которое при N⩾1 дается формулой
˜P(N,b)=P(N,b)1−P(0,b),
является простой модификацией распределения Пуассона:
P(N,b)=1N!e−¯N(b)¯N(b)N
с некоторым параметром ¯N(b). Отличие нашего распределения ˜P(N,b)(24) от пуассоновского (25) состоит только в исключении из него значения N=0: ˜P(0,b)=0, что соответствует исключению дифракционного рассеяния, которому соответствует N=0. Ясно, что при N⩾1 это сводится лишь к введению в (24) дополнительного общего множителя, обеспечивающего правильную нормировку: ∑N=1˜P(N,b)=1.
Расчет среднего числа померонов при заданном b с распределением (24) дает
⟨N(b)⟩=¯N(b)1−P(0,b).
Учитывая, что вероятность σpp(b) недифракционного pp-взаимодействия при заданном фиксированном прицельном параметре b равна вероятности иметь ненулевое число рассеченных померонов, мы имеем
σpp(b)=1−P(0,b)=1−e−¯N(b).
Сравнивая с учетом этого формулы (23) и (26), приходим к выводу, что ¯N(b)∼e−b2/2r20, или, вводя параметр N0, имеем
¯N(b)=N0e−b2/2r20.
Тогда, с учетом формул (25) и (26) среднее число рассеченных померонов при прицельном параметре b определяется выражением
Чтобы получить результат для pp-столкновений без отбора по центральности, необходимо выполнить интегрирование по b. Для этого удобно в плоскости прицельного параметра ввести плотность вероятности недифракционного pp-взаимодействия (см. (19)), нормированную на единицу:
f(b)=σpp(b)σpp,∫f(b)d2b=1.
Тогда, используя (30), можно найти среднее число померонов в недифракционном pp-взаимодействии, усредненное по прицельному параметру:
⟨N⟩=∫⟨N(b)⟩f(b)d2b=∫¯N(b)σppd2b=2πr20N0σpp,
а также соответствующую дисперсию DN≡⟨N2⟩−⟨N⟩2, где
⟨N2⟩=∫⟨N2(b)⟩f(b)d2b=πr20N0(N0+2)σpp=⟨N⟩(1+N02).
Здесь мы использовали формулу (28). Используя эту формулу, мы можем также вычислить в нашей модели полное сечение недифракционного pp-взаимодействия:
σpp=∫σpp(b)d2b=∫[1−e−¯N(b)]d2b=2πr20Φ1(N0),
где
Φm(z)=∞∑k=1(−1)k+1zkk!km,Φ1(z)=∫z0(1−e−t)dtt.
После этого формула (31) принимает окончательный вид:
⟨N⟩=N0Φ1(N0),
т. е. параметр N0 однозначно определяет среднее число рассеченных померонов. Для приведенной дисперсии числа рассеченных померонов ωN согласно (32) и (35) это дает
ωN≡DN⟨N⟩=1+N02−⟨N⟩=1+N02−N0Φ1(N0).
То есть параметр N0, введенный в формуле (28) как коэффициент пропорциональности, однозначно определяет среднее число рассеченных померонов и дисперсию этого числа от события к событию. При этом зависимость от начальной энергии также входит только через этот параметр.
В рамках данной простой модели удается явно вычислить и саму вероятность pN иметь N рассеченых померонов в недифракционном pp-столкновении, усредняя ˜P(N,b)(24) по b при фиксированном N:
pN=∫˜P(N,b)f(b)d2b=1σpp∫P(N,b)d2b,
где мы приняли во внимание соотношение (27). Используя теперь (25), мы находим
где мы использовали формулы (33) и (34). Напомним, что величина σpp, определенная согласно (19), – это сечение недифракционного pp-взаимодействия, поэтому σN имеет смысл вклада в это сечение, возникающего от процесса с N рассечеными померонами.
Формула (43), полученная в рамках изложенного в этом разделе подхода, совпадает с хорошо известными результатами для σN, полученными в рамках квазиэйконального реджевского подхода [24], [25]:
σN=4πλCN[1−e−zN−1∑k=0zkk!],
где
z=2γCλeΔξ,λ=R2+α′ξ,ξ=lnss0.
Здесь s0≃1 ГэВ2, а 1+Δ и α′ – интерсепт и наклон померонной траектории. Параметры γ и R описывают вершину присоединения померона к рассеивающимся адронам. Квазиэйкональный параметр C эффективно учитывает вклад в эту вершину дифракционных процессов.
Совпадение формул (43) и (45) позволяет фиксировать параметры N0 и r0, используемые в нашем подходе, поскольку они однозначно выражаются через известные реджевские параметры:
N0=z=2γCλeΔξ,r0=√2λC,λ=R2+α′ξ.
В наших расчетах при моделировании распределения струн в плоскости прицельного параметра мы использовали следующие значения этих параметров:
которые позволяют описать поведение недифракционного сечения pp-рассеяния при энергиях БАК [26], [27].
5. Формирование струнных кластеров
Для эффективного описания процессов слияния струн и образования струнных кластеров [9]–[11] в поперечной к оси столкновения плоскости мы вводили конечную решетку (сетку) с площадью ячейки порядка поперечного сечения струны. Этот метод был предложен в работе [28] и затем успешно использовался при дальнейших исследованиях этого явления в работах [29]–[31].
В этом подходе предполагается, что происходит слияние в единый струнный кластер всех струн с центрами, попадающими в одну ячейку сетки. При этом в соответствии с моделью слияния струн [32] предполагается следующая зависимость средней множественности заряженных частиц, образующихся от распада кластера, образованного слиянием k струн, в быстротном интервале δy:
⟨μ⟩k=μ(k)0δy=μ0√kδy,
где μ(k)0 и μ(1)0≡μ0 – соответствующие множественности на единицу быстроты от распада кластера и одиночной струны.
Для расчета весовых коэффициентов α(k)(8), входящих в формулу (7), необходимо знать среднее число ⟨n⟩k частиц, рождающихся от всех кластеров, образованных слиянием ровно k струн. С учетом (49) его можно представить следующим образом:
⟨n⟩k=⟨mk⟩μ(k)0δy=⟨mk⟩√kμ0δy,
где ⟨mk⟩ – среднее число кластеров из k струн, образующихся в pp-взаимодействии при заданной энергии и центральности столкновения. Подставляя (50) в (8), находим для весовых коэффициентов
α(k)=⟨mk⟩√k∑∞k=1⟨mk⟩√k.
Из формулы (51) мы видим, что для расчета весовых коэффициентов α(k), входящих в формулу (7), необходимо знать только среднее число ⟨mk⟩ кластеров из k струн, образующихся в pp-взаимодействии при заданной энергии и центральности столкновения. К сожалению, эту часть вычислений пока не удается выполнить аналитически. Поэтому мы осуществили моделирование распределения струн в поперечной плоскости методом Монте-Карло (MK) в соответствии с алгоритмами, представленными в разделе 3. Мы также использовали аналитические результаты раздела 3 для контроля разработанных МК-кодов.
После генерации каждой струнной конфигурации, зная расположение всех центров первичных струн на сетке, мы подсчитали количество кластеров mk, образованных слиянием k струн. Усредняя по конфигурациям, мы нашли ⟨mk⟩, а затем по формуле (51) нашли весовой коэффициент α(k), входящий в формулу (7).
Для расчета сильноинтенсивных переменных Σk(μF,μB) кластеров, образованных слиянием k струн, также входящих в формулу (7) для Σ(nF,nB), мы использовали их выражения (10)–(12) через соответствующие корреляционные функции Λ(k)(Δy). Ранее в работе [22] из анализа разультатов коллаборации ALICE [33] по величине корреляций вперед-назад в pp-столкновениях на БАК [21] было показано, что для корреляционной функции одиночной струны Λ(1)(Δy) с хорошей точностью можно использовать простое экспоненциальное приближение, которое мы будем использовать и для струнных кластеров:
Λ(k)(Δy)=Λ(k)0e−|Δy|/y(k)corr,
где Δy=y1−y2, а y(k)corr – характерная длина корреляции по быстроте для частиц, рожденных данным струнным кластером.
Поскольку согласно (49) быстротная плотность заряженных частиц от распада струнного кластера считается пропорциональной √k, а корреляции имеют место только между соседними сегментами струнного кластера, разумно предположить, что характеристическая корреляционная длина y(k)corr убывает с ростом k как 1/√k. Некоторые дополнительные аргументы в пользу этого предположения были представлены в работе [22]. Поэтому мы предполагаем следующую зависимость параметров корреляционной функции струнного кластера от k:
y(k)corr=y(1)corr√k,Λ(k)0=const.
Отметим, что варианту без учета слияния струн отвечает совсем другой вид этих зависимостей:
˜μ(k)0=kμ(1)0,˜y(k)corr=const,˜Λ(k)0=Λ(1)0k.
Формулы (49) и (53) позволяют вычислить параметры струнного кластера, образованного слиянием k струн, исходя из трех параметров μ(1)0≡μ0, Λ(1)0, y(1)corr, характеризующих одиночную струну (52). После этого по формулам (10)–(12) можно найти Σk(μF,μB), входящие в формулу (7) для Σ(nF,nB). При этом в силу простой зависимости (52) интегралы (11) и (12) вычисляются явно [34].
6. Сравнение с данными эксперимента
Значение параметра μ0 (среднее число заряженных частиц на единицу быстроты от распада одиночной струны) необходимо выбрать равным 0.7 для правильного описания экспериментальной множественности заряженных частиц в центральной области быстрот.
Для фиксации двух оставшихся свободных параметров модели Λ(1)0 и y(1)corr, характеризующих корреляционную функцию одиночной струны (52), нами была рассчитана зависимость переменной Σ(nF,nB) от расстояния по быстроте Δy между окнами наблюдения малой ширины (δy=0.2) для pp-столкновений без отбора по центральности при двух начальных энергиях 0.9 и 7 ТэВ. В MK-реализации алгоритма, описанного в разделе 4, это отвечает случайной генерации вектора прицельного параметра b для каждого события.
Сравнение результатов этих расчетов с предварительными экспериментальными данными коллаборации ALICE для Σ(nF,nB) [19], полученными при анализе рождения заряженных частиц в мягкой области спектра с поперечными импульсами в диапазоне 0.3–1.5 ГэВ/c (см. две нижние кривые на рис. 1), дает следующие значения этих двух параметров:
y(1)corr=2.7,Λ(1)0=0.8.
После этого, уже с полностью фиксированными параметрами модели, мы провели аналогичные вычисления переменной Σ(nF,nB) для случая pp-столкновений без отбора по центральности при начальной энергии 13 ТэВ. Сравнение полученных результатов с экспериментальными данными ALICE [19], полученными при 13 ТэВ (верхняя кривая на рис. 1), показывает, что этот набор параметров также хорошо описывает зависимость наблюдаемой Σ(nF,nB) от быстротного расстояния между центрами окон наблюдения Δy и при начальной энергии 13 ТэВ.
При этом мы учли, что поскольку анализ этих данных при 13 ТэВ в эксперименте ALICE проводился для более широкого диапазона поперечных импульсов 0.2–2 ГэВ/c, который в настоящее время используется в коллаборации ALICE при анализе рождения заряженных частиц в мягкой области спектров, то из-за увеличения общего числа частиц, регистрируемых в этом более широком интервале поперечных импульсов, коэффициент μ0 должен быть увеличен согласно анализу работы [35] в 1.28 раза до значения μ0=0.9.
Из рис. 1 мы видим, что характер зависимости переменной Σ(nF,nB) от расстояния по быстроте между окнами наблюдения Δy в целом соответствует предсказываемому нашей моделью (см. конец раздела 3). Наблюдается рост величины Σ(nF,nB) от значений порядка единицы при малых значениях Δy с некоторой тенденцией к насыщению при больших значениях Δy.
Что касается зависимости от начальной энергии, то сравнивая зависимости при энергиях 0.9, 7 и 13 ТэВ, представленные на рис. 1, мы видим, что в pp-столкновениях при энергиях БАК наблюдается рост наблюдаемой Σ(nF,nB) с увеличением начальной энергии. В нашей модели этот рост объясняется возрастающим с энергией вкладом от образования струнных кластеров с новыми свойствами, состоящих из все большего числа слившихся струн. Отметим, что в модели с одним типом струн, без учета процессов их слияния и образования струнных кластеров, все три кривые на рис. 1 для разных начальных энергий оказались бы совпающими друг с другом.
7. Заключение
В рамках модели, учитывающей слияние кварк-глюонных струн и образование струнных кластеров, для случая pp-рассеяния при высоких энергиях рассчитана сильноинтенсивная переменная Σ(nF,nB), характеризующая корреляции между числами частиц nF и nB, образующихся в двух разнесенных по быстроте интервалах наблюдения. Использование этой переменной для изучения корреляций позволяет минимизировать вклад тривиальных “объемных” флуктуаций и дает возможность получить информацию о свойствах источников (струн и струнных кластеров), образующихся на начальной стадии сильного взаимодействия.
В частности, сравнение полученных теоретических зависимостей для переменной Σ(nF,nB) с предварительными экспериментальными данными [19] коллаборации ALICE на БАК при трех начальных энергиях рр-столкновений позволяет сделать однозначный вывод о наличии эффектов слияния струн и образования струнных кластеров при энергиях БАК уже в рр-взаимодействиях.
При этом, как и ожидалось, вклад струнных кластеров в переменную Σ(nF,nB) растет с ростом начальной энергии. В работе [36] показано, что тот же эффект имеет место и при отборе событий с более центральными pp-столкновениями. Сравнение полученных теоретических зависимостей с экспериментальными данными позволяет также фиксировать параметры, характеризующие кластеры с различным числом слившихся струн.
Благодарности
Авторы благодарны организаторам VII международной конференции “Модели квантовой теории поля” (MQFT-2022).
Конфликт интересов
Авторы заявляют, что у них нет конфликта интересов.
Список литературы
1.
A. Dumitru, F. Gelis, L. McLerran, R. Venugopalan, “Glasma flux tubes and the near side ridge phenomenon at RHIC”, Nucl. Phys. A, 810:1–4 (2008), 91–108
2.
A. B. Kaidalov, “The quark-gluon structure of the pomeron and the rise of inclusive spectra at high energies”, Phys. Lett. B, 116:6 (1982), 459–463
3.
A. Capella, U. Sukhatme, C.-I. Tan, J. T. Thanh Van, “Dual parton model”, Phys. Rep., 236:4–5 (1994), 225–329
4.
A. Casher, J. Kogut, L. Susskind, “Vacuum polarization and the absence of free quarks”, Phys. Rev. D, 10:2 (1974), 732–745
5.
A. Casher, H. Neuberger, S. Nussinov, “Chromoelectric-flux-tube model of particle production”, Phys. Rev. D, 20:1 (1979), 179–188
6.
M. Gyulassy, A. Iwazaki, “Quark and gluon pair production in SU(N) covariant constant fields”, Phys. Lett. B, 165:1–3 (1985), 157–161
7.
F. Bissey, A. I. Signal, D. B. Leinweber, “Comparison of gluon flux-tube distributions for quark-diquark and quark-antiquark hadrons”, Phys. Rev. D, 80:11 (2009), 114506, 6 pp.
8.
P. Cea, L. Cosmai, F. Cuteri, A. Papa, “Flux tubes in the QCD vacuum”, Phys. Rev. D, 95:11 (2017), 114511, 9 pp.
9.
T. S. Biro, H. B. Nielsen, J. Knoll, “Colour rope model for extreme relativistic heavy ion collisions”, Nucl. Phys. B, 245 (1984), 449–468
10.
A. Bialas, W. Czyz, “Chromoelectric flux tubes and the transverse-momentum distribution in high-energy nucleus-nucleus collisions”, Phys. Rev. D, 31:1 (1985), 198–200
11.
M. A. Braun, C. Pajares, “A probabilistic model of interacting strings”, Nucl. Phys. B, 390:2 (1993), 542–558
12.
C. Bierlich, G. Gustafson, L. Lönnblad, A. Tarasov, “Effects of overlapping strings in pp collisions”, JHEP, 2015:3 (2015), 148, 49 pp.
13.
J. Adam, D. Adamova, M. M. Aggarwal et al. [ALICE Collab.], “Enhanced production of multi-strange hadrons in high-multiplicity proton-proton collisions”, Nature Phys., 13 (2017), 535–539
14.
V. N. Kovalenko, A. M. Puchkov, G. A. Feofilov, “Production of strange particles in a multi-pomeron exchange model”, Bull. Russ. Acad. Sci. Phys., 80:8 (2016), 966–969
15.
N. Armesto, D. A. Derkach, G. A. Feofilov, “pt-multiplicity correlations in a multi-Pomeron-exchange model with string collective effects”, Phys. Atom. Nucl., 71:12 (2008), 2087–2095
16.
V. N. Kovalenko, G. A. Feofilov, A. M. Puchkov, F. Valiev, “Multipomeron model with collective effects for high-energy hadron collisions”, Universe, 8:4 (2022), 246, 25 pp.
17.
M. I. Gorenstein, M. Gaździcki, “Strongly intensive quantities”, Phys. Rev. C, 84:1 (2011), 014904, 5 pp.
18.
Е. В. Андронов, “Влияние механизма слияния кварк-глюонных струн на дальние быстротные корреляции и флуктуации”, ТМФ, 185:1 (2015), 28–36
19.
A. Erokhin [ALICE Collab.], “Forward-backward multiplicity correlations with strongly intensive observables in pp collisions”, The VIth International Conference on the Initial Stages of High-Energy Nuclear Collisions (IS2021) (Weizmann Institute of Science, 10–15 January, 2021), Saint Petersburg State University, Laboratory of Ultra-High Energy Physics, St. Petersburg, 2021, Poster/187https://indico.cern.ch/event/854124/contributions/4134683/
20.
С. Н. Белокурова, В. В. Вечернин, “Сильноинтенсивные переменные и дальние корреляции в модели с решеткой в поперечной плоскости”, ТМФ, 200:2 (2019), 195–214
21.
V. Vechernin, “Forward-backward correlations between multiplicities in windows separated in azimuth and rapidity”, Nucl. Phys. A, 939 (2015), 21–45
22.
V. Vechernin, “Short- and long-range rapidity correlations in the model with a lattice in transverse plane”, EPJ Web Conf., 191 (2018), 04011, 8 pp.
23.
V. Vechernin, I. Lakomov, “The dependence of the number of pomerons on the impact parameter and the long-range rapidity correlations in pp collisions”, PoS (Baldin ISHEPP XXI), 2012, 072, 12 pp.
24.
K. A. Ter-Martirosyan, “On the particle multiplicity distributions at high energy”, Phys. Lett. B, 44B:4 (1973), 377–380
25.
A. A. Кайдалов, E. A. Тер-Мартиросян, “Множественное рождение адронов при высоких энергиях в модели кварк-глюонных струн. Сравнение с экспериментом”, ЯФ, 40:1(7) (1984), 211–220
26.
A. Capella, E. G. Ferreiro, “Charged multiplicities in pp and AA collisions at LHC”, Eur. Phys. J. C, 72 (2012), 1936, 6 pp.
27.
J. Bleibel, L. V. Bravina, E. E. Zabrodin, “How many of the scaling trends in pp collisions will be violated at √sNN=14 TeV? Predictions from Monte Carlo quark-gluon string model”, Phys. Rev. D, 93:11 (2016), 114012, 13 pp.
28.
В. В. Вечернин, Р. С. Колеватов, “Дискретный подход к описанию дальных корреляций множественности и pt в модели слияния струн”, Вестн. С.-Петерб. ун-та. Сер. 4. Физика. Химия, 4 (2004), 11–27, arXiv: hep-ph/0305136
29.
M. A. Braun, R. S. Kolevatov, C. Pajares, V. V. Vechernin, “Correlations between multiplicities and average transverse momentum in the percolating color strings approach”, Eur. Phys. J. C, 32:4 (2004), 535–546
30.
M. A. Braun, C. Pajares, V. V. Vechernin, “Anisotropic flows from colour strings: Monte Carlo simulations”, Nucl. Phys. A, 906 (2013), 14–27
31.
M. A. Braun, C. Pajares, V. V. Vechernin, “Ridge from strings”, Eur. Phys. J. A, 51 (2015), 44, 11 pp., arXiv: 1407.4590
32.
M. A. Braun, C. Pajares, “Inplication of percolation of colour strings on multiplicities, correlations and the transverse momentum”, Eur. Phys. J. A, 16:2 (2000), 349–359, arXiv: hep-ph9907332
33.
J. Adam, D. Adamova, M. M. Aggarwal et al. [ALICE Collab.], “Forward-backward multiplicity correlations in pp collisions at √s=0.9, 2.76 and 7 TeV”, JHEP, 2015 (2015), 097, 28 pp., arXiv: 1502.00230
34.
S. Belokurova, “Study of strongly intense quantities and robust variances in multi-particle production at LHC energies”, Phys. Part. Nucl., 53:2 (2022), 154–158
35.
E. Andronov, V. Vechernin, “Strongly intensive observable between multiplicities in two acceptance windows in a string model”, Eur. Phys. J. A, 55 (2019), 14, 12 pp., arXiv: 1808.09770
36.
S. Belokurova, V. Vechernin, “Using a strongly intense observable to study the formation of quark-gluon string clusters in pp collisions at LHC energies”, Symmetry, 14:8 (2022), 1673, 11 pp.
Образец цитирования:
В. В. Вечернин, С. Н. Белокурова, “Сильноинтенсивная переменная в модели высокоэнергетических
pp-взаимодействий с образованием струнных кластеров”, ТМФ, 216:3 (2023), 460–475; Theoret. and Math. Phys., 216:3 (2023), 1299–1312
\RBibitem{VecBel23}
\by В.~В.~Вечернин, С.~Н.~Белокурова
\paper Сильноинтенсивная переменная в модели высокоэнергетических
pp-взаимодействий с образованием струнных кластеров
\jour ТМФ
\yr 2023
\vol 216
\issue 3
\pages 460--475
\mathnet{http://mi.mathnet.ru/tmf10461}
\crossref{https://doi.org/10.4213/tmf10461}
\mathscinet{http://mathscinet.ams.org/mathscinet-getitem?mr=4634826}
\adsnasa{https://adsabs.harvard.edu/cgi-bin/bib_query?2023TMP...216.1299V}
\transl
\jour Theoret. and Math. Phys.
\yr 2023
\vol 216
\issue 3
\pages 1299--1312
\crossref{https://doi.org/10.1134/S0040577923090052}
\scopus{https://www.scopus.com/record/display.url?origin=inward&eid=2-s2.0-85172475862}
Образцы ссылок на эту страницу:
https://www.mathnet.ru/rus/tmf10461
https://doi.org/10.4213/tmf10461
https://www.mathnet.ru/rus/tmf/v216/i3/p460
Эта публикация цитируется в следующих 1 статьяx:
Svetlana Belokurova, Vladimir Vechernin, “Multiplicity Fluctuations in the String Model in pp Collisions at LHC Energies”, Phys. Atom. Nuclei, 87:S3 (2024), S471