text:
\section{Introduction}

Supernova (SN) shocks are generally thought to be the main accelerator of the bulk of
Galactic Cosmic Rays (GCR). Indeed, the power of GCR in the Milky Way is estimated to
10$^{41}$ ergs/s, corresponding to 10-20 \% of the kinetic power of Galactic
supernovae (assuming canonical values of 2 SN per century, each one releasing 10$^{51}$ ergs
of kinetic energy).

The site of the acceleration of GCR remains debatable today, despite more than five decades
of theoretical and observational studies (e.g. Strong et al. 2007
and references therein). Over the years, it has been
suggested that GCR are accelerated in 1) SN remnants (either by the forward or the reverse shock
or both), 2) the interstellar medium (ISM), 3) the winds of massive stars, 4) the interiors
of superbubbles, excavated by the massive star winds and the subsequent SN explosions of an OB
association.

Each one of the proposed sites has its own advantages and shortcomings, regarding the
energetics and/or the composition of accelerated matter. For instance, it has been argued that the hot,
low density environment of a superbubble minmizes radiative losses of SN shocks and
energy losses of accelerated particles, thus allowing the latter to reach substantial energies,
up to the "knee"of the GCR spectrum (e.g. Parizot et al. 2004).
On the other hand, reverse SN shocks running into the SN interior carry insufficient energy
to explain the bulk GCR energetics (Ramaty et al. 1997). Moreover, they should accelerate
$^{59}$Ni, a product of explosive nucleosynthesis which is unstable to e$^-$-capture (with
a lifetime of 10$^4$ yr) and which has not been detected in GCR (Wiedenbeck et al. 1999), while a rapid acceleration
would render it practically stable and thus detectable.

It was realized early on that the {\it elemental} composition of GCR differs significantly
from the one of the ISM. Those differences may provide valuable information
on the origin of GCR particles and, perhaps, on the site - and even the mechanism - of acceleration.
Volatiles behave  differently from refractories: the former display a mass-dependent enrichment
with respect to H, which reaches a factor of 10 for the heaviest of them; the latter are all
overabundant (w.r.t. H) by a factor of 20, while C and O display intermediate overabundances,
by factors of 9 and 5, respectively (e.g. Wiedenbeck 2007 and references therein).

This complex pattern is now thought to result not from ionization effects (as suggested in Cass\'e and
Goret 1978, and further developped by Meyer 1985) but rather from effects related to elemental
condensation temperature (Meyer et al. 1997): refractories are locked in dust grains, which are
sputtered by repeated SN shocks and the released ions are easily picked-up and accelerated (Ellison et al. 1997).
This, quite elaborate, scheme, which builds on earlier ideas by e.g. Cesarsky and Bibring (1981),
accounts quantitatively for most of the observed features of GCR source composition; still, it
leaves unanswered the key issue about the acceleration site of GCR (and how it affects the
composition of accelerated matter).

The most conspicuous feature of GCR source composition   is undoubtely the high isotopic $^{22}$Ne/$^{20}$Ne \ ratio.
Its  value was measured since the late 1970ies (Garcia-Munoz, Simpson and Wefel 1979,
Wiedenbeck and Greiner 1981). The most accurate  measurement today, obtained from
analysis of the CRIS instrument,  leads  to a best estimate  (Binns et al. 2008)
of 0.387 $\pm$ 0.007 (statistical) $\pm$ 0.022 (systematic).
This is 5.3$\pm$0.3 times the value of the ($^{22}$Ne/$^{20}$Ne)$_{\odot}$ \ ratio in the solar wind.
Contrary to the case of the elemental source GCR abundances, which may be affected by various
physico-chemical factors (first ionization potential, condensation temperature, etc.)
isotopic ratios can only be affected by nucleosynthetic processes and thus provide crucial
information on the origin of cosmic ray particles. It should be noticed that up to now there
is now clear evidence for any other GCR isotopic ratio to differ from solar, with the potential exception of
$^{58}$Fe/$^{56}$Fe, which is estimated to 1.5$\pm$0.3 times solar (Binns et al. 2008).

Soon after the discovery of the anomalous GCR $^{22}$Ne/$^{20}$Ne \ ratio, Cass\'e and Paul (1982) suggested that
it could be explained by a mixture of $\sim$2\% of material from the wind of a WC star to 98\% of material
with standard composition. In early He-burning, $^{14}$N (produced through the CNO cycle in the previous
H-burning phase) is transformed almost totally in $^{22}$Ne \ through
$^{14}$N($\alpha,\gamma$)$^{18}$F($\beta^+$)$^{18}$O($\alpha,\gamma$)$^{22}$Ne. He-burning products (like $^{12}$C  and
$^{22}$Ne) are expelled by the stellar winds of  massive stars during their WC phase.
The observed GCR $^{22}$Ne/$^{20}$Ne \ ratio is obtained by assuming dilution of WC material with matter of standard composition.
Subsequent studies put the aforementionned idea on a quantitative basis, with the use of detailed models
of the evolution and nucleosynthesis of massive, mass losing stars (Maeder 1983, Prantzos 1984, Meyer 1985, Prantzos et al. 1987).
In those studies, the acceleration site of GCR was considered as  decoupled from the nucleosynthesis site,
and unrelated to the fraction of admixtured WC material.

Higdon and Lingenfelter (2003) evaluated quantitatively the $^{22}$Ne/$^{20}$Ne \ ratio within a superbubble, created by the
collective action of stellar winds and SN shockwaves. They adopted  stellar wind yields for $^{20}$Ne \ and $^{22}$Ne \ from the models
of Schaller et al. (1992) and SN yields from the models of Woosley and Weaver  (1995). They found that the $^{22}$Ne/$^{20}$Ne \ ratio
in the superbubble decreases with time (since $^{22}$Ne \ from the winds dominates the evolution of $^{22}$Ne/$^{20}$Ne \ at early times) and that
its time average value is compatible with the GCR source $^{22}$Ne/$^{20}$Ne \ inferred from observations. In a subsequent paper, Lingenfelter
and Higdon (2007) recognised that the Schaller et al. (1992) yields of $^{22}$Ne \ were highly overestimated\footnote{The reason
was the excessively high mass loss rates adopted in that work.} and, consequently, "... new detailed calculations of the
expected GCR isotopic ratio are called for...", but they did not attempt such a re-evaluation. In the meantime,
Binns et al. (2005), using updated wind yields of massive stars with rotation (from the Geneva group, see Sec. 2.2), found
good agreement between the observed $^{22}$Ne/$^{20}$Ne \ ratio and an admixture of $\sim$20\% material from WR stars  with 80\% material
of standard composition. According to  Binns et al. (2008), since WR stars are evolutionary products of OB stars,
such an agreement  "...suggests that OB associations within superbubbles are the likely source of at least a substantial fraction
of GCR".

Howewer,  theoretical studies in the past 10 years  are based mostly on the  paradigm of GCR being accelerated in SN remnants, not in
superbubbles, e.g. \cite{Ptuskin05, Berezhko06,Berezhko09,Ptuskin10,Caprioli10,Schure10,Ellison11} \
and references therein. The kinetic energy of the bulk motion of the forward shock of the SN explosion is converted to GCR  energy
through diffusive shock acceleration. The process is highly non-linear and involves the dynamical reaction of both the accelerated
particles and of the magnetic field on the system. Those studies usually take into account the fact that the SN explosion often
occurs within the cavity excavated in the interstellar medium (ISM) by the wind of the massive star prior to the explosion
(\cite{Biermann01}); however, the structure of the circumstellar environment in that case is quite complex
and simplified models are used for its description. Although \cite{Caprioli11b} considered the composition
of GCR (H,He, CNO, MgSiAl, Fe)  resulting from such an acceleration site, none of those studies considered the $^{22}$Ne/$^{20}$Ne \ ratio.

In this work we study the $^{22}$Ne/$^{20}$Ne \ ratio of GCR accelerated by the forward shocks of   SN  explosions, as they run
 through the presupernova winds of  massive stars and through the interstellar medium. We consider the whole mass spectrum
 of massive stars (from $\sim$10 to 120  M$_{\odot}$), including  stars with either small or large mass losses prior to their
 explosions. We consider stellar properties (masses of winds, ejecta, yields etc.) from recent models with mass loss
 and or without rotation (from Hirschi et al. 2005 and Limongi and Chieffi 2006, respectively), the former
 having larger $^{22}$Ne \ enhancements in their winds.
 We adopt a simplified prescription (suggested in Ptuskin and  Zirakashvili 2005 and reformulated in Caprioli 2011)
 to describe the structure of the  circumstellar medium at the time of the explosion and we consider that GCR start being accelerated
 in the Sedov-Taylor (ST)  phase of the SN remnant (see e.g. Ptuskin et al. 2010).  By requiring  the resulting
 IMF averaged $^{22}$Ne/$^{20}$Ne \ ratio to equal the observed one  $R_{Obs}$=($^{22}$Ne/$^{20}$Ne)$_{GCR}$/($^{22}$Ne/$^{20}$Ne)$_{\odot}$=5.3$\pm$0.3
 we are able to constrain the forward shock velocity to  values $>$1900 km/s for rotating stars
 (and to $>$2400 km/s for non rotating ones), i.e.
 we find that GCR are accelerated during the early ST phase, lasting for  a few 10$^ 3$ yr. Assuming that 10\% of the SN
 kinetic energy is converted to GCR, we find that during the acceleration  period a few particles out of a million encountered
 by the forward shock are accelerated.  Finally, we reassess the superbubble paradigm for the origin of GCR, by evaluating
 consistently  the $^{22}$Ne/$^{20}$Ne \ ratio with the aforementioned  stellar yields. We find that it can not be as high as observed, unless
 some extremely favorable assumptions are made (only the early period of the superbubble lifetime considered, no gas left over
 from the formation of the OB association). We conclude that superbubbles cannot be at the origin of the bulk of GCR.

 The plan of the paper is as follows. In Sec. 2 we present the general "set-up" of our model: the adopted stellar models (Sec. 2.2)
 and  their wind yields (Sec. 2.3), the description of the circumstellar environment (Sec. 2.4) and the evolution of the forward shock in the
 ST phase (Sec. 2.5). In Sec. 3 we present our results for the (time-dependent) composition of the accelerated particles, the limits
 imposed on the shock velocity by the observed $^{22}$Ne/$^{20}$Ne \ ratio and
 the efficiency of the particle acceleration. Finally, in Sec. 4 we explore the $^{22}$Ne/$^{20}$Ne \  ratio of GCR,  assumed to be accelerated inside
 a superbubble, and we show that it cannot match the oberved one (unless extreme assumptions are made). The results are
 summarized in Sec. 5.

\section{A toy model for the composition of CR accelerated in massive star winds}

 \subsection{The set-up}
 \label{sub:Toy}

The method adopted here in order to calculate the composition of matter accelerated by a single
SN explosion is schematically illustrated in Fig.~\ref{Fig:SNStruct}. At the end of its
life and at the time of its SN explosion, a star of initial mass
$M_{*}$ \  is left with a mass $M_{Exp}$,   surrounded by a circumstellar shell of mass
$M_{Wind}$=$M_{*}$-$M_{Exp}$,  which has been lost through stellar wind during its prior hydrostatic evolution
After the SN explosion, a mass of ejecta $M_{Ej}$=$M_{Exp}$-$M_{Rem}$ \ (where $M_{Rem}$ \ is the mass of the compact
remmant, neutron star or black hole) expands first within the shell  of mass $M_{Wind}$ \ and then
in the ISM, with the forward shock having initial velocity $\upsilon_0=\sqrt{2 E_0/M_{Ej}}$,
 where $E_0$ is the kinetic energy of the SN explosion.

In the case of stars ending their lives as WR stars, the wind contains both the original
(nuclearly unprocessed) envelope of mass $M_{Env}$, and nuclearly processed layers of mass $M_{Proc}$, enriched
in products of H-burning, and in some cases of He-burning as well. For those stars,
$M_{Wind}$=$M_{Proc}$+$M_{Env}$ \ and $M_{Proc}$ \ is calculated as the difference between the mass of the
nuclearly processed core $M_{HeC}$ \ and the mass at the explosion: $M_{Proc}$=$M_{HeC}$-$M_{Exp}$.
For lower mass stars, exploding as red supergiants, the wind composition results essentially from the 1st dredge-up, i.e. it is a mixture of H-burning products from the
stellar core with the original envelope composition (i.e. mass loss has not uncovered the He-core
at the time of the explosion). The limit between the two classes of stars
depends on their initial mass, mass loss rate and rotational velocity and it is rather poorly
known at present: in general, in models with  no rotation  stars with $M_{*}$$>$32-35 M$_{\odot}$ \
become WR stars e.g. (Heger et al. 2002), while in models
with rotation that limit may be as low as 22 M$_{\odot}$ \ (Meynet and Maeder 2000).

The first phase of the supernova remnant ("free expansion")
takes place at shock velocity $\upsilon\sim const.$ and ends when a mass
$M_{S1}$$\sim$$M_{Ej}$ \ has been swept up in front of the shock wave, at which point the
ST phase sets in. Following Ptuskin et al. (2010), we assume that efficient GCR acceleration
starts at this time, where the situation is energetically most favorable.
In our baseline model we shall consider constant acceleration efficiency ;
time-dependent efficiency
of particle acceleration is the subject of current researches (see Ellison and Bykov 2011,
Drury 2011, and references therein) and will be briefly discussed in Sec. 3.3.

The ST phase proceeds adiabatically, i.e. at $\sim$constant energy and
with decreasing velocity, until
the temperature of the gas engulfed by the shock front  drops to levels allowing a significant fraction
(about 50\%) of the remaining energy to be radiated away. At that time, an amount
of matter $M_{S2}$$>>$$M_{Ej}$ \ has been swept up and the shock enters the "snow-plow" phase. At this
point - and, perhaps, even earlier, during the ST phase - the forward shock is
too weak to accelerate particles to GCR energies any more.

In the aforementioned scenario, GCR are accelerated from a pool of particles with composition
characteristic of the mass $M_{Wind}$ \ early on. Depending on the initial stellar mass, this
composition may be rich in products of H- (and He-) burning. It is progressively diluted
with ambient (first wind - with normal $^{22}$Ne - and then interstellar) gas and at the end of the ST
phase it ressembles closely the one of the ISM.  The  GCR source composition observed on
Earth should correspond to the  average composition
between the early ST phase and some later evolutionary stage of the remnant,
and should result from the whole mass spectrum
of exploding stars, i.e. it should be averaged over a stellar initial mass function (IMF).

 \begin{figure}
   \centering
   \includegraphics[width=8.5cm]{fig1.eps}
   \caption{Schematic representation of a supernova exploding in  the wind  of its
   parent star. The star, of initial mass M$_{*}$ explodes with a mass
   M$_{Exp}$, i.e. it has lost a mass M$_{wind}$=M$_{*}$-M$_{Exp}$. The
   most massive stars become WR stars and  their wind expels not only the H-envelope (of mass M$_{Env}$,
   with composition similar    to the one of the ISM) but also nuclearly processed layers,
   of mass $M_{Proc}$, i.e. M$_{Wind}$=M$_{Env}$+M$_{Proc}$, where M$_{Proc}$=M$_{HeC}$-M$_{Exp}$ and
   M$_{HeC}$ is the mass of the (H-exhausted) He-core. The star leaves a remnant
   (neutron star or black hole) of mass M$_{Rem}$; the mass ejected in the SN explosion
   is M$_{Ej}$=M$_{Exp}$-M$_{Rem}$. Efficient GCR acceleration presumably starts at the beginning of
   the ST phase, when a mass M$_{S1}\sim$M$_{Ej}$ is swept up
   in front of the SN shock wave (which is indicated by arrows).

    }
            \label{Fig:SNStruct}
   \end{figure}

\subsection{Properties of mass losing stars}
\label{sub:StarProp}

   \begin{figure}
   \centering
   \includegraphics[width=0.495\textwidth]{fig2.eps}
   \caption{Input data adopted in this work, for solar metallicity massive stars with mass loss.{\it Left}:
   non rotating models from Limongi and Chieffi (2006); {\it Right:} Rotating models from Hirschi, Meynet and Maeder (2005).
    {\it Top:} Mass of the (H-exhausted) He core
   $M_{HeC}$, of the star at explosion $M_{Exp}$ and of the compact remnant $M_{Rem}$ as a function
    of the initial stellar mass $M_{\bf *}$. {\it Bottom:}
   derived quantities are the mass of the wind $M_{Wind}$=$M_{\bf *}$-$M_{Exp}$, the mass of the
   nuclearly processed layer $M_{Proc}$=$M_{HeC}$-$M_{Rem}$ and of the unprocessed ($^{22}$Ne-normal) envelope
   $M_{Env}$=$M_{Wind}$-$M_{Proc}$ and the mass of the SN ejecta $M_{Ej}$=$M_{Exp}$-$M_{Rem}$
   In all cases, curves are smoothly interpolated between model results.
    }
\label{Fig:SNmasses}
   \end{figure}

We adopt two sets of stellar models in this work. They are calculated
for stars of solar metallicity, and in both cases  the solar mixture of Anders and Grevesse
(1979) is adopted. The corresponding metallicity is Z$_{\odot}$=0.019, substantially larger than
more recent values (Lodders 2003, Asplund et al. 2010) and this difference results
in particular from the reduction in the past decade of the solar abundances of C, N, O and Ne, which
are key elements for the purpose of this work. For obvious consistency reasons, we keep here the
Anders and Grevesse (1979) values, when comparing  our results for GCR to solar ones.

The first set of stellar models
is the one of the Frascati group (Limongi and Chieffi 2006, hereafter LC06). It concerns 15 model
stars between 11 and 120 M$_{\odot}$ \  with mass loss but no rotation. The model includes
all stages of hydrostatic nuclear burning
and simulates the final stellar explosion by imparting an initial velocity
 to a mass coordinate of 1 M$_{\odot}$ \ (i.e. well inside the Fe core of the stars);
the  mass cut (the limit separating the ejecta of mass $M_{Ej}$ \ from the compact remnant of mass $M_{Rem}$) is chosen such as  0.1 M$_{\odot}$ \ of $^{56}$Ni is ejected by the explosion.

The various masses involved in the "toy model" of Sec.~\ref{sub:Toy} and Fig.~\ref{Fig:SNStruct}
are provided in Table 2 of LC06 and are displayed in Fig.~\ref{Fig:SNmasses} ({\it top left}) of this work,
whereas derived quantities are displayed in Fig.~\ref{Fig:SNmasses} ({\it bottom left}).
It can be seen that stars with $M_{*}$$<$20 M$_{\odot}$ \ have lost a negligible
amount of mass prior to the explosion ($M_{Wind}$$<$$M_{Ej}$) and the ST phase starts within the
ambient ISM.
Stars with $M_{*}$$>$30 M$_{\odot}$ \ have $M_{HeC}$$>$$M_{Exp}$, i.e. they explode as WR stars, having
expelled nuclearly processed layers in their winds. But only for the most massive stars
($M_{*}$$>$60 M$_{\odot}$) one has $M_{Ej}$$<$$M_{Proc}$, i.e. in the beginning of the ST phase the shock wave still
expands into material with composition reflecting the one of the nuclearly processed core; for lower
mass stars (30$<$$M_{*}$/M$_{\odot}$$<$60) the ST phase starts when the shock wave encounters material with mixed composition from the core and the envelope, i.e. less enriched in H- and He- burning products.
An interesting feature of those models is that the mass of the ejecta $M_{Ej}$$\sim$10 M$_{\odot}$, is similar
for all of them, leading to similar properties(duration, swept-up mass)
 for their corresponding ST phases.

The second set of stellar models is the one of the Geneva group, calculated by
Hirschi, Meynet and Maeder (2005, hereafter HMM05) and includes both mass loss and rotation for stars with mass between 12 and 60 M$_{\odot}$; it has been complemented with results for stars of 85 and 120 M$_{\odot}$, kindly provided by G. Meynet (private communication). The initial rotational velocity is
$v_{Rot}$=300 km/s on the ZAMS, corrrsponding to an average velocity of 220 km/s on the main sequence,
i.e. close to the average observed value. Non-rotating models with the same physical ingredients (for convection, mass loss, etc.)
have also been calculated, for comparison.
The evolution has been calculated to the end of Si-burning but the final SN explosion was not considered. Instead, an empirical prescription was used to evaluate the mass of the compact remnant.
The masses of our "toy model", as given in Table 2 of HMM05, are displayed in Fig.~\ref{Fig:SNmasses}.
There is an important difference with respect to Table 2 of HMM05: they provide the mass
of He-core at the time of the explosion and for the most massive stars this coincides
with $M_{Exp}$ \ (i.e. the mass left to the star at explosion
is smaller than the maximum extent of the He-core); however, we are interested at the true value of
$M_{HeC}$ \ (since this will determine how much mass of processed material the shock wave will encounter) and
that value is obtained through the detailed results of HMM05 (displaying the wind composition
as a function of time - or of mass left).

Comparing the results of LC06 and HMM05 one sees that rotation increases the mass loss ($M_{Wind}$ \
larger in HMM05), thus leaving the star with a smaller mass at explosion ($M_{Exp}$ \ smaller in HMM05).
Rotation also increases the size of nuclearly processed regions ($M_{HeC}$ \ larger in HMM05), since
matter is rotationally mixed outwards to larger distances
than achieved through convection. In turn, this leads to  larger amounts of processed material
$M_{Proc}$ \ for the HMM05 models.

The aforementioned features of rotating vs non-rotating models, which are explained in detail in e.g.
Maeder and Meynet (2000) are crucial in understanding the differences in the corresponding
wind yields of the stars.

\subsection{The wind composition of massive stars}

LC06 provided (private communication)  {\it yields} $y_i$ of all stable nuclear species, from H to Ge,
included in their models and ejected through the winds of the stars, up to the moment of the explosion.
HMM05 provide (Table 3 in their paper) the {\it net yields} $y_{n,i}$ of the winds of
their models for stable species from
$^3$He to $^{23}$Na, from which the yields can be recovered through
\begin{equation}
y_i(M_*) = y_{n,i}(M_*) + {M_{Wind}} {\rm X_{\odot,i}}
\end{equation}
where the adopted solar values X$_{\odot,i}$ are displayed in Table 1 of HMM05.

The wind yields of a few selected species appear in Fig.~\ref{Fig:windmasses}, for
 the non-rotating models of LC06 and
for both the non-rotating and the rotating models of HMM05. It can be seen that, in general,
there is excellent agreement between the results for non-rotating models  of LC06 and HMM05,  for stars
up to 40 M$_{\odot}$. Their results differ only for the 60 M$_{\odot}$ \ model (and presumably for higher masses
as well) and only for the cases of the He-burning products $^{12}$C, $^{16}$O and $^{22}$Ne.
Since both HMM05 and LC06 use the same prescriptions for mass loss, the reason of that discrepancy could
be the use of a small amount of overshooting in the case of HMM05.

   \begin{figure}
   \centering
   \includegraphics[width=0.49\textwidth]{fig3.eps}
   \caption{Total masses of selected isotopes in the winds of the massive stars, from the
   two sets of models adopted in this work. {Solid curves}  are yields from HMM05 (with rotation)
     and {dotted curves}  from LC06 (no rotation),
      in both cases interpolated between model results.
    The points correspond to results of the HMM05 models with no rotation and,  in general,
     agree well with LC06 results except for the 60 M$_{\odot}$ \ star.
   }
         \label{Fig:windmasses}
   \end{figure}

Rotation has a twofold effect on stellar yields: it increases the size of the nuclearly processed layers
(since it mixes material further than convection alone) and reduces the escape velocity in the stellar equator,
allowing larger amounts of mass to be ejected in the wind. Both effects enhance the wind yields up to some
mass limit; above it, the wind has removed so much mass, that less material is left in the star
to be processed in subsequent stages of the evolution, thus reducing the corresponding yields. This
is the case, for instance, with the He-burning products $^{12}$C \ and $^{22}$Ne, the yields of which decrease above
$\sim$60 M$_{\odot}$ \ in the rotating HMM05 models (see Fig. \ref{Fig:windmasses} and HMM05 for details).

In the following we assume  that the wind interaction
with the ISM has not substantially changed the wind stratification:
 the forward shock will first encounter
the innermost wind layers, containing processed material in the case of  the most massive stars;
later it will encounter the outer wind layers (containing mostly the initial composition), before
running into the ISM.  Fig. 4  displays  the {\it mass integrated composition of the wind} (for a few key metals),
 as encountered
by the forward shock, moving outwards from $M_{Exp}$, in the case of  two rotating model stars of
 25 M$_{\odot}$ \ and 60 M$_{\odot}$ \ (from HMM05). The quantity
 \begin{equation}
 m_i(M) \ = \ \int_{M_{Exp}}^{M_*}X_{Wind,i}(M) \ dM
 \end{equation}
is displayed as a function of mass coordinate $M$, X$_{wind,i}(M)$ being the mass fraction of isotope $i$
in the wind of star of mass $M_{*}$. Obviously, one has $m_i$($M_{Exp}$)=0. and  $m_i$($M_{*}$)=$y_i$($M_{*}$), i.e.
at crossing the last (outermost) wind layer, the forward shock has encountered the totality of the yield $y_i$($M_{*}$).

An inspection of Fig.~\ref{Fig:wind_Struct} shows that:

- in the case of the 25 M$_{\odot}$ \ star (no He-burning products encountered by the shock wave),
 the innermost layers contain most of the $^{14}$N \ produced by the CNO cycle (its mass increases
 slowly in the outermost layers), and little of $^{12}$C \ and $^{16}$O, which are depleted by the CNO cycle
  (their mass increases rapidly in the outermost layers); $^{20}$Ne \ and $^{22}$Ne \ are little affected by H-burning
  and their integrated wind mass increases in a way intermediate between the  cases of $^{14}$N \ and $^{12}$C.

- in the case of the 60 M$_{\odot}$ \ star (He-burning, then H-burning products encountered by the shock wave),
the quasi-totality of the $^{12}$C\ and $^{22}$Ne \ produced by He-burning are encountered in the inner layers
(inside $\sim$25 M$_{\odot}$), while essentially no $^{14}$N \ is left in that region; the majority of $^{14}$N \ is found in the
region 25$<M$/M$_{\odot}$$<$50  while substantial  $^{16}$O \ is found in the unprocessed envelope (beyond 50 M$_{\odot}$).
Finally, $^{20}$Ne \ is little affected everywhere ($\sim$const. mass fraction) and its integrated mass rises regularly with
the wind mass of the star.

   \begin{figure}
   \centering
 \includegraphics[angle=-90,width=0.49\textwidth]{fig4.eps}
 \caption{ Composition encountered by the forward shock wave as it moves
 (from left to right) through the winds of 25 M$_{\odot}$ \ (left) and 60 M$_{\odot}$ (right) rotating stars.
 The shock wave starts at $M_{Exp}$ =10 M$_{\odot}$ and 14.8 M$_{\odot}$, respectively, and the ST phase
 sets in when a mass M$_{S1}$=$M_{Exp}$ \ has been swept up and, presumably, well mixed.
 The quantities displayed correspond to the mass integrated composition of the wind,
 starting at $M_{Exp}$.
 }
         \label{Fig:wind_Struct}
   \end{figure}

If detailed information (i.e. wind composition as function of wind mass coordinate) is available,
the composition of the material swept-up by the forward shock can be calculated in a consistent way. In case only the
wind yields $y_i$ and the various characteristic masses of Fig. 2 are available, one may adopt the following approximations:

(i)  Stars never reaching the WR stage ($M_{*}$$<$22 M$_{\odot}$ \ for rotating models and $M_{*}$$<$34 M$_{\odot}$ for non rotating ones)
have their envelopes fully mixed after the 1st dredge-up episode in the red supergiant phase; the wind composition
is then simply obtained as:
\begin{equation}
X_{wind,i}(M_*)  \ = \  {{y_i(M_*)}\over{M_{Wind}(M_*)}}
\end{equation}

(ii) Stars reaching the WR stage eject first their unprocessed envelope of mass $M_{Env}$ \ and
composition X$_{env,i}$=X$_{ISM,i}$
and then the innermost processed layers of mass $M_{Proc}$ \ and composition X$_{PROC,i}$ such
that
\begin{equation}
y_i(M_*) \ = \ X_{Env,i} M_{env}(M_*) \ + \  X_{Proc,i} M_{Proc}(M_*)
\end{equation}
which allows one to derive  X$_{Proc,i}$, i.e. an average abundance in the processed layers, and obtain thus
an average wind composition profile

\begin{equation}
X_{Wind,i}(M)  \ = \  X_{Proc,i}(M_*)  \ \ {for } M_{Exp}<M< M_{HeC} \\
\end{equation}
\begin{equation}
X_{Wind,i}(M)  \ = \  X_{Env,i}(M_*) \ \ {\rm  for} M_{HeC}<M< M_*
\end{equation}

The method outlined here turns out to provide a good approximation for the detailed composition
profiles of the HMM05 models and has been used here for both sets of LC06 and HMM05 models,
for consistency reasons. Taking ito account all the uncertainties of the models (prescriptions for mass loss,
convective and rotational mixing etc.),  we consider the aforementioned approximation
as fairly satisfactory for the purpose of this work.

\subsection{The circumstellar environment of massive stars}

According to the ideas outlined in Sec. 2.1 particle acceleration to GCR energies
starts when a mass $M_{S1}\sim$$M_{Ej}$ \ has been swept-up by the forward shock.
The composition of that material, at the very beginning of particle acceleration,
is the one of the material located at mass coordinate $A_1$=$M_{Exp}$ + $M_{Ej}$ \
and constitutes obviously an extreme for the abundances of the corresponding elements
(upper limit for those produced in the stellar interior and lower limit for those destroyed, like H).
In the following we shall assume that at the time of the explosion  the stellar wind
has fully kept the stratification of its various layers, as they were progressively leaving the
stellar surface.

 The circumstellar region hit by the forward shock has a complex structure, depending on the
 properties of the exploded star. That structure has been explored in some detail with hydrodynamical models
 by Garcia-Segura et al. (1996a,b) for non-rotating stars. In the case study of a 35 M$_{\odot}$ \ star, they find
 that the rapid wind of the O star excavates a large bubble (radius $\sim$36 pc) of low density (10$^{-3}$ cm$^{-3}$).
 Inside it propagates a dense ($\sim$1 cm$^{-3}$), slow wind - released during the red supergiant
phase - which occupies the innermost few pc (depending on its assumed velocity, of the order of a few tens of km/s).
The subsequent fast ($\sim$10$^3$ km/s) WR wind compresses the RSG wind, and most of the mass of the latter is
found within a thin shell.

 The aforementioned study illustrates the complexity of the situation, but its results can hardly be
 generalized to the whole mass range of
 massive stars (for instance, lower mass stars will not display the fast WR wind).
 In fact, its results cannot even be safely used for the 35 M$_{\odot}$ \ star, since they
 depend so critically on the adopted parameters of the model (mass loss rates and wind velocities for the
 various stages). And they certainly fail to describe the situation for rotating stars, which display
 slow but intense (and not radially symmetric) mass losses on the Main sequence.

In view of these uncertainties, we adopt here a simplified prescription for the
structure of the circumstellar bubble, assuming spherical symmetry
in all cases. We assume that the winds have excavated a bubble of {\it mean density}
$n_0$=0.01 cm$^{-3}$ and, consequently of radius
 \begin{equation}
 R_B \ = \  \left(\frac{3 \ M_W}{4 \pi \ \rho_0}\right)^{1/3}
 \end{equation}
 with $\rho_0$=$n_0 \ m_p$,  $m_p$ being the proton mass.
Inside the bubble, the density profile is $\rho(r) \propto r^{-2}$, i.e. it corresponds to a
 steady stellar wind with mass loss rate ${\dot M_W}$ and velocity $v_W$, which is given by
  \begin{equation}
 \rho(r) \ = \  \frac{\dot M_W}{4 \pi v_W r^2}
 \end{equation}
Our choice of $\rho_0$ automatically fixes the  $\rho(r) $ profile and corresponds to a combination
of ${\dot M_W}$ and $v_W$ values. Obviously, one has
\begin{equation}
M_W \ = \ \int_0^{R_W} \ \rho(r) \ dr
\end{equation}
Outside $R_W$ we assume an ISM with constant density
$\rho_{ISM}$=1 $m_p$ cm$^{-3}$.
Our approach is similar to Caprioli (2011), but we do not consider here
the more complicate case of a WR wind overtaking a RSG wind.

\subsection{Evolution during the Sedov-Taylor (ST) phase}

We follow the propagation of the forward shock first through the wind bubble and then through the ISM
 with a simple model presented in Ptuskin and Zirakashvili (2005) and, in a more concise form, in
 Caprioli (2011, his eqs.
   3.4 to 3.9).  We start the calculation from the "free-expansion"(ejecta dominated) phase, where
 the swept-up mass is smaller than $M_{Ej}$ \ and which can be described by self-similar
 analytical solutions. In the subsequent ST  phase (swept up mass $<$ $M_{Ej}$),  the model is based on the
 "thin shell"approximation (e.g. Ostriker and McKee 1988)
 which assumes that the swept up mass is concentrated in a thin shell behind the shock.
 We solve numerically the time-dependent equations
 for continuity of mass, energy and momentum, to recover shock radius and velocity as a function of time.
 Unlike Caprioli (2011) we assume full adiabaticity
 in the ST phase, i.e. we do not take into account the 10-20\% energy losses of the shock
 through CR acceleration, which would reduce the shock velocity ($\propto E_0^{0.5}$) by less than 10\%.
 The mass swept up in the ST phase inside shock radius $R_S$ is
 \begin{equation}
 M_S(<R_S) \ = \  M_{EJ} \ + \ 4  \pi \int_{R_{S1}}^{R_{S}} \rho(r) r^2  dr
 \end{equation}
We follow the evolution all the way through the  ST  phase, which
ends when a significant fraction of the energy of the cooling remnant is radiated
away (through recombination emission); for a solar mixture this occurs at time
\begin{equation}
t_{S2} \ \sim 4.4 \ 10^4 \ {\rm yr} \  \left({{E_0}\over{10^{51}{\rm erg}}}\right)^{2/9}
\left({{n_{ISM}}\over{{\rm cm^3}}}\right)^{-5/9}
\end{equation}

In the framework of this simple model we are able to calculate the composition
of the material encountered (and presumably accelerated)
by the forward shock as a function of time, or of the swept-up mass:
indeed, the integrated mass of each element swept up by the
forward shock is given by an equation similar to Eq. 2
\begin{equation}
 m_i(M_S) \ = \ \int_{M_{exp}}^{M_S(<R_S)} X_{i}(M_S) \ dM_S
 \end{equation}
 where $X_{i}(M_S)$=X$_{wind,i}$ for $M_S<$$M_{*}$ \ (Sec. 2.3)
 and $X_{i}(M)$=X$_{ISM,i}$ for $M_S>$$M_{*}$ , i.e. when
 the shock propagates in the ISM; the upper limit in the integral is
 given by eq. (10).

Eq. 12 allows one to link the stellar model, i.e. the abundance profiles $X(M)$, to the evolution during the ST
phase  through Eq. 10,  and to the properties of the shock wave.
 Since the mass M$_{S2}$ swept-up in the end of the ST phase is much larger than
 the wind mass in all cases (a few 10$^3$ M$_{\odot}$, compared to $\sim$100 M$_{\odot}$ at most)
for the largest part of the ST phase the swept up material has ISM composition. In order
to obtain significant deviations from the   solar composition, such as the observed $^{22}$Ne/$^{20}$Ne \ ratio, one should assume that
{\it significant acceleration  occurs  only in the early ST phase}, when the
forward shock is stronger and its velocity higher.

 \section{Results}

\subsection{Composition of matter in the ST phase}

   \begin{figure}
   \centering
 \includegraphics[width=0.49\textwidth]{fig5.eps}
 \caption{ Evolution of the SN remnants of two rotating stars of initial masses
 20 M$_{\odot}$ \ ({\it dotted}) and  60 M$_{\odot}$ \ ({\it solid}), respectively. {\it Top}:
 velocity of the forward shock
 and  swept up mass. {\it Middle}: Shock radius and density profile before shock arrival
 (see text).
 {\it Bottom:} Composition of the matter swept up up to time $t$
 by the shock wave,
 for the $^{14}$N/$^{16}$O \ and $^{22}$Ne/$^{20}$Ne \ ratios. In all panels, the {\it thick } portions of the curves indicate
 the period of efficient particle acceleration, i.e. from the beginning of the ST phase and as far as
  $\upsilon_s>\upsilon_{min}$.
 The value of $\upsilon_{min}$ (=1900 km/s for the rotating star models in the figure) is chosen
  so that the IMF averaged theoretical ratio  of $^{22}$Ne/$^{20}$Ne \
 matches the observed one in GCR (see text, Eq. (13)  and Fig.~\ref{Fig:ModelsVsObs}).
 }
         \label{Fig:GCRaccel}
   \end{figure}

Fig. \ref{Fig:GCRaccel} (top) displays the evolution of the velocity $\upsilon_S$
and radius $R_S$ of the forward shock and of the mass $M_S(<R_S)$ swept up
by it, for the cases of a 20 and a 60 M$_{\odot}$ \ rotating star, respectively.
 The density of the unperturbed ISM is taken to be 1 cm$^{-3}$ in all cases.

The similarity of the curves for  $\upsilon_S$, $R_S$ and $M_S(<R_S)$ for the 20 and
60 M$_{\odot}$ \ stars simply reflects the self-similarity of the ST solution. The small differences
in the early ST phase are due to the difference of the ejected mass $M_{Ej}$ \ in the two stars
(6 vs 10 M$_{\odot}$, see Fig. 2), since the other parameters ($E_0$ and $n_{ISM}$) are the same.

The bottom panel of Fig.  \ref{Fig:GCRaccel} displays the evolution of the $^{14}$N/$^{16}$O \ and  $^{22}$Ne/$^{20}$Ne  \
ratios (i.e. the ratios of the corresponding masses of Eq. 12  for each isotope),
 expressed in units of the corresponding solar values. In both cases, the forward shock
first encounters layers with low $^{16}$O \ and high $^{14}$N, resulting in a high $^{14}$N/$^{16}$O \ ratio (in the 60 M$_{\odot}$ \ star,
$^{14}$N \ is depleted from He-burning in the innermost layers, resulting in a slightly smaller
$^{14}$N/$^{16}$O \ ratio than in the 20 M$_{\odot}$ \ case).
Subsequently, in the 20 M$_{\odot}$ \ star the shock moves rapidly through the small remaining
stellar envelope (3.4 M$_{\odot}$) and starts propagating in the ISM, thus decreasing  rapidly its
mass-integrated $^{14}$N/$^{16}$O \ ratio.
In the 60 M$_{\odot}$ \ star, the shock runs through $\sim$20 M$_{\odot}$ \ of processed material, with
a high value of $^{14}$N/$^{16}$O, before getting to the ISM; the corresponding $^{14}$N/$^{16}$O \ ratio decreases
more slowly than in the 20 M$_{\odot}$ \ case.

The evolution of $^{22}$Ne/$^{20}$Ne \ is quite different in the two models. In the 20 M$_{\odot}$ \ star,
no He-burning products are encountered by the shock wave and the $^{22}$Ne/$^{20}$Ne \ ratio
has always its initial (solar) value. A high $^{22}$Ne/$^{20}$Ne \ ratio is initially
encountered in the processed layers
of the 60 M$_{\odot}$ \ star, which is progressively diluted as the shock moves outwards.

   \begin{figure}
   \centering
 \includegraphics[width=0.49\textwidth]{fig6.eps}
 \caption{Abundance ratios of various nuclear species in GCR source normalized to
 the corresponding solar ones,  as a function of the
 initial stellar mass. In all panels, {\it solid curves} correspond to models of
  HMM05 and {\it dotted curves} to models of LC06. {\it Upper,  thin curves}
  are for  GCR  accelerated at the beginning of the ST phase
  and {\it lower,  thick curves } for the time-average at the end of GCR acceleration.
  An average over a Salpeter IMF (and accounting
 for the swept up mass in each case) produces the vertical segments to the right, their top
 point corresponding to the beginning and the bottom one to the end of the GCR acceleration
 phase, respectively (also indicated by {\it filled squares}). These results are
 compared to GCR source abundance
 ratios as derived by ACE data (points at the extreme  right
 with error bars) in Binns et al. (2005). The most significant,
 unaffected by FIP, volatility etc., is the one of $^{22}$Ne/$^{20}$Ne. The end of the GCR
 acceleration phase is assumed to correspond to shock velocities $\upsilon_{min}$ such that
the time and IMF averaged theoretical ratio (squares) of $^{22}$Ne/$^{20}$Ne
 matches the observed one (see text). For the
 set up adopted here we find we find $\upsilon_{min}$=1900 km/s for HMM models and
 $\upsilon_{min}$=2400 km/s for LC06 models.}
         \label{Fig:ModelsVsObs}
   \end{figure}

\subsection{Composition of accelerated particles}

The composition curves of the bottom panel of Fig.  \ref{Fig:GCRaccel} give the time
(or mass) integrated composition encountered by the shock wave and, in consequence,
the composition of the particles that have been accelerated {\it up to that time}.
We assume here that {\it particles are accelerated to GCR energies with the same efficiency
for shock velocities
higher than some critical value $\upsilon_{min}$, which is the same for all stellar masses}.
We determine $\upsilon_{min}$ empirically by requiring that, when averaged over a stellar Initial mass function
(IMF) $\Phi(M_*)$,  the ratio
\begin{equation}
\frac{M_{22}}{M_{20}} = \frac{\int_{10 M_{\odot}}^{120 M_{\odot}} dM_* m_{22}(M_*) \Phi(M_*)}{\int_{10 M_{\odot}}^{120 M_{\odot}} dM_* m_{20}(M_*) \Phi(M_*)} =
R_{Obs} \frac{X_{22,\odot}}{X_{20,\odot}}
\end{equation}
where $R_{Obs}$=5.3$\pm$0.3  is the observationally determined source GCR ratio of $^{22}$Ne/$^{20}$Ne \ in solar units and $m_{22}(M_*)$ and $m_{20}(M_*)$
are calculated from Eq. 12 for stars of mass $M_{*}$ \ and for swept up masses $M_S(\upsilon>\upsilon_{min}$).

We adopt here a Salpeter IMF $\Phi(M_*) \propto M_*^{-X}$ with $X$=2.35.
The results of the procedure appear in Fig. \ref{Fig:ModelsVsObs} for a few selectedabundance ratios and for the models
of both HMM05 ({\it solid curves}) and LC06 ({\it dotted curves}). In all panels, the upper ({\it thin}) curves
correspond to the composition accelerated at the beginning of the ST phase (maximal possible deviations
from solar composition).  It can be seen that stars with mass $<$22 M$_{\odot}$ \ (for HMM05) and 32 M$_{\odot}$ \
(for LC06) display no He-burning products in their accelerated particles.
N is overabundant in lower stellar masses, due to the 1st dredge-up.

The lower ({\it thick}) curves  in all panels of Fig. 6 correspond to composition accelerated {\it up to the end of the
acceleration period} which is assumed to occur for a shock velocity $\upsilon_{min}$.
The corresponding IMF-averaged quantities (between 10 and 120 M$_{\odot}$)
are displayed on the right of the curves: their uppermost point
corresponds to the beginning of the ST phase and the lower one (also indicated with a {\it filled square})
to the end of the acceleration period, i.e. to  $\upsilon_{min}$. The value of $\upsilon_{min}$
is found to be $\sim$ 1900 km/s  for the HMM05 models with rotation and $\sim$20\% higher (2400 km/s)
for the LC06 models without rotation. The reason for that difference is, of course the fact that rotating models
have larger processed layers, requiring  more dilution with circumstellar material.

The top and middle panels of Fig. 6 display the corresponding ratios for $^{12}$C/$^{16}$O \ and $^{14}$N/$^{16}$O, respectively.
In both cases, the IMF averaged ratios are higher than the observed ones in GCR  (which are, in their turn,
higher than solar), by factors of 1.5-2. Unlike $^{22}$Ne/$^{20}$Ne, these are ratios of different elements, having
different atomic properties. GCR source abundances are known to be affected by atomic effects, e.g. First
Ionization Potential (FIP), or, perhaps more plausibly, volatility and mass/charge ratio
(see extensive discussion in Meyer et al. 1997). The analysis of such effects is beyond the scope
of this study. We simply notice here that observations indicate that
refractory elements are relatively more abundant in
GCR sources than volatiles. Meyer et al. (1997) attribute that to the fact that refractories are
locked up in grains, which are sputtered by the shock wave and the released ions are easily picked up
and accelerated. $^{16}$O \ being more refractory than both $^{14}$N \ and $^{12}$C, it is expected that its abundance
in GCR source will be enhanced by that effect, and the corresponding $^{12}$C/$^{16}$O \ and $^{14}$N/$^{16}$O \ ratios in GCR source
will be lower than predicted from stellar nucleosynthesis alone. The fact that in Fig. 6 we obtain
higher than observed ratios for $^{12}$C/$^{16}$O \ and $^{14}$N/$^{16}$O \ is encouraging in that respect, since the aforementioned atomic
effects would lower those ratios, hopefully to their observed values\footnote{One might think that a combined
analysis of the $^{12}$C/$^{16}$O \ and $^{14}$N/$^{16}$O \ ratios of Fig. 6 could help to constrain the atomic processes shaping the GCR source
abundances. However, the abundances of $^{12}$C \ and $^{16}$O \ are affected by the still uncertain  value
of the $^{12}$C($\alpha,\gamma$$^{16}$O \ reaction rate and are unsuitable for such a study.}.

  \begin{figure}
   \centering
 \includegraphics[width=0.49\textwidth]{fig7.eps}
 \caption{Graphical presentation of the results discussed in Sec. 3.3. Particle acceleration
 starts at the beginning of the ST phase, located at mass coordinate $A_1$=$M_{Exp}$ + $M_{Ej}$, i.e. when
 the forward shock (FS, {\it arrows}), launched at $M_{Exp}$,  has swept up a mass $M_{S1}$=$M_{Ej}$. Acceleration stops
 at mass coordinate $A_2$, corresponding in the case dicussed here to a shock velocity of 1900 km/s.
 The mass sampled by the FS between those two regions is $M_{Acc}$=$A_2$-$A_1$. For rotating stars with
 mass $M>$30 M$_{\odot}$, an increasing part of $M_{Acc}$ includes nuclearly processed material ({\it shaded aerea}),
 while for rotating stars with $M<$18 M$_{\odot}$, $M_{ACcc}$ contains only material of ISM (=solar) somposition.
 }
         \label{Fig:AccelvsM}
   \end{figure}

 \subsection{Efficiency of particle acceleration}

 The material of Sec. 3.1 and 3.2 is summarized in Fig. \ref{Fig:AccelvsM} for the rotating models of HMM05. The forward
 shock, launched at mass coordinate $M_{Exp}$, sweeps up a mass $M_{S1}\sim$$M_{Exp}$ \ and then
 starts accelerating particles, at mass coordinate $A_1$=$M_{Exp}$+$M_{S1}$, up to point $A_2$ (where
 its velocity becomes $\upsilon_{min}$).
 For rotating stars with $M<$15 M$_{\odot}$, $A_1$ lies beyond the stellar surface and only ISM is accelerated.
 In stars with 15$<$M$_*$/M$_{\odot}$$<$25, $A_1$ lies beyond the processed/mixed
 interior $M_{HeC}$ \ and the forward shock accelerates first envelope and
 then ISM material. For stars above 35 M$_{\odot}$, the shock first
 accelerates processed material (hatched aerea), then the - $^{22}$Ne \ normal - envelope and then ISM. Finally, for stars with M$>$70 M$_{\odot}$, particle acceleration ends when the shock is still within the massive stellar envelope.

   The mass of circumstellar material
from which particles are accelerated is  $M_{ACC}=A_2 - A_1$ and lies in the range 30-40 M$_{\odot}$.
The mass of particles that have been accelerated is  given by
\begin{equation}
m_{AP} = \sum N_i A_i m_P
\end{equation}
where $N_i$ is the total number of nuclei of species $i$, $A_i$ the corresponding mass number and $m_P$ the proton mass. The number $N_i$, or rather the product $N_i A_i$ can be determined by noticing that
the total energy carried by those accelerated particles is a fraction $f$ of the kinetic energy $E_0$ of the supernova
 \begin{equation}
  \sum N_i A_i \int_0^{\infty} E \ Q(E) \ dE \ = \ f \ E_{0}
 \end{equation}
 where $Q(E)$ is the spectrum of accelerated particles ($A_i$ appears on the left side of Eq. 15 because
  energies are expressed in energy units {\it per nucleon}). The efficiency $f$ of conversion of $E_0$ to
  accelerated particles {\it escaping the supernova} is estimated to $f\sim$0.1, while another
 0.1-0.2 goes to the acceleration of particles which are finally trapped in the SN (Ellison and Bykov 2011).
  The particle spectrum is often described by a power law in momentum

\begin{equation}
Q(E) \ \propto \ \frac{p^{-s}}{\beta} {\rm exp}(-E/E_0)
\end{equation}
where $\beta=v/c$ is the velocity expressed as a
fraction of the light velocity,
$p$ is the particle momentum per nucleon,
the factor $s$ is usually 2$<s<$3 (in the case
of strong shocks) and $E_0$ is a cut-off energy, here taken to be 3 TeV (the results are insensitive
to much higher values).

  Assuming that the spectrum of accelerated particles escaping the SN is given by Eq. 16, one finds
  the efficiency with which particles are accelerated from
  the shocked circumstellar medium, through Eq. 14 to 16:
 \begin{equation}
 W = \frac{m_{AP}}{M_{ACC}}
 \end{equation}
 and for rotating stars it is found to be in the range 3-6 10$^{-6}$ i.e. a few particles out of a million encountered ones are accelerated by the forward shock to GCR energies. In the case of non-rotating stars, the energetics is the same, but the swept-up mass is smaller (by a factor of two, on average, in order to get the observed $^{22}$Ne/$^{20}$Ne \ ratio) and the corresponding efficiency is $W\sim$10$^{-5}$. These estimates constitute
only a gross average, since the efficiency of particle acceleration may depend on several factors, not considered here,
like e.g.  the density of the circumstellar medium or the shock velocity - through a smoothly varying function
$f(v_S)$  instead of the Heavyside function $f(>\upsilon_{min})$=1 and $f(<\upsilon_{min})$=0  considered here -
or the shock radius at the time of acceleration, since particles may subsequently suffer adiabatic cooling before
escaping (e.g. Drury 2011, Bykov and Ellison 2011).
Notice that some of  those effects may have opposite time dependencies.  For instance,
particles accelerated earlier on (at higher shock velocities and presumably with higher efficiencies)  are expected
to suffer more from adiabatic cooling (because they are produced at smaller radii) . These effects require a much more
thorough investigation. {As a first step towards that direction, we also tested the case  where the efficiency of shock acceleration varies
with shock velocity as  $f \propto \upsilon^2$ (Drury 2011). In that case,  material with high   $^{22}$Ne/$^{20}$Ne \ is efficiently accelerated in the inner layers
of the most massive stars. In order to obtain again the observed GCR source $^{22}$Ne/$^{20}$Ne \ ratio, one has to dilute the mixture by allowing acceleration for shock
velocities lower then the reference values found above: we thus find values of $\upsilon_{min}$=1600 km/s for the HMM05 yields and 2150 km/s
for the CL06 yields; the corresponding overall acceleration efficiencies increase then by $\sim$40\% from the reference values given above,
ranging from 4-7 10$^{-6}$ for HMM05 yields to 1-1.4 10$^{-5}$ for CL06 yields.
The reason for obtaining such a small difference (only a few hundreds of km/s) between the non-realistic reference case and the -
perhaps more realistic - case of velocity dependent efficiency, is the adopted unified treatment, which considers acceleration
in both low mass and high mass supernova: the former accelerate almost pure ISM with low $^{22}$Ne/$^{20}$Ne \ and the latter almost pure wind material
with high $^{22}$Ne/$^{20}$Ne. Allowing for lower values of $\upsilon_{min}$ involves a considerably larger amount of ISM processed by
low mass supernovae, since the lower acceleration efficiency ($f \propto \upsilon^2$) is more than compensated  by the  $n_{ISM} 4 \pi r^ 2 dr$
factor.

As approximate as they may be, the results obtained here  through the constrain of the observed GCR source $^{22}$Ne/$^{20}$Ne \ ratio,
indicate clearly that acceleration has to occur for shock velocities larger than $\sim$1500 km/s, and they
may help to improve our understanding of particle acceleration in SN remnants.

\section{GCR cannot be accelerated (mainly) in superbubbles}

The idea that CGR are accelerated mainly in superbubbles has been suggested by Kafatos et al. (1981) and
reassessed by  Higdon et al. (1998).
Massive stars are mainly formed in OB associations and the winds of the most massive of them
initiate the formation of superbubbles, while the subsequent SN explosions power the expansion
of those superbubbles for a few 10$^7$ years. Higdon et al. (1998) argued that
the environment of such superbubbles, enriched with the ejecta of stellar winds and core collapse
SN explosions, provides a composition that compares favourably to the inferred GCR source composition.

One potential problem with that idea is that $^{59}$Ni, a well known product of SN nucleosynthesis,
is absent from GCR arriving on Earth (Wiedenbeck et al. 1999). Since $^{59}$Ni is unstable to electron capure, with a lifetime of
$\sim$10$^5$ y, Higdon et al. (1998) argued that SN explosions occur within the superbubble with a
sufficiently low frequency (less than one SN every 3 10$^5$ y), as to allow for the decay of $^{59}$Ni
between two SN explosions. Prantzos (2005) pointed out that energetic stellar winds can also accelerate particles
and they are not intermittent (like SN) but occur continuously in superbubbles; in that case $^{59}$Ni is
continuously  accelerated to GCR energies and, being unable to capture an electron, it becomes effectively stable
and it should be detectable in GCR. Binns et al. (2008) counter-argued that the period of energetic WR winds
represents a small (albeit not negligible) fraction in the early lifetime of an OB association and therefore only
little (and presumably undetectable) amounts of $^{59}$Ni would be accelerated, thus saving the "superbubble paradigm".

In the meantime, Higdon and Lingenfelter (2003) evaluated the isotopic composition of $^{22}$Ne/$^{20}$Ne \ - the critical abundance anomaly
in GCR -  expected in a superbubble, based on (a very heterogeneous set of)
then available yields of WR stars and SN. They found that the GCR $^{22}$Ne/$^{20}$Ne \ ratio
"...can be easily understood as the result of GCR accelerated primarily out of superbubbles with a mean metallicity
Z$_{SB}$=2.7$\pm$0.4 Z$_{\odot}$..." and that this result "... provides strong, additional evidence for a superbubble origin of GCR".
In a subsequent paper, Linfenfelter and Higdon (2007) explored quantitatively that senario for other GCR abundances, finding
again that results compare favourably to observations.

Despite their apparent sophistication, the aforementioned arguments for a superbubble origin of GCR miss a simple point:
{\it massive stars are the principal source of both $^{20}$Ne \ and $^{22}$Ne \ in the Universe}. As a result, the $^{22}$Ne/$^{20}$Ne \ ratio of a generation
of massive stars (integrated over an IMF and including WR winds and all kinds of SN ejecta) should be solar\footnote{Strictly speaking,
this concerns massive stars of roughly solar metallicity, such as those that contributed to the composition of the solar system; however,
metallicity has evolved very little in the past several billion years in the local volume of $\sim$ 1-2 kpc radius,
where most GCR originate.}. Since  $^{22}$Ne \ is overabundant with respect to $^{20}$Ne \ only in the massive star winds
(and not in the SN ejecta or the ISM),  a popular  exercice in the past 30 years  or so consisted in evaluating the mixing ratio
of WR ejecta with SN ejecta (or with average ISM), as to obtain the observed $^{22}$Ne/$^{20}$Ne \ ratio in GCR.
(see Introduction and references therein).

However,  it is expected that a properly
weighted mixture of the ejecta of massive stars in a superbubble (i.e. including stellar winds and SN ejecta and folded with a
stellar IMF) would produce a solar $^{22}$Ne/$^{20}$Ne \ ratio. In this section, the same exercise is repeated with the yields
of CL05 and HMM05 and by taking into account the corresponding stellar lifetimes. In Fig. \ref{Fig:Ne_yields}
we present again the adopted $^{20}$Ne \ and $^{22}$Ne \ yields of LC06 and HMM05  for the stellar winds (solid and dotted curves, i.e. the same
as in Fig. 3) and also the total  yields  (winds plus SN ejecta, points). The latter are from LC06 for stars with mass loss (but no rotation)
and from Woosley and Weaver (1995, WW95) for stars without mass loss; they are in good overall agreement, except for $^{20}$Ne \ in stars
around 20 M$_{\odot}$.  For the SN ejecta, we consider only stars with M$<$40 M$_{\odot}$,
i.e. we assume that more massive stars eject their $^{20}$Ne \ and $^{22}$Ne \ through their winds and then form
black holes; in that case we maximize the $^{22}$Ne/$^{20}$Ne \ ratio expected from massive stars.

  \begin{figure}
   \centering
 \includegraphics[width=0.49\textwidth]{fig8.eps}
 \caption{Yields   of $^{20}$Ne \ ({\it top}) and $^{22}$Ne \ {\it bottom})   as a function of the
 initial stellar mass.  Wind yields  are  from  HMM05 ({\it solid curves}) and LC06 ({\it dotted curves}), respectively
 (same as in Fig. 3). Total yields are from LC06 ({\it filled squares}) and WW95 ({\it asterisks}).
 }
         \label{Fig:Ne_yields}
   \end{figure}

An inspection of Fig. \ref{Fig:Ne_yields} shows that the total (SN) $^{22}$Ne \ yields of M$<$40 M$_{\odot}$ \ stars are
comparable to the wind $^{22}$Ne \ yields of the most massive ($\sim$100 M$_{\odot}$) stars; however, the corresponding
total $^{20}$Ne \ yields of M$<$40 M$_{\odot}$ \ stars are
at least 10 times larger than the wind $^{20}$Ne \ yields of the most massive ($\sim$100 M$_{\odot}$) stars.
Thus, while the production of $^{22}$Ne \ receives a sizeable contribution from  the most massive
stellar winds (at least for the rotating stars),
the production of $^{20}$Ne \ is totally dominated by the SN ejecta of M$<$40 M$_{\odot}$ \ stars. In a coeval stellar
population, like the one expected in an OB association or a superbubble, the $^{22}$Ne/$^{20}$Ne \ ratio
will evolve then from higher than solar to $\sim$solar values,
as the M$<$40 M$_{\odot}$ \ stars eject their core products a few Myr after
their more massive counterparts.

  \begin{figure}
   \centering
 \includegraphics[width=0.49\textwidth]{fig9.eps}
 \caption{Evolution of composition in a superbubble. {\it Top left:} Ejection rate
 of $^{20}$Ne ({\it dotted}) \ and $^{22}$Ne \ ({\it solid}) in M$_{\odot}$ /Myr per 1 M$_{\odot}$ \ of stars formed (assuming a Salpeter IMF between 0.1 and 120 M$_{\odot}$).
 {\it Top right:} Supernova rate per Myr.
 {\it Middle}: Total mass of $^{20}$Ne \ ({\it dotted})  and $^{22}$Ne \ ({\it solid}) accumulated in the superbubble, in case
 there was no initial gas there (M$_0$=0, {\it left}) and in the presence of an initial mass of gas equal to 1/5 of the
 formed stars (M$_0$=0.20,{\it right}); in the former case $^{22}$Ne \ dominates by mass early on, whereas in the latter the mass of $^{20}$Ne \ dominates
 from the very beginning. {\it Bottom}: Corresponding evolution of the O/H ({\it dashed}) and of the ($^{22}$Ne/$^{20}$Ne)$_{SB}$  ratio  ({\it dotted}) in the SB gas;  the
 ({\it solid} curve represents the ($^{22}$Ne/$^{20}$Ne )$_{CR}$ ratio of particles accelerated {\it up to} that time and has to be compared to the observed GCR source
 ratio ({\it dashed horizontal line}).
 }
         \label{Fig:SBevol}
   \end{figure}

We quantitative illustrate this effect in Fig. \ref{Fig:SBevol}, where we present the evolution of
a stellar population born at time $t$=0 with a Salpeter IMF and total mass of 1 M$_{\odot}$ \ (results can be directly
scaled to masses of  OB associations, while abundances and abundance ratios remain the same).
The top right panel displays the rate of SN explosions; for a total stellar
mass of 10$^4$ M$_{\odot}$, resulting in a $\sim$100 OB stars, there are about 3 SN/Myr, i.e. $\sim$1 SN every 300 000 yr, as
evaluated in Higdon et al. (1998) in order to avoid the problem of $^{59}$Ni acceleration.
The top left panel displays the ejection rates of $^{20}$Ne \ and $^{22}$Ne \ (in M$_{\odot}$/Myr per M$_{\odot}$ \ of stars formed) : $^{22}$Ne \
from stellar winds dominates only for a couple of Myr but  $\sim$5 Myr after the stellar formation, $^{20}$Ne \ from
M$<$40 M$_{\odot}$ explosions dominates.   Notice that, in order to maximise the $^{22}$Ne/$^{20}$Ne \ ratio, we adopted only the
yields of HMM05 for these calculations (i.e. corresponding to the solid curves in Fig. \ref{Fig:Ne_yields} with
only wind yields above M=40 M$_{\odot}$).

The instantaneous ($^{22}$Ne/$^{20}$Ne)$_{SB}$ \ ratio in the superbubble at time $t$ is the ratio of the amounts of $^{20}$Ne \ and $^{22}$Ne \ cumulated up to that time,
which are provided in the middle panels of Fig. \ref{Fig:SBevol}.
In order to evaluate those quantities, one needs to make an assumption about the amount
of pre-existing gas still left inside the superbubble. In the left middle panel of Fig. \ref{Fig:SBevol}
it is assumed that no gas is left after star formation: gas is exclusively supplied by the wind and explosive stellar ejecta.
In that case,  $^{22}$Ne \ dominates early on, but $^{20}$Ne \ takes over after a couple of  Myr. In the right middle panel it is assumed that an amount of gas with solar composition
and equal to 20\% of the mass of formed stars  remained in the superbubble; the stellar ejecta are diluted into it and the superbubble
composition is now dominated always by $^{20}$Ne.

The bottom panels of Fig. \ref{Fig:SBevol} display the corresponding ($^{22}$Ne/$^{20}$Ne)$_{\rm SB}$ \ ratios in the superbubble.
In the case of zero initial gas (left), the instantaneous  ($^{22}$Ne/$^{20}$Ne)$_{\rm SB}$ \ ratio  remains high -
from 18 to 5 times solar -  for a couple of Myr, but soon after stars of M$<$40 M$_{\odot}$ \ start dying its value goes rapidly to solar
(as expected from the nucleosynthesis argument presented in the beginning of this section).
The ($^{22}$Ne/$^{20}$Ne)$_{\rm CR}$ \ ($<t$) ratio of particles  accelerated to CR energies {\it up to time t}  is the time integral
of the instantaneous ($^{22}$Ne/$^{20}$Ne)$_{\rm SB}$ \ ratio in the superbubble
\begin{equation}
\left(\frac{^{22}{\rm Ne}}{^{20}{\rm Ne}}\right)_{\rm CR}(<t) \ = \ \frac{1}{t} \int_0^t  \left(\frac{^{22}{\rm Ne}}{^{20}{\rm Ne}}\right)_{\rm SB} dt'
\end{equation}
where it is assumed that the efficiency of particle acceleration does not vary with time.
 ($^{22}$Ne/$^{20}$Ne)$_{\rm CR}$ \ ($<t$) represents the mean CR composition accelerated by a superbubble which "operates" up to time $t$.

In the bottom left panel it is seen that the  ($^{22}$Ne/$^{20}$Ne)$_{\rm CR}$ \ ($<t$) ratio (in solar units) remains above the observed $R_{Obs}$=5.3$\pm$0.3  for
about 7 Myr, while it declines steadily after that, tending to an asymptotic value of $\sim$2, considerably smaller than observed.
The corresponding metallicity in the superbubble, expressed by the O/H ratio (dashed curve),  is quite high - about 10 times solar -
since it represents pure core collapse SN ejecta; such high metallicity values have never been reported for any astrophysical environment
(except SN remnants).

The bottom right panel of Fig. \ref{Fig:SBevol} displays the corresponding quantities in the more realistic case of a superbubble
endowed with some gas left over from the star formation. In that case, although metallicity  is still high (about 3 times solar
after 7 Myr), the  $^{22}$Ne/$^{20}$Ne \ ratios never get above 3 times solar, either in the superbubble or in the accelerated particles; in fact they
are close to  solar for the largest part of the superbubble lifetime.

The results displayed in Fig.  \ref{Fig:SBevol}  suggest that the observed $R_{Obs}$=5.3$\pm$0.3 ratio in GCR cannot be
obtained from material accelerated in a superbubble: even under the most favorable possible conditions, such as those
adopted here (yields from {\it rotating} massive stars, only wind yields considered above M=40 M$_{\odot}$, no gas left over from star
formation), the resulting $^{22}$Ne/$^{20}$Ne \ ratio remains at high values only in the early evolution  of the superbubble, during a small
fraction of its lifetime. It is interesting to notice in that respect that Binns et al. (2008) suggested that {\it no substantial particle acceleration must occur
during that early period}, in order to avoid the problem of $^{59}$Ni acceleration  from the WR winds; thus, saving
the "superbubble paradigm" for the origin of cosmic rays, requires that acceleration occurs only in the late superbubble evolution,
when the WR winds have essentially stopped and only SN inject intermittently their kinetic energy in the superbubble. However,  we have shown here
that during this period, the ($^{22}$Ne/$^{20}$Ne)$_{CR}$ \ ratio is considerably lower than observed. In fact, under realistic conditions (some initial gas left,
explosive yields also above 40 M$_{\odot}$ \ considered,  at least for some stars), the $^{22}$Ne/$^{20}$Ne \ ratio in a superbubble would never get to values as high
as observed, contrary to claims made in the literature in the past decade. It must be stressed that this conclusion
does not depend on detailed adopted yields, but on a simple argument, namely that the IMF averaged $^{22}$Ne/$^{20}$Ne \ ratio
in a superbubble has to be close to solar during the largest period of the superbubble lifetime.

\section{Summary}

In this work we explore some implications of the idea that cosmic rays
are accelerated by the forward shock in supernova remnants, during their ST phase.
We focus on the chemical composition resulting from such an acceleration and, in particular,
on the $^{22}$Ne/$^{20}$Ne \ ratio, which is the most characteristic feature of the  observed GCR source composition and
is unaffected by atomic effects. For that purpose, we adopt recent models of the nucleosynthesis and
evolution of massive stars with mass loss: those of LC06 with no rotation and those of HMM05 with rotation.

In Sec. 2 we present a detailed summary of the properties of those models and, in particular, of the
chemical composition of their winds, insisting on the fact that rotating models release more
$^{22}$Ne \ in their winds than non-rotating ones. We also present the adopted model for the evolution of a
SN remnant within a stellar wind, based on the ideas of Ptuskin and Zirakashvili (2005) and the equations summarized in Caprioli (2011).

In the framework of the simple model adopted here, we follow the time-dependent composition of  GCR,
accelerated by the forward shock as it runs through either the ISM (in the case
of stars with M$<$25-35 M$_{\odot}$, depending on rotation) or through the stellar wind (in the case of more massive stars).
In fact, during the largest part of the ST phase, the shock runs through ISM and encounters a $\sim$solar composition.
In order to reproduce the observed high value of ($^{22}$Ne/$^{20}$Ne)$_{CR}$    ($R_{Obs}$=5.3$\pm$0.3 in solar units)  after accounting for
the stellar IMF, we have then to assume that acceleration is efficient only during a short early period in the ST phase.
We chose to use the shock velocity as a criterion for efficient acceleration and, based on the aforementioned CR
composition argument, we find that shock velocities larger than $\sim$1900 km/s (for the rotating stellar models) or
2400 km/s ( for the non-rotating ones) are required. This result is obtained by
assuming a step function for the efficiency $f$ of particle acceleration
($f$=0 before the ST phase and after $\upsilon_{min}$, and $f$=1 between the two).
 For the - perhaps, more realistic -
 assumption of a velocity-dependent efficiency  $f\propto \upsilon^2$, we find slightly lower values for
 $\upsilon_{min}$ (1600 km/s for HMM05 yields and 2150 km/s for CL06 yields, respectively).
 In the framework of the adopted models, this corresponds to a circumstellar
mass of  several tens of M$_{\odot}$ \ encountered by the forward shock.
Assuming, furthermore, that 10\% of the SN kinetic energy is used in
acceleration of escaping cosmic rays with standard energy spectra, allows us to evaluate the efficiency of that acceleration:
we find that a few particles out of a million encountered by the forward shock  are accelerated to CR energies.
We also notice that this scheme of GCR acceleration does not suffer from problems related to the absence of unstable
$^{59}$Ni in observed GCR composition: this heavy nucleus is well inside the SN ejecta and is not reached by the
forward shock which accelerates only wind material and ISM.

The aforementioned scenario assumes that even the most massive stars, up to 120 M$_{\odot}$, develop
  strong forward shocks and accelerate the particles of their WR winds. For non-rotating stars, this
  is a rather extreme assumption, since it has been argued (Heger et al. 2003) that non-rotating masssive
  stars of about solar metallicity collapse into black holes, if their
  mass is in the 30-60 M$_{\odot}$ \ range (see their Fig. 1). Notice that stars in the 60-120 M$_{\odot}$ \ range (the most
  important $^{22}$Ne \  producers) end as  black
  holes in their scheme. However,  for slightly higher metallicities - such as those
  encountered in the inner Galactic disk - they find that  only neutron stars are formed,
  because higher stellar mass losses result in a less massive star at explosion. It should be noticed that
  the details of massive star explosions remain poorly understood at present (see e.g. Hanke et al. 2011) and so is the fate of a massive star above 30 M$_{\odot}$ \ (see Fryer et al. 2011 for a recent
  - but certainly not  definitive - assessment).
 The fate of the rotating mass losing stars considered here is even less well known.

  In view of the  aformentioned uncertainties, we feel that the scenario proposed here can be considered as valid
  at present, although future refinements in our understanding of massive star explosions may change it quantitatively
  (and even qualitatively, if it turns out that most masive stars above, say, 50 M$_{\odot}$, end up as black holes).

Finally, we explore the idea that CR are accelerated in superbubbles, in which case their composition results from the
ejecta of both stellar winds and SN explosions. We first notice that simple nucleosynthesis arguments suggest that
the resulting composition, averaged over the stellar IMF, should be very close to solar. We  demonstrate this quantitatively,
with a simple model for the evolving composition of a superbubble, enriched first by the ($^{22}$Ne \ rich) winds of the
most massive stars, then by the ($^{20}$Ne \ rich) SN ejecta of less massive stars. We find that, after a few Myr the superbubble ($^{22}$Ne/$^{20}$Ne)$_{SB}$  \ ratio
tends to solar, and so does the average ($^{22}$Ne/$^{20}$Ne)$_{CR}$ \ ratio in accelerated particles. We conclude that superbubbles cannot provide
the observed high  $R_{Obs}$ value of CR sources and, therefore, are not the main site of CR acceleration.
 On the contrary, SN remnants - including those expanding in the pre-explosion environment of a stellar wind - appear
 as suitable sites of GCR acceleration.

\begin{acknowledgements}

I am grateful to G. Meynet and M. Limongi for providing tables with yields and other properties
of their stellar models and to the referee, L. Drury, for  constructive remarks.

\end{acknowledgements}

redpajama_set_name: RedPajamaArXiv


Generally acceptable to great fit, but not reliable enough to depend on in unfamiliar applications. Quality is fair, but they definitely do not have the life span of most OE hoses. They are generally acceptable for older vehicle with fairly short life expectancy.
These pulleys are very hit and miss on quality, primarily based on bearing quality. The actual pulleys are of good quality regardless of whether a good quality or low quality bearing is installed. They do usually fit acceptably well, so that helps.
The closest thing to a perfect scan tool yet. Amazing capability and data refresh rate. Seamless and easy reprogramming. It does have some quirky PID names and advanced functions.
Full of surprises, both what it can do and what it cannot do. Lots of errors and incomplete tests. But, quite capable for the price. The hardware is of very low quality.


You fill out the Chambliss Judge Sign-Up form by 24 April to let me know which days of the meeting you're available to judge (the sign-up form takes about a minute to complete).
I schedule you to evaluate no more than about five undergraduate and/or graduate student posters, on no more than two days of the meeting (Monday–Wednesday), and only during the morning coffee breaks and/or evening poster sessions.
You judge the posters according to defined criteria. This will take a total of about an hour of your time at the meeting.
That's it! I can promise that our student poster authors will sincerely appreciate your expertise in judging their work, as will I.
For this upcoming competition, I am in particular need of full members to judge; however, if you're an advanced junior member with all coursework completed and you're not entering the competition yourself, you may also sign up to judge — we'll have you judging some of our undergraduate entrants.
Please complete the sign-up form by 24 April.
If you have any questions, please contact me via email. I hope to see you on my judging roster!


text:
Dit artikel geeft een chronologisch overzicht van alle belangrijke gebeurtenissen per dag in de maand mei van het jaar 2010.

Gebeurtenissen

1 mei
 Griekenland, dat in een diepe staatsschuldencrisis verkeert, bereikt overeenstemming met het Internationaal Monetair Fonds en de Europese Unie over maatregelen om een staatsbankroet af te wenden.
 In de Chinese stad Shanghai gaat de Expo 2010 van start.
 Op Times Square in New York wordt een bomaanslag verijdeld.
 Door aanslagen op de Abdalla Shideyemoskee in Mogadishu vallen 39 doden.

2 mei
 FC Twente wordt landskampioen voetbal in Nederland.

3 mei
 In Ukkel blijft de temperatuur steken op 7 graden. Dat is de laagste maximumtemperatuur die ooit op een 3 mei werd opgetekend, aldus weerman Eddy De Mey.
 Neil Robertson wint het WK snooker voor de eerste keer in zijn carrière. In de finale is hij met 18-13 te sterk voor Graeme Dott.

4 mei
 Een man verstoort door te gaan schreeuwen de Nederlandse dodenherdenking op de Dam in Amsterdam - in de ontstane paniek raken 63 mensen gewond.

5 mei
 Tijdens de 24-uursstaking in Athene tegen de miljardenbezuinigingen komen er bij rellen drie mensen om.

6 mei
 In het Verenigd Koninkrijk wordt de Conservatieve Partij van oppositieleider David Cameron de grootste partij bij de parlementsverkiezingen, maar ze haalt geen absolute meerderheid, waardoor voor het eerst sinds decennia een coalitie waarschijnlijk wordt.
 Bij de Mauritiaanse parlementsverkiezingen behaalt de centrumlinkse Alliance de l'Avenir, met onder meer de Parti Travailliste van huidig premier Navin Ramgoolam, een ruime overwinning op de Alliance du Cœur met de marxistische MMM van oppositieleider Paul Bérenger. Ramgoolam kan aanblijven als regeringsleider.
 In Nigeria volgt interim-president Goodluck Jonathan van de centristische PDP partijgenoot Umaru Yar'Adua op als president nadat deze overlijdt aan de gevolgen van pericarditis.

7 mei
 Chili wordt lid van de OESO.

8 mei
 De 93e Ronde van Italië gaat van start in Amsterdam, waar de proloog wordt gewonnen door de Brit Bradley Wiggins. De tweede en derde etappe finishen respectievelijk in Utrecht en Middelburg.
 Tientallen kompels en reddingwerkers komen om het leven bij een mijnramp in Siberië.

9 mei
 De EU-ministers van financiën en het IMF stellen gezamenlijk 720 miljard euro beschikbaar die gebruikt kan worden om landen in de problemen leningen te kunnen verstrekken, en zo een ineenstorting van de euro te voorkomen.
 In Noordrijn-Westfalen, de meest bevolkte deelstaat in Duitsland, verliest de christendemocratische CDU 10% bij de Landdagverkiezingen; het kabinet-Merkel II raakt daardoor tevens zijn meerderheid in de Bondsraad kwijt.

10 mei
 In de Filipijnen kiest men een opvolger voor president Gloria Macapagal-Arroyo, een nieuw Huis van Afgevaardigden en opvolgers voor 12 van de 24 senatoren. Lokaal zijn er tegelijk provinciale-, stadsbestuurs- en gemeenteraadsverkiezingen.
 Bij verschillende aanslagen in de Iraakse hoofdstad Bagdad worden in totaal meer dan honderd mensen gedood en vallen honderden gewonden.
 De Belgische schrijver Bernard Dewulf wint de Libris Literatuur Prijs met zijn novelle Kleine dagen.

11 mei
 De Britse socialistische premier Gordon Brown biedt zijn ontslag aan bij koningin Elizabeth II. David Cameron van de Conservatieve Partij neemt dezelfde dag zijn functie over en zal een kabinet leiden dat een coalitie is van zijn partij met de Liberal Democrats van Nick Clegg.

12 mei
 Afriqiyah Airways-vlucht 771, een internationale passagiersvlucht uit Zuid-Afrika met 104 personen aan boord, stort vlak voor de landing neer bij de Libische hoofdstad Tripoli.
 De Europese Commissie geeft groen licht voor de toetreding van Estland tot de eurozone, waarmee het zo goed als zeker is dat de Baltische staat op 1 januari 2011 als 17e EU-land de euro als munteenheid krijgt.
 Na HFC Haarlem is BV Veendam de tweede Nederlandse profvoetbalclub die in 2010 failliet wordt verklaard, maar dat wordt later teruggedraaid.
 Atlético Madrid verslaat Fulham FC en wordt zo de winnaar van de eerste editie van de UEFA Europa League.

13 mei
 Het half-afzinkbaar platform Aban Pearl zinkt in de Golf van Paria voor de kust van Venezuela. Het was ten tijde van het ongeval niet verbonden met een boorput.
 Aanhangers van oud-president Koermanbek Bakijev bestormen overheidsgebouwen op drie plekken in Zuid-Kirgizië.

14 mei
 Spaceshuttle Atlantis vertrekt voor zijn laatste missie naar het internationaal ruimtestation ISS.

15 mei
 Iran, Turkije en Brazilië sluiten een akkoord over het Iraanse nucleaire programma, waarbij Iran het merendeel van zijn voorraad laag-verrijkt uranium naar Turkije overbrengt.
 Een Antonov An-28 van Blue Wing Airlines met acht inzittenden stort neer neer in de jungle in het oosten van Suriname.
 De Australische Jessica Watson komt aan in haar thuisland nadat ze als jongste persoon ooit solo de wereld rond is gezeild.
 In Antwerpen opent het MAS, of voluit het "Museum aan de Stroom", de deuren. Het is een nieuw museum voor stad, haven en scheepvaart. Het herbergt de collecties van het Volkskundemuseum, het Nationaal Scheepvaartmuseum en deels van het Vleeshuis.

17 mei
 Ten noorden van de Afghaanse hoofdstad Kabul stort Pamir Airways-vlucht 112, met 43 personen aan boord, neer. Het toestel was onderweg van provinciehoofdstad Konduz naar Kabul.
 Bij een helikoptercrash in de Filipijnse stad Lucena komen zes mensen om het leven, onder wie Rafael Nantes, de gouverneur van de provincie Quezon.
 In delen van Europa worden wederom luchthavens gesloten in verband met de aswolk van de vulkaan onder de Eyjafjallajökull-gletsjer in IJsland.

18 mei
 Bij parlementsverkiezingen in de Dominicaanse Republiek behaalt de centrumlinkse PLD van president Leonel Fernández een overweldigende meerderheid. In de senaat gaat van de 32 zetels alleen die die de provincie La Altagracia vertegenwoordigt naar een andere partij, de rechtse PRSC.

19 mei
 In de Thaise hoofdstad Bangkok maakt het leger met geweld een einde aan de wekenlange bezetting van het zakendistrict door de Roodhemden.
 Een internationaal onderzoeksteam concludeert dat Noord-Korea verantwoordelijk is voor het zinken van het Zuid-Koreaanse marineschip Cheonan in maart 2010, waarbij 46 van de 104 bemanningsleden omkwamen.

20 mei
 Uit het Musée d'Art Moderne de la Ville de Paris worden vijf schilderijen gestolen met een gezamenlijke waarde van 100 miljoen euro.
 Een groep Amerikaanse wetenschappers onder leiding van Craig Venter maakt in het blad Science bekend erin geslaagd te zijn "synthetisch leven" te creëren.
 Een ontploffing in een mijn in het noorden van Turkije kost aan circa dertig mijnwerkers het leven.
 Teken-Mohammed-dag, een protestactie voor de vrijheid van meningsuiting, uitgeroepen na de censurering van aflevering 201 van de animatieserie South Park, leidt tot consternatie in de islamitische wereld.
 KV RS Waasland - SK Beveren ontstaat uit een fusie tussen Red Star Waasland en KSK Beveren.

21 mei
 Prinses Máxima reikt, in aanwezigheid van koningin Beatrix en prins Willem-Alexander, de Appeltjes van Oranje uit. De Appeltjes zijn jaarlijkse prijzen van het Oranje Fonds voor bijzondere, innovatieve of succesvolle projecten op sociaal gebied. De prijs bestaat uit een geldbedrag van € 15.000 en een oorkonde.
 Bagdad, de hoofdstad van Irak, wordt getroffen door een aanslag, die het leven kost aan zeker 27 mensen.

22 mei
 Een Boeing 737 van Air India Express stort na een mislukte landing neer nabij de luchthaven van Mangalore in de Indiase staat Karnataka. Het toestel kwam van Dubai in de Verenigde Arabische Emiraten. Zeven van de 166 inzittenden overleven de crash.
 Voetbalclub Inter Milaan wint de Champions League. Ze verslaat in de finale het FC Bayern München van Louis van Gaal met 2-0.

24 mei
 In Ethiopië behaalt het linkse EPRDF van premier Meles Zenawi in alle regio's een ruime overwinning bij fel gecontesteerde parlementsverkiezingen. De Medrekcoalitie van acht oppositiepartijen scoort onder de verwachtingen.
 Bij parlementsverkiezingen in de niet-erkende republiek Nagorno-Karabach (Azerbeidzjan) wordt het Vrij Vaderland van premier Arajik Haroetjoenian met 46% van de stemmen de grootste partij. De Democratische Partij en de Armeens-nationalistische Dashnak worden respectievelijk tweede en derde partij.

25 mei
 In de Jamaicaanse hoofdstad Kingston vallen meer dan 30 doden tijdens gevechten van leger en politie met aanhangers van een ingerekende drugsbaron.
 De Surinaamse parlementsverkiezingen worden gewonnen door de Mega Combinatie van de in Nederland voor drugshandel veroordeelde kandidaat-president Desi Bouterse, het voormalige moederland reageert bezorgd. De coalitie van huidig president Ronald Venetiaan gaat fel achteruit.

26 mei
 Kamla Persad-Bissessar volgt na de overwinning van haar UNC bij parlementsverkiezingen Patrick Manning van de conservatief-socialistische PNM op als premier van Trinidad en Tobago, en wordt zo de eerste vrouwelijke premier van het land.

28 mei
 Maoïstische rebellen laten in de deelstaat West-Bengalen in het noordoosten van India een trein ontsporen, er vallen 145 doden.

29 mei
 Koningin Beatrix en kroonprins Willem-Alexander zijn in de Nieuwe Kerk in Middelburg aanwezig bij de uitreiking van de Four Freedoms Awards. Het Amerikaanse Franklin and Eleanor Roosevelt Institute kent deze onderscheiding toe aan personen die op mondiaal niveau een substantiële, moedige en persoonlijke bijdrage hebben geleverd aan de vrijheid in de wereld. De prijzen worden toegekend aan onder meer het Europese Hof voor de Rechten van de Mens en dr. Asma Jahangir uit Pakistan.
 Bij de Tsjechische parlementsverkiezingen wint de sociaaldemocratische ČSSD nipt van de centrumrechtse ODS van president Václav Klaus.
 Na internationale kritiek verleent de Malawische president Bingu wa Mutharika gratie aan een homoseksueel koppel dat tot een jarenlange celstraf was veroordeeld.
 Lena Meyer-Landrut wint voor Duitsland met het nummer Satellite het Eurovisiesongfestival 2010 met 246 punten. De Belgische vertegenwoordiger Tom Dice eindigt als zesde met "Me and my guitar".

30 mei
 In de Burundese hoofdstad Bujumbura breken rellen uit na het bekendmaken van de resultaten van de gemeenteraadsverkiezingen van vorige week, waarbij de CNDD-FDD van president Pierre Nkurunziza 64% van de stemmen zou hebben gehaald.
 De Rus Denis Kozjoechin wint de Koningin Elisabethwedstrijd voor piano 2010.

31 mei
 Horst Köhler treedt af als bondspresident van Duitsland nadat ophef is ontstaan rond zijn uitspraken over de Duitse militaire inzet in Afghanistan.
 De Israëlische marine valt een konvooi boten met hulpgoederen en activisten van de Free Gaza Movement aan dat op weg was naar de Gazastrook. Er vallen minimaal 19 doden. Zie: Onderschepping van het scheepskonvooi voor Gaza
 Ex-minister van Defensie Juan Manuel Santos van de centrumrechtse PSUN van huidig president Álvaro Uribe wint de eerste ronde van de Colombiaanse presidentsverkiezingen. In een tweede ronde zal hij het opnemen tegen de groene Antanas Mockus, voormalig burgemeester van hoofdstad Bogota.
 Bij een raketaanval van Koerdische rebellen op een legerbasis nabij de stad İskenderun in het zuiden van Turkije komen zes Turkse militairen om het leven.

Overleden

 05
Mei

redpajama_set_name: RedPajamaWikipedia


text:
The art of the Eastern Roman Empire and of its capital Byzantium, or Constantinople. The term denotes more especially those qualities which distinguish this art from that of other countries, or which have caused it to exert an influence upon the art of regions outside of the Eastern Empire. Christian art was dependent for the representation of its new conceptions upon the forms which the time and place of its origin happened to offer. In the beginning, whether at Rome, Ravenna, or Byzantium (Constantinople), it was equally influenced by classical art and by Eastern inclination to allegory. It is a distinguishing characteristic of Constantinople, however, that it was able to maintain a more uniform classical tradition in the face of manifold Oriental influences. These two elements, from the time of Constantine, developed in the Byzantine art more and more of an individual character, though account must also be taken of the friendly intercourse with Western Europe during several hundred years. Beginning with the seventh century, the contrast between the art of the Eastern Empire and that of the Western grew more marked, and Byzantine art underwent a change. It rose to great splendour under the Macedonian emperors (867-1056), then declined up to 1453, and has since existed in the East in a petrified form, so to speak, up to the present time.
The Byzantine Question
In regard to the first period of Byzantine art, which closed either before the reign of Justinian or at the end of the sixth century, scholars differ greatly. Some date Byzantine art proper from the time of Constantine's establishment of his capital. They base this opinion upon certain differences between the art remains of the first period of the Eastern Roman Empire and those of the Western Roman Empire, which differences they maintain are essential. Other scholars hold these peculiarities to be unessential, since they find them here and there in Western countries as well, a fact which the former critics ascribe to Oriental influence. Still other scholars disagree with both views, and distinguish between Oriental art and that specifically Byzantine; that is, between the art of Byzantium, or Constantinople, and that of her dependent provinces, Asia Minor, Syria, Persia, and Egypt. This is a fairly good solution of the "Byzantine question". But as it is difficult to distinguish in detail the combinations of old classic and Christian with Oriental art, we can only group together the principal characteristics of the new style and its materials, with a few examples.
The introduction of Eastern court ceremonial by Constantine was accompanied in the domain of art by the appearances of extraordinary gorgeousness and pomp, expressed, however, with stiffness and formality. The power and pride of the new empire offered the means for great undertakings and gave the impulse to them. The Proconnesian marble, found in the vicinity of the capital, and the stone obtained from other rich quarries provided the material, and, long before this era, the art of working in stone had reached a high state of development, especially in Asia Minor. Moreover, the East had been from ancient times the home of the minor arts. In Constantinople there flourished, along with the art of decorative sculpture, the arts of stone-carving, of working in metal and ivory, of ornamental bronze work, of enamelling, of weaving, and the art of miniature-painting. From classical and ancient Christian art Byzantine genius derived a correct combination of the ideal with truth to nature, harmonious unity along with precision in details, as well as the fondness for mosaics, frescoes, and pictures on panels, in opposition to the dislike of non-Christian and sectarian Orientals to pictorial representation. The iconoclasm of the eighth and ninth centuries wrought great destruction in the domain of art, but these outbreaks were successfully suppressed.
In regard to the influence of the Byzantine style on architecture see BYZANTINE ARCHITECTURE. As to the other arts a few examples may here be given. The church of St. Sophia was adorned in the sixth century with a splendour worthy of Solomon. The interior was sumptuously decorated with mosaics upon a golden background. These mosaics, it is true, with the exception of an "Adoration of Christ by the Emperor", were destroyed, but they were replaced later by others. Some of the walls were ornamented with designs of grape-vines with golden leaves. Pictures of animals decorated the walls of the portico. A silver choir-screen rose above pillars, in the capitals of which medallions of Christ, the Blessed Virgin, saints, and prophets were carved. This is the so-called iconostasis. The altar was of gold inlaid with precious stones; the altar-cloth was of brocaded silk in which were woven pictures of Christ, the prophets and the apostles. The ambo, according to description, was brilliant with gold, silver, precious stones, and ivory. At Parenzo, in Istria, and at Bauit, in Egypt, superb mosaic pictorial ornamentation dating from the sixth century is still preserved. A gold cross decorated with pictures in hammered work was presented by Justin II to the church of St. Peter and is still preserved at the Vatican. A number of ivory book- covers are also still in existence. The illuminated manuscripts of Rossano and Sinope date from the sixth century.
As regards the influence exerted by Byzantine art in the sixth century there can be no doubt that the architecture of Ravenna, though affected by other Eastern influences, strongly reminds us, in its splendid mosaics, of Constantinople. The Proconnesian capitals and other products of decorative art spread even more easily. Like Ravenna, Southern Italy and Gaul came under the influence of the East and Constantinople. Even more specifically Byzantine is African art. In Rome the races of Byzantine art are more difficult to discover than other Oriental influences. In the East itself pictorial art met with opposition, and decorative art came to the forefront. In general, however, after the rise of the Macedonian dynasty the Byzantine style gained the supremacy in all branches of art as well as in architecture. The Byzantine style spread in the East as well as in Northern Italy and Sicily. The numerous mosaic pictures, which are to be found everywhere, still strove to imitate classical models; their symbolism reminds us of the general symbolic tendency of early Christianity, and their form gradually becomes more stiff and fixed. (Painter's Book of Mount Athos.) Purely Oriental, however, was the dislike constantly increasing for sculpture in the round, and the preference for the flat ornamentation in architecture. To the same Oriental influence may be attributed the taste for costly and many-coloured stones and woven fabrics, for goldsmith-work, and enamel. For example, in the treasury of San Marco may be seen Byzantine reliquaries, ivory triptychs, chalices, costly fabrics, and specimens of pictorial art. Some are large and some small, but taken altogether they show how a church of the eleventh century was transformed into a veritable treasure-house. The same taste and the same characteristics of the art of Byzantium (Constantinople) have ever since maintained their supremacy in the East. For further bibliography see BYZANTINE ARCHITECTURE.
KONDAKOFF, Hist. de l'art byzantin considéré principalement dans les miniatures (tr. Paris, 1886-91); IDEM,Les émaux byzantins (Paris, 1892); MOLINIER, Hist. des arts appliqués à l'industrie, I, s.vv. Les ivoires; L'orfèvrerie; TEXIER-PULLAN, Byzantine Architecture (London, 1843-64); LETHABY-SWAINSON, The Church of St. Sophia (London, New York, 1894); FROTHINGHAM, Byzantine Artists in Italy in Am. Journal of Archæology (1894-95); MACPHERSON, The Church of the Holy Sepulchre in Eng. Hist. Rev. (1892), VII; WESTWOOD,Fictile Ivories in the South Kensington Museum (London, 1886); SCHULTZ-BARNSLEY, Phocis and the Dependent Monastery of St. Nicolas-in-the-fields near Skripou (London, 1901); VELUDO, La pala d'oro di S. Marco, Eng. and Fr. trs. (Venice, 1887); CAHIER AND MARTIN, Mélanges d'archéologie (Paris, 1847-56).
G. GIETMANN

redpajama_set_name: RedPajamaCommonCrawl


Portimão is located on the southern coast of Portugal. It was formerly known as Vila Nova de Portimão. The town has a population of about 56 000. Portimão is known for its sporting events. Portimão is well known in the sailing world for being an ideal location. It is famous for its fresh fish restaurants.


OR3D offer CT scanning services. CT is the collection of 2D X-Ray slices that are combined to provide a complete 3D dataset. The CT scanner allows the collection of extremely accurate and high resolution (0.25 micron) scan data.
OR3D offer CT scanning services. CT (Computed Tomography) is the collection of 2D X-Ray slices that are combined to provide a complete 3D dataset. The CT scanner allows the collection of extremely accurate and high resolution (0.25 micron) scan data. This enables the inspection and reverse engineering of very small and highly detailed components or components where surface scanning technologies cannot capture geometry.
This technique uses X-Ray technology enabling NDT (Non Destructive Testing) of part assembles or multi material components. Additional analysis such as porosity checks within castings or plastic injected mouldings offer additional insight.
If you have any questions or think that CT scanning may solve a problem you have, then please contact us.


text:
Home Opinion MLA Report: Prince Albert Northcote
MLA Report: Prince Albert Northcote
Alana Ross
Alana Ross, MLA
As we head into the new year, I send Season's Greetings to all our Greek Orthodox citizens and wish you peace and all the best in 2023.
Saskatchewan's economy continues to lead the nation. Despite global economic uncertainty, private companies continue to invest their dollars in Saskatchewan industries. Significant investments in energy, forestry, technology and agriculture will create jobs and opportunities in Saskatchewan throughout the new year and into the future.
BHP is accelerating the timeline on its $12 billion Jansen potash mine, the largest single investment in Saskatchewan history. This investment brings total capital investment in the province's potash industry to $30 billion since 2007. The mine is expected to be operational in 2027 and will provide over 600 direct jobs at the mine site and at the corporate office in Saskatoon.
Saskatchewan's mining industry also has the unique opportunity to position itself as a leader in lithium production. Global lithium demand is expected to grow significantly in the next decade as it is used in portable electronic devices and electric vehicles. Prairie Lithium and Grounded Lithium have both invested in drilling for lithium in the province. As drilling continues in the province, Saskatchewan can expect growth in the technology, mining, and energy sectors as well as job creation.
Investments by several agricultural companies have helped Saskatchewan maintain and grow its agricultural advantage. Federated Cooperatives Limited is moving forward with its Integrated Agriculture Complex, including a renewable diesel facility and canola crushing plant in partnership with AGT Foods. The renewable diesel plant, which will be the largest of its kind in Canada, represents nearly $2 billion in capital investment, plus 2,500 construction jobs and 150 operating jobs. The entire project will have estimated economic benefits of $4.5 billion.
As part of the Growth Plan goal to double the growth of the province's forestry sector, our government approved timber allocations to support four major forestry projects. Pivot Furniture, a Saskatchewan-based company in Prince Albert, was awarded a timber allocation of 101,000 cubic metres. Between 2021 and 2024 the forestry industry is expected to attract $1 billion in capital investment and create 3,100 jobs. The sector already supports over 8,000 jobs in the province.
The confidence of private investors demonstrates the momentum in our economy as we look to the year ahead. Economic forecasters predict resilience and prosperity for the province in 2023. The future is bright and it's a great time to live, work and raise your family in Saskatchewan.
Saskatchewan families have additional COVID-19 vaccination options for children and teens. As of January 9, children five to eleven are eligible for the Pfizer bivalent vaccine as a booster dose four months following their primary series or four months following their original COVID-19 vaccine booster dose. For children 12 to 17, the Novavax Nuvaxovid is available as a primary series and booster dose. Information on all COVID-19 vaccines approved by Health Canada is available at Saskatchewan.ca/covid19.
If you require assistance with provincial government programs or services, I encourage you to contact my constituency office at #7, 598 15th Street East. We can be reached at 306-763-7677 or by email at panorthcote.mla@sasktel.net.
rosscolumn
Previous articlePatients have right to refuse care says SHA
Next articleSaskatchewan Serious Incident Response Team confirms weapon was recovered during investigation into fatal officer-involved shooting
https://paherald.sk.ca/latest-posts-alana-ross/
The paradoxes of Canada's conservation policies
Millennials Opting Out of Having Children
'All in' cost become too high on WHL trade front
Aging well at home
Announcing corporate bonuses requires tact and empathy
Working outside this winter
Pandering support for Sask Party policies 'a clear and present danger' for the north
A Christmas wish list to fix our health care crisis
Troy Media - December 27, 2022
Canadians are concerned about the mental health of health workers —...
Quoi Media - May 9, 2022
Lyle Karasiuk - January 4, 2023


text: Q: What information does Android's MediaStore use to determine mime type? I'm using the MediaStore to load a list of files on the user's device and using MediaStore.Files.FileColumns.MIME_TYPE in my projection.  I have a file that is a 3gp video which has a file name ending in .3gp, and the data I'm getting back is that its mime type is video/mp4.

redpajama_set_name: RedPajamaStackExchange