In previous work, we developed tools for quantifying the tidal distortion of a black hole's event horizon due to an orbiting companion. These tools use techniques which require large mass ratios (companion mass $\mu$ much smaller than black hole mass $M$), but can be used for arbitrary bound orbits, and for any black hole spin. We also showed how to visualize these distorted black holes by embedding their horizons in a global Euclidean 3-space, ${\mathbb{E}}^3$. Such visualizations illustrate interesting and important information about horizon dynamics. Unfortunately, we could not visualize black holes with spin parameter $a_* > \sqrt{3}/2 \approx 0.866$: such holes cannot be globally embedded into ${\mathbb{E}}^3$. In this paper, we overcome this difficulty by showing how to embed the horizons of tidally distorted Kerr black holes in a hyperbolic 3-space, ${\mathbb{H}}^3$. We use black hole perturbation theory to compute the Gaussian curvatures of tidally distorted event horizons, from which we build a two-dimensional metric of their distorted horizons. We develop a numerical method for embedding the tidally distorted horizons in ${\mathbb{H}}^3$. As an application, we give a sequence of embeddings into ${\mathbb{H}}^3$ of a tidally interacting black hole with spin $a_*=0.9999$. A small amplitude, high frequency oscillation seen in previous work shows up particularly clearly in these embeddings.

The Bondi-van der Burg-Metzner-Sachs (BMS) group is the asymptotic symmetry group of asymptotically flat spacetime. It is infinite dimensional and entails an infinite number of conservation laws. According to the black hole membrane paradigm, null infinity (in asymptotically flat spacetime) and black hole event horizons behave like fluid membranes. The fluid dynamics of the membrane is governed by an infinite set of symmetries and conservation laws. Our main result is to point out that the infinite set of symmetries and conserved charges of the BMS group and the membrane paradigm are the same. This relationship has several consequences. First, it sheds light on the physical interpretation of BMS conservation laws. Second, it generalizes the BMS conservation laws to arbitrary subregions of arbitrary null surfaces. Third, it clarifies the identification of the superrotation subgroup of the BMS group. We briefly comment on the black hole information problem.

Black hole jet power depends on the angular velocity of magnetic field lines, $\Omega_F$. Force-free black hole magnetospheres typically have $\Omega_F/\Omega_H \approx 0.5$, where $\Omega_H$ is the angular velocity of the horizon. We give a streamlined proof of this result using an extension of the classical black hole membrane paradigm. The proof is based on an impedance-matching argument between membranes at the horizon and infinity. Then we consider a general relativistic magnetohydrodynamic simulation of an accreting, spinning black hole and jet. We find that the theory correctly describes the simulation in the jet region. However, the field lines threading the horizon near the equator have much smaller $\Omega_F/\Omega_H$ because the force-free approximation breaks down in the accretion flow.

Numerical simulations indicate that black holes carrying linear momentum and/or orbital momentum can power jets. The jets extract the kinetic energy stored in the black hole's motion. This could provide an important electromagnetic counterpart to gravitational wave searches. We develop the theory underlying these jets. In particular, we derive the analogues of the Penrose process and the Blandford-Znajek jet power prediction for boosted black holes. The jet power we find is $(v/2M)^2 \Phi^2/(4\pi)$, where $v$ is the hole's velocity, $M$ is its mass, and $\Phi$ is the magnetic flux. We show that energy extraction from boosted black holes is conceptually similar to energy extraction from spinning black holes. However, we highlight two key technical differences: in the boosted case, jet power is no longer defined with respect to a Killing vector, and the relevant notion of black hole mass is observer dependent. We derive a new version of the membrane paradigm in which the membrane lives at infinity rather than the horizon and we show that this is useful for interpreting jets from boosted black holes. Our jet power prediction and the assumptions behind it can be tested with future numerical simulations.

We report on the MIT Epoch of Reionization (MITEoR) experiment, a pathfinder low-frequency radio interferometer whose goal is to test technologies that improve the calibration precision and reduce the cost of the high-sensitivity 3D mapping required for 21 cm cosmology. MITEoR accomplishes this by using massive baseline redundancy, which enables both automated precision calibration and correlator cost reduction. We demonstrate and quantify the power and robustness of redundancy for scalability and precision. We find that the calibration parameters precisely describe the effect of the instrument upon our measurements, allowing us to form a model that is consistent with $\chi^2$ per degree of freedom < 1.2 for as much as 80% of the observations. We use these results to develop an optimal estimator of calibration parameters using Wiener filtering, and explore the question of how often and how finely in frequency visibilities must be reliably measured to solve for calibration coefficients. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious Hydrogen Epoch of Reionization Array (HERA) project and other next-generation instruments, which would incorporate many identical or similar technologies.

Spinning black holes tend to expel magnetic fields. In this way they are similar to superconductors. It has been a persistent concern that this black hole "Meissner effect" could quench jet power at high spins. This would make it impossible for the rapidly rotating black holes in Cyg X-1 and GRS 1915+105 to drive Blandford-Znajek jets. We give a simple geometrical argument why fields which become entirely radial near the horizon are not expelled by the Meissner effect and may continue to power jets up to the extremal limit. A simple and natural example is a split-monopole field. We stress that ordinary Blandford-Znajek jets are impossible if the Meissner effect operates and expels the field. Finally, we note that in our general relativistic magnetohydrodynamic simulations of black hole jets, there is no evidence that jets are quenched by the Meissner effect. The simulated jets develop a large split monopole component spontaneously which supports our proposal for how the Meissner effect is evaded and jets from rapidly rotating black holes are powered in nature.

Extremal black holes tend to expel magnetic and electric fields. Fields are unable to reach the horizon because the length of the black hole throat blows up in the extremal limit. The length of the throat is related to the amount of entanglement between modes on either side of the horizon. So it is natural to try to relate the black hole Meissner effect to entanglement. We derive the black hole Meissner effect directly from the low temperature limit of two-point functions in the Hartle-Hawking vacuum. Then we discuss several new examples of the black hole Meissner effect, its applications to astrophysics, and its relationship to gauge invariance.

Using long-duration general relativistic magnetohydrodynamic simulations of radiatively inefficient accretion discs, the energy, momentum and mass outflow rates from such systems are estimated. Outflows occur via two fairly distinct modes: a relativistic jet and a sub-relativistic wind. The jet power depends strongly on the black hole spin and on the magnetic flux at the horizon. Unless these are very small, the energy output in the jet dominates over that in the wind. For a rapidly spinning black hole accreting in the magnetically arrested limit, it is confirmed that jet power exceeds the total rate of accretion of rest mass energy. However, because of strong collimation, the jet probably does not have a significant feedback effect on its immediate surroundings. The power in the wind is more modest and shows a weaker dependence on black hole spin and magnetic flux. Nevertheless, because the wind subtends a large solid angle, it is expected to provide efficient feedback on a wide range of scales inside the host galaxy. Empirical formulae are obtained for the energy and momentum outflow rates in the jet and the wind.

General relativistic magnetohydrodynamic (GRMHD) simulations are providing influential models for black hole spin measurements, gamma ray bursts, and supermassive black hole feedback. Many of these simulations use the same initial condition: a rotating torus of fluid in hydrostatic equilibrium. A persistent concern is that simulation results sometimes depend on arbitrary features of the initial torus. For example, the Bernoulli parameter (which is related to outflows), appears to be controlled by the Bernoulli parameter of the initial torus. In this paper, we give a new equilibrium torus solution and describe two applications for the future. First, it can be used as a more physical initial condition for GRMHD simulations than earlier torus solutions. Second, it can be used in conjunction with earlier torus solutions to isolate the simulation results that depend on initial conditions. We assume axisymmetry, an ideal gas equation of state, constant entropy, and ignore self-gravity. We fix an angular momentum distribution and solve the relativistic Euler equations in the Kerr metric. The Bernoulli parameter, rotation rate, and geometrical thickness of the torus can be adjusted independently. Our torus tends to be more bound and have a larger radial extent than earlier torus solutions. While this paper was in preparation, several GRMHD simulations appeared based on our equilibrium torus. We believe it will continue to provide a more realistic starting point for future simulations.

If spacetime is built out of quantum bits, does the shape of space depend on how the bits are entangled? The ER=EPR conjecture relates the entanglement entropy of a collection of black holes to the cross sectional area of Einstein-Rosen (ER) bridges (or wormholes) connecting them. We show that the geometrical entropy of classical ER bridges satisfies the subadditivity, triangle, strong subadditivity, and CLW inequalities. These are nontrivial properties of entanglement entropy, so this is evidence for ER=EPR. We further show that the entanglement entropy associated to classical ER bridges has nonpositive interaction information. This is not a property of entanglement entropy, in general. For example, the entangled four qubit pure state |GHZ_4>=(|0000>+|1111>)/\sqrt{2} has positive interaction information, so this state cannot be described by a classical ER bridge. Large black holes with massive amounts of entanglement between them can fail to have a classical ER bridge if they are built out of |GHZ_4> states. States with nonpositive interaction information are called monogamous. We conclude that classical ER bridges require monogamous EPR correlations.

Recently it has been observed that the scaling of jet power with black hole spin in galactic X-ray binaries is consistent with the predictions of the Blandford-Znajek (BZ) jet model. These observations motivate us to revisit the BZ model using general relativistic magnetohydrodynamic simulations of magnetized jets from accreting (h/r ~ 0.3), spinning (0 < a_* < 0.98) black holes. We have three main results. First, we quantify the discrepancies between the BZ jet power and our simulations: assuming maximum efficiency and uniform fields on the horizon leads to a ~10% overestimate of jet power, while ignoring the accretion disk leads to a further ~50% overestimate. Simply reducing the standard BZ jet power prediction by 60% gives a good fit to our simulation data. Our second result is to show that the membrane formulation of the BZ model correctly describes the physics underlying simulated jets: torques, dissipation, and electromagnetic fields on the horizon. This provides intuitive yet rigorous pictures for the black hole energy extraction process. Third, we compute the effective resistance of the load region and show that the load and the black hole achieve near perfect impedance matching. Taken together, these results increase our confidence in the BZ model as the correct description of jets observed from astrophysical black holes.

We briefly summarize the method of simulating Sgr A* polarized sub-mm spectra from the accretion flow and fitting the observed spectrum. The dynamical flow model is based on three-dimensional general relativistic magneto hydrodynamic simulations. Fully self-consistent radiative transfer of polarized cyclo-synchrotron emission is performed. We compile a mean sub-mm spectrum of Sgr A* and fit it with the mean simulated spectra. We estimate the ranges of inclination angle theta=42-75deg, mass accretion rate Mdot=(1.4-7.0)*10^{-8}Msun/yr, and electron temperature Te=(3-4)*10^{10}K at 6M. We discuss multiple caveats in dynamical modeling, which must be resolved to make further progress.

Almost all hydrodynamic accretion disk models parametrize viscosity with the dimensionless parameter alpha. There is no detailed model for alpha, so it is usually taken to be a constant. However, global simulations of magnetohydrodynamic disks find that alpha varies with distance from the central object. Also, Newtonian simulations tend to find smaller alpha's than general relativistic simulations. We seek a one-dimensional model for alpha that can reproduce these two observations. We are guided by data from six general relativistic magnetohydrodynamic accretion disk simulations. The variation of alpha in the inner, laminar regions of the flow results from stretching of mean magnetic field lines by the flow. The variation of alpha in the outer, turbulent regions results from the dependence of the magnetorotational instability on the dimensionless shear rate. We give a one-dimensional prescription for alpha(r) that captures these two effects and reproduces the radial variation of alpha observed in the simulations. For thin disks, the prescription simplifies to the formula alpha(r)=0.025[q(r)/1.5]^6, where the shear parameter, q(r), is an analytical function of radius in the Kerr metric. The coefficient and exponent are inferred from our simulations and will change as better simulation data becomes available. We conclude that the alpha-viscosity prescription can be extended to the radially varying alpha's observed in simulations. It is possible that Newtonian simulations find smaller alpha's than general relativistic simulations because the shear parameter is lower in Newtonian flows.

We present results from two long-duration GRMHD simulations of advection-dominated accretion around a non-spinning black hole. The first simulation was designed to avoid significant accumulation of magnetic flux around the black hole. This simulation was run for a time of 200,000GM/c^3 and achieved inflow equilibrium out to a radius \sim90GM/c^2. Even at this relatively large radius, the mass outflow rate \dot{M}_{out} is found to be only 60% of the net mass inflow rate \dot{M}_{BH} into the black hole. The second simulation was designed to achieve substantial magnetic flux accumulation around the black hole in a magnetically arrested disc. This simulation was run for a shorter time of 100,000GM/c^3. Nevertheless, because the mean radial velocity was several times larger than in the first simulation, it reached inflow equilibrium out to a radius \sim170GM/c^2. Here, \dot{M}_{out} becomes equal to \dot{M}_{BH} at r\sim 160GM/c^2. Since the mass outflow rates in the two simulations do not show robust convergence with time, it is likely that the true outflow rates are lower than our estimates. The effect of black hole spin on mass outflow remains to be explored. Neither simulation shows strong evidence for convection, though a complete analysis including the effect of magnetic fields is left for the future.

We obtain estimates of Sgr A* accretion flow and black hole parameters by fitting polarized sub-mm observations with spectra computed using three-dimensional (3D) general relativistic (GR) magnetohydrodynamical (MHD) (GRMHD) simulations. Observations are compiled from averages over many epochs from reports in 29 papers for estimating the mean fluxes Fnu, linear polarization (LP) fractions, circular polarization (CP) fractions, and electric vector position angles (EVPAs). GRMHD simulations are computed with dimensionless spins a_*=0,0.5,0.7,0.9,0.98 over a 20,000M time interval. We perform fully self-consistent GR polarized radiative transfer using our new code to explore the effects of spin a_*, inclination angle \theta, position angle (PA), accretion rate Mdot, and electron temperature Te (Te is reported for radius 6M). By fitting the mean sub-mm fluxes and LP/CP fractions, we obtain estimates for these model parameters and determine the physical effects that could produce polarization signatures. Our best bet model has a_*=0.5, \theta=75deg, PA=115deg, Mdot=4.6*10^{-8}M_Sun/year, and Te=3.1*10^10K at 6M. The sub-mm CP is mainly produced by Faraday conversion as modified by Faraday rotation, and the emission region size at 230GHz is consistent with the VLBI size of 37microas. Across all spins, model parameters are in the ranges \theta=42deg-75deg, Mdot=(1.4-7.0)*10^{-8}M_Sun/year, and Te=(3-4)*10^10K. Polarization is found both to help differentiate models and to introduce new observational constraints on the effects of the magnetic field that might not be fit by accretion models so-far considered.

It is generally thought that the light coming from the inner plunging region of black hole accretion discs contributes negligibly to the disc's overall spectrum, i.e. the plunging fluid is swallowed by the black hole before it has time to radiate. In the standard disc model used to fit X-ray observations of accretion discs, the plunging region is assumed to be perfectly dark. However, numerical simulations that include the full physics of the magnetized flow predict that a small fraction of the disc's total luminosity emanates from the plunging region. We investigate the observational consequences of this neglected inner light. We compute radiative transfer based disc spectra that correspond to 3D general relativistic magnetohydrodynamic simulated discs (which produce light inside their plunging regions). In the context of black hole spin estimation, we find that the neglected inner light only has a modest effect (this bias is less than typical observational systematic errors). For rapidly spinning black holes, we find that the combined emission from the plunging region produces a weak power-law tail at high energies. This indicates that infalling matter is the origin for some of the `coronal' emission observed in the thermal dominant and steep power-law states of X-ray binaries.

We present an analytical solution for thin disk accretion onto a Kerr black hole that extends the standard Novikov-Thorne alpha-disk in three ways: (i) it incorporates nonzero stresses at the inner edge of the disk, (ii) it extends into the plunging region, and (iii) it uses a corrected vertical gravity formula. The free parameters of the model are unchanged. Nonzero boundary stresses are included by replacing the Novikov-Thorne no torque boundary condition with the less strict requirement that the fluid velocity at the innermost stable circular orbit is the sound speed, which numerical models show to be the correct behavior for luminosities below ~30% Eddington. We assume the disk is thin so we can ignore advection. Boundary stresses scale as alpha*h and advection terms scale as h^2 (where h is the disk opening angle (h=H/r)), so the model is self-consistent when h < alpha. We compare our solution with slim disk models and general relativistic magnetohydrodynamic disk simulations. The model may improve the accuracy of black hole spin measurements.

The X-ray spectra of accretion discs of eight stellar-mass black holes have been analyzed to date using the thermal continuum fitting method, and the spectral fits have been used to estimate the spin parameters of the black holes. However, the underlying model used in this method of estimating spin is the general relativistic thin-disc model of Novikov & Thorne, which is only valid for razor-thin discs. We therefore expect errors in the measured values of spin due to inadequacies in the theoretical model. We investigate this issue by computing spectra of numerically calculated models of thin accretion discs around black holes, obtained via three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. We apply the continuum fitting method to these computed spectra to estimate the black hole spins and check how closely the values match the actual spin used in the GRMHD simulations. We find that the error in the dimensionless spin parameter is up to about 0.2 for a non-spinning black hole, depending on the inclination. For black holes with spins of 0.7, 0.9 and 0.98, the errors are up to about 0.1, 0.03 and 0.01 respectively. These errors are comparable to or smaller than those arising from current levels of observational uncertainty. Furthermore, we estimate that the GRMHD simulated discs from which these error estimates are obtained correspond to effective disc luminosities of about 0.4-0.7 Eddington, and that the errors will be smaller for discs with luminosities of 0.3 Eddington or less, which are used in the continuum-fitting method. We thus conclude that use of the Novikov-Thorne thin-disc model does not presently limit the accuracy of the continuum-fitting method of measuring black hole spin.

A typical galaxy is thought to contain tens of millions of stellar-mass black holes, the collapsed remnants of once massive stars, and a single nuclear supermassive black hole. Both classes of black holes accrete gas from their environments. The accreting gas forms a flattened orbiting structure known as an accretion disk. During the past several years, it has become possible to obtain measurements of the spins of the two classes of black holes by modeling the X-ray emission from their accretion disks. Two methods are employed, both of which depend upon identifying the inner radius of the accretion disk with the innermost stable circular orbit (ISCO), whose radius depends only on the mass and spin of the black hole. In the Fe K method, which applies to both classes of black holes, one models the profile of the relativistically-broadened iron line with a special focus on the gravitationally redshifted red wing of the line. In the continuum-fitting method, which has so far only been applied to stellar-mass black holes, one models the thermal X-ray continuum spectrum of the accretion disk. We discuss both methods, with a strong emphasis on the continuum-fitting method and its application to stellar-mass black holes. Spin results for eight stellar-mass black holes are summarized. These data are used to argue that the high spins of at least some of these black holes are natal, and that the presence or absence of relativistic jets in accreting black holes is not entirely determined by the spin of the black hole.

The standard general relativistic model of a razor-thin accretion disk around a black hole, developed by Novikov & Thorne (NT) in 1973, assumes the shear stress vanishes at the radius of the innermost stable circular orbit (ISCO) and that, outside the ISCO, the shear stress is produced by an effective turbulent viscosity. However, astrophysical accretion disks are not razor-thin, it is uncertain whether the shear stress necessarily vanishes at the ISCO, and the magnetic field, which is thought to drive turbulence in disks, may contain large-scale structures that do not behave like a simple local scalar viscosity. We describe three-dimensional general relativistic magnetohydrodynamic simulations of accretion disks around black holes with a range of spin parameters, and we use the simulations to assess the validity of the NT model. Our fiducial initial magnetic field consists of multiple (alternating polarity) poloidal field loops whose shape is roughly isotropic in the disk in order to match the isotropic turbulence expected in the poloidal plane. For a thin disk with an aspect ratio |h/r| ~ 0.07 around a non-spinning black hole, we find a decrease in the accreted specific angular momentum of 2.9% relative to the NT model and an excess luminosity from inside the ISCO of 3.5%. The deviations in the case of spinning black holes are of the same order. In addition, the deviations decrease with decreasing |h/r|. We therefore conclude that magnetized thin accretion disks in x-ray binaries in the thermal/high-soft spectral state ought to be well-described by the NT model, especially at luminosities below 30% of Eddington where we expect a very small disk thickness |h/r| <~ 0.05. We also discuss how the stress and the luminosity inside the ISCO depend on the assumed initial magnetic field geometry. (abridged)

We present the combination of dynamical accretion model based on 3D GRMHD simulations and general relativistic (GR) polarized radiative transfer. We write down the formalism of and perform the GR ray-tracing of cyclo-synchrotron radiation through the model of accretion flow in Sagittarius A*. GR polarimetric imaging is presented as well as the results for spectrum for a probable set of spins and orientations. Precise fitting formulae for Faraday rotation and Faraday conversion coefficients are employed for thermal plasma. The axisymmetic flow pattern and the magnetic field geometry correspond to averaged 3D GRMHD simulations near the black hole, whereas the analytic model was used far from the black hole. The density scaling is found by fitting the sub-mm flux. Spin $a=0.7$ and inclination angle $\theta=0.6$ produce the best fit to sub-mm flux and linear polarization fraction.

We use the Novikov-Thorne thin disk model to fit the thermal continuum X-ray spectra of black hole X-ray binaries, and thereby extract the dimensionless spin parameter a* = a/M of the black hole as a parameter of the fit. We summarize the results obtained to date for six systems and describe work in progress on additional systems. We also describe recent methodological advances, our current efforts to make our analysis software fully available to others, and our theoretical efforts to validate the Novikov-Thorne model.

We present new results on the shock around the southwest radio lobe of Centaurus A using data from the Chandra Very Large Programme observations. The X-ray spectrum of the emission around the outer southwestern edge of the lobe is well described by a single power-law model with Galactic absorption -- thermal models are strongly disfavoured, except in the region closest to the nucleus. We conclude that a significant fraction of the X-ray emission around the southwest part of the lobe is synchrotron, not thermal. We infer that in the region where the shock is strongest and the ambient gas density lowest, the inflation of the lobe is accelerating particles to X-ray synchrotron emitting energies, similar to supernova remnants such as SN1006. This interpretation resolves a problem of our earlier, purely thermal, interpretation for this emission, namely that the density compression across the shock was required to be much larger than the theoretically expected factor of 4. We estimate that the lobe is expanding to the southwest with a velocity of ~2600 km/s, roughly Mach 8 relative to the ambient medium. We discuss the spatial variation of spectral index across the shock region, concluding that our observations constrain gamma_max for the accelerated particles to be 10^8 at the strongest part of the shock, consistent with expectations from diffusive shock acceleration theory. Finally, we consider the implications of these results for the production of ultra-high energy cosmic rays (UHECRs) and TeV emission from Centaurus A, concluding that the shock front region is unlikely to be a significant source of UHECRs, but that TeV emission from this region is expected at levels comparable to current limits at TeV energies, for plausible assumed magnetic field strengths.

By combining data from different published 3-D simulations of Keplerian shearing boxes unstable to the magnetorotational instability (MRI), we highlight tight anti-correlations between the total effective inferred angular momentum transport parameter, $\alpha_{tot}$, its separate Maxwell and Reynolds contributions $\alpha_{mag}$ and $\alpha_{kin}$, and the kinetic to magnetic pressure ratio $\beta$, defined with the initial or saturated (when available) thermal pressure. Plots of $Log (\alpha_{kin}), Log (\alpha_{mag})$, and $Log (\alpha_{tot})$ vs $Log (\beta)$ are well fit by straight lines even as $\alpha_{kin}$, $\alpha_{mag}$,and $\alpha_{tot}$ vary by four orders of magnitude over the simulations included. The ratio $\alpha_{kin}/\alpha_{mag}$ and the product $\alpha_{tot}\beta$ are quite constant and largely independent of the presence or absence of weak mean fields, the choice of initial and boundary conditions, and the resolution. In short, simulations have more strongly constrained the product $\alpha_{tot}\beta$ than $\alpha_{tot}$ itself.

We present results from deeper {\em Chandra} observations of the southwest radio lobe of Centaurus A, first described by Kraft et al. (2003). We find that the sharp X-ray surface brightness discontinuity extends around $\sim$75% of the periphery of the radio lobe, and detect significant temperature jumps in the brightest regions of this discontinuity nearest to the nucleus. This demonstrates that this discontinuity is indeed a strong shock which is the result of an overpressure which has built up in the entire lobe over time. Additionally, we demonstrate that if the mean free path for ions to transfer energy and momentum to the electrons behind the shock is as large as the Spitzer value, the electron and proton temperatures will not have equilibrated along the SW boundary of the radio lobe where the shock is strongest. Thus the proton temperature of the shocked gas could be considerably larger than the observed electron temperature, and the total energy of the outburst correspondingly larger as well. We investigate this using a simple one-dimensional shock model for a two-fluid (proton/electron) plasma. We find that for the thermodynamic parameters of the Cen A shock the electron temperature rises rapidly from $\sim$0.29 keV (the temperature of the ambient ISM) to $\sim$3.5 keV at which point heating from the protons is balanced by adiabatic losses. The proton and electron temperatures do not equilibrate in a timescale less than the age of the lobe. We note that the measured electron temperature of similar features in other nearby powerful radio galaxies in poor environments may considerably underestimate the strength and velocity of the shock.

We investigate the relation between Forbush cosmic ray decrease recovery time and coronal mass ejection transit time between the Sun and Earth. We identify 17 Forbush decreases from ground based neutron count rates between 1978 and 2003 that occur during the same phase in the solar cycle and can be associated with single coronal mass ejections (CMEs) in the SOHO LASCO CME Catalog or previously published reports, and with specific interplanetary coronal mass ejections (ICMEs) crossing the vicinity of Earth. We find an anti-correlation between Forbush recovery times and CME transit time that contradicts the predictions of simple cosmic ray diffusive barrier models. The anti-correlation suggests that the decay rate of ICMEs is anti-correlated with their travel speed. Forbush recovery times range from seven times the transit time for the fastest disturbance to a fifth the Sun-Earth transit time for the slowest. To account for the large range of measured recovery times we infer that the slowest disturbances must decay rapidly with radius whereas the fastest ones must remain strong. The longest recovery times suggest that the fastest disturbances in our sample decayed less rapidly with radius than the ambient solar wind magnetic field strength.