We present the results from an X-ray variability study of IRAS 13224-3809. This is probably the best source for X-ray reverberation studies since it is X-ray bright, extremely variable, and it has been extensively observed with XMM-Newton. We used all the archival XMM-Newton data from the three EPIC cameras (to increase the signal-to-noise) and, given the many observations of the source, we were able to compute the time-lags spectra in three different flux levels/periods. We fitted the time-lags and energy spectra, simultaneously, using a new X-ray reverberation code which computes the time dependent reflection spectra of the disc as a response to an X-ray flash from a point source located on the axis of the black-hole (BH) accretion disc (lamp-post geometry). To the best of our knowledge, this is the first time for an AGN that both time-lags and energy spectra are fitted by a model simultaneously in different flux periods. The model fits in the case when the BH is rapidly rotating are significantly better than the model fits in the case of a Schwarzschild BH. This result strongly favours the hypothesis of a rotating central BH in this source. We also detect significant variations in the height of the X-ray corona. The X-ray height appears to increase from 3-5 gravitational radii when the X-ray luminosity is of the order of 1.5-3 percent of the Eddington limit, up to 10 gravitational radii, when the luminosity doubles.

Reprocessing of primary X-ray emission in the accretion disk of black hole X-ray binaries (BHXRBs) produces the reflection spectrum, with the characteristic Fe K$\alpha$ fluorescence line. Strong low frequency quasi-periodic oscillations (QPOs) are observed from BHXRBs, and the dependence of QPO properties (e.g., phase-lag) on the inclination angle suggests that the observed QPO may be associated with a geometrical effect, e.g., the precession of the X-ray source due to frame-dragging near the spining black hole. Here, in the scenario of Lense-Thirring precession of the X-ray source, we use a Monte-Carlo simulation of radiative transfer to study the irradiation/reflection and the resultant spectral properties including the Fe K$\alpha$ line, as a function of precession phase (time). We found that the reflection fraction, i.e., the ratio of the incident flux towards the disk and the direct flux towards to the observer at infinity, is modulated with the precession phase, which depends on the truncation radius (i.e., the spectral state in the truncated disk model) and the inclination angle. The Fe K$\alpha$ line profile also changes as the primary X-ray source precesses, with the line luminosity and the flux-weighted centroid energy varying with the precession phase. The periodically modulated 2-10 keV continuum flux could apparently lag the line luminosity in phase, if the truncation radius is small enough for Doppler effects due to disk orbital motion to significantly affect the observed radiation.

Accretion disks around stellar-mass black holes (BHs) emit radiation peaking in the soft X-rays when the source is in the thermal state. The emerging photons are polarized and, for symmetry reasons, the polarization integrated over the source is expected to be either parallel or perpendicular to the (projected) disk symmetry axis, because of electron scattering in the disk. However, due to General Relativity effects photon polarization vectors will rotate with respect to their original orientation, by an amount depending on both the BH spin and the observer's inclination. Hence, X-ray polarization measurements may provide important information about strong gravity effects around these sources. Along with the spectral and polarization properties of radiation which reaches directly the observer once emitted from the disk, in this paper we also include the contribution of returning radiation, i.e. photons that are bent by the strong BH gravity to return again on the disk, where they scatter until eventually escaping to infinity. A comparison between our results and those obtained in previous works by different authors show an overall good agreement, despite the use of different code architectures. We finally consider the effects of absorption in the disk material by including more realistic albedo profiles for the disk surface. Our findings in this respect show that considering also the ionization state of the disk may deeply modify the behavior of polarization observables.

We propose the development of X-ray interferometry (XRI), to reveal the universe at high energies with ultra-high spatial resolution. With baselines which can be accommodated on a single spacecraft, XRI can reach 100 $\mu$as resolution at 10 \AA (1.2 keV) and 20 $\mu$as at 2 \AA (6 keV), enabling imaging and imaging-spectroscopy of (for example) X-ray coronae of nearby accreting supermassive black holes (SMBH) and the SMBH `shadow'; SMBH accretion flows and outflows; X-ray binary winds and orbits; stellar coronae within ~100 pc and many exoplanets which transit across them. For sufficiently luminous sources XRI will resolve sub-pc scales across the entire observable universe, revealing accreting binary SMBHs and enabling trigonometric measurements of the Hubble constant with X-ray light echoes from quasars or explosive transients. A multi-spacecraft `constellation' interferometer would resolve well below 1 $\mu$as, enabling SMBH event horizons to be resolved in many active galaxies and the detailed study of the effects of strong field gravity on the dynamics and emission from accreting gas close to the black hole.

The size and geometry of the X-ray emitting corona in AGNs are still not well constrained. Dov\v{c}iak & Done (2016) proposed a method based on calculations assuming a point-like lamp-post corona. To perform more self-consistent calculations of energy spectra of extended coronae, we develop monk, a Monte Carlo radiative transfer code dedicated to calculations of Comptonised spectra in the Kerr spacetime. In monk we assume Klein-Nishina scattering cross section and include all general relativistic effects. We find that for a corona located above the disc, the spectrum is not isotropic, but with harder and less luminous spectra towards observers at lower inclinations, owing to anisotropic illumination of the seed photons. This anisotropy also leads to an underestimated size of the corona if we assume the corona to be a point-like, isotropic source located on the black hole rotation axis, demonstrating the necessity of more self-consistent calculations. We also inspect the effect of motion and geometry of the corona on the emergent spectrum. Finally, we discuss the implication of anisotropic corona emission for the reflection spectrum in AGNs as well as black hole X-ray binaries (BHXRBs). We find that by assuming the corona emission to be isotropic, one may underestimate the soft excess in AGNs and the reflection continuum and iron K fluorescent line flux in BHXRBs.

X-ray reverberation in Active Galactic Nuclei, believed to be the result of the reprocessing of coronal photons by the underlying accretion disc, has allowed us to probe the properties of the inner-most regions of the accretion flow and the central black hole. Our current model (KYNREFREV) computes the time-dependent reflection spectra of the disc as a response to a flash of primary power-law radiation from a point source corona located on the axis of the black hole accretion disc (lamp-post geometry). Full relativistic effects are taken into account. The ionization of the disc is set for each radius according to the amount of the incident primary flux and the density of the accretion disc. We detect wavy residuals around the best-fit reverberation model time lags at high frequencies. This result suggests that the simple lamp-post geometry does not fully explain the X-ray source/disc configuration in Active Galactic Nuclei. There has been a noticeable progress into the development of codes for extended coronae (Wilkins+16, Chainakun & Young 2017, Taylor & Reynolds 2018a,b). Indeed, the model from Chainakun & Young (2017), consisting of two axial point sources illuminating an accretion disc that produce the reverberation lags is able to reproduce the observed time-lag versus frequency spectra. The goal of this paper is to observationally justify the need for an extended corona in order to provide (in the near future) with a mathematical formulation of a model for an extended corona in its simplest form.

We briefly review a fruitful approach to compute a variety of radiation signatures of General Relativity (GR) originating from accretion disks in strong gravity. A set of transfer functions can be pre-computed and then employed to accelerate the ray-traycing computations, generate a series of model spectra, and to fit the model predictions to actual data sets in X-rays. We have been developing this method to examine spectra and light curves and to model the expected polarimetric properties; in particular, to analyse the properties of the "corona--disk-line" geometry with GR effects taken into account. New impetus to this activity is emerging in the anticipation of upcoming missions equipped with the sensitivity to polarimetric properties in X-rays: IXPE (Imaging X-Ray Polarimetry Explorer) and eXTP (enhanced X-ray Timing and Polarimetry mission).

Super-Massive Black Holes reside in galactic nuclei, where they exhibit episodic bright flares due to accretion events. Taking into account relativistic effects, namely, the boosting and lensing of X-ray flares, we further examine the possibility to constraint the mass of the SMBH from the predicted profiles of the observed light curves. To this end, we have studied four bright flares from Sagittarius A*, which exhibit an asymmetric shape consistent with a combination of two intrinsically separate peaks that occur with a specific time delay with respect to each other. We have thus proposed (Karssen et al. 2017, Mon. Not. R. Astron. Soc. 472, 4422) that an interplay of relativistic effects could be responsible for the shape of the observed light curves and we tested the reliability of the method.

We present the first results obtained by the application of the KYNREFREV-reverberation model, which is ready for its use in XSPEC. This model computes the time dependent reflection spectra of the disc as a response to a flash of primary power-law radiation from a point source corona located on the axis of the black hole accretion disc (lamp-post geometry). Full relativistic effects are taken into account. The ionisation of the disc is set for each radius according to the amount of the incident primary flux and the density of the accretion disc. We tested the model by fitting model predictions to the observed time-lag spectra of three Narrow-Line Seyfert 1 galaxies (ARK 564, MCG-6-30-15 and 1H 0707-495), assuming either a rapidly or zero spinning black hole (BH). The time-lags strongly suggest a compact X-ray source, located close to the BH, at a height of approx. 4 gravitational radii. This result does not depend either on the BH spin or the disc ionization. There is no significant statistical difference between the quality of the best-fits in the rapidly and zero spinning BH scenarios in Ark 564 and MCG-6-30-15. But there is an indication that the hypothesis of a non-rotating BH in 1H 0707-495 is not consistent with its time-lag spectrum. Finally, the best-fits to the Ark 564 and 1H 0707-495 data are of rather low quality. We detect wavy-residuals around the best-fit reverberation model time-lags at high frequencies. This result suggests that the simple lamp-post geometry does not fully explain the X-ray source/disc configuration in Active Galactic Nuclei.

We develop a Monte-Carlo code to compute the Compton scattered X-ray flux arising from a hot inner flow which undergoes Lense-Thirring precession. The hot flow intercepts seed photons from an outer truncated thin disk. A fraction of the Comptonized photons will illuminate back the disk and the reflected/reprocessed photons will contribute to the observed spectrum. The total spectrum, including disk thermal emission, hot flow Comptonization, and disk reflection, is modelled within the framework of general relativity, taking light-bending and gravitational redshift into account. The simulations are performed in the context of the Lense-Thirring precession model for the low-frequency quasi-periodic oscillations, so the inner flow is assumed to precess, leading to periodic modulation of the emitted radiation. In this work, we concentrate on the energy-dependent X-ray variability of the model and, in particular, on the evolution of the variability during the spectral transition from hard to soft state, which is implemented by the decrease of the truncation radius of the outer disk towards Innermost Stable Circular Orbit (ISCO). In the hard state where the Comptonizing flow is geometrically thick, the Comptonization is weakly variable with the fractional variability amplitude of $\leq$10\%; in the soft state where the Comptonizing flow is cooled down and thus becomes geometrically thin, and the fractional variability of the Comptonization is highly variable, increasing with photon energy. The fractional variability of the reflection increases with energy, and the reflection emission for low spin is counterintuitively more variable than the one for high spin.

Ultra-Luminous X-ray sources (ULXs) are accreting black holes for which their X-ray properties have been seen to be different to the case of stellar-mass black hole binaries. For most of the cases their intrinsic energy spectra are well described by a cold accretion disc (thermal) plus a curved high-energy emission components. The mass of the black hole (BH) derived from the thermal disc component is usually in the range of 100-10^5 solar masses, which have led to the idea that this can represent strong evidence of the Intermediate Mass Black Holes (IMBH), proposed to exist by theoretical studies but with no firm detection (as a class) so far. Recent theoretical and observational developments are leading towards the idea that these sources are instead compact objects accreting at an unusual super-Eddington regime instead. On the other hand, gravitational waves have been seen to be a useful tool for finding (some of these) IMBHs. We give a brief overview about the recent advent of the discovery of gravitational waves and their relationship with these so far elusive IMBHs.

We address a question whether the observed light curves of X-ray flares originating deep in galactic cores can give us independent constraints on the mass of the central supermassive black hole. To this end we study four brightest flares that have been recorded from Sagittarius A*. They all exhibit an asymmetric shape consistent with a combination of two intrinsically separate peaks that occur at a certain time-delay with respect to each other, and are characterized by their mutual flux ratio and the profile of raising/declining parts. Such asymmetric shapes arise naturally in the scenario of a temporary flash from a source orbiting near a super- massive black hole, at radius of only 10-20 gravitational radii. An interplay of relativistic effects is responsible for the modulation of the observed light curves: Doppler boosting, gravitational redshift, light focusing, and light-travel time delays. We find the flare properties to be in agreement with the simulations (our ray-tracing code sim5lib). The inferred mass for each of the flares comes out in agreement with previous estimates based on orbits of stars; the latter have been observed at radii and over time-scales two orders of magnitude larger than those typical for the X-ray flares, so the two methods are genuinely different. We test the reliability of the method by applying it to another object, namely, the Seyfert I galaxy RE J1034+396.

Ultra-Luminous X-ray sources (ULXs) are accreting black holes for which their X-ray properties have been seen to be different to the case of stellar-mass black hole binaries. For most of the cases their intrinsic energy spectra are well described by a cold accretion disc (thermal) plus a curved high-energy emission components. The mass of the black hole (BH) derived from the thermal disc component is usually in the range of 100-1000 solar masses, which have led to the idea that this might represent strong evidence of the Intermediate Mass Black Holes (IMBH), proposed to exist by theoretical studies but with no firm detection (as a class) so far. Recent theoretical and observational developments are leading towards the idea that these sources are instead stellar-mass BHs accreting at an unusual super-Eddington regime. In this paper we briefly describe the model SLIMULX that can be used in XSPEC for the fit of thermal spectra of slim discs around stellar mass black holes in the super-Eddington regime. This model consistently takes all relativistic effects into account. We present the obtained results from the fit of the X-ray spectra from NGC 5408 X-1.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

The continuum-fitting method is one of the two most advanced methods of determining the black hole spin in accreting X-ray binary systems. There are, however, still some unresolved issues with the underlying disk models. One of them manifests as an apparent decrease in spin for increasing source luminosity. Here, we perform a few simple tests to establish whether outflows from the disk close to the inner radius can address this problem. We employ four different parametric models to describe the wind and compare these to the apparent decrease in spin with luminosity measured in the sources LMC~X-3 and GRS~1915+105. Wind models in which parameters do not explicitly depend on the accretion rate cannot reproduce the spin measurements. Models with mass accretion rate dependent outflows, however, have spectra that emulate the observed ones. The assumption of a wind thus effectively removes the artifact of spin decrease. This solution is not unique; the same conclusion can be obtained with a truncated inner disk model. To distinguish among valid models, high resolution X-ray data and a realistic description of the Comptonization in the wind will be needed.

We analyse and report in detail new near-infrared (1.45 - 2.45 microns) observations of the Dusty S-cluster Object (DSO/G2) during its approach to the black hole at the center of the Galaxy that were carried out with ESO VLT/SINFONI between February and September 2014. Before May 2014 we detect spatially compact Br-gamma and Pa-alpha line emission from the DSO at about 40mas east of SgrA*. The velocity of the source, measured from the red-shifted emission, is 2700+-60 km/s. No blue-shifted emission above the noise level is detected at the position of SgrA* or upstream the presumed orbit. After May we find spatially compact Br-gamma blue-shifted line emission from the DSO at about 30mas west of SgrA* at a velocity of -3320+-60 km/s and no indication for significant red-shifted emission. We do not detect any significant extension of velocity gradient across the source. We find a Br-gamma-line full width at half maximum of 50+-10 Angstroem before and 15+-10 Angstroem after the peribothron transit, i.e. no significant line broadening with respect to last year is observed. Br-gamma line maps show that the bulk of the line emission originates from a region of less than 20mas diameter. This is consistent with a very compact source on an elliptical orbit with a peribothron time passage in 2014.39+-0.14. For the moment, the flaring activity of the black hole in the near-infrared regime has not shown any statistically significant increment. Increased accretion activity of SgrA* may still be upcoming. We discuss details of a source model according to which the DSO is rather a young accreting star than a coreless gas and dust cloud.

We present a statistical analysis of polarized near-infrared (NIR) light from Sgr A*, the radio source associated with the supermassive black hole at the center of the Milky Way. The observations have been carried out using the adaptive optics instrument NACO at the VLT UT4 in the infrared $K_\mathrm{s}$-band from 2004 to 2012. Several polarized flux excursions were observed during these years. Linear polarization at 2.2 $\mu m$, its statistics and time variation, can be used constrain the physical conditions of the accretion process onto this supermassive black hole. With an exponent of about 4 for the number density histogram of fluxes above 5~mJy, the distribution of polarized flux density is closely linked to the single state power-law distribution of the total $K_\mathrm{s}$-band flux densities reported earlier. We find typical polarization degrees of the order of 20\%$\pm$10\% and a preferred polarization angle of 13$^o$$\pm$15$^o$. Simulations show the uncertainties under a total flux density of $\sim 2\,{\rm mJy}$ are probably dominated by observational effects. At higher flux densities there are intrinsic variations of polarization degree and angle within rather well constrained ranges. Since the emission is most likely due to optically thin synchrotron radiation, this preferred polarization angle we find is very likely coupled to the intrinsic orientation of the Sgr A* system i.e. a disk or jet/wind scenario associated with the super massive black hole. If they are indeed linked to structural features of the source the data imply a rather stable geometry and accretion process for the Sgr A* system.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

A fast moving infrared excess source (G2) which is widely interpreted as a core-less gas and dust cloud approaches Sagittarius A* (SgrA*) on a presumably elliptical orbit. VLT K_s-band and Keck K'-band data result in clear continuum identifications and proper motions of this about 19m Dusty S-cluster Object (DSO). In 2002-2007 it is confused with the star S63, but free of confusion again since 2007. Its near-infrared (NIR) colors and a comparison to other sources in the field speak in favor of the DSO being an IR excess star with photospheric continuum emission at 2 microns than a core-less gas and dust cloud. We also find very compact L'-band emission ($<$0.1'') contrasted by the reported extended (0.03'' up to about 0.2'' for the tail) Brgamma emission. The presence of a star will change the expected accretion phenomena, since a stellar Roche lobe may retain a fraction of the material during and after the peri-bothron passage.

The super-massive 4 million solar mass black hole (SMBH) SgrA* shows flare emission from the millimeter to the X-ray domain. The nucleus of the Milky Way has properties (stellar cluster, young stars, molecular gas and an accreting SMBH) that resemble those of currently higher luminous Low Luminosity Active Galactic Nuclei. A detailed analysis of the infrared light curves shows that the flares are probably generated in a single-state process forming a power-law distribution of the flux density. Near-infrared polarimetry shows signatures of strong gravity that are statistically significant against randomly polarized red noise. Details of the emission mechanism are discussed in a synchrotron/self-Compton model. SgrA* also allows to study the interaction of the SMBH with the immediate interstellar and gaseous environment of the central stellar cluster. Through infrared imaging of the central few arcseconds it is possible to study both inflow and outflow phenomena linked to the SgrA* black hole. In this context we also discuss the newly found dusty object that approaches SgrA* and present a comparison between recent Keck and VLT K-band data that clearly supports its detection as a about 19m K'-band continuum source.

We study the observational signatures of the relativistic slim disc of 10 M_sun black hole, in a wide range of mass accretion rate, mdot, dimensionless spin parameter, a_ast, and viewing angle, i. In general, the innermost temperature, T_in increases with the increase of i for a fixed value of mdot and a_ast, due to the Doppler effect. However, for i > 50 and mdot > mdot_turn, T_in starts to decrease with the increase of mdot. This is a result of self-obscuration -- the radiation from the innermost hot part of the disc is blocked by the surrounding cooler part. The value of mdot_turn and the corresponding luminosities depend on a_ast and i. Such obscuration effects cause an interesting behavior on the disc luminosity (L_disc) -- T_in plane for high inclinations. In addition to the standard-disc branch which appears below mdot_turn and which obeys L_disc propto T_in^4 -relation, another branch above mdot_turn, which is nearly horizontal, may be observed at luminosities close to the Eddington luminosity. We show that these features are likely observed in a Galactic X-ray source, GRS 1915+105. We support a high spin parameter (a_ast > 0.9) for GRS 1915+105 since otherwise the high value of T_in and small size of the emitting region (r_in < 1r_S) cannot be explained.

Most of the action in Active Galactic Nuclei (AGN) occurs within a few tens of gravitational radii from the supermassive black hole, where matter in the accretion disk may lose up to almost half of its energy with a copious production of X-rays, emitted via Comptonization of the disk photons by hot electrons in a corona and partly reflected by the accretion disk. Thanks to its large effective area and excellent energy resolution, Athena+ contributions in the understanding of the physics of accretion in AGN will be fundamental - and unique - in many respects. It will allow us to map the disk-corona system - which is crucial to understand the mechanism of energy extraction and the relation of the corona with winds and jets - by studying the time lags between reflected and primary photons. These lags have been recently discovered by XMM-Newton, but only Athena+ will have the sensitivity required to fully exploit this technique. Athena+ will also be able e.g. to determine robustly the spin of the black hole in nearby sources (and to extend these measurements beyond the local Universe), to establish the nature of the soft X-ray components, and to map the circumnuclear matter within the AGN inner parsec with unprecedented details.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

Slim-disk models describe accretion flows at high luminosities, while reducing to the standard thin disk form in the low luminosity limit. We have developed a new spectral model, slimbb, within the framework of XSPEC, which describes fully relativistic slim-disk accretion and includes photon ray-tracing that starts from the disk photosphere, rather than the equatorial plane. We demonstrate the features of this model by applying it to RXTE spectra of the persistent black-hole X-ray binary LMC X-3. LMC X-3 has the virtues of exhibiting large intensity variations while maintaining itself in soft spectral states which are well described using accretion-disk models, making it an ideal candidate to test the aptness of slimbb. Our results demonstrate consistency between the low-luminosity (thin-disk) and high luminosity (slim-disk) regimes. We also show that X-ray continuum-fitting in the high accretion rate regime can powerfully test black-hole accretion disk models.

We study the process of spinning up black holes by accretion from slim disks in a wide range of accretion rates. We show that for super-Eddington accretion rates and low values of the viscosity parameter alpha (\lesssim 0.01) the limiting value of the dimensionless spin parameter a_* can reach values higher than a_*=0.9978 inferred by Thorne (1974) in his seminal study. For \mdot=10\mdot_{Edd} and alpha=0.01 spin equilibrium is reached at a_*=0.9994. We show that the equilibrium spin value depends strongly on the assumed value of alpha. We also prove that for high accretion rates the impact of captured radiation on spin evolution is negligible.

We report on a scheme for incorporating vertical radiative energy transport into a fully relativistic, Kerr-metric model of optically thick, advective, transonic alpha disks. Our code couples the radial and vertical equations of the accretion disk. The flux was computed in the diffusion approximation, and convection is included in the mixing-length approximation. We present the detailed structure of this "two-dimensional" slim-disk model for alpha=0.01. We then calculated the emergent spectra integrated over the disk surface. The values of surface density, radial velocity, and the photospheric height for these models differ by 20%-30% from those obtained in the polytropic, height-averaged slim disk model considered previously. However, the emission profiles and the resulting spectra are quite similar for both types of models. The effective optical depth of the slim disk becomes lower than unity for high values of the alpha parameter and for high accretion rates.

We study oscillations of accretion torus with azimuthal magnetic field. For several lowest-order modes we calculate eigenfrequencies and eigenfunctions and calculate corresponding intensity and polarization light curves using advanced ray-tracing methods.

We calculate optically thick but geometrically thin (and slim) accretion disk models and perform a ray-tracing of photons (in the Kerr geometry) to calculate the observed disk spectra. Previously, it was a common practice to ray-trace photons assuming that they are emitted from the Kerr geometry equatorial plane, z = 0. We show that the spectra calculated with this assumption differ from these calculated under the assumption that photons are emitted from the actual surface of the disc, z = H(r). This implies that a knowledge of the location of the thin disks effective photosphere is relevant for calculating the spectra. In this paper we investigate, in terms of a simple toy model, a possible influence of the (unknown, and therefore ad hoc assumed) vertical dissipation profiles on the vertical structure of the disk and thus on the location of the effective photosphere, and on the observed spectra. For disks with moderate and high mass accretion rates (\dot m>0.01\dot m_C) we find that the photosphere location in the inner, radiation pressure dominated, disk region (where most of the radiation comes from) does not depend on the dissipation profile and therefore emerging disk spectra are insensitive to the choice of the dissipation function. For lower accretion rates the photosphere location depends on the assumed vertical dissipation profile down to the disk inner edge, but the dependence is very weak and thus of minor importance. We conclude that the spectra of optically thick accretion disks around black holes should be calculated with the ray-tracing from the effective photosphere and that, fortunately, the choice of a particular vertical dissipation profile does not substantially influence the calculated spectrum.

A long discussion has been devoted to the issue of clustering of the kHz QPO frequency ratios in the neutron star sources. While the distribution of ratios inferred from an occurrence of a single QPO seems to be consistent with a random walk, the distribution based on simultaneous detections of both peaks indicates a preference of ratios of small integers. Based on the public RXTE data we further investigate this issue for the source 4U 1636-53. Quality factors and rms amplitudes of both the QPOs nearly equal to the points where the frequencies are commensurable, and where the twin QPO detections cluster. We discuss a connection of the clustering with the varying properties of the two QPO modes. Assuming approximative relations for the observed correlations of the QPO properties, we attempt to reproduce the frequency and ratio distributions using a simple model of a random-walk evolution along the observed frequency-frequency correlation. We obtain results which are in qualitative agreement with the observed distributions.

A recently published study on long term evolution of the frequencies of the kilohertz quasi-periodic oscillations (QPOs) in the atoll source 4U 1636-53 concluded that there is no preferred frequency ratio in a distribution of twin QPOs that was inferred from the distribution of a single frequency alone. However, we find that the distribution of the ratio of actually observed pairs of kHz QPO frequencies is peaked close to the 3/2 value, and possibly also close to the 5/4 ratio. To resolve the apparent contradiction between the two studies, we examine in detail the frequency distributions of the lower kHz QPO and the upper kHz QPO detected in our data set. We demonstrate that for each of the two kHz QPOs (the lower or the upper), the frequency distribution in all detections of a QPO differs from the distribution of frequency of the same QPO in the subset of observations where both the kHz QPOs are detected. We conclude that detections of individual QPOs alone should not be used for calculation of the distribution of the frequency ratios.

Double peak kHz QPO frequencies in neutron star sources varies in time by a factor of hundreds Hz while in microquasar sources the frequencies are fixed and located at the line \nu_2 = 1.5 \nu_1 in the frequency-frequency plot. The crucial question in the theory of twin HFQPOs is whether or not those observed in neutron-star systems are essentially different from those observed in black holes. In black hole systems the twin HFQPOs are known to be in a 3:2 ratio for each source. At first sight, this seems not to be the case for neutron stars. For each individual neutron star, the upper and lower kHz QPO frequencies, \nu_2 and \nu_1, are linearly correlated, \nu_2=A \nu_1 + B, with the slope A < 1.5, i.e., the frequencies definitely are not in a 1.5 ratio. In this contribution we show that when considered jointly on a frequency-frequency plot, the data for the twin kHz QPO frequencies in several (as opposed to one) neutron stars uniquely pick out a certain preferred frequency ratio that is equal to 1.5 for the six sources examined so far.

We study strong-gravity effects on modulation of radiation emerging from accreting compact objects as a possible mechanism for flux modulation in QPOs. We construct a toy model of an oscillating torus in the slender approximation assuming thermal bremsstrahlung for the intrinsic emissivity of the medium and we compute observed (predicted) radiation signal including contribution of indirect (higher-order) images and caustics in the Schwarzschild spacetime. We show that the simplest oscillation mode in an accretion flow, axisymmetric up-and-down motion at the meridional epicyclic frequency, may be directly observable when it occurs in the inner parts of accretion flow around black holes. Together with the second oscillation mode, an in-and-out motion at the radial epicyclic frequency, it may then be responsible for the high-frequency modulations of the X-ray flux observed at two distinct frequencies (twin HF-QPOs) in micro-quasars.

We show that a luminous torus in the Schwarzschild metric oscillating along its own axis gives rise to a periodically varying flux of radiation, even though the source of radiation is steady and perfectly axisymmetric. This implies that the simplest oscillation mode in an accretion flow, axisymmetric up-and-down motion at the meridional epicyclic frequency, may be directly observable when it occurs in the inner parts of accretion flow around neutron stars and black holes. The high-frequency modulations of the X-ray flux observed in low-mass X-ray binaries at two frequencies (twin kHz QPOs) could then be a signature of strong gravity both because radial and meridional oscillations have different frequencies in non-Newtonian gravity, and because strong gravitational deflection of light rays causes the flux of radiation to be modulated at the higher frequency.

We find evidence that the two high frequency QPOs in Sco X-1 are, more often than not, approximately in the 2:3 frequency ratio familiar from studies of black hole candidates (e.g., XTE J1550-564, Remillard et al. 2002). This implies that the double kHz QPO phenomenon in neutron stars has its origin in properties of strong-field gravity and has little to do with the rotation of a stellar surface or any magnetic field structure anchored in the star.