Source author record

Alan R. Duffy

Alan R. Duffy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

Gamma-ray emission from the Sagittarius Dwarf Spheroidal galaxy due to millisecond pulsars

The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show via a spatial template analysis that much of the gamma-ray emission associated to the brightest region of substructure -- the so-called cocoon -- is likely due to the Sagittarius dwarf spheroidal (Sgr dSph) galaxy. This large Milky Way satellite is viewed through the Fermi Bubbles from the position of the Solar System. As a tidally and ram-pressure stripped remnant, the Sgr dSph has no on-going star formation, but we nevertheless demonstrate that the dwarf's millisecond pulsar (MSP) population can plausibly supply the gamma-ray signal that our analysis associates to its stellar template. The measured spectrum is naturally explained by inverse Compton scattering of cosmic microwave background photons by high-energy electron-positron pairs injected by MSPs belonging to the Sgr dSph, combined with these objects' magnetospheric emission. This finding plausibly suggests that MSPs produce significant gamma-ray emission amongst old stellar populations, potentially confounding indirect dark matter searches in regions such as the Galactic Centre, the Andromeda galaxy, and other massive Milky Way dwarf spheroidals.

preprint2021arXiv

Fast radio bursts as probes of feedback from active galactic nuclei

Fast Radio Bursts (FRBs) are a promising tool for studying the low-density universe as their dispersion measures (DM) are extremely sensitive probes of electron column density. Active Galactic Nuclei (AGN) inject energy into the intergalactic medium, affecting the DM and its scatter. To determine the effectiveness of FRBs as a probe of AGN feedback, we analysed three different AGN models from the EAGLE simulation series. We measured the mean DM-redshift relation, and the scatter around it, using $2.56 \times 10^8$ sightlines at 131 redshift ($z$) bins between $0 \leq z \leq 3$. While the DM-redshift relation itself is highly robust against different AGN feedback models, significant differences are detected in the scatter around the mean: weaker feedback leads to more scatter. We find $\sim 10^4$ localised FRBs are needed to discriminate between the scatter in standard feedback and stronger, more intermittent feedback models. The number of FRBs required is dependent on the redshift distribution of the detected population. A log-normal redshift distribution at $z=0.5$ requires approximately 50% fewer localised FRBs than a distribution centred at $z=1$. With the Square Kilometre Array expected to detect $>10^3$ FRBs per day, in the future, FRBs will be able to provide constraints on AGN feedback.

preprint2020arXiv

Cosmic Ray Extremely Distributed Observatory

The Cosmic Ray Extremely Distributed Observatory (CREDO) is a newly formed, global collaboration dedicated to observing and studying cosmic rays (CR) and cosmic ray ensembles (CRE): groups of a minimum of two CR with a common primary interaction vertex or the same parent particle. The CREDO program embraces testing known CR and CRE scenarios, and preparing to observe unexpected physics, it is also suitable for multi-messenger and multi-mission applications. Perfectly matched to CREDO capabilities, CRE could be formed both within classical models (e.g. as products of photon-photon interactions), and exotic scenarios (e.g. as results of decay of Super Heavy Dark Matter particles). Their fronts might be significantly extended in space and time, and they might include cosmic rays of energies spanning the whole cosmic ray energy spectrum, with a footprint composed of at least two extensive air showers with correlated arrival directions and arrival times. Since CRE are mostly expected to be spread over large areas and, because of the expected wide energy range of the contributing particles, CRE detection might only be feasible when using available cosmic ray infrastructure collectively, i.e. as a globally extended network of detectors. Thus, with this review article, the CREDO Collaboration invites the astroparticle physics community to actively join or to contribute to the research dedicated to CRE, and in particular to share any cosmic ray data useful for the specific CRE detection strategies.

preprint2020arXiv

Search for ultra-high energy photons through preshower effect with gamma-ray telescopes: Study of CTA-North efficiency

We study the feasibility of detecting preshower initiated by ultra-high energy photons using Monte-Carlo simulations of nearly horizontal air showers for the example of the La Palma site of the Cherenkov Telescope Array. We investigate the efficiency of multivariate analysis in correctly identifying preshower events initiated by 40 EeV photons and cosmic-ray dominated background simulated in the energy range 10 TeV -- 10 EeV. The effective areas for such kind of events are also investigated and event rate predictions related to different ultra-high energy photons production models are presented. While the expected number of preshowers from diffuse emission of UHE photon for 30 hours of observation is estimated around $3.3\times10^{-5}$ based on the upper limits put by the Pierre Auger Observatory, this value is at the level of $2.7\times10^{-4}$ ($5.7\times 10^{-5}$) when considering the upper limits of the Pierre Auger Observatory (Telescope Array) on UHE photon point sources. However, UHE photon emission may undergo possible "boosting" due to gamma-ray burst, increasing the expected number of preshower events up to 0.17 and yielding a minimum required flux of $\sim 0.2$ $\mathrm{km^{-2}yr^{-1}}$ to obtain one preshower event, which is about a factor 10 higher than upper limits put by the Pierre Auger Observatory and Telescope Array (0.034 and 0.019 $\mathrm{km^{-2}yr^{-1}}$, respectively).

preprint2016arXiv

Dark-ages Reionization & Galaxy Formation Simulation I: The dynamical lives of high redshift galaxies

We present the Dark-ages Reionization and Galaxy-formation Observables from Numerical Simulations (DRAGONS) program and Tiamat, the collisionless N-body simulation program upon which DRAGONS is built. The primary trait distinguishing Tiamat from other large simulation programs is its density of outputs at high redshift (100 from z=35 to z=5; roughly one every 10 Myr) enabling the construction of very accurate merger trees at an epoch when galaxy formation is rapid and mergers extremely frequent. We find that the friends-of-friends halo mass function agrees well with the prediction of Watson et al. at high masses, but deviates at low masses, perhaps due to our use of a different halo finder or perhaps indicating a break from "universal" behaviour. We then analyse the dynamical evolution of galaxies during the Epoch of Reionization finding that only a small fraction (~20%) of galactic halos are relaxed. We illustrate this using standard relaxation metrics to establish two dynamical recovery time-scales: i) halos need ~1.5 dynamical times following formation, and ii) ~2 dynamical times following a major (3:1) or minor (10:1) merger to be relaxed. This is remarkably consistent across a wide mass range. Lastly, we use a phase-space halo finder to illustrate that major mergers drive long-lived massive phase-space structures which take many dynamical times to dissipate. This can yield significant differences in the inferred mass build-up of galactic halos and we suggest that care must be taken to ensure a physically meaningful match between the galaxy-formation physics of semi-analytic models and the halo finders supplying their input.

preprint2016arXiv

Dark-ages reionization & galaxy formation simulation IV: UV luminosity functions of high-redshift galaxies

In this paper we present calculations of the UV luminosity function from the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) project, which combines N-body, semi-analytic and semi-numerical modelling designed to study galaxy formation during the Epoch of Reionization. Using galaxy formation physics including supernova feedback, the model naturally reproduces the UV LFs for high-redshift star-forming galaxies from $z{\sim}5$ through to $z{\sim}10$. We investigate the luminosity--star formation rate (SFR) relation, finding that variable SFR histories of galaxies result in a scatter around the median relation of $0.1$--$0.3$ dex depending on UV luminosity. We find close agreement between the model and observationally derived SFR functions. We use our calculated luminosities to investigate the luminosity function below current detection limits, and the ionizing photon budget for reionization. We predict that the slope of the UV LF remains steep below current detection limits and becomes flat at $M_\mathrm{UV}{\gtrsim}{-14}$. We find that $48$ ($17$) per cent of the total UV flux at $z{\sim}6$ ($10$) has been detected above an observational limit of $M_\mathrm{UV}{\sim}{-17}$, and that galaxies fainter than $M_\mathrm{UV}{\sim}{-17}$ are the main source of ionizing photons for reionization. We investigate the luminosity--stellar mass relation, and find a correlation for galaxies with $M_\mathrm{UV}{<}{-14}$ that has the form $M_*{\propto}10^{-0.47M_\mathrm{UV}}$, in good agreement with observations, but which flattens for fainter galaxies. We determine the luminosity--halo mass relation to be $M_\mathrm{vir}{\propto}10^{-0.35M_\mathrm{UV}}$, finding that galaxies with $M_\mathrm{UV}{=}{-20}$ reside in host dark matter haloes of $10^{11.0\pm 0.1}\mathrm{M_\odot}$ at $z{\sim}6$, and that this mass decreases towards high redshift.

preprint2016arXiv

Dark-ages reionization and galaxy formation simulation - III. Modelling galaxy formation and the epoch of reionization

We introduce Meraxes, a new, purpose-built semi-analytic galaxy formation model designed for studying galaxy growth during reionization. Meraxes is the first model of its type to include a temporally and spatially coupled treatment of reionization and is built upon a custom (100 Mpc)$^3$ N-body simulation with high temporal and mass resolution, allowing us to resolve the galaxy and star formation physics relevant to early galaxy formation. Our fiducial model with supernova feedback reproduces the observed optical depth to electron scattering and evolution of the galaxy stellar mass function between $z$=5 and 7, predicting that a broad range of halo masses contribute to reionization. Using a constant escape fraction and global recombination rate, our model is unable to simultaneously match the observed ionizing emissivity at $z{\lesssim}6$. However, the use of an evolving escape fraction of 0.05-0.1 at $z{\sim}6$, increasing towards higher redshift, is able to satisfy these three constraints. We also demonstrate that photoionization suppression of low mass galaxy formation during reionization has only a small effect on the ionization history of the inter-galactic medium. This lack of `self-regulation' arises due to the already efficient quenching of star formation by supernova feedback. It is only in models with gas supply-limited star formation that reionization feedback is effective at regulating galaxy growth. We similarly find that reionization has only a small effect on the stellar mass function, with no observationally detectable imprint at $M_{\rm *}{>}10^{7.5}\,{\rm M_{\odot}}$. However, patchy reionization has significant effects on individual galaxy masses, with variations of factors of 2-3 at $z$=5 that correlate with environment.

preprint2016arXiv

The Theoretical Astrophysical Observatory: Cloud-Based Mock Galaxy Catalogues

We introduce the Theoretical Astrophysical Observatory (TAO), an online virtual laboratory that houses mock observations of galaxy survey data. Such mocks have become an integral part of the modern analysis pipeline. However, building them requires an expert knowledge of galaxy modelling and simulation techniques, significant investment in software development, and access to high performance computing. These requirements make it difficult for a small research team or individual to quickly build a mock catalogue suited to their needs. To address this TAO offers access to multiple cosmological simulations and semi-analytic galaxy formation models from an intuitive and clean web interface. Results can be funnelled through science modules and sent to a dedicated supercomputer for further processing and manipulation. These modules include the ability to (1) construct custom observer light-cones from the simulation data cubes; (2) generate the stellar emission from star formation histories, apply dust extinction, and compute absolute and/or apparent magnitudes; and (3) produce mock images of the sky. All of TAO's features can be accessed without any programming requirements. The modular nature of TAO opens it up for further expansion in the future.

preprint2015arXiv

Dark-ages Reionization & Galaxy Formation Simulation II: Spin and concentration parameters for dark matter haloes during the Epoch of Reionization

We use high resolution N-Body simulations to study the concentration and spin parameters of dark matter haloes in the mass range $10^8\, {\rm M}_{\odot}\, h^{-1} < {\rm M} < 10^{11}\, {\rm M}_{\odot}\, h^{-1}$ and redshifts $5{<}z{<}10$, corresponding to the haloes of galaxies thought to be responsible for reionization. We build a sub-sample of equilibrium haloes and contrast their properties to the full population that also includes unrelaxed systems. Concentrations are calculated by fitting both NFW and Einasto profiles to the spherically-averaged density profiles of individual haloes. After removing haloes that are out-of-equilibrium, we find a $z{>}5$ concentration$-$mass ($c(M)$) relation that is almost flat and well described by a simple power-law for both NFW and Einasto fits. The intrinsic scatter around the mean relation is $Δc_{\rm{vir}}{\sim1}$ (or 20 per cent) at $z=5$. We also find that the analytic model proposed by Ludlow et al. reproduces the mass and redshift-dependence of halo concentrations. Our best-fit Einasto shape parameter, $α$, depends on peak height, $ν$, in a manner that is accurately described by $α{=}0.0070ν^2{+}0.1839$. The distribution of the spin parameter, $λ$, has a weak dependence on equilibrium state; $λ$ peaks at roughly ${\sim}0.033$ for our relaxed sample, and at ${\sim}0.04$ for the full population. The spin--virial mass relation has a mild negative correlation at high redshift.

preprint2015arXiv

The accretion history of dark matter halos I: The physical origin of the universal function

Understanding the universal accretion history of dark matter halos is the first step towards determining the origin of their structure. We use the extended Press-Schechter formalism to derive the halo mass accretion history from the growth rate of initial density perturbations. We show that the halo mass history is well described by an exponential function of redshift in the high-redshift regime. However, in the low-redshift regime the mass history follows a power law because the growth of density perturbations is halted in the dark energy dominated era due to the accelerated expansion of the Universe. We provide an analytic model that follows the expression $M(z) = M_{0}(1+z)^{af(M_{0})}e^{-f(M_{0})z}$, where $M_{0} = M(z=0)$, $a$ depends on cosmology and $f(M_{0})$ depends only on the linear matter power spectrum. The analytic model does not rely on calibration against numerical simulations and is suitable for any cosmology. We compare our model with the latest empirical models for the mass accretion history in the literature and find very good agreement. We provide numerical routines for the model online.

preprint2015arXiv

The accretion history of dark matter halos II: The connections with the mass power spectrum and the density profile

We explore the relation between the structure and mass accretion histories of dark matter halos using a suite of cosmological simulations. We confirm that the formation time, defined as the time when the virial mass of the main progenitor equals the mass enclosed within the scale radius, correlates strongly with concentration. We provide a semi-analytic model for halo mass history that combines analytic relations with fits to simulations. This model has the functional form, $M(z) = M_{0}(1+z)^αe^{βz}$, where the parameters $α$ and $β$ are directly correlated with concentration. We then combine this model for the halo mass history with the analytic relations between $α$, $β$ and the linear power spectrum derived by Correa et al. (2014) to establish the physical link between halo concentration and the initial density perturbation field. Finally, we provide fitting formulas for the halo mass history as well as numerical routines, we derive the accretion rate as a function of halo mass, and we demonstrate how the halo mass history depends on cosmology and the adopted definition of halo mass.

preprint2015arXiv

The accretion history of dark matter halos III: A physical model for the concentration-mass relation

We present a semi-analytic, physically motivated model for dark matter halo concentration as a function of halo mass and redshift. The semi-analytic model combines an analytic model for the halo mass accretion history (MAH), based on extended Press Schechter (EPS) theory, with an empirical relation between concentration and formation time obtained through fits to the results of numerical simulations. Because the semi-analytic model is based on EPS theory, it can be applied to wide ranges in mass, redshift and cosmology. The resulting concentration-mass (c-M) relations are found to agree with the simulations, and because the model applies only to relaxed halos, they do not exhibit the upturn at high masses or high redshifts found by some recent works. We predict a change of slope in the z=0 c-M relation at a mass scale of $10^{11}\rm{M}_{\odot}$. We find that this is due to the change in the functional form of the halo MAH, which goes from being dominated by an exponential (for high-mass halos) to a power-law (for low-mass halos). During the latter phase, the core radius remains approximately constant, and the concentration grows due to the drop of the background density. We also analyse how the c-M relation predicted by this work affects the power produced by dark matter annihilation, finding that at z = 0 the power is two orders of magnitude lower than that obtained from extrapolating best-fitting c-M relations. We provide fits to the c-M relations as well as numerical routines to compute concentrations and MAHs.

preprint2014arXiv

Low-mass galaxy formation and the ionizing photon budget during reionization

We use high-resolution simulations of cosmological volumes to model galaxy formation at high-redshift, with the goal of studying the photon budget for reionization. We demonstrate that galaxy formation models that include a strong, thermally coupled supernovae scheme reproduce current observations of star formation rates and specific star formation rates, both during and after the reionization era. These models produce enough UV photons to sustain reionization at z<8 (z<6) through a significant population of faint, unobserved, galaxies for an assumed escape fraction of 20% (5%). This predicted population is consistent with extrapolation of the faint end of observed UV luminosity functions. We find that heating from a global UV/X-ray background after reionization causes a dip in the total global star formation rate density in galaxies below the current observational threshold. Finally, while the currently observed specific star formation rates are incapable of differentiating between supernovae feedback models, sufficiently deep observations will be able to use this diagnostic in the future to investigate galaxy formation at high redshift.

preprint2014arXiv

Probing the nature of dark energy through galaxy redshift surveys with radio telescopes

Galaxy redshift surveys using optical telescopes have, in combination with other cosmological probes, enabled precision measurements of the nature of dark energy. We show that radio telescopes are rapidly becoming competitive with optical facilities in spectroscopic surveys of large numbers of galaxies. Two breakthroughs are driving this change. Firstly, individual radio telescopes are more efficient at mapping the sky thanks to the large field-of-view of new phased-array feeds. Secondly, ever more dishes can be correlated in a cost-effective manner with rapid increases in computing power. The next decade will see the coming of age of the 21cm radio wavelength as a cosmological probe as first the Pathfinders then, ultimately, the Square Kilometre Array is constructed. The latter will determine precise 3D positions for a billion galaxies, mapping the distribution of matter in the Universe over the last 12 billion years. This radio telescope will be able to constrain the equation of state of dark energy, and its potential evolution, to a precision rivalling that of future optical facilities such as DESI and Euclid.

preprint2012arXiv

Cosmological surveys with the Australian Square Kilometre Array Pathfinder

This is a design study into the capabilities of the Australian Square Kilometre Array Pathfinder in performing a full-sky low redshift neutral hydrogen survey, termed WALLABY, and the potential cosmological constraints one can attain from measurement of the galaxy power spectrum. We find that the full sky survey will likely attain 0.6 million redshifts which, when combined with expected Planck CMB data, will constrain the Dark Energy equation of state to 20%, representing a coming of age for radio observations in creating cosmological constraints.

preprint2012arXiv

Predictions for ASKAP Neutral Hydrogen Surveys

The Australian Square Kilometer Array Pathfinder (ASKAP) will revolutionise our knowledge of gas-rich galaxies in the Universe. Here we present predictions for two proposed extragalactic ASKAP neutral hydrogen (HI) emission-line surveys, based on semi-analytic models applied to cosmological N-body simulations. The ASKAP HI All-Sky Survey, known as WALLABY, is a shallow 3 Pi survey (z = 0 - 0.26) which will probe the mass and dynamics of over 600,000 galaxies. A much deeper small-area HI survey, called DINGO, aims to trace the evolution of HI from z = 0 - 0.43, a cosmological volume of 40 million Mpc^3, detecting potentially 100,000 galaxies. The high-sensitivity 30 antenna ASKAP core (diameter ~2 km) will provide an angular resolution of 30 arcsec (at z=0). Our simulations show that the majority of galaxies detected in WALLABY (87.5%) will be resolved. About 5000 galaxies will be well resolved, i.e. more than five beams (2.5 arcmin) across the major axis, enabling kinematic studies of their gaseous disks. This number would rise to 160,000 galaxies if all 36 ASKAP antennas could be used; the additional six antennas provide baselines up to 6 km, resulting in an angular resolution of 10 arcsec. For DINGO this increased resolution is highly desirable to minimise source confusion; reducing confusion rates from a maximum of 10% of sources at the survey edge to 3%. We estimate that the sources detected by WALLABY and DINGO will span four orders of magnitude in total halo mass (from 10^{11} to 10^{15} Msol) and nearly seven orders of magnitude in stellar mass (from 10^{5} to 10^{12} Msol), allowing us to investigate the process of galaxy formation across the last four billion years.

preprint2011arXiv

Dark matter halo concentrations in the Wilkinson Microwave Anisotropy Probe year 5 cosmology

We use a combination of three large N-body simulations to investigate the dependence of dark matter halo concentrations on halo mass and redshift in the WMAP year 5 cosmology. The median relation between concentration and mass is adequately described by a power-law for halo masses in the range 10^11 - 10^15 Msol/h and redshifts z < 2, regardless of whether the halo density profiles are fit using NFW or Einasto profiles. Compared with recent analyses of the Millennium Simulation, which uses a value of sigma_8 that is higher than allowed by WMAP5, z = 0 halo concentrations are reduced by factors ranging from 23 per cent at 10^11 Msol/h to 16 per cent at 10^14 Msol/h. The predicted concentrations are much lower than inferred from X-ray observations of groups and clusters.

preprint2011arXiv

Modelling neutral hydrogen in galaxies using cosmological hydrodynamical simulations

The characterisation of the atomic and molecular hydrogen content of high-redshift galaxies is a major observational challenge that will be addressed over the coming years with a new generation of radio telescopes. We investigate this important issue by considering the states of hydrogen across a range of structures within high-resolution cosmological hydrodynamical simulations. Additionally, our simulations allow us to investigate the sensitivity of our results to numerical resolution and to sub-grid baryonic physics (especially feedback from supernovae and active galactic nuclei). We find that the most significant uncertainty in modelling the neutral hydrogen distribution arises from our need to model a self-shielding correction in moderate density regions. Future simulations incorporating radiative transfer schemes will be vital to improve on our empirical self-shielding threshold. Irrespective of the exact nature of the threshold we find that while the atomic hydrogen mass function evolves only mildly from redshift two to zero, the molecular hydrogen mass function increases with increasing redshift, especially at the high-mass end. Interestingly, the weak evolution of the neutral hydrogen mass function is insensitive to the feedback scheme utilised, but the opposite is true for the molecular gas, which is more closely associated with the star formation in the simulations.

preprint2010arXiv

Impact of baryon physics on dark matter structures: a detailed simulation study of halo density profiles

The back-reaction of baryons on the dark matter halo density profile is of great interest, not least because it is an important systematic uncertainty when attempting to detect the dark matter. Here, we draw on a large suite of high resolution cosmological hydrodynamical simulations, to systematically investigate this process and its dependence on the baryonic physics associated with galaxy formation. The inclusion of baryons results in significantly more concentrated density profiles if radiative cooling is efficient and feedback is weak. The dark matter halo concentration can in that case increase by as much as 30 (10) per cent on galaxy (cluster) scales. The most significant effects occur in galaxies at high redshift, where there is a strong anti-correlation between the baryon fraction in the halo centre and the inner slope of both the total and the dark matter density profiles. If feedback is weak, isothermal inner profiles form, in agreement with observations of massive, early-type galaxies. However, we find that AGN feedback, or extremely efficient feedback from massive stars, is necessary to match observed stellar fractions in groups and clusters, as well as to keep the maximum circular velocity similar to the virial velocity as observed for disk galaxies. These strong feedback models reduce the baryon fraction in galaxies by a factor of 3 relative to the case with no feedback. The AGN is even capable of reducing the baryon fraction by a factor of 2 in the inner region of group and cluster haloes. This in turn results in inner density profiles which are typically shallower than isothermal and the halo concentrations tend to be lower than in the absence of baryons.

preprint2009arXiv

The physics driving the cosmic star formation history

We investigate the physics driving the cosmic star formation (SF) history using the more than fifty large, cosmological, hydrodynamical simulations that together comprise the OverWhelmingly Large Simulations (OWLS) project. We systematically vary the parameters of the model to determine which physical processes are dominant and which aspects of the model are robust. Generically, we find that SF is limited by the build-up of dark matter haloes at high redshift, reaches a broad maximum at intermediate redshift, then decreases as it is quenched by lower cooling rates in hotter and lower density gas, gas exhaustion, and self-regulated feedback from stars and black holes. The higher redshift SF is therefore mostly determined by the cosmological parameters and to a lesser extent by photo-heating from reionization. The location and height of the peak in the SF history, and the steepness of the decline towards the present, depend on the physics and implementation of stellar and black hole feedback. Mass loss from intermediate-mass stars and metal-line cooling both boost the SF rate at late times. Galaxies form stars in a self-regulated fashion at a rate controlled by the balance between, on the one hand, feedback from massive stars and black holes and, on the other hand, gas cooling and accretion. Paradoxically, the SF rate is highly insensitive to the assumed SF law. This can be understood in terms of self-regulation: if the SF efficiency is changed, then galaxies adjust their gas fractions so as to achieve the same rate of production of massive stars. Self-regulated feedback from accreting black holes is required to match the steep decline in the observed SF rate below redshift two, although more extreme feedback from SF, for example in the form of a top-heavy IMF at high gas pressures, can help.