Suppression of Star Formation in NGC 1266 Part 1

This paper has broken up into 2 parts. Here is part 1. Part 2 is here.

Fundamentally, understanding how galaxies are born, live, and die comes down to understanding the stars sitting inside the galaxy (well, and the dark matter, but that is “beyond the scope of this blog”). That is why for the first half of the blog of this paper, the focus is on how and what we can learn about the stars in a galaxy when we can’t individually resolve stars. Luckily for us, we are sitting directly inside a very good laboratory of the nature of star formation (the Milky Way!) So if we can’t literally count the young stars, how do we calculate a star formation rate?

A schematic of a star-forming region in the Milky Way. Note: the Brick (image to the left) is a cloud with strange star-forming properties, and the Pleiades (stars in the right-hand image) are definitely not a part of the Brick! This image is thus meant as a schematic to explain why we use certain tracers of star formation.

The figure above is a schematic of a star forming region (note: the Brick and the Pleiades are not actually associated). The stars form out of molecular clouds, because it is a region that is both dense and cold. In the region, the youngest and the brightest stars dominate the light. Those stars emit very strongly in the ultraviolet. If we could see all the UV photons in a star-forming region, we could get a good idea of the number of stars that are forming. Alas, because these stars are still in the dense, clumpy “stellar nurseries”, most of those UV photons never make it out. One way that scientists deal with this is make arguments about the geometry of the nursery clouds and count the “escaped” UV photons. At the same time, those UV photons are getting absorbed by the surrounding gas, then re-emitted by the gas that sits at the edge of the star-forming region. Because most of that gas is hydrogen, the emission comes out as hydrogen line emission. This threshold is also where the term H II region comes from. Finally, because these stellar nurseries sit inside of dense clouds, there is also a lot of dust. Unlike the H II region, which emits specific wavelengths of light, the dust doesn’t discriminate the photons it absorbs, and it will re-emit the star-forming light as a blackbody, peaked at the temperature of the dust (up to ~100 Kelvin in the absolute highest star forming objects, like Arp 220). This is why most astronomers use some tracer of the dust emission (like the infrared blackbody emission), the UV emission (directly tracing the young stars, though probably buried under a lot of absorbing material), and the hydrogen lines to infer star formations. Because the dominant energy mechanism in most star-forming galaxies are the young stars, in normal circumstances, these tracers do a pretty good job.

Getting the molecular gas mass requires some similar assumptions. H2 (molecular hydrogen) is the most abundant gas in the star forming nurseries, but is very hard to detect. The reason for this is that H2 lacks a detectable dipole moment, since the molecule is symmetric. It makes direct measurement of H2 in cold clouds really really difficult. Luckily for us, carbon monoxide (CO) has a substantial dipole moment, so when it rotates, it radiates. And that radiation is bright. Carbon and oxygen are also some of the first metals that are formed when stars “turn on”, so it enriches primordial gas fairly quickly. CO requires nearly the same density conditions as H2, so the two molecules are very often in the same region. It was also found within the Milky Way that the luminosity of CO did a very good job tracing the mass of H2 underlying it in giant molecular clouds (you can read about this more from this beautiful astrobite by Adele Plunkett). So, for the most part, in normal galaxies, CO emission can trace the underlying gas mass, like you can use Hα, UV or dust emission to trace the star formation.

Unless you are NGC 1266.

The CS emission and far-IR emission (tracing in part SF)
in NGC 1266 (on top of HST B-band image)

NGC 1266 breaks all of these “reliable rules”. How? Let’s start with the gas (because the star formation is more complicated). The dependence of H2 mass on CO luminosity requires a big assumption: that the gas in your system looks like the gas in the Milky Way. For NGC 1266, it is probable that that is not true. The gas is not distributed in discrete giant molecular clouds (where the entire relation is derived). So, if we can’t use CO, what can we use? Denser gas and gravity. CS (carbon monosulfide) is much less abundant than CO, and require denser conditions to form. This means that you are probably counting the photons that are coming from all the CS molecules, instead of just the ones at the optically thick surface (as is the case with CO). In this paper, we looked at the CS and used that to figure out the H2 mass, assuming a specific abundance of CS to H2, and it gave us a surface density 10^4 M⊙/pc^2, which is close to what we see in mergers, the location of some of the densest molecular gas in the local universe. The other thing we can do is use force balance, measuring the line width in NGC 1266, which traces rotation and dispersion, and calculate how much mass needs to be in the center to balance the forces. This also told us that the surface density was 10^4 M⊙/pc^2. Excellent, we now had a gas mass. We also detected HCN and CS (dense gas tracers) had broad line wings, telling us that this gas was taking part in the outflow. This changed our estimate of the mass outflow rate, from 13 M⊙/yr to 110 M⊙/yr. In a later blog on my Markarian 231 paper which expands upon why the outflow rate is not a depletion rate.

The spectral energy index of NGC 1266, including the different fits to get star formation.

Star formation was more complicated. Many of the simple assumptions that were made on star formation tracers above could not be made for NGC 1266. The Hα emission in NGC 1266 was excited by shocks, the UV emission was partially from older (rather than newborn) stars, and the absorbing column is so high that the UV is impossible to use, and even the dust blackbody could be hiding a bright supermassive black hole. This meant that most of the tracers would give useful limits, but poor estimates. We had to search deeper than the normal star formation tracers, and settled on using a combination of the radio free-free emission (coming from the same H II regions as the Hα, but without the added shock contamination or obscuration) to calculate the star formation rate. This free-free estimate, along with the size of the emitting region were from ALMA observations. Once we had both of these values (star formation rate and gas mass), we looked at how they related to one another. 

The Schmidt-Kennicutt relation with NGC 1266
overlaid. It sits far away from the relation and far
outside the scatter.

There is a reasonable empirical relationship between the density of star formation and the density of molecular gas, known as the Schmidt-Kennicutt relation. This relation works for star-forming normal galaxies, early-type galaxies, and star bursting galaxies. But it does not work for NGC 1266. In fact, the molecular gas in NGC 1266 is under-forming stars by over a factor of 50. This sort of outlier really had not been seen before NGC 1266, and was the main observing result of this paper. 

 

To be continued...

Part 2 of this paper will discuss not only our guess as to what causes this large suppression of star formation in NGC 1266, but also the implications for not only NGC 1266, but this phenomenon in general.

The official published version can be found on NASA ADS.

To get a PDF version made by me, you can download it here.