Spitzer IRS Observations of the Galactic Center: Shocked Gas in the Radio Arc Bubble

We present Spitzer IRS spectra (R ~600, 10 - 38 micron) of 38 positions in the Galactic Center (GC), all at the same Galactic longitude and spanning plus/minus 0.3 degrees in latitude. Our positions include the Arches Cluster, the Arched Filaments, regions near the Quintuplet Cluster, the ``Bubble'' lying along the same line-of-sight as the molecular cloud G0.11-0.11, and the diffuse interstellar gas along the line-of-sight at higher Galactic latitudes. From measurements of the [O IV], [Ne II], [Ne III], [Si II], [S III], [S IV], [Fe II], [Fe III], and H_2 S(0), S(1), and S(2) lines we determine the gas excitation and ionic abundance ratios. The Ne/H and S/H abundance ratios are ~ 1.6 times that of the Orion Nebula. The main source of excitation is photoionization, with the Arches Cluster ionizing the Arched Filaments and the Quintuplet Cluster ionizing the gas nearby and at lower Galactic latitudes including the far side of the Bubble. In addition, strong shocks ionize gas to O^{+3} and destroy dust grains, releasing iron into the gas phase (Fe/H ~ 1.3 times 10^{-6} in the Arched Filaments and Fe/H ~ 8.8 times 10^{-6} in the Bubble). The shock effects are particularly noticeable in the center of the Bubble, but O$^{+3}$ is present in all positions. We suggest that the shocks are due to the winds from the Quintuplet Cluster Wolf-Rayet stars. On the other hand, the H_2 line ratios can be explained with multi-component models of warm molecular gas in photodissociation regions without the need for H_2 production in shocks.Comment: 51 pages, 17 figures To be published in the Astrophysical Journa

The coronal line regions of planetary nebulae NGC6302 and NGC6537: 3-13um grating and echelle spectroscopy

We report on advances in the study of the cores of NGC6302 and NGC6537 using infrared grating and echelle spectroscopy. In NGC6302, emission lines from species spanning a large range of ionization potential, and in particular [SiIX]3.934um, are interpreted using photoionization models (including CLOUDY), which allow us to reestimate the central star's temperature to be about 250000K. All of the detected lines are consistent with this value, except for [AlV] and [AlVI]. Aluminium is found to be depleted to one hundredth of the solar abundance, which provides further evidence for some dust being mixed with the highly ionized gas (with photons harder than 154eV). A similar depletion pattern is observed in NGC6537. Echelle spectroscopy of IR coronal ions in NGC6302 reveals a stratified structure in ionization potential, which confirms photoionization to be the dominant ionization mechanism. The lines are narrow (< 22km/s FWHM), with no evidence of the broad wings found in optical lines from species with similar ionization potentials, such as [NeV]3426A. We note the absence of a hot bubble, or a wind blown bipolar cavity filled with a hot plasma, at least on 1'' and 10km/s scales. We also provide accurate new wavelengths for several of the infrared coronal lines observed with the echelle.Comment: Accepted for publication in MNRA

Stellar Populations and Mass-Loss in M15: A Spitzer Detection of Dust in the Intra-Cluster Medium

We present Spitzer Space Telescope IRAC and MIPS observations of the galactic globular cluster M15 (NGC 7078), one of the most metal-poor clusters with a [Fe/H] = -2.4. Our Spitzer images reveal a population of dusty red giants near the cluster center, a previously detected planetary nebula (PN) designated K648, and a possible detection of the intra-cluster medium (ICM) arising from mass loss episodes from the evolved stellar population. Our analysis suggests 9 (+/-2) x 10^-4 solar masses of dust is present in the core of M15, and this material has accumulated over a period of approximately 10^6 years, a timescale ten times shorter than the last galactic plane crossing event. We also present Spitzer IRS follow up observations of K648, including the detection of the [NeII] 12.81 micron line, and discuss abundances derived from infrared fine structure lines.Comment: Accepted for publication in AJ. 20 pages, 10 figures, 6 tables. Full resolution versions of figures 1, 5, 7, and 8 are available in a PDF version of this manuscript at http://ir.astro.umn.edu/~mboyer/ms_060906.pd

A large-scale R-matrix calculation for electron-impact excitation of the Ne2+^{2+} O-like ion

The five J$\Pi$ levels within a $np^2$ or $np^4$ ground state complex provide an excellent testing ground for the comparison of theoretical line ratios with astrophysically observed values, in addition to providing valuable electron temperature and density diagnostics. The low temperature nature of the line ratios ensure that the theoretically derived values are sensitive to the underlying atomic structure and electron-impact excitation rates. Previous R-matrix calculations for the Ne$^{2+}$ O-like ion exhibit large spurious structure in the cross sections at higher electron energies, which may affect Maxwellian averaged rates even at low temperatures. Furthermore, there is an absence of comprehensive excitation data between the excited states that may provide newer diagnostics to compliment the more established lines discussed in this paper. To resolve these issues, we present both a small scale 56-level Breit-Pauli (BP) calculation and a large-scale 554 levels R-matrix Intermediate Coupling Frame Transformation (ICFT) calculation that extends the scope and validity of earlier JAJOM calculations both in terms of the atomic structure and scattering cross sections. Our results provide a comprehensive electron-impact excitation data set for all transitions to higher $n$ shells. The fundamental atomic data for this O-like ion is subsequently used within a collisional radiative framework to provide the line ratios across a range of electron temperatures and densities of interest in astrophysical observations.Comment: 17 pages, 8 figure

Quantitative analysis of WC stars: Constraints on neon abundances from ISO/SWS spectroscopy

Neon abundances are derived in four Galactic WC stars -- gamma Vel (WR11, WC8+O7.5III), HD156385 (WR90, WC7), HD192103 (WR135, WC8), and WR146 (WC5+O8) - using mid-infrared fine structure lines obtained with ISO/SWS. Stellar parameters for each star are derived using a non-LTE model atmospheric code (Hillier & Miller 1998) together with ultraviolet (IUE), optical (INT, AAT) and infrared (UKIRT, ISO) spectroscopy. In the case of gamma Vel, we adopt results from De Marco et al. (2000), who followed an identical approach. ISO/SWS datasets reveal the [NeIII] 15.5um line in each of our targets, while [NeII] 12.8um, [SIV] 10.5um and [SIII] 18.7um are observed solely in gamma Vel. Using a method updated from Barlow et al. (1988) to account for clumped winds, we derive Ne/He=3-4x10^-3 by number, plus S/He=6x10^-5 for gamma Vel. Neon is highly enriched, such that Ne/S in gamma Vel is eight times higher than cosmic values. However, observed Ne/He ratios are a factor of two times lower than predictions of current evolutionary models of massive stars. An imprecise mass-loss and distance were responsible for the much greater discrepancy in neon content identified by Barlow et al. Our sample of WC5--8 stars span a narrow range in T* (=55--71kK), with no trend towards higher temperature at earlier spectral type, supporting earlier results for a larger sample by Koesterke & Hamann (1995). Stellar luminosities range from 100,000 to 500,000 Lo, while 10^-5.1 < Mdot/(Mo/yr) < 10^-4.5, adopting clumped winds, in which volume filling factors are 10%. In all cases, wind performance numbers are less than 10, significantly lower than recent estimates. Carbon abundances span 0.08 < C/He < 0.25 by number, while oxygen abundances remain poorly constrained.Comment: 16 pages,7 figures accepted for MNRA

The Continuing Slow Decline of AG Pegasi

We analyze optical and ultraviolet observations of the symbiotic binary AG Pegasi acquired during 1992-97. The bolometric luminosity of the hot component declined by a factor of 2-3 from 1980-1985 to 1997. Since 1992, the effective temperature of the hot component may have declined by 10%-20%, but this decline is comparable to the measurement errors. Optical observations of H-beta and He I emission show a clear illumination effect, where high energy photons from the hot component ionize the outer atmosphere of the red giant. Simple illumination models generally account for the magnitude of the optical and ultraviolet emission line fluxes. High ionization emission lines - [Ne V], [Mg V], and [Fe VII] - suggest mechanical heating in the outer portions of the photoionized red giant wind. This emission probably originates in a low density region $\sim$ 30-300 AU from the central binary.Comment: 17 pages, 7 pages, 5 tables; to be published in the Astronomical Journal, July 200

Spitzer Quasar and ULIRG Evolution Study (QUEST). IV. Comparison of 1-Jy Ultraluminous Infrared Galaxies with Palomar-Green Quasars

We report the results from a comprehensive study of 74 ultraluminous infrared galaxies (ULIRGs) and 34 Palomar-Green (PG) quasars within z ~ 0.3$ observed with the Spitzer Infrared Spectrograph (IRS). The contribution of nuclear activity to the bolometric luminosity in these systems is quantified using six independent methods that span a range in wavelength and give consistent results within ~ +/-10-15% on average. The average derived AGN contribution in ULIRGs is ~35-40%, ranging from ~15-35% among "cool" (f_25/f_60 =< 0.2) optically classified HII-like and LINER ULIRGs to ~50 and ~75% among warm Seyfert 2 and Seyfert 1 ULIRGs, respectively. This number exceeds ~80% in PG QSOs. ULIRGs fall in one of three distinct AGN classes: (1) objects with small extinctions and large PAH equivalent widths are highly starburst-dominated; (2) systems with large extinctions and modest PAH equivalent widths have larger AGN contributions, but still tend to be starburst-dominated; and (3) ULIRGs with both small extinctions and small PAH equivalent widths host AGN that are at least as powerful as the starbursts. The AGN contributions in class 2 ULIRGs are more uncertain than in the other objects, and we cannot formally rule out the possibility that these objects represent a physically distinct type of ULIRGs. A morphological trend is seen along the sequence (1)-(2)-(3), in general agreement with the standard ULIRG - QSO evolution scenario and suggestive of a broad peak in extinction during the intermediate stages of merger evolution. However, the scatter in this sequence, implies that black hole accretion, in addition to depending on the merger phase, also has a strong chaotic/random component, as in local AGN. (abridged)Comment: 61 pages, 39 figures, 16 tables, accepted for publication in ApJS, June 2009 issue. Unabbreviated version can be found at http://www.astro.umd.edu/~veilleux/pubs/quest4.pd

Integral cross sections for electron scattering by ground state Ba atoms

We have used the convergent close-coupling method and a unitarized first-order many-body theory to calculate integral cross sections for elastic scattering and momentum transfer, for excitation of the 5d^2 ^1S, 6s6p^1P_1, 6s7p^1P_1, 6s8p^1P_1, 6s5d^1D_2, 5d^2^1D_2, 6s6d^1D_2, 6p5d^1F_3, 6s4f^1F_3, 6p5d^1D_2, 6s6p^3P_{0,1,2}, 6s5d^3D_{1,2,3}, and 6p5d^3D_2 states, for ionization and for total scattering by electron impact on the ground state of barium at incident electron energies from 1 to 1000 eV. These results and all available experimental data have been combined to produce a recommended set of integral cross sections.Comment: 47 pages, 8 tables, 25 figure