summaryrefslogtreecommitdiff
path: root/mathemathica_fwm/RbXMDSSetup.nb
diff options
context:
space:
mode:
Diffstat (limited to 'mathemathica_fwm/RbXMDSSetup.nb')
-rwxr-xr-xmathemathica_fwm/RbXMDSSetup.nb492
1 files changed, 224 insertions, 268 deletions
diff --git a/mathemathica_fwm/RbXMDSSetup.nb b/mathemathica_fwm/RbXMDSSetup.nb
index 21fe46a..801b4c4 100755
--- a/mathemathica_fwm/RbXMDSSetup.nb
+++ b/mathemathica_fwm/RbXMDSSetup.nb
@@ -1389,6 +1389,71 @@ Cell[BoxData[{
Cell[CellGroupData[{
+Cell["\<\
+Eliminate conjugate DM variables\
+\>", "Section"],
+
+Cell[BoxData[{
+ RowBox[{"delvars1", "=",
+ RowBox[{"Intersection", "[",
+ RowBox[{"vars", ",",
+ RowBox[{"Cases", "[",
+ RowBox[{"vars", ",",
+ RowBox[{
+ RowBox[{
+ RowBox[{"DMElementPattern", "[", "]"}], "/;",
+ RowBox[{"!",
+ RowBox[{"OrderedQ", "[",
+ RowBox[{"{",
+ RowBox[{"State2", ",", "State1"}], "}"}], "]"}]}]}], "\[Rule]",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{"State2", ",", "State1"}]], "[", "t", "]"}]}]}], "]"}]}],
+ "]"}]}], "\[IndentingNewLine]",
+ RowBox[{"delpos1", "=",
+ RowBox[{"Flatten", "[",
+ RowBox[{
+ RowBox[{
+ RowBox[{
+ RowBox[{"Position", "[",
+ RowBox[{"vars", ",", "#"}], "]"}], "&"}], "/@", "%"}], ",", "1"}],
+ "]"}]}], "\[IndentingNewLine]",
+ RowBox[{"vars1", "=",
+ RowBox[{"Delete", "[",
+ RowBox[{"vars", ",", "delpos1"}], "]"}]}]}], "Input"],
+
+Cell[BoxData[
+ RowBox[{"delreps1", "=",
+ RowBox[{"delvars1", "/.",
+ RowBox[{
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{"s1_", ",", "s2_"}]], "[", "t", "]"}], "\[Rule]",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{"s1", ",", "s2"}]], "[", "t", "]"}], "\[Rule]",
+ RowBox[{"Conjugate", "[",
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{"s2", ",", "s1"}]], "[", "t", "]"}], "]"}]}],
+ ")"}]}]}]}]], "Input"],
+
+Cell[BoxData[
+ RowBox[{"eqs1", "=",
+ RowBox[{
+ RowBox[{"Delete", "[",
+ RowBox[{"eqs", ",", "delpos1"}], "]"}], "/.", "delreps1"}]}]], "Input"],
+
+Cell[BoxData[
+ RowBox[{"inits1", "=",
+ RowBox[{"Delete", "[",
+ RowBox[{"inits", ",", "delpos1"}], "]"}]}]], "Input"]
+}, Open ]],
+
+Cell[CellGroupData[{
+
Cell["Convert equations to XMDS form", "Section"],
Cell[BoxData[{
@@ -1516,94 +1581,86 @@ Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{
- RowBox[{"DeleteCases", "[",
+ RowBox[{
RowBox[{
RowBox[{
- RowBox[{
- RowBox[{"Join", "[",
- RowBox[{
- RowBox[{"{",
- RowBox[{
- RowBox[{"rt6", "==",
- RowBox[{"N", "[",
- SqrtBox["6"], "]"}]}], ",",
- RowBox[{"rt3", "==",
- RowBox[{"N", "[",
- SqrtBox["3"], "]"}]}], ",",
- RowBox[{"rt2", "==",
- RowBox[{"N", "[",
- SqrtBox["2"], "]"}]}]}], "}"}], ",", "eqs"}], "]"}], "/.",
- RowBox[{"{",
- RowBox[{
+ RowBox[{"Join", "[",
+ RowBox[{
+ RowBox[{"{",
RowBox[{
- SubscriptBox["\[Delta]", "s_"], "\[Rule]",
- RowBox[{
- SubscriptBox["\[Delta]", "s"], "+",
- RowBox[{
- SubscriptBox["k", "s"], "v"}]}]}], ",",
- RowBox[{"\[Delta]\[Omega]", "\[Rule]",
- RowBox[{"\[Delta]\[Omega]", "+",
- RowBox[{"\[Delta]k", " ", "v"}]}]}]}], "}"}]}], "/.",
+ RowBox[{"rt6", "==",
+ RowBox[{"N", "[",
+ SqrtBox["6"], "]"}]}], ",",
+ RowBox[{"rt3", "==",
+ RowBox[{"N", "[",
+ SqrtBox["3"], "]"}]}], ",",
+ RowBox[{"rt2", "==",
+ RowBox[{"N", "[",
+ SqrtBox["2"], "]"}]}]}], "}"}], ",", "eqs1"}], "]"}], "/.",
RowBox[{"{",
RowBox[{
- RowBox[{"NaturalWidth", "\[Rule]", "g"}], ",",
- RowBox[{"HyperfineA", "\[Rule]", "ha"}], ",",
- RowBox[{"HyperfineB", "\[Rule]", "hb"}], ",",
RowBox[{
+ SubscriptBox["\[Delta]", "s_"], "\[Rule]",
RowBox[{
- SuperscriptBox[
- SubscriptBox["\[Rho]",
- RowBox[{"s1_", ",", "s2_"}]], "\[Prime]",
- MultilineFunction->None], "[", "t", "]"}], "\[RuleDelayed]",
- RowBox[{"dr", "[",
- RowBox[{"s1", ",", "s2"}], "]"}]}], ",",
+ SubscriptBox["\[Delta]", "s"], "+",
+ RowBox[{
+ SubscriptBox["k", "s"], "v"}]}]}], ",",
+ RowBox[{"\[Delta]\[Omega]", "\[Rule]",
+ RowBox[{"\[Delta]\[Omega]", "+",
+ RowBox[{"\[Delta]k", " ", "v"}]}]}]}], "}"}]}], "/.",
+ RowBox[{"{",
+ RowBox[{
+ RowBox[{"NaturalWidth", "\[Rule]", "g"}], ",",
+ RowBox[{"HyperfineA", "\[Rule]", "ha"}], ",",
+ RowBox[{"HyperfineB", "\[Rule]", "hb"}], ",",
+ RowBox[{
RowBox[{
- RowBox[{
+ SuperscriptBox[
SubscriptBox["\[Rho]",
- RowBox[{"s1_", ",", "s2_"}]], "[", "t", "]"}], "\[RuleDelayed]",
- RowBox[{"r", "[",
- RowBox[{"s1", ",", "s2"}], "]"}]}], ",",
- RowBox[{"\[Gamma]t", "\[Rule]", "gt"}], ",",
- RowBox[{"\[CapitalOmega]L", "\[Rule]", "WL"}], ",",
- RowBox[{
- SubscriptBox["\[Delta]", "i_"], "\[Rule]",
- RowBox[{"d", "[", "i", "]"}]}], ",",
+ RowBox[{"s1_", ",", "s2_"}]], "\[Prime]",
+ MultilineFunction->None], "[", "t", "]"}], "\[RuleDelayed]",
+ RowBox[{"dr", "[",
+ RowBox[{"s1", ",", "s2"}], "]"}]}], ",",
+ RowBox[{
RowBox[{
- SubscriptBox["k", "i_"], "\[Rule]",
- RowBox[{"k", "[", "i", "]"}]}], ",",
+ SubscriptBox["\[Rho]",
+ RowBox[{"s1_", ",", "s2_"}]], "[", "t", "]"}], "\[RuleDelayed]",
+ RowBox[{"r", "[",
+ RowBox[{"s1", ",", "s2"}], "]"}]}], ",",
+ RowBox[{"\[Gamma]t", "\[Rule]", "gt"}], ",",
+ RowBox[{"\[CapitalOmega]L", "\[Rule]", "WL"}], ",",
+ RowBox[{
+ SubscriptBox["\[Delta]", "i_"], "\[Rule]",
+ RowBox[{"d", "[", "i", "]"}]}], ",",
+ RowBox[{
+ SubscriptBox["k", "i_"], "\[Rule]",
+ RowBox[{"k", "[", "i", "]"}]}], ",",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{"a_.", "+",
+ RowBox[{"\[ImaginaryI]", " ",
+ SubscriptBox["\[Phi]", "i_"]}]}]], "\[Rule]",
RowBox[{
- SuperscriptBox["\[ExponentialE]",
- RowBox[{"a_.", "+",
- RowBox[{"\[ImaginaryI]", " ",
- SubscriptBox["\[Phi]", "i_"]}]}]], "\[Rule]",
+ SuperscriptBox["\[ExponentialE]", "a"],
RowBox[{
- SuperscriptBox["\[ExponentialE]", "a"],
- RowBox[{
- RowBox[{"W", "[", "i", "]"}], "/",
- SubscriptBox["\[CapitalOmega]", "i"]}]}]}], ",",
+ RowBox[{"W", "[", "i", "]"}], "/",
+ SubscriptBox["\[CapitalOmega]", "i"]}]}]}], ",",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{"a_.", "-",
+ RowBox[{"\[ImaginaryI]", " ",
+ SubscriptBox["\[Phi]", "i_"]}]}]], "\[Rule]",
RowBox[{
- SuperscriptBox["\[ExponentialE]",
- RowBox[{"a_.", "-",
- RowBox[{"\[ImaginaryI]", " ",
- SubscriptBox["\[Phi]", "i_"]}]}]], "\[Rule]",
+ SuperscriptBox["\[ExponentialE]", "a"],
RowBox[{
- SuperscriptBox["\[ExponentialE]", "a"],
- RowBox[{
- RowBox[{"Wc", "[", "i", "]"}], "/",
- SubscriptBox["\[CapitalOmega]", "i"]}]}]}]}], "}"}]}], ",",
- RowBox[{
- RowBox[{
- RowBox[{"dr", "[",
- RowBox[{"a_", ",", "b_"}], "]"}], "\[Equal]", "_"}], "/;",
- RowBox[{"!",
- RowBox[{"OrderedQ", "[",
- RowBox[{"{",
- RowBox[{"a", ",", "b"}], "}"}], "]"}]}]}]}], "]"}], "/.",
- "labelrep"}], "/.",
- RowBox[{
- RowBox[{"Power", "[",
- RowBox[{"E", ",", "b_"}], "]"}], "\[Rule]",
- RowBox[{"exp", "[", "b", "]"}]}]}], "/.",
+ RowBox[{"Wc", "[", "i", "]"}], "/",
+ SubscriptBox["\[CapitalOmega]", "i"]}]}]}]}], "}"}]}], "/.",
+ "labelrep"}], "/.",
+ RowBox[{
+ RowBox[{"Power", "[",
+ RowBox[{"E", ",", "b_"}], "]"}], "\[Rule]",
+ RowBox[{"exp", "[", "b", "]"}]}]}], "/.",
+ RowBox[{"Conjugate", "\[Rule]", "conj"}]}], "/.",
RowBox[{
RowBox[{"Complex", "[",
RowBox[{"0", ",", "a_"}], "]"}], "\[Rule]",
@@ -1670,11 +1727,7 @@ Cell[BoxData[{
RowBox[{"\"\<delta\>\"", "<>", "a"}]}], ",",
RowBox[{
RowBox[{"\"\<k(\>\"", "~~", "a_", "~~", "\"\<)\>\""}], "\[RuleDelayed]",
- RowBox[{"\"\<Kvec\>\"", "<>", "a"}]}]}], "}"}]}], "]"}]}], "Input",
- CellChangeTimes->{{3.55456914848632*^9, 3.554569164767181*^9},
- 3.554569287340058*^9, {3.554569321671209*^9, 3.554569324047068*^9}, {
- 3.554569429250296*^9, 3.554569499438798*^9}, {3.554569642022385*^9,
- 3.554569642695839*^9}}],
+ RowBox[{"\"\<Kvec\>\"", "<>", "a"}]}]}], "}"}]}], "]"}]}], "Input"],
Cell[BoxData[
RowBox[{"Export", "[",
@@ -1691,75 +1744,65 @@ Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{
- RowBox[{"DeleteCases", "[",
- RowBox[{
- RowBox[{
- RowBox[{"inits", "/.",
- RowBox[{"{",
- RowBox[{
- RowBox[{
- SubscriptBox["\[Delta]", "s_"], "\[Rule]",
- RowBox[{
- SubscriptBox["\[Delta]", "s"], "+",
- RowBox[{
- SubscriptBox["k", "s"], "v"}]}]}], ",",
- RowBox[{"\[Delta]\[Omega]", "\[Rule]",
- RowBox[{"\[Delta]\[Omega]", "+",
- RowBox[{"\[Delta]k", " ", "v"}]}]}]}], "}"}]}], "/.",
- RowBox[{"{",
+ RowBox[{
+ RowBox[{"inits1", "/.",
+ RowBox[{"{",
+ RowBox[{
RowBox[{
- RowBox[{"NaturalWidth", "\[Rule]", "g"}], ",",
- RowBox[{"HyperfineA", "\[Rule]", "ha"}], ",",
- RowBox[{"HyperfineB", "\[Rule]", "hb"}], ",",
+ SubscriptBox["\[Delta]", "s_"], "\[Rule]",
RowBox[{
+ SubscriptBox["\[Delta]", "s"], "+",
RowBox[{
- SuperscriptBox[
- SubscriptBox["\[Rho]",
- RowBox[{"s1_", ",", "s2_"}]], "\[Prime]",
- MultilineFunction->None], "[", "t_", "]"}], "\[RuleDelayed]",
- RowBox[{"dr", "[",
- RowBox[{"s1", ",", "s2"}], "]"}]}], ",",
- RowBox[{
- RowBox[{
- SubscriptBox["\[Rho]",
- RowBox[{"s1_", ",", "s2_"}]], "[", "t_", "]"}], "\[RuleDelayed]",
- RowBox[{"r", "[",
- RowBox[{"s1", ",", "s2"}], "]"}]}], ",",
- RowBox[{"\[Gamma]t", "\[Rule]", "gt"}], ",",
- RowBox[{
- SubscriptBox["\[Delta]", "i_"], "\[Rule]",
- RowBox[{"d", "[", "i", "]"}]}], ",",
- RowBox[{
- SubscriptBox["k", "i_"], "\[Rule]",
- RowBox[{"k", "[", "i", "]"}]}], ",",
- RowBox[{
- SuperscriptBox["\[ExponentialE]",
- RowBox[{"a_.", "+",
- RowBox[{"\[ImaginaryI]", " ",
- SubscriptBox["\[Phi]", "i_"]}]}]], "\[Rule]",
- RowBox[{
- SuperscriptBox["\[ExponentialE]", "a"],
- RowBox[{
- RowBox[{"W", "[", "i", "]"}], "/",
- SubscriptBox["\[CapitalOmega]", "i"]}]}]}], ",",
- RowBox[{
- SuperscriptBox["\[ExponentialE]",
- RowBox[{"a_.", "-",
- RowBox[{"\[ImaginaryI]", " ",
- SubscriptBox["\[Phi]", "i_"]}]}]], "\[Rule]",
- RowBox[{
- SuperscriptBox["\[ExponentialE]", "a"],
- RowBox[{
- RowBox[{"Wc", "[", "i", "]"}], "/",
- SubscriptBox["\[CapitalOmega]", "i"]}]}]}]}], "}"}]}], ",",
+ SubscriptBox["k", "s"], "v"}]}]}], ",",
+ RowBox[{"\[Delta]\[Omega]", "\[Rule]",
+ RowBox[{"\[Delta]\[Omega]", "+",
+ RowBox[{"\[Delta]k", " ", "v"}]}]}]}], "}"}]}], "/.",
+ RowBox[{"{",
RowBox[{
+ RowBox[{"NaturalWidth", "\[Rule]", "g"}], ",",
+ RowBox[{"HyperfineA", "\[Rule]", "ha"}], ",",
+ RowBox[{"HyperfineB", "\[Rule]", "hb"}], ",",
+ RowBox[{
+ RowBox[{
+ SuperscriptBox[
+ SubscriptBox["\[Rho]",
+ RowBox[{"s1_", ",", "s2_"}]], "\[Prime]",
+ MultilineFunction->None], "[", "t_", "]"}], "\[RuleDelayed]",
+ RowBox[{"dr", "[",
+ RowBox[{"s1", ",", "s2"}], "]"}]}], ",",
RowBox[{
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{"s1_", ",", "s2_"}]], "[", "t_", "]"}], "\[RuleDelayed]",
RowBox[{"r", "[",
- RowBox[{"a_", ",", "b_"}], "]"}], "\[Equal]", "_"}], "/;",
- RowBox[{"!",
- RowBox[{"OrderedQ", "[",
- RowBox[{"{",
- RowBox[{"a", ",", "b"}], "}"}], "]"}]}]}]}], "]"}], "/.",
+ RowBox[{"s1", ",", "s2"}], "]"}]}], ",",
+ RowBox[{"\[Gamma]t", "\[Rule]", "gt"}], ",",
+ RowBox[{
+ SubscriptBox["\[Delta]", "i_"], "\[Rule]",
+ RowBox[{"d", "[", "i", "]"}]}], ",",
+ RowBox[{
+ SubscriptBox["k", "i_"], "\[Rule]",
+ RowBox[{"k", "[", "i", "]"}]}], ",",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{"a_.", "+",
+ RowBox[{"\[ImaginaryI]", " ",
+ SubscriptBox["\[Phi]", "i_"]}]}]], "\[Rule]",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]", "a"],
+ RowBox[{
+ RowBox[{"W", "[", "i", "]"}], "/",
+ SubscriptBox["\[CapitalOmega]", "i"]}]}]}], ",",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{"a_.", "-",
+ RowBox[{"\[ImaginaryI]", " ",
+ SubscriptBox["\[Phi]", "i_"]}]}]], "\[Rule]",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]", "a"],
+ RowBox[{
+ RowBox[{"Wc", "[", "i", "]"}], "/",
+ SubscriptBox["\[CapitalOmega]", "i"]}]}]}]}], "}"}]}], "/.",
"labelrep"}], "/.",
RowBox[{
RowBox[{"Power", "[",
@@ -1839,94 +1882,7 @@ Cell[BoxData[
RowBox[{"ToFileName", "[",
RowBox[{
RowBox[{"NotebookDirectory", "[", "]"}], ",", "\"\<RbInits.txt\>\""}],
- "]"}], ",", "%"}], "]"}]], "Input"],
-
-Cell[BoxData[{
- RowBox[{
- RowBox[{"Cases", "[",
- RowBox[{
- RowBox[{
- RowBox[{"inits", "[",
- RowBox[{"[",
- RowBox[{"All", ",", "1"}], "]"}], "]"}], "/.",
- RowBox[{
- RowBox[{
- SubscriptBox["\[Rho]",
- RowBox[{"s1_", ",", "s2_"}]], "[", "t_", "]"}], "\[RuleDelayed]",
- RowBox[{"r", "[",
- RowBox[{"s1", ",", "s2"}], "]"}]}]}], ",",
- RowBox[{
- RowBox[{"r", "[",
- RowBox[{"a_", ",", "b_"}], "]"}], "/;",
- RowBox[{"!",
- RowBox[{"OrderedQ", "[",
- RowBox[{"{",
- RowBox[{"a", ",", "b"}], "}"}], "]"}]}]}]}], "]"}], "/.",
- "labelrep"}], "\[IndentingNewLine]",
- RowBox[{"\"\<complex \>\"", "<>",
- RowBox[{"StringReplace", "[",
- RowBox[{
- RowBox[{"StringJoin", "[",
- RowBox[{"Riffle", "[",
- RowBox[{
- RowBox[{
- RowBox[{
- RowBox[{"ToString", "@",
- RowBox[{"CForm", "[", "#", "]"}]}], "&"}], "/@", "%"}], ",",
- "\"\<, \>\""}], "]"}], "]"}], ",",
- RowBox[{"{",
- RowBox[{
- RowBox[{"\"\<==\>\"", "\[Rule]", "\"\<=\>\""}], ",",
- RowBox[{
- RowBox[{"Shortest", "[",
- RowBox[{
- "\"\<dr(\\\"\>\"", "~~", "a__", "~~", "\"\<\\\",\\\"\>\"", "~~",
- "b__", "~~", "\"\<\\\")\>\""}], "]"}], ":>",
- RowBox[{"\"\<dr\>\"", "<>", "a", "<>", "b", "<>", "\"\<_dt\>\""}]}],
- ",",
- RowBox[{
- RowBox[{"Shortest", "[",
- RowBox[{
- "\"\<r(\\\"\>\"", "~~", "a__", "~~", "\"\<\\\",\\\"\>\"", "~~", "b__",
- "~~", "\"\<\\\")\>\""}], "]"}], ":>",
- RowBox[{"\"\<r\>\"", "<>", "a", "<>", "b"}]}], ",",
- RowBox[{
- RowBox[{"\"\<ha(\>\"", "~~", "a_", "~~", "\"\<)\>\""}],
- "\[RuleDelayed]",
- RowBox[{"\"\<ha\>\"", "<>", "a"}]}], ",",
- RowBox[{
- RowBox[{"\"\<hb(\>\"", "~~", "a_", "~~", "\"\<)\>\""}],
- "\[RuleDelayed]",
- RowBox[{"\"\<hb\>\"", "<>", "a"}]}], ",",
- RowBox[{
- RowBox[{"\"\<g(\>\"", "~~", "a_", "~~", "\"\<)\>\""}],
- "\[RuleDelayed]",
- RowBox[{"\"\<g\>\"", "<>", "a"}]}], ",",
- RowBox[{
- RowBox[{"\"\<W(\>\"", "~~", "a_", "~~", "\"\<)\>\""}],
- "\[RuleDelayed]",
- RowBox[{"\"\<E\>\"", "<>", "a"}]}], ",",
- RowBox[{
- RowBox[{"\"\<Wc(\>\"", "~~", "a_", "~~", "\"\<)\>\""}],
- "\[RuleDelayed]",
- RowBox[{"\"\<Ec\>\"", "<>", "a"}]}], ",",
- RowBox[{
- RowBox[{"\"\<d(\>\"", "~~", "a_", "~~", "\"\<)\>\""}],
- "\[RuleDelayed]",
- RowBox[{"\"\<d\>\"", "<>", "a"}]}], ",",
- RowBox[{
- RowBox[{"\"\<k(\>\"", "~~", "a_", "~~", "\"\<)\>\""}],
- "\[RuleDelayed]",
- RowBox[{"\"\<k\>\"", "<>", "a"}]}]}], "}"}]}], "]"}], "<>",
- "\"\<;\>\""}]}], "Input"],
-
-Cell[BoxData[
- RowBox[{"Export", "[",
- RowBox[{
- RowBox[{"ToFileName", "[",
- RowBox[{
- RowBox[{"NotebookDirectory", "[", "]"}], ",",
- "\"\<VariableDeclaration.txt\>\""}], "]"}], ",", "%"}], "]"}]], "Input"]
+ "]"}], ",", "%"}], "]"}]], "Input"]
}, Open ]],
Cell[CellGroupData[{
@@ -2050,8 +2006,7 @@ density matrix by\n\t",
StyleBox["d",
FontWeight->"Bold",
FontSlant->"Plain"],
- RowBox[{"m", ",", "n"}]], " "}]}]}]}], TraditionalForm]],
- FormatType->"TraditionalForm"]
+ RowBox[{"m", ",", "n"}]], " "}]}]}]}], TraditionalForm]]]
}], "Text"],
Cell[TextData[{
@@ -2069,8 +2024,7 @@ components of ",
StyleBox["d",
FontWeight->"Bold",
FontSlant->"Plain"],
- RowBox[{"m", ",", "n"}]]}]}], TraditionalForm]],
- FormatType->"TraditionalForm"],
+ RowBox[{"m", ",", "n"}]]}]}], TraditionalForm]]],
" for each of the four transitions, in terms of the reduced dipole matrix \
element ||d||."
}], "Text"],
@@ -2126,8 +2080,7 @@ Cell[TextData[{
Cell[BoxData[
FormBox[
RowBox[{"||", "d",
- SuperscriptBox["||", "2"]}], TraditionalForm]],
- FormatType->"TraditionalForm"],
+ SuperscriptBox["||", "2"]}], TraditionalForm]]],
" in terms of the natural width \[CapitalGamma] of the transition."
}], "Text"],
@@ -2141,17 +2094,19 @@ Cell[BoxData[
RowBox[{
"subsys", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}], ",",
RowBox[{
- RowBox[{"4", "\[Pi]", " ", "\[ImaginaryI]", " ",
- RowBox[{"Energy", "[",
- RowBox[{"Last", "[",
- RowBox[{
- "rmes", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}],
- "]"}], "]"}], " ", "n0", " ",
- RowBox[{
- "rmes", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}], " ",
- RowBox[{
- "polarizationcomponents1", "\[LeftDoubleBracket]", "k",
- "\[RightDoubleBracket]"}]}], "/.", "delreps"}]}], "]"}], "]"}], ",",
+ RowBox[{
+ RowBox[{"4", "\[Pi]", " ", "\[ImaginaryI]", " ",
+ RowBox[{"Energy", "[",
+ RowBox[{"Last", "[",
+ RowBox[{
+ "rmes", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}],
+ "]"}], "]"}], " ", "n0", " ",
+ RowBox[{
+ "rmes", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}], " ",
+ RowBox[{
+ "polarizationcomponents1", "\[LeftDoubleBracket]", "k",
+ "\[RightDoubleBracket]"}]}], "/.", "delreps"}], "/.",
+ "delreps1"}]}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"k", ",",
RowBox[{"Length", "[", "subsys", "]"}]}], "}"}]}], "]"}]}]], "Input"],
@@ -2166,8 +2121,7 @@ Cell[TextData[{
RowBox[{"3",
SubsuperscriptBox["\[Lambda]", "s", "2"],
SubscriptBox["\[CapitalGamma]", "s"], "n"}],
- RowBox[{"4", "\[Pi]"}]]}], TraditionalForm]],
- FormatType->"TraditionalForm"],
+ RowBox[{"4", "\[Pi]"}]]}], TraditionalForm]]],
"."
}], "Text"],
@@ -2264,27 +2218,29 @@ Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{
- RowBox[{"propeqs", "/.",
- RowBox[{"{",
- RowBox[{
- RowBox[{"\[Eta]", "\[Rule]", "eta"}], ",",
- RowBox[{"NaturalWidth", "\[Rule]", "g"}], ",",
+ RowBox[{
+ RowBox[{"propeqs", "/.",
+ RowBox[{"{",
RowBox[{
+ RowBox[{"\[Eta]", "\[Rule]", "eta"}], ",",
+ RowBox[{"NaturalWidth", "\[Rule]", "g"}], ",",
RowBox[{
- SubscriptBox["\[Rho]",
- RowBox[{"s1_", ",", "s2_"}]], "[", "t", "]"}], "\[RuleDelayed]",
- RowBox[{"r", "[",
- RowBox[{"s1", ",", "s2"}], "]"}]}], ",",
- RowBox[{
+ RowBox[{
+ SubscriptBox["\[Rho]",
+ RowBox[{"s1_", ",", "s2_"}]], "[", "t", "]"}], "\[RuleDelayed]",
+ RowBox[{"r", "[",
+ RowBox[{"s1", ",", "s2"}], "]"}]}], ",",
RowBox[{
- SuperscriptBox[
- SubscriptBox["\[CapitalOmega]", "i_"], "\[Prime]",
- MultilineFunction->None], "[", "z", "]"}], "\[Rule]",
- RowBox[{"dW", "[", "i", "]"}]}]}], "}"}]}], "/.", "labelrep"}], "/.",
- RowBox[{
- RowBox[{"Complex", "[",
- RowBox[{"0", ",", "a_"}], "]"}], "\[Rule]",
- RowBox[{"a", " ", "i"}]}]}], "/.",
+ RowBox[{
+ SuperscriptBox[
+ SubscriptBox["\[CapitalOmega]", "i_"], "\[Prime]",
+ MultilineFunction->None], "[", "z", "]"}], "\[Rule]",
+ RowBox[{"dW", "[", "i", "]"}]}]}], "}"}]}], "/.", "labelrep"}], "/.",
+ RowBox[{
+ RowBox[{"Complex", "[",
+ RowBox[{"0", ",", "a_"}], "]"}], "\[Rule]",
+ RowBox[{"a", " ", "i"}]}]}], "/.",
+ RowBox[{"Conjugate", "\[Rule]", "conj"}]}], "/.",
RowBox[{
RowBox[{"a", ":",
RowBox[{"Power", "[",
@@ -2356,7 +2312,7 @@ Cell[BoxData[{
"\"\<RbPropEquations.txt\>\""}], "]"}], ",", "%"}], "]"}]}], "Input"]
}, Open ]]
},
-WindowSize->{974, 898},
+WindowSize->{1331, 890},
WindowMargins->{{0, Automatic}, {Automatic, 3}},
PrivateNotebookOptions->{"FileOutlineCache"->False},
ShowSelection->True,
generated by cgit v1.2.3 (git 2.39.1) at 2024-11-16 09:33:01 +0000