dRho/dt uses conj(E) instead of Ec

* this simplifies the logic in xmds, since otherwise I have to manually
add Ec=conj(E) for appropriate fields.
* This might be a bit slower in computation but should be fine

1 files changed, 13 insertions, 5 deletions

diff --git a/mathemathica_fwm/RbXMDSSetup.nb b/mathemathica_fwm/RbXMDSSetup.nb
index 93b0deb..74e4394 100755
--- a/mathemathica_fwm/RbXMDSSetup.nb
+++ b/mathemathica_fwm/RbXMDSSetup.nb
@@ -727,6 +727,7 @@ Cell[BoxData[
           "\"\<\!\(\*SuperscriptBox[\(\[VeryThinSpace]\), \
 \(2\)]\)\!\(\*SubscriptBox[\(S\), \(1/2\)]\)\>\"", ",", 
            RowBox[{"SingleLetterItalics", "\[Rule]", "False"}]}], "]"}], ",", 
+         
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "2.9"}], ",", "0"}], "}"}]}], "]"}], ",", 
@@ -737,6 +738,7 @@ Cell[BoxData[
           "\"\<\!\(\*SuperscriptBox[\(\[VeryThinSpace]\), \
 \(2\)]\)\!\(\*SubscriptBox[\(P\), \(1/2\)]\)\>\"", ",", 
            RowBox[{"SingleLetterItalics", "\[Rule]", "False"}]}], "]"}], ",", 
+         
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "2.9"}], ",", "1.6"}], "}"}]}], "]"}], ",", 
@@ -747,6 +749,7 @@ Cell[BoxData[
           "\"\<\!\(\*SuperscriptBox[\(\[VeryThinSpace]\), \
 \(2\)]\)\!\(\*SubscriptBox[\(P\), \(3/2\)]\)\>\"", ",", 
            RowBox[{"SingleLetterItalics", "\[Rule]", "False"}]}], "]"}], ",", 
+         
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "2.9"}], ",", "2.7"}], "}"}]}], "]"}]}], "}"}]}], ",", 
@@ -1374,6 +1377,7 @@ Cell[BoxData[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "2.9"}], ",", "2.7"}], "}"}]}], "]"}]}], "}"}]}], ",", 
+   
    RowBox[{"ImageSize", "\[Rule]", 
     RowBox[{"3.5", " ", "72"}]}], ",", 
    RowBox[{"ImagePadding", "\[Rule]", 
@@ -1728,13 +1732,14 @@ Cell[BoxData[{
      RowBox[{
       RowBox[{"\"\<Wc(\>\"", "~~", "a_", "~~", "\"\<)\>\""}], 
       "\[RuleDelayed]", 
-      RowBox[{"\"\<E\>\"", "<>", "a", "<>", "\"\<ac\>\""}]}], ",", 
+      RowBox[{"\"\<conj(E\>\"", "<>", "a", "<>", "\"\<a)\>\""}]}], ",", 
      RowBox[{
       RowBox[{"\"\<d(\>\"", "~~", "a_", "~~", "\"\<)\>\""}], "\[RuleDelayed]", 
       RowBox[{"\"\<delta\>\"", "<>", "a"}]}], ",", 
      RowBox[{
       RowBox[{"\"\<k(\>\"", "~~", "a_", "~~", "\"\<)\>\""}], "\[RuleDelayed]", 
-      RowBox[{"\"\<Kvec\>\"", "<>", "a"}]}]}], "}"}]}], "]"}]}], "Input"],
+      RowBox[{"\"\<Kvec\>\"", "<>", "a"}]}]}], "}"}]}], "]"}]}], "Input",
+ CellChangeTimes->{{3.557696286875835*^9, 3.557696298394461*^9}}],
 
 Cell[BoxData[
  RowBox[{"Export", "[", 
@@ -2095,6 +2100,7 @@ Cell[BoxData[
       RowBox[{
        RowBox[{
        "subsys", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}], ",", 
+       
        RowBox[{
         RowBox[{"4", "\[Pi]", " ", "\[ImaginaryI]", " ", 
          RowBox[{"Energy", "[", 
@@ -2104,6 +2110,7 @@ Cell[BoxData[
            "]"}], "]"}], " ", "n0", " ", 
          RowBox[{
          "rmes", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}], " ", 
+         
          RowBox[{
          "polarizationcomponents1", "\[LeftDoubleBracket]", "k", 
           "\[RightDoubleBracket]"}]}], "/.", "delreps1"}]}], "]"}], "]"}], 
@@ -2250,6 +2257,7 @@ Cell[BoxData[{
      RowBox[{"Power", "[", 
       RowBox[{
        RowBox[{"_", "?", "NumericQ"}], ",", "_"}], "]"}]}], "\[RuleDelayed]", 
+    
     RowBox[{"N", "[", "a", "]"}]}]}], "//", 
   "Simplify"}], "\[IndentingNewLine]", 
  RowBox[{"StringJoin", "[", 
@@ -2316,11 +2324,11 @@ Cell[BoxData[{
      "\"\<RbPropEquations.txt\>\""}], "]"}], ",", "%"}], "]"}]}], "Input"]
 }, Open  ]]
 },
-WindowSize->{1043, 928},
-WindowMargins->{{-1, Automatic}, {Automatic, 2}},
+WindowSize->{956, 1004},
+WindowMargins->{{0, Automatic}, {26, Automatic}},
 PrivateNotebookOptions->{"FileOutlineCache"->False},
 ShowSelection->True,
-FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (October 6, 2011)",
+FrontEndVersion->"8.0 for Linux x86 (32-bit) (February 23, 2011)",
 StyleDefinitions->"Default.nb"
 ]