summaryrefslogtreecommitdiff
path: root/mathemathica_fwm/RbXMDSSetup.nb
diff options
context:
space:
mode:
Diffstat (limited to 'mathemathica_fwm/RbXMDSSetup.nb')
-rwxr-xr-xmathemathica_fwm/RbXMDSSetup.nb1542
1 files changed, 908 insertions, 634 deletions
diff --git a/mathemathica_fwm/RbXMDSSetup.nb b/mathemathica_fwm/RbXMDSSetup.nb
index 70f2769..a8b7155 100755
--- a/mathemathica_fwm/RbXMDSSetup.nb
+++ b/mathemathica_fwm/RbXMDSSetup.nb
@@ -10,10 +10,10 @@
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 145, 7]
-NotebookDataLength[ 606935, 17647]
-NotebookOptionsPosition[ 598405, 17371]
-NotebookOutlinePosition[ 598769, 17387]
-CellTagsIndexPosition[ 598726, 17384]
+NotebookDataLength[ 619915, 17921]
+NotebookOptionsPosition[ 611739, 17657]
+NotebookOutlinePosition[ 612103, 17673]
+CellTagsIndexPosition[ 612060, 17670]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
@@ -801,7 +801,12 @@ Cell[BoxData[
RowBox[{"Expand", "@",
RowBox[{"Hamiltonian", "[",
RowBox[{"system", ",",
- RowBox[{"ElectricField", "\[Rule]", "field"}]}], "]"}]}], "/.", " ",
+ RowBox[{"ElectricField", "\[Rule]", "field"}], ",",
+ RowBox[{"MagneticField", "\[Rule]",
+ RowBox[{"{",
+ RowBox[{"0", ",", "0", ",",
+ RowBox[{"\[CapitalOmega]L", "/", "BohrMagneton"}]}], "}"}]}]}],
+ "]"}]}], "/.", " ",
RowBox[{
RowBox[{
RowBox[{"(",
@@ -819,9 +824,10 @@ Cell[BoxData[
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
+ RowBox[{"\[CapitalOmega]L", "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0", "0",
"0", "0", "0",
RowBox[{
RowBox[{
@@ -915,10 +921,14 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]]}], "0", "0", "0"},
{"0",
+ RowBox[{
+ FractionBox["\[CapitalOmega]L", "2"], "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0",
FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
- "0", "0", "0",
+ RowBox[{
+ SqrtBox["3"], " ", "\[CapitalOmega]L"}], "2"], "0", "0", "0",
RowBox[{
FractionBox[
RowBox[{
@@ -1063,8 +1073,8 @@ Cell[BoxData[
{"0", "0",
FractionBox[
RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
- "0",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0",
+ "\[CapitalOmega]L", "0",
RowBox[{
RowBox[{"-",
FractionBox[
@@ -1249,10 +1259,15 @@ Cell[BoxData[
SubscriptBox["\[Omega]", "4"]}]}], "]"}], " ",
SubscriptBox["\[CapitalOmega]", "4"]}]}], "0"},
{"0", "0", "0",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "2"]}], "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0",
FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
- "0",
+ RowBox[{
+ SqrtBox["3"], " ", "\[CapitalOmega]L"}], "2"], "0",
RowBox[{
RowBox[{"-",
FractionBox[
@@ -1396,9 +1411,11 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]]}]},
{"0", "0", "0", "0",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
+ RowBox[{
+ RowBox[{"-", "\[CapitalOmega]L"}], "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0", "0",
"0",
RowBox[{
RowBox[{
@@ -1492,8 +1509,13 @@ Cell[BoxData[
SubscriptBox["\[CapitalOmega]", "4"]}],
RowBox[{"4", " ",
SqrtBox["6"]}]]}], "0"},
- {"0", "0", "0", "0", "0",
- RowBox[{"-",
+ {"0",
+ FractionBox[
+ RowBox[{
+ SqrtBox["3"], " ", "\[CapitalOmega]L"}], "2"], "0", "0", "0",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "2"]}], "-",
FractionBox[
RowBox[{"5", " ",
RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0",
@@ -1645,7 +1667,7 @@ Cell[BoxData[
SubscriptBox["\[CapitalOmega]", "4"]}],
RowBox[{"8", " ",
SqrtBox["3"]}]]}], "0", "0"},
- {"0", "0", "0", "0", "0", "0",
+ {"0", "0", "\[CapitalOmega]L", "0", "0", "0",
RowBox[{"-",
FractionBox[
RowBox[{"5", " ",
@@ -1831,8 +1853,12 @@ Cell[BoxData[
RowBox[{"t", " ",
SubscriptBox["\[Omega]", "4"]}]}], "]"}], " ",
SubscriptBox["\[CapitalOmega]", "4"]}]}], "0"},
- {"0", "0", "0", "0", "0", "0", "0",
- RowBox[{"-",
+ {"0", "0", "0",
+ FractionBox[
+ RowBox[{
+ SqrtBox["3"], " ", "\[CapitalOmega]L"}], "2"], "0", "0", "0",
+ RowBox[{
+ FractionBox["\[CapitalOmega]L", "2"], "-",
FractionBox[
RowBox[{"5", " ",
RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0",
@@ -2126,6 +2152,8 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]]}], "0",
RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "6"]}], "+",
RowBox[{"Energy", "[", "1", "]"}], "-",
FractionBox[
RowBox[{"5", " ",
@@ -2477,6 +2505,7 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]]}], "0", "0", "0",
RowBox[{
+ FractionBox["\[CapitalOmega]L", "6"], "+",
RowBox[{"Energy", "[", "1", "]"}], "-",
FractionBox[
RowBox[{"5", " ",
@@ -2584,6 +2613,8 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["2"]}]]}], "0", "0", "0", "0", "0",
RowBox[{
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -2727,6 +2758,8 @@ Cell[BoxData[
SubscriptBox["\[Omega]", "4"]}]}], "]"}], " ",
SubscriptBox["\[CapitalOmega]", "4"]}]}], "0", "0", "0", "0", "0",
RowBox[{
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -3062,6 +3095,9 @@ Cell[BoxData[
SubscriptBox["\[CapitalOmega]", "4"]}]}], "0", "0", "0", "0", "0",
"0", "0",
RowBox[{
+ RowBox[{"-",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"]}], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -3170,6 +3206,9 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["2"]}]]}], "0", "0", "0", "0", "0", "0", "0",
RowBox[{
+ RowBox[{"-",
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"]}], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -3216,59 +3255,64 @@ Cell[BoxData[
RowBox[{
RowBox[{"HyperfineA", "[", "_", "]"}], "\[Rule]", ".2"}], ",",
RowBox[{
- RowBox[{"HyperfineB", "[", "_", "]"}], "\[Rule]", ".1"}]}], "}"}]}],
- ",",
+ RowBox[{"HyperfineB", "[", "_", "]"}], "\[Rule]", ".1"}], ",",
+ RowBox[{"\[CapitalOmega]L", "\[Rule]", ".1"}]}], "}"}]}], ",",
RowBox[{"ParityOffset", "\[Rule]", "False"}]}], "]"}]], "Input",
CellID->167259034],
Cell[BoxData[
- GraphicsBox[{{{{},
- LineBox[{{1.6, 0.15000000000000002`}, {2.4,
- 0.15000000000000002`}}]}, {{},
- LineBox[{{0.6, 0.15000000000000002`}, {1.4,
- 0.15000000000000002`}}]}, {{},
+ GraphicsBox[{{{{}, LineBox[{{1.6, 0.25}, {2.4, 0.25}}]}, {{},
+ LineBox[{{0.6, 0.2}, {1.4, 0.2}}]}, {{},
LineBox[{{-0.4, 0.15000000000000002`}, {0.4,
0.15000000000000002`}}]}, {{},
- LineBox[{{-1.4, 0.15000000000000002`}, {-0.6,
- 0.15000000000000002`}}]}, {{},
- LineBox[{{-2.4, 0.15000000000000002`}, {-1.6,
- 0.15000000000000002`}}]}, {{},
- LineBox[{{0.6, -0.25}, {1.4, -0.25}}]}, {{},
- LineBox[{{-0.4, -0.25}, {0.4, -0.25}}]}, {{},
- LineBox[{{-1.4, -0.25}, {-0.6, -0.25}}]}, {{},
- LineBox[{{0.6, 1.75}, {1.4, 1.75}}]}, {{},
+ LineBox[{{-1.4, 0.10000000000000002`}, {-0.6,
+ 0.10000000000000002`}}]}, {{},
+ LineBox[{{-2.4, 0.05000000000000002}, {-1.6,
+ 0.05000000000000002}}]}, {{}, LineBox[{{0.6, -0.3}, {1.4, -0.3}}]}, {{},
+ LineBox[{{-0.4, -0.25}, {0.4, -0.25}}]}, {{},
+ LineBox[{{-1.4, -0.2}, {-0.6, -0.2}}]}, {{},
+ LineBox[{{0.6, 1.7333333333333334`}, {1.4, 1.7333333333333334`}}]}, {{},
LineBox[{{-0.4, 1.75}, {0.4, 1.75}}]}, {{},
- LineBox[{{-1.4, 1.75}, {-0.6, 1.75}}]}, {{},
- LineBox[{{1.6, 3.775}, {2.4, 3.775}}]}, {{},
- LineBox[{{0.6, 3.775}, {1.4, 3.775}}]}, {{},
+ LineBox[{{-1.4, 1.7666666666666666`}, {-0.6,
+ 1.7666666666666666`}}]}, {{},
+ LineBox[{{1.6, 3.9083333333333337`}, {2.4, 3.9083333333333337`}}]}, {{},
+ LineBox[{{0.6, 3.8416666666666663`}, {1.4, 3.8416666666666663`}}]}, {{},
LineBox[{{-0.4, 3.775}, {0.4, 3.775}}]}, {{},
- LineBox[{{-1.4, 3.775}, {-0.6, 3.775}}]}, {{},
- LineBox[{{-2.4, 3.775}, {-1.6, 3.775}}]}}, {{}, {}, {}},
+ LineBox[{{-1.4, 3.708333333333333}, {-0.6, 3.708333333333333}}]}, {{},
+ LineBox[{{-2.4, 3.6416666666666666`}, {-1.6,
+ 3.6416666666666666`}}]}}, {{}, {}, {}},
{Arrowheads[{-0.029508196721311476`, 0.029508196721311476`}],
- ArrowBox[{{1.04, 0.15000000000000002`}, {1.96, 3.775}}],
- ArrowBox[{{0.04000000000000001, 0.15000000000000002`}, {0.96, 1.75}}],
- ArrowBox[{{0.04000000000000001, 0.15000000000000002`}, {0.96, 3.775}}],
- ArrowBox[{{-0.96, 0.15000000000000002`}, {-0.04000000000000001, 1.75}}],
- ArrowBox[{{-0.96, 0.15000000000000002`}, {-0.04000000000000001, 3.775}}],
- ArrowBox[{{-1.96, 0.15000000000000002`}, {-1.04, 1.75}}],
- ArrowBox[{{-1.96, 0.15000000000000002`}, {-1.04, 3.775}}],
- ArrowBox[{{1.04, -0.25}, {1.96, 3.7750000000000004`}}],
- ArrowBox[{{0.04000000000000001, -0.25}, {0.96, 1.75}}],
- ArrowBox[{{0.04000000000000001, -0.25}, {0.96, 3.7750000000000004`}}],
- ArrowBox[{{-0.96, -0.25}, {-0.04000000000000001, 1.75}}],
- ArrowBox[{{-0.96, -0.25}, {-0.04000000000000001, 3.7750000000000004`}}],
- ArrowBox[{{1.96, 0.15000000000000002`}, {1.04, 1.75}}],
- ArrowBox[{{0.96, 0.15000000000000002`}, {0.04000000000000001, 1.75}}],
- ArrowBox[{{0.96, -0.25}, {0.04000000000000001, 1.75}}],
- ArrowBox[{{-0.04000000000000001, 0.15000000000000002`}, {-0.96, 1.75}}],
- ArrowBox[{{-0.04000000000000001, -0.25}, {-0.96, 1.75}}],
- ArrowBox[{{1.96, 0.15000000000000002`}, {1.04, 3.775}}],
- ArrowBox[{{0.96, 0.15000000000000002`}, {0.04000000000000001, 3.775}}],
- ArrowBox[{{0.96, -0.25}, {0.04000000000000001, 3.7750000000000004`}}],
- ArrowBox[{{-0.04000000000000001, 0.15000000000000002`}, {-0.96, 3.775}}],
- ArrowBox[{{-0.04000000000000001, -0.25}, {-0.96, 3.7750000000000004`}}],
- ArrowBox[{{-1.04, 0.15000000000000002`}, {-1.96, 3.775}}],
- ArrowBox[{{-1.04, -0.25}, {-1.96, 3.7750000000000004`}}]},
+ ArrowBox[{{1.04, 0.2}, {1.96, 3.9083333333333337`}}],
+ ArrowBox[{{0.04000000000000001, 0.15000000000000002`}, {0.96,
+ 1.7333333333333334`}}],
+ ArrowBox[{{0.04000000000000001, 0.15000000000000002`}, {0.96,
+ 3.8416666666666663`}}],
+ ArrowBox[{{-0.96, 0.10000000000000002`}, {-0.04000000000000001, 1.75}}],
+ ArrowBox[{{-0.96, 0.10000000000000002`}, {-0.04000000000000001, 3.775}}],
+ ArrowBox[{{-1.96, 0.05000000000000002}, {-1.04, 1.7666666666666666`}}],
+ ArrowBox[{{-1.96, 0.05000000000000002}, {-1.04, 3.708333333333333}}],
+ ArrowBox[{{1., -0.3}, {1., 0.2}}],
+ ArrowBox[{{1.04, -0.3}, {1.96, 3.908333333333334}}],
+ ArrowBox[{{0., -0.25}, {0., 0.15000000000000002`}}],
+ ArrowBox[{{0.04000000000000001, -0.25}, {0.96, 1.7333333333333334`}}],
+ ArrowBox[{{0.04000000000000001, -0.25}, {0.96, 3.841666666666667}}],
+ ArrowBox[{{-1., -0.2}, {-1., 0.10000000000000003`}}],
+ ArrowBox[{{-0.96, -0.2}, {-0.04000000000000001, 1.75}}],
+ ArrowBox[{{-0.96, -0.2}, {-0.04000000000000001, 3.775}}],
+ ArrowBox[{{1.96, 0.25}, {1.04, 1.7333333333333334`}}],
+ ArrowBox[{{0.96, 0.2}, {0.04000000000000001, 1.75}}],
+ ArrowBox[{{0.96, -0.3}, {0.04000000000000001, 1.7499999999999998`}}],
+ ArrowBox[{{-0.04000000000000001, 0.15000000000000002`}, {-0.96,
+ 1.7666666666666666`}}],
+ ArrowBox[{{-0.04000000000000001, -0.25}, {-0.96, 1.7666666666666666`}}],
+ ArrowBox[{{1.96, 0.25}, {1.04, 3.8416666666666663`}}],
+ ArrowBox[{{0.96, 0.2}, {0.04000000000000001, 3.775}}],
+ ArrowBox[{{0.96, -0.3}, {0.04000000000000001, 3.7750000000000004`}}],
+ ArrowBox[{{-0.04000000000000001, 0.15000000000000002`}, {-0.96,
+ 3.708333333333333}}],
+ ArrowBox[{{-0.04000000000000001, -0.25}, {-0.96, 3.708333333333333}}],
+ ArrowBox[{{-1.04, 0.10000000000000002`}, {-1.96, 3.6416666666666666`}}],
+ ArrowBox[{{-1.04, -0.2}, {-1.96, 3.6416666666666666`}}]},
{PointSize[0.0225]}},
ImagePadding->{{2., 2}, {2., 2.}},
ImageSize->244.]], "Output"]
@@ -3666,162 +3710,6 @@ Cell[BoxData[
MatrixForm[BoxForm`e$]]]], "Output"]
}, Open ]],
-Cell[CellGroupData[{
-
-Cell[BoxData[
- RowBox[{"Tr", "[",
- RowBox[{
- RowBox[{"Hamiltonian", "[", "system", "]"}], ",", "List"}], "]"}]], "Input"],
-
-Cell[BoxData[
- RowBox[{"{",
- RowBox[{
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], ",",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], ",",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], ",",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], ",",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], ",",
- RowBox[{"-",
- FractionBox[
- RowBox[{"5", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], ",",
- RowBox[{"-",
- FractionBox[
- RowBox[{"5", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], ",",
- RowBox[{"-",
- FractionBox[
- RowBox[{"5", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], ",",
- RowBox[{
- RowBox[{"Energy", "[", "1", "]"}], "-",
- FractionBox[
- RowBox[{"5", " ",
- RowBox[{"HyperfineA", "[", "1", "]"}]}], "4"]}], ",",
- RowBox[{
- RowBox[{"Energy", "[", "1", "]"}], "-",
- FractionBox[
- RowBox[{"5", " ",
- RowBox[{"HyperfineA", "[", "1", "]"}]}], "4"]}], ",",
- RowBox[{
- RowBox[{"Energy", "[", "1", "]"}], "-",
- FractionBox[
- RowBox[{"5", " ",
- RowBox[{"HyperfineA", "[", "1", "]"}]}], "4"]}], ",",
- RowBox[{
- RowBox[{"Energy", "[", "2", "]"}], "-",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "2", "]"}]}], "4"], "-",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineB", "[", "2", "]"}]}], "4"]}], ",",
- RowBox[{
- RowBox[{"Energy", "[", "2", "]"}], "-",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "2", "]"}]}], "4"], "-",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineB", "[", "2", "]"}]}], "4"]}], ",",
- RowBox[{
- RowBox[{"Energy", "[", "2", "]"}], "-",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "2", "]"}]}], "4"], "-",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineB", "[", "2", "]"}]}], "4"]}], ",",
- RowBox[{
- RowBox[{"Energy", "[", "2", "]"}], "-",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "2", "]"}]}], "4"], "-",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineB", "[", "2", "]"}]}], "4"]}], ",",
- RowBox[{
- RowBox[{"Energy", "[", "2", "]"}], "-",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "2", "]"}]}], "4"], "-",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineB", "[", "2", "]"}]}], "4"]}]}], "}"}]], "Output"]
-}, Open ]],
-
-Cell[CellGroupData[{
-
-Cell[BoxData[
- RowBox[{"SelectState", "[",
- RowBox[{"system", ",", " ",
- RowBox[{
- RowBox[{"Label", "\[Equal]", "1"}], "&&",
- RowBox[{"F", "\[Equal]", "1"}]}]}], "]"}]], "Input"],
-
-Cell[BoxData[
- RowBox[{"AtomicState", "[",
- RowBox[{"1", ",",
- RowBox[{"J", "\[Rule]",
- FractionBox["1", "2"]}], ",",
- RowBox[{"L", "\[Rule]", "1"}], ",",
- RowBox[{"S", "\[Rule]",
- FractionBox["1", "2"]}], ",",
- RowBox[{"NuclearSpin", "\[Rule]",
- FractionBox["3", "2"]}], ",",
- RowBox[{"Parity", "\[Rule]", "Odd"}], ",",
- RowBox[{
- RowBox[{"BranchingRatio", "[", "0", "]"}], "\[Rule]", "1"}], ",",
- RowBox[{"F", "\[Rule]", "1"}], ",",
- RowBox[{"M", "\[Rule]", "1"}]}], "]"}]], "Output"]
-}, Open ]],
-
-Cell[CellGroupData[{
-
-Cell[BoxData[
- RowBox[{
- RowBox[{
- RowBox[{"Hamiltonian", "[",
- RowBox[{"{",
- RowBox[{"SelectState", "[",
- RowBox[{"system", ",", " ",
- RowBox[{
- RowBox[{"Label", "\[Equal]", "1"}], "&&",
- RowBox[{"F", "\[Equal]", "1"}]}]}], "]"}], "}"}], "]"}],
- "\[LeftDoubleBracket]",
- RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "-",
- RowBox[{
- RowBox[{"Hamiltonian", "[",
- RowBox[{"{",
- RowBox[{"SelectState", "[",
- RowBox[{"system", ",", " ",
- RowBox[{
- RowBox[{"Label", "\[Equal]", "0"}], "&&",
- RowBox[{"F", "\[Equal]", "2"}]}]}], "]"}], "}"}], "]"}],
- "\[LeftDoubleBracket]",
- RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}]}]], "Input"],
-
-Cell[BoxData[
- RowBox[{
- RowBox[{"Energy", "[", "1", "]"}], "-",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "-",
- FractionBox[
- RowBox[{"5", " ",
- RowBox[{"HyperfineA", "[", "1", "]"}]}], "4"]}]], "Output"]
-}, Open ]],
-
Cell[TextData[{
"Write ",
Cell[BoxData[
@@ -3893,9 +3781,10 @@ Cell[BoxData[
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
+ RowBox[{"\[CapitalOmega]L", "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0", "0",
"0", "0", "0", "0", "0", "0", "0",
RowBox[{
RowBox[{"-",
@@ -3926,10 +3815,25 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]]}], "0", "0", "0"},
{"0",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
- "0", "0", "0", "0", "0", "0", "0",
+ RowBox[{
+ FractionBox["\[CapitalOmega]L", "2"], "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0",
+ RowBox[{
+ FractionBox["1", "2"], " ",
+ SqrtBox["3"], " ",
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"\[ImaginaryI]", " ", "z", " ",
+ SubscriptBox["k", "1"]}], "-",
+ RowBox[{"\[ImaginaryI]", " ", "z", " ",
+ SubscriptBox["k", "2"]}], "-",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "1"]}], "+",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "2"]}]}]], " ", "\[CapitalOmega]L"}], "0",
+ "0", "0", "0", "0", "0", "0",
RowBox[{
RowBox[{
RowBox[{"-",
@@ -3957,8 +3861,18 @@ Cell[BoxData[
{"0", "0",
FractionBox[
RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
- "0",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"\[ImaginaryI]", " ", "z", " ",
+ SubscriptBox["k", "1"]}], "-",
+ RowBox[{"\[ImaginaryI]", " ", "z", " ",
+ SubscriptBox["k", "2"]}], "-",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "1"]}], "+",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "2"]}]}]], " ", "\[CapitalOmega]L"}], "0",
RowBox[{
RowBox[{"-",
FractionBox[
@@ -4012,10 +3926,25 @@ Cell[BoxData[
SubscriptBox["\[Phi]", "4"]}]], " ",
SubscriptBox["\[CapitalOmega]", "4"]}]}], "0"},
{"0", "0", "0",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
- "0",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "2"]}], "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0",
+ RowBox[{
+ FractionBox["1", "2"], " ",
+ SqrtBox["3"], " ",
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"\[ImaginaryI]", " ", "z", " ",
+ SubscriptBox["k", "1"]}], "-",
+ RowBox[{"\[ImaginaryI]", " ", "z", " ",
+ SubscriptBox["k", "2"]}], "-",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "1"]}], "+",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "2"]}]}]], " ", "\[CapitalOmega]L"}], "0",
RowBox[{
RowBox[{"-",
FractionBox[
@@ -4073,9 +4002,11 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]]}]},
{"0", "0", "0", "0",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
+ RowBox[{
+ RowBox[{"-", "\[CapitalOmega]L"}], "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0", "0",
"0",
RowBox[{
RowBox[{
@@ -4101,12 +4032,28 @@ Cell[BoxData[
RowBox[{"\[ImaginaryI]", " ", "t", " ",
SubscriptBox["\[Omega]", "2"]}]}]], " ",
SubscriptBox["\[CapitalOmega]", "2"]}]}], "0", "0", "0", "0", "0"},
- {"0", "0", "0", "0", "0",
+ {"0",
+ RowBox[{
+ FractionBox["1", "2"], " ",
+ SqrtBox["3"], " ",
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{
+ RowBox[{"-", "\[ImaginaryI]"}], " ", "z", " ",
+ SubscriptBox["k", "1"]}], "+",
+ RowBox[{"\[ImaginaryI]", " ", "z", " ",
+ SubscriptBox["k", "2"]}], "+",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "1"]}], "-",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "2"]}]}]], " ", "\[CapitalOmega]L"}], "0",
+ "0", "0",
RowBox[{
RowBox[{"-",
- FractionBox[
- RowBox[{"5", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "-",
+ FractionBox["\[CapitalOmega]L", "2"]}], "-",
+ FractionBox[
+ RowBox[{"5", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "-",
SubscriptBox["\[Omega]", "1"], "+",
SubscriptBox["\[Omega]", "2"]}], "0", "0", "0", "0", "0", "0", "0",
RowBox[{
@@ -4137,7 +4084,20 @@ Cell[BoxData[
SubscriptBox["\[CapitalOmega]", "4"]}],
RowBox[{"8", " ",
SqrtBox["3"]}]]}], "0", "0"},
- {"0", "0", "0", "0", "0", "0",
+ {"0", "0",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{
+ RowBox[{"-", "\[ImaginaryI]"}], " ", "z", " ",
+ SubscriptBox["k", "1"]}], "+",
+ RowBox[{"\[ImaginaryI]", " ", "z", " ",
+ SubscriptBox["k", "2"]}], "+",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "1"]}], "-",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "2"]}]}]], " ", "\[CapitalOmega]L"}], "0",
+ "0", "0",
RowBox[{
RowBox[{"-",
FractionBox[
@@ -4197,12 +4157,27 @@ Cell[BoxData[
RowBox[{"\[ImaginaryI]", " ", "t", " ",
SubscriptBox["\[Omega]", "2"]}]}]], " ",
SubscriptBox["\[CapitalOmega]", "4"]}]}], "0"},
- {"0", "0", "0", "0", "0", "0", "0",
+ {"0", "0", "0",
RowBox[{
- RowBox[{"-",
- FractionBox[
- RowBox[{"5", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "-",
+ FractionBox["1", "2"], " ",
+ SqrtBox["3"], " ",
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{
+ RowBox[{"-", "\[ImaginaryI]"}], " ", "z", " ",
+ SubscriptBox["k", "1"]}], "+",
+ RowBox[{"\[ImaginaryI]", " ", "z", " ",
+ SubscriptBox["k", "2"]}], "+",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "1"]}], "-",
+ RowBox[{"\[ImaginaryI]", " ", "t", " ",
+ SubscriptBox["\[Omega]", "2"]}]}]], " ", "\[CapitalOmega]L"}], "0",
+ "0", "0",
+ RowBox[{
+ FractionBox["\[CapitalOmega]L", "2"], "-",
+ FractionBox[
+ RowBox[{"5", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "-",
SubscriptBox["\[Omega]", "1"], "+",
SubscriptBox["\[Omega]", "2"]}], "0",
RowBox[{
@@ -4317,6 +4292,8 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]]}], "0",
RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "6"]}], "+",
RowBox[{"Energy", "[", "1", "]"}], "-",
FractionBox[
RowBox[{"5", " ",
@@ -4409,6 +4386,7 @@ Cell[BoxData[
SubscriptBox["\[Omega]", "2"]}]}]], " ",
SubscriptBox["\[CapitalOmega]", "2"]}]}], "0", "0", "0", "0", "0",
RowBox[{
+ FractionBox["\[CapitalOmega]L", "6"], "+",
RowBox[{"Energy", "[", "1", "]"}], "-",
FractionBox[
RowBox[{"5", " ",
@@ -4416,6 +4394,8 @@ Cell[BoxData[
SubscriptBox["\[Omega]", "1"]}], "0", "0", "0", "0", "0"},
{"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0",
RowBox[{
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -4457,6 +4437,8 @@ Cell[BoxData[
SqrtBox["6"]}]]}], "0", "0", "0", "0", "0", "0", "0", "0", "0", "0",
"0",
RowBox[{
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -4578,6 +4560,9 @@ Cell[BoxData[
SubscriptBox["\[CapitalOmega]", "4"]}]}], "0", "0", "0", "0", "0",
"0", "0",
RowBox[{
+ RowBox[{"-",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"]}], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -4644,6 +4629,9 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["2"]}]]}], "0", "0", "0", "0", "0", "0", "0",
RowBox[{
+ RowBox[{"-",
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"]}], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -4724,9 +4712,10 @@ Cell[BoxData[
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
+ RowBox[{"\[CapitalOmega]L", "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0", "0",
"0", "0", "0", "0", "0", "0", "0",
RowBox[{"-",
FractionBox[
@@ -4739,9 +4728,11 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]]}], "0", "0", "0"},
{"0",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
+ RowBox[{
+ FractionBox["\[CapitalOmega]L", "2"], "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0", "0",
"0", "0", "0", "0", "0", "0", "0",
RowBox[{
RowBox[{"-",
@@ -4775,9 +4766,12 @@ Cell[BoxData[
SubscriptBox["\[Phi]", "4"]}]], " ",
SubscriptBox["\[CapitalOmega]", "4"]}], "0"},
{"0", "0", "0",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "2"]}], "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0", "0",
"0",
RowBox[{"-",
FractionBox[
@@ -4800,9 +4794,11 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]]}]},
{"0", "0", "0", "0",
- FractionBox[
- RowBox[{"3", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "0", "0", "0", "0",
+ RowBox[{
+ RowBox[{"-", "\[CapitalOmega]L"}], "+",
+ FractionBox[
+ RowBox[{"3", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "0", "0", "0", "0",
"0",
RowBox[{
RowBox[{"-",
@@ -4815,9 +4811,10 @@ Cell[BoxData[
{"0", "0", "0", "0", "0",
RowBox[{
RowBox[{"-",
- FractionBox[
- RowBox[{"5", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "-",
+ FractionBox["\[CapitalOmega]L", "2"]}], "-",
+ FractionBox[
+ RowBox[{"5", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "-",
SubscriptBox["\[Omega]", "1"], "+",
SubscriptBox["\[Omega]", "2"]}], "0", "0", "0", "0", "0", "0", "0",
FractionBox[
@@ -4855,10 +4852,10 @@ Cell[BoxData[
SubscriptBox["\[CapitalOmega]", "3"]}], "0"},
{"0", "0", "0", "0", "0", "0", "0",
RowBox[{
- RowBox[{"-",
- FractionBox[
- RowBox[{"5", " ",
- RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"]}], "-",
+ FractionBox["\[CapitalOmega]L", "2"], "-",
+ FractionBox[
+ RowBox[{"5", " ",
+ RowBox[{"HyperfineA", "[", "0", "]"}]}], "4"], "-",
SubscriptBox["\[Omega]", "1"], "+",
SubscriptBox["\[Omega]", "2"]}], "0",
FractionBox[
@@ -4898,6 +4895,8 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]], "0",
RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "6"]}], "+",
RowBox[{"Energy", "[", "1", "]"}], "-",
FractionBox[
RowBox[{"5", " ",
@@ -4936,6 +4935,7 @@ Cell[BoxData[
SubscriptBox["\[Phi]", "1"]}]], " ",
SubscriptBox["\[CapitalOmega]", "1"]}], "0", "0", "0", "0", "0",
RowBox[{
+ FractionBox["\[CapitalOmega]L", "6"], "+",
RowBox[{"Energy", "[", "1", "]"}], "-",
FractionBox[
RowBox[{"5", " ",
@@ -4943,6 +4943,8 @@ Cell[BoxData[
SubscriptBox["\[Omega]", "1"]}], "0", "0", "0", "0", "0"},
{"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0",
RowBox[{
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -4965,6 +4967,8 @@ Cell[BoxData[
SqrtBox["6"]}]]}], "0", "0", "0", "0", "0", "0", "0", "0", "0", "0",
"0",
RowBox[{
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -5018,6 +5022,9 @@ Cell[BoxData[
SubscriptBox["\[CapitalOmega]", "3"]}], "0", "0", "0", "0", "0", "0",
"0",
RowBox[{
+ RowBox[{"-",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"]}], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -5047,6 +5054,9 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["2"]}]], "0", "0", "0", "0", "0", "0", "0",
RowBox[{
+ RowBox[{"-",
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"]}], "+",
RowBox[{"Energy", "[", "2", "]"}], "-",
FractionBox[
RowBox[{"3", " ",
@@ -5259,7 +5269,8 @@ Cell[BoxData[
Cell[BoxData[
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
- {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0",
+ {"\[CapitalOmega]L", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0",
+ "0",
RowBox[{"-",
FractionBox[
RowBox[{
@@ -5270,7 +5281,9 @@ Cell[BoxData[
SubscriptBox["\[CapitalOmega]", "4"]}],
RowBox[{"4", " ",
SqrtBox["6"]}]]}], "0", "0", "0"},
- {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0",
+ {"0",
+ FractionBox["\[CapitalOmega]L", "2"], "0", "0", "0", "0", "0", "0", "0",
+ "0", "0", "0", "0",
RowBox[{
RowBox[{"-",
FractionBox["1", "8"]}], " ",
@@ -5298,7 +5311,9 @@ Cell[BoxData[
RowBox[{"-", "\[ImaginaryI]"}], " ",
SubscriptBox["\[Phi]", "4"]}]], " ",
SubscriptBox["\[CapitalOmega]", "4"]}], "0"},
- {"0", "0", "0", "0", "0", "0", "0", "0", "0",
+ {"0", "0", "0",
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "2"]}], "0", "0", "0", "0", "0",
RowBox[{"-",
FractionBox[
RowBox[{
@@ -5319,7 +5334,8 @@ Cell[BoxData[
SubscriptBox["\[CapitalOmega]", "4"]}],
RowBox[{"4", " ",
SqrtBox["6"]}]]}]},
- {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0",
+ {"0", "0", "0", "0",
+ RowBox[{"-", "\[CapitalOmega]L"}], "0", "0", "0", "0", "0",
RowBox[{
RowBox[{"-",
FractionBox["1", "4"]}], " ",
@@ -5331,7 +5347,8 @@ Cell[BoxData[
{"0", "0", "0", "0", "0",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox["\[CapitalOmega]L", "2"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"]}], "0", "0", "0", "0", "0", "0", "0",
FractionBox[
RowBox[{
@@ -5365,8 +5382,8 @@ Cell[BoxData[
SubscriptBox["\[CapitalOmega]", "3"]}], "0"},
{"0", "0", "0", "0", "0", "0", "0",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox["\[CapitalOmega]L", "2"], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"]}], "0",
FractionBox[
RowBox[{
@@ -5404,7 +5421,9 @@ Cell[BoxData[
SubscriptBox["\[CapitalOmega]", "2"]}],
RowBox[{"4", " ",
SqrtBox["6"]}]], "0",
- RowBox[{"-",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "6"]}], "-",
SubscriptBox["\[Delta]", "1"]}], "0", "0", "0", "0", "0", "0", "0"},
{"0", "0", "0",
RowBox[{"-",
@@ -5434,12 +5453,14 @@ Cell[BoxData[
RowBox[{"\[ImaginaryI]", " ",
SubscriptBox["\[Phi]", "1"]}]], " ",
SubscriptBox["\[CapitalOmega]", "1"]}], "0", "0", "0", "0", "0",
- RowBox[{"-",
+ RowBox[{
+ FractionBox["\[CapitalOmega]L", "6"], "-",
SubscriptBox["\[Delta]", "1"]}], "0", "0", "0", "0", "0"},
{"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], "0", "0", "0", "0"},
{
@@ -5454,8 +5475,9 @@ Cell[BoxData[
SqrtBox["6"]}]]}], "0", "0", "0", "0", "0", "0", "0", "0", "0", "0",
"0",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], "0", "0", "0"},
{"0",
@@ -5496,7 +5518,9 @@ Cell[BoxData[
"0",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], "0"},
{"0", "0", "0",
@@ -5519,7 +5543,9 @@ Cell[BoxData[
SqrtBox["2"]}]], "0", "0", "0", "0", "0", "0", "0",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}]}
},
@@ -5570,41 +5596,44 @@ Cell[BoxData[
RowBox[{
RowBox[{"HyperfineA", "[", "_", "]"}], "\[Rule]", ".2"}], ",",
RowBox[{
- RowBox[{"HyperfineB", "[", "_", "]"}], "\[Rule]", ".1"}]}], "}"}]}],
- ",",
+ RowBox[{"HyperfineB", "[", "_", "]"}], "\[Rule]", ".1"}], ",",
+ RowBox[{"\[CapitalOmega]L", "\[Rule]", ".1"}]}], "}"}]}], ",",
RowBox[{"ParityOffset", "\[Rule]", "False"}]}], "]"}]], "Input",
CellID->102096636],
Cell[BoxData[
- GraphicsBox[{{{{}, LineBox[{{1.6, 0.}, {2.4, 0.}}]}, {{},
- LineBox[{{0.6, 0.}, {1.4, 0.}}]}, {{},
+ GraphicsBox[{{{{}, LineBox[{{1.6, 0.1}, {2.4, 0.1}}]}, {{},
+ LineBox[{{0.6, 0.05}, {1.4, 0.05}}]}, {{},
LineBox[{{-0.4, 0.}, {0.4, 0.}}]}, {{},
- LineBox[{{-1.4, 0.}, {-0.6, 0.}}]}, {{},
- LineBox[{{-2.4, 0.}, {-1.6, 0.}}]}, {{},
- LineBox[{{0.6, -0.4}, {1.4, -0.4}}]}, {{},
+ LineBox[{{-1.4, -0.05}, {-0.6, -0.05}}]}, {{},
+ LineBox[{{-2.4, -0.1}, {-1.6, -0.1}}]}, {{},
+ LineBox[{{0.6, -0.45}, {1.4, -0.45}}]}, {{},
LineBox[{{-0.4, -0.4}, {0.4, -0.4}}]}, {{},
- LineBox[{{-1.4, -0.4}, {-0.6, -0.4}}]}, {{},
- LineBox[{{0.6, 1.6}, {1.4, 1.6}}]}, {{},
+ LineBox[{{-1.4, -0.35000000000000003`}, {-0.6, \
+-0.35000000000000003`}}]}, {{},
+ LineBox[{{0.6, 1.5833333333333335`}, {1.4, 1.5833333333333335`}}]}, {{},
LineBox[{{-0.4, 1.6}, {0.4, 1.6}}]}, {{},
- LineBox[{{-1.4, 1.6}, {-0.6, 1.6}}]}, {{},
- LineBox[{{1.6, 3.625}, {2.4, 3.625}}]}, {{},
- LineBox[{{0.6, 3.625}, {1.4, 3.625}}]}, {{},
+ LineBox[{{-1.4, 1.6166666666666667`}, {-0.6,
+ 1.6166666666666667`}}]}, {{},
+ LineBox[{{1.6, 3.7583333333333337`}, {2.4, 3.7583333333333337`}}]}, {{},
+ LineBox[{{0.6, 3.6916666666666664`}, {1.4, 3.6916666666666664`}}]}, {{},
LineBox[{{-0.4, 3.625}, {0.4, 3.625}}]}, {{},
- LineBox[{{-1.4, 3.625}, {-0.6, 3.625}}]}, {{},
- LineBox[{{-2.4, 3.625}, {-1.6, 3.625}}]}}, {{}, {}, {}},
+ LineBox[{{-1.4, 3.558333333333333}, {-0.6, 3.558333333333333}}]}, {{},
+ LineBox[{{-2.4, 3.4916666666666667`}, {-1.6,
+ 3.4916666666666667`}}]}}, {{}, {}, {}},
{Arrowheads[{-0.029508196721311476`, 0.029508196721311476`}],
- ArrowBox[{{0.04000000000000001, 0.}, {0.96, 1.6}}],
- ArrowBox[{{-0.96, 0.}, {-0.04000000000000001, 1.6}}],
- ArrowBox[{{-1.96, 0.}, {-1.04, 1.6}}],
- ArrowBox[{{0.04000000000000001, -0.4}, {0.96, 1.6}}],
- ArrowBox[{{-0.96, -0.4}, {-0.04000000000000001, 1.6}}],
- ArrowBox[{{1.96, 0.}, {1.04, 3.625}}],
- ArrowBox[{{0.96, 0.}, {0.04000000000000001, 3.625}}],
- ArrowBox[{{0.96, -0.4}, {0.04000000000000001, 3.6250000000000004`}}],
- ArrowBox[{{-0.04000000000000001, 0.}, {-0.96, 3.625}}],
- ArrowBox[{{-0.04000000000000001, -0.4}, {-0.96, 3.6250000000000004`}}],
- ArrowBox[{{-1.04, 0.}, {-1.96, 3.625}}],
- ArrowBox[{{-1.04, -0.4}, {-1.96, 3.6250000000000004`}}]},
+ ArrowBox[{{0.04000000000000001, 0.}, {0.96, 1.5833333333333335`}}],
+ ArrowBox[{{-0.96, -0.05}, {-0.04000000000000001, 1.6}}],
+ ArrowBox[{{-1.96, -0.1}, {-1.04, 1.6166666666666667`}}],
+ ArrowBox[{{0.04000000000000001, -0.4}, {0.96, 1.5833333333333335`}}],
+ ArrowBox[{{-0.96, -0.35000000000000003`}, {-0.04000000000000001, 1.6}}],
+ ArrowBox[{{1.96, 0.1}, {1.04, 3.6916666666666664`}}],
+ ArrowBox[{{0.96, 0.05}, {0.04000000000000001, 3.625}}],
+ ArrowBox[{{0.96, -0.45}, {0.04000000000000001, 3.625}}],
+ ArrowBox[{{-0.04000000000000001, 0.}, {-0.96, 3.558333333333333}}],
+ ArrowBox[{{-0.04000000000000001, -0.4}, {-0.96, 3.558333333333333}}],
+ ArrowBox[{{-1.04, -0.05}, {-1.96, 3.4916666666666667`}}],
+ ArrowBox[{{-1.04, -0.35000000000000003`}, {-1.96, 3.4916666666666667`}}]},
{PointSize[0.0225]}},
ImagePadding->{{2., 2}, {2., 2.}},
ImageSize->244.]], "Output"]
@@ -7226,7 +7255,7 @@ Cell[BoxData[
Appearance->Automatic,
ButtonFunction:>OutputSizeLimit`ButtonFunction[
Function[{OutputSizeLimit`Dump`x$},
- TableForm[OutputSizeLimit`Dump`x$]], 27, 23290106009964751349, 5/
+ TableForm[OutputSizeLimit`Dump`x$]], 87, 23290106009964751349, 5/
2],
Enabled->True,
Evaluator->Automatic,
@@ -7245,7 +7274,7 @@ Cell[BoxData[
Appearance->Automatic,
ButtonFunction:>OutputSizeLimit`ButtonFunction[
Function[{OutputSizeLimit`Dump`x$},
- TableForm[OutputSizeLimit`Dump`x$]], 27, 23290106009964751349, 5
+ TableForm[OutputSizeLimit`Dump`x$]], 87, 23290106009964751349, 5
2],
Enabled->True,
Evaluator->Automatic,
@@ -7264,7 +7293,7 @@ Cell[BoxData[
Appearance->Automatic,
ButtonFunction:>OutputSizeLimit`ButtonFunction[
Function[{OutputSizeLimit`Dump`x$},
- TableForm[OutputSizeLimit`Dump`x$]], 27, 23290106009964751349,
+ TableForm[OutputSizeLimit`Dump`x$]], 87, 23290106009964751349,
Infinity],
Enabled->True,
Evaluator->Automatic,
@@ -7309,7 +7338,7 @@ Cell[BoxData[
FrameMargins->5],
Deploy,
DefaultBaseStyle->"Deploy"],
- Out[27]]], "Output"]
+ Out[87]]], "Output"]
}, Open ]],
Cell["\<\
@@ -7438,7 +7467,8 @@ Cell[BoxData[
SubscriptBox["\[Delta]", "_"], "\[Rule]",
RowBox[{"1.", " ",
SuperscriptBox["10", "7"]}]}], ",",
- RowBox[{"\[Gamma]t", "\[Rule]", "10."}]}], "}"}]}], "/.",
+ RowBox[{"\[Gamma]t", "\[Rule]", "10."}], ",",
+ RowBox[{"\[CapitalOmega]L", "\[Rule]", ".1"}]}], "}"}]}], "/.",
"atomicdata"}], "/.",
RowBox[{"Hertz", "\[Rule]", "1"}]}], "/.",
RowBox[{"Mega", "\[Rule]", "1"}]}], "/.",
@@ -8166,7 +8196,7 @@ Cell[BoxData[
Appearance->Automatic,
ButtonFunction:>OutputSizeLimit`ButtonFunction[
Function[{OutputSizeLimit`Dump`x$},
- TableForm[OutputSizeLimit`Dump`x$]], 34, 23290106009964751349, 5/
+ TableForm[OutputSizeLimit`Dump`x$]], 95, 23290106009964751349, 5/
2],
Enabled->True,
Evaluator->Automatic,
@@ -8185,7 +8215,7 @@ Cell[BoxData[
Appearance->Automatic,
ButtonFunction:>OutputSizeLimit`ButtonFunction[
Function[{OutputSizeLimit`Dump`x$},
- TableForm[OutputSizeLimit`Dump`x$]], 34, 23290106009964751349, 5
+ TableForm[OutputSizeLimit`Dump`x$]], 95, 23290106009964751349, 5
2],
Enabled->True,
Evaluator->Automatic,
@@ -8204,7 +8234,7 @@ Cell[BoxData[
Appearance->Automatic,
ButtonFunction:>OutputSizeLimit`ButtonFunction[
Function[{OutputSizeLimit`Dump`x$},
- TableForm[OutputSizeLimit`Dump`x$]], 34, 23290106009964751349,
+ TableForm[OutputSizeLimit`Dump`x$]], 95, 23290106009964751349,
Infinity],
Enabled->True,
Evaluator->Automatic,
@@ -8249,7 +8279,7 @@ Cell[BoxData[
FrameMargins->5],
Deploy,
DefaultBaseStyle->"Deploy"],
- Out[34]]], "Output"]
+ Out[95]]], "Output"]
}, Open ]],
Cell["\<\
@@ -8761,8 +8791,7 @@ Cell[BoxData[
Cell[TextData[{
"Find equations for the steady state by setting time derivatives to zero. We \
-also substitute in numerical values for hyperfine constants and natural \
-widths in units of ",
+also substitute in numerical values for natural widths in units of ",
Cell[BoxData[
FormBox[
SuperscriptBox["10", "6"], TraditionalForm]],
@@ -8893,12 +8922,22 @@ Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",", "2", ",", "2"}], "}"}]}]], "[", "t", "]"}]}],
RowBox[{"4", " ",
- SqrtBox["6"]}]], "-",
+ SqrtBox["6"]}]], "+",
+ RowBox[{"\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",", "2"}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"2", ",", "2", ",", "1"}], "}"}]}]], "[", "t", "]"}]}],
+ "-",
RowBox[{
RowBox[{"(",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -9005,12 +9044,22 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{
- RowBox[{"-",
- RowBox[{"(",
+ FractionBox["1", "2"], " ", "\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
- SubscriptBox["\[Delta]", "2"]}], ")"}]}], " ",
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",", "1"}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"0", ",", "1", ",", "1"}], "}"}]}]], "[", "t", "]"}]}],
+ "-",
+ RowBox[{
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "2"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
+ SubscriptBox["\[Delta]", "2"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
RowBox[{
@@ -9097,6 +9146,16 @@ Cell[BoxData[
RowBox[{"0", ",", "2", ",", "1"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "2", ",", "1"}], "}"}]}]], "[", "t", "]"}]}],
+ "+",
+ RowBox[{
+ FractionBox["1", "2"], " ", "\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",", "1"}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"2", ",", "2", ",", "0"}], "}"}]}]], "[", "t", "]"}]}],
"-",
RowBox[{
RowBox[{"(",
@@ -9362,9 +9421,13 @@ Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",", "2", ",", "0"}], "}"}]}]], "[", "t", "]"}]}],
RowBox[{"4", " ",
- SqrtBox["6"]}]], "+",
+ SqrtBox["6"]}]], "-",
RowBox[{
- SubscriptBox["\[Delta]", "1"], " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "6"]}], "-",
+ SubscriptBox["\[Delta]", "1"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
RowBox[{
@@ -9457,7 +9520,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -9629,11 +9694,22 @@ Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"-",
- RowBox[{"(",
+ FractionBox["1", "2"]}], " ", "\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
- SubscriptBox["\[Delta]", "2"]}], ")"}]}], " ",
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",",
+ RowBox[{"-", "1"}]}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"0", ",", "1", ",",
+ RowBox[{"-", "1"}]}], "}"}]}]], "[", "t", "]"}]}], "-",
+ RowBox[{
+ RowBox[{"(",
+ RowBox[{
+ FractionBox["\[CapitalOmega]L", "2"], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
+ SubscriptBox["\[Delta]", "2"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
RowBox[{
@@ -9759,7 +9835,18 @@ Cell[BoxData[
RowBox[{"0", ",", "2", ",",
RowBox[{"-", "1"}]}], "}"}]}]], "[", "t", "]"}]}],
RowBox[{"4", " ",
- SqrtBox["2"]}]], "+",
+ SqrtBox["2"]}]], "-",
+ RowBox[{
+ FractionBox["1", "2"], " ", "\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",",
+ RowBox[{"-", "1"}]}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"1", ",", "1", ",", "0"}], "}"}]}]], "[", "t", "]"}]}],
+ "+",
RowBox[{
SubscriptBox["\[Delta]", "1"], " ",
RowBox[{
@@ -9861,10 +9948,23 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]], "-",
RowBox[{
+ FractionBox["1", "2"], " ", "\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",",
+ RowBox[{"-", "1"}]}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"2", ",", "2", ",",
+ RowBox[{"-", "2"}]}], "}"}]}]], "[", "t", "]"}]}], "-",
+ RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -10020,9 +10120,22 @@ Cell[BoxData[
RowBox[{"-", "2"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "2", ",",
- RowBox[{"-", "2"}]}], "}"}]}]], "[", "t", "]"}]}], "+",
+ RowBox[{"-", "2"}]}], "}"}]}]], "[", "t", "]"}]}], "-",
+ RowBox[{"\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",",
+ RowBox[{"-", "2"}]}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"1", ",", "1", ",",
+ RowBox[{"-", "1"}]}], "}"}]}]], "[", "t", "]"}]}], "-",
RowBox[{
- SubscriptBox["\[Delta]", "1"], " ",
+ RowBox[{"(",
+ RowBox[{
+ FractionBox["\[CapitalOmega]L", "6"], "-",
+ SubscriptBox["\[Delta]", "1"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
RowBox[{
@@ -10060,10 +10173,22 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"0", ",", "1", ",", "1"}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",", "1"}], "}"}]}]], "[", "t", "]"}]}],
+ "+",
+ RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox["\[CapitalOmega]L", "2"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
@@ -10228,7 +10353,8 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox["\[CapitalOmega]L", "2"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
@@ -10503,9 +10629,13 @@ Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",", "2", ",", "0"}], "}"}]}]], "[", "t", "]"}]}],
RowBox[{"4", " ",
- SqrtBox["6"]}]], "+",
+ SqrtBox["6"]}]], "-",
RowBox[{
- SubscriptBox["\[Delta]", "1"], " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "6"]}], "-",
+ SubscriptBox["\[Delta]", "1"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
RowBox[{
@@ -10626,7 +10756,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -10680,10 +10812,21 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{
+ FractionBox["1", "2"], " ", "\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"0", ",", "1", ",",
+ RowBox[{"-", "1"}]}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",",
+ RowBox[{"-", "1"}]}], "}"}]}]], "[", "t", "]"}]}], "+",
+ RowBox[{
RowBox[{"(",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox["\[CapitalOmega]L", "2"], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
@@ -10943,8 +11086,8 @@ Cell[BoxData[
RowBox[{
RowBox[{"(",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox["\[CapitalOmega]L", "2"], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
@@ -11047,8 +11190,8 @@ Cell[BoxData[
RowBox[{
RowBox[{"(",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox["\[CapitalOmega]L", "2"], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
@@ -11063,7 +11206,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -11157,9 +11302,13 @@ Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",", "2", ",", "0"}], "}"}]}]], "[", "t", "]"}]}],
RowBox[{"4", " ",
- SqrtBox["6"]}]], "-",
+ SqrtBox["6"]}]], "+",
RowBox[{
- SubscriptBox["\[Delta]", "1"], " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "6"]}], "-",
+ SubscriptBox["\[Delta]", "1"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
RowBox[{
@@ -11247,9 +11396,13 @@ Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",", "1", ",", "0"}], "}"}]}]], "[", "t", "]"}]}],
RowBox[{"4", " ",
- SqrtBox["6"]}]], "-",
+ SqrtBox["6"]}]], "+",
RowBox[{
- SubscriptBox["\[Delta]", "1"], " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "6"]}], "-",
+ SubscriptBox["\[Delta]", "1"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
RowBox[{
@@ -11446,9 +11599,13 @@ Cell[BoxData[
RowBox[{"1", ",", "1", ",", "1"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "2", ",", "0"}], "}"}]}]], "[", "t", "]"}]}],
- "-",
+ "+",
RowBox[{
- SubscriptBox["\[Delta]", "1"], " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "6"]}], "-",
+ SubscriptBox["\[Delta]", "1"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
RowBox[{
@@ -11461,7 +11618,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -11533,7 +11692,17 @@ Cell[BoxData[
RowBox[{"0", ",", "2", ",",
RowBox[{"-", "1"}]}], "}"}]}]], "[", "t", "]"}]}],
RowBox[{"4", " ",
- SqrtBox["2"]}]], "-",
+ SqrtBox["2"]}]], "+",
+ RowBox[{
+ FractionBox["1", "2"], " ", "\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"1", ",", "1", ",", "0"}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",",
+ RowBox[{"-", "1"}]}], "}"}]}]], "[", "t", "]"}]}], "-",
RowBox[{
SubscriptBox["\[Delta]", "1"], " ",
RowBox[{
@@ -11643,8 +11812,8 @@ Cell[BoxData[
RowBox[{
RowBox[{"(",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox["\[CapitalOmega]L", "2"], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
@@ -11854,7 +12023,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -11913,9 +12084,22 @@ Cell[BoxData[
RowBox[{"-", "2"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "2", ",",
- RowBox[{"-", "2"}]}], "}"}]}]], "[", "t", "]"}]}], "-",
+ RowBox[{"-", "2"}]}], "}"}]}]], "[", "t", "]"}]}], "+",
+ RowBox[{"\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"1", ",", "1", ",",
+ RowBox[{"-", "1"}]}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",",
+ RowBox[{"-", "2"}]}], "}"}]}]], "[", "t", "]"}]}], "+",
RowBox[{
- SubscriptBox["\[Delta]", "1"], " ",
+ RowBox[{"(",
+ RowBox[{
+ FractionBox["\[CapitalOmega]L", "6"], "-",
+ SubscriptBox["\[Delta]", "1"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
RowBox[{
@@ -12027,12 +12211,22 @@ Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",", "2", ",", "2"}], "}"}]}]], "[", "t", "]"}]}],
RowBox[{"4", " ",
- SqrtBox["6"]}]]}], "+",
+ SqrtBox["6"]}]]}], "-",
+ RowBox[{"\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"2", ",", "2", ",", "1"}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",", "2"}], "}"}]}]], "[", "t", "]"}]}],
+ "+",
RowBox[{
RowBox[{"(",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -12152,6 +12346,16 @@ Cell[BoxData[
RowBox[{"0", ",", "2", ",", "1"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "2", ",", "1"}], "}"}]}]], "[", "t", "]"}]}],
+ "-",
+ RowBox[{
+ FractionBox["1", "2"], " ", "\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"2", ",", "2", ",", "0"}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",", "1"}], "}"}]}]], "[", "t", "]"}]}],
"+",
RowBox[{
RowBox[{"(",
@@ -12238,7 +12442,8 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox["\[CapitalOmega]L", "2"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
@@ -12399,7 +12604,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -12512,7 +12719,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -12635,9 +12844,13 @@ Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",", "2", ",", "0"}], "}"}]}]], "[", "t", "]"}]}],
RowBox[{"4", " ",
- SqrtBox["6"]}]], "+",
+ SqrtBox["6"]}]], "-",
RowBox[{
- SubscriptBox["\[Delta]", "1"], " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["\[CapitalOmega]L", "6"]}], "-",
+ SubscriptBox["\[Delta]", "1"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
RowBox[{
@@ -12651,7 +12864,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -12792,10 +13007,23 @@ Cell[BoxData[
RowBox[{"4", " ",
SqrtBox["6"]}]], "+",
RowBox[{
+ FractionBox["1", "2"], " ", "\[CapitalOmega]L", " ",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"2", ",", "2", ",",
+ RowBox[{"-", "2"}]}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"0", ",", "2", ",",
+ RowBox[{"-", "1"}]}], "}"}]}]], "[", "t", "]"}]}], "+",
+ RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -12900,8 +13128,8 @@ Cell[BoxData[
RowBox[{
RowBox[{"(",
RowBox[{
- RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox["\[CapitalOmega]L", "2"], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"]}], ")"}], " ",
RowBox[{
SubscriptBox["\[Rho]",
@@ -12916,7 +13144,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -13062,7 +13292,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- SubscriptBox["\[Delta]", "1"]}], "+",
+ FractionBox[
+ RowBox[{"4", " ", "\[CapitalOmega]L"}], "3"]}], "-",
+ SubscriptBox["\[Delta]", "1"], "+",
SubscriptBox["\[Delta]", "2"], "-",
SubscriptBox["\[Delta]", "3"]}], ")"}], " ",
RowBox[{
@@ -13389,7 +13621,8 @@ Cell[BoxData[{
SubscriptBox["\[Delta]", "1"], "\[Rule]", "0"}], ",",
RowBox[{
SubscriptBox["\[Delta]", "3"], "\[Rule]", "0"}], ",",
- RowBox[{"\[Gamma]t", "\[Rule]", ".0001"}]}], "}"}]}],
+ RowBox[{"\[Gamma]t", "\[Rule]", ".0001"}], ",",
+ RowBox[{"\[CapitalOmega]L", "\[Rule]", "0"}]}], "}"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{"linabs", "=",
RowBox[{
@@ -13433,9 +13666,7 @@ Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"-",
- RowBox[{
- FractionBox["1",
- SubscriptBox["\[CapitalOmega]", "2"]],
+ FractionBox[
RowBox[{"500.6135764115296`", " ",
RowBox[{"(",
RowBox[{
@@ -13497,10 +13728,9 @@ Cell[BoxData[
RowBox[{"0", ",", "2", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"1", ",", "1", ",", "1"}], "}"}]}]], "[", "t", "]"}],
- "]"}]}]}], ")"}]}]}]}], ",",
- RowBox[{
- FractionBox["1",
- SubscriptBox["\[CapitalOmega]", "2"]],
+ "]"}]}]}], ")"}]}],
+ SubscriptBox["\[CapitalOmega]", "2"]]}], ",",
+ FractionBox[
RowBox[{"250.3067882057648`", " ",
RowBox[{"(",
RowBox[{
@@ -13562,7 +13792,8 @@ Cell[BoxData[
RowBox[{"0", ",", "2", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"1", ",", "1", ",", "1"}], "}"}]}]], "[", "t", "]"}],
- "]"}]}]}], ")"}]}]}]}], "}"}]], "Output"]
+ "]"}]}]}], ")"}]}],
+ SubscriptBox["\[CapitalOmega]", "2"]]}], "}"}]], "Output"]
}, Open ]]
}, Open ]],
@@ -13606,7 +13837,8 @@ Cell[BoxData[{
SubscriptBox["\[Delta]", "1"], "\[Rule]", "0"}], ",",
RowBox[{
SubscriptBox["\[Delta]", "3"], "\[Rule]", "0"}], ",",
- RowBox[{"\[Gamma]t", "\[Rule]", ".1"}]}], "}"}]}],
+ RowBox[{"\[Gamma]t", "\[Rule]", ".1"}], ",",
+ RowBox[{"\[CapitalOmega]L", "\[Rule]", "0."}]}], "}"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"steadyeqs1", "=",
@@ -14016,6 +14248,7 @@ Cell[BoxData[{
RowBox[{"r", "[",
RowBox[{"s1", ",", "s2"}], "]"}]}], ",",
RowBox[{"\[Gamma]t", "\[Rule]", "gt"}], ",",
+ RowBox[{"\[CapitalOmega]L", "\[Rule]", "WL"}], ",",
RowBox[{
SubscriptBox["\[Delta]", "i_"], "\[Rule]",
RowBox[{"d", "[", "i", "]"}]}], ",",
@@ -14117,8 +14350,7 @@ Cell[BoxData[{
"\[RuleDelayed]",
RowBox[{"\"\<Ec\>\"", "<>", "a"}]}], ",",
RowBox[{
- RowBox[{"\"\<d(\>\"", "~~", "a_", "~~", "\"\<)\>\""}], "\[RuleDelayed]",
-
+ RowBox[{"\"\<d(\>\"", "~~", "a_", "~~", "\"\<)\>\""}], "\[RuleDelayed]",
RowBox[{"\"\<d\>\"", "<>", "a"}]}], ",",
RowBox[{
RowBox[{"\"\<k(\>\"", "~~", "a_", "~~", "\"\<)\>\""}], "\[RuleDelayed]",
@@ -14184,20 +14416,23 @@ Cell[BoxData[
RowBox[{"i", " ",
RowBox[{"(",
RowBox[{
+ RowBox[{"WL", " ",
+ RowBox[{"r", "[",
+ RowBox[{"\<\"01\"\>", ",", "\<\"13\"\>"}], "]"}]}], "-",
RowBox[{
- RowBox[{"-",
- RowBox[{"(",
- RowBox[{
- RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
- RowBox[{"d", "[", "2", "]"}], "-",
- RowBox[{"d", "[", "3", "]"}], "-",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "1", "]"}]}], "+",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "2", "]"}]}], "-",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "3", "]"}]}]}], ")"}]}], " ",
+ RowBox[{"(",
+ RowBox[{
+ FractionBox[
+ RowBox[{"2", " ", "WL"}], "3"], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
+ RowBox[{"d", "[", "2", "]"}], "-",
+ RowBox[{"d", "[", "3", "]"}], "-",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "1", "]"}]}], "+",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "2", "]"}]}], "-",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "3", "]"}]}]}], ")"}], " ",
RowBox[{"r", "[",
RowBox[{"\<\"01\"\>", ",", "\<\"13\"\>"}], "]"}]}], "+",
FractionBox[
@@ -14275,19 +14510,22 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{
- RowBox[{"-",
- RowBox[{"(",
- RowBox[{
- RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
- RowBox[{"d", "[", "2", "]"}], "-",
- RowBox[{"d", "[", "3", "]"}], "-",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "1", "]"}]}], "+",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "2", "]"}]}], "-",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "3", "]"}]}]}], ")"}]}], " ",
+ FractionBox["1", "2"], " ", "WL", " ",
+ RowBox[{"r", "[",
+ RowBox[{"\<\"02\"\>", ",", "\<\"14\"\>"}], "]"}]}], "-",
+ RowBox[{
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ RowBox[{"d", "[", "1", "]"}]}], "+",
+ RowBox[{"d", "[", "2", "]"}], "-",
+ RowBox[{"d", "[", "3", "]"}], "-",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "1", "]"}]}], "+",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "2", "]"}]}], "-",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "3", "]"}]}]}], ")"}], " ",
RowBox[{"r", "[",
RowBox[{"\<\"02\"\>", ",", "\<\"14\"\>"}], "]"}]}], "-",
FractionBox[
@@ -14386,11 +14624,14 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{
- RowBox[{"(",
- RowBox[{
- RowBox[{"d", "[", "1", "]"}], "+",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "1", "]"}]}]}], ")"}], " ",
+ RowBox[{"-",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["WL", "6"]}], "-",
+ RowBox[{"d", "[", "1", "]"}], "-",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "1", "]"}]}]}], ")"}]}], " ",
RowBox[{"r", "[",
RowBox[{"\<\"03\"\>", ",", "\<\"09\"\>"}], "]"}]}], "+",
FractionBox[
@@ -14442,7 +14683,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "WL"}], "3"]}], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
RowBox[{"d", "[", "2", "]"}], "-",
RowBox[{"d", "[", "3", "]"}], "-",
RowBox[{"v", " ",
@@ -14556,6 +14799,11 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "WL", " ",
+ RowBox[{"r", "[",
+ RowBox[{"\<\"04\"\>", ",", "\<\"10\"\>"}], "]"}]}], "+",
+ RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"d", "[", "1", "]"}], "+",
@@ -14610,18 +14858,24 @@ Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"-",
- RowBox[{"(",
- RowBox[{
- RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
- RowBox[{"d", "[", "2", "]"}], "-",
- RowBox[{"d", "[", "3", "]"}], "-",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "1", "]"}]}], "+",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "2", "]"}]}], "-",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "3", "]"}]}]}], ")"}]}], " ",
+ FractionBox["1", "2"]}], " ", "WL", " ",
+ RowBox[{"r", "[",
+ RowBox[{"\<\"04\"\>", ",", "\<\"16\"\>"}], "]"}]}], "-",
+ RowBox[{
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox[
+ RowBox[{"4", " ", "WL"}], "3"]}], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
+ RowBox[{"d", "[", "2", "]"}], "-",
+ RowBox[{"d", "[", "3", "]"}], "-",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "1", "]"}]}], "+",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "2", "]"}]}], "-",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "3", "]"}]}]}], ")"}], " ",
RowBox[{"r", "[",
RowBox[{"\<\"04\"\>", ",", "\<\"16\"\>"}], "]"}]}], "-",
FractionBox[
@@ -14706,9 +14960,14 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{
+ RowBox[{"-", "WL"}], " ",
+ RowBox[{"r", "[",
+ RowBox[{"\<\"05\"\>", ",", "\<\"11\"\>"}], "]"}]}], "-",
+ RowBox[{
RowBox[{"(",
RowBox[{
- RowBox[{"d", "[", "1", "]"}], "+",
+ FractionBox["WL", "6"], "-",
+ RowBox[{"d", "[", "1", "]"}], "-",
RowBox[{"v", " ",
RowBox[{"k", "[", "1", "]"}]}]}], ")"}], " ",
RowBox[{"r", "[",
@@ -14735,10 +14994,16 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "WL", " ",
+ RowBox[{"r", "[",
+ RowBox[{"\<\"06\"\>", ",", "\<\"02\"\>"}], "]"}]}], "+",
+ RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
+ FractionBox["WL", "2"]}], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
RowBox[{"d", "[", "2", "]"}], "-",
RowBox[{"v", " ",
RowBox[{"k", "[", "1", "]"}]}], "+",
@@ -14826,7 +15091,8 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
+ FractionBox["WL", "2"]}], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
RowBox[{"d", "[", "2", "]"}], "-",
RowBox[{"v", " ",
RowBox[{"k", "[", "1", "]"}]}], "+",
@@ -14992,11 +15258,14 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{
- RowBox[{"(",
- RowBox[{
- RowBox[{"d", "[", "1", "]"}], "+",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "1", "]"}]}]}], ")"}], " ",
+ RowBox[{"-",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["WL", "6"]}], "-",
+ RowBox[{"d", "[", "1", "]"}], "-",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "1", "]"}]}]}], ")"}]}], " ",
RowBox[{"r", "[",
RowBox[{"\<\"07\"\>", ",", "\<\"09\"\>"}], "]"}]}], "+",
RowBox[{
@@ -15071,7 +15340,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "WL"}], "3"]}], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
RowBox[{"d", "[", "2", "]"}], "-",
RowBox[{"d", "[", "3", "]"}], "-",
RowBox[{"v", " ",
@@ -15115,10 +15386,14 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{
+ FractionBox["1", "2"], " ", "WL", " ",
+ RowBox[{"r", "[",
+ RowBox[{"\<\"08\"\>", ",", "\<\"04\"\>"}], "]"}]}], "+",
+ RowBox[{
RowBox[{"(",
RowBox[{
- RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
+ FractionBox["WL", "2"], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
RowBox[{"d", "[", "2", "]"}], "-",
RowBox[{"v", " ",
RowBox[{"k", "[", "1", "]"}]}], "+",
@@ -15243,8 +15518,8 @@ Cell[BoxData[
RowBox[{
RowBox[{"(",
RowBox[{
- RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
+ FractionBox["WL", "2"], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
RowBox[{"d", "[", "2", "]"}], "-",
RowBox[{"v", " ",
RowBox[{"k", "[", "1", "]"}]}], "+",
@@ -15300,8 +15575,8 @@ Cell[BoxData[
RowBox[{
RowBox[{"(",
RowBox[{
- RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
+ FractionBox["WL", "2"], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
RowBox[{"d", "[", "2", "]"}], "-",
RowBox[{"v", " ",
RowBox[{"k", "[", "1", "]"}]}], "+",
@@ -15313,7 +15588,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
+ FractionBox[
+ RowBox[{"4", " ", "WL"}], "3"]}], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
RowBox[{"d", "[", "2", "]"}], "-",
RowBox[{"d", "[", "3", "]"}], "-",
RowBox[{"v", " ",
@@ -15412,19 +15689,22 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{
- RowBox[{"-",
- RowBox[{"(",
- RowBox[{
- RowBox[{"d", "[", "1", "]"}], "+",
- RowBox[{"v", " ",
- RowBox[{"k", "[", "1", "]"}]}]}], ")"}]}], " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["WL", "6"]}], "-",
+ RowBox[{"d", "[", "1", "]"}], "-",
+ RowBox[{"v", " ",
+ RowBox[{"k", "[", "1", "]"}]}]}], ")"}], " ",
RowBox[{"r", "[",
RowBox[{"\<\"09\"\>", ",", "\<\"15\"\>"}], "]"}]}], "-",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
+ FractionBox[
+ RowBox[{"2", " ", "WL"}], "3"]}], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
RowBox[{"d", "[", "2", "]"}], "-",
RowBox[{"d", "[", "3", "]"}], "-",
RowBox[{"v", " ",
@@ -15533,7 +15813,9 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{
RowBox[{"-",
- RowBox[{"d", "[", "1", "]"}]}], "+",
+ FractionBox[
+ RowBox[{"4", " ", "WL"}], "3"]}], "-",
+ RowBox[{"d", "[", "1", "]"}], "+",
RowBox[{"d", "[", "2", "]"}], "-",
RowBox[{"d", "[", "3", "]"}], "-",
RowBox[{"v", " ",
@@ -15732,69 +16014,73 @@ Cell[BoxData[
Cell[BoxData["\<\"rt6 = 2.449489742783178;\\nrt3 = 1.7320508075688772;\\nrt2 \
= 1.4142135623730951;\\ndr0101_dt = gt/8. - gt*r0101 + (g1*r0909)/2. + \
(g2*r1313)/6. - i*((r0113*E4)/(4.*rt6) - (r1301*Ec4)/(4.*rt6));\\ndr0113_dt = \
-(-(gt*r0113) - (gt + g2)*r0113)/2. - i*(-((-d1 + d2 - d3 - v*k1 + v*k2 - \
-v*k3)*r0113) + (r0101*Ec4)/(4.*rt6) - (r1313*Ec4)/(4.*rt6));\\ndr0202_dt = \
-gt/8. - gt*r0202 + (g1*r0909)/4. + (g1*r1010)/4. + (g2*r1313)/12. + \
-(g2*r1414)/4. - i*((r0214*E4)/8. - (r1402*Ec4)/8.);\\ndr0214_dt = \
-(-(gt*r0214) - (gt + g2)*r0214)/2. - i*(-((-d1 + d2 - d3 - v*k1 + v*k2 - \
-v*k3)*r0214) - (r0206*Ec3)/(8.*rt3) + (r0202*Ec4)/8. - \
-(r1414*Ec4)/8.);\\ndr0303_dt = gt/8. - gt*r0303 + (g1*r0909)/12. + \
-(g1*r1010)/3. + (g1*r1111)/12. + (g2*r1313)/4. + (g2*r1515)/4. - \
-i*((r0309*E1)/(4.*rt6) + (r0315*E4)/8. - (r0903*Ec1)/(4.*rt6) - \
-(r1503*Ec4)/8.);\\ndr0309_dt = (-(gt*r0309) - (gt + g1)*r0309)/2. - i*((d1 + \
-v*k1)*r0309 + (r0303*Ec1)/(4.*rt6) - (r0909*Ec1)/(4.*rt6) - \
-(r0307*Ec2)/(4.*rt6) - (r1509*Ec4)/8.);\\ndr0315_dt = (-(gt*r0315) - (gt + \
-g2)*r0315)/2. - i*(-((-d1 + d2 - d3 - v*k1 + v*k2 - v*k3)*r0315) - \
-(r0915*Ec1)/(4.*rt6) - (r0307*Ec3)/8. + (r0303*Ec4)/8. - \
-(r1515*Ec4)/8.);\\ndr0404_dt = gt/8. - gt*r0404 + (g1*r1010)/4. + \
+(-(gt*r0113) - (gt + g2)*r0113)/2. - i*(WL*r0113 - ((2*WL)/3. - d1 + d2 - d3 \
+- v*k1 + v*k2 - v*k3)*r0113 + (r0101*Ec4)/(4.*rt6) - \
+(r1313*Ec4)/(4.*rt6));\\ndr0202_dt = gt/8. - gt*r0202 + (g1*r0909)/4. + \
+(g1*r1010)/4. + (g2*r1313)/12. + (g2*r1414)/4. - i*((r0214*E4)/8. - \
+(r1402*Ec4)/8.);\\ndr0214_dt = (-(gt*r0214) - (gt + g2)*r0214)/2. - \
+i*((WL*r0214)/2. - (-d1 + d2 - d3 - v*k1 + v*k2 - v*k3)*r0214 - \
+(r0206*Ec3)/(8.*rt3) + (r0202*Ec4)/8. - (r1414*Ec4)/8.);\\ndr0303_dt = gt/8. \
+- gt*r0303 + (g1*r0909)/12. + (g1*r1010)/3. + (g1*r1111)/12. + (g2*r1313)/4. \
++ (g2*r1515)/4. - i*((r0309*E1)/(4.*rt6) + (r0315*E4)/8. - \
+(r0903*Ec1)/(4.*rt6) - (r1503*Ec4)/8.);\\ndr0309_dt = (-(gt*r0309) - (gt + \
+g1)*r0309)/2. - i*(-((-WL/6. - d1 - v*k1)*r0309) + (r0303*Ec1)/(4.*rt6) - \
+(r0909*Ec1)/(4.*rt6) - (r0307*Ec2)/(4.*rt6) - (r1509*Ec4)/8.);\\ndr0315_dt = \
+(-(gt*r0315) - (gt + g2)*r0315)/2. - i*(-(((-2*WL)/3. - d1 + d2 - d3 - v*k1 + \
+v*k2 - v*k3)*r0315) - (r0915*Ec1)/(4.*rt6) - (r0307*Ec3)/8. + (r0303*Ec4)/8. \
+- (r1515*Ec4)/8.);\\ndr0404_dt = gt/8. - gt*r0404 + (g1*r1010)/4. + \
(g1*r1111)/4. + (g2*r1414)/4. + (g2*r1515)/12. + (g2*r1616)/6. - \
i*((r0410*E1)/(4.*rt2) + (r0416*E4)/(4.*rt6) - (r1004*Ec1)/(4.*rt2) - \
(r1604*Ec4)/(4.*rt6));\\ndr0410_dt = (-(gt*r0410) - (gt + g1)*r0410)/2. - \
-i*((d1 + v*k1)*r0410 + (r0404*Ec1)/(4.*rt2) - (r1010*Ec1)/(4.*rt2) - \
-(r0408*Ec2)/(4.*rt6) - (r1610*Ec4)/(4.*rt6));\\ndr0416_dt = (-(gt*r0416) - \
-(gt + g2)*r0416)/2. - i*(-((-d1 + d2 - d3 - v*k1 + v*k2 - v*k3)*r0416) - \
+i*(-(WL*r0410)/2. + (d1 + v*k1)*r0410 + (r0404*Ec1)/(4.*rt2) - \
+(r1010*Ec1)/(4.*rt2) - (r0408*Ec2)/(4.*rt6) - \
+(r1610*Ec4)/(4.*rt6));\\ndr0416_dt = (-(gt*r0416) - (gt + g2)*r0416)/2. - \
+i*(-(WL*r0416)/2. - ((-4*WL)/3. - d1 + d2 - d3 - v*k1 + v*k2 - v*k3)*r0416 - \
(r1016*Ec1)/(4.*rt2) - (r0408*Ec3)/(4.*rt2) + (r0404*Ec4)/(4.*rt6) - \
(r1616*Ec4)/(4.*rt6));\\ndr0505_dt = gt/8. - gt*r0505 + (g1*r1111)/2. + \
(g2*r1515)/6. + (g2*r1616)/3. - i*((r0511*E1)/4. - \
-(r1105*Ec1)/4.);\\ndr0511_dt = (-(gt*r0511) - (gt + g1)*r0511)/2. - i*((d1 + \
-v*k1)*r0511 + (r0505*Ec1)/4. - (r1111*Ec1)/4.);\\ndr0602_dt = -(gt*r0602) - \
-i*((-d1 + d2 - v*k1 + v*k2)*r0602 + (r0614*E4)/8. + \
+(r1105*Ec1)/4.);\\ndr0511_dt = (-(gt*r0511) - (gt + g1)*r0511)/2. - \
+i*(-(WL*r0511) - (WL/6. - d1 - v*k1)*r0511 + (r0505*Ec1)/4. - \
+(r1111*Ec1)/4.);\\ndr0602_dt = -(gt*r0602) - i*(-(WL*r0602)/2. + (-WL/2. - d1 \
++ d2 - v*k1 + v*k2)*r0602 + (r0614*E4)/8. + \
(r1402*Ec3)/(8.*rt3));\\ndr0606_dt = gt/8. - gt*r0606 + (g1*r0909)/12. + \
(g1*r1010)/12. + (g2*r1313)/4. + (g2*r1414)/12. - i*(-(r0614*E3)/(8.*rt3) + \
(r1406*Ec3)/(8.*rt3));\\ndr0614_dt = (-(gt*r0614) - (gt + g2)*r0614)/2. - \
-i*((-d1 + d2 - v*k1 + v*k2)*r0614 - (-d1 + d2 - d3 - v*k1 + v*k2 - \
+i*((-WL/2. - d1 + d2 - v*k1 + v*k2)*r0614 - (-d1 + d2 - d3 - v*k1 + v*k2 - \
v*k3)*r0614 - (r0606*Ec3)/(8.*rt3) + (r1414*Ec3)/(8.*rt3) + \
(r0602*Ec4)/8.);\\ndr0703_dt = -(gt*r0703) - i*((-d1 + d2 - v*k1 + \
v*k2)*r0703 + (r0709*E1)/(4.*rt6) + (r0715*E4)/8. + (r0903*Ec2)/(4.*rt6) + \
(r1503*Ec3)/8.);\\ndr0707_dt = gt/8. - gt*r0707 + (g1*r0909)/12. + \
(g1*r1111)/12. + (g2*r1313)/4. + (g2*r1414)/3. + (g2*r1515)/4. - \
i*(-(r0709*E2)/(4.*rt6) - (r0715*E3)/8. + (r0907*Ec2)/(4.*rt6) + \
-(r1507*Ec3)/8.);\\ndr0709_dt = (-(gt*r0709) - (gt + g1)*r0709)/2. - i*((d1 + \
-v*k1)*r0709 + (-d1 + d2 - v*k1 + v*k2)*r0709 + (r0703*Ec1)/(4.*rt6) - \
-(r0707*Ec2)/(4.*rt6) + (r0909*Ec2)/(4.*rt6) + (r1509*Ec3)/8.);\\ndr0715_dt = \
-(-(gt*r0715) - (gt + g2)*r0715)/2. - i*((-d1 + d2 - v*k1 + v*k2)*r0715 - (-d1 \
-+ d2 - d3 - v*k1 + v*k2 - v*k3)*r0715 + (r0915*Ec2)/(4.*rt6) - (r0707*Ec3)/8. \
-+ (r1515*Ec3)/8. + (r0703*Ec4)/8.);\\ndr0804_dt = -(gt*r0804) - i*((-d1 + d2 \
-- v*k1 + v*k2)*r0804 + (r0810*E1)/(4.*rt2) + (r0816*E4)/(4.*rt6) + \
+(r1507*Ec3)/8.);\\ndr0709_dt = (-(gt*r0709) - (gt + g1)*r0709)/2. - \
+i*(-((-WL/6. - d1 - v*k1)*r0709) + (-d1 + d2 - v*k1 + v*k2)*r0709 + \
+(r0703*Ec1)/(4.*rt6) - (r0707*Ec2)/(4.*rt6) + (r0909*Ec2)/(4.*rt6) + \
+(r1509*Ec3)/8.);\\ndr0715_dt = (-(gt*r0715) - (gt + g2)*r0715)/2. - i*((-d1 + \
+d2 - v*k1 + v*k2)*r0715 - ((-2*WL)/3. - d1 + d2 - d3 - v*k1 + v*k2 - \
+v*k3)*r0715 + (r0915*Ec2)/(4.*rt6) - (r0707*Ec3)/8. + (r1515*Ec3)/8. + \
+(r0703*Ec4)/8.);\\ndr0804_dt = -(gt*r0804) - i*((WL*r0804)/2. + (WL/2. - d1 + \
+d2 - v*k1 + v*k2)*r0804 + (r0810*E1)/(4.*rt2) + (r0816*E4)/(4.*rt6) + \
(r1004*Ec2)/(4.*rt6) + (r1604*Ec3)/(4.*rt2));\\ndr0808_dt = gt/8. - gt*r0808 \
+ (g1*r1010)/12. + (g1*r1111)/12. + (g2*r1414)/12. + (g2*r1515)/4. + \
(g2*r1616)/2. - i*(-(r0810*E2)/(4.*rt6) - (r0816*E3)/(4.*rt2) + \
(r1008*Ec2)/(4.*rt6) + (r1608*Ec3)/(4.*rt2));\\ndr0810_dt = (-(gt*r0810) - \
-(gt + g1)*r0810)/2. - i*((d1 + v*k1)*r0810 + (-d1 + d2 - v*k1 + v*k2)*r0810 + \
-(r0804*Ec1)/(4.*rt2) - (r0808*Ec2)/(4.*rt6) + (r1010*Ec2)/(4.*rt6) + \
-(r1610*Ec3)/(4.*rt2));\\ndr0816_dt = (-(gt*r0816) - (gt + g2)*r0816)/2. - \
-i*((-d1 + d2 - v*k1 + v*k2)*r0816 - (-d1 + d2 - d3 - v*k1 + v*k2 - \
-v*k3)*r0816 + (r1016*Ec2)/(4.*rt6) - (r0808*Ec3)/(4.*rt2) + \
-(r1616*Ec3)/(4.*rt2) + (r0804*Ec4)/(4.*rt6));\\ndr0909_dt = -((gt + \
-g1)*r0909) - i*(-(r0309*E1)/(4.*rt6) + (r0709*E2)/(4.*rt6) + \
-(r0903*Ec1)/(4.*rt6) - (r0907*Ec2)/(4.*rt6));\\ndr0915_dt = (-((gt + \
-g1)*r0915) - (gt + g2)*r0915)/2. - i*(-((d1 + v*k1)*r0915) - (-d1 + d2 - d3 - \
+(gt + g1)*r0810)/2. - i*((d1 + v*k1)*r0810 + (WL/2. - d1 + d2 - v*k1 + \
+v*k2)*r0810 + (r0804*Ec1)/(4.*rt2) - (r0808*Ec2)/(4.*rt6) + \
+(r1010*Ec2)/(4.*rt6) + (r1610*Ec3)/(4.*rt2));\\ndr0816_dt = (-(gt*r0816) - \
+(gt + g2)*r0816)/2. - i*((WL/2. - d1 + d2 - v*k1 + v*k2)*r0816 - ((-4*WL)/3. \
+- d1 + d2 - d3 - v*k1 + v*k2 - v*k3)*r0816 + (r1016*Ec2)/(4.*rt6) - \
+(r0808*Ec3)/(4.*rt2) + (r1616*Ec3)/(4.*rt2) + \
+(r0804*Ec4)/(4.*rt6));\\ndr0909_dt = -((gt + g1)*r0909) - \
+i*(-(r0309*E1)/(4.*rt6) + (r0709*E2)/(4.*rt6) + (r0903*Ec1)/(4.*rt6) - \
+(r0907*Ec2)/(4.*rt6));\\ndr0915_dt = (-((gt + g1)*r0915) - (gt + \
+g2)*r0915)/2. - i*((-WL/6. - d1 - v*k1)*r0915 - ((-2*WL)/3. - d1 + d2 - d3 - \
v*k1 + v*k2 - v*k3)*r0915 - (r0315*E1)/(4.*rt6) + (r0715*E2)/(4.*rt6) - \
(r0907*Ec3)/8. + (r0903*Ec4)/8.);\\ndr1010_dt = -((gt + g1)*r1010) - \
i*(-(r0410*E1)/(4.*rt2) + (r0810*E2)/(4.*rt6) + (r1004*Ec1)/(4.*rt2) - \
(r1008*Ec2)/(4.*rt6));\\ndr1016_dt = (-((gt + g1)*r1016) - (gt + \
-g2)*r1016)/2. - i*(-((d1 + v*k1)*r1016) - (-d1 + d2 - d3 - v*k1 + v*k2 - \
-v*k3)*r1016 - (r0416*E1)/(4.*rt2) + (r0816*E2)/(4.*rt6) - \
+g2)*r1016)/2. - i*(-((d1 + v*k1)*r1016) - ((-4*WL)/3. - d1 + d2 - d3 - v*k1 + \
+v*k2 - v*k3)*r1016 - (r0416*E1)/(4.*rt2) + (r0816*E2)/(4.*rt6) - \
(r1008*Ec3)/(4.*rt2) + (r1004*Ec4)/(4.*rt6));\\ndr1111_dt = -((gt + \
g1)*r1111) - i*(-(r0511*E1)/4. + (r1105*Ec1)/4.);\\ndr1313_dt = -((gt + \
g2)*r1313) - i*(-(r0113*E4)/(4.*rt6) + (r1301*Ec4)/(4.*rt6));\\ndr1414_dt = \
@@ -17369,8 +17655,8 @@ mathemathica_fwm\\\\RbPropEquations.txt\"\>"], "Output"]
}, Open ]]
}, Open ]]
},
-WindowSize->{978, 883},
-WindowMargins->{{60, Automatic}, {11, Automatic}},
+WindowSize->{1104, 893},
+WindowMargins->{{17, Automatic}, {Automatic, 8}},
ShowSelection->True,
FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (October 6, 2011)",
StyleDefinitions->"Default.nb"
@@ -17423,230 +17709,218 @@ Cell[15464, 466, 1150, 43, 29, "Text",
CellID->133602844],
Cell[CellGroupData[{
Cell[16639, 513, 325, 8, 31, "Input"],
-Cell[16967, 523, 478, 12, 50, "Output"]
+Cell[16967, 523, 478, 12, 30, "Output"]
}, Open ]],
Cell[17460, 538, 411, 10, 47, "Text"],
Cell[CellGroupData[{
Cell[17896, 552, 2561, 80, 92, "Input",
CellID->534530029],
-Cell[20460, 634, 4983, 149, 147, "Output"]
+Cell[20460, 634, 4983, 149, 101, "Output"]
}, Open ]],
Cell[25458, 786, 288, 6, 47, "Text",
CellID->462076121],
Cell[CellGroupData[{
-Cell[25771, 796, 621, 19, 31, "Input",
+Cell[25771, 796, 810, 24, 52, "Input",
CellID->494599775],
-Cell[26395, 817, 82830, 2374, 570, "Output"]
+Cell[26584, 822, 84081, 2408, 594, "Output"]
}, Open ]],
-Cell[109240, 3194, 256, 5, 47, "Text",
+Cell[110680, 3233, 256, 5, 47, "Text",
CellID->358620443],
Cell[CellGroupData[{
-Cell[109521, 3203, 643, 18, 52, "Input",
+Cell[110961, 3242, 698, 18, 52, "Input",
CellID->167259034],
-Cell[110167, 3223, 2891, 50, 221, "Output"]
+Cell[111662, 3262, 3230, 55, 230, "Output"]
}, Open ]],
-Cell[113073, 3276, 1513, 54, 68, "Text",
+Cell[114907, 3320, 1513, 54, 50, "Text",
CellID->577766068],
Cell[CellGroupData[{
-Cell[114611, 3334, 956, 31, 52, "Input"],
-Cell[115570, 3367, 997, 34, 30, "Output"],
-Cell[116570, 3403, 1196, 34, 30, "Output"]
-}, Open ]],
-Cell[117781, 3440, 102, 2, 29, "Text"],
-Cell[CellGroupData[{
-Cell[117908, 3446, 447, 12, 31, "Input"],
-Cell[118358, 3460, 7543, 205, 310, "Output"]
-}, Open ]],
-Cell[CellGroupData[{
-Cell[125938, 3670, 125, 3, 31, "Input"],
-Cell[126066, 3675, 2667, 84, 191, "Output"]
-}, Open ]],
-Cell[CellGroupData[{
-Cell[128770, 3764, 190, 5, 31, "Input"],
-Cell[128963, 3771, 527, 14, 47, "Output"]
+Cell[116445, 3378, 956, 31, 52, "Input"],
+Cell[117404, 3411, 997, 34, 30, "Output"],
+Cell[118404, 3447, 1196, 34, 30, "Output"]
}, Open ]],
+Cell[119615, 3484, 102, 2, 29, "Text"],
Cell[CellGroupData[{
-Cell[129527, 3790, 733, 21, 31, "Input"],
-Cell[130263, 3813, 261, 8, 47, "Output"]
+Cell[119742, 3490, 447, 12, 31, "Input"],
+Cell[120192, 3504, 7543, 205, 310, "Output"]
}, Open ]],
-Cell[130539, 3824, 902, 30, 47, "Text"],
+Cell[127750, 3712, 902, 30, 47, "Text"],
Cell[CellGroupData[{
-Cell[131466, 3858, 1142, 31, 52, "Input"],
-Cell[132611, 3891, 29869, 777, 569, "Output"]
+Cell[128677, 3746, 1142, 31, 52, "Input"],
+Cell[129822, 3779, 33911, 877, 570, "Output"]
}, Open ]],
-Cell[162495, 4671, 837, 25, 47, "Text"],
+Cell[163748, 4659, 837, 25, 47, "Text"],
Cell[CellGroupData[{
-Cell[163357, 4700, 765, 20, 38, "Input"],
-Cell[164125, 4722, 12475, 349, 569, "Output"]
+Cell[164610, 4688, 765, 20, 38, "Input"],
+Cell[165378, 4710, 13324, 371, 569, "Output"]
}, Open ]],
-Cell[176615, 5074, 141, 3, 29, "Text"],
+Cell[178717, 5084, 141, 3, 29, "Text"],
Cell[CellGroupData[{
-Cell[176781, 5081, 787, 22, 52, "Input"],
-Cell[177571, 5105, 261, 8, 47, "Output"]
+Cell[178883, 5091, 787, 22, 31, "Input"],
+Cell[179673, 5115, 261, 8, 47, "Output"]
}, Open ]],
Cell[CellGroupData[{
-Cell[177869, 5118, 787, 22, 52, "Input"],
-Cell[178659, 5142, 261, 8, 47, "Output"]
+Cell[179971, 5128, 787, 22, 31, "Input"],
+Cell[180761, 5152, 261, 8, 47, "Output"]
}, Open ]],
Cell[CellGroupData[{
-Cell[178957, 5155, 787, 22, 52, "Input"],
-Cell[179747, 5179, 355, 11, 47, "Output"]
+Cell[181059, 5165, 787, 22, 31, "Input"],
+Cell[181849, 5189, 355, 11, 47, "Output"]
}, Open ]],
-Cell[180117, 5193, 989, 34, 47, "Text"],
+Cell[182219, 5203, 989, 34, 47, "Text"],
Cell[CellGroupData[{
-Cell[181131, 5231, 877, 25, 47, "Input"],
-Cell[182011, 5258, 9961, 278, 559, "Output"]
+Cell[183233, 5241, 877, 25, 47, "Input"],
+Cell[184113, 5268, 10598, 294, 569, "Output"]
}, Open ]],
-Cell[191987, 5539, 92, 2, 29, "Text"],
+Cell[194726, 5565, 92, 2, 29, "Text"],
Cell[CellGroupData[{
-Cell[192104, 5545, 1047, 30, 52, "Input",
+Cell[194843, 5571, 1102, 30, 52, "Input",
CellID->102096636],
-Cell[193154, 5577, 1686, 32, 221, "Output"]
+Cell[195948, 5603, 2044, 35, 230, "Output"]
}, Open ]],
-Cell[194855, 5612, 588, 18, 29, "Text",
+Cell[198007, 5641, 588, 18, 29, "Text",
CellID->610306692],
Cell[CellGroupData[{
-Cell[195468, 5634, 272, 7, 31, "Input",
+Cell[198620, 5663, 272, 7, 31, "Input",
CellID->645617687],
-Cell[195743, 5643, 2716, 58, 266, "Output"]
+Cell[198895, 5672, 2716, 58, 266, "Output"]
}, Open ]],
-Cell[198474, 5704, 460, 13, 29, "Text",
+Cell[201626, 5733, 460, 13, 29, "Text",
CellID->854192725],
Cell[CellGroupData[{
-Cell[198959, 5721, 274, 7, 31, "Input",
+Cell[202111, 5750, 274, 7, 31, "Input",
CellID->465762594],
-Cell[199236, 5730, 49693, 1437, 426, "Output"]
+Cell[202388, 5759, 49693, 1437, 426, "Output"]
}, Open ]],
-Cell[248944, 7170, 69, 1, 29, "Text",
+Cell[252096, 7199, 69, 1, 29, "Text",
CellID->314466782],
Cell[CellGroupData[{
-Cell[249038, 7175, 197, 5, 31, "Input"],
-Cell[249238, 7182, 4464, 129, 126, "Output"]
+Cell[252190, 7204, 197, 5, 31, "Input"],
+Cell[252390, 7211, 4464, 129, 126, "Output"]
}, Open ]],
-Cell[253717, 7314, 232, 4, 29, "Text"],
+Cell[256869, 7343, 232, 4, 29, "Text"],
Cell[CellGroupData[{
-Cell[253974, 7322, 1661, 45, 119, "Input"],
-Cell[255638, 7369, 891, 21, 50, "Output"]
+Cell[257126, 7351, 1661, 45, 119, "Input"],
+Cell[258790, 7398, 891, 21, 50, "Output"]
}, Open ]],
-Cell[256544, 7393, 519, 13, 49, "Text"],
-Cell[257066, 7408, 1704, 44, 120, "Input"],
-Cell[258773, 7454, 105, 3, 29, "Text"],
+Cell[259696, 7422, 519, 13, 49, "Text"],
+Cell[260218, 7437, 1771, 45, 98, "Input"],
+Cell[261992, 7484, 105, 3, 29, "Text"],
Cell[CellGroupData[{
-Cell[258903, 7461, 213, 6, 31, "Input"],
-Cell[259119, 7469, 7417, 209, 240, "Output"]
+Cell[262122, 7491, 213, 6, 31, "Input"],
+Cell[262338, 7499, 7417, 209, 202, "Output"]
}, Open ]],
-Cell[266551, 7681, 93, 2, 29, "Text"],
-Cell[266647, 7685, 187, 6, 31, "Input"],
-Cell[266837, 7693, 80, 2, 29, "Text"],
-Cell[266920, 7697, 216, 8, 31, "Input"],
-Cell[267139, 7707, 88, 2, 29, "Text"],
+Cell[269770, 7711, 93, 2, 29, "Text"],
+Cell[269866, 7715, 187, 6, 31, "Input"],
+Cell[270056, 7723, 80, 2, 29, "Text"],
+Cell[270139, 7727, 216, 8, 31, "Input"],
+Cell[270358, 7737, 88, 2, 29, "Text"],
Cell[CellGroupData[{
-Cell[267252, 7713, 163, 5, 31, "Input"],
-Cell[267418, 7720, 11989, 385, 172, "Output"]
+Cell[270471, 7743, 163, 5, 31, "Input"],
+Cell[270637, 7750, 11989, 385, 152, "Output"]
}, Open ]],
-Cell[279422, 8108, 85, 2, 29, "Text"],
+Cell[282641, 8138, 85, 2, 29, "Text"],
Cell[CellGroupData[{
-Cell[279532, 8114, 190, 6, 31, "Input"],
-Cell[279725, 8122, 4464, 129, 126, "Output"]
+Cell[282751, 8144, 190, 6, 31, "Input"],
+Cell[282944, 8152, 4464, 129, 126, "Output"]
}, Open ]],
-Cell[284204, 8254, 71, 2, 29, "Text"],
+Cell[287423, 8284, 71, 2, 29, "Text"],
Cell[CellGroupData[{
-Cell[284300, 8260, 290, 8, 31, "Input"],
-Cell[284593, 8270, 14287, 488, 297, "Output"]
+Cell[287519, 8290, 290, 8, 31, "Input"],
+Cell[287812, 8300, 14287, 488, 277, "Output"]
}, Open ]],
-Cell[298895, 8761, 323, 9, 29, "Text"],
+Cell[302114, 8791, 297, 8, 29, "Text"],
Cell[CellGroupData[{
-Cell[299243, 8774, 534, 18, 33, "Input"],
-Cell[299780, 8794, 147040, 4357, 3533, "Output"]
+Cell[302436, 8803, 534, 18, 33, "Input"],
+Cell[302973, 8823, 154837, 4560, 3394, "Output"]
}, Open ]],
-Cell[446835, 13154, 1289, 34, 100, "Text"],
+Cell[457825, 13386, 1289, 34, 82, "Text"],
Cell[CellGroupData[{
-Cell[448149, 13192, 869, 26, 52, "Input"],
-Cell[449021, 13220, 4579, 142, 141, "Output"]
+Cell[459139, 13424, 869, 26, 52, "Input"],
+Cell[460011, 13452, 4579, 142, 141, "Output"]
}, Open ]],
-Cell[453615, 13365, 130, 3, 29, "Text"],
+Cell[464605, 13597, 130, 3, 29, "Text"],
Cell[CellGroupData[{
-Cell[453770, 13372, 1170, 34, 52, "Input"],
-Cell[454943, 13408, 189, 5, 49, "Output"]
+Cell[464760, 13604, 1227, 35, 52, "Input"],
+Cell[465990, 13641, 189, 5, 49, "Output"]
}, Open ]],
-Cell[455147, 13416, 154, 3, 29, "Text"],
+Cell[466194, 13649, 154, 3, 29, "Text"],
Cell[CellGroupData[{
-Cell[455326, 13423, 171, 6, 31, "Input"],
-Cell[455500, 13431, 4329, 133, 139, "Output"]
+Cell[466373, 13656, 171, 6, 31, "Input"],
+Cell[466547, 13664, 4286, 131, 125, "Output"]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
-Cell[459878, 13570, 44, 0, 71, "Section"],
-Cell[459925, 13572, 253, 6, 31, "Input"],
-Cell[460181, 13580, 262, 5, 47, "Text"],
+Cell[470882, 13801, 44, 0, 71, "Section"],
+Cell[470929, 13803, 253, 6, 31, "Input"],
+Cell[471185, 13811, 262, 5, 47, "Text"],
Cell[CellGroupData[{
-Cell[460468, 13589, 2077, 62, 112, "Input"],
-Cell[462548, 13653, 1997, 37, 230, "Output"],
-Cell[464548, 13692, 2364, 43, 232, "Output"]
+Cell[471472, 13820, 2135, 63, 112, "Input"],
+Cell[473610, 13885, 1997, 37, 230, "Output"],
+Cell[475610, 13924, 2364, 43, 232, "Output"]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
-Cell[466961, 13741, 49, 0, 71, "Section"],
+Cell[478023, 13973, 49, 0, 71, "Section"],
Cell[CellGroupData[{
-Cell[467035, 13745, 2204, 59, 112, "Input"],
-Cell[469242, 13806, 612, 16, 55, "Output"],
-Cell[469857, 13824, 213, 2, 145, "Output"]
+Cell[478097, 13977, 2204, 59, 112, "Input"],
+Cell[480304, 14038, 612, 16, 33, "Output"],
+Cell[480919, 14056, 213, 2, 145, "Output"]
}, Open ]],
-Cell[470085, 13829, 43, 0, 29, "Text"],
+Cell[481147, 14061, 43, 0, 29, "Text"],
Cell[CellGroupData[{
-Cell[470153, 13833, 884, 26, 52, "Input"],
-Cell[471040, 13861, 1283, 40, 50, "Output"],
-Cell[472326, 13903, 1948, 57, 50, "Output"]
+Cell[481215, 14065, 884, 26, 52, "Input"],
+Cell[482102, 14093, 1283, 40, 50, "Output"],
+Cell[483388, 14135, 1948, 57, 50, "Output"]
}, Open ]],
Cell[CellGroupData[{
-Cell[474311, 13965, 5820, 159, 225, "Input"],
-Cell[480134, 14126, 56587, 1603, 1854, "Output"],
-Cell[536724, 15731, 5304, 74, 1152, "Output"]
+Cell[485373, 14197, 5876, 159, 205, "Input"],
+Cell[491252, 14358, 58538, 1653, 1720, "Output"],
+Cell[549793, 16013, 5569, 78, 1057, "Output"]
}, Open ]],
Cell[CellGroupData[{
-Cell[542065, 15810, 215, 6, 31, "Input"],
-Cell[542283, 15818, 120, 1, 30, "Output"]
+Cell[555399, 16096, 215, 6, 31, "Input"],
+Cell[555617, 16104, 120, 1, 30, "Output"]
}, Open ]],
-Cell[542418, 15822, 61, 2, 29, "Text"],
+Cell[555752, 16108, 61, 2, 29, "Text"],
Cell[CellGroupData[{
-Cell[542504, 15828, 5367, 145, 222, "Input"],
-Cell[547874, 15975, 5875, 158, 219, "Output"],
-Cell[553752, 16135, 722, 9, 962, "Output"]
+Cell[555838, 16114, 5367, 145, 203, "Input"],
+Cell[561208, 16261, 5875, 158, 164, "Output"],
+Cell[567086, 16421, 722, 9, 962, "Output"]
}, Open ]],
Cell[CellGroupData[{
-Cell[554511, 16149, 210, 6, 31, "Input"],
-Cell[554724, 16157, 116, 1, 30, "Output"]
+Cell[567845, 16435, 210, 6, 31, "Input"],
+Cell[568058, 16443, 116, 1, 30, "Output"]
}, Open ]],
-Cell[554855, 16161, 63, 2, 29, "Text"],
-Cell[554921, 16165, 127, 3, 29, "Text"],
+Cell[568189, 16447, 63, 2, 29, "Text"],
+Cell[568255, 16451, 127, 3, 29, "Text"],
Cell[CellGroupData[{
-Cell[555073, 16172, 1162, 35, 52, "Input"],
-Cell[556238, 16209, 18801, 470, 1163, "Output"]
+Cell[568407, 16458, 1162, 35, 52, "Input"],
+Cell[569572, 16495, 18801, 470, 1163, "Output"]
}, Open ]],
-Cell[575054, 16682, 99, 2, 29, "Text"],
-Cell[575156, 16686, 709, 21, 31, "Input"],
-Cell[575868, 16709, 149, 3, 29, "Text"],
+Cell[588388, 16968, 99, 2, 29, "Text"],
+Cell[588490, 16972, 709, 21, 31, "Input"],
+Cell[589202, 16995, 149, 3, 29, "Text"],
Cell[CellGroupData[{
-Cell[576042, 16716, 2094, 53, 181, "Input"],
-Cell[578139, 16771, 5151, 151, 179, "Output"]
+Cell[589376, 17002, 2094, 53, 159, "Input"],
+Cell[591473, 17057, 5151, 151, 137, "Output"]
}, Open ]],
-Cell[583305, 16925, 86, 2, 29, "Text"],
+Cell[596639, 17211, 86, 2, 29, "Text"],
Cell[CellGroupData[{
-Cell[583416, 16931, 916, 26, 33, "Input"],
-Cell[584335, 16959, 674, 23, 30, "Output"]
+Cell[596750, 17217, 916, 26, 33, "Input"],
+Cell[597669, 17245, 674, 23, 30, "Output"]
}, Open ]],
-Cell[585024, 16985, 115, 3, 29, "Text"],
+Cell[598358, 17271, 115, 3, 29, "Text"],
Cell[CellGroupData[{
-Cell[585164, 16992, 721, 21, 52, "Input"],
-Cell[585888, 17015, 5290, 162, 160, "Output"]
+Cell[598498, 17278, 721, 21, 52, "Input"],
+Cell[599222, 17301, 5290, 162, 124, "Output"]
}, Open ]],
-Cell[591193, 17180, 64, 2, 29, "Text"],
+Cell[604527, 17466, 64, 2, 29, "Text"],
Cell[CellGroupData[{
-Cell[591282, 17186, 3489, 94, 194, "Input"],
-Cell[594774, 17282, 2396, 64, 92, "Output"],
-Cell[597173, 17348, 659, 9, 145, "Output"],
-Cell[597835, 17359, 415, 5, 107, "Output"],
-Cell[598253, 17366, 124, 1, 30, "Output"]
+Cell[604616, 17472, 3489, 94, 154, "Input"],
+Cell[608108, 17568, 2396, 64, 92, "Output"],
+Cell[610507, 17634, 659, 9, 107, "Output"],
+Cell[611169, 17645, 415, 5, 107, "Output"],
+Cell[611587, 17652, 124, 1, 30, "Output"]
}, Open ]]
}, Open ]]
}