Shocks, cooling and the origin of star formation rates in spiral galaxies

Understanding star formation is problematic as it originates in the large scale dynamics of a galaxy but occurs on the small scale of an individual star forming event. This paper presents the first numerical simulations to resolve the star formation process on sub-parsec scales, whilst also following the dynamics of the interstellar medium (ISM) on galactic scales. In these models, the warm low density ISM gas flows into the spiral arms where orbit crowding produces the shock formation of dense clouds, held together temporarily by their external pressure. Cooling allows the gas to be compressed to sufficiently high densities that local regions collapse under their own gravity and form stars. The star formation rates follow a Schmidt-Kennicutt \Sigma_{SFR} ~ \Sigma_{gas}^{1.4} type relation with the local surface density of gas while following a linear relation with the cold and dense gas. Cooling is the primary driver of star formation and the star formation rates as it determines the amount of cold gas available for gravitational collapse. The star formation rates found in the simulations are offset to higher values relative to the extragalactic values, implying a constant reduction, such as from feedback or magnetic fields, is likely to be required. Intriguingly, it appears that a spiral or other convergent shock and the accompanying thermal instability can explain how star formation is triggered, generate the physical conditions of molecular clouds and explain why star formation rates are tightly correlated to the gas properties of galaxies.Comment: 13 pages, 12 figures. MNRAS in pres

Detections and Constraints on White Dwarf Variability from Time-Series GALEX Observations

We search for photometric variability in more than 23,000 known and candidate white dwarfs, the largest ultraviolet survey compiled for a single study of white dwarfs. We use gPhoton, a publicly available calibration/reduction pipeline, to generate time-series photometry of white dwarfs observed by GALEX. By implementing a system of weighted metrics, we select sources with variability due to pulsations and eclipses. Although GALEX observations have short baselines (< 30 min), we identify intrinsic variability in sources as faint as Gaia G = 20 mag. With our ranking algorithm, we identify 49 new variable white dwarfs (WDs) in archival GALEX observations. We detect 41 new pulsators: 37 have hydrogen-dominated atmospheres (DAVs), including one possible massive DAV, and four are helium-dominated pulsators (DBVs). We also detect eight new eclipsing systems; five are new discoveries, and three were previously known spectroscopic binaries. We perform synthetic injections of the light curve of WD 1145+017, a system with known transiting debris, to test our ability to recover similar systems. We find that the 3{\sigma} maximum occurrence rate of WD 1145+017-like transiting objects is < 0.5%.Comment: 17 pages, 13 figure

The star-formation history of the universe - an infrared perspective

A simple and versatile parameterized approach to the star formation history allows a quantitative investigation of the constraints from far infrared and submillimetre counts and background intensity measurements. The models include four spectral components: infrared cirrus (emission from interstellar dust), an M82-like starburst, an Arp220-like starburst and an AGN dust torus. The 60 $\mu$m luminosity function is determined for each chosen rate of evolution using the PSCz redshift data for 15000 galaxies. The proportions of each spectral type as a function of 60 $\mu$m luminosity are chosen for consistency with IRAS and SCUBA colour-luminosity relations, and with the fraction of AGN as a function of luminosity found in 12 $\mu$m samples. The luminosity function for each component at any wavelength can then be calculated from the assumed spectral energy distributions. With assumptions about the optical seds corresponding to each component and, for the AGN component, the optical and near infrared counts can be accurately modelled. A good fit to the observed counts at 0.44, 2.2, 15, 60, 90, 175 and 850 $\mu$m can be found with pure luminosity evolution in all 3 cosmological models investigated: $\Omega_o$ = 1, $\Omega_o$ = 0.3 ($\Lambda$ = 0), and $\Omega_o$ = 0.3, $\Lambda$ = 0.7. All 3 models also give an acceptable fit to the integrated background spectrum. Selected predictions of the models, for example redshift distributions for each component at selected wavelengths and fluxes, are shown. The total mass-density of stars generated is consistent with that observed, in all 3 cosmological models.Comment: 20 pages, 25 figures. Accepted for publication in ApJ. Full details of models can be found at http://astro.ic.ac.uk/~mrr/countmodel

The Redshift of GRB 970508

GRB 970508 is the second gamma-ray burst (GRB) for which an optical afterglow has been detected. It is the first GRB for which a distance scale has been determined: absorption and emission features in spectra of the optical afterglow place GRB 970508 at a redshift of z >= 0.835 (Metzger et al. 1997a, 1997b). The lack of a Lyman-alpha forest in these spectra further constrains this redshift to be less than approximately 2.3. I show that the spectrum of the optical afterglow of GRB 970508, once corrected for Galactic absorption, is inconsistent with the relativistic blast-wave model unless a second, redshifted source of extinction is introduced. This second source of extinction may be the yet unobserved host galaxy. I determine its redshift to be z = 1.09^{+0.14}_{-0.41}, which is consistent with the observed redshift of z = 0.835. Redshifts greater than z = 1.40 are ruled out at the 3 sigma confidence level.Comment: Accepted to The Astrophysical Journal (Letters), 10 pages, LaTe

Starburst and AGN activity in ultraluminous infrared galaxies

(Abridged) We examine the power source of 41 local Ultraluminous Infrared Galaxies using archival infrared and optical photometry. We fit the observed Spectral Energy Distributions (SEDs) with starburst and AGN components; each component being drawn from a family of templates. We find all of the sample require a starburst, whereas only half require an AGN. In 90% of the sample the starburst provides over half the IR emission, with a mean fractional luminosity of 82%. When combined with other galaxy samples we find that starburst and AGN luminosities correlate over 6 decades in IR luminosity suggesting that a common factor governs both luminosities, plausibly the gas masses in the nuclear regions. We find that the mid-IR 7.7 micron line-continuum ratio is no indication of the starburst luminosity, or the fractional AGN luminosity, and therefore that this ratio is not a reliable diagnostic of the power source in ULIRGs. We propose that the scatter in the radio-IR correlation in ULIRGs is due to a skewed starburst IMF and/or relic relativistic electrons from a previous starburst, rather than contamination from an obscured AGN. We show that most ULIRGs undergo multiple starbursts during their lifetime, and by inference that mergers between more than two galaxies may be common amongst ULIRGs. Our results support the evolutionary model for ULIRGs proposed by Farrah et al 2001, where they can follow many different evolutionary paths of starburst and AGN activity in transforming merging spiral galaxies into elliptical galaxies, but that most do not go through an optical QSO phase. The lower level of AGN activity in our local sample than in z~1 HLIRGs implies that the two samples are distinct populations. We postulate that different galaxy formation processes at high-z are responsible for this difference.Comment: 24 pages, 8 figures. Accepted for publication in MNRA

Faint Radio Sources and Star Formation History

Faint extragalactic radio sources provide important information about the global history of star formation. Sensitive radio observations of the Hubble Deep Field and other fields have found that sub-mJy radio sources are predominantly associated with star formation activity rather than AGN. Radio observations of star forming galaxies have the advantage of being independent of extinction by dust. We use the FIR-radio correlation to compare the radio and FIR backgrounds, and make several conclusions about the star forming galaxies producing the FIR background. We then use the redshift distribution of faint radio sources to determine the evolution of the radio luminosity function, and thus estimate the star formation density as a function of redshift.Comment: 12 pages, 9 figures, latex using texas.sty, to appear in the CD-ROM Proceedings of the 19th Texas Symposium on Relativistic Astrophysics and Cosmology, held in Paris, France, Dec. 14-18, 1998. Eds.: J. Paul, T. Montmerle, and E. Aubourg (CEA Saclay). No changes to paper, just updated publication info in this commen

Fitting optical source counts with an infrared-defined model: insights into galaxy evolution

In order to improve the fit to the optical counts, the multi-wavelength source count model of Rowan-Robinson (2001) is improved in two ways. First, density evolution is incorporated to correspond to the mergers in hierarchical models of galaxy formation. Secondly, the evolution of dust opacity is taken into account. The luminosity and density evolution parameters and the characteristic opacity of the local universe are allowed to vary freely and the best overall fit to the B band is determined. This parameter set is then translated to other wavebands, showing significant improvement of goodness of fit in the optical, UV and near-IR bands, whilst retaining the goodness of fit in the far-IR and submillimetre bands of the original model. A significant fraction of the improvement of fit is shown to be due to the introduction of optical depth history. The mean value of Av in the local universe for the best-fitting model is 0.4. The goodness of fit to available star formation history data is not significantly altered and the form of the evolution of the luminosity function shown to be consistent with a semi-analytic, forward evolution model.Comment: 12 pages (TeX file), 16 figures(eps format), submitted to MNRAS ; changed in accordance with referee's repor