We perform a suite of two- and three-dimensional magnetohydrodynamic (MHD) simulations with the Athena code of the non-driven Kelvin-Helmholtz instability in the subsonic, weak magnetic field limit. Focusing the analysis on the non-linear turbulent regime, we quantify energy transfer on a scale-by-scale basis and identify the physical mechanisms responsible for energy exchange by developing the diagnostic known as spectral energy transfer function analysis. At late times when the fluid is in a state of MHD turbulence, magnetic tension mediates the dominant mode of energy injection into the magnetic reservoir, whereby turbulent fluid motions twist and stretch the magnetic field lines. This generated magnetic energy turbulently cascades to smaller scales, while being exchanged backwards and forwards with the kinetic energy reservoir, until finally being dissipated. Incorporating explicit dissipation pushes the dissipation scale to larger scales than if the dissipation were entirely numerical. For scales larger than the dissipation scale, we show that the physics of energy transfer in decaying MHD turbulence is robust to numerical effects.

We present the results of a numerical investigation of current-driven instability in magnetized jets. Utilizing the well-tested, relativistic magnetohydrodynamic code Athena, we construct an ensemble of local, co-moving plasma columns in which initial radial force balance is achieved through various combinations of magnetic, pressure, and rotational forces. We then examine the resulting flow morphologies and energetics to determine the degree to which these systems become disrupted, the amount of kinetic energy amplification attained, and the non-linear saturation behaviors. Our most significant finding is that the details of initial force balance have a pronounced effect on the resulting flow morphology. Models in which the initial magnetic field is force-free deform, but do not become disrupted. Systems that achieve initial equilibrium by balancing pressure gradients and/or rotation against magnetic forces, however, tend to shred, mix, and develop turbulence. In all cases, the linear growth of current-driven instabilities is well-represented by analytic models. CDI-driven kinetic energy amplification is slower and saturates at a lower value in force-free models than in those that feature pressure gradients and/or rotation. In rotating columns, we find that magnetized regions undergoing rotational shear are driven toward equipartition between kinetic and magnetic energies. We show that these results are applicable for a large variety of physical parameters, but we caution that algorithmic decisions (such as choice of Riemann solver) can affect the evolution of these systems more than physically motivated parameters.

We have identified the presence of large-scale, low-frequency dynamo cycles in a long-duration, global, magnetohydrodynamic (MHD) simulation of black hole accretion. Such cycles had been seen previously in local shearing box simulations, but we discuss their evolution over 1,500 inner disk orbits of a global pi/4 disk wedge spanning two orders of magnitude in radius and seven scale heights in elevation above/below the disk midplane. The observed cycles manifest themselves as oscillations in azimuthal magnetic field occupying a region that extends into a low-density corona several scale heights above the disk. The cycle frequencies are ten to twenty times lower than the local orbital frequency, making them potentially interesting sources of low-frequency variability when scaled to real astrophysical systems. Furthermore, power spectra derived from the full time series reveal that the cycles manifest themselves at discrete, narrow-band frequencies that often share power across broad radial ranges. We explore possible connections between these simulated cycles and observed low-frequency quasi-periodic oscillations (LFQPOs) in galactic black hole binary systems, finding that dynamo cycles have the appropriate frequencies and are located in a spatial region associated with X-ray emission in real systems. Derived observational proxies, however, fail to feature peaks with RMS amplitudes comparable to LFQPO observations, suggesting that further theoretical work and more sophisticated simulations will be required to form a complete theory of dynamo-driven LFQPOs. Nonetheless, this work clearly illustrates that global MHD dynamos exhibit quasi-periodic behavior on timescales much longer than those derived from test particle considerations.

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.

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.

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 association of an electromagnetic signal with the merger of a pair of supermassive black holes would have many important implications. For example, it would provide new information about gas and magnetic field interactions in dynamical spacetimes as well as a combination of redshift and luminosity distance that would enable precise cosmological tests. A proposal first made by Bode & Phinney (2007) is that because radiation of gravitational waves during the final inspiral and merger of the holes is abrupt and decreases the mass of the central object by a few percent, there will be waves in the disk that can steepen into shocks and thus increase the disk luminosity in a characteristic way. We evaluate this process analytically and numerically. We find that shocks only occur when the fractional mass loss exceeds the half-thickness (h/r) of the disk, hence significant energy release only occurs for geometrically thin disks which are thus at low Eddington ratios. This strongly limits the effective energy release, and in fact our simulations show that the natural variations in disk luminosity are likely to obscure this effect entirely. However, we demonstrate that the reduction of luminosity caused by the retreat of the inner edge of the disk following mass loss is potentially detectable. This decrease occurs even if the disk is geometrically thick, and lasts for a duration on the order of the viscous time of the modified disk. Observationally, the best prospect for detection would be a sensitive future X-ray instrument with a field of view of on the order of a square degree, or possibly a wide-field radio array such as the Square Kilometer Array, if the disk changes produce or interrupt radio emission from a jet.

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.

We examine the evolution and influence of viscosity-induced diskoseismic modes in simulated black hole accretion disks. Understanding the origin and behavior of such oscillations will help us to evaluate their potential role in producing astronomically observed high-frequency quasi-periodic oscillations in accreting black hole binary systems. Our simulated disks are geometrically-thin with a constant half-thickness of five percent the radius of the innermost stable circular orbit. A pseudo-Newtonian potential reproduces the relevant effects of general relativity, and an alpha-model viscosity achieves angular momentum transport and the coupling of orthogonal velocity components in an otherwise ideal hydrodynamic numerical treatment. We find that our simulated viscous disks characteristically develop and maintain trapped global mode oscillations with properties similar to those expected of trapped g-modes and inner p-modes in a narrow range of frequencies just below the maximum radial epicyclic frequency. Although the modes are driven in the inner portion of the disk, they generate waves that propagate at the trapped mode frequencies out to larger disk radii. This finding is contrasted with the results of global magnetohydrodynamic disk simulations, in which such oscillations are not easily identified. Such examples underscore fundamental physical differences between accretion systems driven by the magneto-rotational instability and those for which alpha viscosity serves as a proxy for the physical processes that drive accretion, and we explore potential approaches to the search for diskoseismic modes in full magnetohydrodynamic disks.

We present a set of high-resolution 3D MHD simulations of steady light, supersonic jets, exploring the influence of jet Mach number and the ambient medium on jet propagation and energy deposition over long distances. The results are compared to simple self-similar scaling relations for the morphological evolution of jet-driven structures and to previously published 2D simulations. For this study we simulated the propagation of light jets with internal Mach numbers 3 and 12 to lengths exceeding 100 initial jet radii in both uniform and stratified atmospheres. The propagating jets asymptotically deposit approximately half of their energy flux as thermal energy in the ambient atmosphere, almost independent of jet Mach number or the external density gradient. Nearly one-quarter of the jet total energy flux goes directly into dissipative heating of the ICM, supporting arguments for effective feedback from AGNs to cluster media. The remaining energy resides primarily in the jet and cocoon structures. Despite having different shock distributions and magnetic field features, global trends in energy flow are similar among the different models. As expected the jets advance more rapidly through stratified atmospheres than uniform environments. The asymptotic head velocity in King-type atmospheres shows little or no deceleration. This contrasts with jets in uniform media with heads that are slowed as they propagate. This suggests that the energy deposited by jets of a given length and power depends strongly on the structure of the ambient medium. While our low-Mach jets are more easily disrupted, their cocoons obey evolutionary scaling relations similar to the high-Mach jets.

We consider the effects of the Galactic magnetic field on the propagation of ultra high energy cosmic rays (UHECRs). By employing two methods of trajectory simulation, we investigate the possibility that UHECRs are produced within the Galaxy with paths strongly influenced by the Galactic magnetic field. Such trajectories have the potential to reconcile the existing conflict between proposed local sources and isotropic UHECR arrival directions.