We present Breakthrough Listen's "Exotica" Catalog as the centerpiece of our efforts to expand the diversity of targets surveyed in the Search for Extraterrestrial Intelligence (SETI). As motivation, we introduce the concept of survey breadth, the diversity of objects observed during a program. Several reasons for pursuing a broad program are given, including increasing the chance of a positive result in SETI, commensal astrophysics, and characterizing systematics. The Exotica Catalog is an 865 entry collection of 737 distinct targets intended to include "one of everything" in astronomy. It contains four samples: the Prototype sample, with an archetype of every known major type of non-transient celestial object; the Superlative sample of objects with the most extreme properties; the Anomaly sample of enigmatic targets that are in some way unexplained; and the Control sample with sources not expected to produce positive results. As far as we are aware, this is the first object list in recent times with the purpose of spanning the breadth of astrophysics. We share it with the community in hopes that it can guide treasury surveys and as a general reference work. Accompanying the catalog is extensive discussion of classification of objects and a new classification system for anomalies. We discuss how we intend to proceed with observations in the catalog, contrast it with our extant Exotica efforts, and suggest similar tactics may be applied to other programs.

As it stands today, the search for extraterrestrial intelligence (SETI) is highly dependent on our ability to detect interesting candidate signals, or technosignatures, in radio telescope observations and distinguish these from human radio frequency interference (RFI). Current signal search pipelines look for signals in spectrograms of intensity as a function of time and frequency (which can be thought of as images), but tend to do poorly in identifying multiple signals in a single data frame. This is especially apparent when there are dim signals in the same frame as bright, high signal-to-noise ratio (SNR) signals. In this work, we approach this problem using convolutional neural networks (CNN) to localize signals in synthetic observations resembling data collected by Breakthrough Listen using the Green Bank Telescope. We generate two synthetic datasets, the first with exactly one signal at various SNR levels and the second with exactly two signals, one of which represents RFI. We find that a residual CNN with strided convolutions and using multiple image normalizations as input outperforms a more basic CNN with max pooling trained on inputs with only one normalization. Training each model on a smaller subset of the training data at higher SNR levels results in a significant increase in model performance, reducing root mean square errors by a factor of 5 at an SNR of 25 dB. Although each model produces outliers with significant error, these results demonstrate that using CNNs to analyze signal location is promising, especially in image frames that are crowded with multiple signals.

The discovery of ubiquitous habitable extrasolar planets, combined with revolutionary advances in instrumentation and observational capabilities, has ushered in a renaissance in the search for extra-terrestrial intelligence (SETI). Large scale SETI activities are now underway at numerous international facilities. The Five-hundred-meter Aperture Spherical radio Telescope (FAST) is the largest single-aperture radio telescope in the world, well positioned to conduct sensitive searches for radio emission indicative of exo-intelligence. SETI is one of the five key science goals specified in the original FAST project plan. A collaboration with the Breakthrough Listen Initiative has been initiated in 2016 with a joint statement signed both by Dr. Jun Yan, the then director of the National Astronomical Observatories, Chinese Academy of Sciences (NAOC), and Dr. Peter Worden, the Chairman of the Breakthrough Prize Foundation. In this paper, we highlight some of the unique features of FAST that will allow for novel SETI observations. We identify and describe three different signal types indicative of a technological source, namely, narrow-band, wide-band artificially dispersed, and modulated signals. We here propose observations with FAST to achieve sensitivities never before explored.

Breakthrough Listen (BL) is a ten-year initiative to search for signatures of technologically capable life beyond Earth via radio and optical observations of the local Universe. A core part of the BL program is a comprehensive survey of 1702 nearby stars at radio wavelengths (1-10 GHz). Here, we report on observations with the 64-m CSIRO Parkes radio telescope in New South Wales, Australia, and the 100-m Robert C. Byrd Green Bank radio telescope in West Virginia, USA. Over 2016 January to 2019 March, a sample of 1138 stars was observed at Green Bank using the 1.10-1.90 GHz and 1.80-2.80 GHz receivers, and 189 stars were observed with Parkes over 2.60-3.45 GHz. We searched these data for the presence of engineered signals with Doppler-acceleration drift rates between -4 to 4 Hz/s. Here, we detail our data analysis techniques and provide examples of detected events. After excluding events with characteristics consistent with terrestrial radio interference, we are left with zero candidates. Given the sensitivity of our observations, we can put an upper limit on the power of potential radio transmitters at these frequencies at 2x10^12 W, and 9x10^12 W for GBT and Parkes respectively. These observations constitute the most comprehensive search over 1.10-3.45 GHz for technosignatures to date. All data products, totalling ~219 TB, are available for download as part of the first BL data release (DR1), as described in a companion paper (Lebofsky et. al., 2019)

Low mass X-ray binaries (LMXBs) containing neutron stars are both extremely luminous and compact, emitting up to ~10^6 L_sun within a kilometer-scale boundary layer. This combination allows for easy modulation, motivating X-ray SETI. When X-ray lenses smaller than planets (100 - 1,000 km) magnify the LMXB boundary layer, it brightens by a factor of several thousand for about a second. In addition, there should be occultation events where the neutron star is blocked out. Passive X-ray lenses could require little maintenance and the LMXB light source itself shines for millions of years, serving as an effective beacon for interstellar communication. A very large number of lenses would be needed to ensure detection from all directions, however, and gathering material to construct them could be very difficult. Avoiding collisions between lenses and aiming them pose additional challenges. Both "lens flares" and eclipses of LMXBs are easily detectable in the Galaxy, although they would be rare events, occurring once per decade. Our X-ray instruments could detect the flares to several Mpc, but it is unlikely they would be observing the LMXB during a flare.

The star EPIC 249706694 (HD 139139) was found to exhibit 28 transit-like events over an 87 day period during the Kepler mission's K2 Campaign 15 (Rappaport et al. 2019). These events did not fall into an identifiable pattern and could not be explained by a multitude of transit scenarios explored by the authors. We conduct follow-up observations at C-band frequencies with the Green Bank Telescope as part of the ongoing Breakthrough Listen search for technosignatures. We search for narrow band signals above a signal-to-noise threshold of 10 and with Doppler drift rates within +-5 Hz/s. We detect no evidence of technosignatures from EPIC 249706694 and derive an upper limit for the EIRP (Equivalent Isotropic Radiated Power) of putative transmissions to be 10 TW.

Breakthrough Listen is the most comprehensive and sensitive search for extraterrestrial intelligence (SETI) to date, employing a collection of international observational facilities including both radio and optical telescopes. During the first three years of the Listen program, thousands of targets have been observed with the Green Bank Telescope (GBT), Parkes Telescope and Automated Planet Finder. At GBT and Parkes, observations have been performed ranging from 700 MHz to 26 GHz, with raw data volumes averaging over 1PB / day. A pseudo-real time software spectroscopy suite is used to produce multi-resolution spectrograms amounting to approximately 400 GB hr^-1 GHz^-1 beam^-1. For certain targets, raw baseband voltage data is also preserved. Observations with the Automated Planet Finder produce both 2-dimensional and 1-dimensional high resolution (R~10^5) echelle spectral data. Although the primary purpose of Listen data acquisition is for SETI, a range of secondary science has also been performed with these data, including studies of fast radio bursts. Other current and potential research topics include spectral line studies, searches for certain kinds of dark matter, probes of interstellar scattering, pulsar searches, radio transient searches and investigations of stellar activity. Listen data are also being used in the development of algorithms, including machine learning approaches to modulation scheme classification and outlier detection, that have wide applicability not just for astronomical research but for a broad range of science and engineering. In this paper, we describe the hardware and software pipeline used for collection, reduction, archival, and public dissemination of Listen data. We describe the data formats and tools, and present Breakthrough Listen Data Release 1.0 (BLDR 1.0), a defined set of publicly-available raw and reduced data totalling 1 PB.

The discovery of the ubiquity of habitable extrasolar planets, combined with revolutionary advances in instrumentation and observational capabilities, have ushered in a renaissance in the millenia-old quest to answer our most profound question about the Universe and our place within it - Are we alone? The Breakthrough Listen Initiative, announced in July 2015 as a 10-year 100M USD program, is the most comprehensive effort in history to quantify the distribution of advanced, technologically capable life in the universe. In this white paper, we outline the status of the on-going observing campaign with our primary observing facilities, as well as planned activities with these instruments over the next few years. We also list collaborative facilities which will conduct searches for technosignatures in either primary observing mode, or commensally. We highlight some of the novel analysis techniques we are bringing to bear on multi-petabyte data sets, including machine learning tools we are deploying to search for a broader range of technosignatures than was previously possible.

Glints of light from specular reflection of the Sun are a technosignature of artificial satellites. If extraterrestrial intelligences have left artifacts in the Solar System, these may include flat mirror-like surfaces that also can glint. I describe the characteristics of the resulting flashes. An interplanetary mirror will appear illuminated for several hours, but if it is rotating, its glint may appear as a train of optical pulses. The resulting glints can be very bright, but they will be seen only if the mirror happens to reflect sunlight to the Earth. The detection of large mirrors is limited mainly by the fraction oriented to reflect sunlight toward Earth. I give rough calculations for the expected reach of each exposure of Pan-STARRS1, LSST, and Evryscope for mirror glints. A single exposure of Pan-STARRS1 has an effective reach of 10^-9 - 10^-7 AU^3 for interplanetary mirrors with effective areas of 10 m^2, depending on rotation rate. Over several years, Pan-STARRS1 might accumulate a reach ~10^5 times greater than this, as it tiles the sky and different mirrors enter and exit a favorable geometry.

The Search for Extraterrestrial Intelligence has so far come up negative for Kardashev Type III societies that capture all starlight from a galaxy. One possible reason is that shrouding a star in a megastructure is prohibitively expensive. Most of a galaxy's starlight comes from bright stars, which would require structures even larger than the classical Dyson sphere to enclose. Using a custom spectral synthesis code, I calculate what happens to the spectrum and colors of a galaxy when only stars below a luminosity L_min are cloaked. I find the photometric signatures of galaxies with L_min <= 1 L_sun are minor, especially for blue, galaxies with continuing star formation. Larger luminosity thresholds (>~ 30 L_sun) result in galaxies with unnatural colors and luminosities. Galaxies observed in NIR and galaxies without recent star formation observed at UV-NIR wavelengths become redder than uncloaked galaxies as L_min increases. Recently star-forming galaxies get bluer in UV and blue light when they are cloaked, with colors similar to quasars but very low luminosities. By selecting on color, we may find Type III societies in large photometric surveys. I discuss how different metallicities, ages, and initial mass functions affect the results.

Boyajian's Star (KIC 8462852) undergoes mysterious, irregular eclipses that aren't yet explained. It also appears to have dimmed over a time of several years, possibly decades. I show that Kepler's detection of a phenomenon with a duration of t_anom is only expected if it occurs at a mean rate of >~ 30 Gyr^(-1) (t_anom / 100 yr)^(-1) for each Kepler target and K2 star. If true, the phenomenon occurs hundreds of times during the lifespan of its host stars. Obscuration by the interstellar medium remains a plausible explanation, since it doesn't actually affect the host star. An intervening cloud is consistent with the lack of an observed submillimeter excess but would be abnormally dilute.

It is unclear how frequently life and intelligence arise on planets. I consider a Bayesian prior for the probability P(ETI) that intelligence evolves at a suitable site, with weight distributed evenly over ln(1 - ln P(ETI)). This log log prior can handle a very wide range of P(ETI) values, from 1 to 10^(-10^122), while remaining responsive to evidence about extraterrestrial societies. It is motivated by our uncertainty in the number of conditions that must be fulfilled for intelligence to arise, and it is related to considerations of information, entropy, and state space dimensionality. After setting a lower limit to P(ETI) from the number of possible genome sequences, I calculate a Bayesian confidence of 18% that aliens exist within the observable Universe. With different assumptions about the minimum P(ETI) and the number of times intelligence can appear on a planet, this value falls between 1.4% and 47%. Overall, the prior leans towards our being isolated from extraterrestrial intelligences, but indicates that we should not be confident of this conclusion. I discuss the implications of the prior for the Search for Extraterrestrial Intelligence, concluding that searches for interstellar probes from nearby societies seem relatively effective. I also discuss the possibility of very small probabilities allowed by the prior for the origin of life and the Fermi Paradox, and note that similar priors might be constructed for interesting complex phenomena in general.

[Abridged] Whether technological societies remain small and planet-bound like our own, or ultimately span across galaxies is an open question in the Search for Extraterrestrial Intelligence. Societies that engineer on a galactic scale are classified as Type III on Kardashev's scale. I argue that Type III societies can take the form of blackboxes, entire galaxies veiled in an opaque screen. A blackbox has a temperature that is just above that of the cosmic microwave background. The screen can be made from artificial dust pervading the galaxy. I show that there is enough material in galaxies to build blackboxes if the dust is fashioned into dipole antennas. The thermal emission of a blackbox makes it a bright microwave source. I examine the Planck Catalog of Compact Sources to constrain the abundance of blackboxes. None of the 100 GHz sources has the spectrum expected of a blackbox. The null result rules out shrouded galaxy clusters out to z ~ 1 and shrouded Milky Ways out to (comoving) 700 Mpc. The reach of the results includes 3 million galaxies containing an estimated 300 quadrillion terrestrial planets, as well as tens of thousands of galaxy clusters. Combined with the null results from other searches for Type III societies, I conclude that they are so rare that they basically do not exist within the observable Universe. A hypothesis of "Cosmic Pessimism" is discussed, in which we are alone, our long-term chances for survival are slim, and if we do survive, our future history will be checkered. Our loneliness is suggested by the lack of Type III societies. I discuss the remaining forms of Type III societies not yet well constrained by observation. I argue that the ease of building blackboxes on planetary and Solar System scales may lead, within a few centuries, to environmental catastrophes vastly more devastating than anything we are doing now, boding ill for us.

What is the meaning of the Fermi Paradox -- are we alone or is starfaring rare? Can general relativity be united with quantum mechanics? The searches for answers to these questions could intersect. It is known that an accelerator capable of energizing particles to the Planck scale requires cosmic proportions. The energy required to run a Planck accelerator is also cosmic, of order 100 M_sun c^2 for a hadron collider, because the natural cross section for Planck physics is so tiny. If aliens are interested in fundamental physics, they could resort to cosmic engineering for their experiments. These colliders are detectable through the vast amount of "pollution" they produce, motivating a YeV SETI program. I investigate what kinds of radiation they would emit in a fireball scenario, and the feasibility of detecting YeV radiation at Earth, particularly YeV neutrinos. Although current limits on YeV neutrinos are weak, Kardashev 3 YeV neutrino sources appear to be at least 30--100 Mpc apart on average, if they are long-lived and emit isotropically. I consider the feasibility of much larger YeV neutrino detectors, including an acoustic detection experiment that spans all of Earth's oceans, and instrumenting the entire Kuiper Belt. Any detection of YeV neutrinos implies an extraordinary phenomenon at work, whether artificial and natural. Searches for YeV neutrinos from any source are naturally commensal, so a YeV neutrino SETI program has value beyond SETI itself, particularly in limiting topological defects. I note that the Universe is very faint in all kinds of nonthermal radiation, indicating that cosmic engineering is extremely rare.

The majority of cosmic rays (CRs) generated by star-forming galaxies escape them and enter the intergalactic medium (IGM). Galactic wind termination shocks might also accelerate CRs. I show that the mean pressure of these CRs can reach to within an order of magnitude of the mean Lyman-alpha forest thermal pressure. At z >~ 1, their pressure may have even been dominant. I also demonstrate that, whichever IGM phase the CRs reside in, they contribute significantly to its pressure if its temperature is ~$10^4$ K, as long as pionic and Coulomb losses are negligible. Where CRs end up depends on the structure and strength of intergalactic magnetic fields. I argue that CRs end up at least 30 kpc from their progenitor galaxies. CRs may self-confine in the IGM to the sound speed, generating >~ $10^{-13}$ G magnetic fields. These considerations imply the existence and importance of a nonthermal IGM.

The basic principle of astronomical interferometry is to derive the angular distribution of radiation in the sky from the Fourier transform of the electric field on the ground. What is so special about the Fourier transform? Nothing, it turns out. I consider the possibility of performing other transforms on the electric field with digital technology. The Fractional Fourier Transform (FrFT) is useful for interpreting observations of sources that are close to the interferometer (in the atmosphere for radio interferometers). Essentially, applying the FrFT focuses the array somewhere nearer than infinity. Combined with the other Linear Canonical Transforms, any homogeneous linear optical system with thin elements can be instantiated. The time variation of the electric field can also be decomposed into other bases besides the Fourier modes, which is especially useful for dispersed transients or quick pulses. I discuss why the Fourier basis is so commonly used, and suggest it is partly because most astrophysical sources vary slowly in time.

Imaging Atmosphere Cherenkov Telescopes (IACT) are arrays of very large optical telescopes that are well-suited for rapid photometry of bright sources. I investigate their potential in observing stellar occultations by small objects in the outer Solar System, Transjovian Objects (TJOs). These occultations cast diffraction patterns on the Earth. Current IACT arrays are capable of detecting objects smaller than 100 metres in radius in the Kuiper Belt and 1 km radius out to 5000 AU. The future Cherenkov Telescope Array (CTA) will have even greater capabilities. Because the arrays include several telescopes, they can potentially measure the speeds of TJOs without degeneracies, and the sizes of the TJOs and background stars. I estimate the achievable precision using a Fisher matrix analysis. With CTA, the precisions of these parameter estimations will be as good as a few percent. I consider how often detectable occultations occur by members of different TJO populations, including Centaurs, Kuiper Belt Objects (KBOs), Oort cloud objects, and satellites and Trojans of Uranus and Neptune. The great sensitivity of IACT arrays means that they likely detect KBO occultations once every O(10) hours when looking near the ecliptic. IACTs can also set useful limits on many other TJO populations.

The enhanced star formation in the inner 100 pc of the Galaxy launches a superwind at ~1600 km s$^{-1}$ for M82-like parameters. The ram pressure of the wind is very low compared to more powerful starburst winds. I show that halo gas stops the wind a few kpc from the Galactic Centre. I suggest that the termination shock accelerates cosmic rays, and that the resulting Inverse Compton gamma rays are visible as the Fermi Bubbles. The Bubbles are then wind bubbles, which the starburst can inflate within 10 Myr. They can remain in steady state as long as the starburst lasts. The shock may accelerate PeV electrons and EeV protons. The Bubbles may be analogues of galactic wind termination shocks in the intergalactic medium. I discuss the advantages and problems of this model. I note that any jets from Sgr A* must burrow through the starburst wind bubble before reaching the halo gas, which could affect the early evolution of such jets.

Short lived radionuclides (SLRs) like $^{26}$Al are synthesized by massive stars and are a byproduct of star formation. The abundances of SLRs in the gas of a star-forming galaxy are inversely proportional to the gas consumption time. The rapid evolution of specific star formation rate (SSFR) of normal galaxies implies they had mean SLR abundances ~3--10 times higher at z = 2. During the epoch of Solar system formation, the background SLR abundances of the Galaxy were up to twice as high as at present, if SLR yields from massive stars do not depend on metallicity. If SLRs are homogenized in the gas of galaxies, the high SSFRs of normal galaxies can partly explain the elevated abundance of SLRs like $^{60}$Fe and $^{26}$Al in the early Solar system. Starburst galaxies have much higher SSFRs still, and have enormous mean abundances of $^{26}$Al ($^{26}$Al/$^{27}$Al ~ $10^{-3}$ for Solar metallicity gas). The main uncertainty is whether the SLRs are mixed with the star-forming molecular gas: they could be trapped in hot gas and decay before entering the colder phases, or be blown out by starburst winds. I consider how variability in star-formation rate affects the SLR abundances, and I discuss how SLR transport may differ in these galaxies. The enhanced $^{26}$Al of starbursts might maintain moderate ionization rates ($10^{-18}$ -- $10^{-17}$ s$^{-1}$), possibly dominating ionization in dense clouds not penetrated by cosmic rays. Similar ionization rates would be maintained in protoplanetary discs of starbursts, if the SLRs are well-mixed, and the radiogenic heating of planetesimals would likewise be much higher. In this way, galaxy evolution can affect the geological history of planetary systems.

While star-forming galaxies could be major contributors to the cosmic GeV $\gamma$-ray background, they are expected to be MeV-dim because of the "pion bump" falling off below ~100 MeV. However, there are very few observations of galaxies in the MeV range, and other emission processes could be present. We investigate the MeV background from star-forming galaxies by running one-zone models of cosmic ray populations, including Inverse Compton and bremsstrahlung, as well as nuclear lines (including $^{26}$Al), emission from core-collapse supernovae, and positron annihilation emission, in addition to the pionic emission. We use the Milky Way and M82 as templates of normal and starburst galaxies, and compare our models to radio and GeV--TeV $\gamma$-ray data. We find that (1) higher gas densities in high-z normal galaxies lead to a strong pion bump, (2) starbursts may have significant MeV emission if their magnetic field strengths are low, and (3) cascades can contribute to the MeV emission of starbursts if they emit mainly hadronic $\gamma$-rays. Our fiducial model predicts that most of the unresolved GeV background is from star-forming galaxies, but this prediction is uncertain by an order of magnitude. About ~2% of the claimed 1 MeV background is diffuse emission from star-forming galaxies; we place a firm upper limit of <~10% based on the spectral shape of the background. The star-formation contribution is constrained to be small, because its spectrum is peaked, while the observed background is steeply falling with energy through the MeV-GeV range.

We consider implications of the IceCube excess for hadronuclear (pp) scenarios of neutrino sources such as galaxy clusters/groups and star-forming galaxies. Since the observed neutrino flux is comparable to the diffuse gamma-ray background flux obtained by Fermi, we place new, strong upper limits on the source spectral index, {\Gamma}<2.1-2.2. In addition, the new IceCube data imply that these sources contribute at least 30%-40% of the diffuse gamma-ray background in the 100 GeV range and even ~100% for softer spectra. Our results, which are insensitive to details of the pp source models, are one of the first strong examples of the multimessenger approach combining the measured neutrino and gamma-ray fluxes. The pp origin of the IceCube excess can further be tested by constraining {\Gamma} with sub-PeV neutrino observations, by unveiling the sub-TeV diffuse gamma-ray background and by observing such pp sources with TeV gamma-ray detectors. We also discuss specific pp source models with a multi-PeV neutrino break/cutoff, which are consistent with the current IceCube data.

Recent work has both illuminated and mystified our attempts to understand cosmic rays (CRs) in starburst galaxies. I discuss my new research exploring how CRs interact with the ISM in starbursts. Molecular clouds provide targets for CR protons to produce pionic gamma rays and ionization, but those same losses may shield the cloud interiors. In the densest molecular clouds, gamma rays and Al-26 decay can provide ionization, at rates up to those in Milky Way molecular clouds. I then consider the free-free absorption of low frequency radio emission from starbursts, which I argue arises from many small, discrete H II regions rather than from a "uniform slab" of ionized gas, whereas synchrotron emission arises outside them. Finally, noting that the hot superwind gas phase fills most of the volume of starbursts, I suggest that it has turbulent-driven magnetic fields powered by supernovae, and that this phase is where most synchrotron emission arises. I show how such a scenario could explain the far-infrared radio correlation, in context of my previous work. A big issue is that radio and gamma-ray observations imply CRs also must interact with dense gas. Understanding how this happens requires a more advanced understanding of turbulence and CR propagation.

Frequent supernova explosions in compact starburst regions are a main shaper of these regions' interstellar media (ISM). In most starbursts, the supernova remnants blast open a hot phase that fills the regions and launches a superwind. Denser starbursts are too overpressured for hot wind formation, but supernovae still stir up the ISM. I argue that supernovae power ubiquitous turbulence through each of the starburst ISM phases, including the hot wind, and that a fluctuation dynamo amplifies magnetic fields until they are in equipartition with the turbulence. Supernovae can drive turbulence with speeds of ~1000 km/s in the hot wind and ~20 km/s in the cold molecular gas, depending on the outer scale. I predict magnetic field strengths of 70 muG in the Galactic Center starburst, 300 muG in M82 and NGC 253, and 2 mG in Arp 220's nuclei. The mean magnetic field strengths are a few times stronger in molecular gas than in hot winds, but do not vary strongly with density within a starburst. I explain how the dominance of supernova-driven turbulence leads to near equipartition between the components of starburst ISM. I also consider implications for cosmic ray (CR) diffusion in starbursts. The high amounts of power cascading to small scales could confine CRs very effectively in starbursts, so much that CR transport is dominated by the flow of gas rather than diffusion through the gas. In addition, I discuss the effects of turbulence on X-ray line width, the far-infrared--radio correlation, observed radio polarization, and Faraday rotation measures. Finally, I discuss the many questions raised regarding the physics of turbulence in starbursts.

[Abdriged] The origin of the diffuse hard X-ray (2 - 10 keV) emission from starburst galaxies is a long-standing problem. We suggest that synchrotron emission of 10 - 100 TeV electrons and positrons (e+/-) can contribute to this emission, because starbursts have strong magnetic fields. We consider three sources of e+/- at these energies: (1) primary electrons directly accelerated by supernova remnants; (2) pionic secondary e+/- created by inelastic collisions between CR protons and gas nuclei in the dense ISMs of starbursts; (3) pair e+/- produced between the interactions between 10 - 100 TeV gamma-rays and the intense far-infrared (FIR) radiation fields of starbursts. We create one-zone steady-state models of the CR population in the Galactic Center (R <= 112 pc), NGC 253, M82, and Arp 220's nuclei, assuming a power law injection spectrum for electrons and protons. We compare these models to extant radio and GeV and TeV gamma-ray data for these starbursts, and calculate the diffuse synchrotron X-ray and Inverse Compton (IC) luminosities of these starbursts. If the primary electron spectrum extends to ~PeV energies and has a proton/electron injection ratio similar to the Galactic value, we find that synchrotron contributes 2 - 20% of their unresolved, diffuse hard X-ray emission. Inverse Compton emission is likewise a minority of the unresolved X-ray emission in these starbursts, from 0.1% in the Galactic Center to 10% in Arp 220's nuclei. We also model generic starbursts, including submillimeter galaxies, in the context of the FIR--X-ray relation, finding that up to 2% in the densest starbursts with our fiducial assumptions. Neutrino and TeV gamma-ray data can further constrain the synchrotron X-ray emission of starbursts. Our models do not constrain hard synchrotron X-ray emission from any additional hard components of primary e+/- from sources like pulsars in starbursts.

The low frequency radio emission of starburst galaxies is informative, but it can be absorbed in several ways. Most importantly, starburst galaxies are home to many H II regions, whose free-free absorption blocks low frequency radio waves. These H II regions are discrete objects, but most multiwavelength models of starbursts assume a uniform medium of ionized gas, if they include the absorption at all. I calculate the effective absorption coefficient of H II regions in starbursts, which is ultimately a cross section times the density of H II regions. The cross sections are calculated by assuming that H II regions are Str\"omgren spheres. The coefficient asymptotes to a constant value at low frequencies, because H II regions partially cover the starburst, and are buried part way into the starburst's synchrotron emitting material. Considering Str\"omgren spheres around either OB stars or Super Star Clusters, I demonstrate the method by fitting to the low frequency radio spectrum of M82. I discuss implications of the results for synchrotron spectrum shape, H II region pressure, and free-free emission as a star-formation rate indicator. However, these results are preliminary, and could be affected by systematics. I argue that there is no volume-filling warm ionized medium in starbursts, and that H II regions may be the most important absorption process down to ~10 MHz. Future data at low and high radio frequency will improve our knowledge of the ionized gas.

Equipartition arguments provide an easy way to find a characteristic scale for the magnetic field from radio emission, by assuming the energy densities in cosmic rays and magnetic fields are the same. Yet most of the cosmic ray content in star-forming galaxies is in protons, which are invisible in radio emission. Therefore, the argument needs assumptions about the proton spectrum, typically that of a constant proton/electron ratio. In some environments, particularly starburst galaxies, the reasoning behind these assumptions does not necessarily hold: secondary pionic positrons and electrons may be responsible for most of the radio emission, and strong energy losses can alter the proton/electron ratio. We derive an equipartition expression that should work in a hadronic loss-dominated environment like starburst galaxies. Surprisingly, despite the radically different assumptions from the classical equipartition formula, numerically the results for starburst magnetic fields are similar. We explain this fortuitous coincidence using the energetics of secondary production and energy loss times. We show that these processes cause the proton/electron ratio to be ~100 for GHz-emitting electrons in starbursts.

Cosmic rays are appealing as a source of ionization in starburst galaxies because of the great columns they can penetrate, but in the densest regions of starbursts, they may be stopped by pion production and ionization energy losses. I argue that gamma rays are the source of ionization in the deepest molecular clouds of dense starbursts, creating Gamma-Ray Dominated Regions (GRDRs). Gamma rays are not deflected by magnetic fields, have a luminosity up to ~1/3 that of the injected cosmic rays, and can easily penetrate column depths of ~100 g/cm^2 before being attenuated by gamma-Z pair production. The ionization rates of GRDRs, <~10^-16 s^-1, are much smaller than in cosmic ray dominated regions, but in the most extreme starbursts, they may still reach values comparable to those in Milky Way molecular clouds. The gas temperatures in GRDRs could be likewise low, <~10 K if there is no additional heating from dust or turbulence, while at high densities, the kinetic temperature will approach the dust temperature. The ratio of ambipolar diffusion time to free-fall time inside GRDRs in dense starbursts is expected to be similar to those in Milky Way cores, suggesting star-formation can proceed normally in them. The high columns of GRDRs may be opaque even to millimeter wavelengths, complicating direct studies of them, but I argue that they could appear as molecular line shadows in nearby starbursts with ALMA. Since GRDRs are cold, their Jeans masses are not large, so that star-formation in GRDRs may have a normal or even bottom-heavy initial mass function.

The common approach to compute the cosmic-ray distribution in an starburst galaxy or region is equivalent to assume that at any point within that environment, there is an accelerator inputing cosmic rays at a reduced rate. This rate should be compatible with the overall volume-average injection, given by the total number of accelerators that were active during the starburst age. These assumptions seem reasonable, especially under the supposition of an homogeneous and isotropic distribution of accelerators. However, in this approach the temporal evolution of the superposed spectrum is not explicitly derived; rather, it is essentially assumed ab-initio. Here, we test the validity of this approach by following the temporal evolution and spatial distribution of the superposed cosmic-ray spectrum and compare our results with those from theoretical models that treat the starburst region as a single source. In the calorimetric limit (with no cosmic-ray advection), homogeneity is reached (typically within 20%) across most of the starburst region. However, values of center-to-edge intensity ratios can amount to a factor of several. Differences between the common homogeneous assumption for the cosmic-ray distribution and our models are larger in the case of two-zone geometries, such as a central nucleus with a surrounding disc. We have also found that the decay of the cosmic-ray density following the duration of the starburst process is slow, and even approximately 1 Myr after the burst ends (for a gas density of 35 cm-3) it may still be within an order of magnitude of its peak value. Based on our simulations, it seems that the detection of a relatively hard spectrum up to the highest gamma-ray energies from nearby starburst galaxies favors a relatively small diffusion coefficient (i.e., long diffusion time) in the region where most of the emission originates.

The only invariant speed in special relativity is c; therefore, if some neutrinos travel at even tiny speeds above c, normal special relativity is incomplete and any superluminal speed may be possible. I derive a limit on superluminal neutrino speeds v >> c at high energies by noting that such speeds would increase the size of the neutrino horizon. The increased volume of the Universe visible leads to a brighter astrophysical neutrino background. The nondetection of "guaranteed" neutrino backgrounds from star-forming galaxies and ultrahigh energy cosmic rays (UHECRs) constrains v/c at TeV--ZeV energies. I find that v/c <= 820 at 60 TeV from the nondetection of neutrinos from star-forming galaxies. The nondetection of neutrinos from UHECRs constrains v/c to be less than 2500 at 0.1 EeV in a pessimistic model and less than 4.6 at 4 EeV in an optimistic model. The UHECR neutrino background nondetection is strongly inconsistent with a naive quadratic extrapolation of the OPERA results to EeV energies. The limits apply subject to some caveats, particularly that the expected pionic neutrino backgrounds exist and that neutrinos travel faster than c when they pass the detector. They could be improved substantially as the expected neutrino backgrounds are better understood and with new experimental neutrino background limits. I also point out that extremely subluminal speeds would result in a much smaller neutrino background intensity than expected.

A new class of core-collapse supernovae (SNe) has been discovered in recent years by optical/infrared surveys; these SNe suggest the presence of one or more extremely dense (~10^5-10^11 cm^-3) shells of circumstellar material (CSM) on 10^2-10^4 AU scales. We consider the collisions of the SN ejecta with these massive CSM shells as potential cosmic-ray (CR) accelerators. If ~10% of the SN energy goes into CRs, multi-TeV neutrinos and/or GeV-TeV gamma rays almost simultaneous with the optical/infrared light curves are detectable for SNe at <20-30 Mpc. A new type of coordinated multi-messenger search for such transients of duration ~1 - 10 months is required; these may give important clues to the physical origin of such SNe and to CR acceleration mechanisms.

Imaging Atmospheric Cherenkov Telescopes (IACTs) use large-aperture (~ 10 - 30 m) optical telescopes with arcminute angular resolution to detect TeV gamma-rays in the atmosphere. I show that IACTs are well-suited for optical observations of bright sources (V <= 8 - 10), because these sources are brighter than the sky background. Their advantages are especially great on rapid time-scales. Thus, IACTs are ideal for studying many phenomena optically, including transiting exoplanets and the brightest gamma-ray bursts. In principle, an IACT could achieve millimagnitude photometry of these objects with second-long exposures. I also consider the potential for optical spectroscopy with IACTs, finding that their poor angular resolution limits their usefulness for high spectral resolutions, unless complex instruments are developed. The high photon collection rate of IACTs is potentially useful for precise polarimetry. Finally, I briefly discuss the broader possibilities of extremely large, low resolution telescopes, including a 10" resolution telescope and spaceborne telescopes.

The GeV and TeV emission from M82 and NGC 253 observed by Fermi, HESS, and VERITAS constrains the physics of cosmic rays (CRs) in these dense starbursts. We argue that the gamma rays are predominantly hadronic in origin, as expected by previous studies. The measured fluxes imply that pionic losses are efficient for CR protons in both galaxies: we show that a fraction F_cal ~ 0.2 - 0.4 of the energy injected in high energy primary CR protons is lost to inelastic proton-proton collisions (pion production) before escape, producing gamma rays, neutrinos, and secondary electrons and positrons. We discuss the factor ~2 uncertainties in this estimate, including supernova rate and leptonic contributions to the GeV-TeV emission. We argue that gamma-ray data on ULIRGs like Arp 220 can test whether M82 and NGC 253 are truly calorimetric, and we present upper limits on Arp 220 from the Fermi data. We show that the observed ratio of the GeV to GHz fluxes of the starbursts suggests that non-synchrotron cooling processes are important for cooling the CR electron/positron population. We briefly reconsider previous predictions in light of the gamma-ray detections, including the starburst contribution to the gamma-ray background and CR energy densities. Finally, as a guide for future studies, we list the brightest star-forming galaxies on the sky and present updated predictions for their gamma-ray and neutrino fluxes.

We derive an upper limit to the energy of nuclei accelerated via the Fermi mechanism in any relativistic shockwave, driven by any astrophysical engine. This bound is accessible to current and upcoming ultra-high energy neutrino experiments. Detection of a single neutrino with energy above the upper limit would exclude all sites of shock acceleration, and imply physics beyond the Standard Model. We comment on the possibility that relativistic flows launched by supermassive black hole mergers are the source of the observed ultra-high energy cosmic rays.

Cosmic ray protons generate gamma-rays, neutrinos, and secondary electrons and positrons (e+/-) through pion-producing collisions with gas atoms. Any synchrotron or Inverse Compton (IC) radiation from secondary e+/- is therefore accompanied by pionic gamma-rays. Using the extragalactic gamma-ray background, we constrain the contribution of secondary e+/- to the cosmic radio, X-ray, and soft gamma-ray backgrounds. These bounds depend on the unknown hadronic contribution to the gamma-ray background and the backgrounds' source redshifts. With our assumptions, we find that IC-upscattered light from secondaries is <~ 1/4 of the MeV - GeV gamma-ray background and <~ 10% of the 0.5 keV - 1 MeV background (for sources at a redshift z <~ 10). The low intensity of the observed gamma-ray background is marginally inconsistent with a secondary e+/- origin for the radio background reported by ARCADE at ~3 GHz, unless the magnetic field strength in their sources is milliGauss or greater. These limits on the magnetic field strength are sensitive to uncertainties. However, any contribution to the gamma-ray background from sources not responsible for the ARCADE excess increases the inconsistency.

Primordial black holes (PBHs) are expected to accrete particle dark matter around them to form ultracompact minihalos (UCMHs), if the PBHs themselves are not most of the dark matter. We show that if most dark matter is a thermal relic, then the inner regions of UCMHs around PBHs are highly luminous sources of annihilation products. Flux constraints on gamma rays and neutrinos set strong abundance limits, improving previous limits by orders of magnitude. Assuming enough particle dark matter exists to form UCMHs, we find that Omega_PBH <~ 10^-4 (for m_DM c^2 ~ 100 GeV) for a vast range in PBH mass. We briefly discuss the uncertainties on our limits, including those due to the evolution of the UCMH luminosity as it annihilates.

(Abridged) The far-infrared (FIR) and radio luminosities of star-forming galaxies are linearly correlated over a very wide range in star formation rate, from normal spirals like the Milky Way to the most intense starbursts. Using one-zone models of cosmic ray (CR) injection, cooling, and escape in star-forming galaxies, we attempt to reproduce the observed FIR-radio correlation over its entire span. We show that ~2% of the kinetic energy from supernova explosions must go into primary CR electrons and that ~10 - 20% must go into primary CR protons. Secondary electrons and positrons are likely comparable to or dominate primary electrons in dense starburst galaxies. We discuss the implications of our models for the magnetic field strengths of starbursts, the detectability of starbursts by Fermi, and cosmic ray feedback. Overall, our models indicate that both CR protons and electrons escape from low surface density galaxies, but lose most of their energy before escaping dense starbursts. The FIR-radio correlation is caused by a combination of the efficient cooling of CR electrons (calorimetry) in starbursts and a conspiracy of several factors. For lower surface density galaxies, the decreasing radio emission caused by CR escape is balanced by the decreasing FIR emission caused by the low effective UV dust opacity. In starbursts, bremsstrahlung, ionization, and Inverse Compton cooling decrease the radio emission, but they are countered by secondary electrons/positrons and the decreasing critical synchrotron frequency, which both increase the radio emission. Our conclusions hold for a broad range of variations on our fiducial model.

We construct one-zone steady-state models of cosmic ray (CR) injection, cooling, and escape over the entire dynamic range of the FIR-radio correlation (FRC), from normal galaxies to starbursts, over the redshift interval 0 <= z <= 10. Normal galaxies with low star-formation rates become radio-faint at high z, because Inverse Compton (IC) losses off the CMB cool CR electrons and positrons rapidly, suppressing their nonthermal radio emission. However, we find that this effect occurs at higher redshifts than previously expected, because escape, bremsstrahlung, ionization, and starlight IC losses act to counter this effect and preserve the radio luminosity of galaxies. The radio dimming of star-forming galaxies at high z is not just a simple competition between magnetic field energy density and the CMB energy density; the CMB must also compete with every other loss process. We predict relations for the critical redshift when radio emission is significantly suppressed compared to the z ~ 0 FRC as a function of star-formation rate per unit area. Additionally, we provide a quantitative explanation for the relative radio brightness of some high-z submillimeter galaxies. We show that at fixed star formation rate surface density, galaxies with larger CR scale heights are radio bright with respect to the FRC, because of weaker bremsstrahlung and ionization losses compared to compact starbursts. We predict that these "puffy starbursts" should have steeper radio spectra than compact galaxies with the same star-formation rate surface density. We find that radio bright submillimeter galaxies alone cannot explain the excess radio emission reported by ARCADE2, but they may significantly enhance the diffuse radio background with respect to a naive application of the z ~ 0 FRC.

The Galactic disc is opaque to radio waves from extragalactic sources with frequencies nu less than ~3 MHz. However, radio waves with kHz, Hz, and even lower frequencies may propagate through the intergalactic medium (IGM). I argue that the presence of these waves can be inferred by using the Universe as our detector. I discuss possible sub-MHz sources and set new non-trivial upper limits on the energy density of sub-MHz radio waves in galaxy clusters and the average cosmic background. Limits based on five effects are considered: (1) changes in the expansion of the Universe from the radiation energy density (2) heating of the IGM by free-free absorption; (3) radiation pressure squeezing of IGM clouds by external radio waves; (4) synchrotron heating of electrons in clusters; and (5) Inverse Compton upscattering of sub-MHz radio photons. Any sub-MHz background must have an energy density much smaller than the CMB at frequencies below 1 MHz. The free-free absorption bounds from the Lyman-alpha forest are potentially the strongest, but are highly dependent on the properties of sub-MHz radio scattering in the IGM. I estimate an upper limit of 6 * 10^4 L_sun Mpc^-3 for the emissivity within Lyman-alpha forest clouds in the frequency range 5 - 200 Hz. The sub-MHz energy density in the Coma cluster is constrained to be less than ~10^-15 erg cm^-3. At present, none of the limits is strong enough to rule out a maximal T_b = 10^12 K sub-MHz synchrotron background, but other sources may be constrained with a better knowledge of sub-MHz radio propagation in the IGM.

We consider the physical conditions and origin of the z = 0.0777 absorption system observed in C III, C II, Si III, C IV, O VI, and H I absorption along the line of sight towards the quasar PHL 1811. We analysed the HST/STIS and FUSE spectra of this quasar and compared the results to Cloudy photoionization and collisional ionization models in order to derive densities, temperatures, and metallicities of the absorbing gas. The absorption can be explained by two C III clouds, offset by 35 km/s in velocity, with metallicities of ~one-tenth the solar value. One cloud has a density of order n_H = 1.2 +0.9 -0.5 * 10^-3 cm^-3 (thickness 0.4 +0.3 -0.2 kpc) and produces the observed C II and Si III absorption, while the other has a density of order n_H = 1.2 +0.9 -0.5 * 10^-5 cm^-3 (thickness 80 +70 -40 kpc) and gives rise to the observed weak C IV absorption. Cloud temperatures are ~14,000 +3000 -2000 K and ~34,000 -4000 +2000 K for photoionized models. Although collisionally ionized clouds with T ~ 70,000 K are possible, they are less likely because of the short cooling time-scales involved. Previous studies revealed no luminous galaxy at the absorber's redshift, so it is probably related to tidal debris, ejected material, a dwarf galaxy, or other halo material in a galaxy group. Our models also indicate that one of the two clouds would produce detectable weak Mg II absorption if spectral coverage of that transition existed. We predict what the system would look like at z ~ 1 when the ionizing background radiation was more intense. We find that at z ~ 1 the denser component resembles a C IV absorber. The second C III cloud in this z = 0.0777 absorber may be analogous to a subset of the more diffuse O VI absorbers at higher redshift.

Difference imaging provides a new way to discover gravitationally lensed quasars because few non-lensed sources will show spatially extended, time variable flux. We test the method on lens candidates in the Sloan Digital Sky Survey (SDSS) Supernova Survey region from the SDSS Quasar Lens Search (SQLS) and their surrounding fields. Starting from 20768 sources, including 49 SDSS quasars and 36 candidate lenses/lensed images, we find that 21 sources including 15 SDSS QSOs and 7 candidate lenses/lensed images are non-periodic variable sources. We can measure the spatial structure of the variable flux for 18 of these sources and identify only one as a non-point source. This source does not display the compelling spatial structure of the variable flux of known lensed quasars, so we reject it as a lens candidate. None of the lens candidates from the SQLS survive our cuts. Given our effective survey area of order 0.71 square degrees, this indicates a false positive rate of order one per square degree for themethod. The fraction of quasars not found to be variable and the false positive rate should both fall if we analyze the full, later data releases for the SDSS fields. While application of the method to the SDSS is limited by the resolution, depth, and sampling of the survey, several future surveys such as Pan-STARRS, LSST, and SNAP will avoid these limitations.