From 0ca3e1a4b9212161d1d3fdaa8752e7cf49861413 Mon Sep 17 00:00:00 2001 From: Simon Rochester Date: Sun, 16 Dec 2012 19:10:53 -0800 Subject: Added contour plots of dispersion slope to RbXMDSSetup1.nb. --- mathemathica_fwm/RbXMDSSetup1.nb | 276 ++++++++++++++++++++++++++++++++++++--- 1 file 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[{ + "\"\\"", ",", + "\"\\""}], "}"}]}], ",", + 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[{ + "\"\\"", ",", + "\"\\""}], "}"}]}], ",", + 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" ] -- cgit v1.2.3