Search for: All records

Abstract Star formation is a fundamental, yet poorly understood, process of the Universe. It is important to study how star formation occurs in different galactic environments. Thus, here, in the first of a series of papers, we introduce the Low-metallicity Star Formation (LZ-STAR) survey of the Sh2-284 (hereafter S284) region, which, atZ ∼  0.3–0.5Z_⊙, is one of the lowest-metallicity star-forming regions of our Galaxy. LZ-STAR is a multifacility survey, including observations with JWST, the Atacama Large Millimeter/submillimeter Array (ALMA), Hubble Space Telescope, Chandra, and Gemini. As a starting point, we report JWST and ALMA observations of one of the most massive protostars in the region, S284p1. The observations of shock-excited molecular hydrogen reveal a symmetric, bipolar outflow originating from the protostar, spanning several parsecs, and fully covered by the JWST field of view and ALMA observations of CO(2–1) emission. These allow us to infer that the protostar has maintained a relatively stable orientation of disk accretion over its formation history. The JWST near-infrared continuum observations detect a centrally illuminated bipolar outflow cavity around the protostar, as well as a surrounding cluster of low-mass young stars. We develop new radiative transfer models of massive protostars designed for the low metallicity of S284. Fitting these models to the protostar’s spectral energy distribution implies a current protostellar mass of ∼10M_⊙has formed from an initial ∼100M_⊙core over the last ∼3 × 10⁵yr. Overall, these results indicate that massive stars can form in an ordered manner in low-metallicity, protocluster environments. 

Abstract We present ∼8–40μm SOFIA-FORCAST images of seven regions of “clustered” star formation as part of the SOFIA Massive Star Formation Survey. We identify a total of 34 protostar candidates and build their spectral energy distributions (SEDs). We fit these SEDs with a grid of radiative transfer models based on the turbulent core accretion (TCA) theory to derive key protostellar properties, including initial core mass,M_c, clump environment mass surface density, Σ_cl, and current protostellar mass,m_*. We also carry out empirical graybody (GB) estimation of Σ_cl, which allows a case of restricted SED fitting within the TCA model grid. We also release version 2.0 of the open-source Python packagesedcreator, which is designed to automate the aperture photometry and SED building and fitting process for sources in clustered environments, where flux contamination from close neighbors typically complicates the process. Using these updated methods, SED fitting yields values ofM_c∼ 30–200M_⊙, Σ_cl,SED∼ 0.1–3 g cm⁻², andm_*∼ 4–50M_⊙. The GB fitting yields smaller values of Σ_cl,GB≲ 1 g cm⁻². From these results, we do not find evidence for a critical Σ_clneeded to form massive (≳8M_⊙) stars. However, we do find tentative evidence for a dearth of the most massive (m_*≳ 30M_⊙) protostars in the clustered regions, suggesting a potential impact of environment on the stellar initial mass function. 

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 We present ∼10–40μm SOFIA-FORCAST images of 11isolatedprotostars as part of the SOFIA Massive (SOMA) Star Formation Survey, with this morphological classification based on 37μm imaging. We develop an automated method to define source aperture size using the gradient of its background-subtracted enclosed flux and apply this to build spectral energy distributions (SEDs). We fit the SEDs with radiative transfer models, developed within the framework of turbulent core accretion (TCA) theory, to estimate key protostellar properties. Here, we release the sedcreator python package that carries out these methods. The SEDs are generally well fitted by the TCA models, from which we infer initial core massesM_cranging from 20–430M_⊙, clump mass surface densities Σ_cl∼ 0.3–1.7 g cm⁻², and current protostellar massesm_*∼ 3–50M_⊙. From a uniform analysis of the 40 sources in the full SOMA survey to date, we find that massive protostars form across a wide range of clump mass surface density environments, placing constraints on theories that predict a minimum threshold Σ_clfor massive star formation. However, the upper end of them_*−Σ_cldistribution follows trends predicted by models of internal protostellar feedback that find greater star formation efficiency in higher Σ_clconditions. We also investigate protostellar far-IR variability by comparison with IRAS data, finding no significant variation over an ∼40 yr baseline. 

Abstract 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 ∼30M_⊙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 ∼80M_⊙. 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 ofm_*∼ 40M_⊙forming from a core with initial massM_c∼ 300M_⊙in a clump with mass surface density of Σ_cl∼ 0.8 g cm⁻². 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 ∼1M_⊙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. 

Abstract We present a detailed study of the massive star-forming region G35.2-0.74N with Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm multi-configuration observations. At 0.″2 (440 au) resolution, the continuum emission reveals several dense cores along a filamentary structure, consistent with previous ALMA 0.85 mm observations. At 0.″03 (66 au) resolution, we detect 22 compact sources, most of which are associated with the filament. Four of the sources are associated with compact centimeter continuum emission, and two of these are associated with H30αrecombination line emission. The H30αline kinematics shows the ordered motion of the ionized gas, consistent with disk rotation and/or outflow expansion. We construct models of photoionized regions to simultaneously fit the multiwavelength free–free fluxes and the H30αtotal fluxes. The derived properties suggest the presence of at least three massive young stars with nascent hypercompact Hiiregions. Two of these ionized regions are surrounded by a large rotating structure that feeds two individual disks, revealed by dense gas tracers, such as SO₂, H₂CO, and CH₃OH. In particular, the SO₂emission highlights two spiral structures in one of the disks and probes the faster-rotating inner disks. The¹²CO emission from the general region reveals a complex outflow structure, with at least four outflows identified. The remaining 18 compact sources are expected to be associated with lower-mass protostars forming in the vicinity of the massive stars. We find potential evidence for disk disruption due to dynamic interactions in the inner region of this protocluster. The spatial distribution of the sources suggests a smooth overall radial density gradient without subclustering, but with tentative evidence of primordial mass segregation. 

Context. Hydrogen recombination lines (HRLs) are valuable diagnostics of the physical conditions in ionized regions surrounding high-mass stars. Understanding these lines, including broadening mechanisms and intensity trends, can provide insights into HII region densities, temperatures, and kinematics. Aims. This study aims to investigate the physical properties of ionized gas around massive protostars by analysing the HRLs (Hαand Hβ) in the Q band. Methods. We carried out observations using the Yebes 40m radio telescope in the Q band (30.5–50 GHz) towards six high-mass protostars selected from the SOMA Survey (G45.12+0.13, G45.47+0.05, G28.20−0.05, G35.20−0.74, G19.08−0.29, and G31.28+0.06). The observed line profiles were analysed to assess broadening mechanisms, and electron densities and temperatures were derived. The results were compared with available Q-band data from the TianMa 65-m Radio Telescope (TMRT) that have been reported in the literature, and ALMA Band 1 (35–50 GHz) Science Verification observations towards Orion KL, analysed in this study. Results. A total of eight Hα(n = 51 to 58) and ten Hβ(n = 64 to 73) lines were detected towards G45.12+0.13, G45.47+0.05, and G28.20−0.05; there were no detections in other sources. We derived electron densities of ~1−5 × 10⁶cm⁻³and temperatures of 8000–10 000 K for the sources. However, for Orion KL, we obtained an electron density one order of magnitude lower, while its temperature was found to be more similar. Interestingly, G45.12 and G28.20 show an increasing intensity trend with frequency for both Hαand Hβtransitions, contrary to the decreasing trend observed in Orion KL. Conclusions. The line widths of the detected HRLs indicate contributions from both thermal and dynamical broadening, suggesting the presence of high-temperature ionized gas that is likely kinematically broadened (e.g. due to turbulence, outflows, rapid rotation, or stellar winds). Pressure broadening caused by electron density may also have a minor effect. We discuss different scenarios to explain the measured line widths of the HRLs. The contrasting intensity trends between the sources may reflect variations in local physical conditions or radiative transfer effects, highlighting the need for further investigation through higher-resolution observations and detailed modelling. 

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. 

Context.Massive stars have an impact on their surroundings from early in their formation until the end of their lives. However, very little is known about their formation. Episodic accretion may play a crucial role in the process, but only a handful of observations have reported such events occurring in massive protostars. Aims.We aim to investigate the outburst event from the high-mass star-forming region S255IR where the protostar NIRS3 recently underwent an accretion outburst. We follow the evolution of this source both in photometry and morphology of its surroundings. Methods.We performed near infrared adaptive optics observations on the S255IR central region using the Large Binocular Telescope in theK_sbroadband as well as the H₂and Brγ narrow-band filters with an angular resolution of ~07″.06, close to the diffraction limit. Results.We discovered a new near infrared knot north-east of NIRS3 that we interpret as a jet knot that was ejected during the last accretion outburst and observed in the radio regime as part of a follow-up after the outburst. We measured a mean tangential velocity for this knot of 450 ± 50 km s⁻¹. We analysed the continuum-subtracted images from H₂, which traces jet-shocked emission, and Brγ, which traces scattered light from a combination of accretion activity and UV radiation from the central massive protostar. We observed a significant decrease in flux at the location of NIRS3, withK= 13.48 mag being the absolute minimum in the historic series. Conclusions.Our observations strongly suggest a scenario where the episodic accretion is followed by an episodic ejection response in the near infrared, as was seen in the earlier radio follow-up. The ~2 µm photometry from the past 30 yr suggests that NIRS3 might have undergone another outburst in the late 1980s, making it the first massive protostar with such evidence observed in the near infrared. 

Abstract We present Atacama Large Millimeter/submillimeter Array Band 3 data toward five massive young stellar objects (MYSOs), and investigate relationships between unsaturated carbon-chain species and saturated complex organic molecules (COMs). An HC 5 N ( J = 35–34) line has been detected from three MYSOs, where nitrogen (N)-bearing COMs (CH 2 CHCN and CH 3 CH 2 CN) have been detected. The HC 5 N spatial distributions show compact features and match with a methanol (CH 3 OH) line with an upper-state energy around 300 K, which should trace hot cores. The hot regions are more extended around the MYSOs where N-bearing COMs and HC 5 N have been detected compared to two MYSOs without these molecular lines, while there are no clear differences in the bolometric luminosity and temperature. We run chemical simulations of hot-core models with a warm-up stage, and compare with the observational results. The observed abundances of HC 5 N and COMs show good agreements with the model at the hot-core stage with temperatures above 160 K. These results indicate that carbon-chain chemistry around the MYSOs cannot be reproduced by warm carbon-chain chemistry, and a new type of carbon-chain chemistry occurs in hot regions around MYSOs.