The Crab Nebula's Composition and Precursor Star Mass

Abstract: We present results of new photoionization calculations for investigating gaseous regions that represent potentially expected stages of nuclear processing in the Crab Nebula supernova remnant. In addition to gas resulting from CNO-processing and oxygen-burning, as previously reported, a large component of the nebula appears to be carbon-rich. These results suggest that the precursor star had an initial mass of 9.5 solar masses or more.

\section{INTRODUCTION}

The Crab Nebula (M1~=~NGC1952) is generally recognized as the remnant
of the core-collapse supernova SN1054. Estimates of the precursor star's initial mass have been
in the range $8-12$~M$_{\odot}$
(see Davidson \& Fesen 1985 and references therein).  The visible remnant contains at least $1-2$~M$_{\odot}$
of He-rich line-emitting gas (MacAlpine \& Uomoto 1991)\footnote{Fesen et al. 1997 used global line photometry
to suggest a higher nebular mass, approaching that once postulated by MacAlpine et al. 1989.  However, as
discussed in MacAlpine \& Uomoto 1991, large spatial variations in the He/H ratio cause spuriously large mass
estimates based on global line photometry.  More accurate emitting gas estimates should take into account
spatially resolved line photometry in order to allow for the effects of helium abundance variations on
N(S)/N(H), which is used along with $[$S~II$]$ line intensities to estimate the amount of neutral gas present.},
the neutron star probably has about 1.4~M$_{\odot}$,
and the outer layers of the precursor may have left in a pre-supernova wind (Nomoto et al. 1982)
or a shockwave (e.g., Chevalier 1977; Davidson \& Fesen 1985).

Knowledge of the chemical composition of the remnant is vitally important for a better understanding
of the precursor star's initial mass and associated details of the core-collapse event.  According to
Nomoto et al. (1982), an $8-9.5$~M$_{\odot}$ star would likely have undergone an O-Ne-Mg core collapse
following electron capture, whereas a star more massive than 9.5~M$_{\odot}$ would have developed an Fe-rich
core.  A distinguishing factor between these two possible scenarios would be the nebular
carbon abundance.  As pointed out by Nomoto (1985), for a precursor with initial mass less than about 9.5~M$_{\odot}$,
high N and low C mass fractions (compared with solar) would exist from CNO-cycle processing in the He-rich gas.
On the other hand, for precursor mass $>$~9.5~M$_{\odot}$, helium-burning and nitrogen-processing via
$^{14}$N($\alpha$,$\gamma$)$^{18}$F($\beta$$^{-}$$\nu$)$^{18}$O($\alpha$,$\gamma$)$^{22}$Ne ought to have
produced high C and low N.  In addition, models for a roughly 10~M$_{\odot}$ precursor (Woosley et al. 1980;
Woosley \& Weaver 1986) have suggested the existence of off-center neon and oxygen flashes just prior
to the core collapse, whereas such events may not be expected in lower mass stars (Hillebrandt et al. 1984;
Nomoto 1984).

The necessary chemical composition information should be obtainable.  Because of its young age
and location roughly 180~pc from
the Galactic plane, the Crab Nebula is relatively uncontaminated by interstellar material.
In addition, its line-emitting gas is ionized and heated primarily by locally-generated synchrotron
radiation (rather than shock heating), so the ionization, thermal, and
chemical characteristics of the gas can be accurately analyzed using numerical photoionization codes.

Although there has been a long history of observational and theoretical investigations aimed at understanding
abundances in the Crab Nebula, to some extent the results have been inconclusive or contradictory.
As summarized by Henry (1986), Henry and MacAlpine (1982) and Pequignot and Dennefeld (1983)
compared measured line-intensity data available at the time with photoionization calculations.
In both studies, the deduced C and N mass fractions were roughly solar or below, but there were problems
with interpretation of near-infrared $[$C~I$]$ emission lines (see below).
Davidson et al. (1982) and Blair et al. (1992) also obtained satellite observations of the ultraviolet
C~IV and $[$C~III$]$ lines for certain locations, and they concluded no definitive evidence for a carbon excess,
although the latter study left open the possibility.  More recently, using extensive
long-slit spectroscopy covering much of the nebula, MacAlpine et al. (1996) deduced high nitrogen
mass fractions in some filamants and also high abundances for products of oxygen-burning,
such as sulfur and argon, at other locations.  All of the above results, taken together, have generally been
interpreted as suggestive of a precursor in the $8-9.5$~M$_{\odot}$ range, although high sulfur and argon abundances
would appear to argue for a more massive star.

In this paper, we report on the results of new photoionization computations (summarized by MacAlpine et al. 2007b),
which were developed using new line-intensity measurements
which extend to 1~\micron\ (MacAlpine et al. 2007a, hereinafter Paper~1). These calculations suggest that
strong $[$N~II$]$~$\lambda$$\lambda$6548,6583 emission lines measured at many locations by MacAlpine et al.
(1996, 2007a) do {\it not} necessarily imply high nitrogen abundances, while strong measured $[$C~I$]$~$\lambda$9850
emission {\it can} be interpreted as indicative of high carbon mass fractions.  In addition, high mass fractions
for S and Ar in some locations are confirmed.  These results argue {\it consistently} for a precursor star
initial mass $\gtrsim$~9.5~M$_{\odot}$.

\section{PHOTOIONIZATION COMPUTATIONS AND NUCLEAR PROCESSING REGIONS}

Photoionization models were produced using the numerical code Cloudy, version 06.02
(described in Ferland, G. J. et al. 1998).  Guided by line-intensity measurements and correlations in Paper~1,
we developed plane-parallel computations for a range of physical conditions and chemical abundances
in efforts to investigate gaseous regions that may represent expected stages of nuclear
processing for relevant stellar models (e.g., Woosley \& Weaver 1995).

At least three nuclear processing domains may be inferred from the observations in Paper~1.
Figure~1 of this paper (reproduced from Figure~4 of Paper~1) shows the correlation between
measured, H$\alpha$-normalized, reddening-corrected $[$N~II$]$$\lambda$6583 and $[$S~II$]$$\lambda$6731
line intensities for 180 widely-distributed nebular locations.  It was suggested in Paper~1
that the strongest nitrogen emission may represent gas which has progressed
no further than the CNO-cycle, whereas the apparent majority of measurements (e.g., with normalized $[$N~II$]$$\lambda$6583
over the range 1 to 2.5) could represent material in which some helium-burning has taken place and nitrogen has been
converted into neon.  The strongest sulfur emission was postulated to represent more advanced processing
through oxygen-burning. In this paper, these regions of Figure~1 will be called nuclear processing ``Domains 1, 2, and 3,''
respectively.

\subsection{DOMAIN~2 (HELIUM-BURNING AND NITROGEN DEPLETION)}

We begin by discussing details of our analysis for Domain~2 because it involves the most prevalent nebular gas.
In addition, the relevant photoionization models also help to understand the other domains.

First, we developed a representative ``median observed spectrum'' for Domain~2.  This involved a subset of 19 spectra
for which we have both optical and near-infrared measurements (see Paper~1) and for which
the H$\alpha$-normalized $[$N~II$]$$\lambda$6583 emission lies in the range 1--2.5.  For each emission-line measurement,
we identified a median value and combined those into a ``median spectrum,'' as given in Table~1.
Then we adjusted photoionization model input parameters and chemical abundances in order to obtain
a satisfactory fit to the observed data.

Following considerable experimentation, all models presented here employ the ionizing synchrotron radiation
spectrum given by Davidson \& Fesen (1985), with added points across the unobserved ultraviolet range
as suggested in the documentation which accompanies the Cloudy code.
In addition, we employed a constant hydrogen density of 3000~cm$^{-3}$
and log ionization parameter~=~-3.2. The ionization parameter is defined as
U~=~$\Phi$/cN, where $\Phi$ is the flux of photons more energetic than 1~Rydberg striking the face of a cloud,
c is the speed of light, and N is the total density of hydrogen.  The ionizing spectrum did a good job of
producing typically observed HeII/HeI line ratios (e.g., MacAlpine et al. 1989), the adopted
hydrogen density led to calculated electron densities in S$^+$ zones comparable with those measured
by Fesen \& Kirshner (1982), and the log~U was necessary to obtain a consistent overall
emission line spectrum that matches the observations, in particular the
$[$S~III$]$$\lambda$9069/$[$S~II$]$$\lambda$6731 ratio.

The ionization parameter employed here is almost a factor of 10 lower than that which is necessary
to account for far-ultraviolet C~IV and $[$C~III$]$ lines (see Blair~et~al.~1992).  Henry and MacAlpine (1982)
noted similar ionization parameter differences when comparing photoionization calculations that
could match optical and ultraviolet observations separately. We do not consider this to be a problem.
As shown convincingly by Blair et al., much of the
ultraviolet line emission probably comes from lower-density, more diffuse gas, compared
with the gas in optical line-emitting filaments being analyzed here.

For our Domain~2 model, the helium abundance was set at 89\%
by mass fraction, consistent with some stellar model expectations (Woosley \& Weaver 1995)
for a region with partial helium-burning and also with numerous observations of He~I$\lambda$5876 emission
(e.g., Uomoto \& MacAlpine 1987; MacAlpine et al. 1989; Paper~1).

Figure~2 presents the computed Domain~2 hydrogen and helium ionization fractions as a function of depth into
an emitting gas cloud, with ionizing radiation entering from the left.  We note that the He$^{+2}$
and He$^{+}$ zones do not extend far into the cloud.  Because of the high helium abundance,
ionizing photons more energetic than 54.4~ev (the ionization potential of He$^{+}$)
and 24.6~ev (the ionization potential of He$^{0}$) are absorbed before they can penetrate
a substantial distance.  The hydrogen ionized depth is about 10$^{16}$~cm, which is consistent
with observations for Crab Nebula filaments (see Davidson \& Fesen 1985).

Figure~3 shows the computed electron temperature and density, plotted as a function of distance into the cloud.
The electron density is roughly 8000~cm$^{-3}$ where helium is ionized, drops to 3000~cm$^{-3}$ through
the H$^{+}$ zone, and then falls to 400-800~cm$^{-3}$.  The temperature peaks near
1.2$\times$10$^{4}$~K, is about 6-8$\times$10$^{3}$~K through the H$^{+}$ zone, and then remains
above about 6000~K to the depth plotted.

Figure~4 gives the modeled carbon and nitrogen ionization.  Because the C$^{+}$
ionization potential is comparable with that for He$^{+}$, and the C$^{0}$ ionization potential
is below that for H$^{0}$, the dominant ionization stage for carbon is C$^{+}$,
all the way from the He$^{+}$ zone to where the calculations were stopped at a depth of approximately seven times
the hydrogen ionization zone thickness.
In addition, because the N$^{+}$ and N$^{0}$ ionization potentials are slightly above (within 5~ev of)
those for He$^{0}$ and H$^{0}$, respectively, N$^{+}$ occupies an extensive, very proficient
``$[$N~II$]$ emitting zone,'' in which the temperature and electron density effectively support collisional excitation.

Figure~5 presents an expanded view of the lower 10\% of the vertical scale from Figure~4.
This illustrates how a relatively small, but significant amount of C$^{0}$ exists in the zones
where helium and hydrogen are ionized.  Progressing from left to right in the figure,
there is a small peak in the C$^{0}$ fraction where the electron density is highest due to ionized helium.
Then there is a gradual C$^{0}$ increase as the temperature declines, and finally there is essentially
no C$^{0}$ beyond the hydrogen ionization edge where the electron density becomes very low.
These distribution characteristics for C$^{0}$ are consistent what would be expected from
C$^{+}$$\rightarrow$C$^{0}$ recombination.

Figure~6 shows selected, computed line emissivities for this model, as a function of depth into the cloud.
First, we note that $[$N~II$]$~$\lambda$6584 is ``off the chart'' because of electron collisional
excitation in the warm, high-electron-density $[$N~II$]$ emitting zone.
H~I~$\lambda$6563 illustrates what would be expected from recombination, and $[$C~I$]$~$\lambda$9850
shows a similar distribution.  As discussed above, the latter results from C$^{+}$$\rightarrow$C$^{0}$ recombination
followed by electron collisional excitation.  The $[$O~I$]$~$\lambda$6300 line arises from collisional
excitation of O$^{0}$, and it is primarily important in the neutral hydrogen region because of
charge exchange processes between hydrogen and oxygen.

Previously, high nitrogen mass fractions for many locations were deduced from observations of strong
$[$N~II$]$~$\lambda$$\lambda$6548,6583 emission and models using the Cloudy code in MacAlpine et al. (1996).
However, those models did not adequately take into account very high helium abundance
and its effect on broadening the $[$N~II$]$ emitting zone.
The helium mass fraction over most of the nebula ranges from about 75\% to greater than 90\% (Uomoto \& MacAlpine 1987).
For such helium-rich regions, strong $[$N~II$]$ lines
would {\it not} necessarily imply nitrogen mass fractions greater than solar.

It has been known for decades that $[$C~I$]$~$\lambda$9850 emission is very strong in parts
of the Crab Nebula (Dennefeld \& Andrillat 1981; Henry et al. 1984).
Pequignot \& Dennefeld (1983) tried to match the observed strengths of this line by postulating
a high dielectronic recombination rate, and Henry et al. (1984) postulated that collisional excitation
by H$^{0}$ could play an important role.
Neither group was successful in using $[$C~I$]$~$\lambda$9850 for reliable carbon abundance estimates.
Therefore, difficult measurements
of the ultraviolet $[$C~III$]$ and C~IV emission were attempted as noted previously.
However, production of those lines may not be well understood either, in the sense that model
predictions can be wrong by as much as a factor of two depending on how He~II~L$\alpha$ absorption
is treated (Eastman et al. 1985).  In the present calculations, we understand the production of
$[$C~I$]$~$\lambda$9850 emission as being due to thermal electron collisional excitation.
This $^{1}$D$_{2}$$\rightarrow$$^{3}$P$_{2}$ transition is an analog of the $[$O~III$]$~$\lambda$5007 line.
Making use of effective collision strengths
given by Pequignot \& Aldrovandi (1976), we manually calculated emissivities at various locations
and matched the model-computed emissivities exactly.

Predicted relative emission-line intensities from this model are shown in Table~1, where it may be seen
that there is a reasonable match with the Domain~2 median spectrum for element mass fractions
given in Table~2.  The deduced Domain~2 nitrogen mass fraction is less than half of its solar value.
We also note that the carbon mass fraction required to produce the observed median $[$C~I$]$$\lambda$9850 emission
in Domain~2 is roughly 6 times its solar value, and the deduced oxygen mass fraction is very close to solar.
A surprising result is the sulfur mass fraction necessary to acount for the observed $[$S~II$]$$\lambda$6731
and $[$S~III$]$$\lambda$9069 emission, which is less than a third of its solar value.

\subsection{DOMAIN~1 (STRONGEST NITROGEN EMISSION)}

The photoionization calculations for Domain~1 gas employed the same ionizing spectrum, hydrogen density,
and ionization parameter as discussed above; and the model characteristics are similar to those
in Figures~2-6.  The median measured spectrum and computed line intensities are shown in Table~1.
The measured, normalized $[$N~II$]$$\lambda$6583 values for Domain~1
range from 3.1 to 5.6, and the adopted median is 4.5.  Unfortunately, the subset
of near-infrared data with measurable $[$C~I$]$$\lambda$9850 did not include any nebular locations
with $[$N~II$]$$\lambda$6583~$>$~2.5 (see Figure~10 of Paper~1). Therefore, no $[$S~III$]$ or $[$C~I$]$
line values are given for this domain in Table~1.

As may be seen in Table~2, the helium mass fraction required to match the observed He~1$\lambda$5876 line
is again 89$\%$ for Domain~1, which is consistent with expectations from stellar models in a zone
containing largely CNO-processed material. However, we note that the inferred nitrogen mass fraction is only solar,
because $[$N~II$]$$\lambda$6583 arises so effectively in the $[$N~II$]$ emitting zone discussed above
and illustrated in Figures~4 and 6.  This implies some gas mixing or helium-burning.  The deduced
oxygen mass fraction is solar, while the data indicate that
sulfur is significantly less than solar, as in Domain~2.  The carbon mass fraction input
to the Domain~1 model was {\it assumed} to be 0.5~solar. A lower-than-solar carbon mass fraction
would be consistent for a region in which CNO-processing has taken place.

\subsection{DOMAIN~3 (STRONGEST SULFUR EMISSION)}

The measured spectrum and computed line intensities for Domain~3 are also shown in Table~1.
By combining optical data for locations with measured $[$S~II$]$$\lambda$6731 over the range 3-5.9
(see Figure~1) together with a subset of near-infared data, we derived a representative median
spectrum for this domain.  The best-fit photoionization model was similar to that for Domains~1 and 2,
except it was necessary in this case to lower log~U to -3.5 to match the
$[$S~III$]$$\lambda$9069/$[$S~II$]$$\lambda$6731 ratio.

As seen in Table~2, the adopted helium mass fraction to match observed median He~I~$\lambda$5876
is still 89\% for Domain~3, which would seem to imply gas mixing for regions
in which significant processing beyond helium-burning has taken place.  As expected, the nitrogen mass fraction
is now well below solar.  The deduced carbon and oxygen mass fractions are both about 10 times solar.
We note that the $[$O~I$]$~$\lambda$6300 line emissivity plotted in Figure~6 illustrates that it arises
in a different region from $[$C~I$]$~$\lambda$9850, so the inferred, correlated carbon and oxygen mass fractions
are not the result of high-optical-depth regions in the emitting gas.  This also supports the idea that the
$[$C~I$]$~$\lambda$9850:$[$O~I$]$~$\lambda$6300 line correlations in Figure~9 of Paper~1
result from carbon and oxygen abundances having been produced together by nuclear processing.

The sulfur mass fraction is roughly 4 times solar in this domain where it is hypothesized that
oxygen-burning has taken place.  Further support for oxygen-burning comes from the high inferred argon
mass fraction in this domain, as well as Figure~6 of Paper~1,
which shows a very tight correlation
between $[$Ar~III$]$ and $[$S~II$]$ emission, supporting the idea that sulfur and argon are produced
together by the same process.
Finally, a correlation shown in Figure~8 of Paper~1 is consistent with the hypothesis that iron-peak nuclei
like nickel are produced by ``alpha-rich freeze out'' (see Woosley \& Weaver 1995; Jordan et al. 2003)
in Domain~3 regions containing enhanced silicon-group elements from oxygen-burning.

\section {SUMMARY AND DISCUSSION}

The photoionization calculations presented here, compared with emission-line measurements
from Paper~1, are suggestive of three nuclear processing ``domains'' for gas in the Crab Nebula.

1. Gas with the strongest $[$N~II$]$ lines has roughly solar nitrogen and oxygen mass fractions, probably a carbon
mass fraction less than solar, and significantly less-than-solar sulfur and argon mass fractions.
This domain with 89\% helium may be largely CNO-cycle processed gas in which modest helium-burning and
nitrogen depletion have taken place.

2. The majority of the emitting nebular gas, with $[$N~II$]$$\lambda$6583/H$\alpha$ in the range 1--2.5, also has
$\sim$~89\% helium implied by its median observed spectrum. For this domain, the deduced nitrogen
mass fraction is depleted to less than half its solar value.  In addition, the carbon
mass fraction is about 6 times solar, while oxygen is solar, and sulfur and argon
have approximately a third to half of their solar values.  It is postulated that this domain
represents gas in which significant helium-burning has taken place.

3. The third domain examined here has the highest $[$S~II$]$$\lambda$6731 measured line intensities.
We find that the inferred helium mass fraction is still comparable with the other domains,
and the nitrogen mass fraction is about a fourth of solar.  The carbon and oxygen mass fractions
are roughly 10 times solar, while the sulfur and argon mass fractions are about 4-6 times solar.
This is considered to be a domain in which oxygen-burning has occurred and some mixing has taken place.

Gas resulting from the CNO-cycle and gas with products of oxygen-burning
were noted previously (MacAlpine et al. 1996), where the former was assumed to be the primary
component due to strong observed nitrogen emission throughout most of the nebula.
However, the current computations indicate that strong $[$N~II$]$$\lambda$$\lambda$6548,6583
emission does {\it not} necessarily imply high nitrogen abundance,
whereas strong $[$C~I$]$$\lambda$$\lambda$9823,9850 emission from much of the gas {\it does}
imply high carbon abundance.

The finding of high carbon mass fraction throughout much of the Crab Nebula gas is particularly
significant because it suggests that the supernova precursor star had an initial mass $\gtrsim$~9.5~M$_{\odot}$.
However, it is important to consider whether the near-infrared emission-line measurements
presented in Paper~1, and used here, may be biased toward high $[$C~I$]$$\lambda$9850 emission.
As discussed in Paper~1, the near-infrared data were obtained during non-optimal seeing, and
$[$C~I$]$ measurements were made only where potentially useful $[$S~III$]$ or $[$C~I$]$ line-emission
could be identified in the two-dimensional, near-infrared spectral images. To investigate the possibility
of consequent biasing, we compared the current measurements with $[$C~I$]$$\lambda$$\lambda$9823,9850
emission reported by Henry et al. (1984).
The ranges of measured $\lambda$9850 from that study and this one are comparable, with marginally
lower and higher values in the larger data set used here.  However, we still suspect that the current
median spectrum for Domain~2 may be biased toward stronger $[$C~I$]$ in view of how
the measured positions were selected.  That having been said, we note that the {\it primary} reason
for the current finding of high carbon mass fractions is a photoionization calculation whereby we can derive potentially
useful carbon abundances from $[$C~I{$]$ emission.  Furthermore, even though our $\lambda$9850 {\it median}
values may bias the inferred carbon mass fraction for a domain, the individual measurements
do imply that carbon-rich gas exists at numerous loctions in the nebula.  In any case, a more extensive
set of near-infrared spectra would be very useful for understanding the overall nebular carbon.

We should also reflect here on the fact that the photoionization code gave results for plane-parallel
geometrical situations, whereas line-emitting gas in the Crab Nebula may be immersed in ionizing
synchrotron radiation from all directions.  Therefore, we developed algorithms for converting the
plane-parallel results to what would be expected for convex cylindrical and spherical clouds.
Potentially noteworthy changes in H$\alpha$-normalized lines occurred for only $[$S~II$]$$\lambda$6731
and $[$O~I$]$$\lambda$6300 emission.  In going from plane-parallel to spherical geometry, the computed
$\lambda$6731/H$\alpha$ could be reduced by as much as 30\%, and $\lambda$6300/H$\alpha$
could be reduced by about 40\%.

Because adopted ionization parameters for the photoionization calculations
involved matching the observed $[$S~III$]$$\lambda$9069/$[$S~II$]$$\lambda$6731 ratio, an employed
ionization parameter could be slightly lower for a convex spherical geometry.  However,
we investigated this possibility and found that it would probably not
have a significant effect on the conclusions presented here.  We also considered whether
convex spherical or cylindrical clouds might affect our deduced low sulfur mass fractions for Domains~1 and 2.
By using a spherical gometry, we could increase a sulfur mass fraction by perhaps 30\%,
but it (and that for argon) would still be low in these domains.  Finally, we note that Henry \& MacAlpine (1982)
also reported lower-than-solar sulfur for particular filaments.  We believe the low sulfur mass fractions
in Domains~1 and 2 to be valid results of this investigation.

\acknowledgments

This work was supported by Trinity University and by an endowment from the late
Mr. Gilbert Denman.


Low-temperature specific heat for ferromagnetic and antiferromagnetic orders in CaRu1-xMnxO3

Abstract: Low-temperature specific heat of CaRu1-xMnxO3 was measured to clarify the role of d electrons in ferromagnetic and antiferromagnetic orders observed above x=0.2. Specific heat divided by temperature C_p/T is found to roughly follow a T^2 function, and relatively large magnitudes of electronic specific heat coefficient gamma were obtained in wide x range. In particular, gamma is unchanged from the value at x=0 (84 mJ/K^2 mol) in the paramagnetic state for x<=0.1, but linearly reduced with increasing x above x= 0.2. These features of gamma strongly suggest that itinerant d electrons are tightly coupled with the evolution of magnetic orders in small and intermediate Mn concentrations.

\section{Introduction}
The relationship between anomalous charge transport and magnetism in the vicinity of metal-insulator transition is a longstanding subject in the physics of strongly correlated electron systems. The Ru-based metal CaRuO$_3$ is considered to have a paramagnetic ground state with itinerant characteristics of Ru 4d electrons \cite{rf:Cao1997,rf:Kiyama1998}. In CaRuO$_3$, magnetization is found to be significantly enhanced at low temperatures \cite{rf:Kiyama1999,rf:Felner2000}, while no long-range magnetic order is observed in microscopic measurements such as neutron scattering and NMR \cite{rf:Martinez1995,rf:Mukuda1999,rf:Yoshimura1999}. In addition, electrical resistivity $\rho$ shows an unusually large magnitude at high temperatures, corresponding to a very small mean-free pass of electron conduction comparable to lattice constants \cite{rf:Cao2008}. These features imply that CaRuO$_3$ involves instabilities for both magnetism and charge transport, which are enhanced by strong electronic correlation between itinerant Ru 4d electrons.

In the mixed compound CaRu$_{1-x}$Mn$_x$O$_3$ \cite{rf:Sugiyama1999,rf:Shames2004,rf:Markovich2006,rf:Taniguchi2008,rf:Mizusaki2009,rf:Kawanaka2009}, it is revealed that a ferromagnetic (FM) order evolves above $x=0.2$, and becomes stable in the intermediate $x$ range. Both spontaneous magnetization $M_0$ and Curie temperature $T_{\rm C}$ show maximum values at $x\sim 0.7$ ($M_0\sim 1\mu_{\rm B}/{\rm f.u.}$ and $T_{\rm C}=170\ {\rm K}$), followed by the reductions of them with further increasing $x$. Furthermore, an antiferromagnetic (AFM) spin arrangement with the G-type structure is observed in the FM ordered phase. Finally, at $x=1$ the FM order vanishes, and only the G-type AFM order appears below $T_{\rm N}=120\ {\rm K}$ \cite{rf:Wollan1955}. By doping Mn, the feature of $\rho$ also changes from metallic to less conductive above $x \sim 0.1$, but the magnitude of $\rho$ still remains in small values in the wide $x$ range, that is, smaller than $0.1\ \Omega$ cm even at low temperatures for $x \le 0.8$. At $x=1$, in contrast, $\rho$ dramatically increases with decreasing temperature, whose magnitude at low temperatures exceeds $10^6$ $\Omega\ {\rm cm}$ \cite{rf:Neumeier2000}. These properties indicate that doping Mn into CaRuO$_3$ modifies the itinerant electronic states due to strong correlations between Ru 4d and Mn 3d electrons, and it significantly affects the magnetic and transport properties. To clarify the variation of the d electronic states, we have performed the low-temperature specific-heat measurements for CaRu$_{1-x}$Mn$_x$O$_3$.

\section{Experiment Details}
Polycrystalline samples of CaRu$_{1-x}$Mn$_x$O$_3$ for $0 \le x \le 1$ were prepared by means of the conventional solid-state method. The mixture of appropriate amounts of CaCO$_3$, RuO$_2$ and MnO were first calcined at 800 $^\circ$C for 24 hours. They were shaped into pellets after careful mixing, and then sintered at 1200  $^\circ$C for 48 hours. By means of powder X-ray diffraction technique we checked that all the samples have a distorted perovskite structure (the GaFeO$_3$-type orthorhombic structure) \cite{rf:Sugiyama1999,rf:Shames2004,rf:Markovich2006,rf:Taniguchi2008,rf:Mizusaki2009,rf:Kawanaka2009} without any extrinsic phase. Specific heat $C_p$ was measured between 3 K and 275 K with a thermal-relaxation method using both a commercial system (PPMS: Quantum Design) and a hand-made equipment, where we used plate-shaped samples with the mass of about 4 mg and 40 mg, respectively. The thermal-relaxation curves for all the measurements were well definitive, and the $C_p$ data obtained from both equipments were consistent within the experimental accuracy. Ac-susceptibility was measured in the temperature range of 5-300 K to estimate the FM and AFM transition temperatures. Frequency and amplitude of the applied ac field were 180 Hz and $\sim 0.5\ {\rm Oe}$, respectively.

\section{Results and Discussion}
\begin{figure}[tbp]
\begin{center}
\includegraphics[keepaspectratio,width=0.75\textwidth]{yokofig1.eps}
\end{center}
\vspace{-15pt}
  \caption{
Temperature variations of (a) specific heat and (b) real part of ac susceptibility for CaRu$_{1-x}$Mn$_x$O$_3$. In (a), closed and open arrows indicate the onsets of FM and AFM orders, respectively, estimated from the peak position of the ac susceptibility. The types of magnetic order are determined from the previous reports on magnetization and neutron scattering measurements \cite{rf:Shames2004,rf:Markovich2006,rf:Taniguchi2008,rf:Mizusaki2009,rf:Kawanaka2009}. The inset of (a) shows an enlargement of $C_p$ around the AFM transition temperatures for $x=0.9$ and 1.
}
\end{figure}
Figure 1(a) shows temperature variations of the specific heat $C_p(T)$ for CaRu$_{1-x}$Mn$_x$O$_3$. Clear jumps ascribed to the AFM transition are observed at 122.7 K for $x=1$ and 125.1 K for $x=0.9$ (see inset of Fig. 1(a)), but no distinct anomaly is seen in the $C_p(T)$ curves for $x \le 0.8$. This contrasts with temperature dependence of the real part of ac susceptibility $\chi'(T)$ for $0.2 \le x \le 1$, where the FM and AFM transitions are indicated by the appearance of peaks (Fig.\ 1(b)). The transition temperatures estimated from $\chi'(T)$ are consistent with those determined by the dc magnetization and neutron scattering measurements \cite{rf:Shames2004,rf:Markovich2006,rf:Taniguchi2008,rf:Mizusaki2009,rf:Kawanaka2009}. The peaks seen in $\chi'(T)$ are somewhat broaden probably due to the distribution of the transition temperatures in the sample generated by doping. We consider this effect to be one of the reasons for the reduction of jumps in $C_p(T)$. On the other hand, it is interesting that similar reduction of anomaly is observed in the $C_p(T)$ data for Mn-doped SrRuO$_3$ \cite{rf:Yamaji2010}. SrRuO$_3$ is a ferromagnet with itinerant characteristics of Ru 4d electrons, and $C_p(T)$ exhibits a jump at $T_{\rm C}=162\ {\rm K}$. In the mixed compound SrRu$_{1-x}$Mn$_x$O$_3$, the anomaly at $T_{\rm C}$ in $C_p(T)$ is rapidly suppressed in the FM metallic phase for $x \le 0.3$, but a shoulder-like anomaly again appears in $C_p(T)$ at the AFM transition temperatures in the strongly insulating region for $x\ge 0.4$. As for CaRu$_{1-x}$Mn$_x$O$_3$, the electrical resistivity is much reduced for $x \le 0.8$. The disappearance of anomaly in $C_p(T)$ may thus be related to the small entropy change associated with the magnetic transitions in the small-resistivity region, as well as the broadening of the onset due to the inhomogeneity.

\begin{figure}[tbp]
\begin{center}
\includegraphics[keepaspectratio,width=0.75\textwidth]{yokofig2.eps}
\end{center}
\vspace{-10pt}
  \caption{
(a) Low-temperature specific heat divided by temperature plotted as a function of $T^2$, (b) electronic specific-heat coefficient $\gamma$ and lattice contribution on specific heat $\beta$, and (c) $T_{\rm C}$ and $T_{\rm N}$ estimated from the $\chi'(T)$ data for CaRu$_{1-x}$Mn$_x$O$_3$. The definition of parameters $\gamma$ and $\beta$ is described in the text. The solid lines in (c) are guides to the eye \cite{rf:Shames2004,rf:Markovich2006}.
}
\end{figure}
Displayed in Fig. 2(a) is low-temperature specific heat divided by temperature $C_p/T$ plotted as a function of $T^2$. $C_p/T$ for entire $x$ range are roughly proportional to $T^2$. The finite $C_p/T$ values for $T\to 0$ except at $x=1$ indicate that the Fermi-liquid excitation as well as phonon are mainly responsible for $C_p(T)$. However, $C_p/T$ in the small and intermediate $x$ range shows a very weak deviation from the $T^2$ function with a downward curvature below $T^2\sim 100\ {\rm K}^2$ ($T\sim 10\ {\rm K}$). It may be ascribed to the spin fluctuation generated in the magnetically ordered phase. The $C_p/T$ data are fitted with a function of $\gamma+\beta T^2$ to derive the electronic specific heat coefficient $\gamma$. In Fig. 2(b), we show $x$ variations of $\gamma$ and $\beta$ obtained from the best fit in the temperature range between 3 K and 10 K. For comparison, we also plot the $x-T$ phase diagram obtained from the $\chi'(T)$ data in Fig. 2(c). The $\gamma$ value for $x=0$ is estimated to be 84 mJ/K$^2$ mol, which is consistent with that reported previously \cite{rf:Kikugawa2009}. The magnitude of $\gamma$ does not change in the paramagnetic phase for $x \le 0.1$, but linearly decreases with increasing $x$ in the magnetically ordered region ($x\ge 0.2$). It is highly reduced and becomes $\sim 5\ {\rm mJ/K^2\ mol}$ or less for $x\ge 0.8$, where $T_{\rm C}$ decreases through the maximum at $x\sim 0.7$. On the other hand, $\beta$ is estimated to be 0.19 mJ/K$^4$ mol for $x=0$ that yields the Debye temperature of $\Theta_{\rm D}=370\ {\rm K}$. It is reduced in Mn-rich region and finally becomes 0.085 mJ/K$^4$ mol ($\Theta_{\rm D}=490\ {\rm K}$) at $x=1$. This feature seems to be consistent with the volume contraction in Mn-rich concentrations found by X-ray diffraction measurements \cite{rf:Shames2004,rf:Taniguchi2008}.

In CaRu$_{1-x}$Mn$_x$O$_3$, it is proposed that a double-exchange mechanism between Mn ions plays a crucial role in the evolution of the FM order in Mn-rich region, and in lower $x$ range an interplay of double-exchange and superexchange interactions between d electrons located at Ru and Mn ions may yield the formation of inhomogeneous FM and AFM orders \cite{rf:Shames2004,rf:Markovich2006}. On the other hand, the present $C_p$ investigation revealed that $\gamma$ shows relatively large values in small and intermediate $x$ ranges, while it is strongly reduced in Mn-rich region. In particular, $\gamma$ at $x=0.4$ is estimated to be 38 mJ/K$^2$ mol, which is almost the same magnitude as that for the ferromagnetic metal SrRuO$_3$ (33 mJ/K$^2$ mol). The large $\gamma$ values indicate that d electrons have itinerant characteristics in this $x$ range. Furthermore, we observed that $\gamma$ is nearly constant in the paramagnetic region, and it is then reduced when the FM order occurs above $x=0.2$. Such a variation of $\gamma$ cannot be simply attributed to an increase in the fraction of d electrons localized at Mn ions in the sample, but strongly suggests that the itinerant d electrons originating from both Ru and Mn are tightly coupled with the magnetic orders. This is supported by our recent magnetization and neutron-diffraction measurements for $x=0.4$, where the magnetization shows typical features on fluctuation expected from itinerant electron spins, and the neutron-diffraction experiments revealed the occurrence of the long-range magnetic orders \cite{rf:Kawanaka2009}.

\section{Summary}
We performed $C_p$ measurements for CaRu$_{1-x}$Mn$_x$O$_3$ ($0 \le x \le 1$). Low-temperature $\gamma$ value is estimated to be 84 mJ/K$^2$ mol for pure CaRuO$_3$, and it is unchanged by doping Mn in the paramagnetic region for $x \le 0.1$. In the magnetically ordered phase, $\gamma$ is linearly reduced with increasing $x$ for $0.2 \le x \le 0.6$, and then approaches $\sim 5$ mJ/K$^2$ mol or less for $x \ge 0.8$. The relatively large magnitude of $\gamma$ in the small and intermediate Mn concentrations suggests that d electrons have itinerant characteristics contrary to localized ones expected in Mn-rich concentrations, and they play important roles in the magnetic orderings.

\section*{Acknowledgment}
We thank Y. Nishihara and A. Yamaji for fruitful discussion. This work was carried out by the joint research in ISSP, the University of Tokyo, and partly supported by a Grant-in-Aid for Scientific Research on Innovative Areas ``Heavy Electrons" (No.23102703) from the Ministry of Education, Culture, Sports, Science and Technology of Japan.

\section*{References}

Strongly Correlated Electrons (cond-mat.str-el)


6 W, 1 kHz linewidth, tunable continuous-wave near-infrared laser

Abstract: A modified Coherent 899-21 titanium sapphire laser is injection locked to produce 6-6.5 W of single-frequency light at 852 nm. After closed-loop amplitude control and frequency stabilization to a high-finesse cavity, it delivers 4-4.5 W with <1 kHz linewidth at the output of a single-mode fiber. The laser is tunable from about 700-1000 nm; up to 8 W should be possible at 750-810 nm.

\section{Introduction}
The improvement of tunable, high-power, low linewidth
lasers has been instrumental in advancing atomic physics.
State-of-the-art titanium:sapphire (Ti:sa) lasers span wavelengths
between 700-1000\,nm or more. Typical commercial models use
intracavity etalons for single-mode operation and achieve an
output power around 1.5-2.5\,W. Injection lock
\cite{Tanaka,Cummings,Cha} allows for dispensing with the etalons and
thus reaching a higher power of up to 5\,W so far \cite{Cha2}.
Diode lasers with tapered amplifiers, which are only available in few spectral regions, currently reach 1-2\,W (up to
2.5\,W in $\mu$s pulses \cite{Takase}) at 780\,nm and 0.5\,W near
850\,nm. Because of their lower beam quality and the need for an
optical isolator, however, they rarely deliver more than 0.5\,W
into single-mode fibers. The linewidth of both commercial single
frequency Ti:sapphire and extended-cavity diode lasers is on the
order of a few 100\,kHz. A reduction to the kHz and even Hz level
can be achieved by stabilization to high-finesse optical cavities.
Here, we report an injection-locked Ti:sapphire laser system that
achieves up to 6.5\,W at 850\,nm, up to 4.5\,W after a
single-mode, polarization maintaining (SM-PM) fiber, and $<1\,$kHz
linewidth by stabilization to a high-finesse reference cavity.

\section{Description}
Figure 1 shows the setup. It can be roughly divided into the
injection lock system shown on the left, the frequency
stabilization in the upper right, and the amplitude control shown
in the lower right part of the figure.

\begin{figure}
\centering\epsfig{file=6Wschem.eps, width=0.5\textwidth}
\caption{Setup.}
\end{figure}

The master laser (not shown) is a Coherent 899 with 1.2-1.3\,W
output power when pumped with 10\,W from a Coherent Verdi V-10. It
is frequency stabilized using the internal piezoelectric
transducer (PZT) and Brewster plate as actuators. For this
purpose, we modified the control box as described in Ref.
\cite{Haubrich}. As a frequency reference, we use an
extended-cavity diode laser (ECDL) \cite{tracked}, which is in
turn stabilized to a $D_2$ transition in a cesium vapor cell.
A detuning of 0 to $\pm20$\,GHz can be set by a microwave
synthesizer.

The master passes an acousto-optic modulator, AOM 1 (Fig. 1),
driven by voltage controlled oscillator VCO 1 at about 120\,MHz.
The light is then coupled to an SM-PM fiber, in order to decouple
the alignments of the master and slave. About 0.8\,W are typical
after the fiber.

The slave laser is another Coherent 899. To increase the output
power, we removed the thick and the thin etalon as well as the
Brewster plate, as they are not required for injection locked
operation \cite{birfilter} and replaced the output coupler by one
having 10\% transmission \cite{coupler}. Pumped with 19\,W from a
Coherent Innova 400 argon-ion laser (not shown), we obtain 5.5\,W
output from the free-running slave.

For injection lock, the passive resonance frequency of the slave
laser's cavity must be stabilized to the master. For generating an
error signal without adding sidebands to the radiation, we use the
H\"ansch-Couillaud method \cite{Couillaud}. The master laser's
polarization has a small angle relative to the slave, set by a
half-wave retardation plate. The small angle means that most of
the electric field of the master can add to the slave's field (at
lower master power, an angle of 45$^\circ$ would be used to give a
larger error signal). To generate the error signal, a small
fraction (25\,mW) of the output beam is picked up and its
polarization analyzed.

The error signal is fed back to the piezoelectric transducer (PZT)
that translates the lower fold mirror M3 \cite{remark}. We found
the $\sim 5$\,kHz bandwidth of this path to be sufficient for
locking, but only just. Therefore, we add a fast (about 100\,kHz
bandwidth) feedback channel to VCO 1 (mini-circuits POS-150), in
analogy to \cite{Cummings,Cha}. This makes the lock reliable, at
the price of some broadening of the master laser's linewidth. With
up to 0.8\,W from the master laser, there was no need to optimize
for low master laser powers. We confirmed, however, that the slave
can still be locked with about 200\,mW from the master. Previous
work used master laser powers of 1\% \cite{Tanaka} to 1.5\%
\cite{Cummings,Cha} of the slave, which would be in the range of
50-100\,mW for our power level.

The combined power of master and slave (6\,W typical with a slope efficiency of $\sim 40$\%; 6.75\,W have
been achieved \cite{pumppower}) is deflected by AOM 2 (Crystal
Technology 3080-122) and coupled into a short \cite{SBS} SM-PM
fiber (OFR). In our application, the AOM is
usually driven by 100-$\mu$s pulses; nevertheless, we confirmed
that cw operation for at least a few h is feasible without
destroying the fiber. An adjustable OFR fiber port is used for the spatial mode-matching of the injection and the slave laser cavity: The mode matching is optimized so that the highest coupling efficiency of the combined light into the short fiber is obtained.

Since the injection lock leads to a fixed phase relationship
between the master and the slave, given suitable alignment, the
outputs of both lasers can be added coherently: we obtain 4\,W
typical (4.5\,W best) after the fiber; this means that one of the
AOM or the fiber is at least 81\% efficient. This efficiency is
similar to what we obtain with either the master or the slave
alone, confirming that the outputs of the two lasers do indeed add
coherently, so that the power of the combined lasers into the mode
of the fiber is close to the sum of the individual ones.

The acoustical delay of the AOM is reduced to about $0.5\,\mu$s by
moving the transducer close to the beam with a translation stage;
this adjustment must be made with care, as the very strong
Ti:sapphire beam instantly destroys the AOM if it ever hits the
transducer.

The linewidth of the locked slave is about 500\,kHz full width at
half maximum. To reduce it, we split off a small sample, shift it
by AOM 3 and stabilize it to a Fabry-Perot cavity with a linewidth
of about 25\,kHz (finesse 75,000). The error signal is generated
by the Pound-Drever-Hall method. Phase modulation at 16\,MHz is
applied by an electro-optical modulator (EOM). The light reflected
by the cavity is detected and down-converted to obtain the error
signal. We feed back to AOM3's frequency via VCO 2, with a loop
bandwidth of about 300\,kHz.  VCO 2 also controls the main beam
frequency by driving AOM 2. In effect, the main beam is stabilized
to the cavity independent of the intensity setting of AOM 2. AOM 3 is positioned with a translation stage to match the delays of AOM 2 and AOM 3 so that the noise in the beat note measurement (will be described in the following paragraphs) is minimized. A
slow (several 100\,Hz bandwidth) servo (not shown) zeros the dc
signal to the VCO by feeding back to the cavity via a PZT. This is
used to take out the cavity drift and means that on long time
scales, the slave is locked to the cesium spectroscopy.

To verify the linewidth of the system, a beat note measurement
with an independent low-linewidth reference laser would be ideal.
Such a laser was unfortunately not available to us. As the second
best option, we use the {\em transmitted} light of the cavity as a
reference. This method is as good as an independent laser for
noise frequencies much higher than the cavity linewidth: Such
noise components are filtered out by the cavity and are thus
suppressed in its transmission. Regardless of frequency, the
method measures the laser frequency relative to the cavity
resonance. To minimize fluctuations of the cavity resonance
itself, the cavity is held by viton rings inside a hermetically
sealed (but not evacuated) stainless steel chamber to provide
acoustic and thermal insulation.

For the beat note measurement, we overlap a sample of the power
after the fiber with the transmitted light of the cavity on a
photodetector. The measured spectrum (Fig. 2) reveals a linewidth
of 1.3\,kHz (full width at half maximum), limited by the 1\,kHz
resolution of our spectrum analyzer. Assuming that the linewidths
add geometrically, we conclude that the laser linewidth is
800\,Hz. Frequency modulation sidebands at 24\,kHz are an effect
of mechanical vibrations of the cavity. The wideband spectrum
(Fig. 2, right) shows an about 800\,kHz wide noise pedestal on the
order of -80\,dBc/Hz, the residual of the free-running laser
noise. At frequencies more than 400\,kHz from the carrier, the
noise drops rapidly.

\begin{figure}
\centering\includegraphics[width=7cm]{spectrum.eps} \caption{Beat
measurement between the fiber output and the transmission of the
cavity. Left: Linear scale of the innermost 80\,kHz to show
linewidth. 1\,kHz resolution. Right: 1\,MHz scan shown with
logarithmic scale in dB referred to carrier. 10\,kHz resolution
bandwidth.}
\end{figure}

The laser system as installed is tunable within 780-925\,nm, and
more than 700-1020nm with the short or long wavelength optics
sets. Between 750-810\,nm, the output power of Ti:sapphire lasers
such as the Coherent 899 is typically about 30\% higher than at
850\,nm, so about 8\,W are possible. We have, however, not
ascertained that fiber coupling of this is feasible. If wide
tunability is not required, a considerable cost reduction can be
achieved by using a diode laser (possibly with a tapered
amplifier) as the master laser.

\section{Conclusion}
In summary, we have presented an injection locked laser system.
The master and slave laser's power add coherently to up to
6.75\,W. The system produces power- (up to 4.5\,W) and frequency
stabilized (linewidth below 1\,kHz) and mode-filtered radiation at
the output of a single-mode, polarization maintaining fiber. The
optical frequency is locked to a Cs hyperfine transition with an
offset between 0 and $\pm 20$\,GHz; a high-finesse cavity is used
for linewidth reduction. The amplitude stabilization with its
250\,kHz loop bandwidth allows us to precisely shape light pulses
of arbitrary form. This system is ideally suitable for demanding
applications such as high-order Bragg diffraction \cite{BraggPRL}
in atom interferometry, or as a basis for a powerful, tunable uv
source by frequency multiplication \cite{Mes}.

\subsection{Acknowledgments} H.M. and S.H. thank the Alexander von Humboldt Foundation. This
material is based upon work supported by the National Science
Foundation under Grant No. 0400866 and by the Air Force Office of
Scientific Research under Award Number FA9550-04-1-0040.