\section{\label{sec:level1}Introduction}
This paper discusses the structural characteristics of random packings of planar (2D) granular solids. The science of random granular packing has a very long history and statistical methods have been introduced into this field in the 60's \cite{Bernal1960, Mogami1965}. In 1989, Edwards and co-workers \cite{EdOa89a,EdOa89b} proposed a novel statistical mechanical approach, based on the entropy of configurations that assemblies of grains can take. In this approach, energy is replaced by volume, the Hamiltonian is replaced by a volume function and the temperature, characterising ensemble fluctuations, is replaced by compactivity. This latter quantity has been elusive to measure and quantify from bulk quantities.
Much work followed, aiming both to use this formalism for derivation of macroscopic quantities, as in thermodynamics, and to test it experimentally and numerically  \cite{Knight-etal1995, Richard-etal2005, Song-etal2008}.

Fundamental to this approach is an identification of basic volume elements, which can be used as 'quasi-particles' for the volume function. Blumenfeld and collaborators have proposed that these quasi-particles are {\it quadrons} \cite{BaBl02,BlEd03,BlEd06} - structural elements that tessellate the granular space and can be used to quantify the local structure \cite{BaBl02,Fretal08,Fretal09}. Ball and Blumenfeld \cite{BaBl02} used the quadron description to construct a stress theory of isostatic granular materials.
In contrast to other definitions of volume elements, such as Volonoi-based tesserations \cite{Song-etal2008}, quadrons preserve the connectivity information, since their construction is based on the force-carrying inter-granular contacts. They also allow for an unambiguous tensorial description of every structural element. A further advantage of this description is that, at least in 2D, there are as many quadrons as there are degrees of freedom in the statistical mechanical structural phase space \cite{BlEd03,BlEd06,Bl08a}. This makes it useful for treating these volumes as degrees of freedom in their own right.

It has been shown \cite{Fretal08} that the quadron volume probability density function (PDF), $P(V_q)$, is best understood in terms of its conditional PDF's $P(V_q\mid e)$, where $e$ is the order of the cell associated with the quadron (see Figure \ref{fig:quadron}). Based on geometrical considerations, and supported by numerical results, Frenkel et al. argued that the PDFs $P(V_q\mid e)$ should be independent of inter-granular friction and of the packing process (i.e. the history).
A detailed analysis of this issue and an understanding of the shape of these PDFs from first principles is yet to emerge.

Here, we present a set of numerical experiments on assemblies of polydisperse discs, isotropically packed by the same process to different packing fractions $\phi$ and mean coordination numbers $\bar{z}$.
We investigate in detail the geometrical characterisation of the quadrons and the dependences  of the PDFs $P(V_q)$ and $P(V_q\mid e)$ on $\phi$ and $\bar{z}$.  We examine the decomposition proposed by Frenkel et al. \cite{Fretal08}, focusing on the significance of cell shapes, and we show that the universality suggested by them must be augmented with a mechanics-based consideration.
This is supported by a free cell model, whose results are compared with the numerical experiments.
We conclude with a discussion of the results and the new insight that they provide into structures of random granular packs in mechanical equilibrium.

\section{Probability density function of quadron volumes in statistical mechanics}

As in thermodynamics, one can define the entropy, $S$, of a micro-canonical ensemble of granular assemblies in static equilibrium as the logarithm of  the number of ways, $M$, that a collection of $N$ grains fit into a given volume $V$,

\begin{eqnarray}
S(V,N)=\log M
\label{S}
\end{eqnarray}
where $M$ can be written as
\begin{eqnarray}
M=\sum_i \delta (V-V_i) G_i \to \int \delta (V-W)\ \Theta\ d^{N_{sdf}}\{\eta \}
\label{M}
\end{eqnarray}
Here $\delta$ is the delta function, $G_i$ is the number of configurations (microstates) in mechanical equlibrium having a total volume $V_i$,
$W$ is the volume function - a function that gives the volume in terms of $N_{sdf}$ internal structural degrees of freedom (SDF), $\{\eta \}$, for any grain configuration. $\Theta$ is a product of delta functions, constraining the possible configurations to those specified in the ensemble statistics, e.g. to assemblies in mechanical equilibrium, with a specified mean coordination number $\bar{z}$. From this micro-canonical description, one can readily pass to the canonical one \cite{Edwards_Rafi2007book}.
Assuming a Gibbs ensemble \cite{EdOa89a,EdOa89b}, the occurrence probability of a configuration of volume $V_i$ is

\begin{equation}
P_i = \frac{1}{Z_v} G_i \exp \left( - \frac{ V_i }{X} \right)
\label{Pi}
\end{equation}
where $X = \left( \partial S / \partial V \right)^{-1}$ is a Lagrange multiplier, called compactivity, that characterizes the statistical fluctuations and is analogous to the temperature.
$Z_v$ is a volumetric partition function

\begin{equation}
Z_v = \sum_i e^{-V_i/X} \to \int e^{-W\left(\left\{ \eta \right\}\right)/X} \Theta(\left\{ \eta \right\}) d^{N_{sdf}} \left\{ \eta \right\}
\label{Z}
\end{equation}

To use the  statistical mechanical formalism, one needs a set of elementary entities, whose summed volumes gives the volume of the system. The quadrons, discussed originally in \cite{BaBl02} and adapted to statistical mechanics in \cite{BlEd03,BlEd06}, are structural elements that fit this role nicely. They can be regarded as the 'quasi-particles' of this formalism\cite{BlEd03,BlEd06}. One also needs to identify the SDF. In terms of these, the volume function is  $\sum_q V_q\left(\left\{ \eta \right\}\right)$. It has been shown \cite{BlEd03,BlEd06}  that $N_{sdf}$ is equal to the number of  quadrons in 2D, but is larger in 3D\cite{Bl12}.

An explicit illustration of the quadron in 2D is shown in Figure \ref{fig:quadron}. First, one defines the centroids of grains $g$ and of cells $c$ as the mean position vectors of the contact points around them, respectively. Second, the contact points around every grain are connected to make polygons, whose edges are vectors, $\vec{r}^{gc}$, that circulate the grain in the clockwise direction. Third, one extends vectors $\vec{R}^{gc}$ from the centroids of grains $g$ to the centroids of the cells $c$ that surround them. A quadron is the quadrilateral whose diagonals are the vectors  $\vec{r}^{gc}$ and  $\vec{R}^{gc}$. Since every pair $gc$ corresponds exactly to one quadron, we shall denote a quadron by the index $q$ for  convenience.
The quadron's shape is  quantified by the following structure tensor

\begin{equation}
C^q = \vec{r}^q\otimes\vec{R}^q
\label{QuadTensor}
\end{equation}
and its volume is
\begin{eqnarray}
V_q= \frac{1}{2}\mid\vec{r}^q \times \vec{R}^q\mid = \frac{1}{2} (r_x^q R_y^q - r_y^q R_x^q)
\label{QuadVolume}
\end{eqnarray}

\begin{figure}[!bht]
\begin{minipage}[t]{0.3\textwidth}
\includegraphics[width=1.0\textwidth]{Fig01_quadron.eps}
\caption{Illustration of quadrons (shaded) and cells (yellow) in 2D.}
\label{fig:quadron}
\end{minipage}
\end{figure}

It has been shown that, in 2D, there are as many quadrons as there are degrees of freedom \cite{BlEd03,BlEd06,Bl08a}, $N_q = N_{sdf}$. This makes it possible to use the quadron volumes to span the $N_q$-dimensional phase space, $\vec{V}$.
We define the `density of states'  $g\left(\vec{V} \right)d^{N_q} \vec{V}=g\left(\left\{ V_q \right\}\right)\prod_q{d V_q}$ as the number of states whose volumes are between $\vec{V}$ and $\vec{V}+d\vec{V}$.

The volume partition function can be then written as

\begin{equation}
Z_v = \int e^{-\sum_{j=1}^{N_q} V_q/X} g(\left\{ V_q \right\}) \prod_{j=1}^{N_q} dV_{q,j}
\label{PartFun}
\end{equation}
where the integration is understood as taken only over the systems that satisfy the constraints in $\Theta(\left\{ \eta \right\})$ and $W=\sum_{q=1}^{N_q} V_q $.

Let us now consider a model system - an ideal quadron gas (IQG) - where the quadrons are uncorrelated. The general partition function (\ref{PartFun}) becomes
\begin{equation}
Z^{iqg}_v = \left[ \int e^{-V_q/X} g_1(V_q) dV_q \right]^{N_q}
\label{PartFunIQG}
\end{equation}
where $g_1(V)$ is a one-quadron density of states \cite{BlEd06}.
The propability density function (PDF) of finding a quadron of volume between $V_q$ and $V_q+dV_q$ is
\begin{equation}
P(V_q) = \frac{1}{Z^{iqg}_v} e^{- V_q /X} g_1(V_q)
\label{PDFQuad}
\end{equation}

In the following, we investigate this PDF in our numerical experiments and relate it to the cell structures. The combination of the exponential term and the density of states, whose exact form is unknown, makes exact calculations impossible without additional information. Nevertheless, it is possible to gain insight from such a study, as we show in the following.

\section{Numerical experiments}

Our numerical experiments have been carried out using the Discrete Element Method (DEM) (\cite{Cundall1979,Mat-Chang-2011}).
The method consists of using an incremental time marching scheme, wherein the motions of 2D discs (grains) are each computed using Newton's second law. We postulate a repelling harmonic interaction potential, characterized by normal and tangential spring constants, $k_n$ and $k_s$, activated upon contact and overlap between discs. We set $k_s/k_n=1/4$ in this study.

Our aim is to study polydisperse systems and we have chosen the distribution of the radii of the discs to be log-normal distribution due to its wide use in civil engineering and soil sciences \cite{Mitchell-Soga-2005},
\begin{eqnarray}
P (D) =\frac{1}{\sqrt{2 \pi} \sigma D} \exp \left( \frac{(\ln D - \ln D_0)^2}{2 \sigma^2}   \right)
\end{eqnarray}
where we assign $D_0=1.$ and $\sigma=0.2$, which give for the mean grain size and mode, respectively, $\bar{D}=1.02$ and $D_{mode}=0.961$ (see Figure \ref{fig:size_distrib}).

\begin{figure}[!bht]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig02_size_distrib.eps}
\caption{Grain diameter probability density function.}
\label{fig:size_distrib}
\end{minipage}
\end{figure}

The packing protocol of our systems is as follows.
First, we construct three random packs in a double periodic domain, engineered to be on the verge of jamming. The packs consist of about $21400 \pm 1000$ discs and are made at packing fractions $\phi=0.76, 0.82$ and $0.84$. The small variation in the numbers of particles is due to the different densities of the three initial configurations.
These configurations are then used as initial states for the packing procedure; loose initial state (LIS), intermediated initial states (IIS), and dense initial state (DIS), respectively.

Once an initial state is set, we assign all the particles the same friction coefficient, $\mu$, and apply to the system a slow isotropic stress $\sigma_c$ by changing the periodic length on both sides.
We limit the applied stress to a level that corresponds to an average overlap discs of $\delta=\sigma_c / k_n =10^{-5}$.
No gravitational force is applied and the compression continues until the fluctuations of both grain positions (per mean grain diameter) and inter-granular forces (per mean average contact force) are below very small thresholds - $10^{-9}$ and $10^{-6}$, respectively. This procedure is carried out for each initial state at five different values of the inter-granular friction coefficients: $\mu=0.01, 0.1, 0.2, 0.5$ and $10$, giving altogether 15 assemblies.

Using this procedure, we have computed the packing fractions, the mean coordination numbers, and studied the structures in these systems.
For the determination of $\bar{z}$, we disregard `rattlers', i.e. grains with one or no force-carrying contact. The results are shown in Figure \ref{fig:packing_fraction}.
The upper and lower bounds of the coordination number $\bar{z}_{max}=4$  and $\bar{z}_{min}=3$ correspond to the isostatic states for discs with friction coefficients of $\mu=0$ and $\infty$, respectively.
The three systems with $\mu=0.01$ converge into states with ($\bar{z}$, $\phi$)=(3.941, 0.840),(3.944, 0.842) and (3.943, 0.843), respectively, which are very close to the ideal frictionless jammed state. However, with increasing friction coefficient, the difference between the values of $\bar{z}$ and $\phi$ in the final states of systems started from different initial states increases (see Figure  \ref{fig:packing_fraction}).

\begin{figure}["here"]
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=1.0\textwidth]{Fig03_zg-phi.eps}
\caption{Coordination number vs. packing fraction for 15 systems, generated from three different intial states, LIS, IIS, and DIS. Note the convergence for $\mu\to 0$.}
\label{fig:packing_fraction}
\end{minipage}
\end{figure}

Turning to the analysis of cell structures in these systems, consider the three examples shown in Figures \ref{fig:cell_structure1} - \ref{fig:cell_structure3}, constructed from LIS, with $\mu=10$ (\ref{fig:cell_structure1}) and $0.01$ (\ref{fig:cell_structure2}), and from DIS with $\mu=10$ (\ref{fig:cell_structure3}).
While traces of the initial condition can be detected in both the first two systems, for higher friction one clearly ends up with typically larger cells and a correspondingly lower $\bar{z}$.
Note that the packing fraction in Figures \ref{fig:cell_structure1} and \ref{fig:cell_structure3} are very similar (see also Figure \ref{fig:packing_fraction}), but the cells in  \ref{fig:cell_structure1} are noticeably bigger typically.
Also note the considerable number of rattlers in the large cells, an issue that will be discussed in detail below.
To quantify the structural differences, we next study the quadron volume distribution.

\begin{figure}[!bht]
\begin{minipage}[t]{0.3\textwidth}
\includegraphics[width=1.0\textwidth]{Fig04_cell_loose_mu10.eps}
\caption{Example of an assembly with $\mu=10$, generated from LIS.}
\label{fig:cell_structure1}
\end{minipage}
\end{figure}

\begin{figure}[!bht]
\begin{minipage}[t]{0.3\textwidth}
\includegraphics[width=1.0\textwidth]{Fig05_cell_loose_mu01.eps}
\caption{Example of an assembly with $\mu=0.01$, generated from LIS.}
\label{fig:cell_structure2}
\end{minipage}
\end{figure}

\begin{figure}[!bht]
\begin{minipage}[t]{0.3\textwidth}
\includegraphics[width=1.0\textwidth]{Fig06_cell_dense_mu10.eps}
\caption{Example of an assembly with $\mu=10$, generated from DIS.}
\label{fig:cell_structure3}
\end{minipage}
\end{figure}

\section{Structural characteristics}

\subsection{\label{sec4-1} Quadron volume statistics}

The grain size distribution is the same for all packs, making it convenient to normalize the quadron volumes by the mean grain volume $\bar{V}_g$, $v\equiv V_q/ \bar{V}_g$. The overall PDFs of quadron volumes are plotted in Figure \ref{fig:PDF of quadron} for all friction coefficients and for LIS and DIS.
The corresponding means are plotted against the mean coordination numbers, for all initial states, in Figure \ref{fig:mean quadron volume vs zg}. All the points fall nicely on one curve {\it regardless of initial state}.

\begin{figure}["here"]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig07_vqvg_all.eps}
\caption{PDF of the quadron volumes, normalized by the mean grain volume, $v=V_q/\bar{V}_g$ for all the 10 systems generated from LIS and DIS.}
\label{fig:PDF of quadron}
\end{minipage}
\end{figure}

\begin{figure}["here"]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig08_vqvg-zg.eps}
\caption{Mean quadron volume $\bar{v}$ vs. mean coordination number $\bar{z}$ for all initial states.}
\label{fig:mean quadron volume vs zg}
\end{minipage}
\end{figure}

On close inspection, this very weak dependence, if any, on initial state can also be observed in Figure \ref{fig:PDF of quadron}.
This implies that the dependence of the relation between $\bar{v}$ and $\bar{z}$ on the initial state, seen in Figure  \ref{fig:packing_fraction}, can be made to disappear under the right parameterization. A clue to such a parameterization is the fact that the quadron volumes are computed {\it ignoring the rattlers}.
Following this idea, we can derive the relation between $\bar{v}$ and $\bar{z}$ as follows

\begin{eqnarray}
\phi'=\frac{V_s'}{V}=\frac{\sum^{N_g'} V_g}{\sum^{N_q} V_q}
=\frac{N_g' \bar{V}_g}{N_q \bar{V}_q}
=\frac{1}{\bar{v} \bar{z}}
\label{eq:model03}
\end{eqnarray}

\begin{eqnarray}
\bar{v} = \frac{1}{ \phi' \bar{z} }= \frac{1}{(1-\eta_{r}) \phi \bar{z}}
\label{eq:mod_phi}
\end{eqnarray}
where $\eta_{r}=N_{r}/N$ is the rattlers fraction and $\phi'$ is the packing fraction after removing all the rattlers.
This elegant result suggests that there is indeed a unique relationship between $\bar{z}$ and the packing fraction, but only if the latter is corrected for the rattler fraction.
This is verified by a direct plot of these two quantities, Figure \ref{fig:mean zg vs. mod phi}.
It is also important to note that this relation is independent of initial state as Figure \ref{fig:mean zg vs. mod phi} shows.

\begin{figure}["here"]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig09_zg-phi_mod.eps}
\caption{Mean coordination number $\bar{z}$ vs. the rattlers-free packing fraction $\phi'$.}
\label{fig:mean zg vs. mod phi}
\end{minipage}
\end{figure}

Comparing Figures \ref{fig:mean zg vs. mod phi} and \ref{fig:packing_fraction}, we see that the rattler-free packing fractions are much smaller than the packing fractions commonly reported in the literature.
The removal of rattlers also makes sense on mechanical grounds, since it does not affect the force transmission in our systems (nor in any system experiencing no body forces).

We conclude that the seeming sensitivity of $\bar{z}(\phi)$ to the initial state (e.g. Figure \ref{fig:packing_fraction}), commonly seen in the literature,  stems directly from the variation in $\eta_{r}$. Plotting $\eta_{r}$ as a function of $\bar{z}$ in our systems (Figure \ref{fig:r_rattlers vs. mean zg}), we see that it is the small differences between the curves that gives rise to the observed differences in the conventional plots.

\begin{figure}["here"]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig10_Rrattlers-zg.eps}
\caption{ The rattlers fraction $\eta_{r}$ vs. the mean coordination number $\bar{z}$.}
\label{fig:r_rattlers vs. mean zg}
\end{minipage}
\end{figure}

\begin{figure}["here"]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig11_vqvq_pdf.eps}
\caption{The PDF of the normalized quadron volume $u=V_q/\bar{V}_q$ for  for all the 10 systems generated from LIS and DIS.}
\label{fig:vq_PDF_normalized_by_mean_vq}
\end{minipage}
\end{figure}

It is important to comment here that this result is not really universal. It is correct for the specific packing protocol used here. While we believe that, for any given protocol, the plots of $\bar{z}(\phi')$ would collapse onto one curve, there is no reason that this curve should be universal.  In other words, we expect different protocols to display different $\bar{z}(\phi')$ curves, probably depending on the rate of rattlers generation.
However, these results do provide a universal protocol-independent insight - the relation between the $\bar{z}$ and the packing fraction is directly linked to the mechanical stability of the structure and the way that forces are transmitted. We will discuss this insight more in the concluding section.

Turning to consider the PDFs of the quadron volumes more closely, it is natural to expect the mean quadron volume to increase with $\mu$, simply because the cells get bigger. To get insight into the shapes of the PDFs, we scale the quadron volumes by their means, $u\equiv  V_q / \bar{V}_q=v/\bar{v}$. This simple scaling is sufficient to collapse all the PDFs almost perfectly onto one curve for all the systems, independent of friction and initial state (Figure \ref{fig:vq_PDF_normalized_by_mean_vq}).

Moreover, the collapsed PDF has an exponential tail, which  is a signature of the Boltzmann-like factor of equation \ref{PDFQuad}.
Our best fit to the collapsed curve has the following form
\begin{eqnarray}
P(u) = f(u) e^{\alpha u}
\end{eqnarray}
where $\alpha$ is a slope of the tail and $f(u)$ is a rational function. Comparing this expression with equation \ref{PDFQuad}, we identify
\begin{eqnarray}
X=\bar{V}_q / \alpha
\end{eqnarray}
where we used $u=v/\bar{v}=V_q/\bar{V}_q$.
This implies that the compactivity $X$ is only a function of $\bar{V}_q$, in agreement with the equipartition principle derived recently in \cite{Bl12}.
Intriguingly, combining this observation (or the equipartition principle of \cite{Bl12}) with equation (\ref{eq:mod_phi}), we find a expression relating directly the mean coordination number $\bar{z}$, the compactivity $X$ and the rattler-free packing fraction $\phi'$,

\begin{equation}
\frac{1}{ \phi' \bar{z} X }= \frac{\alpha}{\bar{V}_g}
\label{eq:EOS}
\end{equation}
where the right hand side is a constant that depends only on the grain size distribution.

Thus, the quadron description makes it possible to collapse the statistics of all the systems, given the correct normalisation. As such, it gives better insight into the general characteristics of granular packs - characteristics that are independent of both the inter-granular friction coefficient and the initial state.
In the next section, we explore the reasons for the apparent unified nature of the quadron volumes statistics.

\subsection{\label{sec4-2} Cell order statistics}

To understand better the structural characteristics, we use a recently-proposed decomposition of the quadron volumes into conditional distributions \cite{Fretal08}

\begin{eqnarray}
P\left(v\right)=\sum_e Q(e) P\left(v\mid e\right)
\label{eq_decomposition}
\end{eqnarray}
\\
Here $Q(e)$ is the occurrence probability of cells of order $e$ and $P(v|e)$ is the conditional PDF of the normalized quadron volume, given that it belongs to a cell of order $e$.
$Q(e)$ is essential to the understanding of random granular packing \cite{BlToMaSoon} and it is this PDF that we wish to focus on next.
As expected from Euler topological relation (below) and as can be observed from Figures \ref{fig:cell_structure1}-\ref{fig:cell_structure3}, a lower value of $\bar{z}$ must be accompanied with a higher mean value of $\bar{e}$ and corresponds to a lower packing fraction $\phi'$.

\begin{figure}["here"]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig12_q_e.eps}
\caption{The cell order probability $Q(e)$ for  for the 10 systems generated from LIS and DIS.}
\label{fig:q(e)_dense}
\end{minipage}
\end{figure}

In Figure \ref{fig:q(e)_dense} we plot $Q(e)$ for all the systems generated from LIS and DIS. The plots make evident two points. One is that the higher the inter-granular friction the larger the fraction of high-order cells. The other is  that $Q(e)$ is hardly dependent on the initial state for any $\mu$.

There is a direct relation between the mean cell order, $\bar{e}$, and the mean coordination number $\bar{z}$ and it can be derived from Euler's topological relation for a planar graph,

\begin{equation}
N_V - N_E +N_C = 1
\label{Euler}
\end{equation}
In this relation, $N_V$, $N_E$ and $N_C$ are, respectively, the numbers of the graph's vertices, edges and cells. In this relation we disregard the one large cell making the outside of the system, whose inclusion gives a topological characteristic 2, rather than 1, on the right hand side. This topological characteristic is negligible for sufficiently large systems.
There are several ways to choose the graph in granular assemblies. One is to take the grain centres as vertices and the lines connecting touching grains as edges. This construction is convenient for discs, but can be awkward for grains of arbitrary shapes. Another, more generally applicable, construction is to take the contact points as vertices and the lines connecting contact points around grains as edges. This graph is consistent with the definition of the quadrons (Figure \ref{fig:quadron}).
Both graphs give the same result for packs of $N(\gg1)$ grains

\begin{eqnarray}
\bar{e} = \frac{2 \bar{z}}{ \bar{z} -2} + O\left(\frac{1}{\sqrt{N}}\right)
\label{eq:mean_e}
\end{eqnarray}
where the rightmost term is a (negligible) boundary correction. Figure \ref{fig:Euler's relation} shows that the numerical systems satisfy this relation very well, showing that the sizes of the numerical pack are sufficiently large to ignore finite size effects.

\begin{figure}["here"]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig13_zg-e_relation.eps}
\caption{The mean cell order $\bar{e}$ vs. the mean coordination number $\bar{z}$.}
\label{fig:Euler's relation}
\end{minipage}
\end{figure}

$Q(e)$ is sensitive to the value of $\bar{z}$, which is a function of $\mu$ (see Figure \ref{fig:q(e)_dense}), while $P\left(v\mid e\right)$ has been shown by Frenkel et al \cite{Fretal08}  (see also below) to be hardly dependent on $\mu$. This suggests that there may be a parameterization of $Q(e)$ that collapses all the curves corresponding to the collapse of $P\left(u^q \right)$ in Figure \ref{fig:vq_PDF_normalized_by_mean_vq}.

We find that all the curves collapse if we plot $\bar{e} e Q(e)$ as a function of $e'=(e-\bar{e})/\bar{e}^2$.
We note that the integrals of the curves in Figure \ref{fig:normalized_qe} are exactly 1. The collapsed curve appears to be fitted nicely by a truncated Gaussian

\begin{eqnarray}
\bar{e} e Q(e) =
\frac{\sqrt{2/\pi}}{\sigma_r (1+erf(e'_{min}/\sqrt{2}\sigma_r))}
exp\left(-e'^2/2\sigma_r^2\right)
\label{eq:gaussian}
\end{eqnarray}
with $erf(x)$ being the error function, $e'_{min}=0.076$ and $\sigma_r=0.082$. However, we note small deviation from the Gaussian form at the large $e$-tail, for which we currently have no explanation.

\begin{figure}[!bht]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig14_normalized_qe_gauss.eps}
\caption{The normalized cell probability $\bar{e}e Q(e)$. It is fitted well by a Gaussian form, except for small, but consistent, deviations in the large-$e$ tail (highlighted by a dash lined ellipse).}
\label{fig:normalized_qe}
\end{minipage}
\end{figure}

\subsection{\label{sec4-3} Conditional quadron volume distributions}

To gain further insight into the universal properties of the structure, let us consider in more detail the conditional PDFs $P\left(v\mid e\right)$.
These PDFs were studied by Frenkel et al. \cite{Fretal08}, who argued, on the basis of geometrical considerations, that they should be independent of inter-granular friction. Their argument was based on the observation that, given a collection of $N$ arbitrary grains, the number of ways to arrange $e$ grains into a cell of order $e$  depends only on the grains shapes and not on the inter-granular  friction.
While our results seem to provide a support to this argument, a closer look shows a systematic $\mu$-dependence of $P\left(v\mid e\right)$ for $e> 6$ (Figures \ref{fig:cPDF_e5} and \ref{fig:cPDF_e6}). We can also see this effect in the behaviour of the mean of the conditional quadron volume as a function of $e$, $\bar{v}(e)=\sum_{q\in e}v P\left(v\mid e\right)$, shown in Figure \ref{fig:Mean Vq by cell order}. We include in the figure calculated values of quadron volumes of regular polygonal cells (RPC) of order $e$,

\begin{eqnarray}
v^{RPC}(e)=\frac{V_q^{RPC}(e)}{\bar{V}_g}
=\frac{1}{\pi} cot \frac{\pi}{e}
\label{eq:regular_polygon}
\end{eqnarray}
which is clearly an upper bound for $\bar{v}(e)$.
It is constructive to consider the ratio of the observed mean quadron volume to the regular polygon value, $\gamma(e)\equiv\bar{V}_q(e)/\bar{V}_q^{RPC}(e)$, shown in Figure \ref{fig:vq_vrp}. We observe that $\gamma(e)$, which is always below unity, has a minimum between $e=5$ and $e=6$ and that it decreases with $\mu$. These can be understood as follows. Carrying out a cell shape analysis (to be reported elsewhere), where cells were approximated to lowest order as ellipses, there is a strong correlation between $\gamma(e)$ and the mean aspect ratio of the ellipses for any given $e$.
The implication is that the increase of $\gamma(e)$ with $\mu$ is due to the decrease of mechanical stability of elongated cells as $\mu$ decreases.
The behaviour for any particular $\mu$ depends then on two effects. As $e$ increases, more cell configurations can be realised, allowing for more elongated cells and this accounts for the initial decrease of $\gamma(e)$ with $e$.
However, the more elongated the cell, the less stable it is, leading to fewer elongated cells as $e$ increases.
The competition between these two effects gives rise to the observed minimum.
To test this understanding, we use next an analytical model.

\begin{figure}[!bht]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig15_conditional_vq1.eps}
\caption{The conditional PDF of the quadron volumes for $e \le 5$ collapse nicely for all the systems.}
\label{fig:cPDF_e5}
\end{minipage}
\end{figure}
\begin{figure}[!bht]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig16_conditional_vq2.eps}
\caption{The conditional PDF of the quadron volumes for $e = 6, 8$ collapse, but not as sharply as those for $e \le 5$ (Figure \ref{fig:cPDF_e5}).}
\label{fig:cPDF_e6}
\end{minipage}
\end{figure}
\begin{figure}[!bht]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig17_mean_vq-mean_vg_for_cell.eps}
\caption{The variation of the mean conditional quadron volume, $\bar{v}$, with cell order $e$. The quadron volume of the regular polygon cell  (dashed line) form an upper bound.}
\label{fig:Mean Vq by cell order}
\end{minipage}
\end{figure}

\begin{figure}[!bht]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig18_mean_vq_per_cell.eps}
\caption{The variation of the ratio of the mean quadron volume to regular polygon volume, $\gamma$ ($<1$), as a function of cell order $e$.}
\label{fig:vq_vrp}
\end{minipage}
\end{figure}

\section{A free cell model}

We believe that the geometry-based argument of Frenkel et al \cite{Fretal08}, mentioned above, should be augmented with stability considerations. Indeed, the number of arrangements of $e$ grains, chosen from a given collection of $N$ grains, is independent of $\mu$ if we care only about the geometry. However, not all these shapes are stable under all possible combinations of compressive forces and it is this $\mu$-dependent stability that constraints the possible shapes. We demonstrate this effect, using a 'free cell model' (FCM).
The model consists of a single cell, surrounded by $e$ discs, chosen from the same distribution of sizes as in the numerical simulation. We construct a large number of cells for different values of $\mu$ and analyse their shape distributions.
The cell construction is as follows.
First, a random set of $e$ discs is generated from the size distribution.
We place the first disc at the origin and and another one touching it along the $x$-axis.
We then generate a random angle ($\theta_2$ in Figure \ref{fig:free_cell}) and place the third disc in contact with the second, making sure that it does not touch the first disc.
This process is repeated, making sure that each added disc makes contact only with the previous one,  until we reach the penultimate, $(e-1)$th, disc.
We then determine the size and position of the last disc. If it fits in the remaining opening  then the cell is kept, otherwise it is discarded and we attempt to generate a new cell.
Figure \ref{fig:six_cells} shows examples of free cells, generated by this procedure.
The FCM allows us to generate quickly and efficiently a large number of cells and analyse them statistically.
By comparing this geometry-based model with the DEM results, we expect to identify the effects of mechanical stability for high order cells.

A comparison between the resulting PDFs and those from DEM simulation for all $e=3-8$ are shown in Figure \ref{fig:pdf_vq_freecell}.
The PDFs are similar, but with consistent differences.
The mean $v$ for different $e$ (Figure \ref{fig:vq_vrp})  are lower than of the numerical ones for $e>6$.
This systematic difference is the direct consequence of the cell stability constraint, which has not been considered by Frenkel et al. The DEM cells are closer to regular polygons because they must be stable under compressive forces. The more frictional the grains, the more elongated cells that they can support. This is the reason for the minimum in $\gamma(e)$  decreasing and shifting to higher $e$ with increasing $\mu$.

\begin{figure}[!bht]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig23_free_cell.eps}
\caption{The construction of a cell configuration studied by the free cell model.}
\label{fig:free_cell}
\end{minipage}
\end{figure}

\begin{figure}[!bht]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=0.8\textwidth]{Fig25_six_cells_s.eps}
\caption{ Examples of free cell configurations for $e=6$.}
\label{fig:six_cells}
\end{minipage}
\end{figure}

\begin{figure}[!bht]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig28_vq_free_cell.eps}
\caption{A comparison of the PDF of quadron volumes, generated both with the free cells model and with the DEM numerical experiments.}
\label{fig:pdf_vq_freecell}
\end{minipage}
\end{figure}

\section{Conclusions and Discussions}

To conclude, we have studied in detail the effects of inter-granular friction and initial conditions on the structural characteristics of 2D granular assemblies. We have shown that structure can be analysed quantitatively, using the quadron description\cite{BaBl02,BlEd03}. In particular, we have established that a number of distributions collapse onto single curves, pointing to a `universal'-like structural characteristics that are independent not only of initial state, but also of the friction coefficient. These are significant results, whose implications are discussed at the end of this section.

The raw quadron volume distribution, which is key for the statistical mechanics of granular matter\cite{BlEd03,BlEd06,HiBl12}, changes systematically towards larger values with increasing friction coefficient $\mu$ and hence with decreasing mean coordination number $\bar{z}$. However, when normalised by the mean quadron volume, it collapses nicely onto a single curve. Following the insight proposed by Frenkel et al.\cite{Fretal08}, we have traced this to a similar collapse of the conditional PDFs $P(V_q/\bar{V}_g\mid e)$ for every cell order $e$, except for small, but insightful, deviations for large $e$, which are discussed below.

The common wisdom in the soil mechanics literature is that the relation between the the mean coordination number, $\bar{z}$, and the packing fraction $\phi$ depends on the initial state, as shown in Figure \ref{fig:packing_fraction}.
However, prompted by the observation that the collapse appears once the rattlers are disregarded in the quadron partition of the granular space, we have examined this relation as a function of rattlers-free packing fraction $\phi'$. We discover that, with this correction, the dependence on the initial states disappears and all the relations collapse onto a single curve, $\bar{z}(\phi')$.
This is consistent with our observation that the plot of $\bar{V}_q$ as a function of $\bar{z}$ is also independent of the initial state.
Since rattlers do not support forces in our systems, these results carry a fundamental significance - they suggest that the packing fraction, the mean coordination number, the quadron volume distribution and generally the structure are linked directly to the manner in which stresses are transmitted.
While this gives insight into the structure of the system, it should be remembered that it is difficult to estimate the packing fraction of rattlers in real experiments and it is unclear at this stage how to use this insight for engineering applications.

Considerations of mechanical stability also suggest that the argument proposed by Frenkel at el.\cite{Fretal08}, leading to independence of the conditional quadron volume distributions of inter-granular friction, needs to be modified. Their argument assumes that cell shapes are independent of mechanics and therefore that arranging an $e$-sided cell depends only on the grain shape distribution. This, in turn, eliminates friction as a factor in the occurence probability of cell structures. However, as friction decreases,  elongated cells are less stable mechanically, reducing the number of available cell shapes. Indeed, we can observe a small increase in the occurrence probability of large quadron volumes in $P(V_q/\bar{V}_g\mid e)$ for $e>6$ as $\mu$ increases (Figure \ref{fig:cPDF_e6}).

This effect is seen most clearly when plotting the ratio of the mean quadron volume to that of a regular polygon as a function of $e$ (Figure \ref{fig:vq_vrp}). The initial drop reflects the departure of the cell from a regular polygon due to the increased number of possible geometrical configurations, while the subsequent increase shows the effect of mechanical stability in limiting such a departure. The more pronounced increase for lower $\mu$ is evidence for our argument above.

To  illustrate the effect of mechanical stability, we have constructed a free cell model (FCM) that takes into account only geometric considerations and disregards mechanical stability. Indeed, this model shows the same initial decrease of the above ratio with $e$ and none of the subsequent increase (Figure \ref{fig:vq_vrp}).
Thus, neglecting the mechanical stability effect, the FCM provides a low bound for the mean quadron volume, as illustrated in Figure \ref{fig:VqvsZg}.

\begin{figure}["here"]
\begin{minipage}[t]{0.4\textwidth}
\includegraphics[width=1.0\textwidth]{Fig30_discussion.eps}
\caption{Simple prediction of the mean quadron volume in terms of the mean coordination number.}
\label{fig:VqvsZg}
\end{minipage}
\end{figure}

We have further found that the raw cell order PDF, $Q(e)$, changes systematically with $\mu$, but hardly at all with the initial state.
However, all the PDFs also collapse to a single curve under the  transformations $Q(e)\to\bar{e}eQ(e)$ and $e\to\left(e-\bar{e}\right)/\bar{e}^2$.
We note that the collapse cannot be perfect at the large $e$ tail due to the aforementioned constraint of mechanical stability. This friction-specific behaviour, which is evidenced by Figures \ref{fig:cPDF_e6} and \ref{fig:vq_vrp}, is hardly noticed at all in Figure \ref{fig:normalized_qe}, demonstrating that the effect of mechanical stability is small compared to the geometrical effect discussed by Frenkel et al. \cite{Fretal08}.

Our results are significant for several reasons.
First, they are a step towards a systematic understanding of structural organisation of granular matter in response to a specific packing procedure. In particular, they suggest that the structural characteristics of granular matter are best understood via the statistics of cell configurations.
Second, they make it possible to progress on the statistical mechanics of granular matter, where knowledge of the quadron volumes distribution is key.
Third, the results shed light on the common wisdom in the soil mechanics and civil engineering communities that the initial state affects the final relation between $\bar{z}$ and $\phi$.

Our findings show that this relation is directly linked to mechanical stability of the packing backbone - for any specific pack generation process, the initial state affects the fraction of rattlers, which do not participate in the stress transmission. Once the rattlers are disregarded, there emerges a unique relation between $\bar{z}$ and the corrected rattlers-free packing fraction $\phi'$.
Fourth, our analysis has a significant ramification for the famous packing problem. Different packing procedures lead to different structures. For example, procedures not constrained by mechanical stability would give rise to different $Q(e)$'s than those that are. Specifically, the former would tend to have lower packing fractions due to a higher occurrence probability of large cells, which are mechanically unstable. This reinforces the need for a more accurate definition of the packing problem, which we discuss elsewhere\cite{BlToMaSoon}.

Taken together, these results show that the packing process has a fundamental effect on the structural characteristics while the initial conditions and the inter-granular friction are details that only modify this effect and can be scaled away. The collapsed distributions show how to predict such modifications.
It should be emphasised that the collapses, which we have found, do not imply universality in the traditional sense because the specific forms of the collapsed curves depend on the packing procedure.
For example, it has been demonstrated that, using a different packing protocol, it is possible to generate packs dominated by cells of order 3 \cite{HiBl12}, which cannot be scaled to match the ones we obtained here.  There is in principle an infinite number of possible procedures to generate packs, depending on an astronomically large number of parameters. It would be impossible to map all the possible processes to a manageable parameter space. Each such procedure gives rise to its own characteristic $Q(e)$ and it would be impossible to collapse all of these $Q(e)$ onto one curve.

To fully understand how structural characteristics depend on the packing procedure, one needs a good model of both the geometric effects and the limitations that mechanical stability poses on the cell shape distribution. We are currently developing such a model and will report results in a later publication.
Intriguingly, the relevance of both volumetric and stress effects is reminiscent of the inter-dependence demonstrated recently for the statistical mechanical understanding of granular matter \cite{Bl12} and it would be interesting to find out whether the two issues are related or not.

\nocite{*}


\section{Introduction}
Sensitive radio observations have revealed diffuse emission from the
central regions of some merging galaxy clusters
(e.g. Giovannini et al. 1999, Giovannini \& Feretti 2000,
Kempner \& Sarazin 2001, Govoni et al. 2001a, Bacchi et al. 2003,
Venturi et al. 2007, Venturi et al. 2008, van Weeren et al. 2009,
Rudnick \& Lemmerman 2009, Giovannini et al. 2009, van Weeren et al. 2011).
These radio sources, extended over volumes
of $\sim$1 Mpc$^3$ and called radio halos, are diffuse,
low-surface-brightness ($\simeq$ 1$\mu$Jy~arcsec$^{-2}$ at 1.4 GHz),
and steep-spectrum\footnote{$S(\nu)\propto \nu^{- \alpha}$, with
$\alpha$=spectral index} ($\alpha>1$) synchrotron sources with
no obvious optical counterparts (see e.g. Ferrari et al. 2008,
Feretti et al. 2012).

Radio halos can have quite different length scales, but the largest
halos are the most powerful, in such a way that all these sources may
have similar synchrotron emissivity (Murgia et al. 2009).
Their radio power at 1.4 GHz correlates with the cluster X-ray
luminosity (Feretti 2002), temperature (Liang 1999, Colafrancesco 1999),
mass (Govoni et al. 2001a), and Sunyaev-Zel'dovich effect (Basu 2012)
measurements.

The radio halo morphology is often very similar to the X-ray-emitting
thermal intra-cluster medium (e.g. Govoni et al. 2001b,
Feretti et al. 2001, Keshet \& Loeb 2010, Brown \& Rudnick 2011).
Furthermore, radio halos are preferentially found in clusters that
show evidence of merger activity (e.g. Buote 2001, Schuecker et al. 2001,
Govoni et al. 2004, Cassano et al. 2010, Rossetti et al. 2011), which
suggests a connection between the origin of radio halos
and gravitational processes of cluster formation,
although a one-to-one association between
merging clusters and radio halos is not supported by present observations.

Detailed images of radio halos can provide information on the
cluster magnetic field since the halo brightness fluctuations and the
polarization level are strictly related to the intra-cluster magnetic
field power spectrum (Tribble 1991, Murgia et al. 2004, Govoni et al. 2006).
Vacca et al. (2010) presented a study of the magnetic field power
spectrum in the galaxy cluster A665, which contains a Mpc-scale radio
halo. Their modeling suggests that radio
halos can be effectively polarized, but because of their faintness,
detecting this polarized signal is a very hard task with the current
radio interferometers.  Indeed, radio halos are usually found to be
unpolarized.  Polarized emission from radio halos has been observed so
far only in filaments of the clusters A2255
(Govoni et al. 2005, see also Pizzo et al. 2011)
and MACS J0717+3745 (Bonafede et al. 2009).

The statistics of radio halos remain poor because of their faintness.
To date, only $\sim$40 radio halos are known (see e.g. the recent
compilation by Feretti et al. 2012 and
references therein), thus the discovery of new halos
provide important insights on this class of sources.  As part of an
ongoing program aimed to investigate, in complex X-ray cluster
systems, the presence of halo emission from the sub-clusters in the
merger, we recently discovered the case of the double
merging system, A399
and A401, which both contain a radio halo and can be considered the
only case so far of a double-radio halo system (Murgia et al. 2010).
We also detected an intriguing diffuse radio emission in A781 (Govoni et
al. 2011), an exceptional system characterized by a complex of several
clusters. The archival VLA data set containing the observation of
A781 includes uniform observations for several galaxy clusters.
Thus, to improve the scanty statistics of radio halos, we analyzed the entire
data set and the results of this investigation, complemented by X-ray,
optical, and higher resolution radio data are presented in
this work.

The intrinsic parameters quoted in this paper are computed for a
$\Lambda$CDM cosmology with $H_0$ = 71 km s$^{-1}$Mpc$^{-1}$,
$\Omega_m$ = 0.27, and $\Omega_{\Lambda}$ = 0.73.
Values taken from the literature have been scaled to this cosmology.

\section{VLA Observations and data reduction}

We have analyzed archival observations taken with the VLA at 1.4 GHz
in D configuration (VLA Program AM0469).
This data set contains 41
galaxy clusters observed on average for 15 minutes each.
A brief discussion of the quality of these data is given in the Appendix.

The data were calibrated in phase and amplitude.
Data editing has been made in order to excise RFI.
Since the observation was done in spectral line mode,
data were calibrated in bandpass and data editing
has been made channel by channel.
Surface brightness images were
produced following the standard procedures: Fourier-Transform, Clean,
and Restore implemented in the AIPS task IMAGR.
We used the Multi-scale CLEAN (see e.g. Greisen et al. 2009),
an extension of the classical Clean algorithm,
implemented in the task IMAGR.
We averaged the 2 IFs (and the 7
channels) together in the gridding process under IMAGR.
Self-calibration (phase) was performed to increase the dynamic range
and the sensitivity of the radio images.
The spectral line mode of this data set is not suitable to produce
polarization sensitive images.

We detected evidence of diffuse radio emission in seven clusters
already known to host a radio halo:
A1758, A851, A1995 (Giovannini et al. 2009);
A1351 (Giacintucci et al. 2009, Giovannini et al. 2009);
A773, A2254 (Govoni et al. 2001a);
A2219 (Bacchi et al. 2003, Orr{\'u} et al. 2007).

In addition, we discovered new diffuse sources in
A800, A910, A1550, CL1446+26, and A781.
We analyzed the radio continuum emission and
the spectral index properties between 1.4 and 0.3 GHz of A781
in a separate paper (Govoni et al. 2011). In this paper
we present the properties of the remaining four
clusters (A800, A910, A1550, CL1446+26),
for which only data at 1.4 GHz are available
in the archive.
We note that a peripheral diffuse emission,
classified as a radio relic, was previously detected in CL1446+26
by Giovannini et al. (2009). However, on the basis of the data presented here
the classification of the diffuse emission in this cluster
is more complex than previously stated.

The details of the analyzed radio observations
(pointing, frequency, bandwidth, VLA configuration,
observing time, observing date, and VLA program) are reported in
Table\,\ref{tabl}.

To separate the diffuse radio emission from discrete sources, we produced
with standard procedures images at higher resolution obtained
at 1.4 GHz with the VLA in C configuration for the galaxy
clusters A800, A910, and A1550.
In the case of CL1446+26 we used the VLA data
in C configuration previously analyzed by Giovannini et al. (2009).

The position and the flux density\footnote{All the flux densities given in the tables and in the text have been
primary beam corrected.} of the discrete sources
embedded in the cluster diffuse emission are
given in Table\,\ref{tab2}.
These values have been subtracted to the flux density
obtained by integrating the surface brightness
of the D array data set down to the 3$\sigma$ level to estimate
the residual flux density of the diffuse emission (see Table\,\ref{radio}).
We note, however, the residual flux density associated with
the diffuse cluster emission must be interpreted with caution.
Indeed, a possible variation in the discrete sources flux density,
the slightly different frequency of the two data-sets,
and any absolute calibration error between
them could result in under or over subtraction of the calculated flux density.
In addition, owing to the short exposure time of the archival observations,
some diffuse emission could be missed.
Therefore, a deep follow-up investigation would be necessary
to ensure the recovery of the entire radio flux and to unambiguously
separate the diffuse emission from that
of the unrelated discrete sources.

To ensure that the large-scale diffuse emission
is not caused by the blending of discrete sources, for each cluster
we obtained a total intensity image at 1.4 GHz with
the VLA in D configuration, after subtraction of discrete sources.
We first produced an image of the discrete
sources by using only the longest baselines of the
D configuration data set.
The clean components of the discrete sources were then subtracted
from the original data set in the uv-plane by using the AIPS task UVSUB.
At this points, we calculated the flux density of the residual diffuse emission
by integrating the surface brightness down to the 3$\sigma$ level,
and we checked
its consistency with the values reported in Table\,\ref{radio}.
We found that the agreement is generally good (see the notes on individuals clusters).

\begin{table*}
\caption{Details of the VLA observations.}
\label{tab1}
\centering
\begin{tabular}{lcccccclcr}
\hline\hline
Cluster&  R.A.       & Decl.   &Frequency   & Bandwidth & Conf. & Obs. Time & Date  & Program  \\
       &  (J2000)    & (J2000) & (GHz)      & (MHz)     &       & (min.)   &       &          \\
\hline
A800    & 09 28 22.9 & 37 48 09.0 &  1365/1435 &  7$\times$3.125     &  D & 15    & 1995-March-15 & AM0469  \\
        & 09 25 41.9 & 37 56 50.0 &  1452/1502 &  25                 &  C &  7    & 1984-May-04   & AO0048 \\
A910    & 10 02 57.7 & 67 09 23.0 &  1365/1435 &  7$\times$3.125     &  D & 15    & 1995-March-15 & AM0469  \\
        & 10 03 18.8 & 67 11 02.4 &  1465/1385 &  50                 &  C &  3    & 1993-Jun-28   & AL0297 \\
A1550   & 12 29 19.2 & 47 37 58.0 &  1365/1435 &  7$\times$3.125     &  D & 15    & 1995-March-15 & AM0469 \\
        & 12 28 57.7 & 47 37 58.0 &  1365/1435 &  7$\times$3.125     &  C & 75    & 2003-Jan-03   & AM0702 \\
CL1446+26 & 14 49 28.7 & 26 07 54.1 &  1365/1435 &  7$\times$3.125     &  D & 20    & 1995-March-15 & AM0469     \\
\hline
\multicolumn{9}{l}{\scriptsize Col. 1: Cluster name; Col. 2, Col. 3: Observation pointing
(R.A. J2000, Decl. J2000); Col. 4: Frequency (IF1/IF2); Col. 5: Bandwidth;}\\
\multicolumn{9}{l}{\scriptsize  Col. 6: VLA configuration;
Col. 7: Observing time (minutes); Col. 8: Observing date; Col. 9: VLA program.}\\
\end{tabular}
\label{tabl}
\end{table*}

\section{New diffuse radio emissions}

For each cluster we analyze the results of the VLA observations with a
particular emphasis on the analysis of the radio properties of the
large-scale cluster diffuse emission.

There are only sparse radio, X-ray, and optical information in the
literature for the clusters analyzed in this paper, and nothing is
reported regarding their dynamical state.
In order to test the possible connection between the detection of
diffuse radio emission and the dynamical state of the observed
clusters, we compare our results with the gas and galaxy distribution
of each system.

We investigated the gaseous environment of these clusters
by checking the X-ray emission in the Rosat
All Sky Survey (RASS). The X-ray images in the 0.1$-$2.4\,{\rm keV}
band were background and exposure corrected, and
smoothed with a $\sigma$=45$''$ Gaussian kernel.

We produced density maps of the bi-dimensional galaxy distribution of each observed cluster
on the basis of a multi-scale approach, as described in Ferrari et al. (2005).
In order to avoid as much as possible projection effects, we selected galaxies
lying on the cluster red-sequence and brighter than $L^* + 3$.
For most of the systems of our sample we used the SLOAN Digital Sky Survey
(SDSS) multi-band galaxy
catalogue (Abazajian et al. 2009), applying a color-magnitude selection
based on the $(r - i)$ vs. $r$ and $(i - z)$ vs. $r$ plots (see, e.g., Goto et al. 2002).
Only the field around A910 is not covered by SDSS observations,
therefore we exploited the SuperCOSMOS catalogue (Hambly et al. 2001) and adopted
the $(B - R)$ vs. $B$ and $(B - R)$ vs. $R$ plots for color selection.

In the following, for each cluster we present the radio emission
compared with optical and X-ray images.

\begin{table}
\caption{Information on discrete radio sources.}
\begin{center}
\begin{tabular} {lcccr}
\hline\hline
Cluster & Label   &  R.A.      &  Decl.       &  S$_{\rm 1.4 GHz}$      \\
        &         &  (J2000)   &  (J2000)     &  (mJy)       \\
\hline
A800    &    A    & 09 28 36.1 & 37 47 06     &  3.5$\pm$ 0.3     \\
        &    B    & 09 28 27.5 & 37 45 12     &  2.3$\pm$ 0.2      \\
\hline
A910    &    A    & 10 03 20.8 & 67 09 32     & 25.2$\pm$1.3    \\
        &    B    & 10 03 15.2 & 67 10 20     &  4.9$\pm$0.3   \\
        &    C    & 10 03 14.1 & 67 11 09     &  1.8$\pm$0.2   \\
        &    D    & 10 03 15.3 & 67 12 00     &  0.7$\pm$0.2   \\
        &    E    & 10 03 24.7 & 67 13 22     &  2.8$\pm$0.3    \\
        &    F    & 10 03 32.2 & 67 11 20     &  2.8$\pm$0.3    \\
\hline
A1550   &    A    & 12 29 02.3 & 47 36 56     & 15.8$\pm$0.8    \\
        &    B    & 12 28 58.5 & 47 38 26     &  2.2$\pm$0.2    \\
        &    C    & 12 29 15.7 & 47 37 04     &  0.8$\pm$0.2    \\
\hline
CL1446+26 &    A  & 14 49 29.7 & 26 07 54     & 29.7$\pm$1.5    \\
          &    B  & 14 49 28.2 & 26 08 22     &  2.4$\pm$0.2   \\
          &    C  & 14 49 30.6 & 26 09 11     &  2.2$\pm$0.2   \\
          &    D  & 14 49 33.1 & 26 04 06     & 11.0$\pm$0.5   \\
          &    E  & 14 49 26.7 & 26 05 21     &  3.2$\pm$0.2   \\
          &    F  & 14 49 18.2 & 26 05 30     &  2.5$\pm$0.2   \\
\hline
\multicolumn{5}{l}{\scriptsize Col. 1: Cluster name; Col. 2: Source label;}\\
\multicolumn{5}{l}{\scriptsize Col. 3, 4: Source position (R.A. J2000, Decl. J2000);}\\
\multicolumn{5}{l}{\scriptsize Col. 5: Source flux density (from the C array data set.)}\\
\end{tabular}
\label{tab2}
\end{center}
\end{table}

\begin{figure*}[t]
\centering
\includegraphics[width=16 cm]{Fig1.ps}
\caption{{\bf Abell 800.}
  {\it Top, left}: total intensity radio contours of A800 at 1.4
  GHz with the VLA in D configuration. The image has a FWHM of
  63$^{\prime\prime}$$\times$63$^{\prime\prime}$.  The contour
  levels are drawn at $-$0.35 (green), 0.35 mJy/beam, and the rest are
  spaced by a factor of $\sqrt2$.  The sensitivity (1-$\sigma$) is 0.12
  mJy/beam.  Total intensity radio contours are overlaid on the Rosat
  PSPC X-ray image in the 0.1$-$2.4 keV band, taken from the RASS.
  The X-ray image has been convolved with a Gaussian of
  $\sigma$=45$^{\prime\prime}$.
  {\it Top, right}: total intensity radio contours and grey scale image
  at 1.4 GHz with the VLA in D configuration after subtraction of
  discrete sources.
  {\it Bottom, left}: total intensity radio contours of
  A800 overlaid on the isodensity map of likely cluster members.
  {\it Bottom, right}: zoom of total intensity radio contours in the
  center of A800
  at 1.4 GHz with the VLA in C configuration.  The image has a FWHM
  of 17$^{\prime\prime}$$\times$17$^{\prime\prime}$.  The first contour
  level is drawn at 0.3 mJy/beam, and the rest are spaced by a factor
  $\sqrt{2}$.  The sensitivity (1$\sigma$) is 0.1 mJy/beam.  The
  contours of the radio intensity are overlaid on the optical image
  taken from the SDSS (red plate).}
\label{a800}
\end{figure*}

\subsection{Abell 800}

A800 (RXCJ0928.6+3747) is a low X-ray
luminosity cluster at a redshift z=0.2223 which is part of the
Northern ROSAT All-Sky (NORAS)
Galaxy Cluster Survey by B{\"o}hringer et al. (2000).
Its X-ray luminosity in the 0.1$-$2.4\,{\rm keV}
band is $2.25\times 10^{\rm 44}$\,{\rm erg/sec}.

Low resolution radio contours at 1.4 GHz of A800 are shown in the top-left
and bottom-left panels of Fig.~\ref{a800}.
This image was obtained with the VLA in D
configuration and has been convolved to a resolution of
63$^{\prime\prime}$$\times$63$^{\prime\prime}$.
At the top, the radio contours are overlaid on the cluster X-ray
image taken from the RASS.
At the bottom, the radio contours are overlaid on the isodensity map of
likely cluster members.
The X-ray image shows an elongation in the
northwest-southeast direction. The galaxy distribution has a main clump,
surrounded by several substructures.
The diffuse radio emission, that we classify as a radio halo,
is concentrated within the main galaxy over-densities,
while a clear shift between the galaxies and gas distribution is present.

To separate the diffuse radio emission from discrete sources we
produced an image at higher resolution. In the bottom-right panel of
Fig.~\ref{a800} we present the radio contours of A800 at 1.4 GHz
taken with the VLA in C configuration. This image has been convolved
to a resolution of 17$^{\prime\prime}$$\times$17$^{\prime\prime}$.
The radio contours are overlaid on
the optical image taken from the SDSS.  Two
discrete radio sources (labelled with A and B) are embedded in the
cluster diffuse emission. Source A is located on south of the cluster
X-ray center while it is located in coincidence with the concentration
of likely cluster members. Source B is located on south-west
of the cluster X-ray center at a projected distance of about 2.5$'$ from
source A.

The total flux density is calculated from the D configuration image
after a primary beam correction by integrating the total intensity surface
brightness in the region of the diffuse emission down to the 3$\sigma$
level. The resulting total flux density is estimated to be 16.4$\pm$0.8 mJy.
By subtracting the flux density (see Table\,\ref{tab2})
of the embedded discrete sources A and B, a flux density of
10.6$\pm$0.9 mJy appears to be associated with the low brightness
diffuse emission.
This flux density value corresponds to a radio
power of $P_{\rm 1.4\,GHz} =1.52\times10^{24}$\,{\rm W~Hz}$^{\rm -1}$.

To ensure that the large-scale diffuse emission
is not caused by the blending of discrete sources, in
the top-right panel of Fig.~\ref{a800} we present
the total intensity radio contours at 1.4 GHz with the
VLA in D configuration after subtraction of discrete sources.
In perfect agreement with the value given in Table\,\ref{radio},
the flux density of the diffuse emission calculated on the image
with the point sources subtracted is 10.7$\pm$0.7 mJy.
As measured from the 3$\sigma$ radio
isophote, the overall diffuse emission shows a Largest Linear Size
(LLS) of about 6$'$ ($\simeq 1.28$ Mpc).
However, as pointed out in Murgia et al. (2009), we note that the size
of the diffuse emission calculated from the contour levels should be
considered carefully, since it depends on the sensitivity of the
radio image.

\subsection{Abell 910}

A910 (RXCJ1003.1+6709) is a galaxy cluster at a redshift z=0.2055 which
belongs to the extended sample of the Rosat Brightest Cluster sample
by Ebeling et al. (2000). Its X-ray luminosity in the 0.1$-$2.4\,{\rm keV}
band is $3.32\times 10^{\rm 44}$\,{\rm erg/sec}. This value is in agreement within the
error with that given by B{\"o}hringer et al. (2000).

The radio contours at 1.4 GHz of A910 are shown in the top-left
and bottom-left panels of Fig.~\ref{a910}.
This image was obtained with the VLA in D
configuration and has been convolved to a resolution of
63$^{\prime\prime}$$\times$63$^{\prime\prime}$.
At the top, the radio contours are overlaid on the cluster X-ray
image taken from the RASS.
At the bottom, the radio contours are overlaid on the isodensity map of
likely cluster members.
In A910 at least three X-ray peaks can be identified, possible
indications of a multiple merger.
The galaxy distribution shows several substructures too.
The diffuse radio emission has been
detected in correspondence of a bridge connecting a bright X-ray
region with a fainter clump.

In the bottom-right panel of
Fig.~\ref{a910} we present the radio contours of A910 at 1.4 GHz
taken with the VLA in C configuration. This image has been convolved
to a resolution of 18$^{\prime\prime}$$\times$18$^{\prime\prime}$.
The radio contours are overlaid on
the optical image taken from the Digital Sky Survey.
Several discrete radio sources are embedded in the
cluster diffuse emission.

Note that, in the classical framework of cluster radio source
classification (e.g. Ferrari et al. 2008, Feretti et al. 2012),
the diffuse radio source in A910 could be considered as a relic due to its morphology
and the off-set between the radio and both the intra-cluster medium
and the galaxy density peaks.
Indeed, as shown in top-right panel of Fig.~\ref{a800} where the discrete sources
have been subtracted, the diffuse emission shows an elongated structure
with a LLS of about 5$'$ ($\simeq 1.0$ Mpc) offsetted
from the main gas density and galaxy density peaks of about
135$''$ (450 kpc) and about 200$''$ (670 kpc) respectively.
Therefore, given its location and morphology we classify
this cluster diffuse emission as a relic.

The total flux density of the relic is calculated
from the D configuration image
by integrating the total intensity surface
brightness in the region of the diffuse emission down to the 3$\sigma$
level. The resulting total flux density is estimated to
be 50.3$\pm$2.5 mJy.
By subtracting the flux density of
the discrete sources A, B, C, D, E, and F (see Table\,\ref{tab2}),
a flux density of 12.1$\pm$2.9 mJy appears to be associated
with the low brightness relic emission.
This flux density value corresponds to a radio
power of $P_{\rm 1.4\,GHz} =1.45\times10^{24}$\,{\rm W~Hz}$^{\rm-1}$.
The flux density calculated on the image with the point sources
subtracted is 9.0$\pm$0.6 mJy, consistent within the errors with the
value given above.

\begin{figure*}[t]
\centering
\includegraphics[width=16 cm]{Fig2.ps}
\caption{{\bf Abell 910.}
  {\it Top, left}: total intensity radio contours of A910 at 1.4
  GHz with the VLA in D configuration.
  The image has a FWHM of
  63$^{\prime\prime}$$\times$63$^{\prime\prime}$.  The contour
  levels are drawn at $-$0.45 (green), 0.45 mJy/beam, and the rest are
  spaced by a factor of $\sqrt2$.  The sensitivity (1-$\sigma$) is
  0.15 mJy/beam.  Total intensity radio contours are overlaid on the
  Rosat PSPC X-ray image in the 0.1$-$2.4 keV band, taken from the RASS.
  The X-ray image has been convolved with a Gaussian
  of $\sigma$=45$^{\prime\prime}$.
  {\it Top, right}: total intensity radio contours and grey scale image
  at 1.4 GHz with the VLA in D configuration after subtraction of
  discrete sources.
  {\it Bottom, left}: total intensity radio contours
  of A910 overlaid on the isodensity map of likely cluster members.
  {\it Bottom, right}: zoom of total intensity radio contours in
  the center of A910 at 1.4 GHz with the VLA in C configuration.
  The image has a FWHM of
  18$^{\prime\prime}$$\times$18$^{\prime\prime}$.  The first contour
  level is drawn at 0.3 mJy/beam, and the rest are spaced by a factor
  $\sqrt{2}$.  The sensitivity (1$\sigma$) is 0.1 mJy/beam.  The
  contours of the radio intensity are overlaid on the optical image
  taken from the Digital Sky Survey (red plate).}
\label{a910}
\end{figure*}

\subsection{Abell 1550}

\begin{figure*}[t]
\centering
\includegraphics[width=16 cm]{Fig3.ps}
\caption{{\bf Abell 1550.}
  {\it Top, left}: total intensity radio contours of A1550 at 1.4
  GHz with the VLA in D configuration.
  The image has a FWHM of
  53$^{\prime\prime}$$\times$53$^{\prime\prime}$.  The contour
  levels are drawn at $-$0.2, 0.2 mJy/beam, and the rest are spaced by
  a factor of $\sqrt2$.  The sensitivity (1-$\sigma$) is 0.07
  mJy/beam.  Total intensity radio contours are overlaid on the Rosat
  PSPC X-ray image in the 0.1$-$2.4 keV band, taken from the RASS.
  The X-ray image has been convolved with a Gaussian of
  $\sigma$=45$^{\prime\prime}$.
  {\it Top, right}: total intensity radio contours and grey scale image
  at 1.4 GHz with the VLA in D configuration after subtraction of
  discrete sources.
  {\it Bottom, left}: total intensity radio contours of
  A1550 overlaid on the isodensity map of likely cluster members.
  {\it Bottom, right} : zoom of total intensity radio contours in the center of
  A1550 at 1.4 GHz with the VLA in C configuration. The image has
  a FWHM of 18$^{\prime\prime}$$\times$18$^{\prime\prime}$.
  The first contour
  level is drawn at 0.15 mJy/beam, and the rest are spaced by a factor
  $\sqrt{2}$.  The sensitivity (1$\sigma$) is 0.04 mJy/beam.  The
  contours of the radio intensity are overlaid on the optical image
  taken from the SDSS (red plate).}
\label{a1550}
\end{figure*}

A1550 (RXCJ1229.0+4737) is a galaxy cluster at a redshift z=0.2540 with
an X-ray luminosity in the 0.1$-$2.4\,{\rm keV} band
of $3.51\times 10^{\rm 44}$\,{\rm erg/sec} (B{\"o}hringer
et al. 2000).

Low resolution radio contours at 1.4 GHz of A1550 are shown in the top-left
and bottom-left panels of Fig.~\ref{a1550}.
This image was obtained with the VLA in D
configuration and has been convolved to a resolution of
53$^{\prime\prime}$$\times$53$^{\prime\prime}$.
At the top, the radio contours are overlaid on the cluster X-ray
image taken from the RASS.
At the bottom, the radio contours are overlaid on the isodensity map of
likely cluster members.  The external X-ray emission is elongated in
the east-west direction, while the inner region is extended along
north-south.  The optical galaxy distribution is bimodal but with an
elongation that however does not follow the cluster X-ray emission.
Surrounding the X-ray and the optical center we detected a diffuse
low-surface brightness radio emission, which we classified as a radio halo.
As in the case of A800, the diffuse radio emission seems to be more spatially
correlated to the galaxy than to the gas over-densities.

To separate the diffuse radio emission from discrete sources we
produced an image at higher resolution.
In the bottom-right panel of
Fig.~\ref{a1550} we present the radio contours of A1550 at 1.4 GHz
taken with the VLA in C configuration. This image has
been convolved to a resolution of
18$^{\prime\prime}$$\times$18$^{\prime\prime}$.
The radio contours are overlaid on the optical
image taken from the SDSS. Three discrete radio
sources (labelled as A, B, and C) are clearly embedded in the radio halo. Spots of very faint
emission are also detected in the C configuration data set, but given
that they don't coincide with any discrete source detected in the
FIRST survey (Becker et al. 1995) we consider these emissions associated with the peaks of
the radio halo detected in the D configuration data set.

A total flux density of $\simeq$26.5$\pm$1.3 mJy is calculated
from the D configuration image, by integrating the total
intensity surface brightness in the region of the diffuse emission
down to the 3$\sigma$ level.
By subtracting the flux density
of the discrete sources A, B, and C (see Table ~\ref{tab2}),
a flux density of 7.7$\pm$1.6 mJy appears to be associated
with the low brightness diffuse emission.
This flux density value corresponds to a radio power
of $P_{\rm 1.4\,GHz} = 1.49\times10^{24}$\,{\rm W~Hz}$^{\rm -1}$.

In the top-right panel of Fig.~\ref{a1550} we present
the total intensity radio contours at 1.4 GHz with the
VLA in D configuration after subtraction of discrete sources.
The flux density calculated on the image with the point sources
subtracted is 8.2$\pm$0.5 mJy, while as measured from the 3$\sigma$ radio
isophote, the radio halo shows a LLS of about 6$'$ ($\simeq 1.41$ Mpc).

\subsection{CL 1446+26 - CL 1447+26 - ZwCl 1447.2+2619}

\begin{figure*}
\centering
\includegraphics[width=16 cm]{Fig4.ps}
\caption{{\bf CL 1446+26.}
  {\it Top, left}: total intensity radio contours of CL 1446+26 at
  1.4 GHz with the VLA in D configuration.  The image has a FWHM of
  52$^{\prime\prime}$$\times$52$^{\prime\prime}$.  The contour levels are drawn
  at $-$0.25, 0.25 mJy/beam, and the rest are spaced
  by a factor of $\sqrt2$.  The sensitivity (1-$\sigma$) is 0.08
  mJy/beam.  Total intensity radio contours are overlaid on the Rosat
  PSPC X-ray image in the 0.1$-$2.4 keV band, taken from the RASS.
  The X-ray image has been convolved with a Gaussian of
  $\sigma$=45$^{\prime\prime}$.
  {\it Top, right}: total intensity radio contours and grey scale image
  at 1.4 GHz with the VLA in D configuration after subtraction of
  discrete sources.
 {\it Bottom, left}: total intensity radio contours of
  CL 1446+26 overlaid on the isodensity map of likely cluster members.
 {\it Bottom, right}: zoom of total intensity radio contours in the center
  of CL1446+26 at 1.4 GHz with the VLA in C configuration, image (VLA
  program AO149) taken from Giovannini et al. (2009).  The image has
  a FWHM of 15$''\times 15''$.  The first contour level is
  drawn at 0.15 mJy/beam, and the rest are spaced by a factor
  $\sqrt{2}$.  The sensitivity (1$\sigma$) is 0.05 mJy/beam.  The
  contours of the radio intensity are overlaid on the optical image
  taken from the SDSS (red plate).}
\label{cl1446}
\end{figure*}

CL 1446+26 is a rich galaxy cluster located at z=0.370 with an X-ray
luminosity\footnote{The bolometric X-ray luminosity has been converted to the
0.1$-$2.4 keV band using Table 5
of B{\"o}hringer et al. (2004), assuming an
intra-cluster temperature of T$\simeq$ 5 keV.}
in the 0.1$-$2.4 keV band
of 3.41$\times$10$^{44}$erg/s (Wu et al. 1999).

At radio wavelengths, it has been studied by Owen et al. (1999) and Giovannini et al. (2009).
The presence of a northern radio relic plus several radio galaxies
were pointed out by Giovannini et al. (2009)
(see also radio contours in the bottom-right panel of Fig.~\ref{cl1446}).

Low resolution radio contours at 1.4 GHz of CL 1446+26 are shown in the top-left
and bottom-left panels of Fig.~\ref{cl1446}. This image was obtained with the VLA in D
configuration and has been convolved to a resolution of 52$^{\prime\prime}$$\times$52$^{\prime\prime}$.
At the top, the radio contours are overlaid on the cluster X-ray
image taken from the RASS.
At the bottom, the radio contours are overlaid on the isodensity map of
likely cluster members.
The cluster is characterized by a disturbed morphology both in X-rays
and in optical, with two main clumps, one to the North and one to the South.
Again, the optical and X-ray peaks are not coincident.
As in A800 and A1550, the overall diffuse radio emission seems to follow
more closely the galaxy than the intra-cluster density distribution.
Our lower resolution radio map shows in coincidence with the Northern
clump a radio emission more extended than that
obtained by Giovannini et al. (2009). In particular, this new image
seems to indicate that
the diffuse emission is indeed a central radio halo rather than a peripheral
relic source, as previously stated.
In addition, a new elongated diffuse radio emission
in correspondence of the Southern sub-cluster has been detected.

In the top-right panel of Fig.~\ref{cl1446} we present
the total intensity radio contours at 1.4 GHz with the
VLA in D configuration after subtraction of discrete sources.
Given their location and morphology, in the following we define
the norther emission as a radio halo and the elongated southern emission
as a relic. The radio halo appears to be connected through a low brightness
diffuse bridge to the relic.
However, a deeper follow-up is necessary to clarify the nature of this system.

In the Northern emission a total flux density of 42.0$\pm$2.1 mJy
is calculated from the D configuration image, by integrating the total intensity
surface brightness in the region of the diffuse emission
down to the 3$\sigma$ level.
By subtracting the flux density
of the discrete sources A, B and C (see Table\,\ref{tab2}) a flux density of
7.7$\pm$2.6 mJy appears to be associated to the radio halo.
This flux density value corresponds to a radio
power of $P_{\rm 1.4\,GHz} = 3.57\times10^{24}$\,{\rm W~Hz}$^{\rm -1}$.
While in the Southern emission by subtracting the flux density
of the discrete sources D, E, and F from the total flux density
(22.0$\pm$1.1 mJy), a flux density of
5.3$\pm$1.2 mJy appears to be associated to the relic.
This flux density value corresponds to a radio
power of $P_{\rm 1.4\,GHz} = 2.46\times10^{24}$\,{\rm W~Hz}$^{\rm -1}$.

The LLS is about 4$'$ ($\simeq 1.22$ Mpc) and 5$'$ ($\simeq 1.53$ Mpc)
for the radio halo and relic, respectively.
The flux densities calculated on the image with the point
sources subtracted are consistent within the errors with the
values given in Table\,\ref{radio}, being
 5.3$\pm$0.3 mJy in the halo and 5.5$\pm$0.3 mJy in the relic.

\section{Statistical properties}

\begin{table*}
\caption{Radio information of clusters containing diffuse emission.}
\label{radio}
\centering
\begin{tabular}{lclrrrrrr}
\hline\hline
Name    & z      & kpc/$''$ & Type  &$S_{1.4 GHz}$ & log $P_{1.4 GHz}$ & LLS & Radio$-$X-ray   & Radio$-$Optical \\
        &        &          &       &(mJy)       & (W/Hz)          & (Mpc) & Offset ($''$) & Offset ($''$)  \\
\hline
A800    & 0.2223 &  3.55    & H    &  10.6    &  $1.52\times10^{24}$ & 1.28  &  90 & 5 \\
A910    & 0.2055 &  3.34    & R   &  12.1    &  $1.45\times10^{24}$  & 1.00  & 135 & 200\\
A1550   & 0.2540 &  3.92    & H    &   7.7    &  $1.49\times10^{24}$ & 1.41  & 105 & 5   \\
CL1446+26N & 0.370 & 5.09 & H      &  7.7     &  $3.57\times10^{24}$ & 1.22  & 85  & 40  \\
CL1446+26S & 0.370 & 5.09 & R      &  5.3     &  $2.46\times10^{24}$ & 1.53  & 205 & 250 \\
\hline
\multicolumn{9}{l}{\scriptsize Col. 1: Cluster Name; Col. 2: Redshift; Col. 3: Angular to linear conversion; Col. 4: Type of diffuse radio emission:}\\
\multicolumn{9}{l}{\scriptsize H=halo, R=relic; Col. 5: Diffuse radio flux density at 1.4 GHz; Col. 6: Diffuse radio power at 1.4 GHz;}\\
\multicolumn{9}{l}{\scriptsize Col. 7: Radio Largest Linear Size (LLS); Col. 8: Radio-X-ray offset;}\\
\multicolumn{9}{l}{\scriptsize Col. 9: Radio-Optical offset.}\\
\end{tabular}
\end{table*}

\begin{figure*}
\centering
\includegraphics[width=18 cm]{Fig5.ps}
\caption{
{\it Left}: Monochromatic radio power of radio halos at 1.4 GHz versus
the cluster X-ray luminosity in the 0.1$-$2.4 keV band.
{\it Right}: Monochromatic radio power of halos at 1.4 GHz versus
their largest linear size (LLS) measured at the same frequency.
The data are taken from the recent
compilation by Feretti et al. (2012) and
references therein. Blue dots are classical radio halos while
empty blue dots are outlier clusters. The new halos detected in this work
are indicated with a red triangles.
}
\label{relation}
\end{figure*}

\begin{table}
\caption{Radio-X-ray offsets for a sample of halos.}
\centering
\begin{tabular}{lclrrl}
\hline\hline
Name      & z        & kpc/$''$ & Radio-X-ray    &  LLS &  Ref  \\
          &          &          & Offset ($''$) &  ($'$) &       \\
\hline
A209      &  0.2060  &   3.34   &   20     &  7.0 &  1 \\
A399      &   0.0718 &   1.35   &   135    &  7.0 &  2 \\
A401      &   0.0737 &   1.38   &    70    &  6.3 &  3 \\
A520      &   0.1990 &   3.25   &    20    &  5.7 &  4 \\
A523      &  0.1036  &  1.88    &      185 & 12.0 &  5 \\
A545      &  0.1540  &  2.64    &      50  &  5.6 &  3 \\
A665      &  0.1819  &   3.03   &      50  &   10 &  6 \\
A754      &  0.0542  &   1.04   &     142  & 15.8 & 3  \\
A773      &  0.2170  &   3.48   &      30  &  6.0 & 4  \\
A781      &  0.3004  &   4.42   &      55  &  6.0 & 7  \\
A800      &  0.2223  &   3.55   &      90  &  6.0 & *  \\
A851      &  0.4069  &   5.40   &      75  &  3.3 & 1  \\
A1213     &  0.0469  &   0.91   &     145  &  4.0 & 1  \\
A1351     &  0.3224  &   4.65   &     100  &  3.0 & 1  \\
A1550     &  0.2540  &   3.92   &     105  &  6.0 & *  \\
A1689     &  0.1832  &   3.05   &      15  &  4.0 & 8  \\
A1758a    &  0.2790  &   4.20   &      35  &  6.0 & 1  \\
A1914     &  0.1712  &   2.88   &      45  &  7.4 & 3  \\
A1995     &  0.3186  &   4.61   &      45  &  3.0 & 1  \\
A2034     &  0.1130  &   2.03   &      65  &  5.0 & 1  \\
A2163     &  0.2030  &   3.31   &      30  & 11.5 & 9  \\
A2218     &  0.1756  &   2.94   &      35  &  2.2  & 10  \\
A2219     &  0.2256  &   3.59   &      20  &  8.0 & 3  \\
A2254     &  0.1780  &   2.98   &      45  &  5.0 & 4  \\
A2255     &  0.0806  &   1.50   &      65  & 10.0 & 11  \\
A2294     &  0.1780  &   2.98   &      25  &  3.0 & 1  \\
A2319     &  0.0557  &   1.07   &      35  & 15.9 & 12 \\
A2744     &  0.3080  &   4.50   &      25  &  7.0 & 4  \\
RXCJ1314  &  0.2439  &   3.81   &     100  &  7.0 & 13 \\
CL1446+26 &  0.3700  &   5.09   &      85  &  4.0 & * \\
\hline
\multicolumn{6}{l}{\scriptsize Col. 1: Cluster Name; Col. 2: Redshift;}\\
\multicolumn{6}{l}{\scriptsize Col. 3: Angular to linear conversion; Col. 4: Radio - X-ray offset;}\\
\multicolumn{6}{l}{\scriptsize Col. 5: Radio largest linear size (LLS);}\\
\multicolumn{6}{l}{\scriptsize Col. 6: Reference to the radio halos at 1.4 GHz:}\\
\multicolumn{6}{l}{\scriptsize *=This work, 1=Giovannini et al. (2009),}\\
\multicolumn{6}{l}{\scriptsize 2=Murgia et al. (2010), 3=Bacchi et al. (2003),}\\
\multicolumn{6}{l}{\scriptsize 4=Govoni et al. (2001a), 5=Giovannini et al. (2011),}\\
\multicolumn{6}{l}{\scriptsize 6=Vacca et al. (2010), 7=Govoni et al. (2011),}\\
\multicolumn{6}{l}{\scriptsize 8=Vacca et al. (2011), 9=Feretti et al. (2001),}\\
\multicolumn{6}{l}{\scriptsize 10= Giovannini \& Feretti (2000), 11=Govoni et al. (2005),}\\
\multicolumn{6}{l}{\scriptsize 12=Feretti et al. (1997), 13= Feretti et al. (2005).}\\
\end{tabular}
\label{taboffset}
\end{table}

The basic properties of the new cluster
diffuse radio sources (halos and relics) are reported in Table\,\ref{radio}.

The new relics presented here show properties similar to other published
relic sources (see e.g. the statistical characteristics of radio
relics reviewed by Feretti et al. 2012), and are not further discussed
in the following.

All the clusters presented in this paper and hosting new radio halos
(A800, A1550, CL1446+26) show a clear shift between the galaxies
and gas distribution, as well as complex optical and X-ray morphologies.
Joint numerical and observational analysis indicate that these
are typical features of dynamically disturbed clusters,
with a recent or ongoing major merging event (e.g. Maurogordato et al. 2011).
Deeper X-ray observations and optical spectroscopic follow-up could help in
re-constructing the collision scenario of these perturbed systems.

In Fig.~\ref{relation} we plot the radio power calculated at
1.4 GHz ($P_{\rm 1.4 GHz}$) versus the cluster X-ray
luminosity in the 0.1$-$2.4 keV band ($L_X$) and versus the
Largest Linear Size ($LLS$).
The blue dots refer to radio halos taken from the literature.
Among them the empty dots are outliers, i.e.
on the left panel they refer to the few known
radio halos A523 (Giovannini et al. 2011), A1213 (Giovannini et al. 2009), and
0217+70 (Brown et al. 2011) that are over-luminous in radio with respect
to the cluster X-ray luminosity,
while on the right panel they refer to the few known
radio halos A1213, A1351 (Giovannini et al. 2009),
MACS J0717.5+3745 (Bonafede et
al. 2009) which shows a size smaller than
expected from their radio power.
The new radio halos are indicated with red triangles.
They agree with both the $P_{\rm 1.4 GHz}-L_X$ and $P_{\rm 1.4 GHz}-LLS$ relations
known for the other halos in clusters.

The observed correlation between the radio power of halos and their
size has been firstly found for giant radio halos
by Cassano et al. (2007). In a later analysis Murgia et al. (2009) by
fitting the azimuthally averaged brightness profile with an exponential
found that the radio halo emissivity is remarkably similar from one
halo to the other, despite of quite different length-scales,
in agreement with the $P_{\rm 1.4 GHz}-LLS$ correlation where larger halos show a
higher radio power.

A close similarity between the radio and the X-ray structures
has been found in a number of clusters hosting a radio
halo (see e.g. Govoni et al. 2001b).
This similarity is generally valid for giant and regular halos.
However, more irregular and asymmetric halos have
been found in the literature. In these halos, the radio emission
may show significant displacement from the X-ray emission. Interestingly,
the new halos studied in this paper show a clear offset between
the radio and the X-ray peak (see Table\,\ref {radio}). In addition,
the similarity between the radio and the X-ray morphology
is not clearly present and the radio
emission seems to follow more closely the galaxy distribution rather than
the intra-cluster gas density distribution.

To investigate from a statistical point of view the offset between
the radio halo and the X-ray peak of the cluster emission,
we followed the analysis by Feretti et al. (2011).

In the top panel of Fig.~\ref{offset1} we present the distribution of
the offset between the peak position of radio halo
and the X-ray gas distribution for the new radio halos presented
in this paper together with a sample of cluster at a redshift
$z$$\leq$0.4 previously imaged at 1.4 GHz by our group
(following a uniform data reduction strategy), for which good radio and
X-ray data are available (see Table\,\ref{taboffset} for details).
The mean radio$-$X-ray offset is 182 kpc.
In the bottom panel of Fig.~\ref{offset1} the same distribution
refers to large halos with $LLS > 1$ Mpc (in blue) and small halos
with $LLS < 1$ Mpc (in red), respectively.
Both large and small radio halos can be significantly shifted,
up to hundreds kpc, with respect to the X-ray center.

To highlight radio halos with the most pronounced asymmetric distribution,
the radio$-$X-ray offset normalized to the halo size ($LLS$)
is plotted in Fig. \ref{offset2}.
The mean value of the radio$-$X-ray offset normalized to the
$LLS$ is 0.19.
In the bottom panel of Fig.~\ref{offset2} the same distribution
is considered for large (in blue) and small (in red) halos separately.
The mean value calculated for large halos is 0.15, while the
mean value calculated for small halos is 0.24.
Thus we found that halos can be quite asymmetric with respect
to the X-ray gas distribution, and this becomes more
relevant when halos of smaller size are considered.

\section{Discussion}

Although the properties of the newly discovered, moderate-mass,
radio halos are consistent with the $P_{\rm 1.4 GHz}-L_X$ and $P_{\rm 1.4 GHz}-LLS$
correlations known in the literature,
their radio morphology does not matches exactly the thermal X-ray emission from
the intra-cluster medium.

By analyzing the radio X-ray offset statistically we found
a tendency for smaller halos to show larger displacements.
This result may have important consequences for the origin of radio halos,
suggesting that the morphological properties of radio halos,
may be linked to the cluster dynamical state.

Major cluster mergers can supply energy to the radio emitting
particles, as well as amplify magnetic fields
(e.g. Roettiger et al. 1999, Dolag et al. 1999,
Ricker \& Sarazin 2001, Dolag et al. 2005, Ryu et al. 2008, Vazza et al. 2009,
 Xu et al. 2010, Vazza et al. 2011, Xu et al. 2011).
Since magnetic fields have been found to be ubiquitous in galaxy clusters
(e.g. Kronberg 1994, Carilli \& Taylor 2002, Clarke et al. 2001,
Govoni \& Feretti 2004, Bonafede et al. 2011),
the crucial ingredient for the existence of diffuse radio halos
is the presence of synchrotron electrons with GeV energies.
Continuous acceleration of thermal electrons in galaxy clusters actually
cannot work because the energy gained by the particles
is dissipated to the whole plasma on a timescale much shorter than that of the acceleration
process itself (see e.g. Petrosian 2001, Wolfe \& Melia 2006).
On the other hand, it is believed that
primary relativistic electrons are present in the cluster volume because they were
injected by AGN activity or by star-forming galaxies during the cluster dynamical history.
This population of electrons suffers strong radiation losses mainly
because of synchrotron and inverse Compton emission, thus direct or
stochastic re-acceleration caused by merging shocks and turbulence
is needed to maintain their energy to the level necessary to produce the observed synchrotron
radio emission in the relatively weak intra-cluster magnetic fields
(e.g. Jaffe 1977, Schlickeiser et al. 1987, Petrosian 2001,
Brunetti et al. 2001, Brunetti et al. 2011).
Alternatively, electrons can be continuously injected in the intra-cluster medium as secondary particles by
inelastic nuclear collisions between relativistic protons and nuclei in the thermal ambient intra-cluster medium (e.g. Dennison 1980, Kushnir et al. 2009,
Keshet \& Loeb 2010, En{\ss}lin et al. 2011).
The protons diffuse on large scale because their energy losses are negligible.
They can continuously produce in situ electrons, distributed through the cluster volume.
The possibility that the high
energy electrons, responsible for the synchrotron emission, arise from the decay of
secondary products of the neutralino annihilation in the dark matter halos of galaxy
clusters has also been suggested (Colafrancesco \& Mele 2001).

A detailed discussion of the origin of the non-thermal emission in
galaxy clusters
is beyond the aims of this work, which are in fact mainly observational.
Nevertheless, we can explore, from a qualitative perspective, the reason
for the larger distortions and offsets observed among the smallest halos.
Cluster merger events are expected to release a significant amount of energy
in the intra-cluster medium. This energy is injected on large spatial scales,
and then turbulent cascades may be generated. It is generally believed that
a fraction of this energy is finally delivered in relativistic
particles acceleration and magnetic field amplification.
Thus, in this framework,
the observed offsets could be related to the age of the merger.

In particular, a possible explanation for this behavior can be attributed to
the presence of intra-cluster magnetic field spatial variations on large scales
and/or a non-uniform distribution of relativistic electrons.
Indeed, on the basis of modeling, Vacca et al. (2010) found that if
the intra-cluster magnetic field and relativistic electrons fluctuations
are close to the observing beam scale, the halo results smooth and rounded.
Increasing the magnetic field and the relativistic electrons correlation length
results in a much distorted radio halo morphology and in a significant
offset of the radio halo peak from the cluster center.
Therefore, smaller halos may have a more distorted morphology
because they may be young systems in which the energy is still
on large scales, leading to a non-uniform distribution of relativistic electrons
and/or a magnetic field correlation length
larger than in the more extended, and dynamical older, radio halos.

\begin{figure}
\centering
\includegraphics[width=9 cm]{Fig6.ps}
\caption{
{\it Top}: Offset between the radio and X-ray peaks for a sample
of radio halos (see Table\,\ref{taboffset}).
{\it Bottom}: Offset between the radio and X-ray peaks.
The blue dashed area refers to giant halos ($LLS > 1$ Mpc) and
the red  dashed area refers to small halos ($LLS < 1$ Mpc), respectively.
}
\label{offset1}
\end{figure}

\begin{figure}
\centering
\includegraphics[width=9 cm]{Fig7.ps}
\caption{
{\it Top}: Fractional offset (Offset/$LLS$) between the radio
and X-ray peaks for a sample of radio halos (see Table\,\ref{taboffset}).
{\it Bottom}: Fractional offset (Offset/$LLS$) between the
radio and X-ray peaks.
The blue dashed area refers to giant halos ($LLS > 1$ Mpc) and
the red  dashed area refers to small halos ($LLS < 1$ Mpc), respectively.
}
\label{offset2}
\end{figure}

\section{Conclusion}

In order to investigate the presence of extended diffuse radio sources
in galaxy clusters, we studied the radio continuum emission
in a sample of galaxy clusters by analyzing archival VLA observations
at 1.4 GHz in D configuration.
These observations have been complemented by X-ray, optical, and higher
resolution radio data. The higher resolution radio images ensure
that the detected diffuse emission is real and
not due to a blend of discrete sources.
The point source subtraction is rather relevant for the
detection of diffuse radio emission and the associated flux density.
Therefore, this procedure has been performed carefully by using
different methods.

These data permitted to increase the number of known large
scale diffuse radio emissions.
Indeed, we discovered a new radio halo in the central region of A800 and A1550.
We detected a radio relic in the periphery of A910, and finally,
we revealed both a halo and a relic in CL1446+26.

All the clusters presented in this paper and hosting new radio halos
show a clear shift between the galaxies
and gas distribution as well as complex optical and X-ray morphologies,
typical features of dynamically disturbed clusters.
The radio X-ray similarity is not clearly present and the radio
emission seems to follow more closely the galaxy distribution rather
than the intra-cluster gas
distribution. Interestingly, there is a clear offset between
the radio and the X-ray peak.

By analyzing this offset statistically
we found that radio halos can be quite asymmetric with respect
to the X-ray gas distribution, and this becomes more
relevant when halos of smaller size are considered.
This behavior could be linked to the cluster dynamical state.
Indeed, Vacca et al. (2010) found that a much distorted radio halo
morphology and a significant offset of the radio halo peak from
the cluster center is expected increasing the magnetic field
and the relativistic electrons correlation length in galaxy clusters.
Therefore, smaller halos may have a more distorted morphology
because they may be young systems in which the merger energy is still
on large scales, leading to a non-uniform distribution of relativistic electrons
and/or a magnetic field correlation length
larger than in the more extended, and dynamical older, radio halos.
Thus, in this framework, the observed offsets could be related
to the age of the merger.

\begin{acknowledgements}
  We thank the referee for the suggestions that
  improved the presentation of the paper.
  This research was partially supported by PRIN-INAF2009.  We
  acknowledge financial contribution from the agreement ASI-INAF
  I/009/10/0. CF and CB acknowledge financial support by the ``{\it Agence Nationale de la Recherche}''
  through grant ANR-09-JCJC-0001-01. The National Radio Astronomy Observatory (NRAO) is a
  facility of the National Science Foundation, operated under
  cooperative agreement by Associated Universities, Inc.  Funding for
  the SDSS and SDSS-II has been provided by the Alfred P. Sloan
  Foundation, the Participating Institutions, the National Science
  Foundation, the U.S. Department of Energy, the National Aeronautics
  and Space Administration, the Japanese Monbukagakusho, the Max
  Planck Society, and the Higher Education Funding Council for
  England. The SDSS Web Site is http://www.sdss.org/.  This research
  made use of Montage, funded by the National Aeronautics and Space
  Administration's Earth Science Technology Office, Computational
  Technnologies Project, under Cooperative Agreement Number NCC5-626
  between NASA and the California Institute of Technology. The code is
  maintained by the NASA/IPAC Infrared Science Archive.
\end{acknowledgements}