CC BY
87
УДК 531.9+539.12.01
ПРЕОНОВЫЕ ОБОЛОЧКИ И СТРУКТУРА АТОМА
Трунев Александр Петрович к.ф.-м.н., Ph.D.
Директор, A&E Trounev IT Consulting, Торонто, Канада
В работе рассмотрена модель структуры электронов и кварков, в которой эти частицы представляются состоящими из элементарных частиц преонов. На основе этой модели развита теория электронных оболочек, как продолжение ядерных кварковых оболочек
Ключевые слова: АТОМ, КВАРКИ, НЕЙТРОН, МАГНИТНЫЙ МОМЕНТ, МЕТРИКА, ПРЕОН, ПРОТОН, ЭЛЕКТРОН, ЭЛЕКТРОННЫЕ ОБОЛОЧКИ, ЭНЕРГИЯ СВЯЗИ, ЯДРО
UDC 531.9+539.12.01
PREONS SHELLS AND ATOMIC STRUCTURE
Alexander Trunev Cand.Phys.-Math.Sci., Ph.D.
Director, A&E Trounev IT Consulting, Toronto, Canada
We consider the model of the structure of electrons and quarks, in which these particles are presented consisting of elementary particles preons. From this model, the theory of electron shells, as a continuation of the quark nuclear shells has been proposed
Keywords: ATOM, BINDING ENERGY, QUARKS, MAGNETIC MOMENT, METRIC, PREON, PROTON, NEUTRON, NUCLEI, ELECTRON, ELECTRON SHELL
Введение
Согласно теории атомных и ядерных оболочек [1-4], периодические закономерности в атомах и атомных ядрах объясняются, принципом Паули, который применяется отдельно в отношении электронов и нуклонов, заполняющих ядерные и электронные оболочки соответственно. С другой стороны, в моделях квантовой хромодинамики, которые широко используются для моделирования адронов и атомных ядер [5-8], нуклоны представляются как составные частицы, состоящие из кварков.
В работах [9-11] сформулирована модель метрики адронов, удовлетворяющая основным требованиям физики элементарных частиц и космологии, а также рассмотрена динамика кварков, взаимодействующих с полем Янга-Миллса. Получены результаты по магнитным моментам барионов, согласующиеся с экспериментом с высокой точностью. В работе
[12] рассмотрено применение модели динамики кварков [8-9] к моделированию энергии связи нуклонов в атомных ядрах. Выведено уравнение энергии связи нуклонов в зависимости от содержания кварков. Показано, что кварки в ядрах образуют оболочки, аналогичные
электронным оболочкам. В настоящей работе развита объединенная модель электронных и ядерных оболочек, в которой электроны и кварки представляются как составные частицы, состоящие из элементарных частиц преонов. Показано, что в такой модели наблюдается симметрия ядерных и электронных оболочек.
Основные уравнения модели метрики адронов Рассмотрим центрально-симметричную метрику вида [7-13]
¥ = 1^0 О = - dt2 + e 2vdr2 + dв2 + а2(в^ф2
d2а (1)
= -ка
dв¿
со1 = dt, со2 = evdr, о0 = dв, О = аdj
Здесь = 1 - метрический тензор пространства Минковского сигнатуры (- + + +), к = сот1; - гауссова кривизна квадратичной формы dв2 + а2(в)dj2, Функция V = У(г, ?) определяется путем решения уравнений Янга-Миллса [13]. Всюду, где не оговорено, используется система единиц, в которой константы скорости света и Планка равны единице.
Среди всех решений уравнений Янга-Миллса, в случае метрики (1), есть такое, которое выражается через эллиптическую функцию Вейерштрасса
[13]. В этом случае уравнения модели приводятся к виду [9-12]:
Л„= 2( А2 -к2), еп = Л„ г = t ± г + т„
А = ^/ВДг/^/Й; ^2,2 3), (2)
1 к 1 к
Ь11 =-Ь22 = - А --, Ь33 = Ь44 = - А --, Ь12 = Ь21 = 0.
3 6 6 3
Здесь обозначено: ё2, ё3- инварианты функции Вейерштрасса, причем 2 2 = к л/Т2 ; го - свободный параметр, связанный с выбором начал
координат; Ьу + Ь// - 2(1 Ь у )Лу = Ту - тензор энергии-импульса материи. Отметим, что в этих обозначениях уравнения Эйнштейна имеют вид
Ь/ + Ь/,+ Ь1у = К/ (3)
Ь = 1Ьу ; К// - тензор Риччи.
В метрике (2) можно определить дефект решетки типа пузыря. В области пузыря считаем, что А2 = к2, а во внешней области решение зададим в виде (2), имеем
А2 = к2, еп = 0, Г < г0
а = ъ4Ир(т/ъ4И, ё1, ё 2), еп = Ат, Т|>Г) ^
На границах пузыря непрерывна функция Ли ее первая производная,
к = *Л2р(тй/ЪЛ2, ё1, ё 2), Аг= 0, Г = Г0 (5)
В частном случае решетки с инвариантами заданными в виде ё2 = ^12, ёз = 1, находим первый ноль и соответствующее значение
параметра метрики Г0 = 3 0449983, к = 2Л°38°34 . Отметим, что метрика во внутренней области пузыря является трехмерной, поскольку не содержит радиальной координаты. Действительно, используя уравнения (1) и (4), находим
¥ = ^2 + dв2 + соб2 (л[кв + в0)dj2 (6)
Аналогично строится решение для других корней второго уравнения (5). Все эти решения отличаются только размером пузыря, тогда как значение параметра к не меняется.
Всякий пузырь можно вывернуть наизнанку, просто изменив на противоположные неравенства (4). В этом случае можно определить метрику во внешней области пузыря, используя решение первого уравнения (2), так, чтобы метрика внешнего пространства совпала с метрикой нашей Вселенной [9]. Наконец, третий тип частиц можно
составить как комбинацию двух первых, в результате возникает пузырь, ограниченный оболочкой конечной толщины.
Преобразуем метрику (6) к стандартному виду. Для этого умножим обе части выражения (6) на постоянное число — k и введем новые переменные, отличающиеся от старых переменных на постоянный
множитель Vk, в результате находим
Y ® Y1 = dt2 — dO2 — sin2 Odj (7)
Метрика (7) использовалась для моделирования структуры барионов, в том числе протона и нейтрона [9-10], а также атомных ядер [11-12].
Динамика кварков
Для описания динамики кварков во внутренней области пузыря с метрикой вида (7) рассмотрим систему уравнений Дирака во внешнем поле Янга-Миллса. Отметим, что согласно (2) в метрике (7) тензор энергии импульса является постоянным. Следовательно, будем предполагать, что поле Янга-Миллса во внутренней области пузыря сводится к некоторой совокупности констант. В настоящей модели использованы три константы, а само поле описывается скалярным и векторным потенциалом
К = (fb, Abm )
Кроме того, будем учитывать электромагнитное поле, которое генерируют кварки. Используя результаты работы [14], преобразуем уравнение Дирака к криволинейным координатам (7). Имеем систему уравнений
if‘(Ñm + ЩаьК )¥а = mabVa (8)
Здесь обозначено 7m,qab, A,ya, mab - матрицы Дирака, параметры взаимодействия, векторный потенциал, волновая функция и эффективная
масса поля кварка а входящего в состав частицы Ь соответственно. Матрицы Дирака в метрике (7) имеют вид
■
Г1 0 0 0 " г 0 0 0 -1 е-9"
0 1 0 0 , у9 = 0 0 /е19 0
0 0 -1 0 0 /е -19 0 0
,0 0 0 -1 ч -1 е19 0 0 0 J
■
0
0
б1п в
- е1<р соб в
0 - б1п в е 9 соб вЛ
0 е9 соб в б1п в
есоб в 0 0
- бш в 0 0 J
В этих обозначениях оператор Дирака в метрике (7) можно представить в форме
■Л,
■
■ д , + Увдв+^~д9 Б1П в
Поскольку кварки обладают электрическим зарядом, они генерируют электромагнитное поле, посредством которого взаимодействуют друг с другом. Для описания этого взаимодействия используем уравнения квантовой электродинамики в форме
777 —
ЩаьУаУ^Уа = (д, -^ )А
(9)
«
Здесь а = е /Не - постоянная тонкой структуры, Уа = Уа У ,Уа сопряженный (по Эрмиту) вектор. Таким образом, предполагаем, что токи и заряды кварков суммируются, создавая коллективное поле, с которым кварки взаимодействуют в соответствии с уравнениями (8).
Система уравнений (8)-(9) использовалась для моделирования динамики кварков в случае барионов [10]. В простейшем случае, в котором учитывается только одно электромагнитное поле, модель содержит 15 нелинейных уравнений в частных производных. Для понижения порядка системы представим решение уравнений (8)-(9) в форме
__ -iwt+iLj
r fm Л
f2(d)eij
f3(q)
Kif,(0Yj.
(10)
Здесь L,w - проекция углового момента на выделенную ось и энергия системы соответственно. Система уравнений Дирака для случая представления решения в форме (10), приводится к виду,
fí = (L + VabAb sin в)( f1 COt в + /2) + /2 +
(mab + w - 4ab F b )(f3 Sin в - f4 COS в)
Л' = (L + qab Ab Sin в)(f1 - f2 COt в) - f2 COt в -
(mab +W- 4ab F b )(f3 COS в + f4 Sin в)
= (mab - W + 4abF b )(f1 Sin в - f2 COS в) +
(L + VabAb Sin в)( f3 COt в + /4) + f4 f4 = -(mab -W+ 4abF b )(f1 COS в + f2 Sin в) +
(L + qab Ab Sin в)(f3 - f4 COt в) - f4 COt в
(11)
Здесь предполагается, что Ab
A = A + A™,, Ф, = Ф +Ф
YM
Отметим, что масса и заряд являются индивидуальными для каждого кварка, а момент и энергия всей системы выбираются из условия образования стоячих волн вдоль меридиональной координаты. Вычисляя ток в левой части уравнения (9) и оператор набла в правой части, находим
уравнения, описывающие электродинамическую часть потенциала
ЩьУаУ У а = Щ ab
Ґ 4 Л
(l fi2
V i=1
= -ф''-ф' cotq
(12)
ЩаьУаҐУа = 2°УаЬ (f1f4 - f2f3 )a = -A - АЄ COt в +
A.
sin q
q
Уаї У a
а
а
Здесь по индексу а осуществляется суммирование по всех кваркам, входящим в систему. Таким образом, в случае барионов задача сводится к решению системы из 14 обыкновенных дифференциальных уравнений.
Как известно, электромагнитные свойства элементарных частиц характеризуются электрическим зарядом и магнитным моментом. Поэтому параметры поля Янга-Миллса, фигурирующие в уравнениях (11), должны быть связаны с величиной заряда и магнитного момента системы кварков, которые определяются следующим образом
ж/2
Ґ 4 \
О = I Мчу У У а = 4ж | а^іп вдаЬ £ /2
а
0 V г=1
(13)
2 „
ж/2
Ж / ¿вът2 в£ С!аЬ (/ь/4 - /2/3 )а
0 а
В качестве единицы измерения массы возьмем 1 МэВ, тогда параметры поля Янга-Миллса, векторный потенциал и энергия системы будут выражаться в единицах МэВ. Единицей магнитного момента в этом
случае является mq = 1.0219978 ц,в, где тв - магнетон Бора.
Сомножителем здесь выступает удвоенная масса электрона, выраженная в принятых единицах массы.
а
Модель нуклонов
Влияние векторного потенциала на параметры барионов исследовалось в работе [10]. Было установлено, что масштаб изменения параметров векторного поля Янга-Миллса не превышает 1 МэВ. Следовательно, можно исключить это поле из рассмотрения, заменив его скалярным потенциалом, влияющим на эффективную массу кварков [1112, 15]. Решение системы уравнений (11)-(12) с нулевым векторным
потенциалом Янга-Миллса можно получить в виде ряда по степеням параметра a. Для системы кварков основное состояние с нулевым моментом представляется в стандартном виде (10) с постоянными
функциями ft:
L ° f = fab , f2 = 0 f3 = f4 = gab (14)
В случае (14) система уравнений (11) с нулевым векторным потенциалом приводится к виду:
2gab + (mab - Wab )fab = 0 Wab = ~mab (1 5)
Вычисляя компоненты 4-вектора тока, и используя первое условие нормировки (13), находим
J0 = fib + gi = (1 + mlb )fÍ ,
Jj = 2fabgab Sin 0 = 2m ab fab SÍn 0,
1 (16)
4pJ 0 = 1, fb =--------—
" Ub 4p(1 + mab)
Используем полученные результаты для вычисления магнитных моментов нейтрона и протона. Общие свойства исследуемых нуклонов и кварков представлены в таблицах 1-2. Если предположить, что в составе протона кварки типа u имеют противоположно направленные спины, а в составе нейтрона кварки d имеют противоположно направленные спины, тогда магнитный момент протона зависит от эффективной массы d кварка, а
магнитный момент нейтрона зависит от эффективной массы u кварка. В
этих предположениях находим
2m q
mb / =--------ab 2 , b = n, p; a = u, d. (17)
^b *q 3(1 + ma2b) ^ (17)
В случае протона имеем Mp I ßq = L5544916 x10 3, соответственно
уравнение (17) имеет два корня
mdp = 0.00699556MeV; 142.948MeV . (18)
Для нейтрона ßn / ßq = -1 06479466 х10 3, а эффективная масса u кварка имеет два значения:
mun = 0.0023958MeV; 417.397MeV . (19)
Следовательно, в каждом случае имеем два корня уравнения (17). Один из них соответствует очень малой энергии кварков порядка нескольких кэВ. Второй корень (18) близок к массе заряженного пи-мезона - 139.57018 МэВ, а второй корень (19) близок к величине массы трех пи-мезонов.
Таблица 1. Свойства барионов
Syfflböt Spin C1I5J« Ы3££ RaywiNnniber GFartöi HyperdLafge Isospin QusricCaitait
Р : 1 93B.272Ü3 1 Ei 5594713 1 1 [[DcwiiQuäTk. UpQuaric. UpQuafc]]
F 1 -1 93 ЕЛ72Ü3 -1 i_5 Bi SS4713 -1 г [[DownQuarkBar. UpQuarfcBar. UpQuarkBar}]
Л : 0 &3SLS6536 1 -3.S260Si4i 1 г [[DownQiiaifc. Dour.Qu3ric. UpQuartcj]
JT 1 0 939j5G53G -1 -3.B26QEi4i -1 г [ [Din%T.QuäikBii. Dorvr.QuäikBsi. UpOuarkBar}}
Таблица 2. Свойства кварков
Symbol Spin Charge Xlass BarvonNumber Bottomness Chaim Hyp erch arge Iso Spin Strangeness Topness
u 1 2 22 1 0 0 1 1 0 0
2 j j j ]
ü 1 i 22 1 0 0 1 1 0 0
2 3 3 3 I
d 1 1 5.0 1 0 0 1 1 0 0
2 J J J ]
d 1 1 5.0 1 0 0 1 1 0 0
2 3 3 3 I
Моделирование энергии связи нуклонов в атомных ядрах
Как известно, нуклоны объединяются в атомные ядра под влиянием ядерных сил. Однако сами ядерные силы долгое время оставались загадкой, не смотря на многочисленные феноменологические модели [3-4]. Заметный прогресс в моделировании ядерных сил связан с развитием квантовой хромодинамики и численных моделей нуклонов и легких ядер [5-8, 16].
В работах [11-12] развита модель, которая позволяет объяснить природу ядерных сил динамикой кварков в метрике (7). Предполагается, что ядра состоят из оболочек с метрикой типа (7) и глюонного конденсата, для описания которого используются модели [17-20]. В работе [12] получено уравнение энергии связи нуклонов в атомных ядрах путем разложения по степеням плотности кварков в виде
2 2 2 2
Еь /(N + г) = трпр + тппн - дипи - дапа + дипи гд + гд
2 3 6 (20) + Яиапипагд + Я„р (пипа ) + йпрпр (пипа )
Здесь пр = г /(г + N), nN = N /(г + N) - плотность протонов и
нейтронов; пи = (2г + N)/(г + N), па = (г + 2N)/(г + N) - плотность
двух видов кварков; Qi - параметры взаимодействия кварков, учитывающих вклад кластеров дейтрона и альфа-частиц в энергию связи;
трг, mnN - суммарная масса протонов и нейтронов, входящих в состав
ядра; , - химический потенциал и и ё кварков соответственно; Г -
средний радиус кварковой оболочки.
Для замыкания модели положим
V = 3p3 = Z + N, V =VmU +v0nu3,Vd = V md + dnd2/3 (21)
Первое уравнение (21) описывает зависимость объема и среднего радиуса кварковой (ядерной) оболочки от числа нуклонов, второе и третье уравнения (21) описывают зависимость химического потенциала ферми-газа релятивистских частиц от их плотности [21]. Фигурирующий в этих уравнениях масштабный множитель зададим как среднюю величину максимальной энергии из выражений (18) и (19), следовательно
V0 = (142.948 + 417.397)/2 = 280.173MeV (22)
Модель (20)-(22) получена в следующих предположениях:
1) каждый нуклон в ядре диссоциирует на отдельные кварки, которые образуют кварковые оболочки. Эти оболочки заполняются последовательно, подобно тому, как заполняются электронные оболочки. Поскольку же нуклоны состоят из двух типов кварков, то существуют оболочки двух типов, которые заполняются и и ё кварками соответственно. В этом случае энергия связи, приходящаяся на один нуклон, зависит от концентрации кварков в оболочках и от энергии взаимодействия кварков;
2) кварки каждого типа образуют ферми-газ, обладающий химическим потенциалом как у релятивистских частиц. Отметим, что обычно это предположение относится к самим нуклонам, поэтому стандартная формула энергии связи содержит слагаемое, описывающее кинетическую энергию нуклонов [4, 22-24];
3) при диссоциации масса нуклона расходуется на возбуждение кинетической энергии кварков и на создание связей между кварками;
4) во внутренней области пузыря кварки объединяются в кластеры протонов, нейтронов, ядер дейтрона и альфа-частиц.
Отметим, что выражение (23) позволяет моделировать с высокой точностью энергию связи нуклонов в ядрах изотопов различных элементов [12]. Кластеры ядер легких элементов, кроме водорода и гелия, не учитываются в модели (20), но, очевидно, что их учет сводится только к переопределению параметров модели.
Полагая в уравнении (23) би = = 0 приходим к модели с тремя
параметрами - рис. 1, в которой выделен вклад парного взаимодействия кластеров протонов и нейтронов. Среднее значение параметра
взаимодействия кластеров нуклонов составляет (бпр) = 2 4321^еК , что
близко к массе и кварка. Среднее значение параметра взаимодействия
кварков составляет при этом всего {Яи^ = 0 1847МеУ .
Рис. 1. Зависимость параметров взаимодействия кварков и кластеров нуклонов в ядрах изотопов различных химических элементов от числа протонов.
Структура кварковых оболочек
Правило заполнения кварковых оболочек с учетом образования кластеров может быть выведено из представленной выше динамической модели [12]. Учитывая, что энергия связи кварков относительно мала, а энергия связи кластеров нуклонов приблизительно равна массе кварка, можно предположить, что наиболее вероятной комбинацией является сочетание трех кварков. При этом образуются кластеры, похожие по своим свойствам на нуклоны. Согласно первому уравнению (16) плотность кварков во внутренней области пузыря является постоянной. Следовательно, состояние такой системы определяется только энергией.
Сформулируем правило заполнения оболочек: если две частицы обладают энергией Ei каждый, то вероятность того, что третья частица
обладающая энергией Ei+1 образует с ними кластер, пропорциональна
величине - Ei+1 Ei2 (знак минус обусловлен тем, что энергия связи является отрицательной, тогда как вероятность является положительной величиной). Поскольку статистика кварков определяется распределением Ферми, то в результате приходим к модели, исследованной в работах [2527]:
- E++1 Ef =--------—------------ (23)
i+1 i exp[( Et -gq)/ T) + 1 (23)
Здесь Vq,T,K - энергия, химический потенциал, температура системы и параметр модели соответственно. Все размерные величины в модели (24) имеют размерность МэВ.
На рис. 2 представлена бифуркационная диаграмма модели (23), по которой определяется правило заполнения кварковых оболочек. Мы предполагаем, что вся диаграмма в целом описывает ядерные и
электронные оболочки. Действительно, как следует из данных, приведенных на рис. 2, существует два типа оболочек, которые
соответствуют малой и большой величине параметра K , а также два типа оболочек с малой и большой величиной отношения энергии к температуре при заданной величине параметра K .
Далее заметим, что в случае адиабатического расширения релятивистского газа фермионов выполняется соотношение [21]:
VT 3 = const (24)
Отсюда находим, что T = To ro / rq, где параметры To, ro
характеризуют состояние ядра. Следовательно, уровни энергии
определяются в модели (23) в следующей универсальной форме
К
Е.
exp( - х1 + Ь) + 1 ’
Х,Т 0 Г0 Ь = _тч
Т
(25)
г
Рис. 2. Бифуркационная диаграмма модели (23) иллюстрирующая правило заполнения оболочек в атомах и в ядрах (оболочки выделены синей и красной рамкой соответственно).
Адиабатические процессы связывают электронные и ядерные оболочки путем обмена частицами в реакциях бета-распада. Можно предположить, что процесс обмена частицами в указанных реакциях осуществляется за счет изменения масштаба оболочек при выполнении условия (24). Тогда из (25) находим, что масштаб плотности энергии изменяется как отношении радиуса электронной орбиты к радиусу ядра, что составляет порядка 104. Параметр модели при этом изменяется в
пределах 4.5 < К < 1000 . Следовательно, в случае адиабатических процессов имеем оценку
К = К 1( V гч)к . (26)
Здесь параметры К1, к характеризуют состояние атома.
Таким образом, теория кварковых оболочек позволяет объединить электронные и ядерные оболочки в рамках одной модели.
Теория преонов и симметрия электронных и ядерных оболочек
Возникает вопрос, если кварковые оболочки охватывают и электронные оболочки, то, что тогда представляют собой электроны? Рассмотрим теорию [28-31], согласно которой электроны и кварки состоят из элементарных частиц преонов. Свойства преонов и составных частиц -кварков и электронов, приведены в таблице 3.
Таблица 3. Свойства преонов и составных частиц [28-29]
Частица Символ Спин Заряд Состав
Преон а 1/3
Преон / -2/3
Преон 8 1/3
Антидипреон 0 1/3 {}
Антидипреон 0 -2/3 {а8 }
Антидипреон 0 1/3 {а/}
Кварк и 2/3 {а/ 8 }
Кварк d -1/3 {8/8}
Электрон е ~ -1 {/8 }
Приведенная выше модель кварковых оболочек (24)-(27) и ее бифуркационная диаграмма на рис. 2 не меняется при учете структуры кварков и электронов. Но правило заполнения относится в этом случае к
преонам, а не к кваркам и электронам. В результате такой замены возникает симметрия электронных и кварковых оболочек.
Сформулируем модель энергии связи с учетом наличия структуры электронов и кварков:
1) каждый нуклон в ядре диссоциирует на отдельные кварки, которые распадаются на преоны;
2) преоны каждого типа образуют ферми-газ, обладающий химическим потенциалом как у релятивистских частиц;
3) при диссоциации масса нуклона расходуется на возбуждение кинетической энергии преонов и на создание связей между преонами;
4) во внутренней области пузыря преоны объединяются в кластеры кварков, электронов, протонов, нейтронов, ядер дейтрона, альфа-частиц и других ядер;
5) существует симметрия электронных и ядерных оболочек заключающаяся в последовательности заполнения электронных и протонных оболочек.
С учетом сделанных предположений выражение энергии связи (20) принимает вид:
5
Еъ /(N + I) = трпр + тппы + тепе - ^д1п1 +... (27)
1=1
Здесь пе - плотность электронов, $1, пг - химический потенциал и плотность пяти видов преонов, фигурирующих в составе кварков и электронов. Уравнение (27) можно рассматривать как общую модель атома, с учетом энергии связи ядра и электронных оболочек в совокупности с моделью (25), описывающей правило заполнения ядерных и электронных оболочек. Однако реально вклад преонов в энергию связи нуклонов не отличается от вклада, который дают кварки, согласно уравнению (20). Поэтому сохраняются все результаты, полученные на основе модели кварковых оболочек [12].
Наконец, заметим, что модель преонов является сектором в расширении Стандартной модели на основе теории узлов [30-32]. В этом смысле предложенная выше модель преоновых оболочек является расширением стандартной квантовой теории, объясняющей строение атомов и атомных ядер [1-4]. Как известно, подобие в строении ядерных и электронных оболочек было установлено методами квантовой механики на основе принципа Паули [3-4]. В настоящей работе показано, что существует симметрия в организации ядерных и электронных оболочек, обусловленная наличием единого механизмом заполнения оболочек гипотетическими частицами преонами.
References
1. Зоммерфельд А. Строение атома и спектры. Том 1, 2. М.: ГИТТЛ, 1956.
2. Гейзенберг В. Замечания к теории атомного ядра// УФН (1), 1936.
3. Maria Goeppert-Mayer. On Closed Shells in Nuclei// Phys. Rev. Vol. 74; 1948; Phys. Rev. Vol. 75; 1949.
4. Иваненко Д.Д., Периодическая система химических элементов и атомное ядро // Д.И.Менделеев. Жизнь и труды, АН СССР, М., 1957, с.66-100.
5. Stephen R. Cotanch. THEORETICAL STUDIES OF HADRONS AND NUCLEI/ Rep. 27695-8202, North Carolina State University, Raleigh, North Carolina, 2006.
6. S. Durr, Z. Fodor, J. Frison et all. Ab Initio Determination of Light Hadron Masses// Science, 21 November 2008: Vol. 322, no. 5905 pp. 1224-1227.
7. G. Colangelo et al., Review of lattice results concerning low energy particle physics, Eur. Phys. J. C 71 (2011) 1695, arXiv:1011.4408.
8. H. Leutwyler. Progress in understanding Quantum Chromodynamics/ Schladming Winter School 2012. http://www.leutwyler.itp.unibe.ch/data/2012Schladming.pdf
9. Трунев А.П. Моделирование метрики адронов на основе уравнений Янга-Миллса // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2012. - №10(84). С. 874 - 887. -Режим доступа: http://ej.kubagro.ru/2012/10/pdf/68.pdf, 0,875 у.п.л.
10. Трунев А. П. Динамика кварков в метрике адронов и структура барионов // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2013. - №01(85). С. 525 - 542. -Режим доступа: http://ej.kubagro.ru/2013/01/pdf/42.pdf
11. Трунев А. П. Динамика кварков в метрике барионов и структура ядра // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2013. - №01(85). С. 623 - 636. -Режим доступа: http://ej.kubagro.ru/2013/01/pdf/49.pdf
12. Trunev A.P. QUARK DYNAMICS IN ATOMIC NUCLEI AND QUARK SHELLS// Network electronic scientific journal of the Kuban State Agrarian University (The Journal KubGAU) [electronic resource]. - Krasnodar KubGAU, 2012.
13. Л.Н. Кривоносов, В. А. Лукьянов. Полное решение уравнений Янга-Миллса для центрально-симметричной метрики// Journal of Siberian Federal University, Mathematics & Physics 2011, 4(3), 350-362.
14. V. Dzhunushaliev. Canonical conjugated Dirac equation in a curved space// arXiv:1202.5100, Feb. 25, 2012.
15. J.J.J. Kokkedee. The Quark Model. - W.A. Benjamin Inc., NY-Amsterdam, 1969.
16. S. Quaglioni, P. Navratil, R. Roth, and W. Horiuchi. From nucleons to nuclei to fusion reactions / / arXiv: 1203.0268 [nucl-th]
17. V. Dzhunushaliev. Scalar model of the glueball// Hadronic J. Suppl. 19, 185 (2004); http://arxiv.org/pdf/hep-ph/0312289v4.pdf
18. V. Dzhunushaliev. SU(3) glueball gluon condensate//arXiv:1110.1427 [hep-ph].
19. Трунев А.П. Моделирование массы адронов и энергии возбужденных состояний атомных ядер в модели глюонного конденсата // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2012. - №07(81). С. 545 - 554. - Режим доступа: http://ej.kubagro.ru/2012/07/pdf/40.pdf
20. Alexander Trunev. Hadrons mass spectrum and the gluon thermodynamics//Chaos and Correlation, Nov. 25, 2012, http://chaosandcorrelation.org/Chaos/CR 2 11 2012.pdf
21. Л. Д. Ландау, Е.М. Лифшиц. Теоретическая физика. Т. 5. - Москва, «Наука»,
1976.
22. P. Leboeuf. Regularity and chaos in the nuclear masses/ Lect. Notes Phys. 652, Springer, Berlin Heidelberg 2005, p.245, J. M. Arias and M. Lozano (Eds.).
23. Трунев А.П. Ядерные оболочки и периодический закон Д.И. Менделеева // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2012. - №05(79). С. 414 - 439. -Режим доступа: http://ej.kubagro.ru/2012/05/pdf/29.pdf
24. Трунев А.П. Ядерные оболочки и периодический закон Д.И.Менделеева. Часть
2. // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2012. - №07(81). С. 491 - 514. -Режим доступа: http://ej.kubagro.ru/2012/07/pdf/37.pdf
25. Volov D. B. Specific behavior of one chaotic dynamics near the fine-structure constant//arXiv:1205.6091,05/2012, http://arxiv.org/pdf/1302.3163.pdf
26. Трунев А.П. Бифуркация энергии связи и хаос в атомных ядрах// Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2012. - №05(79). С. 403 - 413. -Режим доступа: http://ej.kubagro.ru/2012/05/pdf/28.pdf
27. Трунев А.П. Хаос и фазовые переходы в атомных ядрах // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2012. - №07(81). С. 515 - 523. - Режим доступа: http://ej.kubagro.ru/2012/07/pdf/38.pdf
28. Jean-Jacques Dugne, Sverker Fredriksson, Johan Hansson, Enrico Predazzi. Preon Trinity - a new model of leptons and quarks// arXiv:hep-ph/9909569v3
29. Sverker Fredriksson. Preon Prophecies by the Standard Model// arXiv:hep-ph/0309213v2
30. Sundance O. Bilson-Thompson, Fotini Markopoulou, Lee Smolin. Quantum gravity and the standard model// arXiv:hep-th/0603022v2
31. Robert J. Finkelstein. The Preon Sector of the SLq(2) (Knot) Model //arXiv:1301.6440v1 [hep-th] 28 Jan 2013
32. Finkelstein R.J. An SLq(2) Extension of the Standard Model// arXiv:1205.1026v3
