Direct Evidence for an Enhancement of Helium in Giant Stars in Omega Centauri

Abstract: The double main sequence identified in the globular cluster Omega Centauri has been interpreted using isochrones to indicate a large variation in the abundance of helium. If true, a helium enhancement carries strong implications for the chemical and stellar evolutionary history of this cluster. However, only indirect measures currently support this conjecture. We report the discovery of a variation in the line strength of the near-infrared He I 10830 A transition in twelve similar red giants in Omega Centauri observed with PHOENIX on Gemini-S. Abundances of these stars derived from Magellan/MIKE spectra taken at Las Campanas Observatory show that the helium transition is not detected in the most metal-poor population ([Fe/H] < -1.8), yet is present in the majority of stars with [Fe/H] >= -1.8. These observations give the first direct evidence for an enhancement of helium in Omega Centauri. The appearance of helium appears better correlated with increased [Al/Fe] and [Na/Fe] abundances than as a function of [Fe/H], giving observational support to the presence of high-temperature H-burning in a prior generation of stars.

\section{Introduction}

Omega Centauri ($\omega$ Cen) is a globular cluster long known
to have a color spread along the red giant branch (Woolley 1966; Geyer
1967) indicative of a spread in the [Fe/H] abundances as originally
identified in cluster RR Lyrae variables (Freeman \&\ Rodgers 1975).
Among subgiants too, four (and possibly five) subgroups of different ages and
metallicities have been identified (Villanova et al.\ 2007; Bellini et al.\ 2010).
The discovery (Anderson 1997; Bedin et al. 2004)  of multiple stellar
populations on the main sequence in a globular cluster, first
identified in $\omega$ Cen was possibly
a singular anomaly. But another globular cluster, NGC 2808, displays multiple branches
on  the  main sequence  (Piotto et al.\ 2007).   These results, based on unevolved dwarf stars,
present firm evidence for multiple episodes of star formation in a
globular cluster, and  upset the long-held
paradigm that a globular cluster consists
of a coeval population of stars (Piotto 2009).

Norris (2004) offered an explanation for the 'double' main sequence in
Omega Cen based on stellar models, as resulting from a
substantially enhanced helium abundance in the 'blue' main sequence
($\Delta$Y=+0.13; Y=0.38), over the primordial abundance
(Y=0.25) in the 'red' sequence.  Subsequently, the blue main-sequence
of stars in $\omega$ Cen was discovered to be more metal-rich
than the red main sequence ([M/H]$_{blue}$=$-$1.26 vs
[M/H]$_{red}$=$-$1.57; Piotto et al.\ 2005). The cluster NGC 2808 shows a similar
effect (Bragaglia et al.\ 2010a).

An enhanced helium abundance would have a number of ramifications.
Stellar models suggest that  excess helium could account for many characteristics of
globular cluster color-magnitude diagrams:
an extended horizontal branch (HB), gaps in the HB,
RR Lyr periods and period distributions, and the ratio
of blue to red HB stars (D'Antona \& Caloi 2004, 2008; Norris 2004;
Catelan 2009; Gratton et al.\ 2010).
It is currently thought that a helium
enhancement and other chemical inhomogeneities might result
from self-enrichment (pollution) of a second generation of cluster stars
by the processed material from the first generation stars. Various sources
for excess helium have been proposed: production by first generation asymptotic branch
stars (AGB) of intermediate mass,  or massive hot stars; and Type II supernovae.
However models
can not uniquely  identify  the source of the  enrichment;
concerns remain because {\it all} of the processed first generation
material is required to produce enrichment, and the material from large numbers of stars
must be homogenized (Norris 2004; Piotto et al. 2005).  Omega Cen
appears to have an anomalously blue HB (D'Cruz et al.\ 2000; Momany
et al.\ 2004) which has been attributed (Norris 2004; Piotto et al.\ 2005)
to the progeny of the blue main sequence.
However,  AGB models (Choi \& Yi 2008) can not produce the
quantity of helium in first generation stars as required for Omega
Cen.

While there is great interest in the possibility of
a helium enhancement in these multiple populations, only
indirect methods have been used to date to infer this enhancement
in globular cluster stars. The observable quantities rely principally on the appearance
of the color magnitude diagrams,  comparison with
stellar isochrones,  and the composition of elements
other than helium, in particular Na and Al, which signal the products of high-temperature
H-burning (Densikov \& Denisenkova 1989; Langer et al.\ 1993) including
CNO, NeNa, and MgAl cycles that lead to helium production (cf. Gratton et al.\ 2001; Gratton et al.\
2004; Bragaglia et al. 2010b).

Helium is a difficult element to detect spectroscopically
in cool stars because high temperatures are generally required to
produce the optical transitions.  Hot stars, such as the HB objects,
where helium might appear, are subject to diffusion and preferential
settling of elements making the determination of abundances ambiguous
(Da Costa et al.\ 1986; Behr et al.\ 1999).  However, the infrared \ion{He}{1}
transition at 10830\AA, has been detected in metal-poor red
giant stars in the field (Dupree et al. 1992, 2009) as well
as in the globular cluster M13 (Smith et al. 2004).
The strength of the helium absorption is independent of [Fe/H] metallicity
in the field-star sample (Dupree et al.\ 2009).  A calculation of the 10830\AA\
profile in the quiet sun with PANDORA (Avrett \& Loeser 2008) demonstrates
that a change in the helium abundance from Y=0.28 to Y=0.4 increases
the equivalent width of the absorption by a factor of 3. Thus a simple measure of
the distribution of equivalent widths in stars within a restricted region  of
temperature and luminosity could reveal a variation
in helium strength. We report the results of such a study here of
stars in the globular cluster $\omega$~Cen.  Because this cluster is believed to have
substantial helium enhancement (Norris 2004), it is an optimum target.

Spectroscopic observations of $\omega$ Cen giants at Gemini-S (for the
near-IR \ion{He}{1} transition)
and Magellan/Clay (for optical spectra) are
presented in Section 2; the derivation of abundances from Magellan spectra is
addressed in Section 3; results are discussed in Section 4, and final
comments in Section 5.

\section{Observations and Reduction}
Selection of target stars in $\omega$ Cen
was guided by the fact that the \ion{He}{1} $\lambda$10830 transition appears
in metal-poor field giants with $T_{eff} \ge 4500K$ where the stellar
chromospheres reach temperatures (T $\sim$ 10,000K)
high enough to produce the helium line (Dupree et al.\
2009).  A restricted region of the red giant branch   was
selected to fit the temperature criterion.   The target region
spanned 0.05 in $B-V$ color and  was limited to a range
of $\sim$0.5 magnitudes in V between 12.95
and 13.48 (Fig. 1).  This selection avoids asymptotic giant branch stars
and minimizes or eliminates any temperature or luminosity dependencies of
the line strength. In addition, horizontal branch stars
are avoided to eliminate objects where elemental diffusion might
occur.   All of the target stars have a membership probability of
98 to 100\% based on proper motion studies (van Leeuwen et al.\ 2000),
and our radial velocity measures, discussed later,
are consistent with membership as well.

\subsection{Gemini-S/PHOENIX Infrared Spectra}

Classical time was awarded (Program GS-2010A-C-7) for 3 nights  on
PHOENIX\footnote{See http://www.noao.edu/ngsc/phoenix/phoenix.html}
mounted on  Gemini-S (Hinkle et al.\ 2003)
to obtain high-resolution spectra of the \ion{He}{1} 10830\AA\ line.
Three clear nights of observations  were obtained
beginning 24 March 2010 UT.
The PHOENIX setup consisted of a slit width of 4 pixels yielding a
spectral resolution of $\sim$50,000. The order-sorting filter J9232 was
selected which spans 1.077--1.089 $\mu$m, and allows access to the
\ion{He}{1} 1.083  $\mu$m line which is shifted by $\sim$+230 km~s$^{-1}$\
(equivalent to 8.3\AA) due to the radial velocity of the $\omega$~Cen cluster.
Standard  procedures were followed by acquiring target spectra using a nodding
mode  (A-B) with a spatial separation of 5 arcsec.
Standard stars and telluric comparison objects were observed with a nodding pattern using a
separation of 8 arcsec to avoid contamination of the targets by a residual
image. Exposure times for most of  the $\omega$~Cen giants
totaled 2 hours, and 1 hour for two of the brightest targets.
To verify the grating setup, a spectrum of a
nearby star with strong P Cygni-type emission at 10830\AA\ was
obtained.      Because standard comparison lamps have a
sparse wavelength pattern in this near-infrared region, we observed a bright
Population I K giant containing many securely identified
narrow absorption lines in order to determine the wavelength scale.
Throughout the night,  spectra
of fast-rotating B stars were obtained at various air masses  to
monitor the strength of the telluric water vapor features.  Flat
exposures and dark exposures were taken at the end of each night.  The
twelve red giants comprising the targets in $\omega$ Cen are listed in Table 1.

A technical issue with PHOENIX became apparent during our run at
Gemini-S. During the periods of good seeing, ($\sim$0.4 arcsec), the
spectrum appeared double in the spatial direction,
leading to slightly degraded image quality.   The source of this doubling is
not now known.\footnote{The doubling phenomenon has been observed at Gemini-S since
PHOENIX began operation on Gemini-S in 2002
(B. Rodgers, private communication). Two possibilities have been
suggested to us: the telescope was slightly out of  focus on the slit
(K. Hinkle, private communication) or an astigmatic collimator whose effects become
substantially more pronounced at longer wavelengths as shown in
PHOENIX spectra obtained at the KPNO 2.1 coud\'e feed in 1998 (C. Kulesa, private
communication).  We assume this was not noticed previously since, in
queue mode, PHOENIX was not likely to be used during optimum seeing
conditions.}
Other than a moderate reduction in the S/N, this issue does not affect
any of our results.

The raw images were corrected for cosmic rays using L.A.Cosmic (van
Dokkum 2001). After this, the nodded image pairs were subtracted from each
other. These subtracted images were divided by a dark-corrected,
normalized flat field, followed by optimal extraction of one-dimensional
spectra. Owing to the paucity of arc lines in the region of the 10830
\AA\ line, the spectra were wavelength calibrated using observations of
strong-lined K giants with known radial velocities. The individual
calibrated spectra were combined with proper variance weighting.
The spectra were normalized by excluding the region between
10833--10844\AA\
where the Doppler-shifted \ion{Si}{1} line ($\lambda$10827.09)
and \ion{He}{1} line ($\lambda$10830.34) occur, and  fitting a second-order
cubic spline to the spectra.  Spectra of the helium line region in the
target giants
are shown in Fig. 2.

Equivalent widths of the \ion{He}{1}
absorption were measured from the unbinned data.  Spectra were not
corrected for telluric absorption, but the equivalent width of a
water vapor line, if blended with helium,
was subtracted from the total helium absorption.
Our spectra of rapidly rotating B stars, taken as standard stars,
display a correlation between the equivalent width of the strongest
water vapor line in this region, $\lambda_{air}$=10832.109\AA\
(Breckinridge \& Hall 1973) and the 10838.034\AA\ line which in some
stars blends with  the helium transition.  A straight line represents
the correlation with a $\pm$2m\AA(1$\sigma$) error of the fit. This correlation
was used to estimate the strength of 10838.034\AA\ in the target stars based on the
observed strength of 10832.109\AA.  Water vapor lines
are marked in Fig. 2 for the $\omega$ Cen giants,  and a sample of the
telluric standard stars is shown in Fig. 3. The water vapor
transition at 10838.034\AA\  varied in strength from 4 to 18 m\AA.
Where helium is not apparent in our targets,
the equivalent width in Table 1 represents an upper limit to the strength of a
non-observed helium line, and generally corresponds to the
equivalent width of the water vapor transition.

Radial velocities of the targets were determined from the high resolution infrared
spectra, based on fitting  the \ion{Si}{1} photospheric absorption
line  with a gaussian profile.  We also measured central wavelengths
of 10-13 \ion{Fe}{2} lines for each star from the MIKE spectra  for the abundance
analysis discussed later, and velocity measures from the individual lines have
a dispersion of less
than 1.7 km~s$^{-1}$. The \ion{Fe}{2}  velocities agree with
the results from \ion{Si}{1} to better than 0.2 km~s$^{-1}$.
Values of the radial velocity are listed in Table 1.
The mean radial velocity of $\omega$ Cen has
a distribution peaked around
$+$233 km~s$^{-1}$\ (Reijns et al.\ 2006). Our targets are clearly
radial velocity members of the cluster confirming the membership
determination from proper motions (van Leeuwen et al.\ 2000). For the nine  targets in common with
Reijns et al.\ (2006), our values  agree to better than 2 km~s$^{-1}$, with
the exception of LEID 54064 where we differ by 8.1 km~s$^{-1}$; this star has the
largest error of 4.7 km~s$^{-1}$\  reported by Reijns et al.\ (2006).

\subsection{Magellan/MIKE spectra}

Echelle spectra of the  PHOENIX targets were obtained at the Magellan/Clay telescope on 30 June - 2
July 2010 (UT) using the MIKE spectrograph (Bernstein et al.\ 2002)
with a 0.70$\times$5 arcsec slit.  The spectral resolution with this
setup ranges from $\sim$26,000 on the blue side ($\lambda\lambda$3350--5000)  to $\sim$36,000 on
the red side ($\lambda\lambda$4900--9300).  Total exposure
times on individual  targets varied from 40 minutes to 1 hour, generally divided between
two exposures.  An IDL pipeline developed by S. Burles, R.
Bernstein, and J. S. Prochaska\footnote{ See
  http://www.lco.cl/telescopes-information/magellan/instruments/mike/idl-tools}
was used to extract the spectra.  This efficient pipeline evaluates
the detector gain, processes
the milky and order-tracing flats, obtains the Th-Ar arc solution,
constructs a slit flat, and extracts a spectrum including sky subtraction.
IRAF procedures
were used to flatten individual echelle orders, and combine multiple exposures.
The signal-to-noise ratio ranged from 100 to 150 in the combined spectra, depending on
position in the echelle blaze. Spectra were  normalized
with  a cubic spline that well fit the  standard white
dwarf LTT 9491. After normalization,
individual spectra were combined.  Lines of \ion{Fe}{2}, \ion{Na}{1},
and \ion{Al}{1}  selected for abundance
determinations occur  in
8 echelle orders drawn from both blue and red sides of MIKE.
The equivalent widths were measured using the IRAF routine {\it splot}
by fitting the profiles with a gaussian curve on the
continuum-normalized spectra. In some cases, where blending was present, the
blended profile was deconvolved using multiple gaussian profiles in
order to extract the equivalent width of the line.

\section{Determination of Stellar Abundances}

Abundances of Fe, Al, and Na  were derived from Magellan/MIKE spectra
of the target stars.
Following the precepts of Kraft and Ivans (2003), we focused on
deriving Fe abundances from \ion{Fe}{2} lines since this ion dominates
in cool giants.  The selection of \ion{Fe}{2} lines
and oscillator strengths was taken from Fulbright (2000).  Of the
26 \ion{Fe}{2} lines listed by Fulbright, we chose 10-13 to be measured
in each star, aiming to find those near the center of the echelle blaze from each
order where the signal-to-noise is highest.

The abundance analysis was restricted to lines with reduced
equivalent width $log$(EW/$\lambda$)$<-4.5$.  In addition, we only used lines with EW $> 9$ m\AA,
unless few stronger or no lines  were available for a given species and star.
Atomic parameters for all lines measured and individual equivalent widths are
listed
in Table 2.

We derived abundances by force-fitting equivalent widths in the 2010
version of MOOG (Sneden 1973\footnote{See
  http://www.as.utexas.edu/\~{}chris/moog.html}).
Effective temperatures for 11
of the stars were estimated using their $V-K$ colors  and the calibration of Alonso
et al.~(1999). The K-magnitudes were taken from the 2MASS All Sky
Catalog
(Cutri et al.\ 2003) which lie within 0.005 mag of the TCS color system
adopted by Alonso et al.~(1999). One star did not have a $K$ magnitude; we used its $B-V$
color instead. A reddening of $E(B-V) = 0.12$ mag was assumed (Harris 1996). These
photometric estimates of effective temperature have small random errors:
0.01 mag in $V-K$ corresponds to a change of only 10K; 0.01 mag in
$E(B-V)$ corresponds similarly to  30K. Systematic effects dominate,
and are estimated at less than 100 K.

Surface gravities were estimated using bolometric corrections from Alonso
et al.~(1999) for the effective temperatures already derived, a distance of
$4.8\pm0.3$ kpc (van de Ven et al.~2006), and an assumed mass of 0.8
$M_{\odot}$. The gravity is insensitive to reasonable uncertainties in
mass and temperature as they enter the calculation only logarithmically.
The largest source of error is the distance, which introduces a $1\sigma$
random error of approximately 0.05 dex.

To obtain chemical abundances, Castelli \& Kurucz (2004) stellar
atmosphere models were employed. These are one-dimensional, plane-parallel models
without convective overshooting. We used $\alpha$-enhanced models
([$\alpha$/Fe] = +0.4) that are generally appropriate for stars in a
globular cluster. Models for the relevant parameters of each star were
derived through interpolation. The effective temperature and surface
gravity were fixed by the photometry as discussed above. We derived the
metallicity [Fe/H] and microturbulence $v_t$ iteratively. To minimize
non-LTE effects, only \ion{Fe}{2} lines were used. The
microturbulence was determined in the usual manner
by requiring that there be no trend of abundance with reduced equivalent
width. The final adopted parameters for each star are given in Table
1.

Abundances for sodium and aluminum were derived using 2--4 weak lines per
species and the stellar parameters estimated from the photometry and the
Fe II lines.  Four lines of \ion{Al}{1} were  selected: $\lambda$6696.03, $\lambda$6698.67,
$\lambda$7835.32, and $\lambda$7836.13.  The \ion{Na}{1} analysis
typically used 5 lines:
$\lambda$4668.57 (in some stars), $\lambda$5682.65, $\lambda$5688.21,
$\lambda$6154.23,
and $\lambda$6160.75. For one star, LEID 63027, we also included the
Na  resonance lines at $\lambda$5889.96 and $\lambda$5895.94 and they
gave results consistent with those of weaker lines.
Wavelengths of the transitions and their {\it gf}-values were taken
from the compilation of Fulbright (2000).  The {\it gf}-values for Al
came from Jacobson et al.\ (2009).  For a subset of the stars,
we could only derive upper limits
to the abundance of aluminum. In these cases, the reported numbers are a
median of the limits for the individual lines.
Abundances are reported relative to solar values. The solar values in
Asplund et al.~(2005) were used for all elements except for iron, for which we assume
log $\epsilon$(Fe) = 7.52 (Sneden et al.~1991).  Fig. 4 shows
a section from the  MIKE spectra of two stars similar in T$_{eff}$, {\it log g},
and [Fe/H] yet differing substantially in Na and helium.

After our analysis was completed, Johnson and Pilachowski (2010) published
abundances for many red giants in Omega Cen based on HYDRA spectra
of lower resolution [R($\lambda$/$\Delta \lambda$)=18,000].
Their sample includes nine of our targets. For the  stars in common,
our values for $T_{eff}$ differ by a mean of 9K, and a maximum of
29K. The average difference in the values of {\it log g} amounts to 0.13.
These authors derived an iron abundance generally based on
a majority of \ion{Fe}{1} lines.  Some stars had no \ion{Fe}{2};
others had one, or at most two lines from \ion{Fe}{2}
available.  Fig. 5 shows the derived abundances for the
stars in common where an mean offset of $-$0.12$\pm$0.12  dex exists between our values
and those from Johnson \& Pilachowski. The maximum deviation of 0.2
dex from this offset is less than the scatter of 0.65 dex among other values
commonly found in the literature (see Fig. 5 of Johnson \& Pilachowski 2010).

\section{Results}

Helium was detected in 5 of the 12 red giants observed (see Fig. 2 and
Table 1), and the
equivalent widths of the lines differ from star to star. Considering the small range
of effective temperature and luminosity in these giants, the
equivalent width variation suggests a similar
variation in the helium abundance.  These helium lines, where
detected, have larger equivalent widths than the helium lines
found (Dupree et al.\ 2009) in metal-poor red giants with
comparable temperatures in the field (Fig. 6).  Additionally, where helium is
not found, the lines in $\omega$ Cen tend to be weaker than in the
field stars.  For the $\omega$ Cen giants,
no correlation is apparent between position in the color-magnitude diagram (CMD)
and the presence of helium or with the iron abundance (Fig. 7).  Our sample
is small, and more targets are necessary to develop a clear pattern
between helium line strength and the iron abundance, if, indeed it
exists.  In the metal-poor field giants, we also found no correlation
of helium strength with [Fe/H]. On the subgiant branch of $\omega$ Cen
also, Villanova et al.~(2007) found poor correlation  between [Fe/H]
and CMD location.  The helium line, when present,
may become stronger in the less luminous objects (Fig. 8).

\subsection{Overall Characteristics of the Presence of Helium}

The current view holds that increased helium abundance causes  the blue main sequence
of $\omega$ Cen (Norris 2004), and the gross characteristics of the
$\omega$~Cen
giants in the PHOENIX sample exhibiting
helium are consistent with both the fraction  and spatial distribution of
the blue  main sequence.      The
fraction of stars displaying helium, 42\%,  is comparable to the
fraction of blue main sequence
stars (25 to 35\%) considered to have enhanced helium and also enhanced
metals (Bedin et al.\  2004; Piotto et al.\ 2005). The spatial
distribution is also of interest.  The blue main
sequence stars have a higher spatial concentration near the cluster core
than the red main sequence (Bellini et al. 2009; Sollima et al.\ 2007). Within the half-mass
radius of the cluster ($\le$ 4.2 arcmin), the ratio of blue to red
main sequence stars ranges from 1.1 to 0.8; at 6 arcminutes
from the center the ratio is still high at 0.5; further out from the
center ($\ge$ 8 arcmin) the ratio falls to 0.4. Taking the
center of the $\omega$ Cen cluster (Anderson \& van der Marel 2010)
as RA = 13:26:47.24, Dec=$-$47:28:46.45 (J2000.0), our sample
of twelve giants is evenly split in radius at a distance of  6 arcminutes from the center.
Of the 6 stars within 6 arcminutes from the center, 3 have
helium or 50\%; this fraction falls to 33\% (2 out of 6 targets)
for the more distant objects. While our sample is small, these
fractions are consistent
with the spatial behavior of the blue:red ratio.
Thus, the percentage
detection and the spatial distribution
of the helium detections are consistent with the blue main sequence population.

The equivalent widths of the detected helium lines are stronger
than the upper limits of the  non-detections by factors of 3 to 10.
Judging from our calculations for the quiet Sun, the amount of this
enhancement could correspond to a change in helium abundance by more than a
factor of two. However, conditions in a red giant chromosphere are certainly
different from a dwarf star.  Outflows can increase the absorption
from the metastable lines in helium (Dupree et al.\ 2009).
Variability  in the   chromospheric structure  could
occur also, although it is not likely  to be large judging by the
lack of H-$\alpha$ wing emission,  and the weak \ion{Ca}{2}
emission reversals in the optical spectra of these giants.  Detailed
non-LTE semi-empirical calculations are underway to extract the
helium abundance, and the chromospheric models must be constrained by
H-$\alpha$ and \ion{Ca}{2} line profiles.

\subsection{Abundances of the Giant Stars}

The relation
between the abundance of iron [Fe/H] and the presence of helium
is shown in Fig. 9.  Stars with [Fe/H] $<$ $-$1.8 do not exhibit
helium; the majority of stars in the range
$-1.8 \leq [Fe/H] \leq -1.4$ display the  helium transition, and the
one metal-rich star in our sample at [Fe/H]=$-$1.15 does not show
helium.\footnote{Our abundance for LEID 54084, [Fe/H]=$-$1.79,
differs from the  Johnson \& Pilchowski (2010) value of
[Fe/H]=$-$1.5. Such a shift would not significantly change the distribution of helium
detections.}  On the main sequence, the metallicity, [M/H] corresponds to
$-$1.26 for the blue main sequence and $-$1.57 for the red main
sequence (Piotto et al.\ 2005).  An empirical transformation from [M/H]
to [Fe/H] was estimated by Villanova et al.\ (2007) from a comparison of
subgiant spectra to spectral syntheses for various values of [M/H].
The best-fitting synthesis for [M/H] was compared to values of [Fe/H]
for many subgiants in $\omega$ Cen. The correction, obtained from the mean difference
between [M/H] and [Fe/H] amounts to $-$0.11 dex applied to values
of [M/H]. For the blue and red main sequences, this leads to
[Fe/H]=$-$1.37 and $-$1.68 respectively. An additional
minor red branch  can be found
on the main sequence (Villanova et al 2007; Anderson and van der Marel
2010)  termed MS-a which has higher [Fe/H] ($-$1.1 to $-$0.6) than the
blue main sequence.  We do not find  a correlation between the
[Fe/H] distribution and the presence of the \ion{He}{1} $\lambda$10830 line.

A different indication of processed material can be found in the light element
abundances of Na and Al which signal high temperature H burning via
the CNO cycle (Denisenkov \& Denisenkova 1989; Langer et al.\ 1993) in a previous
stellar generation (Gratton et al.\ 2001).  In red giants,
the outer convective envelope will never reach the the H-burning shell
where the light elements are produced during evolution on
the red giant branch (Gratton
et al.\ 2004).  For enhanced light elements to be present, they
must have originated in a previous generation of stars.
The abundances of   [Na/Fe]  and [Al/Fe] are
correlated with each other, and clearly enhanced in the stars
displaying \ion{He}{1} lines (Fig. 10). Stars without a helium
detection have low values of [Na/Fe] and principally upper limits
on the [Al/Fe] abundance.   These enhancements of Na and Al  agree with
values from red giants in $\omega$ Cen (Norris \& Da Costa 1995;
Johnson \&
Pilachowski 2010). And the values of [Na/Fe] and [Al/Fe]  for the helium and non-helium giants
are similar to the
abundances found in NGC 2808 for two stars, one each from the blue and
red components of the main sequence respectively (Bragaglia et al.\ 2010a).
Thus the red giants in $\omega$ Cen exhibiting
the abundance pattern expected from
the high temperature H-burning processes are the majority of  red giants
displaying enhanced He lines.

\subsection{Evolutionary Connections}

The population structure in $\omega$ Cen is quite complex.
We now focus on the spectroscopic [Fe/H] values, and consider
how the helium detections fit in  with the continuity of the
color-magnitude diagram (CMD) and the various
population components.  We follow the ideas of Villanova et al.\ (2007),
to interpret the presence of helium in the $\omega$ Cen giants and
diagram these ideas in Fig. 11.
These authors evaluated three factors in presenting an evolution scenario:
(1) morphological continuity on the CMD; (2) population fraction;
and (3) spectroscopic metallicity of the stars.
The bifurcated main sequence (MS) principally represents two values
of [Fe/H], $-$1.68 (spanning 0.2 dex) and $-$1.37,
and a minority population with [Fe/H] = $-$1.1 to $-$0.6.  Detailed abundance study of
the subgiant branch (SGB) reveals  (Villanova et al.\ 2007)  four or five components
which peak at specific values of [Fe/H] and the abundances typically span
0.3 to 0.6 dex (FWHM) each, and range in [Fe/H] from $-$1.7 to
$-$0.6.  These too are marked in
Fig. 11 and denoted by SGB: A, B, C, and D.  Possibly an additional
component bridges the abundance range between B, C and D.  On the
red giant branch (RGB),  Johnson \& Pilachowski (2010) identified 4 components
which they denote (following Sollima et al.\ 2005), RGB: MP, MInt1,
MInt 2+3 and RGB-a where MP and MI stands for metal-poor and
metal-intermediate.  The abbreviation RGB-a denotes the anomalous
red giant branch (Lee et al.\ 1999; Pancino et al.\ 2000). To understand the evolutionary relationships
among these components, Villanova et al.\ (2007) suggest
'cross-over' occurs from the red main sequence (rms) to the
subgiant (SGB) components A, B, and C, and from the
blue main sequence (bms) to the subgiant subpopulations B and C,
and perhaps the intermediate component located between B, C, and D.

From the subgiant branch to the red giant branch, component SGB-A
feeds into the RGB-MP, and stars on SGB-D move to RGB-a. Thus the
absence of helium in giants for [Fe/H]~$<$~$-$1.8 is not
surprising following the route from the MS-rms to RGB-MP. The MS-rms
also overlaps
the metallicity of SGB-B,C and the component percentages suggest that
MS-rms evolves to SGB-B,C, and then to RGB-MInt1 (Villanova et al.\ 2007),
The $\omega$~Cen stars with helium span [Fe/H] = $-$1.8 to $-$1.45
in our sample.
Villanova et al.\ (2007) conclude, based on photometric
and spectroscopic information that the MS-bms evolves through SGB-B,C
which
then, based on the metallicity overlap, contributes to  RGB-MP and
RGB-MInt1,
and perhaps RGB-MInt2+3.  In their abundance
study of 855 red giants in $\omega$ Cen,
Johnson and Pilachowski (2010) found that  O-poor, Na/Al rich
stars occur in three out of the 4 populations of red giants in the cluster,
but not in the RGB-a cohort.  The intermediate population of the [O/Fe] poor
giants (RGB-MInt1 and RGB-MInt 2+3) shares attributes with the blue
main sequence: spatial distribution, metallicity, and percentages,
leading Johnson \& Pilachowski (2010) to suggest that the [O/Fe]
depletion (coupled with the Na/Al enhancement) signals those
giants whose progenitors belonged to the blue main sequence. The mix of
helium detections and non-detections between $-$1.8 and $-$1.45
appears consistent with the crossovers outlined above.

We are left with an apparently anomalous object, LEID 38269, which
is metal-rich, [Fe/H] = $-$1.16 and does not show helium. Its
evolutionary path currently places it on the RGB-MInt2+3, and
its metallicity would seem to prevent an origin on the metal-poor MS-rms
where helium is weak or absent.   However, age dating of
the subgiant population in $\omega$ Cen provides a clue to this object.

In addition to continuity of metallicity values,
Villanova et al.\ (2007)  note that the subgiant branch offers an
optimum region of the CMD to
determine ages by comparison to stellar isochrones.   Interestingly, they find  the broadest
dispersion in stellar ages for stars with  $-1.8 \le [Fe/H] \le -1.5$
which echoes earlier results from studies of the main sequence
turnoff stars by Stanford et al.\ (2006).
Over this metallicity range,
the estimated ages  span almost a factor of two, from 0.6 of the
age of the oldest stars  continuously ranging to the oldest stars in
the cluster.  Thus the mixed appearance
of helium in our targets at this metallicity range appears in
harmony with the span of ages found in subgiants.  At the highest iron abundance
([Fe/H] $\sim$$-$1.1), three subgiants in the Villanova et al.\ (2007) sample have ages
commensurate with those of the most metal-poor stars suggesting a
simple age-metallicity correlation does not exist in $\omega$ Cen.
This might occur as a result of
several star-forming episodes.  Material
from Type II supernovae creates one stellar generation which is
followed by the polluting effects from
longer-lived metal-poor intermediate-mass asymptotic giant branch stars on
a subsequent stellar generation (Johnson \& Pilachowski 2010).
LEID 38269 also does not show enhanced Na and Al,  and
would plausibly  belong to this older population. It
is the coolest and faintest star in our sample, and could form a part of
the 'anomalous' population in $\omega$ Cen.

\subsection{Concluding Comments}

The direct detection of a helium abundance variation and its
correlations with both [Fe/H] and light element abundances
of Al and Na in this small sample of red giants,  supports
the conjecture that helium is enhanced in the blue main sequence
component of the $\omega$~Cen cluster. It appears likely that the
helium-enhanced stars contain the products of high-temperature
hydrogen burning produced by  an earlier generation of stars.
Spectra of additional targets
would be useful to thoroughly understand the relationship between
[Fe/H] and the presence of helium.

The infrared helium line at 10830\AA\ appears to be a powerful
probe of the presence of helium in warm giant stars. Detailed non-LTE
studies of its formation under a variety of conditions are needed to
be able to extend the interpretation to a larger range of temperatures
and luminosities, and also to extract the absolute abundance of
helium.

This research has made use of NASA's Astrophysics Data System.
We thank the US Gemini Office for assistance in preparing Phase II,
and providing  travel support to Gemini-S.   We thank Steve Margheim
and Claudia Winge  for Phase II advice and support during
our classical run at Gemini-S. JS acknowledges support provided
by NASA through Hubble Fellowship grant HST-HF-51237.01,
awarded by the Space Telescope Science Institute, which is
operated by the Association of Universities for Research in Astronomy,
Inc., for NASA, under contract NAS 5-26555.

{\it Facilities:} \facility{Gemini:South
  (PHOENIX)}, \facility{Magellan/Clay (MIKE)}

Solar and Stellar Astrophysics (astro-ph.SR); Astrophysics of Galaxies (astro-ph.GA)


Atomic and Molecular Phases of the Interstellar Medium

Abstract: This review covers four current questions in the behavior of the atomic and molecular interstellar medium. These include whether the atomic gas originates primarily in cold streams or hot flows onto galaxies; what the filling factor of cold gas actually is in galactic regions observationally determined to be completely molecular; whether molecular hydrogen determines or merely traces star formation; and whether gravity or turbulence drives the dynamical motions observed in interstellar clouds, with implications on their star formation properties.

\section{Introduction}

Star formation occurs primarily in atomic and molecular interstellar
gas, although \cite[Keto (2002)]{Keto02} emphasizes that around a massive
enough star, even ionized gas can be accreted. In this review, I
discuss a number of recent advances in our understanding of the
origin and dynamics of that gas and its relation to star formation.

\section{Cold Accretion and its Limits}

The origin of this
gas lies in gas accreted onto galaxies during their formation and
subsequent evolution.
Perspectives on how that gas arrives have changed several times
recently.  Originally, it was thought that all gas accreted into the
halo of a galaxy the size of the Milky Way would be heated to the
virial temperature of over $10^6$~K, followed by gradual cooling and accretion.
However, a decade ago, numerical simulations started to suggest that a major,
possibly even dominant alternative was for dense
streams to accrete from the cosmic web, dense enough to cool quickly
when compressed and thus to never heat to the virial temperature
(e.g. \cite[Fardal et al.\ 2001]{Fardal01}, \cite[Birnboim \& Dekel
2003]{Birnboim03}, \cite[Kere\v{s} et al.\ 2005]{Keres05},
\cite[2009]{Keres09}, \cite[Ocvirk et al.\ 2008]{Ocvirk08},
\cite[Dekel \& Birnboim 2006]{Dekel06}, \cite[Dekel et al.\ 2009]{Dekel09},
\cite[2013]{Dekel13}, \cite[Goerdt et al.\ 2015]{Goerdt15}). However,
\cite[Nelson et al.\ (2013)]{Nelson13}, using the AREPO Voronoi mesh code
(\cite[Springel 2010]{Springel10}), argued that cold streams are
likely only important for small galaxies at high redshift (see Fig.~\ref{Nelson13}), while their
appearance at lower redshifts and for larger galaxies is a numerical
artifact of the smoothed particle hydrodynamics (SPH) method, as exemplified
in that paper by GADGET (\cite[Springel 2005]{Springel05}).
\begin{figure}[bth]
\begin{center}
\includegraphics[width=1.8in,angle=90]{MacLow_Nelson13.eps}
\caption{\label{Nelson13} Comparison of gas accretion rates for hot
  gas (left, with temperature above $T_c$ as shown in legend) and cold
  gas (right) as a function of halo mass, for SPH models (brown) and
  Voronoi mesh models (green).  The two methods agree for small masses
  typical of high redshifts, but the Voronoi mesh code predicts most
  gas will reach high temperatures for large halo sizes, disagreeing
  with the SPH results. Adapted from \cite[Nelson et al.\ (2013)]{Nelson13}.}
\end{center}
\end{figure}

Evidence for continuing accretion of gas in the modern Milky Way can
be gained from observations of high velocity clouds of neutral gas
falling towards the galactic plane at velocities over 70 km~s$^{-1}$.
Material falling in at slower velocities has been shown to have
Galactic chemical abundances, while the high velocity clouds have far
lower abundances, suggesting extragalactic origin, as reviewed by \cite[Putman et al.\
(2012)]{Putman12}. Simple analytic models modeled on the flow of water into a bathtub and
down the drain can guide our understanding of the cycle of gas into the
atomic phase from infalling ionized gas, and then out of it again
through star formation or expulsion of the gas in galactic outflows
(\cite[Bouch\'e et al.\ 2010]{Bouche10}, \cite[Dekel et al.\
2013]{Dekel13}, \cite[Dekel \& Mandelker 2014]{Dekel14}).

\section{Why a Cold, Dense, ISM is Mostly Hot}

The thermodynamic state of the atomic and diffuse molecular gas is controlled
by the balance between optically-thin radiative cooling and photoelectric heating by far ultraviolet
radiation.  The between around 100~K and $10^4$~K the
cooling curve is a sufficiently steep function of temperature to
produce a region of thermal instability (\cite[Field et al.\
1969]{Field69}, \cite[Wolfire et al.\ 1995]{Wolfire95},
\cite[2003]{Wolfire03}). In the case of pressures typical of the Milky
Way, this leads to a two-phase medium (solid line in
Fig.~\ref{Wolfire03}).  When pressures increase (dotted line), only
the cold, dense phase remains thermally stable.
\begin{figure}[bth]
\begin{center}
\includegraphics[width=2.5in,angle=90]{MacLow_Wolfire03.eps}
\caption{\label{Wolfire03} Thermal equilibrium curves at different
  radii in the Milky Way, compared to a typical pressure near the
  Solar circle (solid orange line), and compared to the pressure in the molecular ring or
  starburst galaxies (dashed orange line).  Orange circled points represent
  stable equilibria, where small perturbations of temperature (motions
  perpendicular to the green-dashed line, with increasing
  temperature towards the left) result in restoring
  increases in cooling or heating, as labeled. Purple circled point is
  unstable. Adapted from \cite[Wolfire et al.\ (2003)]{Wolfire03}.}
\end{center}
\end{figure}

\cite[Ostriker et al.\ (2010)]{Ostriker10}
argued that, in fact, this transition leads to self-regulation of star
formation. The pressure range within which a warm phase will be present is
determined by the far ultraviolet radiation field produced by stars.
If the vertical hydrostatic pressure exceeds that range, it will drive
nearly all the gas into the cold phase, promoting star formation, and
thus increasing the far ultraviolet field until a warm medium can be
reestablished. Thus, the combination of thermal and hydrostatic
equilibrium controls,  in the terms of the bathtub model, how
fast the bathtub drains.

While the point that high enough pressures produce a predominantly
cold temperature distribution is widely understood, what is perhaps less appreciated
is the resulting spatial distribution of the gas.  \cite[Gatto et
al.\ (2015)]{Gatto15} and \cite[Li et al.\ (2015)]{MiaoLi15} emphasize
that the filling factor of the cold, dense gas will be small in such a
high-pressure environment.  In regions with supernova explosions
(either from active star formation or even Type Ia supernovae from
older populations) the space in between the cold, dense clumps is
filled with hot, tenuous gas (see Fig.~\ref{Gatto15}), the third phase of \cite[McKee \&
Ostriker (1977)]{McKee77}.
\begin{figure}[bth]
\begin{center}
\includegraphics[width=2in,angle=90]{MacLow_Gatto15.eps}
\caption{\label{Gatto15} Temperature distributions on slices through
  periodic boxes of supernova-driven turbulence.  On the left, supernovae
  explode in random positions, while on the right, half the supernovae
  explode in density peaks, substantially increasing the loss of
  energy to radiative cooling.  On the left, pressures have reached
  high enough levels to eliminate the warm stable phase. Adapted from
  \cite[Gatto et al.\ (2015)]{Gatto15}.}
\end{center}
\end{figure}

Observational tracers of atomic and molecular gas will show that the
mass of the interstellar gas is almost entirely in cold, molecular
material, leading to the description of the region as a fully
molecular interstellar medium.  The hot gas, on the other hand, is
difficult to observe, as its density is low and temperature high
enough that it emits primarily in bremsstrahlung. However, it
still completely dominates the volume filling factor, and needs to be taken into
account in any model of high-pressure, fully molecular, interstellar
gas.  \cite[Joung et al.\ (2009)]{Joung09} demonstrated that following
the Kennicutt-Schmidt law for star formation in stratified disks
results in a three-phase medium, while \cite[Walch et al.\
(2015)]{Walch15} showed that increasing the star formation rate by a
factor of three already leads to a cold interstellar medium filled
with hot gas.

 \section{Molecules are a Consequence of Star Formation}

Observations have demonstrated a linear correlation between the
surface density of star formation and that of molecular hydrogen
(\cite[Rownd \& Young 1999]{Rownd99}, \cite[Wong \& Blitz
2002]{Wong02}, \cite[Bigiel et al.\ 2008]{Bigiel08}, \cite[Leroy et
  al.\ 2008]{Leroy08}).  This raises the question of whether this
  correlation implies a causal relationship.  Is molecular hydrogen
  formation required for star formation to occur?

It has more recently become evident that many tracers of high density
gas correlate just as well with the star formation rate as does
molecular hydrogen.  For example, the tracer molecule HCN was shown by
\cite[Gao \& Solomon (2004)]{Gao04} to also correlate with star
formation rate in both normal spiral galaxies and luminous and
ultraluminous infrared galaxies.  In the solar neighborhood, the
number of young stellar objects in a cloud correlates with the mass of
the cloud contained in regions with a K-band infrared extinction $A_K
> 0.8$ (\cite[Lada et al.\ 2010]{Lada10}).  \cite[Clark \& Glover
  2014]{Clark14} indeed found thresholds for star formation
in numerical simulations that reproduce this observed
threshold. However, the physical process producing this threshold is
{\em not} molecule formation, but rather dust shielding from heating
by the background far ultraviolet radiation.

Indeed \cite[Krumholz et al.\ (2011)]{Krumholz11} and \cite[Glover \&
Clark (2012)]{Glover12} have demonstrated that atomic line cooling is
as effective as molecular cooling in bringing gas down to star forming
temperatures of 10-20~K.  Figure~\ref{Glover12} demonstrates that the same
cold, star-forming gas appears in numerical simulations with full
molecular and atomic cooling; with only atomic cooling, but still including
molecular hydrogen formation; and with all molecular physics suppressed, but still
including dust shielding from far ultraviolet heating.  Only the
removal of dust shielding from the simulation suppresses star
formation.
\begin{figure}[bth]
\begin{center}
\includegraphics[width=4in,angle=90]{MacLow_Glover12.eps}
\caption{\label{Glover12} Distribution of temperature and pressure in
  models of turbulent, atomic gas including both atomic and molecular
  chemistry and cooling processes, as well as photoelectric heating
  modulated by dust shielding.  The bottom panel shows the full model,
  while the upper panels have successively less physics. The first
  panel from the bottom neglects molecular line cooling, the second
  also neglects molecular hydrogent formation, and the top also
  neglects dust shielding.  Gas reaches 10~K and high densities in all
  but the last of these models. Adapted from \cite[Glover \& Clark (2012)]{Glover12}.}
\end{center}
\end{figure}

Thus, it appears that the formation of molecular hydrogen does not control
star formation. Rather, it seems, star formation, or more generally
gravitational collapse of atomic gas, controls molecular hydrogen
formation.

\section{Molecular Cloud Dynamics Driven by Gravity Rather Than
  Turbulence}

Molecular cloud dynamics offer a useful window to their origins.  The
argument has been made that they form by the sweeping up of smaller
clouds in superbubbles or spiral arms (e.g.\ \cite[Tasker \& Tan
2009]{Tasker09}). An alternative proposal is that clouds form from
large scale gravitational collapse, as proposed by
\cite[Hartmann et al.\ (2001)]{Hartmann01} to explain the apparent
correlation of star formation across entire molecular clouds.

Models of the supernova-driven, stratified interstellar medium can be
used to examine these scenarios.  One example of such a model derives
from work by \cite[Joung \& Mac Low (2005)]{Joung05}, \cite[Joung et
al.\ (2009)]{Joung09}, and \cite[Hill et al.\ (2012)]{Hill12}.  It
includes both clustered and random supernovae, a background
gravitational potential from the stellar disk and the halo,
photoelectric heating and equilibrium ionization cooling, and magnetic
fields.  These models swept up gas into filaments that appeared
similar to large molecular clouds.  An example of such a filament is
shown in top panel of Figure~\ref{IbanezMejia15}.

\cite[Ib\'a\~nez-Mej\'{\i}a et al.\ (2015)]{IbanezMejia15}
added the self-gravity of the gas to these models, turning it on after
a full vertical galactic fountain had been set up over the course of
over 200~Myr. This long preparation period also allowed plenty of time
for the dense regions to come into equilibrium with the
supernova-driven turbulence in the surrounding diffuse medium.
However, when the self-gravity was turned on, immediate collapse
ensued, as is shown in the example of Figure~\ref{IbanezMejia15}.  One
immediate and notable result was the
enhancement of the filamentary morphology to the extreme levels
characteristic of observed infrared dark clouds (e.g. \cite[Simon et
al.\ (2006)]{Simon06}).

Within 5 Myr, this collapse lead to the
formation of unresolved dense cores within which rapid star
formation appears inevitable.  As these first models have no mechanism
to simulate the local prompt feedback in ionizing radiation,
stellar winds and so forth from this star formation, we can no longer
follow its evolution after this point.  However it is clear already
that without this local feedback, star formation would reach
unrealistic rates, just as first argued by \cite[Zuckerman \& Palmer
(1974)]{Zuckerman74}.

In order to understand this rapid collapse and fragmentation, it is
necessary to examine the dynamics of these dense
regions. \cite[Ib\'a\~nez-Mej\'{\i}a et al.\ (2015)]{IbanezMejia15}
identified clouds using density contours and demonstrated that these
clouds quite uniformly had velocity dispersions $\sigma <
1$~km~s$^{-1}$, even at radii $R > 30$~pc, quite different from the
\cite[Larson (1981)]{Larson81} size-velocity dispersion relation. This
is shown by the yellow dots in all three panels of
Figure~\ref{IbanezMejia15b}, in contrast to the observed relation
presented as the dashed line, which is a fit to the subset of the
Galactic Ring Survey (\cite[Jackson et al.\ 2006]{Jackson06})
presented in \cite[Heyer et al.\ (2009)]{Heyer09}.  These low velocity
dispersions cannot resist gravitational collapse, and indeed the
clouds do collapse, forming the same sort of filamentary structures
familiar from the free gravitational collapse of the cosmic web.
Although magnetic fields are present in the simulation, they are
unable to effectively resist this collapse.  As they collapse, their
observed velocity dispersions increase to the observed values.

\cite[Ballesteros-Paredes et al.\ (2011)]{BallesterosParedes11}
emphasize that the relationship between surface density, velocity
dispersion, and radius expected from
freely falling material differs by only a factor of $\sqrt{2}$ from
that of material in virial equilibrium, and both appear consistent
with the observations as summarized by \cite[Heyer et
al. (2009)]{Heyer09}.

It is worth noting that in a periodic box simulation of
supernova-driven turbulence, neglecting stratification, but run at
higher resolution, \cite[Padoan et al.\ (2015)]{Padoan15} find
velocity dispersions after 6~Myr that agree with our results.  However
Padoan (priv.\ comm., 2015) finds similar velocity dispersions in
their model even prior to the onset of self-gravity, in contradiction
to our results.  Further investigation will be needed to establish if
the difference is due to their higher resolution or their neglect of
stratification and confinement within the periodic box.

\begin{figure}[bth]
\begin{center}
\includegraphics[width=3.5in,angle=90]{MacLow_IbanezMejia15.eps}
\caption{\label{IbanezMejia15} Vertical projection through a
  sub-region of a supernova-driven, stratified, magnetized model of
  the interstellar medium at times of 0 and 3 Myr after self-gravity
  has been turned on. (The full domain is a $1\times 1\times 20$~kpc
  vertical cut through the galactic disk.)  Prompt collapse occurs,
  forming filaments strongly resembling observed IRDCs, with
  star-forming cores strung along them. Adapted from
  \cite[Ib\'a\~nez-Mej\'{\i}a et al.\ (2015)]{IbanezMejia15}.}
\end{center}
\end{figure}

\begin{figure}[bth]
\begin{center}
\includegraphics[width=1.45in,angle=90]{MacLow_IbanezMejia15b.eps}
\caption{\label{IbanezMejia15b} Radius-velocity dispersion relation
  for clouds drawn from the models prior to the action of
  self-gravity (yellow dots in all three panels), and for clouds at
  times of 2 Myr (left), 4 Myr (center) and 6 Myr (right) after
  self-gravity begins to act.  For comparison, the observed relation
  derived from the data presented by \cite[Heyer et al.\
  (2009)]{Heyer09} is shown by the dashed line. The shaded region
  shows clouds under 10 zones whose velocity dispersion is likely to be strongly suppressed by
  numerical viscosity.  Adapted from
  \cite[Ib\'a\~nez-Mej\'{\i}a et al.\ (2015)]{IbanezMejia15}.}
\end{center}
\end{figure}

\section{Summary}

In this review, I have addressed a few current topics in the theory of
the atomic and molecular interstellar medium.  I have focused on the
following questions, with my current best guesses as to the answers.
\begin{itemize}
  \item Do cold streams contribute to gas accretion onto
  galaxies only at early times and for small galaxies, or far more
  generally as was originally claimed?  The former appears to be the
  result of the numerical simulations most likely to be reliable.
  \item How small is the actual filling factor of dense, molecular gas in
    an interstellar medium when it is observationally determined to be
    completely molecular?  Quite small, with the rest of the volume
    occupied by low-density, hot gas with temperatures exceeding $10^6$~K.
  \item Does molecular hydrogen drive star formation, or rather does
    star formation produce molecular hydrogen and other high density
    tracers as a by-product?  Molecular hydrogen and other molecules
    appear to function far more as a tracer than a cause of star
    formation, particularly since atomic line cooling is nearly as
    effective as molecular line cooling.
  \item Are the observed velocity dispersions in molecular clouds
    caused by supernova-driven turbulence, or by gravitational
    collapse?  Gravitational collapse seems more likely, but this
    question does not yet appear to be settled.
\end{itemize}