Models Of Galaxy Evolution Research Articles

A Review of Recent Observations of Galactic Winds Driven by Star Formation

Galaxy-scale outflows of gas, or galactic winds (GWs), driven by energy from star formation are a pivotal mechanism for regulation of star formation in the current model of galaxy evolution. Observations of this phenomenon have proliferated through the wide application of old techniques on large samples of galaxies, the development of new methods, and advances in telescopes and instrumentation. I review the diverse portfolio of direct observations of stellar GWs since 2010. Maturing measurements of the ionized and neutral gas properties of nearby winds have been joined by exciting new probes of molecular gas and dust. Low-z techniques have been newly applied in large numbers at high z. The explosion of optical and near-infrared 3D imaging spectroscopy has revealed the complex, multiphase structure of nearby GWs. These observations point to stellar GWs being a common feature of rapidly star-forming galaxies throughout at least the second half of cosmic history, and suggest that scaling relationships between outflow and galaxy properties persist over this period. The simple model of a modest-velocity, biconical flow of multiphase gas and dust perpendicular to galaxy disks continues to be a robust descriptor of these flows.

Evolution of the Dust Composition in Damped Lyα Systems

Abstract We present a method for estimating the relative abundances of refractory elements in the interstellar dust of galaxies hosting damped Lyα (DLA) systems. The method requires gas-phase column densities of volatile and refractory elements, obtained from absorption-line spectroscopy, and interstellar abundances of the same elements, predicted by chemical evolution models of DLA galaxies. We applied this method to the sample of DLA systems with measurements of Mg, Si, S, Fe, and Zn column densities. We find that the dust abundance ratios (Si/Fe)d and (Mg/Fe)d decrease by almost two orders of magnitude in the metallicity range between ≃1/100 solar to roughly solar. This decrease is stronger than the well-known decline of α/Fe ratios with metallicity observed in metal-poor stars and galaxies, suggesting the existence of metallicity-dependent mechanisms of dust production. To cast light on these mechanisms we investigated the contributions of different stellar sources and interstellar processes to the galactic cycle of dust. We find that Type II SNe are important contributors to the dust composition at low metallicity ([Fe/H] < −0.6), whereas dust accretion in the interstellar medium appears to be important at higher metallicities, leading to a gradual rise of iron-rich particles, possibly in metal form. To further investigate the nature of the dust, we introduced an idealized model of dust grains based on a mixture of silicates (pyroxenes and olivines) and an iron-rich constituent. The model reproduces the evolutionary trends and suggests that olivines are dominant in silicates, in line with other studies of interstellar dust composition.

Nature versus nurture: what regulates star formation in satellite galaxies?

We use our state-of-the-art Galaxy Evolution and Assembly (GAEA) semi-analytic model to study how and on which time-scales star formation is suppressed in satellite galaxies. Our fiducial stellar feedback model, implementing strong stellar driven outflows, reproduces relatively well the variations of passive fractions as a function of galaxy stellar mass and halo mass measured in the local Universe, as well as the `quenching' time-scales inferred from the data. We show that the same level of agreement can be obtained by using an alternative stellar feedback scheme featuring lower ejection rates at high redshift, and modifying the treatment for hot gas stripping. This scheme over-predicts the number densities of low to intermediate mass galaxies. In addition, a good agreement with the observed passive fractions can be obtained only by assuming that cooling can continue on satellites, at the rate predicted considering halo properties at infall, even after their parent dark matter substructure is stripped below the resolution of the simulation. For our fiducial model, the better agreement with the observed passive fractions can be ascribed to: (i) a larger cold gas fraction of satellites at the time of accretion, and (ii) a lower rate of gas reheating by supernovae explosions and stellar winds with respect to previous versions of our model. Our results suggest that the abundance of passive galaxies with stellar mass larger than ~10^10 Msun is primarily determined by the self-regulation between star formation and stellar feedback, with environmental processes playing a more marginal role.

Starburst galaxies in semi-analytic models of galaxy formation and evolution

We study the shape and evolution of the star formation main sequence in three independently developed semi-analytic models of galaxy formation. We focus, in particular, on the characterization of the model galaxies that are significantly above the main sequence, and that can be identified with galaxies classified as `starburst' in recent observational work. We find that, in all three models considered, star formation triggered by merger events (both minor and major) contribute to only a very small fraction of the cosmic density of star formation. While mergers are associated to bursts of star formation in all models, galaxies that experienced recent merger events are not necessarily significantly above the main sequence. On the other hand, `starburst galaxies' are not necessarily associated with merger episodes, especially at the low-mass end. Galaxies that experienced recent mergers can have relatively low levels of star formation when/if the merger is gas-poor, and galaxies with no recent merger can experience episodes of starbursts due to a combination of large amount of cold gas available from cooling/accretion events and/or small disk radii which increases the cold gas surface density.

Galaxy and Mass Assembly (GAMA): Accurate number densities and environments of massive ultra-compact galaxies at 0.02 < z < 0.3

Context. Massive ultra-compact galaxies (MUGs) are common at z = 2−3, but very rare in the nearby Universe. Simulations predict that the few surviving MUGs should reside in galaxy clusters, whose large relative velocities prevent them from merging, thus maintaining their original properties (namely stellar populations, masses, sizes and dynamical state). Aims. Our goal is to obtain a complete census of the MUG population at 0.02 < z < 0.3, determining the number density, population properties and environment. Methods. We have taken advantage of the high-completeness, large-area spectroscopic GAMA survey, complementing it with deeper imaging from the KiDS and VIKING surveys. We find a set of 22 bona-fide MUGs, defined as having high stellar mass (> 8 × 1010 M⊙) and compact size (Re < 2 kpc). An additional set of seven lower-mass objects (6 × 1010 < M⋆/M⊙ < 8 × 1010) are also potential candidates according to typical mass uncertainties. Results. The comoving number density of MUGs at low redshift (z < 0.3) is constrained at (1.0 ± 0.4)×10−6 Mpc−3, consistent with galaxy evolution models. However, we find a mixed distribution of old and young galaxies, with a quarter of the sample representing (old) relics. MUGs have a predominantly early or swollen disk morphology (Sérsic index 1 < n < 2.5) with high stellar surface densities (⟨Σe⟩∼1010 M⊙ Kpc−2). Interestingly, a large fraction feature close companions – at least in projection – suggesting that many (but not all) reside in the central regions of groups. Halo masses show these galaxies inhabit average-mass groups. Conclusions. As MUGs are found to be almost equally distributed among environments of different masses, their relative fraction is higher in more massive overdensities, matching the expectations that some of these galaxies fell in these regions at early times. However, there must be another channel leading some of these galaxies to an abnormally low merger history because our sample shows a number of objects that do not inhabit particularly dense environments.

The Ages of M31 Star Clusters: Spectral Energy Distribution versus Color–Magnitude Diagram

Abstract It is well-known that fitting Color–Magnitude Diagrams (CMDs) to theoretical isochrones is the main method to determine star cluster ages. However, when CMDs are not available, the Spectral Energy Distribution (SED)-fitting technique is the only other approach, although it suffers the age–metallicity–reddening degeneracy. In this work, we gather the ages, metallicities, and masses of dozens of M31 star clusters from the CMD-fitting with Hubble Space Telescope images from the literature for comparison. We check the reliability of the SED-fitting results with different models, i.e., Bruzual & Charlot model (BC03), Galaxy Evolutionary Synthesis Models (GALEV), and Advanced Stellar Population Synthesis (ASPS) for the simple stellar populations (SSPs) with single stars (ss)-SSP/binary-star (bs)-SSPs models. The photometry bands includes the Galaxy Evolution Explorer GALEX FUV/NUV bands, optical/near-infrared UBVRIJHK bands, as well as the Wide-field Infrared Survey Explorer (WISE) W1/W2 bands. The comparisons show that the SED-fitting ages agree well with the CMD-fitting ages, either with the fixed metallicity or with the free metallicity for both the BC03 and the GALEV model. However, for the ASPS models, it seems that SED-fitting results are systematically older than the CMD ages, especially for the ages log t < 9.0 (year). The fitting also shows that the GALEX FUV/NUV-band are more important than the WISE W1/W2 for constraining the ages, which confirms the previous findings. We also derived the masses of our sample star clusters from the BC03 and GALEV models and it is found that the values agree well with those in the literature.

The Herschel-PACS North Ecliptic Pole Survey

Abstract A detailed analysis of Herschel/Photoconductor Array Camera and Spectrometer (PACS) observations at the North Ecliptic Pole is presented. High-quality maps, covering an area of 0.44 deg2, are produced and then used to derive potential candidate source lists. A rigorous quality-control pipeline has been used to create final legacy catalogues in the PACS Green 100 μm and Red 160 μm bands, containing 1384 and 630 sources respectively. These catalogues reach to more than twice the depth of the current archival Herschel/PACS Point Source Catalogue, detecting 400 and 270 more sources in the short- and long-wavelength bands, respectively. Galaxy source counts are constructed that extend down to flux densities of 6 mJy and 19 mJy (50% completeness) in the Green 100 μm and Red 160 μm bands, respectively. These source counts are consistent with previously published PACS number counts in other fields across the sky. The source counts are then compared with a galaxy evolution model which identifies a population of luminous infrared galaxies as responsible for the bulk of the galaxy evolution over the flux range (5–100 mJy) spanned by the observed counts, contributing approximate fractions of 50% and 60% to the cosmic infrared background at 100 μm and 160 μm, respectively.

The WISSH quasars project

Winds accelerated by active galactic nuclei (AGNs) are invoked in the most successful models of galaxy evolution to explain the observed physical and evolutionary properties of massive galaxies. Winds are expected to deposit energy and momentum into the interstellar medium (ISM), thus regulating both star formation and supermassive black hole (SMBH) growth. We undertook a multiband observing program aimed at obtaining a complete census of winds in a sample of WISE/SDSS selected hyper-luminous (WISSH) quasars (QSOs) at z ≈ 2–4. We analyzed the rest-frame optical (i.e. LBT/LUCI and VLT/SINFONI) and UV (i.e. SDSS) spectra of 18 randomly selected WISSH QSOs to measure the SMBH mass and study the properties of winds both in the narrow line region (NLR) and broad line region (BLR) traced by blueshifted or skewed [OIII] and CIV emission lines, respectively. These WISSH QSOs are powered by SMBH with masses ≳109 M⊙ accreting at 0.4 < λEdd < 3.1. We found the existence of two subpopulations of hyper-luminous QSOs characterized by the presence of outflows at different distances from the SMBH. One population (i.e. [OIII] sources) exhibits powerful [OIII] outflows, a rest-frame equivalent width (REW) of the CIV emission REWCIV ≈ 20–40 Å, and modest CIV velocity shift (vCIVpeak) with respect to the systemic redshift (vCIVpeak <~ 2000 km s−1). The second population (i.e. Weak [OIII] sources), representing ~70% of the analyzed WISSH QSOs, shows weak or absent [OIII] emission and an extremely large blueshifted CIV emission (vCIVpeak up to ~8000 km s−1 and REWCIV <~ 20 Å). We propose two explanations for the observed behavior of the strength of the [OIII] emission in terms of the orientation effects of the line of sight and ionization cone. The dichotomy in the presence of BLR and NLR winds could be likely due to inclination effects considering a polar geometry scenario for the BLR winds. In a few cases these winds are remarkably as powerful as those revealed in the NLR in the [OIII] QSOs (Ėkin ~ 1044−45 erg s−1). We also investigated the dependence of these CIV winds on fundamental AGN parameters such as bolometric luminosity (LBol), Eddington ratio (λEdd), and UV-to-X-ray continuum slope (αOX). We found a strong correlation with LBol and an anti-correlation with αOX whereby the higher the luminosity, the steeper the ionizing continuum described by means of αOX and the larger the blueshift of the CIV emission line. Finally, the observed dependence vCIVpeak ∝ LBol0.28 ± 0.04 is consistent with a radiatively-driven-winds scenario, where a strong UV continuum is necessary to launch the wind and a weakness of the X-rayemission is fundamental to prevent overionization of the wind itself.

On the influence of environment on star-forming galaxies

We use our state-of-the-art semi analytic model for GAlaxy Evolution and Assembly (GAEA), and observational measurements of nearby galaxies to study the influence of the environment on the gas content and gaseous/stellar disc sizes of star-forming galaxies. We analyse the origin of differences between physical properties of satellites and those of their central counterparts, identified by matching the Vmax of their host haloes at the accretion time of the satellites. Our model reproduces nicely the differences between centrals and satellites measured for the HI mass, size of the star-forming region, and stellar radii. In contrast, our model predicts larger differences with respect to data for the molecular gas mass and star formation rate. By analysing the progenitors of central and satellite model galaxies, we find that differences in the gas content arise after accretion, and can be entirely ascribed to the instantaneous stripping of the hot gas reservoir. The suppression of cold gas replenishment via cooling and star formation leads to a reduction of the cold gas and of its density. Therefore, more molecular gas is lost than lower density HI gas, and model satellites have less molecular gas and lower star formation rates than observed satellites. We argue that these disagreements could be largely resolved with the inclusion of a proper treatment for ram-pressure stripping of cold gas and a more gradual stripping of the hot gas reservoir. A more sophisticated treatment of angular momentum exchanges, accounting for the multi-phase nature of the gaseous disc is also required.

On the origin of N in galaxies with galaxy evolution models

AbstractNitrogen is among the most abundant chemical elements in the cosmos, and asymptotic giant branch (AGB) stars are fundamental nucleosynthetic sources of N in galaxies. In this work, we show how the observed N/O versus O/H chemical abundance diagram, both in extragalactic systems and in our own Galaxy, can be used to constrain the nucleosynthetic origin of N in the cosmos. In particular, we review the results of our studies with chemical evolution models, embedded in full cosmological chemodynamical simulations.

Comparing the star formation history of three nearby disk galaxies

AbstractUnderstanding the effect of environment on galaxy formation and evolution is one of the hot topics in extragalactic astronomy. Here we constructed a chemical evolution model of disk galaxies. By comparing the model predictions with the observed profiles, we investigated the star formation history of M33, NGC 300 and NGC 2403. We found that M33 has much longer infall timescale than NGC 300 and NGC 2403, and the star formation process of M33 is still active at later phase. Our results suggested that the cold gas supply of M33 is sufficient in the present-day, which may originate from the HI bridge between M33 and M31. In other words, we argue that the local environment plays an important role on the star formation history of a galaxy, at least for M33.

Metallicity gradients in high-z galaxies: insights from the KLEVER Survey

AbstractWe present reconstructed source plane metallicity maps for a sample of ~30 gravitationally lensed galaxies between 1.2 < z < 2.5, observed in the framework of the KLEVER Survey. Oxygen abundance is derived exploiting a variety of different emission line diagnostics, as provided by the full coverage of the near-infrared bands. The majority of galaxies in our sample present flat radial metallicity gradients, in agreement with galaxy evolution models predicting strong feedback mechanisms in place at these epochs. However, complex patterns as seen in some of our metallicity maps warn against the use of azimuthally-averaged radial gradients as the only observable to constrain chemical evolution models.

The growth and entrainment of cold gas in a hot wind

ABSTRACT Both absorption and emission-line studies show that cold gas around galaxies is commonly outflowing at speeds of several hundred km s−1. This observational fact poses a severe challenge to our theoretical models of galaxy evolution since most feedback mechanisms (e.g. supernovae feedback) accelerate hot gas, and the time-scale it takes to accelerate a blob of cold gas via a hot wind is much larger than the time it takes to destroy the blob. We revisit this long-standing problem using three-dimensional hydrodynamical simulations with radiative cooling. Our results confirm previous findings that cooling is often not efficient enough to prevent the destruction of cold gas. However, we also identify regions of parameter space where the cooling efficiency of the mixed, ‘warm’ gas is sufficiently large to contribute new comoving cold gas, which can significantly exceed the original cold gas mass. This happens whenever, tcool, mix/tcc < 1, where tcool, mix is the cooling time of the mixed warm gas and tcc is the cloud-crushing time. This criterion is always satisfied for a large enough cloud. Cooling ‘focuses’ stripped material on to the tail where mixing takes place and new cold gas forms. A sufficiently large simulation domain is crucial to capturing this behaviour.

Calibrated, cosmological hydrodynamical simulations with variable IMFs I: Method and effect on global galaxy scaling relations

The recently inferred variations in the stellar initial mass function (IMF) among local high-mass early-type galaxies may require a reinterpretation of observations of galaxy populations and may have important consequences for the predictions of models of galaxy formation and evolution. We present a new pair of cosmological, hydrodynamical simulations based on the EAGLE model that self-consistently adopt an IMF that respectively becomes bottom- or top-heavy in high-pressure environments for individual star-forming gas particles. In such models, the excess stellar mass-to-light ($M/L$) ratio with respect to a reference IMF is increased due to an overabundance of low-mass dwarf stars or stellar remnants, respectively. Crucially, both pressure-dependent IMFs have been calibrated to reproduce the observed trends of increasing excess $M/L$ with central stellar velocity dispersion ($\sigma_e$) in early-type galaxies, while maintaining agreement with the observables used to calibrate the EAGLE model, namely the galaxy luminosity function, half-light radii of late-type galaxies, and black hole masses. We find that while the $M/L$ excess is a good measure of the IMF for low-mass slope variations, it depends strongly on the age of the stellar population for high-mass slope variations. The normalization of the [Mg/Fe]$-\sigma_e$ relation is decreased (increased) for bottom- (top-)heavy IMF variations, while the slope is not strongly affected. Bottom-heavy variations have little impact on galaxy metallicities, half-light radii of early-type galaxies, or star formation rates, while top-heavy variations significantly increase these quantities for high-mass galaxies, leading to tension with observations.

The indirect influence of quasars on reionization

ABSTRACT The exact role of quasars during the Epoch of Reionization remains uncertain. With consensus leaning towards quasars producing a negligible amount of ionizing photons, we pose an alternate question: Can quasars indirectly contribute to reionization by allowing ionizing photons from stars to escape more easily? Using the semi-analytic galaxy evolution model to evolve a galaxy population through cosmic time, we construct an idealized scenario in which the escape fraction of stellar ionizing photons (fesc) is boosted following quasar wind events, potentially for several dynamical times. We find that under this scenario, the mean value of fesc as a function of galaxy stellar mass peaks for intermediate mass galaxies. This mass dependence will have consequences for the 21-cm power spectrum, enhancing power at small scales and suppressing it at large scales. This hints that whilst quasars may not directly contribute to the ionizing photon budget, they could influence reionization indirectly by altering the topology of ionized regions.

GAMA/H-ATLAS: the local dust mass function and cosmic density as a function of galaxy type – a benchmark for models of galaxy evolution

ABSTRACT We present the dust mass function (DMF) of 15 750 galaxies with redshift $z$ < 0.1, drawn from the overlapping area of the GAMA and H-ATLAS surveys. The DMF is derived using the density corrected Vmax method, where we estimate Vmax using: (i) the normal photometric selection limit (pVmax) and (ii) a bivariate brightness distribution (BBD) technique, which accounts for two selection effects. We fit the data with a Schechter function, and find $M^{*}=(4.65 \pm 0.18)\times 10^{7}\,h^2_{70}\, \mathrm{ M}_{\odot }$, α = (−1.22 ± 0.01), $\phi ^{*}=(6.26 \pm 0.28)\times 10^{-3}\,h^3_{70}\,\rm Mpc^{-3}\,dex^{-1}$. The resulting dust mass density parameter integrated down to 104 M⊙ is Ωd = (1.11 ± 0.02) × 10−6 which implies the mass fraction of baryons in dust is $f_{m_\mathrm{ b}}=(2.40\pm 0.04)\times 10^{-5}$; cosmic variance adds an extra 7–17 per cent uncertainty to the quoted statistical errors. Our measurements have fewer galaxies with high dust mass than predicted by semi-analytic models. This is because the models include too much dust in high stellar mass galaxies. Conversely, our measurements find more galaxies with high dust mass than predicted by hydrodynamical cosmological simulations. This is likely to be from the long time-scales for grain growth assumed in the models. We calculate DMFs split by galaxy type and find dust mass densities of Ωd = (0.88 ± 0.03) × 10−6 and Ωd = (0.060 ± 0.005) × 10−6 for late types and early types, respectively. Comparing to the equivalent galaxy stellar mass functions (GSMF) we find that the DMF for late types is well matched by the GSMF scaled by (8.07 ± 0.35) × 10−4.

Towards a census of high-redshift dusty galaxies with Herschel

Context. Over the last decade a large number of dusty star-forming galaxies has been discovered up to redshift z = 2 − 3 and recent studies have attempted to push the highly confused Herschel SPIRE surveys beyond that distance. To search for z ≥ 4 galaxies they often consider the sources with fluxes rising from 250 μm to 500 μm (so-called “500 μm-risers”). Herschel surveys offer a unique opportunity to efficiently select a large number of these rare objects, and thus gain insight into the prodigious star-forming activity that takes place in the very distant Universe. Aims. We aim to implement a novel method to obtain a statistical sample of 500 μm-risers and fully evaluate our selection inspecting different models of galaxy evolution. Methods. We consider one of the largest and deepest Herschel surveys, the Herschel Virgo Cluster Survey. We develop a novel selection algorithm which links the source extraction and spectral energy distribution fitting. To fully quantify selection biases we make end-to-end simulations including clustering and lensing. Results. We select 133 500 μm-risers over 55 deg2, imposing the criteria: S500 > S350 > S250, S250 > 13.2 mJy and S500 > 30 mJy. Differential number counts are in fairly good agreement with models, displaying a better match than other existing samples. The estimated fraction of strongly lensed sources is 24+6-5% based on models. Conclusions. We present the faintest sample of 500 μm-risers down to S250 = 13.2 mJy. We show that noise and strong lensing have an important impact on measured counts and redshift distribution of selected sources. We estimate the flux-corrected star formation rate density at 4 < z < 5 with the 500 μm-risers and find it to be close to the total value measured in far-infrared. This indicates that colour selection is not a limiting effect to search for the most massive, dusty z > 4 sources.

Exploring stellar evolution with gravitational-wave observations

Recent detections of gravitational waves from merging binary black holes opened new possibilities to study the evolution of massive stars and black hole formation. In particular, stellar evolution models may be constrained on the basis of the differences in the predicted distribution of black hole masses and redshifts. In this work we propose a framework that combines galaxy and stellar evolution models and use it to predict the detection rates of merging binary black holes for various stellar evolution models. We discuss the prospects of constraining the shape of the time delay distribution of merging binaries using just the observed distribution of chirp masses. Finally, we consider a generic model of primordial black hole formation and discuss the possibility of distinguishing it from stellar-origin black holes.

Post-Newtonian evolution of massive black hole triplets in galactic nuclei – III. A robust lower limit to the nHz stochastic background of gravitational waves

Inspiraling massive black-hole binaries (MBHBs) forming in the aftermath of galaxy mergers are expected to be the loudest gravitational-wave (GW) sources relevant for pulsar-timing arrays (PTAs) at nHz frequencies. The incoherent overlap of signals from a cosmic population of MBHBs gives rise to a stochastic GW background (GWB) with characteristic strain around $h_c\sim10^{-15}$ at a reference frequency of 1 yr$^{-1}$, although uncertainties around this value are large. Current PTAs are piercing into the GW amplitude range predicted by MBHB-population models, but no detection has been reported so far. To assess the future success prospects of PTA experiments, it is therefore important to estimate the minimum GWB level consistent with our current understanding of the formation and evolution of galaxies and massive black holes (MBHs). To this purpose, we couple a semianalytic model of galaxy evolution and an extensive study of the statistical outcome of triple MBH interactions. We show that even in the most pessimistic scenario where all MBHBs stall before entering the GW-dominated regime, triple interactions resulting from subsequent galaxy mergers inevitably drive a considerable fraction of the MBHB population to coalescence. At frequencies relevant for PTA, the resulting GWB is only a factor of 2-to-3 suppressed compared to a fiducial model where binaries are allowed to merge over Gyr timescales. Coupled with current estimates of the expected GWB amplitude range, our findings suggest that the minimum GWB from cosmic MBHBs is unlikely to be lower than $h_c\sim 10^{-16}$ (at $f = 1$ yr$^{-1}$), well within the expected sensitivity of projected PTAs based on future observations with FAST, MeerKAT and SKA.

Red, redder, reddest: SCUBA-2 imaging of colour-selected Herschel sources

High-redshift, luminous, dusty star forming galaxies (DSFGs) constrain the extremity of galaxy formation theories. The most extreme are discovered through follow-up on candidates in large area surveys. Here we present 850 $\mu$m SCUBA-2 follow-up observations of 188 red DSFG candidates from the \textit{Herschel} Multi-tiered Extragalactic Survey (HerMES) Large Mode Survey, covering 274 deg$^2$. We detected 87 per cent with a signal-to-noise ratio $>$ 3 at 850~$\mu$m. We introduce a new method for incorporating the confusion noise in our spectral energy distribution fitting by sampling correlated flux density fluctuations from a confusion limited map. The new 850~$\mu$m data provide a better constraint on the photometric redshifts of the candidates, with photometric redshift errors decreasing from $\sigma_z/(1+z)\approx0.21$ to $0.15$. Comparison spectroscopic redshifts also found little bias ($\langle (z-z_{\rm spec})/(1+z_{\rm spec})\rangle = 0.08 $). The mean photometric redshift is found to be 3.6 with a dispersion of $0.4$ and we identify 21 DSFGs with a high probability of lying at $z > 4$. After simulating our selection effects we find number counts are consistent with phenomenological galaxy evolution models. There is a statistically significant excess of WISE-1 and SDSS sources near our red galaxies, giving a strong indication that lensing may explain some of the apparently extreme objects. Nevertheless, our sample should include examples of galaxies with the highest star formation rates in the Universe ($\gg10^3$ M$_\odot$yr$^{-1}$).