Theory Of Star Formation Research Articles

Kalkayotl 2.0. Bayesian phase-space modelling of star-forming regions, stellar associations, and open clusters

Star-forming regions, stellar associations, and open clusters are fundamental stellar systems where predictions from star-formation theories can be robustly contrasted with observations. We aim to provide the astrophysical community with a free and open-source code to infer the phase-space (i.e. positions and velocities) parameters of stellar systems with lesssim 1000 stars based on Gaia astrometry and possibly observed radial velocities. We upgrade an existing Bayesian hierarchical model and extend it to model 3D (positions) and 6D (positions and velocities) stellar coordinates and system parameters with a flexible variety of statistical models, including a linear velocity field. This velocity field allows for the inference of internal kinematics, including expansion, contraction, and rotation. We extensively validated our statistical models using realistic simulations that mimic the properties of the Gaia Data Release 3. We applied Kalkayotl to beta -Pictoris, the Hyades, and Praesepe, recovering parameter values compatible with those from the literature. In particular, we found an expansion age of $19.1 Myr for beta -Pictoris and rotational signal of $32\! $ for the Hyades and that Praesepe's rotation reported in the literature comes from its periphery. The robust and flexible Bayesian hierarchical model that we make publicly available here represents a step forward in the statistical modelling of stellar systems. The products it delivers, such as expansion, contraction, rotation, and velocity dispersions, can be directly contrasted with predictions from star-formation theories.

Disk Turbulence and Star Formation Regulation in High-z Main-sequence Analog Galaxies

The gas-phase velocity dispersions in disk galaxies, which trace turbulence in the interstellar medium, are observed to increase with lookback time. However, the mechanisms that set this rise in turbulence are observationally poorly constrained. To address this, we combine kiloparsec-scale Atacama Large Millimeter/submillimeter Array observations of CO(3−2) and CO(4−3) with Hubble Space Telescope observations of Hα to characterize the molecular gas and star formation properties of seven local analogs of main-sequence galaxies at z ∼ 1–2, drawn from the DYNAMO sample. Investigating the “molecular gas main sequence” on kiloparsec scales, we find that galaxies in our sample are more gas-rich than local star-forming galaxies at all disk positions. We measure beam-smearing-corrected molecular gas velocity dispersions and relate them to the molecular gas and star formation rate surface densities. Despite being relatively nearby (z ∼ 0.1), DYNAMO galaxies exhibit high velocity dispersions and gas and star formation rate surface densities throughout their disks, when compared to local star-forming samples. Comparing these measurements to predictions from star formation theory, we find very good agreements with the latest feedback-regulated star formation models. However, we find that theories that combine dissipation of gravitational energy from radial gas transport with feedback overestimate the observed molecular gas velocity dispersions.

Simulating Brown Dwarf Observations for Various Mass Functions, Birthrates, and Low-mass Cutoffs

After decades of brown dwarf discovery and follow-up, we can now infer the functional form of the mass distribution within 20 pc, which serves as a constraint on star formation theory at the lowest masses. Unlike objects on the main sequence that have a clear luminosity-to-mass correlation, brown dwarfs lack a correlation between an observable parameter (luminosity, spectral type, or color) and mass. A measurement of the brown dwarf mass function must therefore be procured through proxy measurements and theoretical models. We utilize various assumed forms of the mass function, together with a variety of birthrate functions, low-mass cutoffs, and theoretical evolutionary models, to build predicted forms of the effective temperature distribution. We then determine the best fit of the observed effective temperature distribution to these predictions, which in turn reveals the most likely mass function. We find that a simple power law ( dN/dM∝M−α ) with α ≈ 0.5 is optimal. Additionally, we conclude that the low-mass cutoff for star formation is ≲0.005 M ⊙. We corroborate the findings of Burgasser, which state that the birthrate has a far lesser impact than the mass function on the form of the temperature distribution, but we note that our alternate birthrates tend to favor slightly smaller values of α than the constant birthrate. Our code for simulating these distributions is publicly available. As another use case for this code, we present findings on the width and location of the subdwarf temperature gap by simulating distributions of very old (8–10 Gyr) brown dwarfs.

Absence of High-mass Prestellar Cores in the Orion Giant Molecular Cloud

A fundamental difference between “core-fed” and “clump-fed” star-formation theories lies in the existence or absence of high-mass cores at the prestellar stage. However, only a handful of such cores have been observed. Here, different than previous search in distributed star-formation regions in the Galactic plane, we search for high-mass prestellar cores in the Orion GMC, by observing the seven most massive starless cores selected from previous deep continuum surveys. We present ALMA Atacama Compact Array Band 6 and Band 7 continuum and line observations toward the seven cores, in which we identify nine dense cores at both bands. The derived maximum core mass is less than 11 M ⊙, based on different dust temperatures. We find no high-mass prestellar cores in this sample, aligning with the results of previous surveys, thereby challenging the existence of such cores in Orion. Outside Orion, further detailed studies are needed for remaining high-mass prestellar core candidates to confirm their status as massive, starless cores.

Detection of Two Totally Eclipsing B-type Binaries with Extremely Low Mass Ratios

The detection of O- and B-type stars with extremely low-mass companions is very important for understanding the formation and evolution of binary stars. However, finding them remains a challenge because the low-mass components in such systems contribute such small flux to the total. During our search for pulsations among O- and B-type stars using the Transiting Exoplanet Survey Satellite (TESS) data, we found two short-period and B-type (B9) eclipsing binaries with orbital periods of 1.61613 and 2.37857 days. Photometric solutions of the two close binaries were derived by analyzing the TESS light curves with the Wilson–Devinney method. It is discovered that both of them are detached binaries with extremely low mass ratios of 0.067(2) for TIC 260342097 and 0.140(3) for TIC 209148631. The determined mass ratio indicates that TIC 260342097 is one of the lowest mass ratios among known B-type binary systems. We showed that the two systems have total eclipses with a broad and flat secondary minimum, suggesting that the photometric parameters could be derived reliably. The absolute parameters of the two binaries are estimated and it is found that the secondary components in the two systems are overluminous and oversize when compared with the normal low-mass and cool main-sequence (MS) stars. These findings may imply that the two systems are composed of a B-type MS primary and a cool pre-MS secondary with orbital periods shorter than 2.5 days. They are valuable targets to test theories of binary star formation and evolution.

Metallicity Dependence of Pressure-regulated Feedback-modulated Star Formation in the TIGRESS-NCR Simulation Suite

We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over Z′≡Z/Z⊙∼0.1-3 , gas surface density Σgas ∼ 5–150 M ⊙ pc−2, and stellar surface density Σstar ∼ 1–50 M ⊙ pc−2. The range of emergent star formation rate surface density is ΣSFR ∼ 10−4–0.5 M ⊙ kpc−2 yr−1, and ISM total midplane pressure is P tot/k B = 103–106 cm−3 K, with P tot equal to the ISM weight W . For given Σgas and Σstar, we find ΣSFR∝Z′0.3 . We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to ΣSFR. The thermal feedback yield varies sensitively with both weight and metallicity as ϒth∝W−0.46Z′−0.53 , while the combined turbulent and magnetic feedback yield shows weaker dependence ϒturb+mag∝W−0.22Z′−0.18 . The reduction in ΣSFR at low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved.

A High-resolution Imaging Survey of Massive Young Stellar Objects in the Magellanic Clouds

Constraints on the binary fraction of massive young stellar objects (mYSOs) are important for binary and massive star formation theory. Here, we present speckle imaging of 34 mYSOs located in the Large Magellanic Cloud (1/2 Z ⊙) and Small Magellanic Cloud (∼1/5 Z ⊙), probing projected separations in the 2000 to 20,000 au (at angular scales of 0.″02–0.″2) range, for stars above 8 M ⊙. We find two wide binaries in the Large Magellanic Cloud (from a sample of 23 targets), but none in a sample of 11 in the Small Magellanic Cloud, leading us to adopt a wide binary fraction of 9% ± 5% and <5%, respectively. We rule out a wide binary fraction greater than 35% in the Large Magellanic Cloud and 38% in the Small Magellanic Cloud at the 99% confidence level. This is in contrast to the wide binary fraction of mYSOs in the Milky Way (presumed to be 1 Z ⊙), which within the physical parameter space probed by this study is ∼15%–60% from the literature. We argue that while selection effects could be responsible for the lower binary fraction observed, it is more likely that there are underlying physical mechanisms responsible for the observed properties. This indicates that metallicity and environmental effects may influence the formation of wide binaries among massive stars. Future larger, more statistically significant samples of high-mass systems in low-metallicity environments for comparison to the Milky Way, are essential to confirm or repudiate our claim.

Collapsing molecular clouds with tracer particles – II. Collapse histories

ABSTRACT In order to develop a complete theory of star formation, one essentially needs to know two things: what collapses and how long it takes. This is the second paper in a series, where we query how long a parcel of gas takes to collapse and the process it undergoes. We embed pseudo-Lagrangian tracer particles in simulations of collapsing molecular clouds, identify the particles that end in dense knots, and then examine the collapse history of the gas. We find a nearly universal behaviour of cruise-then-collapse, wherein a core stays at intermediate densities for a significant fraction of its life before finally collapsing. We identify time immediately before each core collapses, $t_{\rm {sing}}$, and examine how it transitions to high density. We find that the time to collapse is uniformly distributed between $0.25 t_{\rm {ff}}$ and the end of the simulation at $\sim\!\! 1 t_{\rm {ff}}$, and that the duration of collapse is universally short, $\Delta t \sim 0.1 t_{\rm {ff}}$, where $t_{\rm {ff}}$ is the free-fall time at the mean density. We describe the collapse in three stages: collection, hardening, and singularity. Collection sweeps low-density gas into moderate density. Hardening brings kinetic and gravitational energies into quasi-equipartition. Singularity is the free-fall collapse, forming an envelope in rough energy balance and central overdensity in $\sim\!\! 0.1 t_{\rm {ff}}$.

A Consistent Explanation for the Unusual Initial Mass Function and Star Formation Rate in the Central Molecular Zone (CMZ)

We examine various physical processes that may explain the shallow high-mass slope of the initial mass function (IMF), as well as the low star formation rate (SFR) in star-forming molecular clouds (MCs) in the Central Molecular Zone (CMZ). We show that the strong tidal field and shear experienced by the CMZ have opposite effects on the collapse of density fluctuations and cannot explain these properties. Similarly, we show that the intense magnetic field in the CMZ provides a negligible pressure support and, for the high densities at play, should not modify the probability density function of the turbulent gas flow, thus affecting negligibly the IMF. However, we show that, in contrast to the MCs in the Galactic disk, the ones in the CMZ experience only one single episode of turbulence cascade. Indeed, their rather short lifetime, due to their high mean densities, is similar to one typical turbulence crossing time. Consequently, according to the Hennebelle–Chabrier theory of star formation, within this “single turbulence cascade episode,” the cloud experiences one single field of turbulence-induced density fluctuations, leading eventually to gravitationally unstable cores. As shown in Hennebelle & Chabrier (2013), this yields a shallower IMF than usual and leads to the correct observed slope for the CMZ star-forming clouds. Similarly, this single large-scale turbulence event within the cloud lifetime yields a 5–6 times lower SFR than under usual conditions, in agreement with the observed values. Therefore, we suggest that this “single turbulence cascade” scenario can explain both the shallow IMF and low SFR of clouds in the CMZ.

TYC 3340-2437-1: A Quadruple System with a Massive Star

Hierarchical massive quadruple systems are ideal laboratories for examining the theories of star formation, dynamical evolution, and stellar evolution. The successive mergers of hierarchical quadruple systems might explain the mass gap between neutron stars and black holes. Looking for light curves of O-type binaries identified by LAMOST, we find a (2+2) quadruple system: TYC 3340-2437-1, located in the stellar bow-shock nebula (SBN). It has a probability of over 99.99% being a quadruple system derived from the surface density of the vicinity stars. Its inner orbital periods are 3.390602(89) days and 2.4378(16) days, respectively, and the total mass is about (11.47 + 5.79) + (5.2 + 2.02) = 24.48 M ⊙. The line-of-sight inclinations of the inner binaries, B1 and B2, are 55.°94 and 78.°2, respectively, indicating that they are not coplanar. Based on observations spanning 34 months and the significance of the astrometric excess noise (D > 2) in Gaia Data Release 3 data, we guess that its outer orbital period might be a few years. If it were true, the quadruple system might form through the disk fragmentation mechanism with outer eccentric greater than zero. This eccentricity could be the cause of both the arc-like feature of the SBN and the noncoplanarity of the inner orbit. The outer orbital period and outer eccentric could be determined with the release of future epoch astrometric data of Gaia.

An Updated Modular Set of Synthetic Spectral Energy Distributions for Young Stellar Objects

Measured properties of young stellar objects (YSOs) are key tools for research into pre-main-sequence stellar evolution. YSO properties are commonly measured by comparing observed radiation to existing grids of template YSO spectral energy distributions (SEDs) calculated by radiative transfer. These grids are often sampled and constructed using simple models of mass assembly/accretion over time. However, because we do not yet have a complete theory of star formation, the choice of model sets the tracked parameters and range of allowed values. By construction, then, the assumed model limits the measurements that can be made using the grid. Radiative transfer models not constrained by specific accretion histories would enable assessment of a wider range of theories. We present an updated version of the Robitaille set (2017) of YSO SEDs, a collection of models with no assumed evolutionary theory. We outline our newly calculated properties: envelope mass, weighted-average dust temperature, disk stability, and circumstellar A V. We also convolve the SEDs with new filters, including JWST, and provide users the ability to perform additional convolutions. We find a correlation between the average temperature and millimeter-wavelength brightness of optically thin dust in our models and discuss its ramifications for mass measurements of pre- and protostellar cores. We also compare the positions of YSOs of different observational classes and evolutionary stages in IR color space and use our models to quantify the extent to which class and stage may be confused due to observational effects. Our updated models are released to the public.

A Radio Counterpart to a Jupiter-mass Binary Object in Orion

Using James Webb Space Telescope near-infrared data of the inner Orion Nebula, Pearson & McCaughrean detected 40 Jupiter-mass binary objects (JuMBOs). These systems are not associated with stars and their components have masses of giant Jupiter-like planets and separations in the plane of the sky of order ∼100 au. The existence of these wide free-floating planetary-mass binaries was unexpected in our current theories of star and planet formation. Here we report the radio continuum (6.1 and 10.0 GHz) Karl G. Jansky Very Large Array detection of a counterpart to JuMBO 24. The radio emission appears to be steady at a level of ∼50 μJy over timescales of days and years. We set an upper limit of ≃15 km s−1 to the velocity of the radio source in the plane of the sky. As in the near-infrared, the radio emission seems to be coming from both components of the binary.

The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE): Core Growth, Cluster Contraction, and Primordial Mass Segregation

The Atacama Large Millimeter/submillimeter Array (ALMA) Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive (M clump ≳ 103 M ⊙), luminous (L bol ≳ 104 L ⊙), and blue-profile (infall signature) clumps by ALMA with resolution of ∼2200–5500 au (median value of 3500 au) at 350 GHz (870 μm). We identified 248 dense cores, including 106 cores showing protostellar signatures and 142 prestellar core candidates. Compared to early stage infrared dark clouds (IRDCs) by ASHES, the core mass and surface density within the ASSEMBLE clumps exhibited a significant increment, suggesting concurrent core accretion during the evolution of the clumps. The maximum mass of prestellar cores was found to be 2 times larger than that in IRDCs, indicating that evolved protoclusters have the potential to harbor massive prestellar cores. The mass relation between clumps and their most massive core (MMCs) is observed in ASSEMBLE but not in IRDCs, which is suggested to be regulated by multiscale mass accretion. The mass correlation between the core clusters and their MMCs has a steeper slope compared to that observed in stellar clusters, which can be due to fragmentation of the MMC and stellar multiplicity. We observe a decrease in core separation and an increase in central concentration as protoclusters evolve. We confirm primordial mass segregation in the ASSEMBLE protoclusters, possibly resulting from gravitational concentration and/or gas accretion.

The multiplicity of massive stars in the Scorpius OB1 association through high-contrast imaging

Context. Despite past efforts, a comprehensive theory of massive star formation is still lacking. One of the most remarkable properties of massive stars is that almost all of them are found in binaries or higher-order multiple systems. Since multiplicity is an important outcome parameter of a star formation process, observations that cover the full companion mass ratio and separation regime are essential to constrain massive star formation theories. Aims. We aim to characterise the multiplicity properties of 20 OB stars (one of which turned out to be a foreground object) in the active star-forming region Sco OB1 in the separation range 0.″15–6″ (∼200−9000 AU), using high-contrast imaging observations. These observations enabled us to reach very large magnitude differences and explore an as of yet uncharted territory of companions around massive stars. Methods. We used VLT/SPHERE to simultaneously observe with IFS and IRDIS, obtaining high-contrast imaging observations that cover a field of view (FoV) of 1.″73 × 1.″73 in YJH bands and 11″ × 12.″5 in K1 and K2 bands, respectively. We extracted low-resolution IFS spectra of candidate companions within 0.″85 and compared them with PHOENIX and ATLAS9 atmosphere models to obtain an estimate of their fundamental parameters. Furthermore, we retrieved an estimate of the mass and age of all sources in the larger IRDIS FoV. The observations reached contrast magnitudes of ΔK1 ∼ 13 on average, so we were able to detect sources at the stellar-substellar boundary. Results. In total, we detect 789 sources, most of which are likely background sources. Thirty objects that are estimated to be in the stellar mass regime have a 20% or lower probability of being spurious associations. We obtain SPHERE companion fractions of 2.3 ± 0.4 and 4.1 ± 0.8 for O- and B-type stars, respectively. Splitting the sample between more massive (&gt; 20 M⊙) and less massive stars (&lt; 20 M⊙), we obtain companion fractions of 2.3 ± 0.4 and 3.9 ± 0.7, respectively. Including all previously detected companions, we find a total multiplicity fraction of 0.89 ± 0.07 for separations in the range of ∼0−12 000 AU. Conclusions. With SPHERE, we are probing an unexplored area in the magnitude versus separation diagram of companions, which is crucial to achieve a complete overview of the multiplicity properties of massive stars and ultimately improve our understanding of massive star formation.

Mother of dragons

Context. Core accretion models of massive star formation require the existence of massive, starless cores within molecular clouds. Yet, only a small number of candidates for such truly massive, monolithic cores are currently known. Aims. Here we analyse a massive core in the well-studied infrared-dark cloud (IRDC) called the ‘dragon cloud’ (also known as G028.37+00.07 or ‘Cloud C’). This core (C2c1) sits at the end of a chain of a roughly equally spaced actively star-forming cores near the center of the IRDC. Methods. We present new high-angular-resolution 1 mm ALMA dust continuum and molecular line observations of the massive core. Results. The high-angular-resolution observations show that this region fragments into two cores, C2c1a and C2c1b, which retain significant background-subtracted masses of 23 M⊙ and 2 M⊙ (31 M⊙ and 6 M⊙ without background subtraction), respectively. The cores do not appear to fragment further on the scales of our highest-angular-resolution images (0.2″, 0.005 pc ~ 1000 AU). We find that these cores are very dense (nH2 &gt; 106 cm−3) and have only trans-sonic non-thermal motions (ℳs ~ 1). Together the mass, density, and internal motions imply a virial parameter of &lt;1, which suggests the cores are gravitationally unstable, unless supported by strong magnetic fields with strengths of ~1–10 mG. From CO line observations, we find that there is tentative evidence for a weak molecular outflow towards the lower-mass core, and yet the more massive core remains devoid of any star formation indicators. Conclusions. We present evidence for the existence of a massive, pre-stellar core, which has implications for theories of massive star formation. This source warrants follow-up higher-angular-resolution observations to further assess its monolithic and pre-stellar nature.

The cosmic waltz of Coma Berenices and Latyshev 2 (Group X)

Context. Open clusters (OCs) are fundamental benchmarks where theories of star formation and stellar evolution can be tested and validated. Coma Berenices (Coma Ber) and Latyshev 2 (Group X) are the second and third OCs closest to the Sun, making them excellent targets to search for low-mass stars and ultra-cool dwarfs. In addition, this pair will experience a flyby in 10–16 Myr, making it a benchmark to test pair interactions of OCs. Aims. We aim to analyse the membership, luminosity, mass, phase-space (i.e. positions and velocities), and energy distributions for Coma Ber and Latyshev 2 and test the hypothesis of the mixing of their populations at the encounter time. Methods. We developed a new phase-space membership methodology and applied it to Gaia data. With the recovered members, we inferred the phase-space, luminosity, and mass distributions using publicly available Bayesian inference codes. Then, with a publicly available orbit integration code and members’ positions and velocities, we integrated their orbits 20 Myr into the future. Results. In Coma Ber, we identified 302 candidate members distributed in the core and tidal tails. The tails are dynamically cold and asymmetrically populated. The stellar system called Group X is made of two structures: the disrupted OC Latyshev 2 (186 candidate members) and a loose stellar association called Mecayotl 1 (146 candidate members), and both of them will fly by Coma Ber in 11.3 ± 0.5 Myr and 14.0 ± 0.6 Myr, respectively, and each other in 8.1 ± 1.3 Myr. Conclusions. We study the dynamical properties of the core and tails of Coma Ber and also confirm the existence of the OC Latyshev 2 and its neighbour stellar association Mecayotl 1. Although these three systems will experience encounters, we find no evidence supporting the mixing of their populations.

The cosmic DANCe of Perseus

Context. Star-forming regions are excellent benchmarks for testing and validating theories of star formation and stellar evolution. The Perseus star-forming region, being one of the youngest (&lt; 10 Myr), closest (280−320 pc), and most studied in the literature, is a fundamental benchmark. Aims. We aim to study the membership, phase-space structure, mass, and energy (kinetic plus potential) distribution of the Perseus star-forming region using public catalogues (Gaia, APOGEE, 2MASS, and Pan-STARRS). Methods. We used Bayesian methodologies that account for extinction to identify the Perseus physical groups in the phase-space, retrieve their candidate members, derive their properties (age, mass, 3D positions, 3D velocities, and energy), and attempt to reconstruct their origin. Results. We identify 1052 candidate members in seven physical groups (one of them new) with ages between 3 and 10 Myr, dynamical super-virial states, and large fractions of energetically unbounded stars. Their mass distributions are broadly compatible with that of Chabrier for masses ≳0.1 M⊙ and do not show hints of over-abundance of low-mass stars in NGC 1333 with respect to IC 348. These groups’ ages, spatial structure, and kinematics are compatible with at least three generations of stars. Future work is still needed to clarify if the formation of the youngest was triggered by the oldest. Conclusions. The exquisite Gaia data complemented with public archives and mined with comprehensive Bayesian methodologies allow us to identify 31% more members than previous studies, discover a new physical group (Gorgophone: 7 Myr, 191 members, and 145 M⊙), and confirm that the spatial, kinematic, and energy distributions of these groups support the hierarchical star formation scenario.

NICMOS Kernel-phase Interferometry. II. Demographics of Nearby Brown Dwarfs

Star formation theories have struggled to reproduce binary brown dwarf population demographics (e.g., frequency, separation, and mass ratio). Kernel-phase interferometry is sensitive to companions at separations inaccessible to classical imaging, enabling tests of these theories at new physical scales below the hydrogen burning limit. We analyze the detections and sensitivity limits from our previous kernel-phase analysis of archival HST/NICMOS surveys of field brown dwarfs. After estimating physical properties of the 105 late-M to T dwarfs using Gaia distances and evolutionary models, we use a Bayesian framework to compare these results to a model companion population defined by log-normal separation and power-law mass-ratio distributions. When correcting for Malmquist bias, we find a companion fraction of and a separation distribution centered at au, smaller and tighter than seen in previous studies. We also find a mass-ratio power-law index that strongly favors equal-mass systems: depending on the assumed age of the field population (0.9–3.1 Gyr). We attribute the change in values to our use of kernel-phase interferometry, which enables us to resolve the peak of the semimajor axis distribution with significant sensitivity to low-mass companions. We confirm the previously seen trends of decreasing binary fraction with decreasing mass and a strong preference for tight and equal-mass systems in the field-age substellar regime; only % of systems are wider than 20 au and % of systems have a mass ratio q < 0.6. We attribute this to turbulent fragmentation setting the initial conditions followed by a brief period of dynamical evolution, removing the widest and lowest-mass companions, before the birth cluster dissolves.

Ultra-diffuse Galaxies as Extreme Star-forming Environments. II. Star Formation and Pressure Balance in H i-rich UDGs

In addition to occupying the extreme, diffuse tail of the dwarf galaxy population, ultra-diffuse galaxies (UDGs) are themselves a key laboratory in which to study star formation in extreme low-density environments. In the second paper of this series, we compare the spatially resolved star formation activity of 22 H i-selected UDGs and 21 “normal” dwarf galaxies within 120 Mpc to predictions within the pressure-regulated, feedback-modulated (PRFM) theory of star formation. To do so, we employ a joint spectral energy distribution fitting method that allows us to estimate star formation rate and stellar mass surface density from UV-optical imaging. We find that the PRFM framework extends successfully to the UDG regime—although the UDGs in our sample show unusually low star formation rate surface densities given their H i content, this low star formation efficiency can be naturally explained by the diffuse structure of the UDGs. In fact, when cast in the PRFM framework, the relationship between midplane pressure and star formation in the UDG sample is in good agreement not only with the “normal” dwarf reference sample, but also with measurements from more massive galaxies. Our results suggest that despite their low star formation efficiencies, the H i-rich UDGs need not be forming stars in an exotic manner. We also find that the UDGs are likely H2 poor compared even to the overall dwarf population.

Isolated Massive Star Formation in G28.20-0.05

We report high-resolution 1.3 mm continuum and molecular line observations of the massive protostar G28.20-0.05 with Atacama Large Millimeter/submillimeter Array. The continuum image reveals a ring-like structure with 2000 au radius, similar to morphology seen in archival 1.3 cm Very Large Array observations. Based on its spectral index and associated H30α emission, this structure mainly traces ionized gas. However, there is evidence for ∼30 M ⊙ of dusty gas near the main millimeter continuum peak on one side of the ring, as well as in adjacent regions within 3000 au. A virial analysis on scales of ∼2000 au from hot core line emission yields a dynamical mass of ∼80 M ⊙. A strong velocity gradient in the H30α emission is evidence for a rotating, ionized disk wind, which drives a larger-scale molecular outflow. An infrared spectral energy distribution (SED) analysis indicates a current protostellar mass of m * ∼ 40 M ⊙ forming from a core with initial mass M c ∼ 300 M ⊙ in a clump with mass surface density of Σcl ∼ 0.8 g cm−2. Thus the SED and other properties of the system can be understood in the context of core accretion models. A structure-finding analysis on the larger-scale continuum image indicates G28.20-0.05 is forming in a relatively isolated environment, with no other concentrated sources, i.e., protostellar cores, above ∼1 M ⊙ found from ∼0.1 to 0.4 pc around the source. This implies that a massive star can form in relative isolation, and the dearth of other protostellar companions within the ∼1 pc environs is a strong constraint on massive star formation theories that predict the presence of a surrounding protocluster.