\section{INTRODUCTION}
\label{sec:intro}

\textsc{Polarbear}\ is designed to measure the B-mode \ac{cmb} polarization signature due to gravitational lensing and search for the possible B-mode signature of inflation. Kermish \textit{et al.}\cite{Kermish_SPIE2012}, also at this conference [8452-47], presents the scientific goals and instrumental overview of \textsc{Polarbear}. In this proceeding, we present the details of the antenna-coupled bolometric focal plane detector array used in the \textsc{Polarbear}\ experiment.

To achieve the high sensitivity and stringent control of systematic instrumental effects required by \textsc{Polarbear}'s scientific goals, there are several important detector design criteria: The individual detectors must couple optical power from the free-space telescope optics to the bolometer with minimal loss. The spectral band of the detectors, centered at 148 GHz, must be designed to utilize the atmospheric window around that frequency. The detectors must be dual-polarization, with diffraction-limited spatial response that are symmetric between polarizations. The bolometer time constants must be such that they can measure the \ac{cmb} structure in the sky as the telescope scans across it. And the bolometer sensitivity must be background-limited, recognising that the statistical variations in the photon flux from the atmosphere are the fundametal limit to detector sensitivity for ground-based observation of the \ac{cmb}.

Section \ref{sec:fparc} of this proceeding goes over the design of the focal plane array of detectors and how it addresses these goals, while Section \ref{sec:fab} gives the details of the detector fabrication in the Berkeley Nanolab. Section \ref{sec:valid} presents the results of tests conducted to validate the focal plane performance. Section \ref{sec:future} concludes with some comments about the future of this technology and of the \textsc{Polarbear}\ experiment.

\section{Focal plane array architecture}
\label{sec:fparc}

   \begin{figure}
   \begin{center}
	\begin{tabular}{cc}
   \includegraphics[width=2.5in]{FocalPlaneArray.eps} & \includegraphics[width=2.4in]{segreCropped.eps} \\
\textbf{(a)} & \textbf{(b)}
\end{tabular}
   \end{center}
   \caption[Array]{\label{fig:arraypicture} \textbf{(a)} A picture of the entire assembled \textsc{Polarbear}\ focal plane. The 6 silver sub-arrays use silicon lenslets with thermoformed \acf{pei} \acf{ar} coatings, while the single white array uses alumina lenslets, a material with similar optical properties to silicon. \textbf{(b)} A picture of a single sub-array. Note how the readout folds back behind the sub-array for a compact, scalable design.}
   \end{figure}

The \textsc{Polarbear}\ array consists of 637 pixels (1274 antenna-coupled bolometers) on 7 hexagonal sub-arrays. Each of these sub-arrays is fabricated on a single 4'' silicon wafer in the Berkeley Nanolab. Pictures of the assembled array and a single assembled sub-array are shown in Figure~ \ref{fig:arraypicture}.

The \textsc{Polarbear}\ detectors integrate several features into a single pixel: polarization-sensitive antennas separate the signal into polarized components and couple these components from free space into superconducting microstrip waveguide, spectral bandpass filtering transmits only the desired frequency band, and superconducting \ac{tes} thermistors on the thermally released bolometers provide the photon-noise limited detection of the transmitted electromagnetic wave. An overview of this architecture is shown in Figures \ref{fig:ellipticalLenslet} and \ref{fig:pixel}.

   \begin{figure}
   \begin{center}
\begin{tabular}{ccc}
\includegraphics[width=2in]{ElliplicalLensletManyRays.eps} & \includegraphics[width=2.0in]{ExtendedHemisphereLabels.eps} & \includegraphics[width=2.0in]{lensExitAngle.eps} \\
\textbf{(a)} & \textbf{(b)} & \textbf{(c)}
\end{tabular}
   \end{center}
   \caption[Refraction at the surface of directly contacting dielectric lenslets]{\label{fig:ellipticalLenslet}Refraction at the surface of a dielectric lenslet for {\bf(a)} an elliptical lenslet and {\bf(b)} an extended hemisphere approximating the elliptical lenslet shape. In the case of the elliptical lenslet, all rays emanating from the center exit the lenslet vertically. In the case of the extended hemisphere, the transmitted angle depends on the extension length and the position on the lenslet surface. \textbf{(c)} Geometric calculation of the transmitted ray angle when exiting the lenslet surface for rays emanating from the center of the planar crossed double-slot dipole. For a true ellipse, this angle would be zero at all points along the radius of the lenslet. The five different lines are for different extension lengths (L) related to the pixel radius as shown in the legend. The vertical lines mark the point after which total internal reflection occurs.}
   \end{figure}

\subsection{Antenna design: lenslet-coupled crossed double-slot dipole}
\label{sec:antennaDesign}

The idea of increasing an antenna's directivity with a directly contacting dielectric lenslet has been used across many applications and spectral regions \cite{RutledgeFirstDielLens,Zmuidzinas_quasi-optical,Filipovic_double_slot,LensletReferenceMap}. This configuration has several benefits:

\begin{itemize}

\item Because there is a difference in near-field dielectric constant between the sides of the planar antenna, the beam is formed more strongly in the direction of the high dielectric constant. In the case of silicon, with a relative dielectric constant of 11.7, the antenna gain is 10 times higher in the forward direction (the direction of the lenslet) than in the backward direction\cite{Filipovic_double_slot}.

\item Standard lithographic techniques produce planar antennas on dielectric substrates. While the gain is preferentially in the direction of the substrate as described above, total internal reflection at the substrate boundry leads the antenna to be sensitive to modes propagating within the planar substrate. An extended hemisphere of dielectric material limits that total internal reflection.

\item The hemispherical dielectric lenslet acts as a beam-forming element, coupling the planar structure's broad beam within the dielectric to a more directed beam in free space---the lenslet acts to magnify the active size of the antenna. This leaves room under the lenslet (away from its center) for other microwave elements, without using the valuable focal plane area coupled to the diffraction-limited field of view of the telescope.

\item Refracting optics have a wide bandwidth of operation. Given that the lenslet medium has low loss, the loss in coupling to free space is dominated by the reflection loss at the surface of the lenslet. This can be mitigated using an \ac{ar} coating. For the \textsc{Polarbear}\ pixel, a single-layer quarter-wavelength \ac{ar} coating is sufficient to reduce the reflection loss over the spectral band of interest. In the future, multichroic antennas coupled to lenslets with multi-layer antireflection coatings will have the capability to measure several spectral bands through a single pixel \cite{Suzuki_SPIE2012,OBrient:Thesis,OBrient2008SPIE}.

\end{itemize}

In the geometric limit, rays from a point source emitter within a dielectric will refract to be parallel if the surface of the dielectric is the far side of an ellipse of the correct eccentricity \cite{Hecht}. Figure~ \ref{fig:ellipticalLenslet}{\bf(a)} shows this refraction for an ellipse of the eccentricity required given the optical index of silicon. For lenslets of high dielectric constant this eccentricity is small, and thus the ellipse can be well-approximated by an extended hemisphere. Figure \ref{fig:ellipticalLenslet}{\bf(b)} shows a cross section of such an extended hemispherical lenslet, and Figure~ \ref{fig:ellipticalLenslet}\textbf{(c)} shows the transmitted angle for such a lenslet as a function of position on the lenslet surface, for several extension lengths. These angles are calculated using geometric optics; to the extent the angles are close to zero, the electromagnetic wave couples to a diffraction-limited beam with waist diameter slightly smaller than the lenslet diameter, just as with a horn antenna.

The nominal \textsc{Polarbear}\ pixel size is $R=1.6\lambda_0$, with $\lambda_0$ the free-space wavelength of the center of the spectral band. For a range of lenslet extension lengths $L$ between $0.36R$ and $0.39R$, the lenslet-coupled double-slot dipole antenna produces a beam with high Gaussisity and directivity \cite{Filipovic_double_slot}. The \textsc{Polarbear}\ dielectric lenslet extensions have some fabrication variability; the above values were chosen as bounds for acceptable extension thickness.

For lenslets of high dielectric constant, the reflection loss at the surface of the lenslet is significant. This reflection loss can be mitigated by employing a quarter-wavelength \ac{ar} coating at the lenslet surface. As the extension length increases, the reflectance at the quarter-wavelength coating increases because of the non-normal incidence of the radiation and the curvature of the surface \cite{Filipovic_double_slot}. A re-optimization of the coating thickness as a function of distance from the center of the lenslet dictates that the coating be slightly thicker at the side of the lenslet, but this optimization was not attempted for the \textsc{Polarbear}\ lenslets, which were \ac{ar} coated with uniform-thickness thermoformed polyetherimide, a plastic of optical index 1.7.

   \begin{figure}
   \begin{center}
	\begin{tabular}{ccc}
\includegraphics[width=1.5in]{antennacad.eps} & \includegraphics[width=2.5in]{pixel_wide.eps} & \includegraphics[width=2.2in]{bolometer_withnumbers.eps} \\
\textbf{(a)} & \textbf{(b)} & \textbf{(c)}
\end{tabular}
   \end{center}
   \caption[Pixel photograph]{\label{fig:pixel}\textbf{(a)} Design drawing of the planar crossed double-slot dipole. Areas in green represent slots in the niobium ground plane. \textbf{(i)} Microstrip feeds are shown at the center of each rectangular slot, with a \textbf{(ii)} via from the microstrip to the ground plane (in purple) internal to each slot. The incoming radiation excites a voltage across the slot and current around the slot, which excites a propagating wave in the microstrip. For the \textsc{Polarbear}\ pixel, $W=33.9~\mu$m ($0.017\lambda_0$), $S=338.7~\mu$m ($0.167\lambda_0$), and $L=677.5~\mu$m ($0.334\lambda_0$). \textbf{(b)} A photograph of the planar structures in one of the \textsc{Polarbear}\ pixels. The structures are labelled in the order that the \ac{cmb} photon travels through them: {\bf(i)} Crossed double-slot dipole, {\bf(ii)} Dolph-Chebychev microstrip transformer {\bf(iii)} microstrip cross-under, {\bf(iv)} cross-under balancing structures, {\bf(v)} microstrip filters, and {\bf(vi)} bolometers. \textbf{(c)} Scanning electron micrograph of a \textsc{Polarbear}\ bolometer, showing \textbf{(i)} the resistive load onto which the observed optical power is dissipated, \textbf{(ii)} the dual-transition \ac{tes} that senses modulation of this power, and \textbf{(iii)} the gold added to the bolometer to increase its heat capacity.}
   \end{figure}

The planar structure coupled to the dielectric lenslet must be able to couple photons efficiently into the on-wafer  waveguide across the desired spectral band, have low cross-polarized response across its spatial beam (at least within the angles that transmit through the aperture stop), and couple efficiently through the lenslet to the diffraction-limited telescope optics. The resonant crossed double-slot dipole used in \textsc{Polarbear}\ was adapted from the antenna presented in Chattopadhyay and Zmuidzinas \cite{doubleSlotAntennaLetter}. This antenna is useful for its relatively low drive impedance, low cross-polarized response, and circular beams \cite{doubleSlotAntennaTest,OBrient:Thesis}.

\subsection{Dolph-Chebychev microstrip transformer}

The antenna described above has a resonant impedance of $30~\Omega$ \cite{doubleSlotAntennaLetter,OBrient:Thesis}. This impedance is achievable with the microstrip waveguide used to feed the antenna, but because of limitations in the resolution and repeatability of the optical lithography used to print the microstrip, wider (lower impedance) microstrip can be fabricated with greater repeatability. For this reason, the band pass filter is designed to couple to a lower input and output impedance of $10~\Omega$. A microstrip transformer is used to transition from the antenna impedance to the filter input impedance. Four of these microstrip transformers are used in each pixel, as shown in Figure~ \ref{fig:pixel}. As a function of the length along the transformer, $l$, the design microstrip impedance, $Z$, is given by \cite{transformerApproxAndTest}

\begin{equation}
Z(l) = Z_{antenna} \cdot e^{ \frac{1}{2} \ln \left ( \frac{Z_{filter}}{Z_{antenna}}\right ) \left [\sin \left (\pi\left ( \frac{l}{L}-\frac{1}{2} \right ) \right )+1 \right ] },
\end{equation}

\noindent where $L$ is the total length of the transformer. The transmittance of that transformer, including the effect of the dielectric loss in the microstrip, is shown in Figure~ \ref{fig:xfmr}.

\begin{figure}[htpb]
\begin{center}
\includegraphics[width=2.5in]{impedancexfmr_nocross.eps}
\end{center}
\caption[Transmittance of the microstrip transformer]{\label{fig:xfmr}Transmitted power ($\left | S_{21} \right | ^2$) as a function of frequency for the microstrip transformer that connects the antenna to the microstrip filter. Also plotted is the transmittance due only to dielectric loss in a line of equal length.}
\end{figure}

\subsection{Microstrip cross-under}
\label{sec:cross-under}

\begin{figure}[htpb]
\begin{center}
\begin{tabular}{cc}
\includegraphics[width=2.0in]{crossundercs.eps} & \includegraphics[width=1.8in]{crossoverWithCaptions.eps} \\
{\bf(a)} & {\bf(b)}
\end{tabular}
\end{center}
\caption[Microstrip cross-under design and simulated performance]{\label{fig:crosssim}Microstrip cross-under {\bf(a)} design and {\bf(b)} simulated performance. {\bf(a)} shows cross-section views from two angles. In the cross-sections, the vertical scale is enhanced by a factor of 10 for clarity. Layers shown are \textbf{(i)} Low-stress silicon nitride, \textbf{(ii)} Nb ground plane, \textbf{(iii)} SiO$_2$ crossunder isolation, \textbf{(iv)} Nb cross-under, \textbf{(v)} SiO$_2$ dielectric, and \textbf{(vi)} Nb microstrip conductor. {\bf(b)} Transmittance, reflection, and crosstalk as simulated in Sonnet.}
\end{figure}

Because of the topology of the antenna, the signals in two of the microstrip lines must cross each other. This is accomplished using a microstrip cross-under structure. The design of this structure is shown in Figure~ \ref{fig:crosssim}{\bf(a)}. The cross-under was simulated using the Sonnet full-wave simulation package\footnote{http://www.sonnetsoftware.com/}. The results of this simulation are shown in figure \ref{fig:crosssim}{\bf(b)}. The reflected power is less than 1\% and the crosstalk power is less than 0.01\%.

\subsection{Bandpass filter}

The bandpass filter is a standard resonant filter design: a shorted quarter-wavelength stub filter \cite{MstripFilterDesignEquations}. Note that the stubs are narrowed where they meet the trunkline. This is necessary for the phase fronts of the wave on the low impedance stub to be perpendicular to the stub length. When the superconducting niobium microstrip is correctly modeled \cite{KerrSonnetCorrection}, an analytical network model of the filter closely reflects the results of a full-wave simulation using Sonnet. Using the analytical model, the parameters of the microstrip filter---center frequency, bandwidth, and steepness of edge (number of poles)---were optimized given the expected atmospheric emission and absorption profile.

\subsection{Bolometer structure}

A bolometer is fundamentally a place where radiative power is deposited that is in thermal contact with a thermistor of heat capacity $C$, all of which has a weak connection to a thermal bath.

In the case of \textsc{Polarbear}, the weak thermal link between the thermistor and the thermal bath is provided by a silicon nitride suspension, where the silicon substrate is etched away beneath the suspension. The scanning electron micrograph in Figure~ \ref{fig:pixel}\textbf{(c)} shows this suspension. The thermal path between the thermistor and the silicon wafer is along the silicon nitride suspension, which is $\sim 1~\mu$m in thickness. The operating power of the bolometer is set by the thermal conductance of this suspension and the thermistor's superconducting transition temperature.

The electrothermal bandwidth of the bolometer $\Delta f_{et}$ is proportional to its small-signal thermal conductance, $g$, the electrothermal loop gain $\mathcal{L}$, and inversely proportional to its heat capacity, $C$ \cite{Gildemeister:Thesis}. Whatever the electrothermal bandwidth, the readout bandwidth for that bolometer must be larger by a factor of $\sim 5$ to stably voltage bias the bolometer\cite{Irwin:etbw}.  The scanning speed of the telescope limits the bandwidth over which there is useful astrophysical signal. The limited readout bandwidth available can best be used multiplexing several bolometers, each with only the necessary electrothermal bandwidth, rather than fewer bolometers of larger $\Delta f_{et}$. Since $g$ is set by the bolometer properties and the desired operating power, $C$ is controlled in fabrication to achieve the desired $\Delta f_{et}$ at reasonable $\mathcal{L}$.

The head capacity, $C$, is controlled by depositing gold on the bolometer and making sure the gold is in good thermal contact with the thermistor \cite{JaredThesis}. Figure~ \ref{fig:pixel}\textbf{(c)} shows the gold used for this purpose. The thermistor material extends underneath the gold, and the gold extends onto the thermistor leads, so as to make the best thermal contact possible without strongly affecting current flow across the thermistor.

\subsubsection{Resistive load}

The optimal resistive load for absorbing power from two balanced $10~\Omega$ microstrip lines has a real impedance of $20~\Omega$ and no imaginary impedance, creating a virtual ground in the center of the load resistor. The \textsc{Polarbear}\ load resistor is fabricated from titanium with a sheet resistance $\sim 10~\Omega/sq$. It's length was adjusted to achieve a measured DC resistance of $20~\Omega$.

\subsubsection{Dual-TES thermistors}
\label{sec:dualtes}

The \textsc{Polarbear}\ bolometers employ a dual-\ac{tes} thermistor \cite{DualTES} to allow for two modes of observation: a low noise mode for cosmological observations and a higher noise mode with larger dynamic range. The low noise mode uses an aluminum/titanium bilayer thermistor with a transition temperature of $T_{c,bilayer} \approx 0.5$ K, while the mode with larger dynamic range employs a bare aluminum thermistor with a transition temperature of $T_{c,Al} = 1.4$ K. The design of the bilayer thermistor is somewhat independent of the Al thermistor. When operating the bolometer with the bilayer thermistor, the Al thermistor is simply a part of the superconducting thermistor leads. The bilayer thermistor operating resistance is $\sim 1~\Omega$, which provides the bolometer responsivity necessary and couples correctly to the readout system.

\section{Focal plane array fabrication}
\label{sec:fab}

A single sub-array consists of a hexagonal device wafer of antenna-coupled bolometers on one side, with its other side in contact with a hexagonal array of lenslets. This section describes the fabrication of both of these items.

\subsection{Lenslet sub-array}

The hemispherical lenslets are produced from semiconductor-grade single-crystal silicon. They are ground into spheres and then ground flat along a plane intersecting a great circle of the sphere\footnote{The commercial grinding process is done by Rayotek Scientific, Inc.}. This leaves a hemisphere to be integrated with a flat silicon spacer and the antenna as described in Section~ \ref{sec:antennaDesign}.

The lenslet sub-array acts both as part of the extension length providing the synthesized elliptical lenslet shape, and as an alignment wafer to properly center each lenslet on its associated crossed double-slot dipole. Circles of diameter slightly ($\sim 40~\mu$m) larger than the diameter of the lenslets are etched into the lenslet wafer using a deep reactive ion etching process \cite{Bosch}. These precise circular depressions serve to locate each lenslet to within $\sim \lambda_0/100$. The entire lenslet wafer is then optically aligned to the device wafer using transmission infrared microscopy.

\subsection{Antenna-coupled bolometer sub-array}

   \begin{table}[htpb]
   \begin{center}

	\begin{tabular}{|c|c|c|c|}
\hline
\textbf{Material} & \textbf{\shortstack{Thickness\\(microns)}} &\textbf{Deposition Process} & \textbf{Etch Process} \\
\hline
Silicon substrate & 500 & Prime-grade wafers & XeF$_2$ vapor \\
Silicon dioxide & 0.05 & Thermal growth & CF$_4$ RIE \\
Low stress nitride & 0.7 & LPCVD & CF$_4$ RIE \\
Niobium & 0.3 & DC magnetron sputter & CF$_4$ / O$_2$ RIE \\
Silicon dioxide & 0.3 & 350\textdegree C PECVD & CHF$_3$ / O$_2$ RIE \\
Niobium & 0.3 & DC magnetron sputter & CF$_4$ / O$_2$ RIE \\
Silicon dioxide & 0.5 & 350\textdegree C PECVD & CHF$_3$ / O$_2$ RIE \\
Niobium & 0.6 & DC magnetron sputter & CF$_4$ / O$_2$ RIE \\
Aluminum/titanium bilayer & 0.04/0.12 & DC magnetron sputter & Premixed wet etch / SF$_6$ RIE \\
Gold & 1.5 & Evaporation & Liftoff in acetone \\
\hline

	\end{tabular}
	\end{center}
	\caption[Layers of the \textsc{Polarbear}\ device wafer]{\label{tab:process}Layers in the \textsc{Polarbear}\ device wafer, along with thickness, deposition process, and etching process. Acronyms used in the table: RIE: Reactive Ion Etch; CVD: Chemical Vapor Deposition; PECVD: Plasma-enhanced CVD; LPCVD: Low-Pressure CVD. }
	\end{table}

The \textsc{Polarbear}\ device wafers require an eight layer process on a silicon substrate. The process is outlined in Table~ \ref{tab:process}, in order of layer deposition. The etching process for each layer is given; the etches for the layer are not in order, because they do not all occur immediately after deposition.

\subsubsection{Substrate}

The silicon substrate forms part of the dielectric contacting lenslet, so its microwave loss properties are important. A \ac{fts} was used to measure loss in commercial semiconductor-grade silicon at the frequencies of interest. The upper bound for loss on this measurement was low enough that it was not a significant source of loss in the pixel. It is important to note that some silicon with a higher impurity content did have significant loss and was unacceptable for this application.

\subsubsection{Photoresist coating and patterning}

Patterning is done in standard positive photoresist with a thickness of 1--4 microns, where the exact recipe depends on the process step. Patterning of the individual pixels is accomplished using a step and repeat 10X optical reduction printer, a GCA wafer stepper. Test pixels, wafer border, and alignment marks are also printed by the wafer stepper. The niobium thermistor leads, which transmit signals from the edge of the wafer to the thermistors, are patterned using a wafer-scale contact printer. After each patterning step, the photoresist undergoes a post-exposure bake and then an automated puddle develop. After development, the wafer is exposed to a 1 minute low-power oxygen plasma to remove any residual photoresist where the pattern was developed.

\subsubsection{Microwave structures}

The continuous niobium ground plane is broken only where necessary for the bolometer release and for the slots of the crossed double-slot dipole. Those slots are connected to superconducting microstrip waveguides that consist of a niobium ground plane, silicon dioxide dielectric, and niobium top conductor. The dielectric loss in the silicon dioxide is the primary source of loss in this superconducting microstrip transmission line \cite{Myers:Thesis}.

The microwave structures, including the cross-under, are constructed of three layers of niobium and two layers of silicon dioxide. Note that the silicon dioxide is not deposited in a furnace. This is important because high temperatures degrade the quality of the niobium as a superconductor \cite{NiobiumOxideTc}. Also note that the thicknesses of the niobium ground plane and microstrip are 300 and 600 nm, respectively. At this thickness, both films are significantly thicker than the penetration depth \cite{NiobiumPenetrationDepth}, which reduces the dependence of the microstrip impedance on film thickness \cite{YassinWithingtonMicrostrip}.

\subsubsection{Bolometer structure}

The \ac{tes} thermistor is an aluminum/titanium bilayer. The layers of this bilayer must be deposited in direct succession without the wafer leaving the vacuum environment, and with an etch before the deposition to remove the native niobium oxide on the niobium thermistor leads. The gold used to increase the heat capacity of the bolometer is deposited using electron beam evaporation onto already patterned photoresist, so that the gold is removed everywhere other than where it is desired by dissolving the photoresist in a lift-off procedure. The \ac{tes} bolometer suspension structure is provided by the Low-stress Silicon Nitride (LSN) film deposited directly on the silicon wafer. This deposition is first, but the etching of the LSN occurs much later in the process, after all of the deposition steps are complete. To release the suspended \ac{tes} structure, trenches are etched into the LSN, and then the silicon substrate underneath the LSN is removed by an isotropic gaseous xenon difluoride ($\mbox{XeF}_2$) etch. The rough surface below the released bolometer (visible in Figure~ \ref{fig:pixel}\textbf{(c)}) is where the silicon has been etched by the $\mbox{XeF}_2$, leaving the suspended bolometer structure.

\section{FOCAL PLANE PERFORMANCE VALIDATION}
\label{sec:valid}

This section reports on several of the tests that were conducted to validate the \textsc{Polarbear}\ detector performance.

\subsection{Pixel beam maps}

Maps of the detectors' spatial response patterns were measured in the laboratory using a modulated thermal source on translation stages. Figure \ref{fig:beammaps} shows the result of that measurement for each polarization of a single pixel. Note that the beam ellipticity is small (best fit elliptical Gaussian beams to this data have ellipticities of 1\%), and that the differential pointing is small, with a best fit value of 0.3\% of the FWHM.

   \begin{figure}[htpb]
   \begin{center}
	\begin{tabular}{cc}
   \includegraphics[width=2.5in]{b_c.eps} & \includegraphics[width=2.5in]{t_c.eps} \\
{\bf(a)} & {\bf(b)}
\end{tabular}
	\end{center}
	\caption[Pixel beam maps]{\label{fig:beammaps}Spatial beam map of the \textsc{Polarbear}\ pixel response for each polarization, as measured using a chopped thermal source on a translation stage.}
	\end{figure}

\subsection{Boresight cross-polar response}

A modulated thermal source was also used to measure the pixel's boresight cross-polar response. The thermal source was placed at the location of peak spatial response, and a lithographed wire grid polarizer on a thin dielectric sheet was placed between the thermal source and the detector. The response of both bolometers in the pixel was recorded as the polarizer was rotated. A model was fit to the data of the form $c_0 +\left ( c_1-c_0 \right ) \cos^2\left (\theta-c_3\right )$. The cross-polar response for each bolometer is defined as $c_0/c_1$, and was found to be less that 1\%. Figure \ref{fig:polpurity}\textbf{(a)} shows one such measurement.

   \begin{figure}[htpb]
   \begin{center}
	\begin{tabular}{cc}
   \includegraphics[width=2.9in]{singlePixelPol.eps} & \includegraphics[width=3in]{SPIESpec.eps} \\
\textbf{(a)} &\textbf{(b)}
\end{tabular}
	\end{center}
	\caption[Polariztion purity]{\label{fig:polpurity} \textbf{(a)} Boresight cross-polar response measurement of the individuel detectors in a \textsc{Polarbear}\ pixel. The figure shows data and a fit sinusoidal model with parameters for amplitude, phase, and cross-polar response. The fit cross-polar response is less than 1\% for each bolometer. \textbf{(b)} The spectral response of each bolometer for a single \textsc{Polarbear}\ pixel, plotted with the design spectral band and a model of the atmospheric transmittance in Chile for 1 mm of precipitable water vapor. The spectra are well-matched to each other and to the design band.}
	\end{figure}

\subsection{Spectral response}

The spectral response of the \textsc{Polarbear}\ detectors was measured in the lab using a \ac{fts}. Repeatable spectral measurements required that the spatial beam of the detector be filled by the output of the spectrometer. This was achieved using a warm UHMWPE lens to match the spectrometer output to the optics of the observation cryostat. Figure \ref{fig:polpurity}\textbf{(b)} shows the spectral band measurement for one of the pixels in the \textsc{Polarbear}\ focal plane array, and Table~ \ref{tab:spectra} enumerates the average spectral band parameters for each of the wafers in the array. Six of the seven average spectral bands are within the goal 5\% of band central location and integrated bandwidth.

\begin{table}
\vspace{6pt}
\begin{center}
\begin{tabular}{c|cc}
\textbf{Wafer} & \textbf{\shortstack{Center\\frequency\\(GHz)}} & \textbf{\shortstack{Integrated\\bandwidth\\(GHz)}} \\
\hline
Design & 148 & 38 \\
8.2.0 & 139 & 31 \\
9.4 & 148 & 36 \\
10.1 & 141 & 37 \\
10.1 & 145 & 38 \\
10.1 & 150 & 35 \\
10.4 & 145 & 36 \\
10.5 & 146 & 36
\end{tabular}
\end{center}
\caption[Spectral properties]{\label{tab:spectra}Per-wafer average measured center frequencies and integrated bandwidths. Design values are 148 GHz and 38 GHz, respectively.}
\end{table}

\subsection{Bolometer yield}

An important factor in total array sensitivity is the number of bolometers that are functional in the deployed array. Of the 1274 antenna-coupled bolometers in the \textsc{Polarbear}\ focal plane, 1015 show nominal optical response to a planet, an optical yield of 80\%. The factors that affect this number are shown in Table~ \ref{tab:yield}.

\subsection{Bolometer optical time constants}

The detector optical time constants were measured under observing conditions using a modulated thermal source that is instrumented through a 6 mm hole in the secondary mirror. A single pole time constant model is a good fit to the detectors' response functions, and the uncertainty in the fit for most detectors is less than 0.2 ms. A histogram of the measured optical time constants, with a peak value around 1.9 ms, is shown in Figure~ \ref{fig:electrothermal}. This is fast enough to measure the changing \ac{cmb} signal as the \textsc{Polarbear}\ telescope scans across it.

   \begin{figure}[htpb]
   \begin{center}
   \includegraphics[width=3.5in]{tau.eps}
	\end{center}
	\caption[Electrothermal properties]{\label{fig:electrothermal} Detector optical time constants as measured during operation. Peak value of 1.9 ms corresponds to an optical bandwidth of 84 Hz, and is fast enough to measure the changing \ac{cmb} signal as the \textsc{Polarbear}\ telescope scans across it.}
	\end{figure}

\begin{table}
\begin{center}
\begin{tabular}{lr}
\textbf{Functional check} & \textbf{Number} \\
\hline
Total fabricated (1274 antenna-coupled, 70 ``dark'') & 1344 \\
Passed warm inspection & 1245 \\
Electrically responsive cold & 1154 \\
Biasable and shows nominal beammaps & 1015
\end{tabular}
\vspace{6pt}
\caption[Bolometer yield]{\label{tab:yield}Bolometer yield in the deployed \textsc{Polarbear}\ focal plane. Of the 1274 antenna-coupled bolometers (70 are not coupled to antennas and are used as checks for spurious signals), $1015/1274=80$\% show nominal beammaps.}
\end{center}
\end{table}

\section{THE FUTURE OF POLARBEAR AND LENSLET-COUPLED DETECTORS}
\label{sec:future}

\textsc{Polarbear}\ is currently observing from Cerro Toco in Northern Chile, which marks the \textsc{Polarbear}\ detector array as the first lenslet-coupled detectors to observe the \ac{cmb}. This is an exciting advance, both because lenslet-coupling proves a path for a multi-band pixel architecture that will be an enabling technology for the next generation of \ac{cmb} experiments\cite{Suzuki_SPIE2012}, and because the current \textsc{Polarbear}\ instrument will improve our understanding of cosmology with the data it collects in the coming years.

\acknowledgments

The \textsc{Polarbear}\ project is funded by the National Science Foundation under grant AST-0618398. Antenna-coupled bolometer development at Berkeley is also funded by NASA under grant NNG06GJ08G. The KEK authors were supported by MEXT KAKENHI Grant Number 21111002. The McGill authors acknowledge funding from the Natural Sciences and Engineering Research Council and Canadian Institute for Advanced Research. MD acknowledges support from an Alfred P. Sloan Research Fellowship and Canada Research Chairs program. All silicon wafer-based technology is fabricated at the UC Berkeley Nanolab.


\section*{\refname\@mkboth{\refname}{\refname}}%
  \list{\@biblabel{\@arabic\c@enumiv}}%
       {\settowidth\labelwidth{\@biblabel{#1}}%
        \leftmargin\labelwidth
        \advance\leftmargin\labelsep

 \setlength\baselineskip{11pt}%
        \@openbib@code
        \usecounter{enumiv}%
        \let\p@enumiv\@empty
        \renewcommand\theenumiv{\@arabic\c@enumiv}}%
  \sloppy
  \clubpenalty4000
  \@clubpenalty\clubpenalty
  \widowpenalty4000%
  \sfcode`\.\@m}
 {\def\@noitemerr
 {\@latex@warning{Empty `thebibliography' environment}}%
\endlist}
\makeatother
\begin{document}
\centerline{{\sl Genshikaku Kenkyu Suppl.} No. 000 (2012)
\begin{center}
\vskip 2mm
{\Large\bf

 Measurements of the photon induced production of $\Lambda$ in the $^2$H$({\gamma},{\Lambda})$X process at threshold energies.

\hspace{-1mm}\footnote{Presented at the International Workshop on Strangeness
Nuclear Physics (SNP12), August 27 - 29, 2012, \\
\hspace*{5mm} Neyagawa, Osaka, Japan.}
}\vspace{5mm}

{
B. Beckford$^a$,  P. Byd$\breve{z}$ovsk$\acute{y}$$^{d}$, A. Chiba$^a$, D. Doi$^a$, T. Fujii$^a$, Y. Fujii$^a$, K. Futatsukawa$^a$,\\ T. Gogami$^a$, O. Hashimoto$^a$,  Y.C. Han$^c$ ,K. Hirose$^b$, S. Hirose $^a$,  R. Honda$^a$, K. Hosomi$^a$, \\T. Ishikawa$^b$,  H. Kanda$^a$, M. Kaneta$^a$, Y. Kaneko$^a$, S. Kato$^a$, D. Kawama$^a$,  C. Kimura$^a$\\ S. Kiyokawa$^a$, T. Koike$^a$,  K. Maeda$^a$, K. Makabe$^a$, M. Matsubara$^a$, K. Miwa$^a$, S. Nagao$^a$,\\ S. N. Nakamura$^a$, A. Okuyama$^a$, K. Shirotori$^a$,  K. Sugihara$^a$, K. Suzuki$^b$, T. Tamae$^b$, H. Tamura$^a$,\\ K. Tsukada$^a$, K. Yagi$^a$,  F. Yamamoto$^a$, T. O. Yamamotoi$^a$, H. Yamazaki$^a$, and Y. Yonemoto$^a$
}\bigskip

{\small
$^a$Department of Physics, Tohoku University,
Sendai 9808578, Japan \\
$^b$Research Center for Electron Photon Science, Tohoku University, Sendai,
982-0826\\
$^c$School f Nuclear Science and Technology, Lanzhou University, Lanzhou,
730000 \\
$^d$Nuclear Physics Institute, $\breve{R}$e$\breve{z}$ near Prague, Czech Republic,
25068 \\

}
\end{center}
\vspace{3mm}

\noindent
{\small \textbf{Abstract}:\quad
An experiment  was carried out with the NKS2+ in 2010 at the Research Center for Electron Photon Science (ELPH), in which tagged photon beams in the range of 0.8 $\le$ $E_{\gamma}$ $\le$ 1.1 GeV were impinged on a liquid $^2$H target positioned at the center of the NKS2+. The produced $\Lambda$ was subsequently detected by the $p{\pi^{-}}$ decay channel.
Integrated cross sections of the $^2$H$({\gamma},{\Lambda})$X in the angular region of 0.9 $\le$ $\cos{\theta}_{\Lambda}^{LAB}$ $\le$ 1.0 was derived and compared with preceding experimental results of the NKS2 collaboration.  In addition, the momentum spectra for two photon energy regions were also procured.  The ${\Lambda}$ angle dependent cross sections as a function of the scattering angle in the laboratory system was additionally deduced. We present the latest results on the excitation function of ${\Lambda}$ photoproduction, the momentum distributions, and polarization.
}%

\section{Introduction}
The photoinduced  production of strangeness has been tirelessly investigated on the proton and most recently the emphasis has veered to the neutron to actualize a comprehensive model capable of describing the reaction mechanism. To facilitate this a deuteron target is used to provide a loosely bound neutron target. Our focus is directed in the threshold energy region, in which it is conceded that the reaction will be significantly less disturbed by higher nucleon resonances, allowing for a simplification in the description of the reaction and thus permitting the study with less uncertainty. The sizable experimental data sets that have been measured for the $^1$H$({\gamma},K^{+}) {\Lambda}$ reaction at facilities such as CLAS and SAPHIR, are not adequate to constrain theoretical models and successfully predict the cross section of the unmeasured channels. The goal of this work is to further investigate strangeness photoproduction in the threshold region by focusing on $\Lambda$ production.

\section{The NKS2+ Spectrometer }
The upgraded large acceptance magnetic spectrometer, NKS2+, is composed of a cryogenic target system located in the center of various detector systems that works in tandem with a photon tagging array.  Using a carbon wire as a target, photon beams are produced via the electromagnetic interaction with the electrons orbiting in the Stretcher Booster Ring (STB) in the form of bremsstrahlung radiation. The STB ring is capable of accelerating electrons injected from a 0.15 GeV Linac injector up to 1.2 GeV.
The photon beam is guided  through a collimator in order to reduce the beam halo, a sweep magnet, and into the NKS2+ and is bombarded on the target. Moving from the inner most position outwards, is the target which is surrounded by a Vertex Drift Chamber (VDC) and an Inner Hodoscope (IH), which acts as the start trigger for time of flight measurements. These detectors are enclosed in a Cylindrical Drift Chamber (CDC). Substantial effort was placed on advancing the study on the neutron as a consequence of the favorable findings of previous experiments with the Neutral Kaon Spectrometer (NKS) and NKS2~\cite{takahashi,watanabe,Kanda,tsukada}; The spectrometer was upgraded by the replacement of the inner detectors in order to further improve upon the acceptance. An experiment has been designed and performed by the NKS2 collaboration was taken at ELPH in 2010.

  \begin{figure}[ht]
  	  \vspace{-.6cm}
    \begin{center}
     \includegraphics[width=10.0cm]{NKS2_upgraded_Schematic_View-eps-converted-to.pdf}
       \vspace{-.4cm}
      \caption{
        The figure shows a schematic perspective the of recently improved NKS2 spectrometer, henceforth referred to as NKS2+.  The inner detector package has been  redesigned and installed and is visible in the zoomed image on the right.
      }
      	  \vspace{-.9cm}
      \label{fig:NKS2_upgraded_detectors}
    \end{center}
  \end{figure}

  \begin{figure}[h!]
\begin{center}
	      \includegraphics[width=7.15cm]{pid_mom_inv_beta.pdf}
	      	\includegraphics[width=8.75cm]{lambda_im_fitted.pdf}
     	 	\caption{The particle momentum multiplied with the particle charge plotted as a function of the inverse velocity is shown on the left.  The proton and pion selection regions are shown n the curved bands. The charged pions are identified by the sign of the momentum. The $p{\pi^{-}}$ invariant mass resolution for the energy bin of 0.9 $\le E_\gamma \le$1.1 GeV (right) was found by fitting as 2.87 $\pm$ 0.19 MeV/c in rms.  }
     	 \label{fig:momentum_versus_inverse_beta}
 \vspace{-.6cm}
\end{center}
\end{figure}

  \section{Data Analysis}
The experiment was performed with tagged photon beams in the range of 0.8 $\le$ $E_{\gamma}$ $\le$ 1.1 GeV, at tagging rate of $1.5 - 2.5$ MHz. The accepted number of events and recorded photons were 0.64 $\times$ 10$^{9} $ and  0.89 $\times$ 10$^{12}$ respectively. The direct approach to measuring the photo production of strangeness is by reconstructing the invariant mass of the produced particles that contain a strange quark. The momentum determination and tracking of the detected particles was performed by the CDC and the particle species was identified by Time Of Flight (TOF) measurements. . The PID capabilities of the NKS2+ is showed in Figure~\ref{fig:momentum_versus_inverse_beta}, where the particle's momentum multiplied with the charge plotted as a function of the inverse velocity is presented. Additionally, an opening angle selection requirement of $-09$ $\le$ cos${\theta}$ $\le$ 0.9 was used in order to reduce the $e^{+}e^{-}$ contamination.  The produced $\Lambda$ was detected in the $p{\pi^{-}}$ decay channel.
Numerous cuts were applied to the obtained raw invariant mass spectrum to resolve and extract a clean $\Lambda$ distribution from the continuum shown on the right in Figure~\ref{fig:momentum_versus_inverse_beta}.

\section{Results and Discussion}
We next proceed to employ various theoretical calculations as a comparison medium for the experimentally yielded outcomes. This initial comparison is made with the momentum dependent differential cross section obtained for the inclusive measurement of ${\Lambda}$ to the Kaon-MAID (KM), Saclay-Lyon A (SLA) and Regge-plus-resonance (RPR-2007) models. The comparison of results with KM and SLA calculations for the addition of amplitudes of the $^2$H$({\gamma},K^{+}) {\Lambda}n$ and $^2$H$({\gamma},K^{0}) {\Lambda}p$ reactions for incident photon beam energy of (top) $0.9 - 1.0$ and (bottom) $1.0 - 1.08$ GeV in the laboratory frame.  The angular regions of  $0.95 - 1.0$ and $0.90 - 0.95$ (cos$\theta_{\Lambda}^{LAB}$)  are shown left to right respectively in the left two columns of Figure~\ref{fig:Comparison_SLA_KM_Lambda_singles}. The $r_{K_{1}K\gamma}$ coupling constant, which is a parameter for the SLA model, is displayed in the figure.  The RPR-2007 comparisons are seen in the right two columns. SLA $r_{K_{1}K\gamma}$ = $-(1.0 -1.5)$ and RPR-2007 adequately reproduce the size and shape of the data, especially in the higher energy bin.
      \begin{figure}[ht]
   	 \begin{center}
	 \includegraphics[width=8.0cm]{SLA_KM_Lambda_singles_SUM_low_comparison.pdf}
	  \includegraphics[width=8.0cm]{Regge_Plus_Lambda_singles_comparison_low_no_error.pdf}
	 \includegraphics[width=8.2cm]{SLA_KM_Lambda_singles_SUM_high_comparison.pdf}
         \includegraphics[width=8.2cm]{Regge_Plus_Lambda_singles_comparison_high_no_error.pdf}
	  \vspace{-.6cm}
      \caption{
Comparison of results with KM and SLA calculations of the semi-inclusive ${\Lambda}$ momentum cross sections for the addition of contributions of the $^2$H$({\gamma},K^{+}) {\Lambda}n$ and $^2$H$({\gamma},K^{0}) {\Lambda}p$ reactions  for photon beam energy of  (top) $0.9 - 1.0$ and (bottom) $1.0 - 1.08$ GeV in the laboratory frame.  The angular integration regions of $0.95 - 1.0$ and $0.90 - 0.95$ (cos$\theta_{\Lambda}^{LAB}$) are shown left to right respectively. The $r_{K_{1}K\gamma}$ parameter for the SLA models is displayed in the figure. The RPR comparisons are given in the right two columns. }
      \label{fig:Comparison_SLA_KM_Lambda_singles}
    \end{center}
  \end{figure}

The experimental results of the inclusive measurement of ${\Lambda}$ on a deuterium target using the NKS2 spectrometer prior to the recent upgrade was successfully carried out and the results were reported~\cite{futatsukawa,Beckford_Baryon10}.
The background subtracted energy dependent cross section of the inclusive measurement of ${\Lambda}$ photoproduction in $^2$H$({\gamma},{\Lambda})$X reaction for the angular integrated region of 0.9$\le$cos${\theta}_{\Lambda}^{LAB}$ $\le$1.0 for the NKS2~\cite{futatsukawa,Beckford_Baryon10} and the latest results, obtained from independent experiments and analysis, are presented in Figure~\ref{fig:total_gd_KpK0_L_SLA_comparison} as the solid circles and open squares respectively. Only statistical errors are shown. The latest results furnished an excitation function that has considerable agreement with the finalized results of the preceeding experimental, emphasizing the uniqueness of the research, and more importantly, the reliability of the data reported by the NKS2/NKS2+ collaboration.
The SLA predictions describe well the data for $r_{K_{1}K\gamma}$ = $-1.4$ shown on the right in Figure~\ref{fig:total_gd_KpK0_L_SLA_comparison}.

\begin{figure}[h!]
\begin{center}
\vspace{-.5cm}
        \includegraphics[width=8.15cm]{EnergyBins_Distributions_Futat_comp_cos_09-1.pdf}
        \includegraphics[width=8.15cm]{total_gd_KpK0_L_KM_SLA_comparison_09-1.pdf}
	\caption{The results of the NKS2~\cite{futatsukawa} measurements of the energy dependent integrated cross section are shown with the NKS2+ results as open squares and black circles respectively (left). Comparison to integrated cross section as a function photon energy predicted by the SLA and KM theoretical models (right).  Where the cross sections of $^2$H$({\gamma},K^{+}) {\Lambda}n$ and $^2$H$({\gamma},K^{0}) {\Lambda}p$ reactions are summed.  }
\vspace{-.3cm}
 \label{fig:total_gd_KpK0_L_SLA_comparison}
\end{center}
\end{figure}

We next report on our measurements of the ${\Lambda}$ recoil polarization observables.
The self analyzing nature of the weak decay of the  ${\Lambda}$ hyperon gives rise to not only parity violation but also an asymmetry.
The polarization,${P_{\Lambda}}$, was computed from the decay asymmetry of the angular distributions and is calculated by Equation~\ref{eq:lambda_polarization}.

  \begin{equation}
	{P_{\Lambda}}  = \frac{2}{{\alpha} }\frac{N_1 - N_2} {N_1 + N_2 }
		\label{eq:lambda_polarization}
  \end{equation}

Where ${\alpha}$ = 0.6421 $\pm$ 0.013, is the decay asymmetry parameter, $N_1$ and $N_2$ are the total number of events where cos${\theta}$ is greater or less than zero respectively.
It this framework cos${\theta}$ denotes the angle that lies between the proton from the $\Lambda$ decay and the normal of the reaction plane.
This asymmetry, typically thought of as the ratio between the difference and the addition of the competing patterns is a consequence of the witnessed interference between the \emph{s} and \emph{p} waves.

  \begin{figure}[ht]
   \begin{center}
    \includegraphics[width=5.25cm]{polarization_momentum_EG_900-1000.pdf}
    \includegraphics[width=5.25cm]{polarization_momentum_EG_950-1050.pdf}
    \includegraphics[width=5.25cm]{polarization_momentum_EG_1000-1100.pdf}
      \caption{${\Lambda}$ recoil polarization as a function of  momentum in the laboratory frame for the energy range of $0.9 - 1.0, 0.95 - 1.05$,  and $1.0 - 1.1$ GeV.
      }
        \vspace{-.5cm}
      \label{fig:polarization_of Lambda_mom_E_ranges}
    \end{center}
  \end{figure}

The momentum dependence of ${\Lambda}$ polarization was determined for three photon energy ranges,  $0.9 - 1.0, 0.95 - 1.05$,  and $1.0 - 1.1$ GeV, in the laboratory frame. The results of polarization of the energy ranges are shown in Figure~\ref{fig:polarization_of Lambda_mom_E_ranges},  left to right respectively.  The polarization in the lower energy bin is distinctively
negative. Whereas, in the other energy bins, the value of the polarization is consistent with zero.

\begin{figure}[h!]
\begin{center}
	\includegraphics[width=8.15cm]{total_k0_k+_exp_data_comp.pdf}
	\includegraphics[width=8.15cm]{total_k0_k+_exp_data_comp_theory.pdf}
	\caption{Comparison of the estimated total $^2$H$({\gamma},K^{0}) {\Lambda}p$ cross section with theoretical predictions and published K$^{+}\Lambda$ and K$^{0}\Lambda$ data using proton and deuteron targets~\cite{SAPHIR_data,SAPHIR_data_98,CLAS_data, futatsukawa_dissertation, bantawa_phd}. Only the statistical errors are shown. }
	\label{fig:total_exp_data_comparison_theory}
	   \vspace{-.8cm}
\end{center}
\end{figure}

\begin{equation}
\sigma_{\mathrm {Total}}^{K^{0}\Lambda} =( \sigma_{\mathrm{Integrated}}^{\gamma d\rightarrow \Lambda X} - \tilde{\sigma}_{\mathrm{Integrated}}^{'\gamma d\rightarrow K^{+}\Lambda n}) \cdot \frac{\tilde{\sigma}_{\mathrm{Total}}^{'\gamma d \rightarrow K^{0}\Lambda p} }{\tilde{\sigma}_{\mathrm{Integrated}}^{'\gamma d \rightarrow K^{0}\Lambda p}}
\label{eq:k0_total_estimation_equation}
\end{equation}

The$^2$H$({\gamma},K^{0}) {\Lambda}p$ total cross section was extrapolated by equation~\ref{eq:k0_total_estimation_equation}.
In the equation, the total and integrated cross of the $\gamma d \rightarrow K^0 \Lambda p$ process are denoted as  $\tilde{\sigma}_{Total}^{'\gamma d \rightarrow K^{0}\Lambda p}$ and $\tilde{\sigma}_{Integrated}^{'\gamma d \rightarrow K^{0}\Lambda p}$. It follows that the integrated cross of the $\gamma d \rightarrow K^+ \Lambda n$ process is designated as $\tilde{\sigma}_{Integrated}^{'\gamma d\rightarrow K^{+}\Lambda n}$. In the above calculation, all integrated cross sections were derived  for the angular integration range of 0.9 $\le$ cos$\theta_{\Lambda}$ $\le$ 1.0 in the laboratory frame of reference.  The values introduced into the calculation for $\tilde{\sigma}_{Total}^{'\gamma d \rightarrow K^{0}\Lambda p}$, $\tilde{\sigma}_{Integrated}^{'\gamma d \rightarrow K^{0}\Lambda p}$ and $ \tilde{\sigma}_{Integrated}^{'\gamma d\rightarrow K^{+}\Lambda n}$ were procured by purely theoretical means.
The systematic error associated with this method  is approximated at less than 15\%. There is a minor model dependence on the extracted total cross section deduced by using either theoretical calculations.

Hence, the estimated total cross section derived with the KM and SLA  $r_{K_{1}K\gamma}$= $-1.5$ can be juxtaposed to experimental data. The published results by SAPHIR~\cite{SAPHIR_data, SAPHIR_data_98} and CLAS~\cite{CLAS_data}, and those taken from the dissertations of K. Bantawa\cite{bantawa_phd} and K. Futatsukawa\cite{futatsukawa_dissertation} are plotted alongside the estimated total cross section of $^2$H$({\gamma},K^{0}) {\Lambda}p$ by the method described in equation~\ref{eq:k0_total_estimation_equation} in Figure~\ref{fig:total_exp_data_comparison_theory}.
The maximum photon energy range of the NKS2+ data analyzed in this work is $E_{\gamma}$= 1.08 GeV. While results of the other collaborations extended into a higher energy regimes. The total $^2$H$({\gamma},K^{0}) {\Lambda}p$ compare favorably well to the measured $^1$H$({\gamma},K^{+}) {\Lambda}$ energy dependent cross sections up to 1.08 GeV, this implies that cross sections are roughly on par in terms of overall magnitude. Despite the general agreement between our semi-inclusive energy dependent cross section to those of the SAPHIR and CLAS results, the MAINZ results exhibits a different trend. Our results over the integrated energy range of $0.9 - 1.08$ GeV has a calculated value of $40 - 20$\% higher than the results reported by their group.

Model predictions for the inclusive $^2$H$({\gamma},K^{+}) {\Lambda}n$ and $^2$H$({\gamma},K^{0}) {\Lambda}p$, calculated in the KM and SLA framework, are compared to the experimental results obtained in the threshold energy region gathered by the NKS2, NKS2+, CLAS, SAPHIR,  and MAINZ collaborations. These comparison are shown in Figure~\ref{fig:total_exp_data_comparison_theory} where the curves for the $^2$H$({\gamma},K^{+}) {\Lambda}$n process predicted by KM and SLA and drawn as solid and dashed black lines.  Calculations of the $^2$H$({\gamma},K^{0}) {\Lambda}p$ for KM and SLA,$r_{K_{1}K_\gamma} $ of $-(1.0 - 2.0)$, are shown as the solid light red, dashed light red, dashed, light blue, and dashed light purple curves respectively. The KM  and SLA predictions of the $^2$H$({\gamma},K^{+}) {\Lambda}$n process are reasonably similar and provides an great description of the data up to $E_{\gamma} \le $1.0 GeV, where after the both predictions slightly under shoot the data between 1.0 $E_{\gamma} \le $1.1 GeV and then return to a consistent agreement.  The general consensus between our experimental results and theoretical predictions of the SLA model is good. The shape and size of the extracted total $^2$H$({\gamma},K^{0}) {\Lambda}p$ cross section is acceptably reproduced for the $r_{K_{1}K_\gamma} $ value around $-1.4$.

\section{Acknowledgements}
 \vspace{-.3cm}
Theoretical calculations for RPR-2007 model were made by P. Vancraeyvald~\cite{RPR_private} of Ghent University.
The research was conducted under the patient and indispensible guidance of the late Professor Osamu Hashimoto.  May he find peace and know that his legacy continues on.

  \vspace{-.2cm}