summaryrefslogtreecommitdiff
path: root/xmds2
diff options
context:
space:
mode:
Diffstat (limited to 'xmds2')
-rw-r--r--xmds2/Shahriar_system/GenerateShahriarSystem.nb502
1 files changed, 278 insertions, 224 deletions
diff --git a/xmds2/Shahriar_system/GenerateShahriarSystem.nb b/xmds2/Shahriar_system/GenerateShahriarSystem.nb
index f5dcb38..9968348 100644
--- a/xmds2/Shahriar_system/GenerateShahriarSystem.nb
+++ b/xmds2/Shahriar_system/GenerateShahriarSystem.nb
@@ -10,10 +10,10 @@
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 145, 7]
-NotebookDataLength[ 63326, 1879]
-NotebookOptionsPosition[ 61101, 1799]
-NotebookOutlinePosition[ 61466, 1815]
-CellTagsIndexPosition[ 61423, 1812]
+NotebookDataLength[ 65931, 1933]
+NotebookOptionsPosition[ 63713, 1853]
+NotebookOutlinePosition[ 64070, 1869]
+CellTagsIndexPosition[ 64027, 1866]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
@@ -50,29 +50,22 @@ frequencies ",
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 ",
+ " and ",
Cell[BoxData[
FormBox[
StyleBox[
- SubscriptBox["\[Omega]", "c"], "InlineMath"], TraditionalForm]]],
- " (near resonant with the ",
+ SubscriptBox["\[Omega]", "3"], "InlineMath"], TraditionalForm]]],
+ "(near resonant with the ",
Cell[BoxData[
StyleBox[
RowBox[{
RowBox[{
RowBox[{
RowBox[{"|", "2"}], "\[RightAngleBracket]"}], "\[Rule]",
- RowBox[{"|", "4"}]}], "\[RightAngleBracket]"}], "InlineMath"]]],
- " transition)."
+ RowBox[{"|", "3"}]}], "\[RightAngleBracket]"}], "InlineMath"]]],
+ " transition)\[AliasDelimiter]."
}], "Text",
+ CellChangeTimes->{{3.524013673755671*^9, 3.524013697746722*^9}},
CellID->525777075],
Cell["\<\
@@ -190,7 +183,8 @@ Cell[BoxData[
RowBox[{"2", ",",
RowBox[{"{",
RowBox[{"Dipole", ",", "1"}], "}"}], ",", "3"}], "]"}]]}], ",", "0",
- ",", "0"}], "}"}]], "Output"]
+ ",", "0"}], "}"}]], "Output",
+ CellChangeTimes->{3.524014012716285*^9}]
}, Open ]],
Cell["\<\
@@ -268,7 +262,8 @@ Cell[BoxData[
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
Function[BoxForm`e$,
- MatrixForm[BoxForm`e$]]]], "Output"]
+ MatrixForm[BoxForm`e$]]]], "Output",
+ CellChangeTimes->{3.524014012961052*^9}]
}, Open ]],
Cell["The level diagram for the system.", "MathCaption",
@@ -289,21 +284,23 @@ Cell[BoxData[
RowBox[{"Energy", "[", "2", "]"}], "\[Rule]",
RowBox[{"-", "1"}]}], ",",
RowBox[{
- RowBox[{"Energy", "[", "4", "]"}], "\[Rule]", ".5"}]}], "}"}]}]}],
+ RowBox[{"Energy", "[", "3", "]"}], "\[Rule]", ".5"}]}], "}"}]}]}],
"]"}]], "Input",
+ CellChangeTimes->{{3.524013707327201*^9, 3.524013712245067*^9}},
CellID->167259034],
Cell[BoxData[
GraphicsBox[{{{{},
LineBox[{{-0.9, -1.5}, {-0.09999999999999998, -1.5}}]}, {{},
- LineBox[{{-0.9, -1}, {-0.09999999999999998, -1}}]}, {{},
- LineBox[{{0.09999999999999998, 0}, {0.9, 0}}]}}, {{}, {}, {}},
+ LineBox[{{-0.9, -1.}, {-0.09999999999999998, -1.}}]}, {{},
+ LineBox[{{0.09999999999999998, 0.}, {0.9, 0.}}]}}, {{}, {}, {}},
{Arrowheads[{-0.07659574468085106, 0.07659574468085106}],
ArrowBox[{{-0.45999999999999996`, -1.5}, {0.45999999999999996`, 0.}}],
- ArrowBox[{{-0.45999999999999996`, -1}, {0.45999999999999996`, 0.}}]},
+ ArrowBox[{{-0.45999999999999996`, -1.}, {0.45999999999999996`, 0.}}]},
{PointSize[0.0225]}},
- ImagePadding->{{2, 2}, {2, 2}},
- ImageSize->94.]], "Output"]
+ ImagePadding->{{2., 2}, {2., 2.}},
+ ImageSize->94.]], "Output",
+ CellChangeTimes->{3.524014013119625*^9}]
}, Open ]],
Cell["\<\
@@ -424,7 +421,8 @@ Cell[BoxData[
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
Function[BoxForm`e$,
- MatrixForm[BoxForm`e$]]]], "Output"]
+ MatrixForm[BoxForm`e$]]]], "Output",
+ CellChangeTimes->{3.52401401316455*^9}]
}, Open ]],
Cell[TextData[{
@@ -491,7 +489,8 @@ Cell[BoxData[
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
Function[BoxForm`e$,
- MatrixForm[BoxForm`e$]]]], "Output"]
+ MatrixForm[BoxForm`e$]]]], "Output",
+ CellChangeTimes->{3.524014013274888*^9}]
}, Open ]],
Cell[TextData[{
@@ -577,7 +576,8 @@ Cell[BoxData[
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
Function[BoxForm`e$,
- MatrixForm[BoxForm`e$]]]], "Output"]
+ MatrixForm[BoxForm`e$]]]], "Output",
+ CellChangeTimes->{3.524014013584079*^9}]
}, Open ]],
Cell["Here are the evolution equations.", "MathCaption",
@@ -601,9 +601,8 @@ Cell[BoxData[
RowBox[{
RowBox[{
SuperscriptBox[
- RowBox[{"(",
- SubscriptBox["\[Rho]",
- RowBox[{"1", ",", "1"}]], ")"}], "\[Prime]",
+ SubscriptBox["\[Rho]",
+ RowBox[{"1", ",", "1"}]], "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{
FractionBox["\[Gamma]t", "2"], "-",
@@ -641,9 +640,8 @@ Cell[BoxData[
RowBox[{
RowBox[{
SuperscriptBox[
- RowBox[{"(",
- SubscriptBox["\[Rho]",
- RowBox[{"1", ",", "2"}]], ")"}], "\[Prime]",
+ SubscriptBox["\[Rho]",
+ RowBox[{"1", ",", "2"}]], "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{
@@ -701,9 +699,8 @@ Cell[BoxData[
RowBox[{
RowBox[{
SuperscriptBox[
- RowBox[{"(",
- SubscriptBox["\[Rho]",
- RowBox[{"1", ",", "3"}]], ")"}], "\[Prime]",
+ SubscriptBox["\[Rho]",
+ RowBox[{"1", ",", "3"}]], "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{
@@ -771,9 +768,8 @@ Cell[BoxData[
RowBox[{
RowBox[{
SuperscriptBox[
- RowBox[{"(",
- SubscriptBox["\[Rho]",
- RowBox[{"2", ",", "1"}]], ")"}], "\[Prime]",
+ SubscriptBox["\[Rho]",
+ RowBox[{"2", ",", "1"}]], "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{
@@ -830,9 +826,8 @@ Cell[BoxData[
RowBox[{
RowBox[{
SuperscriptBox[
- RowBox[{"(",
- SubscriptBox["\[Rho]",
- RowBox[{"2", ",", "2"}]], ")"}], "\[Prime]",
+ SubscriptBox["\[Rho]",
+ RowBox[{"2", ",", "2"}]], "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{
FractionBox["\[Gamma]t", "2"], "-",
@@ -897,9 +892,8 @@ Cell[BoxData[
RowBox[{
RowBox[{
SuperscriptBox[
- RowBox[{"(",
- SubscriptBox["\[Rho]",
- RowBox[{"2", ",", "3"}]], ")"}], "\[Prime]",
+ SubscriptBox["\[Rho]",
+ RowBox[{"2", ",", "3"}]], "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{
@@ -977,9 +971,8 @@ Cell[BoxData[
RowBox[{
RowBox[{
SuperscriptBox[
- RowBox[{"(",
- SubscriptBox["\[Rho]",
- RowBox[{"3", ",", "1"}]], ")"}], "\[Prime]",
+ SubscriptBox["\[Rho]",
+ RowBox[{"3", ",", "1"}]], "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{
@@ -1044,9 +1037,8 @@ Cell[BoxData[
RowBox[{
RowBox[{
SuperscriptBox[
- RowBox[{"(",
- SubscriptBox["\[Rho]",
- RowBox[{"3", ",", "2"}]], ")"}], "\[Prime]",
+ SubscriptBox["\[Rho]",
+ RowBox[{"3", ",", "2"}]], "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{
@@ -1122,9 +1114,8 @@ Cell[BoxData[
RowBox[{
RowBox[{
SuperscriptBox[
- RowBox[{"(",
- SubscriptBox["\[Rho]",
- RowBox[{"3", ",", "3"}]], ")"}], "\[Prime]",
+ SubscriptBox["\[Rho]",
+ RowBox[{"3", ",", "3"}]], "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{"\[Gamma]p", " ",
@@ -1217,7 +1208,8 @@ Cell[BoxData[
Offset[0.2]}, "RowsIndexed" -> {}}],
Column],
Function[BoxForm`e$,
- TableForm[BoxForm`e$]]]], "Output"]
+ TableForm[BoxForm`e$]]]], "Output",
+ CellChangeTimes->{3.524014013741555*^9}]
}, Open ]],
Cell["Convert to c form.", "Text"],
@@ -1252,19 +1244,19 @@ Cell[BoxData[{
FractionBox[
RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Delta]"}], "2"], "+",
RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]}]], " ", "\[Rule]",
- RowBox[{"2",
- RowBox[{"d", "[", "t", "]"}],
+ RowBox[{"2", "d", " ",
RowBox[{"E2", "/", "\[CapitalOmega]2"}]}]}], ",",
RowBox[{
SuperscriptBox["\[ExponentialE]",
RowBox[{
RowBox[{"-",
FractionBox[
- RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Delta]"}], "2"]}],
- "-",
+ RowBox[{GridBox[{
+ {"\[ImaginaryI]"},
+ {"\[Placeholder]"}
+ }], " ", "t", " ", "\[Delta]"}], "2"]}], "-",
RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]}]], " ", "\[Rule]",
- RowBox[{"2",
- RowBox[{"dc", "[", "t", "]"}],
+ RowBox[{"2", "dc", " ",
RowBox[{"E2c", "/", "\[CapitalOmega]2"}]}]}], ",",
RowBox[{
SuperscriptBox["\[ExponentialE]",
@@ -1274,8 +1266,7 @@ Cell[BoxData[{
RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Delta]"}], "2"]}],
"+",
RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]}]], " ", "\[Rule]",
- RowBox[{"2",
- RowBox[{"dc", "[", "t", "]"}],
+ RowBox[{"2", "d", " ",
RowBox[{"E3", "/", "\[CapitalOmega]3"}]}]}], ",",
RowBox[{
SuperscriptBox["\[ExponentialE]",
@@ -1283,8 +1274,7 @@ Cell[BoxData[{
FractionBox[
RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Delta]"}], "2"], "-",
RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]}]], " ", "\[Rule]",
- RowBox[{"2",
- RowBox[{"d", "[", "t", "]"}],
+ RowBox[{"2", "d", " ",
RowBox[{"E3c", "/", "\[CapitalOmega]3"}]}]}], ",",
RowBox[{"\[CapitalGamma]", "\[Rule]",
RowBox[{"2", "G"}]}], ",",
@@ -1349,7 +1339,8 @@ Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{"NotebookDirectory", "[", "]"}], "<>", "\"\<\\\\code.txt\>\""}],
- ",", "%"}], "]"}]}], "Input"],
+ ",", "%"}], "]"}]}], "Input",
+ CellChangeTimes->{{3.524014168318809*^9, 3.524014208362378*^9}}],
Cell[BoxData[
RowBox[{"{",
@@ -1385,13 +1376,9 @@ Cell[BoxData[
RowBox[{"da", " ", "i"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"1", ",", "2"}]]}], "-",
- RowBox[{"i", " ",
+ RowBox[{"d", " ",
RowBox[{"(",
- RowBox[{
- RowBox[{"E2", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E3", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"1", ",", "3"}]]}], "+",
RowBox[{"E1c", " ", "i", " ",
@@ -1405,13 +1392,21 @@ Cell[BoxData[
RowBox[{"-", "E1c"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"1", ",", "1"}]]}], "-",
- RowBox[{"i", " ",
+ RowBox[{
+ FractionBox["1", "2"], " ", "i", " ",
RowBox[{"(",
RowBox[{
- RowBox[{"E3c", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E2c", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"2", " ", "d", " ", "E3c"}], "+",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "i", " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"t", " ", "\[Delta]"}], "+",
+ RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ",
+ "\[CapitalOmega]2"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"1", ",", "2"}]]}], "+",
RowBox[{
@@ -1441,13 +1436,21 @@ Cell[BoxData[
RowBox[{"E1", " ", "i", " ",
SubscriptBox["r",
RowBox[{"2", ",", "3"}]]}], "+",
- RowBox[{"i", " ",
+ RowBox[{
+ FractionBox["1", "2"], " ", "i", " ",
RowBox[{"(",
RowBox[{
- RowBox[{"E3c", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E2c", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"2", " ", "d", " ", "E3c"}], "+",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "i", " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"t", " ", "\[Delta]"}], "+",
+ RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ",
+ "\[CapitalOmega]2"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"3", ",", "1"}]]}]}]}], ",",
RowBox[{
@@ -1457,22 +1460,26 @@ Cell[BoxData[
RowBox[{"2", " ", "gt", " ",
SubscriptBox["r",
RowBox[{"2", ",", "2"}]]}], "-",
- RowBox[{"i", " ",
+ RowBox[{"d", " ",
RowBox[{"(",
- RowBox[{
- RowBox[{"E2", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E3", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"2", ",", "3"}]]}], "+",
- RowBox[{"i", " ",
+ RowBox[{
+ FractionBox["1", "2"], " ", "i", " ",
RowBox[{"(",
RowBox[{
- RowBox[{"E3c", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E2c", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"2", " ", "d", " ", "E3c"}], "+",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "i", " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"t", " ", "\[Delta]"}], "+",
+ RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ",
+ "\[CapitalOmega]2"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"3", ",", "2"}]]}], "+",
RowBox[{"G", " ",
@@ -1486,13 +1493,21 @@ Cell[BoxData[
RowBox[{"-", "E1c"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"2", ",", "1"}]]}], "-",
- RowBox[{"i", " ",
+ RowBox[{
+ FractionBox["1", "2"], " ", "i", " ",
RowBox[{"(",
RowBox[{
- RowBox[{"E3c", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E2c", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"2", " ", "d", " ", "E3c"}], "+",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "i", " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"t", " ", "\[Delta]"}], "+",
+ RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ",
+ "\[CapitalOmega]2"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"2", ",", "2"}]]}], "+",
RowBox[{
@@ -1503,13 +1518,21 @@ Cell[BoxData[
RowBox[{"da", " ", "i"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"2", ",", "3"}]]}], "+",
- RowBox[{"i", " ",
+ RowBox[{
+ FractionBox["1", "2"], " ", "i", " ",
RowBox[{"(",
RowBox[{
- RowBox[{"E3c", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E2c", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"2", " ", "d", " ", "E3c"}], "+",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "i", " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"t", " ", "\[Delta]"}], "+",
+ RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ",
+ "\[CapitalOmega]2"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"3", ",", "3"}]]}]}]}], ",",
RowBox[{
@@ -1519,13 +1542,9 @@ Cell[BoxData[
RowBox[{"E1", " ", "i", " ",
SubscriptBox["r",
RowBox[{"1", ",", "1"}]]}], "+",
- RowBox[{"i", " ",
+ RowBox[{"d", " ",
RowBox[{"(",
- RowBox[{
- RowBox[{"E2", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E3", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"2", ",", "1"}]]}], "+",
RowBox[{
@@ -1546,13 +1565,9 @@ Cell[BoxData[
RowBox[{"E1", " ", "i", " ",
SubscriptBox["r",
RowBox[{"1", ",", "2"}]]}], "+",
- RowBox[{"i", " ",
+ RowBox[{"d", " ",
RowBox[{"(",
- RowBox[{
- RowBox[{"E2", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E3", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"2", ",", "2"}]]}], "+",
RowBox[{
@@ -1563,13 +1578,9 @@ Cell[BoxData[
RowBox[{"da", " ", "i"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"3", ",", "2"}]]}], "-",
- RowBox[{"i", " ",
+ RowBox[{"d", " ",
RowBox[{"(",
- RowBox[{
- RowBox[{"E2", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E3", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"3", ",", "3"}]]}]}]}], ",",
RowBox[{
@@ -1582,32 +1593,38 @@ Cell[BoxData[
RowBox[{"E1", " ", "i", " ",
SubscriptBox["r",
RowBox[{"1", ",", "3"}]]}], "+",
- RowBox[{"i", " ",
+ RowBox[{"d", " ",
RowBox[{"(",
- RowBox[{
- RowBox[{"E2", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E3", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"2", ",", "3"}]]}], "-",
RowBox[{"E1c", " ", "i", " ",
SubscriptBox["r",
RowBox[{"3", ",", "1"}]]}], "-",
- RowBox[{"i", " ",
+ RowBox[{
+ FractionBox["1", "2"], " ", "i", " ",
RowBox[{"(",
RowBox[{
- RowBox[{"E3c", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E2c", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"2", " ", "d", " ", "E3c"}], "+",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "i", " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"t", " ", "\[Delta]"}], "+",
+ RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ",
+ "\[CapitalOmega]2"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"3", ",", "2"}]]}], "-",
RowBox[{"2", " ",
RowBox[{"(",
RowBox[{"G", "+", "gt"}], ")"}], " ",
SubscriptBox["r",
- RowBox[{"3", ",", "3"}]]}]}]}]}], "}"}]], "Output"],
+ RowBox[{"3", ",", "3"}]]}]}]}]}], "}"}]], "Output",
+ CellChangeTimes->{3.524014014687035*^9, 3.524014182827821*^9,
+ 3.524014214903372*^9}],
Cell[BoxData[
RowBox[{"{",
@@ -1643,13 +1660,9 @@ Cell[BoxData[
RowBox[{"da", " ", "i"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"1", ",", "2"}]]}], "-",
- RowBox[{"i", " ",
+ RowBox[{"d", " ",
RowBox[{"(",
- RowBox[{
- RowBox[{"E2", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E3", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"1", ",", "3"}]]}], "+",
RowBox[{"E1c", " ", "i", " ",
@@ -1663,13 +1676,21 @@ Cell[BoxData[
RowBox[{"-", "E1c"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"1", ",", "1"}]]}], "-",
- RowBox[{"i", " ",
+ RowBox[{
+ FractionBox["1", "2"], " ", "i", " ",
RowBox[{"(",
RowBox[{
- RowBox[{"E3c", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E2c", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"2", " ", "d", " ", "E3c"}], "+",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "i", " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"t", " ", "\[Delta]"}], "+",
+ RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ",
+ "\[CapitalOmega]2"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"1", ",", "2"}]]}], "+",
RowBox[{
@@ -1690,22 +1711,26 @@ Cell[BoxData[
RowBox[{"2", " ", "gt", " ",
SubscriptBox["r",
RowBox[{"2", ",", "2"}]]}], "-",
- RowBox[{"i", " ",
+ RowBox[{"d", " ",
RowBox[{"(",
- RowBox[{
- RowBox[{"E2", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E3", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"2", ",", "3"}]]}], "+",
- RowBox[{"i", " ",
+ RowBox[{
+ FractionBox["1", "2"], " ", "i", " ",
RowBox[{"(",
RowBox[{
- RowBox[{"E3c", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E2c", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"2", " ", "d", " ", "E3c"}], "+",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "i", " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"t", " ", "\[Delta]"}], "+",
+ RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ",
+ "\[CapitalOmega]2"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"3", ",", "2"}]]}], "+",
RowBox[{"G", " ",
@@ -1719,13 +1744,21 @@ Cell[BoxData[
RowBox[{"-", "E1c"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"2", ",", "1"}]]}], "-",
- RowBox[{"i", " ",
+ RowBox[{
+ FractionBox["1", "2"], " ", "i", " ",
RowBox[{"(",
RowBox[{
- RowBox[{"E3c", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E2c", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"2", " ", "d", " ", "E3c"}], "+",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "i", " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"t", " ", "\[Delta]"}], "+",
+ RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ",
+ "\[CapitalOmega]2"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"2", ",", "2"}]]}], "+",
RowBox[{
@@ -1736,13 +1769,21 @@ Cell[BoxData[
RowBox[{"da", " ", "i"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"2", ",", "3"}]]}], "+",
- RowBox[{"i", " ",
+ RowBox[{
+ FractionBox["1", "2"], " ", "i", " ",
RowBox[{"(",
RowBox[{
- RowBox[{"E3c", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E2c", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"2", " ", "d", " ", "E3c"}], "+",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "i", " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"t", " ", "\[Delta]"}], "+",
+ RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ",
+ "\[CapitalOmega]2"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"3", ",", "3"}]]}]}]}], ",",
RowBox[{
@@ -1755,52 +1796,65 @@ Cell[BoxData[
RowBox[{"E1", " ", "i", " ",
SubscriptBox["r",
RowBox[{"1", ",", "3"}]]}], "+",
- RowBox[{"i", " ",
+ RowBox[{"d", " ",
RowBox[{"(",
- RowBox[{
- RowBox[{"E2", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E3", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ",
SubscriptBox["r",
RowBox[{"2", ",", "3"}]]}], "-",
RowBox[{"E1c", " ", "i", " ",
SubscriptBox["r",
RowBox[{"3", ",", "1"}]]}], "-",
- RowBox[{"i", " ",
+ RowBox[{
+ FractionBox["1", "2"], " ", "i", " ",
RowBox[{"(",
RowBox[{
- RowBox[{"E3c", " ",
- RowBox[{"d", "[", "t", "]"}]}], "+",
- RowBox[{"E2c", " ",
- RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ",
+ RowBox[{"2", " ", "d", " ", "E3c"}], "+",
+ RowBox[{
+ SuperscriptBox["\[ExponentialE]",
+ RowBox[{
+ RowBox[{"-",
+ FractionBox["1", "2"]}], " ", "i", " ",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{"t", " ", "\[Delta]"}], "+",
+ RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ",
+ "\[CapitalOmega]2"}]}], ")"}], " ",
SubscriptBox["r",
RowBox[{"3", ",", "2"}]]}], "-",
RowBox[{"2", " ",
RowBox[{"(",
RowBox[{"G", "+", "gt"}], ")"}], " ",
SubscriptBox["r",
- RowBox[{"3", ",", "3"}]]}]}]}]}], "}"}]], "Output"],
+ RowBox[{"3", ",", "3"}]]}]}]}]}], "}"}]], "Output",
+ CellChangeTimes->{3.524014014687035*^9, 3.524014182827821*^9,
+ 3.524014214922399*^9}],
Cell[BoxData["\<\"dr11_dt = gt - 2*(gp + gt)*r11 - E1*i*r13 + E1c*i*r31 + \
-G*r33;\\ndr12_dt = (-gp - 2*gt - d1*i + da*i)*r12 - i*(E2*d(t) + \
-E3*dc(t))*r13 + E1c*i*r32;\\ndr13_dt = -(E1c*i*r11) - i*(E3c*d(t) + \
-E2c*dc(t))*r12 + (-G - gp - 2*gt - d1*i)*r13 + E1c*i*r33;\\ndr22_dt = gt - \
-2*gt*r22 - i*(E2*d(t) + E3*dc(t))*r23 + i*(E3c*d(t) + E2c*dc(t))*r32 + G*r33;\
-\\ndr23_dt = -(E1c*i*r21) - i*(E3c*d(t) + E2c*dc(t))*r22 + (-G - 2*gt - \
-da*i)*r23 + i*(E3c*d(t) + E2c*dc(t))*r33;\\ndr33_dt = 2*gp*r11 + E1*i*r13 + \
-i*(E2*d(t) + E3*dc(t))*r23 - E1c*i*r31 - i*(E3c*d(t) + E2c*dc(t))*r32 - 2*(G \
-+ gt)*r33;\\n\"\>"], "Output"],
+G*r33;\\ndr12_dt = (-gp - 2*gt - d1*i + da*i)*r12 - d*(E2 + E3)*i*r13 + \
+E1c*i*r32;\\ndr13_dt = -(E1c*i*r11) - (i*(2*d*E3c + \
+\[CapitalOmega]2/Power(E,(i*(t*\[Delta] + 2*\[Phi]2))/2.))*r12)/2. + (-G - gp \
+- 2*gt - d1*i)*r13 + E1c*i*r33;\\ndr22_dt = gt - 2*gt*r22 - d*(E2 + E3)*i*r23 \
++ (i*(2*d*E3c + \[CapitalOmega]2/Power(E,(i*(t*\[Delta] + \
+2*\[Phi]2))/2.))*r32)/2. + G*r33;\\ndr23_dt = -(E1c*i*r21) - (i*(2*d*E3c + \
+\[CapitalOmega]2/Power(E,(i*(t*\[Delta] + 2*\[Phi]2))/2.))*r22)/2. + (-G - \
+2*gt - da*i)*r23 + (i*(2*d*E3c + \[CapitalOmega]2/Power(E,(i*(t*\[Delta] + 2*\
+\[Phi]2))/2.))*r33)/2.;\\ndr33_dt = 2*gp*r11 + E1*i*r13 + d*(E2 + E3)*i*r23 - \
+E1c*i*r31 - (i*(2*d*E3c + \[CapitalOmega]2/Power(E,(i*(t*\[Delta] + \
+2*\[Phi]2))/2.))*r32)/2. - 2*(G + gt)*r33;\\n\"\>"], "Output",
+ CellChangeTimes->{3.524014014687035*^9, 3.524014182827821*^9,
+ 3.524014214936878*^9}],
-Cell[BoxData["\<\"C:\\\\cygwin\\\\home\\\\Simon\\\\Nresonances\\\\xmds2\\\\\
-Shahriar_system\\\\\\\\code.txt\"\>"], "Output"]
+Cell[BoxData["\<\"/mnt/light_huge_archive/home/evmik/src/my_src/Nresonances/\
+xmds2/Shahriar_system/\\\\code.txt\"\>"], "Output",
+ CellChangeTimes->{3.524014014687035*^9, 3.524014182827821*^9,
+ 3.52401421495504*^9}]
}, Open ]]
}, Open ]]
},
-WindowSize->{1010, 875},
-WindowMargins->{{0, Automatic}, {Automatic, 0}},
+WindowSize->{960, 501},
+WindowMargins->{{Automatic, 0}, {24, Automatic}},
ShowSelection->True,
-FrontEndVersion->"7.0 for Microsoft Windows (64-bit) (February 18, 2009)",
+FrontEndVersion->"8.0 for Linux x86 (32-bit) (February 23, 2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)
@@ -1815,70 +1869,70 @@ CellTagsIndex->{}
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
-Cell[567, 22, 36, 0, 71, "Section"],
-Cell[606, 24, 66, 1, 43, "MathCaption",
+Cell[567, 22, 36, 0, 74, "Section"],
+Cell[606, 24, 66, 1, 44, "MathCaption",
CellID->836781195],
-Cell[675, 27, 85, 2, 31, "Input",
+Cell[675, 27, 85, 2, 30, "Input",
CellID->2058623809],
-Cell[763, 31, 1290, 44, 47, "Text",
+Cell[763, 31, 1128, 37, 52, "Text",
CellID->525777075],
-Cell[2056, 77, 112, 3, 43, "MathCaption",
+Cell[1894, 70, 112, 3, 44, "MathCaption",
CellID->429217524],
-Cell[2171, 82, 1057, 27, 112, "Input",
+Cell[2009, 75, 1057, 27, 126, "Input",
CellID->433132487],
-Cell[3231, 111, 125, 3, 43, "MathCaption",
+Cell[3069, 104, 125, 3, 44, "MathCaption",
CellID->133602844],
Cell[CellGroupData[{
-Cell[3381, 118, 1085, 32, 52, "Input"],
-Cell[4469, 152, 1306, 40, 53, "Output"]
+Cell[3219, 111, 1085, 32, 69, "Input"],
+Cell[4307, 145, 1348, 41, 52, "Output"]
}, Open ]],
-Cell[5790, 195, 166, 4, 43, "MathCaption",
+Cell[5670, 189, 166, 4, 44, "MathCaption",
CellID->462076121],
Cell[CellGroupData[{
-Cell[5981, 203, 508, 15, 31, "Input"],
-Cell[6492, 220, 1755, 50, 71, "Output"]
+Cell[5861, 197, 508, 15, 50, "Input"],
+Cell[6372, 214, 1797, 51, 73, "Output"]
}, Open ]],
-Cell[8262, 273, 76, 1, 43, "MathCaption",
+Cell[8184, 268, 76, 1, 44, "MathCaption",
CellID->358620443],
Cell[CellGroupData[{
-Cell[8363, 278, 464, 15, 31, "Input",
+Cell[8285, 273, 530, 16, 30, "Input",
CellID->167259034],
-Cell[8830, 295, 535, 10, 94, "Output"]
+Cell[8818, 291, 585, 11, 94, "Output"]
}, Open ]],
-Cell[9380, 308, 182, 4, 43, "MathCaption",
+Cell[9418, 305, 182, 4, 44, "MathCaption",
CellID->577766068],
Cell[CellGroupData[{
-Cell[9587, 316, 1379, 38, 72, "Input"],
-Cell[10969, 356, 2462, 70, 117, "Output"]
+Cell[9625, 313, 1379, 38, 69, "Input"],
+Cell[11007, 353, 2503, 71, 112, "Output"]
}, Open ]],
-Cell[13446, 429, 456, 13, 43, "MathCaption",
+Cell[13525, 427, 456, 13, 44, "MathCaption",
CellID->610306692],
Cell[CellGroupData[{
-Cell[13927, 446, 933, 25, 57, "Input",
+Cell[14006, 444, 933, 25, 58, "Input",
CellID->645617687],
-Cell[14863, 473, 717, 20, 71, "Output"]
+Cell[14942, 471, 759, 21, 72, "Output"]
}, Open ]],
-Cell[15595, 496, 390, 12, 43, "MathCaption",
+Cell[15716, 495, 390, 12, 44, "MathCaption",
CellID->854192725],
Cell[CellGroupData[{
-Cell[16010, 512, 1070, 29, 59, "Input",
+Cell[16131, 511, 1070, 29, 58, "Input",
CellID->465762594],
-Cell[17083, 543, 1226, 36, 92, "Output"]
+Cell[17204, 542, 1268, 37, 91, "Output"]
}, Open ]],
-Cell[18324, 582, 76, 1, 43, "MathCaption",
+Cell[18487, 582, 76, 1, 44, "MathCaption",
CellID->314466782],
Cell[CellGroupData[{
-Cell[18425, 587, 250, 7, 31, "Input",
+Cell[18588, 587, 250, 7, 30, "Input",
CellID->298399236],
-Cell[18678, 596, 23733, 623, 315, "Output"]
+Cell[18841, 596, 23487, 615, 298, "Output"]
}, Open ]],
-Cell[42426, 1222, 34, 0, 29, "Text"],
+Cell[42343, 1214, 34, 0, 30, "Text"],
Cell[CellGroupData[{
-Cell[42485, 1226, 4646, 125, 233, "Input"],
-Cell[47134, 1353, 7898, 256, 164, "Output"],
-Cell[55035, 1611, 5277, 171, 107, "Output"],
-Cell[60315, 1784, 630, 8, 145, "Output"],
-Cell[60948, 1794, 125, 1, 30, "Output"]
+Cell[42402, 1218, 4657, 124, 255, "Input"],
+Cell[47062, 1344, 8995, 282, 278, "Output"],
+Cell[56060, 1628, 6444, 201, 205, "Output"],
+Cell[62507, 1831, 958, 13, 210, "Output"],
+Cell[63468, 1846, 217, 3, 29, "Output"]
}, Open ]]
}, Open ]]
}
generated by cgit v1.2.3 (git 2.39.1) at 2024-11-16 20:25:42 +0000