Researcher profile

Robert F. Penna

Robert F. Penna contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

21 published item(s)

preprint2022arXiv

Gravitational Pulse Wave Scattering and the Geroch group

Cylindrical gravitational pulse waves are cylindrical gravitational waves with a pulse profile in the radial direction. The dynamics of cylindrical gravitational pulse waves is governed by a two dimensional integrable sigma model with an infinite dimensional hidden symmetry called the Geroch group. Piran, Safier, and Katz obtained an exact solution for a single pulse wave by making a double analytic continuation of the Kerr metric. Their solution describes a pulse wave that arrives from infinity, bounces off the symmetry axis, and returns to infinity. In this note, we obtain an exact solution for a double pulse wave by making a double analytic continuation of the double Kerr metric. The new solution describes a pair of pulse waves that arrive from infinity, scatter through each other, and return to infinity. The pulses pass through each other and emerge with their shapes intact, as in ordinary soliton scattering. Unlike ordinary soliton scattering, there is no time delay. The metric is free of the conical singularity that appears in the double Kerr metric. In the future, it would be interesting to understand how the Geroch group manifests in the quantum version of gravitational pulse wave scattering.

preprint2016arXiv

BMS invariance and the membrane paradigm

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.

preprint2015arXiv

Black hole jet power from impedance matching

Black hole jet power depends on the angular velocity of magnetic field lines, $Ω_F$. Force-free black hole magnetospheres typically have $Ω_F/Ω_H \approx 0.5$, where $Ω_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 $Ω_F/Ω_H$ because the force-free approximation breaks down in the accretion flow.

preprint2015arXiv

Energy extraction from boosted black holes: Penrose process, jets, and the membrane at infinity

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 Φ^2/(4π)$, where $v$ is the hole's velocity, $M$ is its mass, and $Φ$ 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.

preprint2014arXiv

Black hole Meissner effect and entanglement

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.

preprint2014arXiv

MITEoR: A Scalable Interferometer for Precision 21 cm Cosmology

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

preprint2014arXiv

The Black Hole Meissner Effect and Blandford-Znajek Jets

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 &#34;Meissner effect&#34; 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.

preprint2013arXiv

A new equilibrium torus solution and GRMHD initial conditions

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.

preprint2013arXiv

Are entangled particles connected by wormholes? Support for the ER=EPR conjecture from entropy inequalities

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.

preprint2013arXiv

General Relativistic Magnetohydrodynamic Simulations of Blandford-Znajek Jets and the Membrane Paradigm

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.

preprint2012arXiv

Constraining the Accretion Flow in Sgr A* by General Relativistic Dynamical and Polarized Radiative Modeling

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.

preprint2012arXiv

GRMHD Simulations of Magnetized Advection Dominated Accretion on a Non-Spinning Black Hole: Role of Outflows

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.

preprint2012arXiv

Sagittarius A* Accretion Flow and Black Hole Parameters from General Relativistic Dynamical and Polarized Radiative Modeling

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 θ, 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, θ=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 θ=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.

preprint2012arXiv

The Eye of the Storm: Light from the Inner Plunging Region of Black Hole Accretion Discs

It is generally thought that the light coming from the inner plunging region of black hole accretion discs contributes negligibly to the disc&#39;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&#39;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&#39; emission observed in the thermal dominant and steep power-law states of X-ray binaries.

preprint2012arXiv

The Shakura-Sunyaev Viscosity Prescription with Variable alpha(r)

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&#39;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&#39;s observed in simulations. It is possible that Newtonian simulations find smaller alpha&#39;s than general relativistic simulations because the shear parameter is lower in Newtonian flows.

preprint2011arXiv

Apparent motion of a spherical shell collapsing onto a black hole

We model the collapse of a spherically symmetric, constant density, pressureless shell onto a preexisting black hole. Contrary to a recent claim of Liu and Zhang (2009), we show that an observer at infinity never sees any part of the shell cross the event horizon. The entire shell appears to &#34;freeze&#34; outside the black hole. We find that during intermediate stages of the collapse, the change in apparent size of the inner surface is non-monotonic. However, at late times, the apparent sizes of the inner and outer surfaces of the shell both approach 3\sqrt{3}GM/c^2 asymptotically, where M is the combined mass of the initial black hole and the shell.

preprint2011arXiv

Measuring Black Hole Spin by the Continuum-Fitting Method: Effect of Deviations from the Novikov-Thorne Disc Model

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.

preprint2011arXiv

Measuring the Spins of Accreting Black Holes

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.

preprint2011arXiv

Thin Disk Theory with a Non-Zero Torque Boundary Condition and Comparisons with Simulations

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.

preprint2010arXiv

Simulations of Magnetized Disks Around Black Holes: Effects of Black Hole Spin, Disk Thickness, and Magnetic Field Geometry

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)

preprint2005arXiv

The decay of interplanetary coronal mass ejections and forbush decrease recovery times

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.