(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 33579, 763] NotebookOptionsPosition[ 32681, 728] NotebookOutlinePosition[ 33052, 744] CellTagsIndexPosition[ 33009, 741] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["\<\ Disperson relationship across a bandgap \ \>", "Title", CellChangeTimes->{{3.440930658059*^9, 3.4409306911470003`*^9}, { 3.4409308995109997`*^9, 3.440930910046*^9}}], Cell["\<\ From Optical Imaging and Spectroscopy (Wiley - OSA 2009) By David J. Brady www.opticalimaging.org Figure 4.24\ \>", "Subtitle", CellChangeTimes->{{3.4409307282320004`*^9, 3.440930761152*^9}, { 3.440930913613*^9, 3.440930913829*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"K", "=", RowBox[{"2", " ", RowBox[{"Sin", "[", RowBox[{"Pi", "/", "6"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b", "=", ".01"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"hmf", "=", RowBox[{"8", RowBox[{"b", "/", RowBox[{"K", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"mScale", "=", "31"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"a", "[", "kz_", "]"}], ":=", RowBox[{"1", "-", RowBox[{"kz", "^", "2"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"bmf", "[", "kz_", "]"}], ":=", RowBox[{"4", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", "kz", "]"}], "+", "b"}], ")"}], "/", RowBox[{"K", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"H0off", "[", "kz_", "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"xH", "=", RowBox[{"Transpose", "[", RowBox[{"Join", "[", RowBox[{ RowBox[{"{", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"mScale", "+", "1"}], "}"}]}], "]"}], "}"}], ",", RowBox[{"Transpose", "[", RowBox[{"Join", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"DiagonalMatrix", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"hmf", "^", "2"}], "/", RowBox[{"(", RowBox[{ RowBox[{"16", " ", RowBox[{"n", "^", "2"}]}], "+", RowBox[{"2", " ", "hmf"}], "-", RowBox[{"4", RowBox[{"bmf", "[", "kz", "]"}]}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"n", ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"mScale", "-", "1"}], ")"}]}], "/", "2"}], ",", RowBox[{ RowBox[{"(", RowBox[{"mScale", "-", "1"}], ")"}], "/", "2"}]}], "}"}]}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", "mScale", "}"}]}], "]"}], "}"}]}], "]"}], "]"}]}], "]"}], "]"}]}], ";", RowBox[{"xH", "=", RowBox[{"xH", "+", RowBox[{"Transpose", "[", "xH", "]"}]}]}]}], ")"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"H0", "[", "kz_", "]"}], ":=", RowBox[{"N", "[", RowBox[{"Det", "[", RowBox[{ RowBox[{"IdentityMatrix", "[", RowBox[{"mScale", "+", "1"}], "]"}], "+", RowBox[{"H0off", "[", "kz", "]"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "Input"], Cell[BoxData[{ RowBox[{ RowBox[{"K", "=", RowBox[{"2", " ", RowBox[{"Sin", "[", RowBox[{"Pi", "/", "6"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b", "=", RowBox[{".5", " ", RowBox[{"10", "^", RowBox[{"-", "1"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"mScale", "=", "9"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"a", "[", "kz_", "]"}], ":=", RowBox[{"1", "-", RowBox[{"kz", "^", "2"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Hoff", "[", "kz_", "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"xH", "=", RowBox[{"Transpose", "[", RowBox[{"Join", "[", RowBox[{ RowBox[{"{", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", "mScale", "}"}]}], "]"}], "}"}], ",", RowBox[{"Transpose", "[", RowBox[{"Join", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"DiagonalMatrix", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"b", "/", "2"}], ",", RowBox[{"{", RowBox[{"mScale", "-", "1"}], "}"}]}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"mScale", "-", "1"}], "}"}]}], "]"}], "}"}]}], "]"}], "]"}]}], "]"}], "]"}]}], ";", RowBox[{"xH", "=", RowBox[{"xH", "+", RowBox[{"Transpose", "[", "xH", "]"}]}]}]}], ")"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"H", "[", RowBox[{"kz_", ",", "q_"}], "]"}], ":=", RowBox[{"Det", "[", RowBox[{ RowBox[{"DiagonalMatrix", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"a", "[", "kz", "]"}], "-", RowBox[{ RowBox[{"(", RowBox[{"q", "+", " ", RowBox[{"n", " ", "K"}]}], ")"}], "^", "2"}]}], " ", ",", RowBox[{"{", RowBox[{"n", ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"mScale", "-", "1"}], ")"}]}], "/", "2"}], ",", RowBox[{ RowBox[{"(", RowBox[{"mScale", "-", "1"}], ")"}], "/", "2"}]}], "}"}]}], "]"}], "]"}], "+", RowBox[{"Hoff", "[", "kz", "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"qq", "[", "kz_", "]"}], ":=", RowBox[{"qq", "/.", RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{"H", "[", RowBox[{"kz", ",", "qq"}], "]"}], "\[Equal]", "0"}], ",", "qq"}], "]"}]}]}], ";"}], "\[IndentingNewLine]"}], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"aa", "=", RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{"qq", "[", "kz", "]"}], "\[LeftDoubleBracket]", "11", "\[RightDoubleBracket]"}], ",", RowBox[{"{", RowBox[{"kz", ",", RowBox[{ SqrtBox[ RowBox[{"1", "-", FractionBox[ SuperscriptBox["K", "2"], "4"]}]], "-", "0.09"}], ",", SqrtBox[ RowBox[{"1", "-", FractionBox[ SuperscriptBox["K", "2"], "4"]}]]}], "}"}], ",", RowBox[{"Frame", "\[Rule]", "True"}], ",", RowBox[{"Axes", "\[Rule]", "False"}]}], "]"}]}]], "Input", CellChangeTimes->{3.440931754801*^9}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV0n04lXcYB3AvRayJtNUojiJUYw0p1FenUN5GG1opVEeNsAilksI08lpc ERaT1MzWC53oxTlmWuaY8m7Z7zzPqZ3neWyTXrRlbb/+uK/7+vxzX/d9fW+L HXGbJFoaGhoRtN50S7G7q0OfGnlNMw9zIhUahHCzGGpFymLdKuqVxcc0L1BX +EUVbKb2Vcs7TPrVmH4+rqLDXIWEPI9PtAbUmN8dU3XeTAX5sHdMz6Aav5Uu CIyYr0J4QtBXcb+qwdok+Q/MVaG8Jkq7nqXzkwObbxmosLRQFR36VA1x3+Vn 6n9YWHLrvtHS4vDCoTdsvZpFdITmljBDDiHfldeM9rG4t9Pg4lNzDullrl// KGNxqkej+4gdh6mID0JUDSwSfhIXYw2HiT3lT7zKWFQrmXvX/Tjohf27+FEm CxfJngGrUA4OGx7334pnwcevzuiM5qCKCU7v2cZi2tik7u4UDhbKVasWebPY m2n4SnqCQ7xpy4WrK1iM1R3IE5/hYDAInwQLFn3OP76U1XJwM20tTXqbRX/A Tu2DTRwsNyi4ykkGVtnbjzxu5yAZdTJ8rWQwagfFxl4Ojpnyk7VdDBq2xd1R sRxsFl9/tOs6A/VkRlbiBIfNf4SlxVYzcJnP/7Vai8cPxcPVaTkMhsO/7bhp xGO/V72kbT8DlWFw0nIRj9NlybOCwxhM9/Fvf2jPw6Iza3TEi0HB+FLh6Boe 8v6a5svLGUQ17U4b9OXxcIE0UWrKYKvO08Y9oTzSAyzvlWgzaJLtlP8ZxWOG 2PZ6+5gSNrqWcfUpPLqf70v7/oES3F5xjiibR19k3t5vmpWYs8YzNfsMD5+B yZzxaiWCS0uNltXxuPzCwF96Qgmt8PuH2pp4HNsvlEjjlJjy/CzIqJ2Hp6zA pSZYCQPpR90lvTzs58VrX3RRIvu0qWg2yyNfy3L8krkSC16tMGqfoPfetspc oqmESZfzUIymgNyenjnmLIHkTolh1ywB4XmOjYVtBHddPdQ7zQVUza4cqawg 2J3zn0kldab78lIdap2RG36D1MZ/z3oeU07gkWJ31VckALpr013PEsilc1Od LARcaV0/1HOGoMWJN9ZdJMBXoQzhTxE02Ofj0mIBo647qlq/JPA7ujFeRS3P m5H7HvWYQvu8mbWAF8a1j/edIFgae0D/FLX+ttARURZBXX1E3yEbAU9tdHKT MwiqbR2j/ZYIuBmU8Xz6UYLiRUMl4+8LcP/9dERRAkHo5Tl19nYCenz0O9vi CSwRcCOW2tNpvdGzfQTXtnSMjFGnPBk1Cfyc4H5ho4izF3DKjZFM7SUw0Ci6 pFxO9xPnGxtFEvTndrWIPhSw1tqta5mEoMJUryuMutr3+KjnLoJlzml/PaRu SE0aStxB4B0b6zjsICCs7c+uW9sJsn71vn3fSUDR7BGJOpjAP+qLbqMVAn6X Tj1gggjefSkjAdRrV5jkDn9CUGu8SrubOvdJ4gHZJoIfvK03dDoL2PXW5NRn /gQnB3d8qr9SQK2l4rSPH8HHkZVRG9/Y3CvJ1peAOfZObgd1QZcscmgjgcYN 7fttqwQYza+QTPOgeXuC1XIR8LLvuFy+jiC/N+XZWuq7zW9VHBITmI0/ebeV eomnbfIICNysydabrgJmOIR/usGFYFqjacwrah33/du7VxJ0ikNSXdxoPrZq PT9ngi3bFeek1IE6noa2jgQHi28+urZaQPfAwsfr36f/crbxlwlq0TSvgOil BPrVDS0frBHw84Pm9CO29J66C7Wx1GF3ZmYfsCYobzhXWE+tc3dm1WYrgkfX Sg/z1KlBf6jmLSKwaynabQOaV+Gu9isigmRZzqZI6s9JyNwFZgStHRmra6gL FHO3epsS6CmO2DDU6s1F9lbzCDb1JhmL3AXc0Vf2Zs0hKBuOe72NuizRw3uh IQFL9nBnqY83svYv9Qhm+ryeJaMOzL89cUOT4H9/PQ9T "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], AxesOrigin->{0.78, 0.5}, Frame->True, PlotRange-> NCache[{{0.7760254037844386, Rational[1, 2] 3^Rational[1, 2]}, { 0.5002107373333989, 0.6290175841183598}}, {{0.7760254037844386, 0.8660254037844386}, {0.5002107373333989, 0.6290175841183598}}], PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.440930582343*^9, 3.440931811752*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"bb", "=", RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{"qq", "[", "kz", "]"}], "\[LeftDoubleBracket]", "10", "\[RightDoubleBracket]"}], ",", RowBox[{"{", RowBox[{"kz", ",", SqrtBox[ RowBox[{"1", "-", FractionBox[ SuperscriptBox["K", "2"], "4"]}]], ",", RowBox[{ SqrtBox[ RowBox[{"1", "-", FractionBox[ SuperscriptBox["K", "2"], "4"]}]], "+", "0.09"}]}], "}"}], ",", RowBox[{"Frame", "\[Rule]", "True"}], ",", RowBox[{"Axes", "\[Rule]", "False"}]}], "]"}]}]], "Input", CellChangeTimes->{3.440931776776*^9}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV0Xk41HkcB/BZKascpTwdzrLCFo0jR9RbW5FbeJQmpcNNRK6MbbFJ95Ry tPFU24EoREruos2Vxj0z5kdawu83nTrE2m9/fJ7P8/rncy7fF+bmK8VisXxI /Mhb3R8dMdKnUeLe5XB/koK4b6m7M/HLusGWBQwF7HxnfJp448UNHImYwlVB 46Im4qnhUA+zLgosTvaEjAGNb4czOa6tFPYKI7u3EHvXhX+IbaRQz7EvTybW CfNx8KqlsEKkmVlHzDV/c7G2gkLyri+xM8QRvfd28MooXIjQM7RYQ6M7kto7 u5j0S+WMRhAX2nnFlBdQqCyt8Rom1sivKC66QeGzXOq6VjYNaePT6rvTKUhr VXyYY0jDX0HNaX4ahYUWdL41sZoDr7DvHAW2r6tKKXG/3ZvE/BMUAquWTF02 opEu3yJS+Z1CbId9aSdx2eGczl3xFFJGuSEKxjQCVur8nR9L4bryoCiReNh5 hB8YSUEUmlcVYELDJGa5h28gBRf1dUdNTWlEj6TsHfeksNskxCycOCmMPVrh QSHEPudtHrHM0wBfnhuF1GgpHzUzGi6Gubscncm92pqspc1p7G+W65OyobD2 KEeKb0FDyf8b18aUgsoA91joehqaASGDe5ZScDNhqxYQTwzNGUteTOqlDpWM Ee8OvluZr0zuw7Yf8N9AQ+iLBJYShZeJiy33gcaW9O/lTXMpHNcqeee5kfzP vl29bUqMj34jHOstNOKi7LS3vBKjmdlmtNCJRthkyiufYjH+9ZN7cJdD42FU m63JHjGcWNp5rCAa1pf/yTFVFqNV64zVyTga9TzbIrXn/XhQIrglm0rDY8dt tdQ/+1HeHFNcm0Ej6Enlpz7LfsyZHbnK4DaNC9M88CZEqMlz554vo2G6anFe eqkI53qTQ1c3kPk1I+SjDorwPCJrbVUHjaOfXeQ19EVQyfgv+9gQDVX2oVa7 YSEeC+uipD/SqN6qaH/hqhANp4J+ifuJwYLYL9q2e4QYbtgXqb2AwY1bXzeo LxHiULB11jMNBo7xCsO/dQngYMg10DRgsL1jbNPmNAGsEmUOFK1nYBvEzVZx ECDDKjbc1JGBwqS0JVtGANkkc7VPHAbzB94mDFX2oWnE0yIgmEGSimKaUWQf Jh7uXM2PY+AvZXykVr8PZovu1kSeYGCytPBl9mAvhgTDevMyGcykxmhVZPWi vy5hltctBh2BOToRbr14a3QgQVzGYFvIOYG2dC+CFdltXk8ZDFpqH9St6MF/ UQ41cp0M5FqirQcCevDd1vF42SsGLl1jeeMqPbBZd6d60XsGLXvCD4e96Ebe A6/0OywJOH6N49bx3YjOYTabK0rgZreygLuqGzvnv0ybUpdA6Y8v0QUdXTh9 XcktxkACfdakNovbBd2rFiMjVhLYetuo9ep2QaRTpMpzlOB0Wvavy1o7MblM Ts+AI0GmMy15HdUJP4PGG0mBEgiP24WK1Dsx2vyiViFOgt4xgVV2dQe6O6vv 845LUHjTZPXMvg4I/adt7DMkmBTGiLOm+cg0NtN1JC71ab9yhthzJlzWmXhN mCYnkbg9/XXzNmJHvds9AcSNDc2uXsRiRZd2M+Iircs7A4j1M3+u7ZriI4Uy PZhCfFNjRdYC4k13wl1SiR9VBnvOJmbF5LNPEicOFiz89p2PeAX1j2eI2+uU zw4Qh6+fHXuJuECYmHSPmPNXR+IN4nrfJ4FOxEv85fffIu6OOrZyI3G3ke3m XOKzchgyIXZtqphTQJzZmrlLldjm27VT94ntrww7j0/yMatBGFJGfJ46MI8i ruUpO5cTf84VPOMTW+qemP+YuKq1wLqC+OvH+veVxLP0ZKcLiR/UTPGriTXn ej+6Rhx5yrS0lljmUG7UJWL29vBL9T/29x43PEH8P8YLF1o= "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], AxesOrigin->{0.88, 0.3}, Frame->True, PlotRange-> NCache[{{Rational[1, 2] 3^Rational[1, 2], 0.9560254037844386}, { 0.29654149889889886`, 0.49967653661343114`}}, {{0.8660254037844386, 0.9560254037844386}, {0.29654149889889886`, 0.49967653661343114`}}], PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.440930616848*^9, 3.440931843132*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"cc", "=", RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"Im", "[", RowBox[{ RowBox[{"qq", "[", "kz", "]"}], "\[LeftDoubleBracket]", "10", "\[RightDoubleBracket]"}], "]"}]}], "+", "0.5"}], ",", RowBox[{"{", RowBox[{"kz", ",", RowBox[{ SqrtBox[ RowBox[{"1", "-", FractionBox[ SuperscriptBox["K", "2"], "4"]}]], "-", "0.015"}], ",", RowBox[{ SqrtBox[ RowBox[{"1", "-", FractionBox[ SuperscriptBox["K", "2"], "4"]}]], "+", "0.015"}]}], "}"}], ",", RowBox[{"Frame", "\[Rule]", "True"}], ",", RowBox[{"Axes", "\[Rule]", "False"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Dashing", "[", RowBox[{"{", RowBox[{"0.01`", ",", "0.01`"}], "}"}], "]"}], "}"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4409317987349997`*^9, 3.44093180759*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], Dashing[{0.01, 0.01}], LineBox[CompressedData[" 1:eJxN2Hk0ldv/B3BUVKREhVRICIlQhthkiiRUkiFTKRpIRYUrUSpjpkyRhJBZ KYRMIckU5ZYeIuU5pFuiid/7+a71W/f6x3qtczxnn70/+zMQc3K3OMTBxsZ2 hp2NjfmdbjdzK0GDJmz/+6FImWDhijla//pZl2P4Me1/HV/3/GDcf3ywRE2j 5j82znWXfAmbZht9YeeiiEJ6Bt8oXLUuTVBmIUWWJ/39ewYWOZ1fHMpDkd9R fB/5dWjyVbh7SygvRQavGXZJw9saJLV3LqFI00W/Kk3Ys5ZbdZCPIvnnS7It 4MDl8l/38lMkxvNTzGGY522aZ4EARc67rbngC79ZkhD2dRlFHJz2Hr0O29eN N65bQRE5iyc6FXBP+cfwA0IUWWo8JdcOP7qr3XJEmCLTOhsEh2GPZNJ4ZCVF 6hUTxnm30cRr6zsp01UUseM7meQEe/QGh9SJUkRvQdZlb9ihuc0ySowiMuxv T4bCrg4+krbiFPk+sX37AzgtQlNscC1FwttFJxfo0sRof8dwlSRFTjdZUqvh jyVRgcelKGJdE/pMCZ5sWsUrJE0RycLpNDv43dInAYfXU6Qqot20CG4I5fQY lqXIZ9MLd630aHK1/ZusjwJFegweRB+HU/3ppnmKFKnUYv11EfZrcWwNh6/I W+29B09/TVRP3kQRUV4FDjZ9mvB1F//MUKaI2fN3NlnwzVtNV/1UKeJhojxU Dnselimh4chnV461wbmTZ6at1CjS3qzoNwn/da1jVEEdf98YmKJrQBMzv76J Hg24WmrgHdzWXyUwq0WRXUXuLoLbaTKiOc+wR48iJxTqx2Vh4+JRgY362J8C QW8Cl66pdg6G2/KeBLvAa43zQzYbUMQ0hz+7FCY1RzdHGMLpZbSZEU1UXpk4 yxrDsbMeV41p8mWDVZngLpzPiOaNZPhouDOfM6yu5vu4AD7XFCGcB697O72g B+aueKqpbYZ4lPjntvgOmvAqaL5zMqdIdulQdwUsK7ZZK3E3ReZ2N6uOm9Dk kP083Z59FKlYGsO2exdN4m8k8922x/m0WTQehuO+z1UZhGuv8YX6wr/FrcTF HRDPHBErMuFDDm62aXD3P1flf8CjEyn7bzpSZLzT3zbVjCbqbwbaQpxxPtFH H9LmNKFUDgnrH6bIJX5dj6A9NDHxT72ndgLPG4yWstuP1+f1f288R5EjKrNV Nw/Q5LBzXMfdIIqc1XzX3uCM+KvmW74onCJsAo4ysUdoUvT5s0l/LEWGj/9x bj1Ok1mbsV3nblJEN05jOMeTJs0FZt5v0imy5PQ45eGN/ZD54BCcTRGNlxWV Kb402WpjmMhfQJH+ubLUqQCafKpyUKksgctK3xdexnm2WSzXfYjzEC1xfBmC eGVPG1eppMjmfbul0iJpkm1lFXOumiI+fLzPu2Npclak3i+mliJqiW7lBYl4 ntrplmUNFHHzutVQDy/dVtK9Bk606Gl/DZu85X61Hv65cNvInCTs96pFLZpw pY+QgBU8+NDF5hCsbdd0nC0Z+WWRrHIJbCAqKWZ+kybLkt3rjBopsidzIPif VJqc2t0Z4PQU928vuyLXLZpo16YUucEOc8X6VsJqHF/6POHjTg5y+nCB0Xrh QPjy6nftN+Dy1T+M0+GyuDdCmmmI35YYr0FYOLg3N/g2TbRcw37YNOE8N0/t SYal7KVsnGHJ4eUzhbDGUasHbvAmXUuzPrirPmD/OXjHbPc3uXSaVIrPWsTB vt6dmp3w19jF/C/gd4efvxDJoMnphooWjWY8P4V7uQKcMqr5WQd26Day1YUX harybocp7caPrjAdG6+xh7FQNXsZPDLFs+cY42eFSuaZNNFt4PVNYrwx5sal LJzPgqScSeb5Lp1v42HngaerfzHPT14icQ8+NuAVxtaC988PK+yEz7jbWnMz HrjUvOYuTTao7SwWhQeivX89gkOTZq2N4cEpa/vxbJpMmNhfT4AV4oX05+TQ JFalvzIF9lN9JSMIx2V+GUiHBc/t/a4DJ2/pFcyHTX/uCouFrT3uqNXC5b/1 Krbm4n7xXTD8BEexb1wRcg/rWTwQu+kZ1pM29jsVbpl1/bYZ3rDt3mApHPgy 3ETjGVNv1uf3w8Nf7n3ShWfmSOhvyqPJt6XavXtgN06hU6/hrfM++Z+Bt3HP eSFVgHjdaze/GG7ZRJlrwCb15MB92Nz6cbcprCYXk/8Qdsj27jsDJ+lFa1XD foas4Tr4MefYylb4YeDL3/aFiKdC0dYhWDu32P8UzBMn1zwCP+2M4AiGWQW5 NaNwj5jx/Hz4cN2c1An4a3UV/2+Yl22V0G9Y/neWTHwR1rNplmNpK0UerL2U lwvvfK6XJQBr7nBSqIbNnWb0V8AmiSKbP8AuSc89RWBX1es6ysWIJ8G3VyXh O6fPW7XD/aINeWrwyjGTy/NLaWJ7sqnVBj47wp+1EhYu7blvB/cM9D2Vh6d2 RSbYw5E9Rxbshfl85pk7w3NrgkLT4Pn93/3d4PGoyij1+zTpcxh/6w3vCAss 2Qnr/5V/8hx8N9i42wHWcRJi84GdfF8tC4ZP7PjM8xfce+hbfBe8QiLxQyBc oyqXeuwBTQxC+xPC4VVKX6v9YW9FvVcR8PkN5VQULJS6m+86rCy+XfwRzAp0 PxUNZ3MfzJhXRpNNaglTN2BOTtlGQXh22Jg7AT7I9s8HWXhlSYZQIrx60l/a AnY0t5RMhqP6k3JT4Gdz02ZT4c+vnFqLYE3hhI+3mP3uWj9WDw+FqrWmMc9v KpOn4Z41pwLTYZ+irqItD2kyw95YlAm/yk3sNIbPzJ9xz4I3Zzp+tYN1wt9K 34UnEj8rB8EaKVHXspn1BXE/bIdlfDoW3YPbPlqcHYa92i/lMVbbmaj6Ey5w HjfIg3mXr3+09hFNLiT+7ZIPnzvvcU4VjhmxH2M81F+mthPub7hxvAB+lGVQ 7gXvCBSzKYQleMLPh8AfOU81MQ73eKl+C750yHtjEfyjW+RXKWzgIR/J+KDa wYpmWHnqxijjtpu5Pv1wVEuJVjGzPo6vGl/huY8CQhnfcVH/zVVOkx+R7N2M eZ8FVIrAe8SUlpUw8bex2VcR/rl5uRnj99FLNA3gwyk5QYxNp/f9sYZrJCaK GT+0TX3sDn+I+dDHWPzJB78geH552G/GYevktRLgaweGV5TC01fPzOTBN03G 5Bg7jVdW1cKvTTI0GLdazPXvhTnWL9JjvKVsB2HBC+5KGDC+vTJ6lq0Cr0fR Oox5LvRVL4Mn8hy3MPYeErsgA9s9uSzJeHC7qzaBhaOtFzM2yStk2wPzj7z5 wqyvjG+65gh8/+Kitv+t34sE+MHH1rPSGYf2XdaJYj4v9ewpxlNabexZ8NCd 4q2MndKX1VbA9EzCLLOfrVx2F9vhmuPylYy3HLuzbRjWLXb3ZMyjolTHW0kT icdTz5jz8ko4H7gWvs2h5c544M8TXVXYaqUCD+MHjWb1jrBi9XIlJh5EZeOD vOADdQsfM/ETEvFOLwTeHZGtzdjB6kRDKZx29IMyE28tlfcvNcPqXwIzmPhU Efuj3w8vXP9sMeOFoyGNXI9pcjJqb0cufMa067IInKAcL82YKhY2VITHXQPP 5cD3fbKfWsMthy3mM/dDlJoIdofdk1QMmPtzTU91exA8pVbjx9wvh0VPm/Jg u0strzOYz0sdamaroknnpxz923D+NWMpAThH0cyRub+7vQqDpGC5+GPezP1O 3ulLTOHCxwuiU1qZfMx/Pwm2F20LZvJF58i5pQXwxJp13kw+8e56514Lb5P8 4hDP5K+cXJlPMPemkbVxsMV+3Vubq7FfSbqHmXz1XS/7jxF8SUpDJhJOUlhs Ywe/D8oYYfLdENffy4Lg6rAO81DmvB54hnTAd+r+en8ZTuC/7XWsBusNpX+d h7Vmubr9GU+zDjH5dnD0uGI0HHrCtoXJxzK1aqxHcJ160MXTcKVHhyPXE+SX g+MFx5n9fs628zZ8QnZo5gAc9Mgl5z6cGJa6mKkH6zNauZph3Yo6YaZeePrG 136GE4T0Rfcx9UBWQVWzliZv14p/NIWlr9ivfQX3v59004T73uqLjsESQ60T 6kx8K8mt4qhD/Zrk9lBl8l//9HI52Na5xlqJud8qUQsvwAvpb5Q0bPi+/otU PU2uO68uZOrhtGru+FZYokOjdTGcE36dNoeT7B8N8DD5Rf3AsA8s3dw5zcnU s8ipVy/gW6s9en+ivrpoydR4N6A/mHqzdwAWjFnyOBS2z6jk6mf6gU/fH6XB vF1fi/uYehxbV/IM5jzhMtkFT9K2WWsaaSL50ESmEQ5KiAh/Cvv/SNqdDad/ +2a7ogn7V7z0lhsss6phShaOy5WOcIGL9WOjtOFiySfeTnBtnEqTK1wu661i zfQ7ql5KlTDb7LCTEbzab2qBYzPqObfOJUk4ft7PBzkt6GccrN//jX5JVL7F vBo2tA5g74XvWiayumCn2mVCnXBZlrr4H7h9YZByE9xj5BO26xlNjiQ6iZTC /OG/nb/BOTKXC67B4ctnl2g9p0nIWeX+TS3MfDi23RzOL5EO2QDzeb25cBD+ ZWO2SRqOeflo4hosbbfFfRWcGHO6vRfu/GiezgVnLR2NONmG99ut+9aH/rOG 9yVvxgvMX11XlPzgfzhzeXg60T+7OrIXoj9+Iio8JAq/Kb5y+y4cqX61QgVe 3J61NQ2WP3HYzR6mPJttomC3l2ubiuHSKD/90/BgenLg/i6aDCyvadgMd2qH /8zsRj92UdnzAfr7kvMnP+r0on5IjWncwjxRuXLWdi/82zIz/wZcXxnacQTu jpwVjoB7ZrLKI2CeheGUH/zzYn/oW5hPZGa1DbwtZIfiuVf4/OPaRwTg9qR1 5wtfo/+Z2y5xEfMLq/I195o3+D42zlzG9RSJLX2Q/Z5CP+efrRn/BPd/Y77U kg+oV3OTPitWYf44Zzj2fhTzh25+ZWg58oX+r6g1E5iPy+cFcz+giPGA2QWO SZo0bbJcZFiEeN+08mPrD+QH93FfrnsU2cjJpjA7QxM/y2q5WxkUab57ufrV HBa5dvJd5IcU7H+b1EbhBSyyI2c7yYvDPJOjKaTLyyIPDY4Y+4VRZN6eSWN+ fhbh47TiKQukiE7AU9sdgizyKumR+Blv9PvmG/oOwKvaK07pwqvFovd5wgbB VXf44N81BywS4JAc8cw8L9xntu+GI3Ai21uxoTM4vwsSSoFCLOJUH9prfJoi Iv4BCyqEWWRN05YffzzwfX01HsisYpGZz5lWn10xH6wJ6dCEt7+9b3Qf9q/t Y5nBf5/24fWB588/v9Yb3szNsZQTFo5+GFkH62QUTAsdoQjJVjlqu5pFrKcF Z1VdKHLl5UbR8DUs4mZzINzcCfX5rL9GGly+Zt4ePphv5QvLUnh03rrBdsz3 iQ4nwvrg3lIBy13wPfreT0lRPE9+a6SRA0U62GVeVsNKBy4ObTqAz5eXuPZF jEVMrvQ4vNyP81f/pscpziKe6/iEguEqw/pZYViwwNRHFXZyPHhGDz7lqKGe aIV8FnPH4Qb8U745eP8+iqj/ktiiuZZF2mJ6Ahr3YN5tWjd0RYJF7kxwvDDf RRG5l5MpN+HtZa13v5ti/hlo2F8MWzpKz0uCb/869OJv2IEjQ2twJ0VG5TPL 5dexSGZYZYCbCfrpWMnr3fDZpstXHIwQj05SWqJSLFLV9vpoxzaKZH5fFqcE Rz4NyHKFi6/NHTeAKwPOyrLDLcUDN4/Dzyt8Zjbo4H5wJP8phxfx+Lj+RbDe 9CWP90qziLNxeeZ3DeTfLbMCrrBHZsLFANjz2dgxX5hraKE9N3zta4tIOmzz R51dRB39qu4l3wk4q22PooIq8vXQD42Q9SxSULrCdZUy+tezH6NT4KWqO1/E KCE+eXrpIvhXAMlbCG9XLk16BQ/m77v5RRH1OujEr3UyLGJq/b0qZyNFLgra 7VaDO0pT9gnA4fd25JrAt86Ys3zk8f27pW1OwWwnVz/U24B8IPG+vAbuTIs0 KpdBvnnYsbQbHvgzoMYHj5vUuI3AL9JfWR5aj/ngzE1hXlkWcWk4HTRXGv1/ w77zNrDdTv9RoXXIH/sNOk/A5Y8lr+2VQL87pixzEZ6v190Yuhb7tWxp3114 1oy7nyWG/HGoVe07vDloJPrIavQzzhqK3HIssvG4WJzjKszDjjnSorC1dkqP mQjOw+7KCmP4mBS31iJhigjYTvHaw/TrYKNeQczT1i6cp+GvZx+uv74C/ZOl 3uRNuNR6g2OTAOaJPSWsYnhcIHJqKz/mIwvxoaewrnVhSDgf7oPZ9b/fwPNO GAwULaZIgylb1xd41WPD2ehFFOk3cW/h3MAixjvOlSzmRv027n+yEt46HNxU x4X7ZrTzkQL8//+PljGsLNT/jye3KB8b+49rpO9Jxcr/a8uFySlDG/81iw5Z pqz4r9v4U8W3/8dFGsUbbf/j/wMnUB9X "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], AxesOrigin->{0.855, 0.5}, Frame->True, PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.440931032328*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"aa", ",", "bb", ",", "cc", ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{ SubscriptBox["k", "z"], ",", "q"}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[{{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV0n04lXcYB3AvRayJtNUojiJUYw0p1FenUN5GG1opVEeNsAilksI08lpc ERaT1MzWC53oxTlmWuaY8m7Z7zzPqZ3neWyTXrRlbb/+uK/7+vxzX/d9fW+L HXGbJFoaGhoRtN50S7G7q0OfGnlNMw9zIhUahHCzGGpFymLdKuqVxcc0L1BX +EUVbKb2Vcs7TPrVmH4+rqLDXIWEPI9PtAbUmN8dU3XeTAX5sHdMz6Aav5Uu CIyYr0J4QtBXcb+qwdok+Q/MVaG8Jkq7nqXzkwObbxmosLRQFR36VA1x3+Vn 6n9YWHLrvtHS4vDCoTdsvZpFdITmljBDDiHfldeM9rG4t9Pg4lNzDullrl// KGNxqkej+4gdh6mID0JUDSwSfhIXYw2HiT3lT7zKWFQrmXvX/Tjohf27+FEm CxfJngGrUA4OGx7334pnwcevzuiM5qCKCU7v2cZi2tik7u4UDhbKVasWebPY m2n4SnqCQ7xpy4WrK1iM1R3IE5/hYDAInwQLFn3OP76U1XJwM20tTXqbRX/A Tu2DTRwsNyi4ykkGVtnbjzxu5yAZdTJ8rWQwagfFxl4Ojpnyk7VdDBq2xd1R sRxsFl9/tOs6A/VkRlbiBIfNf4SlxVYzcJnP/7Vai8cPxcPVaTkMhsO/7bhp xGO/V72kbT8DlWFw0nIRj9NlybOCwxhM9/Fvf2jPw6Iza3TEi0HB+FLh6Boe 8v6a5svLGUQ17U4b9OXxcIE0UWrKYKvO08Y9oTzSAyzvlWgzaJLtlP8ZxWOG 2PZ6+5gSNrqWcfUpPLqf70v7/oES3F5xjiibR19k3t5vmpWYs8YzNfsMD5+B yZzxaiWCS0uNltXxuPzCwF96Qgmt8PuH2pp4HNsvlEjjlJjy/CzIqJ2Hp6zA pSZYCQPpR90lvTzs58VrX3RRIvu0qWg2yyNfy3L8krkSC16tMGqfoPfetspc oqmESZfzUIymgNyenjnmLIHkTolh1ywB4XmOjYVtBHddPdQ7zQVUza4cqawg 2J3zn0kldab78lIdap2RG36D1MZ/z3oeU07gkWJ31VckALpr013PEsilc1Od LARcaV0/1HOGoMWJN9ZdJMBXoQzhTxE02Ofj0mIBo647qlq/JPA7ujFeRS3P m5H7HvWYQvu8mbWAF8a1j/edIFgae0D/FLX+ttARURZBXX1E3yEbAU9tdHKT MwiqbR2j/ZYIuBmU8Xz6UYLiRUMl4+8LcP/9dERRAkHo5Tl19nYCenz0O9vi CSwRcCOW2tNpvdGzfQTXtnSMjFGnPBk1Cfyc4H5ho4izF3DKjZFM7SUw0Ci6 pFxO9xPnGxtFEvTndrWIPhSw1tqta5mEoMJUryuMutr3+KjnLoJlzml/PaRu SE0aStxB4B0b6zjsICCs7c+uW9sJsn71vn3fSUDR7BGJOpjAP+qLbqMVAn6X Tj1gggjefSkjAdRrV5jkDn9CUGu8SrubOvdJ4gHZJoIfvK03dDoL2PXW5NRn /gQnB3d8qr9SQK2l4rSPH8HHkZVRG9/Y3CvJ1peAOfZObgd1QZcscmgjgcYN 7fttqwQYza+QTPOgeXuC1XIR8LLvuFy+jiC/N+XZWuq7zW9VHBITmI0/ebeV eomnbfIICNysydabrgJmOIR/usGFYFqjacwrah33/du7VxJ0ikNSXdxoPrZq PT9ngi3bFeek1IE6noa2jgQHi28+urZaQPfAwsfr36f/crbxlwlq0TSvgOil BPrVDS0frBHw84Pm9CO29J66C7Wx1GF3ZmYfsCYobzhXWE+tc3dm1WYrgkfX Sg/z1KlBf6jmLSKwaynabQOaV+Gu9isigmRZzqZI6s9JyNwFZgStHRmra6gL FHO3epsS6CmO2DDU6s1F9lbzCDb1JhmL3AXc0Vf2Zs0hKBuOe72NuizRw3uh IQFL9nBnqY83svYv9Qhm+ryeJaMOzL89cUOT4H9/PQ9T "]]}}, {{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV0Xk41HkcB/BZKascpTwdzrLCFo0jR9RbW5FbeJQmpcNNRK6MbbFJ95Ry tPFU24EoREruos2Vxj0z5kdawu83nTrE2m9/fJ7P8/rncy7fF+bmK8VisXxI /Mhb3R8dMdKnUeLe5XB/koK4b6m7M/HLusGWBQwF7HxnfJp448UNHImYwlVB 46Im4qnhUA+zLgosTvaEjAGNb4czOa6tFPYKI7u3EHvXhX+IbaRQz7EvTybW CfNx8KqlsEKkmVlHzDV/c7G2gkLyri+xM8QRvfd28MooXIjQM7RYQ6M7kto7 u5j0S+WMRhAX2nnFlBdQqCyt8Rom1sivKC66QeGzXOq6VjYNaePT6rvTKUhr VXyYY0jDX0HNaX4ahYUWdL41sZoDr7DvHAW2r6tKKXG/3ZvE/BMUAquWTF02 opEu3yJS+Z1CbId9aSdx2eGczl3xFFJGuSEKxjQCVur8nR9L4bryoCiReNh5 hB8YSUEUmlcVYELDJGa5h28gBRf1dUdNTWlEj6TsHfeksNskxCycOCmMPVrh QSHEPudtHrHM0wBfnhuF1GgpHzUzGi6Gubscncm92pqspc1p7G+W65OyobD2 KEeKb0FDyf8b18aUgsoA91joehqaASGDe5ZScDNhqxYQTwzNGUteTOqlDpWM Ee8OvluZr0zuw7Yf8N9AQ+iLBJYShZeJiy33gcaW9O/lTXMpHNcqeee5kfzP vl29bUqMj34jHOstNOKi7LS3vBKjmdlmtNCJRthkyiufYjH+9ZN7cJdD42FU m63JHjGcWNp5rCAa1pf/yTFVFqNV64zVyTga9TzbIrXn/XhQIrglm0rDY8dt tdQ/+1HeHFNcm0Ej6Enlpz7LfsyZHbnK4DaNC9M88CZEqMlz554vo2G6anFe eqkI53qTQ1c3kPk1I+SjDorwPCJrbVUHjaOfXeQ19EVQyfgv+9gQDVX2oVa7 YSEeC+uipD/SqN6qaH/hqhANp4J+ifuJwYLYL9q2e4QYbtgXqb2AwY1bXzeo LxHiULB11jMNBo7xCsO/dQngYMg10DRgsL1jbNPmNAGsEmUOFK1nYBvEzVZx ECDDKjbc1JGBwqS0JVtGANkkc7VPHAbzB94mDFX2oWnE0yIgmEGSimKaUWQf Jh7uXM2PY+AvZXykVr8PZovu1kSeYGCytPBl9mAvhgTDevMyGcykxmhVZPWi vy5hltctBh2BOToRbr14a3QgQVzGYFvIOYG2dC+CFdltXk8ZDFpqH9St6MF/ UQ41cp0M5FqirQcCevDd1vF42SsGLl1jeeMqPbBZd6d60XsGLXvCD4e96Ebe A6/0OywJOH6N49bx3YjOYTabK0rgZreygLuqGzvnv0ybUpdA6Y8v0QUdXTh9 XcktxkACfdakNovbBd2rFiMjVhLYetuo9ep2QaRTpMpzlOB0Wvavy1o7MblM Ts+AI0GmMy15HdUJP4PGG0mBEgiP24WK1Dsx2vyiViFOgt4xgVV2dQe6O6vv 845LUHjTZPXMvg4I/adt7DMkmBTGiLOm+cg0NtN1JC71ab9yhthzJlzWmXhN mCYnkbg9/XXzNmJHvds9AcSNDc2uXsRiRZd2M+Iircs7A4j1M3+u7ZriI4Uy PZhCfFNjRdYC4k13wl1SiR9VBnvOJmbF5LNPEicOFiz89p2PeAX1j2eI2+uU zw4Qh6+fHXuJuECYmHSPmPNXR+IN4nrfJ4FOxEv85fffIu6OOrZyI3G3ke3m XOKzchgyIXZtqphTQJzZmrlLldjm27VT94ntrww7j0/yMatBGFJGfJ46MI8i ruUpO5cTf84VPOMTW+qemP+YuKq1wLqC+OvH+veVxLP0ZKcLiR/UTPGriTXn ej+6Rhx5yrS0lljmUG7UJWL29vBL9T/29x43PEH8P8YLF1o= "]]}}, {{}, {}, {Hue[0.67, 0.6, 0.6], Dashing[{0.01, 0.01}], LineBox[CompressedData[" 1:eJxN2Hk0ldv/B3BUVKREhVRICIlQhthkiiRUkiFTKRpIRYUrUSpjpkyRhJBZ KYRMIckU5ZYeIuU5pFuiid/7+a71W/f6x3qtczxnn70/+zMQc3K3OMTBxsZ2 hp2NjfmdbjdzK0GDJmz/+6FImWDhijla//pZl2P4Me1/HV/3/GDcf3ywRE2j 5j82znWXfAmbZht9YeeiiEJ6Bt8oXLUuTVBmIUWWJ/39ewYWOZ1fHMpDkd9R fB/5dWjyVbh7SygvRQavGXZJw9saJLV3LqFI00W/Kk3Ys5ZbdZCPIvnnS7It 4MDl8l/38lMkxvNTzGGY522aZ4EARc67rbngC79ZkhD2dRlFHJz2Hr0O29eN N65bQRE5iyc6FXBP+cfwA0IUWWo8JdcOP7qr3XJEmCLTOhsEh2GPZNJ4ZCVF 6hUTxnm30cRr6zsp01UUseM7meQEe/QGh9SJUkRvQdZlb9ihuc0ySowiMuxv T4bCrg4+krbiFPk+sX37AzgtQlNscC1FwttFJxfo0sRof8dwlSRFTjdZUqvh jyVRgcelKGJdE/pMCZ5sWsUrJE0RycLpNDv43dInAYfXU6Qqot20CG4I5fQY lqXIZ9MLd630aHK1/ZusjwJFegweRB+HU/3ppnmKFKnUYv11EfZrcWwNh6/I W+29B09/TVRP3kQRUV4FDjZ9mvB1F//MUKaI2fN3NlnwzVtNV/1UKeJhojxU Dnselimh4chnV461wbmTZ6at1CjS3qzoNwn/da1jVEEdf98YmKJrQBMzv76J Hg24WmrgHdzWXyUwq0WRXUXuLoLbaTKiOc+wR48iJxTqx2Vh4+JRgY362J8C QW8Cl66pdg6G2/KeBLvAa43zQzYbUMQ0hz+7FCY1RzdHGMLpZbSZEU1UXpk4 yxrDsbMeV41p8mWDVZngLpzPiOaNZPhouDOfM6yu5vu4AD7XFCGcB697O72g B+aueKqpbYZ4lPjntvgOmvAqaL5zMqdIdulQdwUsK7ZZK3E3ReZ2N6uOm9Dk kP083Z59FKlYGsO2exdN4m8k8922x/m0WTQehuO+z1UZhGuv8YX6wr/FrcTF HRDPHBErMuFDDm62aXD3P1flf8CjEyn7bzpSZLzT3zbVjCbqbwbaQpxxPtFH H9LmNKFUDgnrH6bIJX5dj6A9NDHxT72ndgLPG4yWstuP1+f1f288R5EjKrNV Nw/Q5LBzXMfdIIqc1XzX3uCM+KvmW74onCJsAo4ysUdoUvT5s0l/LEWGj/9x bj1Ok1mbsV3nblJEN05jOMeTJs0FZt5v0imy5PQ45eGN/ZD54BCcTRGNlxWV Kb402WpjmMhfQJH+ubLUqQCafKpyUKksgctK3xdexnm2WSzXfYjzEC1xfBmC eGVPG1eppMjmfbul0iJpkm1lFXOumiI+fLzPu2Npclak3i+mliJqiW7lBYl4 ntrplmUNFHHzutVQDy/dVtK9Bk606Gl/DZu85X61Hv65cNvInCTs96pFLZpw pY+QgBU8+NDF5hCsbdd0nC0Z+WWRrHIJbCAqKWZ+kybLkt3rjBopsidzIPif VJqc2t0Z4PQU928vuyLXLZpo16YUucEOc8X6VsJqHF/6POHjTg5y+nCB0Xrh QPjy6nftN+Dy1T+M0+GyuDdCmmmI35YYr0FYOLg3N/g2TbRcw37YNOE8N0/t SYal7KVsnGHJ4eUzhbDGUasHbvAmXUuzPrirPmD/OXjHbPc3uXSaVIrPWsTB vt6dmp3w19jF/C/gd4efvxDJoMnphooWjWY8P4V7uQKcMqr5WQd26Day1YUX harybocp7caPrjAdG6+xh7FQNXsZPDLFs+cY42eFSuaZNNFt4PVNYrwx5sal LJzPgqScSeb5Lp1v42HngaerfzHPT14icQ8+NuAVxtaC988PK+yEz7jbWnMz HrjUvOYuTTao7SwWhQeivX89gkOTZq2N4cEpa/vxbJpMmNhfT4AV4oX05+TQ JFalvzIF9lN9JSMIx2V+GUiHBc/t/a4DJ2/pFcyHTX/uCouFrT3uqNXC5b/1 Krbm4n7xXTD8BEexb1wRcg/rWTwQu+kZ1pM29jsVbpl1/bYZ3rDt3mApHPgy 3ETjGVNv1uf3w8Nf7n3ShWfmSOhvyqPJt6XavXtgN06hU6/hrfM++Z+Bt3HP eSFVgHjdaze/GG7ZRJlrwCb15MB92Nz6cbcprCYXk/8Qdsj27jsDJ+lFa1XD foas4Tr4MefYylb4YeDL3/aFiKdC0dYhWDu32P8UzBMn1zwCP+2M4AiGWQW5 NaNwj5jx/Hz4cN2c1An4a3UV/2+Yl22V0G9Y/neWTHwR1rNplmNpK0UerL2U lwvvfK6XJQBr7nBSqIbNnWb0V8AmiSKbP8AuSc89RWBX1es6ysWIJ8G3VyXh O6fPW7XD/aINeWrwyjGTy/NLaWJ7sqnVBj47wp+1EhYu7blvB/cM9D2Vh6d2 RSbYw5E9Rxbshfl85pk7w3NrgkLT4Pn93/3d4PGoyij1+zTpcxh/6w3vCAss 2Qnr/5V/8hx8N9i42wHWcRJi84GdfF8tC4ZP7PjM8xfce+hbfBe8QiLxQyBc oyqXeuwBTQxC+xPC4VVKX6v9YW9FvVcR8PkN5VQULJS6m+86rCy+XfwRzAp0 PxUNZ3MfzJhXRpNNaglTN2BOTtlGQXh22Jg7AT7I9s8HWXhlSYZQIrx60l/a AnY0t5RMhqP6k3JT4Gdz02ZT4c+vnFqLYE3hhI+3mP3uWj9WDw+FqrWmMc9v KpOn4Z41pwLTYZ+irqItD2kyw95YlAm/yk3sNIbPzJ9xz4I3Zzp+tYN1wt9K 34UnEj8rB8EaKVHXspn1BXE/bIdlfDoW3YPbPlqcHYa92i/lMVbbmaj6Ey5w HjfIg3mXr3+09hFNLiT+7ZIPnzvvcU4VjhmxH2M81F+mthPub7hxvAB+lGVQ 7gXvCBSzKYQleMLPh8AfOU81MQ73eKl+C750yHtjEfyjW+RXKWzgIR/J+KDa wYpmWHnqxijjtpu5Pv1wVEuJVjGzPo6vGl/huY8CQhnfcVH/zVVOkx+R7N2M eZ8FVIrAe8SUlpUw8bex2VcR/rl5uRnj99FLNA3gwyk5QYxNp/f9sYZrJCaK GT+0TX3sDn+I+dDHWPzJB78geH552G/GYevktRLgaweGV5TC01fPzOTBN03G 5Bg7jVdW1cKvTTI0GLdazPXvhTnWL9JjvKVsB2HBC+5KGDC+vTJ6lq0Cr0fR Oox5LvRVL4Mn8hy3MPYeErsgA9s9uSzJeHC7qzaBhaOtFzM2yStk2wPzj7z5 wqyvjG+65gh8/+Kitv+t34sE+MHH1rPSGYf2XdaJYj4v9ewpxlNabexZ8NCd 4q2MndKX1VbA9EzCLLOfrVx2F9vhmuPylYy3HLuzbRjWLXb3ZMyjolTHW0kT icdTz5jz8ko4H7gWvs2h5c544M8TXVXYaqUCD+MHjWb1jrBi9XIlJh5EZeOD vOADdQsfM/ETEvFOLwTeHZGtzdjB6kRDKZx29IMyE28tlfcvNcPqXwIzmPhU Efuj3w8vXP9sMeOFoyGNXI9pcjJqb0cufMa067IInKAcL82YKhY2VITHXQPP 5cD3fbKfWsMthy3mM/dDlJoIdofdk1QMmPtzTU91exA8pVbjx9wvh0VPm/Jg u0strzOYz0sdamaroknnpxz923D+NWMpAThH0cyRub+7vQqDpGC5+GPezP1O 3ulLTOHCxwuiU1qZfMx/Pwm2F20LZvJF58i5pQXwxJp13kw+8e56514Lb5P8 4hDP5K+cXJlPMPemkbVxsMV+3Vubq7FfSbqHmXz1XS/7jxF8SUpDJhJOUlhs Ywe/D8oYYfLdENffy4Lg6rAO81DmvB54hnTAd+r+en8ZTuC/7XWsBusNpX+d h7Vmubr9GU+zDjH5dnD0uGI0HHrCtoXJxzK1aqxHcJ160MXTcKVHhyPXE+SX g+MFx5n9fs628zZ8QnZo5gAc9Mgl5z6cGJa6mKkH6zNauZph3Yo6YaZeePrG 136GE4T0Rfcx9UBWQVWzliZv14p/NIWlr9ivfQX3v59004T73uqLjsESQ60T 6kx8K8mt4qhD/Zrk9lBl8l//9HI52Na5xlqJud8qUQsvwAvpb5Q0bPi+/otU PU2uO68uZOrhtGru+FZYokOjdTGcE36dNoeT7B8N8DD5Rf3AsA8s3dw5zcnU s8ipVy/gW6s9en+ivrpoydR4N6A/mHqzdwAWjFnyOBS2z6jk6mf6gU/fH6XB vF1fi/uYehxbV/IM5jzhMtkFT9K2WWsaaSL50ESmEQ5KiAh/Cvv/SNqdDad/ +2a7ogn7V7z0lhsss6phShaOy5WOcIGL9WOjtOFiySfeTnBtnEqTK1wu661i zfQ7ql5KlTDb7LCTEbzab2qBYzPqObfOJUk4ft7PBzkt6GccrN//jX5JVL7F vBo2tA5g74XvWiayumCn2mVCnXBZlrr4H7h9YZByE9xj5BO26xlNjiQ6iZTC /OG/nb/BOTKXC67B4ctnl2g9p0nIWeX+TS3MfDi23RzOL5EO2QDzeb25cBD+ ZWO2SRqOeflo4hosbbfFfRWcGHO6vRfu/GiezgVnLR2NONmG99ut+9aH/rOG 9yVvxgvMX11XlPzgfzhzeXg60T+7OrIXoj9+Iio8JAq/Kb5y+y4cqX61QgVe 3J61NQ2WP3HYzR6mPJttomC3l2ubiuHSKD/90/BgenLg/i6aDCyvadgMd2qH /8zsRj92UdnzAfr7kvMnP+r0on5IjWncwjxRuXLWdi/82zIz/wZcXxnacQTu jpwVjoB7ZrLKI2CeheGUH/zzYn/oW5hPZGa1DbwtZIfiuVf4/OPaRwTg9qR1 5wtfo/+Z2y5xEfMLq/I195o3+D42zlzG9RSJLX2Q/Z5CP+efrRn/BPd/Y77U kg+oV3OTPitWYf44Zzj2fhTzh25+ZWg58oX+r6g1E5iPy+cFcz+giPGA2QWO SZo0bbJcZFiEeN+08mPrD+QH93FfrnsU2cjJpjA7QxM/y2q5WxkUab57ufrV HBa5dvJd5IcU7H+b1EbhBSyyI2c7yYvDPJOjKaTLyyIPDY4Y+4VRZN6eSWN+ fhbh47TiKQukiE7AU9sdgizyKumR+Blv9PvmG/oOwKvaK07pwqvFovd5wgbB VXf44N81BywS4JAc8cw8L9xntu+GI3Ai21uxoTM4vwsSSoFCLOJUH9prfJoi Iv4BCyqEWWRN05YffzzwfX01HsisYpGZz5lWn10xH6wJ6dCEt7+9b3Qf9q/t Y5nBf5/24fWB588/v9Yb3szNsZQTFo5+GFkH62QUTAsdoQjJVjlqu5pFrKcF Z1VdKHLl5UbR8DUs4mZzINzcCfX5rL9GGly+Zt4ePphv5QvLUnh03rrBdsz3 iQ4nwvrg3lIBy13wPfreT0lRPE9+a6SRA0U62GVeVsNKBy4ObTqAz5eXuPZF jEVMrvQ4vNyP81f/pscpziKe6/iEguEqw/pZYViwwNRHFXZyPHhGDz7lqKGe aIV8FnPH4Qb8U745eP8+iqj/ktiiuZZF2mJ6Ahr3YN5tWjd0RYJF7kxwvDDf RRG5l5MpN+HtZa13v5ti/hlo2F8MWzpKz0uCb/869OJv2IEjQ2twJ0VG5TPL 5dexSGZYZYCbCfrpWMnr3fDZpstXHIwQj05SWqJSLFLV9vpoxzaKZH5fFqcE Rz4NyHKFi6/NHTeAKwPOyrLDLcUDN4/Dzyt8Zjbo4H5wJP8phxfx+Lj+RbDe 9CWP90qziLNxeeZ3DeTfLbMCrrBHZsLFANjz2dgxX5hraKE9N3zta4tIOmzz R51dRB39qu4l3wk4q22PooIq8vXQD42Q9SxSULrCdZUy+tezH6NT4KWqO1/E KCE+eXrpIvhXAMlbCG9XLk16BQ/m77v5RRH1OujEr3UyLGJq/b0qZyNFLgra 7VaDO0pT9gnA4fd25JrAt86Ys3zk8f27pW1OwWwnVz/U24B8IPG+vAbuTIs0 KpdBvnnYsbQbHvgzoMYHj5vUuI3AL9JfWR5aj/ngzE1hXlkWcWk4HTRXGv1/ w77zNrDdTv9RoXXIH/sNOk/A5Y8lr+2VQL87pixzEZ6v190Yuhb7tWxp3114 1oy7nyWG/HGoVe07vDloJPrIavQzzhqK3HIssvG4WJzjKszDjjnSorC1dkqP mQjOw+7KCmP4mBS31iJhigjYTvHaw/TrYKNeQczT1i6cp+GvZx+uv74C/ZOl 3uRNuNR6g2OTAOaJPSWsYnhcIHJqKz/mIwvxoaewrnVhSDgf7oPZ9b/fwPNO GAwULaZIgylb1xd41WPD2ehFFOk3cW/h3MAixjvOlSzmRv027n+yEt46HNxU x4X7ZrTzkQL8//+PljGsLNT/jye3KB8b+49rpO9Jxcr/a8uFySlDG/81iw5Z pqz4r9v4U8W3/8dFGsUbbf/j/wMnUB9X "]]}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], AxesOrigin->{0.78, 0.5}, Frame->True, FrameLabel->{ FormBox[ SubscriptBox["k", "z"], TraditionalForm], FormBox["q", TraditionalForm]}, PlotRange->All, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.4409306218050003`*^9, 3.44093170685*^9}] }, Open ]] }, Open ]] }, WindowToolbars->"RulerBar", WindowSize->{804, 583}, WindowMargins->{{Automatic, -872}, {Automatic, 8}}, FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (April 27, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 178, 5, 137, "Title"], Cell[771, 30, 245, 7, 133, "Subtitle"], Cell[1019, 39, 2901, 86, 232, "Input"], Cell[3923, 127, 2770, 86, 212, "Input"], Cell[CellGroupData[{ Cell[6718, 217, 649, 20, 54, "Input"], Cell[7370, 239, 2578, 48, 240, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9985, 292, 649, 20, 54, "Input"], Cell[10637, 314, 2529, 47, 242, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13203, 366, 1033, 31, 82, "Input"], Cell[14239, 399, 6934, 118, 237, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[21210, 522, 269, 7, 31, "Input"], Cell[21482, 531, 11171, 193, 249, "Output"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)