A top-heavy stellar initial mass function in starbursts as an explanation for the high mass-to-light ratios of ultra-compact dwarf galaxies

It has been recently shown that the dynamical V -band mass-to-light ratios of compact stellar systems with masses from 10⁶ to  10⁸ M_⊙  are not consistent with the predictions from simple stellar population models. Top-heavy stellar initial mass functions (IMFs) in these so-called ultra-compact dwarf galaxies (UCDs) offer an attractive explanation for this finding, the stellar remnants and retained stellar envelopes providing the unseen mass. We therefore construct a model which quantifies by how much the IMFs of UCDs would have to deviate in the intermediate- and high-mass range from the canonical IMF in order to account for the enhanced   M / L_V   ratio of the UCDs. The deduced high-mass IMF in the UCDs depends on the age of the UCDs and the number of faint products of stellar evolution retained by them. Assuming that the IMF in the UCDs is a three-part power law equal to the canonical IMF in the low-mass range and taking 20 per cent as a plausible choice for the fraction of the remnants of high-mass stars retained by UCDs, the model suggests the exponent of the high-mass IMF to be ≈1.6 if the UCDs are  13 Gyr  old (i.e. almost as old as the Universe) or ≈1.0 if the UCDs are  7 Gyr  old, in contrast to 2.3 for the Salpeter–Massey IMF. If the IMF was as top heavy as suggested here, the stability of the UCDs might have been threatened by heavy mass loss induced by the radiation and evolution of massive stars. The central densities of UCDs must have been in the range  10⁶ to 10⁷ M_⊙ pc⁻³  when they formed with star formation rates of  10 to 100 M_⊙ yr⁻¹  .     

The initial mass function of the rich young cluster NGC 1818 in the Large Magellanic Cloud

We use deep Hubble Space Telescope photometry of the rich, young (∼20- to 45-Myr old) star cluster NGC 1818 in the Large Magellanic Cloud to derive its stellar mass function (MF) down to  ∼0.15 M_⊙  . This represents the deepest robust MF thus far obtained for a stellar system in an extragalactic, low-metallicity  ([Fe/H]≃−0.4 dex)  environment. Combining our results with the published MF for masses above  1.0 M_⊙  , we obtain a complete present-day MF. This is a good representation of the cluster's initial MF (IMF), particularly at low masses, because our observations are centred on the cluster's uncrowded half-mass radius. Therefore, stellar and dynamical evolution of the cluster will not have affected the low-mass stars significantly. The NGC 1818 IMF is well described by both a lognormal and a broken power-law distribution with slopes of  Γ= 0.46 ± 0.10  and  Γ≃−1.35  (Salpeter-like) for masses in the range from 0.15 to  0.8 M_⊙  and greater than  0.8 M_⊙  , respectively. Within the uncertainties, the NGC 1818 IMF is fully consistent with both the Kroupa solar neighbourhood and the Chabrier lognormal mass distributions.     

A SCUBA survey of Orion – the low-mass end of the core mass function

We have re-analysed all of the Submillimetre Common User Bolometer Array (SCUBA) archive data of the Orion star-forming regions. We have put together all of the data taken at different times by different groups. Consequently, we have constructed the deepest submillimetre maps of these regions ever made. There are four regions that have been mapped: Orion A North and South, and Orion B North and South. We find that two of the regions, Orion A North and Orion B North, have deeper sensitivity and completeness limits, and contain a larger number of sources, so we concentrate on these two. We compare the data with archive data from the Spitzer Space Telescope to determine whether or not a core detected in the submillimetre is pre-stellar in nature. We extract all of the pre-stellar cores from the data and make a histogram of the core masses. This can be compared to the stellar initial mass function (IMF). We find the high-mass core mass function (CMF) follows a roughly Salpeter-like slope, just like the IMF, as seen in previous work. Our deeper maps allow us to see that the CMF turns over at,  ∼1.3 M_⊙  about a factor of 4 higher than our completeness limit. This turnover has never previously been observed, and is only visible here due to our much deeper maps. It mimics the turnover seen in the stellar IMF at  ∼0.1 M_⊙  . The low-mass side of the CMF is a power law with an exponent of, 0.35 ± 0.2 which is consistent with the low-mass slope of the young cluster IMF of 0.3 ± 0.1. This shows that the CMF continues to mimic the shape of the IMF all the way down to the lower completeness limit of these data at  ∼0.3 M_⊙  .     

Optical and submillimetre observations of Bok globules – tracing the magnetic field from low to high density

Although the stellar initial mass function (IMF) has only been directly determined in star clusters, it has been manifoldly applied on galaxy-wide scales. But taking the clustered nature of star formation into account the galaxy-wide IMF is constructed by adding all IMFs of all young star clusters leading to an integrated galactic initial mass function (IGIMF). The IGIMF is top-light compared to the canonical IMF in star clusters and steepens with decreasing total star formation rate (SFR). This discrepancy is marginal for large disc galaxies but becomes significant for Small Magellanic Cloud type galaxies and less massive ones. We here construct IGIMF-based relations between the total far- and near-ultraviolet luminosities of galaxies and the underlying SFR. We make the prediction that the Hα luminosity of star-forming dwarf galaxies decreases faster with decreasing SFR than the ultraviolet (UV) luminosity. This turn-down of the Hα/UV-flux ratio should be evident below total SFRs of  10⁻² M_⊙ yr⁻¹  .     

On the mass of dense star clusters in starburst galaxies from spectrophotometry

The mass of unresolved young star clusters derived from spectrophotometric data may well be off by a factor of 2 or more once the migration of massive stars driven by mass segregation is accounted for. We quantify this effect for a large set of cluster parameters, including variations in the stellar initial mass function (IMF), the intrinsic cluster mass, and mean mass density. Gas-dynamical models coupled with the Cambridge stellar evolution tracks allow us to derive a scheme to recover the real cluster mass given measured half-light radius, one-dimensional velocity dispersion and age. We monitor the evolution with time of the ratio of real to apparent mass through the parameter η. When we compute η for rich star clusters, we find non-monotonic evolution in time when the IMF stretches beyond a critical cut-off mass of  25.5 M_⊙  . We also monitor the rise of colour gradients between the inner and outer volume of clusters: we find trends in time of the stellar IMF power indices overlapping well with those derived for the Large Magellanic Cloud cluster NGC 1818 at an age of 30 Myr. We argue that the core region of massive Antennae clusters should have suffered from much segregation despite their low ages. We apply these results to a cluster mass function, and find that the peak of the mass distribution would appear to observers shifted to lower masses by as much as 0.2 dex. The star formation rate derived for the cluster population is then underestimated by from 20 to 50 per cent.     

Baryonic conversion tree: the global assembly of stars and dark matter in galaxies from the Sloan Digital Sky Survey

Using the spectroscopic sample of the Sloan Digital Sky Survey Data Release 1 (SDSS DR1), we measure how gas was transformed into stars as a function of time and stellar mass: the baryonic conversion tree (BCT). There is a clear correlation between early star formation activity and present-day stellar mass: the more massive galaxies have formed approximately 80 per cent of their stars at   z > 1  , while for the less massive ones the value is only approximately 20 per cent. By comparing the BCT with the dark matter merger tree, we find indications that star formation efficiency at   z > 1  had to be approximately a factor of two higher than today (∼10 per cent) in galaxies with present-day stellar mass larger than  2 × 10¹¹ M_⊙  , if this early star formation occurred in the main progenitor. Therefore, the λ cold dark matter (LCDM) paradigm can accommodate a large number of red objects. On the other hand, in galaxies with present-day stellar mass less than  10¹¹ M_⊙  , efficient star formation seems to have been triggered at   z ∼ 0.2  . We show that there is a characteristic mass  ( M _*∼ 10¹⁰ M_⊙)  for feedback efficiency (or lack of star formation). For galaxies with masses lower than this, feedback (or star formation suppression) is very efficient while for higher masses it is not. The BCT, determined here for the first time, should be an important observable with which to confront theoretical models of galaxy formation.     

Predictions of the extent of self-enrichment in oxygen of giant metal-poor H ii regions

In general, H  ii regions do not show clear signs of self-enrichment in products from massive stars  ( M ≥ 8 M_⊙)  . In order to explore why I modelled the contamination with Wolf–Rayet star ejecta of metal-poor  ( Z = 0.001)  H  ii regions, ionized either by a  10⁶ M_⊙  cluster of coeval stars (cluster 1) or by a cluster resulting from continuous star formation at a rate of  1 M_⊙ yr⁻¹  (cluster 2). The clusters have   Z = 0.001  and a Salpeter initial mass function from 0.1 to  120 M_⊙  . Independent one-dimensional constant density simulations of the emission-line spectra of unenriched H  ii regions were computed at the discrete ages 1, 2, 3, 4 and 5 Myr, with the photoionization code cloudy , using as input, radiative and mechanical stellar feedbacks predicted by the evolutionary synthesis code starburst99 . Each H  ii region was placed at the outer radius of the adiabatically expanding superbubble of Mac Low & McCray. For models with thermal and ionization balance time-scales of less than 1 Myr, and with oxygen emission-line ratios in agreement with observations, the volume of the superbubble and the H  ii region was uniformly and instantaneously polluted with stellar ejecta predicted by starburst99 . I obtained a maximum oxygen abundance enhancement of 0.025 dex, with cluster 1, at 4 Myr. It would be unobservable.     

Single and binary evolution of Population III stars and their supernova explosions

We present stellar evolution calculations for Population III stars for both single- and binary-star evolutions. Our models include 10- and  16.5-M_⊙  single stars and a  10-M_⊙  model star that undergoes an episode of accretion resulting in a final mass of  16.1 M_⊙  . For comparison, we present the evolution of a solar heavy element abundance model. We use the structure from late-stage evolution models to calculate simulated supernova light curves. Light curve comparisons are made between accretion and non-accretion progenitor models, and models for single-star evolution of comparable masses. Where possible, we make comparisons to previous works. Similar investigations have been carried out, but primarily for solar or near-solar heavy metal abundance stars and not including both the evolution and the supernova explosions in one work.     

Modelling the binary progenitor of Supernova 1993J

We have developed a detailed stellar evolution code capable of following the simultaneous evolution of both stars in a binary system, together with their orbital properties. To demonstrate the capabilities of the code, we investigate potential progenitors for the Type IIb Supernova 1993J, which is believed to have been an interacting binary system prior to its primary exploding. We use our detailed binary stellar evolution code to model this system to determine the possible range of primary and secondary masses that could have produced the observed characteristics of this system, with particular reference to the secondary. Using the luminosities and temperatures for both stars (as determined by Maund et al.) and the remaining mass of the hydrogen envelope of the primary at the time of explosion, we find that if mass transfer is 100 per cent efficient, the observations can be reproduced by a system consisting of a  15 M_⊙  primary and a  14 M_⊙  secondary in an orbit with an initial period of 2100 days. With a mass transfer efficiency of 50 per cent, a more massive system consisting of a  17 M_⊙  primary and a  16 M_⊙  secondary in an initial orbit of 2360 days is needed. We also investigate some of the uncertainties in the evolution, including the effects of tidal interaction, convective overshooting and thermohaline mixing.     

Active galactic nuclei jet mass loading and truncation by stellar winds

Active galactic nuclei can produce extremely powerful jets. While tightly collimated, the scale of these jets and the stellar density at galactic centres implies that there will be many jet/star interactions, which can mass load the jet through stellar winds. Previous work employed modest wind mass outflow rates, but this does not apply when mass loading is provided by a small number of high mass-loss stars. We construct a framework for jet mass loading by stellar winds for a broader spectrum of wind mass-loss rates than has previously been considered. Given the observed stellar mass distributions in galactic centres, we find that even highly efficient (0.1 Eddington luminosity) jets from supermassive black holes of masses M _BH≲ 10⁴ M_⊙ are rapidly mass loaded and quenched by stellar winds. For  10⁴ M_⊙ < M _BH < 10⁸ M_⊙  , the quenching length of highly efficient jets is independent of the jet's mechanical luminosity. Stellar wind mass loading is unable to quench efficient jets from more massive engines, but can account for the observed truncation of the inefficient M87 jet, and implies a baryon-dominated composition on scales ≳2 kpc therein even if the jet is initially pair plasma dominated.     

New constraints on turbulent transport in accretion discs

Many objects studied in astronomy follow a power-law distribution function (DF), for example the masses of stars or star clusters. A still used method by which such data is analysed is to generate a histogram and fit a straight line to it. The parameters obtained in this way can be severely biased, and the properties of the underlying DF, such as its shape or a possible upper limit, are difficult to extract. In this work, we review techniques available in the literature and present newly developed (effectively) bias-free estimators for the exponent and the upper limit. Furthermore, we discuss various graphical representations of the data and powerful goodness-of-fit tests to assess the validity of a power law for describing the distribution of data. As an example, we apply the presented methods to the data set of massive stars in R136 and the young star clusters in the Large Magellanic Cloud. For R136 we confirm the result of Koen of a truncated power law with a bias-free estimate for the exponent of  2.20 ± 0.78/2.87 ± 0.98  (where the Salpeter–Massey value is 2.35) and for the upper limit of  143 ± 9/163 ± 9 M_⊙  , depending on the stellar models used. The star clusters in the Large Magellanic Cloud (with ages up to  10^7.5 yr  ) follow a truncated power-law distribution with exponent  1.62 ± 0.06  and upper limit  68 ± 12 × 10³ M_⊙  . Using the graphical data representation, a significant change in the form of the mass function below  10^2.5 M_⊙  can be detected, which is likely caused by incompleteness in the data.     

An iterative filter to reconstruct planetary transit signals in the presence of stellar variability

We report the identification, from a photometric, astrometric and spectroscopic study, of a massive white dwarf member of the nearby, approximately solar metallicity, Coma Berenices open star cluster (Melotte 111). We find the optical to near-infrared energy distribution of WD 1216+260 to be entirely consistent with that of an isolated DA and determine the effective temperature and surface gravity of this object to be   T _eff= 15 739⁺¹⁹⁷₋₁₉₆ K  and  log  g = 8.46^+0.03_−0.02  . We set tight limits on the mass of a putative cool companion,   M ≳ 0.036 M_⊙  (spatially unresolved) and   M ≳ 0.034 M_⊙  (spatially resolved and   a ≲ 2500 au  ). Based on the predictions of CO core, thick H layer evolutionary models we determine the mass and cooling time of WD 1216+260 to be   M _WD= 0.90 ± 0.04 M_⊙  and  τ_cool= 363⁺⁴⁶₋₄₁ Myr  , respectively. For an adopted cluster age of  τ= 500 ± 100 Myr  we infer the mass of its progenitor star to be   M _init= 4.77^+5.37_−0.97 M_⊙  . We briefly discuss this result in the context of the form of the stellar initial mass–final mass relation.     

The low-mass initial mass function in the young cluster NGC 6611

NGC 6611 is the massive young cluster (2–3 Myr) that ionizes the Eagle Nebula. We present very deep photometric observations of the central region of NGC 6611 obtained with the Hubble Space Telescope and the following filters: ACS/WFC F775W and F850LP and NIC2 F110W and F160W, loosely equivalent to ground-based IZJH filters. This survey reaches down to   I ∼ 26 mag  . We construct the initial mass function (IMF) from  ∼1.5 M_⊙  well into the brown dwarf regime (down to  ∼0.02 M_⊙  ). We have detected 30–35 brown dwarf candidates in this sample. The low-mass IMF is combined with a higher-mass IMF constructed from the ground-based catalogue from Oliveira et al. We compare the final IMF with those of well-studied star-forming regions: we find that the IMF of NGC 6611 more closely resembles that of the low-mass star-forming region in Taurus than that of the more massive Orion Nebula Cluster. We conclude that there seems to be no severe environmental effect in the IMF due to the proximity of the massive stars in NGC 6611.     

The host galaxy of a narrow-line Seyfert 1 galaxy, RE J1034+396, with X-ray quasi-periodic oscillations

Using simple stellar population synthesis, we model the bulge stellar contribution in the optical spectrum of a narrow-line Seyfert 1 galaxy, RE J1034+396. We find that its bulge stellar velocity dispersion is  67.7 ± 8 km s⁻¹  . The supermassive black hole (SMBH) mass is about  (1–4) × 10⁶ M_⊙  if it follows the well-known   M _BH–σ_*  relation found in quiescent galaxies. We also derive the SMBH mass from the Hβ second moment, which is consistent with that from its bulge stellar velocity dispersion. The SMBH mass of (1–4)  × 10⁶ M_⊙  implies that the X-ray quasi-periodic oscillation (QPO) of RE J1034+396 can be scaled to a high-frequency QPO at 27–108 Hz found in Galactic black hole binaries with a  10-M_⊙  black hole. With the mass distribution in different age stellar populations, we find that the mean specific star formation rate (SSFR) over the past 0.1 Gyr is  0.0163 ± 0.0011  Gyr⁻¹, the stellar mass in the logarithm is  10.155 ± 0.06  in units of solar mass and the current star formation rate is  0.23 ± 0.016 M_⊙ yr⁻¹  . For RE J1034+396, there is no relation between the Eddington ratio and the SSFR as suggested by Chen et al., despite a larger scatter in their relation. We also suggest that about 7.0 per cent of the total Hα luminosity and 50 per cent of the total [O  ii ] luminosity come from the star formation process.     

Spectroscopic and photometric observations of the eclipsing star UV Leo

High-resolution spectroscopic observations around the Hα line and BVRI photometry of the eclipsing short-period RS CVn star UV Leo are presented. The simultaneous light-curve solution and radial velocity-curve solution led to the following values of the global parameters of the binary: temperatures   T ₁= 6000 ± 100 K  and   T ₂= 5970 ± 20 K  ; masses   M ₁= 0.976 ± 0.067 M_⊙  and   M ₂= 0.931 ± 0.052 M_⊙  ; separation   a = 3.716 ± 0.048 R_⊙  ; orbital inclination     ; radii   R ₁= 1.115 ± 0.052 R_⊙  and   R ₂= 1.078 ± 0.051 R_⊙  ; equatorial velocities   V ₁= 98.8 ± 2.3 km s⁻¹  and   V ₂= 89.6 ± 2.7 km s⁻¹  . These results lead to the conclusion that the two components of UV Leo are slightly oversized for their masses and lie within the main-sequence band on the mass–radius diagram, close to the isochrone 9 × 10¹⁰ yr.     

On the first generation of stars

We argue that the first stars may have spanned the conventional mass range rather than be identified with the very massive objects  (∼100–10³ M_⊙)  favoured by numerical simulations. Specifically, we find that magnetic field generation processes acting in the first protostellar systems suffice to produce fields that exceed the threshold for magneto-rotational instability (MRI) to operate, and thereby allow the MRI dynamo to generate equipartition-amplitude magnetic fields on protostellar mass scales below  ∼50 M_⊙  . Such fields allow primordial star formation to occur at essentially any metallicity by regulating angular momentum transfer, fragmentation, accretion and feedback in much the same way as occurs in conventional molecular clouds.     

On the galaxy stellar mass function, the mass–metallicity relation and the implied baryonic mass function

A comparison between published field galaxy stellar mass functions (GSMFs) shows that the cosmic stellar mass density is in the range 4–8 per cent of the baryon density (assuming  Ω_b= 0.045  ). There remain significant sources of uncertainty for the dust correction and underlying stellar mass-to-light ratio even assuming a reasonable universal stellar initial mass function. We determine the   z < 0.05  GSMF using the New York University Value-Added Galaxy Catalog sample of 49 968 galaxies derived from the Sloan Digital Sky Survey and various estimates of stellar mass. The GSMF shows clear evidence for a low-mass upturn and is fitted with a double Schechter function that has  α₂≃−1.6  . At masses below  ∼10^8.5 M_⊙  , the GSMF may be significantly incomplete because of missing low-surface-brightness galaxies. One interpretation of the stellar mass–metallicity relation is that it is primarily caused by a lower fraction of available baryons converted to stars in low-mass galaxies. Using this principle, we determine a simple relationship between baryonic mass and stellar mass and present an 'implied baryonic mass function'. This function has a faint-end slope,  α₂≃−1.9  . Thus, we find evidence that the slope of the low-mass end of the galaxy mass function could plausibly be as steep as the halo mass function. We illustrate the relationship between halo baryonic mass function → galaxy baryonic mass function → GSMF. This demonstrates the requirement for peak galaxy formation efficiency at baryonic masses  ∼10¹¹ M_⊙  corresponding to a minimum in feedback effects. The baryonic-infall efficiency may have levelled off at lower masses.     

The abundance of brown dwarfs

The amount of mass contained in low-mass objects is investigated anew. Instead of using a mass–luminosity relation to convert a luminosity function to a mass function, I predict the mass–luminosity relation from assumed mass functions and the luminosity functions of Jahreiss & Wielen and Gould, Bahcall & Flynn. Comparison of the resulting mass–luminosity relations with data for binary stars constrains the permissible mass functions. If the mass function is assumed to be a power law, the best-fitting slope lies either side of the critical slope, α =−2, below which the mass in low-mass objects is divergent, depending on the luminosity function adopted. If these power-law mass functions are truncated at 0.001 M_⊙, the contribution to the local density from stars lies between 0.013 and 0.10 M_⊙ pc⁻³ depending on the mass at which the mass function is normalized and the adopted value of α . Recent dynamical estimates of the local mass density rule out stellar mass densities above ∼0.05 M_⊙ pc⁻³. Hence, power laws steeper than α =−2 that extend down to 0.001 M_⊙ are allowed only if one adopts an implausible normalization of the mass function. If the mass function is generalized from a power law to a low-order polynomial in log( M ), the mass in stars with M <0.1 M_⊙ is either negligible or strongly divergent, depending on the order of the polynomial adopted.     

Impact of cosmic rays on Population III star formation

We explore the implications of a possible cosmic-ray (CR) background generated during the first supernova explosions that end the brief lives of massive Population III stars. We show that such a CR background could have significantly influenced the cooling and collapse of primordial gas clouds in minihaloes around redshifts of   z ∼ 15–20  , provided the CR flux was sufficient to yield an ionization rate greater than about 10⁻¹⁹ s⁻¹ near the centre of the minihalo. The presence of CRs with energies  ≲10⁷  eV would indirectly enhance the molecular cooling in these regions, and we estimate that the resulting lower temperatures in these minihaloes would yield a characteristic stellar mass as low as  ∼10 M_⊙  . CRs have a less-pronounced effect on the cooling and collapse of primordial gas clouds inside more massive dark matter haloes with virial masses  ≳10⁸ M_⊙  at the later stages of cosmological structure formation around   z ∼ 10–15  . In these clouds, even without CR flux the molecular abundance is already sufficient to allow cooling to the floor set by the temperature of the cosmic microwave background.