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.     

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.     

Binary formation in stellar clusters

We consider how the tidal potential of a stellar cluster or a dense molecular cloud affects the fragmentation of gravitationally unstable molecular cloud cores. We find that molecular cloud cores which would collapse to form a single star in the absence of tidal shear, can be forced to fragment if they are subjected to tides. This may enhance the frequency of binaries in star-forming regions such as Ophiuchus and the frequency of binaries with separations ≲100 au in the Orion Trapezium Cluster. We also find that clouds which collapse to form binary systems in the absence of a tidal potential will form bound binary systems if exposed to weak tidal shear. However, if the tidal shear is sufficiently strong, even though the cloud still collapses to form two fragments, the fragments are pulled apart while they are forming by the tidal shear and two single stars are formed. This sets an upper limit for the separation of binaries that form near dense molecular clouds or in stellar clusters.     

The brown dwarf desert as a consequence of orbital migration

We show that the dearth of brown dwarfs in short-period orbits around Solar-mass stars – the brown dwarf desert – can be understood as a consequence of inward migration within an evolving protoplanetary disc. Brown dwarf secondaries forming at the same time as the primary star have masses which are comparable to the initial mass of the protoplanetary disc. Subsequent disc evolution leads to inward migration, and destruction of the brown dwarf, via merger with the star. This is in contrast with massive planets, which avoid this fate by forming at a later epoch when the disc is close to being dispersed. Within this model, a brown dwarf desert arises because the mass at the hydrogen-burning limit is coincidentally comparable to the initial disc mass for a Solar mass star. Brown dwarfs should be found in close binaries around very low mass stars, around other brown dwarfs, and around Solar-type stars during the earliest phases of star formation.     

Close binary systems among very low‐mass stars and brown dwarfs

Using Monte Carlo simulations and published radial velocity surveys we have constrained the frequency and separation (a ) distribution of very low‐mass star (VLM) and brown dwarf (BD) binary systems.We find that simple Gaussian extensions of the observed wide binary distribution, with a peak at 4AU and 0.6 < σ _{log(a /AU)} < 1.0, correctly reproduce the observed number of close binary systems, implying a close (a < 2.6 AU) binary frequency of 17–30% and overall frequency of 32–45%. N‐body models of the dynamical decay of unstable protostellar multiple systems are excluded with high confidence because they do not produce enough close binary VLMs/BDs. The large number of close binaries and high overall binary frequency are also completely inconsistent with published smoothed particle hydrodynamical modelling and argue against a dynamical origin for VLMs/BDs. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Enhanced Activity in Close T Tauri Binaries

The number of confirmed and suspected close T Tauri binaries (period days) is increasing. We discuss some systems with enhanced emission line activity and periodic line profile changes. Non-axisymmetric flows of plasma in the region between the circumbinary disk and the stars can be generated through the influence of the secondary component. Such enhanced activity is found around binaries with eccentric as well as circular orbits. We discuss our observations of the T Tauri stars RW Aurigae A and RU Lupi, which may host very close brown dwarf companions. Model simulations indicate that non-axisymmetric flows are generated around close binaries with circumbinary disks, also in systems with circular orbits.     

Spectral atlas of dwarf novae in outburst

Up to now, only a very small number of dwarf novae have been studied during their outburst state (∼30 per cent in the Northern hemisphere). In this paper we present the first comprehensive atlas of outburst spectra of dwarf novae. We study possible correlations between the emission and absorption lines seen in the spectra and some fundamental parameters of the binaries. We find that out of the 48 spectra presented, 12 systems apart from IP Peg show strong He  ii in emission: SS Aur, HL CMa, TU Crt, EM Cyg, SS Cyg, EX Dra, U Gem, HX Peg, GK Per, KT Per, V893 Sco, IY UMa, and seven others less prominently: FO And, V542 Cyg, B  i Ori, TY Psc, VZ Pyx, ER UMa and SS UMi. We conclude that these systems are good targets for finding spiral structure in their accretion discs during outburst if the models of Smak and Ogilvie are correct. This is confirmed by the fact that hints of spiral asymmetries have already been found in the discs of SS Cyg, EX Dra and U Gem.     

On the evolution of giant protoplanets forming in circumbinary discs

We present the results of hydrodynamic simulations of Jovian mass protoplanets that form in circumbinary discs. The simulations follow the orbital evolution of the binary plus protoplanet system acting under their mutual gravitational forces, and forces exerted by the viscous circumbinary disc. The evolution involves the clearing of the inner circumbinary disc initially, so that the binary plus protoplanet system orbits within a low density cavity. Continued interaction between disc and protoplanet causes inward migration of the planet towards the inner binary. Subsequent evolution can take three distinct paths: (i) the protoplanet enters the 4 : 1 mean motion resonance with the binary, but is gravitationally scattered through a close encounter with the secondary star; (ii) the protoplanet enters the 4 : 1 mean motion resonance, the resonance breaks, and the planet remains in a stable orbit just outside the resonance; (iii) when the binary has initial eccentricity   e _bin≥ 0.2  , the disc becomes eccentric, leading to a stalling of the planet migration, and the formation of a stable circumbinary planet.
These results have implications for a number of issues in the study of extrasolar planets. The ejection of protoplanets in close binary systems provides a source of 'free-floating planets', which have been discovered recently. The formation of a large, tidally truncated cavity may provide an observational signature of circumbinary planets during formation. The existence of protoplanets orbiting stably just outside a mean motion resonance (4 : 1) in the simulations indicate that such sites may harbour planets in binary star systems, and these could potentially be observed. Finally, the formation of stable circumbinary planets in eccentric binary systems indicates that circumbinary planets may not be uncommon.     

Thermal stability and nova cycles in permanent superhump systems

Archival data on permanent superhump systems are compiled to test the thermal stability of their accretion discs. We find that their discs are almost certainly thermally stable as expected. This result confirms Osaki's suggestion that permanent superhump systems form a new subclass of cataclysmic variables (CVs), with relatively short orbital periods and high mass-transfer rates. We note that if the high accretion rates estimated in permanent superhump systems represent their mean secular values, then their mass-transfer rates cannot be explained by gravitational radiation, therefore, either magnetic braking should be extrapolated to systems below the period gap or they must have mass-transfer cycles. Alternatively, a new mechanism that removes angular momentum from CVs below the gap should be invoked.
We suggest applying the nova cycle scenarios offered for systems above the period gap to the short orbital period CVs. Permanent superhumps have been observed in the two non-magnetic ex-novae with binary periods below the gap. Their post-nova magnitudes are brighter than their pre-outburst values. In one case (V1974 Cyg) it has been demonstrated that the pre-nova should have been a regular SU UMa system. Thus, it is the first nova whose accretion disc was observed to change its thermal stability. If the superhumps in this system indicate persistent high mass-transfer rates rather than a temporary change induced by irradiation from the hot post-nova white dwarf, it is the first direct evidence for mass-transfer cycles in CVs. The proposed cycles are driven by the nova eruption.     

Are superhumps good measures of the mass ratio for AM CVn systems?

We extend recent work that included the effect of pressure forces to derive the precession rate of eccentric accretion discs in cataclysmic variables to the case of double degenerate systems. We find that the logical scaling of the pressure force in such systems results in predictions of unrealistically high primary masses. Using the prototype AM CVn as a calibrator for the magnitude of the effect, we find that there is no scaling that applies consistently to all the systems in the class. We discuss the reasons for the lack of a superhump period to mass ratio relationship analogous to that known for SU UMa systems and suggest that this is because these secondaries do not have a single valued mass–radius relationship. We highlight the unreliability of mass ratios derived by applying the SU UMa expression to the AM CVn binaries.     

Widely separated binary systems of very low‐mass stars

In this paper we review some recent detections of wide binary brown dwarf systems and discuss them in the context of the multiplicity properties of very low‐mass stars and brown dwarfs. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The standstill luminosity in Z Cam systems

We consider accretion discs in close binary systems. We show that the heating of a disc at the impact point of the accretion stream contributes significantly to the local energy budget at its outer edge. As a result, the thermal balance relation between local accretion rate and surface density (the 'S-curve') changes; the critical mass transfer rate above which no dwarf nova outbursts occur can be up to 40 per cent smaller than without impact heating. Standstills in Z Cam systems thus occur at smaller mass transfer rates than otherwise expected, and are fainter than the peak luminosity during the dwarf nova phase as a result.     

Precessing warped accretion discs in X-ray binaries

We study the radiation-driven warping of accretion discs in the context of X-ray binaries. The latest evolutionary equations are adopted, which extend the classical alpha theory to time-dependent thin discs with non-linear warps. We also develop accurate, analytical expressions for the tidal torque and the radiation torque, including self-shadowing.
We investigate the possible non-linear dynamics of the system within the framework of bifurcation theory. First, we re-examine the stability of an initially flat disc to the Pringle instability. Then we compute directly the branches of non-linear solutions representing steadily precessing discs. Finally, we determine the stability of the non-linear solutions. Each problem involves only ordinary differential equations, allowing a rapid, accurate and well-resolved solution.
We find that radiation-driven warping is probably not a common occurrence in low-mass X-ray binaries. We also find that stable, steadily precessing discs exist for a narrow range of parameters close to the stability limit. This could explain why so few systems show clear, repeatable 'superorbital' variations. The best examples of such systems, Her X-1, SS 433 and LMC X-4, all lie close to the stability limit for a reasonable choice of parameters. Systems far from the stability limit, including Cyg X-2, Cen X-3 and SMC X-1, probably experience quasi-periodic or chaotic variability as first noticed recently by Wijers and Pringle. We show that radiation-driven warping provides a coherent and persuasive framework but that it does not provide a generic explanation for the long-term variabilities in all X-ray binaries.     

The potential for Earth-mass planet formation around brown dwarfs

Recent observations point to the presence of structured dust grains in the discs surrounding young brown dwarfs, thus implying that the first stages of planet formation take place also in the substellar regime. Here, we investigate the potential for planet formation around brown dwarfs and very low-mass stars according to the sequential core accretion model of planet formation. We find that, for a brown dwarf mass 0.05 M_⊙, our models predict a maximum planetary mass of  ∼5   M_⊕  , orbiting with semimajor axis ∼ 1 au. However, we note that the predictions for the mass–semimajor axis distribution are strongly dependent upon the models chosen for the disc surface density profiles and the assumed distribution of disc masses. In particular, if brown dwarf disc masses are of the order of a few Jupiter masses, Earth-mass planets might be relatively frequent, while if typical disc masses are only a fraction of Jupiter mass, we predict that planet formation would be extremely rare in the substellar regime. As the observational constraints on disc profiles, mass dependencies and their distributions are poor in the brown dwarf regime, we advise caution in validating theoretical models only on stars similar to the Sun and emphasize the need for observational data on planetary systems around a wide range of stellar masses. We also find that, unlike the situation around solar-like stars, Type II migration is totally absent from the planet formation process around brown dwarfs, suggesting that any future observations of planets around brown dwarfs would provide a direct measure of the role of other types of migration.     

Evidence for pole switching in the magnetic cataclysmic variable BY Camelopardalis

An analysis of X-ray and optical light curves of the magnetic cataclysmic variable (MCV) BY Cam is presented. This system is one of three MCVs in which the spin period of the white dwarf and the binary orbital period differ by ∼1 per cent. As such these 'BY Cam' stars are important objects with which to probe the field structure of the magnetic white dwarf and ultimately the nature of synchronization of AM Her binaries. We confirm asynchronous rotation of the magnetic white dwarf with respect to the binary. We find evidence that the accretion stream accretes directly on to the white dwarf as in AM Her systems, but further, the stream impacts on to different magnetic poles over the course of the beat period. We present evidence that the optical and hard X-ray light curves modulate in phase, but together they are out of phase with the soft X-ray light curve. We confirm the spin down of the white dwarf which is expected to lead to the synchronization of the spin and orbital periods of BY Cam.     

Numerical simulations of protostellar encounters — II. Coplanar disc–disc encounters

It is expected that an average protostar will undergo at least one impulsive interaction with a neighbouring protostar whilst a large fraction of its mass is still in a massive, extended disc. Such interactions must have a significant impact upon the evolution of the protostars and their discs.   We have carried out a series of simulations of coplanar encounters between two stars, each possessing a massive circumstellar disc, using an SPH code that models gravitational, hydrodynamic and viscous forces. We find that during a coplanar encounter, disc material is swept up into a shock layer between the two interacting stars, and the layer then fragments to produce new protostellar condensations. The truncated remains of the discs may subsequently fragment; and the outer regions of the discs may be thrown off to form circumbinary disc-like structures around the stars. Thus coplanar disc–disc encounters lead efficiently to the formation of multiple star systems and small- N clusters, including substellar objects.     

Stellar-mass black hole binaries as ultraluminous X-ray sources

Ultraluminous X-ray sources (ULXs) with   L _x > 10³⁹ erg s⁻¹  have been discovered in great numbers in external galaxies with ROSAT , Chandra and XMM-Newton . The central question regarding this important class of sources is whether they represent an extension in the luminosity function of binary X-ray sources containing neutron stars and stellar-mass black holes (BHs), or a new class of objects, e.g. systems containing intermediate-mass BHs  (100–1000 M_⊙)  . We have carried out a theoretical study to test whether a large fraction of the ULXs, especially those in galaxies with recent star formation activity, can be explained with binary systems containing stellar-mass BHs. To this end, we have applied a unique set of binary evolution models for BH X-ray binaries, coupled to a binary population synthesis code, to model the ULXs observed in external galaxies. We find that for donor stars with initial masses  ≳10 M_⊙  the mass transfer driven by the normal nuclear evolution of the donor star is sufficient to potentially power most ULXs. This is the case during core hydrogen burning and, to an even more pronounced degree, while the donor star ascends the giant branch, although the latter phases last only ∼5 per cent of the main-sequence phase. We show that with only a modest violation of the Eddington limit, e.g. a factor of ∼10, both the numbers and properties of the majority of the ULXs can be reproduced. One of our conclusions is that if stellar-mass BH binaries account for a significant fraction of ULXs in star-forming galaxies, then the rate of formation of such systems is  ∼3 × 10⁻⁷ yr⁻¹  normalized to a core-collapse supernova rate of 0.01 yr⁻¹.     

Very low mass stars in binaries: a theoretical look

We have undertaken a series of hydrodynamic + N ‐body simulations in order to explore the binary properties of young stars. We find that multiple stars are a natural outcome of collapsing turbulent flows, with a high incidence of N > 2 multiples, specially among the higher mass objects. We find a positive correlation of multiplicity with primary mass and a companion frequency that decreases with age, during the first few Myr after formation. Binary brown dwarfs are rarely formed, in conflict with observations. Brown dwarfs as companions are predominantly found orbiting binaries or triples at large separations. The paucity of ultra low mass and low mass ratio binaries has been investigated further, and we tentatively conclude that their formation is intricately related to an appropriate selection of initial conditions and an accurate modelling of disc accretion and evolution. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Two- and three-dimensional numerical simulations of accretion discs in a close binary system

We perform 2D and 3D numerical simulations of an accretion disc in a close binary system using the simplified flux vector splitting (SFS) finite volume method. In our calculations, the gas is assumed to be ideal with γ =1.01, 1.05, 1.1 and 1.2 . The mass ratio of the mass-losing star to the mass-accreting star is unity. Our results show that spiral shocks are formed on the accretion disc in all cases. In 2D calculations we find that the smaller γ is, the more tightly the spiral winds. We observe this trend in 3D calculations as well in a somewhat weaker sense. Mach numbers in our discs are less than 10. These values are lower than the values in observed accretion discs in close binary systems.
Recently, Steeghs, Harlaftis & Horne found the first convincing evidence for spiral structure in the accretion disc of the eclipsing dwarf nova binary IP Pegasi, using the technique known as Doppler tomography. Although the Mach numbers in present calculations are rather low, we may claim that the spiral structure that we discovered in earlier numerical simulations is now found observationally.     

Photometry of four binary subdwarf B stars and the nature of their unseen companion stars

We present light curves of four binary subdwarf B stars (sdB), Ton 245, Feige 11, PG 1432+159 and PG 1017−086. We also present new spectroscopic data for PG 1017−086 from which we derive its orbital period,   P =0.073 d  , and the mass function,   f _m=0.0010±0.0002 M_⊙.  This is the shortest period for an sdB binary measured to date. The values of P and f _m for the other sdB binaries have been published elsewhere. We are able to exclude the possibility that the unseen companion stars to Ton 245, Feige 11 and PG 1432+159 are main-sequence stars or subgiant stars from the absence of a sinusoidal signal, which would be caused by the irradiation of such a companion star, i.e. they show no reflection effect. The unseen companion stars in these binaries are likely to be white dwarf stars. In contrast, the reflection effect in PG 1017−086 is clearly seen. The lack of eclipses in this binary combined with other data suggests that the companion is a low-mass M-dwarf or, perhaps, a brown dwarf.