# Harvard Astronomy 201b

## ARTICLE: The mass function of dense molecular cores and the origin of the IMF (2007)

In Journal Club, Journal Club 2013, Uncategorized on March 12, 2013 at 9:30 pm

## The mass function of dense molecular cores and the origin of the IMF

by  Joao Alves Marco Lombardi & Charles Lada
Summary by Doug Ferrer

## Introduction

One of the main goals of researching the ISM is understanding the connection between the number and properties of stars and the properties of the surrounding galaxy. We want to be able to look at a stellar population and deduce what sort of material it came from, and the reverse–predict what sort of stars we should expect to form in some region given what we know about its properties. The basics of this connection have been known for a while (eg. Bok 1977). Stars form from the gravitational collapse of dense cores of molecular clouds. Thus the properties of stars are the properties of these dense cores modulated by the physical processes that happen during this collapse.

One of the key items we would like to be able to derive from this understanding of star formation is the stellar initial mass function (IMF)–the number of stars of a particular mass as a function of mass. Understanding how the IMF varies from region to region and across time would be extremely useful in many areas of astrophysics, from cosmology to star clusters. In this paper, Alves et al. attempt to explain the origin of the IMF and provide evidence for this explanation by examining the molecular cloud complex in the Pipe nebula. We will look at some background on the IMF, then review the methods used by Alves et al. and asses the implications for star formation.

## The IMF and its Origins

The IMF describes the probability density for any given star to have a mass M, or equivalently the number of stars in a given region with a mass of M. Early work done by Salpeter (1955) showed that the IMF for relatively high mass stars ( $M > 1 M_{sun}$) follows a power law with the index $\alpha \approx -1.3 - 1.4$. For lower masses, the current consensus is for a break below $1 M_{sun}$, and another break below $.1 M_{sun}$ , with a peak around $.3 M_{sun}$. A cartoon of this sort of IMF is shown in Fig. 1.  The actual underlying distribution may in fact be log-normal. This is consistent with stars being formed from collapsing over-densities caused by supersonic turbulence within molecular clouds (Krumholz 2011). This is not particularly strong evidence, however, as the log-normal distribution can result from any process that depends on the product of many random variables.

Figure 1: An illustration of the accepted form of the IMF. There is a power law region for less than $1 M_{sun}$, and a broad peak for lower masses. Image from Krumholz (2011) produced using Chabrier (2003)

The important implication of these results is that there is a characteristic mass scale for star formation of around $.2-.3 M_{sun}$ . There are two (obvious) ways to explain this:

1. The efficiency of star formation peaks at certain scales.
2. There is a characteristic scale for the clouds that stars form from

There has been quite a lot of theoretical work examining option 1 (see Krumhoz 2011 for a relatively recent, accessible review). There are many different physical processes at play in star formation–turbulent MHD, chemistry, thermodynamics, radiative transfer, and gravitational collapse.  Many of these processes are not separately well understood, and each is occurring in its complex, and highly non-linear regime.  We are not even close to a complete description of the full problem. Thus it would not be at all surprising if there were some mass scale, or even several mass scales that are singled out by the star formation process, even if none is presently known. There is some approximate analytical work showing that feedback from stellar outflows may provide such a scale (eg. Shu et al. 1987) . More recent work (e.g. Price and Bate 2008) has shown that magnetic fields cause significant effects on the structure of collapsing cloud cores in numerical simulations, and may reduce or enhance fragmentation depending on magnetic field strength and mass scale.

Nevertheless, the authors are skeptical of the idea that star-formation is not a scale-free process. Per Larson (2005), they do not believe that turbulence or magnetic fields are likely to be very important for the smallest, densest scale clouds that starts ultimately form from. Supersonic turbulence is quickly dissipated on these scales, and magnetic fields are dissipated through ambipolar diffusion–the decoupling of neutral molecules from the ionic plasma. Thus Larson argues that thermal support is the most important process in small cores, and the Jeans analysis will be approximately correct.

The authors thus turn to option 2. It is clear that if the dense cores of star forming clouds already follow a distribution like the IMF, then there will be no need for option 1 as an explanation. Unfortunately though, the molecular cloud mass function (Fig. 2) does not at first glance show any breaks at low mass and has too shallow of a power law index .  But what if we look at only the smallest, densest clumps?

Figure 2: The cumulative cloud mass function (Ak is proortional to mass) for several different cloud complexes from Lombardi et al. (2008). While this is not directly comparable to the IMF, the important take away is that there are no breaks at low mass.

## Observations

Observations using dense gas emission tracers like $C^{18} O$ and $HC^{13}O$ produces mass distributions more like the stellar IMF (eg. Tachihara et al. 2002, Onishi 2002) . However, there are many systematic uncertainties in emission based analysis. To deal with these issues, the authors instead probed dense cloud mass using wide field extinction mapping (this work).   An extinction map was constructed of the nearby Pipe nebula using the NICER method of Lombardi & Alves (2001), which we have discussed in class. This produced the extinction map shown in Fig. 3 below.

Figure 3: NICER Extinction map of the Pipe nebula. A few dense regions are visible, but the noisey, variable background makes it difficult to seperate out seperate cores in a consistent way.

## Data Analysis

The NICER map of the pipe nebula reveals a complex filamentary structure with very high dynamic range in column density ($~ 10^{21} - 10^{22.5} cm^{-2}$). It is difficult to assign cores to regions in such a data set in a coherent way (Alves et al. 2007) using standard techniques–how do we determine what is a core and what is substructure to a core? Since it is precisely the number and distribution of cores that we are interested in, we cannot use a biased method to identify the cores. To avoid this, the authors used a technique called  multiscale wavelet decomposition. Since the authors do not give much information on how this works, we will give a brief overview following a description of a similar technique from Portier-Fozzani et al. (2001).

### Wavelet Decomposition

Wavelet analysis is a hybrid of Fourier and coordinate space techniques. A wavelet is a function that  has a characteristic frequency, position and spatial extent, like the one in Fig.  4. Thus if we convolve a wavelet of a given frequency and length with a signal, it will tell us how the spatial frequency of the signal varies with position. This is the type of analysis used to produce a spectrogram in audio visualization.

Figure 4: A wavelet. This is one is the product of a sinusoid with a gausian envelope. This choice is a perfect tradeoff between spatial and frequency resolution, but other wavelets may be ideal for some other resolution goal. Note that $\Delta x \Delta k \ge 1/2$

The authors used this technique to separate out dense structures from the background. They performed the analysis using wavelets with spatial scales close to three different length scales separated by factors of 2 (.08 pc, .15 pc, .3 pc). Then they combined the structures identified at each length scale into trees based on spatial overlap, and called only the top level of these trees separate cores. The resulting identified cores are shown in fig. 5.  While the cores are shown as circles in fig. 5, this method does not assume spherical symmetry, and the actual identified core shapes are filamentary, reflecting the qualitative appearance of the image

Figure 5: Dense cores identified in the Pipe nebula (circles) from Alves et al. (2007). The radii of the circles is proportional to the core radius determined by the wavelet decomposition analysis

## Results

After obtaining what they claim to be a mostly complete sample of cores, the authors calculate the mass distribution for them. This is done by converting dust extinction to a dust column density. This gives a dust mass for each clump, which can then be extrapolated to a total mass by assuming a set dust fraction.  The result of this is shown in fig. 6 below. The core mass function they obtain is qualitatively similar in shape to the IMF, but scaled by a factor of ~4 in mass. The analysis is only qualitative, and no statistics are done or functions fit to the data. The authors claim that this result evinces a universal star formation efficiency of ~30%, and that this is good agreement with that calculated analytically by (Shu 2004) and numerical simulations. This is again only a qualitative similarity, however.

We should also note that the IMF is hypothesized to be a log-normal distribution. This sort of distribution can arise out of any process that depends multiplicatively on many independent random factors. Thus the fact that dense cores have a mass function that is a scaled version of the IMF is not necessarily good evidence that they share a simple causal link, in the same way that two variables both being normally distributed does not mean that they are any way related.

Figure 6: The mass function of dense molecular cores (points) and the IMF (solid grey line). The dotted gray line is the IMF with mass argument scaled by a factor of 4. The authors note the qualitative agreement, but do not perform any detailed analysis.

## Conclusions

Based on this, the authors conclude that there is no present need for a favored mass scale in star formation as an explanation of the IMF.  Everything can be explained by the distribution produced by dense cloud cores as they are collapsing.  There are a few points of caution however. This is a study of only one cluster, and the data are analyzed using an opaque algorithm that is not publicly available. Additionally, the distributions are not compared statistically (such as with KS), so we have only qualitative similarity.  It would be interesting to see these results replicated for a different cluster using a more transparent statistical method

## References

Bok, B. J., Sim, M. E., & Hawarden, T. G. 1977, Nature, 266, 145

Krumholz, M. R. 2011, American Institute of Physics Conference Series, 1386, 9

Lombardi, M., Lada, C. J., & Alves, J. 2008, A&A, 489, 143

Krumholz, M. R., & Tan, J. C. 2007, ApJ, 654, 304

Alves, J., Lombardi, M., & Lada, C. J. 2007, A&A, 462, L17

Lombardi, M., Alves, J., & Lada, C. J. 2006, A&A, 454, 781

Larson, R. B. 2005, MNRAS, 359, 211

Shu, F. H., Li, Z.-Y., & Allen, A. 2004, Star Formation in the Interstellar Medium: In Honor of
David Hollenbach, 323, 37

Onishi, T., Mizuno, A., Kawamura, A., Tachihara, K., & Fukui, Y. 2002, ApJ, 575, 950

Tachihara, K., Onishi, T., Mizuno, A., & Fukui, Y. 2002, A&A, 385, 909

Portier-Fozzani, F., Vandame, B., Bijaoui, A., Maucherat, A. J., & EIT Team 2001, Sol. Phys.,
201, 271

Shu, F. H., Adams, F. C., & Lizano, S. 1987, ARA&A, 25, 23

Salpeter, E. E. 1955, ApJ, 121, 161

## ARTICLE: Dark Nebulae, Globules, and Protostars

In Journal Club, Journal Club 2013 on February 19, 2013 at 10:45 pm

### Dark Nebulae, Globules, and Protostars by Bart Bok (1977)

Summary by George Miller

### Introduction

In Bart Bok’s 1977 paper Dark nebulae, globules, and protostars (Bok 1977), largely based on a lecture given upon acceptance of the Astroomical Society of the Pacific’s Bruce Medal, he presents two fundamentally different pictures of star formation. The first, constituting the majority of the paper’s discussion, occurs in large Bok globules which are compact, rounded and remarkably well-defined regions of high-extinction ranging from 3′ to 20′.  The globules show a strong molecular hydrogen and dust component and relatively little signs of higher neutral HI concentrations than its surroundings. In contrast, Bok briefly examines star formation in the Magellanic Clouds which show a vast amount of neutral atomic hydrogen and a comparatively small amount of cosmic dust. In this review, I will summarize a number of key points made by Bok, as well as provide additional information and modern developments since the paper’s original publishing.

### Large Bok Globules

#### A history of observations

In 1908, Barnard drew attention to “a number of very black, small, sharply defined spots or holes” in observations of the emission nebula Messier 8 (Barnard 1908).  39 years later Bok published extensive observations of 16 “globules” present in M8 as well others in $\eta$ Carinae, Sagittarius, Ophiuchus and elsewhere, making initial estimates of their distance, diameter and extinction (Bok & Reilly 1947). He further claimed that these newly coined “globules” were gravitationally contracting clouds present just prior to star formation, comparing them to an “insect’s cocoon” (Bok 1948). As we will see, this bold prediction was confirmed over 40 years later to be correct. Today there over 250 globules known within roughly 500 pc of our sun and, as Bok claims in his 1977 paper, identifying more distant sources is difficult due to their small angular diameter and large number of foreground stars.  There are currently four chief methods of measuring the column density within Bok Globules: extinction mappings of background stars, mm/sub-mm dust continuum emission, absorption measurements of galactic mid-IR background emission, and mapping molecular tracers.  See Figure 1 for a depiction of the first three of these methods.  At the time Bok published his paper in 1977, only extinction mapping and molecular tracer methods were readily available, thus I will primarily discuss these two.  For a more in depth discussion, see Goodman et. al. 2009 and the subsequent AST201b Journal Club review.

Figure 1.  Three methods of determining column density of starless molecular cores or Bok globules. (a) K-band image of Barnard 68 and plot of the $A_K$ as a function of radius from the core.  This method measures the H–K excess, uses the extinction law to convert into $A_V$, and then correlated to the $H_2$ column density from UV line measurements, parameterized by f. (b) 1.2-mm dust continuum emission map and ﬂux versus radius for L1544.  $\kappa_{\nu}$ is the dust opacity per unit gas mass, ρ is the dust density, and m the hydrogen mass (corrected for He). (c) 7-μm ISOCAM image and opacity versus radius for ρ Oph D.  In this method the absorbing opacity is related to the hydrogen column via the dust absorption cross section, $\sigma_{\lambda}$.  Figure taken from Bergin & Tafalla 2007.

#### Measuring photometric extinction

Measuring the photometric absorption, and thus yielding a minimum dust mass, for these globules is itself an arduous process. For globules with $A_v<10$   mag, optical observations with large telescopes can be used to penetrate through the globules and observe the background stars.  Here $A_{\lambda} \equiv m_{\lambda}-m_{\lambda, 0} = 2.5 \, log(\frac{F_{\lambda,0}}{F_{\lambda}})$.  Thus an extinction value of $A_v=10$ mag means the flux is decreased by a factor of $10^4$.  By using proper statistics of the typical star counts and magnitudes seen within a nearby unobstructed field of view, extinction measurements can be made for various regions.  It is important to note that the smaller an area one tries to measure an extinction of, the greater the statistical error (due to a smaller number of background stars).  This is one of the key limitations of extinction mappings.  For the denser cores or more opaque globules with $10 < A_V < 20$ mag, observations in the near infrared are needed (which is relatively simple by today’s standards but not so during Bok’s time). This is further complicated due to imprecisely defined BVRI photometric standard sequences for fainter stars, a problem still present today with various highly-sensitive space telescopes such as the HST. Bok mentions two methods. In the past a Racine (or Pickering) prism was used to produce fainter companion images of known standards, yet as discussed by Christian & Racine 1983 this method can produce important systematic errors. The second, and more widely used, method is to pick an easily accessible progression of faint stars and calibrate all subsequent photographic plates (or ccd images) from this. See Saha et. al. 2005 for a discussion of this problem in regards to the Hubble Space Telescope.

Obtaining an accurate photometric extinction for various regions within the globule is valuable as it leads an estimate of the dust density. Bok reports here from his previous Nature paper (Bok et. al. 1977) that the extinction $A_v$ within the Coalsack Globule 2 varies inversely as the square of distance, thus also implying the dust density varies inversely as the cube of distance from the core.  Modern extinction mappings, as seen in Figure 1(a) of Barnard 68,  show that at as one approaches the central core the extinction vs. distance relation actually flattens out nearly to $r^{-1}$.  This result was a key discovery, for the Bonnor-Ebert (BE) isothermal sphere model predicts a softer power law at small radii.  In his paper, Bok remarks “The smooth density gradient seems to show that Globule 2 is […] an object that reminds one of the polytropic models of stars studied at the turn of the century by Lane (1870) and Emden (1907)”.  It is truly incredible how accurate this assessment was.  The Bonnor-Ebert sphere is a model derived from the Lane-Emden equation for an isothermal, self-gravitating sphere which remains in hydrostatic equilibrium.  Figure 2 displays a modern extinction mapping of Barnard 68 along with the corresponding BE sphere model, showing that the two agree remarkably well.  There are, however, a number of detractors from the BE model applied to Bok globules.  The most obvious is that globules are rarely spherical, implying that some other non-symmetric pressure must be present.  Furthermore, the density gradient between a globule’s core and outer regions often exceeds 14 ($\xi_{max} > 6.5$) as required for a stable BE sphere (Alves, Lada & Lada 2001).

Figure 2.  Dust extinction mapping for Barnard 68 plotted against an isothermal Bonnor-Ebert sphere model.  Figure taken from Alves, Lada & Lada 2001.

#### Using CO as a tracer

Important tracer molecules, such as CO, are used to study the abundance of $H_2$, temperatures and kinematics of these globules. Because the more common $^{12}CO$ isotope tends to become optically thick and saturate in regions of higher column density such as globules, the strength of $^{13}CO$ emission is usually used to indicate the density of H2.  The conversion factor of $N_{H_2} = 5.0 \pm 2.5 \times 10^5 \times N_{13},$ from Dickman 1978 has changed little in over three decades. The column density of $H_2$, combined with its known mass and radius of the globule, can then be used to estimate the globule’s total mass. Furthermore, the correlation of $^{13}CO$ density with photometric extinction, $A_v = 3.7 \times 10^{-16} \times N_{13},$ is another strong indication that $^{13}CO$ emission is an accurate tracer for H$_2$ and dust. Further studies using $C^{17}O$ and $C^{18}O$ have also been used to trace even higher densities when even $^{13}CO$ can become optically thick(Frerking et. al. 1982).  As an example, Figure 3 shows molecular lines from the central region of the high-mass star forming region G24.78+0.08.  In the upper panel we can see the difference between the optically thick $^{12}CO$ and thin $C^{18}O$.  The $^{12}CO$ line shows obvious self-absorption peaks associated with an optically thick regime, and one clearly can not make a Gaussian fit to determine the line intensity.   $^{12}CO$, due to the small dipole moment of its $J=1 \rightarrow 0$ transition and thus ability to thermalize at relatively low densities, is also used to measure the gas temperature within globules. These temperatures usually range from 7K to 15K. Finally, the width of CO lines are used to measure the velocity dispersion within the globule. As Bok states, most velocities range from 0.8 to 1.2 km/s. This motion is often complex and measured excess line-widths beyond their thermal values are usually attributed to turbulence (Bergin & Tafalla 2007). Importantly, the line-width vs. size relationship within molecular clouds first discovered by Barnard 1981 does not extend to their denser cores (which have similar velocity motions as Bok globules).  Instead, a “coherence” radius is seen where the non-thermal component of a linewidth is approximately constant (Goodman et. al. 1998).  In the end, as Bok surmises, the subsonic nature of this turbulence implies it plays a small role compared to thermal motions.

Figure 3.  Spectra taken from the core of the high-mass star forming region G24.78+0.08.  The solid line corresponds to $^{12}CO (1\rightarrow 0)$, $^{12}CO (2\rightarrow 1)$, and $C^{32}S (3\rightarrow 2)$, the dashed line to $^{13}CO (1\rightarrow 0)$$^{13}CO (2\rightarrow 1)$, and $C^{34}S (3\rightarrow 2)$ and the dotted line to $C^{18}O (1\rightarrow 0)$.  From the top panel, one can clearly see the difference between the optically thick, saturated $^{12}CO (1\rightarrow 0)$ line and the optically thin $C^{18}O (1\rightarrow 0)$ transition.  Figure taken from Cesaroni et. al. 2003.

#### The current status of Bok globules

Today, the majority of stars are thought to originate within giant molecular clouds or larger dark cloud complexes, with only a few percent coming from Bok globules. However, the relative simplicity of these globules still make them important objects for studying star formation. While an intense debate rages today regarding the influence of turbulence, magnetic fields, and other factors on star formation in GMCs, these factors are far less important than simple gravitational contraction within Bok globules. The first list of candidate protostars within Bok globules, obtained by co-adding IRAS images, was published in 1990 with the apropos title “Star formation in small globules – Bart Bok was correct” (Yun & Clemens 1990).  To conduct the search, Yun & Clemens first fit a single-temperature modified blackbody model the the IRAS 60 and 100 μm images (after filtering out uncorrelated background emission) to obtain dust temperature and optical depth values.  This result was then used as a map to search for spatially correlated 12 and 25 μm point sources (see Figure 4.).  More evidence of protostar outflows (Yun & Clemens 1992), Herbig-Haro objects due to young-star jets (Reipurth et al. 1992) and the initial stages of protostar collapse (Zhou et. al. 1993) have also been detected within Bok Globules. Over 60 years after Bok’s pronouncement that these globules were “insect cocoons” encompassing the final stages of protostar formation, his hypothesis remains remarkably accurate and validated. It is truly “pleasant indeed that globules are there for all to observe!”

Figure 4.    (a) Contour map of the dust temperature $T_{60/100}$ of the Bok Globule CB60 derived from 60 and 100 μm IRAS images.  (b) 12 μm IRAS image of CB60 after subtracting background emission using median-filtering.  This source is thought to be a young stellar object or protostar located within the globule.  The other 12 μm field sources seen in (b) are thought not to be associated with the globule. Figure taken from Yun & Clemens 1990.

### Magellanic Cloud Star Formation

At the end of his paper, Bok makes a 180 degree turn and discusses the presence of young stars and blue globulars within the Magellanic Clouds. These star formation regions stand in stark contrast to the previously discussed Bok globules; they contain a rich amount of HI and comparatively small traces of dust, they are far larger and more massive, and they form large clusters of stars as opposed to more isolated systems. Much more is known of the star-formation history in the MCs since Bok published this 1977 paper. The youngest star populations in the MCs are found in giant and supergiant shell structures which form filamentary structures throughout the cloud. These shells are thought to form from supernova, ionizing radiation and stellar wind from massive stars which is then swept into the cool, ambient molecular clouds. Further gravitational, thermal and fluid instabilities fragment and coalesce these shells into denser star-forming regions and lead to shell-shell interactions (Dawson et. al. 2013). The initial onset of this new ($\sim$ 125 Myr) star formation is thought to be due to close encounters between the MCs, and is confirmed by large-scale kinematic models (Glatt et al. 2010).

### References

Alves, J. F., Lada, C. J., & Lada, E. A. 2001, Nature, 409, 159
Barnard, E. E. 1908, Astronomische Nachrichten, 177, 231
Bergin, E. A., & Tafalla, M. 2007, ARA&A, 45, 339
Bok, B. J. 1948, Harvard Observatory Monographs, 7, 53
—. 1977, PASP, 89, 597
Bok, B. J., & Reilly, E. F. 1947, ApJ, 105, 255
Bok, B. J., Sim, M. E., & Hawarden, T. G. 1977, Nature, 266, 145
Cesaroni, R., Codella, C., Furuya, R. S., & Testi, L. 2003, A&A, 401, 227
Christian, C. A., & Racine, R. 1983, PASP, 95, 457
Dawson, J. R., McClure-Griffiths, N. M., Wong, T., et al. 2013, ApJ, 763, 56
Dickman, R. L. 1978, ApJS, 37, 407
Frerking, M. A., Langer, W. D., & Wilson, R. W. 1982, ApJ, 262, 590
Glatt, K., Grebel, E. K., & Koch, A. 2010, A&A, 517, A50
Goodman, A. A., Barranco, J. A., Wilner, D. J., & Heyer, M. H. 1998, ApJ, 504, 223
Goodman, A. A., Pineda, J. E., & Schnee, S. L. 2009, ApJ, 692, 91
Reipurth, B., Heathcote, S., & Vrba, F. 1992, A&A, 256, 225
Saha, A., Dolphin, A. E., Thim, F., & Whitmore, B. 2005, PASP, 117, 37
Yun, J. L., & Clemens, D. P. 1990, ApJ, 365, L73
—. 1992, ApJ, 385, L21
Zhou, S., Evans, II, N. J., Koempe, C., & Walmsley, C. M. 1993, ApJ, 404, 232

## ARTICLE: On the Dark Markings in the Sky

In Journal Club, Journal Club 2013 on February 8, 2013 at 2:46 pm

On the Dark Markings in the Sky by Edward E. Barnard (1919)

Summary by Hope Chen

#### Abstract

By examining photographic plates of various regions on the sky, Edward E. Barnard concluded in this paper that what he called “dark markings” were in fact due to the obscuration of nearby nebulae in most cases. This result had a significant impact on the debate regarding the size and the dimension of the Milky Way and also the research of the interstellar medium, particularly work by Vesto Slipher, Heber Curtis and Robert Trumpler. The publication of  Photographic Atlas of Selected Regions of the Milky Way after Barnard’s death, which included many of the regions mentioned in the paper, further provided a new method of doing astronomy research. In this paper and the Atlas, we are also able to see a paradigm very different from that of today.

It is now well-known that the interstellar medium causes extinction of light from background stars. However, think of a time when the infrared imaging was impossible, and the word “photon” meant nothing but a suspicious idea. Back in such a time in the second decade of the twentieth century, Edward Edison Barnard, by looking at hundreds of photographic plates, proposed an insightful idea that “starless” patches of the sky were dark because they are obscured by nearby nebulae. This idea not only built the foundation of the modern concept of the interstellar medium, but also helped astronomers figure out that the Universe extended so much farther beyond the Milky Way.

#### Young Astronomer and His Obsession of the Sky

In 1919, E. E. Barnard published this paper and raised the idea that what he called “dark markings” are mostly obscuration from nebulae close to us. The journey, however, started long before the publication of this paper. Born in Nashville, Tennessee in 1857, Barnard was not able to receive much formal education owing to poverty. His first interest, which became important for his later career, was in photography. He started working as a photographer’s assistant at the age of nine, and the work continued throughout most of his teenage years. He then developed an interest in astronomy, or rather, “star-gazing,” and would go watch the sky almost every night with his own telescope. He took courses in natural sciences at Vanderbilt University and started his professional career as an astronomer at the Lick Observatory in 1888. He helped build the Bruce Photographic Telescope at the Lick Observatory and there he started taking pictures of the sky on photographic plates. He then moved on to his career at the Yerkes Observatory at Chicago University and worked there until his death in 1922. (Introduction of the Atlas, Ref. 2)

Fig. 1 One of the many plates in the Atlas including the region around Rho Ophiuchii, which was constantly mentioned in many of Barnard’s works. (Ref. 2)

Fig. 1 is one of the many plates taken at the Yerkes Observatory. It shows the region near Rho Ophiuchii, which was a region constantly and repetitively visited by Barnard and his telescope. Barnard noted in his description of this plate, “the [luminous] nebula itself is a beautiful object. With its outlying connections and the dark spot in which it is placed and the vacant lanes running to the East from it, … it gives every evidence that it obscures the stars beyond it.” Numerous similar comments spread throughout his descriptions of various regions covered in A Photographic Atlas of Selected Regions of the Milky Way (hereafter, the Atlas). Then finally in his 1919 paper, he concluded, “To me these are all conclusive evidence that masses of obscuring matter exist in space and are readily shown on photographs with the ordinary portrait lenses,” although “what the nature of this matter may be is quite another thing.” The publication of these plates in the Atlas (unfortunately after his death, put together by Miss Mary R. Calvert, who was Barnard’s assistant at the Yerkes Observatory and helped publish many of Barnard’s works after his death) also provided a new way of conducting astronomical research just as the World Wide Telescope does today. The Atlas for the first time allowed researchers to examine the image and the astronomical coordinates along with lists of dominant objects at the same time.

Except quoting Vesto Slipher’s work on spectrometry measurements of these nebulae, most of the evidences in Barnard’s paper seemed rather qualitative than quantitative. So, as of today’s standard, was the “evidence” really conclusive? Again, the question cannot be answered without knowing the limits of astronomical research at the time. Besides an immature understanding of the underlying physics, astronomers in the beginning of the twentieth century were limited by the lack of tools on both the observation and analysis fronts. Photographic plates as those in the Atlas were pretty much the most advanced imaging technique at the time, on which even a quantitative description of “brightness” was not easy, not to mention an estimation of the extinction of these “dark markings.” However, this being said, a very meaningful and somewhat “quantitative” assumption was drawn in Barnard’s paper: the field stars were more or less uniformly distributed. Barnard came to this assumption by looking at many different places, both in the galactic plane and off the plane, and observing the densities of field stars in these regions. Although numbers were not given in the paper, this was inherently similar to a star count study. Eventually, this assumption lead to what Barnard thought as the conclusive evidence of these dark markings being obscuring nebulae instead of “vacancies.” Considering the many technical limits at the time, while the paper might not seem to be scientific in today’s standard, this paper did pose a “conclusion” which was strong enough to sustain many of the more quantitative following examinations.

#### The “Great Debate”

Almost at the same time, perviously mentioned Vesto Slipher (working at the Lowell Observatory) began taking spectroscopic measurements of various clouds and tried to understand the constituents of these clouds. Although limited by the wavelength range and the knowledge of different radiative processes (the molecular transition line emission used largely in the research of the interstellar medium today was not observed until half a century later in 1970, by Robert Wilson, who, on a side note, also discovered the Cosmic Microwave Background), Slipher was able to determine the velocities of clusters by measuring the Doppler shifts and concluded that many of these clusters move at a faster rate than the escape velocity of the Milky Way (Fig. 2). This result, coupled with Barnard’s view of intervening nebulae, revolutionized the notion of the Universe in the 1920s.

Fig. 2 The velocity measurements from spectroscopic observations done by Vesto Slipher. (Ref. 3)

On April 26, 1920 (and in much of the 1920s), the “Great Debate” took place between Harlow Shapley (the Director of Harvard College Observatory at the time) and Curtis Heber (the Lick Observatory, 1902 to 1920). The general debate concerned the dimension of the Universe and the Milky Way, but the basic issue was simply whether distant “spiral nebulae” were small and lay within the Milky Way or whether they were large and independent galaxies. Besides the distance and the velocity measurements, which suffered from large uncertainties due to the technique available at the time, Curtis Heber was able to “win” the debate by claiming that dark lanes in the “Great Andromeda Nebula” resemble local dark clouds as those observed by Barnard (Fig. 3, taken in 1899). The result of the debate then sparked a large amount of work on “extragalactic astronomy” in the next two decades and was treated as the beginning of this particular research field.

Fig. 3 The photographic plate of the “Great Andromeda Nebula” taken in 1988 by Isaac Roberts.

#### The Paper Finally Has a Plot

Then after the first three decades of the twentieth century, astronomers were finally equipped with a relatively more correct view of the Universe, the idea of photons and quantum theory. In 1930, Robert J. Trumpler (the namesake of the Trumpler Award) published his paper about reddening and reconfirmed the existence of local “dark nebulae.” Fig. 4 shows the famous plot in his paper which showed discrepancies between diameter distances and photometric distances of clusters. In the same paper, Trumpler also tried to categorize effects of the ISM on light from background stars, including what he called “selective absorption” or reddening as it is known today. This paper, together with many of Trumpler’s other papers, is one of the first systematic research results in understanding the properties of Barnard’s dark nebulae, which are now known under various names such as clouds, clumps, and filaments, in the interstellar medium.

Fig. 4 Trumpler’s measurements of diameter distances v. photometric distances for various clusters.

#### Moral of the Story

As Alyssa said in class, it is often more beneficial than we thought to understand what astronomers knew and didn’t know at different periods of time and how we came to know what we see as common sense today, not only in the historically interesting sense but also in the sense of better understanding of various ideas. In this paper, Barnard demonstrated a paradigm which we may call unscientific today but made a huge leap into what later became the modern research field of the interstellar medium.

## CHAPTER: How do we know there is an ISM?

In Book Chapter on January 29, 2013 at 8:35 pm

(updated for 2013)

Early astronomers pointed to 3 lines of evidence for the ISM:

• Extinction. The ISM obscures the light from background stars. In 1919, Barnard (JC 2011, 2013) called attention to these “dark markings” on the sky, and put forward the (correct) hypothesis that these were the silhouettes of dark clouds. A good rule of thumb for the amount of extinction present is 1 magnitude of extinction per kpc (for typical, mostly unobscured lines-of-sight).
• Reddening. Even when the ISM doesn’t completely block background starlight, it scatters it. Shorter-wavelength light is preferentially scattered, so stars behind obscuring material appear redder than normal. If a star’s true color is known, its observed color can be used to infer the column density of the ISM between us and the star. Robert Trumpler first used measurements of the apparent “cuspiness” and the brighnesses of star clusters in 1930 to argue for the presence of this effect. Reddening of stars of “known” color is the basis of NICER and related techniques used to map extinction today.
• Stationary Lines. Spectral observations of binary stars show doppler-shifted lines corresponding to the radial velocity of each star. In addition, some of these spectra exhibit stationary (i.e. not doppler-shifted) absorption lines due to stationary material between us and the binary system. Johannes Hartmann first noticed this in 1904 when investigating the spectrum of $\delta$ Orionis: “The calcium line at $\lambda 3934$ [angstroms] exhibits a very peculiar behavior. It is distinguished from all the other lines in this spectrum, first by the fact that it always appears extraordinarily week, but almost perfictly sharp… Closer study on this point now led me to the quite surprising result that the calcium line… does not share in the periodic displacements of the lines caused by the orbital motion of the star”

Helpful References: Good discussion of the history of extinction and reddening, from Michael Richmond.