No Arabic abstract
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.
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 examine new and pre-existing wide-field, continuum-corrected, narrowband images in H$_2$ 1-0 S(1) and Br$gamma$ of three regions of massive star formation: IC 1396, Cygnus OB2, and Carina. These regions contain a variety of globules, pillars, and sheets, so we can quantify how the spatial profiles of emission lines behave in photodissociation regions (PDRs) that differ in their radiation fields and geometries. We have measured 450 spatial profiles of H$_2$ and Br$gamma$ along interfaces between HII regions and PDRs. Br$gamma$ traces photoevaporative flows from the PDRs, and this emission declines more rapidly with distance as the radius of curvature of the interface decreases, in agreement with models. As noted previously, H$_2$ emission peaks deeper into the cloud relative to Br$gamma$, where the molecular gas absorbs far-UV radiation from nearby O-stars. Although PDRs in IC 1396, Cygnus OB2, and Carina experience orders of magnitude different levels of ionizing flux and have markedly differing geometries, all the PDRs have spatial offsets between Br$gamma$ and H$_2$ on the order of $10^{17}$cm. There is a weak negative correlation between the offset size and the intensity of ionizing radiation and a positive correlation with the radius of curvature of the cloud. We can reproduce both the size of the offsets and the dependencies of the offsets on these other variables with simple photoevaporative flow models. Both Br$gamma$ and H$_2$ 1-0 S(1) will undoubtedly be targeted in future JWST observations of PDRs, so this work can serve as a guide to interpreting these images.
Understanding the chemical evolution of young (high-mass) star-forming regions is a central topic in star formation research. Chemistry is employed as a unique tool 1) to investigate the underlying physical processes and 2) to characterize the evolution of the chemical composition. We observed a sample of 59 high-mass star-forming regions at different evolutionary stages varying from the early starless phase of infrared dark clouds to high-mass protostellar objects to hot molecular cores and, finally, ultra-compact HII regions at 1mm and 3mm with the IRAM 30m telescope. We determined their large-scale chemical abundances and found that the chemical composition evolves along with the evolutionary stages. On average, the molecular abundances increase with time. We modeled the chemical evolution, using a 1D physical model where density and temperature vary from stage to stage coupled with an advanced gas-grain chemical model and derived the best-fit chi^2 values of all relevant parameters. A satisfying overall agreement between observed and modeled column densities for most of the molecules was obtained. With the best-fit model we also derived a chemical age for each stage, which gives the timescales for the transformation between two consecutive stages. The best-fit chemical ages are ~10,000 years for the IRDC stage, ~60,000 years for the HMPO stage, ~40,000 years for the HMC stage, and ~10,000 years for the UCHII stage. The total chemical timescale for the entire evolutionary sequence of the high-mass star formation process is on the order of 10^5 years, which is consistent with theoretical estimates. Furthermore, based on the approach of a multiple-line survey of unresolved data, we were able to constrain an intuitive and reasonable physical and chemical model. The results of this study can be used as chemical templates for the different evolutionary stages in high-mass star formation.
How high-mass stars form remains unclear currently. Calculation suggests that the radiation pressure of a forming star can halt spherical infall, preventing its further growth when it reaches 10 M$_{odot}$. Two major theoretical models on the further growth of stellar mass were proposed. One model suggests the mergence of less massive stellar objects, and the other is still through accretion but with the help of disk. Inflow motions are the key evidence of how forming stars further gain mass to build up massive stars. Recent development in technology has boosted the search of inflow motion. A number of high-mass collapse candidates were obtained with single dish observations, mostly showed blue profile. The infalling signatures seem to be more common in regions with developed radiation pressure than in younger cores, which opposes the theoretical prediction and is also very different from that of low mass star formation. Interferometer studies so far confirm such tendency with more obvious blue profile or inverse P Cygni profile. Results seem to favor the accretion model. However, the evolution tendency of the infall motion in massive star forming cores needs to be further explored. Direct evidence for monolithic or competitive collapse processes is still lack. ALMA will enable us to probe more detail of gravity process.
A large number (67) of the compact/ultra-compact H II regions identified in the Coordinated Radio and Infrared Survey for High-Mass Star Formation catalogue were determined to be powered by a Lyman continuum flux in excess of what was expected given their corresponding luminosity. In this study we attempt to reasonably explain the Lyman excess phenomenon in as many of the 67 H II regions as possible through a variety of observational and astrophysical means including new luminosity estimates, new Herschel photometry, new distance determinations, the use of different models for dust and ionized gas covering factors, and the use of different stellar calibrations. This phenomenon has been observed before; however, the objects shown to exhibit this behaviour in the literature have decidedly different physical properties than the regions in our sample, and thus the origin of the excess is not the same. We find that the excess can be reproduced using OB stellar atmosphere models that have been slightly modified in the extreme ultraviolet. Though the exact mechanism producing the excess is still uncertain, we do find that a scaled up magnetospheric accretion model, often used to explain similar emission from T Tauri stars, is unable to match our observations. Our results suggest that the Lyman excess may be associated with younger H II regions, and that it is more commonly found in early B-type stars. Our refined sample of 24 Lyman excess H II regions provides an ideal sample for comparative studies with regular H II regions, and can act as the basis for the further detailed study of individual regions.