Showing 1–50 of 122 results for author: Jones, T W

Hydrodynamic Simulations of a Relativistic Jet Interacting with the Intracluster Medium: Application to Cygnus A

Authors: John A. ZuHone, Paul E. J. Nulsen, Po-Hsun Tseng, Hsi-Yu Schive, Tom W. Jones

Abstract: The Fanaroff-Riley class II radio galaxy Cygnus A hosts jets which produce radio emission, X-ray cavities, cocoon shocks, and X-ray hotspots where the jet interacts with the ICM. Surrounding one hotspot is a peculiar "hole" feature which appears as a deficit in X-ray emission. We use relativistic hydrodynamic simulations of a collimated jet interacting with an inclined interface between lobe and c… ▽ More The Fanaroff-Riley class II radio galaxy Cygnus A hosts jets which produce radio emission, X-ray cavities, cocoon shocks, and X-ray hotspots where the jet interacts with the ICM. Surrounding one hotspot is a peculiar "hole" feature which appears as a deficit in X-ray emission. We use relativistic hydrodynamic simulations of a collimated jet interacting with an inclined interface between lobe and cluster plasma to model the basic processes which may lead to such a feature. We find that the jet reflects off of the interface into a broad, turbulent flow back out into the lobe, which is dominated by gas stripped from the interface at first and from the intracluster medium itself at later times. We produce simple models of X-ray emission from the ICM, the hotspot, and the reflected jet to show that a hole of emission surrounding the hotspot as seen in Cygnus A may be produced by Doppler de-boosting of the emission from the reflected jet as seen by an observer with a sight line nearly along the axis of the outgoing material. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: 16 pages, 11 figures, submitted to the journal Galaxies

Intracluster magnetic filaments and an encounter with a radio jet

Authors: Lawrence Rudnick, Marcus Bruggen, Gianfranco Brunetti, William Cotton, William Forman, Thomas W. Jones, Chris Nolting, Gerrit Schellenberger, Reinout van Weeren

Abstract: Thin synchrotron-emitting filaments are increasingly seen in the intracluster medium (ICM). We present the first example of a direct interaction between a magnetic filament, a radio jet, and a dense ICM clump in the poor cluster Abell 194. This enables the first exploration of the dynamics and possible histories of magnetic fields and cosmic rays in such filaments. Our observations are from the Me… ▽ More Thin synchrotron-emitting filaments are increasingly seen in the intracluster medium (ICM). We present the first example of a direct interaction between a magnetic filament, a radio jet, and a dense ICM clump in the poor cluster Abell 194. This enables the first exploration of the dynamics and possible histories of magnetic fields and cosmic rays in such filaments. Our observations are from the MeerKAT Galaxy Cluster Legacy Survey and the LOFAR Two Metre Sky Survey. Prominent 220 kpc long filaments extend east of radio galaxy 3C40B, with very faint extensions to 300 kpc, and show signs of interaction with its northern jet. They curve around a bend in the jet and intersect the jet in Faraday depth space. The X-ray surface brightness drops across the filaments; this suggests that the relativistic particles and fields contribute significantly to the pressure balance and evacuate the thermal plasma in a $\sim$35 kpc cylinder. We explore whether the relativistic electrons could have streamed along the filaments from 3C40B, and present a plausible alternative whereby magnetized filaments are a) generated by shear motions in the large-scale, post-merger ICM flow, b) stretched by interactions with the jet and flows in the ICM, amplifying the embedded magnetic fields, and c) perfused by re-energized relativistic electrons through betatron-type acceleration or diffusion of turbulently accelerated ICM cosmic ray electrons. We use the Faraday depth measurements to reconstruct some of the 3D structures of the filaments and of 3C40A and B. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: Accepted for publication in ApJ. 32 pages, 28 figures. Each of the last four figures has a corresponding movie file (mpeg), available under the Ancillary Files associated with this preprint

Plasma physics of the intracluster medium

Authors: Matthew W. Kunz, Thomas W. Jones, Irina Zhuravleva

Abstract: This Chapter provides a brief tutorial on some aspects of plasma physics that are fundamental to understanding the dynamics and energetics of the intracluster medium (ICM). The tutorial is split into two parts: one that focuses on the thermal plasma component -- its stability, viscosity, conductivity, and ability to amplify magnetic fields to dynamical strengths via turbulence and other plasma pro… ▽ More This Chapter provides a brief tutorial on some aspects of plasma physics that are fundamental to understanding the dynamics and energetics of the intracluster medium (ICM). The tutorial is split into two parts: one that focuses on the thermal plasma component -- its stability, viscosity, conductivity, and ability to amplify magnetic fields to dynamical strengths via turbulence and other plasma processes; and one that focuses on the non-thermal population of charged particles known as cosmic rays -- their acceleration, re-acceleration, and transport throughout the cluster volume. Observational context is woven throughout the narrative, from constraints on the strength and geometry of intracluster magnetic fields and the effective viscosity of the ICM, to examples of radio halos, radio relics, and cluster shocks that can test theories of particle acceleration. The promise of future X-ray missions to probe intracluster turbulence and discover the impact of small-scale plasma physics, coupled with sensitive, high-resolution radio observations of synchrotron-emitting plasma that reveal the properties of intracluster magnetic fields and particle-acceleration mechanisms, are likely to establish galaxy clusters as the premier cosmic laboratories for deciphering the fundamental physics of hot, dilute plasmas. △ Less

Submitted 5 May, 2022; originally announced May 2022.

Comments: This Chapter will appear in the Section "Galaxy Clusters" (Section Editors: E. Pointecouteau, E. Rasia, A. Simionescu) of the "Handbook of X-ray and Gamma-ray Astrophysics" (Editors in chief: C. Bambi and A. Santangelo)

Exploring the relation between turbulent velocity and density fluctuations in the stratified intracluster medium

Authors: Marco Simonte, Franco Vazza, Fabrizio Brighenti, Marcus Brueggen, Tom W. Jones

Abstract: The dynamics of the intracluster medium (ICM) is affected by turbulence driven by several processes, such as mergers, accretion and feedback from active galactic nuclei. X-ray surface brightness fluctuations have been used to constrain turbulence in galaxy clusters. Here, we use simulations to further investigate the relation between gas density and turbulent velocity fluctuations, with a focus on… ▽ More The dynamics of the intracluster medium (ICM) is affected by turbulence driven by several processes, such as mergers, accretion and feedback from active galactic nuclei. X-ray surface brightness fluctuations have been used to constrain turbulence in galaxy clusters. Here, we use simulations to further investigate the relation between gas density and turbulent velocity fluctuations, with a focus on the effect of the stratification of the ICM. In this work, we studied the turbulence driven by hierarchical accretion by analysing a sample of galaxy clusters simulated with the cosmological code ENZO. We used a fixed scale filtering approach to disentangle laminar from turbulent flows. In dynamically perturbed galaxy clusters, we found a relation between the root mean square of density and velocity fluctuations, albeit with a different slope than previously reported. The Richardson number is a parameter that represents the ratio between turbulence and buoyancy, and we found that this variable has a strong dependence on the filtering scale. However, we could not detect any strong relation between the Richardson number and the logarithmic density fluctuations, in contrast to results by recent and more idealised simulations. In particular, we find a strong effect from radial accretion, which appears to be the main driver for the gas fluctuations. The ubiquitous radial bias in the dynamics of the ICM suggests that homogeneity and isotropy are not always valid assumptions, even if the turbulent spectra follow Kolmogorov's scaling. Finally, we find that the slope of the velocity and density spectra are independent of cluster-centric radii. △ Less

Submitted 1 December, 2021; originally announced December 2021.

Comments: T14 pages, 18 figures

Journal ref: A&A 658, A149 (2022)

Properties of clumps and filaments around galaxy clusters

Authors: M. Angelinelli, S. Ettori, F. Vazza, T. W. Jones

Abstract: We report on the possibility of studying the proprieties of cosmic diffuse baryons by studying self-gravitating clumps and filaments connected to galaxy clusters. While filaments are challenging to detect with X-ray observations, the higher density of clumps makes them visible and a viable tracer to study the thermodynamical proprieties of baryons undergoing accretion along cosmic web filaments on… ▽ More We report on the possibility of studying the proprieties of cosmic diffuse baryons by studying self-gravitating clumps and filaments connected to galaxy clusters. While filaments are challenging to detect with X-ray observations, the higher density of clumps makes them visible and a viable tracer to study the thermodynamical proprieties of baryons undergoing accretion along cosmic web filaments onto galaxy clusters. We developed new algorithms to identify these structures and applied them to a set of non-radiative cosmological simulations of galaxy clusters at high resolution. We find that in those simulated clusters, the density and temperature of clumps are independent of the mass of the cluster where they reside. We detected a positive correlation between the filament temperature and the host cluster mass. The density and temperature of clumps and filaments also tended to correlate. Both the temperature and density decrease moving outward. We observed that clumps are hotter, more massive, and more luminous if identified closer to the cluster center. Especially in the outermost cluster regions (~3*R500,c or beyond), X-ray observations might already have the potential to locate cosmic filaments based on the distribution of clumps and to allow one to study the thermodynamics of diffuse baryons before they are processed by the intracluster medium. △ Less

Submitted 8 September, 2021; v1 submitted 1 February, 2021; originally announced February 2021.

Comments: Accepted to A&A; 19 pages, 25 figures

Journal ref: A&A 653, A171 (2021)

Chandra and XMM-Newton observations of A2256: cold fronts, merger shocks, and constraint on the IC emission

Authors: Chong Ge, Ruo-Yu Liu, Ming Sun, Heng Yu, Lawrence Rudnick, Jean Eilek, Frazer Owen, Sarthak Dasadia, Mariachiara Rossetti, Maxim Markevitch, Tracy E. Clarke, Thomas W. Jones, Simona Ghizzardi, Tiziana Venturi, Alexis Finoguenov, Dominique Eckert

Abstract: We present the results of deep Chandra and XMM-Newton observations of a complex merging galaxy cluster Abell 2256 (A2256) that hosts a spectacular radio relic (RR). The temperature and metallicity maps show clear evidence of a merger between the western subcluster (SC) and the primary cluster (PC). We detect five X-ray surface brightness edges. Three of them near the cluster center are cold fronts… ▽ More We present the results of deep Chandra and XMM-Newton observations of a complex merging galaxy cluster Abell 2256 (A2256) that hosts a spectacular radio relic (RR). The temperature and metallicity maps show clear evidence of a merger between the western subcluster (SC) and the primary cluster (PC). We detect five X-ray surface brightness edges. Three of them near the cluster center are cold fronts (CFs): CF1 is associated with the infalling SC; CF2 is located in the east of the PC; and CF3 is to the west of the PC core. The other two edges at cluster outskirts are shock fronts (SFs): SF1 near the RR in the NW has Mach numbers derived from the temperature and the density jumps, respectively, of $M_T=1.62\pm0.12$ and $M_ρろー=1.23\pm0.06$; SF2 in the SE has $M_T=1.54\pm0.05$ and $M_ρろー=1.16\pm0.13$. In the region of the RR, there is no evidence for the correlation between X-ray and radio substructures, from which we estimate an upper limit for the inverse-Compton emission, and therefore set a lower limit on the magnetic field ($\sim$ 450 kpc from PC center) of $B>1.0\ μみゅー$G for a single power-law electron spectrum or $B>0.4\ μみゅー$G for a broken power-law electron spectrum. We propose a merger scenario including a PC, an SC, and a group. Our merger scenario accounts for the X-ray edges, diffuse radio features, and galaxy kinematics, as well as projection effects. △ Less

Submitted 20 August, 2020; v1 submitted 21 May, 2020; originally announced May 2020.

Comments: 15 pages, 12 figures, MNRAS in press

Shocked Narrow-Angle Tail Radio Galaxies: Simulations and Emissions

Authors: Brian J. O'Neill, T. W. Jones, Chris Nolting, P. J. Mendygral

Abstract: We present a numerical study of the interactions between the elongated AGN outflows representing an evolved, narrow-angle tail (NAT) radio galaxy and planar, transverse ICM shock fronts characteristic of those induced by galaxy cluster mergers (incident Mach numbers 2 - 4). The simulated NAT formation was reported previously in \cite{on19a}. Our simulations utilize a three-dimensional, Eulerian ma… ▽ More We present a numerical study of the interactions between the elongated AGN outflows representing an evolved, narrow-angle tail (NAT) radio galaxy and planar, transverse ICM shock fronts characteristic of those induced by galaxy cluster mergers (incident Mach numbers 2 - 4). The simulated NAT formation was reported previously in \cite{on19a}. Our simulations utilize a three-dimensional, Eulerian magnetohydrodynamic code along with energy-dependent Eulerian transport of passive cosmic ray electrons. Our analysis of the shock/NAT interaction applies a Riemann problem-based theoretical model to interpret complex shock front behavior during passage through the highly heterogeneous structures of the simulated NAT tails. In addition to shock compression, shock-induced vortical motions are observed within the tails that contribute to coherent turbulent dynamo processes that continue to amplify the magnetic fields in the tails well after initial shock compression. We analyze synthetic radio observations spanning the NAT-shock interaction period, and examine the brightness, spectral and polarization properties of our shock-rejuvenated radio tails, as well as the extent to which the pre-shock states of the plasma and particle populations in our tails influence post-shock observations. Finally, we evaluate our findings in the possible context of a physical analogy to our simulated NAT providing the precursor to a cluster ``radio relic'' associated with an impacting ICM shock. △ Less

Submitted 21 October, 2019; v1 submitted 17 October, 2019; originally announced October 2019.

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

Simulated Interactions Between Radio Galaxies and Cluster Shocks -- 2: Jet Axes Orthogonal to Shock Normals

Authors: Chris Nolting, T. W. Jones, Brian O'Neill, P. J. Mendygral

Abstract: We report a 3D MHD simulation study of the interactions between radio galaxies and galaxy-cluster-media shocks in which the incident shock normals are orthogonal to the bipolar AGN jets. Before shock impact, light, supersonic jets inflate lobes (cavities) in a static, uniform ICM. We examine three AGN activity scenarios: 1) continued, steady jet activity; 2) jet source cycled off coincident with s… ▽ More We report a 3D MHD simulation study of the interactions between radio galaxies and galaxy-cluster-media shocks in which the incident shock normals are orthogonal to the bipolar AGN jets. Before shock impact, light, supersonic jets inflate lobes (cavities) in a static, uniform ICM. We examine three AGN activity scenarios: 1) continued, steady jet activity; 2) jet source cycled off coincident with shock/radio lobe impact; 3) jet activity ceased well before shock arrival (a "radio phoenix" scenario). The simulations follow relativistic electrons (CRe) introduced by the jets, enabling synthetic radio synchrotron images and spectra. Such encounters can be decomposed into an abrupt shock transition and a subsequent long term post shock wind. Shock impact disrupts the pre-formed, low density RG cavities into two ring vortices embedded in the post shock wind. Dynamical processes cause the vortex pair to merge as they propagate downwind somewhat faster than the wind itself. When the AGN jets remain active ram pressure bends the jets downwind, generating a narrow angle tail morphology aligned with the axis of the vortex ring. The deflected jets do not significantly alter dynamical evolution of the vortex ring. However, active jets and their associated tails do dominate the synchrotron emission, compromising the observability of the vortex structures. Downwind-directed momentum concentrated by the jets impacts and alters the post-encounter shock. In the "radio phoenix" scenario, no DSA of the fossil electron population is required to account for the observed brightening and flattening of the spectra, adiabatic compression effects are sufficient. △ Less

Submitted 18 September, 2019; originally announced September 2019.

Comments: 23 pages, 11 figures, Accepted for publication in the Astrophysical Journal

A Fresh Look at Narrow-Angle Tail Radio Galaxy Dynamics, Evolution and Emissions

Authors: Brian J. O'Neill, T. W. Jones, Chris Nolting, P. J. Mendygral

Abstract: We present a 3D magnetohydrodynamic (MHD) study of narrow-angle tail (NAT) radio galaxy (RG) dynamics, including passive cosmic ray electrons. We follow evolution of a bipolar-jet RG in a persistent crosswind through hundreds of Myr. We confirm traditional jet-bending models, while noting that our NAT exhibits a transitional morphology reminiscent of wide-angle radio tails. Once deflected, jets re… ▽ More We present a 3D magnetohydrodynamic (MHD) study of narrow-angle tail (NAT) radio galaxy (RG) dynamics, including passive cosmic ray electrons. We follow evolution of a bipolar-jet RG in a persistent crosswind through hundreds of Myr. We confirm traditional jet-bending models, while noting that our NAT exhibits a transitional morphology reminiscent of wide-angle radio tails. Once deflected, jets remain internally stable, but are intermittently disrupted by external disturbances induced by the NAT dynamics itself. The disruptions enhance jet and tail magnetic fields. Disrupted jet plasma is heterogeneously mixed with denser wind plasma, yielding patchy, filamentary tails that grow longer at a rate exceeding the wind speed. Such fast tail extension could, for example, allow NAT tails to overtake extraneous ICM features, such as shocks and shear layers downwind of where the tails first form. Those events, in turn, could generate enhanced radio emissions within the ICM features themselves that do not follow the geometrical extension of the tails past the encounter. Analysis of synthetic radio observations reveals an extended time period once the NAT has developed in which it displays a nearly steady-state morphology with integrated fluxes that are roughly constant, along with a self-similar, curved integrated spectrum. In an appendix, we outline a simple analytic jet trajectory formalism with one adjustable parameter, using it to illustrate explicit trajectories that extend the classic bending model to arbitrary jet-wind orientations. △ Less

Submitted 5 September, 2019; originally announced September 2019.

Comments: 34 pages, 16 figures, accepted for publication in Astrophysical Journal

Turbulent pressure support and hydrostatic mass-bias in the intracluster medium

Authors: M. Angelinelli, F. Vazza, C. Giocoli, S. Ettori, T. W. Jones, G. Brunetti, M. Brüggen, D. Eckert

Abstract: The degree of turbulent pressure support by residual gas motions in galaxy clusters is not well known. Mass modelling of combined X-ray and Sunyaev Zel'dovich observations provides an estimate of turbulent pressure support in the outer regions of several galaxy clusters. Here, we test two different filtering techniques to disentangle bulk from turbulent motions in non-radiative high-resolution cos… ▽ More The degree of turbulent pressure support by residual gas motions in galaxy clusters is not well known. Mass modelling of combined X-ray and Sunyaev Zel'dovich observations provides an estimate of turbulent pressure support in the outer regions of several galaxy clusters. Here, we test two different filtering techniques to disentangle bulk from turbulent motions in non-radiative high-resolution cosmological simulations of galaxy clusters using the cosmological hydro code ENZO. We find that the radial behavior of the ratio of non-thermal pressure to total gas pressure as a function of cluster-centric distance can be described by a simple polynomial function. The typical non-thermal pressure support in the centre of clusters is $\sim$5%, increasing to $\sim$15% in the outskirts, in line with the pressure excess found in recent X-ray observations. While the complex dynamics of the ICM makes it impossible to reconstruct a simple correlation between turbulent motions and hydrostatic bias, we find that a relation between them can be established using the median properties of a sample of objects. Moreover, we estimate the contribution of radial accelerations to the non-thermal pressure support and conclude that it decreases moving outwards from 40% (in the core) to 15% (in the cluster's outskirts). Adding this contribution to one provided by turbulence, we show that it might account for the entire observed hydrostatic bias in the innermost regions of the clusters, and for less than 80% of it at $r > 0.8 r_{200, m}$. △ Less

Submitted 6 April, 2020; v1 submitted 13 May, 2019; originally announced May 2019.

Comments: 20 pages; 21 figures; Substantial Revision; MNRAS in press

Interactions Between Radio Galaxies and Cluster Shocks - 1: Jet Axes Aligned with Shock Normals

Authors: Chris Nolting, T. W. Jones, Brian O'Neill, P. J. Mendygral

Abstract: We report from a study utilizing 3D MHD simulations, including cosmic-ray electrons, of the interactions between radio galaxies (RGs) and dynamically active ICMs. Here we consider interactions involving plane ICM shocks having Mach numbers 2--4 and their normals aligned with steady, active bipolar RG jets penetrating uniform, stationary ICMs. Shock impact disrupts the pre-formed RG jet cocoons int… ▽ More We report from a study utilizing 3D MHD simulations, including cosmic-ray electrons, of the interactions between radio galaxies (RGs) and dynamically active ICMs. Here we consider interactions involving plane ICM shocks having Mach numbers 2--4 and their normals aligned with steady, active bipolar RG jets penetrating uniform, stationary ICMs. Shock impact disrupts the pre-formed RG jet cocoons into ring vortex structures. Sufficiently strong post-shock winds can stop and even reverse the upwind jet, and strip jets to virtually "naked" states, leaving them without a surrounding cocoon. Strong shock-induced vorticity can also disrupt the downwind jet, so that the ring vortex remnant of the cocoons appears ahead of that jet's visible terminus. Magnetic field amplification in the ring vortex can significantly enhance its synchrotron emissions well after the vortex becomes isolated from the RG and its fresh CRe supply. We examine these dynamics and their observable consequences in detail. △ Less

Submitted 11 April, 2019; originally announced April 2019.

Comments: 29 Pages, 15 figures, accepted for publication in the Astrophysical Journal

Plasma 2020 - Intracluster Medium Plasmas

Authors: Damiano Caprioli, Gianfranco Brunetti, Thomas W. Jones, Hyesung Kang, Matthew Kunz, S. Peng Oh, Dongsu Ryu, Irina Zhuravleva, Ellen Zweibel

Abstract: Galaxy clusters are the largest and most massive bound objects resulting from cosmic hierarchical structure formation. Baryons account for somewhat more than 10% of that mass, with roughly 90% of the baryonic matter distributed throughout the clusters as hot ($T>1$ keV), high-$βべーた$, very weakly collisional plasma; the so-called "intracluster medium" (ICM). Cluster mergers, close gravitational encoun… ▽ More Galaxy clusters are the largest and most massive bound objects resulting from cosmic hierarchical structure formation. Baryons account for somewhat more than 10% of that mass, with roughly 90% of the baryonic matter distributed throughout the clusters as hot ($T>1$ keV), high-$βべーた$, very weakly collisional plasma; the so-called "intracluster medium" (ICM). Cluster mergers, close gravitational encounters and accretion, along with violent feedback from galaxies and relativistic jets from active galactic nuclei, drive winds, gravity waves, turbulence and shocks within the ICM. Those dynamics, in turn, generate cluster-scale magnetic fields and accelerate and mediate the transport of high-energy charged particles. Kinetic-scale, collective plasma processes define the basic character and fundamental signatures of these ICM phenomena, which are observed primarily by X-ray and radio astronomers. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: 5 pages, 2 figures, submitted for the 2020 Decadal Assessment of Plasma Science

Radio Galaxy Zoo: The Distortion of Radio Galaxies by Galaxy Clusters

Authors: Avery F. Garon, Lawrence Rudnick, O. Ivy Wong, Tom W. Jones, Jin-Ah Kim, Heinz Andernach, Stanislav S. Shabala, Anna D. Kapińska, Ray P. Norris, Francesco de Gasperin, Jean Tate, Hongming Tang

Abstract: We study the impact of cluster environment on the morphology of a sample of 4304 extended radio galaxies from Radio Galaxy Zoo. A total of 87% of the sample lies within a projected 15 Mpc of an optically identified cluster. Brightest cluster galaxies (BCGs) are more likely than other cluster members to be radio sources, and are also moderately bent. The surface density as a function of separation… ▽ More We study the impact of cluster environment on the morphology of a sample of 4304 extended radio galaxies from Radio Galaxy Zoo. A total of 87% of the sample lies within a projected 15 Mpc of an optically identified cluster. Brightest cluster galaxies (BCGs) are more likely than other cluster members to be radio sources, and are also moderately bent. The surface density as a function of separation from cluster center of non-BCG radio galaxies follows a power law with index $-1.10\pm 0.03$ out to $10~r_{500}$ ($\sim 7~$Mpc), which is steeper than the corresponding distribution for optically selected galaxies. Non-BCG radio galaxies are statistically more bent the closer they are to the cluster center. Within the inner $1.5~r_{500}$ ($\sim 1~$Mpc) of a cluster, non-BCG radio galaxies are statistically more bent in high-mass clusters than in low-mass clusters. Together, we find that non-BCG sources are statistically more bent in environments that exert greater ram pressure. We use the orientation of bent radio galaxies as an indicator of galaxy orbits and find that they are preferentially in radial orbits. Away from clusters, there is a large population of bent radio galaxies, limiting their use as cluster locators; however, they are still located within statistically overdense regions. We investigate the asymmetry in the tail length of sources that have their tails aligned along the radius vector from the cluster center, and find that the length of the inward-pointing tail is weakly suppressed for sources close to the center of the cluster. △ Less

Submitted 5 March, 2019; v1 submitted 16 January, 2019; originally announced January 2019.

Comments: 23 pages, 17 figures, 2 tables. Supplemental data files available in The Astronomical Journal or contact author

Journal ref: AJ 157 126 (2019)

A 1D fluid model of the Centaurus A jet

Authors: Sarka Wykes, Bradford T. Snios, Paul E. J. Nulsen, Ralph P. Kraft, Mark Birkinshaw, Martin J. Hardcastle, Diana M. Worrall, Iain McDonald, Marina Rejkuba, Thomas W. Jones, David J. Stark, William R. Forman, Eileen T. Meyer, Christine Jones

Abstract: We implement a steady, one-dimensional flow model for the X-ray jet of Centaurus A in which entrainment of stellar mass loss is the primary cause of dissipation. Using over 260 ks of new and archival Chandra/ACIS data, we have constrained the temperature, density and pressure distributions of gas in the central regions of the host galaxy of Centaurus A, and so the pressure throughout the length of… ▽ More We implement a steady, one-dimensional flow model for the X-ray jet of Centaurus A in which entrainment of stellar mass loss is the primary cause of dissipation. Using over 260 ks of new and archival Chandra/ACIS data, we have constrained the temperature, density and pressure distributions of gas in the central regions of the host galaxy of Centaurus A, and so the pressure throughout the length of its jet. The model is constrained by the observed profiles of pressure and jet width, and conserves matter and energy, enabling us to estimate jet velocities, and hence all the other flow properties. Invoking realistic stellar populations within the jet, we find that the increase in its momentum flux exceeds the net pressure force on the jet unless only about one half of the total stellar mass loss is entrained. For self-consistent models, the bulk speed only falls modestly, from ~0.67c to ~0.52c over the range of 0.25-5.94 kpc from the nucleus. The sonic Mach number varies between ~5.3 and 3.6 over this range. △ Less

Submitted 25 February, 2019; v1 submitted 11 December, 2018; originally announced December 2018.

Comments: 20 pages, 13 figures. Minor updates. Accepted by MNRAS

WENO-Wombat: Scalable Fifth-Order Constrained-Transport Magnetohydrodynamics for Astrophysical Applications

Authors: J. M. F. Donnert, H. Jang, P. Mendygral, G. Brunetti, D. Ryu, T. W. Jones

Abstract: Due to increase in computing power, high-order Eulerian schemes will likely become instrumental for the simulations of turbulence and magnetic field amplification in astrophysical fluids in the next years. We present the implementation of a fifth order weighted essentially non-oscillatory scheme for constrained-transport magnetohydrodynamics into the code WOMBAT. We establish the correctness of ou… ▽ More Due to increase in computing power, high-order Eulerian schemes will likely become instrumental for the simulations of turbulence and magnetic field amplification in astrophysical fluids in the next years. We present the implementation of a fifth order weighted essentially non-oscillatory scheme for constrained-transport magnetohydrodynamics into the code WOMBAT. We establish the correctness of our implementation with an extensive number tests. We find that the fifth order scheme performs as accurately as a common second order scheme at half the resolution. We argue that for a given solution quality the new scheme is more computationally efficient than lower order schemes in three dimensions. We also establish the performance characteristics of the solver in the WOMBAT framework. Our implementation fully vectorizes using flattened arrays in thread-local memory. It performs at about 0.6 Million zones per second per node on Intel Broadwell. We present scaling tests of the code up to 98 thousand cores on the Cray XC40 machine "Hazel Hen", with a sustained performance of about 5 percent of peak at scale. △ Less

Submitted 11 December, 2018; originally announced December 2018.

Comments: 35 pages, 32 figures

The turbulent pressure support in galaxy clusters revisited

Authors: F. Vazza, M. Angelinelli, T. W. Jones, D. Eckert, Brüggen M., Brunetti G., Gheller C

Abstract: Due to their late formation in cosmic history, clusters of galaxies are not fully in hydrostatic equilibrium and the gravitational pull of their mass at a given radius is expected not to be entirely balanced by the thermal gas pressure. Turbulence may supply additional pressure, and recent (X-ray and SZ) hydrostatic mass reconstructions claim a pressure support of $\sim 5-15\%$ of the total pressu… ▽ More Due to their late formation in cosmic history, clusters of galaxies are not fully in hydrostatic equilibrium and the gravitational pull of their mass at a given radius is expected not to be entirely balanced by the thermal gas pressure. Turbulence may supply additional pressure, and recent (X-ray and SZ) hydrostatic mass reconstructions claim a pressure support of $\sim 5-15\%$ of the total pressure at $R_{\rm 200}$. In this work we show that, after carefully disentangling bulk from small-scale turbulent motions in high-resolution simulations of galaxy clusters, we can constrain which fraction of the gas kinetic energy effectively provides pressure support in the cluster's gravitational potential. While the ubiquitous presence of radial inflows in the cluster can lead to significant bias in the estimate of the non-thermal pressure support, we report that only a part of this energy effectively acts as a source of pressure, providing a support of the order of $\sim 10\%$ of the total pressure at $R_{\rm 200}$. △ Less

Submitted 7 September, 2018; originally announced September 2018.

Comments: 5 pages, 5 pages, accepted, to appear in MNRAS Letters

Local Strain Heterogeneity Influences the Optoelectronic Properties of Halide Perovskites

Authors: Timothy W. Jones, Anna Osherov, Mejd Alsari, Melany Sponseller, Benjamin C. Duck, Young-Kwang Jung, Charles Settens, Farnaz Niroui, Roberto Brenes, Camelia V. Stan, Yao Li, Mojtaba Abdi-Jalebi, Nobumichi Tamura, J. Emyr Macdonald, Manfred Burghammer, Richard H. Friend, Vladimir Bulović, Aron Walsh, Gregory J. Wilson, Samuele Lilliu, Samuel D. Stranks

Abstract: Halide perovskites are promising semiconductors for inexpensive, high-performance optoelectronics. Despite a remarkable defect tolerance compared to conventional semiconductors, perovskite thin films still show substantial microscale heterogeneity in key properties such as luminescence efficiency and device performance. This behavior has been attributed to spatial fluctuations in the population of… ▽ More Halide perovskites are promising semiconductors for inexpensive, high-performance optoelectronics. Despite a remarkable defect tolerance compared to conventional semiconductors, perovskite thin films still show substantial microscale heterogeneity in key properties such as luminescence efficiency and device performance. This behavior has been attributed to spatial fluctuations in the population of sub-bandgap electronic states that act as trap-mediated non-radiative recombination sites. However, the origin of the variations, trap states and extent of the defect tolerance remains a topic of debate, and a precise understanding is critical to the rational design of defect management strategies. By combining scanning X-ray diffraction beamlines at two different synchrotrons with high-resolution transmission electron microscopy, we reveal levels of heterogeneity on the ten-micrometer scale (super-grains) and even ten-nanometer scale (sub-grain domains). We find that local strain is associated with enhanced defect concentrations, and correlations between the local structure and time-resolved photoluminescence reveal that these strain-related defects are the cause of non-radiative recombination. We reveal a direct connection between defect concentrations and non-radiative losses, as well as complex heterogeneity across multiple length scales, shedding new light on the presence and influence of structural defects in halide perovskites. △ Less

Submitted 15 May, 2018; v1 submitted 3 March, 2018; originally announced March 2018.

Acceleration of Cosmic Ray Electrons at Weak Shocks in Galaxy Clusters

Authors: Hyesung Kang, Dongsu Ryu, T. W. Jones

Abstract: According to structure formation simulations, weak shocks with typical Mach number, $M_{\rm s}\lesssim 3$, are expected to form in merging galaxy clusters. The presence of such shocks has been indicated by X-ray and radio observations of many merging clusters. In particular, diffuse radio sources known as radio relics could be explained by synchrotron-emitting electrons accelerated via diffusive s… ▽ More According to structure formation simulations, weak shocks with typical Mach number, $M_{\rm s}\lesssim 3$, are expected to form in merging galaxy clusters. The presence of such shocks has been indicated by X-ray and radio observations of many merging clusters. In particular, diffuse radio sources known as radio relics could be explained by synchrotron-emitting electrons accelerated via diffusive shock acceleration (Fermi I) at quasi-perpendicular shocks. Here we also consider possible roles of stochastic acceleration (Fermi II) by compressive MHD turbulence downstream of the shock. Then we explore a puzzling discrepancy that for some radio relics, the shock Mach number inferred from the radio spectral index is substantially larger than that estimated from X-ray observations. This problem could be understood, if shock surfaces associated with radio relics consist of multiple shocks with different strengths.In that case, X-ray observations tend to pick up the part of shocks with lower Mach numbers and higher kinetic energy flux, while radio emissions come preferentially from the part of shocks with higher Mach numbers and higher cosmic ray (CR) production. We also show that the Fermi I reacceleration model with preexisting fossil electrons supplemented by Fermi II acceleration due to postshock turbulence could reproduce observed profiles of radio flux densities and integrated radio spectra of two giant radio relics. This study demonstrates the CR electrons can be accelerated at collisionless shocks in galaxy clusters just like supernova remnant shock in the interstellar medium and interplanetary shocks in the solar wind. △ Less

Submitted 24 July, 2017; v1 submitted 21 July, 2017; originally announced July 2017.

Comments: 8 pages, 35th International Cosmic Ray Conference, Busan, Korea

Shock Acceleration Model for the Toothbrush Radio Relic

Authors: Hyesung Kang, Dongsu Ryu, T. W. Jones

Abstract: Although many of the observed properties of giant radio relics detected in the outskirts of galaxy clusters can be explained by relativistic electrons accelerated at merger-driven shocks, significant puzzles remain. In the case of the so-called Toothbrush relic, the shock Mach number estimated from X-ray observations ($M_{\rm X}\approx1.2-1.5$) is substantially weaker than that inferred from the r… ▽ More Although many of the observed properties of giant radio relics detected in the outskirts of galaxy clusters can be explained by relativistic electrons accelerated at merger-driven shocks, significant puzzles remain. In the case of the so-called Toothbrush relic, the shock Mach number estimated from X-ray observations ($M_{\rm X}\approx1.2-1.5$) is substantially weaker than that inferred from the radio spectral index ($M_{\rm rad}\approx2.8$).Toward understanding such a discrepancy, we here consider the following diffusive shock acceleration (DSA) models:(1) weak-shock models with $M_{\rm s}\lesssim 2$ and a preexisting population of cosmic-ray electrons (CRe) with a flat energy spectrum,and (2) strong-shock models with $M_{\rm s}\approx3$ and either shock-generated suprathermal electrons or preexisting fossil CRe. We calculate the synchrotron emission from the accelerated CRe, following the time evolution of the electron DSA, and subsequent radiative cooling and postshock turbulent acceleration (TA). We find that both models could reproduce reasonably well the observed integrated radio spectrum of the Toothbrush relic, but the observed broad transverse profile requires the stochastic acceleration by downstream turbulence, which we label "turbulent acceleration" or TA to distinguish it from DSA. Moreover, to account for the almost uniform radio spectral index profile along the length of the relic, the weak-shock models require a preshock region over 400~kpc with a uniform population of preexisting CRe with a high cutoff energy ($\gtrsim 40$ GeV). Due to the short cooling time, it is challenging to explain the origin of such energetic electrons. Therefore, we suggest the strong-shock models with low-energy seed CRe ($\lesssim 150$~MeV) are preferred for the radio observations of this relic. △ Less

Submitted 4 June, 2017; v1 submitted 1 March, 2017; originally announced March 2017.

Comments: Figures 5 and 6 are repalced. Several typos are corrected

WOMBAT: A Scalable and High Performance Astrophysical MHD Code

Authors: Peter Mendygral, Nick Radcliffe, Krishna Kandalla, David Porter, Brian J. O'Neill, Chris Nolting, Paul Edmon, Julius M. F. Donnert, Thomas W. Jones

Abstract: We present a new code for astrophysical magneto-hydrodynamics specifically designed and optimized for high performance and scaling on modern and future supercomputers. We describe a novel hybrid OpenMP/MPI programming model that emerged from a collaboration between Cray, Inc. and the University of Minnesota. This design utilizes MPI-RMA optimized for thread scaling, which allows the code to run ex… ▽ More We present a new code for astrophysical magneto-hydrodynamics specifically designed and optimized for high performance and scaling on modern and future supercomputers. We describe a novel hybrid OpenMP/MPI programming model that emerged from a collaboration between Cray, Inc. and the University of Minnesota. This design utilizes MPI-RMA optimized for thread scaling, which allows the code to run extremely efficiently at very high thread counts ideal for the latest generation of the multi-core and many-core architectures. Such performance characteristics are needed in the era of "exascale" computing. We describe and demonstrate our high-performance design in detail with the intent that it may be used as a model for other, future astrophysical codes intended for applications demanding exceptional performance. △ Less

Submitted 25 January, 2017; originally announced January 2017.

Comments: 28 pages, 23 figures, Accepted for publication in ApJS

Using Collisions of AGN Outflows with ICM Shocks as Dynamical Probes

Authors: T. W. Jones, Chris Nolting, B. J. O'Neill, P. J. Mendygral

Abstract: In this paper we lay out a simple set of relationships connecting the dynamics of fast plasma jets to the dynamical state of their ambient media. The objective is to provide a tool kit that can be used to connect the morphologies of radio AGNs in galaxy clusters to the dynamical state of the local ICM. The formalism is intended to apply to jets whether they are relativistic or non-relativistic. Sp… ▽ More In this paper we lay out a simple set of relationships connecting the dynamics of fast plasma jets to the dynamical state of their ambient media. The objective is to provide a tool kit that can be used to connect the morphologies of radio AGNs in galaxy clusters to the dynamical state of the local ICM. The formalism is intended to apply to jets whether they are relativistic or non-relativistic. Special attention is paid to interactions involving ICM shocks, although the results can be applied more broadly. Our formalism emphasizes the importance of the relative Mach number of the impacting ICM flow and the internal Mach number of the AGN jet in determining how the AGN outflows evolve. △ Less

Submitted 16 December, 2016; originally announced December 2016.

Comments: 11 pages, 5 figures. Submitted to Physics of Plasmas

Turbulence and Vorticity in Galaxy Clusters Generated by Structure Formation

Authors: F. Vazza, T. W. Jones, M. Brüggen, G. Brunetti, C. Gheller, D. Porter, D. Ryu

Abstract: Turbulence is a key ingredient for the evolution of the intracluster medium, whose properties can be predicted with high resolution numerical simulations. We present initial results on the generation of solenoidal and compressive turbulence in the intracluster medium during the formation of a small-size cluster using highly resolved, non-radiative cosmological simulations, with a refined monitorin… ▽ More Turbulence is a key ingredient for the evolution of the intracluster medium, whose properties can be predicted with high resolution numerical simulations. We present initial results on the generation of solenoidal and compressive turbulence in the intracluster medium during the formation of a small-size cluster using highly resolved, non-radiative cosmological simulations, with a refined monitoring in time. In this first of a series of papers, we closely look at one simulated cluster whose formation was distinguished by a merger around $z \sim 0.3$. We separate laminar gas motions, turbulence and shocks with dedicated filtering strategies and distinguish the solenoidal and compressive components of the gas flows using Hodge-Helmholtz decomposition. Solenoidal turbulence dominates the dissipation of turbulent motions ($\sim 95\%$) in the central cluster volume at all epochs. The dissipation via compressive modes is found to be more important ($\sim 30 \%$ of the total) only at large radii ($\geq 0.5 ~r_{\rm vir}$) and close to merger events. We show that enstrophy (vorticity squared) is good proxy of solenoidal turbulence. All terms ruling the evolution of enstrophy (i.e. baroclinic, compressive, stretching and advective terms) are found to be significant, but in amounts that vary with time and location. Two important trends for the growth of enstrophy in our simulation are identified: first, enstrophy is continuously accreted into the cluster from the outside, and most of that accreted enstrophy is generated near the outer accretion shocks by baroclinic and compressive processes. Second, in the cluster interior vortex stretching is dominant, although the other terms also contribute substantially. △ Less

Submitted 12 September, 2016; originally announced September 2016.

Comments: 21 pages, 21 figure. MNRAS accepted, in press

Radio Galaxy Zoo: discovery of a poor cluster through a giant wide-angle tail radio galaxy

Authors: J. K. Banfield, H. Andernach, A. D. Kapinska, L. Rudnick, M. J. Hardcastle, G. Cotter, S. Vaughan, T. W. Jones, I. Heywood, J. D. Wing, O. I. Wong, T. Matorny, I. A. Terentev, A. R. Lopez-Sanchez, R. P. Norris, N. Seymour, S. S. Shabala, K. W. Willett

Abstract: We have discovered a previously unreported poor cluster of galaxies (RGZ-CL J0823.2+0333) through an unusual giant wide-angle tail radio galaxy found in the Radio Galaxy Zoo project. We obtained a spectroscopic redshift of $z=0.0897$ for the E0-type host galaxy, 2MASX J08231289+0333016, leading to M$_r = -22.6$ and a $1.4\,$GHz radio luminosity density of $L_{\rm 1.4} = 5.5\times10^{24}$ W Hz… ▽ More We have discovered a previously unreported poor cluster of galaxies (RGZ-CL J0823.2+0333) through an unusual giant wide-angle tail radio galaxy found in the Radio Galaxy Zoo project. We obtained a spectroscopic redshift of $z=0.0897$ for the E0-type host galaxy, 2MASX J08231289+0333016, leading to M$_r = -22.6$ and a $1.4\,$GHz radio luminosity density of $L_{\rm 1.4} = 5.5\times10^{24}$ W Hz$^{-1}$. These radio and optical luminosities are typical for wide-angle tailed radio galaxies near the borderline between Fanaroff-Riley (FR) classes I and II. The projected largest angular size of $\approx8\,$arcmin corresponds to $800\,$kpc and the full length of the source along the curved jets/trails is $1.1\,$Mpc in projection. X-ray data from the XMM-Newton archive yield an upper limit on the X-ray luminosity of the thermal emission surrounding RGZ J082312.9+033301,at $1.2-2.6\times10^{43}$ erg s$^{-1}$ for assumed intra-cluster medium temperatures of $1.0-5.0\,$keV. Our analysis of the environment surrounding RGZ J082312.9+033301 indicates that RGZ J082312.9+033301 lies within a poor cluster. The observed radio morphology suggests that (a) the host galaxy is moving at a significant velocity with respect to an ambient medium like that of at least a poor cluster, and that (b) the source may have had two ignition events of the active galactic nucleus with $10^7\,$yrs in between. This reinforces the idea that an association between RGZ J082312.9+033301, and the newly discovered poor cluster exists. △ Less

Submitted 16 June, 2016; v1 submitted 15 June, 2016; originally announced June 2016.

Comments: Published in MNRAS

Journal ref: MNRAS, 460, 2376 - 2384 (2016)

The widest-frequency radio relic spectra: observations from 150 MHz to 30 GHz

Authors: Andra Stroe, Timothy Shimwell, Clare Rumsey, Reinout van Weeren, Maja Kierdorf, Julius Donnert, Thomas W. Jones, Huub J. A. Röttgering, Matthias Hoeft, Carmen Rodriguez-Gonzalvez, Jeremy J. Harwood, Richard D. E. Saunders

Abstract: Radio relics are patches of diffuse synchrotron radio emission that trace shock waves. Relics are thought to form when intra-cluster medium electrons are accelerated by cluster merger induced shock waves through the diffusive shock acceleration mechanism. In this paper, we present observations spanning 150 MHz to 30 GHz of the `Sausage' and `Toothbrush' relics from the Giant Metrewave and Westerbo… ▽ More Radio relics are patches of diffuse synchrotron radio emission that trace shock waves. Relics are thought to form when intra-cluster medium electrons are accelerated by cluster merger induced shock waves through the diffusive shock acceleration mechanism. In this paper, we present observations spanning 150 MHz to 30 GHz of the `Sausage' and `Toothbrush' relics from the Giant Metrewave and Westerbork telescopes, the Karl G. Jansky Very Large Array, the Effelsberg telescope, the Arcminute Microkelvin Imager and Combined Array for Research in Millimeter-wave Astronomy. We detect both relics at 30 GHz, where the previous highest frequency detection was at 16 GHz. The integrated radio spectra of both sources clearly steepen above 2 GHz, at the >6$σしぐま$ significance level, supports the spectral steepening previously found in the `Sausage' and the Abell 2256 relic. Our results challenge the widely adopted simple formation mechanism of radio relics and suggest more complicated models have to be developed that, for example, involve re-acceleration of aged seed electrons. △ Less

Submitted 22 October, 2015; originally announced October 2015.

Comments: Accepted for publication in MNRAS; 16 pages, 11 figures;

Vorticity, Shocks and Magnetic Fields in Subsonic, ICM-like Turbulence

Authors: David H. Porter, T. W. Jones, Dongsu Ryu

Abstract: We analyze high resolution simulations of compressible, MHD turbulence with properties resembling conditions in galaxy clusters. The flow is driven to turbulence Mach number $\mathcal{M}_t \sim 1/2$ in an isothermal medium with an initially very weak, uniform seed magnetic field ($βべーた= P_g/P_B = 10^6$). Since cluster turbulence is likely to result from a mix of sheared (solenoidal) and compressive f… ▽ More We analyze high resolution simulations of compressible, MHD turbulence with properties resembling conditions in galaxy clusters. The flow is driven to turbulence Mach number $\mathcal{M}_t \sim 1/2$ in an isothermal medium with an initially very weak, uniform seed magnetic field ($βべーた= P_g/P_B = 10^6$). Since cluster turbulence is likely to result from a mix of sheared (solenoidal) and compressive forcing processes, we examine the distinct turbulence properties for both cases. In one set of simulations velocity forcing is entirely solenoidal ($\nabla\cdot δでるた{\vec u} = 0$), while in the other it is entirely compressive ($\nabla\times δでるた{\vec u} = 0$). Both cases develop a mixture of solenoidal and compressive turbulent motions, since each generates the other. The development of compressive turbulent motions leads to shocks, even when the turbulence is solenoidally forced and subsonic. Shocks, in turn, produce and amplify vorticity, which is especially important in compressively forced turbulence. To clarify those processes we include a pair of appendices that look in detail at vorticity evolution in association with shocks. From our simulation analyses we find that magnetic fields amplified to near saturation levels in predominantly solenoidal turbulence can actually enhance vorticity on small scales by concentrating and stabilizing shear. The properties, evolution rates and relative contributions of the kinetic and magnetic turbulent elements depend strongly on the character of the forcing. Specifically, shocks are stronger, but vorticity evolution and magnetic field amplification are slower and weaker when the turbulence is compressively forced. We identify a simple relation to estimate characteristic shock strengths in terms of the turbulence Mach number and the character of the forcing. Our results will be helpful in understanding flow motions in galaxy clusters. △ Less

Submitted 30 July, 2015; originally announced July 2015.

Comments: 46 pages, 14 figures. Accepted for publication in the Astrophysical Journal

Radio Galaxy Zoo: host galaxies and radio morphologies derived from visual inspection

Authors: J. K. Banfield, O. I. Wong, K. W. Willett, R. P. Norris, L. Rudnick, S. S. Shabala, B. D. Simmons, C. Snyder, A. Garon, N. Seymour, E. Middelberg, H. Andernach, C. J. Lintott, K. Jacob, A. D. Kapinska, M. Y. Mao, K. L. Masters, M. J. Jarvis, K. Schawinski, E. Paget, R. Simpson, H. R. Klockner, S. Bamford, T. Burchell, K. E. Chow , et al. (11 additional authors not shown)

Abstract: We present results from the first twelve months of operation of Radio Galaxy Zoo, which upon completion will enable visual inspection of over 170,000 radio sources to determine the host galaxy of the radio emission and the radio morphology. Radio Galaxy Zoo uses $1.4\,$GHz radio images from both the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) and the Australia Telescope Large Area… ▽ More We present results from the first twelve months of operation of Radio Galaxy Zoo, which upon completion will enable visual inspection of over 170,000 radio sources to determine the host galaxy of the radio emission and the radio morphology. Radio Galaxy Zoo uses $1.4\,$GHz radio images from both the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) and the Australia Telescope Large Area Survey (ATLAS) in combination with mid-infrared images at $3.4\,μみゅー$m from the {\it Wide-field Infrared Survey Explorer} (WISE) and at $3.6\,μみゅー$m from the {\it Spitzer Space Telescope}. We present the early analysis of the WISE mid-infrared colours of the host galaxies. For images in which there is $>\,75\%$ consensus among the Radio Galaxy Zoo cross-identifications, the project participants are as effective as the science experts at identifying the host galaxies. The majority of the identified host galaxies reside in the mid-infrared colour space dominated by elliptical galaxies, quasi-stellar objects (QSOs), and luminous infrared radio galaxies (LIRGs). We also find a distinct population of Radio Galaxy Zoo host galaxies residing in a redder mid-infrared colour space consisting of star-forming galaxies and/or dust-enhanced non star-forming galaxies consistent with a scenario of merger-driven active galactic nuclei (AGN) formation. The completion of the full Radio Galaxy Zoo project will measure the relative populations of these hosts as a function of radio morphology and power while providing an avenue for the identification of rare and extreme radio structures. Currently, we are investigating candidates for radio galaxies with extreme morphologies, such as giant radio galaxies, late-type host galaxies with extended radio emission, and hybrid morphology radio sources. △ Less

Submitted 26 July, 2015; originally announced July 2015.

Comments: Accepted to MNRAS, 17 pages, 13 figures, 1 table. Higher resolution figures available with MNRAS publication or contact author

Injection of $κかっぱ$-like Suprathermal Particles into Diffusive Shock Acceleration

Authors: Hyesung Kang, Vahe Petrosian, Dongsu Ryu, T. W. Jones

Abstract: We consider a phenomenological model for the thermal leakage injection in the diffusive shock acceleration (DSA) process, in which suprathermal protons and electrons near the shock transition zone are assumed to have the so-called $κかっぱ$-distributions produced by interactions of background thermal particles with pre-existing and/or self-excited plasma/MHD waves or turbulence. The $κかっぱ$-distribution has… ▽ More We consider a phenomenological model for the thermal leakage injection in the diffusive shock acceleration (DSA) process, in which suprathermal protons and electrons near the shock transition zone are assumed to have the so-called $κかっぱ$-distributions produced by interactions of background thermal particles with pre-existing and/or self-excited plasma/MHD waves or turbulence. The $κかっぱ$-distribution has a power-law tail, instead of an exponential cutoff, well above the thermal peak momentum. So there are a larger number of potential seed particles with momentum, above that required for participation in the DSA process. As a result, the injection fraction for the $κかっぱ$-distribution depends on the shock Mach number much less severely compared to that for the Maxwellian distribution. Thus, the existence of $κかっぱ$-like suprathermal tails at shocks would ease the problem of extremely low injection fractions, especially for electrons and especially at weak shocks such as those found in the intracluster medium. We suggest that the injection fraction for protons ranges $10^{-4}-10^{-3}$ for a $κかっぱ$-distribution with $10 < κかっぱ_p < 30$ at quasi-parallel shocks, while the injection fraction for electrons becomes $10^{-6}-10^{-5}$ for a $κかっぱ$-distribution with $κかっぱ_e < 2$ at quasi-perpendicular shocks. For such $κかっぱ$ values the ratio of cosmic ray electrons to protons naturally becomes $K_{e/p}\sim 10^{-3}-10^{-2}$, which is required to explain the observed ratio for Galactic cosmic rays. △ Less

Submitted 3 May, 2014; originally announced May 2014.

Comments: 10 pages, 4 figures, submitted to Astrophysical Journal

Filaments in the southern giant lobe of Centaurus A: constraints on nature and origin from modelling and GMRT observations

Authors: Sarka Wykes, Huib T. Intema, Martin J. Hardcastle, Abraham Achterberg, Thomas W. Jones, Helmut Jerjen, Emanuela Orru, Alex Lazarian, Timothy W. Shimwell, Michael W. Wise, Philipp P. Kronberg

Abstract: We present results from imaging of the radio filaments in the southern giant lobe of Centaurus A using data from GMRT observations at 325 and 235 MHz, and outcomes from filament modelling. The observations reveal a rich filamentary structure, largely matching the morphology at 1.4 GHz. We find no clear connection of the filaments to the jet. We seek to constrain the nature and origin of the vertex… ▽ More We present results from imaging of the radio filaments in the southern giant lobe of Centaurus A using data from GMRT observations at 325 and 235 MHz, and outcomes from filament modelling. The observations reveal a rich filamentary structure, largely matching the morphology at 1.4 GHz. We find no clear connection of the filaments to the jet. We seek to constrain the nature and origin of the vertex and vortex filaments associated with the lobe and their role in high-energy particle acceleration. We deduce that these filaments are at most mildly overpressured with respect to the global lobe plasma showing no evidence of large-scale efficient Fermi I-type particle acceleration, and persist for ~ 2-3 Myr. We demonstrate that the dwarf galaxy KK 196 (AM 1318-444) cannot account for the features, and that surface plasma instabilities, the internal sausage mode and radiative instabilities are highly unlikely. An internal tearing instability and the kink mode are allowed within the observational and growth time constraints and could develop in parallel on different physical scales. We interpret the origin of the vertex and vortex filaments in terms of weak shocks from transonic MHD turbulence or from a moderately recent jet activity of the parent AGN, or an interplay of both. △ Less

Submitted 23 May, 2014; v1 submitted 9 April, 2014; originally announced April 2014.

Comments: 17 pages, 4 figures. Minor updates to match accepted version. Accepted by MNRAS

Cosmic rays in galaxy clusters and their non-thermal emission

Authors: G. Brunetti, T. W. Jones

Abstract: Radio observations prove the existence of relativistic particles and magnetic field associated with the intra-cluster-medium (ICM) through the presence of extended synchrotron emission in the form of radio halos and peripheral relics. This observational evidence has fundamental implications on the physics of the ICM. Non-thermal components in galaxy clusters are indeed unique probes of very energe… ▽ More Radio observations prove the existence of relativistic particles and magnetic field associated with the intra-cluster-medium (ICM) through the presence of extended synchrotron emission in the form of radio halos and peripheral relics. This observational evidence has fundamental implications on the physics of the ICM. Non-thermal components in galaxy clusters are indeed unique probes of very energetic processes operating within clusters that drain gravitational and electromagnetic energy into cosmic rays and magnetic fields. These components strongly affect the (micro-)physical properties of the ICM, including viscosity and electrical conductivities, and have also potential consequences on the evolution of clusters themselves. The nature and properties of cosmic rays in galaxy clusters, including the origin of the observed radio emission on cluster-scales, have triggered an active theoretical debate in the last decade. Only recently we can start addressing some of the most important questions in this field, thanks to recent observational advances, both in the radio and at high energies. The properties of cosmic rays and of cluster non-thermal emissions depend on the dynamical state of the ICM, the efficiency of particle acceleration mechanisms in the ICM and on the dynamics of these cosmic rays. In this review we discuss in some detail the acceleration and transport of cosmic rays in galaxy clusters and the most relevant observational milestones that have provided important steps on our understanding of this physics. Finally, looking forward to the possibilities from new generations of observational tools, we focus on what appear to be the most important prospects for the near future from radio and high-energy observations. △ Less

Submitted 29 January, 2014; originally announced January 2014.

Comments: Invited review paper, 74 pages, 21 figures, International Journal of Modern Physics D, in press

Nonthermal Radiation from Supernova Remnants: Effects of Magnetic Field Amplification and Particle Escape

Authors: Hyesung Kang, T. W. Jones, Paul P. Edmon

Abstract: We explore nonlinear effects of wave-particle interactions on the diffusive shock acceleration (DSA) process in Type Ia-like, SNR blast waves, by implementing phenomenological models for magnetic field amplification, Alfv'enic drift, and particle escape in time-dependent numerical simulations of nonlinear DSA. For typical SNR parameters the CR protons can be accelerated to PeV energies only if the… ▽ More We explore nonlinear effects of wave-particle interactions on the diffusive shock acceleration (DSA) process in Type Ia-like, SNR blast waves, by implementing phenomenological models for magnetic field amplification, Alfv'enic drift, and particle escape in time-dependent numerical simulations of nonlinear DSA. For typical SNR parameters the CR protons can be accelerated to PeV energies only if the region of amplified field ahead of the shock is extensive enough to contain the diffusion lengths of the particles of interest. Even with the help of Alfv'enic drift, it remains somewhat challenging to construct a nonlinear DSA model for SNRs in which order of 10 % of the supernova explosion energy is converted to the CR energy and the magnetic field is amplified by a factor of 10 or so in the shock precursor, while, at the same time, the energy spectrum of PeV protons is steeper than E^{-2}. To explore the influence of these physical effects on observed SNR emissions, we also compute resulting radio-to-gamma-ray spectra. Nonthermal emission spectra, especially in X-ray and gamma-ray bands,depend on the time dependent evolution of CR injection process, magnetic field amplification, and particle escape, as well as the shock dynamic evolution. This result comes from the fact that the high energy end of the CR spectrum is composed of the particles that are injected in the very early stages of blast wave evolution. Thus it is crucial to understand better the plasma wave-particle interactions associated with collisionless shocks in detail modeling of nonthermal radiation from SNRs. △ Less

Submitted 30 August, 2013; originally announced August 2013.

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

Production of Z0 bosons in elastic and quasi-elastic ep collisions at HERA

Authors: ZEUS collaboration, H. Abramowicz, I. Abt, L. Adamczyk, M. Adamus, R. Aggarwal, S. Antonelli, P. Antonioli, A. Antonov, M. Arneodo, O. Arslan, V. Aushev, Y. Aushev, O. Bachynska, A. Bamberger, A. N. Barakbaev, G. Barbagli, G. Bari, F. Barreiro, N. Bartosik, D. Bartsch, M. Basile, O. Behnke, J. Behr, U. Behrens , et al. (278 additional authors not shown)

Abstract: The production of Z0 bosons in the reaction ep -> eZ0p*, where p* stands for a proton or a low-mass nucleon resonance, has been studied in ep collisions at HERA using the ZEUS detector. The analysis is based on a data sample collected between 1996 and 2007, amounting to 496 pb-1 of integrated luminosity. The Z0 was measured in the hadronic decay mode. The elasticity of the events was ensured by a… ▽ More The production of Z0 bosons in the reaction ep -> eZ0p*, where p* stands for a proton or a low-mass nucleon resonance, has been studied in ep collisions at HERA using the ZEUS detector. The analysis is based on a data sample collected between 1996 and 2007, amounting to 496 pb-1 of integrated luminosity. The Z0 was measured in the hadronic decay mode. The elasticity of the events was ensured by a cut on eta_max < 3.0, where eta_max is the maximum pseudorapidity of energy deposits in the calorimeter defined with respect to the proton beam direction. A signal was observed at the Z0 mass. The cross section of the reaction ep -> eZ0p* was measured to be sigma(ep -> eZ0p*) = 0.13 +/- 0.06 (stat.) +/- 0.01 (syst.) pb, in agreement with the Standard Model prediction of 0.16 pb. This is the first measurement of Z0 production in ep collisions. △ Less

Submitted 19 October, 2012; originally announced October 2012.

Comments: 13 pages, 3 figures, submitted to PLB

Report number: DESY-12-168

Measurement of high-Q2 neutral current deep inelastic e+p scattering cross sections with a longitudinally polarised positron beam at HERA

Authors: ZEUS Collaboration, H. Abramowicz, I. Abt, L. Adamczyk, M. Adamus, R. Aggarwal, S. Antonelli, P. Antonioli, A. Antonov, M. Arneodo, O. Arslan, V. Aushev, Y. Aushev, O. Bachynska, A. Bamberger, A. N. Barakbaev, G. Barbagli, G. Bari, F. Barreiro, N. Bartosik, D. Bartsch, M. Basile, O. Behnke, J. Behr, U. Behrens , et al. (278 additional authors not shown)

Abstract: Measurements of neutral current cross sections for deep inelastic scattering in e+p collisions at HERA with a longitudinally polarised positron beam are presented. The single-differential cross-sections d(sigma)/dQ2, d(sigma)/dx and d(sigma)/dy and the reduced cross-section were measured in the kinematic region Q2 > 185 GeV2 and y < 0.9, where Q2 is the four-momentum transfer squared, x the Bjorke… ▽ More Measurements of neutral current cross sections for deep inelastic scattering in e+p collisions at HERA with a longitudinally polarised positron beam are presented. The single-differential cross-sections d(sigma)/dQ2, d(sigma)/dx and d(sigma)/dy and the reduced cross-section were measured in the kinematic region Q2 > 185 GeV2 and y < 0.9, where Q2 is the four-momentum transfer squared, x the Bjorken scaling variable, and y the inelasticity of the interaction. The measurements were performed separately for positively and negatively polarised positron beams. The measurements are based on an integrated luminosity of 135.5 pb-1 collected with the ZEUS detector in 2006 and 2007 at a centre-of-mass energy of 318 GeV. The structure functions F3 and F3(gamma)Z were determined by combining the e+p results presented in this paper with previously published e-p neutral current results. The asymmetry parameter A+ is used to demonstrate the parity violation predicted in electroweak interactions. The measurements are well described by the predictions of the Standard Model. △ Less

Submitted 12 May, 2014; v1 submitted 30 August, 2012; originally announced August 2012.

Comments: 46 pages, 10 figures and 15 tables. Version updated after journal review and including full list of authors

Report number: DESY-12-145

Journal ref: Phys Rev D 87 (2013) 052014

Inclusive-jet photoproduction at HERA and determination of alphas

Authors: ZEUS Collaboration, H. Abramowicz, I. Abt, L. Adamczyk, M. Adamus, R. Aggarwal, S. Antonelli, P. Antonioli, A. Antonov, M. Arneodo, V. Aushev, Y. Aushev, O. Bachynska, A. Bamberger, A. N. Barakbaev, G. Barbagli, G. Bari, F. Barreiro, N. Bartosik, D. Bartsch, M. Basile, O. Behnke, J. Behr, U. Behrens, L. Bellagamba , et al. (281 additional authors not shown)

Abstract: Inclusive-jet cross sections have been measured in the reaction ep->e+jet+X for photon virtuality Q2 < 1 GeV2 and gamma-p centre-of-mass energies in the region 142 < W(gamma-p) < 293 GeV with the ZEUS detector at HERA using an integrated luminosity of 300 pb-1. Jets were identified using the kT, anti-kT or SIScone jet algorithms in the laboratory frame. Single-differential cross sections are prese… ▽ More Inclusive-jet cross sections have been measured in the reaction ep->e+jet+X for photon virtuality Q2 < 1 GeV2 and gamma-p centre-of-mass energies in the region 142 < W(gamma-p) < 293 GeV with the ZEUS detector at HERA using an integrated luminosity of 300 pb-1. Jets were identified using the kT, anti-kT or SIScone jet algorithms in the laboratory frame. Single-differential cross sections are presented as functions of the jet transverse energy, ETjet, and pseudorapidity, etajet, for jets with ETjet > 17 GeV and -1 < etajet < 2.5. In addition, measurements of double-differential inclusive-jet cross sections are presented as functions of ETjet in different regions of etajet. Next-to-leading-order QCD calculations give a good description of the measurements, except for jets with low ETjet and high etajet. The influence of non-perturbative effects not related to hadronisation was studied. Measurements of the ratios of cross sections using different jet algorithms are also presented; the measured ratios are well described by calculations including up to O(alphas2) terms. Values of alphas(Mz) were extracted from the measurements and the energy-scale dependence of the coupling was determined. The value of alphas(Mz) extracted from the measurements based on the kT jet algorithm is alphas(Mz) = 0.1206 +0.0023 -0.0022 (exp.) +0.0042 -0.0035 (th.); the results from the anti-kT and SIScone algorithms are compatible with this value and have a similar precision. △ Less

Submitted 28 May, 2012; originally announced May 2012.

Comments: 53 pages, 17 figures

Report number: DESY-12-045

Diffusive Shock Acceleration Simulations of Radio Relics

Authors: Hyesung Kang, Dongsu Ryu, T. W. Jones

Abstract: Recent radio observations have identified a class of structures, so-called radio relics, in clusters of galaxies. The radio emission from these sources is interpreted as synchrotron radiation from GeV electrons gyrating in microG-level magnetic fields. Radio relics, located mostly in the outskirts of clusters, seem to associate with shock waves, especially those developed during mergers. In fact,… ▽ More Recent radio observations have identified a class of structures, so-called radio relics, in clusters of galaxies. The radio emission from these sources is interpreted as synchrotron radiation from GeV electrons gyrating in microG-level magnetic fields. Radio relics, located mostly in the outskirts of clusters, seem to associate with shock waves, especially those developed during mergers. In fact, they seem to be good structures to identify and probe such shocks in intracluster media (ICMs), provided we understand the electron acceleration and re-acceleration at those shocks. In this paper, we describe time-dependent simulations for diffusive shock acceleration at weak shocks that are expected to be found in ICMs. Freshly injected as well as pre-existing populations of cosmic-ray (CR) electrons are considered, and energy losses via synchrotron and inverse Compton are included. We then compare the synchrotron flux and spectral distributions estimated from the simulations with those in two well-observed radio relics in CIZA J2242.8+5301 and ZwCl0008.8+5215. Considering that the CR electron injection is rather inefficient at weak shocks with Mach number M <~ a few, the existence of radio relics could indicate the pre-existing population of low-energy CR electrons in ICMs. The implication of our results on the merger shock scenario of radio relics is discussed. △ Less

Submitted 9 May, 2012; originally announced May 2012.

Comments: Submitted to ApJ; 26 pages with 5 figures. Pdf can be downloaded from http://canopus.cnu.ac.kr/ryu/relics.pdf

Exclusive electroproduction of two pions at HERA

Authors: ZEUS collaboration, H. Abramowicz, I. Abt, L. Adamczyk, M. Adamus, R. Aggarwal, S. Antonelli, P. Antonioli, A. Antonov, M. Arneodo, D. Ashery, V. Aushev, Y. Aushev, O. Bachynska, A. Bamberger, A. N. Barakbaev, G. Barbagli, G. Bari, F. Barreiro, N. Bartosik, D. Bartsch, M. Basile, O. Behnke, J. Behr, U. Behrens , et al. (280 additional authors not shown)

Abstract: The exclusive electroproduction of two pions in the mass range 0.4 < Mπぱいπぱい < 2.5 GeV has been studied with the ZEUS detector at HERA using an integrated luminosity of 82 pb-1. The analysis was carried out in the kinematic range of 2 < Q2 < 80 GeV2, 32 < W < 180 GeV and |t| < 0.6 GeV2, where Q2 is the photon virtuality, W is the photon-proton centre-of-mass energy and t is the squared four-momentum t… ▽ More The exclusive electroproduction of two pions in the mass range 0.4 < Mπぱいπぱい < 2.5 GeV has been studied with the ZEUS detector at HERA using an integrated luminosity of 82 pb-1. The analysis was carried out in the kinematic range of 2 < Q2 < 80 GeV2, 32 < W < 180 GeV and |t| < 0.6 GeV2, where Q2 is the photon virtuality, W is the photon-proton centre-of-mass energy and t is the squared four-momentum transfer at the proton vertex. The two-pion invariant-mass distribution is interpreted in terms of the pion electromagnetic form factor, |F(Mπぱいπぱい)|, assuming that the studied mass range includes the contributions of the ρろー, ρろー' and ρろー" vector-meson states. The masses and widths of the resonances were obtained and the Q2 dependence of the cross-section ratios σしぐま(ρろー' \rightarrow πぱいπぱい)/σしぐま(ρろー) and σしぐま(ρろー" \rightarrow πぱいπぱい)/σしぐま(ρろー) was extracted. The pion form factor obtained in the present analysis is compared to that obtained in e+e- \rightarrow πぱい+πぱい-. △ Less

Submitted 21 November, 2011; originally announced November 2011.

Comments: 24 pages, including 11 figures

Report number: DESY-11-220

Search for single-top production in ep collisions at HERA

Authors: ZEUS Collaboration, H. Abramowicz, I. Abt, L. Adamczyk, M. Adamus, R. Aggarwal, S. Antonelli, P. Antonioli, A. Antonov, M. Arneodo, V. Aushev, Y. Aushev, O. Bachynska, A. Bamberger, A. N. Barakbaev, G. Barbagli, G. Bari, F. Barreiro, N. Bartosik, D. Bartsch, M. Basile, O. Behnke, J. Behr, U . Behrens, L. Bellagamba , et al. (278 additional authors not shown)

Abstract: A search for single-top production, $ep \rightarrow etX$, has been performed with the ZEUS detector at HERA using data corresponding to an integrated luminosity of $0.37\fbi$. No evidence for top production was found, consistent with the expectation from the Standard Model. Limits were computed for single-top production via flavour changing neutral current transitions. The result was combined with… ▽ More A search for single-top production, $ep \rightarrow etX$, has been performed with the ZEUS detector at HERA using data corresponding to an integrated luminosity of $0.37\fbi$. No evidence for top production was found, consistent with the expectation from the Standard Model. Limits were computed for single-top production via flavour changing neutral current transitions. The result was combined with a previous ZEUS result yielding a total luminosity of 0.50fb-1. A 95% credibility level upper limit of 0.13 pb was obtained for the cross section at the centre-of-mass energy of $\sqrt{s}=315\gev$. △ Less

Submitted 4 February, 2012; v1 submitted 16 November, 2011; originally announced November 2011.

Comments: Included approximate NLO and NNLO corrections to the cross section calculation used to derive coupling limit for the top production mediated by gamma exchange. 30 pages, 5 figures, 3 tables. Published in Phys. Lett. B 708 (2012) 27

Report number: DESY-11-204

Journal ref: Phys. Lett. B 708 (2012) 27-36

Scaled momentum distributions for K0s and Lambda/bar Lambda in DIS at HERA

Authors: ZEUS Collaboration, H. Abramowicz, I. Abt, L. Adamczyk, M. Adamus, R. Aggarwal, S. Antonelli, P. Antonioli, A. Antonov, M. Arneodo, V. Aushev, Y. Aushev, O. Bachynska, A. Bamberger, A. N. Barakbaev, G. Barbagli, G. Bari, F. Barreiro, N. Bartosik, D. Bartsch, M. Basile, O. Behnke, J. Behr, U. Behrens, L. Bellagamba , et al. (278 additional authors not shown)

Abstract: Scaled momentum distributions for the strange hadrons K0s and Lambda/bar Lambda were measured in deep inelastic ep scattering with the ZEUS detector at HERA using an integrated luminosity of 330 pb-1. The evolution of these distributions with the photon virtuality, Q2, was studied in the kinematic region 10<Q2<40000 GeV2 and 0.001<x<0.75, where x is the Bjorken scaling variable. Clear scaling viol… ▽ More Scaled momentum distributions for the strange hadrons K0s and Lambda/bar Lambda were measured in deep inelastic ep scattering with the ZEUS detector at HERA using an integrated luminosity of 330 pb-1. The evolution of these distributions with the photon virtuality, Q2, was studied in the kinematic region 10<Q2<40000 GeV2 and 0.001<x<0.75, where x is the Bjorken scaling variable. Clear scaling violations are observed. Predictions based on different approaches to fragmentation were compared to the measurements. Leading-logarithm parton-shower Monte Carlo calculations interfaced to the Lund string fragmentation model describe the data reasonably well in the whole range measured. Next-to-leading-order QCD calculations based on fragmentation functions, FFs, extracted from e+e- data alone, fail to describe the measurements. The calculations based on FFs extracted from a global analysis including e+e-, ep and pp data give an improved description. The measurements presented in this paper have the potential to further constrain the FFs of quarks, anti-quarks and gluons yielding K0s and Lambda/bar Lambda strange hadrons. △ Less

Submitted 19 April, 2012; v1 submitted 15 November, 2011; originally announced November 2011.

Comments: 34 pages, 8 figures

Report number: DESY-11-205

Measurement of the t dependence in exclusive photoproduction of Upsilon(1S) mesons at HERA

Authors: ZEUS collaboration, H. Abramowicz, I. Abt, L. Adamczyk, M. Adamus, R. Aggarwal, S. Antonelli, P. Antonioli, A. Antonov, M. Arneodo, V. Aushev, Y. Aushev, O. Bachynska, A. Bamberger, A. N. Barakbaev, G. Barbagli, G. Bari, F. Barreiro, N. Bartosik, D. Bartsch, M. Basile, O. Behnke, J. Behr, U. Behrens, L. Bellagamba , et al. (278 additional authors not shown)

Abstract: The exclusive photoproduction reaction gamma p -> Upsilon(1S) p has been studied with the ZEUS detector in ep collisions at HERA using an integrated luminosity of 468 pb^-1. The measurement covers the kinematic range 60<W<220 GeV and Q^2<1 GeV^2, where W is the photon-proton centre-of-mass energy and Q^2 is the photon virtuality. The exponential slope, b, of the t dependence of the cross section,… ▽ More The exclusive photoproduction reaction gamma p -> Upsilon(1S) p has been studied with the ZEUS detector in ep collisions at HERA using an integrated luminosity of 468 pb^-1. The measurement covers the kinematic range 60<W<220 GeV and Q^2<1 GeV^2, where W is the photon-proton centre-of-mass energy and Q^2 is the photon virtuality. The exponential slope, b, of the t dependence of the cross section, where t is the squared four-momentum transfer at the proton vertex, has been measured, yielding b = 4.3 +2.0 -1.3 (stat.) +0.5 -0.6 (syst.) GeV^-2. This constitutes the first measurement of the t dependence of the gamma p -> Upsilon(1S) p cross section. △ Less

Submitted 4 February, 2012; v1 submitted 9 November, 2011; originally announced November 2011.

Comments: 15 pages, 4 figures; http://dx.doi.org/10.1016/j.physletb.2012.01.009

Report number: DESY-11-186

Journal ref: Physics Letters B 708 (2012) 14-20

MHD Turbulence Simulation in a Cosmic Structure Context

Authors: T. W. Jones, D. H. Porter, D. Ryu, J. Cho

Abstract: The gaseous media of galaxy clusters and cosmic filaments, which constitute most of the baryonic matter in the universe, is highly dynamic. It is also probably turbulent, although the turbulence properties are poorly known. The gas is highly rarefied, essentially fully ionized plasma. Observational evidence suggests intracluster media (ICMs) are magnetized at some level. There are several possible… ▽ More The gaseous media of galaxy clusters and cosmic filaments, which constitute most of the baryonic matter in the universe, is highly dynamic. It is also probably turbulent, although the turbulence properties are poorly known. The gas is highly rarefied, essentially fully ionized plasma. Observational evidence suggests intracluster media (ICMs) are magnetized at some level. There are several possible origins for ICM seed fields; the observed fields are likely the result of turbulence in the ICM. We are engaged in a simulation study designed to understand in this context how very weak initial magnetic fields evolve in driven turbulence. We find that the magnetic fields eventually evolve towards equipartition levels with the vortical, solenoidal kinetic energy in the turbulence. As they do so the topology of the field structures transition from filamentary forms into ribbon-like structures in which the field orientations are laminated with vorticity structures. △ Less

Submitted 5 August, 2011; originally announced August 2011.

Comments: Twelve pages, 4 figures. Accepted for publication in the SciDAC 2011 Proceedings

Measurement of heavy-quark jet photoproduction at HERA

Authors: ZEUS Collaboration, H. Abramowicz, I. Abt, L. Adamczyk, M. Adamus, R. Aggarwal, S. Antonelli, P. Antonioli, A. Antonov, M. Arneodo, V. Aushev, Y. Aushev, O. Bachynska, A. Bamberger, A. N. Barakbaev, G. Barbagli, G. Bari, F. Barreiro, N. Bartosik, D. Bartsch, M. Basile, O. Behnke, J. Behr, U. Behrens, L. Bellagamba , et al. (287 additional authors not shown)

Abstract: Photoproduction of beauty and charm quarks in events with at least two jets has been measured with the ZEUS detector at HERA using an integrated luminosity of 133 $pb^{-1}$. The fractions of jets containing b and c quarks were extracted using the invariant mass of charged tracks associated with secondary vertices and the decay-length significance of these vertices. Differential cross sections as a… ▽ More Photoproduction of beauty and charm quarks in events with at least two jets has been measured with the ZEUS detector at HERA using an integrated luminosity of 133 $pb^{-1}$. The fractions of jets containing b and c quarks were extracted using the invariant mass of charged tracks associated with secondary vertices and the decay-length significance of these vertices. Differential cross sections as a function of jet transverse momentum, $p_{T}^{\text{jet}}$, and pseudorapidity, $ηいーた^{\text{jet}}$, were measured. The data are compared with previous measurements and are well described by next-to-leading-order QCD predictions. △ Less

Submitted 28 April, 2011; originally announced April 2011.

Comments: 35 pages, 9 figures, 3 tables

Report number: DESY-11-067

Journal ref: Eur.Phys.J.C71:1659,2011

Nonthermal Radiation from Type Ia Supernova Remnants

Authors: Paul P. Edmon, Hyesung Kang, T. W. Jones, Renyi Ma

Abstract: We present calculations of expected continuum emissions from Sedov-Taylor phase Type Ia supernova remnants (SNRs), using the energy spectra of cosmic ray (CR) electrons and protons from nonlinear diffusive shock acceleration (DSA) simulations. A new, general-purpose radiative process code, Cosmicp, was employed to calculate the radiation expected from CR electrons and protons and their secondary p… ▽ More We present calculations of expected continuum emissions from Sedov-Taylor phase Type Ia supernova remnants (SNRs), using the energy spectra of cosmic ray (CR) electrons and protons from nonlinear diffusive shock acceleration (DSA) simulations. A new, general-purpose radiative process code, Cosmicp, was employed to calculate the radiation expected from CR electrons and protons and their secondary products. These radio, X-ray and gamma-ray emissions are generally consistent with current observations of Type Ia SNRs. The emissions from electrons in these models dominate the radio through X-ray bands. Decays of πぱい^0 s from p-p collisions mostly dominate the gamma-ray range, although for a hot, low density ISM case (n_{ISM}=0.003 cm^{-3}), the pion decay contribution is reduced sufficiently to reveal the inverse Compton contribution to TeV gamma-rays. In addition, we present simple scalings for the contributing emission processes to allow a crude exploration of model parameter space, enabling these results to be used more broadly. We also discuss the radial surface brightness profiles expected for these model SNRs in the X-ray and gamma-ray bands. △ Less

Submitted 12 May, 2011; v1 submitted 4 March, 2011; originally announced March 2011.

Comments: 37 pages, 7 figures, accepted in MNRAS

Synthetic Observations of Simulated AGN Jets: X-ray Cavities

Authors: P. J. Mendygral, S. M. O'Neill, T. W. Jones

Abstract: Observations of X-ray cavities formed by powerful jets from AGN in galaxy cluster cores are widely used to estimate the energy output of the AGN. Using methods commonly applied to observations of clusters, we conduct synthetic X-ray observations of 3D MHD simulated jet-ICM interactions to test the reliability of measuring X-ray cavity power. These measurements are derived from empirical estimates… ▽ More Observations of X-ray cavities formed by powerful jets from AGN in galaxy cluster cores are widely used to estimate the energy output of the AGN. Using methods commonly applied to observations of clusters, we conduct synthetic X-ray observations of 3D MHD simulated jet-ICM interactions to test the reliability of measuring X-ray cavity power. These measurements are derived from empirical estimates of the enthalpy content of the cavities and their implicit ages. We explore how such physical factors as jet intermittency and observational conditions such as orientation of the jets with respect to the line of sight impact the reliability of observational measurements of cavity enthalpy and age. An estimate of the errors in these quantities can be made by directly comparing "observationally" derived values with "actual" values from the simulations. In our tests, cavity enthalpy derived from observations was typically within a factor of two of the simulation values. Cavity age and, therefore, cavity power are sensitive to the accuracy of the estimated inclination angle of the jets. Cavity age and power estimates within a factor of two of the actual values are possible given an accurate inclination angle. △ Less

Submitted 7 February, 2011; originally announced February 2011.

Comments: 34 pages, 14 figures, accepted for publication in ApJ

Cluster Turbulence: Simulation Insights

Authors: T. W. Jones, David H. Porter, Dongsu Ryu, Jungyeon Cho

Abstract: Cluster media are dynamical, not static; observational evidence suggests they are turbulent. High-resolution simulations of the intracluster media (ICMs) and of idealized, similar media help us understand the complex physics and astrophysics involved. We present a brief overview of the physics behind ICM turbulence and outline the processes that control its development. High-resolution, compressib… ▽ More Cluster media are dynamical, not static; observational evidence suggests they are turbulent. High-resolution simulations of the intracluster media (ICMs) and of idealized, similar media help us understand the complex physics and astrophysics involved. We present a brief overview of the physics behind ICM turbulence and outline the processes that control its development. High-resolution, compressible, isothermal MHD simulations are used to illustrate important dynamical properties of turbulence that develops in media with initially very weak magnetic fields. The simulations follow the growth of magnetic fields and reproduce the characteristics of turbulence. These results are also compared with full cluster simulations that have examined the properties of ICM turbulence. △ Less

Submitted 20 January, 2011; originally announced January 2011.

Comments: 7 pages, 2 figures. To appear the proceedings of the conference "Non-thermal Phenomena in Colliding Galaxy Clusters" (Nice, 15-18 November, 2010), Memorie della Societa Astronomica Italiana, vol 82, n 3

Measurement of beauty production in deep inelastic scattering at HERA using decays into electrons

Authors: ZEUS collaboration, H. Abramowicz, I. Abt, L. Adamczyk, M. Adamus, R. Aggarwal, S. Antonelli, P. Antonioli, A. Antonov, M. Arneodo, V. Aushev, Y. Aushev, O. Bachynska, A. Bamberger, A. N. Barakbaev, G. Barbagli, G. Bari, F. Barreiro, N. Bartosik, D. Bartsch, M. Basile, O. Behnke, J. Behr, U. Behrens, L. Bellagamba , et al. (289 additional authors not shown)

Abstract: The production of beauty quarks in ep interactions has been studied with the ZEUS detector at HERA for exchanged four-momentum squared Q^2 > 10 GeV^2, using an integrated luminosity of 363 pb^{-1}. The beauty events were identified using electrons from semileptonic b decays with a transverse momentum 0.9 < p_T^e < 8 GeV and pseudorapidity |eta^e| < 1.5. Cross sections for beauty production were me… ▽ More The production of beauty quarks in ep interactions has been studied with the ZEUS detector at HERA for exchanged four-momentum squared Q^2 > 10 GeV^2, using an integrated luminosity of 363 pb^{-1}. The beauty events were identified using electrons from semileptonic b decays with a transverse momentum 0.9 < p_T^e < 8 GeV and pseudorapidity |eta^e| < 1.5. Cross sections for beauty production were measured and compared with next-to-leading-order QCD calculations. The beauty contribution to the proton structure function F_2 was extracted from the double-differential cross section as a function of Bjorken-x and Q^2. △ Less

Submitted 10 March, 2011; v1 submitted 19 January, 2011; originally announced January 2011.

Comments: 34 pages, 8 figures, 4 tables. Version published in EPJ C

Report number: DESY-11-005

Journal ref: Eur.Phys.J.C71:1573,2011

Measurement of beauty production in DIS and F_2^bbbar extraction at ZEUS

Authors: ZEUS collaboration, H. Abramowicz, I. Abt, L. Adamczyk, M. Adamus, R. Aggarwal, S. Antonelli, P. Antonioli, A. Antonov, M. Arneodo, V. Aushev, Y. Aushev, O. Bachynska, A. Bamberger, A. N. Barakbaev, G. Barbagli, G. Bari, F. Barreiro, D. Bartsch, M. Basile, O. Behnke, J. Behr, U. Behrens, L. Bellagamba, A. Bertolin , et al. (289 additional authors not shown)

Abstract: Beauty production in deep inelastic scattering with events in which a muon and a jet are observed in the final state has been measured with the ZEUS detector at HERA using an integrated luminosity of 114 pb^-1. The fraction of events with beauty quarks in the data was determined using the distribution of the transverse momentum of the muon relative to the jet. The cross section for beauty producti… ▽ More Beauty production in deep inelastic scattering with events in which a muon and a jet are observed in the final state has been measured with the ZEUS detector at HERA using an integrated luminosity of 114 pb^-1. The fraction of events with beauty quarks in the data was determined using the distribution of the transverse momentum of the muon relative to the jet. The cross section for beauty production was measured in the kinematic range of photon virtuality, Q^2 > 2 Gev^2, and inelasticity, 0.05 < y < 0.7, with the requirement of a muon and a jet. Total and differential cross sections are presented and compared to QCD predictions. The beauty contribution to the structure function F_2 was extracted and is compared to theoretical predictions. △ Less

Submitted 19 May, 2010; originally announced May 2010.

Comments: 24 pages, 5 figures

Report number: DESY-10-047

Journal ref: Eur.Phys.J.C69:347-360,2010

Comparison of Different Methods for Nonlinear Diffusive Shock Acceleration

Authors: D. Caprioli, Hyesung Kang, A. Vladimirov, T. W. Jones

Abstract: We provide a both qualitative and quantitative comparison among different approaches aimed to solve the problem of non-linear diffusive acceleration of particles at shocks. In particular, we show that state-of-the-art models (numerical, Monte Carlo and semi-analytical), even if based on different physical assumptions and implementations, for typical environmental parameters lead to very consistent… ▽ More We provide a both qualitative and quantitative comparison among different approaches aimed to solve the problem of non-linear diffusive acceleration of particles at shocks. In particular, we show that state-of-the-art models (numerical, Monte Carlo and semi-analytical), even if based on different physical assumptions and implementations, for typical environmental parameters lead to very consistent results in terms of shock hydrodynamics, cosmic ray spectrum and also escaping flux spectrum and anisotropy. Strong points and limits of each approach are also discussed, as a function of the problem one wants to study. △ Less

Submitted 27 July, 2010; v1 submitted 12 May, 2010; originally announced May 2010.

Comments: 26 pages, 4 figures, published version (references updated)

Radio Galaxy NGC 1265 unveils the Accretion Shock onto the Perseus Galaxy Cluster

Authors: Christoph Pfrommer, Tom W. Jones

Abstract: We present a consistent 3D model for the head-tail radio galaxy NGC 1265 that explains the complex radio morphology and spectrum by a past passage of the galaxy and radio bubble through a shock wave. Using analytical solutions to the full Riemann problem and hydrodynamical simulations, we study how this passage transformed the plasma bubble into a toroidal vortex ring. Adiabatic compression of the… ▽ More We present a consistent 3D model for the head-tail radio galaxy NGC 1265 that explains the complex radio morphology and spectrum by a past passage of the galaxy and radio bubble through a shock wave. Using analytical solutions to the full Riemann problem and hydrodynamical simulations, we study how this passage transformed the plasma bubble into a toroidal vortex ring. Adiabatic compression of the aged electron population causes it to be energized and to emit low-surface brightness and steep-spectrum radio emission. The large infall velocity of NGC 1265 and the low Faraday rotation measure values and variance of the jet strongly argue that this transformation was due to the accretion shock onto Perseus situated roughly at R_200. Estimating the volume change of the radio bubble enables inferring a shock Mach number of M = 4.2_{-1.2}^{+0.8}, a density jump of 3.4_{-0.4}^{+0.2}, a temperature jump of 6.3_{-2.7}^{+2.5}, and a pressure jump of 21.5 +/- 10.5 while allowing for uncertainties in the equation of state of the radio plasma and volume of the torus. Extrapolating X-ray profiles, we obtain upper limits on the gas temperature and density in the infalling warm-hot intergalactic medium of kT < 0.4 keV and n < 5e-5 / cm^3. The orientation of the ellipsoidally shaped radio torus in combination with the direction of the galaxy's head and tail in the plane of the sky is impossible to reconcile with projection effects. Instead, this argues for post-shock shear flows that have been caused by curvature in the shock surface with a characteristic radius of 850 kpc. The energy density of the shear flow corresponds to a turbulent-to-thermal energy density of 14%. The shock-injected vorticity might be important in generating and amplifying magnetic fields in galaxy clusters. Future LOFAR observations of head-tail galaxies can be complementary probes of accretion shocks onto galaxy clusters. △ Less

Submitted 8 February, 2011; v1 submitted 20 April, 2010; originally announced April 2010.

Comments: 14 pages, 4 figures, ApJ, in print; v3: typos corrected to match the published version; v2: improved presentation, added 2D numerical simulations and exact solution to the 1D Riemann problem of a shock overrunning a spherical bubble that gets transformed into a vortex ring

Three-Dimensional Simulations of Bi-Directed Magnetohydrodynamic Jets Interacting with Cluster Environments

Authors: S. M. O'Neill, T. W. Jones

Abstract: We report on a series of three-dimensional magnetohydrodynamic simulations of active galactic nucleus (AGN) jet propagation in realistic models of magnetized galaxy clusters. We are primarily interested in the details of energy transfer between jets and the intracluster medium (ICM) to help clarify what role such flows could have in the reheating of cluster cores. Our simulated jets feature a ra… ▽ More We report on a series of three-dimensional magnetohydrodynamic simulations of active galactic nucleus (AGN) jet propagation in realistic models of magnetized galaxy clusters. We are primarily interested in the details of energy transfer between jets and the intracluster medium (ICM) to help clarify what role such flows could have in the reheating of cluster cores. Our simulated jets feature a range of intermittency behaviors, including intermittent jets that periodically switch on and off and one model jet that shuts down completely, naturally creating a relic plume. The ICM into which these jets propagate incorporates tangled magnetic field geometries and density substructure designed to mimic some likely features of real galaxy clusters. We find that our jets are characteristically at least 60% efficient at transferring thermal energy to the ICM. Irreversible heat energy is not uniformly distributed, however, instead residing preferentially in regions very near the jet/cocoon boundaries. While intermittency affects the details of how, when, and where this energy is deposited, all of our models generically fail to heat the cluster cores uniformly. Both the detailed density structure and nominally weak magnetic fields in the ICM play interesting roles in perturbing the flows, particularly when the jets are non-steady. Still, this perturbation is never sufficient to isotropize the jet energy deposition, suggesting that some other ingredient is required for AGN jets to successfully reheat cluster cores. △ Less

Submitted 11 January, 2010; originally announced January 2010.

Comments: 19 pages, 18 figures, Accepted for publication in the Astrophysical Journal

Testing Observational Techniques with 3D MHD Jets in Clusters

Authors: Peter J. Mendygral, Sean M. O'Neill, Tom W. Jones

Abstract: Observations of X-ray cavities formed by powerful jets from AGN in galaxy cluster cores are commonly used to estimate the mechanical luminosity of these sources. We test the reliability of observationally measuring this power with synthetic X-ray observations of 3-D MHD simulations of jets in a galaxy cluster environment. We address the role that factors such as jet intermittency and orientation… ▽ More Observations of X-ray cavities formed by powerful jets from AGN in galaxy cluster cores are commonly used to estimate the mechanical luminosity of these sources. We test the reliability of observationally measuring this power with synthetic X-ray observations of 3-D MHD simulations of jets in a galaxy cluster environment. We address the role that factors such as jet intermittency and orientation of the jets on the sky have on the reliability of observational measurements of cavity enthalpy and age. An estimate of the errors in these quantities can be made by directly comparing ``observationally'' derived values with values from the simulations. In our tests, cavity enthalpy, age and mechanical luminosity derived from observations are within a factor of two of the simulation values. △ Less

Submitted 4 September, 2009; v1 submitted 3 September, 2009; originally announced September 2009.

Comments: 4 pages, 3 figures; to appear in proceedings of The Monster's Fiery Breath: Feedback in Galaxies, Groups, and Clusters (AIP conference series)

Turbulence-induced magnetic fields and the structure of Cosmic Ray modified shocks

Authors: A. Beresnyak, T. W. Jones, A. Lazarian

Abstract: We propose a model for Diffusive Shock Acceleration (DSA) in which stochastic magnetic fields in the shock precursor are generated through purely fluid mechanisms of a so-called small-scale dynamo. This contrasts with previous DSA models that considered magnetic fields amplified through cosmic ray streaming instabilities; i.e., either by way of individual particles resonant scattering in the mag… ▽ More We propose a model for Diffusive Shock Acceleration (DSA) in which stochastic magnetic fields in the shock precursor are generated through purely fluid mechanisms of a so-called small-scale dynamo. This contrasts with previous DSA models that considered magnetic fields amplified through cosmic ray streaming instabilities; i.e., either by way of individual particles resonant scattering in the magnetic fields, or by macroscopic electric currents associated with large-scale cosmic ray streaming. Instead, in our picture, the solenoidal velocity perturbations that are required for the dynamo to work are produced through the interactions of the pressure gradient of the cosmic ray precursor and density perturbations in the inflowing fluid. Our estimates show that this mechanism provides fast growth of magnetic field and is very generic. We argue that for supernovae shocks the mechanism is capable of generating upstream magnetic fields that are sufficiently strong for accelerating cosmic rays up to around 10^16 eV. No action of any other mechanism is necessary. △ Less

Submitted 31 October, 2009; v1 submitted 19 August, 2009; originally announced August 2009.

Comments: 10 pages, 5 figures, ApJ accepted

Journal ref: ApJ, 707, 1541 (2009)