prompt_front: Now, answer the following question based on the text below, only give me the answer and do not output any other words.

Question: What is the proposed approach in this research paper?

Read the following text and answer briefly.

Context 1: \section{Introduction}
\label{sec:introduction}

Probabilistic models have proven to be very useful in a lot of applications in signal processing where signal estimation is needed \cite{rabiner1989tutorial,arulampalam2002tutorial,ji2008bayesian}  Some of their advantages are that 1) they force the designer to specify all the assumptions of the model, 2) they provide a clear separation between the model and the algorithm used to solve it, and 3) they usually provide some measure of uncertainty about the estimation 

On the other hand, adaptive filtering is a standard approach in estimation problems when the input is received as a stream of data that is potentially non-stationary  This approach is widely understood and applied to several problems such as echo cancellation \cite{gilloire1992adaptive}, noise cancellation \cite{nelson1991active}, and channel equalization \cite{falconer2002frequency} 

Although these two approaches share some underlying relations, there are very few connections in the literature  The first important attempt in the signal processing community to relate these two fields was the connection between a linear Gaussian state-space model (i e  Kalman filter) and the RLS filter, by Sayed and Kailath \cite{sayed1994state} and then by Haykin \emph{et al } \cite{haykin1997adaptive}  The RLS adaptive filtering algorithm emerges naturally when one defines a particular state-space model (SSM) and then performs exact inference in that model  This approach was later exploited in \cite{van2012kernel} to design a kernel RLS algorithm based on Gaussian processes 

A first attempt to approximate the LMS filter from a probabilistic perspective was presented in \cite{park2014probabilistic}, focusing on a kernel-based implementation  The algorithm of \cite{park2014probabilistic} makes use of a Maximum a Posteriori (MAP) estimate as an approximation for the predictive step  However, this approximation does not preserve the estimate of the uncertainty in each step, therefore degrading the performance of the algorithm 

In this work, we provide a similar connection between state-space models and least-mean-squares (LMS)  Our approach is based on approximating the posterior distribution with an isotropic Gaussian distribution  We show how the computation of this approximated posterior leads to a linear-complexity algorithm, comparable to the standard LMS  Similar approaches have already been developed for a variety of problems such as channel equalization using recurrent RBF neural networks \cite{cid1994recurrent}, or Bayesian forecasting \cite{harrison1999bayesian}  Here, we show the usefulness of this probabilistic approach for adaptive filtering 

The probabilistic perspective we adopt throughout this work presents two main advantages  Firstly, a novel LMS algorithm with adaptable step size emerges naturally with this approach, making it suitable for both stationary and non-stationary environments  The proposed algorithm has less free parameters than previous LMS algorithms with variable step size \cite{kwong1992variable,aboulnasr1997robust,shin2004variable}, and its parameters are easier to be tuned w r t  these algorithms and standard LMS  Secondly, the use of a probabilistic model provides us with an estimate of the error variance, which is useful in many applications 

Experiments with simulated and real data show the advantages of the presented approach with respect to previous works  However, we remark that the main contribution of this paper is that it opens the door to introduce more Bayesian machine learning techniques, such as variational inference and Monte Carlo sampling methods \cite{barber2012bayesian}, to adaptive filtering \\


\section{Probabilistic Model}

Throughout this work, we assume the observation model to be linear-Gaussian with the following distribution,

\begin{equation}
p(y_k|{\bf w}_k) = \mathcal{N}(y_k;{\bf x}_k^T {\bf w}_k , \sigma_n^2),
\label{eq:mess_eq}
\end{equation}
where  $\sigma_n^2$ is the variance of the observation noise, ${\bf x}_k$ is the regression vector and ${\bf w}_k$ is the parameter vector to be sequentially estimated, both $M$-dimensional column vectors 


In a non-stationary scenario, ${\bf w}_k$ follows a dynamic process  In particular, we consider a diffusion process (random-walk model) with variance $\sigma_d^2$ for this parameter vector:


\begin{equation}
p({\bf w}_k|{\bf w}_{k-1})= \mathcal{N}({\bf w}_k;{\bf w}_{k-1}, \sigma_d^2 {\bf I}),
\label{eq:trans_eq}
\end{equation}
where $\bf I$ denotes the identity matrix  In order to initiate the recursion, we assume the following prior distribution on ${\bf w}_k$

\begin{equation}
p({\bf w}_0)= \mathcal{N}({\bf w}_0;0, \sigma_d^2{\bf I}) \nonumber
\end{equation}

\section{Exact inference in this model: Revisiting the RLS filter}

Given the described probabilistic SSM, we would like to infer the posterior probability distribution $p({\bf w}_k|y_{1:k})$ 
Since all involved distributions are Gaussian, one can perform exact inference, leveraging the probability rules in a straightforward manner  The resulting probability distribution is
\begin{equation}
p({\bf w}_k|y_{1:k}) =  \mathcal{N}({\bf w}_k;{\bf\boldsymbol\mu}_{k}, \boldsymbol\Sigma_{k}), \nonumber
\end{equation}
in which the mean vector ${\bf\boldsymbol\mu}_{k}$ is given by
\begin{equation}
{\bf\boldsymbol\mu}_k = {\bf\boldsymbol\mu}_{k-1} + {\bf K}_k (y_k - {\bf x}_k^T {\bf\boldsymbol\mu}_{k-1}){\bf x}_k, \nonumber
\end{equation}
where we have introduced the auxiliary variable
\begin{equation}
{\bf K}_k = \frac{ \left(\boldsymbol\Sigma_{k-1} + \sigma_d^2 {\bf I}\right)}{{\bf x}_k^T  \left(\boldsymbol\Sigma_{k-1} + \sigma_d^2 {\bf I}\right)  {\bf x}_k + \sigma_n^2}, \nonumber
\end{equation}
and the covariance matrix $\boldsymbol\Sigma_k$ is obtained as
\begin{equation}
\boldsymbol\Sigma_k = \left( {\bf I} -  {\bf K}_k{\bf x}_k {\bf x}_k^T \right) ( \boldsymbol\Sigma_{k-1} +\sigma_d^2), \nonumber
\end{equation}
Note that the mode of $p({\bf w}_k|y_{1:k})$, i e  the maximum-a-posteriori estimate (MAP), coincides with the RLS adaptive rule
\begin{equation}
{{\bf w}}_k^{(RLS)} = {{\bf w}}_{k-1}^{(RLS)} + {\bf K}_k (y_k - {\bf x}_k^T {{\bf w}}_{k-1}^{(RLS)}){\bf x}_k  
\label{eq:prob_rls}
\end{equation}
This rule is similar to the one introduced in \cite{haykin1997adaptive} 

Finally, note that the covariance matrix $\boldsymbol\Sigma_k$ is a measure of the uncertainty of the estimate ${\bf w}_k$ conditioned on the observed data $y_{1:k}$  Nevertheless, for many applications a single scalar summarizing the variance of the estimate could prove to be sufficiently useful  In the next section, we show how such a scalar is obtained naturally when $p({\bf w}_k|y_{1:k})$ is approximated with an isotropic Gaussian distribution  We also show that this approximation leads to an LMS-like estimation 
 


\section{Approximating the posterior distribution: LMS filter }

The proposed approach consists in approximating the posterior distribution $p({\bf w}_k|y_{1:k})$, in general a multivariate Gaussian distribution with a full covariance matrix, by an isotropic spherical Gaussian distribution 

\begin{equation}
\label{eq:aprox_post}
\hat{p}({\bf w}_{k}|y_{1:k})=\mathcal{N}({\bf w}_{k};{\bf \hat{\boldsymbol\mu}}_{k}, \hat{\sigma}_{k}^2 {\bf I} ) 
\end{equation}

In order to estimate the mean and covariance of the approximate distribution $\hat{p}({\bf w}_{k}|y_{1:k})$, we propose to select those that minimize the Kullback-Leibler divergence with respect to the original distribution, i e , 

\begin{equation}
\{\hat{\boldsymbol\mu}_k,\hat{\sigma}_k\}=\arg \displaystyle{  \min_{\hat{\boldsymbol\mu}_k,\hat{\sigma}_k}} \{ D_{KL}\left(p({\bf w}_{k}|y_{1:k}))\| \hat{p}({\bf w}_{k}|y_{1:k})\right) \}  \nonumber
\end{equation}

The derivation of the corresponding minimization problem can be found in Appendix A  In particular, the optimal mean and the covariance are found as
\begin{equation}
{\hat{\boldsymbol\mu}}_{k} = {\boldsymbol\mu}_{k};~~~~~~ \hat{\sigma}_{k}^2 = \frac{{\sf Tr}\{ \boldsymbol\Sigma_k\} }{M} 
\label{eq:sigma_hat}
\end{equation}


We now show that by using \eqref{eq:aprox_post} in the recursive predictive and filtering expressions we obtain an LMS-like adaptive rule  First, let us assume that we have an approximate posterior distribution at $k-1$, $\hat{p}({\bf w}_{k-1}|y_{1:k-1}) =  \mathcal{N}({\bf w}_{k-1};\hat{\bf\boldsymbol\mu}_{k-1}, \hat{\sigma}_{k-1}^2 {\bf I} )$  Since all involved distributions are Gaussian, the predictive distribution
is obtained as %
\begin{eqnarray}
\hat{p}({\bf w}_k|y_{1:k-1}) &=& \int p({\bf w}_k|{\bf w}_{k-1}) \hat{p}({\bf w}_{k-1}|y_{1:k-1}) d{\bf w}_{k-1} \nonumber\\
&=& \mathcal{N}({\bf w}_k;{\bf\boldsymbol\mu}_{k|k-1}, \boldsymbol\Sigma_{k|k-1}), 
\label{eq:approx_pred}
\end{eqnarray}
where the mean vector and covariance matrix are given by
\begin{eqnarray}
\hat{\bf\boldsymbol\mu}_{k|k-1} &=& \hat{\bf\boldsymbol\mu}_{k-1} \nonumber \\
\hat{\boldsymbol\Sigma}_{k|k-1} &=& (\hat{\sigma}_{k-1}^2 + \sigma_d^2 ){\bf I}\nonumber 
\end{eqnarray}

From \eqref{eq:approx_pred}, the posterior distribution at time $k$ can be computed using Bayes' Theorem and standard Gaussian manipulations (see for instance \cite[Ch  4]{murphy2012machine})  Then, we approximate the posterior $p({\bf w}_k|y_{1:k})$ with an isotropic Gaussian,
\begin{equation}
\hat{p}({\bf w}_k|y_{1:k}) =  \mathcal{N}({\bf w}_k ; {\hat{\boldsymbol\mu}}_{k}, \hat{\sigma}_k^2 {\bf I} ),\nonumber
\end{equation}
where 
\begin{eqnarray}
{\hat{\boldsymbol\mu}}_{k} &= & {\hat{\boldsymbol\mu}}_{k-1}+ \frac{ (\hat{\sigma}_{k-1}^2+ \sigma_d^2)  }{(\hat{\sigma}_{k-1}^2+ \sigma_d^2)  \|{\bf x}_k\|^2 + \sigma_n^2} (y_k - {\bf x}_k^T {\hat{\boldsymbol\mu}}_{k-1}){\bf x}_k  \nonumber  \\
&=& {\hat{\boldsymbol\mu}}_{k-1}+ \eta_k (y_k - {\bf x}_k^T {\hat{\boldsymbol\mu}}_{k-1}){\bf x}_k   
\label{eq:prob_lms}
\end{eqnarray}
Note that, instead of a gain matrix ${\bf K}_k$ as in Eq ~\eqref{eq:prob_rls}, we now have a scalar gain $\eta_k$ that operates as a variable step size 


Finally, to obtain the posterior variance, which is our measure of uncertainty, we apply \eqref{eq:sigma_hat} and the trick ${\sf Tr}\{{\bf x}_k{\bf x}_k^T\}= {\bf x}_k^T{\bf x}_k= \|{\bf x}_k \|^2$,

\begin{eqnarray}
\hat{\sigma}_k^2 &=& \frac{{\sf Tr}(\boldsymbol\Sigma_k)}{M} \\
&=& \frac{1}{M}{\sf Tr}\left\{ \left( {\bf I} -  \eta_k {\bf x}_k {\bf x}_k^T \right) (\hat{\sigma}_{k-1}^2 +\sigma_d^2)\right\} \\
&=& \left(1 - \frac{\eta_k \|{\bf x}_k\|^2}{M}\right)(\hat{\sigma}_{k-1}^2 +\sigma_d^2) 
\label{eq:sig_k}
\end{eqnarray}
If MAP estimation is performed, we obtain  an adaptable step-size LMS estimation

\begin{equation}
{\bf w}_{k}^{(LMS)} = {\bf w}_{k-1}^{(LMS)} + \eta_k (y_k - {\bf x}_k^T {\bf w}_{k-1}^{(LMS)}){\bf x}_k, 	
\label{eq:lms}
\end{equation}
with
\begin{equation}
\eta_k = \frac{ (\hat{\sigma}_{k-1}^2+ \sigma_d^2)  }{(\hat{\sigma}_{k-1}^2+ \sigma_d^2)  \|{\bf x}_k\|^2 + \sigma_n^2} \nonumber
\end{equation}
At this point, several interesting remarks can be made:

\begin{itemize}

\item The adaptive rule \eqref{eq:lms} has linear complexity since it does not require us to compute the full matrix $\boldsymbol\Sigma_k$ 

\item For a stationary model, we have $\sigma_d^2=0$ in \eqref{eq:prob_lms} and \eqref{eq:sig_k}  In this case, the algorithm remains valid and both the step size and the error variance, $\hat{\sigma}_{k}$, vanish over time $k$  

\item Finally, the proposed adaptable step-size LMS has only two parameters, $\sigma_d^2$ and $\sigma_n^2$, (and only one, $\sigma_n^2$, in stationary scenarios) in contrast to other variable step-size algorithms \cite{kwong1992variable,aboulnasr1997robust,shin2004variable}  More interestingly, both $\sigma_d^2$ and $\sigma_n^2$ have a clear underlying physical meaning, and they can be estimated in many cases  We will comment more about this in the next section  
\end{itemize}



\section{Experiments}
\label{sec:experiments}

We evaluate the performance of the proposed algorithm in both stationary and tracking experiments  In the first experiment, we estimate a fixed vector ${\bf w}^{o}$ of dimension $M=50$  The entries of the vector are independently and uniformly chosen in the range $[-1,1]$  Then, the vector is normalized so that $\|{\bf w}^o\|=1$  Regressors $\boldsymbol{x}_{k}$ are zero-mean Gaussian vectors with identity covariance matrix  The additive noise variance is such that the SNR is $20$ dB  We compare our algorithm with standard RLS and three other LMS-based algorithms: LMS, NLMS \cite{sayed2008adaptive}, VSS-LMS \cite{shin2004variable} \footnote{The used parameters for each algorithm are: for RLS $\lambda=1$, $\epsilon^{-1}=0 01$; for LMS $\mu=0 01$; for NLMS $\mu=0 5$; and for VSS-LMS $\mu_{max}=1$, $\alpha=0 95$, $C=1e-4$ } The probabilistic LMS algorithm in \cite{park2014probabilistic} is not simulated because it is not suitable for stationary environments 

In stationary environments, the proposed algorithm has only one parameter, $\sigma^2_n$  We simulate both the scenario where we have perfectly knowledge of the amount of noise (probLMS1) and the case where the value $\sigma^2_n$ is $100$ times smaller than the actual value (probLMS2)  The Mean-Square Deviation (${\sf MSD} = {\mathbb E} \| {\bf w}_0 - {\bf w}_k \|^2$), averaged out over $50$ independent simulations, is presented in Fig  \ref{fig:msd_statationary} 



\begin{figure}[htb]
\centering
\begin{minipage}[b]{\linewidth}
  \centering
  \centerline{\includegraphics[width=\textwidth]{results_stationary_MSD}}
\end{minipage}
\caption{Performance in terms of MSD of probabilistic LMS with both optimal (probLMS1) and suboptimal (probLMS2) compared to LMS, NLMS, VS-LMS, and RLS }
\label{fig:msd_statationary}
\end{figure}

The performance of probabilistic LMS is close to RLS (obviously at a much lower computational cost) and largely outperforms previous variable step-size LMS algorithms proposed in the literature  Note that, when the model is stationary, i e  $\sigma^2_d=0$ in \eqref{eq:trans_eq},  both the uncertainty $\hat{\sigma}^2_k$, and the adaptive step size $\eta_k$, vanish over time  This implies that the error tends to zero when $k$ goes to infinity  Fig  \ref{fig:msd_statationary} also shows that the proposed approach is not very sensitive to a bad choice of its only parameter, as demonstrated by the good results of probLMS2, which uses a $\sigma^2_n$ that is $100$ times smaller than the optimal value  


\begin{figure}[htb]
\centering
\begin{minipage}[b]{\linewidth}
  \centering
  \centerline{\includegraphics[width=\textwidth]{fig2_final}}
\end{minipage}
\caption{Real part of one coefficient of the measured and estimated channel in experiment two  The shaded area represents two standard deviations from the prediction {(the mean of the posterior distribution)} }
\label{fig_2}
\end{figure}


\begin{table}[ht]
\begin{footnotesize}
\setlength{\tabcolsep}{2pt}
\def1 5mm{1 5mm}
\begin{center}
\begin{tabular}{|l@{\hspace{1 5mm}}|c@{\hspace{1 5mm}}|c@{\hspace{1 5mm}}|c@{\hspace{1 5mm}}|c@{\hspace{1 5mm}}|c@{\hspace{1 5mm}}|c@{\hspace{1 5mm}}|}
\hline
Method &  LMS &  NLMS & LMS-2013 & VSSNLMS & probLMS & RLS \\
\hline
\hline
MSD (dB) &-28 45 &-21 07 &-14 36 &-26 90 &-28 36 &-25 97\\
\hline                                                                     
\end{tabular}
\end{center}
\caption{Steady-state MSD of the different algorithms for the tracking of a real MISO channel }
\label{tab:table_MSD}
\end{footnotesize}

\end{table}
\newpage
In a second experiment, we test the tracking capabilities of the proposed algorithm with {real} data of a wireless MISO channel acquired in a realistic indoor scenario  More details on the setup can be found in \cite{gutierrez2011frequency}  Fig  \ref{fig_2} shows the real part of one of the channels, and the estimate of the proposed algorithm  The shaded area represents the estimated uncertainty for each prediction, i e  $\hat{\mu}_k\pm2\hat{\sigma}_k$  Since the experimental setup does not allow us to obtain the optimal values for the parameters, we fix these parameters to their values that optimize the steady-state mean square deviation (MSD)  \hbox{Table \ref{tab:table_MSD}} shows this steady-state MSD of the estimate of the MISO channel with different methods  As can be seen, the best tracking performance is obtained by standard LMS and the proposed method  





\section{Conclusions and Opened Extensions}
\label{sec:conclusions}

{We have presented a probabilistic interpretation of the least-mean-square filter  The resulting algorithm is an adaptable step-size LMS that performs well both in stationary and tracking scenarios  Moreover, it has fewer free parameters than previous approaches and these parameters have a clear physical meaning  Finally, as stated in the introduction, one of the advantages of having a probabilistic model is that it is easily extensible:}

\begin{itemize}
\item If, instead of using an isotropic Gaussian distribution in the approximation, we used a Gaussian with diagonal covariance matrix, we would obtain a similar algorithm with different step sizes and measures of uncertainty, for each component of ${\bf w}_k$  Although this model can be more descriptive, it needs more parameters to be tuned, and the parallelism with LMS vanishes 
\item Similarly, if we substitute the transition model of \eqref{eq:trans_eq} by an Ornstein-Uhlenbeck process, 

\begin{equation}
p({\bf w}_k|{\bf w}_{k-1})= \mathcal{N}({\bf w}_k;\lambda {\bf w}_{k-1}, \sigma_d^2), \nonumber
\label{eq:trans_eq_lambda}
\end{equation}
a similar algorithm is obtained but with a forgetting factor $\lambda$ multiplying ${\bf w}_{k-1}^{(LMS)}$ in \eqref{eq:lms}  This algorithm may have improved performance under such a kind of autoregresive dynamics of ${\bf w}_{k}$, though, again, the connection with standard LMS becomes dimmer 

\item As in \cite{park2014probabilistic}, the measurement model \eqref{eq:mess_eq} can be changed to obtain similar adaptive algorithms for classification, ordinal regression, and Dirichlet regression for compositional data  

\item A similar approximation technique could be applied to more complex dynamical models, i e  switching dynamical models \cite{barber2010graphical}  The derivation of efficient adaptive algorithms that explicitly take into account a switch in the dynamics of the parameters of interest is a non-trivial and open problem, though the proposed approach could be useful 

\item Finally, like standard LMS, this algorithm can be kernelized for its application in estimation under non-linear scenarios 

\end{itemize}


\begin{appendices}

\section{KL divergence between a general gaussian distribution and an isotropic gaussian}
\label{sec:kl}

 We want to approximate  $p_{{\bf x}_1}(x) = \mathcal{N}({\bf x}; \boldsymbol\mu_1,\boldsymbol\Sigma_1)$ by $p_{{\bf x}_2}({\bf x}) = \mathcal{N}({\bf x}; \boldsymbol\mu_2,\sigma_2^2 {\bf I})$  In order to do so, we have to compute the parameters of $p_{{\bf x}_2}({\bf x})$, $\boldsymbol\mu_2$ and $\sigma_2^2$, that minimize the following Kullback-Leibler divergence,

\begin{eqnarray}
D_{KL}(p_{{\bf x}_1}\| p_{{\bf x}_2}) &=&\int_{-\infty}^{\infty} p_{{\bf x}_1}({\bf x}) \ln{\frac{p_{{\bf x}_1}({\bf x})}{p_{{\bf x}_2}({\bf x})}}d{\bf x} \nonumber  \\
&= &  \frac{1}{2} \{ -M + {\sf Tr}(\sigma_2^{-2} {\bf I}\cdot \boldsymbol\Sigma_1^{-1})  \nonumber \\
  & &  + (\boldsymbol\mu_2 - \boldsymbol\mu_1 )^T \sigma^{-2}_2{\bf I} (\boldsymbol\mu_2 - \boldsymbol\mu_1 )  \nonumber \\
 & &   +  \ln \frac{{\sigma_2^2}^M}{\det\boldsymbol\Sigma_1} \}   
\label{eq:divergence}
\end{eqnarray}
Using symmetry arguments, we obtain 
\begin{equation}
\boldsymbol\mu_2^{*} =\arg \displaystyle{  \min_{\boldsymbol\mu_2}} \{ D_{KL}(p_{{\bf x}_1}\| p_{{\bf x}_2}) \} = \boldsymbol\mu_1 
\end{equation}
Then, \eqref{eq:divergence} gets simplified into 

\begin{eqnarray}
D_{KL}(p_{{\bf x}_1}\| p_{{\bf x}_2}) = \frac{1}{2}\lbrace { -M + {\sf Tr}(\frac{\boldsymbol\Sigma_1}{\sigma_2^{2}}) + \ln \frac{\sigma_2^{2M}}{\det\boldsymbol\Sigma_1}}\rbrace 
\end{eqnarray}
The variance $\sigma_2^2$ is computed in order to minimize this Kullback-Leibler divergence as

\begin{eqnarray}
\sigma_2^{2*} &=& \arg\min_{\sigma_2^2} D_{KL}(P_{x_1}\| P_{x_2}) \nonumber \\
 &=& \arg\min_{\sigma_2^2}\{ \sigma_2^{-2}{\sf Tr}\{\boldsymbol\Sigma_1\} + M\ln \sigma_2^{2} \}  
\end{eqnarray}
Deriving and making it equal zero leads to

\begin{equation}
\frac{\partial}{\partial \sigma_2^2} \left[ \frac{{\sf Tr}\{\boldsymbol\Sigma_1\}}{\sigma_2^{2}} + M \ln \sigma_2^{2} \right] = \left  {\frac{M}{\sigma_2^{2}}-\frac{{\sf Tr}\{\boldsymbol\Sigma_1\}}{(\sigma_2^{2})^2}}\right|_{\sigma_2^{2}=\sigma_2^{2*}}\left  =0 \right   
\nonumber
\end{equation}
Finally, since the divergence has a single extremum in $R_+$,
\begin{equation}
\sigma_2^{2*} = \frac{{\sf Tr}\{\boldsymbol\Sigma_1\}}{M} 
\end{equation}




\end{appendices}

\vfill
\clearpage

\bibliographystyle{IEEEbib}

Context 2: Paper Info

Title: Compressed quantum error mitigation
Publish Date: 10 May 2023
Author List: Maurits Tepaske (from Physikalisches Institut, Universität Bonn), David Luitz (from Physikalisches Institut, Universität Bonn)

Figure

FIG 3 The out-of-time-ordered correlator C otoc i=L/2,j (t) as a function of the operator position j and time t, for the infinite temperature initial state, for a denoised second-order Trotter supercircuit with Trotter depth Mtrot = 32 and denoiser depth M = 2 We consider evolution times t = 0 5, 1,    , 5, for the periodic L = 14 Heisenberg chain that is affected by two-qubit depolarizing noise with p = 0 01 
FIG  4  The complex eigenvalues λ of the noisy second-order Trotter supercircuit with Mtrot = 16 at time t = 1 (left), the corresponding optimized denoiser with M = 4 (center), and the denoised Trotter supercircuit (right) The Trotter circuit is for a L = 6 Heisenberg model with PBC, and all twoqubit channels are affected by depolarizing noise with p = 0 0046 The unit circle, on which unitary eigenvalues must lie, is shown in black, and the noiseless eigenvalues are shown as blue bars It is evident that the denoiser recovers all the noiseless eigenvalues from the noisy circuit 
FIG  2  The complex eigenvalues λ of the noisy second-order Trotter supercircuit with Mtrot = 16 at time t = 1 (left), the corresponding optimized denoiser with M = 4 (center), and the denoised Trotter supercircuit (right) The Trotter circuit is for a L = 6 Heisenberg model with PBC, and all twoqubit channels are affected by depolarizing noise with p = 0 036 The unit circle, on which unitary eigenvalues must lie, is shown in black, and the noiseless eigenvalues are shown as blue bars It is clear that the denoiser recovers with high accuracy the noiseless eigenvalues from the noisy circuit 
FIG  3  The half-chain channel entanglement entropy S at different two-qubit depolarizing noise strengths p, for a secondorder Trotter supercircuit with Mtrot = 16 and t = 2, for a M = 4 denoiser The Trotter circuit is for a Heisenberg model with PBC of size L = 6 The different curves correspond to the different supercircuits, i e  the noisy supercircuit, the denoiser, the corresponding denoised supercircuit, and the noiseless variant 
FIG  4  The out-of-time-ordered correlator C otoc i=L/2,j (t) as a function of the operator position j and stacked time t, for the infinite temperature initial state, for a denoised secondorder Trotter supercircuit with Trotter depth Mtrot = 32 and denoiser depth M = 2 It is optimized at t = 2 and stacked up to ten times The calculations are for the periodic L = 14 Heisenberg chain that is affected by two-qubit depolarization with p = 0 01 The denoiser is affected by the same noise 
FIG 6 The distribution of the ZZ angle α of M = 2 denoisers (top panels) and M = 8 denoisers (bottom panels), with the lightest color corresponding to the denoiser for the Trotter supercircuit with t = 0 5, and the darkest color with t = 5 As usual, we consider the Heisenberg model on a periodic chain, and second-order Trotter supercircuits with depths Mtrot = 8, 16, 32, 64, which together with the denoiser is affected by a two-qubit depolarizing noise with p = 0 01 The panels are arranged as Mtrot = 8, 16, 32, 64 for top left, top right, bottom left, bottom right, respectively 
FIG  7  The sampling overhead γ of the optimized denoisers from Fig  2 of the main text, with denoiser depths M = 1, 2, 4, 6, 8 and Trotter depths Mtrot = 8, 16, 32, 64 at times t = 0 5, 1,    , 5, for the Heisenberg model on a chain with PBC affected by two-qubit depolarizing noise with p = 0 01 The panels are arranged as Mtrot = 8, 16, 32, 64 for top left, top right, bottom left, bottom right, respectively 
FIG 8 The domain wall magnetization Z dw after evolving a periodic density wall |dw |dw * with the denoised second-order Trotter supercircuits D C from Fig 2of the main text These supercircuits have various Trotter depths Mtrot = 8, 16, 32, 64, denoiser depths M = 1, 2, 4, 6, 8, and evolution times t = 0 5, 1,    , 5, for the periodic L = 14 Heisenberg chain that is affected by two-qubit depolarizing noise of strength p = 0 01 The denoiser is affected by the same noise The non-denoised results are labelled with M = 0 and the noiseless results with p = 0 The panels are arranged as Mtrot = 8, 16, 32, 64 for top left, top right, bottom left, bottom right, respectively We see that the denoiser allows us to recover the noiseless behavior 

abstract

We introduce a quantum error mitigation technique based on probabilistic error cancellation to eliminate errors which have accumulated during the application of a quantum circuit  Our approach is based on applying an optimal "denoiser" after the action of a noisy circuit and can be performed with an arbitrary number of extra gates 
The denoiser is given by an ensemble of circuits distributed with a quasiprobability distribution  For a simple noise model, we show that efficient, local denoisers can be found, and we demonstrate their effectiveness for the digital quantum simulation of the time evolution of simple spin chains  Introduction 
-Quantum information processing has been theoretically shown to hold great promises, and quantum algorithms were developed which can in principle achieve an exponential speed-up over their classical counterparts, both for general purpose computing and quantum simulation   However, present day quantum computing prototypes still suffer from significant noise processes which hinder the execution of many potentially groundbreaking quantum algorithms  
Nontrivial quantum algorithms typically require large sequences of quantum gates, each of which introduces dissipation and hence an overall loss of coherence, eventually rendering the results useless  Until quantum error correction becomes practical, quantum error mitigation seems to be more feasible to increase the accuracy of expectation values 
Here the goal is to induce the (partial) cancellation of errors that stem from noisy quantum gates by extending the circuit corresponding to the desired algorithm with an ensemble of gates , sampled from a quasiprobability distribution  The traditional way to accomplish this is with the gatewise method from , where noise is mitigated by inverting the noise channel of each gate separately, i e  the cancellation of errors is performed for each gate on its own 
Here the local noise channel is approximated in a way such that it can be easily inverted analytically, e g  using Pauli twirling   Gates are then sampled from the inverted noise channel by interpreting it as a quasiprobability distribution  Because in this gate-wise approach every noisy gate has to be modified separately, the sign problem is exponentially large in the number of gates, limiting the practicality of the mitigation 
The success of the gate-wise approach resulted in a large body of work concerning these methods , including extensions for simultaneous mitigation of multiple gates by Pauli-twirling entire layers or variationally constructing a mitigating matrix product operator   In principle, errors during the execution of a circuit can propagate and accumulate 
These propagated errors * david luitz@uni-bonn de ≈ C

C

FIG  1  An example of the quantum error mitigation procedure used in this work for the time evolution of the wave function of a spin chain  The ideal second-order Trotter supercircuit C of depth Mtrot = 1 (light blue) is approximated by applying a denoiser D of depth M = 1 (red) to the noisy Trotter supercircuit C (dark blue) 
Because the denoiser is applied after fully executing the noisy Trotter supercircuit, it represents an approximate inverse of the global noise channel with a precision tunable by the depth of the denoiser  can potentially blow up and lead to large errors for the circuit as a whole   Here we introduce a mitigation technique that takes into account the propagation of errors, can be performed with a tunable number of extra gates, and works for non-Clifford local noise channels since the inversion of the accumulated global noise channel is implicit 
We first execute the targeted noisy circuit completely, letting the noise propagate and accumulate, and only afterwards we apply an extra random circuit sampled from a quasiprobability distribution  We call the corresponding ensemble of random circuits a denoiser, and we construct it such that upon averaging the accumulated errors cancel 
Essentially, the denoiser inverts a global noise channel  Since we will construct it as a local brickwall circuit, following the classical preprocessing approach from , we call this compressed quantum error mitigation  Method  -Due to the inevitable coupling of a quantum processor to its environment, every qubit operation is affected by noise 
Therefore, the simplest technique to minimize the impact of the resulting noise is to minimize the number of operations when performing a quantum algorithm  In we showed that many-body time evolution operators can be efficiently compressed into brick-wall circuits with high fidelity per gate  In this Letter, we consider the noise explicitly by treating quantum operations as (generally non-unitary) quantum channels, corresponding to completely positive and trace preserving (CPTP) maps  
For example, instead of a noiseless two-qubit gate G, which acts on a quantum state |ρ in superoperator form as G|ρ = G⊗G * |ρ , we get the noisy channel G = N G, where the noise channel N implements the two-qubit noise   These channels are used to construct a "supercircuit" C = N G i=1 Gi , consisting of N G channels, which is affected by multi-qubit accumulated noise 
This supercircuit encodes an ensemble of circuits   For simplicity, we assume that the noisy channels Gi in each half brickwall layer are lattice inversion and translation invariant, such that we can construct a denoiser with these properties, limiting the number of variational parameters  The purpose of quantum error mitigation is to modify the ensemble of circuits described by C in a way that we can use it to obtain the noiseless expectation values 
In superoperator language, we do this by following the supercircuit C with a denoiser supercircuit D, such that D C is as close to the noiseless supercircuit C = C ⊗ C * as possible  Here C is the target unitary circuit  Because the noise channel N is non-unitary, hence making the supercircuit C non-unitary, we need to use a non-unitary denoiser to retrieve the unitary C 
We illustrate the mitigation procedure in Fig  , where a denoiser with one layer is used to mitigate errors for a second-order Trotter supercircuit with one layer  This circuit architecture is commonly used to simulate the time evolution of a quantum many-body system, until some time t, with controllable precision , and we will use it to benchmark the denoiser 
In practice, we cannot directly implement a supercircuit, and so we have to utilize its interpretation as an ensemble of circuits  Essentially, after executing a shot of the noisy circuit we sample the denoiser and apply it  The goal is to construct the denoiser in a way that averaging over many of its samples cancels the accumulated errors and gives us a good approximation of the noiseless expectation values 
It should be noted that our approach requires more gate applications on the quantum processor than with the gate-wise scheme, since there each sample from the mitigation quasiprobability distribution can be absorbed into the original circuit, whereas our approach increases the circuit depth  We take this into account by imposing the same noise on the denoiser 
Furthermore, within our scheme, the dimensionality of the quasiprobabilistic mitigating ensemble can be controlled, in contrast to the gate-wise approach where it is equal to the gate count  To facilitate the stochastic interpretation we parameterize each two-qubit denoiser channel G i as a sum of CPTP maps, such that we can sample the terms in this sum and execute the sampled gate on the quantum processor 
Concretely, we use a trace preserv-ing sum of a unitary and a non-unitary channel  For the unitary part we take a two-qubit unitary channel U( φ i ) = U ( φ i ) ⊗ U * ( φ i ), with U ( φ i ) a two-qubit unitary gate parameterized by φ i   For this we take the two-qubit ZZ rotation exp(−iα(σ z ⊗ σ z )) with angle α, which can be obtained from native gates on current hardware , and dress it with four general one-qubit unitaries, only two of which are independent if we want a circuit that is space inversion symmetric around every bond 
The resulting gate has 7 real parameters φ i   For the non-unitary part, which is essential because D has to cancel the non-unitary accumulated noise to obtain the noiseless unitary circuit, we use a general onequbit measurement followed by conditional preparation channel M( , with V a general one-qubit unitary and each κ i a 3-dimensional vector, resulting in a real 9-dimensional ζ i  
This yields the two-qubit correlated measurement M( With these parts we construct the parameterization with coefficients η i ∈ R that satisfy η 0 + η 1 = 1 because G i is trace preserving  Note that here the tensor product symbol corresponds to combining two one-qubit channels to make a two-qubit channel, whereas in most of the paper it is used to link the column and row indices of a density matrix 
We construct the denoiser from the noisy channels Gi = N G i   With this parameterization one denoiser channel has 17 independent real parameters, such that a denoiser of depth M , i e  consisting of M brickwall layers, has 34M real parameters (we use one unique channel per half brickwall layer)  For reference, a general channel has 544M parameters 
To determine the mitigated expectation values we use the full expression where |ρ 0 is the initial state and |1 is the vectorized identity operator on the full Hilbert space  To evaluate this on a quantum processor, we use the stochastic interpretation of (1) to resample   In particular, from each channel (1) we get a unitary with probability p 0 = |η 0 |/γ and a measurement followed by conditional preparation with probability p 1 = |η 1 |/γ 
Here γ = |η 0 | + |η 1 | is the sampling overhead, which characterizes the magnitude of the sign problem from negative η i   For quasiprobability distributions, i e  with γ > 1, every denoiser sample has an extra sign sgn(η) = N G g=1 sgn(η g ), 2  The normalized distance between the denoised Trotter supercircuit D C and the noiseless Trotter supercircuit C (top panels), at evolution times t = 0 5, 1,    , 5, and the twopoint z-spin correlator C zz i=L/2,j=L/2 (t) of a spin on the middle site at times 0 and t (bottom panels), for the infinite temperature initial state 
We consider denoisers with depths M = 1, 2, 4, 6, 8 and second-order Trotter circuits with depths Mtrot = 16, 32, 64  In the top panels we use a Heisenberg chain with L = 8, and in the bottom panels with L = 14, both with periodic boundary conditions  All gates are affected by two-qubit depolarizing noise with p = 0 01
Context 3: Sir Simon William English  (born 30 December 1961) is a New Zealand former National Party politician who served as the 39th prime minister of New Zealand from 2016 to 2017  He had previously served as the 17th deputy prime minister of New Zealand and minister of finance from 2008 to 2016 under John Key and the Fifth National Government 

A farmer and public servant before entering politics, English was elected to the New Zealand Parliament in  as the National Party's candidate in the Wallace electorate  He was elevated to Cabinet in 1996 and in 1999 was made minister of finance, although he served for less than a year due to his party's loss at the 1999 general election  In October 2001, English replaced Jenny Shipley as the leader of the National Party (and consequently as Leader of the Opposition)  He led the party to its worst defeat at the 2002 general election, and as a consequence, in October 2003 he was replaced as leader by Don Brash 

In November 2006, after Brash's resignation, English became deputy leader under John Key  After National's victory at the 2008 general election, he became deputy prime minister and was also made minister of finance for the second time  Under English's direction New Zealand's economy maintained steady growth during National's three terms of government  He became a list-only MP after stepping down as an electorate MP at the 2014 general election 

John Key resigned as leader of the National Party and prime minister in December 2016  English won the resulting leadership election unopposed and was sworn in as prime minister on 12 December 2016  His tenure was only ten months, and included a three-month election campaign  In the 2017 general election, National won the largest number of seats but fell short of a majority  The parties holding the balance of power declined to support the existing government, and English was subsequently replaced as prime minister by Jacinda Ardern, leader of the Labour Party  English initially continued on as Leader of the Opposition, but resigned as leader of the National Party on 27 February 2018 and left parliament two weeks later 

Early life
English was born on 30 December 1961 at Lumsden Maternity Centre in Lumsden  He is the eleventh of twelve children of Mervyn English and Norah (née O'Brien) English  His parents purchased Rosedale, a mixed sheep and cropping farm in Dipton, Southland from Mervyn's uncle, Vincent English, a bachelor, in 1944  English was born in the maternity unit at Lumsden 

English attended St Thomas's School in Winton, then boarded at St  Patrick's College in Upper Hutt, where he became head boy  He played in the first XV of the school's rugby team  English went on to study commerce at the University of Otago, where he was a resident at Selwyn College, and then completed an honours degree in English literature at Victoria University of Wellington 

After finishing his studies, English returned to Dipton and farmed for a few years  From 1987 to 1989, he worked in Wellington as a policy analyst for the New Zealand Treasury, at a time when the free market policies favoured by Labour's finance minister Roger Douglas (known collectively as "Rogernomics") were being implemented 

English joined the National Party in 1980, while at Victoria University  He served for a period as chairman of the Southland branch of the Young Nationals, and became a member of the Wallace electorate committee  After moving to Wellington, he served for periods on the Island Bay and Miramar electorate committees, respectively 

Fourth National Government (1990–1999)

At the 1990 general election, English stood as the National candidate in Wallace, replacing the retiring Derek Angus, and was elected with a large majority  He would hold this seat, renamed Clutha-Southland in 1996, until 2014  He and three other newly elected National MPs (Tony Ryall, Nick Smith, and Roger Sowry) were soon identified as rising stars in New Zealand politics, and at various points were dubbed the "brat pack", the "gang of four", and the "young Turks"  In his first term in parliament, English chaired a select committee into social services  He was made a parliamentary under-secretary in 1993, serving under the Minister of Health 

First period in cabinet (1996–1999)
In early 1996, English was elevated to cabinet by Prime Minister Jim Bolger, becoming the Minister for Crown Health Enterprises and Associate Minister of Education (to Wyatt Creech)  He was 34 at the time, becoming the cabinet's youngest member  After the 1996 general election, the National Party was forced into a coalition with New Zealand First to retain government  In the resulting cabinet reshuffle, English emerged as Minister of Health  However, as a condition of the coalition agreement, NZ First's Neil Kirton (a first-term MP) was made Associate Minister of Health, effectively becoming English's deputy  This arrangement was described in the press as a "shotgun marriage", and there were frequent differences of opinion between the two ministers  After their relationship became unworkable, Kirton was sacked from the role in August 1997, with the agreement of NZ First leader Winston Peters 

As Minister of Health, English was responsible for continuing the reforms to the public health system that National had begun after the 1990 general election  The reforms were unpopular, and health was perceived as one of the government's weaknesses, with the health portfolio consequently being viewed as a challenge  English believed that the unpopularity of the reforms was in part due to a failure in messaging, and encouraged his National colleagues to avoid bureaucratic and money-focused language (such as references to "balance sheets" and "user charges") and instead talk about the improvements to services the government's reforms would bring  He also rejected the idea that public hospitals could be run as commercial enterprises, a view which some of his colleagues had previously promoted 

By early 1997, as dissatisfaction with Bolger's leadership began to grow, English was being touted as a potential successor, along with Jenny Shipley and Doug Graham  His age (35) was viewed as the main impediment to a successful leadership run  National's leadership troubles were resolved in December 1997, when Bolger resigned and Shipley was elected to the leadership unopposed  English had been a supporter of Bolger as leader, but Shipley reappointed him Minister of Health in her new cabinet 

English was promoted to Minister of Finance in a reshuffle in January 1999, a position which was at the time subordinate to the Treasurer, Bill Birch  After a few months, the pair switched positions as part of Birch's transition to retirement, with English assuming the senior portfolio  In early interviews, he emphasised his wish to be seen as a pragmatist rather than an ideologue, and said that the initiatives of some of his predecessors (Roger Douglas's "Rogernomics" and Ruth Richardson's "Ruthanasia") had focused on "fruitless, theoretical debates" when "people just want to see problems solved" 

Opposition (1999–2008)

After the National Party lost the 1999 election to Helen Clark's Labour Party, English continued on in the shadow cabinet as National's spokesperson for finance  He was elected deputy leader of the party in February 2001, following the resignation of Wyatt Creech, with Gerry Brownlee being his unsuccessful opponent 

Leader of the Opposition
In October 2001, after months of speculation, Jenny Shipley resigned as leader of the National Party after being told she no longer had the support of the party caucus  English was elected as her replacement unopposed (with Roger Sowry as his deputy), and consequently became Leader of the Opposition  However, he did not openly organise against Shipley, and according to The Southland Times "there was almost an element of 'aw, shucks, I'll do it then' about Mr English's ascension" 

Aged 39 when he was elected, English became the second-youngest leader in the National Party's history, after Jim McLay (who was 38 when elected in 1984)  He also became only the third Southlander to lead a major New Zealand political party, after Joseph Ward and Adam Hamilton  However, English failed to improve the party's performance  In the 2002 election, National suffered its worst electoral defeat ever, gaining barely more than twenty percent of the vote  English described it as "the worst day of my political life"  Both party insiders and the general public were split as to how much to blame him for the loss, but most of the party believed that English would be able to rebuild National's support 

By late 2003, however, National's performance in opinion polls remained poor  The party had briefly increased its popularity in the year following the election, but by October its support had fallen to levels only slightly better than what it achieved in the last ballot  English also appeared in a boxing match for a charity against entertainer Ted Clarke  This did not boost his polling or that of the National party either, with suggestions that it devalued his image as a serious politician  Don Brash, former governor of the Reserve Bank and a relative newcomer to politics, began to build up support to replace English  On 28 October, Brash gained sufficient backing in Caucus to defeat English in a leadership contest 

Shadow cabinet roles and deputy leader
On 2 November 2003, when Brash changed responsibilities for certain MPs, English became National's spokesman for education, ranked at fifth place in the party's parliamentary hierarchy  He remained in parliament after the 2005 election  In his new shadow education portfolio, English performed strongly, and remained a party favourite despite his election defeat as leader in 2002, eventually being returned to the finance portfolio in August 2004 as deputy spokesman (while still retaining responsibility for education) 

In November 2006, Brash resigned as leader  English was considered as a potential replacement leader (running against John Key) or deputy leader (against incumbent Gerry Brownlee) in the ensuing leadership election  However, a contest was avoided when the MPs agreed a Key/English ticket would run unopposed in a display of party unity  English took over the deputy leadership and the finance portfolio in the Key shadow cabinet 

Fifth National Government (2008–2017)

Deputy Prime Minister and Minister of Finance (2008–2016)

At the 2008 election, English was re-elected by his electorate, winning by a margin of about 15,500 votes  He became Deputy Prime Minister of New Zealand and Minister of Finance in the fifth National Government, being sworn into office on 19 November 2008 and continued to serve in those roles until becoming Prime Minister on 12 December 2014  He was also made Minister of Infrastructure in National's first term of government and Minister responsible for Housing New Zealand Corporation and minister responsible for the New Zealand flag consideration process in its third 

He was comfortably re-elected in Clutha-Southland in the 2011 election but opted to run as a party-list candidate in 2014  

The pairing of John Key as leader of the National Party and English as his deputy has been compared to that of Bob Hawke and Paul Keating (in Australia) and Tony Blair and Gordon Brown (in the UK)



ground_truth: This research paper proposed an approach based on approximating the posterior distribution with an isotropic Gaussian distribution.

prompt_front: Now, answer the following question based on the text below, only give me the answer and do not output any other words.

Question: What size chains were used in the benchmarking?

Read the following text and answer briefly.

Context 2: 
The non-denoised results are labelled with M = 0, and the noiseless values with p = 0  where sgn(η g ) is the sign of the sampled coefficient of the gth channel  γ = 1 means that all signs are positive  Observables Ô p=0 for the noiseless circuit are then approximated by resampling the observables from the denoiser ensemble
where γ = N G g=1 γ g is the overall sampling overhead, with γ g the overhead of the gth gate  Clearly, a large γ implies a large variance of Ô p=0 for a given number of samples, with accurate estimation requiring the cancellation of large signed terms  The number of samples required to resolve this cancellation of signs is bounded by Hoeffding's inequality, which states that a sufficient number of samples to estimate Ô p=0 with error δ at probability 1 − ω is bounded by (2γ 2 /δ 2 ) ln(2/ω)  
Since γ scales exponentially in γ g , it is clear that a denoiser with large M and γ 1 will require many samples  We observed that decompositions with γ > 1 are crucial for an accurate denoiser  Restricting to γ = 1 leads to large infidelity and no improvement upon increasing the number of terms in or the depth M of the denoiser 
Simply put, probabilistic error cancellation of gate noise introduces a sign problem and it is crucial to find optimal parameterizations (1) which minimize γ to make the approach scalable  This issue arises in all high performance error mitigation schemes , because the inverse of a physical noise channel is unphysical and cannot be represented as a positive sum over CPTP maps 
This is clearly visible in the spectra of the denoiser, which lies outside the unit circle (cf  Fig  )  This makes the tunability of the number of gates in each denoiser sample a crucial ingredient, which allows control over the sign problem, because we can freely choose the η i in   For the parametrization (1) of denoiser channels, we try to find a set of parameters for error mitigation by minimizing the normalized Frobenius distance between the noiseless and denoised supercircuits
which bounds the distance of output density matrices and becomes zero for perfect denoising  We carry out the minimization of on a classical processor, using gradient descent with the differential programming algorithm from   Instead of explicitly calculating the accumulated global noise channel and subsequently inverting it, we approximate the noiseless supercircuit C with the denoised supercircuit D C, effectively yielding a circuit representation D of the inverse noise channel 
Results  -To benchmark the denoiser we apply it to the second-order Trotter circuits of the spin-1/2 Heisenberg chain with periodic boundary conditions (PBC) where is the Pauli algebra acting on the local Hilbert space of site i  A second-order Trotter circuit for evolution time t with depth M trot consists of M trot − 1 half brickwall layers with time step t/M trot and two layers with half time step  
We consider circuits that are affected by uniform depolarizing noise with probability p for simplicity, but our approach can be used for any non-Clifford noise  The two-qubit noise channel is which acts on neighboring qubits i and i + 1 and is applied to each Trotter and denoiser gate, and p = 0 01 unless stated otherwise 
We study circuits with depths M trot = 16, 32, 64 for evolution times t = 0 5, 1,    , 5, and denoisers D with depths M = 1, 2, 4, 6, 8  In the top panels of Fig  we show (4) for a chain of size L = 8 as a function of time t  Here it can be seen that even for M trot = 32 a denoiser with M = 1 already improves by roughly an order of magnitude at all considered t 
Depending on M trot and t, further increasing M lowers , with the biggest improvements occurring for high precision Trotter circuits with large depth M trot = 64 and short time t = 0 5, where the Trotter gates are closer to the identity than in the other cases  At the other extreme, for M trot = 16 the improvements are relatively small upon increasing M > 2  In all cases the denoiser works better at early times than at late times, again indicating that it is easier to denoise Trotter gates that are relatively close to the identity 
To probe the accuracy of the denoiser on quantities that do not enter the optimization, as a first test we consider the two-point correlator between spins at different times where we have chosen the infinite temperature initial state, and C(t) is the Trotter supercircuit for time t  In the bottom panels of Fig  we show C zz i=L/2,j=L/2 (t) for the supercircuits from the upper panels, now for a L = 14 chain 
Here we see that at M trot = 16 we can retrieve the noiseless values already with M = 1, but that increasing M trot makes this more difficult  At M trot = 64 we see larger deviations, and improvement upon increasing M is less stable, but nonetheless we are able to mitigate errors to a large extent  As a further test, we compute the out-of-time-ordered correlator (OTOC) ]
In Fig  we show the results for i = L/2, for a Trotter circuit with depth M trot = 32 and a denoiser with depth M = 2  Here we see that a denoiser with M M trot is able to recover the light-cone of correlations, which are otherwise buried by the noise  In the Supplementary Material we consider how the denoiser performs at different noise levels p, and how the denoised supercircuits perform under stacking 
There we also calculate domain wall magnetization dynamics, and show the distribution of the optimized denoiser parameters and the sampling overhead associated to the denoiser as a whole  In Fig  we show the eigenvalues of the noisy supercircuits for a noisy second-order Trotter supercircuit with M trot = 16 at t = 1 (left), the corresponding optimized denoiser with M = 4 (center), and the denoised supercircuit (right) 
The eigenvalues λ of a unitary supercircuit lie on the unit circle, and in the presence of dissipation they are pushed to the center  We see that the spectrum of the denoiser lies outside the unit circle, making it an unphysical channel which cures the effect of the noise on the circuit, such that the spectrum of the denoised circuit is pushed back to the unit circle 
The noiseless eigenvalues are shown as blue bars, making it clear that the denoiser is able to recover the noiseless eigenvalues from the noisy circuit  In the Supplementary Material we show the spectra for a p = 0 036 denoiser, where we observe a clustering of eigenvalues reminiscent of Refs    There we also investigate the channel entropy of the various supercircuits  
Conclusion  -We have introduced a probabilistic error cancellation scheme, where a classically determined denoiser mitigates the accumulated noise of a (generally non-Clifford) local noise channel  The required number of mitigation gates, i e  the dimensionality of the corresponding quasiprobability distribution, is tunable and the parameterization of the corresponding channels provides control over the sign problem that is inherent to probabilistic error cancellation 
We have shown that a denoiser with one layer can already significantly mitigate errors for second-order Trotter circuits with up to 64 layers  This effectiveness of low-depth compressed circuits for denoising, in contrast with the noiseless time evolution operator compression from , can be understood from the non-unitarity of the denoiser channels 
In particu-lar, measurements can have non-local effects, since the measurement of a single qubit can reduce some highly entangled state (e g  a GHZ state) to a product state, whereas in unitary circuits the spreading of correlations forms a light-cone  To optimize a denoiser with convenience at L > 8, the optimization can be formulated in terms of matrix product operators or channels , which is convenient because the circuit calculations leading to the normalized distance and its gradient are easily formulated in terms of tensor contractions and singular value decompositions  
This provides one route to a practical denoiser, which is relevant because the targeted noiseless circuit and the accompanying noisy variant in (4) need to be simulated classically, confining the optimization procedure to limited system sizes with an exact treatment or limited entanglement with tensor networks 
Nonetheless, we can use e g  matrix product operators to calculate (4) for some relatively small t, such that the noiseless and denoised supercircuits in (4) have relatively small entanglement, and then stack the final denoised supercircuit on a quantum processor to generate classically intractable states 
Analogously, we can optimize the channels exactly at some classically tractable size and then execute them on a quantum processor with larger size  Both approaches are limited by the light-cone of many-body correlations, as visualized in Fig  , because finite-size effects appear when the light-cone width becomes comparable with system size 
1  The normalized distance (left) and z spin correlator C zz i=L/2,j=L/2 (right), for a second-order Trotter supercircuit of depth Mtrot = 16 for time t = 1, affected by various twoqubit depolarizing errors p  We compare the values obtained with and without a denoiser, i e  M > 0 and M = 0, to the noiseless values (p = 0) 
The denoiser is affected by the same noise as the Trotter circuit  We consider denoisers with depths M = 1, 2, 4, 6, 8, and we use a L = 8 Heisenberg chain with PBC for the normalized distance, while for the correlator we use L = 14  * david luitz@uni-bonn de to observe that even for larger noise strength p, the local observable C zz improves significantly even with denoisers of depth M = 1 
For large noise strengths, we generally see that the optimization of the denoiser becomes difficult, leading to nonmonotonic behavior as a function of p, presumably because we do not find the global optimum of the denoiser  It is interesting to analyze the spectra of the supercircuits considered in this work 
As mentioned in the main text, the spectrum of the ideal, unitary supercircuit C lies on the unit circle  The comparison to this case is therefore instructive  In the main text, we showed an example of the spectra in Fig  for moderate noise strength  Here, we show additional data for stronger noise p = 0 036 in Fig  for a denoiser with M = 4 layers, optimized to mitigate errors for a second-order Trotter supercircuit with M trot = 16 layers at time t = 1 
The eigenvalues λ of the noisy supercircuit C are clustered close to zero, far away from the unit circle (except for λ = 1), showing that the circuit is strongly affected by the noise  To mitigate the impact of the noise, the denoiser consequently has to renormalize the spectrum strongly  If it accurately represents the inverse of the global noise channel, its spectrum has to lie far outside the unit circle, which is the case 
Interestingly, we observe a clustering of eigenvalues which is reminiscent to the spectra found in   By comparison to these works, we suspect that this is due to the local nature of the denoiser, and warrants further investigation  The right panel of Fig  shows the result of the denoiser, pushing the eigenvalues back to the unit circle, nearly with the exact same distribution along the circle as the noiseless eigenvalues (blue bars) 
Due to the strong noise, this is not achieved perfectly, and it is clear that this cannot work in principle if the global noise channel has a zero eigenvalue  The complexity of an operator can be quantified by its operator entanglement entropy   Here we calculate the half-chain channel entanglement entropy S of the noiseless C, noisy C, denoiser D, and denoised D C supercircuits 
We define S as the entanglement entropy of the state that is related to a supercircuit C via the Choi-Jamio lkowski isomorphism, i e  ψ C = χ C /N , where the process matrix χ ab,cd C = C ac,bd is simply a reshaped supercircuit and N ensures normalization  Then we have S = −Tr [ψ C ln ψ C ]  This entropy measure is a particular instance of the "exchange entropy", which characterizes the information exchange between a quantum system and its environment  
In Fig  we plot the various S for a second-order Trotter circuit with M trot = 16 at t = 2, for a denoiser with M = 4, both affected by two-qubit depolarizing noise with p ∈ [10 −3 , 10 −1 ]  The Trotter circuit is for a Heisenberg model with L = 6 and PBC  We see that at large p, the noise destroys entanglement in the noisy supercircuit, and that the denoiser S increases to correct for this, such that the denoised supercircuit recovers the noiseless S 
Here we investigate how denoised supercircuits perform upon repeated application  We optimize the denoiser for a Trotter supercircuit for a fixed evolution time t  Then, to reach later times, we stack the denoised supercircuit n times to approximate the evolution up to time nt: In Fig  we stack a denoised t = 1 supercircuit up to n = 20 times and calculate the correlation function, defined in the main text, for the middle site 
We consider Trotter depths M trot = 8, 16, 32, 64 and denoiser depths M = 1, 2, 4, 6, 8, for a L = 14 Heisenberg chain with p = 0 01 depolarizing two-qubit noise  The noisy results correspond to M = 0 and the noiseless results to p = 0  In Fig  we calculate the OTOC, defined in the main text, with stacked time evolution for a denoised t = 2 supercircuit with M trot = 32 and M = 2, stacked up to ten times 
We see that the stacked supercircuit performs very well, and the additional precision obtained by using deep denoisers (M = 8) pays off for long evolution times, where we see convergence to the exact result (black dashed lines in Fig  ) as a function of M   FIG    The two-point z-spin correlator C zz i=L/2,j=L/2 (t) of a spin on the middle site at times 0 and t, for the infinite temperature initial state, for denoised second-order Trotter supercircuits that are optimized at evolution time t = 1 and then stacked up to twenty times 
We use Trotter depths Mtrot = 8, 16, 32, 64 and denoiser depths M = 1, 2, 4, 6, 8  The calculations were performed for a periodic Heisenberg model with L = 14 and PBC, affected by two-qubit depolarizing noise with strength p = 0 01, which also affects the denoiser  The non-denoised results are labelled with M = 0, and the noiseless results with p = 0 
The panels are arranged as Mtrot = 8, 16, 32, 64 for top left, top right, bottom left, bottom right, respectively  The costliest and most noise-susceptible operation is the two-qubit ZZ rotation with angle α, which is the foundation of the unitary piece in our channel parameterization, defined in the main text 
For completeness, we here present the α angles of the optimized denoisers  The results are shown in Fig  , which contains histograms for the channel count N G versus α  The histograms are stacked, with the lightest color corresponding to the angles of the denoiser at t = 0 5 and the darkest at t = 5  The top four panels are for a denoiser with M = 2 and the bottom four with M = 8 
We consider M trot = 8, 16, 32, 64  We see that in both cases the distribution widens upon increasing M trot , indicating that the unitary channels start deviating more from the identity  Moreover, while the M = 2 denoisers in all cases except M trot = 64 have ZZ contributions close to the identity, this is clearly not the case for M = 8 
For simplicity, we did not focus on obtaining denoisers with the smallest sampling overhead γ, which is required to minimize the sign problem and hence ease the sampling of mitigated quantities  Instead, we let the optimization freely choose the η i in the denoiser parameterization, as defined in the main text 
In Fig  we show the sampling overhead of the denoisers from Fig  of the main text  We see that for M = 1 and M = 2 the sampling overhead is relatively small and uniform across the different t, whereas for M > 2 the optimization sometimes yields a denoiser with large γ and other times with small γ  This could be related to the difference in α distributions from Fig   
The large fluctuations of γ appears to stem from the difficulty in finding optimal deep denoisers, and our optimization procedure likely only finds a local minimum in these cases  Here C(t) is the Trotter supercircuit for time t  In Fig  we show Z dw for the circuits from Fig 



ground_truth: L = 8 and L = 14.

prompt_front: Now, answer the following question based on the text below, only give me the answer and do not output any other words.

Question: In which electorate was Simon English elected to the New Zealand Parliament?

Read the following text and answer briefly.

Context 3: 

English acceded to the role of Finance Minister in the continuing wake of the financial crisis  In response to New Zealand's rising debt, English made budget deficit-reduction his main priority  His first budget outlined three focuses in New Zealand's financial recovery: "improving the business environment and removing roadblocks to growth; investment in productive infrastructure; and improving the way government works"  One of his first acts was creating the National Infrastructure Unit, charged with formulating a plan for infrastructure projects and investments  He commissioned a government-wide spending review, with an aim to reducing government expenditure—with the exceptions of a two-year stimulus package and long-term increases on infrastructure spending 

In April 2011, the Opposition criticised English for suggesting that New Zealand businesses could use New Zealand's low wages to help it compete with Australia  The National Government campaigned for re-election in 2011 on its economic record  The Government boasted growth for five consecutive quarters up to mid-2010, totalling 1 6% of real GDP 

Strong growth resulted in a surplus of $473 million for the 2015/16 financial year, projected to rise to $8 5 billion by 2020/21  In his 2016 Economic and Fiscal Update address, English stated that reducing debt and tackling the costs of the 2016 Kaikōura earthquake were higher priorities than reducing rates of tax 

Allowances issue
In 2009, the media, including TVNZ and TV3 revealed that English was receiving about NZ$900 a week as part of a living allowance for ministers, to live in his own NZ$1 2 million Wellington home  At the time, English also received $276,200 in his annual salary as Deputy Prime Minister  It was also revealed other ministers with homes in the capital city were also claiming accommodation allowances  On 3 August 2009, Prime Minister John Key started a review of the housing allowances claimed by cabinet ministers  English subsequently paid back $12,000 and only claimed about $24,000 a year in living allowances  The Auditor-General's office said in September 2009 that they were making "preliminary enquiries" into parliamentary housing expenses in response to a letter of complaint from Progressive party leader Jim Anderton  Two days later English stated that he would no longer take up any housing allowance and had paid back all the allowance he had received since the November 2008 election 

Prime Minister (2016–2017)

John Key resigned on 12 December, and endorsed English as his successor in the resulting leadership election  Following the drop-out of both Judith Collins and Jonathan Coleman from the leadership election, English was sworn in as the 39th Prime Minister of New Zealand on 12 December 2016 

English appointed his first cabinet on 18 December  In a reshuffle, he appointed Steven Joyce to succeed him as Finance Minister, while most ministerial portfolios remained the same 

In February 2017, English did not attend Waitangi Day commemorations at the historic treaty grounds, reportedly in response to the Ngāpuhi iwi's decision to stop the Prime Minister from speaking at the marae  Ngāpuhi have protested the Government's negotiation of the Trans Pacific Partnership Agreement (TPPA), which the iwi believe infringes upon Māori sovereignty, and thus does not adhere to the Treaty of Waitangi  English had been invited to attend in an official capacity; his non-attendance was criticised by a Ngāpuhi elder and Opposition leader Andrew Little 

In his first overseas trip as Prime Minister, English travelled to Europe to discuss trade ties, including a prospective New Zealand–European Union free trade agreement  He first travelled to London on 13 January 2017 to meet British Prime Minister Theresa May  Discussing trade relations, English said the two nations were "natural partners" and would "continue to forge ties" after the UK's withdrawal from the EU  He also arranged to meet with London Mayor Sadiq Khan, Belgian Prime Minister Charles Michel and German Chancellor Angela Merkel   In a meeting with Merkel, English received crucial backing from Germany for a trade deal with the EU  On 16 January, English stated that his government would continue to promote TPPA, despite the United States' decision to withdraw from the agreement  He explained that Southeast Asian countries would now be treated as a priority in negotiations—he also asserted that the United States was ceding influence to China by its rejection of the trade pact 

At a press conference at the Beehive on 1 February 2017, English announced that the 2017 general election would be held on 23 September  The Prime Minister later confirmed that his party would approach ACT, United Future and the Māori Party if confidence and supply agreements were required to form a government following the election  In his second cabinet reshuffle on 24 April, English appointed Gerry Brownlee as his new Foreign Affairs Minister; he also promoted Nikki Kaye to the portfolio of Education Minister, and moved Mark Mitchell into the cabinet to become Defence Minister  The reshuffle was perceived as an election preparation 

On 13 February 2017, English welcomed Australian Prime Minister Malcolm Turnbull to Wellington  The two leaders reaffirmed their shared trade agenda, and discussed changes to the Australian citizenship pathway which will affect permanent residents originating from New Zealand 

On 19 June, it was reported that Todd Barclay, who succeeded English as MP for Clutha-Southland, had clandestinely recorded one of his employee's conversations the previous year, and that John Key's leaders' budget was used to pay a confidential settlement after the employee resigned  English admitted that he had been aware of the illegal recording and the settlement, and thus implicated in the scandal 

During the 2017 National campaign launch, English introduced a $379 million social investment package including digital learning academies for high school students, more resources for mathematics, and boosting support for teaching second languages in schools, and maintaining National Standards in the school curriculum  Prime Minister English also sought to defend National's financial management and economic track record and claimed that the opposition Labour Party would raise taxes  Early opinion polling had forecast a poor showing in the election for the Labour Party, but in early August 37-year-old Jacinda Ardern took over as Labour leader and seemingly energised younger voters 

At the 2017 general election, National won the largest share of the party vote (44 4%) and the largest number of seats (56) in the House Representatives  However, National lacked enough seats to govern alone due to two of the party's support partners, the Māori Party and United Future, losing their parliamentary seats  In response, English stated that the party would be entering into talks to form a coalition with New Zealand First  Following talks with the two largest parties, New Zealand First entered a coalition arrangement with the Labour Party  English was succeeded as prime minister by Jacinda Ardern on 26 October 

Opposition (2017–2018)

Leader of the Opposition
English was re-elected as National Party leader on 24 October 2017  At the time of his re-election, English announced his intention to stay on as leader until the next general election  On 13 February 2018, however, he stood down as National Party leader due to personal reasons, and instructed the party to put into motion the processes to elect a new leader  He also retired from Parliament  English's resignation followed weeks of speculation that he would step aside for a new leader  On 27 February, he was succeeded as party leader by Simon Bridges as the result of the leadership election held that day 

Post-premiership 
In 2018, English joined the board of Australian conglomerate, Wesfarmers  English serves in Chairmanships of Mount Cook Alpine Salmon, Impact Lab Ltd and Manawanui Support Ltd  He is also a director of The Instillery, Centre for Independent Studies and The Todd Corporation Limited, and is a member of the Impact Advisory Group of Macquarie Infrastructure and Real Assets 

Political and social views

English is regarded as more socially conservative than his predecessor, John Key  He has stated his opposition to voluntary euthanasia and physician-assisted suicide, same-sex civil unions, and the decriminalisation of prostitution  As Prime Minister he opposed any "liberalisation" of abortion law 

In 2004, English voted against a bill to establish civil unions for both same-sex and opposite-sex couples  In 2005, he voted for the Marriage (Gender Clarification) Amendment Bill, which would have amended the Marriage Act to define marriage as only between a man and a woman  English voted against the Marriage (Definition of Marriage) Amendment Bill, a bill that legalised same-sex marriage in New Zealand  However, in December 2016 he stated, "I'd probably vote differently now on the gay marriage issue  I don't think that gay marriage is a threat to anyone else's marriage" 

In 2009, English voted against the Misuse of Drugs (Medicinal Cannabis) Amendment Bill, a bill aimed at amending the Misuse of Drugs Act so that cannabis could be used for medical purposes 

Personal life 
English met his future wife, Mary Scanlon, at university  She was studying medicine at the time, and became a general practitioner  Both her parents were immigrants, her father being Samoan and her mother Italian, born on the island of Stromboli  They have six children: a daughter and five sons 

English is a practising Roman Catholic, but has stated that he considers his religious beliefs personal and thus separate from politics 

In June 2002, English took part in TV3's Fight For Life, a celebrity boxing fundraiser to raise money for the Yellow Ribbon anti-youth-suicide campaign, influenced by the death of a teenage nephew in 1997  He lost a split decision to former university colleague Ted Clarke 

Honours
In the 2018 Queen's Birthday Honours, English was appointed a Knight Companion of the New Zealand Order of Merit, for services of over 27 years to the State 

See also

List of New Zealand governments
Politics of New Zealand

References

External links

Profile at National Party 
Profile on Parliament nz
Releases and speeches at Beehive govt nz

|-

|-

|-

|-

|-

|-

|-

|-

|-

|-

|-

|-

|-

|-

|-

|-

|-

|-

1961 births
21st-century New Zealand politicians
Candidates in the 2017 New Zealand general election
Deputy Prime Ministers of New Zealand
Leaders of the Opposition (New Zealand)
Living people
Members of the Cabinet of New Zealand
Members of the New Zealand House of Representatives
New Zealand farmers
New Zealand finance ministers
New Zealand list MPs
New Zealand MPs for South Island electorates
New Zealand National Party MPs
New Zealand National Party leaders
New Zealand Roman Catholics
New Zealand people of Irish descent
People educated at St  Patrick's College, Silverstream
People from Dipton, New Zealand
People from Lumsden, New Zealand
Prime Ministers of New Zealand
University of Otago alumni
Victoria University of Wellington alumni
Knights Companion of the New Zealand Order of Merit
New Zealand politicians awarded knighthoods



ground_truth: The Wallace electorate.