Homec4science
Diffusion COSOlver (master) Browse

cosolver/Mathematica_scripts/coulomb_speedfunctions.nb
05a553ae087cmaster

Last Change
Raw File

coulomb_speedfunctions.nb
View Options

(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 11.1' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 27717, 774]
NotebookOptionsPosition[ 25792, 720]
NotebookOutlinePosition[ 26130, 735]
CellTagsIndexPosition[ 26087, 732]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[TextData[{
"Test and Field speed functions \n",
StyleBox["We compute here the closed analytical expressions of the speed \
functions of test and field operators [see Ji & Held, 2006 for details]\n\n\
Verified against the speed functions calculated in COSOlver code !", "Text"]
}], "Title",
CellChangeTimes->{{3.7831627498134527`*^9, 3.78316282626906*^9}, {
3.7838348269063063`*^9, 3.78383483374611*^9}, {3.7850461299935713`*^9,
3.785046141649485*^9}, {3.845459469423867*^9,
3.8454595137435923`*^9}},ExpressionUUID->"32487dfa-0b15-4223-aa2c-\
568b14b50d70"],
Cell[BoxData[
RowBox[{"ClearAll", ";"}]], "Input",
CellChangeTimes->{{3.7851491590993757`*^9, 3.785149159955426*^9},
3.8060326391345882`*^9},ExpressionUUID->"c68f53d3-e92e-4a7b-a10f-\
bc58b6716a3b"],
Cell[CellGroupData[{
Cell["Physical parameters", "Chapter",
CellChangeTimes->{{3.783835321502646*^9, 3.783835331801469*^9},
3.78383577771546*^9},ExpressionUUID->"46e4fc55-f689-420e-8ebb-\
43c6317e4213"],
Cell[BoxData[{
RowBox[{
RowBox[{"\[Sigma]", " ", "=", " ", "1"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Tau]", " ", "=", "1"}], ";"}], "\n",
RowBox[{
RowBox[{"\[Chi]", " ", "=", " ",
RowBox[{"Sqrt", "[",
FractionBox["\[Tau]", "\[Sigma]"], "]"}]}], ";"}], "\n"}], "Input",
CellChangeTimes->{{3.7838386271913443`*^9, 3.783838681513032*^9}, {
3.783838737903454*^9, 3.783838777885253*^9}, {3.783840034319172*^9,
3.783840037217824*^9}, {3.7838404761166058`*^9, 3.783840479011627*^9}, {
3.783840586754912*^9, 3.783840587314732*^9}, {3.783849270958695*^9,
3.783849271035843*^9}, {3.783854630378265*^9, 3.783854641830275*^9}, {
3.7838546743140697`*^9, 3.783854684967579*^9}, {3.7838573022143927`*^9,
3.7838573048687553`*^9}, {3.784954855108571*^9, 3.784954871389327*^9}, {
3.784956171579228*^9, 3.784956191514943*^9}, {3.784956835303032*^9,
3.784956839350119*^9}, {3.784957113706499*^9, 3.7849571266526413`*^9}, {
3.784957198725687*^9, 3.78495721569265*^9}, {3.784960781877899*^9,
3.784960805313635*^9}, {3.78496391250171*^9, 3.784963925890826*^9}, {
3.7849678053626747`*^9, 3.7849678092315817`*^9}, {3.784967856682729*^9,
3.784967859222958*^9}, {3.784968519941201*^9, 3.784968526529997*^9}, {
3.784968765715962*^9, 3.784968769128914*^9}, {3.784968825275635*^9,
3.784968826216567*^9}, {3.784981633498375*^9, 3.784981636215417*^9}, {
3.785043309554091*^9, 3.785043311765821*^9}, {3.7850439560140944`*^9,
3.78504397928918*^9}, {3.785044362130726*^9, 3.785044369926791*^9}, {
3.785045147580247*^9, 3.785045150089452*^9}, {3.785045189358987*^9,
3.785045190496707*^9}, {3.785045745112938*^9, 3.785045748739983*^9}, {
3.785048458591638*^9, 3.785048460919135*^9}, {3.785049669837163*^9,
3.785049672296723*^9}, {3.7851604319389563`*^9, 3.785160433224063*^9}, {
3.785476000445002*^9, 3.785476034414681*^9}, 3.806029330879081*^9, {
3.845459294336763*^9,
3.845459294945724*^9}},ExpressionUUID->"64c6dc69-0b79-4b63-afd3-\
c3339329aa8d"],
Cell[BoxData[""], "Input",
CellChangeTimes->{{3.783835337046235*^9, 3.7838353402875023`*^9}, {
3.783835928226749*^9, 3.783835928298624*^9},
3.783838596328367*^9},ExpressionUUID->"4e16b959-78c1-4e1d-bdd8-\
4fabeb19a910"]
}, Open ]],
Cell[CellGroupData[{
Cell[TextData[{
"Closed Analytical formulations for speed velocity integrals of ",
Cell[BoxData[
FormBox[
RowBox[{
RowBox[{"Erf",
RowBox[{"(", "x", ")"}]}], ":",
SuperscriptBox[
SubscriptBox["e", "ab"], "k"]}], TraditionalForm]],
FormatType->"TraditionalForm",ExpressionUUID->
"9c6d33a3-5e72-433f-a47e-f8ab4ec7036f"],
" and ",
Cell[BoxData[
FormBox[
SubscriptBox[
SuperscriptBox["E", "k"], "ab"], TraditionalForm]],
FormatType->"TraditionalForm",ExpressionUUID->
"3afd57b0-a572-4a9c-baf0-56db420e4cbd"]
}], "Chapter",
CellChangeTimes->{{3.783835375918396*^9, 3.783835411918047*^9}, {
3.783835463669812*^9,
3.783835564933565*^9}},ExpressionUUID->"95727edf-4af9-46c6-848c-\
b2294b3e2ec2"],
Cell[BoxData[
RowBox[{"\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["e",
RowBox[{"ab", ",", "k_"}]], ":=", " ",
FractionBox[
RowBox[{"Factorial", "[",
RowBox[{"k", "-",
RowBox[{"1", "/", "2"}]}], "]"}],
RowBox[{
RowBox[{"Factorial", "[",
RowBox[{
RowBox[{"-", "1"}], "/", "2"}], "]"}], "*",
SuperscriptBox[
RowBox[{"(",
RowBox[{"1", " ", "+", " ",
SuperscriptBox["\[Chi]", "2"]}], ")"}],
RowBox[{"k", "+",
RowBox[{"1", "/", "2"}]}]]}]]}], "\[IndentingNewLine]",
RowBox[{
SubscriptBox["E",
RowBox[{"ab", ",", "k_"}]], ":=", " ",
RowBox[{
FractionBox["\[Chi]", "2"], "*",
RowBox[{
SubsuperscriptBox["\[Sum]",
RowBox[{"j", "=", "0"}], "k"],
RowBox[{
FractionBox[
RowBox[{"Factorial", "[", "k", "]"}],
RowBox[{"Factorial", "[", "j", "]"}]], "*",
SubscriptBox["e",
RowBox[{"ab", ",", "j"}]]}]}]}]}], "\n"}]}]], "Input",
CellChangeTimes->{{3.7838355755556726`*^9, 3.7838355761281843`*^9}, {
3.783838984469953*^9, 3.783838992701725*^9}, {3.784950250517385*^9,
3.784950253596424*^9}, {3.7849508052333307`*^9, 3.784950806632309*^9}, {
3.784957179781451*^9, 3.784957192420298*^9}, {3.784981639967022*^9,
3.784981641023785*^9}, {3.785043350639456*^9, 3.785043351645976*^9}, {
3.785043981970006*^9,
3.785043984130218*^9}},ExpressionUUID->"3c8121a4-1359-4d54-995a-\
f7e8538eca38"]
}, Open ]],
Cell[CellGroupData[{
Cell["Laguerre coefficient", "Chapter",
CellChangeTimes->{{3.783835636820557*^9,
3.7838356414523773`*^9}},ExpressionUUID->"44240ef5-784d-4db2-baa4-\
62f2de21386c"],
Cell[BoxData[
RowBox[{
SubscriptBox["L",
RowBox[{"p_", ",", "j_", ",", "l_", " "}]], ":=", " ",
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{"-", "1"}], ")"}], "l"],
RowBox[{
RowBox[{
RowBox[{
RowBox[{"Factorial", "[",
RowBox[{"p", " ", "+", " ", "j", " ", "+", " ",
RowBox[{"1", "/", "2"}]}], "]"}], "/",
RowBox[{"Factorial", "[",
RowBox[{"j", "-", "l"}], "]"}]}], "/",
RowBox[{"Factorial", "[",
RowBox[{"l", "+", "p", " ", "+",
RowBox[{"1", "/", "2"}]}], "]"}]}], "/",
RowBox[{"Factorial", "[", "l", "]"}]}]}]}]], "Input",
CellChangeTimes->{{3.7838356538583593`*^9,
3.783835720196089*^9}},ExpressionUUID->"3c38910b-7396-4a07-b78a-\
90d351ae55e5"]
}, Open ]],
Cell[CellGroupData[{
Cell[TextData[{
"\n",
StyleBox["Test Speed Function:", "Title"]
}], "Chapter",
CellChangeTimes->{{3.783835348870504*^9, 3.783835355214663*^9}, {
3.783835633228549*^9, 3.783835646220426*^9}, {3.783835729451866*^9,
3.78383576269156*^9},
3.84545932102451*^9},ExpressionUUID->"6ea61762-dbe0-4a42-87e6-\
d055568cb1f9"],
Cell[BoxData[{
RowBox[{
RowBox[{
SubscriptBox["a", "ab0"], "=", " ",
RowBox[{
RowBox[{"-", " ", "4"}], "\[Sigma]", " ",
FractionBox[
RowBox[{"(",
RowBox[{"\[Tau]", "-", "1"}], ")"}], "\[Tau]"]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["a",
RowBox[{"ab1", ",", "p_", ",", "l_"}]], ":=", " ",
RowBox[{
RowBox[{"4",
FractionBox[
RowBox[{"l", " ", "\[Sigma]"}], "\[Tau]"],
RowBox[{"(",
RowBox[{"\[Tau]", " ", "-", "2"}], ")"}]}], " ", "-", " ",
RowBox[{"p",
RowBox[{"(", " ",
RowBox[{"p", " ", "-", " ",
RowBox[{"2", "\[Sigma]"}], " ", "+", " ",
RowBox[{"4", " ",
FractionBox["\[Sigma]", "\[Tau]"]}], " ", "+", " ", "1"}],
")"}]}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["a",
RowBox[{"ab2", ",", "p_", ",", "l_"}]], ":=", " ",
RowBox[{
FractionBox["\[Sigma]", "\[Tau]"],
RowBox[{"(",
RowBox[{
RowBox[{"4",
SuperscriptBox["l", "2"]}], " ", "+", " ",
RowBox[{"4", "l",
RowBox[{"(",
RowBox[{"p", "-", "1"}], ")"}]}], " ", "+", " ",
RowBox[{
FractionBox["3", "2"], "p",
RowBox[{"(",
RowBox[{"p", "-", "1"}], ")"}]}]}], ")"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["a",
RowBox[{"pab0", " "}]], "=", " ",
RowBox[{"4",
RowBox[{"(",
RowBox[{"\[Tau]", "-", "1"}], ")"}],
FractionBox[
RowBox[{"(",
RowBox[{"\[Sigma]", " ", "+", " ", "\[Tau]"}], ")"}], "\[Tau]"]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["a",
RowBox[{"pab1", ",", "p_", ",", "l_"}]], ":=", " ",
RowBox[{
RowBox[{"-", "2"}],
RowBox[{"(",
RowBox[{
RowBox[{"2", "l"}], " ", "+", "p"}], ")"}],
RowBox[{"(",
RowBox[{
RowBox[{
FractionBox["\[Sigma]", "\[Tau]"],
RowBox[{"(",
RowBox[{"\[Tau]", " ", "-", "2"}], ")"}]}], " ", "-", "1"}],
")"}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["a",
RowBox[{"pab2", ",", "p_", ",", "l_"}]], ":=", " ",
RowBox[{
FractionBox[
RowBox[{"-", " ", "\[Sigma]"}], "\[Tau]"],
RowBox[{"(",
RowBox[{
RowBox[{"4",
SuperscriptBox["l", "2"]}], " ", "+", " ",
RowBox[{"4", "l",
RowBox[{"(",
RowBox[{"p", "-", "1"}], ")"}]}], " ", "+", " ",
RowBox[{
FractionBox["3", "2"], "p",
RowBox[{"(",
RowBox[{"p", "-", "1"}], ")"}]}]}], ")"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
SubscriptBox["\[Alpha]",
RowBox[{"Eab", ",", "d_", ",", "j_", ",", "n_", ",", "g_"}]], ":=", " ",
SubscriptBox["E",
RowBox[{"ab", ",",
RowBox[{"d", "+", "j", "-", "n", "+", "g"}]}]]}], "\[IndentingNewLine]",
RowBox[{
SubscriptBox["\[Alpha]",
RowBox[{"eab", ",", "d_", ",", "j_", ",", "n_", ",", "g_"}]], ":=",
SubscriptBox["e",
RowBox[{"ab", ",",
RowBox[{
"d", "+", "j", "-", "n", " ", "+", "g", " ", "+", "1"}]}]]}], "\n",
RowBox[{
SubscriptBox["\[Nu]",
RowBox[{"abT", ",", "p_", ",", "j_", ",", "g_", ",", "d_"}]], ":=",
RowBox[{"4", "*",
RowBox[{
SubsuperscriptBox["\[Sum]",
RowBox[{"i", "=", "0"}], "j"],
RowBox[{
SubscriptBox["L",
RowBox[{"p", ",", "j", ",", "i"}]],
RowBox[{
FractionBox["1",
RowBox[{"Sqrt", "[", "Pi", "]"}]],
RowBox[{"(",
RowBox[{
RowBox[{
SubscriptBox["a", "ab0"], "*",
SubscriptBox["E",
RowBox[{"ab", ",",
RowBox[{"d", "+", "i", "+", "g"}]}]]}], "+", " ",
RowBox[{"\[Chi]", "*",
SubscriptBox["a", "pab0"], "*",
SubscriptBox["e",
RowBox[{"ab", ",",
RowBox[{"d", "+", "i", "+", "g", "+", "1"}]}]]}], " ", "+", " ",
RowBox[{
SubscriptBox["a",
RowBox[{"ab1", ",", "p", ",", "i"}]], "*",
SubscriptBox["E",
RowBox[{"ab", ",",
RowBox[{"d", "+", "i", "+", "g", "-", "1", " "}]}]]}], "+", " ",
RowBox[{"\[Chi]", " ", "*",
SubscriptBox["a",
RowBox[{"pab1", ",", "p", ",", "i"}]], "*",
SubscriptBox["e",
RowBox[{"ab", ",",
RowBox[{"d", "+", "i", "+", "g"}]}]]}], " ", "+", " ",
RowBox[{
SubscriptBox["a",
RowBox[{"ab2", ",", "p", ",", "i"}]],
SubscriptBox["E",
RowBox[{"ab", ",",
RowBox[{"d", "+", "i", "-", "2", "+", "g"}]}]]}], "+", " ",
RowBox[{"\[Chi]", "*",
SubscriptBox["a",
RowBox[{"pab2", ",", "p", ",", "i"}]],
SubscriptBox["e",
RowBox[{"ab", ",",
RowBox[{"d", "+", "i", "+", "g", "-", "1"}]}]]}]}],
")"}]}]}]}]}]}], "\[IndentingNewLine]",
RowBox[{" "}]}], "Input",
CellChangeTimes->{{3.78383578824758*^9, 3.783835919090938*^9}, {
3.7838359701386213`*^9, 3.783836102801579*^9}, {3.783836182617032*^9,
3.783836225192679*^9}, {3.783836421231509*^9, 3.783836511757928*^9}, {
3.78383654688647*^9, 3.78383655784647*^9}, {3.783836593118099*^9,
3.783836806303507*^9}, {3.783836915510477*^9, 3.7838371296370792`*^9}, {
3.783837171692642*^9, 3.783837293657864*^9}, {3.783837614871469*^9,
3.783837662566646*^9}, {3.783837716063058*^9, 3.783837716902781*^9}, {
3.783839537370204*^9, 3.783839541466062*^9}, {3.7838404942600193`*^9,
3.783840504379965*^9}, 3.783851675703457*^9, {3.783851758124722*^9,
3.78385175843633*^9}, {3.783852007140048*^9, 3.78385202055511*^9}, {
3.783852111643689*^9, 3.783852136987213*^9}, {3.78385234441084*^9,
3.783852351873611*^9}, 3.78385242073118*^9, {3.783852698733073*^9,
3.7838527043083572`*^9}, {3.7838528126561823`*^9, 3.783852865318777*^9}, {
3.783852985247779*^9, 3.7838530410319*^9}, {3.7838531163847637`*^9,
3.7838531228839903`*^9}, {3.783853169948185*^9, 3.7838531936194677`*^9}, {
3.783853237572598*^9, 3.783853253379936*^9}, {3.783853302196375*^9,
3.783853402841034*^9}, {3.7838535306440067`*^9, 3.78385361846215*^9}, {
3.783853850193218*^9, 3.783853856574594*^9}, {3.783853905208426*^9,
3.783853922574686*^9}, {3.783853952784111*^9, 3.783853969835998*^9}, {
3.783854191535039*^9, 3.783854191901479*^9}, {3.783854277813509*^9,
3.783854279411846*^9}, {3.783854328788238*^9, 3.7838543577407637`*^9}, {
3.7838544706689873`*^9, 3.7838545148912973`*^9}, {3.78495752929228*^9,
3.78495753194669*^9}, {3.784957700723815*^9, 3.7849577018914433`*^9}, {
3.7849579218219147`*^9, 3.7849579264319277`*^9}, {3.784957979525153*^9,
3.784957987328021*^9}, {3.784958532476543*^9, 3.784958557214192*^9}, {
3.784958934401826*^9, 3.7849589412723017`*^9}, {3.7849591050495787`*^9,
3.784959113455572*^9}, {3.784959640205521*^9, 3.784959644196735*^9}, {
3.784959778182014*^9, 3.7849598427157803`*^9}, {3.785043067396207*^9,
3.785043090762169*^9}, {3.785043937378212*^9, 3.785043942300021*^9}, {
3.785044698832774*^9,
3.785044699550734*^9}},ExpressionUUID->"14b8cbf3-d1ce-4aec-8b6c-\
b2429862c65f"],
Cell["", "Section",
CellChangeTimes->{{3.783840507579067*^9, 3.783840537922998*^9},
3.783841059407319*^9, 3.783866927350333*^9, {3.8454593130805817`*^9,
3.845459314064392*^9}},ExpressionUUID->"164807f6-0ca8-4015-9f35-\
fcd37522d514"]
}, Open ]],
Cell[CellGroupData[{
Cell[TextData[StyleBox["Field Speed Function:", "Title"]], "Chapter",
CellChangeTimes->{{3.783165476266239*^9, 3.783165487706551*^9},
3.783835774387526*^9, {3.845459287256754*^9,
3.845459317672579*^9}},ExpressionUUID->"47ab2a6d-f761-41d5-a332-\
bcd47e5bd9fe"],
Cell[BoxData[{
RowBox[{
RowBox[{
SubscriptBox["b", "ab0"], "=", " ",
RowBox[{"4", "*",
FractionBox[
SuperscriptBox["\[Tau]", "2"], "\[Sigma]"]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
SubscriptBox["b",
RowBox[{"ab1", ",", "p_"}]], ":=", " ",
RowBox[{
RowBox[{"-", "8"}], "*",
FractionBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["p", "2"], " ", "+", " ", "p", " ", "-", "1"}], ")"}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], " ", "-", "1"}], ")"}],
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], " ", "+", "3"}],
")"}]}]]}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["b",
RowBox[{
RowBox[{"+", "ab0"}], ",", "p_", ",", "l_"}]], ":=", " ",
RowBox[{
RowBox[{
RowBox[{"-", "8"}], "*", "p", "*",
RowBox[{"(",
RowBox[{"p", "-", "1"}], ")"}], "*",
FractionBox[
RowBox[{"(",
RowBox[{"l", "+", "1"}], ")"}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], " ", "+", "1"}], ")"}], "*",
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], "-", "1"}], ")"}]}]]}], " ", "-", " ",
RowBox[{"\[Tau]", "*", "8", "*",
FractionBox[
RowBox[{"(",
RowBox[{"1", " ", "+", "p", " ", "-", " ",
RowBox[{"\[Sigma]", "*", "p"}]}], ")"}],
RowBox[{"\[Sigma]", "*",
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], " ", "+", "1"}], ")"}]}]]}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["b",
RowBox[{
RowBox[{"+", "ab1"}], ",", "p_"}]], ":=", " ",
RowBox[{"8", "*",
RowBox[{"(",
RowBox[{"p", "+", "1"}], ")"}], "*",
FractionBox[
RowBox[{"(",
RowBox[{"p", "+", "2"}], ")"}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], " ", "+", "1"}], ")"}], "*",
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], "+", "3"}], ")"}]}]]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["b",
RowBox[{
RowBox[{"-", "ab0"}], ",", "p_", ",", "l_"}]], ":=",
RowBox[{
RowBox[{"4", "*",
RowBox[{"(",
RowBox[{"p", "+", "1"}], ")"}], "*",
RowBox[{"(",
RowBox[{"p", "+", "2"}], ")"}], "*",
FractionBox[
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], "+",
RowBox[{"2", "*", "l"}], "+", "3"}], ")"}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], "+", "1"}], ")"}], "*",
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], " ", "+", "3"}], ")"}]}]]}], " ", "+", " ",
RowBox[{"8", "*",
SuperscriptBox["\[Chi]", "2"], "*",
FractionBox[
RowBox[{"(",
RowBox[{"p", "-", " ",
RowBox[{"\[Sigma]", "*", "p"}], "-", "\[Sigma]"}], ")"}],
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], " ", "+", "1"}], ")"}]]}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["b",
RowBox[{
RowBox[{"-", "ab1"}], ",", "p_"}]], ":=", " ",
RowBox[{
RowBox[{"-", " ", "8"}], "*", "p", "*",
FractionBox[
RowBox[{"(",
RowBox[{"p", "-", "1"}], ")"}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], " ", "+", "1"}], ")"}], "*",
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], " ", "-", "1"}], ")"}]}]]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
SubscriptBox["\[Beta]",
RowBox[{"eab", ",", "g_", ",", "p_", ",", "d_", ",", "n_", ",", "j_"}]], ":=",
RowBox[{
SuperscriptBox["\[Chi]",
RowBox[{"p", "+",
RowBox[{"2", "*", "j"}], "+",
RowBox[{"2", "*", "n"}], " ", "-", "1"}]], "*",
SubscriptBox["e",
RowBox[{"ab", ",",
RowBox[{
"g", " ", "+", " ", "1", " ", "+", " ", "d", "+", "j", "+",
"n"}]}]]}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["\[Beta]",
RowBox[{
RowBox[{"+", "ab"}], ",", "g_", ",", "p_", ",", "d_", ",", "n_", ",",
"j_"}]], ":=", " ",
RowBox[{
SuperscriptBox["\[Chi]",
RowBox[{"p", "+",
RowBox[{"2", "*", "n"}], " ", "-", "1"}]],
RowBox[{
SubsuperscriptBox["\[Sum]",
RowBox[{"i", "=", "0"}], "j"],
FractionBox[
RowBox[{
RowBox[{"Factorial", "[", "j", "]"}], "*",
SuperscriptBox["\[Chi]",
RowBox[{"2", "*", "i"}]], "*",
SubscriptBox["e",
RowBox[{"ab", ",",
RowBox[{"p", " ", "+", "n", "+", "1", "+", "d", "+", "i"}]}]]}],
RowBox[{"2", "*",
RowBox[{"Factorial", "[", "i", "]"}]}]]}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
SubscriptBox["\[Beta]",
RowBox[{
RowBox[{"-", "Eeab"}], ",", "p_", ",", "d_", ",", "n_", ",", "j_"}]], ":=",
" ",
RowBox[{
SuperscriptBox["\[Chi]",
RowBox[{
RowBox[{"2", "*", "n"}], "-", "p", "-", "2"}]], "*",
RowBox[{"(", " ",
RowBox[{
RowBox[{
FractionBox[
RowBox[{"Factorial", "[",
RowBox[{"p", "+", "j", " ", "+",
RowBox[{"1", "/", "2"}]}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{"1", "/", "2"}], "]"}]], "*",
FractionBox[
SubscriptBox["E",
RowBox[{"ab", ",",
RowBox[{"d", "+", "n"}]}]], "2"]}], " ", "-", " ",
RowBox[{
SubsuperscriptBox["\[Sum]",
RowBox[{"i", "=", "0"}],
RowBox[{"p", "+", "j"}]],
RowBox[{
FractionBox[
SuperscriptBox["\[Chi]",
RowBox[{
RowBox[{"2", "*", "i"}], "+", "1"}]],
RowBox[{"Factorial", "[",
RowBox[{"i", " ", "+",
RowBox[{"1", "/", "2"}]}], "]"}]], "*",
RowBox[{"Factorial", "[",
RowBox[{"p", "+", "j", "+",
RowBox[{"1", "/", "2"}]}], "]"}], "*",
FractionBox[
SubscriptBox["e",
RowBox[{"ab", ",",
RowBox[{"d", "+", "n", "+", "i", "+", "1"}]}]], "2"]}]}]}],
")"}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
SubscriptBox["\[Nu]",
RowBox[{"abF", ",", "p_", ",", "j_", ",", "g_", ",", "d_"}]], ":=",
RowBox[{"4",
RowBox[{
SubsuperscriptBox["\[Sum]",
RowBox[{"i", "=", "0"}], "j"],
RowBox[{
SubscriptBox["L",
RowBox[{"p", ",", "j", ",", "i"}]],
RowBox[{
FractionBox["1",
RowBox[{"Sqrt", "[", "Pi", "]"}]],
RowBox[{"(",
RowBox[{
RowBox[{
SubscriptBox["b", "ab0"],
SubscriptBox["\[Beta]",
RowBox[{
"eab", ",", "g", ",", "p", ",", "d", ",", "0", ",", "i"}]]}], " ",
"+", " ",
RowBox[{
SubscriptBox["b",
RowBox[{
RowBox[{"+", "ab0"}], ",", "p", ",", "i"}]],
SubscriptBox["\[Beta]",
RowBox[{
RowBox[{"+", "ab"}], ",", "g", ",", "p", ",", "d", ",", "0", ",",
"i"}]]}], "+",
RowBox[{
SubscriptBox["b",
RowBox[{
RowBox[{"-", "ab0"}], ",", "p", ",", "i"}]],
SubscriptBox["\[Beta]",
RowBox[{
RowBox[{"-", "Eeab"}], ",", "p", ",", "d", ",", "0", ",", "i"}]]}],
" ", "-",
RowBox[{"8", "*",
FractionBox[
RowBox[{
RowBox[{"(",
RowBox[{
SuperscriptBox["p", "2"], " ", "+", " ", "p", " ", "-", "1"}],
")"}],
SubscriptBox["\[Beta]",
RowBox[{
"eab", ",", "g", ",", "p", ",", "d", ",", "1", ",", "i", " "}]],
" "}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], " ", "-", "1"}], ")"}],
RowBox[{"(",
RowBox[{
RowBox[{"2", "*", "p"}], " ", "+", "3"}], ")"}]}]]}], " ", "+",
" ",
RowBox[{
SubscriptBox["b",
RowBox[{
RowBox[{"+", "ab1"}], ",", "p"}]],
SubscriptBox["\[Beta]",
RowBox[{
RowBox[{"+", "ab"}], ",", "g", ",", "p", ",", "d", ",", "1", ",",
"i"}]]}], " ", "+",
RowBox[{
SubscriptBox["b",
RowBox[{
RowBox[{"-", "ab1"}], ",", "p"}]],
SubscriptBox["\[Beta]",
RowBox[{
RowBox[{"-", "Eeab"}], ",", "p", ",", "d", ",", "1", ",",
"i"}]]}]}], ")"}]}]}]}]}]}]}], "Input",
CellChangeTimes->{{3.783162849130638*^9, 3.783162935724229*^9}, {
3.783163953293124*^9, 3.783164021229011*^9}, {3.783164243250895*^9,
3.783164523154057*^9}, {3.783164621129231*^9, 3.783164745159897*^9}, {
3.783164842656008*^9, 3.7831650111659117`*^9}, {3.783165070745102*^9,
3.783165438105013*^9}, {3.7831655978463507`*^9, 3.783165604137204*^9}, {
3.78316567857146*^9, 3.783165689830372*^9}, {3.783165774172901*^9,
3.783165774442171*^9}, 3.7831659730535383`*^9, {3.783166134390232*^9,
3.7831661362749357`*^9}, 3.783166270903467*^9, {3.783166794707615*^9,
3.7831667959715014`*^9}, {3.7831680376816807`*^9, 3.78316804484888*^9}, {
3.783835585486446*^9, 3.7838355916774883`*^9}, {3.783836500735998*^9,
3.7838365076230288`*^9}, {3.783837688207013*^9, 3.7838377677023077`*^9}, {
3.783837900741445*^9, 3.783837941405252*^9}, {3.783838030052445*^9,
3.783838110788252*^9}, {3.7838383142186747`*^9, 3.783838363450362*^9}, {
3.783838466914262*^9, 3.78383846808955*^9}, {3.783854859434825*^9,
3.783854864815896*^9}, {3.783855296363172*^9, 3.783855304573185*^9}, {
3.7838553680025377`*^9, 3.7838554223326693`*^9}, {3.7838557459941263`*^9,
3.7838558283056498`*^9}, {3.783855863667*^9, 3.783855895673208*^9}, {
3.783855960002249*^9, 3.7838559994744377`*^9}, 3.7838560692770243`*^9, {
3.783856158576486*^9, 3.783856163871543*^9}, 3.7838565400137997`*^9,
3.7838566960764112`*^9, {3.78385681285964*^9, 3.783856822259541*^9}, {
3.783856922210801*^9, 3.783856969514193*^9}, {3.783857005442675*^9,
3.7838570104191*^9}, {3.784963885374926*^9, 3.784963902077915*^9}, {
3.7850431065158978`*^9, 3.78504314079982*^9}, {3.78504564662465*^9,
3.785045662329577*^9}, {3.785151669543252*^9,
3.785151671281046*^9}},ExpressionUUID->"6992090f-f457-4931-bbf1-\
83d05cf41736"],
Cell[BoxData[""], "Input",
CellChangeTimes->{{3.783855987466969*^9, 3.783855993493289*^9},
3.78385606719067*^9, {3.78385611335681*^9, 3.783856117316914*^9},
3.783857277266204*^9,
3.845459263854042*^9},ExpressionUUID->"54906d8c-0547-4e30-9b82-\
76ecbcb6401e"]
}, Open ]]
}, Open ]]
},
WindowSize->{1661, 997},
WindowMargins->{{Automatic, -705}, {18, Automatic}},
FrontEndVersion->"11.1 for Linux x86 (64-bit) (April 18, 2017)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)
(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 573, 10, 181, "Title", "ExpressionUUID" -> \
"32487dfa-0b15-4223-aa2c-568b14b50d70"],
Cell[1156, 34, 204, 4, 35, "Input", "ExpressionUUID" -> \
"c68f53d3-e92e-4a7b-a10f-bc58b6716a3b"],
Cell[CellGroupData[{
Cell[1385, 42, 186, 3, 66, "Chapter", "ExpressionUUID" -> \
"46e4fc55-f689-420e-8ebb-43c6317e4213"],
Cell[1574, 47, 2029, 33, 118, "Input", "ExpressionUUID" -> \
"64c6dc69-0b79-4b63-afd3-c3339329aa8d"],
Cell[3606, 82, 227, 4, 33, "Input", "ExpressionUUID" -> \
"4e16b959-78c1-4e1d-bdd8-4fabeb19a910"]
}, Open ]],
Cell[CellGroupData[{
Cell[3870, 91, 738, 22, 113, "Chapter", "ExpressionUUID" -> \
"95727edf-4af9-46c6-848c-b2294b3e2ec2"],
Cell[4611, 115, 1503, 41, 155, "Input", "ExpressionUUID" -> \
"3c8121a4-1359-4d54-995a-f7e8538eca38"]
}, Open ]],
Cell[CellGroupData[{
Cell[6151, 161, 167, 3, 66, "Chapter", "ExpressionUUID" -> \
"44240ef5-784d-4db2-baa4-62f2de21386c"],
Cell[6321, 166, 753, 22, 37, "Input", "ExpressionUUID" -> \
"3c38910b-7396-4a07-b78a-90d351ae55e5"]
}, Open ]],
Cell[CellGroupData[{
Cell[7111, 193, 327, 8, 123, "Chapter", "ExpressionUUID" -> \
"6ea61762-dbe0-4a42-87e6-d055568cb1f9"],
Cell[7441, 203, 7157, 186, 475, "Input", "ExpressionUUID" -> \
"14b8cbf3-d1ce-4aec-8b6c-b2429862c65f"],
Cell[14601, 391, 242, 4, 65, "Section", "ExpressionUUID" -> \
"164807f6-0ca8-4015-9f35-fcd37522d514"]
}, Open ]],
Cell[CellGroupData[{
Cell[14880, 400, 268, 4, 78, "Chapter", "ExpressionUUID" -> \
"47ab2a6d-f761-41d5-a332-bcd47e5bd9fe"],
Cell[15151, 406, 10339, 303, 706, "Input", "ExpressionUUID" -> \
"6992090f-f457-4931-bbf1-83d05cf41736"],
Cell[25493, 711, 271, 5, 33, "Input", "ExpressionUUID" -> \
"54906d8c-0547-4e30-9b82-76ecbcb6401e"]
}, Open ]]
}, Open ]]
}
]
*)
c4science · Help