\section{Introduction}

Jets in active galactic nuclei transfer the energy generated by the central supermassive black hole (SMBH)
to large( $>100$ kpc) distances.
The impact of jets on the environment contributes to the formation and evolution of structures in the
early Universe \citep{croton2006}. The innermost jets (parsec-scales or smaller) of radio-loud quasars are highly relativistic and their observed radiation can be Doppler amplified when observed at small angles to the line of sight.
These jets can be bright at high-energies and thus can provide interesting observational probes of the state of the SMBH activity \citep{begelman84}, but they remain spatially unresolved in X-rays and $\gamma$-rays.

Large scale X-ray jets span distances out to hundreds of kiloparsecs away from the SMBH
and encode the history of SMBH activity during the jet's lifetime (a few Myrs).
Their X-ray emission can be resolved with the \textit{Chandra}\ X-ray Observatory.
The number of such X-ray jets has significantly increased since the launch of \textit{Chandra}\ in 1999, but it is still relatively small in comparison to the number of known quasars.
There are about 100 large scale X-ray jets detected to date and only a few of them have good quality X-ray morphology data \citep{mas11}.
Though a direct connection between SMBH activity and the existence of kpc-scale jets is ambiguous, and the X-ray emission mechanism is not well understood, high-redshift ($z> 3$) jets could potentially establish the dominant energy environment in the early Universe.
Such jets probe the physics of the earliest (first $\sim$2~Gyr of the Universe in the quasars studied) actively accreting SMBH systems and are also interesting for other reasons.
For instance, the ambient medium in these high-redshift galaxies is probably higher than that of lower redshift galaxies \citep[e.g.,][]{dey06} and this may manifest itself in jets with different morphologies, with increased energy dissipation, or with the jets being slower in general than their lower-redshift counterparts.

The X-ray radiation could be attributed to either synchrotron emission
by highly relativistic electrons (Lorentz factors of $\gamma \sim 10^7-10^8$) in relatively strong magnetic fields, or
inverse Compton scattering of the cosmic microwave background (IC/CMB)
photons off the low energy ($\gamma \sim 10^3$) large-scale jet particles \cite[for a review see][]{har06}.  In
the simplest scenario, such models have diverging predictions at high
redshift.  Specifically, we expect a strong redshift dependence in the
X-ray--to--radio- energy flux ratio, ${x/r}\, = \frac{F_{\rm x}}{\nu_{\rm r} f_{\rm r}}$, where the radio fluxes are given at the observed frequency, and the X-ray fluxes are modeled over the energy range $0.5-7$~keV \citep{mas11}.
Typically, ${x/r}\, \propto U_{\rm CMB} \propto~(1+z)^{4}$ for IC/CMB, whereas in
synchrotron models, we do not expect a strong dependence,\footnote{The
energies of the synchrotron emitting electrons are different in the
observed radio and X-ray spectra, $\gamma_{\rm r} \sim 10^3$ vs. $\gamma_{\rm x}
\sim 10^7$ and these electrons may originate in the same population or
two different populations. Therefore, there could be some weak
redshift dependence in the synchrotron model.} ${x/r}\, \propto (1+z)^{0}$. Below we compare the predictions of these
two models for the highest-redshift relativistic jets.

Most \textit{Chandra}\ studies of quasar jets have so far targeted known arcsecond-scale radio jets \citep[e.g.,][]{sam04,mar11}, as most known examples are at $z$ $\stackrel{<}{{}_\sim}$2 \citep{bri84,liu02}.
At the time our program began, there were two high-$z$ quasars with kpc-scale X-ray jets: GB~1508$+$5714 at $z=4.3$ \citep{sie03,yua03,che04} and 1745$+$624 at $z=3.9$ \citep{che06}.
They were observed to have large ${x/r}\,$ values consistent with the IC/CMB model \citep{sch02,che04}, although the small number of high-$z$ detections precluded any definitive statements \citep{kat05,che06}.

We have therefore obtained \textit{Chandra}\ X-ray observations of an additional four high-redshift ($z>3.6$; GB\,1508$+$5714 was previously analyzed) and six intermediate-redshift ($2 \leq z \leq 3$) quasars with known radio jets.
The highest redshift X-ray and radio jet discovered in the sample studied (at $z=4.72$, in GB~1428+4214 (1428+422)) was presented and discussed in detail by \cite{che12}.
New and archival arcsecond-resolution NRAO\footnote{
The National Radio Astronomy Observatory is operated by Associated Universities, Inc.\ under a cooperative agreement with the National Science Foundation.}
imaging observations of these quasars are also presented.

The small number of X-ray photon counts observed from jets relative to their corresponding quasar cores means that detecting X-ray jets is inherently challenging.
Statistically, we must test the hypothesis that a baseline model of the quasar core and a flat background, without a jet, is insufficient to explain the observed data.
We do this test using a multi-scale Bayesian method known as Low Count Image Reconstruction and Analysis\footnote{
LIRA is implemented as a package for the R statistical programming language (r-project.org) that is available for downloading and use at {\tt github.com/astrostat/LIRA}.
} \citep[LIRA;][]{esch04, con07}.
The algorithm models the residual as a multi-scale component, and generates a series of images that capture the emission that may be present in excess of the baseline model.
We can then compute a $p$-value\footnote{
Formally, a $p$-value is the probability that the baseline null hypothesis can generate
 a value  for the test statistic
as large as that which is observed.
In this case, it defines the likelihood that a given
intensity can be obtained under the assumption that the baseline model is the truth.
That is, when the $p$-value is small, the chances that the feature under consideration can be attributed to a fluctuation is small.
This allows us to {\sl reject} the null hypothesis when this probability falls below a pre-defined threshold.
Note, however, that it should never be interpreted as a measure of the probability that the alternate hypothesis is true, nor, if the null cannot be rejected, as a measure of the probability that the null hypothesis is true \citep{wl16}.
}
by generating a series of Monte Carlo simulations of images under the baseline model and fitting each of these simulated images using LIRA.
\cite{stei15} (hereafter Paper~I) show how an upper-bound on the $p$-value can be computed with a small number of MCMC replicates.
We are interested in detecting whether X-ray jets exist in regions where jets were previously observed in the radio band.
In this paper, we will not consider X-ray detections without a corresponding radio emission (such a detection of a jet that was recently reported by \cite{simionescu2016}) when matching our results to the IC/CMB or synchrotron emission model.
We run LIRA to detect jets in pre-defined regions of an X-ray image.
Using the jets detected in X-rays we are able to observe how ${x/r}\,$ is dependent on redshift and whether it matches the predictions of the IC/CMB or synchrotron emission model.

Section 2 describes the sample selection and initial processing of the X-ray and corresponding radio observations.
Section 3 outlines how LIRA is used to find evidence that a jet exists in a region where one  is observed in radio imaging.
Section 4 elaborates on the results of the image analysis methods when applied to the new X-ray observations.
Section 5 gives a final description of our results in context and we summarize our results in Section 6.

\section{Sample Selection and Observations}

Compilations of known radio jets \citep{bri84,liu02} show very few examples with kpc-scale extensions at $z>3$ .
To increase the number of known radio jets, we carried out a VLA survey of the highest redshift (z$>$3.4) quasars and searched for extended radio emission \citep{che05,che08}.
The sample sources are the brightest (\simgt 100 mJy at 1.4 and/or 5\,GHz to facilitate ease of radio mapping), flat-spectrum radio sources catalogued by NED.
The flat spectra is a good proxy for high beaming and our high redshift sample is
representative of a general radio-loud quasars population.

Four flat-spectrum radio jets from the $z>3.4$ radio sample were observed in a \textit{Chandra}\ AO8 program (PI: Cheung) from 2007 Jan - June (ObsIDs: 7871-7874):
1239$+$376 ($z=3.82$),
1754$+$676 ($z=3.60$),
1418$-$064 ($z=3.69$),
1428$+$422 ($z=4.72$),
in addition to one quasar at $z=2.10$ (0833$+$585 ; ObsID 7870)
which had a known long radio jet \citep{mur93}.
The jets have minimum projected lengths of 2.5\hbox{$^{\prime\prime}$}\ and up to $\sim$15\hbox{$^{\prime\prime}$}\ for the $z=2.1$ case (Table~1), so their X-ray counterparts are easily separable from the bright nucleus in the \textit{Chandra}\ images.
The exposures of $\sim 3.8 - 11$ ks were tailored to the radio jet brightness and redshifts.

We also undertook a \textit{Chandra}\, survey of an intermediate redshift ($2 \leq z \leq 3$) sample of five sources in AO10 (PI: Sambruna) from  Jan - May 2009 (ObsIDs: 10307-10311) with exposures set to $\sim 20$ ks.
These radio jets were selected to have lengths greater than 2.5\hbox{$^{\prime\prime}$}\ from the list of \citet{liu02}.
In order to sample a range of jet orientations with respect to our line of sight, the selected targets have a range of radio core to extended flux ratios indicating both core-dominated and lobe-dominated radio sources.
A summary of the basic properties of the sources are recorded in Table~\ref{table-1}.

\subsection{Very Large Array Observations\label{sec:vla}}

The details of the VLA observations for the entire sample are shown in Table~\ref{table-2k}.
These radio jets form the largest sample of the radio jets at high redshift that were observed with {\it{Chandra}}.
We analyzed the archival data and selected the observations with the best astrometry for this project.
Portions of these radio data were presented in \citet{gob11} and \citet{gob14}.

The standard calibration was applied using AIPS \citep{bri94} with scans of primary calibrators, 3C~48 or 3C~286, used to set the flux density scales.
The data was then exported to the Caltech DIFMAP package \citep{she94} for self-calibration and imaging.

The radio jets in 1418$-$064, 1428$+$422, and 1754$+$676 were discovered in VLA 1.4 GHz A-configuration images from Oct - Dec 2004 (program AC755) and discussed in \citet{che05}.
We utilized the archival VLA snapshots of 1239$+$376 \citep{tay96} and 0833$+$585 \citep{mur93}.
The details of the 1428$+$422 observations are described in \citet{che12} but are included in Table~\ref{table-2k} for completeness.

The quasar 1418-064 has a hint of a radio extension in the VLA 1.4 GHz observation thus we obtained new A-array data at 5 GHz (July 27-28, 2007; program S8723) of this object which revealed more detail in the radio jet.

\subsection{\textit{Chandra}\ X-ray Observations.}

All X-ray observations were made with the \textit{Chandra}\ \citep{wei02, sch14} ACIS-S back-illuminated CCD.
The \textit{Chandra}\ observations are listed in Table~\ref{table-chandra}.
We used the nominal aimpoint of the ACIS-S3 chip and a 1/8th subarray mode (0.4s frame time) in order to mitigate pileup of the nucleus.
Roll angles were selected to place a possible charge transfer streak from the X-ray nucleus away from the jet position angles.
We reprocessed the data by using {\tt chandra\_repro} script in CIAO \citep{fru06} and assigned the
calibration available in CALDB v4.6.  We inspected the data for any
possible background flares and concluded that the data were not
affected by flares.  The background level in our short observations
was low and detections of the quasar emission and some jet features
are highly significant. We only use the X-ray events with
energies between 0.5-7\,keV in the spectral and image analysis
described below.

The X-ray nuclei are all clearly detected and the \textit{Chandra}\ positions (peaks) are within the 0.6\hbox{$^{\prime\prime}$}\ (90$\%$ pointing accuracy \citet{wei02}) of the radio positions (Table~\ref{table-1}).
We use the method described below for image analysis and detections of any features outside the core.

All spectral modeling was performed in Sherpa \citep{fre01} CIAO
version v4.6.  We extracted the spectra and created response files for
each observation using CIAO tools.
We used the Nelder-Mead optimization algorithm and Cash likelihood
statistics appropriate for low counts data and fit the spectra in
Sherpa.  We assumed an absorbed power law model for each identified
feature with the absorption column fixed at the appropriate Galactic
value (Colden; \cite{stark1992}).  The photon index of the power law and
the normalization were fit.  For the cores, we also fit a model which
included an additional intrinsic absorption component at the quasar
redshift.

\section{X-ray Image Analysis Methods}

\subsection{Initial Analysis}

For each \textit{Chandra}\ observation we generate an image centered on the quasar rebinned to a pixel size of $0.246$\hbox{$^{\prime\prime}$} (half of the native ACIS pixel size.)
For sources with a radio jet located at $< 7$\hbox{$^{\prime\prime}$}\, from the quasar core we use $64\times 64$ pixel ($15.7$\hbox{$^{\prime\prime}$}$\times 15.7$\hbox{$^{\prime\prime}$}) X-ray images.
For the three sources with larger radio structures we use $128 \times 128$ pixels ($31.5$\hbox{$^{\prime\prime}$}$ \times 31.5$\hbox{$^{\prime\prime}$}) images.

We employ SAOtrace\footnote{CXC Optics: http://cxc.harvard.edu/cal/Hrma/SAOTrace.html} to generate \textit{Chandra}\ point-spread functions (PSFs) for each observation, constructing the PSF with parameters appropriate to the quasar location.
Since X-ray PSFs are energy dependent, the raytraced PSFs include spectral models derived from fitting the spectrum of each quasar.
We use the same pixel size as in the X-ray image when binning the PSF.

We select regions of interest (ROI) for the quasar core and their jets using existing radio images, extracting spectra and response files using {\tt specextract}.
The radio images appear in panel (a) of Figure~\ref{fig:10307} and of Figures~\ref{fig:10308}--\ref{fig:2241} in Appendix B and the ROIs are overlaid on \textit{Chandra}\ ACIS-S images in panel (b) of each figure.
The circular quasar regions are marked `Q' and are centered on the radio core with the intent to encapsulate extended regions around the core, and the elliptical ROIs are associated with jet features and are numbered.
The complementary region is the remainder of the image that excludes both the quasar and jet regions and are marked 'C'.
The number of ROIs vary from one quasar to another.
Spectra are extracted from these regions and fit with an absorbed power-law model.
The results are in Table~\ref{tbl:jets}.
\begin{figure*}[ht]
\centering
\begin{tabular}[b]{@{}p{0.3\textwidth}@{}}
\centering\small (a)
\centering\includegraphics[scale=0.205,trim= 120 0 150 32,clip]{10307_radLab.png}\\
\end{tabular}
\quad
\begin{tabular}[b]{@{}p{0.3\textwidth}@{}}
\centering\small (b)
\centering \includegraphics[scale=0.205,trim= 14 0 10 23,clip]{10307_obs.png}\\
\end{tabular}
\quad
\begin{tabular}[b]{@{}p{0.3\textwidth}@{}}
\centering\small (c)
\centering\includegraphics[scale=0.2,trim= 120 0 170 20,clip]{10307_msLab.png}\\
\end{tabular}\\
\begin{tabular}[b]{@{}p{\textwidth}@{}}
\centering\small (d) \\
\centering\includegraphics[scale=0.5]{10307_fig3.png}
\end{tabular}
\caption{\small Images, regions of interest (ROIs), and results for quasar 0730$+$257 (ObsID 10307) -- (a) The VLA radio image with colors indicating the intensity in janksy per beam. The beam size is shown as an ellipse in the bottom right corner of the panel.-- (b) The observed \textit{Chandra}\ ACIS-S image with colors indicating photon counts. -- (c) The fitted added structure with colors indicating the residual intensity calculated from the posterior mean of LIRA ($\tau_1\Lambda_1$). -- (d.Q)-(d.C) The posterior distributions of $\xi$ in each ROI for the data (blue solid curve), the $50$ simulated replicate images under the null model, which includes the quasar and a background but not a jet (grey solid curves), and the average across them (black dashed curve). The ROIs are shown superimposed on panels (b) and (c). }
\label{fig:10307}
\end{figure*}
\FloatBarrier

\subsection{Image Reconstruction using LIRA}

Detection of an X-ray jet is particularly challenging in part because it is a faint structure located near a much stronger quasar core.
The Bayesian multi-scale fitting method LIRA is well suited to this challenge.
An efficient feature detection method is described in Paper~I as an addendum to the multi-scale fitting.
Section~2 of Paper~I gives full details of the statistical models used by LIRA, which we briefly review here.
We consider square images of $n=2^d{\times}2^d$ pixels, where $d$ is an integer (for typical cases considered here, $d=6$ or $7$).
We denote the counts in each pixel by $\mbox{\boldmath{$y$}} = (y_i, i=1 \ldots n)$, with the total number of photon counts, $N = \sum_i y_i$.
The image is then modeled as the sum of a known {\sl baseline component}, here representing the quasar core and background, and an {\sl added multi-scale component}, here representing the jet.
The counts observed in detector pixel $i$ are modeled as
\begin{equation}\label{eq:detected}
y_i \buildrel{\rm indep}\over\sim  \textrm{Poisson}\left(\sum_{j=1}^n P_{ji}A_j\left(\tau_0\Lambda_{0j} + \tau_1\Lambda_{1j}\right)\right) \,,
\end{equation}
where $(P_{ji}, i=1 \dots n, j=1 \dots n)$ is the PSF, representing the probability that a photon arriving from the direction represented by pixel $j$ is recorded in pixel $i$,
$(A_j, j=1 \dots n)$ is the exposure map, representing the efficiency with which incoming photons at pixel $j$ will be recorded,
$\tau_1$ is the total intensity of the added component,
and $\tau_0$ is the total intensity of the baseline component.
The quantities $\mbox{\boldmath{$\Lambda$}}_k = (\Lambda_{k1},\ldots,\Lambda_{kn})$ (for $k=0,1$) are the proportions of $\tau_k$ in the individual pixels.
Note that this represents a general model that is applied to the full image.
For ROIs smaller than the full image (see below), we first apply the model in Equation~\ref{eq:detected} to the full image, then recompute $\tau_k = \sum_{j{\in}{\rm ROI}} \mu_{kj}$, where $\mu_{kj}=(\tau_k\Lambda_{kj})$ are the inferred intensities in each pixel.

Paper~I describes how LIRA uses a flexible multi-scale model in order to accommodate added structure, such as a jet, in the quasar image.
LIRA combines the likelihood given in Equation~\ref{eq:detected} with prior distributions to formulate a posterior distribution for $\tau_0$, $\tau_1$, and $\Lambda_1$ and uses MCMC to provide a Monte Carlo sample from this posterior distribution.

\subsection{Using LIRA to Quantify Evidence for Extended Jet Emission}

LIRA uses a binned \textit{Chandra}\ image, a ray traced \textit{Chandra}\ PSF, and a composite Gaussian plus constant baseline model representing the quasar core and background respectively.
The MCMC sampler within LIRA returns a sequence of simulated images of the residual multi-scale component, including any non-Poisson deviations from the baseline model, which in the cases studied, would indicate a jet emanating from the X-ray quasar..
If the \textit{Chandra}\ image is consistent with the baseline model (quasar core + background), the fitted image of the multi-scale component contains only random fluctuations.
On the other hand, any X-ray emission not described by the baseline component appears as a systematic deviation from this noise.
Substantial deviations can be taken as evidence for structures beyond the baseline component in the \textit{Chandra}\ image.

Section~3
of Paper~I describes a formal statistical hypothesis test for the presence of structure (e.g., a quasar jet) in the added multi-scale component and how we can compute an upper bound on the $p$-value for this test.
This involves Monte Carlo evaluation of the behavior of the fitted multi-scale component (which describes the deviation from the baseline image) when there is no actual structure in the X-ray image beyond the baseline model (i.e., when there is no jet).
Quantitatively, this involves using a test statistic to efficiently compute an upper bound on the $p$-value.
Here we outline this method; further details can be found in Paper~I.

To start, for each observed image we simulate $50$ replicate images under a baseline-only model using Sherpa's {\tt fake} function.
The baseline model is composed of a 2D Gaussian for the quasar core on top of constant background and set to the best-fit values determined by fitting a 2D model, {\tt psf*gauss2d+const2d}, to the observed X-ray image.
For each quasar, we run LIRA \footnote{We only use LIRA iterations that are obtained after an initial burn in of the MCMC sampler.
In all cases, we obtain 2000 iterations for each run} and discard the first 1000
on the \textit{Chandra}\ images as well as each of its 50 replicate baseline-only simulated images.
The image of the deviation from the baseline when fit to the simulated images contain only random fluctuations, while those fit to the \textit{Chandra}\ images contain estimates of any structure beyond the baseline model.
A jet is considered detected if the signal in the region of interest is stronger than is likely to occur through random fluctuations under the baseline-only model.
Thus, we must compare the fitted deviations from the baseline of the \textit{Chandra}\ observation with the corresponding $50$ fitted deviations of the baseline-only simulations.

This is done quantitatively through a novel test statistic that compares the expected photon counts due to the baseline and added components, i.e., it compares $\tau_0$ with $\tau_1$ in Equation 1.
The test statistic is unusual in that it is a posterior probability, namely,
\begin{equation}
T_c(\mbox{\boldmath{$y$}}) = \Pr( \xi \geq c \mid \mbox{\boldmath{$y$}}),
\end{equation}
given a threshold $c$ where
\begin{equation}
\xi =\tau_1/(\tau_1 + \tau_0)
\end{equation}
is the proportion of expected counts due to the added structure.
If there is no structure in the added component (e.g., no quasar jet), $\xi$ and thus $T_c(\mbox{\boldmath{$y$}})$ are close to zero; as the intensity of the added component grows, $\xi$ and $T_c(\mbox{\boldmath{$y$}})$ tend toward one.
In this regard, $T_c(\mbox{\boldmath{$y$}})$ is a useful statistic for testing for added structure, i.e., for the quasar jet.

\section{Results}

\subsection{Quasars}

The quasar nuclei in our sample were observed by \textit{Chandra}\ with photon counts ranging from 64 to 5500.
We obtained a relatively good spectrum for each quasar core, and fit a power-law with a fixed Galactic absorption model to each spectrum using Sherpa \citep{fre01}.
We also considered an additional intrinsic
absorption component ({\tt
xsphabs(Gal)*xszphabs($z_{qso}$)*powlaw1d}).  The results of our
spectral fits are shown in Table~\ref{tbl:quasar2}.  An intrinsic
absorption is required in four quasars (0730+257, 1218+112, 1834+612, and 1428+422), with absorption columns of
${\approx}10^{22}$~cm$^{-2}$. For the remaining seven quasars we obtained upper limits. The photon indices for the power law
model range between $\Gamma_X=1.35$ and $\Gamma_X=2.26$
with an average of $\Gamma_X = 1.7\pm0.2$, consistent with estimates for large
samples of radio-quiet and radio-loud quasars \citep{kelly2007, young2009, lanzuisi2013}.  All the quasars
are X-ray luminous, with the luminosities in $2-10$\,keV energy band
exceeding $2.5\times 10^{45}$~erg~s$^{-1}$.
The most luminous quasar in our sample (1428$+$422) has a $2-10$\,keV luminosity of $4.0\pm 0.2 \times 10^{47}$~erg~s$^{-1}$.

\subsection{Jets}

We applied the method presented in Section 3 to our sample of 11 high-redshift quasars.
The posterior mean of $\tau_1 \Lambda_1$ is plotted in panel (c) of Figure~\ref{fig:10307} and of Figures \ref{fig:10308} - \ref{fig:2241} in Appendix B.
These represent the fitted multi-scale structure added to the baseline, i.e., the jets.
In order to assess the significance of the X-ray emission in each ROI and claim a detection, we calculate an upper bound on the $p$-value (with an upper tail probability of $\gamma = 0.005$, see section 3.4 of Paper~I).
This requires us to post-process the LIRA MCMC posterior sample to obtain a posterior sample of $\xi$ for each ROI.
Panel (d) of Figure~\ref{fig:10307} and of Figures~\ref{fig:10308}--\ref{fig:2241}  in Appendix B display the smoothed distributions of $\xi$ for each ROI for the data (solid blue line), and the averaged distributions of $\xi$ over the same ROI for the null simulations (dashed black).
The distributions of $\xi$ for each of the $50$ simulated null images are shown as grey lines.
We expect the observed distribution to depart from the distributions based on the null simulations if there is significantly more X-ray emission than expected under the baseline model.

We set a threshold of $p=0.01$ to determine the existence of a significant feature.
That is, if the upper bound on the $p$-value for the distribution of $\xi$ in an ROI is ${\leq}0.01$, we take that to be sufficient evidence for X-ray emission associated with a jet feature in the ROI.
Table~\ref{tbl:sigx} lists the upper bounds on the $p$-value for all the ROIs, showing that 16 features (12 jet features, 3 cores and 1 complementary region) have a $p$-values less than 0.01.
There are an additional 7 ROIs with an upper bound on the $p$-value between 0.01 and 0.02 which we consider marginal detections.
Thus, in total we list 23 significant X-ray features associated with high redshift jets including the corresponding cores and complementary regions.
In all these cases the distributions of $\xi$ based on the \textit{Chandra}\ data are skewed to the right compared to the distributions based on the null simulations.

There are three ROIs associated with the quasar cores (region Q for sources 1508$+$5714, 1428$+$422 and 1834$+$612) for which the upper bound on the $p$-value is lower than our pre-defined $p$-value threshold of $\alpha=0.01$.
This indicates the detection of the X-ray emission in excess of a point source assumed for the spatial model of the quasar core emission.
In general, all the images representing the deviation from the baseline model (panel (c) of Figure~\ref{fig:10307} and of Figures~\ref{fig:10308}--\ref{fig:2241} in Appendix B ) show indications of excess counts in the core ROIs in all the sources; we attribute these to an imperfect knowledge of the true telescope PSF, with the apparent deviations reflecting the uncertainties in the adopted PSF.
However, only in three quasars are they significant, signaling a departure from the baseline model. Note that this also shows that our threshold for detection is stringent enough to avoid detecting spurious features that may arise due to uncertainties in the PSF.

Note that the in the case of 1508$+$5714, a significant secondary structure in the complementary region `C' is seen, implying that not all the X-ray features were captured by the radio-based regions.
This object has a deep \textit{Chandra}\ observation of about 90~ksec, the longest one in our sample, and it is not unusual to find such structure in long observations.

Our X-ray jets detection rate -- 9 sources with detected jets out of 11 with known jets --
is similar to the one reported for the lower redshift sources ($\sim 70\%$ detection rates for exposures of $\sim$5--10 ksec; \cite{sam04,mar11}).

\section{Discussion}

\subsection{X-ray Morphology}

The posterior mean images $\tau_1 \Lambda_1$ representing deviations from the baseline model of the \textit{Chandra}\ data (panel (c) of Figure~\ref{fig:10307} and of Figures~\ref{fig:10308}--\ref{fig:2241} in Appendix B provide a view of the X-ray morphology with the quasar core removed.
However, not all the structures seen in these images are significant, since many features could be attributed to statistical fluctuations.
We have developed a method (Paper~I; see also Section 2 above) to assess the significance of the emission for well-defined ROIs.
Note that these ROIs {\sl must be set prior to the analysis} and cannot be deduced from the LIRA output, since doing so would increase the false detection rate.
We adopt regions based on the locations of the radio jet features, but note that the jet X-ray emission may not always be spatially coincidental with the radio emission \citep{sch02}, nor even have a radio counterpart \citep{jor04, simionescu2016}.
Offsets between the radio and X-ray peaks in the jet features have been reported \cite[e.g.][]{sie2007, worrall2009}.
For instance, in 1754$+$676 (ObsID 7872), the X-ray jet is not detected, but the image showing deviations from the baseline (see panel (c) of Figure~\ref{fig:7872}) suggests the existence of an emission feature between the quasar and the radio jet region.
A longer \textit{Chandra}\ observation is necessary to confirm this emission.
In another source 0805$+$046 (ObsID 10308), the image showing deviations from the baseline (panel (c) of Figure~\ref{fig:10308}) displays considerable emission outside the narrow radio jet, suggesting a more complex X-ray morphology.

We used the complementary regions to assess the possibility of unexpected X-ray emission present outside the pre-defined regions.
In all sources but 1508$+$5714, we do not find a strong indication that such emission is present, though the complementary regions cover a large area and thus statistical tests have relatively low power to detect smaller compact structures.
Future studies of the X-ray morphology in the vicinity of this source is required for understanding the origin of this emission.

The $p$-value upper bound test relies on the test regions being pre-defined.
This is done in order to avoid the loss of power in the test that arises when multiple hypotheses are tested.
We thus take the regions directly from the radio data and do not optimize the regions based on the X-ray data.
This could result in the size of location of the regions to be slightly misaligned, reducing the significance of detection.
For example, ROI 1 in source 1428$+$422 (ObsID 7874) is a marginal detection with an upper bound on the $p$-value of 0.010, but decreasing the region size from 91 to 77 pixels, an arguably better fit for this object in radio, results in an improvement in the $p$-value upper bound to 0.009, which crosses the threshold into a significant detection.
Areas of deviation from the baseline may be difficult to detect in a large region or a region may not encompass all of the relevant area of the image, and the better the region fits an area with deviation from the baseline, the lower the nominal $p$-value is.
But when large numbers of regions are tested, the $p$-value threshold must be reduced correspondingly in order to prevent false claims of detections due to fluctuations.
For instance, if the $\xi$ in 20 regions (say) of radii stepping from 75 to 95 pixels are tested, the appropriate threshold of $\alpha$ must be reduced by a factor of 20, to $\alpha=0.0005$, to maintain the same level of significance.
We emphasize that ROI selection must be consistent across the analysis, and {\sl must be defined before} applying the significance test.
We also note that we apply the test to a total of 47 ROIs, so we expect at most 1 false positive amongst the claimed significant and marginal detections at the significance threshold of $\alpha<0.02$.

As noted above, our current method of region selection depends on
radio data.  In the future, we plan to develop methods that are
independent of the radio (or other wavebands) selections and
autonomously generate regions that adaptively fit the deviations from
the baseline in the LIRA output.  Such a method is needed since the
X-ray emission does not always follow the radio closely.  Naturally,
any such method will trade-off ROI optimization for the statistical power of
the detection routine.

\subsection{Redshift Dependence in Large-Scale X-ray/Radio Emission}

The origin of X-ray jet emission is still under debate.
An early hint at the advantage of studying high-redshift jets came from the $z$=4.3 quasar 1508$+$5714 (Siemiginowska et al.\ 2003; Yuan et al.\ 2003). This quasar has higher X-ray to radio
luminosity ratio (${x/r}\,>$100) than any of its lower-$z$ counterparts (Cheung 2004).
This appears consistent with the $(1+z)^{4}$ amplification in the energy density of Cosmic Microwave Background (CMB):
\begin{equation}
{x/r}\, \propto u_{\rm CMB}/u_{\rm B} \propto (1+z)^{4}(\delta/B)^{2} \,,
\end{equation}
as expected under the IC/CMB model (e.g., Schwartz 2002).

Figure~\ref{fig:fluxRatio} compares the energy flux ratio ${x/r}\, = [F_{\rm x}~{(0.5-7\rm keV)}]/[\nu_{\rm r} f_{\rm r}]$ of the
detected and marginally detected jets across redshift from our
\textit{Chandra}\ sample.  We seek to establish whether or not the energy flux
ratio varies with redshift.  Figure~\ref{fig:density} shows the
posterior distribution from LIRA of the energy flux ratio for each
detected and marginally detected source (blue corresponds to low, and
red to high redshift).  It is visually apparent that there is a
difference in the distributions of the sources with higher redshifts
($z>3$) and those with lower redshifts.  In order to establish a
statistical measure of the significance of this difference, we split
the detected and marginally detected sources into two samples
consisting of the 18 at low redshift ($z<3$; sample $L$) and the 3 at
high redshift ($z>3$; sample $H$).  We then use a hierarchical Gaussian model to
examine whether low and high redshift quasars differ in terms of the
mean and variance of their energy flux ratio.  Appendix A describes a
procedure for evaluating the posterior probability that the difference
between the mean $\log_{10}$ energy flux ratio of the high and low
redshift jets ($\mu_H - \mu_L$) is greater than zero.
Figure~\ref{fig:mustat} shows the distribution of $\mu_H-\mu_L$
calculated from the posterior output.  We find an empirical
probability of $95\%$ that $\mu_H-\mu_L \geq 0$, which is at best
marginal evidence that the observed difference cannot be due to a
statistical fluctuation.  Though highly suggestive, because of the
small number of sources represented in this paper, and given the
disproportionate numbers of jets in the two samples, there is
insufficient evidence to conclude that the mean $\log_{10}$ energy
flux ratio differs between two groups of jets.  More observations at
$z>3$ are required in order to obtain more reliable results.
\begin{figure*}[ht]
\centering
\includegraphics[width=0.9\columnwidth]{fluxPlot_final.png}
\caption{
Ratio of X-ray to radio flux ${x/r}\,$ vs redshift for the detected and marginally detected regions of interest of the jets.
The circles are the energy flux ratios from jets detected in this study.
The diamonds are estimated energy ratios from the marginally detected jets.
The error bars form the $68\%$ interval from the LIRA iterations.
}
\label{fig:fluxRatio}
\end{figure*}

\begin{figure*}[ht]
\centering
\includegraphics[width=0.7\columnwidth]{fluxRatio_density.png}
\caption{
The posterior distribution of the ratio of X-ray to radio fluxes, ${x/r}\,$, for each detected and marginally detected jet.
The energy flux ratio is calculated at every iteration of LIRA.
The color corresponds to redshift.
}
\label{fig:density}
\end{figure*}

\begin{figure*}[ht]
\centering
\includegraphics[width=0.8\columnwidth]{muPlot.png}
\caption{
The difference in the mean X-ray to radio energy flux ratio between high- and medium-redshift quasars, $\mu_H - \mu_L$ at every LIRA iteration.
$\mu_L$ and $\mu_H$ are the average log flux ratio across the lower-redshift ($2<z<3$) and higher-redshift ($z>3$) detected and marginally detected redshift jets.
}
\label{fig:mustat}
\end{figure*}

The difference between the jets in two redshift groups is interesting
because it can also indicate that radio-loud quasars at $z>3$ are
different from low redshift ones. \cite{volonteri2011} hypothesized
that the jets at $z>3$ are systematically slower in comparison to the
jets at $z<3$.  If true, this could affect the energy flux ratio in
the framework of the IC/CMB model, due to a strong dependence of this
ratio on the jet Doppler factor. In this case, assuming that the
comoving jet magnetic field is roughly the same at different
redshifts, the observed increase in the energy flux ratio should be
smaller than that expected from the $(1+z)^4$ scaling. However, the
jet magnetization may evolve with redshift \citep[e.g.][]{singal2013}
which complicates the redshift scaling even further. Our results
indicate that the high redshift jets are different, but more
observations are needed to study the origin of this difference.

\subsection{Quasars at High Redshift}

The X-ray emission of radio loud quasars observed with {\it Chandra}
is unresolved and contained within $< 1.5\hbox{$^{\prime\prime}$}$ circular regions.
This emission could be due to a mixture of at least three components:
a hot corona directly related to the accretion process, a parsec scale
jet, and an unresolved portion of the kpc-scale outflow emitting
X-rays via IC/CMB.  We measured a standard range of photon indices for
the assumed power law model for quasar core spectra and found them
consistent with either process.  However, we detected relatively low
values of photon indices in a few quasars, including two at the
highest redshift (see Table~\ref{tbl:quasar2}).  Lower values of the
photon index are predicted if the jet dominates the X-ray emission.
In this case the beamed jet would make the quasars to appear more
luminous.  We notice that such a trend is present in our small sample
and five sources with $\Gamma_X < 1.6$ are more luminous, with an
average 2-10~keV luminosity of $6.3 \pm 2.5 \times
(10^{46}$erg~s$^{-1})$, than the other six with $\Gamma_X>1.6$ and the
average luminosity of $7.9 \pm 2.5 \times (10^{45}$erg~s$^{-1})$.
Formal correlation tests show that $\Gamma_X$ and ${\log}L_X$ are
indeed correlated, with Pearson's correlation coefficient $\rho=0.61$
($p=0.047$), and Kendall's $\tau=0.53$ ($p=0.024$). Such relation has
been seen in analysis of large samples of radio-loud and radio-quiet
quasars \citep{bechtold1994,
young2009,lanzuisi2013}. \cite{bechtold1994} found a similar trend in
a sample of radio loud quasars and argued that it could be caused by
an increased absorption. On the other hand \cite{young2009} did not
find a significant correlation in a radio-loud subsample of quasars
observed with XMM-Newton. Future studies of large number of radio loud
quasars in X-ray and radio band are necessary for understanding the
presence and origin of this correlation.

In the analysis of the {\it Chandra} images we assumed that the quasar
emission is point-like.  However, we detect signatures of
non-point--like emission in the two highest redshift ($z>4$) quasars
(1428+422 and 1508+5714) and in the one at $z=2.3$ (1834+612) (as evidenced by the fact that the
core component Q is not fully accounted for in the baseline model, see
Table~\ref{tbl:jets} and
Figures~\ref{fig:10311},\ref{fig:7874},\ref{fig:2241}).  This is
unlikely to be due to uncertainties in the shape of the \textit{Chandra}\ PSF,
since the residual core emission is not present in the multi-scale
components for other sources.  We suggest that these are due to
non-negligible contributions from unresolved kpc-scale jets emitting
IC/CMB.  The effect of this residual component on our analysis is
conservative, i.e., the imperfect modeling of the core tends to
increase the strength of the baseline model and systematically dampen
the added multi-scale component.
Thus, our results indicate that the unresolved X-ray cores of radio loud
high redshift quasars can contain significant contributions from
kpc-scale jet emission. This result is in agreement with studies based
on the optical-to-X-ray luminosity ratio \citep{saez2011,wu2013} and
along the expectations from the IC/CMB model.  Such a jet contribution
can potentially bias population studies of quasars and needs to be
taken into account in investigations of radio loud quasars at high
redshift.

\section{Summary}

We have analyzed \textit{Chandra}\ X-ray observations of a sample of eleven $z>2$ quasar radio jets including flat-spectrum and steep-spectrum radio sources.
The quasars were selected based on the known radio jets.
We use the Bayesian multi-scale image reconstruction method, LIRA, to obtain high-quality images outside the quasar core, and assess the significances of X-ray emission features coincident with jets.
We detect X-ray counterparts to radio jets in 9 quasars, including in the highest redshift X-ray and radio jet currently known (GB~1428+4217\ at $z=4.72$, Cheung et al.(2012)).
In particular, 12 radio features are detected at high significance, and an additional 6 at marginal significance.

We find that the ratio of X-ray to radio energy flux may differ
between jets at high ($z>3$) and low ($z<3$) redshift, in accordance
with the IC/CMB mechanism for X-ray emission, or pointing to the
intrinsic differences between low and high redshift quasars.  However,
this difference is subject to a large uncertainty due to small sample
sizes, and more observations are needed in the high redshift regime to
confirm this trend.

\acknowledgments

We remember Dan Harris' passion for understanding jets:  ``After all these years and all these conferences on jets, we still don?t know what jets are made of or how they work" (Harris 2015).
The authors acknowledge Rita Sambruna's contribution as the PI of the Chandra X-ray Observatory program for a sub-sample of quasars analyzed in this work. A.S. thanks Francesco Massaro and Giulia Migliori for discussion. Doug Gobeille participated in the initial selection of the radio jets.
This work was supported in part by the National Science Foundation REU and Department of Defense ASSURE programs under NSF Grant no.\ 1262851 and by the Smithsonian Institution.
The methodological aspects of this project were conducted under the auspices of the CHASC International Astrostatstics Center.
CHASC is supported by NSF grants DMS 1208791, DMS 1209232, DMS 1513546  and 1513492.
DvD acknowledges support from a Wolfson Research Merit Award (WM110023) provided by the British Royal Society and from Marie-Curie Career Integration (FP7-PEOPLE-2012-CIG-321865) and Marie-Skodowska-Curie RISE (H2020-MSCA-RISE-2015-691164) Grants both provided by the European Commission.
This research was supported in part by NASA through contract NAS8-03060 (A.S., D.A.S., V.L.K.) to the \textit{Chandra}\ X-ray Center. It was also supported by the \textit{Chandra}\ Award number GO7-8114 to Brandeis University by the \textit{Chandra}\ X-Ray Observatory Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of NASA under contract NAS8-03060.
L.S. was supported by Polish NSC grant DEC-2012/04/A/ST9/00083.

{\it Facilities:} \facility{CXO}, \facility{VLA}


\section{Introduction}
The Large Hadron Collider (LHC) at CERN, is nowadays at the start of a new program which is expected to run for about other 20 years. It is called High Luminosity LHC~\cite{HL-LHC} (HL-LHC) and it is CERN's number-one priority which aims to increase the number of collisions accumulated in the experiments by a factor of ten from 2024 onwards.

While HL-LHC project is already well defined for the next two decades, CERN has started an exploratory study for a future long-term project based on a new generation of circular colliders with a circumference of about 100 km. The Future Circular Collider (FCC) study~\cite{FCC}  has been undertaken to design a high energy proton-proton machine (FCC-hh), capable of reaching unprecedented energies in the region of 100 TeV, and a high-luminosity e+e- collider (FCC-ee), serving as $Z$, $W$, Higgs and top factory, with luminosities ranging from about $10^{34}$ to $10^{36}$ cm$^{-2}$s$^{-1}$ per collision point as a potential intermediate step towards the realization of the hadron facility. The design of the lepton collider complex will be based on the same infrastructure as the hadron collider.

At high beam intensity, necessary to reach the high luminosity foreseen for FCC-ee, the electromagnetic fields, self-generated by the beam interacting with its immediate surroundings and known as wake fields~\cite{palumbo}, act back on the beam, perturbing the external guiding fields and the beam dynamics. Under unfavorable conditions, the perturbation on the beam further enhances the wake fields; the beam-surroundings interaction then can lead to a reduction of the machine performances and, in some cases, also to instabilities.

The theory of collective beam instabilities induced by the wake fields is a broad subject and it has been assessed over many years by the work of several authors, such as F.~Sacherer~\cite{sacherer}, A.~W.~Chao~\cite{chao}, J.~L.~Laclare~\cite{laclare}, B.~Zotter~\cite{zotter}, C.~Pellegrini~\cite{pellegrini}, M.~Sands~\cite{sands} and others~\cite{yokoya1}.

To simplify the  study of collective effects, in general it is convenient to distinguish between short range wake fields, which influence the single bunch beam dynamics, and long range wake fields, where high quality factor resonant modes excited by a train of bunches can last for many turns exciting, under some conditions, coupled bunch instabilities. In both cases the bunch motion is considered as a sum of coherent oscillation modes perturbed by these wake fields.

In this paper we will focus on the FCC-ee collective effects induced by wake fields. In particular we will first evaluate the wake fields induced by the finite resistivity of the beam vacuum chamber (resistive wall) in section~\ref{s:sec1}. Due to the 100 km of length of the beam pipe, the resistive wall is going to play a fundamental role among the sources of wake fields for this accelerator, and the choice of the pipe geometry (circular or elliptical), material, and dimensions is particularly important.  In section~\ref{s:sec2} we study the collective effects induced by the resistive wall for both the short range (\ref{s:sub1}) and long range wake fields (\ref{s:sub2}), and for both longitudinal and transverse planes. For some instabilities we will resort to the linear theory, while for other cases and for more accurate predictions, we need to use simulation codes.

Section~\ref{s:sec3} is then dedicated to other important sources of wake fields, such as the RF system and tapers, and to the discussions on the choice of some devices, as the synchrotron radiation absorbers, in order to reduce their impact on the total wake field. Finally, concluding remarks and outlook will end the paper.

For reference we report in Table~\ref{tab1} the list of beam parameters we have used for evaluating the effects of wake fields on the beam dynamics.

\begin{table}[!ht]
\caption{Parameter list used to evaluate the beam dynamics effects of wake fields.}
\begin{center}
\begin{tabular}{|l|c|c|c|c|}
\hline
Circumference (km)&  100 & 100 &  100 & 100 \\
\hline
Beam energy (GeV)&  45.6 & 80 &  120 & 175 \\
\hline
Beam current (mA)&  1450 & 152 &  30 & 6.6 \\
\hline
RF frequency (MHz)&  400 & 400 &  400 & 400 \\
\hline
RF voltage (GV)&  0.2 & 0.8 &  3 & 10 \\
\hline
Mom.~compaction (10$^{-5}$) &  0.7 & 0.7 &  0.7 & 0.7 \\
\hline
Bunch length (mm)* &  1.6 & 2.0 &  2.0 & 2.1 \\
\hline
Energy spread (10$^{-3}$)* &  0.37 & 0.65 &  1.0 & 1.4 \\
\hline
Synchrotron tune &  0.025 & 0.037 &  0.056 & 0.075 \\
\hline
Bunches/beam &  90300 & 5260 &  780 & 81 \\
\hline
Bunch population (10$^{11}$)&  0.33 & 0.6 &  0.8 & 1.7 \\
\hline
Betatron tune &  350 & 350 &  350 & 350 \\
\hline
\end{tabular}
\\ * without beamstrahlung (no collision, worst case)
\end{center}
\label{tab1}
\end{table}%

\section{Resistive wall wake fields and impedances}
\label{s:sec1}

The electromagnetic interaction of the beam with the surrounding vacuum chamber, due to its finite resistivity, produces unavoidable wake fields, which, for FCC-ee, result to be of particular importance. If we consider a beam pipe with circular cross section and a single material of infinite thickness, the longitudinal monopolar ($m=0$) coupling impedance is given by~\cite{bane}
\begin{equation}
\label{rw1}
\frac{Z_{||} \left(\omega\right)}{C}=\frac{Z_0 c}{\pi} \frac{1}{\left[1+i\textrm{sgn}\left(\omega\right)\right] 2bc
\sqrt{\frac{\sigma_c Z_0 c}{2 \left| \omega \right|}}-i b^2 \omega}
\end{equation}
and the transverse dipolar ($m=1$) one by
\begin{equation}
\label{rw2}
\frac{Z_{\perp} \left(\omega\right)}{C}=\frac{Z_0 c^2}{\pi} \frac{2}{\left[\textrm{sgn}\left(\omega\right)+i\right] b^3c
\sqrt{2\sigma_c Z_0 c\left| \omega \right|}-i b^4 \omega^2}
\end{equation}
where $C$ is the machine circumference, $Z_0$ the vacuum impedance, $c$ the speed of light, $b$ the pipe radius, and $\sigma_c$ the material conductivity. The above expressions are valid in a frequency range defined by
\begin{equation}
\frac{\chi c} {b} \ll \omega \ll \frac{c \chi^{-1/3}}{b}
\end{equation}
with $\chi=1/\left(Z_0 \sigma_c c b \right)$. The corresponding wake functions are given by~\cite{bane1}
\begin{equation}
\label{eq:long_wake}
\frac{w_{||}\left( z \right)}{C} =\frac{4Z_0 c}{\pi b^2} \left[ \frac{e^{-z/s_0}}{3}
\cos \left( \frac{\sqrt{3}z}{s_0} \right)
-\frac{\sqrt{2}}{\pi} \int_{0}^{\infty} dx \frac{x^2 e^{-z x^2/s_0}}{x^6+8}
\right]
\end{equation}
and
\begin{equation}
w_{\perp}\left( z \right)=\frac{2} {b^2} \frac{dw_{||}\left( z \right)}{dz}
\end{equation}
with $z>0$ and $s_0=\left[2b^2/\left(Z_0 \sigma_c \right) \right]^{1/3}$.

By considering a beam pipe of 35 mm of radius made by copper (conductivity of about 5.9e7 S/m) or aluminium (conductivity of about 3.8e7 S/m), eqs.~(\ref{rw1}) and (\ref{rw2}) are valid in a very large range of frequency.
In addition, it is important to observe that the last term in the denominator of eqs.~(\ref{rw1}) and (\ref{rw2}) is negligible up to high frequencies, giving then the possibility to easily evaluate the scale of the impedance with the pipe radius. Indeed the longitudinal impedance is inversely proportional to the beam pipe radius, and the transverse one to the inverse of the third power of $b$. This scaling can be used to find a compromise for the pipe geometry. By reducing the radius it is possible to reduce the power required for the magnets, but this would increase in particular the coupling impedance and then reduce transverse instability thresholds.

In the following we will use, as a value for the pipe radius, 35 mm, which, as it will be shown in the next section, gives instability thresholds which are quite safe with respect to the nominal beam parameters\footnote{Except for the transverse coupled bunch instability for which, as it will be discussed in section~\ref{s:sub2}, an active feedback system has to be foreseen.}, and represents a good compromise also for the quadrupole magnets' power. Indeed, by using the impedance scaling with radius it is easy to infer scenarios at different pipe radii.

The above equations represent the impedance per unit of length. Due to the fact that the length of the accelerator in this case is of unprecedented value, also the total resistive wall coupling impedance is very high. However it is important to note that the effects of the impedance on the beam dynamics do not scale directly with the machine length. As we will see in the next section, other parameters peculiar of this machine, as the energy and the tunes, are important. In addition we also observe that the above impedances depend on the frequency where they are evaluated. If, for example, we consider the transverse impedance at the revolution frequency, we see that it approximately scales as $C\sqrt{C}$, and not simply with $C$. Some discussions about the scales will be done when dealing with the different instabilities in the next section, when we will see that the mere contribution of the resistive wall to the total wake fields gives important effects in both the longitudinal and transverse planes.

Concerning the choice of the material, we see that, since both the impedances are proportional to the square root of the resistivity, the difference between copper and aluminium gives a factor of about 1.26, which is not considered a critical value.

The above equations, however, are valid only for a single thick layer, for which the skin depth is much smaller than the wall thickness. In Fig.~\ref{skindepth} we show the skin depth for copper and aluminium as a function of frequency. We can see from the plot that at the revolution frequency of about 3 kHz (red line in the figure), the skin depths are about 1.2 and 1.5 mm. For lower frequencies the skin depth is higher. At the moment of writing this paper, the first layer of the vacuum chamber is aluminium with a thickness of 4 mm~\cite{kersevan1}, that means that eqs.~(\ref{rw1}) and (\ref{rw2}) can be used up to very low frequencies, and only at about 100 Hz the difference between one and multiple layers starts to give some differences.

\begin{figure}[!ht]
\center
\includegraphics[width=115mm]{skindepth.pdf}
  \caption{Skin depth as a function of frequency for aluminium (blue) and copper (green). The red line represents the revolution frequency.}
\label{skindepth}
\end{figure}

In Fig.~\ref{multilayer} we show the transverse impedance $Z_{\perp} \left( \omega \right)/C$ as a function of frequency for a circular pipe with radius 35 mm, aluminium, in case of single layer, as given by eq.~(\ref{rw2}), and three layers, with a first layer of aluminium of 4 mm, then 6 mm of dielectric and finally iron with resistivity of $10^{-7}$ $\Omega$m. In the three layers case, the impedance has been evaluated with the code ImpedanceWake2D~\cite{mounet}. As shown in the figure the difference between the two impedances starts to show up at very low frequency. As a conclusion we can say that all the considerations derived from eqs.~(\ref{rw1}) and (\ref{rw2}) are essentially valid also for the multilayer case.

\begin{figure}[!ht]
\center
\includegraphics[width=115mm]{multilayer.pdf}
  \caption{Transverse impedance as a function of frequency for single layer (real part in blue and imaginary part in green) and three layers (real part in red and imaginary part in cyan).}
\label{multilayer}
\end{figure}

For what concerns the pipe geometry, the options are between circular and elliptic cross section. An elliptic chamber with semi-axes 60x35 mm was initially considered for allowing to insert synchrotron radiation absorbers on one side of the vacuum chamber along the major semi-axis. However, even if the elliptic shape gives a reduced impedance in the horizontal plane, it produces an additional quadrupolar transverse impedance which has to be taken into account. In this case the Yokoya form factors~\cite{yokoya2}, shown in Fig.~\ref{yokoya}, give a reduction of the transverse impedance of about 15\%, but an additional quadrupolar impedance which is about 36\% the dipolar one.
\begin{figure}[!ht]
\center
\includegraphics[width=115mm]{yokoya.pdf}
  \caption{Yokoya form factors. $a$ and $b$ are the major and minor semi-axes respectively.}
\label{yokoya}
\end{figure}

Fig.~\ref{RW} shows the total transverse and longitudinal resistive wall impedance as a function of frequency for a circular beam pipe of 35 mm of radius in the three layer case. This impedance is used in the following section for evaluating the resistive wall effect on beam dynamics.
\begin{figure}[!ht]
\center
\includegraphics[width=160mm]{RW_impedance.pdf}
  \caption{Real and imaginary part of transverse (left) and longitudinal (right) impedance of resistive wall as a function of frequency.}
\label{RW}
\end{figure}

\section{Effects of resistive wall on beam dynamics}
\label{s:sec2}

\subsection{Single bunch effects}
\label{s:sub1}

One important effect of the resistive wall on the single bunch dynamics is related to the transverse mode coupling instability, or strong head tail instability~\cite{chao}. In general the transverse motion of the bunch, derived from the Vlasov equation, can be decomposed as a sum of coherent modes of oscillation, called eigenmodes, the coherent frequencies of which depend on the current intensity and on the coupling impedance. If the bunch distribution function can be expressed as a sum of orthogonal polynomials, it is possible to evaluate amplitude and frequencies of the eigenmodes, that could couple together giving rise to an instability. The frequencies of the coherent modes are here calculated with DELPHI~\cite{mounet} code, which considers Laguerre polynomials. In Fig.~\ref{TMCI} we show the real part of the frequency (tune shift) of the first two radial coherent oscillation modes, with the azimuthal number going from -2 to 2, as a function of the bunch population for 45.6 GeV and 80 GeV. As expected, the worst scenario is at the lowest energy, where we find an instability threshold that is a factor of about 6 higher than the nominal bunch population. The higher energy cases, not shown here, give even higher thresholds. In this situation we can see that, if other contributions to the transverse impedance do not exceed the resistive wall, we have a good margin of safety for this kind of instability. However, a more detailed study of transverse mode coupling instability with a more detailed transverse impedance is necessary.
\begin{figure}[!ht]
\center
\includegraphics[width=80mm]{TMCI.pdf}
\includegraphics[width=80mm]{TMCI2.pdf}
   \caption{Real part of the frequency of the first coherent oscillation modes as a function of bunch population for the 45.6 GeV case (left) and the 80 GeV case (right).}
\label{TMCI}
\end{figure}

For what concerns the longitudinal beam dynamics, one main problem caused by the resistive wall is related to the longitudinal potential well distortion and the evaluation of the microwave instability threshold. The microwave instability does not produce a bunch loss, but the consequent longitudinal emittance increase and possible bunch internal oscillations that cannot be counteracted by a feedback system, make the microwave instability an effect that has to be studied with care. In addition to that, there are no reliable analytical expressions that can be used to easily evaluate the instability threshold. For these reasons we have performed a series of simulations by using a tracking code, which we refer here as SBSC~\cite{music}, initially developed to study the longitudinal beam dynamics in the electron storage ring DA$\Phi$NE~\cite{sbunch1}, and successively developed and adapted to other machines~\cite{sbunch2}.

As other single bunch tracking codes taking into account the wake fields effects, in order to reduce the computing time, the code uses macro-particles, divides the bunch into slices, or bins, and evaluates the wake potential, that is the convolution integral of the wake function times the longitudinal distribution function, at the centre of each slice to then interpolate the wake for all the positions of the macro-particles. In Fig.~\ref{long_wakefield} we have represented the wake function given by eq.~(\ref{eq:long_wake}). As we can see, the wake varies considerably in a very short distance. This causes problems to simulations with slices. If we consider, for example, an rms bunch length $\sigma_z$ of about 3-4 mm, and use $\pm 5 \sigma_z$ for the beam dynamics, in order to have about 20 slices in 50 $\mu$m, we need to have several thousands of slices, making the simulation very time consuming.
\begin{figure}[!ht]
\center
\includegraphics[width=120mm]{long_wakefield.pdf}
   \caption{Resistive wall longitudinal wake function used for simulations.}
\label{long_wakefield}
\end{figure}

In Fig.~\ref{long_wakepot_bins} in blue we show the wake potential of a 1.6 mm Gaussian bunch obtained with the tracking code using $2\times 10^4$ slices compared with the analytical calculation (black line) given by the equation~\cite{piwinski}
\begin{equation}
\label{eq:rw}
W_{||}(z)=\int_{-\infty}^{\infty}\lambda (z') w_{||}(z-z')dz'=\frac{cC}{8\sqrt{2}\pi b \sigma_z^{3/2}} \sqrt{\frac{Z_0}{\sigma_c}}F(z/\sigma_c)
\end{equation}
with
\begin{equation}
\label{bessel}
F(x)=|x|^{3/2} e^{-x^2/4}\left(I_{1/4}-I_{-3/4}\pm I_{-1/4} \mp I_{3/4} \right)
\end{equation}
where $I_n$ is the modified Bessel function, the upper signs in eq.~\ref{bessel} are for positive $z$, and $\lambda(z)$ is the longitudinal distribution function.

The oscillations of the wake from slice to slice given by the code are due to the fact that, in order to have a good statistics, in addition to the excessive number of slices, we should also use a very high number of macro-particles.
\begin{figure}[!ht]
\center
\includegraphics[width=120mm]{long_wakepot_bins.pdf}
   \caption{Comparison of the analytical longitudinal wake potential of a 1.6 mm Gaussian bunch (black line) with those given by the simulation code, one using $2\times 10^4$ slices and the wake function as Green function (blue line), and the second one with 400 slices and the wake potential of a 0.15 mm Gaussian bunch as Green function (red line).}
\label{long_wakepot_bins}
\end{figure}

It is possible to overcome this problem by using not the wake function but the wake potential of a very short bunch as Green function for the code. The same figure, with the red line, shows the wake potential given by the same simulation code using only 400 slices and, as Green function, the wake potential of a 0.15 mm Gaussian bunch. As we see, the wake is much smoother due to the fact that the same number of macro-particles are now distributed over 400 slices instead on $2\times 10^4$, and the agreement with eq.~(\ref{eq:rw}) is excellent.

In order to test the code, and to evaluate the effect of the resistive wall on the longitudinal beam dynamics, we have first solved the Ha\"\i ssinski integral equation~\cite{haissinski}, which is able to predict the bunch length and the distortion from a Gaussian distribution for intensities below the microwave instability threshold. The equation can be written as
\begin{equation}
\lambda(z)=\lambda_0 \exp{ \left[ \frac{1} {E_0 \eta \sigma_{\varepsilon 0}^2} \Psi(z) \right]}
\end{equation}
with $\lambda_0$ a normalization constant, $E_0$ the collider energy, $\eta$ the slippage factor, $\sigma_{\varepsilon 0}$ the natural RMS energy spread, and
\begin{equation}
\Psi(z)=\frac{1}{C} \int_{0}^{z} \left[ e V_{RF}(z') - U_0\right] dz'-\frac{e^2 N_p} {C} \int_{0}^z dz' \int_{-\infty}^{z'} \lambda(z'') w_{||}(z'-z'') dz'' \ .
\end{equation}
where $V_{RF}$ represents the total RF voltage, $U_0$ the energy lost per turn due to the synchrotron radiation, and $N_p$ the bunch population.

The bunch shapes for different bunch populations at the lowest energy of 45.6 GeV are shown in Fig.~\ref{haissinski}.
\begin{figure}[!ht]
\center
\includegraphics[width=120mm]{haissinski.pdf}
   \caption{Longitudinal distribution for different bunch population as given by Ha\"\i ssinski equation.}
\label{haissinski}
\end{figure}

The bunch length is about 2.4 - 2.5 mm at the nominal current, but we have to remind that only the resistive wall effect has been taken into account for the moment. For the three shown bunch populations the tracking code gives exactly the same distribution.

In Table~\ref{tab2} we have also reported the bunch length calculated by using the Ha\"\i ssinski equation with the nominal bunch population of Table~\ref{tab1} for all the four energies and due only to the resistive wall effect. We can see that the lowest energy case gives the largest effect, whilst for 120 and 175 GeV the lengthening due to the resistive wall wakefield is almost negligible.

\begin{table}[!ht]
\caption{Bunch length given by Ha\"\i ssinski equation with resistive wall impedance and with the parameters of Table~\ref{tab1}.}
\begin{center}
\begin{tabular}{|l|c|c|c|c|}
\hline
Beam energy (GeV)&  45.6 & 80 &  120 & 175 \\
\hline
Bunch length (mm)&  2.48 & 2.37 &  2.15 & 2.22 \\
\hline
\end{tabular}
\end{center}
\label{tab2}
\end{table}%

The potential well distortion theory described by the Ha\"\i ssinski equation predicts a bunch length increasing with current and a constant energy spread up to a given threshold, called microwave instability threshold, above which also the energy spread increases. In the microwave instability regime, even if the bunch is not lost, it could be characterized by internal turbulent motion which would compromise the machine performances. Several papers have been written to determine the microwave instability threshold~\cite{yunhai}.  In particular, in ref.~\cite{stupakov0}, the microwave instability due to the resistive wall wake fields was analyzed giving a criterion for the threshold evaluation. Applied to the FCC-ee case, it gives a threshold value of $N_p=8.1 \times 10^{10}$, a factor slightly higher than 2 with respect to the nominal bunch population.

This value can be compared with the results of the tracking code. From Fig.~\ref{energyspread}, where we represented the RMS energy spread given by the code as a function of the bunch population, we can see that the energy spread starts to increase at about $8-10\times 10^{10}$. This is in a good agreement with the above analytical estimate.
\begin{figure}[!ht]
\center
\includegraphics[width=120mm]{en_spread.pdf}
   \caption{RMS energy spread as a function of bunch population given by the simulation code with only the resistive wall impedance.}
\label{energyspread}
\end{figure}

As a further check of the tracking code results, a Vlasov-Fokker-Planck solver~\cite{warnok} has also been used, showing that up to a bunch population of $8\times 10^{10}$ the beam is stable and giving the onset of the instability at about $10-12\times 10^{10}$.

Finally, Fig.~\ref{bunchlength} shows the RMS bunch length, obtained with the simulation code, as a function of the bunch population up to an intensity of $2\times 10^{11}$, a value six times higher than the nominal bunch population.
\begin{figure}[!ht]
\center
\includegraphics[width=120mm]{bunch_length.pdf}
   \caption{RMS bunch length as a function of bunch population as given by the simulation code with only the resistive wall impedance.}
\label{bunchlength}
\end{figure}

\subsection{Multi-bunch effects}
\label{s:sub2}

A more critical situation is related to the transverse coupled bunch instability due to the long range transverse wake fields. In this case the study can be performed by considering the motion of the entire beam (not of the single bunch) as a sum of coherent oscillation modes, with coupled bunch modes to be taken into account. By considering the lowest azimuthal mode $m=0$ and a Gaussian bunch, the real part of the coupling impedance can produce stability or instability depending on the sign of the growth rate
\begin{equation}
\label{tcbi}
\alpha_{\mu,\perp}=-\frac{c I} {4 \pi (E_0/e) Q_{\beta}} \sum_{q=-\infty}^{\infty}Re \left[ Z_{\perp} \left( \omega_q \right) \right]
G_{\perp} \left(\sigma_{\tau} \omega'_q \right)
\end{equation}
where $I$ the total beam current, $Q_{\beta}$ the betatron tune, $\sigma_{\tau}$ the RMS bunch length in time, $G_{\perp}$ a form factor which, for our case, is about 1, and
\begin{equation}
\label{eq:modes}
\omega_q = \omega_0 \left(q N_b + \mu + Q_{\beta} \right) \qquad \omega'_q = \omega_q+\omega_0\xi \frac{Q_{\beta}}{\eta}
\end{equation}
with $N_b$ the number of bunches, $\xi$ the chromaticity, and $\omega_0$ the revolution frequency.

In the above equations, $\mu$ represents the $\mu^{th}$ coupled bunch mode, which goes from 0 to $N_b-1$. When $\alpha_{\mu}$ is positive, the corresponding mode is unstable. If we consider, as transverse impedance, the resistive wall one given by eq.~(\ref{rw2}), and ignore the term $-ib^4 \omega^2$, we observe that $Re \left[ Z_{\perp} \left( \omega \right) \right]$ depends on the sign of the frequency $\omega$. Negative frequencies produce unstable modes with an exponential growth given by eq.~(\ref{tcbi}), while positive ones give rise to damped oscillations. In addition to that, the resistive wall impedance grows approximately with the inverse of the square root of the frequency, determining the most dangerous coupled bunch mode when $\omega_q$ is as close to zero as possible. If we consider, as an example, the parameters given by table~\ref{tab1}, and $q=-1$, by denoting with $Q_0$ the integer part of the betatron tune, that is $Q_{\beta} = Q_0+ \nu_{\beta}$, with $\nu_{\beta}$ the fractional part of the tune, which plays a crucial role for this kind of instability, it comes out that the most dangerous coupled bunch mode is $\mu = N_b - Q_{0}-1=89949$, and this mode has its lowest negative frequency at $\omega_q=-\omega_0\left( 1-\nu_{\beta} \right)$.

Fig.~\ref{tbci1} shows the beam spectrum of three coupled bunch modes and the real part of the resistive wall impedance of a circular pipe of aluminium, with radius of 35 mm and three layers, close to zero frequency for two extreme cases of fractional part of the betatron tune, $\nu_{\beta}=0.05$ (left) and $\nu_{\beta}=0.95$ (right), and we see that a smaller fractional tune is preferred to alleviate the transverse coupled bunch instability because the impedance has a lower value. Due to dynamic aperture and  beam-beam issues, and since FCC-ee has 2 interaction points, the fractional tunes are indeed just above the integer~\cite{oide}, and therefore its fractional part is close to zero, mitigating the instability growth rate.
\begin{figure}[!ht]
\center
\includegraphics[width=165mm]{tcbi1.pdf}
  \caption{Coupled bunch spectrum and real part of the resistive wall impedance as a function of frequency around $f=0$ for fractional tune $\nu_{\beta}=0.05$ (left) and $\nu_{\beta}=0.95$ (right).}
\label{tbci1}
\end{figure}

If we consider, as an approximation, not a sum of the impedance over frequency in eq.~(\ref{tcbi}), but the coupling with a single betatron frequency line of the coupled bunch modes, the most dangerous unstable mode has a growth rate given approximately by
\begin{equation}
\alpha_{\perp}=\frac{c I} {4 \pi (E/e) Q_{\beta}} \frac{C}{2\pi b^3}\sqrt{\frac{CZ_0}{\pi |1-\nu_{\beta}|\sigma_c}}
\end{equation}
which, for the best case with $\nu_{\beta}=0.05$, gives a growth rate of about 432.4 s$^{-1}$, corresponding to a rise time of approximately 2.3 ms, that is about 7 machine turns. If the fractional tune increases, the situation worsens. In Fig.~\ref{tbci2} we show the growth rate as a function of the fractional tune for the worst coupled bunch mode (left plot), and as a function of the coupled bunch mode number for the case with $\nu_{\beta}=0.05$ (right plot).
\begin{figure}[!ht]
\center
\includegraphics[width=165mm]{tcbi2.pdf}
  \caption{Growth rate as a function of the fractional tune for the most dangerous coupled bunch mode (left), and as a function of the coupled bunch mode number for the case with $\nu_{\beta}=0.05$ (right).}
\label{tbci2}
\end{figure}

A more precise calculation by considering the sum in eq.~(\ref{tcbi}) and by using the Laguerre polynomials with the DELPHI code confirms the values of the growth rates. As we can see from the right plot of Fig.~\ref{tbci2}, there are several unstable modes, and all of them need to be damped. In addition to that, even if the rise times are in the range of few milliseconds, which are not of particular concern for an accelerator machine, due to the large circumference, the rise times correspond to very few turns, making very challenging the realization of a feedback system. Some schemes that could deal with this problem have been proposed in ref.~\cite{drago}.

In addition to the previous instabilities, for the elliptic vacuum chambers, as we have seen in the previous section, there are  quadrupolar resistive wall wake fields that can produce a substantial tune variation with beam current in the multi-bunch regime. These are particularly important for the lowest energy FCC-ee option (45.6 GeV) with the highest beam current (1450 mA) distributed over a large number of circulating bunches.

For a rectangular vacuum chamber inside dipole magnets, the tune slope is estimated by~\cite{chao2}
\begin{equation}
\frac{d Q_{\beta}}{dI} =\pm \left( \frac{\pi r}{48 Q_{\beta}} \right) \left(\frac{Z_0}{E_0/e} \right)
\left( \frac{R}{b}\right)^2 \left(\frac{L}{C} \right), \qquad r=1+\frac{b^2}{d^2}
\end{equation}
where $R$ is the radius of the machine, $b$ and $d$ the half height and half width of the vacuum chamber, $L/C$ the ring fraction with the rectangular chamber (dipoles).

This equation has reproduced very well the tune shifts observed in the high current lepton factories PEP-II~\cite{chao2} and DA$\Phi$NE~\cite{mikhail}. For the sake of estimates we consider the rectangular vacuum chamber with a cross section of 2bx2d= 70 mm x 120 mm and assume that the dipole vacuum chambers occupy $L/C=0.59$ of the ring, similarly to the HER of PEP-II. In these conditions, for the energy of 45.6 GeV we obtain the tune variation of
\begin{equation}
\frac{d Q_{\beta}}{dI} =\pm  \frac{0.253}{A}
\end{equation}
As a consequence, during the beam injection and up to the nominal value of 1.45 A, the tune variation would be as large as about 0.37. In order to cope with such big tune changes additional feedback systems would be necessary to keep the tune constant as required by the beam-beam interaction, dynamic aperture considerations and other beam dynamics issues.

For what concerns possible longitudinal coupled bunch instabilities excited by HOMs, at this stage it is not possible to quantify their impedance contribution, but we can estimate the maximum shunt impedance that gives a growth rate that can be compensated by the natural radiation damping.

Similarly to the transverse case, starting from the linearized Vlasov equation, and developing any perturbation of the beam distribution as sum of coherent oscillation modes, it is possible to obtain an eigenvalue system representing the coherent frequencies of the modes. By neglecting the coupling between the different azimuthal modes, and by considering only the lowest longitudinal azimuthal mode $m = 1$, it is possible to show that the real part of the HOM impedance can produce stability or instability depending on the sign of the growth rate
\begin{equation}
\alpha_{\mu,||}=\frac{\eta I}{4 \pi \left(E_0/e \right) Q_s} \sum_{q=-\infty}^{\infty} \omega_q Re \left[ Z_{\||} \left( \omega_q \right) \right] G_{||} \left(\sigma_{\tau} \omega_q \right)
\end{equation}
with $Q_s$ the synchrotron tune and $\omega_q = \omega_0(qN_b+\mu+Q_s)$. Stability in this case occurs for negative frequencies because the real part of the longitudinal impedance is always positive, and the worst and simplest unstable case is when the HOM has its resonant angular frequency $\omega_r$ equal to $\omega_q>0$. Similarly to the transverse case, if we consider, as an approximation, not a sum of the impedance over frequency, but the coupling with a single synchrotron frequency line of the coupled bunch modes, the most dangerous unstable mode has a growth rate given approximately by
\begin{equation}
\label{eq:growth_long}
\alpha_{||}=\frac{\eta I}{4 \pi \left(E_0/e \right) Q_s} \omega_r R_s
\end{equation}
with $R_s$ the HOM shunt impedance. Also in this case $G_{||}(x) \simeq 1$, if $f_r \ll 25$ GHz. This growth rate has to be compared with the natural damping rate due to the synchrotron radiation, which, for the lowest energy machine, is about 1320 turns. In Fig.~\ref{HOMs}, we have represented the maximum HOM shunt impedance of eq.~(\ref{eq:growth_long}) as a function of the resonance frequency, such that the corresponding growth rate is exactly balanced by the radiation damping. Of course, similarly to the transverse case, also here a feedback system has to be developed as a further safety knob.

\begin{figure}[!ht]
\center
\includegraphics[width=115mm]{HOMs.pdf}
 \caption{Maximum shunt impedance of a HOM as a function of its resonance frequency, producing a growth rate that is compensated by the natural radiation damping.}
\label{HOMs}
\end{figure}

\section{Other important impedance sources}
\label{s:sec3}

In the previous section, by discussing the effects of the resistive wall, we have seen that its contribution on the beam dynamics is very important, requiring, in some cases, active feedback systems to keep under control beam instabilities.
Due to the above results, we can consider, as a rough indication to accept the impedance of a given device in the design stage, the comparison of its wake potential with the resistive wall one.

Generally, for a high energy e+e- machine, the handling of synchrotron radiation represents one of the most challenging tasks. As a consequence of this, it is necessary to have a sufficient number of cavities to recover the lost energy, and, in addition to that, synchrotron radiation absorbers, necessary to cope with extra heating and eventual background. These devices  could be important sources of impedance, and their evaluation is considered in this section and compared with the resistive wall contribution.

Let us first estimate the impact of the synchrotron radiation absorbers. For FCC-ee a synchrotron radiation absorber will be installed every 4-6 meters, with the purpose of intercepting the radiation that, otherwise, would impact on the beam chamber. Due to their large number, the absorbers represent a very important source of machine impedance.  Initially, an elliptical vacuum pipe, with absorbers on one side of the major semi-axis, as shown in Fig.~\ref{old_absorber}, was proposed. A total of 9228 absorbers was estimated~\cite{kersevan}.
\begin{figure}[!ht]
\center
\includegraphics[width=115mm]{old_absorbers.pdf}
  \caption{Elliptic vacuum chamber with one absorber~\cite{kersevan}.}
\label{old_absorber}
\end{figure}

Of course, the contribution of a single absorber could be neglected, but the effect of all of them is expected to be prohibitively large compared to the resistive wall of a circular beam pipe. In Fig.~\ref{wp_old_abs} we have represented the wake potential of a 4 mm Gaussian bunch given by all the absorbers as obtained by CST Particle studio~\cite{CST} and by using a simpler and more approximated approach that takes into account a rectangular geometry with two absorbers at each side~\cite{stupakov}. In the same figure the wake potential of 10000 absorbers given by CST considering a new alternative design described below is also shown.
\begin{figure}[!ht]
\center
\includegraphics[width=115mm]{wp_old_abs.pdf}
  \caption{Longitudinal wake potential of a 4 mm Gaussian bunch produced by: resistive wall (red line), 9228 absorbers (CST simulations, blue line), 18456 absorbers (two absorbers at each side) with a more approximated approach taking into account a rectangular geometry~\cite{stupakov} (green line), 10000 absorbers for a new alternative design (CST simulation, magenta line).}
\label{wp_old_abs}
\end{figure}

Even if there is a certain degree of uncertainty in the results, and clearly the more approximated approach may overestimate the longitudinal wake potential using two absorbers for symmetry reasons, both approaches demonstrate that the impedance contribution of the absorbers in this condition is unacceptable. In addition to this result, if we further consider the previous comments on the tune shifts produced by the quadrupolar wake fields of the elliptic vacuum chamber, it is clear that the elliptic cross section is not the best solution for FCC-ee from the impedance point of view.

A proposed alternative is the installation of a circular vacuum chamber, with radius 35~mm and winglets on both sides, as the one of Super-KEKB~\cite{SuperKEKB}. Bellows with RF fingers and valves are to be designed with the same profile. Transitions to a circular cross section chamber (without winglets) should be necessary only to connect the pipe to the RF system.

The absorbers are metallic devices shaped like a trapezoid, with a total length of 30 cm, and they are inserted inside the chamber winglets, at about 42.5 mm from the beam axis. Placing slots for vacuum pumps just in front of the absorber allows to efficiently capture the synchrotron radiation and the molecule desorption. The pumping slots have a racetrack profile, length of 100-120~mm and width of 4-6~mm. Behind the slots, a cylindrical volume and a flange will be installed to support a NEG pump~\cite{kersevan}.

Impedance studies of the beam chamber profile with one absorber insertion have been performed using CST. In Fig.~\ref{Absorber}, the geometry of the FCC-ee beam chamber used in CST simulations is shown together with a detail of the absorber inside the beam chamber. Pumping slots and pumps are not included in this preliminary model.
Simulations show that below about 3~GHz the longitudinal impedance is purely inductive, giving, for 10000 elements, a longitudinal broadband impedance $Z/n$ of about 1~m$\Omega$.
\begin{figure}[!ht]
\center
\includegraphics[width=115mm]{Absorber.pdf}
  \caption{3D model of the FCC-ee vacuum chamber with winglets and a synchrotron radiation absorber used for CST simulations.}
\label{Absorber}
\end{figure}

In Fig.~\ref{wpabsorbers}, the wake potential of a 2 and 4 mm bunch length is shown. Even if further analysis is needed, and this first evaluation could overestimate the impedance, by multiplying this wake by the number of absorbers, about 10000 elements, we see that their contribution is much lower than the resistive wall one, as shown in Fig.~\ref{wp_old_abs}. About the transverse contribution to the impedance of a single absorber, this is so low that, up to now, we did not manage to obtain reliable results.
\begin{figure}[!ht]
\center
\includegraphics[width=115mm]{wpabsorbers.pdf}
  \caption{Wake potential of a single abosorber for 2 and 4 mm RMS bunch length from CST code.}
\label{wpabsorbers}
\end{figure}

Another important source of impedance is represented by the RF system. Several options are under investigation~\cite{calaga}, and in our study 400 MHz cavities with a single cell, as the one shown in Fig.~\ref{cavity} left side, have been considered. The wake potentials for 2 and 4 mm RMS bunch length, as given by the ABCI code~\cite{ABCI}, are shown on the right side of the same figure.
\begin{figure}[!ht]
\center
\includegraphics[width=80mm]{cell_shape.pdf}
\includegraphics[width=80mm]{cav_wake.pdf}
  \caption{400 MHz single cell cavity used in ABCI (left) and wake potentials for 2 and 4 mm RMS bunch length, as given by the ABCI code (red and blue lines), and by eq.~\ref{eq:cavity} (green and black).}
\label{cavity}
\end{figure}

In order to check the ABCI results, in the same plot we have also represented the analytical wake potentials, for the same bunch lengths, of a cavity with attached tubes at high frequencies, given by the expression~\cite{palumbo}
\begin{equation}
\label{eq:cavity}
W(x)=\tilde{W} |x|^{1/4}e^{-x}\left[I_{-1/4}(x)+\textrm{sign}(z)I_{1/4}(x) \right]
\end{equation}
with $x=[z/(2\sigma_z)]^2$. As we can see from the figure, the analytical curves are almost  superimposed to those given by ABCI.

By using the single cell, a group of 4 cavities are put together as shown in Fig.~\ref{rfsystem}, and connected to other two groups by two tapers. A total number of 100 cells has been considered as a global contribution to the impedance model, and, as a consequence, 25 double tapers. Of course the impedance produced by the tapers strongly depends on their length, which we have considered here to be of 500 mm.
\begin{figure}[!ht]
\center
\includegraphics[width=120mm]{rfsystem.pdf}
\caption{Sketch, not in scale, of a group of 4 cells with tapers at each end.}
\label{rfsystem}
\end{figure}

The wake potential of a single double taper (in and out, considered independent) in this condition is shown in Fig.~\ref{tapers}. The results have been obtained with ABCI and another electromagnetic code, ECHO~\cite{echo,stupakov1}, and the agreement can be considered satisfactory.
\begin{figure}[!ht]
\center
\includegraphics[width=120mm]{tapers.pdf}
\caption{Wake potentials of a double taper for 2 and 4 mm RMS bunch length, as given by the ABCI code (red and blue lines), and by the ECHO code (green and black).}
\label{tapers}
\end{figure}

Since the RF system (including the connecting tapers) and the absorbers are very important sources of impedance, in particular these last ones for their high number, in Fig.~\ref{waketot} we have represented the total longitudinal wake potential of these sources for 2 and 4 mm bunch length compared with the resistive wall contribution. In the figure the red and the blue curves have to be compared with the green and the black one respectively. We can see a factor of about 3 between the resistive wall contribution and the other sources, indicating that, up to now, the resistive wall remains the main source of impedance and wake fields.
\begin{figure}[!ht]
\center
\includegraphics[width=120mm]{waketot.pdf}
\caption{Total wake potentials for 2 and 4 mm RMS bunch length due to the RF system and the absorbers (green and black lines) compared with the resistive wall contribution (red and blue lines).}
\label{waketot}
\end{figure}

Of course, also the transverse contribution of the previous devices has to be taken into account. However, since their longitudinal contribution is quite smaller than the resistive wall, even if a more careful evaluation has to be carried out, we do not expect surprising results. There are several other sources of impedance, such as bellows, RF fingers, BPMs and other devices for diagnostics, and the criterion for their design should be such as to generate an impedance much smaller with respect to the resistive wall part.

Possible trapped modes in the interaction region deserve special studies, and work on other collective effects, such as the fast ion and the electron cloud instabilities, is in progress.

\section{Conclusions and outlook}

In this paper we have discussed single beam collective effects in FCC-ee due to the beam coupling impedance. In particular we focused our study primarily on the resistive wall effects because this is, up to now, the main source of impedance.

We have found that, in the single bunch case, the transverse mode coupling instability is about a factor 6 higher than the nominal bunch population at the lowest energy (45.6 GeV), and even higher for the other machine energies. Also the microwave instability due to the resistive wall has a margin of safety of about 2.4 with respect to the bunch population. However these numbers will be reduced as other sources of impedance will be taken into account.

Regarding the multi-bunch effect, we have concluded that the resistive wall transverse coupled bunch instability has to be counteracted by a feedback system, which requires innovative ideas for its design. Also the quadrupolar resistive wall wake fields, due to the elliptic vacuum chamber, produce a very high tune variation, which would not exist for the circular geometry. For the longitudinal case, at this stage, it is not possible to evaluate the characteristics of trapped HOMs, but an estimate of the maximum allowed shunt impedance as a function of the resonant frequency has been given.

In addition to the assessment of the resistive wall effects, we have started the evaluation of the impedance budget for other devices, with the goal of designing them in order to give an impact on the beam dynamics much smaller than that of the resistive wall itself. With an accelerator of 100 km of length, this is a long work, and the strategy is to identify the most important sources of impedance. We have started with the synchrotron radiation absorbers and the RF system. For the former, it resulted that the initial proposal could not be accepted, and, consequently, a new geometry of the beam pipe has been presented. Both contributions together amount to about one third of the resistive wall one in the longitudinal plane. Their evaluation in the transverse plane is in progress.

Other machine devices are going to be studied with the goal of comparing their contribution to the resistive wall impedance. Also two-stream instabilities, such as the fast ion and the electron cloud, must be determined.

\section*{Acknowledgements}

We acknowledge many helpful and stimulating discussions with R. Calaga, R. Kersevan, K. Oide, E. Shaposhnikova, G. Stupakov, D. Zhou, and F. Zimmermann.


\section*{Acknowledgments}
This work has been supported by MIUR project P.O.N.~RICERCA E COMPETITIVITA' 2007-2013 - Avviso n.~713/Ric.~del 29/10/2010, Titolo II - ``Sviluppo/Potenziamento di DAT e di LPP'' (project n.~PON02-00576-3333585). M.D., A.G. and F.V.P. acknowledge support by Istituto Nazionale di Fisica Nucleare (INFN) through the project QUANTUM. V.T. acknowledges the support, during the first part of this project, of the German Space Agency DLR with funds provided by the Federal Ministry of Economics and Technology (BMWi) under grant no. DLR 50 WM 1556.

\section*{Bibliography}