Magnetic fields play vital roles in intracluster media (ICMs), but estimating their strengths and distributions from observations is a major challenge. Faraday rotation measures (RMs) are widely applied to this task, so it is critical to understand inherent uncertainties in RM analysis. In this paper, we seek to characterize those uncertainties given the types of information available today, independent of the specific technique used. We conduct synthetic RM observations through the ICM of a galaxy cluster drawn from an MHD cosmological simulation in which the magnetic field is known. We analyze the synthetic RM observations using an analytical formalism based on commonly used model assumptions allowing us to relate model physical variables to outcome uncertainties. Despite the simplicity of some assumptions, and unknown physical parameters, we are able to extract an approximate magnitude of the central magnetic field within an apparently irreducible uncertain factor approximately 3. Principal, largely irreducible, uncertainties come from the unknown depth along the line of sight of embedded polarized sources, the lack of robust coherence lengths from area-constrained polarization sampling, and the unknown scaling between ICM electron density and magnetic field strength. The RM-estimated central magnetic field strengths span more than an order of magnitude including the full range of synthetic experiments.

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.

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.

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.

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.

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.

In galaxy clusters, modern radio interferometers observe non-thermal radio sources with unprecedented spatial and spectral resolution. For the first time, the new data allows to infer the structure of the intra-cluster magnetic fields on small scales via Faraday tomography. This leap forward demands new numerical models for the amplification of magnetic fields in cosmic structure formation - the cosmological magnetic dynamo. Here we present a novel numerical approach to astrophyiscal MHD simulations aimed to resolve this small-scale dynamo in future cosmological simulations. As a first step, we implement a fifth order WENO scheme in the new code WOMBAT. We show that this scheme doubles the effective resolution of the simulation and is thus less expensive than common second order schemes. WOMBAT uses a novel approach to parallelization and load balancing developed in collaboration with performance engineers at Cray Inc. This will allow us scale simulation to the exaflop regime and achieve kpc resolution in future cosmological simulations of galaxy clusters. Here we demonstrate the excellent scaling properties of the code and argue that resolved simulations of the cosmological small scale dynamo within the whole virial radius are possible in the next years.

Deep learning is a promising tool to determine the physical model that describes our universe. To handle the considerable computational cost of this problem, we present CosmoFlow: a highly scalable deep learning application built on top of the TensorFlow framework. CosmoFlow uses efficient implementations of 3D convolution and pooling primitives, together with improvements in threading for many element-wise operations, to improve training performance on Intel(C) Xeon Phi(TM) processors. We also utilize the Cray PE Machine Learning Plugin for efficient scaling to multiple nodes. We demonstrate fully synchronous data-parallel training on 8192 nodes of Cori with 77% parallel efficiency, achieving 3.5 Pflop/s sustained performance. To our knowledge, this is the first large-scale science application of the TensorFlow framework at supercomputer scale with fully-synchronous training. These enhancements enable us to process large 3D dark matter distribution and predict the cosmological parameters $\Omega_M$, $\sigma_8$ and n$_s$ with unprecedented accuracy.

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.

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.

Supernova remnants (SNRs) are believed to accelerate particles up to high energies through the mechanism of diffusive shock acceleration (DSA). Except for direct plasma simulations, all modeling efforts must rely on a given form of the diffusion coefficient, a key parameter that embodies the interactions of energetic charged particles with the magnetic turbulence. The so-called Bohm limit is commonly employed. In this paper we revisit the question of acceleration at perpendicular shocks, by employing a realistic model of perpendicular diffusion. Our coefficient reduces to a power-law in momentum for low momenta (of index $\alpha$), but becomes independent of the particle momentum at high momenta (reaching a constant value $\kappa_{\infty}$ above some characteristic momentum $p_{\rm c}$). We first provide simple analytical expressions of the maximum momentum that can be reached at a given time with this coefficient. Then we perform time-dependent numerical simulations to investigate the shape of the particle distribution that can be obtained when the particle pressure back-reacts on the flow. We observe that, for a given index $\alpha$ and injection level, the shock modifications are similar for different possible values of $p_{\rm c}$, whereas the particle spectra differ markedly. Of particular interest, low values of $p_{\rm c}$ tend to remove the concavity once thought to be typical of non-linear DSA, and result in steep spectra, as required by recent high-energy observations of Galactic SNRs.

We present a pair of 3-d magnetohydrodynamical simulations of intermittent jets from a central active galactic nucleus (AGN) in a galaxy cluster extracted from a high resolution cosmological simulation. The selected cluster was chosen as an apparently relatively relaxed system, not having undergone a major merger in almost 7 Gyr. Despite this characterization and history, the intra-cluster medium (ICM) contains quite active "weather". We explore the effects of this ICM weather on the morphological evolution of the AGN jets and lobes. The orientation of the jets is different in the two simulations so that they probe different aspects of the ICM structure and dynamics. We find that even for this cluster that can be characterized as relaxed by an observational standard, the large-scale, bulk ICM motions can significantly distort the jets and lobes. Synthetic X-ray observations of the simulations show that the jets produce complex cavity systems, while synthetic radio observations reveal bending of the jets and lobes similar to wide-angle tail (WAT) radio sources. The jets are cycled on and off with a 26 Myr period using a 50% duty cycle. This leads to morphological features similar to those in "double-double" radio galaxies. While the jet and ICM magnetic fields are generally too weak in the simulations to play a major role in the dynamics, Maxwell stresses can still become locally significant.

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.

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.

The high-resolution X-ray spectroscopy made possible by the 1999 deployment of the Chandra X-ray Observatory has revolutionized our understanding of stellar X-ray emission. Many puzzles remain, though, particularly regarding the mechanisms of X-ray emission from OB stars. Although numerous individual stars have been observed in high-resolution, realizing the full scientific potential of these observations will necessitate studying the high-resolution Chandra dataset as a whole. To facilitate the rapid comparison and characterization of stellar spectra, we have compiled a uniformly processed database of all stars observed with the Chandra High Energy Transmission Grating (HETG). This database, known as X-Atlas, is accessible through a web interface with searching, data retrieval, and interactive plotting capabilities. For each target, X-Atlas also features predictions of the low-resolution ACIS spectra convolved from the HETG data for comparison with stellar sources in archival ACIS images. Preliminary analyses of the hardness ratios, quantiles, and spectral fits derived from the predicted ACIS spectra reveal systematic differences between the high-mass and low-mass stars in the atlas and offer evidence for at least two distinct classes of high-mass stars. A high degree of X-ray variability is also seen in both high and low-mass stars, including Capella, long thought to exhibit minimal variability. X-Atlas contains over 130 observations of approximately 25 high-mass stars and 40 low-mass stars and will be updated as additional stellar HETG observations become public. The atlas has recently expanded to non-stellar point sources, and Low Energy Transmission Grating (LETG) observations are currently being added as well.