VV 655, a dwarf irregular galaxy with HI tidal debris, is a companion to the lenticular luminous infrared galaxy (LIRG) NGC 4418. NGC 4418 stands out among nearby LIRGs due to its dense central concentration of molecular gas and the dusty, bi-polar structures along its minor axis suggestive of a wind driven by a central starburst and possible nuclear activity. We seek to understand the consequences of the ongoing minor interaction between VV 655 and NGC 4418 for the evolution of the LIRG, including the origin of the gas supply responsible for its unusual nuclear properties. We investigate the structural, kinematic, and chemical properties of VV 655 and NGC 4418 by analyzing archival imaging data and optical spectroscopic observations from the SDSS-III and new spectra from SALT-RSS. We characterize their gas-phase metal abundances and spatially resolved, ionized gas kinematics, and demonstrate that the gas-phase metallicity in NGC 4418 significantly exceeds that in VV 655. No kinematic disturbances in the ionized gas are observed along the minor axis of NGC 4418, but we see evidence for ionized gas outflows from VV 655 that may increase the cross-section for gas stripping in grazing collisions. A faint, asymmetric outer arm is detected in NGC 4418 of the type normally associated with galaxy-galaxy interactions. The simplest model suggests that the minor interaction between VV 655 and NGC 4418 produced the unusual nuclear properties of the LIRG via tidal torquing of the interstellar medium of NGC 4418 rather than through a significant gas transfer event. In addition to inducing a central concentration of gas in NGC 4418, this interaction also produced an enhanced star formation rate and an outer tidal arm in the LIRG. The VV 655-NGC 4418 system offers an example of the potential for minor collisions to alter the evolutionary pathways of giant galaxies.

Recent analyses of the gamma-ray spectrum from the ultra-luminous infrared galaxy Arp 220 have revealed a discrepancy in the cosmic ray energy injection rates derived from the gamma-rays versus the radio emission. While the observed radio emission is consistent with the star formation rate inferred from infrared observations, a significantly higher cosmic ray population is necessary to accurately model the measured gamma-ray flux. To resolve this discrepancy between the radio and gamma-ray observations, we find that we must increase the cosmic ray energy injection rate and account for an infrared optical depth greater than unity. Raising the energy injection rate naturally raises the total gamma-ray flux but also raises the radio flux unless there is also an increase in the energy loss rate for cosmic ray leptons. A optically thick medium results in an increase in energy losses via inverse Compton for cosmic ray leptons and preserves agreement with submillimeter, millimeter, and infrared wavelength observations.

The Cygnus X region contains giant molecular cloud complexes and populous associates of massive young stars. The discovery of spatially extended, hard gamma-ray emission in Cygnus X by both Milagro and Fermi indicates that Cygnus X is also a potential source of high-energy Galactic neutrinos. Here, we adapt our single-zone model for cosmic ray interactions in the central molecular zones of starburst galaxies for use in Cygnus X. We calculate the potential neutrino flux corresponding to the hard gamma-ray emission from the "Cygnus Cocoon" and to the soft, diffuse interstellar gamma-ray emission. We check our results by comparing the corresponding gamma-ray emission against the Fermi interstellar emission model and Milagro, ARGO-YBJ, and HAWC observations. In comparing our results against a recent IceCube analysis and the current sensitivity limits, we find that neutrino emission from the Cocoon has a large enough flux that it could plausibly be detected, provided hadronic interactions are occurring at sufficiently high energies. High-energy neutrinos from Cygnus X would provide direct evidence for the presence of as yet unidentified PeV energy accelerators in the Galactic disk.

Extragalactic cosmic ray populations are important diagnostic tools for tracking the distribution of energy in nuclei and for distinguishing between activity powered by star formation versus active galactic nuclei (AGNs). Here, we compare different diagnostics of the cosmic ray populations of the nuclei of Arp 220 based on radio synchrotron observations and the recent gamma-ray detection. We find the gamma-ray and radio emission to be incompatible; a joint solution requires at minimum a factor of 4 - 8 times more energy coming from supernovae and a factor of 40 - 70 more mass in molecular gas than is observed. We conclude that this excess of gamma-ray flux in comparison to all other diagnostics of star-forming activity indicates that there is an AGN present that is providing the extra cosmic rays, likely in the western nucleus.

We analyse new observations with the International Low Frequency Array (LOFAR) telescope, and archival data from the Multi-Element Radio Linked Interferometer Network (MERLIN) and the Karl G. Jansky Very Large Array (VLA). We model the spatially resolved radio spectrum of Arp 220 from 150 MHz to 33 GHz. We present an image of Arp 220 at 150 MHz with resolution $0.65''\times0.35''$, sensitivity 0.15 mJy beam$^{-1}$, and integrated flux density $394\pm59$ mJy. More than 80% of the detected flux comes from extended ($6''\approx$2.2 kpc) steep spectrum ($\alpha=-0.7$) emission, likely from star formation in the molecular disk surrounding the two nuclei. We find elongated features extending $0.3''$ (110 pc) and $0.9''$ (330 pc) from the eastern and western nucleus respectively, which we interpret as evidence for outflows. The extent of radio emission requires acceleration of cosmic rays far outside the nuclei. We find that a simple three component model can explain most of the observed radio spectrum of the galaxy. When accounting for absorption at 1.4 GHz, Arp 220 follows the FIR/radio correlation with $q=2.36$, and we estimate a star formation rate of 220 M$_\odot\text{yr}^{-1}$. We derive thermal fractions at 1 GHz of less than 1% for the nuclei, which indicates that a major part of the UV-photons are absorbed by dust. International LOFAR observations shows great promise to detect steep spectrum outflows and probe regions of thermal absorption. However, in LIRGs the emission detected at 150 MHz does not necessarily come from the main regions of star formation. This implies that high spatial resolution is crucial for accurate estimates of star formation rates for such galaxies at 150 MHz.

The energy densities in magnetic fields and cosmic rays (CRs) in galaxies are often assumed to be in equipartition, allowing for an indirect estimate of the magnetic field strength from the observed radio synchrotron spectrum. However, both primary and secondary CRs contribute to the synchrotron spectrum, and the CR electrons also loose energy via bremsstrahlung and inverse Compton. While classical equipartition formulae avoid these intricacies, there have been recent revisions that account for the extreme conditions in starbursts. Yet, the application of the equipartition formula to starburst environments also presupposes that timescales are long enough to reach equilibrium. Here, we test equipartition in the central molecular zones (CMZs) of nearby starburst galaxies by modeling the observed gamma-ray spectra, which provide a direct measure of the CR energy density, and the radio spectra, which provide a probe of the magnetic field strength. We find that in starbursts, the magnetic field energy density is significantly larger than the CR energy density, demonstrating that the equipartition argument is frequently invalid for CMZs.

The cores of Arp 220, the closest ultra-luminous infrared starburst galaxy, provide an opportunity to study interactions of cosmic rays under extreme conditions. In this paper, we model the populations of cosmic rays produced by supernovae in the central molecular zones of both starburst nuclei. We find that ~65 - 100% of cosmic rays are absorbed in these regions due to their huge molecular gas contents, and thus, the nuclei of Arp 220 nearly complete proton calorimeters. As the cosmic ray protons collide with the interstellar medium, they produce secondary electrons that are also contained within the system and radiate synchrotron emission. Using results from chi-squared tests between the model and the observed radio spectral energy distribution, we predict the emergent gamma-ray and high-energy neutrino spectra and find the magnetic field to be at milligauss levels. Because of the extremely intense far-infrared radiation fields, the gamma-ray spectrum steepens significantly at TeV energies due to gamma-gamma absorption.

The conditions in the Galactic Center are often compared with those in starburst systems, which contain higher supernova rates, stronger magnetic fields, more intense radiation fields, and larger amounts of dense molecular gas than in our own Galactic disk. Interactions between such an augmented interstellar medium and cosmic rays result in brighter radio and gamma-ray emission. Here, we test how well the comparisons between the Galactic Center and starburst galaxies hold by applying a model for cosmic ray interactions to the Galactic Center to predict the resulting gamma-ray emission. The model only partially explains the observed gamma-ray and radio emission. The model for the gamma-ray spectrum agrees with the data at TeV energies but not at GeV energies. Additionally, as the fits of the model to the radio and gamma-ray spectra require significant differences in the optimal wind speed and magnetic field strength, we find that the single-zone model alone cannot account for the observed emission from the Galactic Center. Our model is improved by including a soft, additional cosmic-ray population. We assess such a cosmic ray population and its potential sources and find that a cosmic-ray electron spectrum is energetically favored over a cosmic-ray proton spectrum.

The Milky Way appears is a typical barred spiral, and comparisons can be made between its nuclear region and those of structurally similar nearby spirals. Maffei 2, M83, IC 342 and NGC 253 are nearby systems whose nuclear region properties contrast with those of the Milky Way. Stellar masses derived from NIR photometery, molecular gas masses and star formation rates allow us to assess the evolutionary states of this set of nuclear regions. These data suggest similarities between nuclear regions in terms of their stellar content while highlighting significant differences in current star formation rates. In particular current star formation rates appear to cover a larger range than expected based on the molecular gas masses. This behavior is consistent with nuclear region star formation experiencing episodic variations. Under this hypothesis the Milky Way's nuclear region currently may be in a low star formation rate phase.

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.

How cosmic rays sample the multi-phase interstellar medium (ISM) in starburst galaxies has important implications for many science goals, including evaluating the cosmic ray calorimeter model for these systems, predicting their neutrino fluxes, and modeling their winds. Here, we use Monte Carlo simulations to study cosmic ray sampling of a simple, two-phase ISM under conditions similar to those of the prototypical starburst galaxy M82. The assumption that cosmic rays sample the mean density of the ISM in the starburst region is assessed over a multi-dimensional parameter space where we vary the number of molecular clouds, the galactic wind speed, the extent to which the magnetic field is tangled, and the cosmic ray injection mechanism. We evaluate the ratio of the emissivity from pion production in molecular clouds to the emissivity that would be observed if the cosmic rays sampled the mean density, and seek areas of parameter space where this ratio differs significantly from unity. The assumption that cosmic rays sample the mean density holds over much of parameter space; however, this assumption begins to break down for high cloud density, injection close to the clouds, and a very tangled magnetic field. We conclude by evaluating the extent to which our simulated starburst region behaves as a proton calorimeter and constructing the time-dependent spectrum of a burst of cosmic rays.

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.