(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 117904, 3446] NotebookOptionsPosition[ 115170, 3348] NotebookOutlinePosition[ 115527, 3364] CellTagsIndexPosition[ 115484, 3361] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["setup the system", "Section"], Cell["This loads the package.", "MathCaption", CellID->836781195], Cell[BoxData[ RowBox[{"<<", "AtomicDensityMatrix`"}]], "Input", CellID->2058623809], Cell[TextData[{ "We define an atomic system consisting of two even-parity lower states and \ two odd-parity upper states. We apply a light field with components at \ frequencies ", Cell[BoxData[ FormBox[ StyleBox[ SubscriptBox["\[Omega]", "1"], "InlineMath"], TraditionalForm]]], " (near resonant with the ", Cell[BoxData[ StyleBox[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"|", "1"}], "\[RightAngleBracket]"}], "\[Rule]", RowBox[{"|", "3"}]}], "\[RightAngleBracket]"}], "InlineMath"]]], " transition), ", Cell[BoxData[ FormBox[ StyleBox[ SubscriptBox["\[Omega]", "2"], "InlineMath"], TraditionalForm]]], " (near resonant with the ", Cell[BoxData[ StyleBox[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"|", "2"}], "\[RightAngleBracket]"}], "\[Rule]", RowBox[{"|", "3"}]}], "\[RightAngleBracket]"}], "InlineMath"]]], " transition), and ", Cell[BoxData[ FormBox[ StyleBox[ SubscriptBox["\[Omega]", "c"], "InlineMath"], TraditionalForm]]], " (near resonant with the ", Cell[BoxData[ StyleBox[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"|", "2"}], "\[RightAngleBracket]"}], "\[Rule]", RowBox[{"|", "4"}]}], "\[RightAngleBracket]"}], "InlineMath"]]], " transition)." }], "Text", CellID->525777075], Cell["Define the atomic system.", "MathCaption", CellID->429217524], Cell[BoxData[ RowBox[{ RowBox[{"system", "=", RowBox[{"Sublevels", "[", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"AtomicState", "[", RowBox[{"1", ",", RowBox[{"J", "\[Rule]", "0"}], ",", RowBox[{"L", "\[Rule]", "0"}], ",", RowBox[{"S", "\[Rule]", "0"}], ",", RowBox[{"NaturalWidth", "\[Rule]", "0"}], ",", RowBox[{"Energy", "\[Rule]", "0"}], ",", RowBox[{ RowBox[{"Polarizability", "[", "0", "]"}], "\[Rule]", "0"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"AtomicState", "[", RowBox[{"2", ",", RowBox[{"J", "\[Rule]", "1"}], ",", RowBox[{"L", "\[Rule]", "1"}], ",", RowBox[{"S", "\[Rule]", "0"}], ",", RowBox[{"NaturalWidth", "\[Rule]", "\[CapitalGamma]"}], ",", RowBox[{"Energy", "\[Rule]", "\[Omega]0"}], ",", RowBox[{ RowBox[{"Polarizability", "[", "0", "]"}], "\[Rule]", RowBox[{ RowBox[{"-", "2"}], "/", "3"}]}], ",", RowBox[{ RowBox[{"Polarizability", "[", "2", "]"}], "\[Rule]", RowBox[{"2", "/", "3"}]}], ",", RowBox[{ RowBox[{"BranchingRatio", "[", "1", "]"}], "\[Rule]", "1"}]}], "]"}]}], "\[IndentingNewLine]", "}"}], "]"}]}], ";"}]], "Input", CellID->433132487], Cell["Define the optical field ", "MathCaption", CellID->133602844], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"field", "=", RowBox[{ FractionBox[ RowBox[{"2", SqrtBox["6"]}], RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}]], RowBox[{"{", RowBox[{ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", RowBox[{"(", RowBox[{"\[Phi]x", "-", RowBox[{"t", " ", "\[Omega]"}]}], ")"}]}]], " ", "\[CapitalOmega]Rx"}], " ", ",", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", RowBox[{"(", RowBox[{"\[Phi]y", "-", RowBox[{"t", " ", "\[Omega]"}]}], ")"}]}]], " ", "\[CapitalOmega]Ry"}], " ", ",", "0"}], "}"}]}]}]], "Input", CellID->26742303], Cell[BoxData[ RowBox[{"{", RowBox[{ FractionBox[ RowBox[{"2", " ", SqrtBox["6"], " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", RowBox[{"(", RowBox[{"\[Phi]x", "-", RowBox[{"t", " ", "\[Omega]"}]}], ")"}]}]], " ", "\[CapitalOmega]Rx"}], RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}]], ",", FractionBox[ RowBox[{"2", " ", SqrtBox["6"], " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", RowBox[{"(", RowBox[{"\[Phi]y", "-", RowBox[{"t", " ", "\[Omega]"}]}], ")"}]}]], " ", "\[CapitalOmega]Ry"}], RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}]], ",", "0"}], "}"}]], "Output"] }, Open ]], Cell["The Hamiltonian for the system subject to the optical field.", \ "MathCaption", CellID->462076121], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"H", "=", RowBox[{ RowBox[{ RowBox[{"Expand", "@", RowBox[{"Hamiltonian", "[", RowBox[{"system", ",", RowBox[{"ElectricField", "\[Rule]", "field"}], ",", RowBox[{"MagneticField", "\[Rule]", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"\[CapitalOmega]rf", " ", RowBox[{"Sin", "[", RowBox[{"\[Omega]rf", " ", "t"}], "]"}]}], ",", "0", ",", "\[CapitalOmega]L"}], "}"}], "/", "BohrMagneton"}]}]}], "]"}]}], "+", RowBox[{"Hamiltonian", "[", RowBox[{"system", ",", RowBox[{"ElectricField", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0", ",", SqrtBox["\[CapitalDelta]0"]}], "}"}]}], ",", RowBox[{"Interaction", "\[Rule]", RowBox[{"{", "Polarizability", "}"}]}]}], "]"}]}], "//", "FullSimplify"}]}], "]"}]], "Input", CellID->494599775], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", RowBox[{ RowBox[{"-", "2"}], " ", RowBox[{"(", RowBox[{ RowBox[{"\[CapitalOmega]Rx", " ", RowBox[{"Cos", "[", RowBox[{"\[Phi]x", "-", RowBox[{"t", " ", "\[Omega]"}]}], "]"}]}], "+", RowBox[{"\[ImaginaryI]", " ", "\[CapitalOmega]Ry", " ", RowBox[{"Cos", "[", RowBox[{"\[Phi]y", "-", RowBox[{"t", " ", "\[Omega]"}]}], "]"}]}]}], ")"}]}], "0", RowBox[{ RowBox[{"2", " ", "\[CapitalOmega]Rx", " ", RowBox[{"Cos", "[", RowBox[{"\[Phi]x", "-", RowBox[{"t", " ", "\[Omega]"}]}], "]"}]}], "-", RowBox[{"2", " ", "\[ImaginaryI]", " ", "\[CapitalOmega]Ry", " ", RowBox[{"Cos", "[", RowBox[{"\[Phi]y", "-", RowBox[{"t", " ", "\[Omega]"}]}], "]"}]}]}]}, { RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "\[CapitalOmega]Rx", " ", RowBox[{"Cos", "[", RowBox[{"\[Phi]x", "-", RowBox[{"t", " ", "\[Omega]"}]}], "]"}]}], "+", RowBox[{"2", " ", "\[ImaginaryI]", " ", "\[CapitalOmega]Ry", " ", RowBox[{"Cos", "[", RowBox[{"\[Phi]y", "-", RowBox[{"t", " ", "\[Omega]"}]}], "]"}]}]}], RowBox[{"\[Omega]0", "+", "\[CapitalOmega]L"}], FractionBox[ RowBox[{"\[CapitalOmega]rf", " ", RowBox[{"Sin", "[", RowBox[{"t", " ", "\[Omega]rf"}], "]"}]}], SqrtBox["2"]], "0"}, {"0", FractionBox[ RowBox[{"\[CapitalOmega]rf", " ", RowBox[{"Sin", "[", RowBox[{"t", " ", "\[Omega]rf"}], "]"}]}], SqrtBox["2"]], RowBox[{"\[CapitalDelta]0", "+", "\[Omega]0"}], FractionBox[ RowBox[{"\[CapitalOmega]rf", " ", RowBox[{"Sin", "[", RowBox[{"t", " ", "\[Omega]rf"}], "]"}]}], SqrtBox["2"]]}, { RowBox[{ RowBox[{"2", " ", "\[CapitalOmega]Rx", " ", RowBox[{"Cos", "[", RowBox[{"\[Phi]x", "-", RowBox[{"t", " ", "\[Omega]"}]}], "]"}]}], "+", RowBox[{"2", " ", "\[ImaginaryI]", " ", "\[CapitalOmega]Ry", " ", RowBox[{"Cos", "[", RowBox[{"\[Phi]y", "-", RowBox[{"t", " ", "\[Omega]"}]}], "]"}]}]}], "0", FractionBox[ RowBox[{"\[CapitalOmega]rf", " ", RowBox[{"Sin", "[", RowBox[{"t", " ", "\[Omega]rf"}], "]"}]}], SqrtBox["2"]], RowBox[{"\[Omega]0", "-", "\[CapitalOmega]L"}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output"] }, Open ]], Cell["The level diagram for the system.", "MathCaption", CellID->358620443], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"LevelDiagram", "[", RowBox[{"system", ",", RowBox[{"H", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Omega]0", "\[Rule]", "1"}], ",", RowBox[{"\[CapitalOmega]rf", "\[Rule]", ".5"}], ",", RowBox[{"\[CapitalOmega]L", "\[Rule]", ".3"}], ",", RowBox[{"\[CapitalDelta]0", "\[Rule]", ".1"}]}], "}"}]}]}], "]"}]], "Input", CellID->167259034], Cell[BoxData[ GraphicsBox[{{{{}, LineBox[{{-0.4, 0}, {0.4, 0}}]}, {{}, LineBox[{{0.6, 1.3}, {1.4, 1.3}}]}, {{}, LineBox[{{-0.4, 1.1}, {0.4, 1.1}}]}, {{}, LineBox[{{-1.4, 0.7}, {-0.6, 0.7}}]}}, {{}, {}, {}}, {Arrowheads[{-0.05, 0.05}], ArrowBox[{{0.04000000000000001, 0}, {0.96, 1.3}}], ArrowBox[{{0.04000000000000001, 1.1}, {0.96, 1.3}}], ArrowBox[{{-0.04000000000000001, 0}, {-0.96, 0.7}}], ArrowBox[{{-0.96, 0.7}, {-0.04000000000000001, 1.1}}]}, {PointSize[0.0225]}}, ImagePadding->{{2, 2}, {2, 2}}, ImageSize->144.]], "Output"] }, Open ]], Cell["Apply the rotating-wave approximation to the Hamiltonian.", \ "MathCaption", CellID->577766068], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(", RowBox[{"Hrwa", "=", RowBox[{ RowBox[{"RotatingWaveApproximation", "[", RowBox[{"system", ",", "H", ",", "\[Omega]"}], "]"}], "/.", RowBox[{"\[Omega]", "\[Rule]", RowBox[{"\[Omega]0", "+", "\[CapitalDelta]"}]}]}]}], ")"}], "//", "MatrixForm"}]], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", RowBox[{ RowBox[{ RowBox[{"-", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]]}], " ", "\[CapitalOmega]Rx"}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry"}]}], "0", RowBox[{ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx"}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry"}]}]}, { RowBox[{ RowBox[{ RowBox[{"-", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]]}], " ", "\[CapitalOmega]Rx"}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry"}]}], RowBox[{ RowBox[{"-", "\[CapitalDelta]"}], "+", "\[CapitalOmega]L"}], RowBox[{ FractionBox[ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf"}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf"}], RowBox[{"2", " ", SqrtBox["2"]}]]}], "0"}, {"0", RowBox[{ FractionBox[ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf"}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf"}], RowBox[{"2", " ", SqrtBox["2"]}]]}], RowBox[{ RowBox[{"-", "\[CapitalDelta]"}], "+", "\[CapitalDelta]0"}], RowBox[{ FractionBox[ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf"}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf"}], RowBox[{"2", " ", SqrtBox["2"]}]]}]}, { RowBox[{ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx"}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry"}]}], "0", RowBox[{ FractionBox[ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf"}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf"}], RowBox[{"2", " ", SqrtBox["2"]}]]}], RowBox[{ RowBox[{"-", "\[CapitalDelta]"}], "-", "\[CapitalOmega]L"}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output"] }, Open ]], Cell[TextData[{ Cell[BoxData[ ButtonBox["IntrinsicRelaxation", BaseStyle->"Link", ButtonData->"paclet:AtomicDensityMatrix/ref/IntrinsicRelaxation"]]], " and ", Cell[BoxData[ ButtonBox["TransitRelaxation", BaseStyle->"Link", ButtonData->"paclet:AtomicDensityMatrix/ref/TransitRelaxation"]]], " supply the relaxation matrices." }], "MathCaption", CellID->610306692], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"relax", "=", RowBox[{ RowBox[{"IntrinsicRelaxation", "[", "system", "]"}], "+", RowBox[{"TransitRelaxation", "[", RowBox[{"system", ",", "\[Gamma]t"}], "]"}]}]}], "]"}]], "Input", CellID->645617687], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"\[Gamma]t", "0", "0", "0"}, {"0", RowBox[{"\[CapitalGamma]", "+", "\[Gamma]t"}], "0", "0"}, {"0", "0", RowBox[{"\[CapitalGamma]", "+", "\[Gamma]t"}], "0"}, {"0", "0", "0", RowBox[{"\[CapitalGamma]", "+", "\[Gamma]t"}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output"] }, Open ]], Cell["Remove explict time dependence from the density matrix.", "MathCaption", CellID->690131918], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"SetOptions", "[", RowBox[{"DensityMatrix", ",", RowBox[{"TimeDependence", "\[Rule]", "False"}], ",", RowBox[{"ComplexExpandVariables", "\[Rule]", "False"}]}], "]"}]], "Input", CellID->718931880], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"TimeDependence", "\[Rule]", "False"}], ",", RowBox[{"Representation", "\[Rule]", "Zeeman"}], ",", RowBox[{"DMSymbol", "\[Rule]", "\[Rho]"}], ",", RowBox[{"Label", "\[Rule]", "None"}], ",", RowBox[{"ComplexExpandVariables", "\[Rule]", "False"}], ",", RowBox[{"TimeVariable", "\[Rule]", "t"}]}], "}"}]], "Output", ImageSize->{432, 33}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}] }, Open ]], Cell[TextData[{ Cell[BoxData[ ButtonBox["OpticalRepopulation", BaseStyle->"Link", ButtonData->"paclet:AtomicDensityMatrix/ref/OpticalRepopulation"]]], " and ", Cell[BoxData[ ButtonBox["TransitRepopulation", BaseStyle->"Link", ButtonData->"paclet:AtomicDensityMatrix/ref/TransitRepopulation"]]], " supply the repopulation matrices." }], "MathCaption", CellID->854192725], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"repop", "=", RowBox[{ RowBox[{ RowBox[{"OpticalRepopulation", "[", "system", "]"}], "+", RowBox[{"TransitRepopulation", "[", RowBox[{"system", ",", "\[Gamma]t"}], "]"}]}], "/.", RowBox[{ RowBox[{"BranchingRatio", "[", RowBox[{"a_", ",", "b_"}], "]"}], "\[Rule]", SubscriptBox["R", RowBox[{"a", ",", "b"}]]}]}]}], "]"}]], "Input", CellID->465762594], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"\[Gamma]t", "+", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "+", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}]}], "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output"] }, Open ]], Cell["Here are the evolution equations.", "MathCaption", CellID->314466782], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"TableForm", "[", RowBox[{ RowBox[{"eqs", "=", RowBox[{ RowBox[{"LiouvilleEquation", "[", RowBox[{"system", ",", "Hrwa", ",", "relax", ",", "repop"}], "]"}], "//", "Expand"}]}], ",", RowBox[{"TableHeadings", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"DMVariables", "[", "system", "]"}], ",", "None"}], "}"}]}]}], "]"}]], "Input", CellID->298399236], Cell[BoxData[ TagBox[ TagBox[GridBox[{ { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{"\[Gamma]t", "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "+", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{ FractionBox["1", "2"], " ", "\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", "\[CapitalOmega]L", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]]}], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{ FractionBox["1", "2"], " ", "\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]0", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{ FractionBox["1", "2"], " ", "\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", "\[CapitalOmega]L", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{ FractionBox["1", "2"], " ", "\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", "\[CapitalOmega]L", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]0", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", "\[CapitalOmega]L", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"2", " ", "\[ImaginaryI]", " ", "\[CapitalOmega]L", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{ FractionBox["1", "2"], " ", "\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]0", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]0", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", "\[CapitalOmega]L", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]0", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", "\[CapitalOmega]L", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{ FractionBox["1", "2"], " ", "\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", "\[CapitalOmega]L", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "+", RowBox[{"2", " ", "\[ImaginaryI]", " ", "\[CapitalOmega]L", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", "\[CapitalDelta]0", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", "\[CapitalOmega]L", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]]}]}]}, { TagBox[ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]], HoldForm], RowBox[{"0", "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[CapitalOmega]Rx", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[CapitalOmega]Ry", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], "-", RowBox[{"\[CapitalGamma]", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", RowBox[{"\[Gamma]t", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Omega]rf"}]], " ", "\[CapitalOmega]rf", " ", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}]}]]}], RowBox[{"2", " ", SqrtBox["2"]}]]}]}]} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxDividers->{ "Columns" -> {False, {True}, False}, "ColumnsIndexed" -> {}, "Rows" -> {{False}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.5599999999999999]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], OutputFormsDump`HeadedColumn], Function[BoxForm`e$, TableForm[BoxForm`e$, TableHeadings -> {{ Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {1, 0}, {1, 0}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {1, 0}, {2, 1}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {1, 0}, {2, 0}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {1, 0}, {2, -1}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, 1}, {1, 0}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, 1}, {2, 1}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, 1}, {2, 0}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, 1}, {2, -1}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, 0}, {1, 0}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, 0}, {2, 1}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, 0}, {2, 0}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, 0}, {2, -1}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, -1}, {1, 0}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, -1}, {2, 1}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, -1}, {2, 0}], Subscript[AtomicDensityMatrix`DensityMatrix`\[Rho], {2, -1}, {2, -1}]}, None}]]]], "Output"] }, Open ]], Cell["Convert to c form.", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"Collect", "[", RowBox[{ RowBox[{ RowBox[{"eqs", "[", RowBox[{"[", RowBox[{"All", ",", "2"}], "]"}], "]"}], "/.", RowBox[{"{", RowBox[{ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]x"}]], " ", "\[Rule]", RowBox[{"Ex", "/", "\[CapitalOmega]Rx"}]}], ",", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]x"}]], " ", "\[Rule]", RowBox[{"Exc", "/", "\[CapitalOmega]Rx"}]}], ",", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "\[Phi]y"}]], " ", "\[Rule]", RowBox[{"Ey", "/", "\[CapitalOmega]Ry"}]}], ",", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]y"}]], " ", "\[Rule]", RowBox[{"Eyc", "/", "\[CapitalOmega]Ry"}]}], ",", " ", RowBox[{"\[Gamma]t", "\[Rule]", "gt"}], ",", RowBox[{"\[CapitalGamma]", "\[Rule]", "g0"}], ",", RowBox[{"\[CapitalOmega]L", "\[Rule]", "deltaL"}], ",", " ", RowBox[{"\[CapitalDelta]0", "\[Rule]", RowBox[{"delta0z", "[", "z", "]"}]}], ",", RowBox[{"\[CapitalDelta]", "\[Rule]", "delta"}]}], "}"}]}], ",", RowBox[{"DMElementPattern", "[", "]"}], ",", "FullSimplify"}], "]"}], "/.", " ", RowBox[{ RowBox[{"Sin", "[", RowBox[{"t", " ", "\[Omega]rf"}], "]"}], "\[Rule]", RowBox[{ SqrtBox["2"], RowBox[{"rf", "/", " ", "\[CapitalOmega]rf"}]}]}]}], "/.", RowBox[{"\[Rho]", "\[Rule]", "r"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"MapThread", "[", RowBox[{"Equal", ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"DMVariables", "[", "system", "]"}], "/.", RowBox[{"\[Rho]", "\[Rule]", "dr"}]}], ",", "%"}], "}"}]}], "]"}], "/.", RowBox[{ RowBox[{"Complex", "[", RowBox[{"0", ",", "a_"}], "]"}], "\[Rule]", RowBox[{"a", " ", "i"}]}]}], "/.", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Label", "[", RowBox[{"system", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", RowBox[{"M", "[", RowBox[{"system", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}], "\[Rule]", "i"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "system", "]"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"DeleteCases", "[", RowBox[{"%", ",", RowBox[{ RowBox[{ SubscriptBox["dr", RowBox[{"a_", ",", "b_"}]], "\[Equal]", "_"}], "/;", RowBox[{"b", "<", "a"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"StringJoin", "[", RowBox[{ RowBox[{ RowBox[{"StringReplace", "[", RowBox[{ RowBox[{ RowBox[{"ToString", "@", RowBox[{"CForm", "[", "#", "]"}]}], "<>", "\"\<;\\n\>\""}], ",", RowBox[{"{", RowBox[{ RowBox[{"\"\<==\>\"", "\[Rule]", "\"\<=\>\""}], ",", RowBox[{ RowBox[{ "\"\\"", "~~", "a_", "~~", "\"\<,\>\"", "~~", "b_", "~~", "\"\<)\>\""}], ":>", RowBox[{"\"\\"", "<>", "a", "<>", "b", "<>", "\"\<_dt\>\""}]}], ",", RowBox[{ RowBox[{ "\"\\"", "~~", "r_", "~~", "\"\<,\>\"", "~~", "a_", "~~", "\"\<,\>\"", "~~", "b_", "~~", "\"\<)\>\""}], "\[RuleDelayed]", RowBox[{"r", "<>", "a", "<>", "b"}]}]}], "}"}]}], "]"}], "&"}], "/@", "%"}], "]"}], "\[IndentingNewLine]", RowBox[{"Export", "[", RowBox[{ RowBox[{"ToFileName", "[", RowBox[{ RowBox[{"NotebookDirectory", "[", "]"}], ",", "\"\\""}], "]"}], ",", "%"}], "]"}]}], "Input", CellChangeTimes->{{3.523190032024437*^9, 3.523190039895821*^9}}], Cell[BoxData["\<\"dr11_dt = gt - gt*r11 + (-Ey - Ex*i)*r12 + (-Ey + Ex*i)*r14 \ + (-Eyc + Exc*i)*r21 + g0*r22 + g0*r33 + (-Eyc - Exc*i)*r41 + \ g0*r44;\\ndr12_dt = (Eyc - Exc*i)*r11 + (-g0/2. - gt - delta*i + i*\ \[CapitalOmega]L)*r12 + i*rf*r13 + (-Eyc + Exc*i)*r22 + (-Eyc - \ Exc*i)*r42;\\ndr13_dt = i*rf*r12 + (-g0/2. - gt - delta*i + i*delta0z(z))*r13 \ + i*rf*r14 + (-Eyc + Exc*i)*r23 + (-Eyc - Exc*i)*r43;\\ndr14_dt = (Eyc + \ Exc*i)*r11 + i*rf*r13 + (-g0/2. - gt - delta*i - i*\[CapitalOmega]L)*r14 + \ (-Eyc + Exc*i)*r24 + (-Eyc - Exc*i)*r44;\\ndr22_dt = (Ey + Ex*i)*r12 + (Eyc - \ Exc*i)*r21 + (-g0 - gt)*r22 + i*rf*r23 - i*rf*r32;\\ndr23_dt = (Ey + \ Ex*i)*r13 + i*rf*r22 + (-g0 - gt - i*\[CapitalOmega]L + i*delta0z(z))*r23 + \ i*rf*r24 - i*rf*r33;\\ndr24_dt = (Ey + Ex*i)*r14 + (Eyc + Exc*i)*r21 + \ i*rf*r23 + (-g0 - gt - 2*i*\[CapitalOmega]L)*r24 - i*rf*r34;\\ndr33_dt = \ -(i*rf*r23) + i*rf*r32 + (-g0 - gt)*r33 + i*rf*r34 - i*rf*r43;\\ndr34_dt = \ -(i*rf*r24) + (Eyc + Exc*i)*r31 + i*rf*r33 + (-g0 - gt - i*\[CapitalOmega]L - \ i*delta0z(z))*r34 - i*rf*r44;\\ndr44_dt = (Ey - Ex*i)*r14 - i*rf*r34 + (Eyc + \ Exc*i)*r41 + i*rf*r43 + (-g0 - gt)*r44;\\n\"\>"], "Output"], Cell[BoxData["\<\"C:\\\\cygwin\\\\home\\\\Simon\\\\Nresonances\\\\xmds2\\\\\ Genas_system\\\\code.txt\"\>"], "Output"] }, Open ]], Cell["Find medium polarization", "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"TensorForm", "[", RowBox[{"dcart", "=", RowBox[{"Most", "@", RowBox[{"ToCartesian", "@", RowBox[{"WignerEckart", "[", RowBox[{"system", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}]}], "]"}]}]}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", FractionBox[ RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}], SqrtBox["6"]], "0", RowBox[{"-", FractionBox[ RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}], SqrtBox["6"]]}]}, { FractionBox[ RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}], SqrtBox["6"]], "0", "0", "0"}, {"0", "0", "0", "0"}, { RowBox[{"-", FractionBox[ RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}], SqrtBox["6"]]}], "0", "0", "0"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]], ",", TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", FractionBox[ RowBox[{"\[ImaginaryI]", " ", RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}]}], SqrtBox["6"]], "0", FractionBox[ RowBox[{"\[ImaginaryI]", " ", RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}]}], SqrtBox["6"]]}, { RowBox[{"-", FractionBox[ RowBox[{"\[ImaginaryI]", " ", RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}]}], SqrtBox["6"]]}], "0", "0", "0"}, {"0", "0", "0", "0"}, { RowBox[{"-", FractionBox[ RowBox[{"\[ImaginaryI]", " ", RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}]}], SqrtBox["6"]]}], "0", "0", "0"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{ "4", "\[Pi]", " ", "\[ImaginaryI]", " ", "\[Omega]0", " ", "n", " ", FractionBox[ RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}], RowBox[{"2", SqrtBox["6"]}]], " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"DensityMatrix", "[", "system", "]"}], "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], RowBox[{"#", "[", RowBox[{"[", RowBox[{"j", ",", "i"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "system", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "i"}], "}"}]}], "]"}]}], "&"}], "/@", "dcart"}], "//", "Simplify"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"ExpandDipoleRME", "[", RowBox[{"system", ",", "%"}], " ", "]"}], "/.", RowBox[{"\[Omega]0", "\[Rule]", RowBox[{"2", RowBox[{"\[Pi]", "/", "\[Lambda]"}]}]}]}], "//", "Simplify"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"%", " ", RowBox[{"eta", "/", RowBox[{"(", FractionBox[ RowBox[{"3", " ", "n", " ", "\[CapitalGamma]", " ", SuperscriptBox["\[Lambda]", "2"], " "}], RowBox[{"16", " ", "\[Pi]"}]], ")"}]}]}], "/.", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Label", "[", RowBox[{"system", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", RowBox[{"M", "[", RowBox[{"system", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}], "\[Rule]", "i"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "system", "]"}]}], "}"}]}], "]"}]}], "/.", RowBox[{"\[Rho]", "\[Rule]", "r"}]}]}], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", FractionBox["1", "3"]}], " ", "\[ImaginaryI]", " ", "n", " ", "\[Pi]", " ", "\[Omega]0", " ", SuperscriptBox[ RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}], "2"], " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]], "-", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], ")"}]}], ",", RowBox[{ RowBox[{"-", FractionBox["1", "3"]}], " ", "n", " ", "\[Pi]", " ", "\[Omega]0", " ", SuperscriptBox[ RowBox[{"ReducedME", "[", RowBox[{"1", ",", RowBox[{"{", RowBox[{"Dipole", ",", "1"}], "}"}], ",", "2"}], "]"}], "2"], " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]], "+", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], ")"}]}]}], "}"}]], "Output"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"3", " ", "\[ImaginaryI]", " ", "n", " ", "\[CapitalGamma]", " ", SuperscriptBox["\[Lambda]", "2"], " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]], "-", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], ")"}]}], RowBox[{"16", " ", "\[Pi]"}]]}], ",", RowBox[{"-", FractionBox[ RowBox[{"3", " ", "n", " ", "\[CapitalGamma]", " ", SuperscriptBox["\[Lambda]", "2"], " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]], "+", SubscriptBox["\[Rho]", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}]]}], ")"}]}], RowBox[{"16", " ", "\[Pi]"}]]}]}], "}"}]], "Output"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "eta", " ", RowBox[{"(", RowBox[{ RowBox[{"-", SubscriptBox["r", RowBox[{"2", ",", "1"}]]}], "+", SubscriptBox["r", RowBox[{"4", ",", "1"}]]}], ")"}]}], ",", RowBox[{ RowBox[{"-", "eta"}], " ", RowBox[{"(", RowBox[{ SubscriptBox["r", RowBox[{"2", ",", "1"}]], "+", SubscriptBox["r", RowBox[{"4", ",", "1"}]]}], ")"}]}]}], "}"}]], "Output"] }, Open ]] }, Open ]] }, WindowSize->{984, 1031}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, ShowSelection->True, FrontEndVersion->"8.0 for Linux x86 (32-bit) (February 23, 2011)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[567, 22, 35, 0, 74, "Section"], Cell[605, 24, 66, 1, 44, "MathCaption", CellID->836781195], Cell[674, 27, 85, 2, 30, "Input", CellID->2058623809], Cell[762, 31, 1290, 44, 71, "Text", CellID->525777075], Cell[2055, 77, 68, 1, 44, "MathCaption", CellID->429217524], Cell[2126, 80, 1346, 33, 107, "Input", CellID->433132487], Cell[3475, 115, 68, 1, 44, "MathCaption", CellID->133602844], Cell[CellGroupData[{ Cell[3568, 120, 818, 26, 53, "Input", CellID->26742303], Cell[4389, 148, 922, 29, 53, "Output"] }, Open ]], Cell[5326, 180, 105, 2, 44, "MathCaption", CellID->462076121], Cell[CellGroupData[{ Cell[5456, 186, 987, 27, 86, "Input", CellID->494599775], Cell[6446, 215, 3061, 82, 113, "Output"] }, Open ]], Cell[9522, 300, 76, 1, 44, "MathCaption", CellID->358620443], Cell[CellGroupData[{ Cell[9623, 305, 401, 11, 30, "Input", CellID->167259034], Cell[10027, 318, 581, 12, 87, "Output"] }, Open ]], Cell[10623, 333, 102, 2, 44, "MathCaption", CellID->577766068], Cell[CellGroupData[{ Cell[10750, 339, 330, 9, 30, "Input"], Cell[11083, 350, 4689, 131, 123, "Output"] }, Open ]], Cell[15787, 484, 384, 12, 44, "MathCaption", CellID->610306692], Cell[CellGroupData[{ Cell[16196, 500, 272, 7, 30, "Input", CellID->645617687], Cell[16471, 509, 814, 22, 72, "Output"] }, Open ]], Cell[17300, 534, 98, 1, 44, "MathCaption", CellID->690131918], Cell[CellGroupData[{ Cell[17423, 539, 230, 5, 30, "Input", CellID->718931880], Cell[17656, 546, 471, 11, 50, "Output"] }, Open ]], Cell[18142, 560, 390, 12, 44, "MathCaption", CellID->854192725], Cell[CellGroupData[{ Cell[18557, 576, 459, 13, 50, "Input", CellID->465762594], Cell[19019, 591, 1463, 43, 75, "Output"] }, Open ]], Cell[20497, 637, 76, 1, 44, "MathCaption", CellID->314466782], Cell[CellGroupData[{ Cell[20598, 642, 424, 13, 50, "Input", CellID->298399236], Cell[21025, 657, 79447, 2252, 558, "Output"] }, Open ]], Cell[100487, 2912, 34, 0, 30, "Text"], Cell[CellGroupData[{ Cell[100546, 2916, 4144, 113, 237, "Input"], Cell[104693, 3031, 1184, 15, 297, "Output"], Cell[105880, 3048, 118, 1, 30, "Output"] }, Open ]], Cell[106013, 3052, 40, 0, 30, "Text"], Cell[CellGroupData[{ Cell[106078, 3056, 288, 8, 30, "Input"], Cell[106369, 3066, 3317, 96, 229, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[109723, 3167, 1976, 61, 153, "Input"], Cell[111702, 3230, 1530, 49, 83, "Output"], Cell[113235, 3281, 1375, 42, 52, "Output"], Cell[114613, 3325, 529, 19, 30, "Output"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)