Search

It has been advanced, on experimental (P.-M. Robitaille, IEEE Trans. Plasma Sci. 2003, v. 31(6), 1263-1267) and theoretical (P.M. Robitaille, Progr. Phys. 2006, v.2, 22-23) grounds, that blackbody radiation is not universal and remains closely linked to the emission of graphite and soot. In order to strengthen such claims, a conceptual analysis of the proofs for universality is presented. This treatment reveals that Gustav Robert Kirchhoff has not properly considered the combined effects of absorption, reflection, and the directional nature of emission in real materials. In one instance, this leads to an unintended movement away from thermal equilibrium within cavities. Using equilibrium arguments, it is demonstrated that the radiation within perfectly reflecting or arbitrary cavities does not necessarily correspond to that emitted by a blackbody.

In this contribution, I review the applications and potential limitations of the spectral energy distribution fitting tool that I have developed, with a strong emphasis on the limits to which this tool can be used to improve our understanding of massive star formation. I discuss why our current grid of models cannot be used to distinguish between the several competing theories of massive star formation. I also discuss stellar mass determinations, artificial correlations between parameters in the grid of models, multiplicity, confusion, dust assumptions, and unique fits. I briefly review the improvements we intend to carry out for our next grid of models, which will eliminate many of these limitations. Finally, I show examples of applications of this tool to massive young stars., Comment: 8 Pages, 3 Figures, to appear in the proceedings for "Massive Star Formation: Observations confront Theory 2007"

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., Comment: 10 pages; 3 figures; 3 tables

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., Comment: 47 pages, 12 Figures, 3 Tables, accepted for publication in AJ

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., Comment: 13 pages, 5 figures (3 in color), 2 tables. Accepted for publication in ApJ

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., Comment: 16p preprint - ApJ accepted

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., Comment: Accepted by ApJ, 46 pages, 12 figures

(Abridged) Spiral galaxies often have extended outflows that permeate beyond the region of the disk. Such outflows have been seen both in starburst galaxies, actively star forming galaxies and galaxies with an AGN. In the latter galaxies it is unknown whether the large-scale outflows are driven by star formation activity or purely by the active nucleus. The aim of our investigation is to study the frequency of extended minor-axis outflows in edge-on Seyfert galaxies to investigate the role of the AGN, the circumnuclear environment and star formation activity within the disk regions, and their importance for IGM enrichment on large scales. Narrowband imaging in two different ionizational stages (H-alpha and [OIII]) was performed to attempt a discrimination between processes associated with the active nucleus and those connected to star forming activity within the disk. The H-alpha morphology of the Seyfert galaxies is usually complex, but only in three out of 14 galaxies did we find evidence for minor axis disk outflows. At the sensitivity of our observations [OIII] emission is generally detected only in the nuclear region. Overall, our results show that extraplanar emission of similar brightness and extent as in the previously known cases of NGC3079 and NGC4388 is not common in Seyfert galaxies of otherwise similar properties. Comparison with our previous results shows that for nearby edge-on spiral galaxies star formation may be a more powerful mechanism for producing DIG than AGN activity. While in general AGN activity undoubtedly plays some role in driving minor-axis outflows, this probably requires higher AGN luminosities than are encountered in our small distance-limited sample., Comment: accepted for publication in Astronomy & Astrophysics, 36 pages, 25 figures, 5 tables (Figs. 6-24 and Tables 2 + 5 in the online part), Article with Full resolution images is available at http://www.stsci.edu/~marel/psgzdir/seyfert_outflows.ps.gz

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 ., Comment: 54 pages, 7 figures, 12 tables. Accepted for publication in ApJS. Preprint with full resolution figures available 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%., Comment: Accepted to ApJ

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. [...]., Comment: 69 pages, 28 figures, Accepted for publication in ApJS. Preprint with full resolution figures available at http://www.astro.wisc.edu/protostars/

We present numerical simulations of the passage of clumpy gas through a galactic spiral shock and the subsequent formation of giant molecular clouds (GMCs). The spiral shock forms dense clouds while dissipating kinetic energy, producing regions that are locally gravitationally bound and collapse to form stars. The effect of the clumpiness of gas as it passes through the shock is to generate chaotic internal motions in the gas. The kinematics of these motions are found to agree with the observed velocity-dispersion/size relation found in star-forming regions. In contrast to the standard picture where continuously driven turbulence generates the density inhomogeneities in star-forming clouds, we find here that it is the clumpiness of the interstellar gas that produces the chaotic motions as it passes through the spiral shock and initiates the star formation process. The velocity dispersion can be understood as being due to the random mass loading of clumps as they converge in the spiral shock. In this model there is no need for any internal or external continuous driving mechanism for the 'turbulence'. The coupling of the clouds' internal kinematics to their externally triggered formation removes the need for the clouds to be self-gravitating. Indeed, while clearly some parts of the clouds are self-gravitating and able to form stars, most of the molecular material remains gravitationally unbound. This can provide a simple explanation for the low efficiency of star formation., Comment: 9 pages, 6 figures, accepted for publication in MNRAS

Through the formulation of his law of thermal emission, Kirchhoff conferred upon blackbody radiation the quality of universality [G.Kirchhoff, Annalen der Physik 109, 275 (1860)]. Consequently, modern physics holds that such radiation is independent of the nature and shape of the emitted object. Recently, Kirchhoff's experimental work and theoretical conclusions have been reconsidered [P.M.L. Robitaille, IEEE Transactions on Plasma Science 31(6), 1263 (2003). In this work, Einstein's derivation of the Planckian relation is reexamined. It is demonstrated that claims of universality in blackbody radiation are invalid., Comment: Presented at the New England Section Meeting of the APS on April 2, 2005 at MIT

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., Comment: 10 pages, 8 figures, to appear in the proceedings of IAU Symp 227, Massive Star Birth: A Crossroads of Astrophysics, (Cesaroni R., Churchwell E., Felli M., Walmsley C. editors)

In this work, a liquid model of the Sun is presented wherein the entire solar mass is viewed as a high density/high energy plasma. This model challenges our current understanding of the densities associated with the internal layers of the Sun, advocating a relatively constant density, almost independent of radial position. The incompressible nature of liquids is advanced to prevent solar collapse from gravitational forces. The liquid plasma model of the Sun is a non-equilibrium approach, where nuclear reactions are free to occur throughout the solar mass. The primary means of addressing internal heat transfer are convection and conduction. As a result of the convective processes on the solar surface, the liquid model brings into question the established temperature of the solar photosphere by highlighting a violation of Kirchhoff's law of thermal emission. Along these lines, the model emphasizes that radiative emission is a surface phenomenon. Strong evidence is provided that the Sun is a high density/high energy liquid plasma. This evidence is based on our knowledge of Planckian thermal emission and condensed matter, including the existence of pressure ionization and liquid metallic hydrogen at high temperatures and pressures. The equations of magnetohydrodynamics are invoked as the proper vehicle for the understanding od stellar convection and structure. Prior to introducing the liquid plasma model, the historic and scientific justifications for the gaseous model of the Sun are reviewed and the gaseous equations of state are also discussed., Comment: corrected typos in abstract

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., Comment: 13 pages, 13 figures, accepted for publication in MNRAS

The trend towards ever larger telescopes and more advanced adaptive optics systems is driving the need for deformable mirrors with a large number of low cost actuators. Liquid mirrors have long been recognized a potential low cost alternative to conventional solid mirrors. By using a water or oil based ferrofluid we are able to benefit from a stronger magnetic response than is found in magnetic liquid metal amalgams and avoid the difficulty of passing a uniform current through a liquid. Depositing a thin silver colloid known as a metal liquid-like film (MELLF) on the ferrofluid surface solves the problem of low reflectivity of pure ferrofluids. This combination provides a liquid optical surface that can be precisely shaped in a magnetic field. We present experimental results obtained with a prototype deformable liquid mirror based on this combination., Comment: 9 pages. To appear in SPIE meeting in Hawaii August 2002. Astronomical Telescopes and Instrumentation

This paper concerns the enumeration of rotation-type and congruence-type convex polyominoes on the square lattice. These can be defined as orbits of the groups C4, of rotations, and D4, of symmetries of the square acting on (translation- type) polyominoes. In virtue of Burnside's Lemma, it is sufficient to enumerate the various symmetry classes (fixed points) of polyominoes defined by the elements of C4 and D4. Using the Temperley--Bousquet-Melou methodology, we solve this problem and provide explicit or recursive formulas for their generating functions according to width, height and area. We also enumerate the class of asymmetric convex polyominoes, using Moebius inversion, and prove that their number is asymptotically equivalent to the number of convex polyominoes, a fact which is empirically evident., Comment: LaTeX 2e, 31 pages with 22 figures. Accepted in Advances in Applied Mathematics