\section{Introduction}

The group ${\rm{SL}}(n,\mathbb Z)$ acts faithfully by linear transformations on $\mathbb R^n$ and hence on the sphere at infinity $\S^{n-1}$,
but if $n\ge 3$ then it cannot act on  lower-dimensional spheres or Euclidean spaces.
Indeed,  in \cite{BV}, Bridson and Vogtmann\begin{footnote}{The proof of a
similar theorem announced earlier  \cite{Parwani} is not valid.}\end{footnote} proved that if $n\ge 3$ and $d<n$, then
${\rm{SL}}(n,\mathbb Z)$ cannot act non-trivially by homeomorphisms on any contractible manifold of dimension $d$, nor on any homology sphere
of dimension $d-1$. This is an immediate corollary of a more general theorem proved in \cite{BV}: if $n\ge 3$ then
${\rm{SAut}}(F_n)$, the unique subgroup of index $2$ in the automorphism group of the free group $F_n$, cannot act
non-trivially by homeomorphisms
 on any $\mathbb Z/2$-acyclic $\mathbb Z/2$-homology manifold of dimension less than $n$, nor any
 generalized $\mathbb Z/2$-homology sphere of dimension less than
$n-1$.
The proof revolved around the elementary abelian $2$-groups in ${\rm{SAut}}(F_n)$. Using Smith theory, one
argues that one of a small number of involutions have to act trivially; a detailed analysis of the quotients
of ${\rm{SAut}}(F_n)$ by these involutions then allows one to conclude that the entire group has to act trivially.

Smith theory applies for any prime $p$. When $p$ is odd there is a stronger restriction on the dimension of a sphere on which an elementary abelian $p$-subgroup can act.  If $n$ is even, then ${\rm{SAut}}(F_n)$ contains an elementary
3-group of rank $n/2$, and by exploiting this one can give a considerably shorter proof of the fact that
${\rm{SAut}}(F_{2m})$  (and hence ${\rm{SL}}(2m,\mathbb Z)$) cannot act on an acyclic manifold of dimension less than $2m$
or a $\mathbb Z/3$-homology sphere of dimension less than $2m-1$ (see \cite{BV}, Theorem 4.9).
The purpose of the
present note is to point out that the strategy of this shorter proof, combined with  ideas of
Zimmermann \cite{Z1},
can be used to establish similar theorems for many  arithmetic groups. In most cases, this will give only a lower bound on the dimension
of spheres and acyclic spaces where the
group can act. But in good cases, where one has a suitably large elementary abelian $p$-subgroup for some $p$,
the bounds that one obtains can be sharp.

One such case is that of the symplectic groups. In this case,  we prove that for $n>2$ any action of ${\rm Sp}(2n,\mathbb Z)$ on  a generalized $\mathbb Z/3$-homology sphere of dimension  less than $2n-1$ or a $\mathbb Z/3$-acyclic $\mathbb Z/3$-homology manifold of dimension less than  $2n$ must be trivial. If  $n=2$ any such action factors through the abelianization of ${\rm{Sp}}(4,\mathbb Z)$, which is $\mathbb Z/2$.
The dimension bound is sharp,
since ${\rm{Sp}}(2n,\mathbb Z)$ does act (linearly) on $\mathbb R^{2n}$ and $\S^{2n-1}$.
This result, with a weaker conclusion for $n=2$, was obtained independently by Zimmermann (see \cite{Z}).

Our investigations are motivated by the following general question:
\begin{question} Given an affine algebraic group scheme $\mathbb G$ defined over $\mathbb Z$ and a ring $\mathfrak o_S$ of $S$-integers  in a number field $k,$ what is the smallest dimensional sphere (or Euclidean space) on which $\mathbb G(\mathfrak o_S)$ can act? In particular, are there examples of actions below the
first dimension in which  $\mathbb G(\mathfrak o_S)$ has a non-trivial linear representation over $\mathbb R$?
\end{question}

The bounds that we establish depend on the torsion in the groups $\mathbb G(\mathfrak o_S)$ and are summarized in Theorem~\ref{t:perfect}.
 In Section \ref{s:sp} we show how everything works in the concrete setting of the symplectic group ${\rm{Sp}}(2n,\mathbb Z)$.  In Section \ref{s:lemma} we point out that the basic features of the argument work for much more general arithmetic  groups, and a slightly weaker hypothesis may be used.  We emphasize here, however, that this general argument does not eliminate the need for special arguments in low-dimensional cases, even for the symplectic group.

In Section \ref{s:lattices} we discuss conditions on linear algebraic groups $G$ and rings of $S$-integers $\mathfrak o_S$ that are sufficient to ensure that
the lattices $\Gamma=G(\mathfrak o_S)$ have the features needed for our general argument. The lower bounds on the dimension in which such groups $\Gamma$
have interesting actions depend on the existence of elementary $p$-subgroups; the larger
the rank of such a subgroup, the better the bounds one obtains.
We provide some elementary examples of such subgroups, but we leave open the question of how
to identify the largest such subgroup in general.

\medskip

\noindent{\em{Acknowledgements.}} We thank Alex Lubotzky, Gopal Prasad and Alan Reid  for their helpful comments concerning the
material in Section 4. Most particularly, we thank Dan Segal for his notes on this material, from which we borrowed heavily. We also
thank the Institute Mittag-Leffler (Djursholm, Sweden) for its hospitality during the preparation of this manuscript.

Tragically, the second author did not survive to see this project completed. He is sorely missed for many reasons.
Any deficiencies in the final version of this paper are the responsibility of the first and third authors alone.

\section{The integral symplectic group}\label{s:sp}
The symplectic group ${\rm{Sp}}(2n,\mathbb Z)$ is the subgroup of ${\rm{GL}}(2n,\mathbb Z)$  consisting of matrices which preserve the standard symplectic form.  With respect to a symplectic basis $$\{a_1,\ldots,a_n,b_1,\ldots, b_n\}$$
for $V=\mathbb Z^{2n}$, the form is represented by the matrix
$$
J=\begin{pmatrix}
      0 & I \\
      -I&0\\
 \end{pmatrix}
$$
 and an integer matrix $M$ lies in   ${\rm{Sp}}(2n,\mathbb Z)$ if and only if  $^tMJM=J$.

 There is a subgroup $T$ of ${\rm{Sp}}(2n,\mathbb Z)$ isomorphic to $(\mathbb Z/3)^n$, generated by the symplectic transformations
 $$R_i\colon
 \begin{cases}
a_i\mapsto -b_i  \\
  b_i\mapsto a_i-b_i\\
  a_k\mapsto a_k \hbox{ for } k\neq i\\
  b_k\mapsto b_k \hbox{ for } k\neq i\\
  \end{cases}
  $$
  If $V_i$ is the subspace with basis $\{a_i,b_i\}$, then $R_i$ is the identity on $V^\perp_i\subset\mathbb R^{2n}$  and the matrix of $R_i$ restricted to  $V_i$   is
$
\begin{pmatrix}
      0& 1\\
      -1&{-1}\\
 \end{pmatrix}.
$
We remark that $T$ is also a subgroup of ${\rm{SL}}(2n,\mathbb Z),$ and is the same subgroup as was  used in \cite{BV}, Lemma  3.2.

We now recall the statement of Theorem 4.7 of  \cite{BV}  for odd primes $p$, which is proved using Smith Theory.

\begin{theorem}\label{Smith}  Let $p$ be an odd prime.   If $m < 2d-1$, then $(\mathbb Z/p)^d$ cannot act effectively
by homeomorphisms on a generalized  $\mathbb Z/p$-homology sphere of dimension  $m$ or a $\mathbb Z/p$-acyclic $\mathbb Z/p$-homology manifold of dimension $m+1$.
\end{theorem}

For $T\cong (\mathbb Z/3)^n$ this says:

\begin{corollary}\label{trivial} Whenever $T$ acts by homeomorphisms on a generalized $\mathbb Z/3$-homology sphere of dimension less than  $2n-1$ or a
$\mathbb Z/3$-acyclic $\mathbb Z/3$-homology manifold of dimension less than $2n$,  some non-trivial element of $T$ acts trivially.
\end{corollary}

We shall use the structure of ${\rm{Sp}}(2n,\mathbb Z)$ to deduce Theorem \ref{t:sp} from Corollary \ref{trivial}.
The case $n=2$ is special: for $n\ge 3$ the group
${\rm{Sp}}(2n,\mathbb Z)$ is perfect, but ${\rm{Sp}}(4,\mathbb Z)$ maps onto ${\rm{Sp}}(4,2)$, which is
isomorphic to the symmetric group $\Sigma_6$ and hence maps onto $\mathbb Z/2$, which acts non-trivially on the
line and the $0$-sphere.

\begin{theorem}\label{t:sp} Let $X$ be a generalized $\mathbb Z/3$-homology sphere of dimension  less than $2n-1$ or a $\mathbb Z/3$-acyclic $\mathbb Z/3$-homology manifold of dimension less than  $2n$.
If $n\geq 3$, then any action of ${\rm{Sp}}(2n,\mathbb Z)$ by homeomorphisms on $X$ is trivial. If $n=2$, then every
action of ${\rm{Sp}}(2n,\mathbb Z)$ on $X$ factors through the abelianization of ${\rm{Sp}}(4,\mathbb Z)$, which is  $\mathbb Z/2$.
\end{theorem}

\begin{proof}
Let  $\rho\colon{\rm{Sp}}(2n,\mathbb Z)\to{\rm{Homeo}}(X)$ be an action of ${\rm{Sp}}(2n,\mathbb Z)$ on $X$ and
let  $T\cong (\mathbb Z/3)^n$ be the subgroup defined above. By Corollary~\ref{trivial}, there is a non-trivial element $t\in T$ in the kernel of $\rho$.

The center of ${\rm{Sp}}(2n,\mathbb Z)$ is $\pm I$, and is not contained in $T$.  By the Margulis normal subgroup theorem, any normal subgroup of ${\rm{Sp}}(2n,\mathbb Z)$ is either central or has finite index.  Since the kernel of $\rho$ contains a non-central element, it must have finite index.
${\rm{Sp}}(2n,\mathbb Z)$ has the congruence subgroup property, so  the kernel of  $\rho$   contains a principal congruence subgroup $\Gamma_m$ for some $m$.
Thus our action factors through the quotient ${\rm{Sp}}(2n,\mathbb Z)/\Gamma_m\cong {\rm{Sp}}(2n,\mathbb Z/m)$, and it suffices to show that any action of ${\rm{Sp}}(2n,\mathbb Z/m)$ is trivial.

Now ${\rm{Sp}}(2n,\mathbb Z/m)\cong \prod_j{\rm{Sp}}(2n,\mathbb Z/{p_i^{e_i}})$, where $p_i$ runs over the prime divisors of $m$ (see, e.g., \cite{N} for an elementary proof of this).  If an action of ${\rm{Sp}}(2n,\mathbb Z/m)$
 is non-trivial, then its restriction to some ${\rm{Sp}}(2n,\mathbb Z/{p_i^{e_i}})$ must be non-trivial. So it is enough to prove
the theorem with
${\rm{Sp}}(2n,\mathbb Z)$ replaced by ${\rm{Sp}}(2n,\mathbb Z/{p_i^{e_i}})$.

So now we assume that $\rho$ is an action of ${\rm{Sp}}(2n,\mathbb Z/p^e)$  on $X$, and we let
$Q\leq{\rm{Homeo}}(X)$ denote its image.  $T$ maps injectively to ${\rm{Sp}}(m,\mathbb Z/p^e)$, so we
identify it with its image. As above,  some non-trivial $t\in T$  acts trivially on $X$.

The kernel of the natural map from ${\rm{Sp}}(2n,\mathbb Z/p^e)$ onto ${\rm{Sp}}(2n,\mathbb Z/p)$ is a $p$-group.
Except in the case $n=p=2$ the quotient ${\rm{PSp}}(2n,p)={\rm{Sp}}(2n,\mathbb Z/p)/\{\pm I\}$ is simple, so the composition factors of ${\rm{Sp}}(2n,\mathbb Z/p^e)$ are ${\rm{PSp}}(2n,p)$, $\mathbb Z/2$ and $\mathbb Z/p$.  Furthermore (except for ${\rm Sp}(4,2)$) ${\rm{Sp}}(2n,\mathbb Z/p^e)$  is perfect; thus if $Q$ is non-trivial it must map onto ${\rm PSp}(2n,p)$.  But any two maps ${\rm{Sp}}(2n,\mathbb Z/p^e)\twoheadrightarrow {\rm{PSp}}(2n,p)$ have the same kernel, which does {\it not} contain $t$.   This contradicts the fact that the image of $t$ in $Q$ is trivial.

If $n=p=2$, then we must deal with the possibility that the image of
$\rho$ is a non-trivial $2$-group. A convenient way to do this is to note that the following
element of order $5$ in ${\rm{Sp}}(4,\mathbb Z)$ has normal closure of index $2$:
$$
\begin{pmatrix}
      0& 1& 0& 0\\
      -1& 0& 1& {-1}\\
      0& 1& 0& 1\\
     0& 0& -1& 1\\
 \end{pmatrix}
$$
It follows that the only non-trivial
finite 2-group onto which ${\rm{Sp}}(4,\mathbb Z)$ can map is $\mathbb Z/2$.

(This matrix is the image of an element $\mu$ of order $10$ in the mapping class group ${\rm{Mod}}_2$ of a genus 2 surface, under the natural map ${\rm{Mod}}_2\to {\rm{Sp}}(4,\mathbb Z)$.  We refer the reader to \cite{HK} for a  geometric argument  that  the normal closure of  $\mu$ in ${\rm{Mod}}_2$ has index $2$.)
\end{proof}

\section{A General Argument}\label{s:lemma}

The argument for ${\rm{Sp}}(2n,\mathbb Z)$ used several features of this group which are common in higher-rank arithmetic groups.   ${\rm{Sp}}(2n)$ is an example of a group scheme defined over $\mathbb Z$, and it  is convenient to use this language to phrase a general question about actions of such groups on spheres and Euclidean spaces.

Recall that an {\it affine algebraic group scheme $\mathbb G$ defined over $\mathbb Z$} is a functor from the category of rings to groups which is {\it represented} by some finitely-generated ring $A$; this means that  there is a natural bijection between $\mathbb G(R)$ and ${\rm{Hom}}(A,R)$ for any ring $R$.
As a simple example, consider the functor $\mathbb G={\rm{SL}}_2$; it is an affine group scheme over $\mathbb Z$ represented by the ring   $A=\mathbb Z[x_1,x_2,x_3,x_4]/\langle x_1x_4-x_2x_3-1\rangle$.  For a  very readable introduction to affine group schemes, see \cite{Waterhouse}.

\begin{question} Given an affine algebraic group scheme $\mathbb G$ defined over $\mathbb Z$ and a ring $\mathfrak o_S$ of $S$-integers  in a number field $k,$ what is the smallest dimensional sphere (or Euclidean space) on which $\mathbb G(\mathfrak o_S)$ can act? In particular, are there examples of actions below the
first dimension in which  $\mathbb G(\mathfrak o_S)$ has a non-trivial linear representation over $\mathbb R$?
\end{question}

The answer to this question will depend on both the group scheme $\mathbb G$ and on  $\mathfrak o_S$.  In our answer for $\mathbb G={\rm{Sp}}(2n)$ and $\mathfrak o_S=\mathbb Z$ we used Smith theory to find a non-central element in ${\rm{Sp}}(2n,\mathbb Z)$ which acted trivially; this allowed us to reduce the problem to  actions of ${\rm{Sp}}(2n,\mathbb Z/p^e)$.  We then began over, with a Smith theory argument for ${\rm{Sp}}(2n,\mathbb Z/p^e)$.  Instead of starting over, we could have hypothesized the existence of one element $t\in{\rm{Sp}}(2n,\mathbb Z)$ which acts trivially and which projects to a non-central element in each ${\rm{Sp}}(2n,\mathbb Z/p^e)$.  A slightly more complicated argument, which we give below, can then be used to show that  $t$ normally generates ${\rm{Sp}}(2n,\mathbb Z)$ so the whole action is trivial.  The advantage of this  new  argument is that it applies whenever   $\Gamma=\mathbb G(\mathfrak o_S)$   enjoys the following properties:

\begin{enumerate}

\item {\bf  Normal Subgroup Property.}
Every normal subgroup of $\Gamma$ is either central or has finite index

\item  {\bf Weak Congruence Subgroup Property.}  For every finite-index normal subgroup $K<\Gamma$,
the pre-image in $\Gamma$ of the centre of $\Gamma/K$
 contains a principal congruence subgroup $\Gamma(\mathfrak a)=\ker (\Gamma\to \mathbb G(\mathfrak o_S/\mathfrak a))$.

\item  {\bf Very Strong Approximation.}  The map $\Gamma\to \mathbb G(\mathfrak o_S/\mathfrak a)$ is surjective for all ideals $\mathfrak a$.

\item {\bf Quasi-Solvable Quotients.}
For every prime ideal $\P$ in $\mathfrak o_S$ the quotient of $\mathbb G(\mathfrak o_S/\P)$ by its centre is centreless with only solvable proper quotients
(e.g. simple).

\end{enumerate}

\begin{proposition}\label{Normal}  Let $\Gamma$ satisfy properties (1)-(4) above and
suppose that $\Gamma$ contains an element  $t$ whose image in $\mathbb G(\mathfrak o_S/\P)$ is non-central for every prime ideal $\P$.
Then the quotient of $\Gamma$ by the normal closure of $t$ is finite and solvable. In particular, if $\Gamma$
is perfect then the normal closure of $t$ is equal to $\Gamma$.
\end{proposition}

\begin{proof} Let $K$ be the subgroup of $\Gamma$ normally generated by $t$,
let  $N$ be the pre-image in $\Gamma$ of the centre of $\Gamma/K$, and let $f\colon \Gamma\twoheadrightarrow  \Gamma/N := Q$
be the quotient map.

Since $N$ contains the non-central element $t$, it has finite index, by $(1)$.
 The Congruence Subgroup Property (2) tells us that  $N$
 must contain a congruence subgroup $\G({\mathfrak a})=\ker (\Gamma\to \mathbb G(\mathfrak o_S/\mathfrak a))$, and   (3) tells us that $\Gamma\to \mathbb G(\mathfrak o_S/\mathfrak a)$ is onto. Thus $f$ induces an epimorphism
$$\overline f\colon  \mathbb G(\mathfrak o_S/\mathfrak a)\twoheadrightarrow Q.$$

Since $\mathfrak o_S$ is a Dedekind domain, $\mathfrak a$ factors uniquely as a product of powers of prime ideals,
say $\mathfrak a=\prod \P_i^{e_i}$. By the Chinese Remainder Theorem,
 $\mathbb G(\mathfrak o_S/\mathfrak a)= \prod \mathbb G(\mathfrak o_S/\P_i^{e_i})$. We will be done if we can argue that the restriction of $\overline f$
 to each $\mathbb G(\mathfrak o_S/\P_i^{e_i})$ has solvable image.

An elementary calculation shows that the kernel  of the natural map
$\mathbb G(\mathfrak o_S/\P^e)\to\mathbb G(\mathfrak o_S/\P)$ is a $p$-group, where $p$ is the characteristic of $\mathfrak o_S/\P$. Thus the kernel of $$\prod \mathbb G(\mathfrak o_S/\P_i^{e_i})
\to \prod \mathbb G(\mathfrak o_S/\P_i)$$ is a nilpotent group, which we denote $P$.

Now consider the following commutative diagram of epimorphisms:
$$
\begin{matrix}
\prod\mathbb G(\mathfrak o_S/\P_i^{e_i})&\xrightarrow{\overline{f}}&Q\\
\downarrow &&\downarrow\\
\prod \mathbb G(\mathfrak o_S/\P_i)& \xrightarrow{} & Q/\overline{f}(P)
\end{matrix}
$$
 Let $(t_1,\dots,t_n)$ denote the image of $t$ in $\prod \mathbb G(\mathfrak o_S/\P_i)$. By hypothesis, each $t_i$ is non-central
 in $\mathbb G(\mathfrak o_S/\P_i)$.
On the other hand, the image of $t$ in $Q$ is trivial, by hypothesis, so the image of each $t_i$ in
$Q/\overline{f}(P)$ is central by the virtue of the following elementary observation:

\smallskip
{\em{Let $\phi:G_1\times G_2\to A$ be an epimorphism of groups. If $\phi(g_1,g_2)=1$
then $\phi(g_1,1)$ and $\phi(1,g_2)$ are central in $A$.}}

\smallskip

Property (4) now implies that  the image of $\mathbb G(\mathfrak o_S/\P_i)$ in $Q/\overline{f}(P)$ is  solvable,
and hence $Q/\overline{f}(P)$, a product of commuting solvable groups, is solvable.
Since $\overline{f}(P)$ is nilpotent, we conclude that $Q$ is solvable and hence so is
its central extension $\Gamma/K$, as required.
\end{proof}

We combine this proposition with Theorem~\ref{Smith} to obtain restrictions on the possible actions of $\Gamma$.

\begin{theorem}\label{t:perfect}
Let $p$ be an odd prime,   let  $\Gamma=\mathbb G(\mathfrak o_S)$, as above, be a perfect
group with properties (1)-(4), and consider an action
of $\Gamma$ on a generalized $m$-dimensional $\mathbb Z/p$-sphere or  $(m + 1)$-dimensional   $\mathbb Z/p$-acyclic homology manifold over $\mathbb Z/p$. Suppose that $\Gamma$ contains an elementary abelian $p$-group $T$ of rank $\lfloor m/2\rfloor+1$ whose projection to $\mathbb G(\mathfrak o_S/\P)$ intersects the centre trivially for each
prime ideal $\P$ in $\mathfrak o_S$.  Then the action is trivial.
\end{theorem}

\begin{proof}  Theorem~\ref{Smith}  gives us an element $t\in T\subset\Gamma$ which acts trivially.  Then Proposition~\ref{Normal} allows us to conclude that the entire group $\Gamma$ acts trivially.
\end{proof}

The subgroup $T$ used in our previous proof for ${\rm{Sp}}(2n,\mathbb Z)$ satisfies the hypothesis of Theorem~\ref{t:perfect}, so this gives a slightly different proof for
${\rm{Sp}}(2n,\mathbb Z)$ when $n\ge 3$.

\section{Arithmetic groups}\label{s:lattices}

 Let $G$ be an algebraic group defined over $\mathbb Q$.  We fix an embedding $G\hookrightarrow {\rm{GL}}_n$.  This means that we identify $G$ with a set of invertible matrices whose entries satisfy specific polynomial equations with coefficients in $\mathbb Q.$
 These equations make sense for any field $k$ containing $\mathbb Q$, so we can define $G(k)$ to be ${\rm{GL}}(n,k)\cap G$.  If $k$ is an algebraic number field with ring of integers $\mathfrak o$ and $S$ a set of valuations,  we may then define $G(\mathfrak o_S)$ to be ${\rm{GL}}(n,\mathfrak o_S)\cap G(k).$

 In this section we discuss conditions on  $G, k$ and $S$  which will guarantee that $\Gamma=G(\mathfrak o_S)$ has the properties we need for Proposition~\ref{Normal}.    We assume throughout that $G$ is connected (in the Zariski topology), simply connected (has no proper \'etale covers), and absolutely simple (i.e. simple when considered as an algebraic group over the algebraic closure of the defining field). Examples of such groups include the classical special linear and symplectic groups, as well as commutator subgroups of the orthogonal groups.  We also assume that  $\Gamma$ is infinite throughout.

\subsection{Normal subgroups} The first property in the list   from Section~\ref{s:lemma}  is the Normal Subgroup Property, which was established for a large class of algebraic groups by Margulis \cite{Margulis}.  In particular, under the conditions fixed above, the Normal Subgroup Property holds for $\Gamma$ as long as the dimension of a maximal $k$-split torus in $G$ is at least 2 ($rank_k(G)\geq 2$) (\cite{Margulis},  quoted in \cite{PR} as Theorem 9.9).

\subsection{Congruence Subgroups}
Sufficient conditions for ensuring Property (2) (the Congruence Subgroup property) are given at length in \cite{PR}, chapter 9.  Our standing assumptions on $G$ are sufficient together with:
\begin{itemize}
\item $G$ has $k$-rank at least 2,
\item if $G$ is one of $B_2$ or $G_2$ then $S$ contains all prime divisors of $2$ with residue degree 1.
\end{itemize}

 \subsection{Strong Approximation}
The third property on our list  is the most problematic.     For $G={\rm{SL}}_n$ the statement is classical and can be found, for example, in Bourbaki (\cite{Bourbaki}, Chapter VII, section 2, part 4 on Dedekind domains).  The proof uses the fact that ${\rm{SL}}(n,\mathfrak o_S)$ and ${\rm{SL}}(n,\mathfrak o_S/\mathfrak a)$  are generated by elementary matrices $I+\lambda E_{i,j}$ together with the Chinese Remainder Theorem.

For most linear algebraic groups $G$ one cannot expect $\Gamma\to G(\mathfrak o_S/\mathfrak a)$ to be surjective for all $\mathfrak a$.  One reason is that the equations defining $G$  may change drastically when reduced modulo a prime, or even become completely trivial; in fact we should write $G^{(\mathfrak a)}(\mathfrak o_S/\mathfrak a)$ to indicate that we are considering the defining equations for $G$ modulo $\mathfrak a$ in the image.  However, in our proof we cannot afford to exclude any primes because we have no control over which ideal $\mathfrak a$ will be given to us by the  Congruence Subgroup Property.    It turns out that it is enough to  assume that
\begin{itemize}
\item For every prime $\P$ in $\mathfrak o_S$, the map $G(\mathfrak o_\P)\to G^{(\P)}(\mathfrak o/\P)$ is surjective.
\end{itemize}
(Here $\mathfrak o_\P$ is the localization of $\mathfrak o$ at the prime $\P$, i.e. the ring obtained from $\mathfrak o$ by inverting all primes except $\P$, and $\mathfrak o/\P=\mathfrak o_S/\P=\mathfrak o_\P/\P$.)  This assumption is satisfied, for example, by all of the Chevalley groups (see \cite{Carter}).

Under this assumption, Strong Approximation implies the needed surjectivity.  This argument goes as follows. The Strong Approximation as stated in \cite{PR}, Theorem 7.12, can be rephrased as:

\begin{theorem}\label{thm:SA}   Let $\P_{1},\ldots,\P_{r}$ be primes not in $S$.
Let $e_{i}\in\mathbb{N}$ and put
\[
K_{i}=\left\{  g\in G(\mathfrak o_{\P_{i}})\mid g\equiv I_{n}~(\operatorname{mod}\P_{i}^{e_{i}})\right\}  .
\]
Embed $\Gamma$ diagonally in $G(\mathfrak o_{\P_{1}})\times
\cdots\times G(\mathfrak o_{\P_{r}})$. Then
\[
\Gamma\cdot(K_{1}\times\cdots\times K_{r})=G(\mathfrak o_{\P_{1} })\times\cdots\times G(\mathfrak o_{\P_{r}}).
\]
\end{theorem}

When $G={\rm{SL}}_n$, $\mathfrak o=\mathbb Z$, $S=\emptyset$ and we have only one prime $p$,
the theorem reduces to the statement ${\rm{SL}}(n,\mathbb Z)\cdot K={\rm{SL}}(n,\mathbb Z_p)$, i.e. every determinant 1 matrix with entries in $\mathbb Z_p$ can be written as a product of a matrix with entries in $\mathbb Z$ times a matrix  congruent to the identity mod $p^e$ for any natural number $e$.

Now let $\mathfrak{a}=\prod\P_{i}^{e_{i}}$ be an arbitrary ideal of $\mathfrak o_S$. It follows that the image of $\Gamma$ in ${\rm{GL}}_n(\mathfrak o/\mathfrak{a})$ is the same as the image of $G(\mathfrak o_{\P_{1}})\times\cdots\times G(\mathfrak o_{\P_{r}})$.   Using the Chinese Remainder theorem and the assumption that the maps $G(\mathfrak o_\P)\to G^{(\P)}(\mathfrak o/\P)$ are surjective, one shows that this latter image is just $G(\mathfrak o/\mathfrak{a})$. (This includes showing
 that the surjectivity of $G(\mathfrak o_\P)\to G^{(\P)}(\mathfrak o/\P)$ implies the surjectivity of  $G(\mathfrak o_\P)\to G^{(\P)}(\mathfrak o/\P^e)$, which
is explained in Lemma 5 on page 393 of \cite{LS} in the case $\mathfrak o_\P=\mathbb Z_p$.)

\subsection{Simplicity}
If $G$ is a Chevalley group and $G\hookrightarrow GL_n$ is the standard representation given by the Chevalley basis (see \cite{Carter2}, Chapter 4), then $G^{(\P)}(\mathfrak o/\P)$ is almost always simple modulo its centre and perfect \cite{Carter}.  The only exceptions are the groups in the following list:
$$ A_1(\mathbb F_2), A_1(\mathbb F_3), ^2\!A_2(\mathbb F_2), B_2(\mathbb F_2), ^2\!B_2(\mathbb F_2), G_2(\mathbb F_2), ^2\!F_4(\mathbb F_2),
^2\!G_2(\mathbb F_3).$$

\subsection{Some good groups}
All of our conditions are satisfied if
$G$ is a Chevalley group of rank at least $2$, $G\hookrightarrow GL_n$ is the standard representation, and if $G$ is  $B_2$ or $G_2$ then $S$ contains all prime divisors of $2$ with residue degree $1$.
One can extract from the references quoted above that there are also many other groups with the properties we require.

\section{Special linear and symplectic groups} Let $k$ be a number field.
In this section we restrict to special linear and symplectic groups over rings $\mathfrak o_S$ of $S$-integers.  For $\mathfrak o_S=\mathbb Z$, of course, we have already established optimal bounds for actions on homology spheres and acyclic homology manifolds.

\subsection{Linear actions.}     Let $k$ and $\mathfrak o_S$ be as above. Theorem~\ref{Smith} tells us that
whenever it acts on $\S^{n-2}$ or $\mathbb R^{n-1}$,  the subgroup ${\rm{SL}}(n,\mathbb Z)\leq {\rm{SL}}(n,\mathfrak o_S)$ must act trivially.  Therefore the kernel of the action contains elements which are not central, and are not central mod $\mathfrak p$ for any prime $\mathfrak p$.  Thus  the entire action is trivial by Proposition~\ref{Normal}.  The same observation shows that ${\rm{Sp}}(2n,\mathfrak o_S)$ cannot act on $\mathbb R^{2n-1}$ or $\S^{2n-2}$.

If $k$ is totally real, then ${\rm{SL}}(n,\mathfrak o_S)$ and ${\rm{Sp}}(2n,\mathfrak o_S)$ embed in ${\rm{SL}}(n, \mathbb R)$ and ${\rm{SL}}(2n, \mathbb R)$, respectively, so the above
observations provide a complete answer to Question \ref{q:smallest} in these cases.
But if  $k$ is not totally real,   then we have only that
${\rm{SL}}(n,\mathfrak o_S)$ acts on $\mathbb C^n$ (real dimension $2n$) and $\S^{2n-1}$; and similarly,
 ${\rm{Sp}}(2n,\mathfrak o_S)$ acts non-trivially on $\mathbb R^{4n}$  and $\S^{4n-1}$. Thus, in this case, we have not determined the least dimension of a
 contractible manifold or sphere on which these groups act non-trivially.

\subsection{Roots of unity.}  Suppose that $\mathfrak o$ contains a  $p$-th root of unity  $\eta$ for some odd prime $p$.  The subgroup
$D$ of ${\rm{SL}}(n,\mathfrak o_S)$ generated by diagonal matrices with powers of $\eta$ on the diagonal is a copy of $(\mathbb Z/p)^{n-1}$.   By Theorem~\ref{Smith},
whenever ${\rm{SL}}(n,\mathfrak o_S)$ acts on $\mathbb R^{2n-3}$ or $\S^{2n-4}$, some element of $D$ acts trivially.
$D$ intersects the centre of  ${\rm{SL}}(n,\mathfrak o_S)$ if and only if $n$ is divisible by $p$.  Moreover,
if $p$ does not divide $n$, then no element of $D$ become central modulo any prime $\mathfrak p$.

If $p$ does not divide $n$,  we conclude that any action of ${\rm{SL}}(n,\mathfrak o_S)$ on a Euclidean space of dimension less than $2n-2$ is trivial, i.e. ${\rm{SL}}(n,\mathfrak o_S)$ can't act on $\mathbb R^{2n-3}$.    Since there are actions on $\mathbb R^{2n}$ and $\S^{2n-1}$, this leaves only two dimensions in which there may or may not be actions.

If $p$ divides $n$,  then ${\rm{SL}}(n,\mathfrak o_p)$ contains a non-central copy of $(\mathbb Z/p)^{n-2}$, so ${\rm{SL}}(n,\mathfrak o_p)$ can't act on $\mathbb R^{2n-5}$.

\section{Questions}

The Smith theory techniques underlying this note rely on the presence of  $p$-torsion.  To determine the scope  of their
applicability, one needs to address the following natural question:

\begin{question}\label{q:smallest}  Given an affine group scheme $\mathbb G$, a ring of $S$-integers $\mathfrak o_S$, and a rational prime $p$,
 what is the rank of the largest elementary abelian $p$-group in $\mathbb G(\mathfrak o_S)$?
\end{question}

There are obvious bounds coming from
the observation that the generators of any such elementary abelian $p$-subgroup are all simultaneously diagonalizable over $\mathbb C$ (as can be seen easily from considering their Jordan forms), so the rank is at most the dimension   of a maximal torus in $G(\mathbb C)$.  For example, in ${\rm{SL}}_n$, it is at most $n-1$.

Smith theory techniques tell us  nothing in the absence of torsion, and little is known about how torsion-free lattices in semi-simple Lie groups
of higher rank might act on contractible manifolds or spheres. Such problems fit naturally into the context of
the Zimmer Programme \cite{zimmProg}, but the techniques developed until now are not well-adapted to actions by homeomorphisms.

\bibliographystyle{siam}


\section{Introduction}
Cosmic Microwave Background (CMB) probes have recently broadened our knowledge of primordial acoustic oscillations to small angular scales,
extending previous measurements of temperature power spectrum (Wilkinson Microwave Anisotropy Probe 7, \cite{Komatsu2011}) up to $l \sim 3000$
with first compelling evidence of CMB lensing from the South Pole Telescope (SPT, \cite{k11}) and Atacama Cosmology Telescope (ACT, \cite{act_2011}).
The impact of small-scale CMB measurements and gravitational lensing on cosmology
is relevant \cite{lewis_challinor_2006} and can be used to constrain cosmological parameters and to address one of the major issues of present cosmology,
 that is to say the nature of dark energy  \citep{verde_spergel_2002, giovi_etal_2003, 2004PhRvD..70b3515A, acquaviva_baccigalupi_2006, hu_etal_2006, Sherwin:2011gv}.
The simplest framework for dark energy models considers dark energy as a
cosmological constant $\Lambda$, contributing to about $74\%$ of the total
energy density in the universe and providing late time cosmic acceleration,
while a Cold Dark Matter represents about $21\%$ ($\Lambda$CDM model). Though theoretically in good agreement with present
observations, a cosmological constant is somewhat unpleaseantly
affected by the coincidence and fine-tuning problems which seem unavoidable in
such a framework. Many alternative models have been proposed, though it is fair to say
that so far no one  completely avoids these problems. Some
encouraging arguments have been put forward in the framework of dynamical dark
energy models, where a scalar field (quintessence or cosmon) rolls down a
 suitable potential \cite{wetterich_1988, ratra_peebles_1988} possibly
interacting with dark matter \cite{amendola_2000, pettorino_baccigalupi_2008}
or gravity \cite{Matarrese:2004xa, Perrotta:2002sw} and therefore modifying the growth of structure. Usually, one of the features of
such dynamical dark energy models is to have a non-negligible amount of dark
energy at early times. The amount of early dark energy (\emph{early} referring
to the time of decoupling) influences CMB peaks in various ways and can be
strongly constrained when including small scale measurements, as shown for instance in Refs. \cite{calabrese_etal_2011,
reichardt_etal_2011}.

In this paper, we consider the case of coupled dark energy models, in which dark
matter particles feel an interaction, additional to gravity, mediated by the
dark energy scalar field.
Such an interaction introduces effectively a coupling between the evolution of
the dark energy scalar field and dark matter particles. In this sense, this
class of models is both an example in which a non-negligible amount of early
dark energy is present as well as a typical scenario of modified gravity
theories. When seen in the Jordan frame, a coupling between matter and dark
energy can be reformulated in terms of scalar-tensor theories (or $f(R)$  models).
This is exactly true when the contribution of baryons is neglected. In the
Einstein frame, it is common use to neglect a coupling to baryon within coupled
dark energy models, and consider only dark energy - dark matter interactions.
 Alternatively, in the Jordan frame, scalar-tensor theories ($f(R)$ models) require
some sort of screening mechanism (like chameleon \cite{khoury_etal_2004,
Hui:2009kc, 2012arXiv1204.3906U, 2011arXiv1102.5278D} or symmetrons
\cite{2011PhRvD..84j3521H})  that protects the dark energy scalar field and its
mass within high density regions, so that local solar system constraints are
satisfied.

The strength of the coupling affects CMB  in several ways, changing the
amplitude, the position of the peaks as well as contributing to the Late
Integrated Sachs-Wolfe (ISW) effect (manifest at large length scales) and to
gravitational lensing (appearing at small length scales in the temperature
spectrum). Moreover, the coupling is degenerate with the amount of cold dark
matter $\Omega_{c}$, the spectral index $n$, the Hubble parameter $H(z)$ (see
\cite{amendola_etal_2012} for a review).
After recalling the effects of the coupling on CMB, we use a Monte Carlo
analysis to constrain the coupling combining WMAP and SPT  real data.
Furthermore we extend our analysis to forecasting the constraints that Planck
data are expected to put on the coupling parameter, combined with mock SPT
data.

This paper is organized as follows. In Section II we recall the main features of
coupled dark energy (CDE) cosmologies. In section III we recall effects of the
coupling on the CMB spectrum and describe the methods used, both with regard to
the implementation of the numerical code and the data used for this paper. In
Section IV we derive
the constraints from existing data for several different runs, including effects
of the effective relativistic degree of freedom $N_{eff}$, CMB-lensing, curvature and
massive neutrinos. Here we also forecast the  constraining capability in
presence of the forthcoming Planck data, joined with SPT mock data. Finally, in Section V we derive our
conclusions.

 \section{Coupled Dark Energy} \label{cde}

 We consider the case in which an interaction is present between dark energy and
dark matter, as illustrated in
 \cite{Kodama:1985bj, amendola_2000, Amendola:2003wa, amendola_etal_2003,
amendola_quercellini_2003, pettorino_baccigalupi_2008}. Such cosmologies have to
be seen within the framework of modified gravity, since effectively the
gravitational interaction perceived by dark matter particles is modified with respect
to standard General Relativity.
 Coupled dark energy cosmologies considered here are described by the
lagrangian:
 \begin{equation}}\newcommand{\ee}{\end{equation}}\newcommand{\bea}{\begin{eqnarray}}\newcommand{\eea}{\end{eqnarray}}\newcommand{\beann}{\begin{eqnarray*}}\newcommand{\eeann}{\end{eqnarray*}}\newcommand{\nn}{\nonumber}\newcommand{\drv}{{\partial_r \delta v}}\newcommand{\dtv}{{\partial_t \delta v}}\newcommand{\at}[2]{\left. #2 \right|_{#1}}\newcommand{\divvec}[1]{\mathbf{\nabla\hspace{-.3em}\cdot\hspace{-.2em} #1}}\newcommand{\scalprod}[2]{\mathbf{#1 \hspace{-.2em}\cdot\hspace{-.2em} #2}}\newcommand{\exvalue}[1]{\langle #1 \rangle \label{L_phi} {\cal L} =
 -\frac{1}{2}\partial^\mu \phi \partial_\mu \phi - U(\phi) -
 m(\phi)\bar{\psi}\psi + {\cal L}_{\rm kin}[\psi] \,, \ee in which the mass of
 matter fields $\psi$ is a function of the scalar field $\phi
 $. We consider the case in which the DE is only coupled
 to CDM (hereafter denoted with a subscript $c$, while the subscript $b$ will denote baryons). In this
case, the coupling is not affected by tests on the equivalence principle and
solar system constraints and can therefore be stronger than the one with
baryons.
 The choice
 $m(\phi)$ specifies the coupling and as a consequence the quantity
 $Q_{(\phi) \mu}$ via the expression: \begin{equation}}\newcommand{\ee}{\end{equation}}\newcommand{\bea}{\begin{eqnarray}}\newcommand{\eea}{\end{eqnarray}}\newcommand{\beann}{\begin{eqnarray*}}\newcommand{\eeann}{\end{eqnarray*}}\newcommand{\nn}{\nonumber}\newcommand{\drv}{{\partial_r \delta v}}\newcommand{\dtv}{{\partial_t \delta v}}\newcommand{\at}[2]{\left. #2 \right|_{#1}}\newcommand{\divvec}[1]{\mathbf{\nabla\hspace{-.3em}\cdot\hspace{-.2em} #1}}\newcommand{\scalprod}[2]{\mathbf{#1 \hspace{-.2em}\cdot\hspace{-.2em} #2}}\newcommand{\exvalue}[1]{\langle #1 \rangle Q_{(\phi) \mu} = \frac{\partial
   \ln{m(\phi)}}{\partial \phi} \rho_c \, \partial_\mu \phi . \ee

$Q_{(\phi) \mu} $   acts as a source term in the Bianchi identities:
     \begin{align}
T_{\nu:\mu}^{\mu} = Q_{(\phi) \mu}
\label{tensor_conserv_alpha}
\end{align}

If no other species is involved in the coupling then $Q_{(c) \mu} = - Q_{(\phi)
\mu}$.
 Various choices of couplings have been investigated in literature, including
constant $\beta $
\cite{Wetterich1994,Amendola1999,Mangano_etal_2003,Amendola2004,Koivisto:2005nr,
Guo:2007zk,Pettorino:2008ez,quercellini_etal_2008,Quartin:2008px,
Valiviita:2009nu} and varying couplings \cite{ Baldi:2010vv}.
 For a constant coupling, typical values of $\beta$ presently allowed by
observations (within current CMB
 data) are within the range $0\le \beta < 0.06$ (at 95\% CL for a constant
coupling and an exponential potential) \cite{amendola_etal_2003,
Amendola:2003wa, amendola_quercellini_2003, Bean:2008ac}, or possibly more
\cite{lavacca2009, Kristiansen:2009yx} if neutrinos are taken into account or
for more realistic time-dependent choices of the coupling.
 Analysis of the
models and constraint on these couplings have been obtained in several other
ways, including spherical collapse (\cite{Wintergerst:2010ui, Mainini:2006zj}
and references therein), higher-order expansions with the time renormalizazion group \cite{Saracco:2009df},
$N$-body simulations \cite{Baldi:2010td,Baldi_etal_2010,Baldi:2010vv} and
effects on supernovae, CMB and cross-correlation of CMB and LSS
\cite{amendola_etal_2003,amendola_quercellini_2003,Bean:2008ac,lavacca2009,
Kristiansen:2009yx,Mainini:2010ng,DeBernardis:2011iw,Xia:2009zzb,martinelli}
together with Fisher matrix forecasts analysis combining power spectrum and Baryonic Acoustic Oscillations
measurements as expected by the Euclid satellite \footnote{http://www.euclid-ec.org/ and http://sci.esa.int/euclid} and CMB as expected from Planck
\cite{amendola_etal_2012}.

 The zero-component of equation (\ref{tensor_conserv_alpha}) gives the
 conservation equations for the energy densities of each species:
 \bea
 \label{rho_conserv_eq_phi} \rho_{\phi}' &=& -3 {\cal{H}} \rho_{\phi} (1 + w_\phi) - Q_{(\phi)0} \,\,\,\, , \\
 \label{rho_conserv_eq_c} \rho_{c}' &=& -3 {\cal{H}} \rho_{c} + Q_{(\phi)0} \pp
 \nonumber \eea Here we have treated each component as a fluid with
 ${T^\nu}_{(\alpha)\mu} = (\rho_\alpha + p_\alpha) u_\mu u^\nu + p_\alpha
 \delta^\nu_\mu$, where $u_\mu = (-a, 0, 0, 0)$ is the fluid 4-velocity and
 $w_\alpha \equiv p_\alpha/\rho_\alpha$ is the equation of state. Primes denote derivative with respect to conformal time $\tau$.  The class of models considered here corresponds to the choice: \begin{equation}}\newcommand{\ee}{\end{equation}}\newcommand{\bea}{\begin{eqnarray}}\newcommand{\eea}{\end{eqnarray}}\newcommand{\beann}{\begin{eqnarray*}}\newcommand{\eeann}{\end{eqnarray*}}\newcommand{\nn}{\nonumber}\newcommand{\drv}{{\partial_r \delta v}}\newcommand{\dtv}{{\partial_t \delta v}}\newcommand{\at}[2]{\left. #2 \right|_{#1}}\newcommand{\divvec}[1]{\mathbf{\nabla\hspace{-.3em}\cdot\hspace{-.2em} #1}}\newcommand{\scalprod}[2]{\mathbf{#1 \hspace{-.2em}\cdot\hspace{-.2em} #2}}\newcommand{\exvalue}[1]{\langle #1 \rangle \label{coupling_const}
 m(\phi) = m_0 e^{-\beta \frac{\phi}{M}} \,, \ee with the coupling term
 equal to \begin{equation}}\newcommand{\ee}{\end{equation}}\newcommand{\bea}{\begin{eqnarray}}\newcommand{\eea}{\end{eqnarray}}\newcommand{\beann}{\begin{eqnarray*}}\newcommand{\eeann}{\end{eqnarray*}}\newcommand{\nn}{\nonumber}\newcommand{\drv}{{\partial_r \delta v}}\newcommand{\dtv}{{\partial_t \delta v}}\newcommand{\at}[2]{\left. #2 \right|_{#1}}\newcommand{\divvec}[1]{\mathbf{\nabla\hspace{-.3em}\cdot\hspace{-.2em} #1}}\newcommand{\scalprod}[2]{\mathbf{#1 \hspace{-.2em}\cdot\hspace{-.2em} #2}}\newcommand{\exvalue}[1]{\langle #1 \rangle Q_{(\phi)0} = - \frac{\beta}{M} \rho_c \phi' \,. \ee
 Equivalently, the scalar field evolves according to the Klein-Gordon equation, which now includes an extra term that depends on CDM energy density:
 \begin{equation}}\newcommand{\ee}{\end{equation}}\newcommand{\bea}{\begin{eqnarray}}\newcommand{\eea}{\end{eqnarray}}\newcommand{\beann}{\begin{eqnarray*}}\newcommand{\eeann}{\end{eqnarray*}}\newcommand{\nn}{\nonumber}\newcommand{\drv}{{\partial_r \delta v}}\newcommand{\dtv}{{\partial_t \delta v}}\newcommand{\at}[2]{\left. #2 \right|_{#1}}\newcommand{\divvec}[1]{\mathbf{\nabla\hspace{-.3em}\cdot\hspace{-.2em} #1}}\newcommand{\scalprod}[2]{\mathbf{#1 \hspace{-.2em}\cdot\hspace{-.2em} #2}}\newcommand{\exvalue}[1]{\langle #1 \rangle \label{kg} \phi'' + 2{\cal H} \phi' + a^2 \frac{dV}{d \phi} = a^2 \beta
 \rho_{c} \,\, . \ee
 Throughout this paper we choose an inverse power law potential defined as:
 \begin{equation}}\newcommand{\ee}{\end{equation}}\newcommand{\bea}{\begin{eqnarray}}\newcommand{\eea}{\end{eqnarray}}\newcommand{\beann}{\begin{eqnarray*}}\newcommand{\eeann}{\end{eqnarray*}}\newcommand{\nn}{\nonumber}\newcommand{\drv}{{\partial_r \delta v}}\newcommand{\dtv}{{\partial_t \delta v}}\newcommand{\at}[2]{\left. #2 \right|_{#1}}\newcommand{\divvec}[1]{\mathbf{\nabla\hspace{-.3em}\cdot\hspace{-.2em} #1}}\newcommand{\scalprod}[2]{\mathbf{#1 \hspace{-.2em}\cdot\hspace{-.2em} #2}}\newcommand{\exvalue}[1]{\langle #1 \rangle \label{potential}
  V = V_0 \phi^{-\sigma}
 \ee
 with $\sigma$ and $V_0$ constants.

\section{The coupling effect on the CMB power spectrum}
As discussed in \cite{amendola_etal_2012}, the coupling has two main effects on the CMB: 1) it shifts the position
of the acoustic peaks to larger $\ell$'s due to the increase in the
last scattering surface distance (sometimes called projection effect,
\citep{Pettorino:2008ez} and references therein); 2) it reduces the
ratio of baryons to dark matter at decoupling with respect to its
present value, since coupled dark matter dilute faster than in an
uncoupled model. Both effects are clearly visible in Fig. (\ref{fig:CMB-for-three})
for various values of $\beta$.

\begin{figure}
\includegraphics[scale=1]{plots_cq/lensedcl_wmap7.eps}
\caption{\label{fig:CMB-for-three}CMB TT temperature spectra for
three values of $\beta$. Data are taken from WMAP7 \cite{Larson:2010gs}.}
\end{figure}

In Fig.\ref{fig:CMB-for-three_timeslsquared_lensing} the effect of the coupling on the CMB is more evident. The figure shows the quantity $C_l(l+1)l^3/(2\pi)$
(i.e. the usual TT spectrum plotted in Fig.(\ref{fig:CMB-for-three}) multiplied by $l^2$), as suggested for example in \cite{Hlozek:2011pc}.
 In Fig.\ref{fig:CMB-for-three_timeslsquared_lensing}  we also show the effect
of including CMB-lensing on the TT spectrum for the same value of the coupling
constant: the unlensed TT spectrum (blue dot dashed) clearly differs from the
corresponding lensed one (solid blue), an effect which is larger at small scales
(large l). The position of the peaks remains invariant but the amplitude is
larger and, for $\ell\lsim 2000$,  the throats appear more pronounced than the peaks.

\begin{figure}
\includegraphics[scale=1]{plots_cq/cq_cl}
   \includegraphics[width=13.cm]{plots_cq/cq_lensingpotential.eps}%
\caption{\label{fig:CMB-for-three_timeslsquared_lensing}CMB TT temperature spectra (top panel) and dimensionless lensing potential $l(l+1)C^{\phi \phi}/(2\pi)$  versus multipole $l$ for three values of the coupling $\beta$ .}
\end{figure}

The purpose of our analysis is to use recent CMB observations reaching multipoles up to $l \sim 3000$ to
constrain dynamical dark energy models in which a coupling is present between
the quintessence scalar field (or cosmon, seen as the mediator of a fifth force
between dark matter particles) and dark matter.
In order to do so, we proceed as follows.

\subsection{Theoretical spectra}
Theoretical CMB and lensing power spectra have been produced using the code IDEA
(Interacting Dark Energy Anisotropies) based on CAMB \cite{Lewis:1999bs} and
able to include dynamical dark energy, Early Dark Energy parameterizations (not
included in this analysis) as well as  interacting dark energy models.
In order to include the coupling, both background and linear perturbations have
been modified following Refs. \cite{Amendola1999, Pettorino:2008ez}.
The output has been compared to an independent
code \cite{amendola_quercellini_2003} that is built on CMBFAST and
the agreement was better than 1$\%$. The difficulty in the implementation relies
on the fact that the initial conditions cannot be obtained analytically as in
simple dark energy parameterizations (early dark energy or $(w_0, w_a)$):
instead, they must be found by trial and error, through an iterative routine
that finds the initial conditions required to get the desired present values of
the cosmological parameters.

We have then performed a Monte Carlo analysis integrating IDEA within COSMOMC
\cite{cosmomc_lewis_bridle_2002} comparing our theoretical predictions with the
data presented in the next subsection.

We recall that the CMB coming from the last scattering surface (LSS) is bent by gravitational
structures on the path towards us; this effect is called CMB-lensing
\citep{Bartelmann:1999yn, lewis_challinor_2006}.
The standard deviation of the deflection angle is of the order of 2 arcminutes, which would correspond to
small scales and $l > 3000$ multipoles, where CMB peaks are already damped by
photon diffusion. However, deflection angles are correlated with each other over
degree scales, so that lensing can have an important effect on the scales of
the primary acoustic peaks, mainly smoothing them and transferring power to
larger multipoles. Recently, CMB lensing detection has been claimed by several groups analyzing
total intensity CMB anisotropies \citep{act_2011, k11}.

CMB lensing is a probe of the expansion rate and naturally depends on the growth of perturbation and on
the gravitational potentials; since dark energy affects both aspects, CMB
lensing represents a way to discriminate among dynamical Dark Energy models and
$\Lambda$CDM, with promising results \citep{verde_spergel_2002, giovi_etal_2003, 2004PhRvD..70b3515A, acquaviva_baccigalupi_2006, hu_etal_2006, Sherwin:2011gv}.
In particular, the CMB weak lensing theory in generalized cosmologies has been outlined in \cite{2004PhRvD..70b3515A}.
The difference between the lensed and unlensed curves in Fig.\ref{fig:CMB-for-three} shows the typical effects from CMB lensing. The
acoustic peaks are smeared because of the correlation between different scales induced by lensing, and for
the same reason a fraction of power is transferred to the angular domain corresponding to the damping tail, therefore
dominating that part of the spectrum. Earlier works \cite{acquaviva_baccigalupi_2006} have pointed out how the
lensing is most relevant in particular in early dark energy models, as it injects power at the onset
of cosmic acceleration, $z\simeq 1\pm 0.5$ constraining the dark energy abundance in the corresponding epoch.
On these lines, the analysis in \citep{Sherwin:2011gv} has shown that the inclusion of lensing data
promote the CMB alone to be a probe of the existence of dark energy,
breaking geometrical degeneracies associated to the pure CMB anisotropies at
last scattering.

Also, the lensing depends on time, combining information from decoupling (from
the last scattering surface of the CMB) and $z < 5$ (when large scale structures
formed); recent studies
\cite{carbone_etal_2009, Teyssier_etal_2009} focus on implementing and investigating simulations of
CMB lensing through cosmological structures in N-body simulations.

During matter dominated era (MDE), the potentials encountered along the way are constant in the linear regime and the gradient of the potential causes a total deflection angle given by:
\begin{equation}
\alpha = -2 \int_0^{\chi_{*}} d \chi \frac{f_K(\chi_* - \chi)}{f_K(\chi_*)} \nabla_{\perp} \Psi (\chi \hat{n}; \tau_0 - \chi)\ ,
\label{alpha_lensing}
\end{equation}
where $\chi_*$ is the conformal distance of the source acting as a lens, $\Psi$ is its gravitational potential, $\eta_0 - \chi$ is the conformal time at which the CMB photon was at position $\chi \hat{n}$.

One can also define the gravitational lensing potential
\begin{equation}
\psi(\hat{n}) \equiv -2 \int_0^{\chi_*} d \chi \frac{f_K(\chi_*-\chi)}{f_K(\chi_*)f_K(\chi)} \Psi(\chi \hat{n}; \eta_0 - \chi)\ .
\end{equation}
The lensed CMB temperature $\tilde{T}_{\hat{n}}$ in a direction $\hat{n}$ is
given by the unlensed temperature in the deflected direction $\tilde{T}(\hat{n})
= T(\hat{n}') = T(\hat{n} + \alpha)$ where at lowest order the deflection angle
$\alpha = \nabla \psi$ is just the gradient of the lensing potential.
Expanding the lensing potential into spherical harmonics, one can define also
the angular power spectrum $C_l^\psi$ corresponding to the lensing potential,
defined as $<\psi_{lm} {\psi^*}_{l'm'} > = \delta_{ll'}
\delta_{mm'} C_l^{\psi}$; the latter (multiplied by $[l(l+1)]^2$) is shown in
Fig.\ref{fig:CMB-for-three_timeslsquared_lensing} (lower panel) for different values of the coupling. As we can see
from the plot, the CMB lensing potential mainly gives contribution to large
scales up to $l \sim 1000$ or less. However, the lensed CMB temperature power
spectrum  depends on the convolution between the lensing potential and the
unlensed temperature spectrum (see \cite{lewis_challinor_2006} for more details)
whose effect is of several percent at $l > 1000$, thus being important when
estimating the spectrum up to small scales of $l \sim 3000$, as for the data we consider in the following.

\subsection{Observed data}
We have compared theoretical predictions of CMB lensed spectra with two
datasets. The first includes WMAP7 temperature spectra \cite{Komatsu2011}. The
second includes the recently released power spectrum data from SPT \cite{k11}.
Together with the ATC \cite{act_2011}, SPT
\cite{k11} has recently shown evidence of CMB-lensing with an enhancement in the
CMB temperature power spectrum up to $l \sim 3000$. We do not use ACT data in
this analysis.

In order to use SPT data, we have installed the likelihood provided by the SPT
team \cite{k11} on their SPT website
\footnote{http://pole.uchicago.edu/public/data/keisler11/index.html} and
integrated it with the recommended version of COSMOMC
\footnote{http://cosmologist.info/cosmomc} (August 2011). We have then
implemented IDEA on this version.
Care has to be used whenever small multipoles in the range $2000 - 3000$ are
used, due to several sources of foregrounds active on those scales.
In particular, whenever SPT data are included in the analysis we also
marginalize over the three nuisance parameters described in \cite{k11}: two of
them refer to Poisson point sources and clustered point sources; the third one
adds power from the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects, which
are also an example of secondary CMB anisotropies.
These effects are relevant when small scales ($l \gsim 2000$) are included and
have to be taken into account whenever SPT data are used.
In Fig.(\ref{fig:cq_data}) and Fig.(\ref{fig:cq_data_l2}) we show the superposition of WMAP7 and SPT data, together with the best fit of the combined set.

\begin{figure*}
\centering
   \includegraphics[width=15.cm]{plots_cq/cq_data.eps}%
     \caption{\small WMAP7 and SPT data used for this work. SZ effects are not included in the best fit line (solid black) (see fig.5 of \cite{k11} for the best fit including SZ). The best fit from WMAP+SPT, for a $\Lambda$CDM without foregrounds is also shown (dash-dotted black line). This plot is similar to fig.5 of \cite{k11} (without SZ in the best fit); we reproduced it here for convenience and because it will allow us to neglect SZ when combining Planck with mock SPT data, in section V. All foregrounds are instead included when SPT is combined with WMAP7.}
\label{fig:cq_data}
\end{figure*}
\normalsize

\begin{figure*}
\centering
   \includegraphics[width=15.cm]{plots_cq/cq_data_l2.eps}%
     \caption{\small WMAP7 and SPT data used for this work, with an extra factor of $\ell^2$. SZ effects are not included in the best fit line (solid black) (see fig.5 of \cite{k11} for the best fit including SZ). The best fit from WMAP+SPT, for a $\Lambda$CDM without foregrounds is also shown (dash-dotted black line).}
\label{fig:cq_data_l2}
\end{figure*}
\normalsize

The baseline set of parameters includes $\Theta = {\Omega_b h^2, \Omega_{c} h^2, \Omega_{de},
\theta_s, log{\cal{A}}, n_s, \tau}$. As in other papers, these parameters represent the fractional abundances of the
various species, as well as the amplitude and shape of the primordial power spectrum, and the reionization optical depth;
as stressed already, we also consider the three nuisance parameters when SPT is included; when we  impose spatial flatness,
the present dark energy density  $\Omega_{de}$ becomes a derived parameter;
in addition, when CQ is included two more parameters are added: $\beta$ and
$\sigma$. Again, $\beta$ represents the coupling between dark matter particles in eq.\ref{coupling_const},
$\sigma$ is the parameter in the scalar field potential  (\ref{potential}) that drives the long
range interaction. Different runs are illustrated in Tab.\ref{tab:runs}. The Helium abundance $Y_{He}$ is derived following BBN consistency (see \cite{k11} for details).

 \begin{center}
\begin{table}
\begin{tabular}{cccccccc}
Run & \begin{minipage}{100pt}
\center
CMB-Lensing \\
\end{minipage} &
\begin{minipage}{50pt}
\center
WMAP7
\end{minipage} &
\begin{minipage}{50pt}
\center
SPT
\end{minipage} &
\begin{minipage}{50pt}
\center
Planck
\end{minipage} &
\begin{minipage}{50pt}
\center
Parameters
\end{minipage} &
\\
\hline
$ \emph{cq1}  $ & \checkmark & \checkmark & \checkmark & X & baseline + $\sigma$ + $\beta$ \\
$ \emph{cq2}  $ & \checkmark & \checkmark & X & X &baseline + $\sigma$ + $\beta$\\
$ \emph{cq4}  $ & \checkmark & \checkmark & \checkmark &X & baseline + $\sigma$ + $\beta$ + $N_{\nu}$\\
$ \emph{cq1NL}  $ & X & \checkmark & \checkmark &X & baseline + $\sigma$ + $\beta$ \\
$ \emph{cq1K}  $ &  \checkmark & \checkmark & \checkmark &X & baseline + $\sigma$ + $\beta$ + curvature \\
$ \emph{cq1AL}  $ &  \checkmark & \checkmark & \checkmark &X & baseline + $\sigma$ + $\beta$ + $A_l$ \\
$ \emph{cq1}\nu $ &  \checkmark & \checkmark & \checkmark &X & baseline + $\sigma$ + $\beta$ + $f_\nu$ \\
$ \emph{cq1hst} $ &  \checkmark & \checkmark & \checkmark &X & baseline + $\sigma$ + $\beta$ + HST + BAO + SNae\\
$ \emph{cq1Pl} $ &  \checkmark & X & \checkmark (mock) & \checkmark (mock) & baseline + $\sigma$ + $\beta$ + HST + BAO + SNae\\
\hline \end{tabular}
\caption{COSMOMC Monte Carlo simulation runs described in this paper; they all refer to CQ models.} \label{tab:runs}
\end{table}
\par\end{center}

\normalsize

\section{Results} \label{sec:results}
As a first step, we have  performed a run using WMAP7 and SPT with a
$\Lambda$CDM model, and we find results compatible with \cite{k11}.
Note for the following, that in \cite{k11} the authors report the mean values of
each parameter, together with its standard deviation. We instead report the best
fit values and the marginalized errors at 68$\%$ and 95$\%$ confidence level. We
now describe results from the runs illustrated in Tab.\ref{tab:runs}.

\subsection{Baseline plus $\beta$ and $\sigma$}
The first two Monte Carlo runs we describe (\emph{cq1} and \emph{cq2}) use the baseline
set of parameters $\Theta \equiv \{\Omega_b h^2, \Omega_{c} h^2, \theta_s,
log{\cal{A}}, n_s, \tau \}$; in addition, two more parameters are added to account
for the coupling ($\beta$) and the dark energy scalar field potential
($\sigma$).
Run \emph{cq1} compares theoretical spectra with WMAP7+SPT data; run  \emph{cq2} includes
WMAP7 data only. Results are shown in Tab.\ref{tab:cq1cq2_bestfit}, where we report best fit values with 1-sigma (68\%) errors on various parameters. For $\beta$ we also
report upper 1 and 2 $\sigma$ limits.
When WMAP7 only is considered, the coupling is constrained to be $\beta <
0.078 (0.14)$ at 1 (2) $\sigma$, while including SPT data, and therefore small scales and large
multipoles, the coupling is constrained down to $< 0.063 (0.11)$.

In the past, \cite{amendola_etal_2003} found $\beta_A < 0.16 $ at 95\% c.l.
(note that our definition of $\beta$ is $\beta = \sqrt(2/3) \beta_A = 0.13$ at 2
$\sigma$) using
COBE, Boomerang, Maxima, DASI and fixing the optical depth $\tau$.
\cite{amendola_quercellini_2003} found $\beta < 0.061 (0.11)$ at 1 (2) $\sigma$
using WMAP1 data and $\beta < 0.11 (0.16)$ for pre-WMAP data (using a different
set than \cite{amendola_etal_2003}).
Ref. \cite{Bean:2008ac} found C < 0.034 (0.066) (in their notation $C \equiv
\beta_A \equiv \sqrt{2/3} \beta$ so that their result is equivalent to $\beta <
0.028 (0.054)$) using WMAP5, SNLS, HST, LRG, SDSS.
See also \cite{Xia:2009zzb} (with a different definition of the coupling) and
\cite{lavacca2009} (where massive neutrinos were included).
These constraints are not easily comparable to ours since in these papers different
priors have been used and/or some parameters have been kept fixed.
Overall, however, we find an agreement with the more recent  constraints to within a factor
of 50\% at most.

The 2D confidence contours are plotted in Fig.\ref{fig:like_cont_baseline_cq1_cq2}.
Here we show a selection of the most interesting likelihood contours vs the
coupling $\beta$.
In Fig.\ref{fig:like_cont_1D_cq1_cq2} we also show 1D likelihood contours,
comparing results from WMAP7+ SPT with WMAP7 only. Note that there is no
dependence of cosmological parameters from $\sigma$, as expected since $\sigma$ only affects late time cosmology; the
range in $\sigma$ was therefore safely chosen to be between 0.13 and 0.5, small enough to get
reasonable speed for the runs. The value of $w$ is arbitrary and approximately related to
$\sigma$ via the expression: $w = - 2 / (\sigma + 2)$; the interval chosen for
$\sigma$ is such that $w$ still assumes reasonable values, at least smaller than
-0.8. As it appears clearly from Fig. \ref{fig:like_cont_1D_cq1_cq2}, CMB is practically insensitive
to $\sigma$ or $w$ within the range we consider.

 \begin{center}
\begin{table}
\begin{tabular}{lllll}
\textbf{Best fit values} & \textbf{for coupled quintessence} & \\
\hline
\begin{minipage}{30pt}
\flushleft
Parameter
\\
\end{minipage} &
\begin{minipage}{150pt}
\flushleft
cq1 (WMAP7 + SPT)
\end{minipage} &
\begin{minipage}{150pt}
\flushleft
cq2 (WMAP7)
\end{minipage} &
\\
\\
\hline
$  \Omega_b h^2 $ & $0.022^{ +0.0007}_{-0.00013}$           & $0.023^{+0.00044}_{-0.00070}$		\\
$ \Omega_{c} h^2  $ &$0.11^{+0.0022}_{-0.011}$                                  & $0.11^{+0.0019}_{-0.016}$ \\
      $ \theta_s$ & $1.04^{+0.0024}_{-0.00072}$                    & $1.04^{+0.0027}_{0.0025} $ \\
   $ \tau$ & $0.091^{+0.0013}_{-0.012} $                             & $0.089^{0.0076}_{0.0072}$\\
   $ n_s$ & $0.96^{+0.019}_{-0.0056} $                                                     & $0.97^{+0.021}_{-0.012}$ \\
   $w$ &  $-0.88^{+0.080}_{-0.12}  $                                                                                & $-0.97^{+0.17}_{-0.03}$ \\
   $\beta$ &  $0.012^{+0.050}_{-0.012}$                          & $0.0066^{0.071}_{-0.0066}$\\
      $\beta$ &  $< 0.063 (0.11) $                          & $< 0.078 (0.14)$\\
    $ \sigma$ & $0.22^{+0.28}_{-0.090}$                                              & $0.13^{+0.37}_{-0.0048}$ \\

   $ \Omega_{de}$ & $0.72^{+0.076}_{-0.012}$                                  & $0.72^{+0.093}_{-0.016} $ \\
   $ Age/Gyr$ & $13.8^{+0.012}_{-0.39} $                                              & $13.8^{+0.07}_{-0.5}$ \\
   $ z_{re}$ & $10.9^{+0.63}_{-1.8} $                          & $10.6^{+1.2}_{-1.2}$ \\
   $ H_0$ & $68.7^{+8.4}_{-0.98}$                             & $ 69.3^{+10.7}_{-1.7} $  \\
            $ D_{3000}^{SZ}$ & $4.0^{+2.1}_{-4.0}$                    & - \\
                  $ D_{3000}^{PS}$ & $21.5^{+1.7}_{-3.7}$                    &- \\
                           $ D_{3000}^{CL}$ & $5.3^{+1.9}_{-2.4} $                    & - \\
   -$Log$(Like) &$ 3756$                      							& $3737 $  \\
\hline
\end{tabular}
\caption{Best fit values and 1-$\sigma$ errors comparing runs \emph{cq1} and \emph{cq2}. Both runs include coupling; cq1 uses WMAP7 + SPT while \emph{cq2} uses WMAP7 only. For $\beta$ we also write in brackets the value of the 2$\sigma$ marginalized error.}
 \label{tab:cq1cq2_bestfit}
\end{table}
\par\end{center}
\normalsize

 \begin{figure*}
 \centering
 \includegraphics[width=18.cm]{plots_cq/compare_cq1_v10_cq2_v10_showmore_andro.eps}
      \caption{\small Confidence contours for the cosmological parameters for coupled quintessence models. We compare runs \emph{cq1} (red) and \emph{cq2} (white). The light blue asterisks mark the best fit points for $cq1$.  1-sigma and 2-sigma contours are shown.}
 \label{fig:like_cont_baseline_cq1_cq2}
 \end{figure*}
 \normalsize

 \begin{figure*}
 \centering
    \includegraphics[width=17.cm]{plots_cq/compare1D_cq1_v10_cq2_v10_sptwmap_andro_more.eps}%
      \caption{\small 1D likelihoods for the cosmological parameters for coupled quintessence models. We compare runs  \emph{cq1} (solid blue) and  \emph{cq2} (dashed red). As expected, in coupled quintessence models, no dependence is seen on the value of $\sigma$ and $w$, which are arbitrary. We restrict the analysis to a reasonable range in $w$ and a range in $\sigma$ small enough to speed up the iterative routine that finds initial conditions.}
 \label{fig:like_cont_1D_cq1_cq2}
 \end{figure*}
 \normalsize

\subsection{Effective number of relativistic species $N_{eff}$}
The effective number of relativistic species before recombination, usually
denoted by $N_{eff}$
is higher than the number of relativistic neutrino species (3.046) due to
photons produced in electron-positron annihilation at the end of neutrino
freeze-out \cite{dicus_etal_1982,lopez_etal_1999, mangano_etal_2005,
2011JCAP...03..035M}.
Using WMAP7 + SPT data for a $\Lambda$CDM, $N_{eff} = 3.046$ was found to be
preferred over zero relativistic species ($N_{eff} = 0$) \cite{k11}.
If $N_{eff}$ is left free to vary and marginalized over, its best fit value can
be even larger: \cite{Komatsu2011} found $N_{eff} > 2.7$ at $95\%$ CL using
WMAP7 alone; \cite{dunkley_etal_2010}  found $N_{eff} = 5.3 \pm 1.3$ using WMAP7
+ ACT; \cite{k11} found $N_{eff} = 3.85 \pm 0.62$ using WMAP7 + SPT.
If relativistic species are present, the expansion rate
during radiation dominated era increases \cite{hu_white_1996,
Bashinsky_Seljak_2004, 2011arXiv1104.2333H,  k11}.
We have redone the analysis in the case of coupled quintessence, to check whether the number of relativistic degrees of freedom is degenerate with the coupling.
The effect of marginalizing over $N_{eff}$ on the coupling is shown in
Tab.\ref{tab:cq1cq4_bestfit}. In Fig.\ref{fig:like_cont_neff} we plot the
likelihood contours for a selection of parameters, vs the coupling $\beta$,
comparing different runs.

Contours are larger when we allow $N_{eff}$ to vary, but the range in $\beta$ is
not affected considerably ($\beta < 0.074 (0.12)$ instead of $\beta < 0.063 (0.11)$). When a
coupling is present, we find that the best fit for the number of relativistic
species is given by $N_{eff} = 3.84^{+0.74}_{-0.49}$ when using WMAP and SPT, similar
to the value mentioned before and evaluated in absence of a coupling.  The
allowed range for $Y_{He}$, obtained through BBN consistency, increases a lot when $N_{eff}$ is free to vary
\cite{Bashinsky_Seljak_2004}. $N_{eff}$ is degenerate with dark matter
and the spectral index, which in turn are weakly degenerate with $\beta$, though no
direct degeneracy appears between $\beta$ and $N_{eff}$, as shown in Fig.\ref{fig:like_cont_neff_yhe}.

It is interesting to see (Fig.\ref{fig:like_cont_neff}) that when we allow for an effective number of relativistic degrees of freedom, marginalizing over $N_{eff}$, the coupling from WMAP7+SPT data increases to a best fit value of $\beta \sim 0.03$. A larger value of $N_{eff}$ (best fit $\sim 3.8)$ favors larger couplings between dark matter and dark energy as well as values of the spectral index closer to -1 ($n_s \sim 0.99$).

\begin{center}
\begin{table}
\begin{tabular}{lllll}
\textbf{Best fit values} & \textbf{for coupled quintessence} & WMAP7 + SPT\\
\hline
\begin{minipage}{30pt}
\flushleft
Parameter
\\
\end{minipage} &
\begin{minipage}{150pt}
\flushleft
 \emph{cq1} (baseline + $\beta$+$\sigma$)
\end{minipage} &
\begin{minipage}{150pt}
\flushleft
 \emph{cq4} (baseline + $\beta$+$\sigma$+$N_{eff}$)
\end{minipage} &
\\
\\
\hline
$  \Omega_b h^2 $ & $0.022^{ +0.0007}_{-0.00013}$           & $0.023^{+0.00054}_{-0.00054}$		\\
$ \Omega_{c} h^2  $ &$0.11^{+0.0022}_{-0.011}$                                  & $0.12^{+0.0076}_{-0.016}$ \\
      $ \theta_s$ & $1.04^{+0.0024}_{-0.00072}$                    & $1.04^{+0.0012}_{-0.0025} $ \\
   $ \tau$ & $0.091^{+0.0013}_{-0.012} $                             & $0.086^{+0.012}_{-0.0030}$\\
   $ n_s$ & $0.96^{+0.019}_{-0.0056} $                                                     & $0.99^{+0.029}_{-0.014}$ \\
   $w$ &  $-0.88^{+0.080}_{-0.12}  $                                                                                & $-0.89^{+0.090}_{-0.11}$ \\
   $\beta$ &  $0.012^{+0.050}_{-0.012}$                          & $0.032^{+0.042}_{-0.032}$\\
      $\beta$ &  $< 0.063 (0.11) $                          & $< 0.074 (0.12)$\\
    $ \sigma$ & $0.22^{+0.28}_{-0.090}$                                              & $0.14^{+0.36}_{-0.11}$ \\
   $ N_{eff}$ & -  & 						$3.84^{+0.74}_{-0.49}$  \\

   $ \Omega_{de}$ & $0.72^{+0.076}_{-0.012}$                                  & $0.74^{+0.060}_{-0.031} $ \\
   $ Age/Gyr$ & $13.8^{+0.012}_{-0.39} $                                              & $13.0^{+0.40}_{-0.77}$ \\
   $ z_{re}$ & $10.9^{+0.63}_{-1.8} $                          & $10.7^{+1.7}_{-1.0}$ \\
   $ H_0$ & $68.7^{+8.4}_{-0.98}$                             & $ 75.6^{+9.9}_{-4.3} $  \\
            $ D_{3000}^{SZ}$ & $4.0^{+2.1}_{-4.0}$                    & $6.7^{+2.2}_{-3.3}$ \\
                  $ D_{3000}^{PS}$ & $21.5^{+1.7}_{-3.7}$                    & $20.2^{+3.0}_{-2.5}$  \\
                           $ D_{3000}^{CL}$ & $5.3^{+1.9}_{-2.4} $                    &   $4.4^{+3.2}_{-1.2} $   \\
   -$Log$(Like) &$ 3756$                      							& $3756 $  \\
\hline
\hline
\end{tabular}
\caption{Best fit values and 1-$\sigma$ errors comparing runs  \emph{cq1} and  \emph{cq4}. Both runs include coupling and use WMAP7 + SPT; in addition,  \emph{cq4} marginalizes over $N_{eff}$. For $\beta$ we also write in brackets the value of the 2$\sigma$ marginalized error.}
 \label{tab:cq1cq4_bestfit}
\end{table}
\par\end{center}
\normalsize

 \begin{figure*}
 \centering
    \includegraphics[width=18.cm]{plots_cq/compare_cq1cq2cq4_showfour_andro.eps}%
      \caption{\small Confidence contours for the cosmological parameters for coupled quintessence models. We compare runs  \emph{cq1} (blue),  \emph{cq2} (green) and \emph{cq4} (yellow). The light blue asterisks mark the best fit points for $cq1$ while the pink asterisks mark the best fit points for $cq4$. 1-sigma and 2-sigma contours are shown.}
 \label{fig:like_cont_neff}
 \end{figure*}
 \normalsize

 \begin{figure*}
 \centering
    \includegraphics[width=18.cm]{plots_cq/plot2D_single_neff.eps}%
      \caption{\small Confidence contours for a choice of cosmological parameters vs $N_{eff}$ within coupled quintessence models for run  \emph{cq4}. 1-sigma and 2-sigma contours are shown..  $N_{eff}$ is degenerate with $\Omega_{DM}h^2$, $n_s$, which are in turn degenerate with $\beta$. No direct degeneracy can however be seen between $\beta$ and $N_{eff}$. The pink asterisks mark the best fit for $cq4$.}
 \label{fig:like_cont_neff_yhe}
 \end{figure*}
 \normalsize

\subsection{CMB lensing}

In order to test the effect of CMB-lensing on coupled quintessence, we have
redone a run  \emph{cq1} (WMAP + SPT) without lensing.
The presence of a constant coupling doesn't seem to be very much affected by
lensing in the TT CMB spectra, as we can see in Fig.(\ref{fig:like_cont_NL})
where we compare run  \emph{cq1} with run  \emph{cq1NL}. If no lensing is included, the bound on
$\beta$ is slightly (but not significantly) larger: $\beta < 0.068 (0.13)$.

Furthermore, similarly to \cite{k11} one can rescale the lensing potential power
spectrum by a factor $A_L$:
\begin{equation}
C_{l}^{\phi\phi} \rightarrow A_L C_l^{\phi \phi}
\end{equation}
All runs discussed so far fix $A_L = 1$. In order to test the effect of lensing
we also performed a run ($cq1AL$) in which we vary $A_L$ and marginalize over it. The
$A_L$ parameter was found to be $A_L = 0.94 \pm 0.15$ when using WMAP + SPT with
a $\Lambda$CDM model \cite{reichardt_etal_2011} (see also \cite{2009ApJ...694.1200R,
calabrese_etal_2008, das_etal_2011a, act_2011} for different datasets). When a
coupling between dark matter and dark energy is included, we find that the best
fit for $A_L$ is $A_L = 0.86^{+0.34}_{-0.12}$, still compatible with one. The bound on
the coupling is of the same order as in the case in which $A_L$
is fixed: $\beta < 0.063 (0.11)$. In other words, we don't gain much marginalizing over
$A_L$ instead of fixing it to one, given that $A_L$ best fit is very close and
fully compatible with one. Though with the data considered here the effect is
not significant, $A_L$ is also correlated with dark matter and $n_s$, which in
turn are correlated with $\beta$, as shown in fig.\ref{fig:like_cont_al}.

 \begin{figure*}
 \centering
    \includegraphics[width=18.cm]{plots_cq/compare_cq1_v10_cq1_v10_nolensing_showmore_andro.eps}%
      \caption{\small Likelihood contours for cosmological parameters in presence of a coupling for WMAP+SPT (run  \emph{cq1}, including lensing, red contours), as compared to the same run done without CMB lensing (white contours). Light blue asterisks mark the best fit points for $cq1$.}
 \label{fig:like_cont_NL}
 \end{figure*}
 \normalsize

 \begin{figure*}
 \centering
    \includegraphics[width=15.cm]{plots_cq/plot2D_single_al.eps}
      \caption{\small Likelihood contours for cosmological parameters when $A_L$ is allowed to vary.}
 \label{fig:like_cont_al}
 \end{figure*}
 \normalsize

\subsection{Curvature}
If we release the constraint of a flat universe and allow for curvature and
coupling (run $cq1K$), we get $\Omega_K = -0.0068^{+0.0092}_{-0.036}$, which is compatible with a flat
Universe. In this case, the constraint on $\beta$ is slightly less restrictive,
$\beta < 0.071 (0.13)$ but the bound on $\beta$ is already stringent enough not
to be affected so much by the uncertainty on curvature. Contours are shown in
Fig.\ref{fig:like_cont_curvature} where they are compared to run  \emph{cq1}, in which a flat
universe was assumed. We also show in Fig.\ref{fig:like_cont_omegaK} how
curvature is degenerate, as expected, with the Hubble parameter, though no
direct degeneracy is seen between $\Omega_K$ and the coupling $\beta$.

 \begin{figure*}
 \centering
    \includegraphics[width=18.cm]{plots_cq/compare_cq1_v10_cq1_v10_curvature_showmore_andro.eps}
      \caption{\small Likelihood contours for cosmological parameters of run  \emph{cq1} (red contours) as compared to run  \emph{cq1}, in which $\Omega_K$ is allowed to vary (white contours). Light blue asterisks mark the best fit points of $cq1$ while yellow asterisks mark the best fit points for $cq1K$.}
 \label{fig:like_cont_curvature}
 \end{figure*}
 \normalsize

 \begin{figure*}
 \centering
    \includegraphics[width=15.cm]{plots_cq/plot2D_single_K.eps}
      \caption{\small Likelihood contours for cosmological parameters when $\Omega_K$ is allowed to vary. Yellow asterisks mark best fit points for run $cq1K$.}
 \label{fig:like_cont_omegaK}
 \end{figure*}
 \normalsize

\subsection{Massive neutrinos}
Up to now, we have fixed the relative fraction of massive neutrinos $f_\nu$ to
zero. We now consider run  \emph{cq1}$\nu$ in which we also allow for a non zero
fraction of massive neutrinos and marginalize over $f_\nu$. In this case the
range allowed for the coupling is $\beta < 0.084 (0.14)$, slightly bigger than when
using massless neutrinos, as expected \cite{lavacca2009}. We update the results of \cite{lavacca2009}
 using both WMAP7 and SPT data. The degeneracy between massive neutrinos
and $\beta$ is clearly shown in Fig.\ref{fig:like_cont_fnu}. The best fit value for $f_\nu$ is $f_\nu = 0.065^{+0.017}_{-0.065}$

 \begin{figure*}
 \centering
    \includegraphics[width=15.cm]{plots_cq/plot2D_single_fnu.eps}
      \caption{\small Likelihood contours for cosmological parameters when $f_\nu$ is allowed to vary. Pink asterisks mark best fit points for run $cq1\nu$.}
 \label{fig:like_cont_fnu}
 \end{figure*}
 \normalsize

\subsection{Combining WMAP7 and SPT with HST, BAO and Supernovae Ia data}

As discussed earlier on in this paper and in \cite{amendola_etal_2012}, the
coupling $\beta$ is degenerate with the Hubble parameter. In order to
investigate the effect of this degeneracy on the constraints from data, we did
another run ($cq1hst$) in which we combined the data used for $cq1$ (WMAP7 and
SPT) to also included baryon acoustic oscillations (BAO) \cite{percival_etal_2010}, Hubble Space Telescope constraints on $H_0$ (HST) \cite{riess_etal_2011} and Supernovae Ia (SNae) data \cite{union2_2010} as from COSMOMC (Aug 2011).
Likelihood contours for this case are shown in Fig.{\ref{fig:like_cont_hstv2}}.
Best fit values and errors are shown in Tab.\ref{tab:cq1mock_bestfit}, left
column. There seems to be an interesting preference for a non zero coupling,
though the values are clearly still compatible with zero at 1$\sigma$.
This peak comes mostly from a slight tension between the Hubble parameter HST result ($h=0.738\pm 0.024 $) and our
 WMAP7+SPT best fit for $\beta=0$ ($ h=0.685\pm 0.025 $ ).
Notice however that even for $\beta=0$ we are not in an exact $\Lambda$CDM
since in our model $w$  is close, but not exactly equal, to -1.
It is interesting then to test whether the  forecoming data from Planck can confirm or reject  this non zero coupling.
This we do in the next section.

 \begin{figure*}
 \centering
    \includegraphics[width=15.cm]{plots_cq/plot2D_single_hstv2.eps}
      \caption{\small Likelihood contours for cosmological parameters when data from HST, BAO and Supernovae Ia are included. Blue asterisks mark best fit points for run $cq1hst$.}
 \label{fig:like_cont_hstv2}
 \end{figure*}
 \normalsize

\subsection{Combining Planck and SPT mock data}
As a further analysis, we forecast the effect that Planck data would have on the
coupling parameter, when combined with the power measurement of $SPT$ (run
$cq1Pl$).
Since Planck data are not yet available, we produce a set of mock data
\cite{planck2011}. We have therefore implemented FutureCMB
\cite{perrotto_etal_2006_futurecmb} in our modified version of COSMOMC and,
using as fiducial power spectrum a
$\Lambda$CDM model
$D_l^{th}$, with $D_l = l(l+1)/2\pi  C_l$ (black dot-dashed in
Fig.\ref{fig:mock}).
We have then produced an SPT mock spectrum as follows:
\begin{description}
\item[ ]  we have added to the same fiducial model used to generate Planck mock
data, the effect of poisson sources (PS) (blue dotted) using the value for the
nuisance parameters ${{D_l}^{PS}}_{3000} = 18.1 \mu K^2$ with a dependence from
the momentum $\propto (l/3000)^2$  \cite{k11};
\item[ ] we have added to $D_{l}^{th}$ + PS the effect of clustered sources (CL)
(light blue dashed line), using the value for the nuisance parameters
${{D_l}^{CL}}_{3000} = 3.5 \mu K^2$ with a dependence from the momentum $\propto
(l/3000)^0.8$ for all $l$. This is not entirely correct as the dependence is
slightly different for $l > 1500$ and $l < 1500$ but differences are thought to
be small \cite{k11}; we obtain the $\tilde{D}_l \equiv D_l^{th} + PS + CL$. We
neglect here the effect of SZ.
\item[ ] we have then convolved the $\tilde{D}_l$ for the SPT window functions
and created an SPT mock data with $D_l$ given by the convolved $\tilde{D}_l$ and
errors given by SPT data from \cite{k11}. The final SPT mock data is plotted in
Fig.\ref{fig:mock}) (red dotted line).
\end{description}

 \begin{figure*}
 \centering
    \includegraphics[width=18.cm]{plots_cq/cl_mockspt.eps}%
      \caption{\small SPT mock data as obtained from the fiducial $\Lambda$CDM spectrum.}
 \label{fig:mock}
 \end{figure*}
 \normalsize

\begin{center}
\begin{table}
\begin{tabular}{lllll}
\textbf{Best fit values} & \textbf{for coupled quintessence}\\
\hline
\begin{minipage}{30pt}
\flushleft
Parameter
\\
\end{minipage} &
\begin{minipage}{150pt}
\flushleft
 \emph{cq1hst} (baseline + $\beta$+$\sigma$)\\ real data
\end{minipage} &
\begin{minipage}{150pt}
\flushleft
$cq1Pl$  \\ forecasts
\end{minipage} &
\\
\\
\hline
$  \Omega_b h^2  $ & $0.022^{+00044}_{-0.00040}$           & $0.023^{+0.0001}_{-0.00024} $		\\
$ \Omega_{c} h^2  $ & $0.11^{+0.0022}_{-0.0039}$                                  & $0.11^{+0.0013}_{-0.0012}$ \\
      $ \theta_s$ & $1.04^{+0.0016}_{-0.0014} $                    & $1.05^{+0.00008}_{-0.0006}$ \\
   $ \tau$ & $0.084^{+0.0061}_{-0.0067} $                             & $0.090^{+0.00011}_{-0.0039}$\\
   $ n_s$ & $0.96^{+0.010}_{-0.0086} $                                                     & $0.98^{+0.0024}_{-0.0062} $ \\
   $w$ &  $-0.90^{+0.10}_{-0.098} $                                                                                & $-1.0^{+0.1}_{-0.003}$ \\
   $\beta$ &  $ 0.041^{+0.0092}_{-0.041}$                          & $ 0.0035^{+0.0089}_{-0.0035} $\\
   $\beta$ &  $ < 0.050 (0.074)$                          & $ < 0.012 (0.030) $\\
    $ \sigma$ & $0.19^{+0.31}_{-0.059}$                                              & $0.48^{+0.015}_{-0.11}$ \\

   $ \Omega_{de}$ & $0.74^{+0.012}_{-0.015} $                                  & $0.73^{+0.019}_{-0.000056} $ \\
   $ Age/Gyr$ & $13.7^{+0.082}_{-0.093} $                                              & $13.6^{+0.005}_{-0.08}$ \\
   $ z_{re}$ & $10.3^{+1.1}_{-1.2} $                          & $10.6^{+0.2}_{-0.5} $ \\
   $ H_0$ & $71.0^{+1.0}_{1.7}$                             & $ 69.6^{+2.4}_{-0.13} $  \\
         $ D_{3000}^{SZ}$ & $5.0^{+1.2}_{-5.0} $                    & $2.1^{+2.8}_{-1.1} $ \\
                  $ D_{3000}^{PS}$ & $20.6^{+1.9}_{-3.4} $                    & $15.6^{+2.6}_{-2.3} $ \\
                           $ D_{3000}^{CL}$ & $5.1^{+3.0}_{-1.5} $                    & $2.9^{+1.0}_{-2.8} $ \\
   -$Log$(Like) &$ 4024$                      							& $660 $  \\
\hline
\hline
\end{tabular}
\caption{Best fit values and 1-$\sigma$ errors comparing runs  \emph{cq1hst} and
mock $cq1Pl$ including SPT + PLANCK, plus HST, BAO, SN data. For $\beta$ we also
write in brackets the value of the 2$\sigma$ marginalized error.}
 \label{tab:cq1mock_bestfit}
\end{table}
\par\end{center}
\normalsize

For this case, we also include BAO, HST and Supernovae Ia data as from Cosmomc (August 2011 version),
to break the degeneracy between the coupling and the Hubble parameter. Some representative
2D confidence regions are in Fig. \ref{fig:like_cont_curvature_mock}.
We then find that $\beta < 0.012 (0.030)$; though stronger than WMAP7+SPT, this is still a pessimistic bound,
since it includes Planck but still considers SPT data with errors released by
\cite{k11}; by the time Planck data will be available, better SPT (or ACT) data
may have been released. In Table \ref{tab:cq1mock_bestfit} we report the best
fit values, together with the left and right errors at $68 \%$ and $95\%$ CL
around the best fit. The left column refers to run $cq1hst$, done with real data
(WMAP+SPT+HST+BAO+SNae); the right column shows the forecasted values for
Planck+SPT mock data, plus SNae, HST and BAO.

 \begin{figure*}
 \centering
    \includegraphics[height=12.cm]{plots_cq/plot2D_single_planckspt_try3.eps}%
      \caption{\small Confidence contours for Planck + SPT mock. Blue asterisks mark best fit points for run $cq1Pl$.}
 \label{fig:like_cont_curvature_mock}
 \end{figure*}
 \normalsize

We finally plot in Fig. \ref{fig:finalbeta} a comparison between current data and
future observations,  marginalizing over all parameters except $\beta$.
It is clear  that Planck data will reach a precision sufficient  to tell whether
the peak in $\beta$ is a real detection or just a fluke.

 \begin{figure*}
 \centering
    \includegraphics[width=10.cm]{plots_cq/compare1D_cq1v10_cq1Plv12.eps}%
      \caption{\small 1D likelihood for the coupling $\beta$ from runs $cq1$, $cq1hst$, $cq1Pl$. We recall that $cq1$ and $cq1hst$ are based on real data while $cq1Pl$ estimates the forecasted constraints from mock data, around a fiducial $\Lambda$CDM cosmology.}
 \label{fig:finalbeta}
 \end{figure*}
 \normalsize

\section{Conclusions} \label{conclusions}
We have considered the possibility that the evolution of dark matter and dark energy might be
connected by a constant coupling, of the type illustrated in
\cite{amendola_2000, pettorino_baccigalupi_2008}. We have used current CMB data
from WMAP7 and SPT to constrain the coupling parameter $\beta$. We find that
$\beta$ is constrained to be less than 0.063 (0.11) at $68\%$ ($95\%$) C.L. when SPT data are
included, with respect to $\beta < 0.078 (0.14)$ coming from WMAP7 only.
We have done a number of tests to check whether this bound depends on the degeneracy with other parameters (lensing, curvature, massive neutrinos, $N_{eff},$ HST/BAO/SNae data).
If the effective number of relativistic degrees of freedom $N_{eff}$ is allowed to
vary, no much gain is obtained on $\beta$, which still needs to be $\beta <
0.074 (0.12)$.
We have further considered the effect of CMB-lensing, both with a run which
includes no lensing and by marginalizing over $A_L$, a parameter which encodes
the rescaling of the lensing power spectrum. $A_L$ is slightly degenerate with
$n_s, \Omega_{DM}h^2$ which in turn are degenerate with $\beta$, though no
direct degeneracy is seen between $\beta$ and $A_L$.
If the assumption of a flat universe is released, constraints on $\beta$ weaken
back almost to the level of constraints given by WMAP only (flat universe), with $\beta <
0.071 (0.13)$. Degeneracy with massive neutrinos widens the coupling constraints to be $\beta < 0.084 (0.14)$ when we marginalize over the fraction of massive neutrino species $f_{\nu}$.
We conclude that the bound on $\beta$ from current data is already strong enough to be quite stable with respect to a better knowledge of other parameters and to all cases considered.

When WMAP+SPT are considered (run $cq1$), the best fit value for $\beta$, though
still fully compatible with zero, has a best fit of $\beta =
0.012^{+0.050}_{-0.012}$. It is interesting to see that when we allow for an
effective number of relativistic degrees of freedom, marginalizing over
$N_{eff}$ the coupling from WMAP7+SPT data increases to a best fit value of
$\beta \sim 0.03$. A larger value of $N_{eff}$ favors larger couplings between
dark matter and dark energy and values of the spectral index closer to $1$.
Including SPT data does not improve significantly constraints on the coupling
$\beta$.
Inclusion of additional priors from HST, BAO and SNae moves the best fit to $\beta=0.041$, again still compatible
with zero at 1$\sigma$.
We forecast that the inclusion of Planck data will be able to
pin down the coupling to about $1\%$ and therefore detect whether the small non-zero
coupling present in current data is washed away  with more data.

\begin{acknowledgments}
V.P. is supported by the Marie Curie IEF, Project DEMO - Dark Energy Models and Observations. Support was given to V.P. and C.B. by the Italian Space Agency through the ASI contracts Euclid-IC (I/031/10/0). V.P. thanks Joanna Dunkley for useful discussion and advices. We thank Christof Wetterich for helpful comments.
C.B. also acknowledges support from the PD51 INFN initiative.
L.A. acknowledges support from DFG - Deutsche Forschungsgemeinschaft through the project TRR33 "The Dark Universe".
\end{acknowledgments}