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

Let $\mu(n)$ denote the Moebius arithmetic function and $M(x)$ its sum
function. Von Mangoldt first established the estimate $M(x) = o(x)$,
essentially going through the steps used in proving the Prime Number
Theorem (PNT).  A few years later, Landau showed by relatively simple
real variable arguments that this and several other estimates followed
from the PNT \cite[\S150]{landau}.  Similarly, these relations imply each
other and the PNT; thus they are said to be ``equivalent'' to the PNT.

In this article we consider an analog of the $M$ bound for Beurling
generalized (g-)numbers.  A Beurling g-prime system is an unbounded
sequence of real numbers ${\mathcal{P}} =\{p_i\}$ satisfying $1 < p_1 \le p_2
\le \ldots\,$, and the multiplicative semigroup generated by $\mathcal{P}$ and
1 is called the associated collection $\mathcal{N}$ of g-integers
(cf.~\cite{BD69}, \cite{diamond-zhangbook}, \cite{MV07}).  The
counting function $N(x)$ of $\mathcal{N}$ is the number of elements of $\mathcal{N}$
not exceeding $x$.  The g-Chebyshev function of $\mathcal{N}$ is
\[
\psi(x) = \sum_{p_i^{\alpha_i} \le x} \log p_i.
\]
(Here and below we give our Beurling functions the same names used in
classical number theory.)

The PNT-related assertions are somewhat different
for g-numbers: not all implications between the several corresponding
assertions hold unconditionally, see e.g.~\cite{d-vPNTequiv2016},
\cite[Chap.~14]{diamond-zhangbook}, \cite{zhang1987}.  Here we study
the Beurling version of the assertion $M(x) = o(x)$ and show this can
be deduced (a) from the PNT under an O-boundedness condition
$N(x)=O(x)$ or (b) from a Chebyshev-type upper bound assuming
$N(x)\sim ax$ (for $a>0$) and the integral condition \eqref{L1cond}
(below).  At the end, we give examples in which $M$ estimates fail.

A word about the Beurling version of $M(x)$.  The characteristic
property of the Moebius function is its being the (multiplicative)
convolution inverse of the $1$ function.  For g-numbers, we define the
measure $\mathrm{d}M$ as the convolution inverse of $\mathrm{d}N$; by familiar Mellin
transform properties (cf. \cite{BD04}, \cite{diamond-zhangbook},
\cite{MV07}),
\[
 \int_{1-}^\infty u^{-s} \mathrm{d}M(u) = 1/\zeta(s) \ \, {\rm with} \ \,
\zeta(s) = \zeta_{\mathcal{N}}(s) = \int_{1-}^\infty u^{-s} \mathrm{d}N(u),
\]
the Beurling number version of the Riemann zeta function. Note that
these definitions make sense even when $\mathrm{d}N$ is not discrete or when
factorization into primes is not unique.
\smallskip

We thank Wen-Bin Zhang for his useful comments.

\section{PNT hypothesis}  \label{sec:PNT}

In this section we show that the PNT together with the ``O-density''
condition $N(x) \ll x$ implies $M(x) = o(x)$.
\begin{theorem}  \label{th:M}
Let $\mathcal{N}$ be a g-number system for which the PNT holds and $N(x) \ll
x$.   Then $M(x) = o(x)$.
\end{theorem}

Other sufficient conditions for an $M$ estimate are known,
e.g.~\cite[ Prop.~14.10]{diamond-zhangbook}, if the PNT holds for a
g-number system and the integer counting function satisfies the
logarithmic density condition
\begin{equation}
\label{log-density}
\int_{1^{-}}^{x}\frac{\mathrm{d}N(t)}{t}\sim a\log x,
\end{equation}
then $M(x) = o(x)$.  The present result differs from the other in that
log-density and O-density are conditions that do not imply one another;
also, the proofs are very different.

The key to our argument is the following relation.
\begin{lemma} \label{lem:key}
Under the hypotheses of the theorem,
\begin{equation}  \label{eq:main-ident}
\frac{M(x)}{x} = \frac {-1}{\log x}\int_1^x \!\frac{M(t)}{t^2}\,
\mathrm{d}t + o(1).
\end{equation}
\end{lemma}

\begin{proof}[Proof of the lemma]
A variant of Chebyshev's identity for primes reads
\[
L\mathrm{d}M =  -\mathrm{d}M * \mathrm{d}\psi,
\]
where $L$ is the operator of multiplication by $\log t$ and $*$ is
multiplicative convolution.  This can be verified (in the classical or
in the Beurling case) via a Mellin transform.  Note that this
transform carries convolutions into pointwise products and the $L$
operator into differentiation.  The equivalent Mellin formula is the
identity
\[
\Big\{\frac{1}{\zeta(s)} \Big\}' = -\frac{1}{\zeta(s)} \cdot
\frac {\zeta'(s)}{\zeta(s)}.
\]

Now add and subtract the term $(\delta_1+\mathrm{d}t)*\mathrm{d}M$ in
the variant of the
Chebyshev relation, with $\delta_1$ the point mass at 1 and $\mathrm{d}t$
the Lebesgue measure on $(1,\, \infty)$.  Integrating, we find
\begin{equation}  \label{eq:Cheb'}
\int_{1-}^x L\,\mathrm{d}M =  \int_{1-}^x\mathrm{d}M *
(\delta_1+\mathrm{d}t - \mathrm{d}\psi) - \int_{1-}^x\mathrm{d}M
* (\delta_1+\mathrm{d}t).
\end{equation}

Integrating by parts the left side of the last formula, we get
\[
M(x) \log x - \int_1^x \frac{M(t)}{t}\, \mathrm{d}t = M(x) \log x + O(x),
\]
since $|\mathrm{d}M| \le \mathrm{d}N$, and hence $|M(t)| \le N(t) \ll
t$.  If we evaluate the convolution integrals by the iterated integral
formula
\[
 \int_{1-}^x\mathrm{d}A * \mathrm{d}B = \int\!\!\!\int_{st \le x}\!
 \mathrm{d}A(s)\, \mathrm{d}B(t) = \int_{1-}^xA(x/t)\, \mathrm{d}B(t)
\]
and use the PNT, we find the first term on the right side of
\eqref{eq:Cheb'} to be
\[
\int_{1-}^x \{x/t-\psi(x/t)\}\,\mathrm{d}M(t) =
\int_{1-}^x o(x/t)\,\mathrm{d}N(t) = o(x \log x).
\]
Also, the last term of \eqref{eq:Cheb'} is, upon integrating by parts,
\[
\int_{1-}^x \frac xt \, \mathrm{d}M(t) = M(x) + x\!\int_1^x
\!\frac{M(t)}{t^2}\,\mathrm{d}t.
\]

Finally, we combine the bounds for the terms of \eqref{eq:Cheb'},
divide through by $x \log x$, and note that $M(x) \ll x$ to obtain
\eqref{eq:main-ident}.
\end{proof}

\begin{proof}[Proof of Theorem \ref{th:M}]
First, we can assume that $M(x)$ has an infinite number of sign
changes.  Otherwise, there is some number $z$ such that $M(x)$ is of
one sign for all $x \ge z$.  By \eqref{eq:main-ident}, as $x \to \infty$,
\[
\frac{M(x)}{x}
= \frac {-1}{\log x}\, \int_1^z
\!\frac{O(t)}{t^2}\,\mathrm{d}t -\frac{1}{\log x}\int_z^x
\!\frac{M(t)}{t^2}\,\mathrm{d}t  + o(1).
\]
Thus we have
\[
\frac{M(x)}{x}
+ \frac{1}{\log x}\int_z^x \!\frac{M(t)}{t^2}\,\mathrm{d}t
=  \frac {O(\log z)}{\log x} +o(1) = o(1).
\]
Since $M(x)$ and the integral are of the same sign, $M(x)/x \to 0$
as $x \to \infty$, and this case is done.

Now suppose that $M(x)$, which we regard as a right-continuous
function, changes sign at $x$.  We show that $M(x)/x = o(1)$.  If
$M(x)/x = 0$, there is nothing more to say here.  If, on the other
hand, $M(x) > 0$, then there is a number $y \in (x-1,\ x)$ with $M(y)
\le 0$.  If we apply \eqref{eq:main-ident} again, we find
\begin{align*}
\frac{M(x)}{x}
&= \frac {-\log y}{\log x}\, \frac{1}{\log y}\int_1^y
\!\frac{M(t)}{t^2}\,\mathrm{d}t -\frac{1}{\log x}\int_y^x
\!\frac{M(t)}{t^2}\,\mathrm{d}t  + o(1) \\
&= \frac {\log y}{\log x}\,\frac{M(y)}{y} +
\frac{1}{\log x} \int_y^x \frac{O(t)\,\mathrm{d}t}{t^2} +o(1).
\end{align*}
Thus $M(x)/x - M(y)/y  = o(1)$, and since $M(x)/x> 0 \ge M(y)/y$,
each is $o(1)$.  A similar story holds if $M(x) < 0$.

Finally, suppose that $M(t)$ changes sign at $t=y$ (so that
$M(y)/y = o(1)$) and  $M(t)$ is of one sign
for $y < t \le z$.  By yet another application of
\eqref{eq:main-ident},  we find for any $x \in (y,\, z]$
\[
\frac{M(x)}{x}
= \frac {-\log y}{\log x}\, \frac{1}{\log y}\int_1^y
\!\frac{M(t)}{t^2}\,\mathrm{d}t -\frac{1}{\log x}\int_y^x
\!\frac{M(t)}{t^2}\,\mathrm{d}t  + o(1)
\]
or
\[
\frac{M(x)}{x} + \frac{1}{\log x}\int_y^x \!\frac{M(t)}{t^2}\,\mathrm{d}t
= \frac {\log y}{\log x}\,\frac{M(y)}{y} +o(1) = o(1)
\]
as $y \to \infty$.  Since $M(x)/x$ and the integral are of the same
sign, it follows that $M(x)/x = o(1)$.
\end{proof}

\section{Chebyshev hypothesis}  \label{sec:cheb}

What happens if the PNT hypothesis of the last theorem is weakened to
just a Chebyshev upper bound? In the examples of Section \ref{Sect
  examples} we show that even two-sided Chebyshev estimates by
themselves are not strong enough to ensure that $M(x) = o(x)$ holds.
Furthermore, this bound could fail even if, in addition to Chebyshev
estimates, one also assumes that $N$ satisfies both
(\ref{log-density}) and $N(x)\ll x$, see Example \ref{ex1C-M}. We
shall show, however, that if the regularity hypothesis on $N$ is
slightly augmented, then one can indeed deduce the desired $M$ bound.

\begin{theorem} \label{thCheby-moebius} Suppose that a g-number
 system satisfies the following conditions: \begin{enumerate}
 \item [(a)] a Chebyshev upper bound, that is,
\begin{equation}
\label{chebyeq1} \limsup_{x\to\infty} \frac{\psi(x)}{x} < \infty,
\end{equation}
\item [(b)]  for some positive constant $a, \ N(x) \sim ax$,
\item [(c)] for some $\b \in (0,1/2)$ and all $\sigma \in (1,2)$
\begin{equation}
\label{L1cond} \int_{1}^{\infty} \frac{|N(x)-ax|}{x^{\sigma+1}}\:
\mathrm{d}x \ll (\sigma-1)^{-\b}\,.
\end{equation}
 \end{enumerate}
Then, $M(x)=o(x)$ holds.
\end{theorem}

As a simple consequence we have an improvement of a result of
W.-B.~Zhang \cite[Cor.~2.5]{zhang1987}:
\begin{corollary}  If a g-number system satisfies \eqref{chebyeq1} and
\[
N(x) - ax \ll x \log ^{-\gamma} x
\]
for some $\gamma > 1/2$, then $M(x)=o(x)$ holds.
\end{corollary}

 W.-B. Zhang had conjectured\footnote{Oral communication to the authors.} that a Chebyshev
 bound along with the $L^{1}$ bound
\begin{equation}
\label{strongL1cond} \int_{1}^{\infty} \frac{|N(x)-ax|}{x^{2}}
\mathrm{d}x <\infty
\end{equation}
imply $M(x) = o(x)$.  Naturally, (\ref{strongL1cond}) is included in
(\ref{L1cond}) and $N(x)\sim ax$.

It is easy to show that $M(x) = o(x)$ is always implied by
\begin{equation} \label{eqlandausharprelation}
 m(x) =\int_{1^{-}}^{x}\frac{\mathrm{d}M(u)}{u}=

 o(1),
\end{equation}
but the converse implication is not true in general
\cite{diamond-zhangbook}.  On the other hand, it has recently been
shown \cite[Thm. 2]{d-vPNTequiv2016} that the assertions become
equivalent under the additional hypothesis (\ref{strongL1cond}).  So, Theorem
\ref{thCheby-moebius} and the quoted result yield at once:

\begin{corollary} \label{corCheby-moebius} Under (\ref{chebyeq1}) and
  (\ref{strongL1cond}), relation (\ref{eqlandausharprelation})
  holds as well.
\end{corollary}

We can also strengthen another result of W.-B. Zhang
\cite[Thm. 2.3]{zhang1987}.

\begin{corollary} \label{cor2Cheby-moebius} The condition
(\ref{strongL1cond})
and
\begin{equation}\label{avcondpointwiseeq}
\int_{1}^{x}\frac{(N(t)-at)\log t}{t}\:\mathrm{d}t\ll x
\end{equation}
imply (\ref{eqlandausharprelation}).
\end{corollary}
\begin{proof} Chebyshev bounds are known to hold under the hypotheses
  (\ref{strongL1cond}) and (\ref{avcondpointwiseeq})
  \cite[Thm. 11.1]{diamond-zhangbook}. The rest follows from Corollary
  \ref{corCheby-moebius}.
\end{proof}

We shall prove Theorem \ref{thCheby-moebius} using several lemmas.
Our method is inspired by W.-B. Zhang's proof of a H{a}l\'{a}sz-type
theorem for Beurling primes \cite{zhang1987}.  Our first step is to
replace $M(x)=o(x)$ by an equivalent asymptotic relation.

\begin{lemma}\label{lem1C-M} Suppose that the integer counting
  function of a g-number system has a positive density, i.e.,
  $N(x)\sim a x$ for some $a>0$. Then, $M(x)=o(x)$ if and only if
\begin{equation}
\label{MoebiusAv} f(x)=\int_{1}^{x}\left(\int_{1}^{u}\log t\:
  \mathrm{d}M(t) \right)\frac{\mathrm{du}}{u}=o(x\log x).
\end{equation}
\end{lemma}
\begin{proof}The direct implication is trivial. For the converse, we
  set
\[
g(x)=\int_{1^{-}}^{x} \log t\, (\mathrm{d}N(t)+\mathrm{d}M(t)).
\]
Notice that
\[
\mathrm{d}N+\mathrm{d}M=2\sum_{n=0}^{\infty}\,
\mathrm{d}\Pi^{\ast 2n}/(2n)!
\]
is a non-negative measure, so that $g$ is non-decreasing. The
hypotheses (\ref{MoebiusAv}) and $N(x)\sim ax$ give
$$
\int_{1}^{x}\frac{g(u)}{u}\,\mathrm{d}u\sim  a x\log x.
$$
Hence,
$$
\int_{1}^{x}g(u)\,\mathrm{d}u= x \int_{1}^{x}\frac{g(u)}{u} \,\mathrm{d}u
- \int_{1}^{x}\!\left(\int_{1}^{t}\frac{g(u)}{u} \,\mathrm{d}u\right)
\mathrm{d}t \sim \frac{a}{2} x^{2}\log x.
$$
Since $g$ is a non-decreasing function, we see by a simple
differencing argument (see e.g. \cite[p. 34]{korevaarbook}) that
$g(x)\sim a x\log x$. But also we have $\int_{1^{-}}^{x}\log u\,
\mathrm{d}N(u)\sim a x\log x $; consequently, $\int_{1^{-}}^{x}\log
u\, \mathrm{d}M(u) = o( x\log x)$.  Integration by parts now yields
$$
\log x\, M(x)= \int_{1}^{x}\log u\, \mathrm{d}M(u)
+\int_{1}^{x} \frac{M(u)}{u}\, \mathrm{d}u = o(x\log x)+O(x),
$$
and the result then follows by dividing by $\log x$.
\end{proof}

The next lemma provides a crucial analytic estimate.
\begin{lemma}\label{lem3CC-M} Suppose that $N(x)\sim a x$, with $a>0$. Then,
\begin{equation}
\label{invzeta}
\frac{1}{\zeta(\sigma+it)}=o\left(\frac{1}{\sigma-1}\right), \quad
\sigma \to 1^{+},
\end{equation}
uniformly for $t$ on compact intervals.
\end{lemma}
\begin{proof} We first show that (\ref{invzeta}) holds pointwise,
i.e., for each fixed $t\in \mathbb{R}$, without the uniformity
requirement. If $t=0$ this is clear because $\zeta(\sigma)\sim
a/(\sigma-1)$ and thus $1/\zeta(\sigma)=o(1)$. Note that $N(x)\sim a x$
implies
$$
\zeta(s)=\frac{a}{s-1}+o\left(\frac{|s|}{\sigma-1}\right)
$$
uniformly. If $t\neq 0$, we obtain, $\zeta(\sigma+2it)=o_{t}(1/(\sigma-1))$.
Applying the 3-4-1 inequality, we conclude that
$$
1\leq |\zeta(\sigma)|^{3}|\zeta(\sigma+it)|^{4} |\zeta(\sigma+2it)|
=|\zeta(\sigma+it)|^{4} \,o_{t}\!\left(\frac{1}{(\sigma-1)^{4}}\right),
$$
which shows our claim. Equivalently, we have
$$
\exp\left(-\int_{1}^{\infty}x^{-\sigma}(1+\cos(t\log x))
\,\mathrm{d}\Pi(x)\right)
=\frac{1}{|\zeta(\sigma)\zeta(\sigma+it)|}=o_{t}(1).
$$
The left-hand side of the last formula is a net of continuous
monotone functions in the variable $\sigma$ that tend pointwise to 0 as $\sigma\to 1^{+}$;
Dini's theorem then asserts that it must also converge to 0 uniformly
for $t$ on compact sets, as required.
\end{proof}

The next lemma is very simple but useful.

\begin{lemma}\label{lem3C-M}Let $F$ be a right continuous function of
  local bounded variation with support in $[1,\infty)$ satisfying the
  bound $F(x)=O(x)$. Set $\widehat{F}(s)=\int_{1^{-}}^{\infty}x^{-s}
  \mathrm{d}F(x)$, $\Re s>1$. Then,
$$
\int_{\Re s=\sigma} \left| \frac{\widehat{F}(s)}{s}\right|^{2}
|\mathrm{d}s|\ll \frac{1}{\sigma-1} \,.
$$
\end{lemma}
\begin{proof}
Indeed, by the Plancherel theorem,
$$
\int_{\Re s=\sigma}  \left| \frac{\widehat{F}(s)}{s}\right|^{2}
|\mathrm{d}s|=2\pi\! \int_{0}^{\infty}\!\! e^{-2\sigma x}
|F(e^{x})|^{2}\,\mathrm{d}x\ll \int_{0}^{\infty}\!\! e^{-2(\sigma-1) x}
\,\mathrm{d}x= \frac{1}{2(\sigma-1)}\,. \qedhere
$$
\end{proof}

If a g-number system satisfies a Chebyshev upper bound, then Lemma
\ref{lem3C-M} implies
\begin{equation}
\label{L2Chevy}
\int_{\Re s=\sigma} \left|\frac{\zeta'(s)}{s\zeta(s)}\right|^{2}
|\mathrm{d}s|\ll \frac{1}{\sigma-1}.
\end{equation}
Similarly, $N(x)=O(x)$ yields
\begin{equation}
\label{L2zeta}
\int_{\Re s=\sigma} \left|\frac{\zeta(s)}{s}\right|^{2} |\mathrm{d}s|
\ll \frac{1}{\sigma-1}.
\end{equation}

Also, we shall need the following version of the Wiener-Wintner theorem
\cite{bateman-diamond2000,montgomeryBook1970}.
\begin{lemma}\label{lem4C-M} Let $F_{1}$ and $F_{2}$ be right
continuous functions of local bounded variation with support in
$[1,\infty)$. Suppose that their Mellin-Stieltjes transforms
$\widehat{F}_{j}(s)=\int_{1^{-}}^{\infty}x^{-s} \mathrm{d}F_{j}(x)$
are convergent on $\Re s>\alpha$, that $F_{1}$ is non-decreasing, and
$|\mathrm{d}F_{2}|\leq \mathrm{d}F_{1}$. Then, for all
$b\in\mathbb{R}$, $c>0$, and $\sigma>\alpha$,
$$
\int_{b}^{b+c}|\widehat{F}_{2}(\sigma+it)|^{2}\,\mathrm{d}t
\leq 2 \int_{-c}^{c}|\widehat{F}_{1}(\sigma+it)|^{2}\,\mathrm{d}t.
$$
\end{lemma}

We are ready to present the proof of Theorem \ref{thCheby-moebius}.
\begin{proof}[Proof of Theorem \ref{thCheby-moebius}]
In view of Lemma
  \ref{lem1C-M}, it suffices to show (\ref{MoebiusAv}).
The Mellin-Stieltjes transform of the function $f$ is
$$
-\frac{1}{s}\left(\frac{1}{\zeta(s)}\right)'.
$$
Given $x > 1$, it is convenient to set $\ss=1+1/\log x$.  By the
Perron inversion formula, the Cauchy-Schwarz inequality, and
(\ref{L2Chevy}), we have
\begin{align*}
\nonumber
\frac{|f(x)|}{x}&=\frac{1}{2\pi}\left|\int_{\Re s=\ss}
\frac{x^{s-1}\zeta'(s)}{s^{2}\zeta^{2}(s)}\, \mathrm{d}s
\right|\leq \frac{e}{2\pi} \int_{\Re s=\ss}
\left|\frac{\zeta'(s)}{s^{2}\zeta^{2}(s)}\right| |\mathrm{d}s|  \\
&
\ll \log^{1/2}x \left( \int_{\Re s=\ss} \left|
\frac{1}{s\zeta(s)}\right|^{2} |\mathrm{d}s| \right)^{1/2}.
\end{align*}

Next, we take a large number $\lambda$, arbitrary but fixed. We
split the integration line $\{\Re s=\ss\}$ of the last integral into two
parts, $\{\ss+it:\: |t|\geq \lambda\} $ and $\{\ss+it:\: |t|\leq
\lambda\} $, and we denote the corresponding integrals over these sets
as $I_{1}(x)$ and $I_{2}(x)$ respectively, so that
\begin{equation}
\label{eq1proofM}
\frac{|f(x)|}{x}\ll \big((I_{1}(x))^{1/2}+(I_{2}(x))^{1/2}\big) \log^{1/2}x.
\end{equation}

To estimate $I_{1}(x)$, we apply Lemma \ref{lem4C-M} to
$|\mathrm{d}M|\leq \mathrm{d}N$ and employ  (\ref{L2zeta}),
$$
\left(\int_{-\lambda -m-1}^{-\lambda-m} +\int_{\lambda+m}^{\lambda+m+1}
\right)\left|\frac{1}{\zeta(\ss+it)}\right|^{2} \mathrm{d}t\leq 4
\int_{-1}^{1} \left|\zeta(\ss+it)\right|^{2} \mathrm{d}t \ll \log
x.
$$
Therefore, we have the bound
\begin{equation}
\label{eq2proofM}
I_{1}(x)\ll \log x \sum_{m=0}^{\infty} \frac{1}{1+m^{2}
+\lambda^{2}}\ll \frac{\log x}{\lambda}.
\end{equation}

To deal with $I_{2}(x)$, we need to derive further properties of the
zeta function. Using the hypothesis \eqref{L1cond}, we find that
\begin{align*}
\zeta(s)-\frac{a}{s-1}&= s \!\int_{1}^{\infty}\!\!
x^{-s}\frac{N(x)-ax}{x}\:\mathrm{d}x +a \\
&\ll |s|\!\int_{1}^{\infty}\!\! x^{-\sigma}\frac{|N(x)-ax|}{x}\,
\mathrm{d}x =O_t\!\left(\frac{1}{(\sigma-1)^{\beta}}\right).
\end{align*}
Hence, we obtain
\begin{equation}
\label{eq3proofM}
\zeta(\ss+it)= \frac{a}{\ss-1+it}
+ O(\log^{\beta} x)
\end{equation}
for some number $\beta \in (0,\,1/2)$, uniformly for $t$ on
compact sets. We are ready to estimate $I_{2}(x)$. Set
$\eta=(1-2\beta)/(1-\beta)$ and note that $\eta \in (0,\, 1)$.  Then,
using Lemma~\ref{lem3CC-M},
$$
I_{2}(x)\leq \int_{-\lambda}^{\lambda} \left|
\frac{1}{\zeta(\ss+it)}\right|^{2}\mathrm{d}t\leq
\left(\int_{-\lambda}^{\lambda} \left|
\frac{1}{\zeta(\ss+it)}\right|^{2-\eta}\mathrm{d}t\right)
 o_{\lambda}(\log^{\eta}x).
$$
On the other hand, applying  Lemma \ref{lem4C-M} to
$\mathrm{d}F_{1}=\exp^{\ast}(-(1-\eta/2)\,\mathrm{d}\Pi)$ and
$\mathrm{d}F_{2}=\exp^{\ast}((1-\eta/2)\,\mathrm{d}\Pi)$
and using  (\ref{eq3proofM}), we find
\begin{align*}
\int_{-\lambda}^{\lambda}\left|\frac{1}{\zeta(\ss+it)}\right|^{2-\eta}\mathrm{d}t
&
\leq 4\int_{-\lambda}^{\lambda} \left|\zeta(\ss+it)\right|^{2-\eta}\mathrm{d}t
\\
&
 \ll  \int_{-\lambda}^{\lambda} \frac{\mathrm{d}t}{ ((\ss-1)^{2}+t^{2})^{1-\eta/2}}
+\lambda\log^{(2-\eta)\beta}x\ll \lambda \log^{1-\eta}x,
\end{align*}
which implies $I_{2}(x)= o_{\lambda}(\log
x)$. Inserting this and the bound (\ref{eq2proofM}) into
(\ref{eq1proofM}), we arrive at
$
|f(x)|/(x\log x)\ll \lambda^{-1/2} +o_{\lambda}(1).$ Taking first the limit
superior as $x\to\infty$ and then $\lambda\to\infty$, we have shown
that
$$
\lim _{x\to\infty} \frac{f(x)}{x\log x}=0.
$$
By Lemma~\ref{lem1C-M}, $M(x) = o(x)$.
This completes the proof of Theorem \ref{thCheby-moebius}.
\end{proof}

\section{Three examples}\label{Sect examples}
The examples of this section center on the importance of
$N(x)$ being close to $ax$ in Theorem \ref{thCheby-moebius}.  In
the first example, $N(x)/x$ has excessive wobble and in the second
one, excessive growth; in both cases $M(x) = o(x)$ fails. The third
example shows that condition~\eqref{L1cond} is not sufficient to
insure the convergence of $N(x)/x$, whence the introduction of this
hypothesis.

In preparation for treating the first two examples, we give a
necessary condition for $M(x) = o(x)$. An analytic function $G(s)$ on
the half-plane $\{s:\:\Re s>\alpha\}$ is said to have a right-hand zero of
order $\beta>0$ at $s=it_{0}+\alpha$ if $\lim_{\sigma\to\alpha^{+}}
(\sigma -\alpha)^{-\beta}G(it_{0}+\sigma)$ exists and is non-zero.
Our examples violate the following necessary condition:
\begin{lemma}
\label{lexC-M}
If $M(x)=o(x)$, then $\zeta(s)$ does not have any right-hand zero of
order $\geq 1$ on $\{\Re s=1\}$.
\end{lemma}
\begin{proof} We must have
$$
\frac{1}{\zeta(\sigma+it)}=s\!\int_{1}^{\infty} x^{-it-\sigma}
\frac{M(x)}{x}\:\mathrm{d}x= o\!\left(\frac{1}{\sigma-1}\right), \quad
\sigma\to 1^{+},
$$
uniformly for $t$ on compact intervals.
\end{proof}

\begin{example}
\label{ex1C-M}
We consider
$$
\Pi (x)=\sum_{2^{k+1/2}\leq x} \frac{2^{k+1/2}}{k}\,.
$$
This satisfies the Chebyshev bounds: we have
\[
\psi(x)=  2^{\lfloor \log x/\log 2 + 1/2\rfloor +1/2}\log 2 +O(x/\log x).
\]
Thus
$$
\liminf_{x\to\infty} \frac{\psi(x)}{x}=\log 2 \quad \mbox{and} \quad
\limsup_{x\to\infty} \frac{\psi(x)}{x}=2 \log 2.
$$
Further, the zeta function of this g-number system can be
explicitly computed:
$$
\log \zeta(s)= 2^{-(s-1)/2}\sum_{k=1}^{\infty} \frac{2^{-k(s-1)}}{k} = - 2^{-(s-1)/2}\log(1-2^{-(s-1)}).
$$
We conclude that
$$\zeta(\sigma)\sim \frac{1}{(\sigma-1)\log 2}, \quad \sigma\to1^{+};$$
therefore, by the Hardy-Littlewood-Karamata Tauberian Theorem
\cite{diamond-zhangbook}, \cite{korevaarbook}, $N$ has logarithmic density
\begin{equation}
\label{log-densityex1}
\int_{1^{-}}^{x} \frac{\mathrm{d}N(u)}{u}\sim \frac{\log x}{\log 2}.
\end{equation}
Furthermore, $\zeta(s)$  has infinitely many right-hand zeros of order
1 at the points $s=1\pm i 2\pi (2n+1) /\log 2$, $n\in\mathbb{N}$, because
\[
\zeta\left(\sigma \pm i\frac{2\pi(2n+1)}{\log 2}\right)= \frac{1}{
  \zeta (\sigma)}\sim (\sigma-1)\log 2.
\]
It follows, by Lemma \ref{lexC-M}, that
$$M(x)=\Omega(x).$$

To show the wobble of $F(x) = N(x)/x$, we apply an idea of Ingham
\cite{ingham1942} that is based on a finite form of the Wiener-Ikehara
method.  We use (essentially) the result given in \cite[Thm.~11.12]{BD04}.
For $\Re s > 0$,
\[
\int_1^\infty x^{-s-1} F(x)\,\mathrm{d}x = \frac{\zeta(s+1)}{s+1} = G(s).
\]
The discontinuities of $G(s)$ on the line segment $(-8\pi i / \log 2,
8\pi i / \log 2)$, which provide a measure of the wobble, occur at
$s=0$ and $s = \pm 4\pi i / \log 2$.

We analyze the behavior of $G$ near these points.  Let $s$ be a
complex number with $\Re s \ge 0$ (to avoid any logarithmic fuss)
and $0 <|s| \le 1/2$ (so $\log s \ll \log |1/s|$ is valid).
 For $n = 0, \, \pm 1$, a small calculation shows that
\[
\zeta(1+ s + 4 \pi n i/\log 2) = \frac{1}{s \log 2}
+ O(\log |1/s|).
\]
For $n = -1,\,0, \, 1$ set
\[
\gamma_n =  4 \pi n/\log 2, \quad \  \a_n = 1/(\log 2+ 4 \pi in).
\]
Take $T$ a number between $4\pi/\log 2$ and $8\pi/\log 2$, e.g.~$T=36$,
and set
\[
G^*(s) =  \sum_{-1\le n \le 1}\frac{\a_{n}}{s-i \gamma_n}, \quad
F_T^*(u) = \sum_{-1\le n \le 1}\!\!\a_{n}\Big(1 - \frac{|\gamma_{n}|}{T}
\,\Big)\, e^{i \gamma_nu}.
\]

Now $G - G^*$ has a continuation to the closed strip $\{s\!: \sigma
\ge 0, \, |t| \le T\}$ as a function that is continuous save
logarithmic singularities at ${\gamma}_{-1},\,{\gamma}_0,\,{\gamma}_1$.  In particular,
$G - G^*$ is integrable on the imaginary segment $(-iT, \, iT)$.  The
result in \cite{BD04} is stated for an extension that is continuous at
all points of such an interval, but integrability is a sufficient
condition for the result to hold.

We find that
\[
\liminf_{u \to \infty} F(u) \le \inf_u F_T^*(u) < \sup_uF_T^*(u)
\le \limsup_{u \to \infty} F(u).
\]
By a little algebra,
\begin{align*}
\sup_uF_T^*(u) &= \frac 1{\log 2} + \frac 2{|\log 2 - 4 \pi i|}\Big(1 -
\frac{4 \pi}{36 \log 2}\Big) > 1.52, \\
\inf_uF_T^*(u) &= \frac 1{\log 2} - \frac 2{|\log 2 - 4 \pi i|}\Big(1 -
\frac{4 \pi}{36 \log 2}\Big) < 1.37.
\end{align*}
Thus $N(x)/x$ has no asymptote as $x \to \infty$.

It is interesting to note that $\zeta$ is $4\pi i/\log 2$ periodic.
So $\zeta(s)$ has an analytic continuation to $\{s\!: \Re s=1,\,
s\neq1+4n\pi i/\log 2,\: n\in\mathbb{N}\}$. One could show a larger
oscillation by a more elaborate analysis exploiting additional
singularities.
\medskip

Finally, we discuss $m(x)=\int_{1^{-}}^{x}u^{-1} \mathrm{d}M(u)$ for
this example.  Since $M(x)=\Omega(x)$, we necessarily have
$m(x)=\Omega(1)$. We will prove that
\begin{equation}
\label{m O-bound}
m(x)=O(1).
\end{equation}
This shows that, in general, having Chebyshev bounds, log-density
(\ref{log-density}), $N(x)\ll x $, and (\ref{m O-bound}) together do
not suffice to deduce the estimate $M(x)=o(x)$.

To prove (\ref{m O-bound}), we first need to improve (\ref{log-densityex1}) to
\begin{equation}
\label{2log-densityex1}
\int_{1^{-}}^{x} \frac{\mathrm{d}N(u)}{u}= \frac{\log x}{\log 2}
-\frac{\log \log x}{2}+O(1).
\end{equation}
 This estimate can be shown by applying a Tauberian theorem of
 Ingham-Fatou-Riesz type \cite{debruyne-vindasCT}
 (cf.~\cite{ingham1935,korevaarbook}). In fact, the Laplace transform
 of the non-decreasing function $\tau_{1}(x)=\int_{1^{-}}^{e^{x}}
 u^{-1} \mathrm{d}N(u)$ is analytic on $\Re s>0$ and
\begin{align*}
G(s)&=\mathcal{L}\{\tau_{1}; s\}- \frac{1}{s^{2}\log 2}
+\frac{\log (1/s)}{2s}=\frac{\zeta(s+1)}{s}-\frac{1}{s^{2}\log 2}
+\frac{\log (1/s)}{2s} \\
&
= \frac{1+\log\log2}{2s}+O(\log^{2} |1/s|), \quad |s|<1/2,
\end{align*}
as a small computation shows. In the terminology of
\cite{debruyne-vindasCT}, $G(s)$ has local pseudo\-measure boundary
behavior on the imaginary segment $(-i/2,i/2)$; in fact, its boundary
value on that segment is the sum of a pseudomeasure and a locally
integrable function.  One then deduces (\ref{2log-densityex1})
directly from \cite[Thm. 3.7]{debruyne-vindasCT}.
Similarly, we use the fact that $(s\,\zeta(s+1))^{-1}$ has a continuous
extension to the same imaginary segment and we apply the same
Tauberian result to the non-decreasing function
$\tau_{2}(x)=\int_{1^{-}}^{e^{x}} u^{-1}
(\mathrm{d}M(u)+\mathrm{d}N(u))$, whose Laplace transform is
$\mathcal{L}\{\tau_{2}; s\}= \mathcal{L}\{\tau_{1}; s\}+(s\zeta(s+1))^{-1}.$
The conclusion is again the asymptotic formula $\tau_{2}(\log x)= \log
x/\log 2-(\log \log x)/2+O(1)$.  One then obtains (\ref{m O-bound}) upon
subtracting  (\ref{2log-densityex1}) from this formula.
\end{example}

\begin{example}
\label{ex2C-M} As a second example, we consider a modification of the
Beurling-Diamond examples from \cite{beurling1937,diamond1970} (see
also \cite{d-s-v}), namely, the continuous prime measure
$\mathrm{d}\Pi_{B}$ given by
$$
\Pi_{B}(x)= \int_{1}^{x}\frac{1-\cos (\log u)}{\log u}\: \mathrm{d}u
$$
and the discrete g-prime system
$$q_{k}=\Pi^{-1}_{B}(k), \quad k=1,2,\dots,$$
with g-prime and g-integer counting functions $\pi_{D}(x)$ and $N_{D}(x)$.

We study here the continuous prime measure
$\mathrm{d}\Pi_{C}=2\mathrm{d}\Pi_{B}$ and the discrete $g$-primes
formed by taking each $q_{k}$ twice, that is, the g-prime system
$$\mathcal{P}= \{q_{1},\, q_{1},\, q_{2},\, q_{2},\, q_{3},\,
q_{3},\,\dots\}.$$
The associated number-theoretic functions will be denoted as
$\pi_{\mathcal{P}}(x)$, $\Pi_{\mathcal{P}}(x)$, $N_{\mathcal{P}}(x)$,
$M_{\mathcal{P}}(x)$, and $\zeta_{\mathcal{P}}(s)$, and those
corresponding to $\mathrm{d}\Pi_{C}$ we denote by $N_{C}(x)$, $M_{C}(x)$, and
$\zeta_{C}(s)$.

It is easy to verify that
\begin{equation*}
\Pi_{C}(x)=  \frac{x}{\log x}\left(2-\sqrt{2}\cos \left(
\log x-\frac{\pi}{4}\right)\right) +O\left(\frac{x}{\log^{2} x}\right)
\end{equation*}
and, since $\pi_{\mathcal{P}}(x)=2\pi_{D}(x)= 2\lfloor\Pi_{B}(x)\rfloor= \Pi_{C}(x)+ O(1)$,
\begin{equation}
\label{asympdiamondexeq1}
\pi_{P}(x)=\frac{x}{\log x}\left(2-\sqrt{2}\cos \left(
\log x-\frac{\pi}{4}\right)\right)  +O\left(\frac{x}{\log^{2}
x}\right) ,
\end{equation}
whence both $\Pi_{C}$ and $\pi_{P}$ satisfy lower and upper Chebyshev
bounds. The zeta function $\zeta_{C}(s)$ can be explicitly computed:
$$
\log\zeta_{C}(s)= -2\log(s-1)+\log (s-1-i)+\log(s-1+i),
$$
and so
$$
\zeta_{C}(s)= \frac{(s-1)^{2}+1}{(s-1)^{2}}= 1+\frac{1}{(s-1)^{2}}\,.
$$

Now $\zeta_{C}(s)$ has right-hand zeros of order 1 located at $1\pm i$,
and so Lemma \ref{lexC-M} implies that $M_{C}(x)=\Omega(x)$. For the
discrete g-number system, we have
$$
\zeta_{\mathcal{P}}(s)=\zeta_{C}(s)\exp(G(s)),
$$
where $G(s)=\int_{1}^{\infty}x^{-s}\mathrm{d}(\Pi_{\mathcal{P}}-\Pi_{C})(x)$ is
analytic on $\Re s>1/2$ because
$$\Pi_{\mathcal{P}}(x)=\pi_{\mathcal{P}}(x)+O(\sqrt x)=\Pi_{C}(x)+O(\sqrt x).$$
Thus, the same argument yields $M_{\mathcal{P}}(x)=\Omega(x)$.

Note that $\zeta_{C}$ is the Mellin transform of the measure
$\mathrm{d}N_{C}= \delta_{1}+\log u\, \mathrm{d}u$, and therefore we
have the exact formula
$$
N_{C}(x)= x\log x- x+2, \quad x\geq 1.
$$
We now show that $N_{\mathcal{P}}$ satisfies a lower bound of a
similar type.   It is well known \cite{d-s-v,diamond1970} that
$N_{D}(x)= cx +O(x\log^{-3/2}x)$
with $c>0$. Thus $N_D(x) \ge c'x$ for some $c'>0$ and all
$x\ge 1$ and so
\begin{align*}
N_{\mathcal{P}}(x) &= \int_{1}^x \mathrm{d}N_D*\mathrm{d}N_D
= \int_{1}^x N_D(x/t)\,\mathrm{d}N_D(t)
\ge \int_{1}^x c'x/t \,\mathrm{d}N_D(t) \\
&\ge c'\!\!\int_{1}^x N_D(x/t)\,\mathrm{d}t
\ge c'^{\, 2}\!\!\int_{1}^x x/t \,\mathrm{d}t = c'^{\, 2}x \log x \ne O(x).
\end{align*}

Applying the Dirichlet hyperbola method, one can actually obtain the
sharper asymptotic estimate
\begin{equation}
\label{eqasympex2}
N_{\mathcal{P}}(x)=c^{2} x\log x+ bx+ O\left(\frac{x}{\log^{1/2}x}\right)
\end{equation}
with certain constants $b, \, c \in\mathbb{R}$. We leave the
verification of (\ref{eqasympex2})  to the reader.
\end{example}

\begin{example}  \label{ex:3C-M}
There exists a non-decreasing function $N$ on $[1, \, \infty)$ for which
\[
\int_{1}^{\infty} \frac{|N(x)-x|}{x^{\sigma + 1}}\,
\mathrm{d}x \ll (\sigma -1)^{-1/3}\, \tag{\ref{L1cond} bis}.
\]
holds for $1 < \sigma < 2$, but  $\limsup_{x \to \infty} N(x)/x = \infty$.
\end{example}

To see this, set $f(n) = e^{e^n}$ and take
\[
N(x) = \begin{cases}
x, \ x \ge 1, & x \not\in \displaystyle\bigcup_{m=1}^\infty \big[f(m),
\,e^{m/3}f(m)\big] \\
e^{n/3}f(n), & f(n) \le x \le e^{n/3}f(n).
\end{cases}
\]
Clearly, $N\big(e^{e^n}\big)/ e^{e^n} = e^{n/3} \to \infty.$

On the other hand, for each $n$ we have
\begin{align*}
\int_{f(n)}^{\,e^{n/3}f(n)} \frac{|N(x)-x|}{x^{\sigma + 1}}\,\mathrm{d}x
&<e^{n/3}f(n)\int_{f(n)}^{\,e^{n/3}f(n)} \frac{\mathrm{d}x}{x^{\sigma + 1}} \\
&< e^{n/3}f(n)\int_{f(n)}^\infty \frac{\mathrm{d}x}{x^{\sigma + 1}}
< e^{n/3} f(n)^{-(\sigma -1)}.
\end{align*}
Instead of summing $ e^{n/3} f(n)^{-(\sigma -1)}$, we calculate the
corresponding integral:
\[
\int_0^\infty e^{u/3} e^{-(\sigma -1) e^u} \mathrm{d}u
= \int_0^\infty v^{1/3} e^{-(\sigma -1) v} \frac{\mathrm{d}v}{v} =\Gamma(1/3)
(\sigma -1)^{-1/3}.
\]
Also, $\{e^{n/3} f(n)^{-(\sigma -1)}\}$ is a unimodal sequence whose maximal
term is of size at most $(3e(\sigma -1))^{-1/3} \ll (\sigma-1)^{-1/3}$.
Thus (\ref{L1cond} bis) holds.


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

Jets launched by supermassive black holes are presently being evaluated as important mechanisms for galaxy evolution, affecting multiple phases of the interstellar medium (ISM). By depositing energy in the hot and tenuous gas, jets are proposed to delay the inflow of intergalactic gas into ellipticals and moderate their growth \citep{bower06, croton06}. By depositing energy in the cold and dense gas, jets are proposed to influence the star-formation properties of disks \citep{wagner11, gaibler12}. Observational evidence that the impact of a jet upon dense clouds initiates molecular gas winds has been successfully gathered \citep[e.g.,][]{morganti05, rupke11, sakamoto14, garcia-burillo14, dasyra15, aalto16}.

In this letter, we provide observational evidence that the molecular gas accelerated or impacted by a jet is more highly excited and less optically thick than the rest of the ISM. We prove this for the galaxy IC5063, an elliptical with a molecular disk in its center. The jet launched from the black hole of IC5063 propagates through the disk, nearly parallel to its plane. Once the jet encounters clouds, it scatters gas, forming atomic and molecular winds \citep{oosterloo00,morganti07,morganti13,tadhunter14}. Using ultra-deep ESO VLT near infrared (NIR) data, we identified at least four discrete starting points of ionized and warm molecular gas winds near the jet trail \citep{dasyra15}. This constitutes the highest number of known jet-driven winds in a galaxy. The winds sweep the galaxy's inner kpc$^2$, allowing for a spatially-resolved study of the accelerated gas properties. Atacama Large Millimeter Array (ALMA) data have revealed the presence of cold molecular gas in most of these winds \citep{morganti15}.

\section{The ALMA data}
\label{sec:data}
New, band 8 data of IC5063 were taken with ALMA on May 18 2016 for the program 2015.1.00420.S (PI Combes). Forty two antennas observed the galaxy for 39 minutes on-source. Three 1.875\,GHz-bandwidth spectral windows were employed. Two of the spectral windows were partly overlapping to cover the $\pm$1000\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ range around \hbox{$\rm CO(4$-$3)$} . The third spectral window was used for duplication purposes in the $\pm$300\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ range: its role was to ensure that no significant continuum offset exists between the former windows due to flux and bandpass calibration uncertainties. The data were reduced within the Common Astronomy Software Applications (CASA) environment, version 4.6.0. The calibration script was executed as delivered by ESO. Pallas, J1924-2914, and J2056-4714 were respectively used as flux, bandpass, and phase calibrators. The flux calibration accuracy was better than 5\%, as cross-checked with J1924-2914, which is not variable as J2056-4714 is. To construct images from visibilities, we chose a 25\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ spectral bin and a 0.06\arcsec\ pixel scale. The synthesized beam size was 0.4\arcsec $\times$0.4\arcsec .

To study the gas excitation, we also used \hbox{$\rm CO(2$-$1)$}\ data from the program 2012.1.00435.S \citep{morganti15}. Thirty one antennas observed IC5063 in band 6 for 30 mins on-source. J2056-4714 was used as bandpass and phase calibrator, Neptune as flux calibrator. Four 1.875\,GHz-bandwidth spectral windows were employed, with one of them covering the $\pm$1000\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ range around \hbox{$\rm CO(2$-$1)$}. The flux calibration was performed in CASA 4.2.2. Imaging was carried out in CASA 4.6.0 (because of continuum subtraction improvements). The cube was created at a 25\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ spectral bin and at a 0.1\arcsec\ pixel scale. The synthesized beam was 0.56\arcsec $\times$0.53\arcsec .

All images were reconstructed using Briggs robust=0.5 weighting, and primary-beam corrected. No important flux losses were considered for them: the {\it uv} plane coverage was high and the area under examination  was small (3-5\arcsec\ in individual channels) compared to the primary beam in both bands (28\arcsec,\,14\arcsec ). To compare them, we re-gridded the \hbox{$\rm CO(2$-$1)$}\ cube to the pixel scale of the \hbox{$\rm CO(4$-$3)$}\ cube, and we smoothed by Gaussian convolution the \hbox{$\rm CO(4$-$3)$}\ cube to the resolution of the \hbox{$\rm CO(2$-$1)$}\ beam. In the common resolution, the noise is 0.26 mJy/beam for \hbox{$\rm CO(2$-$1)$} , and 3.4 mJy/beam for \hbox{$\rm CO(4$-$3)$} .

\section{Results: gas excitation differences in the wind and in the ambient medium}
\label{sec:results}

Images of the spatially-resolved line emission are shown in Fig.~\ref{fig:basics}, following the subtraction of the (also spatially-resolved) continuum. \hbox{$\rm CO(2$-$1)$}\ is bright along the jet propagation axis and in spirals arms, forming the S-shaped feature identified by \citet{morganti15}. Contrarily,  \hbox{$\rm CO(4$-$3)$}\ is mainly detected along the jet axis. In both lines, the maximum emission is detected near the north-west radio lobe, the south-east radio lobe, and the radio core. Gaussian modeling of the emission in both lines indicates that the projected rotational velocity reaches $\sim$200\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ within 600\,pc from the nucleus. High  \hbox{$\rm CO(4$-$3)$}\  velocity dispersion (either real or caused by the projection of several kinematic components) is seen north of the nucleus, near the north-west radio lobe, and along the biconical outflow detected by \citet{dasyra15} in NIR \,\hbox{$\rm{H_ 2}$}\ and \,\hbox{[\ion{Fe}{II}}]\ data. The latter outflow's starting point is mid-way between the nucleus and the north-west lobe, and its axis is oriented perpendicularly to the jet trail.

The total flux detected in \hbox{$\rm CO(4$-$3)$}\ is 118.0\,$\pm$\,5.0\, Jy\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$} , whereas\, that detected\, in \hbox{$\rm CO(2$-$1)$}\ is 70.2 $\pm$ 4.7 Jy\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$} .  The \hbox{$\rm CO(2$-$1)$}\ flux is 30\% lower than that measured from single-dish, APEX data. The difference is attributed to uncertainties in the continuum determination in the APEX data, in which the blue CO wing was covering most of the available band \citep{morganti13, morganti15}.

The overall  \hbox{$\rm CO(4$-$3)$} /\hbox{$\rm CO(2$-$1)$}\ flux ratio, \hbox{$\rm F_{43}/F_{21}$} , is indicative of a sub-thermally excited gas reservoir (see Section~\ref{sec:discussion}). However, when integrated within a beam's radius from the nucleus (at 20:52:02.37 -57:04:07.6) and the radio lobes (at 20:52:02.15 -57:04:06.7 and  20:52:02.54 -57:04:08.3), the flux in \hbox{$\rm CO(4$-$3)$}\  is $\sim$4 times higher than that in \hbox{$\rm CO(2$-$1)$}\ (Fig.~\ref{fig:spectra}).

To derive the fraction of the flux originating from the jet accelerated or impacted gas, we modeled the disk emission and subtracted it from the observed spectrum in each pixel. A Gaussian fitting algorithm was employed, prefered over a two-component Gaussian fitting algorithm for its lower errors. The disk velocity $V$ was tracing the peak of the line emission, never exceeding $\pm$250\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ in the inner 1\,kpc$^2$. The velocity dispersion $\sigma$ was constrained to values below 50\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$} : this is the average $\sigma$ for both lines in regions that are located more than a beam away from the nucleus and the radio lobes. In the residual cubes, the integrated flux is 12.3($\pm$2.5) Jy\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ for \hbox{$\rm CO(2$-$1)$} , and 59.7($\pm$9.2) Jy\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ for \hbox{$\rm CO(4$-$3)$} . This result indicates that $\sim$half of the reservoir probed by \hbox{$\rm CO(4$-$3)$}\ is kinematically or excitationally affected by the jet. The \hbox{$\rm F_{43}/F_{21}$}\ ratio is 4.9($\pm$)1.7 for this gas,\,and 1.0($\pm$)0.4 for the rest of the gas in the disk. For the fastest gas in the wind, which is moving with -600$<$V$<$-400\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ and which is seen on the north-west radio lobe, \hbox{$\rm F_{43}/F_{21}$}\ is 8.0($\pm$2.4).

Regions with high CO excitation are shown in a \hbox{$\rm F_{43}/F_{21}$}\ map that we constructed after collapsing the disk-subtracted cubes in the \hbox{-600}$<$V$<$400\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ range (Fig. ~\ref{fig:residual_maps}). These regions often coincide with wind starting points known from NIR \,\hbox{$\rm{H_ 2}$}\ and \,\hbox{[\ion{Fe}{II}}]\ data \citep{tadhunter14,dasyra15}. Values $\ge$8 are seen north of the north-west radio lobe, where the \,\hbox{$\rm{H_ 2}$}\ is also highly excited \citep{dasyra15}, and where diffuse 17.8 GHz emission is detected \citep{morganti07}. Values approaching 16 are seen in blobs that nearly coincide with the biconical outflow axis perpendicular to the jet \citep[see the line defined by R1 and R2 in Fig. 8 of ][]{dasyra15}. Within 200\,pc from the nucleus, the wind associated with \hbox{$\rm F_{43}/F_{21}$}$\ge$6 regions was seen in \hbox{$\rm CO(2$-$1)$}\ by \citet{morganti15}. East of the south-east radio lobe, high-velocity Fe ions and highly excited CO molecules co-exist (Fig. ~\ref{fig:residual_maps}). Counter-rotating gas is detected near the \hbox{$\rm F_{43}/F_{21}$} $\sim$7 region north of the same radio lobe in both CO lines (Fig.~\ref{fig:ratio_maps}).
All the above-mentioned regions are also seen in \hbox{$\rm F_{43}/F_{21}$}\ maps that we created from the cubes prior to the disk subtraction (Fig.~\ref{fig:ratio_maps}).

\section{Discussion: gas mass in the wind}
\label{sec:discussion}

Our analysis indicates the presence of at least two distinct molecular gas components: a low-excitation component associated with clouds in the disk, and a high-excitation component associated with gas in the jet-driven winds or near jet-cloud impact points. To be bright in \hbox{$\rm CO(4$-$3)$}, which has a critical density of 1.7$\times$10$^5$ cm$^{-3}$, both components must contain clouds of high volume density. However, while the low excitation component is optically thick, the high excitation component is (partly) optically thin: \hbox{$\rm F_{43}/F_{21}$}\ values $>$4 can only be consistent with emission from optically thin gas at high temperatures.

From the antenna temperature definition, assuming that the molecular gas is in local thermodynamic equilibrium (LTE) and that its excitation temperature \,\hbox{$T_{\rm ex}$}\ greatly exceeds that of the cosmic microwave background (see below), we obtain analytic solutions for \hbox{$\rm F_{43}/F_{21}$}\ (with fluxes in units of Jy\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$} ). For optically thick gas, \hbox{$\rm F_{43}/F_{21}$}\ = ${ \rm (\nu_{43}/\nu_{21})^3 [e^{(h \nu_{21} /kT_{ex})}-1]/[e^{(h \nu_{43} /kT_{ex})}-1]}$. An upper limit of 4 is reached at the Rayleigh-Jeans limit (as expected from the same line brightness temperature). For optically thin gas, \hbox{$\rm F_{43}/F_{21}$}\ = ${\rm (\nu_{43}/\nu_{21})^4 e^{-7h \nu_{10} /kT_{ex}}}$. An asymptotic value of 16 is obtained for  \,\hbox{$T_{\rm ex}$} $>$$>$7h$\nu_{10}$/k $\sim$ 39 K. From this relation, we find that \,\hbox{$T_{\rm ex}$} =32K, if all of the gas kinematically or excitationally affected by the jet (with \hbox{$\rm F_{43}/F_{21}$} =5) is optically thin. For the fastest gas in the wind (with \hbox{$\rm F_{43}/F_{21}$} =8),
\,\hbox{$T_{\rm ex}$}\ is 56K. Near the biconical outflow base, where \hbox{$\rm F_{43}/F_{21}$}\ is 11, \,\hbox{$T_{\rm ex}$}\ is 100K. The temperatures of the gas in the wind can thus be quite high for LTE, and they can be even higher for non-LTE.

It is unlikely that this result reflects infrared (IR) pumping of the populations of quantum states with high rotational number J. The molecules should see an IR source of $>$160K with a large filling factor \citep{carroll81} in all pertinent regions of Figs.~\ref{fig:residual_maps} and \ref{fig:ratio_maps}. Infrared pumping has been deemed unimportant for HCN in the wind of Mrk231 \citep{aalto12}. The high \,\hbox{$T_{\rm ex}$}\ is instead related to the creation of dense fragments once the jet impinges on molecular clouds. The contribution of the jet to the dissipation and heating of the clouds could thus have a negative impact on the overall star formation of IC5063. Still, it doesn't rule out a local enhancement of star formation in compressed clouds.

Given the above bimodality, we determined the gas mass in the wind and in the disk using the appropriate CO-luminosity-to-\,\hbox{$\rm{H_ 2}$} -mass conversion factor, $\alpha_{CO}$, for optically thin and thick emission, respectively. For the optically thick gas, we use Eq. 3 of \citet{solomon97} with $\alpha_{CO}$=4.6\,\hbox{$M_{\odot}$} /(K\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\,pc$^2$), as the far infrared luminosity of IC5063 is 2$\times$10$^{10}$\,\hbox{$L_{\odot}$} .  For \hbox{$\rm F_{21}/F_{10}$}=4, a mass of 4.0($\pm$0.5)$\times$10$^8$\,\hbox{$M_{\odot}$}\ is obtained for the disk. For the optically thin gas, \hbox{$\rm F_{21}/F_{10}$} = $ {\rm (\nu_{21}/\nu_{10})^4 e^{-2h \nu_{10} /kT_{ex}}}$ and $\alpha_{CO}=0.26\,(\,\hbox{$T_{\rm ex}$} /30)\,e^{5.53/\,\hbox{$T_{\rm ex}$}-0.184}$ \,\hbox{$M_{\odot}$} / (K\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\,pc$^2$), for a Galactic CO/\,\hbox{$\rm{H_ 2}$}\ abundance \citep{bolatto13}. \hbox{For $\rm F_{21}$=12.3($\pm$2.5) Jy\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ and \hbox{$\rm F_{43}/F_{21}$}=5, then $\alpha_{CO}$} \hbox{is 0.27\,\hbox{$M_{\odot}$}/\,(K\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\,pc$^2$) and the wind mass  is 1.7($\pm$0.3)}  \hbox{$\times$10$^6$\,\hbox{$M_{\odot}$} . For \hbox{$\rm F_{43}/F_{21}$}=8, $\alpha_{CO}$ is 0.44\,\hbox{$M_{\odot}$} /\,(K\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\,pc$^2$)} and the wind mass  is 2.4($\pm$0.5)$\times$10$^6$\,\hbox{$M_{\odot}$} . An order-of-magnitude difference in $\alpha_{CO}$ between the disk and the wind was already assumed by \citet{morganti15}, but the use of an optically-thick \hbox{$\rm F_{21}/F_{10}$} value led to a wind mass of 1.9$\times$10$^7$\,\hbox{$M_{\odot}$} . Our current knowledge that the wind is partly optically thin leads us to revise its molecular content down: 2$\times$10$^6$\,\hbox{$M_{\odot}$}\ is its lower mass limit, if all of its gas is optically thin. This is comparable to the atomic wind mass \citep[3.6$\times$10$^6$\,\hbox{$M_{\odot}$} ; ][]{morganti07}. The actual molecular wind mass will be derived when all data, including new CO (1-0) and (3-2) observations, are analyzed assuming non-LTE conditions: when both high-temperature gas and low-temperature gas (contributing to low J lines) are used to fit the CO spectral line energy distribution (SLED). Still, our simplified LTE analysis proves that molecular wind masses can be overestimated \citep[contrary to what is known from Mrk231;][]{cicone12}.

\section{Summary}
\label{sec:conclusions}

We obtained new ALMA \hbox{$\rm CO(4$-$3)$}\ data and combined them with previous \hbox{$\rm CO(2$-$1)$}\ data to study the molecular gas excitation in the jet-driven winds of IC5063. We found that \hbox{$\rm CO(4$-$3)$}\ is mainly detected along the jet. Half of its emission comes from gas that was impacted by the jet, while only one sixth of the \hbox{$\rm CO(2$-$1)$}\ emission comes from the same gas. In this gas, the average \hbox{$\rm F_{43}/F_{21}$}\ ratio is 5. In regions with a strong jet-cloud interaction, the ratio reaches 12 in the original cubes and approaches 16 in the disk-subtracted cubes. Contrarily, the same ratio for the rest of the gas in the disk is 1. We thus detect at least two gas components: one optically thick related to the disk, and another optically thin related to the winds. Because the outflowing molecular gas is partly optically thin, its mass is lower than that previously thought, as low as 2$\times$10$^6$\,\hbox{$M_{\odot}$} . This is comparable to the mass of the outflowing HI.  Excitation temperatures of 30-100K are common for the outflowing molecules. The skewing of the CO SLED toward higher J numbers can facilitate the detection of winds at high z.

\begin{acknowledgements}
K. D. acknowledges support by the European Commission through a Marie Curie Intra-European Fellowship (PIEF-GA-2013-627195;  'BHs shaping galaxies') awarded under the 7th Framework Programme. R. M. gratefully acknowledges support from the European Research Council under the European Union's 7th Framework Programme, ERC Advanced Grant RADIOLIFE-320745.
This paper makes use of the ALMA data 2015.1.00420.S and 2012.1.00435.S.  ALMA is a partnership of ESO (representing its member states), NSF (USA), NINS (Japan), together with NRC (Canada), NSC and
ASIAA (Taiwan) and KASI (Republic of Korea), in cooperation with the Republic of Chile.  The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ.
\end{acknowledgements}

\begin{figure*}[h!]
\begin{center}
\includegraphics[width=19cm]{basics.eps} \\
\caption{ {\it Left:} images of CO (4-3), (2-1), and their nearby continua, extracted from the optimal-resolution cubes (i.e., 0.4\arcsec$\times$0.4\arcsec\ beam in band 8 and 0.56\arcsec$\times$0.53\arcsec\ beam in band 6). The line fluxes in mJy/beam were averaged over the -100$<$V$<$100\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ range. Crosses indicate the locations of the nucleus and the two radio lobes, as seen in the \hbox{$\rm CO(2$-$1)$}\ continuum (see Sect.~\ref{sec:data} for the coordinates). {\it Right:} maps of the amplitude, velocity, and velocity dispersion of a Gaussian line-profile fit. All panels are centered at the nucleus, at 20:52:02.37 -57:04:07.6. They have a field of view (FOV) of 7\arcsec$\times$7\arcsec .}
\label{fig:basics}
\end{center}
\end{figure*}

\begin{figure*}[h!]
\begin{center}
\includegraphics[width=0.78\columnwidth]{spectra.eps}
\caption{ Integrated \hbox{$\rm CO(4$-$3)$}\ and \hbox{$\rm CO(2$-$1)$}\ spectra, extracted from the cubes at the common spatial resolution. The \hbox{$\rm CO(2$-$1)$}\ flux is multiplied by a factor of four to facilitate comparison with the \hbox{$\rm CO(2$-$1)$}\ flux.  {\it Upper left:} Emission integrated over the whole disk, that is within an ellipse with a 12\arcsec\ major axis, 4\arcsec\ minor axis, and \hbox{-65$^{\rm o}$} position angle.  {\it Upper right and bottom:} Emission integrated within a beam (i.e., within a radius of 0.28\arcsec ) from the nucleus, and from the north and south radio lobes. }
\label{fig:spectra}
\end{center}
\end{figure*}

\begin{figure*}[h!]
\begin{center}
\includegraphics[width=0.9\textwidth]{model_cube_ratio.eps} \\
\caption{ {\it Left, middle:} residual images from the subtraction of a disk model from the data (see Section~\ref{sec:results}). Their fluxes, integrated in the -600$<$V$<$400\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ range, are in units of Jy\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ per pixel. The common beam is 0.56\arcsec$\times$0.53\arcsec . Pixels with signal-to-noise (S/N)\,$<$5 have been masked. {\it Right:} \hbox{$\rm F_{43}/F_{21}$}\ of the residual images. Contours of the outflowing-to-ambient \,\hbox{[\ion{Fe}{II}}]\ emission from \citet{dasyra15} are overplotted, in steps of 0.2 starting from 0.1 for the fast winds near the lobes, and steps of 0.33 starting from 0.33 for the intermediate-velocity (component) of the other winds. The FOV and center are as in Fig.~\ref{fig:basics}.}
\label{fig:residual_maps}
\end{center}
\end{figure*}

\begin{figure*}[h!]
\begin{center}
\includegraphics[width=0.9\textwidth]{cube_ratios.eps} \\
\caption{ CO (4-3) and (2-1) images and \hbox{$\rm F_{43}/F_{21}$}\ ratio, integrated in 200\,\hbox{\hbox{km}\,\hbox{s}$^{-1}$}\ bins. All panels were constructed directly from the data, and, unlike Fig.~\ref{fig:residual_maps}, they contain the disk emission. The image units, beam, FOV, and center are as in Fig.~\ref{fig:residual_maps}. Pixels with S/N$<$3 have been masked. The contours in the last panel indicate where the \hbox{$\rm CO(4$-$3)$}\ emission reaches S/N of 5,\,10,\,and 20.}
\label{fig:ratio_maps}
\end{center}
\end{figure*}