CC BY
54
УДК 539.1
РОЖДЕНИЕ БОЗОНОВ ХИГГСА В РЕДЖЕВСКОМ ПРЕДЕЛЕ КВАНТОВОЙ ХРОМОДИНАМИКИ НА БАК
М. А. Нефедов, В. А. Салеев
Самарский государственный университет,
Россия, 443011, Самара, ул. Академика Павлова, 1.
E-mails: nef edovma@gmail. com, saleev@samsu. ru
В теории реджезованных партонов Л. Н. Липатова рассматривается рождение скалярного бозона Хиггса Н Стандартной модели и псевдоскалярного бозона Хиггса А минимальной суперсимметричной Стандартной модели в процессах слияния реджезованных глюонов (RR —s- Н, А) при энергиях Большого адронного коллайдера. Расчеты выполнены в лидирующем приближении по константе сильного взаимодействия as с учетом вкладов от t- и b-кварков в петлевые интегралы. Рассчитаны спектры по поперечному импульсу и зависимость полных сечений рождения от масс бозонов. Полученные результаты согласуются с предсказаниями, полученными ранее в коллинеарной партонной модели с учетом следующих поправок по as и эффектов суммирования вкладов, усиленных большими логарифмами, во всех порядках теории возмущений.
Ключевые слова: бозон Хиггса, Стандартная модель, минимальная суперсим-метричная Стандартная модель, квантовая хромодинамика, теория реджезованных партонов.
Введение. Исследование процессов рождения бозонов Хиггса в Стандартной модели (СМ) и минимальной суперсимметричной Стандартной модели (МССМ) на Большом адронном коллайдере (БАК) рассматривается как центральная задача физики высоких энергий текущего десятилетия. Обнаружение на БАК стандартного скалярного бозона Хиггса будет блестящим подтверждением СМ и очередным триумфом квантовой калибровочной теории поля. Другая важная проблема, которую ожидается решить в экспериментах на БАК, это получение ответа на вопрос о месте суперсимметричных обобщений СМ. Первые указания на существование скалярной частицы с массой около 125 ГэВ были получены в 2012 году Коллаборациями ATLAS [1] и CMS [2]. Превышение сигнала над фоном составляет около 4 стандартных отклонений, что позволяет надеяться на достоверность полученных результатов. В этой связи ещё более актуальным становится моделирование процессов рождения бозонов Хиггса на БАК в теории реджезованных партонов (ТРП), которая позволяет учесть нетривиальные эффекты, связанные с особенностями квантовой хромодинамики (КХД) при высоких и сверхвысоких энергиях.
В СМ основной вклад в сечения рождения хиггсовского бозона с массой порядка 100 4- 500 ГэВ в протон-протонных взаимодейсвиях при энергии протонов в системе центра масс на БАК y/s = 7 4- 14 ТэВ дает процесс глюон-глюонного слияния g + g —> Н, идущий через кварковую петлю. Наибольший вклад в этот процесс дает петля t-кварка, т. к. вершина взаимодействия
Максим Александрович Нефедов, аспирант, каф. общей и теоретической физики. Владимир Анатольевич Салеев (д.ф.-м.н., проф.), профессор, каф. математического моделирования в механике.
бозона Хиггса СМ пропорциональна массе кварка в петле:
.т 1
&н = г—, V =
где Ср — константа Ферми. Величина V ~ 246 ГэВ характеризует энергетический масштаб спонтанного нарушения электрослабой симметрии. Следующий по величине вклад дает процесс рождения глюонами хиггсовского бозона совместно с парой «кварк-антикварк», в основном с парой 6-кварков (ЪЪН), затем идут процессы ассоциативного рождения бозона Хиггса с тяжелыми калибровочными бозонами (\¥Н, ZH) и парой £-кварков (ЫН).
В связи с экспериментальными условиями, предполагающими, что на начальном этапе поиска бозона Хиггса на БАК статистика будет недостаточной для измерения дифференциальных распределений, имеющиеся расчеты в коллинеарной партонной модели (КПМ) в основном ориентированы на определение полного сечения рождения бозона Хиггса в СМ, которое в настоящее время известно как в следующем за лидирующим приближении по а3 (СЛП) с учётом конечной массы топ кварка [5], так и в следующем за СЛП КПМ в пределе т* —> оо [6]. Для расчета распределений по поперечному импульсу (рт) используется широкий класс генераторов Монте—Карло [7] в сочетании с процедурой суммирования больших логарифмических поправок в КПМ [7-9]. Сочетание этих методов позволяет получать стабильные предсказания для ^-распределений, которых достаточно для целей феноменологии СМ. Однако в феноменологии МССМ мы вынуждены иметь дело с большим числом свободных параметров и делать предсказания при их варьировании в широкой области. Это приводит к существенным неопределенностям при использовании упомянутых выше подходов.
В МССМ спонтанное нарушение электрослабой симметрии обеспечивается введением двух дублетов хиггсовских полей, что приводит к возникновению двух скалярных (Н и /г), одного псевдоскалярного (А) и двух заряженных (Н±) бозонов Хиггса. Генерация масс верхних и нижних фермионов в Би(2)-дублете обеспечивается разными дублетами хиггсовских полей с различными вакуумными средними для этих дублетов 1)( и В результате естественным параметром модели является отношение этих вакуумных средних:
Вершины взаимодействия псевдоскалярного бозона Хиггса МССМ с ферми-онами имеют вид
V V
Экспериментальные данные [4] дают следующие ограничения на свободные параметры МССМ: 10 < tg(/3) < 30 и Ма > 100 ГэВ; при таких значениях параметров вклады петель £- и 6-кварков в амплитуду рождения бозона Хиггса А будут иметь одинаковый порядок величины. С другой стороны, неколлинеарные поправки к петле 6-кварка могут быть значительны, т. к. поперечные импульсы реджезованных глюонов имеют величину порядка массы 6-кварка, но много меньше массы £-кварка. В КПМ последнее обстоятельство
приводит к необходимости непростого учета старших поправок теории возмущений, в то время как в ТРП основная часть этих поправок учитывается уже в лидирующем приближении.
1. Амплитуды рождения скалярного и псевдоскалярного бозонов Хиггса
в ТРП. При высоких энергиях основной вклад в инклюзивные сечения рождения частиц дают процессы в мультиреджевской кинематике. В этих условиях начальные £-канальные партоны реджезуются, что является фундаментальным свойством всех калибровочных теорий поля. В работах Л. Н. Липатова сформулирована ТРП как обобщение КХД [11]. Лагранжиан ТРП содержит поля реджезованных глюонов и кварков, а их взаимодействие с Янг—Милл-совскими глюонами и кварками ПМ осуществляется при помощи индуцированных вершин [12]. В рассматриваемом процессе важна вершина перехода реджезованного глюона в Янг—Миллсовский К —>■ д:
где ??+" = 2Р»//Б, п~» = 2Р£/у/Б, (п+)2 = (и")2 = 0, (п+п~) = 2, 2(РіР2) =
где <71,2 = Жі,2-Рі,2 + <7і,2т — 4-импульсы реджезованных партонов, дт = = (0^т,0), = q2 2т- В ТРП амплитуда рождения скалярного бозона
Хиггса в слиянии двух реджезованых глюонов (упускаем константу связи и цветовую структуру) имеет вид
где т — масса кварка, а Мн — масса бозона Хиггса. Скалярные произведения, через которые выражается амплитуда, записываются в таком виде:
(<7і + <72 )2 — Мн, 2 (<71(72) — Мн + + і‘2, хіх'ів — Мн + р т — М^,
где рт = qlт + q2T• В рамках метода размерной регуляризации и процедуры Пассарино—Вельтмана [13] данная амплитуда может быть сведена к линейной комбинации базисных однопетлевых интегралов Во и Со с коэффициентами, зависящими от 1\, £2- Выражения для этих коэффициентов являются
5, Р\,2 — 4-импульсы сталкивающхся протонов.
Диаграммы ТРП в лидирующем приближении для процесса ПК —> Н приведены на рис. 1.
Необходимо отметить, что для согласования правил Фейнмана в ТРП с правилами Фейнмана в стандартной КХД амплитуды в ТРП домножается на нормировочный множитель
а
б
Рис. 1
АА = ,
4^1*2
ХіХ-ів
X
{/
4 Тг [7»(д + т)У(д - <12 + т)(д + д\ + т)} ^ (д2 — т2)((д — д2)2 — т2)((д + д\)2 — т2)
+ (91 ** 92, ^ ** V) | ,
(1)
весьма громоздкими, но они могут быть упрощены, если ввести новые функции:
D±(ti,t2) = Mfj ± + t2) + — t2)2.
Тогда выражение для амплитуды принимает вид As = ~tb m2; m2>m2)+
+ 8В ■ Во(Мн,т2,m2) - 4C(ti,t2) ■ B0(-ti,m2,m2)-
- 4C(t2, t{) ■ B0(-t2, m2, m2) + 32D}, (2)
где коэффициенты А, В, С, D имеют следующий вид:
А = Му (£>_ -D+- mfj) х х (SMfj (3D_ + bD+) - 2D_D+ + 3D2_ + 64m2D+Mf1 - D\ + 48M|) +
+ 4£>+М|(-8М^ (L>_ + 3D+) + 2D_D+ + D2_ +
+ QAm2D+M2H - 3D\ + 16M|), В = Mjj (M| (24Mfj (D_ + L>+) + 2£>_£>+ - 3D2 + - 48M|) -
-AD+M2h (D_ + 3D+ - 4M^)) , C(t\,t2) = {D- — D+ + 4M^(t2 — £i)) (3M2 (_D_ — 4M#) (M^ — (£2 — ti)) + +D+ (M| (5M| - 4(t2 - ti)) + ШАН + 3M|(t2 - ti))) ,
= М|£>+ (М|(£>_ - 4M^) + D+(4Mfj - M|)) .
Амплитуда (2) не имеет ультрафиолетовых расходимостей, сингулярные члены входящие в определение функции До, сокращаются. На стадии чис-
ленных расчетов мы используем представление функции Со через комплексные дилогарифмы Ъ\2(г), полученное в работе [14]. Квадрат модуля амплитуды рождения СМ Хиггса в ЛП ТРП с учетом вклада 6-кварка и опущенных ранее цветовых множителей и констант связи имеет вид
______ Tj-2 2
|Asm|2 = -j=GFa2s\mtAs(m = mt) + mbAs{m = mt)\ .
После перехода к пределу t\t2 —> 0 в (3) и усреднения по углу между поперечными импульсами партонов квадрат модуля амплитуды в ТРП переходит в соответствующий квадрат модуля амплитуды в КПМ [15]:
И1 + (1-тжп)|2,
Рассмотрим рождение псевдоскалярного бозона Хиггса МССМ. Амплитуда записывается аналогично (1):
л ___ ( 1) Х\Х2Р1цР2и дци / , , д^2 \
Р3 ~ (2тт)4 у/Г1л/Ь ЛР^~1ъ ~*2’ А’ Ш)’
/ _/ м2 !' ^^№(9 + т)У{д - д2 + т)7б(д + 91 + т)] ,
1’ ] а д (<?2_т2)(((?_(?2)2_т2)(((? + (?1)2_т2) +
+(<71 ^ <72,/л о г/).
Можно показать, что квадрат модуля амплитуды рождения бозона Хиггса А в ЛП ТРП приводится к виду
______ 2 а2М4
\Аа\2 = ^2(Ф) \™ъСо(М%, —*1, —*2, гп2ъ, т2ъ) tg(/?)+
+т^С0(м1, -и, ^2) т?, т?, т42) ctg(/5)|2 , где ф — угол между векторами qlг и q2т• В коллинеарном пределе получаем
2 дж4
а%М'
\ЛС/М\2 = Ш(п) tg(/5) + г*/(г4) й£(/3)|2,
где = Ат2ь/Мц, что совпадает с известным выражением в КПМ [16].
2. Сечения рождения бозона Хиггса СМ на БАК. Дифференциальное по поперечному импульсу (рт) и быстроте (у) сечение рождения бозона Хиггса в слиянии реджезованных глюонов согласно теореме факторизации в ТРП записывается в виде [17]:
ГСО Г 27Г
ск/йрт = М1.1 о ,10 >
где
XI,2 = ~/^е±У’ Мт = \!МН +РТ> 12=Рт + и-‘}ртл/Г1СОи(ф1)]
ф\— угол между ql7’ и р^, а Ф(х, £, (л2) — глюонная неинтегрированная пар-тонная функция распределения (НПФР). В наших расчетах мы будем использовать глюонную НПФР, полученную в работах Кимбера, Мартина и Рыски-на [19], которая соответствует коллинеарной ПФР МИ8Т2002 [18]. В ТРП мы учитываем лидирующие логарифмические поправки вида 1п(///Л<дс£>) и 1п(1/ж), в то же время из расчётов в СЛП КПМ известно, что нелогарифмические поправки также дают большой вклад в сечение. Эти вклады обусловлены петлевыми глюонными поправками к вершине рождения Хиггса и могут быть эффективно учтены при помощи умножения сечения на А'-фактор:
Г7Г
К = ехр у-САа3{11цЦ ,
где С а = Л'с = 3. Данный Л'-фактор получается из аналитического продолжения формфактора Судакова во временпподобную область переданных импульсов (см., например, [20] и имеющиеся там ссылки).
При вычислении сечений мы используем следующий выбор масштабов перенормировки и факторизации: /хд = цр = £Мт, где для центральной кривой на графиках соответствует прямая £ = 1, а коридор теоретической неопределённое™ вносимой выбором масштаба, строится путём варьирования 1/2 < £ < 2.
На рис. 2 приведено сравнение наших предсказаний (сплошная гистограмма и сплошная заливка) для ру-спектра бозонов Хиггса при л/е = 7 ТэВ и Мн = 120 ГэВ с учетом вклада 6-кварка с результатами СЛП КПМ методом Монте—Карло, взятыми из работы [7]. Полоса неопределенности с вертикальной штриховкой соответствует результату генератора Монте—Карло РуЬМа с программным кодом РСМНЕС [10]. Полоса неопределенности с горизонтальной штриховкой показывает результат расчета в СЛП КПМ в программе НдТ [9].
Рис. 2. показывает, что расчет в ЛП ТРП удовлетворительно воспроизводит результаты СЛП КПМ, особенно если принять во внимание степень неопределенности последних. Результат ЛП ТРП менее устойчив по отношению к варьированию масштабов перенормировки и факторизации, что, впрочем, характерно для расчетов в лидирующем приближении.
Предсказания ЛП ТРП и СЛП КПМ для зависимости полного сечения рождения бозона Хиггса СМ от его массы в наиболее интересной для эксперимента области 120 < Мн < 130 ГэВ отличаются не более чем на 25 %, что является хорошим аргументом в пользу применения расчетов в ЛП ТРП для оценки различных дифференциальных сечений. В частности мы можем сделать предсказание для ру-спектра бозона Хиггса с массой 126 ГэВ при у/е = 14 ТэВ, который может быть измерен после модернизации БАК до этой энергии (рис. 3).
(кг/брг, пб/ГэВ йст/фг, пб/ГэВ
10° С---------------------------------------------------------- г---------------------
Рис. 2 Рис. 3
3. Сечения рождения псевдоскалярного бозона Хиггса МССМ. На рис. 4 приведено сравнение наших предсказаний для полного сечения рождения бозона Хиггса МССМ А в ЛП ТРП при 1^(/3) = 3 (кривая 1) и 1^(/3) = 30
cr(RR —> А), пб
МА, ГэВ
(кривая 3) с соответствующими результатами, полученными в СЛП КПМ [16] (кривые 2 и 4). Видно, что при tg(/3) = 3 ЛП ТРП хорошо воспроизводит результат СЛП КПМ практически для всех значений массы бозона, в то время как для t.g(/3) = 30 при больших значениях массы наш подход предсказывает несколько меньшее сечение, чем СЛП КПМ. Однако кривая КПМ лежит в пределах полосы теоретической неопределенности наших расчетов, что позволяет надеяться на разумное согласие предсказаний и для спектров по поперечному импульсу с результатами КПМ.
Рис. 4 Результаты расчетов дифференци-
ального эффективного сечения рождения псевдоскалярного бозона Хиггса МССМ в процессе слияния редже-зованных глюонов в ЛП ТРП для л/s = 14 ТэВ при значениях Ма = 126, 200, 500, 1000 ГэВ приведены на рис. 5 (а, б, в, г, соответственно). Сплошная кривая соответствует значению t.g(/3) = 3, а штриховая кривая — t.g(/3) = 30. Из рисунков видно, что спектры при большом и малом значениях t.g(/3) имеют разный наклон. В случае измерения данных спектров на эксперименте этот факт можно использовать для более точного определения значения этого параметра.
da/dpr, пб/ГэВ
da/dpr, пб/ГэВ
da/dpr, пб/ГэВ
рт, ГэВ
Рис. 5
Также можно заключить, что при увеличении Ма сечение, соответствующее tg(/5) = 30, становится меньше, чем сечение для tg(/5) = 3. Это связано с влиянием интерференции t и 6-кварковых вкладов, которая становится деструктивной при Ма TTlt b в силу свойств функции Со.
Заключение. Эксперименты на БАК в тэвной области энергий показывают, что доминирующую роль в процессах рождения частиц в центральной области по быстроте играют подпроцессы в мультиреджевской кинематике, для описания которых необходимо использовать ТРП вместо КПМ. Результаты этой работы, а также наших расчетов спектров чармониев [21], 6-кварковых струй [22], адронных струй и прямых фотонов [23] подтверждают это наблюдение.
Работа выполнена при поддержке РФФИ (проекты №№ 11-02-00769-а, 12-02-31701-мол-а).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. ATLAS Collaboration, С. Aad, et. al, “Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC” // Phys. Lett. B., 2012. Vol. 716, no. 1. Pp. 1-29, arXiv: 1207.7214 [hep-ex].
2. CMS Collaboration, S. Chatrchyan, et. al., “Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC” // Phys. Lett. B., 2012. Vol. 716, no. 1. Pp. 30-61, arXiv: 1207.7235 [hep-ex].
3. J. Baglio, A. Djouadi, R. M. Codbole, “The apparent excess in the Higgs to di-photon rate at the LHC: New Physics or QCD uncertainties?” // Phys. Lett. B, 2012. Vol. 716, no. 1. Pp. 203-207, arXiv: 1207.1451 [hep-ph],
4. ATLAS Collaboration, C. Aad, et. al., “Search for neutral MSSM Higgs bosons decaying to t+t~ pairs in proton-proton collisions at y/s = 7 TeV with the ATLAS detector” // Phys. Lett. B, 2011. Vol. 705, no. 3. Pp. 174-192, arXiv: 1107.5003 [hep-ex],
5. M. Spira, A. Djouadi, D. Craudenz, P. M. Zerwas, “Higgs boson production at the LHC” // Nucl. Phys. B., 1995. Vol. 453, no. 1-2. Pp. 17-82, arXiv: hep-ph/9504378.
6. C. Anastasiou, K. Melnikov, “Higgs boson production at hadron colliders in NNLO QCD” // Nucl. Phys. B., 2002. no. 1-2. Pp. 220-256, arXiv: hep-ph/0207004.
7. LHC Higgs Cross Section Working Group, S. Dittmaier, et. al., Handbook of LHC Higgs Cross Sections: 2. Differential Distributions: to be submitted to CERN Report, 2012. 275 pp., arXiv: hep-ph/1201.3084.
8. Yu. L. Dokshitzer, D. I. Dyakonov, S. I. Troyan, “Hard processes in quantum chromodynamics” // Phys. Rep., 1980. Vol. 58, no. 5. Pp. 269-395.
9. C. Bozzi, S. Catani, D. de Florian, M. Crazzini, “The qr spectrum of the Higgs boson at the LHC in QCD perturbation theory” // Phys. Lett. B, 2003. Vol. 564, no. 1-2. Pp. 65-72, arXiv: hep-ph/0302104.
10. S. Frixione, P. Nason, C. Oleari, “Matching NLO QCD computations with parton shower simulations: the POWHEG method” // JHEP, 2007. Vol. 2007, no. 11, 070.
11. L. N. Lipatov, “Gauge invariant effective action for high energy processes in QCD” // Nucl. Phys. B., 1995. Vol. 452, no. 1-2. Pp. 369-397, arXiv: hep-ph/9502308.
12. E. N. Antonov, L. N. Lipatov, E. A. Kuraev, I. O. Cherednikov, “Feynman rules for effective Regge action”// Nucl. Phys. B., 2005. Vol. 721, no. 1-3. Pp. 111-135, arXiv: hep-ph/0411185.
13. C. Passarino, M. Veltman, “One-loop corrections for e+e~ annihilation into уГ in the Weinberg model” // Nucl. Phys. B., 1979. Vol. 160, no. 1. Pp. 151-207.
14. О. Л. Веретин, О. В. Теряев, Аксиальная аномалия при произвольных виртуальностях: Препринт ОИЯИ № Р2-94-483. Дубна: ОИЯИ, 1994. 12 с. [О. L. Veretin, О. V. Teryaev, Axial anomaly at the arbitrary external momenta: JINR Preprint P2-94-483. Dubna: JINR, 1994. 12 pp.]
15. A. Djouadi, “The anatomy of electroweak symmetry breaking: Tome I: The Higgs boson in the Standard Model”// Phys. Rept., 2008. Vol. 457, no. 1-4. Pp. 1-216, arXiv: hep-ph/0503172.
16. A. Djouadi, “The anatomy of electroweak symmetry breaking Tome II: The Higgs bosons in the Minimal Supersymmetric Model” // Phys. Rept., 2008. Vol. 459, no. 1-6. Pp. 1-241, arXiv: hep-ph/0503173.
17. B. L. Ioffe, V. S. Fadin, L. N. Lipatov, Quantum Chromodynamics, Perturbative and Nonperturbative Aspects. Cambridge, UK: Cambridge University Press, 2010. 596 pp.
18. A. D. Martin, R. G. Roberts, W. J. Stirling, R. S. Thorne, “NNLO global parton analysis” // Phys. Lett. B, 2002. Vol. 531, no. 3-4. Pp. 216-224, arXiv: hep-ph/0201127.
19. M. A. Kimber, A. D. Martin, M. G. Ryskin, “Unintegrated parton distributions” // Phys. Rev. D, 2001. Vol. 63, no. 11, 114027. 10 pp., arXiv: hep-ph/0101348.
20. G. Watt, A. D. Martin, M. G. Ryskin, “Unintegrated parton distributions and electroweak boson production at hadron colliders” // Phys. Rev. D, 2004. Vol. 70, no. 1, 014012. 9 pp., arXiv: hep-ph/0309096.
21. V. A. Saleev, M. A. Nefedov, A. V. Shipilova, “Prompt J/ф production in the Regge limit of QCD: From the Tevatron to the LHC” // Phys. Rev. D, 2012. Vol. 85, no. 7, 074013. 9 pp., arXiv: 1201.3464 [hep-ph].
22. V. A. Saleev, A. V. Shipilova, “Inclusive 6-jet and 66-dijet production at the LHC via Reggeized gluons”// Phys. Rev. D, 2012. Vol. 86, no. 3, 034032. 9 pp., arXiv: 1201.4640 [hep-ph].
23. B. A. Kniehl, V. A. Saleev, A. V. Shipilova, E. V. Yatsenko, “Single jet and prompt-photon inclusive production with multi-Regge kinematics: From Tevatron to LHC” // Phys. Rev. D, 2011. Vol. 84, no. 7, 074017. 8 pp.
Поступила в редакцию 15/XI/2012; в окончательном варианте — 27/III/2013.
MSC: 81V05, 81Т60
HIGGS BOSON PRODUCTION AT THE LHC IN THE FRAMEWORK OF REGGE LIMIT OF QUANTUM CHROMODYNAMICS
M. A. Nefedov, V. A. Saleev
Samara State University,
1, Academician Pavlov St., Samara, 443011, Russia.
E-mails: nef edovma@gmail. com, saleev@samsu. ru
In the framework of the L. N. Lipatov Reggeized Partons Theory we consider production of scalar Higgs boson H of the Standard Model and pseudoscalar Higgs boson A of the Minimal Supersymmetric Standard Model in the Reggeized gluon fusion at the energy range of the Large Hadron Collider. Working in the leading order approximation in the strong coupling constant as we take into account contributions from, t- and b-quark in the loop integrals. We have calculated transverse momentum boson spectra and total cross sections as a function of boson mass. The obtained results agree with predictions which have been obtained, earlier in the collinear Parton Model with, the next to leading order in as corrections and, all-order resummation of large logarithmi-cal corrections.
Key words: Higgs boson, Standard Model, Minimal Supersymmetric Standard Model, Quantum, Chromodynamics, Theory of Reggeized, Partons.
Original article submitted 15/XI/2012; revision submitted 27/111/2013.
Maxim A. Nefedov, Postgraduate Student, Dept, of General and Theoretical Physics.. Vladimir A. Saleev (Dr. Sci. (Phys. & Math.)), Professor, Dept, of Mathematical Modeling in Mechanics..
