summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
Diffstat
-rw-r--r--mathemathica_fwm/RbXMDSSetup1.nb276
1 files changed, 255 insertions, 21 deletions
diff --git a/mathemathica_fwm/RbXMDSSetup1.nb b/mathemathica_fwm/RbXMDSSetup1.nb
index 8c6e618..58a451b 100644
--- a/mathemathica_fwm/RbXMDSSetup1.nb
+++ b/mathemathica_fwm/RbXMDSSetup1.nb
@@ -749,7 +749,6 @@ Cell[BoxData[
"\"\<\!\(\*SuperscriptBox[\(\[VeryThinSpace]\), \
\(2\)]\)\!\(\*SubscriptBox[\(S\), \(1/2\)]\)\>\"", ",",
RowBox[{"SingleLetterItalics", "\[Rule]", "False"}]}], "]"}], ",",
-
RowBox[{"{",
RowBox[{
RowBox[{"-", "2.9"}], ",", "0"}], "}"}]}], "]"}], ",",
@@ -760,7 +759,6 @@ Cell[BoxData[
"\"\<\!\(\*SuperscriptBox[\(\[VeryThinSpace]\), \
\(2\)]\)\!\(\*SubscriptBox[\(P\), \(1/2\)]\)\>\"", ",",
RowBox[{"SingleLetterItalics", "\[Rule]", "False"}]}], "]"}], ",",
-
RowBox[{"{",
RowBox[{
RowBox[{"-", "2.9"}], ",", "1.6"}], "}"}]}], "]"}], ",",
@@ -771,7 +769,6 @@ Cell[BoxData[
"\"\<\!\(\*SuperscriptBox[\(\[VeryThinSpace]\), \
\(2\)]\)\!\(\*SubscriptBox[\(P\), \(3/2\)]\)\>\"", ",",
RowBox[{"SingleLetterItalics", "\[Rule]", "False"}]}], "]"}], ",",
-
RowBox[{"{",
RowBox[{
RowBox[{"-", "2.9"}], ",", "2.7"}], "}"}]}], "]"}]}], "}"}]}], ",",
@@ -979,10 +976,8 @@ Cell[BoxData[
RowBox[{"\[Gamma]t", "\[Rule]", "10."}], ",",
RowBox[{
SubscriptBox["\[CapitalOmega]L", "x"], "\[Rule]", ".1"}], ",",
-
RowBox[{
SubscriptBox["\[CapitalOmega]L", "y"], "\[Rule]", ".2"}], ",",
-
RowBox[{
SubscriptBox["\[CapitalOmega]L", "z"], "\[Rule]", ".3"}]}],
"}"}]}], "/.", "atomicdata"}], "/.",
@@ -1286,7 +1281,6 @@ Cell[BoxData[
SuperscriptBox[
RowBox[{
SubscriptBox["\[CapitalOmega]", "Im"], "[", "z", "]"}], "2"], "+",
-
SuperscriptBox[
RowBox[{
SubscriptBox["\[CapitalOmega]", "Re"], "[", "z", "]"}], "2"]}]],
@@ -1370,7 +1364,7 @@ Cell[BoxData[{
RowBox[{
RowBox[{
SubscriptBox["\[CapitalOmega]",
- RowBox[{"1", ",", "p"}]], "\[Rule]", "100."}], ",",
+ RowBox[{"1", ",", "m"}]], "\[Rule]", "100."}], ",",
RowBox[{
SubscriptBox["\[CapitalOmega]",
RowBox[{"2", ",", "p"}]], "\[Rule]", ".01"}], ",",
@@ -1405,11 +1399,10 @@ Cell[BoxData[{
RowBox[{
RowBox[{"{",
RowBox[{"d", ",", "obs"}], "}"}], "/.",
- RowBox[{"Chop", "@",
- RowBox[{
- RowBox[{"NSolve", "[",
- RowBox[{"steadyeqs1", ",", "vars"}], "]"}], "\[LeftDoubleBracket]",
- "1", "\[RightDoubleBracket]"}]}]}], "/.", "params"}], ",",
+ RowBox[{
+ RowBox[{"NSolve", "[",
+ RowBox[{"steadyeqs1", ",", "vars"}], "]"}], "\[LeftDoubleBracket]",
+ "1", "\[RightDoubleBracket]"}]}], "/.", "params"}], ",",
RowBox[{"{",
RowBox[{"d", ",",
RowBox[{"-", "200"}], ",", "200", ",", "2"}], "}"}]}], "]"}]}],
@@ -1565,7 +1558,6 @@ Cell[BoxData[{
RowBox[{"{",
RowBox[{
RowBox[{"-", "2.9"}], ",", "2.7"}], "}"}]}], "]"}]}], "}"}]}], ",",
-
RowBox[{"ImageSize", "\[Rule]",
RowBox[{"3.5", " ", "72"}]}], ",",
RowBox[{"ImagePadding", "\[Rule]",
@@ -1688,7 +1680,6 @@ Cell[BoxData[{
RowBox[{"{",
RowBox[{
RowBox[{"-", "2.9"}], ",", "2.7"}], "}"}]}], "]"}]}], "}"}]}], ",",
-
RowBox[{"ImageSize", "\[Rule]",
RowBox[{"3.5", " ", "72"}]}], ",",
RowBox[{"ImagePadding", "\[Rule]",
@@ -1699,7 +1690,251 @@ Cell[BoxData[{
RowBox[{"{",
RowBox[{"10", ",", "10"}], "}"}]}], "}"}]}], ",",
RowBox[{"PopulationStyle", "\[Rule]",
- RowBox[{"PointSize", "[", ".018", "]"}]}]}], "]"}]}], "Input"]
+ RowBox[{"PointSize", "[", ".018", "]"}]}]}], "]"}]}], "Input"],
+
+Cell[BoxData[{
+ RowBox[{
+ RowBox[{"table", "=",
+ RowBox[{"Flatten", "[",
+ RowBox[{
+ RowBox[{"Table", "[", "\[IndentingNewLine]",
+ RowBox[{
+ RowBox[{
+ RowBox[{"params", "=",
+ RowBox[{"{",
+ RowBox[{
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]",
+ RowBox[{"1", ",", "m"}]], "\[Rule]",
+ SuperscriptBox["10", "e1m"]}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]",
+ RowBox[{"2", ",", "p"}]], "\[Rule]", ".001"}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]",
+ RowBox[{"3", ",", "m"}]], "\[Rule]",
+ SuperscriptBox["10", "e3m"]}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]",
+ RowBox[{"_", ",", "_"}]], "\[Rule]", "0."}], ",",
+ RowBox[{
+ SubscriptBox["\[Delta]", "2"], "\[Rule]", "d"}], ",",
+ RowBox[{
+ SubscriptBox["\[Delta]", "1"], "\[Rule]", "0"}], ",",
+ RowBox[{
+ SubscriptBox["\[Delta]", "3"], "\[Rule]", "0"}], ",",
+ RowBox[{"\[Gamma]t", "\[Rule]", ".1"}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]L", "x"], "\[Rule]", "0"}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]L", "y"], "\[Rule]", "0"}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]L", "z"], "\[Rule]", "0"}]}],
+ "}"}]}], ";", "\[IndentingNewLine]",
+ RowBox[{"steadyeqs1", "=",
+ RowBox[{"steadyeqs", "/.", "params"}]}], ";", "\[IndentingNewLine]",
+ RowBox[{"points", "=",
+ RowBox[{"Table", "[",
+ RowBox[{
+ RowBox[{
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"d", ",", "obs"}], "}"}], "/.",
+ RowBox[{
+ RowBox[{"NSolve", "[",
+ RowBox[{"steadyeqs1", ",", "vars"}], "]"}],
+ "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}], "/.",
+ "params"}], ",",
+ RowBox[{"{",
+ RowBox[{"d", ",",
+ RowBox[{
+ RowBox[{"-", "1"}], " ",
+ SuperscriptBox["10",
+ RowBox[{"-", "4"}]]}], ",",
+ RowBox[{"1", " ",
+ SuperscriptBox["10",
+ RowBox[{"-", "4"}]]}], ",",
+ RowBox[{"2", " ",
+ SuperscriptBox["10",
+ RowBox[{"-", "4"}]]}]}], "}"}]}], "]"}]}], ";",
+ "\[IndentingNewLine]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"e1m", ",", "e3m", ",",
+ RowBox[{
+ RowBox[{"Subtract", "@@",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{
+ RowBox[{"d\[Phi]", "[",
+ RowBox[{
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]",
+ RowBox[{"2", ",", "p"}]], "/.", "params"}], ",", "#"}],
+ "]"}], "&"}], "/@",
+ RowBox[{"points", "[",
+ RowBox[{"[",
+ RowBox[{"All", ",", "2", ",", "1"}], "]"}], "]"}]}], ")"}]}],
+ "/",
+ RowBox[{"Subtract", "@@",
+ RowBox[{"points", "[",
+ RowBox[{"[",
+ RowBox[{"All", ",", "1"}], "]"}], "]"}]}]}]}], "}"}], "/.",
+ "params"}]}], ",",
+ RowBox[{"{",
+ RowBox[{"e1m", ",",
+ RowBox[{"-", "1."}], ",", "3.", ",", ".2"}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"e3m", ",",
+ RowBox[{"-", "1."}], ",", "3.", ",", ".2"}], "}"}]}], "]"}], ",",
+ "1"}], "]"}]}], ";"}], "\n",
+ RowBox[{"ListContourPlot", "[",
+ RowBox[{"table", ",",
+ RowBox[{"PlotRange", "\[Rule]", "All"}], ",",
+ RowBox[{"ContourLabels", "\[Rule]", "All"}], ",",
+ RowBox[{"FrameLabel", "\[Rule]",
+ RowBox[{"{",
+ RowBox[{
+ "\"\<log(\!\(\*SubscriptBox[\(\[CapitalOmega]\), \(1, m\)]\) [Mrad/sec])\
+\>\"", ",",
+ "\"\<log(\!\(\*SubscriptBox[\(\[CapitalOmega]\), \(3, m\)]\) \
+[Mrad/sec])\>\""}], "}"}]}], ",",
+ RowBox[{"BaseStyle", "\[Rule]", "14"}], ",",
+ RowBox[{"Contours", "\[Rule]",
+ RowBox[{"Join", "[",
+ RowBox[{
+ RowBox[{"Range", "[",
+ RowBox[{
+ RowBox[{"-", "5"}], ",", "0", ",", ".5"}], "]"}], ",",
+ RowBox[{"10", "^",
+ RowBox[{"Range", "[",
+ RowBox[{
+ RowBox[{"-", "4."}], ",",
+ RowBox[{"-", "2."}], ",", "1"}], "]"}]}]}], "]"}]}]}],
+ "]"}]}], "Input"],
+
+Cell[BoxData[
+ RowBox[{
+ RowBox[{"table", "=",
+ RowBox[{"Flatten", "[",
+ RowBox[{
+ RowBox[{"Table", "[", "\[IndentingNewLine]",
+ RowBox[{
+ RowBox[{
+ RowBox[{"params", "=",
+ RowBox[{"{",
+ RowBox[{
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]",
+ RowBox[{"1", ",", "m"}]], "\[Rule]",
+ SuperscriptBox["10", "e1m"]}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]",
+ RowBox[{"2", ",", "p"}]], "\[Rule]", ".001"}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]",
+ RowBox[{"3", ",", "m"}]], "\[Rule]",
+ SuperscriptBox["10", "e3m"]}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]",
+ RowBox[{"_", ",", "_"}]], "\[Rule]", "0."}], ",",
+ RowBox[{
+ SubscriptBox["\[Delta]", "2"], "\[Rule]", "d"}], ",",
+ RowBox[{
+ SubscriptBox["\[Delta]", "1"], "\[Rule]", "0"}], ",",
+ RowBox[{
+ SubscriptBox["\[Delta]", "3"], "\[Rule]", "0"}], ",",
+ RowBox[{"\[Gamma]t", "\[Rule]", "1."}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]L", "x"], "\[Rule]", "0"}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]L", "y"], "\[Rule]", "0"}], ",",
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]L", "z"], "\[Rule]", "0"}]}],
+ "}"}]}], ";", "\[IndentingNewLine]",
+ RowBox[{"steadyeqs1", "=",
+ RowBox[{"steadyeqs", "/.", "params"}]}], ";", "\[IndentingNewLine]",
+ RowBox[{"points", "=",
+ RowBox[{"Table", "[",
+ RowBox[{
+ RowBox[{
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"d", ",", "obs"}], "}"}], "/.",
+ RowBox[{
+ RowBox[{"NSolve", "[",
+ RowBox[{"steadyeqs1", ",", "vars"}], "]"}],
+ "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}], "/.",
+ "params"}], ",",
+ RowBox[{"{",
+ RowBox[{"d", ",",
+ RowBox[{
+ RowBox[{"-", "1"}], " ",
+ SuperscriptBox["10",
+ RowBox[{"-", "4"}]]}], ",",
+ RowBox[{"1", " ",
+ SuperscriptBox["10",
+ RowBox[{"-", "4"}]]}], ",",
+ RowBox[{"2", " ",
+ SuperscriptBox["10",
+ RowBox[{"-", "4"}]]}]}], "}"}]}], "]"}]}], ";",
+ "\[IndentingNewLine]",
+ RowBox[{
+ RowBox[{"{",
+ RowBox[{"e1m", ",", "e3m", ",",
+ RowBox[{
+ RowBox[{"Subtract", "@@",
+ RowBox[{"(",
+ RowBox[{
+ RowBox[{
+ RowBox[{"d\[Phi]", "[",
+ RowBox[{
+ RowBox[{
+ SubscriptBox["\[CapitalOmega]",
+ RowBox[{"2", ",", "p"}]], "/.", "params"}], ",", "#"}],
+ "]"}], "&"}], "/@",
+ RowBox[{"points", "[",
+ RowBox[{"[",
+ RowBox[{"All", ",", "2", ",", "1"}], "]"}], "]"}]}], ")"}]}],
+ "/",
+ RowBox[{"Subtract", "@@",
+ RowBox[{"points", "[",
+ RowBox[{"[",
+ RowBox[{"All", ",", "1"}], "]"}], "]"}]}]}]}], "}"}], "/.",
+ "params"}]}], ",",
+ RowBox[{"{",
+ RowBox[{"e1m", ",",
+ RowBox[{"-", "1."}], ",", "3.", ",", ".2"}], "}"}], ",",
+ RowBox[{"{",
+ RowBox[{"e3m", ",",
+ RowBox[{"-", "1."}], ",", "3.", ",", ".2"}], "}"}]}], "]"}], ",",
+ "1"}], "]"}]}], ";"}]], "Input"],
+
+Cell[BoxData[
+ RowBox[{"ListContourPlot", "[",
+ RowBox[{"table", ",",
+ RowBox[{"PlotRange", "\[Rule]", "All"}], ",",
+ RowBox[{"ContourLabels", "\[Rule]", "All"}], ",",
+ RowBox[{"FrameLabel", "\[Rule]",
+ RowBox[{"{",
+ RowBox[{
+ "\"\<log(\!\(\*SubscriptBox[\(\[CapitalOmega]\), \(1, m\)]\) [Mrad/sec])\
+\>\"", ",",
+ "\"\<log(\!\(\*SubscriptBox[\(\[CapitalOmega]\), \(3, m\)]\) \
+[Mrad/sec])\>\""}], "}"}]}], ",",
+ RowBox[{"BaseStyle", "\[Rule]", "14"}], ",",
+ RowBox[{"Contours", "\[Rule]",
+ RowBox[{"Join", "[",
+ RowBox[{
+ RowBox[{"Range", "[",
+ RowBox[{
+ RowBox[{"-", "5"}], ",", "0", ",", ".05"}], "]"}], ",",
+ RowBox[{"10", "^",
+ RowBox[{"Range", "[",
+ RowBox[{
+ RowBox[{"-", "3."}], ",",
+ RowBox[{"-", "2."}], ",", "1"}], "]"}]}]}], "]"}]}]}],
+ "]"}]], "Input"]
}, Open ]],
Cell[CellGroupData[{
@@ -2518,7 +2753,6 @@ Cell[BoxData[
RowBox[{
RowBox[{
"subsys", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}], ",",
-
RowBox[{
RowBox[{"4", "\[Pi]", " ", "\[ImaginaryI]", " ",
RowBox[{"Energy", "[",
@@ -2528,7 +2762,6 @@ Cell[BoxData[
"]"}], "]"}], " ", "n0", " ",
RowBox[{
"rmes", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}], " ",
-
RowBox[{
"polarizationcomponents1", "\[LeftDoubleBracket]", "k",
"\[RightDoubleBracket]"}]}], "/.", "delreps1"}]}], "]"}], "]"}],
@@ -2730,7 +2963,8 @@ Cell[BoxData[{
RowBox[{
SubscriptBox["\[CapitalOmega]",
RowBox[{"j_", ",", "m"}]], "[", "z", "]"}], "\[Rule]",
- RowBox[{"Em", "[", "j", "]"}]}]}], "}"}]}], "/.", "labelrep"}], "/.",
+ RowBox[{"Em", "[", "j", "]"}]}]}], "}"}]}], "/.", "labelrep"}], "/.",
+
RowBox[{
RowBox[{"Complex", "[",
RowBox[{"0", ",", "a_"}], "]"}], "\[Rule]",
@@ -2845,11 +3079,11 @@ Cell[BoxData[{
3.56319258550938*^9}}]
}, Open ]]
},
-WindowSize->{956, 1004},
-WindowMargins->{{Automatic, 0}, {26, Automatic}},
+WindowSize->{838, 902},
+WindowMargins->{{-2, Automatic}, {Automatic, 0}},
PrivateNotebookOptions->{"FileOutlineCache"->False},
ShowSelection->True,
-FrontEndVersion->"8.0 for Linux x86 (32-bit) (February 23, 2011)",
+FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (October 6, 2011)",
StyleDefinitions->"Default.nb"
]
generated by cgit v1.2.3 (git 2.39.1) at 2024-11-16 07:32:37 +0000