\section{Introduction.}

The properties of confined spins in semiconductor nanostructures are
interesting because of their expected important role in quantum computing
and spintronic devices. In particular, the hole spin attracts much
attention because of its enhanced coherence time. In systems that
efficiently couple to optical
fields, like quantum wells and self-assembled quantum dots, an
interesting possibility is to study the spin dynamics using optical
spectroscopy tools. A method which is particularly suited for this
purpose is the time-resolved Kerr (or Faraday) rotation (TRKR, TRFR)
\cite{syperek07,kugler09,korn10}, where the evolution of the occupations of the hole
and trion Zeeman states is traced by investigating the rotation of the
polarization plane of the reflected or transmitted probe pulse.

In a recent work
\cite{machnikowski10}, we proposed a general description for the
dynamics of a confined hole-trion system in a tilted magnetic
field. We showed that a TRKR experiment provides rich
information about the rates of spin precession and decoherence (both
longitudinal and transverse). In particular, the optical response at
slightly tilted magnetic fields contains contributions related both to
the hole spin relaxation (longitudinal decoherence) with respect
to the hole spin quantization axis, as well as the dephasing of spin
coherences (transverse decoherence) with respect to this axis. Under
favorable experimental conditions (a sufficient separation of time scales), the
corresponding two decoherence rates can be deduced from a single run
of an experiment, although in a realistic system the latter may be
convoluted with (or even dominated by) the inhomogeneous dephasing.

Although the TRKR or TRFR method may be the most obvious choice for the
investigation of the spin dynamics the spin precession
of the trion and hole in a tilted field will lead to transitions between optically active
and inactive states which should be manifested in any form of the
optical response. In particular, signatures of the spin-related
dynamics should be visible in the four wave mixing (FWM) nonlinear
spectroscopy. As the FWM spectroscopy is one of the most widely used
methods in the investigation of semiconductors and their
nanostructures \cite{borri01,vagov04}
it may
be interesting to study how the spin precession and decoherence affect
the FWM response and what information on the spin-related kinetics can
be extracted from this kind of experiments.

In this contribution, we study the FWM response of the same system of
confined holes as discussed in our previous work
\cite{machnikowski10}. We show
that the FWM signal is also affected by the spin dynamics in a way
that allows one to extract the Larmor frequencies of electrons and
holes. The available information on the decoherence rates
is less specific than that contained in the TRKR response and is
available only as long as the hole spin dephasing times are not much longer
than the life time of the optical coherence. On the other hand, the
spins undergo partial refocusing in the two-pulse optical echo
experiment which leads to the appearance of components in the optical
response that depend on the inherent spin dephasing rates but are
insensitive to the inhomogeneous distribution of the Larmor
frequencies. In this way, the FWM method may be a valuable
experimental tool to study the spin dynamics, complementary
to the TRKR or TRFR techniques.

The paper is organized as follows. In Sec.~\ref{sec:model}, we define
the model under study. In Sec. \ref{sec:evol} we derive the FWM
response of the system. Next, in Sec.~\ref{sec:discusion}, we discuss
the result and compare the spin-related information contained in the
FWM signal to that available from the TRKR response.

\section{Model.}
\label{sec:model}

We study a system similar to that discussed in our previous work
\cite{machnikowski10}, consisting of an ensemble of quantum dots or
trapping centers in quantum wells with one confined hole in each of the
dots. The system is excited with a pump pulse at the time $t=-\tau$ and a
probe pulse at $t=0$. Unlike in our previous work, we are now
interested in the FWM response from the system. Experimentally, the
relevant third-order response is isolated by choosing the
appropriate excitation and detection directions (see
Ref.~\cite{machnikowski10} for a discussion). In the modeling, we
assume the pump and probe pulses to have phases $\phi_{1}$ and
$\phi_{2}$, respectively, and calculate the terms in the response of an
inhomogeneously broadened system that carry the phase
$2\phi_{2}-\phi_{1}$. Out of many possible configurations of the
polarizations of the excitation and detection, we choose the
circularly co-polarized one, with both pulses having the $\sigma_{+}$
polarizations and the detection being performed at the same
polarization.
The system is placed in a magnetic field
tilted by the angle $\vartheta$ to the normal, which defines the trion spin
quantization. The hole spin is quantized
along the axis at an angle $\varphi$ from the normal to the
sample, defined by the components of the highly anisotropic hole
Land\'e tensor. We will assume that the system geometry is very close
to the Voigt configuration, with the magnetic field only slightly
tilted from the system plane, so that $\cos\varphi\approx 0$ and
$\sin\varphi\approx 1$.

The system dynamics is a composition of many processes. For the optical
hole-trion transition, we assume the radiative recombination rate
$\gamma_{1}$ and the additional pure dephasing rate $\gamma_{0}$. The
hole and the trion spin precess in the magnetic field with the Larmor frequencies
$\wh$ and $\wt$, respectively. The longitudinal relaxation rate for
the hole spin is
$T_{1}^{-1}=\kappa_{+}+\kappa_{-}$, where $\kappa_{-},\kappa_{+}$ are the rates
for the spin flip transitions to the lower and to the upper Zeeman
level, related by the detailed balance condition
$\kappa_{-}/\kappa_{+}=\exp(\hbar\wh/k_{\mathrm{B}}T)$, where
$k_{\mathrm{B}}$ is the Boltzmann constant and $T$ is the
temperature. The hole spin transverse relaxation rate is
$T_{2}^{-1}=T_{1}^{-1}/2+\kappa_{0}$, where $\kappa_{0}$ is the
additional pure dephasing rate. We assume that the electron spin dephasing is slow
compared to the trion lifetime so that the electron spin coherence is
dominated by the latter.

The simulations presented in this contribution are performed for
$T=4$~K and for the magnetic field of 7~T tilted at an angle
$\pi/2-\vartheta=4^{\circ}$ from the system plane. For the in-plane and
perpendicular components of the hole g-tensor equal to $0.04$ and
$0.6$, respectively, this yields the hole Larmor frequency $\wh=0.036$~ps$^{-1}$ and
the hole spin quantization axis oriented at $\varphi=44^{\circ}$ from the
normal direction. The trion Larmor frequency is 0.16~ps$^{-1}$
(corresponding to the electron g-factor of 0.26). We will assume the
decoherence times $1/\gamma_{1}=1.2$~ns and $T_{1}=2$~ns and no
additional pure dephasing effects ($\gamma_{0}=\kappa_{0}=0$).

\section{The FWM response.}
\label{sec:evol}

Our theoretical analysis is based on the method developed in
Ref.~\cite{machnikowski10}: the system evolution is studied in the
density matrix formalism, with the optical and spin-related dephasing
included via a Lindblad dissipator in the evolution equation. In the
limit of coherent ultrafast optical driving, this approach yields
analytical formulas for the single system dynamics, which allows one to
perform averaging over an inhomogeneous distribution of various
parameters in the ensemble.

For a $\sigma_{+}$ probe, the only element of the density matrix
just after the probe pulse that carries the $e^{2i\phi_{2}}$ phase
dependence is
\begin{equation*}
\rho_{31}(0^{+})=\rho_{13}(0^{-})e^{2i\phi_{2}}
\sin^{2}\frac{\alpha_{2}}{2},
\end{equation*}
where $\phi_{2}$ is the phase of the probe pulse, $\alpha_{2}$ is its
area, and $0^{-}$ denotes the time instant just before the arrival of
the pulse. The evolution of the system state is then calculated using
the Master equation in the Lindblad form, like in Ref.~\cite{machnikowski10}. The
$\sigma_{+}$ interband coherence at a time $t>0$ is
\begin{equation*}
\rho_{31}(t)=\sum_{\pm,\pm}
d_{\pm,\pm}e^{\lambda_{\pm,\pm}t}\rho_{13}(0^{+})
e^{-iEt/\hbar},
\end{equation*}
where the amplitudes and the exponents are given by
\begin{eqnarray*}
d_{1}=\frac{\cos\varphi+ 1}{4}, & \quad &
\lambda_{1}=\frac{2i\wt -(2i\wh -\beta)}{4}, \\
d_{2}=\frac{\cos\varphi- 1}{4}, & \quad &
\lambda_{2}=\frac{2i\wt +(2i\wh -\beta)}{4}, \\
d_{3}=-\frac{\cos\varphi+ 1}{4}, & \quad &
\lambda_{3}=\frac{-2i\wt -(2i\wh -\beta)}{4}, \\
d_{4}=-\frac{\cos\varphi- 1}{4}, & \quad &
\lambda_{4}=\frac{-2i\wt +(2i\wh -\beta)}{4}. \\
\end{eqnarray*}
Here $E$ is the interband transition energy (at zero field) and
the dephasing constants are
\begin{equation*}
\Gamma=4\gamma_{0}+2\gamma_{1}+\kappa_{0}+\kappa_{+}+\kappa_{-}
\end{equation*}
and
\begin{equation*}
\beta=(\kappa_{-}-\kappa_{+})\left[
(\cos\varphi-1)\sin^{2}\varphi +1\right].
\end{equation*}

Now, we need to find $\rho_{13}(0^{-})$. For a $\sigma_{+}$--polarized
pump, just after the pump pulse (at $t=-\tau^{+}$),
the only non-zero element linear in the pump amplitude is
$\rho_{31}$. By
solving the Master equation one finds at $t=0^{-}$
\begin{equation}
\label{ro13}
\rho_{13}(0^{-})=\sum_{i}
c_{i}e^{\lambda_{i}^{*}\tau}e^{iE\tau/\hbar},
\end{equation}
where
\begin{equation*}
c_{i}=\rho_{13}(-\tau^{+})d_{i}.
\end{equation*}
Thus, at the time $t$ one finds
\begin{eqnarray*}
\rho_{31}(t) & = & \frac{1}{2}\sin\alpha_{1}\sin^{2}\frac{\alpha_{2}}{2}
e^{i(2\phi_{2}-\phi_{1})} \\
&&e^{-iE(t-\tau)/\hbar}
\sum_{ij} d_{i}d_{j}e^{\lambda_{i}t+\lambda_{j}^{*}\tau}.
\end{eqnarray*}
When averaging over the distribution of interband energies $E$ in a
typical QD ensemble, the exponent produces a narrow echo peak around
$t=\tau$. Since the spin-related evolution (frequencies $\lambda_{i}$)
is slow on this time scale, we can put $t=\tau$ under the
summation. Then the magnitude of the time-integrated response is
proportional to
\begin{eqnarray}
\lefteqn{\mathrm{FWM}  \sim  \left|
\sum_{i}d_{i}e^{\lambda_{i}\tau}\right|^{2}} \nonumber \\
&& =  \frac{1}{8}(1+\cos\omega_{\mathrm{t}}\tau)
e^{-\Gamma\tau/2}\Bigl[
(\cos\varphi-1)^{2}e^{-\beta\tau/2} \nonumber \\
&&\;\;\;\;+(\cos\varphi+1)^{2}e^{\beta\tau/2}
+2(1-\cos^{2}\varphi)\cos\omega_{\mathrm{h}}\tau \Bigr].
\label{fwm}
\end{eqnarray}
In order to take into account also the inhomogeneity of the Land\'e
tensors at different trapping centers this result should be averaged
over a certain distribution of the Larmor frequencies $\wh$ and $\wt$,
which we will assume to be Gaussian and characterized by the variances
$\sigma_{\mathrm{h}}^{2}$, $\sigma_{\mathrm{t}}^{2}$. We assume here
that the distribution of Land\'e tensors is uncorrelated to the
spectral positions of the trion transitions.

\begin{figure}[t]
\includegraphics*[width=\linewidth]{homo}
\caption{The FWM signal from a system without inhomogeneous broadening
  of the Larmor frequencies. (a) The
  time-integrated FWM signal as a function of the delay time; (b) The
  real part of the Fourier transform of the time-resolved signal.}
\label{fig:homo}
\end{figure}

In Fig.~\ref{fig:homo}(a), we show the calculated time-integrated FWM
signal in the absence of inhomogeneous distribution of the Larmor
frequencies. As follows from Eq.~\eqref{fwm}, oscillations at various
frequencies are present in the optical signal, corresponding to
combinations of the hole and trion frequencies. A more transparent
picture is obtained after performing a Fourier-transform of this
signal, as shown in Fig.~\ref{fig:homo}(b). Here one can see the zero-frequency
line, as well as lines at the frequencies $\pm\wt,\pm\wh$, and
$\pm\wt\pm\wh$.

\begin{figure}[t]
\includegraphics*[width=\linewidth]{inhomo}
\caption{The FWM signal from a system with inhomogeneously broadened
  Larmor frequencies. (a) The
  time-integrated FWM signal as a function of the delay time; (b) The
  real part of the Fourier transform of the time-resolved signal.}
\label{fig:inhomo}
\end{figure}

In a real system, the optical response is affected by the inhomogeneous
distribution of the relevant parameters. In the present case, the
Larmor frequencies are usually not identical for each hole-trion
system due to variations of the $g$-factors in the nanostructure. In
Fig.~\ref{fig:inhomo} we show the results of a simulation performed
for the same parameters as in the previous case but with the
additional effect of inhomogeneous broadening of the Larmor frequency,
which is assumed to be equal to 20\% of their average values. Due to
this inhomogeneity effect, the oscillations in the time integrated response
are damped and the signal is dominated by the monotonic decay (except
for short delay times). Correspondingly, all the non-zero frequency
peaks in the Fourier spectrum are strongly broadened. However, the
central peak remains unaffected.

\section{Discussion and conclusions.}
\label{sec:discusion}

An interesting point in the discussion presented above is the presence
of a zero-frequency
component which produces a central peak in the Fourier transform,
composed of two Lorentzian contributions: one with the width
$(\Gamma+\beta)/2$ and another one with $(\Gamma-\beta)/2$.
Interestingly, as follows from Eq.~\eqref{ro13}, there are no
zero-frequency components in the polarization evolution between the
pulses. Therefore, the non-oscillatory part of the FWM response can be
interpreted as a result of partial refocusing of the Larmor precession
by the probe pulse.

The presence of this zero-frequency component, which is insensitive to
the inhomogeneous distribution of the $g$-factor, opens a possibility of
extracting useful information on the spin-related system parameters.
If the trion lifetime is sufficiently long compared to the hole spin
decoherence times ($\beta$ not too small compared to $\Gamma$) then
the two Lorentzian components of the central line, with the widths
$\Gamma\pm\beta$ can be separated by fitting, from which the values of
$\Gamma$ and $\beta$ can be deduced. The value of $\beta$, along with
the detailed balance relation, allows one to extract the rates
$\kappa_{\pm}$, hence the longitudinal decoherence time. On the other
hand, $\Gamma$ involves both the spin-related rates $\kappa_{\pm,0}$
and the optical
lifetime and dephasing rates $\gamma_{0,1}$. The latter, however, can
often be deduced independently for a given system, which can allow one
to find also the value of $\kappa_{0}$ and to calculate the
intrinsic transverse spin decoherence time $T_{2}$.

This should be compared with the information available from a TRKR
experiment \cite{machnikowski10}. There, the hole-related response
originates from the occupations of the hole Zeeman sublevels after the
trion recombination. Therefore, the TRKR experiment allows one to
extract the hole spin-related information even if the exciton lifetime
and coherence time are very short, as was indeed the case in the
experiments \cite{syperek07,kugler09}. However, as the information on
the transverse spin dephasing in that experiment is deduced from the
spin precession signal, it is convoluted with the inhomogeneous effect
which may even completely dominate the decay of coherent precession
signal. Therefore, in a TRKR experiment essentially only the
inhomogeneous transverse decoherence time $T_{2}^{*}$ is available.

Thus, we have shown that a FWM experiment performed on an ensemble of
quantum dots or other trapping centers doped with excess holes can provide
useful information on the properties of hole spin precession and
decoherence. Although this information is only available under
favorable conditions related to the various decoherence rates in the
system (long exciton coherence) it is insensitive to inhomogeneous
effects, in particular to a variation of $g$-factors. Therefore, we
conclude that the FWM experiments may be a useful tool to extract
spin-related information which is complementary to that available from
the TRKR study.

\section*{Acknowledgments.}

This work was supported in part by the
Alexander von Humboldt Foundation within a Research Group Linkage
Grant. P.M. acknowledges support from the TEAM programme of the Foundation for
Polish Science, co-financed from the European Regional Development
Fund.

\providecommand{\WileyBibTextsc}{}
\let\textsc\WileyBibTextsc
\providecommand{\othercit}{}
\providecommand{\jr}[1]{#1}
\providecommand{\etal}{~et~al.}


\section{Introduction}

The boron to carbon nuclear ratio is very important for understanding the
cosmic ray propagation in the Galaxy.
Boron and carbon are respectively a pure secondary and a pure primary cosmic ray species.
Since B is directly produced by C spallation in the interstellar medium, their relative abundance is
sensitive to the traversed matter from the source, and the B/C energy distribution is an
interesting quantity to study cosmic rays diffusion and interactions in the galactic medium.
A precise measurement of the B/C ratio constraints the propagation parameters and, as a
consequence, provides a better understanding of the astrophysical processes at the sources.
The AMS-01 measurement of the B/C ratio is presented here in the kinetic energy range 0.4-19 GeV/n.

\section{The AMS-01 detector}

The AMS-01 spectrometer was composed of a cylindrical permanent magnet, a silicon microstrip tracker,
time of flight (TOF) scintillator planes, a Cerenkov counter and anti-coincidence counters \cite{ref_AMS01}.
The magnet (inner diameter $1.1 m$) provided a central dipole
field with an analyzing power $BL^2=0.14 Tm^2$. Six layers of double-sided silicon microstrip tracker
measured the trajectory of charged particles with an accuracy
of $10 \mu m$ ($30 \mu m$) in the (non)bending coordinates, as well as providing multiple energy loss measurements.
The TOF system had two orthogonal planes at each end of the magnet, and measured the particle
transit time with an accuracy of $\sim$120 psec; they also yielded energy loss measurement up to $|Z|$=2.
A layer of anti-coincidence scintillation counters lined the inner surface of the magnet.
A thin carbon layer were used as a shield to absorb low energy particles.
\begin{figure}[!htbp]
\centering
  \includegraphics[width=0.45\textwidth]{icrc0182_fig01.eps}
\caption{The AMS-01 detector}\label{ams01}
\end{figure}

The detector response was simulated under a \texttt{GEANT3}-based software
\cite{ref_GEANT}. The effects of energy loss, multiple scattering, nuclear interactions and decays
were included in the code, as well as efficiencies, resolutions and reconstruction algorithms.

\section{Event reconstruction and selection}
The main physical characteristics of a cosmic ray particle traversing the detector
are the arrival direction ($\theta$,$\phi$), the particle identity and its
kinetic energy or rigidity. These quantities were reconstructed combining the independent
measurements provided by the various sub-detectors.
The particle rigidity $R=pc/Ze$ was given by the deflection of the particle
trajectory reconstructed from hits in at least four out of the six tracker planes.
The velocity $\beta=v/c$ was determined from the transit time in
the TOF along the track lenght.
The particle charge $|Z|$  was obtained by the analysis of the multiple measurements of energy loss,
as described in detail in the next section.
\\

In order to reject poorly reconstructed particles, a set of quality cuts was applied to the sample:
\begin{enumerate}

\item particles were required to be measured by fully efficient TOF scintillator paddles; this rejects events
passing through a single TOF paddle providing unreliable timing information
\item tracks fitted with high ${\chi}^{2}$ were removed according to a rigidity dependent cut-off value on the ${\chi}^{2}$;
\item tracks with poor \emph{half rigidity agreement}\footnote{For a given track, the half rigidities $R_1$ and $R_2$ are two different measurements of the particle rigidity; $R_1$ ($R_2)$ was obtained with hits from the upper (lower) silicon planes.} were also removed, according to rigidity-dependent criteria;
\item consistency between independent $\beta$ and R measurements was required,
being $\beta$ (particle velocity) and R (rigidity) obtained by two distinct subdetectors;
\item the absence of hits was required close to the track extrapolation on missing layers,
namely on tracker layers with no hit belonging to the reconstructed track.

\end{enumerate}
For the present work we considered data collected during these two periods:
(a) \emph{MIR-docking phase}, 4 days while the shuttle was docked with the MIR space station, and
(b) \emph{post-docking phase}, 2 days pointing at fixed directions ($0^0$, $20^0$ and $45^0$ with respect to the zenith).

Data collected while AMS-01 was passing in the region of the South Atlantic Anomaly were excluded
(latitude: $-47^0 \div -4^0$, longitude: $-85^0 \div -5^0$)
and only downward-going particles (in the AMS-01 reference frame)
were considered, within a restricted acceptance of $32^\circ$.
Approximately 75,000 nuclei from Lithium to Oxygen were selected with the
above mentioned criteria.

\section{Particle identification}

A nuclear species is defined by its charge Z. In this analysis
the charge of each species is determined by the energy losses recorded in the
silicon layers, due to the limited dynamical response of the
TOF scintillators.
The ionization energy generated by a charged particle in a silicon microstrip detector is
collected by a \emph{cluster} of adjacent strips. Tracker clusters were recognized online
and then re-processed with the reconstruction software.
A multi-step procedure of normalization of the cluster amplitudes was performed.
The method accounts for saturation effects, electronics response, particle inclination and velocity dependence
of the energy loss.
The charge identification algorithm, applied to the corrected signals, was based on the \emph{maximum likelihood} method
which determines the best Z value corresponding to the maximum value of the \emph{log-likelihood function}:
\begin{equation}\label{equa1}
L(Z) = \log_{10}\{\prod_{k=1}^{k=6}P_{Z}^{k}(x_{k},\beta)\}
\end{equation}
The k-index runs over the silicon layers, $x_{k}$ are the corrected
cluster amplitudes as observed on k-th layer.
The $P_{Z}^{k}(x_{k},\beta)$ functions have to be viewed as the probabilities of a given charge
Z with velocity $\beta$ to produce a signal $x_{k}$ on the k-th layers. These probability density
functions were estimated from a clean reference sample of flight data. For a fixed $\beta$,
they have gaussian-like shape; their dependence on $\beta$ is shown on Fig.\ref{dedxvsbeta}.
Due to inefficiencies in charge collection,
some energy losses may produce charge responses which do not carry reliable
information on the particle charge. For this reason, not all
the six layer clusters were used to determine Z.
Firstly, a selection of the tracker clusters was done according to (1) reasonable strip occupancy levels
(clusters containing noisy or dead channels were removed) and (2) cluster morphology (double-peaked clusters
and single-strip clusters were removed).
Finally the \emph{best} set of three clusters given by
$\Omega=\{k_1,k_2,k_3\}$ was recognized. $\Omega$ was provided by the maximum likelihood method as a
parameter as well as Z.
Hence, the actual \emph{log-likelihood} function used was:
\begin{equation}\label{equa2}
L(Z,\Omega) = \log_{10}\{P_{Z}^{k_1}\cdot P_{Z}^{k_2}\cdot P_{Z}^{k_3} \}
\end{equation}
where Z ranges from 3 (lithium) to 8 (oxygen) and $\Omega$ runs over all the
$\{k_1,k_2,k_3\}$ 3-fold combinations of the tracker signals.
It should be noted that, although the final Z estimation was provided by three tracker clusters,
all the selected clusters were processed in Eq.\ref{equa2}.

\begin{figure}[!h]
\centering
  \includegraphics[width=0.45\textwidth]{icrc0182_fig02.eps}
\caption{
Signal amplitude (ADC) of the mean energy loss in the silicon tracker vs velocity.
Nuclear families fall into distinct charge bands.
MPVs of the $P_{Z}^{k}(x_{k},\beta)$ functions are superimposed for Z=3 to 8 (dashed lines).
}\label{dedxvsbeta}
\end{figure}

\begin{figure}[!htbp]
\centering
  \includegraphics[width=0.45\textwidth]{icrc0182_fig03.eps}
\caption{
Charge histograms in the B and C region.
The signal amplitudes of Fig.\ref{dedxvsbeta} are here equalized to $\beta\equiv 1$.
Boron and Carbon families fall in distinct charge peaks.
}\label{dedx_bc}
\end{figure}

\section{Charge contamination}

To obtain the differential fluxes impinging the detector, the influence of charge contamination
must be accounted for.
Each nucleus of charge $3\leq Z\leq 8$ produces a charge estimation $\hat{Z}$ which is
related to its \emph{true} impinging charge $Z$ by using a set of coefficients $F_{Z}^{\hat{Z}}$.

The 6$\times$6 matrix $\|F\|$ is diagonally-dominated, and each off-diagonal element
$F_{Z}^{\hat{Z}}$ represents the probability of a nuclear species $Z$ to be
mis-identified as $\hat{Z}$ due to interactions in the detector material and fluctuations
of the energy loss:
\begin{equation}\label{equa3}
F_{Z}^{\hat{Z}}=P(\hat{Z}|Z)
\end{equation}

Two different contributions produce a charge migration $Z\rightarrow \hat{Z}$ :
\\

\begin{enumerate}
  \item after interacting in TOF material, an incoming nucleus Z may fragment and physically turn into $\hat{Z}$.
These events were typically removed by the anticoincidence veto. However, a fraction of them produces
a clean track on the tracker, passing trigger and selection. Since $\hat{Z}$$<$$Z$, the corresponding matrix is triangular.
  \item due to energy loss fluctuations and broadening of the charge response, the overlapping
of the measured $dE/dx$ distributions produces an additional $Z\rightarrow \hat{Z}$ migration probability.
 This effect is tipically symmetric and in most cases $\hat{Z}=Z\pm 1$.
\\

\end{enumerate}

Both the effects were studied with a Monte Carlo simulation where nuclear interactions and
detector responses were reproduced. An upper limit to the contribution (2) was also estimated by the data.
In principle, by the exact knowledge of $F_{Z}^{\hat{Z}}$ coefficients, it would be possible to correct the
measured abundances $\hat{Z}$ to get the corrected ones (e.g. inverting the $\|F\|$ matrix).
Since the overall distortion is estimated to be $2\% \div 6\%$ for the ratio B/C (i.e. smaller than
the statistical uncertainties), this effect was included in the systematic errors.

\section{Flux calculation}

Detector counts are relate to the differential energy flux as a function of the kinetic
energy per nucleon $E_{kn}$ according to the following relation:

\begin{equation}\label{equa4}
\Phi^{Z}(E_{kn})=\frac{N^{Z}(E_{kn})}{A^{Z}(E_{kn})\cdot T^{Z}(E_{kn})\cdot \Delta E_{kn}}
\end{equation}
where $N^{Z}(E_{kn})$ is the number of the particles detected after selection in the energy bin around $E_{kn}$,
$A^{Z}(E_{kn})$ is the detector acceptance in the selected angular range, $T^{Z}(E_{kn})$ is the exposure time and $\Delta E_{kn}$ is
the width of the energy window around $E_{kn}$.

The acceptance $A^{Z}(E_{kn})$ is the convolution of the geometrical factor with the energy dependent
efficiency $\epsilon$, which includes
(1) trigger efficiency, (2) reconstruction efficiency and (3) selection efficiency.
A Monte Carlo estimation of the acceptance was performed according to the formula \cite{ref_Sullivan}:
\begin{equation}\label{equa5}
A^{Z}(E_{kn})=A^{0}\frac{N_{D}^{Z}(E_{kn})}{N_{G}^{Z}(E_{kn})}
\end{equation}
where $N_{G}^{Z}(E_{kn})$ is the number of Z-charged nuclei generated in the energy bin $E_{kn}$
and $N_{D}^{Z}(E_{kn})$ the number of those selected.
Trajectories were generated from an ideal plane placed above the detector,
and $A^0=47.78$ $m^2sr$ is the top plane geometrical acceptance.
To calculate $A^{B}$ and $A^{C}$, $2\cdot 10^{8}$ boron and carbon trajectories were generated
in $2\div 250$ GeV/c momentum range, according to an isotropic distribution and a
power law momentum distribution with spectral index $\gamma=-1$.

The geomagnetic cut-off, varying in the 1-15 GV interval along the AMS-01 orbit, introduces different
distortions of the measured energy spectrum for different nuclear species.
The vertical cut-off was calculated in the
eccentric dipole approximation \cite{ref_Smart}:
\begin{equation}\label{equa6}
R_{vc}=15.0\cos^{4}(\lambda)/R^{2}
\end{equation}
where $\lambda$ is the geomagnetic latitude and R is the distance from the dipole center in terrestrial radii units.
To remove any distortion due to geomagnetic field, only particles with rigidity greater
than $1.2 R_{vc}$ were accepted.
This means that, for a given species Z in the energy window $\Delta E_{kn}$, only some geomagnetic regions
were accessible by AMS-01 during the orbit.
The exposure time $T^{Z}(E_{kn})$  is then defined as the overall time spent by AMS-01 in the regions
above cut-off for a certain species Z in the energy window around $E_{kn}$.
These $T^{Z}$ were also corrected for the trigger deadtime fraction.
\\

The B/C ratio is then derived as the ratio of the differential fluxes

\begin{equation}\label{equa7}
B/C = \Phi^{B}(E_{kn})/\Phi^{C}(E_{kn})
\end{equation}
Two distinct measurements were made:
(a) in the $0.4\div 2.0$ GeV/n range, the $E_{kn}$ was directly obtained from $\beta$ measurement;
(b) between $2.0$ GeV/n and $19.0$ GeV/n the $E_{kn}$ was calculated by the rigidity from tracker.
In the latter case, since the conversion $R\rightarrow E_{kn}$ requires the particle mass,
an assumption on the B and C isotopic compositions is needed\footnote{
Rigidity is the quantity directly measured by the tracker; it was translated into kinetic energy per nucleon according to the relation: $E_{kn}(R)=\sqrt{R^{2}(Z/A)^{2}+M_{p}^{2}}-M_{p}$, being $M_p$= nucleon mass and A= mass number.}.
According to various measurements and predictions (e.g. \cite{ref_SMILI, ref_Voyager, ref_GALPROP}),
the cosmic ray boron is composed of a non-negligible fraction of $^{10}B$. We considered B as a $^{10}B+^{11}B$
isotopic mixture and assumed a $^{10}B$ isotopic fraction as
\begin{equation}\label{equa8}
Y_{B}=^{10}B/(^{10}B+^{11}B)=0.35\pm0.15
\end{equation}
whereas the carbon flux is assumed to be composed of pure $^{12}C$. The uncertainty on $Y_{B}$ is consequently
transposed as additional errors on the B/C ratio ($3\%\div5\%$ relative errors).

\section{Results and discussion}
The results for B/C ratio are presented in Fig.\ref{bc_ratio}.
The ratio B/C is calculated in ten energy intervals between 0.4 and 19.0 GeV/n.
The energy range was decided according to restrictions due to signal saturations
(at lower energies) and requirement of reasonable statistics (at higher energies).
The binning is uniform in logarithm scale; for each data point the reference value
$\langle E_{kn} \rangle$ is defined as the geometric mean of the bin margins.
The results are compared with HEAO-3-C2 data \cite{ref_HEAO} collected
from October 1979 and June 1980.
Both the detectors operated in space during periods of similar solar activity and polarity;
altitude and inclination of the orbit were also comparable.
The two measurements are in substantial agreement.
A theoretical prediction, calculated with the \texttt{GALPROP} numerical software \cite{ref_GALPROP}
(conventional diffusion model modulated with $\Phi=450 MV$ \cite{ref_Gleeson}) is also reported.
The dashed line is referred to a diffusion model (D. Maurin \etal \cite{ref_Maurin}) where
the acceleration and diffusion processes are considered in terms of a rigidity-dependent diffusion coefficient
\begin{equation}\label{equa9}
 K(R)=K_{0}\beta \left( \frac{R}{1GV} \right) ^{\delta}
\end{equation}
and a power-law distribution $R^{-\alpha}$ is assumed at the source;
$\delta=0.6$ and $\alpha=2.2$ are used; the solar modulation is described in the force-field
approximation using $\Phi=500MV$, consistent with the period of the AMS-01 flight.

Both statistical and systematic uncertainties are represented in the Fig.4.
Most of the systematic uncertainties arising from the detector response cancel into the B-C ratio. Differences in the trigger and selection efficiencies of the two species are expected from the Z dependence of the delta-rays production and fragmentation effects in the detector material. Deep investigation on these effects has been undertaken with GEANT4 and FLUKA \cite{ref_FLUKA} based simulations with the aim of obtaining additional estimations of acceptances and their related systematics.

Finally, since a considerable amount of light nuclei were recognized in the same
energy range, one may expect to extend this analysis to other nuclear species.

\begin{figure}[htbp]
\centering
  \includegraphics[width=0.55\textwidth]{icrc0182_fig04.eps}
\caption{B/C ratio from this work (filled circles), from HEAO-3-C2 experiment (open diamonds) and from two diffusion models (solid and dashed lines).}\label{bc_ratio}
\end{figure}

\section{Aknowledgements}
The support of INFN, Italy, ETH-Zurich, the University of
Geneva, the Chinese Academy of Sciences, Academia Sinica
and National Central University, Taiwan, the RWTH Aachen,
Germany, the University of Turku, the University of Technology
of Helsinki, Finland, the US DOE and MIT, CIEMAT,
Spain, LIP, Portugal and IN2P3, France, is gratefully acknowledged.
The success of the first AMS mission is due to many individuals
and organizations outside of the collaboration. The support
of NASA was vital in the inception, development and operation
of the experiment. Support from the Max-Planck Institute for
Extraterrestrial Physics, from the space agencies of Germany
(DLR), Italy (ASI), France (CNES) and China and from CSIST,
Taiwan also played important roles in the success of AMS.
\\


\section{Introduction}

MWC 349A is a peculiar B[e] star which happens to be the brightest continuum
source of its kind at cm wavelengths (e.g.\ Braes et al.\ 1972; Tafoya et al.\ 2004).
It is located in the direction of the Cyg OB2 association, at a distance which has
been estimated between 1.2 kpc (Cohen et al.\ 1985) and 1.7 kpc (Hofmann et al.\
2002). MWC349A may have a B0 III companion (named MWC 349B), located
$2\rlap.{''}4$ to its west  (Merrill et al.\ 1932; Cohen et al.\ 1985). The cm radio
continuum emission from MWC 349A can be successfully described in terms of
an isothermal wind photoevaporating from the surface of a neutral Keplerian disk
surrounding a young massive star (Lugo et al. 2004).

MWC 349A also exhibits hydrogen recombination line emission detectable from
the radio (Altenhoff et al.\ 1981) to the optical (Hartmann et al.\ 1980). The
intensity of the recombination lines at cm wavelengths was found to be consistent
with nearly thermal, LTE conditions (e.g.\ Altenhoff et al.\ 1981; Escalante et al.\
1989; Rodr\'\i guez \& Bastian 1994). However, in a remarkable discovery,
Mart\'\i n-Pintado et al.\ (1989a) found that the mm recombination lines towards
MWC 349A were emitting in a maser mode, and were therefore tens of times
stronger than expected for thermal conditions. This discovery was rapidly confirmed,
and it is now known that all the hydrogen recombination lines at mm, sub-mm, far-IR
and mid-IR wavelengths (from H39$\alpha$ at about 100 GHz to H7$\alpha$ at
about 16,000 GHz) are masing (Gordon et al.\ 2001; Thum et al.\ 1994a; Strelnitski et
al.\ 1996a, 1996b; Thum et al.\ 1998). These maser lines are also known to be time
variable (Mart\'\i n-Pintado et al.\ 1989b; Gordon et al.\ 2001). To this date, MWC
349A remains the only source known to show such high-gain natural hydrogen
maser emission. The line amplification has been investigated theoretically by Thum et
al.\ (1994b) and Ponomarev et al.\ (1994), and is believed to arise in the inner region
of a nearly edge-on disk where the eastern region is redshifted with respect to the
western one (Planesas et al.\ 1992; Weintroub et al.\ 2008).

In an unexpected turn of events, Mart\'\i n-Pintado et al.\ (1993) reported that the
H66$\alpha$ recombination line at 22 GHz was about 5 times stronger than expected
from thermal emission, and argued that this line must, therefore, also be masing. This
was surprising and implied that the maser effect could still be present at wavelengths
as large as 1.3 cm. The observations reported by Mart\'\i n-Pintado et al.\ (1993),
however, were taken with the Very Large Array (VLA) at the time when it only had
limited spectral capabilities. In particular, the observations were made with a total
velocity coverage of only about 160 km s$^{-1}$, barely enough to include the entire
width of the H66$\alpha$ line. This restriction could have led to an improper continuum
subtraction and to an overestimate of the line intensity. In this {\it Letter}, we take advantage
of the greatly enhanced spectral capabilities of the Expanded Very Large Array (EVLA) to
reobserve this important source and obtain accurate H66$\alpha$ line parameters.

\section{Observations}

The data were collected with the Expanded Very Large Array in its most compact (D)
configuration on June 16, 2010 while 22 antennas were available in the array. Five
of these antennas had to be flagged because they consistently delivered low visibility
amplitudes or because they exhibited significant phase jumps. As a consequence,
the final data products were constructed from observations collected on 17 antennas
only. The correlator was set up to deliver 256 spectral channels, each 500 kHz wide,
centered on the rest frequency of the H66$\alpha$ line at 22.364 GHz. This
corresponds to a total velocity coverage in excess of 1,700 km s$^{-1}$ (about 10 times
the velocity coverage of Mart\'\i n-Pintado et al.\ 1993), and a spectral resolution of
about 6.7 km s$^{-1}$. The data were calibrated with the Common Astronomy Software
Applications (CASA) package following standard procedures. The bandpass
was measured using observations of J1743--0350 and J0319$+$4130 (3C84) obtained
at the beginning and the end of the observations. The flux scale was set using an
observation of J0137$+$3309 (3C48). Finally, the amplitude and phase gains were
monitored using observations of J2015$+$3710 (with a bootstrapped flux density of 2.3 Jy)
obtained roughly every 10 minutes throughout the 3 hours of observations.

Both the continuum and the H66$\alpha$ line were clearly detected in the observations,
so two data sets were created from the calibrated visibilities. A continuum data set
was generated by averaging 137 line-free spectral channels (corresponding to
an equivalent bandwidth of about 70 MHz). Simultaneously, a  continuum-free
spectral line data set was created by subtracting a constant continuum (measured over
the same 137 line-free channels) from the initial data. Both data sets were then
imaged using robust weighting and 0.5 arcsecond pixels. The resulting angular
resolution was \msec{3}{35} $\times$ \msec{2}{07}; P.A.\ --60.7$^\circ$. To improve
the overall calibration, the continuum data set was self-calibrated --in phase
only-- and the resulting gain corrections were applied both to the continuum
and the continuum-free spectral line data sets which were then re-imaged.

\section{Results}

The continuum image shows the presence of a strong point-like
source at the expected position of MWC 349A ($\alpha_{J2000.0}$ =
\dechms{20}{32}{45}{52}; $\delta_{J2000.0}$ = \decdms{40}{39}{36}{62}). A 2D
Gaussian fit to this source yields an integrated flux of 411 mJy with a formal error of
a fraction of a mJy. This error, however, does not account for the systematic uncertainty
on the absolute flux calibration. We estimate this systematic uncertainty to be about
10\%, and conclude that the 22.3 GHz flux of MWC 349A is $F_\nu$ = 411 $\pm$ 41
mJy. This value is in good agreement with the flux obtained in previous 1.3 cm observations
of MWC 349A. For instance, Tafoya et al.\ (2004) reported a 22.4 GHz flux density of
446.2 $\pm$ 44.8 mJy. We note, also, that changes in the morphology of MWC 349A at
cm wavelengths are known to have occurred (Rodr\'{\i}guez et al.\ 2007), so moderate
variations in its flux density might also be present.

In the continuum-free spectral data set, emission is clearly detected in more than 40
contiguous spectral channels around the frequency of the H66$\alpha$ line. This
emission remains unresolved at the angular resolution of our observations, so a
single-pixel spectrum at the peak position of the source (Figure 1) provides a complete
description of the line emission. This spectrum shows a very prominent double-peaked
feature surrounded by a 1,200 km s$^{-1}$ emission-free spectral baseline.
This baseline is remarkably linear but not flat. Instead, it shows a slope which most certainly
results from the combined effect of the true positive spectral index of MWC 349A (e.g.\
Tafoya et al.\ 2004) and the uncorrected negative spectral indices of the flux and
bandpass calibrators.\footnote{In the future, CASA will be able to take into account the
spectral dependences of the calibrators, and it will become possible to measure directly
the spectral index of the target from observations such as those presented here.}
A spectrum with a linear baseline removed  is shown in Figure 1b. The VLA observations
reported by Mart\'\i n-Pintado et al.\  (1993) only covered the spectral bandwidth shown
by the thick horizontal bracket on Figure 1b. It is clear that this spectral window barely contains
the most prominent spectral line detected here, and does not include at all the weaker
feature seen at more negative velocities. Since observations with such a limited velocity
coverage would contain no channels free of line emission, it would clearly be nearly
impossible to properly subtract a continuum from them.

\begin{figure}[!ht]
\begin{center}
\includegraphics[width=0.65\linewidth,angle=0]{f1.eps}
\end{center}
\caption{(a) Raw spectrum of MWC 349A centered on the rest frequency of the H66$\alpha$
line. The best fit to the slope in the spectral baseline with a straight line is shown as a solid
black line. (b) Baseline-subtracted spectrum.
The velocity coverage of the VLA data presented by  Mart\'\i n-Pintado et al.\ (1993) is shown
as a horizontal bracket below the spectral lines. The blue solid line shows a fit with two unconstrained
Gaussian components. The red solid line shows a similar fit but with the velocity separation
between the two components constrained to be 122.2 km s$^{-1}$, the expected velocity
separation between the H66$\alpha$ and the He66$\alpha$ lines. }
\end{figure}

A simultaneous fit to the spectrum with two Gaussian components (blue solid line in Figure 1b) provides
the parameters listed in Table 1. The separation in velocity between the two spectral features
is 130.4 $\pm$ 4.4 km s$^{-1}$, very similar to the separation (122.2 km s$^{-1}$) expected
between the hydrogen and helium 66$\alpha$ recombination lines. Indeed, a fit where the
separation between the two spectral components is constrained to be 122.2 km s$^{-1}$
(red solid line in Figure 1b) is very nearly as good as a fit where the separation can vary freely
(the parameters corresponding to that constrained fit are also given in Table 1). We conclude
that the weaker spectral feature seen at more negative velocities in our spectrum is the
He66$\alpha$ line. We note that Thum et al.\ (1992) also identified a blue-shifted feature in
their H41$\alpha$ spectrum of MWC 349A, and favored an interpretation where that additional
component is the He41$\alpha$ line. The present results clearly lend strong support to that
interpretation. Interestingly, however, the He41$\alpha$ line detected by Thum et al.\
(1992) appears to be at a somewhat more negative radial velocity than its hydrogen counterpart
($V_{lsr}$ = --31 km s$^{-1}$ vs.\ $V_{lsr}$ = $+$1 km s$^{-1}$). If our unconstrained fit is taken at
face value, the He66$\alpha$ line is also at a more negative radial velocity than the H66$\alpha$
line ($V_{lsr}$ = --15.2 $\pm$ 4.3 km s$^{-1}$ vs.\ $V_{lsr}$ = --7.0 $\pm$ 0.7 km s$^{-1}$), but
the difference is significantly less than for the 41$\alpha$ lines. We will come back to that point in
Section 4.3. We note, finally, that the H66$\alpha$ line also shows a hint of excess emission on
the redshifted side that may be related to the winds predicted in the models of Avalos \& Lizano
(in preparation).

\begin{table}[!t]
\caption{Spectral line parameters}
\begin{center}
\begin{tabular}{lrrr}
\hline
\hline
 & \multicolumn{1}{c}{$S_\nu$ (mJy)} & \multicolumn{1}{c}{$V_{lsr}$ (km s$^{-1}$)} & \multicolumn{1}{c}{$\Delta V$ (km s$^{-1}$)\tablenotemark{a}}\\%
\hline
\multicolumn{1}{l}{Unconstrained fit:}\\%
Feature 1 & 24.9 $\pm$ 0.3 & --7.0 $\pm$ 0.7 & 90.0 $\pm$ 1.7 \\%
Feature 2 & 4.5 $\pm$ 0.5 & --137.4 $\pm$ 4.3 & 61.3 $\pm$ 9.9 \\%
\multicolumn{1}{l}{Constrained fit:}\\%
H66$\alpha$ & 25.0 $\pm$ 0.3 & --7.1 $\pm$ 0.7 & 88.7 $\pm$ 1.6 \\%
He66$\alpha$ & 4.4 $\pm$ 0.4 & \multicolumn{1}{c}{$''$} & 62.8 $\pm$ 8.2 \\%
\hline
\hline
\tablenotetext{a}{$\Delta V$ is the full width at half maximum.}
\end{tabular}
\end{center}
\end{table}

\section{Discussion}

\subsection{Excitation of the H66$\alpha$ line}

The H66$\alpha$ line to continuum ratio reported by Mart\'\i n-Pintado et al.\ (1993)
is typically about 20\%. In contrast, the line peak temperature found here (25 mJy) is
only about 6\% of the underlying continuum. A somewhat more robust line to
continuum ratio of about 5\% is obtained by averaging over the central 90 km
s$^{-1}$ of the line. In LTE conditions and adopting a full width at half power of
90 km s$^{-1}$ and an electron temperature of 8,000 K, the H66$\alpha$ line is
predicted to represent about 5\% of the underlying continuum (Altenhoff et al.\
1981; Rodr\'\i guez 1982; Mart\'\i n-Pintado 2002) in excellent agreement with the
observed values. We conclude that the H66$\alpha$ line in MWC 349A is largely
thermally excited and is not strongly masing as had been suggested by
Mart\'\i n-Pintado et al.\  (1993).

One could argue, however, that our results and those of Mart\'\i n-Pintado et al.\ (1993)
are not necessarily incompatible because the two sets of observations might trace
different regions of MWC 349A. In particular, the present observations were obtained
in the most compact (D) configuration of the EVLA and are, therefore, sensitive to even
fairly extended emission. If it existed, such extended emission would have been filtered out in the
observations reported by Mart\'\i n-Pintado et al.\  (1993), because they were collected in
the most extended (A) configuration of the VLA. Could our observations be dominated
by extended, thermally excited gas, and those of Mart\'\i n-Pintado et al.\ (1993) by denser,
more compact masing regions?

This appears not to be the case for the following reasons. At a wavelength of 1.3 cm,
the (E)VLA in the A configuration is sensitive to any structure more compact than about
2.5 arcseconds\footnote{http://evlaguides.nrao.edu/index.php?title=Category:Status}.
Our observations, on the other hand, show that the line emission
detected here originates in a region more compact than about 1.5 arcseconds. Thus,
there are no spatial scales contributing to the emission reported here and to which the
observations by Mart\'\i n-Pintado et al.\  (1993) were not sensitive. {\em The two sets of
observations trace essentially the same gas.} Moreover, at 1.3 cm, the size of MWC
349A is about \msec{0}{2} $\times$  \msec{0}{2} (e.g.\ Tafoya et al.\  2004;  Rodr\'{\i}guez
et al.\ 2007; Mart\'\i n-Pintado et al.\ 1993). Since the angular resolution of the (E)VLA
at 1.3 cm in the A configuration is about \msec{0}{1}, MWC 349A only contains 4 independent
resolution elements at that wavelength. Thus, the four spectra reported by Mart\'\i n-Pintado
et al.\  (1993) in their Figure 2 trace most of the extent of the source, and the continuum emission
integrated over these four positions will be about 300 mJy (the total continuum emission
estimated by Mart\'\i n-Pintado et al.\ 1993). According to the spectra shown in Figure 2
of Mart\'\i n-Pintado et al.\ (1993), the mean intensity of the H66$\alpha$ line over the
velocity range from $-$60 to $+$50 km s$^{-1}$ is about 20\% of the underlying
continuum. The line intensity integrated over these 4 positions must therefore be
about 60 mJy. As a consequence, the integrated line flux measured over any larger area should
be at least 60 mJy in that velocity range.\footnote{Recall that flux densities measured
in Jy do not get diluted: if a 1 Jy source is present at the center of a given field, the total
flux density taken over any area containing the source will be 1 Jy, independently of
how large the chosen area might be.} Yet, in our spectrum (corresponding to a
\msec{3}{35} $\times$ \msec{2}{07} beam), the flux density is systematically smaller
than 25 mJy; averaged over the $-$60 to $+$50 km s$^{-1}$ velocity range, it is 18
mJy. We conclude that Mart\'\i n-Pintado et al.\  (1993) have overestimated the line
intensity by a factor 3 to 4.

We can only speculate about the exact reasons which led Mart\'\i n-Pintado et al.\
(1993) to over-estimate the H66$\alpha$ line intensity in MWC 349A. In principle,
the lack of line-free channels in their observations should have led them to over-estimate
the continuum flux (since the channels on the edge of the bandpass used to calculate
the continuum still contained line emission) and to under-estimate the line to continuum
ratio (rather than the other way around). Instead, a more subtle effect, also related
to the narrow bandwidth of their observations might have been the culprit. The velocity
channels on the edge of the VLA bandpass are well known to be very noisy. A
continuum measurement based on these channels may have seriously underestimated
the continuum flux density, and  adversely affected the line to continuum ratio
measurement.

\subsection{Helium abundance and electron temperature}

Our results show that the excitation of the H66$\alpha$ line
very nearly corresponds to LTE conditions. In this situation, the line behaves
as if it were approximately optically thin (see Rodr\'{\i}guez 1982). Given its
similar excitation and lower abundance, the He66$\alpha$ line will also
behave as if it were approximately optically thin, so the integrated line intensity
ratio provides a direct measure of the ionized helium over ionized hydrogen
abundance ratio $y^+$. From the unconstrained fit  (Section 3), we obtain
$y^+ = 0.12 \pm 0.02$, in good agreement with the value (0.10 $\pm$
0.01) found by Thum et al.\ (1992) for the 41$\alpha$ lines.

As a consequence of its relative abundance, the opacity of the helium lines is
expected to be about 10 times less than that of the hydrogen lines, so the helium
lines are not anticipated to be masing (even at mm wavelengths). This presumably
explains why the helium lines associated with strongly masing hydrogen lines at
mm, sub-mm, far-IR and mid-IR wavelengths have not been detected in existing
observations. Since the hydrogen lines are tens of times stronger than thermal,
the helium to hydrogen line intensity ratio will typically be less than 1\% at mm
wavelengths. The noise level of the mm observations was almost certainly adapted
to detect the hydrogen lines, and was insufficient to reach the much weaker helium
lines.

The high quality of the data presented here allows us to estimate the electron
temperature of the ionized gas in MWC 349A. Following Rodr\'\i guez et al.\
(2009), and assuming LTE conditions, the electron temperature, $T_e^*$, of
an ionized, partially optically-thick outflow is given by:

\begin{equation}\Biggl[{{T_e^*} \over {K}}\Biggr] = \Biggl[{4700 \over (1 + y^+)}
{{S_C} \over {S_L}} \biggl({{\nu_L}
\over {GHz}} \biggr)^{1.1}
 \biggl({{\Delta v} \over
{km~s^{-1}}}\biggr)^{-1} \Biggr]^{0.87}, \end{equation}

\noindent where $\nu_L$ is the line frequency, $S_C$ is the continuum flux
density, $S_L$ is the hydrogen peak line flux density, and $\Delta v$ is the hydrogen
FWHM line width. Using the helium abundance estimated above and the line
and continuum parameters measured in Section 3, we obtain $T_e^* = 6,300 \pm 600$ K.
This value is consistent with those determined from other radio observations
(White \& Becker 1985; Cohen et al.\ 1985; Escalante et al.\ 1989).

\subsection{Relative kinematics of the hydrogen and helium lines}

The kinematics of the ionized gas associated with MWC 349A presents several
problems that require additional research. For instance, the direction of rotation
determined from the cm thermal lines (eastern side blueshifted and western side
redshifted; Rodr\'\i guez \& Bastian 1994) is opposite to that found in the mm maser
lines (Planesas et al.\ 1992; Weintroub et al.\ 2008). Also, while the thermal radio
recombination lines at cm wavelengths peak at negative velocities (typically
$V_{lsr}$ $\sim$ --20 km s$^{-1}$; e.g.\ Rodr\'\i guez \& Bastian 1994), the centroid
of maser mm radio recombination lines is at positive velocities ($V_{lsr}$ $\sim$
$+$10 km s$^{-1}$; e.g.\ Weintroub et al.\ 2008).

In our unconstrained fit, the LSR velocity of the He66$\alpha$ line appears to
be --15.2 $\pm$ 4.3 km s$^{-1}$ whereas the LSR velocity of the He41$\alpha$ line
reported by Thum et al.\ (1992) is --31 $\pm$ 8 km s$^{-1}$. Thus, on average,
the Helium lines appear to peak at $V_{lsr}$ $\sim$ --20 km s$^{-1}$, similar to
the velocity of the thermal cm lines of hydrogen. This strongly suggests that the helium
lines originate in the same gas as the thermal cm hydrogen lines.

The H66$\alpha$ line detected here peaks at a velocity
(--7.0 $\pm$ 0.7 km s$^{-1}$ -- Section 2), intermediate between the $+$10 km s$^{-1}$
masing mm lines, and the --20 km s$^{-1}$ thermal cm lines.
A plausible (albeit somewhat {\it ad hoc}) explanation is that the H66$\alpha$
line is largely dominated by a thermal component at --20 km s$^{-1}$ but also contains a
moderate masing component around $+$10 km s$^{-1}$. The reversed situation might
occur with the H41$\alpha$ line at 92 GHz reported by Thum et al.\ (1992) which peaks
at $+$1.0 $\pm$ 0.9 km s$^{-1}$. In that second case, the emission might be dominated
by the masing line at $+$10 km s$^{-1}$ but with a moderate thermal contribution. It would
be interesting to analyze further the relative kinematics of the hydrogen and helium lines,
for instance by obtaining high resolution EVLA observations of the 53$\alpha$ and
66$\alpha$ lines at 0.7 and 1.3 cm, respectively.

\section{Conclusions and perspectives}

We have presented new EVLA observations of the 22.3 GHz continuum and the
H66$\alpha$ line emissions toward the well-studied Galactic emission-line star MWC 349A.
The continuum flux density (411 $\pm$ 41 mJy) is in good agreement with previous
determinations. The intensity of the H66$\alpha$ line (25 mJy) and the line-to-continuum
flux density ratio (5\%), on the other hand, are those expected for local thermodynamic
equilibrium. This shows that the H66$\alpha$ line excitation is largely thermal, and that
the line is not strongly masing as had been proposed by Mart\'\i n-Pintado et al.\  (1993).
The intermediate radial velocity of the H66$\alpha$ line, however, might suggest the
existence of a moderate maser component. The He66$\alpha$ recombination line is also
detected in our observations, and shows that the ionized helium to ionized hydrogen
abundance ratio in MWC 349A is $y^+$ = 0.12 $\pm$ 0.02. Because of their shared
kinematics, it is likely that the ionized helium gas and the thermally excited ionized
hydrogen gas are well mixed. The electron temperature of that gas is 6,300 $\pm$ 600 K.

High-resolution EVLA observations of this and other recombination lines (particularly
the H53$\alpha$ at 0.7 cm) might help resolve some of the pending issues concerning
the kinematics of the gas in MWC 349A. On a more general note, we would like to stress
that the present observations demonstrate the enormous potential of the EVLA to observe
wide spectral lines. When previous VLA observations had to rely on {\it ad hoc} and
often very uncertain procedures to remove a spectral baseline, the present observations
provide a very clean identification of the line and extremely reliable determinations of
the line parameters.

\acknowledgments We thank Clemens Thum and an anonymous referee for
thoughtful comments on our manuscript, and Vivek Dhawan, Lorant Sjouwerman,
Rick Perley, George Moellenbrock, and Michael Rupen for their help and advices
at various stage of the data calibration. We acknowledge the support of DGAPA,
UNAM, and of CONACyT (M\'exico). LL is indebted to the Guggenheim Memorial
Foundation for financial support.


\section{Introduction}

Let $M$ be a closed (that is, compact and without boundary) Riemannian manifold
with fundamental group $\pi_1(M)$. A residually finite group $G$, a
surjective homomorphism $\pi_1(M)\twoheadrightarrow G$ and a nested sequence of
finite index normal subgroups of $G$ with trivial intersection gives rise
to a tower of  finite sheeted regular covers of $M$; conversely, every tower of
finite sheeted regular
covers arises in this manner. In summary, writing $G_0=G$ and $M_0=M$,
we have:
\begin{equation*}
\label{setup}\tag{$*$}
\begin{minipage}[l]{0.6\linewidth}
  $G_0 \trianglerighteqslant G_1 \trianglerighteqslant G_2 \trianglerighteqslant \cdots$,\;\;\hbox{ with }
  $\displaystyle{\bigcap_{n=0}^{\infty}} G_n= \{\, 1 \,\}$,\;\;\\[2pt]
  $\pi_1(M_n)\twoheadrightarrow G_n, \hbox{ and finite groups } \Gamma_n := G/G_n$,
\end{minipage}
\begin{minipage}[r]{0.4\linewidth}
  \xymatrix{%
     \ar@{..>}[d] &&  \ar@{..>}[d] \\
   \Gamma_2 \ar[d] && M_2 \ar[d] \\
   \Gamma_1 \ar[d] && M_1 \ar[d] \\
   \Gamma_0 && M_0.  }
\end{minipage}
\end{equation*}

In the context of spectral geometry of towers of covers one studies the asymptotic behavior
of the first non-zero eigenvalues $\lambda_1(M_n)$ of the Laplacian,
that is, of the Laplace-Beltrami operator of the individual
Riemannian manifolds $M_n$.  In particular, the following questions are classical:
\begin{ilist}
  \item Does there exist a tower with $\lambda_1(M_n)\geqslant c>0$ uniformly over $n$?
  \item Does there exist a tower with $\lambda_1(M_n)\to 0$ as $n\to\infty$?
\end{ilist}

In this note we are concerned with (b). The earliest positive result on this
question is due to Randol, who studied the case of cyclic covers using
the trace formula \cite{MR0400316}. Subsequent results of Brooks
\cite{MR840402,MR799796} and Burger \cite{MR832070} were obtained by
relating the
eigenvalues $\lambda_1(M_n)$ to combinatorial properties of the Cayley graphs of
finite groups of deck transformations $\Gamma_n$.  Similar results are due to Sunada~\cite{MR782558}.

In all cases, the method to build a tower of covers satisfying (b)
rests on choosing an {\it amenable\/} group $G$ for the construction~(\ref{setup}). Our main result is that it is possible to obtain such a
tower when $G$ is the free group on two generators.
In the statement, $H^{(2)}$ denotes the subgroup of the discrete group
$H$ generated by the squares of its elements.

\begin{thmm}
  Let $M$ be a closed Riemannian manifold, whose fundamental group
  admits a virtual\footnote{A virtual homomorphism is a homomorphism of a finite-index subgroup.} surjective homomorphism onto
  the free group
   of rank 2.  Taking the nested sequence of
  subgroups in \textup{(}\ref{setup}\textup{)}
  to be the sequence of iterated squares in the free group
\begin{equation*}
\label{squares}\tag{$**$}
 G_0=\mathbb{F}_2,\quad G_1 = \mathbb{F}_2^{(2)}, \quad G_2 =(\mathbb{F}_2^{(2)})_{\strut}^{(2)},
       \quad G_3 =((\mathbb{F}_2^{(2)})_{\strut}^{(2)})_{\strut}^{(2)},\quad \cdots
\end{equation*}
  we obtain a tower of covers of
  $M$ for which $\lambda_1(M_n)\to 0$ as $n\to\infty$.   This tower is not
  obtainable up to uniform quasi-isometry \textup{(}or even uniform
  coarse equivalence\textup{)} by the construction
  \textup{(}\ref{setup}\textup{)} with an amenable $G$.
\end{thmm}

\noindent
Observe here that each $G_n$
is normal, even characteristic, in $\mathbb{F}_2$.

The hypothesis of the theorem means that the fundamental group is
{\it large\/} (the terminology is due to Gromov~\cite{gromov}). It
applies to many hyperbolic manifolds~\cite{l}, in particular, to a
closed orientable surface of genus at least two -- the fundamental
group of such a manifold surjects onto $\mathbb{F}_2$.

We conclude the introduction by remarking that in more modern
terminology the classical problems above concerning the construction
(\ref{setup}) can be rephrased in terms of {\it Property\/ $\tau$}:
(a) asks for $G$ to have Property $\tau$ with respect to the family of
subgroups $(G_n)_{n\geqslant 0}$, whereas (b) asks, after perhaps passing to
a subsequence, for $G$ to {\it not\/} have Property $\tau$ with
respect to the $(G_n)_{n\geqslant 0}$.  This is explained in the work of
Burger and Brooks cited above.  Thus, the first assertion in the
theorem is essentially equivalent to the assertion that $\mathbb{F}_2$ does
not have Property $\tau$ with respect to the subgroups appearing in
(\ref{squares}).  For the definition and relevant facts about Property
$\tau$ see \cite{lubotzky}.

\section{Eigenvalues}

A {\it graph\/} is a collection of {\it vertices\/} and {\it edges\/}.
With a small number of exceptions, we permit neither multiple edges
nor loops, so that an edge is uniquely determined by its incident
vertices.  Our graphs are unoriented.
The {\it Cheeger constant\/} of a finite graph $\Gamma$ is
\begin{equation}
\label{cheeger}
  h(\Gamma) = \inf \frac{\# E(A,B)}{\min\{\, \# A, \# B \,\}},
\end{equation}
where the infimum is taken over all decompositions of the vertex set
of $\Gamma$ as a disjoint union $A\sqcup B$ and where, for such a
decomposition, $E(A,B)$ denotes the set of edges with one incident
vertex in $A$ and the other in $B$.

We shall make use of the following result of Brooks which, in the
notation of (\ref{setup}), relates the eigenvalues of the $M_n$ to the
Cheeger constants of the Cayley graphs of the $\Gamma_n$
computed  with respect to the canonical images of generators of $G$ and denoted, by an abuse of notation, again by $\Gamma_n$.
  We shall require only the forward implication, which is the content of
\cite[Lemma~1]{MR840402}.

\begin{thmm}[Brooks]
  In the notation of
  \textup{(}\ref{setup}\textup{)} we have
  $h(\Gamma_n)\to 0$  precisely when $\lambda_1(M_n)\to 0$.  \qed
\end{thmm}

Thus, the first statement in theorem of the introduction is reduced to
the following:

\begin{prop}\label{prop:sq}
  Let $G=\mathbb{F}_2$ be the free group of rank 2.  Consider the tower of
  iterated squares \textup{(}\ref{squares}\textup{)} and the
  corresponding quotients:
\begin{equation*}
    \Gamma_0 = \{\, 1 \,\} \leftarrow \Gamma_1 = \mathbb{F}_2/\mathbb{F}_2^{(2)}
      \leftarrow \Gamma_2 = \mathbb{F}_2/(\mathbb{F}_2^{(2)})_{\strut}^{(2)} \leftarrow \cdots
  \end{equation*}
  Abusing notation, view each $\Gamma_n$ as a Cayley graph
  with respect to the images of the standard free generators of $\mathbb{F}_2$. Then we
  have $h(\Gamma_n)\to 0$ as $n\to\infty$.
\end{prop}

In preparation for the proof we recall the construction of the {\it
$\mathbb{Z}/2$-homology cover\/} of a finite graph $\Sigma$.  Fix a maximal
tree $T$ in $\Sigma$ and let $e_1,\dots,e_r$ be the edges of
$\Sigma$ {\it not\/} in $T$.  The vertex and edge sets of the
$\mathbb{Z}/2$-homology cover $\widetilde\Sigma$ are
\begin{equation*}
  \widetilde V = V \times \oplus_1^r \mathbb{Z}/2, \quad
  \widetilde E = E \times \oplus_1^r \mathbb{Z}/2,
\end{equation*}
where $E$ and $V$ denote the vertex and edge sets of $\Sigma$. Let
$e\in E$ and let $v$, $w\in V$ be the vertices incident with $e$.
Consider the edge $(e,\alpha)\in\widetilde E$. Incidence is defined
in two cases:
\begin{equation*}
  \text{$(e,\alpha)$ contains} \begin{cases}
          \text{$(v,\alpha)$ and $(w,\alpha)$},
                &\text{when $e$ belongs to the maximal tree $T$} \\
          \text{$(v,\alpha)$ and $(w,\alpha+\overline{e}_j)$},
                &\text{when $e=e_j$, for some $1\leqslant j \leqslant r$.}
                          \end{cases}
\end{equation*}
Here $\overline{e}_j=(\dots,1,\dots)$ is the standard basis vector
with a single $1$ in the $j$-the position and $0$'s elsewhere.
Strictly speaking, when defining incidence it is necessary to {\it
  direct\/} the edges $e_j$. It is quickly verified however that,
while the edges are parameterized in a different manner, the
underlying {\it undirected\/} graph is independent of the choice. We
shall not dwell on this aspect.

\begin{rem}
  The construction given here of the $\mathbb{Z}/2$-homology cover is a
  special case of the classical construction of a finite sheeted regular cover of
  $\Sigma$ corresponding to
  a given normal subgroup of finite index in  $\pi_1(\Sigma)$, see, for example,~\cite[Ch. 2]{stillwell}.
  Indeed, with $e_1,\dots,e_r$ as above, and after directing each
  $e_j$, we identify
\begin{equation*}
 \pi_1(\Sigma) \cong \mathbb{F}_r = \langle\; e_1,\dots,e_r \;\rangle.
\end{equation*}
Then the cover corresponding to the kernel of the epimorphism
\begin{equation*}
      \pi_1(\Sigma)\cong \mathbb{F}_r \twoheadrightarrow {\mathbb{F}_r}/{\mathbb{F}_r^{(2)}}\cong
\oplus_1^r \mathbb{Z}/2\quad \hbox{ defined by } \quad
e_j\mapsto\overline{e}_j,
 \end{equation*}
is the $\mathbb{Z}/2$-homology cover.
\end{rem}

\begin{lem}\label{lem:ch}
  Let $\Sigma$ be a finite graph, with vertex set $V$; let
  $\widetilde\Sigma$ be its
  $\mathbb{Z}/2$-homology cover.  We have
  \begin{equation*}
    h(\widetilde\Sigma)\leqslant \frac{2}{\# V}.
 \end{equation*}
\end{lem}

\begin{proof}
  We employ the notation introduced above for $\widetilde\Sigma$.
  We shall exhibit a
  decomposition of the vertex set $\widetilde V = A \sqcup B$ for which the
  quotient in  (\ref{cheeger}) is bounded by $2/\# V$.  Let
  \begin{equation*}
    A = \{\, (v,\alpha)\in \widetilde V: \alpha=(*,\dots,*,0) \;\}, \quad
    B = \{\, (w,\beta)\in \widetilde V : \beta=(*,\dots,*,1) \;\},
  \end{equation*}
each of which contains exactly  $2^{r-1}\, \# V$ vertices. The edges
in $\widetilde E$ with one vertex in $A$ and the other in $B$ are
exactly those of the form $(e_r,\gamma)$, for arbitrary
$\gamma\in\oplus_1^r \mathbb{Z}/2$; thus $E(A,B)$ contains exactly $2^r$
edges.
\end{proof}

\begin{proof}[Proof of Proposition~\ref{prop:sq}]
  The Cayley graph $\Gamma_n$ is the $n$-th iterated $\mathbb{Z}/2$-homology
  cover of the ``figure $8$''.  Since the number of vertices in
  $\Gamma_n$ tends to infinity, the result follows from the previous
  lemma.
\end{proof}

\begin{rem}
  A more detailed analysis gives information on the {\it rate\/} of
  the convergence $h(\Gamma_n)\to 0$.  Indeed, let $V_n$ be the set
  of vertices and $E_n$ the set of edges of (the Cayley graph of)
  $\Gamma_n$.  We have
  \begin{equation*}
    \frac{\# V_{n+1}}{\# V_{n}} =    \frac{\# E_{n+1}}{\# E_{n}}
         = 2^{\rk \pi_1(\Gamma_{n})}.
  \end{equation*}
  Now, the rank
  of the fundamental group $\pi_1(\Gamma_{n})$ is
  the number of edges
  {\it not\/} belonging to a fixed maximal tree in $\Gamma_{n}$.
  Since $\# E_n = 2 \cdot\# V_n$, the rank of $\pi_1(\Gamma_{n})$ is
  $\# V_n+1$.
  Thus, we get the recursive formula
\begin{equation*}
  \# V_{n+1} = \# V_n \cdot 2^{\# V_n +1}.
\end{equation*}
In particular, $\# V_n$ grows faster than an iterated exponential and,
according to the previous lemma, the Cheeger constant
$h(\Gamma_{n+1})$ decays as the reciprocal of $\# V_n$.
\end{rem}

\section{Non uniform coarse equivalence}

We shall now show that the tower constructed in the previous section
cannot be duplicated beginning with an amenable group in
(\ref{setup}), thus completing the proof of the theorem in the
introduction.

 Two families $(X_n)_{n\geqslant 0}$ and $(Y_n)_{n\geqslant 0}$ of metric spaces are {\it
  uniformly   quasi-isometric\/} if there exist functions $f_n:X_n\to Y_n$ and
constants $C\geqslant 1$ and $D\geqslant 0$ such that for all
$x,y\in X_n$ and $z\in Y_n$, we have
\begin{itemize}
  \item $C^{-1} d(x,y) - D \leqslant d(f_n(x),f_n(y)) \leqslant C d(x,y) + D$,
    \item $d(z, f_n(X_n))\leqslant D$.
\end{itemize}
The families $(X_n)_{n\geqslant 0}$ and $(Y_n)_{n\geqslant 0}$ are {\it uniformly coarsely equivalent\/}
if there exist functions $f_n:X_n\to Y_n$ with the following two
properties:
\begin{itemize}
  \item $\forall A\; \exists B$ such that
           $\forall n\;\forall x,y\in X_n$ we have
           $d(x,y)\leqslant A\ \Rightarrow\ d(f_n(x),f_n(y))\leqslant B$,
  \item $\forall A\; \exists B$ such that
           $\forall n\;\forall x,y\in X_n$ we have
           $d(x,y)\geqslant B\ \Rightarrow\ d(f_n(x),f_n(y))\geqslant A$.
\end{itemize}
If two families are uniformly quasi-isometric then they are uniformly
coarsely equivalent.  Observe that these notions apply to individual
spaces, which we regard as trivial families containing a single space.
We say, for example, two spaces are coarsely equivalent.

\begin{prop}[The Uniform \v{S}varc-Milnor Lem\-ma]\label{prop:uqi}
  Continue with the notation of
  \textup{(}\ref{setup}\textup{)}.  Equip each $\Gamma_n$ with the
  word metric associated to a fixed finite
  generating set for $G$; equip each $M_n$ with the path metric
  associated to its Riemannian structure.  The families $(\Gamma_n)_{n\geqslant 0}$ and
  $(M_n)_{n\geqslant 0}$ are uniformly quasi-isometric.
\end{prop}

\begin{rem}
  In the situation of (\ref{setup}) the group $G$ is indeed
  finitely generated.  Further, the statement in the proposition is
  independent of the choice of generators for $G$.
\end{rem}

\begin{proof}[Proof of Proposition~\ref{prop:uqi}]
  The result follows from the \v{S}varc-Milnor Lem\-ma \cite[Prop.~I.8.19]{bh},
  observing that the inherent quasi-isometry constants (see the proof of the Lemma) depend
  only on the diameter of a fundamental domain for the action. In
  detail, $\Gamma_n$ is the group of deck transformations of the cover
  $M_n$ of $M$, whereas $G$ is the group of deck transformations
  of the cover corresponding to the kernel of the surjective
  homomorphism $\pi_1(M) \twoheadrightarrow G$.  Further, the image in
  $M_n$ of a bounded fundamental
  domain for the action of $G$ is a
  fundamental domain for the action of $\Gamma_n$, of no greater diameter.
\end{proof}

Thus, the second statement in the theorem of the introduction is
reduced to the following:

\begin{prop}\label{prop:nqe}
  Consider the tower of iterated squares
  \textup{(}\ref{squares}\textup{)} of the free
  group $\mathbb{F}_2$ and the corresponding quotients
  \begin{equation*}
    \Gamma_0 = \{\, 1 \,\} \leftarrow \Gamma_1 = \mathbb{F}_2/\mathbb{F}_2^{(2)}
      \leftarrow \Gamma_2 = \mathbb{F}_2/(\mathbb{F}_2^{(2)})_{\strut}^{(2)} \leftarrow \cdots
  \end{equation*}
Then the family $(\Gamma_n)_{n\geqslant 0}$ is not uniformly coarsely equivalent to any family of
quotients of an amenable group.
\end{prop}

Let $G$ be a finitely generated discrete group, and let $\ell$ be the
word length associated to a fixed finite and symmetric set of
generators. Of the many equivalent definitions of amenabilitiy we
shall work with {\it Reiter's condition\/} -- $G$ is {\it amenable\/}
if for every $\varepsilon>0$ and for every $R>0$ there exists a
finitely supported
$\xi\in\ell^1(G)$ such that $\xi\geqslant 0$, $\|\xi\|=1$ and
\begin{equation}
\label{amen}
   \ell(g)\leqslant R \Rightarrow \| g\cdot \xi - \xi \| <\varepsilon,
\end{equation}
where the action of $G$ on $\ell^1(G)$ is defined by $g\cdot \xi(h)=\xi(g^{-1}h).$

Our main tool to prove Proposition~\ref{prop:nqe} is the use of Property A, a weak form of amenability, introduced
by Yu in the context of the Baum-Connes conjecture in topology~\cite{Yu}.

Let $X$ be a discrete metric space of {\it bounded geometry\/} -- that
is, the number of points in a ball of fixed radius is bounded, the
bound depending only on the radius of the ball and not on its center.
Of the many equivalent definitions of Property $A$ we choose
the one most closely related to Reiter's condition --  $X$ has {\it  Property
$A$\/} if for every $\varepsilon>0$ and $R>0$ there exists an $S>0$ and
for each $x\in X$ a function $\xi_x\in \ell^1(X)$ such that $\xi_x\geqslant
0$, $\|\xi_x\|=1$ and \begin{align*}
  d(x,y)\leqslant R &\Rightarrow \| \xi_x -\xi_y \| < \varepsilon, \\
  \xi_x(y) \neq 0 &\Rightarrow d(x,y)\leqslant S.
\end{align*}
The analogy with amenability being clear, we say that a metric space
having Property $A$ is {\it coarsely amenable\/} whereas one not
having Property $A$ is {\it coarsely non-amenable}.

Finally, a metric space $X$ is the {\it coarse union\/} of its
subspaces $X_n$ if $X=\sqcup X_n$ (disjoint union), and if
$d(X_n,X_m)\to \infty$ as $n+m\to\infty$. If the $X_n$ are metric
spaces each having finite diameter, then there exists a metric space
$X$ which is the coarse union of (isometric copies of) the
$X_n$. Further, any
two such unions are coarsely equivalent.  Moreover, if $Y$ is the coarse
union of the $Y_n$ then $X$ and $Y$ are coarsely equivalent
when the $X_n$ and $Y_n$ are uniformly coarsely equivalent.

We require the following slight
generalization of \cite[Prop. 11.39]{roe}.  We include a proof which is
both different from other proofs in the literature and
convenient for our result.

\begin{prop}\label{prop:cu}
  Let $G$ be a finitely generated amenable group.
  Every quotient of $G$ is amenable; the coarse union of any family of finite
  quotients of $G$ is coarsely amenable.
\end{prop}

\begin{proof}
  Let $H$ be a quotient of $G$ and identify $H$ with a set of cosets
  $\{\, gN \,\}$.  Fix a finite and symmetric set of
  generators for $G$ and equip $G$ with the associated word length;
  equip $H$ with the word length associated to the induced generators.
  With these conventions
  \begin{equation*}
    \ell_H(x)\leqslant R\ \Longleftrightarrow\
       \text{$\exists \; g\in x$ such that $\ell_G(g)\leqslant R$}
  \end{equation*}
  and, in particular, the map $G\twoheadrightarrow H$ is contractive. Given
  $\varepsilon>0$ and $R>0$, obtain $\xi\in\ell^1(G)$ as in
  (\ref{amen}).  Define
  \begin{equation}
\label{eta}
    \eta(x) = \sum_{g\in x} \xi(g),
  \end{equation}
  so that $\eta\geqslant 0$ and $\| \eta\|=1$.  Further, when $z\in H$ has
  length at most $R$ we obtain $g\in G$ of length at most $R$ such
  that $z=gN$.  We then calculate
  \begin{equation*}
  \| z\cdot \eta -\eta \| = \sum_{x\in H} | \eta(g^{-1}x)-\eta(x)|
       \leqslant \sum_{x\in H} \sum_{h\in x} | \xi(g^{-1}h)-\xi(h)|
       = \| g\cdot \xi-\xi\| <\varepsilon.
\end{equation*}
We conclude that $H$ is amenable.

When dealing with a coarse union the essential observation is that, in
the previous argument, if $\xi$ is supported on the elements of length
at most $S$ then the same is true of $\eta$. Thus, let $H_n$ be a
family of quotients of $G$, each equipped with a length function as
above, and let $X$ be a coarse union of the $H_n$. Given
$\varepsilon>0$ and $R>0$ proceed as above -- obtain a Reiter function
$\xi$ for $G$ and define $\eta_n$ as in (\ref{eta}). For $x\in X$
define
\begin{equation*}
  \xi_x = \begin{cases}
        \chi_N, &\text{$x\in H_n$, $n\leqslant N$}\\
        x\cdot \eta_n, & \text{$x\in H_n$, $n>N$},
  \end{cases}
\end{equation*}
where $N$ is chosen large enough so that for $n>N$ the distance
between $H_n$ and any other $H_m$ is at least $R$; $\chi_N$ is the
characteristic function of $H_1\cup\dots\cup H_N$.  Finally, choose
$S$ larger than the diameter of $H_1\cup\dots\cup H_N$ and large
enough so that $\xi$ is supported on elements of length at most $S$ in
$G$.

The required properties are easily verified.
\end{proof}

\begin{proof}[Proof of Proposition~\ref{prop:nqe}]
  The  iterated squares are proper characteristic subgroups of the free group, hence,
by Levi's theorem~\cite[Ch.I, Prop.~3.3]{ls},    they have trivial intersection,
  $\cap \mathbb{F}_2^{(2)\ldots(2)}= \{\, 1 \,\}.$
     Thus, the coarse union of the metric
  spaces $\Gamma_n$ is an example of a coarsely non-amenable {\it box
    space\/}. See \cite[Def. 11.24 and Prop. 11.39]{roe}.  (This statement is the converse of the
  previous proposition, and can also be proved by modifying the above
  argument.) This finishes the proof as coarse amenability is invariant under coarse equivalence,
  see, for example,~\cite[Prop. 4.2]{Tu}.
\end{proof}

We conclude with two remarks.  First, we have used a very crude
invariant from coarse geometry to distinguish towers constructed
from the sequence of iterated squares (\ref{squares}) from those
constructed beginning with an amenable group in (\ref{setup}) -- the
former are coarsely non-amenable while the latter are coarsely
amenable.  More refined invariants would be needed to establish the
existence of coarsely inequivalent towers constructed as in
(\ref{setup}) from a given non-amenable group.

Second, our construction involving the iterated squares
(\ref{squares}) is particular to the free group.  It would be
interesting to remove the hypothesis of `largeness' from our theorem.

\bibliographystyle{plain}


\section{Introduction}
Understanding the phenomena close to quantum critical points (QCPs)
is one of the major challenges in condensed matter physics. A QCP
emerges when the characteristic temperature of a system, e.g. the
N\'{e}el temperature, is suppressed to zero by variation of a
non-thermal parameter $r$ like pressure, magnetic field or doping.
Upon pressure-tuning a system towards a QCP, the Gr\"{u}neisen ratio
$\Gamma$ is expected to display singular behavior, namely a
divergence towards zero-temperature, while by tuning a system by
magnetic field towards quantum criticality $\Gamma_{\rm H}$, as
defined below, diverges.\cite{ZhuLijun:UnidGp,GarstM:SigctG} Both
Gr\"{u}neisen parameters consist of ratios of the control-parameter
and temperature derivatives of the entropy, which is accumulated
close to the QCP ($S$, $\alpha$, $C$ and $M$ denote the entropy,
volume thermal expansion, specific heat and magnetization,
respectively):
\begin{eqnarray}
&&\Gamma=-\frac{1}{V_mT}\frac{(\partial S/\partial p)_T}{(\partial S/\partial T)_p}= \frac{\alpha}{C_p}\\
&&\Gamma_{\rm H}=-\frac{1}{T}\frac{(\partial S/\partial
H)_T}{(\partial S/\partial T)_H}=-\frac{(\partial M/\partial
T)}{C_H}=\frac{1}{T}\left.\frac{dT}{dH}\right|_S \nonumber
\end{eqnarray}

Indeed, both $\Gamma$ and $\Gamma_H$ have been experimentally shown
to diverge in materials near
QCPs.\cite{KuchlerR:DivtGr,KuchlerR:Grurdq,Lorenz07,ISI:000263389500045}
Since the magnetic field, in contrast to pressure or doping, can
easily be tuned continuously in-situ, the magnetic Gr\"uneisen ratio
is ideally suited to investigate the nature of QCPs. So far, the
magnetic Gr\"uneisen ratio close to QCPs has been determined only
indirectly by calculation from magnetization and specific heat
data\cite{ISI:000263389500045}. Note, that $\Gamma_H$ can be
determined by a single measurement of the magnetocaloric (MCE)
effect, $(dT/dH)_S$, as shown above. Very recently, the
field-dependence of the entropy close to a QCP has been calculated
from highly non-adiabatic (quasi isothermal) $T(H)$ traces,
separating the effect of heat flow from/to the bath and the
magnetocaloric effect of the sample, using independent measurements
of the thermal conductance between sample and
bath.\cite{ARostScience09} However, a more direct measurement of the
MCE under truly adiabatic conditions is highly desired for
comparison with theoretical scenarios. Below, we show that the
alternating-field technique is suitable for this purpose.
Previously, the spin-flip transition in insulating GdVO$_4$ has been
investigated using such technique down to 0.5~K.\cite{Fischer} In
this study not much care has been taken on the absolute values of
the MCE and eddy current heating plays no role, since the material
is an electrical insulator. Below, we report the application of the
low-frequency alternating-field technique to clean metals and down
to mK temperatures. We carefully investigate the effect of eddy
current heating. Furthermore, we prove by comparison with
magnetization and specific heat measurements, that accurate absolute
values of the adiabatic MCE under are obtained.

\section{Experimental Setup}
The definition of the MCE implies that a change of the sample
temperature $\Delta T$ caused by a small change of the magnetic
field $\Delta H$ under adiabatic conditions needs to be detected.
Below, we demonstrate how the MCE could be detected directly by
applying a low-frequency alternating field, $H_0\sin(2\pi ft)$, and
detecting the temperature oscillation of the sample with the same
frequency, $T_0\sin[2\pi f(t+\delta t)]$ ($\delta t$ is a small
delay due to the thermal resistance between the sample and
thermometer).

Our setup for generating the alternating field is shown in Fig.~1. A
SR830 Lock-In amplifier (Stanford Research Systems) is used to
produce a current with low frequency (typically 0.01$\sim$0.1\,Hz,
the choice of the frequency will be discussed later). The current is
amplified by a Kepco bipolar power amplifier up to 8\,A and sent
through a modulation coil (8\,A current corresponds to 25\,mT) which
is inserted to the main superconducting magnet and cooled within the
liquid 4-He. The current is determined by measuring the voltage
across a shunt resistance (4\,m$\Omega$) with a Keithley 195A
multimeter. The modulation field is superimposed to the desired
static magnetic field of interest which is generated by a large
superconducting magnet. Note, that a finite static magnetic field is
required since the MCE depends on the magnetization of the system
(cf. eq. 1).

 \begin{figure}[t]
 \includegraphics{fieldgeneration.eps}
 \caption{\label{}Generation and control of the alternating magnetic
  field superposed to the large magnetic field of the
  main superconducting magnet.}
 \end{figure}

Figure 2 shows the setup to detect the sample temperature under
quasi-adiabatic conditions. The sample is glued on a thin sapphire
plate and is cooled or heated through the weak thermal link (thin
CuNi wire). If $\tau$ denotes the thermal relaxation through the
link, the frequency of the field modulation must be adjusted such
that $1/f\ll\tau$ in order to obtain quasi-adiabatic conditions.
Such adjustment is easily be performed at several different
temperatures by taking temperature traces with differing frequency
and calculation of the MCE (see below). If the MCE is frequency
independent, the measurements are performed under quasi-adiabatic
conditions. The RuO$_2$ resistive thermometer is glued on the sample
and its leads are made from superconducting filaments to avoid an
additional heat leak. The sample temperature is measured using a
LS370 resistance bridge (Lake Shore). The heat sink is weakly
coupled to the mixing chamber and its temperature is PID controlled
with a heater within the field-compensated region of the large
superconducting magnet. In a typical measurement of the temperature
dependence of the MCE, we slowly sweep the temperature of the heat
sink leading to a slow sweep of the average sample temperature. Due
to the alternating magnetic field, the sample temperature oscillates
through its average value. The amplitude of this oscillation gives
the MEC $\Gamma_H$, while the average value determines $T$.

 \begin{figure}
 \includegraphics{setup.eps}
 \caption{\label{}Schematic view of the platform used for
 alternating field magnetocaloric effect measurements under
 quasi-adiabatic conditions.}
 \end{figure}

\section{Application to strongly correlated materials}

Figure~3 shows a typical example of a temperature sweep during an
alternating field MCE measurement. We have chosen YbRh$_2$Si$_2$ as
this is a system close to a QCP, at which a pronounced divergence of
the magnetic Gr\"uneisen ratio has been
observed.\cite{ISI:000263389500045} The data discussed here are
obtained on a slightly Fe-doped single crystal of 20.1\,mg
mass in a shape of plate with a thickness of 0.25mm (residual resistivity ratio $\sim$12\cite{future_pub}). A relatively large piece of single crystal was chosen to determine precisely its magnetization and specific heat, which we will discuss later. The field was applied within the easy magnetic plane of the
system (parallel to the plane of the sample plate), which is perpendicular to the tetragonal $c$-axis. The
static field is set to 50\,mT which is the critical field for
undoped YbRh$_2$Si$_2$ and the alternating field has an amplitude of
4.5\,mT with frequency of 0.1\,Hz. A clear oscillation of the sample
temperature due to the MCE is resolved. The inset of Fig.~3(a)
displays the temperature variation on a larger time scale such that
individual oscillations could not be resolved. Here the thickness of
the line is a measure of the oscillation amplitude. Indeed the
thickness of the $T(t)$ trace increases as temperature is decreased,
indicating a larger temperature oscillation $dT/dH$ with decreasing
temperature. This is a clear evidence of quantum criticality in this
material, since $dT/dH$ is expected to decrease linearly with
decreasing temperature for a non-critical material with constant
magnetic Gr\"uneisen ratio $\Gamma_H=const$.

 \begin{figure}
 \includegraphics{dtdh2.eps}
 \caption{\label{}(Color online) Sample temperature (a) and magnetic field (b) vs time for a measurement
 on doped YbRh$_2$Si$_2$. The red solid and dotted lines show a fit to $T(t)=a_0+a_1 t+a_2 t^2+T_0\sin(2\pi f(t-\delta t))$ and its background term, $a_0+a_1 t+a_2 t^2$, respectively. The inset displays the temperature variation of the sample over
  a longer time scale.}
 \end{figure}

In order to quantify the MCE at a given average temperature
$T_{av}$, we chose time intervals of five periods (5$\times 1/f$
seconds in time) and determine the amplitude of the temperature
oscillation by fitting the temperature trace to $T(t)=a_0+a_1 t+a_2
t^2+T_0\sin(2\pi f(t-\delta t))$ (red solid line in Fig.~3(a)),
where $f$ is the frequency of the field modulation
$H(t)=H_{av}+H_0\sin(2\pi ft)$. The second order polynomial in time
($a_0+a_1 t+a_2 t^2$, cf. red dotted line in Fig.~3(a)) is added to
the fitting function to take into account the background temperature
drift due to the slow sweep of the heat sink temperature. The
amplitude of the temperature oscillation $T_0$ divided by the field
modulation amplitude $H_0$ gives the MCE. This value is associated
to the average temperature $T_{av}$ given by the background
temperature at the center of the time interval, $T_{av}=a_0+a_1
(t_0+2.5/f)+a_2 (t_0+2.5/f)^2$, where $t_0$ is the time of the first
data point for fitting. The magnetic Gr\"uneisen ratio is then
calculated by $\Gamma_H(T_{av})=1/T_{av}\times T_0/H_0$. By shifting
the fitting interval across the temperature trace a dense and
continuous curve of $\Gamma_H$ is obtained, as shown in Fig. 4. In
this measurement from 3~K down to 0.3~K, the temperature of the heat
sink was slowly decreased, using PID control, and below 0.3\,K the
heater for the heat sink was turned off, since the thermal
conductance of the thermal link is very small at low temperatures,
and thus, the decrease of the average sample temperature is very
slow. We have also performed a measurement with half value of the
frequency ($f$=0.05\,Hz) and the same $H_0$ in the entire
temperature range from 3 to 60\,mK, using the same cooling rate of
the heat sink, and confirmed that the result is exactly the same.
This proves that the measurement is taken under effective adiabatic
conditions ($1/f\ll\tau$; i.e., heat flow between sample and heat
sink within one period is negligible.). The data are compared with
$\Gamma_H(T)$ calculated using Eq. 1 from magnetization and specific
heat data~\cite{future_pub} taken at the same field of 0.05~T and
measured using the same piece of single crystal. The excellent
agreement between the two independent ways to obtain the magnetic
Gr\"uneisen ratio prove that this technique is best-suited to
quantitatively obtain the true adiabatic MCE. Furthermore, the
comparison shows the extraordinarily high resolution of the
alternating field technique.

 \begin{figure}
 \includegraphics{compare_yrs_fe.eps}
 \caption{\label{}Magnetic Gr\"{u}neisen ratio for Yb(Rh$_{0.95}$Fe$_{0.05}$)$_2$Si$_2$ as
  a function of temperature. The open circles are determined by the alternating field
   technique while the solid squares are
   calculated from magnetization and specific heat data using $-(dM/dT)/C$.}
 \end{figure}

At last we discuss the problem of Joule heating caused by the field
modulation. If the sample is very clean and the mean free path of
the electrons is rather high, the frequency and the amplitude of
alternating field have to be low to avoid Joule heating due to eddy
currents. Figure~5 shows data obtained in the normal state of a
high-quality single crystal of the clean heavy-fermion
superconductor CeCoIn$_5$ in magnetic field applied parallel to
c-axis. Note here that the material has a very large electronic mean
free path of 810\AA\cite{movshovich:prl-01} and alternating-field is
applied perpendicular to the plane of a platelet sample
(1.3$\times$1.1$\times$0.33\,mm$^3$) with a twice larger frequency
of 0.2\,Hz than the one used for Fe-doped YbRh$_2$Si$_2$. This
condition causes a severe problem with eddy current heating. As
clearly seen in Fig.~5, the temperature of the sample oscillates
with twice the frequency of the field modulation which is due to
Joule heating. Eddy currents induced by the variation of a magnetic
field are proportional to the time derivative of the field,
$dH/dt=2\pi f H_0\cos(2\pi ft)$. Since eddy current heating is
proportional to the square of the current, $(dH/dt)^2\sim
f^2H_0^2[1+\cos(2\pi 2ft)]$, it results in an oscillation with $2f$
frequency. Note that the heating effect has two components; constant
and 2$f$-oscillating ones. Figure~5 also shows that as soon as the
alternating field is turned on, the average sample temperature
raises immediately due to the constant heating component, followed
by the 2$f$-oscillation. We have also verified that the amplitude of
$2f$ component increases rapidly with amplitude $H_0$ and frequency
of the alternating field, whereas the $1f$ component due to the
intrinsic MCE remains constant. The eddy current heating can very
effectively be suppressed by a reduction of the cross-section of the
sample perpendicular to the field axis and as well as by reducing
the frequency and/or amplitude of the field modulation.

 \begin{figure}
 \includegraphics{2omega2.eps}
 \caption{\label{} Sample temperature (a) and magnetic field (b) vs time for a high-quality
 single crystal of CeCoIn$_5$ for $H\parallel$[001]. Note that the temperature oscillation has twice
 the frequency of the field modulation.}
 \end{figure}

At an elevated temperature, where the electrical conductivity is
smaller, and therefore, the Eddy current heating is less severe, a
mixture of intrinsic 1$f$ and extrinsic 2$f$ oscillations can be
found, as shown in Fig.~6. The temperature as a function of time is
well fitted by a function with two oscillating components,
$T(t)=a_0+a_1 t+a_2 t^2+T_0\sin(2\pi f(t-\delta
t))+T_0^{2f}\cos(2\pi 2f(t-\delta t_{2f}))$ (red line in Fig.~6). We
deduced $T_0$, using several different frequencies but the same
$H_0$ and found a constant $T_0$. This indicates that the obtained
$T_0$ is still valid under the influence of eddy current heating.

 \begin{figure}
 \includegraphics{mix_osci.eps}
 \caption{\label{}(Color online) Mixture of 1$f$- and 2$f$-oscillations of sample temperature of CeCoIn$_5$ under the same condition as Fig.~5. The red line is a fit to $T(t)=a_0+a_1 t+a_2 t^2+T_0\sin(2\pi f(t-\delta t))+T_0^{2f}\cos(2\pi 2f(t-\delta t_{2f}))$}
 \end{figure}

\section{Summary}
We developed a high-resolution low-temperature alternating-field
technique for quasi-adiabatic measurements of the MCE and magnetic
Gr\"uneisen ratio for clean metallic samples down to mK temperatures
and in large magnetic fields. The performance check has revealed
excellent agreement with the magnetic Gr\"uneisen ratio calculated
by magnetization and specific heat measurements. The main advantages
of the new technique compared to the latter are that it is much
faster and easier and that its resolution is much higher. The
capability of continuous temperature sweeps makes possible thorough
investigations of the magnetic Gr\"uneisen ratio in the entire
$H$-$T$ phase space, which is of particular interest for systems
close to magnetic-field induced QCPs.

\section{Acknowledgements}
We are grateful to M. Brando, M. Lang, R. Movshovic, A. Steppke, K.
Winzer and B. Wolf for valuable conversations. We thank H. S. Jeevan
and E. D. Bauer for providing the
Yb(Rh$_{0.95}$Fe$_{0.05}$)$_2$Si$_2$ and CeCoIn$_5$ single crystals,
respectively. Work supported by the German Science Foundation
through SFB 602 TPA19 and FOR 960 (Quantum phase transitions).