Dynamics and structure of an aging binary colloidal glass

Abstract: We study aging in a colloidal suspension consisting of micron-sized particles in a liquid. This system is made glassy by increasing the particle concentration. We observe samples composed of particles of two sizes, with a size ratio of 1:2.1 and a volume fraction ratio 1:6, using fast laser scanning confocal microscopy. This technique yields real-time, three-dimensional movies deep inside the colloidal glass. Specifically, we look at how the size, motion and structural organization of the particles relate to the overall aging of the glass. Particles move in spatially heterogeneous cooperative groups. These mobile regions tend to be richer in small particles, and these small particles facilitate the motion of nearby particles of both sizes.

\section{Introduction}

A liquid can be crystallized by slowly lowering
its temperature below $T_{m}$, the crystal melting
temperature. On the other hand a rapid temperature quench
to below $T_{g}$, the glass transition temperature, yields
an amorphous solid with interesting non-equilibrium properties
\cite{sillescu99,gotze92,EdigerAngellNagelJPhysChem1996,DawsonCurrOpColloidInterfSci2002,StillingerScience1995,AngellScience1995,angell00}.
More specifically, the glassy dynamics slow down dramatically and do
so on macroscopically large time scales.  This phenomenon is known
as aging and has been seen in a variety of glass-forming materials
\cite{hodge95,vanmegenpre98,CianciWeeks2006,CianciWeeksSSC2006,angell00,nowak98,cipelletti02,cipelletti03,cipelletti03b,cipelletti05,leheny04,bonn07,munch03,schweizer07a}.
However, the microscopic mechanisms by which a system ages are
still unclear \cite{ngai07pm}.

Dense colloidal suspensions are good models to study glassy systems.
These are composed of small micron-sized solid particles in
a liquid.  In these systems, the glass transition is reached
by increasing the volume fraction, and a glass is formed once
the volume fraction is greater than $\phi_{g} \approx$ $0.58$
\cite{PuseyvanMegenNature1986}.
Previous work looked at aging in colloidal
glasses composed of particles of one size
\cite{CourtlandWeeksJPhysCondMat2003,CianciWeeks2006}.
A small ($5\%$) size polydispersity was used to
inhibit crystallization on experimental time scales
\cite{KegelLangmuir2000,citeulike:2423268}.  In this manuscript, we
study aging in a binary colloidal system.  To an extent the binary
system facilitates comparison with prior simulations
\cite{KobBarratPRL1997,KobTartagliaJPhysCondMat2000}, but our
larger size ratio ($2.1 : 1$) also highlights the role of the
different sizes in the aging process.

In particular, we find that smaller particles are more mobile
than expected. These significantly larger movements are
key components of the aging process since, in our system,
the structural arrangement of the colloidal particles
completely determines the ``age'' \cite{CianciWeeks2006}.
Particles moving to new positions change the structure and
allow the system to become ``older''.  Thus, we conjecture
that the small, unusually mobile particles are a key component
of the aging process, and perhaps facilitate other, more
subtle structural changes of the slower moving large particles
\cite{hoffman92,mewis94,vanmegen01}.  This is reflected in the
tendency for the neighbors of small particles to have higher
mobility on average, and gives insight into the structural
origin of the aging process.

\section{Experimental Methods}
\label{methods}

\begin{figure}[b]
\centering
\includegraphics[width=7cm]{picforpaper.ps}
\caption{Confocal micrograph of the binary colloidal glass studied
in this experiment. The image was taken $20 \mu m$ inside the
sample.  The scale bar represents $5 \mu m$.
}
\label{pic}
\end{figure}

We use two sizes of poly(methyl methacrylate)
(PMMA) particles that are dyed with Rhodamine 6G
\cite{DinsmoreWeitzApplOptics2001,AntlCollSurf1986}, as shown
in Fig.~\ref{pic}.  These particles are suspended in a mixture
of $85\%$ cyclohexylbromide and $15\%$ decalin by weight.
This mixture closely matches both the density and refractive
index of the particles \cite{DinsmoreWeitzApplOptics2001}.
The small particles have a radius of $a_{S} = 0.56 \pm .05 \mu m$
and the large particles have a radius of $a_{L} = 1.17 \pm .05 \mu
m$; the error bars reflect the uncertainty of the mean radius
of each particle species, and additionally each species has a
polydispersity of $\sim 5$\%.  The number ratio of small to large
particles is $1.48 : 1$, but given the large size ratio $1 : 2.1$,
the volume fraction ratio is $1 : 6$.  The total volume fraction
is $\phi_{S+L}\approx 0.62$.  We do not see any demixing of the
two colloidal species, perhaps because the sample is glassy and
thus particle motion is difficult.

In individual {\it monodisperse} samples, the two particle
species would each be expected to have a glass transition at
$\phi_g \approx 0.58$.  For volume fractions slightly less
than $\phi_g$, a monodisperse sample should form crystals in
a finite time, and for samples with volume fractions slightly
above $\phi_g$, previous colloidal samples similar to ours
only find crystals which nucleate at flat boundaries and
which only grow into the bulk of the sample a short distance
\cite{WeeksWeitzScience2000,PuseyvanMegenNature1986,henderson96}.
It is the lack of homogeneous crystal nucleation that is used to
define $\phi_g$ for these monodisperse samples, although in practice
it is probable that a slight polydispersity of the particles is
needed to truly prevent crystallization \cite{henderson96,auer01b}.
Given that the two monodisperse colloidal samples should have the
same $\phi_g$, it is worth noting that our binary sample is not
analogous to mixtures of two polymer glasses with different glass
transition temperatures $T_g$, which are known to have interesting
behaviors \cite{Cowie1989,Cowie2005}.  In practice, for our sample,
we cannot use the lack of crystallization to define $\phi_g$ as
our binary sample does not crystallize at any volume fraction.
Instead, we regard our sample as glassy as the dynamics do not
equilibrate, but change with age for as long as we have observed,
as will be discussed in detail in Sec.~\ref{results}.

We use fast laser scanning confocal microscopy which yields clear
images deep inside our dense samples \cite{prasad07a}.  Despite the
high density, the two colloidal species can be easily discerned.  We
acquire three-dimensional scans of our sample yielding a $53 \times
53 \times 15$~$\mu$m$^3$ observation volume.  The acquisition time
for one 3D image ($\sim 8$~s) is several orders of magnitude faster
than the diffusion time for the particles at this volume fraction.
We continuously follow the sample's aging over $1.5$ hours at a rate
of three images per minute.  This allows us to track $\sim 5000$
particles in three dimensions over the course of the experiment
\cite{CrockerGrierJCollInterfSci1996,DinsmoreWeitzApplOptics2001}.

To study aging, we first must initialize aging in the sample.
There are two approaches to initialization in any glassy
system \cite{mckenna03}.  The first is to quench a sample from a
liquid state.  In a conventional glass, this is done by rapidly
lowering the temperature or increasing the pressure.  In a hard
sphere or colloidal glass, this would be done by starting with a
sample with a volume fraction $\phi < \phi_{g}$ and then rapidly
increase the volume fraction to above $\phi_{g}$.  Currently,
there is no experimental technique that allows for this method
in a sample of hard PMMA colloids.  The second method for
initializing aging involves shear rejuvenation:  by applying
a large shear to a system (above any yield stress), in some
cases the aging appears to be reinitialized.  This is the method
used in this study and previous studies on aging with colloids
\cite{CourtlandWeeksJPhysCondMat2003,CianciWeeks2006,cloitre00,bonn02,viasnoff02}.
There are questions as to whether these two methods lead to the
same glass \cite{mckenna03}, but for our current study, we stress
that only the second of these methods is experimentally available.
Experimentally, this leads to a reproducible initial state
\cite{CourtlandWeeksJPhysCondMat2003} and the dynamics appear
similar to a conventional glass aging via a rapid quench from a
liquid state, as will be shown below.

To shear rejuvenate the sample, we begin each experiment by pulling
a small metal wire that is present in the samples with a hand
held magnet, thus shear melting the glass and re-initializing the
aging process.   We begin taking data, and thus define our initial
time $t_{w}=0$, within 20~s after we finish pulling the wire.
In practice, transient flows continue for a few seconds after
cessation of the stirring, and these flows likely still shear the
glass, resulting in some uncertainty of the initial time.  However,
the 20~s delay between ending the stirring and starting the data
acquisition is long enough to allow for these flows to decrease
quite significantly, and we see no signs of shearing even at the
start of data acquisition.  Thus there is at most an uncertainty
of 20~s for our choice of $t_w=0$, and in practice the results
shown in this work are for time scales much longer than 20~s.

\section{Results}
\label{results}

\begin{figure}[b]
\centering
\includegraphics[width=7.0cm]{lynchweeks.fig1d.eps}
\caption{Aging mean squared displacement for large and
  small particles of a binary colloidal glass at
  $\phi\approx0.62$. The five curves represent five different
  ages of the sample.  The geometric mean ages of the curves,
  from top to bottom, are $\bar{t}_{w} = 8.7, 15, 26, 45$, and
  $76$ min; see Table~\protect\ref{table1} for details.
  The mean square displacements are normalized relative to each
  particle size, that is, for the small particles we plot $\langle
  \Delta r^2 \rangle / a^2_S$ and for the large particles we plot
  $\langle \Delta r^2 \rangle / a^2_L$.}
\label{msd}
\end{figure}

Aging in glasses is defined as the dependence of the dynamical
properties of the system on the time elapsed since vitrification. To
quantify these changing dynamics, we calculate the mean
squared displacement (MSD), defining displacements $\Delta \vec
r_i(t_w,\Delta t) = \vec{r}(t_w+\Delta t) - \vec{r}(t_w)$, and the
MSD as $\langle \Delta r^2 \rangle_{i,t_w}$.  The angle brackets
indicate an average over all particles $i$ (of a given particle
size), and over all times $t_w$; the MSD is thus a function of the
lag time $\Delta t$.  In practice, we wish to understand the $t_w$
dependence of the MSD as well.  Thus, while we continuously take
data for 1.5~h, we divide the data into five temporal windows,
and calculate the MSD restricting the average over times $t_w$
within each window.  The windows are formed logarithmically based
on the sample age $t_{w}$ and each temporal window is
characterized by its geometric mean age $\bar{t}_w$;
see Table~\ref{table1} for details.
For each window, the MSD is calculated for lag times $\Delta t$
up to $\sim 0.5 \bar{t}_w$, and it is important to recognize that
for the largest lag times in each window, the data are aging
during the observation window.  Thus, care must be taken in
interpreting the MSD curves.  Nonetheless, this method provides a way to
characterize the dynamics for a given age $\bar{t}_w$.

Figure~\ref{msd} shows the MSDs of the glass over a total period
of 1.5 h for the large and small particles.  At short time
scales, the motion of a particle is expected to be diffusive as
each particle explores its local environment without noticing
its neighbors. However, at these high densities, particles are
crowded, and this diffusive regime can only be observed on time
scales shorter than our image acquisition rate and is therefore
not visible in Fig.~\ref{msd}. Instead, we observe a plateau for
most of our time scales, indicating that a particle's motion
is slower and more inhibited than normal diffusion.  This is
because particles are trapped in cages formed by their neighbors
\cite{weeksweitzprl02}.  At long time scales, the MSDs
show an upturn, with the time scale of this upturn increasing with
the age of the sample.

\begin{table}
\caption{\label{table1}
We divide the data into five temporal windows, each window
beginning at $t_w = t_1$ and ending at $t_w = t_2$, with values
$t_1$ and $t_2$ listed below.  For each window, $t_2 = 1.7 t_1$
and the ``age'' $\bar t_w$ is defined as the geometric mean of $t_1$
and $t_2$.  All times are in minutes.
}
\begin{ruledtabular}
\begin{tabular}{ccc}
$\bar t_w$ & $t_1$ & $t_2$ \\
\hline
8.7  &  6.7  & 11.3\\
15  &  11.3  & 20\\
26  &  20    & 34\\
45  &  34    & 58\\
76  &  58    & 100\\
\end{tabular}
\end{ruledtabular}
\end{table}

Both the large and small particles' dynamics are aging similarly,
showing upturns at similar lag times.  Clearly in our experiments,
the aging of each species is strongly coupled to that of the other.
Were we comparing two monodisperse samples with different particle
sizes $a_S$ and $a_L$, the horizontal axis of Fig.~\ref{msd}
would need to be rescaled by the diffusion time scale for each
particle, $\tau_D = 3 \pi \eta a^3 / k_B T$, where $\eta$ is the
solvent viscosity, $k_B$ is Boltzmann's constant, and $T$ is the
absolute temperature.  $\tau_D$ thus depends sensitively on particle
radius $a$, and the ratio of diffusion time scales for our particle
sizes is $(a_L/a_S)^3 \approx 9$.  However, in our data, clearly
the upturns in the MSD curves occur at quite similar time scales
(Fig.~\ref{msd}), rather than differing by an order of magnitude.
These essentially similar upturns indicate the system has a single
time scale for structural rearrangements of both species, and this
time scale grows as the sample ages.  As noted in
Sec.~\ref{methods}, this is perhaps reasonable for a binary
colloidal mixture but stands in contrast to mixtures of polymers
with differing glass transition temperatures
\cite{Cowie1989,Cowie2005}.  In fact, it suggests that in
our experiment the relevant time scale determining structural
rearrangements is $\bar t_w$ (which is the same for both species)
and not the diffusion time scale $\tau_D$ (which differs by a
factor of 12).

In Fig.~\ref{msd} we nondimensionalized $\langle \Delta
r^2 \rangle$ to better compare the behavior of large and
small colloids.  The height of the plateau of the MSD curve
is related to the cage size \cite{weeksweitzprl02}, which would be
the same nondimensional size for two monodisperse samples.
In our binary sample, we observe the nondimensional cage
size is larger for small particles. For example, at $\bar{t}_w =
76$~min and using a time scale $\Delta t=10$~min, we find
$(\langle \Delta r_L^2 \rangle)^{1/2} = 0.25 \mu m = 0.21 a_L$, and
$(\langle \Delta r_S^2 \rangle)^{1/2} = 0.43 \mu m = 0.76 a_S$.  Thus,
the small particles are more mobile in their cages, both in an
absolute sense, and especially so when taking into account their
smaller size.

\begin{figure}[b]
\centering
\includegraphics[width=7.0cm]{isf.ps}
\caption{Self intermediate scattering function for each
  particle species.  Each group of five curves represent five
  different ages of the sample.  The geometric mean ages of the
  curves, from bottom to top, are $\bar{t}_{w} = 8.7, 15, 26,
  45$, and $76$ min; see Table~\protect\ref{table1} for details.
  Each scattering function is computed for the wave vector $k_{\rm
  max}$ which maximizes the structure factor for that particle
  species.  The specific values are $k_{\rm max}=5.61$~$\mu$m$^{-1}$
  for the small particles and $k_{\rm max}=3.14$~$\mu$m$^{-1}$
  for the large particles.
  }
\label{scattering}
\end{figure}

A further way to characterize the dynamics is to examine the
self intermediate scattering function for the two species.
This is calculated as
\begin{equation}
F_s(\vec k,\Delta t) = \langle \exp(-i \vec k \cdot \Delta \vec
r_i)_{i,t_w}
\rangle
\end{equation}
where the angle brackets again indicate an average over starting
times $t_w$ and over all particles of a given size, similar to
the definition of the mean square displacement.  In practice, we
choose $| \vec k |$ from the maximum of the static structure factor
for each species, and additionally we average over orientations
of $\vec k$.  The results are shown for the five time windows
in Fig.~\ref{scattering}.  $F_s(\Delta t)$ has a downturn for
both species around similar values of $\Delta t$ for each age
$\bar t_w$.  The overall magnitude of $F_s(\Delta t)$ is larger for
the large particles, reflecting that they move less, as shown in
Fig.~\ref{msd}.  Thus the two particle sizes clearly have different
albeit connected dynamics.  As with the MSD curves, note that at
the largest time scales $\Delta t$ shown in Fig.~\ref{scattering},
the sample is aging over that time scale.  Given that the sample
ages before complete decay can be seen, it is difficult to
characterize the decay of $F_s(\Delta t)$, but we note that the
downturn seen at longer lag times is linked to the small
irreversible structural rearrangements of the sample.

\begin{figure}[t]
\centering
\includegraphics[width=9.0cm]{correctcolortime60percent30ppt.ps}
\caption{
(Color online)
A narrow slice of a three-dimensional sample with width
$\Delta z = 4$~$\mu$m at $t_{w} = 20$~min.
Particles with the $30\%$ greatest displacements $\langle \Delta
r \rangle \geq 0.375$~$\mu$m for $\Delta
t = 10$ min are colored with color indicating the direction
of displacement projected onto the $xy$-plane, where $0^\circ$
indicates motion to the right and $90^\circ$ indicates motion
directly upwards.  When viewed in grayscale, lighter shades indicate
motion upwards, and darker shades indicate motion downwards.
}
\label{matlab}
\end{figure}

We wish to investigate the relationship between structure
and dynamics, that is, how a particle's local environment
influences its motion and the motion of its neighbors
\cite{CianciWeeks2006,CourtlandWeeksJPhysCondMat2003}.
Prior work indicated that the particles rearrange in
cooperative groups \cite{CourtlandWeeksJPhysCondMat2003},
similar to what is seen in supercooled liquids
\cite{EdigerAngellNagelJPhysChem1996,WeeksWeitzScience2000,KegelvanScience2000,DonatiGlotzerPRL1998,harrowell04,sillescu02,berthier04,yamamoto98,ediger00}
and granular materials \cite{Mehta2008,GoldmanSwinney2006}.
A starting point to look for cooperativity
is to ask where the mobile particles are located
\cite{DonatiGlotzerPRL1998,berthier04,CourtlandWeeksJPhysCondMat2003};
Fig.~\ref{matlab} highlights the most mobile particles at
$t_w = 20$~min, using $\Delta t=10$~min.
The coloring indicates the direction of
motion.  Several features are seen.  First, mobile particles
are clustered, as was reported in monodisperse samples
\cite{CourtlandWeeksJPhysCondMat2003}.  A mobile particle is
likely to have mobile neighbors.  Second, particles of similar
colors are likely to be neighbors, indicating that groups of
mobile particles are moving in similar directions.  For example,
the magnified region in Fig.~\ref{matlab} shows a group of particles
all moving approximately up and to the left (green and blue
colors).
Third, occasionally
neighboring particles have significantly different directions of
motion, similar to the ``mixing'' particles seen in supercooled
fluids \cite{weeksweitzprl02}.  This tendency is enhanced
for the small species, which can often move in directions different
from their larger neighbors.  These three observations highlight
that indeed the motion of a particle is related to the motion
of its neighbors, as discussed before.  It is these cooperative
motions that are responsible for the upturn in the MSD curves
(Fig.~\ref{msd}).

\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{lynchweeks.fig6b.eps}
\caption{
(Color online)
Trajectories of two neighboring
particles:  a large particle (red/dark gray) and a small particle
(green/light gray).  The particles have large displacements at
$t_w=20$~min, using $\Delta t=10$~min.  The vertical dashed
lines indicate this time range $t_w,t_w+\Delta t$.
}
\label{traj}
\end{figure}

The trajectories of two mobile particles are shown in
Fig.~\ref{traj}, further demonstrating the cooperative motion.
These are neighboring particles taken from those shown in
Fig.~\ref{matlab}.  One is a large particle (red/dark gray) and
the other is small (green/light gray).  At $t_w=20$~min they both
begin to move in similar directions, with the motion essentially
complete by $t_w=30$~min.  Subsequent to the rearrangement, these
two particles remain caged for the next 30 minutes (further data
not shown).  Figure \ref{traj} shows that these large displacements
happen slowly in our aging glass; while the particles each move
$\sim 1$~$\mu$m, this takes the whole time interval $\Delta
t=10$~min.  Occasionally some particles move more rapidly, and
some move more slowly; similar results to Fig.~\ref{matlab} can
be seen using other choices of $\Delta t$, and indeed motions of
other durations can be seen in the trajectories of Fig.~\ref{traj}
at other times.

\begin{figure}[t]
\centering
\includegraphics[width=8.0cm]{lynchweeks.fig3d.eps}
\caption{Probability of a large particle (solid line) and a small
  particle (dashed line) having a certain number of (a) large
  neighbors and a certain number of (b) small neighbors.  These
  distributions do not change as the sample ages.}
\label{prob}
\end{figure}

\begin{figure}[b]
\centering
\includegraphics[width=7.0cm]{lynchweeks.fig5b.eps}
\caption{Pair correlation functions
  $g_{SS}(r)$ (solid line), $g_{LL}(r)$ (dotted line), and $g_{SL}(r)$
  (dashed line).  The arrow at $r=2.8$~$\mu$m indicates the cutoff used to define nearest neighbors.
  These functions do not change as the sample ages.}
\label{gr}
\end{figure}

These observations do not yet show
how the local structural environment influences the motion.
To quantify the structure of the local environment of a
particle, in Fig.~\ref{prob} we plot the probability of a
particle having $N_{L}$ large neighbors (a) or $N_{S}$ small
neighbors (b). The probabilities are computed separately for
large and small reference particles. We define a neighbor as a
particle within 2.8~$\mu$m from the reference particle being
considered.  This cut-off distance is set by the first minimum
of the pair correlation function for large particles, as shown
in Fig.~\ref{gr}.  However, the trends seen in the distributions
of Fig.~\ref{prob} are not sensitive to small variations
of this value.  Consider Fig.~\ref{prob}(a) which shows the
probability of having $N_{L}$ neighbors for both species of
particles. The two distributions are very similar
indicating that the neighborhoods of both large and
small reference particles contain similar numbers of large
particles. Figure~\ref{prob}(b) shows a similar plot counting
the likelihood of having $N_{S}$ small neighbors for both
species. Again we find comparable distributions, showing that
the sample is statistically homogeneous and there is no size
segregation.  The most probable neighborhood [the maximum of
the probability distribution $P(N_L,N_S)$] is $N_L=7, N_S=10$
for large particles, and $N_L=6, N_S=10$ for small particles.

\begin{figure}[b]
\centering
\includegraphics[width=8.0cm]{lynchweeks.fig4d.eps}
\caption{Average particle mobility versus the number and type of
  neighbors the particle has (circles) for $\Delta t =$ 10 min.
  The solid line is proportional to the probability of having
  that number and type of neighbors.  (a) Large particles with
  large neighbors, (b) large particles with small neighbors, (c)
  small particles with large neighbors, (d) and small particles
  with small neighbors.  The trends in these graphs do not change
  as the sample ages or with choice of $\Delta t$.  The curves are
  truncated where the probabilities of having that many neighbors
  drops below 3\%.}
\label{mob}
\end{figure}

The probability distributions in Fig.~\ref{prob} do not depend
on sample age.  This is in agreement with previous work on
aging monodisperse glasses where other local measures of
structure did not show any aging \cite{CianciWeeks2006}.  However,
in that study the structure showed some correlation to the
dynamics. To establish whether this still holds in a binary
suspension, we first calculate the mobility of each particle.
To do this, we choose a fixed lag time $\Delta t = 10$~min,
and calculate $\Delta r^2(t_w,\Delta t)$ for each particle
at each time $t_w$.
The results that follow do not depend on this choice of $\Delta t$.
We then average the mobility over all
particles of a given size with the same number of large or
small neighbors, and plot these results in Fig.~\ref{mob}.
We observe that, on average, having more small neighbors
allows a reference particle to be more mobile.  Conversely,
a local environment rich in large neighbors tends to inhibit
mobility. This seems to hold whether the reference particle is
large or small, as can be seen by comparing the top and bottom
rows of Fig.~\ref{mob}.  Therefore, in this binary sample,
not only are the small particles more mobile on average,
but they also facilitate the mobility of their neighbors
independently of their size.  Of course, this is only true on
average; these results vary quite a lot from particle to particle
\cite{conrad05}.

This is similar to previous observations of binary systems
using rheology \cite{hoffman92,mewis94} and light scattering
\cite{vanmegen01}.  These studies found that binary suspensions
have lower viscosities and faster microscopic motion compared
to monodisperse samples of equal total volume fraction.  This was
attributed to the fact that binary suspensions are capable of
packing to higher volume fractions; thus at a given volume
fraction, there is more free volume for a binary suspension as
compared to a monodisperse suspension.  In particular, a glassy
binary suspension has faster dynamics and can even become
liquid-like with a sufficiently high population of small
particles \cite{vanmegen01}.  For our experiment, we have
purposely chosen the overall volume fraction to ensure that the
sample still has glassy behavior.  Nonetheless, Fig.~\ref{mob} shows
that the
small particles ``lubricate'' the motion of the large ones, as
proposed to interpret the earlier experiments \cite{vanmegen01}.

\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{lynchweeks.fig2e.eps}
\caption{$\langle{\Delta x^{2}} \rangle$ (solid) and $\gamma =
  \sqrt{\langle x^{4} \rangle / 3} $ (dotted) as a function of sample
  age for large particles and small particles as
  indicated.  The displacements are defined using $\Delta
  t = 10$ min.  The solid lines are fit to power laws (dashed).
  Inset:  Same data plotted on a log-log plot; the upper line
  corresponds to the small particles, and the lower line to
  the large particles.
}
\label{x2}
\end{figure}

To further investigate the role the small particles have in the aging
of the sample, we examine the distribution of particle
displacements over the same lag time of $10$ minutes.  Figure~\ref{x2}
shows $\langle \Delta x^2(t_w) \rangle$ (solid lines), which is the second
moment of the distribution of displacements, as a function of sample
age, $t_{w}$, for large particles (top) and small particles (bottom).
Note that here the angle brackets $\langle \cdot \rangle$ indicate an
average over all particles, but not over times $t_w$.  The data
decrease, reflecting the slowing of motion already seen in
Fig.~\ref{msd}.  The small fluctuations in each graph indicate
periods of extra activity (of the sort pictured in
Fig.~\ref{matlab}), and have been seen previously
\cite{CourtlandWeeksJPhysCondMat2003,cipelletti03,cipelletti05,castillo07}.

The non-Gaussianity of the sample can be determined by
comparing these solid curves to the dotted curves, which
represent $\gamma = \sqrt{\langle \Delta x^{4} \rangle / 3}$.
The ratio $\gamma / \langle \Delta x^{2} \rangle = 1$ for a
Gaussian distribution. In our sample, this ratio is $\approx
1.2 $ for the large particles and $\approx 1.6$ for the
small particles.  Both of these ratios are greater than $1$
indicating that the dynamics are anomalous: there is an
excess of large steps when compared to a Gaussian distribution
with the same width. This result is typical in a monodisperse
dense suspension \cite{WeeksWeitzScience2000}. Interestingly,
this ratio is slightly greater for the small particles,
supporting our assertion that the small particles move in
an even more heterogeneous way. This suggests that
the subset of small particles with unusually high mobilities
dominate the dynamics.  Note that we use $\Delta t=10$~min
because the experiment includes time scales $t_w \ll \Delta
t$, $t_w \sim \Delta t$, and $t_w \gg \Delta t$, and the
non-Gaussian behavior of the small particles is present at all of
these time scales.

To characterize the slowing of the dynamics, we fit the large $t_w$
decay of $\langle x^2 \rangle(t_w)$ in Fig.~\ref{x2} to a power
law decay of the form $\langle \Delta x^{2} \rangle \sim t^{-b}$
where $b_{L} =0.64 \pm 0.05$ for the large particles and $b_{S}
= 0.60 \pm 0.05$ for the small particles; these fits are shown
as dashed lines in the figure.  The inset of Fig.~\ref{x2} shows
a log-log plot with the same data and highlights the power-law
decay at large $t_w$.  We note that it is possible and perhaps
likely that this power-law decay is a transient effect that could
disappear at substantially longer $t_w$ \cite{mckenna03}; for our
data, we regard this as simply one way to characterize the behavior.
The values of the exponents $b$ vary from experiment to experiment
and are in the range of 0.60 - 0.80 for this sample.  Within our
uncertainty, the large and small particles decay almost at the
same rates, as suggested by the similar shapes of the MSD curves
(Fig.~\ref{msd}).  However, comparisons between different data
sets always find $b_L > b_S$ by a small amount.  A study of aging
in a monodisperse colloidal sample found a range of power law decay
exponents of $0.05 < b < 0.5$ \cite{CourtlandWeeksJPhysCondMat2003}.
This prior experiment studied colloidal particles with radius $a
\approx 1.18$~$\mu$m which is intermediate to the particle radii
used in our experiment.  However, the power law decay exponent
range found in Ref.~\cite{CourtlandWeeksJPhysCondMat2003} is
less than both exponents found in this study; apparently this
binary suspension ages ``faster''.  Given our observations that
the smaller particles facilitate the motion of their neighbors,
perhaps the presence of these small particles results in this faster
aging as compared with Ref.~\cite{CourtlandWeeksJPhysCondMat2003}.

\section{Conclusion}

We study aging in a binary colloidal glass using confocal
microscopy and distinguish between the small and large
particles within the sample.  These two particle species
slow down with age similarly to each other, but the smaller
particles appear to be more important in the aging process.
These small particles are more mobile, more dynamically
heterogeneous, and facilitate the motion of other small
particles and the motion of the large particles.  Furthermore,
significant motions within the sample (presumably responsible
for the aging) take place in cooperative groups of mobile
particles.  From a practical viewpoint, these results suggest
that the aging of a sample could be influenced by controlling
the size ratio and number ratio of the two species.  Another
likelihood is that a poorly-mixed sample would have spatially
heterogeneous aging, depending on the local composition.

\acknowledgements

We thank E.~Baker, G.~L.~Hunter, and C.~B.~Roth for useful
conversations.  This work was funded by the National Science
Foundation under Grant No.~DMR-0239109.  The work of J.~M.~L.~was
also supported by a SIRE Independent Research Grant from Emory
University.

Soft Condensed Matter (cond-mat.soft); Disordered Systems and Neural Networks (cond-mat.dis-nn)


Magnetic Fields in Clusters: Theory vs. Observations

Abstract: It is now well established that the plasma in galaxy clusters is magnetized. In some cases, at least, the field is strong enough to be dynamically important. Perhaps this is the time to move past simple detection experiments, and to work towards a general understanding of the strength, structure and maintenance of the cluster field.

\section{Introduction}

I am not aware of any astrophysical plasma in which the magnetic field
is clearly unimportant.  Space plasmas, which we can observe {\it in
situ}, are commonly magnetized; in some cases the magnetic field
and associated currents dominate the plasma dynamics.  The
interstellar medium (ISM) in our galaxy is magnetized, at
approximately an equipartion level.  The diffuse plasma in radio
galaxies is obviously magnetized, probably at dynamically important
levels.

Until recently, magnetic fields have generally been ignored in
dynamical models of the intracluster plasma medium (ICM).  However,
this is changing as data accumulates -- and the discussion at this
conference is a case in point.  It seems very likely that the ICM is
magnetized in every cluster -- just as we believe the interstellar
medium (ISM) is
magnetized in every galaxy.  I suspect the dynamics of the magnetized
ICM are in many ways similar to the dynamics of the magnetized ISM.
The magnetized plasmas in both systems are probably subject to strong
effects of  cluster, or galactic, weather ({\it e.g.}, Binney 1998,
Zweibel \& Heiles 1997).  That is, seed fields are probably provided
to the ISM by stellar ejecta; and to the ICM by ejecta from both
active and ``normal'' galaxies (Kronberg {\it et al} 1999, V\"olk \&
Aronyan 1999).  The seed fields are enhanced and shaped by turbulence,
flows and dissipation, a complex process which results in the fields
we see at present.

Up to now observational work on cluster fields has mainly been aimed
at detection:  can we find  any evidence of magnetic fields in any
clusters?  The answer to that seems clearly to be ``yes'', as this
meeting has demonstrated.  It may be time to change our focus.  We
know that the ICM is magnetized;  we should concentrate now
exploring the cluster field in more detail, both observationally and
theoretically.  In this paper I will try to convey a rather personal
impression of what we know about magnetized plasmas from other sources
(laboratory and space data, simulations, theory), and to consider how
this information can be used to understand magnetic fields in
clusters.

\section{What Does the Data Tell Us?}

We have two ways to detect magnetic fields in the cluster plasma.
We can detect diffuse synchrotron emission (from radio haloes and
relics), or we can detect Faraday rotation from individual radio
galaxies, either behind the cluster or embedded in it.

\subsection{Overview of the Observations}

We have three possible observational approaches to measuring the
magnetic field in a cluster of galaxies.

{\bf Radio Haloes and Relics}.  Some clusters show diffuse synchrotron
emission which is associated with the ICM rather than with discrete
radio galaxies within the  clusters ({\it e.g.} Feretti \& Giovannini
1996, Feretti 1999).   This emission can exist
throughout the cluster (haloes), or be peripheral (relics).  In either
case, this is direct evidence of  intracluster B fields.
Minimum-energy or minimum-pressure arguments, sometimes combined with
X-ray inverse Compton detections or limits, suggest $B \>\hbox{\lower.25em\hbox{$\buildrel <\over\sim$}}\> 1 \mu$G is
typical on large scales (comparable to the size of the cluster core,
or the cluster itself).

{\bf Faraday Rotation of Background Sources}.  Comparison of Faraday
signals from cluster and non-cluster sources allows a statistical test
of cluster magnetic fields.  Current studies do show an excess rotation
measure from sources in or behind rich clusters (Kim {\it et al} 1991,
Clarke 1999).  Modelling this data in terms of a ``typical'' ICM
density and scale, and a cell size $\sim 10$ kpc, again suggests a
microgauss-level field is typical of the cluster on moderate to large
scales.

{\bf Faraday Rotation of Embedded Sources}.  Another approach
uses high-resolution images of radio sources within the cluster.  As
this is the closest to my own research, I will say a bit more about
it. Radio imaging gives us our only  evidence on the structure of the
magnetic field.  Most studies have been of cluster-center sources, but a
few studies of non-central sources now exist ({\it e.g.} Feretti {\it
et al} 1999).  In Eilek \& Owen (1999), we present new Faraday data on
two large, tailed radio sources (3C75, 3C465) in the centers of two
clusters (A400, A2634 respectively) which have only
weak cooling cores.  We also collect published results of other imaged
Faraday data on cluster-center sources (which have studied smaller
radio sources in strong cooling cores).
These results are summarized in Table 1, which lists
the radio source, its radial extent, the
characteristic scale of its measured rotation measure, the mass
accretion rate inferred from simple cooling-flow analysis, the mean
magnetic field along the line of sight inferred from Faraday data, and
the ratio of magnetic pressure to gas pressure.  This data allows us
to make a general statement:  the magnetic field is dynamically
significant in the central regions of these clusters.  The typical central
field strength $\sim 10-30 \mu$G.  Comparison to pressures derived
from X-ray data shows that $B^2 / 8 \pi \sim n k_B T$ is typical for
the cluster-center sample.  We can also use the images to estimate
coherence lengths for the rotation measure.  The general trend we see,
is that cthe luster-core rotation measure has a characteristic length
scale $\sim$ 10 kpc.  This is our only direct information so far on
the coherence length of the cluster field.

\begin{table}
\begin{tabular}{lccccc}
\multicolumn{6}{c}{{\bf Table 1}. Cluster-Center Magnetic
Fields$^{\rm  a}$}\\
\hline\hline
Source &$ L_{RS}$ & $L_{RM}$ & $\dot M$ & $\langle B_{\parallel} \rangle$
     & $ p_B / p_g$ \\
  & kpc  & kpc & $M_{\sun}$/yr & $\mu$G
\\
\hline
M87              &  3  & 1  & 10  & 50 & 0.7
\\
A2199            & 4  & 3 & 150 & 60 & 5.0
\\
A1795            & 7 & 3  & 400 & 10 & .09
\\
A2052            & 8 & 2  & 90  & 40 & 4.5
\\
A2029            &  10 & 3 & 370 & 75 & 13
\\
A4059            &  10 & 3 &  120 & 25 & 1.8
\\
Hydra A          &  50 & 5 & 600 & 60 & 13
\\
Cyg A            & 70 & 10  & 200 & 40 & 4.5
\\
A 400            &  100 & 10 & $< 1$ & 7 & 0.7
\\
A 2634           &   140& 20  & $< 1$ & 5 & 0.4
\\
 \hline\hline
\end{tabular}
\begin{list}{}{}
\item[$^{\rm a}$]From Eilek \& Owen 1999, and references therein.
\end{list}
 \end{table}

\subsection{What is the Nature of the Cluster Field?}

Collecting all the evidence, I think there is a strong case for
magnetic fields being common in all clusters at (at least) modest
levels.  Both the radio haloes and the Faraday rotation of background
sources are more or less consistent with a microgauss level field
existing throughout most clusters.\footnote{Why don't we detect radio
haloes from all clusters? I suspect it is a question of sensitivity.
The ICM is very likely magnetized in nearly all
clusters. The galaxy ejecta that carry seed fields should also carry
``seed'' cosmic rays. Could it be the case, then, that we are seeing
only the brightest few haloes with our present observations?}  A
field at this level is not dynamically important at typical ICM densities
and temperatures.  We also know, however, that stronger fields (tens
of $\mu$G) exist in cluster centers, where we detect them in front of
central radio galaxies.  These fields are dynamically important.
There are also hints of stronger fields associated with non-central
radio galaxies.  We do not yet know whether these stronger fields are
typical of the diffuse ICM in all clusters, or whether they are
amplified by energy input from the radio galaxies.

What is the next step?
We would like to use these observations to answer fundamental
questions.  What is the strength and structure of the field?  What is
its role in, and relation to, the dynamics of the cluster gas?  The
observations we can make are both tantalizing and frustrating.  They
cannot answer these qustions directly, but can provide important
clues.  We should look to other information that we have on magnetized plasmas.

\section{There Are No Uniform Magnetic Fields}

One is often tempted to interpret an observation in terms of the
simplest possible model.  For the case of magnetic fields, the
very simplest model is a uniform field throught the cluster; the next
simplest is  the old idea of random, Gaussian ``cells'' of field.
We should be very cautious here.  These seems to be no plasma
(whether astrophysical or simulated) in which the field is truly
uniform, nor in which the field structure can be described by simple,
space-filling cells.  The magnetic field in a magnetized plasma is
more likely to be bunched into elongated, high-field regions:  we can
call  these ``filaments'' or ``flux ropes''.

\subsection{Evidence for Filamented Fields}

The evidence comes from several areas.

{\bf Simulations of MHD Turbulence} find that intermittency is common:
the magnetic field is organized into localized, often elongated,
high-field regions {\it e.g.}, Menguzzi {\it et al} 1981, Nordlund
{\it et al} 1992, Miller {\it etal} 1996, Kinney {\it et al} 1995).
This appears similar to the small-scale vortex intermittency seen in
hydrodynamical turbulence;  vorticity in fluid flows sometimes behaves
similarly  to magnetic fields in MHD flows.   The characteristic
intermittency scale is  thought to be similar to the turbulent
dissipation scale (but this conculsion may still be sensitive to
numerical resolution).

{\bf Space Plasma Data} is measured {\it in situ}.  These measurements
again find direct evidence that the field is often bunched into flux
ropes. These are apparent in the flux transfer events that are part of
patchy reconnection at planetary magnetopauses, and they have been
detected in interplanetary space, associated with coronal mass
ejections. (The compilation edited by Russell {\it et al.\/}\ 1990 is a good
source of references to this area.)  In addition, optical and X-ray
images of the sun make it clear that solar plasma is filamented;  this
is almost certainly due to magnetic fields.\footnote{In many ways the
sun is a poor analogy to cluster magnetic fields, due to the very
heavy ball of gas in which the field lines are ``footed''.  However I
suspect that the field dynamics on small scales are probably similar,
although one must take into consideration the more ordered turbulent,
convective flow that exists on the sun.}

{\bf The Magnetized ISM} displays rich substructure and
filamentation. Dramatic examples include the filamentation apparent in
HI images (Hartmann \& Burton 1997), or recent radio images of the
galactic center (Lang {\it et al} 1999).  Once again, we do not have
direct evidence that the filaments are magnetically created, but such
seems very likely to be the case.

{\bf Extragalactic Radio Sources} commonly show filamentation in their
synchrotron emission.  A particularly beautiful example is the new
image of the M87 halo (Owen, Eilek \& Kassim 1999, also Owen 1999).
These are almost certainly high-field regions; polarization data when
available shows that such filaments are aligned with the magnetic
field.

\subsection{Consequences of Filamented Fields}

Inhomogeneous fields affect ``simple'' interpretation of observations
in several ways; we must be cautions.

 {\bf Minimum pressure or equipartion} estimates are well known
to be sensitive to the inhomogeneity of the source.  In following the
usual applications (for instance, that of Pacholczyk 1970), we must
remember that two filling factors are needed.  One is the fraction of
the volume filled by relativistic particles, which may of course be
less than unity.  The second is the fraction of the volume filled by
relativistic particles that is {\it also} filled with magnetic field
(denoted by $\phi$ in the standard notation).  Incomplete volume
filling leads to a {\it higher} magnetic field value, derived from
these standard arguments, than would be the case for a homogeneous
source.

 Interpretation of {\bf inverse Compton} (IC) X-ray detections or
limits is also sensitive to the degree of mixing of the relativistic
electrons and the magnetic field.  Incomplete mixing (such as magnetic
filamentation) will lead to a lower ratio of synchrotron to IC
emission than would apply to a homogeneous source. Conversely, an
observation of a fixed IC/synchrotron ratio requires a higher magnetic
field value if the source is incompletely filled.

{\bf Particle aging} due to synchrotron losses is also sensitive
to magnetic field inhomogeneity.  Particles lose energy more rapidly,
and radiate at higher frequencies, in a strong field.  One consequence
of this is that the high-frequency decay of the synchrotron spectrum
can be dramatically slowed down if the particles spend much of the
time in weak-field regions (Eilek, Melrose \& Walker 1997).  Applied
to interpretation of data, one must be aware that radio sources can be
much older than they look.

Finally, interpretation of {\bf synchrotron spectra} must be
done with care if the magnetic field is inhomogeneous.  The usual
analysis assumes a uniform field, and interprets the photon spectrum
as a direct measure of the particle energy distribution.  In reality,
however, the photon spectrum is a convolution of the field and
particle distributions.  One example is that a power-law distribution
of magnetic field strength, combined with a peaked particle
energy distribution, can produce a power-law photon spectrum
(Eilek \& Arendt 1996).  It is quite likely that observed spectra are
telling us as much about the field distribution as about the
particles.

\subsection{Nature of the Filaments}

What can we say about the small-scale structure of these magnetic
filaments?  I think our best evidence comes from space plasmas, where
{\it in situ} field measurements are possible.  The filaments tend to
be structured as a {\it flux rope}.  That is, a twisted (helical)
magnetic field structure, in which the field lies along the axis in
the center of the tube, and becomes increasingly helical going away
from the axis.   Such flux ropes can be confined by external pressure,
or can be partly or fully self-confining, depending on the exact
structure of the field.  The latter is particularly nice, as it allows
strong fields to be self-organized into filamentary structures without
the need for external pressure confinement.

What is the origin of these structures?  Fluid flow offers analogs.
Consider is a tornado:  a strong vortex line which
forms due to the combination of shear and convergeance in the local
atmosphere. As vorticity and magnetic field are often useful analogs of
each other,  one would expect magnetic flux ropes to arise from
similarly favorable flow fields.  (One MHD example might be the solar
convective
zone, where converging flows accumulate the magnetic field into
strong, localized, low-density filaments.   Add a bit of
circulation and
one can imagine a twisted flux rope arising.) One can also look to
space plasma analogs.  Flux rope formation associated with flux
transfer events is thought to be due to plasma instabilities,
particularly reconnection, again in the presence of some twist.
Tearing instabilities lead to magnetic island
formation in two-dimensional problems (such as the neutral sheet in a
reconnection situation), and thus connect to helical flux ropes in three
dimensions.

In other domains, MHD simulations  find that passive magnetic fields
(weak compared to flow energy density) in sheared flows collect into
filaments (such as the radio galaxy simulations of  Clarke
1996).   Laboratory evidence shows that strong fields (compared to
plasma energy densities) can organize themselves into force-free
structures ({\it e.g.} Taylor 1986); in cylindrical geometry these have
the general flux-rope structure described above.  Thus, it seems that
we can generally think of magnetic fields as filamented, probably on a
ranges of scales, whether the field is weak or strong.

\section{Large Scale Structure of the Magnetic Field}

The existence of small-scale filaments does not preclude an organized
large-scale structure.  An obvious example is the sun:  the surface
field is strongly inhomogeneous and filamented, and yet shows a simple
dipolar structure when measured on larger scales.  What can we say
about the large-scale structure of the cluster field?

\subsection{Observational Evidence}

We have little direct information as yet.  We know that the magnitude
of the rotation measure seems to drop going away from cluster center;
this is apparent in the statistical data of Clarke (1999), and also in
the two large-scale sources measured by Eilek \& Owen (1999).  We also
know that large-scale radio haloes have a finite size, which is no
larger than the cluster size.

Both of these observations suggest that the field strength decreases
with distance from the cluster center (except possibly clusters with
in offset radio relics, which may have an enhanced magnetic field near
their periphery).  However, both observations must be interpreted with
care.  The observed decay of rotation measure must be due in part to
the decrease of gas density, as well as to the field structure.
Similarly, the structure of a radio halo is determined by the
relativistic particle density as well as the field strength.
I am not aware of any work at this point which can separate out the
effects of the field and the particles in these observations.

\subsection{Theoretical Arguments}

We can make one general statement:  cluster fields must involve a
dynamo.  For our purposes, I define a dynamo (interpreted
loosely) as anything that plasma flows can do to a magnetic field.
The most common application involves the inductive enhancement of a
weak seed field to interesting levels.
(This comes from the $\bf{v} \times \bf{B}$ term in the induction
equation, given for instance in Garasi \& Eilek 1999)  As a general rule,
flows can amplify a field up to a state of dynamic balance, $B^2 \>\hbox{\lower.25em\hbox{$\buildrel <\over\sim$}}\>
4 \pi \rho v^2$.  However, the
rate of growth and the final field level depend on details of the flow
field, particularly turbulent and dissipative effects. In addition,
some situations (such as flux ropes in the sun) result in the field on
small scales coming into balance with the static pressure:  $B^2
/ 8 \pi \sim n k T$.   This results in a higher field that that
produced by the usual dynamic balance estimate, if the flows
controlling the dynamo are subsonic. Figure 1  illustrates some of the
possible branch points in turbulent and/or advective dynamos.

\begin{figure}
 \resizebox{\hsize}{!}{\includegraphics{Eilek1_fig1.ps}}
\caption[]{Schematic flow chart for a cluster dynamo.  The key point
is that plasma flows can support both large-scale and small-scale
magnetic fields;  the nature of the resultant field depends, of
course, on the detailed properties of the flow.  Beck {\it et al.\/}\ (1996)
has a similar figure and a detailed discussion of galactic dynamos.}
\end{figure}

What can dynamo theory tell us about the cluster field?  We can, in
principle, inquire about the strength and structure of the field.

\subsubsection{Field Strength}

Several authors have suggested that turbulent dynamos maintain a low
level magnetic field throughout the cluster ({\it e.g.} Ruzmaikin {\it
et al} 1989, Goldman \& Rephaeli 1991, De Young 1992). These models
generally address only the local field strength (from the energy
balance relation given above), and assume the field is disordered with
small-scale tangling.  The field strength depends directly on the
turbulent velocity assumed:  $v_t = 400$km/s gives $B = 6 \mu$G if the
plasma density is $10^{-3}$cm$^{-3}$.

An important related question is what drives the turbulence that these
models require to maintain the field.  Earlier work argued that galaxy
driving in a relaxed cluster could maintain interesting levels of
turbulence, and magnetic field, throughout
the cluster. However, more recent work disagrees.  There now seems to
be a general agreement that galaxy driving cannot maintain the strong
turbulence, (probably transonic;  {\it e.g.} Baum 1992, or Keel {\it
et al} 1996) and strong magnetic fields, seen in cluster cores.  It
may not even be able to maintain a microgauss level field.  I have done
calculations for turbulent driving in a rich, Coma-type cluster, and
find turbulent velocities of a few tens of km/s, which maintain fields
at only a fraction of a microgauss.

This difficulty leaves one with two options.  One can argue that the
stronger fields discussed at this meeting are atypical, perhaps
maintained only locally by energy input from an embedded radio
galaxy.  Alternatively, ongoing evolution of the cluster gravitational
potential (otherwise known as cluster mergers) can provide a stronger
energy input, in principle maintaining stronger turbulent and field
levels.  Norman \& Bryan (1998) indeed found turbulent velocities
between 1/4 and 1/2 of the virial velocity were maintained for
significant times in one set of merger simulations.

\subsubsection{Field Structure}

Less work has been done on the spatial structure of the magnetic
field.

Several authors have noted that the compression associated with
convergent cooling flows will amplify an initial seed field towards
the cluster center.  For a field which remains purely turbulent and
simply advected, one would expect $B \propto \rho^{3/2}$ (this would
be the case if fluid turbulence or reconnection maintained a locally
turbulent field during the inflow).  An alternative picture is that of
Soker \& Sarazin (1990), who proposed that the inflow stretches an
initally disordered field radially.  They find the field varies as $B
\propto 1 / r^2$ in this case.  Such arguments also assume a weak,
passive field; when the field becomes strong it will affect the flow
and break these scaling laws.  Numerical simulations (such as Dolag,
Bartlemann \& Lesch 1999, Garasi \& Eilek 1999, Roettinger 1999) will
be required to determine the final answer here.

What field structure td simple turbulent dynamo models predict?
One expects that the turbulently maintained field will be strong where
the turbulence in strong.  Locally strong turbulence might be the
reason why some haloes, or relics, do not share the symmetry of the
cluster in which they sit.  I have not seen this addressed in
published work; again, ongoing numerical simulations have the
potential to address this.

The behavior of turbulent dynamos has the potential to be much richer
than these simple arguments.  In particular, laboratory plasmas show a
capacity for {\it self-organization}.  If the system has the proper
initial asymmetry (helicity), the magnetic field can undergo a
 ``relaxation'' process to a final state which displays more order on
large scales than the initial state had (Taylor, 1986).  I have
considered one possible final state for clusters, namely, nested
toroidal flux
tubes.  This solution predicts that the magnetic field increases
outward from a central minimum, and then decays as $B \propto 1/r$.
It is not yet clear whether these formal solutions (derived from
laboratory considerations) apply to clusters which are still in a
non-steady state.  We must keep in mind, however, that dynamo effects
in clusters may be more complex than has been considered so far.

\section{Impact of the Field on the ICM}

To summarize:  the  observations suggest that microgauss fields
 extend throughout the cluster plasma.  Larger fields exist around
embedded radio sources;  it is not yet clear if they are
representative of the general diffuse ICM.  What are the consequences
of a such magnetic fields for the larger questions about the clusters?

\subsection{Microphysics}

There has been quite a lot of discussion in the literature regarding
the effect of magnetic fields on transport processes in the ICM.  This
has come up most particularly in regard to thermal conduction (or its
suppression) in cooling-flow models;  and it is also relevant to
questions of relativistic particle diffusion in radio haloes.  The
detailed workings of particle transport depend on the detailed
structure of the magnetic field, particularly on small scales.  The
literature tends to split into two camps:  those working with a
spectral representation of the turbulent magnetic field ({\it e.g.}, Jokipii
1966),  and those working with a chaotic representation ({\it e.g.}, Chandran
\& Cowley 1998).   The essence of both approaches involves determining
a characteristic, coherence, length for the turbulent field (if,
indeed, one exists).  It may be here that observations (such as the
structure of Faraday rotation images) can make connection with the
theory.

\subsection{Structure}

A significant magnetic field in the center of the cluster will of
course affect the hydrostatic balance, and thus the structure, of the
ICM.  In a (quasi) static
situation the field can provide an important part of the pressure
support.  This may be part of the answer to the discrepancy found
between cluster masses determined from gravitational lensing, and from
X-ray imaging (Miralda-Escud\'e \& Babul, 1995;  Loeb \& Mao 1994).
Just how general this result is does not seem clear yet, as only a few
clusters have been so analyzed both in X-rays and with lensing. It is also
worth noting that X-ray mass determinations will also be misleading if
the ICM is itself dynamic rather than in hydrostatic equilibrium ({\it e.g.},
Squires {\it et al.\/}, 1996).   The theoretical approach to this question will
probably lean heavily on numerical simulations of cluster evolution
({\it e.g.}, Garasi \& Eilek 1999, Roettiger 1999).

\subsection{Dynamics}

A significant central magnetic field will also impact the dynamics of
the ICM.  In general, a magnetic field providing pressure support to
gas at the bottom of a potential well is unstable; it is energetically
favorable for the dense gas above to fall, and the magnetized, lower
density gas to rise.  On small scales, one would expect turbulence or
convection to be driven by magnetic bouyancy ({\it cf.}\ Balbus and Soker,
1989, also Zoabi, Soker \& Regev 1996, for applications to clusters).
On larger scales, one would expect a magnetically supported cluster
core to be subject to Parker-type instabilities (Parker, 1966).  In
either regime, the end effect should be the same -- we should expect
gas in the central region of a magnetized cluster to be disordered,
turbulent, and anything but static.

\begin{acknowledgements}
I am grateful to many of my colleagues for
conversations and collaborations relating to magnetic fields in
clusters.  It is a particular pleasure to thank  Frazer Owen, Chris
Garasi, Ellen Zweibel, Greg Taylor, and
Tracy Clarke for recent discussions. However, any errors or
misconceptions herein are my fault, not theirs.  This work was
partially supported by NSF grant AST-9720263.
\end{acknowledgements}


Anomalous Pressure Dependence of the Superconducting Transition Temperature in TlNi$_2$Se$_{2-x}$S$_x$

Abstract: We report the pressure dependence of the superconducting transition temperature, $T_c$, in TlNi$_2$Se$_{2-x}$S$_x$ detected via the AC susceptibility method. The pressure-temperature phase diagram constructed for TlNi$_{2}$Se$_{2}$, TlNi$_{2}$S$_{2}$ and TlNi$_{2}$SeS exhibits two unexpected features: (a) a sudden collapse of the superconducting state at moderate pressure for all three compositions and (b) a dome-shaped pressure dependence of $T_c$ for TlNi$_{2}$SeS. These results point to the nontrivial role of S substitution and its subtle interplay with applied pressure, as well as novel superconducting properties of the TlNi$_2$Se$_{2-x}$S$_x$ system.

\section{Introduction}

The phase diagram of iron-based superconductors often feature a superconducting region in the proximity of antiferromagnetism. The multi-orbital nature of the system and the nesting between well-separated electron and hole Fermi surfaces provide a framework for the discussion of spin-fluctuation-mediated superconductivity as well as spin-density-wave (SDW)  type antiferromagnetism \cite{Mazin08, Kuroki08}.  The 122-type Fe-pnictide family of materials offers a prominent example of these effects \cite{Ishida09, Paglione10, Johnston10, Shibauchi14}, where the application of pressure or chemical substitution can tune the system away from antiferromagnetism and towards superconductivity.

The situation in the 122-type nickel-based systems is however rather different. The stoichiometric parent compounds already exhibit superconductivity, with BaNi$_2$P$_2$ \cite{Mine08}, BaNi$_2$As$_2$ \cite{Ronning08}, SrNi$_2$P$_2$ \cite{Ronning09} and SrNi$_2$As$_2$ \cite{Bauer08} having comparatively low $T_c$'s of 3~K, 0.7~K, 1.4~K and 0.6~K respectively. The structural transition seen in the Fe-based compounds is of first order in the Ni materials, and may even be absent altogether \cite{Bauer08,Hirai12}. Importantly, magnetic ordering is universally absent in these superconducting materials even with tuning, for example in chemically tuned BaNi$_2$(As$_{1-x}$P$_x$)$_2$ (0 $\leq$ x $\leq$ 0.13) \cite{Kudo12} and pressure tuned BaNi$_2$As$_2$ (up to 27~kbar) \cite{Park10}.

These observations, taken together with studies of electron-phonon coupling \cite{Subedi08} and Fermi surfaces \cite{Shein09,Ideta14,Zhou11,Terashima09} suggest that the superconductivity in the Ni-based systems is likely of the conventional type, although there remain further compounds to be explored. Understanding the phase diagrams and Fermiology of these materials is an important step in understanding the complex physics of the Fe- and Ni-based pnictides and chalcogenides as a whole, and may help shed light on how the Fe-based pnictides achieve such high $T_c$'s.

TlNi$_{2}$Se$_{2}$\ is a relatively new and unstudied member of this family. It superconducts below 3.7~K, and crystallizes in a tetragonal ThCr$_2$Si$_2$-type structure. The normal state is a Pauli paramagnetic metal involving unusually heavy electrons with an effective mass of around 14--20 $m_e$ \cite{Wang13a}. Thermal conductivity measurement suggests that TlNi$_{2}$Se$_{2}$\ possesses multiple, nodeless superconducting gaps \cite{Hong14}. The material can be tuned by replacing Se with S, and Wang \etal\ have reported a smooth but non-monotonic variation of $T_c$ as a function of sulphur concentration $x$ in the isostructural series TlNi$_2$Se$_{2-x}$S$_x$ \cite{Wang13b}. Magnetic susceptibility and electrical resistivity measurements in the normal state up to 300~K did not detect any signature of additional phase transition for the entire substitution series, in contract to the closely related system KNi$_2$S$_2$ where a number of structural transitions were observed in the normal state \cite{Neilson13}.

In TlNi$_2$Se$_{2-x}$S$_x$ sulphur substitution is expected to introduce chemical pressure into the system, due to the smaller ionic radius of sulphur. At the same time this process introduces disorder into the system, as quantified by a significant reduction of the residual resistivity ratio (RRR=$\rho_{300{\rm K}}/\rho_{4{\rm K}}$) from $\sim100$ for $x=0$ to less than 10 for $x=1$ and $x=2$. Clearly, the $x$-dependence of $T_c$ not only reflects the effect of chemical pressure, but also the pair-breaking effects of disorder \cite{Wang13b}. The application of hydrostatic pressure on the other hand separates these two effects, and allows us to probe directly the intrinsic pressure dependence of $T_c$.

Single crystals of TlNi$_{2}$Se$_{2}$, TlNi$_{2}$SeS\ and TlNi$_{2}$S$_{2}$\ were synthesized using the self-flux method as described elsewhere \cite{Wang13a, Wang13b}. To track the superconducting transition under pressure, we implemented a two-coil technique in a Moissanite anvil cell, in which a 140-turn modulation coil was wound around the Moissanite anvil and a 10-turn pickup coil was placed inside the gasket hole together with the crystal \cite{Alireza03, Goh08, Klintberg10}. Owing to the advantageous volume filling factor of this setup, typically $\sim$30\%, the signal from the superconducting transition is clear, enabling us to follow the evolution of the superconducting state under pressure accurately. A Helium-3 dipper provided the low temperature environment. To avoid heating, we used a small modulation current of 1~mA, which gave a modulation field of around 1~Oe with a modulation frequency of 1.1~kHz. Glycerin was used as the pressure medium to provide hydrostatic pressure \cite{Osakabe08}, and the pressure was determined via ruby fluorescence spectroscopy.

\begin{figure}[!t]\centering
      \resizebox{8.5cm}{!}{
              \includegraphics{Figure1.pdf}}
              \caption{\label{fig1} (Color online) Selected AC susceptibility traces showing the superconducting transitions for (a) TlNi$_{2}$Se$_{2}$, (b) TlNi$_{2}$SeS\ and (c) TlNi$_{2}$S$_{2}$. The dotted lines in (a) illustrates our determination of the superconducting transition temperature.}
\end{figure}

Fig. \ref{fig1}(a) shows the temperature dependence of the normalised pick-up voltage for TlNi$_{2}$Se$_{2}$. This voltage is proportional to the AC susceptibility of the sample, and the transition to the superconducting state is clearly marked by the drop in the pick-up voltage. At ambient pressure, the transition temperature $T_c$, defined as the onset of the transition, is $3.5$~K. This agrees well with the reported value of $3.7$~K determined by resistivity and heat capacity \cite{Wang13a}. With applied pressure, $p$, $T_c$ decreases with an initial slope of d$T_c/$d$p\sim-59$~mK/kbar. This trend continues to about 22~kbar, followed by a rapid suppression of the superconducting state (c.f. Fig. \ref{fig2}a).

In the isostructural substitution series TlNi$_2$Se$_{2-x}$S$_x$, $T_c$ evolves smoothly as a function of the sulphur content $x$. As shown in the inset to Fig. \ref{fig2}(c), $T_c$ first decreases for $x<1$, and remains more or less constant between $x=1$ and $x=1.6$. For $x>1.6$, $T_c(x)$ shows a gentle positive slope. The replacement of Se by S is accompanied by a monotonic decrease of the lattice parameters $c$ and $a$ \cite{Wang13b}: moving from $x=0$ to $x=2$ results in a $\sim$2\% (5\%) reduction in the lattice parameter $a$ ($c$). If we consider the substitution of S as only providing a chemical pressure to the system, our observation of the collapse of the superconducting state in TlNi$_{2}$Se$_{2}$\ at $\sim$22~kbar is unexpected.

\begin{figure}[!t]\centering
       \resizebox{9cm}{!}{
              \includegraphics{Figure2.pdf}}
              \caption{\label{fig2} (Color online) Temperature-pressure phase diagram for (a) TlNi$_{2}$Se$_{2}$, (b) TlNi$_{2}$SeS\ and (c) TlNi$_{2}$S$_{2}$. For each composition, we conducted two sets of measurements to check the reproducibility of the results. Closed and open symbols denoted the first and the second runs, respectively. For TlNi$_{2}$Se$_{2}$, the numbers next to the closed circles denote the sequence of the measurements. For the other runs, the measurements were taken in the order of increasing pressures. The inset to (c) shows the $x$ dependence of $T_c$ at ambient pressure in TlNi$_2$Se$_{2-x}$S$_x$ \cite{Wang13b}}
\end{figure}

\begin{figure*}[!t]\centering
       \resizebox{17cm}{!}{
              \includegraphics{Figure3.pdf}}
              \caption{\label{fig3} (Color online) (a) Electronic band structure of TlNi$_{2}$Se$_{2}$\ along selected high symmetry directions, giving rise to four Fermi surface sheets as shown in (b) -- (e). Two of the four Fermi surface sheets of TlNi$_{2}$S$_{2}$\ are shown in (f) and (g). For the illustrations of the Fermi surface sheets, the $\Gamma$-point is located at the centre, and the $Z$-point is at the top of the Brillouin zone directly above $\Gamma$. The colored side of each Fermi surface corresponds to occupied states.}
\end{figure*}
The pressure dependence of $T_c$ in TlNi$_{2}$S$_{2}$\ and TlNi$_{2}$SeS\ reveals further unusual features. Figs. \ref{fig1}(b) and \ref{fig1}(c) display the corresponding AC susceptibility data for each of these materials,  allowing us to construct the $T$~--~$p$ phase diagrams as summarized in Fig. \ref{fig2}. If the effect of S substitution were to provide only chemical pressure, the application of physical pressure on fully S-substituted compound, TlNi$_{2}$S$_{2}$, would extend the trend of $T_c(x)$ near $x=2$, i.e. $T_c$ would be enhanced under pressure. Instead, $T_c$ decreases rapidly under pressure with a large initial slope of $-68$~mK/kbar. An even more surprising evolution of $T_c(p)$ is observed in TlNi$_{2}$SeS: $T_c$ first \textit{increases} with an initial slope
of $+27$~mK/kbar, reaches a maximum value of $\sim$1.3~K before it plummets, thus forming a dome-shaped $T_c(p)$ frequently observed in many correlated electron systems.

To check the reproducibility of our results, we conducted two sets of measurements using different samples for each composition. For TlNi$_{2}$Se$_{2}$\ and TlNi$_{2}$SeS, the agreement between sets of measurements is excellent. For TlNi$_{2}$S$_{2}$, the second set (Sample 2) gives an overall higher $T_c(p)$, which we attribute to a slight difference in sample quality. Importantly, the curvature of $T_c(p)$ is identical for both sets of data. In fact, these two sets of data can be brought to a broad agreement by including a relative offset of $\sim$0.4~K. These measurements demonstrate the robustness of our data, and place confidence on our observation of unusual, multi-dome $T_c(p)$ across the TlNi$_2$Se$_{2-x}$S$_x$ series.

A striking similarity across these $T$~--~$p$ phase diagrams concerns the abrupt disappearance of the superconducting state under high pressure. It might be tempting to attribute the disappearance to the onset of pressure inhomogeneity, however we find this unlikely due to previous work using a similar experimental technique. A recent study performed by some of us on BaFe$_2$(As,P)$_2$ up to 70~kbar \cite{Klintberg10} used the same pressure transmitting fluid, same type of anvil cell and similar gasket thickness. However, $T_c(p)$ varied smoothly in BaFe$_2$(As,P)$_2$ up to the highest pressure achieved. Moreover, in the present study the disappearance of superconductivity occurs at different pressures in samples with different sulphur content -- the pressure at which superconductivity collapses is $\sim$25~kbar, 40~kbar and 22~kbar in TlNi$_{2}$Se$_{2}$, TlNi$_{2}$SeS\ and TlNi$_{2}$S$_{2}$, respectively. As a consequence of these facts, it is difficult to associate the disappearance with pressure inhomogeneity.

To gain further insight, we examine the electronic structure of the stoichiometric end compounds TlNi$_{2}$Se$_{2}$\ and TlNi$_{2}$S$_{2}$. Bandstructure calculations were carried out with the WIEN2k \cite{Wien} package, which is based on density functional theory in the local density approximation. Exchange correlations were approximated with the Perdew-Enzerhoff approximation. Due to the rather heavy cores, additional variational steps accounting for spin-orbit coupling and relativistic local orbits were included. The calculations were performed with a resolution of 100,000 points in the first Brillouin Zone.

We compute the lattice parameters $a$ and $c$ for the fully relaxed structure, and we obtain $a=3.90$~\AA\ and $c=13.56$~\AA\ for TlNi$_{2}$Se$_{2}$, in excellent agreement with the experimental values $a=3.87$~\AA\ and $c=13.43$~\AA\ \cite{Wang13b}. For TlNi$_{2}$S$_{2}$, the calculated (experimental \cite{Wang13b}) lattice parameters are $a=3.81$~\AA\ ($a=3.79$~\AA) and $c=13.00$~\AA\ ($c=12.77$~\AA). Fig. \ref{fig3}(a) shows the dispersion relation of TlNi$_{2}$Se$_{2}$\ along selected high symmetry directions. Four bands cross the Fermi level giving rise to Fermi surface sheets A--D shown in Figs. \ref{fig3}(b)--(e). For TlNi$_{2}$S$_{2}$, the electronic structure is very similar to that of TlNi$_{2}$Se$_{2}$. In particular, Band A, Band B and Band D assume a very similar Fermi surface topology as the TlNi$_{2}$Se$_{2}$\ counterparts. However, it is important to note that the Fermi surface sheet associated with Band C becomes more rounded around the $Z$-point in TlNi$_{2}$S$_{2}$. Since the lattice parameters of TlNi$_2$Se$_{2-x}$S$_x$ evolve smoothly with $x$, it is thus very reasonable to assume that the Fermi surfaces of the system would evolve smoothly with $x$. This implies that, for TlNi$_{2}$SeS, the Fermi surface sheets associated with Band C would be more (less) rounded at $Z$ than the corresponding sheets in TlNi$_{2}$Se$_{2}$\ (TlNi$_{2}$S$_{2}$).

The progressive rounding of the Fermi surface sheet associated with Band C will have strong influence on the interband nesting probability with the Fermi surface sheets associated with Band D. The nesting of the well separated hole and electron sheets not only provides an intuitive picture to understand a spin-density-wave type antiferromagnetism, but also present a framework for the discussion of
spin-fluctuation mediated superconductivity \cite{Mazin08, Kuroki08}. Our calculations thus suggest the tuning of the strength of this nesting condition by varying the sulphur concentration or applying pressures.

The existence of multiple $T_c$ domes has been revealed in several iron-based systems, notably the $T_c(p)$ in FeSe$_{1-x}$ \cite{Miyoshi09, Bendele10} and the $T_c(x)$ in LaFeAsO$_{1-x}$H$_x$ \cite{Iimura12, Hiraishi14}. In these systems, the mechanism of superconductivity was unclear when the multi-dome $T_c$ was first reported \cite{Miyoshi09, Iimura12}. However, subsequent sensitive microscopic measurements detected the existence of antiferromagnetism bordering the superconducting phase \cite{Bendele10, Hiraishi14}, thereby allowing a unified treatment of superconductivity within the framework of spin fluctuation model. Our discovery of the multi-dome $T_c(p)$ in the TlNi$_2$Se$_{2-x}$S$_x$ bears striking resemblance to the phase diagrams of FeSe$_{1-x}$ and LaFeAsO$_{1-x}$H$_x$. With the picture of the Fermi surface nesting from our electronic structure calculations, it is urgently needed to investigate the magnetism of this system under pressure and over a wide range sulphur concentrations.

The clearly different phase diagrams produced through tuning by pressure and chemical substitution in the TlNi$_{2}$Se$_{2}$\ system could be accounted for in a number of ways. It may be for instance that the inequivalence arises from the sensitivity of the superconducting state to disorder induced pair-breaking. The most likely pairing scenario for TlNi$_{2}$Se$_{2}$\ is a multi-gapped nodeless s$\pm$ wave state, supported by recent thermal conductivity measurements \cite{Hong14}. In this state $T_c$ is expected to be rapidly suppressed with disorder \cite{Onari09}, a result of a breakdown of Anderson's theorem arising from interband scattering between sign-reversing Fermi surface sheets. The disorder introduced through isoelectronic substitution of S for Se could be felt in this manner. A second possibility is that there exists an as-yet-unknown structural or magnetic phase transition that occurs under high pressure. Our results motivate further work on the high pressure structural properties and magnetic properties of this family.

\textbf{Acknowledgment.} We acknowledge funding support from the EPSRC, Trinity College (Cambridge), CUHK (Startup Grant, Direct Grant No. 4053071), UGC Hong Kong (ECS/24300214), Cusanuswerk and the Royal Society. The work in ZJU and HNU was supported by the Natural Science Foundation of China (Grants No. 11374261 and No. 11204059), the Ministry of Science and Technology of China (National Basic Research Program No. 2011CBA00103, No. 2012CB821404, and No. 2015CB921004).

Superconductivity (cond-mat.supr-con); Strongly Correlated Electrons (cond-mat.str-el)


Nodeless superconductivity in Ir$_{1-x}$Pt$_x$Te$_2$ with strong spin-orbital coupling

Abstract: The thermal conductivity $\kappa$ of superconductor Ir$_{1-x}$Pt$_{x}$Te$_2$ ($x$ = 0.05) single crystal with strong spin-orbital coupling was measured down to 50 mK. The residual linear term $\kappa_0/T$ is negligible in zero magnetic field. In low magnetic field, $\kappa_0/T$ shows a slow field dependence. These results demonstrate that the superconducting gap of Ir$_{1-x}$Pt$_{x}$Te$_2$ is nodeless, and the pairing symmetry is likely conventional s-wave, despite the existence of strong spin-orbital coupling and a quantum critical point.

\section{Introduction}

The effect of strong spin-orbital coupling (SOC) on
superconductivity has recently attracted much attention. One example
is the topological superconductor, such as candidate
Cu$_x$Bi$_2$Se$_3$ in which Cu atoms are intercalated into
topological insulator Bi$_2$Se$_3$ with strong SOC. \cite{YSHor}
Novel superconducting state was claimed in Cu$_x$Bi$_2$Se$_3$ by the
point-contact spectra and superfluid density measurements.
\cite{SSasaki,MKriener1} Another example is the noncentrosymmetric
superconductor, such as Li$_2$Pt$_3$B in which the spatial inversion
symmetry is broken. \cite{PBadica} The strong SOC in Li$_2$Pt$_3$B
gives large spin-triplet pairing component and produces line nodes
in the superconducting gap. \cite{HQYuan,MNishiyama}

More recently, superconductivity was discovered in the layered
compound IrTe$_2$ by Pd intercalation (Pd$_x$IrTe$_2$),
\cite{JJYang} Pd substitution (Ir$_{1-x}$Pd$_x$Te$_2$),
\cite{JJYang} Pt substitution (Ir$_{1-x}$Pt$_x$Te$_2$), \cite{SPyon}
or Cu intercalation (Cu$_x$IrTe$_2$). \cite{MKamitani} Since the SOC
is proportional to $Z^4$, where $Z$ is the atomic number, the
superconductivity in doped IrTe$_2$ must associated with strong SOC
due to the large $Z$.

Furthermore, the parent compound IrTe$_2$ exhibits an intriguing
structural phase transition from a high-temperature trigonal to a
low-temperature monoclinic phase near 270 K. \cite{JJYang} Initially
it was related to a charge-density-wave (CDW) induced by Ir $5d$
$t_{2g}$ orbitals, \cite{JJYang,DOotsuki} however, later no CDW gap
was detected from the optical spectroscopy \cite{AFFang} and
angle-resolved photoemission spectroscopy (ARPES)
\cite{DOotsukiARPES} measurements. At this moment, the origin of the
transition is still under hot debate, with proposals such as crystal
field effect from Te $5p$ orbital splitting, \cite{AFFang} the
depolymerization-polymerization of anionic Te bonds, \cite{YSOhPRL}
and Ir $5d$ orbital order. \cite{HBCao} With increasing the doping
level $x$, the structure transition is gradually suppressed and
superconductivity emerges, showing a dome-like phase diagram with
the maximum $T_c$ of 3 K near $x \approx$ 0.04. \cite{JJYang,SPyon}
Such a phase diagram of doped IrTe$_2$ is reminiscent of high-$T_c$
cuprates and some heavy fermion superconductors, in which
superconductivity appears close to a magnetic quantum critical point
(QCP). This means that there likely exists a QCP under the
superconducting dome of doped IrTe$_2$, and the superconductivity
may be unconventional. \cite{PMonthoux} Therefore it is of great
interest to investigate whether there is novel superconducting state
in doped IrTe$_2$.

The ultra-low-temperature thermal conductivity measurement is a bulk
tool to study the gap structure of superconductors.
\cite{Shakeripour} The existence of a finite residual linear term
$\kappa_0/T$ in zero field is usually considered as the signature of
nodal superconducting gap. Further information of nodal gap, gap
anisotropy, or multiple gaps may be obtained from the field
dependence of $\kappa_0/T$. \cite{Shakeripour} Previously,
single-gap $s$-wave superconductivity near the QCP of CDW has
been clearly shown in Cu$_x$TiSe$_2$ by thermal conductivity
measurements. \cite{SYLi1}

In this paper, we probe the superconducting gap structure of
Ir$_{1-x}$Pt$_{x}$Te$_2$ ($x$ = 0.05) single crystal by measuring
the thermal conductivity $\kappa$ down to 50 mK. The residual linear
term $\kappa_0/T$ is negligible in zero magnetic field. The field
dependence of $\kappa_0/T$ is slow at low field, unlike that of a
nodal superconductor. Both results suggest nodeless superconducting
gap in Ir$_{1-x}$Pt$_{x}$Te$_2$.

\section{Experimental}

Single crystals of Ir$_{1-x}$Pt$_x$Te$_2$ were grown via self-flux
method. \cite{AFFang} The dc magnetic susceptibility was measured by
using a SQUID (MPMS, Quantum Design). The heat capacity measurement
was preformed in a physical property measurement system (PPMS,
Quantum Design) via the relaxation method. The
Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single crystal was cut to a rectangular
shape of dimensions 2.0 $\times$ 0.55 mm$^2$ in the $ab$ plane and
20 $\mu$m thickness along the $c$ axis. Four silver wires were
attached to the sample surface with silver paint, which were used
for both in-plane resistivity and thermal conductivity measurements.
The contacts are metallic with typical resistance 15 m$\Omega$ at 2
K. In-plane thermal conductivity was measured in a dilution
refrigerator, using a standard four-wire steady-state method with
two RuO$_2$ chip thermometers, calibrated {\it in situ} against a
reference RuO$_2$ thermometer. Magnetic fields were applied along
the $c$ axis and perpendicular to the heat current. To ensure a
homogeneous field distribution in the sample, all fields for
resistivity and thermal conductivity measurements were applied at
temperature above $T_c$.

\begin{figure}
\includegraphics[clip,width=7cm]{Fig1.eps}
\caption{(Color online). (a) Low-temperature magnetic susceptibility
of Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single crystal. The measurements
were performed in an applied field of $H$ = 20 Oe parallel to the
$ab$ plane with zero-field-cooled (ZFC) and field-cooled (FC)
processes. (b) Specific heat $C/T$ as a function of $T^2$ for
Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single crystals. The $C/T$ versus $T^2$
shows a linear dependence above $T_c$, giving the electronic
specific-heat coefficient $\gamma$ = 9.20 mJ/mol K$^2$. }
\end{figure}

\begin{figure}
\includegraphics[clip,width=9cm]{Fig2.eps}
\caption{(Color online). (a) Temperature dependence of the
resistivity $\rho$ for Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single crystal.
The data between 3.5 and 31 K can be fitted to $\rho(T) = \rho_0
+AT^n$, as shown by the solid line, with $\rho_0$ = 4.34 $\pm$ 0.002
$\mu\Omega$cm and $n$ = 3.01 $\pm$ 0.03. (b) The resistivity $\rho$
as a function of $T^{3.01}$ below 31 K. The solid line is the
fitting curve in (a). (c) Low-temperature resistivity of
Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single crystal in magnetic fields up to
0.75 T. (d) Temperature dependence of the upper critical field
$H_{c2}$, defined at the point $\rho$ dropping to zero on the curves
in (c). The dashed line is a guide to the eye, which points to
$H_{c2}(0) \approx$ 0.09 T.}
\end{figure}

\section{Results and Discussion}

Figure 1(a) presents the dc magnetic susceptibility of
Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single crystal. It was measured in
magnetic field $H$ = 20 Oe parallel to the $ab$ plane, with
zero-field-cooled (ZFC) and field-cooled (FC) processes. Sharp
superconducting transition with $T_c \approx$ 3.0 K and about 100\%
shielding volume fraction were observed for the ZFC process,
suggesting the homogeneous bulk superconductivity in our sample.

Heat capacity was measured on three pieces of
Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single crystals, with a total mass of
8.8 mg. It is plotted in Fig. 1(b), as $C/T$ versus $T^2$. Above
$T_c$, the data can be well fitted by $C/T = \gamma +\beta T^2$,
giving the electronic specific-heat coefficient $\gamma$ = 9.20
mJ/mol K$^2$. The significant jump was observed at $T_c \approx$ 3.0
K, which also indicates the high quality of our single crystals.

Figure 2(a) shows the resistivity of Ir$_{0.95}$Pt$_{0.05}$Te$_2$
single crystal in zero field. No resistivity anomaly is observed
above $T_c$, suggesting that no structural transition occurs in this
sample near optimal doping. The data between 3.5 and 31 K can be
fitted to $\rho(T) = \rho_0 +AT^n$, with $\rho_0$ = 4.34 $\pm$ 0.002
$\mu\Omega$cm and $n$ = 3.01 $\pm$ 0.03. To see more clearly, Fig.
2(b) plots $\rho$ versus $T^{3.01}$ and the solid line represents
the fitting curve data in Fig. 2(a). Such a temperature dependence
of $\rho(T) \sim T^n$ with $n \approx 2.8$ has been observed in
Ir$_{1-x}$Pt$_x$Te$_2$ polycrystal and attributed to phonon-assisted
interband scattering, \cite{SPyon} as in TiSe$_2$.
\cite{AFKusmartseva}

Previously, the upper critical field $H_{c2}(0) \approx$ 0.17 T has
been determined for Ir$_{0.96}$Pt$_{0.04}$Te$_2$ polycrystal by
resistivity measurements. \cite{SPyon} In order to obtain the
$H_{c2}(0)$ of our Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single crystal, we
also measure its resistivity with magnetic field parallel to the $c$
axis up to $H$ = 0.75 T and down to 50 mK, shown in Fig. 2(c). In
zero field, the resistivity drops to zero at 2.94 K with a narrow
transition width of 0.06 K. The superconducting transition is
gradually suppressed in magnetic fields. In Fig. 2(d), we plot the
temperature dependence of $H_{c2}$, where $T_c$ is defined at the
resistivity $\rho$ dropping to zero on the curves in Fig. 2(c). The
dashed line is a guide to the eye, which points to $H_{c2}(0)
\approx$ 0.09 T. This value is only about half of that in
Ir$_{0.96}$Pt$_{0.04}$Te$_2$ polycrystal. \cite{SPyon} It is not
surprising since in polycrystal, due to the random grain
orientation, the $H_{c2}$ represents the maximum value for all field
configurations, presumably $H \parallel ab$.

\begin{figure}
\includegraphics[clip,width=7.5cm]{Fig3.eps}
\caption{(Color online). Low-temperature thermal conductivity of
Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single crystal in zero and magnetic
fields. The solid lines are fits to $\kappa/T = a + bT^{\alpha-1}$.
The dash line is the normal-state Wiedemann-Franz law expectation
$L_0$/$\rho_0$, with the Lorenz number $L_0$ = 2.45 $\times$
10$^{-8}$ W$\Omega$K$^{-2}$ and $\rho_0$ = 4.34 $\mu\Omega$cm.}
\end{figure}

\begin{figure}
\includegraphics[clip,width=7.73cm]{Fig4.eps}
\caption{(Color online). Normalized residual linear term
$\kappa_0/T$ of Ir$_{0.95}$Pt$_{0.05}$Te$_2$ as a function of
$H/H_{c2}$. For comparison, similar data are shown for the clean
$s$-wave superconductor Nb, \cite{Lowell} the dirty $s$-wave
superconducting alloy InBi, \cite{Willis} an overdoped $d$-wave
cuprate superconductor Tl-2201, \cite{Proust} and the single-gap
$s$-wave superconductor Cu$_{0.06}$TiSe$_2$. \cite{SYLi1}}
\end{figure}

The thermal conductivities of Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single
crystal in zero and magnetic fields up to $H$ = 0.08 T are plotted
in Fig. 3, as $\kappa / T$ versus $T$. We fit all the curves to
$\kappa/T$ = $a + bT^{\alpha-1}$, in which the two terms $aT$ and
$bT^{\alpha}$ represent contributions from electrons and phonons,
respectively. \cite{Sutherland1,SYLi2} The power $\alpha$ of the
second term contributed by phonons is typically between 2 and 3 for
single crystals, due to the specular reflections of phonons at the
boundary. \cite{Sutherland1,SYLi2} In zero field, the fitting gives
a residual linear term $\kappa_0/T \equiv a$ = -7 $\pm$ 18 $\mu$W
K$^{-2}$ cm$^{-1}$, with $\alpha$ = 2.91 $\pm$ 0.04. Since our
experimental error bar is 5 $\mu$W K$^{-2}$cm$^{-1}$, the
$\kappa_0/T$ in zero filed is essentially zero, comparing to the
normal-state Wiedemann-Franz law expectation $L_0$/$\rho_0$ = 5.65
mW K$^{-2}$ cm$^{-1}$, with the Lorenz number $L_0$ = 2.45 $\times$
10$^{-8}$ W$\Omega$K$^{-2}$ and $\rho_0$ = 4.34 $\mu\Omega$cm.

For $s$-wave nodeless superconductors, there are no fermionic
quasiparticles to conduct heat when $T \rightarrow 0$ since all
electrons become Cooper pairs. Therefore, there is no residual
linear term of $\kappa_0/T$, as seen in Nb \cite{Lowell} and InBi.
\cite{Willis} However, a finite $\kappa_0/T$ in zero field is
usually observed in a superconductor with nodal gap, coming from the
nodal quasiparticles. \cite{Shakeripour} For example, $\kappa_0/T$ =
1.41 mW K$^{-2}$ cm$^{-1}$ for the overdoped cuprate
Tl$_2$Ba$_2$CuO$_{6+\delta}$ (Tl-2201), a $d$-wave superconductor
with $T_c$ = 15 K, \cite{Proust} and $\kappa_0/T$ = 17 mW K$^{-2}$
cm$^{-1}$ for the ruthenate Sr$_2$RuO$_4$, a $p$-wave superconductor
with $T_c$ = 1.5 K. \cite{Suzuki} The negligible $\kappa_0/T$ of our
Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single crystal in zero field suggests a
nodeless superconducting gap.

The field dependence of $\kappa_0/T$ will give more information of
the superconducting gap structure. For a clean type-II $s$-wave
superconductor with isotropic gap, $\kappa_0/T$ should grow
exponentially with field (above $H_{c1}$), as in Nb. \cite{Lowell}
In the case of nodal superconductor, $\kappa_0/T$ increases rapidly
($\sim H^{1/2}$) in low field due to the Volovik effect,
\cite{Volovik} as in Tl-2201. \cite{Proust}

In Fig. 3, all the curves in magnetic field are roughly linear.
Therefore we fit these curves to $\kappa/T$ = $a + bT^{\alpha-1}$
with $\alpha$ fixed to 2. Previously $\alpha \approx$ 2.2 was found
in Cu$_{0.06}$TiSe$_2$, \cite{SYLi1} and recently $\alpha \approx$ 2
has been observed in some iron-based superconductors such as
BaFe$_{1.9}$Ni$_{0.1}$As$_2$, \cite{LDing} KFe$_2$As$_2$,
\cite{JKDong} and Ba(Fe$_{1-x}$Ru$_x$)$_2$As$_2$ single crystals.
\cite{XQiu} Note that in the field of $H$ = 0.068 T and 0.08 T,
$\kappa_0/T$ = 5.53 $\pm$ 0.03 and 5.49 $\pm$ 0.05 mW K$^{-2}$
cm$^{-1}$ were obtained, respectively. Both values are close to the
normal-state Wiedemann-Franz law expectation $L_0$/$\rho_0$ = 5.65
mW K$^{-2}$ cm$^{-1}$. We take $H$ = 0.068 T as its bulk
$H_{c2}(0)$. A slightly different $H_{c2}(0)$ does not affect our
discussion on the field dependence of $\kappa_0/T$ below.

In Fig. 4, the normalized $\kappa_0/T$ of
Ir$_{0.95}$Pt$_{0.05}$Te$_2$ single crystal is plotted as a function
of $H/H_{c2}$. For comparison, we also plot the data of the clean
$s$-wave superconductor Nb, \cite{Lowell} the dirty $s$-wave
superconducting alloy InBi, \cite{Willis} the $d$-wave cuprate
superconductor Tl-2201, \cite{Proust} and the $s$-wave
superconductor Cu$_{0.06}$TiSe$_2$. \cite{SYLi1} From Fig. 4, the
$\kappa_0/T$ of Ir$_{0.95}$Pt$_{0.05}$Te$_2$ shows a field
dependence similar to that of $s$-wave superconductors InBi and
Cu$_{0.06}$TiSe$_2$. This further supports that
Ir$_{0.95}$Pt$_{0.05}$Te$_2$ is a nodeless superconductor.

Nodeless gap does not essentially mean conventional
superconductivity. For example, the nodeless gap observed in
optimally doped iron-based superconductor may be unconventional
$s_{\pm}$-wave resulting from antiferromagnetic spin fluctuations.
\cite{Hirschfeld} However, here in Ir$_{1-x}$Pt$_x$Te$_2$, the
nodeless gap is unlikely $s_{\pm}$-wave. One obvious reason is that
Ir$_{1-x}$Pt$_x$Te$_2$ does not have that kind of multiple Fermi
surfaces as in iron-based superconductors. \cite{JJYang,Hirschfeld}
Therefore, the pairing symmetry in Ir$_{1-x}$Pt$_x$Te$_2$ is likely
conventional $s$-wave. In this sense, the appearance of
superconductivity may have little relationship with quantum
fluctuations near QCP. This is not surprising since the single-gap
$s$-wave superconductivity in Cu$_x$TiSe$_2$ is also not related to
quantum fluctuations of CDW.

From the aspect of strong SOC, an unconventional odd-parity pairing
state was claimed for the topological superconductor candidate
Cu$_x$Bi$_2$Se$_3$. \cite{SSasaki,MKriener1} Both nodeless gap or
gap with point nodes are allowed for the odd-parity superconducting
state. \cite{SSasaki} In this context, more experiments such as
point-contact spectra are needed to completely exclude novel
superconductivity in Ir$_{1-x}$Pt$_x$Te$_2$.

\section{Summary}

In summary, we investigate the superconducting gap structure of
Ir$_{1-x}$Pt$_{x}$Te$_2$ ($x$ = 0.05) single crystal by thermal
conductivity measurements. The $\kappa_0/T$ in zero field is
negligible, and the field dependence of $\kappa_0/T$ is slow at low
field. Both of them suggest nodeless superconductivity in
Ir$_{1-x}$Pt$_{x}$Te$_2$. The pairing symmetry is likely
conventional $s$-wave, although the odd-parity superconducting state
can not be completely excluded from our measurements.

\begin{center}
{\bf ACKNOWLEDGEMENTS}
\end{center}
We thank J. J. Yang, Y. S. Oh, and S-W. Cheong for providing
Ir$_{0.96}$Pt$_{0.04}$Te$_2$ polycrystal and nominal
Ir$_{0.8}$Pt$_{0.2}$Te$_2$ single crystal to initialize this study.
This work is supported by the Natural Science Foundation of China,
the Ministry of Science and Technology of China (National Basic
Research Program No: 2009CB929203 and 2012CB821402), and the Program
for Professor of Special Appointment (Eastern Scholar) at Shanghai
Institutions of Higher Learning. \\

$^*$ E-mail: shiyan$\_$li@fudan.edu.cn