We present results from a survey of molecular hydrogen emission from a sample of Starburst and Seyfert galaxies carried out with the Infrared Space Observatory (ISO). Pure rotational H_2 emission has been detected in a number of extragalactic objects and a variety of environments. A number of transitions from S(7) to S(0) are detected in both Starbursts and Seyferts. Using excitation diagrams we derive temperatures and masses of the ``warm'' molecular hydrogen. We find that the temperature of the ``warm'' gas is similar in Starbursts and Seyferts (those Seyferts for which we have firm detections of the S(0) line) with a value of around T~150 K. This ``warm'' gas accounts for as much as 10% of the total galactic mass (as probed by CO molecular observations) in Starbursts. The fraction of ``warm'' gas is overall higher in Seyferts, ranging between 2--35%. We then investigate the origin of the warm H_2 emission. Comparison with published theoretical models and Galactic templates implies that although emission from photodissociation regions (PDR) alone could explain the emission from Starbursts and Seyferts, most likely a combination of PDR, shock emission and gas heated by X-rays (mostly for the Seyferts) is responsible for H_2 excitation in extragalactic environments. Finally, we find that although PAH and H_2 line emission correlate well in Starbursts and the large scale emission in AGN, H_2 emission is much stronger compared to PAH emission in cases where a ``pure'' AGN dominates the energy output.

We present new observations of the central regions of the starburst galaxy M82 consisting of near-infrared H- and K-band integral field spectroscopy from the MPE 3D instrument, and 2.4-45 micron spectroscopy from the Short Wavelength Spectrometer (SWS) on board ISO. These measurements are used, together with data from the literature, to re-examine the issue of extinction, determine the physical conditions of the interstellar medium, and characterize the composition of the stellar populations. We find that purely foreground extinction cannot reproduce the global relative intensities of H recombination lines from optical to radio wavelengths. A good fit is provided by a homogeneous mixture of dust and sources, and with a visual extinction of 52 mag. The SWS data provide evidence for deviations from commonly assumed extinction laws between 3 and 10 micron. The fine-structure lines of Ne, Ar, and S detected with SWS imply an electron density of about 300 cm^{-3}, and abundance ratios Ne/H and Ar/H nearly solar and S/H one-fourth solar. The excitation of the ionized gas indicates an average temperature for the OB stars of 37400 K, with little spatial variation. A random distribution of gas clouds and ionizing clusters and an ionization parameter of about 10^{-2.3} represent well the star-forming regions. From population synthesis and the mass-to-K-light ratio, we conclude that the near-infrared continuum emission across the starburst regions is dominated by red supergiants with average temperatures ranging from 3600 K to 4500 K, and roughly solar metallicity. Our data rule out significant contributions from older metal-rich giants in the central few tens of parsecs of M82.

We present new ISO-SWS data for a sample of 27 starburst galaxies, and with these data examine the issues of formation and evolution of the most massive stars in starburst galaxies. Using starburst models which incorporate time evolution, new stellar atmosphere models for massive stars, and a starburst model geometry derived from observations of the prototypical starburst M82, we model the integrated mid-infrared line ratio [NeIII](15.6 microns)/[NeII](12.8 microns). This line ratio is sensitive to the hardness of the stellar energy distribution and therefore to the most massive stars present. We conclude from our models, with consideration of recent determinations of the stellar census in local, high-mass star forming regions, that the [NeIII]/[NeII] ratios we measure are consistent with the formation of massive (~50-100 solar mass) stars in most starbursts. In this framework, the low nebular excitation inferred from the measured line ratios can be attributed to aging effects. By including estimates of the ratio of infrared-to-Lyman continuum luminosity for the galaxies in our sample, we further find that most starbursts are relatively short-lived (1-10 million years), only a few O-star lifetimes. We discuss a possible cause of such short events: the effectiveness of stellar winds and supernovae in destroying the starburst environment.

We present ISO observations of the S(0), S(1), S(2), and S(3) rotational lines of molecular hydrogen from two molecular clouds near the Galactic Center (GC). We have also measured continuum dust emission at infrared wavelengths with ISO and the rotational radio lines J=1-0 of 13CO and C18O and J=2-1 of C18O with the IRAM-30m telescope. Using the dust continuum spectra and the CO lines we derive a total visual extinction of 15-20 magnitudes toward these GC clouds. After correcting the H2 data for extinction, the gas temperatures are \~250 K and the column densities of warm gas are ~2e21 cm-2. This is the first direct measure of the H2 column densities of the warm component; with this, we estimate an NH3 abundance in the warm gas of ~ 2e-7. The column density of warm gas is, at least, a factor of 100 larger than the corresponding column densities derived from the warm dust. The observed ortho-to-para ratio (OTPR) is ~ 1, clearly below the local thermodynamical equilibrium (LTE) OTPR for gas at 250 K of ~3. Low velocity shocks (~ 10 km/s) are the most likely explanation for the column densities of warm gas and dust and the non-LTE H2 OTPR.

We present an inventory of mid-infrared spectral features detected in high resolution 2.4 to 45um ISO-SWS spectra of the galaxies M82, NGC253, Circinus, NGC1068, and the 30Doradus region of the LMC. We discuss their identifications and highlight possible relations between these features and the physical state of the interstellar medium in galaxies. In the 5--13 um region, emission from unidentified infrared bands (UIBs), usually ascribed to aromatic molecules, and apparent silicate absorption dominate the spectrum. We demonstrate the danger of overestimating the amount of silicate absorption by comparing the spectrum of M82 to the (absorption free) spectrum of the reflection nebula NGC7023. No strong silicate absorption is present in M82. The (very small grain) dust continuum under the UIB emission in our starburst templates can be modeled by a simple power law, starting at wavelengths between 8 and 9um. The spectra are valuable templates for future mid-infrared missions. We smooth our data to simulate low resolution spectra as obtained with ISOCAM-CVF, ISOPHOT-S, and in the future with the low resolution mode of SIRTF-IRS, and use our high spectral resolution information to highlight possible identification problems at low resolving power that are caused by coincidences of lines and features.

We report on observations of SN 1987A with ISO SWS/LWS made 9 - 11 years after the explosion. No emission from the supernova was seen. In particular, the upper limits on the fluxes of [Fe I] 24.05 microns and [Fe II] 25.99 microns on day 3999 are roughly 1.1 Jy and 1.4 Jy, respectively. Assuming a homogeneous distribution of Ti-44 inside 2000 km/s, we have made theoretical models to estimate the mass of ejected Ti-44. Assessing various uncertainties of the model, we obtain an upper limit of 1.5\EE{-4} solar masses. The implications of this are discussed. The LWS data display continuum emission as well as nebular lines of [O I], [C II] and [O III] from neighboring photoexcited regions in the LMC. The [O III] lines indicate an electron density of 120 +/- 75 cm^{-3}, and the continuum can be explained by dust with a temperature of roughly 37 K. A second dust component with approximately 10 K may also be present.

We have observed the S(0), S(1), S(2), S(3), S(4) and S(5) rotational lines of molecular hydrogen (H2) towards the peak of the photodissociation region (PDR) associated with the reflection nebula NGC 7023. The observed H2 line ratios show that they arise in warm gas with kinetic temperatures ~300 - 700 K. However, the data cannot be fitted by an ortho- to para- (OTP) ratio of 3. An OTP ratio in the range ~1.5 - 2 is necessary to explain our observations. This is the first detection of a non-equilibrium OTP ratio measured from the H2 pure-rotational lines in a PDR. The existence of a dynamical PDR is discussed as the most likely explanation for this low OTP ratio.

We report the detection of faint emission in the high-excitation [OIV] 25.90um line in a number of starburst galaxies, from observations obtained with the Short Wavelength Spectrometer (SWS) on board ISO. Further observations of M82 spatially resolve the [OIV] emitting region. Detection of this line in starbursts is surprising since it is not produced in measurable quantities in HII regions around hot main-sequence stars, the dominant energy source of starburst galaxies. We discuss various models for the formation of this line. [OIV] that is spatially resolved by ISO cannot originate in a weak AGN and must be due to very hot stars or ionizing shocks related to the starburst activity. For low-excitation starbursts like M82, shocks are the most plausible source of [OIV] emission.

We present an ISO SWS and ISOPHOT-S, mid-infrared spectroscopic survey of 15 ultra-luminous IRAS galaxies. We combine the survey results with a detailed case study, based on near-IR and mm imaging spectroscopy, of one of the sample galaxies (UGC 5101). We compare the near- and mid-IR characteristics of these ultra-luminous galaxies to ISO and literature data of thirty starburst and active galactic nuclei (AGN), template galaxies. We find that 1) 70-80% of the ultra-luminous IRAS galaxies in our sample are predominantly powered by recently formed massive stars. 20-30% are powered by a central AGN. These conclusions are based on a new infrared 'diagnostic diagram' involving the ratio of high to low excitation mid-IR emission lines on the one hand, and on the strength of the 7.7um PAH feature on the other hand. 2) at least half of the sources probably have simultaneously an active nucleus and starburst activity in a 1-2 kpc diameter circum-nuclear disk/ring. 3) the mid-infrared emitting regions are highly obscured. After correction for these extinctions, we estimate that the star forming regions in ULIRGs have ages between 10^7 and 10^8 years, similar to but somewhat larger than those found in lower luminosity starburst galaxies. 4) in the sample we have studied there is no obvious trend for the AGN component to dominate in the most compact, and thus most advanced mergers. Instead, at any given time during the merger evolution, the time dependent compression of the circum-nuclear interstellar gas, the accretion rate onto the central black hole and the associated radiation efficiency may determine whether star formation or AGN activity dominates the luminosity of the system.