(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[      9721,        311]
NotebookOptionsPosition[      8751,        274]
NotebookOutlinePosition[      9093,        289]
CellTagsIndexPosition[      9050,        286]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
double check the integral done in diracPotential.tex, and perform the first \
Fourier transformation.\
\>", "Text",
 CellChangeTimes->{{3.6569733057557535`*^9, 3.6569733169497957`*^9}, {
  3.6569761518663044`*^9, 3.6569761590057106`*^9}}],

Cell[BoxData[{
 RowBox[{"ClearAll", "[", " ", 
  RowBox[{
  "kappa", ",", " ", "k", ",", " ", "m", ",", " ", "a", ",", " ", "t", ",", 
   " ", "ee"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"$Assumptions", " ", "=", " ", 
   RowBox[{
    RowBox[{"kappa", " ", "\[Element]", " ", "Reals"}], "  ", "&&", " ", 
    RowBox[{"kappa", " ", ">", " ", "0"}], " ", "&&", " ", 
    RowBox[{"\[HBar]", " ", "\[Element]", " ", "Reals"}], " ", "&&", " ", 
    RowBox[{"\[HBar]", " ", ">", " ", "0"}], "&&", " ", 
    RowBox[{"k", " ", "\[Element]", " ", "Reals"}], "  ", "&&", " ", 
    RowBox[{"k", " ", ">", " ", "0"}], " ", "&&", " ", 
    RowBox[{"m", " ", "\[Element]", " ", "Reals"}], "  ", "&&", " ", 
    RowBox[{"m", " ", ">", " ", "0"}], " ", "&&", " ", 
    RowBox[{"t", " ", "\[Element]", " ", "Reals"}], "  ", "&&", " ", 
    RowBox[{"t", " ", ">", " ", "0"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.6569729905185766`*^9, 3.656973120631898*^9}, {
  3.6569731523235383`*^9, 3.6569731589795065`*^9}, {3.6569732692144012`*^9, 
  3.656973289265952*^9}, {3.6569734211467867`*^9, 3.656973424536827*^9}, {
  3.6569735480448956`*^9, 3.65697360571821*^9}, {3.6569755088083477`*^9, 
  3.65697552114665*^9}, {3.6569756398723*^9, 3.656975651738186*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Integrate", "[", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"(", " ", 
     RowBox[{
      RowBox[{"Sqrt", "[", "kappa", "]"}], " ", 
      RowBox[{"E", "^", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "kappa"}], " ", 
         RowBox[{"Abs", "[", "x", "]"}]}], " ", ")"}]}]}], ")"}], "^", "2"}], 
   ",", " ", 
   RowBox[{"{", 
    RowBox[{"x", ",", " ", 
     RowBox[{"-", "Infinity"}], ",", " ", "Infinity"}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{3.656973612314127*^9}],

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

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"a", "[", "k_", "]"}], " ", "=", " ", 
  RowBox[{
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{"kappa", "/", "2"}], "/", "Pi"}], "]"}], " ", 
   RowBox[{"Integrate", "[", " ", 
    RowBox[{
     RowBox[{"E", "^", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "I"}], " ", "k", " ", "x"}], " ", "-", " ", 
        RowBox[{"kappa", " ", 
         RowBox[{"Abs", "[", "x", "]"}]}]}], ")"}]}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"x", ",", " ", 
       RowBox[{"-", "Infinity"}], ",", " ", "Infinity"}], "}"}]}], 
    "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ee", "[", "kappa_", "]"}], " ", "=", " ", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"kappa", " ", "\[HBar]"}], ")"}], "^", "2"}]}], "/", 
   RowBox[{"(", 
    RowBox[{"2", "m"}], ")"}]}]}]}], "Input",
 CellChangeTimes->{{3.656973107313674*^9, 3.6569731086338587`*^9}, {
   3.656973431146762*^9, 3.6569735257158575`*^9}, 3.656974724340966*^9, {
   3.656975533113385*^9, 3.656975566220161*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["kappa", 
    RowBox[{"3", "/", "2"}]], " ", 
   SqrtBox[
    FractionBox["2", "\[Pi]"]]}], 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"k", "-", 
     RowBox[{"\[ImaginaryI]", " ", "kappa"}]}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"k", "+", 
     RowBox[{"\[ImaginaryI]", " ", "kappa"}]}], ")"}]}]]], "Output",
 CellChangeTimes->{
  3.6569735654925566`*^9, 3.6569736207759285`*^9, 3.6569747366289806`*^9, {
   3.656975559244889*^9, 3.6569755675013113`*^9}, 3.6569757734084024`*^9}],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox[
   RowBox[{
    SuperscriptBox["kappa", "2"], " ", 
    SuperscriptBox["\[HBar]", "2"]}], 
   RowBox[{"2", " ", "m"}]]}]], "Output",
 CellChangeTimes->{
  3.6569735654925566`*^9, 3.6569736207759285`*^9, 3.6569747366289806`*^9, {
   3.656975559244889*^9, 3.6569755675013113`*^9}, 3.656975773412405*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"a", "[", "k", "]"}], " ", "/.", " ", 
    RowBox[{"kappa", " ", "\[Rule]", " ", 
     RowBox[{"(*", 
      RowBox[{"m", " ", 
       RowBox[{
        SubscriptBox["V", "0"], "/", 
        RowBox[{"\[HBar]", "^", "2"}]}]}], "*)"}], " ", "\[Kappa]"}]}], " ", "//",
    " ", "Simplify"}], " ", "//", " ", 
  "TraditionalForm"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"a", "[", "k", "]"}], " ", "/.", " ", 
    RowBox[{"kappa", " ", "\[Rule]", " ", 
     RowBox[{"m", " ", 
      RowBox[{
       SubscriptBox["V", "0"], "/", 
       RowBox[{"\[HBar]", "^", "2"}]}]}]}]}], " ", "//", " ", "Simplify"}], 
  " ", "//", " ", "TraditionalForm"}]}], "Input",
 CellChangeTimes->{{3.656973634880628*^9, 3.6569736529875293`*^9}, {
  3.6569740436975603`*^9, 3.6569740637816296`*^9}, {3.6569741931870136`*^9, 
  3.656974199790921*^9}, {3.6569749931393147`*^9, 3.656975000892543*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{
    SqrtBox[
     FractionBox["2", "\[Pi]"]], " ", 
    SuperscriptBox["\[Kappa]", 
     RowBox[{"3", "/", "2"}]]}], 
   RowBox[{
    SuperscriptBox["\[Kappa]", "2"], "+", 
    SuperscriptBox["k", "2"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{3.656974742365814*^9, 3.656975001999484*^9}]
}, Open  ]],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{
    SqrtBox[
     FractionBox["2", "\[Pi]"]], "  ", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"m", " ", 
        SubscriptBox["V", "0"]}], ")"}], 
      RowBox[{"3", "/", "2"}]], "/", "\[HBar]"}]}], 
   RowBox[{
    RowBox[{
     SuperscriptBox["k", "2"], " ", 
     SuperscriptBox["\[HBar]", "2"]}], "+", 
    RowBox[{
     SuperscriptBox["m", "2"], " ", 
     RowBox[{
      SubsuperscriptBox["V", "0", "2"], "/", 
      SuperscriptBox["\[HBar]", "2"]}]}]}]], TraditionalForm]], "Input",
 CellChangeTimes->{{3.656975039297426*^9, 3.6569750676022367`*^9}}],

Cell[BoxData[
 FormBox["", TraditionalForm]], "Input",
 CellChangeTimes->{{3.656974551462776*^9, 3.65697457769448*^9}, 
   3.6569747312351356`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"psi", "[", 
   RowBox[{"x_", ",", " ", "t_"}], "]"}], " ", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"1", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{"2", " ", "Pi"}], "]"}]}], ")"}], 
    RowBox[{"Integrate", "[", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"Exp", "[", " ", 
        RowBox[{
         RowBox[{"I", " ", "k", " ", "x"}], " ", "-", " ", 
         RowBox[{"I", " ", "\[HBar]", " ", 
          RowBox[{"k", "^", "2"}], " ", 
          RowBox[{"t", "/", 
           RowBox[{"(", 
            RowBox[{"2", "m"}], ")"}]}]}]}], "]"}], " ", 
       RowBox[{"a", "[", "k", "]"}]}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"k", ",", " ", 
        RowBox[{"-", "Infinity"}], ",", " ", "Infinity"}], "}"}]}], "]"}]}], 
   " ", "//", " ", "Simplify"}]}]], "Input",
 CellChangeTimes->{{3.6569755738359904`*^9, 3.656975623397522*^9}, {
  3.6569756810659313`*^9, 3.656975741537611*^9}, {3.656975862315708*^9, 
  3.656975900748431*^9}}],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.6569757445543156`*^9, 3.6569757963778133`*^9, 
  3.656976711215311*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"psi", "[", 
   RowBox[{"x_", ",", " ", "t_"}], "]"}], " ", "=", 
  RowBox[{"FourierTransform", "[", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"Exp", "[", " ", 
      RowBox[{
       RowBox[{"-", " ", "I"}], " ", "\[HBar]", " ", 
       RowBox[{"k", "^", "2"}], " ", 
       RowBox[{"t", "/", 
        RowBox[{"(", 
         RowBox[{"2", "m"}], ")"}]}]}], "]"}], " ", 
     RowBox[{"a", "[", "k", "]"}]}], ",", " ", "k", ",", "x"}], "]"}], 
  " "}]], "Input",
 CellChangeTimes->{{3.6569767209140234`*^9, 3.6569767446013403`*^9}, {
  3.6569767805287533`*^9, 3.6569768124358835`*^9}}],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.656977795384096*^9}]
}, Open  ]]
},
WindowSize->{707, 636},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 248, 5, 66, "Text"],
Cell[808, 27, 1258, 22, 116, "Input"],
Cell[CellGroupData[{
Cell[2091, 53, 522, 17, 66, "Input"],
Cell[2616, 72, 70, 1, 40, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2723, 78, 1073, 31, 116, "Input"],
Cell[3799, 111, 541, 16, 96, "Output"],
Cell[4343, 129, 345, 9, 63, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4725, 143, 960, 25, 91, "Input"],
Cell[5688, 170, 349, 11, 113, "Output"]
}, Open  ]],
Cell[6052, 184, 636, 21, 109, InheritFromParent],
Cell[6691, 207, 147, 3, 38, InheritFromParent],
Cell[CellGroupData[{
Cell[6863, 214, 1006, 28, 116, "Input"],
Cell[7872, 244, 128, 2, 40, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8037, 251, 618, 17, 66, "Input"],
Cell[8658, 270, 77, 1, 72, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)