Мезоны

[BJIaquMup]

Hy a теперь вишенка на тортике.

Если по чесноку, то в доказательство моей системы масс частиц я очень много ожидал от мезонов. Но мало получил.
Здесь хоть вроде бы и всё просто: кварк возводится в степень кварка поочерёдно, вариаций и всего две (в отличие от барионов). Кроме того здесь явно прорисовывается минимум. Но говорю, что толком оправдать ни один мезон так и не смог. Кроме ud. И то кое-как с грехом пополам привязался к единице: V=1.
Да, получился минимум с массой где-то около 90% от ожидаемого результата. Строго говоря, это малоутешительно.

Но тут, в очередной раз перебирая вариации q^(1/q) подумалось: "а почему бы не поменять местами V и A?"
Странно, мне почему-то казалось, что один хрен, что в лоб, что по лбу. Поменял. Оказалось — нефига подобного.
И интересен результат.
Вот прога:

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

to4 = 50; pa3 = 200;
S = 0.04125643197`50;
ew = 30.3662`50;
ev = 0.51099891`50;
V = 0.04125643197`50;
(* === === === === === === === === === === *)
x := N[((1 + x0) + S)^(1/(1 - y0)), to4];
(* === === === === === === === === === === *)
x0 = 0; y0 = 0;
ko := 0;
(* Quark *)
u := x^(x^(V^x));
d := x^(x^Y);
(*Meson*)
(*p[2] := d^u; даёт микроскопический вклад в массу *)
p[1] := u^d;
(*Calculate*)
x1 = 0.001`50;
x2 = 0.12`50;
b1 = 0.1`50;
b7 = 0.15`50;
Do[
    pb = b7 - b1; g = pb/7;
    b2 = b1 + g; b6 = b7 - g;
    Y = b1;
    ko = 0;
    For[i = 1, i < 2,
      res = FindMinimum[p[i], {x, x1, x2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ax1 = ko;
    Y = b1 + g;
    ko = 0;
    For[i = 1, i < 2,
      res = FindMinimum[p[i], {x, x1, x2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ax2 = ko;
    Y = b1 + 2*g;
    ko = 0;
    For[i = 1, i < 2,
      res = FindMinimum[p[i], {x, x1, x2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ax3 = ko;
    Y = b7 - 2*g;
    ko = 0;
    For[i = 1, i < 2,
      res = FindMinimum[p[i], {x, x1, x2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ax5 = ko;
    Y = b7 - g;
    ko = 0;
    For[i = 1, i < 2,
      res = FindMinimum[p[i], {x, x1, x2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ax6 = ko;
    Y = b7;
    ko = 0;
    For[i = 1, i < 2,
      res = FindMinimum[p[i], {x, x1, x2}];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]/ew; ko = ko + k1;
      i++];
    ax7 = ko;
    If[And[ax1 < ax2, ax2 < ax3], {b1 = b2}];
    If[And[ax7 < ax6, ax6 < ax5], {b7 = b6}],
    {pa3}];
o = ax3*ev;
Print["===   Это наглость!   ==="]
Print["Pi +- = 139.57018 MeV experiment"];
Print["Macca = ", o, " MeV"];
Print["========================="];
Print["V = ", V];
Print["A = ", Y];

Всё очень просто. Берём V = S и вычисляем. Находим число A. Оно в данном случае 0.11019... но это неважно.
Фишка в том, что это не минимум уже, а МАКСИМУМ!
И он равен массе пиона.        
Результат       == 139.813 MeV
Эксперимент  == 139.57018 MeV

Напомню, что
ew — это масса электрона в абсолютных единицах,
ev — это собственно масса электрона в MeV'ах.

Вот так, имея на руках полученное значение S, как результат игры с тройкой электрически заряженных лептонов, мы очень просто получаем массу Пи-мезона. Вот просто. С лёту. Блин! Я oxyeваю.

[BJIaquMup]

Вот чего получается пока что.

Всё-таки эта хренотень с равенством числа S и числа V — это весьма презабавная штуковина. Она выдаёт массу пиона. Только и всего-то. Нету ни минимума никакого, нихрена. Просто совпадение, и всё.
Пошарился в других мезонах — пока нuxy9l интересного.
Но есть такая хреновина, как неточное совпадение с единицей, когда ожидается, что при V = 1, масса пиона должна где-то барахтаться в районе экспериментального значения массы пиона. Да она там и барахтается, в районе 90% от экспериментальной.

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

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

pa3 = 200; vw = 37;
s = 0.04124014`50;
mes = 139.57018`50;
mso = 0.00035`50;
x := N[((1 + x0) + s)^(1/(1 - y0)), 50];
V = s;
x0 = 0; y0 = 0; ko := 0;
(* Quark *)
u := x^(x^(V^x));
d := x^(x^Y);
(*Meson*)
(*q := d^u; даёт небольшой вклад в массу *)
p := u^d;
q := d^u;
x1 = 0.001`50;
x2 = 0.2`50;
b1 = 0.000001`50;
b7 = 0.15`50;
Do[
    pb = b7 - b1; g = pb/7; b2 = b1 + g; b6 = b7 - g;
    Y = b1;
    ko =
    0; res = FindMinimum[p, {x, x1, x2}, AccuracyGoal -> Automatic, \
PrecisionGoal -> Automatic, WorkingPrecision -> vw];
    x0 = x /. Last[res]; y0 = First[res]; k1 = p; ko = ko + k1;
    ax1 = ko;
    Y = b1 + g;
    ko = 0; res = FindMinimum[p, {x, x1, x2}, AccuracyGoal -> Automatic,
    PrecisionGoal -> Automatic, WorkingPrecision -> vw];
    x0 = x /. Last[res]; y0 = First[res]; k1 = p; ko = ko + k1;
    ax2 = ko;
    Y = b1 + 2*g;
    ko = 0; res = FindMinimum[p, {x, x1, x2}, AccuracyGoal -> Automatic,
    PrecisionGoal -> Automatic, WorkingPrecision -> vw];
    x0 = x /. Last[res]; y0 = First[res]; k1 = p; ko = ko + k1;
    ax3 = ko;
    Y = b7 - 2*g;
    ko = 0; res = FindMinimum[p, {x, x1, x2},
    AccuracyGoal -> Automatic, PrecisionGoal ->
    Automatic, WorkingPrecision -> vw];
    x0 = x /. Last[res]; y0 = First[res]; k1 = p; ko = ko + k1;
    ax5 = ko;
    Y = b7 - g;
    ko = 0; res = FindMinimum[p, {x, x1, x2}, AccuracyGoal ->
        Automatic, PrecisionGoal -> Automatic, WorkingPrecision -> vw];
    x0 = x /. Last[res]; y0 = First[res]; k1 = p; ko = ko + k1;
    ax6 = ko;
    Y = b7;
    ko = 0; res = FindMinimum[p, {x, x1, x2}, AccuracyGoal ->
    Automatic, PrecisionGoal -> Automatic, WorkingPrecision -> vw];
    x0 = x /. Last[res]; y0 = First[res]; k1 = p; ko = ko + k1;
    ax7 = ko;
    If[And[ax1 < ax2, ax2 < ax3], {b1 = b2}];
    If[And[ax7 < ax6, ax6 < ax5], {b7 = b6}],
    {pa3}];
o = ax3;
Print["A = ", Y];
W = o/mes;
Print["W = ", W];
x1 = 0.00001`50;
x2 = 0.0002`50;
(*Plot[q, {x, x1, x2}, PlotStyle -> Hue[.6], ImageSize -> {400, 400}];*)
res = FindMinimum[q, {x, x1,
      x2}, AccuracyGoal -> Automatic, PrecisionGoal -> Automatic,
        WorkingPrecision -> vw];
x0 = x /. Last[res];
y0 = First[res];
Print["x = ", x0];
Print["y = ", y0];
Print["Macca Abs = ", q];
Print["Macca MeV = ", q/W];
o = o + q/W;
Print["q = ", q];
Print["q = ", q/W];
Print["========================="];
Print["Pi +  = ", mes + mso, " MeV experiment"];
Print["Macca = ", o/W, " MeV"];
Print["Pi -  = ", mes - mso, " MeV experiment"];
Print["========================="];
Print["Macca = ", o, " Absolut"];
Print["========================="];


[BJIaquMup]

Мы можем перейти к старой программе поиска минимума массы пиона, при V приближающемся к 1, но его не достигающем.
Вот эту добавку к  единице я никак не мог вменяемо объяснить. Окуда она, cyka, появляется? А без неё нет точности. Есть только эти 90%.
А вот теперь мы и показываем, ОТКУДА эта "добавочка" к единице. Вот этой вот программой:

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

to4 = 50; pa3 = 100; vw = 37;
s = 0.04124014`50;
W = 59.4106547958649980539418004451885112`50;
mes = 139.57018`50;
V = 1.0016054972`50;
(*_______________________*)
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
(*_______________________*)
x0 = 0; y0 = 0;
ko := 0;
(*Quark*)
ui := x^(x^(Y^x));
di := x^(x^V);
(*Meson*)
p[1] := di^ui;
p[2] := ui^di;
(*Calculate*)
(* ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^*)
x1 = 0.001`50;
x2 = 0.12`50;
b1 = 0.000001`50;
b7 = 0.0001`50;
Do[pb = b7 - b1; g = pb/7;
    b2 = b1 + g; b6 = b7 - g;
    Y = b1;
    ko = 0;
    For[i = 1, i < 3, res = FindMinimum[p[i], {x,
    x1, x2}, AccuracyGoal -> Automatic, PrecisionGoal ->
     Automatic, WorkingPrecision -> vw];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]; ko = ko + k1;
      i++];
    ax1 = ko;
    Y = b1 + g;
    ko = 0;
    For[i = 1, i < 3, res = FindMinimum[p[i], {x, x1, x2}, AccuracyGoal ->
 Automatic, PrecisionGoal -> Automatic, WorkingPrecision -> vw];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]; ko = ko + k1;
      i++];
    ax2 = ko;
    Y = b1 + 2*g;
    ko = 0;
    For[i = 1, i < 3, res = FindMinimum[p[i], {x, x1, x2}, AccuracyGoal ->
    Automatic, PrecisionGoal -> Automatic, WorkingPrecision -> vw];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]; ko = ko + k1;
      i++];
    ax3 = ko;
    Y = b7 - 2*g;
    ko = 0;
    For[i = 1, i < 3, res = FindMinimum[p[i], {
          x, x1, x2}, AccuracyGoal -> Automatic, PrecisionGoal ->
      Automatic, WorkingPrecision -> vw];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]; ko = ko + k1;
      i++];
    ax5 = ko;
    Y = b7 - g;
    ko = 0;
    For[i = 1, i < 3, res = FindMinimum[p[i], {x,
    x1, x2}, AccuracyGoal -> Automatic, PrecisionGoal -> Automatic,
    WorkingPrecision -> vw];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]; ko = ko + k1;
      i++];
    ax6 = ko;
    Y = b7;
    ko = 0;
    For[i = 1, i < 3, res = FindMinimum[p[i], {x, x1,
          x2}, AccuracyGoal -> Automatic, PrecisionGoal -> Automatic, \
WorkingPrecision -> vw];
      x0 = x /. Last[res];
      y0 = First[res];
      k1 = p[i]; ko = ko + k1;
      i++];
    ax7 = ko;
    If[And[ax1 > ax2, ax2 > ax3], {b1 = b2}];
    If[And[ax7 > ax6, ax6 > ax5], {b7 = b6}], {pa3}];
o = ax3/W;
Print["======================="]
Print["Pi +- = 139.57018"];
Print["Macca = ", o, " MeV"];
Print["========================"];
Print["V = ", V];
Print["A = ", Y];
Print[ax3];

В результате, V = 1.0016054972... , т.е. "добавочка" равна 0.0016054972... .

Вообще, пока это как бы не значит ничего. С дальнейшими мезонами пока глухо, как в Т-55.
Но тут два ништяка:
1. Куда более простая привязка всего лишь к одному мезону.
2. Возможность надыбать какую-то чисто математическую подоплёку ко всему этому 6ляqcmвy.