^ Гocnoдa физики, где барион SCC?

[Ost]

Ну не знаю. Значит я ошибся.
Вообще-то надо увеличивая параметр "o", уменьшать параметр "t", выдерживая слияние максимума и минимума.
Тогда эта точка слияния имеет минимальное значение по (No+N2)/2
 

При o = 0.6611;
       t = 0.774*1.0009; минимума тоже не наблюдается.
Точка слияния может наблюдаться при разных вариантах выбора o и t. Однозначности в решении нет.

[BJIaquMup]

Код: [Выделить]

to4 = 50;
s = 0.04168723670021172737215360293164504721700486232504;
ew = 30.1567307802595803785344923081668874541094417852214;
ev = 0.51099891000000;
No = 0.083;
N2 = 0.095;
r = 0;
o = 0.675;(* + *)
t = 0.633;(* - *)
q = (No + N2)/2;
Print["A = ", r + q];
Print["B = ", o + q];
Print["Y = ", t + q];
(* _______________________ *)
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
(* _______________________ *)
kv = 20;
kv1 = kv + 1;
x0 = 0; y0 = 0; ko := 0;
(* u = c, d = s *)
di := x^(x^(x^(x^B)));
ai := x^(x^(x^(x^(A^x))));
ui := x^(x^(x^(x^(Y^x))));
(* Proton *)
p[1] := ai^(ui^(1/di));
p[2] := ai^(di^(1/ui));
p[3] := (ui^(1/di))^ai;
p[4] := (di^(1/ui))^ai;
p[5] := (ui^di)^(1/ai);
p[6] := (di^ui)^(1/ai);
p[7] := ai^(1/(ui^di));
p[8] := ai^(1/(di^ui));
p[9] := ai^(1/(ui^(1/di)));
p[10] := ai^(1/(di^(1/ui)));
p[11] := (ui^(1/di))^(1/ai);
p[12] := (di^(1/ui))^(1/ai);
(* End Proton *)
(* ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ *)
cnucok := Array[g, kv];
m1 = 0.001; m2 = 0.4;
dd = (N2 - No)/kv;
For[j = 1, j < kv1,
    ko = 0;
    No = No + dd;
    A = No + r;
    B = No + o;
    Y = No + t;
    For[i = 1, i < 13,
      res = FindMinimum[p[i], {x, m1, m2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ko = ko*ev/4;
    g[j] = ko;
    j++];
ListPlot[cnucok, PlotJoined -> True, PlotStyle -> Hue[.6]];
Это для
o = 0.675;(* + *)
t = 0.633;(* - *)

А вот это

Код: [Выделить]

to4 = 50;
s = 0.04168723670021172737215360293164504721700486232504;
ew = 30.1567307802595803785344923081668874541094417852214;
ev = 0.510998910000000000000000000000000000000000000000;
No = 0.08300000000000000000000000000000000000000000000;
N2 = 0.096000000000000000000000000000000000000000000000;
r = 0.000000000000000000000000000000000000000000000000;;
o = 0.6497000000000000000000000000000000000000000000;(* - *)
t = 0.9290000000000000000000000000000000000000000000;(* + *)
q = (No + N2)/2;
Print["A = ", r + q];
Print["B = ", o + q];
Print["Y = ", t + q];
(* _______________________ *)
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
(* _______________________ *)
kv = 15;
kv1 = kv + 1;
(* u = c, d = s *)
di := x^(x^(x^(x^B)));
ai := x^(x^(x^(x^(A^x))));
ui := x^(x^(x^(x^(Y^x))));
(* Proton *)
p[1] := ai^(ui^(1/di));
p[2] := ai^(di^(1/ui));
p[3] := (ui^(1/di))^ai;
p[4] := (di^(1/ui))^ai;
p[5] := (ui^di)^(1/ai);
p[6] := (di^ui)^(1/ai);
p[7] := ai^(1/(ui^di));
p[8] := ai^(1/(di^ui));
p[9] := ai^(1/(ui^(1/di)));
p[10] := ai^(1/(di^(1/ui)));
p[11] := (ui^(1/di))^(1/ai);
p[12] := (di^(1/ui))^(1/ai);
(* End Proton *)
(* ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ *)
cnucok := Array[g, kv];
m1 = 0.010000000000000000000000000000000000000000000;
m2 = 0.400000000000000000000000000000000000000000000;
dd = (N2 - No)/kv;
For[j = 1, j < kv1,
    ko = 0;
    No = No + dd;
    A = No + r;
    B = No + o;
    Y = No + t;
    For[i = 1, i < 13,
      res = FindMinimum[p[i], {x, m1, m2}, AccuracyGoal -> 28,
    PrecisionGoal -> 20, WorkingPrecision -> to4];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ko = ko*ev;
    g[j] = ko;
    j++];
ListPlot[cnucok, PlotJoined -> True, PlotStyle -> Hue[.6]];
для
o = 0.6497000000000000000000000000000000000000000000;(* - *)
t = 0.9290000000000000000000000000000000000000000000;(* + *)

Смотрите значения эти для No и N2
и значения для той программы, которую дал ранее.

[Ost]

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

[BJIaquMup]

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

Так я же и говорю, по параметрам No & N2.

[Ost]

Так я же и говорю, по параметрам No & N2.

Если судить по программе, то параметр No в сумме с  (o, t, r),  определяет начальное значение параметров A, B, Y.
A = No + r;
B = No + o;
Y = No + t;

Можно было просто записать
A += dd;
B += dd;
Y += dd;
Задав начальные A, B, Y и выбросив (No, o, t, r) как лишние параметры.

Фактически Вы не ответили на мой вопрос "И как применить такое разнообразие для вычисления массы?"
Так как выбор No имеет смысл только в сумме с параметрами (o, t, r), а  N2 просто участвует в формировании
диапазона вычисления.
A ... An
B ... Bn
Y ... Yn.
n=kv.

В контексте сказанного, можно вопрос перефразировать.
Какие начальные значения A, B, Y из разных вариантов слияния определяют массу?
 

[BJIaquMup]

Фактически Вы не ответили на мой вопрос "И как применить такое разнообразие для вычисления массы?"

Ну, если уж вы так заинтересовались, то могу присоветовать вам всё же начать с более простого по исполнению. А именно, с мезонов. И с мезонов наиболее простых - ud.
Без этого я не смогу ответить на ваш вопрос.

Здесь же вы можете играться, как хотите. Мне без разницы. Я ни на чём не настаиваю. Возможно, я ошибаюсь. Вот вы и поправите.
Давайте, всё же, посмотрим, как это делается для мезонов? Там гораздо проще. Но только в другой теме.

[BJIaquMup]

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

Код: [Выделить]

s = 0.0416872367;
ew = 30.1567498951; ev = 0.51099891000000;
(* _______________________ *)
x := N[((1 + x0) + s)^(1/(1 - y0)), 20];
(* _______________________ *)
A = 0.083;
B = A + 0.6611;
Y = A + 0.774;
ko := 0;
(* d = s *)
di := x^(x^(x^(x^B)));
ai := x^(x^(x^(x^(A^x))));
ci := x^(x^(x^(x^(Y^x))));
(* Proton *)
p[1] := ai^(ci^(1/di));
p[2] := ai^(di^(1/ci));
p[3] := (ci^(1/di))^ai;
p[4] := (di^(1/ci))^ai;
p[5] := (ci^di)^(1/ai);
p[6] := (di^ci)^(1/ai);
p[7] := ai^(1/(ci^di));
p[8] := ai^(1/(di^ci));
p[9] := ai^(1/(ci^(1/di)));
p[10] := ai^(1/(di^(1/ci)));
p[11] := (ci^(1/di))^(1/ai);
p[12] := (di^(1/ci))^(1/ai);
(* End Proton *)
(* ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ *)
pa3 = 0;
m1 = 0.05; m2 = 0.7;
For[i = 1, i < 13,
    pa3++;
    (*Plot[p[i], {x, m1, m2}];*)
    res = FindMinimum[p[i], {x, m1, m2}];
    x0 = x /. Last[res];
    y0 = First[res];
    Print[pa3, " =  ", x0, "   ", y0];
    k1 = p[i]/ew; ko = ko + k1;
    Print[k1*ev, " MeV"];
    i++];
o = ko*ev;
Plot[p[1], {x, m1, m2}];
Plot[p[2], {x, m1, m2}];
Print["========================"];
Print["s(cs)  = ", o/4, " MeV"];


[BJIaquMup]

Немножко конкретизировал массу \(\Omega^+_{cc} = SCC\) бариона.

Очень тяжело вычислять, трудоёмко, но попробовал 4 ГэВ и 15 ГэВ.
Минимум слияния экстремумов заметно меньше на 15 ГэВ, чем на предсказанных (CASTRO) 4 ГэВ.

Но я однозначно неправ со старыми вычислениями. Они наивны. Вот здесь, на этой странице http://privaloff.narod.ru/cern/barion.html практически, вся таблица есть наивняк. Соответствует только в очень грубом приближении.

Но на протоне получается очень интересно. Там идёт соединение двух экстремумов. Вот на этом соединении и определяется масса бариона.
В принципе, барион SCC является тем же протоном, только состоящим из кварков второго поколения. Так что для него легче всего определиться точно.

Кстати, о наивной модели...  Чего-то мне помнится, что сам CASTRO как-то говаривал что-то о наивной модели...  Причём, не на форуме, а в своём научно-популярном трэде. Это относилось к трёхкварковому представлению барионов.

[BJIaquMup]

Уточнил массу ещё неоткрытого \(\Omega^+_{cc}\) SCC бариона.
Она находится приблизительно в пределах от 19 до 26 ГэВ.

Ранее вычисленное значение до 15 ГэВ неверно, т.к. вычислено без уточняющих процедур поиска минимумов (по умолчанию, там в вольфрамовской математике до 10 знаков, а это мало).
Не знаю, может и ошибся, но старая схема (а она заключается в одинаковых параметрах для всех кварков) не верна однозначно.

В этом свете вариант бариона \(\Omega^{++}_{ccc}\) мне не вычислить. А он наиболее интересен.
Здесь, с протоноподобными, просто повезло, что выявилось схождение экстремумов.

[BJIaquMup]

p.s.
Важнейшая поправка по протону и бариону SCC

Признаю себя ослом и посыпаю голову пеплом. В общем, я оказался неправ и барион SCC выколотить можно. Масса его немного не такая, какая предсказана Стандартной Моделью. И, тем не менее, крушения гипотезы нет. Скорее, наоборот.

Там очень хитровыделанная система выползла. Оказывается, в этой системе есть минимум на слиянии двух экстремумов (минимума и макимума). (А схема слияния минимума и максимума ровно такая же, как и на протоне). В принципе, барион SCC (имеется ввиду барион со спином 1/2, а не 3/2) является этаким "двойником" протона. Если составные кварки протона d- и u- из одного поколения кварков, то s- и c- кварки - из следующего.

Так получилось, что минимум на слиянии двух экстремумов (максимума и минимума) я заметил вначале на барионное SCC. Но такое явление вылезло и на протоне. (На протоне не сразу это заметишь). Теперь можно, ну, если не дезавуировать, то сильно-сильно поправить страницу "Вернуть протон" .

А вот теперь уже привязка к протону железная, стопроцентная. Здесь уже нельзя сослаться на какие-то «с потолка взятые» параметры. Здесь математика в чистом виде. Да-да, масса протона получается именно из чистейшей математики. Если принять в качестве гипотезы формулу для тройки лептонов: см. http://privaloff.narod.ru/ . Именно, что никаких параметров. Поскольку эти параметры, они сами выходят из минимума на слиянии двух экстремумов.
Только здесь не надо путать с минимумом на пи-мезоне, где минимум именно массы. Здесь же, на протоне, минимум на столкновении экстремумов. А значение массы вычисляется уже из этого минимума.

Ещё раз, значит. Надо учесть, что в случае барионов суммируются все возможные степенные связи. А их для протона всего 12 штук. (Для случая, например, бариона ddd, можно ограничиться 6-ю, а для барионов, например, udc, их 24). Поэтому, следует применять коэффициент при суммировании вариаций. Для протона он равен 1/2 (как и для бариона SCC). Коэффициент необходим, потому что в самом простом варианте, когда участвует только один аромат кварка, этот коэффициент равен 1.
Для протона, если коэффициент равен 1/2, то значение массы – 802.275 MeV
Для бариона \(\Omega^+_{cc} = SCC\) значение массы ожидается равным 53.63 GeV

[BJIaquMup]

Вот программа для вычисления значения числа слияния экстремумов для малой массы протона

Код: [Выделить]

to4 = 28;
s = 0.0416872367002117273721536029316450472170;
ew = 30.15673078025958037853449230816688745410;
w = 0.0260499593922058169319422986287300000000;
ev = 0.510998910000000000000000000000000000000;
No = 0.8;
N2 = 0.82;
Print["Ecmb CMEщение!"];
t = 0;
o = -0.04;
r = 0.0603;
Print[(No + N2)/2];
(* _______________________ *)
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
(* _______________________ *)
x0 = 0; y0 = 0;
ko := 0;
(* Quark *)
di := x^(x^B);
ai := x^(x^(A^x));
ui := x^(x^(Y^x));
(* Proton *)
p[1] := ai^(ui^(1/di));
p[2] := ai^(di^(1/ui));
p[3] := (ui^(1/di))^ai;
p[4] := (di^(1/ui))^ai;
p[5] := (ui^di)^(1/ai);
p[6] := (di^ui)^(1/ai);
p[7] := ai^(1/(ui^di));
p[8] := ai^(1/(di^ui));
p[9] := ai^(1/(ui^(1/di)));
p[10] := ai^(1/(di^(1/ui)));
p[11] := (ui^(1/di))^(1/ai);
p[12] := (di^(1/ui))^(1/ai);
(* End Proton *)
(* ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ *)
kvo = 50;
kv1 = kvo + 1;
cnucok := Array[g, kvo];
m1 = 0.100000000000000000000000000000000000;
m2 = 0.700000000000000000000000000000000000;
dd = (N2 - No)/kvo;
For[j = 1, j < kv1,
    ko = 0;
    No = No + dd;
    A = No + r;
    B = No + o;
    Y = No + t;
    For[i = 1, i < 13,
      (*res = FindMinimum[p[i], {x, m1, m2}, AccuracyGoal -> 28,
    PrecisionGoal -> 22, WorkingPrecision -> to4];*)
      res = FindMinimum[p[i], {x, m1, m2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ko = ko*ev;
    g[j] = ko/2;
    j++];
ListPlot[cnucok, PlotJoined -> True, PlotStyle -> Hue[.6]];
Print["Pro  = 938.272036"];
Это число здесь равно 0.81

А вот для большой массы

Код: [Выделить]

to4 = 28;
s = 0.0416872367002117273721536029316450472170;
ew = 30.15673078025958037853449230816688745410;
w = 0.0260499593922058169319422986287300000000;
ev = 0.510998910000000000000000000000000000000;
No = 0.80;
N2 = 0.822;
Print["Ecmb CMEщение!"];
t = 0;
o = 0.118;
r = 0.136;
Print[(No + N2)/2];
(* _______________________ *)
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
(* _______________________ *)
x0 = 0; y0 = 0;
ko := 0;
(* Quark *)
di := x^(x^B);
ai := x^(x^(A^x));
ui := x^(x^(Y^x));
(* Proton *)
p[1] := ai^(ui^(1/di));
p[2] := ai^(di^(1/ui));
p[3] := (ui^(1/di))^ai;
p[4] := (di^(1/ui))^ai;
p[5] := (ui^di)^(1/ai);
p[6] := (di^ui)^(1/ai);
p[7] := ai^(1/(ui^di));
p[8] := ai^(1/(di^ui));
p[9] := ai^(1/(ui^(1/di)));
p[10] := ai^(1/(di^(1/ui)));
p[11] := (ui^(1/di))^(1/ai);
p[12] := (di^(1/ui))^(1/ai);
(* End Proton *)
(* ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ *)
kvo = 50;
kv1 = kvo + 1;
cnucok := Array[g, kvo];
m1 = 0.100000000000000000000000000000000000;
m2 = 0.700000000000000000000000000000000000;
dd = (N2 - No)/kvo;
For[j = 1, j < kv1,
    ko = 0;
    No = No + dd;
    A = No + r;
    B = No + o;
    Y = No + t;
    For[i = 1, i < 13,
      (*res = FindMinimum[p[i], {x, m1, m2}, AccuracyGoal -> 28,
    PrecisionGoal -> 22, WorkingPrecision -> to4];*)
      res = FindMinimum[p[i], {x, m1, m2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ko = ko*ev;
    g[j] = ko/2;
    j++];
ListPlot[cnucok, PlotJoined -> True, PlotStyle -> Hue[.6]];
Print["Pro  = 938.272036"];
Здесь значение 0.811

[BJIaquMup]

А вот программа для собственно минимума по протону

Код: [Выделить]

to4 = 28;
s = 0.0416872367002117273721536029316450472170;
ew = 30.15673078025958037853449230816688745410;
w = 0.0260499593922058169319422986287300000000;
ev = 0.510998910000000000000000000000000000000;
h = 0.8734230000000000000000000000000000000000;
No = 0.798850000000000000000000000000000000000;
N2 = 0.799100000000000000000000000000000000000;
Print["Ecmb CMEщение!"];
t = 0;
o = 0.0734230000000000000000000000000000000000;
r = 0.1015274300000000000000000000000000000000;
Print[(No + N2)/2];
(* _______________________ *)
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
(* _______________________ *)
x0 = 0; y0 = 0;
ko := 0;
(* Quark *)
di := x^(x^B);
ai := x^(x^(A^x));
ui := x^(x^(Y^x));
(* Proton *)
p[1] := ai^(ui^(1/di));
p[2] := ai^(di^(1/ui));
p[3] := (ui^(1/di))^ai;
p[4] := (di^(1/ui))^ai;
p[5] := (ui^di)^(1/ai);
p[6] := (di^ui)^(1/ai);
p[7] := ai^(1/(ui^di));
p[8] := ai^(1/(di^ui));
p[9] := ai^(1/(ui^(1/di)));
p[10] := ai^(1/(di^(1/ui)));
p[11] := (ui^(1/di))^(1/ai);
p[12] := (di^(1/ui))^(1/ai);
(* End Proton *)
(* ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ *)
kvo = 50;
kv1 = kvo + 1;
cnucok := Array[g, kvo];
m1 = 0.100000000000000000000000000000000000;
m2 = 0.700000000000000000000000000000000000;
dd = (N2 - No)/kvo;
For[j = 1, j < kv1,
    ko = 0;
    No = No + dd;
    A = No + r;
    B = No + o;
    Y = No + t;
    For[i = 1, i < 13,
      res = FindMinimum[p[i], {x, m1, m2}, AccuracyGoal -> 28, PrecisionGoal \
-> 22, WorkingPrecision -> to4];
      (*res = FindMinimum[p[i], {x, m1, m2}];*)
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ko = ko*ev;
    g[j] = ko/2;
    j++];
ListPlot[cnucok, PlotJoined -> True, PlotStyle -> Hue[.6]];
Print["Pro  = 938.272036"];
Этот минимум равен 0.7989750
и соответствует значению массы 802.275 MeV

Можно хихикать, что не совпадает. Посмотрим, как у вас совпадёт SCC с 4 GeV

[BJIaquMup]

Для вариаций кварков из разных поколений такая закономерность не прослеживается. Например, для бариона c(dc) узел экстремумов есть, но минимум не прослеживается. А для u(su) таковой узел вообще не виден.
Это я имею ввиду протоноподобные барионы. А есть куча дельта-подобных барионов, с тем же допустим кварковым составом. Только d(cc) и s(uu). (Это барионы со спином 3/2).
Вычислительной работы тут море. Особенно, если учесть смешанные барионы с тремя различными ароматами кварков.

По идее, можно вычислить и гипотетический барион t(bt). Но программа даёт переполнение: слишком мала мощность моего слабого компьютера, приобретённого ещё добелоленточного движения. 

[VPD]

Для вариаций кварков из разных поколений такая закономерность не прослеживается. Например, для бариона c(dc) узел экстремумов есть, но минимум не прослеживается. А для u(su) таковой узел вообще не виден.
Это я имею ввиду протоноподобные барионы. А есть куча дельта-подобных барионов, с тем же допустим кварковым составом. Только d(cc) и s(uu). (Это барионы со спином 3/2).
Вычислительной работы тут море. Особенно, если учесть смешанные барионы с тремя различными ароматами кварков.

По идее, можно вычислить и гипотетический барион t(bt). Но программа даёт переполнение: слишком мала мощность моего слабого компьютера, приобретённого ещё добелоленточного движения. 

Может быть это, что Вы ищите?
http://www.dailytechinfo.org/news/9364-uchenye-cern-obnaruzhili-novuyu-dvazhdy-ocharovannuyu-elementarnuyu-chasticu.html

[BJIaquMup]

Может быть это, что Вы ищите?
http://www.dailytechinfo.org/news/9364-uchenye-cern-obnaruzhili-novuyu-dvazhdy-ocharovannuyu-elementarnuyu-chasticu.html

He, это не то. Это говорится вот про это
http://pdg.lbl.gov/2013/listings/rpp2013-list-xicc-plus.pdf
Про \(\Omega^-_{cc}\) там ни слова.

[BJIaquMup]

По идее, можно вычислить и гипотетический барион t(bt). Но программа даёт переполнение: слишком мала мощность моего слабого компьютера, приобретённого ещё добелоленточного движения. 

Кстати, просто ради антиресу, хто-нить, из штюдентов, не попробовали бы на своих мощных, купленных на папины деньги, компьютерах проверить, протестировать вариацию t(bt)? По образу и подобию u(du) и c(sc). A? 

[BJIaquMup]

Код: [Выделить]

to4 = 28;
s = 0.0412564319700000000000000000000000000000;
ew = 30.36623013531664800000000000000000000000;
ev = 0.510998910000000000000000000000000000000;
(*  *)
No = 0.08820000000000000000000000000000000000000000;
N2 = 0.0890000000000000000000000000000000000000000;
(*  *)
ba = 0;
(*  *)
on = 0.6583;
tw = 0.80971;
Print[(No + N2)/2];
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
x0 = 0; y0 = 0; ko := 0;
(*Quark*)
base := x^(x^(x^(x^(Ba^x))));
one := x^(x^(x^(x^Q1)));
two := x^(x^(x^(x^(q2^x))));
(*Barion*)
p[1] := base^(two^(1/one));
p[2] := base^(one^(1/two));
p[3] := (two^(1/one))^base;
p[4] := (one^(1/two))^base;
p[5] := (two^one)^(1/base);
p[6] := (one^two)^(1/base);
p[7] := base^(1/(two^one));
p[8] := base^(1/(one^two));
p[9] := base^(1/(two^(1/one)));
p[10] := base^(1/(one^(1/two)));
p[11] := (two^(1/one))^(1/base);
p[12] := (one^(1/two))^(1/base);
(*End Barion*)
kvo = 50;
kv1 = kvo + 1;
cnucok := Array[g, kvo];
m1 = 0.05;
m2 = 0.2;
n1 = 0.01;
n2 = 0.7;
dd = (N2 - No)/kvo;
For[j = 1, j < kv1,
    ko = 0;
    No = No + dd;
    Ba = No + ba;
    Q1 = No + on;
    q2 = No + tw;
    For[i = 1, i < 3,
      res = FindMinimum[p[i], {x, m1, m2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    For[i = 3, i <
         13,(*res = FindMinimum[p[i], {x, n1, n2}, AccuracyGoal -> 28,
      PrecisionGoal -> 22, WorkingPrecision -> to4];*)
      res = FindMinimum[p[i], {x, n1, n2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ko = ko*ev;
    g[j] = ko/2;
    j++];
ListPlot[cnucok,
            PlotJoined -> True, PlotStyle -> Hue[.6], ImageSize -> {400, 400}]
              ;

Поправка невелика, но результат тем не менее = 40481.4 MeV

Но есть куча вопросов. (В теме "Барионы" http://bolshoyforum.com/forum/index.php?topic=559752.0 ).

[BJIaquMup]

По бариону SCC получены важные поправки

Здесь почему-то даже проще, чем в случае с протоном.
Фишка с минимумом по оси абцисс проходит по всем трём вариациям.

(Подробнее, см. Тему "Барионы").

Результат для первой вариации, приблизительно 21.5 GeV. То же и для второй.
Для третьей - что-то около 5 МэВ, поэтому его можно практически не учитывать.
Окончательно, масса бариона SCC = 43 GeV.

Вот программы.

Для первых двух вариаций.

Код: [Выделить]

to4 = 28;
s = 0.0412564319700000000000000000000000000000;
ew = 30.36623013531664800000000000000000000000;
ev = 0.510998910000000000000000000000000000000;
(*  *)
No = 0.08;
N2 = 0.095;
(*  *)
on = 0;
(*  *)
tw = 0.587;
tr = 0.6535;
Print[(No + N2)/2];
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
x0 = 0; y0 = 0; ko := 0;
(*Quark*)
one := x^(x^(x^(x^(q1^x))));
two := x^(x^(x^(x^(q2^x))));
three := x^(x^(x^(x^q3)));
(*Barion*)
p[1] := one^(two^(1/three));
p[2] := one^(three^(1/two));
p[3] := two^(one^(1/three));
p[4] := two^(three^(1/one));
(*End Barion*)
Print["  "];
kvo = 50;
kv1 = kvo + 1;
cnucok := Array[g, kvo];
m1 = 0.050000000000000000000000000000000000;
m2 = 0.400000000000000000000000000000000000;
dd = (N2 - No)/kvo;
For[j = 1, j < kv1,
    ko = 0;
    No = No + dd;
    q1 = No + on;
    q2 = No + tw;
    q3 = No + tr;
    For[i = 1, i <
 5,(*res =
     FindMinimum[p[i], {x, m1, m2}, AccuracyGoal -> 28, PrecisionGoal -> 22,
      WorkingPrecision -> to4];*)
      res = FindMinimum[p[i], {x, m1, m2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ko = ko*ev/4;
    g[j] = ko;
    j++];
ListPlot[cnucok, PlotJoined -> True, PlotStyle -> Hue[.6], ImageSize -> {400,
    400}];

Для третьей вариации

Код: [Выделить]

to4 = 28;
s = 0.0412564319700000000000000000000000000000;
ew = 30.36623013531664800000000000000000000000;
ev = 0.510998910000000000000000000000000000000;
(*  *)
No = 0.76800000000000000000000000000000000000;
N2 = 0.77700000000000000000000000000000000000;
(*  *)
tr = 0;
(*  *)
on = -0.60000000000000000000000000000000000000;
tw = -0.62712000000000000000000000000000000000;
Print[(No + N2)/2];
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
x0 = 0; y0 = 0; ko := 0;
(*Quark*)
one := x^(x^(x^(x^(q1^x))));
two := x^(x^(x^(x^(q2^x))));
three := x^(x^(x^(x^q3)));
(*Barion*)
p[1] := one^(two^(1/three));
p[2] := one^(three^(1/two));
p[3] := two^(one^(1/three));
p[4] := two^(three^(1/one));
(*End Barion*)
Print["  "];
kvo = 50;
kv1 = kvo + 1;
cnucok := Array[g, kvo];
m1 = 0.050000000000000000000000000000000000;
m2 = 0.700000000000000000000000000000000000;
dd = (N2 - No)/kvo;
For[j = 1, j < kv1,
    ko = 0;
    No = No + dd;
    q1 = No + on;
    q2 = No + tw;
    q3 = No + tr;
    For[i = 1, i <
       5, res = FindMinimum[p[
        i], {x, m1, m2}, AccuracyGoal -> 28,
          PrecisionGoal -> 22, WorkingPrecision -> to4];
      (*res = FindMinimum[p[i], {x, m1, m2}];*)
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ko = ko*ev/4;
    g[j] = ko;
    j++];
ListPlot[cnucok, PlotJoined -> True, PlotStyle -> Hue[.6], ImageSize -> {400,
      400}];
Print["Macca  = ", 805.0 + 133.7, " MeV"];


[BJIaquMup]

Гы-гы-гы...  А ведь я нашел эту самую \(\Omega^+_{cc}\). Эту самую SCC. И причём, с кастриным значением!
 &/   
Причём, решение наипростейшее. (Как проглядел?  )
Уже не помню, сколько говорил CASTRO... Кажись, не более 4 ГэВ

Так вот. Получается, что  масса этой самой омеги около 3.668 GeV (3668 MeV).
 
Могу посчитать точнее.

  Вообще-то, это явилось для меня неожиданным. И ничего хорошего в этом нету.
Забавно конечно получится, если это значение с большой точностью совпадёт с экспериментальным. Но если оно действительно совпадёт, то тогда это полностью перечеркнёт всё, касательно барионов по гипотезе.

[BJIaquMup]

Ну вот, маленько подобрался к значению массы для гипотетического T(BT)-бариона. 

По образу и подобию протона и \( \Omega^+_{cc} \)-бариона.

Та же самая тильда. Пока выяснил 3-ю диспозицию. Она равна \( 8.17... * 10^7 \) MeV

RAM в 8 Гиг не позволяет точно выяснить, сколько в двух первых диспозициях (равных друг другу). Предположительно, около 10²⁰⁵ МэВ

Код: [Выделить]

to4 = 28;
s = 0.0415473362267507`35;
ew = 30.40588968977591515`35;
ev = 0.510998910`35;
(**)
No = 0.865`35;
N2 = 0.89`35;
(**)
tr = 0;
(**)
on = -0.71695`35;
tw = on;
Print[(No + N2)/2];
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
x0 = 0; y0 = 0; ko := 0;
(*Quark*)
one := x^(x^(x^(x^(x^(x^(q1^x))))));
two := x^(x^(x^(x^(x^(x^(q2^x))))));
three := x^(x^(x^(x^(x^(x^q3)))));
(*Barion*)
p[1] := one^(two^(1/three));
p[2] := one^(three^(1/two));
p[3] := two^(one^(1/three));
p[4] := two^(three^(1/one));
(*End Barion*)
Print["  "];
kvo = 50;
kv1 = kvo + 1;
cnucok := Array[g, kvo];
m1 = 0.0005`35;
m2 = 0.23`35;
dd = (N2 - No)/kvo;
For[j = 1, j < kv1, ko = 0;
    No = No + dd;
    q1 = No + on;
    q2 = No + tw;
    q3 = No + tr;
    For[i = 1, i < 5, res = FindMinimum[p[i], {x, m1, m2},
    AccuracyGoal -> Automatic, PrecisionGoal ->
Automatic, WorkingPrecision -> to4];
      (*res = FindMinimum[p[i], {x, m1, m2}];*)
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ko = ko*ev/4;
    g[j] = ko;
    j++];
ListPlot[cnucok, PlotJoined -> True, PlotStyle -> Hue[.6], ImageSize -> {400,
 400}];
Print["Macca  = ", 805.0 + 133.7, " MeV это для протона"];
Print["Здесь вычислен 3-ий вариант для массы T(BT)."];
Print[" Первые 2 не известны"];

Ну, в общем, для экспериментаторов, которые всё же сдуру захотят поискать  \( \Omega^+_{cc} \)-барион, цифирь всё та же == 53 ГэВ. +- 10 ГэВ. "Плюс бездетность, минус подоходные" (c)