Wide-orbit binary evolution
[Papers in prep]
The last two decades of observational advances at long wavelengths and in the infrared have painted an increasingly detailed picture of the star formation process. Stars collapse out of dense cores within filamentary structures in molecular clouds. These filaments are generated by the interplay between supersonic turbulence, self-gravity, and magnetic fields. When these cores collapse to form stars, a majority become binary or higher order multiple systems rather than single stars (Duchene and Kraus 2013). Unlike the initial mass function, which changes little with star forming environment (Offner et al. 2013), the frequency and properties of binary systems do vary from one star cluster to the next (Duchene and Kraus 2013). Thus, binarity appears to be a more sensitive probe of the star formation process. An understanding of binary star formation is required both to explain the properties of the majority of star systems and to gain deeper insight into the properties of the progenitor star forming environments. One of the most fundamental unanswered questions is: why do some filaments and cores produce multiple bound objects, while others produce only single stars?
Initial Conditions matter
Simulations of low-mass star formation must use initial conditions that are consistent with the observed physical properties of star-forming clouds. As first shown by Larson (1981) and confirmed by many studies, the non-thermal motions within the molecular cloud obey a linewidth-size relation, which relates the non-thermal velocity dispersion of the gas and the size of the cloud. We initialize a 512^3 unigrid domain with uniform gas properties and stir the gas according to the prescription of Dubinski (1995). This method drives the gas in Fourier space by using a flat power spectrum on the large scales. The only difference between the simulations is the strength of the magnetic field, values constituting magnetic field energies that exceed, equal, and are less than the turbulent and gravitational energies of the gas.
Once the gas is stirred, gravity is turned on and the simulation commences.
Collapse —> Stars
Gravity condenses the regions of enhanced density resulting from the turbulent stirring phase, and eventually these regions enter a phase of collapse. The collapse proceeds until we initialize protostars as individual “sink particles” that can roam around the domain. The number of stars that form is a function of the magnetic field strength: stronger magnetic fields impede collapse and reduce the star formation efficiency.
However, the multiplicity fractions between the two simulations are a weak function of the magnetic field: all contain relatively the same number of stars in higher-order systems. The magnetic field appears to be a weak trigger of multiplicity, at least for binaries formed from turbulent fragmentation.
Star Formation is a dynamical process
Turbulent fragmentation through gravitational collapse results in a clustered system of particles along filaments. These particles accrete from the gas and also gravitationally interact with each other. Sometimes, these systems get a little too close and disrupt. Here is an example of a massive star (center of frame) that roams near the vicinity of a quadruple star system (upper-right) that we will show below. Interestingly, while the disk system gets destroyed, the inner binaries remain bound to each other.
Wide-orbit binaries migrate into tighter systems
Binary Systems Forming Through Turbulent Fragmentation
Observations of nearby star forming regions reveal two clusters of binary systems: short-period binaries with separations of ~100 AU and wide-orbit binaries with separations around ~10000 AU. The smaller orbits could potentially form from disk fragmentation or from orbital migration. While our simulations lack the resolution to capture disk fragmentation, we reproduce this bimodal distribution through turbulent fragmentation and migration alone. Short-period binaries are formed from turbulent fragmentation and migration and are composed mostly of older protostars, whereas wide-orbit binaries tend to be younger stars that have just condensed from their core’s gas. Wide orbit binaries can also be maintained through dynamical interactions, but many of those interactions tend to disrupt the binary entirely. The age differentiation between wide-orbit and short-period binaries has observational support.