We present MeerKAT 1000 MHz and 1400 MHz observations of a bright radio galaxy in the southern hemisphere, ESO 137-006. The galaxy lies at the centre of the massive and merging Norma galaxy cluster. ...The MeerKAT continuum images (rms ∼0.02 mJy beam
−1
at ∼10″ resolution) reveal new features that have never been seen in a radio galaxy before: collimated synchrotron threads of yet unknown origin, which link the extended and bent radio lobes of ESO 137-006. The most prominent of these threads stretches in projection for about 80 kpc and is about 1 kpc in width. The radio spectrum of the threads is steep, with a spectral index of up to
α
≃ 2 between 1000 and 1400 MHz.
We report a tight linear relation between the H i circular velocity measured at 6 R
e and the stellar velocity dispersion measured within 1 R
e for a sample of 16 early-type galaxies with stellar ...mass between 1010 and 1011 M⊙. The key difference from previous studies is that we only use spatially resolved v
circ(H i) measurements obtained at large radius for a sizeable sample of objects. We can therefore link a kinematical tracer of the gravitational potential in the dark-matter dominated outer regions of galaxies with one in the inner regions, where baryons control the distribution of mass. We find that v
circ(H i)= 1.33 σe with an observed scatter of just 12 per cent. This indicates a strong coupling between luminous and dark matter from the inner- to the outer regions of early-type galaxies, analogous to the situation in spirals and dwarf irregulars. The v
circ(H i)–σe relation is shallower than those based on v
circ measurements obtained from stellar kinematics and modelling at smaller radius, implying that v
circ declines with radius – as in bulge-dominated spirals. Indeed, the value of v
circ(H i) is typically 25 per cent lower than the maximum v
circ derived at ∼0.2 R
e from dynamical models. Under the assumption of power-law total density profiles ρ ∝ r
−γ, our data imply an average logarithmic slope 〈γ〉 = 2.18 ± 0.03 across the sample, with a scatter of 0.11 around this value. The average slope and scatter agree with recent results obtained from stellar kinematics alone for a different sample of early-type galaxies.
We present new observations of Fornax A taken at ∼1 GHz with the MeerKAT telescope and at ∼6 GHz with the Sardinia Radio Telescope (SRT). The sensitive (noise ∼16
μ
Jy beam
−1
), high-resolution ...(≲10″) MeerKAT images show that the lobes of Fornax A have a double-shell morphology, where dense filaments are embedded in a diffuse and extended cocoon. We study the spectral properties of these components by combining the MeerKAT and SRT observations with archival data between 84 MHz and 217 GHz. For the first time, we show that multiple episodes of nuclear activity must have formed the extended radio lobes. The modelling of the radio spectrum suggests that the last episode of injection of relativistic particles into the lobes started ∼24 Myr ago and stopped 12 Myr ago. More recently (∼3 Myr ago), a less powerful and short (≲1 Myr) phase of nuclear activity generated the central jets. Currently, the core may be in a new active phase. It appears that Fornax A is rapidly flickering. The dense environment around Fornax A has lead to a complex recent merger history for this galaxy, including mergers spanning a range of gas contents and mass ratios, as shown by the analysis of the galaxy’s stellar- and cold-gas phases. This complex recent history may be the cause of the rapid, recurrent nuclear activity of Fornax A.
We present a systematic study of the extraplanar gas (EPG) in a sample of 15 nearby late-type galaxies at intermediate inclinations using publicly available, deep interferometric H I data from the ...Hydrogen Accretion in LOcal GAlaxieS (HALOGAS) survey. For each system we masked the H I emission coming from the regularly rotating disc and used synthetic datacubes to model the leftover “anomalous” H I flux. Our model consists of a smooth, axisymmetric thick component described by three structural and four kinematical parameters, which are fit to the data via a Markov chain Monte Carlo (MCMC) based Bayesian method. We find that extraplanar H I is nearly ubiquitous in disc galaxies as we fail to detect it in only two of the systems with the poorest spatial resolution. The EPG component encloses ∼5−25% of the total H I mass with a mean value of 14%, and has a typical thickness of a few kpc which is incompatible with expectations based on hydrostatic equilibrium models. The EPG kinematics is remarkably similar throughout the sample, and consists of a lagging rotation with typical vertical gradients of ∼ − 10 km s−1 kpc−1, a velocity dispersion of 15−30 km s−1, and, for most galaxies, a global inflow in both the vertical and radial directions with speeds of 20−30 km s−1. The EPG H I masses are in excellent agreement with predictions from simple models of the galactic fountain that are powered by stellar feedback. The combined effect of photo-ionisation and interaction of the fountain material with the circumgalactic medium can qualitatively explain the kinematics of the EPG, but dynamical models of the galactic fountain are required to fully test this framework.
We present MeerKAT observations of neutral hydrogen gas (H I) in the nearby merger remnant NGC 1316 (Fornax A), the brightest member of a galaxy group which is falling into the Fornax cluster. We ...find H I on a variety of scales, from the galaxy centre to its large-scale environment. For the first time we detect H I at large radii (70–150 kpc in projection), mostly distributed on two long tails associated with the galaxy. Gas in the tails dominates the H I mass of NGC 1316: 7 × 108 M⊙– 14 times more than in previous observations. The total H I mass is comparable to the amount of neutral gas found inside the stellar body, mostly in molecular form. The H I tails are associated with faint optical tidal features thought to be the remnant of a galaxy merger occurred a few billion years ago. They demonstrate that the merger was gas-rich. During the merger, tidal forces pulled some gas and stars out to large radii, where we now detect them in the form of optical tails and, thanks to our new data, H I tails; while torques caused the remaining gas to flow towards the centre of the remnant, where it was converted into molecular gas and fuelled the starburst revealed by the galaxy’s stellar populations. Several of the observed properties of NGC 1316 can be reproduced by a ∼10:1 merger between a dominant, gas-poor early-type galaxy and a smaller, gas-rich spiral occurred 1–3 Gyr ago, likely followed by subsequent accretion of satellite galaxies.
Kinematical parametrizations of disc galaxies, employing emission line observations, are indispensable tools for studying the formation and evolution of galaxies. Future large-scale H i surveys will ...resolve the discs of many thousands of galaxies, allowing a statistical analysis of their disc and halo kinematics, mass distribution and dark matter content. Here, we present an automated procedure which fits tilted-ring models to H i data cubes of individual, well-resolved galaxies. The method builds on the 3D Tilted Ring Fitting Code (TiRiFiC) and is called Fully Automated TiRiFiC (FAT). To assess the accuracy of the code, we apply it to a set of 52 artificial galaxies and 25 real galaxies from the Local Volume H i Survey (LVHIS). Using LVHIS data, we compare our 3D modelling to the 2D modelling methods DiskFit and ROTCUR. A conservative result is that FAT accurately models the kinematics and the morphologies of galaxies with an extent of eight beams across the major axis in the inclination range 20°–90° without the need for priors such as disc inclination. When comparing to 2D methods we find that velocity fields cannot be used to determine inclinations in galaxies that are marginally resolved. We conclude that with the current code tilted-ring models can be produced in a fully automated fashion. This will be essential for future H i surveys, with the Square Kilometre Array and its pathfinders, which will allow us to model the gas kinematics of many thousands of well-resolved galaxies. Performance studies of FAT close to our conservative limits, as well as the introduction of more parametrized models will open up the possibility to study even less resolved galaxies.
Abstract
We present the full panchromatic afterglow light-curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By ...compiling all archival data and reprocessing a subset of it, we have evaluated the impact of differences in data processing or flux determination methods used by different groups and attempted to mitigate these differences to provide a more uniform data set. Simple power-law fits to the uniform afterglow light curve indicate a
t
0.86±0.04
rise, a
t
−1.92±0.12
decline, and a peak occurring at 155 ± 4 days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index (−0.584 ± 0.002) across all epochs. This gives a precise and updated estimate of the electron power-law index,
p
= 2.168 ± 0.004. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized interstellar medium density, <0.01 cm
−3
, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having a magnetic field strength between 10
10.4
and 10
16
G. Our fits to the afterglow data using an analytical model that includes Very Long Baseline Interferometry proper motion from Mooley et al., and a structured jet model that ignores the proper motion, indicates that the proper-motion measurement needs to be considered when seeking an accurate estimate of the viewing angle.
Abstract
We report the detection of stimulated hydrogen radio recombination line (RRL) emission from ionized gas in a
z
= 0.89 galaxy using 580–1670 MHz observations from the MeerKAT Absorption Line ...Survey. The RRL emission originates in a galaxy that intercepts and strongly lenses the radio blazar PKS 1830−211 (
z
= 2.5). This is the second detection of RRLs outside of the local Universe and the first clearly associated with hydrogen. We detect effective H144
α
(and H163
α
) transitions at observed frequencies of 1156 (798) MHz by stacking 17 (27) RRLs with 21
σ
(14
σ
) significance. The RRL emission contains two main velocity components and is coincident in velocity with H
i
21 cm and OH 18 cm absorption. We use the RRL spectral line energy distribution and a Bayesian analysis to constrain the density (
n
e
) and the volume-averaged path length (
ℓ
) of the ionized gas. We determine
log
(
n
e
)
=
2.0
−
0.7
+
1.0
cm
−3
and
log
(
ℓ
)
=
−
0.7
−
1.1
+
1.1
pc toward the northeast (NE) lensed image, likely tracing the diffuse thermal phase of the ionized ISM in a thin disk. Toward the southwest (SW) lensed image, we determine
log
(
n
e
)
=
3.2
−
1.0
+
0.4
cm
−3
and
log
(
ℓ
)
=
−
2.7
−
0.2
+
1.8
pc, tracing gas that is more reminiscent of H scii regions. We estimate a star formation (surface density) rate of Σ
SFR
∼ 0.6
M
⊙
yr
−1
kpc
−2
or SFR ∼ 50
M
⊙
yr
−1
, consistent with a star-forming main-sequence galaxy of
M
⋆
∼ 10
11
M
⊙
. The discovery presented here opens up the possibility of studying ionized gas at high redshifts using RRL observations from current and future (e.g., SKA and ngVLA) radio facilities.
We present MeerKAT neutral hydrogen (H
I
) observations of the Fornax A group, which is likely falling into the Fornax cluster for the first time. Our H
I
image is sensitive to 1.4 × 10
19
atoms cm
...−2
over 44.1 km s
−1
, where we detect H
I
in 10 galaxies and a total of (1.12 ± 0.02) × 10
9
M
⊙
of H
I
in the intra-group medium (IGM). We search for signs of pre-processing in the 12 group galaxies with confirmed optical redshifts that reside within the sensitivity limit of our H
I
image. There are 9 galaxies that show evidence of pre-processing and we classify each galaxy into their respective pre-processing category, according to their H
I
morphology and gas (atomic and molecular) scaling relations. Galaxies that have not yet experienced pre-processing have extended H
I
discs and a high H
I
content with a H
2
-to-H
I
ratio that is an order of magnitude lower than the median for their stellar mass. Galaxies that are currently being pre-processed display H
I
tails, truncated H
I
discs with typical gas fractions, and H
2
-to-H
I
ratios. Galaxies in the advanced stages of pre-processing are the most H
I
deficient. If there is any H
I
, they have lost their outer H
I
disc and efficiently converted their H
I
to H
2
, resulting in H
2
-to-H
I
ratios that are an order of magnitude higher than the median for their stellar mass. The central, massive galaxy in our group (NGC 1316) underwent a 10:1 merger ∼2 Gyr ago and ejected 6.6−11.2 × 10
8
M
⊙
of H
I
, which we detect as clouds and streams in the IGM, some of which form coherent structures up to ∼220 kpc in length. We also detect giant (∼100 kpc) ionised hydrogen (H
α
) filaments in the IGM, likely from cool gas being removed (and subsequently ionised) from an in-falling satellite. The H
α
filaments are situated within the hot halo of NGC 1316 and there are localised regions that contain H
I
. We speculate that the H
α
and multiphase gas is supported by magnetic pressure (possibly assisted by the NGC 1316 AGN), such that the hot gas can condense and form H
I
that survives in the hot halo for cosmological timescales.
NGC 4945 is one of the nearest (D ≈ 3.8 Mpc; 1″ ≈ 19 pc) starburst galaxies. To investigate the structure, dynamics, and composition of the dense nuclear gas of this galaxy, ALMA band 3 (λ ≈ 3−4 mm) ...observations were carried out with ≈2″ resolution. Three HCN and two HCO+ isotopologues, CS, C3H2, SiO, HCO, and CH3C2H were measured. Spectral line imaging demonstrates the presence of a rotating nuclear disk of projected size 10″ × 2″ reaching out to a galactocentric radius of r ≈ 100 pc with position angle PA = 45° ± 2°, inclination i = 75° ± 2° and an unresolved bright central core of size ≲2″. The continuum source, representing mostly free-free radiation from star forming regions, is more compact than the nuclear disk by a linear factor of two but shows the same position angle and is centered 0.′′39 ± 0.′′14 northeast of the nuclear accretion disk defined by H2O maser emission. Near the systemic velocity but outside the nuclear disk, both HCN J = 1 → 0 and CS J = 2 → 1 delineate molecular arms of length ≳15″ (≳285 pc) on opposite sides of the dynamical center. These are connected by a (deprojected) ≈ 0.6 kpc sized molecular bridge, likely a dense gaseous bar seen almost ends-on, shifting gas from the front and back side into the nuclear disk. Modeling this nuclear disk located farther inside (r ≲100 pc) with tilted rings provides a good fit by inferring a coplanar outflow reaching a characteristic deprojected velocity of ≈50 km s−1. All our molecular lines, with the notable exception of CH3 C2H, show significant absorption near the systemic velocity (≈571 km s−1), within the range ≈500–660 km s−1. Apparently, only molecular transitions with low critical H2 density (ncrit ≲ 104 cm−3) do not show absorption. The velocity field of the nuclear disk, derived from CH3 C2H, provides evidence for rigid rotation in the inner few arcseconds and a dynamical mass of Mtot = (2.1 ± 0.2) × 108 M⊙ inside a galactocentric radius of 2.′′45 (≈45 pc), with a significantly flattened rotation curve farther out. Velocity integrated line intensity maps with most pronounced absorption show molecular peak positions up to ≈1.′′5 (≈30 pc) southwest of the continuum peak, presumably due to absorption, which appears to be most severe slightly northeast of the nuclear maser disk. A nitrogen isotope ratio of 14N/15N ≈ 200–450 is estimated. This range of values is much higher then previously reported on a tentative basis. Therefore, because 15N is less abundant than expected, the question for strong 15N enrichment by massive star ejecta in starbursts still remains to be settled.