We report on the detection of very-high-energy (VHE; $E > 100$ GeV) $\gamma$-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multi-wavelength observations with Fermi/LAT, XRT and UVOT on board the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum towards the VHE $\gamma$-ray regime, we deduce a 95% confidence level upper limit on the unknown redshift of KUV 00311-1938 of z < 0.98, and of PKS 1440-389 of z < 0.53. When combined with previous spectroscopy results the redshift of KUV 00311-1938 is constrained to $0.51 \leq z < 0.98$ and for PKS 1440-389 to $0.14 \lessapprox z < 0.53$.

Gamma-ray bursts (GRBs) are brief flashes of gamma rays, considered to be the most energetic explosive phenomena in the Universe. The emission from GRBs comprises a short (typically tens of seconds) and bright prompt emission, followed by a much longer afterglow phase. During the afterglow phase, the shocked outflow -- produced by the interaction between the ejected matter and the circumburst medium -- slows down, and a gradual decrease in brightness is observed. GRBs typically emit most of their energy via gamma-rays with energies in the kiloelectronvolt-to-megaelectronvolt range, but a few photons with energies of tens of gigaelectronvolts have been detected by space-based instruments. However, the origins of such high-energy (above one gigaelectronvolt) photons and the presence of very-high-energy (more than 100 gigaelectronvolts) emission have remained elussive. Here we report observations of very-high-energy emission in the bright GRB 180720B deep in the GRB afterglow -ten hours after the end of the prompt emission phase, when the X-ray flux had already decayed by four orders of magnitude. Two possible explanations exist for the observed radiation: inverse Compton emission and synchrotron emission of ultrarelativistic electrons. Our observations show that the energy fluxes in the X-ray and gamma-ray range and their photon indices remain comparable to each other throughout the afterglow. This discovery places distinct constraints on the GRB environment for both emission mechanisms, with the inverse Compton explanation alleviating the particle energy requirements for the emission observed at late times. The late timing of this detection has consequences for the future observations of GRBs at the highest energies.

Flat-spectrum radio-quasars (FSRQs) are rarely detected at very-high-energies (VHE; E>100 GeV) due to their low-frequency-peaked SEDs. At present, only 6 FSRQs are known to emit VHE photons, representing only 7% of the VHE extragalactic catalog. Following the detection of MeV-GeV gamma-ray flaring activity from the FSRQ PKS 0736+017 (z=0.189) with Fermi, the H.E.S.S. array of Cherenkov telescopes triggered ToO observations on February 18, 2015, with the goal of studying the gamma-ray emission in the VHE band. H.E.S.S. ToO observations were carried out during the nights of February 18, 19, 21, and 24, 2015. Together with Fermi-LAT, the multi-wavelength coverage of the flare includes Swift observations in soft-X-rays and optical/UV, and optical monitoring (photometry and spectro-polarimetry) by the Steward Observatory, the ATOM, the KAIT and the ASAS-SN telescope. VHE emission from PKS 0736+017 was detected with H.E.S.S. during the night of February 19, 2015, only. Fermi data indicate the presence of a gamma-ray flare, peaking at the time of the H.E.S.S. detection, with a flux doubling time-scale of around six hours. The gamma-ray flare was accompanied by at least a 1 mag brightening of the non-thermal optical continuum. No simultaneous observations at longer wavelengths are available for the night of the H.E.S.S. detection. The gamma-ray observations with H.E.S.S. and Fermi are used to put constraints on the location of the gamma-ray emitting region during the flare: it is constrained to be just outside the radius of the broad-line-region with a bulk Lorentz factor $\simeq 20$, or at the level of the radius of the dusty torus with Gamma > 60. PKS 0736+017 is the seventh FSRQ known to emit VHE photons and, at z=0.189, is the nearest so far. The location of the gamma-ray emitting region during the flare can be tightly constrained thanks to opacity, variability, and collimation arguments.

The flat spectrum radio quasar 3C 279 is known to exhibit pronounced variability in the high-energy ($100\,$MeV$<E<100\,$GeV) $\gamma$-ray band, which is continuously monitored with Fermi-LAT. During two periods of high activity in April 2014 and June 2015 Target-of-Opportunity observations were undertaken with H.E.S.S. in the very-high-energy (VHE, $E>100\,$GeV) $\gamma$-ray domain. While the observation in 2014 provides an upper limit, the observation in 2015 results in a signal with $8.7\,\sigma$ significance above an energy threshold of $66\,$GeV. No VHE variability has been detected during the 2015 observations. The VHE photon spectrum is soft and described by a power-law index of $4.2\pm 0.3$. The H.E.S.S. data along with a detailed and contemporaneous multiwavelength data set provide constraints on the physical parameters of the emission region. The minimum distance of the emission region from the central black hole is estimated using two plausible geometries of the broad-line region and three potential intrinsic spectra. The emission region is confidently placed at $r\gtrsim 1.7\times10^{17}\,$cm from the black hole, i.e., beyond the assumed distance of the broad-line region. Time-dependent leptonic and lepto-hadronic one-zone models are used to describe the evolution of the 2015 flare. Neither model can fully reproduce the observations, despite testing various parameter sets. Furthermore, the H.E.S.S. data are used to derive constraints on Lorentz invariance violation given the large redshift of 3C 279.

PKS 1830-211 is a known macrolensed quasar located at a redshift of z=2.5. Its high-energy gamma-ray emission has been detected with the Fermi-LAT instrument and evidence for lensing was obtained by several authors from its high-energy data. Observations of PKS 1830-211 were taken with the H.E.S.S. array of Imaging Atmospheric Cherenkov Telescopes in August 2014, following a flare alert by the Fermi- LAT collaboration. The H.E.S.S observations were aimed at detecting a gamma-ray flare delayed by 20-27 days from the alert flare, as expected from observations at other wavelengths. More than twelve hours of good quality data were taken with an analysis threshold of $\sim67$ GeV. The significance of a potential signal is computed as a function of the date as well as the average significance over the whole period. Data are compared to simultaneous observations by Fermi-LAT. No photon excess or significant signal is detected. An upper limit on PKS 1830-211 flux above 67 GeV is computed and compared to the extrapolation of the Fermi-LAT flare spectrum.

We study the oscillations of an axisymmetric, viscous, radiative, general relativistic hydrodynamical simulation of a geometrically thin disk around a non-rotating, $6.62\,M_\odot$ black hole. The numerical setup is initialized with a Novikov-Thorne, gas-pressure-dominated accretion disk, with an initial mass accretion rate of $\dot{m} = 0.01\,L_\mathrm{Edd}/c^2$ (where $L_\mathrm{Edd}$ is the Eddington luminosity and $c$ is the speed of light). Viscosity is treated with the $\alpha$-prescription. The simulation was evolved for about $1000$ Keplerian orbital periods at three Schwarzschild radii (ISCO radius). Power density spectra of the radial and vertical fluid velocity components, the total (gas $+$ radiation) midplane pressure, and the vertical component of radiative flux from the photosphere, all reveal strong power at the local breathing oscillation frequency. The first, second and third harmonics of the breathing oscillation are also clearly seen in the data. We quantify the properties of these oscillations by extracting eigenfunctions of the radial and vertical velocity components and total pressure. This confirms that these oscillations are associated with breathing motion.

We present results from two-dimensional, general relativistic, viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks accreting onto stellar mass Schwarzschild black holes. We consider cases on both the gas- and radiation-pressure-dominated branches of the thermal equilibrium curve, with mass accretion rates spanning the range from $\dot{M} = 0.01 L_\mathrm{Edd}/c^2$ to $10 L_\mathrm{Edd}/c^2$. The simulations directly test the stability of this standard disk model on the different branches. We find clear evidence of thermal instability for all radiation-pressure-dominated disks, resulting universally in the vertical collapse of the disks, which in some cases then settle onto the stable, gas-pressure-dominated branch. Although these results are consistent with decades-old theoretical predictions, they appear to be in conflict with available observational data from black hole X-ray binaries. We also find evidence for a radiation-pressure-driven instability that breaks the unstable disks up into alternating rings of high and low surface density on a timescale comparable to the thermal collapse. Since radiation is included self-consistently in the simulations, we are able to calculate lightcurves and power density spectra (PDS). For the most part, we measure radiative efficiencies (ratio of luminosity to mass accretion rate) close to 6%, as expected for a non-rotating black hole. The PDS appear as broken power laws, with a break typically around 100 Hz. There is no evidence of significant excess power at any frequencies, i.e. no quasi-periodic oscillations are observed.

Index of H.E.S.S. conference proceedings to the 35th ICRC, Busan, Korea

Very high-energy gamma-rays (VHE, E>100 GeV) propagating over cosmological distances can interact with the low-energy photons of the extragalactic background light (EBL) and produce electron-positron pairs. The transparency of the universe to VHE gamma-rays is then directly related to the spectral energy distribution (SED) of the EBL. The observation of features in the VHE energy spectra of extragalactic sources allows the EBL to be measured, which otherwise is very difficult to determine. An EBL-model independent measurement of the EBL SED with the H.E.S.S. array of Cherenkov telescopes is presented. It is obtained by extracting the EBL absorption signal from the reanalysis of high-quality spectra of blazars. From H.E.S.S. data alone the EBL signature is detected at a significance of 9.5 sigma, and the intensity of the EBL obtained in different spectral bands is presented together with the associated gamma-ray horizon.

The diffuse very high-energy (VHE, > 100 GeV) gamma-ray emission observed in the central 200 pc of the Milky Way by H.E.S.S. was found to follow the dense matter distribution in the Central Molecular Zone (CMZ) up to a longitudinal distance of about 130 pc to the Galactic Centre (GC), where the flux rapidly decreases. Recent sensitive H.E.S.S. analysis revealed that the cosmic-ray (CR) density profile drops with the distance to the centre, making data compatible with a steady cosmic PeVatron at the GC. In this paper, we extend this analysis to obtain for the first time a detailed characterisation of the correlation with matter and to search for additional features and individual gamma-ray sources in the inner 200 pc. Taking advantage of 250 hours of H.E.S.S. data and improved analysis techniques we perform a detailed morphology study of the diffuse VHE emission observed from the GC ridge and reconstruct its total spectrum. We show that the emission correlated with dense matter covers the full CMZ and that its flux is about half the total diffuse emission flux. We also detect some emission at higher latitude likely produced by hadronic collisions of CRs in less dense regions of the GC interstellar medium. We detect an additional emission component centred on the GC and extending over about 15 pc that is consistent with the existence of a strong CR density gradient and confirms the presence of a CR accelerator at the very centre of our Galaxy. We show that the spectrum of the full ridge diffuse emission is compatible with the one previously derived from the central regions, suggesting that a single population of particles fills the entire CMZ. Finally, we report the discovery of a VHE gamma-ray source near the GC radio arc and argue that it is produced by the pulsar wind nebula candidate G0.13-0.11.

We performed axisymmetric, grid-based, ideal magnetohydrodynamic (MHD) simulations of oscillating cusp-filling tori orbiting a non-rotating neutron star. A pseudo-Newtonian potential was used to construct the constant angular momentum tori in equilibrium. The inner edge of the torus is terminated by a "cusp" in the effective potential. The initial motion of the model tori was perturbed with uniform sub-sonic vertical and diagonal velocity fields. As the configuration evolved in time, we measured the mass accretion rate on the neutron star surface and obtained the power spectrum. The prominent mode of oscillation in the cusp torus is the radial epicyclic mode. It would appear that vertical oscillations are suppressed by the presence of the cusp. From our analysis it follows that the mass accretion rate carries a modulation imprint of the oscillating torus, and hence so does the boundary layer luminosity.

We have modeled hydrodynamical tori in the Klu{\'z}niak-Lee pseudo-Newtonian potential. The tori in equilibrium were perturbed with uniform sub-sonic velocity fields: vertical, radial and diagonal respectively, and allowed to evolve in time. We identify the eigenmodes corresponding to those of slender tori. The results of our simulations are relevant to the investigation of high-frequency quasi-periodic oscillations observed in stellar-mass black hole binaries.

We consider the three currently known pulsing ultraluminous X--ray sources (PULXs). We show that in one of them the observed spinup rate requires super--Eddington accretion rates at the magnetospheric radius, even if magnetar--strength fields are assumed. In the two other systems a normal--strength neutron star field implies super--Eddington accretion at the magnetosphere. Adopting super--Eddington mass transfer as the defining characteristic of ULX systems, we find the parameters required for self--consistent simultaneous fits of the luminosities and spinup rates of the three pulsed systems. These imply near--equality between their magnetospheric radii $R_M$ and the spherization radii $R_{\rm sph}$ where radiation pressure becomes important and drives mass loss from the accretion disc. We interpret this near--equality as a necessary condition for the systems to appear as pulsed, since if it is violated the pulse fraction is small. We show that as a consequence all PULXs must have spinup rates $\dot\nu \gtrsim 10^{-10}\, {\rm s^{-2}}$, an order of magnitude higher than in any other pulsing neutron--star binaries. The fairly tight conditions required for ULXs to show pulsing support our earlier suggestion that many unpulsed ULX systems must actually contain neutron stars rather than black holes.

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.

We performed axisymmetric hydrodynamical simulations of oscillating tori orbiting a non-rotating black hole. The tori in equilibrium were constructed with a constant distribution of angular momentum in a pseudo-Newtonian potential (Klu{\'z}niak-Lee). Motions of the torus were triggered by adding sub-sonic velocity fields: radial, vertical and diagonal to the tori in equilibrium. As the perturbed tori evolved in time, we measured $L_{2}$ norm of density and obtained the power spectrum of $L_{2}$ norm which manifested eigenfrequencies of tori modes. The most prominent modes of oscillation excited in the torus by a quasi-random perturbation are the breathing mode and the radial and vertical epicyclic modes. The radial and the plus modes, as well as the vertical and the breathing modes will have frequencies in an approximate 3:2 ratio if the torus is several Schwarzschild radii away from the innermost stable circular orbit. Results of our simulations may be of interest in the context of high-frequency quasi-periodic oscillations (HF QPOs) observed in stellar-mass black hole binaries, as well as in supermassive black holes.

Aims: We aim for an understanding of the morphological and spectral properties of the supernova remnant RCW~86 and for insights into the production mechanism leading to the RCW~86 very high-energy gamma-ray emission. Methods: We analyzed High Energy Spectroscopic System data that had increased sensitivity compared to the observations presented in the RCW~86 H.E.S.S. discovery publication. Studies of the morphological correlation between the 0.5-1~keV X-ray band, the 2-5~keV X-ray band, radio, and gamma-ray emissions have been performed as well as broadband modeling of the spectral energy distribution with two different emission models. Results:We present the first conclusive evidence that the TeV gamma-ray emission region is shell-like based on our morphological studies. The comparison with 2-5~keV X-ray data reveals a correlation with the 0.4-50~TeV gamma-ray emission.The spectrum of RCW~86 is best described by a power law with an exponential cutoff at $E_{cut}=(3.5\pm 1.2_{stat})$ TeV and a spectral index of $\Gamma$~$1.6\pm 0.2$. A static leptonic one-zone model adequately describes the measured spectral energy distribution of RCW~86, with the resultant total kinetic energy of the electrons above 1 GeV being equivalent to $\sim$0.1\% of the initial kinetic energy of a Type I a supernova explosion. When using a hadronic model, a magnetic field of $B$~100$\mu$G is needed to represent the measured data. Although this is comparable to formerly published estimates, a standard E$^{-2}$ spectrum for the proton distribution cannot describe the gamma-ray data. Instead, a spectral index of $\Gamma_p$~1.7 would be required, which implies that ~$7\times 10^{49}/n_{cm^{-3}}$erg has been transferred into high-energy protons with the effective density $n_{cm^{-3}}=n/ 1$ cm^-3. This is about 10\% of the kinetic energy of a typical Type Ia supernova under the assumption of a density of 1~cm^-3.

Estimations of black hole spin in the three Galactic microquasars GRS 1915+105, GRO J1655-40, and XTE J1550-564 have been carried out based on spectral and timing X-ray measurements and various theoretical concepts. Among others, a non-linear resonance between axisymmetric epicyclic oscillation modes of an accretion disc around a Kerr black hole has been considered as a model for the observed high-frequency quasi-periodic oscillations (HF QPOs). Estimates of spin predicted by this model have been derived based on the geodesic approximation of the accreted fluid motion. Here we assume accretion flow described by the model of a pressure-supported torus and carry out related corrections to the mass-spin estimates. We find that for dimensionless black hole spin a<0.9, the resonant eigenfrequencies are very close to those calculated for the geodesic motion. Their values slightly grow with increasing torus thickness. These findings agree well with results of a previous study carried out in the pseudo-Newtonian approximation. The situation becomes different for a>0.9, in which case the resonant eigenfrequencies rapidly decrease as the torus thickness increases. We conclude that the assumed non-geodesic effects shift the lower limit of the spin, implied for the three microquasars by the epicyclic model and independently measured masses, from a~0.7 to a~0.6. Their consideration furthermore confirms compatibility of the model with the rapid spin of GRS 1915+105 and provides highly testable predictions of the QPO frequencies. Individual sources with a moderate spin (a<0.9) should exhibit a smaller spread of the measured 3:2 QPO frequencies than sources with a near-extreme spin (a~1). This should be further examined using the large amount of high-resolution data expected to become available with the next generation of X-ray instruments, such as the proposed Large Observatory for X-ray Timing (LOFT).

If the ultraluminous source (ULX) M82 X-2 sustains its measured spin-up value of $\dot \nu= 10^{-10}\,{\rm s^{-2}}$, it will become a millisecond pulsar in less than $10^5\,$ years. The observed (isotropic) luminosity of $10^{40}\,$ erg/s also supports the notion that the neutron star will spin up to a millisecond period upon accreting about $0.1\,{\rm M_\odot}$---the reported hard X-ray luminosity of this ULX, together with the spin-up value, implies torques consistent with the accretion disk extending down to the vicinity of the stellar surface, as expected for low values of the stellar dipole magnetic field ($B\lesssim 10^9\,$G). This suggests a new channel of millisecond pulsar formation---in high-mass X-ray binaries (HMXBs)---and may have implications for studies of gravitational waves, and possibly for the formation of low-mass black holes through accretion-induced collapse.

We numerically investigate the thermally unstable accretion disks around black holes. We adopt an evolutionary viscous stress equation to replace the standard alpha-prescription based on the results of two MHD simulations. We find a kind of interesting oscillations on some running models in limit-cycle outburst state. The oscillations arise near the inner boundary and propagate radially outwards. We deem that they are the trapped $p$-mode oscillations excited by sonic-point instability. We directly integrate the local radiation cooling fluxes to construct the mimic bolometric light-curve. We find a series of overtones beside the fundamental harmonic on the power spectra of mimic light-curves. The frequency of the fundamental harmonic is very close to the maximum epicyclic frequency of the disk and the frequency ratio of the fundamental harmonic and overtones is a regular integer series. We suggest that the code for ray-tracing calculation must be time-dependent in virtual observation and point out the robustness of the black hole spin measurement with high frequency QPOs.

Our goal is to perform global simulations of thin accretion discs around compact bodies like neutron stars with dipolar magnetic profile and black holes by exploiting the facilities provided by state-of-the-art grid-based, high resolution shock capturing (HRSC) and finite volume codes. We have used the Godunov-type code PLUTO to simulate a thin disc around a compact object prescribed with a pseudo-Newtonian potential in a purely hydrodynamical (HD) regime, with numerical viscosity as a first step towards achieving our goal as mentioned above.

We present eigenfrequencies and eigenfunctions of trapped acoustic-inertial oscillations of thin accretion disks for a Schwarzschild black hole and a rapidly rotating Newtonian star (a Maclaurin spheroid). The results are derived in the formalism of Nowak and Wagoner (1991) with the assumption that the oscillatory motion is parallel to the midplane of the disk. The first four radial modes for each of five azimuthal modes $m = 0$ through $m = 4$ are presented. The frequencies and wavefunctions of the lowest modes may be accurately approximated by the Airy function.

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.

Kilohertz QPOs can be used as a probe of the inner regions of accretion disks in compact stars and hence also of the properties of the central object. Most models of kHz QPOs involve epicyclic frequencies to explain their origin. We compute the epicyclic frequencies of nearly circular orbits around rotating strange quark stars. The MIT bag model is used to model the equation of state of quark matter and the uniformly rotating stellar configurations are computed in full general relativity. The vertical epicyclic frequency and the related nodal precession rate of inclined orbits are very sensitive to the oblateness of the rotating star. For slowly rotating stellar models of moderate and high mass strange stars, the sense of the nodal precession changes at a certain rotation rate. At lower stellar rotation rates the orbital nodal precession is prograde, as it is in the Kerr metric, while at higher rotation rates the precession is retrograde, as it is for Maclaurin spheroids. Thus, qualitatively, the orbits around rapidly rotating strange quark stars are affected more strongly by the effects of stellar oblateness than by the effects of general relativity. We show that epicyclic and orbital frequencies calculated numerically for small mass strange stars are in very good agreement with analytical formulae for Maclaurin spheroids.

Context. Some microquasars exhibit millisecond quasi-periodic oscillations (QPO) that are likely related to phenomena occuring in the immediate vicinity of the central black hole. Oscillations of accretion tori have been proposed to model these QPOs. Aims. Here, we aim at determining the observable spectral signature of slender accretion tori surrounding Kerr black holes. We analyze the impact of the inclination and spin parameters on the power spectra. Methods. Ray-traced power spectra of slender tori oscillation modes are computed in the Kerr metric. Results. We show that the power spectral densities of oscillating tori are very sensitive to the inclination and spin parameters. This strong dependency of the temporal spectra on inclination and spin may lead to observable constraints of these parameters. Conclusions. This work goes a step further in the analysis of the oscillating torus QPO model. It is part of a long-term study that will ultimately lead to comparison with observed data.

Accretion disks are ubiquitous in the universe and it is generally accepted that magnetic fields play a pivotal role in accretion-disk physics. The spin history of millisecond pulsars, which are usually classified as magnetized neutron stars spun up by an accretion disk, depends sensitively on the magnetic field structure, and yet highly idealized models from the 80s are still being used for calculating the magnetic field components. We present a possible way of improving the currently used models with a semi-analytic approach. The resulting magnetic field profile of both the poloidal and the toroidal component can be very different from the one suggested previously. This might dramatically change our picture of which parts of the disk tend to spin the star up or down.

Optically thin coronae around neutron stars suffering an X-ray burst can be ejected as a result of rapid increase in stellar luminosity. In general relativity (GR), radiation pressure from the central luminous star counteracts gravitational attraction more strongly than in Newtonian physics. However, motion near the neutron star is very effectively impeded by the radiation field. We discuss coronal ejection in a general relativistic calculation of the motion of a test particle in a spherically symmetric radiation field. At every radial distance from the star larger than that of the ISCO, and any initial luminosity of the star, there exists a luminosity change which leads to coronal ejection. The luminosity required to eject from the system the inner parts of the optically thin neutron-star corona is very high in the presence of radiation drag and always close to the Eddington luminosity. Outer parts of the corona, at a distance of ~20 $R_G$ or more, will be ejected by a sub-Eddington outburst. Mildly fluctuating luminosity will lead to dissipation in the plasma and may explain the observed X-ray temperatures of coronae in low mass X-ray binaries (LMXBs). At large radial distances from the star ($3\cdot 10^3 R_G$ or more) the results do not depend on whether or not Poynting-Robertson drag is included in the calculation.

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.

In some galaxy clusters powerful AGN have blown bubbles with cluster scale extent into the ambient medium. The main pressure support of these bubbles is not known to date, but cosmic rays are a viable possibility. For such a scenario copious gamma-ray emission is expected as a tracer of cosmic rays from these systems. Hydra A, the closest galaxy cluster hosting a cluster scale AGN outburst, located at a redshift of 0.0538, is investigated for being a gamma-ray emitter with the High Energy Stereoscopic System (H.E.S.S.) array and the Fermi Large Area Telescope (Fermi-LAT). Data obtained in 20.2 hours of dedicated H.E.S.S. observations and 38 months of Fermi-LAT data, gathered by its usual all-sky scanning mode, have been analyzed to search for a gamma-ray signal. No signal has been found in either data set. Upper limits on the gamma-ray flux are derived and are compared to models. These are the first limits on gamma-ray emission ever presented for galaxy clusters hosting cluster scale AGN outbursts. The non-detection of Hydra A in gamma-rays has important implications on the particle populations and physical conditions inside the bubbles in this system. For the case of bubbles mainly supported by hadronic cosmic rays, the most favorable scenario, that involves full mixing between cosmic rays and embedding medium, can be excluded. However, hadronic cosmic rays still remain a viable pressure support agent to sustain the bubbles against the thermal pressure of the ambient medium. The largest population of highly-energetic electrons which are relevant for inverse-Compton gamma-ray production is found in the youngest inner lobes of Hydra A. The limit on the inverse-Compton gamma-ray flux excludes a magnetic field below half of the equipartition value of 16 muG in the inner lobes.

The two high frequency quasi periodic oscillations (HFQPOs) recently reported in the black hole candidate IGR J17091-3624 by Altamirano and Belloni (2012) are in a 5:2 frequency ratio (164 Hz to 66 Hz). This ratio is strongly suggestive of period doubling and nonlinear resonance analogous to phenomena known in RV Tauri-type pulsating stars (and recently discovered also in oscillations of RR Lyrae-type and of BL Herculis-type variables). An interpretation of the frequency ratio in terms of nonlinear interactions and a comparison with the HFQPOs reported in GRS 1915+105 may imply a mass of about 6 solar masses for the black hole in IGR J17091- 3624.

Several newly discovered very-high-energy (VHE; E > 100 GeV) gamma-ray sources in the Galaxy are thought to be associated with energetic pulsars. Among them, middle-aged (> 1E+4 yr) systems exhibit large centre-filled VHE nebulae, offset from the pulsar position, which result from the complex relationship between the pulsar wind and the surrounding medium, and reflect the past evolution of the pulsar. Imaging Atmospheric Cherenkov Telescopes (IACTs) have been successful in revealing extended emission from these sources in the VHE regime. Together with radio and X-ray observations, this observational window allows one to probe the energetics and magnetic field inside these large-scale nebulae. H.E.S.S., with its large field of view, angular resolution of < 0.1deg and unprecedented sensitivity, has been used to discover a large population of such VHE sources. In this paper, the H.E.S.S. data from the continuation of the Galactic Plane Survey (-80deg < l < 60deg, |b| < 3deg), together with the existing multi-wavelength observations, are used. A new VHE gamma-ray source was discovered at R.A. (J2000) = 13h56m00s, Dec. (J2000) = -64d30m00s with a 2' statistical error in each coordinate, namely HESS J1356-645. The source is extended, with an intrinsic Gaussian width of (0.20 +/- 0.02)deg. Its integrated energy flux between 1 and 10 TeV of 8E-12 erg cm-2 s-1 represents ~ 11% of the Crab Nebula flux in the same energy band. The energy spectrum between 1 and 20 TeV is well described by a power law dN/dE ~ E-Gamma with photon index Gamma = 2.2 +/- 0.2stat +/- 0.2sys. The inspection of archival radio images at three frequencies and the analysis of X-ray data from ROSAT/PSPC and XMM-Newton/MOS reveal the presence of faint non-thermal diffuse emission coincident with HESS J1356-645. HESS J1356-645 is most likely associated with the young and energetic pulsar PSR J1357-6429 (Abridged)

The H.E.S.S. very-high-energy (VHE, E > 0.1 TeV) gamma-ray telescope system has discovered a new source, HESS J1747-248. The measured integral flux is (1.2 +/- 0.3) \times 10^-12 cm-2 s-1 above 440 GeV for a power-law photon spectral index of 2.5 +/- 0.3 stat +/- 0.2 sys. The VHE gamma-ray source is located in the close vicinity of the Galactic globular cluster Terzan 5 and extends beyond the H.E.S.S. point spread function (0.07 degree). The probability of a chance coincidence with Terzan 5 and an unrelated VHE source is quite low (~ 10^-4). With the largest population of identified millisecond pulsars (msPSRs), a very high core stellar density and the brightest GeV range flux as measured by Fermi-LAT, Terzan 5 stands out among Galactic globular clusters. The properties of the VHE source are briefly discussed in the context of potential emission mechanisms, notably in relation to msPSRs. Interpretation of the available data accommodates several possible origins for this VHE gamma-ray source, although none of them offers a satisfying explanation of its peculiar morphology.

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.

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.

The recent discovery of the radio shell-type supernova remnant (SNR), G353.6-0.7, in spatial coincidence with the unidentified TeV source HESS J1731-347 has motivated further observations of the source with the High Energy Stereoscopic System (H.E.S.S.) Cherenkov telescope array to test a possible association of the gamma-ray emission with the SNR. With a total of 59 hours of observation, representing about four times the initial exposure available in the discovery paper of HESS J1731-347, the gamma-ray morphology is investigated and compared with the radio morphology. An estimate of the distance is derived by comparing the interstellar absorption derived from X-rays and the one obtained from 12CO and HI observations. The deeper gamma-ray observation of the source has revealed a large shell-type structure with similar position and extension (r~0.25{\deg}) as the radio SNR, thus confirming their association. By accounting for the H.E.S.S. angular resolution and projection effects within a simple shell model, the radial profile is compatible with a thin, spatially unresolved, rim. Together with RX J1713.7-3946, RX J0852.0-4622 and SN 1006, HESS J1731-347 is now the fourth SNR with a significant shell morphology at TeV energies. The derived lower limit on the distance of the SNR of 3.2 kpc is used together with radio and X-ray data to discuss the possible origin of the gamma-ray emission, either via inverse Compton scattering of electrons or the decay of neutral pions resulting from proton-proton interaction.

Observations of the globular clusters NGC 6388 and M 15 were carried out by the H.E.S.S. array of Cherenkov telescopes for a live time of 27.2 and 15.2 hours respectively. No gamma-ray signal is found at the nominal target position of NGC 6388 and M 15. In the primordial formation scenario, globular clusters are formed in a dark matter halo and dark matter could still be present in the baryon-dominated environment of globular clusters. This opens the possibility of observing a dark matter self-annihilation signal. The dark matter content of the globular clusters NGC 6388 and M 15 is modelled taking into account the astrophysical processes that can be expected to influence the dark matter distribution during the evolution of the globular cluster: the adiabatic contraction of dark matter by baryons, the adiabatic growth of a black hole in the dark matter halo and the kinetic heating of dark matter by stars. 95% confidence level exclusion limits on the dark matter particle velocity-weighted annihilation cross section are derived for these dark matter haloes. In the TeV range, the limits on the velocity-weighted annihilation cross section are derived at the 10-25 cm3 s-1 level and a few 10-24 cm3 s-1 for NGC 6388 and M 15 respectively.

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.

Observations of very high energy gamma-rays from blazars provide information about acceleration mechanisms occurring in their innermost regions. Studies of variability in these objects allow a better understanding of the mechanisms at play. To investigate the spectral and temporal variability of VHE (>100 GeV) gamma-rays of the well-known high-frequency-peaked BL Lac object PKS 2155-304 with the H.E.S.S. imaging atmospheric Cherenkov telescopes over a wide range of flux states. Data collected from 2005 to 2007 are analyzed. Spectra are derived on time scales ranging from 3 years to 4 minutes. Light curve variability is studied through doubling timescales and structure functions, and is compared with red noise process simulations. The source is found to be in a low state from 2005 to 2007, except for a set of exceptional flares which occurred in July 2006. The quiescent state of the source is characterized by an associated mean flux level of 4.32 +/-0.09 x 10^-11 cm^-2 s^-1 above 200 GeV, or approximately 15% of the Crab Nebula, and a power law photon index of 3.53 +/-0.06. During the flares of July 2006, doubling timescales of ~2 min are found. The spectral index variation is examined over two orders of magnitude in flux, yielding different behaviour at low and high fluxes,which is a new phenomenon in VHE gamma-ray emitting blazars. The variability amplitude characterized by the fractional r.m.s. is strongly energy-dependent and is proportional to E^(0.19 +/- 0.01). The light curve r.m.s. correlates with the flux. This is the signature of a multiplicative process which can be accounted for as a red noise with a Fourier index of ~2. This unique data set shows evidence for a low level gamma-ray emission state from PKS 2155-304, which possibly has a different origin than the outbursts. The discovery of the light curve lognormal behaviour might be an indicator ..

AIMS: The properties of the broad-band emission from the high-frequency peaked BL Lac H 2356-309 (z=0.165) are investigated. METHODS: Very High Energy (VHE; E > 100 GeV) observations of H 2356-309 were performed with the High Energy Stereoscopic System (HESS) from 2004 through 2007. Simultaneous optical/UV and X-ray observations were made with the XMM-Newton satellite on June 12/13 and June 14/15, 2005. NRT radio observations were also contemporaneously performed in 2005. ATOM optical monitoring observations were also made in 2007. RESULTS: A strong VHE signal, ~13 sigma total, was detected by HESS after the four years HESS observations (116.8 hrs live time). The integral flux above 240 GeV is I(>240 GeV) = (3.06 +- 0.26 {stat} +- 0.61 {syst}) x 10^{-12} cm^{-2} s^{-1}, corresponding to ~1.6% of the flux observed from the Crab Nebula. A time-averaged energy spectrum is measured from 200 GeV to 2 TeV and is characterized by a power law (photon index of Gamma = 3.06 +- 0.15 {stat} +- 0.10 {syst}). Significant small-amplitude variations in the VHE flux from H 2356-309 are seen on time scales of months and years, but not on shorter time scales. No evidence for any variations in the VHE spectral slope are found within these data. The XMM-Newton X-ray measurements show a historically low X-ray state, characterized by a hard, broken-power-law spectrum on both nights. CONCLUSIONS: The broad-band spectral energy distribution (SED) of the blazar can be adequately fit using a simple one-zone synchrotron self-Compton (SSC) model. In the SSC scenario, higher VHE fluxes could be expected in the future since the observed X-ray flux is at a historically low level.

Thin viscous Keplerian accretion disks are considered asymptotically stable, even though they can show significant dynamic activity on short timescales. In this paper the dynamics of non-axisymmetric hydrodynamical disturbances of disks are investigated analytically building upon the steady state three-dimensional structure and evolution of axisymmetric perturbations explored in previous work. Assuming a polytropic equation of state solutions are found by means of an asymptotic expansion in the small parameter measuring the ratio of the disk thickness to characteristic radius. In-depth analysis shows that every perturbation that disturbs the radial velocity induces significant transient growth in the (acoustic) energy of the evolving disturbance. This effect is most evident in the density and vertical velocity. The transient growth observed is tied to the non-separable nature of the solutions where, in particular, pattern evolution is controlled by a similarity variable composed of the radial coordinate and time. This leads to growing winding perturbations that display successive radial peaks and troughs. We argue that these transient non-axisymmetric structures may precipitate secondary instabilities which, consequently, may be a critical element for a new alternative picture of turbulence arousal in non-magnetized astrophysical disks.

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.

We explore some properties of twin kilohertz quasiperiodic oscillations (QPOs) in a simple toy-model consisting of two oscillation modes coupled by a general nonlinear force. We examine resonant effects by slowly varying the values of the tunable, and nearly commensurable, eigenfrequencies. The behavior of the actual oscillation frequencies and amplitudes during a slow transition through the 3:2 resonance is examined in detail and it is shown that both are significantly affected by the nonlinearities in the governing equations. In particular, the amplitudes of oscillations reflect a resonant exchange of energy between the modes, as a result the initially weaker mode may become dominant after the transition. We note that a qualitatively similar behavior has been recently reported in several neutron star sources by Torok (2008, arXiv:0812.4751), who found that the difference of amplitudes in neutron star twin peak QPOs changes sign as the observed frequency ratio of the QPOs passes through the value 3:2.

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.

We compute the properties of a geometrically thin, steady accretion disk surrounding a central rotating, magnetized star. The magnetosphere is assumed to entrain the disk over a wide range of radii. The model is simplified in that we adopt two (alternate) ad hoc, but plausible, expressions for the azimuthal component of the magnetic field as a function of radial distance. We find a solution for the angular velocity profile tending to corotation close to the central star, and smoothly matching a Keplerian curve at a radius where the viscous stress vanishes. The value of this ''transition'' radius is nearly the same for both of our adopted B-field models. We then solve analytically for the torques on the central star and for the disk luminosity due to gravity and magnetic torques. When expressed in a dimensionless form, the resulting quantities depend on one parameter alone, the ratio of the transition radius to the corotation radius. For rapid rotators, the accretion disk may be powered mostly by spin-down of the central star. These results are independent of the viscosity prescription in the disk. We also solve for the disk structure for the special case of an optically thick alpha disk. Our results are applicable to a range of astrophysical systems including accreting neutron stars, intermediate polar cataclysmic variables, and T Tauri systems.

We study epicyclic oscillations of fluid tori around black holes (in the Paczynski-Wiita potential) and derive exact analytic expressions for their radial and vertical eigenfrequencies nu_r and nu_z to second-order accuracy in the width of the torus. We prove that pressure effects make the eigenfrequencies smaller than those for free particles. However, the particular ratio nu_z/nu_r=3/2, which is important for the theory of high-frequency quasi-periodic oscillations (QPOs), occurs when the fluid tori epicyclic frequencies nu_r and nu_z are about 15% higher than the ones corresponding to free particles. Our results therefore suggest that previous estimates of black hole spins from QPOs have produced values that are too high.

We consider a new version of the twin paradox. The twins move along the same circular free photon path around the Schwarzschild center. In this case, despite their different velocities, all twins have the same non-zero acceleration. On the circular photon path, the symmetry between the twins situations is broken not by acceleration (as it is in the case of the classic twin paradox), but by the existence of an absolute standard of rest (timelike Killing vector). The twin with the higher velocity with respect to the standard of rest is younger on reunion. This closely resembles the case of periodic motions in compact (non-trivial topology) 3-D space recently considered in the context of the twin paradox by Barrow and Levin, except that there accelerations of all twins were equal to zero, and that in the case considered here, the 3-D space has trivial topology.

Kilohertz Quasi-Periodic Oscillations (QPOs) have been detected in many accreting X-ray binaries. It has been suggested that the highest QPO frequencies observed in the modulation of the X-ray flux reflect a non-linear resonance between two modes of accreting disk oscillation. This hypothesis implies certain very general predictions, several of which have been borne out by observations. Some of these follow from properties of non-linear oscillators, while the others are specific to oscillations of fluid in strong gravity. A 3:2 resonant ratio of frequencies can be clearly recognized in the black-hole as well as in the neutron-star QPO data.

Fluids in external gravity may oscillate with frequencies characteristic of the epicyclic motions of test particles. We explicitly demonstrate that global oscillations of a slender, perfect fluid torus around a Kerr black hole admit incompressible vertical and radial epicyclic modes. Our results may be directly relevant to one of the most puzzling astrophysical phenomena -- high (hundreds of hertz) frequency quasiperiodic oscillations (QPOs) detected in X-ray fluxes from several black hole sources. Such QPOs are pairs of stable frequencies in the 3/2 ratio. It seems that they originate a few gravitational radii away from the black hole and thus observations of them have the potential to become an accurate probe of super-strong gravity.

Some of the more promising ideas about the origin of the high frequency variability (kHz QPOs) in the observed X-ray emissions of low-mass X-ray binaries are contrasted with less promising ones.

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.

Recent observations, reported by Warner and Woudt, of Dwarf Nova Oscillations (DNOs) exhibiting frequency drift, period doubling, and 1:2:3 harmonic structure, can be understood as disc oscillations that are excited by perturbations at the spin frequency of the white dwarf or of its equatorial layers. Similar quasi-periodic disc oscillations in black hole low-mass X-ray binary (LMXB) transients in a 2:3 frequency ratio show no evidence of frequency drift and correspond to two separate modes of disc oscillation excited by an internal resonance. Just as no effects of general relativity play a role in white dwarf DNOs, no stellar surface or magnetic field effects need be invoked to explain the black hole QPOs.

We have carried out a systematic study of the properties of the kilo-Hertz quasi-periodic oscillations (QPO) observed in the X-ray emission of the neutron star low-mass X-ray binary 4U1608-52, using archival data obtained with the Rossi X-ray Timing Explorer. We have investigated the quality factor, Q, of the oscillations (defined as the ratio of the frequency of the QPO peak to its full width at half maximum). In order to minimise the effect of long-term frequency drifts, power spectra were computed over the shortest times permitted by the data statistics. We show that the high Q of ~200 reported by Berger et al. (1996) for the lower frequency kilo-Hz QPO in one of their observations is by no means exceptional, as we observe a mean Q value in excess of 150 in 14 out of the 21 observations analysed and Q can remain above 200 for thousands of seconds. The frequency of the QPO varies over the wide range 560--890 Hz and we find a systematic trend for the coherence time of the QPO, estimated as tau=Q /(pi nu), to increase with the frequency, up to a maximum level at ~ 800 Hz, beyond which it appears to decrease, at frequencies where the QPO weakens. There is a more complex relationship between tau and the QPO root mean squared amplitude (RMS), in which positive and negative correlations can be found. A higher-frequency QPO, revealed by correcting for the frequency drift of the 560-890 Hz one, has a much lower Q (~10) which does not follow the same pattern. We discuss these results in the framework of competing QPO models and show that those involving clumps orbiting within or above the accretion disk are ruled out.

We present and discuss a short and simple derivation of orbital epicyclic frequencies for circular geodesic orbits in stationary and axially symmetric spacetimes. Such spacetimes include as special cases analytically known black hole Kerr and Schwarzschild spacetimes, as well as the analytic Hartle-Thorne spacetime and all numerically constructed spacetimes relevant for rotating neutron stars. Our derivation follows directly from energy and angular momentum conservation and it uses the concept of the effective potential. It has never been published, except for a few special cases, but it has already become a part of the common knowledge in the field.

The origin of quasi periodic modulations of flux in the kilohertz range (kHz QPOs), observed in low-mass X-ray binaries, is usually assumed to be physically distinct from that of the ``normal branch oscillations'' (NBOs) in the Z-sources. We show that a low-frequency modulation of the kHz QPOs is a natural consequence of the non-linear relativistic resonance suggested previously to explain the properties of the high-frequency twin peaks. The theoretical results discussed here are reminiscent of the 6 Hz variations of frequency and amplitude of the kHz QPOs reported by Yu, van der Klis and Jonker (2001).

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 have performed numerical simulations of a radially perturbed "accretion" torus around a black hole or neutron star and find that the torus performs radial and vertical motions at the appropriate epicyclic frequencies. We find clear evidence that vertical motions are excited in a non-linear resonance when the applied perturbation is periodic in time. The strongest resonant response occurs when the frequency difference of the two oscillations is equal to one-half the forcing frequency, precisely as recently observed in the accreting pulsar, SAX J1808.4-3658, where the observed kHz QPO peak separation is half the spin frequency of 401 Hz.

Recently, twin-peak QPOs have been observed in a 3:2 ratio for three Galactic black-hole microquasars with frequencies that have been shown to scale as 1/M, as expected for general relativisitic motion near a black hole. It may be possible to extend this result to distinguish between the following two disparate models that have been proposed for the puzzling ultraluminous X-ray sources (ULXs): (1) an intermediate-mass black hole M ~1000 solar mass) radiating very near the Eddington limit and (2) a conventional black hole (M ~ 10 solar mass) accreting at a highly super-Eddington rate with its emission beamed along the rotation axis. We suggest that it may be possible to distinguish between these models by detecting the counterpart of a Galactic twin-peak QPO in a ULX: the expected frequency for the intermediate-mass black hole model is only about 1 Hz, whereas, for the conventional black hole model the expected frequency would be the ~100 Hz value observed for the Galactic microquasars.

The presence of a kHz frequency in LMXBs has been expected from scaling laws, by analogy with the QPO phenomenon in HMXB X-ray pulsars. Interpretation of the two kHz frequencies, observed in accreting neutron stars, in terms of non-linear resonance in strong-field gravity led to the prediction of twin QPOs in black hole systems, in a definite frequency ratio (such as 2/3). The imprint of a subharmonic of the 401 Hz rotation rate in the frequencies of the QPOs detected in the accreting millisecond pulsar is at once a signature of non-linear resonance and of coupling between accretion disk modes and the neutron star spin.

Many Galactic black hole and neutron star sources in low X-ray mass binaries show QPOs (quasi periodic oscillations) in their observed X-ray fluxes, i.e. peaks in the Fourier variability power spectra. Pairs of twin peaks are observed, and in black-hole systems their frequencies (upp), (down) are in rational ratios. For example, in all four microquasars with twin peaks observed, (upp):(down) = 3:2. The rational ratios have been postulated in a model that explained twin peak QPOs as a non-linear resonance between modes of accretion disk oscillations. For those microquasars where the mass of the X-ray source is known, we determine the black-hole spin, following from the observed QPO frequencies within various models of resonance.

In all the microquasars with two hHz QPOs, the ratio of the frequencies is 3:2, supporting our suggestion that a non-linear resonance between two modes of oscillation in the accretion disk plays a role in exciting the observed modulations of the X-ray flux. We discuss the evidence in favor of this interpretation, and we relate the black hole spin to the frequencies expected for various types of resonances that may occur in nearly Keplerian disks in strong gravity. For those microquasars where the mass of the central X-ray source is known, the black hole spin can be deduced from a comparison of the observed and expected frequencies.

Twin quasi-periodic millisecond modulations of the X-ray flux (kHz QPOs) have recently been reported from an accreting 2.5 ms X-ray pulsar. We identify modes of disk oscillations whose frequencies are in agreement with the observed ones when the rotating neutron star is modeled with realistic equations of state. The frequency difference of the twin QPOs, equal to about one half of the neutron-star spin rate, clearly indicates that resonant oscillations of the accretion disk have been observed. Similar non-linear resonances may also be spontaneously excited in the accretion disk. The two QPO frequencies in the pulsar system are close to a 5:7 ratio and this suggests a link with the QPOs in black hole systems, where frequency ratios of 2:3 and 3:5 have been reported.

Millisecond variability detected in low-mass X-ray binaries (LMXBs) probably reflects motions of accreting fluid in the strong-field gravity of neutron stars or neutron stars. Parametric resonance between two oscillators with unequal eigenfrequencies provides a natural explanation for the 2:3 frequency ratio observed in some black-hole systems, and probably also in the brightest neutron stars. The two oscillators likely correspond to the meridional and radial epicyclic motions.

We have explored the ideas that parametric resonance affects nearly geodesic motion around a black hole or a neutron star, and that it may be relevant to the high frequency (twin) quasi-periodic oscillations occurring in some low-mass X-ray binaries. We have assumed the particles or fluid elements of an accretion disc to be subject to an isotropic perturbation of a hypothetical but rather general form. We find that the parametric resonance is indeed excited close to the radius where epicyclic frequencies of radial and meridional oscillations are in a 2:3 ratio. The location and frequencies of the highest amplitude excitation vary with the strength of the perturbation. These results agree with actual frequency ratios of twin kHz QPOs that have been reported in some black hole candidates, and they may be consistent also with correlation of the twin peaks in Sco X-1.

In 3-d SPH simulations of the coalescence of a quark star with a pseudo-Newtonian black hole all of the quark matter is quickly accreted by the black hole. The Madsen-Caldwell-Friedman argument against the existence of quark stars may need to be re-examined.

Recently, several ways of obtaining observational proof of the existence of black-hole horizons have been proposed. We argue here that such proof is fundamentally impossible: observations can provide arguments, sometimes very strong ones, in favour of the existence of the event horizon, but they cannot prove it. This applies also to future observations, which will trace very accurately the details of the spacetime metric of a body suspected of being a black hole.

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.

The X-ray emission of RXJ1856.5-3754 has been found to coincide to unprecedented accuracy with that of a blackbody, of radius 5.8+-0.9 km for the measured parallax distance of 140 pc (Burwitz et al. 2001, Drake et al. 2002). If the emission is uniform over the whole surface of a non-rotating star, the mass of the star cannot exceed 0.75+-0.12 solar mass regardless of its composition. If the compact object is a quark star described by the MIT-bag equation of state (a ``strange star''), the mass is no more than 0.3 solar mass. Comparably small masses are also obtained for the X-ray bursters Aql X-1 and KS1731-260 for some fits to their spectra.

We apply Feynman's principle, ``The same equations have the same solutions'', to Kepler's problem and show that Newton's dynamics in a properly curved 3-D space is identical with that described by Einstein's theory in the 3-D optical geometry of Schwarzschild's spacetime. For this reason, rather unexpectedly, Newton's formulae for Kepler's problem, in the case of nearly circular motion in a static, spherically spherical gravitational potential accurately describe strong field general relativistic effects, in particular vanishing of the radial epicyclic frequency at the marginally stable orbit.

The spin-up of a neutron star crucially depends on the maximum orbital frequency around it, as do a host of other high energy accretion phenomena in low mass X-ray binaries, including quasi-periodic oscillations (QPOs) in the X-ray flux. We compare the maximum orbital frequencies for MIT-bag quark stars and for neutron stars modeled with the FPS equation of state. The results are based on relativistic calculations of constant baryon sequences of uniformly rotating strange star models, and are presented as a function of the stellar rotational frequency. The marginally stable orbit is present outside quark stars for a wide range of parameters, but outside the neutron stars it is present only for the highest values of mass. This allows a discrimination between quark stars and neutron stars in the resonance theory of kHz QPOs.

Non-linear acoustic coupling of modes in accretion disks allows a parametric resonance between epicyclic motions. In black hole disks, such a resonance first occurs, and is strongest, when the radial and vertical epicyclic frequencies are approximately in a 2:3 ratio, in agreement with the 300 and 450 Hz frequencies reported in GRO J1655-40 and the 184 and 276 Hz frequencies reported in XTE J1550-564. The first overtone in the same resonance is in a 5:3 ratio with the fundamental, in agreement with the 69.2 and 41.5 Hz frequencies reported in GRS 1915+105. The narrow width and the frequency stability of the corresponding QPOs follow from general properties of parametric resonance.

We compute and present the distribution in mass of single and binary neutron stars, strange stars, and black holes. The calculations were performed using a stellar population synthesis code. We follow evolution of massive single stars as well as binaries with high mass primaries. The final product of the latter evolution can be either a binary composed of a white dwarf and a compact object (neutron star, black hole, strange star), two compact objects in a binary, or two single stars if the system was disrupted. We find in binaries a population of black holes which are more massive than single black holes which are a product of either binary or single evolution. We also find that if quark stars exist at all, their population can be as large as the population of black holes.

We present Newtonian three-dimensional hydrodynamical simulations of the merger of quark stars with black holes. The initial conditions correspond to non-spinning stars in Keplerian orbits, the code includes gravitational radiation reaction in the quadrupole approximation for point masses. We find that the quark star is disrupted, forming transient accretion structures around the black hole, but 0.03 of the original stellar mass survives the initial encounter and remains in an elongated orbit as a rapidly rotating quark starlet, in all cases. No resolvable amount of mass is dynamically ejected during the encounters--the black hole eventually accretes 99.99% +/- 0.01% of the quark matter initially present.

Observations of kHz quasi-periodic oscillations (QPOs) in the X-ray fluxes of low-mass X-ray binaries (LMXBs) have been used in attempts to constrain the external metric of the compact members of these binaries, as well as their masses and the equation of state of matter at supranuclear denisties. We compute the maximum orbital frequency of stable circular motion around uniformly rotating strange stars described by the MIT bag model. The calculations are performed for both normal and supramassive constant baryon mass sequences of strange stars rotating at all possible rates. We find the lower limits on the maximum orbital frequency and discuss them for a range of masses and for all rotational frequencies allowed in the model considered. We show that for slowly and moderately rotating strange stars the maximum value of orbital frequency can be a good indicator of the mass of the compact object. However, for rapidly rotating strange stars the same value of orbital frequency in the innermost stable circular orbit is obtained for stars with masses ranging from that of a planetoid to about three solar masses. At sufficiently high rotation rates of the strange star, the rotational period alone constrains the stellar mass to a surprisingly narrow range.

Heyl (2000) made an interesting suggestion that the observed shifts in QPO frequency in type I X-ray bursts could be influenced by the same geometrical effect of strong gravity as the one that causes centrifugal force reversal discovered by Abramowicz and Lasota (1974). However, his main result contains a sign error. Here we derive the correct formula and conclude that constraints on the M(R) relation for neutron stars deduced from the rotational-modulation model of QPO frequency shifts are of no practical interest because the correct formula implies a weak condition R* > 1.3 Rs, where Rs is the Schwarzschild radius. We also argue against the relevance of the rotational-modulation model to the observed frequency modulations.

We note that the recently discovered 450 Hz frequency in the X-ray flux of the black hole candidate GRO J1655-40 is in a 3:2 ratio to the previously known 300 Hz frequency of quasi-periodic oscillations (QPO) in the same source. If the origin of high frequency QPOs in black hole systems is a resonance between orbital and epicyclic motion of accreting matter, as suggested previously, the angular momentum of the black hole can be accurately determined, given its mass. We find that the dimensionless angular momentum is in the range $0.2<j<0.65$ if the mass is in the (corresponding) range of 5.5 to 7.9 solar masses.

We explain the origin of the puzzling high frequency peaks (QPOs) in the variability power spectra of accreting neutron stars and black holes as a non-linear 1:2 or 1:3 resonance between orbital and radial epicyclic motion. These resonances are present because the gravitational field deviates strongly from a Newtonian 1/r potential. Our theory agrees with the recently reported observations of two QPOs, at 300 Hz and 450 Hz, in the black hole candidate GRO J1655-40.

We present the first numerical results on the binary coalescence of a quark star with a black hole, obtained with a 3-D Newtonian smooth particle hydro (SPH) code. The star is initially represented by 17,000 particles modeling a self-gravitating fluid with the equation of state P= (rho - rho0)c**2/3, and the black hole by a point mass with an absorbing boundary at the Schwarzschild radius. As in similar calculations carried out for a stiff polytrope, the stellar core survives the initial episode of mass transfer, but here an accretion disk is clearly formed as well.

We show that in the limit of low masses (M<0.1M solar), the orbital and rotational frequencies in exact relativistic numerical models of rotating, axially symmetric, quark stars coincide with those for Maclaurin spheroids. In particular, when the eccentricity of the spheroid exceeds a critical value (e>0.83458318), circular orbits in the equatorial plane are unstable for a range of orbital radii outside the stellar surface. The orbital frequency in the marginally stable orbit, and all other orbits, around Maclaurin spheroids goes to zero in the limit e=1. Maximum orbital frequencies in accreting Galactic X-ray sources may be a poor indication of the source mass.

It has been suggested that observations of quasi periodic oscillations (QPOs) in the X-ray fluxes from low mass X-ray binaries can be used to constrain the mass of the compact object and the equation of state of its matter. A specific suggestion that the kHz QPO frequency saturates at the maximum orbital frequency has been widely considered. We examine rotating compact stars described by a new model of strange quark matter (Dey et al. 1998). We calculate the maximum orbital frequencies for both normal and supramassive strange stars described by the Dey model, and present these frequencies for sequences of equilibrium models with constant baryon mass. The maximum orbital frequencies for these compact objects are always higher than the kHz QPO frequencies observed to date.

After reviewing the general-relativistic ``gap'' model of accretion, I discuss its relation to the high frequency quasi-periodic oscillations observed in low-mass X-ray binaries. The ``300'' Hz frequency seen in some X-ray bursts may be a relativistic signature of keplerian rotation of the neutron star.

The basic observed properties of neutron stars are reviewed. I suggest that neutron stars in low-mass X-ray binaries are the best of all known sites for testing strong-field effects of general relativity.

We compute and present the distribution in mass of single and binary neutron stars, strange stars, and black holes. The calculations were performed using a stellar population synthesis code. We follow all phases of single and binary evolution, starting from a ZAMS binary and ending in the creation of one compact object (neutron star, black hole, strange star) and a white dwarf, or two compact objects (single or binary). We assume that neutron stars are formed in the collapse of iron/nickel cores in the mass range M0 < M < M1, quark stars in the range M1 < M < M2, and black holes for core masses M > M2 and find that the population of quark stars can easily be as large as the population of black holes, even if there is only a small mass window for their formation.

An analytic solution is presented to the three-dimensional problem of steady axisymmetric fluid flow through an accretion disk. The solution has been obtained through a systematic expansion in the small parameter epsilon =H/R (the ratio of disk thickness to its radial dimension) of the equations of viscous hydrodynamics. The equation of state was assumed to be polytropic. For all values alpha< 0.685 of the viscosity parameter, we find significant backflow in the midplane of the disk occuring at all radii larger than a certain value; however, in the inner regions of the disk the fluid always flows toward the accreting object. The region of backflow is separated from the region of inflow by a surface flaring outwards from a circular locus of stagnation points situated in the midplane of the disk.

One of a family previously proposed ``central engines'' for cosmic gamma-ray burst sources (Klu\'zniak & Ruderman 1998) is considered in some detail. A steadily accreting $10^6$ Gauss magnetic white dwarf should ultimately collapse to a strongly differentially rotating, millisecond-rotation-period neutron star for a wide range of steady accretion rates and initial masses if the accreting white dwarf has an evolved O-Ne-Mg composition. A similar neutron star could also result from an initial C-O white dwarf but only for more constrained accretion rates. Because the collapsing white dwarf begins as a $\gamma=4/3$ polytrope, the final neutron star's spin-rate increases strongly with cylindrical radius. A stable wind-up of the neutron star's poloidal magnetic field then produces buoyant magnetic toroids which grow, break loose, rise, and partly penetrate the neutron star surface to form a transient, $B\approx 10^{17}$ G millisecond spin-period pulsar with a powerful pulsar wind (Usov 1992). This pulsar wind emission is then rapidly suppressed by the surface shear motion from the strong stellar differential rotation. This wind-up and transient pulsar formation can occur at other times on different cylinders and/or repeat on the same one, with (re-)wind up and surface penetration time scales hugely longer han the neutron star's millisecond spin period. In this way, differential rotation both opens and closes the doors which allow neutron star spin-energy to be emitted in powerful bursts of pulsar wind. Predictions of this model compare favorably to needed central engine properties of gamma-ray burst sources (total energy, duration, sub-burst fluctuations and time scales, variability among burst events, and baryon loading).

It has been suggested that the frequency in the co-rotating innermost stable circular orbit (ISCO) about a compact stellar remnant can be determined through X-ray observations of low-mass X-ray binaries, and that its value can be used to constrain the equation of state of ultradense matter. Upon constructing numerical models of rapidly rotating strange (quark) stars in general relativity, we find that for stars rotating at the equatorial mass-shedding limit, the ISCO is indeed above the stellar surface, for a wide range of central energy densities at a height equal to 11% of the circumferential stellar radius, which scales inversely with the square root of the energy density, of self-bound quark matter at zero presure. In contrast to static stars, the ISCO frequencies of rapidly rotating strange stars can be as low as 0.9 kHz for a 1.3 solar mass strange star. Hence, the presence of strange stars in low-mass X-ray binaries cannot be excluded on the basis of the currently observed frequencies of kHz QPOs, such as the cut-off frequency of 1066 Hz in 4U 1820-30.

If Lorentz invariance is broken at an energy scale Eq, as has recently been suggested in the context of attempts to quantize gravity, the kinematics of photon-photon collisions would be profoundly affected at lower energies. Specifically, electron-positron pair creation on soft photons may be forbidden at photon energies as low as 30 TeV times square root of (Eq/10**17 GeV) and the Universe would then be transparent to high energy photons. The proposition that Lorentz invariance is broken may be falsified by the techniques of TeV astronomy.

We present a numerical study of the hydrodynamics in the final stages of inspiral of a black hole-neutron star binary, when the binary separation becomes comparable to the stellar radius. We use a Newtonian three-dimensional Smooth Particle Hydrodynamics (SPH) code, and model the neutron star with a soft (adiabatic index Gamma=5/3) polytropic equation of state and the black hole as a Newtonian point mass which accretes matter via an absorbing boundary at the Schwarzschild radius. Our initial conditions correspond to tidally locked binaries in equilibrium, and we have explored configurations with different values of the mass ratio q=Mns/Mbh, ranging from q=1 to q=0.1. The dynamical evolution is followed for approximately 23 ms, and in every case studied here we find that the neutron star is tidally disrupted on a dynamical timescale, forming a dense torus around the black hole that contains a few tenths of a solar mass. A nearly baryon-free axis is present in the system throughout the coalescence, and only modest beaming of a fireball that could give rise to a gamma-ray burst would be sufficient to avoid excessive baryon contamination. We find that some mass (on the order of 0.01 to 0.001 solar masses) may be dynamically ejected from the system, and could thus contribute substantially to the amount of observed r-process material in the galaxy. We calculate the gravitational radiation waveforms and luminosity in the quadrupole approximation.

We present a numerical study of the hydrodynamics in the final stages of inspiral in a black hole-neutron star binary, when the separation becomes comparable to the stellar radius. We use a Newtonian three-dimensional Smooth Particle Hydrodynamics (SPH) code, and model the neutron star with a soft (adiabatic index Gamma=5/3) polytropic equation of state and the black hole as a Newtonian point mass which accretes matter via an absorbing boundary at the Schwarzschild radius. Our initial conditions correspond to tidally locked binaries in equilibirium. The dynamical evolution is followed for approximately 23 ms, and in every case for Gamma=5/3 we find that the neutron star is tidally disrupted on a dynamical timescale, forming a dense torus around the black hole that contains a few tenths of a solar mass. A nearly baryon-free axis is present in the system throughout the coalescence, a fireball expanding along that axis would avoid excessive baryon contamination and could give rise to a modestly beamed gamma-ray burst.

We used a newtonian SPH (smooth-particle hydrodynamics) code to follow the final stages of evolution of a coalescing binary system of a neutron star and a black hole. We find that the outcome of the "merger" is very sensitive to the equation of state describing the neutron star. A neutron star with a soft equation of state (polytrope with index Gamma=5/3) is completely disrupted, and a fairly large and long-lived accretion torus is formed.

The kilohertz quasi periodic oscillations (QPOs) discovered in several low mass X-ray binaries (LMXBs) by the Rossi X-ray Timing Explorer (XTE) are thought to occur at the orbital frequency in accretion discs whose inner edge corresponds to the innermost (marginally) stable orbit allowed by general relativity. These ideas have been applied to constrain the equation of state (e.o.s.) of the central neutron star. Here we discuss another possibility, that the central object is a strange star, and show how kHz QPOs constrain the e.o.s. of strange matter.

We present a detailed study of the hydrodynamical interactions in a Newtonian black hole-neutron star binary during the last stages of inspiral. We consider close binaries which are tidally locked, use a stiff equation of state (with an adiabatic index Gamma=3) throughout, and explore the effect of different initial mass ratios on the evolution of the system. We calculate the gravitational radiation signal in the quadrupole approximation. Our calculations are carried out using a Smooth Particle Hydrodynamics (SPH) code.

If the ordinary neutrinos oscillate into a sterile flavor in a manner consistent with the Super-Kamiokande data on the zenith-angle dependence of atmospheric mu-neutrino flux, an energy sufficient to power a typical cosmic gamma-ray burst (GRB) (about 10^{52} erg) can be carried by sterile neutrinos away from the source and deposited in a region relatively free of baryons. Hence, ultra-relativistic bulk motion (required by the theory of and observations of GRBs and their afterglows) can easily be achieved in the vicinity of plausible sources of GRBs. Oscillations between sterile and ordinary neutrinos would thus provide a solution to the ``baryon-loading problem'' in the theory of GRBs.

GRBs are thought to arise in relativistic blast-wave shocks at distances of 10 to 1000 AU from the point where the explosive energy is initially released. To account for the observed duration and variability of the gamma-ray emission in most GRBs, a central engine powering the shocks must remain active for several seconds to many minutes but must strongly fluctuate in its output on much shorter timescales. We show how a neutron star differentially rotating at millisecond periods (DROMP) could be such an engine. A magnetized DROMP would repeatedly wind up toroidal magnetic fields to about 10**17 G and only release the corresponding magnetic energy, when each buoyant magnetic field torus floats up to, and breaks through, the stellar surface. The resulting rapid sub-bursts, separated by relatively quiescent phases, repeat until the kinetic energy of differential rotation is exhausted by these events. Calculated values of the energy released and of the various timescales are in agreement with observations of GRBs. The baryon loading in each sub-burst may also be consistent with theoretical requirements for a blast wave capable of giving the X-ray, optical and radio afterglows recently observed from cosmological distances. DROMPs could be created in several kinds of astrophysical events; some of these would be expected to occur at about the observed GRB rate. The requisite differential rotation could be imparted to neutron stars as they are born or at the end of their existence: some DROMPs may be created close to star forming regions while others may arise far from galaxies.

If the observed millisecond variability in the X-ray flux of several neutron-star low-mass X-ray binaries (LMXBs) is interpreted within a general-relativistic framework (Kluzniak, Michelson \& Wagoner 1990) extant at the time of discovery, severe constraints can be placed on the equation of state (e.o.s.) of matter at supranuclear densities. The reported maximum frequency (1.14 +- 0.01 kHz) of quasiperiodic oscillations observed in sources as diverse as Sco X-1 and 4U 1728-34 would imply that the neutron star masses in these LMXBs are M > 1.9 M_solar, and hence many equations of state would be excluded. Among the very few still viable equations of state are the e.o.s. of Phandaripande and Smith (1975), and e.o.s. AV14 + UVII of Wiringa, Fiks \& Fabrocini (1988).

We present the results of Newtonian hydrodynamic simulations of the coalescence of a binary consisting of a black hole with a neutron star. The calculations show that for a wide range of initial conditions the core of the neutron star survives the initial mass transfer episode. We therefore identify black hole-neutron star binaries as the astrophysical production site of low mass neutron stars unstable to explosion. The relevance of the simulations to the theory of gamma-ray bursts is also discussed.

We present the results of hydrodynamic (SPH) simulations showing the coalescence of a black hole with a neutron star to be a promising theoretical source of short duration gamma-ray bursts. The favorable features of the process include rapid onset, millisecond variability, a duration much longer than the dynamical timescale, and a range of outcomes sufficient to allow variety in the properties of individual gamma-ray bursts. Interestingly, the process of coalescence differs rather markedly from past predictions.

Angular momentum loss via the emission of gravitational waves must eventually drive compact binaries containing black holes and/or neutron stars to coalesce. The resulting events are primary candidate sources for detectors such as VIRGO and LIGO. We present calculations of gravitational radiation waveforms and luminosities for the coalescence of a black hole-neutron star binary, performed in the quadrupole approximation using a Newtonian smooth particle hydrodynamics code. We discuss the dependence of the waveforms and the total emitted luminosity as well as the final configuration of the system on the initial mass ratio and the degree of tidal locking.