Interferometric visibilities, reflecting the complex correlations between signals recorded at antennas in an interferometric array, carry information about the angular structure of a distant source. While unknown antenna gains in both amplitude and phase can prevent direct interpretation of these measurements, certain combinations of visibilities called closure phases and closure amplitudes are independent of antenna gains and provide a convenient set of robust observables. However, these closure quantities have subtle noise properties and are generally both linearly and statistically dependent. These complications have obstructed the proper use of closure quantities in interferometric analysis, and they have obscured the relationship between analysis with closure quantities and other analysis techniques such as self calibration. We review the statistics of closure quantities, noting common pitfalls that arise when approaching low signal-to-noise due to the nonlinear propagation of statistical errors. We then develop a strategy for isolating and fitting to the independent degrees of freedom captured by the closure quantities through explicit construction of linearly independent sets of quantities along with their noise covariance in the Gaussian limit, valid for moderate signal-to-noise, and we demonstrate that model fits have biased posteriors when this covariance is ignored. Finally, we introduce a unified procedure for fitting to both closure information and partially calibrated visibilities, and we demonstrate both analytically and numerically the direct equivalence of inference based on closure quantities to that based on self calibration of complex visibilities with unconstrained antenna gains.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

The Event Horizon Telescope image of the supermassive black hole in the galaxy M87 is dominated by a bright, unresolved ring. General relativity predicts that embedded within this image lies a thin "photon ring," which is composed of an infinite sequence of self-similar subrings that are indexed by the number of photon orbits around the black hole. The subrings approach the edge of the black hole "shadow," becoming exponentially narrower but weaker with increasing orbit number, with seemingly negligible contributions from high order subrings. Here, we show that these subrings produce strong and universal signatures on long interferometric baselines. These signatures offer the possibility of precise measurements of black hole mass and spin, as well as tests of general relativity, using only a sparse interferometric array.

Context. The Event Horizon Telescope (EHT) collaboration recently obtained first images of the surroundings of the supermassive compact object M87* at the center of the galaxy M87. Aims. We want to develop a simple analytic disk model for the accretion flow of M87*. Compared to general-relativistic magnetohydrodynamic (GRMHD) models, it has the advantage of being independent of the turbulent character of the flow, and controlled by only few easy-to-interpret, physically meaningful parameters. We want to use this model to predict the image of M87* assuming that it is either a Kerr black hole, or an alternative compact object. Methods. We compute the synchrotron emission from the disk model and propagate the resulting light rays to the far-away observer by means of relativistic ray tracing. Such computations are performed assuming different spacetimes (Kerr, Minkowski, non-rotating ultracompact star, rotating boson star or Lamy spinning wormhole). We perform numerical fits of these models to the EHT data. Results. We discuss the highly-lensed features of Kerr images and show that they are intrinsically linked to the accretion-flow properties, and not only to gravitation. This fact is illustrated by the notion of secondary ring that we introduce. Our model of spinning Kerr black hole predicts mass and orientation consistent with the EHT interpretation. The non-Kerr images result in similar quality of the numerical fits and may appear very similar to Kerr images, once blurred to the EHT resolution. This implies that a strong test of the Kerr spacetime may be out of reach with the current data. We notice that future developments of the EHT could alter this situation. Conclusions. Our results show the importance of studying alternatives to the Kerr spacetime in order to be able to test the Kerr paradigm unambiguously.

High-resolution imaging of supermassive black holes is now possible, with new applications to testing general relativity and horizon-scale accretion and relativistic jet formation processes. Over the coming decade, the EHT will propose to add new strategically placed VLBI elements operating at 1.3mm and 0.87mm wavelength. In parallel, development of next-generation backend instrumentation, coupled with high throughput correlation architectures, will boost sensitivity, allowing the new stations to be of modest collecting area while still improving imaging fidelity and angular resolution. The goal of these efforts is to move from imaging static horizon scale structure to dynamic reconstructions that capture the processes of accretion and jet launching in near real time.

In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerful confirmation of Einstein's theory of gravity, or general relativity, was made possible by assembling a global network of radio telescopes operating at millimetre wavelengths that for the first time included the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA as an anchor station has enabled a giant leap forward by increasing the sensitivity limits of the EHT by an order of magnitude, effectively turning it into an imaging array. The published image demonstrates that it is now possible to directly study the event horizon shadows of SMBHs via electromagnetic radiation, thereby transforming this elusive frontier from a mathematical concept into an astrophysical reality. The expansion of the array over the next few years will include new stations on different continents - and eventually satellites in space. This will provide progressively sharper and higher-fidelity images of SMBH candidates, and potentially even movies of the hot plasma orbiting around SMBHs. These improvements will shed light on the processes of black hole accretion and jet formation on event-horizon scales, thereby enabling more precise tests of general relativity in the truly strong field regime.

We report on a new class of solutions of black hole accretion disks that we have found through three-dimensional, global, radiative magnetohydrodynamic simulations in general relativity. It combines features of the canonical thin, slim and thick disk models but differs in crucial respects from each of them. We expect these new solutions to provide a more realistic description of black hole disks than the slim disk model. We are presenting a disk solution for a non-spinning black hole at a sub-Eddington mass accretion rate, $\dot M=0.6\,\dot M_{\rm Edd}$. By the density scale-height measure the disk appears to be thin, having a high density core near the equatorial plane of height $h_{\rho} \sim 0.1 \,r$, but most of the inflow occurs through a highly advective, turbulent, optically thick, Keplerian region that sandwiches the core and has a substantial geometrical thickness comparable to the radius, $H \sim r$. The accreting fluid is supported above the midplane in large part by the magnetic field, with the gas and radiation to magnetic pressure ratio $\beta \sim 1$, this makes the disk thermally stable, even though the radiation pressure strongly dominates over gas pressure. A significant part of the radiation emerging from the disk is captured by the black hole, so the disk is less luminous than a thin disk would be at the same accretion rate.

Operating 1.5 million km from Earth at the Sun-Earth L2 Lagrange point, the Origins Space Telescope equipped with a slightly modified version of its HERO heterodyne instrument could function as a uniquely valuable node in a VLBI network. The unprecedented angular resolution resulting from the combination of Origins with existing ground-based millimeter/submillimeter telescope arrays would increase the number of spatially resolvable black holes by a factor of a million, permit the study of these black holes across all of cosmic history, and enable new tests of general relativity by unveiling the photon ring substructure in the nearest black holes.

The Event Horizon Telescope (EHT) recently produced the first horizon-scale image of a supermassive black hole. Expanding the array to include a 3-meter space telescope operating at >200 GHz enables mass measurements of many black holes, movies of black hole accretion flows, and new tests of general relativity that are impossible from the ground.

We present the design and implementation of an automated data calibration and reduction pipeline for very-long-baseline interferometric (VLBI) observations taken at millimeter wavelengths. These short radio-wavelengths provide the best imaging resolution available from ground-based VLBI networks such as the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA), but require specialized processing due to the strong effects from atmospheric opacity and turbulence as well as the heterogeneous nature of existing global arrays. The pipeline builds upon a calibration suite (HOPS) originally designed for precision geodetic VLBI. To support the reduction of data for astronomical observations, we have developed an additional framework for global phase and amplitude calibration which provides output in a standard data format for astronomical imaging and analysis. We test the pipeline on observations taken at 3.5 mm (86 GHz) by the GMVA joined by the phased Atacama Large Millimeter/submillimeter Array in April 2017, and demonstrate the benefits from the specialized processing of high frequency VLBI data with respect to classical analysis techniques.

We investigate three-dimensional relativistic trajectories of test particles in the spacetime of a slowly rotating compact star, under the combined influence of gravity and a strong, near-Eddington radiation field. While in the static case a spherically symmetric shell of matter suspended above the stellar surface can be formed at the location of radial equilibrium of effective forces, the same is not true for a rotating star. In the latter case the symmetry is broken by the interplay between motion in the non-static spacetime and the influence of strong radiation drag forces, pushing particles towards the equatorial plane. As a result an expanding spherical shell of matter ejected from the neutron star surface collapses on a short timescale into a single stable equatorial ring supported by radiation. These findings have implications for the geometry of optically thin outflows during luminous neutron star bursts.

Neutron stars with near-Eddington observable luminosities were shown to harbor levitating atmospheres, suspended above their surface. We report a new method to simultaneously measure the mass and radius of a neutron star based on oscillations of such atmospheres. In this paper, we present an analytic derivation of a family of relativistic, oscillatory, spherically symmetric eigenmodes of the optically and geometrically thin levitating atmospheres, including the damping effects induced by the radiation drag. We discover characteristic maxima in the frequencies of the damped oscillations and show that using the frequency maxima, one can estimate mass and radius of the neutron star, given the observed frequency and the corresponding luminosity of the star during the X-ray burst. Thus, our model provides a new way to probe the stellar parameters. We also show that the ratio of any two undamped eigenfrequencies depends only on the adiabatic index of the atmosphere, while for the damped eigenfrequencies, this ratio varies with the luminosity. The damping coefficient is independent of the mode number of the oscillations. Signatures of these atmospheres' dynamics will be reflected in the source's X-ray light curves.

(Abridged) The CASA software suite, can now reduce very long baseline interferometry (VLBI) data with the recent addition of a fringe fitter. Here, we present the Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD), which is an open-source VLBI calibration and imaging pipeline built on top of the CASA framework. The pipeline is capable of reducing data from different VLBI arrays. It can be run non-interactively after only a few non-default input parameters are set and delivers high-quality calibrated data. CPU scalability based on a message-passing interface (MPI) implementation ensures that large bandwidth data from future arrays can be processed within reasonable computing times. Phase calibration is done with a Schwab-Cotton fringe fit algorithm. For the calibration of residual atmospheric effects, optimal solution intervals are determined based on the signal-to-noise ratio (S/N) of the data for each scan. Different solution intervals can be set for different antennas in the same scan to increase the number of detections in the low S/N regime. These novel techniques allow rPICARD to calibrate data from different arrays, including high-frequency and low-sensitivity arrays. The amplitude calibration is based on standard telescope metadata, and a robust algorithm can solve for atmospheric opacity attenuation in the high-frequency regime. Standard CASA tasks are used for CLEAN imaging and self-calibration. In this work we demonstrate the capabilities of rPICARD by calibrating and imaging 7 mm VLBA data of the central radio source in the M87 galaxy. The reconstructed jet image reveals a complex collimation profile and edge-brightened structure. A potential counter-jet is detected that has 10 % of the brightness of the approaching jet. This constrains jet speeds close to the radio core to about half the speed of light for small inclination angles.

We present a detailed description of the special procedures for calibration and quality assurance of Atacama Large Millimeter/submillimeter Array (ALMA) observations in Very Long Baseline Interferometry (VLBI) mode. These procedures are required to turn the phased ALMA array into a fully calibrated VLBI station. As an illustration of these methodologies, we present full-polarization observations carried out with ALMA as a phased array at 3mm (Band 3) and 1.3mm (Band 6) as part of Cycle-4. These are the first VLBI science observations conducted with ALMA and were obtained during a 2017 VLBI campaign in concert with other telescopes worldwide as part of the Global mm-VLBI Array (GMVA, April 1-3) and the Event Horizon Telescope (EHT, April 5-11) in ALMA Bands 3 and 6, respectively.

Context. The properties of the accretion flow surrounding the supermassive central black hole of the Galaxy, Sgr A*, will be scrutinized by the new-generation instrument GRAVITY and the Event Horizon Telescope (EHT). Developing fast, robust, and simple models of such flows is thus important and very timely. Aims. We want to model the quiescent emission of Sgr A* from radio to mid-infrared by considering a magnetized compact torus and an extended jet. Results. We find perfect spectral fit both for face-on and edge-on views. These best fits give parameters values very close to that found by the most recent numerical simulations, which are much more complex than our model. The intrinsic radio size of Sgr A* is found to be in reasonable agreement with the centimetric observed constraints. Our best-fit infrared spectral index is in perfect agreement with the latest constraints. Our emission region at 1.3 mm, although larger than the Doeleman et al. (2008) Gaussian best-fit, does contain bright features at the 40 microarcsec scale. EHT-reconstructed images show that torus/jet-specific features persist after the reconstruction procedure, and that these features are sensitive to inclination. Conclusions. The main interest of our model is to give a simple and fast model of the quiescent state of Sgr A*, which gives extremely similar results as compared to state-of-the-art numerical simulations. Our model is easy to use and we publish all the material necessary to reproduce our spectra and images, so that anyone interested can use our results rather straightforwardly. We hope that such a public tool can be useful in the context of the recent and near-future GRAVITY and EHT results. Our model can in particular be easily used to test alternative compact objects models, or alternative gravity theories. The main limitation of our model is that we do not yet treat the X-ray emission.

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI.

If primordial black holes with masses of $10^{25}\,\mbox{g}\gtrsim m \gtrsim 10^{17}\,\mbox{g}$ constitute a non-negligible fraction of the galactic dark-matter haloes, their existence should have observable consequences: they necessarily collide with galactic neutron stars, nest in their centers and accrete the dense matter, eventually converting them to neutron-star mass black holes while releasing the neutron-star magnetic field energy. Such processes may explain the fast radio bursts phenomenology, in particular their millisecond durations, large luminosities ${\sim}10^{43}$ erg/s, high rate of occurrence $\gtrsim 1000/\mbox{day}$, as well as high brightness temperatures, polarized emission and Faraday rotation. Longer than the dynamical timescale of the Bondi-like accretion for light primordial black holes allows for the repeating fast radio bursts. This explanation follows naturally from (assumed) existence of the dark matter primordial black holes and requires no additional unusual phenomena, in particular no unacceptably large magnetic fields of neutron stars. In our model, the observed rate of fast radio bursts throughout the Universe follows from the presently known number of neutron stars in the Galaxy.

We present a numerical method which evolves a two-temperature, magnetized, radiative, accretion flow around a black hole, within the framework of general relativistic radiation magnetohydrodynamics. As implemented in the code KORAL, the gas consists of two sub-components -- ions and electrons -- which share the same dynamics but experience independent, relativistically consistent, thermodynamical evolution. The electrons and ions are heated independently according to a standard prescription from the literature for magnetohydrodynamical turbulent dissipation. Energy exchange between the particle species via Coulomb collisions is included. In addition, electrons gain and lose energy and momentum by absorbing and emitting synchrotron and bremsstrahlung radiation, and through Compton scattering. All evolution equations are handled within a fully covariant framework in the relativistic fixed-metric spacetime of the black hole. Numerical results are presented for five models of low luminosity black hole accretion. In the case of a model with a mass accretion rate $\dot{M}\sim10^{-8} \dot M_{\rm Edd}$, we find that radiation has a negligible effect on either the dynamics or the thermodynamics of the accreting gas. In contrast, a model with a larger $\dot{M}\sim 4\times 10^{-4} \dot M_{\rm Edd}$ behaves very differently. The accreting gas is much cooler and the flow is geometrically less thick, though it is not quite a thin accretion disk.

We present an extension to the general relativistic radiation magnetohydrodynamic code HARMRAD to account for emission and absorption by thermal cyclo-synchrotron, double Compton, bremsstrahlung, low-temperature OPAL opacities as well as Thomson and Compton scattering. We approximate the radiation field as a Bose-Einstein distribution and evolve it using the radiation number-energy-momentum conservation equations in order to track photon hardening. We perform various simulations to study how these extensions affect the radiative properties of magnetically-arrested disks accreting at Eddington to super-Eddington rates. We find that double Compton dominates bremsstrahlung in the disk within a radius of $r\sim 15r_g$ (gravitational radii) at a hundred times the Eddington accretion rate, and within smaller radii at lower accretion rates. Double Compton and cyclo-synchrotron regulate radiation and gas temperatures in the corona, while cyclo-synchrotron regulates temperatures in the jet. Interestingly, as the accretion rate drops to Eddington, an optically thin corona develops whose gas temperature of $T\sim 10^9$K is $\sim 100$ times higher than the disk's black body temperature. Our results show the importance of double Compton and synchrotron in super-Eddington disks, magnetized coronae, and jets.

We construct models of static, spherically symmetric shells supported by the radiation flux of a luminous neutron star in the Schwarzschild metric. The atmospheres are disconnected from the star and levitate above its surface. Gas pressure and density inversion appear in the inner region of these atmospheres, which is a purely relativistic phenomenon. We account for the scattering opacity dependence on temperature and utilize the relativistic M1 closure scheme for the radiation tensor, hence allowing for a fully GR-consistent treatment of the photon flux and radiation tensor anisotropy. In this way we are able to address atmospheres of both large and moderate/low optical depths with the same set of equations. We discuss properties of the levitating atmospheres and find that they may indeed be optically thick, with the distance between star surface and the photosphere expanding as luminosity increases. These results may be relevant for the photosphereric radius expansion X-ray bursts.

Polish doughnuts (PDs) are geometrically thick disks that rotate with super-Keplerian velocities in their innermost parts, and whose long and narrow funnels along rotation axes collimate the emerging radiation into beams. In this paper we construct an extremal family of PDs that maximize both geometrical thickness and radiative efficiency. We then derive upper limits for these quantities and subsequently for the related ability to collimate radiation. PDs with such extreme properties may explain the observed properties of the ultraluminous X-ray sources without the need for the black hole masses to exceed ~ 10 solar masses. However, we show that strong advective cooling, which is expected to be one of the dominant cooling mechanisms in accretion flows with super-Eddington accretion rates, tends to reduce the geometrical thickness and luminosity of PDs substantially. We also show that the beamed radiation emerging from the PD funnels corresponds to "isotropic" luminosities that linearly scale with the mass accretion rate, and do not obey the familiar and well-known logarithmic relation.

In general relativity static gaseous atmospheres may be in hydrostatic balance in the absence of a supporting stellar surface, provided that the luminosity is close to the Eddington value. We construct analytic models of optically thin, spherically symmetric shells supported by the radiation pressure of a luminous central body in the Schwarzschild metric. Opacity is assumed to be dominated by Thomson scattering. The inner parts of the atmospheres, where the luminosity locally has supercritical values, are characterized by a density and pressure inversion. The atmospheres are convectively and Rayleigh-Taylor stable, and there is no outflow of gas.

We investigate the local stability of strongly magnetized relativistic tori orbiting Kerr black holes, for the case of a purely toroidal magnetic field topology. Our approach encompasses both solving the full dispersion relation numerically and solving a simplified analytic treatment. Both approaches indicate that such tori are subject to an unstable non-axisymmetric magnetorotational mode, regardless of such background properties as the gas-to-magnetic pressure ratio and specific angular momentum distribution. We demonstrate that our modal analysis matches well with what is seen in global, three-dimensional, general relativistic magnetohydrodynamic simulations.

Context: An impulsive increase in luminosity by one half or more of the Eddington value will lead to ejection of all optically thin plasma from Keplerian orbits around the radiating star, if gravity is Newtonian and the Poynting-Robertson drag is neglected. Radiation drag may bring some particles down to the stellar surface. On the other hand, general relativistic calculations show that gravity may be balanced by a sufficiently intense radiation field at a certain distance from the star. Aims: We investigate the motion of test particles around highly luminous stars to determine conditions under which plasma may be ejected from the system. Results: In Einstein's gravity, if the outburst is close to the Eddington luminosity, all test particles orbiting outside an "escape sphere" will be ejected from the system, while all others will be captured from their orbits onto the surface of another sphere, which is well above the stellar surface, and may even be outside the escape sphere, depending on the value of luminosity. Radiation drag will bring all the captured particles to rest on this "Eddington capture sphere," where they will remain suspended in an equilibrium state as long as the local flux of radiation does not change and remains at the effective Eddington value.

We show that results of a simple dynamical gedanken experiment interpreted according to standard Newton's gravitational theory, may reveal that three-dimensional space is curved. The experiment may be used to reconstruct the curved geometry of space, i.e. its non-Euclidean metric. The perihelion of Mercury advance and the light bending calculated from the Poisson equation and the equation of motion in the curved geometry have the correct (observed) values. Independently, we also show that Newtonian gravity theory may be enhanced to incorporate the curvature of three dimensional space by adding an extra equation which links the Ricci scalar with the density of matter. Like in Einstein's general relativity, matter is the source of curvature. In the spherically symmetric (vacuum) case, the metric of space 3gik that follows from this extra equation agrees, to the expected accuracy, with the metric measured by the Newtonian gedanken experiment mentioned above.

We present a toy model of mildly super-Eddington, optically thin accretion onto a compact star in the Schwarzschild metric, which predicts periodic variations of luminosity when matter is supplied to the system at a constant accretion rate. These are related to the periodic appearance and disappearance of the Eddington Capture Sphere. In the model the frequency is found to vary inversely with the luminosity. If the input accretion rate varies (strictly) periodically, the luminosity variation is quasi-periodic, and the quality factor is inversely proportional to the relative amplitude of mass accretion fluctuations, with its largest value approximately Q= 1/(10 |delta Mdot/Mdot|) attained in oscillations at about 1 to 2 kHz frequencies for a 2 solar mass star.

Test particles infalling from infinity onto a compact spherical star with a mildly super-Eddington luminosity at its surface are typically trapped on the "Eddington Capture Sphere" and do not reach the surface of the star. The presence of a sphere on which radiation pressure balances gravity for static particles was first discovered some twenty five years ago. Subsequently, it was shown to be a capture sphere for particles in radial motion, and more recently also for particles in non-radial motion, in which the Poynting-Robertson radiation drag efficiently removes the orbital angular momentum of the particles, reducing it to zero. Here we develop this idea further, showing that "levitation" on the Eddington sphere (above the stellar surface) is a state of stable equilibrium, and discuss its implications for Hoyle-Lyttleton accretion onto a luminous star. When the Eddington sphere is present, the cross-section of a compact star for actual accretion is typically less than the geometrical cross-section (pi Rsquared), direct infall onto the stellar surface only being possible for relativistic particles, with the required minimum particle velocity at infinity typically ~1/2 the speed of light. We further show that particles on typical trajectories in the vicinity of the stellar surface will also be trapped on the Eddington Capture Sphere.

We derive and investigate the dispersion relation for accretion disks with retarded or advanced heating. We follow the alpha-prescription but allow for a time offset (\tau) between heating and pressure perturbations, as well as for a diminished response of heating to pressure variations. We study in detail solutions of the dispersion relation for disks with radiation-pressure fraction 1 - \beta . For \tau <0 (delayed heating) the number and sign of real solutions for the growth rate depend on the values of the time lag and the ratio of heating response to pressure perturbations, \xi . If the delay is larger than a critical value (e.g., if \Omega \tau <-125 for \alpha =0.1, \beta =0 and \xi =1) two real solutions exist, which are both negative. These results imply that retarded heating may stabilize radiation-pressure dominated accretion disks.