We present multi-wavelength images observed with SOFIA-FORCAST from ~10 to 40 $\mu$m of 14 protostars, selected as intermediate-mass protostar candidates, as part of the SOFIA Massive (SOMA) Star Formation Survey. We build protostellar spectral energy distributions (SEDs) with the SOFIA observations, together with archival data from Spitzer, Herschel and IRAS. We then fit the SEDs with radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, to estimate key properties of the protostars. The SEDs generally indicate the validity of these RT models down to intermediate-mass and/or early-stage protostars. The protostars analyzed so far in the SOMA survey span a range of luminosities from ~$10^{2}$ to ~$10^{6} L_{\odot}$, current protostellar masses from ~0.5 to ~30 $M_{\odot}$, and ambient clump mass surface densities, $\Sigma_{\rm cl}$ of 0.1 - 3 g cm$^{-2}$. A wide range of evolutionary states of the individual protostars and of the protocluster environments are also probed. The 19 to 37 $\mu$m spectral index of the sources correlates with outflow cavity opening angle, ratio of this angle to viewing angle, and evolutionary stage. We have also added a sample of ~50 protostellar sources identified from within Infrared Dark Clouds and expected to be at the earliest stages of their evolution. With this global sample, most of the evolutionary stages of high- and intermediate-mass protostars are probed. From the best fitting models of the protostars, there is no evidence of a threshold value of protocluster clump mass surface density being needed to form protostars up to ~25 $M_\odot$. However, to form more massive protostars, there is tentative evidence that $\Sigma_{\rm cl}$ needs to be at least 1 g cm$^{-2}$. We discuss how this is consistent with expectations from core accretion models that include internal feedback from the forming massive star.

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

We study a very young star-forming region in the outer Galaxy that is the most concentrated source of outflows in the Spitzer Space Telescope GLIMPSE360 survey. This region, dubbed CMa-l224, is located in the Canis Major OB1 association. CMa-l224 is relatively faint in the mid-infrared, but it shines brightly at the far-infrared wavelengths as revealed by the Herschel Space Observatory data from the Hi-GAL survey. Using the 3.6 and 4.5 $\mu$m data from the Spitzer/GLIMPSE360 survey, combined with the JHK$_s$ 2MASS and the 70-500 $\mu$m Herschel/Hi-GAL data, we develop a young stellar object (YSO) selection criteria based on color-color cuts and fitting of the YSO candidates' spectral energy distributions with YSO 2D radiative transfer models. We identify 293 YSO candidates and estimate physical parameters for 210 sources well-fit with YSO models. We select an additional 47 sources with GLIMPSE360-only photometry as `possible YSO candidates'. The vast majority of these sources are associated with high H$_2$ column density regions and are good targets for follow-up studies. The distribution of YSO candidates at different evolutionary stages with respect to Herschel filaments supports the idea that stars are formed in the filaments and become more dispersed with time. Both the supernova-induced and spontaneous star formation scenarios are plausible in the environmental context of CMa-l224. However, our results indicate that a spontaneous gravitational collapse of filaments is a more likely scenario. The methods developed for CMa-l224 can be used for larger regions in the Galactic plane where the same set of photometry is available.

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

We present radiation transfer models of rotating young stellar objects (YSOs) with hotspots in their atmospheres, inner disk warps and other 3-D effects in the nearby circumstellar environment. Our models are based on the geometry expected from the magneto-accretion theory, where material moving inward in the disk flows along magnetic field lines to the star and creates stellar hotspots upon impact. Due to rotation of the star and magnetosphere, the disk is variably illuminated. We compare our model light curves to data from the Spitzer YSOVAR project (Morales-Calderon et al. 2014, Cody et al. 2014) to determine if these processes can explain the variability observed at optical and mid-infrared wavelengths in young stars. We focus on those variables exhibiting "dipper" behavior that may be periodic, quasi-periodic, or aperiodic. We find that the stellar hotspot size and temperature affects the optical and near-infrared light curves, while the shape and vertical extent of the inner disk warp affects the mid-IR light curve variations. Clumpy disk distributions with non-uniform fractal density structure produce more stochastic light curves. We conclude that the magneto-accretion theory is consistent with certain aspects of the multi-wavelength photometric variability exhibited by low-mass YSOs. More detailed modeling of individual sources can be used to better determine the stellar hotspot and inner disk geometries of particular sources.

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.

It has been suggested that the gaps and cavities recently discovered in transitional disks are opened by planets. To explore this scenario, we combine two-dimensional two fluid (gas + particle) hydrodynamical calculations with three-dimensional Monte Carlo Radiative Transfer simulations, and study the observational signatures of gaps opened by one or several planets, making qualitative comparisons with observations. We find that a single planet as small as 0.2 MJ can produce a deep gap at millimeter (mm) wavelengths and almost no features at near-infrared (NIR) wavelengths, while multiple planets can open up a few *10 AU wide common gap at both wavelengths. Both the contrast ratio of the gaps and the wavelength dependence of the gap sizes are broadly consistent with data. We also confirm previous results that NIR gap sizes may be smaller than mm gap sizes due to dust-gas coupling and radiative transfer effects. When viewed at a moderate inclination angle, a physically circular on-centered gap could appear to be off-centered from the star due to shadowing. Planet-induced spiral arms are more apparent at NIR than at mm wavelengths. Overall, our results suggest that the planet-opening-gap scenario is a promising way to explain the origin of the transitional disks. Finally, inspired by the recent ALMA release of the image of the HL Tau disk, we show that multiple narrow gaps, well separated by bright rings, can be opened by 0.2 MJ planets soon after their formation in a relatively massive disk.

We present the first results from the SAGE-Var program, a follow on to the Spitzer legacy program Surveying the Agents of Galaxy Evolution (SAGE; Meixner, et al. 2006). We obtained 4 epochs of photometry at 3.6 & 4.5 microns covering the bar of the Large Magellanic Cloud (LMC) and the central region of the Small Magellanic Cloud (SMC) in order to probe the variability of extremely red sources missed by variability surveys conducted at shorter wavelengths, and to provide additional epochs of observation for known variables. Our 6 total epochs of observations allow us to probe infrared variability on 15 different timescales ranging from ~20 days to ~5 years. Out of a full catalog of 1,717,554 (LMC) and 457,760 (SMC) objects, we find 10 (LMC) and 6 (SMC) large amplitude AGB variables without optically measured variability owing to circumstellar dust obscuration. The catalog also contains multiple observations of known AGB variables, type I and II Cepheids, eclipsing variables, R CrB stars and young stellar objects which will be discussed in following papers. Here we present infrared Period-Luminosity (PL) relations for classical Cepheids in the Magellanic Clouds, as well as improved PL relationships for AGB stars pulsating in the fundamental mode using mean magnitudes constructed from 6 epochs of observations.

We present near-IR and far-UV observations of the pre-transitional (gapped) disk in HD 169142 using NASA's Infrared Telescope Facility and Hubble Space Telescope. The combination of our data along with existing data sets into the broadband spectral energy distribution reveals variability of up to 45% between ~1.5-10 {\mu}m over a maximum timescale of 10 years. All observations known to us separate into two distinct states corresponding to a high near-IR state in the pre-2000 epoch and a low state in the post-2000 epoch, indicating activity within the <1 AU region of the disk. Through analysis of the Pa {\beta} and Br {\gamma} lines in our data we derive a mass accretion rate in May 2013 of (1.5 - 2.7) x 10^-9 Msun/yr. We present a theoretical modeling analysis of the disk in HD 169142 using Monte-Carlo radiative transfer simulation software to explore the conditions and perhaps signs of planetary formation in our collection of 24 years of observations. We find that shifting the outer edge (r = 0.3 AU) of the inner disk by 0.05 AU toward the star (in simulation of accretion and/or sculpting by forming planets) successfully reproduces the shift in NIR flux. We establish that the ~40-70 AU dark ring imaged in the NIR by Quanz et al. (2013) and Momose et al. (2013) and at 7 mm by Osorio et al. (2014) may be reproduced with a 30% scaled density profile throughout the region, strengthening the link to this structure being dynamically cleared by one or more planetary mass bodies.

Among the most fascinating and hotly-debated areas in contemporary astrophysics are the means by which planetary systems are assembled from the large rotating disks of gas and dust which attend a stellar birth. Although important work has already been, and is still being done both in theory and observation, a full understanding of the physics of planet formation can only be achieved by opening observational windows able to directly witness the process in action. The key requirement is then to probe planet-forming systems at the natural spatial scales over which material is being assembled. By definition, this is the so-called Hill Sphere which delineates the region of influence of a gravitating body within its surrounding environment. The Planet Formation Imager project (PFI) has crystallized around this challenging goal: to deliver resolved images of Hill-Sphere-sized structures within candidate planet-hosting disks in the nearest star-forming regions. In this contribution we outline the primary science case of PFI. For this purpose, we briefly review our knowledge about the planet-formation process and discuss recent observational results that have been obtained on the class of transition disks. Spectro-photometric and multi-wavelength interferometric studies of these systems revealed the presence of extended gaps and complex density inhomogeneities that might be triggered by orbiting planets. We present detailed 3-D radiation-hydrodynamic simulations of disks with single and multiple embedded planets, from which we compute synthetic images at near-infrared, mid-infrared, far-infrared, and sub-millimeter wavelengths, enabling a direct comparison of the signatures that are detectable with PFI and complementary facilities such as ALMA. From these simulations, we derive some preliminary specifications that will guide the array design and technology roadmap of the facility.

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30-day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes, and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" having discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with H$\alpha$ emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk.

The Magellanic Bridge is the nearest low-metallicity, tidally stripped environment, offering a unique high-resolution view of physical conditions in merging and forming galaxies. In this paper we present analysis of candidate massive young stellar objects (YSOs), i.e., {\it in situ, current} massive star formation (MSF) in the Bridge using {\it Spitzer} mid-IR and complementary optical and near-IR photometry. While we definitely find YSOs in the Bridge, the most massive are $\sim10 M_\odot$, $\ll45 M_\odot$ found in the Large Magellanic Cloud (LMC). The intensity of MSF in the Bridge also appears decreasing, as the most massive YSOs are less massive than those formed in the past. To investigate environmental effects on MSF, we have compared properties of massive YSOs in the Bridge to those in the LMC. First, YSOs in the Bridge are apparently less embedded than in the LMC: 81% of Bridge YSOs show optical counterparts, compared to only 56% of LMC sources with the same range of mass, circumstellar dust mass, and line-of-sight extinction. Circumstellar envelopes are evidently more porous or clumpy in the Bridge's low-metallicity environment. Second, we have used whole samples of YSOs in the LMC and the Bridge to estimate the probability of finding YSOs at a given \hi\ column density, N(HI). We found that the LMC has $\sim3\times$ higher probability than the Bridge for N(HI) $>10\times10^{20}$ cm$^{-2}$, but the trend reverses at lower N(HI). Investigating whether this lower efficiency relative to HI is due to less efficient molecular cloud formation, or less efficient cloud collapse, or both, will require sensitive molecular gas observations.

Based on more than four weeks of continuous high cadence photometric monitoring of several hundred members of the young cluster NGC 2264 with two space telescopes, NASA's Spitzer and the CNES CoRoT (Convection, Rotation, and planetary Transits), we provide high quality, multi-wavelength light curves for young stellar objects (YSOs) whose optical variability is dominated by short duration flux bursts, which we infer are due to enhanced mass accretion rates. These light curves show many brief -- several hour to one day -- brightenings at optical and near-infrared (IR) wavelengths with amplitudes generally in the range 5-50% of the quiescent value. Typically, a dozen or more of these bursts occur in a thirty day period. We demonstrate that stars exhibiting this type of variability have large ultraviolet (UV) excesses and dominate the portion of the u-g vs. g-r color-color diagram with the largest UV excesses. These stars also have large Halpha equivalent widths, and either centrally peaked, lumpy Halpha emission profiles or profiles with blue-shifted absorption dips associated with disk or stellar winds. Light curves of this type have been predicted for stars whose accretion is dominated by Rayleigh-Taylor instabilities at the boundary between their magnetosphere and inner circumstellar disk, or where magneto-rotational instabilities modulate the accretion rate from the inner disk. Amongst the stars with the largest UV excesses or largest Halpha equivalent widths, light curves with this type of variability greatly outnumber light curves with relatively smooth sinusoidal variations associated with long-lived hot spots. We provide quantitative statistics for the average duration and strength of the accretion bursts and for the fraction of the accretion luminosity associated with these bursts.

Massive young stellar objects (YSOs), like low-mass YSOs, are thought to be surrounded by optically thick envelopes and/or discs and are observed to have associated regions that produce polarized light at near-infrared wavelengths. These polarized regions are thought to be lower-density outflows along the polar axes of the YSO envelopes. Using the 0.2 arcsec spatial resolution of the Near-Infrared Camera and Multi-Object Spectrometer on the Hubble Space Telescope we are examining the structure of the envelopes and outflow regions of massive YSOs in star-forming regions within a few kpc of the Sun. Here we report on 2 micron polarimetry of Mon R2-IRS3, S140-IRS1, and AFGL 2591. All three sources contain YSOs with highly-polarized monopolar outflows, with Mon R2-IRS3 containing at least two YSOs in a small cluster. The central stars of all four YSOs are also polarized, with position angles perpendicular to the directions of the outflows. We infer that this polarization is due to scattering and absorption by aligned grains. We have modelled our observations of S140-IRS1 and AFGL 2591 as light scattered and absorbed both by spherical grains and by elongated grains that are aligned by magnetic fields. Models that best reproduce the observations have a substantial toroidal component to the magnetic field in the equatorial plane. Moreover, the toroidal magnetic field in the model that best fits AFGL 2591 extends a large fraction of the height of the model cavity, which is 10^5 au. We conclude that the massive YSOs in this study all show evidence of the presence of a substantial toroidal magnetic field.

We have updated our publicly available dust radiative transfer code (HOCHUNK3D) to include new emission processes and various 3-D geometries appropriate for forming stars. The 3-D geometries include warps and spirals in disks, accretion hotspots on the central star, fractal clumping density enhancements, and misaligned inner disks. Additional axisymmetric (2-D) features include gaps in disks and envelopes, "puffed-up inner rims" in disks, multiple bipolar cavity walls, and iteration of disk vertical structure assuming hydrostatic equilibrium. We include the option for simple power-law envelope geometry, which combined with fractal clumping, and bipolar cavities, can be used to model evolved stars as well as protostars. We include non-thermal emission from PAHs and very small grains, and external illumination from the interstellar radiation field. The grid structure was modified to allow multiple dust species in each cell; based on this, a simple prescription is implemented to model dust stratification. We describe these features in detail, and show example calculations of each. Some of the more interesting results include the following: 1) Outflow cavities may be more clumpy than infalling envelopes. 2) PAH emission in high-mass stars may be a better indicator of evolutionary stage than the broadband SED slope; and related to this, 3) externally illuminated clumps and high-mass stars in optically thin clouds can masquerade as YSOs. 4) Our hydrostatic equilibrium models suggest that dust settling is likely ubiquitous in T Tauri disks, in agreement with previous observations.

We characterize the completeness of point source lists from Spitzer Space Telescope surveys in the four Infrared Array Camera (IRAC) bandpasses, emphasizing the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) programs (GLIMPSE I, II, 3D, 360; Deep GLIMPSE) and their resulting point source Catalogs and Archives. The analysis separately addresses effects of incompleteness resulting from high diffuse background emission and incompleteness resulting from point source confusion (i.e., crowding). An artificial star addition and extraction analysis demonstrates that completeness is strongly dependent on local background brightness and structure, with high-surface-brightness regions suffering up to five magnitudes of reduced sensitivity to point sources. This effect is most pronounced at the IRAC 5.8 and 8.0 microns bands where UV-excited PAH emission produces bright, complex structures (photodissociation regions; PDRs). With regard to diffuse background effects, we provide the completeness as a function of stellar magnitude and diffuse background level in graphical and tabular formats. These data are suitable for estimating completeness in the low-source-density limit in any of the four IRAC bands in GLIMPSE Catalogs and Archives and some other Spitzer IRAC programs that employ similar observational strategies and are processed by the GLIMPSE pipeline. Point source incompleteness is primarily a consequence of structure in the diffuse background emission rather than photon noise. With regard to source confusion in the high-source-density regions of the Galactic Plane, we provide figures illustrating the 90% completeness levels as a function of point source density at each band. (Slightly abridged)

ALMA observations of the Galactic center with spatial resolution $2.61"\times0.97"$ resulted in the detection of 11 SiO (5-4) clumps of molecular gas within 0.6pc (15$"$) of Sgr A*, interior to the 2-pc circumnuclear molecular ring. The three SiO (5-4) clumps closest to Sgr A* show the largest central velocities, $\sim150$ \kms, and broadest asymmetric linewidths with full width zero intensity (FWZI) $\sim110-147$ \kms. The remaining clumps, distributed mainly to the NE of the ionized mini-spiral, have narrow FWZI ($\sim18-56$ \kms). Using CARMA SiO (2-1) data, LVG modeling of the the SiO line ratios for the broad velocity clumps, constrains the column density N(SiO) $\sim10^{14}$ cm$^{-2}$, and the H$_2$ gas density n$_{\rm H_2}=(3-9)\times10^5$ cm$^{-3}$ for an assumed kinetic temperature 100-200K. The SiO clumps are interpreted as highly embedded protostellar outflows, signifying an early stage of massive star formation near Sgr A* in the last $10^4-10^5$ years. Support for this interpretation is provided by the SiO (5-4) line luminosities and velocity widths which lie in the range measured for protostellar outflows in star forming regions in the Galaxy. Furthermore, SED modeling of stellar sources shows two YSO candidates near SiO clumps, supporting in-situ star formation near Sgr A*. We discuss the nature of star formation where the gravitational potential of the black hole dominates. In particular, we suggest that external radiative pressure exerted on self-shielded molecular clouds enhances the gas density, before the gas cloud become gravitationally unstable near Sgr A*. Alternatively, collisions between clumps in the ring may trigger gravitational collapse.

We present the first near infrared (NIR) spatially resolved images of the circumstellar transitional disk around SR21. These images were obtained with the Subaru HiCIAO camera, adaptive optics and the polarized differential imaging (PDI) technique. We resolve the disk in scattered light at H-band for stellocentric 0.1"<r<0.6" (12<r<75AU). We compare our results with previously published spatially-resolved 880 micron continuum Submillimeter Array (SMA) images that show an inner r<36AU cavity in SR21. Radiative transfer models reveal that the large disk depletion factor invoked to explain SR21's sub-mm cavity cannot be "universal" for all grain sizes. Even significantly more moderate depletions (delta=0.1, 0.01 relative to an undepleted disk) than those that reproduce the sub-mm cavity (delta~10^-6) are inconsistent with our H-band images when they are assumed to carry over to small grains, suggesting that surface grains scattering in the NIR either survive or are generated by whatever mechanism is clearing the disk midplane. In fact, the radial polarized intensity profile of our H-band observations is smooth and steeply inwardly-increasing (r^-3), with no evidence of a break at the 36AU sub-mm cavity wall. We hypothesize that this profile is dominated by an optically thin disk envelope or atmosphere component. We also discuss the compatibility of our data with the previously postulated existence of a sub-stellar companion to SR21 at r~10-20AU, and find that we can neither exclude nor verify this scenario. This study demonstrates the power of multiwavelength imaging of transitional disks to inform modeling efforts, including the debate over precisely what physical mechanism is responsible for clearing these disks of their large midplane grains.

We present 30 and 40 micron imaging of the massive protostar G35.20-0.74 with SOFIA-FORCAST. The high surface density of the natal core around the protostar leads to high extinction, even at these relatively long wavelengths, causing the observed flux to be dominated by that emerging from the near-facing outflow cavity. However, emission from the far-facing cavity is still clearly detected. We combine these results with fluxes from the near-infrared to mm to construct a spectral energy distribution (SED). For isotropic emission the bolometric luminosity would be 3.3x10^4 Lsun. We perform radiative transfer modeling of a protostar forming by ordered, symmetric collapse from a massive core bounded by a clump with high mass surface density, Sigma_cl. To fit the SED requires protostellar masses ~20-34 Msun depending on the outflow cavity opening angle (35 - 50 degrees), and Sigma_cl ~ 0.4-1 g cm-2. After accounting for the foreground extinction and the flashlight effect, the true bolometric luminosity is ~ (0.7-2.2)x10^5 Lsun. One of these models also has excellent agreement with the observed intensity profiles along the outflow axis at 10, 18, 31 and 37 microns. Overall our results support a model of massive star formation involving the relatively ordered, symmetric collapse of a massive, dense core and the launching bipolar outflows that clear low density cavities. Thus a unified model may apply for the formation of both low and high mass stars.

We present a self-consistent three-dimensional Monte-Carlo radiative transfer model of the stellar and dust emission in the Milky-Way, and have computed synthetic observations of the 3.6 to 100 microns emission in the Galactic mid-plane. In order to compare the model to observations, we use the GLIMPSE, MIPSGAL, and IRAS surveys to construct total emission spectra, as well as longitude and latitude profiles for the emission. The distribution of stars and dust is taken from the SKY model, and the dust emissivities includes an approximation of the emission from polycyclic aromatic hydrocarbons in addition to thermal emission. The model emission is in broad agreement with the observations, but a few modifications are needed to obtain a good fit. Firstly, by adjusting the model to include two major and two minor spiral arms rather than four equal spiral arms, the fit to the longitude profiles for |l|>30 degrees can be improved. Secondly, introducing a deficit in the dust distribution in the inner Galaxy results in a better fit to the shape of the IRAS longitude profiles at 60 and 100 microns. With these modifications, the model fits the observed profiles well, although it systematically under-estimates the 5.8 and 8.0 microns fluxes. One way to resolve this discrepancy is to increase the abundance of PAH molecules by 50% compared to the original model, although we note that changes to the dust distribution or radiation field may provide alternative solutions. Finally, we use the model to quantify which stellar populations contribute the most to the heating of different dust types, and which stellar populations and dust types contribute the most to the emission at different wavelengths.

Through detailed radiative transfer modeling, we present a disk+cavity model to simultaneously explain both the SED and Subaru H-band polarized light imaging for the pre-transitional protoplanetary disk PDS 70. Particularly, we are able to match not only the radial dependence, but also the absolute scale, of the surface brightness of the scattered light. Our disk model has a cavity 65 AU in radius, which is heavily depleted of sub-micron-sized dust grains, and a small residual inner disk which produces a weak but still optically thick NIR excess in the SED. To explain the contrast of the cavity edge in the Subaru image, a factor of ~1000 depletion for the sub-micron-sized dust inside the cavity is required. The total dust mass of the disk may be on the order of 1e-4 M_sun, only weakly constrained due to the lack of long wavelength observations and the uncertainties in the dust model. The scale height of the sub-micron-sized dust is ~6 AU at the cavity edge, and the cavity wall is optically thick in the vertical direction at H-band. PDS 70 is not a member of the class of (pre-)transitional disks identified by Dong et al. (2012), whose members only show evidence of the cavity in the millimeter-sized dust but not the sub-micron-sized dust in resolved images. The two classes of (pre-)transitional disks may form through different mechanisms, or they may just be at different evolution stages in the disk clearing process.

We have obtained time series observations of the Orion Nebula Cluster at 70 microns and 160 microns from the Herschel/PACS Photometer. This represents the first wide-field far-infrared photometric monitoring of a young star forming region. The acquired 35'x35' maps show complex extended structures, with unprecedented details, that trace the interaction between the molecular gas and the young hot stars. We detect 43 protostars, most of which are situated along the integral-shaped filament extending from the Orion nebula, through OMC2 and to OMC3. We present high-reliability light curves for some of these objects using the first six epochs of our observing program spread over 6 weeks. We find amplitude variations in excess of 20% for a fraction of the detected protostars over periods as short as a few weeks. This is inconsistent with the dynamical time-scales of cool far-IR emitting material that orbits at hundreds of AU from the protostar, and it suggests that the mechanism(s) responsible for the observed variability originates from the inner region of the protostars, likely driven by variable mass accretion.

Transitional circumstellar disks around young stellar objects have a distinctive infrared deficit around 10 microns in their Spectral Energy Distributions (SED), recently measured by the Spitzer Infrared Spectrograph (IRS), suggesting dust depletion in the inner regions. These disks have been confirmed to have giant central cavities by imaging of the submillimeter (sub-mm) continuum emission using the Submillimeter Array (SMA). However, the polarized near-infrared scattered light images for most objects in a systematic IRS/SMA cross sample, obtained by HiCIAO on the Subaru telescope, show no evidence for the cavity, in clear contrast with SMA and Spitzer observations. Radiative transfer modeling indicates that many of these scattered light images are consistent with a smooth spatial distribution for micron-sized grains, with little discontinuity in the surface density of the micron-sized grains at the cavity edge. Here we present a generic disk model that can simultaneously account for the general features in IRS, SMA, and Subaru observations. Particularly, the scattered light images for this model are computed, which agree with the general trend seen in Subaru data. Decoupling between the spatial distributions of the micron-sized dust and mm-sized dust inside the cavity is suggested by the model, which, if confirmed, necessitates a mechanism, such as dust filtration, for differentiating the small and big dust in the cavity clearing process. Our model also suggests an inwardly increasing gas-to-dust-ratio in the inner disk, and different spatial distributions for the small dust inside and outside the cavity, echoing the predictions in grain coagulation and growth models.

In 2007, R Coronae Borealis (R CrB) went into an historically deep and long decline. In this state, the dust acts like a natural coronagraph at visible wavelengths, allowing faint nebulosity around the star to be seen. Imaging has been obtained from 0.5 to 500 micron with Gemini/GMOS, HST/WFPC2, Spitzer/MIPS, and Herschel/SPIRE. Several of the structures around R CrB are cometary globules caused by wind from the star streaming past dense blobs. The estimated dust mass of the knots is consistent with their being responsible for the R CrB declines if they form along the line of sight to the star. In addition, there is a large diffuse shell extending up to 4 pc away from the star containing cool 25 K dust that is detected all the way out to 500 micron. The SED of R CrB can be well fit by a 150 AU disk surrounded by a very large diffuse envelope which corresponds to the size of the observed nebulosity. The total masses of the disk and envelope are 10^-4 and 2 M(Sun), respectively, assuming a gas-to-dust ratio of 100. The evidence pointing toward a white-dwarf merger or a final-helium-shell flash origin for R CrB is contradictory. The shell and the cometary knots are consistent with a fossil planetary nebula. Along with the fact that R CrB shows significant Lithium in its atmosphere, this supports the final-helium-shell flash. However, the relatively high inferred mass of R CrB and its high fluorine abundance support a white-dwarf merger.

This is a collection of articles that were originally published in the Newsletter of the American Astronomical Society (AAS) between May 2008 and September 2011 by the Committee on Employment. Authors representing a wide range of career paths tell their stories and provide insight and advice that is relevant to success in various job sectors. Although all of these articles are available individually from the AAS archives, we are posting the complete collection here to make them more accessible as a resource for the astronomy community. The collection includes the following articles: (1) Changing Priorities: the Hard Money Wild Card, (2) Beyond Ivory Towers, (3) Astronomers Working in Public Outreach, (4) Bush-Whacking a Career Trail, (5) Science Communication as a Press Officer, (6) Jobs in Industry, (7) Back to School: A Ph.D. Enters the Classroom, (8) Working at a Soft-Money Institute, (9) Balancing Research and Service at NOAO, (10) Succeeding in a Large Research Collaboration, and (11) Preparing for the College Teaching Job Market. The final published versions of the articles can be found at http://aas.org/career/

The Small Magellanic Cloud (SMC) provides a unique laboratory for the study of the lifecycle of dust given its low metallicity (~1/5 solar) and relative proximity (~60 kpc). This motivated the SAGE-SMC (Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud) Spitzer Legacy program with the specific goals of studying the amount and type of dust in the present interstellar medium, the sources of dust in the winds of evolved stars, and how much dust is consumed in star formation. This program mapped the full SMC (30 sq. deg.) including the Body, Wing, and Tail in 7 bands from 3.6 to 160 micron using the IRAC and MIPS instruments on the Spitzer Space Telescope. The data were reduced, mosaicked, and the point sources measured using customized routines specific for large surveys. We have made the resulting mosaics and point source catalogs available to the community. The infrared colors of the SMC are compared to those of other nearby galaxies and the 8 micron/24 micron ratio is somewhat lower and the 70 micron/160 micron ratio is somewhat higher than the average. The global infrared spectral energy distribution shows that the SMC has ~3X lower aromatic emission/PAH (polycyclic aromatic hydrocarbon) abundances compared to most nearby galaxies. Infrared color-magnitude diagrams are given illustrating the distribution of different asymptotic giant branch stars and the locations of young stellar objects. Finally, the average spectral energy distribution (SED) of HII/star formation regions is compared to the equivalent Large Magellanic Cloud average HII/star formation region SED. These preliminary results are expanded in detail in companion papers.

We have used the Spitzer Space Telescope Infrared Spectrograph (IRS) to observe the 5-37 micron thermal emission of comet 73P/Schwassmann-Wachmann 3 (SW3), components B and C. We obtained low spectral resolution (R ~ 100) data over the entire wavelength interval, along with images at 16 and 22 micron. These observations provided an unprecedented opportunity to study nearly pristine material from the surface and what was until recently the interior of an ecliptic comet - cometary surface having experienced only two prior perihelion passages, and including material that was totally fresh. The spectra were modeled using a variety of mineral types including both amorphous and crystalline components. We find that the degree of silicate crystallinity, ~ 35%, is somewhat lower than most other comets with strong emission features, while its abundance of amorphous carbon is higher. Both suggest that SW3 is among the most chemically primitive solar system objects yet studied in detail, and that it formed earlier or farther from the sun than the bulk of the comets studied so far. The similar dust compositions of the two fragments suggests that these are not mineralogically heterogeneous, but rather uniform throughout their volumes. Atomic abundances derived from the spectral models indicates a depletion of O compared to solar photospheric values, despite the inclusion of water ice and gas in the models. Atomic C may be solar or slightly sub-solar, but its abundance is complicated by the potential contribution of spectrally featureless mineral species to the portion of the spectra most sensitive to the derication of the C abundance. We find a relatively high bolometric albedo, ~ 0.13 for the dust, considering the large amount of dark carbonaceous material, but consistent with the presence of abundant small particles and strong emission features.

We investigate the infrared (IR) properties of cool, evolved stars in the Small Magellanic Cloud (SMC), including the red giant branch (RGB) stars and the dust-producing red supergiant (RSG) and asymptotic giant branch (AGB) stars using observations from the Spitzer Space Telescope Legacy program entitled: "Surveying the Agents of Galaxy Evolution in the Tidally-stripped, Low Metallicity SMC", or SAGE-SMC. The survey includes, for the first time, full spatial coverage of the SMC bar, wing, and tail regions at infrared (IR) wavelengths (3.6 - 160 microns). We identify evolved stars using a combination of near-IR and mid-IR photometry and point out a new feature in the mid-IR color-magnitude diagram that may be due to particularly dusty O-rich AGB stars. We find that the RSG and AGB stars each contribute ~20% of the global SMC flux (extended + point-source) at 3.6 microns, which emphasizes the importance of both stellar types to the integrated flux of distant metal-poor galaxies. The equivalent SAGE survey of the higher-metallicity Large Magellanic Cloud (SAGE-LMC) allows us to explore the influence of metallicity on dust production. We find that the SMC RSG stars are less likely to produce a large amount of dust (as indicated by the [3.6]-[8] color). There is a higher fraction of carbon-rich stars in the SMC, and these stars appear to able to reach colors as red as their LMC counterparts, indicating that C-rich dust forms efficiently in both galaxies. A preliminary estimate of the dust production in AGB and RSG stars reveals that the extreme C-rich AGB stars dominate the dust input in both galaxies, and that the O-rich stars may play a larger role in the LMC than in the SMC.

We present new Spitzer/IRS spectra of two hot R Coronae Borealis (RCB) stars, one in the Galaxy,V348 Sgr, and one lying in the LMC, HV 2671. These two objects may constitute a link between the RCB stars and the late Wolf-Rayet ([WCL]) class of central stars of planetary nebula (CSPNe) such as CPD -56 8032 that has little or no hydrogen in their atmospheres. HV 2671 and V348 Sgr are members of a rare subclass that has significantly higher effective temperatures than most RCB stars, but sharing the traits of hydrogen deficiency and dust formation that define the cooler RCB stars. The [WC] CSPNe star, CPD -56 8032, displays evidence for dual-dust chemistry showing both PAHs and crystalline silicates in its mid-IR spectrum. HV 2671 shows strong PAH emission but shows no sign of having crystalline silicates. The spectrum of V348 Sgr is very different from those of CPD -56 8032 and HV 2671. The PAH emission seen strongly in the other two stars is not present. Instead, the spectrum is dominated by a broad emission centered at about 8.2 micron. The mid-IR spectrum of CPD -56 8032 shows emission features that may be associated with C60. The other two stars do not show evidence for C60. HV 2671 has also been detected by Herschel/PACS and SPIRE. V348 Sgr and CPD -56 8032 have been detected by AKARI/FIS. These data were combined with Spitzer, IRAS, 2MASS and other photometry to produce their spectral energy distributions from the visible to the far-IR. Monte Carlo radiative transfer modeling was used to study the circumstellar dust around these stars. HV 2671 and CPD -56 8032 require both a flared inner disk with warm dust and an extended diffuse envelope with cold dust to to fit their SEDs. The SED of V348 Sgr can be fit with a much smaller disk and envelope.

I outline methods for calculating the solution of Monte Carlo Radiative Transfer (MCRT) in scattering, absorption and emission processes of dust and gas, including polarization. I provide a bibliography of relevant papers on methods with astrophysical applications.

We present initial results from time series imaging at infrared wavelengths of 0.9 sq. degrees in the Orion Nebula Cluster (ONC). During Fall 2009 we obtained 81 epochs of Spitzer 3.6 and 4.5 micron data over 40 consecutive days. We extracted light curves with ~3% photometric accuracy for ~2000 ONC members ranging from several solar masses down to well below the hydrogen burning mass limit. For many of the stars, we also have time-series photometry obtained at optical (Ic) and/or near-infrared (JKs) wavelengths. Our data set can be mined to determine stellar rotation periods, identify new pre-main-sequence (PMS) eclipsing binaries, search for new substellar Orion members, and help better determine the frequency of circumstellar disks as a function of stellar mass in the ONC. Our primary focus is the unique ability of 3.6 & 4.5 micron variability information to improve our understanding of inner disk processes and structure in the Class I and II young stellar objects (YSOs). In this paper, we provide a brief overview of the YSOVAR Orion data obtained in Fall 2009, and we highlight our light curves for AA-Tau analogs - YSOs with narrow dips in flux, most probably due to disk density structures passing through our line of sight. Detailed follow-up observations are needed in order to better quantify the nature of the obscuring bodies and what this implies for the structure of the inner disks of YSOs.

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.

With high-resolution infrared data becoming available that can probe the formation of high-mass stellar clusters for the first time, models that make testable predictions of these objects are necessary. We utilize a three-dimensional radiative transfer code, including a hierarchically clumped medium, to study the earliest stages of super star cluster evolution. We explore a range of parameter space in geometric sequences that mimic the evolution of an embedded super star cluster. The inclusion of a hierarchically clumped medium can make the envelope porous, in accordance with previous models and supporting observational evidence. The infrared luminosity inferred from observations can differ by a factor of two from the true value in the clumpiest envelopes depending on the viewing angle. The infrared spectral energy distribution also varies with viewing angle for clumpy envelopes, creating a range in possible observable infrared colors and magnitudes, silicate feature depths and dust continua. General observable features of cluster evolution differ between envelopes that are relatively opaque or transparent to mid-infrared photons. The [70]-[160] color can be used to determine star formation efficiency; the Spitzer IRAC/MIPS [8.0]-[24] color is able to constrain Rin and Rout values; and the IRAC [3.6]-[5.8] color is sensitive to the fraction of the dust distributed in clumps. Finally, in a comparison of these models to data of ultracompact HII regions, we find good agreement, suggesting that these models are physically relevant, and will provide useful diagnostic ability for datasets of resolved, embedded SSCs with the advent of high-resolution infrared telescopes like JWST.

Fundamental parameters and time-evolution of mass loss are investigated for post-main-sequence stars in the Galactic globular cluster 47 Tucanae (NGC 104). This is accomplished by fitting spectral energy distributions (SEDs) to existing optical and infrared photometry and spectroscopy, to produce a true Hertzsprung--Russell diagram. We confirm the cluster's distance as 4611 (+213, -200) pc and age as 12 +/- 1 Gyr. Horizontal branch models appear to confirm that no more RGB mass loss occurs in 47 Tuc than in the more-metal-poor omega Centauri, though difficulties arise due to inconsistencies between the models. Using our SEDs, we identify those stars which exhibit infrared excess, finding excess only among the brightest giants: dusty mass loss begins at a luminosity of ~ 1000 Lsun, becoming ubiquitous above 2000 Lsun. Recent claims of dust production around lower-luminosity giants cannot be reproduced, despite using the same archival Spitzer imagery.

We present the photometric catalogs for the star-forming cluster NGC 602 in the wing of the Small Magellanic Cloud covering a range of wavelengths from optical HST/ACS (F555W, F814W) and SMARTS/ANDICAM (V, I) to infrared (Spitzer/IRAC 3.6, 4.5, 5.8, and 8 micron and MIPS 24 micron). Combining this with IRSF (InfraRed Survey Facility) near-infrared photometry (J, H, Ks), we compare the young main sequence (MS) and pre-main sequence (PMS) populations prominent in the optical with the current young stellar object (YSO) populations revealed by the infrared (IR). We analyze the MS and PMS population with isochrones in color-magnitude diagrams to derive ages and masses. The optical data reveal ~565 PMS candidates, low mass Stage III YSOs. We characterize ~40 YSOs by fitting their spectral energy distributions (SEDs) to a grid of models (Robitaille et al. 2007) to derive luminosities, masses and evolutionary phase (Stage I-III). The higher resolution HST images reveal that ~70% of the YSO candidates are either multiples or protoclusters. For YSOs and PMS sources found in common, we find a consistency in the masses derived. We use the YSO mass function to derive a present-day star-formation rate of ~0.2-1.0 Msun/yr/kpc^2, similar to the rate derived from the optical star formation history suggesting a constant star formation rate for this region. We demonstrate a progression of star formation from the optical star cluster center to the edge of the star forming dust cloud. We derive lifetimes of a few 10^5 years for the YSO Stages I and II.

A group of four compact HII regions associated with the well-known 50 km/s molecular cloud is the closest site of on-going star formation to the dynamical center of the Galaxy, at a projected distance of ~6 pc. We present a study of ionized gas based on the [NeII] (12.8 micron) line, as well as multi-frequency radio continuum, HST Pa alpha and Spitzer IRAC observations of the most compact member of the HII group, Sgr A East HII D. The radio continuum image at 6cm shows that this source breaks up into two equally bright ionized features, D1 and D2. The SED of the D source is consistent with it being due to a 25\pm3 solar mass, star with a luminosity of 8\pm3x10^4 solar luminosity. The inferred mass, effective temperature of the UV source and the ionization rate are compatible with a young O9-B0 star. The ionized features D1 and D2 are considered to be ionized by UV radiation collimated by an accretion disk. We consider that the central massive star photoevaporates its circumstellar disk on a timescale of 3x10^4 years giving a mass flux ~3x10^{-5} solar mass per yr and producing the ionized material in D1 and D2 expanding in an inhomogeneous medium. The ionized gas kinematics, as traced by the [Ne II] emission, is difficult to interpret, but it could be explained by the interaction of a bipolar jet with surrounding gas along with what appears to to be a conical wall of lower velocity gas. The other HII regions, Sgr A East A-C, have morphologies and kinematics that more closely resemble cometary flows seen in other compact HII regions, where gas moves along a paraboloidal surface formed by the interaction of a stellar wind with a molecular cloud.

We present the classification of 197 point sources observed with the Infrared Spectrograph in the SAGE-Spec Legacy program on the Spitzer Space Telescope. We introduce a decision-tree method of object classification based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, which is used to classify the SAGE-Spec sample of point sources. The decision tree has a broad application to mid-infrared spectroscopic surveys, where supporting photometry and variability information are available. We use these classifications to make deductions about the stellar populations of the Large Magellanic Cloud and the success of photometric classification methods. We find 90 asymptotic giant branch (AGB) stars, 29 young stellar objects, 23 post-AGB objects, 19 red supergiants, eight stellar photospheres, seven background galaxies, seven planetary nebulae, two HII regions and 12 other objects, seven of which remain unclassified.

The LMC star, SSTISAGE1C J050756.44-703453.9, was first noticed during a survey of EROS-2 lightcurves for stars with large irregular brightness variations typical of the R Coronae Borealis (RCB) class. However, the visible spectrum showing emission lines including the Balmer and Paschen series as well as many Fe II lines is emphatically not that of an RCB star. This star has all of the characteristics of a typical UX Ori star. It has a spectral type of approximately A2 and has excited an H II region in its vicinity. However, if it is an LMC member, then it is very luminous for a Herbig Ae/Be star. It shows irregular drops in brightness of up to 2 mag, and displays the reddening and "blueing" typical of this class of stars. Its spectrum, showing a combination of emission and absorption lines, is typical of a UX Ori star that is in a decline caused by obscuration from the circumstellar dust. SSTISAGE1C J050756.44-703453.9 has a strong IR excess and significant emission is present out to 500 micron. Monte Carlo radiative transfer modeling of the SED requires that SSTISAGE1C J050756.44-703453.9 has both a dusty disk as well as a large extended diffuse envelope to fit both the mid- and far-IR dust emission. This star is a new member of the UX Ori subclass of the Herbig Ae/Be stars and only the second such star to be discovered in the LMC.

IRAS04325+2402 (herafter IRAS04325) is a complex protostellar system hosting two young stellar objects (AB and C in the following) at a separation of 1250AU. Here we present new deep Gemini imaging and spectroscopy for the system covering the wavelength regime from 1-12mu as well as Sub-Millimeter Array interferometry at 870mu, in combination with Spitzer and literature data. Object AB is a low-mass star with a disk/envelope system and an outflow cavity, which is prominently seen in infrared images. Object C, previously suspected to be a brown dwarf, is likely a very low mass star, with an effective temperature of ~3400K. It features an edge-on disk and an elongated envelope, and shows strong indications for accretion and ejection activity. Both objects are likely to drive parsec-scale molecular outflows. The two objects are embedded in an isolated, dense molecular cloud core. High extinction, lack of X-ray emission, and relatively high bolometric luminosity argue for a very young age below 1Myr. The disk/outflow systems of AB and C are misaligned by ~60deg against each other and by 80 and 40deg against the orbital plane of the binary. The system might be a good case for primordial misalignment, as opposed to misalignment caused by dynamical interactions, because the outflow direction is constant and the realignment timescale is likely larger than the system age. This favours turbulent fragmentation, rather than rotational fragmentation, as the formation scenario. We show that the spectral energy distributions and images for the two objects can be reproduced with radiative transfer models for disk/envelope systems. Our analysis provides reassurance in the established paradigm for the structure and early evolution of YSOs, but stresses the importance of developing 3D models with sophisticated dust chemistry (abridged).

We present a catalog of 5324 massive stars in the Small Magellanic Cloud (SMC), with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 3654 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer, SAGE-SMC survey database, for which we present uniform photometry from 0.3-24 um in the UBVIJHKs+IRAC+MIPS24 bands. We compare the color magnitude diagrams and color-color diagrams to those of the Large Magellanic Cloud (LMC), finding that the brightest infrared sources in the SMC are also the red supergiants, supergiant B[e] (sgB[e]) stars, luminous blue variables, and Wolf-Rayet stars, with the latter exhibiting less infrared excess, the red supergiants being less dusty and the sgB[e] stars being on average less luminous. Among the objects detected at 24 um are a few very luminous hypergiants, 4 B-type stars with peculiar, flat spectral energy distributions, and all 3 known luminous blue variables. We detect a distinct Be star sequence, displaced to the red, and suggest a novel method of confirming Be star candidates photometrically. We find a higher fraction of Oe and Be stars among O and early-B stars in the SMC, respectively, when compared to the LMC, and that the SMC Be stars occur at higher luminosities. We estimate mass-loss rates for the red supergiants, confirming the correlation with luminosity even at the metallicity of the SMC. Finally, we confirm the new class of stars displaying composite A & F type spectra, the sgB[e] nature of 2dFS1804 and find the F0 supergiant 2dFS3528 to be a candidate luminous blue variable with cold dust.

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory's PACS and SPIRE photometry cameras, we imaged a 2x8 square degree strip through the LMC, at a position angle of ~22.5 degrees as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 {\mu}m and SPIRE 250, 350 and 500 {\mu}m. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 {\mu}m emission is in excess of the models by 6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_PAH) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J=1-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65(+15,-18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287(+25,-42) that is more consistent with a GDR inferred from extinction.

We demonstrate the unique capabilities of Herschel to study very young luminous extragalactic young stellar objects (YSOs) by analyzing a central strip of the Large Magellanic Cloud obtained through the HERITAGE Science Demonstration Program. We combine PACS 100 and 160, and SPIRE 250, 350, and 500 microns photometry with 2MASS (1.25-2.17 microns) and Spitzer IRAC and MIPS (3.6-70 microns) to construct complete spectral energy distributions (SEDs) of compact sources. From these, we identify 207 candidate embedded YSOs in the observed region, ~40% never-before identified. We discuss their position in far-infrared color-magnitude space, comparing with previously studied, spectroscopically confirmed YSOs and maser emission. All have red colors indicating massive cool envelopes and great youth. We analyze four example YSOs, determining their physical properties by fitting their SEDs with radiative transfer models. Fitting full SEDs including the Herschel data requires us to increase the size and mass of envelopes included in the models. This implies higher accretion rates (greater than or equal to 0.0001 M_sun/yr), in agreement with previous outflow studies of high-mass protostars. Our results show that Herschel provides reliable longwave SEDs of large samples of high-mass YSOs; discovers the youngest YSOs whose SEDs peak in Herschel bands; and constrains the physical properties and evolutionary stages of YSOs more precisely than was previously possible.

HII regions are the birth places of stars, and as such they provide the best measure of current star formation rates (SFRs) in galaxies. The close proximity of the Magellanic Clouds allows us to probe the nature of these star forming regions at small spatial scales. We aim to determine the monochromatic IR band that most accurately traces the bolometric IR flux (TIR), which can then be used to estimate an obscured SFR. We present the spatial analysis, via aperture/annulus photometry, of 16 LMC and 16 SMC HII region complexes using the Spitzer IRAC and MIPS bands. UV rocket data and SHASSA H-alpha data are also included. We find that nearly all of the LMC and SMC HII region SEDs peak around 70um, from ~10 to ~400 pc from the central sources. As a result, the sizes of HII regions as probed by 70um is approximately equal to the sizes as probed by TIR (about 70 pc in radius); the radial profile of the 70um flux, normalized by TIR, is constant at all radii (70um ~ 0.45 TIR); the 1-sigma standard deviation of the 70um fluxes, normalized by TIR, is a lower fraction of the mean (0.05 to 0.12 out to ~220 pc) than the normalized 8, 24, and 160um normalized fluxes (0.12 to 0.52); and these results are invariant between the LMC and SMC. From these results, we argue that 70um is the most suitable IR band to use as a monochromatic obscured star formation indicator because it most accurately reproduces the TIR of HII regions in the LMC and SMC and over large spatial scales. We also explore the general trends of the 8, 24, 70, and 160um bands in the LMC and SMC HII region SEDs, radial surface brightness profiles, sizes, and normalized (by TIR) radial flux profiles. We derive an obscured SFR equation that is modified from the literature to use 70um luminosity, SFR(Mo/yr) = 9.7(0.7)x10^{-44} L(70)(ergs/s), which is applicable from 10 to 300 pc distance from the center of an HII region.

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.

The SAGE-Spec Spitzer Legacy program is a spectroscopic follow-up to the SAGE-LMC photometric survey of the Large Magellanic Cloud carried out with the Spitzer Space Telescope. We present an overview of SAGE-Spec and some of its first results. The SAGE-Spec program aims to study the life cycle of gas and dust in the Large Magellanic Cloud, and to provide information essential to the classification of the point sources observed in the earlier SAGE-LMC photometric survey. We acquired 224.6 hours of observations using the InfraRed Spectrograph and the SED mode of the Multiband Imaging Photometer for Spitzer. The SAGE-Spec data, along with archival Spitzer spectroscopy of objects in the Large Magellanic Cloud, are reduced and delivered to the community. We discuss the observing strategy, the specific data reduction pipelines applied and the dissemination of data products to the scientific community. Initial science results include the first detection of an extragalactic "21 um" feature towards an evolved star and elucidation of the nature of disks around RV Tauri stars in the Large Magellanic Cloud. Towards some young stars, ice features are observed in absorption. We also serendipitously observed a background quasar, at a redshift of z~0.14, which appears to be host-less.

Using Spitzer IRAC observations from the SAGE-SMC Legacy program and archived Spitzer IRAC data, we investigate dust production in 47 Tuc, a nearby massive Galactic globular cluster. A previous study detected infrared excess, indicative of circumstellar dust, in a large population of stars in 47 Tuc, spanning the entire Red Giant Branch (RGB). We show that those results suffered from effects caused by stellar blending and imaging artifacts and that it is likely that no stars below about 1 mag from the tip of the RGB are producing dust. The only stars that appear to harbor dust are variable stars, which are also the coolest and most luminous stars in the cluster.

We present initial results from a population synthesis model aimed at determining the star formation rate of the Milky-Way. We find that a total star formation rate of 0.68 to 1.45 Msun/yr is able to reproduce the observed number of young stellar objects in the Spitzer/IRAC GLIMPSE survey of the Galactic plane, assuming simple prescriptions for the 3D Galactic distributions of YSOs and interstellar dust, and using model SEDs to predict the brightness and color of the synthetic YSOs at different wavelengths. This is the first Galaxy-wide measurement derived from pre-main-sequence objects themselves, rather than global observables such as the total radio continuum, Halpha, or FIR flux. The value obtained is slightly lower than, but generally consistent with previously determined values. We will extend this method in the future to fit the brightness, color, and angular distribution of YSOs, and simultaneously make use of multiple surveys, to place constraints on the input assumptions, and reduce uncertainties in the star formation rate estimate. Ultimately, this will be one of the most accurate methods for determining the Galactic star formation rate, as it makes use of stars of all masses (limited only by sensitivity) rather than solely massive stars or indirect tracers of massive stars.

We combine near-infrared (2MASS) and mid-infrared (Spitzer-IRAC) photometry to characterize the IR extinction law (1.2-8 microns) over nearly 150 degrees of contiguous Milky Way midplane longitude. The relative extinctions in 5 passbands across these wavelength and longitude ranges are derived by calculating color excess ratios for G and K giant red clump stars in contiguous midplane regions and deriving the wavelength dependence of extinction in each one. Strong, monotonic variations in the extinction law shape are found as a function of angle from the Galactic center, symmetric on either side of it. These longitudinal variations persist even when dense interstellar regions, known a priori to have a shallower extinction curve, are removed. The increasingly steep extinction curves towards the outer Galaxy indicate a steady decrease in the absolute-to-selective extinction ratio (R_V) and in the mean dust grain size at greater Galactocentric angles. We note an increasing strength of the 8 micron extinction inflection at high Galactocentric angles and, using theoretical dust models, show that this behavior is consistent with the trend in R_V. Along several lines of sight where the solution is most feasible, A_lambda/A_Ks as a function of Galactic radius is estimated and shown to have a Galactic radial dependence. Our analyses suggest that the observed relationship between extinction curve shape and Galactic longitude is due to an intrinsic dependence of the extinction law on Galactocentric radius.

New 12CO J=4-3 and 13CO J=3-2 observations of the N159 region in the Large Magellanic Cloud have been made. The 12CO J=4-3 distribution is separated into three clumps. These new measurements toward the three clumps are used in coupled calculations of molecular rotational excitation and line radiation transfer, along with other transitions of the 12CO as well as the isotope transitions of 13CO. The temperatures and densities are determined to be ~70-80K and ~3x10^3 cm-3 in N159W and N159E and ~30K and ~1.6x10^3 cm-3 in N159S. These results are compared with the star formation activity. The N159E clump is associated with embedded cluster(s) as observed at 24 micron and the derived high temperature is explained as due to the heating by these sources. The N159E clump is likely responsible for a dark lane in a large HII region by the dust extinction. The N159W clump is associated with embedded clusters mainly toward the eastern edge of the clump only. These clusters show offsets of 20"-40" from the 12CO J=4-3 peak and are probably responsible for heating indicated by the derived high temperature. The N159W clump exhibits no sign of star formation toward the 12CO J=4-3 peak position and its western region. We suggest that the N159W peak represents a pre-star-cluster core of ~105M_sol which deserves further detailed studies. Note that recent star formation took place between N159W and N159E as indicated by several star clusters and HII regions, while the natal molecular gas toward the stars have already been dissipated by the ionization and stellar winds of the OB stars. The N159S clump shows little sign of star formation as is consistent with the lower temperature and somewhat lower density. The N159S clump is also a candidate for future star formation.

We investigate dust production and stellar mass loss in the Galactic globular cluster NGC 362. Due to its close proximity to the Small Magellanic Cloud (SMC), NGC 362 was imaged with the IRAC and MIPS cameras onboard the Spitzer Space Telescope as part of the Surveying the Agents of Galaxy Evolution (SAGE-SMC) Spitzer Legacy program. We detect several cluster members near the tip of the Red Giant Branch that exhibit infrared excesses indicative of circumstellar dust and find that dust is not present in measurable quantities in stars below the tip of the Red Giant Branch. We modeled the spectral energy distribution (SED) of the stars with the strongest IR excess and find a total cluster dust mass-loss rate of 3.0(+2.0/-1.2) x 10^-9 solar masses per year, corresponding to a gas mass-loss rate of 8.6(+5.6/-3.4) x 10^-6 solar masses per year, assuming [Fe/H] = -1.16. This mass loss is in addition to any dust-less mass loss that is certainly occurring within the cluster. The two most extreme stars, variables V2 and V16, contribute up to 45% of the total cluster dust-traced mass loss. The SEDs of the more moderate stars indicate the presence of silicate dust, as expected for low-mass, low-metallicity stars. Surprisingly, the SED shapes of the stars with the strongest mass-loss rates appear to require the presence of amorphous carbon dust, possibly in combination with silicate dust, despite their oxygen-rich nature. These results corroborate our previous findings in omega Centauri.

We present a catalog of 1750 massive stars in the Large Magellanic Cloud, with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 1268 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer SAGE survey database, for which we present uniform photometry from 0.3-24 microns in the UBVIJHKs+IRAC+MIPS24 bands. The resulting infrared color-magnitude diagrams illustrate that the supergiant B[e], red supergiant and luminous blue variable (LBV) stars are among the brightest infrared point sources in the Large Magellanic Cloud, due to their intrinsic brightness, and at longer wavelengths, due to dust. We detect infrared excesses due to free-free emission among ~900 OB stars, which correlate with luminosity class. We confirm the presence of dust around 10 supergiant B[e] stars, finding the shape of their spectral energy distributions (SEDs) to be very similar, in contrast to the variety of SED shapes among the spectrally variable LBVs. The similar luminosities of B[e] supergiants (log L/Lo>=4) and the rare, dusty progenitors of the new class of optical transients (e.g. SN 2008S and NGC 300 OT), plus the fact that dust is present in both types of objects, suggests a common origin for them. We find the infrared colors for Wolf-Rayet stars to be independent of spectral type and their SEDs to be flatter than what models predict. The results of this study provide the first comprehensive roadmap for interpreting luminous, massive, resolved stellar populations in nearby galaxies at infrared wavelengths.

We have used Spitzer/IRAC to conduct a photometric monitoring program of the IC1396A dark globule in order to study the mid-IR (3.6 - 8 micron) variability of the heavily embedded Young Stellar Objects (YSOs) present in that area. We obtained light curves covering a 14 day timespan with a twice daily cadence for 69 YSOs, and continuous light curves with approximately 12 second cadence over 7 hours for 38 YSOs. Typical accuracies for our relative photometry were 1-2% for the long timespan data and a few mmag, corresponding to less than 0.5%, for the 7 hour continuous "staring-mode" data. More than half of the YSOs showed detectable variability, with amplitudes from ~0.05 mag to ~0.2 mag. About thirty percent of the YSOs showed quasi-sinusoidal light curve shapes with apparent periods from 5-12 days and light curve amplitudes approximately independent of wavelength over the IRAC bandpasses. We have constructed models which simulate the time dependent spectral energy distributions of Class I and I I YSOs in order to attempt to explain these light curves. Based on these models, the apparently periodic light curves are best explained by YSO models where one or two high latitude photospheric spots heat the inner wall of the circumstellar disk, and where we view the disk at fairly large inclination angle. Disk inhomogeneities, such as increasing the height where the accretion funnel flows to the stellar hotspot, enhances the light curve modulations. The other YSOs in our sample show a range of light curve shapes, some of which are probably due to varying accretion rate or disk shadowing events. One star, IC1396A-47, shows a 3.5 hour periodic light curve; this object may be a PMS Delta Scuti star.

The central kpc of the Milky Way might be expected to differ significantly from the rest of the Galaxy with regard to gas dynamics and the formation of YSOs. We probe this possibility with mid-infrared observations obtained with IRAC and MIPS on Spitzer and with MSX. We use color-color diagrams and SED fits to explore the nature of YSO candidates (including objects with 4.5 micron excesses possibly due to molecular emission). There is an asymmetry in the distribution of the candidate YSOs, which tend to be found at negative Galactic longitudes; this behavior contrasts with that of the molecular gas, approximately 2/3 of which is at positive longitudes. The small scale height of these objects suggests that they are within the Galactic center region and are dynamically young. They lie between two layers of infrared dark clouds and may have originated from these clouds. We identify new sites for this recent star formation. The methanol masers appear to be associated with young, embedded YSOs characterized by 4.5 micron excesses. We use the SEDs of these sources to estimate their physical characteristics. Within the central 400x50 pc (|l|<1.3\degr and |b|<10') the star formation rate based on the identification of Stage I evolutionary phase of YSO candidates is about 0.14 solar mass/yr. We suggest that a recent burst of star formation took place within the last 10^5 years. This suggestion is also consistent with estimates of star formation rates within the last ~10^7 years showing a peak around 10^5 years ago. Lastly, we find that the Schmidt-Kennicutt Law applies well in the central 400 pc of the Galaxy. This implies that star formation does not appear to be dramatically affected by the extreme physical conditions in the Galactic center region.

We present high spatial resolution mid-IR images of the ring of UCHII regions in W49A obtained at Gemini North, allowing us to identify the driving source of its powerful H2O maser outflow. These data also confirm our previous report that several radio sources in the ring are undetected in the mid-IR because they are embedded deep inside the cloud core. We locate the source of the water maser outflow at the position of the compact mid-IR peak of source G (source G:IRS1). This IR source is not coincident with any identified compact radio continuum source, but is coincident with a hot molecular core, so we propose that G:IRS1 is a hot core driving an outflow analogous to the wide-angle bipolar outflow in OMC-1. G:IRS1 is at the origin of a larger bipolar cavity and CO outflow. The water maser outflow is orthogonal to the bipolar CO cavity, so the masers probably reside near its waist in the cavity walls. Models of the IR emission require a massive protostar of 45Msun, 3e5Lsun, and an effective envelope accretion rate of 1e-3Msun/yr. Feedback from the central star could potentially drive the H2O maser outflow, but it has insufficient radiative momentum to have driven the large-scale CO outflow, requiring that this massive star had an active accretion disk over the past 10^4 yr. Combined with the spatialy resolved morphology in IR images, G:IRS1 in W49 provides compelling evidence for a massive protostar that formed by accreting from a disk, accompanied by a bipolar outflow.

Massive young stellar objects (YSOs), like low-mass YSOs, appear to be surrounded by optically thick envelopes and/or disks and have regions, often bipolar, that are seen in polarized scattered light at near-infrared wavelengths. We are using the 0.2'' spatial resolution of NICMOS on Hubble Space Telescope to examine the structure of the disks and outflow regions of massive YSOs in star-forming regions within a few kpc of the Sun. Here we report on 2 micron polarimetry of NGC 6334 V and S255 IRS1. NGC 6334 V consists of a double-lobed bright reflection nebula seen against a dark region, probably an optically thick molecular cloud. Our polarization measurements show that the illuminating star lies ~ 2'' south of the line connecting the two lobes; we do not detect this star at 2 micron, but there are a small radio source and a mid-infrared source at this location. S255 IRS1 consists of two YSOs (NIRS1 and NIRS3) with overlapping scattered light lobes and luminosities corresponding to early B stars. Included in IRS1 is a cluster of stars from whose polarization we determine the local magnetic field direction. Neither YSO has its scattered light lobes aligned with this magnetic field. The line connecting the scattered light lobes of NIRS1 is twisted symmetrically around the star; the best explanation is that the star is part of a close binary and the outflow axis of NIRS1 is precessing as a result of non-coplanar disk and orbit. The star NIRS3 is also offset from the line connecting its two scattered light lobes. We suggest that all three YSOs show evidence of episodic ejection of material as they accrete from dense, optically thick envelopes.

We present empirical relations describing excess emission from evolved stars in the Large Magellanic Cloud (LMC) using data from the SAGE (Surveying the Agents of a Galaxy's Evolution) survey which includes the IRAC 3.6, 4.5, 5.8 and 8.0 \mu m and MIPS 24, 70 and 160 \mu m bands. We combine the SAGE data with the Two Micron All Sky Survey (2MASS; J, H and Ks) and the optical Magellanic Cloud Photometric Survey (MCPS; U, B, V and I) point source catalogs to create complete spectral energy distributions (SEDs) of the asymptotic giant branch (AGB) star candidates in the LMC. AGB star outflows are among the main producers of dust in a galaxy, and this mass loss results in an excess in the fluxes observed in the 8 and 24 \mic m bands. We identify oxygen-rich, carbon-rich and extreme AGB star populations in our sample based on their 2MASS and IRAC colors. We calculate excesses fluxes in the mid-IR bands by comparison of the SEDs with model photospheres. We find about 16,000 O-rich, 6300 C-rich and 1000 extreme sources with reliable 8 \mu m excesses, and about 4500 O-rich, 5300 C-rich and 960 extreme sources with reliable 24 \mic m excesses. The excesses are in the range 0.1 mJy--5 Jy. The 8 and 24 \mic m excesses for all three types of AGB candidates show a general increasing trend with luminosity. The dust color temperature derived from the ratio of the 8 and 24 \mic m excesses decreases with an increase in excess, while the 24 \mic m optical depth increases with excess. The extreme AGB candidates are the major contributors to the mass loss, and we estimate the total AGB mass-loss return to the LMC to be (5.9--13)\times 10^{-3} Msun yr^{-1}.

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.

(abridged) We study photoelectric heating throughout the Large Magellanic Cloud. We quantify the importance of the [CII] cooling line and the photoelectric heating process of various environments in the LMC and investigate which parameters control the extent of photoelectric heating. We use the BICE [CII] map and the Spitzer/SAGE infrared maps. We examine the spatial variations in the efficiency of photoelectric heating: photoelectric heating rate over power absorbed by grains. We correlate the photoelectric heating efficiency and the emission from various dust constituents and study the variations as a function of H\alpha emission, dust temperatures, and the total infrared luminosity. From this we estimate radiation field, gas temperature, and electron density. We find systematic variations in photoelectric efficiency. The highest efficiencies are found in the diffuse medium, while the lowest coincide with bright star-forming regions (~1.4 times lower). The [CII] line emission constitutes 1.32% of the far infrared luminosity across the whole of the LMC. We find correlations between the [CII] emission and ratios of the mid infrared and far infrared bands, which comprise various dust constituents. The correlations are interpreted in light of the spatial variations of the dust abundance and by the local environmental conditions that affect the dust emission properties. As a function of the total infrared surface brightness, S_{TIR}, the [CII] surface brightness can be described as: S_{[CII]}=1.25 S_{TIR}^{0.69} [10^{-3} erg s^{-1} cm^{-2} sr^{-1}]. The [CII] emission is well-correlation with the 8 micrometer emission, suggesting that the polycyclic aromatic hydrocarbons play a dominant role in the photoelectric heating process.

The Tail region of the Small Magellanic Cloud (SMC) was imaged using the MIPS instrument on the Spitzer Space Telescope as part of the SAGE-SMC Spitzer Legacy. Diffuse infrared emission from dust was detected in all the MIPS bands. The Tail gas-to-dust ratio was measured to be 1200 +/- 350 using the MIPS observations combined with existing IRAS and HI observations. This gas-to-dust ratio is higher than the expected 500-800 from the known Tail metallicity indicating possible destruction of dust grains. Two cluster regions in the Tail were resolved into multiple sources in the MIPS observations and local gas-to-dust ratios were measured to be ~440 and ~250 suggests dust formation and/or significant amounts of ionized gas in these regions. These results support the interpretation that the SMC Tail is a tidal tail recently stripped from the SMC that includes gas, dust, and young stars.

We present initial results and source lists of variable sources in the Large Magellanic Cloud (LMC) for which we detect thermal infrared variability from the SAGE (Surveying the Agents of a Galaxy's Evolution) survey, which had 2 epochs of photometry separated by three months. The SAGE survey mapped a 7 degree by 7 degree region of the LMC using the IRAC and the MIPS instruments on board Spitzer. Variable sources are identified using a combination of the IRAC 3.6, 4.5, 5.8, 8.0 \micron bands and the MIPS 24 \micron bands. An error-weighted flux difference between the two epochs is used to assess the variability. Of the ~ 3 million sources detected at both epochs we find ~ 2,000 variable sources for which we provide electronic catalogs. Most of the variable sources can be classified as asymptotic giant branch (AGB) stars. A large fraction (> 66%) of the extreme AGB stars are variable and only smaller fractions of carbon-rich (6.1%) and oxygen-rich (2.0%) stars are detected as variable. We also detect a population of variable young stellar object candidates.

We review the history and structure of star formation in the Taurus-Auriga dark clouds. Our discussion includes a summary of the macroscopic cloud properties, the population of single and binary pre-main sequence stars, the properties of jets and outflows, and detailed summaries of selected individual objects. We include comprehensive tables of dark clouds, young stars, and jets in the clouds.

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.

The spectral energy distributions (SEDs), spanning the mid-infrared to millimeter wavelengths, of a sample of 13 high-mass protostellar objects (HMPOs) were studied using a large archive of 2-D axisymmetric radiative transfer models. Measurements from the Spitzer GLIMPSE and MIPSGAL surveys and the MSX survey were used in addition to our own surveys at millimeter and submillimeter wavelengths to construct the SEDs, which were then fit to the archive of models. These models assumed that stars of all masses form via accretion and allowed us to make estimates for the masses, luminosities and envelope accretion rates for the HMPOs. The models fit the observed SEDs well. The implied envelope accretion rates are high, $\approx 10^{-2.5} \msun/yr$, consistent with the accretion-based scenario of massive star formation. With the fitted accretion rates and with mass estimates of up to $\sim 20 \msun$ for these objects, it appears plausible that stars with stellar masses $M_{\ast} > 20 \msun$ can form via accretion.

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 have modeled optical to far infrared images, photometry, and spectroscopy of the object known as Gomez's Hamburger. We reproduce the images and spectrum with an edge-on disk of mass 0.3M_sun and radius 1600AU, surrounding an A0 III star at a distance of 280pc. Our mass estimate is in excellent agreement with recent CO observations. However, our distance determination is more than an order of magnitude smaller than previous analyses which inaccurately interpreted the optical spectrum. To accurately model the infrared spectrum we have extended our Monte Carlo radiation transfer codes to include emission from polycyclic aromatic hydrocarbon (PAH) molecules and very small grains (VSG). We do this using pre-computed PAH/VSG emissivity files for a wide range of values of the mean intensity of the exciting radiation field. When Monte Carlo energy packets are absorbed by PAHs/VSGs we reprocess them to other wavelengths by sampling from the emissivity files, thus simulating the absorption and re-emission process without reproducing lengthy computations of statistical equilibrium, excitation and de-excitation in the complex many level molecules. Using emissivity lookup tables in our Monte Carlo codes gives the flexibility to use the latest grain physics calculations of PAH/VSG emissivity and opacity that are being continually updated in the light of higher resolution infrared spectra. We find our approach gives a good representation of the observed PAH spectrum from the disk of Gomez's Hamburger. Our models also indicate the PAHs/VSGs in the disk have a larger scaleheight than larger radiative equilibrium grains, providing evidence for dust coagulation and settling to the midplane.

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 present three-dimensional, self-consistent radiative transfer solutions with a new Monte Carlo radiative equilibrium code. The code, RADISHE ($\bf{RAD}$iative transfer $\bf{I}$n $\bf{S}$moothed particle $\bf{H}$ydrodynamics and $\bf{E}$ulerian codes), can be applied to calculate the emergent spectral energy distributions (SEDs) and broadband images from optical to millimeter wavelengths of arbitrary density geometries with distributed sources of radiation. One of the primary uses of this code has been to interface with hydrodynamical codes to calculate emergent SEDs along a simulation time sequence. The primary methodological focus of this paper is on the radiative equilibrium temperature calculation. We find that an iterative calculation of the temperature, which takes as the Monte Carlo estimator for the mean free intensity the sum of photon flight paths, is significantly faster than relaxation temperature calculation methods, particularly when large numbers of grid cells are required, i.e., in modeling three-dimensional geometries such as the dust envelopes of turbulent massive protostellar cores or infrared bright galaxies. We present simulated color-color plots for infrared bright galaxies at a range of redshifts, and unfold these plots as color vs the fractional AGN luminosity, to demonstrate that $\it{Herschel}$ will be able to effectively discriminate between submillimeter galaxies where the energy source is dominated by AGN and those where star formation dominates. [abridged]

We present IRAC and MIPS images and photometry of a sample of previously known planetary nebulae (PNe) from the SAGE survey of the Large Magellanic Cloud (LMC) performed with the Spitzer Space Telescope. Of the 233 known PNe in the survey field, 185 objects were detected in at least two of the IRAC bands, and 161 detected in the MIPS 24 micron images. Color-color and color-magnitude diagrams are presented using several combinations of IRAC, MIPS, and 2MASS magnitudes. The location of an individual PN in the color-color diagrams is seen to depend on the relative contributions of the spectral components which include molecular hydrogen, polycyclic aromatic hydrocarbons (PAHs), infrared forbidden line emission from the ionized gas, warm dust continuum, and emission directly from the central star. The sample of LMC PNe is compared to a number of Galactic PNe and found to not significantly differ in their position in color-color space. We also explore the potential value of IR PNe luminosity functions (LFs) in the LMC. IRAC LFs appear to follow the same functional form as the well-established [O III] LFs although there are several PNe with observed IR magnitudes brighter than the cut-offs in these LFs.

We present a list of 552 sources with suspected variability, based on a comparison of mid-infrared photometry from the GLIMPSE I and MSX surveys, which were carried out nearly a decade apart. We were careful to address issues such as the difference in resolution and sensitivity between the two surveys, as well as the differences in the spectral responses of the instruments. We selected only sources where the IRAC 8.0 and MSX 8.28 micron fluxes differ by more than a factor of two, in order to minimize contamination from sources where the difference in fluxes at 8 micron is due to a strong 10 micron silicate feature. We present a subset of 40 sources for which additional evidence suggests variability, using 2MASS and MIPSGAL data. Based on a comparison with the variability flags in the IRAS and MSX Point-Source Catalogs we estimate that at least a quarter of the 552 sources, and at least half of the 40 sources are truly variable. In addition, we tentatively confirm the variability of one source using multi-epoch IRAS LRS spectra. We suggest that most of the sources in our list are likely to be Asymptotic Giant Branch stars.

We use Spitzer Space Telescope observations from the Spitzer Survey of the Small Magellanic Cloud (S3MC) to study the young stellar content of N66, the largest and brightest HII region in the SMC. In addition to large numbers of normal stars, we detect a significant population of bright, red infrared sources that we identify as likely to be young stellar objects (YSOs). We use spectral energy distribution (SED) fits to classify objects as ordinary (main sequence or red giant) stars, asymptotic giant branch stars, background galaxies, and YSOs. This represents the first large-scale attempt at blind source classification based on Spitzer SEDs in another galaxy. We firmly identify at least 61 YSOs, with another 50 probable YSOs; only one embedded protostar in the SMC was reported in the literature prior to the S3MC. We present color selection criteria that can be used to identify a relatively clean sample of YSOs with IRAC photometry. Our fitted SEDs indicate that the infrared-bright YSOs in N66 have stellar masses ranging from 2 Msun to 17 Msun, and that approximately half of the objects are Stage II protostars, with the remaining YSOs roughly evenly divided between Stage I and Stage III sources. We find evidence for primordial mass segregation in the HII region, with the most massive YSOs being preferentially closer to the center than lower-mass objects. Despite the low metallicity and dust content of the SMC, the observable properties of the YSOs appear consistent with those in the Milky Way. Although the YSOs are heavily concentrated within the optically bright central region of N66, there is ongoing star formation throughout the complex and we place a lower limit on the star formation rate of 3.2 x 10^-3 Msun/yr over the last ~1 Myr.

We present new Spitzer photometry of the Eagle Nebula (M16, containing the optical cluster NGC 6611) combined with near-infrared photometry from 2MASS. We use dust radiative transfer models, mid-infrared and near-infrared color-color analysis, and mid-infrared spectral indices to analyze point source spectral energy distributions, select candidate young stellar objects (YSOs), and constrain their mass and evolutionary state. Comparison of the different protostellar selection methods shows that mid-infrared methods are consistent, but as has been known for some time, near-infrared-only analysis misses some young objects. We reveal more than 400 protostellar candidates, including one massive young stellar object (YSO) that has not been previously highlighted. The YSO distribution supports a picture of distributed low-level star formation, with no strong evidence of triggered star formation in the ``pillars''. We confirm the youth of NGC 6611 by a large fraction of infrared-excess sources, and reveal a younger cluster of YSOs in the nearby molecular cloud. Analysis of the YSO clustering properties shows a possible imprint of the molecular cloud's Jeans length. Multiwavelength mid-IR imaging thus allows us to analyze the protostellar population, to measure the dust temperature and column density, and to relate these in a consistent picture of star formation in M16.

We compare H-alpha, radio continuum, and Spitzer Space Telescope (SST) images of 58 planetary nebulae (PNe) recently discovered by the Macquarie-AAO-Strasbo- urg H-alpha PN Project (MASH) of the SuperCOSMOS H-alpha Survey. Using InfraRed Array Camera (IRAC) data we define the IR colors of PNe and demonstrate good isolation between these colors and those of many other types of astronomical object. The only substantive contamination of PNe in the color-color plane we illustrate is due to YSOs. However, this ambiguity is readily resolved by the unique optical characteristics of PNe and their environs. We also examine the relationships between optical and MIR morphologies from 3.6 to 8.0um and explore the ratio of mid-infrared (MIR) to radio nebular fluxes, which is a valuable discriminant between thermal and nonthermal emission. MASH emphasizes late evolutionary stages of PNe compared with previous catalogs, enabling study of the changes in MIR and radio flux that attend the aging process. Spatially integrated MIR energy distributions were constructed for all MASH PNe observed by the GLIMPSE Legacy Project, using the H-alpha morphologies to establish the dimensions for the calculations of the Midcourse Space Experiment (MSX), IRAC, and radio continuum (from the Molonglo Observatory Synthesis Telescope and the Very Large Array) flux densities. The ratio of IRAC 8.0-um to MSX 8.3-um flux densities provides a measure of the absolute diffuse calibration of IRAC at 8.0 um. We independently confirm the aperture correction factor to be applied to IRAC at 8.0um to align it with the diffuse calibration of MSX. The result agrees with the recommendations of the Spitzer Science Center and with results from a parallel study of HII regions. These PNe probe the diffuse calibration of IRAC on a spatial scale of 9-77 arcsec.

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.

OH(1720 MHz) and methanol masers are now recognized to be excellent probes of the interactions of supernova remnants with molecular clouds and tracers of massive star formation, respectively. To better understand the nature of star formation activity in the central region of the Galaxy, we have used these two classes of masers combined with the IRAC and MIPS data to study prominent sites of ongoing star formation in the nuclear disk. The nuclear disk is characterized by massive GMCs with elevated gas temperatures, compared to their dust temperatures. We note an association between methanol masers and a class of mid-infrared ``green sources''. These highly embedded YSOs show enhanced 4.5micron emission due to excited molecular lines. The distribution of methanol masers and supernova remnants suggest a low efficiency of star formation (with the exception of Sgr B2), which we believe is due to an enhanced flux of cosmic ray electrons impacting molecular clouds in the nuclear disk. We also highlight the importance of cosmic rays in their ability to heat molecular clouds, and thus increase the gas temperature.

We present a method to analyze the spectral energy distributions (SEDs) of young stellar objects (YSOs). Our approach is to fit data with pre-computed 2-D radiation transfer models spanning a large region of parameter space. This allows us to determine not only a single set of physical parameter values but the entire range of values consistent with the multi-wavelength observations of a given source. In this way we hope to avoid any over-interpretation when modeling a set of data. We have constructed spectral energy distributions from optical to sub-mm wavelengths, including new Spitzer IRAC and MIPS photometry, for 30 young and spatially resolved sources in the Taurus-Auriga star-forming region. We demonstrate fitting model SEDs to these sources, and find that we correctly identify the evolutionary stage and physical parameters found from previous independent studies, such as disk mass, disk accretion rate, and stellar temperature. We also explore how fluxes at various wavelengths help to constrain physical parameters, and show examples of degeneracies that can occur when fitting SEDs. A web-based version of this tool is available to the community at http://www.astro.wisc.edu/protostars .

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 report on the discovery of two spiral galaxies located behind the southern Milky Way, within the least explored region of the Great Attractor. They lie at 317, -0.5 deg galactic, where obscuration from Milky Way stars and dust exceeds 13 to 15 mag of visual extinction. The galaxies were the most prominent of a set identified using mid-infrared images of the low-latitude (|b| < 1 deg) Spitzer Legacy program GLIMPSE. Follow-up HI radio observations reveal that both galaxies have redshifts that place them squarely in the Norma Wall of galaxies, which appears to extend diagonally across the Galactic Plane from Norma in the south to Centaurus & Vela in the north. We report on the near-infrared, mid-infrared and radio properties of these newly discovered galaxies, and discuss their context in the larger view of the Great Attractor. The work presented here demonstrates that mid-infrared surveys open up a new window to study galaxies in the Zone of Avoidance.

We investigate the diffuse absolute calibration of the InfraRed Array Camera on the Spitzer Space Telescope at 8.0microns using a sample of 43 HII regions with a wide range of morphologies near GLON=312deg. For each region we carefully measure sky-subtracted,point-source- subtracted, areally-integrated IRAC 8.0-micron fluxes and compare these with Midcourse Space eXperiment (MSX) 8.3-micron images at two different spatial resolutions, and with radio continuum maps. We determine an accurate median ratio of IRAC 8.0-micron/MSX\8.3-micron fluxes, of 1.55+/-0.15. From robust spectral energy distributions of these regions we conclude that the present 8.0-micron diffuse calibration of the SST is 36% too high compared with the MSX validated calibration, perhaps due to scattered light inside the camera. This is an independent confirmation of the result derived for the diffuse calibration of IRAC by the Spitzer Science Center (SSC). From regression analyses we find that 843-MHz radio fluxes of HII regions and mid-infrared (MIR) fluxes are linearly related for MSX at 8.3-microns and Spitzer at 8.0 microns, confirming the earlier MSX result by Cohen & Green. The median ratio of MIR/843-MHz diffuse continuum fluxes is 600 times smaller in nonthermal than thermal regions, making it a sharp discriminant. The ratios are largely independent of morphology up to a size of ~24 arcsec. We provide homogeneous radio and MIR morphologies for all sources. MIR morphology is not uniquely related to radio structure. Compact regions may have MIR filaments and/or diffuse haloes, perhaps infrared counter- parts to weakly ionized radio haloes found around compact HII regions. We offer two IRAC colour-colour plots as quantitative diagnostics of diffuse HII regions.

We present the results of a search for new members of the Taurus star-forming region using the Infrared Array Camera (IRAC) aboard the Spitzer Space Telescope}. With IRAC images of 29.7 deg^2 of Taurus at 3.6, 4.5, 5.8, and 8.0 um, we have identified sources with red mid-infrared colors indicative of disk-bearing objects and have obtained optical and infrared spectra of 23 of these candidate members. Through this work, we have discovered 13 new members of Taurus, two of which have spectral types later than M6 and thus are likely to be brown dwarfs according to the theoretical evolutionary models of Chabrier and Baraffe. This survey indicates that the previous census of Taurus has a completeness of ~80% for members with disks. The new members that we have found do not significantly modify the previously measured distributions of Taurus members as a function of position, mass, and extinction. For instance, we find no evidence for a population of highly reddened brown dwarfs (A_K~2) that has been missed by previous optical and near-infrared surveys, which suggests that brown dwarf disks are not significantly more flared than disks around stars. In addition to the new members, we also present IRAC photometry for the 149 previously known members that appear within this survey, which includes 27 objects later than M6.

We present a grid of radiation transfer models of axisymmetric young stellar objects (YSOs), covering a wide range of stellar masses (from 0.1Msun to 50Msun) and evolutionary stages (from the early envelope infall stage to the late disk-only stage). The grid consists of 20,000 YSO models, with spectral energy distributions (SEDs) and polarization spectra computed at ten viewing angles for each model, resulting in a total of 200,000 SEDs. [...]. These models are publicly available on a dedicated WWW server: http://www.astro.wisc.edu/protostars/ . In this paper we summarize the main features of our models, as well as the range of parameters explored. [...]. We examine the dependence of the spectral indices of the model SEDs on envelope accretion rate and disk mass. In addition, we show variations of spectral indices with stellar temperature, disk inner radius, and disk flaring power for a subset of disk-only models. We also examine how changing the wavelength range of data used to calculate spectral indices affects their values. We show sample color-color plots of the entire grid as well as simulated clusters at various distances with typical {\it Spitzer Space Telescope} sensitivities. We find that young embedded sources generally occupy a large region of color-color space due to inclination and stellar temperature effects. Disk sources occupy a smaller region of color-color space, but overlap substantially with the region occupied by embedded sources, especially in the near- and mid-IR. We identify regions in color-color space where our models indicate that only sources at a given evolutionary stage should lie. [...].

Color-magnitude diagrams (CMDs) are presented for the Spitzer SAGE (Surveying the Agents of a Galaxy's Evolution) survey of the Large Magellanic Cloud (LMC). IRAC and MIPS 24 um epoch one data are presented. These data represent the deepest, widest mid-infrared CMDs of their kind ever produced in the LMC. Combined with the 2MASS survey, the diagrams are used to delineate the evolved stellar populations in the Large Magellanic Cloud as well as Galactic foreground and extragalactic background populations. Some 32000 evolved stars brighter than the tip of the red giant branch are identified. Of these, approximately 17500 are classified as oxygen-rich, 7000 carbon-rich, and another 1200 as ``extreme'' asymptotic giant branch (AGB) stars. Brighter members of the latter group have been called ``obscured'' AGB stars in the literature owing to their dusty circumstellar envelopes. A large number (1200) of luminous oxygen--rich AGB stars/M supergiants are also identified. Finally, there is strong evidence from the 24 um MIPS channel that previously unexplored, lower luminosity oxygen-rich AGB stars contribute significantly to the mass loss budget of the LMC (1200 such sources are identified).

Diffuse emission in the mid-infrared shows a wealth of structure, that lends itself to high-resolution structure analysis of the interstellar gas. A large part of the emission comes from polycyclic aromatic hydrocarbons, excited by nearby ultra-violet sources. Can the observed diffuse emission structure be interpreted as column density structure? We discuss this question with the help of a set of model molecular clouds bathed in the radiation field of a nearby O-star. The correlation strength between column density and ``observed'' flux density strongly depends on the absolute volume density range in the region. Shadowing and irradiation effects may completely alter the appearance of an object. Irradiation introduces additional small-scale structure and it can generate structures resembling shells around HII-regions in objects that do not possess any shell-like structures whatsoever. Nevertheless, structural information about the underlying interstellar medium can be retrieved. In the more diffuse regime ($n({HI})\lesssim 100$cm$^{-3}$), flux density maps may be used to trace the 3D density structure of the cloud via density gradients. Thus, while caution definitely is in order, mid-infrared surveys such as GLIMPSE will provide quantitative insight into the turbulent structure of the interstellar medium.

We are performing a uniform and unbiased, ~7x7 degrees imaging survey of the Large Magellanic Cloud (LMC), using the IRAC and MIPS instruments on board the Spitzer Space Telescope in order to survey the agents of a galaxy's evolution (SAGE), the interstellar medium (ISM) and stars in the LMC. The detection of diffuse ISM with column densities >1.2x10^21 H cm^-2 permits detailed studies of dust processes in the ISM. SAGE's point source sensitivity enables a complete census of newly formed stars with masses >3 solar masses that will determine the current star formation rate in the LMC. SAGE's detection of evolved stars with mass loss rates >1x10^-8 solar masses per year will quantify the rate at which evolved stars inject mass into the ISM of the LMC. The observing strategy includes two epochs in 2005, separated by three months, that both mitigate instrumental artifacts and constrain source variability. The SAGE data are non-proprietary. The data processing includes IRAC and MIPS pipelines and a database for mining the point source catalogs, which will be released to the community in support of Spitzer proposal cycles 4 and 5. We present initial results on the epoch 1 data with a special focus on the N79 and N83 region. The SAGE epoch 1 point source catalog has ~4 million sources. The point source counts are highest for the IRAC 3.6 microns band and decrease dramatically towards longer wavelengths consistent with the fact that stars dominate the point source catalogs and that the dusty objects, e.g. young stellar objects and dusty evolved stars that detected at the longer wavelengths, are rare in comparison. We outline a strategy for identifying foreground MW stars, that may comprise as much as 18% of the source list, and background galaxies, that may comprise ~12% of the source list.

Using Spitzer Space Telescope and Hubble Space Telescope observations of the superbubble N51D, we have identified three young stellar objects (YSOs) in dust globules, and made the first detection of a Herbig-Haro object outside the Galaxy. The spectral energy distributions of these YSOs suggest young massive stars with disk, envelope, and outflow cavities. The interstellar conditions are used to assess whether the star formation was spontaneous or induced by external pressure.

Using Infrared Array Camera (IRAC) images at 3.6, 4.5, 5.8, and 8 microns from the GLIMPSE Legacy science program on the Spitzer Space Telescope, we searched for infrared counterparts to the 95 known supernova remnants that are located within galactic longitudes 65>|l|>10 degrees and latitudes |b|<1 degree. Eighteen infrared counterparts were detected. Many other supernova remnants could have significant infrared emission but are in portions of the Milky Way too confused to allow separation from bright HII regions and pervasive mid-infrared emission from atomic and molecular clouds along the line of sight. Infrared emission from supernova remnants originates from synchrotron emission, shock-heated dust, atomic fine-structure lines, and molecular lines. The detected remnants are G11.2-0.3, Kes 69, G22.7-0.2, 3C 391, W 44, 3C 396, 3C 397, W 49B, G54.4-0.3, Kes 17, Kes 20A, RCW 103, G344.7-0.1, G346.6-0.2, CTB 37A, G348.5-0.0, and G349.7+0.2. The infrared colors suggest emission from molecular lines (9 remnants), fine-structure lines (3), and PAH (4), or a combination; some remnants feature multiple colors in different regions. None of the remnants are dominated by synchrotron radiation at mid-infrared wavelengths. The IRAC-detected sample emphasizes remnants interacting with relatively dense gas, for which most of the shock cooling occurs through molecular or ionic lines in the mid-infrared.

We use 3-D radiative transfer models to show the effects of clumpy circumstellar material on the observed infrared colors of high mass stars embedded in molecular clouds. We highlight differences between 3-D clumpy and 1-D smooth models which can affect the interpretation of data. We discuss several important properties of the emergent spectral energy distribution (SED): More near-infrared light (scattered and direct from the central source) can escape than in smooth 1-D models. The near- and mid-infrared SED of the same object can vary significantly with viewing angle, depending on the clump geometry along the sightline. Even the wavelength-integrated flux can vary with angle by more than a factor of two. Objects with the same average circumstellar dust distribution can have very different near-and mid-IR SEDs depending on the clump geometry and the proximity of the most massive clump to the central source. Although clumpiness can cause similar objects to have very different SEDs, there are some observable trends. Near- and mid-infrared colors are sensitive to the weighted average distance of clumps from the central source and to the magnitude of clumpy density variations (smooth-to-clumpy ratio). Far-infrared emission remains a robust measure of the total dust mass. We present simulated SEDs, colors, and images for 2MASS and Spitzer filters. We compare to observations of some UCHII regions and find that 3-D clumpy models fit better than smooth models. In particular, clumpy models with fractal dimensions in the range 2.3-2.8, smooth to clumpy ratios of <50%, and density distributions with shallow average radial density profiles fit the SEDs best.

We study the effects of an annular gap induced by an embedded proto-planet on disk scattered light images and the infrared spectral energy distribution. We find that the outer edge of a gap is brighter in the scattered light images than a similar location in a gap-free disk. The stellar radiation that would have been scattered by material within in the gap is instead scattered by the disk wall at the outer edge of the gap, producing a bright ring surrounding the dark gap in the images. Given sufficient resolution, such gaps can be detected by the presence of this bright ring in scattered light images. A gap in a disk also changes the shape of the SED. Radiation that would have been absorbed by material in the gap is instead reprocessed by the outer gap wall. This leads to a decrease in the SED at wavelengths corresponding to the temperature at the radius of the missing gap material and a corresponding flux increase at longer wavelengths corresponding to the temperature of the outer wall. We note that, unlike an inner hole in the disk, the presence of an annular gap does not change the bolometric IR flux; it simply redistributes the radiation, previously produced by material within the gap, to longer wavelengths. This implies that the changes in the SED generally will be smaller for gaps than holes. Although it will be difficult on the basis of the SED alone to distinguish between the presence of a gap and other physical effects, the level of changes can be sufficiently large to be measurable with current instruments (e.g., Spitzer).

The GLIMPSE (Galactic Legacy Mid-Plane Survey Extraordinaire) Point Source Catalog of ~ 30 million mid-infrared sources towards the inner Galaxy, 10 < |l| < 65 degrees and |b| < 1 degree, was used to determine the distribution of stars in Galactic longitude, latitude, and apparent magnitude. The counts versus longitude can be approximated by the modified Bessel function N=N_0*(l/l_0)*K_1(l/l_0), where l_0 is insensitive to limiting magnitude, band choice, and side of Galactic center: l_0= 17-30 degrees with a best fit value in the the 4.5 micron band of l_0=24 +/- 4 degrees. Modeling the source distribution as an exponential disk yields a radial scale length of H= 3.9 +/- 0.6 kpc. There is a pronounced north-south asymmetry in source counts for |l| < 30 degrees, with ~ 25% more stars in the north. For l=10-30 degrees, there is a strong enhancement of stars of m= 11.5-13.5 mag. A linear bar passing through the Galactic center with half-length R_bar=4.4 +/- 0.5 kpc, tilted by phi=44 +/- 10 degrees to the Sun-Galactic Center line, provides the simplest interpretation of this data. We examine the possibility that enhanced source counts at l=26-28 degrees, 31.5-34 degrees, and 306-309 degrees are related to Galactic spiral structure. Total source counts are depressed in regions where the counts of red objects (m_K-m_[8.0] >3) peak. In these areas, the counts are reduced by extinction due to molecular gas and/or high diffuse backgrounds associated with star formation.

2-D and 3-D radiation transfer models of forming stars generally produce bluer 1-10 micron colors than 1-D models of the same evolutionary state and envelope mass. Therefore, 1-D models of the shortwave radiation will generally estimate a lower envelope mass and later evolutionary state than multidimensional models. 1-D models are probably reasonable for very young sources, or longwave analysis (wavelengths > 100 microns). In our 3-D models of high-mass stars in clumpy molecular clouds, we find no correlation between the depth of the 10 micron silicate feature and the longwave (> 100 micron) SED (which sets the envelope mass), even when the average optical extinction of the envelope is >100 magnitudes. This is in agreement with the observations of Faison et al. (1998) of several UltraCompact HII (UCHII) regions, suggesting that many of these sources are more evolved than embedded protostars. We have calculated a large grid of 2-D models and find substantial overlap between different evolutionary states in the mid-IR color-color diagrams. We have developed a model fitter to work in conjunction with the grid to analyze large datasets. This grid and fitter will be expanded and tested in 2005 and released to the public in 2006.

Spitzer IRAC 3.6-8 micron photometry obtained as part of the GLIMPSE survey has revealed mid-infrared excesses for 33 field stars with known spectral types in a 1.2 sq. degree field centered on the southern Galactic HII region RCW49. These stars comprise a subset of 184 stars with known spectral classification, most of which were pre-selected to have unusually red IR colors. We propose that the mid-IR excesses are caused by circumstellar dust disks that are either very late remnants of stellar formation or debris disks generated by planet formation. Of these 33 stars, 29 appear to be main-sequence stars based on optical spectral classifications. Five of the 29 main-sequence stars are O or B stars with excesses that can be plausibly explained by thermal bremsstrahlung emission, and four are post main-sequence stars. The lone O star is an O4V((f)) at a spectrophotometric distance of 3233+ 540- 535 pc and may be the earliest member of the Westerlund 2 cluster. Of the remaining 24 main-sequence stars, 18 have SEDs that are consistent with hot dusty debris disks, a possible signature of planet formation. Modeling the excesses as blackbodies demonstrates that the blackbody components have fractional bolometric disk-to-star luminosity ratios, LIR/L*, ranging from 10^-3 to 10^-2 with temperatures ranging from 220 to 820 K. The inferred temperatures are more consistent with asteroid belts rather than the cooler temperatures expected for Kuiper belts. Mid-IR excesses are found in all spectral types from late B to early K.

We present a new planetary nebula, first identified in images from the Australia Telescope Compact Array, although not recognized at that time. Recent observations with the Spitzer Space Telescope during the GLIMPSE Legacy program have rediscovered the object. The high-resolution radio and infrared images enable the identification of the central star or its wind, the recognition of the radio emission as thermal, and the probable presence of polycylic aromatic hydrocarbons in and around the source. These lead to the conclusion that G313.3+00.3 is a planetary nebula. This object is of particular interest because it was discovered solely through radio and mid-infrared imaging, without any optical (or near-infrared) confirmation, and acts as a proof of concept for the discovery of many more highly extinguished planetary nebulae. G313.3+00.3 is well-resolved by both the instruments with which it was identified, and suffers extreme reddening due to its location in the Scutum-Crux spiral arm.

We present Monte Carlo scattered light models of a warped disk that reproduce the observed photopolarimetric variability of the classical T Tauri star AA Tau. For a system inclination of 75 deg and using an analytic description for a warped inner disk, we find that the shape and amplitude of the photopolarimetric variability are reproduced with a warp which occults the star, located at 0.07 au, amplitude 0.016 au, extending over radial and azimuthal ranges 0.0084 au and 145 deg. We also show a time sequence of high spatial resolution scattered light images, showing a dark shadow cast by the warp sweeping round the disk. Using a modified smooth particle hydrodynamics code, we find that a stellar dipole magnetic field of strength 5.2 kG, inclined at 30 deg to the stellar rotation axis can reproduce the required disk warping to explain AA Tau's photopolarimetric variability.

We present observations of four T Tauri stars using long baseline infrared interferometry from the Palomar Testbed Interferometer. The target sources, T Tau N, SU Aur, RY Tau and DR Tau, are all known to be surrounded by dusty circumstellar disks. The observations directly trace the inner regions (< 1 AU) of the disk and can be used to constrain the physical properties of this material. For three of the sources observed, the infrared emission is clearly resolved. We first use geometric models to characterize the emission region size, which ranges from 0.04 to 0.3 AU in radius. We then use Monte Carlo radiation transfer models of accretion disks to jointly model the spectral energy distribution and the interferometric observations with disk models including accretion and scattering. With these models, we are able to reproduce the data set with extended emission arising from structures larger than 10 milliarcseconds contributing less than 6% of the K band emission, consistent with there being little or no envelope remaining for these Class II sources (d log(lambda*F_lambda)/d log(lambda) ~ -2 to 0 in the infrared). The radiation transfer models have inner radii for the dust similar to the geometric models; however, for RY Tau emission from gas within the inner dust radius contributes significantly to the model flux and visibility at infrared wavelengths. The main conclusion of our modeling is that emission from inner gas disks (between the magnetic truncation radius and the dust destruction radius) can be a significant component in the inner disk flux for sources with large inner dust radii.

We determine and tabulate A(lambda)/A(K), the wavelength dependence of interstellar extinction, in the Galactic plane for 1.25um<lambda<8.0um along two lines of sight: l~42 and l~284. The first is a relatively quiescent and unremarkable region; the second contains the giant HII region RCW 49 as well as a ``field'' region unrelated to the cluster and nebulosity. Areas near these Galactic longitudes were imaged at J, H, and K bands by 2MASS and at 3--8um by Spitzer for the GLIMPSE Legacy program. We measure the mean values of the color excess ratios (A(lambda)-A(K))/(A(J)-A(K)) directly from the color distributions of observed stars. The extinction ratio between two of the filters, e.g. A(J)/A(K), is required to calculate A(lambda)/A(K) from those measured ratios. We use the apparent JHK magnitudes of giant stars along our two sightlines, and fit the reddening as a function of magnitude (distance) to determine A(J)/kpc, A(K)/kpc, and A(J)/A(K). Our values of A(lambda)/A(K) show a flattening across the 3--8um wavelength range, roughly consistent with the Lutz 1996 extinction measurements derived for the sightline toward the Galactic center.

Spitzer Space Telescope imaging as part of the Galactic Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) reveals a previously unidentified low-latitude rich star cluster near l=31.3 degrees, b=-0.1 degrees. Near-infrared JHK' photometry from the Wyoming Infrared Observatory indicates an extinction of A_V ~ 15+/-3 mag for cluster members. Analysis of 13CO features along the same sightline suggests a probable kinematic distance of 3.1 -- 5.2 kpc. The new cluster has an angular diameter of ~1-2 pc, a total magnitude m_{K_0}=2.1 corrected for extinction, and a luminosity of M_K ~ -10.3 at 3.1 kpc. In contrast to young massive Galactic clusters with ages <100 Myr, the new cluster has no significant radio emission. Comparison to theoretical K-band luminosity functions indicates an age of at least several Gyr and a mass of at least 10^5 solar masses. Unlike known old open clusters, this new cluster lies in the inner Galaxy at R_{GC} ~ 6.1 kpc. We designate this object ``GLIMPSE-C01'' and classify it as a Milky Way globular cluster passing through the Galactic disk. We also identify a region of star formation and fan-shaped outflows from young stellar objects in the same field as the cluster. The cluster's passage through the Galactic molecular layer may have triggered this star formation activity.

We model how the mid-infrared colors of Young Stellar Objects (YSOs) vary with stellar temperature. The spectral energy distribution (SED) of each object has contributions from thermal emission of circumstellar dust, from direct stellar photospheric emission, and from scattered stellar emission. We first isolate the effects of stellar contributions (direct+scattered) to the SED using homologous "Class I" models, where the shape of the thermal contribution to the SED remains constant. Stellar light contributes relatively more to the MIR emission of YSOs with lower temperature stars (Tstar ~4000 K) because the emission peak wavelength of the star is closer to that of the thermal radiation. Next we consider realistic Class 0, I, and II source models and find that the other dominant effect of varying stellar temperature on YSO SEDs is that of the inner disk wall: In high-Tstar models, the dust destruction radius is much further out with a consequently larger inner disk wall that contributes relatively more to the 2-10 um flux. This effect partially offsets that of the stellar contribution leading to varying behaviors of the 2-10 um flux. Our modeled MIR protostellar colors have implications for interpretations of Spitzer IRAC observations of star formation regions: It is commonly assumed that the slope of the SED at 1-10 um is directly related to evolutionary state. We show that inclination effects, aperture size, scattered light, and stellar temperature cause a broad spread in the colors of a source at a single evolutionary state. Color magnitude diagrams can help sort out these effects by separating sources with different Tstar based on their different brightness (for sources at the same distance).

GLIMPSE imaging using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope indicates that star formation is ongoing in the RCW 49 giant H II region. A photometric comparison of the sources in RCW 49 to a similar area to its north finds that at least 300 stars brighter than 13th magnitude in band [3.6] have infrared excesses inconsistent with reddening due to foreground extinction. These are likely young stellar objects (YSOs) more massive than 2.5 M_sun suggesting that thousands more low-mass stars are forming in this cloud. Some of the YSOs are massive (B stars) and therefore very young, suggesting that a new generation of star formation is occurring, possibly triggered by stellar winds and shocks generated by the older (2-3 Myr) central massive cluster. The Spitzer IRAC camera has proven to be ideally suited for distinguishing young stars from field stars, and the GLIMPSE survey of the Galactic Plane will likely find thousands of new star formation regions.

We present synthetic spectra for circumstellar disks that are heated by radiation from a central brown dwarf. Under the assumption of vertical hydrostatic equilibrium, our models yield scaleheights for brown dwarf disks in excess of three times those derived for classical T Tauri (CTTS) disks. If the near-IR excess emission observed from brown dwarfs is indeed due to circumstellar disks, then the large scaleheights we find could have a significant impact on the optical and near-IR detectability of such systems. Our radiation transfer calculations show that such highly flared disks around brown dwarfs will result in a large fraction of obscured sources due to extinction of direct starlight by the disk over a wide range of sightlines. The obscured fraction for a 'typical' CTTS is less than 20%. We show that the obscured fraction for brown dwarfs may be double that for CTTS, but this depends on stellar and disk mass. We also comment on possible confusion in identifying brown dwarfs via color-magnitude diagrams: edge-on CTTS display similar colors and magnitudes as a face-on brown dwarf plus disk systems.

We present images of the young stellar object ASR 41 in the NGC 1333 star forming region at the wavelengths of H_alpha and [SII] and in the I, J, H, and K-bands. ASR 41 has the near-infrared morphology of an edge-on disk object, but appears an order of magnitude larger than typical systems of this kind. We also present detailed models of the scattering and radiative transfer in systems consisting of a young star surrounded by a proto-planetary disk, and the whole system being embedded in either an infalling envelope or a uniform molecular cloud. The best fit to the observed morphology can be achieved with a disk of approx. 200 AU diameter, immersed in a low density cloud. The low cloud density is necessary to stay below the sub-mm flux upper limits and to preserve the shadow cast by the disk via single scattering. The results demonstrate that ASR 41 is probably not inherently different from typical edge-on disk objects, and that its large apparent size is due to the shadow of a much smaller disk being projected into the surrounding dusty molecular material.

We present model spectral energy distributions, colors, polarization, and images for an evolutionary sequence of a low-mass protostar from the early collapse stage (Class 0) to the remnant disk stage (Class III). We find a substantial overlap in colors and SEDs between protostars embedded in envelopes (Class 0-I) and T Tauri disks (Class II), especially at mid-IR wavelengths. Edge-on Class I-II sources show double-peaked spectral energy distributions, with a short-wavelength hump due to scattered light and the long-wavelength hump due to thermal emission. These are the bluest sources in mid-IR color-color diagrams. Since Class 0 and I sources are diffuse, the size of the aperture over which fluxes are integrated has a substantial effect on the computed colors, with larger aperture results showing significantly bluer colors. This causes overlap in color-color diagrams between all evolutionary states, especially in the mid-IR. However the near-IR polarization of the Class 0 sources is much higher than the Class I-II sources, providing a means to separate these evolutionary states. We varied the grain properties in the circumstellar envelope, allowing for larger grains in the disk midplane and smaller in the envelope. We find that grain growth in disks of Class I sources can be detected at wavelengths greater than 100 $\mu$m. Our image calculations predict that the diffuse emission from edge-on Class I and II sources should be detectable in the mid-IR with the Space Infrared Telescope Facility (SIRTF) in nearby star forming regions (out to several hundred parsecs).

We present the results of a Monte Carlo radiative transfer code that calculates the polarization produced by multiple Thomson scattering and variable absorptive opacity in a circumstellar disk around one component of a close detached binary system. We consider in detail the polarization variations over the binary cycle that result from the disk's illumination by the external star and by its own volume emission. We identify key features of these polarization phase curves and investigate their behavior as functions of optical depth, albedo, and inclination for geometrically thin and thick disks. The polarization due to disk self-illumination is sensitive to the internal optical characteristics of the disk, while the polarization arising from external illumination is mainly sensitive to the disk's geometrical thickness. With appropriate flux weighting, these results, combined with those for an internally illuminated disk, allow simulation of the polarization signature from an arbitrary binary-disk system.

GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire), a SIRTF Legacy Science Program, will be a fully sampled, confusion-limited infrared survey of the inner two-thirds of the Galactic disk with a pixel resolution of \~1.2" using the Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 microns. The survey will cover Galactic latitudes |b| <1 degree and longitudes |l|=10 to 65 degrees (both sides of the Galactic center). The survey area contains the outer ends of the Galactic bar, the Galactic molecular ring, and the inner spiral arms. The GLIMPSE team will process these data to produce a point source catalog, a point source data archive, and a set of mosaicked images. We summarize our observing strategy, give details of our data products, and summarize some of the principal science questions that will be addressed using GLIMPSE data. Up-to-date documentation, survey progress, and information on complementary datasets are available on the GLIMPSE web site: www.astro.wisc.edu/glimpse.

We present 2-D radiation transfer models of Class I Protostars and show the effect of including more realistic geometries on the resulting spectral energy distributions and images. We begin with a rotationally flattened infalling envelope as our comparison model, and add a flared disk and bipolar cavity. The disk affects the spectral energy distribution most strongly at edge-on inclinations, causing a broad dip at about 10 um (independent of the silicate feature) due to high extinction and low scattering albedo in this wavelength region. The bipolar cavities allow more direct stellar+disk radiation to emerge into polar directions, and more scattering radiation to emerge into all directions. The wavelength-integrated flux, often interpreted as luminosity, varies with viewing angle, with pole-on viewing angles seeing 2-4 times as much flux as edge-on, depending on geometry. Thus, observational estimates of luminosity should take into account the inclination of a source. The envelopes with cavities are significantly bluer in near-IR and mid-IR color-color plots than those without cavities. Using 1-D models to interpret Class I sources with bipolar cavities would lead to an underestimate of envelope mass and an overestimate of the implied evolutionary state. We compute images at near-, mid-, and far-IR wavelengths. We find that the mid-IR colors and images are sensitive to scattering albedo, and that the flared disk shadows the midplane on large size scales at all wavelengths plotted. Finally, our models produce polarization spectra which can be used to diagnose dust properties, such as albedo variations due to grain growth. Our results of polarization across the 3.1 um ice feature agree well with observations for ice mantles covering 5% of the radius of the grains.

[Abridged] We present ISO far-infrared (IR) observations of 21 hard X-ray selected AGN from the HEAO-1 A2 sample. We compare the far-IR to X-ray spectral energy distributions (SEDs) of this sample with various radio and optically selected AGN samples. The hard-X-ray selected sample shows a wider range of optical/UV shapes extending to redder near-IR colors. The bluer objects are Seyfert 1s, while the redder AGN are mostly intermediate or type 2 Seyferts. This is consistent with a modified unification model in which the amount of obscuring material increases with viewing angle and may be clumpy. Such a scenario, already suggested by differing optical/near-IR spectroscopic and X-ray AGN classifications, allows for different amounts of obscuration of the continuum emission in different wavebands and of the broad emission line region which results in a mixture of behaviors for AGN with similar optical emission line classifications. The resulting limits on the column density of obscuring material through which we are viewing the redder AGN are 100 times lower than for the standard optically thick torus models. The resulting decrease in optical depth of the obscuring material allows the AGN to heat more dust at larger radial distances. We show that an AGN-heated, flared, dusty disk with mass 10^9 solar and size of few hundred pc is able to generate optical-far-IR SEDs which reproduce the wide range of SEDs present in our sample with no need for an additional starburst component to generate the long-wavelength, cooler part of the IR continuum.

The unusual spectral energy distribution (SED) of the classical T Tauri star GM Aurigae provides evidence for the presence of an inner disc hole extending to several au. Using a combination of hydrodynamical simulations and Monte Carlo radiative transport, we investigate whether the observed SED is consistent with the inner hole being created and maintained by an orbiting planet. We show that an ~ 2 M_Jupiter planet, orbiting at 2.5 au in a disc with mass 0.047 M_sun and radius 300 au, provides a good match both to the SED and to CO observations which constrain the velocity field in the disc. A range of planet masses is allowed by current data, but could in principle be distinguished with further observations between 3 and ~ 20 microns. Future high precision astrometric instruments should also be able to detect the motion of the central star due to an orbiting Jupiter mass planet. We argue that the small number of T Tauri stars with SEDs resembling that of GM Aur is broadly consistent with the expected statistics of embedded migrating planets.

By radiation transfer models with a realistic power spectra of the projected density distributions, we show that the optical properties of grains are poorly constrained by observations of reflection nebulae. The ISM is known to be hierarchically clumped from a variety of observations (molecules, H I, far-infrared). Our models assume the albedo and phase parameter of the dust, the radial optical depth of the sphere averaged over all directions, and random distributions of the dust within the sphere. The outputs are the stellar extinction, optical depth, and flux of scattered light as seen from various viewing angles. Observations provide the extinction and scattered flux from a particular direction. Hierarchical geometry has a large effect on the flux of scattered light emerging from a nebula for a particular extinction of the exciting star. There is a very large spread in both scattered fluxes and extinctions for any distribution of dust. Consequently, an observed stellar extinction and scattered flux can be fitted by a wide range of albedos. With hierarchical geometry it is not completely safe to determine even relative optical constants from multiwavelength observations of the same reflection nebula. The geometry effectively changes with wavelength as the opacity of the clumps varies. Limits on the implications of observing the same object in various wavelengths are discussed briefly. Henry (2002) uses a recipe to determine the scattered flux from a star with a given extinction. It is claimed to be independent of the geometry. It provides considerably more scattering than our models, probably leading to an underestimate of the grain albedos from the UV Diffuse Galactic Light.

We investigate the observational signatures of a straightforward evolutionary scenario for protoplanetary disks in which the disk mass of small (50 micron) particles decreases homologously with time, but the disk structure and stellar parameters do not change. Our goal is to identify optimal infrared spectral indicators of the existence of disks, their structure, and mass evolution that may be tested with the upcoming SIRTF mission. We present simulated spectral energy distributions and colors over a wide range of masses. The SED is most sensitive to disk mass in the far-IR and longer wavelengths, which is already known from millimeter and radio observations. As the disk mass decreases, the excess emission of the disk over the stellar photosphere diminishes more rapidly at the longest rather than at short wavelengths. At near-infrared wavelengths, the disk remains optically thick to stellar radiation over a wide range of disk mass, resulting in a slower decline in the SED in this spectral regime. Therefore, near-IR excesses (K-L) provide a robust means of detecting disks in star clusters down to 1E-7 solar masses, while the far-IR excess probes the disk mass. Reducing the disk mass results in a clear progression in color-color diagrams with low mass disks displaying the bluest colors. We interpret color-color diagrams for Taurus-Auriga sources in the context of decreasing disk mass.

We present spectral energy distribution (SED) models for the edge-on classical T Tauri star HH30 IRS that indicate dust grains have grown to larger than 50 microns within its circumstellar disk. The disk geometry and inclination are known from previous modeling of multiwavelength Hubble Space Telescope images and we use the SED to constrain the dust size distribution. Model spectra are shown for different circumstellar dust models: a standard ISM mixture and larger grain models. As compared to ISM grains, the larger dust grain models have a shallower wavelength dependent opacity. Models with the larger dust grains provide a good match to the currently available data, but mid and far-IR observations are required to more tightly constrain the dust size distribution. The accretion luminosity in our models is L_acc<0.2 L_star corresponding to an accretion rate of 4E-9M_sun/yr. Dust size distributions that are simple power-law extensions (i.e., no exponential cutoff) yield acceptable fits to the optical/near-IR but too much emission at mm wavelengths and require larger disk masses. Such a simple size distribution would not be expected in an environment such as the disk of HH30 IRS, particularly over such a large range in grain sizes. However, its ability to adequately characterize the grain populations may be determined from more complete observational sampling of the SED in the mid to far-IR.

Ground based imaging, imaging polarimetry, and recent Hubble Space Telescope WFPC2 and NICMOS images of protostars have revealed very complex scattered light patterns that cannot be entirely explained by 2-D radiation transfer models. We present here for the first time radiation transfer models of T Tau and IRAS 04016+2610 that are fully 3-D, with the aim of investigating the effects on image morphology of multiple illuminating sources and infalling envelopes that have been shaped by multiple outflows. For T Tau we have constructed scattered light models where the illumination of the surrounding envelope is by a binary with each source surrounded by its own small circumstellar disk or envelope. We find that the asymmetries in the WFPC2 image of T Tau can be reproduced if the disks in the binary system are misaligned, consistent with a recently discovered bipolar outflow believed to originate from the secondary. For IRAS 04016+2610 we find that the observed scattered light pattern can be reproduced by scattering in an envelope with cavities carved by two sets of bipolar outflows, suggestive of an embedded binary system.

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.

We report results of VRI photometric monitoring of the T Tauri star plus disk system HH30. We find that HH30 is highly variable over timescales of a few days with Delta V~1.5 mag, Delta I~1.1 mag. Furthermore we find hints of periodicity with periodograms indicating possible periods of 11.6 and 19.8 days. The VRI photometry is available through the anonymous ftp service. We model the variability with Monte Carlo radiation transfer simulations for a spotted star plus disk system and find that the large variability is best reproduced with a single hot spot and circumstellar grains that are larger than typical interstellar grains. The apparent existence of a single hot spot and the need for larger, more forward throwing grains is consistent with previous modeling of HST imagery.

The close-in extrasolar giant planets [CEGPs], $\ltorder$ 0.05 AU from their parent stars, may have a large component of optically reflected light. We present theoretical optical photometric light curves and polarization curves for the CEGP systems, from reflected planetary light. Different particle sizes of three condensates are considered. In the most reflective case, the variability is $\approx 100$ micromagnitudes, which will be easily detectable by the upcoming satellite missions MOST, COROT, and MONS, and possibly from the ground in the near future. The least reflective case is caused by small, highly absorbing grains such as solid Fe, with variation of much less than one micromagnitude. Polarization for all cases is lower than current detectability limits. We also discuss the temperature-pressure profiles and resulting emergent spectra of the CEGP atmospheres. We discuss the observational results of Tau Boo b by Cameron et al. (1999) and Charbonneau et al. (1999) in context of our model results. The predictions - the shape and magnitude of the light curves and polarization curves - are highly dependent on the size and type of condensates present in the planetary atmosphere.

We describe planetesimal accretion calculations for the dusty ring observed in the nearby A0 star HR 4796A. Models with initial masses of 10-20 times the minimum mass solar nebula produce a ring of width 7-15 AU and height 0.3-0.6 AU at 70 AU in roughly 10 Myr. The ring has a radial optical depth of 1. These results agree with limits derived from infrared images and from the excess infrared luminosity.