(************** Content-type: application/mathematica ************** Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 40572, 986]*) (*NotebookOutlinePosition[ 41239, 1009]*) (* CellTagsIndexPosition[ 41195, 1005]*) (*WindowFrame->Normal*) Notebook[{ Cell["\<\ Quantum Simulations for a single Qubit System \ \>", "Subtitle"], Cell[BoxData[{ \(<< Graphics`Graphics`\), "\[IndentingNewLine]", \(<< Statistics`DataManipulation`\)}], "Input"], Cell[BoxData[{\(Off[General::spell1]\), "\n", \(Off[ General::spell]\), "\n", \(up = {1, 0};\), "\n", RowBox[{\(down = {0, 1};\), "\n"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ StyleBox[ RowBox[{\(U[\[Theta]_, \[Phi]_]\), "=", TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(Cos[\[Theta]\/2]\), \(\(-\[ExponentialE]\^\(\(-\ \[ImaginaryI]\)\ \[Phi]\)\)\ Sin[\[Theta]\/2]\)}, {\(\[ExponentialE]\^\(\[ImaginaryI]\ \[Phi]\)\ Sin[\ \[Theta]\/2]\), \(Cos[\[Theta]\/2]\)} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]}], Background->GrayLevel[1]], StyleBox[";", Background->None]}], StyleBox["\n", Background->None]}], "\[IndentingNewLine]", StyleBox[\(upX = U[Pi/2, 0] . up;\), Background->None], "\n", StyleBox[\(\(downX = \(-U[Pi/2, 0] . down\);\)\(\n\) \), Background->None], "\[IndentingNewLine]", StyleBox[\(upY = U[Pi/2, Pi/2] . up;\), Background->None], "\n", StyleBox[\(downY = I\ U[Pi/2, Pi/2] . down;\), Background->None], "\n", RowBox[{ StyleBox[\(\[Psi] = U[Pi/3, Pi/8] . up;\), Background->None], "\[IndentingNewLine]", StyleBox["\[IndentingNewLine]", Background->None]}], "\n", RowBox[{\(measurementSz[\[Psi]in_] := Module[{p, s}, p = \((Abs[Conjugate[up] . \[Psi]in])\)^2\ ; Which[s = Random[]; s \[LessEqual] p, {1, "\\>", \[Psi]out = up}\ , s > p, {\(-1\), \ "\\>", \[Psi]out = down}]]\), "\[IndentingNewLine]"}], "\n", RowBox[{\(measurementSx[\[Psi]in_] := Module[{p, s}, p = \((Abs[Conjugate[upX] . \[Psi]in])\)^2\ ; Which[s = Random[]; s \[LessEqual] p, {1, "\\>", \[Psi]out = upX}\ , s > p, {\(-1\), \ "\\>", \ \[Psi]out = downX}]]\), "\[IndentingNewLine]"}], "\n", \(measurementSy[\[Psi]in_] := Module[{p, s}, p = \((Abs[Conjugate[upY] . \[Psi]in])\)^2\ ; Which[s = Random[]; s \[LessEqual] p, {1, "\\>", \[Psi]out = upY}\ , s > p, {\(-1\), \ "\\>", \[Psi]out \ = downY}]]\), "\[IndentingNewLine]", "Null"}], "Input"], Cell[BoxData[{ StyleBox[\(The\ system\ is\ in\ a\ quantum\ state\ \ | \(\(\[Psi]\)\(>\)\ \), \ can\ we\ figure\ which\ state\ the\ system\ is\), FontSize->16], "\[IndentingNewLine]", RowBox[{ StyleBox[\(in\ using\ our\ three\ measuring\ devices\ Sz\), FontSize->16], StyleBox[",", FontSize->16], StyleBox["Sx", FontSize->16], StyleBox[",", FontSize->16], RowBox[{ StyleBox["Sy", FontSize->16], StyleBox["?", FontSize->16], StyleBox["\[IndentingNewLine]", "Text"]}]}]}], "Input", Evaluatable->False], Cell[CellGroupData[{ Cell[BoxData[{ \(measurementSz[\[Psi]]\), "\[IndentingNewLine]", \(measurementSx[\[Psi]]\), "\[IndentingNewLine]", \(measurementSy[\[Psi]]\), "\[IndentingNewLine]", \(\)}], "Input"], Cell[BoxData[ \({1, "Collapsed into |+>", {1, 0}}\)], "Output"], Cell[BoxData[ \({1, "Collapsed into |X \[UpArrow] >", {1\/\@2, 1\/\@2}}\)], "Output"], Cell[BoxData[ \({\(-1\), "Collapsed into |Y \[DownArrow]>", {\(-\(1\/\@2\)\), \ \[ImaginaryI]\/\@2}}\)], "Output"] }, Open ]], Cell[BoxData[{ \(\(samplesize = 1000;\)\), "\[IndentingNewLine]", \(\(results = Table[measurementSx[\[Psi]], {i, 1, samplesize}];\)\)}], "Input"], Cell[BoxData[ \(\(data = Table[results[\([i, 1]\)], {i, 1, samplesize}];\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Frequencies[data]\)], "Input"], Cell[BoxData[ \({{95, \(-1\)}, {905, 1}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(BarChart[%]\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations -0.181818 0.454545 0.0147151 0.000898632 [ [.27273 .00222 -6 -9 ] [.27273 .00222 6 0 ] [.72727 .00222 -3 -9 ] [.72727 .00222 3 0 ] [-0.0125 .10458 -18 -4.5 ] [-0.0125 .10458 0 4.5 ] [-0.0125 .19444 -18 -4.5 ] [-0.0125 .19444 0 4.5 ] [-0.0125 .2843 -18 -4.5 ] [-0.0125 .2843 0 4.5 ] [-0.0125 .37417 -18 -4.5 ] [-0.0125 .37417 0 4.5 ] [-0.0125 .46403 -18 -4.5 ] [-0.0125 .46403 0 4.5 ] [-0.0125 .55389 -18 -4.5 ] [-0.0125 .55389 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .27273 .01472 m .27273 .02097 L s [(-1)] .27273 .00222 0 1 Mshowa .72727 .01472 m .72727 .02097 L s [(1)] .72727 .00222 0 1 Mshowa 0 .01472 m 1 .01472 L s 0 .10458 m .00625 .10458 L s [(100)] -0.0125 .10458 1 0 Mshowa 0 .19444 m .00625 .19444 L s [(200)] -0.0125 .19444 1 0 Mshowa 0 .2843 m .00625 .2843 L s [(300)] -0.0125 .2843 1 0 Mshowa 0 .37417 m .00625 .37417 L s [(400)] -0.0125 .37417 1 0 Mshowa 0 .46403 m .00625 .46403 L s [(500)] -0.0125 .46403 1 0 Mshowa 0 .55389 m .00625 .55389 L s [(600)] -0.0125 .55389 1 0 Mshowa .125 Mabswid 0 .03269 m .00375 .03269 L s 0 .05066 m .00375 .05066 L s 0 .06863 m .00375 .06863 L s 0 .08661 m .00375 .08661 L s 0 .12255 m .00375 .12255 L s 0 .14052 m .00375 .14052 L s 0 .1585 m .00375 .1585 L s 0 .17647 m .00375 .17647 L s 0 .21241 m .00375 .21241 L s 0 .23039 m .00375 .23039 L s 0 .24836 m .00375 .24836 L s 0 .26633 m .00375 .26633 L s 0 .30228 m .00375 .30228 L s 0 .32025 m .00375 .32025 L s 0 .33822 m .00375 .33822 L s 0 .3562 m .00375 .3562 L s 0 .39214 m .00375 .39214 L s 0 .41011 m .00375 .41011 L s 0 .42809 m .00375 .42809 L s 0 .44606 m .00375 .44606 L s 0 .482 m .00375 .482 L s 0 .49998 m .00375 .49998 L s 0 .51795 m .00375 .51795 L s 0 .53592 m .00375 .53592 L s 0 .57187 m .00375 .57187 L s 0 .58984 m .00375 .58984 L s 0 .60781 m .00375 .60781 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 1 0 0 r .09091 .01472 m .09091 .32474 L .45455 .32474 L .45455 .01472 L F 0 g .5 Mabswid .09091 .01472 m .09091 .32474 L .45455 .32474 L .45455 .01472 L .09091 .01472 L s 1 0 0 r .54545 .01472 m .54545 .60332 L .90909 .60332 L .90909 .01472 L F 0 g .54545 .01472 m .54545 .60332 L .90909 .60332 L .90909 .01472 L .54545 .01472 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg0?ooo`030000003o0000o`0005d0o`0000<0 00000?ooo`3oool05`3oool001D0oooo00<000000?ooo`3oool0SP3oool00`000000o`000?l0001M 0?l000030000003oool0oooo01L0oooo000E0?ooo`800000S`3oool00`000000o`000?l0001M0?l0 00030000003oool0oooo01L0oooo000E0?ooo`030000003oool0oooo08h0oooo00<000000?l0003o 0000G@3o00000`000000oooo0?ooo`0G0?ooo`005@3oool00`000000oooo0?ooo`2>0?ooo`030000 003o0000o`0005d0o`0000<000000?ooo`3oool05`3oool001D0oooo00<000000?ooo`3oool0SP3o ool00`000000o`000?l0001M0?l000030000003oool0oooo01L0oooo000E0?ooo`030000003oool0 oooo08h0oooo00<000000?l0003o0000G@3o00000`000000oooo0?ooo`0G0?ooo`005@3oool20000 08l0oooo00<000000?l0003o0000G@3o00000`000000oooo0?ooo`0G0?ooo`005@3oool00`000000 oooo0?ooo`2>0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3oool001D0oooo 00<000000?ooo`3oool0SP3oool00`000000o`000?l0001M0?l000030000003oool0oooo01L0oooo 00030?ooo`030000003oool0oooo0080oooo0P0000040?ooo`8000001@3oool00`000000oooo0?oo o`2>0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3oool000<0oooo00D00000 0?ooo`3oool0oooo000000020?ooo`040000003oool0oooo00000080oooo00<000000?ooo`3oool0 0P3oool00`000000oooo0?ooo`2>0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool0 5`3oool000D000000P3oool010000000oooo0?ooo`0000020?ooo`040000003oool0oooo000000@0 oooo0P00002?0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3oool000040000 003oool0oooo000000<0oooo00@000000?ooo`3oool000000P3oool010000000oooo0?ooo`000004 0?ooo`030000003oool0oooo08h0oooo00<000000?l0003o0000G@3o00000`000000oooo0?ooo`0G 0?ooo`0000@0oooo0000003oool000000`3oool010000000oooo0?ooo`0000020?ooo`040000003o ool0oooo000000@0oooo00<000000?ooo`3oool0SP3oool00`000000o`000?l0001M0?l000030000 003oool0oooo01L0oooo00020?ooo`8000000`3oool010000000oooo0?ooo`0000020?ooo`040000 003oool0oooo000000@0oooo00<000000?ooo`3oool0SP3oool00`000000o`000?l0001M0?l00003 0000003oool0oooo01L0oooo00030?ooo`030000003oool0oooo0080oooo0P0000040?ooo`800000 1@3oool2000008l0oooo00<000000?l0003o0000G@3o00000`000000oooo0?ooo`0G0?ooo`005@3o ool00`000000oooo0?ooo`2>0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3o ool001D0oooo00<000000?ooo`3oool0SP3oool00`000000o`000?l0001M0?l000030000003oool0 oooo01L0oooo000E0?ooo`030000003oool0oooo08h0oooo00<000000?l0003o0000G@3o00000`00 0000oooo0?ooo`0G0?ooo`005@3oool00`000000oooo0?ooo`2>0?ooo`030000003o0000o`0005d0 o`0000<000000?ooo`3oool05`3oool001D0oooo0P00002?0?ooo`030000003o0000o`0005d0o`00 00<000000?ooo`3oool05`3oool001D0oooo00<000000?ooo`3oool0SP3oool00`000000o`000?l0 001M0?l000030000003oool0oooo01L0oooo000E0?ooo`030000003oool0oooo08h0oooo00<00000 0?l0003o0000G@3o00000`000000oooo0?ooo`0G0?ooo`005@3oool00`000000oooo0?ooo`2>0?oo o`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3oool001D0oooo00<000000?ooo`3o ool0SP3oool00`000000o`000?l0001M0?l000030000003oool0oooo01L0oooo000E0?ooo`800000 S`3oool00`000000o`000?l0001M0?l000030000003oool0oooo01L0oooo000E0?ooo`030000003o ool0oooo08h0oooo00<000000?l0003o0000G@3o00000`000000oooo0?ooo`0G0?ooo`005@3oool0 0`000000oooo0?ooo`2>0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3oool0 01D0oooo00<000000?ooo`3oool0SP3oool00`000000o`000?l0001M0?l000030000003oool0oooo 01L0oooo000E0?ooo`030000003oool0oooo08h0oooo00<000000?l0003o0000G@3o00000`000000 oooo0?ooo`0G0?ooo`005@3oool2000008l0oooo00<000000?l0003o0000G@3o00000`000000oooo 0?ooo`0G0?ooo`005@3oool00`000000oooo0?ooo`2>0?ooo`030000003o0000o`0005d0o`0000<0 00000?ooo`3oool05`3oool000040?ooo`0000000000000000@0oooo0P0000040?ooo`8000001@3o ool00`000000oooo0?ooo`2>0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3o ool000@0oooo00@000000?ooo`3oool000000P3oool010000000oooo0?ooo`0000020?ooo`030000 003oool0oooo0080oooo00<000000?ooo`3oool0SP3oool00`000000o`000?l0001M0?l000030000 003oool0oooo01L0oooo00040?ooo`040000003oool0oooo00000080oooo00@000000?ooo`3oool0 00000P3oool00`000000oooo0?ooo`020?ooo`800000S`3oool00`000000o`000?l0001M0?l00003 0000003oool0oooo01L0oooo0000103oool000000000000000030?ooo`040000003oool0oooo0000 0080oooo00@000000?ooo`3oool00000103oool00`000000oooo0?ooo`2>0?ooo`030000003o0000 o`0005d0o`0000<000000?ooo`3oool05`3oool000030?ooo`000000oooo00@0oooo00@000000?oo o`3oool000000P3oool010000000oooo0?ooo`0000040?ooo`030000003oool0oooo08h0oooo00<0 00000?l0003o0000G@3o00000`000000oooo0?ooo`0G0?ooo`0000<0oooo0000003oool0103oool0 10000000oooo0?ooo`0000020?ooo`040000003oool0oooo000000@0oooo00<000000?ooo`3oool0 SP3oool00`000000o`000?l0001M0?l000030000003oool0oooo01L0oooo00000`3oool000000000 0002000000<0oooo0P0000040?ooo`8000001@3oool00`000000oooo0?ooo`2>0?ooo`030000003o 0000o`0005d0o`0000<000000?ooo`3oool05`3oool001D0oooo0P00002?0?ooo`030000003o0000 o`0005d0o`0000<000000?ooo`3oool05`3oool001D0oooo00<000000?ooo`3oool0SP3oool00`00 0000o`000?l0001M0?l000030000003oool0oooo01L0oooo000E0?ooo`030000003oool0oooo08h0 oooo00<000000?l0003o0000G@3o00000`000000oooo0?ooo`0G0?ooo`005@3oool00`000000oooo 0?ooo`2>0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3oool001D0oooo00<0 00000?ooo`3oool0SP3oool00`000000o`000?l0001M0?l000030000003oool0oooo01L0oooo000E 0?ooo`800000S`3oool00`000000o`000?l0001M0?l000030000003oool0oooo01L0oooo000E0?oo o`030000003oool0oooo08h0oooo00<000000?l0003o0000G@3o00000`000000oooo0?ooo`0G0?oo o`005@3oool00`000000oooo0?ooo`2>0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3o ool05`3oool001D0oooo00<000000?ooo`3oool0SP3oool00`000000o`000?l0001M0?l000030000 003oool0oooo01L0oooo000E0?ooo`030000003oool0oooo08h0oooo00<000000?l0003o0000G@3o 00000`000000oooo0?ooo`0G0?ooo`005@3oool2000008l0oooo00<000000?l0003o0000G@3o0000 0`000000oooo0?ooo`0G0?ooo`005@3oool00`000000oooo0?ooo`2>0?ooo`030000003o0000o`00 05d0o`0000<000000?ooo`3oool05`3oool001D0oooo00<000000?ooo`3oool0SP3oool00`000000 o`000?l0001M0?l000030000003oool0oooo01L0oooo000E0?ooo`030000003oool0oooo08h0oooo 00<000000?l0003o0000G@3o00000`000000oooo0?ooo`0G0?ooo`005@3oool00`000000oooo0?oo o`2>0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3oool001D0oooo0P00002? 0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3oool001D0oooo00<000000?oo o`3oool0SP3oool00`000000o`000?l0001M0?l000030000003oool0oooo01L0oooo00020?ooo`80 0000103oool2000000@0oooo0P0000050?ooo`030000003oool0oooo08h0oooo00<000000?l0003o 0000G@3o00000`000000oooo0?ooo`0G0?ooo`0000D0oooo0000003oool0oooo000000020?ooo`04 0000003oool0oooo00000080oooo00@000000?ooo`3oool00000103oool00`000000oooo0?ooo`2> 0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3oool000050?ooo`000000oooo 0?ooo`0000000P3oool010000000oooo0?ooo`0000020?ooo`040000003oool0oooo000000@0oooo 0P00002?0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool05`3oool000040?ooo`00 00000000000000<0oooo00@000000?ooo`3oool000000P3oool010000000oooo0?ooo`0000040?oo o`030000003oool0oooo08h0oooo00<000000?l0003o0000G@3o00000`000000oooo0?ooo`0G0?oo o`0000<0oooo0000003oool0103oool010000000oooo0?ooo`0000020?ooo`040000003oool0oooo 000000@0oooo00<000000?ooo`3oool0SP3oool00`000000o`000?l0001M0?l000030000003oool0 oooo01L0oooo00000`3oool000000?ooo`040?ooo`040000003oool0oooo00000080oooo00@00000 0?ooo`3oool00000103oool00`000000oooo0?ooo`2>0?ooo`030000003o0000o`0005d0o`0000<0 00000?ooo`3oool05`3oool00080oooo0`0000030?ooo`800000103oool2000000D0oooo00<00000 0?ooo`3oool0SP3oool00`000000o`000?l0001M0?l000030000003oool0oooo01L0oooo000E0?oo o`800000S`3oool00`000000o`000?l0001M0?l000030000003oool0oooo01L0oooo000E0?ooo`03 0000003oool0oooo08h0oooo00<000000?l0003o0000G@3o00000`000000oooo0?ooo`0G0?ooo`00 5@3oool00`000000oooo0?ooo`2>0?ooo`030000003o0000o`0005d0o`0000<000000?ooo`3oool0 5`3oool001D0oooo00<000000?ooo`3oool0SP3oool00`000000o`000?l0001M0?l000030000003o ool0oooo01L0oooo000E0?ooo`030000003oool0oooo08h0oooo00<000000?l0003o0000G@3o0000 0`000000oooo0?ooo`0G0?ooo`005@3oool2000008l0oooo00<000000?l0003o0000G@3o00000`00 0000oooo0?ooo`0G0?ooo`005@3oool00`000000oooo0?ooo`2>0?ooo`030000003o0000o`0005d0 o`0000<000000?ooo`3oool05`3oool001D0oooo00<000000?ooo`3oool0SP3oool00`000000o`00 0?l0001M0?l000030000003oool0oooo01L0oooo000E0?ooo`030000003oool0oooo08h0ooooH@00 000I0?ooo`005@3oool00`000000oooo0?ooo`3o0?ooo`T0oooo000E0?ooo`800000o`3oool:0?oo o`005@3oool00`000000oooo0?ooo`3o0?ooo`T0oooo000E0?ooo`030000003oool0oooo0?l0oooo 2@3oool00?l0oooo8@3oool00?l0oooo8@3oool00001\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {0.223435, -60.6727, \ 0.00828076, 4.18856}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ \(mean = Sum[results[\([i, 1]\)], {i, 1, samplesize}]/ samplesize\), "\[IndentingNewLine]", \(N[%]\)}], "Input"], Cell[BoxData[ \(31\/100\)], "Output"], Cell[BoxData[ \(0.31`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\[Psi]\)], "Input"], Cell[BoxData[ \({\@3\/2, 1\/2\ \[ExponentialE]\^\(\(\[ImaginaryI]\ \[Pi]\)\/8\)}\)], "Output"] }, Open ]], Cell[BoxData[ \(\[Sigma] = mean\ of\ the\ deviation\ from\ mean\ \ \@\((\(\(<\)\(\(x\)\(-\)\)\(<\)\ \(x\)\)\(>>\))\)\^2\)], "Input", Evaluatable->False], Cell[CellGroupData[{ Cell[BoxData[ \(sigma = Sum[Sqrt[\((results[\([i, 1]\)] - mean)\)^2], {i, 1, samplesize}]/ samplesize\)], "Input"], Cell[BoxData[ \(9039\/10000\)], "Output"] }, Open ]], Cell[BoxData[{ \(Lets\ also\ measure\ the\ results\ of\ operators; \ measurementSx[\[Psi]]\), "\[IndentingNewLine]", \(\(measurementSy[\[Psi]]\ and\ perform\ the\ necessary\ \ statitics;\)\)}], "Input", Evaluatable->False], Cell[CellGroupData[{ Cell[BoxData[{ \(\(\("\";\)\(\ \[IndentingNewLine]\) \)\), "\[IndentingNewLine]", \(\(Mup = Outer[Times, up, Conjugate[up]];\)\), "\[IndentingNewLine]", \(\(Mupdagger = Transpose[Conjugate[Mup]];\)\), "\[IndentingNewLine]", \(MatrixForm[Mup]\), "\[IndentingNewLine]", \(\(Mdown = Outer[Times, down, Conjugate[down]];\)\), "\[IndentingNewLine]", \(\(Mdowndagger = Transpose[Conjugate[Mdown]];\)\), "\[IndentingNewLine]", \(MatrixForm[Mdown]\)}], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "0"}, {"0", "0"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0"}, {"0", "1"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[BoxData[ \(\(\*"\"\< \ p(a)=<\[Psi]|\!\(\(M\_+\)\^\[Dagger]\)\!\(\(M\_+\)\)|\[Psi]>\>\"";\)\)], \ "Input", Evaluatable->False], Cell[BoxData[ \(thus; \ probabilty\ for\ get\ Sz\ pointing\ up\ is\)], "Input", Evaluatable->False], Cell[CellGroupData[{ Cell[BoxData[ \(Conjugate[\[Psi]] . Mupdagger . Mup . \[Psi]\)], "Input"], Cell[BoxData[ \(3\/4\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ \(\(\("\< probability for get Sz pointing down is\>";\)\(\ \[IndentingNewLine]\) \)\), "\[IndentingNewLine]", \(Conjugate[\[Psi]] . Mdowndagger . Mdown . \[Psi]\)}], "Input"], Cell[BoxData[ \(1\/4\)], "Output"] }, Open ]], Cell[BoxData[ \(\("\";\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(MupY = Outer[Times, upY, Conjugate[upY]];\)\), "\n", \(\(MupYdagger = Transpose[Conjugate[MupY]];\)\), "\n", \(\(MdownY = Outer[Times, downY, Conjugate[downY]];\)\), "\n", \(\(MdownYdagger = Transpose[Conjugate[MdownY]];\)\), "\n", \(MatrixForm[MupY]\), "\[IndentingNewLine]", \(MatrixForm[MdownY]\)}], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(1\/2\), \(-\(\[ImaginaryI]\/2\)\)}, {\(\[ImaginaryI]\/2\), \(1\/2\)} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(1\/2\), \(\[ImaginaryI]\/2\)}, {\(-\(\[ImaginaryI]\/2\)\), \(1\/2\)} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(N[Conjugate[\[Psi]] . MupYdagger . MupY . \[Psi]]\)], "Input"], Cell[BoxData[ \(\(\(0.6657067870177957`\)\(\[InvisibleSpace]\)\) + 2.7755575615628914`*^-17\ \[ImaginaryI]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(N[Conjugate[\[Psi]] . MdownYdagger . MdownY . \[Psi]]\)], "Input"], Cell[BoxData[ \(\(\(0.3342932129822041`\)\(\[InvisibleSpace]\)\) + 0.`\ \[ImaginaryI]\)], "Output"] }, Open ]], Cell[BoxData[{ \(From\ the\ definition\ of\ the\ mean\ value\ of\ a\ random\ variable\), \ "\[IndentingNewLine]", \(\(\(\(<\)\(X\)\(>\)\(\ \)\(\[Congruent]\)\(\ \)\(\(x\_1\) p \((x\_1)\) + \ \(x\_2\) p \((x\_2)\) + \)\) ... \) . \ \ x\_n p \((x\_n)\)\[IndentingNewLine]\), "\[IndentingNewLine]", \(\(p \((x\_1)\) + \ p \((x\_2)\) + ... \) . \ \ p \((x\_n)\) = 1\)}], "Input", Evaluatable->False], Cell[BoxData[{ \(Thus\ \ \[IndentingNewLine]\), "\[IndentingNewLine]", \(\(\(<\)\(SpinupZ\)\(>\)\)\ = \ \(+1\)\ p \((\(+1\))\)\ + \ \((\(-1\))\ \)\ p \((\(-1\))\)\), "\[IndentingNewLine]", \(\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \(\(=\)\(\ \)\(\ \[Sum]\+m m\ \ < \[Psi] | \(M\_m\^\[Dagger]\) M\_m | \[Psi] >= \[Sum]\+m m < \[Psi] | M\_m | \[Psi] >= \[Sum]\+m\(\(<\)\(\[Psi]\)\) | m\ \(M\_+\) | \[Psi] > \[IndentingNewLine]now\ \ \ \(\[Sum]\+m m\ \ \ M\_m\) = M\)\)\), "\[IndentingNewLine]", \(\(\(<\)\(SpinupZ\)\(>\)\)\ \ = \ \ \(\(<\)\(\[Psi]\)\) | \[Sigma]z | \ \(\(\[Psi]\)\(>\)\)\)}], "Input", Evaluatable->False], Cell[BoxData[{ \(\(\[Sigma]z = {{1, 0}, {0, \(-1\)}};\)\), "\[IndentingNewLine]", \(\(\[Sigma]y = {{0, \(-I\)}, {I, 0}};\)\), "\[IndentingNewLine]", \(\(\[Sigma]x = {{0, 1}, {1, 0}};\)\)}], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Conjugate[\[Psi]] . \[Sigma]z . \[Psi]\)], "Input"], Cell[BoxData[ \(1\/2\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(N[Conjugate[\[Psi]] . \[Sigma]y . \[Psi]]\)], "Input"], Cell[BoxData[ \(\(\(0.3314135740355918`\)\(\[InvisibleSpace]\)\) + 0.`\ \[ImaginaryI]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(N[Conjugate[\[Psi]] . \[Sigma]x . \[Psi]]\)], "Input"], Cell[BoxData[ \(\(\(0.8001031451912655`\)\(\[InvisibleSpace]\)\) + 0.`\ \[ImaginaryI]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(psiguess = Sqrt[3/4] up + Sqrt[1/4] down\)], "Input"], Cell[BoxData[ \({\@3\/2, 1\/2}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Conjugate[psiguess] . \[Sigma]z . psiguess\)], "Input"], Cell[BoxData[ \(1\/2\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Conjugate[psiguess] . \[Sigma]y . psiguess\)], "Input"], Cell[BoxData[ \(0\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(N[Conjugate[psiguess] . \[Sigma]x . psiguess]\)], "Input"], Cell[BoxData[ \(0.8660254037844386`\)], "Output"] }, Open ]], Cell[BoxData[{ \(Lets\ Build\ a\ new\ measuring\ device\ using\ a\ linear\ combination\ \ of\ the\), "\[IndentingNewLine]", \(X\ Y\ Z\ devices\ defined\ above\), "\[IndentingNewLine]", \(\[Sigma]n = \[Sigma]z\ Cos[\[Theta]] + \[Sigma]y\ Sin[\[Theta]]\ Sin[\ \[Phi]] + \ \[Sigma]x\ Sin[\[Theta]]\ Cos[\[Phi]]\)}], "Input", Evaluatable->False], Cell[BoxData[ \(newdevice[\[Theta]_, \[Phi]_, \[Psi]in_] := Module[{p, s}, p = \((Abs[Conjugate[U[\[Theta], \[Phi]] . up] . \ \[Psi]in])\)^2\ ; Which[s = Random[]; s \[LessEqual] p, {1, "\\>", \[Psi]out = U[\[Theta], \[Phi]] . up}\ , s > p, {\(-1\), \ "\\>", \ \[Psi]out = U[\[Theta], \[Phi]] . down}]]\)], "Input"], Cell[BoxData[{ \(\[IndentingNewLine]\(\(Homwork;\)\(\[IndentingNewLine]\) \)\), "\[IndentingNewLine]", \(Try\ various\ values\ for\ \[Theta]\ \ \[Phi]\ ; \ see\ if\ you\ find\ \ \ the\ mystery\ state\ | \(\(\[Psi]\)\(>\)\(\ \[IndentingNewLine]\)\)\)}], "Input", Evaluatable->False], Cell[BoxData[{ \(Suppose\ the\ quantum\ state\ is\ \), "\[IndentingNewLine]", \(\(\(\(\(|\)\(\(\[Psi]\)\(>\)\)\)\ = \ 1/\@2 \((\(\(|\)\(+\(\(>\)\(\ \)\(-\)\)\)\(\ \)\(|\)\) \[Rule] > )\)\)\ \(\[IndentingNewLine]\) \)\), "\[IndentingNewLine]", \(\(\(Calulate\ the\ uncertainties\)\(\ \)\(\[IndentingNewLine]\) \)\), "\[IndentingNewLine]", \(\[CapitalDelta]\ \[Sigma]z, \ \[CapitalDelta]\ \[Sigma]y\ , \ \ \[CapitalDelta]\ \[Sigma]z, \[IndentingNewLine]\)}], "Input", Evaluatable->False], Cell[BoxData[{ \(On\ the\ same\ state, \ consider\ the\ following\ cascade\ \(measurement\ : \[IndentingNewLine]\ \[IndentingNewLine]first\ \ measure\ \[Sigma]z\), \ then\ immediately\ after\ \[Sigma]y, \ followed\ by\ \[Sigma]z\ again, \[IndentingNewLine]\[IndentingNewLine]Do\ \ the\ numerical\ experiment\ and\ compare\ the\ probabilties\ obtained, \ and\ mean\ value\ for\ this\), "\[IndentingNewLine]", \(measurement . \ Can\ you\ predict\ these\ results\ without\ the\ \ \(\(simulation\)\(?\)\(\ \ \)\)\)}], "Input", Evaluatable->False] }, FrontEndVersion->"4.1 for X", ScreenRectangle->{{0, 1280}, {0, 1024}}, WindowSize->{789, 931}, WindowMargins->{{233, Automatic}, {Automatic, 22}}, StyleDefinitions -> "Classroom.nb" ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[1705, 50, 78, 3, 66, "Subtitle"], Cell[1786, 55, 121, 2, 63, "Input"], Cell[1910, 59, 2526, 56, 505, "Input"], Cell[4439, 117, 638, 20, 100, "Input", Evaluatable->False], Cell[CellGroupData[{ Cell[5102, 141, 200, 4, 95, "Input"], Cell[5305, 147, 67, 1, 47, "Output"], Cell[5375, 150, 89, 1, 68, "Output"], Cell[5467, 153, 125, 3, 68, "Output"] }, Open ]], Cell[5607, 159, 163, 3, 63, "Input"], Cell[5773, 164, 91, 1, 47, "Input"], Cell[CellGroupData[{ Cell[5889, 169, 50, 1, 47, "Input"], Cell[5942, 172, 58, 1, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6037, 178, 44, 1, 47, "Input"], Cell[6084, 181, 25100, 484, 186, 2703, 203, "GraphicsData", "PostScript", \ "Graphics"], Cell[31187, 667, 130, 3, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[31354, 675, 151, 4, 63, "Input"], Cell[31508, 681, 41, 1, 64, "Output"], Cell[31552, 684, 39, 1, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[31628, 690, 39, 1, 47, "Input"], Cell[31670, 693, 105, 2, 68, "Output"] }, Open ]], Cell[31790, 698, 165, 4, 56, "Input", Evaluatable->False], Cell[CellGroupData[{ Cell[31980, 706, 133, 3, 47, "Input"], Cell[32116, 711, 45, 1, 64, "Output"] }, Open ]], Cell[32176, 715, 238, 5, 63, "Input", Evaluatable->False], Cell[CellGroupData[{ Cell[32439, 724, 557, 11, 159, "Input"], Cell[32999, 737, 227, 7, 81, "Output"], Cell[33229, 746, 227, 7, 81, "Output"] }, Open ]], Cell[33471, 756, 137, 4, 49, "Input", Evaluatable->False], Cell[33611, 762, 105, 2, 47, "Input", Evaluatable->False], Cell[CellGroupData[{ Cell[33741, 768, 77, 1, 47, "Input"], Cell[33821, 771, 38, 1, 64, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[33896, 777, 202, 4, 79, "Input"], Cell[34101, 783, 38, 1, 64, "Output"] }, Open ]], Cell[34154, 787, 93, 1, 47, "Input"], Cell[CellGroupData[{ Cell[34272, 792, 359, 6, 127, "Input"], Cell[34634, 800, 276, 7, 97, "Output"], Cell[34913, 809, 276, 7, 97, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[35226, 821, 82, 1, 47, "Input"], Cell[35311, 824, 131, 2, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[35479, 831, 86, 1, 47, "Input"], Cell[35568, 834, 110, 2, 47, "Output"] }, Open ]], Cell[35693, 839, 461, 9, 95, "Input", Evaluatable->False], Cell[36157, 850, 710, 12, 173, "Input", Evaluatable->False], Cell[36870, 864, 212, 3, 79, "Input"], Cell[CellGroupData[{ Cell[37107, 871, 71, 1, 47, "Input"], Cell[37181, 874, 38, 1, 64, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[37256, 880, 74, 1, 47, "Input"], Cell[37333, 883, 110, 2, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[37480, 890, 74, 1, 47, "Input"], Cell[37557, 893, 110, 2, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[37704, 900, 73, 1, 47, "Input"], Cell[37780, 903, 48, 1, 68, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[37865, 909, 75, 1, 47, "Input"], Cell[37943, 912, 38, 1, 64, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[38018, 918, 75, 1, 47, "Input"], Cell[38096, 921, 35, 1, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[38168, 927, 78, 1, 47, "Input"], Cell[38249, 930, 53, 1, 47, "Output"] }, Open ]], Cell[38317, 934, 357, 6, 79, "Input", Evaluatable->False], Cell[38677, 942, 474, 9, 79, "Input"], Cell[39154, 953, 302, 6, 111, "Input", Evaluatable->False], Cell[39459, 961, 523, 10, 150, "Input", Evaluatable->False], Cell[39985, 973, 583, 11, 127, "Input", Evaluatable->False] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)