Protostellar Disks Research Articles

Revealing the dynamics of Class 0 protostellar discs with ALMA

We present synthetic ALMA observations of Keplerian, protostellar discs in the Class 0 stage studying the emission of molecular tracers like $^{13}$CO, C$^{18}$O, HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, and H$_2$CO. We model the emission of discs around low- and intermediate-mass protostars. We show that under optimal observing conditions ALMA is able to detect the discs already in the earliest stage of protostellar evolution, although the emission is often concentrated to the innermost 50 AU. Therefore, a resolution of a few 0.1" might be too low to detect Keplerian discs around Class 0 objects. We also demonstrate that under optimal conditions for edge-on discs Keplerian rotation signatures are recognisable, from which protostellar masses can be inferred. For this we here introduce a new approach, which allows us to determine protostellar masses with higher fidelity than before. Furthermore, we show that it is possible to reveal Keplerian rotation even for strongly inclined discs and that ALMA should be able to detect possible signs of fragmentation in face-on discs. In order to give some guidance for future ALMA observations, we investigate the influence of varying observing conditions and source distances. We show that it is possible to probe Keplerian rotation in inclined discs with an observing time of 2 h and a resolution of 0.1", even in the case of moderate weather conditions. Furthermore, we demonstrate that under optimal conditions, Keplerian discs around intermediate-mass protostars should be detectable up to kpc-distances.

On the mechanism of self gravitating Rossby interfacial waves in proto-stellar accretion discs

The dynamical response of edge waves under the influence of self-gravity is examined in an idealised two-dimensional model of a proto-stellar disc, characterised in steady state as a rotating vertically infinite cylinder of fluid with constant density except for a single density interface at some radius . The fluid in basic state is prescribed to rotate with a Keplerian profile modified by some additional azimuthal sheared flow. A linear analysis shows that there are two azimuthally propagating edge waves, kin to the familiar Rossby waves and surface gravity waves in terrestrial studies, which move opposite to one another with respect to the local basic state rotation rate at the interface. Instability only occurs if the radial pressure gradient is opposite to that of the density jump (unstably stratified) where self-gravity acts as a wave stabiliser irrespective of the stratification of the system. The propagation properties of the waves are discussed in detail in the language of vorticity edge waves. The roles of both Boussinesq and non-Boussinesq effects upon the stability and propagation of these waves with and without the inclusion of self-gravity are then quantified. The dynamics involved with self-gravity non-Boussinesq effect is shown to be a source of vorticity production where there is a jump in the basic state density In addition, self-gravity also alters the dynamics via the radial main pressure gradient, which is a Boussinesq effect. Further applications of these mechanical insights are presented in the conclusion including the ways in which multiple density jumps or gaps may or may not be stable.

Supermassive star formation via episodic accretion: protostellar disc instability and radiative feedback efficiency

The formation of SMSs is a potential pathway to seed SMBHs in the early universe. A critical issue for forming SMSs is stellar UV feedback, which may limit the stellar mass growth via accretion. In this paper we study the evolution of an accreting SMS and its UV emissivity under conditions of realistic variable accretion from a self-gravitating circumstellar disc. First we conduct a 2D hydrodynamical simulation to follow the long-term protostellar accretion until the stellar mass exceeds $10^4~M_\odot$. The disc fragments due to gravitational instability, creating a number of small clumps that rapidly migrate inward to fall onto the star. The resulting accretion history is thus highly time-dependent: short episodic accretion bursts are followed by longer, relative quiescent phases. We show that the circumstellar disc for the so-called direct collapse model is more unstable and generates greater variability over shorter timescales than normal Pop III cases. We conduct a post-process stellar evolution calculation using the obtained accretion history. Our results show that, regardless of the strong variability of the accretion rates, the stellar radius monotonically increases with almost constant effective temperature at $T_{\rm eff} \simeq 5000$ K as the stellar mass increases. The resulting UV feedback is too weak to hinder mass accretion due to the low flux of stellar UV photons, thus verifying our implicit assumption of no stellar feedback during the hydrodynamic simulations. The insensitivity of stellar evolution to variable accretion is attributed to the fact that typical timescales of variability, $\lesssim 10^3$ years, are too short to affect the stellar structure. We argue that this evolution will continue until the SMS eventually collapses to produce a massive black hole by the general relativistic instability after the stellar mass reaches $\gtrsim 10^5~M_\odot$.

MAGNETIC FIELDS IN EARLY PROTOSTELLAR DISK FORMATION

ABSTRACT We consider formation of accretion disks from a realistically turbulent molecular gas using 3D MHD simulations. In particular, we analyze the effect of the fast turbulent reconnection described by the Lazarian & Vishniac model for the removal of magnetic flux from a disk. With our numerical simulations we demonstrate how the fast reconnection enables protostellar disk formation resolving the so-called “magnetic braking catastrophe.” In particular, we provide a detailed study of the dynamics of a 0.5 M ⊙ protostar and the formation of its disk for up to several thousands years. We measure the evolution of the mass, angular momentum, magnetic field, and turbulence around the star. We consider effects of two processes that strongly affect the magnetic transfer of angular momentum, both of which are based on turbulent reconnection: the first, “reconnection diffusion,” removes the magnetic flux from the disk; the other involves the change of the magnetic field's topology, but does not change the absolute value of the magnetic flux through the disk. We demonstrate that for the first mechanism, turbulence causes a magnetic flux transport outward from the inner disk to the ambient medium, thus decreasing the coupling of the disk to the ambient material. A similar effect is achieved through the change of the magnetic field's topology from a split monopole configuration to a dipole configuration. We explore how both mechanisms prevent the catastrophic loss of disk angular momentum and compare both above turbulent reconnection mechanisms with alternative mechanisms from the literature.

Self-gravitational instability of dusty plasma with dissipative effects

The influence of resistivity and viscosity effects on Jeans self-gravitational instability of dusty plasma is studied. The governing equations are constructed using magnetohydrodynamic model for dusty plasma. The general dispersion relation and Jeans criteria for instability are obtained employing the plane wave solutions on linearized perturbation equations. The Jeans criterion is found to be unaffected due to resistivity and viscosity parameter. The dispersion relation is analyzed for various directions of propagation with respect to magnetic field. The numerical calculations have been performed to show the effects of resistivity and viscosity on the growth rate of Jeans instability. It is observed that the presence of resistivity has destabilizing influence while viscosity has stabilizing effects on the growth rate of the Jeans instability. The analytical results are verified by the numerical analysis using both real and normalized astrophysical parameters. The present work is applicable to various space and astrophysical plasma systems, particularly for molecular clouds, protostellar disks, interstellar & circumstellar clouds and the ionosphere of the earth, etc.

EFFECTS OF DYNAMICAL EVOLUTION OF GIANT PLANETS ON THE DELIVERY OF ATMOPHILE ELEMENTS DURING TERRESTRIAL PLANET FORMATION

ABSTRACT Recent observations started revealing the compositions of protostellar disks and planets beyond the solar system. In this paper, we explore how the compositions of terrestrial planets are affected by the dynamical evolution of giant planets. We estimate the initial compositions of the building blocks of these rocky planets by using a simple condensation model, and numerically study the compositions of planets formed in a few different formation models of the solar system. We find that the abundances of refractory and moderately volatile elements are nearly independent of formation models, and that all the models could reproduce the abundances of these elements of the Earth. The abundances of atmophile elements, on the other hand, depend on the scattering rate of icy planetesimals into the inner disk, as well as the mixing rate of the inner planetesimal disk. For the classical formation model, neither of these mechanisms are efficient and the accretion of atmophile elements during the final assembly of terrestrial planets appears to be difficult. For the Grand Tack model, both of these mechanisms are efficient, which leads to a relatively uniform accretion of atmophile elements in the inner disk. It is also possible to have a “hybrid” scenario where the mixing is not very efficient but the scattering is efficient. The abundances of atmophile elements in this case increase with orbital radii. Such a scenario may occur in some of the extrasolar planetary systems, which are not accompanied by giant planets or those without strong perturbations from giants. We also confirm that the Grand Tack scenario leads to the distribution of asteroid analogues where rocky planetesimals tend to exist interior to icy ones, and show that their overall compositions are consistent with S-type and C-type chondrites, respectively.

THE VLA NASCENT DISK AND MULTIPLICITY SURVEY: FIRST LOOK AT RESOLVED CANDIDATE DISKS AROUND CLASS 0 AND I PROTOSTARS IN THE PERSEUS MOLECULAR CLOUD

ABSTRACT We present the first dust emission results toward a sample of seven protostellar disk candidates around Class 0 and I sources in the Perseus molecular cloud from the VLA Nascent Disk and Multiplicity (VANDAM) survey with ∼0.″05 or 12 AU resolution. To examine the surface brightness profiles of these sources, we fit the Ka-band 8 mm dust-continuum data in the u, v-plane to a simple, parametrized model based on the Shakura–Sunyaev disk model. The candidate disks are well-fit by a model with a disk-shaped profile and have masses consistent with known Class 0 and I disks. The inner-disk surface densities of the VANDAM candidate disks have shallower density profiles compared to disks around more evolved Class II systems. The best-fit model radii of the seven early-result candidate disks are R c > 10 AU; at 8 mm, the radii reflect lower limits on the disk size since dust continuum emission is tied to grain size and large grains radially drift inwards. These relatively large disks, if confirmed kinematically, are inconsistent with theoretical models where the disk size is limited by strong magnetic braking to <10 AU at early times.

Magnetic fields in primordial accretion disks

Magnetic fields are considered as a vital ingredient of contemporary star formation, and may have been important during the formation of the first stars in the presence of an efficient amplification mechanism. Initial seed fields are provided via plasma fluctuations, and are subsequently amplified by the small-scale dynamo, leading to a strong tangled magnetic field. Here we explore how the magnetic field provided by the small-scale dynamo is further amplified via the $\alpha-\Omega$ dynamo in a protostellar disk and assess its implications. For this purpose, we consider two characteristic cases, a typical Pop.~III star with $10$~M$_\odot$ and an accretion rate of $10^{-3}$~M$_\odot$~yr$^{-1}$, and a supermassive star with $10^5$~M$_\odot$ and an accretion rate of $10^{-1}$~M$_\odot$~yr$^{-1}$. For the $10$~M$_\odot$ Pop.~III star, we find that coherent magnetic fields can be produced on scales of at least $100$~AU, which are sufficient to drive a jet with a luminosity of $100$~L$_\odot$ and a mass outflow rate of $10^{-3.7}$~M$_\odot$~yr$^{-1}$. For the supermassive star, the dynamical timescales in its environment are even shorter, implying smaller orbital timescales and an efficient magnetization out to at least $1000$~AU. The jet luminosity corresponds to $\sim10^{6.0}$~L$_\odot$, and a mass outflow rate of $10^{-2.1}$~M$_\odot$~yr$^{-1}$. We expect that the feedback from the supermassive star can have a relevant impact on its host galaxy.

The Evolution of Self-Gravitating Accretion Discs

AbstractIt is quite likely that self-gravity will play an important role in the evolution of accretion discs, in particular, those around young stars, and those around supermassive black holes. We summarise, here, our current understanding of the evolution of such discs, focussing more on discs in young stellar system, than on discs in active galactic nuclei. We consider the conditions under which such discs may fragment to form bound objects, and when they might, instead, be expected to settle into a quasi-steady, self-regulated state. We also discuss how this understanding may depend on the mass of the disc relative to the mass of the central object, and how it might depend on the presence of external irradiation. Additionally, we consider whether or not fragmentation might be stochastic, where we might expect it to occur in an actual protostellar disc, and if there is any evidence for fragmentation actually playing a role in the formation of planetary-mass bodies. Although there are still a number of outstanding issue, such as the convergence of simulations of self-gravitating discs, whether or not there is more than one mode of fragmentation, and quite what role self-gravitating discs may play in the planet-formation process, our general understanding of these systems seems quite robust.

A Theoretical Perspective on the Formation and Fragmentation of Protostellar Discs

AbstractWe discuss the factors influencing the formation and gravitational fragmentation of protostellar discs. We start with a review of how observations of prestellar cores can be analysed statistically to yield plausible initial conditions for simulations of their subsequent collapse. Simulations based on these initial conditions show that, despite the low levels of turbulence in prestellar cores, they deliver primary protostars and associated discs which are routinely subject to stochastic impulsive perturbations; consequently misalignment of the spins and orbits of protostars are common. Also, the simulations produce protostars that collectively have a mass function and binary statistics matching those observed in nearby star-formation regions, but only if a significant fraction of the turbulent energy in the core is solenoidal, and accretion onto the primary protostar is episodic with a duty cycle $\stackrel{>}{\sim }3\,000\,{\rm yr}$. Under this circumstance, a core typically spawns between 4 and 5 protostars, with high efficiency, and the lower mass protostars are mainly formed by disc fragmentation. The requirement that a proto-fragment in a disc lose thermal energy on a dynamical timescale dictates that there is a sweet spot for disc fragmentation at radii $70\,{\rm AU}\stackrel{<}{\sim }R\stackrel{<}{\sim }100\,{\rm AU}$ and temperatures $10\,{\rm K}\stackrel{<}{\sim }T\stackrel{<}{\sim }20\,{\rm K}$, and this might explain the brown dwarf desert.

Formation, disruption and energy output of Population III X-ray binaries

The first astrophysical objects shaped the cosmic environment by reionizing and heating the intergalactic medium (IGM). In particular, X-rays are very efficient at heating the IGM before it became completely ionized, an effect that can be measured through the 21 cm line of neutral hydrogen. High-mass X-ray binaries (HMXBs), known to be prolific X-ray sources in star-forming galaxies at lower redshifts, are prime candidates for driving the thermal evolution of the IGM at redshifts $z > 20$. Despite their importance, the formation efficiency of HMXBs from the first stellar populations is not well understood---as such, their collective X-ray emission and the subsequent imprint on the 21 cm signature are usually evaluated using free parameters. Using $N$-body simulations, we estimate the rate of HMXB formation via mutual gravitational interactions of nascent, small groups of the first stars (Pop III stars). We run two sets of calculations: one in which stars form in small groups of five in nearly Keplerian initial orbits, and another in which two such groups collide (expected from mergers of host protogalaxies). We find that HMXBs form at a rate of one per $~10^{4}~{\rm M}_{\odot}$ in newly born stars, and that they emit with a power of $\sim 10^{41} ~{\rm erg}~{\rm s} ^{-1}$ in the $2-10$ keV band per solar mass per year of star formation. This value is a factor $\sim 10^{2}$ larger than what is observed in star forming galaxies at lower redshifts; the X-ray production from early HMXBs would have been even more copious, if they also formed $in$ $situ$ or via migration in protostellar disks. Combining our results with earlier studies suggests that early HMXBs were highly effective at heating the IGM and leaving a strong 21 cm signature. We discuss broader implications of our results, such as the rate of long GRBs from PopIII stars and the direct collapse black hole formation.

Disc fragmentation rarely forms planetary-mass objects

It is now reasonably clear that disc fragmentation can only operate in the outer parts of protostellar discs ($r > 50$ au). It is also expected that any object that forms via disc fragmentation will have an initial mass greater than that of Jupiter. However, whether or not such a process actually operates, or can play a significant role in the formation of planetary-mass objects, is still unclear. We do have a few examples of directly imaged objects that may have formed in this way, but we have yet to constrain how often disc fragmentation may actually form such objects. What we want to consider here is whether or not we can constrain the likely population of planetary-mass objects formed via disc fragmentation by considering how a population of objects at large radii ($a > 50$) au - if they do exist - would evolve under perturbations from more distant stellar companions. We find that there is a specific region of parameter space to which such objects would be scattered and show that the known exoplanets in that region have properties more consistent with that of the bulk exoplanet population, than with having been formed via disc fragmentation at large radii. Along with the scarcity of directly-imaged objects at large radii, our results provide a similar, but independent, constraint on the frequency of objects formed via disc fragmentation.

SELF-GRAVITATIONAL FORCE CALCULATION OF SECOND-ORDER ACCURACY FOR INFINITESIMALLY THIN GASEOUS DISKS IN POLAR COORDINATES

Investigating the evolution of disk galaxies and the dynamics of proto-stellar disks can involve the use of both a hydrodynamical and a Poisson solver. These systems are usually approximated as infinitesimally thin disks using two- dimensional Cartesian or polar coordinates. In Cartesian coordinates, the calcu- lations of the hydrodynamics and self-gravitational forces are relatively straight- forward for attaining second order accuracy. However, in polar coordinates, a second order calculation of self-gravitational forces is required for matching the second order accuracy of hydrodynamical schemes. We present a direct algorithm for calculating self-gravitational forces with second order accuracy without artifi- cial boundary conditions. The Poisson integral in polar coordinates is expressed in a convolution form and the corresponding numerical complexity is nearly lin- ear using a fast Fourier transform. Examples with analytic solutions are used to verify that the truncated error of this algorithm is of second order. The kernel integral around the singularity is applied to modify the particle method. The use of a softening length is avoided and the accuracy of the particle method is significantly improved.

Spiral-driven accretion in protoplanetary discs

We numerically investigate the dynamics of a 2D non-magnetised protoplanetary disc surrounded by an inflow coming from an external envelope. We find that the accretion shock between the disc and the inflow is unstable, leading to the generation of large-amplitude spiral density waves. These spiral waves propagate over long distances, down to radii at least ten times smaller than the accretion shock radius. We measure spiral-driven outward angular momentum transport with 1e-4 < alpha < 1e-2 for an inflow accretion rate Mout>1e-8 Msun/yr. We conclude that the interaction of the disc with its envelope leads to long-lived spiral density waves and radial angular momentum transport with rates that cannot be neglected in young non-magnetised protostellar discs.

THE MIGRATION OF GAS GIANT PLANETS IN GRAVITATIONALLY UNSTABLE DISKS

Planets form in the discs of gas and dust that surround young stars. It is not known whether gas giant planets on wide orbits form the same way as Jupiter or by fragmentation of gravitationally unstable discs. Here we show that a giant planet, which has formed in the outer regions of a protostellar disc, initially migrates fast towards the central star (migration timescale ~10,000 yr) while accreting gas from the disc. However, in contrast with previous studies, we find that the planet eventually opens up a gap in the disc and the migration is essentially halted. At the same time, accretion-powered radiative feedback from the planet, significantly limits its mass growth, keeping it within the planetary mass regime (i.e. below the deuterium burning limit) at least for the initial stages of disc evolution. Giant planets may therefore be able to survive on wide orbits despite their initial fast inward migration, shaping the environment in which terrestrial planets that may harbour life form.

NATURE’S STARSHIPS. II. SIMULATING THE SYNTHESIS OF AMINO ACIDS IN METEORITE PARENT BODIES

Carbonaceous chondrite meteorites are known for having high water and organic material contents, including amino acids. Here we address the origin of amino acids in the warm interiors of their parent bodies (planetesimals) within a few million years of their formation, and connect this with the astrochemistry of their natal protostellar disks. We compute both the total amino acid abundance pattern as well as the relative frequencies of amino acids within the CM2 (e.g. Murchison) and CR2 chondrite subclasses based on Strecker reactions within these bodies. We match the relative frequencies to well within an order of magnitude among both CM2 and CR2 meteorites for parent body temperatures $<$ 200$^{\circ}$C. These temperatures agree with 3D models of young planetesimal interiors. We find theoretical abundances of approximately 7x10$^5$ parts-per-billion (ppb), which is in agreement with the average observed abundance in CR2 meteorites of 4$\pm$7x10$^5$, but an order of magnitude higher than the average observed abundance in CM2 meteorites of 2$\pm$2x10$^4$. We find that the production of hydroxy acids could be favoured over the production of amino acids within certain meteorite parent bodies (e.g. CI1, CM2) but not others (e.g. CR2). This could be due to the relatively lower NH$_3$ abundances within CI1 and CM2 meteorite parent bodies, which leads to less amino acid synthesis. We also find that the water content in planetesimals is likely to be the main cause of variance between carbonaceous chondrites of the same subclass. We propose that amino acid abundances are primarily dependent on the ammonia and water content of planetesimals that are formed in chemically distinct regions within their natal protostellar disks.

Particle splitting in smoothed particle hydrodynamics based on Voronoi diagram

We present a novel method for particle splitting in smoothed particle hydrodynamics simulations. Our method utilizes the Voronoi diagram for a given particle set to determine the position of fine daughter particles. We perform several test simulations to compare our method with a conventional splitting method in which the daughter particles are placed isotropically over the local smoothing length. We show that, with our method, the density deviation after splitting is reduced by a factor of about two compared with the conventional method. Splitting would smooth out the anisotropic density structure if the daughters are distributed isotropically, but our scheme allows the daughter particles to trace the original density distribution with length scales of the mean separation of their parent. We apply the particle splitting to simulations of the primordial gas cloud collapse. The thermal evolution is accurately followed to the hydrogen number density of 10^12 /cc. With the effective mass resolution of ~10^-4 Msun after the multi-step particle splitting, the protostellar disk structure is well resolved. We conclude that the method offers an efficient way to simulate the evolution of an interstellar gas and the formation of stars.

The Composition of the Protosolar Disk and the Formation Conditions for Comets

Conditions in the protosolar nebula have left their mark in the composition of cometary volatiles, thought to be some of the most pristine material in the solar system. Cometary compositions represent the end point of processing that began in the parent molecular cloud core and continued through the collapse of that core to form the protosun and the solar nebula, and finally during the evolution of the solar nebula itself as the cometary bodies were accreting. Disentangling the effects of the various epochs on the final composition of a comet is complicated. But comets are not the only source of information about the solar nebula. Protostellar disks around young stars similar to the protosun provide a way of investigating the evolution of disks similar to the solar nebula while they are in the process of evolving to form their own solar systems. In this way we can learn about the physical and chemical conditions under which comets formed, and about the types of dynamical processing that shaped the solar system we see today. This paper summarizes some recent contributions to our understanding of both cometary volatiles and the composition, structure and evolution of protostellar disks.

Dust trapping by spiral arms in gravitationally unstable protostellar discs

In this paper we discuss the influence of gravitational instabilities in massive protostellar discs on the dynamics of dust grains. Starting from a Smoothed Particle Hydrodynamics (SPH) simulation, we have computed the evolution of the dust in a quasi-static gas density structure typical of self-gravitating disc. For different grain size distributions we have investigated the capability of spiral arms to trap particles. We have run 3D radiative transfer simulations in order to construct maps of the expected emission at (sub-)millimetre and near-infrared wavelengths. Finally, we have simulated realistic observations of our disc models at (sub-)millimetre and near-infrared wavelengths as they may appear with the Atacama Large Millimetre/sub-millimetre Array (ALMA) and the High-Contrast Coronographic Imager for Adaptive Optics (HiCIAO) in order to investigate whether there are observational signatures of the spiral structure. We find that the pressure inhomogeites induced by gravitational instabilities produce a non-negligible dynamical effect on centimetre sized particles leading to significant overdensities in spiral arms. We also find that the spiral structure is readily detectable by ALMA over a wide range of (sub-)millimetre wavelengths and by HiCIAO in near-infrared scattered light for non-face-on discs located in the Ophiucus star-forming region. In addition, we find clear spatial spectral index variations across the disc, revealing that the dust trapping produces a migration of large grains that can be potentially investigated through multi-wavelenghts observations in the (sub-)millimetric. Therefore, the spiral arms observed to date in protoplanetary disc might be interpreted as density waves induced by the development of gravitational instabilities.

The doubly eclipsing quintuple low-mass star system 1SWASP J093010.78+533859.5

Our discovery of 1SWASP J093010.78+533859.5 as a probable doubly eclipsing quadruple system containing a contact binary with P~0.23 d and a detached binary with P~1.31 d was announced in 2013. Subsequently Koo et al. confirmed the detached binary spectroscopically and identified a fifth set of static spectral lines at its location, corresponding to a further non-eclipsing component of the system. Here we present new spectroscopic and photometric observations, allowing confirmation of the contact binary and improved modelling of all four eclipsing components. The detached binary is found to contain components of masses 0.837(8) and 0.674(7) M_sol, with radii of 0.832(18) and 0.669(18) R_sol and effective temperatures of 5185(-20,+25) and 4325(-15,+20) K respectively, the contact system has masses 0.86(2) and 0.341(11) M_sol, radii of 0.79(4) and 0.52(5) R_sol respectively, and a common T_eff of 4700(50) K. The fifth star is of similar temperature and spectral type to the primaries in the two binaries. Long-term photometric observations indicate the presence of a spot on one component of the detached binary, moving at an apparent rate of approximately one rotation every two years. Both binaries have consistent system velocities around -11 to -12 km/s, which match the average radial velocity of the fifth star, consistent distance estimates for both subsystems of d=78(3) and d=73(4) pc are also found, and (with some further assumptions) of d=83(9) pc for the fifth star. These findings strongly support the claim that both binaries (and very probably all five stars) are gravitationally bound in a single system. The consistent angles of inclination found for the two binaries (88.2(3) and 86(4) degrees) may also indicate that they originally formed by fragmentation (~9-10 Gyr ago) from a single protostellar disk and subsequently remained in the same orbital plane.