the discovery of a primordial magnetic field ( pmf ) would have a profound impact on our understanding of the early universe @xcite and would help explain the origin of the observed magnetic fields in galaxies and clusters @xcite .
several observational probes are currently being investigated and the cosmic microwave background ( cmb ) is a promising tool to discover and study the pmf on the largest cosmic scales .    a pmf can affect the thermal distribution and the polarization of the cmb .
the most relevant cmb signature depends on the form of the magnetic field spectrum and also on the level of instrumental noise for different observations .
for instance , current cmb bounds of a few ng @xcite on a scale - invariant pmf @xcite derive from temperature anisotropies sourced by the magnetic stress - energy .
comparable bounds were recently obtained in @xcite using the lyman-@xmath0 forest spectra @xcite .
as the stress - energy is quadratic in the magnetic field strength , improving this bound by a factor of @xmath1 would require a @xmath2 fold improvement in the accuracy of the spectra .
on the other hand , faraday rotation ( fr ) of cmb polarization , being linear in the magnetic field strength , offers an alternative way to improve cmb bounds on a scale - invariant pmf by an order of magnitude @xcite , and possibly more , with the next generation cmb experiments .
even current cmb data is close to providing competitive bounds on scale - invariant pmf through their fr signature @xcite .
in contrast , a pmf produced causally in an aftermath of the electroweak or qcd phase transition @xcite would have a blue spectrum @xcite with most of the power concentrated near a small cutoff scale set by the plasma conductivity at recombination .
it has been suggested @xcite that such small - scale fields can appreciably alter the recombination history via enhancement of small scale baryonic inhomogeneities .
consequently , the strongest cmb constraint would come from an overall shift in the distance to last scattering . according to @xcite
, upcoming cmb experiments can rule out causally produced fields with a comoving strength larger than @xmath3 .
given the promise of fr to significantly improve the bounds on scale - invariant pmf @xcite , it is important to assess the strength of the rotation induced by the magnetic fields in our own galaxy .
a rotation measure ( rm ) map of the milky way has recently been assembled by oppermann et al @xcite based on an extensive catalog of fr of compact extragalactic polarized radio sources .
they also presented the angular power spectrum of rm and find that it is very close to a scale- invariant spectrum at @xmath4 .
1 of @xcite shows that the root - mean - square rm away from the galactic plane is about @xmath5-@xmath6 .
additional insight into galactic rm on smaller scales can be gained from haverkorn et al @xcite who studied stokes parameters generated by diffuse polarized sources residing inside the milky way .
their power spectra generally agree with the scale - invariance of the large scale tail of the power spectrum of @xcite , but also show evidence of a break in the spectrum , indicating a power law suppression of the rm power on small scales .
such a break in the rm spectrum is also present in the model of minter and spangler @xcite who , by fitting to a small set of rm data from extragalactic sources , argued that the spectrum should be set by kolmogorov turbulence on small scales .
the approximately scale - invariant shape of the galactic rm spectrum can obscure fr constraints on scale - invariant pmf . based on calculations in @xcite
, one can estimate that a typical galactic rm of @xmath6 corresponds to an _ effective _ , or `` energy equivalent '' ( see eq .
( [ beff ] ) for the definition ) pmf strength of @xmath7 ng .
this is well below the bounds from planck @xcite obtained using non - fr diagnostics and , as we confirm in this paper , below the levels detectable by planck via the mode - coupling statistics induced by fr .
however , as shown in @xcite , future cmb experiments capable of detecting the weak lensing b - mode will be able to constrain a scale - invariant pmf of @xmath8 ng strength at @xmath9 ghz . operating at lower frequencies and combining fr information from several channels may further improve the constraints .
clearly , with such high sensitivities the contribution of the galactic rm to the total fr signal will become important .    in this paper
, we investigate the imprint of the galactic rm on cmb observables , and its impact on detectability of the pmf via the eb and tb mode - coupling correlations .
we start by introducing the necessary concepts and reviewing the known galactic rm measurements in sec .
[ frcmb ] . in sec .
[ sec : detection ] , we estimate detectability of the galactic rm by upcoming and future cmb experiments and forecast future bounds on the scale - invariant pmf under various assumptions . we conclude with a summary in sec .  [
sec : summary ] .
a cmb experiment measures stokes parameters in different directions on the sky , with parameters @xmath10 and @xmath11 quantifying linear polarization . if cmb photons pass through ionized regions permeated by magnetic fields , the direction of linear polarization is rotated by an angle @xcite @xmath12 where @xmath13 is the direction along the line of sight , @xmath14 is the differential optical depth , @xmath15 is the observed wavelength of the radiation , @xmath16 is the `` comoving '' magnetic field strength ( the physical field strength scales with the expansion as @xmath17 ) and @xmath18 is the comoving length element along the photon trajectory . the rotation measure , @xmath19 , is a frequency independent quantity used to describe the strength of fr . under the rotation of the polarization vector ,
the two stokes parameters transform as @xmath20 \includegraphics[height=0.46\textwidth]{clbb } \caption{the cmb b - mode spectrum from faraday rotation sourced at $ 30$~ghz   by a scale - invariant primordial magnetic field of 1 ng strength ( red solid ) , by the full sky   galactic magnetic field ( blue dot ) , and by the galactic field with planck 's sky   mask ( $ f_{\rm sky}=0.6 $ ) .
the black short - dash line is the input e - mode spectrum ,   the black dash - dot line   is the contribution from inflationary gravitational waves with $ r=0.1 $ , while the black   long - dash line is the expected contribution from gravitational lensing by large scale   structure . } \label{fig : cl } \end{figure }    the quadratic estimator method is based on the assumption of statistical isotropy of primordial perturbations .
this implies statistical independence of primordial $ e_{l m}$ and $ b_{l m}$ for different $ lm$ pairs , { \it e.~g . }
$ \langle e_{l m}^ * e_{l'm ' } \rangle =   \delta_{ll ' } \delta_{mm ' } c_l^{ee}$ for the primordial e - mode .
furthermore , if primordial fields are gaussian , all of their correlation functions can be expressed in terms of the power spectra .
fr introduces correlations between unequal $ lm$ and generates connected four - point correlations .
the corresponding non - gaussian signal can be used to extract the rotation angle .
namely , given a cmb polarization map , one constructs quantities such as \cite{kamionkowski:2008fp , gluscevic:2009 mm } \be { \hat d}_{ll'}^{lm,{\rm map}}={4\pi \over ( 2l+1)(2l'+1 ) } \sum_{mm ' } b_{lm}^{\rm map}e_{l'm'}^{\rm map * } \xi_{lml'm'}^{lm } \label{eq : dll } \ee which is the minimum variance unbiased estimator for $ d_{ll'}^{lm } = 2 \alpha_{lm}c_l^{ee}h_{ll'}^l$ as shown   in \cite{pullen:2007tu}. then , each pair of $ l$ and $ l'$ provides an estimate of $ \alpha_{lm}$ via \be [ { \hat \alpha}_{lm}]_{ll ' } = { { \hat d}_{ll'}^{lm,{\rm map } } \over 2c_l^{ee}h_{ll'}^l } \ . \label{alphallpr } \ee the minimum variance estimator $ { \hat \alpha}_{lm}$ is obtained from an appropriately weighted sum over estimators for each $ ll'$ pair .
if all such pairs were statistically independent , the weighting would be given by the inverse variance of each estimator .
however , one has to account for the correlation   between the $ ll'$ and $ l'l$ pairs and a detailed derivation of the appropriate weighting can be found in \cite{kamionkowski:2008fp}. quantities such as $ { \hat d}_{ll'}^{lm,{\rm map}}$ can also be constructed from products of t and b or , more generally , all possible quadratic combinations , { \it i.e. } \{eb , be , tb , bt , te , et , ee\}. one can construct an estimator $ { \hat \alpha}_{lm}$ that utilizes all these combinations while accounting for the covariance between them \cite{gluscevic:2009mm}. in this paper , we opt to consider eb / be and tb / bt separately , since one of them is typically   much more informative than other combinations . for further details of the method the reader
is referred to   \cite{kamionkowski:2008fp , gluscevic:2009mm}.    we will assume a `` stochastic '' primordial magnetic field such that its value at any given location is described by a probability distribution function .
this implies that one can only   predict statistical properties of fr , such as the rotation power spectrum ,   $ c_l^{\alpha \alpha}$. the estimator for the power spectrum directly follows   from the estimator of $ \alpha_{lm}$ and derives from four - point correlations ,   such as ebeb . in the next section
, we will examine detectability of the rotation   spectrum caused by the galactic and primordial magnetic fields .
we note that the form of the quadratic estimator allows for contributions from the   monopole ( $ l=0 $ ) and dipole ( $ l=1 $ ) of the fr field which , in principle , should not be   ignored .
however , the monopole is generally not expected for fr , since it would imply   a non - zero magnetic charge enclosed by the cmb surface , while the dipole ,   corresponding to a uniform magnetic field , is strongly constrained for the pmf and is   relatively unimportant for the galactic rm .
in addition to mode - coupling correlations of eb and tb type , fr also contributes   to the b - mode polarization spectrum , $ c_l^{bb}$. in fig.~\ref{fig : cl } we show the b- mode spectrum at $ 30 $ ghz due to fr by a pmf with a scale - invariant spectrum and   effective field strength of $ 1~{\rm ng}$.the b - mode   spectrum due to fr by the galactic field is also shown for two choices of sky cuts .
for   reference , we show the $ e$ mode auto - correlation spectrum which acts as a source   for the fr b modes , the b modes from inflationary gravitational waves with $ r=0.1 $ ,   and the expected contribution from weak lensing ( wl ) . as can be seen from fig.~ \ref{fig : cl } , the shape of the fr induced b mode spectrum largely mimics that of the e- mode , and the galactic contribution has a very similar shape to that of the scale- invariant pmf . a more detailed discussion of the fr induced b - mode spectrum can be   found in \cite{pogosian:2011qv}. as shown in sec.~\ref{sec : detection } , the   fr signatures of the galactic or scale - invariant pmf are always less visible in   the b - mode spectrum   compared to the eb quadratic estimator , with the latter promising to provide the   tightest constraints .
we note that for causally generated pmf , which have blue spectra   with most of the power concentrated at the dissipation scale , the signal - to - noise ratio   ( snr ) for detecting fr in bb is higher than that in the quadratic estimators   \cite{yadav:2012uz } , although the resulting constraints are not competitive with   other probes of   causal fields . in
what follows , we will focus on fr signatures of scale - invariant pmf   only .
\subsection{a scale - invariant primordial magnetic field } \label{pmf }    assuming statistical homogeneity and isotropy , the   magnetic field correlation function in fourier space is~\cite{1975mit .. bookr .... m } \begin{eqnarray } \langle { b}^ * _ i ( { \bm k } ) { b } _ j ( { \bm k } ' ) \rangle & = & \left [ \left(\delta _ { ij } - { k_i k_j \over k^2 } \right )    s ( k )   + i \epsilon_{ijl } { k_l \over k } a(k ) \right ] \nonumber \\ & & \hskip 2 cm \times ( 2 \pi)^3 \delta^3 ( { \bm k } -{\bm k } ' )   \end{eqnarray } where repeated indices are summed over , $ b_i({\bm k})$ denotes the fourier transform of the magnetic field at wave vector $ { \bm k}$ , $ k=|{\bm k}|$ , and $ s(k)$ and $ a(k)$ are the symmetric and   antisymmetric ( helical ) parts of the magnetic field power spectrum .     the helical part of the power spectrum does not play a role in faraday rotation ;   only the symmetric part of the spectrum is relevant for us . on scales larger than the inertial scale ,
{ \it i.e. } $ k < k_i$ , $ s(k)$ will fall off as a power law . on scales   smaller than a dissipative scale , { \it i.e. }
$ k > k_d$ , the spectrum will get sharply cut - off .
a second power - law behavior is possible for $ k_i < k <
k_d$.   however , we shall assume that $ k_i \approx k_d$ , so that the magnetic field spectrum is given by a single power law behavior at all scales larger than the   dissipation scale .
these features can be summarized by \cite{pogosian:2011qv } \begin{eqnarray } s(k ) = \begin{cases } \omega_{b\gamma } { \rho}_\gamma                \frac{32\pi^3 n}{k_i^3 }      \left ( \frac{k}{k_i } \right ) ^{2n-3 } , & \mbox{$0<k < k_{\rm diss}$}\\ 0 , & \mbox{$k_{\rm diss } < k$ }         \end{cases } \label{eq : pb } \end{eqnarray } where $ \omega_{b\gamma}$ is the ratio of cosmological magnetic and photon   energy densities , and $ \rho_\gamma$ is the photon energy density .    for convenience , we define an `` effective '' , or `` energy equivalent '' , magnetic field strength \begin{equation } b_{\rm eff } \equiv ( 8\pi \omega_{b\gamma } \rho_\gamma ) ^{1/2 }   = 3.25\times 10^{-6}\sqrt{\omega_{b\gamma } } ~{\rm gauss}. \label{beff}\ ] ] a uniform magnetic field of strength @xmath21 will have the same energy density as a stochastic magnetic field with spectrum in eq .
( [ eq : pb ] ) .
scale invariance for @xmath22 corresponds to @xmath23 , in which case the expressions in eq .
( [ eq : pb ] ) are not well defined .
one reason is that the energy density diverges logarithmically for a strictly scale invariant field .
we get around this issue by only considering approximately scale invariant magnetic fields for which @xmath24 is small .
note that , for scale invariant fields , @xmath21 is equal to the field strength smoothed over a scale @xmath25 ( if it is larger than the dissipation scale ) , or @xmath26 , frequently used in the literature ( _ e.g. _ in @xcite ) .
as we have already noted , the galactic rm spectrum is approximately scale invariant .
the rm induced by a scale invariant pmf is also scale invariant @xcite .
hence the amplitudes of the galactic and pmf rotation measure spectra can be related directly for a scale invariant pmf .    for a scale - invariant angular spectrum @xmath27 ,
the quantity @xmath28 is constant .
hence , one can define an amplitude @xmath29 for the range of @xmath30 over which the rm spectrum is approximately scale - invariant .
one can see from fig .
[ fig : rmcl ] that this holds for the galactic @xmath19 for @xmath31 .
similarly , one can define an amplitude a_= the relation between the faraday rotation angle in radians and the rm is given by @xmath32 where @xmath33 .
for example , a characteristic galactic rm of @xmath34 at @xmath35 gives a rotation angle of @xmath36 rad , or @xmath37 .
the amplitudes of @xmath38 in radians and @xmath39 in @xmath40 are therefore also related via a_= 10 ^ -4 ^-2_30 a_rm  m^2  .
[ alpharm ]    to gain some intuition into the relative importance of the galactic rm prior to more detail forecasts in sec  [ detect - pmf ] , let us estimate the @xmath21 of scale - invariant pmf that gives rm comparable to that of our galaxy . for this , we can compare the amplitude @xmath38 from pmf to the amplitude for galactic spectra at the same frequency .
the relation between pmf and @xmath38 involves integration along the line of sight ( eq .  ( [ alpha - fr ] ) ) . rather than trying to find an approximate analytical solution ,
we use exact spectra found numerically in @xcite . from fig .  1 of ref
.  @xcite , one can infer that for a scale invariant pmf with @xmath41 , @xmath42 rad@xmath43 at @xmath35 .
we can then use eqs .
( [ alpharm ] ) and eq .
( [ beff ] ) to relate @xmath21 to @xmath39 : b_eff 0.021  a_rm  ng  m^2 rad  .
[ eq : beffrm ] from fig .
[ fig : rmcl ] , we can see that @xmath44 for the galactic rm spectrum obtained without any sky cuts .
it follows from eq .
( [ eq : beffrm ] ) that this corresponds to a scale - invariant pmf with @xmath45 .
such pmf are well below levels detectable by planck through their fr signal , but can be detected at high significance by future sub - orbital and space - based experiments @xcite .
we can also estimate the peak value of the b - mode spectrum , @xmath46 , in ( @xmath47k)@xmath43 , for a given @xmath39 and frequency . from fig .
[ fig : cl ] , we see that @xmath48k)@xmath49 for a scale invariant pmf of @xmath50  ng strength . using this fact along with eq .
( [ eq : beffrm ] ) , we can deduce p_b ^2 ( k)^2 ^2 ( k)^2  .
thus , for a typical galactic @xmath39 of @xmath51 rad / m@xmath43 , at @xmath51  ghz , we should expect a peak in the b - mode spectrum of about @xmath52 ( @xmath47k)@xmath43 , which agrees with the blue dotted line in fig .
[ fig : cl ] corresponding to the unmasked galactic rm . in practice ,
the galactic plane is masked in all cmb experiments , _
e.g. _ the mask applied to planck s @xmath51  ghz map only uses @xmath7 of the sky .
this significantly reduces the amplitude of the galactic rm spectrum and , as a consequence , the b - mode spectrum .
this is shown as a blue dashed line in fig .
[ fig : cl ] .
one can see that the galactic b - mode is at least one , and more likely two , orders of magnitude below the weak lensing signal , so there is almost no hope of seeing it even at a frequency as low as @xmath51  ghz .
going a factor of @xmath53 higher in frequency , which is more typical in cmb studies , results in another factor of @xmath54 decrease in the b - mode spectrum .
one of the conclusions of this section is that the galactic rm is practically invisible in the bb auto - correlation .
hence , we should look at the mode - coupling estimators , as they are known to be much more sensitive to fr @xcite .
as shown in @xcite , future experiments can easily probe @xmath55 strength pmf fields with the mode - coupling estimators .       sky cut .
the bottom panel also shows lines of symmetric cuts corresponding to @xmath56 and @xmath57 .
, scaledwidth=85.0% ]    , obtained from the rm map of oppermann et al @xcite with different cuts . shown
are the rm spectra corresponding to , from top to bottom , a scale - invariant pmf of @xmath50 ng , galaxy with no sky cut , with a mask used by planck for their @xmath58 ghz map , a planck mask for the @xmath51 ghz map , and symmetric cuts corresponding to @xmath56 and @xmath57 .
, scaledwidth=40.0% ]    a rotation measure ( rm ) map of milky way has been put together by oppermann et al @xcite based on an extensive catalog of faraday rotation data of compact extragalactic polarized radio sources .
the rm is strongest along the galactic plane , but remains significant at the poles .
however , the latitude dependence of the rms rm disappears at latitudes higher than about @xmath59 , where the dependence flattens out at a rms value of about @xmath60 rad / m@xmath43 .
oppermann et al @xcite also presented the angular power spectrum of rm calculated after factoring out the latitude dependence .
they find that @xmath61 at @xmath62 , which is very close to a scale - invariant spectrum .
some insight into galactic rm on smaller scales is provided by haverkorn et al @xcite .
they studied stokes parameters , generated by diffuse polarized sources residing inside the galaxy , at several frequencies around 350 mhz .
they focused on two separate patches away from the galactic plane , each covering about @xmath63 sq .
the rotation measures in both were less than @xmath60 rad / m@xmath43 .
however , since they are looking at sources inside the galaxy , they only see a fraction of the total rm along each line of sight .
thus , their numbers can not be directly compared to rm of @xcite .
still their power spectrum gives a sense for the scale dependence of the total rm , and generally agrees with the flat nature of the large scale tail of the power spectrum of @xcite . at a scale of about @xmath64 ( @xmath65 ) , they see evidence of a break in the spectrum , indicating a power law suppression of the rm power on small scales .
such a break in the spectrum is also present in the model of minter and spangler @xcite , whose analysis is based on a small set of rm data from extragalactic sources .
they argued that on small scales the spectrum should be due to kolmogorov turbulence .
the corresponding scaling of @xmath66 would be @xmath67 . to summarize , based on the evidence in @xcite , if we were to measure rm in a few degree patches near galactic poles
, we would expect a rm spectrum that is roughly scale invariant at @xmath68 turning into a @xmath67 tail on smaller scales . as we will show in the next section ,
the quadratic estimator of the fr angle is largely uninformative for @xmath69 .
hence , we do not need a rm map with a resolution higher than that of oppermann et al @xcite .
note that the cmb polarization even in a small patch around the poles is affected by very long wavelength modes of variations of the rm . in terms of eq .
( [ eq : blm ] ) , lowest multipoles of @xmath70 contribute significantly to the much higher multipoles of @xmath71 .
thus , to predict the observability of the fr in cmb in a small patch one can not simply use the high @xmath72 portion of the rm spectrum .
we should use the full rm spectrum obtained from the rm after applying a mask normally used to block the galactic plane .    in fig .
[ fig : rmmap ] we show the rm map for our galaxy obtained from the faraday rotation of extragalactic sources in units of rad / m@xmath73 by oppermann et al @xcite .
the top panel is the full map .
the mid - panel shows the uncertainty in rm in rad / m@xmath73 .
most of the signal and uncertainty come from the galactic disk .
the bottom panel shows the rm map with part of the sky removed with a mask used for planck s mask @xmath51 ghz map @xcite .
the bottom panel also shows lines of symmetric latitude cuts corresponding to sky fractions of @xmath57 and @xmath74 .    in fig .
[ fig : rmcl ] we present the power spectra @xmath75 of the galactic rm under different symmetric sky cuts and two planck mask cuts : one for @xmath51 ghz with @xmath76 , and one for @xmath58 ghz with @xmath77 .
the figure also shows the rm spectrum from a scale - invariant pmf with @xmath78 ng . to calculate the galactic rm spectra ,
we first use healpix  @xcite to find _ lm = d  w ( )  r ( )  y^*_lm ( ) , [ tildealm ] where @xmath79 is the rm field , and @xmath80 is the mask which is @xmath81 if a pixel is masked and @xmath50 if it is not .
we then evaluate c^rm_l = 1 f_sky(2l+1 ) _ m =-
l^l a^*_lm a_lm  , [ eq : clrm ] which are the so - called `` pseudo-@xmath82 '' rescaled by @xmath83 .    the commonly used procedure for estimating @xmath82 s from a partial sky is detailed in @xcite .
it assumes statistical isotropy and involves calculating a conversion matrix that relates `` pseudo-@xmath82 '' , obtained from a partial sky , to the `` full - sky '' @xmath82 . in eq .
( [ eq : clrm ] ) , as a crude approximation , we ignore the mode - coupling and replace the matrix with a rescaling by the @xmath83 factor . note that the galactic rm map is not isotropic and the `` full - sky '' spectrum reconstructed from a given patch is not the same as the actual spectrum evaluated from the full - sky .
thus , our estimates in sec .
[ sec : detection ] of the detectability of the galactic rm spectrum , and of the detectability of the pmf with the galactic rm as a foreground , are specific to particular masks .
the error bars in fig .
[ fig : rmcl ] receive contributions from the uncertainty in the measurements of the rm map and from the cosmic variance .
we are interested in answering two questions : 1)_how well can cmb polarization experiments measure the galactic rm spectrum ? _ 2)_what bounds on a scale - invariant pmf can be placed from measurements of fr ? _ both contributions , the galactic fr and the primordial fr , contribute in the same way to the quadratic estimators discussed in sec .
[ basics ] .
let @xmath84 be the estimator of the rotation angle that one reconstructs from linear combinations of quadratic terms @xmath85 .
the variance in @xmath84 , for a statistically isotropic rotation angle , is defined as @xmath86 with @xmath87 given by @xcite ^2_l = c_l^ + n_l^eb  , [ eq : variance ] where @xmath88 is the rotation power spectrum , and @xmath89 is the part of the variance that does not depend on the rotation angle and contains contributions from variances of individual estimators @xmath90 ( see eq .
( [ eq : dll ] ) ) , while accounting for their covariance .
it is given by @xcite ( n_l^eb)^-1 = 1 _  ( 2 + 1)(2+1)[c^ee_ h^l_]^2 _
^ee c_^bb  , [ eq : ebnoise ] where @xmath91even , and ^xx_^prim_+f_dl[c^xx]^wl_+i^xx _  , [ clvariance ] is the measured spectrum , with @xmath92 standing for either ee or bb , that includes the primordial contribution , @xmath93^{\rm prim}$ ] , the weak lensing contribution , @xmath93^{\rm wl}$ ] , and the noise term associated with the experiment , @xmath94 .
we allow for the possibility of partial subtraction of the weak lensing contribution by introducing a de - lensing fraction , @xmath95 , in ( [ clvariance ] ) .
the efficiency of de - lensing depends on the method used to de - lens @xcite , as well as the noise and the resolution parameters of the experiment . in @xcite
it was found that the quadratic estimator method of de - lensing @xcite , which is similar but orthogonal @xcite to the one we use to extract the fr rotation , can reduce the lensing contribution to bb by up to a factor of @xmath96 as it reaches the white noise limit . on the other hand , the so - called iterative method @xcite can further reduce the lensing contribution by a factor of @xmath60 or larger . in our analysis we have only considered two cases  one with no de - lensing ( @xmath97 ) and one with a factor of @xmath9 reduction of the lensing b - mode ( @xmath98 ) .
we adopt the latter optimistic assumption with the aim of illustrating the relative significance of lensing contamination for an experiment of given noise and resolution .
note that the numerator in eq .
( [ eq : ebnoise ] ) contains the unlensed e - mode spectrum , while the denominator includes the lensing contribution to e- and b - mode spectra .
predictions for both are readily obtained from camb @xcite for a given cosmological model .
the instrumental noise term , @xmath94 , accounts for the detector noise as well as the suppression of the spectrum on scales smaller than the width of the beam @xcite : i^ee_=i^bb_= ^2_p ( ^2 ^2_fwhm / 82 )  , [ eq : polnoise ] where @xmath99 quantifies detector noise associated with measurements of cmb polarization and @xmath100 is the full width at half maximum of the gaussian beam .    in the case of tb , the expression for the variance analogous to eq .
( [ eq : ebnoise ] ) is @xcite ( n_l^tb)^-1 = 1 _  ( 2 + 1)(2+1)[c^te_ h^l_]^2 _ ^tt c_^bb .
[ eq : tbnoise ]    eq .
( [ variance - alm ] ) assumes that the distribution of the rotation angle is isotropic , and the angular brackets denote an ensemble average over many realizations of @xmath0 as well as realizations of the primordial cmb sky .
this is a well - motivated assumption for the pmf and for the galactic rm near galactic poles .
away from the poles , the galactic rm shows a strong dependence on the galactic latitude @xcite .
this means that the forecasted variance can differ from the actual uncertainty for larger values of @xmath83 , although we expect it to be of the correct order of magnitude . for the galactic rm map , we take the rotation angular spectrum to be @xmath101 , with @xmath102 calculated from the rm map of @xcite after applying sky cuts , as plotted in fig .  [
fig : rmcl ] and explained in sec .
[ galaxy ] .
we will consider several sets of experimental parameters corresponding to ongoing , future and hypothetical experiments .
namely , we consider planck s @xmath51 ghz lfi and @xmath9 ghz hfi channels based on actual performance characteristics @xcite , polarbear @xcite at @xmath103 ghz with parameters compiled in @xcite , quiet phase ii @xcite at @xmath104 ghz using the parameters compiled in @xcite , @xmath51 , @xmath105 , @xmath58 and @xmath9 ghz channels of a proposed cmbpol satellite @xcite , as well as an optimistic hypothetical sub - orbital and space - based experiments at @xmath51 and @xmath103 ghz .
the assumed sky coverage ( @xmath83 ) , resolution ( @xmath100 ) , and instrumental noise ( @xmath99 ) parameters are listed in tables  [ table : galaxy ] and [ table : pmf ] .
for an optimistic @xmath51 ghz sub - orbital experiment with ( solid blue ) and without ( blue dot ) de - lensing , as well as for a hypothetical future @xmath51 ghz space probe with ( red short dash ) and without ( red long dash ) de - lensing ( @xmath98 ) .
the spiky nature of the plot is due to the spikiness of the rm spectra , as seen in fig .
[ fig : rmcl ] . experimental parameters assumed in this plot are given in the text and in table [ table : galaxy].,scaledwidth=45.0% ]    [ cols="^,^,^,^,^,^,^",options="header " , ]     -1 cm    the snr of the detection of the primordial fr rotation angle is given by ( s n ) ^2_xb = _ l=1^l_max f_sky ( 2l+1 ) [ c_l^,pmf]^2 2[c_l^,pmf + f_dgc_l^,g+n^xb_l]^2  .
[ eq : ebsnrp ] where @xmath106 stands for either e or t , depending on whether one uses the eb or tb estimator , @xmath107 is the primordial contribution , while @xmath108 is the `` de - galaxing '' factor , _
i.e. _ the fraction of the galactic rotation angle spectrum that is known from other sources and can be subtracted .    as in the case of the galactic rm ,
the snr receives a non - zero contribution from all multipoles @xmath30 , including the smallest , even for a small sky coverage . in the case
when a large fraction of the sky is available , there is an optimal cut of the sky , @xmath109 , for which the s / n is maximal .
this is because the galactic rm contribution is weaker when large parts of the sky are cut , but there are more available modes when more of the sky is kept . in
what follows , we will state when values of @xmath109 are different from @xmath83 .    in fig .
[ fig : dsnpmg ] , we plot contributions per @xmath110 to the net snr in detection of the fr from a scale - invariant pmf with @xmath111 ng for the same two hypothetical experiments considered in fig .
[ fig : dsngal ] .
in addition to the four cases considered in fig .  [
fig : dsngal ] , the green dot - dash line shows the case when the galactic contribution is partially subtracted for the space - based probe with @xmath112 .
a few observations can be made from this plot .
first is that both experiments can detect a @xmath8 ng pmf at high significance , which constitutes a big improvement on prospects of detecting a pmf in cmb via non - fr signatures @xcite .
second is that de - lensing makes a significant difference for both the space - based and sub - orbital experiments .
third is that subtracting the galactic rm moderately improves the detection at multipoles in the @xmath113 range .
finally , we can see that most of the information comes from @xmath114 , meaning that future cmb experiments can reconstruct fr maps corresponding to pmf of @xmath8 ng strength up to a degree resolution .    in table
[ table : pmf ] we present our forecasts for @xmath115 bounds on @xmath21 that can be expected from various ongoing and future cmb experiments .
we present results obtained from the eb estimator , because they lead to the strongest bounds , except for the case of planck s @xmath51 ghz channel , for which the tb estimator gives better constraints . for each experiment
, we check the effect of de - lensing ( + dl ) with @xmath98 as well as further partial subtraction of the galactic rm ( + dl+dg ) with @xmath112 .
the optimal sky cuts are indicated when relevant .
one of the interesting facts one can extract from this table is that very competitive bounds can be placed by polarbear and quiet , while even @xmath116 fields can be constrained in principle by future experiments .    to show the dependence of the pmf bound on the instrumental noise and resolution , we plot contours of constant @xmath115 bounds on @xmath21 in fig .
[ fig : contour ] .
the three panels show the cases with and without de - lensing and with an additional partial subtraction of the galactic rm .
an optimal sky cut was used for each set of parameters .
we have studied the detectability of a scale - invariant pmf due to fr of the cmb using polarization correlators in the presence of the foreground imposed by the milky way magnetic field .
we have found that the galactic rm is not a serious barrier to observe a scale - invariant pmf and that the eb quadratic correlator is in general more powerful than the bb corelator .
we have also studied the detectability of the pmf as a function of the sky coverage .
we find that suborbital experiments can be almost as effective as space borne experiments , and the obstruction caused by the galaxy is relatively weak if the observed patch is near the poles . also , as the mode - coupling correlations of cmb are mostly sourced by the largest scale features ( low @xmath30 ) of the rotation measure , a full sky cmb map is not necessary to access the pmf on the largest scales .
our results on the sensitivities of various experiments to the galactic magnetic field are summarized in table i. similarly in table ii we show the effective magnetic field strength that can be constrained by those experiments , taking into account the lensing of the cmb and the fr due to the milky way . the dependence of the observable pmf to the parameters of the experiment are shown in fig .  [
fig : contour ] .
here we have focussed on the observability of cmb fr at a single frequency .
cross - correlating polarization maps at multiple frequencies with comparable sensitivity to fr can further boost the significance of detection .
we will address this possibility in a future publication .
also , frequency dependence of polarization can be used to separate out fr induced b - modes from those induced by gravitational waves .
we thank the two anonymous referees for their careful reading of the initially submitted version and constructive suggestions that helped to improve the paper .
we are grateful to niels oppermann and collaborators @xcite for making their rotation measure maps publicly available at @xcite .
some of the results in this paper have been derived using the healpix @xcite package .
lp and tv benefited from previous collaborations with amit yadav .
we thank steven gratton , yin - zhe ma , phil mauskopf and vlad stolyarov for discussions and helpful comments .
sd is supported by a sese exploration postdoctoral fellowship ; lp is supported by an nserc discovery grant ; tv is supported by the doe at arizona state university .
d.  grasso and h.  r.  rubinstein , phys .
rept .   * 348 * , 163 ( 2001 ) [ astro - ph/0009061 ] .
l.  m.  widrow , rev .
phys .   * 74 * , 775 ( 2002 ) [ astro - ph/0207240 ] .
p.  a.  r.  ade _ et al .
_ [ planck collaboration ] , arxiv:1303.5076 [ astro-ph.co ] .
d.  paoletti , f.  finelli and f.  paci , mon .  not .
* 396 * ( 2009 ) 523 [ arxiv:0811.0230 [ astro - ph ] ] .
d.  paoletti and f.  finelli , arxiv:1208.2625 [ astro-ph.co ] . m.  s.  turner and l.  m.  widrow , phys .
d * 37 * , 2743 ( 1988 ) .
b.  ratra , astrophys .  j.   * 391 * , l1 ( 1992 ) .
t.  kahniashvili , y.  maravin , a.  natarajan , n.  battaglia and a.  g.  tevzadze , astrophys .
j.   * 770 * , 47 ( 2013 ) [ arxiv:1211.2769 [ astro-ph.co ] ] .
r.  a.  c.  croft , d.  h.  weinberg , m.  bolte , s.  burles , l.  hernquist , n.  katz , d.  kirkman and d.  tytler , astrophys .  j.   * 581 * , 20 ( 2002 ) [ astro - ph/0012324 ] .
t.  kahniashvili , y.  maravin and a.  kosowsky , phys .
d * 80 * , 023009 ( 2009 ) [ arxiv:0806.1876 [ astro - ph ] ] .
a.  yadav , l.  pogosian and t.  vachaspati , phys .
d * 86 * , 123009 ( 2012 ) [ arxiv:1207.3356 [ astro-ph.co ] ] .
t.  vachaspati , phys .
b * 265 * , 258 ( 1991 ) .
r.  durrer and c.  caprini , jcap * 0311 * , 010 ( 2003 ) [ astro - ph/0305059 ] .
k.  jedamzik and g.  sigl , phys .
d * 83 * , 103005 ( 2011 ) [ arxiv:1012.4794 [ astro-ph.co ] ] .
k.  jedamzik and t.  abel , arxiv:1108.2517 [ astro-ph.co ] .
n.  oppermann , h.  junklewitz , g.  robbers , m.  r.  bell , t.  a.  ensslin , a.  bonafede , r.  braun and j.  c.  brown _ et al .
_ , arxiv:1111.6186 [ astro-ph.ga ] .
m.  haverkorn , p.  katgert and a.  g.  de bruyn , astron .
astrophys .   * 403 * , 1045 ( 2003 ) [ astro - ph/0303644 ] .
l.  pogosian , a.  p.  s.  yadav , y.  -f .
ng and t.  vachaspati , phys .
d * 84 * , 043530 ( 2011 ) [ erratum - ibid .
d * 84 * , 089903 ( 2011 ) ] [ arxiv:1106.1438 [ astro-ph.co ] ] .
a.  kosowsky and a.  loeb , astrophys .
j.   * 469 * , 1 ( 1996 ) [ astro - ph/9601055 ] .
d.  d.  harari , j.  d.  hayward and m.  zaldarriaga , phys .
d * 55 * , 1841 ( 1997 ) [ astro - ph/9608098 ] .
m.  kamionkowski , phys .
lett .   * 102 * , 111302 ( 2009 ) [ arxiv:0810.1286 [ astro - ph ] ] .
s.  yadav , r.  biswas , 1 , m.  su and m.  zaldarriaga , phys .
d * 79 * , 123009 ( 2009 ) [ arxiv:0902.4466 [ astro-ph.co ] ] .
v.  gluscevic , m.  kamionkowski and a.  cooray , phys .
d * 80 * , 023510 ( 2009 ) [ arxiv:0905.1687 [ astro-ph.co ] ] . v.  gluscevic , d.  hanson , m.  kamionkowski and c.  m.  hirata , arxiv:1206.5546 [ astro-ph.co ] .
s.  yadav , m.  shimon and b.  g.  keating , phys .
d * 86 * , 083002 ( 2012 ) [ arxiv:1207.6640 [ astro-ph.co ] ] .
w.  hu and t.  okamoto , astrophys .
j.   * 574 * , 566 ( 2002 ) [ astro - ph/0111606 ] .
a.  r.  pullen and m.  kamionkowski , phys .
d * 76 * , 103529 ( 2007 ) [ arxiv:0709.1144 [ astro - ph ] ] .
monin a.  s. , iaglom a.  m. , statistical fluid mechanics : mechanics of turbulence . volume 2 , cambridge , mass . , mit press , ( 1975 ) .
k.  m.  gorski , e.  hivon , a.  j.  banday , b.  d.  wandelt , f.  k.  hansen , m.  reinecke and m.  bartelman , astrophys .
j.   * 622 * , 759 ( 2005 ) [ astro - ph/0409513 ] .
e.  hivon , k.  m.  gorski , c.  b.  netterfield , b.  p.  crill , s.  prunet and f.  hansen , astrophys .  j.   * 567 * , 2 ( 2002 ) [ astro - ph/0105302 ] . c.  m.  hirata and u.  seljak , phys .
d * 67 * , 043001 ( 2003 ) [ astro - ph/0209489 ] .
u.  seljak and c.  m.  hirata , phys .
d * 69 * , 043005 ( 2004 ) [ astro - ph/0310163 ] .
http://camb.info/    a.  lewis , a.  challinor and a.  lasenby , astrophys .
j.   * 538 * , 473 ( 2000 ) [ arxiv : astro - ph/9911177 ] . j.  r.  bond , g.  efstathiou and m.  tegmark , mon .  not .
astron .  soc .
* 291 * , l33 ( 1997 ) [ astro - ph/9702100 ] .
a.  r.  ade _ et al .
_ [ planck collaboration ] , arxiv:1303.5062 [ astro-ph.co ] .
b.  keating , s.  moyerman , d.  boettger , j.  edwards , g.  fuller , f.  matsuda , n.  miller and h.  paar _ et al .
_ , arxiv:1110.2101 [ astro-ph.co ] . y.  -z .
ma , w.  zhao and m.  l.  brown , jcap * 1010 * , 007 ( 2010 ) [ arxiv:1007.2396 [ astro-ph.co ] ] .
d.  samtleben [ quiet collaboration ] , nuovo cim .
b * 122 * , 1353 ( 2007 ) [ arxiv:0802.2657 [ astro - ph ] ] .
d.  baumann _ et al .
_ [ cmbpol study team collaboration ] , aip conf .
proc .   *
1141 * , 10 ( 2009 ) [ arxiv:0811.3919 [ astro - ph ] ] .


the problem of reconstructing binary matrices from their row and column sums is a classical inverse task in combinatorics ; see @xcite .
even though the term was introduced much later it can be seen as a root of _ discrete tomography _ ; see @xcite , and the collected editions @xcite , @xcite .
of particular relevance for the present paper are two well - known observations : ( 1 ) the question of consistency of the data i.e. , the question whether there exists a matrix whose row and column sums coincide with the given data , can be solved in polynomial time .
( 2 ) typically , its row and column sums do not determine the underlying matrix uniquely .
see @xcite for characterizations of the ` rare ' cases of uniqueness .    in the present paper
we address the second issue by adding additional window constraints , specifying ( or giving bounds on ) how many points there are in certain minors of the matrix .
these constraints come up naturally in _ dynamic discrete tomography _ @xcite .
an application of particular relevance in physics is that of particle tracking ; see @xcite , @xcite . here
, the positions of particles over time are to be reconstructed from two ( sometimes more ) high speed camera images .
window constraints are a natural way of modeling additional physical information for instance on the speed of the particles ; see @xcite and the background literature given there for further information .
here we are taking a complexity theoretical view commemorating observation ( 1 ) .
in fact , when adding various kinds of window constraints , we are interested in the boundary between polynomial - time solvability and @xmath1-hardness . intuitively speaking ,
we ask which kinds of additional constraints can be added without imposing a significant extra computational effort . or , phrased differently , what is the computational price to pay for reducing the number of solutions by utilizing additional window information .
we will focus on the effect of three different parameters : @xmath2 corresponds to the _ size _ of the window , @xmath3 is the _ number of @xmath4 s _ in the nonzero minors , and @xmath5 specifies the allowed positions of @xmath4 s in the windows , referred to as the _ pattern _ of the window .
the choice for these parameters will specify the given problem rec@xmath6 which is formally introduced in section  [ sec - results ] . hence @xmath7 @xmath8 and @xmath5
are given beforehand ( i.e. , are not part of the input ) .
as it will turn out , the problem exhibits various unexpected complexity jumps .    omitting technical details ( which are all given in section [ sec - results ] )
some of these jumps can be summarized as follows ; see table [ table : sumtable1 ] .
for @xmath9 the problems are in @xmath10 regardless on how the other parameters are set ; see theorem  [ thm : tract]([thm : tract : i ] ) .
for @xmath11 it depends on @xmath3 and @xmath5 whether the problems are in @xmath10 or @xmath1-hard ; see theorems  [ thm : tract]([thm : tract : ii]),([thm : tract : iii ] ) and  [ thm : np ] . for @xmath11 ,
some values of @xmath3 render the problems tractable while others make them @xmath1-hard ; see theorems  [ thm : tract]([thm : tract : ii ] ) and  [ thm : np]([thm : np : ii ] ) even if the patterns are not restricted at all . adding a pattern constraint may turn an otherwise @xmath1-hard problem into a polynomial time solvable problem ; see theorems  [ thm : np]([thm : np : ii ] ) and  [ thm : tract]([thm : tract : iii ] ) .
the reverse complexity jump , however , can also be observed ; see theorems  [ thm : tract]([thm : tract : ii ] ) and  [ thm : np]([thm : np : i ] ) .
the present paper is organized as follows : section [ sec - results ] will introduce our notation and state the main results . the proofs of our tractability results are given in section [ sec - tractability ] while the @xmath1-hardness results are proved in section [ sec - intractability ] .
section  [ sect : finalremarks ] contains some final remarks .
in particular , it provides a ( potentially also quite interesting ) extension of our main problem and collects related implications of our main results .
we begin with some standard notation .
let @xmath12 @xmath13 and @xmath14 denote the set of integers , natural numbers , and non - negative integers , respectively
. for @xmath15 set @xmath16 , @xmath17 @xmath18:=\{1,\dots , k\}$ ] , and @xmath18_0:=[k]\cup\{0\}.$ ] the support of a vector @xmath19 is defined as @xmath20 . with @xmath21
we denote the all - ones vector of the corresponding dimension .
the cardinality of a finite set @xmath22 is denoted by  @xmath23 we will use the notational convention that specific settings of variables and parameters are signified by a superscript @xmath24    in this paper we often refer to points ( or variables ) @xmath25 where @xmath26 are points of the integer grid  @xmath27 .
then , @xmath28 and @xmath29 denote the @xmath30- and @xmath31-coordinates of @xmath26 , respectively . in the grid @xmath32\times[n],$ ] the set @xmath33 $ ] and @xmath32\times\{j\}$ ] is called _
column  @xmath34 _ and _ row  @xmath35 _ respectively . for @xmath36
we refer to @xmath37 $ ] and @xmath32\times\{(j-1)k+1,\dots , jk\}$ ] as _ vertical strip  @xmath34 ( of width  @xmath2 ) _ and _ horizontal strip  @xmath38 ( of width  @xmath2 ) _ , respectively .
of course , the standard matrix notation can be obtained from our notation by the coordinate transformation @xmath39    generally , any subset of the @xmath32\times[n]$ ] grid is called a _
window_. windows of the form @xmath40\times[c , d])\cap\mathbb{z}^2,$ ] with @xmath41 and @xmath42 @xmath43 are called _ boxes . _ defining for any @xmath36 and @xmath44 the set of _ ( lower - left ) corner points _
@xmath45\times [ n])\cap(k\mathbb{n}_0 + 1)^2,$ ] we call any box @xmath46 ^ 2 $ ] with @xmath47 a _ block _ ; see figure  [ fig : ill ] for an illustration .
( here we depict the structure both as point set and as pixels . in the following
we restrict the figures to pixel images , which we find more intuitive . )
the blocks form a partition of @xmath32\times[n],$ ] i.e. , @xmath48\times[n].$ ] in the main part of this paper , we consider such non - overlapping blocks , which play also a role in super - resolution imaging @xcite . in section  [
sect : finalremarks ] we consider also other windows , which may be positioned at other places than those defined by @xmath49    next we introduce three patterns that we will study in detail since they exhibit already the general complexity jump behavior we are particular interested in .
the first pattern @xmath50 is unconstrained , i.e. , does not pose any additional restrictions on the positions of @xmath4 s . the second pattern @xmath51 forces all elements in the @xmath52 block to be @xmath53 except possibly for the two entries in the lower - left and upper - right corner .
the third pattern @xmath54 excludes all patterns that admit more than one  @xmath4 in each row of the @xmath52 block ( see also figure  [ fig : c1 ] ) . here
are the formal definitions
.    for @xmath55 let @xmath56_0 ^ 2}$ ] denote the power set of @xmath57_0 ^ 2.$ ] then we set @xmath58_0 ^ 2},\\ p(k,1)&:=\{\{(0,0)\},\{(k-1,k-1)\}\},\\ p(k,2)&:=\{m\in2^{[k-1]_0
^ 2 } \::\ : |m\cap\left([k-1]_0\times\{j\}\right)|\leq 1 \textnormal { for all } j\in [ k]_0.\}\end{aligned}\ ] ]    further , for @xmath59 and @xmath60,j\in[n]}$ ] we set @xmath61    for @xmath62 and @xmath63 we define now the following problem rec@xmath64    rec@xmath65 +    @xmath66 & \textnormal{task:}\quad & & \omit \rlap{find   $ \xi_{p , q}\in\{0,1\},$ \:$(p , q)\in[m]\times[n],$ with } \\
\displaystyle & & & \sum_{\mathclap{p\in[m ] } } \xi_{p , q}&&=r_q , & &    q\in[n ] , & & \textnormal{(row sums ) } \\
\displaystyle & & & \sum_{q\in[n ] } \xi_{p , q}&&=c_p , & &   p \in [ m ] , & & \textnormal{(column sums)}\\    \displaystyle & & & \sum_{\mathclap{(p , q)\in b_k(i , j ) } } \xi_{p , q } & & \leq v(i , j ) , & &   ( i , j)\in c(m , n , k ) , \quad & & \textnormal{(block constraints)}\\ \displaystyle & & & \textnormal{pat}_k(x , i , j)&&\in p(k , t),\qquad & &   ( i , j)\in c(m , n , k),\qquad & & \textnormal{(pattern constraints),}\\[1.2ex ] & & & \omit \rlap{or decide that no such solution exists.}\\\end{aligned}\ ] ]    in other words , we ask for @xmath67-solutions that satisfy given row and column sums , _ block constraints _ of the form @xmath68 with given @xmath69 and _ pattern constraints _ that restrict the potential locations of the 1 s in each block .
we remark that the @xmath70 @xmath71 can be viewed as some part of prior knowledge . in particle tracking
@xcite , for instance , the @xmath72 may reflect prior knowledge about physically meaningful particle trajectories ; see also @xcite .
our main results show that the computational complexity of @xmath73 may change drastically when  @xmath7  @xmath8 or  @xmath5 is varied .
[ thm : tract ]    a.   rec@xmath74 for any @xmath75 and @xmath76 [ thm : tract : i ] b.   rec@xmath77 for any @xmath78 [ thm : tract : ii ] c.   rec@xmath79 for any @xmath80 and @xmath81 [ thm : tract : iii ]    [ thm : np ]    a.   rec@xmath82-hard for any
@xmath78 [ thm : np
: i ] b.   rec@xmath83-hard for any @xmath78 [ thm : np : ii ]    the most notable changes are summarized in table  [ table : sumtable1 ] .
some of these changes may at first glance seem somewhat counterintuitive . for instance , restricting the solution space via pattern constraints turns the @xmath0-hard problem rec@xmath84 @xmath85 into the polynomial time solvable problem rec@xmath86 figure  [ fig : c1 ] depicts the possible types of blocks in rec@xmath87 conversely , additional pattern constraints convert the tractable problem rec@xmath88 into the @xmath0-hard problem rec@xmath89 @xmath90    ccc & @xmath91 & @xmath0-hard + varying @xmath2 & rec@xmath92 & rec@xmath84 + & rec@xmath93 & rec@xmath89 + varying @xmath3 & rec@xmath94 & rec@xmath84 + varying @xmath5 & rec@xmath95 & rec@xmath84 + & rec@xmath88 & rec@xmath89 +   +    -entries ( black pixels ) in each block of any solution of a @xmath96 instance.,scaledwidth=50.0% ]
this section contains the proofs of our tractability results stated in theorem  [ thm : tract ] .
clearly , a polynomial - time algorithm for rec@xmath97 can be given along the classical lines for deciding whether given row and column sums of a binary matrix are consistent . in the following we give one of the various arguments .
the problem obviously reduces to that of reconstructing a @xmath67-matrix from given row and column sums with some of the matrix entries fixed to @xmath98 it is well known that the coefficient matrix is totally unimodular , see , e.g. , @xcite , ( * ? ? ?
the problem can thus be solved efficiently , e.g. , by linear programming .
next we turn to the proof of theorem  [ thm : tract]([thm : tract : ii ] ) . for @xmath99
let @xmath100)\cap i|\qquad \textnormal{and}\qquad \rho_j(i):=|([i]\times\{j\})\cap i|\ ] ] denote the number of blocks  @xmath101 with @xmath102 @xmath103 and @xmath104 that lie in the same vertical and horizontal strip , respectively , as the block @xmath105 moreover , for @xmath106 we set @xmath107\times[n],\\ \pi_x(i)&:=\{i\in[m]:\exists j\in[n]:(i , j)\in i\},\\ \pi_y(i)&:=\{j \in [ n]:\exists i \in [ n]:(i , j)\in i\}.\end{aligned}\ ] ] note that @xmath108 and @xmath109 denotes the projection of  @xmath110 onto the first and second coordinate , respectively .
we use a ( slightly generalized version of a ) result of @xcite on the following problem .
dr@xmath111 +    @xmath112_0 , & & \textnormal{(row sum measurem . ) }   \\ & & & \omit \rlap{$\displaystyle c_{i+l}\in \mathbb{n}_0,$ } & & & & i\in\pi_x(i ) , \:\ : l\in [ k-1]_0 , & & \textnormal{(column sum measurem.)}\\[1.2ex ] & \textnormal{task:}\quad & & \omit \rlap{find   $ \xi_{p , q}\in\{0,1\ } , \:(p , q)\in g(i)$ with } \\
\displaystyle & & & \sum_{\mathclap{p:(p , j)\in g(i)}}\xi_{p , j+l}&&=r_{j+l } , & &   j\in\pi_y(i ) , \:\ : l\in [ k-1]_0,&&\textnormal{(row sums ) } \\
\displaystyle & & & \sum_{\mathclap{q:(i , q)\in g(i)}}\xi_{i+l , q}&&=c_{i+l } , & &   i\in\pi_x(i ) , \:\ : l\in [ k-1]_0,&&\textnormal{(column sums)}\\   \displaystyle & & & \sum_{\mathclap{(p , q)\in b_k(i , j)}}\xi_{p , q}&&=1,&\:\ : & ( i , j)\in i , & & \textnormal{(block constraints),}\\ & & & \omit \rlap{or decide that no such solution exists.}\end{aligned}\ ] ]    [ prop : mono ]    a.   [ mono2 ] an instance @xmath113 of dr@xmath111 is feasible if , and only if , for every @xmath114 we have @xmath115 b.   [ mono1 ] dr@xmath116    as a service to the reader we remark that a solution @xmath117 @xmath118 for a given instance of dr@xmath111 is obtained by setting for every @xmath114 and @xmath119 @xmath120 an illustration is given in figure  [ fig : example1 ] ( taken from @xcite ) .     with @xmath121 in gray color .
( right ) solution defined by  .,scaledwidth=45.0% ]    now we provide a polynomial - time algorithm for rec@xmath122    let @xmath123 denote an instance of rec@xmath122 we proceed in two steps .    first , we solve the following reconstruction - from - row - and - column - sums instance of finding @xmath124 @xmath71 satisfying the constraints @xmath125\cap(k\mathbb{n}_0 + 1),\\   \sum_{j:(i , j)\in i } \eta_{i , j}&=\sum_{l=0}^{k-1}c_{i+l } ,   & & i\in [ m]\cap(k\mathbb{n}_0 + 1),\\   \eta_{i , j}&=0 , & & ( i , j)\in \{(a , b)\in c(m , n , k):v(a , b)=0\}. \end{aligned}\ ] ] if no solution exists , we report infeasibility of @xmath126 the idea behind solving   first is that in this way we locate boxes that we want to contain a  1 .    in a second step , if is feasible , we invoke dr@xmath111 from proposition  [ prop : mono ] for the instance @xmath127 where @xmath128 and @xmath129 if there is no solution we report infeasibility of @xmath130 otherwise , we return the solution of dr@xmath131 it remains to be shown that this solves  @xmath113 correctly .    to this end , suppose that @xmath113 has a solution @xmath132 @xmath133\times[n].$ ] then ,   is also feasible ( a solution is given by @xmath134 @xmath47 ) and @xmath135 satisfies the corresponding instance of dr@xmath131 conversely , let @xmath136 @xmath71 denote a solution to   and let @xmath110 be defined as in  . for every @xmath114 we have clearly @xmath115 by proposition  [ prop : mono ]
there is thus a solution to dr@xmath131 this solution satisfies all row and column sum constraints of  @xmath113 by definition .
further since @xmath137 by definition of  , there are 1 s in the solution only in those blocks @xmath138 @xmath71 for which @xmath139 holds .
the block constraints are thus also satisfied , which proves the claim .
the integer linear program   can be solved in polynomial time as the coefficient matrix is well known to be totally unimodular ; see , e.g. , @xcite , ( * ? ? ?
the set @xmath110 in   is determined in polynomial time , and dr@xmath140 by proposition  [ prop : mono]([mono1 ] ) .
hence , rec@xmath141    next we turn to the proof of theorem  [ thm : tract]([thm : tract : iii ] ) .
it relies on the following result .
[ prop : totallyunimodular ] let @xmath142 be the node - edge incidence matrix of a ( simple ) bipartite graph , and let @xmath143 denote a binary @xmath144 matrix with the following properties .
a.   for every @xmath145 $ ] there exists an index @xmath146 $ ] with @xmath147 ; [ lem1:prop2 ] b.   the supports of any two distinct vectors @xmath148 and @xmath149 of @xmath150 with @xmath151 for some @xmath152 $ ] do not intersect .
[ lem1:prop3 ]    then , @xmath153 is totally unimodular .
it suffices to prove the assertion under the additional assumption that latexmath:[\[\label{eq : h }    all vectors @xmath149 whose support is a singleton can be added later since appending any subset of rows of the @xmath155-identity matrix to a totally unimodular matrix yields again a totally unimodular matrix .
as the matrix @xmath156 is the node - edge incidence matrix of a bipartite graph , let @xmath157 denote a corresponding partition of ( the indices of ) the vertices of the graph i.e. , a partition of @xmath158 $ ] .
clearly , for every @xmath159 with @xmath160 we have @xmath161    it suffices to prove that every collection of rows of @xmath162 can be split into two parts such that the sum of the rows in one part minus the sum of the rows in the other part is a vector with entries in @xmath163 ; see e.g. ( * ? ? ?
* theorem  19.3 ) . hence , let @xmath164 $ ] and @xmath165 $ ] denote such a collection of row indices of  @xmath156 and  @xmath166 respectively . for @xmath167
$ ] we set @xmath168 of course , the sets @xmath169 and @xmath170 form a partition of @xmath171 i.e. , @xmath172 also , it follows from  ( i ) , ( ii ) ,  ,  , and the fact that the graph is simple that @xmath173 i.e. , the sets constitute a partition of @xmath174 .    by  ( i ) ,  ( ii ) ,  ,   and   we have @xmath175 hence @xmath176 similarly , @xmath177 therefore , @xmath178 i.e. , the collection @xmath110 and @xmath174 of rows of @xmath156 and @xmath166 respectively , is split into two parts such that the sum of the rows in one part minus the sum of the rows in the other part is a vector with entries in @xmath163 .
now we present a polynomial - time algorithm for rec@xmath179 when @xmath180 .
we claim that the problem can be formulated as that of finding an integer solution to @xmath181 with a totally unimodular matrix  @xmath182 and integral right - hand side  @xmath183 showing that the problem is solvable in polynomial time ; see @xcite , ( * ? ? ?
16.2 ) .    assembling the variables
@xmath184 for @xmath185\times[n]$ ] into an @xmath186-dimensional vector and rephrasing the row and column sum constraints in matrix form @xmath187 where @xmath188 contains the  @xmath189 s and  @xmath190 contains the  @xmath191 s ,  @xmath156 is the node - edge incidence matrix of a bipartite graph , and hence totally unimodular .    as the pattern constraints ensure that no row of a block contains two 1 s , and since @xmath192 or @xmath193 we
can ( equivalently ) replace the block and pattern constraints by the _ box constraints _ @xmath194_0,\ ] ] where @xmath195\times\{j+l\}).$ ]    we can rephrase the box constraints in matrix form @xmath196 with a binary matrix  @xmath150 and a vector  @xmath197 containing the respective right - hand sides .    since each box @xmath198 is contained in a row of @xmath32\times [ n]$ ]
the condition in lemma  [ prop : totallyunimodular ] is satisfied .
also , any two of these boxes are disjoint since the blocks are disjoint .
hence the condition in lemma  [ prop : totallyunimodular ] is satisfied .
therefore , by lemma  [ prop : totallyunimodular ] , the matrix @xmath199 is also totally unimodular . appending @xmath200 and the identity matrix @xmath201 yields again a totally unimodular matrix ,
hence @xmath202 is totally unimodular .
since the corresponding right - hand side vector  @xmath203 is integral , the instance can be solved by linear programming .
this section contains the proof of theorem [ thm : np ] .
we begin with the @xmath0-hardness of @xmath204 for @xmath80 .
it suffices to show the result for @xmath205 since the @xmath0-hardness for larger  @xmath2 can be inferred from that for @xmath205 by setting the row and column sums @xmath206 and @xmath207 to zero for every @xmath47 and @xmath208 see @xcite .
so , let @xmath209 we use a transformation from the following problem 3-color tomography , which is @xmath0-hard by  @xcite .
( note that the third color is just `` blank . '' )    3-color tomography +    @xmath210.$}\\[1.2ex ] & \textnormal{task:}\quad & & \omit \rlap { find $ \xi_{p , q}^{(a)}\in\{0,1\},\:\:$ $ ( p , q)\in[m]\times[n],$ $ a \in [ 2],$ with } \\
\displaystyle & & & \sum_{p\in[m ] } \xi^{(a)}_{p , q}&&=r_q^{(a)},\qquad & &   q\in[n ] , a\in[2],&&\textnormal{(row sums ) } \\
\displaystyle & & & \sum_{q\in[n ] } \xi^{(a)}_{p , q}&&=c_p^{(a ) } , & &   p\in[m ] , a\in[2],&&\textnormal{(columns sums ) } \\
\displaystyle & & & \xi_{p , q}^{(1)}+\xi_{p , q}^{(2)}&&\leq1 , & &   ( p , q)\in [ m]\times[n]\qquad&&\textnormal{(disjointness condition),}\\[1.2ex ] & & & \omit \rlap{or decide that no such solution exists.}\end{aligned}\ ] ]    let @xmath211 denote an instance of 3-color tomography . we are setting up an instance  @xmath127 of our reconstruction problem by defining the row and column sums @xmath212 and @xmath213 , respectively .
the block and pattern constraints ensure that there are only three possibilities for setting the 1 s in each block .
they are shown in figure  [ fig:3colors ] .
instances.,scaledwidth=30.0% ]    note that the types can be viewed as representing three colors , which are counted by different row and column sums .
the row and column sums with odd indices count type  1 blocks , the other sums count type  2 blocks .    setting @xmath214,\:\ : & a\in[2 ] , \\
\displaystyle c_{2p+a-2}&:=c^{(a)}_p,\quad & & p\in[m],\:\ : & a\in[2],\end{aligned}\ ] ] and using the correspondences @xmath215,\\ \displaystyle { \xi_{p , q}^{(1)}}^*={\xi_{p , q}^{(2)}}^*=0 \quad & & \leftrightarrow \quad ( \xi_{2p-1,2q-1}^*,\xi_{2p-1,2q}^*,\xi_{2p,2q-1}^*,\xi_{2p,2q}^ * ) \textnormal { is of type~3},&&\end{aligned}\ ] ] we conclude that @xmath113 admits a solution if , and only if , @xmath127 admits a solution .
next we turn to the @xmath1-hardness of rec@xmath84 for @xmath216 .
again it suffices to prove the result for  @xmath209    the general structure of the proof follows that of the @xmath0-hardness result under data uncertainty given in @xcite . however , some adaptations and additions are required as will be detailed below .    as in  @xcite
we give a transformation from the @xmath0-hard problem 1-in-3-sat @xcite , which asks for a satisfying truth assignment that sets exactly one literal true in each clause of a given boolean formula in conjunctive normal form where all clauses contain three literals ( involving three different variables ) .    for a given instance of 1-in-3-sat a `` circuit board ''
is constructed that contains an initializer , several connectors , and clause chips .
the general structure of the circuit board from the proof of  @xcite and its modification , which contains additional connectors and initializers , are outlined in figure  [ fig : p2layoutleq ] .
the initializers contain , for every variable @xmath217 @xmath218,$ ] so - called _
@xmath219-chips _ , while the connectors contain _ @xmath220-chips . _
the clause chips are more complex as they consist of two _ collectors _ , two _ verifiers _ , and a _
transmitter_. the initializer holds a truth assignment for the variables @xmath221 of the given instance  @xmath113 of 1-in-3-sat .
the boolean values true and false are encoded in @xmath219-chips , @xmath218,$ ] by the type  1 and  2 blocks , respectively , that are shown in figure  [ fig : p2types ] .
-chips.,scaledwidth=17.0% ]    the proof in  @xcite shows how the truth assignment is transmitted through the circuit board , and how the verifier chips indeed check the feasibility of the given 1-in-3-sat instance .
figure  [ fig : p2example2 ] depicts a specific problem instance and a corresponding solution for the original construction of  @xcite .
the proof in  @xcite employs three different types of block constraints . for every block @xmath222 @xmath223 there is one of the following constraints : @xmath224 the names in parenthesis are those from @xcite ( even though  @xmath225 may seem more natural here in the third case ) .    as for binary variables ,
the constraint @xmath226 is equivalent to @xmath227 it suffices to adapt the construction in such a way that no @xmath228-block constraints are required .
in fact , the @xmath228-block constraints in the original construction are only employed for the @xmath220-chips , @xmath218,$ ] contained in the connectors .    for each @xmath220-chip , @xmath218,$ ] in a connector , we replace the @xmath228-block constraints by a @xmath229-block constraint . then , between any two clause chips we insert an additional copy of a connector and an initializer each .
these copies are placed as indicated in figure  [ fig : p2layoutleq](b ) .
the additional row and column sums are prescribed again to values  @xmath4 and  @xmath230 in alternation .
the new copies ensure that each @xmath220-chip , @xmath218,$ ] is contained in a horizontal or vertical strip that contains a @xmath219-chip of an initializer .
the corresponding row or column sums have values  @xmath4 and  @xmath231 and since we have @xmath229-block constraints for both chips , we can distribute the @xmath232 ones only in such a way that each of the two chips contains exactly two 1 s . in other words ,
each @xmath220-chip , @xmath218,$ ] is required to contain exactly two 1 s .
figure  [ fig : change ] shows the solution of our correspondingly adapted circuit board from  figure  [ fig : p2example2 ] .     a solution @xmath233 is shown ( non - zero components  @xmath234 are depicted as black pixels ) , representing the boolean solution @xmath235 , scaledwidth=60.0% ]    with this adaptation the proof from  @xcite carries over .
as a service to the reader we explain how the truth assignment is transmitted in our particular example ; see figure  [ fig : transmission ] .
the formal proof that , in particular , this transmission and , in general , the complete transformation works as needed , follows exactly as in  @xcite .
( for the example shown in  figure  [ fig : change ] ) .
, scaledwidth=60.0% ]    suppose the @xmath236-chip of the initializer in the lower right corner of the circuit board shown in figure  [ fig : transmission ] is of type  1 .
the truth assignment is vertically transmitted from the initializer to a connector , which negates the boolean value as the x - rays are set to  @xmath4 and  @xmath230 in alternation .
then , this boolean value is horizontally transmitted to an initializer ( again negated , so that we have the original values back ) , then vertically to a connector ( again negated ) . from here
the boolean value is transmitted horizontally to the first clause chip , where it enters and exits , yet again negated , through a collector . in the part between the two collectors , the boolean values for the literals appearing in the first clause
take a path through the verifiers that ensure that exactly one literal is set to true ; the values for the other literals are just transmitted ( again with negation ) .
this process continues in a similar way for the remaining clauses .
we end up with feasible solution for our instance of rec@xmath84 if , and only if , there is a satisfying truth assignment for our given instance of 1-in-3-sat .
let us first point out that the question of inferring information about an otherwise unknown object from observations `` under the microscope , '' i.e. , through certain windows , is fundamental for a wide range of applications . for corresponding results on the scanning of binary or integer matrices
; see e.g.  @xcite . in our terminology ,
@xcite studies the task of reconstructing binary matrices from their number of  @xmath4 s in each windows of a fixed size .
in addition to the patterns @xmath247 , and @xmath54 we consider , for @xmath216 , a fourth one @xmath248_0 ^ 2 } \::\ : |m\cap\left([k-1]_0\times\{j\}\right)|\geq k-1 \textnormal { for all } j\in [ k]_0\}.\ ] ] note that @xmath249 can be viewed as the colored - inverted version of @xmath250        @xmath254 & \textnormal{task:}\quad & & \omit \rlap{find   $ \xi_{p , q}\in\{0,1\},$ \:$(p , q)\in[m]\times[n],$ with } \\
\displaystyle & & & \sum_{\mathclap{p\in[m ] } } \xi_{p , q}&&=r_q , & &    q\in[n ] , & & \textnormal{(row sums ) } \\
\displaystyle & & & \sum_{q\in[n ] } \xi_{p , q}&&=c_p , & &   p \in [ m ] , & & \textnormal{(column sums)}\\    \displaystyle & & & \sum_{\mathclap{(p , q)\in w_{k}(i , j ) } } \xi_{p , q } & & \sim_{i , j } v(i , j ) , & &   ( i , j)\in c(m , n , k , l ) , \quad & & \textnormal{(window constraints)}\\ \displaystyle & & & \textnormal{pat}_{k}(x , i , j)&&\in p(k , t ) , \qquad \qquad \qquad & &   ( i , j)\in c(m , n , k , l),\qquad   & & \textnormal{(pattern constraints),}\\[1.2ex ] & & & \omit \rlap{or decide that no such solution exists.}\end{aligned}\ ] ]    tables  [ tablep01 ] and  [ tablenp01 ] below list tractability and intractability results for @xmath253 , respectively , which are simple corollaries to our main theorems [ thm : tract ] and [ thm : np ] .
in fact , they are either just reinterpretations of these results or rely , in addition , on one or two of the following simple transformation principles ( t1 ) , ( t2 ) , or ( t3 ) .
( t1 ) is the trivial observation that rather than reconstructing the ones in an image , we may reconstruct the zeros instead .
more precisely , this process of `` color inversion '' works as follows
. let @xmath255 denote the number of positions in the window @xmath256 .
then , we associate a given instance @xmath113 of @xmath253 with the instance @xmath127 of @xmath257 where @xmath258 @xmath259 and @xmath260 , @xmath261 @xmath262,$ ] @xmath263 @xmath264,$ ] @xmath265 and , for all @xmath266 , @xmath267 is replaced by @xmath268 clearly , the problems @xmath253 and @xmath257 lie in the same complexity class .
the transformation ( t2 ) adds empty rows and columns to an instance @xmath113 of @xmath269 to extend @xmath0-hardness results for @xmath205 to higher values of @xmath2 .
more precisely , let @xmath270 @xmath271 and set for @xmath272\times[n])\cap c(m , n,2)$ ] and @xmath273 $ ] @xmath274,\\ 0&:&\textnormal{otherwise,}\end{array}\right.\\ c'_{\frac{k}{2}(i-1)+l}&:=\left\{\begin{array}{lll}c_{i+l}&:&l\in[2],\\ 0&:&\textnormal{otherwise,}\end{array}\right.\\ ( \sim'_{i , j},v'(k(i-1)/2 + 1,k(j-1)/2 + 1)))&:=(\sim_{i , j},v(i , j)).\end{aligned}\ ] ] this defines an instance @xmath127 of @xmath275 with row sums @xmath276 column sums @xmath277 and block constraints @xmath278 @xmath279 clearly , the instance  @xmath113 admits a solution if , and only if ,  @xmath127 admits a solution ( by filling / extracting the  @xmath280-blocks of  @xmath113 into / from the  @xmath52-blocks of @xmath127 ) .
as this is a polynomial - time transformation , @xmath0-hardness of @xmath269 implies @xmath0-hardness of @xmath281    in the final transformation ( t3 ) additional rows and columns are filled by ones . for a given instance @xmath113 of @xmath282 and @xmath80 , let @xmath283 , @xmath284 , and set for @xmath272\times[n])\cap c(m , n,2)$ ] and @xmath273 $ ] @xmath274,\\ m'&:&\textnormal{otherwise,}\end{array}\right.\\ c'_{\frac{k}{2}(i-1)+l}&:=\left\{\begin{array}{lll}c_{i+l}&:&l\in[2],\\ n'&:&\textnormal{otherwise,}\end{array}\right.\\ ( \sim'_{i , j},v'(k(i-1)/2 + 1,k(j-1)/2 + 1)))&:=(\sim_{i , j},k^2+v(i , j)-4).\end{aligned}\ ] ]
this defines an instance @xmath127 of @xmath285 with row sums @xmath276 column sums @xmath277 block constraints @xmath278 @xmath286 @xmath287 @xmath288 and @xmath289 clearly , the instance  @xmath113 admits a solution if , and only if ,  @xmath127 admits a solution ( by filling / extracting the  @xmath280-blocks of  @xmath113 into / from the  @xmath52-blocks of @xmath127 ) . as this is a polynomial - time transformation , @xmath0-hardness of @xmath269 implies @xmath0-hardness of @xmath290    c@3exc@3exc@3exc@3exc@3exc@3exl@3ex@xmath2&@xmath291&@xmath292&@xmath293&@xmath5&@xmath174&reference + 1&@xmath14&@xmath294&@xmath294&@xmath295_0$]&@xmath296&thm .
[ thm : tract]([thm : tract : i ] ) + 2&@xmath297&@xmath298&@xmath299_0$]&0&@xmath300&@xcite + @xmath301&@xmath302&@xmath303&@xmath304&@xmath53&@xmath305&thm .
[ thm : tract]([thm : tract : ii ] ) + @xmath301&@xmath297&@xmath306&@xmath304&@xmath53&@xmath305&thm .
[ thm : tract]([thm : tract : ii ] ) & ( t1 ) + @xmath301&@xmath307&@xmath294&@xmath297&@xmath308&@xmath305&thm .
[ thm : tract]([thm : tract : iii ] ) + @xmath301&@xmath294&@xmath309_0$]&@xmath303&@xmath230&@xmath305&thm .
[ thm : tract]([thm : tract : iii ] ) & ( t1 ) +   +    c@4exc@4exc@4exc@4exc@4exc@4exl@4ex@xmath2&@xmath291&@xmath292&@xmath293&@xmath5&@xmath174&reference + @xmath301&@xmath294&@xmath297&@xmath310&0&@xmath296&@xcite [ tablenp01:14 ] + @xmath301&@xmath302&@xmath294&@xmath294&1&@xmath305&thm .
[ thm : np]([thm : np : i ] ) + @xmath301&@xmath311&@xmath294&@xmath294&0&@xmath305&thm .
[ thm : np]([thm : np : ii ] ) [ tablenp01:17 ] + @xmath301&@xmath312&@xmath294&@xmath294&0&@xmath305&thm .
[ thm : np]([thm : np : ii ] ) & ( t3 ) + @xmath301&@xmath294&@xmath313&@xmath294&0&@xmath305&thm .
[ thm : np]([thm : np : ii ] ) & ( t1 ) + @xmath301&@xmath294&@xmath314&@xmath294&0&@xmath305&thm .
[ thm : np]([thm : np : ii ] ) & ( t1 ) & ( t2 ) +   +    let us finally point out , that the tables  [ tablep01 ] and  [ tablenp01 ] can be largely extended . in particular , we can consider reconstruction problems for different types of windows which are even allowed to vary for each @xmath315\times[n]$ ] . while some of the complexity results in this more general setting may seem rather marginal
, other may be interesting for certain applications ; see @xcite . as this would go far beyond the scope of the present paper
we refrain , however , from carrying the results to the extremes here .
p.  gritzmann . on the reconstruction of finite lattice sets from their x - rays . in _ discrete geometry for computer imagery ( eds .
: e.  ahronovitz and c.  fiorio ) , dcgi97 , lecture notes on computer science 1347 , springer _ , pages 1932 , 1997 .
c.  kisielowski , p.  schwander , f.  h. baumann , m.  seibt , y.  kim , and a.  ourmazd .
an approach to quantitative high - resolution transmission electron microscopy of crystalline materials .
, 58(2):131155 , 1995 .
p.  schwander , c.  kisielowski , f.  h. baumann , y.  kim , and a.  ourmazd .
mapping projected potential , interfacial roughness , and composition in general crystalline solids by quantitative transmission electron microscopy .
, 71(25):41504153 , 1993 .    j.  zhu , j.  gao , a.  ehn , m.  aldn , z.  li , d.  moseev , y.  kusano , m.  salewski , a.  alpers , p.  gritzmann , and m.  schwenk .
measurements of 3d slip velocities and plasma column lengths of a gliding arc discharge . , 106(4):04410114 , 2015 .