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.

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.

X-ray burst spectra have long been used to estimate neutron star masses and radii. These estimates assumed that burst spectra are accurately described by the model atmosphere spectra developed over the last three decades. We compared RXTE data from a superburst with these spectra and found that the spectra predicted by previously published model atmospheres are strongly inconsistent with these high-precision measurements. In contrast, a simple Bose-Einstein spectrum is fully consistent with the data, as are recently published model atmosphere spectra. We discuss the implications of our results for determinations of neutron star masses and radii via constraints on their surface gravity and redshift, as originally suggested by Majczyna and Madej.

It has been known for nearly three decades that the energy spectra of thermonuclear X-ray bursts are often well-fit by Planck functions with temperatures so high that they imply a super-Eddington radiative flux at the emitting surface, even during portions of bursts when there is no evidence of photospheric radius expansion. This apparent inconsistency is usually set aside by assuming that the flux is actually sub-Eddington and that the fitted temperature is so high because the spectrum has been distorted by the energy-dependent opacity of the atmosphere. Here we show that the spectra predicted by currently available conventional atmosphere models appear incompatible with the highest-precision measurements of burst spectra made using the Rossi X-ray Timing Explorer, such as during the 4U 1820-30 superburst and a long burst from GX 17+2. In contrast, these measurements are well-fit by Bose-Einstein spectra with high temperatures and modest chemical potentials. Such spectra are very similar to Planck spectra. They imply surface radiative fluxes more than a factor of three larger than the Eddington flux. We find that segments of many other bursts from many sources are well-fit by similar Bose-Einstein spectra, suggesting that the radiative flux at the emitting surface also exceeds the Eddington flux during these segments. We suggest that burst spectra can closely approximate Bose-Einstein spectra and have fluxes that exceed the Eddington flux because they are formed by Comptonization in an extended, low-density radiating gas supported by the outward radiation force and confined by a tangled magnetic field.

We investigate further a model of the accreting millisecond X-ray pulsars we proposed earlier. In this model, the X-ray-emitting regions of these pulsars are near their spin axes but move. This is to be expected if the magnetic poles of these stars are close to their spin axes, so that accreting gas is channeled there. As the accretion rate and the structure of the inner disk vary, gas is channeled along different field lines to different locations on the stellar surface, causing the X-ray-emitting areas to move. We show that this "nearly aligned moving spot model" can explain many properties of the accreting millisecond X-ray pulsars, including their generally low oscillation amplitudes and nearly sinusoidal waveforms; the variability of their pulse amplitudes, shapes, and phases; the correlations in this variability; and the similarity of the accretion- and nuclear-powered pulse shapes and phases in some. It may also explain why accretion-powered millisecond pulsars are difficult to detect, why some are intermittent, and why all detected so far are transients. This model can be tested by comparing with observations the waveform changes it predicts, including the changes with accretion rate.

We have shown previously that many of the properties of persistent accretion-powered millisecond pulsars can be understood if their X-ray emitting areas are near their spin axes and move as the accretion rate and structure of the inner disk vary. Here we show that this "nearly aligned moving spot model" may also explain the intermittent accretion-powered pulsations that have been detected in three weakly magnetic accreting neutron stars. We show that movement of the emitting area from very close to the spin axis to about 10 degrees away can increase the fractional rms amplitude from less than about 0.5 percent, which is usually undetectable with current instruments, to a few percent, which is easily detectable. The second harmonic of the spin frequency usually would not be detected, in agreement with observations. The model produces intermittently detectable oscillations for a range of emitting area sizes and beaming patterns, stellar masses and radii, and viewing directions. Intermittent oscillations are more likely in stars that are more compact. In addition to explaining the sudden appearance of accretion-powered millisecond oscillations in some neutron stars with millisecond spin periods, the model explains why accretion-powered millisecond oscillations are relatively rare and predicts that the persistent accretion-powered millisecond oscillations of other stars may become undetectable for brief intervals. It suggests why millisecond oscillations are frequently detected during the X-ray bursts of some neutron stars but not others and suggests mechanisms that could explain the occasional temporal association of intermittent accretion-powered oscillations with thermonuclear X-ray bursts.

Detection of paired kilohertz quasi-periodic oscillations (kHz QPOs) in the X-ray emission of a compact object is compelling evidence that the object is an accreting neutron star. In many neutron stars, the stellar spin rate is equal or roughly equal to Delta-nu, the frequency separation of the QPO pair, or to 2Delta-nu. Hence, if the mechanism that produces the kilohertz QPOs is similar in all stars, measurement of Delta-nu can provide an estimate of the star's spin rate. The involvement of the stellar spin in producing Delta-nu indicates that the magnetic fields of these stars are dynamically important. We focus here on the implications of the paired kHz QPOs recently discovered in the low-mass X-ray binary (LMXB) system Cir X-1 (Boutloukos et al. 2006). The kHz QPOs discovered in Cir X-1 are generally similar to those seen in other stars, establishing that the compact object in the Cir X-1 system is a neutron star. However, the frequency nu-u of its upper kHz QPO is up to a factor of three smaller than is typical, and Delta-nu varies by about a factor 2 (167 Hz, the largest variation so far observed). Periodic oscillations have not yet been detected from Cir X-1, so its spin rate has not yet been measured directly. The low values of nu-u and the large variation of Delta-nu challenge current models of the generation of kHz QPOs. Improving our understanding of Cir X-1 will improve our knowledge of the spin rates and magnetic fields of all neutron stars.

Using the Rossi X-ray Timing Explorer (RossiXTE), astronomers have discovered that disk-accreting neutron stars with weak magnetic fields produce three distinct types of high-frequency X-ray oscillations. These oscillations are powered by release of the binding energy of matter falling into the strong gravitational field of the star or by the sudden nuclear burning of matter that has accumulated in the outermost layers of the star. The frequencies of the oscillations reflect the orbital frequencies of gas deep in the gravitational field of the star and/or the spin frequency of the star. These oscillations can therefore be used to explore fundamental physics, such as strong-field gravity and the properties of matter under extreme conditions, and important astrophysical questions, such as the formation and evolution of millisecond pulsars. Observations using RossiXTE have shown that some two dozen neutron stars in low-mass X-ray binary systems have the spin rates and magnetic fields required to become millisecond radio-emitting pulsars when accretion ceases, but that few have spin rates above about 600 Hz. The properties of these stars show that the paucity of spin rates greater than 600 Hz is due in part to the magnetic braking component of the accretion torque and to the limited amount of angular momentum that can be accreted in such systems. Further study will show whether braking by gravitational radiation is also a factor. Analysis of the kilohertz oscillations has provided the first evidence for the existence of the innermost stable circular orbit around dense relativistic stars that is predicted by strong-field general relativity. It has also greatly narrowed the possible descriptions of ultradense matter.

We discovered kHz QPOs in 80 archived RXTE observations from the peculiar low-mass X-ray binary (LMXB) Circinus X-1. In 11 cases these appear in pairs in the frequency range of ~230 Hz to ~500 Hz for the upper kHz QPO and ~56 Hz to ~225 Hz for the lower kHz QPO. Their correlation with each other, which is similar to that of frequencies of kHz QPO pairs in other LMXBs containing a neutron star, and their variation by a factor two confirm that the central object is a neutron star. These are the lowest frequencies of kHz QPO pairs discovered so far and extend the above correlation over a frequency range of factor four. In this new frequency range the frequency difference of the two kHz QPOs increases monotonically by more than ~170 Hz with increasing kHz QPO frequency, challenging theoretical models.

We use perturbation theory and the relativistic Cowling approximation to numerically compute characteristic oscillation modes of rapidly rotating relativistic stars which consist of a perfect fluid obeying a polytropic equation of state. We present a code that allows the computation of modes of arbitrary order. We focus here on the overall distribution of frequencies. As expected, we find an infinite pressure mode spectrum extending to infinite frequency. In addition we obtain an infinite number of inertial mode solutions confined to a finite, well-defined frequency range which depends on the compactness and the rotation frequency of the star. For non-axisymmetric modes we observe how this range is shifted with respect to the axisymmetric ones, moving towards negative frequencies and thus making all m>2 modes unstable. We discuss whether our results indicate that the star's spectrum must have a continuous part, as opposed to simply containing an infinite number of discrete modes.

We report the discovery with the RXTE/PCA of twin kHz QPOs in Cir X-1. Eleven cases of simultaneous double QPOs occurred, with significances of up to 6.3 and 5.5 sigma and centroid frequencies ranging between approximately 56-225 and 230-500 Hz for the two QPO peaks, respectively, i.e., for the most part at frequencies well below those of other sources. The QPO properties clearly indicate that these double peaks are the kHz QPOs known from low magnetic field neutron stars, and not black-hole high-frequency QPOs, confirming that Cir X-1 is a neutron star. The kHz QPO peak separation varies over a wide range, \~175-340 Hz, and increases with QPO frequency. This is contrary to what is seen in other sources but agrees with predictions of the relativistic precession model and Alfven wave models; beat-frequency models require modification to accommodate this. In other observations single kHz QPOs can be seen down to frequencies as low as ~12 Hz, as well as a strong low-frequency (LF) QPO between 1 and 30 Hz. The relations between the frequencies of the kHz QPOs and the LF QPO are in good agreement with those found previously in Z sources, confirming that Cir X-1 may be a peculiar Z source. We suggest that the low frequencies of the kHz QPOs in Cir X-1 and to a lesser extent in (other) Z sources might be due to a relatively stronger radial inflow to the neutron star than in other kHz QPO sources.

If the cluster formation rate is constant and the disruption time of clusters depends on their initial mass as t_dis= t_4 x (M_i/10^4 Msun)^gamma, where M_i is the initial mass of a cluster, t_4 is the disruption time of a cluster of M_i=10^4 Msun and gamma is an exponent, then the values of t_4 and gamma can be derived in a very simple way from the age and mass histograms of large homogeneous samples of clusters with reliable age and/or mass determinations. We demonstrate the method and derive the values of t_4 and gamma from cluster samples in selected fields of four galaxies: M51, M33, SMC and the solar neighbourhood. The values of the exponent gamma are the same within their uncertainties in the four galaxies with gamma = 0.62 +- 0.06. However, the disruption time t_4 of clusters of 10^4 Msun is very different in the different locations. It is shortest in the inner region of M51 (t_4 = 40 Myr) and highest in the SMC (t_4=8 Gyr).

We present a new method to derive the cluster disruption time in selected regions of galaxies from the mass or age distribution of magnitude-limited cluster samples. If the disruption time of clusters in a region of a galaxy depends on their initial mass as t_4 x (M_cluster/10^4 M_sun)^gamma and if the cluster formation rate is constant, then the mass and age distributions of the observed clusters will each show two powerlaw relations. The values of t_4 and gamma can be derived from these relations. We used this method to derive the cluster disruption time in specific regions in four galaxies: the inner region of M51, a region of M33, the SMC and the solar neighbourhood. The values of gamma are the same in the four galaxies within the uncertainty and the mean value is gamma= 0.62 +- 0.06. However the disruption time t_4 of a cluster of 10^4 M_sun is very different in the different galaxies. The clusters in the SMC have the longest disruption time, t_4 = 8 Gyr, and the clusters at 1 to 3 kpc from the nucleus of M51 have the shortest disruption time of t_4 = 0.04 Gyr. The disruption time of clusters 1 to 5 kpc from the nucleus of M33 is t_4 = 0.13 Gyr and for clusters within 1 kpc from the Sun we find t_4 = 1.0 Gyr.