Abstract
We describe ipole, a new public ray-tracing code for covariant, polarized radiative transport. The code extends the ibothros scheme for covariant, unpolarized transport using two ...representations of the polarized radiation field: In the coordinate frame, it parallel transports the coherency tensor; in the frame of the plasma it evolves the Stokes parameters under emission, absorption, and Faraday conversion. The transport step is implemented to be as spacetime- and coordinate- independent as possible. The emission, absorption, and Faraday conversion step is implemented using an analytic solution to the polarized transport equation with constant coefficients. As a result, ipole is stable, efficient, and produces a physically reasonable solution even for a step with high optical depth and Faraday depth. We show that the code matches analytic results in flat space, and that it produces results that converge to those produced by Dexter's grtrans
polarized transport code on a complicated model problem. We expect ipole will mainly find applications in modelling Event Horizon Telescope sources, but it may also be useful in other relativistic transport problems such as modelling for the IXPE mission.
Simple assumptions made regarding electron thermodynamics often limit the extent to which general relativistic magnetohydrodynamic (GRMHD) simulations can be applied to observations of low-luminosity ...accreting black holes. We present, implement, and test a model that self-consistently evolves an entropy equation for the electrons and takes into account the effects of spatially varying electron heating and relativistic anisotropic thermal conduction along magnetic field lines. We neglect the backreaction of electron pressure on the dynamics of the accretion flow. Our model is appropriate for systems accreting at ...10... of the Eddington accretion rate, so radiative cooling by electrons can be neglected. It can be extended to higher accretion rates in the future by including electron cooling and proton-electron Coulomb collisions. We present a suite of tests showing that our method recovers the correct solution for electron heating under a range of circumstances, including strong shocks and driven turbulence. Our initial applications to axisymmetric simulations of accreting black holes show that (1) physically motivated electron heating rates that depend on the local magnetic field strength yield electron temperature distributions significantly different from the constant electron-to-proton temperature ratios assumed in previous work, with higher electron temperatures concentrated in the coronal region between the disc and the jet; (2) electron thermal conduction significantly modifies the electron temperature in the inner regions of black hole accretion flows if the effective electron mean free path is larger than the local scaleheight of the disc (at least for the initial conditions and magnetic field configurations we study). The methods developed in this work are important for producing more realistic predictions for the emission from accreting black holes such as Sagittarius A* and M87; these applications will be explored in future work. (ProQuest: ... denotes formulae/symbols omitted.)
Abstract
We calculate the radiative properties of Sagittarius A* – spectral energy distribution, variability and radio-infrared images – using the first 3D, physically motivated black hole accretion ...models that directly evolve the electron thermodynamics in general relativistic MHD simulations. These models reproduce the coupled disc-jet structure for the emission favoured by previous phenomenological analytic and numerical works. More specifically, we find that the low frequency radio emission is dominated by emission from a polar outflow while the emission above 100 GHz is dominated by the inner region of the accretion disc. The latter produces time variable near-infrared (NIR) and X-ray emission, with frequent flaring events (including IR flares without corresponding X-ray flares and IR flares with weak X-ray flares). The photon ring is clearly visible at 230 GHz and 2 μm, which is encouraging for future horizon-scale observations. We also show that anisotropic electron thermal conduction along magnetic field lines has a negligible effect on the radiative properties of our model. We conclude by noting limitations of our current generation of first-principles models, particularly that the outflow is closer to adiabatic than isothermal and thus underpredicts the low frequency radio emission.
ABSTRACT We present bhlight, a numerical scheme for solving the equations of general relativistic radiation magnetohydrodynamics using a direct Monte Carlo solution of the frequency-dependent ...radiative transport equation. bhlight is designed to evolve black hole accretion flows at intermediate accretion rate, in the regime between the classical radiatively efficient disk and the radiatively inefficient accretion flow (RIAF), in which global radiative effects play a sub-dominant but non-negligible role in disk dynamics. We describe the governing equations, numerical method, idiosyncrasies of our implementation, and a suite of test and convergence results. We also describe example applications to radiative Bondi accretion and to a slowly accreting Kerr black hole in axisymmetry.
Measurements of stellar orbits provide compelling evidence that the compact radio source Sagittarius A* at the Galactic Centre is a black hole four million times the mass of the Sun. With the ...exception of modest X-ray and infrared flares, Sgr A* is surprisingly faint, suggesting that the accretion rate and radiation efficiency near the event horizon are currently very low. Here we report the presence of a dense gas cloud approximately three times the mass of Earth that is falling into the accretion zone of Sgr A*. Our observations tightly constrain the cloud's orbit to be highly eccentric, with an innermost radius of approach of only ∼3,100 times the event horizon that will be reached in 2013. Over the past three years the cloud has begun to disrupt, probably mainly through tidal shearing arising from the black hole's gravitational force. The cloud's dynamic evolution and radiation in the next few years will probe the properties of the accretion flow and the feeding processes of the supermassive black hole. The kilo-electronvolt X-ray emission of Sgr A* may brighten significantly when the cloud reaches pericentre. There may also be a giant radiation flare several years from now if the cloud breaks up and its fragments feed gas into the central accretion zone.
Full text
Available for:
DOBA, IJS, IZUM, KILJ, KISLJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Conservative numerical schemes for general relativistic magnetohydrodynamics (GRMHD) require a method for transforming between "conserved" variables such as momentum and energy density and ..."primitive" variables such as rest-mass density, internal energy, and components of the four-velocity. The forward transformation (primitive to conserved) has a closed-form solution, but the inverse transformation (conserved to primitive) requires the solution of a set of five nonlinear equations. Here we discuss the mathematical properties of the inverse transformation and present six numerical methods for performing the inversion. The first method solves the full set of five nonlinear equations directly using a Newton-Raphson scheme and a guess from the previous time step. The other methods reduce the five nonlinear equations to either one or two nonlinear equations that are solved numerically. Comparisons between the methods are made using a survey over phase space, a two-dimensional explosion problem, and a general relativistic MHD accretion disk simulation. The run time of the methods is also examined. Code implementing the schemes is available with the electronic edition of the article.
The Moon is believed to have formed in the aftermath of a giant impact between a planetary-mass body and the proto-Earth. In a typical giant impact scenario, a disk of vapor, liquid, and solid debris ...forms around the proto-Earth and-after possibly decades of evolution-condenses to form the Moon. Using state-of-the-art numerical simulations, we investigate the dynamical effects of magnetic fields on the Moon-forming giant impact. We show that turbulence generated by the collision itself, shear in the boundary layer between the post-impact debris field and the proto-Earth, and turbulence in the vapor component of the disk amplify the field to dynamically significant strengths. Magnetically driven turbulence promotes angular momentum transport in the protolunar disk. Debris material is accreted onto the proto-Earth, making Moon formation less efficient, while the disk is forced to spread to larger radii, cooling at its outer edge. Magnetic fields speed the evolution of the vapor component of the protolunar disk and hasten the formation of the Moon.
Numerical simulations of self-gravitating flows evolve a momentum equation and an energy equation that account for accelerations and gravitational energy releases due to a time-dependent ...gravitational potential. In this work, we implement a fully conservative numerical algorithm for self-gravitating flows, using source terms, in the astrophysical magnetohydrodynamics framework Athena++. We demonstrate that properly evaluated source terms are conservative when they are equivalent to the divergence of a corresponding "gravity flux" (i.e., a gravitational stress tensor or a gravitational energy flux). We provide test problems that demonstrate several advantages of the source-term-based algorithm, including second-order convergence and round-off error total momentum and total energy conservation. The fully conservative scheme suppresses anomalous accelerations that arise when applying a common numerical discretization of the gravitational stress tensor that does not guarantee curl-free gravity.
ABSTRACT The leading theory for the formation of Earth's Moon invokes a collision between a Mars-sized body and the proto-Earth to produce a disk of orbiting material that later condenses to form the ...Moon. We show that the disk opacity is large, and cooling is therefore inefficient ( ). In this regime, angular momentum transport in the disk leads to steady heating unless . Following earlier work by Charnoz and Michaut, and Carballido et al., we show that once the disk is completely vaporized it is well coupled to the magnetic field. We consider a scenario in which turbulence driven by magnetic fields leads to a brief, hot phase where the disk is geometrically thick, with strong turbulent mixing. The disk cools by spreading until it decouples from the field. We point out that approximately half the accretion energy is dissipated in the boundary layer where the disk meets the Earth's surface. This creates high entropy material close to the Earth, driving convection and mixing. Finally, a hot magnetized disk could drive bipolar outflows that remove mass and angular momentum from the Earth-Moon system.
Some active galactic nuclei, microquasars, and gamma-ray bursts may be powered by the electromagnetic braking of a rapidly rotating black hole. We investigate this possibility via axisymmetric ...numerical simulations of a black hole surrounded by a magnetized plasma. The plasma is described by the equations of general relativistic magnetohydrodynamics, and the effects of radiation are neglected. The evolution is followed for 2000GM/c super(3), and the computational domain extends from inside the event horizon to typically 40GM/c super(2). We compare our results to two analytic steady state models, including the force-free magnetosphere of Blandford & Znajek. Along the way we present a self-contained rederivation of the Blandford-Znajek model in Kerr- Schild (horizon penetrating) coordinates. We find that (1) low-density polar regions of the numerical models agree well with the Blandford-Znajek model, (2) many of our models have an outward Poynting flux on the horizon in the Kerr- Schild frame, (3) none of our models have a net outward energy flux on the horizon, and (4) one of our models, in which the initial disk has net magnetic flux, shows a net outward angular momentum flux on the horizon. We conclude with a discussion of the limitations of our model, astrophysical implications, and problems to be addressed by future numerical experiments.