Protobinary Systems Research Articles

Multiple chemical tracers finally unveil the intricate NGC 1333 IRAS 4A outflow system. FAUST XVI

ABSTRACT The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS 4A protobinary system. This analysis is based on ALMA observations of SiO(5–4), H2CO(30, 3–20, 3), and HDCO(41, 4–31, 3) with a spatial resolution of ∼150 au. Leveraging an astrochemical approach involving the use of diverse tracers beyond traditional ones has enabled the identification of novel features and a comprehensive understanding of the broader outflow dynamics. Our analysis reveals the presence of two jets in the redshifted emission, emanating from IRAS 4A1 and IRAS 4A2, respectively. Furthermore, we identify four distinct outflows in the region for the first time, with each protostar, 4A1 and 4A2, contributing to two of them. We characterize the morphology and orientation of each outflow, challenging previous suggestions of bends in their trajectories. The outflow cavities of IRAS 4A1 exhibit extensions of 10 and 13 arcsec with position angles (PA) of 0° and -12°, respectively, while those of IRAS 4A2 are more extended, spanning 18 and 25 arcsec with PAs of 29° and 26°. We propose that the misalignment of the cavities is due to a jet precession in each protostar, a notion supported by the observation that the more extended cavities of the same source exhibit lower velocities, indicating they may stem from older ejection events.

Observations of high-order multiplicity in a high-mass stellar protocluster

The dominant mechanism forming multiple stellar systems in the high-mass regime (M* ≳ 8 M⊙) remained unknown because direct imaging of multiple protostellar systems at early phases of high-mass star formation is very challenging. High-mass stars are expected to form in clustered environments containing binaries and higher-order multiplicity systems. So far only a few high-mass protobinary systems, and no definitive higher-order multiples, have been detected. Here we report the discovery of one quintuple, one quadruple, one triple and four binary protostellar systems simultaneously forming in a single high-mass protocluster, G333.23–0.06, using Atacama Large Millimeter/submillimeter Array high-resolution observations. We present a new example of a group of gravitationally bound binary and higher-order multiples during their early formation phases in a protocluster. This provides the clearest direct measurement of the initial configuration of primordial high-order multiple systems, with implications for the in situ multiplicity and its origin. We find that the binary and higher-order multiple systems, and their parent cores, show no obvious sign of disk-like kinematic structure. We conclude that the observed fragmentation into binary and higher-order multiple systems can be explained by core fragmentation, indicating its crucial role in establishing the multiplicity during high-mass star cluster formation.

The impact of the initial core temperature on protostellar disc fragmentation

ABSTRACT Ground-based and satellite observations have revealed dust temperatures as low as ∼5−7 K in the centre of low-mass, pre-stellar cloud cores, where star formation takes place. However, external heating may rise the outer core temperatures up to ∼15−20 K. Such low temperatures at the centre of pre-stellar cores are a key factor to constrain the conditions that lead to the formation of gravitationally bound protostellar systems as was recently captured by highly-resolved Atacama large millimeter/submillimeter array observations. Here, we report consistent smoothed particle hydrodynamics collapse calculations of cold cores that demonstrate the formation of close protobinary systems via small-scale fragmentation of a gravitationally unstable protostellar disc. The results indicate that mean binary separations, of tens of astronomical units, are a consequence of disc fragmentation in cold pre-stellar cores. Cloud cores initially with temperatures ≤6 K and a low amplitude (a = 0.1), m = 2 density perturbation formed close protobinaries that were followed deep into the non-isothermal collapse for several orbital periods and appeared to survive as independent stellar entities. At temperatures ≥7 K disc fragmentation is no longer observed and the calculations terminate with the formation of a wide protobinary, which may occasionally be accompanied by small substellar objects emerging by fragmentation of the circumbinary disc. When the perturbation amplitude is raised to a = 0.25, disc fragmentation occurs again only in cores with initial temperatures ≤6 K. Therefore, increasing the perturbation amplitude does not necessarily imply that there will be disc fragmentation at higher core temperatures.

Turbulence and its connection to episodic accretion in binary YSOs

ABSTRACT We report signatures of episodic accretion in young stellar objects (YSOs) that emerge in protobinary configurations in a gravoturbulent gas collapse. We find in most of these protobinary systems strong accretion bursts between the two companions with a recurrence time-scale of about 1 kyr. The accretion rate on to the secondary star typically exceeds that on to the primary with a peak value of 2 × 10−2 M⊙ yr−1 for the former and 6 × 10−3 M⊙ yr−1 for the latter. We propose that the secondary companion, which remains more active in its episodes of accretion bursts, especially for the gas cores with subsonic velocity dispersion, may provide observational opportunities to find traces of episodic accretion in the surrounding gas of the embedded YSOs that are in a binary configuration. Also, protostars evolving as single objects in the same environment show fewer accretion bursts and all together a more steady mass growth history. The prestellar cores with subsonic velocity dispersion exhibit an order of magnitude more intense accretion bursts than in the case of cores with supersonic velocity dispersions. The latter shows the formation of some of the protobinaries in which the primary acts as a more actively accreting companion. This can support these binaries to become systems of extreme mass ratio. Moreover, the YSOs in binary configurations with small semimajor axis a ≈ 50 au and high mass ratio q > 0.7 support phases of intense episodic accretion. The eccentricity, however, seems to play no significant role in the occurrence of accretion bursts.

Misaligned Twin Molecular Outflows from the Class 0 Protostellar Binary System VLA 1623A Unveiled by ALMA

Abstract We present the results of ALMA observations toward the low-mass Class 0 binary system VLA 1623Aab in the Ophiuchus molecular cloud in 12CO, 13CO, and C18O(2–1) lines. Our 12CO (J = 2–1) data reveal that the VLA 1623 outflow consists of twin spatially overlapped outflows/jets. The redshifted northwestern jet exhibits three cycles of wiggle with a spatial period of 1360 ± 10 au, corresponding to a time period of 180 yr. The wiggle-like structure is also found in the position–velocity (PV) diagram, showing an amplitude in the velocity of about 0.9 km s−1. Both the period and velocity amplitude of the wiggle are roughly consistent with those expected from the binary parameters, i.e., the orbital period (460 ± 20 yr) and the Keplerian velocity (2.2 km s−1). Our 13CO and C18O images show a dense gas nature in the two centimeter/millimeter sources, VLA 1623B and W, and its relation to the outflows, and strongly support the previous interpretation that both are shocked cloudlets. The driving sources of the twin molecular outflows are, therefore, likely to be the VLA 1623Aab binary. The outflow axes of the two molecular outflows are estimated to be inclined by 70° to each other across the plane of sky, implying that protostellar disks are also misaligned by . Such nature together with a small binary separation of 34 au in one of the youngest protobinary systems seems difficult to explain by disk fragmentation in quiescent environments. Other effects such as turbulence probably play roles.

Twin Jets and Close Binary Formation

Abstract The formation of a close binary system is investigated using a three-dimensional resistive magnetohydrodynamic simulation. Starting from a prestellar cloud, the cloud evolution is calculated until ∼400 yr after protostar formation. Fragmentation occurs in the gravitationally collapsing cloud, and two fragments evolve into protostars. The protostars orbit each other and a protobinary system appears. A wide-angle low-velocity outflow emerges from the circumbinary streams that enclose two protostars, while each protostar episodically drives high-velocity jets. Thus, the two high-velocity jets are surrounded by the low-velocity circumbinary outflow. The speed of the jets exceeds ≳100 . Although the jets have a collimated structure, they are swung back on the small scale and are tangled at the large scale due to the binary orbital motion. A circumstellar disk also appears around each protostar. In the early main accretion phase, the binary orbit is complicated, while the binary separation is within <30 au. For the first time, all the characteristics of protobinary systems recently observed with large telescopes are reproduced in a numerical simulation.

Gas flow and accretion via spiral streamers and circumstellar disks in a young binary protostar.

The majority of stars are part of gravitationally bound stellar systems, such as binaries. Observations of protobinary systems constrain the conditions that lead to stellar multiplicity and subsequent orbital evolution. We report high-angular resolution observations of the circumbinary disk around [BHB2007] 11, a young binary protostar system. The two protostars are embedded in circumstellar disks that have radii of 2 to 3 astronomical units and probably contain a few Jupiter masses. These systems are surrounded by a complex structure of filaments connecting to the larger circumbinary disk. We also observe accretion and radio jets associated with the protobinary system. The accretion is preferentially onto the lower-mass protostar, consistent with theoretical predictions.

The formation of protostellar binaries in primordial minihalos

The first stars are known to form in primordial gas, either in minihalos with about $10^6$~M$_\odot$ or so-called atomic cooling halos of about $10^8$~M$_\odot$. Simulations have shown that gravitational collapse and disk formation in primordial gas yield dense stellar clusters. In this paper, we focus particularly on the formation of protostellar binary systems, and aim to quantify their properties during the early stage of their evolution. For this purpose, we combine the smoothed particle hydrodynamics code GRADSPH with the astrochemistry package KROME. The GRADSPH-KROME framework is employed to investigate the collapse of primordial clouds in the high-density regime, exploring the fragmentation process and the formation of binary systems. We observe a strong dependence of fragmentation on the strength of the turbulent Mach number $\mathcal{M}$ and the rotational support parameter $\beta{}$. Rotating clouds show significant fragmentation, and have produced several Pop.~III proto-binary systems. We report maximum and minimum mass accretion rates of $2.31 \times 10^{-1}$~M$_{\odot}$ yr$^{-1}$ and $2.18\times 10^{-4}$~M$_{\odot}$ yr$^{-1}$. The mass spectrum of the individual Pop III proto-binary components ranges from $0.88$~M$_{\odot}$ to $31.96$~M$_{\odot}$ and has a sensitive dependence on the Mach number $\mathcal{M}$ as well as on the rotational parameter $\beta{}$. We also report a range from $\sim0.01$ to $\sim1$ for the mass ratio of our proto-binary systems.

A SUB-ARCSECOND SURVEY TOWARD CLASS 0 PROTOSTARS IN PERSEUS: SEARCHING FOR SIGNATURES OF PROTOSTELLAR DISKS

We present a CARMA 1.3 mm continuum survey toward 9 Class 0 protostars in the Perseus molecular cloud at $\sim$0.3$^{\prime\prime}$ (70 AU) resolution. This study approximately doubles the number of Class 0 protostars observed with spatial resolutions $<$ 100 AU at millimeter wavelengths, enabling the presence of protostellar disks and proto-binary systems to be probed. We detect flattened structures with radii $>$ 100 AU around 2 sources (L1448 IRS2 and Per-emb-14) and these sources may be strong disk candidates. Marginally-resolved structures with position angles within 30$^{\circ}$ of perpendicular to the outflow are found toward 3 protostars (L1448 IRS3C, IRAS 03282+3035, and L1448C) and are considered disk candidates. Two others (L1448 IRS3B and IRAS 03292+3039) have resolved structure, possibly indicative of massive inner envelopes or disks; L1448 IRS3B also has a companion separated by 0.9$^{\prime\prime}$ ($\sim$210 AU). IC348-MMS does not have well-resolved structure and the candidate first hydrostatic core L1451-MMS is marginally resolved on 1$^{\prime\prime}$ scales. The strong disk candidate sources were followed-up with C$^{18}$O ($J=2\rightarrow1$) observations, detecting velocity gradients consistent with rotation, but it is unclear if the rotation is Keplerian. We compare the observed visibility amplitudes to radiative transfer models, finding that visibility amplitude ratios suggest a compact component (possibly a disk) is necessary for 5 of 9 Class 0 sources; envelopes alone may explain the other 4 systems. We conclude that there is evidence for the formation of large disks in the Class 0 phase with a range of radii and masses dependent upon their initial formation conditions.

SMA OBSERVATIONS OF CLASS 0 PROTOSTARS: A HIGH ANGULAR RESOLUTION SURVEY OF PROTOSTELLAR BINARY SYSTEMS

We present high angular resolution 1.3 mm and 850 um dust continuum data obtained with the Submillimeter Array toward 33 Class 0 protostars in nearby clouds (distance < 500 pc), which represents so far the largest survey toward protostellar binary/multiple systems. The median angular resolution in the survey is 2.5 arcsec, while the median linear resolution is approximately 600 AU. Compact dust continuum emission is observed from all sources in the sample. Twenty-one sources in the sample show signatures of binarity/multiplicity, with separations ranging from 50 to 5000 AU. The numbers of singles, binaries, triples, and quadruples in the sample are 12, 14, 5, and 2, respectively. The derived multiplicity frequency (MF) and companion star fraction (CSF) for Class 0 protostars are 0.64+/-0.08 and 0.91+/-0.05, respectively, with no correction for completeness. The derived MF and CSF in this survey are approximately two times higher than the values found in the binary surveys toward Class I YSOs, and approximately three (for MF) and four (for CSF) times larger than the values found among MS stars, with a similar range of separations. Furthermore, the observed fraction of high order multiple systems to binary systems in Class 0 protostars (0.50+/-0.09) is also larger than the fractions found in Class I YSOs (0.31+/-0.07) and MS stars (< 0.2). These results suggest that binary properties evolve as protostars evolve, as predicted by numerical simulations. The distribution of separations for Class 0 protostellar binary systems shows a general trend in which companion star fraction increases with decreasing companion separation. We find that 67%+/-8% of the protobinary systems have circumstellar mass ratios below 0.5, implying that unequal-mass systems are preferred in the process of binary star formation. We suggest an empirical sequential fragmentation picture for binary star formation.

Toward understanding the formation of multiple systems

The formation process of binary stars and multiple systems is poorly understood. Here, we seek to determine the typical outcome of protostellar collapse and to constrain models of binary formation by core fragmentation during collapse, using high-resolution millimeter continuum imaging of very young (Class 0) protostars observed at the beginning of the main accretion phase. We carried out a pilot high-resolution study of 5 Class 0 objects, using the most extended (A) configuration of the IRAM PdBI at 1.3 mm, which allow us to probe the multiplicity of Class 0 protostars down to separations a ~50 AU and circumstellar mass ratios q ~0.07. We show that our PdBI observations revealed only wide (>1500 AU) protobinary systems and/or outflow-generated features. When combined with previous millimeter interferometric observations of Class 0 protostars, our pilot PdBI study tentatively suggests that the binary fraction in the ~ 75-1000 AU range increases from the Class 0 to the Class I stage. It also seems to argue against purely hydrodynamic models of binary star formation. We briefly discuss possible alternative scenarios to reconcile the low multiplicity rate of Class 0 protostars on small scales with the higher binary fraction observed at later evolutionary stages.

Orbital Motions in Binary Protostellar Systems

AbstractUsing high-resolution (~ 01), multi-epoch Very Large Array observations, we have detected orbital motions in several low-luminosity protobinary systems in the Taurus and ρ Ophiuchus molecular complexes. The masses obtained from Kepler’s third law are of the order of 0.5 to 2 M⊙, as expected for such low-mass protostars. The relatively large bolometric luminosities of these young systems corroborates the notion that protostars obtain most of their luminosity from accretion and not from nuclear reactions. In addition, in one of the sources studied (a multiple system in Taurus), a low-mass young star has shown a drastic change in its orbit after a close approach with another component of the system, presumed to be a double star. The large proper motion achieved by this low mass protostar (20 km s−1), suggests an ejection from the system.

Orbital Proper Motions in the Protobinary System L1527/IRAS 04368+2557?

Using high angular resolution 7 mm Very Large Array (VLA) observations, we show that the low-mass protostellar object IRAS 043682557 in L1527 is composed of two radio sources separated in projection by 25 AU. One of the two components is extended in the direction perpendicular to the outflow powered by IRAS 043682557 and is, most likely, a compact accretion disk ( AU), similar to those found in L1551 R ∼ 20 IRS 5 by Rodro ´guez et al. As in L1551 IRS 5, the disk found here is small compared to those around T Tauri stars. Tidal interactions with nearby companions provide a natural way of truncating disks, and we argue that the two millimeter-wavelength sources in L1527 trace a compact binary system, where the disk surrounding one of the components has been truncated by the tidal influence of the other. A comparison between observations obtained in 1996 and 2002 reveals large proper motions, which can be only partly attributed to the overall largescale motion of the region in which IRAS 043682557 is located. The remaining “residual” proper motions might trace the orbital motion of the binary and would suggest a total mass for the system larger than 0.2 M, and likely of the order of 0.5–2 M,. This mass (only a small fraction of which is in the disk) is of the same order as that of the extended surrounding envelope of gas and dust traced by far-infrared and submillimeter observations ( M,), implying that the stars that will eventually form out of IRAS 043682557 would M ∼ 0.5 env have already acquired a significant fraction of their final masses. It is worth noting that multiepoch VLA studies of nearby protobinary systems similar to that presented here could provide direct mass estimates in most nearby star-forming sites. Combined with submillimeter observations of the surrounding envelopes, this would provide a more reliable measure of the evolutionary status of binary protostars. Subject headings: astrometry — binaries: general — ISM: jets and outflows — stars: formation

Observations of Protobinary Systems

High resolution images at millimeter wavelengths are providing new insights into the formation of binary and multiple star systems. These wavelengths are particularly useful in studying the earliest stages of multiplicity because they trace the bulk material distribution in the circumstellar environment and can penetrate 1000's of magnitudes of visual extinction. Current millimeter wavelength observations are finding a high incidence of multiplicity among young systems and that multiplicity begins at birth, or before. While the statistics are poor, the types of systems found (independent envelope, common envelope, and common disk systems) follow theoretical ideas about binary and multiple star formation. Systems can be identified which exhibit the characteristics of prompt initial collapse, central fragmentation during collapse, and fragmentation in high angular momentum scenarios. Expansion of this work to more systems and to more detailed studies of the structure and kinematics of individual systems will provide valuable insights into the formation of multiple systems.

Accretion and the Properties of Protobinary Systems

We discuss the evolution of the properties of a protobinary system which forms via fragmentation within a collapsing molecular cloud core and grows to its final mass by accreting from the gaseous envelope. We predict how the mass ratio distributions and the masses of the circumstellar and circumbinary discs depend on the binary's separation, the primary's mass, and the properties of the pre-collapse molecular cloud core. We also discuss how the properties of binaries may depend on whether they were formed in isolated molecular cloud cores or in stellar clusters.

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.

Modelling accretion in protobinary systems

A method for following fragmentation simulations further in time using smoothed particle hydrodynamics (SPH) is presented. In a normal SPH simulation of the collapse and fragmentation of a molecular cloud, high-density regions of gas that form protostars are represented by many particles with small separations. These high-density regions require small time steps, limiting the time for which the simulation can be followed. Thus, the end result of the fragmentation can never be definitively ascertained, and comparisons between cloud fragmentation calculations and the observed characteristics of stellar systems cannot be made. In this paper, each high-density region is replaced by a single, non-gaseous particle, with appropriate boundary conditions, which contains all the mass in the region and accretes any infalling mass. This enables the evolution of the cloud and the resulting protostars to be followed for many orbits or until most of the original cloud mass has been accreted.

DYNAMICS OF BINARY AND PLANETARY SYSTEM INTERACTION WITH DISKS: ECCENTRICITY CHANGES

ABSTRACT Protostellar and protoplanetary systems, as well as merging galactic nuclei, often interact tidally and resonantly with the astrophysical disks via gravity. Underlying our understanding of the formation processes of stars, planets and some galaxies is a dynamical theory of such interactsion. Its main goals are to determine the geometry of the binary-disk system and, through the torque calculations, the rate of change of orbital elements of the components. We present some recent developments in this field concentrating on eccentricity driving mechanisms in protoplanetary and protobinary systems. In those two types of systems the result of the interaction is opposite. A small body embedded in a disk suffers a decrease_ of orbital eccentricity, whereas newly formed binary stars surrounded by protostellar disks may undergo a significant orbital evolution increasing_ their eccentricities.