Star-forming Environments Research Articles

Radiation Hydrodynamic Simulations of Dust-driven Winds

Abstract We study dusty winds driven by radiation pressure in the atmosphere of a rapidly star-forming environment. We apply the variable Eddington tensor algorithm to re-examine the two-dimensional radiation hydrodynamic problem of a column of gas that is accelerated by a constant infrared radiation flux. In the absence of gravity, the system is primarily characterized by the initial optical depth of the gas. We perform several runs with different initial optical depths and resolutions. We find that the gas spreads out along the vertical direction, as its mean velocity and velocity dispersion increase. In contrast to previous work using the flux-limited diffusion algorithm, we find little evolution in the trapping factor. The momentum coupling between radiation and gas in the absence of gravity is similar to that with gravity. For Eddington ratio increasing with the height in the system, the momentum transfer from the radiation to the gas is not merely , but amplified by a factor of , where is the integrated infrared optical depth through the system, and , decreasing with the optical depth. We apply our results to the atmosphere of galaxies and conclude that radiation pressure may be an important mechanism for driving winds in the most rapidly star-forming galaxies and starbursts.

The Molecular Gas Environment in the 20 km s−1 Cloud in the Central Molecular Zone

Abstract We recently reported a population of protostellar candidates in the 20 km s−1 cloud in the Central Molecular Zone of the Milky Way, traced by H2O masers in gravitationally bound dense cores. In this paper, we report molecular line studies with high angular resolution (∼3″) of the environment of star formation in this cloud. Maps of various molecular line transitions as well as the continuum at 1.3 mm are obtained using the Submillimeter Array. Five NH3 inversion lines and the 1.3 cm continuum are observed with the Karl G. Jansky Very Large Array. The interferometric observations are complemented with single-dish data. We find that the CH3OH, SO, and HNCO lines, which are usually shock tracers, are better correlated spatially with the compact dust emission from dense cores among the detected lines. These lines also show enhancement in intensities with respect to SiO intensities toward the compact dust emission, suggesting the presence of slow shocks or hot cores in these regions. We find gas temperatures of ≳100 K at 0.1 pc scales based on RADEX modeling of the H2CO and NH3 lines. Although no strong correlations between temperatures and linewidths/H2O maser luminosities are found, in high-angular-resolution maps we note several candidate shock-heated regions offset from any dense cores, as well as signatures of localized heating by protostars in several dense cores. Our findings suggest that at 0.1 pc scales in this cloud star formation and strong turbulence may together affect the chemistry and temperature of the molecular gas.

High Dense Gas Fraction in a Gas-rich Star-forming Galaxy at z = 1.2∗

Abstract We report observations of dense molecular gas in the star-forming galaxy EGS 13004291 (z = 1.197) using the Plateau de Bure Interferometer. We tentatively detect HCN and HNC emission when stacked together at significance, yielding line luminosities of K km s−1 pc2 and K km s−1 pc2, respectively. We also set 3σ upper limits of <7–8 ×109 K km s−1 pc2 on the , ), and HC3N(J = 20 → 19) line luminosities. We serendipitously detect CO emission from two sources at and in the same field of view. We also detect CO( ) emission in EGS 13004291, showing that the excitation in the previously detected CO( ) line is subthermal ( ). We find a line luminosity ratio of / = 0.17 ± 0.07, as an indicator of the dense gas fraction. This is consistent with the median ratio observed in galaxies ( / = 0.16 ± 0.07) and nearby ULIRGs ( / = 0.13 ± 0.03), but higher than that in local spirals ( / = 0.04 ± 0.02). Although EGS 13004291 lies significantly above the galaxy main sequence at , we do not find an elevated star formation efficiency (traced by / ) as in local starbursts, but a value consistent with main-sequence galaxies. The enhanced dense gas fraction, the subthermal gas excitation, and the lower than expected star formation efficiency of the dense molecular gas in EGS 13004291 suggest that different star formation properties may prevail in high-z starbursts. Thus, using / as a simple recipe to measure the star formation efficiency may be insufficient to describe the underlying mechanisms in dense star-forming environments inside the large gas reservoirs.

The star formation environment of the FU Ori type star V582 Aur

We have studied the environment of the FU Ori type star V582 Aur. Our aim is to explore the star-forming region associated with this young eruptive star. Using slitless spectroscopy we searched for H alpha emission stars within a field of 11.5arcmin \times 11.5arcmin, centred on V582 Aur. Based on UKIDSS and Spitzer Space Telescope data we further selected infrared-excess young stellar object candidates. In all, we identified 68 candidate low-mass young stars, 16 of which exhibited H alpha emission in the slitless spectroscopic images. The colour-magnitude diagram of the selected objects, based on IPHAS data, suggests that they are low-mass pre-main-sequence stars associated with the Aur OB 1 association, located at a distance of 1.3 kpc from the Sun. The bright-rimmed globules in the local environment of V582 Aur probably belong to the dark cloud LDN~1516. Our results suggest that star formation in these globules might have been triggered by the radiation field of a few hot members of Aur OB 1. The bolometric luminosity of V582 Aur, based on archival photometric data and on the adopted distance, is 150-320 Lsun.

STAR FORMATION AND SUPERCLUSTER ENVIRONMENT OF 107 NEARBY GALAXY CLUSTERS

ABSTRACT We analyze the relationship between star formation (SF), substructure, and supercluster environment in a sample of 107 nearby galaxy clusters using data from the Sloan Digital Sky Survey. Previous works have investigated the relationships between SF and cluster substructure, and cluster substructure and supercluster environment, but definitive conclusions relating all three of these variables has remained elusive. We find an inverse relationship between cluster SF fraction (f SF) and supercluster environment density, calculated using the Galaxy luminosity density field at a smoothing length of 8 h −1 Mpc (D8). The slope of f SF versus D8 is −0.008 ± 0.002. The f SF of clusters located in low-density large-scale environments, 0.244 ± 0.011, is higher than for clusters located in high-density supercluster cores, 0.202 ± 0.014. We also divide superclusters, according to their morphology, into filament- and spider-type systems. The inverse relationship between cluster f SF and large-scale density is dominated by filament- rather than spider-type superclusters. In high-density cores of superclusters, we find a higher f SF in spider-type superclusters, 0.229 ± 0.016, than in filament-type superclusters, 0.166 ± 0.019. Using principal component analysis, we confirm these results and the direct correlation between cluster substructure and SF. These results indicate that cluster SF is affected by both the dynamical age of the cluster (younger systems exhibit higher amounts of SF); the large-scale density of the supercluster environment (high-density core regions exhibit lower amounts of SF); and supercluster morphology (spider-type superclusters exhibit higher amounts of SF at high densities).

Gemini/GRACES spectroscopy of stars in Tri II

The chemical abundance ratios and radial velocities for two stars in the recently discovered Triangulum II faint dwarf galaxy have been determined from high resolution, medium signal-to-noise ratio spectra from the Gemini-GRACES facility. These stars have stellar parameters and metallicities similar to those derived from their photometry and medium-resolution Ca II triplet spectra, and supports that Triangulum II has a metallicity spread consistent with chemical evolution in a dwarf galaxy. The elemental abundances show that both stars have typical calcium abundances and barium upper limits for their metallicities, but low magnesium and sodium. This chemical composition resembles some stars in dwarf galaxies, attributed to inhomogeneous mixing in a low star formation environment, and/or yields from only a few supernova events. One of our targets (Star40) has an enhancement in potassium, and resembles some stars in the unusual outer halo star cluster, NGC 2419. Our other target (Star46) appears to be a binary based on a change in its radial velocity (Delta v(rad) = 24.5 +/- 2.1 km/s). This is consistent with variations found in binary stars in other dwarf galaxies. While this serves as a reminder of the high binary fraction in these ultra faint dwarf galaxies, this particular object has had little impact on the previous determination of the velocity dispersion in Triangulum II.

A PORTRAIT OF COLD GAS IN GALAXIES AT 60 pc RESOLUTION AND A SIMPLE METHOD TO TEST HYPOTHESES THAT LINK SMALL-SCALE ISM STRUCTURE TO GALAXY-SCALE PROCESSES

ABSTRACT The cloud-scale density, velocity dispersion, and gravitational boundedness of the interstellar medium (ISM) vary within and among galaxies. In turbulent models, these properties play key roles in the ability of gas to form stars. New high-fidelity, high-resolution surveys offer the prospect to measure these quantities across galaxies. We present a simple approach to make such measurements and to test hypotheses that link small-scale gas structure to star formation and galactic environment. Our calculations capture the key physics of the Larson scaling relations, and we show good correspondence between our approach and a traditional “cloud properties” treatment. However, we argue that our method is preferable in many cases because of its simple, reproducible characterization of all emission. Using, low-J 12CO data from recent surveys, we characterize the molecular ISM at 60 pc resolution in the Antennae, the Large Magellanic Cloud (LMC), M31, M33, M51, and M74. We report the distributions of surface density, velocity dispersion, and gravitational boundedness at 60 pc scales and show galaxy-to-galaxy and intragalaxy variations in each. The distribution of flux as a function of surface density appears roughly lognormal with a 1σ width of ∼0.3 dex, though the center of this distribution varies from galaxy to galaxy. The 60 pc resolution line width and molecular gas surface density correlate well, which is a fundamental behavior expected for virialized or free-falling gas. Varying the measurement scale for the LMC and M31, we show that the molecular ISM has higher surface densities, lower line widths, and more self-gravity at smaller scales.

HERSCHEL/HIFI SPECTRAL MAPPING OF C+, CH+, AND CH IN ORION BN/KL: THE PREVAILING ROLE OF ULTRAVIOLET IRRADIATION IN CH+ FORMATION

ABSTRACT The CH+ ion is a key species in the initial steps of interstellar carbon chemistry. Its formation in diverse environments where it is observed is not well understood, however, because the main production pathway is so endothermic (4280 K) that it is unlikely to proceed at the typical temperatures of molecular clouds. We investigate the formation of this highly reactive molecule with the first velocity-resolved spectral mapping of the CH+ J = 1−0, 2−1 rotational transitions, three sets of CH Λ-doubled triplet lines, 12C+ and 13C+ , and CH3OH 835 GHz E-symmetry Q-branch transitions, obtained with Herschel/HIFI over a region of ≈12 arcmin2 centered on the Orion BN/KL source. We present the spatial morphologies and kinematics, cloud boundary conditions, excitation temperatures, column densities, and 12C+ optical depths. Emission from all of C+, CH+, and CH is indicated to arise in the diluted gas, outside the explosive, dense BN/KL outflow. Our models show that UV irradiation provides favorable conditions for steady-state production of CH+ in this environment. Surprisingly, no spatial or kinematic correspondences of the observed species are found with H2 S(1) emission tracing shocked gas in the outflow. We propose that C+ is being consumed by rapid production of CO to explain the lack of both C+ and CH+ in the outflow. Hence, in star-forming environments containing sources of shocks and strong UV radiation, a description of the conditions leading to CH+ formation and excitation is incomplete without including the important—possibly dominant—role of UV irradiation.

Mass estimates for very cold (<8 K) gas in molecular cloud cores

The mass of prestellar cores is an essential ingredient to understand the onset of star formation in the core. The low level of emission from cold dust may keep parts of it hidden from observation. We aim to determine the fraction of core mass in the temperature range < 8 K that can be expected for typical low- and high-mass star formation regions. We calculate the dust temperature within standard spherically symmetric prestellar cores for a grid of density powerlaws in the outer core regions, core masses, and variation in the external multi-component radiation field. The dust is assumed to be composed of amorphous silicate and carbon, and variations of its optical properties are discussed. As measure for the distribution of cores and clumps, we use core mass functions derived for various environments. In view of the high densities in very cold central regions, dust and gas temperatures are assumed to be equal. We find that the fraction of mass with temperatures < 8 K in typical low- and high-mass cores is < 20\%. It is possible to obtain higher fractions of very cold gas by placing intermediate- or high-mass cores in a typical low-mass star formation environment. We show that the mass uncertainty arising from FIR to mm modeling of very cold dust emission is smaller than the mass uncertainty due to the unknown dust opacities.

PANCHROMATIC HUBBLE ANDROMEDA TREASURY. XVI. STAR CLUSTER FORMATION EFFICIENCY AND THE CLUSTERED FRACTION OF YOUNG STARS

ABSTRACT We use the Panchromatic Hubble Andromeda Treasury survey data set to perform spatially resolved measurements of star cluster formation efficiency (Γ), the fraction of stellar mass formed in long-lived star clusters. We use robust star formation history and cluster parameter constraints, obtained through color–magnitude diagram analysis of resolved stellar populations, to study Andromeda’s cluster and field populations over the last ∼300 Myr. We measure Γ of 4%–8% for young, 10–100 Myr-old populations in M31. We find that cluster formation efficiency varies systematically across the M31 disk, consistent with variations in mid-plane pressure. These Γ measurements expand the range of well-studied galactic environments, providing precise constraints in an H i-dominated, low-intensity star formation environment. Spatially resolved results from M31 are broadly consistent with previous trends observed on galaxy-integrated scales, where Γ increases with increasing star formation rate surface density (ΣSFR). However, we can explain observed scatter in the relation and attain better agreement between observations and theoretical models if we account for environmental variations in gas depletion time (τ dep) when modeling Γ, accounting for the qualitative shift in star formation behavior when transitioning from a H2-dominated to a H i-dominated interstellar medium. We also demonstrate that Γ measurements in high ΣSFR starburst systems are well-explained by τ dep-dependent fiducial Γ models.

LATE POP III STAR FORMATION DURING THE EPOCH OF REIONIZATION: RESULTS FROM THE RENAISSANCE SIMULATIONS

ABSTRACT We present results on the formation of Population III (Pop III) stars at redshift 7.6 from the Renaissance Simulations, a suite of extremely high-resolution and physics-rich radiation transport hydrodynamics cosmological adaptive-mesh refinement simulations of high-redshift galaxy formation performed on the Blue Waters supercomputer. In a survey volume of about 220 comoving Mpc3, we found 14 Pop III galaxies with recent star formation. The surprisingly late formation of Pop III stars is possible due to two factors: (i) the metal enrichment process is local and slow, leaving plenty of pristine gas to exist in the vast volume; and (ii) strong Lyman–Werner radiation from vigorous metal-enriched star formation in early galaxies suppresses Pop III formation in (“not so”) small primordial halos with mass less than ∼3 × 107 M ⊙. We quantify the properties of these Pop III galaxies and their Pop III star formation environments. We look for analogs to the recently discovered luminous Ly α emitter CR7, which has been interpreted as a Pop III star cluster within or near a metal-enriched star-forming galaxy. We find and discuss a system similar to this in some respects, however, the Pop III star cluster is far less massive and luminous than CR7 is inferred to be.

Molecular cloud formation in high-shear, magnetized colliding flows

The colliding flows (CF) model is a well-supported mechanism for generating molecular clouds. However, to-date most CF simulations have focused on the formation of clouds in the normal-shock layer between head-on colliding flows. We performed simulations of magnetized colliding flows that instead meet at an oblique-shock layer. Oblique shocks generate shear in the post-shock environment, and this shear creates inhospitable environments for star formation. As the degree of shear increases (i.e. the obliquity of the shock increases), we find that it takes longer for sink particles to form, they form in lower numbers, and they tend to be less massive. With regard to magnetic fields, we find that even a weak field stalls gravitational collapse within forming clouds. Additionally, an initially oblique collision interface tends to reorient over time in the presence of a magnetic field, so that it becomes normal to the oncoming flows. This was demonstrated by our most oblique shock interface, which became fully normal by the end of the simulation.

Low-end mass function of the Quintuplet cluster

The Quintuplet and Arches clusters, which were formed in the harsh environment of the Galactic Center (GC) a few million years ago, have been excellent targets for studying the effects of a star-forming environment on the initial mass function (IMF). In order to estimate the shape of the low-end IMF of the Arches cluster, Shin & Kim devised a novel photometric method that utilizes pixel intensity histograms (PIHs) of the observed images. Here, we apply the PIH method to the Quintuplet cluster and estimate the shape of its low-end IMF below the magnitude of completeness limit as set by conventional photometry. We found that the low-end IMF of the Quintuplet is consistent with that found for the Arches cluster--Kroupa MF, with a significant number of low-mass stars below 1 $M_{\odot}$. We conclude that the most likely IMFs of the Quintuplet and the Arches clusters are not too different from the IMFs found in the Galactic disc. We also find that the observed PIHs and stellar number density profiles of both clusters are best reproduced when the clusters are assumed to be at three-dimensional distances of approximately 100 pc from the GC.

Resolving the chemical substructure of Orion-KL(Corrigendum)

The Kleinmann-Low nebula in Orion (Orion-KL) is the nearest example of a high-mass star-forming environment. For the first time, we complemented 1.3 mm Submillimeter Array (SMA) interferometric line survey with IRAM 30 m single-dish observations of the Orion-KL region. Covering a 4 GHz bandwidth in total, this survey contains over 160 emission lines from 20 species (25 isotopologues), including 11 complex organic molecules (COMs). At a spatial resolution of 1200 AU, the continuum substructures are resolved. Extracting the spectra from individual substructures and providing the intensity-integrated distribution map for each species, we studied the small-scale chemical variations in this region. Our main results are: (1) We identify lines from the low-abundance COMs CH3COCH3 and CH3CH2OH, as well as tentatively detect CH3CHO and long carbon-chains C6H and HC7N. (2) We find that while most COMs are segregated by type, peaking either towards the hot core (e.g., N-bearing species) or the compact ridge (e.g., O-bearing species like HCOOCH3 and CH3OCH3), while the distributions of others do not follow this segregated structure (e.g., CH3CH2OH, CH3OH, CH3COCH3). (3) We find a second velocity component of HNCO, SO2, 34SO2, and SO lines, which may be associated with a strong shock event in the low-velocity outflow. (4) Temperatures and molecular abundances show large gradients between central condensations and the outflow regions, illustrating a transition between hot molecular core and shock-chemistry dominated regimes. Our observations of spatially resolved chemical variations in Orion-KL provide the nearest reference source for hot molecular core and outflow chemistry, which will be an important example for interpreting the chemistry of more distant HMSFRs.

Dust as interstellar catalyst

Context. Interstellar dust particles, which represent 1% of the total mass, are recognized to be very powerful interstellar catalysts in star-forming regions. The presence of dust can have a strong impact on the chemical composition of molecular clouds. While observations show that many species that formed onto dust grains populate the gas phase, the process that transforms solid state into gas phase remains unclear. Aims. The aim of this paper is to consider the chemical desorption process, i.e. the process that releases solid species into the gas phase, in astrochemical models. These models allow determining the chemical composition of star-forming environments with an accurate treatment of the solid-phase chemistry. Methods. In paper I we derived a formula based on experimental studies with which we quantified the efficiencies of the chemical desorption process. Here we extend these results to astrophysical conditions. Results. The simulations of astrophysical environments show that the abundances of gas-phase methanol and H2O2 increase by four orders of magnitude, whereas gas-phase H2CO and HO2 increase by one order of magnitude when the chemical desorption process is taken into account. The composition of the ices strongly varies when the chemical desorption is considered or neglected. Conclusions. We show that the chemical desorption process, which directly transforms solid species into gas-phase species, is very efficient for many reactions. Applied to astrophysical environments such as Rho Oph A, we show that the chemical desorption efficiencies derived in this study reproduce the abundances of observed gas-phase methanol, HO2, and H2O2, and that the presence of these molecules in the gas shows the last signs of the evolution of a cloud before the frost.

Open cluster Dolidze 25: Stellar parameters and the metallicity in the Galactic anticentre

The young open cluster Dolidze 25, in the direction of the Galactic Anticentre, has been attributed a very low metallicity, with typical abundances between $-0.5$ and $-0.7$ dex below solar. We intend to derive accurate cluster parameters and accurate stellar abundances for some of its members. We have obtained a large sample of intermediate- and high-resolution spectra for stars in and around Dolidze 25. We used the FASTWIND code to generate stellar atmosphere models to fit the observed spectra. We derive stellar parameters for a large number of OB stars in the area, and abundances of oxygen and silicon for a number of stars with spectral types around B0. We measure low abundances in stars of Dolidze 25. For the three stars with spectral types around B0, we find $0.3$ dex (Si) and $0.5$ dex (O) below the values typical in the solar neighbourhood. These values, even though not as low as those given previously, confirm Dolidze 25 and the surrounding H II region Sh2-284 as the most metal-poor star-forming environment known in the Milky Way. We derive a distance $4.5\pm0.3\:$kpc to the cluster ($r_{\textrm{G}}\approx12.3\:$kpc). The cluster cannot be older than $\sim3\:$Myr, and likely is not much younger. One star in its immediate vicinity, sharing the same distance, has Si and O abundances at most $0.15\:$dex below solar. The low abundances measured in Dolidze 25 are compatible with currently accepted values for the slope of the Galactic metallicity gradient, if we take into account that variations of at least $\pm0.15\:$dex are observed at a given radius. The area traditionally identified as Dolidze 25 is only a small part of a much larger star-forming region that comprises the whole dust shell associated with Sh2-284 and very likely several other smaller H II regions in its vicinity.

Probing star formation in the dense environments of z ∼ 1 lensing haloes aligned with dusty star-forming galaxies detected with the South Pole Telescope

We probe star formation in the environments of massive (similar to 10(13) M-circle dot) dark matter haloes at redshifts of z similar to 1. This star formation is linked to a submillimetre clustering signal which we detect in maps of the Planck High Frequency Instrument that are stacked at the positions of a sample of high redshift (z \textgreater 2) strongly lensed dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope (SPT) 2500 deg(2) survey. The clustering signal has submillimetre colours which are consistent with the mean redshift of the foreground lensing haloes (z similar to 1). We report a mean excess of star formation rate (SFR) compared to the field, of (2700 +/- 700) M-circle dot yr(-1) from all galaxies contributing to this clustering signal within a radius of 3.5 arcmin from the SPT DSFGs. The magnitude of the Planck excess is in broad agreement with predictions of a current model of the cosmic infrared background. The model predicts that 80 per cent of the excess emission measured by Planck originates from galaxies lying in the neighbouring haloes of the lensing halo. Using Herschel maps of the same fields, we find a clear excess, relative to the field, of individual sources which contribute to the Planck excess. The mean excess SFR compared to the field is measured to be (370 +/- 40) M-circle dot yr(-1) per resolved, clustered source. Our findings suggest that the environments around these massive z similar to 1 lensing haloes host intense star formation out to about 2 Mpc. The flux enhancement due to clustering should also be considered when measuring flux densities of galaxies in Planck data.

RECONSIDERING THE EFFECTS OF LOCAL STAR FORMATION ON TYPE Ia SUPERNOVA COSMOLOGY

Recent studies found a correlation with $\sim$3 sigma significance between the local star formation measured by GALEX in Type Ia supernova (SN Ia) host galaxies and the distances or dispersions derived from these SNe. We search for these effects by using data from recent cosmological analyses to greatly increase the SN Ia sample; we include 179 GALEX-imaged SN Ia hosts with distances from the JLA and Pan-STARRS SN Ia cosmology samples and 157 GALEX-imaged SN Ia hosts with distances from the Riess et al. (2011) H$_0$ measurement. We find little evidence that SNe Ia in locally star-forming environments are fainter after light curve correction than SNe Ia in locally passive environments. We find a difference of only 0.000$\pm$0.018 (stat+sys) mag for SNe fit with SALT2 and 0.029$\pm$0.027 (stat+sys) mag for SNe fit with MLCS2k2 (R$_V$ = 2.5), which suggests that proposed changes to recent measurements of H$_0$ and w are not significant and numerically smaller than the parameter measurement uncertainties. We find the greatly reduced significance of these distance modulus differences compared to Rigault et al. (2013) and Rigault et al. (2015) result from two improvements with fairly equal effects, our larger sample size and the use of JLA and Riess et al. (2011) sample selection criteria. Without these improvements, we recover the results of Rigault et al. (2015). We find that both populations have more similar dispersion in distance than found by Rigault et al. (2013), Rigault et al. (2015), and Kelly et al. (2015), with slightly smaller dispersion for locally passive SNe Ia fit with MLCS, the opposite of the effect seen by Rigault et al. (2015) and Kelly et al. (2015). We caution that measuring local SNe Ia environments in the future may require a higher-resolution instrument than GALEX and that SN sample selection has a significant effect on local star formation biases.

The formation of submillimetre-bright galaxies from gas infall over a billion years.

Submillimetre-bright galaxies at high redshift are the most luminous, heavily star-forming galaxies in the Universe and are characterized by prodigious emission in the far-infrared, with a flux of at least five millijanskys at a wavelength of 850 micrometres. They reside in haloes with masses about 10(13) times that of the Sun, have low gas fractions compared to main-sequence disks at a comparable redshift, trace complex environments and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy that simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter haloes have increasing rates of star formation that peak at collective rates of about 500-1,000 solar masses per year at redshifts of two to three, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fuelled in part by the infall of a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a lifetime of nearly a billion years, our simulations show that the submillimetre-bright phase of high-redshift galaxies is prolonged and associated with significant mass buildup in early-Universe proto-clusters, and that many submillimetre-bright galaxies are composed of numerous unresolved components (for which there is some observational evidence).

THE INFLUENCE OF GALAXY SURFACE BRIGHTNESS ON THE MASS–METALLICITY RELATION

We study the effect of surface brightness on the mass-metallicity relation using nearby galaxies whose gas content and metallicity profiles are available. Previous studies using fiber spectra indicated that lower surface brightness galaxies have systematically lower metallicity for their stellar mass, but the results were uncertain because of aperture effect. With stellar masses and surface brightnesses measured at WISE W1 and W2 bands, we re-investigate the surface brightness dependence with spatially-resolved metallicity profiles and find the similar result. We further demonstrate that the systematical difference cannot be explained by the gas content of galaxies. For two galaxies with similar stellar and gas masses, the one with lower surface brightness tends to have lower metallicity. Using chemical evolution models, we investigate the inflow and outflow properties of galaxies of different masses and surface brightnesses. We find that, on average, high mass galaxies have lower inflow and outflow rates relative to star formation rate. On the other hand, lower surface brightness galaxies experience stronger inflow than higher surface brightness galaxies of similar mass. The surface brightness effect is more significant for low mass galaxies. We discuss implications on the different inflow properties between low and high surface brightness galaxies, including star formation efficiency, environment and mass assembly history.