Haojing Yan

aut viam inveniam aut faciam

Dusty Galaxies

ULIRGs at z > 6

The detection of CO gas in the SDSS quasar SDSS J114816.64+525150.3 at z = 6.42 (Walter et al. 2003) was a very important discovery because it showed that metals were already abundant as early as when the universe had not yet finished the hydrogen reionization. Metal can only be produced by stars. This means that the universe must have already begun active star formation much earlier than z = 6.42 and been producing/dispersing dust effectively so that this quasar host galaxy could have such a huge amount of CO gas. In 2017, the timeline was moved further back by the detection of [CII] line (λ=158μm) in the z = 7.54 quasar J1342+0928 (Venemans et al. 2017). At this redshift, the age of the universe was only 700 Myr.

Along a parallel path, metals also manifests themselves in the epoch of reionization through ULIRGs at z > 6.3. ULIRGs are dusty starbursts. To form dust, metals are necessary. Riechers et al. (2013) reported an ULIRG at z = 6.34, which was selected using the Herschel SPIRE data in the (Spitzer's) First Look Survey (FLS) field as a "500 μm riser". Strandet et al. (2017) discovered a mm-galaxy, SPT-S J031132-5823.4, at z = 6.900. Even after taking the possible gravitational lensing effect into account, both objects still have SFR > 100 Msolar/yr as inferred from their rest-frame IR luminosity. In other words, these are extreme starbursts completely embedded by dust that was produced by the previous generation of stars; and these previous generation of stars were possibly also formed in starbursts so that there could be enough metals to form dust so quickly.

In this context, finding ULIRGs at z > 6 is extremely interesting to me because they would throw light on the formation of first stars and first galaxies. In Yan et al. (2020, ApJS, 249, 1), we presented a systematic search for very high-z ULIRGs using the 106.5 deg2 HerMES data and the 2.98 deg2 S2CLS data falling in HerMES. The HerMES data alone resulted in 629 500 μm risers, while the S2CLS+HerMES data resulted in 95 "SPIRE-dropouts". For various reasons, neither method could provide a clean selection of very high-z ULIRGs, and those genuine z > 6 objects are likely only a small fraction among the candidates. As compared to other similar selections in the literature, we moved one step further by using deep radio data to flag those that are the most promising candidates at z > 6.3. This is based on the fact that high-z ULIRGs should be radio-weak. Follow-up analysis on our sample is underway.

Yan, Ma, Huang, & Fan (2020, ApJS, 249, 1)

(U)LIRGs : Supernova Fireworks

(U)LIRGs must have high rates of supernova because they have high SFR. A number of repeated IR and radio observations on a few local (U)LIRGs have shown that they indeed have high occurrence of supernovae and/or supernova remnants. However, this is yet to be demonstrated beyond the local universe. The difficulty is two-fold. First, one would need to do this in the rest-frame IR or radio to circumvent the severe dust extinction. Second, one would need high resolution to discern the SNe. Both become exceedingly difficult beyond the local universe. In Yan et al. (2018, ApJ, 867, 21), we presented indirect evidence by using the data in the Spitzer IRAC Dark Field (IDF), whose 3.6/4.5 μm observations spanned 14+ years on a nearly bi-weekly basis. Herschel/SPIRE also observed this area (as it is the region that has the lowest IR background), so we have (U)LIRGs by SPIRE detections. IRAC did not have the spatial resolution to discern SNe individually; however, the outburst of SN would change the brightness of the host galaxy. Therefore, by searching for 3.6/4.5 μm variability of the (U)LIRG hosts, the high rate of SNe could be revealed. We indeed found two very promising cases, whose IRAC Ch1/2 light curves are shown below in the left panel (the middle two objects). Their "net" variations after subtracting the host light are shown in the right panel, which are consistent with the summation of multiple SNe overlapping in time.

A JWST monitoring program will be able to resolve such SNe in (U)LIRGs, which will open up a new window of studying high-z SNe. The IDF is the best field for this purpose because it is in the JWST's CVZ and thus is the easiest to schedule observations.

Using a series of toy model of SN light curves, the features ("bumps") in the net variation can be reasonably explained.

(Yan, Ma, Beacon, & Runge 2018, ApJ, 867, 21)

"Red" but not "Dead" : W4BCGs

Elliptical galaxies are thought to be dominated by old stellar populations ("red") and have long ceased their star formation ("dead"). But it's known that there are counter-examples. Some of the most interesting counter-examples are found among the brightest cluster galaxies (BCGs). The vast majority of the BCG family are indeed red-and-dead, but a very small minority are not: they reside in the so-called "cool-core" (CC) clusters, and they can have SFRs from a few to ~ 100 Msolar/yr fueled by "cooling flows".

While exploring the WISE data archive, we have found an interesting population, which we dubbed the "W4BCGs". Runge & Yan (2018, ApJ, 853, 47; which was James' thesis paper) presented the analysis on a sample of 389 such objects (the largest sample of its kind): they are brightest cluster galaxies (BCGs) selected from the SDSS data, so in theory they should be some of the best representatives of local elliptical galaxies (i.e., high-mass and and long being "dead"); but these are bright in the WISE W4 band (22 μm) and not like anything passively evolving. They are IR-luminous, with the median IR luminosity well in the LIRG regime; in fact, > 27% of them qualify as ULIRGs. While some of them are AGNs, the majority are not (as shown the BPT diagnostics) and they should owe their strong IR emissions to dust-enshrouded star formation (the latter is confirmed by a small sub-sample that has Herschel SPIRE data). Ten of the W4BCGs have Chandra X-ray data, and it seems that seven of them might be in CC clusters. However, we show that in most cases (5 out of these 7) the mass deposition rates from cooling flows cannot account for their observed SFRs. This casts doubt on the theory that cooling flows are the cause of the star formation in non-quiescent BCGs.

Runge & Yan (2018, ApJ, 853, 47)

ULIRG Quasars : starburst and AGN have nothing to do with each other

"AGN feedback" is an often-used mechanism in theoretical works to turn off star formation in galaxies so that they don't grow too big to be matched with observations. On the other hand, there have also been suggestions that this "feedback" could be "positive" so that the AGN could somehow trigger/enhance star formation. Quasars detected in sub-mm/mm are used as examples for the latter. But really?

In Ma & Yan (2015, ApJ, 811, 58), we assembled a large sample of quasars that have strong far-IR emissions and studied the co-evolution of starburst and AGN. This sample has 354 quasars spanning 0.14 ≤ z ≤ 4.7, and was derived based on the SDSS quasar sample and the far-IR data from all the Herschel wide-field surveys. This paper, while being only part of Zhiyuan's final thesis, is 25 pages long and contains a lot of interesting analysis and implications. If I were to highlight, two points are the most important in terms of the "big picture": (1) the amplitude of the starburst in the quasar host galaxies (indicated by the IR luminosity) has nothing to do with the AGN strength (in terms of quasar absolute magnitude, black hole mass, or X-ray luminosity), and (2) the fraction of IR-luminous quasars peaks at around 1.5 ≲ z ≲ 2.5, which is different from the distribution of the parent SDSS quasar sample. See figures below. Both suggest that starburst has nothing to do with AGN if we see them co-exist in galaxies.

(IR luminosity of IR-qusars is not correlated with AGN strength)

(Redshift distribution of IR-luminous quasars peaks at 1.5 ≲ z ≲ 2.5, which is different from that of the parent SDSS sample)

Ma & Yan (2015, ApJ, 811, 58)

Dusty starbursting regions in ULIRGs must be of small sizes (r ≲ 2 kpc)

Inspired by one of the preliminary conclusions in Ma & Yan (2015, ApJ, 811, 58), we carried out a further analysis of the LIR-T relation (which was first observed in SMGs). As here we only cared about LIR, we used a much enlarged far-IR quasar sample. The details are given in Yan & Ma (2016, ApJ, 820, L16). Basically, an object following the LIR-T relation is best understood using the equivalent of the Stefan-Boltzmann law for an MBB with an effective radius of Reff. The very existence of the LIR-T relation means that the size of the dusty starbursting regions in an ULIRG must be small (with Reff ≲ 2 kpc). While strictly speaking Reff is not exactly the physical size, it is a good proxy. Outside of the region outlined by the LIR-T relation is the "zone of avoidance" in terms of minimum star formation rate density of the starbursting region ("minSFRD" as we call it).

Yan & Ma (2016, ApJ, 820, L16)