Infrared to millimetre photometry of ultra-luminous IR galaxies: new evidence favouring a 3-stage dust model

Infrared to millimetre spectral energy distributions have been obtained for 41 bright ultra-luminous infrared galaxies. The observations were carried out with ISOPHOT between 10 and 200 micron and supplemented for 16 sources with SCUBA at 450 and 850 micron and with SEST at 1.3 mm. In addition, seven sources were observed at 1.2 and 2.2 $\mu$m with the 2.2 m telescope on Calar Alto. These new SEDs represent the most complete set of infrared photometric templates obtained so far on ULIRGs in the local universe.Comment: 23 pages, 11 figures, accepted for publication in Astronomy & Astrophysic

HD 152246 - a new high-mass triple system and its basic properties

Analyses of multi-epoch, high-resolution (R ~ 50.000) optical spectra of the O-type star HD 152246 (O9 IV according to the most recent classification), complemented by a limited number of earlier published radial velocities, led to the finding that the object is a hierarchical triple system, where a close inner pair (Ba-Bb) with a slightly eccentric orbit (e = 0.11) and a period of 6.0049 days revolves in a 470-day highly eccentric orbit (e = 0.865) with another massive and brighter component A. The mass ratio of the inner system must be low since we were unable to find any traces of the secondary spectrum. The mass ratio A/(Ba+Bb) is 0.89. The outer system has recently been resolved using long-baseline interferometry on three occasions. The interferometry confirms the spectroscopic results and specifies elements of the system. Our orbital solutions, including the combined radial-velocity and interferometric solution indicate an orbital inclination of the outer orbit of 112{\deg} and stellar masses of 20.4 and 22.8 solar masses. We also disentangled the spectra of components A and Ba and compare them to synthetic spectra from two independent programmes, TLUSTY and FASTWIND. In either case, the fit was not satisfactory and we postpone a better determination of the system properties for a future study, after obtaining observations during the periastron passage of the outer orbit (the nearest chance being March 2015). For the moment, we can only conclude that component A is an O9 IV star with v*sin(i) = 210 +\- 10 km/s and effective temperature of 33000 +\- 500 K, while component Ba is an O9 V object with v*sin(i) = 65 +/- 3 km/s and T_eff = 33600 +\- 600 K.Comment: 9 pages, 6 figures, accepted for publication in Astronomy and Astrophysic

3-D Models of Embedded High-Mass Stars: Effects of a Clumpy Circumstellar Medium

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.Comment: accepted to ApJ; version with full-res figures: http://www.astro.virginia.edu/~ri3e/clumpy3d.pd

Infall models of Class 0 protostars

We have carried out radiative transfer calculations of infalling, dusty envelopes surrounding embedded protostars to understand the observed properties of the recently identified ``Class 0'' sources. To match the far-infrared peaks in the spectral energy distributions of objects such as the prototype Class 0 source VLA 1623, pure collapse models require mass infall rates \sim10^{-4}\msunyr$^{-1}$. The radial intensity distributions predicted by such infall models are inconsistent with observations of VLA 1623 at sub-mm wavelengths, in agreement with the results of Andre et al. (1993) who found a density profile of $\rho \propto r^{-1/2}$ rather than the expected $\rho \propto r^{-3/2}$ gradient. To resolve this conflict, while still invoking infall to produce the outflow source at the center of VLA 1623, we suggest that the observed sub-mm intensity distribution is the sum of two components: an inner infall zone, plus an outer, more nearly constant-density region. This explanation of the observations requires that roughly half the total mass observed within 2000 AU radius of the source lies in a region external to the infall zone. The column densities for this external region are comparable to those found in the larger Oph A cloud within which VLA 1623 is embedded. The extreme environments of Class 0 sources lead us to suggest an alternative or additional interpretation of these objects: rather than simply concluding with Andre et al. that Class 0 objects only represent the earliest phases of protostellar collapse, and ultimately evolve into older ``Class I'' protostars, we suggest that many Class 0 sources could be the protostars of very dense regions. (Shortened)Comment: 22 pages, including 3 PostScript figures, accepted for publication in The Astrophysical Journa

Generalized Quasilinear Approximation: Application to Zonal Jets

Quasilinear theory is often utilized to approximate the dynamics of fluids exhibiting significant interactions between mean flows and eddies. We present a generalization of quasilinear theory to include dynamic mode interactions on the large scales. This generalized quasilinear (GQL) approximation is achieved by separating the state variables into large and small zonal scales via a spectral filter rather than by a decomposition into a formal mean and fluctuations. Nonlinear interactions involving only small zonal scales are then removed. The approximation is conservative and allows for scattering of energy between small-scale modes via the large scale (through nonlocal spectral interactions). We evaluate GQL for the paradigmatic problems of the driving of large-scale jets on a spherical surface and on the beta plane and show that it is accurate even for a small number of large-scale modes. As GQL is formally linear in the small zonal scales, it allows for the closure of the system and can be utilized in direct statistical simulation schemes that have proved an attractive alternative to direct numerical simulation for many geophysical and astrophysical problems

High Mass Star Formation. II. The Mass Function of Submillimeter Clumps in M17

We have mapped an approximately 5.5 by 5.5 pc portion of the M17 massive star-forming region in both 850 and 450 micron dust continuum emission using the Submillimeter Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). The maps reveal more than 100 dusty clumps with deconvolved linear sizes of 0.05--0.2 pc and masses of 0.8--120 solar masses, most of which are not associated with known mid-infrared point sources. Fitting the clump mass function with a double power law gives a mean power law exponent of alpha_high = -2.4 +/- 0.3 for the high-mass power law, consistent with the exponent of the Salpeter stellar mass function. We show that a lognormal clump mass distribution with a peak at about 4 solar masses produces as good a fit to the clump mass function as does a double power law. This 4 solar mass peak mass is well above the peak masses of both the stellar initial mass function and the mass function of clumps in low-mass star-forming regions. Despite the difference in intrinsic mass scale, the shape of the M17 clump mass function appears to be consistent with the shape of the core mass function in low-mass star-forming regions. Thus, we suggest that the clump mass function in high-mass star-forming regions may be a scaled-up version of that in low-mass regions, instead of its extension to higher masses.Comment: 33 pages, 6 figures, 3 tables. Accepted for publication in the Astrophysical Journa

Properties of the Narrow-Line Region in Seyfert Galaxies

We study the narrow-line region (NLR) of six Seyfert-1 and six Seyfert-2 galaxies by means of spatially resolved optical spectroscopy and photoionization modelling. From spatially resolved spectral diagnostics, we find a transition between the AGN-excited NLR and the surrounding star-forming regions, allowing us to determine the NLR size independent of stellar contamination. CLOUDY photoionization models show that the observed transition represents a true difference in ionization source and cannot be explained by variations of physical parameters. The electron density and ionization parameter decrease with radius indicating that the NLR is photoionized by the central source only. The velocity field suggests a disky NLR gas distribution

Structure and Thermodynamics of the Mixed Alkali Alanates

The thermodynamics and structural properties of the hexahydride alanates (M2M'AlH6) with the elpasolite structure have been investigated. A series of mixed alkali alanates (Na2LiAlH6, K2LiAlH6 and K2NaAlH6) were synthesized and found to reversibly absorb and desorb hydrogen without the need for a catalyst. Pressure-composition isotherms were measured to investigate the thermodynamics of the absorption and desorption reactions with hydrogen. Isotherms for catalyzed (4 mol% TiCl3) and uncatalyzed Na2LiAlH6 exhibited an increase in kinetics, but no change in the bulk thermodynamics with the addition of a dopant. A structural analysis using synchrotron x-ray diffraction showed that these compounds favor the Fm-3m space group with the smaller ion (M') occupying an octahedral site. These results demonstrate that appropriate cation substitutions can be used to stabilize or destabilize the material and may provide an avenue to improving the unfavorable thermodynamics of a number of materials with promising gravimetric hydrogen densities.Comment: 6 pages, 7 figures,3 tables, submitted to PR

Freshly Formed Dust in the Cassiopeia A Supernova Remnant as Revealed by the Spitzer Space Telescope

We performed Spitzer Infrared Spectrograph mapping observations covering nearly the entire extent of the Cassiopeia A supernova remnant (SNR), producing mid-infrared (5.5-35 micron) spectra every 5-10". Gas lines of Ar, Ne, O, Si, S and Fe, and dust continua were strong for most positions. We identify three distinct ejecta dust populations based on their continuum shapes. The dominant dust continuum shape exhibits a strong peak at 21 micron. A line-free map of 21 micron-peak dust made from the 19-23 micron range closely resembles the [Ar II], [O IV], and [Ne II] ejecta-line maps implying that dust is freshly formed in the ejecta. Spectral fitting implies the presence of SiO2, Mg protosilicates, and FeO grains in these regions. The second dust type exhibits a rising continuum up to 21 micron and then flattens thereafter. This ``weak 21 micron'' dust is likely composed of Al2O3 and C grains. The third dust continuum shape is featureless with a gently rising spectrum and is likely composed of MgSiO3 and either Al2O3 or Fe grains. Using the least massive composition for each of the three dust classes yields a total mass of 0.02 Msun. Using the most-massive composition yields a total mass of 0.054 Msun. The primary uncertainty in the total dust mass stems from the selection of the dust composition necessary for fitting the featureless dust as well as 70 micron flux. The freshly formed dust mass derived from Cas A is sufficient from SNe to explain the lower limit on the dust masses in high redshift galaxies.Comment: 8 figures: Accepted for the publication in Ap