2018-04-14 22:36

8 files changed, 697 insertions, 443 deletions

diff --git a/MX2/chapter.tex b/MX2/chapter.tex
index 6e5bdb1..98c37cc 100644
--- a/MX2/chapter.tex
+++ b/MX2/chapter.tex
@@ -91,26 +91,11 @@ frequencies $\omega_1$ and $\omega_2$.  %
 The $\omega_2$ beam is split into two beams, denoted by $\omega_2$ and $\omega_2^\prime$.  %
 These three beams are focused onto the MoS\textsubscript{2} thin film at angles, creating an output
 beam in the phase-matched direction
-$\mathbf{k}_{\mathrm{out}}=\mathbf{k}_1-\mathbf{k}_2+\mathbf{k}_{2^\prime}$ where $\mathbf{k}$ is the wave
-vector for each beam and the subscripts label the excitation frequencies.  %
+$\mathbf{k}_{\mathrm{out}}=\mathbf{k}_1-\mathbf{k}_2+\mathbf{k}_{2^\prime}$ where $\mathbf{k}$ is
+the wave vector for each beam and the subscripts label the excitation frequencies.  %
 Multidimensional spectra result from measuring the output intensity dependence on frequency and
 delay times.  %
 
-\begin{figure}
-  \includegraphics[width=0.5\textwidth]{MX2/01}
-  \caption[CMDS tutorial]{
-    (a) Example delays of the $\omega_1$, $\omega_2$, and $\omega_{2^\prime}$ excitation pulses.
-    (b) Dependence of the output intensity on the $\tau_{22^\prime}$ and $\tau_{21}$ time delays
-    for $\omega_1=\omega_2$.
-    The solid lines define the regions for the six different time orderings of the $\omega_1$,
-    $\omega_2$, and $\omega_{2^\prime}$ excitation pulses.
-    We have developed a convention for numbering these time orderings, as shown.
-    (c) Diagram of the band structure of MoS\textsubscript{2} at the $K$ point.
-    The A and B exciton transitions are shown. (d) Two dimensional frequency-frequency plot
-    labeling two diagonal and cross-peak features for the A and B excitons.}
-  \label{fig:Czech01}
-\end{figure}
-
 \autoref{fig:Czech01} introduces our conventions for representing multidimensional spectra.  %
 \autoref{fig:Czech01}b,d are simulated data.  %
 \autoref{fig:Czech01}a shows one of the six time orderings of the three excitation pulses where
@@ -155,14 +140,23 @@ The intensity of the cross-peaks depends on the importance of state filling and
 relaxation of hot A excitons as well as the presence of interband population trnasfer of the A and
 B exciton states.  %
 
-\section{Methods}  % ==============================================================================
-
 \begin{figure}
-  \includegraphics[width=\textwidth]{MX2/S1}
-  \caption{Schemiatic of the synthetic setup used for Mo thin film sulfidation reactions.}
-  \label{fig:CzechS1}
+  \includegraphics[width=0.5\textwidth]{MX2/01}
+  \caption[Introduction to CMDS.]{
+    (a) Example delays of the $\omega_1$, $\omega_2$, and $\omega_{2^\prime}$ excitation pulses.
+    (b) Dependence of the output intensity on the $\tau_{22^\prime}$ and $\tau_{21}$ time delays
+    for $\omega_1=\omega_2$.
+    The solid lines define the regions for the six different time orderings of the $\omega_1$,
+    $\omega_2$, and $\omega_{2^\prime}$ excitation pulses.
+    We have developed a convention for numbering these time orderings, as shown.
+    (c) Diagram of the band structure of MoS\textsubscript{2} at the $K$ point.
+    The A and B exciton transitions are shown. (d) Two dimensional frequency-frequency plot
+    labeling two diagonal and cross-peak features for the A and B excitons.}
+  \label{fig:Czech01}
 \end{figure}
 
+\section{Methods}  % ==============================================================================
+
 MoS\textsubscript{2} thin films were prepared \textit{via} a Mo film sulfidation reaction, similar
 to methods reported by \textcite{LaskarMasihhurR2013a}.  %
 A 1 nm amount of Mo (Kurt J. Lesker, 99.95\%) metal was electron-beam evaporated onto a fused
@@ -194,15 +188,11 @@ TEM experiments were performed on a FEI Titan aberration corrected (S)TEM under
 voltage.  %
 
 \begin{figure}
-  \includegraphics[width=\textwidth]{MX2/10}
-  \caption[Mask and epi vs transmissive.]{
-    (a) Mask.
-    (b) 2D delay spectra at the BB diagonal ($\omega_1=\omega_2\approx1.95$ eV) for transmissive
-    and reflective geometries.
-    Transmissive signal is a mixture of MoS\textsubscript{2} signal and a large amount of driven
-    signal from the substrate that only appears in the pulse overlap region. Reflective signal is
-    representative of the pure MoS\textsubscript{2} response.}
-  \label{fig:Czech10}
+  \includegraphics[width=\textwidth]{MX2/S1}
+  \caption[Schematic of sulfidation reaction chamber.]{
+    Schemiatic of the synthetic setup used for Mo thin film sulfidation reactions.
+  }
+  \label{fig:CzechS1}
 \end{figure}
 
 The coherent multidimensional spectroscopy system used a 35 fs seed pulse, centered at 800 nm and
@@ -216,12 +206,6 @@ Signal and idler were not filtered out, but played no role due to their low phot
 Pulse $\omega_2$ was split into pulses labeled $\omega_2$ and $\omega_{2^\prime}$ to create a total
 of three excitation pulses.  %
 
-\begin{figure}
-  \includegraphics[width=\textwidth]{MX2/S4}
-  \caption{OPA outputs at each color explored.}
-  \label{fig:CzechS4}
-\end{figure}
-
 In this experiment we use motorized OPAs which allow us to set the output color in software.  %
 OPA1 and OPA2 were used to create the $\omega_1$ and $\omega_2$ frequencies, respectively.  %
 In \autoref{fig:CzechS4} we compare the spectral envelope generated by the OPA at each set
@@ -237,12 +221,6 @@ focused onto the sample surface by a 1 meter focal length spherical mirror in a
 geometry to form a 630, 580, and 580 $\mu$m FWHM spot sizes for $\omega_1$, $\omega_2$, and
 $\omega_{2^\prime}$, respectively.  %
 
-\begin{figure}
-  \includegraphics[width=\textwidth]{MX2/S5}
-  \caption{Spectral delay correction.}
-  \label{fig:CzechS5}
-\end{figure}
-
 \autoref{fig:CzechS5} represents delay corrections applied for each OPA.  %
 The corrections were experimentally determined using driven FWM output from fused silica.  %
 Corrections were approximately linear against photon energy, in agreement with the normal dispersion
@@ -251,6 +229,18 @@ OPA2 required a relatively small correction along $\tau_{22^\prime}$ (middle sub
 for any dispersion experienced differently between the two split beams.  %
 OPA1 was not split and therefore needed no such correction.  %
 
+\begin{figure}
+  \includegraphics[width=\textwidth]{MX2/S4}
+  \caption{OPA outputs at each color explored.}
+  \label{fig:CzechS4}
+\end{figure}
+
+\begin{figure}
+  \includegraphics[width=\textwidth]{MX2/S5}
+  \caption{Spectral delay correction.}
+  \label{fig:CzechS5}
+\end{figure}
+
 \autoref{fig:Czech10}a represents the to-scale mask that defines our distorted BOXCARS
 configuration.  %
 Relative to the center of the BOXCARS mask (small black dot), $\omega_1$, $\omega_2$, and
@@ -276,6 +266,38 @@ signal id is the geometry chosen for this experiment.  %
 This discrimination between a film and the substrate was also seen in reflective and transmissive
 CARS microscopy experiments. \cite{VolkmerAndreas2001a}  %
 
+
+\begin{figure}
+  \includegraphics[width=\textwidth]{MX2/10}
+  \caption[Mask and epi vs transmissive.]{
+    (a) Mask.
+    (b) 2D delay spectra at the BB diagonal ($\omega_1=\omega_2\approx1.95$ eV) for transmissive
+    and reflective geometries.
+    Transmissive signal is a mixture of MoS\textsubscript{2} signal and a large amount of driven
+    signal from the substrate that only appears in the pulse overlap region. Reflective signal is
+    representative of the pure MoS\textsubscript{2} response.}
+  \label{fig:Czech10}
+\end{figure}
+
+Once measured, the FWM signal was sent through a four-stage workup process to create the data set
+shown here.  %
+This workup procedure is visualized in \autoref{fig:Czech11}.  %
+We use a chopper and boxcar in active background subtraction mode (averaging 100 laser shots) to
+extract the FWM signal from $\omega_1$ and $\omega_2$ scatter.  %
+We collect this differential signal (\autoref{fig:Czech11}b) in software with an additional 50
+shots of averaging.  %
+In post-process we subtract $\omega_2$ scatter and smooth the data using a 2D Kaiser window.  %
+Finally, we represent the homodyne collected data as (sig)$^{1/2}$ to make the dynamics and line
+widths comparable to heterodyne-collected techniques like absorbance and pump-probe spectra.  %
+Throughout this work, zero signal on the color bar is set to agree with the average rather than the
+minimum of noise.  %
+Values below zero due to measurement uncertainty underflow the color bar and are plotted in
+white.  %
+This is especially evident in lots such as +120 fs in \autoref{fig:Czech08}, where there is no real
+signal.  %
+IPython \cite{PerezFernando2007a} and matplotlib \cite{HunterJohnD2007a} were important for data
+processing and plotting in this work.
+
 \begin{figure}
   \includegraphics[width=\textwidth]{MX2/11}
   \caption[MoS\textsubscript{2} post processing.]{
@@ -303,42 +325,8 @@ CARS microscopy experiments. \cite{VolkmerAndreas2001a}  %
   \label{fig:Czech11}
 \end{figure}
 
-Once measured, the FWM signal was sent through a four-stage workup process to create the data set
-shown here.  %
-This workup procedure is visualized in \autoref{fig:Czech11}.  %
-We use a chopper and boxcar in active background subtraction mode (averaging 100 laser shots) to
-extract the FWM signal from $\omega_1$ and $\omega_2$ scatter.  %
-We collect this differential signal (\autoref{fig:Czech11}b) in software with an additional 50
-shots of averaging.  %
-In post-process we subtract $\omega_2$ scatter and smooth the data using a 2D Kaiser window.  %
-Finally, we represent the homodyne collected data as (sig)$^{1/2}$ to make the dynamics and line
-widths comparable to heterodyne-collected techniques like absorbance and pump-probe spectra.  %
-Throughout this work, zero signal on the color bar is set to agree with the average rather than the
-minimum of noise.  %
-Values below zero due to measurement uncertainty underflow the color bar and are plotted in
-white.  %
-This is especially evident in lots such as +120 fs in \autoref{fig:Czech08}, where there is no real
-signal.  %
-IPython \cite{PerezFernando2007a} and matplotlib \cite{HunterJohnD2007a} were important for data
-processing and plotting in this work.
-
 \section{Results and discussion}  % ===============================================================
 
-\begin{figure}
-  \includegraphics[width=0.75\textwidth]{MX2/02}
-  \caption[Few-layer MoS\textsubscript{2} thin film characterization.]{
-    Characterization of the few-layer MoS\textsubscript{2} film studied in this work.
-    Optical images of the
-    MoS\textsubscript{2} thin film on fused silica substrate in (a) transmission and (b)
-    reflection.
-    (c) Raman spectrum of the $E_{2g}^1$ and $A_{1g}$ vibrational modes.
-    (d) High-resolution TEM image and its corresponding FFT shown in the inset.
-    (e) Absorption (blue), photoluminescence (green), Gaussian fits to the A and B excitons, along
-    with the residules betwen the fits and absorbance (dotted), A and B exciton centers (dotted)
-    and representative excitation pulse shape (red).}
-  \label{fig:Czech02}
-\end{figure}
-
 The few-layer MoS\textsubscript{2} thin film sample studied in this work was prepared on a
 transparent fused silica substrate by a simple sufidation reaction of a Mo thin film using a
 procedure modified from a recent report. \cite{LaskarMasihhurR2013a}  %
@@ -356,15 +344,6 @@ corresponds to approximately four monolayers.  %
 and B excitonic line shapes that were extracted from the absorption spectrum. A representative
 excitation pulse profile is also shown in red for comparison.  %
 
-\begin{figure}
-  \includegraphics[width=\textwidth]{MX2/S3}
-  \caption[MoS\textsubscript{2} absorbance.]{Extraction of excitonic features from absorbance
-    spectrum. (a) Second derivative spectra of absorbance (black) and fit second derivative
-    spectrum (green). Gaussian fit parameters are shown in the inset table. (b) Absorption curve
-    (black), Gaussian fits (blue and red), and remainder (black dotted).}
-  \label{fig:CzechS3}
-\end{figure}
-
 Extracting the exciton absorbance spectrum is complicated by the large ``rising background'' signal
 from other MoS\textsubscript{2} bands.  %
 With this in mind, we fit the second derivative absorption spectrum to a sum of two second
@@ -374,6 +353,30 @@ the amplitude) of the remainder between the fit and the absorption spectrum.  %
 The fit parameters can be found in the inset table in \autoref{fig:CzechS3}.  %
 The Gaussians themselves and the remainder can be found in \autoref{fig:CzechS3}.  %
 
+\begin{figure}
+  \includegraphics[width=0.75\textwidth]{MX2/02}
+  \caption[Few-layer MoS\textsubscript{2} thin film characterization.]{
+    Characterization of the few-layer MoS\textsubscript{2} film studied in this work.
+    Optical images of the
+    MoS\textsubscript{2} thin film on fused silica substrate in (a) transmission and (b)
+    reflection.
+    (c) Raman spectrum of the $E_{2g}^1$ and $A_{1g}$ vibrational modes.
+    (d) High-resolution TEM image and its corresponding FFT shown in the inset.
+    (e) Absorption (blue), photoluminescence (green), Gaussian fits to the A and B excitons, along
+    with the residules betwen the fits and absorbance (dotted), A and B exciton centers (dotted)
+    and representative excitation pulse shape (red).}
+  \label{fig:Czech02}
+\end{figure}
+
+\begin{figure}
+  \includegraphics[width=\textwidth]{MX2/S3}
+  \caption[MoS\textsubscript{2} absorbance.]{Extraction of excitonic features from absorbance
+    spectrum. (a) Second derivative spectra of absorbance (black) and fit second derivative
+    spectrum (green). Gaussian fit parameters are shown in the inset table. (b) Absorption curve
+    (black), Gaussian fits (blue and red), and remainder (black dotted).}
+  \label{fig:CzechS3}
+\end{figure}
+
 The multiresonant CMDS experiment uses $\approx$70 fs excitation pulses created by two
 independently tunable optical parametric amplifiers (OPAs).  %
 Automated delay stages and neutral density filters set the excitation time delays over all values
@@ -387,23 +390,9 @@ In order to compare the FWM spectra with the absorption spectrum, the signal has
 the square root of the measured FWM signal since FWM depends quadratically on the sample
 concentration and path length.  %
 
-\begin{figure}
-  \includegraphics[width=\textwidth]{MX2/03}
-  \caption[MoS\textsubscript{2} frequency-frequency slices.]{2D frequency-frequency spectra of the
-    MoS\textsubscript{2} sample in the epi configuration. In all spectra $\tau_{22^\prime}=0$ fs,
-    while $\tau_{21}$ is designated in the bottom-right corner of each spectral panel. The color
-    bar defines the square root of the intensity normalized to the most intense feature in the
-    series of spectra. The integration of the signal onto the $\hbar\omega_1=\hbar\omega_m$ and
-    $\hbar\omega_2$ axes are represented ans the blue curves in the top and right side plots,
-    respectively. The side plots also contain the absorbance spectrum (black line) to aid
-    intepretation of the dynamics of the integrated 2D signals. The dashed lines mark the centers
-    of the A and B excitons, as designated from the absorption spectrum.}
-  \label{fig:Czech03}
-\end{figure}
-
 The main set of data presented in this work is an $\omega_1\omega_2\tau_{21}$ ``movie'' with
 $\tau_{22\prime}=0$. 
-\autoref{fig:Czec03} shows representative 2D frequency-frequency slices from this movie at
+\autoref{fig:Czech03} shows representative 2D frequency-frequency slices from this movie at
 increasingly negative $\tau_{21}$ times.  %
 Each 2D frequency spectrum contains side plots along both axes that compare the absorbance spectrum
 (black) to the projection of the integrated signal onto the axis (blue).  %
@@ -414,24 +403,19 @@ Instead, the signal amplitude increases toward bluer $\omega_2$ values.  %
 The decrease in FWM above 2.05 eV is caused by a drop in the $\omega_2$ OPA power.
 
 \begin{figure}
-  \includegraphics[width=0.75\textwidth]{MX2/04}
-  \caption[MoS\textsubscript{2} $\omega_1$ Wigner progression.]{Mixed $\omega_1$---$\tau_{21}$
-    time---frequency representations of the 3D data set at five ascending $\omega_2$ excitation
-    frequencies (solid black lines) showing the impact of the $\omega_2$ excitation frequency on
-    the $\omega_1$ spectral line shape as a function of time. The A and B exciton energies are
-    marked as dashed lines within each spectrum.}
-  \label{fig:Czech04}
+  \includegraphics[width=\textwidth]{MX2/03}
+  \caption[MoS\textsubscript{2} frequency-frequency slices.]{2D frequency-frequency spectra of the
+    MoS\textsubscript{2} sample in the epi configuration. In all spectra $\tau_{22^\prime}=0$ fs,
+    while $\tau_{21}$ is designated in the bottom-right corner of each spectral panel. The color
+    bar defines the square root of the intensity normalized to the most intense feature in the
+    series of spectra. The integration of the signal onto the $\hbar\omega_1=\hbar\omega_m$ and
+    $\hbar\omega_2$ axes are represented ans the blue curves in the top and right side plots,
+    respectively. The side plots also contain the absorbance spectrum (black line) to aid
+    intepretation of the dynamics of the integrated 2D signals. The dashed lines mark the centers
+    of the A and B excitons, as designated from the absorption spectrum.}
+  \label{fig:Czech03}
 \end{figure}
 
-\begin{figure}
-  \includegraphics[width=0.75\textwidth]{MX2/05}
-  \caption[MoS\textsubscript{2} $\omega_2$ Wigner progression.]{Mixed $\omega_2$---$\tau_{21}$
-    time---frequency representations of the 3D data set at five ascending $\omega_1$ probe
-    frequencies (solid black lines) showing the impact of the $\omega_1$ excitation frequency on
-    the $\omega_2$ spectral line shape as a function of time. The A and B exciton energies are
-    marked as dashed lines within each spectrum.}
-  \label{fig:Czech05}
-\end{figure}
 
 Figures \ref{fig:Czech04} and \ref{fig:Czech05} show representative 2D frequency-delay slices from
 this movie, where the absicissa is the $\omega_1$ or $\omega_2$ frequency, respectively, the
@@ -441,9 +425,9 @@ The color bar is normalized to the brightest feature in each subplot.  %
 This normalization allows comparison of the time dependence of the line shapes, positions, and
 relative signal amplitudes along the $\omega_1$ or $\omega_2$ axis directly.  %
 
-Each subplot in \autoref{fig:Czech04} is similar to published pump-probe, transient absorption, multidimensionaland
-transient reflection experiments that have measured the electronc dynamics of the A and B excitons.
-\cite{XiaoDi2012a, FangHui2014a, KumarNardeep2013a, NieZhaogang2014a, SunDezheng2014a,
+Each subplot in \autoref{fig:Czech04} is similar to published pump-probe, transient absorption,
+and transient reflection experiments that have measured the electronc dynamics of the A and B
+excitons. \cite{XiaoDi2012a, FangHui2014a, KumarNardeep2013a, NieZhaogang2014a, SunDezheng2014a,
   SimSangwan2013a, MakKinFai2012a, ThomallaMarkus2006a}  %
 These previous experiments measure relaxation dynamics on the same $\approx$400-600 fs time scale
 that is characteristic of \autoref{fig:Czech04}.  %
@@ -459,15 +443,26 @@ Both the line shapes and the dynamics of the spectral features are very similar.
 features do not depend strongly on the $\omega_1$ frequency.
 
 \begin{figure}
-  \includegraphics[width=0.5\textwidth]{MX2/06}
-  \caption[Pathway V, VI liouville pathways.]{Liouville pathways for \autoref{fig:Czech04}. gg and
-    ee designate ground- and excited-state populations, the eg, 2e,e, and e$^\prime$+e,e represent
-    the excitonic and biexcitonic output coherences, and the arrows are labeled with the
-    frequencies or population transfer responsible for the transitions. e and e$^\prime$ represent
-    either A or B excitonic states.}
-  \label{fig:Czech06}
+  \includegraphics[width=0.75\textwidth]{MX2/04}
+  \caption[MoS\textsubscript{2} $\omega_1$ Wigner progression.]{Mixed $\omega_1$---$\tau_{21}$
+    time---frequency representations of the 3D data set at five ascending $\omega_2$ excitation
+    frequencies (solid black lines) showing the impact of the $\omega_2$ excitation frequency on
+    the $\omega_1$ spectral line shape as a function of time. The A and B exciton energies are
+    marked as dashed lines within each spectrum.}
+  \label{fig:Czech04}
+\end{figure}
+
+\begin{figure}
+  \includegraphics[width=0.75\textwidth]{MX2/05}
+  \caption[MoS\textsubscript{2} $\omega_2$ Wigner progression.]{Mixed $\omega_2$---$\tau_{21}$
+    time---frequency representations of the 3D data set at five ascending $\omega_1$ probe
+    frequencies (solid black lines) showing the impact of the $\omega_1$ excitation frequency on
+    the $\omega_2$ spectral line shape as a function of time. The A and B exciton energies are
+    marked as dashed lines within each spectrum.}
+  \label{fig:Czech05}
 \end{figure}
 
+
 The spectral features in Figures \ref{fig:Czech03}, \ref{fig:Czech04} and \ref{fig:Czech05} depend
 on the quantum mechanical interference effects caused by the different pathways.  %
 \autoref{fig:Czech06} shows all of the Liouville pathways required to understand the spectral
@@ -506,6 +501,18 @@ A quantitative treatment of the cancellation effects between the GSB, SE, and ES
 knowledge of the transition moments and state degeneracies and is beyond the scope of this paper.
 \cite{WongCathyY2011a}  %
 
+\begin{figure}
+  \includegraphics[width=0.5\textwidth]{MX2/06}
+  \caption[Liouville pathways in time orderings V, VI.]{
+    Liouville pathways for \autoref{fig:Czech04}. gg and
+    ee designate ground- and excited-state populations, the eg, 2e,e, and e$^\prime$+e,e represent
+    the excitonic and biexcitonic output coherences, and the arrows are labeled with the
+    frequencies or population transfer responsible for the transitions. e and e$^\prime$ represent
+    either A or B excitonic states.
+  }
+  \label{fig:Czech06}
+\end{figure}
+
 The most important characteristic of the experimental spectra is the contrast between the absence
 of well-resolved excitonic features that depend on $\omega_2$ in Figures \ref{fig:Czech03} and
 \ref{fig:Czech05} and the well-defined excitonic features that depend on $\omega_1$.  %
@@ -557,17 +564,6 @@ The B/A ratio is higher when $\omega_2$ is resonant with the B excitonic transit
 $\omega_2$ is lower than the A exciton frequency (the top subplot).  %
 If population transfer of holes from the B to A valence bands occurred during temporal overlap, the
 B/A ratio would be independent of pump frequency at $\tau_21<0$.
-
-\begin{figure}
-  \includegraphics[width=\textwidth]{MX2/07}
-  \caption[MoS\textsubscript{2} transients.]{
-    Transients taken at the different $\omega_1$ and $\omega_2$ frequencies indicated by the
-    colored markers on the 2D spectrum.
-    The dynamics are assigned to a 680 fs fast time constant (black solid line) and a slow time
-    constant represented as an unchanging offset over this timescale (black dashed line).}
-  \label{fig:Czech07}
-\end{figure}
-
 \autoref{fig:Czech07} shows the delay transients at the different frequencies shown in the 2D
 spectrum.  %
 The colors of the dots on the 2D frequency-frequency spectrum match the colors of the
@@ -583,11 +579,13 @@ The 680 fs decay is similar to previously published pump-probe and transient abs
 experiments. \cite{NieZhaogang2014a, SunDezheng2014a, DochertyCallumJ2014a}  %
 
 \begin{figure}
-  \includegraphics[width=\textwidth]{MX2/08}
-  \caption[MoS\textsubscript{2} frequency-frequency slices near pulse overlap.]{2D
-    frequency-frquency spectra near zero $\tau_{21}$ delay times. The signal amplitude is
-    normalized to the brightest features in each spectrum.}
-  \label{fig:Czech08}
+  \includegraphics[width=\textwidth]{MX2/07}
+  \caption[MoS\textsubscript{2} transients.]{
+    Transients taken at the different $\omega_1$ and $\omega_2$ frequencies indicated by the
+    colored markers on the 2D spectrum.
+    The dynamics are assigned to a 680 fs fast time constant (black solid line) and a slow time
+    constant represented as an unchanging offset over this timescale (black dashed line).}
+  \label{fig:Czech07}
 \end{figure}
 
 The spectral features change quantitatively for delay times near temporal overlap.  %
@@ -600,14 +598,6 @@ The spectra also develop more diagonal character as the delay time moves from ne
 values.  %
 The AB cross-peak is also a strong feature in the spectrum at early times.  %
 
-\begin{figure}
-  \includegraphics[width=0.5\textwidth]{MX2/09}
-  \caption[Pathways I, III Liouville pathways.]{Liouville pathways for the $\omega_1$, $\omega_2$,
-    and $\omega_{2^\prime}$ time ordering of pulse interactions. e and e$^\prime$ represent either
-    A or B excitonic states.}
-  \label{fig:Czech09}
-\end{figure}
-
 The pulse overlap region is complicated by the multiple Liouville pathways that must be
 considered.  %
 Additionally, interference between scattered light from the $\omega_1$ excitation beam and the
@@ -634,6 +624,23 @@ More positie values of $\tau_{21}$ emphasize the \autoref{fig:Czech09} pathways
 \autoref{fig:Czech06} pathways, accounting for the increasing percentage of diagonal character at
 increasingly positive delays.  %
 
+\begin{figure}
+  \includegraphics[width=\textwidth]{MX2/08}
+  \caption[MoS\textsubscript{2} frequency-frequency slices near pulse overlap.]{2D
+    frequency-frquency spectra near zero $\tau_{21}$ delay times. The signal amplitude is
+    normalized to the brightest features in each spectrum.}
+  \label{fig:Czech08}
+\end{figure}
+
+\begin{figure}
+  \includegraphics[width=0.5\textwidth]{MX2/09}
+  \caption[Liouville pathways in time orderings I, III.]{
+    Liouville pathways for the $\omega_1$, $\omega_2$,
+    and $\omega_{2^\prime}$ time ordering of pulse interactions. e and e$^\prime$ represent either
+    A or B excitonic states.}
+  \label{fig:Czech09}
+\end{figure}
+
 \section{Conclusions}  % ==========================================================================
 
 This paper presents the first coherent multidimensional spectroscopy of MoS\textsubscript{2} thin
@@ -650,7 +657,7 @@ These spectra are complementary to previous work because they allow a direct com
 initially excited excitonic states and the states creating the final output coherence.  %
 The spectra show that the same hot A and B exciton continuum states that are observed in the
 absorption spectrum also dominate the CMDS excitation spectra.  %
-They also show that rapid, <70 fs intraband relaxation occurs to create the band-edge A and B
+They also show that rapid, $<70$ fs intraband relaxation occurs to create the band-edge A and B
 excitonic features observed in the CMSD spectrum.  %
 The relative intensity of the diagonal peak features depends on the relative absorption strength of
 the A and B excitons.  %
@@ -669,4 +676,26 @@ complex MoS\textsubscript{2} and other TMDC heterostructures with quantum-state
 The frequency domain based multiresonant CMDS methods described in this paper will play a central
 role in these measurements.  %
 They use longer, independently tunable pulses that provide state-selective excitation over a wide
-spectral range without the requirement for interferometric stability.  %
\ No newline at end of file
+spectral range without the requirement for interferometric stability.  %
+
+\section{Errata}  % ===============================================================================
+
+The following is an errata for the of this chapter, published in November 2015.
+\cite{CzechKyleJonathan2015a}  %
+
+\begin{ditemize}
+  \item Reference 13 is identical to reference 9.
+  \item In the last paragraph of the introduction the sentence ``The experimental spectra differ
+    from the simple 2D spectrum shown in Figure 1d and those of earlier CMDS experiments with model
+    systems'' appears. This sentence cites references 6 through 10. Instead, it should cite
+    references 15 through 20.
+  \item In the last paragraph begining on page 12148, the text ``Automated delay stages and neutral
+    density filters set the excitation time delays over all values of $\tau_{21}$ with
+    $\tau_{22}=0$''    appears. For the second $\tau$, the subscript should read $22^\prime$, not
+    $22$.
+  \item Caption of Figure 5 reads, in part: "showing the impact of the $\omega_1$ excitation
+    frequency on the $\omega_1$ spectral line shape". This should instead read "showing the impact
+    of the $\omega_1$ excitation frequency on the $\omega_2$ spectral line shape". The subscript on
+    the last $\omega$ should be a 2 and not a 1.
+  \item Figure 6 e$^\prime$+e,e$^\prime$ should read e$^\prime$+e,e and vice versa.
+\end{ditemize}
\ No newline at end of file
diff --git a/PbSe_global_analysis/chapter.tex b/PbSe_global_analysis/chapter.tex
index 9805c7c..001f0ed 100644
--- a/PbSe_global_analysis/chapter.tex
+++ b/PbSe_global_analysis/chapter.tex
@@ -23,8 +23,8 @@ The fits of the combined dataset reveal and quantify the presence of continuum t
 \section{Introduction}  % =========================================================================
 
 Lead chalcogenide nanocrystals are among the simplest manifestations of quantum
-confinement \cite{Wise2000} and provide a foundation for the rational design of nano-engineered
-photovoltaic materials.  %
+confinement \cite{WiseFrankW2000a} and provide a foundation for the rational design of
+nano-engineered photovoltaic materials.  %
 The time and frequency resolution capabilities of the different types of ultrafast pump-probe
 methods have provided the most detailed understanding of quantum dot (QD) photophysics.  %
 Transient absorption (TA) studies have dominated the literature.  %
@@ -37,18 +37,18 @@ Although the real component is less important for photovoltaic performance, it i
 indicator of underlying structure and dynamics.  %
 In practice, having both real and imaginary components is often helpful.  
 For example, the fully-phased response is crucial for correctly interpreting spectroscopy when
-interfaces are important, which is common in evaluation of materials. \cite{Price2015, Yang2015,
-  Yang2017}  %
+interfaces are important, which is common in evaluation of materials. \cite{PriceMichaelB2015a,
+  YangYe2015a, YangYe2017a}  %
 The real and imaginary responses are directly related by the Kramers-Kronig relation, but it is
 experimentally difficult to measure the ultrafast response over the range of frequencies required
 for a Hilbert transform.  % 
 Interferometric methods, such as two-dimensional eletronic spectroscopy (2DES), can resolve both
-components, but they are demanding methods and not commonly used.  %
-% note that they often use TA to phase spectra
-
+components.  %
+In this work we report a strategy for resolving fully phased spectra using frequency domain
+coherent multidimensional spectroscopy.  %
 
 Transient grating (TG) is a pump-probe method closely related to TA.
-Figures \ref{fig:tg_vs_ta} illustrates both methods. 
+\autoref{psg:fig:tg_vs_ta} illustrates both methods. 
 In TG, two pulsed and independently tunable excitation fields, $E_1$ and $E_2$, are incident on a
 sample.  %
 The TG experiment modulates the optical properties of the sample by creating a population grating
@@ -60,14 +60,14 @@ modulates the ground and excited state populations with a chopper.  %
 TA can be seen as a special case of a TG experiment in which the grating fringes
 become infinitely spaced ($\vec{k}_2-\vec{k}_{2^\prime} \rightarrow \vec{0}$)
 and, instead of being diffracted, the nonlinear field overlaps and interferes with the probe
-beam.  %
+beam (self heterodyne).  %
 
 Like TA, TG does not fully characterize the non-linear response.  %
 Both imaginary and real parts of the refractive index spatially modulate in the TG experiment. 
 The diffracted probe is sensitive only to the total grating contrast (the response
 \textit{amplitude}), and not the phase relationships of the grating.  %
 Since both techniques are sensitive to different components of the non-linear response, however,
-the combination of both TA and TG can solve the fully-phased response.  %
+\emph{the combination of both TA and TG can solve the fully-phased response.}  %
 
 Here we report the results of dual 2DTA-2DTG experiments of PbSe quantum dots at the 1S exciton
 transition.  %
@@ -86,48 +86,42 @@ and dynamics.  %
 
 \begin{figure}
 	\includegraphics[width=\linewidth]{"PbSe_global_analysis/ta_vs_tg"}
-	\caption{The similarities between transient grating and transient absorption measurements. 
+	\caption[Similarities between transient grating and transient absorption measurements.]{
+    The similarities between transient grating and transient absorption measurements. 
 		Both signals are derived from creating a population difference in the sample. 
-			(a) A transient grating experiment crosses two pump beams of the same optical frequency
-      ($E_2$, $E_{2^\prime}$) to create an intensity grating roughly perpendicular to the direction
-      of propagation. 
-			(b) The intensity grating consequently spatially modulates the balance of ground state and
-      excited state in the sample.  
-				The probe beam ($E_1$) is diffracted, and the diffracted intensity is measured. 
-			In transient absorption (c), the probe creates a monolithic population difference, which
-      changes the attenuation the probe beam experiences through the sample. 
-			(d) The pump is modulated by a chopper, which facilitates measurement of the population
-      difference.}
+		The intensity grating consequently spatially modulates the balance of ground state and
+    excited state in the sample.  
+		The probe beam ($E_1$) is diffracted, and the diffracted intensity is measured. 
+		In transient absorption, the pump creates a monolithic population difference, which
+    changes the attenuation the probe beam experiences through the sample.
+    TA is sensitive to only the imaginary portion of the nonlinearity, while TG measures the total
+    (magnitude) response.  %
+  }
 	\label{psg:fig:tg_vs_ta}
 \end{figure}
 
 \section{Theory}  % ===============================================================================
 
-[FIGURE]
-
-The optical non-linearity of near-bandgap QD excitons has been extensively investigated. [CITE]  %
+The optical non-linearity of near-bandgap QD excitons has been extensively investigated.  %
 The response derives largely from state-filling and depends strongly on the exciton occupancy of
 the dots.  %
 In a PbSe quantum dot, the 1S peak is composed of transitions between 8 1S electrons and 8 1S
-holes. \cite{Kang1997}  %
-Figure \ref{fig:model_system}a shows the ground state configuration for a PbSe quantum dot.  %
+holes. \cite{KangInuk1997a}  %
+\autoref{psg:fig:model_system}a shows the ground state configuration for a PbSe quantum dot.  %
 The 8-fold degeneracy means there are $8 \times 8 = 64$ states in the single-exciton ($|1\rangle$)
 manifold and $7 \times 7 = 49$ states in the biexciton ($|2\rangle$) manifold, so $1/4$ of optical
 transitions are lost upon single exciton creation.  %
 
-Figure \ref{fig:model_system} shows the model system used in this study and the parameters that
+Figure \ref{psg:fig:model_system} shows the model system used in this study and the parameters that
 control the third-order response.  %
 
 We assign all electron-hole transitions the same dipole moment, $\mu_{eh}$, so that the total
 cross-section between manifolds, $N_i \mu_{eh}^2$, is determined by the number optically active
 transitions available, $N_i$.  %
-% BJT: state more correctly about what we are doing--there is the assumption that all dipoles are
-% the same, and there is the observable that cross-sections correspond to the number of optically
-% active transitions.
-Although this assumption has come under scrutiny \cite{Karki2013,Gdor2015} it remains valid for the
-perturbative fluence used in this study.  %
+Although this assumption has come under scrutiny \cite{KarkiKhadgaJ2013a, GdorItay2015b} it remains
+valid for the perturbative fluence used in this study.  %
 This model corresponds to the weakly interacting boson model, used to describe the four-wave mixing
-response of quantum wells,\cite{Svirko1999} in the limit of small quantum well area.  %
+response of quantum wells, \cite{SvirkoP1999a} in the limit of small quantum well area.  %
 
 With this excitonic structure, we now describe the resulting non-linear polarization.  %
 We restrict ourselves to field-matter interactions in which the pump, $E_2$, precedes the probe,
@@ -140,77 +134,93 @@ considered: $\text{Tr}\left[ \rho \right]  = (1-\bar{n})\rho_{00} + \bar{n} \rho
 $\bar{n}\ll 1$ is the (average) fractional conversion of population.  %
 In this limit, $\bar{n} = \frac{\sigma_0(\omega_2)}{\hbar \omega_2} \int{I_2(t) dt}$, where
 $I_2(t)$ is the pump intensity and $\sigma_0$ is the ground state absorptive cross-section.  %
-%In TA, $I_2(t) = \frac{nc}{8\pi}\left| E_2(t) \right|^2$, while in TG, the average pump intensity is $I_2(t) = \frac{nc}{4\pi}\left| E_2(t) \right|^2$
-%The pump induces a non-equilibrium population difference closely approximated by the Poisson distribution; in the case of low fluence ($\hbar \omega_2 \ll \sigma_{QD} \int{I_2 dt}$), 
-%The expected population conversion is $\langle n \rangle = \frac{\sigma_0(\omega_2)}{\hbar \omega_2} \int{I_2(t) dt}$, where $ I_2(t) = \frac{nc}{8\pi}\left| E_2(t) \right|^2 $ and $\sigma_0$ is the ground state absorptive cross-section.  
 For a Gaussian temporal profile of standard deviation $\Delta_t$, $I_2(t) = I_{\text{2,peak}} \exp
 \left( -t^2 / 	2\Delta_t^2 \right)$, the exciton population is
-\begin{equation}\label{psg:eq:n}
+\begin{equation} \label{psg:eqn:n}
 	\bar{n} = \frac{\sigma_0(\omega_2)\sqrt{2\pi}}{\hbar \omega_2} \Delta_t I_{\text{2,peak}}.
 \end{equation}
 When the probe interrogates this ensemble; each population will interact linearly:
-\begin{equation}\label{psg:eq:ptot}
+\begin{equation} \label{psg:eqn:ptot}
 \begin{split}
-	P_{\text{tot}} &= \left( 1-\bar{n} \right)\chi_0^{(1)}(\omega_1)E_1 + \bar{n} \chi_1^{(1)} (\omega_1)E_1 \\ 
+	P_{\text{tot}} &= \left( 1-\bar{n} \right)\chi_0^{(1)}(\omega_1)E_1 +
+  \bar{n} \chi_1^{(1)} (\omega_1)E_1 \\ 
 		&= \chi_0^{(1)}E_1 + \bar{n}\left( \chi_1^{(1)}(\omega_1) - \chi_0^{(1)}(\omega_1) \right)E_1.
 \end{split}
 \end{equation}
-Here $\chi_0^{(1)}$ ($\chi_1^{(1)}$) denotes the linear susceptibility of the pure state $|0\rangle$ ($|1\rangle$). The third-order field scales as $I_2 E_1$, so from Equation \ref{eq:ptot}
-\begin{equation}\label{psg:eq:chi3}
+Here $\chi_0^{(1)}$ ($\chi_1^{(1)}$) denotes the linear susceptibility of the pure state
+$|0\rangle$ ($|1\rangle$). The third-order field scales as $I_2 E_1$, so from
+\autoref{psg:eqn:ptot}
+\begin{equation} \label{psg:eqn:chi3}
 	\chi^{(3)} \propto \sigma_0 (\omega_2) \left( \chi_1^{(1)} - \chi_0^{(1)} \right).
 \end{equation}
-This expression accounts for the familiar population-level pathways such as excited state absorption/emission and ground state depletion. 
-Conforming the linear susceptibilities to our model, the non-linear portion of Equation \ref{eq:chi3} can be written as:
+This expression accounts for the familiar population-level pathways such as excited state
+absorption/emission and ground state depletion.  % 
+Conforming the linear susceptibilities to our model, the non-linear portion of
+\autoref{psg:eqn:chi3} can be written as:
 \begin{gather}
 	\chi^{(3)} \propto \mu_{eh}^4\text{Im}\left[L_0(\omega_2)\right]\left[ SL_1(\omega) -
-    L_0(\omega_1) \right], \label{psg:eq:chi3_lorentz} \\
+    L_0(\omega_1) \right], \label{psg:eqn:chi3_lorentz} \\
 	L_0(\omega) = \frac{1}{\omega - \omega_{10} + i\Gamma_{10}} ,\\
 	L_1(\omega) = \frac{1}{\omega - (\omega_{10} - \epsilon) + i\xi\Gamma_{10}} ,
 \end{gather}
 where $S = \frac{N_2}{N_1 + 1}$. The denominator of $S$ is larger than $N_1$ due to the
 contribution of stimulated emission; this contribution is often neglected.  % 
-From Equation \ref{eq:chi3_lorentz} we can see that a finite response can result from three
+From \autoref{psg:eqn:chi3_lorentz} we can see that a finite response can result from three
 conditions: $S\neq 1$, $\xi \neq 1$, and/or $\epsilon \neq 0$.  %
 The first inequality is the model's manifestation of state-filling, $S < 1$.  %
 If we assume that all 64 ground state transitions are optically active, then $S = 0.75$. 
 The second condition is met by exciton-induced dephasing (EID), $\xi > 1$, 
-% EID has also been attributed to stark splitting of exciton states
 and the third from the net attractive Coulombic coupling of excitons, $ \epsilon > 0 $. 
-The finite bandwidth of the monochromator can be accounted for by convolving equation
-\ref{eq:chi3_lorentz} with the monochromator instrumental function.  %
+The finite bandwidth of the monochromator can be accounted for by convolving
+\autoref{psg:eqn:chi3_lorentz} with the monochromator instrumental function.  %
 
-\subsection{The Bleach Nonlinearity}  % -----------------------------------------------------------
+\begin{figure}
+	\includegraphics[width=\linewidth]{"PbSe_global_analysis/model_system"}
+	\caption[Model system for the 1S band of PbSe quantum dots.]{
+		Model system for the 1S band of PbSe quantum dots. 
+		(a) The ground state shown in the electron-hole basis. 
+		All electrons (holes) are in the valence (conduction) band. 
+		There are two electrons and holes in each of the four degenerate $L$ points. 
+		(b) The excitonic basis and the transitions accessible in this experiment. 
+		The arrows illustrate the available absorptive or emissive transitions that take place in the
+    $\chi^{(3)}$ experiment, and are labeled by parameters that control the cross-sectional
+    strength (arrow width qualitatively indicates transition strength).
+	}
+	\label{psg:fig:model_system}
+\end{figure}
+
+\subsection{The bleach nonlinearity}  % -----------------------------------------------------------
 
 The bleach of the 1S band is the most studied nonlinear signature of the PbSe quantum dots. 
 Experiments keep track of the bleach factor, $\phi$, which is the proportionality factor that
 relates the relative change in the absorption coefficient at the exciton resonance,
 $\alpha_\text{FWM}(\omega_{10})$, with the average exciton occupation:  %
-\begin{equation}\label{psg:eq:bleach_factor}
+\begin{equation}\label{psg:eqn:bleach_factor}
 	\frac{\alpha_\text{FWM}(\omega_{10})}{\alpha_0(\omega_{10})} = -\phi \bar{n}
 \end{equation}
 where $\alpha_0$ is the linear absorption coefficient.  
 If QDs are completely bleached by the creation of a single exciton, then $\phi = 1$; if QDs are
 unperturbed by the exciton, then $\phi=0$.  %
 For PbSe, 1S-resonant values of $\phi=0.25$ and $0.125$ have been reported in literature
-\cite{Gdor2013a, Schaller2003, Nootz2011, Omari2012, Geiregat2014}, each with supporting theories
-on how state-filling should behave in an 8-fold degenerate system.  %
-Inspection of Equation \ref{eq:chi3_lorentz} shows that $\phi = \frac{\text{Im} \left[ L_0 - SL_1
+\cite{Gdor2013a, Schaller2003, NootzGero2011a, OmariAbdoulghafar2012a, GeiregatPieter2014a}, each
+with supporting theories on how state-filling should behave in an 8-fold degenerate system.  %
+Inspection of \autoref{psg:eqn:chi3_lorentz} shows that $\phi = \frac{\text{Im} \left[ L_0 - SL_1
   \right]}{L_0}$; if the 1S nonlinearity is dominated by state-filling ($\xi=1$ and $\epsilon=0$),
 then the bleach fraction has perfect correspondence with the change in the number of optically
 active states: $\phi = 1-S$.  %
 Because Coulombic shifts and EID act to decrease the resonant absorption of the excited state, we
 have the strict relation $\phi \geq 1-S$.  %
 
-More recently, a bleach factor metric has been adopted\cite{Trinh2008, Trinh2013} as the
+More recently, a bleach factor metric has been adopted \cite{Trinh2008, Trinh2013} as the
 proportionality between the spectrally integrated probe and the carrier concentration:  %
-\begin{equation}\label{psg:eq:bleach_factor_int}
+\begin{equation} \label{psg:eqn:bleach_factor_int}
 	\frac{\int{\alpha_\text{FWM}(\omega) d\omega}}{\int{\alpha_0(\omega)d\omega}} =
   -\phi_{\text{int}} \bar{n}.
 \end{equation}
 This metric is a more robust description of state filling, because it is unaffected by Coulomb
-shifts or EID: Equation \ref{eq:chi3_lorentz} gives $\phi_{\text{int}}=1-S$ regardless of $\xi$ and
+shifts or EID: \autoref{psg:eqn:chi3_lorentz} gives $\phi_{\text{int}}=1-S$ regardless of $\xi$ and
 $\epsilon$.  %
-An experimental value of $\phi_{\text{int}}=0.25$ has been reported\cite{Trinh2013} which
+An experimental value of $\phi_{\text{int}}=0.25$ has been reported \cite{Trinh2013} which
 consequently supports the measurement of $\phi = 0.25$.  %
 
 \subsection{TG/TA scaling}  % ---------------------------------------------------------------------
@@ -229,16 +239,15 @@ For TG, the polarization modulated in the phase-matched direction is given by
 	P_{\text{TG}}(z) = E_2^* E_{2^\prime} e^{-\alpha_0(\omega_2)z} E_1 e^{-\alpha_0(\omega_1)z / 2} \chi^{(3)}
 \end{equation}
 The TG electric field propagation can be solved using the slowly varying envelope approximation,
-which yields an output intensity of \cite{Carlson1989}
+which yields an output intensity of \cite{CarlsonRogerJohn1989a}
 \begin{gather}
 	I_{\text{TG}} \propto \omega_1^2 M_{\text{TG}}^2 I_1 I_2 I_{2^\prime} \left| \chi^{(3)} \right|^2 \ell^2 \\
 	M_{\text{TG}}(\omega_1, \omega_2) = e^{-\alpha_0(\omega_1)\ell / 2}
   \frac{1-e^{-\alpha_0(\omega_2)\ell}}{\alpha_0(\omega_2)\ell}.
 \end{gather}
-%$A_0(\omega)= \alpha_0 \ell \ln 10$ is the absorbance spectrum of the sample with path length $\ell$. 
 This motivates the following metric for TG:
 \begin{equation}
-\begin{split} \label{psg:eq:S_TG}
+\begin{split} \label{psg:eqn:S_TG}
 	S_{\text{TG}} &\equiv \frac{1}{M_{\text{TG}}\omega_1}\sqrt{\frac{N_{\text{TG}}}{I_1 I_2^2}} \\
 	&\propto \left| \chi^{(3)}\right|
 \end{split}
@@ -274,7 +283,7 @@ the change in the absorption $T = I_1 e^{-\alpha_{\text{tot}}\ell}$.  %
 We quantify the non-linear absorption by $\alpha_\text{FWM}=1/\ell \ln \frac{T-T_0}{T_0}= \alpha_0
 - \alpha_{\text{tot}}$, which can now be written as  %
 \begin{gather}
-	\alpha_\text{FWM}(\omega_1; \omega_2) = \frac{32\pi^2 \omega_1}{nc^2}M_{\text{TA}}(\omega_2) I_2 \text{Im} \left[\chi^{(3)}\right], \label{psg:eq:alpha_fwm} \\
+	\alpha_\text{FWM}(\omega_1; \omega_2) = \frac{32\pi^2 \omega_1}{nc^2}M_{\text{TA}}(\omega_2) I_2 \text{Im} \left[\chi^{(3)}\right], \label{psg:eqn:alpha_fwm} \\
 	M_{\text{TA}}(\omega_2) = \frac{1-e^{-\alpha_0(\omega_2)\ell}}{\alpha_0(\omega_2) \ell}.
 \end{gather}
 Similar to TG propagation, $\alpha_\text{FWM}$ suffers from absorption effects that distort the
@@ -285,8 +294,8 @@ automatically.  %
 Note that $M_{\text{TA}}$ simply accounts for the fraction of pump light absorbed in the sample,
 and consequently is closely related to the average exciton occupation across the entire path length
 of the sample, $\bar{n}_\ell = \ell^{-1} \int_0^\ell \bar{n}(z)dz$, which can be written using
-Equation \ref{eq:n} as:  %
-\begin{equation} \label{psg:eq:n_tot}
+\autoref{psg:eqn:n} as:  %
+\begin{equation} \label{psg:eqn:n_tot}
 	\bar{n}_\ell = \frac{\sqrt{2} \sigma_0 (\omega_2) \Delta_t}{\hbar \omega_2} I_{\text{2,peak}} M_{\text{TA}}
 \end{equation}
 We define an experimental metric that isolates the $\chi^{(3)}$ tensor:
@@ -300,43 +309,45 @@ For the general $\chi^{(3)}$ response, the imaginary and real components of the
 complicated relations owing to the causality of all three laser interactions.  %
 For the pump-probe time-ordered processes, the probe causality is separable from the pump
 excitation event, which makes the causality relation of the pump and probe separable.
-\cite{Hutchings1992}  %
+\cite{HutchingsDC1992a}  %
 The causality relations of the probe spectrum are the famous Kramers-Kronig equations that relate
 ground state absorption to the index of refraction.  %
-This relation is foundational in analysis of TG \cite{Hogemann1996} and TA measurements.
-% DK: need better citations for this
+This relation is foundational in analysis of TG \cite{HogemannClaudia1996a} and TA measurements.  %
 
 Theoretically, TA probe spectra alone could be transformed to generate the real spectrum. 
 In practice, such a transform is difficult because the spectral breadth needed to accurately
 calculate the integral is experimentally difficult to achieve.  % 
 When both the absolute value and the imaginary projection of $\chi^{(3)}$ are known, however, the
 real part can also be defined by the much simpler relation:  %
-\begin{equation} \label{psg:eq:chi_real}
+\begin{equation} \label{psg:eqn:chi_real}
 	\text{Re} \left[ \chi^{(3)} \right] =
   \pm \sqrt{\left| \chi^{(3)} \right|^2 - \text{Im} \left[ \chi^{(3)} \right]^2}
 \end{equation}
 % DK: concluding sentence
 
-\subsection{The Absorptive Third-Order Susceptibility}  % -----------------------------------------
+\subsection{The absorptive third-order susceptibility}  % -----------------------------------------
 
 Though the bleach factor is defined within the context of absorptive measurements, it can be
 converted into the form of a third-order susceptibility as well.  %
-Equations \ref{eq:ptot} and \ref{eq:bleach_factor} motivate alternative expressions for
+Equations \ref{psg:eqn:ptot} and \ref{psg:eqn:bleach_factor} motivate alternative expressions for
 differential absorptivity of the probe:  %
-\begin{equation}\label{psg:eq:alpha_fwm_to_bleach1}
+\begin{equation} \label{psg:eqn:alpha_fwm_to_bleach1}
 \begin{split}
 	\alpha_\text{FWM} &= \bar{n} \left( \alpha_1(\omega_1) - \alpha_0(\omega_1) \right) \\
 		& =-\phi \bar{n} \alpha_0(\omega_1).
 \end{split}
 \end{equation}
-Substituting Equation \ref{eq:n_tot} into \ref{eq:alpha_fwm_to_bleach1}, we can write the non-linear absorption as
+Substituting \autoref{psg:eqn:n_tot} into \autoref{psg:eqn:alpha_fwm_to_bleach1}, we can write the
+non-linear absorption as 
 \begin{equation}
-	\alpha_\text{FWM}(\omega_1) = -\phi \frac{\sqrt{2\pi}}{\hbar \omega_2} \alpha_0(\omega_1)\alpha_0(\omega_2) \Delta_t I_{\text{2,peak}} M_{\text{TA}}
+	\alpha_\text{FWM}(\omega_1) =
+  -\phi \frac{\sqrt{2\pi}}{\hbar \omega_2} \alpha_0(\omega_1)\alpha_0(\omega_2)
+  \Delta_t I_{\text{2,peak}} M_{\text{TA}}
 \end{equation}
-By Equation \ref{eq:alpha_fwm}, $\chi^{(3)}$ and the molecular hyperpolarizability, $\gamma_{\text{QD}}^{(3)} = \chi_{\text{QD}}^{(3)} / F^4 N_{\text{QD}}$, where $F$ is the local field factor, can also be related to the bleach fraction:
-\begin{gather}\label{psg:eq:chi3_state_filling}
+By \autoref{psg:eqn:alpha_fwm}, $\chi^{(3)}$ and the molecular hyperpolarizability, $\gamma_{\text{QD}}^{(3)} = \chi_{\text{QD}}^{(3)} / F^4 N_{\text{QD}}$, where $F$ is the local field factor, can also be related to the bleach fraction:
+\begin{gather} \label{psg:eqn:chi3_state_filling}
 	\text{Im}\left[ \chi^{(3)}_{\text{QD}} \right] = - \phi \frac{\sqrt{2\pi} n c^2 \Delta_t N_{\text{QD}}}{32\pi^2 \hbar \omega_1 \omega_2} \sigma_0(\omega_2)\sigma_0(\omega_1) \\
-	\text{Im}\left[ \gamma^{(3)}_{\text{QD}} \right] = - \phi \frac{\sqrt{2\pi} n c^2 \Delta_t }{32\pi^2 \hbar \omega_1 \omega_2 F^4} \sigma_0(\omega_2)\sigma_0(\omega_1). \label{psg:eq:gamma3_state_filling} 
+	\text{Im}\left[ \gamma^{(3)}_{\text{QD}} \right] = - \phi \frac{\sqrt{2\pi} n c^2 \Delta_t }{32\pi^2 \hbar \omega_1 \omega_2 F^4} \sigma_0(\omega_2)\sigma_0(\omega_1). \label{psg:eqn:gamma3_state_filling} 
 \end{gather}
 Because this formula only predicts the imaginary component of the signal, its magnitude gives an
 approximate lower limit for the peak susceptibility and hyperpolarizability.  %
@@ -346,17 +357,17 @@ Absorptive cross-sections have been experimentally determined for PbSe QDs.  %
 \section{Methods}  % ==============================================================================
 
 Quantum dot samples used in this study were synthesized using the hot injection
-method. \cite{Wehrenberg2002}  %
+method. \cite{WehrenbergBrianL2002a}  %
 Samples were kept in a glovebox after synthesis and exposure to visible and UV light was minimized.  
 These conditions preserved the dots for several months.  
 Two samples, Batch A and Batch B, are presented in this study, in an effort to show the robustness
 of the results.  %
-Properties of their optical characterization are shown in Table \ref{tab:QD_abs}.  
+Properties of their optical characterization are shown in \autoref{psg:tab:QD_abs}.  
 The 1S band of Batch A is broader than Batch B, an effect which is usually attributed to a wider
 size distribution and therefore greater inhomogeneous broadening.  %
 
 The experimental system for the TG experiment has been previously
-explained. \cite{Kohler2014, Czech2015}   %
+explained. \cite{KohlerDanielDavid2014a, CzechKyleJonathan2015a}   %
 Briefly, two independently tunable OPAs are used to make pulses $E_1$ and $E_2$ with colors
 $\omega_1$ and $\omega_2$.   %
 The third beam, $E_{2^\prime}$, is split off from $E_2$. The TG experiment utilized here uses
@@ -374,7 +385,7 @@ The signal, $N_{\text{TG}}$, was detected with an InSb photodiode. Reflective ne
 filters (Inconel) limit the pulse fluence to avoid multi-photon absorption.  % 
 To control for frequency-dependent changes in pulse arrival time due to the OPAs and the neutral
 density, a calibration table was established to assign a correct zero delay for each color
-combination (see supporting information for more details).  %
+combination.  %
 
 The TA experiments were designed to minimally change the TG experimental conditions.  %
 The $E_{2^\prime}$ beam was blocked and signal in the $\vec{k}_1$ direction was measured.  %
@@ -382,16 +393,8 @@ $E_2$ was chopped and the differential signal and the average signal were measur
 and $T$ needed to compute $\Delta A$.  %  
 Just as in TG experiments, the excitation frequencies were scanned while the monochromator was
 locked at $\omega_m=\omega_1$.  %
-% DK: perhaps leave this part out
-%Finally, fluence studies resonant with the 1S band were performed to test for indications of
-%intensity-dependent relaxation.   
-%These studies showed no indication of accelerated Auger recombination rates (see supporting info).
 
-\begin{table}[]
-	\centering
-	\caption{Batch Parameters extracted from absorption spectra.  $\langle d \rangle$: average QD
-    diameter, as inferred by the 1S transition energy.}
-	\label{psg:tab:QD_abs}
+\begin{table}
 	\begin{tabular}{l|cc}
 		 													  & A    & B \\
 		\hline
@@ -400,17 +403,22 @@ locked at $\omega_m=\omega_1$.  %
 		$ \langle d \rangle  \left(\text{nm}\right)$ 		  & 4    & 4.8 \\
 		$ \sigma_0 \left( \times 10^{16} \text{cm}^2 \right)$ & 1.7  & 2.9
 	\end{tabular}
+	\caption[Batch parameters extracted from absorption spectra.]{
+    Batch Parameters extracted from absorption spectra.
+    $\langle d \rangle$: average QD diameter, as inferred by the 1S transition energy.
+  }
+	\label{psg:tab:QD_abs}
 \end{table}
 
 \section{Results}  % ==============================================================================
 
-\subsection{Pump-Probe 3D acquisitions for TA and TG}  % ------------------------------------------
+\subsection{Pump-probe 3D acquisitions for TA and TG}  % ------------------------------------------
 
 For both samples, 2D spectra were collected for increments along the population rise time.  % 
 For these acquisitions, concentrated samples ($\text{OD}_{\text{1S}} \sim 0.6, 0.8$) were used to
 minimize contributions from non-resonant background.  % 
 Both samples maintained constant signal amplitude for at least hundreds of picoseconds after initial excitation, indicating multiexcitons and trapping were negligible effects in these studies. 
-The TA and TG results for both batches are shown in Figure \ref{fig:movies}. For $T<0$ (probe
+The TA and TG results for both batches are shown in \autoref{psg:fig:movies}. For $T<0$ (probe
 arrives before pump), both collections show spectral line-narrowing in the anti-diagonal
 direction.  % 
 This highly correlated line shape is indicative of an inhomogeneous distribution, but the
@@ -419,8 +427,8 @@ as the pump, the typical pump-probe pathways are suppressed and more unconventio
 probe-pump and pump-probe-pump pulse orderings are enhanced.  %
 Such pathways exhibit resonant enhancement when $\omega_1=\omega_2$, even in the absence of
 inhomogeneity.  %
-The pulse overlap effect is well-understood in both TA\cite{BritoCruz1988} and TG\cite{Kohler2017}
-experiments.  % 
+The pulse overlap effect is well-understood in both TA \cite{BritoCruzCH1988a} and TG
+\cite{KohlerDanielDavid2017a} experiments.  % 
 
 After the initial excitation rise time ($T > 50$ fs), the signal reaches a maximum, followed by a slight loss of signal ($\sim 10\%$) over the course of ~150 fs, after which the signal converges to a line shape that remains static over the dynamic range of our experiment ($200$ ps). 
 This signal loss occurs in both samples in both TA and TG; in TA measurements, the loss of
@@ -445,7 +453,7 @@ Our spectra show that the 2D line shape of the 1S exciton is significantly disto
 contributions from hot carrier excitation just above the 1S state.  %
 These hot carriers arise from transitions between the 1S and 1P resonances, which have been
 attributed to either the “rising edge” of the continuum or the pseudo-forbidden 1S-1P exciton
-transition \cite{Schins2009, Peterson2007}.  %
+transition \cite{SchinsJuleon2009a, PetersonJJ2007a}.  %
 Contributions from these hot carriers distort the 1S 2D line shape for $\omega_2 >
 \omega_{\text{1S}}$, resulting in a bleach feature centered at $\omega_1=\omega_{\text{1S}}$ and
 containing bleach contributions from the unresolved ensemble.  %
@@ -478,12 +486,14 @@ The difference in prominence shows that the redshifted feature is primarily refr
 character.  %
 
 \begin{figure}
-	\includegraphics[scale=0.5]{"PbSe_global_analysis/movies_combined"}
-	\caption{$S_{\text{TG}}$ (left) and $S_{\text{TA}}$ 2D spectra (see colorbar
-    labels) of Batch A (top) and Batch B (bottom) as a function of T delay. The
-    colors of each 2D spectrum are normalized to the global maximum of the 3D
-    acquisition, while the contour lines are normalized to each particular 2D
-    spectrum.}
+	\includegraphics[width=\textwidth]{"PbSe_global_analysis/movies_combined"}
+	\caption[TA and TG 2D spectra as function of delay.]{
+    $S_{\text{TG}}$ (left) and $S_{\text{TA}}$ 2D spectra (see colorbar labels) of Batch A (top)
+    and Batch B (bottom) as a function of T delay.
+    The colors of each 2D spectrum are normalized to the global maximum of the 3D acquisition,
+    while the contour lines are normalized to each particular 2D spectrum.
+    The vertical axis is pump energy, horizontal axis probe energy.  %
+  }
 	\label{psg:fig:movies}
 \end{figure}
 
@@ -492,25 +502,38 @@ character.  %
 \subsection{Comparison of TA and TG line shapes}  % -----------------------------------------------
 
 We first attempted simple fits on a subset of the data to reduce the parameter complexity. 
-We consider our experimental data at $T=120$ fs to remove contributions from probe-pump and pump-probe-pump time-ordered processes. 
-By further restricting our considerations to a single probe slice ($\omega_2 = \omega_{\text{1S}}$), we discriminate against inhomogeneous broadening and thus neglect ensemble effects for initial considerations. 
-We fit our probe spectrum with Equation \ref{eq:chi3_lorentz} along with the added treatment of convolving the response with our monochromator instrumental function. 
-Proper accounting of all time-ordered pathways, finite pulse bandwidth, and inhomogeneity are treated later.
+We consider our experimental data at $T=120$ fs to remove contributions from probe-pump and
+pump-probe-pump time-ordered processes.  %
+By further restricting our considerations to a single probe slice ($\omega_2 =
+\omega_{\text{1S}}$), we discriminate against inhomogeneous broadening and thus neglect ensemble
+effects for initial considerations.  %
+We fit our probe spectrum with \autoref{psg:eqn:chi3_lorentz} along with the added treatment of
+convolving the response with our monochromator instrumental function.   %
+Proper accounting of all time-ordered pathways, finite pulse bandwidth, and inhomogeneity are
+treated later.  %
+
+We find that the TA spectra are more sensitive to the model parameters than TG, and that the
+parameter interplay necessary to reproduce the spectra can be easily described.  %
+We note three features of the TA spectra that are crucial to reproduce in simulation: 1. the net
+bleach, 2. the photon energy of the bleach feature minimum is blue of the 1S absorption peak, 3.
+the probe spectrum decays quickly to zero away from the bleach resonance, leaving a minimal ESA
+feature to the red.  %
+These features are consistent with the vast majority of published TA spectra of the 1S exciton,
+\cite{Trinh2013, SchinsJuleon2009a, GesueleF2012a, Gdor2013a, KraatzIngvarT2014a,
+  DeGeyterBram2012a} and can only be reproduced when all three of our nonlinearities
+(state-filling, excitation-induced broadening, and Coulombic coupling) are invoked (exposition of
+this result is found in supporting information).  %
+The extracted fit parameters are listed in \autoref{psg:tab:fit1}.
 
 \begin{figure}
 	\includegraphics[scale=0.5]{"PbSe_global_analysis/kramers_kronig"}
-	\caption{Kramers-Kronig analysis of TA spectra compared with TG spectra.}
+	\caption[Kramers-Kronig analysis of TA, TG spectra.]{
+    Kramers-Kronig analysis of TA spectra compared with TG spectra.
+  }
 	\label{psg:fig:kramers_kronig}
 \end{figure}
 
-We find that the TA spectra are more sensitive to the model parameters than TG, and that the parameter interplay necessary to reproduce the spectra can be easily described. 
-We note three features of the TA spectra that are crucial to reproduce in simulation: (1) the net bleach; (2) the photon energy of the bleach feature minimum is blue of the 1S absorption peak; (3) the probe spectrum decays quickly to zero away from the bleach resonance, leaving a minimal ESA feature to the red. 
-These features are consistent with the vast majority of published TA spectra of the 1S exciton,\cite{Trinh2013,Schins2009,Gesuele2012,Gdor2013a,Kraatz2014,DeGeyter2012} and can only be reproduced when all three of our nonlinearities (state-filling, excitation-induced broadening, and Coulombic coupling) are invoked (exposition of this result is found in supporting information). 
-The extracted fit parameters are listed in Table \ref{tab:fit1}.
-
 \begin{table}[]
-	\centering
-	\caption{Parameters used in fitting experimental probe slices using Equation \ref{eq:chi3_lorentz}; $S=0.75$, $\omega_2 = \omega_\text{1S}$.}
 	\label{psg:tab:fit1}
 	\begin{tabular}{l|cc}
 		& \multicolumn{2}{l}{Batch} \\
@@ -520,69 +543,88 @@ The extracted fit parameters are listed in Table \ref{tab:fit1}.
 		$ \Gamma_{10} \left(\text{cm}^{-1}\right)$ 			   & 380  & 200         \\
 		$ \xi $ 								   			   & 1.35 & 1.39        
 	\end{tabular}
+	\caption[Parameters used in fitting probe slices.]{
+    Parameters used in fitting experimental probe slices using \autoref{psg:eqn:chi3_lorentz};
+    $S=0.75$, $\omega_2 = \omega_\text{1S}$.
+  }
 \end{table}
 
-With the 1S band TA well-characterized, we now consider applying the fitted parameters to the TG signal.
+With the 1S band TA well-characterized, we now consider applying the fitted parameters to the TG
+signal.  %
 Transferring this simulation to the TG data poses technical challenges. 
 A critical factor is appropriately scaling the TG signals relative to TA signals.
-The TA simulation can define the scaling through the Kramers-Kronig transform, which gives the transient refraction. 
-The computed transient refraction is unique to within an arbitrary offset; for a single resonant TA/TG feature, the transient refraction offset is zero.  
+The TA simulation can define the scaling through the Kramers-Kronig transform, which gives the
+transient refraction.  %
+The computed transient refraction is unique to within an arbitrary offset; for a single resonant
+TA/TG feature, the transient refraction offset is zero.  %  
 We take the offset to be zero now and address this assumption later.
-The transient refraction (Figure \ref{fig:cw_sim1}, third column) shows highly dispersive character with a node near resonance. 
+The transient refraction shows highly dispersive character with a node near resonance.   %
 This means that there is a point in our spectrum at which $\left| \chi^{(3)} \right| = \left| \text{Im}\left[ \chi^{(3)} \right] \right|$.  
 Of course, we also have the constraint $ \left| \chi^{(3)} \right| \geq \left| \text{Im} \left[ \chi^{(3)} \right]  \right|$ for every probe color. 
 These two constraints uniquely determine the appropriate scaling factor as the minimum scalar $c_0$ that satisfies $c_0 S_{\text{TG}} \geq \left| S_{\text{TA}} \right|$ for all probe colors. 
-%Such a scaling of the experimental data is consistent with our TA fit because the peak TA component is nearly equal to the peak TG amplitude (when the arbitrary offset of the KK-transform is zero).
 
 As we alluded, the arbitrary offset of the Kramers-Kronig transform deserves special consideration.  
-%A single TA resonance should not cause an offset in the transient reflection spectrum, but it is conceivable that states outside our spectral range are strongly coupled to the 1S band and produce strong refractive signals at these colors.  
-%While the peaked TA line shape might seem to imply a dispersively shaped $\text{Re} \left[ \chi^{(3)} \right]$ with a node near the bleach center, this is not guaranteed by the Kramers-Kronig relations.  
 The physical origin for this offset would be coupling between the 1S band and states outside our spectral range. 
 If the coupling is sufficiently strong, the $\text{Re}\left[ \chi^{(3)} \right] $ offset may be large enough to remove the node, invalidating the minimum scaling factor method.  
 We believe such a large offset is not viable for several reasons.
 From a physical standpoint, it seems very unlikely a non-resonant state would have coupling stronger coupling to the 1S band than the 1S band itself.  % DK: elaborate?
-Furthermore, TA studies with larger spectral ranges see no evidence of coupling features with strength comparable to the 1S bleach.\cite{Gdor2013a,Trinh2013}  
-Finally, a sufficiently large refractive offset would produce a dip in the $\left| \chi^{(3)} \right|$ line shape near the FHWM points; such features are definitively absent in the TG spectra.  
-
-% DK:  a more direct topic sentence might be nice--jump to the fact that equation 4 is inadequate
-%While we have confidence in relating the TG and TA measurements using the minimum scaling factor, we find two inconsistencies when fitting Equation \ref{eq:chi3_lorentz} with experiment. 
-%Firstly, the TA simulation gives a resonant bleach factor that is much greater than that predicted by state-filling alone: for instance, with $1-S=0.25$, we see $\phi>0.5$ for both batches (Figure \ref{fig:cw_sim1}, first column). 
-%While our parameters successfully recreate the features of the TA line shapes, the simulation grossly overshoots the magnitude of the non-linearity. 
-While we have confidence in relating the TG and TA measurements using the minimum scaling factor, Equation \ref{eq:chi3_lorentz} fails to accurately reproduce the TG spectrum (Figure \ref{fig:cw_sim1}, third column). 
-The errors are systematic: in both batches, our simulation misses the characteristic red skew of our experimental TG and instead skews signal to the blue. 
-Based on the excellent agreement with $S_{\text{TA}}$ (Figure \ref{fig:cw_sim1}, second column), it follows that the chief source of error in our simulation lies in $\text{Re} \left[ \chi^{(3)} \right]$.  
-The experimentally-consistent $\text{Re} \left[ \chi^{(3)} \right]$ can be calculated from Equation \ref{eq:chi_real} (Figure \ref{fig:cw_sim1}, fourth column). 
-The dark green curve highlights which of the two roots of Equation \ref{eq:chi_real} is closest to our simulation.
-The choice of a negative non-linear refractive index value on resonance is in agreement with z-scan measurements.\cite{Moreels2006} 
-The discrepancy between the experimental and simulated real components is well-approximated by a constant offset.
-
-\begin{figure}
-	\includegraphics[width=\textwidth]{"CW_sim2"}
-	\caption{Top row: Global fits of $S_\text{TA}$ (blue), $S_\text{TG}$ (red), and the associated
-    real projection (green) using Equation \ref{eq:offset_fit}.
-    Light colors indicate the simulations and the darker lines indicate the experimental data.
-    Bottom row:  Final simulated absorption spectra for the excited state and the ground state.
-  }
-	\label{psg:fig:cw_sim2}
-\end{figure}
+Furthermore, TA studies with larger spectral ranges see no evidence of coupling features with
+strength comparable to the 1S bleach. \cite{Gdor2013a, Trinh2013}  
+Finally, a sufficiently large refractive offset would produce a dip in the $\left| \chi^{(3)}
+\right|$ line shape near the FHWM points; such features are definitively absent in the TG
+spectra.  %
+
+While we have confidence in relating the TG and TA measurements using the minimum scaling factor,
+\autoref{psg:eqn:chi3_lorentz} fails to accurately reproduce the TG spectrum.  %
+The errors are systematic: in both batches, our simulation misses the characteristic red skew of
+our experimental TG and instead skews signal to the blue.  %
+Based on the excellent agreement with $S_{\text{TA}}$, it follows that the chief source of error in
+our simulation lies in $\text{Re} \left[ \chi^{(3)} \right]$.  %
+The experimentally-consistent $\text{Re} \left[ \chi^{(3)} \right]$ can be calculated from
+\autoref{psg:eqn:chi_real}.  %
+The dark green curve highlights which of the two roots of \autoref{psg:eqn:chi_real} is closest to
+our simulation.  %
+The choice of a negative non-linear refractive index value on resonance is in agreement with z-scan
+measurements. \cite{Moreels2006} 
+The discrepancy between the experimental and simulated real components is well-approximated by a
+constant offset.  %
 
 The presence of this offset forced a re-evaluation of the model.
-By revising our simulation to include a complex-valued offset term, $\Delta = |\Delta| e^{i\theta}$, so that
-\begin{equation} \label{psg:eq:offset_fit}
-	\chi^{(3)} = \text{Im} \left[ L_0(\omega_2) \right] \left[ SL_1(\omega_1) - L_0(\omega_1) + \Delta \right],
+By revising our simulation to include a complex-valued offset term, $\Delta = |\Delta|
+e^{i\theta}$, so that
+\begin{equation} \label{psg:eqn:offset_fit}
+	\chi^{(3)} = \text{Im} \left[ L_0(\omega_2) \right] \left[ SL_1(\omega_1) -
+    L_0(\omega_1) + \Delta \right],
 \end{equation}
 the discrepancy between $S_{\text{TA}}$ and $S_{\text{TG}}$ can be resolved. 
-It was found, however, that minimizing error between Equation \ref{eq:offset_fit} and the two datasets alone does not confine all variables uniquely. 
-Specifically, the effects of the EID non-linearity ($\xi$) on the line shapes were highly correlated with those of the broadband photoinduced absorption (PA) term, $\text{Im}\left[ \Delta \right]$, so that getting a unique parameter combination was not possible. 
+It was found, however, that minimizing error between \autoref{psg:eqn:offset_fit} and the two datasets
+alone does not confine all variables uniquely.  %
+Specifically, the effects of the EID non-linearity ($\xi$) on the line shapes were highly
+correlated with those of the broadband photoinduced absorption (PA) term, $\text{Im}\left[ \Delta
+\right]$, so that getting a unique parameter combination was not possible.  % 
 The fitting routine was robust, however, when the resonant bleach magnitude was pinned to the
 state-filling: $\phi \approx 1-S$.  %
 Robustness here is defined as the ability to permute the fitting parameter order when minimizing
 the residual. For example, $\tau_{10}$ can be fit either before or after $\Delta$ is fit without
 significantly changing the resulting parameters.  %
-The resulting parameters are shown in Table \ref{tab:fit2}, and the results of the fit are shown in Figure \ref{fig:cw_sim2}. 
+The resulting parameters are shown in \autoref{psg:tab:fit2}, and the results of the fit are shown
+in \autoref{psg:fig:cw_sim2}.  %
 As both $\phi=0.25$ and $\phi_{\text{int}}=0.25$ have been measured, this added constraint has a reasonable precedence. 
-As mentioned earlier, EID and Coulombic coupling prevent this equality (as in equation \ref{eq:chi3_lorentz}), but PA can counteract these terms and keep the resonant bleach near $1-S$. 
-In addition, since EID and PA are correlated, decreasing the resonant bleach also decreases the need for EID to influence the line shape, which results in a reduced EID non-linearity in this fit (compare $\xi$ in Table \ref{tab:fit1} and Table \ref{tab:fit2}).
+As mentioned earlier, EID and Coulombic coupling prevent this equality (as in
+\autoref{psg:eqn:chi3_lorentz}), but PA can counteract these terms and keep the resonant bleach
+near $1-S$.  %
+In addition, since EID and PA are correlated, decreasing the resonant bleach also decreases the
+need for EID to influence the line shape, which results in a reduced EID non-linearity in this fit
+(compare $\xi$ in \autoref{psg:tab:fit1} and \autoref{psg:tab:fit2}).  %
+
+\begin{figure}
+	\includegraphics[width=0.5\textwidth]{"PbSe_global_analysis/driven_complete"}
+	\caption[Fits of probe traces.]{
+    Top row:  Final simulated absorption spectra for the excited state and the ground state.
+    Bottom row: magnitude (purple), imagninary (blue) and real (red) projections and fits.  %
+  }
+	\label{psg:fig:cw_sim2}
+\end{figure}
 
 \begin{table}[]
 	\begin{tabular}{l|cc}
@@ -594,86 +636,110 @@ In addition, since EID and PA are correlated, decreasing the resonant bleach als
 		$ \left|\Delta \right| / \text{Im}\left[ L_0(\omega_\text{1S}) \right] $ &  0.07 (0.06) & 0.06 (0.06)   \\
 		$ \theta \left( \text{deg} \right)$ 					 			     & 151 (156)    & 146 (148)
 	\end{tabular}
-  \caption{
+  \caption[Parameters extracted by global fits.]{
     Parameters of the simulated $\chi^{(3)}$ response extracted by global fits of TA and TG
-    at $T=120$ fs using Equation \ref{eq:offset_fit} and with $S=0.75$.
+    at $T=120$ fs using Equation \ref{psg:eqn:offset_fit} and with $S=0.75$.
     Numbers in parentheses refer to fits at $T=300$ fs.  %
   }
 	\label{psg:tab:fit2}
 \end{table}
 
-The nature of $\Delta$ may be associated with resonant and non-resonant transitions from the 1S excitonic state to the continuum of intraband states involving the electron and/or hole. 
-The magnitude and phase of this contribution would then depend on the ensemble average from all transitions. 
+The nature of $\Delta$ may be associated with resonant and non-resonant transitions from the 1S
+excitonic state to the continuum of intraband states involving the electron and/or hole.  % 
+The magnitude and phase of this contribution would then depend on the ensemble average from all
+transitions.  %
 This contribution has been identified in previous TA studies.
-DeGeyter et.al. isolated a net absorption at sub-bandgap probe
-frequencies.\cite{DeGeyter2012}
-Geigerat et.al. found an absorptive contribution was needed to explain the fluence dependence of the 1S-resonant bleach.\cite{Geiregat2014}  
-The absorptive contribution was $5\%$ of the 1S absorbance, a value that agrees with this work (see Table \ref{tab:fit2}). 
-Our data unifies both observations by showing that additional contribution persists at both bandgap and sub-bandgap frequencies. 
+\textcite{DeGeyterBram2012a} isolated a net absorption at sub-bandgap probe frequencies.  %
+\textcite{GeiregatPieter2014a} found an absorptive contribution was needed to explain the fluence
+dependence of the 1S-resonant bleach.  %
+The absorptive contribution was $5\%$ of the 1S absorbance, a value that agrees with this work (see
+\autoref{psg:tab:fit2}). 
+Our data unifies both observations by showing that additional contribution persists at both bandgap
+and sub-bandgap frequencies.  %
 In addition, our data provides the spectral phase of the contribution. 
-It also shows that the red skew of the TG line shape is very sensitive to the relative importance of the 1S resonance and the additional contribution.
-
-There may also be a relationship to studies on CdSe quantum dots where an additional broad ESA feature was observed for $\omega_1 < \omega_{\text{1S}}$. 
-The feature had separate narrow and broad components.
-The narrow component closest to the band edge bleach corresponded to the Coulombically shifted biexciton transition. 
-Since the broad component correlated with inadequate surface passivation, it was attributed to the surface inducing ESA transitions to the broad band of continuum states that would normally be forbidden. 
-In addition to creating additional ESA transitions, it also created a short-lived transient that was similar to the transients attributed to multiexciton relaxation and multiexcion generation.
-
-\subsection{Determination of State Filling Factor}  % ---------------------------------------------
-
-% Given the 
-%We measured the peak $chi^{(3)}$ hyperpolarizability of Batch B via standard additions (SI) and found good agreement with the TA hyperpolarizability from Equation \ref{eq:gamma3_state_filling} (for $\phi=0.25$).
-%A offset that would remove the node would create a difference between the peak $\chi^{(3)}$ values measured from TA and TG.  
-%To check this possibility we measured the absolute susceptibility of the TG response and compared it to the susceptibility due to TA. 
-Our results show that the peak susceptibility is almost entirely imaginary, which means we can attribute our TG peak $\left|\chi^{(3)}\right|$ susceptibility to the TA state-filling factor.  
-A standard addition method was used to extract the peak TG hyperpolarizability of $\left| \gamma^{(3)} \right|=7\pm 2 \cdot 10^{-31} \text{cm}^6 \ \text{erg}^{-1}$ for Batch A, while Equation \ref{eq:gamma3_state_filling} predicts a peak imaginary hyperpolarizability of $\text{Im}\left[ \gamma^{(3)} \right] =2.1\phi \cdot 10^{-30} \text{cm}^6 \ \text{erg}^{-1}$. 
-The two values are similar for $\phi=0.25$, while $\phi=0.125$ gives an imaginary component that is much smaller than the total susceptibility.  
-We conclude that only the $\phi=0.25$ bleach factor is consistent with our observations. % DK: also discuss what S likely is
-
-\subsection{Inhomogeneity and the Pulse Overlap Response}  % --------------------------------------
-
-Our parameter extraction above gives plausible parameters to explain the observed photophysics of a small slice of our multidimensional data. 
+It also shows that the red skew of the TG line shape is very sensitive to the relative importance
+of the 1S resonance and the additional contribution.  %
+
+There may also be a relationship to studies on CdSe quantum dots where an additional broad ESA
+feature was observed for $\omega_1 < \omega_{\text{1S}}$.  % 
+The feature had separate narrow and broad components.  %
+The narrow component closest to the band edge bleach corresponded to the Coulombically shifted
+biexciton transition.  %
+Since the broad component correlated with inadequate surface passivation, it was attributed to the
+surface inducing ESA transitions to the broad band of continuum states that would normally be
+forbidden.  % 
+In addition to creating additional ESA transitions, it also created a short-lived transient that
+was similar to the transients attributed to multiexciton relaxation and multiexcion generation.  %
+
+\subsection{Determination of state filling factor}  % ---------------------------------------------
+
+Our results show that the peak susceptibility is almost entirely imaginary, which means we can
+attribute our TG peak $\left|\chi^{(3)}\right|$ susceptibility to the TA state-filling factor.  %
+A standard addition method was used to extract the peak TG hyperpolarizability of $\left|
+  \gamma^{(3)} \right|=7\pm 2 \cdot 10^{-31} \text{cm}^6 \ \text{erg}^{-1}$ for Batch A, while
+Equation \ref{psg:eqn:gamma3_state_filling} predicts a peak imaginary hyperpolarizability of
+$\text{Im}\left[ \gamma^{(3)} \right] =2.1\phi \cdot 10^{-30} \text{cm}^6 \ \text{erg}^{-1}$.  %
+The two values are similar for $\phi=0.25$, while $\phi=0.125$ gives an imaginary component that is
+much smaller than the total susceptibility.  %
+We conclude that only the $\phi=0.25$ bleach factor is consistent with our observations.  %
+
+\subsection{Inhomogeneity and the pulse overlap response}  % --------------------------------------
+
+Our parameter extraction above gives plausible parameters to explain the observed photophysics of a
+small slice of our multidimensional data.  %
 We now apply a more rigorous simulation of the model system to address the entire dataset and consider the broader experimental space. 
 This rigorous simulation is meant to account for the complex signals that arise at temporal pulse overlap, the pulsed nature of our excitation beams, and sample inhomogeneity. 
-We calculate signal through numerical integration techniques. % DK: cite paper 1
+We calculate signal through numerical integration techniques. \cite{KohlerDanielDavid2017a}  %
 The homogeneous and inhomogeneous broadening were constrained to compensate each other so that the
 total ensemble line shape was kept constant and equal to that extracted from absorption
-measurements (Table \ref{tab:QD_abs}).  %
+measurements (\autoref{psg:tab:QD_abs}).  %
 For a Lorentzian of FWHM $2\Gamma_{10}$ and a Gaussian line shape of standard deviation
 $\sigma_{\text{inhom}}$, the resulting Voigt line shape has a FWHM well-approximated by
 $\text{FWHM}_{\text{tot}} \left[ \text{cm}^{-1} \right] \approx 5672 \Gamma_{10}\left[
   \text{fs}^{-1} \right] + \sqrt{2298 \Gamma_{10}\left[ \text{fs}^{-1} \right] + 8 \ln 2
-  \left(\sigma_{\text{inhom}}\left[ \text{cm}^{-1} \right]\right) }$. \cite{Olivero1977}  %
+  \left(\sigma_{\text{inhom}}\left[ \text{cm}^{-1} \right]\right) }$. \cite{OliveroJJ1977a}  %
+
+Rather than a complete global fit of all parameters, we accept the non-linear parameters extracted
+earlier (\autoref{psg:tab:fit2}) and optimize the inhomogeneous-homogeneous ratio only, using the
+ellipticity of the 2D peak shape \cite{OkumuraKo1999a} at late population times as the figure of
+merit.  %
+\autoref{psg:tab:fit3} shows the degree of inhomogeneity in the sample from matching the ellipticity
+of the peak shape.  %
+As expected, Batch B fits to a smaller Gaussian distribution of transition energies than Batch A,
+but it also exhibits a longer homogeneous dephasing time, suggesting different system-bath coupling
+strengths for both samples.  %
+Previous research on the coherent dynamics of PbSe 1S exciton feature noted homogeneous lifetimes
+are a significant source of broadening on the 1S exciton; \cite{KohlerDanielDavid2014a} our results
+demonstrate that the relationship between exciton size distribution and 1S exciton linewidth is
+further complicated by sample-dependent system-bath coupling.  %
+
+The results of this final simulation are compared with the experimental data in Figure
+\autoref{psg:fig:nise_fits}.  %
+It is important to note that the simulations get many details of the rise-time spectra correct. 
+Particularly, the transition from narrow, diagonal line shape to a broad, less diagonal line shape
+is reproduced very well in both TA and TG simulations.  %
+Such behavior is expected for responses from excitonic peaks of material systems; the rise time
+behavior for such systems was studied in detail previously. \cite{KohlerDanielDavid2017a}  % 
+Because these simulations do not account for hot-exciton creation from the pump, simulations differ
+from experiment increasingly as the pump becomes bluer than the 1S center.  %
 
 \begin{table}[]
-	\centering
-	\caption{}
-	\label{psg:tab:fit3}
 	\begin{tabular}{l|cc}
 		Batch & A          & B         \\
 		\hline
 		$ \Gamma_{10} \left( \text{cm}^{-1}  \right) $ 					& 220 & 130    \\
 		$\text{FWHM}_\text{inhom} \left( \text{cm}^{-1} \right)$		& 520 & 360
 	\end{tabular}
+  \caption[Homogeneous and inhomogeneous linewidths.]{
+    Homogeneous and inhomogeneous linewidths extracted by global analysis using numerical
+    integration.  %
+  }
+	\label{psg:tab:fit3}
 \end{table}
 
-Rather than a complete global fit of all parameters, we accept the non-linear parameters extracted
-earlier (Table \ref{tab:fit2}) and optimize the inhomogeneous-homogeneous ratio only, using the
-ellipticity of the 2D peak shape\cite{Okumura1999} at late population times as the figure of
-merit.  %
-Table \ref{tab:fit3} shows the degree of inhomogeneity in the sample from matching the ellipticity
-of the peak shape.  %
-As expected, Batch B fits to a smaller Gaussian distribution of transition energies than Batch A,
-but it also exhibits a longer homogeneous dephasing time, suggesting different system-bath coupling
-strengths for both samples.  %
-Previous research on the coherent dynamics of PbSe 1S exciton feature noted homogeneous lifetimes
-are a significant source of broadening on the 1S exciton;\cite{Kohler2014} our results demonstrate
-that the relationship between exciton size distribution and 1S exciton linewidth is further
-complicated by sample-dependent system-bath coupling.  %
-
 \begin{figure}
 	\includegraphics[width=\linewidth]{"PbSe_global_analysis/movies_fitted"}
-	\caption{
+	\caption[Global simulation.]{
 		Global simulation using numerical integration and comparison with experiment. 
 		Batches A (left block) and B (right block) are shown, with the TG experimental (top), the
     simulated TG (2nd row), the experimental TA (3rd row), and the simulated TA (bottom row) data. 
@@ -683,13 +749,6 @@ complicated by sample-dependent system-bath coupling.  %
 	\label{psg:fig:nise_fits}
 \end{figure}
 
-The results of this final simulation are compared with the experimental data in Figure
-\ref{fig:nise_fits}.  %
-It is important to note that the simulations get many details of the rise-time spectra correct. 
-Particularly, the transition from narrow, diagonal line shape to a broad, less diagonal line shape is reproduced very well in both TA and TG simulations. 
-Such behavior is expected for responses from excitonic peaks of material systems; the rise time behavior for such systems was studied in detail previously.\cite{Kohler2017} 
-Because these simulations do not account for hot-exciton creation from the pump, simulations differ from experiment increasingly as the pump becomes bluer than the 1S center.
-
 \section{Conclusion}  % ===========================================================================
 
 By combining TA and TG measurements, we have described the complex third-order, 2D susceptibility
diff --git a/PbSe_susceptibility/chapter.tex b/PbSe_susceptibility/chapter.tex
index 8bacebf..c28b300 100644
--- a/PbSe_susceptibility/chapter.tex
+++ b/PbSe_susceptibility/chapter.tex
@@ -664,7 +664,8 @@ Our method gives agreement with the $\phi=0.25$ bleach factor.  %
 		$ |B| / |\gamma_\text{sol}| $ & $1.3 \cdot 10^5$ & $1.6 \cdot 10^5$ & 1.3 \\
 		$ |\gamma_{\nu1}| / |\gamma_\text{sol}| $ & $0.25 \pm 0.04 $ & $5.1$ & 20.4 \\
 	\end{tabular}
-  \caption{Non-linear parameters relative to CCl$_4$ hyperpolarizability.  $\gamma_{\nu 1}$:
+  \caption[Non-linear parameters relative to carbon tetrachloride hyperpolarizability.]{
+    Non-linear parameters relative to CCl$_4$ hyperpolarizability.  $\gamma_{\nu 1}$:
     hyperpolarizability of the $\nu_1$ Raman transition.}
 	\label{pss:tab:gamma_ratio}
 \end{table}
@@ -683,7 +684,8 @@ Our method gives agreement with the $\phi=0.25$ bleach factor.  %
 		& 4 & 7 & 120 \\
 		$\phi$ (\autoref{pss:eq:gamma_to_phi}) & 0.3 & 0.6 & 0.15 \\
 	\end{tabular}
-  \caption{Comparison of these measurements with PbX measurements in literature.
+  \caption[Comparison with PbX literature measurements.]{
+    Comparison of these measurements with PbX measurements in literature.
     $\gamma_\text{raman}^{(3)}$ refers to the 465 $\text{cm}^{-1}$ (symmetric stretch) mode of
     $\text{CCl}_4$}
 	\label{pss:tab:litcompare}
diff --git a/bibliography.bib b/bibliography.bib
index b6f050e..a9fbb04 100644
--- a/bibliography.bib
+++ b/bibliography.bib
@@ -1,14 +1,10 @@
-@article{HendersonGiles1994a,

-  author =       {Giles Henderson and Robert C. Rittenhouse and John C. Wright and Jon L. Holmes},

-  title =        {How a Photon is Created or Absorbed},

-  journal =      {Journal of Chemical Education},

-  volume =       71,

-  number =       4,

-  pages =        300,

-  year =         1994,

-  doi =          {10.1021/ed071p300.2},

-  month =        {apr},

-  publisher =    {American Chemical Society ({ACS})},

+@article{OliveroJJ1977a,

+  author =       {Olivero, J.J. and Longbothum, R.L.},

+  title =        {{Empirical fits to the voigt line width: A brief review}},

+  journal =      {J. Quant. Spectrosc. Radiat. Transfer},

+  volume =       17,

+  pages =        {233--236},

+  year =         1977,

 }

 

 @article{AartsmaThijsJ1976a,

@@ -552,8 +548,6 @@
   doi =          {10.1021/nn9001616},

   isbn =         {1936-0851},

   issn =         {1936-086X},

-  keywords =     {Composition,Molar extinction coefficient,PbSe,Semiconductor nanocrystal,Size

-                  dependence,informative,size-exciton correlation},

   month =        {jun},

   pmid =         19435305,

 }

@@ -588,22 +582,22 @@
 }

 

 @article{DeBoeijWimP1998a,

-author = {de Boeij, Wim P and Pshenichnikov, Maxim S and Wiersma, Douwe A},

-doi = {10.1146/annurev.physchem.49.1.99},

-isbn = {0066-426X},

-issn = {0066-426X},

-journal = {Annual Review of Physical Chemistry},

-keywords = {correlation function,ultrafast spectroscopy},

-month = {oct},

-number = {1},

-pages = {99--123},

-pmid = {15012426},

-title = {{ULTRAFAST SOLVATION DYNAMICS EXPLORED BY FEMTOSECOND PHOTON ECHO SPECTROSCOPIES}},

-volume = {49},

-year = {1998}

-}

-

-@article{DeGeyterBram2012a,

+  author =       {de Boeij, Wim P and Pshenichnikov, Maxim S and Wiersma, Douwe A},

+  title =        {{ULTRAFAST SOLVATION DYNAMICS EXPLORED BY FEMTOSECOND PHOTON ECHO

+                  SPECTROSCOPIES}},

+  journal =      {Annual Review of Physical Chemistry},

+  volume =       49,

+  number =       1,

+  pages =        {99--123},

+  year =         1998,

+  doi =          {10.1146/annurev.physchem.49.1.99},

+  isbn =         {0066-426X},

+  issn =         {0066-426X},

+  month =        {oct},

+  pmid =         15012426,

+}

+

+@article{gDeGeyterBram2012a,

   author =       {{De Geyter}, Bram and Geiregat, Pieter and Gao, Yunan and {Ten Cate}, Sybren and

                   Houtepen, Arjan J. and Schins, Juleon M. and {Van Thourhout}, Dries and

                   Siebbeles, Laurens D A and Hens, Zeger},

@@ -710,9 +704,6 @@ year = {1998}
   year =         2008,

   doi =          {10.1016/j.chemphys.2008.02.050},

   issn =         03010104,

-  keywords =     {Analytical spectroscopy,Methylene spectroscopy,Multidimensional infrared

-                  spectroscopy,Picosecond infrared spectroscopy,Polymer science,Protein

-                  spectroscopy,Time resolved infrared spectroscopy,Vibrational spectroscopy},

   month =        {jun},

 }

 

@@ -839,8 +830,6 @@ year = {1998}
   year =         1998,

   doi =          {10.1098/rsta.1998.0172},

   issn =         {1364-503X},

-  keywords =     {femtosecond spectroscopy,nonlinear spectroscopy,photon echoes,solvation

-                  dynamics,spectral broadening},

   month =        {feb},

 }

 

@@ -865,7 +854,6 @@ year = {1998}
   year =         1992,

   doi =          {10.1063/1.463565},

   issn =         {0021-9606},

-  keywords =     {molecular-dynamics high-resolution sodium vapor li},

   month =        {jul},

 }

 

@@ -960,6 +948,33 @@ year = {1998}
   month =        {dec},

 }

 

+@article{GdorItay2015a,

+  author =       {Gdor, Itay and Shapiro, Arthur and Yang, Chunfan and Yanover,

+                  Diana and Lifshitz, Efrat and Ruhman, Sanford},

+  title =        {{Three-pulse femtosecond spectroscopy of PbSe nanocrystals: 1S

+                  bleach nonlinearity and sub-band-edge excited-state absorption

+                  assignment}},

+  journal =      {ACS Nano},

+  volume =       9,

+  number =       2,

+  pages =        {2138--2147},

+  year =         2015,

+  doi =          {10.1021/nn5074868},

+}

+

+@article{GdorItay2015b,

+  author =       {Gdor, Itay and Yanover, Dianna and Yang, Chunfan and Shapiro, Arthur and

+                  Lifshitz, Efrat and Ruhman, Sanford},

+  title =        {{Three Pulse Femtosecond Spectroscopy of PbSe Nano-Crystals ; 1S Bleach

+                  Nonlinearity and Sub Band-Edge Excited State Absorption Assignment Three Pulse

+                  Femtosecond Spectroscopy of PbSe nano-Crystals ; 1S Bleach Nonlinearity and Sub

+                  Band-edge Excited State A}},

+  journal =      {ACS Nano},

+  pages =        {2--5},

+  year =         2015,

+  doi =          {10.1021/nn5074868},

+}

+

 @article{GeiregatPieter2014a,

   author =       {Geiregat, Pieter and Houtepen, Arjan J. and Justo, Yolanda and Grozema, Ferdinand

                   C. and {Van Thourhout}, Dries and Hens, Zeger},

@@ -1020,6 +1035,17 @@ year = {1998}
   pmid =         19449854,

 }

 

+@article{GesueleF2012a,

+  author =       {Gesuele, F and Sfeir, M Y and Murray, C B and Heinz, T F and Wong, C W},

+  title =        {{Biexcitonic Effects in Excited Ultrafast Supercontinuum S{\aa}ectroscopy of

+                  Carrier Multiplication and Biexctionic Effects in Excited States of PbS Quantum

+                  Dots}},

+  journal =      {Nano letters},

+  volume =       12,

+  pages =        2658,

+  year =         2012,

+}

+

 @misc{GitHub,

   note =         {Accessed: 2018-04-04},

   title =        {The world's leading software development platform · GitHub},

@@ -1177,6 +1203,19 @@ year = {1998}
   issn =         00796727,

 }

 

+@article{HendersonGiles1994a,

+  author =       {Giles Henderson and Robert C. Rittenhouse and John C. Wright and Jon L. Holmes},

+  title =        {How a Photon is Created or Absorbed},

+  journal =      {Journal of Chemical Education},

+  volume =       71,

+  number =       4,

+  pages =        300,

+  year =         1994,

+  doi =          {10.1021/ed071p300.2},

+  month =        {apr},

+  publisher =    {American Chemical Society ({ACS})},

+}

+

 @article{HoPP1979a,

   author =       {Ho, P P and Alfano, R R},

   title =        {{Optical Kerr effect in liquids}},

@@ -1190,6 +1229,19 @@ year = {1998}
   month =        {nov},

 }

 

+@article{HogemannClaudia1996a,

+  author =       {H{\"{o}}gemann, Claudia and Pauchard, Marc and Vauthey, Eric},

+  title =        {{Picosecond transient grating spectroscopy: The nature of the diffracted

+                  spectrum}},

+  journal =      {Review of Scientific Instruments},

+  volume =       67,

+  number =       10,

+  pages =        {3449--3453},

+  year =         1996,

+  doi =          {10.1063/1.1147157},

+  month =        {oct},

+}

+

 @article{HorakErikH2018a,

   author =       {Horak, Erik H. and Rea, Morgan T. and Heylman, Kevin D. and Gelbwaser-Klimovsky,

                   David and Saikin, Semion K. and Thompson, Blaise J. and Kohler, Daniel D. and

@@ -1217,6 +1269,17 @@ year = {1998}
   publisher =    {Institute of Electrical and Electronics Engineers ({IEEE})},

 }

 

+@article{HutchingsDC1992a,

+  author =       {Hutchings, D C and Sheik-Bahae, M and Hagan, D J and {Van Stryland}, E W},

+  title =        {{Kramers-Kronig relations in nonlinear optics}},

+  journal =      {Optical and Quantum Electronics},

+  volume =       24,

+  number =       1,

+  pages =        {1--30},

+  year =         1992,

+  doi =          {10.1007/BF01234275},

+}

+

 @article{HuygensChristiaan1913a,

   author =       {Huygens, Christiaan},

   title =        {Treatise on Light},

@@ -1251,7 +1314,6 @@ year = {1998}
   year =         2001,

   doi =          {10.1016/S0301-0104(01)00233-6},

   issn =         03010104,

-  keywords =     {inertial memory},

   mendeley-tags ={inertial memory},

   month =        {may},

 }

@@ -1347,6 +1409,32 @@ year = {1998}
   month =        {jan},

 }

 

+@article{KangInuk1997a,

+  author =       {Kang, Inuk and Wise, Frank W.},

+  title =        {{Electronic structure and optical properties of PbS and PbSe quantum dots}},

+  journal =      {Journal of the Optical Society of America B},

+  volume =       14,

+  number =       7,

+  pages =        1632,

+  year =         1997,

+  doi =          {10.1364/JOSAB.14.001632},

+  month =        {jul},

+}

+

+@article{KarkiKhadgaJ2013a,

+  author =       {Karki, Khadga J and Ma, Fei and Zheng, Kaibo and Zidek, Karel and Mousa,

+                  Abdelrazek and Abdellah, Mohamed and Messing, Maria E and Wallenberg, L Reine and

+                  Yartsev, Arkadi and Pullerits, T{\~{o}}nu},

+  title =        {{Multiple exciton generation in nano-crystals revisited: consistent calculation

+                  of the yield based on pump-probe spectroscopy.}},

+  journal =      {Scientific reports},

+  volume =       3,

+  pages =        2287,

+  year =         2013,

+  doi =          {10.1038/srep02287},

+  month =        {jan},

+}

+

 @article{KearnsNicholasM2017a,

   author =       {Nicholas M. Kearns and Randy D. Mehlenbacher and Andrew C. Jones and Martin T.

                   Zanni},

@@ -1435,6 +1523,15 @@ year = {1998}
   publisher =    {{AIP} Publishing},

 }

 

+@article{KraatzIngvarT2014a,

+  author =       {Kraatz, Ingvar T and Booth, Matthew and Whitaker, Benjamin J and Nix, Michael G D

+                  and Critchley, Kevin},

+  title =        {{Sub-Bandgap Emission and Intraband Defect-Related Excited-State Dynamics in

+                  Colloidal CuInS 2 /ZnS Quantum Dots Revealed by Femtosecond Pump − Dump − Probe

+                  Spectroscopy}},

+  year =         {2014},

+}

+

 @article{KumarNardeep2013a,

   author =       {Kumar, Nardeep and He, Jiaqi and He, Dawei and Wang, Yongsheng

                   and Zhao, Hui},

@@ -1637,7 +1734,6 @@ year = {1998}
   publisher =    {American Chemical Society ({ACS})},

 }

 

-

 @article{MaimanTheodore1960a,

 	author = {Maiman, Theodore},

 	title = {Stimulated Optical Radiation in Ruby},

@@ -1735,6 +1831,7 @@ year = {1998}
   month =        {may},

 }

 

+

 @article{McClainBrianL2004a,

   author =       {McClain, Brian L. and Finkelstein, Ilya J. and Fayer, M. D.},

   title =        {Vibrational echo experiments on red blood cells: Comparison of the dynamics of

@@ -1828,7 +1925,6 @@ year = {1998}
   year =         2007,

   doi =          {10.1021/cm071410q},

   issn =         {0897-4756},

-  keywords =     {reference,size-exciton correlation},

   mendeley-tags ={reference},

   month =        {dec},

 }

@@ -1854,6 +1950,7 @@ year = {1998}
   isbn =         9780195092783,

   series =       {Oxford series in optical and imaging sciences},

 }

+

 @article{MukamelShaul2000a,

   author =       {Mukamel, Shaul},

   title =        {Multidimensional Femtosecond Correlation Spectroscopies of Electronic and

@@ -1961,10 +2058,8 @@ year = {1998}
   doi =          {10.1016/S0009-2614(99)01173-2},

   isbn =         {0009-2614},

   issn =         00092614,

-  keywords =     {3rd-order,dynamics,glasses,phase,probes,spectroscopy},

   month =        {dec},

 }

-

 @book{OliphantTravisE2006a,

   author =       {Travis E Oliphant},

   title =        {A guide to NumPy},

@@ -2190,6 +2285,13 @@ year = {1998}
   month =        {feb},

 }

 

+@article{PetersonJJ2007a,

+  author =       {Peterson, JJ and Huang, Libai and Delerue, C. and Allan, Guy},

+  title =        {{Uncovering forbidden optical transitions in PbSe nanocrystals}},

+  journal =      {nano Letters},

+  year =         2007,

+}

+

 @article{PollardW1992a,

   author =       {Pollard, W.},

   title =        {{Analysis of Femtosecond Dynamic Absorption Spectra of Nonstationary States}},

@@ -2226,6 +2328,21 @@ year = {1998}
   publisher =    {{ACM} Press},

 }

 

+@article{PriceMichaelB2015a,

+  author =       {Price, Michael B. and Butkus, Justinas and Jellicoe, Tom C. and Sadhanala, Aditya

+                  and Briane, Anouk and Halpert, Jonathan E. and Broch, Katharina and Hodgkiss,

+                  Justin M. and Friend, Richard H. and Deschler, Felix},

+  title =        {{Hot-carrier cooling and photoinduced refractive index changes in

+                  organic–inorganic lead halide perovskites}},

+  journal =      {Nature Communications},

+  volume =       6,

+  number =       {May},

+  pages =        8420,

+  year =         2015,

+  doi =          {10.1038/ncomms9420},

+  publisher =    {Nature Publishing Group},

+}

+

 @misc{PyPI,

   note =         {Accessed: 2018-04-04},

   title =        {PyPI - the Python Package Index},

@@ -2269,7 +2386,6 @@ year = {1998}
   year =         2008,

   doi =          {10.1364/JOSAB.25.001738},

   issn =         {0740-3224},

-  keywords =     {1904400,1904710,OCIS codes},

   month =        {oct},

 }

 

@@ -2358,7 +2474,6 @@ year = {1998}
   year =         2008,

   doi =          {10.1016/j.chemphys.2007.12.020},

   issn =         03010104,

-  keywords =     {cdse,laser,nanocrystal,non-linear response,photo echo,quantum dot,ultrafast},

   month =        {jun},

 }

 

@@ -2371,8 +2486,6 @@ year = {1998}
   year =         2008,

   doi =          {10.1109/MS.2008.84},

   ISSN =         {0740-7459},

-  keywords =     {Professional end-user developer; scientific application software; software

-                  development; software testing; empirical study},

   month =        07,

 }

 

@@ -2391,6 +2504,20 @@ year = {1998}
   month =        {dec},

 }

 

+@article{SchinsJuleon2009a,

+  author =       {Schins, Juleon and Trinh, M. and Houtepen, Arjan and Siebbeles, Laurens},

+  title =        {{Probing formally forbidden optical transitions in PbSe nanocrystals by time- and

+                  energy-resolved transient absorption spectroscopy}},

+  journal =      {Physical Review B},

+  volume =       80,

+  number =       3,

+  pages =        035323,

+  year =         2009,

+  doi =          {10.1103/PhysRevB.80.035323},

+  issn =         {1098-0121},

+  month =        {jul},

+}

+

 @article{SchweigertIgorV2008a,

   author =       {Schweigert, Igor V. and Mukamel, Shaul},

   title =        {{Simulating multidimensional optical wave-mixing signals with finite-pulse

@@ -2617,6 +2744,21 @@ year = {1998}
   isbn =         {0824742435},

 }

 

+@article{SvirkoP1999a,

+  author =       {{P. Svirko}, Yuri and Shirane, Masayuki and Suzuura, Hidekatsu and

+                  Kuwata-Gonokami, Makoto},

+  title =        {{Four-Wave Mixing Theory at the Excitonic Resonance: Weakly Interacting Boson

+                  Model}},

+  journal =      {Journal of the Physical Society of Japan},

+  volume =       68,

+  number =       2,

+  pages =        {674--682},

+  year =         1999,

+  doi =          {10.1143/JPSJ.68.674},

+  issn =         {0031-9015},

+  month =        {feb},

+}

+

 @article{Thalhammer1983,

   author =       {Thalhammer, M and Penzkofer, A},

   title =        {{Measurement of third-order nonlinear susceptibilities by non-phase matched

@@ -2827,6 +2969,19 @@ year = {1998}
   month =        {oct},

 }

 

+@article{WehrenbergBrianL2002a,

+  author =       {Wehrenberg, Brian L. and Wang, Congjun and Guyot-Sionnest, Philippe},

+  title =        {{Interband and Intraband Optical Studies of PbSe Colloidal Quantum Dots}},

+  journal =      {The Journal of Physical Chemistry B},

+  volume =       106,

+  number =       41,

+  pages =        {10634--10640},

+  year =         2002,

+  doi =          {10.1021/jp021187e},

+  issn =         {1520-6106},

+  month =        {oct},

+}

+

 @article{Weiner1991,

   author =       {Weiner, A M and Wiederrecht, Gary P and Nelson, Keith A and Leaird, D E},

   title =        {{Femtosecond multiple-pulse impulsive stimulated Raman scattering spectroscopy}},

@@ -2880,8 +3035,6 @@ year = {1998}
   author =       {G. Wilson},

   doi =          {10.1109/MCSE.2006.122},

   journal =      {Computing in Science {\&} Engineering},

-  keywords =     {physics education; computation in undergraduate education; continuing education;

-                  software engineering},

   month =        11,

   pages =        {66-69},

   title =        {Software Carpentry: Getting Scientists to Write Better Code by Making Them More

@@ -2905,6 +3058,18 @@ year = {1998}
   publisher =    {Public Library of Science ({PLoS})},

 }

 

+@article{WiseFrankW2000a,

+  author =       {Wise, Frank W.},

+  title =        {{Lead salt quantum dots: the limit of strong quantum confinement.}},

+  journal =      {Accounts of chemical research},

+  volume =       33,

+  number =       11,

+  pages =        {773--80},

+  year =         2000,

+  month =        {nov},

+  pmid =         11087314,

+}

+

 @article{WongCathyY2011a,

   author =       {Cathy Y. Wong and Gregory D. Scholes},

   title =        {Using two-dimensional photon echo spectroscopy to probe the fine structure of the

@@ -2977,7 +3142,6 @@ year = {1998}
   month =        {oct},

 }

 

-

 @article{WrightJohnCurtis2011a,

   author =       {John C. Wright},

   title =        {Multiresonant Coherent Multidimensional Spectroscopy},

@@ -3005,6 +3169,7 @@ year = {1998}
   publisher =    {Elsevier B.V.},

 }

 

+

 @article{XiaoDi2012a,

   author =       {Di Xiao and Gui-Bin Liu and Wanxiang Feng and Xiaodong Xu and

                   Wang Yao},

@@ -3044,6 +3209,31 @@ year = {1998}
 	doi = {10.1038/nphys2942},

 }

 

+@article{YangYe2015a,

+  author =       {Yang, Ye and Yan, Yong and Yang, Mengjin and Choi, Sukgeun and Zhu, Kai and

+                  Luther, Joseph M and Beard, Matthew C},

+  title =        {{Low surface recombination velocity in}},

+  journal =      {Nature Communications},

+  volume =       6,

+  pages =        {1--6},

+  year =         2015,

+  doi =          {10.1038/ncomms8961},

+  publisher =    {Nature Publishing Group},

+}

+

+@article{YangYe2017a,

+  author =       {Yang, Ye and Yang, Mengjin and Moore, David?T. and Yan, Yong and Miller, Elisa M.

+                  and Zhu, Kai and Beard, Matthew?C.},

+  title =        {{Top and bottom surfaces limit carrier lifetime in lead iodide perovskite films}},

+  journal =      {Nature Energy},

+  volume =       2,

+  number =       2,

+  pages =        16207,

+  year =         2017,

+  doi =          {10.1038/nenergy.2016.207},

+  month =        {jan},

+}

+

 @article{YeeTK1978a,

   author =       {Yee, T. K. and Gustafson, T. K.},

   title =        {Diagrammatic analysis of the density operator for nonlinear optical calculations:

diff --git a/dissertation.tex b/dissertation.tex
index 72b351e..9c0d90a 100644
--- a/dissertation.tex
+++ b/dissertation.tex
@@ -89,11 +89,10 @@ This dissertation is approved by the following members of the Final Oral Committ
 \part{Appendix} \label{prt:appendix}

 \begin{appendix}

 % AFTER DEFENSE \include{public/chapter}

-\include{procedures/chapter}

-\include{hardware/chapter}

-\include{irf/chapter}

-\include{quantitative_ta/chapter}

-\include{tg/chapter}

+\include{procedures/chapter}  % prc

+\include{hardware/chapter}  % hdw

+\include{irf/chapter}  % irf

+\include{quantitative_ta/chapter}  % qta

 \include{opa_phase/chapter}  % rpa

 %\include{errata/chapter}

 % AFTER DEFENSE \include{colophon/chapter}

diff --git a/errata/chapter.tex b/errata/chapter.tex
deleted file mode 100644
index 159c713..0000000
--- a/errata/chapter.tex
+++ /dev/null
@@ -1,25 +0,0 @@
-\chapter{Errata}
-
-In this appendix I list the errors I am currently aware of in my publications and associated
-work.  %
-
-\section{Czech 2015}  % ---------------------------------------------------------------------------
-
-The following is an errata for \textcite{CzechKyleJonathan2015a}, published in November 2015.
-
-\begin{itemize}
-  \item Reference 13 is identical to reference 9.
-  \item In the last paragraph of the introduction the sentence ``The experimental spectra differ
-    from the simple 2D spectrum shown in Figure 1d and those of earlier CMDS experiments with model
-    systems'' appears. This sentence cites references 6 through 10. Instead, it should cite
-    references 15 through 20.
-  \item In the last paragraph begining on page 12148, the text ``Automated delay stages and neutral
-    density filters set the excitation time delays over all values of $\tau_{21}$ with
-    $\tau_{22}=0$''    appears. For the second $\tau$, the subscript should read $22^\prime$, not
-    $22$.
-  \item Caption of Figure 5 reads, in part: "showing the impact of the $\omega_1$ excitation
-    frequency on the $\omega_1$ spectral line shape". This should instead read "showing the impact
-    of the $\omega_1$ excitation frequency on the $\omega_2$ spectral line shape". The subscript on
-    the last $\omega$ should be a 2 and not a 1.
-  \item Figure 6 e$^\prime$+e,e$^\prime$ should read e$^\prime$+e,e and vice versa.
-\end{itemize}
\ No newline at end of file
diff --git a/introduction/chapter.tex b/introduction/chapter.tex
index fe43cb4..f6adc46 100644
--- a/introduction/chapter.tex
+++ b/introduction/chapter.tex
@@ -205,7 +205,7 @@ In \autoref{cha:mx2} we report the first MR-CMDS study performed on a molybdenum
 film.  %

 MoS\textsubscript{2} is a member of a class of materials called transition metal dichacolgenides

 which have recently attracted a large amount of attention for their unique photophysical

-properties.  %

+properties. \cite{MakKinFai2010a, XuXiaodong2014a, XiaoDi2012a}  %

 These thin films have relatively low optical density and are highly scattering, making them

 particularly challenging for MR-CMDS experiments.  %

 We employed several strategies to overcome these challenges, and performed three dimensional

diff --git a/todo.org b/todo.org
index 45021e4..de66d3f 100644
--- a/todo.org
+++ b/todo.org
@@ -56,7 +56,7 @@
 * TODO generalizability                                                 :opa:
 * TODO future directions                                                :opa:
 * TODO figure captions                                                  :opa:
-* DONE insert content from main                                         :pss:
+* DONE insert content from main              0:00:00                            :pss:
   CLOSED: [2018-04-12 Thu 21:04]
 * IDEA insert content from SI                                           :pss:
 * DONE insert content from main                                         :psg: