Source author record

Monika Moscibrodzka

Monika Moscibrodzka appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

23 published item(s)

preprint2022arXiv

MeqSilhouette v2: Spectrally-resolved polarimetric synthetic data generation for the Event Horizon Telescope

We present MeqSilhouette v2.0 (MeqSv2), a fully polarimetric, time-and frequency-resolved synthetic data generation software for simulating millimetre (mm) wavelength very long baseline interferometry (VLBI) observations with heterogeneous arrays. Synthetic data are a critical component in understanding real observations, testing calibration and imaging algorithms, and predicting performance metrics of existing or proposed sites. MeqSv2 applies physics-based instrumental and atmospheric signal corruptions constrained by empirically-derived site and station parameters to the data. The new version is capable of applying instrumental polarization effects and various other spectrally-resolved effects using the Radio Interferometry Measurement Equation (RIME) formalism and produces synthetic data compatible with calibration pipelines designed to process real data. We demonstrate the various corruption capabilities of MeqSv2 using different arrays, with a focus on the effect of complex bandpass gains on closure quantities for the EHT at 230 GHz. We validate the frequency-dependent polarization leakage implementation by performing polarization self-calibration of synthetic EHT data using PolSolve. We also note the potential applications for cm-wavelength VLBI array analysis and design and future directions.

preprint2022arXiv

Millimeter light curves of Sagittarius A* observed during the 2017 Event Horizon Telescope campaign

The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very-long-baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare, and exhibit strongly enhanced variability. All of the light curves are consistent with a red noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 min and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 h. The spectral energy distribution is flat at 220 GHz and there are no time-lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines.

preprint2022arXiv

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm, ipole, and igrmonty, paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.

preprint2022arXiv

Resolving the inner parsec of the blazar J1924-2914 with the Event Horizon Telescope

The blazar J1924-2914 is a primary Event Horizon Telescope (EHT) calibrator for the Galactic Center's black hole Sagittarius A*. Here we present the first total and linearly polarized intensity images of this source obtained with the unprecedented 20 $μ$as resolution of the EHT. J1924-2914 is a very compact flat-spectrum radio source with strong optical variability and polarization. In April 2017 the source was observed quasi-simultaneously with the EHT (April 5-11), the Global Millimeter VLBI Array (April 3), and the Very Long Baseline Array (April 28), giving a novel view of the source at four observing frequencies, 230, 86, 8.7, and 2.3 GHz. These observations probe jet properties from the subparsec to 100-parsec scales. We combine the multi-frequency images of J1924-2914 to study the source morphology. We find that the jet exhibits a characteristic bending, with a gradual clockwise rotation of the jet projected position angle of about 90 degrees between 2.3 and 230 GHz. Linearly polarized intensity images of J1924-2914 with the extremely fine resolution of the EHT provide evidence for ordered toroidal magnetic fields in the blazar compact core.

preprint2022arXiv

The Photon Ring in M87*

We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.

preprint2020arXiv

On the prospects of imaging Sagittarius A* from space

Very Long Baseline Interferometry (VLBI) at sub-millimeter waves has the potential to image the shadow of the black hole in the Galactic Center, Sagittarius A* (Sgr A*), and thereby test basic predictions of the theory of general relativity. We investigate the imaging prospects of a new Space VLBI mission concept. The setup consists of two satellites in polar or equatorial circular Medium-Earth Orbits with slightly different radii, resulting in a dense spiral-shaped uv-coverage with long baselines, allowing for extremely high-resolution and high-fidelity imaging of radio sources. We simulate observations of a general relativistic magnetohydrodynamics model of Sgr A* for this configuration with noise calculated from model system parameters. After gridding the $uv$-plane and averaging visibilities accumulated over multiple months of integration, images of Sgr A* with a resolution of up to 4 $μ$as could be reconstructed, allowing for stronger tests of general relativity and accretion models than with ground-based VLBI.

preprint2019arXiv

General relativistic polarized radiative transfer with inverse Compton scatterings

We present {\tt radpol} - a numerical scheme for integrating multifrequency polarized radiative transfer equations along rays propagating in a curved spacetime. The scheme includes radiative processes such as synchrotron emission, absorption, Faraday rotation and conversion, and, for the first time, relativistic Compton scatterings including effects of light polarization. The scheme is fully covariant and is applicable to model radio-$γ$-ray emission and its polarization from, e.g., relativistic jets and accretion flows onto black holes and other exotic objects described in alternative metric theories and modeled semi-analytically or with time-dependent magnetohydrodynamical simulations. We perform a few tests to validate the implemented numerical algorithms that handle light polarization in curved spacetime. We demonstrate application of the scheme to model broadband emission spectra from a relativistically hot, geometrically thick coronal-like inflow around a supermassive black hole where the disk model is realized in a general relativistic magnetohydrodynamical simulation.

preprint2019arXiv

Scale-invariance of black hole accretion: modeling emission from a black hole X-ray binary with relativistic accretion flow simulations

We model non-thermal emission spectrum of the extremely sub-Eddington X-ray binary system A0620-00. It is believed that this non-thermal emission is produced by a radiatively inefficient "quiescent" accretion onto a stellar-mass black hole present in the system. We post-process GRMHD simulations with multiwavelength, fully polarized, relativistic radiative transfer calculations to predict broadband spectra and emission polarization levels for a range of electron models and accretion rates. We find that a model with strong coupling of electrons and ions in the accretion disk and accretion rate of only $\dot{M}=3\times10^{-13} [M_\odot yr^{-1}]$ is able to recover the observed X-ray spectral slope as well as the excess of linear polarization detected in the source in NIR/optical bands. Our models constrain the spectral properties of a putative relativistic jet produced in this system. In addition we show that the magnetized winds from our hot accretion flow carry away a small fraction of the orbital angular momentum of the binary unable to explain the observed rapid orbital decay of the system. Similar to the present GRMHD simulations are often used to explain emission from sub-Eddington supermassvie black holes in Sgr A* or M87; the present simulations allow us to test whether some aspects of the quiescent black hole accretion are scale invariant.

preprint2016arXiv

General relativistic magnetohydrodynamical simulations of the jet in M87

(abridged) The connection between black hole, accretion disk, and radio jet can be best constrained by fitting models to observations of nearby low luminosity galactic nuclei, in particular the well studied sources Sgr~A* and M87. There has been considerable progress in modeling the central engine of active galactic nuclei by an accreting supermassive black hole coupled to a relativistic plasma jet. However, can a single model be applied to a range of black hole masses and accretion rates? Here we want to compare the latest three-dimensional numerical model, originally developed for Sgr A* in the center of the Milky Way, to radio observations of the much more powerful and more massive black hole in M87. We postprocess three-dimensional GRMHD models of a jet-producing radiatively inefficient accretion flow around a spinning black hole using relativistic radiative transfer and ray-tracing to produce model spectra and images. As a key new ingredient to these models, we allow the proton-electron coupling in these simulations depend on the magnetic properties of the plasma. We find that the radio emission in M87 is well described by a combination of a two-temperature accretion flow and a hot single-temperature jet. The model fits the basic observed characteristics of the M87 radio core. The best fit model has a mass-accretion rate of Mdot approx 9x10^{-3} MSUN/YR and a total jet power of P_j \sim 10^{43} erg/s. Emission at 1.3mm is produced by the counter jet close to the event horizon. Its characteristic crescent shape surrounding the black hole shadow could be resolved by future millimeter-wave VLBI experiments. The model was successfully derived from one for the supermassive black hole in center of the Milky Way by appropriately scaling mass and accretion rate. This suggests the possibility that this model could also apply to a larger range of low-luminosity black holes.

preprint2016arXiv

The Proper Motion of the Galactic Center Pulsar Relative to Sagittarius A*

We measure the proper motion of the pulsar PSR J1745-2900 relative to the Galactic Center massive black hole, Sgr A*, using the Very Long Baseline Array (VLBA). The pulsar has a transverse velocity of 236 +/- 11 km s^-1 at position angle 22 +/- 2 deg East of North at a projected separation of 0.097 pc from Sgr A*. Given the unknown radial velocity, this transverse velocity measurement does not conclusively prove that the pulsar is bound to Sgr A*; however, the probability of chance alignment is very small. We do show that the velocity and position is consistent with a bound orbit originating in the clockwise disk of massive stars orbiting Sgr A* and a natal velocity kick of <~ 500 km s^-1. An origin among the isotropic stellar cluster is possible but less probable. If the pulsar remains radio-bright, multi-year astrometry of PSR J1745-2900 can detect its acceleration and determine the full three-dimensional orbit. We also demonstrate that PSR J1745-2900 exhibits the same angular broadening as Sgr A* over a wavelength range of 3.6 cm to 0.7 cm, further confirming that the two sources share the same interstellar scattering properties. Finally, we place the first limits on the presence of a wavelength-dependent shift in the position of Sgr A*, i.e., the core shift, one of the expected properties of optically-thick jet emission. Our results for PSR J1745-2900 support the hypothesis that Galactic Center pulsars will originate from the stellar disk and deepens the mystery regarding the small number of detected Galactic Center pulsars.

preprint2015arXiv

ALMA and VLA measurements of frequency-dependent time lags in Sagittarius A*: evidence for a relativistic outflow

Radio and mm-wavelength observations of Sagittarius A* (Sgr A*), the radio source associated with the supermassive black hole at the center of our Galaxy, show that it behaves as a partially self-absorbed synchrotron-emitting source. The measured size of Sgr A* shows that the mm-wavelength emission comes from a small region and consists of the inner accretion flow and a possible collimated outflow. Existing observations of Sgr A* have revealed a time lag between light curves at 43 GHz and 22 GHz, which is consistent with a rapidly expanding plasma flow and supports the presence of a collimated outflow from the environment of an accreting black hole. Here we wish to measure simultaneous frequency-dependent time lags in the light curves of Sgr A* across a broad frequency range to constrain direction and speed of the radio-emitting plasma in the vicinity of the black hole. Light curves of Sgr A* were taken in May 2012 using ALMA at 100 GHz using the VLA at 48, 39, 37, 27, 25.5, and 19 GHz. As a result of elevation limits and the longitude difference between the stations, the usable overlap in the light curves is approximately four hours. Although Sgr A* was in a relatively quiet phase, the high sensitivity of ALMA and the VLA allowed us to detect and fit maxima of an observed minor flare where flux density varied by ~10%. The fitted times of flux density maxima at frequencies from 100 GHz to 19 GHz, as well as a cross-correlation analysis, reveal a simple frequency-dependent time lag relation where maxima at higher frequencies lead those at lower frequencies. Taking the observed size-frequency relation of Sgr A* into account, these time lags suggest a moderately relativistic (lower estimates: 0.5c for two-sided, 0.77c for one-sided) collimated outflow.

preprint2014arXiv

Observational appearance of inefficient accretion flows and jets in 3D GRMHD simulations: Application to Sgr~A*

Radiatively inefficient accretion flows (RIAFs) are believed to power supermassive black holes (SMBH) in the underluminous cores of galaxies. Such black holes are typically accompanied by flat-spectrum radio cores indicating the presence of moderately relativistic jets. One of the best constrained RIAFs is associated with the SMBH in the Galactic center, Sgr A*. Since the plasma in RIAFs is only weakly collisional, the dynamics and the radiative properties of these systems are very uncertain. Here we want to study the impact of varying electron temperature on the appearance of accretion flows and jets. Using 3-D GRMHD accretion flow simulations, we use ray tracing methods to predict spectra and radio images of RIAFs allowing for different electron heating mechanisms in the in- and outflowing parts of the simulations. We find that small changes in the electron temperature can result in dramatic differences in the relative dominance of jets and accretion flows. Application to Sgr A* shows that radio spectrum and size of this source can be well reproduced with a model where electrons are more efficiently heated in the jet. The X-ray emission is sensitive to the electron heating mechanism in the jets and disk and therefore X-ray observations put strong constraints on electron temperatures and geometry of the accretion flow and jet. For Sgr A*, the jet model also predicts a significant frequency-dependent core shift which could place independent constraints on the model once measured accurately. We conclude that more sophisticated models for electron distribution functions are crucial for constraining GRMHD simulations with actual observations. For Sgr A*, the radio appearance may well be dominated by the outflowing plasma. Nonetheless, at the highest radio frequencies, the shadow of the event horizon should still be detectable with future Very Long Baseline Interferometric observations.

preprint2013arXiv

Accretion and outflow from a magnetized, neutrino cooled torus around the gamma ray burst central engine

We calculate the structure and short-term evolution of a gamma ray burst central engine in the form of a turbulent torus accreting onto a stellar mass black hole. Our models apply to the short gamma ray burst events, in which a remnant torus forms after the neutron star-black hole or a double neutron star merger and is subsequently accreted. We study the 2-dimensional, relativistic models and concentrate on the effects of black hole and flow parameters as well as the neutrino cooling. We compare the resulting structure and neutrino emission to the results of our previous 1-dimensional simulations. We find that the neutrino cooled torus launches a powerful mass outflow, which contributes to the total neutrino luminosity and mass loss from the system. The neutrino luminosity may exceed the Blandford-Znajek luminosity of the polar jets and the subsequent annihilation of neutrino-antineutrino pairs will provide an additional source of power to the GRB emission.

preprint2012arXiv

The Galactic Center Weather Forecast

In accretion-based models for Sgr A* the X-ray, infrared, and millimeter emission arise in a hot, geometrically thick accretion flow close to the black hole. The spectrum and size of the source depend on the black hole mass accretion rate $\dot{M}$. Since Gillessen et al. have recently discovered a cloud moving toward Sgr A* that will arrive in summer 2013, $\dot{M}$ may increase from its present value $\dot{M}_0$. We therefore reconsider the "best-bet" accretion model of Moscibrodzka et al., which is based on a general relativistic MHD flow model and fully relativistic radiative transfer, for a range of $\dot{M}$. We find that for modest increases in $\dot{M}$ the characteristic ring of emission due to the photon orbit becomes brighter, more extended, and easier to detect by the planned Event Horizon Telescope submm VLBI experiment. If $\dot{M} \gtrsim 8 \dot{M}_0$ this "silhouette of the black hole will be hidden beneath the synchrotron photosphere at 230 GHz, and for $\dot{M} \gtrsim 16 \dot{M}_0$ the silhouette is hidden at 345 GHz. We also find that for $\dot{M} > 2 \dot{M}_0$ the near-horizon accretion flow becomes a persistent X-ray and mid-infrared source, and in the near-infrared Sgr A* will acquire a persistent component that is brighter than currently observed flares.

preprint2011arXiv

Accretion and outflow from a magnetized, neutrino cooled torus in the gamma ray burst central engine

Gamma Ray Bursts (GRB) are the extremely energetic transient events, visible from the most distant parts of the Universe. They are most likely powered by accretion on the hyper-Eddington rates that proceeds onto a newly born stellar mass black hole. This central engine gives rise to the most powerful, high Lorentz factor jets that are responsible for energetic gamma ray emission. We investigate the accretion flow evolution in GRB central engine, using the 2D MHD simulations in General Relativity. We compute the structure and evolution of the extremely hot and dense torus accreting onto the fast spinning black hole, which launches the magnetized jets. We calculate the chemical structure of the disk and account for neutrino cooling. Our preliminary runs apply to the short GRB case (remnant torus accreted after NS-NS or NS-BH merger). We estimate the neutrino luminosity of such an event for chosen disk and central BH mass

preprint2011arXiv

Pair Production in Low Luminosity Galactic Nuclei

Electron-positron pairs may be produced near accreting black holes by a variety of physical processes, and the resulting pair plasma may be accelerated and collimated into a relativistic jet. Here we use a self-consistent dynamical and radiative model to investigate pair production by γγcollisions in weakly radiative accretion flows around a black hole of mass M and accretion rate \dot{M}. Our flow model is drawn from general relativistic magnetohydrodynamic simulations, and our radiation field is computed by a Monte Carlo transport scheme assuming the electron distribution function is thermal. We argue that the pair production rate scales as r^{-6} M^{-1} \dot{M}^{6}. We confirm this numerically and calibrate the scaling relation. This relation is self-consistent in a wedge in M, \dot{M} parameter space. If \dot{M} is too low the implied pair density over the poles of the black hole is below the Goldreich-Julian density and γγpair production is relatively unimportant; if \dot{M} is too high the models are radiatively efficient. We also argue that for a power-law spectrum the pair production rate should scale with the observables L_X \equiv X-ray luminosity and M as L_X^2 M^{-4}. We confirm this numerically and argue that this relation likely holds even for radiatively efficient flows. The pair production rates are sensitive to black hole spin and to the ion-electron temperature ratio which are fixed in this exploratory calculation. We finish with a brief discussion of the implications for Sgr A* and M87.

preprint2010arXiv

Probing the accretion disk - jet connection via instabilities in the inner accretion flow. From microquasars to quasars

We present various instability mechanisms in the accreting black hole systems which might indicate at the connection between the accretion disk and jet. The jets observed in microquasars can have a persistent or blobby morphology. Correlated with the accretion luminosity, this might provide a link to the cyclic outbursts of the disk. Such duty-cycle type of behavior on short timescales results from the thermal instability caused by the radiation pressure domination. The same type of instability may explain the cyclic radioactivity of the supermassive black hole systems. The somewhat longer timescales are characteristic for the instability caused by the partial hydrogen ionization. The distortions of the jet direction and complex morphology of the sources can be caused by precession of the disk-jet axis.

preprint2010arXiv

The Jet in the Galactic Center: An Ideal Laboratory for Magnetohydrodynamics and General Relativity

In this paper we review and discuss some of the intriguing properties of the Galactic Center supermassive black hole candidate Sgr A*. Of all possible black hole sources, the event horizon of Sgr A*, subtends the largest angular scale on the sky. It is therefore a prime candidate to study and image plasma processes in strong gravity and it even allows imaging of the shadow cast by the event horizon. Recent mm-wave VLBI and radio timing observations as well as numerical GRMHD simulations now have provided several breakthroughs that put Sgr A* back into the focus. Firstly, VLBI observations have now measured the intrinsic size of Sgr A* at multiple frequencies, where the highest frequency measurements have approached the scale of the black hole shadow. Moreover, measurements of the radio variability show a clear time lag between 22 GHz and 43 GHz. The combination of size and timing measurements, allows one to actually measure the flow speed and direction of magnetized plasma at some tens of Schwarzschild radii. This data strongly support a moderately relativistic outflow, consistent with an accelerating jet model. This is compared to recent GRMHD simulation that show the presence of a moderately relativistic outflow coupled to an accretion flow Sgr A*. Further VLBI and timing observations coupled to simulations have the potential to map out the velocity profile from 5-40 Schwarzschild radii and to provide a first glimpse at the appearance of a jet-disk system near the event horizon. Future submm-VLBI experiments would even be able to directly image those processes in strong gravity and directly confirm the presence of an event horizon.

preprint2009arXiv

grmonty: a Monte Carlo Code for Relativistic Radiative Transport

We describe a Monte Carlo radiative transport code intended for calculating spectra of hot, optically thin plasmas in full general relativity. The version we describe here is designed to model hot accretion flows in the Kerr metric and therefore incorporates synchrotron emission and absorption, and Compton scattering. The code can be readily generalized, however, to account for other radiative processes and an arbitrary spacetime. We describe a suite of test problems, and demonstrate the expected $N^{-1/2}$ convergence rate, where $N$ is the number of Monte Carlo samples. Finally we illustrate the capabilities of the code with a model calculation, a spectrum of the slowly accreting black hole Sgr A* based on data provided by a numerical general relativistic MHD model of the accreting plasma.

preprint2009arXiv

Magnetized Accretion Flows: Effects of Gas Pressure

We study how axisymmetric magnetohydrodynamical (MHD) accretion flows depend on gamma adiabatic index in the polytropic equation of state. This work is an extension of Moscibrodzka & Proga (2008), where we investigated the gamma dependence of 2-D Bondi-like accretion flows in the hydrodynamical (HD) limit. Our main goal is to study if simulations for various gamma can give us insights into to the problem of various modes of accretion observed in several types of accretion systems such as black hole binaries (BHB), active galactic nuclei (AGN), and gamma-ray bursts (GRBs). We find that for gamma >~ 4/3, the fast rotating flow forms a thick torus that is supported by rotation and gas pressure. As shown before for gamma=5/3, such a torus produces a strong, persistent bipolar outflow that can significantly reduce the polar funnel accretion of a slowly rotating flow. For low gamma, close to 1, the torus is thin and is supported by rotation. The thin torus produces an unsteady outflow which is too weak to propagate throughout the polar funnel inflow. Compared to their HD counterparts, the MHD simulations show that the magnetized torus can produce an outflow and does not exhibit regular oscillations. Generally, our simulations demonstrate how the torus thickness affects the outflow production. They also support the notion that the geometrical thickness of the torus correlates with the power of the torus outflow. Our results, applied to observations, suggest that the torus ability to radiatively cool and become thin can correspond to a suppression of a jet as observed in the BHB during a transition from a hard/low to soft/high spectral state and a transition from a quiescent to hard/low state in AGN.

preprint2009arXiv

Radiative Models of Sgr A* from GRMHD Simulations

Using flow models based on axisymmetric general relativistic magnetohydrodynamics (GRMHD) simulations, we construct radiative models for sgr A*. Spectral energy distributions that include the effects of thermal synchrotron emission and absorption, and Compton scattering, are calculated using a Monte Carlo technique. Images are calculated using a ray-tracing scheme. All models are scaled so that the 230 GHz flux density is 3.4 Jy. The key model parameters are the dimensionless black hole spin a*, the inclination i, and the ion-to-electron temperature ratio Ti/Te. We find that: (1) models with Ti/Te=1 are inconsistent with the observed submillimeter spectral slope; (2) the X-ray flux is a strongly increasing function of a*; (3) the X-ray flux is a strongly increasing function of i; (4) 230 GHz image size is a complicated function of i, a*, and Ti/Te, but the Ti/Te = 10 models are generally large and at most marginally consistent with the 230 GHz VLBI data; (5) for models with Ti/Te=10 and i=85 deg the event horizon is cloaked behind a synchrotron photosphere at 230 GHz and will not be seen by VLBI, but these models overproduce NIR and X-ray flux; (6) in all models whose SEDs are consistent with observations the event horizon is uncloaked at 230 GHz; (7) the models that are most consistent with the observations have a* \sim 0.9. We finish with a discussion of the limitations of our model and prospects for future improvements.

preprint2009arXiv

Time Evolution of the 3-D Accretion Flows: Effects of the Adiabatic Index and Outer Boundary Condition

We study a slightly rotating accretion flow onto a black hole, using the fully three dimensional (3-D)numerical simulations. We consider hydrodynamics of an inviscid flow, assuming a spherically symmetric density distribution at the outer boundary and a small, latitude-dependent angular momentum. We investigate the role of the adiabatic index and gas temperature, and the flow behaviour due to non-axisymmetric effects. Our 3-D simulations confirm axisymmetric results: the material that has too much angular momentum to be accreted forms a thick torus near the equator and the mass accretion rate is lower than the Bondi rate. In our previous study of the 3-D accretion flows, for gamma=5/3, we found that the inner torus precessed, even for axisymmetric conditions at large radii. The present study shows that the inner torus precesses also for other values of the adiabatic index: gamma=4/3, 1.2 and 1.01. However, the time for the precession to set increases with decreasing gamma. In particular, for gamma=1.01 we find that depending on the outer boundary conditions, the torus may shrink substantially due to the strong inflow of the non-rotating matter and the precession will have insufficient time to develop. On the other hand, if the torus is supplied by the continuous inflow of the rotating material from the outer radii, its inner parts will eventually tilt and precess, as it was for the larger gamma's.

preprint2006arXiv

The pattern of accretion flow onto Sgr A*

The material accreting onto Sgr A* most probably comes from the nearby stars. We analyze the pattern of this flow at distances of a fraction of a parsec and we argue that the net angular momentum of this material is low but non-negligible, and the initially supersonic disk accretion changes into subsonic flow with constant angular momentum. Next we estimate the flow parameters at a distance $R_{BHL}$ from the black hole and we argue that for the plausible parameter range the accretion flow is non-stationary. The inflow becomes supersonic at distance of $\sim 10^4 R_g$ but the solution does not continue below the horizon and the material piles up forming a torus, or a ring, at a distance of a few up to tens of Schwarzchild radii. Such a torus is known to be unstable and may explain strong variability of the flow in Sgr A*. Our considerations show that the temporary formation of such a torus seems to be unavoidable. Our best fitting model predicts a rather large accretion rate of around $4 \cdot 10^{-6} M_{\odot}/yr$ directly on Sgr A*. We argue that magnetic fields in the flow are tangled and this allows our model to overcome the disagreement with the Faraday rotation limits.