Predicting the properties of binary stellar systems: the evolution of accreting protobinary systems

We investigate the formation of binary stellar systems. We consider a model where a 'seed' protobinary system forms, via fragmentation, within a collapsing molecular cloud core and evolves to its final mass by accreting material from an infalling gaseous envelope. This accretion alters the mass ratio and orbit of the binary, and is largely responsible for forming the circumstellar and/or circumbinary discs.
Given this model for binary formation, we predict the properties of binary systems and how they depend on the initial conditions within the molecular cloud core. We predict that there should be a continuous trend such that closer binaries are more likely to have equal-mass components and are more likely to have circumbinary discs than wider systems. Comparing our results with observations, we find that the observed mass-ratio distributions of binaries and the frequency of circumbinary discs as a function of separation are most easily reproduced if the progenitor molecular cloud cores have radial density profiles between uniform and 1/ r (e.g., Gaussian) with near-uniform rotation. This is in good agreement with the observed properties of pre-stellar cores. Conversely, we find that the observed properties of binaries cannot be reproduced if the cloud cores are in solid-body rotation and have initial density profiles which are strongly centrally condensed. Finally, in agreement with the radial-velocity searches for extrasolar planets, we find that it is very difficult to form a brown dwarf companion to a solar-type star with a separation ≲10 au, but that the frequency of brown dwarf companions should increase with larger separations or lower mass primaries.     

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.     

Binary formation with different metallicities: dependence on initial conditions

The fragmentation process in collapsing clouds with various metallicities is studied using three-dimensional nested-grid hydrodynamics. Initial clouds are specified by three parameters: cloud metallicity, initial rotation energy and initial cloud shape. For different combinations of these parameters, we calculate 480 models in total and study cloud evolution, fragmentation conditions, orbital separation and binary frequency. For the cloud to fragment during collapse, the initial angular momentum must be higher than a threshold value, which decreases with decreasing metallicity. Although the exact fragmentation conditions depend also on the initial cloud shape, this dependence is only modest. Our results indicate a higher binary frequency in lower metallicity gas. In particular, with the same median rotation parameter as in the solar neighbourhood, a majority of stars are born as members of binary/multiple systems for  <10⁻⁴ Z_⊙  . With initial mass  <0.1 M_⊙  , if fragments are ejected in embryo from the host clouds by multibody interaction, they evolve to substellar-mass objects. This provides a formation channel for low-mass stars in zero- or low-metallicity environments.     

How well do starlab and nbody4 compare? I. Simple models

N -body simulations are widely used to simulate the dynamical evolution of a variety of systems, among them star clusters. Much of our understanding of their evolution rests on the results of such direct N -body simulations. They provide insight in the structural evolution of star clusters, as well as into the occurrence of stellar exotica. Although the major pure N -body codes starlab/kira and nbody4 are widely used for a range of applications, there is no thorough comparison study yet.
Here, we thoroughly compare basic quantities as derived from simulations performed either with starlab/kira or nbody4 .
We construct a large number of star cluster models for various stellar mass function settings (but without stellar/binary evolution, primordial binaries, external tidal fields, etc.), evolve them in parallel with starlab/kira and nbody4 , analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison, we develop a bootstrap algorithm for functional dependencies.
We find an overall excellent agreement between the codes, both for the clusters' structural and energy parameters as well as for the properties of the dynamically created binaries. However, we identify small differences, like in the energy conservation before core collapse and the energies of escaping stars, which deserve further studies.
Our results reassure the comparability and the possibility to combine results from these two major N -body codes, at least for the purely dynamical models (i.e. without stellar/binary evolution) we performed. Further detailed comparison studies for more complex systems, e.g. including stellar/binary evolution, are required.     

Fragmentation in molecular clouds and its connection to the IMF

We present an analysis of star-forming gas cores in a smooth particle hydrodynamics simulation of a giant molecular cloud. We identify cores using their deep potential wells. This yields a smoother distribution with clearer boundaries than density. Additionally, this gives an indication of future collapse, as bound potential cores (p-cores) represent the earliest stages of fragmentation in molecular clouds. We find that the mass function of the p-cores resembles the stellar initial mass function and the observed clump mass function, although p-core masses  (∼0.7 M_⊙)  are smaller than typical density clumps. The bound p-cores are generally subsonic, have internal substructure and are only quasi-spherical. We see no evidence of massive bound cores supported by turbulence. We trace the evolution of the p-cores forward in time, and investigate the connection between the original p-core mass and the stellar mass that formed from it. We find that there is a poor correlation, with considerable scatter suggesting accretion on to the core is dependent on more factors than just the initial core mass. During the accretion process the p-cores accrete from beyond the region first bound, highlighting the importance of the core environment to its subsequent evolution.     

Theoretical models for wide binary evolution

Summary Wide binary evolution under the influence of stars, giant molecular clouds, and the galactic tidal field is reviewed. Recent results show that binary lifetimes ata≳0.1 pc are determined by an interplay between relatively gentle evolution due to stellar encounters and catastrophic collisions with cloud subclumps. It is emphasized that binary semimajor axis distributions may not be understood from lifetimes alone. To determine the theoretical distribution binary ‘birthrate functions’ must be convolved with probability distributions for the evolution of an ensemble of binaries from given initial separations. Simple models show no sharp breaks below the tidal cut-off imposed by the Galaxy.     

An analysis of v sin (i) correlations in early-type binaries

We use information on the   v sin ( i )  values of early-type binaries in order to search for correlations which may constrain the relative orientation of the stellar spin axes in binary systems. We find correlations in the case of close binaries which suggest that tidal synchronization is effective for binaries whose separation exceeds the stellar radius by more than an order of magnitude, in line with the theoretical predictions of Goldreich & Nicholson and the previous observational analysis of Giuricin et al. In the case of wide binaries, the   v sin ( i )  values are not well correlated, which requires that the magnitude of the spin speeds is not tightly correlated. Under this assumption, we then find that the data provide no significant constraints on the degree of alignment of spin axes. The data are therefore compatible with scenarios (such as disc fragmentation or capture) which differ widely in the expected degree of spin alignment.     

ECHA J0843.3–7905: Discovery of an ''old'' classical T Tauri star in the η Chamaeleontis cluster

Cold, dense clouds of gas have been proposed to explain the dark matter in Galactic haloes, and have also been invoked in the Galactic disc as an explanation for the excess faint submillimetre sources detected by SCUBA. Even if their dust-to-gas ratio is only a small percentage of that in conventional gas clouds, these dense systems would be opaque to visible radiation. We examine the possibility that the data sets of microlensing experiments searching for massive compact halo objects can also be used to search for occultation signatures by such clouds. We compute the rate and time-scale distribution of stellar transits by clouds in the Galactic disc and halo. We find that, for cloud parameters typically advocated by theoretical models, thousands of transit events should already exist within microlensing survey data sets. We examine the seasonal modulation in the rate caused by the Earth's orbital motion and find it provides an excellent probe of whether detected clouds are of disc or halo origin.     

Dense and diffuse gas in dynamically active clouds

We investigate the chemical and observational implications of repetitive transient dense core formation in molecular clouds. We allow a transient density fluctuation to form and disperse over a period of 1 Myr, tracing its chemical evolution. We then allow the same gas immediately to undergo further such formation and dispersion cycles. The chemistry of the dense gas in subsequent cycles is similar to that of the first, and a limit cycle is reached quickly (2–3 cycles). Enhancement of hydrocarbon abundances during a specific period of evolution is the strongest indicator of previous dynamical history. The molecular content of the diffuse background gas in the molecular cloud is expected to be strongly enhanced by the core formation and dispersion process. Such enhancement may remain for as long as 0.5 Myr. The frequency of repetitive core formation should strongly determine the level of background molecular enhancement.
We also convolve the emission from a synthesized dark cloud, comprised of ensembles of transient dense cores. We find that the dynamical history of the gas, and therefore the chemical state of the diffuse intercore medium, may be determined if a sufficient sample of cores is present in an ensemble. Molecular ratios of key hydrocarbons with SO and SO₂ are crucial to this distinction. Only surveys with great enough angular resolution to resolve individual cores, or very small groupings, are expected to show evidence of repetitive dynamical processing. The existence of non-equilibrium chemistry in the diffuse background may have implications for the initial conditions used in chemical models. Observed variations in the chemistries of diffuse and translucent regions may be explained by lines of sight which intersect a number of molecular cloud cores in various stages of evolution.     

On the formation of massive stars

We present a model for the formation of massive ( M ≳10 M⊙) stars through accretion-induced collisions in the cores of embedded dense stellar clusters. This model circumvents the problem of accreting on to a star whose luminosity is sufficient to reverse the infall of gas. Instead, the central core of the cluster accretes from the surrounding gas, thereby decreasing its radius until collisions between individual components become sufficient. These components are, in general, intermediate-mass stars that have formed through accretion on to low-mass protostars. Once a sufficiently massive star has formed to expel the remaining gas, the cluster expands in accordance with this loss of mass, halting further collisions. This process implies a critical stellar density for the formation of massive stars, and a high rate of binaries formed by tidal capture.     

Close encounters in young stellar clusters: implications for planetary systems in the solar neighbourhood

The stars that populate the solar neighbourhood were formed in stellar clusters. Through N -body simulations of these clusters, we measure the rate of close encounters between stars. By monitoring the interaction histories of each star, we investigate the singleton fraction in the solar neighbourhood. A singleton is a star which formed as a single star, has never experienced any close encounters with other stars or binaries, or undergone an exchange encounter with a binary. We find that, of the stars which formed as single stars, a significant fraction is not singletons once the clusters have dispersed. If some of these stars had planetary systems, with properties similar to those of the Solar System, the planets' orbits may have been perturbed by the effects of close encounters with other stars or the effects of a companion star within a binary. Such perturbations can lead to strong planet–planet interactions which eject several planets, leaving the remaining planets on eccentric orbits. Some of the single stars exchange into binaries. Most of these binaries are broken up via subsequent interactions within the cluster, but some remain intact beyond the lifetime of the cluster. The properties of these binaries are similar to those of the observed binary systems containing extrasolar planets. Thus, dynamical processes in young stellar clusters will alter significantly any population of Solar System-like planetary systems. In addition, beginning with a population of planetary systems exactly resembling the Solar System around single stars, dynamical encounters in young stellar clusters may produce at least some of the extrasolar planetary systems observed in the solar neighbourhood.     

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.     

On the binarity of Herbig Ae/Be stars

We present high-resolution spectro-astrometry of a sample of 28 Herbig Ae/Be and three F-type pre-main-sequence stars. The spectro-astrometry, which is essentially the study of unresolved features in long-slit spectra, is shown from both empirical and simulated data to be capable of detecting binary companions that are fainter by up to 6 mag at separations larger than ∼0.1 arcsec. The nine targets that were previously known to be binary are all detected. In addition, we report the discovery of six new binaries and present five further possible binaries. The resulting binary fraction is 68 ± 11 per cent. This overall binary fraction is the largest reported for any observed sample of Herbig Ae/Be stars, presumably because of the exquisite sensitivity of spectro-astrometry for detecting binary systems. The data hint that the binary frequency of the Herbig Be stars is larger than that of the Herbig Ae stars. The Appendix presents model simulations to assess the capabilities of spectro-astrometry and reinforces the empirical findings. Most spectro-astrometric signatures in this sample of Herbig Ae/Be stars can be explained by the presence of a binary system. Two objects, HD 87643 and Z CMa, display evidence for asymmetric outflows. Finally, the position angles of the binary systems have been compared with available orientations of the circumprimary disc and these appear to be coplanar. The alignment between the circumprimary discs and the binary systems strongly suggests that the formation of binaries with intermediate-mass primaries is due to fragmentation as the alternative, stellar capture, does not naturally predict aligned discs. The alignment extends to the most massive B-type stars in our sample. This leads us to conclude that formation mechanisms that do result in massive stars, but predict random angles between the binaries and the circumprimary discs, such as stellar collisions, are also ruled out for the same reason.     

The approach to collapse of molecular clouds

The dense molecular cloud cores that form stars, like other self-gravitating objects, undergo bulk oscillations. Just at the point of gravitational instability, their fundamental oscillation mode has zero frequency. We study, using perturbation theory, the evolution of a spherical cloud that possesses such a frozen mode. We find that the cloud undergoes a prolonged epoch of subsonic, accelerating contraction. This slow contraction occurs whether the cloud is initially inflated or compressed by the oscillation. The subsonic motion described here could underlie the spectral infall signature observed in many starless dense cores.     

The proper-motion signal of unresolved binaries in the Hipparcos catalogue

We present an investigation of the differences between quasi-instantaneous stellar proper motions from the Hipparcos catalogue and long-term proper motions determined by combining Hipparcos and the Astrographic Catalogue. Our study is based on a sample of about 12000 stars of visual magnitude from 7 to 10 in two declination zones on the northern and equatorial sky. The distribution of the proper-motion differences shows an excess of large deviations. This is caused by the influence of orbital motion of unresolved binary systems. The proper-motion deviations provide statistical evidence for 360 astrometric binaries in the investigated zones, corresponding to about 2400 such binaries in the entire Hipparcos catalogue, in addition to those already known. In order to check whether the observed deviations are compatible with standard assumptions on the basic parameters of binary stars, we model the impact of orbital motion on the observed proper motions in a Monte Carlo simulation. We show that the simulation yields an acceptable approximation of the observations, if a binary frequency between 70% and 100% is assumed, i.e.if most of the stars in the sample are assumed to have a companion. Thus Hipparcos astrometric binaries confirm that the frequency of non-single stars among field stars is very high. We also investigate the influence of the mass function for the secondary component on the result of the simulation. A constant mass function and mass functions with moderate increase towards low masses lead to results, which are compatible with the observed proper-motion effects. A high preponderance of very-low-mass or substellar companions as produced, for example, by a M^—1 power law is not in agreement with the frequency of proper-motion deviations in our sample of stars.     

Comments on Gravoturbulent Star Formation

Understanding the star formation process is central to much of modern astrophysics. Stellar birth is intimately linked to the dynamical behavior of the parental gas cloud. Gravoturbulent fragmentation determines where and when protostellar cores form, and how they contract and grow in mass via accretion from the surrounding cloud material to build up stars. Supersonic turbulence can provide support against gravitational collapse on global scales, whereas at the same time it produces localized density enhancements that allow for collapse on small scales. The efficiency and timescale of stellar birth in Galactic molecular clouds strongly depend on the properties of the interstellar turbulent velocity field, with slow, inefficient, isolated star formation being a hallmark of turbulent support, and fast, efficient, clustered star formation occurring in its absence.     

The evolution of binary fractions in globular clusters

We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code ( startrack ) and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing three- and four-body interactions ( fewbody ). We find that the combination of stellar evolution and dynamical interactions (binary–single and binary–binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster such as 47 Tuc, assuming an initial binary fraction of 100 per cent, is only ∼ 5–10 per cent. We show that this is in good agreement with recent Hubble Space Telescope observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modelling of globular clusters, suggesting that 'realistic models' should incorporate much larger initial binary fractions than has usually been the case in the past.     

Star formation and extinct radioactivities

A supernova trigger for star formation, such as was discussed by Cameron and Truran a few years ago, lacks observational confirmation, and the nucleosynthesis arguments for suggesting its relevance to the formation of the solar system no longer seem compelling. Observations of star formation currently point toward nearly isolated individual events taking place in the interiors of dense molecular clouds, in which cores are formed in the clouds and collapse to form stars. If the formation of the solar system proceeded in this manner, then a viable theory of cloud core formation must provide circumstances in which it is not uncommon for the cloud cores to be formed with a significant supply of fresh short-lived radioactivities. A review is made of the evidence for the existence of now extinct radioactivities in primitive solar system material and an examination is made of the implications for the early stages of formation of the Sun and solar system. The characteristics of possible disturbances in dense molecular clouds which can initiate the formation of cloud cores is discussed,and in particular those disturbances which can produce fresh radioactivities are considered. A red giant star of roughly one solar mass on the asymptotic giant branch appears to have been the best candidate to account for the short-lived extinct radioactivities in the early solar system. Star formation itself can generate disturbances through the T Tauri stellar winds and bipolar outflows, and it appears likely that several such disturbances in a chain of star-forming events were needed before the radioactive material was diluted enough to be consistent with the abundances seen in early solar system material.     

A turbulent MHD model for molecular clouds and a new method of accretion on to star-forming cores

We describe the results of a sequence of simulations of gravitational collapse in a turbulent magnetized region. The parameters are chosen to be representative of molecular cloud material. We find that several protostellar cores and filamentary structures of higher than average density form. The filaments inter connect the high-density cores. Furthermore, the magnetic field strengths are found to correlate positively with the density, in agreement with recent observations. We make synthetic channel maps of the simulations, and show that material accreting on to the cores is channelled along the magnetized filamentary structures. This is compared with recent observations of S106, and shown to be consistent with these data. We postulate that this mechanism of accretion along filaments may provide a means for molecular cloud cores to grow to the point where they become gravitationally unstable and collapse to form stars.     

Monte Carlo simulations of star clusters – II. Tidally limited, multimass systems with stellar evolution

A revision of Stodółkiewicz's Monte Carlo code is used to simulate evolution of large star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. A survey of the evolution of N -body systems influenced by the tidal field of a parent galaxy and by stellar evolution is presented. The process of energy generation is realized by means of appropriately modified versions of Spitzer's and Mikkola's formulae for the interaction cross-section between binaries and field stars and binaries themselves. The results presented are in good agreement with theoretical expectations and the results of other methods (Fokker–Planck, Monte Carlo and N -body). The initial rapid mass loss, resulting from stellar evolution of the most massive stars, causes expansion of the whole cluster and eventually leads to the disruption of less bound systems ( W ₀=3). Models with larger W ₀ survive this phase of evolution and then undergo core collapse and subsequent post-collapse expansion, like isolated models. The expansion phase is eventually reversed when tidal limitation becomes important. The results presented are the first major step in the direction of simulating evolution of real globular clusters by means of the Monte Carlo method.