Accretion in stellar clusters and the initial mass function

We present a simple physical mechanism that can account for the observed stellar mass spectrum for masses M ∗≳0.5 M_⊙ . The model depends solely on the competitive accretion that occurs in stellar clusters where each star's accretion rate depends on the local gas density and the square of the accretion radius. In a stellar cluster, there are two different regimes depending on whether the gas or the stars dominate the gravitational potential. When the cluster is dominated by cold gas, the accretion radius is given by a tidal-lobe radius. This occurs as the cluster collapses towards a ρ  ∝  R ⁻² distribution. Accretion in this regime results in a mass spectrum with an asymptotic limit of γ =−3/2 (where Salpeter is γ =−2.35) . Once the stars dominate the potential and are virialized, which occurs first in the cluster core, the accretion radius is the Bondi–Hoyle radius. The resultant mass spectrum has an asymptotic limit of γ =−2 with slightly steeper slopes ( γ ≈−2.5) if the stars are already mass-segregated. Simulations of accretion on to clusters containing 1000 stars show that, as expected, the low-mass stars accumulate the majority of their masses during the gas-dominated phase whereas the high-mass stars accumulate the majority of their masses during the stellar-dominated phase. This results in a mass spectrum with a relatively shallow γ ≈3/2 power law for low-mass stars and a steeper power law for high-mass stars −2.5≲ γ ≤−2 . This competitive accretion model also results in a mass-segregated cluster.     

Non-emission-line young stars of intermediate mass

We present optical spectra of four intermediate-mass candidate young stellar objects that have often been classified as Herbig Ae/Be stars. Typical Herbig Ae/Be emission features are not present in the spectra of these stars. Three of them, HD 36917, HD 36982 and HD 37062, are members of the young Orion nebula cluster (ONC). This association constrains their ages to be ≲1 Myr. The lack of appreciable near-infrared excess in them suggests the absence of hot dust close to the central star. However, they do possess significant amounts of cold and extended dust as revealed by the large excess emission observed at far-infrared wavelengths. The fractional infrared luminosities  ( L _ir/ L _★)  and the dust masses computed from IRAS fluxes are systematically lower than those found for Herbig Ae/Be stars but higher than those for Vega-like stars. These stars may thus represent the youngest examples of the Vega phenomenon known so far. In contrast, the other star in our sample, HD 58647, is more likely to be a classical Be star, as is evident from the low   L _ir/ L _★  , the scarcity of circumstellar dust, the low polarization, the presence of H α emission and near-infrared excess, and the far-infrared spectral energy distribution consistent with free–free emission similar to other well-known classical Be stars.     

Analysing the atolls: X-ray spectral transitions of accreting neutron stars

We systematically analyse all the available X-ray spectra of disc accreting neutron stars (atolls and millisecond pulsars) from the RXTE data base. We show that while all these have similar spectral evolution as a function of mass accretion rate, there are also subtle differences. There are two different types of hard/soft transition, those where the spectrum softens at all energies, leading to a diagonal track on a colour–colour diagram, and those where only the higher energy spectrum softens, giving a vertical track. The luminosity at which the transition occurs is correlated with this spectral behaviour, with the vertical transition at   L / L _Edd∼ 0.02  while the diagonal one is at ∼0.1. Superimposed on this is the well-known hysteresis effect, but we show that classic, large-scale hysteresis occurs only in the outbursting sources, indicating that its origin is in the dramatic rate of change of mass accretion rate during the disc instability. We show that the long-term mass accretion rate correlates with the transition behaviour, and speculate that this is due to the magnetic field being able to emerge from the neutron star surface for low average mass accretion rates. While this is not strong enough to collimate the flow except in the millisecond pulsars, its presence may affect the inner accretion flow by changing the properties of the jet.     

The impact of a supernova explosion in a very massive binary

We consider the effect of a supernova (SN) explosion in a very massive binary that is expected to form in a portion of Population III stars with the mass higher than  100 M_⊙  . In a Population III binary system, a more massive star can result in the formation of a black hole (BH) and a surrounding accretion disc. Such BH accretion could be a significant source of the cosmic reionization in the early Universe. However, a less massive companion star evolves belatedly and eventually undergoes a SN explosion, so that the accretion disc around a BH might be blown off in a lifetime of companion star. In this paper, we explore the dynamical impact of a SN explosion on an accretion disc around a massive BH, and elucidate whether the BH accretion disc is totally demolished or not. For the purpose, we perform three-dimensional hydrodynamic simulations of a very massive binary system, where we assume a BH of  10³ M_⊙  that results from a direct collapse of a very massive star and a companion star of  100 M_⊙  that undergoes a SN explosion. We calculate the remaining mass of a BH accretion disc as a function of time. As a result, it is found that a significant portion of gas disc can survive through three-dimensional geometrical effects even after the SN explosion of a companion star. Even if the SN explosion energy is higher by two orders of magnitude than the binding energy of gas disc, about a half of disc can be left over. The results imply that the Population III BH accretion disc can be a long-lived luminous source, and therefore could be an important ionizing source in the early Universe.     

Numerical simulations of protostellar encounters — I. Star–disc encounters

It appears that most stars are born in clusters, and that at birth most stars have circumstellar discs which are comparable in size to the separations between the stars. Interactions between neighbouring stars and discs are therefore likely to play a key role in determining disc lifetimes, stellar masses, and the separations and eccentricities of binary orbits. Such interactions may also cause fragmentation of the discs, thereby triggering the formation of additional stars.   We have carried out a series of simulations of star–disc interactions using an SPH code which treats self-gravity, hydrodynamic and viscous forces. We find that interactions between discs and stars provide a mechanism for removing energy from, or adding energy to, the orbits of the stars, and for truncating the discs. However, capture during such encounters is unlikely to be an important binary formation mechanism.   A more significant consequence of such encounters is that they can trigger fragmentation of the disc, via tidally and compressionally induced gravitational instabilities, leading to the formation of additional stars and substellar objects. When the disc spins and stellar orbits are randomly oriented, encounters lead to the formation of new companions to the original star in 20 per cent of encounters. If most encounters are prograde and coplanar, as suggested by simulations of dynamically triggered star formation, then new companions are formed in approximately 50 per cent of encounters.     

Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

Infrared spectroscopy of cataclysmic variables – III. Dwarf novae below the period gap and nova-like variables

We present K -band spectra of the short-period dwarf novae YZ Cnc, LY Hya, BK Lyn, T Leo, SW UMa and WZ Sge, the nova-like variables DW UMa, V1315 Aql, RW Tri, VY Scl, UU Aqr and GP Com, and a series of field dwarf stars with spectral types ranging from K2 to M6.
The spectra of the dwarf novae are dominated by emission lines of H  i and He  i . The large velocity and equivalent widths of these lines, in conjunction with the fact that the lines are double-peaked in the highest inclination systems, indicate an accretion disc origin. In the case of YZ Cnc and T Leo, for which we obtained time-resolved data covering a complete orbital cycle, the emission lines show modulations in their equivalent widths that are most probably associated with the bright spot (the region where the gas stream collides with the accretion disc). There are no clear detections of the secondary star in any of the dwarf novae below the period gap, yielding upper limits of 10–30 per cent for the contribution of the secondary star to the observed K -band flux. In conjunction with the K -band magnitudes of the dwarf novae, we use the derived secondary star contributions to calculate lower limits to the distances to these systems.
The spectra of the nova-like variables are dominated by broad, single-peaked emission lines of H  i and He  i – even the eclipsing systems we observed do not show the double-peaked profiles predicted by standard accretion disc theory. With the exception of RW Tri, which exhibits Na  i , Ca  i and ¹²CO absorption features consistent with a M0V secondary contributing 65 per cent of the observed K -band flux, we find no evidence for the secondary star in any of the nova-like variables. The implications of this result are discussed.     

On the effects of the stellar magnetic field on the structure of T Tauri accretion discs

The structure of accretion discs around magnetic T Tauri stars is calculated numerically using a particle hydrodynamical code, in which magnetic interaction is included in the framework of King's diamagnetic blob accretion model. Setting up the calculation so as to simulate the density structure of a quasi-steady disc in the equatorial plane of a T Tauri star, we find that the central star's magnetic field typically produces a central hole in the disc and spreads out the surface density distribution. We argue that this result suggets a promising mechanism for explaining the unusual flatness (IR excess) of T Tauri accretion disc spectra.     

Accretion disc dynamics in extreme mass ratio compact binaries

An analysis is presented of a numerical investigation of the dynamics and geometry of accretion discs in binary systems with mass ratios   q = M ₂/ M ₁ < 0.1  , applicable to ultracompact X-ray binaries, AM CVn stars and very short period cataclysmic variables. The steady-state geometry of the disc in the binary reference frame is found to be quite different from that expected at higher mass ratios. For   q ∼ 0.1  , the disc takes on the usual elliptical shape, with the major axis aligned perpendicular to the line of centres of the two stars. However, at smaller mass ratios the elliptical gaseous orbits in the outer regions of the disc are rotated in the binary plane. The angle of rotation increases with gas temperature, but is found to vary inversely with q . At   q = 0.01  , the major axis of these orbits is aligned almost parallel to the line of centres of the two stars. These effects may be responsible for the similar disc structure inferred from Doppler tomography of the AM CVn star GP Com, which has   q = 0.02  . The steady-state geometry at low mass ratios is not predicted by an inviscid, restricted three-body model of gaseous orbits; it is related to the effects of tidal-viscous truncation of the disc near the Roche lobe boundary. Since the disc geometry can be inferred observationally for some systems, it is proposed that this may offer a useful diagnostic for the determination of mass ratios in ultracompact binaries.     

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.     

The component masses of the cataclysmic variable V347 Puppis

We present time-resolved spectroscopy and photometry of the double-lined eclipsing cataclysmic variable V347 Pup (=LB 1800). There is evidence of irradiation on the inner hemisphere of the secondary star, which we correct for using a model to give a secondary-star radial velocity of   K _R= 198 ± 5 km s⁻¹  . The rotational velocity of the secondary star in V347 Pup is found to be   v sin  i = 131 ± 5 km s⁻¹  and the system inclination is   i = 840 ± 23  . From these parameters we obtain masses of   M ₁= 0.63 ± 0.04 M_⊙  for the white dwarf primary and   M ₂= 0.52 ± 0.06 M_⊙  for the M0.5V secondary star, giving a mass ratio of   q = 0.83 ± 0.05  . On the basis of the component masses, and the spectral type and radius of the secondary star in V347 Pup, we find tentative evidence for an evolved companion. V347 Pup shows many of the characteristics of the SW Sex stars, exhibiting single-peaked emission lines, high-velocity S-wave components and phase-offsets in the radial velocity curve. We find spiral arms in the accretion disc of V347 Pup and measure the disc radius to be close to the maximum allowed in a pressureless disc.     

Global 3‐D solar‐type star‐disc dynamo systems: I. MHD modeling

We have carried out global three‐dimensional magnetohydrodynamic simulations of the star‐disc interaction region around a young solar‐type star. The magnetic field is generated and maintained by dynamos in the star as well as in the disc. The developing mass flows possess non‐periodic time‐variable azimuthal structure and are controlled by the nonaxisymmetric magnetic fields. Since the stellar field drives a strong stellar wind, accretion is anti‐correlated with the stellar field strength and disc matter is spiraling onto the star at low latitudes, both contrary to the generally assumed accretion picture. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The accretion of brown dwarfs and planets by giant stars – II. Solar-mass stars on the red giant branch

This paper extends our previous study of planet/brown dwarf accretion by giant stars to solar-mass stars located on the red giant branch. The model assumes that the planet is dissipated at the bottom of the convective envelope of the giant star. The evolution of the giant is then followed in detail. We analyse the effects of different accretion rates and different initial conditions. The computations indicate that the accretion process is accompanied by a substantial expansion of the star, and, in the case of high accretion rates, hot bottom burning can be activated. The possible observational signatures that accompany the engulfing of a planet are also extensively investigated. They include the ejection of a shell and a subsequent phase of IR emission, an increase in the ⁷Li surface abundance and a potential stellar metallicity enrichment, spin-up of the star because of the deposition of orbital angular momentum, the possible generation of magnetic fields and the related X-ray activity caused by the development of shear at the base of the convective envelope, and the effects on the morphology of the horizontal branch in globular clusters. We propose that the IR excess and high Li abundance observed in 4–8 per cent of the G and K giants originate from the accretion of a giant planet, a brown dwarf or a very low-mass star.     

Boundary layer on the surface of a neutron star

In an attempt to model the accretion on to a neutron star in low-mass X-ray binaries, we present 2D hydrodynamical models of the gas flow in close vicinity of the stellar surface. First, we consider a gas pressure-dominated case, assuming that the star is non-rotating. For the stellar mass we take   M _star= 1.4 × 10⁻² M_⊙  and for the gas temperature   T = 5 × 10⁶ K  . Our results are qualitatively different in the case of a realistic neutron star mass and a realistic gas temperature of T ≃ 10⁸ K, when the radiation pressure dominates. We show that to get the stationary solution in a latter case, the star most probably has to rotate with the considerable velocity.     

X-rays from cusps of compact remnants near galactic centres

Compact remnants – stellar mass black holes and neutron stars formed in the inner few parsec of galactic centres are predicted to sink into the central parsec due to dynamical friction on low-mass stars, forming a high concentration cusp. Same physical region may also contain very high-density molecular clouds and accretion discs that are needed to fuel supermassive black hole (SMBH) activity. Here we estimate gas capture rates on to the cusp of stellar remnants, and the resulting X-ray luminosity, as a function of the accretion disc mass. At low disc masses, most compact objects are too dim to be observable, whereas in the high disc case most of them are accreting at their Eddington rates. We find that for low accretion disc masses, compact remnant cusps may be more luminous than the central SMBHs. This 'diffuse' emission may be of importance for local moderately bright active galactic nuclei (AGNs), especially low-luminosity AGNs. We also briefly discuss how this expected emission can be used to put constraints on the black hole cusp near our Galactic Centre.     

Influence of the r-mode instability on hypercritically accreting neutron stars

We have investigated the influence of the r-mode instability on hypercritically accreting neutron stars in close binary systems during their common envelope phases, based on the scenario proposed by Brown et al. On the one hand, neutron stars are heated by the accreted matter at the stellar surface, but on the other hand they are also cooled down by the neutrino radiation. At the same time, the accreted matter transports its angular momentum and mass to the star. We have studied the evolution of the stellar mass, temperature and rotational frequency.
The gravitational-wave-driven instability of the r-mode oscillation strongly suppresses spinning up of the star, the final rotational frequency of which is well below the mass-shedding limit, in fact typically as low as 10 per cent of that of the mass-shedding state. On a very short time-scale the rotational frequency tends to approach a certain constant value and saturates there, as long as the amount of accreted mass does not exceed a certain limit to collapse to a black hole. This implies that a similar mechanism of gravitational radiation to that in the so-called 'Wagoner star' may work in this process. The star is spun up by accretion until the angular momentum loss by gravitational radiation balances the accretion torque. The time-integrated dimensionless strain of the radiated gravitational wave may be large enough to be detectable by gravitational wave detectors such as LIGO II.     

The stellar composition of the star formation region CMa R1 – II. Spectroscopic and photometric observations of nine young stars

We present new high- and low-resolution spectroscopic and photometric data of nine members of the young association CMa R1. All the stars have circumstellar dust at some distance, as could be expected from their association with reflection nebulosity. Four stars (HD 52721, HD 53367, LkH α  220 and LkH α  218) show H α emission and we argue that they are Herbig Be stars with discs. Our photometric and spectroscopic observations of these stars reveal new characteristics of their variability. We present first interpretations of the variability of HD 52721, HD 53367 and the two LkH α stars in terms of a partially eclipsing binary, a magnetic activity cycle and circumstellar dust variations, respectively. The remaining five stars show no clear indications of H α emission in their spectra, although their spectral types and ages are comparable with those of HD 52721 and HD 53367. This indicates that the presence of a disc around a star in CMa R1 may depend on the environment of the star. In particular we find that all H α emission stars are located at or outside the arc-shaped border of the H  ii region, which suggests that the stars inside the arc have lost their discs through evaporation by UV photons from nearby O stars, or from the nearby (<25 pc) supernova, about 1 Myr ago.     

General relativistic spectra of accretion discs around rapidly rotating neutron stars: effect of light bending

We present computed spectra, as seen by a distant observer, from the accretion disc around a rapidly rotating neutron star. Our calculations are carried out in a fully general relativistic framework, with an exact treatment of rotation. We take into account the Doppler shift, gravitational redshift and light-bending effects in order to compute the observed spectrum. We find that light bending significantly modifies the high-energy part of the spectrum. Computed spectra for slowly rotating neutron stars are also presented. These results would be important for modelling the observed X-ray spectra of low-mass X-ray binaries containing fast-spinning neutron stars.     

Periodicities in the high-mass X-ray binary system RX J0146.9+6121/LS I+61°235

The high-mass X-ray binary RX J0146.9+6121, with optical counterpart LS I+61°235 (V831 Cas), is an intriguing system on the outskirts of the open cluster NGC 663. It contains the slowest Be type X-ray pulsar known with a pulse period of around 1400 s and, primarily from the study of variation in the emission line profile of Hα, it is known to have a Be decretion disc with a one-armed density wave period of approximately 1240 d. Here we present the results of an extensive photometric campaign, supplemented with optical spectroscopy, aimed at measuring short time-scale periodicities. We find three significant periodicities in the photometric data at, in order of statistical significance, 0.34, 0.67 and 0.10 d. We give arguments to support the interpretation that the 0.34 and 0.10 d periods could be due to stellar oscillations of the B-type primary star and that the 0.67 d period is the spin period of the Be star with a spin axis inclination of  23⁺¹⁰₋₈  degrees. We measured a systemic velocity of  −37.0 ± 4.3 km s⁻¹  confirming that LS I+61°235 has a high probability of membership in the young cluster NGC 663 from which the system's age can be estimated as 20–25 Myr. From archival RXTE All Sky Monitor (ASM) data we further find 'super' X-ray outbursts roughly every 450 d. If these super outbursts are caused by the alignment of the compact star with the one-armed decretion disc enhancement, then the orbital period is approximately 330 d.