The ages of pre-main-sequence stars

The position of pre-main-sequence or protostars in the Hertzsprung–Russell diagram is often used to determine their mass and age by comparison with pre-main-sequence evolution tracks. On the assumption that the stellar models are accurate, we demonstrate that, if the metallicity is known, the mass obtained is a good estimate. However, the age determination can be very misleading, because it is significantly (generally different by a factor of 2 to 5) dependent on the accretion rate and, for ages less than about 10⁶ yr, the initial state of the star. We present a number of accreting protostellar tracks that can be used to determine age if the initial conditions can be determined and the underlying accretion rate has been constant in the past. Because of the balance established between the Kelvin–Helmholtz, contraction time-scale and the accretion time-scale, a pre-main-sequence star remembers its accretion history. Knowledge of the current accretion rate, together with an HR-diagram position, gives information about the rate of accretion in the past, but does not necessarily improve any age estimate. We do not claim that ages obtained by comparison with these particular accreting tracks are likely to be any more reliable than those from comparisons with non-accreting tracks. Instead, we stress the unreliability of any such comparisons, and use the disparities between various tracks to estimate the likely errors in age and mass estimates. We also show how a set of coeval accreting objects do not appear coeval when compared with non-accreting tracks. Instead, accreting pre-main-sequence stars of around a solar mass are likely to appear older than those of either smaller or larger mass.     

Variable accretion and emission from the stellar winds in the Galactic Centre

We present numerical simulations of stellar wind dynamics in the central parsec of the Galactic Centre, studying in particular the accretion of gas on to Sgr A*, the supermassive black hole. Unlike our previous work, here we use state-of-the-art observational data on orbits and wind properties of individual wind-producing stars. Since wind velocities were revised upwards and non-zero eccentricities were considered, our new simulations show fewer clumps of cold gas and no conspicuous disc-like structure. The accretion rate is dominated by a few close 'slow-wind stars' ( v _w≤ 750 km s⁻¹), and is consistent with the Bondi estimate, but variable on time-scales of tens to hundreds of years. This variability is due to the stochastic infall of cold clumps of gas, as in earlier simulations, and to the eccentric orbits of stars. The present models fail to explain the high luminosity of Sgr A* a few hundred years ago implied by Integral observations, but we argue that the accretion of a cold clump with a small impact parameter could have caused it. Finally, we show the possibility of constraining the total mass-loss rate of the 'slow-wind stars' using near infrared observations of gas in the central few arcseconds.     

Magnetospheric accretion and pre-main-sequence stellar masses

We present a method of determining lower limits on the masses of pre-main-sequence (PMS) stars and so constraining the PMS evolutionary tracks. This method uses the redshifted absorption feature observed in some emission-line profiles of T Tauri stars, indicative of infall. The maximum velocity of the accreting material measures the potential energy at the stellar surface, which, combined with an observational determination of the stellar radius, yields the stellar mass. This estimate is a lower limit owing to uncertainties in the geometry and projection effects. Using available data, we show that the computed lower limits can be larger than the masses derived from PMS evolutionary tracks for M   0.5 M. Our analysis also supports the notion that accretion streams do not impact near the stellar poles but probably hit the stellar surface at moderate latitudes.     

Accretion and dynamical interactions in small-N star-forming clusters: N= 5

We present results from high-resolution hydrodynamical simulations that explore the effects of small-scale clustering in star-forming regions. A large ensemble of small- N clusters with five stellar seeds have been modelled and the resulting properties of stars and brown dwarfs statistically derived and compared with observational data.
Close dynamical interactions between the protostars and competitive accretion driven by the cloud collapse are shown to produce a distribution of final masses that is bimodal, with most of the mass residing in the binary components. When convolved with a suitable core mass function, the final distribution of masses resembles the observed initial mass function, in both the stellar and substellar regimes. Binaries and single stars are found to constitute two kinematically distinct populations, with about half of the singles attaining velocities ≥2 km s⁻¹, which might deprive low-mass star-forming regions of their lightest members in a few crossing times. The eccentricity distribution of binaries and multiples is found to follow a distribution similar to that of observed long-period (uncircularized) binaries.
The results obtained support a mechanism in which a significant fraction of brown dwarfs form under similar circumstances as those of normal stars but are ejected from the common envelope of unstable multiple systems before their masses exceed the hydrogen burning limit. We predict that many close binary stars should have wide brown dwarf companions. Brown dwarfs, and, in general, very low-mass stars, would be rare as pure binary companions. The binary fraction should be a decreasing function of primary mass, with low-mass or substellar primaries being scarce. Where such binaries exist, they are expected either to be close enough (semimajor axis ∼10 au) to survive strong interactions with more massive binaries or to be born in very small molecular cloud cores.     

A submillimetre survey for protostellar accretion discs using the JCMT–CSO interferometer

We present a submillimetre continuum survey for accretion discs around seven embedded protostars in the Perseus and Serpens molecular clouds. Observations were made at frequencies between 339 and 357 GHz using the James Clerk Maxwell Telescope–Caltech Submillimeter Observatory single-baseline interferometer on Mauna Kea, Hawaii. All the objects in our survey show compact dust emission on scales ≲1 arcsec, assumed to arise in a circumstellar accretion disc. We compare the properties of this compact component with evolutionary indicators, such as the ratio of compact to extended emission, and bolometric temperature. We find that discs of mass ∼0.01 M_⊙ have formed by the Class 0 stage, and that similar mass discs are observed in Class I and Class II sources. A trend is observed whereby the ratio of compact to extended emission in our sources increases from Class 0 to Class II sources. For three of the objects in the survey, NGC 1333 IRAS2:CR1 and SVS13 in Perseus, and FIRS1 in Serpens, the signal-to-noise ratio is sufficient to allow us to model the brightness distributions with elliptical Gaussian and power-law disc models. The Gaussian fits give semimajor half-power radii of approximately 90 to 140 au, at the assumed distance of 350 pc to the Perseus and Serpens clouds.     

SPH simulations of negative (nodal) superhumps: a parametric study

Negative superhumps in cataclysmic variable systems result when the accretion disc is tilted with respect to the orbital plane. The line of nodes of the tilted disc precesses slowly in the retrograde direction, resulting in a photometric signal with a period slightly less than the orbital period. We use the method of smoothed particle hydrodynamics to simulate a series of models of differing mass ratio and effective viscosity to determine the retrograde precession period and superhump period deficit  ɛ₋  as a function of system mass ratio q . We tabulate our results and present fits to both  ɛ₋  and  ɛ₊  versus q , as well as compare the numerical results with those compiled from the literature of negative superhump observations. One surprising result is that while we find negative superhumps most clearly in simulations with an accretion stream present, we also find evidence for negative superhumps in simulations in which we shut off the mass transfer stream completely, indicating that the origin of the photometric signal is more complicated than previously believed.     

The anticorrelation between the hard X-ray photon index and the Eddington ratio in low-luminosity active galactic nuclei

We find a significant anticorrelation between the hard X-ray photon index Γ and the Eddington ratio   L _bol/ L _Edd  for a sample of low-ionization nuclear emission-line regions and local Seyfert galaxies, compiled from literature with Chandra or XMM–Newton observations. This result is in contrast with the positive correlation found in luminous active galactic nuclei (AGN), while it is similar to that of X-ray binaries (XRBs) in the low/hard state. Our result is qualitatively consistent with the spectra produced from advection-dominated accretion flows (ADAFs). It implies that the X-ray emission of low-luminosity active galactic nuclei (LLAGN) may originate from the Comptonization process in ADAF, and the accretion process in LLAGN may be similar to that of XRBs in the low/hard state, which is different from that in luminous AGN.     

Optical spectroscopy and Doppler tomography of Cygnus X-2

We present phase resolved optical spectroscopy and Doppler tomography of V1341 Cygni, the optical counterpart to the neutron star low-mass X-ray binary (LMXB) Cygnus X-2 (Cyg X-2). We derive a radial velocity (RV) curve for the secondary star, finding a projected RV semi-amplitude of   K ₂= 79 ± 3 km s⁻¹  , leading to a mass function of  0.51 ± 0.06 M_⊙, ∼30  per cent lower than the previous estimate. We tentatively attribute the lower value of K ₂ (compared to that obtained by other authors) to variations in the X-ray irradiation of the secondary star at different epochs of observations. The limited phase coverage and/or longer timebase of previous observations may also contribute to the difference in K ₂. Our value for the mass function implies a primary mass of  1.5 ± 0.3 M_⊙  , somewhat lower than previous dynamical estimates, but consistent with the value found by analysis of type-I X-ray bursts from this system. Our Doppler tomography of the broad He  ii λ4686 line reveals that most of the emission from this line is produced on the irradiated face of the donor star, with little emission from the accretion disc. In contrast, the Doppler tomogram of the N  iii λ4640.64 Bowen blend line shows bright emission from near the gas stream/accretion disc impact region, with fainter emission from the gas stream and secondary star. This is the first LMXB for which the Bowen blend is dominated by emission from the gas stream/accretion disc impact region, without comparable emission from the secondary star. This has implications for the interpretation of Bowen blend Doppler tomograms of other LMXBs for which the ephemeris may not be accurately known.     

Searching for proto‐brown dwarfs: Extending near IR spectroscopy of protostars below the hydrogen burning limit

Recent observations of nearby star forming regions have offered evidence that young brown dwarfs undergo a period of mass accretion analogous to the T Tauri phase observed in young stars. Brown dwarf analogs to stellar protostars, however, have yet to be definitively observed. These young, accreting objects would shed light on the nature of the dominant brown dwarf formation process, as well as provide ideal laboratories to investigate the dependence of the accretion mechanism on protostellar mass. Recent near infrared surveys have identified candidate proto‐brown dwarfs and characterized low mass protostars in nearby star forming regions. These techniques allow near infrared spectra to diagnose the effective temperature, accretion luminosity, magnetic field strength and rotation velocity of young low mass stars across the stellar/substellar boundary. The lowest mass proto‐brown dwarfs (M < 40 M_Jup), however, will prove challenging to observe given current near IR observational capabilities. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Rotation and activity in the solar-type stars of NGC 2547

We present high-resolution spectroscopy of a sample of 24 solar-type stars in the young (15–40 Myr), open cluster NGC 2547. We use our spectra to confirm cluster membership in 23 of these stars, to determine projected equatorial velocities and chromospheric activity, and to search for the presence of accretion discs. We find examples of both fast ( v _e sin  i >50 km s⁻¹) and slow ( v _e sin  i <10 km s⁻¹) rotators, but no evidence for active accretion in any of the sample. The distribution of projected rotation velocities is indistinguishable from the slightly older IC 2391 and IC 2602 clusters, implying similar initial angular momentum distributions and circumstellar disc lifetimes. The presence of very slow rotators indicates either that long (10–40 Myr) disc lifetimes or internal differential rotation are needed, or that NGC 2547 (and IC 2391/2602) were born with more slowly rotating stars than are presently seen in even younger clusters and associations. The solar-type stars in NGC 2547 follow a similar rotation–activity relationship to that seen in older clusters. X-ray activity increases until a saturation level is reached for v _e sin  i >15–20 km s⁻¹. We are unable to explain why this saturation level, of log( L _x L _bol)≃−3.3, is a factor of 2 lower than in other clusters, but rule out anomalously slow rotation rates or uncertainties in X-ray flux calculations.     

High-resolution simulations of clump–clump collisions using SPH with particle splitting

We investigate, by means of numerical simulations, the phenomenology of star formation triggered by low-velocity collisions between low-mass molecular clumps. The simulations are performed using a smoothed particle hydrodynamics code which satisfies the Jeans condition by invoking on-the-fly particle splitting.
Clumps are modelled as stable truncated (non-singular) isothermal, i.e. Bonnor–Ebert, spheres. Collisions are characterized by M ₀ (clump mass), b (offset parameter, i.e. ratio of impact parameter to clump radius) and     (Mach number, i.e. ratio of collision velocity to effective post-shock sound speed). The gas subscribes to a barotropic equation of state, which is intended to capture (i) the scaling of pre-collision internal velocity dispersion with clump mass, (ii) post-shock radiative cooling and (iii) adiabatic heating in optically thick protostellar fragments.
The efficiency of star formation is found to vary between 10 and 30 per cent in the different collisions studied and it appears to increase with decreasing M ₀, and/or decreasing b , and/or increasing     . For   b < 0.5  collisions produce shock-compressed layers which fragment into filaments. Protostellar objects then condense out of the filaments and accrete from them. The resulting accretion rates are high,     , for the first     . The densities in the filaments,     , are sufficient that they could be mapped in NH₃ or CS line radiation, in nearby star formation regions.     

The X-ray and extreme-ultraviolet flux evolution of SS Cygni throughout outburst

We present the most complete multiwavelength coverage of any dwarf nova outburst: simultaneous optical, Extreme Ultraviolet Explorer and Rossi X-ray Timing Explorer observations of SS Cygni throughout a narrow asymmetric outburst. Our data show that the high-energy outburst begins in the X-ray waveband 0.9–1.4 d after the beginning of the optical rise and 0.6 d before the extreme-ultraviolet rise. The X-ray flux drops suddenly, immediately before the extreme-ultraviolet flux rise, supporting the view that both components arise in the boundary layer between the accretion disc and white dwarf surface. The early rise of the X-ray flux shows that the propagation time of the outburst heating wave may have been previously overestimated.
The transitions between X-ray and extreme-ultraviolet dominated emission are accompanied by intense variability in the X-ray flux, with time-scales of minutes. As detailed by Mauche & Robinson, dwarf nova oscillations are detected throughout the extreme-ultraviolet outburst, but we find they are absent from the X-ray light curve.
X-ray and extreme-ultraviolet luminosities imply accretion rates of  3 × 10¹⁵ g s⁻¹  in quiescence,  1 × 10¹⁶ g s⁻¹  when the boundary layer becomes optically thick, and  ∼10¹⁸ g s⁻¹  at the peak of the outburst. The quiescent accretion rate is two and a half orders of magnitude higher than predicted by the standard disc instability model, and we suggest this may be because the inner accretion disc in SS Cyg is in a permanent outburst state.     

The violent past of Cygnus X-2

Cygnus X-2 appears to be the descendant of an intermediate-mass X-ray binary (IMXB). Using Mazzitelli's stellar code we compute detailed evolutionary sequences for the system and find that its prehistory is sensitive to stellar input parameters, in particular the amount of core overshooting during the main-sequence phase. With standard assumptions for convective overshooting a case B mass transfer starting with a 3.5-M_⊙ donor star is the most likely evolutionary solution for Cygnus X-2. This makes the currently observed state rather short-lived, of order 3 Myr, and requires a formation rate > 10⁻⁷–10⁻⁶ yr⁻¹ of such systems in the Galaxy. Our calculations show that neutron star IMXBs with initially more massive donors (≳4 M_⊙) encounter a delayed dynamical instability; they are unlikely to survive this rapid mass transfer phase. We determine limits for the age and initial parameters of Cygnus X-2 and calculate possible dynamical orbits of the system in a realistic Galactic potential, given its observed radial velocity. We find trajectories which are consistent with a progenitor binary on a circular orbit in the Galactic plane inside the solar circle that received a kick velocity ≤200 km s⁻¹ at the birth of the neutron star. The simulations suggest that about 7 per cent of IMXBs receiving an arbitrary kick velocity from a standard kick velocity spectrum would end up in an orbit similar to Cygnus X-2, while about 10 per cent of them reach yet larger Galactocentric distances.     

The origin of the initial mass function and its dependence on the mean Jeans mass in molecular clouds

We investigate the dependence of stellar properties on the mean thermal Jeans mass in molecular clouds. We compare the results from the two largest hydrodynamical simulations of star formation to resolve the fragmentation process down to the opacity limit, the first of which was reported by Bate, Bonnell & Bromm. The initial conditions of the two calculations are identical except for the radii of the clouds, which are chosen so that the mean densities and mean thermal Jeans masses of the clouds differ by factors of 9 and 3, respectively.
We find that the denser cloud, with the lower mean thermal Jeans mass, produces a higher proportion of brown dwarfs and has a lower characteristic (median) mass of the stars and brown dwarfs. This dependence of the initial mass function (IMF) on the density of the cloud may explain the observation that the Taurus star-forming region appears to be deficient in brown dwarfs when compared with the Orion Trapezium cluster. The new calculation also produces wide binaries (separations >20 au), one of which is a wide binary brown dwarf system.
Based on the hydrodynamical calculations, we develop a simple accretion/ejection model for the origin of the IMF. In the model, all stars and brown dwarfs begin with the same mass (set by the opacity limit for fragmentation) and grow in mass until their accretion is terminated stochastically by their ejection from the cloud through dynamically interactions. The model predicts that the main variation of the IMF in different star-forming environments should be in the location of the peak (due to variations in the mean thermal Jeans mass of the cloud) and in the substellar regime. However, the slope of the IMF at high masses may depend on the dispersion in the accretion rates of protostars.     

The aperiodic broad-band X-ray variability of Cygnus X-3

We study the soft X-ray variability of Cygnus X-3. By combining data from the All-Sky Monitor and Proportional Counter Array instruments on the RXTE satellite with EXOSAT /Medium Energy (ME) detector observations, we are able to analyse the power density spectrum (PDS) of the source from 10⁻⁹ to 0.1 Hz, thus covering time-scales from seconds to years. As the data on the longer time-scales are unevenly sampled, we combine traditional power spectral techniques with simulations to analyse the variability in this range. The PDS at higher frequencies  (≳10⁻³ Hz)  are for the first time compared for all states of this source. We find that it is for all states well described by a power law, with index  ∼−2  in the soft states and a tendency for a less steep power law in the hard state. At longer time-scales, we study the effect of the state transitions on the PDS, and find that the variability below  ∼10⁻⁷ Hz  is dominated by the transitions. Furthermore, we find no correlation between the length of a high/soft-state episode and the time since the previous high/soft state. On intermediate time-scales, we find evidence for a break in the PDS at time-scales of the order of the orbital period. This may be interpreted as evidence for the existence of a tidal resonance in the accretion disc around the compact object, and constraining the mass ratio to   M ₂/ M ₁≲ 0.3  .     

Resistive relaxation of a magnetically confined mountain on an accreting neutron star

Three-dimensional numerical magnetohydrodynamic (MHD) simulations are performed to investigate how a magnetically confined mountain on an accreting neutron star relaxes resistively. No evidence is found for non-ideal MHD instabilities on a short time-scale, such as the resistive ballooning mode or the tearing mode. Instead, the mountain relaxes gradually as matter is transported across magnetic surfaces on the diffusion time-scale, which evaluates to  τ_I∼ 10⁵–10⁸ yr  (depending on the conductivity of the neutron star crust) for an accreted mass of   M _a= 1.2 × 10⁻⁴ M_⊙  . The magnetic dipole moment simultaneously re-emerges as the screening currents dissipate over  τ_I  . For non-axisymmetric mountains, ohmic dissipation tends to restore axisymmetry by magnetic reconnection at a filamentary neutral sheet in the equatorial plane. Ideal-MHD oscillations on the Alfvén time-scale, which can be excited by external influences, such as variations in the accretion torque, compress the magnetic field and hence decrease  τ_I  by one order of magnitude relative to its standard value (as computed for the static configuration). The implications of long-lived mountains for gravitational wave emission from low-mass X-ray binaries are briefly explored.     

Inflow and outflow from the accretion disc of the microquasar SS 433: UKIRT spectroscopy

A succession of near-infrared (near-IR) spectroscopic observations, taken nightly throughout an entire cycle of SS 433's orbit, reveal (i) the persistent signature of SS 433's accretion disc, having a rotation speed of  ∼500 km s⁻¹  , (ii) the presence of circumbinary disc recently discovered at optical wavelengths by Blundell, Bowler & Schmidtobreick (2008) and (iii) a much faster outflow than has previously been measured for the disc wind, with a terminal velocity of  ∼1500 km s⁻¹  . The increased wind terminal velocity results in a mass-loss rate of  ∼10⁻⁴ M_⊙ yr⁻¹  . These, together with the newly (upwardly) determined masses for the components of the SS 433 system, result in an accurate diagnosis of the extent to which SS 433 has super-Eddington flows. Our observations imply that the size of the companion star is comparable with the semiminor axis of the orbit which is given by     , where e is the eccentricity. Our relatively spectral resolution at these near-IR wavelengths has enabled us to deconstruct the different components that comprise the Brackett-γ (Brγ) line in this binary system, and their physical origins. With this line being dominated throughout our series of observations by the disc wind, and the accretion disc itself being only a minority (∼15 per cent) contribution, we caution against use of the unresolved Brγ line intensity as an 'accretion signature' in X-ray binaries or microquasars in any quantitative way.     

On the effect of overshooting as predicted by the modelling of the pre-main-sequence evolution of a 2 M⊙ star

We discuss the effects of convective overshooting in the pre-main-sequence (PMS) evolution of intermediate-mass stars, by analysing in detail the early evolution towards the main sequence of a  2 M_⊙  stellar model. These effects can be extremely important in the end of the PMS, when the abundances in CNO elements approach the equilibrium in the centre. We provide a possible physical explanation on why a moderate amount of overshooting produces, as the star approaches the zero-age main-sequence, an extra loop in the evolutionary tracks on the Hertzsprung–Russell diagram.
An interesting feature is that there is a very well defined amount of overshooting (for a given stellar mass and chemical composition) beyond which a loop is produced. For smaller amounts of overshooting such a loop does not take place and the evolutionary tracks are similar to those found in the literature. The amount of overshooting needed to produce the loop decreases with stellar mass.
We discuss the underlining physical reasons for the behaviour predicted by the evolution models and argue that it provides a crucial observational test for convective overshooting in the core of intermediate-mass stars.     

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.