Source author record

J. B. R. Oonk

J. B. R. Oonk appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

46works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2022arXiv

Deep sub-arcsecond widefield imaging of the Lockman Hole field at 144 MHz

High quality low-frequency radio surveys have the promise of advancing our understanding of many important topics in astrophysics, including the life cycle of active galactic nuclei (AGN), particle acceleration processes in jets, the history of star formation, and exoplanet magnetospheres. Currently leading low-frequency surveys reach an angular resolution of a few arcseconds. However, this resolution is not yet sufficient to study the more compact and distant sources in detail. Sub-arcsecond resolution is therefore the next milestone in advancing these fields. The biggest challenge at low radio frequencies is the ionosphere. If not adequately corrected for, ionospheric seeing blurs the images to arcsecond or even arcminute scales. Additionally, the required image size to map the degree-scale field of view of low-frequency radio telescopes at this resolution is far greater than what typical soft- and hardware is currently capable of handling. Here we present for the first time (to the best of our knowledge) widefield sub-arcsecond imaging at low radio frequencies. We derive ionospheric corrections in a few dozen individual directions and apply those during imaging efficiently using a recently developed imaging algorithm (arXiv:1407.1943, arXiv:1909.07226). We demonstrate our method by applying it to an eight hour observation of the International LOw Frequency ARray (LOFAR) Telescope (ILT) (arXiv:1305.3550). Doing so we have made a sensitive $7.4\ \mathrm{deg}^2$ $144\ \mathrm{MHz}$ map at a resolution of $0.3''$ reaching $25\ μ\mathrm{Jy\ beam}^{-1}$ near the phase centre. The estimated $250,000$ core hours used to produce this image, fit comfortably in the budget of available computing facilities. This result will enable future mapping of the entire northern low-frequency sky at sub-arcsecond resolution.

preprint2022arXiv

Filamentary structures of ionized gas in Cygnus X

Ionized gas probes the influence of massive stars on their environment. The Cygnus X region (d~1.5 kpc) is one of the most massive star forming complexes in our Galaxy, in which the Cyg OB2 association (age of 3-5 Myr and stellar mass $2 \times 10^{4}$ M$_{\odot}$) has a dominant influence. We observe the Cygnus X region at 148 MHz using the Low Frequency Array (LOFAR) and take into account short-spacing information during image deconvolution. Together with data from the Canadian Galactic Plane Survey, we investigate the morphology, distribution, and physical conditions of low-density ionized gas in a $4^{\circ} \times 4^{\circ}$ (100 pc $\times$ 100 pc) region at a resolution of 2' (0.9 pc). The Galactic radio emission in the region analyzed is almost entirely thermal (free-free) at 148 MHz, with emission measures of $10^3 < EM~{\rm[pc~cm^{-6}]} < 10^6$. As filamentary structure is a prominent feature of the emission, we use DisPerSE and FilChap to identify filamentary ridges and characterize their radial ($EM$) profiles. The distribution of radial profiles has a characteristic width of 4.3 pc and a power-law distribution ($β= -1.8 \pm 0.1$) in peak $EM$ down to our completeness limit of 4200 pc cm$^{-6}$. The electron densities of the filamentary structure range from $10 < n_e~{\rm[cm^{-3}]} < 400$ with a median value of 35 cm$^{-3}$, remarkably similar to [N II] surveys of ionized gas. Cyg OB2 may ionize at most two-thirds of the total ionized gas and the ionized gas in filaments. More than half of the filamentary structures are likely photoevaporating surfaces flowing into a surrounding diffuse (~5 cm$^{-3}$) medium. However, this is likely not the case for all ionized gas ridges. A characteristic width in the distribution of ionized gas points to the stellar winds of Cyg OB2 creating a fraction of the ionized filaments through swept-up ionized gas or dissipated turbulence.

preprint2022arXiv

The LOFAR Two-metre Sky Survey -- V. Second data release

In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$μ$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$μ$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data.

preprint2021arXiv

A parsec-scale faint jet in the nearby changing-look Seyfert galaxy Mrk 590

Broad Balmer emission lines in active galactic nuclei (AGN) may display dramatic changes in amplitude, even disappearance and re-appearance in some sources. As a nearby galaxy at a redshift of z = 0.0264, Mrk 590 suffered such a cycle of Seyfert type changes between 2006 and 2017. Over the last fifty years, Mrk 590 also underwent a powerful continuum outburst and a slow fading from X-rays to radio wavelengths with a peak bolometric luminosity reaching about ten per cent of the Eddington luminosity. To track its past accretion and ejection activity, we performed very long baseline interferometry (VLBI) observations with the European VLBI Network (EVN) at 1.6 GHz in 2015. The EVN observations reveal a faint (~1.7 mJy) radio jet extending up to ~2.8 mas (projected scale ~1.4 pc) toward north, and probably resulting from the very intensive AGN activity. To date, such a parsec-scale jet is rarely seen in the known changing-look AGN. The finding of the faint jet provides further strong support for variable accretion as the origin of the type changes in Mrk 590.

preprint2021arXiv

Sub-arcsecond imaging with the International LOFAR Telescope I. Foundational calibration strategy and pipeline

[abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $μ$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $μ$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit.

preprint2020arXiv

A molecular absorption line survey toward the AGN of Hydra-A

We present Atacama Large Millimeter/submillimeter Array observations of the brightest cluster galaxy Hydra-A, a nearby ($z=0.054$) giant elliptical galaxy with powerful and extended radio jets. The observations reveal CO(1-0), CO(2-1), $^{13}$CO(2-1), CN(2-1), SiO(5-4), HCO$^{+}$(1-0), HCO$^{+}$(2-1), HCN(1-0), HCN(2-1), HNC(1-0) and H$_{2}$CO(3-2) absorption lines against the galaxy's bright and compact active galactic nucleus. These absorption features are due to at least 12 individual molecular clouds which lie close to the centre of the galaxy and have velocities of approximately $-50$ to $+10$ km/s relative to its recession velocity, where positive values correspond to inward motion. The absorption profiles are evidence of a clumpy interstellar medium within brightest cluster galaxies composed of clouds with similar column densities, velocity dispersions and excitation temperatures to those found at radii of several kpc in the Milky Way. We also show potential variation in a $\sim 10$ km/s wide section of the absorption profile over a two year timescale, most likely caused by relativistic motions in the hot spots of the continuum source which change the background illumination of the absorbing clouds.

preprint2020arXiv

LOFAR 144-MHz follow-up observations of GW170817

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$σ$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $α^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations.

preprint2020arXiv

Pioneering the Exascale era with Astronomy

SURF is the collaborative ICT organisation for Dutch education & research and coordinates the national e-infrastructure (surf.nl). To accelerate scientific discovery, SURF invests in, operates and explores high-end IT solutions for and with researchers in the Netherlands. In this ADASS 2019 contribution we present our latest developments in high performance and high throughput cloud computing. These developments are particularly relevant for Astronomy, as this science domain: (i) generates large (Petabyte sized) data collections, (ii) uses rapid release and deployment schemes for their software, (iii) requires flexible and interactive test and staging environments and (iv) needs to execute complex workflows on diverse data structures. We highlight our new OpenStack-based, cloud infrastructure layer and focus on Spider, our new addition to SURF's high throughput data processing platforms.

preprint2020arXiv

Tied-array holography with LOFAR

A radio interferometer uses time delays to maximize its response to radiation coming from a particular direction. These time delays compensate for differences in the time of arrival of the wavefront at the different elements of the interferometer, and for delays in the instrument's signal chain. If the radio interferometer is operated as a phased array (tied array), the time delays cannot be accounted for after an observation, so they must be determined in advance. Our aim is to characterize the time delays between the stations in the core of the LOw Frequency ARray (LOFAR). We used radio holography to determine the time delays for the core stations of LOFAR (innermost 3.5 km). Using the multibeaming capability of LOFAR we map the voltage beam faster than with a raster scan, while simultaneously calibrating the observed beam continuously. For short radio holographic observations (60 s and 600 s) of 3C196, 3C147, and 3C48 we are able to derive time delays with errors of less than one nanosecond. After applying the derived time delays to the beamformer, the beam shows residuals of less than $20\%$ with respect to the theoretical beam shape. Tied-array holography could be a way towards semi-real-time beam calibration for the Square Kilometer Array.

preprint2019arXiv

Searching for the largest bound atoms in space

(abridged) Radio recombination lines (RRLs) at frequencies $ν$ < 250 MHz trace the cold, diffuse phase of the ISM. Next generation low frequency interferometers, such as LOFAR, MWA and the future SKA, with unprecedented sensitivity, resolution, and large fractional bandwidths, are enabling the exploration of the extragalactic RRL universe. We observed the radio quasar 3C 190 (z~1.2) with the LOFAR HBA. In reducing this data for spectroscopic analysis, we have placed special emphasis on bandpass calibration. We devised cross-correlation techniques to significantly identify the presence of RRLs in a low frequency spectrum. We demonstrate the utility of this method by applying it to existing low-frequency spectra of Cassiopeia A and M 82, and to the new observations of 3C 190. RRLs have been detected in the foreground of 3C 190 at z = 1.12355 (assuming a carbon origin), owing to the first detection of RRLs outside of the local universe (first reported in Emig et al. 2019). Towards the Galactic supernova remnant Cas A, we uncover three new detections: (1) C$ε$-transitions ($Δ$n = 5) for the first time at low radio frequencies, (2) H$α$-transitions at 64 MHz with a FWHM of 3.1 km/s, the most narrow and one of the lowest frequency detections of hydrogen to date, and (3) C$α$ at v$_{LSR}$ = 0 km/s in the frequency range 55-78 MHz for the first time. Additionally we recover C$α$, C$β$, C$γ$, and C$δ$ from the -47 km/s and -38 km/s components. In the nearby starburst galaxy, M 82, we do not find a significant feature. Our current searches for RRLs in LOFAR observations are limited to narrow (< 100 km/s) features, owing to the relatively small number of channels available for continuum estimation. Future strategies making use of larger contiguous frequency coverage would aid calibration to deeper sensitivities and broader features.

preprint2016arXiv

ALMA observations of cold molecular gas filaments trailing rising radio bubbles in PKS0745-191

We present ALMA observations of the CO(1-0) and CO(3-2) line emission tracing filaments of cold molecular gas in the central galaxy of the cluster PKS0745-191. The total molecular gas mass of 4.6 +/- 0.3 x 10^9 solar masses, assuming a Galactic X_{CO} factor, is divided roughly equally between three filaments each extending radially 3-5 kpc from the galaxy centre. The emission peak is located in the SE filament roughly 1 arcsec (2 kpc) from the nucleus. The velocities of the molecular clouds in the filaments are low, lying within +/-100 km/s of the galaxy's systemic velocity. Their FWHMs are less than 150 km/s, which is significantly below the stellar velocity dispersion. Although the molecular mass of each filament is comparable to a rich spiral galaxy, such low velocities show that the filaments are transient and the clouds would disperse on <10^7 yr timescales unless supported, likely by the indirect effect of magnetic fields. The velocity structure is inconsistent with a merger origin or gravitational free-fall of cooling gas in this massive central galaxy. If the molecular clouds originated in gas cooling even a few kpc from their current locations their velocities would exceed those observed. Instead, the projection of the N and SE filaments underneath X-ray cavities suggests they formed in the updraft behind bubbles buoyantly rising through the cluster atmosphere. Direct uplift of the dense gas by the radio bubbles appears to require an implausibly high coupling efficiency. The filaments are coincident with low temperature X-ray gas, bright optical line emission and dust lanes indicating that the molecular gas could have formed from lifted warmer gas that cooled in situ.

preprint2016arXiv

ALMA reveals optically thin, highly excited CO gas in the jet-driven winds of the galaxy IC5063

Using CO (4-3) and (2-1) Atacama Large Millimeter Array (ALMA) data, we prove that the molecular gas in the jet-driven winds of the galaxy IC5063 is more highly excited than the rest of the molecular gas in the disk of the same galaxy. On average, the CO(4-3)/CO(2-1) flux ratio is 1 for the disk and 5 for the jet accelerated or impacted gas. Spatially-resolved maps reveal that in regions associated with winds, the CO(4-3)/CO(2-1) flux ratio significantly exceeds the upper limit of 4 for optically thick gas. It frequently takes values between 5 and 11, and it occasionally further approaches the upper limit of 16 for optically thin gas. Excitation temperatures of 30-100 K are common for the molecules in these regions. If all of the outflowing molecular gas is optically thin, at 30-50 K, then its mass is 2*10^6 M_sun. This lower mass limit is an order of magnitude below the mass derived from the CO(2-1) flux in the case of optically thick gas. Molecular winds can thus be less massive, but more easily detectable at high z than they were previously thought to be.

preprint2016arXiv

Carbon and hydrogen radio recombination lines from the cold clouds towards Cassiopeia A

We use the Low Frequency Array to perform a systematic high spectral resolution investigation of the low-frequency 33-78 MHz spectrum along the line of sight to Cassiopeia A. We complement this with a 304-386 MHz Westerbork Synthesis Radio telescope observation. In this first paper we focus on the carbon radio recombination lines. We detect Cn$α$ lines at -47 and -38 km s$^{-1}$ in absorption for quantum numbers n=438-584 and in emission for n=257-278 with high signal to noise. These lines are associated with cold clouds in the Perseus spiral arm component. Hn$α$ lines are detected in emission for n=257-278. In addition, we also detect Cn$α$ lines at 0 km s$^{-1}$ associated with the Orion arm. We analyze the optical depth of these transitions and their line width. Our models show that the carbon line components in the Perseus arm are best fit with an electron temperature 85 K and an electron density 0.04 cm$^{-3}$ and can be constrained to within 15\%. The electron pressure is constrained to within 20\%. We argue that much of these carbon radio recombination lines arise in the CO-dark surface layers of molecular clouds where most of the carbon is ionized but hydrogen has made the transition from atomic to molecular. The hydrogen lines are clearly associated with the carbon line emitting clouds, but the low-frequency upperlimits indicate that they likely do not trace the same gas. Combining the hydrogen and carbon results we arrive at a firm lower limit to the cosmic ray ionization rate of 2.5$\times$10$^{-18}$ s$^{-1}$, but the actual value is likely much larger.

preprint2016arXiv

Embedded star formation in the extended narrow line region of Centaurus A: extreme mixing observed by MUSE

We present a detailed study of the complex ionization structure in a small (~250 pc) extended narrow line region (ENLR) cloud near Centaurus A using the Multi Unit Spectroscopic Explorer. This cloud is located in the so-called outer filament of ionized gas (about 15 kpc from the nucleus) where jet-induced star formation has been suggested to occur by different studies. We find that, despite the small size, a mixture of ionization mechanisms is operating, resulting in considerable complexity in the spatial ionization structure. The area includes two H II regions where star formation is occurring and another location where star formation must have ceased very recently. Interestingly, the extreme Balmer decrement of one of the star forming regions (H_alpha/H_beta~6) indicates that it is still heavily embedded in its natal cocoon of gas and dust. At all three locations a continuum counterpart is found with spectra matching those of O/B stars local to Centaurus A. The H II regions are embedded in a larger gas complex which is photoionized by the radiation of the central active galactic nucleus (AGN), but the O/B stars affect the spatial ionization pattern in the ENLR cloud very locally. In particular, in the surroundings of the youngest star forming region, we can isolate a tight mixing sequence in the diagnostic diagram going from gas with ionization due to a pure stellar continuum to gas only photoionized by the AGN. These results emphasize the complexity and the mixture of processes occurring in star forming regions under the influence of an AGN radiation. This is relevant for our understanding of AGN-induced star formation suggested to occur in a number of objects, including this region of Centaurus A. They also illustrate that these young stars influence the gas over only a limited region.

preprint2016arXiv

LOFAR VLBI Studies at 55 MHz of 4C 43.15, a z=2.4 Radio Galaxy

The correlation between radio spectral index and redshift has been exploited to discover high redshift radio galaxies, but its underlying cause is unclear. It is crucial to characterise the particle acceleration and loss mechanisms in high redshift radio galaxies to understand why their radio spectral indices are steeper than their local counterparts. Low frequency information on scales of $\sim$1 arcsec are necessary to determine the internal spectral index variation. In this paper we present the first spatially resolved studies at frequencies below 100 MHz of the $z = 2.4$ radio galaxy 4C 43.15 which was selected based on its ultra-steep spectral index ($α< -1$; $S_ν \sim ν^α$ ) between 365 MHz and 1.4 GHz. Using the International Low Frequency Array (LOFAR) Low Band Antenna we achieve sub-arcsecond imaging resolution at 55 MHz with VLBI techniques. Our study reveals low-frequency radio emission extended along the jet axis, which connects the two lobes. The integrated spectral index for frequencies $<$ 500 MHz is -0.83. The lobes have integrated spectral indices of -1.31$\pm$0.03 and -1.75$\pm$0.01 for frequencies $\geq$1.4 GHz, implying a break frequency between 500 MHz and 1.4 GHz. These spectral properties are similar to those of local radio galaxies. We conclude that the initially measured ultra-steep spectral index is due to a combination of the steepening spectrum at high frequencies with a break at intermediate frequencies.

preprint2016arXiv

Molecular Gas Along a Bright H-alpha Filament in 2A 0335+096 Revealed by ALMA

We present ALMA CO(1-0) and CO(3-2) observations of the brightest cluster galaxy (BCG) in the 2A 0335+096 galaxy cluster (z = 0.0346). The total molecular gas mass of (1.13+/-0.15) x 10^9 M_sun is divided into two components: a nuclear region and a 7 kpc long dusty filament. The central molecular gas component accounts for (3.2+/-0.4) x 10^8 M_sun of the total supply of cold gas. Instead of forming a rotationally-supported ring or disk, it is composed of two distinct, blueshifted clumps south of the nucleus and a series of low-significance redshifted clumps extending toward a nearby companion galaxy. The velocity of the redshifted clouds increases with radius to a value consistent with the companion galaxy, suggesting that an interaction between these galaxies <20 Myr ago disrupted a pre-existing molecular gas reservoir within the BCG. Most of the molecular gas, (7.8+/-0.9) x 10^8 M_sun, is located in the filament. The CO emission is co-spatial with a 10^4 K emission-line nebula and soft X-rays from 0.5 keV gas, indicating that the molecular gas has cooled out of the intracluster medium over a period of 25-100 Myr. The filament trails an X-ray cavity, suggesting that the gas has cooled from low entropy gas that has been lifted out of the cluster core and become thermally unstable. We are unable to distinguish between inflow and outflow along the filament with the present data. Cloud velocities along the filament are consistent with gravitational free-fall near the plane of the sky, although their increasing blueshifts with radius are consistent with outflow.

preprint2016arXiv

The warm molecular hydrogen of PKS B1718-649: feeding a newly born radio AGN

We present new SINFONI VLT observations of molecular hydrogen (H2) in the central regions (< 2.5 kpc) of the youngest and closest radio source PKS B1718-649. We study the distribution of the H2 traced by the 1-0 S(1) ro-vibrational line, revealing a double disk structure with the kinematics of both disks characterised by rotation. An outer disk (r > 650 pc) is aligned with other components of the galaxy (atomic hydrogen, stars, dust), while the inner disk (r< 600 pc) is perpendicular to it and is polar with respect to the stellar distribution. However, in the innermost 75 pc, the data show the presence of H2 gas redshifted with respect to the rotating inner disk ($Δv\,$+150 km/s) which may trace gas falling into the super-massive black hole associated with the central radio source. Along the same line of sight, earlier observations had shown the presence in the central regions of PKS B1718-649 of clouds of atomic hydrogen with similar unsettled kinematics. The range of velocities and mass of these unsettled clouds of HI and H2 suggest they may be actively contributing in fuelling the central newly-born radio source.

preprint2015arXiv

Cold gas and the disruptive effect of a young radio jet

Newly born and young radio sources are in a delicate phase of their life. Their jets are fighting their way through the surrounding gaseous medium, strongly experiencing this interaction while, at the same time, impacting and affecting the interstellar medium (ISM). Here we present the results from two studies of HI (in absorption) and molecular gas illustrating what can be learned from these phases of the gas. We first describe a statistical study with the WSRT. The study shows that the young radio sources not only have an higher detection rate of HI, but also systematically broader and more asymmetric HI profiles, most of them blueshifted. This supports the idea that we are looking at young radio jets making their way through the surrounding ISM, which also appears to be, on average, richer in gas than in evolved radio sources. Signatures of the impact of the jet are seen in the kinematics of the gas. However, even among the young sources, we identify a population that remains undetected in HI even after stacking their profiles. Orientation effects can only partly explain the result. These objects either are genuinely gas-poor or have different conditions of the medium, e.g. higher spin temperature. We further present the ALMA study of molecular gas in IC5063 to trace in detail the jet impacting the ISM. The kinematics of the cold, molecular gas co-spatial with the radio plasma shows this process in action. The ALMA data reveal a fast outflow of molecular gas extending along the entire radio jet (~1 kpc), with the highest outflow velocities at the location of the brighter hot-spot. We propose a scenario where the radio jet is expanding into a clumpy medium, interacting directly with the clouds and inflating a cocoon that drives a lateral outflow into the ISM.

preprint2015arXiv

PKS B1718-649: an HI and H2 perspective on the birth of a compact radio source

We present neutral hydrogen (HI) and warm molecular hydrogen (H2) observations of the young (10^2 years) radio galaxy PKS B1718-649. We study the morphology and the kinematics of both gas components, focusing, in particular, on their properties in relation to the triggering of the radio activity. The regular kinematics of the large scale HI disk, seen in emission, suggests that an interaction event occurred too long ago to be responsible for the recent triggering of the radio activity. In absorption, we detect two absorption lines along the narrow line of sight of the compact (r<2 pc) radio source. The lines trace two clouds with opposite radial motions. These may represent a population of clouds in the very inner regions of the galaxy, which may be involved in triggering the radio activity. The warm molecular hydrogen (H2 1-0 S(1) ro-vibrational line) in the innermost kilo-parsec of the galaxy appears to be distributed in a circum-nuclear disk following the regular kinematics of the HI and of the stellar component. An exception to this behaviour arises only in the very centre, where a highly dispersed component is detected. These particular HI and H2 features suggest that a strong interplay between the radio source and the surrounding ISM is on-going. The physical properties of the cold gas in the proximity of the radio source may regulate the accretion recently triggered in this AGN.

preprint2015arXiv

Radio Jets Clearing theWay Through a Galaxy: Watching Feedback in Action in the Seyfert galaxy IC 5063

High-resolution (0.5 arcsec) CO(2-1) observations performed with the Atacama Large Millimetre/submillimetre Array have been used to trace the kinematics of the molecular gas in the Seyfert 2 galaxy{IC~5063}. Although one of the most radio-loud Seyfert galaxy, IC~5063 is a relatively weak radio source (P_1.4GHz = 3 x 10^23 W Hz^-1). The data reveal that the kinematics of the gas is very complex. A fast outflow of molecular gas extends along the entire radio jet (~ 1 kpc), with the highest outflow velocities about 0.5 kpc from the nucleus, at the location of the brighter hot-spot in the W lobe. All the observed characteristics can be described by a scenario of a radio plasma jet expanding into a clumpy medium, interacting directly with the clouds and inflating a cocoon that drives a lateral outflow into the interstellar medium. This suggests that most of the observed cold molecular outflow is due to fast cooling of the gas after the passage of a shock and that it is the end product of the cooling process.

preprint2015arXiv

The fast molecular outflow in the Seyfert galaxy IC5063 as seen by ALMA

We use high-resolution (0.5 arcsec) CO(2-1) observations performed with ALMA to trace the kinematics of the molecular gas in the Seyfert 2 galaxy IC5063. A fast outflow of molecular gas extends along the entire radio jet, with the highest outflow velocities about 0.5kpc from the nucleus, at the location of the brighter hot-spot in the W lobe. The data show that a massive, fast outflow with velocities up to 650 km/s of cold molecular gas is present, in addition to one detected earlier in warm H2, HI and ionised gas. Both the central AGN and the radio jet could energetically drive the outflow. However, the characteristics of the outflowing gas point to the radio jet being the main driver. This is important, because IC5063, although one of the most powerful Seyfert galaxies, is a relatively weak radio source (P = 3x10^23 W/Hz). All the observed characteristics can be described by a scenario of a radio plasma jet expanding into a clumpy medium, interacting directly with the clouds and inflating a cocoon that drives a lateral outflow into the interstellar medium. This model is consistent with results obtained by recent simulations such as those of Wagner et al.. A stronger, direct interaction between the jet and a gas cloud is present at the location of the brighter W lobe. Even assuming the most conservative values for the conversion factor CO-to-H2, the mass of the outflowing gas is between 1.9 and 4.8x10^7 Msun. These amounts are much larger than those of the outflow of warm gas (molecular and ionized) and somewhat larger than of the HI outflow. This suggests that most of the observed cold molecular outflow is due to fast cooling after being shocked. This gas is the end product of the cooling process. Our CO observations demonstrate that fast outflows of molecular gas can be driven by relativistic jets.

preprint2015arXiv

The jet-ISM interaction in the Outer Filament of Centaurus A

The interaction between the radio plasma ejected by the active nucleus of a galaxy and the surrounding medium is a key process that can have a strong impact on the interstellar medium of the galaxy and hence on galaxy evolution. The closest laboratory where we can observe and investigate this phenomenon is the radio galaxy Centaurus A. About 15 kpc north-east of this galaxy, a particularly complex region is found: the so-called Outer Filament where jet-cloud interactions have been proposed to occur. We investigate the presence of signatures of jet-ISM interaction by a detailed study of the kinematics of the ionized gas, expanding on previous results obtained from the HI. We observed two regions of the outer filament with VLT/VIMOS in the IFU observing mode. Emission from Hbeta and [OIII]4959,5007Å is detected in both pointings. We found two distinct kinematical components of ionized gas that well match the kinematics of the nearby HI cloud. One component follows the regular kinematics of the rotating gas while the second shows similar velocities to those of the nearby HI component thought to be disturbed by an interaction with the radio jet. We suggest that the ionized and atomic gas are part of the same dynamical gas structure originating as result of the merger that shaped Centaurus A and which is regularly rotating around Centaurus A as proposed by other authors. The gas (ionized and HI) with anomalous velocities is tracing the interaction of the Large-Scale radio Jet with the ISM, suggesting that, although poorly collimated as structure, the jet is still active. However, we can exclude that a strong shock is driving the ionization of the gas. It is likely that a combination of jet entrainment and photoionization by the UV continuum from the central engine is needed in order to explain both the ionization and the kinematics of the gas in the Outer Filament.

preprint2015arXiv

The Physics of the Cold Neutral Medium: Low-frequency Carbon Radio Recombination Lines with the Square Kilometre Array

The Square Kilometre Array (SKA) will transform our understanding of the role of the cold, atomic gas in galaxy evolution. The interstellar medium (ISM) is the repository of stellar ejecta and the birthsite of new stars and, hence, a key factor in the evolution of galaxies over cosmic time. Cold, diffuse, atomic clouds are a key component of the ISM, but so far this phase has been difficult to study, because its main tracer, the HI 21 cm line, does not constrain the basic physical information of the gas (e.g., temperature, density) well. The SKA opens up the opportunity to study this component of the ISM through a complementary tracer in the form of low-frequency (<350 MHz) carbon radio recombination lines (CRRL). These CRRLs provide a sensitive probe of the physical conditions in cold, diffuse clouds. The superb sensitivity, large field of view, frequency resolution and coverage of the SKA allows for efficient surveys of the sky, that will revolutionize the field of low-frequency recombination line studies. By observing these lines with the SKA we will be able determine the thermal balance, chemical enrichment, and ionization rate of the cold, atomic medium from degree-scales down to scales corresponding to individual clouds and filaments in our Galaxy, the Magellanic Clouds and beyond. Furthermore, being sensitive only to the cold, atomic gas, observations of low-frequency CRRLs with the SKA will aid in disentangling the warm and cold constituents of the HI 21 cm emission.

preprint2014arXiv

A $10^{10}$ Solar Mass Flow of Molecular Gas in the Abell 1835 Brightest Cluster Galaxy

We report ALMA Early Science observations of the Abell 1835 brightest cluster galaxy (BCG) in the CO (3-2) and CO (1-0) emission lines. We detect $5\times 10^{10}~\rm M_\odot$ of molecular gas within 10 kpc of the BCG. Its ensemble velocity profile width of $\sim 130 ~\rm km~s^{-1}$ FWHM is too narrow for the molecular cloud sto be supported in the galaxy by dynamic pressure. The gas may instead be supported in a rotating, turbulent disk oriented nearly face-on. Roughly $10^{10}~\rm M_\odot$ of molecular gas is projected $3-10 ~\rm kpc$ to the north-west and to the east of the nucleus with line of sight velocities lying between $-250 ~\rm km~s^{-1}$ to $+480 ~\rm km~s^{-1}$ with respect to the systemic velocity. The high velocity gas may be either inflowing or outflowing. However, the absence of high velocity gas toward the nucleus that would be expected in a steady inflow, and its bipolar distribution on either side of the nucleus, are more naturally explained as outflow. Star formation and radiation from the AGN are both incapable of driving an outflow of this magnitude. If so, the molecular outflow may be associated a hot outflow on larger scales reported by Kirkpatrick and colleagues. The molecular gas flow rate of approximately $200~\rm M_\odot ~yr^{-1}$ is comparable to the star formation rate of $100-180~\rm M_\odot ~yr^{-1}$ in the central disk. How radio bubbles would lift dense molecular gas in their updrafts, how much gas will be lost to the BCG, and how much will return to fuel future star formation and AGN activity are poorly understood. Our results imply that radio-mechanical (radio mode) feedback not only heats hot atmospheres surrounding elliptical galaxies and BCGs, it is able to sweep higher density molecular gas away from their centers.

preprint2014arXiv

Discovery of Carbon Radio Recombination Lines in absorption towards Cygnus~A

We present the first detection of carbon radio recombination line absorption along the line of sight to Cygnus A. The observations were carried out with the LOw Frequency ARray in the 33 to 57 MHz range. These low frequency radio observations provide us with a new line of sight to study the diffuse, neutral gas in our Galaxy. To our knowledge this is the first time that foreground Milky Way recombination line absorption has been observed against a bright extragalactic background source. By stacking 48 carbon $α$ lines in the observed frequency range we detect carbon absorption with a signal-to-noise ratio of about 5. The average carbon absorption has a peak optical depth of 2$\times$10$^{-4}$, a line width of 10 km s$^{-1}$ and a velocity of +4 km s$^{-1}$ with respect to the local standard of rest. The associated gas is found to have an electron temperature $T_{e}\sim$ 110 K and density $n_{e}\sim$ 0.06 cm$^{-3}$. These properties imply that the observed carbon $α$ absorption likely arises in the cold neutral medium of the Orion arm of the Milky Way. Hydrogen and helium lines were not detected to a 3$σ$ peak optical depth limit of 1.5$\times$10$^{-4}$ for a 4 km s$^{-1}$ channel width. Radio recombination lines associated with Cygnus A itself were also searched for, but are not detected. We set a 3$σ$ upper limit of 1.5$\times$10$^{-4}$ for the peak optical depth of these lines for a 4 km s$^{-1}$ channel width.

preprint2014arXiv

Discovery of Carbon Radio Recombination Lines in M82

Carbon radio recombination lines (RRLs) at low frequencies (<=500 MHz) trace the cold, diffuse phase of the interstellar medium, which is otherwise difficult to observe. We present the detection of carbon RRLs in absorption in M82 with LOFAR in the frequency range of 48-64 MHz. This is the first extragalactic detection of RRLs from a species other than hydrogen, and below 1 GHz. Since the carbon RRLs are not detected individually, we cross-correlated the observed spectrum with a template spectrum of carbon RRLs to determine a radial velocity of 219 +- 9 km/s . Using this radial velocity, we stack 22 carbon-α transitions from quantum levels n = 468-508 to achieve an 8.5 sigma detection. The absorption line profile exhibits a narrow feature with peak optical depth of 0.003 and FWHM of 31 km/s. Closer inspection suggests that the narrow feature is superimposed on a broad, shallow component. The total line profile appears to be correlated with the 21 cm H I line profile reconstructed from H I absorption in the direction of supernova remnants in the nucleus. The narrow width and centroid velocity of the feature suggests that it is associated with the nuclear starburst region. It is therefore likely that the carbon RRLs are associated with cold atomic gas in the direction of the nucleus of M82.

preprint2014arXiv

Exact Bound-Bound Gaunt Factor Values for Quantum Levels up to n=2000

Comparison of observations of radio recombination lines in the interstellar medium with theoretical models can be used to constrain electron temperature and density of the gas. An important component of the models is spontaneous transition rates between bound levels. Calculating these rates relies on accurate bound-bound oscillator strengths, which can be cast in terms of the Gaunt factor. The Gaunt factor contains terminating hypergeometric functions that cannot be calculated with sufficient accuracy for high quantum levels ($n \gtrsim 50$) by standard machine-precision methods. Methods to overcome the accuracy problem have been developed, which include asymptotic expansions and recursion relations. These methods, used in astrophysical models to calculate oscillator strengths, can introduce errors, sometimes up to as much as $\sim 8$ per cent. Detections of radio recombination lines with the new Low Frequency Array (LOFAR) has prompted an examination of theoretical models of the interstellar medium. We revisit the calculation of the Gaunt factor, employing modern arbitrary-precision computational methods to tabulate the Gaunt factor for transitions up to quantum level $n=2000$, sufficient to model low frequency Carbon radio recombination lines. The calculations provide a relative error of $\sim3\times 10^{-4}$ when compared to more detailed calculations including relativistic corrections. Our values for the Gaunt factor are provided for download in a tabular format to be used for a wide range of applications.

preprint2014arXiv

Extended, Dusty Star Formation Fueled by a Residual Cooling Flow in the Cluster of Galaxies Sersic 159-03

While the cooling of the hot intra-cluster medium (ICM) in the cores of galaxy clusters is mostly counteracted by heating from the central active galactic nucleus (AGN), the balance is not perfect. This can lead to residual cooling flows and low-level star formation, the physics of which is not well understood. Here we present a detailed study of the residual cooling flow in the center of the low mass galaxy cluster Sersic 159-03 (A S1101) using far-ultraviolet imaging from the Hubble Space Telescope and far-infrared (FIR) spectroscopy and photometry from the Herschel space observatory, along with a wealth of archival data. We detect extended emission at UV, FIR, and [CII], indicating a star formation rate of ~1-3 Msun/yr, depending on the indicator and assumptions made. The most recently formed stars appear spatially coincident with the lowest entropy ICM. We speculate that this low-entropy gas has been displaced by the central AGN ~7.5 kpc north of the cD galaxy. These data demonstrate that the displacement of the cooling core from the direct vicinity of the central AGN can temporarily break the feedback cycle and lead to cooling and star formation that is offset from the center of the galaxy. We find an abundance (~10^7 Msun) of cold (20K) dust in the center of the cluster and a second FIR peak ~30kpc to the north of the central galaxy. If confirmed to be associated with the cooling filaments, this would be the most extended complex of dust yet found in a cool core cluster.

preprint2014arXiv

Jet acceleration of the fast molecular outflows in the Seyfert galaxy IC5063

Massive outflows driven by active galactic nuclei (AGN) are widely recognised to play a key role in the evolution of galaxies, heating the ambient gas, expelling it from the nuclear regions, and thereby affecting the star formation histories of the galaxy bulges. It has been proposed that the powerful jets of relativistic particles launched by some AGN can both accelerate and heat the molecular gas, which often dominates the mass budgets of the outflows. However, clear evidence for this mechanism in the form of detailed associations between the molecular gas kinematics and features in the radio-emitting jets has been lacking. Here we show that the warm molecular hydrogen gas in the western radio lobe of the Seyfert galaxy IC5063 is moving at high velocities - up to 600 km/s - relative to the galaxy disk. This suggests that the molecules have been accelerated by fast shocks driven into the interstellar medium (ISM) by the expanding radio jets. These results demonstrate the general feasibility of accelerating molecular outflows in fast shocks driven by AGN.

preprint2014arXiv

LOFAR low-band antenna observations of the 3C295 and Bootes fields: source counts and ultra-steep spectrum sources

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($α< -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Boötes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

preprint2014arXiv

The origin of cold gas in giant elliptical galaxies and its role in fueling radio-mode AGN feedback

We present a multi-wavelength study of the interstellar medium in eight nearby, X-ray and optically bright, giant elliptical galaxies. Using Herschel PACS, we map the cold gas in the lines of [CII], [OI], and [OIb]. Additionally, we present Ha+[NII] imaging of warm ionized gas with the SOAR telescope, and a study of the hot X-ray emitting plasma with Chandra. All systems with extended Ha emission in our sample (6/8 galaxies) display significant [CII] line emission indicating the presence of cold gas. This emission is co-spatial with the Ha+[NII] emitting nebulae and the lowest entropy X-ray emitting plasma. The entropy profiles of the hot galactic atmospheres show a clear dichotomy, with the systems displaying extended emission line nebulae having lower entropies beyond r~1 kpc than the cold-gas-poor systems. We show that while the hot atmospheres of the cold-gas-poor galaxies are thermally stable outside of their innermost cores, the atmospheres of the cold-gas-rich systems are prone to cooling instabilities. This result indicates that the cold gas is produced chiefly by thermally unstable cooling from the hot phase. We show that cooling instabilities may develop more easily in rotating systems and discuss an alternative condition for thermal instability for this case. The hot atmospheres of cold-gas-rich galaxies display disturbed morphologies indicating that the accretion of clumpy multiphase gas in these systems may result in variable power output of the AGN jets, potentially triggering sporadic, larger outbursts. In the two cold-gas-poor, X-ray morphologically relaxed galaxies of our sample, NGC 1399 and NGC 4472, powerful AGN outbursts may have destroyed or removed most of the cold gas from the cores, allowing the jets to propagate and deposit most of their energy further out, increasing the entropy of the hot galactic atmospheres and leaving their cores relatively undisturbed.

preprint2013arXiv

A Ten Billion Solar Mass Outflow of Molecular Gas Launched by Radio Bubbles in the Abell 1835 Brightest Cluster Galaxy

We report ALMA Early Science observations of the Abell 1835 brightest cluster galaxy (BCG) in the CO (3-2) and CO (1-0) emission lines. We detect 5E10 solar masses of molecular gas within 10 kpc of the BCG. Its velocity width of ~130 km/s FWHM is too narrow to be supported by dynamical pressure. The gas may instead be supported in a rotating, turbulent disk oriented nearly face-on. The disk is forming stars at a rate of 100-180 solar masses per year. Roughly 1E10 solar masses of molecular gas is projected 3-10 kpc to the north-west and to the east of the nucleus with line of sight velocities lying between -250 km/s to +480 km/s with respect to the systemic velocity. Although inflow cannot be ruled out, the rising velocity gradient with radius is consistent with a broad, bipolar outflow driven by radio jets or buoyantly rising X-ray cavities. The molecular outflow may be associated with an outflow of hot gas in Abell 1835 seen on larger scales. Molecular gas is flowing out of the BCG at a rate of approximately 200 solar masses per year, which is comparable to its star formation rate. How radio bubbles lift dense molecular gas in their updrafts, how much gas will be lost to the BCG, and how much will return to fuel future star formation and AGN activity are poorly understood. Our results imply that radio-mechanical (radio mode) feedback not only heats hot atmospheres surrounding elliptical galaxies and BCGs, it is able to sweep higher density molecular gas away from their centers.

preprint2013arXiv

Cold gas dynamics in Hydra-A: evidence for a rotating disk

We present multi-frequency observations of the radio galaxy Hydra-A (3C218) located in the core of a massive, X-ray luminous galaxy cluster. IFU spectroscopy is used to trace the kinematics of the ionised and warm molecular hydrogen which are consistent with a ~ 5 kpc rotating disc. Broad, double-peaked lines of CO(2-1), [CII]157 $μ$m and [OI]63 $μ$m are detected. We estimate the mass of the cold gas within the disc to be M$_{gas}$ = 2.3 $\pm$ 0.3 x 10$^9$ M$_{\odot}$. These observations demonstrate that the complex line profiles found in the cold atomic and molecular gas are related to the rotating disc or ring of gas. Finally, an HST image of the galaxy shows that this gas disc contains a substantial mass of dust. The large gas mass, SFR and kinematics are consistent with the levels of gas cooling from the ICM. We conclude that the cold gas originates from the continual quiescent accumulation of cooled ICM gas. The rotation is in a plane perpendicular to the projected orientation of the radio jets and ICM cavities hinting at a possible connection between the kpc-scale cooling gas and the accretion of material onto the black hole. We discuss the implications of these observations for models of cold accretion, AGN feedback and cooling flows.

preprint2013arXiv

LOFAR detections of low-frequency radio recombination lines towards Cassiopeia A

Cassiopeia A was observed using the Low-Band Antennas of the LOw Frequency ARray (LOFAR) with high spectral resolution. This allowed a search for radio recombination lines (RRLs) along the line-of-sight to this source. Five carbon-alpha RRLs were detected in absorption between 40 and 50 MHz with a signal-to-noise ratio of > 5 from two independent LOFAR datasets. The derived line velocities (v_LSR ~ -50 km/s) and integrated optical depths (~ 13 s^-1) of the RRLs in our spectra, extracted over the whole supernova remnant, are consistent within each LOFAR dataset and with those previously reported. For the first time, we are able to extract spectra against the brightest hotspot of the remnant at frequencies below 330 MHz. These spectra show significantly higher (15-80 %) integrated optical depths, indicating that there is small-scale angular structure on the order of ~1 pc in the absorbing gas distribution over the face of the remnant. We also place an upper limit of 3 x 10^-4 on the peak optical depths of hydrogen and helium RRLs. These results demonstrate that LOFAR has the desired spectral stability and sensitivity to study faint recombination lines in the decameter band.

preprint2013arXiv

Massive molecular gas flows in the Abell 1664 brightest cluster galaxy

We report ALMA Early Science CO(1-0) and CO(3-2) observations of the brightest cluster galaxy (BCG) in Abell 1664. The BCG contains 1.1x10^{10} solar masses of molecular gas divided roughly equally between two distinct velocity systems: one from -250 to +250 km/s centred on the BCG's systemic velocity and a high velocity system blueshifted by 570 km/s with respect to the systemic velocity. The BCG's systemic component shows a smooth velocity gradient across the BCG center with velocity proportional to radius suggestive of solid body rotation about the nucleus. However, the mass and velocity structure are highly asymmetric and there is little star formation coincident with a putative disk. It may be an inflow of gas that will settle into a disk over several 10^8 yr. The high velocity system consists of two gas clumps, each ~2 kpc across, located to the north and southeast of the nucleus. Each has a line of sight velocity spread of 250-300 km/s. The velocity of the gas in the high velocity system tends to increase towards the BCG center and could signify a massive high velocity flow onto the nucleus. However, the velocity gradient is not smooth and these structures are also coincident with low optical-UV surface brightness regions, which could indicate dust extinction associated with each clump. If so, the high velocity gas would be projected in front of the BCG and moving toward us along the line of sight in a massive outflow most likely driven by the AGN. A merger origin is unlikely but cannot be ruled out.

preprint2013arXiv

The nature of filamentary cold gas in the core of the Virgo Cluster

We present a multi-wavelength study of the emission-line nebulae located southeast of the nucleus of M87, the central dominant galaxy of the Virgo Cluster. We report the detection of far-infrared (FIR) [CII] line emission from the nebulae using observations made with Herschel PACS. The infrared line emission is extended and cospatial with optical Hα+[NII], far-ultraviolet CIV lines, and soft X-ray emission. The filamentary nebulae evidently contain multi-phase material spanning a temperature range of at least 5 orders of magnitude, from ~100 K to ~10^7 K. This material has most likely been uplifted by the AGN from the center of M87. The thermal pressure of the 10^4 K phase appears to be significantly lower than that of the surrounding hot intra-cluster medium (ICM) indicating the presence of additional turbulent and magnetic pressure in the filaments. If the turbulence in the filaments is subsonic then the magnetic field strength required to balance the pressure of the surrounding ICM is B~30-70 μG. The spectral properties of the soft X-ray emission from the filaments indicate that it is due to thermal plasma with kT~0.5-1 keV, which is cooling by mixing with the cold gas and/or radiatively. Charge exchange can be ruled out as a significant source of soft X-rays. Both cooling and mixing scenarios predict gas with a range of temperatures. This is at first glance inconsistent with the apparent lack of X-ray emitting gas with kT<0.5 keV. However, we show that the missing very soft X-ray emission could be absorbed by the cold gas in the filaments with an integrated hydrogen column density of ~1.6x10^21 cm^-2, providing a natural explanation for the apparent temperature floor to the X-ray emission at kT~0.5 keV. The FIR through ultra-violet line emission is most likely primarily powered by the ICM particles penetrating the cold gas following a shearing induced mixing process.

preprint2012arXiv

A Giant Metrewave Radio Telescope/Chandra view of IRAS 09104+4109: A type 2 QSO in a cooling flow

IRAS 09104+4109 is a rare example of a dust enshrouded type 2 QSO in the centre of a cool-core galaxy cluster. Previous observations of this z=0.44 system showed that as well as powering the hyper-luminous infrared emission of the cluster-central galaxy, the QSO is associated with a double-lobed radio source. However, the steep radio spectral index and misalignment between the jets and ionised optical emission suggested that the orientation of the QSO had recently changed. We use a combination of new, multi-band Giant Metrewave Radio Telescope observations and archival radio data to confirm that the jets are no longer powered by the QSO, and estimate their age to be 120-160 Myr. This is in agreement with the ~70-200 Myr age previously estimated for star-formation in the galaxy. Previously unpublished Very Long Baseline Array data reveal a 200 pc scale double radio source in the galaxy core which is more closely aligned with the current QSO axis and may represent a more recent period of jet activity. These results suggest that the realignment of the QSO, the cessation of jet activity, and the onset of rapid star-formation may have been caused by a gas-rich galaxy merger. A Chandra X-ray observation confirms the presence of cavities associated with the radio jets, and we estimate the energy required to inflate them to be ~7.7x10^60 erg. The mechanical power of the jets is sufficient to balance radiative cooling in the cluster, provided they are efficiently coupled to the intra-cluster medium (ICM). We find no evidence of direct radiative heating and conclude that the QSO either lacks the radiative luminosity to heat the ICM, or that it requires longer than 100-200 Myr to significantly impact its environment. [Abridged]

preprint2012arXiv

Multiphase Signatures of AGN Feedback in Abell 2597

We present new Chandra X-ray observations of the brightest cluster galaxy (BCG) in the cool core cluster Abell 2597. The data reveal an extensive kpc-scale X-ray cavity network as well as a 15 kpc filament of soft-excess gas exhibiting strong spatial correlation with archival VLA radio data. In addition to several possible scenarios, multiwavelength evidence may suggest that the filament is associated with multiphase (10^3 - 10^7 K) gas that has been entrained and dredged-up by the propagating radio source. Stemming from a full spectral analysis, we also present profiles and 2D spectral maps of modeled X-ray temperature, entropy, pressure, and metal abundance. The maps reveal an arc of hot gas which in projection borders the inner edge of a large X-ray cavity. Although limited by strong caveats, we suggest that the hot arc may be (a) due to a compressed rim of cold gas pushed outward by the radio bubble or (b) morphologically and energetically consistent with cavity-driven active galactic nucleus (AGN) heating models invoked to quench cooling flows, in which the enthalpy of a buoyant X-ray cavity is locally thermalized as ambient gas rushes to refill its wake. If confirmed, this would be the first observational evidence for this model.

preprint2012arXiv

Residual Cooling and Persistent Star Formation amid AGN Feedback in Abell 2597

New Chandra X-ray and Herschel FIR observations enable a multiwavelength study of active galactic nucleus (AGN) heating and intracluster medium (ICM) cooling in the brightest cluster galaxy of Abell 2597. The new Chandra observations reveal the central < 30 kiloparsec X-ray cavity network to be more extensive than previously thought, and associated with enough enthalpy to theoretically inhibit the inferred classical cooling flow. Nevertheless, we present new evidence, consistent with previous results, that a moderately strong residual cooling flow is persisting at 4%-8% of the classically predicted rates in a spatially structured manner amid the feedback-driven excavation of the X-ray cavity network. New Herschel observations are used to estimate warm and cold dust masses, a lower-limit gas-to-dust ratio, and a star formation rate consistent with previous measurements. The cooling time profile of the ambient X-ray atmosphere is used to map the locations of the observational star formation entropy threshold as well as the theoretical thermal instability threshold. Both lie just outside the < 30 kpc central region permeated by X-ray cavities, and star formation as well as ionized and molecular gas lie interior to both. The young stars are distributed in an elongated region that is aligned with the radio lobes, and their estimated ages are both younger and older than the X-ray cavity network, suggesting both jet-triggered as well as persistent star formation over the current AGN feedback episode. Bright X-ray knots that are coincident with extended Ly-alpha and FUV continuum filaments motivate a discussion of structured cooling from the ambient hot atmosphere along a projected axis that is perpendicular to X-ray cavity and radio axis. We conclude that the cooling ICM is the dominant contributor of the cold gas reservoir fueling star formation and AGN activity in the Abell 2597 BCG.

preprint2012arXiv

The "toothbrush-relic": evidence for a coherent linear 2-Mpc scale shock wave in a massive merging galaxy cluster?

Some merging galaxy clusters host diffuse extended radio emission, so-called radio halos and relics. Here we present observations between 147 MHz and 4.9 GHz of a new radio-selected galaxy cluster 1RXS J0603.3+4214 (z=0.225). The cluster is also detected as an extended X-ray source in the RASS. It hosts a large bright 1.9 Mpc radio relic, an elongated ~2 Mpc radio halo, and two smaller radio relics. The large radio relic has a peculiar linear morphology. For this relic we observe a clear spectral index gradient, in the direction towards the cluster center. We performed Rotation Measure (RM) Synthesis between 1.2 and 1.7 GHz. The results suggest that for the west part of the large relic some of the Faraday rotation is caused by ICM and is not only due to galactic foregrounds. We also carried out a detailed spectral analysis of this radio relic and created radio color-color diagrams. We find (i) an injection spectral index of -0.6 to -0.7, (ii) steepening spectral index and increasing spectral curvature in the post-shock region, and (iii) an overall power-law spectrum between 74 MHz and 4.9 GHz with α=-1.10 \pm 0.02. Mixing of emission in the beam from regions with different spectral ages is probably the dominant factor that determines the shape of the radio spectra. Changes in the magnetic field, total electron content, or adiabatic gains/losses do not play a major role. A model in which particles are (re)accelerated in a first order Fermi process at the front of the relic provides the best match to the observed spectra. We speculate that in the post-shock region particles are re-accelerated by merger induced turbulence to form the radio halo as the relic and halo are connected. The 1RXS J0603.3+4214 merger is probably more complex than the "simple'" binary merger events that are thought to give rise to symmetric double radio relics.

preprint2011arXiv

Far Ultraviolet Emission in the A2597 and A2204 Brightest Cluster Galaxies

We use the Hubble Space Telescope ACS/SBC and Very Large Telescope FORS cameras to observe the Brightest Cluster Galaxies in Abell 2597 and Abell 2204 in the far-ultraviolet (FUV) F150LP and optical U, B, V, R, I Bessel filters. The FUV and U band emission is enhanced in bright, filamentary structures surrounding the BCG nuclei. These filaments can be traced out to 20 kpc from the nuclei in the FUV. Excess FUV and U band light is determined by removing emission due to the underlying old stellar population and mapped with 1 arcsec spatial resolution over the central 20 kpc regions of both galaxies. We find the FUV and U excess emission to be spatially coincident and a stellar interpretation requires the existence of a significant amount of 10000-50000 K stars. Correcting for nebular continuum emission and dust intrinsic to the BCG further increases the FUV to U band emission ratio and implies that stars alone may not suffice to explain the observations. However, lack of detailed information on the gas and dust distribution and extinction law in these systems prevents us from ruling out a purely stellar origin. Non-stellar processes, such as the central AGN, Scattering, Synchrotron and Bremsstrahlung emission are investigated and found to not be able to explain the FUV and U band measurements in A2597. Contributions from non-thermal processes not treated here should be investigated. Comparing the FUV emission to the optical H-alpha line emitting nebula shows good agreement on kpc-scales in both A2597 and A2204. In concordance with an earlier investigation by O'Dea et al. (2004) we find that O-stars can account for the ionising photons necessary to explain the observed H-alpha line emission.

preprint2010arXiv

Herschel observations of FIR emission lines in brightest cluster galaxies

The question of how much gas cools in the cores of clusters of galaxies has been the focus of many, multiwavelength studies in the past 30 years. In this letter we present the first detections of the strongest atomic cooling lines, [C II], [O I] and [N I] in two strong cooling flow clusters, A1068 and A2597, using Herschel PACS. These spectra indicate that the substantial mass of cold molecular gas (>10^9 Mo) known to be present in these systems is being irradiated by intense UV radiation, most probably from young stars. The line widths of these FIR lines indicate that they share dynamics similar but not identical to other ionised and molecular gas traced by optical, near-infrared and CO lines. The relative brightness of the FIR lines compared to CO and FIR luminosity is consistent with other star-forming galaxies indicating that the properties of the molecular gas clouds in cluster cores and the stars they form are not unusual. These results provide additional evidence for a reservoir of cold gas that is fed by the cooling of gas in the cores of the most compact clusters and provide important diagnostics of the temperature and density of the dense clouds this gas resides in.

preprint2010arXiv

Herschel photometry of brightest cluster galaxies in cooling flow clusters

The dust destruction timescales in the cores of clusters of galaxies are relatively short given their high central gas densities. However, substantial mid-infrared and sub-mm emission has been detected in many brightest cluster galaxies. In this letter we present Herschel PACS and SPIRE photometry of the brightest cluster galaxy in three strong cooling flow clusters, A1068, A2597 and Zw3146. This photometry indicates that a substantial mass of cold dust is present (>3 x 10^7 Mo) at temperatures significantly lower (20-28K) than previously thought based on limited MIR and/or sub-mm results. The mass and temperature of the dust appear to match those of the cold gas traced by CO with a gas-to-dust ratio of 80-120.

preprint2010arXiv

Sub-millimeter to centimeter excess emission from the Magellanic Clouds. I. Global spectral energy distribution

In order to reconstruct the global SEDs of the Magellanic Clouds over eight decades in spectral range, we combined literature flux densities representing the entire LMC and SMC respectively, and complemented these with maps extracted from the WMAP and COBE databases covering the missing the 23--90 GHz (13--3.2 mm) and the poorly sampled 1.25--250 THz (240--1.25 micron). We have discovered a pronounced excess of emission from both Magellanic Clouds, but especially the SMC, at millimeter and sub-millimeter wavelengths. We also determined accurate thermal radio fluxes and very low global extinctions for both LMC and SMC. Possible explanations are briefly considered but as long as the nature of the excess emission is unknown, the total dust masses and gas-to-dust ratios of the Magellanic Clouds cannot reliably be determined.

preprint2010arXiv

The Distribution and Condition of the Warm Molecular Gas in Abell 2597 and Sersic 159-03

We have used the SINFONI integral field spectrograph to map the near-infrared K-band emission lines of molecular and ionised hydrogen in the central regions of two cool core galaxy clusters, Abell 2597 and Sersic 159-03. Gas is detected out to 20 kpc from the nuclei of the brightest cluster galaxies and found to be distributed in clumps and filaments around it. The ionised and molecular gas phases trace each other closely in extent and dynamical state. Both gas phases show signs of interaction with the active nucleus. Within the nuclear regions the kinetic luminosity of this gas is found to be somewhat smaller than the current radio luminosity. Outside the nuclear region the gas has a low velocity dispersion and shows smooth velocity gradients. There is no strong correlation between the intensity of the molecular and ionised gas emission and either the radio or X-ray emission. The molecular gas in Abell 2597 and Sersic 159-03 is well described by a gas in local thermal equilibrium (LTE) with a single excitation temperature T_exc ~ 2300 K. The emission line ratios do not vary strongly as function of position, with the exception of the nuclear regions where the ionised to molecular gas ratio is found decrease. These constant line ratios imply a single source of heating and excitation for both gas phases.

preprint2009arXiv

The discovery of diffuse steep spectrum sources in Abell 2256

Context: Hierarchical galaxy formation models indicate that during their lifetime galaxy clusters undergo several mergers. Here we report on the discovery of three diffuse radio sources in the periphery of Abell 2256, using the GMRT. Aims: The aim of the observations was to search for diffuse ultra-steep spectrum radio sources within the galaxy cluster Abell 2256. Methods: We have carried out GMRT 325 MHz radio continuum observations of Abell 2256. V, R and I band images of the cluster were taken with the 4.2m WHT. Results: We have discovered three diffuse elongated radio sources located about 1 Mpc from the cluster center. Two are located to the west of the cluster center, and one to the southeast. The sources have a measured physical extent of 170, 140 and 240 kpc, respectively. The two western sources are also visible in deep low-resolution 115-165 MHz WSRT images, although they are blended into a single source. For the combined emission of the blended source we find an extreme spectral index of -2.05\pm 0.14 between 140 and 351 MHz. The extremely steep spectral index suggests these two sources are most likely the result of adiabatic compression of fossil radio plasma due to merger shocks. Conclusions: The discovery of the steep spectrum sources implies the existence of a population of faint diffuse radio sources in (merging) clusters with such steep spectra that they have gone unnoticed in higher frequency (\gtrsim 1 GHz) observations. An exciting possibility therefore is that such sources will determine the general appearance of clusters in low-frequency high resolution radio maps as will be produced by for example LOFAR or LWA.