Added contour plots of dispersion slope to RbXMDSSetup1.nb.

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"
 ]