We offer to the star formation community the third installment of the Massive Star-forming Regions (MSFRs) Omnibus X-ray Catalog (MOXC3), a compilation of X-ray point sources detected in 50 archival Chandra/ACIS observations of 14 Galactic MSFRs and surrounding fields. The MOXC3 MSFRs are NGC 2264, NGC 6193, RCW 108-IR, Aur OB1, DR15, NGC 6231, Berkeley 87, NGC 6357, AFGL 4029, h Per (NGC 869), NGC 281, Onsala 2S, G305, and RCW 49 (Wd 2); they have distances of 0.7 kpc to 4.2 kpc. Most exhibit clumped or clustered young stellar populations; several contain at least two distinct massive young stellar clusters. The total MOXC3 catalog includes 27,923 X-ray point sources. We take great care to identify even the faintest X-ray point sources across these fields. This allows us to remove this point source light, revealing diffuse X-ray structures that pervade and surround MSFRs, often generated by hot plasmas from massive star feedback. As we found in MOXC1 and MOXC2, diffuse X-ray emission is traceable in all MOXC3 MSFRs; here we perform spectral fitting to investigate the origins of selected diffuse regions. Once again, MOXC3 shows the value of high spatial resolution X-ray studies of MSFRs enabled by Chandra.

We present a novel infrared spectral energy distribution (SED) modeling methodology that uses likelihood-based weighting of the model fitting results to construct probabilistic H-R diagrams (pHRD) for X-ray identified, intermediate-mass (2-8 $M_{\odot}$), pre-main sequence young stellar populations. This methodology is designed specifically for application to young stellar populations suffering strong, differential extinction ($\Delta A_V > 10$ mag), typical of Galactic massive star-forming regions. We pilot this technique in the Carina Nebula Complex (CNC) by modeling the 1-8 $\mu$m SEDs of 2269 likely stellar members that exhibit no excess emission from circumstellar dust disks at 4.5 $\mu$m or shorter wavelengths. A subset of ${\sim}100$ intermediate-mass stars in the lightly-obscured Trumpler 14 and 16 clusters have available spectroscopic $T_{\rm eff}$, measured from the Gaia-ESO survey. We correctly identify the stellar temperature in 70% of cases, and the aggregate pHRD for all sources returns the same peak in the stellar age distribution as obtained using the spectroscopic $T_{\rm eff}$. The SED model parameter distributions of stellar mass and evolutionary age reveal significant variation in the duration of star formation among four large-scale stellar overdensities within the CNC and a large distributed stellar population. Star formation began ${\sim}10$ Myr ago and continues to the present day, with the star formation rate peaking ${<}3$ Myr ago when the massive Trumpler 14 and 16 clusters formed. We make public the set of 100,000 SED models generated from standard pre-main sequence evolutionary tracks and our custom software package for generating pHRDs and mass-age distributions from the SED fitting results.

Second only to initial mass, the rate of wind-driven mass loss determines the final mass of a massive star and the nature of its remnant. Motivated by the need to reconcile observational values and theory, we use a recently vetted technique to analyze the mass-loss rates in a sample of OB stars that generate bowshock nebulae. We measure peculiar velocities from new Gaia parallax and proper motion data and their spectral types from new optical and infrared spectroscopy. For our sample of 67 central stars in morphologically selected bowshocks nebulae, 67 are OB stars. The median peculiar velocity is 11 km/s, significantly smaller than classical `runaway star' velocities. Mass-loss rates for these O and early B stars agree with recently lowered theoretical predictions, ranging from ~10^-7 Msun/yr for mid-O dwarfs to 10^-9 Msun/yr for late-O dwarfs---a factor of about 2.7 lower than the often-used Vink et al. (2001) formulation. Our results provide the first observational mass-loss rates for B0--B3 dwarfs and giants---10^-9 to 10^-8 Msun/yr. We find evidence for an increase in the mass-loss rates below a critical effective temperature, consistent with predictions of the bi-stability phenomenon in the range Teff=19,000--27,000 K. The sample exhibits a correlation between modified wind momentum and luminosity, consistent in slope but lower by 0.43 dex in magnitude compared to canonical wind-luminosity relations. We identify a small subset of objects deviating most significantly from theoretical expectations as probable radiation-driven bow wave nebulae by virtue of their low stellar-to-nebular luminosity ratios. For these, the inferred mass-loss rates must be regarded as upper limits.

Citizen science has helped astronomers comb through large data sets to identify patterns and objects that are not easily found through automated processes. The Milky Way Project (MWP), a citizen science initiative on the Zooniverse platform, presents internet users with infrared (IR) images from Spitzer Space Telescope Galactic plane surveys. MWP volunteers make classification drawings on the images to identify targeted classes of astronomical objects. We present the MWP second data release (DR2) and an updated data reduction pipeline written in Python. We aggregate ${\sim}3$ million classifications made by MWP volunteers during the years 2012-2017 to produce the DR2 catalogue, which contains 2600 IR bubbles and 599 candidate bow-shock driving stars. The reliability of bubble identifications, as assessed by comparison to visual identifications by trained experts and scoring by a machine-learning algorithm, is found to be a significant improvement over DR1. We assess the reliability of IR bow shocks via comparison to expert identifications and the colours of candidate bow-shock driving stars in the 2MASS point-source catalogue. We hence identify highly-reliable subsets of 1394 DR2 bubbles and 453 bow-shock driving stars. Uncertainties on object coordinates and bubble size/shape parameters are included in the DR2 catalog. Compared with DR1, the DR2 bubbles catalogue provides more accurate shapes and sizes. The DR2 catalogue identifies 311 new bow shock driving star candidates, including three associated with the giant HII regions NGC 3603 and RCW 49.

It is usually assumed that the stellar initial mass function (IMF) takes a universal form and that there exists a direct mapping between this and the distribution of natal core masses (the core mass function, CMF). The IMF and CMF have been best characterized in the Solar neighborhood. Beyond 500~pc from the Sun, in diverse environments where metallicity varies and massive star feedback may dominate, the IMF has been measured only incompletely and imprecisely, while the CMF has hardly been measured at all. In order to establish if the IMF and CMF are indeed universal and related to each other, it is necessary to: 1) perform multi-wavelength large-scale imaging and spectroscopic surveys of different environments across the Galaxy; 2) require an angular resolution of < 0.1'' in the optical/near-IR for stars and < 5'' in the far-IR for cores; 3) achieve far-IR sensitivities to probe 0.1~Msol cores at 2--3 kpc.

We present a stacking analysis of 2.61 Msec of archival Chandra observations of stellar wind bow shocks. We place an upper limit on the X-ray luminosity of IR-detected bow shocks of $<2\times10^{29}$ erg s$^{-1}$, a more stringent constraint than has been found in previous archival studies and dedicated observing campaigns of nearby bow shocks. We compare the X-ray luminosities and $L_X/L_{\rm bol}$ ratios of bow shock driving stars to those of other OB stars within the Chandra field of view. Driving stars are, on average, of later spectral type than the "field of view" OB stars, and we do not observe any unambiguously high $L_X/L_{\rm bol}$ ratios indicative of magnetic stars in our sample. We additionally asses the feasibility of detecting X-rays from stellar wind bow shocks with the proposed Lynx X-ray Observatory. If the X-ray flux originating from the bow shocks is just below our Chandra detection limit, the nearest bow shock in our sample (at $\sim$0.4 kpc with an absorbing column of $\sim10^{21}$ cm$^{-2}$) should be observable with Lynx in exposure times on the order of $\sim$100 kiloseconds.

We present a multiwavelength study of 28 Galactic massive star-forming H II regions. For 17 of these regions, we present new distance measurements based on Gaia DR2 parallaxes. By fitting a multicomponent dust, blackbody, and power-law continuum model to the 3.6 $\mu$m through 10 mm spectral energy distributions, we find that ${\sim}34$% of Lyman continuum photons emitted by massive stars are absorbed by dust before contributing to the ionization of H II regions, while ${\sim}68$% of the stellar bolometric luminosity is absorbed and reprocessed by dust in the H II regions and surrounding photodissociation regions. The most luminous, infrared-bright regions that fully sample the upper stellar initial mass function (ionizing photon rates $N_C \ge 10^{50}~{\rm s}^{-1}$ and dust-processed $L_{\rm TIR}\ge 10^{6.8}$ L$_{\odot}$) have on average higher percentages of absorbed Lyman continuum photons ($\sim$51%) and reprocessed starlight ($\sim$82%) compared to less luminous regions. Luminous H II regions show lower average PAH fractions than less luminous regions, implying that the strong radiation fields from early-type massive stars are efficient at destroying PAH molecules. On average, the monochromatic luminosities at 8, 24, and 70 $\mu$m combined carry 94% of the dust-reprocessed $L_{\rm TIR}$. $L_{70}$ captures ${\sim}52$% of $L_{\rm TIR}$, and is therefore the preferred choice to infer the bolometric luminosity of dusty star-forming regions. We calibrate SFRs based on $L_{24}$ and $L_{70}$ against the Lyman continuum photon rates of the massive stars in each region. Standard extragalactic calibrations of monochromatic SFRs based on population synthesis models are generally consistent with our values.

Photometric detections of dust circumstellar disks around pre-main sequence (PMS) stars, coupled with estimates of stellar ages, provide constraints on the time available for planet formation. Most previous studies on disk longevity, starting with Haisch, Lada & Lada (2001), use star samples from PMS clusters but do not consider datasets with homogeneous photometric sensitivities and/or ages placed on a uniform timescale. Here we conduct the largest study to date of the longevity of inner dust disks using X-ray and 1--8 micrometre infrared photometry from the MYStIX and SFiNCs projects for 69 young clusters in 32 nearby star-forming regions with ages t<=5 Myr. Cluster ages are derived by combining the empirical AgeJX method with PMS evolutionary models, which treat dynamo-generated magnetic fields in different ways. Leveraging X-ray data to identify disk-free objects, we impose similar stellar mass sensitivity limits for disk-bearing and disk-free YSOs while extending the analysis to stellar masses as low as M~0.1 Mo. We find that the disk longevity estimates are strongly affected by the choice of PMS evolutionary model. Assuming a disk fraction of 100% at zero age, the inferred disk half-life changes significantly, from t1/2 ~ 1.3--2 Myr to t1/2 ~ 3.5 Myr when switching from non-magnetic to magnetic PMS models. In addition, we find no statistically significant evidence that disk fraction varies with stellar mass within the first few Myr of life for stars with masses <2 Mo, but our samples may not be complete for more massive stars. The effects of initial disk fraction and star-forming environment are also explored.

The Carina Nebula is an active star forming region in the southern sky that is of particular interest due to the presence of a large number of massive stars in a wide array of evolutionary stages. Here we present the results of the spectroscopic analysis of 82 B-type stars and 33 O-type stars that were observed in 2013 and 2014. For 82 B-type stars without line blending, we fit model spectra from the Tlusty BSTAR2006 grid to the observed profiles of H{\gamma} and He {\lambda}{\lambda} 4026, 4388, and 4471 to measure the effective temperatures, surface gravities, and projected rotational velocities. We also measure the masses, ages, radii, bolometric luminosities, and distances of these stars. From the radial velocities measured in our sample, we find 31 single lined spectroscopic binary candidates. We find a high dispersion of radial velocities among our sample stars, and we argue that the Carina Nebula stellar population has not yet relaxed and become virialized.

The rate at which massive stars eject mass in stellar winds significantly influences their evolutionary path. Cosmic rates of nucleosynthesis, explosive stellar phenomena, and compact object genesis depend on this poorly known facet of stellar evolution. We employ an unexploited observational technique for measuring the mass-loss rates of O- and early-B stars. Our approach, which has no adjustable parameters, uses the principle of pressure equilibrium between the stellar wind and the ambient interstellar medium for a high-velocity star generating an infrared bowshock nebula. Results for twenty bowshock-generating stars show good agreement with two sets of theoretical predictions for O5--O9.5 main-sequence stars, yielding $\dot M=$1.3$\times$10$^{-6}$ to 2$\times$10$^{-9}$ solar masses per year. Although $\dot M$ values derived for this sample are smaller than theoretical expectations by a factor of about two, this discrepancy is greatly reduced compared to canonical mass-loss methods. Bowshock-derived mass-loss rates are factors of ten smaller than H$\alpha$-based measurements (uncorrected for clumping) for similar stellar types and are nearly an order of magnitude larger than P$^{4+}$ and some other UV absorption-line-based diagnostics. Ambient interstellar densities of at least several cm$^{-3}$ appear to be required for formation of a prominent infrared bowshock nebula. $\dot M$ measurements for early-B stars are not yet compelling owing to the small number in our sample and the lack of clear theoretical predictions in the regime of lower stellar luminosities. These results may constitute a partial resolution of the extant "weak-wind problem" for late-O stars. The technique shows promise for determining mass-loss rates in the weak-wind regime.

We present the second installment of the Massive Star-forming Regions (MSFRs) Omnibus X-ray Catalog (MOXC2), a compilation of X-ray point sources detected in Chandra/ACIS observations of 16 Galactic MSFRs and surrounding fields. MOXC2 includes 13 ACIS mosaics, three containing a pair of unrelated MSFRs at different distances, with a total catalog of 18,396 point sources. The MSFRs sampled range over distances of 1.3 kpc to 6 kpc and populations varying from single massive protostars to the most massive Young Massive Cluster known in the Galaxy. By carefully detecting and removing X-ray point sources down to the faintest statistically-significant limit, we facilitate the study of the remaining unresolved X-ray emission. Through comparison with mid-infrared images that trace photon-dominated regions and ionization fronts, we see that the unresolved X-ray emission is due primarily to hot plasmas threading these MSFRs, the result of feedback from the winds and supernovae of massive stars. The 16 MSFRs studied in MOXC2 more than double the MOXC1 sample, broadening the parameter space of ACIS MSFR explorations and expanding Chandra's substantial contribution to contemporary star formation science.

Arcuate infrared nebulae are ubiquitous throughout the Galactic Plane and are candidates for partial shells, bubbles, or bowshocks produced by massive runaway stars. We tabulate infrared photometry for 709 such objects using images from the Spitzer Space Telescope (SST), Wide-Field Infrared Explorer (WISE), and Herschel Space Observatory (HSO). Of the 709 objects identified at 24 or 22 microns, 422 are detected at the HSO 70 micron bandpass. Of these, only 39 are detected at HSO 160 microns. The 70 micron peak surface brightnesses are 0.5 to 2.5 Jy/square arcminute. Color temperatures calculated from the 24 micron to 70 micron ratios range from 80 K to 400 K. Color temperatures from 70 micron to 160 micron ratios are systematically lower, 40 K to 200 K. Both of these temperature are, on average, 75% higher than the nominal temperatures derived by assuming that dust is in steady-state radiative equilibrium. This may be evidence of stellar wind bowshocks sweeping up and heating --- possibly fragmenting but not destroying --- interstellar dust. Infrared luminosity correlates with standoff distance, R_0, as predicted by published hydrodynamical models. Infrared spectral energy distributions are consistent with interstellar dust exposed to a either single radiant energy density, U=10^3 to 10^5 (in more than half of the objects) or a range of radiant energy densities U_min=25 to U_max=10^3 to 10^5 times the mean interstellar value for the remainder. Hence, the central OB stars dominate the energetics, making these enticing laboratories for testing dust models in constrained radiation environments. SEDs are consistent with PAH fractions q_PAH <1% in most objects.

Massive, O and early B-type (OB) stars remain incompletely catalogued in the nearby Galaxy due to high extinction, bright visible and infrared nebular emission in H II regions, and high field star contamination. These difficulties are alleviated by restricting the search to stars with X-ray emission. Using the X-ray point sources from the Massive Young star-forming complex Study in Infrared and X-rays (MYStIX) survey of OB-dominated regions, we identify 98 MYStIX candidate OB (MOBc) stars by fitting their 1-8 micron spectral energy distributions (SEDs) with reddened stellar atmosphere models. We identify 27 additional MOBc stars based on JHK photometry of X-ray stars lacking SED fitting. These candidate OB stars indicate that the current census of stars earlier than B1, taken across the 18 MYStIX regions studied, is less than 50% complete. We also fit the SEDs of 239 previously-published OB stars to measure interstellar extinction and bolometric luminosities, revealing six candidate massive binary systems and five candidate O-type (super)giants. As expected, candidate OB stars have systematically higher extinction than previously-published OB stars. Notable results for individual regions include: identification of the OB population of a recently discovered massive cluster in NGC 6357; an older OB association in the M17 complex; and new massive luminous O stars near the Trifid Nebula. In several relatively poorly-studied regions (RCW 38, NGC 6334, NGC 6357, Trifid, and NGC 3576), the OB populations may increase by factors of >2.

The Carina star-forming region is one of the largest in the Galaxy, and its massive star population is still being unveiled. The large number of stars combined with high, and highly variable, interstellar extinction makes it inherently difficult to find OB stars in this type of young region. We present the results of a spectroscopic campaign to study the massive star population of the Carina Nebula, with the primary goal to confirm or reject previously identified Carina OB star candidates. A total of 141 known O- and B-type stars and 94 candidates were observed, of which 73 candidates had a high enough signal-to-noise ratio to classify. We find 23 new OB stars within the Carina Nebula, a 32% confirmation rate. One of the new OB stars has blended spectra and is suspected to be a double-lined spectroscopic binary (SB2). We also reclassify the spectral types of the known OB stars and discover nine new SB2s among this population. Finally, we discuss the spatial distribution of these new OB stars relative to known structures in the Carina NSebula.

The Massive Young Star Forming Complex in Infrared and X-ray (MYStIX) project provides a new census on stellar members of massive star forming regions within 4 kpc. Here the MYStIX Infrared Excess catalog (MIRES) and Chandra-based X-ray photometric catalogs are mined to obtain high-quality samples of Class I protostars using criteria designed to reduce extragalactic and Galactic field star contamination. A total of 1,109 MYStIX Candidate Protostars (MCPs) are found in 14 star forming regions. Most are selected from protoplanetary disk infrared excess emission, but 20% are found from their ultrahard X-ray spectra from heavily absorbed magnetospheric flare emission. Two-thirds of the MCP sample is newly reported here. The resulting samples are strongly spatially associated with molecular cores and filaments on Herschel far-infrared maps. This spatial agreement and other evidence indicate that the MCP sample has high reliability with relatively few 'false positives' from contaminating populations. But the limited sensitivity and sparse overlap among the infrared and X-ray subsamples indicate that the sample is very incomplete with many 'false negatives'. Maps, tables, and source descriptions are provided to guide further study of star formation in these regions. In particular, the nature of ultrahard X-ray protostellar candidates without known infrared counterparts needs to be eludicated.

We identify 709 arc-shaped mid-infrared nebula in 24 micron Spitzer Space Telescope or 22 micron Wide Field Infrared Explorer surveys of the Galactic Plane as probable dusty interstellar bowshocks powered by early-type stars. About 20% are visible at 8 microns or shorter mid-infrared wavelengths as well. The vast majority (660) have no previous identification in the literature. These extended infrared sources are strongly concentrated near Galactic mid-Plane with an angular scale height of ~0.6 degrees. All host a symmetrically placed star implicated as the source of a stellar wind sweeping up interstellar material. These are candidate "runaway" stars potentially having high velocities in the reference frame of the local medium. Among the 286 objects with measured proper motions, we find an unambiguous excess having velocity vectors aligned with the infrared morphology --- kinematic evidence that many of these are "runaway" stars with large peculiar motions responsible for the bowshock signature. We discuss a population of "in-situ" bowshocks (103 objects) that face giant HII regions where the relative motions between the star and ISM may be caused by bulk outflows from an overpressured bubble. We also identify 58 objects that face 8 micron bright-rimmed clouds and apparently constitute a sub-class of in-situ bowshocks where the stellar wind interacts with a photo-evaporative flow from an eroding molecular cloud interface (i.e., "PEF bowshocks"). Orientations of the arcuate nebulae exhibit a correlation over small angular scales, indicating that external influences such as HII regions are responsible for producing some bowshock nebulae. However, the vast majority of this sample appear to be isolated (499 objects) from obvious external influences.

We present a catalog of 840 X-ray sources and first results from a 100 ks Chandra X-ray Observatory imaging study of the filamentary infrared dark cloud G014.225$-$00.506, which forms the central regions of a larger cloud complex known as the M17 southwest extension (M17 SWex). In addition to the rich population of protostars and young stellar objects with dusty circumstellar disks revealed by Spitzer Space Telescope archival data, we discover a population of X-ray-emitting, intermediate-mass pre--main-sequence stars (IMPS) that lack infrared excess emission from circumstellar disks. We model the infrared spectral energy distributions of this source population to measure its mass function and place new constraints on the inner dust disk destruction timescales for 2-8 $M_{\odot}$ stars. We also place a lower limit on the star formation rate (SFR) and find that it is quite high ($\dot{M}\ge 0.007~M_{\odot}$ yr$^{-1}$), equivalent to several Orion Nebula Clusters in G14.225$-$0.506 alone, and likely accelerating. The cloud complex has not produced a population of massive, O-type stars commensurate with its SFR. This absence of very massive (${\ge}20~M_{\odot}$) stars suggests that either (1) M17 SWex is an example of a distributed mode of star formation that will produce a large OB association dominated by intermediate-mass stars but relatively few massive clusters, or (2) the massive cores are still in the process of accreting sufficient mass to form massive clusters hosting O stars.

We present a statistical study of the distribution and physical properties of cold dense material in and around the inner Galactic Plane near infrared bubbles as catalogued by the Milky Way Project citizen scientists. Using data from the ATLASGAL 870 um survey, we show that 48 +/- 2% of all cold clumps in the studied survey region (|l| <= 65 degrees, |b| <= 1 degree) are found in close proximity to a bubble, and 25 +/- 2% appear directly projected towards a bubble rim. A two-point correlation analysis confirms the strong correlation of massive cold clumps with expanding bubbles. It shows an overdensity of clumps along bubble rims that grows with increasing bubble size, which shows how interstellar medium material is reordered on large scales by bubble expansion around regions of massive star formation. The highest column density clumps appear resistent to the expansion, remaining overdense towards the bubbles' interior rather than being swept up by the expanding edge. Spectroscopic observations in ammonia show that cold dust clumps near bubbles appear to be denser, hotter and more turbulent than those in the field, offering circumstantial evidence that bubble-associated clumps are more likely to be forming stars. These observed differences in physical conditions persist for beyond the region of the bubble rims.

Nearby Type III (galaxy-spanning) Kardashev supercivilizations would have high mid-infrared (MIR) luminosities. We have used the Wide-field Infrared Survey Explorer (WISE) to survey ~$1 \times 10^5$ galaxies for extreme MIR emission, $10^3$ times more galaxies than the only previous such search. We have calibrated the WISE All-sky Catalog pipeline products to improve its photometry for extended sources. We present 563 extended sources with $|b| \ge 10$ and red MIR colors, having visually vetted them to remove artifacts. No galaxies in our sample host an alien civilization reprocessing more than 85% of its starlight into the MIR, and only 50 galaxies, including Arp 220, have MIR luminosities consistent with >50% reprocessing. Ninety of these (likely) extragalactic sources have little literature presence; in most cases they are likely barely resolved galaxies or pairs of galaxies undergoing large amounts of star formation. Five are new to science and deserve further study. The Be star 48 Librae sits within a MIR nebula, and we suggest that it may be creating dust. WISE, 2MASS, and Spitzer imagery shows that IRAS 04287+6444 is consistent with a previously unnoticed, heavily extinguished cluster of young stellar objects. We identify five "passive" (i.e. red) spiral galaxies with unusually high MIR and low NUV luminosity. We search a set of optically "dark" HI galaxies for MIR emission, and find none. These 90 poorly understood sources and five anomalous passive spirals deserve follow-up via both SETI and conventional astrophysics.

We report the results of infrared photometry on 39 extended red objects (EROs) in the Carina Nebula, observed with the Spitzer Space Telescope. Most EROs are identified by bright, extended 8.0 um emission, which ranges from 10'' to 40'' in size, but our sample also includes 4 EROs identified by extended 24 um emission. Of particular interest are nine EROs associated with late O or early B-type stars and characterized by arc-shaped morphology, suggesting dusty, stellar-wind bow shocks. These objects are preferentially oriented towards the central regions of the Carina Nebula, suggesting that these bow shocks are generally produced by the interactions of OB winds with the bulk expansion of the H II region rather than high proper motion. We identify preferred regions of mid-infrared color space occupied by our bow shock candidates, which also contain bow shock candidates in M17 and RCW 49 but are well-separated from polycyclic aromatic hydrocarbon emission or circumstellar discs. Color cuts identify an additional 12 marginally-resolved bow shock candidates, 10 of which are also associated with known late O or early B stars. H II region expansion velocities derived from bow shock candidate standoff distances are ~10 km/s, assuming typical H II region gas densities, comparable to expansion velocities derived from bow shocks in M17 and RCW 49. One candidate bow shock provides direct evidence of physical interaction between the massive stellar winds originating in the Trumpler 15 and Trumpler 14 clusters, supporting the conclusion that both clusters are at similar heliocentric distances.

We motivate the \^G infrared search for extraterrestrial civilizations with large energy supplies. We discuss some philosophical difficulties of SETI, and how communication SETI circumvents them. We review "Dysonian SETI", the search for artifacts of alien civilizations, and find that it is highly complementary to traditional communication SETI; the two together might succeed where either one, alone, has not. We discuss the argument of Hart (1975) that spacefaring life in the Milky Way should be either galaxy-spanning or non-existent, and examine a portion of his argument that we dub the "monocultural fallacy". We discuss some rebuttals to Hart that invoke sustainability and predict long Galaxy colonization timescales. We find that the maximum Galaxy colonization timescale is actually much shorter than previous work has found ($< 10^9$ yr), and that many "sustainability" counter-arguments to Hart's thesis suffer from the monocultural fallacy. We extend Hart's argument to alien energy supplies, and argue that detectably large energy supplies can plausibly be expected to exist because life has potential for exponential growth until checked by resource or other limitations, and intelligence implies the ability to overcome such limitations. As such, if Hart's thesis is correct then searches for large alien civilizations in other galaxies may be fruitful; if it is incorrect, then searches for civilizations within the Milky Way are more likely to succeed than Hart argued. We review some past Dysonian SETI efforts, and discuss the promise of new mid-infrared surveys, such as that of WISE.

We describe the framework and strategy of the \^G infrared search for extraterrestrial civilizations with large energy supplies, which will use the wide-field infrared surveys of WISE and Spitzer to search for these civilizations' waste heat. We develop a formalism for translating mid-infrared photometry into quantitative upper limits on extraterrestrial energy supplies. We discuss the likely sources of false positives, how dust can and will contaminate our search, and prospects for distinguishing dust from alien waste heat. We argue that galaxy-spanning civilizations may be easier to distinguish from natural sources than circumstellar civilizations (i.e., Dyson spheres), although Gaia will significantly improve our capability to identify the latter. We present a "zeroth order" null result of our search based on the WISE all-sky catalog: we show, for the first time, that Kardashev Type III civilizations (as Kardashev originally defined them) are very rare in the local universe. More sophisticated searches can extend our methodology to smaller waste heat luminosities, and potentially entirely rule out (or detect) both Kardashev Type III civilizations and new physics that allows for unlimited "free" energy generation.

The clusters of young stars in massive star-forming regions show a wide range of sizes, morphologies, and numbers of stars. Their highly subclustered structures are revealed by the MYStIX project's sample of 31,754 young stars in nearby sites of star formation (regions at distances <3.6 kpc that contain at least one O-type star.) In 17 of the regions surveyed by MYStIX, we identify subclusters of young stars using finite mixture models -- collections of isothermal ellipsoids that model individual subclusters. Maximum likelihood estimation is used to estimate the model parameters, and the Akaike Information Criterion is used to determine the number of subclusters. This procedure often successfully finds famous subclusters, such as the BN/KL complex behind the Orion Nebula Cluster and the KW-object complex in M17. A catalog of 142 subclusters is presented, with 1 to 20 subclusters per region. The subcluster core radius distribution for this sample is peaked at 0.17 pc with a standard deviation of 0.43 dex, and subcluster core radius is negatively correlated with gas/dust absorption of the stars -- a possible age effect. Based on the morphological arrangements of subclusters, we identify four classes of spatial structure: long chains of subclusters, clumpy structures, isolated clusters with a core-halo structure, and isolated clusters well fit by a single isothermal ellipsoid.

We present the Massive Star-forming Regions (MSFRs) Omnibus X-ray Catalog (MOXC), a compendium of X-ray point sources from {\em Chandra}/ACIS observations of a selection of MSFRs across the Galaxy, plus 30 Doradus in the Large Magellanic Cloud. MOXC consists of 20,623 X-ray point sources from 12 MSFRs with distances ranging from 1.7 kpc to 50 kpc. Additionally, we show the morphology of the unresolved X-ray emission that remains after the catalogued X-ray point sources are excised from the ACIS data, in the context of \Spitzer\ and {\em WISE} observations that trace the bubbles, ionization fronts, and photon-dominated regions that characterize MSFRs. In previous work, we have found that this unresolved X-ray emission is dominated by hot plasma from massive star wind shocks. This diffuse X-ray emission is found in every MOXC MSFR, clearly demonstrating that massive star feedback (and the several-million-degree plasmas that it generates) is an integral component of MSFR physics.

A major impediment to understanding star formation in massive star forming regions (MSFRs) is the absence of a reliable stellar chronometer to unravel their complex star formation histories. We present a new estimation of stellar ages using a new method that employs near-infrared (NIR) and X-ray photometry, AgeJX. Stellar masses are derived from X-ray luminosities using the Lx - Mass relation from the Taurus cloud. J-band luminosities are compared to mass-dependent pre-main-sequence evolutionary models to estimate ages. AgeJX is sensitive to a wide range of evolutionary stages, from disk-bearing stars embedded in a cloud to widely dispersed older pre-main sequence stars. The MYStIX (Massive Young Star-Forming Complex Study in Infrared and X-ray) project characterizes 20 OB-dominated MSFRs using X-ray, mid-infrared, and NIR catalogs. The AgeJX method has been applied to 5525 out of 31,784 MYStIX Probable Complex Members. We provide a homogeneous set of median ages for over a hundred subclusters in 15 MSFRs; median subcluster ages range between 0.5 Myr and 5 Myr. The important science result is the discovery of age gradients across MYStIX regions. The wide MSFR age distribution appears as spatially segregated structures with different ages. The AgeJX ages are youngest in obscured locations in molecular clouds, intermediate in revealed stellar clusters, and oldest in distributed populations. The NIR color index J-H, a surrogate measure of extinction, can serve as an approximate age predictor for young embedded clusters.

MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula.

Spitzer IRAC observations and stellar photometric catalogs are presented for the Massive Young Star-Forming Complex Study in the Infrared and X-ray (MYStIX). MYStIX is a multiwavelength census of young stellar members of twenty nearby (d < 4 kpc), Galactic, star-forming regions that contain at least one O star. All regions have data available from the Spitzer Space Telescope, consisting of GLIMPSE or other published catalogs for eleven regions and results of our own photometric analysis of archival data for the remaining nine regions. This paper seeks to construct deep and reliable catalogs of sources from the Spitzer images. Mid-infrared study of these regions faces challenges of crowding and high nebulosity. Our new catalogs typically contain fainter sources than existing Spitzer studies, which improves the match rate to Chandra X-ray sources that are likely to be young stars, but increases the possibility of spurious point-source detections, especially peaks in the nebulosity. IRAC color-color diagrams help distinguish spurious detections of nebular PAH emission from the infrared excess associated with dusty disks around young stars. The distributions of sources on the mid-infrared color-magnitude and color-color diagrams reflect differences between MYStIX regions, including astrophysical effects such as stellar ages and disk evolution.

The MYStIX project (Massive Young Star-Forming Complex Study in Infrared and X-rays) provides a comparative study of 20 Galactic massive star-forming complexes (d = 0.4 to 3.6 kpc). Probable stellar members in each target complex are identified using X-ray and/or infrared data via two pathways: (1) X-ray detections of young/massive stars with coronal activity/strong winds; or (2) infrared excess (IRE) selection of young stellar objects (YSOs) with circumstellar disks and/or protostellar envelopes. We present the methodology for the second pathway, using Spitzer/IRAC, 2MASS, and UKIRT imaging and photometry. Although IRE selection of YSOs is welltrodden territory, MYStIX presents unique challenges. We combine IR spectral energy distribution (SED) fitting with IR color cuts and spatial clustering analysis to identify IRE sources and isolate probable YSO members in each MYStIX target field from the myriad types of contaminating sources that can resemble YSOs: extragalactic sources, evolved stars, nebular knots, and even unassociated foreground/background YSOs. Applying our methodology consistently across 18 of the target complexes, we produce the MYStIX IRE Source (MIRES) Catalog comprising 20,719 sources, including 8686 probable stellar members of the MYStIX target complexes. We also classify the SEDs of 9365 IR counterparts to MYStIX X-ray sources to assist the first pathway, the identification of X-ray detected stellar members. The MIRES catalog provides a foundation for follow-up studies of diverse phenomena related to massive star cluster formation, including protostellar outflows, circumstellar disks, and sequential star formation triggered by massive star feedback processes.

The Massive Young star-forming Complex Study in Infrared and X-rays (MYStIX) project requires samples of young stars that are likely members of 20 nearby Galactic massive star-forming regions. Membership is inferred from statistical classification of X-ray sources, from detection of a robust infrared excess that is best explained by circumstellar dust in a disk or infalling envelope, and from published spectral types that are unlikely to be found among field stars. We present the MYStIX membership lists here, and describe in detail the statistical classification of X-ray sources via a \Naive Bayes Classi

Observations of the spatial distributions of young stars in star-forming regions can be linked to the theory of clustered star formation using spatial statistical methods. The MYStIX project provides rich samples of young stars from the nearest high-mass star-forming regions. Maps of stellar surface density reveal diverse structure and subclustering. Young stellar clusters and subclusters are fit with isothermal spheres and ellipsoids using the Bayesian Information Criterion to estimate the number of subclusters. Clustering is also investigated using Cartwright and Whitworth's Q statistic and the inhomogeneous two-point correlation function. Mass segregation is detected in several cases, in both centrally concentrated and fractally structured star clusters, but a few clusters are not mass segregated.

The Milky Way Project citizen science initiative recently increased the number of known infrared bubbles in the inner Galactic plane by an order of magnitude compared to previous studies. We present a detailed statistical analysis of this dataset with the Red MSX Source catalog of massive young stellar sources to investigate the association of these bubbles with massive star formation. We particularly address the question of massive triggered star formation near infrared bubbles. We find a strong positional correlation of massive young stellar objects (MYSOs) and H II regions with Milky Way Project bubbles at separations of < 2 bubble radii. As bubble sizes increase, a statistically significant overdensity of massive young sources emerges in the region of the bubble rims, possibly indicating the occurrence of triggered star formation. Based on numbers of bubble-associated RMS sources we find that 67+/-3% of MYSOs and (ultra)compact H II regions appear associated with a bubble. We estimate that approximately 22+/-2% of massive young stars may have formed as a result of feedback from expanding H II regions. Using MYSO-bubble correlations, we serendipitously recovered the location of the recently discovered massive cluster Mercer 81, suggesting the potential of such analyses for discovery of heavily extincted distant clusters.

Massive stars produce copious quantities of ultraviolet radiation beyond the Lyman limit, photoionizing the interstellar medium (ISM) and producing H II regions. As strong sources of recombination- and forbidden-line emission, infrared continuum, and thermal (free-free) radio continuum, H II regions serve as readily-observable beacons of massive star formation in the Milky Way and external galaxies. Along with supernovae, H II regions are dominant sources of feedback in star-forming galaxies, injecting radiative and mechanical luminosity into the ISM. H II regions may prove more important than supernovae as triggers of star formation through localized compression of cold cloud cores. In this review, I give a broad overview of the structure and time-evolution of H II regions, emphasizing complications to the theoretical picture revealed by multiwavelength observations. I discuss a recent controversy surrounding the dominant feedback mechanism in 30 Doradus, the most luminous H II region in the Local Group. I summarize the first results from the Milky Way Project (MWP), which has produced a new catalog of several thousand candidate Galactic H II regions by enlisting >35,000 "citizen scientists" to search Spitzer Space Telescope survey images for bubble-shaped structures. The MWP and similar large catalogs enable empirical studies of Galactic H II region evolution across the full range of luminosities and statistical studies of triggered star formation.

We present a new catalogue of 5,106 infrared bubbles created through visual classification via the online citizen science website 'The Milky Way Project'. Bubbles in the new catalogue have been independently measured by at least 5 individuals, producing consensus parameters for their position, radius, thickness, eccentricity and position angle. Citizen scientists - volunteers recruited online and taking part in this research - have independently rediscovered the locations of at least 86% of three widely-used catalogues of bubbles and H ii regions whilst finding an order of magnitude more objects. 29% of the Milky Way Project catalogue bubbles lie on the rim of a larger bubble, or have smaller bubbles located within them, opening up the possibility of better statistical studies of triggered star formation. Also outlined is the creation of a 'heat map' of star-formation activity in the Galactic plane. This online resource provides a crowd-sourced map of bubbles and arcs in the Milky Way, and will enable better statistical analysis of Galactic star-formation sites.

The star formation rate (SFR) of the Milky Way remains poorly known, with often-quoted values ranging from 1 to 10 solar masses per year. This situation persists despite the potential for the Milky Way to serve as the ultimate SFR calibrator for external galaxies. We show that various estimates for the Galactic SFR are consistent with one another once they have been normalized to the same initial mass function (IMF) and massive star models, converging to 1.9 +/- 0.4 M_sun/yr. However, standard SFR diagnostics are vulnerable to systematics founded in the use of indirect observational tracers sensitive only to high-mass stars. We find that absolute SFRs measured using resolved low/intermediate-mass stellar populations in Galactic H II regions are systematically higher by factors of ~2-3 as compared with calibrations for SFRs measured from mid-IR and radio emission. We discuss some potential explanations for this discrepancy and conclude that it could be allayed if (1) the power-law slope of the IMF for intermediate-mass (1.5 M_sun < m < 5 M_sun) stars were steeper than the Salpeter slope, or (2) a correction factor was applied to the extragalactic 24 micron SFR calibrations to account for the duration of star formation in individual mid-IR-bright H II regions relative to the lifetimes of O stars. Finally, we present some approaches for testing if a Galactic SFR of ~2 M_sun/yr is consistent with what we would measure if we could view the Milky Way as external observers. Using luminous radio supernova remnants and X-ray point sources, we find that the Milky Way deviates from expectations at the 1-3 sigma level, hinting that perhaps the Galactic SFR is overestimated or extragalactic SFRs need to be revised upwards.

We present a catalog of 1439 young stellar objects (YSOs) spanning the 1.42 deg^2 field surveyed by the Chandra Carina Complex Project (CCCP), which includes the major ionizing clusters and the most active sites of ongoing star formation within the Great Nebula in Carina. Candidate YSOs were identified via infrared (IR) excess emission from dusty circumstellar disks and envelopes, using data from the Spitzer Space Telescope Vela--Carina survey and the Two-Micron All Sky Survey. We model the 1--24 /mu m IR spectral energy distributions of the YSOs to constrain physical properties. Our Pan-Carina YSO Catalog (PCYC) is dominated by intermediate-mass (2 Msun < m < 10 Msun) objects with disks, including Herbig Ae/Be stars and their less evolved progenitors. The PCYC provides a valuable complementary dataset to the CCCP X-ray source catalogs, identifying 1029 YSOs in Carina with no X-ray detection. We also catalog 410 YSOs with X-ray counterparts, including 62 candidate protostars. Candidate protostars with X-ray detections tend to be more evolved than those without. In most cases, X-ray emission apparently originating from intermediate-mass, disk-dominated YSOs is consistent with the presence of low-mass companions, but we also find that X-ray emission correlates with cooler stellar photospheres and higher disk masses. We suggest that intermediate-mass YSOs produce X-rays during their early pre-main sequence evolution, perhaps driven by magnetic dynamo activity during the convective atmosphere phase, but this emission dies off as the stars approach the main sequence. Extrapolating over the stellar initial mass function scaled to the PCYC population, we predict a total population of >2x10^4 YSOs and a present-day star formation rate (SFR) of >0.008 Msun/yr. The global SFR in the Carina Nebula, averaged over the past ~5 Myr, has been approximately constant.

The Chandra Carina Complex contains 200 known O- and B type stars. The Chandra survey detected 68 of the 70 O stars and 61 of 127 known B0-B3 stars. We have assembled a publicly available optical/X-ray database to identify OB stars that depart from the canonical Lx/Lbol relation, or whose average X-ray temperatures exceed 1 keV. Among the single O stars with high kT we identify two candidate magnetically confined wind shock sources: Tr16-22, O8.5 V, and LS 1865, O8.5 V((f)). The O4 III(fc) star HD 93250 exhibits strong, hard, variable X-rays, suggesting it may be a massive binary with a period of >30 days. The visual O2 If* binary HD 93129A shows soft 0.6 keV and hard 1.9 keV emission components, suggesting embedded wind shocks close to the O2 If* Aa primary, and colliding wind shocks between Aa and Ab. Of the 11 known O-type spectroscopic binaries, the long orbital-period systems HD 93343, HD 93403 and QZ Car have higher shock temperatures than short-period systems such as HD 93205 and FO 15. Although the X-rays from most B stars may be produced in the coronae of unseen, low-mass pre-main-sequence companions, a dozen B stars with high Lx cannot be explained by a distribution of unseen companions. One of these, SS73 24 in the Treasure Chest cluster, is a new candidate Herbig Be star.

The distribution of young stars found in the Chandra Carina Complex Project (CCCP) is examined for clustering structure. X-ray surveys are advantageous for identifying young stellar populations compared to optical and infrared surveys in suffering less contamination from nebular emission and Galactic field stars. The analysis is based on smoothed maps of a spatially complete subsample of about 3000 brighter X-ray sources classified as Carina members, and about 10,000 stars from the full CCCP sample. The principal known clusters are recovered, and some additional smaller groups are identified. No rich embedded clusters are present, although a number of sparse groups are found. The CCCP reveals considerable complexity in clustering properties. The Trumpler 14 and 15 clusters have rich stellar populations in unimodal, centrally concentrated structures several parsecs across. Non-spherical internal structure is seen, and large-scale low surface density distributions surround these rich clusters. Trumpler 16, in contrast, is comprised of several smaller clusters within a circular boundary. Collinder 228 is a third type of cluster which extends over tens of parsecs with many sparse compact groups likely arising from triggered star formation processes. A widely dispersed, but highly populous, distribution of X-ray stars across the about 50 pc CCCP mosaic supports a model of past generations of star formation in the region. Collinder 234, a group of massive stars without an associated cluster of pre-main sequence stars, may be part of this dispersed population.

We report the results of a new survey of massive, OB stars throughout the Carina Nebula using the X-ray point source catalog provided by the Chandra Carina Complex Project (CCCP) in conjunction with infrared (IR) photometry from the Two Micron All-Sky Survey and the Spitzer Space Telescope Vela--Carina survey. Mid-IR photometry is relatively unaffected by extinction, hence it provides strong constraints on the luminosities of OB stars, assuming that their association with the Carina Nebula, and hence their distance, is confirmed. We fit model stellar atmospheres to the optical (UBV) and IR spectral energy distributions (SEDs) of 182 OB stars with known spectral types and measure the bolometric luminosity and extinction for each star. We find that the extinction law measured toward the OB stars has two components: Av=1--1.5 mag produced by foreground dust with a ratio of total-to-selective absorption Rv=3.1 plus a contribution from local dust with Rv>4.0 in the Carina molecular clouds that increases as Av increases. Using X-ray emission as a strong indicator of association with Carina, we identify 94 candidate OB stars with Lbol\geq10^4 Lsun by fitting their IR SEDs. If the candidate OB stars are eventually confirmed by follow-up spectroscopic observations, the number of cataloged OB stars in the Carina Nebula will increase by ~50%. Correcting for incompleteness due to OB stars falling below the Lbol cutoff or the CCCP detection limit, these results potentially double the size of the young massive stellar population.

The Chandra Carina Complex Project (CCCP) provides a sensitive X-ray survey of a nearby starburst region over >1 square degree in extent. Thousands of faint X-ray sources are found, many concentrated into rich young stellar clusters. However, significant contamination from unrelated Galactic and extragalactic sources is present in the X-ray catalog. We describe the use of a naive Bayes classifier to assign membership probabilities to individual sources, based on source location, X-ray properties, and visual/infrared properties. For the particular membership decision rule adopted, 75% of CCCP sources are classified as members, 11% are classified as contaminants, and 14% remain unclassified. The resulting sample of stars likely to be Carina members is used in several other studies, which appear in a Special Issue of the ApJS devoted to the CCCP.

We describe detailed simulations of X-ray-emitting populations to evaluate the levels of contamination by both Galactic and extragalactic X-ray sources unrelated to a star-forming region under study. For Galactic contaminations, we consider contribution from main-sequence stars and giants (not including cataclysmic variables and other classes of accretion-driven X-ray binary systems) as they make the dominant contribution at the position of the Carina Nebula. The simulations take into consideration a variety of technical factors involving a Galactic population synthesis model, stellar X-ray luminosity functions, Chandra telescope response, source detection methodology, and possible spatial variations in the X-ray background and absorption through molecular clouds. When applied to the 1.42 square-degree field of the Chandra Carina Complex Project (CCCP), the simulations predict ~5000 contaminating sources (1 source per square arcminute of the survey), evenly distributed across the field. The results of the simulations are further employed in a companion CCCP study to assign membership probabilities to individual sources.

The Great Nebula in Carina provides an exceptional view into the violent massive star formation and feedback that typifies giant HII regions and starburst galaxies. We have mapped the Carina star-forming complex in X-rays, using archival Chandra data and a mosaic of 20 new 60ks pointings using the Chandra X-ray Observatory's Advanced CCD Imaging Spectrometer, as a testbed for understanding recent and ongoing star formation and to probe Carina's regions of bright diffuse X-ray emission. This study has yielded a catalog of properties of >14,000 X-ray point sources; >9800 of them have multiwavelength counterparts. Using Chandra's unsurpassed X-ray spatial resolution, we have separated these point sources from the extensive, spatially-complex diffuse emission that pervades the region; X-ray properties of this diffuse emission suggest that it traces feedback from Carina's massive stars. In this introductory paper, we motivate the survey design, describe the Chandra observations, and present some simple results, providing a foundation for the 15 papers that follow in this Special Issue and that present detailed catalogs, methods, and science results.

We report the first results of imaging the Carina Nebula with Spitzer/IRAC, providing a catalog of point sources and YSOs based on SED fits. We discuss several aspects of the extended emission, including dust pillars that result when a clumpy molecular cloud is shredded by massive star feedback. There are few "extended green objects" (EGOs) normally taken as signposts of outflow activity, and none of the HH jets detected optically are seen as EGOs. A population of "extended red objects" tends to be found around OB stars, some with clear bow-shocks. These are dusty shocks where stellar winds collide with flows off nearby clouds. Finally, the relative distributions of O stars and subclusters of YSOs as compared to dust pillars shows that while some YSOs are located within pillars, many more stars and YSOs reside just outside pillar heads. We suggest that pillars are transient phenomena, part of a continuous outwardly propagating wave of star formation driven by massive star feedback. As pillars are destroyed, they leave newly formed stars in their wake, which are then subsumed into the young OB association. Altogether, the current generation of YSOs shows no strong deviation from a normal IMF. The number of YSOs suggests a roughly constant star-formation rate over the past 3Myr, implying that star formation in pillars constitutes an important mechanism to construct unbound OB associations. Accelerated pillars may give birth to O-type stars that, after several Myr, could appear to have formed in isolation.

Through analysis of archival images and photometry from the Spitzer GLIMPSE and MIPSGAL surveys combined with 2MASS and MSX data, we have identified 488 candidate young stellar objects (YSOs) in the giant molecular cloud M17 SWex, which extends ~50 pc southwest from the prominent Galactic H II region M17. Our sample includes >200 YSOs with masses >3 Msun that will become B-type stars on the main sequence. Extrapolating over the stellar initial mass function (IMF), we find that M17 SWex contains >1.3x10^4 young stars, representing a proto-OB association. The YSO mass function is significantly steeper than the Salpeter IMF, and early O stars are conspicuously absent from M17 SWex. Assuming M17 SWex will form an OB association with a Salpeter IMF, these results reveal the combined effects of (1) more rapid circumstellar disk evolution in more massive YSOs and (2) delayed onset of massive star formation.

We present 737 candidate Young Stellar Objects (YSOs) near the W51 Giant Molecular Cloud (GMC) over an area of 1.25 deg x 1.00 deg selected from Spitzer Space Telescope data. We use spectral energy distribution (SED) fits to identify YSOs and distinguish them from main-sequence or red giant stars, asymptotic giant branch stars, and background galaxies. Based on extinction of each YSO, we separate a total of 437 YSOs associated with the W51 region from the possible foreground sources. We identify 69 highly embedded Stage 0/I candidate YSOs in our field with masses > 5 Msun (corresponding to mid-to early-B main-sequence spectral types) 46 of which are located in the central active star forming regions of W51A and W51B. From the YSOs associated with W51, we find evidence for mass segregation showing that the most massive YSOs are concentrated on the W51 HII region complex. We find a variation in the spatial distribution of the mass function (MF) of YSOs in the mass range between 5 Msun and 18 Msun. The derived slopes of the MF are -1.26 and -2.36 in the active star-forming region and the outer region, respectively. The variation of the MF for YSOs embedded in the molecular cloud implies that the distribution of stellar masses in clusters depends on the local conditions in the parent molecular cloud.

M17 is one of the youngest and most massive nearby star-formation regions in the Galaxy. It features a bright H II region erupting as a blister from the side of a giant molecular cloud (GMC). Combining photometry from the Spitzer GLIMPSE survey with complementary infrared (IR) surveys, we identify candidate young stellar objects (YSOs) throughout a 1.5 deg x 1 deg field that includes the M17 complex. The long sightline through the Galaxy behind M17 creates significant contamination in our YSO sample from unassociated sources with similar IR colors. Removing contaminants, we produce a highly-reliable catalog of 96 candidate YSOs with a high probability of association with the M17 complex. We fit model spectral energy distributions to these sources and constrain their physical properties. Extrapolating the mass function of 62 intermediate-mass YSOs (M >3 Msun), we estimate that >1000 stars are in the process of forming in the extended outer regions of M17. From IR survey images from IRAS and GLIMPSE, we find that M17 lies on the rim of a large shell structure ~0.5 deg in diameter (~20 pc at 2.1 kpc). We present new maps of CO and 13CO (J=2-1) emission, which show that the shell is a coherent, kinematic structure associated with M17 at v = 19 km/s. The shell is an extended bubble outlining the photodissociation region of a faint, diffuse H II region several Myr old. We provide evidence that massive star formation has been triggered by the expansion of the bubble. The formation of the massive cluster ionizing the M17 H II region itself may have been similarly triggered. We conclude that the star formation history in the extended environment of M17 has been punctuated by successive waves of massive star formation propagating through a GMC complex.

We present an analysis of late-O/early-B-powered, parsec-sized bubbles and associated star-formation using 2MASS, GLIMPSE, MIPSGAL and MAGPIS surveys. Three bubbles were selected from the Churchwell et al. (2007) catalog. We confirm that the structure identified in Watson et al. (2008) holds in less energetic bubbles, i.e. a PDR, identified by 8 um emission due to PAHs surrounds hot dust, identified by 24 um emission and ionized gas, identified by 20 cm continuum. We estimate the dynamical age of two bubbles by comparing bubble sizes to numerical models of Hosokawa & Inutsuka (2006). We also identify and analyze candidate young stellar objects (YSOs) using SED fitting and identify sites of possible triggered star-formation. Lastly, we identify likely ionizing sources for two sources based on SED fitting.

Using images from the Spitzer GLIMPSE Legacy survey, we have identified more than 300 extended 4.5 micron sources (abbreviated EGO, Extended Green Object, for the common coding of the [4.5] band as green in 3-color composite IRAC images). We present a catalog of these EGOs, including integrated flux density measurements at 3.6, 4.5, 5.8, 8.0, and 24 microns from the GLIMPSE and MIPSGAL surveys. The average angular separation between a source in our sample and the nearest IRAS point source is >1 arcminute. The majority of EGOs are associated with infrared dark clouds (IRDCs), and where high-resolution 6.7 GHz methanol maser surveys overlap the GLIMPSE coverage, EGOs and 6.7 GHz methanol masers are strongly correlated. Extended 4.5 micron emission is thought to trace shocked molecular gas in protostellar outflows; the association of EGOs with IRDCs and 6.7 GHz methanol masers suggests that extended 4.5 micron emission may pinpoint outflows specifically from massive protostars. The mid-infrared colors of EGOs lie in regions of color-color space occupied by young protostars still embedded in infalling envelopes.

We present a highly reliable flux-limited census of 18,949 point sources in the Galactic mid-plane that have intrinsically red mid-infrared colors. These sources were selected from the Spitzer Space Telescope GLIMPSE I and II surveys of 274 deg^2 of the Galactic mid-plane, and consist mostly of high- and intermediate- mass young stellar objects (YSOs) and asymptotic giant branch (AGB) stars. The selection criteria were carefully chosen to minimize the effects of position- dependent sensitivity, saturation, and confusion. The distribution of sources on the sky and their location in IRAC and MIPS 24 microns color-magnitude and color-color space are presented. Using this large sample, we find that YSOs and AGB stars can be mostly separated by simple color-magnitude selection criteria into approximately 50-70% of YSOs and 30-50% of AGB stars. Planetary nebulae and background galaxies together represent at most 2-3% of all the red sources. 1,004 red sources in the GLIMPSE II region, mostly AGB stars with high mass- loss rates, show significant (>0.3 mag) variability at 4.5 and/or 8.0 microns. With over 11,000 likely YSOs and over 7,000 likely AGB stars, this is to date the largest uniform census of AGB stars and high- and intermediate mass YSOs in the Milky-Way Galaxy.

We report the discovery of six infrared stellar-wind bowshocks in the Galactic massive star formation regions M17 and RCW49 from Spitzer GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire) images. The InfraRed Array Camera (IRAC) on the Spitzer Space Telescope clearly resolves the arc-shaped emission produced by the bowshocks. We combine Two Micron All-Sky Survey (2MASS), Spitzer, MSX, and IRAS observations to obtain the spectral energy distributions (SEDs) of the bowshocks and their individual driving stars. We use the stellar SEDs to estimate the spectral types of the three newly-identified O stars in RCW49 and one previously undiscovered O star in M17. One of the bowshocks in RCW49 reveals the presence of a large-scale flow of gas escaping the H II region at a few 10^2 km/s. Radiation-transfer modeling of the steep rise in the SED of this bowshock toward longer mid-infrared wavelengths indicates that the emission is coming principally from dust heated by the star driving the shock. The other 5 bowshocks occur where the stellar winds of O stars sweep up dust in the expanding H II regions.

We present an analysis of wind-blown, parsec-sized, mid-infrared bubbles and associated star-formation using GLIMPSE/IRAC, MIPSGAL/MIPS and MAGPIS/VLA surveys. Three bubbles from the Churchwell et al. (2006) catalog were selected. The relative distribution of the ionized gas (based on 20 cm emission), PAH emission (based on 8 um, 5.8 um and lack of 4.5 um emission) and hot dust (24 um emission) are compared. At the center of each bubble there is a region containing ionized gas and hot dust, surrounded by PAHs. We identify the likely source(s) of the stellar wind and ionizing flux producing each bubble based upon SED fitting to numerical hot stellar photosphere models. Candidate YSOs are also identified using SED fitting, including several sites of possible triggered star formation.

We have observed the central region of the IR-dark cloud filament associated with IRAS 18507+0121 at millimeter wavelengths in CO(1-0), 13CO(1-0), and C18O(1-0) line emission and with the Spitzer Space Telescope at mid-IR wavelengths. Five massive outflows from two cloud cores were discovered. Three outflows are centered on or near an Ultracompact HII region (G34.4+0.23) while the remaining two outflows originate from the millimeter core G34.4+0.23 MM. Modeling of the SEDs of the mid-IR sources identified 31 young stellar objects in the filament with a combined stellar mass of ~127 +/- 27 Msun. An additional 22 sources were identified as probable cluster members based on the presence of strong 24 micron emission. The total star formation efficiency in the G34.4 cloud filament is estimated to be ~7% while the massive and intermediate mass star formation efficiency in the entire cloud filament is estimated to be roughly 2%. A comparison of the gravitational binding energy with the outflow kinetic energy suggests that the compact core containing G34.4+0.23 MM is being destroyed by its molecular outflows whereas the outflows associated with more massive core surrounding the G34.4 UC HII region are not likely to totally disrupt the cloud. Additionally, a qualitative evaluation of the region appears to suggest that stars in this region may have formed in two stages: first lower mass stars formed and then, a few Myrs later, the more massive stars began to form.

We have identified 230 Tycho-2 Spectral Catalog stars that exhibit 8 micron mid-infrared extraphotospheric excesses in the MidCourse Space Experiment (MSX) and Spitzer Space Telescope Galactic Legacy MidPlane Survey Extraordinaire (GLIMPSE) surveys. Of these, 183 are either OB stars earlier than B8 in which the excess plausibly arises from a thermal bremsstrahlung component or evolved stars in which the excess may be explained by an atmospheric dust component. The remaining 47 stars have spectral classifications B8 or later and appear to be main sequence or late pre-main-sequence objects harboring circumstellar disks. Six of the 47 stars exhibit multiple signatures characteristic of pre-main-sequence circumstellar disks, including emission lines, near-infrared K-band excesses, and X-ray emission. Approximately one-third of the remaining 41 sources have emission lines suggesting relative youth. Of the 25 GLIMPSE stars with SST data at >24 microns, 20 also show an excess at 24 microns. Three additional objects have 24 micron upper limits consistent with possible excesses, and two objects have photospheric measurements at 24 microns. Six MSX sources had a measurement at wavelengths >8 microns. We modeled the excesses in 26 stars having two or more measurements in excess of the expected photospheres as single-component blackbodies. We determine probable disk temperatures and fractional infrared luminosities in the range 191 < T < 787 and 3.9x10^-4 < L_IR/L_* < 2.7x10^-1. We estimate a lower limit on the fraction of Tycho-2 Spectral Catalog main-sequence stars having mid-IR, but not near-IR, excesses to be 1.0+-0.3%.

We present spectroscopic observations of the short-period cataclysmic variable SW Ursa Majoris, obtained by the Far Ultraviolet Spectroscopic Explorer (FUSE) satellite while the system was in quiescence. The data include the resonance lines of O VI at 1031.91 and 1037.61 A. These lines are present in emission, and they exhibit both narrow (~ 150 km/s) and broad (~ 2000 km/s) components. The narrow O VI emission lines exhibit unusual double-peaked and redshifted profiles. We attribute the source of this emission to a cooling flow onto the surface of the white dwarf primary. The broad O VI emission most likely originates in a thin, photoionized surface layer on the accretion disk. We searched for emission from H_2 at 1050 and 1100 A, motivated by the expectation that the bulk of the quiescent accretion disk is in the form of cool, molecular gas. If H_2 is present, then our limits on the fluxes of the H_2 lines are consistent with the presence of a surface layer of atomic H that shields the interior of the disk. These results may indicate that accretion operates primarily in the surface layers of the disk in SW UMa. We also investigate the far-UV continuum of SW UMa and place an upper limit of 15,000 K on the effective temperature of the white dwarf.

We present Hubble Space Telescope (HST) Near-Infrared Camera and Multi-object Spectrometer (NICMOS) observations of the reflection nebulosity associated with the T Tauri star HH 30. The images show the scattered light pattern characteristic of a highly inclined, optically thick disk with a prominent dustlane whose width decreases with increasing wavelength. The reflected nebulosity exhibits a lateral asymmetry in the upper lobe on the opposite side to that reported in previously published Wide Field Planetary Camera 2 (WFPC2) images. The radiation transfer model which most closely reproduces the data has a flared accretion disk with dust grains larger than standard interstellar medium grains by a factor of approximately 2.1. A single hotspot on the stellar surface provides the necessary asymmetry to fit the images and is consistent with previous modeling of the light curve and images. Photometric analysis results in an estimated extinction of Av>~80; however, since the photometry measures only scattered light rather than direct stellar flux, this a lower limit. The radiative transfer models require an extinction of Av = 7,900.