Search

[Abridged] The relation between nuclear star clusters (NSCs) and the growth of the central SMBHs, as well as their connection to the properties of the host galaxies, is crucial for understanding the evolution of galaxies. Recent observations have revealed that about 10 per cent of nucleated galaxies host hybrid nuclei, consisting of both NSCs and accreting SMBHs that power active galactic nuclei (AGN). Motivated by the potential of the recently published multi-wavelength data sets from LeMMINGs survey, here we present the most thorough investigation to date of the incidence of hybrid nuclei in a large sample of 100 nearby nucleated galaxies (10 E, 25 S0, 63 S, and 2 Irr), covering a wide range in stellar mass ($M_{*,\rm gal} \sim 10^{8.7}-10^{12}~\rm M_{sun}$). We identify the nuclei and derive their properties by performing detailed 1D and 2D multi-component decompositions of the optical and near-infrared $HST$ stellar light distributions of the galaxies using S\'ersic and core-S\'ersic models. Our AGN diagnostics are based on homogeneously derived nuclear 1.5 GHz $e$-MERLIN radio, $Chandra$ X-ray (0.3--10 keV) and optical emission-line data. We determine the nucleation fraction ($f_{\rm nuc} $) as the relative incidence of nuclei across the LeMMINGs $HST$ sample and find $f_{\rm nuc} =~ $100/149 (= 67 $\pm$ 7 per cent), confirming previous work, with a peak value of 49/56~(= $88 \pm 13$ per cent) at bulge masses $M_{*,\rm bulge} \sim 10^{9.4}$- $10^{10.8}~\rm M_{sun}$. We identify 30 nucleated LeMMINGs galaxies that are optically active, radio-detected and X-ray luminous ($L_{X} > 10^{39}$ erg s$^{-1}$). This indicates that our nucleated sample has a lower limit $\sim$ 30 per cent occupancy of hybrid nuclei, which is a function of $M_{*,\rm bulge}$ and $M_{*,\rm gal}$. We find that hybrid nuclei have a number density of $(1.5 \pm 0.4) \times 10^{-5}$ Mpc$^{-3}$., Comment: Accepted for publication in MNRAS (21 pages, 10 figures, 9 tables)

Multiwavelength studies indicate that nuclear activity and bulge properties are closely related, but the details remain unclear. To study this further, we combine $Hubble~Space~Telescope$ bulge structural and photometric properties with 1.5 GHz, $e$-MERLIN nuclear radio continuum data from the LeMMINGs survey for a large sample of 173 `active' galaxies (LINERs and Seyferts) and `inactive' galaxies (H IIs and absorption line galaxies, ALGs). Dividing our sample into active and inactive, they define distinct (radio core luminosity)$-$(bulge mass), L_R,core-M_*,bulge, relations, with a mass turnover at M_*, bulge ~ 10^(9.8 +- 0.3) M_sun (supermassive black hole mass M_BH ~ 10^(6.8 +- 0.3) M_sun), which marks the transition from AGN-dominated nuclear radio emission in more massive bulges to that mainly driven by stellar processes in low-mass bulges. None of our 10/173 bulgeless galaxies host an AGN. The AGN fraction increases with increasing M_*, bulge such that f_optical_AGN $\propto$ M_*,bulge^(0.24 +- 0.06) and f_radio_AGN $\propto$ M_*,bulge^(0.24 +- 0.05). Between M_*,bulge ~ 10^8.5 and 10^11.3 M_sun, f_optical_AGN steadily rises from 15 +- 4 to 80 +- 5 per cent. We find that at fixed bulge mass, the radio loudness, nuclear radio activity and the (optical and radio) AGN fraction exhibit no dependence on environment. Radio-loud hosts preferentially possess an early-type morphology than radio-quiet hosts, the two types are however indistinguishable in terms of bulge S\'ersic index and ellipticity, while results on the bulge inner logarithmic profile slope are inconclusive. We finally discuss the importance of bulge mass in determining the AGN triggering processes, including potential implications for the nuclear radio emission in nearby galaxies., Comment: 27 pages, 15 figures, 4 tables, accepted for publication in MNRAS

[Abridged] We use high-resolution $HST$ imaging and $e$-MERLIN 1.5-GHz observations of galaxy cores from the LeMMINGs survey to investigate the relation between optical structural properties and nuclear radio emission for a large sample of galaxies. We perform accurate, multi-component decompositions of new surface brightness profiles extracted from $HST$ images for 163 LeMMINGs galaxies and fit up to six galaxy components (e.g., bulges, discs, AGN, bars, rings, spiral arms, and nuclear star clusters) simultaneously with S\'ersic and/or core-S\'ersic models. By adding such decomposition data for 10 LeMMINGs galaxies from our past work, the final sample of 173 nearby galaxies (102 Ss, 42 S0s, 23 Es plus 6 Irr) with bulge stellar mass (typically) M_*, bulge ~ 10^6-10^12.5 M_sun, encompasses all optical spectral classes (LINER, Seyfert, ALG and H II). We show that the bulge mass can be significantly overestimated in many galaxies when components such as bars, rings and spirals are not included in the fits. We additionally implement a Monte Carlo method to determine errors on bulge, disc and other fitted structural parameters. Moving (in the opposite direction) across the Hubble sequence, i.e., from the irregular to elliptical galaxies, we confirm that bulges become larger, more prominent and round. Such bulge dominance is associated with a brighter radio core luminosity. We also find that the radio detection fraction increases with bulge mass. At M_*,bulge > 10^11 M_sun, the radio detection fraction is 77%, declining to 24% for M_bulge < 10^10 M_sun. Furthermore, we observe core-S\'ersic bulges tend to be systematically round and to possess high radio core luminosities and boxy-distorted or pure elliptical isophotes., Comment: 102 pages, 26 figures, 10 tables, accepted for publication in A&A

All 280 of the statistically-complete Palomar sample of nearby (<120 Mpc) galaxies dec > 20 degrees have been observed at 1.5 GHz as part of the LeMMINGs e-MERLIN legacy survey. Here, we present Chandra X-ray observations of the nuclei of 213 of these galaxies, including a statistically-complete sub-set of 113 galaxies in the declination range 40 degrees to 65 degrees. We observed galaxies of all optical spectral types, including 'active' galaxies (e.g., LINERs and Seyferts) and 'inactive' galaxies like HII galaxies and absorption line galaxies (ALG). The X-ray flux limit of our survey is 1.65$\times$10$^{-14}$~erg s$^{-1}$ cm$^{-2}$ (0.3$-$10 keV). We detect X-ray emission coincident within 2-arcsec of the nucleus in 150/213 galaxies, including 13/14 Seyferts, 68/77 LINERs, 13/22 ALGs and 56/100 HII galaxies, but cannot completely rule out contamination from non-AGN processes in sources with nuclear luminosities <10$^{39}$ erg s$^{-1}$. We construct an X-ray Luminosity function (XLF) and find that the local galaxy XLF, when including all AGN types, can be represented as a single power-law of slope $-0.54 \pm 0.06$. The Eddington ratio of the Seyferts is usually 2-4 decades higher than that of the LINERs, ALGs and HII galaxies, which are mostly detected with Eddington ratios <10$^{-3}$. Using [O III] line measurements and BH masses from the literature, we show that LINERs, HII galaxies and ALGs follow similar correlations to low luminosities, suggesting that some 'inactive' galaxies may harbour AGN., Comment: 20 pages + 42 pages of online supplementary material 15 figures, 7 tables (additional 150 figures and 6 tables in online supplementary material) Accepted for publication in MNRAS, 15th November 2021

What determines the nuclear radio emission in local galaxies? We combine optical [O III] line emission, robust black hole (BH) mass estimates, and high-resolution e-MERLIN 1.5-GHz data, from the LeMMINGs survey, of a statistically-complete sample of 280 nearby, optically active (LINER and Seyfert) and inactive HII and Absorption line galaxies [ALG]) galaxies. Using [O III] luminosity ($L_{\rm [O~III]}$) as a proxy for the accretion power, local galaxies follow distinct sequences in the optical-radio planes of BH activity, which suggest different origins of the nuclear radio emission for the optical classes. The 1.5-GHz radio luminosity of their parsec-scale cores ($L_{\rm core}$) is found to scale with BH mass ($M_{\rm BH}$) and [O~III] luminosity. Below $M_{\rm BH} \sim$10$^{6.5}$ M$_{\odot}$, stellar processes from non-jetted HII galaxies dominate with $L_{\rm core} \propto M_{\rm BH}^{0.61\pm0.33}$ and $L_{\rm core} \propto L_{\rm [O~III]}^{0.79\pm0.30}$. Above $M_{\rm BH} \sim$10$^{6.5}$ M$_{\odot}$, accretion-driven processes dominate with $L_{\rm core} \propto M_{\rm BH}^{1.5-1.65}$ and $L_{\rm core} \propto L_{\rm [O~III]}^{0.99-1.31}$ for active galaxies: radio-quiet/loud LINERs, Seyferts and jetted HII galaxies always display (although low) signatures of radio-emitting BH activity, with $L_{\rm 1.5\, GHz}\gtrsim$10$^{19.8}$ W Hz$^{-1}$ and $M_{\rm BH}\gtrsim10^{7}$ M$_{\odot}$, on a broad range of Eddington-scaled accretion rates ($\dot{m}$). Radio-quiet and radio-loud LINERs are powered by low-$\dot{m}$ discs launching sub-relativistic and relativistic jets, respectively. Low-power slow jets and disc/corona winds from moderately high to high-$\dot{m}$ discs account for the compact and edge-brightened jets of Seyferts, respectively. Jetted HII galaxies may host weakly active BHs. Fuel-starved BHs and recurrent activity account for ALG properties. [abridged], Comment: Accpeted for publication on MNRAS (27 pages, 9 Figures, 3 Tables, Appendix A)

We assemble a large set of 2-10 GHz radio flux density measurements and upper limits of 294 different supernovae (SNe), from the literature and our own and archival data. Only 31% of the SNe were detected. We characterize the SN lightcurves near the peak using a two-parameter model, with $t_{\rm pk}$ being the time to rise to a peak and $L_{\rm pk}$ the spectral luminosity at that peak. Over all SNe in our sample at $D<100$ Mpc, we find that $t_{\rm pk} = 10^{1.7\pm0.9}$ d, and that $L_{\rm pk} = 10^{25.5\pm1.6}$ erg s$^{-1}$ Hz$^{-1}$, and therefore that generally, 50% of SNe will have $L_{\rm pk} < 10^{25.5}$ erg s$^{-1}$ Hz$^{-1}$. These $L_{\rm pk}$ values are ~30 times lower than those for only detected SNe. Types I b/c and II (excluding IIn's) have similar mean values of $L_{\rm pk}$ but the former have a wider range, whereas Type IIn SNe have ~10 times higher values with $L_{\rm pk} = 10^{26.5\pm1.1}$ erg s$^{-1}$ Hz$^{-1}$. As for $t_{\rm pk}$, Type I b/c have $t_{\rm pk}$ of only $10^{1.1\pm0.5}$ d while Type II have $t_{\rm pk} = 10^{1.6\pm1.0}$ and Type IIn the longest timescales with $t_{\rm pk} = 10^{3.1\pm0.7}$ d. We also estimate the distribution of progenitor mass-loss rates, $\dot M$, and find the mean and standard deviation of log$_{10}(\dot M/$Msol) yr$^{-1}$ are $-5.4\pm1.2$ (assuming $v_{\rm wind}=1000$ km s$^{-1}$) for Type I~b/c SNe, and $-6.9\pm1.4$ (assuming $v_{\rm wind} = 10$ km s$^{-1}$ for Type II SNe excluding Type IIn., Comment: Accepted for publication in the Astrophysical Journal 15 Figures, 4 Tables; Full version of Table 1 in ancillary files. Minor revisions only from version 1

We present the second data release of high-resolution ($\leq0.2$ arcsec) 1.5-GHz radio images of 177 nearby galaxies from the Palomar sample, observed with the e-MERLIN array, as part of the LeMMINGs (Legacy e-MERLIN Multi-band Imaging of Nearby Galaxy Sample) survey. Together with the 103 targets of the first LeMMINGs data release, this represents a complete sample of 280 local active (LINER and Seyfert) and inactive galaxies HII galaxies and Absorption Line Galaxies, ALG). This large program is the deepest radio survey of the local Universe, $\gtrsim$10$^{17.6}$ W Hz$^{-1}$, regardless of the host and nuclear type: we detect radio emission $\gtrsim$0.25 mJy beam$^{-1}$ for 125/280 galaxies (44.6 per cent) with sizes of typically $\lesssim$100 pc. Of those 125, 106 targets show a core which coincides within 1.2 arcsec with the optical nucleus. Although we observed mostly cores, around one third of the detected galaxies features jetted morphologies. The detected radio core luminosities of the sample range between $\sim$10$^{34}$ and 10$^{40}$ erg s$^{-1}$. LINERs and Seyferts are the most luminous sources, whereas HII galaxies are the least. LINERs show FRI-like core-brightened radio structures, while Seyferts reveal the highest fraction of symmetric morphologies. The majority of HII galaxies have single radio core or complex extended structures, which probably conceal a nuclear starburst and/or a weak active nucleus (seven of them show clear jets). ALGs, which are typically found in evolved ellipticals, although the least numerous, exhibit on average the most luminous radio structures, similar to LINERs., Comment: Accepted for publication on MNRAS (48 pages, 7 figures, 10 tables, Appendix)

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5\mu$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1., Comment: 21 pages, 12 colour figures, 3 tables -- accepted for publication in MNRAS

A key characteristic of many active galactic nuclei (AGN) is their variability, but its origin is poorly understood, especially in the radio domain. Williams et al. (2017) reported a ~50 per cent increase in peak flux density of the AGN in the Seyfert galaxy NGC 4151 at 1.5 GHz with the e-MERLIN array. We present new high resolution e-MERLIN observations at 5 GHz and compare these to archival MERLIN observations to investigate the reported variability. Our new observations allow us to probe the nuclear region at a factor three times higher-resolution than the previous e-MERLIN study. We separate the core component, C4, into three separate components: C4W, C4E and X. The AGN is thought to reside in component C4W, but this component has remained constant between epochs within uncertainties. However, we find that the Eastern-most component, C4E, has increased in peak flux density from 19.35$\pm$1.10 to 37.09$\pm$1.86 mJy/beam, representing a 8.2 sigma increase on the MERLIN observations. We attribute this peak flux density increase to continued interaction between the jet and the emission line region (ELR), observed for the first time in a low-luminosity AGN such as NGC 4151. We identify discrete resolved components at 5 GHz along the jet axis, which we interpret as areas of jet-ELR interaction., Comment: 9 pages, 3 figures, 1 table, accepted to MNRAS

Submillimetre-selected galaxies (SMGs) at high redshift ($z$ $\sim$ 2) are potential host galaxies of active galactic nuclei (AGN). If the local Universe is a good guide, $\sim$ 50$\%$ of the obscured AGN amongst the SMG population could be missed even in the deepest X-ray surveys. Radio observations are insensitive to obscuration; therefore, very long baseline interferometry (VLBI) can be used as a tool to identify AGN in obscured systems. A well-established upper limit to the brightness temperature of 10$^5$ K exists in star-forming systems, thus VLBI observations can distinguish AGN from star-forming systems via brightness temperature measurements. We present 1.6 GHz European VLBI Network (EVN) observations of four SMGs (with measured redshifts) to search for evidence of compact radio components associated with AGN cores. For two of the sources, e-MERLIN images are also presented. Out of the four SMGs observed, we detect one source, J123555.14, that has an integrated EVN flux density of 201 $\pm$ 15.2 $\mu$Jy, corresponding to a brightness temperature of 5.2 $\pm$ 0.7 $\times$ 10$^5$ K. We therefore identify that the radio emission from J123555.14 is associated with an AGN. We do not detect compact radio emission from a possible AGN in the remaining sources (J123600.10, J131225.73, and J163650.43). In the case of J131225.73, this is particularly surprising, and the data suggest that this may be an extended, jet-dominated AGN that is resolved by VLBI. Since the morphology of the faint radio source population is still largely unknown at these scales, it is possible that with a $\sim$ 10 mas resolution, VLBI misses (or resolves) many radio AGN extended on kiloparsec scales., Comment: 11 pages, 5 figures

We studied the nearby, interacting galaxy NGC 5195 (M51b) in the radio, optical and X-ray bands. We mapped the extended, low-surface-brightness features of its radio-continuum emission; determined the energy content of its complex structure of shock-ionized gas; constrained the current activity level of its supermassive nuclear black hole. In particular, we combined data from the European Very Long Baseline Interferometry Network (~1-pc scale), from our new e-MERLIN observations (~10-pc scale), and from the Very Large Array (~100-1000-pc scale), to obtain a global picture of energy injection in this galaxy. We put an upper limit to the luminosity of the (undetected) flat-spectrum radio core. We find steep-spectrum, extended emission within 10 pc of the nuclear position, consistent with optically-thin synchrotron emission from nuclear star formation or from an outflow powered by an active galactic nucleus (AGN). A linear spur of radio emission juts out of the nuclear source towards the kpc-scale arcs (detected in radio, Halpha and X-ray bands). From the size, shock velocity, and Balmer line luminosity of the kpc-scale bubble, we estimate that it was inflated by a long-term-average mechanical power ~3-6 x 10^{41} erg/s over the last 3-6 Myr. This is an order of magnitude more power than can be provided by the current level of star formation, and by the current accretion power of the supermassive black hole. We argue that a jet-inflated bubble scenario associated with previous episodes of AGN activity is the most likely explanation for the kpc-scale structures., Comment: 16 pages, 7 figures. Accepted by MNRAS on 2017 Feb 12

We present the first data release of high-resolution ($\leq0.2$ arcsec) 1.5-GHz radio images of 103 nearby galaxies from the Palomar sample, observed with the eMERLIN array, as part of the LeMMINGs survey. This sample includes galaxies which are active (LINER and Seyfert) and quiescent (HII galaxies and Absorption line galaxies, ALG), which are reclassified based upon revised emission-line diagrams. We detect radio emission $\gtrsim$ 0.2 mJy for 47/103 galaxies (22/34 for LINERS, 4/4 for Seyferts, 16/51 for HII galaxies and 5/14 for ALGs) with radio sizes typically of $\lesssim$100 pc. We identify the radio core position within the radio structures for 41 sources. Half of the sample shows jetted morphologies. The remaining half shows single radio cores or complex morphologies. LINERs show radio structures more core-brightened than Seyferts. Radio luminosities of the sample range from 10$^{32}$ to 10$^{40}$ erg s$^{-1}$: LINERs and HII galaxies show the highest and the lowest radio powers respectively, while ALGs and Seyferts have intermediate luminosities. We find that radio core luminosities correlate with black hole (BH) mass down to $\sim$10$^{7}$ M$_{\odot}$, but a break emerges at lower masses. Using [O III] line luminosity as a proxy for the accretion luminosity, active nuclei and jetted HII galaxies follow an optical fundamental plane of BH activity, suggesting a common disc-jet relationship. In conclusion, LINER nuclei are the scaled-down version of FR I radio galaxies; Seyferts show less collimated jets; HII galaxies may host weak active BHs and/or nuclear star-forming cores; and recurrent BH activity may account for ALG properties., Comment: accepted for publication on MNRAS (45 pages, 10 figures with reduced image quality)

V959 Mon (Nova Mon 2012) was first detected by the Fermi Large Area Telescope in June 2012, as a transient gamma-ray source. Subsequent optical observations showed that this gamma-ray emission was due to a classical nova explosion. Multi-frequency observations of V959 Mon with the VLA between June and September 2012 revealed dramatic brightening, and a spectrum that steepened with increasing frequency. High resolution radio images of V959 Mon using e-MERLIN are presented here, at six epochs between September 2012 and February 2014 which show morphological evolution of the source. While early e-MERLIN observations of V959 Mon show an east-west elonga- tion in the ejecta morphology, subsequent observations suggest that the ejecta become elongated in the north-south direction. Our high-resolution observations of this sur- prising evolution in the structure of V959 Mon can assist us in further understanding the behaviour and morphology of nova ejecta.

We present high-sensitivity eMERLIN radio images of the Seyfert galaxy NGC 4151 at 1.5 GHz. We compare the new eMERLIN images to those from archival MERLIN observations in 1993 to determine the change in jet morphology in the 22 years between observations. We report an increase by almost a factor of 2 in the peak flux density of the central core component, C4, thought to host the black hole, but a probable decrease in some other components, possibly due to adiabatic expansion. The core flux increase indicates an AGN which is currently active and feeding the jet. We detect no significant motion in 22 years between C4 and the component C3, which is unresolved in the eMERLIN image. We present a spectral index image made within the 512 MHz band of the 1.5 GHz observations. The spectrum of the core, C4, is flatter than that of other components further out in the jet. We use HST emission line images (H$\alpha$, [O III] and [O II]) to study the connection between the jet and the emission line region. Based on the changing emission line ratios away from the core and comparison with the eMERLIN radio jet, we conclude that photoionisation from the central AGN is responsible for the observed emission line properties further than 4" (360 pc) from the core, C4. Within this region, several evidences (radio-line co-spatiality, low [O III]/H$\alpha$ and estimated fast shocks) suggest additional ionisation from the jet., Comment: 13 pages, 8 figures, 2 tables. Accepted to MNRAS: MN-17-2603-MJ.R1

We present high-resolution e-MERLIN radio continuum maps of the Dwarf Irregular galaxy IC10 at 1.5GHz and 5GHz. We detect 11 compact sources at 1.5GHz, 5 of which have complementary detections at 5GHz. We classify 3 extended sources as compact H{\sc ii} regions within IC10, 5 sources as contaminating background galaxies and identify 3 sources which require additional observations to classify. We do not expect that any of these 3 sources are Supernova Remnants as they will likely be resolved out at the assumed distance of IC10 (0.7Mpc). We correct integrated flux densities of IC10 from the literature for contamination by unrelated background sources and obtain updated flux density measurements of $354\pm11$\,mJy at 1.5GHz and $199\pm9$\,mJy at 4.85GHz. The background contamination does not contribute significantly to the overall radio emission from IC10, so previous analysis concerning its integrated radio properties remain valid., Comment: 15 pages, 7 Figures, 5 Tables, Accepted for publication in MNRAS

We present a radio light curve of supernova (SN) 2014C taken with the Arcminute Microkelvin Imager (AMI) Large Array at 15.7 GHz. Optical observations presented by Milisavljevic et al. demonstrated that SN 2014C metamorphosed from a stripped-envelope Type Ib SN into a strongly interacting Type IIn SN within 1 year. The AMI light curve clearly shows two distinct radio peaks, the second being a factor of 4 times more luminous than the first peak. This double bump morphology indicates two distinct phases of mass-loss from the progenitor star with the transition between density regimes occurring at 100-200 days. This reinforces the interpretation that SN 2014C exploded in a low density region before encountering a dense Hydrogen-rich shell of circumstellar material that was likely ejected by the progenitor prior to the explosion. The AMI flux measurements of the first light curve bump are the only reported observations taken within ~50 to ~125 days post-explosion, before the blast-wave encountered the Hydrogen shell. Simplistic synchrotron self-absorption (SSA) and free-free absorption (FFA) modelling suggest that some physical properties of SN 2014C, such as the mass-loss rate, are consistent with the properties of other Type Ibc and IIn SNe. However, our single frequency data does not allow us to distinguish between these two models, which implies they are likely too simplistic to describe the complex environment surrounding this event. Lastly, we present the precise radio location of SN 2014C obtained with eMERLIN, which will be useful for future VLBI observations of the SN., Comment: 19 pages, 11 figures, 3 tables, accepted for publication in MNRAS

We present results from sensitive, multi-epoch NuSTAR observations of the late-type star-forming galaxy M83 (d=4.6 Mpc), which is the first investigation to spatially resolve the hard (E>10 keV) X-ray emission of this galaxy. The nuclear region and ~ 20 off-nuclear point sources, including a previously discovered ultraluminous X-ray (ULX) source, are detected in our NuSTAR observations. The X-ray hardnesses and luminosities of the majority of the point sources are consistent with hard X-ray sources resolved in the starburst galaxy NGC 253. We infer that the hard X-ray emission is most likely dominated by intermediate accretion state black hole binaries and neutron star low-mass X-ray binaries (Z-sources). We construct the X-ray binary luminosity function (XLF) in the NuSTAR band for an extragalactic environment for the first time. The M83 XLF has a steeper XLF than the X-ray binary XLF in NGC 253, consistent with previous measurements by Chandra at softer X-ray energies. The NuSTAR integrated galaxy spectrum of M83 drops quickly above 10 keV, which is also seen in the starburst galaxies NGC253, NGC 3310 and NGC 3256. The NuSTAR observations constrain any AGN to be either highly obscured or to have an extremely low luminosity of $_{\sim}^<$10$^{38}$ erg/s (10-30 keV), implying it is emitting at a very low Eddington ratio. An X-ray point source consistent with the location of the nuclear star cluster with an X-ray luminosity of a few times 10$^{38}$ erg/s may be a low-luminosity AGN but is more consistent with being an X-ray binary., Comment: Accepted for publication in ApJ (25 pages, 17 figures)

The SKA will be a transformational instrument in the study of our local Universe. In particular, by virtue of its high sensitivity (both to point sources and diffuse low surface brightness emission), angular resolution and the frequency ranges covered, the SKA will undertake a very wide range of astrophysical research in the field of nearby galaxies. By surveying vast numbers of nearby galaxies of all types with $\mu$Jy sensitivity and sub-arcsecond angular resolutions at radio wavelengths, the SKA will provide the cornerstone of our understanding of star-formation and accretion activity in the local Universe. In this chapter we outline the key continuum and molecular line science areas where the SKA, both during phase-1 and when it becomes the full SKA, will have a significant scientific impact., Comment: 23 pages, 4 figures, to appear as part of 'Continuum Science' in Proceedings 'Advancing Astrophysics with the Square Kilometre Array' PoS(AASKA14)070

We report on simultaneous observations of the local starburst system Arp 299 with NuSTAR and Chandra, which provides the first resolved images of this galaxy up to energies of ~ 45 keV. Fitting the 3-40 keV spectrum reveals a column density of $N_{\rm H}$ ~ 4 x10^{24} cm^{-2}, characteristic of a Compton-thick AGN, and a 10-30 keV luminosity of 1.2x 10^{43} ergs s^{-1}. The hard X-rays detected by NuSTAR above 10 keV are centered on the western nucleus, Arp 299-B, which previous X-ray observations have shown to be the primary source of neutral Fe-K emission. Other X-ray sources, including Arp 299-A, the eastern nucleus which is also thought to harbor an AGN, as well as X-ray binaries, contribute $\lesssim 10%$ to the 10-20 keV emission from the Arp 299 system. The lack of significant emission above 10 keV other than that attributed to Arp 299-B suggests that: a) any AGN in Arp 299-A must be heavily obscured ($N_{\rm H}$ > 10^{24} cm^{-2}) or have a much lower luminosity than Arp 299-B and b) the extranuclear X-ray binaries have spectra that cut-off above ~10 keV. Such soft spectra are characteristic of ultraluminous X-ray (ULX) sources observed to date by NuSTAR., Comment: 9 pages; accepted for publication in Astrophysical Journal

We report deep EVN and eMERLIN observations of the Type Ia SN 2014J in the nearby galaxy M 82. Our observations represent, together with JVLA observations of SNe 2011fe and 2014J, the most sensitive radio studies of Type Ia SNe ever. By combining data and a proper modeling of the radio emission, we constrain the mass-loss rate from the progenitor system of SN 2014J to $\dot{M} \lesssim 7.0\times 10^{-10}\, {\rm M_{\odot}\, yr^{-1}}$ (3-$\sigma$; for a wind speed of $100\, {\rm km s^{-1}}$). If the medium around the supernova is uniform, then $n_{\rm ISM} \lesssim 1.3 {\rm cm^3}$ (3-$\sigma$), which is the most stringent limit for the (uniform) density around a Type Ia SN. Our deep upper limits favor a double-degenerate (DD) scenario--involving two WD stars--for the progenitor system of SN 2014J, as such systems have less circumstellar gas than our upper limits. By contrast, most single-degenerate (SD) scenarios, i.e., the wide family of progenitor systems where a red giant, main-sequence, or sub-giant star donates mass to a exploding WD, are ruled out by our observations. Our estimates on the limits to the gas density surrounding SN 2011fe, using the flux density limits from Chomiuk et al. (2012), agree well with their results. Although we discuss possibilities for a SD scenario to pass observational tests, as well as uncertainties in the modeling of the radio emission, the evidence from SNe 2011fe and 2014J points in the direction of a DD scenario for both., Comment: 10 pages, 5 figures. Accepted for publication on The Astrophysical Journal

We present results from three nearly simultaneous NuSTAR and Chandra monitoring observations between 2012 Sep 2 and 2012 Nov 16 of local star-forming galaxy NGC 253. The 3-40 keV NuSTAR intensity of the inner 20 arcsec (~400 pc) nuclear region varied by a factor of ~2 across the three monitoring observations. The Chandra data reveal that the nuclear region contains three bright X-ray sources, including a luminous (L2-10 keV ~ few x 10^39 erg/s) point source ~1 arcsec from the dynamical center of the galaxy (within the 3sigma positional uncertainty of the dynamical center); this source drives the overall variability of the nuclear region at energies >3 keV. We make use of the variability to measure the spectra of this single hard X-ray source when it was in bright states. The spectra are well described by an absorbed (NH ~ 1.6 x 10^23 cm^-2) broken power-law model with spectral slopes and break energies that are typical of ultraluminous X-ray sources (ULXs), but not AGN. A previous Chandra observation in 2003 showed a hard X-ray point source of similar luminosity to the 2012 source that was also near the dynamical center (~0.4 arcsec); however, this source was offset from the 2012 source position by ~1 arcsec. We show that the probability of the 2003 and 2012 hard X-ray sources being unrelated is >>99.99% based on the Chandra spatial localizations. Interestingly, the Chandra spectrum of the 2003 source (3-8 keV) is shallower in slope than that of the 2012 hard X-ray source. Its proximity to the dynamical center and harder Chandra spectrum indicate that the 2003 source is a better AGN candidate than any of the sources detected in our 2012 campaign; however, we were unable to rule out a ULX nature for this source. Future NuSTAR and Chandra monitoring would be well equipped to break the degeneracy between the AGN and ULX nature of the 2003 source, if again caught in a high state., Comment: Accepted for publication in ApJ (7 pages, 5 figures)

This paper presents the results of the 2009-2010 monitoring sessions of the starburst galaxy M82, obtained with the Multi-Element Radio-Linked Interferometer Network (MERLIN) at 5GHz and e-MERLIN at 6GHz. Combining several 5GHz MERLIN epochs to form a map with 33.0 uJy/bm noise level, 52 discrete sources, mostly supernova remnants and HII regions, are identified. These include three objects which were not detected in the 2002 5GHz MERLIN monitoring session: supernova SN2008iz, the transient source 43.78+59.3, and a new supernova remnant shell. Flux density variations, in the long (1981 to 2010), medium (2002 to 2010) and short (2009 to 2010) term, are investigated. We find that flux densities of SNRs in M82 stay constant in most of the sample (~95%). In addition, aside from SN2008iz and the well-known variable source 41.95+57.5, two sources display short and medium term variations over the period 2009-2010. These sources being among the most compact SNR in M82, these flux density variations could be due to changes in the circumstellar and interstellar medium in which the shocks travel., Comment: 15 pages, 10 figures

The high star formation rates of luminous infrared galaxies (LIRGs) make them ideal places for core-collapse supernova (CCSN) searches. At radio frequencies, free from dust extinction, it is possible to detect compact components within the innermost LIRG nuclear regions, such as SNe and SN remnants, as well as AGN buried deep in the LIRG nuclei. We studied the LIRG IC883 aiming at: (i) investigating its (circum-)nuclear regions using the e-EVN at 5GHz, and e-MERLIN at 6.9GHz, complemented by archival VLBI data; (ii) detecting at radio frequencies the two recently reported circumnuclear SNe 2010cu and 2011hi, which were discovered by near-IR (NIR) adaptive optics observations of IC883; and (iii) further investigating the nature of SN2011hi at NIR by means of observations with Gemini-North. The circumnuclear regions traced by e-MERLIN at 6.9GHz have an extension of ~1kpc, and show a striking double-sided structure, which very likely corresponds to a warped rotating ring, in agreement with previous studies. Our e-EVN observations at 5GHz and complementary archival VLBI data at 5GHz and 8.4GHz, reveal the presence of various milliarcsec compact components in the nucleus of IC883. A single compact source, an AGN candidate, dominates the emission at both nuclear and circumnuclear scales, as imaged with the e-EVN and e-MERLIN, respectively. The other milliarcsec components are very suggestive of ongoing nuclear CCSN activity. Our e-EVN observations also resulted in upper limits to the radio luminosity of the two SNe in IC883 recently discovered at NIR. We refine the classification of SN2011hi as a Type IIP SN according to our latest Gemini-North epoch from 2012, in agreement with a low-luminosity radio SN nature. We estimate a CCSN rate lower limit of 1.1_{-0.6}^{+1.3} yr^{-1} for the entire galaxy, based on three nuclear radio SNe and the circumnuclear SNe 2010cu and 2011hi. (abridged), Comment: 9 pages, 5 figures and 2 tables. Accepted for publication in A&A

We report on the VLBI detection of supernova SN2011dh at 22GHz using a subset of the EVN array. The observations took place 14 days after the discovery of the supernova, thus resulting in a VLBI image of the youngest radio-loud supernova ever. We provide revised coordinates for the supernova with milli-arcsecond precision, linked to the ICRF. The recovered flux density is a factor 2 below the EVLA flux density reported by other authors at the same frequency and epoch of our observations. This discrepancy could be due to extended emission detected with the EVLA or to calibration problems in the VLBI and/or EVLA observations., Comment: Letter. Accepted in A&A

From a combination of MERLIN (Multi-Element Radio-Linked Interferometer Network) and global VLBI (Very Long Baseline Interferometry) observations of the starburst galaxy M82, images of 36 discrete sources at resolutions ranging from ~3 to ~80 mas at 1.7 GHz are presented. Of these 36 sources, 32 are identified as supernova remnants, 2 are HII regions, and 3 remain unclassified. Sizes, flux densities and radio brightnesses are given for all of the detected sources. Additionally, global VLBI only data from this project are used to image four of the most compact radio sources. These data provide a fifth epoch of VLBI observations of these sources, covering a 19-yr time-line. In particular, the continued expansion of one of the youngest supernova remnants, 43.31+59.3 is discussed. The deceleration parameter is a power-law index used to represent the time evolution of the size of a supernova remnant. For the source 43.31+59.3, a lower limit to the deceleration parameter is calculated to be 0.53+/-0.06, based on a lower limit of the age of this source., Comment: 31 pages, 12 figures, 7 tables

A faint new radio source has been detected in the nuclear region of the starburst galaxy M82 using MERLIN radio observations designed to monitor the flux density evolution of the recent bright supernova SN2008iz. This new source was initially identified in observations made between 1-5th May 2009 but had not been present in observations made one week earlier, or in any previous observations of M82. In this paper we report the discovery of this new source and monitoring of its evolution over its first 9 months of existence. The true nature of this new source remains unclear, and we discuss whether this source may be an unusual and faint supernova, a supermassive blackhole associated with the nucleus of M82, or intriguingly the first detection of radio emission from an extragalactic microquasar., Comment: 6 pages, 3 figures (2 colour), MNRAS letters accepted

The class of type Ic supernovae have drawn increasing attention since 1998 owing to their sparse association (only four so far) with long duration gamma-ray bursts. Although both phenomena originate from the core collapse of a massive star, supernovae emit mostly at optical wavelengths, whereas GRBs emit mostly in soft gamma-rays or hard X-rays. Though the GRB central engine generates ultra-relativistic jets, which beam the early emission into a narrow cone, no relativistic outflows have hitherto been found in type Ib/c supernovae explosions, despite theoretical expectations and searches. Here we report radio (interferometric) observations that reveal a mildly relativistic expansion in a nearby type Ic supernova, SN 2007gr. Using two observational epochs 60 days apart, we detect expansion of the source and establish a conservative lower limit for the average apparent expansion velocity of 0.6c. Independently, a second mildly relativistic supernova has been reported. Contrary to the radio data, optical observations of SN 2007gr indicate a typical type Ic supernova with ejecta velocities ~6000 km/s, much lower than in GRB-associated supernovae. We conclude that in SN 2007gr a small fraction of the ejecta produced a low-energy mildly relativistic bipolar radio jet, while the bulk of the ejecta were slower and, as shown by optical spectro-polarimetry, mildly aspherical., Comment: Nature, 10 pages (including supplementary material), 2 figures

Using the VLA, a series of high velocity resolution observations have been made of the M82 starburst at 1.6 GHz. These observations follow up on previous studies of the main line OH maser emission in the central kiloparsec of this starburst region, but with far greater velocity resolution, showing significant velocity structure in some of the maser spots for the first time. A total of thirteen masers were detected, including all but one of the previously known sources. While some of these masers are still unresolved in velocity, these new results clearly show velocity structure in spectra from several of the maser regions. Position-velocity plots show good agreement with the distribution of H{\sc i} including interesting velocity structure on the blueward feature in the west of the starburst which traces the velocity distribution seen in the ionised gas., Comment: MNRAS in press. 15 pages, 9 figures

We report on phase-referenced 23 GHz Very-Long-Baseline-Interferometry (VLBI) observations of the type IIb supernova SN 2008ax, made with the Very Long Baseline Array (VLBA) on 2 April 2008 (33 days after explosion). These observations resulted in a marginal detection of the supernova. The total flux density recovered from our VLBI image is 0.8$\pm$0.3 mJy (one standard deviation). As it appears, the structure may be interpreted as either a core-jet or a double source. However, the supernova structure could be somewhat confused with a possible close by noise peak. In such a case, the recovered flux density would decrease to 0.48$\pm$0.12 mJy, compatible with the flux densities measured with the VLA at epochs close in time to our VLBI observations. The lowest average expansion velocities derived from our observations are $(1.90 \pm 0.30) \times 10^5$ km s$^{-1}$ (case of a double source) and $(5.2 \pm 1.3) \times 10^4$ km s$^{-1}$ (taking the weaker source component as a spurious, close by, noise peak, which is the more likely interpretation). These velocities are 7.3 and 2 times higher, respectively, than the maximum ejecta velocity inferred from optical-line observations., Comment: 4 pages, 3 figures. Accepted in A&A on 24/03/2009

The results of an extremely deep, 8-day long observation of the central kpc of the nearby starburst galaxy M82 using MERLIN (Multi-Element Radio Linked Interferometer Network) at 5 GHz are presented. The 17E-06 Jy/beam, rms noise level in the naturally weighted image make it the most sensitive high resolution radio image of M82 made to date. Over 50 discrete sources are detected, the majority of which are supernova remnants, but with 13 identified as HII regions. Sizes, flux densities and radio brightnesses are given for all of the detected sources, which are all well resolved with a majority showing shell or partial shell structures. Those sources within the sample which are supernova remnants have diameters ranging from 0.3 to 6.7 pc, with a mean size of 2.9 pc. From a comparison with previous MERLIN 5 GHz observations made in July 1992, which gives a 9.75 year timeline, it has been possible to measure the expansion velocities of ten of the more compact sources, eight of which have not been measured before. These derived expansion velocities range between 2200 and 10500 km/s., Comment: 34 pages, 10 figures. Accepted by MNRAS

A study of the distribution of OH gas in the central region of the nearby active starburst galaxy M82 has confirmed two previously known bright masers and revealed several new main line masers. Three of these are seen only at 1665 MHz, one is detected only at 1667 MHz, while the rest are detected in both lines. Observations covering both the 1665 and 1667 MHz lines, conducted with both the Very Large Array (VLA) and the Multi-Element Radio Linked Interferometer Network (MERLIN), have been used to accurately measure the positions and velocities of these features. This has allowed a comparison with catalogued continuum features in the starburst such as HII regions and supernova remnants, as well as known water and satellite line OH masers. Most of the main line masers appear to be associated with known HII regions although the two detected only at 1665 MHz are seen along the same line of sight as known supernova remnants., Comment: MNRAS accepted. 16 pages, 13 figures

We report on 8.4GHz Very Long Baseline Interferometry (VLBI) observations of the type II-P supernova SN2004et in the spiral galaxy NGC 6946, made on 20 February 2005 (151 days after explosion). The Very Large Array (VLA) flux density was 1.23$\pm$0.07 mJy, corresponding to an isotropic luminosity at 8.4GHz of (4.45$\pm$0.3)$\times10^{25}$ erg s$^{-1}$ Hz$^{-1}$ and a brightness temperature of (1.3$\pm$0.3)$\times10^{8}$ K. We also provide an improved source position, accurate to about 0.5 mas in each coordinate. The VLBI image shows a clear asymmetry. From model fitting of the size of the radio emission, we estimate a minimum expansion velocity of 15,700$\pm$2,000 km s$^{-1}$. This velocity is more than twice the expected mean expansion velocity estimated from a synchrotron self-absorbed emission model, thus suggesting that synchrotron self-absorption is not relevant for this supernova. With the benefit of an optical spectrum obtained 12 days after explosion, we favor an emission model which consists of two hot spots on an underlying expanding shell of width comparable to that of SN 1993J., Comment: Accepted for publication in A&A (22/05/07)

We present the results of a second epoch of 18cm global Very Long-Baseline Interferometry (VLBI) observations, taken on 23 February 2001, of the central kiloparsec of the nearby starburst galaxy Messier 82. These observations further investigate the structural and flux evolution of the most compact radio sources in the central region of M82. The two most compact radio objects in M82 have been investigated (41.95+575 and 43.31+592). Using this recent epoch of data in comparison with our previous global VLBI observations and two earlier epochs of European VLBI Network observations we measure expansion velocities in the range of 1500-2000km/s for 41.95+575, and 9000-11000km/s for 43.31+592 using various independent methods. In each case the measured remnant expansion velocities are significantly larger than the canonical expansion velocity (500km/s) of supernova remnants within M82 predicted from theoretical models. In this paper we discuss the implications of these measured expansion velocities with respect to the high density environment that the SNR are expected to reside in within the centre of the M82 starburst., Comment: Accepted for publication in MNRAS, 9 pages, 8 figures

Supernova 2004dj in the nearby spiral galaxy NGC2403 was detected optically in July 2004. Peaking at a magnitude of 11.2, this is the brightest supernova detected for several years. Here we present Multi-Element Radio Linked Interferometer Network (MERLIN) observations of this source, made over a four month period, which give a position of R.A. = 07h37m17.044s, Dec =+65deg35'57.84" (J2000.0). We also present a well-sampled 5 GHz light curve covering the period from 5 August to 2 December 2004. With the exception of the unusual and very close SN 1987A, these observations represent the first detailed radio light curve for the prompt emission from a Type II-P supernova., Comment: (1) Jodrell Bank Observatory (2) University of Valencia (3) University of Sheffield 6 pages, 1 figure. To appear in ApJ letters

Several new main line OH masers have been detected in the nearby starburst galaxy M82. Eight masers have been detected to 5 sigma, six of which are new detections. Observations covering both the 1665 and 1667 MHz lines with both the Very Large Array (VLA) and the Multi-Element Radio Linked Interferometer Network (MERLIN) have been used to accurately measure the positions and velocities of these features. Following analysis of the data, another six objects below 5 sigma, but with velocities consistent with the distribution inferred from the more certain detections, have been detected. These are classified as possible detections., Comment: 2 pages. 2 figures. Proceedings of the 7th European VLBI Network Symposium held in Toledo, Spain on October 12-15, 2004. Editors: R. Bachiller, F. Colomer, J.-F. Desmurs, P. de Vicente (Observatorio Astronomico Nacional), p. 151-152. Needs evn2004.cls

As part of an on-going radio supernova monitoring program, we have discovered a variable, compact steep spectrum radio source ~65 arcsec (~4 kpc) from the centre of the starburst galaxy NGC 3310. If the source is at the distance of NGC 3310, then its 5 GHz luminosity is ~3 x 10^{19} WHz^-1. The source luminosity, together with its variability characteristics, compact structure (<17 mas) and its association with a group of HII regions, leads us to propose that it is a previously uncatalogued type II radio supernova. A search of archival data also shows an associated X-ray source with a luminosity similar to known radio supernova., Comment: 5 pages, 3 figures. Accepted by MNRAS

A new two-image gravitational lens system has been discovered as a result of the Cosmic Lens All-Sky Survey (CLASS). Radio observations with the VLA, MERLIN and the VLBA at increasingly higher resolutions all show two components with a flux density ratio of ~7:1 and a separation of 1.34". Both components are compact and have the same spectral index. Followup observations made with the VLA at 8.4 GHz show evidence of a feature to the south-east of the brighter component and a corresponding extension of the weaker component to the north-west. Optical observations with the WHT show ~1.7" extended emission aligned in approximately the same direction as the separation between the radio components with an R-band magnitude of 21.8 +/- 0.4., Comment: 6 pages, 5 figures, accepted by MNRAS