Giant pulses of pulsar radio emission

Abstract: Review report of giant pulses of pulsar radio emission, based on our detections of four new pulsars with giant pulses, and the comparative analysis of the previously known pulsars with giant pulses, including the Crab pulsar and millisecond pulsar PSR B1937+21.

\section{Introduction}

Giant pulses (GPs)- a short duration outbursts- are a special form
of pulsar radio emission.

This is a most striking phenomenon of the pulsar radio emission.
For normal pulsars, the intensity of single pulses varies by no
more then one order of magnitude (Hesse \& Wielebinski 1974,
Ritchings 1976).  GPs peak flux densities can exceed hundreds and
thousands of times the peak flux density of regular pulses. GPs
are the brightest sources of the radio emission  observed in the
Universe.

This  rare phenomenon has been detected only in 11 pulsars among
more than 1\,500 known ones.
 The firsts two of them to be discovered were Crab pulsar
  PSR B0531+21  (Staelin \& Reifenstein 1968; Argyle \& Gower 1972;
  Gower \& Argyle 1972) and the millisecond pulsar PSR B1937+21
  (Wolsczcan, Cordes \& Stinebring 1984), and for a long time
  remain little investigated.

For over 20 years only these two pulsars were known to emit GPs.
The last 5 years were marked by a systematic search of GPs that
benefit a fast progress in the detection of nine new pulsar with
GPs. They are PSR B1821-24 (Romani \& Johnston 2001), PSR B1112+50
(Ershov \& Kuzmin 2003), PSR B0540-69 (Johnston \& Romani 2003),
PSR B0031-07 (Kuzmin, Ershov \& Losovsky 2004), PSR J0218+42
(Joshi et al. 2004), PSR B1957+20 (Joshi et al. 2004), PSR
J1752+2359 (Ershov \& Kuzmin 2005), PSR J1823-3021A (Knight,
Bailes et al. 2005) and PSR B0656+14 (Kuzmin \& Ershov 2006).

\section{Observations }
GPs are observed as a rare very strong pulses of pulsar radio
emission standing out of the noise background and underlying
ordinary individual pulses. An example of one GP of the Crab
pulsar observed at frequency 111 \,MHz inside of 150 pulsar
periods is demonstrated in Fig.1.
\begin{figure}
   \resizebox{\hsize}{!}{\includegraphics{kuzmin_fig1.eps}}
      \caption{ An example of a GP of the Crab pulsar inside o
      f 150 normal pulses masked by noise}
      \label{Fig1}
   \end{figure}

GPs are distinguished by a number of characteristic features.
They possess very large excess of the flux density and energy of
radio emission relative to an average pulse (AP), the energy
distribution of GPs has a  power-law statistic, have a short pulse
time-scale and occur in a narrow-phase window of an AP (except of
the Crab pulsar for which GPs can occur anywhere within the AP).

\subsection{Flux density, energy, spectra and intensity distribution.}
GPs are the most extremal phenomena of pulsar radio emission.
Their peak flux densities can exceed hundreds and thousands of
times the peak flux density of an AP. An example of a GP in comparison with an
AP of the millisecond pulsar PSR B1937+21 is shown in Fig.2. The
peak flux densitiy of GPs in this pulsar at 1650\,MHz ranges up to
65 \,kJy and exceeds the peak flux density of an AP up to factor
of $3 \times 10^5$ (Soglasnov et al. 2004). An energy of this GP
$E = S_{peak} \times \tau $, where $\tau$ is the pulse duration,
exceeds the energy of the AP by a factor of 60.

\begin{figure}
\centerline{\psfig{file=kuzmin_fig2.eps,width=7.0cm,clip=} }
\caption{ An example of a GP (solid line) in comparison
with an AP (dotted line) of the pulsar PSR B1937+21. \it (After Soglasnov et al. 2005).
\rm Note different y-scales for a GP (left) and an AP (right).
\label{Fig2}}
\end{figure}
The peak flux densities of GPs from the Crab Nebula pulsar PSR
B0531+21 at 2228\,MHz exceeds the peak flux density of an AP  up
to factor of $5 \times 10^5$. An energy $E = S_{peak} \times \tau
$ of the GPs exceeds the energy of the AP by a factor of 80
(Kostyuk et al., 2003).

GPs of the Crab pulsar have been detected in a very wide frequency
range from 23 \,MHz (Popov et al. 2006) up to 15 GHz (Hankins et al. 2000).
Only small part of GPs occurs simultaneously in different frequency
ranges. Radio spectra of these GPs were studied by simultaneous
multi-frequency observations.
Sallmen et al. (1999) observations of GPs in pulsar PSR
B0531+21 at two frequencies 1.4 and 0.6 \,GHz show that the GPs
spectral indices fall between -2.2 and -4.9, which may be compared
to the AP value for this pulsar -3.0. Popov et al.(2006)
observations of this pulsar at three frequencies 600, 111 and 43
\,MHz reveal that the GPs spectral indices fall between -1.6 and
-3.1 with mean value -2.7, that also may be compared to the AP
value for this pulsar. Simultaneous two-frequency observations of
GPs from PSR B1937+21 at 2210-2250 and 1414-1446 MHz (Popov \&
Stappers 2003) don't reveal any GPs which occur simultaneously in
both frequency ranges. They conclude that radio spectra of
detected GPs are limited in frequency at a scale of about $\Delta
\nu /\nu < 0.5$. Simultaneous two-frequency observations of GPs in
PSR B0031-07 at 111 and 40 MHz reveal only a few common GPs
(Kuzmin et al. 05).

Kinkhabwala et al. (2000) have estimated the average spectral
properties of the GPs emission of the pulsar PSR B1937+21. Using
the top eight GPs at three frequencies 430, 1420 and 2380 \,MHz,
they find a somewhat steeper slope of -3.1 for the GPs spectrum,
compared the -2.6 slope for the normal emission spectrum of this
pulsar.

Another distinguishing characteristic of pulsars with GPs, as
demonstrated in Fig.3, is their two-mode pulse intensity
distribution. At low intensities of ordinary pulses, the pulse
strength distribution is Gaussian one, but above a certain
threshold the pulse strength of GPs is roughly power-law
distributed (Argyle \& Gower 1972, Cognard et al. 1996) .

\begin{figure}
\centerline{\psfig{file=kuzmin_fig3.eps,width=6.5cm,clip=} }
\caption{\textbf{a)} Cumulative distribution of GPs intensity of pulsar
PSR B1937+21 measured in units of the mean pulse and interpulse (after Cognars et al. 1996)
\label{Fig3}}
\end{figure}
\subsection{Size and brightness temperature of an emitting region.}

The next peculiarity of GPs, as demonstrated in Fig.2 for PSR
B1937+21, is a much shorter time-scale of GP as compared with an
AP. The width of  this GP is $\tau \le 15$ \,ns (Soglasnov et al.
2004).

Hankins et al. (2003) found in pulsar PSR B0531+21 the pulse
structure as short as 2 \,ns. If one interprets this pulse
duration in terms of the maximum possible size of emitting region
$l \le c \times \tau$, where $c$ is the speed of a light, the
time-scale
$\tau=$2 \,ns corresponds to a light-travel size of an
emitting body $l$ of only 60 \,cm, the smallest object ever
detected outside our solar system.

The brightness temperature of the GPs is
$$
  T_{\rm B} = S_{peak} \lambda^2 /2k\Omega~,
$$
where $S_{peak}$ is the peak flux density, $\lambda$ is the radio
wavelength, $k$ is the Boltzmann's constant, and  $\Omega \simeq
(l/d)^2$ is the solid angle of the radio emission region.

For observed in Crab pulsar at frequency 5.5\,GHz peak flux
density of $10^3$\,Jy  with $\tau = 2$\,ns (Hankins et al. 2003)
the brightness temperature of GPs is as high as $10^{37}$ \,K.

Soglasnov et al. (2004) have proved that the majority of GPs from
the millisecond pulsar PSR B1937+21 are shorter than 15 \,ns. At a
frequency of 1.65 \,GHz a peak flux density of $65 \times 10^3$
\,Jy with $\tau = 15$ \,ns show that the brightness temperature of
GPs in this pulsar is $T_B \ge 5 \times 10^{39}$ \,K, the
brightest radio emission observed in Universe.

However, these evaluations of  size  and brightness temperature of
an emitting body are not unambiguous.  Gil \& Melikidze (2004)
argued  that the apparent duration of the observed impulse
$\tau_{obs}$ may be shorter $^{1)}$ than the duration of the
emitted one as $\tau_{rad}$ as $\tau_{obs} = \tau_{rad} \times
\gamma^{-2}$. For Lorentz factor $\gamma \approx 100$, the
time-scale of pulse structure for B1937+21 will transformed from
observed $\tau_{obs}$= 2 ns to emitted $\tau_{rad}$ =20 $\mu s$.
\footnote{Relativistic shortening of a pulse duration were claimed
by Smith (1970) and Zheleznyakov (1971).} As well as $T_B \propto
\Omega^{-1} \propto \tau^{-2}$, the brightness temperature will be
reduced to $T_B \approx 10^{31}$ \,K.

This aspect needs a further refinement.

The duration of GPs may be masked by the scatter broadening of a
pulse due to multi-path propagation in the interstellar plasma.

GPs are clustered around a small phase window (except of Crab
pulsar, for which GPs can occur anywhere within the AP). An
example of such clustering for pulsar PSR B0031-07 is shown in the
bottom of Fig.4.

\begin{figure}
\centerline{\psfig{file=kuzmin_fig4.eps,width=6.5cm,clip=} }
\caption{\textbf{(top)} GP in pulsar PSR B0656+14 at 111  MHz (bold line)
       together with the AP (dotted line). The plot of AP is presented on a
       500 times larger scale and flux densities of GP and AP are shown
       separately on the left and right sides of the "y" scales. \textbf{(bottom)} The
       phases of the observed GPs relative to the center of the
       AP. N is the number of GP.
       (after Kuzmin \& Ershov ,2006)
       \label{Fig4}}
\end{figure}

The top of Fig4 shows the GP together with the AP for the pulsar
PSR B0656+14 (Kuzmin \& Ershov 2006). The peak flux density of this
GP exceeds the peak flux density of AP by factor of 630.

Kuzmin \& Ershov (2004) two frequency observations of GPs in the
pulsar PSR B0031-07 reveal that some GPs have two components and
that the spacing of these components at 40 \,MHz is larger than at
111 \,MHz. This is similar to a larger width of the AP at lower
frequency, which is interpreted as a divergence of magnetic field
lines in the radius-to-frequency mapping in the hollow cone model
of pulsar radio emission. This suggests that  GPs from this pulsar
originate in the polar region rather than in the light cylinder
region.

One may suggests that there are two classes of GPs, one associated
with high-energy emission from the outer gap, the other associated
with the polar radio emission. GPs of PSR B0031--07 will be of the
second class.

\rm Table 1 summarizes the comparative data for known pulsars with
GPs. Here PSR is a pulsar name, \it P \rm is the pulsar period,
$B_{LC}$ is a magnetic field strength at the light cylinder, \it
Freq \rm is the observations frequency, $S_{GP}/S_{AP}$ is an
excess of the peak flux density of a strongest GP over the peak
flux density of an AP, $T_B$ is a brightness temperature of the
strongest GP radio emission and $E_{GP}/E_{GP}$ is an excess of
the energy of the strongest GP over the energy of an AP.

\begin{table*}[t]
  \caption[]{ Properties of the Giant Pulses}

  \label{Tab:publ-works}
  \begin{center}\begin{tabular}{clclclclc}
  \hline\noalign{\smallskip}
PSR&\it P \rm&$log B_{LC}$&\it Freq \rm& $S_{GP}/S_{AP}$ &$T_B$& $E_{GP}/E_{GP}$&References\\
    &ms &G & GHz  &             &  K   &  &  \\
      \hline\noalign{\smallskip}
 B0031-07&943&6.9&0.04&  400        &$\ge10^{28}$             & 75&Kuzmin \& Ershov 2004 \\
         & & &0.11&  120        &$\ge10^{26}$             &  8    &Kuzmin et al. 2004 \\
 J0218+42&2.3&$3.2\times 10^5$&0.61     &                 &1.3& 51&Joshi et al. 2004 \\
 B0531+21&33.4&$9.8\times 10^5$&0.14     & 300               &&   &Argyle \& Gower 1972\\
         & & &0.59&  $6\times 10^4$&$\ge10^{36}$   &  75          &Kostyuk et al. 2003 \\
         & & &2.23&  $5\times 10^5$&$\ge10^{34}$  &  80 &          Kostyuk et al. 2003 \\
         & & &5.50&                &$\ge10^{37}$  &        &       Hankins et al. 2003 \\
 B0540-69&50.5&$3.5\times 10^5$&1.38        &$5\times 10^3$& & &   Johnston \& Romani 2003 \\
 B0656+14&385&770 &0.11&  600   &$ \ge10^{26}$     & 110&          Kuzmin \& Ershov 2006 \\
 B1112+50&1656&4.1 &0.11&   80      &$ \ge10^{26}$      & 10&      Ershov \& Kuzmin 2003 \\
 J1752+2359&409&71 &0.11&  260      &$ \ge10^{28}$      &200&      Ershov \& Kuzmin 2005 \\
 B1821-24&3.0&$7.2\times 10^5$ &1.51&           &         & 81&    Romani \& Johnston 2001 \\
 J1823-3021A&5.4&$2.5\times 10^5$&0.68& 680  &             & 64&   Knight et al. 2005 \\
         & & &1.40&1700 &      &      28  &                        Knight et al. 2005 \\
 B1937+21&1.5&$9.8\times 10^5$& 0.11& 600&$\ge 10^{35}$  & 65&     Kuzmin \& Losovsky \\
         & & &1.65&  &$\ge 5\times 10^{39}$ &    60          &     Soglasnov et al. 2004 \\
 B1957+20&1.6&$3.8\times 10^5$ &0.61      &     &         &129&    Joshi et al. 2004\\
  \noalign{\smallskip}\hline
  \end{tabular}\end{center}
\end{table*}

\section{Are giant pulses inherent to some special property of pulsars ?}
The first detected pulsars with GPs PSR B0531+21 and PSR B1937+21
were among a small group of pulsars with highest magnetic field at
the light cylinder $B_{LC}=10^4 -10^5$\,G. This gave a rise to the
suggestion that GPs are inherent to pulsars with very strong
magnetic field at the light cylinder and the search of GPs was
oriented on those pulsars. As a result GPs were detected in five
other pulsars with very strong magnetic field at the light
cylinder: PSR B1821-24 (Romani \& Johnston 2001), PSR B0540-69
(Johnston \& Romani 2003), PSR J0218+42 (Joshi et al. 2004), PSR
B1957+20 (Joshi et al. 2004) and PSR J1823-3021 (Knight et al.
2005).

However detection of the GPs in four pulsar with ordinary magnetic
field at the light cylinder: PSR B0031-07 (Kuzmin, Ershov \&
Losovsky 2004), PSR B0656+14 (Kuzmin \& Ershov 2006), PSR
B1112+50 (Ershov \& Kuzmin 2003) and PSR J1752+2359 (Ershov \&
Kuzmin 2005) have revealed that GPs can occur in ordinary pulsars
too.

Johnson \& Romany (2004) claimed that GPs may be associated with
pulsars which possess a high energy X-ray emission. But the
alignment of X-ray pulses with GPs is observed only for four
objects among eleven known pulsars with GPs and is not a proper
indication for GPs.

Knight et all. (2005) argued that GPs may be indicated by large
rotation loss luminosity $\dot E \propto P^{-3} \dot P$. But in
fact the rotation loss luminosity of the known 11 pulsars with GPs
are differ by 6 orders of magnitude and  is not a proper
indication for GPs also.

\rm GPs occur in various types of pulsars in a wide range of
periods $P=1.5-1600$ \,ms and magnetic field at the light cylinder
$\log B_{LC}=4 -10^6$\,G over a wide range of radio frequencies.

\section{Summary}
\label{sect:summary}

Giant pulses  is a special  form of  pulsar radio emission, that
is characterized by very large excess of flux density and energy
of radio emission relative to an  average pulse, the power-law
statistic of the energy distribution , giant pulses occur in a
narrow-phase window of an average pulse and have a short pulse
time-scale.

The flux density of giant pulses increases over the flux density
of an average pulse $S_{GP}/S_{AP}$ up to  $5 \times 10^5$ .

Giant pulses energy excess over the energy of an average pulse is
$E_{GP}/E_{AP}$ = 50 - 200 and is nearly the same for different
magnetic field at the light cylinder, pulsar periods and
frequencies.

A light-travel size of an emitting body indicate the smallest
object ever detected outside our solar system.

Giant pulses are the brightest sources of radio emission  among
the known astronomical objects.

Giant pulses exist in various types of pulsars in a wide range of
periods , magnetic field at the light cylinder  and broad
frequency range.

\begin{acknowledgements}
This work was supported in part be the Russian Foundation of Basic
Research, Project N 05-02-16415.
\end{acknowledgements}

Giant pulses of pulsar radio emission


Determination of the phase diagram of the electron doped superconductor Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$

Abstract: Systematic measurements of the resistivity, heat capacity, susceptibility and Hall coefficient are presented for single crystal samples of the electron-doped superconductor Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$. These data delineate an $x-T$ phase diagram in which the single magnetic/structural phase transition that is observed for undoped BaFe$_2$As$_2$ at 134 K apparently splits into two distinct phase transitions, both of which are rapidly suppressed with increasing Co concentration. Superconductivity emerges for Co concentrations above $x \sim 0.025$, and appears to coexist with the broken symmetry state for an appreciable range of doping, up to $x \sim 0.06$. The optimal superconducting transition temperature appears to coincide with the Co concentration at which the magnetic/structural phase transitions are totally suppressed, at least within the resolution provided by the finite step size between crystals prepared with different doping levels. Superconductivity is observed for a further range of Co concentrations, before being completely suppressed for $x \sim 0.018$ and above. The form of this $x-T$ phase diagram is suggestive of an association between superconductivity and a quantum critical point arising from suppression of the magnetic and/or structural phase transitions.

\section{Introduction}

BaFe$_2$As$_2$ is a prototypical member of the new family of Fe-pnictides that play host to high temperature superconductivity. The stoichiometric compound suffers a coupled structural and antiferromagnetic transition at $\sim$ 140 K \cite{rotter}. Suppression of the broken-symmetry state by either applied pressure \cite{alireza} or chemical doping \cite{rotter-2,mandrus}  results in superconductivity, and one of the key questions associated with this entire class of material is the role played by spin fuctuations associated with the incipient tendency towards magnetism \cite{Norman, Mazin}.

Initial work on chemical substitution in BaFe$_2$As$_2$ focused on the hole-doped material Ba$_{1-x}$K$_x$Fe$_2$As$_2$ \cite{rotter-2}. In this case, it has been established that the structural/magnetic transition is totally suppressed for K concentrations above $x = 0.4$ \cite{xhchen-2}, whereas superconductivity appears for a wide range of concentrations from $x$ = 0.3 to 1 \cite{xhchen-2}. More recently, it has been shown that Co-substitution (i.e. Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$) also suppresses the magnetic and structural transitions, eventually leading to superconductivity \cite{mandrus}. Co-doped BaFe$_2$As$_2$ is especially attractive for two reasons. First, from a fundamental point of view, this material allows an exploration of the interrelation between the structural/magnetic phase transitions with superconductivity in the context of an electron-doped system. Secondly, from a more practical perspective, crystal growth using cobalt does not suffer from the inherent difficulties associated with using potassium. Cobalt has neither a large vapor pressure nor attacks the quartz tubing used to encapsulate the growths. It is thus possible to obtain more homogeneous crystals and more reproducible synthesis conditions than is otherwise feasible for potassium-doped samples.

In this paper we present results of resistivity, Hall coefficient, heat capacity and susceptibility measurements of single crystal samples of Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$. We find that the structural/magnetic phase transition which occurs at 134 K in BaFe$_2$As$_2$ is rapidly suppressed, and splits in to two successive phase transitions with increasing Co concentration. The signatures of these two distinct phase transitions are relatively sharp, and it is unlikely that this effect is associated with phase separation due to a variation in Co concentration across an individual sample.  Superconductivity appears for Co concentrations $0.025 < x < 0.18$, with a maximum $T_c$ for a Co concentration of $x \sim 0.06$, coincident with the concentration at which the structural/magnetic phase transitions are suppressed below T$_c$. The apparent coexistence of superconductivity with the broken-symmetry state on the underdoped side of the phase diagram is unlikely to be due to phase separation due to Co inhomogeneity, and mirrors what is seen for Ba$_{1-x}$K$_x$Fe$_2$As$_2$ \cite{xhchen-2}. These measurements establish an essential symmetry between electron and hole doping of BaFe$_2$As$_2$, in which the appearance of superconductivity in both cases is clearly associated with suppression of the magnetic/structural phase transitions. Nevertheless, there are some clear differences that emerge from this study. Specifically, Co-substitution is found to be significantly more effective at suppressing the structural/magnetic phase transitions than K-substitution, and the superconducting ``dome'' extends over a much smaller range of doping. The resulting phase diagram is highly suggestive that superconductivity is intimately related to the presence of a quantum critical point associated with suppression the structural and/or magnetic phase transitions.

\section{Results and Discussion}

Single crystals of Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$ were grown from a self flux using similar conditions to published methods \cite{xhchen,mandrus}. Ba and FeAs in the molar ration 1:4 with additional Co were placed in an alumina crucible and sealed in evacuated quartz tubes. The mixture was heated to 1150C, and held for 24 hours, before slowly cooling to 1000 C, at which temperature the remaining flux was decanted using a centrifuge. The crystals have a plate-like morphology, with the $c$-axis perpendicular to the plane of the plates, and grow up to several millimeters on a side. The Co concentration was measured by electron microprobe analysis (EMPA) using undoped BaFe$_2$As$_2$ and elemental Co as standards. Measurements were made at several locations on each sample. Measured values were close to the nominal melt composition in all cases, and the variation in Co concentration across individual samples was typically characterized by a standard deviation of $0.15\%$. For instance, the Co concentration corresponding to optimal doping is $x$ = 0.061 $\pm$ 0.002.

\begin{figure}[tbh]
\includegraphics[width=8.5cm]{RvsT}
\caption{\label{RvsT} Temperature dependence of the in-plane ($\rho_{ab}$) resistivity of Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$. Data are shown normalized by the room temperature resistivity, $\rho_{ab}$(300K) , to remove uncertainty in estimates of the absolute value due to geometric factors. Successive data sets are offset vertically by 0.1 for clarity (except for $x = 0$ data, which is offset by 1.4). Values of the Co-concentration $x$ are listed in the legend, and were determined by microprobe analysis for all concentrations except the specific sample with $x = 0.070$, for which $x$ was estimated based on the observed linear relation between nominal and actual $x$.}
\end{figure}

\begin{figure}[tbh]
\includegraphics[width=8.5cm]{CP}
\caption{\label{CP}(a) Heat capacity data for four representative Co concentrations ($x$ = 0, 0.016, 0.025, 0.036 and 0.061) showing the suppression of the structural/magnetic phase transitions with increasing Co concentration. (b) Superconducting anomaly for the sample with $x$ = 0.061 shown as $\Delta C/T$ vs. $T$, where $\Delta C = C_S - C_N$ and $C_N$ is estimated as described in the text from the heat capacity measured in an applied field of 14 T. Data are shown for zero field, and for applied fields of 1, 3 and 5 T along the $c$-axis.}
\end{figure}

The in-plane electrical resistivity ($\rho_{ab}$) was measured with a standard four-probe configuration using a Quantum Design Physical Properties Measurement System (PPMS), and the results are summarized in Fig. \ref{RvsT}. Data are shown normalized by the room temperature resistivity, $\rho_{ab}$(300K), to remove uncertainty in estimates of the absolute value due to geometric factors  \footnote{ Absolute values of $\rho_{ab}(300K)$ for three representative Co concentrations are 0.676 $m \Omega cm$ ($x$=0, undoped), 0.906 $m \Omega cm$ ($x$=0.061, optimally doped), and 0.696 $m\Omega cm$ ($x$=0.018, completely overdoped)}. The anomaly associated with the magnetic/structural phase transition in undoped BaFe$_2$As$_2$ is clearly visible at 134 K. This feature changes from a sharp drop in the resistivity below the transition temperature for BaFe$_2$As$_2$, to an abrupt upturn for Co-doped samples \footnote{The origin of the significant difference between the resistivity anomaly associated with the structural/magnetic transitions in undoped BaFe$_2$As$_2$ and Co-doped BaFe$_2$As$_2$ is unclear. Recent optical conductivity measurements indicate that a substantial fraction of the FS of BaFe$_2$As$_2$ is gapped below the SDW transition \cite{leo}. In this case, we can speculate that there is a subtle balance between effects associated with the reduction in the number of carriers and changes in the scattering. We note that the resistivity anomaly associated with the magnetic/structural phase transitions observed for intermediate Co concentrations is reminiscent of that which is observed for BaFe$_2$As$_2$ grown from a Sn flux \cite{canfield}, possibly indicating that Sn substitution also affects the scattering in a similarly sensitive manner. }. The anomaly is clearly suppressed in temperature with increasing $x$. For samples with $x > 0.035$ superconductivity is evident from a sharp drop to zero in the resistivity, while the anomaly associated with the structural/magnetic phase transition is still observable in the normal state. As $x$ increases beyond this value, the superconducting critical temperature T$_c$\ (defined here as the midpoint of the resistive transition) increases, whereas the magnetic/structural anomaly continues to be suppressed. Finally T$_c$\ reaches the highest value of 24 K in the sample with $x = 0.061$ (``optimal" doping), for which the normal state anomaly is no longer observable. Further increasing $x$ beyond this value results in an eventual suppression of superconductivity (``overdoping"). For the sample with $x = 0.184$, T$_c$\ is below our instrumental base temperature of 1.8 K.

The suppression of the magnetic/structural phase transitions was also followed by heat capacity, measured for single crystal samples by a relaxation technique using a Quantum Design PPMS. Results are shown in Fig. \ref{CP}(a) for four representative Co concentrations. For undoped BaFe$_2$As$_2$, a sharp peak at 134 K marks the first order structural/magnetic phase transition as reported previously \cite{dong}. With increasing $x$ both the transition temperature and the magnitude of the peak are rapidly suppressed, until the anomaly becomes almost unobservable for $x > 0.036$. In order to provide an estimate of the phonon background, the heat capacity was also measured for an optimally doped sample (i.e. $x$ = 0.061, which is the smallest Co concentration for which the magnetic/structural phase transition is totally suppressed). These data were taken to lower temperature, revealing the superconducting anomaly. Data for this particular Co-concentration are shown in Fig. \ref{CP}(b) as $\Delta C/T$ vs. $T$, where $\Delta C$ is the difference between the heat capacity in the superconducting state $C_S$ and the normal state $C_N$. The latter quantity has been estimated over the plotted temperature range from measurements made in a field of 14 T applied along the $c$-axis. These data indicate $\Delta C/T_c$ $\approx$ 23 (mJ/K$^2$ mol) in zero field.

\begin{figure}[tbh]
\includegraphics[width=8.5cm]{zoom}
\caption{\label{zoom}Resistivity (panels a-c) and heat capacity (panels d-f) data for samples with $x$ = 0, 0.016 and 0.025 over a narrower temperature interval, showing the splitting of the single phase transition observed for $x$ = 0 in to two distinct phase transitions for $x > 0$.}
\end{figure}

Careful inspection of the resistivity and heat capacity data shown in Fig.s \ref{RvsT} and \ref{CP} reveal that the single anomaly associated with the magnetic/structural phase transition which is seen for undoped BaFe$_2$As$_2$ splits in to two distinct features for the Co-doped samples. In fig \ref{zoom} we show the resistivity and heat capacity anomaly near these transitions for the undoped and two lightly doped samples. It is evident that the single sharp downturn in BaFe$_2$As$_2$ becomes two successive sharp changes in slope in the resistivity of the lightly doped Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$.  As for the heat capacity, the single sharp first-order peak observed for undoped BaFe$_2$As$_2$, becomes two distinct features, resembling a second order-like step followed by a lower temperature peak, reminiscent of a broadened first-order transition.

To further investigate the apparent splitting of the structural/magnetic phase transitions, magnetic susceptibility and Hall coefficient data were also measured. Results are shown in Fig \ref{G3} for two representative Co concentrations. Significantly, both measurements indicate the presence of two distinct phase transitions for $x > 0$ evidenced by distinct breaks in the slope of each quantity. In the case of the susceptibility, this effect is more clearly apparent as two peaks in the derivative $d(\chi T)/dT$ \footnote{Susceptibility data were smoothed using a linear interpolation and adjacent points averaging before taking the derivative.}. These data can be compared with the heat capacity measurements (right axis of the upper panels), where we have explicitly subtracted the normal state heat capacity of the optimal doped Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$ with $x = 0.061$ as described above, to provide an estimate of the contribution to the heat capacity arising from the two phase transitions. The background subtracted heat capacity is in excellent agreement with the derivative $d(\chi T)/dT$ of the susceptibility, which is proportional to the heat capacity near a second-order phase transition\cite{fisher-1962}. These data are also in excellent agreement with the derivative of the resistivity $d\rho/dT$. In short, all four measurements indicate unambiguously the presence of two distinct phase transitions. It is very unlikely that this splitting arises from Co inhomogeneity, which in it's simplest form would cause a continuous broadening rather than a distinct splitting of the two phase transitions. Indeed, the presence of Co inhomogeneity can be inferred from the broadening of each of these two phase transitions (Fig. \ref{CPDR}), but this effect remains less than the actual splitting.

\begin{figure*}[tbh]
\includegraphics[width=17cm]{G3}
\caption{\label{G3} Comparison of heat capacity, susceptibility, resistivity and Hall coefficient data for samples with (a) $x = 0$, (b) $x = 0.016$ and (c) $x = 0.025$. In the upper panels, the heat capacity (right axis) is shown as $\Delta C$, having subtracted the background phonon contribution as described in the main text, together with the in-plane susceptibility $\chi_{ab}$ (left axis) and the derivative $d(\chi T)/dT$ (arb. units). In the lower panels, absolute values of the in-plane resistivity $\rho_{ab}$ (left axis) are shown together with the Hall coefficient $R_H$ (right axis) and the derivative -$d\rho/dT$ (arb units). Vertical dotted lines mark successive phase transitions, evident in all four physical properties.}
\end{figure*}

\begin{figure}[tbh]
\includegraphics[width=8.5cm]{CPDR}
\caption{\label{CPDR}Progression of the two successive phase transitions, marked as $\alpha$ and $\beta$, as a function of Co concentration as observed in (a) heat capacity and (b) the derivative of the resistivity -$d\rho/dT$. Heat capacity data are shown as $\Delta C$, having subtracted the background phonon contribution as described in the main text. Data for $x$ = 0.025 and 0.036 have been multiplied by factors of 3 and 5 respectively for clarity.}
\end{figure}

Although the magnitude of the heat capacity anomalies associated with these two phase transitions is rapidly suppressed for $x > 0.036$, nevertheless the close correspondence between the heat capacity and the resistivity derivative can be used to follow these phase transitions to even higher Co concentrations. Data for Co concentrations up to $x$ = 0.051 are shown in Figure \ref{CPDR}. As can be seen, both transitions, marked as $\alpha$ and $\beta$, are broadened with increasing Co concentration, but the two specific features which can be associated with the critical temperatures for $x$ = 0.016, 0.025 and 0.036 remain visible in the resistivity derivative even for $x$ = 0.045 and 0.051. For larger values of the Co concentration, no phase transitions are evident above $T_c$.

Based on the measurements described above, we can establish a composition-temperature ($x-T$) phase diagram for Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$, shown in Fig. \ref{Phase}. Phase boundaries delineate the successive phase transitions $T_{\alpha}$, $T_{\beta}$ and $T_c$ as described above. In analogy to the oxy-pnictide LaFeAsO \cite{Dai-1} it is tempting to speculate that the upper transition $T_{\alpha}$ corresponds to the structural distortion, while the lower transition $T_{\beta}$ corresponds to the SDW transition, but the current data are not sufficient to make this association and further scattering and/or microscopic measurements will be necessary to determine the nature of the two transitions.

\begin{figure}[tbh]
\includegraphics[width=8.5cm]{Phase}
\caption{\label{Phase}Phase diagram for Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$ showing the suppression of the two successive phase transitions $\alpha$ and $\beta$ with increasing Co concentration, and the eventual superconducting (SC) "dome". Data points for $T_{\alpha}$ and $T_{\beta}$ were obtained from heat capacity, resistivity, Hall coefficient and susceptibility data for $x$ = 0, 0.016, 0.025 and 0.036, and from resistivity data alone for $x$ = 0.045 and 0.051. Superconducting $T_c$ values were obtained from resistivity data. Error bars indicate 10 and 90$\%$ of the resistive transition. Uncertainty in the Co concentration corresponds to less than the width of the data markers.}
\end{figure}

The apparent coexistence of superconductivity with the broken symmetry SDW state on the underdoped side of the phase diagram raises the question of whether there is macroscopic phase separation associated with a variation in the Co concentration. For length scales above one micron this can be ruled out based on the microprobe analysis. Specifically, the standard deviation in the cobalt doping is much smaller than the range across which we observed this coexistence. From a materials viewpoint it would also be surprising to find Co inhomogeneity on shorter length scales given the large solid solubility. This does not, however, rule out the intriguing possibility of a spontaneous phase separation in an otherwise homogeneous material. In this regard susceptibility measurements provide some insight.

\begin{figure}[tbh]
\includegraphics[width=8.5cm]{fczfc}
\caption{\label{fczfc} Temperature dependece of susceptibility following zero-field cooling (zfc) and field cooling (fc) cycles, for both single crystal and powdered samples for three representative Co concentrations corresponding to $x = 0.045$ (underdoped), $x = 0.061$ (optimally doped) and $x = 0.109$ (overdoped). Measurements were made in a field of 50 Oe applied parallel to the $ab$-plane of the single crystals.}
\end{figure}

Susceptibility data were taken for representative Co concentrations on the underdoped and overdoped side of the phase diagram, as well as for an optimally doped sample. These measurements were made for applied fields of 50 Oe, oriented in the $ab$-plane, for both zero-field cooling and field cooling cycles. The measurements were then repeated for each sample after grinding to a fine powder to examine the role of screening currents. The results of these measurements are shown in Fig \ref{fczfc}. First we consider the single crystal data (black data points). Field cooled values are relatively small for each Co concentration, indicative of substantial flux trapping. However, zero-field cooled values are large for both underdoped, optimally doped and overdoped cases indicating the presence of significant screening currents. These values do not change substantially when the single crystals are ground to fine powders (red data points) for the optimally doped and overdoped cases, but significantly the underdoped material has a substantially smaller value for the zfc susceptibility. It is clear that the underdoped material is still a bulk superconductor. However, the ability of the underdoped material to support screening currents is significantly reduced by powdering. This effect might be related to the carrier concentration. Specifically, the superfluid density will presumably be smaller for the underdoped material given the coexistence with the SDW state, leading to larger values for the penetration depth. In this case, it is possible that some fraction of the grains in the finely powdered sample will have a size smaller than the penetration depth, reducing the zfc susceptibility. Equally, it is possible that powdering causes significant damage to the surface of the resulting grains. The asymmetry of the superconducting dome might then lead to a stronger suppression of superconductivity for underdoped relative to overdoped samples. Most intriguingly, the lower zfc susceptibility might also reflect the presence of intrinsic inhomogeneity on the underdoped side of the phase diagram, because the percolating shielding currents do not
encompass the entire sample. The current measurements are not sufficient to determine the origin of this rather dramatic effect on the susceptibility. Nevertheless, it is clear that there is a substantial volume fraction even in the underdoped region of the phase diagram, but that the superconductivity is to some extent more fragile, at least to the effect of powdering.

The phase diagram shown in Fig. \ref{Phase} is highly suggestive of a scenario in which the superconductivity is intimately related to the presence of a quantum critical point associated with the eventual suppression of the structural and/or magnetic phase transitions. Clearly this needs to be investigated in greater detail. Nevertheless, this observation raises the question of whether a similar mechanism is at work in the hole doped analog Ba$_{1-x}$K$_x$Fe$_2$As$_2$, and more broadly for all of the superconducting Fe-pnictides. In the specific case of BaFe$_2$As$_2$, Co-doping is clearly more effective at suppressing the magnetic/structural transitions than K-doping, but it remains to be determined whether this is related to an asymmetry in the electronic structure or to other possibilities.

During preparation of this manuscript we became aware of two similar studies of the effects of Co-substitution in BaFe$_2$As$_2$ recently posted on the ArXiv.org preprint server. In the first of these, Ni \textit{et. al.} map out a similar phase diagram, also finding a splitting of the magnetic/structural phase transition \cite{canfield-2}. These results broadly corroborate our findings. In the second, Ning \textit{et. al.} use NMR measurements to follow spin dynamics for representative underdoped and overdoped samples, also surmising the presence of a quantum critical point \cite{mandrus-2}.

\section{Conclusion}
In conclusion, we have determined the phase diagram for the electron-doped superconductor Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$. We find that the single structural/magnetic phase transition that occurs in BaFe$_2$As$_2$ splits with Co-doping, although the nature of the two split transitions remains to be determined. Critical temperatures associated with both transitions are progressively reduced with increasing Co concentration, and are suppressed below $T_c$ coincident with optimal doping. The superconducting "dome" extends over a limited range of Co concentrations, from $x$ = 0.025 to $x$ = 0.18. These observations clearly delineate regions of the $x-T$ phase diagram for which further experiments have the potential to reveal the interplay between the structural/magnetic and superconducting phase transitions.

\section{Acknowledgements}
The authors thank R. E. Jones for technical assistance with EMPA measurements, and T. H. Geballe, M. R. Beasley and L. DeGiorgi for helpful discussions. This work is supported by the DOE, Office of Basic Energy Sciences, under contract no. DE-AC02-76SF00515.

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