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.     

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)     

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 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.     

On the migration of a system of protoplanets

The evolution of a system consisting of a protoplanetary disc with two embedded Jupiter-sized planets is studied numerically. The disc is assumed to be flat and non-self-gravitating; this is modelled by the planar (two-dimensional) Navier–Stokes equations. The mutual gravitational interaction of the planets and the star, and the gravitational torques of the disc acting on the planets and the central star are included. The planets have an initial mass of one Jupiter mass M _Jup each, and the radial distances from the star are one and two semimajor axes of Jupiter, respectively.
During the evolution a joint wide annular gap is created by the planets. Both planets increase their mass owing to accretion of gas from the disc: after about 2500 orbital periods of the inner planet it has reached a mass of 2.3  M _Jup, while the outer planet has reached a mass of 3.2  M _Jup. The net gravitational torques exerted by the disc on the planets result in an inward migration of the outer planet on time-scales comparable to the viscous evolution time of the disc. The semimajor axis of the inner planet remains constant as there is very little gas left in its vicinity to induce any migration. When the distance of close approach eventually becomes smaller than the mutual Hill radius, the eccentricities increase strongly and the system may become unstable.
If disc depletion occurs rapidly enough before the planets come too close to each other, a stable system similar to our own Solar system may remain. Otherwise the orbits may become unstable and produce systems like υ And.     

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.     

Predictions for the frequency and orbital radii of massive extrasolar planets

We investigate the migration of massive extrasolar planets caused by gravitational interaction with a viscous protoplanetary disc. We show that a model in which planets form at 5 au at a constant rate, before migrating, leads to a predicted distribution of planets that is a steeply rising function of log( a ), where a is the orbital radius. Between 1 and 3 au, the expected number of planets per logarithmic interval in a roughly doubles. We demonstrate that, once selection effects are accounted for, this is consistent with current data, and then extrapolate the observed planet fraction to masses and radii that are inaccessible to current observations. In total, approximately 15 per cent of stars targeted by existing radial velocity searches are predicted to possess planets with masses  0.3< M _p sin( i )<10 M _J  and radii  0.1< a <5 au  . A third of these planets (around 5 per cent of the target stars) lie at the radii most amenable to detection via microlensing. A further  5–10  per cent of stars could have planets at radii of  5< a <8 au  that have migrated outwards. We discuss the probability of forming a system (akin to the Solar system) in which significant radial migration of the most massive planet does not occur. Approximately  10–15  per cent of systems with a surviving massive planet are estimated to fall into this class. Finally, we note that a smaller fraction of low-mass planets than high-mass planets is expected to survive without being consumed by the star. The initial mass function for planets is thus predicted to rise more steeply towards small masses than the observed mass function.     

A source of high-velocity white dwarfs

We investigate whether the recently observed population of high-velocity white dwarfs can be derived from a population of binaries residing initially within the thin disc of the Galaxy. In particular, we consider binaries where the primary is sufficiently massive to explode as a Type II supernova. A large fraction of such binaries are broken up when the primary then explodes as a supernova, owing to the combined effects of the mass loss from the primary and the kick received by the neutron star on its formation. For binaries where the primary evolves to fill its Roche lobe, mass transfer from the primary leads to the onset of a common envelope phase during which the secondary and the core of the primary spiral together as the envelope is ejected. Such binaries are the progenitors of X-ray binaries if they are not broken up when the primary explodes. For those systems that are broken up, a large number of the secondaries receive kick velocities ∼100–200 km s⁻¹ and subsequently evolve into white dwarfs. We compute trajectories within the Galactic potential for this population of stars and relate the birth rate of these stars over the entire Galaxy to those seen locally with high velocities relative to the local standard of rest (LSR) . We show that for a reasonable set of assumptions concerning the Galactic supernova rate and the binary population, our model produces a local number density of high-velocity white dwarfs compatible with that inferred from observations. We therefore propose that a population of white dwarfs originating in the thin disc may make a significant contribution to the observed population of high-velocity white dwarfs.     

The migration and growth of protoplanets in protostellar discs

We investigate the gravitational interaction of a Jovian-mass protoplanet with a gaseous disc with aspect ratio and kinematic viscosity expected for the protoplanetary disc from which it formed. Different disc surface density distributions are investigated. We focus on the tidal interaction with the disc with the consequent gap formation and orbital migration of the protoplanet. Non-linear two-dimensional hydrodynamic simulations are employed using three independent numerical codes.
A principal result is that the direction of the orbital migration is always inwards and such that the protoplanet reaches the central star in a near-circular orbit after a characteristic viscous time‐scale of ∼10⁴ initial orbital periods. This is found to be independent of whether the protoplanet is allowed to accrete mass or not. Inward migration is helped by the disappearance of the inner disc, and therefore the positive torque it would exert, because of accretion on to the central star. Maximally accreting protoplanets reach about 4 Jovian masses on reaching the neighbourhood of the central star. Our results indicate that a realistic upper limit for the masses of closely orbiting giant planets is ∼5 Jupiter masses, if they originate in protoplanetary discs similar to the minimum-mass solar nebula. This is because of the reduced accretion rates obtained for planets of increasing mass.
Assuming that some process such as termination of the inner disc through a magnetospheric cavity stops the migration, the range of masses estimated for a number of close orbiting giant planets as well as their inward orbital migration can be accounted for by consideration of disc–protoplanet interactions during the late stages of giant planet formation.     

A General Relativistic Calculation of Boundary Layer and Disc Luminosity for Accreting Non-magnetic Neutron Stars in Rapid Rotation

We calculate the disc and boundary layer luminosities for accreting rapidly rotating neutron stars with low magnetic fields in a fully general relativistic manner. Rotation increases the disc luminosity and decreases the boundary layer luminosity. A rapid rotation of the neutron star substantially modifies these quantities as compared with the static limit. For a neutron star rotating close to the centrifugal mass shed limit, the total luminosity has contribution only from the extended disc. For such maximal rotation rates, we find that well before the maximum stable gravitational mass configuration is reached, there exists a limiting central density, for which particles in the innermost stable orbit will be more tightly bound than those at the surface of the neutron star. We also calculate the angular velocity profiles of particles in Keplerian orbits around the rapidly rotating neutron star. The results are illustrated for a representative set of equation of state models of neutron star matter.     

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.     

GD 552: a cataclysmic variable with a brown dwarf companion?

GD 552 is a high proper motion star with the strong, double-peaked emission lines characteristic of the dwarf nova class of cataclysmic variable (CV) star, and yet no outburst has been detected during the past 12 yr of monitoring. We present spectroscopy taken with the aim of detecting emission from the mass donor in this system. We fail to do so at a level which allows us to rule out the presence of a near-main-sequence star donor. Given GD 552's orbital period of 103 min, this suggests that it is either a system that has evolved through the ∼80-minute orbital period minimum of CV stars and now has a brown dwarf mass donor, or that has formed with a brown dwarf donor in the first place. This model explains the low observed orbital velocity of the white dwarf and GD 552's low luminosity. It is also consistent with the absence of outbursts from the system.     

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.     

Hotspots and a clumpy disc: variability of brown dwarfs and stars in the young σ Ori cluster

The properties of accretion discs around stars and brown dwarfs in the σ Ori cluster (age 3 Myr) are studied based on near-infrared (IR) time series photometry supported by mid-IR spectral energy distributions (SEDs). We monitor ∼30 young low-mass sources over eight nights in the J and K band using the duPont telescope at Las Campanas. We find three objects showing variability with J -band amplitudes  ≥0.5 mag  ; five additional objects exhibit low-level variations. All three highly variable sources have been previously identified as highly variable; thus, we establish the long-term nature of their flux changes. The light curves contain periodic components with time-scales of  ∼0.5–8 d  , but have additional irregular variations superimposed – the characteristic behaviour for classical T Tauri stars. Based on the colour variability, we conclude that hotspots are the dominant cause of variations in two objects (#19 and #33), including one likely brown dwarf, with spot temperatures in the range of 6000–7000 K. For the third one (#2), a brown dwarf or very low-mass star, inhomogeneities at the inner edge of the disc are the likely origin of variability. Based on mid-IR data from Spitzer , we confirm that the three highly variable sources are surrounded by circum-(sub)-stellar discs. They show typical SEDs for T Tauri-like objects. Using SED models, we infer an enhanced scaleheight in the disc for the object #2, which favours the detection of disc inhomogeneities in light curves and is thus consistent with the information from variability. In the σ Ori cluster, about every fifth accreting low-mass object shows persistent high-level photometric variability. We demonstrate that estimates for fundamental parameters in such objects can be significantly improved by determining the extent and origin of the variations.     

Launching of conical winds and axial jets from the disc–magnetosphere boundary: axisymmetric and 3D simulations

We investigate the launching of outflows from the disc–magnetosphere boundary of slowly and rapidly rotating magnetized stars using axisymmetric and exploratory 3D magnetohydrodynamic simulations. We find long-lasting outflows in the following cases. (1) In the case of slowly rotating stars , a new type of outflow, a conical wind , is found and studied in simulations. The conical winds appear in cases where the magnetic flux of the star is bunched up by the disc into an X-type configuration. The winds have the shape of a thin conical shell with a half-opening angle  θ∼ 30°–40°  . About 10–30 per cent of the disc matter flows from the inner disc into the conical winds. The conical winds may be responsible for episodic as well as long-lasting outflows in different types of stars. There is also a low-density, higher velocity component (a jet) in the region inside the conical wind. (2) In the case of rapidly rotating stars (the 'propeller regime'), a two-component outflow is observed. One component is similar to the conical winds. A significant fraction of the disc matter may be ejected into the winds. The second component is a high-velocity, low-density magnetically dominated axial jet where matter flows along the opened polar field lines of the star. The jet has a mass flux of about 10 per cent of that of the conical wind, but its energy flux (dominantly magnetic) can be larger than the energy flux of the conical wind. The jet's angular momentum flux (also dominantly magnetic) causes the star to spin down rapidly. Propeller-driven outflows may be responsible for the jets in protostars and for their rapid spin-down. The jet is collimated by the magnetic force while the conical winds are only weakly collimated in the simulation region. Exploratory 3D simulations show that conical winds are axisymmetric about the rotational axis (of the star and the disc), even when the dipole field of the star is significantly misaligned.     

Numerical simulations of type III planetary migration – II. Inward migration of massive planets

We present a numerical study of rapid, so-called type III migration for Jupiter-sized planets embedded in a protoplanetary disc. We limit ourselves to the case of inward migration, and study in detail its evolution and physics, concentrating on the structure of the corotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using the results of global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. The initial conditions are chosen to satisfy the condition for rapid inward migration. We find that type III migration can be divided into two regimes, fast and slow. The structure of the co-orbital region, mass accumulation rate and migration behaviour differ between these two regimes. All our simulations show a transition from the fast to the slow regime, ending type III migration well before reaching the star. The stopping radius is found to be larger for more massive planets and less massive discs. A sharp density drop is also found to be an efficient stopping mechanism. In the fast migration regime the migration rate and induced eccentricity are lower for less massive discs, but almost do not depend on planet mass. Eccentricity is damped on the migration time-scale.     

The Monitor project: searching for occultations in young open clusters

The Monitor project is a photometric monitoring survey of nine young (1–200 Myr) clusters in the solar neighbourhood to search for eclipses by very low mass stars and brown dwarfs and for planetary transits in the light curves of cluster members. It began in the autumn of 2004 and uses several 2- to 4-m telescopes worldwide. We aim to calibrate the relation between age, mass, radius and where possible luminosity, from the K dwarf to the planet regime, in an age range where constraints on evolutionary models are currently very scarce. Any detection of an exoplanet in one of our youngest targets (≲10 Myr) would also provide important constraints on planet formation and migration time-scales and their relation to protoplanetary disc lifetimes. Finally, we will use the light curves of cluster members to study rotation and flaring in low-mass pre-main-sequence stars.
The present paper details the motivation, science goals and observing strategy of the survey. We present a method to estimate the sensitivity and number of detections expected in each cluster, using a simple semi-analytic approach which takes into account the characteristics of the cluster and photometric observations, using (tunable) best-guess assumptions for the incidence and parameter distribution of putative companions, and we incorporate the limits imposed by radial velocity follow-up from medium and large telescopes. We use these calculations to show that the survey as a whole can be expected to detect over 100 young low and very low mass eclipsing binaries, and ∼3 transiting planets with radial velocity signatures detectable with currently available facilities.     

Mass transfer in tidally unstable compact binaries

The 2001 outburst of WZ Sagittae has shown the most compelling evidence yet for an enhancement of the mass-transfer rate from the donor star during a dwarf nova outburst in the form of hotspot brightening. I show that, even in this extreme case, the brightening can be attributed to tidal heating near the interaction point of an accretion stream with the expanding edge of an eccentric accretion disc, with no need at all for an increase in the mass-transfer rate. Furthermore, I confirm previous suggestions that an increase in mass-transfer rate through the stream damps any eccentricity in an accretion disc and suppresses the appearance of superhumps, in contradiction to observations. Tidal heating is expected to be most significant in systems with small mass ratios. It follows that systems like WZ Sagittae – which has a tiny mass ratio – are those most likely to show a brightening in the hotspot region.     

S 235 B explained: an accreting Herbig Be star surrounded by reflection nebulosity

The intent of this study is to determine the nature of the star and associated nebulosity S 235 B, which are located in a region of active star formation still heavily obscured by the parent molecular cloud. Low-resolution  ( R = 400)  long-slit spectra of the star and nebulosity, and medium-  ( R = 1800)  and high-resolution  ( R = 60 000)  spectra of the central star are presented along with the results of Fabry–Perot interferometric imaging of the entire region. Based on the long-slit and Fabry–Perot observations, the nebulosity appears to be entirely reflective in nature, with the stellar component S 235 B★ providing most of the illuminating flux. The stellar source itself is classified here as a B1V star, with emission-line profiles indicative of an accretion disc. S 235 B★ thus belongs to the relatively rare class of early-type Hebrig Be stars. Based on the intensity of the reflected component, it is concluded that the accretion disc must be viewed nearly edge-on. Estimates of the accretion rate of S 235 B★ from the width of the Hα profile at 10 per cent of maximum intensity, a method which has been used lately for T Tauri stars and Brown Dwarfs, appear to be inconsistent with the mass outflow rate and accretion rate implied from previous infrared observations by Felli et al., suggesting this empirical law does not extend to higher masses.     

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.