The IceCube neutrino discovery was punctuated by three showers with $E_\nu$ ~ 1-2 PeV. Interest is intense in possible fluxes at higher energies, though a marked deficit of $E_\nu$ ~ 6 PeV Glashow resonance events implies a spectrum that is soft and/or cutoff below ~few PeV. However, IceCube recently reported a through-going track event depositing 2.6 $\pm$ 0.3 PeV. A muon depositing so much energy can imply $E_{\nu_\mu} \gtrsim$ 10 PeV. We show that extending the soft $E_\nu^{-2.6}$ spectral fit from TeV-PeV data is unlikely to yield such an event. Alternatively, a tau can deposit this much energy, though requiring $E_{\nu_\tau}$ ~10x higher. We find that either scenario hints at a new flux, with the hierarchy of $\nu_\mu$ and $\nu_\tau$ energies suggesting a window into astrophysical neutrinos at $E_\nu$ ~ 100 PeV if a tau. We address implications, including for ultrahigh-energy cosmic-ray and neutrino origins.

Any interpretation of the astrophysical neutrinos discovered by IceCube must accommodate a variety of multimessenger constraints. We address implications of these neutrinos being produced in transient sources, principally if buried within supernovae so that gamma rays are absorbed by the star. This would alleviate tension with the isotropic Fermi GeV background that >10 TeV neutrinos rival in detected energy flux. We find that IceCube data constrain transient properties, implying buried GeV-TeV electromagnetic emission near or exceeding canonical SN explosion energies of ~10^51 erg, indicative of an origin within superluminous SNe. TeV neutrino bursts with dozens of IceCube events -- which would be of great use for understanding r-process nucleosynthesis and more -- may be just around the corner if they are a primary component of the flux.

The complex interplay of processes at the Galactic Center is at the heart of numerous past, present, and (likely) future mysteries. We aim at a more complete understanding of how spectra extending to >10 TeV result. We first construct a simplified model to account for the peculiar energy and angular dependence of the intense central parsec photon field. This allows for calculating anisotropic inverse Compton scattering and mapping gamma-ray extinction due to gamma gamma -> e^+ e^- attenuation. Coupling these with a method for evolving electron spectra, we examine several clear and present excesses, including the diffuse hard X-rays seen by NuSTAR and GeV gamma rays by Fermi. We address further applications to cosmic rays, dark matter, neutrinos, and gamma rays from the Center and beyond.

The Galactic Center (GC) has been long known to host gamma-ray emission detected to >10 TeV. HESS data now points to two plausible origins: the supermassive black hole (perhaps with >PeV cosmic rays and neutrinos) or high-energy electrons from the putative X-ray pulsar wind nebula G359.95-0.04 observed by Chandra and NuSTAR. We show that if the magnetic field experienced by PWN electrons is near the several mG ambient field strength suggested by radio observations of the nearby GC magnetar SGR J1745-29, synchrotron losses constrain the TeV gamma-ray output to be far below the data. Accounting for the peculiar geometry of GC infrared emission, we also find that the requisite TeV flux could be reached if the PWN is ~1 pc from Sgr A* and the magnetic field is two orders of magnitude weaker, a scenario that we discuss in relation to recent data and theoretical developments. Otherwise, Sgr A* is left, which would then be a PeV link to other AGN.

Gamma-ray observations have shown pulsars to be efficient converters of rotational energy into GeV photons and it is of wide-ranging interest to determine their contribution to the gamma-ray background. We arrive at flux predictions from both the young (<~ Myr) and millisecond (~Gyr) Galactic pulsar populations. We find that unresolved pulsars can yield both a significant fraction of the excess GeV gamma rays near the Galactic Center and an inverse Compton flux in the inner kpc similar to that inferred by Fermi. We compare models of the young pulsar population and millisecond pulsar population to constraints from gamma-ray and radio observations. Overall, we find that the young pulsars should outnumber millisecond pulsars as unassociated gamma-ray point sources in this region. The number of young radio pulsars discovered near the Galactic Center is in agreement with our model of the young pulsar population. Deeper radio observations at higher latitudes can constrain the total gamma-ray emission from both young and millisecond pulsars from the inner galaxy. While this is a step towards better understanding of pulsars, cosmic rays in the Milky Way, and searches for dark matter, we also discuss a few interesting puzzles that arise from the underlying physics of pulsar emission and evolution.

Calculations of the cosmic rate of core collapses, and the associated neutrino flux, commonly assume that a fixed fraction of massive stars collapse to black holes. We argue that recent results suggest that this fraction instead increases with redshift. With relatively more stars vanishing as "unnovae" in the distant universe, the detectability of the cosmic MeV neutrino background is improved due to their hotter neutrino spectrum, and expectations for supernova surveys are reduced. We conclude that neutrino detectors, after the flux from normal SNe is isolated via either improved modeling or the next Galactic SN, can probe the conditions and history of black hole formation.

The IceCube neutrino discovery presents an opportunity to answer long-standing questions in high-energy astrophysics. For their own sake and relations to other processes, it is important to understand neutrinos arising from the Milky Way, which should have an accompanying flux of gamma rays. Examining Fermi TeV data, and applying other constraints up to >1 PeV, it appears implausible that the Galactic fraction of the IceCube flux is large, though could be present at some level. We address Sgr A*, where the TeV-PeV neutrinos may outrun gamma rays due to gamma-gamma opacity, and further implications, including dark matter and cosmic-ray electrons.

The steep spectrum of neutrinos measured by IceCube extending from >1 PeV down to ~10 TeV has an energy flux now encroaching on the Fermi isotropic GeV background. We examine several implications starting from source energetics requirements for neutrino production. We show how the environment of extragalactic nuclei can extinguish ~10-100 TeV gamma rays and convert their energy to X-rays for plausible conditions of infrared luminosity and magnetic field, so that the Fermi background is not overwhelmed by cascades. We address a variety of scenarios, such as for acceleration by supermassive black holes and hadronic scenarios, and observations that may help elucidate the neutrinos' shadowy origins.

Studies of Fermi data indicate an excess of GeV gamma rays around the Galactic center (GC), possibly due to dark matter. We show that young gamma-ray pulsars can yield a similar signal. First, a high concentration of GC supernovae naturally leads to a population of kicked pulsars symmetric about the GC. Second, while very-young pulsars with soft spectra reside near the Galactic plane, pulsars with spectra that have hardened with age accumulate at larger angles. This combination, including unresolved foreground pulsars, traces the morphology and spectrum of the Excess.

The IceCube experiment recently detected the first flux of high-energy neutrinos in excess of atmospheric backgrounds. We examine whether these neutrinos originate from within the same extragalactic sources as ultrahigh-energy cosmic rays. Starting from rather general assumptions about spectra and flavors, we find that producing a neutrino flux at the requisite level through pion photoproduction leads to a flux of protons well below the cosmic-ray data at ~10^18 eV, where the composition is light, unless pions/muons cool before decaying. This suggests a dominant class of accelerator that allows for cosmic rays to escape without significant neutrino yields.

Neutrinos and gravitational waves are the only direct probes of the inner dynamics of a stellar core collapse. They are also the first signals to arrive from a supernova and, if detected, establish the moment when the shock wave is formed that unbinds the stellar envelope and later initiates the optical display upon reaching the stellar surface with a burst of UV and X-ray photons, the shock breakout (SBO). We discuss how neutrino observations can be used to trigger searches to detect the elusive SBO event. Observation of the SBO would provide several important constraints on progenitor structure and the explosion, including the shock propagation time (the duration between the neutrino burst and SBO), an observable that is important in distinguishing progenitor types. Our estimates suggest that next generation neutrino detectors could exploit the overdensity of nearby SNe to provide several such triggers per decade, more than an order of magnitude improvement over the present.

The metallicity of a star strongly affects both its evolution and the properties of the stellar remnant that results from its demise. It is generally accepted that stars with initial masses below ~8 M_sun leave behind white dwarfs and that some sub-population of these lead to Type Ia supernovae. However, it is often tacitly assumed that metallicity has no effect on the rate of SNe Ia. We propose that a consequence of the effects of metallicity is to significantly increase the SN Ia rate in lower-metallicity galaxies, in contrast to previous expectations. This is because lower-metallicity stars leave behind higher-mass white dwarfs, which should be easier to bring to explosion. We first model SN Ia rates in relation to galaxy masses and ages alone, finding that the elevation in the rate of SNe Ia in lower-mass galaxies measured by LOSS is readily explained. However, we then see that models incorporating this effect of metallicity agree just as well. Using the same parameters to estimate the cosmic SN Ia rate, we again find good agreement with data up to z~2. We suggest that this degeneracy warrants more detailed examination of host galaxy metallicities. We discuss additional implications, including for hosts of high-z SNe Ia, the SN Ia delay time distribution, super-Chandrasekhar SNe, and cosmology.

Observations of high-z galaxies and gamma-ray bursts now allow for empirical studies during reionization. However, even deep surveys see only the brightest galaxies at any epoch and must extrapolate to arbitrary lower limits to estimate the total rate of star formation. We first argue that the galaxy populations seen in LBG surveys yield a GRB rate at z > 8 that is an order of magnitude lower than observed. We find that integrating the inferred UV luminosity functions down to M_UV ~ -10 brings LBG- and GRB-inferred SFR density values into agreement up to z ~ 8. GRBs, however, favor a far larger amount of as yet unseen star formation at z > 9. We suggest that the SFR density may only slowly decline out to z ~ 11, in accord with WMAP and Planck reionization results, and that GRBs may be useful in measuring the scale of this multitude of dwarf galaxies.

Isotropy is a key assumption in many models of cosmic-ray electrons and positrons. We find that simulation results imply a critical energy of ~10-1000 GeV above which electrons and positrons can spend their entire lives in streams threading magnetic fields, due to energy losses. This would restrict the number of electron/positron sources contributing at Earth, likely leading to smooth electron and positron spectra, as is observed. For positrons, this could be as few as one, with an enhanced flux that would ease energetics concerns of a pulsar origin of the positron excess, or even zero, bringing dark matter into play. We conclude that ideas about electron/positron propagation based on either isotropic diffusion or turbulent fields must be changed.

The ultrahigh-energy cosmic-ray anisotropies discovered by the Pierre Auger Observatory give the potential to finally address both the particles' origins and properties of the nearby extragalactic magnetic field (EGMF). We examine the implications of the excess of ~ 10^20 eV events around the nearby radio galaxy Centaurus A. We find that, if Cen A is the source of these cosmic rays, the angular distribution of events constrains the EGMF strength within several Mpc of the Milky Way to > 20 nG for an assumed primary proton composition. Our conclusions suggest that either the observed excess is a statistical anomaly or the local EGMF is stronger then conventionally thought. We discuss the implications of this field, including UHECR scattering from more distant sources, time delays from transient sources, and the possibility of using magnetic lensing signatures to attain tighter constraints.

The legacy of solar neutrinos suggests that large neutrino detectors should be sited underground. However, to instead go underwater bypasses the need to move mountains, allowing much larger water Cherenkov detectors. We show that reaching a detector mass scale of ~5 Megatons, the size of the proposed Deep-TITAND, would permit observations of neutrino "mini-bursts" from supernovae in nearby galaxies on a roughly yearly basis, and we develop the immediate qualitative and quantitative consequences. Importantly, these mini-bursts would be detected over backgrounds without the need for optical evidence of the supernova, guaranteeing the beginning of time-domain MeV neutrino astronomy. The ability to identify, to the second, every core collapse in the local Universe would allow a continuous "death watch" of all stars within ~5 Mpc, making practical many previously-impossible tasks in probing rare outcomes and refining coordination of multi-wavelength/multi-particle observations and analysis. These include the abilities to promptly detect otherwise-invisible prompt black hole formation, provide advance warning for supernova shock-breakout searches, define tight time windows for gravitational-wave searches, and identify "supernova impostors" by the non-detection of neutrinos. Observations of many supernovae, even with low numbers of detected neutrinos, will help answer questions about supernovae that cannot be resolved with a single high-statistics event in the Milky Way.

Due to their higher concentrations and small internal velocities, Milky Way subhalos can be at least as important as the smooth halo in accounting for the GeV positron excess via dark matter annihilation. After showing how this can be achieved in various scenarios, including in Sommerfeld models, we demonstrate that, in this case, the diffuse inverse-Compton emission resulting from electrons and positrons produced in substructure leads to a nearly-isotropic signal close to the level of the isotropic GeV gamma-ray background seen by Fermi. Moreover, we show that HESS cosmic-ray electron measurements can be used to constrain multi-TeV internal bremsstrahlung gamma rays arising from annihilation to charged leptons.

At energies above a few TeV, no measurements of the cosmic-ray electron spectrum exist yet. By considering the similarity of air showers induced by electrons and gamma rays as seen by ground-based arrays, we use published limits on isotropic gamma-ray fluxes to place first constraints on the >10 TeV electron spectrum. We demonstrate that, due the proximity of known sources, the flux of such electrons (and positrons) can be large. We show how these smoothly connect to lower-energy positrons measured by PAMELA and relate to exciting new indications from Fermi.

We argue that the recent groundbreaking discovery by Badenes et al. (2009) of a nearby (~50 pc) white dwarf-neutron star (or black hole) binary (SDSS 1257+5428) with a merger timescale ~500 Myr implies that such systems are common; we estimate that there are of order 10^6 in the Galaxy. Although subject to large uncertainties, the nominal derived merger rate is ~5 x 10^-4 per yr in the Milky Way, just ~3-6 and ~20-40 times less than the Type Ia and core-collapse supernova (SN) rates, respectively. This implies that the merger rate is ~0.5-1 x 10^4 per Gpc^3 per yr in the local universe, ~5000-10000 times more than the observed (beaming-uncorrected) long-duration gamma-ray burst (GRB) rate. We estimate the lower limit on the rate in the Galaxy to be >2.5 x 10^-5 per yr at 95% confidence. We briefly discuss the implications of this finding for the census of long- and short-duration GRBs and their progenitors, the frequency of tight binary companions to Type Ib/c SN progenitors, the origin of ultra-high energy cosmic rays (UHECRs), the formation of rapidly rotating neutron stars and ~2-3 M_sun black holes, the census of faint Ia-like SNe, as well as for upcoming and current transient surveys (e.g., LOSS, PTF, LSST), and for high- (LIGO) and low-frequency (LISA) gravitational wave searches.

High-redshift gamma-ray bursts (GRBs) offer an extraordinary opportunity to study aspects of the early Universe, including the cosmic star formation rate (SFR). Motivated by the two recent highest-z GRBs, GRB 080913 at z = 6.7 and GRB 090423 at z = 8.1, and more than four years of Swift observations, we first confirm that the GRB rate does not trace the SFR in an unbiased way. Correcting for this, we find that the implied SFR to beyond z = 8 is consistent with LBG-based measurements after accounting for unseen galaxies at the faint end of the UV luminosity function. We show that this provides support for the integrated star formation in the range 6 < z < 8 to have been alone sufficient to reionize the Universe.

The fraction of ionizing photons that escape their host galaxies and so are able to ionize hydrogen in the inter-galactic medium (IGM) is a critical parameter in analyses of the reionization era and early galaxy formation. Studies of the reionization history normally suffer from a degeneracy between the unknown values for the efficiency with which high redshift galaxies turn mass into stars and the escape fraction of ionizing photons. Recent gamma-ray burst (GRB) measurements of the star formation rate density during reionization provide the first opportunity to break this degeneracy. We confront a semi-analytic model for reionization with the GRB-derived star formation rate, as well as observations of the Ly-alpha forest and the CMB. Assuming that UV photons produced in star-forming galaxies dominate the reionization process, we show that the escape fraction of ionizing photons from high redshift galaxies is ~5% [log(f_esc)=-1.35+/-0.15 (68%)] for our fiducial model. This value is reasonably stable against uncertainties in the modeling, including the implementation of radiative feedback, the possibility of an evolving escape fraction, and the unknown shape of the IMF, which in sum contribute ~0.2 dex of additional systematic uncertainty on the value of escape fraction.

The Geminga pulsar has long been one of the most intriguing MeV-GeV gamma-ray point sources. We examine the implications of the recent Milagro detection of extended, multi-TeV gamma-ray emission from Geminga, finding that this reveals the existence of an ancient, powerful cosmic-ray accelerator that can plausibly account for the multi-GeV positron excess that has evaded explanation. We explore a number of testable predictions for gamma-ray and electron/positron experiments (up to ~100 TeV) that can confirm the first "direct" detection of a cosmic-ray source.

The progenitors of SN 2008S and the 2008 transient in NGC300 were dust-enshrouded, with extremely red mid-infrared (MIR) colors and relatively low luminosities. The transients were optically faint (-13 < M_V < -15) compared to normal core-collapse supernovae (SNe), and their spectra exhibited narrow emission lines. These events are unique among transient-progenitor pairs and hence constitute a new class. Whether they are true SNe or bright massive-star eruptions, we argue that their rate is ~20% of the SN rate. This fact is remarkable in light of the observation that a very small fraction of all massive stars have the MIR colors of the SN 2008S and NGC300 progenitors, as we show using MIR and optical luminosity, color, and variability properties of massive stars in M33. We find that the fraction of massive stars with colors consistent with these progenitors is 1/10000. In fact, only < 10 similar objects exist in M33 - all of which lie at the luminous red extremum of the AGB sequence. That these transients are relatively common with respect to SNe, while their progenitors are rare compared to the massive star population, implies that the dust-enshrouded phase is a short-lived phase in the lives of many massive stars. This shrouded epoch can occur only in the last ~10^4 yr before explosion, be it death or merely eruption. We discuss the implications of this finding for the evolution and census of ``low-mass'' massive stars (8-12 Msun), and we connect it with theoretical discussions of electron-capture SNe and the explosive birth of white dwarfs. A systematic census with (warm) Spitzer of galaxies in the local universe for analogous progenitors would significantly improve our knowledge of this channel to massive stellar explosions, and potentially to others with obscured progenitors. (Abridged)

While the high-z frontier of star formation rate (SFR) studies has advanced rapidly, direct measurements beyond z ~ 4 remain difficult, as shown by significant disagreements among different results. Gamma-ray bursts, owing to their brightness and association with massive stars, offer hope of clarifying this situation, provided that the GRB rate can be properly related to the SFR. The Swift GRB data reveal an increasing evolution in the GRB rate relative to the SFR at intermediate z; taking this into account, we use the highest-z GRB data to make a new determination of the SFR at z = 4-7. Our results exceed the lowest direct SFR measurements, and imply that no steep drop exists in the SFR up to at least z ~ 6.5. We discuss the implications of our result for cosmic reionization, the efficiency of the universe in producing stellar-mass black holes, and ``GRB feedback'' in star-forming hosts.

The association of long gamma-ray bursts with supernovae naturally suggests that the cosmic GRB rate should trace the star formation history. Finding otherwise would provide important clues concerning these rare, curious phenomena. Using a new estimate of Swift GRB energetics to construct a sample of 36 luminous GRBs with redshifts in the range z=0-4, we find evidence of enhanced evolution in the GRB rate, with ~4 times as many GRBs observed at z~4 than expected from star formation measurements. This direct and empirical demonstration of needed additional evolution is a new result. It is consistent with theoretical expectations from metallicity effects, but other causes remain possible, and we consider them systematically.

A number of theories, spanning a wide range of mass scales, predict dark matter candidates that have lifetimes much longer than the age of the universe, yet may produce a significant flux of gamma rays in their decays today. We constrain such late decaying dark matter scenarios model-independently by utilizing gamma-ray line emission limits from the Galactic Center region obtained with the SPI spectrometer on INTEGRAL, and the determination of the isotropic diffuse photon background by SPI, COMPTEL and EGRET observations. We show that no more than ~5% of the unexplained MeV background can be produced by late dark matter decays either in the Galactic halo or cosmological sources.

We report the discovery of the progenitor of the recent type IIn SN 2008S in the nearby galaxy NGC 6946. Surprisingly, it was not found in deep, pre-explosion optical images of its host galaxy taken with the Large Binocular Telescope, but only through examination of archival Spitzer mid-IR data. A source coincident with the SN 2008S position is clearly detected in the 4.5, 5.8, and 8.0 micron IRAC bands, showing no evident variability in the three years prior to the explosion, yet is undetected at 3.6 and 24 micron. The distinct presence of ~440 K dust, along with stringent LBT limits on the optical fluxes, suggests that the progenitor of SN 2008S was engulfed in a shroud of its own dust. The inferred luminosity of 3.5x10^4 Lsun implies a modest mass of ~10 Msun. We conclude that objects like SN 2008S are not exclusively associated with the deaths or outbursts of very massive eta Carinae-like objects. This conclusion holds based solely on the optical flux limits even if our identification of the progenitor with the mid-IR source is incorrect.

Extragalactic transient searches have historically been limited to looking for the appearance of new sources such as supernovae. It is now possible to carry out a new kind of survey that will do the opposite, that is, search for the disappearance of massive stars. This will entail the systematic observation of galaxies within a distance of 10 Mpc in order to watch ~10^6 supergiants. Reaching this critical number ensures that something will occur yearly, since these massive stars must end their lives with a core collapse within ~10^6 years. Using deep imaging and image subtraction it is possible to determine the fates of these stars whether they end with a bang (supernova) or a whimper (fall out of sight). Such a survey would place completely new limits on the total rate of all core collapses, which is critical for determining the validity of supernova models. It would also determine the properties of supernova progenitors, better characterize poorly understood optical transients, such as eta Carina-like mass ejections, find and characterize large numbers of Cepheids, luminous blue variables and eclipsing binaries, and allow the discovery of any new phenomena that inhabit this relatively unexplored parameter space.

Gamma-ray bursts, which are among the most violent events in the universe, are one of the few viable candidates to produce ultrahigh energy cosmic rays. Recently, observations have revealed that GRBs generally originate from metal-poor galaxies and do not directly trace cosmic star formation, as might have been assumed from their association with core-collapse supernovae. Several implications follow from these findings. The redshift distribution of observed GRBs is expected to peak at higher redshift (compared to cosmic star formation), which is supported by the mean redshift of the Swift GRB sample, <z>~3. If GRBs are, in fact, the source of the observed UHECR, then cosmic-ray production would evolve with redshift in a stronger fashion than has been previously suggested. This necessarily leads, through the GZK process, to an enhancement in the flux of cosmogenic neutrinos, providing a near-term approach for testing the gamma-ray burst-cosmic ray connection with ongoing and proposed UHE neutrino experiments.

Recent Milagro observations of the Cygnus region have revealed both diffuse TeV gamma-ray emission and a bright and extended TeV source, MGRO J2019+37, which seems to lack an obvious counterpart at other wavelengths. Additional study of this curious object also promises to provide important clues concerning one of the Milky Way's most active environments. We point out some of the principal facts involved by following three modes of attack. First, to gain insight into this mysterious source, we consider its relation to known objects in both the Cygnus region and the rest of the Galaxy. Second, we find that a simple hadronic model can easily accommodate Milagro's flux measurement (which is at a single energy), as well as other existing observations spanning nearly seven orders of magnitude in gamma-ray energy. Third, since a hadronic gamma-ray spectrum necessitates an accompanying TeV neutrino flux, we show that IceCube observations may provide the first direct evidence of a Galactic cosmic-ray accelerator.

Recent observations, particularly from the HESS Collaboration, have revealed rich Galactic populations of TeV gamma-ray sources, including a collection unseen in other wavelengths. Many of these gamma-ray spectra are well measured up to ~10 TeV, where low statistics make observations by air Cerenkov telescopes difficult. To understand these mysterious sources, especially at much higher energies--where a cutoff should eventually appear--new techniques are needed. We point out the following: (1) For a number of sources, it is very likely that pions, and hence TeV neutrinos, are produced; (2) As a general point, neutrinos should be a better probe of the highest energies than gamma rays, due to increasing detector efficiency; and (3) For several specific sources, the detection prospects for km^3 neutrino telescopes are very good, about 1-10 events/year, with low atmospheric neutrino background rates above reasonable energy thresholds. Such signal rates, as small as they may seem, will allow neutrino telescopes to powerfully discriminate between models for the Galactic TeV sources, with important consequences for our understanding of cosmic-ray production.