Activity and quiescence in galaxies at redshifts 1.4<z<3.5. The role of the starburst temperature

Abstract: We investigate 'activity' and 'quiescence' in galaxies at relatively high redshifts by modelling the line(and continuum) spectra of each object. The models account consistently for photoionization and shocks. We claim that the starburst effective temperature, the flux from an AGN, and the shock velocity are critical to activity. The results confirm that two sample galaxies show intense starburst activity with temperatures reaching Ts=2x10^5K and shock velocities Vs> 250 km/s, while for the remaining galaxies of our sample the models show quiescent star formation with Ts< 7x10^4K. A Seyfert 2 like AGN is proposed in one galaxy. The O/H relative abundances derived by the detailed modelling of the spectra are nearly solar for all the sample galaxies, in contrast to those obtained by direct methods.

\section{Introduction}

In  previous papers (Contini 2014, hereafter Paper I, Contini 2013, and references therein) we have
investigated the physical conditions and element abundances in galaxies   throughout  redshifts
0.001$<$z$<$3.4   by   detailed modelling of the emission line ratios.
We have found  that both active galactic nuclei (AGN) and starbursts (SB) coexist in  almost all galaxies and that
collisional process   should be accounted for.

In this paper we present  modelling  results  of
 new spectra observed from selected  galaxies at 1.4$<$z$<$3.5 (Table 1), in particular those
with enough lines (e.g. of  oxygen in different ionization levels)  that  constrain the models  (Paper I).
We cross-check the models, when possible, by the upper limits of significant lines,
by reproduction of the observed spectral energy distribution (SED) of the continuum,
by the morphological evidence, etc.

\begin{table}
\caption{Description of the observed galaxies}
\tiny{
\begin{tabular}{lccccccccccc} \hline  \hline
\ object    & z     & galaxy type   & telescope   &   ref \\  \hline
\ CDFS-695  & 2.225 & luminous  red  & Gemini NIRS&1 \\
\ 1255-0    & 2.1865&quiescent & Gemini NIRS&2   \\
\ A1835     & 2.073 &low-mass SF&Hale tripelspec&3 \\
\ A773      & 2.303 &low-mass SF&Hale tripelspec&3 \\
\ MACS0717  &2.55  &low-mass SF&Hale tripelspec&3 \\
\ WISP     &1.444-2.315&line selected      & Mag. Baade (FIRE)&4\\
\ J1000+0221S&3.417 & SF dwarf &CANDLES team& 5 \\ \hline
\end{tabular}}

1: van Dokkum et al (2005);2:  Kriek et al (2006);3: Belli et al (2013);4: Masters et al (2014);5: Amor\'{i}n et al. (2014)
\end{table}

\begin{figure}
\includegraphics[width=10.2cm]{fig1.eps}
\caption{
Black open squares: the observation data. Black: AGN(V=100 (dash-dot), 200 (dot), 300 (solid), 500 km s$^{-1}$ (dash);
x: $F$=0, log$F$=9,10,11,11.5,12. SB:red (T=10$^4$K), blue (T=5 10$^4$K), magenta (T=10$^5$K),
x:V=100,300,500$\rm km\, s^{-1}$, $U$=0.01 (dot), 0.1 (solid), 1 (dash), 10 (dot-dash)
}
\end{figure}

We associate the `activity' in galaxies (e.g. CDFS-695,  van Dokkum et al 2005) to
 an AGN, a SB, and to shocks on the basis of  a strong wind throughout the galaxy, while we refer to
`quiescent' galaxies e.g. 1255-0,   which  was classified  as a quiescent ultra-dense object by Kriek et al. (2009)
because of emitting  weak lines. Those imply a low luminosity AGN and
 a low star formation rate (SFR = 1-3 M$_{\odot}$ yr$^{-1}$)  compared with models of  stellar population
synthesis SEDs.
In this paper we investigate activity and quiescence  in
 the galaxies described in Table 1,
considering  the  SB effective temperature ($\rm T_{*}$~) and ionization parameter ($U$),
 the AGN flux intensity ($F$), and  the shock velocity ($\rm V_{s}$~) as critical parameters.

 In Sect. 2 we   present  calculations and modelling of the galaxy spectra.
Results are presented and discussed in Sect. 3.  Concluding remarks follow in Sect. 4.

\begin{table}
\caption{Modelling  the  line ratios to H${\beta}$~=1 observed  by  van Dokkum et al. (2005) and Kriek et al. (2009)}
\tiny{
\begin{tabular}{lccccccccccc} \hline  \hline
\ line &            Obs$^1$  &     M$_{AGN}$ &      M$_{SB}$             & Obs$^2$             & m$_{AGN}$          & m$_{SB}$     \\ \hline
\ z    &            2.225    &   -           &     -       &         2.1865                    &     -              &  -   \\
\ Ly$\alpha$~         &        -           &       41.87       &       53.97           &    -                &  30.39             &27.7       \\
\ NV 1240     &          -           &     0.007      &     0.072              &    -                &  0.31              &1.27            \\
\ CIV 1550    &         -           &      0.53        &      0.67              &    -                &  0.8               & 2.64           \\
\ HeII 1640   &         -           &      3.43        &      3.32              &    -                &  0.89              &0.16             \\
\ [OII] 3727+ &      6.13            &    6.7          &    6.4                 & $>$5.33             & 6.                 &5.             \\
\ [NeIII] 3869+&     2.17             &    1.76        &    2.5                 &      -              &       -            &   -         \\
\ [OIII] 5007+&      10.2           &       10          &   11.                 & $<$1.36             &1.48                &1.4          \\
\ [OI] 6300+  &     $<$0.26          &      0.7        &    0.9                 &    -                &    -               &      -      \\
\ [NII] 6548+  &     2.46            &  3.2           &  3.6                    & $>$5.26                &4.7                 & 5.          \\
\  H${\alpha}$~         & 3.                  &  3.            & 3.                     & 3.07                 & 3.                & 3.           \\
\ [SII] 6718  &      0.5             &    0.5       &     0.29                 & $>$1.8             & 1.1                & 1.1            \\
\ [SII] 6731  &      0.45            &    0.52      &     0.37                 & $>$0.67            &1.1                 & 1.7           \\
\ H${\beta}$~$_{obs}^3$ &  2.9 10$^{-17}$     &    -        &         -                & $<$4.2 10$^{-18}$  &     -              &      -            \\
\ H${\beta}$~$_{mod}^3$   &      -            &   0.059     &    0.051                 &    -               &0.013               & 0.009             \\
\ $\rm V_{s}$~ ($\rm km\, s^{-1}$)    &        -           &    100       &     250                  &   -                 &180                 &150             \\
\ $\rm n_{0}$ ($\rm cm^{-3}$)    &        -           &    150       &     130                  &   -                 &100                 & 300            \\
\ $F$ $^3$      &        -           &    40        &     -                    &-                    & 1.6               &  -           \\
\ $\rm T_{*}$~ (10$^4$ K) &       -           &    -        &     20                    &-                    &-                   & 7             \\
\ $U$            &       -          &     -        &      0.03                 &-                    &-                 & 0.007     \\
\ $D$ (10$^{19}$) cm&     -          &     2         &      1                   &-                    &1.                & 0.6    \\
\ N/H (10$^{-4}$)     &  -          &     0.8       &      0.4                  &-                   & 1.0               & 1.2  \\
\ O/H (10$^{-4}$)     &  -          &     6.6       &      5.6                  &-                   &5.6                 &6.6   \\
\  S/H (10$^{-4}$)&     -           &     0.06      &      0.06                 &-                   &0.07               &0.06    \\ \hline
\end{tabular}}

$^1$ van Dokkum et al. (2005);
$^2$ Kriek et al. (2009);
$^3$ in $\rm erg\, cm^{-2}\, s^{-1}$ ;
$^4$ in  10$^9$ photons cm$^{-2}$ s$^{-1}$ eV$^{-1}$ at the Lyman limit

\end{table}

\section{The models}

For the calculations of  both the  line and continuum spectra we adopted the code
  {\sc suma}\footnote{http://wise-obs.tau.ac.il/$\sim$marcel/suma/index.htm} (see also Contini  et al. 2012), which
accounts for the coupled effect of
photoionization from a radiation source and shocks.
The input parameters related with shocks are  the  shock velocity $\rm V_{s}$~, the atomic preshock density $\rm n_{0}$, and
the preshock magnetic field $\rm B_{0}$.
 $\rm V_{s}$~ determines the maximum
temperature T of the gas downstream, near the shock front. The cooling rate in the recombination zone  depends
on the square density of the gas and on T in  single slabs of gas.
We adopt  $\rm B_{0}$=10$^{-4}$ Gauss.

The input parameter  that represents the radiation field  in AGNs is the power-law
flux  from the active centre $F$  in number of photons cm$^{-2}$ s$^{-1}$ eV$^{-1}$ at the Lyman limit.
The spectral indices are $\alpha_{UV}$=-1.5
and $\alpha_X$=-0.7.
For SBs the radiation is a   black body (bb).
The input parameters are the  effective temperature of the  starburst $\rm T_{*}$~
and  the ionization parameter $U$
(in  number of photons per  number  of electrons at the nebula).

Additional input parameters are the geometrical thickness of the emitting nebula $D$ which
 determines  whether the model is radiation-  or matter-bound, the
  abundances of  He, C, N, O, Ne, Mg, Si, S, A, and Fe relative to H and  the dust-to-gas
($d/g$) ratio which  affect  the cooling rate.
The input parameters are all correlated throughout the Rankine-Hugoniot equations which regulate the
shock dynamics at the shock-front and downstream, and  in the ionization equations  that are used to calculate the
ion  fraction abundances for each element in each slab of gas downstream.
Uncertainty in the calculation results is in average $<$10\%.

The primary radiation flux
and the secondary (diffuse) radiation flux,
emitted by the gas slabs heated at relatively high temperatures by the shock,
are calculated by radiative transfer throughout the clouds.
The  line ratios  calculated by a power-law flux  or by   bb   are different, so
we calculate  AGN and SB dominated models separately for all the galaxies.
There is no  degeneracy   for the results of  [OIII]/H${\beta}$~ and  [OII]/H${\beta}$~  calculated  by SB models
as functions of   $\rm T_{*}$~  or of $U$, because the line intensities, which
depend on the atomic configuration of the  corresponding ions,
 show different trends with $\rm T_{*}$~ and $U$.

Briefly,
to reproduce the data by  model calculations, we  first choose the input parameter range from the grids
of models (Fig. 1) calculated
by Contini \& Viegas (2001a,b, hereafter GRID1 and  GRID2, for AGN and SB, respectively).
Then we refine the fit (Tables 2-6)
calculating  new  grids  suitable to each spectrum.
We cross-check the models resulting from  the line spectrum fit  by reproducing consistently
 the  data   in the continuum SED diagrams (Fig. 2), where
 for each model    one line  represents the
bremsstrahlung (from the gas) and  another, in the IR,  the reprocessed
radiation  from dust. An initial grain radius of 0.2 $\mu$m~ is adopted.
The bb radiation  fluxes from the SB  and   from  the background old population stars
also appear in  Fig. 2 diagrams.

The results of detailed modelling  are different from those obtained by direct methods
because we calculate  the emitted line fluxes  integrating throughout the  downstream recombination region
where the  shock determines the profile of  electron temperature and density.

\begin{table*}
\caption{Modelling the  line ratios to H${\beta}$~=1 observed   by  Belli et al. (2013)}
\tiny{
\begin{tabular}{lccccccccc} \hline  \hline
\ line       & obs (A1835) & MB1$_{SB}$& MB1$_{AGN}$& obs(A773)$^1$ & MB2$_{SB}$& MB2$_{AGN}$& obs(MACS0717)&MB3$_{SB}$& MB3$_{AGN}$\\ \hline
\ z          &2.073        &     -     &     -      & 2.303         &     -     &       -    & 2.55         &      -   &     -      \\
\ [OII]3727+   & 2.17       & 2.3       & 2.2       & 4.65          & 4.6       & 4.8        & 2.27         & 2.2      & 2.3 \\
\ [NeIII]3870& 0.35        & 0.38      & 0.6        & $<$4.3        & 1         & 0.85       & $<$ 1.07     & 0.29     & 0.5 \\
\ [[OIII]5007+& $<$3.5     & 3.45      & 3.47       & 6.5           & 6.4       & 6.7        & 2.67         & 2.8      & 2.6 \\
\ H${\alpha}$~        & 2.73        & 2.96      & 3.25       & 3.            & 3.6       & 3.4        & 2.9          & 2.96     & 3.16 \\
\  [NII] 6384& 0.33        &0.386       &0.39       & 0.59          & 0.47      & 0.6        & 0.7          & 0.57     & 0.59 \\
\ H${\beta}$~$_{obs}^2$ &16.6$\pm$3.1 10$^{-17}$&-          & -         &8.3$\pm$3.5 10$^{-17}$   &  -        & -       & 8.5$\pm$2.9 10$^{-17}$ &-      & -    \\
\ H${\beta}$~$_{mod}^2$ & -           &0.067      & 0.033      & -             &0.85       & 0.035      & -            & 0.017    & 0.013 \\
\ $\rm V_{s}$~ ($\rm km\, s^{-1}$) & -           &130        & 100        & -             & 100       & 100        &  -           & 140      & 100  \\
\ $\rm n_{0}$ ($\rm cm^{-3}$) & -           &100        & 200        & -             & 200       & 140        & -            & 100      & 80  \\
\ $\rm T_{*}$~ (10$^4$ K) & -       & 7         & -          & -             & 20        & -          & -            & 6.3      & -   \\
\ $U$           &  -       & 0.02      & -          & -             & 0.8       & -          & -            & 0.02     & -   \\
\ $F$ $^3$      & -        &  -        & 5.5        & -             & -         & 10          & -            & -        & 2 \\
\ $D$ (10$^{19}$cm)  &   -  &   1      & 1          & -             & 6         & 1          &  -           & 1        & 1.5 \\
\ N/H (10$^{-4}$)     &   -  & 0.3      & 0.3        & -             & 0.3       & 0.3        &  -           & 0.3      & 0.3 \\
\ O/H (10$^{-4}$)     &   -  & 6.6      & 5.5        & -             & 7.6       & 6.6        &  -           & 6.6      & 5.6 \\ \hline

\end{tabular}}

$^1$ reddening corrected;
$^2$  in   $\rm erg\, cm^{-2}\, s^{-1}$ ;
$^3$  in 10$^{9}$ photons cm$^{-2}$ s$^{-1}$ eV$^{-1}$ at the Lyman limit

\end{table*}

\section{Results and discussion}

\subsection{The physical conditions and the element abundances}

The set of  input parameters   which  yields the  best  fit  of calculated to observed data
 determines the  model for each galaxy.
It represents averaged conditions  for each object  when the   observations  cover
the whole galaxy.

Generally,
we have found that  $\rm V_{s}$~, $\rm n_{0}$, and $D$  result within the ranges calculated for a large number
of objects (Paper I, Fig. 5), while  the
 SB temperatures (Tables 2-3-5-6)  $\rm T_{*}$~=2.5 10$^5$K  for
CDFS-695 (van Dokkum et al)  and A773 (Belli et al) are outstanding.
$\rm T_{*}$~ $\sim$7 10$^4$ K is found  in the quiescent galaxy 1255-0 (Kriek et al),  in
 WISP 170-106 (Masters et al.), and in A1835 (Belli et al.) and
 $\rm T_{*}$~=6.3 10$^4$K  in MCS 0717 (Belli et al.).

The  results  suggest that a  strong starburst activity  occurred   in
the two  sample galaxies  with maximum $\rm T_{*}$~,  while quiescence is more adapted to
the other galaxies with $\rm T_{*}$~ $\leq$7 10$^4$ K  similar to  the starburst temperatures  calculated for most
of the objects in Paper I and in Contini (2013).
 $\rm T_{*}$~=2 10$^5$ K  accompanied by a relatively high $\rm V_{s}$~ in  SBs  is  suitable to young stars
and it may indicate a revival of activity in these galaxies.
Such high temperatures were found, for example, in  the  recurrent nova T Pyxidis
 (Contini \& Prialnik 1997), in the  symbiotic star  HM Sge at outburst (Angeloni et al. 2007), etc.
The two  `active' galaxies  also  show a relatively high $U$ ($\sim$ 1)
and  correspond to SFR$\sim$60 and 400 M$_{\odot}$ yr$^{-1}$, respectively, which  are derived from the H${\alpha}$~ luminosities
(Kennicutt 1998).

 The modelling presented in Tables 2, 3 and 6 suggests that   AGN and   SB  coexist (see also Fig. 1)
because the line ratios in each spectrum are satisfactorily fitted (within the observed errors) by
both AGN and SB models.
The same  was found  throughout many local merger galaxies, due to
mutual  feedbacks.
The   AGN fluxes  calculated for the present  sample are similar to those found in
low luminosity AGN (LLAGN), except for  J1000+0221S (Amor\'{i}n et al. 2014)
which   suits a Seyfert 2 galaxy.

We have found  O/H  close to solar  (6.6 10$^{-4}$, Allen 1976) for all the objects (Fig. 3).
van Dokkum et al refer to diagnostic diagrams calculated by solar O/H and Kriek et al.
mention  solar or supersolar  O/H.
N/H modelled for Kriek et al. quiescent galaxy is similar to  the  maximum  (N/H=1.4 10$^{-4}$)
calculated in local galaxies (e.g. NGC 7212, Contini et al. 2012, fig. 14).
 Masters et al  suggest an average O/H=2.18 10$^{-4}$ for their sample  objects, while our models
 show O/H  slightly higher than solar (Allen 1976), but within  the  solar
O/H= 8.5 10$^{-4}$ given by Anders \& Grevesse (1989)(cf Paper I, table 1).
Calculated N/H are generally lower than solar ($\sim$ 10$^{-4}$), in agreement with the results
obtained for galaxies in the same z range.
For J1000+0.221S  we  found solar O/H and N/H, while
Amor\'{i}n et al. derived O/H=2.75$^{2.88}_{2.6}$ 10$^{-5}$
by the strong-line and calibration method presented by Maiolino et al. (2008) (see Paper I).

\begin{figure}
\centering
\includegraphics[width=4.4cm]{fig2a.eps}
\includegraphics[width=4.4cm]{fig2b.eps}
\caption{SED of Kriek et al (2009, fig. 2) and Belli et al (2009, fig. 1) MACS0717 galaxies.
Black circles : the data;
blue solid lines : AGN models;
red dashed lines : SB models;
magenta dot-dashed lines : bb radiation calculated for $\rm T_{*}$~ ;  black dotted : bb for $\rm T_{*}$~=10$^4$K}
\end{figure}

\begin{figure}
\centering
\includegraphics[width=6.0cm]{fig3.eps}
\caption{Direct method (d) versus model (m) results for O/H}
\end{figure}

\subsection{Notes on individual galaxies}

{\bf CDFS-695 }. We  cross-checked the results obtained by   modelling the line spectrum (Table 2, col.3 and 4)
calculating the distance r of the emitting cloud from the AGN in the NLR.
The H${\beta}$~ flux observed at Earth (2.9 10$^{-17}$ $\rm erg\, cm^{-2}\, s^{-1}$)  after reddening correction (Osterbrock 1974, eq 7.6, fig. 7.1)
results H${\beta}$~ $_{obs}$= 6.10$^{-16}$ $\rm erg\, cm^{-2}\, s^{-1}$. Combining it with the H${\beta}$~ flux calculated at the cloud (H${\beta}$~$_{mod}$) :
\noindent
 H${\beta}$~$_{obs}$ d$^2$ = H${\beta}$~$_{mod}$ r$^2$   (where d is the distance to Earth)
 and adopting a filling factor $\sim$ 1 throughout the NLR,
 r results $\sim$ 1.8 kpc, in rough  agreement with van Dokkum et al who suggested the presence of  shocks
on the basis of
the line ratios extending to  $\sim$ 10 kpc from the nucleus.
The distance of the emitting clouds from the SB region,  calculated in the same way,
results r$<$ 1.9  kpc,  adopting a filling factor $<$ 1.
The SB region does not  cover the entire galaxy.
SB have typical radii of 100-1000 pc (1-10\% of the size of their 'host' galaxies).

{\bf 1255-0}.
We constrained the  models (Table 2)
of this  `quiescent' galaxy  (Kriek et al 2009)
by the  relatively low [OIII]5007+/[OII]3727+,
 [OII]/H${\beta}$~, and [OIII]+/H${\beta}$~ (Table 2, col. 5).
The  modelling  results appear in Table 2, col. 6 and 7 (m$_{AGN}$, m$_{SB}$).
The relatively low calculated  H${\alpha}$~ flux ($\sim$  3 H${\beta}$~$_{mod}$) partly  explains the low observed H${\alpha}$~
flux (1.29$^{+0.16}_{-0.18}$  $\times$ 10$^{-17}$ $\rm erg\, cm^{-2}\, s^{-1}$) and
luminosity L$_{H\alpha}$~(0.46$^{+0.06}_{-0.06}$ $\times$ 10$^{-42}$ erg s$^{-1}$)  confirming
 the low SFR   observed by Kriek et al.

In Fig. 2 (left diagram)   we  model the  continuum SED   presented
by Kriek et al. (2009, fig. 2). The data are normalized to
unity at 1 $\mu$m~.
 The  bremsstrahlung calculated
by model m$_{AGN}$ reproduces the data  for $\nu$ between 10$^{14}$ and 10$^{15}$ Hz.
The  Planck function corresponding to 7 10$^4$ K
 peaks in the UV where the flux is  strongly absorbed by the Galaxy.
A few datapoints in the SED between 10$^{15}$ and 10$^{16}$ Hz fit
the   bremsstrahlung calculated by m$_{SB}$ with perhaps some contribution from
 the bb flux from the SB.
Models m$_{AGN}$  and m$_{SB}$ were calculated by    dust-to-gas ratios (by number)
of 10$^{-15}$ and  10$^{-14.3}$, respectively.

{\bf A1835, A773, and MACS0717}. The results for the galaxies
 from the Belli et al. sample are presented in Table 3.
Although the observed line ratios are typical of SBs, we have  searched for
AGNs in each object.
The results   show that  AGN could be present (Fig. 1).
 $\rm T_{*}$~  is relatively high, in particular for A773 (2 10$^5$ K).
The emitting clouds are  geometrically thick ($D$$\geq$ 10$^{19}$ cm)  which  seems unusual for AGNs
in  this z range ($D$$\leq$ 10$^{16}$ cm, Contini 2013, fig. 3).
The models selected for MACS0717 (Table 3)  are cross-checked by fitting  the  SED in Fig. 2.
The observed photometry
follows the bremsstrahlung maximum. The peak frequency depends on $\rm V_{s}$~ and the peak intensity on $F$.
The old star background population  contributes to
the MACS0717 SED by a bb flux corresponding to a temperature of 10$^4$ K.

\begin{table}
\centering
\caption{Modelling  the  line ratios (H${\beta}$~=1) observed by Masters et al (2014)}
\tiny{
\begin{tabular}{cccccccccccccccc} \hline  \hline
\ WISP   &   z   &  [OII]       &  [OIII]      & H${\alpha}$~   & [NII]     & H${\beta}$~$_{obs}$  $^1$ \\
\        &       &  3727+       & 5007+        &       & 6583      &                 \\ \hline
\ 173-205& 1.444 & 1.19         & 6.27         & 2.97  & -         & 7.5(0.4) 10$^{-18}$\\
\  m1    &   -   & 1.1          & 6.3          & 3.03  & 0.19      &-  \\
\ 9-73   & 1.454 & 1.88         & 5.3          &4.     & 0.51      & 31.6(1.3) 10$^{-18}$\\\
\  m2    & -     & 2.0          & 5.4          & 3.2   & 0.40      &-  \\
\ 25-53  & 1.486 & 2.13         & 2.55         & 2.8   & 0.51      & 31.0(0.9) 10$^{-18}$\\\
\  m3    & -     & 2.0          & 2.4          & 3.    & 0.42      & -\\
\ 46-75  & 1.504 & 2.8          & 1.36         & 3.05  & 0.45      & 6.6(1.1) 10$^{-18}$\\\
\  m4    & -     & 2.7          & 1.46         & 3.    & 0.37      &- \\
\ 170-106& 2.165 & 3.57         & 8.64         & 4.    & $<$0.35   & 3.7(0.4) 10$^{-18}$\\\
\  m5    & -     & 3.           & 8.32         & 2.9   & 0.13      &- \\
\ 64-210 & 2.177 & 1.49         & 6.76         & 4.2   & 0.75      & 20.1(1.4) 10$^{-18}$\\\
\  m6    & -     & 1.4          & 6.3          & 2.9   & 0.43     &- \\
\ 27-95  & 2.192 & 2.69         & 4.87         & 4.17  & 0.52     & 18.6(0.7) 10$^{-18}$\\\
\  m7    & -     & 2.67         & 4.57         & 3.1   & 0.4      &-\\
\ 90-58  & 2.212 & 2.15         & 2.65         & 3.79  & 0.54     & 8.7(0.3) 10$^{-18}$\\\
\  m8    & -     & 2.1          & 2.65         & 3.    & 0.5      & -\\
\ 56-210 & 2.304 & 1.55         & 6.25         & 2.8   & -        & 2.7(0.4) 10$^{-18}$\\\
\ m9     &  -    & 1.6          & 6.67         & 3.1   & -        &-\\
\ 206-261& 2.315 & 1.08         & 6.64         & 4.1   & $<$0.15  & 7.4(0.5) 10$^{-18}$\\\
\ m10    & -     & 1.1          & 6.67         & 3.    & 0.13 \\ \hline
\end{tabular}}

$^1$ H${\beta}$~ observed flux in  $\rm erg\, cm^{-2}\, s^{-1}$. The uncertainty appears in parenthesis

\end{table}

\begin{table}
\centering
\caption{Gas physical conditions and relative abundances calculated for the Masters et al. (2014) sample}
\tiny{
\begin{tabular}{lcccccccl} \hline  \hline
\  mod &  $\rm V_{s}$~ & $\rm n_{0}$  & $U$  & $\rm T_{*}$~       & $D$          & N/H        & O/H      & H${\beta}$~$_{mod}$   \\
\         & $\rm km\, s^{-1}$ & $\rm cm^{-3}$ & -    & 10$^4$ K  & unit$^1$      & 10$^{-4}$  & 10$^{-4}$ &unit$^2$       \\ \hline
\   m1    & 180  & 160  & 0.07 & 4.8       & 0.6          & 0.8        & 7.3      & 0.0036            \\
\   m2    & 170  & 120  & 0.06 & 4.2       & 0.5          & 1.5        & 6.8      & 0.001             \\
\   m3    & 140  & 150  & 0.03 & 4.0       & 1.0          & 0.7        & 7.       & 0.0046          \\
\  m4     & 140  & 150  & 0.01 & 4.5       & 1.0          & 0.3        & 6.       & 0.0043          \\
\   m5    & 320  & 200  & 0.12 & 7.5       & 3.0          & 0.1        & 8.5      & 0.1            \\
\   m6    & 250  & 210  & 0.1  & 5.6       & 2.0          & 0.7        & 7.5      & 0.065           \\
\   m7    & 230  & 100  & 0.05& 4.        & 1.0          & 0.85       & 7.6      & 0.0013         \\
\   m8    & 200  & 100  & 0.04 & 4.        & 2.0          & 0.8        & 6.4      & 0.0042        \\
\   m9    & 180  & 160  & 0.07 & 4.5       & 0.4          & 0.8        & 7.3      & 0.0018         \\
\   m10   & 180  & 160  & 0.07 & 4.8       & 0.6          & 0.6        & 7.6      & 0.0035          \\ \hline
\end{tabular}}

 $^1$ 10$^{16}$ cm ;
 $^2$  $\rm erg\, cm^{-2}\, s^{-1}$ ;

\end{table}

{\bf The WISP galaxies}.
The modelling  of  line ratios  from  the Masters et al sample of galaxies (Table 1)
with  prominent H${\alpha}$~ and/or [OIII]
lines, is shown  in Table 4.
 We  selected the galaxies corresponding to 2.8 $<$H${\alpha}$~/H${\beta}$~$<$ 4.2,
i.e. those which do not need  reddening correction
(Osterbrock 1974).  Modelling results  appear in Table 5.
 WISP 170-106  shows a relatively high  shock velocity  (320 $\rm km\, s^{-1}$)  and a relatively high $U$ (0.12) for $\rm T_{*}$~=7.5 10$^4$ K
 suggesting  a  fading activity.

{\bf J1000+0221S}.  Amor\'{i}n et al. (2014) report the spectrum observed from  this strongly
gravitationally lensed galaxy.
The  observed line ratios  appear in Table 6.
They are  similar to those   calculated by models presented in  GRID1 (table 4, mod 36) for AGNs and in GRID2
(table 1, mod 6  and table 2, mod 7) for SBs.
The AGN flux $F$ =3 10$^{11}$
photons cm$^{-2}$ s$^{-1}$ eV$^{-1}$ at the Lyman limit, is suitable to  Seyfert 2 AGN and
is higher by  $\sim$7-200  than  in the other galaxies  of the present sample.
 Amor\'{i}n et al. define this object as a  metal-poor star forming galaxy  based on the
low O/H and low $U$ deduced by empirical +photoionization method.
Our  models  show different results because
of the coupled effect of photoionization and shocks. This would  lead to
a different  mass-metallicity- SFR relation.

\begin{table}
\centering
\caption{Comparison of model results with the data for J1000+0221S
(Amor\'{i}n et al. 2014)}
\tiny{
\begin{tabular}{lcccccccc} \hline  \hline
\                        & obs        & M$_{AGN}$ & M1$_{SB}$ & M2$_{SB}$ \\ \hline
\        z               & 3.417      &     -     &     -     &       -     \\
\ [OII]3727              &$<$0.3      & 0.3        &$<$0.5    &0.43     \\
\ [NeIII]3869           & $<$0.2     & 0.2        & 0.25     & 0.4    \\
\ [OIII]5007+4959        & 5.9$\pm$2.6& 4.7        &5.5       &8.45     \\
\ H${\beta}$~                    & 1          & 1          & 1        & 1      \\
\ H${\beta}$~$_{mod}$ $^1$       & -          & 21.9       &0.13      & 2.58   \\
\ $\rm V_{s}$~ ($\rm km\, s^{-1}$)             & -          & 200        &100        & 200  \\
\ $\rm n_{0}$  ($\rm cm^{-3}$)            & -          & 200        &100        & 200   \\
\ $F$  $^2$                & -          & 3          & -         & -     \\
\ $\rm T_{*}$~(10$^4$K)           & -          & -          & 5         & 5    \\
\     $U$                & -          & -          &0.1         & 1    \\
\ $D$ (10$^{19}$ cm)    & -          & 1          & 1         & 1    \\
\   O/H (10$^{-4}$)      &  -           & 6.6        & 6.6      & 6.6  \\
\ N/H  (10$^{-4}$)       & -          &0.9          &  0.9     & 0.9  \\ \hline
\end{tabular}}

$^1$ in $\rm erg\, cm^{-2}\, s^{-1}$ ;
$^2$ in 10$^{11}$ photons cm$^{-2}$ s$^{-1}$ eV$^{-1}$ at the Lyman limit

\end{table}

\section{Concluding remarks}

The  high  starburst temperature  ($\rm T_{*}$~$\sim$ 2 10$^5$ K)
suggesting activity  in  the two galaxies CDFS-695 (van Dokkum et al. 2005) and  A773 (Belli et al 2013)
is not seen in the  other galaxies of the modelled samples.
Activity in terms of a  Seyfert 2  AGN  flux is found in
the J1000+0221S galaxy (Amor\'{i}n et al 2014).
The ionization parameter  $U$ calculated for the galaxy bulk
 shows  larger scattering than $F$
because $U$ depends on the hot source radiation flux  as well as on  the ISM  conditions
So it is less significant to activity than
$F$ which depends  mainly on the intensity of the  non-thermal radiation
flux from the AGN.

By detailed modelling of the spectra we have found  O/H  close to solar  (Allen 1976,
Anders \& Grevesse 1989).
N/H are generally lower than solar, in agreement with the results
obtained for galaxies at those redshifts.

\section*{Aknowledgements}
I am grateful to the referee for useful comments.

\section*{References}

\def\par\noindent\hangindent 18pt{\par\noindent\hangindent 18pt}
\tiny{
\par\noindent\hangindent 18pt Allen, C.W. 1976 Astrophysical Quantities, London: Athlone (3rd edition)
\par\noindent\hangindent 18pt Amor\'{i}n, R. et al. 2014, ApJL, 788, L4
\par\noindent\hangindent 18pt Anders E., Grevesse N. 1989, Geochim. Cosmochim. Acta, 53, 197
\par\noindent\hangindent 18pt Angeloni, R., Contini, M., Ciroi, S., Rafanelli, P. 2007, AJ, 134, 205
\par\noindent\hangindent 18pt Belli, S., Jones, T., Ellis, R.S., Richard, J. 2013 ApJ, 772, 141
\par\noindent\hangindent 18pt Contini, M. 2014, A\&A, 564, 19, Paper I
\par\noindent\hangindent 18pt Contini, M. 2013, arXiv1310.5447
\par\noindent\hangindent 18pt Contini, M., Cracco, V., Ciroi, S., La Mura, G. 2012, A\&A, 545, 72
\par\noindent\hangindent 18pt Contini, M. \& Prialnik, D. 1997, ApJ, 475, 803
\par\noindent\hangindent 18pt Contini, M. \& Viegas, S.M.A. 2001a, ApJS, 132, 211, GRID1
\par\noindent\hangindent 18pt Contini, M. \& Viegas, S.M.A. 2001b, ApJS, 137, 75, GRID2
\par\noindent\hangindent 18pt Kennicutt, R.C. 1998, ARA\&A, 36, 189
\par\noindent\hangindent 18pt Kriek, M., van Dokkum, P.G., Labbe', I., Franx, M., Illingworth, G.D.,
Marchesini, D., Quadri, R.F. 2009, ApJ, 700, 221
\par\noindent\hangindent 18pt Kriek, M. et al. 2006, ApJ, 649, L71
\par\noindent\hangindent 18pt Maiolino, R. et al. 2008, A\&A, 488, 463
\par\noindent\hangindent 18pt Masters, D. et al. 2014,  ApJ, 785, 153
\par\noindent\hangindent 18pt Osterbrock, D.E. 1974 in  Astrophysics of Gaseous Nebulae. W.H.Freeman and Company,
San Francisco, Ed.
\par\noindent\hangindent 18pt van Dokkum, P.G., Kriek, M., Rodgers, B., Franx, M., Puxley, P. 2005, ApJ, 622, L16
}
\end{document}


The r-Process in Metal Poor Stars and Black Hole Formation

Abstract: Nucleosynthesis of heavy nuclei in metal-poor stars is generally thought to occur via the r-process because the r-process is a primary process that would have operated early in the Galaxy's history. This idea is strongly supported by the fact that the abundance pattern in many metal-poor stars matches well the inferred Solar r-process abundance pattern in the mass range between the second and third r-process abundance peaks. Nevertheless, a significant number of metal-poor stars do not share this standard r-process template. In this Letter we suggest that the nuclides observed in many of these stars are produced by the r-process, but that it is prevented from running to completion in more massive stars by collapse to black holes before the r-process is completed, creating a "truncated r-process," or "tr-process." We find that the observed fraction of tr-process stars is qualitatively what one would expect from the initial mass function, and that an apparent sharp truncation observed at around mass 160 could result from a combination of collapses to black holes and the difficulty of observing the higher mass rare earths. We test the tr-process hypothesis with r-process calculations that are terminated before all r-process trajectories have been ejected. We find qualitative agreement between observation and theory when black hole collapse and observational realities are taken into account.

\section{Introduction}
The r-process has been understood for many years to (1) synthesize
half the nuclides heavier than iron, (2) synthesize all of the nuclides
heavier than $^{209}$Bi, and (3) be primary in the sense that its nucleosynthesis
does not appear to depend on preexisting nuclides
\citep{b2fh,woosley94,wallerstein97,farouqi09}. In this latter
context, its production of the heavier r-process nuclides in many metal poor stars
appears to be very uniform \citep{sneden03a}, and it also produces relative
r-process abundances that are essentially the same as those found in more modern stars.
A standard interpretation of the r-process suggests that it occurs in the neutrino
wind that emanates from core-collapse supernovae
\citep{woosley94, takahashi94,farouqi09}, although that interpretation is not
without its issues \citep{qian96, meyer98, hudepohl10a,hudepohl10b, fischer10, roberts10}.
Even with these difficulties,
though, there are a sufficient number of uncertainties, e.g., the possibility of
sterile neutrinos \citep{mclaughlin99, mclaughlin03}, that this r-process site might turn out to be viable.

At the same time, other metal poor stars exhibit some of the features of r-process
nucleosynthesis, but their abundances represent a surprisingly poor match with
the ``standard'' abundance template \citep{aoki00, honda07, roederer10a},
perhaps most frequently identified as that observed in
CS-22892-052 \citep{sneden03a}. A recent paper \citep{roederer10a}
has summarized the situation that exists for metal-poor stars. The data in that paper seem
to suggest that a distribution of r-process abundances exists in metal poor stars,
with some resembling the standard abundance set, but with a significant fraction
of stars having abundances that do not match the standard template. These latter
stars appear to favor the lighter r-process nuclides at varying levels, and many
seem to have abundance patterns that terminate around Dy, that is, around mass 160.

In this Letter we point out that the abundance patterns observed for the stars that
do not fit the standard r-process template could be produced by stars that are
sufficiently massive that their core collapse first produce neutron stars
but that the infall that occurs following the formation of the neutron star subsequently
results in collapse to a black hole - so called ``fallback supernovae.''
This class of stars spans a mass range from roughly 25 to 40 solar
masses \citep{heger92} for low-metallicity stars.
When the neutron star collapses to a black
hole the r-process enrichment of the interstellar medium would cease,
terminating either when the r-processed regions were swallowed by the black hole,
or when the electron antineutrinos fell below the event horizon \citep{sasaqui05}.
Thus, this truncated r-process, or tr-process, nucleosynthesis would terminate at
different stages of that process,
depending on the precise time at which the black hole prevented
further r-process production or emission of nuclides into the interstellar medium.
We suggest that the delayed collapse to the black hole,
combined with another effect, namely, the difficulties in observing the higher mass
rare earths, could produce the cutoff in the r-process distributions observed around 160 u.

The general scenario we envision merely assumes that mass elements that predominantly produce the
lighter r-process nuclei are ejected before the mass elements that produce the heavier r-process nuclides.
Any setting within a core-collapse supernova that satisfies this condition would allow for
a truncated r-process.  To evaluate this idea quantitatively, however, we apply the
neutrino-driven wind model for the r-process.  Although, as noted above, this model
is not without issues; it is plausible and well discussed in the literature and makes a good setting for discussing
the truncated r-process.

\section{Neutrino Wind Model of the r-Process}
In the neutrino wind scenario of the r-process, neutrinos from the nascent neutron
star heat material at the neutron star surface and drive it away in a wind.  As a wind
element leaves the neutron star surface, it expands and cools.  Nucleons in the
wind element assemble into heavier nuclei, which serve as the seeds for the subsequent
r-process \citep{woosley94}. Importantly, the abundance pattern for the seed nuclei
is well-described by a quasi-equilibrium, which, for the entropy and neutron richness
characteristic of neutrino-driven wind envirnoments, has a relatively sharp
abundance peak near A$\sim$100.  As the temperature in the wind element falls further,
charged-particle reactions freeze out.  If a high neutron density accompanies
the seed nuclei, neutron captures and beta decays successively promote
the seed nuclei to higher mass.  This upwards flow in mass slows
at the neutron closed shells at 82 and then 126 neutrons, and
thereby produces the r-process abundance peaks at A$\sim$130 and 195 u
\citep{woosley94, wanajo02}.

For what follows, it is essential to note that the neutrino-driven wind leaves the
star over a period of several seconds.  In the standard scenario
\citep{meyer92, woosley94, takahashi94}, the wind elements that leave the star early
have lower entropy and a lower degree of neutron richness.  As the neutron star
evolves, the wind elements that leave later have higher entropy and a greater
degree of neutron richness.  It is in these latter wind elements that the heaviest
nuclei are made.  The lower-entropy, less neutron rich wind elements that leave
the neutron star earlier produce ligher nuclei.  If the neutron star collapses to a
black hole after the lower-entropy elements leave the star but before the higher-entropy
do, then the r-process has been truncated, and the abundance pattern will be dominated
by lighter r-process nuclei.

Current supernova models do not naturally produce the parameters
required to produce a successful r-process, most notably the entropy, but
also the requisite neutron richness.  While it may in fact be the case that
core-collapse supernovae are not the site of the r-process, it may also be that
current supernova models are not sufficiently detailed to provide an adequate
description of the r-process conditions in core-collapse events.
For example, results from multidimensional hydrodynamics calculations suggest that the
instabilities resulting from those calculations may ultimately be shown to produce
the entropy required for making the heaviest r-process nuclei \citep{burrows06}.
Interestingly, these instabilities may generate the neutron star
kicks \citep{burrows06, guilet10} that give rise to high space velocities of some
pulsars \citep{arzou02}.
It will also be necessary
to include all of the detailed neutrino physics \citep{burrows02, roberts10, liebendoerfer05}
to adequately characterize the neutron richness in matter expelled from the supernova.

Future work will determine whether more advanced core-collapse supernova models will provide the
conditions needed for the r-process.  For the purposes of the present paper,
we will simply assume that
neutrino-driven winds in core-collapse supernovae are at least one of the sites of
r-process nucleosynthesis.

\section{Model Calculations}
To study the truncated r-process in the neutrino-driven wind,
we applied the basic idea behind the study of \citet{woosley94},
which assumed that the r-process occurred in the neutrino wind from a core collapse supernova. In that
study, a succession of 40 ``trajectories'' (that is, thin shell wind elements), all assumed to
have originated deep within the (assumed spherically symmetric) star, but
having initially different conditions of density,
temperature, entropy, and electron fraction, were emitted into the
interstellar medium, thus contributing
to the total r-process nucleosynthesis. The bubble
was evolving in time, so that the conditions under which the individual trajectories were
processed changed with the identity
of the trajectory. We assumed that the different
trajectories were emitted
from the star successively, but ceased to be emitted when the collapse to the black hole
occurred. This would be consistent with \citet{woosley94}, who assumed successive
emissions of the trajectories to generate what turned out to be a good representation of
the Solar  r-process abundances.

For our r-process calculations, we used a network code
based on libnucnet, a library of C codes for storing and managing nuclear reaction
networks \citep{meyer07}.  The nuclear and reaction data for the calculations were
taken from the JINA reaclib database \citep{cyburt10}.
We performed calculations for trajectories
24 - 40 in the \citet{woosley94} hydrodynamics model.
For each trajectory, reaction network calculations were
performed for T$_9<$2.5 using initial abundances derived from the \citet{woosley94}
results.  Our calculations were simplified by assuming an initial abundance of massive
nuclei from a single nucleus heavy seed with a mass A equal to the average mass at
T$_9$=2.5 and an atomic number derived from the average mass number, the Y$_e$, and
the neutron and alpha mass fractions at T$_9$=2.5 in \citet{woosley94}.  We justify
this by noting that the heavy-nuclide distribution is typically sharply peaked near T$_9$=2.5.
An adiabatic expansion was assumed for each trajectory in the nucleosynthesis, with entropy
constant within a trajectory but varying between trajectories, again consistent with
the approach of \citet{woosley94} for times at which the temperature T$_9<$2.5.
The material was taken to expand at constant velocity on a timescale consistent with
that derived from \citet{woosley94}.
For each trajectory the calculation was continued
until neutron capture reactions had ceased and the abundance distribution versus
mass number had frozen out. Our representation of the
full r-process abundances, shown as the dots in Figure \ref{woo_traj},
is not as good as that
achieved by \citet{woosley94},
but our calculations do produce the basic features of the
full r-process, namely the mass 130 and 195u peaks.

The simulation of \citet{woosley94} that produced a good r-process representation
began with trajectory number 24 and summed the nucleosynthesis yields from the remaining
16 elements since all of these trajectories would have been ejected into the interstellar
medium. We also began with trajectory 24, and
performed a mass weighted sum of the nucleosynthesis from the trajectories up to some
higher number trajectory to observe the resulting nucleosynthesis when the trajectories
beyond our maximum trajectory were assumed to be consumed by the collapse to a black
hole.
This produced a result similar to that of \citet{woosley94} when all trajectories
from 24 to 40 were included. The results are shown in Figure \ref{woo_traj}. There it can be seen that
truncating the r-process at increasing trajectory number does terminate the r-process nucleosynthesis
at increasingly higher nuclear mass. Note that although the curve representing trajectories
24 through 31 does reach the mass 195u peak, the abundance produced in that calculation
is nearly two orders of magnitude below that of the full r-process, which would render the
higher mass nuclides difficult to observe.
The abundances for that calculation would therefore appear observationally
to terminate at a mass of about 140 u.

In Figure \ref{ele_abun}  we compare several tr-process calculations with the derived elemental
abundance pattern in the metal poor halo star HD 122563. This low-metallicity star ([Fe/H]
= -2.7) is deficient in the heavy neutron-capture elements (Ba and heavier) relative to the
light neutron-capture elements (Sr through Cd) when compared with the scaled Solar
r-process pattern. The HD 122563 abundances are a very poor match to the scaled Solar
s-process pattern. Its abundance pattern matches the scaled Solar r-process pattern
better up to an atomic mass of about 70 (mass of $\sim$174 u), but even this fit is
unsatisfying \citep{sneden83, honda06}. Stars like
HD 122563 may be candidates for enrichment by the tr-process. Figure \ref{ele_abun} demonstrates
that the tr-process predictions, while far from a perfect match to the individual
abundances, can reproduce the relative downward trend in abundance with increasing
atomic number seen in some metal poor stars. More exploration of tr-process calculations
is obviously needed, but the general trend is encouraging.

\section{Probability of Occurence of tr-Process Stars}
In principle, a test of the tr-process model would be provided by a large set of
data for metal poor stars that spans the masses of the nuclides produced in the
r-process from mass $\sim$80 to the heaviest observable nuclides that the r-process
synthesizes,  usually thorium. Unfortunately this is challenging due to the difficulty
of observing the higher mass rare-earth nuclides, starting at mass $\sim$160. Thus all
r-process events that terminate at masses between 160 and lead would appear to terminate
at mass 160, producing the effect observed by \citet{roederer10a}.

Any standard initial mass function (IMF) reflects the fact that less-massive stars are
more numerous than more-massive stars.  Assuming that stars from 8 to 40 solar
masses may produce comparable amounts of r-process material, we can use the IMF
to estimate the fraction of stars whose r-process may be truncated by collapse to a
black hole.  We adopt the \citet{salpeter54} IMF, for which dN/dm $\propto$ m$^{-2.35}$
for massive stars, where m is the stellar mass in units of the solar mass, and N is the
number of stars of a given mass per unit volume.  The ratio of the number of stars that
would be expected to collapse to black holes (25-40 solar masses) to those expected to
collapse to neutron stars (8-25 solar masses) is 0.13 in a well-sampled IMF.

While a number of r-process rich stars (here taken to mean stars with [Eu/Fe] $>$ +1.0)
have been discovered and studied in detail in the last two decades, these stars constitute
a relatively minor fraction of all metal-poor stars.  More unbiased samples
\citep{mcwilliam95, barklem05} find that stars with [Eu/Fe] $>$ +1.0 comprise $<$10\%
of all stars with [Fe/H] $<$ -2.0.  Figure 11 of \citet{roederer10a} suggests that stars
with [Eu/Fe] $<$ 0 are candidates for enrichment by a truncated r-process.  Difficulties
in detecting Eu in metal-poor stars with [Eu/Fe] $<$ 0 make it even more challenging to
estimate how numerous these stars are.  Ba is more easily detected and may be used to
represent Eu and other heavy elements in stars lacking s-process enrichment.  From the
Ba abundances in the large survey of \citet{barklem05} we estimate a lower limit of
$\sim$55\% of metal-poor stars as candidates for enrichment by a truncated r-process.
This is significantly higher than the 13\% derived above
assuming a Salpeter IMF, yet Figure \ref{ele_abun} demonstrates that a
tr-process may be one way to produce the heavy nuclides observed in
some metal-poor stars. The discrepancy in these fractions could indicate that additional
nucleosynthesis channels, such as those proposed by, e.g.,
\citet{qian08}, \citet{farouqi09}, \citet{wanajo11}, or \citet{nakamura11},
together with the tr-process may contribute to the chemical enrichment of the early
Galaxy.
Given that massive stars with short lifetimes will dominate the chemical
enrichment at early times, though, the tr-process may very well have been a major source of
heavy nuclei at these epochs.

\section{Other r-Process Issues}
As a possible additional benefit of the tr-process, we note that the stars that were
truncated at mass points 28 or less produced no nuclides in the mass 130 peak or beyond.
This would have the effect of boosting the yields of the lighter r-process nuclides
relative to the heavier ones, filling in the mass 110 to 120 u region. It has been known
for some time that r-process calculations that produce a mass 195 u peak underproduce
the nuclides in the 110-120 mass region. The r-process nuclides that were even
lighter than those were underproduced by an even larger factor, but this would
also be consistent with distributed cutoffs of the r-process nucleosynthesis. Thus,
again guided by the IMF, but along with the mass dependent cutoff times, the tr-process
might potentially solve this troublesome aspect of current r-process analyses
\citep{kratz07,farouqi09}.
Conclusions in this regard, however, must
await more detailed analyses.

An obvious test of the tr-process model could occur from renewed effort to observe the
higher mass rare earth elements in very metal poor stars with low [Eu/Fe] ratios. If
elements in that mass region can be observed, and the tr-process is correct, additional
observations would be expected to map out the black hole collapse times over that region.
While no doubt challenging, a relationship between observations and black hole collapse
times would be of great interest.
\acknowledgments
RNB's work has been supported under the auspices of the Lawrence Livermore
National Security, LLC (LLNL) under Contract No. DE-AC52-07NA27344; MF's
by National Science Foundation grant PHY-0855013; BSM's by NASA grant
NNX10AH78G; YM's by the NEXT Program of JSPS and CSTP (GR098);
TK's by Grants-in-Aid for Scientific Research of JSPS (20244035),
Scientific Research on Innovative Area of MEXT (20105004), and Heiwa Nakajima
Foundation; and IUR's by the Carnegie Institution of Washington through the
Carnegie Observatories Fellowship. This is document LLNL-JRNL-491647.

Solar and Stellar Astrophysics (astro-ph.SR); Nuclear Theory (nucl-th)