trunk/multipole/ChargeTensor.nb

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[         0,          0]
NotebookDataLength[     19818,        604]
NotebookOptionsPosition[     18290,        543]
NotebookOutlinePosition[     18645,        559]
CellTagsIndexPosition[     18602,        556]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"rc", "=."}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Kappa]", "=."}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", 
    RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", 
   RowBox[{
    RowBox[{"Erfc", "[", 
     RowBox[{"\[Kappa]", " ", "r"}], "]"}], "/", "r"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"g", "[", 
    RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "2"}], "*", "\[Kappa]"}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"E", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"r", "^", "2"}], "*", 
          RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*", 
       RowBox[{"Sqrt", "[", "Pi", "]"}], "*", "r"}], ")"}]}], "-", 
    RowBox[{
     RowBox[{"Erfc", "[", 
      RowBox[{"r", "*", "\[Kappa]"}], "]"}], "/", 
     RowBox[{"r", "^", "2"}]}]}]}], " ", ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"h", "[", 
    RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"4", "*", "\[Kappa]"}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"E", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"r", "^", "2"}], "*", 
          RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*", 
       RowBox[{"Sqrt", "[", "Pi", "]"}], "*", 
       RowBox[{"r", "^", "2"}]}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"4", "*", 
       RowBox[{"\[Kappa]", "^", "3"}]}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"E", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"r", "^", "2"}], "*", 
          RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*", 
       RowBox[{"Sqrt", "[", "Pi", "]"}]}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"2", "*", 
       RowBox[{"Erfc", "[", 
        RowBox[{"r", "*", "\[Kappa]"}], "]"}]}], ")"}], "/", 
     RowBox[{"r", "^", "3"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"f0", "[", 
    RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", 
   RowBox[{
    RowBox[{"f", "[", 
     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ", 
    RowBox[{"f", "[", 
     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"r", "-", "rc"}], ")"}], 
     RowBox[{"g", "[", 
      RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"g0", "[", 
    RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", 
   RowBox[{
    RowBox[{"g", "[", 
     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ", 
    RowBox[{"g", "[", 
     RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"h0", "[", 
    RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", 
   RowBox[{"h", "[", 
    RowBox[{"r", ",", "\[Kappa]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vprime", "[", 
    RowBox[{"r_", ",", "method_"}], "]"}], ":=", 
   RowBox[{"Switch", "[", 
    RowBox[{"method", ",", "\[IndentingNewLine]", "BARE", ",", 
     RowBox[{"g", "[", 
      RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", 
     "TSF", ",", " ", 
     RowBox[{"g0", "[", 
      RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", "\[IndentingNewLine]", "SP",
      ",", " ", 
     RowBox[{"g", "[", 
      RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", "\[IndentingNewLine]", 
     "GSF", ",", " ", 
     RowBox[{
      RowBox[{"g", "[", 
       RowBox[{"r", ",", "\[Kappa]"}], "]"}], " ", "-", " ", 
      RowBox[{"g", "[", 
       RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}], "\[IndentingNewLine]", 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vprimeprime", "[", 
    RowBox[{"r_", ",", " ", "method_"}], "]"}], ":=", " ", 
   RowBox[{"Switch", "[", 
    RowBox[{"method", ",", " ", "\[IndentingNewLine]", "BARE", ",", 
     RowBox[{"h", "[", 
      RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", 
     "TSF", ",", " ", 
     RowBox[{"h0", "[", 
      RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", 
     "SP", ",", " ", 
     RowBox[{"h", "[", 
      RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", 
     "GSF", ",", " ", 
     RowBox[{"h", "[", 
      RowBox[{"r", ",", "\[Kappa]"}], "]"}]}], "\[IndentingNewLine]", "]"}]}],
   ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"R", "[", "\[Alpha]_", "]"}], " ", ":=", 
   RowBox[{"r", " ", 
    RowBox[{"Switch", "[", 
     RowBox[{"\[Alpha]", ",", " ", "\[IndentingNewLine]", "1", ",", 
      RowBox[{
       RowBox[{"Sin", "[", "\[Theta]", "]"}], 
       RowBox[{"Cos", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]",
       "2", ",", 
      RowBox[{
       RowBox[{"Sin", "[", "\[Theta]", "]"}], 
       RowBox[{"Sin", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]",
       "3", ",", 
      RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"\[Alpha]_", ",", "\[Beta]_", ",", " ", "method_"}], "]"}], " ", ":=",
    " ", 
   RowBox[{
    RowBox[{
     RowBox[{"KroneckerDelta", "[", 
      RowBox[{"\[Alpha]", ",", "\[Beta]"}], "]"}], "*", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"vprime", "[", 
        RowBox[{"r", ",", " ", "method"}], "]"}], "/", "r"}], ")"}]}], "+", 
    " ", 
    RowBox[{
     RowBox[{"R", "[", "\[Alpha]", "]"}], 
     RowBox[{"R", "[", "\[Beta]", "]"}], 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"vprimeprime", "[", 
         RowBox[{"r", ",", "method"}], "]"}], " ", "-", " ", 
        RowBox[{
         RowBox[{"vprime", "[", 
          RowBox[{"r", ",", "method"}], "]"}], "/", "r"}]}], ")"}], "/", 
      RowBox[{"r", "^", "2"}]}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"TBare", " ", "=", " ", 
   RowBox[{"MatrixForm", "[", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Integrate", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"r", "^", "2"}], 
           RowBox[{"Integrate", "[", 
            RowBox[{
             RowBox[{"Integrate", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"Sin", "[", "\[Theta]", "]"}], 
                RowBox[{"T", "[", 
                 RowBox[{"i", ",", "j", ",", "BARE"}], "]"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], 
             ",", 
             RowBox[{"{", 
              RowBox[{"\[Phi]", ",", "0", ",", 
               RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", 
          RowBox[{"{", 
           RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"TSP", " ", "=", " ", 
   RowBox[{"MatrixForm", "[", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Integrate", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"r", "^", "2"}], 
           RowBox[{"Integrate", "[", 
            RowBox[{
             RowBox[{"Integrate", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"Sin", "[", "\[Theta]", "]"}], 
                RowBox[{"T", "[", 
                 RowBox[{"i", ",", "j", ",", "SP"}], "]"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], 
             ",", 
             RowBox[{"{", 
              RowBox[{"\[Phi]", ",", "0", ",", 
               RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", 
          RowBox[{"{", 
           RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"TGSF", " ", "=", " ", 
   RowBox[{"MatrixForm", "[", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Integrate", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"r", "^", "2"}], 
           RowBox[{"Integrate", "[", 
            RowBox[{
             RowBox[{"Integrate", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"Sin", "[", "\[Theta]", "]"}], 
                RowBox[{"T", "[", 
                 RowBox[{"i", ",", "j", ",", "GSF"}], "]"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], 
             ",", 
             RowBox[{"{", 
              RowBox[{"\[Phi]", ",", "0", ",", 
               RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", 
          RowBox[{"{", 
           RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"TTSF", " ", "=", " ", 
   RowBox[{"MatrixForm", "[", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Integrate", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"r", "^", "2"}], 
           RowBox[{"Integrate", "[", 
            RowBox[{
             RowBox[{"Integrate", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"Sin", "[", "\[Theta]", "]"}], 
                RowBox[{"T", "[", 
                 RowBox[{"i", ",", "j", ",", "TSF"}], "]"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], 
             ",", 
             RowBox[{"{", 
              RowBox[{"\[Phi]", ",", "0", ",", 
               RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", 
          RowBox[{"{", 
           RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"QBARE", " ", "=", 
   RowBox[{"3", 
    RowBox[{
     RowBox[{"Part", "[", 
      RowBox[{"TBare", ",", "1", ",", "1", ",", "1"}], "]"}], "/", 
     RowBox[{"(", 
      RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"QTSP", " ", "=", " ", 
   RowBox[{"3", 
    RowBox[{
     RowBox[{"Part", "[", 
      RowBox[{"TSP", ",", "1", ",", "1", ",", "1"}], "]"}], "/", 
     RowBox[{"(", 
      RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"QGSF", " ", "=", "  ", 
   RowBox[{"3", 
    RowBox[{
     RowBox[{"Part", "[", 
      RowBox[{"TGSF", ",", "1", ",", "1", ",", "1"}], "]"}], "/", 
     RowBox[{"(", 
      RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"QTSF", " ", "=", " ", 
   RowBox[{"3", 
    RowBox[{
     RowBox[{"Part", "[", 
      RowBox[{"TTSF", ",", "1", ",", "1", ",", "1"}], "]"}], "/", 
     RowBox[{"(", 
      RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.613731558054296*^9, 3.613732200209785*^9}, 
   3.6137322781282454`*^9, {3.61373233169027*^9, 3.613732397354877*^9}, {
   3.613732488348338*^9, 3.6137327713906393`*^9}, {3.613732810835518*^9, 
   3.613732821957967*^9}, {3.613732862763562*^9, 3.613732913141625*^9}, {
   3.61373295595298*^9, 3.613733106904414*^9}, {3.613733177027177*^9, 
   3.613733199740512*^9}, {3.61373326352699*^9, 3.6137332727080936`*^9}, {
   3.613733425731758*^9, 3.6137335719815893`*^9}, {3.613733605778666*^9, 
   3.613733648768111*^9}, {3.613733776819376*^9, 3.6137338150569887`*^9}, {
   3.6137338699918118`*^9, 3.61373393057332*^9}, {3.613733985426153*^9, 
   3.613733985628289*^9}, {3.613734054134171*^9, 3.6137340588514023`*^9}, {
   3.6137371511302147`*^9, 3.613737188118826*^9}, {3.613815304031159*^9, 
   3.613815306899302*^9}, {3.613815415274508*^9, 3.613815450160325*^9}, {
   3.613815535717202*^9, 3.613815604127789*^9}, {3.613815665097225*^9, 
   3.613815673383039*^9}, {3.613816023942026*^9, 3.6138160527230864`*^9}, {
   3.613825536655637*^9, 3.613825555245915*^9}, {3.613825660601698*^9, 
   3.6138256692518787`*^9}, {3.613825968874196*^9, 3.613826000721004*^9}, {
   3.6138263115515537`*^9, 3.613826312408416*^9}, {3.613826919678602*^9, 
   3.613826920319913*^9}, {3.613827178556149*^9, 3.6138272197490597`*^9}, {
   3.6138410090125504`*^9, 3.613841090576209*^9}, {3.613841142597171*^9, 
   3.613841191624243*^9}, {3.6138412274048777`*^9, 3.613841229484283*^9}, {
   3.6138413897003736`*^9, 3.613841438224633*^9}, {3.61384150925097*^9, 
   3.6138415133903227`*^9}, {3.613841547354535*^9, 3.613841573983226*^9}, {
   3.613898318732862*^9, 3.613898326606717*^9}, {3.614594975804741*^9, 
   3.614595013097396*^9}, {3.61459504604961*^9, 3.614595067267371*^9}}],

Cell[CellGroupData[{

Cell[BoxData["QBARE"], "Input",
 CellChangeTimes->{{3.61389845863691*^9, 3.6138984602794447`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"3", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       FractionBox["8", "3"]}], " ", 
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{
        RowBox[{"-", 
         SuperscriptBox["rc", "2"]}], " ", 
        SuperscriptBox["\[Kappa]", "2"]}]], " ", 
      SqrtBox["\[Pi]"], " ", "rc", " ", "\[Kappa]"}], "+", 
     RowBox[{
      FractionBox["4", "3"], " ", "\[Pi]", " ", 
      RowBox[{"Erf", "[", 
       RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]}], ")"}]}], 
  RowBox[{"4", " ", "\[Pi]"}]]], "Output",
 CellChangeTimes->{3.61389846058971*^9, 3.61459514866705*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["QTSP"], "Input",
 CellChangeTimes->{{3.613898461774685*^9, 3.61389846255107*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"3", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       FractionBox["8", "3"]}], " ", 
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{
        RowBox[{"-", 
         SuperscriptBox["rc", "2"]}], " ", 
        SuperscriptBox["\[Kappa]", "2"]}]], " ", 
      SqrtBox["\[Pi]"], " ", "rc", " ", "\[Kappa]"}], "+", 
     RowBox[{
      FractionBox["4", "3"], " ", "\[Pi]", " ", 
      RowBox[{"Erf", "[", 
       RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]}], ")"}]}], 
  RowBox[{"4", " ", "\[Pi]"}]]], "Output",
 CellChangeTimes->{3.6138984628297167`*^9, 3.614595151448133*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["QGSF"], "Input",
 CellChangeTimes->{{3.613898464699469*^9, 3.613898465498917*^9}}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{
  3.613898465771669*^9, {3.614595155790647*^9, 3.61459516007798*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["QTSF"], "Input",
 CellChangeTimes->{{3.613898466555105*^9, 3.613898468075191*^9}}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{3.613898468428026*^9, 3.614595163440785*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", "QTSP", "]"}]], "Input",
 CellChangeTimes->{{3.613898477301433*^9, 3.613898481802129*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   FractionBox[
    RowBox[{"2", " ", 
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{
       RowBox[{"-", 
        SuperscriptBox["rc", "2"]}], " ", 
       SuperscriptBox["\[Kappa]", "2"]}]], " ", "rc", " ", "\[Kappa]"}], 
    SqrtBox["\[Pi]"]]}], "+", 
  RowBox[{"Erf", "[", 
   RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]], "Output",
 CellChangeTimes->{3.613898482310692*^9, 3.614595171545216*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", "QBARE", "]"}]], "Input",
 CellChangeTimes->{{3.6145951775304317`*^9, 3.6145951844023333`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   FractionBox[
    RowBox[{"2", " ", 
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{
       RowBox[{"-", 
        SuperscriptBox["rc", "2"]}], " ", 
       SuperscriptBox["\[Kappa]", "2"]}]], " ", "rc", " ", "\[Kappa]"}], 
    SqrtBox["\[Pi]"]]}], "+", 
  RowBox[{"Erf", "[", 
   RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]], "Output",
 CellChangeTimes->{3.614595184646675*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Limit", "[", 
  RowBox[{"QBARE", ",", 
   RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.6145952058026037`*^9, 3.6145952197375603`*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{3.6145952203235397`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Limit", "[", 
  RowBox[{"QTSP", ",", 
   RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.614595380663261*^9, 3.614595396297574*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.614595392877947*^9, 3.614595396693039*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Limit", "[", 
  RowBox[{"QGSF", ",", 
   RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.6145953985545464`*^9, 3.614595405611649*^9}}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{3.614595405937716*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Limit", "[", 
  RowBox[{"QTSF", ",", 
   RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.6145954076554337`*^9, 3.614595415455311*^9}}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{3.6145954158135767`*^9}]
}, Open  ]]
},
WindowSize->{884, 1337},
WindowMargins->{{Automatic, 292}, {Automatic, 0}},
FrontEndVersion->"10.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (June 27, \
2014)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[400, 13, 13403, 363, 811, "Input"],
Cell[CellGroupData[{
Cell[13828, 380, 98, 1, 28, "Input"],
Cell[13929, 383, 639, 19, 58, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14605, 407, 95, 1, 28, "Input"],
Cell[14703, 410, 643, 19, 58, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15383, 434, 96, 1, 28, "Input"],
Cell[15482, 437, 118, 2, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15637, 444, 96, 1, 28, "Input"],
Cell[15736, 447, 92, 1, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15865, 453, 134, 2, 28, "Input"],
Cell[16002, 457, 445, 13, 58, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16484, 475, 139, 2, 28, "Input"],
Cell[16626, 479, 423, 13, 58, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17086, 497, 192, 4, 28, "Input"],
Cell[17281, 503, 72, 1, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17390, 509, 187, 4, 28, "Input"],
Cell[17580, 515, 94, 1, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17711, 521, 189, 4, 28, "Input"],
Cell[17903, 527, 70, 1, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18010, 533, 189, 4, 28, "Input"],
Cell[18202, 539, 72, 1, 28, "Output"]
}, Open  ]]
}
]
*)

Revision:	4195
Committed:	Thu Jul 17 18:24:41 2014 UTC (9 years, 11 months ago) by gezelter
File size:	19819 byte(s)
Log Message:	Most recent version including dielectric stuff
#	User	Rev	Content
1	gezelter	4195	(CacheID: 234)
2			(* Internal cache information:
3			NotebookFileLineBreakTest
4			NotebookFileLineBreakTest
5			NotebookDataPosition[ 0, 0]
6			NotebookDataLength[ 19818, 604]
7			NotebookOptionsPosition[ 18290, 543]
8			NotebookOutlinePosition[ 18645, 559]
9			CellTagsIndexPosition[ 18602, 556]
10			WindowFrame->Normal*)
11
12			(* Beginning of Notebook Content *)
13			Notebook[{
14			Cell[BoxData[{
15			RowBox[{
16			RowBox[{"rc", "=."}], ";"}], "\[IndentingNewLine]",
17			RowBox[{
18			RowBox[{"\[Kappa]", "=."}], ";"}], "\[IndentingNewLine]",
19			RowBox[{
20			RowBox[{
21			RowBox[{"f", "[",
22			RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ",
23			RowBox[{
24			RowBox[{"Erfc", "[",
25			RowBox[{"\[Kappa]", " ", "r"}], "]"}], "/", "r"}]}],
26			";"}], "\[IndentingNewLine]",
27			RowBox[{
28			RowBox[{
29			RowBox[{"g", "[",
30			RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ",
31			RowBox[{
32			RowBox[{
33			RowBox[{"(",
34			RowBox[{
35			RowBox[{"-", "2"}], "*", "\[Kappa]"}], ")"}], "/",
36			RowBox[{"(",
37			RowBox[{
38			RowBox[{"E", "^",
39			RowBox[{"(",
40			RowBox[{
41			RowBox[{"r", "^", "2"}], "*",
42			RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
43			RowBox[{"Sqrt", "[", "Pi", "]"}], "*", "r"}], ")"}]}], "-",
44			RowBox[{
45			RowBox[{"Erfc", "[",
46			RowBox[{"r", "*", "\[Kappa]"}], "]"}], "/",
47			RowBox[{"r", "^", "2"}]}]}]}], " ", ";"}], "\[IndentingNewLine]",
48			RowBox[{
49			RowBox[{
50			RowBox[{"h", "[",
51			RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=",
52			RowBox[{
53			RowBox[{
54			RowBox[{"(",
55			RowBox[{"4", "*", "\[Kappa]"}], ")"}], "/",
56			RowBox[{"(",
57			RowBox[{
58			RowBox[{"E", "^",
59			RowBox[{"(",
60			RowBox[{
61			RowBox[{"r", "^", "2"}], "*",
62			RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
63			RowBox[{"Sqrt", "[", "Pi", "]"}], "*",
64			RowBox[{"r", "^", "2"}]}], ")"}]}], "+",
65			RowBox[{
66			RowBox[{"(",
67			RowBox[{"4", "*",
68			RowBox[{"\[Kappa]", "^", "3"}]}], ")"}], "/",
69			RowBox[{"(",
70			RowBox[{
71			RowBox[{"E", "^",
72			RowBox[{"(",
73			RowBox[{
74			RowBox[{"r", "^", "2"}], "*",
75			RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
76			RowBox[{"Sqrt", "[", "Pi", "]"}]}], ")"}]}], "+",
77			RowBox[{
78			RowBox[{"(",
79			RowBox[{"2", "*",
80			RowBox[{"Erfc", "[",
81			RowBox[{"r", "*", "\[Kappa]"}], "]"}]}], ")"}], "/",
82			RowBox[{"r", "^", "3"}]}]}]}], ";"}], "\[IndentingNewLine]",
83			RowBox[{
84			RowBox[{
85			RowBox[{"f0", "[",
86			RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ",
87			RowBox[{
88			RowBox[{"f", "[",
89			RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ",
90			RowBox[{"f", "[",
91			RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-",
92			RowBox[{
93			RowBox[{"(",
94			RowBox[{"r", "-", "rc"}], ")"}],
95			RowBox[{"g", "[",
96			RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}]}], " ",
97			";"}], "\[IndentingNewLine]",
98			RowBox[{
99			RowBox[{
100			RowBox[{"g0", "[",
101			RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ",
102			RowBox[{
103			RowBox[{"g", "[",
104			RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ",
105			RowBox[{"g", "[",
106			RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]",
107			RowBox[{
108			RowBox[{
109			RowBox[{"h0", "[",
110			RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=",
111			RowBox[{"h", "[",
112			RowBox[{"r", ",", "\[Kappa]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
113			RowBox[{
114			RowBox[{
115			RowBox[{"vprime", "[",
116			RowBox[{"r_", ",", "method_"}], "]"}], ":=",
117			RowBox[{"Switch", "[",
118			RowBox[{"method", ",", "\[IndentingNewLine]", "BARE", ",",
119			RowBox[{"g", "[",
120			RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]",
121			"TSF", ",", " ",
122			RowBox[{"g0", "[",
123			RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", "\[IndentingNewLine]", "SP",
124			",", " ",
125			RowBox[{"g", "[",
126			RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", "\[IndentingNewLine]",
127			"GSF", ",", " ",
128			RowBox[{
129			RowBox[{"g", "[",
130			RowBox[{"r", ",", "\[Kappa]"}], "]"}], " ", "-", " ",
131			RowBox[{"g", "[",
132			RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}], "\[IndentingNewLine]",
133			"]"}]}], ";"}], "\[IndentingNewLine]",
134			RowBox[{
135			RowBox[{
136			RowBox[{"vprimeprime", "[",
137			RowBox[{"r_", ",", " ", "method_"}], "]"}], ":=", " ",
138			RowBox[{"Switch", "[",
139			RowBox[{"method", ",", " ", "\[IndentingNewLine]", "BARE", ",",
140			RowBox[{"h", "[",
141			RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]",
142			"TSF", ",", " ",
143			RowBox[{"h0", "[",
144			RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]",
145			"SP", ",", " ",
146			RowBox[{"h", "[",
147			RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]",
148			"GSF", ",", " ",
149			RowBox[{"h", "[",
150			RowBox[{"r", ",", "\[Kappa]"}], "]"}]}], "\[IndentingNewLine]", "]"}]}],
151			";"}], "\[IndentingNewLine]",
152			RowBox[{
153			RowBox[{
154			RowBox[{"R", "[", "\[Alpha]_", "]"}], " ", ":=",
155			RowBox[{"r", " ",
156			RowBox[{"Switch", "[",
157			RowBox[{"\[Alpha]", ",", " ", "\[IndentingNewLine]", "1", ",",
158			RowBox[{
159			RowBox[{"Sin", "[", "\[Theta]", "]"}],
160			RowBox[{"Cos", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]",
161			"2", ",",
162			RowBox[{
163			RowBox[{"Sin", "[", "\[Theta]", "]"}],
164			RowBox[{"Sin", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]",
165			"3", ",",
166			RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "]"}]}]}],
167			";"}], "\[IndentingNewLine]",
168			RowBox[{
169			RowBox[{
170			RowBox[{"T", "[",
171			RowBox[{"\[Alpha]_", ",", "\[Beta]_", ",", " ", "method_"}], "]"}], " ", ":=",
172			" ",
173			RowBox[{
174			RowBox[{
175			RowBox[{"KroneckerDelta", "[",
176			RowBox[{"\[Alpha]", ",", "\[Beta]"}], "]"}], "*", " ",
177			RowBox[{"(",
178			RowBox[{
179			RowBox[{"vprime", "[",
180			RowBox[{"r", ",", " ", "method"}], "]"}], "/", "r"}], ")"}]}], "+",
181			" ",
182			RowBox[{
183			RowBox[{"R", "[", "\[Alpha]", "]"}],
184			RowBox[{"R", "[", "\[Beta]", "]"}],
185			RowBox[{
186			RowBox[{"(",
187			RowBox[{
188			RowBox[{"vprimeprime", "[",
189			RowBox[{"r", ",", "method"}], "]"}], " ", "-", " ",
190			RowBox[{
191			RowBox[{"vprime", "[",
192			RowBox[{"r", ",", "method"}], "]"}], "/", "r"}]}], ")"}], "/",
193			RowBox[{"r", "^", "2"}]}]}]}]}], ";"}], "\[IndentingNewLine]",
194			RowBox[{
195			RowBox[{"TBare", " ", "=", " ",
196			RowBox[{"MatrixForm", "[",
197			RowBox[{"Table", "[",
198			RowBox[{
199			RowBox[{"Table", "[",
200			RowBox[{
201			RowBox[{"Integrate", "[",
202			RowBox[{
203			RowBox[{
204			RowBox[{"r", "^", "2"}],
205			RowBox[{"Integrate", "[",
206			RowBox[{
207			RowBox[{"Integrate", "[",
208			RowBox[{
209			RowBox[{
210			RowBox[{"Sin", "[", "\[Theta]", "]"}],
211			RowBox[{"T", "[",
212			RowBox[{"i", ",", "j", ",", "BARE"}], "]"}]}], ",",
213			RowBox[{"{",
214			RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
215			",",
216			RowBox[{"{",
217			RowBox[{"\[Phi]", ",", "0", ",",
218			RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
219			RowBox[{"{",
220			RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
221			RowBox[{"{",
222			RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
223			RowBox[{"{",
224			RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
225			";"}], "\[IndentingNewLine]",
226			RowBox[{
227			RowBox[{"TSP", " ", "=", " ",
228			RowBox[{"MatrixForm", "[",
229			RowBox[{"Table", "[",
230			RowBox[{
231			RowBox[{"Table", "[",
232			RowBox[{
233			RowBox[{"Integrate", "[",
234			RowBox[{
235			RowBox[{
236			RowBox[{"r", "^", "2"}],
237			RowBox[{"Integrate", "[",
238			RowBox[{
239			RowBox[{"Integrate", "[",
240			RowBox[{
241			RowBox[{
242			RowBox[{"Sin", "[", "\[Theta]", "]"}],
243			RowBox[{"T", "[",
244			RowBox[{"i", ",", "j", ",", "SP"}], "]"}]}], ",",
245			RowBox[{"{",
246			RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
247			",",
248			RowBox[{"{",
249			RowBox[{"\[Phi]", ",", "0", ",",
250			RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
251			RowBox[{"{",
252			RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
253			RowBox[{"{",
254			RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
255			RowBox[{"{",
256			RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
257			";"}], "\[IndentingNewLine]",
258			RowBox[{
259			RowBox[{"TGSF", " ", "=", " ",
260			RowBox[{"MatrixForm", "[",
261			RowBox[{"Table", "[",
262			RowBox[{
263			RowBox[{"Table", "[",
264			RowBox[{
265			RowBox[{"Integrate", "[",
266			RowBox[{
267			RowBox[{
268			RowBox[{"r", "^", "2"}],
269			RowBox[{"Integrate", "[",
270			RowBox[{
271			RowBox[{"Integrate", "[",
272			RowBox[{
273			RowBox[{
274			RowBox[{"Sin", "[", "\[Theta]", "]"}],
275			RowBox[{"T", "[",
276			RowBox[{"i", ",", "j", ",", "GSF"}], "]"}]}], ",",
277			RowBox[{"{",
278			RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
279			",",
280			RowBox[{"{",
281			RowBox[{"\[Phi]", ",", "0", ",",
282			RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
283			RowBox[{"{",
284			RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
285			RowBox[{"{",
286			RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
287			RowBox[{"{",
288			RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
289			";"}], "\[IndentingNewLine]",
290			RowBox[{
291			RowBox[{"TTSF", " ", "=", " ",
292			RowBox[{"MatrixForm", "[",
293			RowBox[{"Table", "[",
294			RowBox[{
295			RowBox[{"Table", "[",
296			RowBox[{
297			RowBox[{"Integrate", "[",
298			RowBox[{
299			RowBox[{
300			RowBox[{"r", "^", "2"}],
301			RowBox[{"Integrate", "[",
302			RowBox[{
303			RowBox[{"Integrate", "[",
304			RowBox[{
305			RowBox[{
306			RowBox[{"Sin", "[", "\[Theta]", "]"}],
307			RowBox[{"T", "[",
308			RowBox[{"i", ",", "j", ",", "TSF"}], "]"}]}], ",",
309			RowBox[{"{",
310			RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
311			",",
312			RowBox[{"{",
313			RowBox[{"\[Phi]", ",", "0", ",",
314			RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
315			RowBox[{"{",
316			RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
317			RowBox[{"{",
318			RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
319			RowBox[{"{",
320			RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
321			";"}], "\[IndentingNewLine]",
322			RowBox[{
323			RowBox[{"QBARE", " ", "=",
324			RowBox[{"3",
325			RowBox[{
326			RowBox[{"Part", "[",
327			RowBox[{"TBare", ",", "1", ",", "1", ",", "1"}], "]"}], "/",
328			RowBox[{"(",
329			RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]",
330			RowBox[{
331			RowBox[{"QTSP", " ", "=", " ",
332			RowBox[{"3",
333			RowBox[{
334			RowBox[{"Part", "[",
335			RowBox[{"TSP", ",", "1", ",", "1", ",", "1"}], "]"}], "/",
336			RowBox[{"(",
337			RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]",
338			RowBox[{
339			RowBox[{"QGSF", " ", "=", " ",
340			RowBox[{"3",
341			RowBox[{
342			RowBox[{"Part", "[",
343			RowBox[{"TGSF", ",", "1", ",", "1", ",", "1"}], "]"}], "/",
344			RowBox[{"(",
345			RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]",
346			RowBox[{
347			RowBox[{"QTSF", " ", "=", " ",
348			RowBox[{"3",
349			RowBox[{
350			RowBox[{"Part", "[",
351			RowBox[{"TTSF", ",", "1", ",", "1", ",", "1"}], "]"}], "/",
352			RowBox[{"(",
353			RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}]}], "Input",
354			CellChangeTimes->{{3.613731558054296^9, 3.613732200209785^9},
355			3.6137322781282454`^9, {3.61373233169027^9, 3.613732397354877*^9}, {
356			3.613732488348338^9, 3.6137327713906393`^9}, {3.613732810835518*^9,
357			3.613732821957967^9}, {3.613732862763562^9, 3.613732913141625*^9}, {
358			3.61373295595298^9, 3.613733106904414^9}, {3.613733177027177*^9,
359			3.613733199740512^9}, {3.61373326352699^9, 3.6137332727080936`*^9}, {
360			3.613733425731758^9, 3.6137335719815893`^9}, {3.613733605778666*^9,
361			3.613733648768111^9}, {3.613733776819376^9, 3.6137338150569887`*^9}, {
362			3.6137338699918118`^9, 3.61373393057332^9}, {3.613733985426153*^9,
363			3.613733985628289^9}, {3.613734054134171^9, 3.6137340588514023`*^9}, {
364			3.6137371511302147`^9, 3.613737188118826^9}, {3.613815304031159*^9,
365			3.613815306899302^9}, {3.613815415274508^9, 3.613815450160325*^9}, {
366			3.613815535717202^9, 3.613815604127789^9}, {3.613815665097225*^9,
367			3.613815673383039^9}, {3.613816023942026^9, 3.6138160527230864`*^9}, {
368			3.613825536655637^9, 3.613825555245915^9}, {3.613825660601698*^9,
369			3.6138256692518787`^9}, {3.613825968874196^9, 3.613826000721004*^9}, {
370			3.6138263115515537`^9, 3.613826312408416^9}, {3.613826919678602*^9,
371			3.613826920319913^9}, {3.613827178556149^9, 3.6138272197490597`*^9}, {
372			3.6138410090125504`^9, 3.613841090576209^9}, {3.613841142597171*^9,
373			3.613841191624243^9}, {3.6138412274048777`^9, 3.613841229484283*^9}, {
374			3.6138413897003736`^9, 3.613841438224633^9}, {3.61384150925097*^9,
375			3.6138415133903227`^9}, {3.613841547354535^9, 3.613841573983226*^9}, {
376			3.613898318732862^9, 3.613898326606717^9}, {3.614594975804741*^9,
377			3.614595013097396^9}, {3.61459504604961^9, 3.614595067267371*^9}}],
378
379			Cell[CellGroupData[{
380
381			Cell[BoxData["QBARE"], "Input",
382			CellChangeTimes->{{3.61389845863691^9, 3.6138984602794447`^9}}],
383
384			Cell[BoxData[
385			FractionBox[
386			RowBox[{"3", " ",
387			RowBox[{"(",
388			RowBox[{
389			RowBox[{
390			RowBox[{"-",
391			FractionBox["8", "3"]}], " ",
392			SuperscriptBox["\[ExponentialE]",
393			RowBox[{
394			RowBox[{"-",
395			SuperscriptBox["rc", "2"]}], " ",
396			SuperscriptBox["\[Kappa]", "2"]}]], " ",
397			SqrtBox["\[Pi]"], " ", "rc", " ", "\[Kappa]"}], "+",
398			RowBox[{
399			FractionBox["4", "3"], " ", "\[Pi]", " ",
400			RowBox[{"Erf", "[",
401			RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]}], ")"}]}],
402			RowBox[{"4", " ", "\[Pi]"}]]], "Output",
403			CellChangeTimes->{3.61389846058971^9, 3.61459514866705^9}]
404			}, Open ]],
405
406			Cell[CellGroupData[{
407
408			Cell[BoxData["QTSP"], "Input",
409			CellChangeTimes->{{3.613898461774685^9, 3.61389846255107^9}}],
410
411			Cell[BoxData[
412			FractionBox[
413			RowBox[{"3", " ",
414			RowBox[{"(",
415			RowBox[{
416			RowBox[{
417			RowBox[{"-",
418			FractionBox["8", "3"]}], " ",
419			SuperscriptBox["\[ExponentialE]",
420			RowBox[{
421			RowBox[{"-",
422			SuperscriptBox["rc", "2"]}], " ",
423			SuperscriptBox["\[Kappa]", "2"]}]], " ",
424			SqrtBox["\[Pi]"], " ", "rc", " ", "\[Kappa]"}], "+",
425			RowBox[{
426			FractionBox["4", "3"], " ", "\[Pi]", " ",
427			RowBox[{"Erf", "[",
428			RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]}], ")"}]}],
429			RowBox[{"4", " ", "\[Pi]"}]]], "Output",
430			CellChangeTimes->{3.6138984628297167`^9, 3.614595151448133^9}]
431			}, Open ]],
432
433			Cell[CellGroupData[{
434
435			Cell[BoxData["QGSF"], "Input",
436			CellChangeTimes->{{3.613898464699469^9, 3.613898465498917^9}}],
437
438			Cell[BoxData["1"], "Output",
439			CellChangeTimes->{
440			3.613898465771669^9, {3.614595155790647^9, 3.61459516007798*^9}}]
441			}, Open ]],
442
443			Cell[CellGroupData[{
444
445			Cell[BoxData["QTSF"], "Input",
446			CellChangeTimes->{{3.613898466555105^9, 3.613898468075191^9}}],
447
448			Cell[BoxData["1"], "Output",
449			CellChangeTimes->{3.613898468428026^9, 3.614595163440785^9}]
450			}, Open ]],
451
452			Cell[CellGroupData[{
453
454			Cell[BoxData[
455			RowBox[{"FullSimplify", "[", "QTSP", "]"}]], "Input",
456			CellChangeTimes->{{3.613898477301433^9, 3.613898481802129^9}}],
457
458			Cell[BoxData[
459			RowBox[{
460			RowBox[{"-",
461			FractionBox[
462			RowBox[{"2", " ",
463			SuperscriptBox["\[ExponentialE]",
464			RowBox[{
465			RowBox[{"-",
466			SuperscriptBox["rc", "2"]}], " ",
467			SuperscriptBox["\[Kappa]", "2"]}]], " ", "rc", " ", "\[Kappa]"}],
468			SqrtBox["\[Pi]"]]}], "+",
469			RowBox[{"Erf", "[",
470			RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]], "Output",
471			CellChangeTimes->{3.613898482310692^9, 3.614595171545216^9}]
472			}, Open ]],
473
474			Cell[CellGroupData[{
475
476			Cell[BoxData[
477			RowBox[{"FullSimplify", "[", "QBARE", "]"}]], "Input",
478			CellChangeTimes->{{3.6145951775304317`^9, 3.6145951844023333`^9}}],
479
480			Cell[BoxData[
481			RowBox[{
482			RowBox[{"-",
483			FractionBox[
484			RowBox[{"2", " ",
485			SuperscriptBox["\[ExponentialE]",
486			RowBox[{
487			RowBox[{"-",
488			SuperscriptBox["rc", "2"]}], " ",
489			SuperscriptBox["\[Kappa]", "2"]}]], " ", "rc", " ", "\[Kappa]"}],
490			SqrtBox["\[Pi]"]]}], "+",
491			RowBox[{"Erf", "[",
492			RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]], "Output",
493			CellChangeTimes->{3.614595184646675*^9}]
494			}, Open ]],
495
496			Cell[CellGroupData[{
497
498			Cell[BoxData[
499			RowBox[{"Limit", "[",
500			RowBox[{"QBARE", ",",
501			RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
502			CellChangeTimes->{{3.6145952058026037`^9, 3.6145952197375603`^9}}],
503
504			Cell[BoxData["0"], "Output",
505			CellChangeTimes->{3.6145952203235397`*^9}]
506			}, Open ]],
507
508			Cell[CellGroupData[{
509
510			Cell[BoxData[
511			RowBox[{"Limit", "[",
512			RowBox[{"QTSP", ",",
513			RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
514			CellChangeTimes->{{3.614595380663261^9, 3.614595396297574^9}}],
515
516			Cell[BoxData["0"], "Output",
517			CellChangeTimes->{{3.614595392877947^9, 3.614595396693039^9}}]
518			}, Open ]],
519
520			Cell[CellGroupData[{
521
522			Cell[BoxData[
523			RowBox[{"Limit", "[",
524			RowBox[{"QGSF", ",",
525			RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
526			CellChangeTimes->{{3.6145953985545464`^9, 3.614595405611649^9}}],
527
528			Cell[BoxData["1"], "Output",
529			CellChangeTimes->{3.614595405937716*^9}]
530			}, Open ]],
531
532			Cell[CellGroupData[{
533
534			Cell[BoxData[
535			RowBox[{"Limit", "[",
536			RowBox[{"QTSF", ",",
537			RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
538			CellChangeTimes->{{3.6145954076554337`^9, 3.614595415455311^9}}],
539
540			Cell[BoxData["1"], "Output",
541			CellChangeTimes->{3.6145954158135767`*^9}]
542			}, Open ]]
543			},
544			WindowSize->{884, 1337},
545			WindowMargins->{{Automatic, 292}, {Automatic, 0}},
546			FrontEndVersion->"10.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (June 27, \
547			2014)",
548			StyleDefinitions->"Default.nb"
549			]
550			(* End of Notebook Content *)
551
552			(* Internal cache information *)
553			(*CellTagsOutline
554			CellTagsIndex->{}
555			*)
556			(*CellTagsIndex
557			CellTagsIndex->{}
558			*)
559			(*NotebookFileOutline
560			Notebook[{
561			Cell[400, 13, 13403, 363, 811, "Input"],
562			Cell[CellGroupData[{
563			Cell[13828, 380, 98, 1, 28, "Input"],
564			Cell[13929, 383, 639, 19, 58, "Output"]
565			}, Open ]],
566			Cell[CellGroupData[{
567			Cell[14605, 407, 95, 1, 28, "Input"],
568			Cell[14703, 410, 643, 19, 58, "Output"]
569			}, Open ]],
570			Cell[CellGroupData[{
571			Cell[15383, 434, 96, 1, 28, "Input"],
572			Cell[15482, 437, 118, 2, 28, "Output"]
573			}, Open ]],
574			Cell[CellGroupData[{
575			Cell[15637, 444, 96, 1, 28, "Input"],
576			Cell[15736, 447, 92, 1, 28, "Output"]
577			}, Open ]],
578			Cell[CellGroupData[{
579			Cell[15865, 453, 134, 2, 28, "Input"],
580			Cell[16002, 457, 445, 13, 58, "Output"]
581			}, Open ]],
582			Cell[CellGroupData[{
583			Cell[16484, 475, 139, 2, 28, "Input"],
584			Cell[16626, 479, 423, 13, 58, "Output"]
585			}, Open ]],
586			Cell[CellGroupData[{
587			Cell[17086, 497, 192, 4, 28, "Input"],
588			Cell[17281, 503, 72, 1, 28, "Output"]
589			}, Open ]],
590			Cell[CellGroupData[{
591			Cell[17390, 509, 187, 4, 28, "Input"],
592			Cell[17580, 515, 94, 1, 28, "Output"]
593			}, Open ]],
594			Cell[CellGroupData[{
595			Cell[17711, 521, 189, 4, 28, "Input"],
596			Cell[17903, 527, 70, 1, 28, "Output"]
597			}, Open ]],
598			Cell[CellGroupData[{
599			Cell[18010, 533, 189, 4, 28, "Input"],
600			Cell[18202, 539, 72, 1, 28, "Output"]
601			}, Open ]]
602			}
603			]
604			*)
605