The sources of ultrahigh-energy cosmic rays (UHECRs) have been difficult to catch. It was recently pointed out that while sources of UHECR protons exhibit anisotropy patterns that become denser and compressed with rising energy, nucleus-emitting-sources give rise to a cepa stratis (onion-like) structure with layers that become more distant from the source position with rising energy. The peculiar shape of the hot spots from nucleus-accelerators is steered by the competition between energy loss during propagation and deflection on the Galactic magnetic field (GMF). Here, we run a full-blown simulation study to accurately characterize the deflections of UHECR nuclei in the GMF. We show that while the cepa stratis structure provides a global description of anisotropy patterns produced by UHECR nuclei en route to Earth, the hot spots are elongated depending on their location in the sky due to the regular structure of the GMF. We demonstrate that with a high-statistics sample at the high-energy-end of the spectrum, like the one to be collected by NASA's POEMMA mission, the energy dependence of the hot-spot contours could become a useful observable to identify the nuclear composition of UHECRs. This new method to determine the nature of the particle species is complementary to those using observables of extensive air showers, and therefore is unaffected by the large systematic uncertainties of hadronic interaction models.

Diffuse gamma-ray emission from interstellar clouds results largely from cosmic ray (CR) proton collisions with ambient gas, regardless of the gas state, temperature, or dust properties of the cloud. The interstellar medium is predominantly transparent to both CRs and gamma-rays, so GeV emission is a unique probe of the total gas column density. The gamma-ray emissivity of a cloud of known column density is then a measure of the impinging CR population and may be used to map the kpc-scale CR distribution in the Galaxy. To this end, we test a number of commonly used column density tracers to evaluate their effectiveness in modeling the GeV emission from the relatively quiescent, nearby $\rho$ Ophiuchi molecular cloud. We confirm that both \hi\ and an appropriate H$_2$ tracer are required to reproduce the total gas column densities probed by diffuse gamma-ray emisison. We find that the optical depth at 353 GHz $\tau_{353}$ from Planck reproduces the gamma-ray data best overall based on the test statistic across the entire region of interest, but near infrared stellar extinction also performs very well, with smaller spatial residuals in the densest parts of the cloud.

We analyze the gamma-ray emission from 9 high latitude, translucent molecular clouds taken with the Fermi Large Area Telescope (LAT) between 250 MeV and 10 GeV. Observations of gamma-rays allow us to probe the density and spectrum of cosmic rays in the solar neighborhood. The clouds studied lie within $\sim\!$270 pc from the Sun and are selected from the Planck all-sky CO map. Gamma-rays in this energy range mostly result from cosmic ray interactions with the interstellar medium, which is traced with three components: HI, CO, and dark gas. Every cloud is detected and shows significant, extended gamma-ray emission from molecular gas. The gamma-ray emission is dominated by the CO-emitting gas in some clouds, but by the CO-dark gas in others. The average emissivity and gamma-ray power law index from HI above 1 GeV shows no evidence of a systematic variation. The CO-to-H$_2$ conversion factor shows no variation between clouds over this small spatial range, but shows significant variations within each cloud. The average CO-to-H$_2$ conversion factor suggests that the CO-dark gas is molecular as opposed to optically thick HI.

The physical relationship between the far-infrared and radio fluxes of star forming galaxies has yet to be definitively determined. The favored interpretation, the "calorimeter model," requires that supernova generated cosmic ray (CR) electrons cool rapidly via synchrotron radiation. However, this cooling should steepen their radio spectra beyond what is observed, and so enhanced ionization losses at low energies from high gas densities are also required. Further, evaluating the minimum energy magnetic field strength with the traditional scaling of the synchrotron flux may underestimate the true value in massive starbursts if their magnetic energy density is comparable to the hydrostatic pressure of their disks. Gamma-ray spectra of starburst galaxies, combined with radio data, provide a less ambiguous estimate of these physical properties in starburst nuclei. While the radio flux is most sensitive to the magnetic field, the GeV gamma-ray spectrum normalization depends primarily on gas density. To this end, spectra above 100 MeV were constructed for two nearby starburst galaxies, NGC 253 and M82, using Fermi data. Their nuclear radio and far-infrared spectra from the literature are compared to new models of the steady-state CR distributions expected from starburst galaxies. Models with high magnetic fields, favoring galaxy calorimetry, are overall better fits to the observations. These solutions also imply relatively high densities and CR ionization rates, consistent with molecular cloud studies.

We explore the cosmic evolution of radio AGN with low radio powers (L_1.4GHz < 5\times10^25 W/Hz) out to z=1.3 using to-date the largest sample of ~600 low luminosity radio AGN at intermediate redshift drawn from the VLA-COSMOS survey. We derive the radio luminosity function for these AGN, and its evolution with cosmic time assuming two extreme cases: i) pure luminosity and ii) pure density evolution. The former and latter yield L_*\propto(1+z)^(0.8+/-0.1), and Phi*\propto (1+z)^(1.1+/-0.1), respectively, both implying a fairly modest change in properties of low radio-power AGN since z=1.3. We show that this is in stark contrast with the evolution of powerful (L_1.4GHz > 5\times10^25 W/Hz) radio AGN over the same cosmic time interval, constrained using the 3CRR, 6CE, and 7CRS radio surveys by Willott et al. (2001). We demonstrate that this can be explained through differences in black hole fueling and triggering mechanisms, and a dichotomy in host galaxy properties of weak and powerful AGN. Our findings suggest that high and low radio-power AGN activity is triggered in different stages during the formation of massive red galaxies. We show that weak radio AGN occur in the most massive galaxies already at z~1, and they may significantly contribute to the heating of their surrounding medium and thus inhibit gas accretion onto their host galaxies, as recently suggested for the `radio mode' in cosmological models.

We derive the cosmic star formation history (CSFH) out to z=1.3 using a sample of ~350 radio-selected star-forming galaxies, a far larger sample than in previous, similar studies. We attempt to differentiate between radio emission from AGN and star-forming galaxies, and determine an evolving 1.4 GHz luminosity function based on these VLA-COSMOS star forming galaxies. We precisely measure the high-luminosity end of the star forming galaxy luminosity function (SFR>100 M_Sol/yr; equivalent to ULIRGs) out to z=1.3, finding a somewhat slower evolution than previously derived from mid-infrared data. We find that more stars are forming in luminous starbursts at high redshift. We use extrapolations based on the local radio galaxy luminosity function; assuming pure luminosity evolution, we derive $L_* \propto (1+z)^{2.1 \pm 0.2}$ or $L_* \propto (1+z)^{2.5 \pm 0.1}$, depending on the choice of the local radio galaxy luminosity function. Thus, our radio-derived results independently confirm the ~1 order of magnitude decline in the CSFH since z~1.

We examine the faint-end slope of the rest-frame V-band luminosity function (LF), with respect to galaxy spectral type, of field galaxies with redshift z<0.5, using a sample of 80,820 galaxies with photometric redshifts in the Cosmic Evolution Survey (COSMOS) field. For all galaxy spectral types combined, the LF slope, alpha, ranges from -1.24 to -1.12, from the lowest redshift bin to the highest. In the lowest redshift bin (0.02<z<0.1), where the magnitude limit is M(V) ~ -13, the slope ranges from ~ -1.1 for galaxies with early-type spectral energy distributions (SEDs), to ~ -1.9 for galaxies with low-extinction starburst SEDs. In each galaxy SED category (Ell, Sbc, Scd/Irr, and starburst), the faint-end slopes grow shallower with increasing redshift; in the highest redshift bin (0.4<z<0.5), the slope is ~ -0.5 and ~ -1.3 for early-types and starbursts respectively. The steepness of alpha at lower redshift could be qualitatively explained by large numbers of faint dwarf galaxies, perhaps of low surface brightness, which are not detected at higher redshifts.

We present 5.2" x 2.6" resolution interferometry of CO J=1-0 emission from the starburst galaxy NGC 253. The high spatial resolution of these new data, in combination with recent high resolution maps of 13CO, HCN and near-infrared emission, allow us for the first time to link unambiguously the gas properties in the central starburst of NGC 253 with its bar dynamics. We confirm that the star formation results from bar-driven gas flows as seen in "twin peaks" galaxies. Two distinct kinematic features are evident from the CO map and position-velocity diagram: a group of clouds rotating as a solid body about the kinematic center of the galaxy, and a more extended gas component associated with the near-infrared bar. We model the line intensities of CO, HCN and 13CO to infer the physical conditions of the gas in the nucleus of NGC 253. The results indicate increased volume densities around the radio nucleus in a twin-peaks morphology. Compared with the CO kinematics, the gas densities appear highest near the radius of a likely inner Linblad resonance, and slightly lead the bar minor axis. This result is similar to observations of the face-on, twin-peaks galaxy NGC 6951, and is consistent with models of starburst generation due to gas inflow along a bar.

We present spectra of the extended 12CO and 13CO J=1-0 emission along the major axes of 17 nearby galaxies. Spatial variations in the ratio of CO and 13CO integrated intensities, R, are found in nearly every galaxy observed. There is an overall variation in R of 20-40% from the inner 2 kpc to the disk. Roughly a third of the survey galaxies have such gradients in R detected above the 2-sigma confidence level. Though some galaxies show a lower central value of R, on average R inside 2 kpc is 10-30% higher than R outside of 2 kpc. An increase in the CO/13CO intensity ratio from disk to nucleus may imply that the conversion factor between CO intensity and H2 column density, X, is lower in galactic nuclei. Also variations in physical conditions, most notably the gas kinetic temperature, affect both R and X. Abundance variations probably do not cause the gradient in R, though we do not rule out a decrease in effective cloud column densities in galactic nuclei possibly caused by destructive starburst superwinds. A modest rise in temperature (less than a factor of 2 or 3) from outside a 2 kpc radius towards the nucleus can easily account for the observed gradient. These results support previous work implying that X is lower in the center of the Milky Way and probably most galactic nuclei. Therefore calculating H2 masses using the standard Galactic X-factor, especially within the central few kpc of galaxies, overestimates the true mass by factors of a few. The standard X-factor still appears to be appropriate for galactic disks.

We present a search for high energy gamma-ray emission from 9 nearby starburst galaxies and M31 with the EGRET instrument aboard CGRO. Though the diffuse gamma-ray emission from starburst galaxies was suspected to be detectable, we find no emission from NGC 253, M82 nor from the average of all 9 galaxies. The 2 sigma upper limit for the EGRET flux above 100 MeV for the averaged survey observations is 1.8 x 10-8 ph cm-2 s-1. From a model of the expected radio and gamma-ray emission, we find that the magnetic field in the nuclei of these galaxies is > 25 micro Gauss, and the ratio of proton and electron densities is < 400. The EGRET limits indicate that the rate of massive star formation in the survey galaxies is only about an order of magnitude higher than in the Milky Way. The upper limit to the gamma-ray flux above 100 MeV for M31 is 1.6 x 10-8 ph cm-2 s-1. At the distance of M31, the Milky Way flux would be over twice this value, indicating higher gamma-ray emissivities in our Galaxy. Therefore, since the supernova rate of the Milky Way is higher than in M31, our null detection of M31 supports the theory of the supernova origin of cosmic rays in galaxies.