Microstructures and Mechanical Properties of Dense Particle Gels: Microstructural Characterization

Abstract: The macroscopic mechanical properties of densely packed coagulated colloidal particle gels strongly depend on the local arrangement of the powder particles on length scales of a few particle diameters. Heterogeneous microstructures exhibit up to one order of magnitude higher elastic properties and yield strengths than their homogeneous counterparts. The microstructures of these gels are analyzed by the straight path method quantifying quasi-linear particle arrangements of particles. They show similar characteristics than force chains bearing the mechanical load in granular material. Applying this concept to gels revealed that heterogeneous colloidal microstructures show a significantly higher straight paths density and exhibit longer straight paths than their homogeneous counterparts.

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

Mechanical properties of coagulated colloids are important in many
technical areas. Sediments~\cite{Yun_2007}, ceramic pastes and
suspensions~\cite{Feng_2000}, pharmaceutical formulations such as
cr\`{e}mes and emulsions~\cite{Adeyeye_2002} and some
food~\cite{Marti_2005} are examples, for which it is desirable to
control the mechanical properties.

Recently, an internal gelation method (DCC = direct coagulation
casting)~\cite{Gauckler_1999,Tervoort_2004} was developed to
process electrostatically stabilized colloidal suspensions to
coagulated particle gels. Thereby it was found that this method
permits to control the gels' microstructures for volume fractions
ranging from 0.2 to 0.6. The method allows for an \it in situ\rm,
i.e., undisturbed destabilization of the colloidal suspension by
either changing the pH of the solution ($\Delta$pH-method)
resulting in a ``homogeneous'' microstructure or by increasing the
ionic strength ($\Delta$I-method) leading to a ``heterogeneous''
microstructure~\cite{Wyss_2004}.

\begin{figure}
\includegraphics[width=\linewidth]{Fig1.eps}
\caption{Microstructure of coagulated silica suspensions at a
volume fraction of 0.4 formed by the $\Delta$pH- (left) and the
$\Delta$I-method (right) using the DCC process (particle diameter
0.525~$\mu$m). Micrographs obtained by cryogenic scanning electron
microscopy~\cite{Wyss_2004}.} \label{fig:load_Fig1}
\end{figure}

Figure~\ref{fig:load_Fig1} shows cryogenic scanning electron
microscopy (cryo-SEM) pictures of both microstructures. The
particles in the $\Delta$pH-coagulated gel present a highly
homogeneous microstructure whereas in the $\Delta$I-system
inhomogeneities on the length scale of a few particle diameters
are observed. This qualitative observation was quantified by the
three-dimensional pair-correlation function calculated from stereo
cryo-SEM images~\cite{Wyss_2004} and by \it in situ \rm performed
diffusing wave spectroscopy (DWS) experiments during
destabilization~\cite{Wyss_2001}. The structural differences
correspond to differences in the heterogeneity: the first peak of
the pair-correlation function at $r/d = 1$ ($r$ being the distance
between the particles and $d$ the particle diameter) is higher for
the $\Delta$I-system indicating locally denser regions with higher
average coordination number. Also, local maxima are present in the
pair-correlation function for the $\Delta$I-system at $r/d \approx
1.4$ and at $r/d \approx 1.6$ corresponding to characteristic
peaks in hexagonally packed particle arrangements. During the
$\Delta$I-destabilization, particle rearrangements lead to larger
pores and thus to less homogeneous
microstructures~\cite{Wyss_2001}, whereas during the
$\Delta$pH-destabilization, the initially stabilized, liquid-like
microstructure, is ``frozen'', which results in a more homogeneous
microstructure.

Rheological and uniaxial compression experiments on coagulated
colloidal structures obtained by
DCC~\cite{Balzer_1999,Wyss_1_2005} showed that colloids with
heterogeneous microstructures have significantly higher elastic
moduli and yield strengths than their homogeneous counterparts.
Measuring the dynamics of the colloids during destabilization by
DWS confirmed these findings quantitatively using a model for the
storage modulus proposed by Krall and
Weitz~\cite{Wyss_2001,Krall_1998}.

Alkali-swellable polymer (ASP) particles were used in order to
introduce heterogeneities during the pH-destabilization. Small
amounts of 80~nm ASP particles were admixed to the powders under
acidic conditions. These particles swell upon changing pH during
the internal gelling reaction of the DCC process and enfold to
0.7~$\mu$m size producing less homogeneous microstructures. Those
samples with swollen polymer particles showed much higher
mechanical properties in comparison to samples without swellable
polymers and hence more homogenous
microstructures~\cite{Hesselbarth_2000}. In particular, they
present the same high mechanical properties as samples with
heterogeneous microstructures produced by the
$\Delta$I-method~\cite{Hesselbarth_2000}.

Examples of the rheologically measured elastic properties of
alumina particle suspensions (average particle diameter $d_0 =
0.4$~$\mu$m) for different volume fractions destabilized by the
$\Delta$pH- and the $\Delta$I-method, respectively, are shown in
Fig.~\ref{fig:load_Fig2}. Almost one order of magnitude higher
elastic plateau storage moduli are measured for heterogeneous
microstructures than for those with homogenous
microstructures~\cite{Wyss_1_2005}.

\begin{figure}
\includegraphics[width=\linewidth]{Fig2.eps}
\caption{Elastic plateau storage modulus $G'_p$ of alumina
particle suspensions (average particle diameter $d_0 =
0.4$~$\mu$m) in dependence of the volume fraction formed by the
$\Delta$pH- and the $\Delta$I-method of the DCC
process~\cite{Wyss_1_2005}.} \label{fig:load_Fig2}
\end{figure}

In summary, strong evidence is given that the differences in
macroscopic mechanical properties of coagulated particle
suspensions are controlled by the differences in microstructure.
An open question is now how these microstructural differences on
particle length scales can have such a dramatic influence on the
mechanical properties. Therefore, a concise quantitative analysis
of microstructures of colloidal particle systems is needed.

In preceding works, various characterization methods, such as the
radial pair-correlation function~\cite{Huetter_2000}, the bond
angle distribution function~\cite{Huetter_2000}, the triangle
distribution function~\cite{Huetter_2000} and the Minkowski
functionals in conjunction with the parallel-body
technique~\cite{Huetter_2003}, were applied to sets of
microstructures generated by Brownian dynamics (BD)
simulations~\cite{Huetter_1999}. These simulations were used to
study the coagulation dynamics and the evolving microstructures in
dense colloidal suspensions and the resulting microstructures
agree well with experiments~\cite{Wyss_2002}.

While the pair-correlation function permits to quantify the amount
of structural rearrangement during the
coagulation~\cite{Huetter_2000}, the analysis using the Minkowski
functionals in conjunction with the parallel-body technique
supplies additional information on the structure's morphology
resolving microstructural differences on a length scale limited by
the largest pore size~\cite{Huetter_2003}. The bond angle
distribution function and the triangle distribution functions are
useful means to examine the local building blocks of the particle
network~\cite{Huetter_2000}. Particular features, as, for example,
peaks in the respective distribution function have successfully
been correlated to the structure's heterogeneity and porosity. The
same conclusion is valid for the Minkowski
functionals~\cite{Huetter_2003}.

All these four methods are good means to compare structures in
terms of their heterogeneity. However, these structural
descriptions do not unambiguously help to understand why more
heterogeneous colloidal structures possess stronger mechanical
properties, as they do not adequately capture the microstructural
characteristics that are responsible for the mechanical properties
of these particle systems.

It is well known from granular matter physics that particulate
systems under mechanical stress carry load via chains of
contacting particles, termed as force chains, as observed in
experiments on granular materials~\cite{Drescher_1972} and in
simulations~\cite{Cundall_1979}. Granular materials, as, for
example, sand piles, that contain large amounts of particles
arranged in chains of contacting particles are expected to possess
higher mechanical properties than those in which the particles
have to rearrange upon applied external load in order to form such
chains. Hence it is plausible that these chains of contacting
particles control the mechanical properties.

The aim of this paper is to adequately analyze the microstructure
of dense colloids and to link microstructural differences to
differences in macroscopic mechanical properties upon applied
mechanical stress. Thereby, the focus is on the identification of
densely packed regions and chains of contacting particles in these
particle networks that may correlate differences in mechanical
properties with differences in microstructures. Our approach is to
analyze the distribution of densely packed regions in each
microstructure using the common neighbor distribution and the
dihedral angle distribution function. Additionally, a new method
called the straight path method, characterizing the quasi-linear
arrangement of particle chains, is introduced. The evaluation of
these characterization methods is performed with regard to their
ability to distinguish quantitatively between homogeneous and
heterogeneous microstructures obtained from
BD-simulations~\cite{Huetter_1999} mimicking well those found in
experiments~\cite{Wyss_2002}.

\section{\label{sec:MatAndMeth}Materials and Methods}

\subsection{Structure Generation}

The homogeneous and the heterogeneous microstructures are fully
coagulated colloidal suspensions obtained from Brownian dynamics
(BD) simulations~\cite{Huetter_1999}. The
Derjaguin-Landau-Verweg-Overbeek (DLVO) theory~\cite{Russel_1989}
was used to describe the particle-particle interaction, where
pair-wise particle potential interactions are assumed given by the
sum of the van der Waals attraction $V^{vdw}$
(Eq.~(\ref{eq:load_Vvdw})) and the electrostatic double layer
repulsion $V^{el}$ (Eq.~(\ref{eq:load_Vel})). Thus, $V^{dlvo} =
V^{vdw} + V^{el}$ with

 \begin{equation}
    V^{vdw}(r) = - \frac{A_H}{12}\left[\frac{d^2}{r^2-d^2} + \frac{d^2}{r^2} + 2 \ln \left(\frac{r^2-d^2}{r^2}\right)\right]
 \label{eq:load_Vvdw}
 \end{equation}

\noindent and

 \begin{equation}
    V^{el}(r) = \pi \epsilon_r \epsilon_{\mathrm{0}}
    \left[ \frac{4 k_b T}{z_V e} \tanh \left( \frac{z_V e}{4 k_b T} \Psi_{\mathrm{0}} \right)
    \right]^2
    d \exp \left( -\kappa \{ r-d\} \right),
 \label{eq:load_Vel}
 \end{equation}

\noindent respectively. The DLVO parameters are summarized in
Table~\ref{tab:load_BDparameter} and the potential curves for
various values of the surface potential $\Psi_{\mathrm{0}}$ are
shown in Fig.~\ref{fig:load_Fig3}.

\begin{table}
\caption{Potential Parameters for the Coagulating Suspensions}
\label{tab:load_BDparameter}
\begin{center}
\begin{tabular}{lll}
\hline\noalign{\smallskip}
Parameter & Symbol & Value  \\
\noalign{\smallskip}\hline\noalign{\smallskip}
Hamaker constant of Al$_2$O$_3$ in water  & $A_H$ & $4.76 \times 10^{-20}$ J\\
Particle diameter  & $d$ & 0.5~$\mu$m \\
Relative dielectric constant of water & $\epsilon_{r}$ & 81 \\
Absolute temperature  & $T$ & 293 K \\
Valency of ions  & $z_V$ & 1 \\
Inverse Debye screening length & $\kappa$ & 10$^{8}$ m$^{-1}$ \\
\noalign{\smallskip}\hline
\end{tabular}
\end{center}
\end{table}

\begin{figure}
  \includegraphics[width=\linewidth]{Fig3.eps}
\caption{DLVO interaction for different surface potentials: (a)
pure van der Waals attraction ($\Delta$pH-destabilization), (b)
$\Delta E = 0 k_B T$ and (c) $\Delta E = 5.65 k_B T$
($\Delta$I-destabilization).} \label{fig:load_Fig3}
\end{figure}

For $\Psi_{\mathrm{0}} = 0$ mV, the electrostatic double layer
repulsion is zero and the inter-particle potential is only given
by the attractive van der Waals potential. For $\Psi_{\mathrm{0}}
= 12$ mV a secondary minimum appears and $\Delta$E denotes the
energy barrier between the local maximum and the secondary
minimum. For $\Psi_{\mathrm{0}} = 15$ mV, a repulsive barrier of
$\Delta$E = 5.65 $k_B T$ exists and the secondary minimum is found
at $r = 1.08 d$. The model further contains the frictional Stokes'
drag force and a random Brownian force caused by the suspending
liquid.

We analyze microstructures generated using $\Psi_{\mathrm{0}} = 0$
mV, which correspond to the $\Delta$pH destabilization method and
$\Psi_{\mathrm{0}} = 15$ mV, corresponding to the $\Delta$I
destabilization method. In the following, the
$\Delta$pH-microstructures are labeled as homogeneous and the
$\Delta$I-microstructures as heterogeneous. All structures have
the same volume fraction of 0.4 and contain 8000 spherical
particles forming one percolating cluster. The particles have a
diameter of 0.5~$\mu$m.

\subsection{Structure Characterization Methods}

\subsubsection{Common Neighbor and Dihedral Angle Distribution
Function}

The common neighbor analysis~\cite{Clarke_1993} considers pairs of
particles, referred to as configurations, and determines the
number of particles that are in contact with both particles of the
configuration. In two dimensions, configurations can have at most
two common neighbors. In three dimensions, configurations can have
up to five common neighbors. The number of configurations with n
common neighbors is denoted $CN_n$ with $n=1,\ldots,5$. In this
study, absolute numbers of configurations are compared as the
microstructures analyzed have identical volume fraction, the same
number of particles and equal particle diameters.

\begin{figure}
  \includegraphics[width=\linewidth]{Fig4.eps}
\caption{Configuration (grey particles) with two common neighbors
(white particles) separated by a distance $r$ and having a
dihedral angle $\xi$. Perspective view (a) and top view (b).}
\label{fig:load_Fig4}
\end{figure}

The common neighbor distribution describes the short range
arrangement of particles. In particular, $CN_0$ is the number of
pairs of contacting particles with no common neighbor. $CN_1$
counts the number of configurations with exactly one neighbor
common to both particles, forming equilateral triangles. For $n$ =
2 the four particles of the configuration form a regular
tetrahedron if the two common neighbors are in contact, and a
generalized tetrahedron, having one longer edge, otherwise.

Configurations with two or more common neighbors ($CN_2,\ldots,
CN_5$) can be further characterized by the dihedral angle
distribution. This distribution analyses the arrangement of the
neighboring particles around the two particles of the
configuration, which form regular triangles with each of their
common neighbors. For $n \geq 2$ there are at least two common
neighbors. In this case, the dihedral angle $\xi$ is defined as
the angle between two planes spanned by triangles belonging to the
same configuration as is shown in Fig.~\ref{fig:load_Fig4}. $\xi$
is measured between one triangle and the triangles formed with
each other common neighbor of the original couple of particles.
Chemistry uses this method to characterize molecular
structures~\cite{Wielopolski_1986}. In the field of granular
matter, the dihedral angle distribution has been used to study the
structural organization and correlations in very large packings of
monodispersed spherical particles~\cite{Aste_2005,Aste_2006}.

For two common neighbors in contact with each other, the four
particles form a regular tetrahedron and the dihedral angle is
70.5$^{\circ}$. This angle constitutes the lower limit of the
dihedral angle distribution. The upper limit is 180$^{\circ}$ when
the four particles are in a plane and form a square.

Formally, the dihedral angle $\xi$ between two common neighbors of
a configuration having a separation distance $r$ between each
other is defined by $\xi = 2 \arcsin \left( \frac{r}{\sqrt{3} d}
\right)$. The dihedral angle distribution is given by $p(\xi +
\delta \xi / 2) = \delta N_{\xi} (\xi, \xi + \delta \xi)$ where
$\delta N_{\xi} (\xi, \xi + \delta \xi)$ is the number of
generalized (open) tetrahedra with a dihedral angle between $\xi$
and $\xi + \delta \xi$~\cite{Aste_2006}.

\subsubsection{Straight Path Method}

A straight path is a quasi-linear chain of contacting particles.
The shortest chain consists of three particles. In order for three
particles to form a straight path they have to fulfill the
following condition: The angle $\theta$ between the vector
connecting the first two particles (P$_1$, P$_2$) and the vector
connecting the second two particles (P$_2$, P$_3$) must be smaller
than a chosen threshold angle $\theta_c$. This is exemplified in
Fig.~\ref{fig:load_Fig5}.

\begin{figure}
  \includegraphics[width=\linewidth]{Fig5.eps}
\caption{Chain of three particles forming a straight path.}
\label{fig:load_Fig5}
\end{figure}

Analogously, a straight path of length $l$ is one of length $l-1$
that has another particle attached to one end and whose connection
vector fulfills the condition to stay in a direction within a cone
of angle $\theta_c$ from the vector connecting the previous two
particles. The distribution of straight path lengths is obtained
by determining the number of straight paths $SP_l$ of length
$l=3,\ldots,l_{max}$, with $l_{max}$ the longest path under
consideration. Absolute numbers are calculated and compared in the
distribution of $SP_l$.

The angle $\theta_c$ determines the quasi-linearity and specifies
the maximally allowed deflection from an absolute linear
arrangement of the particles. The dependence of the straight path
statistic on the angle $\theta_c$ is presented as part of the
results section and is used in order to choose a suitable angle
$\theta_c$. The distribution of straight paths for the homogeneous
and heterogeneous microstructures is compared. The average path
length $SP_{mean} = \Sigma_{l\geq 3} l SP_l / \Sigma_{l\geq 3}
SP_l$ and the number of paths longer than a chosen length $l_0$ is
determined. The latter is defined by $SP_{l\geq l_0} =
\Sigma_{l\geq l_0} SP_l$ ($l_0$ = 4 and 5) and permits to quantify
differences in the number of straight paths with longer lengths,
neglecting paths shorter than $l_0$ particles.

\section{Results and Discussion}
\label{ResAndDisc}

\subsection{Common Neighbor Analysis}

The number of particle configurations $CN_n$ is higher for the
heterogeneous than for the homogeneous microstructure for
$n=0,\ldots,3$ (Fig.~\ref{fig:load_Fig6}). This can easily be
understood as the mean coordination number of a particle in the
heterogeneous microstructure is with 5.2 roughly 10\% larger than
in the homogeneous microstructures with 4.7. When normalized by
the mean coordination number, the distributions collapse,
indicating that the relative distribution of common neighbor
configurations is identical for both microstructures. Thus, the
common neighbor distribution does not capture in a suitable way
the structural differences between these microstructures.

\begin{figure}
  \includegraphics[width=\linewidth]{Fig6.eps}
\caption{Common neighbor distribution for the homogeneous and
heterogeneous microstructures: number of configurations $CN_n$ in
dependence of the number of common neighbors $n$.}
\label{fig:load_Fig6}
\end{figure}

\subsection{Dihedral Angle Distribution}

The dihedral angle distributions of both microstructures show a
characteristic peak at 70.5$^{\circ}$ (Fig.~\ref{fig:load_Fig7})
corresponding to configurations that form regular tetrahedra with
two common neighbors. The peak height is the same for both
microstructures. The heterogeneous microstructure has a peak at
77$^{\circ}$, which is missing for the homogeneous one. This angle
corresponds to a separation distance between the common neighbors
of $1.08 d$ which is the distance between particles trapped in the
secondary minimum (Fig.~\ref{fig:load_Fig3}).

\begin{figure}
  \includegraphics[width=\linewidth]{Fig7.eps}
\caption{Dihedral angle distribution for the homogeneous and
heterogeneous microstructures. For clarity purposes, only every
fourth data point is indicated by a symbol.} \label{fig:load_Fig7}
\end{figure}

At $\xi = 109.5^{\circ}$ the homogeneous microstructure has a peak
that corresponds to a common neighbor separation of $1.41 d
\approx \sqrt{2}d$, which is absent in the case of the
heterogeneous microstructure. This separation distance corresponds
to two configurations forming a square and sharing two common
neighbors. An additional characteristic peak is found at
141$^{\circ}$ having roughly the same height for both
microstructures and corresponding to two juxtaposed tetrahedra.
The peak at 148$^{\circ}$ is solely present for the heterogeneous
microstructure and corresponds to one regular tetrahedron
juxtaposed to one with a distance of $1.08 d$ between the common
neighbors.

In conclusion, even if there are sizable differences between
homogeneous and heterogeneous structures, the absolute number of
regular tetrahedra is virtually equal for both microstructures and
the analysis does not reveal any structural difference that could
explain the better mechanical properties of the heterogeneous
microstructures. Thus, the common neighbor and the dihedral angle
distributions are not suited to link the microstructural
differences to the different mechanical properties.

\subsection{Straight Path Analysis}

The straight path analysis depends on the threshold angle
$\theta_c$, hence, a suitable value $\theta_c$ has to be
determined. For small angles $\theta_c$ only a small number of
straight paths are expected because only a few particle chains
satisfy the criterion to be a straight path. For large angles
$\theta_c$ many straight paths are found but the physical meaning
of a straight path is lost as the path becomes less ``straight''
and thus loses its load bearing capacity. Also for increasing
$\theta_c$ the ratio of the number of paths with equal lengths
between the heterogeneous and the homogeneous microstructure
decreases. Figure~\ref{fig:load_Fig8} shows the number of straight
path longer or equal to five particles for the homogeneous
($SP_{l\geq 5}^{HO}$) and for the heterogeneous microstructure
($SP_{l\geq 5}^{HE}$) as well as their ratio ($SP_{l\geq 5}^{HE} :
SP_{l\geq 5}^{HO}$) in dependence of the cone angle $\theta_c$
ranging from 10$^{\circ}$ to 40$^{\circ}$. The lines serve as
guide for the eyes.

\begin{figure}
  \includegraphics[width=\linewidth]{Fig8.eps}
\caption{Number of straight paths for both microstructures longer
or equal to five ($SP_{l=5}^{HE}$ and $SP_{l=5}^{HO}$, open
symbols) and their ratio ($SP_{l=5}^{HE} : SP_{l=5}^{HO}$,
triangles) as a function of the threshold angle $\theta_c$.}
\label{fig:load_Fig8}
\end{figure}

For $\theta_c$ between 14$^{\circ}$ and 18$^{\circ}$ the ratio
$SP_{l\geq 5}^{HE} : SP_{l\geq 5}^{HO}$ shows a maximum. However,
below 22$^{\circ}$ the number of straight paths is quite small
($SP_{l\geq 5}<100$) for both microstructures. Toward larger
angles the absolute number of straight paths is increasing while
the ratio $SP_{l\geq 5}^{HE} : SP_{l\geq l5}^{HO}$ is decreasing.
The angle $\theta_c = 30^{\circ}$ was chosen as this angle gives
reasonable values for both the ratio ($SP_{l\geq 5}^{HE} :
SP_{l\geq 5}^{HO} \approx 3$) and absolute number of paths
($SP_{l\geq 5}^{HE}= 451$ and $SP_{l\geq 5}^{HO}= 163$).

\begin{figure}
  \includegraphics[width=\linewidth]{Fig9.eps}
\caption{Number of straight paths $SP_l$ in dependence of the
straight path length $l_{SP}$ for the homogeneous and
heterogeneous microstructures (squares and circles, respectively)
and exponential fits (dashed and solid line, respectively).}
\label{fig:load_Fig9}
\end{figure}

For both microstructures the distribution of straight paths shows
an exponential decrease with straight path length $l_{SP}$
(Fig.~\ref{fig:load_Fig9}). Both microstructures have
approximately the same number of straight paths of length three.
Toward longer path lengths the microstructures present a diverging
behavior: the heterogeneous microstructure has more paths of
longer length than the homogeneous microstructure. The
distributions were fitted by an exponential law using $SP_l
\propto \exp (-\lambda l)$, with $\lambda$ the characteristic
inverse path length. A remarkably good fit of the data is obtained
for both microstructures indicated by the values close to one for
the correlation coefficients R$^2$. The higher exponent $\lambda$
found for the homogeneous microstructure demonstrates the more
rapid decrease of the number of straight paths toward longer path
lengths than in the case of the heterogeneous microstructure.

The values for the average path lengths ($SP_{mean}$) and the
numbers of paths longer or equal to four ($SP_{l\geq 4}$) and five
($SP_{l\geq 5}$) are summarized in Table~\ref{tab:load_SP}. The
average path length for the heterogeneous microstructure is only
5\% longer than for the homogeneous microstructure because
$SP_{l=3}$ dominates the statistics and is almost identical for
both structures. On the other hand, in the heterogeneous
microstructure longer chains $SP_{l\geq 4}$ and $SP_{l\geq 5}$ are
found 1.8 and 2.8 times more often, respectively, than in the
homogeneous one.

\begin{table}
\caption{Average path length $SP_{mean}$ and numbers of paths
longer or equal to four ($SP_{l\geq 4}$) and five ($SP_{l\geq 5}$)
for both microstructures} \label{tab:load_SP}
\begin{center}
\begin{tabular}{lll}
\hline\noalign{\smallskip}
 & Homogeneous & Heterogeneous \\
\noalign{\smallskip}\hline\noalign{\smallskip}
$SP_{mean}$ & 3.28 & 3.45 \\
$SP_{l\geq 4}$ & 779 & 1385 \\
$SP_{l\geq 5}$ & 163 & 451 \\
\noalign{\smallskip}\hline
\end{tabular}
\end{center}
\end{table}

To conclude, the number of straight paths in the homogeneous and
heterogeneous microstructures follows an exponential distribution
and a higher characteristic inverse path length $\lambda$ is found
for the homogeneous microstructure. The heterogeneous
microstructure contains significantly more straight paths of
longer lengths than the homogeneous one: in absolute numbers twice
as many paths with a length $\geq 4$ and three times as many
having a length $\geq 5$ are observed. These are very significant
differences between the two microstructures.

Our results suggest that the differences in the straight path
distributions are characteristic for the better mechanical
properties of the heterogeneous microstructures. Indeed, in
granular materials load is transmitted via force chains that are
quasi-linear substructures of the particle network. This has been
experimentally verified using photo-elasticity
techniques~\cite{Drescher_1972} and computationally modelled using
the discrete element method~\cite{Cundall_1979}. Also, the
computational study of the force chain network in 2D granular
assemblies of polydispersed particles subjected to indentation by
a rigid flat punch has shown that the force chain length follows
an exponential distribution~\cite{Peters_2005}.

\section{Conclusions}
\label{Conclusions}

Different chemical pathways during the internal destabilization of
colloidal suspensions lead to wet ceramic green bodies with
drastically different mechanical properties. The microstructures
of these destabilized colloids reveal differences in the
heterogeneity of the particle arrangements on a length scale in
the order of a few particle diameters. Various microstructural
characterization methods, i.e., the common neighbor distribution,
the dihedral angle distribution and the straight path method, have
been applied in order to analyze homogeneous and heterogeneous
microstructures generated by Brownian dynamics simulations, which
were shown to agree well with the experimentally determined
microstructures of such coagulated colloidal particle systems.

Toward our goal, to establish a correlation between the
microstructural differences and the differences in macroscopic
mechanical properties, it was found that the common neighbor and
the dihedral angle distributions do not discriminate enough the
differences between the homogeneous and heterogeneous
microstructures and can therefore hardly account for the large
differences in mechanical properties. The newly introduced
straight path method reveals significantly more straight paths of
longer lengths for the heterogeneous than for the homogenous
microstructure: twice as many straight paths of length $\geq 4$
particles and three times as many of length $\geq 5$ particles are
found in the heterogeneous microstructure. These differences in
straight path number and length seem suitable to characterize and
to differentiate between these structures and these differences
are considered to be sufficiently large to account for the
differences in mechanical properties. The straight paths may
capture best the characteristic microstructural features which are
relevant for the mechanical properties. The number of straight
paths was found to follow an exponential distribution just like
the distribution of the force chain lengths in mechanically loaded
granular material~\cite{Peters_2005}. Additionally, the
quasi-linear structure of the straight paths seems to correspond
to the geometrical shape of force
chains~\cite{Drescher_1972,Cundall_1979} that are well known to
determine the load bearing capacity of granular matter. These
findings are encouraging for an attempt to establish a correlation
between defined microstructural features and the differences in
mechanical properties of destabilized colloids.

Further quantitative structural analyses and computational
modelling using the discrete element method are in progress to
study the force chains in ceramic green bodies and to correlate
them with the straight paths of their microstructures.

\begin{acknowledgments}
The authors would like to express their gratitude to Markus
H\"utter for providing the colloidal microstructure data from the
Brownian dynamics simulations.
\end{acknowledgments}

\newpage

Soft Condensed Matter (cond-mat.soft); Materials Science (cond-mat.mtrl-sci)


Plerionic Supernova Remnants

Abstract: Plerions represent ideal laboratories for the search for neutron stars, the study of their relativistic winds, and their interaction with their surrounding supernova ejecta and/or the interstellar medium. As well, they are widely believed to represent efficient engines for particle acceleration up to the knee of the cosmic ray spectrum (at about 1E15 eV). Multi-wavelength observations from the radio to the highest TeV energies, combined with modelling, have opened a new window to study these objects, and particularly shed light on their intrinsic properties, diversity, and evolution. High-resolution X-ray observations are further revealing the structure and sites for shock acceleration. The missing shells in the majority of these objects remain puzzling, and the presence of plerions around highly magnetized neutron stars is still questionable. I review the current status and statistics of observations of plerionic supernova remnants (SNRs), highlighting combined radio and X-ray observations of a growing class of atypical, non Crab-like, plerionic SNRs in our Galaxy. I will also briefly describe the latest developments to our high-energy SNRs catalogue recently released to the community, and finally highlight the key questions to be addressed in this field with future high-energy missions, including Astro-H in the very near future.

\section{Plerionic Supernova Remnants}

Supernova Remnants (SNRs) are generally classified under three categories: Shells, Filled-centre, or Composites\footnote{see Dave Green's catalogue, http://www.mrao.cam.ac.uk/surveys/snrs/,
for a compilation and classification of the Galactic Supernova Remnants based on radio studies.}.
Filled-centre SNRs have been referred to in the literature  as `plerions' or  Pulsar Wind Nebulae (PWNe), and Composite-type SNRs are classified as either Plerionic Composites (those powered by plerions) or thermal Composites (also referred to as `Mixed-Morphology' SNRs).
The term `plerion' originates from the ancient greek word pleres ($\pi \lambda \eta \rho \eta s$) \cite{weilerpanagia78}
which means `filled'; the centrally filled morphology arising from the non-thermal emission powered by a neutron star.

Typified by the Crab nebula, plerions are bubbles of relativistic particles inflated by the pulsar's relativistic wind.
As such they are  `calorimeters' for the study of neutron stars and, in the absence of a pulsar detection, excellent pathfinders for pulsar discovery.
They allow us to probe the interaction between the neutron stars' relativistic winds and their surroundings,
which consist of the supernova ejecta at the earlier stages of SNR evolution, the shocked interstellar matter or supernova blast wave as the neutron star encounters the expanding SNR shell,
and the interstellar medium (ISM) after the neutron star has left its hosting SNR.
Plerions seed the Galaxy with highly energetic particles and strong magnetic fields, and are thus considered efficient engines for cosmic ray acceleration up to TeV energies.

 Plerions are observed from radio to the highest energy gamma-rays, and  are {\it generally} characterized by the following properties:
1) an increase in the brightness towards the centre, unlike shells which have a limb-brightened morphology;
2) their X-ray size is generally shorter than
their radio and optical size, due to the smaller synchrotron lifetimes of the
higher energy photons;
3) a flat radio spectral index, $\alpha$ = 0--0.3 ($S_{\nu}$$\sim$$\nu^{\-\alpha}$) with
 a steeper index in the X-ray due to synchrotron and radiation losses and with an average photon index
$\Gamma$~=~$\alpha$+1$\sim$2; and
4) they are highly polarized with a preferentially radial magnetic
field whose strength varies from a few $\mu$G -- mG.

Modern X-ray telescopes, combined with numerical simulations, have greatly advanced our understanding of plerions \cite{gaenslerslane2006, KP2008, KP2010, bucciantini11}. In particular, they have allowed us to
1) zoom in on the interaction of the relativistic pulsar wind with its surroundings, showing evidence of torus and jet-like structures successfully reproduced in simulations,
and 2) map the photon index across the plerion, answering questions related to the energetics of the powering engine and the wind magnetization parameter at the termination shock.
Correlating the X-ray with radio and gamma-ray studies sheds light on the magnetic field configuration and strength, the conditions of the confining ambient medium, the age and the distance to the object.
While their radio to X-ray emission arises from synchrotron radiation from relativistic particles injected by the neutron star.
 their gamma-ray emission arises, in the leptonic model, from Inverse Compton scattering (with the seed photons coming from the synchrotron emitting relativistic particles in the pulsar wind, or from the cosmic microwave background, dust and starlight), and in the hadronic model from pion decay.

\subsection{Evolution of Plerions}
The evolution of a plerion in an SNR has been studied by several authors, including \cite{blondin2001, bucciantini2003, vanderswaluw2004}.
Initially plerions expand in the freely moving ejecta of the SNR. As the SNR evolves from the free expansion phase into the Sedov-Taylor phase, a reverse shock develops that runs back into the expanding ejecta and can eventually reach the central plerion
on a timescale of $\sim$10 kyr. This is followed by a crushing and re-expansion phases of the plerion.
This action displaces the plerion from the pulsar, leaving a `relic' nebula which can still be detected in the radio.
The `crushing' phase is subject to Rayleigh-Taylor instabilities that result in the mixing of thermal (ejecta) and non-thermal (PWN) fluids.
Furthermore, asymmetries in the position of the plerion relative to the pulsar and explosion site arise from asymmetries in the surrounding ISM.
Meanwhile the continuous injection of relativistic particles from the neutron star forms an X-ray nebula surrounding it and that is offset from the relic nebula.
The pulsar's motion inside the SNR further causes this wind bubble of freshly injected particles to move with it through the SNR interior.
When its motion becomes supersonic, it forms a bow-shock nebula.

The best and most studied example for a remnant displaying both a bow-shock and relic nebulae interacting with an asymmetrical reverse shock  is the Vela-X region in the $\sim$11-kyr old SNR Vela.
The complex morphology and properties of the nebula in the radio \cite{dodson2003}, X-rays \cite{MO95, pavlov2003, lamassa08} and gamma-rays \cite{aharonian2006} suggest
that a reverse shock has displaced the original plerion (Vela X) from the pulsar.
While Vela X is believed to be a relic nebula, the Vela pulsar has created a new nebula in its vicinity that clearly has a bow-shock morphology \cite{blondin2001, CR2011}.
Another example of an offset plerion likely interacting with the reverse shock is in the Composite-type SNR G327.1$-$1.1 \cite{temim2009},
 estimated to have an age of $\sim$18~kyr and recently discovered in TeV gamma-rays with H.E.S.S. \cite{acero12}.
The overall morphology is also similar to that seen in other offset TeV PWNe, such as HESS J1825$-$139 \cite{aharonian2006b}.

 \section{Statistics and the first X-ray plus gamma-ray SNR catalogue}

As of October 2012, we know of 309 Galactic SNRs; see \cite{FSSH2012} for a
an up-to-date (and first of its kind) catalogue\footnote{http://www.physics.umanitoba.ca/snr/SNRcat} of X-ray and gamma-ray observations
of all Galactic SNRs. Our high-energy catalogue complements Dave Green's catalogue based on radio studies.
 Out of these 309 remnants, $\sim$90 SNRs contain plerions or plerionic candidates,  about
30 of which lack shells (including the Crab). These will be referred to as `naked' PWNe.
 Out of the 70 Galactic TeV known sources, about half are identified as PWNe or PWN candidates.
 In the majority of the latter class, the TeV emission is found to be offset from the X-ray peak suggestive of `relic' PWNe (like for Vela-X).
 The TeV detections are particularly useful for studying low-magnetic field plerions, since these could be X-ray faint but gamma-ray bright.

\section{Growing Class of `unusual' Plerions}

For over three decades, our view of plerions has been shaped by the study
of the young and bright Crab nebula. However, there are several plerions whose properties
differ from the Crab, hinting to a diversity of PWNe which include many previously missed evolved objects.
In the following, I highlight
some interesting cases that we have been studying using coordinated radio and X-ray observations,
and that need to be further investigated in gamma-rays and with deeper X-ray observations.

\begin{figure}
 \includegraphics[height=.5\textheight]{da495_g76_v3_compressed.eps}
 \caption{Top: The radio images of the plerionic SNRs DA~495 (left) and G76.9+1.1 (right), with the arrow pointing to the location of the X-ray source discovered with $Chandra$.
 The bottom panels show the $Chandra$ images of the central regions of the SNRs highlighting the discovery of compact X-ray nebulae, inside the depression seen in the radio, surrounding pulsar candidates.
 For  DA~495, the X-ray nebula is $\sim$40$^{\prime\prime}$ in diameter surrounding the point-like source CXOU J202221.68+384214.
 For G76.9+1.0, the X-ray image has a scale size of 1$^{\prime}$$\times$1$^{\prime}$ and the central ellipse shows the compact nebula with a size of 8$^{\prime\prime}$.3$\times$5$^{\prime\prime}$.1
 surrounding PSR~J2022+3842.
Adapted from figures and data published in \cite{Arzoumanian2008}, \cite{Kothes2008}, and \cite{Arzoumanian2011}.}
 \end{figure}

\subsection{DA~495 (G65.7+1.2) and G76.9+1.0}

The SNRs DA~495 (G65.7+1.2) and G76.9+1.0 (see Fig.~1) are unusual due to the following reasons:
1) unlike shell-type SNRs, they are not limb-brightened, but their non-thermal flux falls off away from centre, implying a neutron star as a central engine;
2) their radio spectrum has a very steep spectral index, $\alpha$ $\sim$ 0.6, which is more typical of Shell-type remnants;
3) they are large and have a low-frequency spectral break in the radio; and
4) in polarized emission, they have a strong axisymmetry and a double-lobed morphology.

\begin{table}[tbh]
\begin{tabular}{lrr|rr}
\hline
 Component & DA495 Nebula & pulsar candidate & G76.9+1.0 Nebula & PSR~J2022+3842 \\
\hline
$N_{\rm H}$ & $\leq$0.4$\times$10$^{22}$~cm$^{-2}$ & ... & 1.6$\pm$0.3$\times$10$^{22}$~cm$^{-2}$ & ... \\
$\Gamma$/kT (MK) & 1.6 & 2.5 & 1.4 & 1.0$\pm$0.2\\
$L_x$\tablenote{in units of erg~s$^{-1}$, and at a distance of 1~kpc for DA~495 and 10 kpc for G76.9+1.0} & 3.3$\times$10$^{31}$ & $\sim$10$^{31}$ & 5.6$\times$10$^{32}$  & 7$\times$10$^{33}$ \\
\hline
\end{tabular}
\caption{Summary of the properties of the plerion and pulsar candidate/pulsar in the SNRs DA~495 and G76.9+1.0 \cite{Arzoumanian2008, Arzoumanian2011}. All models are power-law except for the point source in DA~495 which is fitted with a blackbody model.}
\end{table}

While their nature as SNRs has been questioned in the past, $Chandra$ observations have allowed the discovery of compact non-thermal nebulae surrounding point-like sources, sitting in the heart of the radio depression. The morphology and spectral properties of both the nebulae and point-like sources (see Table~1) left no doubt about their nature as plerions. Combined X-ray and radio studies have further suggested that these are evolved PWNe with a very low X-ray to radio luminosity and size ratios, suggestive of a high nebular magnetic field~\cite{Kothes2008}.
Follow-up radio pulsation searches failed to detect a pulsar in DA~495, but led to the discovery of a 24 ms pulsar (PSR J2022+3842)  in G76.9+1.0, a result that was confirmed in X-rays with the Rossi X-ray Timing Explorer \cite{Arzoumanian2011}.
While DA~495 hasn't been yet detected in gamma-rays, there has been recently a reported Fermi detection in G76.9+1.0 \cite{nolan12}.

\subsection{CTB~87 (G74.9 +1.2) and G63.7+1.1}

\begin{figure}[tbh]
 \includegraphics[height=.35\textheight]{ctb87figv3_compressed.eps}
\caption{Left Panel: The \textit{Chandra} greyscale X-ray image of CTB~87 with the Canadian Galactic Plane Survey radio contours overlaid (in red). The top right panel shows the RGB energy image of the extended diffuse nebula
with a cometary morphology powered by a pulsar candidate. The bottom panel shows a more compact nebula surrounding the putative pulsar with evidence of faint jet-like and torus-like structures. See \cite{MSHK2012} for details.}
\end{figure}

CTB~87 and G63.7+1.1 are two plerionic candidates, previously unexplored in X-rays, but that also have unusual characteristics in radio wavelengths although with a morphology different from the plerions discussed in the previous section.
They both have a centrally filled morphology (as for plerions) but  they are also unusual in that they have 1) a steep radio spectral index ($\alpha$~=~0.4--0.5), 2) a low-frequency spectral break in the radio, and 3) they are unusually large implying evolved PWNe.
Our $Chandra$ and \textit{XMM-Newton} studies unveiled faint X-ray PWNe in both objects, surrounding pulsar candidates, and that are clearly offset from the peak of the radio emission (papers in preparation). In CTB~87 shown in Fig.~2, we find evidence of both a compact nebula with a torus/jet-like structure and a diffuse nebula with a cometary-like morphology \cite{MSHK2012}. The radio and X-ray study suggests that CTB~87 represents an evolved ($\sim$10--20~kyr old) plerion where the radio emission peak corresponds to a relic PWN, similar to what's observed in G327.1--1.1.
Recent gamma-ray observations with VERITAS led to the detection of this source \cite{aliu11}.

\subsection{Plerions around highly magnetized neutron stars?}

\begin{figure}[tbh]
 \includegraphics[height=.3\textheight]{g292j1119_compressed.eps}
  \caption{ \emph{Right}: The compact PWN and jet discovered with \textit{Chandra} surrounding the high-magnetic field pulsar J1119$-$6127~\cite{sshkumar08} in the SNR G292.2$-$0.5 shown in the
  left panel. \emph{Left:}   \textit{XMM-Newton} image and contours in greyscale, with the \textit{Chandra} contours overlaid in blue (darker contours). The red contours show the optical/DSS emission
  illustrating the presence of dark clouds in the \textit{eastern} field of G292.2$-$0.5. The purple cross in the \textit{western} field refers to the H.E.S.S. detection in the SNR interior \cite{djannati09}. Figure adapted from figures published in \cite{kumarsshgonz12} and \cite{sshkumar08} and best seen in colour in the electronic version.}
\end{figure}

The question on whether highly magnetized neutron stars power plerions remains to be answered. So far among the known 23 magnetars plus 7 X-ray detected high-magnetic radio field pulsars,
there are only 4 objects confirmed to power PWNe discovered or detected in X-rays:
the high-magnetic field pulsars PSR J1119$-$6127 in the SNR G292.2$-$0.5~\cite{sshkumar08, gonzssh03} (see Fig.~3) and
PSR J1846$-$0258 in the SNR Kes~75~\cite{kumarssh08, gavriil08, ng08},
PSR J1819-1458 (classified as a RRAT or Rotating Radio Transient) \cite{rea09}, and the magnetar Swift J1834.9$-$0846~\cite{younes12}.
The PWN claimed to be associated with the anomalous X-ray pulsar 1E 1547.0$-$5408 \cite{bambavink09} was recently proposed to be a dust scattering halo \cite{olausen11}.
Among these, Kes 75 and G292.2$-$0.5 have been detected in gamma-rays with H.E.S.S.\footnote{http://www.mpi-hd.mpg.de/hfm/HESS/pages/home/som/2008/10/}$^,$\footnote{http://www.mpi-hd.mpg.de/hfm/HESS/pages/home/som/2009/11/}.
In G292.2$-$0.5, the H.E.S.S. source was located in the western field of the SNR, between the pulsar and the SNR shell where hard non-thermal X-ray emission was also detected (Fig.~3).
While the origin of the gamma-ray emission is not yet understood, it could be attributed to an offset plerion, where the nebula has been displaced after being crushed by an asymmetric reverse shock
\cite{kumarsshgonz12, djannati09}.

\begin{figure}
 \includegraphics[height=.3\textheight,width=1.0\textwidth]{g21g327ima_compressed.eps}
  \caption{\emph{Left}: The young plerionic SNR G21.5$-$0.9 with $Chandra$ showing evidence for thermal emission from `knots' outside the main PWN component, and a limb-brightened shell to the east. Image credit: \cite{MSSH2010}. \emph{Right}: The composite $Chandra$ (in blue) and radio (in red) image of the Composite-type SNR G327.1--1.1 harbouring an offset plerion recently detected in gamma-rays. Image credit: \cite{temim2009}.}
\end{figure}

\section{Conclusions and future prospects}

The synergy between radio and high-energy observations has been crucial to understanding the emission mechanism in plerions and in addressing
their apparent diversity. In particular, radio plus X-ray observations are allowing the discovery of a new class of plerions which seem to have unusual properties,
but that likely represent evolved objects.
Some of the outstanding questions in this field include:
\begin{itemize}
\item{Where is the shell around the dozens of `naked' plerions? Is the absence of a shell caused by the nature of the progenitor star creating a wind-blown bubble in which the plerion is expanding?
Are they born in lower energy supernova explosions?}
\item{Do all high magnetic field pulsars and magnetars power nebulae?  What powers their X-ray emission: rotation, magnetism, both, or some other mechanism?}
\end{itemize}

Deep observations with current observatories as well as targeted multi-wavelength
observations will address these, and other, questions related to plerion physics and evolution.
In X-rays, we look forward to the launch of the Astro-H mission \cite{takahashi2010} which will have an unprecedented spectral resolution and sensitivity, as well as a broadband coverage in the 0.5--600 keV energy range.
This will allow the search for the long sought SNR shells and thermal plasma in the `naked' plerions, and for the thermal emission expected in evolved plerions crushed by the reverse shock (see Fig.~4).
The unique broadband coverage will also measure their broadband X-ray spectra and, in conjunction with radio and gamma-ray studies, help constrain their magnetic field and emission mechanism (leptonic versus hadronic).
Finally, the synergy between X-ray and gamma-ray observations with current and the new generation observatories (including Fermi,  H.E.S.S. II,  VERITAS, and the Cerenkov Telescope Array in the future)
will allow us to complete the census of Plerionic SNRs in our Galaxy and beyond.

\begin{theacknowledgments}
 The author acknowledges the support of the Canada Research Chairs program, the Natural Sciences and Engineering Research Council of Canada,
the Canadian Institute for Theoretical Astrophysics, the Canadian Space Agency, and the Canada Foundation for Innovation.
She is grateful to the organizers for the invitation and the opportunity to attend a stimulating and well organized conference
in a nice city. She also acknowledges the many contributions to the results presented in this review by her
collaborators particularly Zaven Arzoumanian, Roland Kothes and Tom Landecker; group members
including Heather Matheson, Gilles Ferrand, and Harsha Kumar; as well as the entire active and productive SNR/PWN community.
\end{theacknowledgments}

\bibliographystyle{aipprocl}

Astrophysics of Galaxies (astro-ph.GA); High Energy Astrophysical Phenomena (astro-ph.HE)

Plerionic Supernova Remnants


Intermediate Scale Structures in BL Lac objects

Abstract: The parsec-scale total intensity structures of BL Lac objects share many characteristics with quasars: e.g. well-collimated jets and superluminal motion. However, the jets of BL Lac objects appear to fade much more quickly than those in quasars and on VLA scales many BL Lacs have structures comparable to those of low luminosity (FR I) radio galaxies, in which kpc-scale flow speeds are non-relativistic.The region between parsec and kiloparsec scales is therefore one of transition in BL Lac objects. We have carried out VLBI observations of a small sample of BL Lac objects at frequencies between 2.3 GHz and 327 MHz in order to investigate this transition region. Preliminary results for two of these sources are presented here.Coherent structures are detected in these sources out to a distance of several tens of parsecs. We find evidence for significant changes in jet structure on scales of a few tens of milliarcseconds, both in terms of the orientation and collimation of the jets.

\section{Introduction}
BL~Lac objects are a sub-class of extragalactic radio sources exhibiting highly
variable polarized radio emission and are distinguished from optically
violently variable quasars by their weak optical line emission (equivalent
width less than 5 \AA, \citealp{stickel93}).

The parsec (pc) and kpc scale jets of BL~Lac objects are often misaligned and
the kpc scale emission is typically more diffuse. This indicates that the jet
has become decollimated between the milliarcsecond (mas) and arcsec scales.
There is evidence in some BL~Lac objects for significant total intensity and
polarized flux density on scales intermediate to the mas and arcsec resolution
provided by VLBI and VLA observations respectively \citep{gabuzda94}. Further,
the pc scale jets of BL~Lac objects exhibit inferred magnetic fields which tend
to be perpendicular to the direction of the jet motion (e.g.
\citealp{gabuzda99} and references therein), while on arcsecond scales the
magnetic field is often aligned with the jet \citep{kollgaard92}. The region
over which these changes take place lies on scales between those accessible to
high-frequency ($> 5$~GHz) VLBI and more conventional interferometers such as
the VLA, and has therefore been little explored to date.

In order to investigate these changes of structure with scale, we have carried
out observations at frequencies of 2.3~GHz, 1.7~GHz, 610~MHz and 327~MHz with
the VLBA. At the lower frequencies more diffuse jet emission should become more
dominant allowing us to trace the development of the jet on larger scales.

Linear sizes and luminosities have been calculated using $H_0 = 75$ km s$^{-1}$ Mpc$^{-1}$, $q_0=0.05$.

\section{Observations}
The observations took place on 13th December 1999 using all ten elements of the
VLBA in a dual polarization mode at a recording rate of 128 Mb s$^{-1}$. The
data were correlated at the VLBA correlator in Socorro, New Mexico. The
observing frequencies were 2.3~GHz, 1.7~GHz, 610~MHz and 327~MHz. The 2.3~GHz
and 1.7~GHz data were processed in full polarization mode, but the polarization
maps are not presented here. The a priori calibration and self-calibration of
the data were performed using the standard programmes in the NRAO Astronomical
Image Processing Software (\texttt{AIPS}) package.

In total, a small sample of 5 objects were observed, maps for two of which are
presented here.

\section{Results}
The observed sources were inferred to have significant flux density on
intermediate scales from a comparison of their flux densities in previous VLA
and VLBI observations at a frequency of 5~GHz. The expected 5~GHz intermediate
scale flux densities of the two sources are presented in
Table~\ref{table:isflux}. For both of these sources we have been able to image
significant emission on larger scales than has been reported in previous
observations at higher frequencies.

\begin{table}
\caption[]{Intermediate scale (I.S.) flux densities inferred from
previous 5~GHz VLA and VLBI images \citep{gabuzda94,gabuzda92}.}

\label{table:isflux}
$$
\begin{array}{p{0.5\linewidth}p{0.4\linewidth}}
\hline
\noalign{\smallskip}
Source      &  5~GHz I.S. flux density (Jy) \\
\noalign{\smallskip}
\hline
\noalign{\smallskip}
1219+285  &   0.22    \\
0735+178  &   0.87    \\
\noalign{\smallskip}
\hline
\end{array}
$$
\end{table}

\subsection{1219+285}
1219+285 lies in the centre of an elliptical galaxy at a redshift of 0.102
\citep{weistrop85}. This source is essentially unresolved to the VLA D array
\citep{kollgaard92}. The VLBI jet extends to the east of the core and a number
of superluminal components have been identified \citep{gabuzda94}.

The images in Figures~\ref{fig:1219+285_2.3ghz} to \ref{fig:1219+285_0.3ghz}
show a large diffuse extension towards the south roughly perpendicular to the
known VLBI jet. There is an apparent change in jet direction visible in the
2.3~GHz data. It is possible that the jet is changing direction from an
easterly to a southerly direction and finally in a westerly direction, while
also becoming more diffuse. If the jet is curving towards the  line of sight
then it is also possible that this more diffuse emission from the southern part
of the source represents emission from the end of a highly projected jet.
However, it is also possible that there is a qualitative change in the jet
structure and that the jet itself is becoming diffuse and decollimated. It is
also of note that there is a region of enhanced brightness and relatively flat
spectral index (perhaps resulting from free-free absorption of a steep
spetctrum jet, see Figure~\ref{fig:1219+285_specind}) in approximately the
position that the jet direction apparently begins to change from an easterly to
a southerly direction. This may be indicative of an interaction between the jet
and the surrounding medium which results in the jet changing direction. Further
analysis of the associated polarization data should provide more information on
this point.

The expected intermediate scale flux density from comparison of VLA and VLBI
images is 0.22~Jy at 5~GHz. The 2.3~GHz image presented here detects
approximately 0.11~Jy outside the region which was observed in the previous
5~GHz VLBI observations \citep{gabuzda94}. Even allowing for spectral index
effects, this represents a large fraction of the missing 5~GHz flux density,
although it also suggests that there is significant emission on still larger
scales than these observations are sensitive to. It is interesting that no
significant new emission becomes visible in the lower frequency 610 and
327~MHz images when compared with the 1.7 and 2.3~GHz images. The sources
appear essentially unresolved at these lower frequencies.

At the redshift of this object, 1 mas $ \sim 1.71$ pc. The luminosity detected
in the new, extended structure in this observation is $\sim 2.2 \times 10^{24}$
W Hz$^{-1}$.

\begin{figure}
\centering
\includegraphics[width=6.5cm]{creynolds_fig1.ps}
   \caption{2.3~GHz image of 1219+285. The peak flux density is 0.38 Jy/beam
and the contour levels are $-0.36$, 0.36, 0.72, 1.44 ... 92.2 \%.
           }
      \label{fig:1219+285_2.3ghz}
\end{figure}

\begin{figure}
\centering
\includegraphics[width=6.5cm]{creynolds_fig2.ps}
   \caption{1.7~GHz image of 1219+285. The peak flux density is 0.38 Jy/beam
and the contour levels are $-0.36$, 0.36, 0.72, 1.44 ... 92.2 \%.
           }
      \label{fig:1219+285_1.7ghz}
\end{figure}

\begin{figure}
\centering
\includegraphics[width=6.5cm]{creynolds_fig3.ps}
   \caption{610~MHz image of 1219+285. The peak flux density is 0.61 Jy/beam
and the contour levels are $-0.8$, 0.8, 1.6, 3.2 ... 51.2 \%.
           }
      \label{fig:1219+285_0.6ghz}
\end{figure}

\begin{figure}
\centering
\includegraphics[width=6.5cm]{creynolds_fig4.ps}
   \caption{327~MHz image of 1219+285. The peak flux density is 0.66 Jy/beam
and the contour levels are $-0.45$, 0.45, 0.9, 1.8 ... 57.6 \%.
           }
      \label{fig:1219+285_0.3ghz}
\end{figure}

\begin{figure}
\centering
\includegraphics[width=6.5cm]{creynolds_fig5.ps}
   \caption{The contours are of the 2.3~GHz image shown in
Figure~\ref{fig:1219+285_1.7ghz}. The grey scale is the spectral index
measured between 2.3 and 1.7~GHz.
           }
      \label{fig:1219+285_specind}
\end{figure}

\subsection{0735+178}
No emission lines have yet been identified in the optical spectrum of 0735+178,
though a lower limit to the redshift of 0.424 is provided by a pair of
absorption lines \citep{carswell74}. VLBI images of this source show a compact
core and a jet that extends to the northeast (e.g. \citealp{gabuzda92}) while
arcsecond scale VLA images reveal a point source
\citep{ulvestad83,kollgaard92}. 8, 22 and 43~GHz VLBI images reveal a jet out to
approximately 10~mas with two sharp apparent bends of 90$^{\circ}$ within 2~mas
of the core. It also appears that there was a dramatic change in the trajectory
of the jet within the first few mas between 1992 and 1995
\citep{gomez99,gomez01b}. This could be interpreted either as a precessing jet
nozzle or more likely as a result of pressure gradients in the external medium
through which the jet propagates \citep{gomez01b}.

The images presented here show, for the first time, coherent structures in this
source on scales of several tens of mas. The sharp twists close to the core
reported at higher frequencies are not seen here owing to the low resolution,
but the jet appears to turn from its initial easterly direction towards the
south at a distance of about 30~mas from the core. The beginning of this bend
is apparent in the 2.3~GHz, 1.7~GHz and 610~MHz images. In the 327~MHz image
there is a suggestion of emission 100~mas due south of the core which may
indicate that the jet continues on its southerly path after the initial turn at
30~mas from the core, although the detection is not strong enough to confirm
this.

The flux density in the previously undetected structures, beyond 10~mas from
the core, is 0.122 Jy, which, after accounting for the spectral index effects,
represents a very large fraction of the predicted intermediate scale flux at
5~GHz of $\sim 0.090$~Jy \citep{gabuzda94}.

At the lower limit for the redshift of this object, 1 mas $ \sim 4.9$ pc. The
luminosity detected in the new, extended, structure in this observation is
therefore $> 6.1 \times 10^{25}$ W Hz$^{-1}$.

\begin{figure}
\centering
\includegraphics[width=6.5cm]{creynolds_fig6.ps}
   \caption{2.3~GHz image of 0735+178. The peak flux density is 0.78 Jy/beam
and the contour levels are $-0.2$, 0.2, 0.4, 0.8 ... 51.2 \%.
           }
      \label{fig:0735+178_2.3ghz}
\end{figure}

\begin{figure}
\centering
\includegraphics[width=6.5cm]{creynolds_fig7.ps}
   \caption{1.7~GHz image of 0735+178. The peak flux density is 0.79 Jy/beam
and the contour levels are $-0.2$, 0.2, 0.4, 0.8 ... 51.2 \%.
           }
      \label{fig:0735+178_1.7ghz}
\end{figure}

\begin{figure}
\centering
\includegraphics[width=6.5cm]{creynolds_fig8.ps}
   \caption{610~MHz image of 0735+178. The peak flux density is 0.65 Jy/beam
and the contour levels are $-1.5$, 1.5, 3.0, 6.0 ... 96.0 \%.
           }
      \label{fig:0735+178_0.6ghz}
\end{figure}

\begin{figure}
\centering
\includegraphics[width=6.5cm]{creynolds_fig9.ps}
   \caption{327~MHz image of 0735+178. The peak flux density is 0.43 Jy/beam
and the contour levels are $-1.2$, 1.2, 2.4, 4.8 ... 76.8 \%.
           }
      \label{fig:0735+178_0.3ghz}
\end{figure}

\section{Discussion}
It has been suggested that BL~Lac objects are low luminosity (FR\,I) radio
galaxies orientated such that Doppler boosting enhances their apparent
brightness (e.g. \citealt{browne83}). If this model is valid, the extended
structure in BL~Lacs should have a similar radio power to that found in the
extended structure of FR\,I radio galaxies. The FR\,I/FR\,II luminosity
division is $\sim 10^{24.5}$ -- $10^{25}$ W Hz$^{-1}$ at 1.4 GHz
\citep{bridle84b}. Allowing for a spectral index of $\sim -0.7$, this would
correspond to a 2.3~GHz division at approximately $\sim 2.0 \times 10^{24}$ --
$6.3 \times 10^{24}$ W Hz$^{-1}$.

For the two objects presented here the detected intermediate scale flux density
was a large fraction of that expected from the discrepancy between previous
5~GHz VLA and VLBI images. In the case of 0735+178, this may well account for
almost all of the intermediate scale ``missing" flux. The intermediate scale
flux detected for 1219+285 is somewhat less than the predicted intermediate
scale flux which means that more flux may lie on larger scales than are probed
in this observation.

The luminosity of the intermediate scale emission detected in 1219+285 is $\sim
2.2 \times 10^{24}$ W Hz$^{-1}$. This emission appears to be rather diffuse,
and might be associated with a plume-like structure formed by the disruption of
the well-collimated parsec-scale jet seen closer to the core. Even if all of
this emission is associated with a diffuse unbeamed component, it is at the
lower end of the transition region between FR\,I and FR\,II radio sources, and
so does not pose a problem for associating the diffuse emission in this galaxy
with FR\,I structure.

The intermediate scale flux density detected in 0735+178 appears to be
associated with a well-collimated jet, and is therefore likely to be
significantly Doppler beamed. For this reason, it is probably not significant
that the observed intermediate scale luminosity is above the FR\,I limit.

\section{Conclusions}
\begin{enumerate}
\item
In the two sources presented here, coherent structures extend out to
projected distances of $> 30$~pc from the core before disruption.
\item
In one case, 1219+285, diffuse -- possibly plume-like -- emission was
detected.
\item
The luminosity detected in the intermediate scale structure of 1219+285
presented here was found to be at the lower end of the
FR\,I/FR\,II luminosity division. 0735+178 had a luminosity in the
extended structure above this division, but as the emission appears to be
associated with a jet rather than a diffuse component this does not necessarily
imply that the source contains FR\,II extended structure.
\end{enumerate}

\begin{acknowledgements}
The National Radio Astronomy Observatory is operated by Associated Universities
Incorporated, under cooperative agreement with the NSF.
\end{acknowledgements}