Showing papers in "Monthly Notices of the Royal Astronomical Society: Letters in 2015"

Raising the bar: new constraints on the Hubble parameter with cosmic chronometers at z ∼ 2

Abstract: One of the most compelling tasks of modern cosmology is to constrain the expansion history of the Universe, since this measurement can give insights on the nature of dark energy and help to estimate cosmological parameters. In this letter are presented two new measurements of the Hubble parameter H(z) obtained with the cosmic chronometer method up to z ∼ 2. Taking advantage of near-infrared spectroscopy of the few very massive and passive galaxies observed at z > 1.4 available in literature, the differential evolution of this population is estimated and calibrated with different stellar population synthesis models to constrain H(z), including in the final error budget all possible sources of systematic uncertainties (star formation history, stellar metallicity, model dependences). This analysis is able to extend significantly the redshift range coverage with respect to present-day constraints, crossing for the first time the limit at z ∼ 1.75. The new H(z) data are used to estimate the gain in accuracy on cosmological parameters with respect to previous measurements in two cosmological models, finding a small but detectable improvement (∼5 per cent) in particular on Ω_M and w_0. Finally, a simulation of a Euclid-like survey has been performed to forecast the expected improvement with future data. The provided constraints have been obtained just with the cosmic chronometers approach, without any additional data, and the results show the high potentiality of this method to constrain the expansion history of the Universe at these redshifts.

On planet formation in HL Tau

Abstract: We explain the axisymmetric gaps seen in recent long-baseline observations of the HL Tau protoplanetary disc with the Atacama Large Millimetre/Submillimetre Array (ALMA) as being due to the different response of gas and dust to embedded planets in protoplanetary discs. We perform global, three dimensional dusty smoothed particle hydrodynamics calculations of multiple planets embedded in dust/gas discs which successfully reproduce most of the structures seen in the ALMA image. We find a best match to the observations using three embedded planets with masses of 0.2, 0.27 and 0.55 MJ in the three main gaps observed by ALMA, though there remain uncertainties in the exact planet masses from the disc model.

Formation of complex organic molecules in cold objects: the role of gas-phase reactions

Abstract: While astrochemical models are successful in reproducing many of the observed interstellar species, they have been struggling to explain the observed abundances of complex organic molecules. Current models tend to privilege grain surface over gas phase chemistry in their formation. One key assumption of those models is that radicals trapped in the grain mantles gain mobility and react on lukewarm (>30 K) dust grains. Thus, the recent detections of methyl formate (MF) and dimethyl ether (DME) in cold objects represent a challenge and may clarify the respective role of grain surface and gas phase chemistry. We propose here a new model to form DME and MF with gas phase reactions in cold environments, where DME is the precursor of MF via an efficient reaction overlooked by previous models. Furthermore, methoxy, a precursor of DME, is also synthetized in the gas phase from methanol, which is desorbed by a non-thermal process from the ices. Our new model reproduces fairy well the observations towards L1544. It also explains, in a natural way, the observed correlation between DME and MF. We conclude that gas phase reactions are major actors in the formation of MF, DME and methoxy in cold gas. This challenges the exclusive role of grain-surface chemistry and favours a combined grain-gas chemistry.

polychord: nested sampling for cosmology

Abstract: PolyChord is a novel nested sampling algorithm tailored for high dimensional parameter spaces. In addition, it can fully exploit a hierarchy of parameter speeds such as is found in CosmoMC and CAMB. It utilises slice sampling at each iteration to sample within the hard likelihood constraint of nested sampling. It can identify and evolve separate modes of a posterior semi-independently and is parallelised using openMPI. PolyChord is available for download at: this http URL

Gas-phase formation of the prebiotic molecule formamide: insights from new quantum computations

Abstract: New insights into the formation of interstellar formamide, a species of great relevance in prebiotic chemistry, are provided by electronic structure and kinetic calculations for the reaction NH2 + H2CO -> NH2CHO + H. Contrarily to what previously suggested, this reaction is essentially barrierless and can, therefore, occur under the low temperature conditions of interstellar objects thus providing a facile formation route of formamide. The rate coefficient parameters for the reaction channel leading to NH2CHO + H have been calculated to be A = 2.6x10^{-12} cm^3 s^{-1}, beta = -2.1 and gamma = 26.9 K in the range of temperatures 10-300 K. Including these new kinetic data in a refined astrochemical model, we show that the proposed mechanism can well reproduce the abundances of formamide observed in two very different interstellar objects: the cold envelope of the Sun-like protostar IRAS16293-2422 and the molecular shock L1157-B2. Therefore, the major conclusion of this Letter is that there is no need to invoke grain-surface chemistry to explain the presence of formamide provided that its precursors, NH2 and H2CO, are available in the gas-phase.

ALMA and Herschel reveal that X-ray-selected AGN and main-sequence galaxies have different star formation rate distributions

Abstract: Using deep Herschel and ALMA observations, we investigate the star formation rate (SFR) distributions of X-ray-selected active galactic nucleus (AGN) host galaxies at 0.5 < z < 1.5 and 1.5 < z < 4, comparing them to that of normal, star-forming (i.e. `main-sequence', or MS) galaxies. We find that 34-55 per cent of AGNs in our sample have SFRs at least a factor of 2 below that of the average MS galaxy, compared to ≈15 per cent of all MS galaxies, suggesting significantly different SFR distributions. Indeed, when both are modelled as lognormal distributions, the mass and redshift-normalized SFR distributions of X-ray AGNs are roughly twice as broad, and peak ≈0.4 dex lower, than that of MS galaxies. However, like MS galaxies, the normalized SFR distribution of AGNs in our sample appears not to evolve with redshift. Despite X-ray AGNs and MS galaxies having different SFR distributions, the linear-mean SFR of AGNs derived from our distributions is remarkably consistent with that of MS galaxies, and thus with previous results derived from stacked Herschel data. This apparent contradiction is due to the linear-mean SFR being biased by bright outliers, and thus does not necessarily represent a true characterization of the typical SFR of X-ray AGNs.

The dust mass in z > 6 normal star-forming galaxies

Abstract: We interpret recent ALMA observations of z > 6 normal star forming galaxies by means of a semi-numerical method, which couples the output of a cosmological hydrodynamical simulation with a chemical evolution model which accounts for the contribution to dust enrichment from supernovae, asymptotic giant branch stars and grain growth in the interstellar medium. We find that while stellar sources dominate the dust mass of small galaxies, the higher level of metal enrichment experienced by galaxies with Mstar > 10 9 M⊙ allows efficient grain growth, which provides the dominant contribution to the dust mass. Even assuming maximally efficientsupernova dust production, the observed dust mass of the z = 7.5 galaxy A1689-zD1 requires very efficient grain growth. This, in turn, implies that in this galaxy the average density of the cold and dense gas, where grain growth occurs, is comparable to that inferred from observations of QSO host galaxies at similar redshifts. Although plausible, the upper limits on the dust continuum emission of galaxies at 6.5 < z < 7.5 show that these conditions must not apply to the bulk of the high redshift galaxy population.

Measuring line-of-sight-dependent Fourier-space clustering using FFTs

Abstract: Observed galaxy clustering exhibits local transverse statistical isotropy around the lineof-sight (LOS). The variation of the LOS across a galaxy survey complicates the measurement of the observed clustering as a function of the angle to the LOS, as fast Fourier transforms (FFTs) based on Cartesian grids, cannot individually allow for this. Recent advances in methodology for calculating LOS-dependent clustering in Fourier space include the realization that power spectrum LOS-dependent moments can be constructed from sums over galaxies, based on approximating the LOS to each pair of galaxies by the LOS to one of them. We show that we can implement this method using multiple FFTs, each measuring the LOS-weighted clustering along different axes. The N logN nature of FFTs means that the computational speed-up is a factor of > 1000 compared with summing over galaxies. This development should be beneficial for future projects such as DESI and Euclid which will provide an order of magnitude more galaxies than current surveys.

X-ray coherent pulsations during a sub-luminous accretion disc state of the transitional millisecond pulsar XSS J12270−4859

Abstract: We present the first detection of X-ray coherent pulsations from the transitional millisecond pulsar XSS J12270-4859, while it was in a sub-luminous accretion disk state characterized by a 0.5-10 keV luminosity of 5E33 erg/s (assuming a distance of 1.4 kpc). Pulsations were observed by XMM-Newton at an rms amplitude of (7.7 +/- 0.5)% with a second harmonic stronger than the the fundamental frequency, and were detected when the source is neither flaring nor dipping. The most likely interpretation of this detection is that matter from the accretion disk was channelled by the neutron star magnetosphere and accreted onto its polar caps. According to standard disk accretion theory, for pulsations to be observed the mass in-flow rate in the disk was likely larger than the amount of plasma actually reaching the neutron star surface; an outflow launched by the fast rotating magnetosphere then probably took place, in agreement with the observed broad-band spectral energy distribution. We also report about the non-detection of X-ray pulsations during a recent observation performed while the source behaved as a rotationally-powered radio pulsar.

Separate universe simulations

Abstract: The large-scale statistics of observables such as the galaxy density are chiey determined by their dependence on the local coarse-grained matter density. This dependence can be measured directly and eciently in N-body simulations by using the fact that a uniform density perturbation with respect to some ducial background cosmology is equivalent to modifying the background and including curvature, i.e., by simulating a \separate universe". We derive this mapping to fully non-linear order, and provide a step-by-step description of how to perform and analyse the separate universe simulations. This technique can be applied to a wide range of observables. As an example, we calculate the response of the non-linear matter power spectrum to long-wavelength density perturbations, which corresponds to the angle-averaged squeezed limit of the matter bispectrum and higher n-point functions. Using only a modest simulation volume, we obtain results with percent-level precision over a wide range of scales.

Binary black hole accretion during inspiral and merger

Abstract: We present the results of 2D, moving mesh, viscous hydrodynamical simulations of accretion onto merging supermassive black hole (SMBH) binaries. We include viscous heating, shock heating, and radiative cooling, and simulate the transition from the "pre-decoupling" epoch, where the inspiral timescale is longer than the viscous timescale, to the "post-decoupling" epoch, where the inspiral timescale is shorter than the viscous timescale. We find that there is no abrupt halt to the accretion at decoupling, but rather the accretion shows a slow decay, with significant accretion well after the expected decoupling. Moreover, we find that the luminosity in X-rays is significantly higher prior to the merger, as orbital energy from the SMBH binary is converted to heat via strong shocks inside the cavity, and radiated away. Following the merger, the cavity refills viscously and the accretion rate relaxes to the Shakura-Sunyaev value, while the X-ray luminosity drops as the shocks quickly dissipate.

Recombination energy in double white dwarf formation

Abstract: In this Letter we investigate the role of recombination energy during a common envelope event. We confirm that taking this energy into account helps to avoid the formation of the circumbinary envelope commonly found in previous studies. For the first time, we can model a complete common envelope event, with a clean compact double white dwarf binary system formed at the end. The resulting binary orbit is almost perfectly circular. In addition to considering recombination energy, we also show that between 1/4 and 1/2 of the released orbital energy is taken away by the ejected material. We apply this new method to the case of the double-white dwarf system WD 1101+364, and we find that the progenitor system at the start of the common envelope event consisted of a $\sim1.5M_\odot$ red giant star in a $\sim 30$ day orbit with a white dwarf companion.

Characteristic signatures in the thermal emission from accreting binary black holes

Abstract: We present the results of a calculation of the thermal spectrum from a 2D, moving mesh, high-accuracy, viscous hydrodynamical simulation of an accreting supermassive black hole binary. We include viscous heating, shock heating, and radiative cooling, evolving for longer than a viscous time so that we reach a quasi-steady accretion state. In agreement with previous work, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We also find that emission from the shock-heated mini-disks and accretion streams prevents any deficit in high-energy emission that may be expected inside the circumbinary cavity, and instead leads to a characteristic brightening of the spectrum beginning in soft X-rays.

Type Ia supernovae with bimodal explosions are common – possible smoking gun for direct collisions of white dwarfs

Abstract: We discover clear doubly-peaked line profiles in 3 out of ~20 type Ia supernovae (SNe Ia) with high-quality nebular-phase spectra. The profiles are consistently present in three well-separated Co/Fe emission features. The two peaks are respectively blue-shifted and red-shifted relative to the host galaxies and are separated by ~5000 km/s. The doubly-peaked profiles directly reflect a bi-modal velocity distribution of the radioactive Ni56 in the ejecta that powers the emission of these SNe. Due to their random orientations, only a fraction of SNe with intrinsically bi-modal velocity distributions will appear as doubly-peaked spectra. Therefore SNe with intrinsic bi-modality are likely common, especially among the SNe in the low-luminosity part on the Philips relation (\Delta m15(B) >~ 1.3; ~40% of all SNe Ia). Such bi-modality is naturally expected from direct collisions of white dwarfs (WDs) due to the detonation of both WDs and is demonstrated in a 3D 0.64 M_Sun-0.64 M_Sun WD collision simulation. In the future, with a large sample of nebular spectra and a comprehensive set of numerical simulations, the collision model can be unambiguously tested as the primary channel for type Ia SNe, and the distribution of nebular line profiles will either be a smoking gun or rule it out.

Dark matter-induced collapse of neutron stars: a possible link between fast radio bursts and the missing pulsar problem

Abstract: Fast radio bursts (FRBs) are an emerging class of short and bright radio transients whose sources remain enigmatic. Within the Galactic Centre, the non-detection of pulsars within the inner ∼10 pc has created a missing pulsar problem that has intensified with time. With all reserve, we advance the notion that the two problems could be linked by a common solution: the collapse of neutron stars (NS) due to capture and sedimentation of dark matter (DM) within their cores. Bramante & Linden showed that certain DM properties allow for rapid NS collapse within the high DM density environments near galactic centres while permitting NS survival elsewhere. Each DM-induced collapse could generate an FRB as the NS magnetosphere is suddenly expelled. This scenario could explain several features of FRBs: their short time scales, large energies, locally produced scattering tails, and high event rates. We predict that FRBs are localized to galactic centres, and that our own galactic centre harbours a large population of NS-mass (M ∼ 1.4 M⊙) black holes. The DM-induced collapse scenario is intrinsically unlikely because it can only occur in a small region of allowable DM parameter space. However, if observed to occur, it would place tight constraints on DM properties.

Einstein's legacy in galaxy surveys

Abstract: Non-Gaussianity in the primordial fluctuations that seeded structure formation produces a signal in the galaxy power spectrum on very large scales. This signal contains vital information about the primordial Universe, but it is very challenging to extract, because of cosmic variance and large-scale systematics-especially after the Planck experiment has already ruled out a large amplitude for the signal. Whilst cosmic variance and experimental systematics can be alleviated by the multitracer method, we here address another systematic-introduced by not using the correct relativistic analysis of the power spectrum on very large scales. In order to reduce the errors on f NL , we need to include measurements on the largest possible scales. Failure to include the relativistic effects on these scales can introduce significant bias in the best-fit value of f NL from future galaxy surveys.

Peanuts, brezels and bananas: food for thought on the orbital structure of the Galactic bulge

Abstract: Recent observations have discovered the presence of a Box/Peanut or X-shape structure in the Galactic bulge. Such Box/Peanut structures are common in external disc galaxies, and are well-known in N-body simulations where they form following the buckling instability of a bar. From studies of analytical potentials and N-body models it has been claimed in the past that Box/Peanut bulges are supported by "bananas", or x1v1 orbits. We present here a set of N-body models where instead the peanut bulge is mainly supported by brezel-like orbits, allowing strong peanuts to form with short extent relative to the bar length. This shows that stars in the X-shape do not necessarily stream along banana orbits which follow the arms of the X-shape. The brezel orbits are also found to be the main orbital component supporting the peanut shape in our recent Made-to-Measure dynamical models of the Galactic bulge. We also show that in these models the fraction of stellar orbits that contribute to the X-structure account for 40-45% of the stellar mass.

AGN flickering and chaotic accretion

Abstract: Observational arguments suggest that the growth phases of the supermassive black holes in active galactic nuclei have a characteristic timescale $\sim 10^5$ yr. We show that this is the timescale expected in the chaotic accretion picture of black hole feeding, because of the effect of self-gravity in limiting the mass of any accretion disc feeding event.

Black hole jets without large-scale net magnetic flux

Abstract: We propose a scenario for launching relativistic jets from rotating black holes, in which small-scale magnetic flux loops, sustained by disc turbulence, are forced to inflate and open by differential rotation between the black hole and the accretion flow. This mechanism does not require a large-scale net magnetic flux in the accreting plasma. Estimates suggest that the process could operate effectively in many systems, and particularly naturally and efficiently when the accretion flow is retrograde. We present the results of general-relativistic force-free electrodynamic simulations demonstrating the time evolution of the black hole's magnetosphere, the cyclic formation of jets, and the effect of magnetic reconnection. The jets are highly variable on timescales ~ 10-10^3 r_ g/c, where r_g is the black hole's gravitational radius. The reconnecting current sheets observed in the simulations may be responsible for the hard X-ray emission from accreting black holes.

Detection of two power-law tails in the probability distribution functions of massive GMCs

Abstract: We report the novel detection of complex high-column density tails in the probability distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, MonR2, NGC6334), obtained from dust emission observed with Herschel. The low column density range can be fit with a lognormal distribution. A first power-law tail starts above an extinction (Av) of ~6-14. It has a slope of alpha=1.3-2 for the rho~r^-alpha profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above Av~40, 60, and 140, we detect an excess that can be fitted by a flatter power law tail with alpha>2. It correlates with the central regions of the cloud (ridges/hubs) of size ~1 pc and densities above 10^4 cm^-3. This excess may be caused by physical processes that slow down collapse and reduce the flow of mass towards higher densities. Possible are: 1. rotation, which introduces an angular momentum barrier, 2. increasing optical depth and weaker cooling, 3. magnetic fields, 4. geometrical effects, and 5. protostellar feedback. The excess/second power-law tail is closely linked to high-mass star-formation though it does not imply a universal column density threshold for the formation of (high-mass) stars.

Photon underproduction crisis: Are QSOs sufficient to resolve it?

Abstract: We investigate the recent claim of ‘photon underproduction crisis’ by Kollmeier et al. (2014) which suggests that the known sources of ultra-violet (UV) radiation may not be sufficient to generate the inferred Hi photoionization rate ( HI) in the low redshift inter-galactic medium. Using the updated QSO emissivities from the recent studies and our radiative transfer code developed to estimate the UV background, we show that the QSO contributions to HI is higher by a factor �2 as compared to the previous estimates. Using self-consistently computed combinations of star formation rate density and dust attenuation, we show that a typical UV escape fraction of 4% from star forming galaxies should be sufficient to explain the inferred HI by Kollmeier et al.(2014). Interestingly, we find that the contribution from QSOs alone can explain the recently inferred HI by Shull et al.(2015) which used the same observational data but different simulation. Therefore, we conclude that the crisis is not as severe as it was perceived before and there seems no need to look for alternate explanations such as low luminosity hidden QSOs or decaying dark matter particles.

ALMA imaging of SDP.81 – I. A pixelated reconstruction of the far-infrared continuum emission

Abstract: We present a sub-50 parsec scale analysis of the gravitational lens system SDP.81 at redshift 3.042 using Atacama Large Millimetre/submillimetre Array science verification data. We model both the mass distribution of the gravitational lensing galaxy and the pixelated surface brightness distribution of the background source using a novel Bayesian technique that fits the data directly in visibility space. We find the 1 and 1.3 mm dust emission to be magnified by a factor of μtot = 17.6 ± 0.4, giving an intrinsic total star formation rate of 315 ± 60 M⊙ yr-1 and a dust mass of 6.4 ± 1.5 × 108 M⊙. The reconstructed dust emission is found to be non-uniform, but composed of multiple regions that are heated by both diffuse and strongly clumped star formation. The highest surface brightness region is a ˜1.9 × 0.7 kpc disc-like structure, whose small extent is consistent with a potential size-bias in gravitationally lensed starbursts. Although surrounded by extended star formation, with a density of 20-30 ± 10 M⊙ yr-1 kpc-2, the disc contains three compact regions with densities that peak between 120 and 190 ± 20 M⊙ yr-1 kpc-2. Such star formation rate densities are below what is expected for Eddington-limited star formation by a radiation pressure supported starburst. There is also a tentative variation in the spectral slope of the different star-forming regions, which is likely due to a change in the dust temperature and/or opacity across the source.

Refined critical balance in strong Alfvénic turbulence

Abstract: We present numerical evidence that in strong Alfvenic turbulence, the critical balance principle---equality of the nonlinear decorrelation and linear propagation times---is scale invariant, in the sense that the probability distribution of the ratio of these times is independent of scale. This result only holds if the local alignment of the Elsasser fields is taken into account in calculating the nonlinear time. At any given scale, the degree of alignment is found to increase with fluctuation amplitude, supporting the idea that the cause of alignment is mutual dynamical shearing of Elsasser fields. The scale-invariance of critical balance (while all other quantities of interest are strongly intermittent, i.e., have scale-dependent distributions) suggests that it is the most robust of the scaling principles used to describe Alfvenic turbulence. The quality afforded by situ fluctuation measurements in the solar wind allows for direct verification of this fundamental principle.

Fast cold gas in hot AGN outflows

Abstract: Observations of the emission from spatially extended cold gas around bright high-redshift quasars (QSOs) reveal surprisingly large velocity widths exceeding 2000kms −1 , out to projected distances as large as 30kpc. The high-velocity widths have been interpreted as the signature of powerful AGN-driven outflows. Naively, these findings appear in tension with hydrodynamicmodelsinwhichAGN-drivenoutflowsareenergy-drivenandthusveryhotwith typical temperatures T 10 6−7 K. Using the moving-mesh code AREPO, we perform ‘zoom-in’ cosmologicalsimulationsofa z ∼6QSOanditsenvironment,followingblackholegrowthand feedback via energy-driven outflows. In the simulations, the QSO host galaxy is surrounded by a clumpy circumgalactic medium pre-enriched with metals due to supernovae-driven galactic outflows. As a result, part of the AGN-driven hot outflowing gas can cool radiatively, leading to large amounts (10 9 M� ) of cold gas comoving with the hot bipolar outflow. This results in velocity widths of spatially extended cold gas similar to those observed. We caution, however, thatgasinflows,randommotionsinthedeeppotentialwelloftheQSOhostgalaxyandcooling of supernovae-driven winds contribute significantly to the large velocity width of the cold gas in the simulations, complicating the interpretation of observational data.

The reversal of the SF–density relation in a massive, X-ray-selected galaxy cluster at z = 1.58: results from Herschel

Abstract: Dusty, star-forming galaxies have a critical role in the formation and evolution of massive galaxies in the Universe. Using deep far-infrared imaging in the range 100-500um obtained with the Herschel telescope, we investigate the dust-obscured star formation in the galaxy cluster XDCP J0044.0-2033 at z=1.58, the most massive cluster at z >1.5, with a measured mass M200= 4.7x10$^{14}$ Msun. We perform an analysis of the spectral energy distributions (SEDs) of 12 cluster members (5 spectroscopically confirmed) detected with >3$\sigma$ significance in the PACS maps, all ULIRGs. The individual star formation rates (SFRs) lie in the range 155-824 Ms/yr, with dust temperatures of 24$\pm$35 K. We measure a strikingly high amount of star formation (SF) in the cluster core, SFR ( 1875$\pm$158 Ms/yr, 4x higher than the amount of star formation in the cluster outskirts. This scenario is unprecedented in a galaxy cluster, showing for the first time a reversal of the SF-density relation at z~1.6 in a massive cluster.

An ultraluminous nascent millisecond pulsar

Abstract: If the ultraluminous source (ULX) M82 X-2 sustains its measured spin-up value of $\dot u= 10^{-10}\,{\rm s^{-2}}$, it will become a millisecond pulsar in less than $10^5\,$ years. The observed (isotropic) luminosity of $10^{40}\,$ erg/s also supports the notion that the neutron star will spin up to a millisecond period upon accreting about $0.1\,{\rm M_\odot}$---the reported hard X-ray luminosity of this ULX, together with the spin-up value, implies torques consistent with the accretion disk extending down to the vicinity of the stellar surface, as expected for low values of the stellar dipole magnetic field ($B\lesssim 10^9\,$G). This suggests a new channel of millisecond pulsar formation---in high-mass X-ray binaries (HMXBs)---and may have implications for studies of gravitational waves, and possibly for the formation of low-mass black holes through accretion-induced collapse.

Black hole accretion versus star formation rate: theory confronts observations

Abstract: We use a suite of hydrodynamical simulations of galaxy mergers to compare star formation rate (SFR) and black hole accretion rate (BHAR) for galaxies before the interaction ('stochastic' phase), during the `merger' proper, lasting ~0.2-0.3 Gyr, and in the `remnant' phase. We calculate the bi-variate distribution of SFR and BHAR and define the regions in the SFR-BHAR plane that the three phases occupy. No strong correlation between BHAR and galaxy-wide SFR is found. A possible exception are galaxies with the highest SFR and the highest BHAR. We also bin the data in the same way used in several observational studies, by either measuring the mean SFR for AGN in different luminosity bins, or the mean BHAR for galaxies in bins of SFR. We find that the apparent contradiction or SFR versus BHAR for observed samples of AGN and star forming galaxies is actually caused by binning effects. The two types of samples use different projections of the full bi-variate distribution, and the full information would lead to unambiguous interpretation. We also find that a galaxy can be classified as AGN-dominated up to 1.5 Gyr after the merger-driven starburst took place. Our study is consistent with the suggestion that most low-luminosity AGN hosts do not show morphological disturbances.

Galaxy pairs in the Sloan Digital Sky Survey – XII. The fuelling mechanism of low-excitation radio-loud AGN

Abstract: We investigate whether the fuelling of low excitation radio galaxies (LERGs) is linked to major galaxy interactions. Our study utilizes a sample of 10,800 spectroscopic galaxy pairs and 97 post-mergers selected from the Sloan Digital Sky Survey with matches to multi-wavelength datasets. The LERG fraction amongst interacting galaxies is a factor of 3.5 higher than that of a control sample matched in local galaxy density, redshift and stellar mass. However, the LERG excess in pairs does not depend on projected separation and remains elevated out to at least 500 h 1 70 kpc, suggesting that major mergers are not their main fuelling channel. In order to identify the primary fuelling mechanism of LERGs, we compile samples of control galaxies that are matched in various host galaxy and environmental properties. The LERG excess is reduced, but not completely removed, when halo mass or D4000 are included in the matching parameters. However, when both Mhalo and D4000 are matched, there is no LERG excess and the 1.4 GHz luminosities (which trace jet mechanical power) are consistent between the pairs and control. In contrast, the excess of optical and mid-IR selected AGN in galaxy pairs is unchanged when the additional matching parameters are implemented. Our results suggest that whilst major interactions may trigger optically and mid-IR selected AGN, the gas which fuels the LERGs has two secular origins: one associated with the large scale environment, such as accretion from the surrounding medium or minor mergers, plus an internal stellar mechanism, such as winds from evolved stars.

PSR J1738+0333: the first millisecond pulsar + pulsating white dwarf binary.

Abstract: We report the discovery of the first millisecond pulsar with a pulsating white dwarf companion. Following the recent discoveries of pulsations in extremely low-mass (ELM, <0.3 Msol) white dwarfs (WDs), we targeted ELM WD companions to two millisecond pulsars with high-speed Gemini photometry. We find significant optical variability in PSR J1738+0333 with periods between roughly 1790-3060 s, consistent in timescale with theoretical and empirical observations of pulsations in 0.17 Msol He-core ELM WDs. We additionally put stringent limits on a lack of variability in PSR J1909-3744, showing this ELM WD is not variable to <0.1 per cent amplitude. Thanks to the accurate distance and radius estimates from radio timing measurements, PSR J1738+0333 becomes a benchmark for low-mass, pulsating WDs. Future, more extensive time-series photometry of this system offers an unprecedented opportunity to constrain the physical parameters (including the cooling age) and interior structure of this ELM WD, and in turn, the mass and spin-down age of its pulsar companion.

The colours of satellite galaxies in the Illustris simulation

Abstract: Observationally, the fraction of blue satellite galaxies decreases steeply with host halo mass, and their radial distribution around central galaxies is significantly shallower in massive (M* ≥ 1011 M⊙) than in Milky Way-like systems. Theoretical models, based primarily on semi-analytical techniques, have had a long-standing problem with reproducing these trends, instead predicting too few blue satellites in general but also estimating a radial distribution that is too shallow, regardless of primary mass. In this Letter, we use the Illustris cosmological simulation to study the properties of satellite galaxies around isolated primaries. For the first time, we find good agreement between theory and observations. We identify the main source of this success relative to earlier work to be a consequence of the large gas contents of satellites at infall, a factor ∼5–10 times larger than in semi-analytical models. Because of their relatively large gas reservoirs, satellites can continue to form stars long after infall, with a typical time-scale for star-formation to be quenched ∼2 Gyr in groups but more than ∼5 Gyr for satellites around Milky Way-like primaries. The gas contents we infer are consistent with z = 0 observations of H I gas in galaxies, although we find large discrepancies among reported values in the literature. A testable prediction of our model is that the gas-to-stellar mass ratio of satellite progenitors should vary only weakly with cosmic time.