The galaxies M82, NGC 253, NGC 1068, and NGC 4945 have been detected in gamma-rays by Fermi. Previously, we developed and tested a model for cosmic ray interactions in the starburst galaxy M82. Now, we aim to explore the differences between starburst
and active galactic nuclei (AGN) environments by applying our self-consistent model to the starburst galaxy NGC 253 and the Seyfert galaxy NGC 1068. Assuming constant cosmic-ray acceleration efficiency by supernova remnants with Milky-Way parameters, we calculate the cosmic-ray proton and primary and secondary electron/positron populations, predict the radio and gamma-ray spectra, and compare with published measurements. We find that our models easily fits the observed gamma-ray spectrum for NGC 253 while constraining the cosmic ray source spectral index and acceleration efficiency. However, we encountered difficultly modeling the observed radio data and constraining the speed of the galactic wind and the magnetic field strength, unless the gas mass is less than currently preferred values. Additionally, our starburst model consistently underestimates the observed gamma-ray flux and overestimates the radio flux for NGC 1068; these issues would be resolved if the AGN is the primary source of gamma-rays. We discuss the implications of these results and make predictions for the neutrino fluxes for both galaxies.
We construct a family of models for the evolution of energetic particles in the starburst galaxy M82 and compare them to observations to test the calorimeter assumption that all cosmic ray energy is radiated in the starburst region. Assuming constant
cosmic ray acceleration efficiency with Milky Way parameters, we calculate the cosmic-ray proton and primary and secondary electron/positron populations as a function of energy. Cosmic rays are injected with Galactic energy distributions and electron-to-proton ratio via type II supernovae at the observed rate of 0.07/yr. From the cosmic ray spectra, we predict the radio synchrotron and gamma-ray spectra. To more accurately model the radio spectrum, we incorporate a multiphase interstellar medium in the starburst region of M82. Our model interstellar medium is highly fragmented with compact dense molecular clouds and dense photoionized gas, both embedded in a hot, low density medium in overall pressure equilibrium. The spectra predicted by this one-zone model are compared to the observed radio and gamma-ray spectra of M82. Chi-squared tests are used with radio and gamma-ray observations and a range of model predictions to find the best-fit parameters. The best-fit model yields constraints on key parameters in the starburst zone of M82, including a magnetic field strength of ~250 mu G and a wind advection speed in the range of 300-700 km/s. We find that M82 is a good electron calorimeter but not an ideal cosmic-ray proton calorimeter and discuss the implications of our results for the astrophysics of the far infrared-radio correlation in starburst galaxies.
We present deep WIYN H_alpha SparsePak and DensePak spatially-resolved optical spectroscopy of the dwarf irregular starburst galaxy NGC 1140. The different spatial resolutions and coverage of the two sets of observations have allowed us to investigat
e the properties and kinematics of the warm ionized gas within both the central regions of the galaxy and the inner halo. We find that the position angle of the H_alpha rotation axis for the main body of the galaxy is consistent with the HI rotation axis at PA = 39 deg, but that the ionized gas in the central 20x20 arcsecs (~2x2 kpc) is kinematically decoupled from the rest of the system, and rotates at a PA approximately perpendicular to that of the main body of the galaxy at +40 deg. We find no evidence of coherent large-scale galactic outflows. Instead multiple narrow emission line components seen within a radius of ~1-1.5 kpc, and high [SII]/H_alpha ratios found beyond ~2 kpc implying a strong contribution from shocks, suggest that the intense star formation is driving material outwards from the main star forming zone in the form of a series of interacting superbubbles/shells. A broad component (100<FWHM<230 km/s) to the H_alpha line is identified throughout galaxy disk out to >2 kpc. Based on recent work, we conclude that it is produced in turbulent mixing layers on the surfaces of cool gas knots embedded within the ISM, set up by the feedback from young massive star clusters. Our data suggest a physical limit to the radius where the broad emission line component is significant, and we propose that this limit marks a significant transition point in the development of the galactic outflow, where turbulent motion becomes less dominant. This mirrors what has recently been found in another similar irregular starburst galaxy NGC 1569.
NGC1275, the central galaxy in the Perseus cluster, is the host of gigantic hot bipolar bubbles inflated by AGN jets observed in the radio as Perseus A. It presents a spectacular $H{alpha}$-emitting nebulosity surrounding NGC1275, with loops and fila
ments of gas extending to over 50 kpc. The origin of the filaments is still unknown, but probably correlates with the mechanism responsible for the giant buoyant bubbles. We present 2.5 and 3-dimensional MHD simulations of the central region of the cluster in which turbulent energy, possibly triggered by star formation and supernovae (SNe) explosions is introduced. The simulations reveal that the turbulence injected by massive stars could be responsible for the nearly isotropic distribution of filaments and loops that drag magnetic fields upward as indicated by recent observations. Weak shell-like shock fronts propagating into the ICM with velocities of 100-500 km/s are found, also resembling the observations. The isotropic outflow momentum of the turbulence slows the infall of the intracluster medium, thus limiting further starburst activity in NGC1275. As the turbulence is subsonic over most of the simulated volume, the turbulent kinetic energy is not efficiently converted into heat and additional heating is required to suppress the cooling flow at the core of the cluster. Simulations combining the MHD turbulence with the AGN outflow can reproduce the temperature radial profile observed around NGC1275. While the AGN mechanism is the main heating source, the supernovae are crucial to isotropize the energy distribution.
(Abridged) In this second paper of the series, we present the results from optical Gemini-North GMOS-IFU and WIYN DensePak IFU spectroscopic observations of the starburst and inner wind zones of M82, with a focus on the state of the T~10^4 K ionized
interstellar medium. Our electron density maps show peaks of a few 1000 cm-3, local small spatial-scale variations, and a fall-off in the minor axis direction. We discuss the implications of these results with regards to the conditions/locations that may favour the escape of individual cluster winds. Our findings imply that the starburst environment is highly fragmented into a range of clouds from small/dense clumps with low filling factors (<1pc, n_e>10^4 cm-3) to larger filling factor, less dense gas. The near-constant state of the ionization state of the ~10^4 K gas throughout the starburst can be explained as a consequence of the small cloud sizes, which allow the gas conditions to respond quickly to any changes. We have examined in more detail both the broad (FWHM 150-350 km/s) line component found in Paper I that we associated with emission from turbulent mixing layers on the gas clouds, and the discrete outflow channel identified within the inner wind. The channel appears as a coherent, expanding cylindrical structure of length >120 pc and and width 35-50 pc and the walls maintain an approximately constant (but subsonic) expansion velocity of ~60 km/s. We use the channel to examine further the relationship between the narrow and broad component emitting gas within the inner wind. Within the starburst energy injection zone, we find that turbulent motions (as traced by the broad component) appear to play an increasing role with height.
Like other starburst galaxies, M82 hosts compact, massive young star clusters that are interesting both in their own right and as benchmarks for population synthesis models. Can spectral synthesis models at resolutions around 1000 adequately reproduc
e the near-IR spectral features and the energy distribution of these clusters between 0.8 and 2.4 microns? How do the derived cluster properties compare with previous results from optical studies? We analyse the spectra of 5 massive clusters in M82, using data acquired with the spectrograph SpeX on the InfraRed Telescope Facility (NASA/IRTF) and a new population synthesis tool with a highly improved near-IR extension, based on a recent collection of empirical and theoretical spectra of red supergiant stars. We obtain excellent fits across the near-IR with models at quasi-solar metallicity and a solar neighbourhood extinction law. Spectroscopy breaks a strong degeneracy between age and extinction in the near-IR colours in the red supergiant-dominated phase of evolution. The estimated near-IR ages cluster between 9 and 30 Myr, i.e. the ages at which the molecular bands due to luminous red supergiants are strongest in the current models. They do not always agree with optical spectroscopic ages. Adding optical data sometimes leads to the rejection of the solar neighbourhood extinction law. This is not surprising considering small-scale structure around the clusters, but it has no significant effect on the near-IR based spectroscopic ages. [abridged]
We have obtained optical spectra of 29 early-type (E/S0) galaxies that hosted type Ia supernovae (SNe Ia). We have measured absorption-line strengths and compared them to a grid of models to extract the relations between the supernova properties and
the luminosity-weighted age/composition of the host galaxies. The same analysis was applied to a large number of early-type field galaxies selected from the SDSS spectroscopic survey. We find no difference in the age and abundance distributions between the field galaxies and the SN Ia host galaxies. We do find a strong correlation suggesting that SNe Ia in galaxies whose populations have a characteristic age greater than 5 Gyr are ~ 1 mag fainter at V(max) than those found in galaxies with younger populations. However, the data cannot discriminate between a smooth relation connecting age and supernova luminosity or two populations of SN Ia progenitors. We find that SN Ia distance residuals in the Hubble diagram are correlated with host-galaxy metal abundance, consistent with the predictions of Timmes, Brown & Truran (2003). The data show that high iron abundance galaxies host less-luminous supernovae. We thus conclude that the time since progenitor formation primarily determines the radioactive Ni production while progenitor metal abundance has a weaker influence on peak luminosity, but one not fully corrected by light-curve shape and color fitters. Assuming no selection effects in discovering SNe Ia in local early-type galaxies, we find a higher specific SN Ia rate in E/S0 galaxies with ages below 3 Gyr than in older hosts. The higher rate and brighter luminosities seen in the youngest E/S0 hosts may be a result of recent star formation and represents a tail of the prompt SN Ia progenitors.
