Detached Shell Carbon Stars: Tracing Thermal Pulses on the Asymptotic Giant Branch

We consider whether the subset of carbon-rich asymptotic giant branch (AGB) stars that exhibit detached, expanding circumstellar shells may reveal the past histories of these stars as having undergone helium shell flashes (thermal pulses) on the AGB. We exploit newly available Gaia parallaxes and photometry, along with archival infrared photometry, to obtain refined estimates of the luminosities of all (12) known detached shell carbon stars. We examine the relationship between these luminosities and the estimated dynamical ages (ejection times) of the detached shells associated with the 12 stars, which range from ∼1000 to ∼30,000 yr. When arranged according to detached shell dynamical age, the (implied) luminosity evolution of the known detached shell carbon stars closely follows the predicted “light curves” of individual thermal pulses obtained from models of AGB stars. The comparison between data and models suggests that detached shell carbon stars are descended from ∼2.5 to 4.0 M ⊙ progenitors. We conclude that detached shell carbon stars may serve as effective tracers of the luminosity evolution of AGB thermal pulses.

The evolution of the low-frequency radio AGN population to z ≃ 1.5 in the ELAIS N1 field

ABSTRACT We study the cosmic evolution of radio sources out to z ≃ 1.5 using a GMRT 610 MHz survey covering ∼1.86 deg2 of the ELAIS N1 field with a minimum/median rms noise 7.1/19.5 μJy beam−1 and an angular resolution of 6 arcsec. We classify sources as star forming galaxies (SFGs), radio-quiet (RQ) and radio-loud (RL) Active Galactic Nuclei (AGNs) using a combination of multiwavelength diagnostics and find evidence in support of the radio emission in SFGs and RQ AGN arising from star formation, rather than AGN-related processes. At high luminosities, however, both SFGs and RQ AGN display a radio excess when comparing radio and infrared star formation rates. The vast majority of our sample lie along the $\rm {SFR - \mathit{ M}_{\star }}$ ‘main sequence’ at all redshifts when using infrared star formation rates. We derive the 610 MHz radio luminosity function for the total AGN population, constraining its evolution via continuous models of pure density and pure luminosity evolution with $\rm {\Phi ^{\star }\, \propto \, (\, 1+\, \mathit{ z})^{(2.25\pm 0.38)-(0.63\pm 0.35)z}}$ and $\rm {\mathit{ L}_{610\, MHz}\, \propto \, (\, 1+\, \mathit{ z})^{(3.45\pm 0.53)-(0.55\pm 0.29)\mathit{ z}}}$, respectively. For our RQ and RL AGN, we find a fairly mild evolution with redshift best fitted by pure luminosity evolution with $\rm {\mathit{ L}_{610\, MHz}\, \propto \, (\, 1+\, \mathit{ z})^{(2.81\pm 0.43)-(0.57\pm 0.30)\mathit{ z}}}$ for RQ AGN and $\rm {\mathit{ L}_{610\, MHz}\, \propto \, (\, 1+\, \mathit{ z})^{(3.58\pm 0.54)-(0.56\pm 0.29)\mathit{ z}}}$ for RL AGN. The 610 MHz radio AGN population thus comprises two differently evolving populations whose radio emission is mostly SF-driven or AGN-driven, respectively.

Predicted rates of merging neutron stars in galaxies

ABSTRACT We compute rates of merging neutron stars (MNS) in different galaxies, as well as the cosmic MNS rate in different cosmological scenarios. Our aim is to provide predictions of kilonova rates for future observations both at low and high redshift. In the adopted galaxy models, the production of r-process elements either by MNS or core-collapse supernovae is taken into account. To compute the MNS rates, we adopt either a constant total time delay for merging (10 Myr) or a distribution function of such delays. We conclude (i) the observed present time MNS rate in our Galaxy is well reproduced either with a constant time delay or a distribution function ∝t−1. (ii) The [Eu/Fe] versus [Fe/H] relation can be well reproduced with only MNS, if the time delay is short and constant. If a distribution function of delays is adopted, core-collapse supernovae are also required. (iii) The present time cosmic MNS rate can be well reproduced in several cosmological scenarios. (iv) Spiral galaxies are the major contributors to the cosmic MNS at all redshifts in hierarchical scenarios. In the pure luminosity evolution scenario, the spirals are the major contributors locally, whereas at high redshift ellipticals dominate. (v) The predicted cosmic MNS rate well agrees with the cosmic rate of short gamma-ray bursts, if the distribution function of delays is adopted in a cosmological hierarchical scenario observationally derived. (vi) Future observations of kilonovae in ellipticals will allow us to disentangle among constant or a distribution of time delays and among different cosmological scenarios.

The XXL Survey

We studied a sample of 274 radio and X-ray selected quasars (XQSOs) detected in the COSMOS and XXL-S radio surveys at 3 GHz and 2.1 GHz, respectively. This sample was identified by adopting a conservative threshold in X-ray luminosity, LX [2−10 keV] ≥ 1044 erg s−1, selecting only the most powerful quasars. A number of previous studies on the origin of radio emission in type-1 quasars have focused on the radio loudness distributions, some claiming to have found evidence for bimodality, pointing toward the existence of two physically different mechanisms for the radio emission. Using available multiwavelength data, we examined various criteria for the selection of radio-loud (RL) and radio-quiet (RQ) XQSOs and found that the number of RL/RQ XQSOs changes significantly depending on the chosen criterion. This discrepancy arises due to the different criteria tracing different physical processes and due to the fact that our sample was selected from flux-limited radio and X-ray surveys. Another approach to study the origin of radio emission in XQSOs is via their radio luminosity functions (RLF). We constructed the XQSO 1.4 GHz RLFs in six redshift bins at 0.5 ≤ z ≤ 3.75. The lower-1.4 GHz luminosity end shows a higher normalization than expected only from AGN contribution in all studied redshift bins. We found that the so-called “bump” is mostly dominated by emission due to star-forming processes within the host galaxies of XQSOs. As expected, AGN-related radio emission is the dominant contribution at the higher-luminosity end of RLF. To study the evolution of the XQSO RLF, we used a combination of analytic forms from the literature to constrain the “bump” due to star formation and the higher-luminosity AGN part of the RLF. We defined two 1.4 GHz luminosity thresholds, Lth, SF and Lth, AGN, below and above which more than 80% of sources contributing to the RLF are dominated by star formation and AGN-related activity, respectively. The two thresholds evolve with redshift, which is most likely driven by the strong evolution of star formation rates of the XQSO host galaxies. We found that both the lower and higher luminosity ends evolve significantly in density, while their luminosity evolution parameters are consistent with being constant. We found that the lower-luminosity end evolves both in density and luminosity, while the higher-luminosity end evolves significantly only in density. Our results expose the dichotomy of the origin of radio emission: while the higher-luminosity end of the XQSO RLF is dominated by AGN activity, the lower-luminosity end is dominated by the star formation-related processes.

The luminosity evolution of nova shells

Over the last decade, nova shells have been discovered around a small number of cataclysmic variables that had not been known to be post-novae, while other searches around much larger samples have been mostly unsuccessful. This raises the question about how long such shells are detectable after the eruption and whether this time limit depends on the characteristics of the nova. So far, there has been only one comprehensive study of the luminosity evolution of nova shells, undertaken almost two decades ago. Here, we present a re-analysis of the Hα and [O III] flux data from that study, determining the luminosities while also taking into account newly available distances and extinction values, and including additional luminosity data of “ancient” nova shells. We compare the long-term behaviour with respect to nova speed class and light curve type. We find that, in general, the luminosity as a function of time can be described as consisting of three phases: an initial shallow logarithmic decline or constant behaviour, followed by a logarithmic main decline phase, with a possible return to a shallow decline or constancy at very late stages. The luminosity evolution in the first two phases is likely to be dominated by the expansion of the shell and the corresponding changes in volume and density, while for the older nova shells, the interaction with the interstellar medium comes into play. The slope of the main decline is very similar for almost all groups for a given emission line, but it is significantly steeper for [O III], compared to Hα, which we attribute to the more efficient cooling provided by the forbidden lines. The recurrent novae are among the notable exceptions, along with the plateau light curve type novae and the nova V838 Her. We speculate that this is due to the presence of denser material, possibly in the form of remnants from previous nova eruptions, or of planetary nebulae, which might also explain some of the brighter ancient nova shells. While there is no significant difference in the formal quality of the fits to the decline when grouped according to light curve type or to speed class, the former presents less systematic scatter. It is also found to be advantageous in identifying points that would otherwise distort the general behaviour. As a by-product of our study, we revised the identification of all novae included in our investigation with sources in the Gaia Data Release 2 catalogue.

SN 2018hti: a nearby superluminous supernova discovered in a metal-poor galaxy

ABSTRACT SN 2018hti is a Type I superluminous supernova (SLSN I) with an absolute g-band magnitude of −22.2 at maximum brightness, discovered by the Asteroid Terrestrial-impact Last Alert System in a metal-poor galaxy at a redshift of 0.0612. We present extensive photometric and spectroscopic observations of this supernova, covering the phases from ∼−35 d to more than +340 d from the r-band maximum. Combining our BVgri-band photometry with Swift UVOT optical/ultraviolet photometry, we calculated the peak luminosity as ∼3.5 × 1044 erg s−1. Modelling the observed light curve reveals that the luminosity evolution of SN 2018hti can be produced by an ejecta mass of 5.8 M⊙ and a magnetar with a magnetic field of B = 1.8 × 1013 G having an initial spin period of P0 = 1.8 ms. Based on such a magnetar-powered scenario and a larger sample, a correlation between the spin of the magnetar and the kinetic energy of the ejecta can be inferred for most SLSNe I, suggesting a self-consistent scenario. Like for other SLSNe I, the host galaxy of SN 2018hti is found to be relatively faint (Mg = −17.75 mag) and of low metallicity (Z = 0.3 Z⊙), with a star formation rate of 0.3 M⊙ yr−1. According to simulation results of single-star evolution, SN 2018hti could originate from a massive, metal-poor star with a zero-age main sequence (ZAMS) mass of 25–40 M⊙, or from a less massive rotating star with MZAMS ≈ 16–25 M⊙. For the case of a binary system, its progenitor could also be a star with $M_\mathrm{ZAMS} \gtrsim 25\, \mathrm{ M}_\odot$.

Infrared luminosity functions and dust mass functions in the EAGLE simulation

ABSTRACT We present infrared luminosity functions and dust mass functions for the EAGLE cosmological simulation, based on synthetic multiwavelength observations generated with the SKIRT radiative transfer code. In the local Universe, we reproduce the observed infrared luminosity and dust mass functions very well. Some minor discrepancies are encountered, mainly in the high luminosity regime, where the EAGLE-SKIRT luminosity functions mildly but systematically underestimate the observed ones. The agreement between the EAGLE-SKIRT infrared luminosity functions and the observed ones gradually worsens with increasing lookback time. Fitting modified Schechter functions to the EAGLE-SKIRT luminosity and dust mass functions at different redshifts up to z = 1, we find that the evolution is compatible with pure luminosity/mass evolution. The evolution is relatively mild: within this redshift range, we find an evolution of L⋆,250 ∝ (1 + z)1.68, L⋆,TIR ∝ (1 + z)2.51 and M⋆,dust ∝ (1 + z)0.83 for the characteristic luminosity/mass. For the luminosity/mass density we find ε250 ∝ (1 + z)1.62, εTIR ∝ (1 + z)2.35, and ρdust ∝ (1 + z)0.80, respectively. The mild evolution of the dust mass density is in relatively good agreement with observations, but the slow evolution of the infrared luminosity underestimates the observed luminosity evolution significantly. We argue that these differences can be attributed to increasing limitations in the radiative transfer treatment due to increasingly poorer resolution, combined with a slower than observed evolution of the SFR density in the EAGLE simulation and the lack of AGN emission in our EAGLE-SKIRT post-processing recipe.