No Arabic abstract
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.
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.
We present molecular line and 1.4 mm continuum observations towards five massive star forming regions at arcsecond resolution using the Submillimeter Array (SMA). We find that the warm molecular gas surrounding each HII region (as traced by SO_2 and OCS) appears to be undergoing bulk rotation. From the molecular line emission and thermal component of the continuum emission, we independently derived gas masses for each region which are consistent with each other. From the free-free component of the continuum emission we estimate the minimum stellar mass required to power the HII region and find that this mass, when added to the derived gas mass, is a significant fraction of the dynamical mass for that region.
We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle-2 observations of the HBC 494 molecular outflow and envelope. HBC 494 is an FU Ori-like object embedded in the Orion A cloud and is associated with the reflection nebulae Re50 and Re50N. We use $^{12}$CO, $^{13}$CO and C$^{18}$O spectral line data to independently describe the outflow and envelope structures associated with HBC 494. The moment-1 map of the $^{12}$CO emission shows the widest outflow cavities in a Class I object known to date (opening angle $sim$ 150$^{^{circ}}$). The morphology of the wide outflow is likely to be due to the interaction between winds originating in the inner disc and the surrounding envelope. The low-velocity blue- and red-shifted $^{13}$CO and C$^{18}$O emission trace the rotation and infall motion of the circumstellar envelope. Using molecular line data and adopting standard methods for correcting optical depth effects, we estimate its kinematic properties, including an outflow mass on the order of 10$^{-1}$ M$_{odot}$. Considering the large estimated outflow mass for HBC 494, our results support recent theoretical work suggesting that wind-driven processes might dominate the evolution of protoplanetary discs via energetic outflows.
Given the rarity of young O star candidates, compact HII regions embedded in dense molecular cores continue to serve as potential sites to peer into the details of high-mass star formation. To uncover the ionizing sources of the most luminous and compact HII regions embedded in the RCW106 and RCW122 giant molecular clouds, known to be relatively nearby (2-4 kpc) and isolated, thus providing an opportunity to examine spatial scales of a few hundred to a thousand AU in size. High spatial resolution (0.3), mid-infrared spectra (R=350), including the fine structure lines [ArIII] and [NeII], were obtained for four luminous compact HII regions, embedded inside the dense cores within the RCW106 and RCW122 molecular cloud complexes. At this resolution, these targets reveal point-like sources surrounded by nebulosity of different morphologies, uncovering details at spatial dimensions of <1000AU. The point-like sources display [ArIII] and [NeII] lines - the ratios of which are used to estimate the temperature of the embedded sources. The derived temperatures are indicative of mid-late O type objects for all the sources with [ArIII] emission. Previously known characteristics of these targets from the literature, including evidence of disk or accretion suggest that the identified sources may grow more to become early-type O stars by the end of the star formation process.
Accurately determining the age of H,{ii} regions and the stars they host is as important as it is challenging. Historically the most popular method has been isochrone fitting to Hertzsprung-Russell Diagrams or Colour-Magnitude Diagrams. Here we introduce a different method for age determination using BPASS and hoki. We infer the most likely ages of the regions D118 and D119 NGC 300 to be log(age/years)=6.86$^{0.05}_{0.06}$ and we also deduce stellar mass and number counts by comparison with the BPASS models. We compare how our binary and single star models perform and find that the latter are unable to predict 20 per cent ($pm$ 10 per cent) of our sample. We also discuss how results obtained from isochrone fitting would differ. We conclude that ages could be underestimated by ~0.2 dex and that the limitations of the isochrone method is not solely due to the lack of binary stars. We propose that the method presented here is more reliable and more widely applicable since it can be used on smaller samples. Alongside this study, we release new hoki features to allow easy implementation of this method.