Evolution of Class 0 protostars: models versus observations

The rates at which mass accumulates into protostellar cores can now be predicted in numerical simulations. Our purpose here is to develop methods to compare the statistical properties of the predicted protostars with the observable parameters. This requires (1) an evolutionary scheme to convert numerically derived mass accretion rates into evolutionary tracks and (2) a technique to compare the tracks to the observed statistics of protostars. Here, we use a 3D Kolmogorov–Smirnov test to quantitatively compare model evolutionary tracks and observations of Class 0 protostars.
We find that the wide range of accretion functions and time-scales associated with gravoturbulent simulations naturally overcome difficulties associated with schemes that use a fixed accretion pattern. This implies that the location of a protostar on an evolutionary track does not precisely determine the present age or final accrued mass. Rather, we find that predictions of the final mass for protostars from observed   T _bol– L _bol  values are uncertain by a factor of 2 and that the bolometric temperature is not always a reliable measure of the evolutionary stage. Furthermore, we constrain several parameters of the evolutionary scheme and estimate a lifetime of Class 0 sources of  2–6 × 10⁴ yr  , which is related to the local free-fall time and thus to the local density at the onset of the collapse. Models with Mach numbers smaller than six are found to best explain the observational data. Generally, only a probability of 70 per cent was found that our models explain the current observations. This is caused by not well-understood selection effects in the observational sample and the simplified assumptions in the models.     

Simultaneous X-ray/optical observations of GX 9+9 (4U 1728−16)

We report on the results of the first simultaneous X-ray ( RXTE ) and optical [South African Astronomical Observatory (SAAO)] observations of the luminous low-mass X-ray binary (LMXB) GX 9+9 in 1999 August. The high-speed optical photometry revealed an orbital period of 4.1958 h and confirmed previous observations, but with greater precision. No X-ray modulation was found at the orbital period. On shorter time-scales, a possible 1.4-h variability was found in the optical light curves which might be related to the MHz quasi-periodic oscillations seen in other LMXBs. We do not find any significant X-ray/optical correlation in the light curves. In X-rays, the colour–colour and hardness-intensity diagrams indicate that the source shows characteristics of an atoll source in the upper banana state, with a correlation between intensity and spectral hardness. Time-resolved X-ray spectroscopy suggests that two-component spectral models give a reasonable fit to the X-ray emission. Such models consist of a blackbody component which can be interpreted as the emission from an optically thick accretion disc or an optically thick boundary layer, and a hard Comptonized component for an extended corona.     

Evolution of the phase-space density of dark matter haloes and mixing effects in merger events

Cosmological N -body simulations were performed to study the evolution of the phase-space density   Q =ρ/σ³  of dark matter haloes. No significant differences in the scale relations   Q ∝σ^−2.1  or   Q ∝ M ^−0.82  are seen for the 'cold' or 'warm' dark matter models. The follow-up of individual haloes from   z = 10  up to the present time indicate the existence of two main evolutionary phases: an early and fast one  (10 > z > 6.5)  , in which Q decreases on the average by a factor of 40 as a consequence of the randomization of bulk motions, and a late and long one  (6.5 > z ≥ 0)  , in which Q decreases by a factor of 20 because of mixing induced by merger events. The study of these haloes has also evidenced that rapid and positive variations of the velocity dispersion, induced by merger episodes, are related to a fast decrease of the phase-space density Q .     

Substellar companions and isolated planetary-mass objects from protostellar disc fragmentation

Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time-scale that is short compared with the orbital period. We use a combination of hydrodynamic simulations and N -body orbit integrations to study the long-term evolution of a fragmenting disc with an initial mass ratio to the star of   M _disc/ M _*= 0.1  . For a disc that is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar-mass star) up to  ≈0.01 M_⊙  . Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 – in which the companions are close to or beyond the deuterium burning limit – appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ∼1 au should differ from that of stars with lower mass planetary companions.     

Non-linear bending waves in Keplerian accretion discs

The non-linear dynamics of a warped accretion disc is investigated in the important case of a thin Keplerian disc with negligible viscosity and self-gravity. A one-dimensional evolutionary equation is formally derived that describes the primary non-linear and dispersive effects on propagating bending waves other than parametric instabilities. It has the form of a derivative non-linear Schrödinger (DNLS) equation with coefficients that are obtained explicitly for a particular model of a disc. The properties of this equation are analysed in some detail and illustrative numerical solutions are presented. The non-linear and dispersive effects both depend on the compressibility of the gas through its adiabatic index Γ. In the physically realistic case Γ < 3, non-linearity does not lead to the steepening of bending waves but instead enhances their linear dispersion. In the opposite case Γ > 3, non-linearity leads to wave steepening and solitary waves are supported. The effects of a small effective viscosity, which may suppress parametric instabilities, are also considered. This analysis may provide a useful point of comparison between theory and numerical simulations of warped accretion discs.     

A transient relativistic radio jet from Cygnus X-1

We report the first observation of a transient relativistic jet from the canonical black hole candidate, Cygnus X-1, obtained with the Multi-Element Radio-Linked Interferometer Network (MERLIN). The jet was observed in only one of six epochs of MERLIN imaging of the source during a phase of repeated X-ray spectral transitions in 2004 Jan–Feb, and this epoch corresponded to the softest 1.5–12 keV X-ray spectrum. With only a single epoch revealing the jet, we cannot formally constrain its velocity. Nevertheless, several lines of reasoning suggest that the jet was probably launched 0.5–4.0 d before this brightening, corresponding to projected velocities of  0.2 c ≲ v _app≲ 1.6 c   , and an intrinsic velocity of  ≳0.3 c   . We also report the occurrence of a major radio flare from Cyg X-1, reaching a flux density of ∼120 mJy at 15 GHz, and yet not associated with any resolvable radio emission, despite a concerted effort with MERLIN. We discuss the resolved jet in terms of the recently proposed 'unified model' for the disc–jet coupling in black hole X-ray binaries, and tentatively identify the 'jet line' for Cyg X-1. The source is consistent with the model in the sense that a steady jet appears to persist initially when the X-ray spectrum starts softening, and that once the spectral softening is complete the core radio emission is suppressed and transient ejecta/shock observed. However, there are some anomalies, and Cyg X-1 clearly does not behave like a normal black hole transient in progressing to the canonical soft/thermal state once the ejection event has happened.