Showing 1–50 of 80 results for author: Bicknell, G V

Probing the role of self-gravity in clouds impacted by AGN-driven winds

Authors: Ankush Mandal, Dipanjan Mukherjee, Christoph Federrath, Geoffrey V. Bicknell, Nicole P. H. Nesvadba, Andrea Mignone

Abstract: The impact of winds and jet-inflated bubbles driven by active galactic nuclei (AGN) are believed to significantly affect the host galaxy's interstellar medium (ISM) and regulate star formation. To explore this scenario, we perform a suite of hydrodynamic simulations to model the interaction between turbulent star-forming clouds and highly pressurised AGN-driven outflows, focusing on the effects of… ▽ More The impact of winds and jet-inflated bubbles driven by active galactic nuclei (AGN) are believed to significantly affect the host galaxy's interstellar medium (ISM) and regulate star formation. To explore this scenario, we perform a suite of hydrodynamic simulations to model the interaction between turbulent star-forming clouds and highly pressurised AGN-driven outflows, focusing on the effects of self-gravity. Our results demonstrate that the cloudlets fragmented by the wind can become gravitationally bound, significantly increasing their survival time. While external pressurisation leads to a global collapse of the clouds in cases of weaker winds ($10^{42}-10^{43}~{\rm erg~s^{-1}}$), higher-power winds ($10^{44}-10^{45}~{\rm erg~s^{-1}}$) disperse the gas and cause localised collapse of the cloudlets. We also demonstrate that a kinetic energy-dominated wind is more efficient in accelerating and dispersing the gas than a thermal wind with the same power. The interaction can give rise to multi-phase outflows with velocities ranging from a few 100 to several 1000~${\rm km\,s^{-1}}$. The mass outflow rates are tightly correlated with the wind power, which we explain by an ablation-based mass-loss model. Moreover, the velocity dispersion and the virial parameter of the cloud material can increase by up to one order of magnitude through the effect of the wind. Even though the wind can suppress or quench star formation for about 1 Myr during the initial interaction, a substantial number of gravitationally bound dense cloudlets manage to shield themselves from the wind's influence and subsequently undergo rapid gravitational collapse, leading to an enhanced star formation rate (SFR). △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: Accepted for publication in MNRAS. Accepted 2024 May 16. Received 2024 May 16; in original form 2024 January 18

Modelling observable signatures of jet-ISM interaction: thermal emission and gas kinematics

Authors: Moun Meenakshi, Dipanjan Mukherjee, Alexander Y. Wagner, Nicole P. H. Nesvadba, Geoffrey V. Bicknell, Raffaella Morganti, Reinier M. J. Janssen, Ralph S. Sutherland, Ankush Mandal

Abstract: Relativistic jets are believed to have a substantial impact on the gas dynamics and evolution of the interstellar medium (ISM) of their host galaxies. In this paper, we aim to draw a link between the simulations and the observable signatures of jet-ISM interactions by analyzing the emission morphology and gas kinematics resulting from jet-induced shocks in simulated disc and spherical systems. We… ▽ More Relativistic jets are believed to have a substantial impact on the gas dynamics and evolution of the interstellar medium (ISM) of their host galaxies. In this paper, we aim to draw a link between the simulations and the observable signatures of jet-ISM interactions by analyzing the emission morphology and gas kinematics resulting from jet-induced shocks in simulated disc and spherical systems. We find that the jet-induced laterally expanding forward shock of the energy bubble sweeping through the ISM causes large-scale outflows, creating shocked emission and high-velocity dispersion in the entire nuclear regions ($\sim2$ kpcs) of their hosts. The jetted systems exhibit larger velocity widths (> 800 km/s), broader Position-Velocity maps and distorted symmetry in the disc's projected velocities than systems without a jet. We also investigate the above quantities at different inclination angles of the observer with respect to the galaxy. Jets inclined to the gas disc of its host are found to be confined for longer times, and consequently couple more strongly with the disc gas. This results in prominent shocked emission and high-velocity widths, not only along the jet's path, but also in the regions perpendicular to them. Strong interaction of the jet with a gas disc can also distort its morphology. However, after the jets escape their initial confinement, the jet-disc coupling is weakened, thereby lowering the shocked emission and velocity widths. △ Less

Submitted 19 September, 2022; v1 submitted 19 March, 2022; originally announced March 2022.

Comments: Matches the Published version

The extent of ionization in simulations of radio-loud AGNs impacting kpc gas discs

Authors: Moun Meenakshi, Dipanjan Mukherjee, Alexander Y. Wagner, Nicole P. H. Nesvadba, Raffaella Morganti, Reinier M. J. Janssen, Geoffrey V. Bicknell

Abstract: We use the results of relativistic hydrodynamic simulations of jet-ISM interactions in a galaxy with a radio-loud AGN to quantify the extent of ionization in the central few kpcs of the gaseous galactic disc. We perform post-process radiative transfer of AGN radiation through the simulated gaseous jet-perturbed disc to estimate the extent of photo-ionization by the AGN with an incident luminosity… ▽ More We use the results of relativistic hydrodynamic simulations of jet-ISM interactions in a galaxy with a radio-loud AGN to quantify the extent of ionization in the central few kpcs of the gaseous galactic disc. We perform post-process radiative transfer of AGN radiation through the simulated gaseous jet-perturbed disc to estimate the extent of photo-ionization by the AGN with an incident luminosity of $10^{45}~\mathrm{erg\,s^{-1}}$. We also map the gas that is collisionally ionized due to shocks driven by the jet. The analysis was carried out for simulations with similar jet power ($10^{45}~\mathrm{erg\,s^{-1}}$) but different jet orientations with respect to the gas disc. We find that the shocks from the jets can ionize a significant fraction (up to 33$\%$) of dense gas ($n>100\,\mathrm{cm^{-3}}$) in the disc, and that the jets clear out the central regions of gas for AGN radiation to penetrate to larger distances in the disc. Jets inclined towards the disc plane couple more strongly with the ISM and ionize a larger fraction of gas in the disc as compared to the vertical jet. However, similar to previous studies, we find that the AGN radiation is quickly absorbed by the outer layers of dense clouds in the disc, and is not able to substantially ionize the disc on a global scale. Thus, compared to jet-ISM interactions, we expect that photo-ionization by the AGN radiation only weakly affects the star-formation activity in the central regions of the galactic disc ($\lesssim 1$ kpc), although the jet-induced shocks can spread farther out. △ Less

Submitted 8 February, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: Published in MNRAS

Tracing the Ionization Structure of the Shocked Filaments of NGC 6240

Authors: Anne M. Medling, Lisa J. Kewley, Daniela Calzetti, George C. Privon, Kirsten Larson, Jeffrey A. Rich, Lee Armus, Mark G. Allen, Geoffrey V. Bicknell, Tanio Díaz-Santos, Timothy M. Heckman, Claus Leitherer, Claire E. Max, David S. N. Rupke, Ezequiel Treister, Hugo Messias, Alexander Y. Wagner

Abstract: We study the ionization and excitation structure of the interstellar medium in the late-stage gas-rich galaxy merger NGC 6240 using a suite of emission line maps at $\sim$25 pc resolution from the Hubble Space Telescope, Keck NIRC2 with Adaptive Optics, and ALMA. NGC 6240 hosts a superwind driven by intense star formation and/or one or both of two active nuclei; the outflows produce bubbles and fi… ▽ More We study the ionization and excitation structure of the interstellar medium in the late-stage gas-rich galaxy merger NGC 6240 using a suite of emission line maps at $\sim$25 pc resolution from the Hubble Space Telescope, Keck NIRC2 with Adaptive Optics, and ALMA. NGC 6240 hosts a superwind driven by intense star formation and/or one or both of two active nuclei; the outflows produce bubbles and filaments seen in shock tracers from warm molecular gas (H$_2$ 2.12$μ$m) to optical ionized gas ([O III], [N II], [S II], [O I]) and hot plasma (Fe XXV). In the most distinct bubble, we see a clear shock front traced by high [O III]/H$β$ and [O III]/[O I]. Cool molecular gas (CO(2-1)) is only present near the base of the bubble, towards the nuclei launching the outflow. We interpret the lack of molecular gas outside the bubble to mean that the shock front is not responsible for dissociating molecular gas, and conclude that the molecular clouds are partly shielded and either entrained briefly in the outflow, or left undisturbed while the hot wind flows around them. Elsewhere in the galaxy, shock-excited H$_2$ extends at least $\sim$4 kpc from the nuclei, tracing molecular gas even warmer than that between the nuclei, where the two galaxies' interstellar media are colliding. A ridgeline of high [O III]/H$β$ emission along the eastern arm aligns with the south nucleus' stellar disk minor axis; optical integral field spectroscopy from WiFeS suggests this highly ionized gas is centered at systemic velocity and likely photoionized by direct line-of-sight to the south AGN. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 27 pages, 18 figures; accepted for publication in ApJ

Resolved simulations of jet-ISM interaction: Implications for gas dynamics and star formation

Authors: Dipanjan Mukherjee, Geoffrey V. Bicknell, Alexander Y. Wagner

Abstract: Relativistic jets can interact with the ambient gas distribution of the host galaxy, before breaking out to larger scales. In the past decade several studies have simulated jet-driven outflows to understand how they affect the nearby environment, and over what spatial and temporal scales such interactions occur. The simulations are able to capture the interaction of the jets with the turbulent clu… ▽ More Relativistic jets can interact with the ambient gas distribution of the host galaxy, before breaking out to larger scales. In the past decade several studies have simulated jet-driven outflows to understand how they affect the nearby environment, and over what spatial and temporal scales such interactions occur. The simulations are able to capture the interaction of the jets with the turbulent clumpy interstellar medium and the resultant energetics of the gas. In this review we summarise the results of such recent studies and discuss their implications on the evolution of the dynamics of the gas distribution and the star formation rate. △ Less

Submitted 22 October, 2021; originally announced October 2021.

Comments: To appear in Astronomische Nachrichten, proceedings paper from "6th Workshop on Compact Steep Spectrum and GHz-Peaked Spectrum Radio Sources"

Impact of relativistic jets on the star formation rate: a turbulence-regulated framework

Authors: Ankush Mandal, Dipanjan Mukherjee, Christoph Federrath, Nicole P. H. Nesvadba, Geoffrey V. Bicknell, Alexander Y. Wagner, Moun Meenakshi

Abstract: We apply a turbulence-regulated model of star formation to calculate the star formation rate (SFR) of dense star-forming clouds in simulations of jet-ISM interactions. The method isolates individual clumps and accounts for the impact of virial parameter and Mach number of the clumps on the star formation activity. This improves upon other estimates of the SFR in simulations of jet--ISM interaction… ▽ More We apply a turbulence-regulated model of star formation to calculate the star formation rate (SFR) of dense star-forming clouds in simulations of jet-ISM interactions. The method isolates individual clumps and accounts for the impact of virial parameter and Mach number of the clumps on the star formation activity. This improves upon other estimates of the SFR in simulations of jet--ISM interactions, which are often solely based on local gas density, neglecting the impact of turbulence. We apply this framework to the results of a suite of jet-ISM interaction simulations to study how the jet regulates the SFR both globally and on the scale of individual star-forming clouds. We find that the jet strongly affects the multi-phase ISM in the galaxy, inducing turbulence and increasing the velocity dispersion within the clouds. This causes a global reduction in the SFR compared to a simulation without a jet. The shocks driven into clouds by the jet also compress the gas to higher densities, resulting in local enhancements of the SFR. However, the velocity dispersion in such clouds is also comparably high, which results in a lower SFR than would be observed in galaxies with similar gas mass surface densities and without powerful radio jets. We thus show that both local negative and positive jet feedback can occur in a single system during a single jet event, and that the star-formation rate in the ISM varies in a complicated manner that depends on the strength of the jet-ISM coupling and the jet break-out time-scale. △ Less

Submitted 26 October, 2021; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: Published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 4, December 2021, Pages 4738-4757

Tracing the Milky Way's Vestigial Nuclear Jet

Authors: Gerald Cecil, Alexander Y. Wagner, Joss Bland-Hawthorn, Geoffrey V. Bicknell, Dipanjan Mukherjee

Abstract: MeerKAT radio continuum and XMM-Newton X-ray images have recently revealed a spectacular bipolar channel at the Galactic Center that spans several degrees ($\sim$0.5 kpc). An intermittent jet likely formed this channel and is consistent with earlier evidence of a sustained, Seyfert-level outburst fueled by black-hole accretion onto Sgr A$^*$ several Myr ago. Therefore, to trace a now weak jet that… ▽ More MeerKAT radio continuum and XMM-Newton X-ray images have recently revealed a spectacular bipolar channel at the Galactic Center that spans several degrees ($\sim$0.5 kpc). An intermittent jet likely formed this channel and is consistent with earlier evidence of a sustained, Seyfert-level outburst fueled by black-hole accretion onto Sgr A$^*$ several Myr ago. Therefore, to trace a now weak jet that perhaps penetrated, deflected, and percolated along multiple paths through the interstellar medium, relevant interactions are identified and quantified in archival X-ray images, Hubble Space Telescope Paschen $α$ images and Atacama Large Millimeter/submillimeter Array millimeter-wave spectra, and new SOAR telescope IR spectra. Hydrodynamical simulations are used to show how a currently weak jet can explain these structures and inflate the ROSAT/eROSITA X-ray and Fermi $γ$-ray bubbles that extend $\pm75°$ from the Galactic plane. Thus, our Galactic outflow has features in common with energetic, jet-driven structures in the prototypical Seyfert galaxy NGC 1068. △ Less

Submitted 7 December, 2021; v1 submitted 2 September, 2021; originally announced September 2021.

Comments: 30 pages, 32 figures in main text, 17 interactive figure and control Javascript in /anc subdirectory. Please retrieve paper with full resolution figures and interactive figure at http://physics.unc.edu/~cecil/MWjet; paper text and figures revised to agree with published version in the ApJ except somehow they managed to completely botch Fig 1 so the correct version is in this archive

Unravelling the enigmatic ISM conditions in Minkowski's Object

Authors: Henry R. M. Zovaro, Robert Sharp, Nicole P. H. Nesvadba, Lisa Kewley, Ralph Sutherland, Philip Taylor, Brent Groves, Alexander Y. Wagner, Dipanjan Mukherjee, Geoffrey V. Bicknell

Abstract: Local examples of jet-induced star formation lend valuable insight into its significance in galaxy evolution and can provide important observational constraints for theoretical models of positive feedback. Using optical integral field spectroscopy, we present an analysis of the ISM conditions in Minkowski's Object ($z = 0.0189$), a peculiar star-forming dwarf galaxy located in the path of a radio… ▽ More Local examples of jet-induced star formation lend valuable insight into its significance in galaxy evolution and can provide important observational constraints for theoretical models of positive feedback. Using optical integral field spectroscopy, we present an analysis of the ISM conditions in Minkowski's Object ($z = 0.0189$), a peculiar star-forming dwarf galaxy located in the path of a radio jet from the galaxy NGC 541. Full spectral fitting with PPXF indicates that Minkowski's Object primarily consists of a young stellar population $\sim 10$ Myr old, confirming that the bulk of the object's stellar mass formed during a recent jet interaction. Minkowski's Object exhibits line ratios largely consistent with star formation, although there is evidence for a low level ($\lesssim 15$ per cent) of contamination from a non-stellar ionising source. Strong-line diagnostics reveal a significant variation in the gas-phase metallicity within the object, with $\log\left( \rm O / H \right) + 12$ varying by $\sim 0.5$ dex, which cannot be explained by in-situ star formation, an enriched outflow from the jet, or enrichment of gas in the stellar bridge between NGC 541 and NGC 545/547. We hypothesise that Minkowski's Object either (a) was formed as a result of jet-induced star formation in pre-existing gas clumps in the stellar bridge, or (b) is a gas-rich dwarf galaxy that is experiencing an elevation in its star formation rate due to a jet interaction, and will eventually redden and fade, becoming an ultra-diffuse galaxy as it is processed by the cluster. △ Less

Submitted 6 October, 2020; originally announced October 2020.

Comments: 22 pages, 18 figures. Accepted to the Monthly Notices of the Royal Astronomical Society

Gas, dust, and star formation in the positive AGN feedback candidate 4C 41.17 at z=3.8

Authors: N. P. H. Nesvadba, G. V. Bicknell, D. Mukherjee, A. Y. Wagner

Abstract: We present new, spatially resolved [CI]1-0, [CI]2-1, CO(7-6), and dust continuum observations of 4C 41.17 at $z=3.8$ obtained with the IRAM NOEMA interferometer. This is one of the best-studied radio galaxies in this epoch and is arguably the best candidate of jet-triggered star formation at high redshift currently known in the literature. 4C 41.17 shows a narrow ridge of dust continuum extending… ▽ More We present new, spatially resolved [CI]1-0, [CI]2-1, CO(7-6), and dust continuum observations of 4C 41.17 at $z=3.8$ obtained with the IRAM NOEMA interferometer. This is one of the best-studied radio galaxies in this epoch and is arguably the best candidate of jet-triggered star formation at high redshift currently known in the literature. 4C 41.17 shows a narrow ridge of dust continuum extending over 15 kpc near the radio jet axis. Line emission is found within the galaxy in the region with signatures of positive feedback. Using the [CI]1-0 line as a molecular gas tracer, and multifrequency observations of the far-infrared dust heated by star formation, we find a total gas mass of $7.6\times 10^{10}$ M$_{\odot}$, which is somewhat greater than that previously found from CO(4-3). The gas mass surface density of $10^3$ M$_{\odot}$ yr$^{-1}$ pc$^{-2}$ and the star formation rate surface density of 10 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$ were derived over the 12 kpc$\times$8 kpc area, where signatures of positive feedback have previously been found. These densities are comparable to those in other populations of massive, dusty star-forming galaxies in this redshift range, suggesting that the jet does not currently enhance the efficiency with which stars form from the gas. This is consistent with expectations from simulations, whereby radio jets may facilitate the onset of star formation in galaxies without boosting its efficiency over longer timescales, in particular after the jet has broken out of the interstellar medium, as is the case in 4C 41.17. △ Less

Submitted 18 June, 2020; originally announced June 2020.

Comments: A&A in press

Searching for signs of jet-driven negative feedback in the nearby radio galaxy UGC 05771

Authors: Henry R. M. Zovaro, Nicole P. H. Nesvadba, Robert Sharp, Geoffrey V. Bicknell, Brent Groves, Dipanjan Mukherjee, Alexander Y. Wagner

Abstract: Hydrodynamical simulations predict that the jets of young radio sources can inhibit star formation in their host galaxies by injecting heat and turbulence into the interstellar medium (ISM). To investigate jet-ISM interactions in a galaxy with a young radio source, we have carried out a multi-wavelength study of the $z = 0.025$ Compact Steep Spectrum radio source hosted by the early-type galaxy UG… ▽ More Hydrodynamical simulations predict that the jets of young radio sources can inhibit star formation in their host galaxies by injecting heat and turbulence into the interstellar medium (ISM). To investigate jet-ISM interactions in a galaxy with a young radio source, we have carried out a multi-wavelength study of the $z = 0.025$ Compact Steep Spectrum radio source hosted by the early-type galaxy UGC 05771. Using Keck/OSIRIS observations, we detected H\textsubscript{2} 1--0 S(1) and [Fe \textsc{ii}] emission at radii of 100s of pc, which traces shocked molecular and ionised gas being accelerated outwards by the jets to low velocities, creating a `stalling wind'. At kpc radii, we detected shocked ionised gas using observations from the CALIFA survey, covering an area much larger than the pc-scale radio source. We found that existing interferometric radio observations fail to recover a large fraction of the source's total flux, indicating the likely existence of jet plasma on kpc scales, which is consistent with the extent of shocked gas in the host galaxy. To investigate the star formation efficiency in UGC 05771, we obtained IRAM CO observations to analyse the molecular gas properties. We found that UGC 05771 sits below the Kennicutt-Schmidt relation, although we were unable to definitively conclude if direct interactions from the jets are inhibiting star formation. This result shows that jets may be important in regulating star formation in the host galaxies of compact radio sources. △ Less

Submitted 31 August, 2019; originally announced September 2019.

Comments: 19 pages, 15 figures. Accepted to the Monthly Notices of the Royal Astronomical Society

Jets blowing bubbles in the young radio galaxy 4C 31.04

Authors: Henry R. M. Zovaro, Robert Sharp, Nicole P. H. Nesvadba, Geoffrey V. Bicknell, Dipanjan Mukherjee, Alexander Y. Wagner, Brent Groves, Shreyam Krishna

Abstract: We report the discovery of shocked molecular and ionized gas resulting from jet-driven feedback in the compact radio galaxy 4C 31.04 using near-IR imaging spectroscopy. 4C 31.04 is a $\sim 100$ pc double-lobed Compact Steep Spectrum source believed to be a very young AGN. It is hosted by a giant elliptical with a $\sim 10^{9}~\rm M_\odot$ multi-phase gaseous circumnuclear disc. We used high spatia… ▽ More We report the discovery of shocked molecular and ionized gas resulting from jet-driven feedback in the compact radio galaxy 4C 31.04 using near-IR imaging spectroscopy. 4C 31.04 is a $\sim 100$ pc double-lobed Compact Steep Spectrum source believed to be a very young AGN. It is hosted by a giant elliptical with a $\sim 10^{9}~\rm M_\odot$ multi-phase gaseous circumnuclear disc. We used high spatial resolution, adaptive optics-assisted $H$- and $K$-band integral field Gemini/NIFS observations to probe (1) the warm ($\sim 10^3~\rm K$) molecular gas phase, traced by ro-vibrational transitions of $\rm H_2$, and (2), the warm ionized medium, traced by the [Fe II]$_{1.644~\rm μm}$ line. The [Fe II] emission traces shocked gas ejected from the disc plane by a jet-blown bubble $300-400~\rm pc$ in diameter, whilst the $\rm H_2$ emission traces shock-excited molecular gas in the interior $\sim 1~\rm kpc$ of the circumnuclear disc. Hydrodynamical modelling shows that the apparent discrepancy between the extent of the shocked gas and the radio emission can occur when the brightest regions of the synchrotron-emitting plasma are temporarily halted by dense clumps, whilst less bright plasma can percolate through the porous ISM and form an energy-driven bubble that expands freely out of the disc plane. This bubble is filled with low surface-brightness plasma not visible in existing VLBI observations of 4C 31.04 due to insufficient sensitivity. Additional radial flows of jet plasma may percolate to $\sim \rm kpc$ radii in the circumnuclear disc, driving shocks and accelerating clouds of gas, giving rise to the $\rm H_2$ emission. △ Less

Submitted 20 February, 2019; v1 submitted 21 November, 2018; originally announced November 2018.

Comments: 18 pages, 10 figures. Accepted to the Monthly Notices of the Royal Astronomical Society

How frequent are close supermassive binary black holes in powerful jet sources?

Authors: Martin G. H. Krause, Stanislav S. Shabala, Martin J. Hardcastle, Geoffrey V. Bicknell, Hans Böhringer, Gayoung Chon, Mohammad A. Nawaz, Marc Sarzi, Alexander Y. Wagner

Abstract: Supermassive black hole binaries may be detectable by an upcoming suite of gravitational wave experiments. Their binary nature can also be revealed by radio jets via a short-period precession driven by the orbital motion as well as the geodetic precession at typically longer periods. We have investigated Karl G. Jansky Very Large Array (VLA) and MERLIN radio maps of powerful jet sources for morpho… ▽ More Supermassive black hole binaries may be detectable by an upcoming suite of gravitational wave experiments. Their binary nature can also be revealed by radio jets via a short-period precession driven by the orbital motion as well as the geodetic precession at typically longer periods. We have investigated Karl G. Jansky Very Large Array (VLA) and MERLIN radio maps of powerful jet sources for morphological evidence of geodetic precession. For perhaps the best studied source, Cygnus A, we find strong evidence for geodetic precession. Projection effects can enhance precession features, for which we find indications in strongly projected sources. For a complete sample of 33 3CR radio sources we find strong evidence for jet precession in 24 cases (73 per cent). The morphology of the radio maps suggests that the precession periods are of the order of 10^6 - 10^7 yr. We consider different explanations for the morphological features and conclude that geodetic precession is the best explanation. The frequently observed gradual jet angle changes in samples of powerful blazars can be explained by orbital motion. Both observations can be explained simultaneously by postulating that a high fraction of powerful radio sources have sub-parsec supermassive black hole binaries. We consider complementary evidence and discuss if any jetted supermassive black hole with some indication of precession could be detected as individual gravitational wave source in the near future. This appears unlikely, with the possible exception of M87. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 24 pages, 36 figures, accepted by MNRAS

Relativistic jet feedback III: feedback on gas disks

Authors: Dipanjan Mukherjee, Geoffrey V. Bicknell, Alexander Y. Wagner, Ralph S. Sutherland, Joseph Silk

Abstract: We study the interactions of a relativistic jet with a dense turbulent gaseous disk of radius $\sim 2$ kpc. We have performed a suite of simulations with different mean density, jet power and orientation. Our results show that: A) the relativistic jet couples strongly with the gas in the inner kpc, creating a cavity and launching outflows. B) The high pressure bubble inflated by the jet and its ba… ▽ More We study the interactions of a relativistic jet with a dense turbulent gaseous disk of radius $\sim 2$ kpc. We have performed a suite of simulations with different mean density, jet power and orientation. Our results show that: A) the relativistic jet couples strongly with the gas in the inner kpc, creating a cavity and launching outflows. B) The high pressure bubble inflated by the jet and its back-flow compresses the disk at the outer edges, driving inflows within the disk. C) Jets inclined towards the disk couple more strongly with the disk and launch sub-relativistic, wide-angle outflows along the minor axis. D) Shocks driven directly by the jet and the jet-driven energy bubble raise the velocity dispersion throughout the disk by several times its initial value. E) Compression by the jet-driven shocks can enhance the star formation rate in the disk, especially in a ring-like geometry close to the axis. However, enhanced turbulent dispersion in the disk also leads to quenching of star formation. Whether positive or negative feedback dominates depends on jet power, ISM density, jet orientation with respect to the disc, and the time-scale under consideration. Qualitatively, our simulations compare favourably with kinematic and morphological signatures of several observed galaxies such as NGC 1052, NGC 3079, 3C 326 and 3C 293. △ Less

Submitted 3 July, 2018; v1 submitted 22 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS

The jet-ISM interactions in IC 5063

Authors: Dipanjan Mukherjee, Alex Y. Wagner, Geoffrey V. Bicknell, Raffaella Morganti, Tom Oosterloo, Nicole Nesvadba, Ralph S. Sutherland

Abstract: The interstellar medium of the radio galaxy IC 5063 is highly perturbed by an AGN jet expanding in the gaseous disc of the galaxy. We model this interaction with relativistic hydrodynamic simulations and multiphase initial conditions for the interstellar medium and compare the results with recent observations. As the jets flood through the inter-cloud channels of the disc, they ablate, accelerate,… ▽ More The interstellar medium of the radio galaxy IC 5063 is highly perturbed by an AGN jet expanding in the gaseous disc of the galaxy. We model this interaction with relativistic hydrodynamic simulations and multiphase initial conditions for the interstellar medium and compare the results with recent observations. As the jets flood through the inter-cloud channels of the disc, they ablate, accelerate, and disperse clouds to velocities exceeding $400 \mbox{km s}^{-1}$. Clouds are also destroyed or displaced in bulk from the central regions of the galaxy. Our models with jet powers of $10^{44} \mbox{erg s}^{-1}$ and $10^{45} \mbox{erg s}^{-1}$ are capable of reproducing many of the features seen in the position-velocity diagram published in Morganti et al. (2015) and confirm the notion that the jet is responsible for the strongly perturbed gas dynamics seen in the ionized, neutral, and molecular gas phases. In our simulations, we also see strong venting of the jet plasma perpendicular to the disc, which entrains clumps and diffuse filaments into the halo of the galaxy. Our simulations are the first 3D hydrodynamic simulations of the jet and ISM of IC 5063. △ Less

Submitted 21 January, 2018; originally announced January 2018.

Comments: Accepted for publication in MNRAS

Relativistic jet feedback II: Relationship to gigahertz peak spectrum and compact steep spectrum radio galaxies

Authors: Geoffrey V. Bicknell, Dipanjan Mukherjee, Alexander Y. Wagner, Ralph S. Sutherland, Nicole P. H. Nesvadba

Abstract: We propose that Gigahertz Peak Spectrum (GPS) and Compact Steep Spectrum (CSS) radio sources are the signposts of relativistic jet feedback in evolving galaxies. Our simulations of relativistic jets interacting with a warm, inhomogeneous medium, utilize cloud densities and velocity dispersions in the range derived from optical observations, show that free-free absorption can account for the… ▽ More We propose that Gigahertz Peak Spectrum (GPS) and Compact Steep Spectrum (CSS) radio sources are the signposts of relativistic jet feedback in evolving galaxies. Our simulations of relativistic jets interacting with a warm, inhomogeneous medium, utilize cloud densities and velocity dispersions in the range derived from optical observations, show that free-free absorption can account for the $\sim \rm GHz$ peak frequencies and low frequency power laws inferred from the radio observations. These new computational models replace the power-law model for the free-free optical depth in the \citep{bicknell97a} model by a more fundamental model involving disrupted log-normal distributions of warm gas. One feature of our new models is that at early stages, the low frequency spectrum is steep but progressively flattens as a result of a broader distribution of optical depths, suggesting that the steep low frequency spectra discovered by \citet{callingham17a} may possibly be attributed to young sources. We also investigate the inverse correlation between peak frequency and size and find that the initial location on this correlation is determined by the average density of the warm ISM. The simulated sources track this correlation initially but eventually fall below it, indicating the need for a more extended ISM than presently modelled. GPS and CSS sources can potentially provide new insights into the phenomenon of AGN feedback since their peak frequencies and spectra are indicative of the density, turbulent structure and distribution of gas in the host galaxy. △ Less

Submitted 15 January, 2018; originally announced January 2018.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Centrifugally driven winds from protostellar accretion discs. I - Formulation and initial results

Authors: C. A. Nolan, R. Salmeron, C. Federrath, G. V. Bicknell, R. S. Sutherland

Abstract: Protostellar discs play an important role in star formation, acting as the primary mass reservoir for accretion onto young stars and regulating the extent to which angular momentum and gas is released back into stellar nurseries through the launching of powerful disc winds. In this study, we explore how disc structure relates to the properties of the wind-launching region, mapping out the regions… ▽ More Protostellar discs play an important role in star formation, acting as the primary mass reservoir for accretion onto young stars and regulating the extent to which angular momentum and gas is released back into stellar nurseries through the launching of powerful disc winds. In this study, we explore how disc structure relates to the properties of the wind-launching region, mapping out the regions of protostellar discs where wind launching could be viable. We combine a series of 1.5D semi-analytic, steady-state, vertical disc-wind solutions into a radially extended 1+1.5D model, incorporating all three diffusion mechanisms (Ohm, Hall and ambipolar). We observe that the majority of mass outflow via disc winds occurs over a radial width of a fraction of an astronomical unit, with outflow rates attenuating rapidly on either side. We also find that the mass accretion rate, magnetic field strength and surface density profile each have significant effects on both the location of the wind-launching region and the ejection/accretion ratio $\dot{M}_{\rm out}/\dot{M}_{\rm in}$. Increasing either the accretion rate or the magnetic field strength corresponds to a shift of the wind-launching region to smaller radii and a decrease in $\dot{M}_{\rm out}/\dot{M}_{\rm in}$, while increasing the surface density corresponds to launching regions at larger radii with increased $\dot{M}_{\rm out}/\dot{M}_{\rm in}$. Finally, we discover a class of disc winds containing an ineffective launching configuration at intermediate radii, leading to two radially separated regions of wind launching and diminished $\dot{M}_{\rm out}/\dot{M}_{\rm in}$. We find that the wind locations and ejection/accretion ratio are consistent with current observational and theoretical estimates. △ Less

Submitted 4 July, 2017; originally announced July 2017.

Comments: 19 pages, 11 figures, accepted for publication in MNRAS

Relativistic jet feedback in high-redshift galaxies I: Dynamics

Authors: Dipanjan Mukherjee, Geoffrey V. Bicknell, Ralph S. Sutherland, A. Y. Wagner

Abstract: We present the results of three dimensional relativistic hydrodynamic simulations of interaction of AGN jets with a dense turbulent two-phase interstellar medium, which would be typical of high redshift galaxies. We describe the effect of the jet on the evolution of the density of the turbulent ISM. The jet driven energy bubble affects the gas to distances up to several kiloparsecs from the inject… ▽ More We present the results of three dimensional relativistic hydrodynamic simulations of interaction of AGN jets with a dense turbulent two-phase interstellar medium, which would be typical of high redshift galaxies. We describe the effect of the jet on the evolution of the density of the turbulent ISM. The jet driven energy bubble affects the gas to distances up to several kiloparsecs from the injection region. The shocks resulting from such interactions create a multi-phase ISM and radial outflows. One of the striking result of this work is that low power jets (P_jet < 10^{43} erg/s) although less efficient in accelerating clouds, are trapped in the ISM for a longer time and hence affect the ISM over a larger volume. Jets of higher power drill through with relative ease. Although the relativistic jets launch strong outflows, there is little net mass ejection to very large distances, supporting a galactic fountain scenario for local feedback. △ Less

Submitted 3 June, 2016; originally announced June 2016.

Comments: Accepted for publication in MNRAS

Jet-Intracluster Medium interaction in Hydra A. II The Effect of Jet Precession

Authors: M. A. Nawaz, G. V. Bicknell, A. Y. Wagner, R. S. Sutherland, B. R. McNamara

Abstract: We present three dimensional relativistic hydrodynamical simulations of a precessing jet interacting with the intracluster medium and compare the simulated jet structure with the observed structure of the Hydra A northern jet. For the simulations, we use jet parameters obtained in the parameter space study of the first paper in this series and probe different values for the precession period and p… ▽ More We present three dimensional relativistic hydrodynamical simulations of a precessing jet interacting with the intracluster medium and compare the simulated jet structure with the observed structure of the Hydra A northern jet. For the simulations, we use jet parameters obtained in the parameter space study of the first paper in this series and probe different values for the precession period and precession angle. We find that for a precession period P = 1 Myr and a precession angle = 20 degree the model reproduces i) the curvature of the jet, ii) the correct number of bright knots within 20 kpc at approximately correct locations, and iii) the turbulent transition of the jet to a plume. The Mach number of the advancing bow shock = 1.85 is indicative of gentle cluster atmosphere heating during the early stages of the AGN's activity. △ Less

Submitted 9 February, 2016; originally announced February 2016.

Comments: 15 pages, 11 figures, 2 tables

Galaxy-scale AGN Feedback - Theory

Authors: A. Y. Wagner, G. V. Bicknell, M. Umemura, R. S. Sutherland, J. Silk

Abstract: Powerful relativistic jets in radio galaxies are capable of driving strong outflows but also inducing star-formation by pressure-triggering collapse of dense clouds. We review theoretical work on negative and positive active galactic nuclei feedback, discussing insights gained from recent hydrodynamical simulations of jet-driven feedback on galaxy scales that are applicable to compact radio source… ▽ More Powerful relativistic jets in radio galaxies are capable of driving strong outflows but also inducing star-formation by pressure-triggering collapse of dense clouds. We review theoretical work on negative and positive active galactic nuclei feedback, discussing insights gained from recent hydrodynamical simulations of jet-driven feedback on galaxy scales that are applicable to compact radio sources. The simulations show that the efficiency of feedback and the relative importance of negative and positive feedback depends strongly on interstellar medium properties, especially the column depth and spatial distribution of clouds. Negative feedback is most effective if clouds are distributed spherically and individual clouds have small column depths, while positive feedback is most effective if clouds are predominantly in a disc-like configuration. △ Less

Submitted 13 October, 2015; originally announced October 2015.

Comments: 8 pages, 6 figures, Proceedings from the 5th Workshop on Compact Steep Spectrum and GHz-peaked Spectrum Radio Sources

Turbulent mixing layers in supersonic protostellar outflows, with application to DG Tauri

Authors: Marc C. White, Geoffrey V. Bicknell, Ralph S. Sutherland, Raquel Salmeron, Peter J. McGregor

Abstract: Turbulent entrainment processes may play an important role in the outflows from young stellar objects at all stages of their evolution. In particular, lateral entrainment of ambient material by high-velocity, well-collimated protostellar jets may be the cause of the multiple emission-line velocity components observed in the microjet-scale outflows driven by classical T Tauri stars. Intermediate-ve… ▽ More Turbulent entrainment processes may play an important role in the outflows from young stellar objects at all stages of their evolution. In particular, lateral entrainment of ambient material by high-velocity, well-collimated protostellar jets may be the cause of the multiple emission-line velocity components observed in the microjet-scale outflows driven by classical T Tauri stars. Intermediate-velocity outflow components may be emitted by a turbulent, shock- excited mixing layer along the boundaries of the jet. We present a formalism for describing such a mixing layer based on Reynolds decomposition of quantities measuring fundamental properties of the gas. In this model, the molecular wind from large disc radii provides a continual supply of material for entrainment. We calculate the total stress profile in the mixing layer, which allows us to estimate the dissipation of turbulent energy, and hence the luminosity of the layer. We utilize MAPPINGS IV shock models to determine the fraction of total emission that occurs in [Fe II] 1.644 μm line emission in order to facilitate comparison to previous observations of the young stellar object DG Tauri. Our model accurately estimates the luminosity and changes in mass outflow rate of the intermediate-velocity component of the DG Tau approaching outflow. Therefore, we propose that this component represents a turbulent mixing layer surrounding the well-collimated jet in this object. Finally, we compare and contrast our model to previous work in the field. △ Less

Submitted 5 October, 2015; originally announced October 2015.

Comments: 18 pages, 13 figures, accepted for publication in MNRAS

A Unified Model of the Fermi Bubbles, Microwave Haze, and Polarized Radio Lobes: Reverse Shocks in the Galactic Center's Giant Outflows

Authors: Roland M. Crocker, Geoffrey V. Bicknell, Andrew M. Taylor, Ettore Carretti

Abstract: The Galactic Center's giant outflows are manifest in three different, non-thermal phenomena: i) the hard-spectrum, γ-ray `Fermi Bubbles' emanating from the nucleus and extending to |b| ~ 50 degrees; ii) the hard-spectrum, total-intensity microwave (~ 20-40 GHz) `Haze' extending to |b| ~ 35 degrees in the lower reaches of the Fermi Bubbles; and iii) the steep spectrum, polarized, `S-PASS' radio (~… ▽ More The Galactic Center's giant outflows are manifest in three different, non-thermal phenomena: i) the hard-spectrum, γ-ray `Fermi Bubbles' emanating from the nucleus and extending to |b| ~ 50 degrees; ii) the hard-spectrum, total-intensity microwave (~ 20-40 GHz) `Haze' extending to |b| ~ 35 degrees in the lower reaches of the Fermi Bubbles; and iii) the steep spectrum, polarized, `S-PASS' radio (~ 2-20 GHz) Lobes that envelop the Bubbles and extend to |b| ~ 60 degrees. We find that the nuclear outflows inflate a genuine bubble in each Galactic hemisphere that has the classical structure, working outwards, of reverse shock, contact discontinuity, and forward shock. Expanding into the finite pressure of the halo and given appreciable cooling and gravitational losses, the contact discontinuity of each bubble is now expanding only very slowly. We find observational signatures in both hemispheres of giant, reverse shocks at heights of ~ 1 kpc above the nucleus; their presence ultimately explains all three of the non-thermal phenomena mentioned above. Synchrotron emission from shock-reaccelerated cosmic-ray electrons explains the spectrum, morphology, and vertical extent of the microwave Haze and the polarized radio Lobes. Collisions between shock-reaccelerated hadrons and denser gas in cooling condensations that form inside the contact discontinuity account for most of the Bubbles' γ-ray emissivity. △ Less

Submitted 17 June, 2015; v1 submitted 23 December, 2014; originally announced December 2014.

Comments: Accepted for publication in the Astrophysical Journal. 14 figures. Qualitative results unchanged from vers. 1 but calculations made more robust and text clarified

Journal ref: ApJ, 808, 107 (2015)

Multi-dimensional simulations of the expanding supernova remnant of SN 1987A

Authors: T. M. Potter, L. Staveley-Smith, B. Reville, C. -Y. Ng, G. V. Bicknell, R. S. Sutherland, A. Y. Wagner

Abstract: The expanding remnant from SN 1987A is an excellent laboratory for investigating the physics of supernovae explosions. There are still a large number of outstanding questions, such the reason for the asymmetric radio morphology, the structure of the pre-supernova environment, and the efficiency of particle acceleration at the supernova shock. We explore these questions using three-dimensional simu… ▽ More The expanding remnant from SN 1987A is an excellent laboratory for investigating the physics of supernovae explosions. There are still a large number of outstanding questions, such the reason for the asymmetric radio morphology, the structure of the pre-supernova environment, and the efficiency of particle acceleration at the supernova shock. We explore these questions using three-dimensional simulations of the expanding remnant between days 820 and 10,000 after the supernova. We combine a hydrodynamical simulation with semi-analytic treatments of diffusive shock acceleration and magnetic field amplification to derive radio emission as part of an inverse problem. Simulations show that an asymmetric explosion, combined with magnetic field amplification at the expanding shock, is able to replicate the persistent one-sided radio morphology of the remnant. We use an asymmetric Truelove & McKee progenitor with an envelope mass of $10 M_{\sun}$ and an energy of $1.5 \times 10^{44} J$. A termination shock in the progenitor's stellar wind at a distance of $0\farcs43-0\farcs51$ provides a good fit to the turn on of radio emission around day 1200. For the H\textsc{ii} region, a minimum distance of $0\farcs63\pm0\farcs01$ and maximum particle number density of $(7.11\pm1.78) \times 10^7$ m$^{-3}$ produces a good fit to the evolving average radius and velocity of the expanding shocks from day 2000 to day 7000 after explosion. The model predicts a noticeable reduction, and possibly a temporary reversal, in the asymmetric radio morphology of the remnant after day 7000, when the forward shock left the eastern lobe of the equatorial ring. △ Less

Submitted 14 September, 2014; originally announced September 2014.

Comments: 31 pages, 19 figures

Jet-Intracluster Medium interaction in Hydra A. I Estimates of jet velocity from inner knots

Authors: M. A. Nawaz, A. Y. Wagner, G. V. Bicknell, R. S. Sutherland, B. R. McNamara

Abstract: We present the first stage of an investigation of the interactions of the jets in the radio galaxy Hydra A with the intracluster medium. We consider the jet kinetic power, the galaxy and cluster atmosphere, and the inner structure of the radio source. Analysing radio observations of the inner lobes of Hydra A by Taylor et al. (1990) we confirm the jet power estimates of about 1e45 ergs/s derived b… ▽ More We present the first stage of an investigation of the interactions of the jets in the radio galaxy Hydra A with the intracluster medium. We consider the jet kinetic power, the galaxy and cluster atmosphere, and the inner structure of the radio source. Analysing radio observations of the inner lobes of Hydra A by Taylor et al. (1990) we confirm the jet power estimates of about 1e45 ergs/s derived by Wise et al. (2007) from dynamical analysis of the X-ray cavities. With this result and a model for the galaxy halo, we explore the jet-intracluster medium interactions occurring on a scale of 10 kpc using two-dimensional, axisymmetric, relativistic pure hydrodynamic simulations. A key feature is that we identify the three bright knots in the northern jet as biconical reconfinement shocks, which result when an over pressured jet starts to come into equilibrium with the galactic atmosphere. Through an extensive parameter space study we deduce that the jet velocity is approximately 0.8 c at a distance 0.5 kpc from the black hole. The combined constraints of jet power, the observed jet radius profile along the jet, and the estimated jet pressure and jet velocity imply a value of the jet density parameter approximately 13 for the northern jet. We show that for a jet velocity = 0.8c and angle between the jet and the line of sight = 42 deg, an intrinsic asymmetry in the emissivity of the northern and southern jet is required for a consistent brightness ratio approximately 7 estimated from the 6cm VLA image of Hydra A. △ Less

Submitted 19 August, 2014; originally announced August 2014.

Comments: 17 pages, 10 figures, 4 tables

Multi-epoch Sub-arcsecond [Fe II] Spectroimaging of the DG Tau Outflows with NIFS. II. On the Nature of the Bipolar Outflow Asymmetry

Authors: Marc C White, Geoffrey V Bicknell, Peter J McGregor, Raquel Salmeron

Abstract: The origin of bipolar outflow asymmetry in young stellar objects (YSOs) remains poorly understood. It may be due to an intrinsically asymmetric outflow launch mechanism, or it may be caused by the effects of the ambient medium surrounding the YSO. Answering this question is an important step in understanding outflow launching. We have investigated the bipolar outflows driven by the T Tauri star DG… ▽ More The origin of bipolar outflow asymmetry in young stellar objects (YSOs) remains poorly understood. It may be due to an intrinsically asymmetric outflow launch mechanism, or it may be caused by the effects of the ambient medium surrounding the YSO. Answering this question is an important step in understanding outflow launching. We have investigated the bipolar outflows driven by the T Tauri star DG Tauri on scales of hundreds of AU, using the Near-infrared Integral Field Spectrograph (NIFS) on Gemini North. The approaching outflow consists of a well-collimated jet, nested within a lower-velocity disc wind. The receding outflow is composed of a single-component bubble-like structure. We analyse the kinemat- ics of the receding outflow using kinetic models, and determine that it is a quasi-stationary bubble with an expanding internal velocity field. We propose that this bubble forms because the receding counterjet from DG Tau is obstructed by a clumpy ambient medium above the circumstellar disc surface, based on similarities between this structure and those found in the modeling of active galactic nuclei outflows. We find evidence of interaction between the obscured counterjet and clumpy ambient material, which we attribute to the large molecular envelope around the DG Tau system. An analytical model of a momentum-driven bubble is shown to be consistent with our interpretation. We conclude that the bipolar outflow from DG Tau is intrinsically symmetric, and the observed asymmetries are due to environmental effects. This mechanism can potentially be used to explain the observed bipolar asymmetries in other YSO outflows. △ Less

Submitted 23 April, 2014; originally announced April 2014.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Multi-epoch Sub-arcsecond [Fe II] Spectroimaging of the DG Tau Outflows with NIFS. I. First data epoch

Authors: Marc C. White, Peter J. McGregor, Geoffrey V. Bicknell, Raquel Salmeron, Tracy L. Beck

Abstract: Investigating the outflows emanating from young stellar objects (YSOs) on sub-arcsecond scales provides important clues to the nature of the underlying accretion-ejection process occurring near the central protostar. We have investigated the structures and kinematics of the outflows driven by the YSO DG Tauri, using the Near-infrared Integral Field Spectrograph (NIFS) on Gemini North. The blueshif… ▽ More Investigating the outflows emanating from young stellar objects (YSOs) on sub-arcsecond scales provides important clues to the nature of the underlying accretion-ejection process occurring near the central protostar. We have investigated the structures and kinematics of the outflows driven by the YSO DG Tauri, using the Near-infrared Integral Field Spectrograph (NIFS) on Gemini North. The blueshifted outflow shows two distinct components in [Fe II] 1.644 micron emission, which are separated using multi-component line fitting. A stationary recollimation shock is observed, in agreement with previous X-ray and FUV observations. The presence of this shock indicates that the innermost streamlines of the high-velocity component are launched at a very small radius, 0.01-0.15 AU, from the central star. The jet accelerates and expands downstream of the recollimation shock; the 'acceleration' is likely a sign of velocity variations in the jet. No evidence of rotation is found, and we compare this non-detection to previous counter-claims. Moving jet knots, likely the result of the jet velocity variations, are observed. One of these knots moves more slowly than previously observed knots, and the knot ejection interval appears to be non-periodic. An intermediate-velocity component surrounds this central jet, and is interpreted as the result of a turbulent mixing layer along the jet boundaries. Such lateral entrainment requires the presence of a magnetic field of strength a few mG or less at hundreds of AU above the disc surface, which is argued to be a reasonable proposition. In H2 1-0 S(1) 2.1218 micron emission, a wide-angle, intermediate-velocity blueshifted outflow is observed. Both outflows are consistent with being launched by a magnetocentrifugal disc wind, although an X-wind origin for the high-velocity jet cannot be ruled out. △ Less

Submitted 2 April, 2014; originally announced April 2014.

Comments: 28 pages, 19 figures, accepted for publication in MNRAS

Steady-State Hadronic Gamma-Ray Emission from 100-Myr-Old Fermi Bubbles

Authors: Roland M. Crocker, Geoffrey V. Bicknell, Ettore Carretti, Alex S. Hill, Ralph S. Sutherland

Abstract: The Fermi Bubbles are enigmatic γ-ray features of the Galactic bulge. Both putative activity (within $\sim$ few $\times$ Myr) connected to the Galactic center super-massive black hole and, alternatively, nuclear star formation have been claimed as the energising source of the Bubbles. Likewise, both inverse-Compton emission by non-thermal electrons (`leptonic' models) and collisions between non-th… ▽ More The Fermi Bubbles are enigmatic γ-ray features of the Galactic bulge. Both putative activity (within $\sim$ few $\times$ Myr) connected to the Galactic center super-massive black hole and, alternatively, nuclear star formation have been claimed as the energising source of the Bubbles. Likewise, both inverse-Compton emission by non-thermal electrons (`leptonic' models) and collisions between non-thermal protons and gas (`hadronic' models) have been advanced as the process supplying the Bubbles' γ-ray emission. An issue for any steady state hadronic model is that the very low density of the Bubbles' plasma seems to require that they accumulate protons over a multi-Gyr timescale, much longer than other natural timescales occurring in the problem. Here we present a hadronic model where the timescale for generating the Bubbles' hadronic γ-ray emission is $\sim$ few $\times 10^8$ years. Our model invokes collapse of the Bubbles' thermally-unstable plasma, leading to an accumulation of cosmic rays and magnetic field into localised, warm ($\sim 10^4$ K), and likely filamentary condensations of higher density gas. Under the condition that these filaments are supported by non-thermal pressure, we can predict the hadronic emission from the Bubbles to be $L_γ\simeq 2 \times 10^{37}$ erg/s $\ \dot{M}_\mathrm{in}/(0.1 \ M_{Sun}/$ year $) \ T_\mathrm{FB}^2/(3.5 \times 10^7 K) ^2 M_{fil}/M_{pls}$ ; precisely their observed luminosity (normalizing to the star-formation-driven mass flux into the Bubbles and their measured plasma temperature and adopting the further result that the mass in the filaments, $M_{fil}$ is approximately equal to that of the Bubbles' plasma, $M_{pls}$). △ Less

Submitted 17 August, 2014; v1 submitted 2 December, 2013; originally announced December 2013.

Comments: Accepted for publication in ApJ Letters. Shortened from ver.1 with an important modification to equation 2; qualitative results the same; vers 3: title fixed

Journal ref: ApJ, 791, L20 (2014)

Global simulations of magnetorotational turbulence I: convergence and the quasi-steady state

Authors: E. R. Parkin, G. V. Bicknell

Abstract: Magnetorotational turbulence provides a viable mechanism for angular momentum transport in accretion disks. We present global, three dimensional (3D), MHD accretion disk simulations that investigate the dependence of the turbulent stresses on resolution. Convergence in the time-and-volume-averaged stress-to-gas-pressure ratio, at a value of $\sim0.04$, is found for a model with radial, vertical, a… ▽ More Magnetorotational turbulence provides a viable mechanism for angular momentum transport in accretion disks. We present global, three dimensional (3D), MHD accretion disk simulations that investigate the dependence of the turbulent stresses on resolution. Convergence in the time-and-volume-averaged stress-to-gas-pressure ratio, at a value of $\sim0.04$, is found for a model with radial, vertical, and azimuthal resolution of 12-51, 27, and 12.5 cells per scale-height (the simulation mesh is such that cells per scale-height varies in the radial direction). A control volume analysis is performed on the main body of the disk (|z|<2H) to examine the production and removal of magnetic energy. Maxwell stresses in combination with the mean disk rotation are mainly responsible for magnetic energy production, whereas turbulent dissipation (facilitated by numerical resistivity) predominantly removes magnetic energy from the disk. Re-casting the magnetic energy equation in terms of the power injected by Maxwell stresses on the boundaries of, and by Lorentz forces within, the control volume highlights the importance of the boundary conditions (of the control volume). The different convergence properties of shearing-box and global accretion disk simulations can be readily understood on the basis of choice of boundary conditions and the magnetic field configuration. Periodic boundary conditions restrict the establishment of large-scale gradients in the magnetic field, limiting the power that can be delivered to the disk by Lorentz forces and by stresses at the surfaces. The factor of three lower resolution required for convergence in turbulent stresses for our global disk models compared to stratified shearing-boxes is explained by this finding. (Abridged) △ Less

Submitted 8 July, 2013; v1 submitted 5 June, 2013; originally announced June 2013.

Comments: 20 pages, 11 figures. (Submitted to MNRAS on 21st March 2013)

The Kinematics and Ionization of Nuclear Gas Clouds in Centaurus A

Authors: Geoffrey V. Bicknell, Ralph S. Sutherland, Nadine Neumayer

Abstract: Neumayer et al. established the existence of a blue-shifted cloud in the core of Centaurus A, within a few parsecs of the nucleus and close to the radio jet. We propose that the cloud has been impacted by the jet, and that it is in the foreground of the jet, accounting for its blue-shifted emission on the Southern side of the nucleus. We consider both shock excitation and photoionization models fo… ▽ More Neumayer et al. established the existence of a blue-shifted cloud in the core of Centaurus A, within a few parsecs of the nucleus and close to the radio jet. We propose that the cloud has been impacted by the jet, and that it is in the foreground of the jet, accounting for its blue-shifted emission on the Southern side of the nucleus. We consider both shock excitation and photoionization models for the excitation of the cloud. Shock models do not account for the [SiVI] and [CaVIII] emission line fluxes. However, X-ray observations indicate a source of ionizing photons in the core of Centaurus A; photoionization by the inferred flux incident on the cloud can account for the fluxes in these lines relative to Brackett-gamma. The power-law slope of the ionizing continuum matches that inferred from synchrotron models of the X-rays. The logarithm of the ionization parameter is -1.9, typical of that in Seyfert galaxies and consistent with the value proposed for dusty ionized plasmas. The model cloud density depends upon the Lorentz factor of the blazar and the inclination of our line of sight to the jet axis. For acute inclinations, the inferred density is consistent with expected cloud densities. However, for moderate inclinations of the jet to the line of sight, high Lorentz factors imply cloud densities in excess of $10^5 cm^{-3}$ and very low filling factors, suggesting that models of the gamma ray emission should incorporate jet Lorentz factors $\lesssim 5$. △ Less

Submitted 5 March, 2013; originally announced March 2013.

Comments: Published online by the Astrophysical Journal

Journal ref: ApJ, 2013, 766, 36 http://stacks.iop.org/0004-637X/766/36

Equilibrium disks, MRI mode excitation, and steady state turbulence in global accretion disk simulations

Authors: E. R. Parkin, G. V. Bicknell

Abstract: Global three dimensional magnetohydrodynamic (MHD) simulations of turbulent accretion disks are presented which start from fully equilibrium initial conditions in which the magnetic forces are accounted for and the induction equation is satisfied. The local linear theory of the magnetorotational instability (MRI) is used as a predictor of the growth of magnetic field perturbations in the global si… ▽ More Global three dimensional magnetohydrodynamic (MHD) simulations of turbulent accretion disks are presented which start from fully equilibrium initial conditions in which the magnetic forces are accounted for and the induction equation is satisfied. The local linear theory of the magnetorotational instability (MRI) is used as a predictor of the growth of magnetic field perturbations in the global simulations. The linear growth estimates and global simulations diverge when non-linear motions - perhaps triggered by the onset of turbulence - upset the velocity perturbations used to excite the MRI. The saturated state is found to be independent of the initially excited MRI mode, showing that once the disk has expelled the initially net flux field and settled into quasi-periodic oscillations in the toroidal magnetic flux, the dynamo cycle regulates the global saturation stress level. Furthermore, time-averaged measures of converged turbulence, such as the ratio of magnetic energies, are found to be in agreement with previous works. In particular, the globally averaged stress normalized to the gas pressure, <α_{\rm P}> = 0.034, with notably higher values achieved for simulations with higher azimuthal resolution. Supplementary tests are performed using different numerical algorithms and resolutions. Convergence with resolution during the initial linear MRI growth phase is found for 23-35 cells per scaleheight (in the vertical direction). △ Less

Submitted 6 December, 2012; originally announced December 2012.

Comments: 19 pages, 13 figures, accepted for ApJ

Ultra Fast Outflows: Galaxy-Scale Active Galactic Nucleus Feedback

Authors: A. Y. Wagner, M. Umemura, G. V. Bicknell

Abstract: We show, using global 3D grid-based hydrodynamical simulations, that Ultra Fast Outflows (UFOs) from Active Galactic Nuclei (AGN) result in considerable feedback of energy and momentum into the interstellar medium (ISM) of the host galaxy. The AGN wind interacts strongly with the inhomogeneous, two-phase ISM consisting of dense clouds embedded in a tenuous hot hydrostatic medium. The outflow flood… ▽ More We show, using global 3D grid-based hydrodynamical simulations, that Ultra Fast Outflows (UFOs) from Active Galactic Nuclei (AGN) result in considerable feedback of energy and momentum into the interstellar medium (ISM) of the host galaxy. The AGN wind interacts strongly with the inhomogeneous, two-phase ISM consisting of dense clouds embedded in a tenuous hot hydrostatic medium. The outflow floods through the inter-cloud channels, sweeps up the hot ISM, and ablates and disperses the dense clouds. The momentum of the UFO is primarily transferred to the dense clouds via the ram pressure in the channel flow, and the wind-blown bubble evolves in the energy-driven regime. Any dependence on UFO opening angle disappears after the first interaction with obstructing clouds. On kpc scales, therefore, feedback by UFOs operates similarly to feedback by relativistic AGN jets. Negative feedback is significantly stronger if clouds are distributed spherically, rather than in a disc. In the latter case the turbulent backflow of the wind drives mass inflow toward the central black hole. Considering the common occurrence of UFOs in AGN, they are likely to be important in the cosmological feedback cycles of galaxy formation. △ Less

Submitted 1 January, 2013; v1 submitted 25 November, 2012; originally announced November 2012.

Comments: 5 pages, 4 figures

Periodic structure in the Mpc-scale jet of PKS 0637-752

Authors: L. E. H. Godfrey, J. E. J. Lovell, S. Burke-Spolaor, R. Ekers, G. V. Bicknell, M. Birkinshaw, D. M. Worrall, D. L. Jauncey, D. A. Schwartz, H. L. Marshall, J. Gelbord, E. S. Perlman

Abstract: We present 18 GHz Australia Telescope Compact Array imaging of the Mpc-scale quasar jet PKS 0637-752 with angular resolution ~0.58 arcseconds. We draw attention to a spectacular train of quasi-periodic knots along the inner 11 arcseconds of the jet, with average separation ~1.1 arcsec (7.6 kpc projected). We consider two classes of model to explain the periodic knots: those that involve a static p… ▽ More We present 18 GHz Australia Telescope Compact Array imaging of the Mpc-scale quasar jet PKS 0637-752 with angular resolution ~0.58 arcseconds. We draw attention to a spectacular train of quasi-periodic knots along the inner 11 arcseconds of the jet, with average separation ~1.1 arcsec (7.6 kpc projected). We consider two classes of model to explain the periodic knots: those that involve a static pattern through which the jet plasma travels (e.g. stationary shocks); and those that involve modulation of the jet engine. Interpreting the knots as re-confinement shocks implies the jet kinetic power Q ~ 10^{46} erg/s, but the constant knot separation along the jet is not expected in a realistic external density profile. For models involving modulation of the jet engine, we find that the required modulation period is 2 x 10^3 yr < τ< 3 x 10^5 yr. The lower end of this range is applicable if the jet remains highly relativistic on kpc-scales, as implied by the IC/CMB model of jet X-ray emission. We suggest that the quasi-periodic jet structure in PKS 0637-752 may be analogous to the quasi-periodic jet modulation seen in the microquasar GRS 1915+105, believed to result from limit cycle behaviour in an unstable accretion disk. If variations in the accretion rate are driven by a binary black hole, the predicted orbital radius is 0.7 < a < 30 pc, which corresponds to a maximum angular separation of ~0.1 - 5 mas. △ Less

Submitted 20 September, 2012; originally announced September 2012.

Comments: 5 Pages, 2 Figures

A Multi-Wavelength Study of the Jet, Lobes and Core of the Quasar PKS 2101-490

Authors: L. E. H. Godfrey, G. V. Bicknell, J. E. J. Lovell, D. L. Jauncey, J. Gelbord, D. A. Schwartz, E. S. Perlman, H. L. Marshall, M. Birkinshaw, D. M. Worrall, M. Georganopoulos, D. W. Murphy

Abstract: We present a detailed study of the X-ray, optical and radio emission from the jet, lobes and core of the quasar PKS 2101-490 as revealed by new Chandra, HST and ATCA images. We extract the radio to X-ray spectral energy distributions from seven regions of the 13 arcsecond jet, and model the jet X-ray emission in terms of Doppler beamed inverse Compton scattering of the cosmic microwave background… ▽ More We present a detailed study of the X-ray, optical and radio emission from the jet, lobes and core of the quasar PKS 2101-490 as revealed by new Chandra, HST and ATCA images. We extract the radio to X-ray spectral energy distributions from seven regions of the 13 arcsecond jet, and model the jet X-ray emission in terms of Doppler beamed inverse Compton scattering of the cosmic microwave background (IC/CMB) for a jet in a state of equipartition between particle and magnetic field energy densities. This model implies that the jet remains highly relativistic hundreds of kpc from the nucleus, with a bulk Lorentz factor Gamma ~ 6 and magnetic field of order 30 microGauss. We detect an apparent radiative cooling break in the synchrotron spectrum of one of the jet knots, and are able to interpret this in terms of a standard one-zone continuous injection model, based on jet parameters derived from the IC/CMB model. However, we note apparent substructure in the bright optical knot in one of the HST bands. We confront the IC/CMB model with independent estimates of the jet power, and find that the IC/CMB model jet power is consistent with the independent estimates, provided that the minimum electron Lorentz factor gamma_min > 50, and the knots are significantly longer than the jet width, as implied by de-projection of the observed knot lengths. △ Less

Submitted 31 July, 2012; originally announced August 2012.

Comments: 16 pages, 10 figures, 6 tables

Driving Outflows with Relativistic Jets and the Dependence of AGN Feedback Efficiency on ISM Inhomogeneity

Authors: A. Y. Wagner, G. V. Bicknell, M. Umemura

Abstract: We examine the detailed physics of the feedback mechanism by relativistic AGN jets interacting with a two-phase fractal interstellar medium in the kpc-scale core of galaxies using 29 3D grid-based hydrodynamical simulations. The feedback efficiency, as measured by the amount of cloud-dispersal generated by the jet-ISM interactions, is sensitive to the maximum size of clouds in the fractal cloud di… ▽ More We examine the detailed physics of the feedback mechanism by relativistic AGN jets interacting with a two-phase fractal interstellar medium in the kpc-scale core of galaxies using 29 3D grid-based hydrodynamical simulations. The feedback efficiency, as measured by the amount of cloud-dispersal generated by the jet-ISM interactions, is sensitive to the maximum size of clouds in the fractal cloud distribution but not to their volume filling factor. Feedback ceases to be efficient for Eddington ratios P_jet/L_edd<10^-4, although systems with large cloud complexes ~50 pc require jets of Eddington ratio in excess of 10^-2 to disperse the clouds appreciably. Based on measurements of the bubble expansion rates in our simulations we argue that sub-grid AGN prescriptions resulting in negative feedback in cosmological simulations without a multi-phase treatment of the ISM are good approximations if the volume filling factor of warm phase material is less than 0.1 and the cloud complexes are smaller than ~25 pc. We find that the acceleration of the dense embedded clouds is provided by the ram pressure of the high velocity flow through the porous channels of the warm phase, flow that has fully entrained the shocked hot-phase gas it has swept up, and is additionally mass-loaded by ablated cloud material. This mechanism transfers 10% to 40% of the jet energy to the cold and warm gas, accelerating it within 10 to 100 Myr to velocities that match those observed in a range of high and low redshift radio galaxies hosting powerful radio jets. △ Less

Submitted 9 August, 2012; v1 submitted 1 May, 2012; originally announced May 2012.

Comments: 26 pages, 12 figures. Accepted by the Astrophysical Journal on August 9, 2012

Searches for very high energy gamma rays from blazars with CANGAROO-III telescope in 2005-2009

Authors: Y. Mizumura, J. Kushida, K. Nishijima, G. V. Bicknell, R. W. Clay, P. G. Edwards, S. Gunji, S. Hara, S. Hayashi, S. Kabuki, F. Kajino, A. Kawachi, T. Kifune, R. Kiuchi, K. Kodani, Y. Matsubara, T. Mizukami, Y. Mizumoto, M. Mori, H. Muraishi, T. Naito, M. Ohishi, V. Stamatescu, D. L. Swaby, T. Tanimori , et al. (6 additional authors not shown)

Abstract: We have searched for very high energy (VHE) gamma rays from four blazars using the CANGAROO-III imaging atmospheric Cherenkov telescope. We report the results of the observations of H 2356-309, PKS 2155-304, PKS 0537-441, and 3C 279, performed from 2005 to 2009, applying a new analysis to suppress the effects of the position dependence of Cherenkov images in the field of view. No significant VHE g… ▽ More We have searched for very high energy (VHE) gamma rays from four blazars using the CANGAROO-III imaging atmospheric Cherenkov telescope. We report the results of the observations of H 2356-309, PKS 2155-304, PKS 0537-441, and 3C 279, performed from 2005 to 2009, applying a new analysis to suppress the effects of the position dependence of Cherenkov images in the field of view. No significant VHE gamma ray emission was detected from any of the four blazars. The GeV gamma-ray spectra of these objects were obtained by analyzing Fermi/LAT archival data. Non-simultaneous wide range (radio to VHE gamma-ray bands) spectral energy distributions (SEDs) including CANGAROO-III upper limits, GeV gamma-ray spectra, and archival data are discussed using a one-zone synchrotron self-Compton (SSC) model in combination with a external Compton (EC) radiation. The HBLs (H 2356-309 and PKS 2155-304) can be explained by a simple SSC model, and PKS 0537-441 and 3C 279 are well modeled by a combination of SSC and EC model. We find a consistency with the blazar sequence in terms of strength of magnetic field and component size. △ Less

Submitted 19 February, 2012; v1 submitted 9 February, 2012; originally announced February 2012.

Comments: 11 pages, 8 figures, Accepted for publication in Astroparticle Physics

Journal ref: Astroparticle Physics 35 (2012), pp. 563-572

CANGAROO-III observation of TeV gamma rays from the unidentified gamma-ray source HESS J1614-518

Authors: T. Mizukami, H. Kubo, T. Yoshida, T. Nakamori, R. Enomoto, T. Tanimori, M. Akimoto, G. V. Bicknell, R. W. Clay, P. G. Edwards, S. Gunji, S. Hara, T. Hara, S. Hayashi, H. Ishioka, S. Kabuki, F. Kajino, H. Katagiri, A. Kawachi, T. Kifune, R. Kiuchi, T. Kunisawa, J. Kushida, T. Matoba, Y. Matsubara , et al. (24 additional authors not shown)

Abstract: We report the detection, with the CANGAROO-III imaging atmospheric Cherenkov telescope array, of a very high energy gamma-ray signal from the unidentified gamma-ray source HESS J1614-518, which was discovered in the H.E.S.S. Galactic plane survey. Diffuse gamma-ray emission was detected above 760 GeV at the 8.9 sigma level during an effective exposure of 54 hr from 2008 May to August. The spectrum… ▽ More We report the detection, with the CANGAROO-III imaging atmospheric Cherenkov telescope array, of a very high energy gamma-ray signal from the unidentified gamma-ray source HESS J1614-518, which was discovered in the H.E.S.S. Galactic plane survey. Diffuse gamma-ray emission was detected above 760 GeV at the 8.9 sigma level during an effective exposure of 54 hr from 2008 May to August. The spectrum can be represented by a power-law: 8.2+-2.2_{stat}+-2.5_{sys}x10^{-12}x (E/1TeV)^{-Gamma} cm^{-2} s^{-1} TeV^{-1} with a photon index Gamma of 2.4+-0.3_{stat}+-0.2_{sys}, which is compatible with that of the H.E.S.S. observations. By combining our result with multi-wavelength data, we discuss the possible counterparts for HESS J1614-518 and consider radiation mechanisms based on hadronic and leptonic processes for a supernova remnant, stellar winds from massive stars, and a pulsar wind nebula. Although a leptonic origin from a pulsar wind nebula driven by an unknown pulsar remains possible, hadronic-origin emission from an unknown supernova remnant is preferred. △ Less

Submitted 21 July, 2011; v1 submitted 20 July, 2011; originally announced July 2011.

Comments: 9 pages, 7 figures, accepted for publication in ApJ

Relativistic Jet Feedback in Evolving Galaxies

Authors: Alexander Y. Wagner, Geoffrey V. Bicknell

Abstract: Over cosmic time, galaxies grow through the hierarchical merging of smaller galaxies. However, the bright region of the galaxy luminosity function is incompatible with the simplest version of hierarchical merging, and it is believed that feedback from the central black hole in the host galaxies reduces the number of bright galaxies and regulates the co-evolution of black hole and host galaxy. Nume… ▽ More Over cosmic time, galaxies grow through the hierarchical merging of smaller galaxies. However, the bright region of the galaxy luminosity function is incompatible with the simplest version of hierarchical merging, and it is believed that feedback from the central black hole in the host galaxies reduces the number of bright galaxies and regulates the co-evolution of black hole and host galaxy. Numerous simulations of galaxy evolution have attempted to include the physical effects of such feedback with a resolution usually exceeding a kiloparsec. However, interactions between jets and the interstellar medium involve processes occurring on less than kiloparsec scales. In order to further the understanding of processes occurring on such scales, we present a suite of simulations of relativistic jets interacting with a fractal two-phase interstellar medium with a resolution of two parsecs and a largest scale of one kiloparsec. The transfer of energy and momentum to the interstellar medium is considerable, and we find that jets with powers in the range of 10^43-10^46 erg s^-1 can inhibit star formation through the dispersal of dense gas in the galaxy core. We determine the effectiveness of this process as a function of the ratio of the jet power to the Eddington luminosity of the black hole, the pressure of the interstellar medium and the porosity of the dense gas. △ Less

Submitted 6 December, 2010; originally announced December 2010.

Comments: 11 pages, 6 figures, 1 table

Optical IFU Observations of the Brightest Cluster Galaxy NGC 4696: The Case for a Minor Merger and Shock-excited Filaments

Authors: C. L. Farage, P. J. McGregor, M. A. Dopita, G. V. Bicknell

Abstract: We present deep optical integral-field spectroscopic observations of the nearby (z ~ 0.01) brightest cluster galaxy NGC 4696 in the core of the Centaurus Cluster, made with the Wide Field Spectrograph (WiFeS) on the ANU 2.3m telescope at Siding Spring Observatory. We investigate the morphology, kinematics, and excitation of the emission-line filaments and discuss these in the context of a model of… ▽ More We present deep optical integral-field spectroscopic observations of the nearby (z ~ 0.01) brightest cluster galaxy NGC 4696 in the core of the Centaurus Cluster, made with the Wide Field Spectrograph (WiFeS) on the ANU 2.3m telescope at Siding Spring Observatory. We investigate the morphology, kinematics, and excitation of the emission-line filaments and discuss these in the context of a model of a minor merger. We suggest that the emission-line filaments in this object have their origin in the accretion of a gas-rich galaxy and that they are excited by v ~100-200 km/s shocks driven into the cool filament gas by the ram pressure of the transonic passage of the merging system through the hot halo gas of NGC 4696. △ Less

Submitted 15 September, 2010; originally announced September 2010.

Comments: 19 pages, 15 figures, accepted for publication in The Astrophysical Journal

Search for VHE gamma rays from SS433/W50 with the CANGAROO-II telescope

Authors: Sei. Hayashi, F. Kajino, T. Naito, A. Asahara, G. V. Bicknell, R. W. Clay, Y. Doi, P. G. Edwards, R. Enomoto, S. Gunji, S. Hara, T. Hara, T. Hattori, C. Itoh, S. Kabuki, H. Katagiri, A. Kawachi, T. Kifune, L. T. Ksenofontov, H. Kubo, T. Kurihara, R. Kurosaka, J. Kushida, Y. Matsubara, Y. Miyashita , et al. (26 additional authors not shown)

Abstract: SS433, located at the center of the supernova remnant W50, is a close proximity binary system consisting of a compact star and a normal star. Jets of material are directed outwards from the vicinity of the compact star symmetrically to the east and west. Non-thermal hard X-ray emission is detected from lobes lying on both sides. Shock accelerated electrons are expected to generate sub-TeV gamma… ▽ More SS433, located at the center of the supernova remnant W50, is a close proximity binary system consisting of a compact star and a normal star. Jets of material are directed outwards from the vicinity of the compact star symmetrically to the east and west. Non-thermal hard X-ray emission is detected from lobes lying on both sides. Shock accelerated electrons are expected to generate sub-TeV gamma rays through the inverse-Compton process in the lobes. Observations of the western X-ray lobe region of SS433/W50 system have been performed to detect sub-TeV gamma-rays using the 10m CANGAROO-II telescope in August and September, 2001, and July and September, 2002. The total observation times are 85.2 hours for ON source, and 80.8 hours for OFF source data. No significant excess of sub-TeV gamma rays has been found at 3 regions of the western X-ray lobe of SS433/W50 system. We have derived 99% confidence level upper limits to the fluxes of gamma rays and have set constraints on the strengths of the magnetic fields assuming the synchrotron/inverse-Compton model for the wide energy range of photon spectrum from radio to TeV. The derived lower limits are 4.3 microgauss for the center of the brightest X-ray emission region and 6.3 microgauss for the far end from SS433 in the western X-ray lobe. In addition, we suggest that the spot-like X-ray emission may provide a major contribution to the hardest X-ray spectrum in the lobe. △ Less

Submitted 1 September, 2009; originally announced September 2009.

Comments: 7 pages, 8 figures, to be published in Astroparticle Physics

Journal ref: Astropart.Phys.32:112-119,2009

CANGAROO-III search for TeV Gamma-rays from two clusters of galaxies

Authors: R. Kiuchi, M. Mori, G. V. Bicknell, R. W. Clay, P. G. Edwards, R. Enomoto, S. Gunji, S. Hara, T. Hara, T. Hattori, S. Hayashi, Y. Higashi, Y. Hirai, K. Inoue, C. Itoh, S. Kabuki, F. Kajino, H. Katagiri, A. Kawachi, T. Kifune, H. Kubo, J. Kushida, Y. Matsubara, T. Mizukami, Y. Mizumoto , et al. (25 additional authors not shown)

Abstract: Because accretion and merger shocks in clusters of galaxies may accelerate particles to high energies, clusters are candidate sites for the origin of ultra-high-energy (UHE) cosmic-rays. A prediction was presented for gamma-ray emission from a cluster of galaxies at a detectable level with the current generation of imaging atmospheric Cherenkov telescopes. The gamma-ray emission was produced via… ▽ More Because accretion and merger shocks in clusters of galaxies may accelerate particles to high energies, clusters are candidate sites for the origin of ultra-high-energy (UHE) cosmic-rays. A prediction was presented for gamma-ray emission from a cluster of galaxies at a detectable level with the current generation of imaging atmospheric Cherenkov telescopes. The gamma-ray emission was produced via inverse Compton upscattering of cosmic microwave background (CMB) photons by electron-positron pairs generated by collisions of UHE cosmic rays in the cluster. We observed two clusters of galaxies, Abell 3667 and Abell 4038, searching for very-high-energy gamma-ray emission with the CANGAROO-III atmospheric Cherenkov telescope system in 2006. The analysis showed no significant excess around these clusters, yielding upper limits on the gamma-ray emission. From a comparison of the upper limit for the north-west radio relic region of Abell 3667 with a model prediction, we derive a lower limit for the magnetic field of the region of ~0.1 micro G. This shows the potential of gamma-ray observations in studies of the cluster environment. We also discuss the flux upper limit from cluster center regions using a model of gamma-ray emission from neutral pions produced in hadronic collisions of cosmic-ray protons with the intra-cluster medium (ICM). The derived upper limit of the cosmic-ray energy density within this framework is an order of magnitude higher than that of our Galaxy. △ Less

Submitted 23 August, 2009; originally announced August 2009.

Comments: 7 pages, 6 figures, Accepted in ApJ

Accretion Discs in Blazars

Authors: E. J. D. Jolley, Z. Kuncic, G. V. Bicknell, S. Wagner

Abstract: The characteristic properties of blazars (rapid variability, strong polarization, high brightness) are widely attributed to a powerful relativistic jet oriented close to our line of sight. Despite the spectral energy distributions (SEDs) being strongly jet-dominated, a "big blue bump" has been recently detected in sources known as flat spectrum radio quasars (FSRQs). These new data provide a uni… ▽ More The characteristic properties of blazars (rapid variability, strong polarization, high brightness) are widely attributed to a powerful relativistic jet oriented close to our line of sight. Despite the spectral energy distributions (SEDs) being strongly jet-dominated, a "big blue bump" has been recently detected in sources known as flat spectrum radio quasars (FSRQs). These new data provide a unique opportunity to observationally test coupled jet-disc accretion models in these extreme sources. In particular, as energy and angular momentum can be extracted by a jet magnetically coupled to the accretion disc, the thermal disc emission spectrum may be modified from that predicted by the standard model for disc accretion. We compare the theoretically predicted jet-modified accretion disc spectra against the new observations of the "big blue bump" in FSRQs. We find mass accretion rates that are higher, typically by a factor of two, than predicted by standard accretion disc theory. Furthermore, our results predict that the high redshift blazars PKS 0836+710, PKS 2149-307, B2 0743+25 and PKS 0537-286 may be predominantly powered by a low or moderate spin (a < 0.6) black hole with high mass accretion rates mdot_a ~ 50 - 200 msol/yr, while 3C 273 harbours a rapidly spinning black hole (a = 0.97) with mdot_a ~ 20 msol/yr. We also find that the black hole masses in these high redshift sources must be > 5 * 10^9 msol. △ Less

Submitted 19 August, 2009; v1 submitted 17 August, 2009; originally announced August 2009.

Comments: Accepted for publication (17 August 2009) in MNRAS

CANGAROO-III Observation of TeV Gamma Rays from the vicinity of PSR B1 706-44

Authors: R. Enomoto, J. Kushida, T. Nakamori, T. Kifune, G. V. Bicknell, R. W. Clay, P. G. Edwards, S. Gunji, S. Hara, T. Hara, T. Hattori, S. Hayashi, Y. Higashi, Y. Hirai, K. Inoue, H. Ishioka, S. Kabuki, F. Kajino, H. Katagiri, A. Kawachi, R. Kiuchi, H. Kubo, T. Kunisawa, T. Matoba, Y. Matsubara , et al. (28 additional authors not shown)

Abstract: Observation by the CANGAROO-III stereoscopic system of the Imaging Cherenkov Telescope has detected extended emission of TeV gamma rays in the vicinity of the pulsar PSR B1706$-$44. The strength of the signal observed as gamma-ray-like events varies when we apply different ways of emulating background events. The reason for such uncertainties is argued in relevance to gamma-rays embedded in the… ▽ More Observation by the CANGAROO-III stereoscopic system of the Imaging Cherenkov Telescope has detected extended emission of TeV gamma rays in the vicinity of the pulsar PSR B1706$-$44. The strength of the signal observed as gamma-ray-like events varies when we apply different ways of emulating background events. The reason for such uncertainties is argued in relevance to gamma-rays embedded in the "off-source data", that is, unknown sources and diffuse emission in the Galactic plane, namely, the existence of a complex structure of TeV gamma-ray emission around PSR B1706$-$44. △ Less

Submitted 12 August, 2009; originally announced August 2009.

Comments: 10 pages, 13 figures, to be published in ApJ

Journal ref: Astrophysical Journal 703 (2009) 1725-1733

Starburst-Driven Galactic Winds: Filament Formation and Emission Processes

Authors: J. L. Cooper, G. V. Bicknell, R. S. Sutherland, J. Bland-Hawthorn

Abstract: We have performed a series of three-dimensional simulations of the interaction of a supersonic wind with a non-spherical radiative cloud. These simulations are motivated by our recent three-dimensional model of a starburst-driven galactic wind interacting with an inhomogeneous disk, which show that an optically emitting filament can be formed by the break-up and acceleration of a cloud into a su… ▽ More We have performed a series of three-dimensional simulations of the interaction of a supersonic wind with a non-spherical radiative cloud. These simulations are motivated by our recent three-dimensional model of a starburst-driven galactic wind interacting with an inhomogeneous disk, which show that an optically emitting filament can be formed by the break-up and acceleration of a cloud into a supersonic wind. In this study we consider the evolution of a cloud with two different geometries (fractal and spherical) and investigate the importance of radiative cooling on the cloud's survival. We have also undertaken a comprehensive resolution study in order to ascertain the effect of the assumed numerical resolution on the results. We find that the ability of the cloud to radiate heat is crucial for its survival. While an adiabatic cloud is destroyed over a short period of time, a radiative cloud is broken up via the Kelvin-Helmholtz instability into numerous small, dense cloudlets, which are drawn into the flow to form a filamentary structure. The degree of fragmentation is highly dependent on the resolution of the simulation, with the number of cloudlets formed increasing as the Kelvin-Helmholtz instability is better resolved. Nevertheless, there is a clear qualitative trend, with the filamentary structure still persistent at high resolution. We confirm the mechanism behind the formation of the H-alpha emitting filaments found in our global simulations of a starburst-driven wind. Based on our resolution study, we conclude that bow shocks around accelerated gas clouds, and their interaction, are the main source of the soft X-ray emission observed in these galactic-scale winds. [ABRIDGED] △ Less

Submitted 23 July, 2009; originally announced July 2009.

Comments: Accepted to ApJ, 39 pages, 21 figures, movie file can obtained at http://www.mso.anu.edu.au/~jcooper/movie/halpha.mov

Enhanced MHD transport in astrophysical accretion flows: turbulence, winds and jets

Authors: Peter B. Dobbie, Zdenka Kuncic, Geoffrey V. Bicknell, Raquel Salmeron

Abstract: Astrophysical accretion is arguably the most prevalent physical process in the Universe; it occurs during the birth and death of individual stars and plays a pivotal role in the evolution of entire galaxies. Accretion onto a black hole, in particular, is also the most efficient mechanism known in nature, converting up to 40% of accreting rest mass energy into spectacular forms such as high-energ… ▽ More Astrophysical accretion is arguably the most prevalent physical process in the Universe; it occurs during the birth and death of individual stars and plays a pivotal role in the evolution of entire galaxies. Accretion onto a black hole, in particular, is also the most efficient mechanism known in nature, converting up to 40% of accreting rest mass energy into spectacular forms such as high-energy (X-ray and gamma-ray) emission and relativistic jets. Whilst magnetic fields are thought to be ultimately responsible for these phenomena, our understanding of the microphysics of MHD turbulence in accretion flows as well as large-scale MHD outflows remains far from complete. We present a new theoretical model for astrophysical disk accretion which considers enhanced vertical transport of momentum and energy by MHD winds and jets, as well as transport resulting from MHD turbulence. We also describe new global, 3D simulations that we are currently developing to investigate the extent to which non-ideal MHD effects may explain how small-scale, turbulent fields (generated by the magnetorotational instability -- MRI) might evolve into large-scale, ordered fields that produce a magnetized corona and/or jets where the highest energy phenomena necessarily originate. △ Less

Submitted 8 December, 2009; v1 submitted 19 March, 2009; originally announced March 2009.

Comments: 8 pages, 2 figures. Minor revision, published version: Proc 14th International Congress on Plasma Physics, Fukuoka, Japan, Sep 2008

Journal ref: Plas. Fus. Res., 4, 017 (2009)

New Global 3D MHD Simulations of Black Hole Disk Accretion and Outflows

Authors: Peter B. Dobbie, Zdenka Kuncic, Geoffrey V. Bicknell, Raquel Salmeron

Abstract: It is widely accepted that quasars and other active galactic nuclei (AGN) are powered by accretion of matter onto a central supermassive black hole. While numerical simulations have demonstrated the importance of magnetic fields in generating the turbulence believed necessary for accretion, so far they have not produced the high mass accretion rates required to explain the most powerful sources.… ▽ More It is widely accepted that quasars and other active galactic nuclei (AGN) are powered by accretion of matter onto a central supermassive black hole. While numerical simulations have demonstrated the importance of magnetic fields in generating the turbulence believed necessary for accretion, so far they have not produced the high mass accretion rates required to explain the most powerful sources. We describe new global 3D simulations we are developing to assess the importance of radiation and non-ideal MHD in generating magnetized outflows that can enhance the overall rates of angular momentum transport and mass accretion. △ Less

Submitted 2 February, 2009; originally announced February 2009.

Comments: 2 pages, including 1 colour figure. To appear in proceedings of IAU Symposium 259: "Cosmic Magnetic Fields: From Planets, To Stars and Galaxies", Tenerife, Nov 2008

A Multi-Wavelength Study of the High Surface Brightness Hotspot in PKS1421-490

Authors: L. E. H. Godfrey, G. V. Bicknell, J. E. J. Lovell, D. L. Jauncey, J. Gelbord, D. A. Schwartz, H. L. Marshall, M. Birkinshaw, M. Georganopoulos, D. W. Murphy, E. S. Perlman, D. M. Worrall

Abstract: Long Baseline Array imaging of the z=0.663 broad line radio galaxy PKS1421-490 reveals a 400 pc diameter high surface brightness hotspot at a projected distance of approximately 40kpc from the active galactic nucleus. The isotropic X-ray luminosity of the hotspot, L_{2-10 keV} = 3 10^{44} ergs/s, is comparable to the isotropic X-ray luminosity of the entire X-ray jet of PKS0637-752, and the peak… ▽ More Long Baseline Array imaging of the z=0.663 broad line radio galaxy PKS1421-490 reveals a 400 pc diameter high surface brightness hotspot at a projected distance of approximately 40kpc from the active galactic nucleus. The isotropic X-ray luminosity of the hotspot, L_{2-10 keV} = 3 10^{44} ergs/s, is comparable to the isotropic X-ray luminosity of the entire X-ray jet of PKS0637-752, and the peak radio surface brightness is hundreds of times greater than that of the brightest hotspot in Cygnus A. We model the radio to X-ray spectral energy distribution using a one-zone synchrotron self Compton model with a near equipartition magnetic field strength of 3 mG. There is a strong brightness asymmetry between the approaching and receding hotspots and the hot spot spectrum remains flat (alpha ~ 0.5) well beyond the predicted cooling break for a 3 mG magnetic field, indicating that the hotspot emission may be Doppler beamed. A high plasma velocity beyond the terminal jet shock could be the result of a dynamically important magnetic field in the jet. There is a change in the slope of the hotspot radio spectrum at GHz frequencies from alpha~0.5 to alpha<0.2, which we model by incorporating a cut-off in the electron energy distribution at gamma_{min} ~ 650, with higher values implied if the hotspot emission is Doppler beamed. We show that a sharp decrease in the electron number density below a Lorentz factor of 650 would arise from the dissipation of bulk kinetic energy in an electron/proton jet with a Lorentz factor Gamma_{jet} ~ 5. △ Less

Submitted 22 January, 2009; originally announced January 2009.

Comments: Accepted for publication in the Astrophysical Journal. 19 pages. 9 figures

Journal ref: Astrophys.J.695:707-723,2009

MRC B0319-454: Probing the large-scale structure with a giant radio galaxy

Authors: V. Safouris, R. Subrahmanyan, G. V. Bicknell, L. Saripalli

Abstract: We present an investigation of the relationships between the radio properties of a giant radio galaxy MRC B0319-454 and the surrounding galaxy distribution with the aim of examining the influence of intergalactic gas and gravity associated with the large-scale structure on the evolution in the radio morphology. Our new radio continuum observations of the radio source, with high surface brightnes… ▽ More We present an investigation of the relationships between the radio properties of a giant radio galaxy MRC B0319-454 and the surrounding galaxy distribution with the aim of examining the influence of intergalactic gas and gravity associated with the large-scale structure on the evolution in the radio morphology. Our new radio continuum observations of the radio source, with high surface brightness sensitivity, images the asymmetries in the megaparsec-scale radio structure in total intensity and polarization. We compare these with the 3-D galaxy distribution derived from galaxy redshift surveys. Galaxy density gradients are observed along and perpendicular to the radio axis: the large-scale structure is consistent with a model wherein the galaxies trace the ambient intergalactic gas and the evolution of the radio structures are ram-pressure limited by this associated gas. Additionally, we have modeled the off-axis evolution of the south-west radio lobe as deflection of a buoyant jet backflow by a transverse gravitational field: the model is plausible if entrainment is small. The case study presented here is a demonstration that giant radio galaxies may be useful probes of the warm-hot intergalactic medium believed to be associated with moderately over dense galaxy distributions. △ Less

Submitted 10 December, 2008; originally announced December 2008.

Comments: 27 pages, 15 figures, accepted for publication in MNRAS

Three-Dimensional Simulations of a Starburst-Driven Galactic Wind

Authors: Jackie L. Cooper, Geoffrey V. Bicknell, Ralph S. Sutherland, Joss Bland-Hawthorn

Abstract: We have performed a series of three-dimensional simulations of a starburst-driven wind in an inhomogeneous interstellar medium. The introduction of an inhomogeneous disk leads to differences in the formation of a wind, most noticeably the absence of the ``blow-out'' effect seen in homogeneous models. A wind forms from a series of small bubbles that propagate into the tenuous gas between dense cl… ▽ More We have performed a series of three-dimensional simulations of a starburst-driven wind in an inhomogeneous interstellar medium. The introduction of an inhomogeneous disk leads to differences in the formation of a wind, most noticeably the absence of the ``blow-out'' effect seen in homogeneous models. A wind forms from a series of small bubbles that propagate into the tenuous gas between dense clouds in the disk. These bubbles merge and follow the path of least resistance out of the disk, before flowing freely into the halo. Filaments are formed from disk gas that is broken up and accelerated into the outflow. These filaments are distributed throughout a biconical structure within a more spherically distributed hot wind. The distribution of the inhomogeneous interstellar medium in the disk is important in determining the morphology of this wind, as well as the distribution of the filaments. While higher resolution simulations are required in order to ascertain the importance of mixing processes, we find that soft X-ray emission arises from gas that has been mass-loaded from clouds in the disk, as well as from bow shocks upstream of clouds, driven into the flow by the ram pressure of the wind, and the interaction between these shocks. △ Less

Submitted 29 October, 2007; originally announced October 2007.

Comments: 37 pages, 16 figures, mpg movie can be obtained at http://www.mso.anu.edu.au/~jcooper/movie/video16.mpg, accepted for publication in ApJ

Interaction of Jets with the ISM of Radio Galaxies

Authors: R. S. Sutherland, G. V. Bicknell

Abstract: We present three dimensional simulations of the interaction of a light hypersonic jet with an inhomogeneous thermal and turbulently supported disk in an elliptical galaxy, including Radio and multi-band X-ray visualisations. These simulations are applicable to the GPS/CSS phase of some extragalactic radio sources. We identify four generic phases in the evolution of such a jet with the interste… ▽ More We present three dimensional simulations of the interaction of a light hypersonic jet with an inhomogeneous thermal and turbulently supported disk in an elliptical galaxy, including Radio and multi-band X-ray visualisations. These simulations are applicable to the GPS/CSS phase of some extragalactic radio sources. We identify four generic phases in the evolution of such a jet with the interstellar medium. The first is a `flood and channel'' phase, dominated by complex jet interactions with the dense cloudy medium close to the nucleus. A spherical, energy driven, bubble phase follows, where the bubble is larger than the disk scale, but the jet remains fully disrupted close to the nucleus. Then in a rapid, jet break--out phase, the jet breaks free of the last obstructing dense clouds, becomes collimated and pierces the more or less spherical bubble. In the final classical phase, the jet propagates in a momentum-dominated fashion similar to jets in single component hot haloes, leading to the classical jet -- cocoon -- bow-shock structure. △ Less

Submitted 25 July, 2007; originally announced July 2007.

Comments: 22 pages, 10 figures, PDFLaTeX, aastex macros, graphicx package. Proceedings of the 5th Stromlo Symposium: Disks, Winds, and Jets - from Planets to Quasars. Accepted, to be published in Astronomy and Space Physics

Journal ref: Astrophys.Space Sci.311:293-303,2007

Interactions of a Light Hypersonic Jet with a Non-Uniform Interstellar Medium

Authors: R. S. Sutherland, G. V. Bicknell

Abstract: We present three dimensional simulations of the interaction of a light hypersonic jet with an inhomogeneous thermal and turbulently supported disk in an elliptical galaxy. We model the jet as a light, supersonic non-relativistic flow with parameters selected to be consistent with a relativistic jet with kinetic power just above the FR1/FR2 break. We identify four generic phases in the evolutio… ▽ More We present three dimensional simulations of the interaction of a light hypersonic jet with an inhomogeneous thermal and turbulently supported disk in an elliptical galaxy. We model the jet as a light, supersonic non-relativistic flow with parameters selected to be consistent with a relativistic jet with kinetic power just above the FR1/FR2 break. We identify four generic phases in the evolution of such a jet with the inhomogeneous interstellar medium: 1) an initial ``flood and channel'' phase, where progress is characterized by high pressure gas finding changing weak points in the ISM, flowing through channels that form and re-form over time, 2) a spherical, energy-driven bubble phase, were the bubble is larger than the disk scale, but the jet remains fully disrupted close to the nucleus, 3) a rapid, jet break--out phase the where jet breaks free of the last dense clouds, becomes collimated and pierces the spherical bubble, and 4) a classical phase, the jet propagates in a momentum-dominated fashion leading to the classical jet + cocoon + bow-shock structure. Mass transport in the simulations is investigated, and we propose a model for the morphology and component proper motions in the well-studied Compact Symmetric Object 4C31.04. △ Less

Submitted 25 July, 2007; originally announced July 2007.

Comments: 66 pages, 22 figures, PDFLaTeX, aastex macros, graphicx and amssymb packages, Accepted, to be published 2007 ApJS

Towards a New Standard Theory for Astrophysical Disk Accretion

Authors: Zdenka Kuncic, Geoffrey V. Bicknell

Abstract: We briefly review recent developments in black hole accretion disk theory, placing new emphasis on the vital role played by magnetohydrodynamic (MHD) stresses in transporting angular momentum. The apparent universality of accretion-related outflow phenomena is a strong indicator that vertical transport of angular momentum by large-scale MHD torques is important and may even dominate radial trans… ▽ More We briefly review recent developments in black hole accretion disk theory, placing new emphasis on the vital role played by magnetohydrodynamic (MHD) stresses in transporting angular momentum. The apparent universality of accretion-related outflow phenomena is a strong indicator that vertical transport of angular momentum by large-scale MHD torques is important and may even dominate radial transport by small-scale MHD turbulence. This leads to an enhanced overall rate of angular momentum transport and allows accretion of matter to proceed at an interesting rate. Furthermore, we argue that when vertical transport is important, the radial structure of the accretion disk is modified and this affects the disk emission spectrum. We present a simple model demonstrating that energetic, magnetically-driven outflows give rise to a disk spectrum that is dimmer and redder than a standard accretion disk accreting at the same rate. We briefly discuss the implications of this key result for accreting black holes in different astrophysical systems. △ Less

Submitted 12 July, 2007; originally announced July 2007.

Comments: Accepted for publication as brief review in Mod. Phys. Let. A

Journal ref: Mod.Phys.Lett.A22:1685-1700,2007