(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 4.0, MathReader 4.0, or any compatible application. 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[ 78502, 2351]*) (*NotebookOutlinePosition[ 114445, 3543]*) (* CellTagsIndexPosition[ 114030, 3525]*) (*WindowFrame->Normal*) Notebook[{ Cell[TextData[StyleBox["\"Superlumial\" phenomena in linear optics", \ "Title"]], "Text", PageWidth->WindowWidth, TextJustification->0], Cell["\<\ S. M. Rochester and D. Budker UC Berkeley Physics Dept. Sept. 9, 2002\ \>", "Author"], Cell[CellGroupData[{ Cell["Introduction", "Subtitle"], Cell[CellGroupData[{ Cell["Background", "Section 1"], Cell[TextData[{ "Recently, both the ", ButtonBox["scientific", ButtonData:>{ URL[ "http://focus.aps.org/v3/st37.html"], None}, ButtonStyle->"Hyperlink"], " and popular (e.g., ", ButtonBox["NYT", ButtonData:>{ URL[ "http://www.nytimes.com/library/national/science/021899sci-slow-light.\ html"], None}, ButtonStyle->"Hyperlink"], " and ", ButtonBox["NPR", ButtonData:>{ URL[ "http://search.npr.org/cf/cmn/segment_display.cfm?segID=79717"], None}, ButtonStyle->"Hyperlink"], ") media have paid much attention to optical phenomena dubbed \"slow\" and \ \"fast\" light. When a pulse of weak probe light propagates through a \ nonlinear medium in the presence of a strong drive field, the light pulse \ can, under certain circumstances, be transmitted without much distortion, and \ the ", StyleBox["apparent", FontSlant->"Italic"], " group velocity can be very slow (\[Tilde]1 ", Cell[BoxData[ FormBox[ RowBox[{ StyleBox["m", FontSlant->"Plain"], "\[VeryThinSpace]\[VeryThinSpace]", SuperscriptBox[ StyleBox["s", FontSlant->"Plain"], \(-1\)]}], TraditionalForm]]], " in some experiments), faster than ", Cell[BoxData[ \(TraditionalForm\`c\)]], ", or negative (see recent comprehensive reviews by ", ButtonBox["Matsko et. al., 2001", ButtonData:>"Mat2001", ButtonStyle->"Hyperlink"], "; ", ButtonBox["Boyd and Gauthier, 2002", ButtonData:>"Boy2002", ButtonStyle->"Hyperlink"], "). Although the recent studies of these phenomena have been in the \ context of ", StyleBox["nonlinear", FontSlant->"Italic"], " optics, in which at least two light fields interact with an atomic \ medium, these effects can be understood by considering the simpler case of ", StyleBox["linear", FontSlant->"Italic"], " optics." }], "Text", CellMargins->{{Inherited, 0}, {Inherited, Inherited}}, CellSize->{7, Inherited}], Cell[TextData[{ "In this tutorial, we use considerations from linear optics to show how the \ variation in group velocity comes about. Although the physics is quite \ straightforward, and has been worked out in detail both theoretically (", ButtonBox["Garrett and McCumber, 1970", ButtonData:>"Gar70", ButtonStyle->"Hyperlink"], ") and experimentally (", ButtonBox["Chu and Wong, 1982", ButtonData:>"Chu82", ButtonStyle->"Hyperlink"], "), we attempt to remove any remaining mystery from the subject by \ discussing it in very basic language and illustrating each step with graphics \ derived directly from the discussion. (See also ", ButtonBox["Sprangle et. al., 2001", ButtonData:>"Spr2001", ButtonStyle->"Hyperlink"], "; ", ButtonBox["McDonald, 2001", ButtonData:>"McD2001", ButtonStyle->"Hyperlink"], "; ", ButtonBox["Dogariu et. al. 2001", ButtonData:>"Dog2001", ButtonStyle->"Hyperlink"], ")." }], "Text"] }, Open ]], Cell[CellGroupData[{ Cell["Approach", "Section 1"], Cell[TextData[{ "We consider an optical medium consisting of an atomic vapor, with light \ whose frequency is tuned near an isolated resonance. We introduce the complex \ refractive index, and first describe the propagation of a monochromatic light \ wave through a slab of the medium. The imaginary part of the refractive index \ leads to light absorption, while the real part affects the light's phase. \ Next, we consider a smooth (Gaussian) pulse incident on the slab of our \ atomic medium. We decompose the pulse into Fourier components, and analyze \ the propagation of each of these components through the medium separately. \ This is justified in the case of linear medium whose optical properties are \ independent of the presence of light. Finally, we recombine the resulting \ Fourier components into a \"new\" pulse, and analyze its spatial and temporal \ distribution. Clearly, this discussion does not contain any physical or \ mathematical concepts beyond those found in any introductory physics course. \ Nevertheless, this will be sufficient to illustrate all of our ", ButtonBox["conclusions", ButtonData:>"conclusions", ButtonStyle->"Hyperlink"], "." }], "Text"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Mathematica Setup", "Subtitle", InitializationCell->True], Cell["Turn off annoying warnings", "Text", InitializationCell->True], Cell[BoxData[{ \(Off[General::"\"]\), "\[IndentingNewLine]", \(Off[General::"\"]\), "\[IndentingNewLine]", \(Off[Solve::"\"]\)}], "InputOnly", InitializationCell->True], Cell["Load packages we'll use", "Text", InitializationCell->True], Cell[BoxData[{ \(\(Needs["\"];\)\), "\[IndentingNewLine]", \(\(Needs["\"];\)\)}], "InputOnly", InitializationCell->True], Cell["Set options for plotting", "Text", InitializationCell->True], Cell[BoxData[ RowBox[{ RowBox[{"SetOptions", "[", RowBox[{ "Plot", ",", \(PlotRange \[Rule] All\), ",", \(Frame \[Rule] True\), ",", \(Axes \[Rule] {True, False}\), ",", \(ImageSize \[Rule] {6\ 72, \ Automatic}\), ",", \(Ticks \[Rule] None\), ",", \(FrameTicks \[Rule] None\), ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\\\"Italic\\\"]\)\>\"", ",", "\"\\\\"Italic\\\"]\)\>\""}], "}"}]}]}], "]"}], ";"}]], "InputOnly", InitializationCell->True], Cell[BoxData[ \(\(SetOptions[ListPlot, PlotRange \[Rule] All, Frame \[Rule] True, Axes \[Rule] {True, False}, PlotJoined \[Rule] True, ImageSize \[Rule] {6\ 72, \ Automatic}, Ticks \[Rule] None, FrameTicks \[Rule] None];\)\)], "InputOnly", InitializationCell->True], Cell[BoxData[ \(\($TextStyle\ = \ {FontFamily -> "\", FontSize \[Rule] 20};\)\)], "InputOnly", InitializationCell->True], Cell["Set ComplexExpand to behave the way we like", "Text", InitializationCell->True], Cell[BoxData[ \(\(SetOptions[ComplexExpand, TargetFunctions \[Rule] {Re, Im}];\)\)], "InputOnly", InitializationCell->True] }, Closed]], Cell[CellGroupData[{ Cell["Discussion of light propagation in linear media", "Subtitle"], Cell["Monochromatic waves, phase velocity", "Section 1"], Cell[CellGroupData[{ Cell["Definition of the light field", "Section"], Cell[TextData[{ "Any plane light wave propagating in the ", Cell[BoxData[ \(TraditionalForm\`z\)]], " direction can be decomposed into a number of sinusoidal waves. We use the \ complex notation in which each field is described by an amplitude times the \ real part of" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(EField[ k_, \[Omega]_] = \[ExponentialE]\^\(\[ImaginaryI] \((k\ z\ - \ \ \[Omega]\ t)\)\)\)], "Input", InitializationCell->True], Cell[BoxData[ \(\[ExponentialE]\^\(\[ImaginaryI]\ \((k\ z - t\ \[Omega])\)\)\)], \ "Output"] }, Open ]], Cell["where", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(k[n_, \[Omega]_] = \(n\ \[Omega]\)\/c\)], "Input", InitializationCell->True], Cell[BoxData[ \(\(n\ \[Omega]\)\/c\)], "Output"] }, Open ]], Cell[TextData[{ "is the wavenumber, which depends on the light frequency ", Cell[BoxData[ \(TraditionalForm\`\[Omega]\)]], " and the index of refraction of the medium ", Cell[BoxData[ \(TraditionalForm\`n(\[Omega])\)]], ". ", Cell[BoxData[ \(TraditionalForm\`c\)]], " is the speed of light in vacuum. The real field is given by" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(ReEField[k_, \[Omega]_] = ComplexExpand[Re[EField[k, \[Omega]]], k]\)], "Input", InitializationCell->True], Cell[BoxData[ \(\[ExponentialE]\^\(\(-z\)\ Im[k]\)\ Cos[ t\ \[Omega] - z\ Re[k]]\)], "Output"] }, Open ]], Cell[BoxData[ \(\(Table[ Plot[ReEField[2 \[Pi], 2 \[Pi]], {z, 0, 8}, PlotRange \[Rule] {{0, 8}, {\(-1.25\), 1.25}}], {t, 0, .9, .1}];\)\)], "Input"], Cell[TextData[{ StyleBox["A monochromatic wave ", FontColor->GrayLevel[0]], "\n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Phase velocity \[LongDash] general medium\ \>", "Section"], Cell[TextData[{ "The ", StyleBox["phase velocity", FontSlant->"Italic"], " is the velocity of a point of constant phase of the electric field, for \ example, a zero crossing. To find the phase velocity we solve for ", Cell[BoxData[ \(TraditionalForm\`z\)]], " such that ", Cell[BoxData[ \(TraditionalForm\`Re[E(z)] = 0\)]], " as a function of time:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(zConstPhase = \(Solve[ReEField[k, \[Omega]] \[Equal] 0, z]\)\[LeftDoubleBracket]1\[RightDoubleBracket]\)], "Input", InitializationCell->True], Cell[BoxData[ \({z \[Rule] \(\(-\(\[Pi]\/2\)\) + t\ \[Omega]\)\/Re[k]}\)], "Output"] }, Open ]], Cell[TextData[{ "The phase velocity is given by ", Cell[BoxData[ \(TraditionalForm\`\(d\ z\)\/\(d\ t\)\)]] }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Vph[k_, \[Omega]_] = D[z /. zConstPhase, t]\)], "Input", InitializationCell->True], Cell[BoxData[ \(\[Omega]\/Re[k]\)], "Output"] }, Open ]], Cell["In terms of the index of refraction", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(VphOfN[n_] = ComplexExpand[Vph[k[n, \[Omega]], \[Omega]], n]\)], "Input",\ InitializationCell->True], Cell[BoxData[ \(c\/Re[n]\)], "Output"] }, Open ]], Cell[TextData[{ "In a vacuum, ", Cell[BoxData[ \(TraditionalForm\`n = 1\)]], " and the phase velocity is " }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(VphOfN[1]\)], "Input"], Cell[BoxData[ \(c\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Wave propagation through linear atomic media", "Section"], Cell[TextData[{ "In a linear medium near an atomic resonance of frequency ", Cell[BoxData[ \(TraditionalForm\`\[Omega]0\)]], " and width ", Cell[BoxData[ \(TraditionalForm\`\[CapitalGamma]\)]], ", the index of refraction is given by (see, for example, ", ButtonBox["Lipson, 1995, Ch. 13.3.2", ButtonData:>"Lip95", ButtonStyle->"Hyperlink"], "; ", ButtonBox["Griffiths, 1999, Ch. 8.4", ButtonData:>"Gri99", ButtonStyle->"Hyperlink"], ")" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(nMedium[\[Omega]_] = 1 - \(\(\[Alpha]\ c\)\/\(\[Omega] - \[Omega]0 + \[ImaginaryI]\ \ \[CapitalGamma]/2\)\) \[CapitalGamma]\/\(4\ \[Omega]0\)\)], "Input", InitializationCell->True], Cell[BoxData[ \(1 - \(c\ \[Alpha]\ \[CapitalGamma]\)\/\(4\ \((\(\[ImaginaryI]\ \ \[CapitalGamma]\)\/2 + \[Omega] - \[Omega]0)\)\ \[Omega]0\)\)], "Output"] }, Open ]], Cell[TextData[{ "where ", Cell[BoxData[ \(TraditionalForm\`\[Alpha]\)]], " is the absorption coefficient on resonance. The real part of the index of \ refraction, responsible for phase shifts, is given by" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(RenMedium[\[Omega]_] = Collect[ComplexExpand[ Re[nMedium[\[Omega]]]] /. \(c\ \[Alpha]\ \[CapitalGamma]\)\/\(\(\ \ \)\((\[CapitalGamma]\^2\/4 + \((\[Omega] - \[Omega]0)\)\^2)\)\) \[Rule] 4\ \[Xi]\ , \[Xi]]\)], "Input", InitializationCell->True], Cell[BoxData[ \(1 + \[Xi]\ \((1 - \[Omega]\/\[Omega]0)\)\)], "Output"] }, Open ]], Cell[TextData[{ "where we have defined the parameter ", Cell[BoxData[ \(TraditionalForm\`\[Xi]\)]], " to be" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\[Xi][\[Omega]_] = \(c\ \[Alpha]\ \[CapitalGamma]\)\/\(\(\ \)\(4 \((\ \[CapitalGamma]\^2\/4 + \((\[Omega] - \[Omega]0)\)\^2)\)\)\)\)], "Input", InitializationCell->True], Cell[BoxData[ \(\(c\ \[Alpha]\ \[CapitalGamma]\)\/\(4\ \((\[CapitalGamma]\^2\/4 + \((\ \[Omega] - \[Omega]0)\)\^2)\)\)\)], "Output"] }, Open ]], Cell[TextData[{ "In a physical medium, the real and the imaginary parts of the refractive \ index are related through the Kramers - Kronig relations (", ButtonBox["Landau et. al., 1995, Sec. 82", ButtonData:>"Lan95", ButtonStyle->"Hyperlink"], "). In a simulation, we can easily create an artificial medium that does \ not posess this property; however, our results immediately become \ nonphysical, and causality is violated." }], "Commentary"], Cell["\<\ The following function plots the real and imaginary parts of the index of \ refraction vs. frequency:\ \>", "Text"], Cell[BoxData[ RowBox[{\(IndexPlot[Parameters_, opts___]\), ":=", RowBox[{"Plot", "[", RowBox[{\(Evaluate[\(Join[{Re[nMedium[\[Omega]]] - 1, Im[nMedium[\[Omega]]]}, If[\(PlotEnvelope /. {opts}\) /. PlotEnvelope \[Rule] False, {\[CapitalGamma]p/ 7\ Envelope[\[CapitalDelta]1\ \[CapitalGamma]]}, \ {}]] /. \[Omega] \[Rule] \[Omega]0 + \[CapitalDelta]1\ \[CapitalGamma]\) /. Parameters]\), ",", \({\[CapitalDelta]1, \(-4\)\ , 4}\), ",", \(PlotStyle \[Rule] {RGBColor[1, 0, 0], RGBColor[0, 1, 0], RGBColor[0, 0, 1]}\), ",", \(DisplayFunction \[Rule] $DisplayFunction\), ",", \(FrameLabel \[Rule] {"\", "\"}\), ",", \(FrameTicks \[Rule] {True, False, False, False}\), ",", RowBox[{"PlotLegend", "\[Rule]", RowBox[{"{", RowBox[{"\"\\\\"Italic\\\"]\) - 1\>\"", ",", "\"\\\\"Italic\\\"]\)\>\"", ",", \(StyleForm["\", FontSize \[Rule] 12]\)}], "}"}]}], ",", \(LegendPosition \[Rule] { .23, .2}\), ",", \(LegendSize \[Rule] { .6, .25}\), ",", \(LegendShadow \[Rule] { .015, \(- .015\)}\), ",", \(Axes \[Rule] True\)}], "]"}]}]], "InputOnly", InitializationCell->True], Cell[BoxData[ \(\(IndexPlot[{c \[Rule] 1, \[Alpha] \[Rule] 1, \[CapitalGamma] \[Rule] 2, \[Omega]0 \[Rule] 2 \[Pi]}];\)\)], "InputOutdent"], Cell[TextData[StyleBox["The real and imaginary parts of the index of \ refraction in a medium as a function of frequency.", FontColor->GrayLevel[0]]], "CaptionOutdent"], Cell["The phase velocity in the medium is given by", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(VphMedium[\[Omega]_] = Vph[k[RenMedium[\[Omega]], \[Omega]], \[Omega]] // ComplexExpand\)], "Input", InitializationCell->True], Cell[BoxData[ \(c\/\(1 + \[Xi]\ \((1 - \[Omega]\/\[Omega]0)\)\)\)], "Output"] }, Open ]], Cell[TextData[{ "Since ", Cell[BoxData[ \(TraditionalForm\`\[Xi]\)]], " is always positive, we see that the phase velocity is ", Cell[BoxData[ \(TraditionalForm\`c\)]], " on resonance, less than ", Cell[BoxData[ \(TraditionalForm\`\(\(\ \)\(c\)\)\)]], " below resonance, and ", StyleBox["greater", FontSlant->"Italic"], " than ", Cell[BoxData[ \(TraditionalForm\`c\)]], " above resonance. To illustrate this, we can plot the light electric field \ as it traverses a medium. The wavenumber in the medium is given by" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(kMedium[\[Omega]_] = k[nMedium[\[Omega]], \[Omega]] // Simplify\)], "Input", InitializationCell->True], Cell[BoxData[ \(\[Omega]\/c - \(\[Alpha]\ \[CapitalGamma]\ \[Omega]\)\/\(4\ \((\(\ \[ImaginaryI]\ \[CapitalGamma]\)\/2 + \[Omega] - \[Omega]0)\)\ \ \[Omega]0\)\)], "Output"] }, Open ]], Cell["and the electric field is given by", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(EFieldMedium[z_, t_, \[Omega]_] = EField[kMedium[\[Omega]], \[Omega]] // Simplify\)], "Input", InitializationCell->True], Cell[BoxData[ \(\[ExponentialE]\^\(\(-\(1\/4\)\)\ \[ImaginaryI]\ \[Omega]\ \((4\ t - \ \(4\ z\)\/c + \(z\ \[Alpha]\ \[CapitalGamma]\)\/\(\((\(\[ImaginaryI]\ \ \[CapitalGamma]\)\/2 + \[Omega] - \[Omega]0)\)\ \[Omega]0\))\)\)\)], "Output"] }, Open ]], Cell[TextData[{ "The light intensity is proportional to ", Cell[BoxData[ \(TraditionalForm\`Abs[E\^2]\)]], ", where Abs stands for absolute value. We can verify that the intensity \ has the correct dependence on the absorption coefficient:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Abs[EFieldMedium[z, t, \[Omega]0]\^2] // ComplexExpand\) // PowerExpand\)], "Input"], Cell[BoxData[ \(\[ExponentialE]\^\(\(-z\)\ \[Alpha]\)\)], "Output"] }, Open ]], Cell["\<\ The following function produces an animation of a monochromatic wave \ traversing a medium. \ \>", "Text"], Cell[BoxData[ \(PlotMonochromatic[{z0_, z1_, z2_, z3_}, {tMin_, tMax_, tStep_}, Parameters_] := \((\[IndentingNewLine]SetOptions[Plot, MaxBend \[Rule] 6, DisplayFunction \[Rule] Identity]; \[IndentingNewLine]\[IndentingNewLine]FieldOptions = PlotStyle \[Rule] {AbsoluteThickness[ 2]}; \[IndentingNewLine]GhostFieldOptions = PlotStyle \[Rule] {AbsoluteThickness[1.5], AbsoluteDashing[{2, 5}]}; \[IndentingNewLine]\[IndentingNewLine]Region1Field[ z_, t_] = \(EFieldMedium[z, t, \[Omega]] /. {\[Alpha] \[Rule] 0}\) /. Parameters; \[IndentingNewLine]Region2Field[z_, t_] = Region1Field[z1, 0] EFieldMedium[z - z1, t, \[Omega]] /. Parameters; Region3Field[z_, t_] = \(Region2Field[z2, 0] EFieldMedium[z - z2, t, \[Omega]] /. {\[Alpha] \[Rule] 0}\) /. Parameters; \[IndentingNewLine]\[IndentingNewLine]Table[\ \[IndentingNewLine]Show[\[IndentingNewLine]Plot[ Re[Region1Field[z, t]], {z, z1, z3}, Evaluate[GhostFieldOptions]], \[IndentingNewLine]Plot[ Re[Region1Field[z, t]], {z, z0, z1}, Evaluate[FieldOptions]], \[IndentingNewLine]Plot[ Re[Region2Field[z, t]], {z, z1, z2}, Evaluate[FieldOptions]], \[IndentingNewLine]Plot[ Re[Region3Field[z, t]], {z, z2, z3}, Evaluate[ FieldOptions]], \[IndentingNewLine]DisplayFunction \[Rule] \ $DisplayFunction, Epilog \[Rule] {AbsoluteThickness[3], Line[{{z1, \(-1.4\)}, {z1, 1.4}, {z2, 1.4}, {z2, \(-1.4\)}, {z1, \(-1.4\)}}]}, PlotRange \[Rule] {{z0 - 0.000001, z3}, {\(-1.5\), 1.5}}\[IndentingNewLine]], \[IndentingNewLine]{t, tMin, tMax - tStep, tStep}]; \[IndentingNewLine]\[IndentingNewLine]SetOptions[Plot, MaxBend \[Rule] 10, DisplayFunction \[Rule] \ $DisplayFunction];\[IndentingNewLine])\)\)], "InputOnly", InitializationCell->True], Cell[TextData[{ "Here is a plot of the electric field for a light wave tuned above \ resonance. The phase of the wave at the output is advanced relative to what \ it would have been, had there been no medium (dotted line). The amplitude of \ the wave is reduced because of absorption (imaginary part of the index of \ refraction). Try changing the parameters and re-evaluating. Tuning the light \ frequency closer to resonance reduces the phase shift and increases \ absorption. Tuning below resonance (i.e. making ", Cell[BoxData[ \(TraditionalForm\`\(\(\[Omega] < \[Omega]0\)\()\)\)\)]], ", one sees that the phase of the wave at the output is retarded relative \ to what it would have been, had there been no medium (dotted line) because \ the phase velocity is ", Cell[BoxData[ \(TraditionalForm\`\(\(<\)\(c\)\)\)]], " in this case." }], "Text"], Cell[BoxData[ \(PlotMonochromatic[{0, 3, 5, 8}, {0. , 0.95, .05}, {c \[Rule] 1, \[Alpha] \[Rule] 3. , \[Omega] \[Rule] 2.1 \[Pi], \[CapitalGamma] \[Rule] 1, \[Omega]0 \[Rule] 2 \[Pi]}]\)], "Input"], Cell[TextData[{ StyleBox["A monochromatic wave propagating through a linear medium.", FontColor->GrayLevel[0]], "\n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"], Cell["", "Text"] }, Open ]], Cell["Light pulses, group velocity", "Section 1"], Cell[CellGroupData[{ Cell["Constructing a light pulse", "Section"], Cell[CellGroupData[{ Cell["Introduction", "Subsection"], Cell[TextData[{ "To construct a Gaussian pulse of width ", Cell[BoxData[ \(TraditionalForm\`\[CapitalGamma]p\)]], ", we superpose waves with amplitudes given by " }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Envelope[\[CapitalDelta]_] = \(\(1\/\(\[CapitalGamma]p \@\( 2 \ \[Pi]\)\)\) \(\[ExponentialE]\^\(-\(\[CapitalDelta]\^2\/\(2 \[CapitalGamma]p\ \^2\)\)\)\)\(\ \)\)\)], "Input", InitializationCell->True], Cell[BoxData[ \(\[ExponentialE]\^\(-\(\[CapitalDelta]\^2\/\(2\ \ \[CapitalGamma]p\^2\)\)\)\/\(\@\(2\ \[Pi]\)\ \[CapitalGamma]p\)\)], "Output"] }, Open ]], Cell[TextData[{ "where ", Cell[BoxData[ \(TraditionalForm\`\[CapitalDelta]\)]], " is the detuning from the central frequency \[Omega], and the first factor \ provides overall normalization. The more waves of different frequencies we \ add, the more the resultant envelope looks like a pulse, as we will now \ show." }], "Text"] }, Open ]], Cell[CellGroupData[{ Cell["Discrete frequencies", "Subsection"], Cell[CellGroupData[{ Cell["Two Frequencies", "Subsubsection"], Cell[TextData[{ "For illustration, let us look at waves traveling in vacuum, in which ", Cell[BoxData[ \(TraditionalForm\`n = 1\)]], ". A \"pulse\" with two frequencies is given by" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(TwoFreq = \(Sum[ Envelope[\[CapitalDelta]] EField[k[ 1, \[Omega] + \[CapitalDelta]], \[Omega] + \ \[CapitalDelta]], {\[CapitalDelta], \(-3\) \[Pi]/2, 3 \[Pi]/2, 3 \[Pi]}] // ComplexExpand\) // FullSimplify\)], "Input"], Cell[BoxData[ \(\(\[ExponentialE]\^\(\(-\(\(9\ \[Pi]\^2\)\/\(8\ \[CapitalGamma]p\^2\)\)\ \) - \[ImaginaryI]\ \((t - z\/c)\)\ \[Omega]\)\ \@\(2\/\[Pi]\)\ Cos[\(3\ \ \[Pi]\ \((c\ t - z)\)\)\/\(2\ c\)]\)\/\[CapitalGamma]p\)], "Output"] }, Open ]], Cell["The envelope is given by the absolute value of the field", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(TwoFreqEnvelope = \(\(Abs[TwoFreq] // ComplexExpand\) // Simplify\) // PowerExpand\)], "Input"], Cell[BoxData[ \(\(\[ExponentialE]\^\(-\(\(9\ \[Pi]\^2\)\/\(8\ \[CapitalGamma]p\^2\)\)\)\ \ \@\(2\/\[Pi]\)\ Cos[\(3\ \[Pi]\ \((c\ t - z)\)\)\/\(2\ c\)]\)\/\ \[CapitalGamma]p\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Re[TwoFreq /. c \[Rule] 1] // ComplexExpand\)], "Input"], Cell[BoxData[ \(\(\[ExponentialE]\^\(-\(\(9\ \[Pi]\^2\)\/\(8\ \[CapitalGamma]p\^2\)\)\)\ \ \@\(2\/\[Pi]\)\ Cos[3\/2\ \[Pi]\ \((t - z)\)]\ Cos[\((\(-t\) + z)\)\ \ \[Omega]]\)\/\[CapitalGamma]p\)], "Output"] }, Open ]], Cell["\<\ Here is the plot. The sinusoidal envelope is an indication of the pulse-like \ behavior that will come when we add waves of more frequencies.\ \>", "Text"], Cell[BoxData[ \(\(Table[ Plot[Evaluate[ Re[TwoFreq /. {c \[Rule] 1, \[Omega] \[Rule] 15 \[Pi], \[CapitalGamma]p \[Rule] \[Pi]}]], {z, \(-8\)/\ \((2 \[Pi])\), 8/\((2 \[Pi])\)}, PlotRange \[Rule] {{\(-8\)/\((2 \[Pi])\), 8/\((2 \[Pi])\)}, {\(- .1\), .1}}, PlotPoints \[Rule] 70], {t, 0, 4/3 - 1/9, 1/9}];\)\)], "Input"], Cell[TextData[{ "Superposition of two monochromatic waves \n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"] }, Open ]], Cell[CellGroupData[{ Cell["Three frequencies", "Subsubsection"], Cell["Superposing waves of three frequencies, we get", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(ThreeFreq = \(Sum[ Envelope[\[CapitalDelta]] EField[k[ 1, \[Omega] + \[CapitalDelta]], \[Omega] + \ \[CapitalDelta]], {\[CapitalDelta], \(-\[Pi]\), \[Pi], \[Pi]}] // ComplexExpand\) // FullSimplify\)], "Input"], Cell[BoxData[ \(\(\[ExponentialE]\^\(\(-\(\[Pi]\^2\/\(2\ \[CapitalGamma]p\^2\)\)\) - \ \[ImaginaryI]\ \((t - z\/c)\)\ \[Omega]\)\ \ \((\[ExponentialE]\^\(\[Pi]\^2\/\(2\ \[CapitalGamma]p\^2\)\) + 2\ Cos[\[Pi]\ \ \((t - z\/c)\)])\)\)\/\(\@\(2\ \[Pi]\)\ \[CapitalGamma]p\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(Abs[ThreeFreq] // ComplexExpand\) // Simplify\) // PowerExpand\)], "Input"], Cell[BoxData[ \(\(\[ExponentialE]\^\(-\(\[Pi]\^2\/\(2\ \[CapitalGamma]p\^2\)\)\)\ \((\ \[ExponentialE]\^\(\[Pi]\^2\/\(2\ \[CapitalGamma]p\^2\)\) + 2\ Cos[\[Pi]\ \ \((t - z\/c)\)])\)\)\/\(\@\(2\ \[Pi]\)\ \[CapitalGamma]p\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(Re[ThreeFreq] // ComplexExpand\) // FullSimplify\) // PowerExpand\)], "Input"], Cell[BoxData[ \(\(\[ExponentialE]\^\(-\(\[Pi]\^2\/\(2\ \[CapitalGamma]p\^2\)\)\)\ \((\ \[ExponentialE]\^\(\[Pi]\^2\/\(2\ \[CapitalGamma]p\^2\)\) + 2\ Cos[\[Pi]\ \ \((t - z\/c)\)])\)\ Cos[\((\(-t\) + z\/c)\)\ \[Omega]]\)\/\(\@\(2\ \[Pi]\)\ \ \[CapitalGamma]p\)\)], "Output"] }, Open ]], Cell[BoxData[ \(\(Table[ Plot[Evaluate[ Re[ThreeFreq /. {c \[Rule] 1, \[Omega] \[Rule] 20 \[Pi], \[CapitalGamma]p \[Rule] \[Pi]}]], {z, \ \(-12\)/\((2 \[Pi])\), 12/\((2 \[Pi])\)}, PlotRange \[Rule] {{\(-12\)/\((2 \[Pi])\), 12/\((2 \[Pi])\)}, {\(- .3\), .3}}, PlotPoints \[Rule] 70], {t, 0, 2 - 1/9, 1/9}];\)\)], "Input"], Cell[TextData[{ "Superposition of three monochromatic waves \n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"] }, Open ]], Cell[CellGroupData[{ Cell["Four frequencies", "Subsubsection"], Cell["With four frequencies:", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(FourFreq = Sum[Envelope[\[CapitalDelta]] EField[k[ 1, \[Omega] + \[CapitalDelta]], \[Omega] + \[CapitalDelta]], {\ \[CapitalDelta], \(-3\) \[Pi]/2, 3 \[Pi]/2, \[Pi]}]\)], "Input"], Cell[BoxData[ \(\[ExponentialE]\^\(\(-\(\(9\ \[Pi]\^2\)\/\(8\ \[CapitalGamma]p\^2\)\)\) \ + \[ImaginaryI]\ \((\(-t\)\ \((\(-\(\(3\ \[Pi]\)\/2\)\) + \[Omega])\) + \(z\ \ \((\(-\(\(3\ \[Pi]\)\/2\)\) + \[Omega])\)\)\/c)\)\)\/\(\@\(2\ \[Pi]\)\ \ \[CapitalGamma]p\) + \[ExponentialE]\^\(\(-\(\[Pi]\^2\/\(8\ \ \[CapitalGamma]p\^2\)\)\) + \[ImaginaryI]\ \((\(-t\)\ \((\(-\(\[Pi]\/2\)\) + \ \[Omega])\) + \(z\ \((\(-\(\[Pi]\/2\)\) + \[Omega])\)\)\/c)\)\)\/\(\@\(2\ \ \[Pi]\)\ \[CapitalGamma]p\) + \[ExponentialE]\^\(\(-\(\[Pi]\^2\/\(8\ \ \[CapitalGamma]p\^2\)\)\) + \[ImaginaryI]\ \((\(-t\)\ \((\[Pi]\/2 + \[Omega])\ \) + \(z\ \((\[Pi]\/2 + \[Omega])\)\)\/c)\)\)\/\(\@\(2\ \[Pi]\)\ \ \[CapitalGamma]p\) + \[ExponentialE]\^\(\(-\(\(9\ \[Pi]\^2\)\/\(8\ \ \[CapitalGamma]p\^2\)\)\) + \[ImaginaryI]\ \((\(-t\)\ \((\(3\ \[Pi]\)\/2 + \ \[Omega])\) + \(z\ \((\(3\ \[Pi]\)\/2 + \[Omega])\)\)\/c)\)\)\/\(\@\(2\ \[Pi]\ \)\ \[CapitalGamma]p\)\)], "Output"] }, Open ]], Cell[BoxData[ \(\(Table[ Plot[Evaluate[ Re[FourFreq /. {c \[Rule] 1, \[Omega] \[Rule] 20 \[Pi], \[CapitalGamma]p \[Rule] \[Pi]}]], {z, \ \(-12\)/\((2 \[Pi])\), 12/\((2 \[Pi])\)}, PlotRange \[Rule] {{\(-12\)/\((2 \[Pi])\), 12/\((2 \[Pi])\)}, {\(- .5\), .5}}, PlotPoints \[Rule] 70], {t, 0, 6/3 - 1/9, 1/9}];\)\)], "Input"], Cell[TextData[{ "Superposition of four monochromatic waves \n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"], Cell["\<\ Note that we see periodic \"pulse trains\" because the frequencies that we \ use have small least common multiples. We could vary the frequencies to \ remove this periodic behavior, at the expense of destroying the smooth \ \"pulse-like\" appearance. Or, by adding more and more frequencies, we can \ obtain nonperiodic pulses with smooth shape, as shown in the next section.\ \>", "Text"], Cell["\<\ Many short-pulse lasers operate in the so-called mode-locked regime, emitting \ a continuous train of pulses separated by a fixed time interval, the round \ trip time of light in the laser's cavity. The spectrum of such a laser \ corresponds to a \"frequency comb,\" a series of sharp uniformly spaced peaks \ extending over a spectal range corresponding to the reciprocal of the \ duration of an individual pulse, so the shorter the pulse, the broader the \ \"comb.\"\ \>", "Commentary"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Continuous frequency - Gaussian pulse", "Subsection"], Cell["\<\ A real pulse contains every frequency within its bandwidth. Integration over \ all frequencies yields the formula\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(GaussianPulse[z_, t_] = \((Integrate[ Envelope[\[CapitalDelta]] EField[k[ 1, \[Omega]0 + \[CapitalDelta]], \[Omega]0 + \ \[CapitalDelta]], {\[CapitalDelta], \(-\[Infinity]\), \[Infinity]}, Assumptions \[Rule] \[CapitalGamma]p > 0] // FullSimplify)\)\)], "Input", InitializationCell->True], Cell[BoxData[ \(\[ExponentialE]\^\(-\(\(\((c\ t - z)\)\ \((\((c\ t - z)\)\ \ \[CapitalGamma]p\^2 + 2\ \[ImaginaryI]\ c\ \[Omega]0)\)\)\/\(2\ \ c\^2\)\)\)\)], "Output"] }, Open ]], Cell["The envelope is, not surprisingly, a Gaussian:", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(Abs[GaussianPulse[z, t]] // ComplexExpand\) // Simplify\) // PowerExpand\)], "Input"], Cell[BoxData[ \(\[ExponentialE]\^\(-\(\(\((\(-c\)\ t + z)\)\^2\ \[CapitalGamma]p\^2\)\/\ \(2\ c\^2\)\)\)\)], "Output"] }, Open ]], Cell[BoxData[ \(\(Table[ Plot[Evaluate[ Re[GaussianPulse[z, t] /. {c \[Rule] 1, \[CapitalGamma]p \[Rule] \[Pi], \[Omega]0 \[Rule] 20 \[Pi]}]], {z, \(-12\)/\((2 \[Pi])\), 12/\((2 \[Pi])\)}, PlotRange \[Rule] {{\(-12\)/\((2 \[Pi])\), 12/\((2 \[Pi])\)}, {\(-1.5\), 1.5}}, PlotPoints \[Rule] 70], {t, 0, \[Pi], \[Pi]/16}];\)\)], "Input"], Cell[TextData[{ "A Gaussian pulse\n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Group velocity \[LongDash] general transparent medium\ \>", "Section"], Cell[" Consider two waves with slightly different frequencies", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Table[\[Omega]\_q = \[Omega]c + \(\((\(-1\))\)\^q\) \[Delta]/2, {q, 1, 2}] // TableForm\)], "Input", InitializationCell->True], Cell[BoxData[ InterpretationBox[GridBox[{ {\(\(-\(\[Delta]\/2\)\) + \[Omega]c\)}, {\(\[Delta]\/2 + \[Omega]c\)} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], TableForm[ { Plus[ Times[ Rational[ -1, 2], \[Delta]], \[Omega]c], Plus[ Times[ Rational[ 1, 2], \[Delta]], \[Omega]c]}]]], "Output"] }, Open ]], Cell[TextData[{ "Since the difference ", Cell[BoxData[ \(TraditionalForm\`\[Delta]\)]], " between the freqencies ", Cell[BoxData[ \(TraditionalForm\`\[Omega]\_1\)]], " and ", Cell[BoxData[ \(TraditionalForm\`\[Omega]\_2\)]], " is small, we can approximate the frequency dependence of ", Cell[BoxData[ \(TraditionalForm\`k\)]], " by a series about the central frequency ", Cell[BoxData[ \(TraditionalForm\`\[Omega]c\)]], ":" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(kSeries[\[Omega]_] = Series[k[\[Omega]], {\[Omega], \[Omega]c, 1}]\)], "Input", InitializationCell->True], Cell[BoxData[ InterpretationBox[ RowBox[{\(k[\[Omega]c]\), "+", RowBox[{ RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}], " ", \((\[Omega] - \[Omega]c)\)}], "+", InterpretationBox[\(O[\[Omega] - \[Omega]c]\^2\), SeriesData[ \[Omega], \[Omega]c, {}, 0, 2, 1]]}], SeriesData[ \[Omega], \[Omega]c, { k[ \[Omega]c], Derivative[ 1][ k][ \[Omega]c]}, 0, 2, 1]]], "Output"] }, Open ]], Cell[TextData[{ "Here ", Cell[BoxData[ FormBox[ RowBox[{\(\(k'\)[\[Omega]c]\), "=", RowBox[{\(\(d\ k\)\/\(d\ \[Omega]\)\), SubscriptBox[ AdjustmentBox["|", BoxMargins->{{-0.333333, 0.333333}, {0, 0}}], "\[Omega]c"]}]}], TraditionalForm]]], ". The wavenumbers of the two waves are given in this approximation by" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Table[k\_q = kSeries[\[Omega]\_q] // Normal, {q, 1, 2}] // TableForm\)], "Input", InitializationCell->True], Cell[BoxData[ InterpretationBox[GridBox[{ { RowBox[{\(k[\[Omega]c]\), "-", RowBox[{\(1\/2\), " ", "\[Delta]", " ", RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}]}]}]}, { RowBox[{\(k[\[Omega]c]\), "+", RowBox[{\(1\/2\), " ", "\[Delta]", " ", RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}]}]}]} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], TableForm[ { Plus[ k[ \[Omega]c], Times[ Rational[ -1, 2], \[Delta], Derivative[ 1][ k][ \[Omega]c]]], Plus[ k[ \[Omega]c], Times[ Rational[ 1, 2], \[Delta], Derivative[ 1][ k][ \[Omega]c]]]}]]], "Output"] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{"Evaluate", "[", RowBox[{\({Re[EField[k\_1, \[Omega]\_1]], Re[EField[k\_2, \[Omega]\_2]]}\), "/.", RowBox[{"{", RowBox[{\(\[Omega]c \[Rule] 6\), ",", \(k[\[Omega]c] \[Rule] 6\), ",", " ", RowBox[{ RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}], "\[Rule]", "1"}], ",", \(\[Delta] \[Rule] 1\)}], "}"}]}], "]"}], ",", \({z, 0, 2 \[Pi]}\), ",", \(PlotStyle \[Rule] {RGBColor[0, 0, 0], RGBColor[1, 0, 0]}\), ",", \(PlotRange \[Rule] {{0, 2 \[Pi]}, {\(-2.5\), 2.5}}\), ",", \(PlotPoints \[Rule] 30\)}], "]"}], ",", \({t, 0, 2 \[Pi] \((17/20)\), 2 \[Pi]/30}\)}], "]"}], ";"}]], "Input"], Cell[TextData[{ "Two monochromatic waves with different frequencies and wavenumbers\n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"], Cell["\<\ We see that the two waves are \"in phase\" at some regions of space and time, \ and are out of phase in others. Now let us look at the field corresponding to \ the sum of the two waves given by\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(EField[k\_1, \[Omega]\_1] + EField[k\_2, \[Omega]\_2]\)], "Input", InitializationCell->True], Cell[BoxData[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", RowBox[{"(", RowBox[{\(\(-t\)\ \((\(-\(\[Delta]\/2\)\) + \[Omega]c)\)\), "+", RowBox[{"z", " ", RowBox[{"(", RowBox[{\(k[\[Omega]c]\), "-", RowBox[{\(1\/2\), " ", "\[Delta]", " ", RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}]}]}], ")"}]}]}], ")"}]}]], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", RowBox[{"(", RowBox[{\(\(-t\)\ \((\[Delta]\/2 + \[Omega]c)\)\), "+", RowBox[{"z", " ", RowBox[{"(", RowBox[{\(k[\[Omega]c]\), "+", RowBox[{\(1\/2\), " ", "\[Delta]", " ", RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}]}]}], ")"}]}]}], ")"}]}]]}]], "Output"] }, Open ]], Cell["Simplifying:", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(SumField = % // FullSimplify\)], "Input", InitializationCell->True], Cell[BoxData[ RowBox[{ "2", " ", \(\[ExponentialE]\^\(\(-\[ImaginaryI]\)\ \((t\ \[Omega]c - z\ k[\[Omega]c])\)\)\), " ", RowBox[{"Cos", "[", RowBox[{\(1\/2\), " ", "\[Delta]", " ", RowBox[{"(", RowBox[{"t", "-", RowBox[{"z", " ", RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}]}]}], ")"}]}], "]"}]}]], "Output"] }, Open ]], Cell["The real field is given by", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Re[SumField] // ComplexExpand\) // PowerExpand\)], "Input"], Cell[BoxData[ RowBox[{"2", " ", \(Cos[t\ \[Omega]c - z\ k[\[Omega]c]]\), " ", RowBox[{"Cos", "[", RowBox[{\(1\/2\), " ", "\[Delta]", " ", RowBox[{"(", RowBox[{"t", "-", RowBox[{"z", " ", RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}]}]}], ")"}]}], "]"}]}]], "Output"] }, Open ]], Cell["and the envelope is given by", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(SumFieldEnv = \(Abs[SumField] // ComplexExpand\) // PowerExpand\)], "Input", InitializationCell->True], Cell[BoxData[ RowBox[{"2", " ", RowBox[{"Cos", "[", RowBox[{\(1\/2\), " ", "\[Delta]", " ", RowBox[{"(", RowBox[{"t", "-", RowBox[{"z", " ", RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}]}]}], ")"}]}], "]"}]}]], "Output"] }, Open ]], Cell[TextData[{ "We assume here that ", Cell[BoxData[ \(TraditionalForm\`k\)]], " is real, i.e. that the medium is transparent." }], "Text"], Cell[BoxData[ RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{"Evaluate", "[", RowBox[{\({Re[SumField], SumFieldEnv, \(-SumFieldEnv\)}\), "/.", RowBox[{"{", RowBox[{\(\[Omega]c \[Rule] 6\), ",", \(k[\[Omega]c] \[Rule] 4\), ",", " ", RowBox[{ RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}], "\[Rule]", "1"}], ",", \(\[Delta] \[Rule] 1. \)}], "}"}]}], "]"}], ",", \({z, 0, 8 \[Pi]}\), ",", \(PlotRange \[Rule] {{0, 8 \[Pi]}, {\(-2.5\), 2.5}}\), ",", \(PlotPoints \[Rule] 50\), ",", \(PlotStyle \[Rule] {RGBColor[0, 0, 0], RGBColor[1, 0, 0], RGBColor[1, 0, 0]}\)}], "]"}], ",", \({t, 0, 4 \[Pi] \((19/20)\), 4 \[Pi]/20}\)}], "]"}], ";"}]], "Input"], Cell[TextData[{ "Two superposed monochromatic waves in a medium, plotted with their \ envelope. Note that since the group velocity does not equal the phase \ velocity, a phase difference between the light field and the envelope is \ acumulated in time.\n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"], Cell[TextData[{ "To determine the velocity of the pulse, we consider the motion of a given \ point on the envelope, for example, a zero crossing: ", Cell[BoxData[ \(TraditionalForm\`z[t]\)]], " such that the ", Cell[BoxData[ \(TraditionalForm\`\(\(|\)\(E[t]\)\(|\)\) = 0\)]], ":" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(zConstantEnv = \(Solve[SumFieldEnv \[Equal] 0, z]\)\[LeftDoubleBracket]1, 1\[RightDoubleBracket]\)], "Input", InitializationCell->True], Cell[BoxData[ RowBox[{"z", "\[Rule]", FractionBox[\(\(-\[Pi]\) + t\ \[Delta]\), RowBox[{"\[Delta]", " ", RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}]}]]}]], "Output"] }, Open ]], Cell[TextData[{ "The derivative of ", Cell[BoxData[ \(TraditionalForm\`z\)]], " with respect to time is the group velocity" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Vg[k_] = D[z /. zConstantEnv, t]\)], "Input", InitializationCell->True], Cell[BoxData[ FractionBox["1", RowBox[{ SuperscriptBox["k", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}]]], "Output"] }, Open ]], Cell[TextData[{ "as a function of ", Cell[BoxData[ \(TraditionalForm\`\[Omega]\)]], ", ", Cell[BoxData[ \(TraditionalForm\`k\)]], " is given by " }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(kOf\[Omega][\[Omega]_] = k[n[\[Omega]], \[Omega]]\)], "Input", InitializationCell->True], Cell[BoxData[ \(\(\[Omega]\ n[\[Omega]]\)\/c\)], "Output"] }, Open ]], Cell["\<\ so the group velocity is given in terms of the index of refraction by\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(VgOfn[n_] = Vg[kOf\[Omega]] // Simplify\)], "Input", InitializationCell->True], Cell[BoxData[ FractionBox["c", RowBox[{\(n[\[Omega]c]\), "+", RowBox[{"\[Omega]c", " ", RowBox[{ SuperscriptBox["n", "\[Prime]", MultilineFunction->None], "[", "\[Omega]c", "]"}]}]}]]], "Output"] }, Open ]], Cell["\<\ This derivation is for the case of a transparent medium. As is pointed out in \ most optics texts, one has to be careful when applying the concept of group \ velocity to a medium with strong absorption, because pulse reshaping will \ occur as the pulse travels through the medium. We will look at this more \ closely next.\ \>", "Text"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Pulse propagation in linear atomic media\ \>", "Section"], Cell[CellGroupData[{ Cell["Nontransparent media", "Subsection"], Cell["\<\ We would like to determine when our formula for group velocity in transparent \ media will apply to nontransparent media. For a Gaussian pulse, the group \ velocity refers to the velocity of the peak of the pulse. Although inside an \ absorptive medium the pulse shape is distorted spatially, we can still \ measure the group velocity by comparing the times at which the pulse reaches \ a maximum at different points inside the medium. Our formula for the group \ velocity should hold as long as there is not too much distortion of the pulse \ shape, which would render the concept of a \"group\" velocity (the velocity \ at which the entire pulse travels) meaningless. Distortion is caused when \ resonant spectral components are absorbed much more strongly than \ off-resonant ones.\ \>", "Text"], Cell[TextData[{ "Absorption is governed by the imaginary part of ", Cell[BoxData[ \(TraditionalForm\`k\)]], ":" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(ImkMedium[\[Omega]_] = \(\[ImaginaryI]\ Im[kMedium[\[Omega]]] // ComplexExpand\) // FullSimplify\)], "Input"], Cell[BoxData[ \(\(\[ImaginaryI]\ \[Alpha]\ \[CapitalGamma]\^2\ \[Omega]\)\/\(2\ \((\ \[CapitalGamma]\^2 + 4\ \((\[Omega] - \[Omega]0)\)\^2)\)\ \[Omega]0\)\)], \ "Output"] }, Open ]], Cell["\<\ Let's consider a Gaussian input pulse. The spectrum of the input intensity is \ given by\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Abs[Envelope[\[CapitalDelta]]]\^2 // ComplexExpand\)], "Input"], Cell[BoxData[ \(\[ExponentialE]\^\(-\(\[CapitalDelta]\^2\/\[CapitalGamma]p\^2\)\)\/\(2\ \ \[Pi]\ \[CapitalGamma]p\^2\)\)], "Output"] }, Open ]], Cell["\<\ The spectrum of the output intensity is given by\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Transmission[\[CapitalDelta]_] = \(Abs[Envelope[\[CapitalDelta]] \ EField[ImkMedium[\[CapitalDelta] + \[Omega]0], \[Omega]] /. \[CapitalDelta] + \ \[Omega]0 \[Rule] \[Omega]0]\^2 // ComplexExpand\) // PowerExpand\)], "Input"], Cell[BoxData[ \(\[ExponentialE]\^\(\(-\(\[CapitalDelta]\^2\/\[CapitalGamma]p\^2\)\) - \ \(z\ \[Alpha]\ \[CapitalGamma]\^2\)\/\(\[CapitalGamma]\^2 + 4\ \ \[CapitalDelta]\^2\)\)\/\(2\ \[Pi]\ \[CapitalGamma]p\^2\)\)], "Output"] }, Open ]], Cell[TextData[{ "We assume here that ", Cell[BoxData[ \(TraditionalForm\`\[Omega] \[TildeEqual] \[Omega]0\)]], "." }], "Text"], Cell[TextData[{ "A Gaussian spectrum has a maxium only in the center. If we require this to \ be true for the output spectrum, we can obtain a condition for avoiding pulse \ reshaping. Setting the derivative with respect to ", Cell[BoxData[ \(TraditionalForm\`\[CapitalDelta]\)]], " equal to zero, we can find the extrema:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Extrema = Solve[D[Transmission[\[CapitalDelta]], \[CapitalDelta]] \[Equal] 0, \[CapitalDelta]]\)], "Input"], Cell[BoxData[ \({{\[CapitalDelta] \[Rule] 0}, {\[CapitalDelta] \[Rule] \(-\@\(\(-\(\[CapitalGamma]\^2\/4\)\) \ - 1\/2\ \@z\ \@\[Alpha]\ \[CapitalGamma]\ \[CapitalGamma]p\)\)}, {\ \[CapitalDelta] \[Rule] \@\(\(-\(\[CapitalGamma]\^2\/4\)\) - 1\/2\ \@z\ \@\ \[Alpha]\ \[CapitalGamma]\ \[CapitalGamma]p\)}, {\[CapitalDelta] \[Rule] \ \(-\@\(\(-\(\[CapitalGamma]\^2\/4\)\) + 1\/2\ \@z\ \@\[Alpha]\ \ \[CapitalGamma]\ \[CapitalGamma]p\)\)}, {\[CapitalDelta] \[Rule] \@\(\(-\(\ \[CapitalGamma]\^2\/4\)\) + 1\/2\ \@z\ \@\[Alpha]\ \[CapitalGamma]\ \ \[CapitalGamma]p\)}}\)], "Output"] }, Open ]], Cell[TextData[{ "The first solution is zero, and the second and third are imaginary. To \ find the condition for excluding the last two (i.e., requiring them to be \ imaginary), we solve for ", Cell[BoxData[ \(TraditionalForm\`\[CapitalDelta] = 0\)]], " in the fourth solution:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Solve[\[CapitalDelta] \[Equal] 0 /. Extrema\[LeftDoubleBracket]4\[RightDoubleBracket], \[Alpha]]\)], \ "Input"], Cell[BoxData[ \({{\[Alpha] \[Rule] \[CapitalGamma]\^2\/\(4\ z\ \ \[CapitalGamma]p\^2\)}}\)], "Output"] }, Open ]], Cell[TextData[{ "So we need ", Cell[BoxData[ \(TraditionalForm\`\[Alpha]\ z << \ \(\[CapitalGamma]\^2\)/\[CapitalGamma]p\^2\)]], " to avoid significant pulse reshaping. If this condition is satisfied, the \ formula we derived for group velocity in a transparent medium will remain \ approximately correct." }], "Text"], Cell[BoxData[ \(\(IndexPlot[{c \[Rule] 1, \[Alpha] \[Rule] 1, \[CapitalGamma] \[Rule] 1, \[Omega]0 \[Rule] 2 \[Pi], \[CapitalGamma]p \[Rule] .1}, PlotEnvelope \[Rule] True];\)\)], "InputOutdent"], Cell[TextData[StyleBox["The real and imaginary parts of the index of \ refraction in a medium as a function of frequency, plotted together with the \ Gaussian pulse spectrum of width \[CapitalGamma]p.", FontColor->GrayLevel[0]]], "CaptionOutdent"] }, Open ]], Cell[CellGroupData[{ Cell["Group velocity", "Subsection"], Cell["\<\ We now apply the expression for group velocity that we derived for \ transparent media to a linear atomic medium near resonance. Substituting our \ formula for the real part of the index of refraction in a atomic medium into \ the expression for group velocity, we obtain for the group velocity in a \ linear atomic medium near resonance\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(VgOfn[RenMedium]\)], "Input", InitializationCell->True], Cell[BoxData[ \(c\/\(1 - \(\[Xi]\ \[Omega]c\)\/\[Omega]0 + \[Xi]\ \((1 - \[Omega]c\/\ \[Omega]0)\)\)\)], "Output"] }, Open ]], Cell[TextData[{ "Using our definition of ", Cell[BoxData[ \(TraditionalForm\`\[Xi]\)]], " and assuming that we are directly on resonance gives" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(VgMedium = \(% /. \[Xi] \[Rule] \[Xi][\[Omega]c]\) /. \[Omega]c \[Rule] \ \[Omega]0\)], "Input", InitializationCell->True], Cell[BoxData[ \(c\/\(1 - \(c\ \[Alpha]\)\/\[CapitalGamma]\)\)], "Output"] }, Open ]], Cell[TextData[{ "We can see that on resonance, when ", Cell[BoxData[ \(TraditionalForm\`\[Alpha]\ c < \[CapitalGamma]\)]], " the group velocity is greater than ", Cell[BoxData[ \(TraditionalForm\`c\)]], " and when ", Cell[BoxData[ \(TraditionalForm\`\[Alpha]\ c > \[CapitalGamma]\)]], " the group velocity is negative. While at first glance, these results may \ seem unphysical (especially the second, which implies that the pulse moves \ backwards in the medium!) they actually have simple, consistent \ interpretations. We illustrate this below." }], "Text"], Cell["\<\ It is important to note that there is a change in the spatial width of the \ pulse upon propagation through an interface between vacuum and the medium -- \ the pulse compresses or expands depending on the value of the group velocity. \ Clearly, while the pulse is going through the interface it is severely \ reshaped; however, we may still discuss the behaviour of the pulse when most \ of its spatial extent is contained within the medium or one of the two vacuum \ regions -- before and after the medium.\ \>", "Commentary"] }, Open ]], Cell[CellGroupData[{ Cell["Setup for plotting", "Subsection"], Cell[TextData[{ "There is a fast-oscillating phase ", Cell[BoxData[ \(\[ExponentialE]\^\(\(-\[ImaginaryI]\)\ \[Omega]0 \((t - z)\)\)\)]], " in the electric field that we don't need to include if we are plotting \ the envelope. We'll divide it out of the expression for the E-field:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Map[Factor, ExpandAll[ EFieldMedium[z, t, \[Omega]0 + \[CapitalDelta]] \[ExponentialE]\^\(\[ImaginaryI]\ \ \[Omega]0 \((t - z/c)\)\)], {2}]\)], "Input", InitializationCell->True], Cell[BoxData[ \(\[ExponentialE]\^\(\(-\(\(z\ \[Alpha]\ \[CapitalGamma]\)\/\(2\ \((\ \[CapitalGamma] - 2\ \[ImaginaryI]\ \[CapitalDelta])\)\)\)\) - \[ImaginaryI]\ \ t\ \[CapitalDelta] + \(\[ImaginaryI]\ z\ \[CapitalDelta]\)\/c - \(z\ \[Alpha]\ \ \[CapitalGamma]\ \[CapitalDelta]\)\/\(2\ \((\[CapitalGamma] - 2\ \ \[ImaginaryI]\ \[CapitalDelta])\)\ \[Omega]0\)\)\)], "Output"] }, Open ]], Cell[TextData[{ "Assuming that ", Cell[BoxData[ \(TraditionalForm\`\[CapitalGamma] << \[Omega]0\)]], ", we can set the term containing ", Cell[BoxData[ \(TraditionalForm\`\[CapitalGamma]/\[Omega]0\)]], " to zero:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(EFieldEnv[z_, t_, \[CapitalDelta]_] = % /. \[CapitalGamma]/\[Omega]0 \[Rule] 0\)], "Input", InitializationCell->True], Cell[BoxData[ \(\[ExponentialE]\^\(\(-\(\(z\ \[Alpha]\ \[CapitalGamma]\)\/\(2\ \((\ \[CapitalGamma] - 2\ \[ImaginaryI]\ \[CapitalDelta])\)\)\)\) - \[ImaginaryI]\ \ t\ \[CapitalDelta] + \(\[ImaginaryI]\ z\ \[CapitalDelta]\)\/c\)\)], "Output"] }, Open ]], Cell["This function plots a pulse traversing a medium", "Text"], Cell[BoxData[ \(PlotField[FieldFunction_, {z1_, z2_}, EnvelopePlotOptions_] := \[IndentingNewLine]Plot[ Abs[FieldFunction], {z, z1, z2}, Evaluate[ Join[EnvelopePlotOptions, {DisplayFunction \[Rule] Identity, PlotDivision \[Rule] 2, PlotPoints \[Rule] 10}]]]\)], "Input", InitializationCell->True], Cell[BoxData[ \(NPlotField[FieldFunction_, {z1_, z2_}, IntRange_, EnvelopePlotOptions_] := \[IndentingNewLine]Plot[ Abs[NIntegrate[ FieldFunction, {\[CapitalDelta], \(-IntRange\), IntRange}, PrecisionGoal \[Rule] 3, WorkingPrecision \[Rule] 12, MaxRecursion \[Rule] 12, MinRecursion \[Rule] 2]], {z, z1, z2}, Evaluate[ Join[EnvelopePlotOptions, {DisplayFunction \[Rule] Identity, PlotDivision \[Rule] 2, PlotPoints \[Rule] 10}]]]\)], "Input", InitializationCell->True], Cell[BoxData[ RowBox[{\(PlotPacket[{z0_, z1_, z2_, z3_}, {tMin_, tMax_, tStep_}, Parameters_]\), ":=", RowBox[{"(", "\[IndentingNewLine]", StyleBox[\( (*\ Define\ styles\ for\ plotting\ *) \), FontColor->GrayLevel[0.500008]], "\[IndentingNewLine]", RowBox[{\(FieldOptions = {PlotStyle \[Rule] {AbsoluteThickness[2]}}\), ";", "\[IndentingNewLine]", \(EnvelopeOptions = {PlotStyle \[Rule] \ {AbsoluteThickness[2], RGBColor[1, 0, 0]}}\), ";", "\[IndentingNewLine]", \(GhostFieldOptions = {PlotStyle \[Rule] \ {AbsoluteThickness[1.5], AbsoluteDashing[{2, 5}]}}\), ";", "\[IndentingNewLine]", \(GhostEnvelopeOptions = {PlotStyle \[Rule] \ {AbsoluteThickness[1.5], RGBColor[1, 0, 0], AbsoluteDashing[{2, 5}]}}\), ";", "\[IndentingNewLine]", StyleBox[\( (*\ Define\ the\ field\ for\ each\ phase\ *) \), FontColor->GrayLevel[0.500008]], "\[IndentingNewLine]", \(NoMedField[z_, t_] = \(GaussianPulse[z, t] /. \[Omega]0 \[Rule] 0\) /. Parameters\), ";", "\[IndentingNewLine]", \(Region1Field[z_, t_, \[CapitalDelta]1_] = \(EFieldEnv[z, t, \[CapitalDelta]1] /. {\[Alpha] \[Rule] 0}\) /. Parameters\), ";", "\[IndentingNewLine]", \(Region2Field[z_, t_, \[CapitalDelta]1_] = Region1Field[z1, 0, \[CapitalDelta]1] EFieldEnv[z - z1, t, \[CapitalDelta]1] /. Parameters\), ";", "\[IndentingNewLine]", \(Region3Field[z_, t_, \[CapitalDelta]1_] = \(Region2Field[z2, 0, \[CapitalDelta]1] EFieldEnv[z - z2, t, \[CapitalDelta]1] /. {\[Alpha] \[Rule] 0}\) /. Parameters\), ";", "\[IndentingNewLine]", \(R3Mag = \(Region3Mag /. Parameters\) /. Region3Mag \[Rule] 0\), ";", "\[IndentingNewLine]", \(Env[\[CapitalDelta]1_] = Envelope[\[CapitalDelta]1] /. Parameters\), ";", "\[IndentingNewLine]", \(IntRange = 4\ \[CapitalGamma]p /. Parameters\), ";", "\[IndentingNewLine]", StyleBox[\( (*\ Draw\ plots\ for\ each\ region\ *) \), FontColor->GrayLevel[0.500008]], "\[IndentingNewLine]", \(Table[\[IndentingNewLine]Show[\ \[IndentingNewLine]PlotField[NoMedField[z, t], {z0, z1}, EnvelopeOptions], \[IndentingNewLine]PlotField[ NoMedField[z, t], {z1, z3}, GhostEnvelopeOptions], \[IndentingNewLine]NPlotField[ 10^\((R3Mag \((z - z1)\)/\((z2 - z1)\))\)\ Env[\[CapitalDelta]] Region2Field[z, t, \[CapitalDelta]], {z1, z2}, IntRange, EnvelopeOptions], \[IndentingNewLine]NPlotField[ 10\^R3Mag\ Env[\[CapitalDelta]] Region3Field[z, t, \[CapitalDelta]], {z2, z3}, IntRange, EnvelopeOptions], \[IndentingNewLine]DisplayFunction \[Rule] \ $DisplayFunction, Axes \[Rule] {True, False}, Epilog \[Rule] {AbsoluteThickness[3], Line[{{z1, \(-1.4\)}, {z1, 1.4}, {z2, 1.4}, {z2, \(-1.4\)}, {z1, \(-1.4\)}}], If[R3Mag \[NotEqual] 0, {Text[\\!\\("\<\[Times]10\>" \\^ \ R3Mag\\), {z2 + .2 \((z3 - z2)\), 1.3}]}, {}]}, PlotRange \[Rule] {{z0 - 0.00001, z3}, {0, 1.5}}\[IndentingNewLine]], \[IndentingNewLine]{t, tMin, tMax, tStep}]\), ";"}], "\[IndentingNewLine]", ")"}]}]], "InputOnly", InitializationCell->True] }, Open ]], Cell[CellGroupData[{ Cell["Plots of pulse propagation", "Subsection"], Cell[CellGroupData[{ Cell["Superluminal pulse", "Subsubsection"], Cell[TextData[{ "We start with the case in which ", Cell[BoxData[ \(TraditionalForm\`\[Alpha]\ c < \[CapitalGamma]\)]], ". Here ", Cell[BoxData[ \(TraditionalForm\`\[Alpha]\ c/\[CapitalGamma] = 0.5\)]], ", and the group velocity is " }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(VgMedium /. {\[Alpha] \[Rule] 1/\((2 c)\), \[CapitalGamma] \[Rule] 1}\)], "Input"], Cell[BoxData[ \(2\ c\)], "Output"] }, Open ]], Cell[TextData[{ "For all the plots below, we set the light frequency to resonance. We see \ that even though the pulse is partly absorbed (i.e. its amplitude at the \ output is smaller than that at the input), the output pulse has approximately \ the same shape as the input. The peak of the output pulse is advanced \ relative to where the peak would have been hand there been no medium (dashed \ line), i.e., the pulse appears to travel faster than ", Cell[BoxData[ \(TraditionalForm\`c\)]], " inside the medium. No photons travel faster than ", Cell[BoxData[ \(TraditionalForm\`c\)]], " here, rather, the trailing edge of the pulse is absorbed more strongly \ than the leading edge creating the ", "illusion", " of superluminal propagation. " }], "Text"], Cell[BoxData[ \(PlotPacket[{0. , 25. , 35. , 60. }, {\(-30. \), 90. , 3. }, {c \[Rule] 1, \[Alpha] \[Rule] .5, \[CapitalGamma]p \[Rule] .1, \ \[CapitalGamma] \[Rule] 1. }]\)], "Input"], Cell[TextData[{ StyleBox["A Gaussian pulse propagating through a linear medium.", FontColor->GrayLevel[0]], "\n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"] }, Open ]], Cell[CellGroupData[{ Cell["Negative group velocity", "Subsubsection"], Cell[TextData[{ "When ", Cell[BoxData[ \(TraditionalForm\`\[Alpha]\ c < \[CapitalGamma]\)]], ", the group velocity is negative:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(VgMedium /. {\[Alpha] \[Rule] 6/c, \[CapitalGamma] \[Rule] 2}\)], "Input"], Cell[BoxData[ \(\(-\(c\/2\)\)\)], "Output"] }, Open ]], Cell["\<\ This implies that the peak of the output pulse will exit the medium before \ the peak of the input pulse enters it:\ \>", "Text"], Cell[BoxData[ \(PlotPacket[{0. , 25. , 35. , 60. }, {\(-30. \), 90. , 3. }, {c \[Rule] 1, \[Alpha] \[Rule] 6. , \[CapitalGamma]p \[Rule] .1, \[CapitalGamma] \[Rule] 2. , Region3Mag \[Rule] 8}]\)], "Input"], Cell[TextData[{ StyleBox["A Gaussian pulse propagating through a linear medium.", FontColor->GrayLevel[0]], "\n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"], Cell[TextData[{ "The absorption is now so large that we have to scale the amplitude of the \ output pulse by ", Cell[BoxData[ \(TraditionalForm\`10\^8\)]], " in order to see what is going on. Inside the medium the scaling is ramped \ up exponentially so that the plot remains smooth." }], "Text"], Cell[TextData[{ "By adjusting the parameters, we can actually see the pulse moving \ backwards in the medium, although the ", Cell[BoxData[ \(TraditionalForm\`10\^34\)]], "scaling required to observe the output makes this situation rather \ unphysical. We can also see that the pulse is compressed inside the medium." }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(VgMedium /. {\[Alpha] \[Rule] 16/c, \[CapitalGamma] \[Rule] 4}\)], "Input"], Cell[BoxData[ \(\(-\(c\/3\)\)\)], "Output"] }, Open ]], Cell[BoxData[ \(PlotPacket[{0. , 25. , 35. , 60. }, {\(-30. \), 90. , 3. }, {c \[Rule] 1, \[Alpha] \[Rule] 16. , \[CapitalGamma]p \[Rule] .1, \[CapitalGamma] \[Rule] 4. , Region3Mag \[Rule] 34}]\)], "Input"], Cell[TextData[{ StyleBox["A Gaussian pulse propagating through a linear medium.", FontColor->GrayLevel[0]], "\n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"], Cell[TextData[{ "Actually, the spatial pulse-like shape in the medium is an artifact of our \ scaling. The electric field at any given time is, in fact, a monotonically \ decreasing function of ", Cell[BoxData[ \(TraditionalForm\`z\)]], ". However, if we look at the value of the field at a given point as a \ function of time, we can observe pulse-like behavior even inside the medium." }], "Text"], Cell["\<\ Negative group velocities can be observed in experiments in linear optics. A \ situation much closer to the common experimental setup than the above case is \ one in which the pulse advance is much smaller than the total pulse width. \ Here the pulse is much wider than in the previous examples, but we still see \ that the output pulse peaks at the output before the input pulse peaks at the \ input. \ \>", "Text"], Cell[BoxData[ \(PlotPacket[{0. , 25. , 35. , 60. }, {\(-3. \)\ 10\^10, 3. \ 10\^10, 3. \ 10\^9}, {c \[Rule] 1, \[Alpha] \[Rule] 1. , \[CapitalGamma]p \[Rule] 10\^\(-10\), \[CapitalGamma] \[Rule] 10\^\(-9\), Region3Mag \[Rule] 2}]\)], "Input"], Cell[TextData[{ StyleBox["A Gaussian pulse propagating through a linear medium.", FontColor->GrayLevel[0]], "\n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Non-Gaussian pulse", "Section 1"], Cell[TextData[{ "The pulse shape is not restricted to a simple Gaussian, as long as the \ pulse shape is analytic, the entire pulse spectrum is within the region of \ linear anomalous dispersion, and ", Cell[BoxData[ \(TraditionalForm\`\[Alpha]\ z << \ \(\[CapitalGamma]\^2\)/\[CapitalGamma]p\^2\)]], ", as derived above. Although complex pulse shapes can pass through a \ medium with group velocity greater than ", Cell[BoxData[ \(TraditionalForm\`c\)]], ", information does not travel faster than the speed of light, since the \ requirement that the pulse shape is analytic means that all the information \ in the pulse is contained in the front \"tail\" of the pulse. The fact that \ the peak of the pulse moves faster than ", Cell[BoxData[ \(TraditionalForm\`c\)]], " does not mean that the information in the pulse does." }], "Text"], Cell["\<\ We illustrate the propagation of a non-Gaussian pulse with a double-peaked \ pulse (the sum of two Gaussians).\ \>", "Text"], Cell[CellGroupData[{ Cell["Setup for plotting", "Subsection"], Cell[BoxData[ RowBox[{\(Plot2Packets[{z0_, z1_, z2_, z3_}, {tMin_, tMax_, tStep_}, Parameters_]\), ":=", RowBox[{"(", "\[IndentingNewLine]", StyleBox[\( (*\ Define\ styles\ for\ plotting\ *) \), FontColor->GrayLevel[0.500008]], "\[IndentingNewLine]", RowBox[{\(FieldOptions = {PlotStyle \[Rule] {AbsoluteThickness[2]}}\), ";", "\[IndentingNewLine]", \(EnvelopeOptions = {PlotStyle \[Rule] \ {AbsoluteThickness[2], RGBColor[1, 0, 0]}}\), ";", "\[IndentingNewLine]", \(GhostFieldOptions = {PlotStyle \[Rule] \ {AbsoluteThickness[1.5], AbsoluteDashing[{2, 5}]}}\), ";", "\[IndentingNewLine]", \(GhostEnvelopeOptions = {PlotStyle \[Rule] \ {AbsoluteThickness[1.5], RGBColor[1, 0, 0], AbsoluteDashing[{2, 5}]}}\), ";", "\[IndentingNewLine]", StyleBox[\( (*\ Define\ the\ field\ for\ each\ phase\ *) \), FontColor->GrayLevel[0.500008]], "\[IndentingNewLine]", \(NoMedField[z_, t_] = \(GaussianPulse[z, t] + GaussianPulse[z, t + ts] /. \[Omega]0 \[Rule] 0\) /. Parameters\), ";", "\[IndentingNewLine]", \(Region1Field[z_, t_, \[CapitalDelta]1_] = \(EFieldEnv[z, t, \[CapitalDelta]1] + EFieldEnv[z, t + ts, \[CapitalDelta]1] /. {\[Alpha] \[Rule] 0}\) /. Parameters\), ";", "\[IndentingNewLine]", \(Region2Field[z_, t_, \[CapitalDelta]1_] = Region1Field[z1, 0, \[CapitalDelta]1] EFieldEnv[z - z1, t, \[CapitalDelta]1] /. Parameters\), ";", "\[IndentingNewLine]", \(Region3Field[z_, t_, \[CapitalDelta]1_] = \(Region2Field[z2, 0, \[CapitalDelta]1] EFieldEnv[z - z2, t, \[CapitalDelta]1] /. {\[Alpha] \[Rule] 0}\) /. Parameters\), ";", "\[IndentingNewLine]", \(R3Mag = \(Region3Mag /. Parameters\) /. Region3Mag \[Rule] 0\), ";", "\[IndentingNewLine]", \(Env[\[CapitalDelta]1_] = Envelope[\[CapitalDelta]1] /. Parameters\), ";", "\[IndentingNewLine]", \(IntRange = 4\ \[CapitalGamma]p /. Parameters\), ";", "\[IndentingNewLine]", StyleBox[\( (*\ Draw\ plots\ for\ each\ region\ *) \), FontColor->GrayLevel[0.500008]], "\[IndentingNewLine]", \(Table[\[IndentingNewLine]Show[\ \[IndentingNewLine]PlotField[NoMedField[z, t], {z0, z1}, EnvelopeOptions], \[IndentingNewLine]PlotField[ NoMedField[z, t], {z1, z3}, GhostEnvelopeOptions], \[IndentingNewLine]NPlotField[ 10^\((R3Mag \((z - z1)\)/\((z2 - z1)\))\)\ Env[\[CapitalDelta]] Region2Field[z, t, \[CapitalDelta]], {z1, z2}, IntRange, EnvelopeOptions], \[IndentingNewLine]NPlotField[ 10\^R3Mag\ Env[\[CapitalDelta]] Region3Field[z, t, \[CapitalDelta]], {z2, z3}, IntRange, EnvelopeOptions], \[IndentingNewLine]DisplayFunction \[Rule] \ $DisplayFunction, Axes \[Rule] {True, False}, Epilog \[Rule] {AbsoluteThickness[3], Line[{{z1, \(-1.4\)}, {z1, 1.4}, {z2, 1.4}, {z2, \(-1.4\)}, {z1, \(-1.4\)}}], If[R3Mag \[NotEqual] 0, {Text[\\!\\("\<\[Times]10\>" \\^ \ R3Mag\\), {z2 + .2 \((z3 - z2)\), 1.3}]}, {}]}, PlotRange \[Rule] {{z0 - 0.00001, z3}, {0, 1.5}}\[IndentingNewLine]], \[IndentingNewLine]{t, tMin, tMax, tStep}]\), ";"}], "\[IndentingNewLine]", ")"}]}]], "InputOnly", InitializationCell->True] }, Open ]], Cell[CellGroupData[{ Cell["Plots", "Subsection"], Cell["\<\ Here we plot a double-peaked pulse propagating through a linear medium. The \ entire pulse propagates at the same group velocity, and exits the medium \ advanced relative to where it would have been, had there been no medium.\ \>", "Text"], Cell[BoxData[ \(Plot2Packets[{0. , 25. , 35. , 60. }, {\(-30. \), 90. , 3. }, {c \[Rule] 1, \[Alpha] \[Rule] .5, \[CapitalGamma]p \[Rule] .1, \ \[CapitalGamma] \[Rule] 1. , ts \[Rule] 25. , Region3Mag \[Rule] 1}]\)], "Input"], Cell[TextData[{ StyleBox["A double-peaked pulse propagating through a linear medium.", FontColor->GrayLevel[0]], "\n", StyleBox["(double-click plot to view as a function of time)", FontColor->RGBColor[0.972671, 0.308598, 0.273442]] }], "Caption"] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Conclusion", "Subtitle", CellTags->"conclusions"], Cell["We hope that this tutorial has demonstrated that:", "Text"], Cell["\<\ All the \"unusual\" light propagation effects in this system follow directly \ from the form of the complex refractive index and elementary Fourier \ analysis.\ \>", "Text", CellDingbat->"\[FivePointedStar]"], Cell[TextData[{ "Even though it may appear as if the peak of the light pulse is moving \ faster than ", Cell[BoxData[ \(TraditionalForm\`c\)]], " (under certain conditions), nothing is actually traveling that fast. The \ \"superluminal\" effects arise because the back of the pulse is absorbed more \ than the front, shifting the apparent peak of the pulse forward." }], "Text", CellDingbat->"\[FivePointedStar]"], Cell[TextData[{ "The usual formula for the group velocity \n\t", Cell[BoxData[ FormBox[GridBox[{ {\(v\_g\), "=", RowBox[{"Re", FractionBox[ StyleBox[ RowBox[{ StyleBox["d", FontSlant->"Italic"], "\[Omega]"}]], \(d\ k\)]}]}, {" ", "=", \(c\/\(Re\ n\ + \ \[Omega]\ Re\ \(d\ n\)\/\(d\ \[Omega]\)\ \)\)} }], TraditionalForm]]], ",\nwhere ", Cell[BoxData[ \(TraditionalForm\`\[Omega]\)]], " is the light frequency, ", Cell[BoxData[ \(TraditionalForm\`k\)]], " is the complex wavenumber, and ", Cell[BoxData[ \(TraditionalForm\`n\)]], " is the complex refractive index, remains valid even in the cases where \ such velocity becomes infinite or negative, with all the seemingly strange \ consequencies to which this leads. " }], "Text", CellDingbat->"\[FivePointedStar]"], Cell[CellGroupData[{ Cell["Acknowledgement and disclaimer", "Section 1"], Cell[TextData[{ StyleBox["We would like to acknowldge helpful discussions with M. \ Zolotorev. Creation of this tutorial was supported by a UC Berkeley Classroom \ Technologies Grant (2002). The material is based in part upon work supported \ by the ", FontSize->12], ButtonBox["National Science Foundation", ButtonData:>{ URL[ "http://www.nsf.gov/"], None}, ButtonStyle->"Hyperlink"], StyleBox[" (NSF) under Grant No. PHY-9733479 and by the ", FontSize->12], ButtonBox["Office of Naval Research", ButtonData:>{ URL[ "http://www.onr.navy.mil/"], None}, ButtonStyle->"Hyperlink"], StyleBox[" (ONR) under Grant No. N00014-97-1-0214. Any opinions, findings, \ conclusions or recomendations expressed in this material are those of the \ authors and do not necessarily reflect the views of the NSF or the ONR. \n", FontSize->12] }], "Text"] }, Open ]], Cell[CellGroupData[{ Cell["References", "Section 1"], Cell[TextData[{ StyleBox["R.W. Boyd and D.J. Gauthier, in ", FontFamily->"Arial", FontSize->12], StyleBox["Progress in Optics", FontFamily->"Arial", FontSize->12, FontSlant->"Italic"], StyleBox[", edited by E. Wolf (Elsevier, Amsterdam, 2002), Vol. 43.", FontFamily->"Arial", FontSize->12] }], "Reference", CellTags->"Boy2002"], Cell[TextData[{ StyleBox["R.Y. Chiao, in ", FontFamily->"Arial", FontSize->12], StyleBox["Amazing light: a volume dedicated to Charles Hard Townes on his \ 80th birthday", FontFamily->"Arial", FontSize->12, FontSlant->"Italic"], StyleBox[", edited by R.Y. Chiao (Springer, New York, 1996).", FontFamily->"Arial", FontSize->12] }], "Reference", CellTags->"Chi96"], Cell[TextData[{ StyleBox["S. Chu and S. Wong, Physical Review Letters ", FontFamily->"Arial", FontSize->12], StyleBox["48", FontFamily->"Arial", FontSize->12, FontWeight->"Bold"], StyleBox[", 738 (1982).", FontFamily->"Arial", FontSize->12], "" }], "Reference", CellTags->"Chu82"], Cell[TextData[{ StyleBox["A. Dogariu, A. Kuzmich, and L.J. Wang, Physical Review A (Atomic, \ Molecular, and Optical Physics) ", FontFamily->"Arial", FontSize->12], StyleBox["63", FontFamily->"Arial", FontSize->12, FontWeight->"Bold"], StyleBox[", 053806/1 (2001).", FontFamily->"Arial", FontSize->12] }], "Reference", CellTags->"Dog2001"], Cell[TextData[{ StyleBox["C.G.B. Garrett and D.E. McCumber, Physical Review A (General \ Physics) ", FontFamily->"Arial", FontSize->12], StyleBox["1", FontFamily->"Arial", FontSize->12, FontWeight->"Bold"], StyleBox[", 305 (1970).", FontFamily->"Arial", FontSize->12] }], "Reference", CellTags->"Gar70"], Cell[TextData[{ StyleBox["D. J. Griffiths, ", FontFamily->"Arial", FontSize->12], StyleBox["Introduction to electrodynamics", FontFamily->"Arial", FontSize->12, FontSlant->"Italic"], StyleBox[" (Prentice Hall, Upper Saddle River, N.J., 1999).", FontFamily->"Arial", FontSize->12] }], "Reference", CellTags->"Gri99"], Cell[TextData[{ StyleBox["A. Kuzmich, A. Dogariu, L.J. Wang, et al., Physical Review \ Letters ", FontFamily->"Arial", FontSize->12], StyleBox["86", FontFamily->"Arial", FontSize->12, FontWeight->"Bold"], StyleBox[", 3925 (2001).", FontFamily->"Arial", FontSize->12] }], "Reference", CellTags->"Kuz2001"], Cell[TextData[{ StyleBox["L.D. Landau, E.M. Lifshits, and L.P. Pitaevskii, ", FontFamily->"Arial", FontSize->12], StyleBox["Electrodynamics of continuous media", FontFamily->"Arial", FontSize->12, FontSlant->"Italic"], StyleBox[" (Butterworth-Heinemann, Oxford England, 1995).", FontFamily->"Arial", FontSize->12] }], "Reference", CellTags->"Lan95"], Cell[TextData[{ StyleBox["S.G. Lipson, H. Lipson, and D.S. Tannhauser, ", FontFamily->"Arial", FontSize->12], StyleBox["Optical physics", FontFamily->"Arial", FontSize->12, FontSlant->"Italic"], StyleBox[" (Cambridge University Press, Cambridge 1995).", FontFamily->"Arial", FontSize->12] }], "Reference", CellTags->"Lip95"], Cell[TextData[{ StyleBox["A. B. Matsko, O. A. Kocharovskaya, Y. Rostovtsev, et al., in ", FontFamily->"Arial", FontSize->12], StyleBox["Advances in Atomic, Molecular, and Optical Physics", FontFamily->"Arial", FontSize->12, FontSlant->"Italic"], StyleBox[", edited by B. Bederson (Academic Press, New York, 2001), Vol. \ 46, p. 191.", FontFamily->"Arial", FontSize->12] }], "Reference", CellTags->"Mat2001"], Cell[TextData[{ StyleBox["K.T. McDonald, American Journal of Physics ", FontFamily->"Arial", FontSize->12], StyleBox["69", FontFamily->"Arial", FontSize->12, FontWeight->"Bold"], StyleBox[", 607 (2001).", FontFamily->"Arial", FontSize->12] }], "Reference", CellTags->"McD2001"], Cell[TextData[{ StyleBox["P. Sprangle, J.R. Penano, and B. Hafizi, Physical Review E \ (Statistical, Nonlinear, and Soft Matter Physics) ", FontFamily->"Arial", FontSize->12], StyleBox["64", FontFamily->"Arial", FontSize->12, FontWeight->"Bold"], StyleBox[", 026504/1 (2001).", FontFamily->"Arial", FontSize->12] }], "Reference", CellTags->"Spr2001"] }, Open ]] }, Open ]] }, FrontEndVersion->"4.0 for Microsoft Windows", ScreenRectangle->{{0, 1024}, {0, 695}}, AutoGeneratedPackage->None, WindowSize->{1014, 668}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, StyleDefinitions -> Notebook[{ Cell[CellGroupData[{ Cell["Style Definitions", "Subtitle"], Cell["\<\ Modify the definitions below to change the default appearance of all cells in \ a given style. Make modifications to any definition using commands in the Format menu.\ \>", "Text"], Cell[CellGroupData[{ Cell["Style Environment Names", "Section"], Cell[StyleData[All, "Working"], PageWidth->WindowWidth, CellBracketOptions->{"Color"->RGBColor[0.4, 0, 0.239994]}, ScriptMinSize->9, Background->RGBColor[1, 0.960006, 0.900008]], Cell[StyleData[All, "Presentation"], PageWidth->WindowWidth, CellBracketOptions->{"Color"->RGBColor[0.744411, 0.47155, 0.371542]}, ScriptMinSize->9], Cell[StyleData[All, "Printout"], PageWidth->PaperWidth, ShowCellLabel->False, ImageSize->{200, 200}, PrivateFontOptions->{"FontType"->"Outline"}] }, Closed]], Cell[CellGroupData[{ Cell["Notebook Options", "Section"], Cell["\<\ The options defined for the style below will be used at the Notebook level.\ \>", "Text"], Cell[StyleData["Notebook"], PageHeaders->{{Cell[ TextData[ { CounterBox[ "Page"]}], "PageNumber"], None, Cell[ TextData[ { ValueBox[ "FileName"]}], "Header"]}, {Cell[ TextData[ { ValueBox[ "FileName"]}], "Header"], None, Cell[ TextData[ { CounterBox[ "Page"]}], "PageNumber"]}}, CellFrameLabelMargins->6, StyleMenuListing->None] }, Closed]], Cell[CellGroupData[{ Cell["Styles for Headings", "Section"], Cell[CellGroupData[{ Cell[StyleData["Title"], CellFrame->{{0, 0}, {0, 0.25}}, CellMargins->{{18, 10}, {0, 20}}, CellGroupingRules->{"TitleGrouping", 0}, PageBreakBelow->False, CellFrameMargins->6, CellFrameColor->GrayLevel[0], InputAutoReplacements->{"TeX"->StyleBox[ RowBox[ {"T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "LaTeX"->StyleBox[ RowBox[ {"L", StyleBox[ AdjustmentBox[ "A", BoxMargins -> {{-0.36, -0.1}, {0, -0}}, BoxBaselineShift -> -0.2], FontSize -> Smaller], "T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "mma"->"Mathematica", "Mma"->"Mathematica", "MMA"->"Mathematica"}, LineSpacing->{0.95, 0}, CounterAssignments->{{"Section", 0}, {"Equation", 0}, {"Figure", 0}}, FontFamily->"Helvetica", FontSize->36, FontWeight->"Bold", FontColor->RGBColor[0.4, 0, 0.239994]], Cell[StyleData["Title", "Presentation"], CellMargins->{{20, 10}, {2, 20}}, CellFrameMargins->5, FontColor->RGBColor[0.744411, 0.47155, 0.371542]], Cell[StyleData["Title", "Printout"], CellMargins->{{18, 30}, {0, 0}}, CellFrameMargins->4, FontSize->24] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Subtitle"], CellMargins->{{20, 30}, {2, 10}}, CellGroupingRules->{"TitleGrouping", 10}, PageBreakBelow->False, CellFrameMargins->{{0, 4}, {8, 4}}, InputAutoReplacements->{"TeX"->StyleBox[ RowBox[ {"T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "LaTeX"->StyleBox[ RowBox[ {"L", StyleBox[ AdjustmentBox[ "A", BoxMargins -> {{-0.36, -0.1}, {0, -0}}, BoxBaselineShift -> -0.2], FontSize -> Smaller], "T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "mma"->"Mathematica", "Mma"->"Mathematica", "MMA"->"Mathematica"}, LineSpacing->{1, 0}, CounterAssignments->{{"Section", 0}, {"Equation", 0}, {"Figure", 0}}, FontFamily->"Helvetica", FontSize->24], Cell[StyleData["Subtitle", "Presentation"], CellMargins->{{20, 10}, {2, 10}}], Cell[StyleData["Subtitle", "Printout"], CellMargins->{{18, 30}, {0, 10}}, FontSize->18] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Section 1"], CellFrame->{{0, 0}, {0, 3}}, CellMargins->{{18, 10}, {10, 30}}, CellGroupingRules->{"SectionGrouping", 40}, PageBreakBelow->False, CellFrameMargins->{{0, 4}, {8, 4}}, CellFrameColor->RGBColor[0.4, 0, 0.239994], InputAutoReplacements->{"TeX"->StyleBox[ RowBox[ {"T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "LaTeX"->StyleBox[ RowBox[ {"L", StyleBox[ AdjustmentBox[ "A", BoxMargins -> {{-0.36, -0.1}, {0, -0}}, BoxBaselineShift -> -0.2], FontSize -> Smaller], "T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "mma"->"Mathematica", "Mma"->"Mathematica", "MMA"->"Mathematica"}, CounterIncrements->"Section", CounterAssignments->{{"Subsection", 0}, {"Subsubsection", 0}}, FontFamily->"Helvetica", FontSize->18, FontWeight->"Bold"], Cell[StyleData["Section 1", "Presentation"], CellMargins->{{20, 10}, {6, 30}}, CellFrameMargins->5], Cell[StyleData["Section 1", "Printout"], CellMargins->{{18, 30}, {4, 30}}, CellFrameMargins->4, FontSize->14] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Section"], CellFrame->{{0, 0}, {0, 0.25}}, CellMargins->{{18, 10}, {10, 30}}, CellGroupingRules->{"SectionGrouping", 40}, PageBreakBelow->False, CellFrameMargins->4, CellFrameColor->RGBColor[0.4, 0, 0.239994], InputAutoReplacements->{"TeX"->StyleBox[ RowBox[ {"T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "LaTeX"->StyleBox[ RowBox[ {"L", StyleBox[ AdjustmentBox[ "A", BoxMargins -> {{-0.36, -0.1}, {0, -0}}, BoxBaselineShift -> -0.2], FontSize -> Smaller], "T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "mma"->"Mathematica", "Mma"->"Mathematica", "MMA"->"Mathematica"}, CounterIncrements->"Section", CounterAssignments->{{"Subsection", 0}, {"Subsubsection", 0}}, FontFamily->"Helvetica", FontSize->18, FontWeight->"Bold"], Cell[StyleData["Section", "Presentation"], CellMargins->{{20, 10}, {0, 30}}], Cell[StyleData["Section", "Printout"], CellMargins->{{18, 30}, {0, 30}}, CellFrameMargins->5, FontSize->14] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Subsection"], CellDingbat->"\[FilledSquare]", CellMargins->{{38, 30}, {2, 20}}, CellGroupingRules->{"SectionGrouping", 50}, PageBreakBelow->False, InputAutoReplacements->{"TeX"->StyleBox[ RowBox[ {"T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "LaTeX"->StyleBox[ RowBox[ {"L", StyleBox[ AdjustmentBox[ "A", BoxMargins -> {{-0.36, -0.1}, {0, -0}}, BoxBaselineShift -> -0.2], FontSize -> Smaller], "T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "mma"->"Mathematica", "Mma"->"Mathematica", "MMA"->"Mathematica"}, CounterIncrements->"Subsection", CounterAssignments->{{"Subsubsection", 0}}, FontFamily->"Times", FontSize->14, FontWeight->"Bold"], Cell[StyleData["Subsection", "Presentation"], CellMargins->{{35, 30}, {0, 20}}], Cell[StyleData["Subsection", "Printout"], CellMargins->{{18, 30}, {0, 10}}, FontSize->12] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Subsubsection"], CellDingbat->"\[FilledSmallSquare]", CellMargins->{{55, 30}, {4, 10}}, CellGroupingRules->{"SectionGrouping", 60}, PageBreakBelow->False, InputAutoReplacements->{"TeX"->StyleBox[ RowBox[ {"T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "LaTeX"->StyleBox[ RowBox[ {"L", StyleBox[ AdjustmentBox[ "A", BoxMargins -> {{-0.36, -0.1}, {0, -0}}, BoxBaselineShift -> -0.2], FontSize -> Smaller], "T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "mma"->"Mathematica", "Mma"->"Mathematica", "MMA"->"Mathematica"}, CounterIncrements->"Subsubsection", FontFamily->"Times", FontSize->12, FontWeight->"Bold"], Cell[StyleData["Subsubsection", "Presentation"], CellMargins->{{31, 30}, {0, 12}}], Cell[StyleData["Subsubsection", "Printout"], CellMargins->{{18, 30}, {0, 12}}, FontSize->10] }, Closed]] }, Open ]], Cell[CellGroupData[{ Cell["Styles for Body Text", "Section"], Cell[CellGroupData[{ Cell[StyleData["Text"], PageWidth->550, CellMargins->{{55, 0}, {6, 6}}, CellSize->{7, Inherited}, InputAutoReplacements->{"TeX"->StyleBox[ RowBox[ {"T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "LaTeX"->StyleBox[ RowBox[ {"L", StyleBox[ AdjustmentBox[ "A", BoxMargins -> {{-0.36, -0.1}, {0, -0}}, BoxBaselineShift -> -0.2], FontSize -> Smaller], "T", AdjustmentBox[ "E", BoxMargins -> {{-0.075, -0.085}, {0, 0}}, BoxBaselineShift -> 0.5], "X"}]], "mma"->"Mathematica", "Mma"->"Mathematica", "MMA"->"Mathematica"}, TextJustification->1, Hyphenation->True, LineSpacing->{1, 2}, FontFamily->"Times", FontSize->14], Cell[StyleData["Text", "Presentation"], CellMargins->{{20, 10}, {6, 6}}, TextAlignment->Left, TextJustification->0, LineSpacing->{1.3, 0}], Cell[StyleData["Text", "Printout"], CellMargins->{{18, 4}, {4, 4}}, LineSpacing->{1, 3}] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Commentary"], PageWidth->550, CellMargins->{{120, 0}, {2, 6}}, TextAlignment->Left, TextJustification->1, Hyphenation->True, LineSpacing->{1, 2}, FontFamily->"Helvetica", FontSize->10, FontColor->RGBColor[0, 0, 0.4], Background->RGBColor[1, 0.749996, 0.749996]], Cell[StyleData["Commentary", "Presentation"], CellMargins->{{60, 30}, {2, 6}}, TextAlignment->Left, TextJustification->1, LineSpacing->{1.3, 0}, FontSize->12], Cell[StyleData["Commentary", "Printout"], CellMargins->{{18, 30}, {3, 0}}, TextAlignment->Left, LineSpacing->{1, 3}, FontFamily->"Times", FontSize->10] }, Closed]] }, Open ]], Cell[CellGroupData[{ Cell["Styles for Input/Output", "Section"], Cell["\<\ The cells in this section define styles used for input and output to the \ kernel. Be careful when modifying, renaming, or removing these styles, \ because the front end associates special meanings with these style names.\ \>", "Text"], Cell[CellGroupData[{ Cell[StyleData["Input"], CellFrame->{{1, 1}, {0, 1}}, CellMargins->{{55, 10}, {0, 0}}, Evaluatable->True, CellGroupingRules->"InputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, CellLabelPositioning->Automatic, CellLabelMargins->{{23, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultInputFormatType, HyphenationOptions->{"HyphenationCharacter"->"\[Continuation]"}, LanguageCategory->"Formula", FormatType->InputForm, ShowStringCharacters->True, NumberMarks->True, LinebreakAdjustments->{0.85, 2, 10, 0, 1}, FontSize->14, FontWeight->"Bold", Background->GrayLevel[0.966674]], Cell[StyleData["Input", "Presentation"], CellMargins->{{60, 10}, {0, 10}}, Background->GrayLevel[0.850004]], Cell[StyleData["Input", "Printout"], CellMargins->{{55, 10}, {0, 10}}, LinebreakAdjustments->{0.85, 2, 10, 1, 1}, FontSize->10, Background->GrayLevel[0.850004]] }, Closed]], Cell[StyleData["InlineInput"], Evaluatable->True, CellGroupingRules->"InputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, DefaultFormatType->DefaultInputFormatType, HyphenationOptions->{"HyphenationCharacter"->"\[Continuation]"}, AutoItalicWords->{}, FormatType->InputForm, ShowStringCharacters->True, NumberMarks->True, CounterIncrements->"Input", FontWeight->"Bold"], Cell[CellGroupData[{ Cell[StyleData["Output"], CellFrame->{{1, 1}, {1, 0}}, CellMargins->{{55, 10}, {15, 0}}, CellEditDuplicate->True, CellGroupingRules->"OutputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, CellLabelMargins->{{23, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultOutputFormatType, HyphenationOptions->{"HyphenationCharacter"->"\[Continuation]"}, LanguageCategory->"Formula", FormatType->InputForm, FontSize->12, Background->GrayLevel[0.850004]], Cell[StyleData["Output", "Presentation"], CellMargins->{{60, Inherited}, {10, 0}}], Cell[StyleData["Output", "Printout"], CellMargins->{{55, Inherited}, {10, 0}}, FontSize->10] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["InputOnly"], CellFrame->1, CellMargins->{{55, 10}, {15, 0}}, Evaluatable->True, CellGroupingRules->"InputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, CellLabelPositioning->Automatic, CellLabelMargins->{{23, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultInputFormatType, HyphenationOptions->{"HyphenationCharacter"->"\[Continuation]"}, LanguageCategory->"Formula", FormatType->InputForm, ShowStringCharacters->True, NumberMarks->True, LinebreakAdjustments->{0.85, 2, 10, 0, 1}, FontSize->12, FontWeight->"Bold", Background->GrayLevel[0.966674]], Cell[StyleData["InputOnly", "Presentation"], CellMargins->{{60, Inherited}, {10, 10}}], Cell[StyleData["InputOnly", "Printout"], CellMargins->{{55, Inherited}, {10, 10}}, LinebreakAdjustments->{0.85, 2, 10, 1, 1}, FontSize->10] }, Closed]], Cell[StyleData["InputOutdent"], CellFrame->{{1, 1}, {0, 1}}, CellMargins->{{55, 3}, {0, 0}}, Evaluatable->True, CellGroupingRules->"InputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, CellLabelPositioning->Automatic, CellLabelMargins->{{23, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultInputFormatType, HyphenationOptions->{"HyphenationCharacter"->"\[Continuation]"}, LanguageCategory->"Formula", FormatType->InputForm, ShowStringCharacters->True, NumberMarks->True, LinebreakAdjustments->{0.85, 2, 10, 0, 1}, FontSize->12, FontWeight->"Bold", Background->GrayLevel[0.966674]], Cell[RawData["\<\ Cell[StyleData[\"InputOnly\", \"Presentation\"], CellMargins->{{60, Inherited}, {10, 10}}]\ \>"], CellMargins->{{60, Inherited}, {10, 10}}], Cell[RawData["\<\ Cell[StyleData[\"InputOnly\", \"Printout\"], CellMargins->{{55, Inherited}, {10, 10}}, LinebreakAdjustments->{0.85, 2, 10, 1, 1}, FontSize->10]\ \>"], CellMargins->{{55, Inherited}, {10, 10}}, LinebreakAdjustments->{0.85, 2, 10, 1, 1}, FontSize->10], Cell[CellGroupData[{ Cell[StyleData["Message"], CellFrame->{{1, 1}, {0, 0}}, CellDingbat->"\[LongDash]", CellMargins->{{55, 10}, {0, 0}}, CellGroupingRules->"OutputGrouping", PageBreakWithin->False, GroupPageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, ShowCellLabel->False, CellLabelMargins->{{23, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultOutputFormatType, HyphenationOptions->{"HyphenationCharacter"->"\[Continuation]"}, FormatType->InputForm, StyleMenuListing->None, FontFamily->"Helvetica", FontSize->10, FontSlant->"Oblique", FontColor->GrayLevel[0], Background->GrayLevel[1]], Cell[StyleData["Message", "Presentation"], CellMargins->{{60, Inherited}, {Inherited, Inherited}}, LineSpacing->{1, 0}], Cell[StyleData["Message", "Printout"], CellMargins->{{55, Inherited}, {0, 0}}, FontSize->8] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Print"], CellMargins->{{55, 26}, {1, 6}}, CellGroupingRules->"OutputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, ShowCellLabel->False, CellLabelMargins->{{23, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultOutputFormatType, TextAlignment->Left, HyphenationOptions->{"HyphenationCharacter"->"\[Continuation]"}, FormatType->InputForm, StyleMenuListing->None, Background->GrayLevel[0.700008]], Cell[StyleData["Print", "Presentation"], CellMargins->{{60, Inherited}, {10, 2}}], Cell[StyleData["Print", "Printout"], CellMargins->{{54, Inherited}, {2, 6}}, FontSize->10] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Graphics"], CellFrame->{{1, 1}, {0, 0}}, CellMargins->{{55, 3}, {0, 0}}, CellGroupingRules->"GraphicsGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, ShowCellLabel->False, DefaultFormatType->DefaultOutputFormatType, FormatType->InputForm, ImageMargins->{{35, Inherited}, {Inherited, 0}}, StyleMenuListing->None, Background->GrayLevel[0.850004]], Cell[StyleData["Graphics", "Presentation"], CellMargins->{{60, Inherited}, {0, 0}}, ImageMargins->{{10, 10}, {10, 10}}], Cell[StyleData["Graphics", "Printout"], CellMargins->{{55, Inherited}, {0, 0}}, ImageSize->{0.0625, 0.0625}] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["CellLabel"], StyleMenuListing->None, FontFamily->"Helvetica", FontSize->10, FontSlant->"Oblique", FontColor->RGBColor[0.6, 0, 0.6]], Cell[StyleData["CellLabel", "Presentation"], CellMargins->{{18, Inherited}, {Inherited, Inherited}}], Cell[StyleData["CellLabel", "Printout"], CellMargins->{{0, Inherited}, {Inherited, Inherited}}, FontSize->8] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Unique Styles", "Section"], Cell[CellGroupData[{ Cell[StyleData["Author"], CellMargins->{{20, 30}, {45, 5}}, CellGroupingRules->{"TitleGrouping", 20}, PageBreakBelow->False, CellFrameMargins->{{0, 4}, {8, 4}}, LineSpacing->{1, 0}, CounterAssignments->{{"Section", 0}, {"Equation", 0}, {"Figure", 0}}, FontFamily->"Helvetica", FontSize->14, FontSlant->"Italic"], Cell[StyleData["Author", "Presentation"], CellMargins->{{20, 30}, {45, 10}}], Cell[StyleData["Author", "Printout"], CellMargins->{{18, 30}, {45, 5}}] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Abstract"], CellFrame->False, CellMargins->{{18, 140}, {4, 30}}, Hyphenation->True, LineSpacing->{0.9, 0}, FontFamily->"Times", FontSize->12], Cell[StyleData["Abstract", "Presentation"], CellFrame->True, CellMargins->{{20, 10}, {Inherited, 30}}], Cell[StyleData["Abstract", "Printout"], LineSpacing->{1, 2}, FontSize->10] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Caption"], PageWidth->WindowWidth, CellFrame->{{1, 1}, {1, 0}}, CellMargins->{{55, 17}, {5, 0}}, PageBreakAbove->False, Hyphenation->True, FontFamily->"Helvetica", FontSize->9, FontColor->GrayLevel[0], Background->GrayLevel[0.850004]], Cell[StyleData["Caption", "Presentation"], PageWidth->PaperWidth, CellMargins->{{60, 65}, {6, 4}}, FontSize->10], Cell[StyleData["Caption", "Printout"], PageWidth->PaperWidth, CellMargins->{{55, 55}, {5, 4}}, LineSpacing->{1, 2}, FontSize->8] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["CaptionOutdent"], PageWidth->WindowWidth, CellFrame->{{1, 1}, {1, 0}}, CellMargins->{{55, 10}, {5, 0}}, PageBreakAbove->False, Hyphenation->True, FontFamily->"Helvetica", FontSize->9, FontColor->GrayLevel[0], Background->GrayLevel[0.850004]], Cell[StyleData["Caption", "Presentation"], PageWidth->WindowWidth, CellMargins->{{60, 65}, {6, 4}}, FontSize->10], Cell[StyleData["Caption", "Printout"], PageWidth->WindowWidth, CellMargins->{{55, 55}, {5, 4}}, LineSpacing->{1, 2}, FontSize->8] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Reference"], CellMargins->{{24, 40}, {6, 6}}, TextJustification->1, Hyphenation->True, LineSpacing->{1, 0}, FontFamily->"Times"], Cell[StyleData["Reference", "Presentation"], CellMargins->{{20, 40}, {Inherited, 6}}, TextJustification->0, LineSpacing->{1, 4}, FontSize->12], Cell[StyleData["Reference", "Printout"], CellMargins->{{18, 4}, {4, 4}}, FontSize->9] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["PictureGroup"], CellFrame->{{1, 1}, {0, 0}}, CellMargins->{{55, 20}, {0, 0}}, CellGroupingRules->"GraphicsGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, ShowCellLabel->False, ImageMargins->{{35, Inherited}, {Inherited, 0}}, StyleMenuListing->None, Background->GrayLevel[0.850004]], Cell[StyleData["PictureGroup", "Presentation"], CellMargins->{{60, Inherited}, {0, 0}}, ImageMargins->{{10, 10}, {10, 10}}], Cell[StyleData["PictureGroup", "Printout"], CellMargins->{{55, Inherited}, {0, 0}}, ImageSize->{0.0625, 0.0625}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Hyperlink Styles", "Section"], Cell["\<\ The cells below define styles useful for making hypertext ButtonBoxes. The \ \"Hyperlink\" style is for links within the same Notebook, or between \ Notebooks.\ \>", "Text"], Cell[CellGroupData[{ Cell[StyleData["Hyperlink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`NotebookLocate[ #2]}]&), Active->True, ButtonNote->ButtonData}], Cell[StyleData["Hyperlink", "Presentation"], FontSize->16], Cell[StyleData["Hyperlink", "Printout"], FontSize->10, FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Closed]], Cell["\<\ The following styles are for linking automatically to the on-line help \ system.\ \>", "Text"], Cell[CellGroupData[{ Cell[StyleData["MainBookLink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`HelpBrowserLookup[ "MainBook", #]}]&), Active->True, ButtonFrame->"None"}], Cell[StyleData["MainBookLink", "Presentation"], FontSize->16], Cell[StyleData["MainBookLink", "Printout"], FontSize->10, FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["AddOnsLink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontFamily->"Courier", FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`HelpBrowserLookup[ "AddOns", #]}]&), Active->True, ButtonFrame->"None"}], Cell[StyleData["AddOnsLink", "Presentation"], FontSize->16], Cell[StyleData["AddOnsLink", "Printout"], FontSize->10, FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["RefGuideLink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontFamily->"Courier", FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`HelpBrowserLookup[ "RefGuide", #]}]&), Active->True, ButtonFrame->"None"}], Cell[StyleData["RefGuideLink", "Presentation"], FontSize->16], Cell[StyleData["RefGuideLink", "Printout"], FontSize->10, FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["GettingStartedLink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`HelpBrowserLookup[ "GettingStarted", #]}]&), Active->True, ButtonFrame->"None"}], Cell[StyleData["GettingStartedLink", "Presentation"], FontSize->16], Cell[StyleData["GettingStartedLink", "Printout"], FontSize->10, FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["OtherInformationLink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`HelpBrowserLookup[ "OtherInformation", #]}]&), Active->True, ButtonFrame->"None"}], Cell[StyleData["OtherInformationLink", "Presentation"], FontSize->16], Cell[StyleData["OtherInformationLink", "Printout"], FontSize->10, FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Palette Styles", "Section"], Cell["\<\ The cells below define styles that define standard ButtonFunctions, for use \ in palette buttons.\ \>", "Text"], Cell[StyleData["Paste"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`NotebookApply[ FrontEnd`InputNotebook[ ], #, After]}]&)}], Cell[StyleData["Evaluate"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`NotebookApply[ FrontEnd`InputNotebook[ ], #, All], SelectionEvaluate[ FrontEnd`InputNotebook[ ], All]}]&)}], Cell[StyleData["EvaluateCell"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`NotebookApply[ FrontEnd`InputNotebook[ ], #, All], FrontEnd`SelectionMove[ FrontEnd`InputNotebook[ ], All, Cell, 1], FrontEnd`SelectionEvaluateCreateCell[ FrontEnd`InputNotebook[ ], All]}]&)}], Cell[StyleData["CopyEvaluate"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`SelectionCreateCell[ FrontEnd`InputNotebook[ ], All], FrontEnd`NotebookApply[ FrontEnd`InputNotebook[ ], #, All], FrontEnd`SelectionEvaluate[ FrontEnd`InputNotebook[ ], All]}]&)}], Cell[StyleData["CopyEvaluateCell"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`SelectionCreateCell[ FrontEnd`InputNotebook[ ], All], FrontEnd`NotebookApply[ FrontEnd`InputNotebook[ ], #, All], FrontEnd`SelectionEvaluateCreateCell[ FrontEnd`InputNotebook[ ], All]}]&)}] }, Closed]], Cell[CellGroupData[{ Cell["Formulas and Programming", "Section"], Cell[CellGroupData[{ Cell[StyleData["DisplayFormula"], CellMargins->{{55, Inherited}, {Inherited, Inherited}}, CellHorizontalScrolling->True, DefaultFormatType->DefaultInputFormatType, HyphenationOptions->{"HyphenationCharacter"->"\[Continuation]"}, LanguageCategory->"Formula", ScriptLevel->0, SingleLetterItalics->True, StyleMenuListing->None, UnderoverscriptBoxOptions->{LimitsPositioning->True}], Cell[StyleData["DisplayFormula", "Presentation"], CellMargins->{{60, 10}, {Inherited, 6}}, LineSpacing->{1, 5}], Cell[StyleData["DisplayFormula", "Printout"], CellMargins->{{18, 4}, {4, 4}}] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["ChemicalFormula"], CellMargins->{{55, 10}, {Inherited, 0}}, CellHorizontalScrolling->True, DefaultFormatType->DefaultInputFormatType, TextJustification->1, HyphenationOptions->{"HyphenationCharacter"->"\[Continuation]"}, LanguageCategory->"Formula", AutoSpacing->False, ScriptLevel->1, ScriptBaselineShifts->{0.6, Automatic}, SingleLetterItalics->False, ZeroWidthTimes->True], Cell[StyleData["ChemicalFormula", "Presentation"], CellMargins->{{60, 10}, {Inherited, 6}}], Cell[StyleData["ChemicalFormula", "Printout"], CellMargins->{{18, 4}, {4, 4}}, LineSpacing->{1, 3}, FontSize->10] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Program"], CellMargins->{{55, 10}, {Inherited, 0}}, CellHorizontalScrolling->True, LanguageCategory->"Formula", FontFamily->"Courier"], Cell[StyleData["Program", "Presentation"], CellMargins->{{60, 10}, {Inherited, 6}}], Cell[StyleData["Program", "Printout"], CellMargins->{{18, 4}, {4, 4}}, LineSpacing->{1, 3}, FontSize->9.5] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Styles for Automatic Numbering", "Section"], Cell["\<\ The following styles are useful for numbered equations, figures, etc. They \ automatically give the cell a FrameLabel containing a reference to a \ particular counter, and also increment that counter.\ \>", "Text"], Cell[CellGroupData[{ Cell[StyleData["NumberedEquation"], CellMargins->{{55, 85}, {Inherited, Inherited}}, CellFrameLabels->{{None, Cell[ TextData[ {"(", CounterBox[ "NumberedEquation"], ")"}]]}, {None, None}}, DefaultFormatType->DefaultInputFormatType, HyphenationOptions->{"HyphenationCharacter"->"\[Continuation]"}, CounterIncrements->"NumberedEquation", FormatTypeAutoConvert->False, FontFamily->"Times"], Cell[StyleData["NumberedEquation", "Presentation"], CellMargins->{{60, 10}, {Inherited, 6}}, LineSpacing->{1, 0}], Cell[StyleData["NumberedEquation", "Printout"], CellMargins->{{18, 4}, {4, 4}}, FontSize->8] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["NumberedFigure"], CellMargins->{{55, 95}, {Inherited, Inherited}}, CellFrameLabels->{{None, None}, {Cell[ TextData[ {"Figure ", CounterBox[ "NumberedFigure"]}], FontWeight -> "Bold"], None}}, CounterIncrements->"NumberedFigure", FormatTypeAutoConvert->False, FontFamily->"Times"], Cell[StyleData["NumberedFigure", "Presentation"], CellMargins->{{60, 80}, {Inherited, 6}}, LineSpacing->{1, 0}], Cell[StyleData["NumberedFigure", "Printout"], CellMargins->{{18, 80}, {4, 4}}, FontSize->8] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["NumberedTable"], CellMargins->{{55, 95}, {Inherited, Inherited}}, CellFrameLabels->{{None, None}, {Cell[ TextData[ {"Table ", CounterBox[ "NumberedTable"]}], FontWeight -> "Bold"], None}}, TextAlignment->Center, CounterIncrements->"NumberedTable", FormatTypeAutoConvert->False, FontFamily->"Times"], Cell[StyleData["NumberedTable", "Presentation"], CellMargins->{{60, 80}, {Inherited, 6}}, LineSpacing->{1, 0}], Cell[StyleData["NumberedTable", "Printout"], CellMargins->{{18, 80}, {4, 4}}, FontSize->8] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Styles for Headers and Footers", "Section"], Cell[StyleData["Header"], StyleMenuListing->None, FontFamily->"Helvetica", FontSize->9], Cell[StyleData["PageNumber"], StyleMenuListing->None, FontFamily->"Helvetica", FontSize->9, FontWeight->"Bold"], Cell[StyleData["Footer"], TextAlignment->Center, StyleMenuListing->None, FontFamily->"Helvetica", FontSize->9] }, Closed]] }, Open ]] }] ] (*********************************************************************** 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->{ "conclusions"->{ Cell[71278, 2090, 57, 1, 39, "Subtitle", CellTags->"conclusions"]}, "Boy2002"->{ Cell[74070, 2174, 363, 12, 27, "Reference", CellTags->"Boy2002"]}, "Chi96"->{ Cell[74436, 2188, 398, 13, 27, "Reference", CellTags->"Chi96"]}, "Chu82"->{ Cell[74837, 2203, 318, 13, 27, "Reference", CellTags->"Chu82"]}, "Dog2001"->{ Cell[75158, 2218, 375, 13, 27, "Reference", CellTags->"Dog2001"]}, "Gar70"->{ Cell[75536, 2233, 339, 13, 27, "Reference", CellTags->"Gar70"]}, "Gri99"->{ Cell[75878, 2248, 351, 12, 27, "Reference", CellTags->"Gri99"]}, "Kuz2001"->{ Cell[76232, 2262, 340, 13, 27, "Reference", CellTags->"Kuz2001"]}, "Lan95"->{ Cell[76575, 2277, 385, 12, 27, "Reference", CellTags->"Lan95"]}, "Lip95"->{ Cell[76963, 2291, 360, 12, 27, "Reference", CellTags->"Lip95"]}, "Mat2001"->{ Cell[77326, 2305, 443, 13, 42, "Reference", CellTags->"Mat2001"]}, "McD2001"->{ Cell[77772, 2320, 313, 12, 27, "Reference", CellTags->"McD2001"]}, "Spr2001"->{ Cell[78088, 2334, 386, 13, 27, "Reference", CellTags->"Spr2001"]} } *) (*CellTagsIndex CellTagsIndex->{ {"conclusions", 112847, 3482}, {"Boy2002", 112940, 3485}, {"Chi96", 113030, 3488}, {"Chu82", 113118, 3491}, {"Dog2001", 113208, 3494}, {"Gar70", 113298, 3497}, {"Gri99", 113386, 3500}, {"Kuz2001", 113476, 3503}, {"Lan95", 113566, 3506}, {"Lip95", 113654, 3509}, {"Mat2001", 113744, 3512}, {"McD2001", 113836, 3515}, {"Spr2001", 113928, 3518} } *) (*NotebookFileOutline Notebook[{ Cell[1717, 49, 139, 3, 53, "Text"], Cell[1859, 54, 95, 4, 98, "Author"], Cell[CellGroupData[{ Cell[1979, 62, 32, 0, 39, "Subtitle"], Cell[CellGroupData[{ Cell[2036, 66, 31, 0, 70, "Section 1"], Cell[2070, 68, 2002, 57, 192, "Text"], Cell[4075, 127, 973, 26, 138, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[5085, 158, 29, 0, 70, "Section 1"], Cell[5117, 160, 1196, 19, 246, "Text"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[6362, 185, 65, 1, 39, "Subtitle", InitializationCell->True], Cell[6430, 188, 70, 1, 30, "Text", InitializationCell->True], Cell[6503, 191, 208, 4, 91, "InputOnly", InitializationCell->True], Cell[6714, 197, 67, 1, 30, "Text", InitializationCell->True], Cell[6784, 200, 177, 3, 71, "InputOnly", InitializationCell->True], Cell[6964, 205, 68, 1, 30, "Text", InitializationCell->True], Cell[7035, 208, 705, 16, 71, "InputOnly", InitializationCell->True], Cell[7743, 226, 298, 5, 71, "InputOnly", InitializationCell->True], Cell[8044, 233, 156, 3, 51, "InputOnly", InitializationCell->True], Cell[8203, 238, 87, 1, 30, "Text", InitializationCell->True], Cell[8293, 241, 138, 3, 51, "InputOnly", InitializationCell->True] }, Closed]], Cell[CellGroupData[{ Cell[8468, 249, 67, 0, 33, "Subtitle"], Cell[8538, 251, 56, 0, 70, "Section 1"], Cell[CellGroupData[{ Cell[8619, 255, 48, 0, 67, "Section"], Cell[8670, 257, 293, 7, 66, "Text"], Cell[CellGroupData[{ Cell[8988, 268, 159, 4, 36, "Input", InitializationCell->True], Cell[9150, 274, 96, 2, 50, "Output"] }, Open ]], Cell[9261, 279, 21, 0, 30, "Text"], Cell[CellGroupData[{ Cell[9307, 283, 98, 2, 46, "Input", InitializationCell->True], Cell[9408, 287, 52, 1, 60, "Output"] }, Open ]], Cell[9475, 291, 372, 11, 48, "Text"], Cell[CellGroupData[{ Cell[9872, 306, 135, 3, 35, "Input", InitializationCell->True], Cell[10010, 311, 106, 2, 49, "Output"] }, Open ]], Cell[10131, 316, 186, 4, 35, "Input"], Cell[10320, 322, 223, 6, 50, "Caption"] }, Open ]], Cell[CellGroupData[{ Cell[10580, 333, 68, 2, 67, "Section"], Cell[10651, 337, 391, 12, 66, "Text"], Cell[CellGroupData[{ Cell[11067, 353, 174, 3, 35, "Input", InitializationCell->True], Cell[11244, 358, 88, 1, 71, "Output"] }, Open ]], Cell[11347, 362, 128, 4, 34, "Text"], Cell[CellGroupData[{ Cell[11500, 370, 104, 2, 35, "Input", InitializationCell->True], Cell[11607, 374, 49, 1, 62, "Output"] }, Open ]], Cell[11671, 378, 51, 0, 30, "Text"], Cell[CellGroupData[{ Cell[11747, 382, 123, 3, 35, "Input", InitializationCell->True], Cell[11873, 387, 42, 1, 62, "Output"] }, Open ]], Cell[11930, 391, 130, 5, 30, "Text"], Cell[CellGroupData[{ Cell[12085, 400, 42, 1, 35, "Input"], Cell[12130, 403, 35, 1, 49, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[12214, 410, 63, 0, 67, "Section"], Cell[12280, 412, 494, 16, 66, "Text"], Cell[CellGroupData[{ Cell[12799, 432, 209, 4, 49, "Input", InitializationCell->True], Cell[13011, 438, 158, 2, 72, "Output"] }, Open ]], Cell[13184, 443, 229, 6, 48, "Text"], Cell[CellGroupData[{ Cell[13438, 453, 297, 6, 59, "Input", InitializationCell->True], Cell[13738, 461, 74, 1, 60, "Output"] }, Open ]], Cell[13827, 465, 132, 5, 30, "Text"], Cell[CellGroupData[{ Cell[13984, 474, 191, 3, 59, "Input", InitializationCell->True], Cell[14178, 479, 136, 2, 73, "Output"] }, Open ]], Cell[14329, 484, 459, 9, 84, "Commentary"], Cell[14791, 495, 125, 3, 48, "Text"], Cell[14919, 500, 1581, 29, 191, "InputOnly", InitializationCell->True], Cell[16503, 531, 157, 2, 35, "InputOutdent"], Cell[16663, 535, 170, 2, 36, "CaptionOutdent"], Cell[16836, 539, 60, 0, 30, "Text"], Cell[CellGroupData[{ Cell[16921, 543, 164, 4, 35, "Input", InitializationCell->True], Cell[17088, 549, 81, 1, 69, "Output"] }, Open ]], Cell[17184, 553, 577, 18, 84, "Text"], Cell[CellGroupData[{ Cell[17786, 575, 131, 3, 35, "Input", InitializationCell->True], Cell[17920, 580, 177, 3, 72, "Output"] }, Open ]], Cell[18112, 586, 50, 0, 30, "Text"], Cell[CellGroupData[{ Cell[18187, 590, 149, 3, 35, "Input", InitializationCell->True], Cell[18339, 595, 241, 3, 75, "Output"] }, Open ]], Cell[18595, 601, 264, 6, 66, "Text"], Cell[CellGroupData[{ Cell[18884, 611, 113, 2, 35, "Input"], Cell[19000, 615, 71, 1, 49, "Output"] }, Open ]], Cell[19086, 619, 116, 3, 48, "Text"], Cell[19205, 624, 2181, 41, 511, "InputOnly", InitializationCell->True], Cell[21389, 667, 871, 16, 156, "Text"], Cell[22263, 685, 235, 4, 35, "Input"], Cell[22501, 691, 259, 6, 50, "Caption"], Cell[22763, 699, 16, 0, 30, "Text"] }, Open ]], Cell[22794, 702, 49, 0, 70, "Section 1"], Cell[CellGroupData[{ Cell[22868, 706, 45, 0, 67, "Section"], Cell[CellGroupData[{ Cell[22938, 710, 34, 0, 41, "Subsection"], Cell[22975, 712, 186, 5, 48, "Text"], Cell[CellGroupData[{ Cell[23186, 721, 224, 4, 63, "Input", InitializationCell->True], Cell[23413, 727, 146, 2, 84, "Output"] }, Open ]], Cell[23574, 732, 339, 8, 66, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[23950, 745, 42, 0, 41, "Subsection"], Cell[CellGroupData[{ Cell[24017, 749, 40, 0, 31, "Subsubsection"], Cell[24060, 751, 205, 5, 48, "Text"], Cell[CellGroupData[{ Cell[24290, 760, 296, 6, 35, "Input"], Cell[24589, 768, 236, 3, 92, "Output"] }, Open ]], Cell[24840, 774, 72, 0, 30, "Text"], Cell[CellGroupData[{ Cell[24937, 778, 124, 2, 35, "Input"], Cell[25064, 782, 190, 3, 92, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[25291, 790, 76, 1, 35, "Input"], Cell[25370, 793, 209, 3, 92, "Output"] }, Open ]], Cell[25594, 799, 165, 3, 48, "Text"], Cell[25762, 804, 406, 8, 55, "Input"], Cell[26171, 814, 198, 4, 50, "Caption"] }, Open ]], Cell[CellGroupData[{ Cell[26406, 823, 42, 0, 31, "Subsubsection"], Cell[26451, 825, 62, 0, 30, "Text"], Cell[CellGroupData[{ Cell[26538, 829, 286, 6, 35, "Input"], Cell[26827, 837, 287, 4, 100, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[27151, 846, 106, 2, 35, "Input"], Cell[27260, 850, 239, 3, 100, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[27536, 858, 109, 2, 35, "Input"], Cell[27648, 862, 277, 4, 100, "Output"] }, Open ]], Cell[27940, 869, 410, 8, 55, "Input"], Cell[28353, 879, 200, 4, 50, "Caption"] }, Open ]], Cell[CellGroupData[{ Cell[28590, 888, 41, 0, 31, "Subsubsection"], Cell[28634, 890, 38, 0, 30, "Text"], Cell[CellGroupData[{ Cell[28697, 894, 232, 5, 35, "Input"], Cell[28932, 901, 941, 13, 94, "Output"] }, Open ]], Cell[29888, 917, 411, 8, 55, "Input"], Cell[30302, 927, 199, 4, 50, "Caption"], Cell[30504, 933, 399, 6, 102, "Text"], Cell[30906, 941, 498, 8, 99, "Commentary"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[31453, 955, 59, 0, 41, "Subsection"], Cell[31515, 957, 137, 3, 48, "Text"], Cell[CellGroupData[{ Cell[31677, 964, 388, 9, 55, "Input", InitializationCell->True], Cell[32068, 975, 171, 3, 63, "Output"] }, Open ]], Cell[32254, 981, 62, 0, 30, "Text"], Cell[CellGroupData[{ Cell[32341, 985, 116, 2, 35, "Input"], Cell[32460, 989, 122, 2, 62, "Output"] }, Open ]], Cell[32597, 994, 475, 10, 55, "Input"], Cell[33075, 1006, 173, 4, 50, "Caption"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[33297, 1016, 80, 2, 67, "Section"], Cell[33380, 1020, 71, 0, 30, "Text"], Cell[CellGroupData[{ Cell[33476, 1024, 159, 3, 35, "Input", InitializationCell->True], Cell[33638, 1029, 472, 15, 79, "Output"] }, Open ]], Cell[34125, 1047, 487, 17, 48, "Text"], Cell[CellGroupData[{ Cell[34637, 1068, 134, 3, 35, "Input", InitializationCell->True], Cell[34774, 1073, 517, 12, 49, "Output"] }, Open ]], Cell[35306, 1088, 412, 11, 52, "Text"], Cell[CellGroupData[{ Cell[35743, 1103, 136, 3, 35, "Input", InitializationCell->True], Cell[35882, 1108, 1025, 29, 79, "Output"] }, Open ]], Cell[36922, 1140, 1143, 25, 75, "Input"], Cell[38068, 1167, 223, 4, 50, "Caption"], Cell[38294, 1173, 217, 4, 66, "Text"], Cell[CellGroupData[{ Cell[38536, 1181, 114, 2, 35, "Input", InitializationCell->True], Cell[38653, 1185, 1176, 25, 55, "Output"] }, Open ]], Cell[39844, 1213, 28, 0, 30, "Text"], Cell[CellGroupData[{ Cell[39897, 1217, 89, 2, 35, "Input", InitializationCell->True], Cell[39989, 1221, 502, 12, 62, "Output"] }, Open ]], Cell[40506, 1236, 42, 0, 30, "Text"], Cell[CellGroupData[{ Cell[40573, 1240, 81, 1, 35, "Input"], Cell[40657, 1243, 438, 10, 62, "Output"] }, Open ]], Cell[41110, 1256, 44, 0, 30, "Text"], Cell[CellGroupData[{ Cell[41179, 1260, 133, 3, 35, "Input", InitializationCell->True], Cell[41315, 1265, 392, 10, 62, "Output"] }, Open ]], Cell[41722, 1278, 152, 5, 30, "Text"], Cell[41877, 1285, 1145, 25, 75, "Input"], Cell[43025, 1312, 392, 7, 64, "Caption"], Cell[43420, 1321, 313, 9, 48, "Text"], Cell[CellGroupData[{ Cell[43758, 1334, 168, 3, 35, "Input", InitializationCell->True], Cell[43929, 1339, 282, 7, 64, "Output"] }, Open ]], Cell[44226, 1349, 147, 5, 30, "Text"], Cell[CellGroupData[{ Cell[44398, 1358, 93, 2, 35, "Input", InitializationCell->True], Cell[44494, 1362, 163, 4, 64, "Output"] }, Open ]], Cell[44672, 1369, 179, 8, 30, "Text"], Cell[CellGroupData[{ Cell[44876, 1381, 110, 2, 35, "Input", InitializationCell->True], Cell[44989, 1385, 62, 1, 62, "Output"] }, Open ]], Cell[45066, 1389, 93, 2, 30, "Text"], Cell[CellGroupData[{ Cell[45184, 1395, 100, 2, 35, "Input", InitializationCell->True], Cell[45287, 1399, 265, 7, 62, "Output"] }, Open ]], Cell[45567, 1409, 346, 6, 84, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[45950, 1420, 67, 2, 67, "Section"], Cell[CellGroupData[{ Cell[46042, 1426, 42, 0, 41, "Subsection"], Cell[46087, 1428, 808, 12, 192, "Text"], Cell[46898, 1442, 135, 5, 30, "Text"], Cell[CellGroupData[{ Cell[47058, 1451, 140, 2, 35, "Input"], Cell[47201, 1455, 174, 3, 67, "Output"] }, Open ]], Cell[47390, 1461, 112, 3, 30, "Text"], Cell[CellGroupData[{ Cell[47527, 1468, 83, 1, 35, "Input"], Cell[47613, 1471, 136, 2, 81, "Output"] }, Open ]], Cell[47764, 1476, 72, 2, 30, "Text"], Cell[CellGroupData[{ Cell[47861, 1482, 245, 3, 35, "Input"], Cell[48109, 1487, 228, 3, 81, "Output"] }, Open ]], Cell[48352, 1493, 138, 5, 30, "Text"], Cell[48493, 1500, 347, 7, 66, "Text"], Cell[CellGroupData[{ Cell[48865, 1511, 145, 3, 35, "Input"], Cell[49013, 1516, 590, 9, 70, "Output"] }, Open ]], Cell[49618, 1528, 302, 7, 66, "Text"], Cell[CellGroupData[{ Cell[49945, 1539, 140, 3, 35, "Input"], Cell[50088, 1544, 106, 2, 67, "Output"] }, Open ]], Cell[50209, 1549, 330, 8, 66, "Text"], Cell[50542, 1559, 223, 3, 35, "InputOutdent"], Cell[50768, 1564, 249, 3, 36, "CaptionOutdent"] }, Open ]], Cell[CellGroupData[{ Cell[51054, 1572, 36, 0, 41, "Subsection"], Cell[51093, 1574, 361, 6, 84, "Text"], Cell[CellGroupData[{ Cell[51479, 1584, 77, 2, 35, "Input", InitializationCell->True], Cell[51559, 1588, 118, 2, 69, "Output"] }, Open ]], Cell[51692, 1593, 167, 5, 30, "Text"], Cell[CellGroupData[{ Cell[51884, 1602, 144, 3, 35, "Input", InitializationCell->True], Cell[52031, 1607, 77, 1, 68, "Output"] }, Open ]], Cell[52123, 1611, 592, 14, 102, "Text"], Cell[52718, 1627, 537, 8, 114, "Commentary"] }, Open ]], Cell[CellGroupData[{ Cell[53292, 1640, 40, 0, 41, "Subsection"], Cell[53335, 1642, 306, 6, 66, "Text"], Cell[CellGroupData[{ Cell[53666, 1652, 228, 6, 36, "Input", InitializationCell->True], Cell[53897, 1660, 378, 5, 59, "Output"] }, Open ]], Cell[54290, 1668, 247, 8, 30, "Text"], Cell[CellGroupData[{ Cell[54562, 1680, 159, 4, 35, "Input", InitializationCell->True], Cell[54724, 1686, 245, 3, 59, "Output"] }, Open ]], Cell[54984, 1692, 63, 0, 30, "Text"], Cell[55050, 1694, 353, 7, 75, "Input", InitializationCell->True], Cell[55406, 1703, 563, 10, 95, "Input", InitializationCell->True], Cell[55972, 1715, 3725, 63, 578, "InputOnly", InitializationCell->True] }, Open ]], Cell[CellGroupData[{ Cell[59734, 1783, 48, 0, 41, "Subsection"], Cell[CellGroupData[{ Cell[59807, 1787, 43, 0, 31, "Subsubsection"], Cell[59853, 1789, 267, 8, 30, "Text"], Cell[CellGroupData[{ Cell[60145, 1801, 114, 2, 35, "Input"], Cell[60262, 1805, 38, 1, 49, "Output"] }, Open ]], Cell[60315, 1809, 783, 16, 156, "Text"], Cell[61101, 1827, 210, 4, 35, "Input"], Cell[61314, 1833, 255, 6, 50, "Caption"] }, Open ]], Cell[CellGroupData[{ Cell[61606, 1844, 48, 0, 31, "Subsubsection"], Cell[61657, 1846, 152, 5, 30, "Text"], Cell[CellGroupData[{ Cell[61834, 1855, 105, 2, 35, "Input"], Cell[61942, 1859, 47, 1, 60, "Output"] }, Ope