The W3 GMC is a prime target for investigating the formation of high-mass stars and clusters. This second study of W3 within the HOBYS Key Program provides a comparative analysis of subfields within W3 to further constrain the processes leading to th
e observed structures and stellar population. Probability density functions (PDFs) and cumulative mass distributions (CMDs) were created from dust column density maps, quantified as extinction Av. The shape of the PDF, typically represented with a lognormal function at low Av breaking to a power-law tail at high Av, is influenced by various processes including turbulence and self-gravity. The breaks can also be identified, often more readily, in the CMDs. The PDF break from lognormal (Av(SF)= 6-10 mag) appears to shift to higher Av by stellar feedback, so that high-mass star-forming regions tend to have higher PDF breaks. A second break at Av > 50 mag traces structures formed or influenced by a dynamic process. Because such a process has been suggested to drive high-mass star formation in W3, this second break might then identify regions with potential for hosting high-mass stars/clusters. Stellar feedback appears to be a major mechanism driving the local evolution and state of regions within W3. A high initial star formation efficiency in a dense medium could result in a self-enhancing process, leading to more compression and favourable star-formation conditions (e.g., colliding flows), a richer stellar content, and massive stars. This scenario would be compatible with the convergent constructive feedback model introduced in our previous Herschel study.
This paper introduces the data cubes from GHIGLS, deep Green Bank Telescope surveys of the 21-cm line emission of HI in 37 targeted fields at intermediate Galactic latitude. The GHIGLS fields together cover over 1000 square degrees at 9.55 spatial re
solution. The HI spectra have an effective velocity resolution about 1.0 km/s and cover at least -450 < v < +250 km/s. GHIGLS highlights that even at intermediate Galactic latitude the interstellar medium is very complex. Spatial structure of the HI is quantified through power spectra of maps of the column density, NHI. For our featured representative field, centered on the North Ecliptic Pole, the scaling exponents in power-law representations of the power spectra of NHI maps for low, intermediate, and high velocity gas components (LVC, IVC, and HVC) are -2.86 +/- 0.04, -2.69 +/- 0.04, and -2.59 +/- 0.07, respectively. After Gaussian decomposition of the line profiles, NHI maps were also made corresponding to the narrow-line and broad-line components in the LVC range; for the narrow-line map the exponent is -1.9 +/- 0.1, reflecting more small scale structure in the cold neutral medium (CNM). There is evidence that filamentary structure in the HI CNM is oriented parallel to the Galactic magnetic field. The power spectrum analysis also offers insight into the various contributions to uncertainty in the data. The effect of 21-cm line opacity on the GHIGLS NHI maps is estimated.
We use three dimensional hydrodynamical simulations to show that a highly misaligned accretion disk around one component of a binary system can exhibit global Kozai-Lidov cycles, where the inclination and eccentricity of the disk are interchanged per
iodically. This has important implications for accreting systems on all scales, for example, the formation of planets and satellites in circumstellar and circumplanetary disks, outbursts in X-ray binary systems and accretion on to supermassive black holes.
Be/X-ray binary systems exhibit both periodic (Type I) X-ray outbursts and giant (Type II) outbursts, whose origin has remained elusive. We suggest that Type II X-ray outbursts occur when a highly misaligned decretion disk around the Be star becomes
eccentric, allowing the compact object companion to capture a large amount of material at periastron. Using 3D smoothed particle hydrodynamics simulations we model the long term evolution of a representative Be/X-ray binary system. We find that periodic (Type I) X-ray outbursts occur when the neutron star is close to periastron for all disk inclinations. Type II outbursts occur for large misalignment angles and are associated with eccentricity growth that occurs on a timescale of about 10 orbital periods. Mass capture from the eccentric decretion disk results in an accretion disk around the neutron star whose estimated viscous time is long enough to explain the extended duration of Type II outbursts. Previous studies suggested that the outbursts are caused by a warped disk but our results suggest that this is not sufficient, the disk must be both highly misaligned and eccentric to initiate a Type II accretion event.
Protoplanetary disks are likely to be threaded by a weak net flux of vertical magnetic field that is a remnant of the much larger fluxes present in molecular cloud cores. If this flux is approximately conserved its dynamical importance will increase
as mass is accreted, initially by stimulating magnetorotational disk turbulence and subsequently by enabling wind angular momentum loss. We use fits to numerical simulations of ambipolar dominated disk turbulence to construct simplified one dimensional evolution models for weakly magnetized protoplanetary disks. We show that the late onset of significant angular momentum loss in a wind can give rise to two timescale disk evolution in which a long phase of viscous evolution precedes rapid dispersal as the wind becomes dominant. The wide dispersion in disk lifetimes could therefore be due to varying initial levels of net flux. Magnetohydrodynamic (MHD) wind triggered dispersal differs from photoevaporative dispersal in predicting mass loss from small (less that 1 AU) scales, where thermal winds are suppressed. Our specific models are based on a limited set of simulations that remain uncertain, but qualitatively similar evolution appears likely if mass is lost from disks more quickly than flux, and if MHD winds become important as the plasma beta decreases.
A t by n random matrix A is formed by sampling n independent random column vectors, each containing t components. The random Gram matrix of size n, G_n, contains the dot products between all pairs of column vectors in the randomly generated matrix A;
that is, G_n = transpose(A) A. The matrix G_n has characteristic roots coinciding with the singular values of A. Furthermore, the sequences det(G_i) and per(G_i) (for i = 0, 1, ..., n) are factors that comprise the expected coefficients of the characteristic and permanental polynomials of G_n. We prove theorems that relate the generating functions and recursions for the traces of matrix powers, expected characteristic coefficients, expected determinants E(det(G_n)), and expected permanents E(per(G_n)) in terms of each other. Using the derived recursions, we exhibit the efficient computation of the expected determinant and expected permanent of a random Gram matrix G_n, formed according to any underlying distribution. These theoretical results may be used both to speed up numerical algorithms and to investigate the numerical properties of the expected characteristic and permanental coefficients of any matrix comprised of independently sampled columns.
The W3 GMC is a prime target for the study of the early stages of high-mass star formation. We have used Herschel data from the HOBYS key program to produce and analyze column density and temperature maps. Two preliminary catalogs were produced by ex
tracting sources from the column density map and from Herschel maps convolved to the 500 micron resolution. Herschel reveals that among the compact sources (FWHM<0.45 pc), W3 East, W3 West, and W3 (OH) are the most massive and luminous and have the highest column density. Considering the unique properties of W3 East and W3 West, the only clumps with on-going high-mass star formation, we suggest a convergent constructive feedback scenario to account for the formation of a cluster with decreasing age and increasing system/source mass toward the innermost regions. This process, which relies on feedback by high-mass stars to ensure the availability of material during cluster formation, could also lead to the creation of an environment suitable for the formation of Trapezium-like systems. In common with other scenarios proposed in other HOBYS studies, our results indicate that an active/dynamic process aiding in the accumulation, compression, and confinement of material is a critical feature of the high-mass star/cluster formation, distinguishing it from classical low-mass star formation. The environmental conditions and availability of triggers determine the form in which this process occurs, implying that high-mass star/cluster formation could arise from a range of scenarios: from large scale convergence of turbulent flows, to convergent constructive feedback or mergers of filaments.
The submillimeter opacity of dust in the diffuse Galactic interstellar medium (ISM) has been quantified using a pixel-by-pixel correlation of images of continuum emission with a proxy for column density. We used three BLAST bands at 250, 350, and 500
mu m and one IRAS at 100 mu m. The proxy is the near-infrared color excess, E(J-Ks), obtained from 2MASS. Based on observations of stars, we show how well this color excess is correlated with the total hydrogen column density for regions of moderate extinction. The ratio of emission to column density, the emissivity, is then known from the correlations, as a function of frequency. The spectral distribution of this emissivity can be fit by a modified blackbody, whence the characteristic dust temperature T and the desired opacity sigma_e(1200) at 1200 GHz can be obtained. We have analyzed 14 regions near the Galactic plane toward the Vela molecular cloud, mostly selected to avoid regions of high column density (N_H > 10^{22} cm^-2) and small enough to ensure a uniform T. We find sigma_e(1200) is typically 2 to 4 x 10^{-25} cm^2/H and thus about 2 to 4 times larger than the average value in the local high Galactic latitude diffuse atomic ISM. This is strong evidence for grain evolution. There is a range in total power per H nucleon absorbed (re-radiated) by the dust, reflecting changes in the interstellar radiation field and/or the dust absorption opacity. These changes affect the equilibrium T, which is typically 15 K, colder than at high latitudes. Our analysis extends, to higher opacity and lower T, the trend of increasing opacity with decreasing T that was found at high latitudes. The recognition of changes in the emission opacity raises a cautionary flag because all column densities deduced from dust emission maps, and the masses of compact structures within them, depend inversely on the value adopted.
In this work we have carried out an in-depth analysis of the young stellar content in the W3 GMC. The YSO population was identified and classified in the IRAC/MIPS color-magnitude space according to the `Class scheme and compared to other classificat
ions based on intrinsic properties. Class 0/I and II candidates were also compared to low/intermediate-mass pre-main-sequence stars selected through their colors and magnitudes in 2MASS. We find that a reliable color/magnitude selection of low-mass PMS stars in the infrared requires prior knowledge of the protostar population, while intermediate mass objects can be more reliably identified. By means of the MST algorithm and our YSO spatial distribution and age maps we investigated the YSO groups and the star formation history in W3. We find signatures of clustered and distributed star formation in both triggered and quiescent environments. The central/western parts of the GMC are dominated by large scale turbulence likely powered by isolated bursts of star formation that triggered secondary star formation events. Star formation in the eastern high density layer also shows signs of extended periods of star formation. While our findings support triggering as a key factor for inducing and enhancing some of the major star forming activity in the HDL (e.g., W3 Main/W3(OH)), we argue that some degree of quiescent or spontaneous star formation is required to explain the observed YSO population. Our results also support previous studies claiming an spontaneous origin for the isolated massive star(s) powering KR 140.
We present a power spectrum analysis of the Herschel-SPIRE observations of the Polaris flare, a high Galactic latitude cirrus cloud midway between the diffuse and molecular phases. The SPIRE images of the Polaris flare reveal for the first time the s
tructure of the diffuse interstellar medium down to 0.01 parsec over a 10 square degrees region. These exceptional observations highlight the highly filamentary and clumpy structure of the interstellar medium even in diffuse regions of the map. The power spectrum analysis shows that the structure of the interstellar medium is well described by a single power law with an exponent of -2.7 +- 0.1 at all scales from 30 to 8 degrees. That the power spectrum slope of the dust emission is constant down to the SPIRE angular resolution is an indication that the inertial range of turbulence extends down to the 0.01 pc scale. The power spectrum analysis also allows the identification of a Poissonian component at sub-arcminute scales in agreement with predictions of the cosmic infrared background level at SPIRE wavelengths. Finally, the comparison of the SPIRE and IRAS 100 micron data of the Polaris flare clearly assesses the capability of SPIRE in maping diffuse emission over large areas.
