The collapse of a massive star with low angular momentum content is commonly argued to result in the formation of a black hole without an accompanying bright transient. Our goal in this Letter is to understand the flow in and around a newly-formed black hole, involving accretion and rotation, via general relativistic hydrodynamics simulations aimed at studying the conditions under which infalling material can accrete without forming a centrifugally supported structure and, as a result, generate no effective feedback. If the feedback from the black hole is, on the other hand, significant, the collapse would be halted and we suggest that the event is likely to be followed by a bright transient. We find that feedback is only efficient if the specific angular momentum of the infalling material at the innermost stable circular orbit exceeds that of geodesic circular flow at that radius by at least $\approx 20\%$. We use the results of our simulations to constrain the maximal stellar rotation rates of the disappearing massive progenitors PHL293B-LBV and N6946-BH1, and to provide an estimate of the overall rate of disappearing massive stars. We find that about a few percent of single O-type stars with measured rotational velocities are expected to spin below the critical value before collapse and are thus predicted to vanish without a trace.

We prescribe a method to study the effects of self-gravity of accretion disk around a black hole associated with long Gamma Ray Bursts (GRBs) in an evolving background Kerr metric. This is an extension to our previous work where we presented possible constraints for the final masses and spins of these astrophysical black holes. Incorporating the self-force of the accreting cloud around the black hole is a very important aspect due to the transient nature of the event, in which a huge amount of mass is accreted and changes the fundamental black hole parameters i.e. its mass and spin, during the process. Understanding of the GRBs engine is important because they are possible sources of high-energy particles and gravitational waves as most of the energy released from the dynamical evolution is in the form of gravitational radiation. Here, we describe the analytical framework we developed to employ in our numerical model. The numerical studies are planned for the future work.

We present the results of simulations done with the code HARM-COOL developed in the CTP PAS Warsaw research group over the years 2017-2019. It is based in the original GR MHD scheme proposed by Gammie et al. (2003) for the simulation of Active Galactic Nucleus, but now it has been suited for the engine of a short Gamma Ray Burst event. We compute time-dependent evolution of a black hole accretion disk, in two-dimensional, axisymmetric scheme. The code includes neutrino cooling and accounts for nuclear structure of dense, degenerate matter. Free protons, neutrons, and electron-positron pairs form a neutron-rich, magnetically driven outflow that provides site for subsequent r-process nucleosynthesis. Here the heavy elements up to the Uranium and Gold are synthesized and may contribute to the chemical enrichment of the circum-burst medium. Their radio-active decay will give signal in lower energies in a timescale of weeks-months after the GRB prompt phase. In addition, the magnetic fields are responsible for the launching of ultra-relativistic jets along the rotation axis of the central black hole, according to the well-known Blandford-Znajek mechanism. These jets are sites of variable high energy emission in gamma rays. We find that the magnetic field and the black hole spin account for the observed variability timescales and jet energetics.

Using axisymetric general relativistic magnetohydrodynamics simulations we study evolution of accretion torus around black hole endowed with different initial magnetic field configurations. Due to accretion of material onto black hole, parabolic magnetic field will develop in accretion torus funnel around vertical axis, for any initial magnetic field configuration.

We present the relativistic 2D and 3D GRMHD simulation of axisymmetric, inviscid, hydrodynamic accretion flows in a fixed Kerr black hole gravitational field. The flow is having low angular momentum with respect to Keplerian one. A relativistic fluid where its bulk velocity is comparable to the speed of light, flowing in the accretion disk very close to the horizon should be described by adiabtic index: 4=3 < g < 5=3 .The time dependent evolution of shock position and respective effect on mass accretion rate and oscillation frequency with varying adiabatic index has been studied. Here we present some of the results for adiabatic index = 1.4 in a 2D and 3D model.

We investigate the collapsar scenario, in which, according to Woosley (1993) and Paczynski (1998), the long gamma ray bursts (LGRBs) can be explained. The total energetics of explosions observed in the gamma ray band are consistent here with the total binding energy of a pogenitor star. However, the details characteristic for various events, such as their individual duration times, lightcurve profiles, variability, and connection with supernovae of different spectral properties, are still subject of many studies. In our scenario, the evolved progenitor star is collapsing onto a compact core, which forms a spinning black hole. The process of accretion via rotationally supported torus powers the ejection of relativistic jets. Here, the rotational energy of the black hole is presumably transported to the remote jet, which is mediated by magnetic fields. The rotation of pre-supernova star is a key property of the model. In our study, we investigate different types of collapsing stars: magnetized and non-magnetized. We probe the distributions of angular momentum inside the collapsar and the process of spinning up the black hole. Both Kerr parameter, and increasing mass of the black hole, put a constraint for the existence of rotationally supported torus inside the collapsar's envelope.

We compute the evolution of a quasi-spherical, slowly rotating accretion flow around a black hole, whose mass and spin evolve adequately to the mass-energy transfer through the horizon. Our model is relevant for the central engine driving a long gamma ray burst, that originates from the collapse of a massive star. Our results show how much mass and spin a newly formed black hole should possess during collapsar to launch long GRB.

I present the publicly available code GLADIS (GLobal Accretion Disk Instability Simulation) developed in my reserach group over the years 2002-2017. It can be freely downloaded and modified by the users via the link from the Astrophysics Source Code Library. The software computes time-dependent evolution of a black hole accretion disk, in one-dimensional, axisymmetric, vertically integrated scheme. The main applications are to explain the variability of accretion disks that can be subject to radiation-pressure instability. The phenomenon is relevant for fast variable microquasars, as well as for a class of changing-look AGN.

Short gamma-ray bursts require a rotating black hole, surrounded by a magnetized relativistic accretion disk, such as the one formed by coalescing binary neutron stars or neutron star - black hole systems. The accretion onto a Kerr black hole is the mechanism of launching a baryon-free relativistic jet. An additional uncollimated outflow, consisting of sub-relativistic neutron-rich material which becomes unbound by thermal, magnetic and viscous forces, is responsible for blue and red kilonova. We explore the formation, composition and geometry of the secondary outflow by means of simulating accretion disks with relativistic magneto-hydrodynamics and employing realistic nuclear equation of state. We calculate the nucleosynthetic r-process yields by sampling the outflow with a dense set of tracer particles. Nuclear heating from the residual r-process radioactivities in the freshly synthesized nuclei is expected to power a red kilonova, contributing independently from the dynamical ejecta component, launched at the time of merger, and neutron-poor broad polar outflow, launched from the surface of the hypermassive neutron star by neutrino wind. Our simulations show that both magnetisation of the disk and high black hole spin are able to launch fast wind outflows ($v/c\sim 0.11-0.23$) with a broad range of electron fraction $Y_{\rm e}\sim 0.1-0.4$, and help explain the multiple components observed in the kilonova lightcurves. The total mass loss from the post-merger disk via unbound outflows is between 2\% and 17\% of the initial disk mass.

We study the role of adiabatic index in determining the critical points in the transonic low angular momentum accretion flow onto a black hole. We present the general relativistic 2D hydrodynamic simulations of axisymmetric, inviscid accretion flows in a fixed Kerr black hole gravitational field. A relativistic fluid where its bulk velocity is comparable to the speed of light, flowing in the accretion disk very close to the horizon can be described by an adiabatic index of 4/3 < {\gamma} < 5/3. The time dependent evolution of the shock position and respective effects on mass accretion rate and oscillation frequency with varying adiabatic index is discussed in the context of the observed microquasars.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Short gamma ray bursts are presumably results of binary neutron star mergers, which lead to the formation of a stellar mass black hole, surrounded by a remnant matter. The strong magnetic fields help collimate jets of plasma, launched along the axis of the black hole rotation. We study the structure and evolution of the accreting plasma in the short GRBs and we model the formation of the base of a relativistic, Poynting-dominated jets. Our numerical models are based on the general relativistic MHD, axisymmetric simulations. We discuss the origin of variability in the GRB jet emission, which timescales are related to the action of the magneto-rotational instability in the accreting plasma. We also estimate the value of a maximum achievable Lorentz factor in the jets produced by our simulations, and reached at the large distances, where the gamma ray emission is produced.

We compute the evolution of a quasi-spherical, slowly rotating accretion flow around a black hole, whose mass and spin evolve adequately to the mass-energy transfer through the horizon. Our model is relevant for the central engine driving a long gamma ray burst, that originates from the collapse of a massive star. The computations of a GRB engine in a dynamically evolving spacetime metric are important specifically due to the transient nature of the event, in which a huge amount of mass is accreted and changes the fundamental black hole parameters, its mass and spin, during the process. We discuss the results in the context of angular momentum magnitude of the collapsing star. We also study the possible formation and evolution of shocks in the envelope, which may temporarily affect accretion. Our results are important for the limitations on the mass and spin range of black holes detected independently by electromagnetic observations of GRBs and gravitational waves. We speculate on the possible constraints for the final masses and spins of these astrophysical black holes. It is shown that the most massive BHs were rahter not formed in a powerful GRB explosion, if the cores of their progenitors were only weakly rotating.

The X-ray Integral Field Unit (X-IFU) is the high resolution X-ray spectrometer of the ESA Athena X-ray observatory. Over a field of view of 5' equivalent diameter, it will deliver X-ray spectra from 0.2 to 12 keV with a spectral resolution of 2.5 eV up to 7 keV on ~5 arcsecond pixels. The X-IFU is based on a large format array of super-conducting molybdenum-gold Transition Edge Sensors cooled at about 90 mK, each coupled with an absorber made of gold and bismuth with a pitch of 249 microns. A cryogenic anti-coincidence detector located underneath the prime TES array enables the non X-ray background to be reduced. A bath temperature of about 50 mK is obtained by a series of mechanical coolers combining 15K Pulse Tubes, 4K and 2K Joule-Thomson coolers which pre-cool a sub Kelvin cooler made of a 3He sorption cooler coupled with an Adiabatic Demagnetization Refrigerator. Frequency domain multiplexing enables to read out 40 pixels in one single channel. A photon interacting with an absorber leads to a current pulse, amplified by the readout electronics and whose shape is reconstructed on board to recover its energy with high accuracy. The defocusing capability offered by the Athena movable mirror assembly enables the X-IFU to observe the brightest X-ray sources of the sky (up to Crab-like intensities) by spreading the telescope point spread function over hundreds of pixels. Thus the X-IFU delivers low pile-up, high throughput (>50%), and typically 10 eV spectral resolution at 1 Crab intensities, i.e. a factor of 10 or more better than Silicon based X-ray detectors. In this paper, the current X-IFU baseline is presented, together with an assessment of its anticipated performance in terms of spectral resolution, background, and count rate capability. The X-IFU baseline configuration will be subject to a preliminary requirement review that is scheduled at the end of 2018.

We model the structure and evolution of black hole accretion disks, and their neighboring regions, using numerical simulations. The numerics is governed by the equations of general relativistic magneto-hydrodynamics (GRMHD). In particular, such disks and outflows can be found at the base of very energetic ultra-relativistic jets produced by cosmic explosions, so called gamma-ray bursts (GRBs). Another, more persistent type of the jet phenomena, are blazars, emitted from the centers of galaxies. Long-lasting, detailed computations are essential to properly determine the physics of these explosions, and confront the theoretical models with any potential observables. From the point of view of numerical methods and computational techniques, three ingredients need to be considered. First, the numerical scheme must work in a conservative manner, which is achieved by solving a set of non-linear equations at each time-step, to advance the conserved quantities from one time step to the next. Second, the efficiency of computations intrinsically depends on the code parallelization methods, which may use various techniques. Third, the analysis of results is possible via the post-processing of the computed time-dependent physical quantities, and visualization of the flow properties. This is done via implementing various packages and libraries that are standardized in the field of computational astrophysics and supported by community developers. In the present paper, we discuss the physical picture of the cosmic sources which are modeled using numerical framework. We also describe several technical issues, in the particular context of our own experience with the performance of the GRMHD code which we develop. We also present a suite of performance tests, done on the High-Performance Computer cluster (HPC) in the Center for Mathematical Modeling of the Warsaw University.

The collapsar model was proposed to explain the long-duration gamma-ray bursts (GRBs), while the short GRBs are associated with the mergers of compact objects. In the first case, mainly the energetics of the events is consistent with the proposed progenitor models, while the duration, time variability, as well as the afterglow emission may shed some light on the detailed properties of the collapsing massive stars. In the latter case, the recent discovery of the binary neutron star (NS-NS) merger in the gravitational wave observation made by LIGO (GW170817), and the detection of associated electromagnetic counterparts, for the first time gave a direct proof of the NS-NS merger being a progenitor of a short GRB. In general, all GRBs are believed to be powered by accretion through a rotationally supported torus, or by fast rotation of a compact object. For long ones, the rotation of the progenitor star is a key property in order to support accretion over relatively long activity periods, and also to sustain the rotation of the black hole itself. The latter is responsible for ejection of the relativistic jets, which are powered due to the extraction of the BH rotational energy, mitigated by the accretion torus and magnetic fields. The jets must break through the stellar envelope though, which poses a question on the efficiency of this process. Similar mechanisms of powering the jet ejection may act in short GRBs, which in this case may freely propagate through the interstellar medium. The power of the jets launched from the rotating black hole is at first associated mostly with the magnetic Poynting flux, and then at large distances it is transferred to the kinetic and finally radiative energy of the expanding shells.

We study the structure and evolution of the accreting plasma in gamma ray burst central engines. The models are based on the general relativistic MHD simulations. The nuclear equation of state adequate for dense and degenerate plasma, is incorporated to the numerical scheme instead of a simple polytropic EOS. Plasma is cooled by neutrinos and energy is extracted from the rotating BH by magnetic fields. We discuss our results in the frame of the observable GRBs, and speculate about the origin of variability in the GRB jet emission. We also discuss the possible contribution from the heavy elements synthesized in winds from GRB engines to the kilonova emission.

We perform 1D/2D/3D relativistic hydrodynamical simulations of accretion flows with low angular momentum, filling the gap between spherically symmetric Bondi accretion and disc-like accretion flows. Scenarios with different directional distributions of angular momentum of falling matter and varying values of key parameters such as spin of central black hole, energy and angular momentum of matter are considered. In some of the scenarios the shock front is formed. We identify ranges of parameters for which the shock after formation moves towards or outwards the central black hole or the long lasting oscillating shock is observed. The frequencies of oscillations of shock positions which can cause flaring in mass accretion rate are extracted. The results are scalable with mass of central black hole and can be compared to the quasi-periodic oscillations of selected microquasars (such as GRS 1915+105, XTE J1550-564 or IGR J17091-3624), as well as to the supermassive black holes in the centres of weakly active galaxies, such as Sgr $A^{*}$.

The thermal stability of accretion disk and the possibility to see a limit-cycle behaviour strongly depends on the ability of the disk plasma to cool down. Various processes connected with radiation-matter interaction appearing in hot accretion disk plasma contribute to opacity. For the case of geometrically thin and optically thick accretion disk, we can estimate the influence of several different components of function \kappa, given by the Roseland mean. In the case of high temperatures, the electron Thomson scattering is dominant. At lower temperatures atomic processes become important. The slope d log \kappa/d log T can have locally stabilizing or destabilizing effect on the disk. Although the local MHD simulation postulate the stabilizing influence of the atomic processes, only the global time-dependent model can reveal the global disk stability range estimation. This is due to global diffusive nature of that processes. In this paper, using previously tested GLADIS code with modified prescription of the viscous dissipation, we examine the stabilizing effect of the Iron Opacity Bump.

Black holes (BHs) surrounded by accretion disks are present in the Universe at different scales of masses, from microquasars up to the active galactic nuclei (AGNs). The current picture remains ad hoc due to the complexity of the magnetic field action. In addition, accretion disks at high Eddington rates can be radiation-pressure dominated and, according to some of the heating prescriptions, thermally unstable. The observational verification of their resulting variability patterns may shed light on both the role of radiation pressure and magnetic fields. We compute the structure and time evolution of an accretion disk, using the code GLADIS (which models the global accretion disk instability). We supplement this model with a modified viscosity prescription, which can to some extent describe the magnetisation of the disk. We study the results for a large grid of models, to cover the whole parameter space, and we derive conclusions separately for different scales of black hole masses, which are characteristic for various types of cosmic sources. We show the dependencies between the flare or outburst duration, its amplitude, and period, on the accretion rate and viscosity scaling. We show that if the heating rate in the accretion disk grows more rapidly with the total pressure and temperature, the instability results in longer and sharper flares. In general, we confirm that the disks around the supermassive black holes are more radiation-pressure dominated and present relatively brighter bursts. Our method can also be used as an independent tool for the black hole mass determination, which we confront now for the intermediate black hole in the source HLX-1.

We calculate the structure and evolution of a gamma ray burst central engine where an accreting torus has formed around the newly born black hole. We study the general relativistic,MHD models and we self-consistently incorporate the nuclear equation of state. The latter accounts for the degeneracy of relativistic electrons, protons, and neutrons, and is used in the dynamical simulation, instead of a standard polytropic -law. The EOS provides the conditions for the nuclear pressure in the function of density and temperature, which evolve with time according to the conservative MHD scheme. We analyze the structure of the torus and outflowing winds, and compute the neutrino flux emitted through the nuclear reactions balance in the dense and hot matter. We also estimate the rate of transfer of the black hole rotational energy to the bipolar jets. Finally, we elaborate on the nucleosynthesis of heavy elements in the accretion flow and the wind, through computations of the thermonuclear reaction network. We discuss the possible signatures of the radioactive elements decay in the accretion flow. We suggest that further detailed modeling of the accretion flow in GRB engine, together with its microphysics, may be a valuable tool to constrain the black hole mass and spin. It can be complementary to the gravitational wave analysis, if the waves are detected with an electromagnetic counterpart.

We study the structure and evolution of the hyperaccreting disks and outflows in the gamma ray bursts central engines. The torus around a stellar mass black hole is composed of free nucleons, Helium, electron-positron pairs, and is cooled by neutrino emission. Accretion of matter powers the relativistic jets, responsible for the gamma ray prompt emission. The significant number density of neutrons in the disk and outflowing material will cause subsequent formation of heavier nuclei. We study the process of nucleosynthesis and its possible observational consequences. We also apply our scenario to the recent observation of the gravitational wave signal, detected on September 14th, 2015 by the two Advanced LIGO detectors, and related to an inspiral and merger of a binary black hole system. A gamma ray burst that could possibly be related with the GW150914 event was observed by the Fermi satellite. It had a duration of about 1 second and appeared about 0.4 seconds after the gravitational-wave signal. We propose that a collapsing massive star and a black hole in a close binary could lead to the event. The gamma ray burst was powered by a weak neutrino flux produced in the star remnant's matter. Low spin and kick velocity of the merged black hole are reproduced in our simulations. Coincident gravitational-wave emission originates from the merger of the collapsed core and the companion black hole.

The hyperluminous X-ray source (HLX-1, the peak X-ray luminosity $\sim 10^{42}\rm erg\ s^{-1}$) near the spiral galaxy ESO 243-49 is possibly the best candidate for intermediate mass black hole (IMBH), which underwent recurrent outbursts with a period of $\sim 400$ days. The physical reason for this quasi-periodic variability is still unclear. We explore the possibility of radiation-pressure instability in accretion disk by modeling the light curve of HLX-1, and find that it can roughly reproduce the duration, period and amplitude of the recurrent outbursts HLX-1 with an IMBH of ~10^5Msun. Our result provides a possible mechanism to explain the recurrent outbursts in HLX-1. We further find a universal correlation between the outburst duration and the bolometric luminosity for the BH sources with a very broad mass range (e.g., X-ray binaries, XRBs, HLX-1 and active galactic nuclei, AGNs), which is roughly consistent with the prediction of radiation-pressure instability of the accretion disk. These results imply that "heartbeat" oscillations triggered by radiation-pressure instability may appears in different-scale BH systems.

We consider the gamma ray burst (GRB) central engine that is powered by the collapse of a massive rotating star or compact binary merger. The engine is a hot and dense accretion disk, which is composed of free nucleons, electron-positron pairs, and helium, and cooled by neutrino emission. A significant number density of neutrons in the inner disk body provide conditions for neutron rich plasma in the GRB outflows or jets. Helium is synthesized in the inner disk if the accretion rate is large, and heavy nuclei are also formed in the outer disk at distances above 150-250 $r_{g}$ from the black hole. We study the process of nucleosynthesis in the GRB engine, depending on its physical properties. The GRB central engine is hydrodynamically modeled in the frame of a dense and hot disk which accretes with a high rate (up to 1 Solar mass per second) onto a maximally spinning, stellar mass black hole. The synthesis of heavy nuclei up to germanium and gallium is then followed by the nuclear reaction network. The accretion at high rate onto a Kerr black hole feeds the engine activity and establishes conditions for efficient synthesis of heavy nuclei in the disk. These processes may have important observational implications for the jet deceleration process and heavy elements observed in the spectra of GRB afterglows.

During the bright outburst in 2011, the black hole candidate IGR J17091-3624 exhibited strong quasi-periodic flare-like events (on timescales of tens of seconds) in some characteristic states, the so-called heartbeat state. From the theoretical point of view, these oscillations may be modeled by the process of accretion disk instability, driven by the dominant radiation pressure and enhanced heating of the plasma. Although the mean accretion rate in this source is probably below the Eddington limit, the oscillations will still have large amplitudes. As the observations show, the source can exhibit strong wind outflow during the soft state. This wind may help to partially or even completely stabilize the heartbeat. Using our hydrodynamical code GLADIS, we modeled the evolution of an accretion disk responsible for X-ray emission of the source. We accounted for a variable wind outflow from the disk surface. We examined the data archive from the Chandra and XMM-Newton satellites to find the observed limitations on the wind physical properties, such as its velocity and ionization state. We also investigated the long-term evolution of this source, which lasted over about 600 days of observations, using the data collected by the Swift and RXTE satellites. During this long period, the oscillations pattern and the observable wind properties changed systematically. We found that this source probably exhibits observable outbursts of appropriate timescales and amplitudes as a result of the disk instability. Our model requires a substantial wind component to explain the proper variability pattern, and even complete suppression of flares in some states. The wind mass-loss rate extracted from the data agrees quantitatively well with our scenario.

The microquasar IGR J17091-3624, which is the recently discovered analogue of the well known source GRS 1915+105, exhibits quasi-periodic outbursts, with a period of 5-70 seconds, and regular amplitudes, referred to as "heartbeat state". We argue that these states are plausibly explained by accretion disk instability, driven by the dominant radiation pressure. Using our GLobal Accretion DIsk Simulation hydrodynamical code, we model these outbursts quantitatively. We also find a correlation between the presence of massive outflows launched from the accretion disk and the stabilization of its oscillations. We verify the theoretical predictions with the available timing and spectral observations. Furthermore, we postulate that the underlying non-linear differential equations that govern the evolution of an accretion disk are responsible for the variability pattern of several other microquasars, including XTE J1550-564, GX 339-4, and GRO J1655-40. This is based on the signatures of deterministic chaos in the observed lightcurves of these sources, which we found using the recurrence analysis method. We discuss these results in the frame of the accretion disk instability model.

The X-ray Integral Field Unit (X-IFU) on board the Advanced Telescope for High-ENergy Astrophysics (Athena) will provide spatially resolved high-resolution X-ray spectroscopy from 0.2 to 12 keV, with 5 arc second pixels over a field of view of 5 arc minute equivalent diameter and a spectral resolution of 2.5 eV up to 7 keV. In this paper, we first review the core scientific objectives of Athena, driving the main performance parameters of the X-IFU, namely the spectral resolution, the field of view, the effective area, the count rate capabilities, the instrumental background. We also illustrate the breakthrough potential of the X-IFU for some observatory science goals. Then we briefly describe the X-IFU design as defined at the time of the mission consolidation review concluded in May 2016, and report on its predicted performance. Finally, we discuss some options to improve the instrument performance while not increasing its complexity and resource demands (e.g. count rate capability, spectral resolution). The X-IFU will be provided by an international consortium led by France, The Netherlands and Italy, with further ESA member state contributions from Belgium, Finland, Germany, Poland, Spain, Switzerland and two international partners from the United States and Japan.

Data from the Fermi Gamma-ray Burst Monitor satellite observatory suggested that the recently discovered gravitational wave source, a pair of two coalescing black holes, was related to a gamma-ray burst. The observed high-energy electromagnetic radiation (above 50 keV) originated from a weak transient source and lasted for about 1 second. Its localization is consistent with the direction to GW150914. We speculate about the possible scenario for the formation of a gamma-ray burst accompanied by the gravitational-wave signal. Our model invokes a tight binary system consisting of a massive star and a black hole which leads to the triggering of a collapse of the star's nucleus, the formation of a second black hole, and finally to the binary black hole merger. For the most-likely configuration of the binary spin vectors with respect to the orbital angular momentum in the GW150914 event, the recoil speed (kick velocity) acquired by the final black hole through gravitational wave emission is of the order of a few hundred km/s and this might be sufficient to get it closer to the envelope of surrounding material and capture a small fraction of matter from the remnant of the host star. The gamma-ray burst is produced by the accretion of this remnant matter onto the final black hole. The moderate spin of the final black hole suggests that the gamma-ray burst jet is powered by weak neutrino emission rather than the Blandford-Znajek mechanism, and hence explains the low power available for the observed GRB signal.

We study the X-ray emission of the microquasar XTE J1550-564 and analyze the properties of its light curves using the recurrence analysis method. The indicators for non-linear dynamics of the accretion flow are found in the very high state and soft state of this source. The significance of deterministic variability depends on the energy band. We discuss the non-linear dynamics of the accretion flow in the context of the disc-corona geometry and propagating oscillations in the accretion flow.

We describe the methods and results of numerical simulations of coalescing black holes. The simulation in dynamical spacetime covers the inspiral, merger, and ringdown phases. We analyze the emission of gravitational waves and properties of a black hole being the merger product. We discuss the results in the context of astrophysical environment of black holes that exist in the Universe.

The high energy radiation emitted by black hole X-ray binaries originates in an accretion disk, hence the variability of the lightcurves mirrors the dynamics of the disc. We study the time evolution of the emitted flux in order to find evidences, that low dimensional non-linear equations govern the accretion flow. Here we test the capabilities of our novel method to find chaotic behaviour on the two numerical time series describing the motion of a test particle around a black hole surrounded by a thin massive disc, one being regular and the other one chaotic.

We study the low angular momentum accretion of matter onto Schwarzschild or Kerr black hole using fully relativistic numerical simulations. Our aim is to include the rotation of both the gas and the black hole and the presence of the magnetic field into the picture. We compare resulting structure of transonic flow with results of 1D pseudo-Newtonian semi-analytical computations of non-magnetized flow.

The black hole candidates exhibit fast variability of their X-ray emission on a wide range of timescales. The short, coherent variations, with frequencies above 1 Hz, are referred to as quasi-periodic oscillations, and are generally explained by resonant effects in the black hole accretion flow. The purely stochastic variability that occurs due to turbulent conditions in the plasma, is quantified by the power density spectra and appears practically in all types of sources and their spectral states. The specific kind of quasi-periodic flares is expected, when the global structure of the accretion flow, governed by the nonlinear hydrodynamics, induces fluctuations around a fixed point solution. These conditions, which occur at high accretion rates, lead to the variability in the sense of deterministic chaos. We study the nonlinear behaviour of X-ray sources using the recurrence analysis method. We estimate quantitatively the indications for deterministic chaos, such as the Renyi's entropy, for the observed time series, and we compare them with the surrogate data. Using the data collected by RXTE satellite, we reveal the oscillations pattern and the observable properties of six black hole systems. We confirm that both the well known microquasar GRS 1915+105, as well as its recently discovered analogue, IGR J17091-3624, exhibit variability characteristic to deterministic chaotic system. We conclude therefore that the underlying nature of the process must be intrinsically connected in these sources with the accretion flow instability, that leads to the limit cycle oscillations around a fixed point. Furthermore, we studied four other black hole candidates with high accretion rates. For GX 339-4, GRO J1655-40 and XTE J1550-564 we obtained significant results for some observations, whereas for XTE J1650-500 no such data set was found.

Hardly any of the observed black hole accretion disks in X-Ray binaries and active galaxies shows constant flux. When the local stochastic variations of the disk occur at specific regions where a resonant behaviour takes place, there appear the Quasi-Periodic Oscillations (QPOs). If the global structure of the flow and its non-linear hydrodynamics affects the fluctuations, the variability is chaotic in the sense of deterministic chaos. Our aim is to solve a problem of the stochastic versus deterministic nature of the black hole binaries vari- ability. We use both observational and analytic methods. We use the recurrence analysis and we study the occurence of long diagonal lines in the recurrence plot of observed data series and compare it to the sur- rogate series. We analyze here the data of two X-Ray binaries - XTE J1550-564, and GX 339-4 observed by Rossi X-ray Timing Explorer. In these sources, the non-linear variability is expected because of the global conditions (such as the mean accretion rate) leading to the pos- sible instability of an accretion disk. The thermal-viscous instability and fluctuations around the fixed-point solution occurs at high accretion rate, when the radiation pressure gives dominant contribution to the stress tensor.

We present the results of a study of the low-mass X-ray binary 4U 1636-536. We have performed temporal analysis of all available RXTE/ASM, Swift/BAT and MAXI data. We have confirmed the previously discovered quasi-periodicity of ~45 d present during ~2004, however we found it continued to 2006. At other epochs, the quasi-periodicity is only transient, and the quasi-period, if present, drifts. We have then applied a time-dependent accretion disc model to the interval with the significant X-ray quasi-periodicity. For our best model, the period and the amplitude of the theoretical light curve agree well with that observed. The modelled quasi-periodicity is due to the hydrogen thermal-ionization instability occurring in outer regions of the accretion disc. The model parameters are the average mass accretion rate (estimated from the light curves), and the accretion disc viscosity parameters, for the hot and cold phases. Our best model gives relatively low values of viscosity parameter for cold phase 0.01 and for hot phase 0.03.

The ultrarelativistic jets responsible for prompt and afterglow emission in gamma ray bursts are presumably driven by a central engine that consists of a dense accretion disk around a spinning black hole. We consider such engine, composed of free nucleons, electron-positron pairs, Helium nuclei, and cooled by neutrino emission. A significant number density of neutrons in the disk provide conditions for neutron rich plasma in the outflows and jets. Heavy nuclei are also formed in the accretion flow, at the distances 150-250 gravitational radii from the black hole. We study the process of nucleosynthesis in the GRB engine, depending on its physical properties. Our results may have important observational implications for the jet deceleration process and heavy elements observed in the spectra of GRB afterglows.

The X-ray weakness of broad absorption line (BAL) quasars in comparison to non-BAL objects is possibly caused by the absorption of X-ray emission by the shielding material near the equatorial plane. On the other hand, the radio-loud BALQSOs are more X-ray loud than the radio-quiet ones. This suggests that part of the X-ray emission may arise from the radio jet. To investigate this possibility, we modelled the nuclear spectra of two radio-loud BALQSOs. We focus on the emission from the very centres of these two objects. The source of emission was approximated by a single, homogeneous component that produces synchrotron and inverse-Compton (IC) radiation. The simplicity of the model allowed us to estimate the basic physical parameters of the emitting regions, using a universal analytic approach. Such methods have already been used in blazars. For the first time we propose this solution for quasars. In addition, we modelled the radiation spectra of the accretion disk and its corona to compare them with the jets' spectra. We find that in the case of luminous 3C270.1, the nuclear X-ray continuum is dominated by the non-thermal, IC emission from the innermost parts of the radio jet. However, the radio core of the lobe-dominated PG 1004+130 is probably too weak to produce significant part of the observed X-ray emission. A large contribution from the X-ray emitting accretion disk and corona is produced in our model. However, it then exceeds the observed flux. Because the large intrinsic absorption was postulated recently by the NuSTAR observations, we propose that the disk-corona component may still account for the X-rays produced in this source. This part of the spectrum must nevertheless be dominated by the X-ray jet. We conclude that the jet-linked X-ray emission is present in strong and weak radio sources, but its fraction seems to scale with the radio jet power.

We derive the conditions for shock formation in a quasi-spherical, slightly rotating flows. We verify the results of semi-analytical, stationary calculations with the time evolution studied by numerical hydro-simulations, and we study the oscillations of the shock position. We also study the behaviour of flows with varying specific angular momentum, where the 'hysteresis' type of loop is found when passing through the multiple sonic points region. Our results are in agreement with the timescales and shapes of the luminosity flares observed in Sgr A*. These models may also be applicable for the Galactic stellar mass black holes, like GX 339-4 or GRS 1915+105, where periodic oscillations of X-ray luminosity are detected.

We consider a scenario for the longest duration gamma ray bursts, resulting from the collapse of a massive rotating star in a close binary system with a companion black hole. The primary black hole born during the core collapse is first being spun up and increases its mass during the fallback of the stellar envelope just after its birth. As the companion black hole enters the outer envelope, it provides an additional angular momentum to the gas. After the infall and spiral-in towards the primary, the two black holes merge inside the circumbinary disk. The second episode of mass accretion and high final spin of the post-merger black hole prolongs the gamma ray burst central engine activity. The observed events should have two distinct peaks in the electromagnetic signal, separated by the gravitational wave emission. The gravitational recoil of the burst engine is also possible.

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.

We consider a scenario for the longest duration gamma ray bursts, resulting from the collapse of a massive star in a close binary system with a companion black hole. The primary black hole born during the core collapse is spun up and increases its mass during the fallback of the stellar envelope. The companion black hole provides an additional angular momentum to the envelope, which ultimately makes the core BH spinning with a high Kerr parameter. After the infall and spiral-in, the two black holes merge inside the circumbinary disk. The second episode of mass accretion and final, even larger spin of the post-merger black hole prolongs the gamma ray burst central engine activity. The observed events should have two distinct peaks in the electromagnetic signal, separated by the gravitational wave emission. The gravitational recoil of the burst engine is also possible.

We discuss the issues of stability of accretion disks that may undergo the limit-cycle oscillations due to the two main types of thermal-viscous instabilities. These are induced either by the domination of radiation pressure in the innermost regions close to the central black hole, or by the partial ionization of hydrogen in the zone of appropriate temperatures. These physical processes may lead to the intermittent activity in AGN on timescales between hundreds and millions of years. We list a number of observational facts that support the idea of the cyclic activity in high accretion rate sources. We conclude however that the observed features of quasars may provide only indirect signatures of the underlying instabilities. Also, the support from the sources with stellar mass black holes, whose variability timescales are observationally feasible, is limited to a few cases of the microquasars. Therefore we consider a number of plausible mechanisms of stabilization of the limit cycle oscillations in high accretion rate accretion disks. The newly found is the stabilizing effect of the stochastic viscosity fluctuations.

Gamma-ray catalogs contain a considerable amount of unidentified sources. Many of these are located out of the Galactic plane and therefore may have extragalactic origin. Here we assume that the formation of massive black holes in galactic nuclei proceeds through a quasi-star stage and consider the possibility of jet production by such objects. Those jets would be the sources of collimated synchrotron and Compton emission, extending from radio to gamma rays. The expected lifetimes of quasi-stars are of the order of million of years while the jet luminosities, somewhat smaller than that of quasar jets, are sufficient to account for the unidentified gamma-ray sources. The jet emission dominates over the thermal emission of a quasi-star in all energy bands, except when the jet is not directed towards an observer. The predicted synchrotron emission peaks in the IR band, with the flux close to the limits of the available IR all sky surveys. The ratio of the $\gamma$-ray flux to the IR flux is found to be very large ($\sim 60$), much larger than in BL Lac objects but reached by some radio-loud quasars. On the other hand, radio-loud quasars show broad emission lines while no such lines are expected from quasi-stars. Therefore the differentiation between various scenarios accounting for the unidentified gamma-ray sources will be possible at the basis of the photometry and spectroscopy of the IR/optical counterparts.

The standard thin accretion disk model has been successfully used to explain the soft X-ray spectra of Galactic black hole systems and perhaps the UV emission of Active Galactic Nuclei. However, radiation pressure dominated disks are known to be viscously unstable and should produce large amplitude oscillations that are typically not observed. Instead, these sources exhibit stochastic variability which may naturally arise due to viscous fluctuations in a turbulent disk. Here we investigate whether these aperiodic viscous fluctuations can stabilize the inner radiation pressure dominated disks and hence maybe the answer to a forty year old problem in accretion disk theory.

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

We present the state of current research of Very Short Gamma Ray Bursts (VSGRBs) from seven GRB detectors. We found that VSGRBs form distinct class of GRBs, which in our opinion, in most cases can originate from the evaporating Primordial Black Holes (PBHs). Arguments supporting our opinion: 1. GRBs with time duration (T90) < 100 ms form distinct class: VSGRBs. 2. We observe significant anisotropy in the galactic angular distribution of BATSE VSGRB events. 3. V/Vmax distribution for BATSE VSGRB events indicates the local distance production. 4. VSGBBs have more energetic {\gamma}-ray burst than other GRBs with longer duration (KONUS). 5. We observe small number of afterglows in SWIFT VSGRB sample (25%), in contrast with the noticeable afterglow frequency in SGRB sample (78%). 6. Time profile of rising part BATSE VSGRBs is in agreement with the evaporation PBH model.

We report the high resolution radio observations and their analysis of a radio-loud compact steep spectrum (CSS) quasar FIRST J164311.3+315618, one of the members of a binary system. The second component of the system is a radio-quiet AGN. The projected separation of this pair is 2.3 arcsec (15 kpc) and it is one of the known smallest separation binary quasars. The multi-band images of this binary system made with the Hubble Space Telescope showed that the host galaxy of the radio-loud quasar is highly disturbed. The radio observations presented here were made with the multi-element radio linked interferometer network (MERLIN) at 1.66 GHz and 5 GHz. We show that the radio morphology of FIRST J164311.3+315618 is complex on both frequencies and exhibits four components, which indicate on the intermittent activity with a possible rapid change of the jet direction and/or restart of the jet due to the interaction with the companion. The radio components that are no longer powered by the jet can quickly fade away. We suggest that this makes the potential distortions of the radio structure to be short-lived phenomena. On the other hand, our numerical simulations show that the influence of the companion can lead to the prolonged current and future activity. FIRST J164311.3+315618 is an unusual and statistically very rare low redshift binary quasar in which probably the first close encounter is just taking place.

We discuss two important instability mechanisms that may lead to the limit-cycle oscillations of the luminosity of the accretion disks around compact objects: ionization instability and radiation-pressure instability. Ionization instability is well established as a mechanism of X-ray novae eruptions in black hole binary systems but its applicability to AGN is still problematic. Radiation pressure theory has still very weak observational background in any of these sources. In the present paper we attempt to confront the parameter space of these instabilities with the observational data. At the basis of this simple survey of sources properties we argue that the radiation pressure instability is likely to be present in several Galactic sources with the Eddington ratios above 0.15, and in AGN with the Eddington ratio above 0.025. Our results favor the parameterization of the viscosity through the geometrical mean of the radiation and gas pressure both in Galactic sources and AGN. More examples of the quasi-regular outbursts in the timescales of 100 seconds in Galactic sources, and hundreds of years in AGN are needed to formulate firm conclusions. We also show that the disk sizes in the X-ray novae are consistent with the ionization instability. This instability may also considerably influence the lifetime cycle and overall complexity in the supermassive black hole environment.

We consider the radiation pressure instability operating on short timescales 10^3 - 10^6 years in the accretion disk around a supermassive black hole as the origin of the intermittent activity of radio sources. We test whether this instability can be responsible for short ages (<10^4 years) of Compact Steep Spectrum sources measured by hot spots propagation velocities in VLBI observations and statistical overabundance of Gigahertz Peaked Spectrum sources.The implied timescales are consistent with the observed ages of the sources. We aslo discuss possible implications of the intermittent activity on the complex morphology of radio sources, such as the quasar 1045+352, dominated by a knotty jet showing several bends. It is possible that we are whitnessing an ongoing jet precession in this source due to internal instabilities within the jet flow.

Both types of long and short gamma ray bursts involve a stage of a hyper-Eddington accretion of hot and dense plasma torus onto a newly born black hole. The prompt gamma ray emission originates in jets at some distance from this 'central engine' and in most events is rapidly variable, having a form of spikes and subpulses. This indicates at the variable nature of the engine itself, for which a plausible mechanism is an internal instability in the accreting flow. We solve numerically the structure and evolution of the neutrino-cooled torus. We take into account the detailed treatment of the microphysics in the nuclear equation of state that includes the neutrino trapping effect. The models are calculated for both Schwarzschild and Kerr black holes. We find that for sufficiently large accretion rates (> 10 Msun/s for non-rotating black hole, and >1 Msun/s for rotating black hole, depending on its spin), the inner regions of the disk become opaque, while the helium nuclei are being photodissociated. The sudden change of pressure in this region leads to the development of a viscous and thermal instability, and the neutrino pressure acts similarly to the radiation pressure in sub-Eddington disks. In the case of rapidly rotating black holes, the instability is enhanced and appears for much lower accretion rates. We also find the important and possibly further destabilizing role of the energy transfer from the rotating black hole to the torus via the magnetic coupling.

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.

We present new more sensitive high-resolution radio observations of a compact broad absorption line (BAL) quasar, 1045+352, made with the EVN+MERLIN at 5 GHz. They allowed us to trace the connection between the arcsecond structure and the radio core of the quasar. The radio morphology of 1045+352 is dominated by a knotty jet showing several bends. We discuss possible scenarios that could explain such a complex morphology: galaxy merger, accretion disk instability, precession of the jet and jet-cloud interactions. It is possible that we are witnessing an ongoing jet precession in this source due to internal instabilities within the jet flow, however, a dense environment detected in the submillimeter band and an outflowing material suggested by the X-ray absorption could strongly interact with the jet. It is difficult to establish the orientation between the jet axis and the observer in 1045+352 because of the complex structure. Nevertheless taking into account the most recent inner radio structure we conclude that the radio jet is oriented close to the line of sight which can mean that the opening angle of the accretion disk wind can be large in this source. We also suggest that there is no direct correlation between the jet-observer orientation and the possibility of observing BALs.

We discuss a very peculiar subgroup of gamma-ray bursts among the BATSE sources. These bursts are very short ($T_{90} \le $0.1 s), hard, and came predominantly from a restricted direction of the sky (close to the Galactic anti-center). We analyze their arrival times and possible correlations, as well as the profiles of individual bursts. We find no peculiarities in the arrival times of Very Short Bursts (VSBs) despite their highly non-uniform spatial distribution. There is no dependence in the burst shapes on location. Bursts coming both from the burst-enhancement Galactic Anticenter region and from all other directions show considerable dispersion in their rise and fall times. Significant fraction of VSBs have multiple peaks despite their extremely short duration. Burst time properties are most likely to be consistent with two origin mechanisms: either with binary NS-NS mergers with low total masses passing through a phase of hypermassive neutron star, or with evaporation of the primordial black holes in the scenario of no photosphere formation.

Large scale X-ray jets that extend to >100 kpc distances from the host galaxy indicate the importance of jets interactions with the environment on many different physical scales. Morphology of X-ray clusters indicate that the radio-jet activity of a cD galaxy is intermittent. This intermittency might be a result of a feedback and/or interactions between galaxies within the cluster. Here we consider the radiation pressure instability operating on short timescales (<10^5 years) as the origin of the intermittent behaviour. We test whether this instability can be responsible for short ages (< 10^4 years) of Compact Symmetric Objects measured by hot spots propagation velocities in VLBI observations. We model the accretion disk evolution and constrain model parameters that may explain the observed compact radio structures and over-abundance of GPS sources. We also describe effects of consequent outbursts.

We report on the third phase of our study of the neutrino-cooled hyperaccreting torus around a black hole that powers the jet in Gamma Ray Bursts. We focus on the influence of the black hole spin on the properties of the torus. The structure of a stationary torus around the Kerr black hole is solved numerically. We take into account the detailed treatment of the microphysics in the nuclear equation of state that includes the neutrino trapping effect. We find, that in the case of rapidly rotating black holes, the thermal instability discussed in our previous work is enhanced and develops for much lower accretion rates. We also find the important role of the energy transfer from the rotating black hole to the torus, via the magnetic coupling.

We present new results on X-ray properties of radio loud broad absorption line (BAL) quasars and focus on broad-band spectral properties of a high ionization BAL (HiBAL) compact steep spectrum (CSS) radio-loud quasar 1045+352. This HiBAL quasar has a very complex radio morphology indicating either strong interactions between a radio jet and the surrounding interstellar medium or a possible re-start of the jet activity. We detected 1045+352 quasar in a short 5 ksec Chandra ACIS-S observation. We applied theoretical models to explain spectral energy distribution (SED) of 1045+352 and argue that non-thermal, inverse-Compton emission from the innermost parts of the radio jet can account for a large fraction of the observed X-ray emission. In our analysis we also consider a scenario in which the observed X-ray emission from radio-loud BAL quasars can be a sum of inverse-Compton jet X-ray emission and optically thin corona X-ray emission. We compiled a sample of radio-loud BAL quasars that were observed in X-rays to date and report no correlation between their X-ray and radio luminosity. However, the radio-loud BAL quasars show a large range of X-ray luminosities and absorption columns. This is consistent with the results obtained earlier for radio-quiet BAL quasars and may indicate an orientation effect in BAL quasars or more complex dependence between X-ray emission, radio emission and an orientation based on the radio morphology.

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.

We associate the existence of short-lived compact radio sources with the intermittent activity of the central engine caused by a radiation pressure instability within an accretion disk. Such objects may constitute a numerous sub-class of Giga-Hertz Peaked Spectrum sources, in accordance with the population studies of radio-loud active galaxies, as well as detailed investigations of their radio morphologies. We perform the model computations assuming the viscosity parametrization as proportional to a geometrical mean of the total and gas pressure. The implied timescales are consistent with the observed ages of the sources. The duration of an active phase for a moderate accretion rate is short enough (< 10^3-10^4 years) that the ejecta are confined within the host galaxy and thus these sources cannot evolve into large size radio galaxies unless they are close to the Eddington limit.

In the frame of the collapsar model for long gamma ray bursts (GRBs), we investigate the formation of a torus around a spinning BH and we check what rotational properties a progenitor star must have in order to sustain torus accretion over relatively long activity periods. We also study the time evolution of the BH spin parameter. We take into account the coupling between BH mass, its spin parameter and the critical specific angular momentum of accreting gas, needed for the torus to form. The large BH spin reduces the critical angular momentum which in turn can increase the GRB duration with respect to the Schwarzschild BH case. We quantify this effect and estimate the GRB durations in three cases: when a hyper accreting torus operates or a BH spins very fast or both. We show under what conditions a given progenitor star produces a burst that can last as short as several seconds and as long as several hundred of seconds. Our models indicate that it is possible for a single collapse to produce three kinds of jets: (1) a very short, lasting between a fraction of a second and a few seconds, 'precursor' jet, powered only by a hyper accreting torus before the BH spins up, (2) an 'early' jet, lasting several tens of seconds and powered by both hyper accretion and BH rotation, and (3) a 'late' jet, powered only by the spinning BH.

We report on the fourth phase of our study of slightly rotating accretion flows onto black holes. The main new element of this study is that we used fully three dimensional (3-D) numerical simulations. We consider hydrodynamics of inviscid accretion flows. We assume a spherically symmetric density distribution at the outer boundary, but brake the flow symmetry by introducing a small, latitude-dependent angular momentum. We also consider cases where angular momentum at large radii is latitude- and azimuth-dependent. For the latitude-dependent angular momentum, 3-D simulations confirm axisymmetric results: the material that has too much angular momentum to be accreted forms a thick torus near the equator. Consequently, accretion proceeds only through the polar funnel, and the mass accretion rate through the funnel is constrained by the size and shape of the torus, not by the outer conditions. In 3-D simulations, we found that the torus precesses, even for axisymmetric conditions at large radii. For the latitude and azimuth-dependent angular momentum, the non-rotating gas near the equator can also significantly affect the evolution of the rotating gas. In particular, it may prevent the formation of a proper torus (i.e. its closing, in the azimuthal direction). In such models, the mass accretion rate is only slightly less than the corresponding Bondi rate.

We select a sample of 10 radio-quiet quasars with confirmed intranight optical variability and with available X-ray data. We compare the variability properties and the broad band spectral constraints to the predictions of intranight variability by three models: (i) irradiation of an accretion disk by a variable X-ray flux (ii) an accretion disk instability (iii) the presence of a weak blazar component. We concluded that the third model, e.g. the blazar component model, is the most promising if we adopt a cannonball model for the jet variable emission. In this case, the probability of detecting the intranight variability is within 20-80%, depending on the ratio of the disk to the jet optical luminosity. Variable X-ray irradiation mechanism is also possible but only under additional requirement: either the source should have a very narrow Hbeta line or occasional extremely strong flares should appear at very large disk radii.

The collapsar model is the most promising scenario to explain the huge release of energy associated with long duration gamma-ray-bursts (GRBs). Within this scenario GRBs are believed to be powered by accretion through a rotationally support torus or by fast rotation of a compact object. In both cases then, rotation of the progenitor star is one of the key properties because it must be high enough for the torus to form, the compact object to rotate very fast, or both. Here, we check what rotational properties a progenitor star must have in order to sustain torus accretion over relatively long activity periods as observed in most GRBs. We show that simple, often cited, estimates of the total mass available for torus formation and consequently the duration of a GRB are only upper limits. We revise these estimates by taking into account the long term effect that as the compact object accretes the minimum specific angular momentum needed for torus formation increases. This in turn leads to a smaller fraction of the stellar envelope that can form a torus. We demostrate that this effect can lead to a significant, an order of magnidute, reduction of the total energy and overall duration of a GRB event. This of course can be mitigated by assuming that the progenitor star rotates faster then we assumed. However, our assumed rotation is already high compared to observational and theoretical constraints. We also discuss implications of our result.

We show the locations of the SWIFT short hard bursts (SHB) with afterglows on the galactic map and compare with the VSB BATSE events. As we have pointed out before, there is an excess of events in the galactic map of BATSE VSB events. We not that none of VSB SWIFT era events fall into this cluster. More SWIFT events are needed to check this claim. We also report a new study with KONUS data of the VSB sample with an average energy above 90 keV showing a clear excess of events below 100 ms duration (T90) that have large mean energy protons. We suggest that VSB themselves consist of two subclasses: a fraction of events have peculiar distribution properties and have no detectable counterparts, as might be expected for exotic sources such as primordial black holes. We show how GLAST could add key new information to the study of VSB bursts and could help test the black hole concept.

We investigate the properties and evolution of accretion tori formed after the coalescence of two compact objects. At these extreme densities and temperatures, the accreting torus is cooled mainly by neutrino emission produced primarily by electron and positron capture on nucleons (beta reactions). We solve for the disc structure and its time evolution by introducing a detailed treatment of the equation of state which includes photodisintegration of helium, the condition of beta-equilibrium, and neutrino opacities. We self-consistently calculate the chemical equilibrium in the gas consisting of helium, free protons, neutrons and electron-positron pairs and compute the chemical potentials of the species, as well as the electron fraction throughout the disc. We find that, for sufficiently large accretion rates (> 10 solar masses per second), the inner regions of the disk become opaque and develop a viscous and thermal instability. The identification of this instability might be relevant for GRB observations.

We develop a two-flow model of accretion onto a black hole which incorporates the effect of the local magneto-rotational instability. The flow consists of an accretion disk and an accreting corona, and the local dynamo affects the disk/corona mass exchange. The model is aimed to explain the power spectrum density of the sources in their soft, disk-dominated states. The local perturbations of the magnetic field in the disk are described as in King et al. (2004) and Mayer and Pringle (2005), but the time-dependent local magnetic field is assumed to affect the local supply of the material to the corona. Since the viscous timescale in the corona is much shorter than in the disk, the local perturbations are not smeared in the corona. Simple analytical estimates and full time-dependent computations of the disk-corona system are performed. The accreting corona model can reproduce the broad power spectra of Soft State X-ray binaries and AGN. The model, however, predicts that (i) sources undergoing radiation pressure instability (GRS 1915+105) should have systematically steeper power spectra than other sources, (ii) AGN should have systematically steeper power spectra than GBH, even if their disks are described using viscosity proportional to the gas pressure. More precise measurements of power spectra of Soft State sources are clearly needed.

We develop a model of time delays between the continuum bands in the Narrow Line Seyfert 1 galaxy Mrk 335 to explain the observed delays measured in this source. We consider two geometries: an accretion disk with fully ionized warm absorber of considerable optical depth, located close to the symmetry axis, and an accretion disk with a hot corona. Both media lead to significant disk irradiation but the disk/corona geometry gives lower values of the time delays. Only the disk/corona models give results consistent with measurements of Sergeev et al., and a low value of the disk inclination is favored. The presence of an optically thick, fully ionized outflow is ruled out at the 2-sigma level.

We show the locations of the SWIFT short hard bursts (SHB) with afterglows on the Galactic map and compare with the very short bursts (VSB) BATSE events. As we have pointed out before, there is an excess of events in the Galactic map of BATSE VSB events. We note, that none of VSB SWIFT era events fall into this cluster. More SWIFT events are needed to check this claim. We also report a new study with KONUS data of the VSB sample with an average energy above 90 keV showing a clear excess of events below 100 ms duration (T_90) that have large mean energy photons. We suggest that VSB themselves consists of two subclasses: a fraction of events have peculiar distribution properties and have no detectable counter parts, as might be expected for exotic sources such as Primordial Black Holes.

The modified version of a bullet model for gamma ray bursts is studied. The central engine of the source produces multiple sub-jets that are contained within a cone. The emission of photons in the source frame of a sub-jet either takes part in an infinitesimally thin shell, or during its expansion for a finite time. The analysis of the observed profiles of GRBs taken by BATSE leads us to the conclusion that the latter possibility is much more favored. We also study the statistical distribution of GRBs, in the context of their bimodality of durations, taking into account the detector's capability of observing the signal above a certain flux limit. The model with shells emitting for a finite time is able to reproduce only one class of bursts, short or long, depending on the adopted physical parameters. Therefore we suggest that the GRB bimodality is intrinsically connected with two separate classes of sources.

Our recent studies of Very Short Gamma Ray Bursts (VSB) of duration T_{90} \le 100 ms have indicated a significant angular asymmetry and a uniform V/V_{max} distribution from the BATSE data. Here we update these studies, and we extend our research to events observed with KONUS satellite which gives a new insight into the spectra not possible with the BATSE data. KONUS has observed 18 events with T_{90} \le 100 ms duration. These events display considerable numbers of photons with energies above 1 MeV and have photons above 5 MeV in all cases. These appear to be some of the most energetic photons observed from any classes of GRB to date.

At the extreme densities and temperatures typical of the central engine of GRBs, the accreting torus is cooled mainly by advection and by neutrino emission. The latter process is dominated by electron and positron capture onto nucleons ($\beta$ reactions). We calculate the reaction rates and the nuclear composition of matter, assuming that the torus consists of helium, eletron-positron pairs, free neutrons and protons. After determining the equation of state and solving for the disk structure for a given initial accretion rate, we subsequently follow its time evolution. We find that, for accretion rates of the order of $10 M_{\odot}$/s, likely typical for the early stages of the accretion event, the disk becomes unstable, giving rise to variable energy output. This instability may play an important role for producing internal shocks.

The hard X-ray lightcurves exhibit delays of $\sim 1$ s with respect to the soft X-ray lightcurves when the microquasar GRS 1915+105 is in the state of frequent, regular outbursts (states $\rho$ and $\kappa$ of Belloni et al. 2000). Such outbursts are supposed to be driven by the radiation pressure instability of the inner disc parts. The hard X-ray delays are then caused by the time needed for the adjustment of the corona to changing conditions in the underlying disc. We support this claim by the computation of the time evolution of the disc, including a non-stationary evaporation of the disc and mass exchange with the corona.

Rapid, hyper-Eddington accretion is likely to power the central engines of gamma-ray bursts (GRBs). In the extreme conditions of densities and temperatures the accreting torus is cooled by neutrino emission rather than by radiation. Another important cooling mechanism is the advection of energy into the central black hole. We compute the time evolution of a neutrino-dominated disc that proceeds during the burst and investigate the changes in its density and temperature. The discrimination between short and long bursts is made on the basis of the different rates of material inflow to the outer parts of the disc, thus favoring the binary merger scenario for the short GRBs and the collapsar scenario for the long ones. Within the context of the collapsar model, we also study the evolution of the photon luminosity of the remnant disc up to times of the order of 1 day, and discuss its implications for the production of emission lines in GRB spectra.

We consider effects of the MHD turbulence on the viscosity during the evolution of the thermal-viscous ionization instability in the standard $\alpha$-accretion disks. We consider the possibility that the accretion onto a supermassive black hole proceeds through an outer standard accretion disk and inner, radiatively inefficient and advection dominated flow. In this scenario we follow the time evolution of the accretion disk in which the viscosity parameter $\alpha$ is constant throughout the whole instability cycle, as implied by the strength of MHD turbulence. We conclude that the hydrogen ionization instability is a promising mechanism to explain the intermittent activity in AGN.

The time dependent evolution of the accretion disk around black hole is computed. The classical description of the $\alpha$-viscosity is adopted so the evolution is driven by the instability operating in the innermost radiation-pressure dominated part of the accretion disk. We assume that the optically thick disk always extends down to the marginally stable orbit so it is never evacuated completely. We include the effect of the advection, coronal dissipation and vertical outflow. We show that the presence of the corona and/or the outflow reduce the amplitude of the outburst. If only about half of the energy is dissipated in the disk (with the other half dissipated in the corona and carried away by the outflow) the outburst amplitude and duration are consistent with observations of the microquasar GRS 1915+105. Viscous evolution explains in a natural way the lack of direct transitions from the state C to the state B in color-color diagram of this source. Further reduction of the fraction of energy dissipated in the optically thick disk switches off the outbursts which may explain why they are not seen in all high accretion rate sources being in the Very High State.

We observed 17 optically-selected, radio-quiet high-redshift quasars with the Chandra Observatory ACIS, and detected 16 of them. The quasars have redshift between 3.70 and 6.28 and include the highest redshift quasars known. When compared to low-redshift quasars observed with ROSAT, these high redshift quasars are significantly more X-ray quiet. We also find that the X-ray spectral index of the high redshift objects is flatter than the average at lower redshift. These trends confirm the predictions of models where the accretion flow is described by a cold, optically-thick accretion disk surrounded by a hot, optically thin corona, provided the viscosity parameter alpha >= 0.02. The high redshift quasars have supermassive black holes with masses ~10^{10} M_{sun}, and are accreting material at ~0.1 the Eddington limit. We detect 10 X-ray photons from the z=6.28 quasar SDS 1030+0524, which may have a Gunn-Peterson trough and be near the redshift of reionization of the intergalactic medium. The X-ray data place an upper limit on the optical depth of the intergalactic medium tau(IGM) < 10^6, compared to the lower limit from the spectrum of Lyalpha and Lybeta, which implies tau(IGM) > 20.

We propose a theoretical model to explain the spectrum of the quasar PG1211+143 emitted in the Optical/X-ray bands. In particular, we suggest that the inner accretion disk may develop a warm, optically thick skin, which produces a profound emission feature observed in the soft X-ray band. This is well modelled with the Comptonized black body emission. The same warm, mildly ionized medium may also be responsible for the hard X-ray reflection and the presence of the iron $K_{\alpha}$ line. However, in our model it still remains an open question, whether the seed photons for Comptonization come from the cold accretion disc or from the hotter plasma. High resolution spectroscopy available through the Chandra and XMM data may provide now an independent test of the physical conditions in the Comptonizing and reflecting warm skin.

The central engine of the the Narrow Line Seyfert 1 galaxies is being probed. We use the ASCA and RXTE data to model the X-ray primary continuum as well as the reflected component and iron $K\alpha$ line. Since these are strongly coupled, we obtain independent measurements of the disc ionization level and the orientation dependent reflection amplitude. Using the available Optical/UV data we also estimate the black hole masses and the $L/L_{Edd}$ ratios, which are probably related to Boroson & Green eigenvector~1.

We present the study of the emission properties of the quasar PG1211+143, which belongs to the class of Narrow Line Seyfert 1 galaxies. On the basis of observational data analyzed by us and collected from the literature, we study the temporal and spectral variability of the source in the optical/UV/X-ray bands and we propose a model that explains the spectrum emitted in this broad energy range. In this model, the intrinsic emission originating in the warm skin of the accretion disk is responsible for the spectral component that is dominant in the softest X-ray range. The shape of reflected spectrum as well as Fe K line detected in hard X-rays require the reflecting medium to be mildly ionized (xi~500). We identify this reflector with the warm skin of the disk and we show that the heating of the skin is consistent with the classical alpha P_{tot} prescription, while alpha P_{gas} option is at least two orders of magnitude too low to provide the required heating. We find that the mass of the central black hole is relatively small (M_BH~10^7- 10^8 Msun, which is consistent with the Broad Line Region mapping results and characteristic for NLS1 class.

Physical mechanism responsible for high viscosity in accretion disks is still under debate. Parameterization of the viscous stress as $\alpha P$ proved to be a successful representation of this mechanism in the outer parts of the disk, explaining the dwarf novae and X-ray novae outbursts as due to ionization instability. We show that this parameterization can be also adopted in the innermost part of the disk where the adoption of the $\alpha$-viscosity law implies the presence of the instability in the radiation pressure dominated region. We study the time evolution of such disks. We show that the time-dependent behavior of GRS 1915+105 can be well reproduced if $\alpha$-viscosity disk model is calculated accurately (with proper numerical coefficients in vertically averaged equations and with advection included), and if the model is supplemented with (i) moderate corona dissipating 50% of energy (ii) jet carrying luminosity-dependent fraction of energy. These necessary modifications in the form of the presence of a corona and a jet are well justified observationally. The model predicts outbursts at luminosity larger than 0.16$\dot M_{Edd}$, as required, correct outburst timescales and amplitudes, including the effect of increasing outburst timescale with mean luminosity. This result strongly suggests that the $\alpha$-viscosity law is a good description of the actual mechanism responsible for angular momentum transfer also in the innermost, radiation pressure dominated part of the disk around a black hole.

The standard accretion discs are known to be thermally and viscously unstable over a certain range of temperatures. In the inner disc regions there may develop radiation pressure driven instability, which is possibly related to the rapid variability detected in AGNs in the UV range. In the outer disc develops the ionization instability, similar to that in the cataclysmic variables, but operating on much longer timescales. Due to this process the spectrum of the accretion disc differs from that of a stationary one.

We present the results of both analytical and numerical calculations of the amplitude of the reflection component in X-ray spectra of galactic black hole systems. We take into account the anisotropy of Compton scattering and the systematic relativistic bulk motion of the hot plasma. In the case of single scattering approximation the reflection from the disc surface is significantly enhanced due to the anisotropy of Compton scattering. On the other hand the calculations of multiple scattering obtained using the Monte Carlo method show that the anisotropy effect is much weaker in that case. Therefore, the enhanced backscattered flux may affect the observed spectra only if the disc surface is highly ionized, which reduces the absorption in the energy band corresponding to the first Compton scattering.

We present the results of spectral analysis of the ASCA data for the Narrow-Line Seyfert 1 galaxy (NLS1) Ton S180 and simultaneous ASCA and RXTE data modelling for the NLS1 Ark 564. We model both the primary and reflected continuum as well as the iron K alpha line, the energy of which depends on the ionization state of the reprocessor. We show that the reprocessing matter is mildly ionized, and we find the soft to hard luminosity ratio to be about 2.5. The accretion rate approximately corresponds to the Eddington limit value.

We propose a model of the accretion flow onto a black hole consisting of the accretion disk with an accreting two-temperature corona. The model is based on assumptions about the radiative and conductive energy exchange between the two phases and the pressure equilibrium. The complete model is determined by the mass, the accretion rate, and the viscosity parameter. We present the radial dependencies of parameters of such a two-phase flow, with advection in the corona and the disk/corona mass exchange due to evaporation/condensation included, and we determine the transition radius from a two-phase disk/corona accretion to a single-phase optically thin flow (ADAF) in the innermost part of the disk as a function of accretion rate. We identify the NLS1 galaxies with objects accreting at a rate close to the Eddington accretion rate. The strong variability of these objects may be related to the limit cycle behaviour expected in this luminosity range, as the disk, unstable due to the dominance by the radiation pressure, oscillates between the two stable branches: the advection-dominated optically thick branch and the evaporation branch.

We analyse the radiation pressure instability of the standard, optically thisc disc and propose it as an explanation of the obseved variability of the superluminal microquasar GRS 1915+105. The quantitave agreement with the observed outburst amplitudes is obtained under the assumption that substantial part of the energy is carried away by the jet.

We study in detail the properties of the accreting corona model for active galactic nuclei and galactic black holes. In this model the fraction of the energy liberated in the corona at a given radius is calculated from the global parameters of the model (mass of the central object, accretion rate and viscosity parameter) and it appears to be a strong function of the radius. The model predicts the relative decrease of the coronal hard X-ray emission with an increase of the accretion rate. The presented description of disc/corona interaction forms a basis for further studies of disc disruption mechanism.

We consider the role of advection in the two-temperature accreting corona with an underlying optically thick disc. The properties of coronal solutions depend significantly on the description of advection. Local parameterization of advection by a constant coeficient $\delta$ replacing the radial derivatives lead to complex topology of solutions, similar to some extent to other advection-dominated accretion flow solutions. One, radiatively cooled branch exists for low accertion rates. For higher accretion rates two solutions exist in a broad range of radii: one is radiatively cooled and the other one is advection-dominated. With further increase of accretion rate the radial extensions of the two solutions shrink and no solutions are found above certain critical value. However, these trends change if the local parameterization of advection is replaced by proper radial derivatives computed iteratively from the model. Only one, radiatively cooled solution remains, and it exists even for high accretion rates. The advection-dominated branch disappears during the iteration process which means that a self-consistently described advection-dominated flow cannot exist in the presence of an underlying cold disc.