Search for: All records

ABSTRACT In their early, formative stages star clusters can undergo rapid dynamical evolution leading to strong gravitational interactions and ejection of “runaway” stars at high velocities. While O/B runaway stars have been well studied, lower-mass runaways are so far very poorly characterized, even though they are expected to be much more common. We carried out spectroscopic observations with MAG2-MIKE to follow-up 27 high priority candidate runaways consistent with having been ejected from the Orion Nebula Cluster (ONC) $$\gt 2.5$$ Myr ago, based on Gaia astrometry. We derive spectroscopic youth indicators (Li and H $$\alpha$$) and radial velocities, enabling detection of bona fide runaway stars via signatures of youth and 3D traceback. We successfully confirmed 11 of the candidates as low-mass Young Stellar Objects (YSOs) on the basis of our spectroscopic criteria and derived radial velocities (RVs) with which we performed 3D traceback analysis. Three of these confirmed YSOs have kinematic ejection ages $$\gt 4\:$$ Myr, with the oldest being 4.7 Myr. Assuming that these stars indeed formed in the ONC and were then ejected, this yields an estimate for the overall formation time of the ONC to be at least $$\sim 5\:$$ Myr, i.e. about 10 free-fall times, and with a mean star formation efficiency per free-fall time of $$\bar{\epsilon }_{\rm ff}\lesssim 0.05$$. These results favour a scenario of slow, quasi-equilibrium star cluster formation, regulated by magnetic fields and/or protostellar outflow feedback. 

Context.Most stars form in clusters or associations, but only a small number of these groups are expected to remain bound for longer than a few megayears. Once star formation has ended and the molecular gas around young stellar objects has been expelled via feedback processes, most initially bound young clusters lose the majority of their binding mass and begin to disperse into the Galactic field. Aims.This process can be investigated by analysing the structure and kinematic trends in nearby young clusters, particularly by analyzing the trend of expansion, which is a tell-tale sign that a cluster is no longer gravitationally bound and dispersing into the field. Methods.We combinedGaiaDR3 five-parameter astrometry with calibrated RVs for members of the nearby young clusterλOri (Collinder 69). Results.We characterised the plane-of-sky substructure of the cluster using theQ-parameter and the angular dispersion parameter. We find evidence that the cluster contains a significant substructure but that this is preferentially located away from the central cluster core, which is smooth and likely remains bound. We found strong evidence for expansion inλOri in the plane of sky by using a number of metrics, but we also found that the trends are asymmetric at the 5σsignificance level, with the maximum rate of expansion being directed nearly parallel to the Galactic plane. We subsequently inverted the maximum rate of expansion of 0.144_−0.003^+0.003kms⁻¹pc⁻¹to give an expansion timescale of 6.944_−0.142^+0.148Myr, which is slightly larger than the typical literature age estimates for the cluster. We also found asymmetry in the velocity dispersion as well as signatures of cluster rotation, and we calculated the kinematic ages for individual cluster members by tracing their motion back in time to their closest approach to the cluster centre. 

Context. Evidence that the chemical characteristics around low- and high-mass protostars are similar has been found: notably, a variety of carbon-chain species and complex organic molecules (COMs) form around both types. On the other hand, the chemical compositions around intermediate-mass (IM) protostars (2M_⊙<m_*< 8M_⊙) have not been studied with large samples. In particular, it is unclear the extent to which carbon-chain species form around them. Aims. We aim to obtain the chemical compositions of a sample of IM protostars, focusing particularly on carbon-chain species. We also aim to derive the rotational temperatures of HC₅N to confirm whether carbon-chain species are formed in the warm gas around these stars. Methods. We conducted Q-band (31.5–50 GHz) line survey observations toward 11 mainly IM protostars with the Yebes 40 m radio telescope. The target protostars were selected from a subsample of the source list of the SOFIA Massive Star Formation project. Assuming local thermodynamic equilibrium, we derived the column densities of the detected molecules and the rotational temperatures of HC₅N and CH₃OH. Results. Nine carbon-chain species (HC₃N, HC₅N, C₃H, C₄Hlinear-H₂CCC,cyclic-C₃H₂, CCS, C₃S, and CH₃CCH), three COMs (CH₃OH, CH₃CHO, and CH₃CN), H₂CCO, HNCO, and four simple sulfur-bearing species (¹³CS, C³⁴S, HCS⁺, and H₂CS) are detected. The rotational temperatures of HC₅N are derived to be ~20–30 K in three IM protostars (Cepheus E, HH288, and IRAS 20293+3952). The rotational temperatures of CH₃OH are derived in five IM sources and found to be similar to those of HC₅N. Conclusions. The rotational temperatures of HC₅N around the three IM protostars are very similar to those around low- and high-mass protostars. These results indicate that carbon-chain molecules are formed in lukewarm gas (~20–30 K) around IM protostars via the warm carbon-chain chemistry process. Thus, carbon-chain formation occurs ubiquitously in the warm gas around protostars across a wide range of stellar masses. Carbon-chain molecules and COMs coexist around most of the target IM protostars, which is similar to the situation for low- and high-mass protostars. In summary, the chemical characteristics around protostars are the same in the low-, intermediate- and high-mass regimes. 

Abstract The “tuning-fork” (TF) analysis of CO and Hαemission has been used to estimate the lifetimes of molecular clouds in nearby galaxies. With simple model calculations, we show that this analysis does not necessarily estimate cloud lifetimes, but instead captures a duration of the cloud evolutionary cycle, from dormant to star-forming, and then back to a dormant phase. We adopt a hypothetical setup in which molecular clouds (e.g., traced in CO) live forever and form stars (e.g., Hiiregions) at some frequency, which then drift away from the clouds. The TF analysis still returns a timescale for the immortal clouds. This model requires drifting motion to separate the newborn stars from the clouds, and we discuss its origin. We also discuss the physical origin of the characteristic spatial separation term in the TF analysis and a bias due to systematic error in the determination of the reference timescale. 

Abstract Magnetic fields may play a crucial role in setting the initial conditions of massive star and star cluster formation. To investigate this, we report SOFIA-HAWC+ 214μm observations of polarized thermal dust emission and high-resolution GBT-Argus C¹⁸O(1-0) observations toward the massive Infrared Dark Cloud (IRDC) G28.37+0.07. Considering the local dispersion ofB-field orientations, we produce a map of the B-field strength of the IRDC, which exhibits values between ∼0.03 and 1 mG based on a refined Davis–Chandrasekhar–Fermi method proposed by Skalidis & Tassis. Comparing to a map of inferred density, the IRDC exhibits aB–nrelation with a power-law index of 0.51 ± 0.02, which is consistent with a scenario of magnetically regulated anisotropic collapse. Consideration of the mass-to-flux ratio map indicates that magnetic fields are dynamically important in most regions of the IRDC. A virial analysis of a sample of massive, dense cores in the IRDC, including evaluation of magnetic and kinetic internal and surface terms, indicates consistency with virial equilibrium, sub-Alfvénic conditions, and a dominant role forB-fields in regulating collapse. A clear alignment of magnetic field morphology with the direction of the steepest column density gradient is also detected. However, there is no preferred orientation of protostellar outflow directions with theB-field. Overall, these results indicate that magnetic fields play a crucial role in regulating massive star and star cluster formation, and therefore they need to be accounted for in theoretical models of these processes. 

Abstract We study the astrochemical diagnostics of the isolated massive protostar G28.20-0.05. We analyze data from Atacama Large Millimeter/submillimeter Array 1.3 mm observations with a resolution of 0.″2 (∼1000 au). We detect emission from a wealth of species, including oxygen-bearing (e.g., H₂CO, CH₃OH, CH₃OCH₃), sulfur-bearing (SO₂, H₂S), and nitrogen-bearing (e.g., HNCO, NH₂CHO, C₂H₃CN, C₂H₅CN) molecules. We discuss their spatial distributions, physical conditions, correlation between different species, and possible chemical origins. In the central region near the protostar, we identify three hot molecular cores (HMCs). HMC1 is part of a millimeter continuum ring-like structure, is closest in projection to the protostar, has the highest temperature of ∼300 K, and shows the most line-rich spectra. HMC2 is on the other side of the ring, has a temperature of ∼250 K, and is of intermediate chemical complexity. HMC3 is further away, ∼3000 au in projection, cooler (∼70 K), and is the least line-rich. The three HMCs have similar mass surface densities (∼10 g cm⁻²), number densities (n_H∼ 10⁹cm⁻³), and masses of a few solar masses. The total gas mass in the cores and in the region out to 3000 au is ∼25M_⊙, which is comparable to that of the central protostar. Based on spatial distributions of peak line intensities as a function of excitation energy, we infer that the HMCs are externally heated by the protostar. We estimate column densities and abundances of the detected species and discuss the implications for hot core astrochemistry. 

ABSTRACT Compression in giant molecular cloud (GMC) collisions is a promising mechanism to trigger the formation of massive star clusters and OB associations. We simulate colliding and non-colliding magnetized GMCs and examine the properties of pre-stellar cores, selected from projected mass surface density maps, including after synthetic ALMA observations. We then examine core properties, including mass, size, density, velocity, velocity dispersion, temperature, and magnetic field strength. After 4 Myr, ∼1000 cores have formed in the GMC collision, and the high-mass end of the core mass function (CMF) can be fit by a power-law dN/dlogM ∝ M−α with α ≃ 0.7, i.e. relatively top heavy compared to a Salpeter mass function. Depending on how cores are identified, a break in the power law can appear around a few $$\times 10\, \mathrm{M}_\odot$$. The non-colliding GMCs form fewer cores with a CMF with α ≃ 0.8–1.2, i.e. closer to the Salpeter index. We compare the properties of these CMFs to those of several observed samples of cores. Considering other properties, cores formed from colliding clouds are typically warmer, have more disturbed internal kinematics, and are more likely to be gravitational unbound, than cores formed from non-colliding GMCs. The dynamical state of the protocluster of cores formed in the GMC–GMC collision is intrinsically subvirial but can appear to be supervirial if the total mass measurement is affected by observations that miss mass on large scales or at low densities. 

Abstract We adopt the deep learning methodcasi-3d(convolutional approach to structure identification-3D) to infer the orientation of magnetic fields in sub-/trans-Alfvénic turbulent clouds from molecular line emission. We carry out magnetohydrodynamic simulations with different magnetic field strengths and use these to generate synthetic observations. We apply the 3D radiation transfer coderadmc-3dto model¹²CO and¹³CO (J = 1−0) line emission from the simulated clouds and then train acasi-3dmodel on these line emission data cubes to predict magnetic field morphology at the pixel level. The trainedcasi-3dmodel is able to infer magnetic field directions with a low error (≲10° for sub-Alfvénic samples and ≲30° for trans-Alfvénic samples). We further test the performance ofcasi-3don a real sub-/trans- Alfvénic region in Taurus. Thecasi-3dprediction is consistent with the magnetic field direction inferred from Planck dust polarization measurements. We use our developed methods to produce a new magnetic field map of Taurus that has a three times higher angular resolution than the Planck map. 

Context.Traditionally, supersonic turbulence is considered to be one of the most likely mechanisms slowing the gravitational collapse in dense clumps, thereby enabling the formation of massive stars. However, several recent studies have raised differing points of view based on observations carried out with sufficiently high spatial and spectral resolution. These studies call for a re-evaluation of the role turbulence plays in massive star-forming regions. Aims.Our aim is to study the gas properties, especially the turbulence, in a sample of massive star-forming regions with sufficient spatial and spectral resolution, which can both resolve the core fragmentation and the thermal line width. Methods.We observed NH₃metastable lines with the Very Large Array (VLA) to assess the intrinsic turbulence. Results.Analysis of the turbulence distribution histogram for 32 identified NH₃cores reveals the presence of three distinct components. Furthermore, our results suggest that (1) sub- and transonic turbulence is a prevalent (21 of 32) feature of massive star-forming regions and those cold regions are at early evolutionary stage. This investigation indicates that turbulence alone is insufficient to provide the necessary internal pressure required for massive star formation, necessitating further exploration of alternative candidates; and (2) studies of seven multi-core systems indicate that the cores within each system mainly share similar gas properties and masses. However, two of the systems are characterized by the presence of exceptionally cold and dense cores that are situated at the spatial center of each system. Our findings support the hub-filament model as an explanation for this observed distribution. 

Abstract We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high- and low-density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations.