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.

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.

We test the DARKexp model for relaxed, self-gravitating, collisionless systems against equilibrium dark matter halos from the Millennium-II simulation. While limited tests of DARKexp against simulations and observations have been carried out elsewhere, this is the first time the testing is done with a large sample of simulated halos spanning a factor of ~ 50 in mass, and using independent fits to density and energy distributions. We show that DARKexp, a one shape parameter family, provides very good fits to the shapes of density profiles, \rho(r), and differential energy distributions, N(E), of individual simulated halos. The best fit shape parameter $\phi_{0}$ obtained from the two types of fits are correlated, though with scatter. Our most important conclusions come from \rho(r) and N(E) that have been averaged over many halos. These show that the bulk of the deviations between DARKexp and individual Millennium-II halos come from halo-to-halo fluctuations, likely driven by substructure, and other density perturbations. The average \rho(r) and N(E) are quite smooth and follow DARKexp very closely. The only deviation that remains after averaging is small, and located at most bound energies for N(E) and smallest radii for \rho(r). Since the deviation is confined to 3-4 smoothing lengths, and is larger for low mass halos, it is likely due to numerical resolution effects.