BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

T. W. Jones

T. W. Jones appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.HEastro-phastro-ph.GAastro-ph.SRphysics.flu-dyn

Catalog footprint

What is connected

36works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

36 published item(s)

preprint2020arXiv

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

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}$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Shocked Narrow-Angle Tail Radio Galaxies: Simulations and Emissions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Turbulence and Vorticity in Galaxy Clusters Generated by Structure Formation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmic rays in galaxy clusters and their non-thermal emission

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Diffusive Shock Acceleration Simulations of Radio Relics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Inclusive-jet photoproduction at HERA and determination of alphas

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Search for single-top production in ep collisions at HERA

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$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cluster Turbulence: Simulation Insights

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Exclusive electroproduction of two pions at HERA

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 π+π-.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measurement of heavy-quark jet photoproduction at HERA

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

MHD Turbulence Simulation in a Cosmic Structure Context

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nonthermal Radiation from Type Ia Supernova Remnants

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Synthetic Observations of Simulated AGN Jets: X-ray Cavities

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Comparison of Different Methods for Nonlinear Diffusive Shock Acceleration

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Self-Similar Evolution of Cosmic-Ray Modified Shocks: The Cosmic-Ray Spectrum

We use kinetic simulations of diffusive shock acceleration (DSA) to study the time-dependent evolution of plane, quasi-parallel, cosmic-ray (CR) modified shocks. Thermal leakage injection of low energy CRs and finite Alfvén wave propagation and dissipation are included. Bohm diffusion as well as the diffusion with the power-law momentum dependence are modeled. As long as the acceleration time scale to relativistic energies is much shorter than the dynamical evolution time scale of the shocks, the precursor and subshock transition approach the time-asymptotic state, which depends on the shock sonic and Alfvénic Mach numbers and the CR injection efficiency. For the diffusion models we employ, the shock precursor structure evolves in an approximately self-similar fashion, depending only on the similarity variable, x/(u_s t). During this self-similar stage, the CR distribution at the subshock maintains a characteristic form as it evolves: the sum of two power-laws with the slopes determined by the subshock and total compression ratios with an exponential cutoff at the highest accelerated momentum, p_{max}(t). Based on the results of the DSA simulations spanning a range of Mach numbers, we suggest functional forms for the shock structure parameters, from which the aforementioned form of CR spectrum can be constructed. These analytic forms may represent approximate solutions to the DSA problem for astrophysical shocks during the self-similar evolutionary stage as well as during the steady-state stage if p_{max} is fixed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Three-dimensional Magnetohydrodynamic Simulations of Buoyant Bubbles in Galaxy Clusters

We report results of 3D MHD simulations of the dynamics of buoyant bubbles in magnetized galaxy cluster media. The simulations are three dimensional extensions of two dimensional calculations reported by Jones & De Young (2005). Initially spherical bubbles and briefly inflated spherical bubbles all with radii a few times smaller than the intracluster medium (ICM) scale height were followed as they rose through several ICM scale heights. Such bubbles quickly evolve into a toroidal form that, in the absence of magnetic influences, is stable against fragmentation in our simulations. This ring formation results from (commonly used) initial conditions that cause ICM material below the bubbles to drive upwards through the bubble, creating a vortex ring; that is, hydrostatic bubbles develop into "smoke rings", if they are initially not very much smaller or very much larger than the ICM scale height. Even modest ICM magnetic fields with beta = P_gas/P_mag ~ 10^3 can influence the dynamics of the bubbles, provided the fields are not tangled on scales comparable to or smaller than the size of the bubbles. Quasi-uniform, horizontal fields with initial beta ~ 10^2 bifurcated our bubbles before they rose more than about a scale height of the ICM, and substantially weaker fields produced clear distortions. On the other hand, tangled magnetic fields with similar, modest strengths are generally less easily amplified by the bubble motions and are thus less influential in bubble evolution. Inclusion of a comparably strong, tangled magnetic field inside the initial bubbles had little effect on our bubble evolution, since those fields were quickly diminished through expansion of the bubble and reconnection of the initial field.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Acceleration of Cosmic Rays at Large Scale Cosmic Shocks in the Universe

Cosmological hydrodynamic simulations of large scale structure in the universe have shown that accretion shocks and merger shocks form due to flow motions associated with the gravitational collapse of nonlinear structures. Estimated speed and curvature radius of these shocks could be as large as a few 1000 km/s and several Mpc, respectively. According to the diffusive shock acceleration theory, populations of cosmic-ray particles can be injected and accelerated to very high energy by astrophysical shocks in tenuous plasmas. In order to explore the cosmic ray acceleration at the cosmic shocks, we have performed nonlinear numerical simulations of cosmic ray (CR) modified shocks with the newly developed CRASH (Cosmic Ray Amr SHock) numerical code. We adopted the Bohm diffusion model for CRs, based on the hypothesis that strong Alfvén waves are self-generated by streaming CRs. The shock formation simulation includes a plasma-physics-based ``injection'' model that transfers a small proportion of the thermal proton flux through the shock into low energy CRs for acceleration there. We found that, for strong accretion shocks, CRs can absorb most of shock kinetic energy and the accretion shock speed is reduced up to 20 %, compared to pure gas dynamic shocks. For merger shocks with small Mach numbers, however, the energy transfer to CRs is only about 10-20 % with an associated CR particle fraction of $10^{-3}$. Nonlinear feedback due to the CR pressure is insignificant in the latter shocks. Although detailed results depend on models for the particle diffusion and injection, these calculations show that cosmic shocks in large scale structure could provide acceleration sites of extragalactic cosmic rays of the highest energy.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

On the Exchange of Kinetic and Magnetic Energy Between Clouds and the Interstellar Medium

We investigate, through 2D MHD numerical simulations, the interaction of a uniform magnetic field oblique to a moving interstellar cloud. In particular we explore the transformation of cloud kinetic energy into magnetic energy as a result of field line stretching. Some previous simulations have emphasized the possible dynamical importance of a ``magnetic shield'' formed around clouds when the magnetic field is perpendicular to the cloud motion (Jones et al. 1996, Miniati et al. 1998). It was not clear, however, how dependent those findings were to the assumed field configuration and cloud properties. To expand our understanding of this effect, we examine several new cases by varing the magnetic field orientation angle with respect to the cloud motion (θ), the cloud-background density contrast, and the cloud Mach number. We show that in 2D and with θlarge enough, the magnetic field tension can become dominant in the dynamics of the motion of high density contrast, low Mach number clouds. In such cases a significant fraction of cloud kinetic energy can be transformed into magnetic energy with the magnetic pressure at the cloud nose exceeding the ram pressure of the impinging flow. We derive a characteristic timescale for this process of energy ``conversion''. We find also that unless the cloud motion is highly aligned to the magnetic field, reconnection through tearing mode instabilities in the cloud wake limit the formation of a strong flux rope feature following the cloud. Finally we attempt to interpret some observational properties of the magnetic field in view of our results.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

The MHD Kelvin-Helmholtz Instability III: The Role of Sheared Magnetic Field in Planar Flows

We have carried out simulations of the nonlinear evolution of the magnetohydrodynamic (MHD) Kelvin-Helmholtz (KH) instability for compressible fluids in $2\frac{1}{2}$-dimensions, extending our previous work by Frank et al (1996) and Jones \etal (1997). In the present work we have simulated flows in the x-y plane in which a ``sheared'' magnetic field of uniform strength ``smoothly'' rotates across a thin velocity shear layer from the z direction to the x direction, aligned with the flow field. We focus on dynamical evolution of fluid features, kinetic energy dissipation, and mixing of the fluid between the two layers, considering their dependence on magnetic field strength for this geometry. The introduction of magnetic shear can allow a Cat's Eye-like vortex to form, even when the field is stronger than the nominal linear instability limit given above. For strong fields that vortex is asymmetric with respect to the preliminary shear layer, however, so the subsequent dissipation is enhanced over the uniform field cases of comparable field strength. In fact, so long as the magnetic field achieves some level of dynamical importance during an eddy turnover time, the asymmetries introduced through the magnetic shear will increase flow complexity, and, with that, dissipation and mixing. The degree of the fluid mixing between the two layers is strongly influenced by the magnetic field strength. Mixing of the fluid is most effective when the vortex is disrupted by magnetic tension during transient reconnection, through local chaotic behavior that follows.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Hydrodynamics of Cloud Collisions in 2D: The Fate of Clouds in a Multi-phase Medium

We have studied head-on collisions between equal-mass, mildly supersonic (Mach number 1.5) HI clouds, in a standard Two-phase ISM (T_cl = 74 K, n_cl = 22 cm^-3, χ= 100). We explore the role of various factors, including the radiative cooling parameter η= t_rad/t_coll (t_coll=R_c/v_c), evolutionary modifications on the cloud structure (by colliding clouds ``evolved'' through independent motion within the intercloud medium (WIM)), and the symmetry of the problem (by colliding initially identical clouds, evolved to different ages before impact). The presence of bow shocks and ram pressure from material in the cloud wake, developed during such evolution through the WIM, significantly alters these interactions with respect to the standard case of non-evolved clouds. In general, in our adiabatic collisions the clouds are disrupted and convert their gas into a few low density contrast clumps. By contrast, for symmetric radiative cases we find that the two clouds coalesce, with almost all the initial kinetic energy radiated away. On the other hand, for both adiabatic and radiative collisions, asymmetric collisions have a much greater tendency to disrupt the two clouds. Fragmentation of the clouds may occur, and instabilities are in general enhanced. In addition, radiative cooling is less efficient in our asymmetric interactions, so that those parts of the clouds that initially seem to merge are more likely to re-expand and fade into the WIM. Since the majority of real cloud collisions should be asymmetric for one reason or another, we conclude that most gasdynamical diffuse cloud collisions will be disruptive, at least in the absence of significant self-gravity or of a significant magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Diffusive Shock Acceleration Simulations: Comparison with Particle Methods and Bow Shock Measurements

Direct comparisons of diffusive particle acceleration numerical simulations have been made against Monte Carlo and hybrid plasma simulations by Ellison {\it et. al.} (1993) and against observations at the earth's bow shock presented by Ellison {\it et. al.} (1990). Toward this end we have introduced a new numerical scheme for injection of cosmic-ray particles out of the thermal plasma, modeled by way of the diffusive scattering process itself; that is, the diffusion and acceleration across the shock front of particles out of the suprathermal tail of the Maxwellian distribution. We find that all of these computational methods (diffusion-advection, two-fluid, Monte Carlo and hybrid) are in substantial agreement on the issues they can simultaneously address, so that the essential physics of diffusive particle acceleration is adequately contained within each. This is despite the fact that each makes what appear to be very different assumptions or approximations.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for Multi-Dimensional Flow

We present for astrophysical use a multi-dimensional numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference method on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. Multiple spatial dimensions are treated through a Strang-type operator splitting. The constraint of a divergence-free field is enforced exactly by calculating a correction via a gauge transformation in each time step. Results from two-dimensional shock tube tests show that the code captures correctly discontinuities in all three MHD waves families as well as contact discontinuities. The numerical viscosities and resistivity in the code, which are useful in order to understand simulations involving turbulent flows, are estimated through the decay of two-dimensional linear waves. Finally, the robustness of the code in two-dimensions is demonstrated through calculations of the Kelvin-Helmholtz instability and the Orszag-Tang vortex.

0 citations0 reviewsNo code linkedOpen access