Unresolved emission and ionized gas in the bulge of M31

We study the origin of unresolved X-ray emission from the bulge of M31 based on archival Chandra and XMM–Newton observations. We demonstrate that three different components are present. (i) Broad-band emission from a large number of faint sources – mainly accreting white dwarfs and active binaries, associated with the old stellar population, similar to the Galactic ridge X-ray emission of the Milky Way. The X-ray to K -band luminosity ratios are compatible with those for the Milky Way and for M32; in the 2–10 keV band, the ratio is  (3.6 ± 0.2) × 10²⁷ erg s⁻¹ L⁻¹_⊙  . (ii) Soft emission from ionized gas with a temperature of about ∼300 eV and a mass of  ∼2 × 10⁶ M_⊙  . The gas distribution is significantly extended along the minor axis of the galaxy, suggesting that it may be outflowing in the direction perpendicular to the galactic disc. The mass and energy supply from evolved stars and Type Ia supernovae is sufficient to sustain the outflow. We also detect a shadow cast on the gas emission by spiral arms and the 10-kpc star-forming ring, confirming significant extent of the gas in the 'vertical' direction. (iii) Hard extended emission from spiral arms, most likely associated with young stellar objects and young stars located in the star-forming regions. The   L _X/SFR  (star formation rate) ratio equals  ∼9 × 10³⁸ (erg s⁻¹)(M_⊙ yr⁻¹)⁻¹  , which is about ∼1/3 of the high-mass X-ray binary contribution, determined earlier from Chandra observations of other nearby galaxies.     

The Jeans mass and the origin of the knee in the IMF

We use numerical simulations of the fragmentation of a  1000 M_⊙  molecular cloud and the formation of a stellar cluster to study how the initial conditions for star formation affect the resulting initial mass function (IMF). In particular, we are interested in the relation between the thermal Jeans mass in a cloud and the knee of the IMF, i.e. the mass separating the region with a flat IMF slope from that typified by a steeper, Salpeter-like, slope. In three isothermal simulations with   M _Jeans= 1, 2  and  5 M_⊙  , the number of stars formed, at comparable dynamical times, scales roughly with the number of initial Jeans masses in the cloud. The mean stellar mass also increases (though less than linearly) with the initial Jeans mass in the cloud. It is found that the IMF in each case displays a prominent knee, located roughly at the mass scale of the initial Jeans mass. Thus clouds with higher initial Jeans masses produce IMFs which are shallow to higher masses. This implies that a universal IMF requires a physical mechanism that sets the Jeans mass to be near  1 M_⊙  . Simulations including a barotropic equation of state as suggested by Larson, with cooling at low densities followed by gentle heating at higher densities, are able to produce realistic IMFs with the knee located at  ≈1 M_⊙  , even with an initial   M _Jeans= 5 M_⊙  . We therefore suggest that the observed universality of the IMF in the local Universe does not require any fine tuning of the initial conditions in star forming clouds but is instead imprinted by details of the cooling physics of the collapsing gas.     

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.     

Quantifying the coexistence of massive black holes and dense nuclear star clusters

In large spheroidal stellar systems, such as elliptical galaxies, one invariably finds a  10⁶–10⁹ M_⊙  supermassive black hole at their centre. In contrast, within dwarf elliptical galaxies one predominantly observes a  10⁵–10⁷ M_⊙  nuclear star cluster. To date, few galaxies have been found with both types of nuclei coexisting and even less have had the masses determined for both central components. Here, we identify one dozen galaxies housing nuclear star clusters and supermassive black holes whose masses have been measured. This doubles the known number of such hermaphrodite nuclei – which are expected to be fruitful sources of gravitational radiation. Over the host spheroid (stellar) mass range  10⁸–10¹¹ M_⊙  , we find that a galaxy's nucleus-to-spheroid (baryon) mass ratio is not a constant value but decreases from a few per cent to ∼0.3 per cent such that  log[( M _BH+ M _NC)/ M _sph]=−(0.39 ± 0.07) log[ M _sph/10¹⁰ M_⊙]− (2.18 ± 0.07)  . Once dry merging commences and the nuclear star clusters disappear, this ratio is expected to become a constant value.
As a byproduct of our investigation, we have found that the projected flux from resolved nuclear star clusters is well approximated with Sérsic functions having a range of indices from ∼0.5 to ∼3, the latter index describing the Milky Way's nuclear star cluster.     

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.     

Cosmological origin of the lowest metallicity halo stars

We explore the predictions of the standard hierarchical clustering scenario of galaxy formation, regarding the numbers and metallicities of PopIII stars that are likely to be found within our Galaxy today. By PopIII we refer to stars formed at large redshift ( z >4), with low metallicities ([ Z /Z_⊙]<−2.5) and in small systems (total mass ≲ 2×10⁸ M_⊙) that are extremely sensitive to stellar feedback, and which through a prescribed merging history end up becoming part of the Milky Way today. An analytic, extended Press–Schechter formalism is used to obtain the mass functions of haloes which will host PopIII stars at a given redshift, and which will end up in Milky Way sized systems today. Each of these is modelled as a mini-galaxy, with a detailed treatment of the dark halo structure, angular momentum distribution, final gas temperature and disc instabilities, all of which determine the fraction of the baryons that are subject to star formation. The use of new primordial metallicity stellar evolutionary models allows us to trace the history of the stars formed, and give accurate estimates of their expected numbers today and their location in L /L_⊙ versus T /K Hertzsprung–Russell (HR) diagrams. A first comparison with observational data suggests that the initial mass function (IMF) of the first stars was increasingly high-mass weighted towards high redshifts, levelling off at z ≳9 at a characteristic stellar mass scale m _s=10–15 M_⊙.     

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.     

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.     

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.     

The census of nuclear activity of late-type galaxies in the Virgo cluster

The first spectroscopic census of active galactic nuclei (AGNs) associated with late-type galaxies in the Virgo cluster was carried out by observing 213 out of a complete set of 237 galaxies more massive than   M _dyn > 10^8.5 M_⊙  . Among them, 77 are classified as AGNs [including 21 transition objects, 47 low-ionization nuclear emission regions (LINERs) and nine Seyferts] and comprise 32 per cent of the late-type galaxies in Virgo. Due to spectroscopic incompleteness, at most 21 AGNs are missed in the survey, so that the fraction would increase up to 41 per cent. Using corollary near-infrared observations that enable us to estimate galaxy dynamical masses, it is found that AGNs are hosted exclusively in massive galaxies, i.e.   M _dyn≳ 10¹⁰ M_⊙  . Their frequency increases steeply with the dynamical mass from zero at   M _dyn≈ 10^9.5 M_⊙  to virtually 1 at   M _dyn > 10^11.5 M_⊙  . These frequencies are consistent with those of low-luminosity AGNs found in the general field by the Sloan Digital Sky Survey. Massive galaxies that harbour AGNs commonly show conspicuous r -band star-like nuclear enhancements. Conversely, they often, but not necessarily, contain massive bulges. A few well-known AGNs (e.g. M61, M100, NGC 4535) are found in massive Sc galaxies with little or no bulge. The AGN fraction seems to be only marginally sensitive to galaxy environment. We infer the black hole masses using the known scaling relations of quiescent black holes. No black holes lighter than  ∼10⁶ M_⊙  are found active in our sample.     

Angular momentum and clustering properties of early dark matter haloes

In this paper, we study the angular momentum properties of simulated dark matter haloes at high redshifts that likely host the first stars in the Universe. Calculating the spin distributions of these  10⁶– 10⁷ M_⊙  haloes in redshift slices from   z = 15  to 6, we find that they are well fit by a lognormal distribution as is found for lower redshift and more massive haloes in earlier work. We find that both the mean value of the spin and dispersion are largely unchanged with redshift for all haloes. Our key result is that subsamples of low- and high-spin, 10⁶ and  10⁷ M_⊙  , haloes show difference in clustering strength. In both mass bins, higher spin haloes are more strongly clustered in concordance with a tidal torquing picture for the growth of angular momentum in dark matter haloes in the cold dark matter paradigm.     

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⁻¹  .     

A ultraluminous X-ray source associated with a cloud collision in M 99

The Sc galaxy M 99 in the Virgo Cluster has been strongly affected by tidal interactions and recent close encounters, responsible for an asymmetric spiral pattern and a high star formation rate. Our XMM–Newton study shows that the inner disc is dominated by hot plasma at kT ≈ 0.30 keV, with a total X-ray luminosity of ≈10⁴¹ erg s⁻¹ in the 0.3–12 keV band. At the outskirts of the galaxy, away from the main star-forming regions, there is an ultraluminous X-ray source (ULX) with an X-ray luminosity of ≈2 × 10⁴⁰ erg s⁻¹ and a hard spectrum well fitted by a power law of photon index Γ≈ 1.7. This source is close to the location where a massive H  i cloud appears to be falling on to the M 99 disc at a relative speed of >100 km s⁻¹. We suggest that there may be a direct physical link between fast cloud collisions and the formation of bright ULXs, which may be powered by accreting black holes with masses ∼100 M_⊙. External collisions may trigger large-scale dynamical collapses of protoclusters, leading to the formation of very massive (≳200 M_⊙) stellar progenitors; we argue that such stars may later collapse into massive black holes if their metal abundance is sufficiently low.     

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.     

Masses for the Local Group and the Milky Way

We use the very large Millennium Simulation of the concordance Λ cold dark matter cosmogony to calibrate the bias and error distribution of Timing Argument estimators of the masses of the Local Group and of the Milky Way. From a large number of isolated spiral–spiral pairs similar to the Milky Way/Andromeda system, we find the interquartile range of the ratio of timing mass to true mass to be a factor of 1.8, while the 5 and 95 per cent points of the distribution of this ratio are separated by a factor of 5.7. Here, we define true mass as the sum of the 'virial' masses, M ₂₀₀, of the two dominant galaxies. For present best values of the distance and approach velocity of Andromeda, this leads to a median likelihood estimate of the true mass of the Local Group of  5.27 × 10¹² M_⊙  or  log  M _LG/M_⊙= 12.72  , with an interquartile range of [12.58, 12.83] and a 5–95 per cent range of [12.26, 13.01]. Thus, a 95 per cent lower confidence limit on the true mass of the Local Group is  1.81 × 10¹² M_⊙  . A timing estimate of the Milky Way's mass based on the large recession velocity observed for the distant satellite Leo I works equally well, although with larger systematic uncertainties. It gives an estimated virial mass for the Milky Way of  2.43 × 10¹² M_⊙  with a 95 per cent lower confidence limit of  0.80 × 10¹² M_⊙  .     

Magnetically collimated jets with high mass flux

We present HST /WFPC2 observations of UGC 4483, an irregular galaxy in the M81/NGC 2403 complex. Stellar photometry was carried out with HSTphot, and is complete to V ≃26.0 and I ≃24.7. We measure the red giant branch tip at I =23.56±0.10, and calculate a distance modulus of μ ₀=27.53±0.12 (corresponding to a distance of 3.2±0.2 Mpc), placing UGC 4483 within the NGC 2403 subgroup. We were able to measure properties of a previously known young star cluster in UGC 4483, finding integrated magnitudes of V =18.66±0.21 and I =18.54±0.10 for the stellar contribution (integrated light minus H α and [O  iii ] contribution), corresponding to an age of ∼10–15 Myr and an initial mass of ∼10⁴ M_⊙. This is consistent with the properties of the cluster's brightest stars, which were resolved in the data for the first time. Finally, a numerical analysis of the galaxy's stellar content yields a roughly constant star formation rate of 1.3×10⁻³ M_⊙ yr⁻¹ and mean metallicity of [Fe/H]=−1.3 dex from 15 Gyr ago to the present.     

The tidal disruption rate in dense galactic cusps containing a supermassive binary black hole

We consider the problem of tidal disruption of stars in the centre of a galaxy containing a supermassive binary black hole with unequal masses. We assume that over the separation distance between the black holes, the gravitational potential is dominated by the more massive primary black hole. Also, we assume that the number density of stars is concentric with the primary black hole and has a power-law cusp. We show that the bulk of stars with a small angular-momentum component normal to the black hole binary orbit can reach a small value of total angular momentum through secular evolution in the gravitational field of the binary, and hence they can be tidally disrupted by the larger black hole. This effect is analogous to the so-called Kozai effect well known in celestial mechanics. We develop an analytical theory for the secular evolution of the stellar orbits and calculate the rate of tidal disruption. We compare our analytical theory with a simple numerical model and find very good agreement.
Our results show that for a primary black hole mass of  ∼10⁶–10⁷ M_⊙  , the black hole mass-ratio   q > 10⁻²  , cusp size ∼1 pc, the tidal disruption rate can be as large as  ∼10⁻²–1 M_⊙ yr⁻¹  . This is at least 10²–10⁴ times larger than estimated for the case of a single supermassive black hole. The duration of the phase of enhanced tidal disruption is determined by the dynamical-friction time-scale, and it is rather short: ∼10⁵ yr. The dependence of the tidal disruption rate on the mass ratio, and on the size of the cusp, is also discussed.     

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.