Search

Evidence for similar chemical characteristics around low- and high-mass protostars has been found: in particular, a variety of carbon-chain species and complex organic molecules (COMs) are formed around them. On the other hand, the chemical compositions around intermediate-mass (IM; $2 M_{\odot} < m_* <8 M_{\odot}$) protostars have not been studied with large samples. In particular, it is unclear the extent to which carbon-chain species are formed around them. We aim to obtain the chemical compositions, particularly focusing on carbon-chain species, towards a sample of IM protostars. We have conducted Q-band (31.5-50 GHz) line survey observations towards eleven mainly intermediate-mass protostars with the Yebes 40 m radio telescope. The target protostars were selected from a sub-sample of the source list of the SOFIA Massive (SOMA) Star Formation project. Nine carbon-chain species (HC$_3$N, HC$_5$N, C$_3$H, C$_4$H, $linear-$H$_2$CCC, $cyclic-$C$_3$H$_2$, CCS, C$_3$S, and CH$_3$CCH), three COMs (CH$_3$OH, CH$_3$CHO, and CH$_3$CN), H$_2$CCO, HNCO, and four simple sulfur (S)-bearing species ($^{13}$CS, C$^{34}$S, HCS$^+$, H$_2$CS) have been detected. The rotational temperatures of HC$_5$N are derived to be $\sim20-30$ K in three IM protostars and they are very similar compared to those around low- and high-mass protostars. These results indicate that carbon-chain molecules are formed in lukewarm ($\sim20-30$ K) gas around the IM protostars by the Warm Carbon-Chain Chemistry (WCCC) process. 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 in low- and high-mass protostars. The chemical characteristics around protostars are common in the low-, intermediate- and high-mass regimes., Comment: Accepted for publication in the Astronomy and Astrophysics (A&A)

Filamentary structures are ubiquitously found in high-mass star-forming clouds. To investigate the relationship between filaments and star formation, we carry out the INFANT (INvestigations of massive Filaments ANd sTar formation) survey, a multi-scale, multi-wavelength survey of massive filamentary clouds with ALMA band 3/band 6 and VLA K band. In this first paper, we present the ALMA band 6 continuum observations toward a sample of 8 high-mass star forming filaments. We covered each target with approximately rectangular mosaic field of view with two 12-m array configurations, achieving an angular resolution of $\sim$0.6" (2700 AU at 4.5 kpc) and a continuum rms of $\sim$0.1 mJy/beam ($\sim$0.06 Msun in gas mass assuming 15 K). We identify cores using the getsf and astrodendro and find the former is more robust in terms of both identification and measuring flux densities. We identify in total 183 dense cores (15--36 cores in each cloud) and classify their star formation states via outflow and warm gas tracers. The protostellar cores are statistically more massive than the prestellar cores, possibly indicating further accretion onto cores after formation of protostars. For the high-mass end ($M_\text{core}$ $>$ 1.5 Msun) of the core mass function (CMF) we derive a power-law index of $-$1.15 $\pm$ 0.12 for the whole sample, and $-$1.70 $\pm$ 0.25 for the prestellar population. We also find a steepening trend in CMF with cloud evolution ($-$0.89 $\pm$ 0.15 for the young group v.s. $-$1.44 $\pm$ 0.25 for the evolved group) and discuss its implication for cluster formation., Comment: 25 pages, 8 figures, accepted for ApJ

Massive protostars launch accretion-powered, magnetically-collimated outflows, which play crucial roles in the dynamics and diagnostics of the star formation process. Here we calculate the shock heating and resulting free-free radio emission in numerical models of outflows of massive star formation within the framework of the Turbulent Core Accretion model. We post-process 3D magneto-hydrodynamic simulation snapshots of a protostellar disk wind interacts with an infalling core envelope, and calculate shock temperatures, ionization fractions, and radio free-free emission. We find heating up to ~10^7 K and near complete ionization in shocks at the interface between the outflow cavity and infalling envelope. However, line-of-sight averaged ionization fractions peak around ~10%, in agreement with values reported from observations of massive protostar G35.20-0.74N. By calculating radio continuum fluxes and spectra, we compare our models with observed samples of massive protostars. We find our fiducial models produce radio luminosities similar to those seen from low and intermediate-mass protostars that are thought to be powered by shock ionization. Comparing to more massive protostars, we find our model radio luminosities are ~10 to 100 times less luminous. We discuss how this apparent discrepancy either reflects aspects of our modeling related to the treatment of cooling of the post-shock gas or a dominant contribution in the observed systems from photoionization. Finally, our models exhibit 10-year radio flux variability of ~5%, especially in the inner 1000 au region, comparable to observed levels in some hyper-compact HII regions., Comment: Submitted to ApJ

To test theoretical models of massive star formation it is important to compare their predictions with observed systems. To this end, we conduct CO molecular line radiative transfer post-processing of 3D magneto-hydrodynamic (MHD) simulations of various stages in the evolutionary sequence of a massive protostellar core, including its infall envelope and disk wind outflow. Synthetic position-position-velocity (PPV) cubes of various transitions of CO, 13CO, and C18O emission are generated. We also carry out simulated Atacama Large Millimeter/submillimeter Array (ALMA) observations of this emission. We compare the mass, momentum and kinetic energy estimates obtained from molecular lines to the true values, finding that the mass and momentum estimates can have uncertainties of up to a factor of four. However, the kinetic energy estimated from molecular lines is more significantly underestimated. Additionally, we compare the mass outflow rate and momentum outflow rate obtained from the synthetic spectra with the true values. Finally, we compare the synthetic spectra with real examples of ALMA-observed protostars and determine the best fitting protostellar masses and outflow inclination angles. We then calculate the mass outflow rate and momentum outflow rate for these sources, finding that both rates agree with theoretical protostellar evolutionary tracks., Comment: ApJ Accepted

We report the first extragalactic detection of the higher-order millimeter hydrogen recombination lines ($\Delta n>2$). The $\gamma$-, $\epsilon$-, and $\eta$-transitions have been detected toward the millimeter continuum source N105-1A in the star-forming region N105 in the Large Magellanic Cloud (LMC) with the Atacama Large Millimeter/submillimeter Array (ALMA). We use the H40$\alpha$ line, the brightest of the detected recombination lines (H40$\alpha$, H36$\beta$, H50$\beta$, H41$\gamma$, H57$\gamma$, H49$\epsilon$, H53$\eta$, and H54$\eta$), and/or the 3 mm free-free continuum emission to determine the physical parameters of N105-1A (the electron temperature, emission measure, electron density, and size) and study ionized gas kinematics. We compare the physical properties of N105-1A to a large sample of Galactic compact and ultracompact (UC) H II regions and conclude that N105-1A is similar to the most luminous ($L>10^5$ $L_{\odot}$) UC H II regions in the Galaxy. N105-1A is ionized by an O5.5 V star, it is deeply embedded in its natal molecular clump, and likely associated with a (proto)cluster. We incorporate high-resolution molecular line data including CS, SO, SO$_2$, and CH$_3$OH ($\sim$0.12 pc), and HCO$^{+}$ and CO ($\sim$0.087 pc) to explore the molecular environment of N105-1A. Based on the CO data, we find evidence for a cloud-cloud collision that likely triggered star formation in the region. We find no clear outflow signatures, but the presence of filaments and streamers indicates on-going accretion onto the clump hosting the UC H II region. Sulfur chemistry in N105-1A is consistent with the accretion shock model predictions., Comment: 51 pages, 30 figures, 2 tables (including appendices); accepted for publication in The Astrophysical Journal (ApJ)

The destiny of complex organic molecules (COMs) in star-forming regions is interlinked with various evolutionary phases. Therefore, identifying these species in diversified environments of identical star-forming regions would help to comprehend their physical and chemical heritage. We identified multiple COMs utilizing the Large Program `Astrochemical Surveys At IRAM' (ASAI) data, dedicated to chemical surveys in Sun-like star-forming regions with the IRAM 30 m telescope. It was an unbiased survey in the millimetre regime, covering the prestellar core, protostar, outflow region, and protoplanetary disk phase. Here, we have reported some transitions of seven COMs, namely, methanol (CH3OH), acetaldehyde (CH3CHO), methyl formate (CH3OCHO), ethanol (C2H5OH), propynal (HCCCHO), dimethyl ether (CH3OCH3), and methyl cyanide (CH3CN) in some sources L1544, B1-b, IRAS4A, and SVS13A. We found a trend among these species from the derived abundances using the rotational diagram method and MCMC fit. We have found that the abundances of all of the COMs, except for HCCCHO, increase from the L1544 (prestellar core) and peaks at IRAS16293-2422 (class 0 phase). It is noticed that the abundance of these molecules correlate with the luminosity of the sources. The obtained trend is also visible from the previous interferometric observations and considering the beam dilution effect., Comment: 44 pages, 25 figures, and 12 tables. Accepted for the publication in the Astrophysical Journal

Recent millimeter/sub-millimeter facilities have revealed the physical properties of filamentary molecular clouds in relation to high-mass star formation. A uniform survey of the nearest, face-on star-forming galaxy, the Large Magellanic Cloud (LMC), complements the Galactic knowledge. We present ALMA survey data with a spatial resolution of $\sim$0.1 pc in the 0.87 mm continuum and HCO$^{+}$(4-3) emission toward 30 protostellar objects with luminosities of 10$^4$-10$^{5.5}$ $L_{\odot}$ in the LMC. The spatial distributions of the HCO$^{+}$(4-3) line and thermal dust emission are well correlated, indicating that the line effectively traces dense, filamentary gas with an H$_2$ volume density of $\gtrsim$10$^5$ cm$^{-3}$ and a line mass of $\sim$10$^3$-10$^{4}$ $M_{\odot}$ pc$^{-1}$. Furthermore, we obtain an increase in the velocity linewidths of filamentary clouds, which follows a power-law dependence on their H$_2$ column densities with an exponent of $\sim$0.5. This trend is consistent with observations toward filamentary clouds in nearby star-forming regions withiin $ \lesssim$1 kpc from us and suggests enhanced internal turbulence within the filaments owing to surrounding gas accretion. Among the 30 sources, we find that 14 are associated with hub-filamentary structures, and these complex structures predominantly appear in protostellar luminosities exceeding $\sim$5 $\times$10$^4$ $L_{\odot}$. The hub-filament systems tend to appear in the latest stages of their natal cloud evolution, often linked to prominent H$\;${\sc ii} regions and numerous stellar clusters. Our preliminary statistics suggest that the massive filaments accompanied by hub-type complex features may be a necessary intermediate product in forming extremely luminous high-mass stellar systems capable of ultimately dispersing the parent cloud., Comment: 21 pages, 8 figures, 4 tables, accepted for publication in ApJ

Young massive clusters (YMCs) are the most massive star clusters forming in nearby galaxies and are thought to be a young analogue to the globular clusters. Understanding the formation process of YMCs leads to looking into very efficient star formation in high-redshift galaxies suggested by recent JWST observations. We investigate possible observational signatures of their formation stage, particularly when the mass of a cluster is increasing via accretion from a natal molecular cloud. To this end, we study the broad-band continuum emission from ionized gas and dust enshrouding YMCs, whose formation is followed by recent radiation-hydrodynamics simulations. We perform post-process radiative transfer calculations using simulation snapshots and find characteristic spectral features at radio and far-infrared frequencies. We show that a striking feature is long-lasting, strong free-free emission from a $\sim$ 10pc-scale HII region with a large emission measure of $\gtrsim 10^7 \mathrm{cm}^{-6} \ \mathrm{pc}$, corresponding to the mean electron density of $\gtrsim 10^3~\mathrm{cm}^{-3}$. There is a turnover feature below $\sim$ 10 GHz, a signature of the optically-thick free-free emission, often found in Galactic ultra-compact HII regions. These features come from the peculiar YMC formation process, where the cluster's gravity effectively traps photoionized gas for a long duration and enables continuous star formation within the cluster. Such large and dense HII regions show distinct distribution on the density-size diagram, apart from the standard sequence of Galactic HII regions. This is consistent with the observational trend inferred for extragalactic HII regions associated with YMCs., Comment: 12 pages, 10 figures, accepted for publication in MNRAS

We report the first detection of hot molecular cores in the Small Magellanic Cloud, a nearby dwarf galaxy with 0.2 solar metallicity. We observed two high-mass young stellar objects in the SMC with ALMA, and detected emission lines of CO, HCO+, H13CO+, SiO, H2CO, CH3OH, SO, and SO2. Compact hot-core regions are traced by SO2, whose spatial extent is about 0.1 pc, and the gas temperature is higher than 100 K based on the rotation diagram analysis. In contrast, CH3OH, a classical hot-core tracer, is dominated by extended (0.2-0.3 pc) components in both sources, and the gas temperature is estimated to be 39+-8 K for one source. Protostellar outflows are also detected from both sources as high-velocity components of CO. The metallicity-scaled abundances of SO2 in hot cores are comparable among the SMC, LMC, and Galactic sources, suggesting that the chemical reactions leading to SO2 formation would be regulated by elemental abundances. On the other hand, CH3OH shows a large abundance variation within SMC and LMC hot cores. The diversity in the initial condition of star formation (e.g., degree of shielding, local radiation field strength) may lead to the large abundance variation of organic molecules in hot cores. This work, in conjunction with previous hot-core studies in the LMC and outer/inner Galaxy, suggests that the formation of a hot core would be a common phenomenon during high-mass star formation across the metallicity range of 0.2-1 solar metallicity. High-excitation SO2 lines will be a useful hot-core tracer in the low-metallicity environments of the SMC and LMC., Comment: Accepted for publication in ApJL, 17 pages, 8 figures, 4 tables. arXiv admin note: text overlap with arXiv:2109.11123

We study the astrochemical diagnostics of the isolated massive protostar G28.20-0.05. We analyze data from ALMA 1.3~mm observations with resolution of 0.2 arcsec ($\sim$1,000 au). We detect emission from a wealth of species, including oxygen-bearing (e.g., $\rm{H_2CO}$, $\rm{CH_3OH}$, $\rm{CH_3OCH_3}$), sulfur-bearing (SO$_2$, H$_2$S) and nitrogen-bearing (e.g., HNCO, NH$_2$CHO, C$_2$H$_3$CN, C$_2$H$_5$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 mm continuum ring-like structure, is closest in projection to the protostar, has the highest temperature of $\sim300\:$K, and shows the most line-rich spectra. HMC2 is on the other side of the ring, has a temperature of $\sim250\:$K, and is of intermediate chemical complexity. HMC3 is further away, $\sim3,000\:$au in projection, cooler ($\sim70\:$K) and is the least line-rich. The three HMCs have similar mass surface densities ($\sim10\:{\rm{g\:cm}}^{-2}$), number densities ($n_{\rm H}\sim10^9\:{\rm{cm}}^{-3}$) and masses of a few $M_\odot$. The total gas mass in the cores and in the region out to $3,000\:$au is $\sim 25\:M_\odot$, 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., Comment: Accepted for publication in ApJ

Star formation is ubiquitously associated with the ejection of accretion-powered outflows that carve bipolar cavities through the infalling envelope. This feedback is expected to be important for regulating the efficiency of star formation from a natal pre-stellar core. These low-extinction outflow cavities greatly affect the appearance of a protostar by allowing the escape of shorter wavelength photons. Doppler-shifted CO line emission from outflows is also often the most prominent manifestation of deeply embedded early-stage star formation. Here, we present 3D magneto-hydrodynamic simulations of a disk wind outflow from a protostar forming from an initially $60\:M_\odot$ core embedded in a high pressure environment typical of massive star-forming regions. We simulate the growth of the protostar from $m_*=1\:M_\odot$ to $26\:M_\odot$ over a period of $\sim$100,000 years. The outflow quickly excavates a cavity with half opening angle of $\sim10^\circ$ through the core. This angle remains relatively constant until the star reaches $4\:M_\odot$. It then grows steadily in time, reaching a value of $\sim 50^\circ$ by the end of the simulation. We estimate a lower limit to the star formation efficiency (SFE) of 0.43. However, accounting for continued accretion from a massive disk and residual infall envelope, we estimate that the final SFE may be as high as $\sim0.7$. We examine observable properties of the outflow, especially the evolution of the cavity opening angle, total mass and momentum flux, and velocity distributions of the outflowing gas, and compare with the massive protostars G35.20-0.74N and G339.88-1.26 observed by ALMA, yielding constraints on their intrinsic properties., Comment: Accepted for publication in ApJ

Molecular lines tracing the orbital motion of gas in a well-defined disk are valuable tools for inferring both the properties of the disk and the star it surrounds. Lines that arise only from a disk, and not also from the surrounding molecular cloud core that birthed the star or from the outflow it drives, are rare. Several such emission lines have recently been discovered in one example case, those from NaCl and KCl salt molecules. We studied a sample of 23 candidate high-mass young stellar objects (HMYSOs) in 17 high-mass star-forming regions to determine how frequently emission from these species is detected. We present five new detections of water, NaCl, KCl, PN, and SiS from the innermost regions around the objects, bringing the total number of known briny disk candidates to nine. Their kinematic structure is generally disk-like, though we are unable to determine whether they arise from a disk or outflow in the sources with new detections. We demonstrate that these species are spatially coincident in a few resolved cases and show that they are generally detected together, suggesting a common origin or excitation mechanism. We also show that several disks around HMYSOs clearly do not exhibit emission in these species. Salty disks are therefore neither particularly rare in high-mass disks, nor are they ubiquitous., Comment: accepted to ApJ

We present a detailed study of the massive star-forming region G35.2-0.74N with 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{\alpha} recombination line emission. The H30{\alpha} line kinematics show ordered motion of the ionized gas, consistent with disk rotation and/or outflow expansion. We construct models of photoionized regions to simultaneously fit the multi-wavelength free-free fluxes and the H30{\alpha} total fluxes. The derived properties suggest the presence of at least three massive young stars with nascent hypercompact Hii regions. Two of these ionized regions are surrounded by a large rotating structure that feeds two individual disks, revealed by dense gas tracers, such as SO2, H2CO, and CH3OH. In particular, the SO2 emission highlights two spiral structures in one of the disks and probes the faster-rotating inner disks. The 12CO 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 dynamical 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., Comment: 46 pages, 21 figures, 5 tables. Accepted to ApJ

Protostellar outflows are one of the most outstanding features of star formation. Observational studies over the last several decades have successfully demonstrated that outflows are ubiquitously associated with low- and high-mass protostars in the solar-metallicity Galactic conditions. However, the environmental dependence of protostellar outflow properties is still poorly understood, particularly in the low-metallicity regime. Here we report the first detection of a molecular outflow in the Small Magellanic Cloud with 0.2 $Z_{\odot}$, using Atacama Large Millimeter/submillimeter Array observations at a spatial resolution of 0.1 pc toward the massive protostar Y246. The bipolar outflow is nicely illustrated by high-velocity wings of CO(3-2) emission at $\gtrsim$15 km s$^{-1}$. The evaluated properties of the outflow (momentum, mechanical force, etc.) are consistent with those of the Galactic counterparts. Our results suggest that the molecular outflows, i.e., the guidepost of the disk accretion at the small scale, might be universally associated with protostars across the metallicity range of $\sim$0.2-1 $Z_{\odot}$., Comment: 7 pages, 2 figures, Accepted for publication in ApJL

Cosmic metallicity evolution possibly creates the diversity of star formation modes at different epochs. Gravitational fragmentation of circumstellar discs provides an important formation channel of multiple star systems, including close binaries. We here study the nature of disc fragmentation, systematically performing a suite of two-dimensional radiation-hydrodynamic simulations, in a broad range of metallicities, from the primordial to the solar values. In particular, we follow relatively long-term disc evolution over 15 kyr after the disc formation, incorporating the effect of heating by the protostellar irradiation. Our results show that the disc fragmentation occurs at all metallicities $1$--$0$ $Z_{\odot}$, yielding self-gravitating clumps. Physical properties of the clumps, such as their number and mass distributions, change with the metallicity due to different gas thermal evolution. For instance, the number of clumps is the largest for the intermediate metallicity range of $10^{-2}$--$10^{-5}$ $Z_{\odot}$, where the dust cooling is effective exclusively in a dense part of the disc and causes the fragmentation of spiral arms. The disc fragmentation is more modest for $1$--$0.1$ $Z_{\odot}$ thanks to the disc stabilization by the stellar irradiation. Such metallicity dependence agrees with the observed trend that the close binary fraction increases with decreasing metallicity in the range of $1$--$10^{-3}$ $Z_{\odot}$., Comment: 17 pages, 16 figures, submitted to MNRAS

We present $\sim10-40\,\mu$m SOFIA-FORCAST images of 11 isolated protostars as part of the SOFIA Massive (SOMA) Star Formation Survey, with this morphological classification based on 37 $\mu$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 masses $M_c$ ranging from $20-430\:M_\odot$, clump mass surface densities $\Sigma_{\rm cl}\sim0.3-1.7\:{\rm{g\:cm}}^{-2}$ and current protostellar masses $m_*\sim3-50\:M_\odot$. 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 $\Sigma_{\rm cl}$ for massive star formation. However, the upper end of the $m_*-\Sigma_{\rm cl}$ distribution follows trends predicted by models of internal protostellar feedback that find greater star formation efficiency in higher $\Sigma_{\rm cl}$ conditions. We also investigate protostellar far-IR variability by comparison with IRAS data, finding no significant variation over an $\sim$40 year baseline., Comment: Accepted for publication in ApJ. Published version including the proofs

Massive dense clumps in the Large Magellanic Cloud can be an important laboratory to explore the formation of populous clusters. We report multiscale ALMA observations of the N159W-North clump, which is the most CO-intense region in the galaxy. High-resolution CO isotope and 1.3 mm continuum observations with an angular resolution of $\sim$0."25($\sim$0.07 pc) revealed more than five protostellar sources with CO outflows within the main ridge clump. One of the thermal continuum sources, MMS-2, shows especially massive/dense nature whose total H$_2$ mass and peak column density are $\sim$10$^{4}$ $M_{\odot}$ and $\sim$10$^{24}$ cm$^{-2}$, respectively, and harbors massive ($\sim$100 $M_{\odot}$) starless core candidates identified as its internal substructures. The main ridge containing this source can be categorized as one of the most massive protocluster systems in the Local Group. The CO high-resolution observations found several distinct filamentary clouds extending southward from the star-forming spots. The CO (1-0) data set with a larger field of view reveals a conical-shaped, $\sim$30 pc long complex extending toward the northern direction. These features indicate that a large-scale gas compression event may have produced the massive star-forming complex. Based on the striking similarity between the N159W-North complex and the previously reported other two high-mass star-forming clouds in the nearby regions, we propose a $"$teardrops inflow model$"$ that explains the synchronized, extreme star formation across $>$50 pc, including one of the most massive protocluster clumps in the Local Group., Comment: 25 pages, 14 figures, 6 tables, Accepted for publication in ApJ

We have analyzed Atacama Large Millimeter/submillimeter Array Band 6 data of the hyper-compact H$_{\rm {II}}$ region G24.78+0.08 A1 (G24 HC H$_{\rm {II}}$) and report the detection of vibrationally-excited lines of HC$_{3}$N ($v_{7}=2$, $J=24-23$). The spatial distribution and kinematics of a vibrationally-excited line of HC$_{3}$N ($v_{7}=2$, $J=24-23$, $l=2e$) are found to be similar to the CH$_{3}$CN vibrationally-excited line ($v_{8}=1$), which indicates that the HC$_{3}$N emission is tracing the disk around the G24 HC H$_{\rm {II}}$ region previously identified by the CH$_{3}$CN lines. We derive the $^{13}$CH$_{3}$CN/HC$^{13}$CCN abundance ratios around G24 and compare them to the CH$_{3}$CN/HC$_{3}$N abundance ratios in disks around Herbig Ae and T Tauri stars. The $^{13}$CH$_{3}$CN/HC$^{13}$CCN ratios around G24 ($\sim 3.0-3.5$) are higher than the CH$_{3}$CN/HC$_{3}$N ratios in the other disks ($\sim 0.03-0.11$) by more than one order of magnitude. The higher CH$_{3}$CN/HC$_{3}$N ratios around G24 suggest that the thermal desorption of CH$_{3}$CN in the hot dense gas and efficient destruction of HC$_{3}$N in the region irradiated by the strong UV radiation are occurring. Our results indicate that the vibrationally-excited HC$_{3}$N lines can be used as a disk tracer of massive protostars at the HC H$_{\rm {II}}$ region stage, and the combination of these nitrile species will provide information of not only chemistry but also physical conditions of the disk structures., Comment: 21 pages, 13 figures, 5 tables, Accepted by The Astrophysical Journal

We report high-resolution 1.3~mm continuum and molecular line observations of the massive protostar G28.20-0.05 with ALMA. The continuum image reveals a ring-like structure with 2,000~au radius, similar to morphology seen in archival 1.3~cm VLA observations. Based on its spectral index and associated H$30\alpha$ emission, this structure mainly traces ionised gas. However, there is evidence for $\sim30$~M$_{\odot}$ of dusty gas near the main mm continuum peak on one side of the ring, as well as in adjacent regions within 3,000~au. A virial analysis on scales of $\sim$2,000~au from hot core line emission yields a dynamical mass of $\sim80\:M_\odot$. A strong velocity gradient in the H$30\alpha$ emission is evidence for a rotating, ionized disk wind, which drives a larger-scale molecular outflow. An infrared SED analysis indicates a current protostellar mass of $m_*\sim40\:M_\odot$ forming from a core with initial mass $M_c\sim300\:M_\odot$ in a clump with mass surface density of $\Sigma_{\rm cl}\sim 0.8\:{\rm g\:cm}^{-2}$. Thus the SED and other properties of the system can be understood in the context of core accretion models. 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 $\sim 1\:M_\odot$ found from $\sim$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 $\sim1$~pc environs is a strong constraint on massive star formation theories that predict the presence of a surrounding protocluster., Comment: Accepted to ApJ

Multi-epoch narrow-band HST images of the bipolar Hii region Sh2-106 reveal highly supersonic nebular proper motions which increase with projected distance from the massive young stellar object S106~IR, reaching over ~30 mas/year (~150 km/s at D=1.09 kpc) at a projected separation of ~1.4' (0.44 pc) from S106~IR. We propose that S106~IR experienced a $\sim10^{47}$ erg explosion ~3,500 years ago. The explosion may be the result of a major accretion burst, a recent encounter with another star, or a consequence of the interaction of a companion with the bloated photosphere of S106~IR as it grew from ~10 through ~15 Solar masses at a high accretion rate. Near-IR images reveal fingers of molecular hydrogen emission pointing away from S106~IR and an asymmetric photon-dominated region surrounding the ionized nebula. Radio continuum and Brackett-gamma emission reveal a C-shaped bend in the plasma, either indicating motion of S106~IR toward the east, or deflection of plasma toward the west by the surrounding cloud. The Hii region bends around a ~1' diameter dark bay west of S106~IR that may be shielded from direct illumination by a dense molecular clump. Herbig-Haro (HH) and Molecular Hydrogen Objects (MHOs) tracing outflows powered by stars in the Sh2-106 proto-cluster such as the Class 0 source S106 FIR are discussed., Comment: 47 pages, 24 figures

The Next Generation Very Large Array (ngVLA) has excellent capabilities to unveil various dynamical and chemical processes in massive star formation at the unexplored innermost regions. Based on the recent observations of ALMA/VLA as well as theoretical predictions, we propose several intriguing topics in massive star formation from the perspective of the ngVLA. In the disk scale of $\lesssim$ 100 au around massive protostars, dust grains are expected to be destructed/sublimated because the physical conditions of temperature, shocks, and radiation are much more intense than those in the envelopes, which are typically observed as hot cores. The high sensitivity and resolution of the ngVLA will enable us to detect the gaseous refractories released by dust destruction, e.g., SiO, NaCl, and AlO, which trace disk kinematics and give new insights into the metallic elements in star-forming regions, i.e., astromineralogy. The multi-epoch survey by the ngVLA will provide demographics of forming massive multiples with separations of $\lesssim$ 10 au with their proper motion. Combining with observations of refractory molecular lines and hydrogen recombination lines, we can reproduce the three-dimensional orbital motions of massive proto-binaries. Moreover, the 1-mas resolution of the ngVLA could possibly take the first-ever picture of the photospheric surface of an accreting protostar, if it is bloated to the au scale by the high accretion rates of mass and thermal energy., Comment: 7 pages, 4 figures, Published in the ngVLA-Japan Memo Series

We present three-dimensional astrochemical simulations and synthetic observations of magnetised, turbulent, self-gravitating molecular clouds. We explore various galactic interstellar medium environments, including cosmic-ray ionization rates in the range of $\zeta_{\rm CR}=10^{-17}$-$10^{-14}\,{\rm s}^{-1}$, far-UV intensities in the range of $G_0=1$-$10^3$ and metallicities in the range of $Z=0.1$-$2\,{\rm Z}_{\odot}$. The simulations also probe a range of densities and levels of turbulence, including cases where the gas has undergone recent compression due to cloud-cloud collisions. We examine: i) the column densities of carbon species across the cycle of CII, CI and CO, along with OI, in relation to the HI-to-H$_2$ transition; ii) the velocity-integrated emission of [CII]~$158\mu$m, [$^{13}$CII]~$158\mu$m, [CI]~$609\mu$m and $370\mu$m, [OI]~$63\mu$m and $146\mu$m, and of the first ten $^{12}$CO rotational transitions; iii) the corresponding Spectral Line Energy Distributions; iv) the usage of [CII] and [OI]~$63\mu$m to describe the dynamical state of the clouds; v) the behavior of the most commonly used ratios between transitions of CO and [CI]; and vi) the conversion factors for using CO and CI as H$_2$-gas tracers. We find that enhanced cosmic-ray energy densities enhance all aforementioned line intensities. At low metallicities, the emission of [CII] is well connected with the H$_2$ column, making it a promising new H$_2$ tracer in metal-poor environments. The conversion factors of $X_{\rm CO}$ and $X_{\rm CI}$ depend on metallicity and the cosmic-ray ionization rate, but not on FUV intensity. In the era of ALMA, SOFIA and the forthcoming CCAT-prime telescope, our results can be used to understand better the behaviour of systems in a wide range of galactic and extragalactic environments., Comment: 33 pages, 24 figures, 3 tables. Submitted to MNRAS. Comments welcome!

We systematically perform hydrodynamics simulations of 20 km s^-1 converging flows of the warm neutral medium (WNM) to calculate the formation of the cold neutral medium (CNM), especially focusing on the mean properties of the multiphase interstellar medium (ISM), such as the average shock front position and the mean density on a 10 pc scale. Our results show that the convergence in those mean properties requires 0.02 pc spatial resolution that resolves the cooling length of the thermally unstable neutral medium (UNM) to follow the dynamical condensation from the WNM to CNM. We also find that two distinct post-shock states appear in the mean properties depending on the amplitude of the upstream WNM density fluctuation (= sqrt()/rho_0). When the amplitude > 10 %, the interaction between shocks and density inhomogeneity leads to a strong driving of the post-shock turbulence of > 3 km s^-1, which dominates the energy budget in the shock-compressed layer. The turbulence prevents the dynamical condensation by cooling and the following CNM formation, and the CNM mass fraction remains as ~ 45 %. In contrast, when the amplitude <= 10 %, the shock fronts maintain an almost straight geometry and CNM formation efficiently proceeds, resulting in the CNM mass fraction of ~ 70 %. The velocity dispersion is limited to the thermal-instability mediated level of 2 - 3 km s^-1 and the layer is supported by both turbulent and thermal energy equally. We also propose an effective equation of state that models the multiphase ISM formed by the WNM converging flow as a one-phase ISM., Comment: 21 pages, 13 figures, 3 tables, submitted to the ApJ, included the referee comments

Through the magnetic braking and the launching of protostellar outflows, magnetic fields play a major role in the regulation of angular momentum in star formation, which directly impacts the formation and evolution of protoplanetary disks and binary systems. The aim of this paper is to quantify those phenomena in the presence of non-ideal magnetohydrodynamics effects, namely the Ohmic and ambipola r diffusion. We perform three-dimensional simulations of protostellar collapses varying the mass of the prestellar dense core, the thermal support (the $\alpha$ ratio) and the dust grain size-distribu tion. The mass mostly influences the magnetic braking in the pseudo-disk, while the thermal support impacts the accretion rate and hence the properties of the disk. Removing the grains smaller than 0. 1 $\mu$m in the Mathis, Rumpl, Nordsieck (MRN) distribution enhances the ambipolar diffusion coefficient. Similarly to previous studies, we find that this change in the distribution reduces the magnet ic braking with an impact on the disk. The outflow is also significantly weakened. In either case, the magnetic braking largely dominates the outflow as a process to remove the angular momentum from t he disk. Finally, we report a large ionic precursor to the outflow with velocities of several km s$^{-1}$, which may be observable., Comment: 20 pages, 16 figures, Accepted in ApJ

We report results of 0.05"-resolution observations toward the O-type proto-binary system IRAS 16547-4247 with the Atacama Large Millimeter/submillimeter Array (ALMA). We present dynamical and chemical structures of the circumbinary disk, circumstellar disks, outflows and jets, illustrated by multi-wavelength continuum and various molecular lines. In particular, we detect sodium chloride, silicon compounds, and vibrationally-excited water lines as probes of the individual protostellar disks at a scale of 100 au. These are complementary to typical hot-core molecules tracing the circumbinary structures on a 1000-au scale. The H2O line tracing inner-disks has an upper-state energy of Eu/k>3000K, indicating a high temperature of the disks. On the other hand, despite the detected transitions of NaCl, SiO, and SiS not necessarily having high upper-state energies, they are enhanced only in the vicinity of the protostars. We interpret that these molecules are the products of dust destruction, which only happens in the inner disks. This is the second detection of alkali metal halide in protostellar systems after the case of the disk of Orion Source I, and also one of few massive protostellar disks associated with high-energy transition water and silicon compounds. These new results suggest these "hot-disk" lines may be common in innermost disks around massive protostars, and have great potential for future research of massive star formation. We also tentatively find that the twin disks are counter-rotating, which might give a hint of the origin of the massive proto-binary system IRAS 16547-4247., Comment: 15 pages, 5 figures, 2 appendix figures. Published in ApJ Letters

We present $\sim10-40\,\mu$m SOFIA-FORCAST images of 14 intermediate-mass protostar candidates as part of the SOFIA Massive (SOMA) Star Formation Survey. We build spectral energy distributions (SEDs), also utilizing archival Spitzer, Herschel and IRAS data. We then fit the SEDs with radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, to estimate key protostellar properties. With the addition of these intermediate-mass sources, SOMA protostars span luminosities from $\sim10^{2}-10^{6}\:L_{\odot}$, current protostellar masses from $\sim0.5-30\:M_{\odot}$ and ambient clump mass surface densities, $\Sigma_{\rm cl}$ from $0.1-3\:{\rm{g\:cm}^{-2}}$. A wide range of evolutionary states of the individual protostars and of the protocluster environments are also probed. We have also considered about 50 protostars identified in Infrared Dark Clouds and expected to be at the earliest stages of their evolution. With this global sample, most of the evolutionary stages of high- and intermediate-mass protostars are probed. From the best fitting models, there is no evidence of a threshold value of protocluster clump mass surface density being needed to form protostars up to $\sim25\:M_\odot$. However, to form more massive protostars, there is tentative evidence that $\Sigma_{\rm{cl}}$ needs to be $\gtrsim1\:{\rm{g\,cm}}^{-2}$. We discuss how this is consistent with expectations from core accretion models that include internal feedback from the forming massive star., Comment: 40 pages, 21 figures, 4 tables, accepted to ApJ

We present the results of 0.1-pc-scale observations in 250 GHz and 350GHz towards a newly-discovered hot molecular core in a nearby low-metallicity galaxy, the Large Magellanic Cloud (LMC), with the Atacama Large Millimeter/submillimeter Array. A variety of C/N/O/Si/S-bearing molecules are detected towards the high-mass young stellar object, ST16. A rotating protostellar envelope is for the first time detected outside our Galaxy by SO2 and 34SO lines. An outflow cavity is traced by CCH and CN. The isotope abundance of sulfur in the source is estimated to be 32S/34S = 17 and 32S/33S = 53 based on SO, SO2, and CS isotopologues, suggesting that both 34S and 33S are overabundant in the LMC. Rotation diagram analyses show that the source is associated with hot gas (>100K) traced by high-excitation lines of CH3OH and SO2, as well as warm gas (~50K) traced by CH3OH, SO2, 34SO, OCS, CH3CN lines. A comparison of molecular abundances between LMC and Galactic hot cores suggests that organic molecules (e.g., CH3OH, a classical hot core tracer) show a large abundance variation in low metallicity, where the present source is classified into an organic-poor hot core. Our astrochemical simulations suggest that different grain temperature during the initial ice-forming stage would contribute to the chemical differentiation. In contrast, SO2 shows similar abundances within all the known LMC hot cores and the typical abundance roughly scales with the LMC's metallicity. Nitrogen-bearing molecules are generally less abundant in LMC hot cores, except for NO. The present results suggest that chemical compositions of hot cores do not always simply scale with the metallicity., Comment: Accepted for publication in ApJ, 30 pages, 14 figures, 10 tables

Massive protostars generate strong radiation feedback, which may help set the mass they achieve by the end of the accretion process. Studying such feedback is therefore crucial for understanding the formation of massive stars. We report the discovery of a photoionized bipolar outflow towards the massive protostar G45.47+0.05 using high-resolution observations at 1.3 mm with the Atacama Large Millimeter/Submillimeter Array (ALMA) and at 7 mm with the Karl G. Jansky Very Large Array (VLA). By modeling the free-free continuum, the ionized outflow is found to be a photoevaporation flow with an electron temperature of 10,000 K and an electron number density of ~1.5x10^7 cm^-3 at the center, launched from a disk of radius of 110 au. H30alpha hydrogen recombination line emission shows strong maser amplification, with G45 being one of very few sources to show such millimeter recombination line masers. The mass of the driving source is estimated to be 30-50 Msun based on the derived ionizing photon rate, or 30-40 Msun based on the H30alpha kinematics. The kinematics of the photoevaporated material is dominated by rotation close to the disk plane, while accelerated to outflowing motion above the disk plane. The mass loss rate of the photoevaporation outflow is estimated to be ~(2-3.5)x10^-5 Msun/yr. We also found hints of a possible jet embedded inside the wide-angle ionized outflow with non-thermal emissions. The possible co-existence of a jet and a massive photoevaporation outflow suggests that, in spite of the strong photoionization feedback, accretion is still on-going., Comment: Accepted to ApJL. 16 pages, 5 figures, 3 appendix figures

Almost all massive stars have bound stellar companions, existing in binaries or higher-order multiples. While binarity is theorized to be an essential feature of how massive stars form, essentially all information about such properties is derived from observations of already formed stars, whose orbital properties may have evolved since birth. Little is known about binarity during formation stages. Here we report high angular resolution observations of 1.3 mm continuum and H30alpha recombination line emission, which reveal a massive protobinary with apparent separation of 180 au at the center of the massive star-forming region IRAS07299-1651. From the line-of-sight velocity difference of 9.5 km/s of the two protostars, the binary is estimated to have a minimum total mass of 18 solar masses, consistent with several other metrics, and maximum period of 570 years, assuming a circular orbit. The H30alpha line from the primary protostar shows kinematics consistent with rotation along a ring of radius of 12 au. The observations indicate that disk fragmentation at several hundred au may have formed the binary, and much smaller disks are feeding the individual protostars., Comment: Published in Nature Astronomy. This is author's version. Full article is available here (https://rdcu.be/brENk). 47 pages, 10 figures, including methods and supplementary information

We present multi-wavelength images observed with SOFIA-FORCAST from $\sim$10 to 40 $\mu$m of seven high luminosity massive protostars, as part of the SOFIA Massive (SOMA) Star Formation Survey. Source morphologies at these wavelengths appear to be influenced by outflow cavities and extinction from dense gas surrounding the protostars. Using these images, we build spectral energy distributions (SEDs) of the protostars, also including archival data from Spitzer, Herschel and other facilities. Radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, are then fit to the SEDs to estimate key properties of the protostars. Considering the best five models fit to each source, the protostars have masses $m_{*} \sim 12-64 \: M_{\odot}$ accreting at rates of $\dot{m}_{*} \sim 10^{-4}-10^{-3} \: M_{\odot} \: \rm yr^{-1}$ inside cores of initial masses $M_{c} \sim 100-500 \: M_{\odot}$ embedded in clumps with mass surface densities $\Sigma_{\rm cl} \sim 0.1-3 \: \rm g \: cm^{-2}$ and span a luminosity range of $10^{4} -10^{6} \: L_{\odot}$. Compared with the first eight protostars in Paper I, the sources analyzed here are more luminous, and thus likely to be more massive protostars. They are often in a clustered environment or have a companion protostar relatively nearby. From the range of parameter space of the models, we do not see any evidence that $\Sigma_{\rm cl}$ needs to be high to form these massive stars. For most sources the RT models provide reasonable fits to the SEDs, though the cold clump material often influences the long wavelength fitting. However, for sources in very clustered environments, the model SEDs may not be such a good description of the data, indicating potential limitations of the models for these regions., Comment: 30 pages, 19 figures, Accepted for publication in ApJ

We report molecular line observations of the massive protostellar source G339.88-1.26 with the Atacama Large Millimeter/Submillimeter Array. The observations reveal a highly collimated SiO jet extending from the 1.3 mm continuum source, which connects to a slightly wider but still highly collimated CO outflow. Rotational features perpendicular to the outflow axis are detected in various molecular emissions, including SiO, SO2, H2S, CH3OH, and H2CO emissions. Based on their spatial distributions and kinematics, we find that they trace different parts of the envelope-disk system. The SiO emission traces the disk and inner envelope in addition to the jet. The CH3OH and H2CO emissions mostly trace the infalling-rotating envelope, and are enhanced around the transition region between envelope and disk, i.e., the centrifugal barrier. The SO2 and H2S emissions are enhanced around the centrifugal barrier, and also trace the outer part of the disk. Envelope kinematics are consistent with rotating-infalling motion, while those of the disk are consistent with Keplerian rotation. The radius and velocity of the centrifugal barrier are estimated to be about 530 au and 6 km/s, leading to a central mass of about 11 solar masses, consistent with estimates based on spectral energy distribution fitting. These results indicate that an ordered transition from an infalling-rotating envelope to a Keplerian disk through a centrifugal barrier, accompanied by changes of types of molecular line emissions, is a valid description of this massive protostellar source. This implies that at least some massive stars form in a similar way as low-mass stars via Core Accretion., Comment: Accepted to ApJ

We perform a sequence of 3D magnetohydrodynamic (MHD) simulations of the outflow-core interaction for a massive protostar forming via collapse of an initial cloud core of $60~{M_\odot}$. This allows us to characterize the properties of disk wind driven outflows from massive protostars, which can allow testing of different massive star formation theories. It also enables us to assess quantitatively the impact of outflow feedback on protostellar core morphology and overall star formation efficiency. We find that the opening angle of the flow increases with increasing protostellar mass, in agreement with a simple semi-analytic model. Once the protostar reaches $\sim24~{M_\odot}$ the outflow's opening angle is so wide that it has blown away most of the envelope, thereby nearly ending its own accretion. We thus find an overall star formation efficiency of $\sim50\%$, similar to that expected from low-mass protostellar cores. Our simulation results therefore indicate that the MHD disk wind outflow is the dominant feedback mechanism for helping to shape the stellar initial mass function from a given prestellar core mass function., Comment: Accepted for publication in ApJ

We conduct a theoretical study of the formation of massive stars over a wide range of metallicities from 1e-5 to 1Zsun and evaluate the star formation efficiencies (SFEs) from prestellar cloud cores taking into account multiple feedback processes. Unlike for simple spherical accretion, in the case of disk accretion feedback processes do not set upper limits on stellar masses. At solar metallicity, launching of magneto-centrifugally-driven outflows is the dominant feedback process to set SFEs, while radiation pressure, which has been regarded to be pivotal, has only minor contribution even in the formation of over-100Msun stars. Photoevaporation becomes significant in over-20Msun star formation at low metallicities of <1e-2Zsun, where dust absorption of ionizing photons is inefficient. We conclude that if initial prestellar core properties are similar, then massive stars are rarer in extremely metal-poor environments of 1e-5 - 1e-3Zsun. Our results give new insight into the high-mass end of the initial mass function and its potential variation with galactic and cosmological environments., Comment: 13 pages, 9 figures, accepted for publication in The Astrophysical Journal

In this series of papers, we model the formation and evolution of the photoionized region and its observational signatures during massive star formation. Here we focus on the early break out of the photoionized region into the outflow cavity. Using results of 3-D magnetohydrodynamic-outflow simulations and protostellar evolution calculations, we perform post-processing radiative-transfer. The photoionized region first appears at a protostellar mass of 10Msun in our fiducial model, and is confined to within 10-100AU by the dense inner outflow, similar to some observed very small hypercompact HII regions. Since the ionizing luminosity of the massive protostar increases dramatically as Kelvin-Helmholz (KH) contraction proceeds, the photoionized region breaks out to the entire outflow region in <10,000yr. Accordingly, the radio free-free emission brightens significantly in this stage. In our fiducial model, the radio luminosity at 10 GHz changes from 0.1 mJy kpc2 at m=11Msun to 100 mJy kpc2 at 16Msun, while the infrared luminosity increases by less than a factor of two. The radio spectral index also changes in the break-out phase from the optically thick value of 2 to the partially optically thin value of 0.6. Additionally, we demonstrate that short-timescale variation in free-free flux would be induced by an accretion burst. The outflow density is enhanced in the accretion burst phase, which leads to a smaller ionized region and weaker free-free emission. The radio luminosity may decrease by one order of magnitude during such bursts, while the infrared luminosity is much less affected, since internal protostellar luminosity dominates over accretion luminosity after KH contraction starts. Such variability may be observable on timescales as short 10-100 yr, if accretion bursts are driven by disk instabilities., Comment: 9 pages, 5 figures, accepted for publication in ApJ

We present calculations of molecular, atomic and ionic line emission from simulations of giant molecular cloud (GMC) collisions. We post-process snapshots of the magneto-hydrodynamical simulations presented in an earlier paper in this series by Wu et al. (2017) of colliding and non-colliding GMCs. Using photodissociation region (PDR) chemistry and radiative transfer we calculate the level populations and emission properties of $^{12}$CO $J=1-0$, [CI] $^3{\rm P}_1\rightarrow{^3{\rm P}}_0$ at $609\,\mu$m, [CII] $158\,\mu$m and [OI] $^3{\rm P}_1\rightarrow{^3{\rm P}}_0$ transition at $63\,\mu$m. From integrated intensity emission maps and position-velocity diagrams, we find that fine-structure lines, particularly the [CII] $158\,\mu$m, can be used as a diagnostic tracer for cloud-cloud collision activity. These results hold even in more evolved systems in which the collision signature in molecular lines has been diminished., Comment: 10 pages, 7 figures, accepted for publication in ApJ, comments welcome

We study feedback during massive star formation using semi-analytic methods, considering the effects of disk winds, radiation pressure, photoevaporation and stellar winds, while following protostellar evolution in collapsing massive gas cores. We find that disk winds are the dominant feedback mechanism setting star formation efficiencies (SFEs) from initial cores of ~0.3-0.5. However, radiation pressure is also significant to widen the outflow cavity causing reductions of SFE compared to the disk-wind only case, especially for >100Msun star formation at clump mass surface densities Sigma<0.3g/cm2. Photoevaporation is of relatively minor importance due to dust attenuation of ionizing photons. Stellar winds have even smaller effects during the accretion stage. For core masses Mc~10-1000Msun and Sigma~0.1-3g/cm2, we find the overall SFE to be 0.31(Rc/0.1pc)^{-0.39}, potentially a useful sub-grid star-formation model in simulations that can resolve pre-stellar core radii, Rc=0.057(Mc/60Msun)^{1/2}(Sigma/g/cm2)^{-1/2}pc. The decline of SFE with Mc is gradual with no evidence for a maximum stellar-mass set by feedback processes up to stellar masses of ~300Msun. We thus conclude that the observed truncation of the high-mass end of the IMF is shaped mostly by the pre-stellar core mass function or internal stellar processes. To form massive stars with the observed maximum masses of ~150-300Msun, initial core masses need to be >500-1000Msun. We also apply our feedback model to zero-metallicity primordial star formation, showing that, in the absence of dust, photoevaporation staunches accretion at ~50Msun. Our model implies radiative feedback is most significant at metallicities ~10^{-2}Zsun, since both radiation pressure and photoevaporation are effective in this regime., Comment: 18 pages, 7 figures, accepted to ApJ

We present an overview and first results of the Stratospheric Observatory For Infrared Astronomy Massive (SOMA) Star Formation Survey, which is using the FORCAST instrument to image massive protostars from $\sim10$--$40\:\rm{\mu}\rm{m}$. These wavelengths trace thermal emission from warm dust, which in Core Accretion models mainly emerges from the inner regions of protostellar outflow cavities. Dust in dense core envelopes also imprints characteristic extinction patterns at these wavelengths, causing intensity peaks to shift along the outflow axis and profiles to become more symmetric at longer wavelengths. We present observational results for the first eight protostars in the survey, i.e., multiwavelength images, including some ancillary ground-based MIR observations and archival {\it{Spitzer}} and {\it{Herschel}} data. These images generally show extended MIR/FIR emission along directions consistent with those of known outflows and with shorter wavelength peak flux positions displaced from the protostar along the blueshifted, near-facing sides, thus confirming qualitative predictions of Core Accretion models. We then compile spectral energy distributions and use these to derive protostellar properties by fitting theoretical radiative transfer models. Zhang and Tan models, based on the Turbulent Core Model of McKee and Tan, imply the sources have protostellar masses $m_*\sim10$--50$\:M_\odot$ accreting at $\sim10^{-4}$--$10^{-3}\:M_\odot\:{\rm{yr}}^{-1}$ inside cores of initial masses $M_c\sim30$--500$\:M_\odot$ embedded in clumps with mass surface densities $\Sigma_{\rm{cl}}\sim0.1$--3$\:{\rm{g\:cm}^{-2}}$. Fitting Robitaille et al. models typically leads to slightly higher protostellar masses, but with disk accretion rates $\sim100\times$ smaller. We discuss reasons for these differences and overall implications of these first survey results for massive star formation theories., Comment: Accepted to ApJ, 32 pages

We present an evolutionary sequence of models of the photoionized disk-wind outflow around forming massive stars based on the Core Accretion model. The outflow is expected to be the first structure to be ionized by the protostar and can confine the expansion of the HII region, especially in lateral directions in the plane of the accretion disk. The ionizing luminosity increases as Kelvin-Helmholz contraction proceeds, and the HII region is formed when the stellar mass reaches ~10-20Msun depending on the initial cloud core properties. Although some part of outer disk surface remains neutral due to shielding by the inner disk and the disk wind, almost the whole of the outflow is ionized in 1e3-1e4 yr after initial HII region formation. Having calculated the extent and temperature structure of the HII region within the immediate protostellar environment, we then make predictions for the strength of its free-free continuum and recombination line emission. The free-free radio emission from the ionized outflow has a flux density of ~(20-200)x(nu/10GHz)^p mJy for a source at a distance of 1 kpc with a spectral index p~0.4-0.7, and the apparent size is typically ~500AU at 10GHz. The H40alpha line profile has a width of about 100km/s. These properties of our model are consistent with observed radio winds and jets around forming massive protostars., Comment: 23 pages, 19 figures, accepted for publication in ApJ

Fragmentation of protostellar disks controls the growth of protostars and plays a key role in determining the final mass of newborn stars. In this paper, we investigate the structure and gravitational stability of the protostellar disks in the full metallicity range between zero and the solar value. Using the mass-accretion rates evaluated from the thermal evolution in the preceding collapse phase of the pre-stellar cores, we calculate disk structures and their evolution in the framework of the standard steady disks. Overall, with higher metallicity, more efficient cooling results in the lower accretion rate and lower temperature inside the disk: at zero metallicity, the accretion rate is ~ 1e-3Msun/yr and the disk temperature is ~ 1000 K, while at solar metallicity, ~ 1e-6Msun/yr and 10 K. Despite the large difference in these values, the zero- and solar-metallicity disks have similar stability properties: the Toomre parameter for the gravitational stability, which can be written using the ratio of temperatures in the disk and in the envelope as Q ~ (T_disk/T_env)^3/2, is > 1, i.e., marginally stable. At intermediate metallicities of 1e-5--1e-3Zsun, however, the disks are found to be strongly unstable with Q ~ 0.1--1 since dust cooling, which is effective only in the disks due to their high density (> 1e10 cm^-3), makes the temperature in the disks lower than that in the envelopes. This indicates that masses of the individual stars formed as a result of the protostellar disk fragmentation can be significantly smaller than their parent core in this metallicity range. The typical stellar mass in this case would be a few Msun, which is consistent with the observationally suggested mass-scale of extremely metal-poor stars., Comment: 15 pages, 8 figures, accepted for publication in MNRAS