(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 17940, 529] NotebookOptionsPosition[ 16674, 483] NotebookOutlinePosition[ 17118, 500] CellTagsIndexPosition[ 17075, 497] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"$Assumptions", "=", RowBox[{"{", RowBox[{ RowBox[{"k", ">", "0"}], ",", " ", RowBox[{"\[Kappa]", ">", "0"}]}], "}"}]}]], "Input", CellChangeTimes->{{3.4797313014238*^9, 3.4797313297866*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"k", ">", "0"}], ",", RowBox[{"\[Kappa]", ">", "0"}]}], "}"}]], "Output", CellChangeTimes->{3.4797313314714003`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"Step", " ", "Potential"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.4797387338026*^9, 3.4797387390442*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Solve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"A", "+", "B"}], "==", "C"}], ",", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"Continuity", " ", "of", " ", "\[Psi]", " ", "at", " ", "x"}], " ", "=", " ", "0"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[ImaginaryI]", " ", "k", " ", "A"}], " ", "-", " ", RowBox[{"\[ImaginaryI]", " ", "k", " ", "B"}]}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"-", "C"}], " ", "\[Kappa]"}]}]}], "}"}], ",", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"Continuity", " ", "of", " ", RowBox[{"d\[Psi]", "/", "dx"}], " ", "at", " ", "x"}], " ", "=", " ", "0"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{"B", ",", "C"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.4797313391466*^9, 3.4797314534166*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"B", "\[Rule]", RowBox[{ RowBox[{"-", "A"}], "+", FractionBox[ RowBox[{"2", " ", "A", " ", "k"}], RowBox[{"k", "+", RowBox[{"\[ImaginaryI]", " ", "\[Kappa]"}]}]]}]}], ",", RowBox[{"C", "\[Rule]", FractionBox[ RowBox[{"2", " ", "A", " ", "k"}], RowBox[{"k", "+", RowBox[{"\[ImaginaryI]", " ", "\[Kappa]"}]}]]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.479731459407*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Simplify", "[", RowBox[{"Together", "[", RowBox[{"Solve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"A", "+", "B"}], "==", "C"}], ",", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{ "Continuity", " ", "of", " ", "\[Psi]", " ", "at", " ", "x"}], " ", "=", " ", "0"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[ImaginaryI]", " ", "k", " ", "A"}], " ", "-", " ", RowBox[{"\[ImaginaryI]", " ", "k", " ", "B"}]}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"-", "C"}], " ", "\[Kappa]"}]}]}], "}"}], ",", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"Continuity", " ", "of", " ", RowBox[{"d\[Psi]", "/", "dx"}], " ", "at", " ", "x"}], " ", "=", " ", "0"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{"B", ",", "C"}], "}"}]}], "]"}], "]"}], "]"}]], "Input", CellChangeTimes->{{3.4797314964414*^9, 3.4797315210894003`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"B", "\[Rule]", FractionBox[ RowBox[{"A", " ", RowBox[{"(", RowBox[{"k", "-", RowBox[{"\[ImaginaryI]", " ", "\[Kappa]"}]}], ")"}]}], RowBox[{"k", "+", RowBox[{"\[ImaginaryI]", " ", "\[Kappa]"}]}]]}], ",", RowBox[{"C", "\[Rule]", FractionBox[ RowBox[{"2", " ", "A", " ", "k"}], RowBox[{"k", "+", RowBox[{"\[ImaginaryI]", " ", "\[Kappa]"}]}]]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.4797315011837997`*^9, 3.4797315216198*^9}}] }, Open ]], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", "Tunneling", " ", "*)"}]}]], "Input", CellChangeTimes->{{3.4797387447226*^9, 3.4797387498394003`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Clear", "[", RowBox[{"B", ",", "D"}], "]"}]], "Input", CellChangeTimes->{{3.4797328727826*^9, 3.4797328905978003`*^9}, { 3.4797330416526003`*^9, 3.4797330420738*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"Clear", "::", "\<\"wrsym\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Symbol \\\\\\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\(D\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" is \ Protected.\\\\\\\"\\\", \\\"MT\\\"]\\) \ \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/message/General/wrsym\\\", \ ButtonNote -> \\\"Clear::wrsym\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{{3.4797328843733997`*^9, 3.4797328913154*^9}, 3.4797330426354*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Solve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"A", "+", "B"}], "\[Equal]", RowBox[{"C", "+", "D"}]}], ",", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"cont", " ", "of", " ", "\[Psi]", " ", "at", " ", "x"}], " ", "=", " ", "0"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[ImaginaryI]", " ", "k", " ", "A"}], " ", "-", " ", RowBox[{"\[ImaginaryI]", " ", "k", " ", "B"}]}], " ", "\[Equal]", " ", RowBox[{ RowBox[{ RowBox[{"-", "\[Kappa]"}], " ", "C"}], "+", RowBox[{"\[Kappa]", " ", "D"}]}]}], ",", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"cont", " ", "of", " ", RowBox[{"d\[Psi]", "/", "dx"}], " ", "at", " ", "x"}], " ", "=", " ", "0"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"C", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[Kappa]"}], " ", "a"}]]}], " ", "+", " ", RowBox[{"D", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[Kappa]", " ", "a"}]]}]}], "==", RowBox[{"F", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "k", " ", "a"}]]}]}], ",", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"cont", " ", "of", " ", "\[Psi]", " ", "at", " ", "x"}], " ", "=", " ", "a"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", "\[Kappa]"}], " ", "C", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[Kappa]"}], " ", "a"}]]}], "+", " ", RowBox[{"\[Kappa]", " ", "D", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[Kappa]", " ", "a"}]]}]}], "==", " ", RowBox[{"\[ImaginaryI]", " ", "k", " ", "F", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "k", " ", "a"}]]}]}]}], "}"}], ",", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"cont", " ", "of", " ", RowBox[{"d\[Psi]", "/", "dx"}], " ", "at", " ", "x"}], " ", "=", " ", "a"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{"B", ",", " ", "C", ",", " ", "D", ",", " ", "F"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.4797329246526003`*^9, 3.4797329635122004`*^9}, { 3.4797330170982*^9, 3.4797332027538*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"F", "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"4", " ", "A", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "a", " ", "k"}], "+", RowBox[{"a", " ", "\[Kappa]"}]}]], " ", "k", " ", "\[Kappa]"}], RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SuperscriptBox["k", "2"]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["k", "2"]}], "-", RowBox[{"2", " ", "k", " ", "\[Kappa]"}], "-", RowBox[{"2", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", "k", " ", "\[Kappa]"}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[Kappa]", "2"]}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["\[Kappa]", "2"]}]}]]}]}], ",", RowBox[{"B", "\[Rule]", FractionBox[ RowBox[{"A", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]]}], ")"}], " ", RowBox[{"(", RowBox[{ SuperscriptBox["k", "2"], "+", SuperscriptBox["\[Kappa]", "2"]}], ")"}]}], RowBox[{ RowBox[{"-", SuperscriptBox["k", "2"]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["k", "2"]}], "+", RowBox[{"2", " ", "\[ImaginaryI]", " ", "k", " ", "\[Kappa]"}], "+", RowBox[{"2", " ", "\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", "k", " ", "\[Kappa]"}], "+", SuperscriptBox["\[Kappa]", "2"], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["\[Kappa]", "2"]}]}]]}], ",", RowBox[{"C", "\[Rule]", FractionBox[ RowBox[{"2", " ", "A", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", "k", " ", RowBox[{"(", RowBox[{"k", "+", RowBox[{"\[ImaginaryI]", " ", "\[Kappa]"}]}], ")"}]}], RowBox[{ RowBox[{"-", SuperscriptBox["k", "2"]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["k", "2"]}], "+", RowBox[{"2", " ", "\[ImaginaryI]", " ", "k", " ", "\[Kappa]"}], "+", RowBox[{"2", " ", "\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", "k", " ", "\[Kappa]"}], "+", SuperscriptBox["\[Kappa]", "2"], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["\[Kappa]", "2"]}]}]]}], ",", RowBox[{"D", "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"2", " ", "A", " ", "k", " ", RowBox[{"(", RowBox[{"k", "-", RowBox[{"\[ImaginaryI]", " ", "\[Kappa]"}]}], ")"}]}], RowBox[{ RowBox[{"-", SuperscriptBox["k", "2"]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["k", "2"]}], "+", RowBox[{"2", " ", "\[ImaginaryI]", " ", "k", " ", "\[Kappa]"}], "+", RowBox[{"2", " ", "\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", "k", " ", "\[Kappa]"}], "+", SuperscriptBox["\[Kappa]", "2"], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["\[Kappa]", "2"]}]}]]}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.479733203487*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"F", "=", RowBox[{"-", FractionBox[ RowBox[{"4", " ", "A", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "a", " ", "k"}], "+", RowBox[{"a", " ", "\[Kappa]"}]}]], " ", "k", " ", "\[Kappa]"}], RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SuperscriptBox["k", "2"]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["k", "2"]}], "-", RowBox[{"2", " ", "k", " ", "\[Kappa]"}], "-", RowBox[{"2", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", "k", " ", "\[Kappa]"}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[Kappa]", "2"]}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["\[Kappa]", "2"]}]}]]}]}]], "Input", CellChangeTimes->{{3.4797333397374*^9, 3.479733350205*^9}}], Cell[BoxData[ RowBox[{"-", FractionBox[ RowBox[{"4", " ", "A", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "a", " ", "k"}], "+", RowBox[{"a", " ", "\[Kappa]"}]}]], " ", "k", " ", "\[Kappa]"}], RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SuperscriptBox["k", "2"]}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["k", "2"]}], "-", RowBox[{"2", " ", "k", " ", "\[Kappa]"}], "-", RowBox[{"2", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", "k", " ", "\[Kappa]"}], "+", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[Kappa]", "2"]}], "-", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]], " ", SuperscriptBox["\[Kappa]", "2"]}]}]]}]], "Output", CellChangeTimes->{3.4797333511878*^9}] }, Open ]], Cell[BoxData[ RowBox[{"Simplify", "[", RowBox[{"F", "/", "A"}], "]"}]], "Input", CellChangeTimes->{{3.4797333564762*^9, 3.4797333673806*^9}}], Cell[BoxData[ RowBox[{ FractionBox[ RowBox[{"4", " ", "\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"a", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "k"}], "+", "\[Kappa]"}], ")"}]}]], " ", "k", " ", "\[Kappa]"}], RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]]}], ")"}], " ", SuperscriptBox["k", "2"]}], "+", RowBox[{"2", " ", "\[ImaginaryI]", " ", RowBox[{"(", RowBox[{"1", "+", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]]}], ")"}], " ", "k", " ", "\[Kappa]"}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "a", " ", "\[Kappa]"}]]}], ")"}], " ", SuperscriptBox["\[Kappa]", "2"]}]}]], "\[AliasDelimiter]"}]], "Input", CellChangeTimes->{{3.4797374722692003`*^9, 3.4797374723004*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"FullSimplify", " ", "[", RowBox[{"A", "/", "F"}], "]"}]], "Input", CellChangeTimes->{{3.4797333779574003`*^9, 3.4797333890646*^9}}], Cell[BoxData[ FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", "a", " ", "k"}]], " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "k", " ", "\[Kappa]", " ", RowBox[{"Cosh", "[", RowBox[{"a", " ", "\[Kappa]"}], "]"}]}], "-", RowBox[{"\[ImaginaryI]", " ", RowBox[{"(", RowBox[{"k", "-", "\[Kappa]"}], ")"}], " ", RowBox[{"(", RowBox[{"k", "+", "\[Kappa]"}], ")"}], " ", RowBox[{"Sinh", "[", RowBox[{"a", " ", "\[Kappa]"}], "]"}]}]}], ")"}]}], RowBox[{"2", " ", "k", " ", "\[Kappa]"}]]], "Output", CellChangeTimes->{{3.4797333840726*^9, 3.4797333912642*^9}}] }, Open ]] }, WindowSize->{1008, 647}, WindowMargins->{{0, Automatic}, {Automatic, 31}}, Magnification:>FEPrivate`If[ FEPrivate`Equal[FEPrivate`$VersionNumber, 6.], 1.5, 1.5 Inherited], FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (February 18, 2009)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[567, 22, 231, 6, 43, "Input"], Cell[801, 30, 172, 5, 42, "Output"] }, Open ]], Cell[988, 38, 157, 3, 43, "Input"], Cell[CellGroupData[{ Cell[1170, 45, 934, 24, 98, "Input"], Cell[2107, 71, 523, 17, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2667, 93, 1072, 28, 98, "Input"], Cell[3742, 123, 594, 18, 64, "Output"] }, Open ]], Cell[4351, 144, 194, 3, 98, "Input"], Cell[CellGroupData[{ Cell[4570, 151, 196, 4, 43, "Input"], Cell[4769, 157, 613, 12, 35, "Message"] }, Open ]], Cell[CellGroupData[{ Cell[5419, 174, 2476, 63, 160, "Input"], Cell[7898, 239, 4362, 110, 235, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12297, 354, 1137, 29, 76, "Input"], Cell[13437, 385, 1068, 28, 69, "Output"] }, Open ]], Cell[14520, 416, 146, 3, 43, "Input"], Cell[14669, 421, 1119, 32, 82, InheritFromParent], Cell[CellGroupData[{ Cell[15813, 457, 159, 3, 43, "Input"], Cell[15975, 462, 683, 18, 68, "Output"] }, Open ]] } ] *) (* End of internal cache information *)