full_reportE_appendices.tex

\appendix

\clearpage

\iftoggle{debuggingaidsappendix}{
\section{Debugging aids}

If the $\chi^2$ plots in Figures \ref{bestEstimateColorBoxPlot}-\ref{invertedBestEstimateColorBoxPlot} indicate that the data is poorly fit by GST, the next step is to begin ``debugging'' the experiments and/or the fit.  Most commonly, a poor fit is due to non-Markovian behavior.  However, there are many kinds of non-Markovian behavior.  The most straightforward occurs when the gate set fluctuates over time, or when there is other time correlation in the experiments (e.g., due to memory effects).  However, another possibility that must be considered is that repeated gate operations cause changes in the system, e.g. heating it up (as is seen in 2-qubit trapped-ion gates) so that the data from long gate sequences is simply chaotic and inconsistent with shorter experiments.

Figure \ref{directLSGSTChi2BoxPlot} provides a test (albeit currently an unreliable one) for such an effect.  Like Figures \ref{bestEstimateColorBoxPlot}-\ref{invertedBestEstimateColorBoxPlot}, it displays per-experiment $\chi^2$ values -- but \emph{not} for any single gate set.  Instead, this \emph{direct GST} analysis treats each base sequence as an independent process (\emph{not} as a product of many gates), and analyzes it using LSGST together with the individual gates (which are necessary to model the effect of the fiducial sequences that precede and follow the base sequence being analyzed).  The resulting direct GST estimate is then used to assign probabilities for the corresponding experiments.  

This analysis decouples the various germs and base sequences from each other.  Unlike a standard GST analysis, it does not model different base sequences as being generated from the same gates.  Therefore, this analysis \emph{should} be fairly consistent (i.e., lots of gray/white squares and few or no red ones) even if Figures \ref{bestEstimateColorBoxPlot}-\ref{invertedBestEstimateColorBoxPlot} indicate severe fit problems.
%However, beware one caveat.  Because this analysis uses \emph{only} LGST, it is not aware of $\chi^2$ (unlike the full GST analysis, which explicitly seeks to minimize $\chi^2$).  In particular, it can and often does assign slightly negative probabilities to experiments with frequencies close to zero.  This produces implausibly high $\chi^2$ values for those experiments.  Therefore, we recommend ignoring sporadic red squares in this analysis.

\begin{figure}
\begin{center}
\includegraphics[width=6.50in,height=8.00in,keepaspectratio]{full_reportE_files/directLongSeqGSTChi2Boxes.pdf}
\caption{\pdftooltip{{\color{blue}\texttt{.figures['directLongSeqGSTColorBoxPlot']}}\quad}{Access this information in pyGSTi via<ResultsObject>.figures['directLongSeqGSTColorBoxPlot']}\textbf{$\chi^2$ values for ``direct GST'' fit.}  This plot indicates how well direct GST analysis of each base sequence fits the observed data.  By decoupling different base sequences, this analysis largely avoids making any assumptions about Markovianity, and therefore serves as a sanity check on the full GST analysis.  See text for details.\label{directLSGSTChi2BoxPlot}}
\end{center}
\end{figure}

Figure \ref{directLSGSTDeviationBoxPlot} uses direct GST analysis in another way, but to detect similar effects.   The goal here is to to compare two different predictions for each base sequence:  (1) the one given by the overall GST fit; and (2) the one obtained by direct GST on that base sequence, as in Figure \ref{directLSGSTChi2BoxPlot}.  The analysis shown in Figure \ref{directLSGSTDeviationBoxPlot} does so by computing each process's fidelity with the closest unitary.  This is a measure of non-unitary decoherence.  Then, it plots the \emph{difference} between the values of this ``unitarity'' for (1) the overall GST fit, and (2) the direct GST estimate.  The direct GST estimate is not very precise, but it is extremely reliable, because it is not influenced by any data that are not \emph{directly} connected with that base sequence.  Thus, this serves as another sanity check and debugging aid.  Large absolute values indicate that the full GST fit is signicantly inconsistent with the data (complementing the analysis in Figure \ref{bestEstimateColorBoxPlot}).  Positive numbers indicate that full GST has \emph{overestimated} the amount of decoherence, while negative numbers indicate that it \emph{underestimated} it.

\begin{figure}
\begin{center}
\includegraphics[width=4.00in,height=5.00in,keepaspectratio]{full_reportE_files/directLongSeqGSTDeviationBoxes.pdf}
\caption{\pdftooltip{{\color{blue}\texttt{.figures['directLongSeqGSTDeviationColorBoxPlot']}}\quad}{Access this information in pyGSTi via<ResultsObject>.figures['directLongSeqGSTDeviationColorBoxPlot']}\textbf{Inconsistency of unitarity between GST and direct GST.}  This plot shows, for each base sequence, the increase in ``upper bound of fidelity with unitary'' (see Table \ref{bestGatesetClosestUnitaryTable}) when using the direct-GST result for a gate sequence instead of the process given by the best gate set.\label{directLSGSTDeviationBoxPlot}}
\end{center}
\end{figure}

Finally, for various reasons including the diagnosis of non-Markovian behavior, it is often useful to have a direct and reliable estimate of the per-gate incoherent error rate observed in each base sequence.  This is shown in Figure \ref{smallEigvalErrRateColorBoxPlot}.  To obtain these numbers, direct GST is performed on each base sequence.  Then, each resulting process matrices is diagonalized, and the smallest eigenvalue (corresponding the most rapid loss of information/coherence/polarization) is extracted.  This eigenvalue is then raised to the $1/L$ power, where $L$ is the length of the base sequence, to estimate the rate of decoherence per gate, and subtracted from 1 (to convert it to an error rate).  These numbers become much more reliable towards the right-hand side of the plot, because errors in LGST become far less significant for these long sequences.  Large changes in this direct-GST error rate as $L$ is changed (on any given row of the table) are a ``smoking gun'' for non-Markovian decoherence -- especially when the error rate \emph{decreases} with increasing $L$.  Ideally, all numbers in a given row should be the same.

\begin{figure}
\begin{center}
\includegraphics[width=4.00in,height=5.00in,keepaspectratio]{full_reportE_files/smallEigvalErrRateBoxes.pdf}
\caption{\pdftooltip{{\color{blue}\texttt{.figures['smallEigvalErrRateColorBoxPlot']}}\quad}{Access this information in pyGSTi via<ResultsObject>.figures['smallEigvalErrRateColorBoxPlot']}\textbf{Per-gate error rates, extrapolated from the smallest magnitude eigenvalue of the direct GST gate matrix.}  See text for further details.\label{smallEigvalErrRateColorBoxPlot}}
\end{center}
\end{figure}

%\begin{figure}
%\begin{center}
%\XXputfield{bestEstimateColorBoxPlot_hist}{Histogram of best gate set chi2}
%\caption{\XXputfield{tt_bestEstimateColorBoxPlot_hist}{}Histogram plot of best gate set $\chi^2$ values.}
%\end{center}
%\end{figure}
}{}

\iftoggle{gaugeoptappendix}{
\section{Best gate set in different gauges\label{appendix_gauge_opts}}

In this appendix, we report the non-gauge-invariant quantities from the main text for different gauge choices.  In each section, GST's best gate set estimate is gauge optimized according to a different objective function, as specified within the section.

\subsection{Gauge optimization to the target gate set}
The gauge is chosen to minimize the Frobenius distance between the estimate and the target gates, with equal weight given to the gate and SPAM parameters.

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
Operator & Matrix & Hilbert-Schmidt vector (Pauli basis) \\ \hline
$\rho_{0}$ & $ \left(\!\!\begin{array}{cc}
0.961 & 0.001e^{-i0.13} \\ 
0.001e^{i0.13} & 0.047
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.712 \\ 
0.001 \\ 
1\e{-4} \\ 
0.646
 \end{array} $
 \\ \hline
$E_{0}$ & $ \left(\!\!\begin{array}{cc}
0.048 & 0.001e^{i1.79} \\ 
0.001e^{-i1.79} & 0.962
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.714 \\ 
-4\e{-4} \\ 
-0.002 \\ 
-0.646
 \end{array} $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the SPAM operations}.  Compare to Table \ref{targetSpamTable}.\label{bestTargetGatesetSpamTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
 & $E_{0}$ & $E_C$ \\ \hline
$\rho_{0}$ & 0.091 & 0.909 \\ \hline
\end{tabular}

\caption{\textbf{GST estimate of SPAM probabilities}.  Computed by taking the dot products of vectors in Table \ref{bestTargetGatesetSpamTable}.\label{bestTargetGatesetSpamParametersTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|}
\hline
Gate & Superoperator (Pauli basis) \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
1 & 3\e{-4} & 4\e{-4} & -3\e{-4} \\ 
-0.001 & 0.945 & -0.006 & -0.009 \\ 
4\e{-4} & 0.001 & 0.947 & 0.004 \\ 
0.002 & 2\e{-5} & 0.001 & 0.949
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
1 & 9\e{-5} & 0.007 & -0.007 \\ 
4\e{-4} & 0.951 & -0.005 & -0.005 \\ 
0.01 & 0.001 & 2\e{-4} & -0.949 \\ 
0.011 & 0.001 & 0.949 & -0.002
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
1 & -0.007 & -4\e{-5} & -0.007 \\ 
-0.01 & 5\e{-5} & -0.004 & 0.949 \\ 
3\e{-4} & 0.001 & 0.951 & 0.001 \\ 
0.01 & -0.949 & -0.006 & -0.001
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the logic gate operations}.  Compare to Table \ref{targetGatesTable}.\label{bestTargetGatesetGatesTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|}
\hline
Gate & \begin{tabular}{c}Process\\Infidelity\end{tabular} & \begin{tabular}{c}$\nicefrac{1}{2}$ Trace\\Distance\end{tabular} & $\nicefrac{1}{2}$ $\Diamond$-Norm \\ \hline
Gi & 0.04 & 0.04 & 0.04 \\ \hline
Gx & 0.038 & 0.039 & 0.039 \\ \hline
Gy & 0.038 & 0.038 & 0.039 \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|}
\hline
Gate & Error Generator \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
1\e{-5} & 3\e{-4} & 4\e{-4} & -3\e{-4} \\ 
-0.001 & -0.057 & -0.007 & -0.009 \\ 
5\e{-4} & 0.001 & -0.054 & 0.004 \\ 
0.002 & 2\e{-5} & 0.001 & -0.053
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
-1\e{-5} & 8\e{-5} & 0.011 & 5\e{-4} \\ 
4\e{-4} & -0.051 & 9\e{-5} & -0.008 \\ 
0.017 & 0.001 & -0.05 & -0.001 \\ 
4\e{-4} & 3\e{-6} & 0.001 & -0.054
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
2\e{-5} & -0.011 & -6\e{-5} & -0.001 \\ 
-0.016 & -0.052 & 0.001 & 6\e{-5} \\ 
3\e{-4} & 0.002 & -0.05 & 6\e{-5} \\ 
-3\e{-4} & -0.001 & -0.008 & -0.053
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{Comparison of GST estimated gates to target gates}.  This table presents, for each of the gates, three different measures of distance or discrepancy from the GST estimate to the ideal target operation.  See text for more detail.  The column labeled ``Error Generator'' gives the Lindbladian $\mathbb{L}$ that describes \emph{how} the gate is failing to match the target.  This error generator is defined by the equation $\hat{G} = G_{\mathrm{target}}e^{\mathbb{L}}$. \label{bestTargetGatesetVsTargetTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|}
\hline
Gate & Eigenvalues & Fixed pt & Rotn. axis & Diag. decay & Off-diag. decay \\ \hline
Gi & $ \begin{array}{c}
1 \\ 
0.949 \\ 
0.946e^{i0.00} \\ 
0.946e^{-i0.00}
 \end{array} $
 & $ \begin{array}{c}
0.998 \\ 
-0.022 \\ 
0.01 \\ 
0.032
 \end{array} $
 & $ \begin{array}{c}
0.007 \\ 
-0.934 \\ 
0.149 \\ 
0.325
 \end{array} $
 & 0.051 & 0.054 \\ \hline
Gx & $ \begin{array}{c}
0.949e^{i1.57} \\ 
0.949e^{-i1.57} \\ 
1 \\ 
0.951
 \end{array} $
 & $ \begin{array}{c}
1 \\ 
0.006 \\ 
1\e{-4} \\ 
0.011
 \end{array} $
 & $ \begin{array}{c}
0.002 \\ 
-1 \\ 
2\e{-7} \\ 
-0.001
 \end{array} $
 & 0.049 & 0.051 \\ \hline
Gy & $ \begin{array}{c}
0.949e^{i1.57} \\ 
0.949e^{-i1.57} \\ 
1 \\ 
0.951
 \end{array} $
 & $ \begin{array}{c}
1 \\ 
-0.001 \\ 
0.007 \\ 
0.01
 \end{array} $
 & $ \begin{array}{c}
-6\e{-6} \\ 
0.005 \\ 
-1 \\ 
0.001
 \end{array} $
 & 0.049 & 0.051 \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|c|c|c|}
\hline
\multirow{2}{*}{Gate} & \multirow{2}{*}{Angle} & \multicolumn{3}{c|}{Angle between Rotation Axes} \\ \cline{3-5}
 & & Gi & Gx & Gy \\ \hline
Gi & 0.001$\pi$ &  & 0.117$\pi$ & 0.549$\pi$ \\ \hline
Gx & 0.5$\pi$ & 0.117$\pi$ &  & 0.502$\pi$ \\ \hline
Gy & 0.5$\pi$ & 0.549$\pi$ & 0.502$\pi$ &  \\ \hline
\end{tabular}

\caption{\textbf{Eigen-decomposition of estimated gates}.  Each estimated gate is described in terms of: (1) the eigenvalues of the superoperator; (2) the gate's fixed point (as a vector in $\mathcal{B}(\mathcal{H})$, in the Pauli basis); (3)  the axis around which it rotates, as a vector in $\mathcal{B}(\mathcal{H})$; (4) the angle of the rotation that it applies; (5) the decay rate along the axis of rotation (``diagonal decay''); and (6) the decay rate perpendicular to the axis of rotation (``off-diagonal decay'').  ``--'' indicates that the decomposition failed or couldn't be interpreted. \label{bestTargetGatesetDecompTable}}
\end{center}
\end{table}


\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|}
\hline
Gate & \begin{tabular}{c}Process\\Infidelity\end{tabular} & \begin{tabular}{c}$\nicefrac{1}{2}$ Trace\\Distance\end{tabular} & \begin{tabular}{c}Rotation\\Axis\end{tabular} & \begin{tabular}{c}Rotation\\Angle\end{tabular} & Sanity \checkmark \\ \hline
Gi & 0.04 & 0.04 & $ \begin{array}{c}
0 \\ 
0.261 \\ 
0.741 \\ 
-0.618
 \end{array} $
 & 0.002$\pi$ & 5\e{-6} \\ \hline
Gx & 0.038 & 0.039 & $ \begin{array}{c}
0 \\ 
1 \\ 
-0.003 \\ 
3\e{-4}
 \end{array} $
 & 0.5$\pi$ & 0.001 \\ \hline
Gy & 0.038 & 0.038 & $ \begin{array}{c}
0 \\ 
-0.003 \\ 
1 \\ 
4\e{-4}
 \end{array} $
 & 0.5$\pi$ & 0.001 \\ \hline
\end{tabular}

\caption{Information pertaining to the closest unitary gate to each of the estimated gates.\label{bestTargetGatesetClosestUnitaryTable}}
\end{center}
\end{table}


\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
Gate & Choi matrix (Pauli basis) & Eigenvalues \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
0.96 & 0.001e^{-i1.75} & 0.002e^{-i1.48} & 0.002e^{i1.39} \\ 
0.001e^{i1.75} & 0.012 & 0.001e^{-i2.79} & 0.002e^{i3.13} \\ 
0.002e^{i1.48} & 0.001e^{i2.79} & 0.013 & 0.001e^{i0.26} \\ 
0.002e^{-i1.39} & 0.002e^{-i3.13} & 0.001e^{-i0.26} & 0.014
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.011 \\ 
0.013 \\ 
0.017 \\ 
0.96
 \end{array} $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
0.487 & 0.475e^{i1.57} & 0.005e^{-i0.32} & 0.002e^{i0.96} \\ 
0.475e^{-i1.57} & 0.488 & 0.004e^{-i1.83} & 0.001e^{i2.57} \\ 
0.005e^{i0.32} & 0.004e^{i1.83} & 0.013 & 7\e{-5}e^{-i1.77} \\ 
0.002e^{-i0.96} & 0.001e^{-i2.57} & 7\e{-5}e^{i1.77} & 0.012
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.01 \\ 
0.013 \\ 
0.015 \\ 
0.962
 \end{array} $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
0.488 & 0.005e^{-i2.74} & 0.474e^{i1.57} & 0.002e^{i1.17} \\ 
0.005e^{i2.74} & 0.012 & 0.004e^{-i1.75} & 1\e{-4}e^{i2.26} \\ 
0.474e^{-i1.57} & 0.004e^{i1.75} & 0.488 & 0.001e^{i2.44} \\ 
0.002e^{-i1.17} & 1\e{-4}e^{-i2.26} & 0.001e^{-i2.44} & 0.012
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.01 \\ 
0.012 \\ 
0.015 \\ 
0.962
 \end{array} $
 \\ \hline
\end{tabular}

\caption{\textbf{Choi matrix representation of the GST estimated gate set}.  This table lists Choi representations of the estimated gates, and their eigenvalues.  Unitary gates have a spectrum $(1,0,0\ldots)$, just like pure quantum states.  Negative eigenvalues are non-physical, and may represent either statistical fluctuations or violations of the CPTP model used by GST.\label{bestTargetGatesetChoiTable}}
\end{center}
\end{table}


\clearpage

\subsection{Gauge optimization to the target SPAM}
The gauge is chosen to minimize the Frobenius distance between the estimate and the target gates, with 99\% weight given to the SPAM parameters, 1\% to the gate parameters.

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
Operator & Matrix & Hilbert-Schmidt vector (Pauli basis) \\ \hline
$\rho_{0}$ & $ \left(\!\!\begin{array}{cc}
0.983 & 9\e{-6}e^{i1.05} \\ 
9\e{-6}e^{-i1.05} & 0.037
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.721 \\ 
6\e{-6} \\ 
-1\e{-5} \\ 
0.669
 \end{array} $
 \\ \hline
$E_{0}$ & $ \left(\!\!\begin{array}{cc}
0.055 & 2\e{-5}e^{i1.19} \\ 
2\e{-5}e^{-i1.19} & 1.002
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.748 \\ 
1\e{-5} \\ 
-3\e{-5} \\ 
-0.669
 \end{array} $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the SPAM operations}.  Compare to Table \ref{targetSpamTable}.\label{bestTargetSpamGatesetSpamTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
 & $E_{0}$ & $E_C$ \\ \hline
$\rho_{0}$ & 0.091 & 0.909 \\ \hline
\end{tabular}

\caption{\textbf{GST estimate of SPAM probabilities}.  Computed by taking the dot products of vectors in Table \ref{bestTargetSpamGatesetSpamTable}.\label{bestTargetSpamGatesetSpamParametersTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|}
\hline
Gate & Superoperator (Pauli basis) \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
1 & 4\e{-4} & 4\e{-4} & -4\e{-4} \\ 
-2\e{-4} & 0.945 & -0.006 & -0.009 \\ 
2\e{-4} & 0.001 & 0.947 & 0.004 \\ 
0.005 & 5\e{-5} & 0.001 & 0.948
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
1.001 & 1\e{-4} & 0.009 & -0.01 \\ 
5\e{-4} & 0.951 & -0.007 & -0.001 \\ 
0.067 & 3\e{-4} & 0.004 & -0.947 \\ 
0.072 & 0.002 & 0.952 & -0.007
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
1.001 & -0.009 & -5\e{-5} & -0.01 \\ 
-0.066 & 0.003 & -0.009 & 0.947 \\ 
2\e{-4} & 0.002 & 0.951 & 4\e{-4} \\ 
0.072 & -0.952 & -0.003 & -0.004
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the logic gate operations}.  Compare to Table \ref{targetGatesTable}.\label{bestTargetSpamGatesetGatesTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|}
\hline
Gate & \begin{tabular}{c}Process\\Infidelity\end{tabular} & \begin{tabular}{c}$\nicefrac{1}{2}$ Trace\\Distance\end{tabular} & $\nicefrac{1}{2}$ $\Diamond$-Norm \\ \hline
Gi & 0.04 & 0.04 & 0.04 \\ \hline
Gx & 0.038 & 0.06 & 0.072 \\ \hline
Gy & 0.037 & 0.059 & 0.073 \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|}
\hline
Gate & Error Generator \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
4\e{-5} & 4\e{-4} & 4\e{-4} & -4\e{-4} \\ 
-2\e{-4} & -0.057 & -0.007 & -0.009 \\ 
2\e{-4} & 0.001 & -0.054 & 0.004 \\ 
0.005 & 5\e{-5} & 0.001 & -0.053
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
3\e{-5} & 1\e{-4} & 0.015 & -4\e{-4} \\ 
0.001 & -0.051 & -0.005 & -0.007 \\ 
0.112 & 0.001 & -0.044 & 0.002 \\ 
0.005 & 0.001 & 0.005 & -0.061
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
1\e{-4} & -0.015 & -1\e{-4} & -0.001 \\ 
-0.111 & -0.047 & -0.004 & -0.003 \\ 
4\e{-4} & 0.002 & -0.05 & -0.001 \\ 
0.006 & -0.005 & -0.01 & -0.058
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{Comparison of GST estimated gates to target gates}.  This table presents, for each of the gates, three different measures of distance or discrepancy from the GST estimate to the ideal target operation.  See text for more detail.  The column labeled ``Error Generator'' gives the Lindbladian $\mathbb{L}$ that describes \emph{how} the gate is failing to match the target.  This error generator is defined by the equation $\hat{G} = G_{\mathrm{target}}e^{\mathbb{L}}$. \label{bestTargetSpamGatesetVsTargetTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|}
\hline
Gate & Eigenvalues & Fixed pt & Rotn. axis & Diag. decay & Off-diag. decay \\ \hline
Gi & $ \begin{array}{c}
1 \\ 
0.949 \\ 
0.946e^{i0.00} \\ 
0.946e^{-i0.00}
 \end{array} $
 & $ \begin{array}{c}
0.991 \\ 
-0.02 \\ 
0.009 \\ 
0.091
 \end{array} $
 & $ \begin{array}{c}
-0.008 \\ 
0.933 \\ 
-0.149 \\ 
-0.326
 \end{array} $
 & 0.051 & 0.054 \\ \hline
Gx & $ \begin{array}{c}
0.949e^{i1.57} \\ 
0.949e^{-i1.57} \\ 
1 \\ 
0.951
 \end{array} $
 & $ \begin{array}{c}
0.995 \\ 
0.008 \\ 
-4\e{-4} \\ 
0.071
 \end{array} $
 & $ \begin{array}{c}
0.002 \\ 
-1 \\ 
0.001 \\ 
-0.001
 \end{array} $
 & 0.049 & 0.051 \\ \hline
Gy & $ \begin{array}{c}
0.949e^{i1.57} \\ 
0.949e^{-i1.57} \\ 
1 \\ 
0.951
 \end{array} $
 & $ \begin{array}{c}
0.995 \\ 
0.001 \\ 
0.006 \\ 
0.071
 \end{array} $
 & $ \begin{array}{c}
-0.001 \\ 
0.006 \\ 
-1 \\ 
-0.003
 \end{array} $
 & 0.049 & 0.051 \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|c|c|c|}
\hline
\multirow{2}{*}{Gate} & \multirow{2}{*}{Angle} & \multicolumn{3}{c|}{Angle between Rotation Axes} \\ \cline{3-5}
 & & Gi & Gx & Gy \\ \hline
Gi & 0.001$\pi$ &  & 0.883$\pi$ & 0.45$\pi$ \\ \hline
Gx & 0.5$\pi$ & 0.883$\pi$ &  & 0.502$\pi$ \\ \hline
Gy & 0.5$\pi$ & 0.45$\pi$ & 0.502$\pi$ &  \\ \hline
\end{tabular}

\caption{\textbf{Eigen-decomposition of estimated gates}.  Each estimated gate is described in terms of: (1) the eigenvalues of the superoperator; (2) the gate's fixed point (as a vector in $\mathcal{B}(\mathcal{H})$, in the Pauli basis); (3)  the axis around which it rotates, as a vector in $\mathcal{B}(\mathcal{H})$; (4) the angle of the rotation that it applies; (5) the decay rate along the axis of rotation (``diagonal decay''); and (6) the decay rate perpendicular to the axis of rotation (``off-diagonal decay'').  ``--'' indicates that the decomposition failed or couldn't be interpreted. \label{bestTargetSpamGatesetDecompTable}}
\end{center}
\end{table}


\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|}
\hline
Gate & \begin{tabular}{c}Process\\Infidelity\end{tabular} & \begin{tabular}{c}$\nicefrac{1}{2}$ Trace\\Distance\end{tabular} & \begin{tabular}{c}Rotation\\Axis\end{tabular} & \begin{tabular}{c}Rotation\\Angle\end{tabular} & Sanity \checkmark \\ \hline
Gi & 0.04 & 0.04 & $ \begin{array}{c}
0 \\ 
0.26 \\ 
0.743 \\ 
-0.616
 \end{array} $
 & 0.002$\pi$ & 3\e{-5} \\ \hline
Gx & 0.038 & 0.06 & $ \begin{array}{c}
0 \\ 
-1 \\ 
0.003 \\ 
-0.002
 \end{array} $
 & 0.5$\pi$ & 0.023 \\ \hline
Gy & 0.037 & 0.059 & $ \begin{array}{c}
0 \\ 
-0.003 \\ 
1 \\ 
0.002
 \end{array} $
 & 0.5$\pi$ & 0.022 \\ \hline
\end{tabular}

\caption{Information pertaining to the closest unitary gate to each of the estimated gates.\label{bestTargetSpamGatesetClosestUnitaryTable}}
\end{center}
\end{table}


\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
Gate & Choi matrix (Pauli basis) & Eigenvalues \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
0.96 & 0.001e^{-i1.53} & 0.002e^{-i1.51} & 0.002e^{i1.04} \\ 
0.001e^{i1.53} & 0.012 & 0.002e^{-i2.36} & 0.002e^{-i3.12} \\ 
0.002e^{i1.51} & 0.002e^{i2.36} & 0.013 & 0.001e^{i0.13} \\ 
0.002e^{-i1.04} & 0.002e^{i3.12} & 0.001e^{-i0.13} & 0.014
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.01 \\ 
0.013 \\ 
0.017 \\ 
0.96
 \end{array} $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
0.487 & 0.475e^{i1.57} & 0.019e^{-i0.04} & 0.016e^{i0.12} \\ 
0.475e^{-i1.57} & 0.489 & 0.021e^{-i1.66} & 0.014e^{i1.56} \\ 
0.019e^{i0.04} & 0.021e^{i1.66} & 0.015 & 0.001e^{-i0.08} \\ 
0.016e^{-i0.12} & 0.014e^{-i1.56} & 0.001e^{i0.08} & 0.01
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
-0.01 \\ 
0.014 \\ 
0.033 \\ 
0.963
 \end{array} $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
0.488 & 0.019e^{-i3.09} & 0.475e^{i1.57} & 0.016e^{i0.16} \\ 
0.019e^{i3.09} & 0.014 & 0.021e^{-i1.66} & 0.001e^{i3.08} \\ 
0.475e^{-i1.57} & 0.021e^{i1.66} & 0.488 & 0.014e^{i1.62} \\ 
0.016e^{-i0.16} & 0.001e^{-i3.08} & 0.014e^{-i1.62} & 0.011
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
-0.009 \\ 
0.013 \\ 
0.033 \\ 
0.963
 \end{array} $
 \\ \hline
\end{tabular}

\caption{\textbf{Choi matrix representation of the GST estimated gate set}.  This table lists Choi representations of the estimated gates, and their eigenvalues.  Unitary gates have a spectrum $(1,0,0\ldots)$, just like pure quantum states.  Negative eigenvalues are non-physical, and may represent either statistical fluctuations or violations of the CPTP model used by GST.\label{bestTargetSpamGatesetChoiTable}}
\end{center}
\end{table}

\clearpage

\subsection{Gauge optimization to the target gates:}
The gauge is chosen to minimize the Frobenius distance between the estimate and the target gates, with 99\% weight given to the gate parameters, 1\% to the SPAM parameters.


\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
Operator & Matrix & Hilbert-Schmidt vector (Pauli basis) \\ \hline
$\rho_{0}$ & $ \left(\!\!\begin{array}{cc}
0.99 & 0.001e^{-i0.15} \\ 
0.001e^{i0.15} & 0.01
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.707 \\ 
0.002 \\ 
3\e{-4} \\ 
0.693
 \end{array} $
 \\ \hline
$E_{0}$ & $ \left(\!\!\begin{array}{cc}
0.083 & 0.002e^{i1.84} \\ 
0.002e^{-i1.84} & 0.918
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.708 \\ 
-0.001 \\ 
-0.002 \\ 
-0.591
 \end{array} $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the SPAM operations}.  Compare to Table \ref{targetSpamTable}.\label{bestTargetGatesGatesetSpamTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
 & $E_{0}$ & $E_C$ \\ \hline
$\rho_{0}$ & 0.091 & 0.909 \\ \hline
\end{tabular}

\caption{\textbf{GST estimate of SPAM probabilities}.  Computed by taking the dot products of vectors in Table \ref{bestTargetGatesGatesetSpamTable}.\label{bestTargetGatesGatesetSpamParametersTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|}
\hline
Gate & Superoperator (Pauli basis) \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
1 & 3\e{-4} & 3\e{-4} & 8\e{-5} \\ 
-0.001 & 0.945 & -0.006 & -0.009 \\ 
0.001 & 0.001 & 0.947 & 0.004 \\ 
0.001 & 2\e{-5} & 0.001 & 0.949
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
1 & 7\e{-5} & -5\e{-5} & 2\e{-5} \\ 
4\e{-4} & 0.951 & -0.005 & -0.005 \\ 
2\e{-4} & 3\e{-4} & 1\e{-4} & -0.949 \\ 
2\e{-4} & 3\e{-4} & 0.949 & -0.002
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
1 & -2\e{-5} & -3\e{-5} & 5\e{-5} \\ 
2\e{-4} & 6\e{-5} & -0.004 & 0.949 \\ 
3\e{-4} & 0.001 & 0.951 & 0.001 \\ 
-2\e{-4} & -0.949 & -0.006 & -0.001
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the logic gate operations}.  Compare to Table \ref{targetGatesTable}.\label{bestTargetGatesGatesetGatesTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|}
\hline
Gate & \begin{tabular}{c}Process\\Infidelity\end{tabular} & \begin{tabular}{c}$\nicefrac{1}{2}$ Trace\\Distance\end{tabular} & $\nicefrac{1}{2}$ $\Diamond$-Norm \\ \hline
Gi & 0.04 & 0.04 & 0.04 \\ \hline
Gx & 0.038 & 0.038 & 0.038 \\ \hline
Gy & 0.038 & 0.038 & 0.038 \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|}
\hline
Gate & Error Generator \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
5\e{-7} & 3\e{-4} & 3\e{-4} & 8\e{-5} \\ 
-0.001 & -0.057 & -0.007 & -0.009 \\ 
0.001 & 0.001 & -0.054 & 0.004 \\ 
0.001 & 2\e{-5} & 0.001 & -0.053
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
-8\e{-6} & 7\e{-5} & -5\e{-5} & -2\e{-5} \\ 
4\e{-4} & -0.051 & 5\e{-4} & -0.008 \\ 
3\e{-4} & 5\e{-4} & -0.05 & -0.001 \\ 
3\e{-5} & -1\e{-6} & 0.001 & -0.054
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
2\e{-5} & 2\e{-5} & -3\e{-5} & 6\e{-5} \\ 
4\e{-4} & -0.052 & 0.002 & 1\e{-4} \\ 
3\e{-4} & 0.002 & -0.05 & 6\e{-5} \\ 
1\e{-5} & -4\e{-4} & -0.008 & -0.053
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{Comparison of GST estimated gates to target gates}.  This table presents, for each of the gates, three different measures of distance or discrepancy from the GST estimate to the ideal target operation.  See text for more detail.  The column labeled ``Error Generator'' gives the Lindbladian $\mathbb{L}$ that describes \emph{how} the gate is failing to match the target.  This error generator is defined by the equation $\hat{G} = G_{\mathrm{target}}e^{\mathbb{L}}$. \label{bestTargetGatesGatesetVsTargetTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|}
\hline
Gate & Eigenvalues & Fixed pt & Rotn. axis & Diag. decay & Off-diag. decay \\ \hline
Gi & $ \begin{array}{c}
1 \\ 
0.949 \\ 
0.946e^{i0.00} \\ 
0.946e^{-i0.00}
 \end{array} $
 & $ \begin{array}{c}
0.999 \\ 
-0.024 \\ 
0.011 \\ 
0.023
 \end{array} $
 & $ \begin{array}{c}
0.004 \\ 
-0.934 \\ 
0.149 \\ 
0.325
 \end{array} $
 & 0.051 & 0.054 \\ \hline
Gx & $ \begin{array}{c}
0.949e^{i1.57} \\ 
0.949e^{-i1.57} \\ 
1 \\ 
0.951
 \end{array} $
 & $ \begin{array}{c}
1 \\ 
0.008 \\ 
-1\e{-5} \\ 
2\e{-4}
 \end{array} $
 & $ \begin{array}{c}
0.001 \\ 
-1 \\ 
-1\e{-6} \\ 
-3\e{-4}
 \end{array} $
 & 0.049 & 0.051 \\ \hline
Gy & $ \begin{array}{c}
1 \\ 
0.951 \\ 
0.949e^{i1.57} \\ 
0.949e^{-i1.57}
 \end{array} $
 & $ \begin{array}{c}
1 \\ 
-2\e{-5} \\ 
0.007 \\ 
-3\e{-4}
 \end{array} $
 & $ \begin{array}{c}
-0.001 \\ 
0.005 \\ 
-1 \\ 
0.001
 \end{array} $
 & 0.049 & 0.051 \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|c|c|c|}
\hline
\multirow{2}{*}{Gate} & \multirow{2}{*}{Angle} & \multicolumn{3}{c|}{Angle between Rotation Axes} \\ \cline{3-5}
 & & Gi & Gx & Gy \\ \hline
Gi & 0.001$\pi$ &  & 0.117$\pi$ & 0.549$\pi$ \\ \hline
Gx & 0.5$\pi$ & 0.117$\pi$ &  & 0.502$\pi$ \\ \hline
Gy & 0.5$\pi$ & 0.549$\pi$ & 0.502$\pi$ &  \\ \hline
\end{tabular}

\caption{\textbf{Eigen-decomposition of estimated gates}.  Each estimated gate is described in terms of: (1) the eigenvalues of the superoperator; (2) the gate's fixed point (as a vector in $\mathcal{B}(\mathcal{H})$, in the Pauli basis); (3)  the axis around which it rotates, as a vector in $\mathcal{B}(\mathcal{H})$; (4) the angle of the rotation that it applies; (5) the decay rate along the axis of rotation (``diagonal decay''); and (6) the decay rate perpendicular to the axis of rotation (``off-diagonal decay'').  ``--'' indicates that the decomposition failed or couldn't be interpreted. \label{bestTargetGatesGatesetDecompTable}}
\end{center}
\end{table}


\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|}
\hline
Gate & \begin{tabular}{c}Process\\Infidelity\end{tabular} & \begin{tabular}{c}$\nicefrac{1}{2}$ Trace\\Distance\end{tabular} & \begin{tabular}{c}Rotation\\Axis\end{tabular} & \begin{tabular}{c}Rotation\\Angle\end{tabular} & Sanity \checkmark \\ \hline
Gi & 0.04 & 0.04 & $ \begin{array}{c}
0 \\ 
-0.261 \\ 
-0.741 \\ 
0.618
 \end{array} $
 & 0.002$\pi$ & 5\e{-6} \\ \hline
Gx & 0.038 & 0.038 & $ \begin{array}{c}
0 \\ 
-1 \\ 
0.003 \\ 
-3\e{-6}
 \end{array} $
 & 0.5$\pi$ & 5\e{-7} \\ \hline
Gy & 0.038 & 0.038 & $ \begin{array}{c}
0 \\ 
0.003 \\ 
-1 \\ 
7\e{-8}
 \end{array} $
 & 0.5$\pi$ & 4\e{-7} \\ \hline
\end{tabular}

\caption{Information pertaining to the closest unitary gate to each of the estimated gates.\label{bestTargetGatesGatesetClosestUnitaryTable}}
\end{center}
\end{table}


\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
Gate & Choi matrix (Pauli basis) & Eigenvalues \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
0.96 & 0.001e^{-i1.80} & 0.002e^{-i1.48} & 0.002e^{i1.40} \\ 
0.001e^{i1.80} & 0.012 & 0.001e^{-i2.93} & 0.002e^{i3.12} \\ 
0.002e^{i1.48} & 0.001e^{i2.93} & 0.013 & 0.001e^{i0.28} \\ 
0.002e^{-i1.40} & 0.002e^{-i3.12} & 0.001e^{-i0.28} & 0.014
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.011 \\ 
0.013 \\ 
0.017 \\ 
0.96
 \end{array} $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
0.487 & 0.475e^{i1.57} & 0.001e^{-i1.55} & 0.001e^{i1.52} \\ 
0.475e^{-i1.57} & 0.488 & 0.001e^{-i3.10} & 0.001e^{i3.10} \\ 
0.001e^{i1.55} & 0.001e^{i3.10} & 0.013 & 9\e{-5}e^{-i2.02} \\ 
0.001e^{-i1.52} & 0.001e^{-i3.10} & 9\e{-5}e^{i2.02} & 0.012
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.011 \\ 
0.013 \\ 
0.014 \\ 
0.962
 \end{array} $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
0.488 & 0.002e^{-i1.54} & 0.474e^{i1.57} & 0.001e^{i1.61} \\ 
0.002e^{i1.54} & 0.012 & 0.001e^{i3.05} & 1\e{-4}e^{i2.00} \\ 
0.474e^{-i1.57} & 0.001e^{-i3.05} & 0.488 & 0.001e^{-i3.09} \\ 
0.001e^{-i1.61} & 1\e{-4}e^{-i2.00} & 0.001e^{i3.09} & 0.012
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.011 \\ 
0.012 \\ 
0.015 \\ 
0.962
 \end{array} $
 \\ \hline
\end{tabular}

\caption{\textbf{Choi matrix representation of the GST estimated gate set}.  This table lists Choi representations of the estimated gates, and their eigenvalues.  Unitary gates have a spectrum $(1,0,0\ldots)$, just like pure quantum states.  Negative eigenvalues are non-physical, and may represent either statistical fluctuations or violations of the CPTP model used by GST.\label{bestTargetGatesGatesetChoiTable}}
\end{center}
\end{table}

\clearpage

\subsection{Gauge optimization to TP}
The gauge is chosen to make the gate set as trace-preserving as possible.  This is done by minimizing the sum of the squared euclidian distance between the first row of each estimated gate matrix and the vector $(1,0,\cdots 0)$.  The Frobenius distance between the estimate and the target gate set, weighted by $10^{-4}$, is added to the aforementioned sum of distances is to give the final objective function.  Ideally, a perfectly TP gate set will be found and the Frobenius distance term causes the optimization to choose the gate set closest to the target gates that is also in TP.  If a perfectly TP gate set cannot be gauge-optimized to, then the resulting gate set compromises between being TP and being close to the target gate set, with the intent that the TP penalty term dominates.

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
Operator & Matrix & Hilbert-Schmidt vector (Pauli basis) \\ \hline
$\rho_{0}$ & $ \left(\!\!\begin{array}{cc}
0.99 & 0.001e^{-i0.15} \\ 
0.001e^{i0.15} & 0.01
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.707 \\ 
0.002 \\ 
3\e{-4} \\ 
0.693
 \end{array} $
 \\ \hline
$E_{0}$ & $ \left(\!\!\begin{array}{cc}
0.083 & 0.002e^{i1.84} \\ 
0.002e^{-i1.84} & 0.918
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.708 \\ 
-0.001 \\ 
-0.002 \\ 
-0.591
 \end{array} $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the SPAM operations}.  Compare to Table \ref{targetSpamTable}.\label{bestTPGatesetSpamTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
 & $E_{0}$ & $E_C$ \\ \hline
$\rho_{0}$ & 0.091 & 0.909 \\ \hline
\end{tabular}

\caption{\textbf{GST estimate of SPAM probabilities}.  Computed by taking the dot products of vectors in Table \ref{bestTPGatesetSpamTable}.\label{bestTPGatesetSpamParametersTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|}
\hline
Gate & Superoperator (Pauli basis) \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
1 & 3\e{-4} & 3\e{-4} & 8\e{-5} \\ 
-0.001 & 0.945 & -0.006 & -0.009 \\ 
0.001 & 0.001 & 0.947 & 0.004 \\ 
0.001 & 2\e{-5} & 0.001 & 0.949
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
1 & 7\e{-5} & -3\e{-5} & 6\e{-6} \\ 
4\e{-4} & 0.951 & -0.005 & -0.005 \\ 
2\e{-4} & 3\e{-4} & 1\e{-4} & -0.949 \\ 
2\e{-4} & 3\e{-4} & 0.949 & -0.002
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
1 & -1\e{-6} & -2\e{-5} & 0 \\ 
3\e{-4} & 7\e{-5} & -0.004 & 0.949 \\ 
3\e{-4} & 0.001 & 0.951 & 0.001 \\ 
-2\e{-4} & -0.949 & -0.006 & -0.001
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the logic gate operations}.  Compare to Table \ref{targetGatesTable}.\label{bestTPGatesetGatesTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|}
\hline
Gate & \begin{tabular}{c}Process\\Infidelity\end{tabular} & \begin{tabular}{c}$\nicefrac{1}{2}$ Trace\\Distance\end{tabular} & $\nicefrac{1}{2}$ $\Diamond$-Norm \\ \hline
Gi & 0.04 & 0.04 & 0.04 \\ \hline
Gx & 0.038 & 0.038 & 0.038 \\ \hline
Gy & 0.038 & 0.038 & 0.038 \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|}
\hline
Gate & Error Generator \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
5\e{-7} & 3\e{-4} & 3\e{-4} & 8\e{-5} \\ 
-0.001 & -0.057 & -0.007 & -0.009 \\ 
0.001 & 0.001 & -0.054 & 0.004 \\ 
0.001 & 2\e{-5} & 0.001 & -0.053
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
-8\e{-6} & 7\e{-5} & -3\e{-5} & -2\e{-5} \\ 
4\e{-4} & -0.051 & 5\e{-4} & -0.008 \\ 
3\e{-4} & 5\e{-4} & -0.05 & -0.001 \\ 
3\e{-5} & 3\e{-6} & 0.001 & -0.054
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
2\e{-5} & -9\e{-7} & -3\e{-5} & 9\e{-7} \\ 
4\e{-4} & -0.052 & 0.002 & 8\e{-5} \\ 
3\e{-4} & 0.002 & -0.05 & 6\e{-5} \\ 
4\e{-5} & -4\e{-4} & -0.008 & -0.053
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{Comparison of GST estimated gates to target gates}.  This table presents, for each of the gates, three different measures of distance or discrepancy from the GST estimate to the ideal target operation.  See text for more detail.  The column labeled ``Error Generator'' gives the Lindbladian $\mathbb{L}$ that describes \emph{how} the gate is failing to match the target.  This error generator is defined by the equation $\hat{G} = G_{\mathrm{target}}e^{\mathbb{L}}$. \label{bestTPGatesetVsTargetTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|}
\hline
Gate & Eigenvalues & Fixed pt & Rotn. axis & Diag. decay & Off-diag. decay \\ \hline
Gi & $ \begin{array}{c}
1 \\ 
0.949 \\ 
0.946e^{i0.00} \\ 
0.946e^{-i0.00}
 \end{array} $
 & $ \begin{array}{c}
0.999 \\ 
-0.024 \\ 
0.011 \\ 
0.023
 \end{array} $
 & $ \begin{array}{c}
0.004 \\ 
-0.934 \\ 
0.149 \\ 
0.325
 \end{array} $
 & 0.051 & 0.054 \\ \hline
Gx & $ \begin{array}{c}
0.949e^{i1.57} \\ 
0.949e^{-i1.57} \\ 
1 \\ 
0.951
 \end{array} $
 & $ \begin{array}{c}
1 \\ 
0.008 \\ 
-1\e{-5} \\ 
2\e{-4}
 \end{array} $
 & $ \begin{array}{c}
0.001 \\ 
-1 \\ 
1\e{-6} \\ 
-3\e{-4}
 \end{array} $
 & 0.049 & 0.051 \\ \hline
Gy & $ \begin{array}{c}
1 \\ 
0.951 \\ 
0.949e^{i1.57} \\ 
0.949e^{-i1.57}
 \end{array} $
 & $ \begin{array}{c}
1 \\ 
-2\e{-6} \\ 
0.007 \\ 
-2\e{-4}
 \end{array} $
 & $ \begin{array}{c}
-0.001 \\ 
0.005 \\ 
-1 \\ 
0.001
 \end{array} $
 & 0.049 & 0.051 \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|c|c|c|}
\hline
\multirow{2}{*}{Gate} & \multirow{2}{*}{Angle} & \multicolumn{3}{c|}{Angle between Rotation Axes} \\ \cline{3-5}
 & & Gi & Gx & Gy \\ \hline
Gi & 0.001$\pi$ &  & 0.117$\pi$ & 0.549$\pi$ \\ \hline
Gx & 0.5$\pi$ & 0.117$\pi$ &  & 0.502$\pi$ \\ \hline
Gy & 0.5$\pi$ & 0.549$\pi$ & 0.502$\pi$ &  \\ \hline
\end{tabular}

\caption{\textbf{Eigen-decomposition of estimated gates}.  Each estimated gate is described in terms of: (1) the eigenvalues of the superoperator; (2) the gate's fixed point (as a vector in $\mathcal{B}(\mathcal{H})$, in the Pauli basis); (3)  the axis around which it rotates, as a vector in $\mathcal{B}(\mathcal{H})$; (4) the angle of the rotation that it applies; (5) the decay rate along the axis of rotation (``diagonal decay''); and (6) the decay rate perpendicular to the axis of rotation (``off-diagonal decay'').  ``--'' indicates that the decomposition failed or couldn't be interpreted. \label{bestTPGatesetDecompTable}}
\end{center}
\end{table}


\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|}
\hline
Gate & \begin{tabular}{c}Process\\Infidelity\end{tabular} & \begin{tabular}{c}$\nicefrac{1}{2}$ Trace\\Distance\end{tabular} & \begin{tabular}{c}Rotation\\Axis\end{tabular} & \begin{tabular}{c}Rotation\\Angle\end{tabular} & Sanity \checkmark \\ \hline
Gi & 0.04 & 0.04 & $ \begin{array}{c}
0 \\ 
0.261 \\ 
0.741 \\ 
-0.618
 \end{array} $
 & 0.002$\pi$ & 5\e{-6} \\ \hline
Gx & 0.038 & 0.038 & $ \begin{array}{c}
0 \\ 
1 \\ 
-0.003 \\ 
-2\e{-6}
 \end{array} $
 & 0.5$\pi$ & 5\e{-7} \\ \hline
Gy & 0.038 & 0.038 & $ \begin{array}{c}
0 \\ 
-0.003 \\ 
1 \\ 
3\e{-7}
 \end{array} $
 & 0.5$\pi$ & 4\e{-7} \\ \hline
\end{tabular}

\caption{Information pertaining to the closest unitary gate to each of the estimated gates.\label{bestTPGatesetClosestUnitaryTable}}
\end{center}
\end{table}


\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
Gate & Choi matrix (Pauli basis) & Eigenvalues \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
0.96 & 0.001e^{-i1.80} & 0.002e^{-i1.48} & 0.002e^{i1.40} \\ 
0.001e^{i1.80} & 0.012 & 0.001e^{-i2.93} & 0.002e^{i3.12} \\ 
0.002e^{i1.48} & 0.001e^{i2.93} & 0.013 & 0.001e^{i0.28} \\ 
0.002e^{-i1.40} & 0.002e^{-i3.12} & 0.001e^{-i0.28} & 0.014
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.011 \\ 
0.013 \\ 
0.017 \\ 
0.96
 \end{array} $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
0.487 & 0.475e^{i1.57} & 0.001e^{-i1.54} & 0.001e^{i1.52} \\ 
0.475e^{-i1.57} & 0.488 & 0.001e^{-i3.09} & 0.001e^{i3.10} \\ 
0.001e^{i1.54} & 0.001e^{i3.09} & 0.013 & 8\e{-5}e^{-i1.97} \\ 
0.001e^{-i1.52} & 0.001e^{-i3.10} & 8\e{-5}e^{i1.97} & 0.012
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.011 \\ 
0.013 \\ 
0.014 \\ 
0.962
 \end{array} $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
0.488 & 0.002e^{-i1.54} & 0.474e^{i1.57} & 0.001e^{i1.61} \\ 
0.002e^{i1.54} & 0.012 & 0.001e^{i3.07} & 1\e{-4}e^{i2.04} \\ 
0.474e^{-i1.57} & 0.001e^{-i3.07} & 0.488 & 0.001e^{-i3.09} \\ 
0.001e^{-i1.61} & 1\e{-4}e^{-i2.04} & 0.001e^{i3.09} & 0.012
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.011 \\ 
0.012 \\ 
0.015 \\ 
0.962
 \end{array} $
 \\ \hline
\end{tabular}

\caption{\textbf{Choi matrix representation of the GST estimated gate set}.  This table lists Choi representations of the estimated gates, and their eigenvalues.  Unitary gates have a spectrum $(1,0,0\ldots)$, just like pure quantum states.  Negative eigenvalues are non-physical, and may represent either statistical fluctuations or violations of the CPTP model used by GST.\label{bestTPGatesetChoiTable}}
\end{center}
\end{table}

\clearpage

\subsection{Gauge optimization to CPTP}
The gauge is chosen to minimize a quantity indicating the ``non-CPTP-ness'' of the gate set, along with a slight preference for being close to the target gate set.  In particular, we first optimize the gate set to TP (as described above), and then minimize the logarithm of a CP ``penalty term'' which we take as the sum of
\begin{enumerate}
\item the absolute values of any negative eigenvalues possessed by a gate's Choi matrix
\item the absolute value of any negative eigenvalues of a state preparation density matrix
\item the absolute difference between the trace of a state preparation density matrix and one
\item the distance between an eigenvalue of an effect operator and the interval $[0,1]$.
\end{enumerate}
within the space of TP-preserving gauge transformations.  This CP penalty term is added to the Frobenius distance between the estimate and the target gate set, weighted by $10^{-4}$ to give the final objective function.

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
Operator & Matrix & Hilbert-Schmidt vector (Pauli basis) \\ \hline
$\rho_{0}$ & $ \left(\!\!\begin{array}{cc}
0.99 & 0.001e^{-i0.16} \\ 
0.001e^{i0.16} & 0.01
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.707 \\ 
0.002 \\ 
3\e{-4} \\ 
0.693
 \end{array} $
 \\ \hline
$E_{0}$ & $ \left(\!\!\begin{array}{cc}
0.083 & 0.002e^{i1.84} \\ 
0.002e^{-i1.84} & 0.918
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.708 \\ 
-0.001 \\ 
-0.002 \\ 
-0.591
 \end{array} $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the SPAM operations}.  Compare to Table \ref{targetSpamTable}.\label{bestCPTPGatesetSpamTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
 & $E_{0}$ & $E_C$ \\ \hline
$\rho_{0}$ & 0.091 & 0.909 \\ \hline
\end{tabular}

\caption{\textbf{GST estimate of SPAM probabilities}.  Computed by taking the dot products of vectors in Table \ref{bestCPTPGatesetSpamTable}.\label{bestCPTPGatesetSpamParametersTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|}
\hline
Gate & Superoperator (Pauli basis) \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
1 & 3\e{-4} & 3\e{-4} & 8\e{-5} \\ 
-0.001 & 0.945 & -0.006 & -0.009 \\ 
0.001 & 0.001 & 0.947 & 0.004 \\ 
0.001 & 2\e{-5} & 0.001 & 0.949
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
1 & 7\e{-5} & 2\e{-9} & 2\e{-6} \\ 
4\e{-4} & 0.951 & -0.005 & -0.005 \\ 
2\e{-4} & 3\e{-4} & 4\e{-5} & -0.949 \\ 
2\e{-4} & 3\e{-4} & 0.949 & -0.002
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
1 & -3\e{-5} & -2\e{-5} & -2\e{-6} \\ 
2\e{-4} & 5\e{-5} & -0.004 & 0.949 \\ 
3\e{-4} & 0.001 & 0.951 & 0.001 \\ 
-2\e{-4} & -0.949 & -0.006 & -5\e{-4}
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the logic gate operations}.  Compare to Table \ref{targetGatesTable}.\label{bestCPTPGatesetGatesTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|}
\hline
Gate & \begin{tabular}{c}Process\\Infidelity\end{tabular} & \begin{tabular}{c}$\nicefrac{1}{2}$ Trace\\Distance\end{tabular} & $\nicefrac{1}{2}$ $\Diamond$-Norm \\ \hline
Gi & 0.04 & 0.04 & 0.04 \\ \hline
Gx & 0.038 & 0.038 & 0.038 \\ \hline
Gy & 0.038 & 0.038 & 0.038 \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|}
\hline
Gate & Error Generator \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
5\e{-7} & 3\e{-4} & 3\e{-4} & 8\e{-5} \\ 
-0.001 & -0.057 & -0.007 & -0.009 \\ 
0.001 & 0.001 & -0.054 & 0.004 \\ 
0.001 & 2\e{-5} & 0.001 & -0.053
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
-8\e{-6} & 7\e{-5} & -2\e{-6} & 2\e{-6} \\ 
4\e{-4} & -0.051 & 5\e{-4} & -0.009 \\ 
4\e{-4} & 5\e{-4} & -0.051 & -0.001 \\ 
5\e{-6} & -4\e{-5} & 0.001 & -0.054
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
2\e{-5} & -3\e{-5} & -3\e{-5} & 2\e{-5} \\ 
3\e{-4} & -0.052 & 0.002 & -5\e{-5} \\ 
3\e{-4} & 0.002 & -0.05 & 1\e{-4} \\ 
-3\e{-5} & -0.001 & -0.008 & -0.053
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{Comparison of GST estimated gates to target gates}.  This table presents, for each of the gates, three different measures of distance or discrepancy from the GST estimate to the ideal target operation.  See text for more detail.  The column labeled ``Error Generator'' gives the Lindbladian $\mathbb{L}$ that describes \emph{how} the gate is failing to match the target.  This error generator is defined by the equation $\hat{G} = G_{\mathrm{target}}e^{\mathbb{L}}$. \label{bestCPTPGatesetVsTargetTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|}
\hline
Gate & Eigenvalues & Fixed pt & Rotn. axis & Diag. decay & Off-diag. decay \\ \hline
Gi & $ \begin{array}{c}
1 \\ 
0.949 \\ 
0.946e^{i0.00} \\ 
0.946e^{-i0.00}
 \end{array} $
 & $ \begin{array}{c}
0.999 \\ 
-0.024 \\ 
0.011 \\ 
0.023
 \end{array} $
 & $ \begin{array}{c}
0.004 \\ 
-0.934 \\ 
0.149 \\ 
0.325
 \end{array} $
 & 0.051 & 0.054 \\ \hline
Gx & $ \begin{array}{c}
0.949e^{i1.57} \\ 
0.949e^{-i1.57} \\ 
1 \\ 
0.951
 \end{array} $
 & $ \begin{array}{c}
1 \\ 
0.008 \\ 
5\e{-6} \\ 
2\e{-4}
 \end{array} $
 & $ \begin{array}{c}
0.001 \\ 
-1 \\ 
-3\e{-5} \\ 
-3\e{-4}
 \end{array} $
 & 0.049 & 0.051 \\ \hline
Gy & $ \begin{array}{c}
1 \\ 
0.951 \\ 
0.949e^{i1.57} \\ 
0.949e^{-i1.57}
 \end{array} $
 & $ \begin{array}{c}
1 \\ 
-5\e{-5} \\ 
0.007 \\ 
-2\e{-4}
 \end{array} $
 & $ \begin{array}{c}
-5\e{-4} \\ 
0.005 \\ 
-1 \\ 
0.001
 \end{array} $
 & 0.049 & 0.051 \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|c|c|c|}
\hline
\multirow{2}{*}{Gate} & \multirow{2}{*}{Angle} & \multicolumn{3}{c|}{Angle between Rotation Axes} \\ \cline{3-5}
 & & Gi & Gx & Gy \\ \hline
Gi & 0.001$\pi$ &  & 0.117$\pi$ & 0.549$\pi$ \\ \hline
Gx & 0.5$\pi$ & 0.117$\pi$ &  & 0.502$\pi$ \\ \hline
Gy & 0.5$\pi$ & 0.549$\pi$ & 0.502$\pi$ &  \\ \hline
\end{tabular}

\caption{\textbf{Eigen-decomposition of estimated gates}.  Each estimated gate is described in terms of: (1) the eigenvalues of the superoperator; (2) the gate's fixed point (as a vector in $\mathcal{B}(\mathcal{H})$, in the Pauli basis); (3)  the axis around which it rotates, as a vector in $\mathcal{B}(\mathcal{H})$; (4) the angle of the rotation that it applies; (5) the decay rate along the axis of rotation (``diagonal decay''); and (6) the decay rate perpendicular to the axis of rotation (``off-diagonal decay'').  ``--'' indicates that the decomposition failed or couldn't be interpreted. \label{bestCPTPGatesetDecompTable}}
\end{center}
\end{table}


\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|}
\hline
Gate & \begin{tabular}{c}Process\\Infidelity\end{tabular} & \begin{tabular}{c}$\nicefrac{1}{2}$ Trace\\Distance\end{tabular} & \begin{tabular}{c}Rotation\\Axis\end{tabular} & \begin{tabular}{c}Rotation\\Angle\end{tabular} & Sanity \checkmark \\ \hline
Gi & 0.04 & 0.04 & $ \begin{array}{c}
0 \\ 
0.261 \\ 
0.741 \\ 
-0.619
 \end{array} $
 & 0.002$\pi$ & 5\e{-6} \\ \hline
Gx & 0.038 & 0.038 & $ \begin{array}{c}
0 \\ 
-1 \\ 
0.003 \\ 
-5\e{-6}
 \end{array} $
 & 0.5$\pi$ & 5\e{-7} \\ \hline
Gy & 0.038 & 0.038 & $ \begin{array}{c}
0 \\ 
-0.003 \\ 
1 \\ 
-3\e{-5}
 \end{array} $
 & 0.5$\pi$ & 3\e{-7} \\ \hline
\end{tabular}

\caption{Information pertaining to the closest unitary gate to each of the estimated gates.\label{bestCPTPGatesetClosestUnitaryTable}}
\end{center}
\end{table}


\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
Gate & Choi matrix (Pauli basis) & Eigenvalues \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
0.96 & 0.001e^{-i1.80} & 0.002e^{-i1.48} & 0.002e^{i1.40} \\ 
0.001e^{i1.80} & 0.012 & 0.001e^{-i2.93} & 0.002e^{i3.12} \\ 
0.002e^{i1.48} & 0.001e^{i2.93} & 0.013 & 0.001e^{i0.28} \\ 
0.002e^{-i1.40} & 0.002e^{-i3.12} & 0.001e^{-i0.28} & 0.014
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.011 \\ 
0.013 \\ 
0.017 \\ 
0.96
 \end{array} $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
0.487 & 0.475e^{i1.57} & 0.001e^{-i1.53} & 0.001e^{i1.52} \\ 
0.475e^{-i1.57} & 0.488 & 0.001e^{-i3.09} & 0.001e^{i3.10} \\ 
0.001e^{i1.53} & 0.001e^{i3.09} & 0.013 & 8\e{-5}e^{-i1.67} \\ 
0.001e^{-i1.52} & 0.001e^{-i3.10} & 8\e{-5}e^{i1.67} & 0.012
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.011 \\ 
0.013 \\ 
0.014 \\ 
0.962
 \end{array} $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
0.488 & 0.002e^{-i1.55} & 0.474e^{i1.57} & 0.001e^{i1.62} \\ 
0.002e^{i1.55} & 0.012 & 0.001e^{i3.07} & 1\e{-4}e^{i2.33} \\ 
0.474e^{-i1.57} & 0.001e^{-i3.07} & 0.488 & 0.001e^{-i3.10} \\ 
0.001e^{-i1.62} & 1\e{-4}e^{-i2.33} & 0.001e^{i3.10} & 0.012
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.011 \\ 
0.012 \\ 
0.015 \\ 
0.962
 \end{array} $
 \\ \hline
\end{tabular}

\caption{\textbf{Choi matrix representation of the GST estimated gate set}.  This table lists Choi representations of the estimated gates, and their eigenvalues.  Unitary gates have a spectrum $(1,0,0\ldots)$, just like pure quantum states.  Negative eigenvalues are non-physical, and may represent either statistical fluctuations or violations of the CPTP model used by GST.\label{bestCPTPGatesetChoiTable}}
\end{center}
\end{table}

\clearpage
}{}

\iftoggle{pixelplotsappendix}{
\section{Pixel plots for intermediate gate sets\label{appendix_chi2_pixelplots}}

This appendix contains $\chi^2$ pixel-plots for the intermediate estimates computed by GST on dataset ``MyDefaultReportName'' using a subset of the available data.  Each estimate was computed using only base sequences of length up to $L$, for $L=1,2,4,8,\ldots$.  These plots are \emph{not} identical to subsets of Figure \ref{bestEstimateColorBoxPlot}, because the estimated gate set is different.  Since these estimates only attempt to fit a subset of the data, they generally do much better at fitting that subset.  (They do very badly, in general, at predicting longer sequences; for this reason the $\chi^2$ values of those out-of-sample sequences are not shown).  See main text and caption of Figure \ref{bestEstimateColorBoxPlot} for more details.


\begin{figure}
\begin{center}
\includegraphics[width=2.60in,height=2.91in,keepaspectratio]{full_reportE_files/L1_Chi2Boxes.pdf}
\caption{Box plot of iteration 1 (L=1) gateset $\chi^2$ values.\label{L1GatesetChi2BoxPlot}}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[width=3.90in,height=2.91in,keepaspectratio]{full_reportE_files/L2_Chi2Boxes.pdf}
\caption{Box plot of iteration 2 (L=2) gateset $\chi^2$ values.\label{L2GatesetChi2BoxPlot}}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[width=5.20in,height=7.27in,keepaspectratio]{full_reportE_files/L3_Chi2Boxes.pdf}
\caption{Box plot of iteration 3 (L=4) gateset $\chi^2$ values.\label{L3GatesetChi2BoxPlot}}
\end{center}
\end{figure}

}{}

\iftoggle{whackamoleappendix}{
\section{Whack-a-mole plots for select gate sequences\label{appendix_whack_a_mole}}

This appendix shows whimsically named ``Whack-a-mole'' plots.  The point of these plots is to answer the question ``Why is GST failing to fit the data in the block corresponding to base sequence XXX?''  This question usually comes up after examining Figure \ref{bestEstimateColorBoxPlot} and identifying a particularly bad block.  Usually, the answer is ``Because GST was also trying to fit other data, in block YYY.''  It's often useful to know specifically \emph{which} experiments are the stumbling block -- i.e., which experiments would be even more badly fit by the estimate \emph{if} we demanded that the data in base sequence XXX be fit better.

This question can be answered by examining the derivatives of the aggregated $\chi^2$ function with respect to  the gate set, evaluated at the stationary point (local minimum) corresponding to the GST estimate.  By doing so, it is possible to simulate what would happen if, by demanding a better fit to the data in block XXX, we attempted to ``whack'' the large $\chi^2$ values in block XXX.  Unsurprisingly, whacking the XXX mole causes another one to pop up somewhere else (the GST fit is a local minimum of $\chi^2$, and hopefully a global one too, so there is no way to improve aggregate $\chi^2$).  Where this occurs can provide some insight into what is stopping GST from fitting the data better (and therefore into what's wrong with the gates).

In the Whack-a-mole plots shown below, each of the longest base sequences is independently ``whacked'' -- i.e., the analysis attempts to reduce the $\chi^2$ for that particular base sequence.  The plots answer the question, ``If the $\chi^2$ for this base sequence was forced downward by 10 units, how much would the $\chi^2$ for other experiments have to rise?''  Thus, the number on the ``whacked'' block is always $-10$.  Blocks with negative values in this analysis are correlated tightly with the whacked block (improving the fit on the whacked block improves their fit as well).  Blocks with large positive values are incompatible with the whacked block (improving the fit on the whacked block makes the fit on such a block worse).  These plots can thus be used to identify inconsistencies in the data.


\begin{figure}
\begin{center}
\includegraphics[width=6.50in,height=8.00in,keepaspectratio]{full_reportE_files/whackGxMoleBoxes.pdf}
\caption{Whack-a-$\chi^2$-mole box plot for $\mathrm{Gx}^{8}$.  Hitting with hammer of weight 10.0.\label{WhackGxMoleBoxPlot}}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[width=6.50in,height=8.00in,keepaspectratio]{full_reportE_files/whackGyMoleBoxes.pdf}
\caption{Whack-a-$\chi^2$-mole box plot for $\mathrm{Gy}^{8}$.  Hitting with hammer of weight 10.0.\label{WhackGyMoleBoxPlot}}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[width=6.50in,height=8.00in,keepaspectratio]{full_reportE_files/whackGiMoleBoxes.pdf}
\caption{Whack-a-$\chi^2$-mole box plot for $\mathrm{Gi}^{8}$.  Hitting with hammer of weight 10.0.\label{WhackGiMoleBoxPlot}}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[width=4in,height=5in,keepaspectratio]{full_reportE_files/whackGxMoleBoxesSummed.pdf}
\caption{Whack-a-$\chi^2$-mole box plot for $\mathrm{Gx}^{8}$, summed over fiducial matrix.  Hitting with hammer of weight 10.0.\label{WhackGxMoleBoxPlotSummed}}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[width=4in,height=5in,keepaspectratio]{full_reportE_files/whackGyMoleBoxesSummed.pdf}
\caption{Whack-a-$\chi^2$-mole box plot for $\mathrm{Gy}^{8}$, summed over fiducial matrix.  Hitting with hammer of weight 10.0.\label{WhackGyMoleBoxPlotSummed}}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[width=4in,height=5in,keepaspectratio]{full_reportE_files/whackGiMoleBoxesSummed.pdf}
\caption{Whack-a-$\chi^2$-mole box plot for $\mathrm{Gi}^{8}$, summed over fiducial matrix.  Hitting with hammer of weight 10.0.\label{WhackGiMoleBoxPlotSummed}}
\end{center}
\end{figure}

}{}