Here we investigate the possibility that some globular clusters (GCs) harbor intermediate mass black holes (BH) in their centers and are also embedded in a low-mass dark matter (DM) halo. Up to date, there is no evidence on whether or not GCs have DM in their constitution. For standard cold DM cosmology, it is expected that GCs form with their own DM halos. Other studies investigate the possibility that GCs were initially embedded in massive DM halos that evolved during the cluster lifetime. An additional intriguing question is related to the existence of intermediate mass black holes (IMBH) in the of GCs. The determination of whether GCs hold IMBHs would be able to answer important questions about GCs formation and the circumstances that gave rise to the IMBHs. DM & IMBH in the context of GCs are interesting subjects to be studied and we will perform such studies here, assuming the coexistence of both of them in some GCs. Having such information, we perform the study possible DM signals from GCs. One important subject to be studied is the DM density profile. In the inner regions of GCs, the DM density profile is still an open question of fundamental importance to DM studies, specially for the study of radio and $\gamma$-ray emission from DM annihilation in such regions (abridged).

Until recently, it was assumed that the diffuse extragalactic background light (EBL) in the FUV has a sharp and complete cut-off at energies above the Lyman limit. However, newly discovered extreme emission line starburst galaxies (EELGs) and other galaxies with substantial escape fractions of EUV photons in the Lyman continuum (LyC) region are increasingly common at high redshifts. The recent observations indicate that many more LyC photons at high redshifts escape from such galaxies into intergalactic space than was previously suspected. They are probably the major cause of reionization of the intergalactic medium. In this paper we estimate the contribution of these hitherto unconsidered EUV photons from EELGs at all redshifts to the (EBL), and their subsequent effect on the absorption of extragalactic high energy $\gamma$-rays. For the more distant $\gamma$-ray sources, particularly at $z \ge 3$, the intergalactic opacity above a few GeV is significantly higher than previously estimated. We calculate the results of this increased opacity on observed $\gamma$-ray spectra, which produces a high-energy turnover starting at lower energies than previously thought, and a gradual spectral steepening that may also be observable.

In this paper we explore the potential astrophysical signatures of dark matter (DM) annihilations in ultracompact dwarf galaxies (UCDs) considering two of the richest known galaxy clusters within 100 million light-years, nominally, Virgo and Fornax. Fornax UCD3 is the most luminous UCD and M59 UCD3 is the most massive UCD. With the detection of a 3.5 million solar mass black hole (BH) in Fornax UCD3, we carefully model several DM enhanced profiles scenarios, considering both the presence of the supermassive black hole (SMBH) and DM. For Fornax UCD3, the comparison of the stellar and dynamical masses suggests that there is little content of DM in UCDs. M59 UCD3 did not receive the same attention in simulations as Fornax UCD3, but deep radio imaging and X-ray observations were performed for M59 UCD3 and can be used to place limits in DM content of these UCDs. We work with an average estimative of dark matter content considering the Salpeter and Kroupa mass functions. We model Fornax UCD3 and M59 UCD3 to have a DM content that is the average of these mass functions. We then analyze the constraints for Fornax and M59 UCD3 coming from gamma-ray and radio sources considering in our simulations, a dark matter particle with mass between $10-34\,{\rm GeV}$. In the absence of strong $\gamma$-ray signatures, we show that synchrotron emission from electrons and positrons produced by DM annihilations can be very sensitive to indirect DM search. We find that DM parameters can be significantly constrained at radio frequencies and the spike profiles play an interesting rule in order to deep study the enhancements of DM \& BH interactions in ultracompact galaxies.

We give theoretical predictions for the radio emission of a dark matter candidate annihilating into 2-lepton and 4-lepton final states. We then compare our results with the known radio measurements of the sky temperature as a function of the frequency. In particular, we calculate the radio emission for some dark matter candidates annihilating into intermediate bosons that subsequently decay into a 4-lepton channel with a thermal annihilation cross-section. We show that within the range of frequencies from $20\,{\rm MHz}$ to $5\,{\rm GHz}$, this channel can produce a stronger signature than direct annihilation into leptons.

The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to observe cosmic neutrinos (CNs) above 20 PeV and ultra-high energy cosmic rays (UHECRs) above 20 EeV over the full sky. The POEMMA mission calls for two identical satellites flying in loose formation, each comprised of a 4-meter wide field-of-view (45 degrees) Schmidt photometer. The hybrid focal surface includes a fast (1 ${\mu}$s) ultraviolet camera for fluorescence observations and an ultrafast (10 ns) optical camera for Cherenkov observations. POEMMA will provide new multi-messenger windows onto the most energetic events in the universe, enabling the study of new astrophysics and particle physics at these otherwise inaccessible energies.

The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is a NASA Astrophysics probe-class mission designed to observe ultra-high energy cosmic rays (UHECRs) and cosmic neutrinos from space. Astro2020 APC white paper: Medium-class Space Particle Astrophysics Project.

Recent radio surveys have discovered a large number of low luminosity core dominated radio galaxies that are much more abundant than those at higher luminosities. These objects will be too faint in gamma-rays to be detected individually by Fermi. Nevertheless, they may contribute significantly to the unresolved extragalactic gamma-ray background. We consider here the possible contribution of these core dominated radio galaxies to the diffuse extragalactic gamma-ray background. Using published data available for all 45 of the radiogalaxies listed as detected counterparts in the Fermi FL8Y source list update to the 3FGL catalog, we have searched for radio maps which can resolve the core flux from the total source flux. Using high resolution radio maps we were able to obtain core fluxes for virtually every source. We then derived a relation between core radio flux and gamma-ray flux that we extrapolated to sources with low radio luminosities that are known to be highly core dominated. We then employed a very recent determination of the luminosity function for core dominated radio galaxies in order to obtain the contribution of all possible gamma-ray emitting radio galaxies to the unresolved extragalactic gamma-ray background. We find this contribution to be a possibly non-negligible, 4% - 18% of the background.

Observations of high energy neutrinos, both in the laboratory and from cosmic sources, can be a useful probe in searching for new physics. Such observations can provide sensitive tests of Lorentz invariance violation (LIV), which may be a the result of quantum gravity physics (QG). We review some observationally testable consequences of LIV using effective field theory (EFT) formalism. To do this, one can postulate the existence of additional small LIV terms in free particle Lagrangians, suppressed by powers of the Planck mass. The observational consequences of such terms are then examined. In particular, one can place limits on a class of non-renormalizable, mass dimension five and six Lorentz invariance violating operators that may be the result of QG.

We consider here three dark matter models with the gauge symmetry of the standard model plus an additional local $U(1)_D$ factor. One model is truly secluded and the other two models begin flipped, but end up secluded. All of these models include one dark fermion and one vector boson that gains mass via the Stueckelberg mechanism. We show that the would be flipped models provide an example dark matter composed of "almost least interacting particles" (ALIPs). Such particles are therefore compatible with the constraints obtained from both laboratory measurements and astrophysical observations.

We discuss some of the tests of Lorentz symmetry made possible by astrophysical observations of ultrahigh energy cosmic rays, gamma-rays, and neutrinos. These are among the most sensitive tests of Lorentz symmetry violation because they are the highest energy phenomena known to man.

The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) mission is being designed to establish charged-particle astronomy with ultra-high energy cosmic rays (UHECRs) and to observe cosmogenic tau neutrinos (CTNs). The study of UHECRs and CTNs from space will yield orders-of-magnitude increase in statistics of observed UHECRs at the highest energies, and the observation of the cosmogenic flux of neutrinos for a range of UHECR models. These observations should solve the long-standing puzzle of the origin of the highest energy particles ever observed, providing a new window onto the most energetic environments and events in the Universe, while studying particle interactions well beyond accelerator energies. The discovery of CTNs will help solve the puzzle of the origin of UHECRs and begin a new field of Astroparticle Physics with the study of neutrino properties at ultra-high energies.

We have previously calculated the intergalactic background light (IBL) as a function of redshift in the far ultraviolet to near infrared range, based purely on data from deep galaxy surveys. Here we utilize similar methods to determine the mid- and far infrared IBL out to a wavelength of 850 microns. Our approach enables us to constrain the range of photon densities, based on the uncertainties from observationally determined luminosity densities and colors. By also including the effect of the 2.7 K cosmic background photons, we determine 68% confidence upper and lower limits on the opacity of the universe to gamma-rays up to PeV energies. Our direct results on the IBL are consistent with those from complimentary gamma-ray analyses using observations from the Fermi $\gamma$-ray space telescope and the H.E.S.S. air Cherenkov telescope. Thus, we find no evidence of previously suggested processes for the modification of gamma-ray spectra other than that of absorption by pair production alone.

The discovery of extraterrestrial very-high-energy neutrinos by the IceCube collaboration has launched a quest for the identification of their astrophysical sources. Gamma-ray blazars have been predicted to yield a cumulative neutrino signal exceeding the atmospheric background above energies of 100 TeV, assuming that both the neutrinos and the gamma-ray photons are produced by accelerated protons in relativistic jets. Since the background spectrum falls steeply with increasing energy, the individual events with the clearest signature of being of an extraterrestrial origin are those at PeV energies. Inside the large positional-uncertainty fields of the first two PeV neutrinos detected by IceCube, the integrated emission of the blazar population has a sufficiently high electromagnetic flux to explain the detected IceCube events, but fluences of individual objects are too low to make an unambiguous source association. Here, we report that a major outburst of the blazar PKS B1424-418 occurred in temporal and positional coincidence with the third PeV-energy neutrino event (IC35) detected by IceCube. Based on an analysis of the full sample of gamma-ray blazars in the IC35 field and assuming a photo-hadronic emission model, we show that the long-term average gamma-ray emission of blazars as a class is in agreement with both the measured all-sky flux of PeV neutrinos and the spectral slope of the IceCube signal. The outburst of PKS B1424-418 has provided an energy output high enough to explain the observed PeV event, indicative of a direct physical association.

We propose to investigate a secluded WIMP dark matter model consisting of neutral fermions as the dark matter candidate and a Proca-Wentzel (PW) field as a mediator. In the model that we consider here, dark matter WIMPs interact with standard model (SM) particles only through the PW field of ~ MeV -- multi-GeV mass particles. The interactions occur via an U(1)' mediator, V_{\mu}', which couples to the SM by kinetic mixing with U(1) hypercharge bosons, B_{\mu}. One important difference between our model and other such models in the literature is the absence of an extra singlet scalar, so that the parameter with dimension of mass M^2_V is not related to a spontaneous symmetry breaking. This QED based model is also renormalizable. The mass scale of the mediator and the absence of the singlet scalar can lead to interesting astrophysical signatures. The dominant annihilation channels are different from those usually considered in previous work. We show that the GeV-energy gamma-ray excess in the galactic center region, as derived from Fermi-LAT Gamma-ray Space Telescope data, can be attributed to such secluded dark matter WIMPs, given parameters of the model that are consistent with the cosmological dark matter density. Secluded WIMP models are also consistent with suggested upper limits on the DM contribution to the cosmic-ray antiproton flux.

High energy cosmic neutrino observations provide a sensitive test of Lorentz invariance violation, which may be a consequence of quantum gravity theories. We consider a class of non-renormalizable, Lorentz invariance violating operators that arise in an effective field theory description of Lorentz invariance violation in the neutrino sector inspired by Planck-scale physics and quantum gravity models. We assume a conservative generic scenario for the redshift distribution of extragalactic neutrino sources and employ Monte Carlo techniques to describe superluminal neutrino propagation, treating kinematically allowed energy losses of superluminal neutrinos caused by both vacuum pair emission and neutrino splitting. We consider EFTs with both non-renormalizable CPT-odd and non-renormalizable CPT-even operator dominance. We then compare the spectra derived using our Monte Carlo calculations in both cases with the spectrum observed by IceCube in order to determine the implications of our results regarding Planck-scale physics. We find that if the drop off in the neutrino flux above ~2 PeV is caused by Planck scale physics, rather than by a limiting energy in the source emission, a potentially significant pileup effect would be produced just below the drop off energy in the case of CPT-even operator dominance. However, such a clear drop off effect would not be observed if the CPT-odd, CPT-violating term dominates.

The IceCube observation of cosmic neutrinos with $E_{\nu} > 60$ TeV, most of which are likely of extragalactic origin, allows one to severely constrain Lorentz invariance violation (LIV) in the neutrino sector, allowing for the possible existence of superluminal neutrinos. The subsequent neutrino energy loss by vacuum $e^+e^-$ pair emission (VPE) is strongly dependent on the strength of LIV. In this paper we explore the physics and cosmology of superluminal neutrino propagation. We consider a conservative scenario for the redshift distribution of neutrino sources. Then by propagating a generic neutrino spectrum, using Monte Carlo techniques to take account of energy losses from both VPE and redshifting, we obtain the best present constraints on LIV parameters involving neutrinos. We find that $\delta_{\nu e} = \delta_{\nu} - \delta_e \le 5.2 \times 10^{-21}$. Taking $\delta_e \le 5 \times 10^{-21}$, we then obtain an upper limit on the superluminal velocity fraction for neutrinos alone of $1.0 \times 10^{-20}$. Interestingly, by taking $\delta_{\nu e} = 5.2 \times 10^{-21}$, we obtain a cutoff in the predicted neutrino spectrum above 2 PeV that is consistent with the lack of observed neutrinos at those energies, and particularly at the Glashow resonance energy of 6.3 PeV. Thus, such a cutoff could be the result of neutrinos being slightly superluminal, with $\delta_{\nu}$ being $(0.5 \ {\rm to} \ 1.0) \times 10^{-20}$.

We describe the science motivation and development of a pair production telescope for medium-energy gamma-ray polarimetry. Our instrument concept, the Advanced Energetic Pair Telescope (AdEPT), takes advantage of the Three-Dimensional Track Imager, a low-density gaseous time projection chamber, to achieve angular resolution within a factor of two of the pair production kinematics limit (~0.6 deg at 70 MeV), continuum sensitivity comparable with the Fermi-LAT front detector (<3x10-6 MeV cm-2 s-1 at 70 MeV), and minimum detectable polarization less than 10% for a 10 millicrab source in 106 seconds.

We extend our previous model-independent determination of the intergalactic background light (IBL), based purely on galaxy survey data, out to a wavelength of 5 microns. Our approach enables us to constrain the range of photon densities, based on the uncertainties from observationally determined luminosity densities and colors. We further determine a 68% confidence upper and lower limit on the opacity of the universe to gamma-rays up to energies of 1.6/(1+z) TeV. A comparison of our lower limit redshift-dependent opacity curves to the opacity limits derived from the results of both ground-based air Cherenkov telescope and Fermi-LAT observations of PKS 1424+240 allows us to place a new upper limit on the redshift of this source, independent of IBL modeling.

The observation of two PeV-scale neutrino events reported by Ice Cube can, in principle, allows one to place constraints on Lorentz invariance violation (LIV) in the neutrino sector. After first arguing that at least one of the IceCube events was of extragalactic origin, I derive an upper limit for {\it the difference} between putative superluminal neutrino and electron velocities of $\le \sim 5.6 \times 10^{-19}$ in units where $c = 1$, confirming that the observed PeV neutrinos could have reached Earth from extragalactic sources. I further derive a new constraint on the superluminal electron velocity, obtained from the observation of synchrotron radiation in the Crab Nebula flare of September, 2010. The inference that the $>$ 1 GeV $\gamma$-rays from synchrotron emission in the flare were produced by electrons of energy up to $\sim 5.1$ PeV indicates the non-occurrence of vacuum \'{C}erenkov radiation by these electrons. This implies a new, strong constraint on superluminal electron velocities $\delta_e \le \sim 5 \times 10^{-21}$. It immediately follows that one then obtains an upper limit on the superluminal neutrino velocity {\it alone} of $\delta_{\nu} \le \sim 5.6 \times 10^{-19}$, many orders of magnitude better than the time-of-flight constraint from the SN1987A neutrino burst. However, if the electrons are {\it subluminal} the constraint on $|\delta_e| \le \sim 8 \times 10^{-17}$, obtained from the Crab Nebula $\gamma$-ray spectrum, places a weaker constraint on superluminal neutrino velocity of $\delta_{\nu} \le \sim 8 \times 10^{-17}$.

Report of the CF6 Working Group at Snowmass 2013. Topics addressed include ultra-high energy cosmic rays, neutrinos, gamma rays, baryogenesis, and experiments probing the fundamental nature of spacetime.

Some Quantum Gravity (QG) theories allow for a violation of Lorentz invariance (LIV), manifesting as a dependence of the velocity of light in vacuum on its energy. If such a dependence exists, then photons of different energies emitted together by a distant source will arrive at the Earth at different times. High-energy (GeV) transient emissions from distant astrophysical sources such as Gamma-ray Bursts (GRBs) and Active Galaxy Nuclei can be used to search for and constrain LIV. The Fermi collaboration has previously analyzed two GRBs in order to put constraints on the dispersion parameter in vacuum, and on the energy scale at which QG effects causing LIV may arise. We used three different methods on four bright GRBs observed by the Fermi-LAT to get more stringent and robust constraints. No delays have been detected and strong limits on the QG energy scale are derived: for linear dispersion we set tight constraints placing the QG energy scale above the Planck mass; a quadratic leading LIV effect is also constrained.

I show that the high energy neutrino flux predicted to arise from AGN cores can explain the PeV neutrinos detected by Ice Cube without conflicting with the constraints from the observed extragalactic cosmic ray and gamma-ray backgrounds.

High-energy astrophysics observations provide the best possibilities to detect a very small violation of Lorentz invariance, such as may be related to the structure of space-time near the Planck scale. I discuss the possible signatures of Lorentz invariance violation that can be manifested by observing the spectra, polarization, and timing of gamma-rays from active galactic nuclei and gamma-ray bursts. Other sensitive tests are provided by observations of the spectra of ultrahigh-energy cosmic rays and very high-energy neutrinos. I also discuss a new time-of-flight analysis of observations of GRB 090510 by the Fermi gamma-ray Space Telescope. These results, based on high-energy astrophysical observations, have fundamental implications for space-time physics and quantum gravity models.

We analyze the MeV/GeV emission from four bright Gamma-Ray Bursts (GRBs) observed by the Fermi-Large Area Telescope to produce robust, stringent constraints on a dependence of the speed of light in vacuo on the photon energy (vacuum dispersion), a form of Lorentz invariance violation (LIV) allowed by some Quantum Gravity (QG) theories. First, we use three different and complementary techniques to constrain the total degree of dispersion observed in the data. Additionally, using a maximally conservative set of assumptions on possible source-intrinsic spectral-evolution effects, we constrain any vacuum dispersion solely attributed to LIV. We then derive limits on the "QG energy scale" (the energy scale that LIV-inducing QG effects become important, E_QG) and the coefficients of the Standard Model Extension. For the subluminal case (where high energy photons propagate more slowly than lower energy photons) and without taking into account any source-intrinsic dispersion, our most stringent limits (at 95% CL) are obtained from GRB090510 and are E_{QG,1}>7.6 times the Planck energy (E_Pl) and E_{QG,2}>1.3 x 10^11 GeV for linear and quadratic leading order LIV-induced vacuum dispersion, respectively. These limits improve the latest constraints by Fermi and H.E.S.S. by a factor of ~2. Our results disfavor any class of models requiring E_{QG,1} \lesssim E_Pl.

We discuss the prospects of doing tests of Lorentz invariance with gamma-rays observed with present and future ground based gamma-ray observatories.

Secondary nuclear production physics is receiving increased attention given the high-quality measurements of the gamma-ray emissivity of local interstellar gas between ~50 MeV and ~40 GeV, obtained with the Large Area Telescope on board the Fermi space observatory. More than 90% of the gas-related emissivity above 1 GeV is attributed to gamma-rays from the decay of neutral pions formed in collisions between cosmic rays and interstellar matter, with lepton-induced processes becoming increasingly important below 1 GeV. The elementary kinematics of neutral pion production and decay are re-examined in light of two physics questions: does isobaric production follow a scaling behavior? and what is the minimum proton kinetic energy needed to make a gamma-ray of a certain energy formed through intermediate pi0 production? The emissivity spectrum will allow the interstellar cosmic-ray spectrum to be determined reliably, providing a reference for origin and propagation studies as well as input to solar modulation models. A method for such an analysis and illustrative results are presented.

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

I present the results of a new approach to the intensity and photon density spectrum of the intergalactic background light as a function of redshift using observational data obtained in many different wavelength bands from local to deep galaxy surveys. This enables an empirical determination of both the EBL and its observationally based uncertainties. Using these results one can place 68% confidence upper and lower limits on the opacity of the universe to {\gamma}-rays, free of the theoretical assumptions that were needed for past calculations. I compare our results with measurements of the extragalactic background light, upper limits obtained from observations made by the Fermi Gamma-ray Space Telescope, and new observationally based results from Fermi and H.E.S.S. using recent analyses of blazar spectra.

We calculate the intensity and photon spectrum of the intergalactic background light (IBL) as a function of redshift using an approach based on observational data obtained in many different wavelength bands from local to deep galaxy surveys. This allows us to obtain an empirical determination of the IBL and to quantify its observationally based uncertainties. Using our results on the IBL, we then place 68% confidence upper and lower limits on the opacity of the universe to gamma-rays, free of the theoretical assumptions that were needed for past calculations. We compare our results with measurements of the extragalactic background light and upper limits obtained from observations made by the Fermi Gamma-ray Space Telescope.

Constraints on possible Lorentz invariance violation (LIV) to first order in $E/M_{\rm Planck}$ for photons in the framework of effective field theory (EFT) are discussed, taking cosmological factors into account. Then, using the reported detection of polarized soft $\gamma$-ray emission from the $\gamma$-ray burst GRB041219a that is indicative of an absence of vacuum birefringence, together with a very recent improved method for estimating the redshift of the burst, we derive constraints on the dimension 5 Lorentz violating modification to the Lagrangian of an effective local QFT for QED. Our new constraints are more than five orders of magnitude better than recent constraints from observations of the Crab Nebula. We obtain the upper limit on the Lorentz violating dimension 5 EFT parameter $|\xi|$ of $2.4 \times 10^{-15}$, corresponding to a constraint on the dimension 5 standard model extension parameter, $k^{(5)}_{(V)00} \le 4.2 \times 10^{-34}$ GeV$^{-1}$.

We present new theoretical estimates of the relative contributions of unresolved blazars and star-forming galaxies to the extragalactic gamma-ray background (EGB) and discuss constraints on the contributions from alternative mechanisms such as dark matter annihilation and truly diffuse gamma-ray production. We find that the Fermi source count data do not rule out a scenario in which the EGB is dominated by emission from unresolved blazars, though unresolved star-forming galaxies may also contribute significantly to the background, within order-of-magnitude uncertainties. In addition, we find that the spectrum of the unresolved star-forming galaxy contribution cannot explain the EGB spectrum found by EGRET at energies between 50 and 200 MeV, whereas the spectrum of unresolved FSRQs, when accounting for the energy-dependent effects of source confusion, could be consistent with the combined spectrum of the low-energy EGRET EGB measurements and the Fermi-LAT EGB measurements.

We have previously shown that a very small amount of Lorentz invariance violation (LIV), which suppresses photomeson interactions of ultrahigh energy cosmic rays (UHECRs) with cosmic background radiation (CBR) photons, can produce a spectrum of cosmic rays that is consistent with that currently observed by the Pierre Auger Observatory (PAO) and HiRes experiments. Here, we calculate the corresponding flux of high energy neutrinos generated by the propagation of UHECR protons through the CBR in the presence of LIV. We find that LIV produces a reduction in the flux of the highest energy neutrinos and a reduction in the energy of the peak of the neutrino energy flux spectrum, both depending on the strength of the LIV. Thus, observations of the UHE neutrino spectrum provide a clear test for the existence and amount of LIV at the highest energies. We further discuss the ability of current and future proposed detectors make such observations.

We derive a relation for the steepening of blazar gamma-ray spectra between the multi-GeV Fermi energy range and the TeV energy range observed by atmospheric Cerenkov telescopes. The change in spectral index is produced by two effects: (1) an intrinsic steepening, independent of redshift, owing to the properties of emission and absorption in the source, and (2) a redshift-dependent steepening produced by intergalactic pair production interactions of blazar gamma-rays with low energy photons of the intergalactic background light (IBL). Given this relation, with good enough data on the mean gamma-ray SED of TeV selected BL Lacs, the redshift evolution of the IBL can, in principle, be determined independently of stellar evolution models. We apply our relation to the results of new Fermi observations of TeV selected blazars.

The topic of Lorentz invariance violation is a fundamental question in physics that has taken on particular interest in theoretical explorations of quantum gravity scenarios. I discuss various gamma-ray observations that give limits on predicted potential effects of Lorentz invariance violation. Among these are spectral data from ground based observations of the multi-TeV gamma-rays from nearby AGN, INTEGRAL detections of polarized soft gamma-rays from the vicinity of the Crab pulsar, Fermi Gamma Ray Space Telescope studies of photon propagation timing from gamma-ray bursts, and Auger data on the spectrum of ultrahigh energy cosmic rays. These results can be used to seriously constrain or rule out some models involving Planck scale physics. Possible implications of these limits for quantum gravity and Planck scale physics will be discussed.

Ultrahigh energy cosmic rays that produce giant extensive showers of charged particles and photons when they interact in the Earth's atmosphere provide a unique tool to search for new physics. Of particular interest is the possibility of detecting a very small violation of Lorentz invariance such as may be related to the structure of space-time near the Planck scale of $\sim 10^{-35}$m. We discuss here the possible signature of Lorentz invariance violation on the spectrum of ultrahigh energy cosmic rays as compared with present observations of giant air showers. We also discuss the possibilities of using more sensitive detection techniques to improve searches for Lorentz invariance violation in the future. Using the latest data from the Pierre Auger Observatory, we derive a best fit to the LIV parameter of $3.0^{+1.5}_{-3.0} \times 10^{-23}$, corresponding to an upper limit of $4.5 \times 10^{-23}$ at a proton Lorentz factor of $\sim 2 \times 10^{11}$. This result has fundamental implications for quantum gravity models.

A fundamental question that can be answered in the next decade is: WHAT IS THE ORIGIN OF THE HIGHEST ENERGY COSMIC PARTICLES? The discovery of the sources of the highest energy cosmic rays will reveal the workings of the most energetic astrophysical environments in the recent universe. Candidate sources range from the birth of compact objects to explosions related to gamma-ray bursts or generated around supermassive black holes in active galactic nuclei. In addition to beginning a new era of high-energy astrophysics, the study of ultra-high energy cosmic rays will constrain the structure of the Galactic and extragalactic magnetic fields. The propagation of these particles from source to Earth also probes the cosmic background radiation and gives insight into particle interactions at orders of magnitude higher energy than can be achieved in terrestrial laboratories. Next generation observatories designed to study the highest energy cosmic rays will have unprecedented sensitivity to ultra-high energy photons and neutrinos, which will further illuminate the workings of the universe at the most extreme energies. For this challenge to be met during the 2010-2020 decade, a significant increase in the integrated exposure to cosmic rays above 6 1019 eV will be necessary. The technical capabilities for answering this open question are at hand and the time is ripe for exploring Charged Particle Astronomy.

There has been much interest in possible violations of Lorentz invariance, particularly motivated by quantum gravity theories. It has been suggested that a small amount of Lorentz invariance violation (LIV) could turn off photomeson interactions of ultrahigh energy cosmic rays (UHECRs) with photons of the cosmic background radiation and thereby eliminate the resulting sharp steepening in the spectrum of the highest energy CRs predicted by Greisen Zatsepin and Kuzmin (GZK). Recent measurements of the UHECR spectrum reported by the HiRes and Auger collaborations, however, indicate the presence of the GZK effect. We present the results of a detailed calculation of the modification of the UHECR spectrum caused by LIV using the formalism of Coleman and Glashow. We then compare these results with the experimental UHECR data from Auger and HiRes. Based on these data, we find a best fit amount of LIV of $4.5^{+1.5}_{-4.5} \times 10^{-23}$,consistent with an upper limit of $6 \times 10^{-23}$. This possible amount of LIV can lead to a recovery of the cosmic ray spectrum at higher energies than presently observed. Such an LIV recovery effect can be tested observationally using future detectors.

The MAGIC collaboration has recently reported the detection of the strong gamma-ray blazar 3C279 during a 1-2 day flare. They have used their spectral observations to draw conclusions regarding upper limits on the opacity of the Universe to high energy gamma-rays and, by implication, upper limits on the extragalactic mid-infrared background radiation. In this paper we examine the effect of gamma-ray absorption by the extragalactic infrared radiation on intrinsic spectra for this blazar and compare our results with the observational data on 3C279. We find agreement with our previous results, contrary to the recent assertion of the MAGIC group that the Universe is more transparent to \gray s than our calculations indicate. Our analysis indicates that in the energy range between ~80 and ~500 GeV, 3C279 has a best-fit intrinsic spectrum with a spectral index ~1.78 using our fast evolution model and ~2.19 using our baseline model. However, we also find that spectral indices in the range of 0.0 to 3.0 are almost as equally acceptable as the best fit spectral indices. Assuming the same intrinsic spectral index for this flare as for the 1991 flare from 3C279 observed by EGRET, viz., 2.02, which lies between our best fit indices, we estimate that the MAGIC flare was ~3 times brighter than the EGRET flare observed 15 years earlier.

Analysis of data from the EGRET gamma-ray detector on the Compton Gamma Ray Observatory indicated an anomaly in the form of an excess diffuse galactic flux at GeV energies over that which was theoretically predicted. Various explanations for this anomaly have been put forth, including the invocation of supersymmetric dark matter annihilation. We reexamine these explanations here, including a new discussion of the possible systematic errors in the sensitivity determination of the EGRET detector. We conclude that the most likely explanation of the EGRET ``GeV anomaly'' was an error in the estimation of the of the EGRET sensitivity at energies above ~1 GeV. We give reasons why such a situation could have occurred. We find evidence from our new all-sky analysis which is inconsistent with the assumption that the anomaly can be a signal of supersymmetric dark matter annihilation. We also reconfirm the original results of the EGRET team on the extragalactic gamma-ray background spectrum. There are important implications of our analysis for the upcoming Gamma Ray Large Area Telescope (GLAST) mission.

We examine the effect of gamma-ray absorption by the extragalactic infrared radiation on intrinsic spectra predicted for 1ES0229+200 and compare our results with the observational data. We find agreement with our previous results on the shape of the IR spectral energy distribution (SED), contrary to the recent assertion of the HESS group. Our analysis indicates that 1ES0229+200 has a very hard intrinsic spectrum with a spectral index between 1.1 +/- 0.3 and 1.5 +/- 0.3 in the energy range between ~0.5 TeV and ~15 TeV. Under the assumptions that (1) the SED models of Stecker, Malkan & Scully (2006) are reasonable as derived from numerous detailed IR observations, and (2) spectral indexes in the range 1 < \Gamma < 1.5 have been shown to be obtainable from relativistic shock acceleration under the astrophysical conditions extant in blazar flares (Stecker, Baring & Summerlin 2007), the fits to the observations of 1ES0229+200 using our previous IR SEDs are consistent with both the IR and gamma-ray observations. Our analysis presents evidence indicating that the energy spectrum of relativistic particles in 1ES0229+200 is produced by relativistic shock acceleration, producing an intrinsic gamma-ray spectrum with index 1 < \Gamma < 1.5 and with no evidence of a peak in the SED up to energies ~15 TeV.

Ongoing deep surveys of galaxy luminosity functions, spectral energy distributions and backwards evolution models of star formation rates can be used to calculate the past history of intergalactic photon densities and, from them, the present and past optical depth of the universe. This procedure can be reversed by looking for sharp cutoffs in the spectra of extragalactic gamma-ray sources at high redshifts in the multi-GeV energy range with GLAST. Determining the cutoff energies of sources with known redshifts and little intrinsic absorption can enable a more precise determination of the past intergalactic photon density and thus allow a better determination of the past history of the total star formation rate, including that from galaxies too faint to be observed. Conversely, observations of sharp high energy cutoffs in the gamma-ray spectra of blazars at unknown redshifts can be used instead of spectral lines to give a measure of their redshifts. Also, given a knowledge photon densities, one can derive the intrinsic gamma-ray spectra and luminosities of blazars over a range of redshifts and look for possible trends in blazar evolution. I present some evidence hinting that TeV blazars with flatter spectra have higher intrinsic TeV gamma-ray luminosities and that there may be a correlation of flatness and luminosity with redshift. (abridged)

The observed spectra of blazars, their intrinsic emission, and the underlying populations of radiating particles are intimately related. The use of these sources as probes of the extragalactic infrared background, a prospect propelled by recent advances in TeV-band telescopes, soon to be augmented by observations by NASA's upcoming Gamma-Ray Large Area Space Telescope (GLAST), has been a topic of great recent interest. Here, it is demonstrated that if particles in blazar jets are accelerated at relativistic shocks, then gamma-ray spectra with indices less than 1.5 can be produced. This, in turn, loosens the upper limits on the near infrared extragalactic background radiation previously proposed. We also show evidence hinting that TeV blazars with flatter spectra have higher intrinsic TeV gamma-ray luminosities and we indicate that there may be a correlation of flatness and luminosity with redshift.

Stecker, Malkan and Scully, have shown how ongoing deep surveys of galaxy luminosity functions, spectral energy distributions and backwards evolution models of star formation rates can be used to calculate the past history of intergalactic photon densities for energies from 0.03 eV to the Lyman limit at 13.6 eV and for redshifts out to 6 (called here the intergalactic background light or IBL). From these calculations of the IBL at various redshifts, they predict the present and past optical depth of the universe to high energy gamma-rays owing to interactions with photons of the IBL and the 2.7 K CMB. We discuss here how this proceedure can be reversed by looking for sharp cutoffs in the spectra of extragalactic gamma-ray sources such as blazars at high redshifts in the multi-GeV energy range with GLAST. By determining the cutoff energies of sources with known redshifts, we can refine our determination of the IBL photon densities in the past, i.e., the "archeo-IBL", and therefore get a better measure of the past history of the total star formation rate. Conversely, observations of sharp high energy cutoffs in the gamma-ray spectra of sources at unknown redshifts can be used instead of spectral lines to give a measure of their redshifts.

Table 1 in our paper, ApJ 648, 774 (2006) entitled "Intergalactic Photon Spectra from the Far IR to the UV Lyman Limit for 0 < z < 6 and the Optical Depth of the Universe to High Energy Gamma-Rays" had erroneous numbers for the coefficients fitting the parametric form for the optical depth of the universe to gamma-rays. The correct values for these parameters as described in the original text are given here in a corrected table for various redshifts for the baseline model (upper row) and fast evolution (lower row) for each individual redshift. The parametric approximation is good for optical depths between 0.01 and 100 and for gamma-ray energies up to ~2 TeV for all redshifts but also for energies up to ~10 TeV for redshifts less than 1.

The determination of the densities of intergalactic photons from the FIR to the UV produced by stellar emission and dust reradiation at various redshifts can provide an independent measure of the star formation history of the universe. High energy gamma-rays can annihilate with FIR-UV photons to produce electron-positron pairs which result in high-end absorption cutoffs in the gamma-ray spectra of extragalactic sources. Future measurements of such absorption in the spectra of extragalactic high energy gamma-ray sources at higher redshifts from detectors such as the (soon to be launched) GLAST space telescope can be used to determine intergalactic photon densities in the distant past, thereby shedding light on the history of star formation and galaxy evolution. (abridged),

We derive a new and user friendly simple analytic approximation for determining the effect of intergalactic absorption in the energy range 0.2-2 TeV and the redshift range 0.05-0.4. In these ranges, the form of the absorption coeeficient is approximately logarithmic in energy. The effect of this energy dependence is to steepen intrinsic source spectra such that a source with an approximate power-law intrinsic spectrum in this energy range with spectral index $\Gamma_{s}$ is steepened to a power-law with an observed spectral index $\Gamma_{o} = $\Gamma_{s} + $\Delta \Gamma (z)$ where $\Delta \Gamma (z)$ is a linear function of z in the redshift range 0.05-0.4. We apply this approximation to the spectra of seven TeV blazars.

Loeb and Waxman have argued that high energy neutrinos from the decay of pions produced in interactions of cosmic rays with interstellar gas in starburst galaxies would be produced with a large enough flux to be observable. Here we obtain an upper limit to the diffuse neutrino flux from starburst galaxies which is a factor of $\sim$5 lower than the flux which they predict. Compared with predicted fluxes from other extragalactic high energy neutrino sources, starburst neutrinos with $\sim$ PeV energies would have a flux considerably below that predicted for AGN models. We also estimate an upper limit for the diffuse GeV $\gamma$-ray flux from starbust galaxies to be $\cal{O}$$(10^{-2})$ of the observed $\gamma$-ray background, much less than the component from unresolved blazars.

Loeb and Waxman have argued that high energy neutrinos from the decay of pions produced in interactions of cosmic rays with interstellar gas in starburst galaxies would be produced with a large enough flux to be observable. Their model is reexamined here and we obtain an upper limit to the diffuse neutrino flux from starburst galaxies. The upper limit obtained here is a factor of ~5 lower than the flux which they predict. Our predicted neutrino flux would be below the atmospheric neutrino foreground flux at energies below \~300 TeV and therefore would be unobservable. PeV neutrinos from starburst galaxies are also unlikely to be detected. Compared with predicted fluxes from other extragalactic high energy neutrino sources, starburst neutrinos with ~PeV energies would have a flux considerably below that predicted for AGN models. We also estimate an upper limit for the diffuse GeV gamma-ray flux from starburst galaxies to be about two orders of magnitude below the observed gamma-ray background, much less than the background from blazars and more than an order of magnitude below that calculated by Thompson et al.

(ABRIDGED) The Gamma-ray Large Area Space Telescope (GLAST) will measure the spectra of distant extragalactic sources of high energy gamma-rays. GLAST can look for energy dependent propagation effects from such sources as a signal of Lorentz invariance violation (LIV). Such sources should also exhibit high energy spectral cutoffs from pair production interactions with low energy photons. The properties of such cutoffs can also be used to test LIV. Detectors to measure gamma-ray polarization can look for the depolarizing effect of space-time birefingence predicted by loop quantum gravity. A spaceborne detector array looking down on Earth to study extensive air showers produced by ultrahigh energy cosmic rays can study their spectral properties and look for a possible deviation from the predicted GZK effect as another signal of LIV.

We calculate the intergalactic photon density as a function of both energy and redshift for 0 < z < 6 for photon energies from .003 eV to the Lyman limit cutoff at 13.6 eV in a Lambda-CDM universe with $\Omega_{\Lambda} = 0.7$ and $\Omega_{m} = 0.3$. Our galaxy evolution model gives results which are consistent with Spitzer deep number counts and the spectral energy distribution of the extragalactic background radiation. We use our photon density results to extend previous work on the absorption of high energy gamma-rays in intergalactic space owing to interactions with low energy photons and the 2.7 K cosmic background radiation. We calculate the optical depth of the universe, tau, for gamma-rays having energies from 4 GeV to 100 TeV emitted by sources at redshifts from ~0 to 5. We also give an analytic fit with numerical coefficients for approximating $\tau(E_{\gamma}, z)$. As an example of the application of our results, we calculate the absorbed spectrum of the blazar PKS 2155-304 at z = 0.117 and compare it with the spectrum observed by the H.E.S.S. air Cherenkov gamma-ray telescope array.

We calculate the flux of neutrinos generated by the propagation of ultra-high energy nuclei over cosmological distances. The propagation takes into account the interactions with cosmic background radiations including the CMB and the most recent estimates of higher energy (infra-red, optical, and ultra violet) backgrounds. We assume that the composition of ultra-high energy cosmic rays (UHECRs) at the source is the same as the observed one at low energies. This assumption fits well the present data at the highest energies. We compare the cosmogenic neutrino flux from mixed composition sources to pure proton sources. We find that the neutrino flux in the mixed composition case has a high energy peak, mainly due to photopion production off CMB photons, of similar shape and amplitude to the proton case. At low energies both composition cases have significant neutrino flux with a peak around 1014.5 eV due to the higher energy backgrounds. The mixed composition case induces a higher flux of neutrinos at energies below 1013 eV due to the neutron decay component that extends down to low energies. Detection of diffuse neutrino fluxes at ultra high energies can strongly constrain the source distribution of UHECR whereas fluxes at lower energies could be used to constrain confinement of VHE and UHE cosmic rays if combined with composition analysis from cosmic ray experiments.

We discuss the production of cosmogenic neutrinos on extragalactic infrared photons in a model of its cosmological evolution. The relative importance of these infrared photons as a target for proton interactions is significant, especially in the case of steep injection spectra of the ultrahigh energy cosmic rays. For an E$^{-2.5}$ cosmic ray injection spectrum, for example, the event rate of neutrinos of energy above 1 PeV is more than doubled.

Taking into account new physics and astronomy developments I give a revised high energy neutrino flux for the AGN core model of Stecker, Done, Salamon and Sommers

Ultrahigh-energy cosmic rays passing through the central region of the Galaxy interact with starlight and the infrared photons. Both nuclei and protons generate secondary fluxes of photons and neutrinos on their passage through the central region. We compute the fluxes of these secondary particles, the observations of which can be used to improve one's understanding of origin and composition of ultrahigh-energy comic rays, especially if the violation of the Greisen--Zatespin--Kuzmin cutoff is confirmed by the future data.

The H.E.S.S. collaboration recently performed the accurate measurements of the gamma-ray spectra above 200 GeV for two BL Lac-type AGN in the Southern sky: PKS 2155-304 and PKS 2005-489. The TeV spectrum of the BL Lac object 1ES 2344+514 has also been recently measured by the Whipple collaboration. Using the results of phenomenological calculations of the SED of the EBL we have extracted the intrinsic gamma-ray spectra of these blazars. Consequently we have shown that these spectra can be fitted well by means of a synchrotron self Compton model.

We propose a method for estimating the composition, i.e. the relative amounts of leptons and protons, of extragalactic jets which exhibit X-ray bright knots in their kpc scale jets. The method relies on measuring, or setting upper limits on, the component of the Cosmic Microwave Background (CMB) radiation that is bulk-Comptonized by cold electrons in the relativistically flowing jet. These measurements, along with modeling of the broadband knot emission that constrain the bulk Lorentz factor of the jets, can yield estimates of the jet power carried by protons and leptons. We provide an explicit calculation of the spectrum of the bulk-Comptonized (BC) CMB component and apply these results to PKS 0637--752 and 3C 273, two superluminal quasars with Chandra-detected large scale jets. What makes these sources particularly suited for such a procedure is the absence of significant non-thermal jet emission in the `bridge', the region between the core and the first bright jet knot, which guarantees that most of the electrons are cold there, leaving the BC scattered CMB radiation as the only significant source of photons in this region. At lambda=3.6-8.0 microns the most likely band to observe the BC scattered CMB emission, the Spitzer angular resolution (~ 1''-3) is considerably smaller than the `bridges' of these jets (~10''), making it possible to both measure and resolve this emission.

Owing to their isotropy, it is generally believed that ultrahigh energy cosmic rays (UHECRs) are extragalactic in origin. It is then expected that interactions of these cosmic rays with photons of the cosmic background radiation (CBR) should produce a drastic reduction in their flux above and energy of about $5 \times 10^{19}$ eV (50 EeV), the so-called ``GZK effect''. At present, the existence of this effect is uncertain owing to conflicting observational data and small number statistics. We show here that a small amount of Lorentz invariance violation (LIV), which could turn off photomeson interactions of UHECRs with the CBR, could explain the UHECR spectrum as measured by {\it AGASA} which shows an excess of UHECRs at energies above 100 EeV. If new results from the {\it Auger} array agree with the {\it AGASA} spectrum, this may be interpreted as evidence for a small amount of LIV. If, on the other hand, the new results are consistent with the {\it HiRes} results favoring a GZK effect, this would place severe constraints on LIV and, by implication, on some Planck scale quantum gravity models. We also discuss the power requirements needed to explain the UHECR spectrum for a range of assumptions, including source evolution and LIV and show that in all cases our results disfavor a $\gamma$-ray burst origin for the UHECRs.

We propose a method for estimating the matter content of quasar jets which exhibit Chandra - detected knots in their kpc scale jets. The method relies on measuring the component of the Cosmic Microwave Background (CMB) radiation that is bulk-Comptonized (BC) by the cold electrons in the relativistically flowing jet. We apply our method to PKS 0637 - 752, a superluminal quasar with an one -- sided Chandra - detected large scale jet. What makes this source particularly suited for such a procedure is the absence of significant non-thermal jet emission in the 'bridge', the region between the core and the first bright knot, guaranteeing that most of the electrons in the bridge are cold, and leaving the BC scattered CMB radiation as the only significant source of photons in this region. At lambda=3.6-8.0 microns, the most likely band for the BC scattered emission to appear, the angular resolution of Spitzer (~ 1''- 3'') is considerably smaller than the bridge of PKS 0637 -- 752 (~8''), making it possible to both measure and resolve this emission.

Multi-TeV spectra of Mkn 501 and other BL Lac objects exhibit the high energy cutoffs predicted to be the result of intergalactic annihilation interactions, primarily with IR photons having the flux level determined by various astrophysical observations. After correcting for such intergalactic absorption, these spectra can be explained within the framework of synchrotron self-Compton emission models. Stecker and Glashow have shown that the existence of this annihilation via electron-positron pair production puts strong constriants on Lorentz invariance violation. Such constriants have important implications for some quantum gravity and large extra dimension models. A much smaller amount of Lorentz invariance violation has potential implications for understanding the spectra of ultrahigh energy cosmic rays.

The goal of the Orbiting Wide-field Light-collectors (OWL) mission is to study the origin and physics of the highest energy particles known in nature, the ultrahigh energy cosmic rays (UHECRs). The OWL mission consists of telescopes with UV sensitive cameras on two satellites operating in tandem to view in stereo the development of the giant particle showers induced in the Earth's atmosphere by UHECRs. This paper discusses the characteristics of the OWL mission.

Constraints on possible Lorentz symmetry violation (LV) of order E/M_Planck for electrons and photons in the framework of effective field theory (EFT) are discussed. Using (i) the report of polarized MeV emission from GRB021206 and (ii) the absence of vacuum \v{C}erenkov radiation from synchrotron electrons in the Crab nebula, we improve previous bounds by 10^-10 and 10^-2 respectively. We also show that the LV parameters for positrons and electrons are different, discuss electron helicity decay, and investigate investigate how prior constraints are modified by the relations between LV parameters implied by EFT.

We consider the advantages of and the problems associated with hypotheses to explain the origin of ulthrahigh energy cosmic rays (UHECR: E > 10 EeV) and the "trans GZK" cosmic rays (TGZK: E > 100 EeV), both through "old physics" (acceleration in cosmic sources) and "new physics" (new particles, topological defects, fat neutrino cross sections, Lorentz invariance violation).

Cosmic rays have been observed up to energies $10^8$ times larger than those of the best particle accelerators. Studies of astrophysical particles (hadrons, neutrinos and photons) at their highest observed energies have implications for fundamental physics as well as astrophysics. Thus, the cosmic high energy frontier is the nexus to new particle physics. This overview discusses recent advances being made in the physics and astrophysics of cosmic rays and cosmic gamma-rays at the highest observed energies as well as the related physics and astrophysics of very high energy cosmic neutrinos. These topics touch on questions of grand unification, violation of Lorentz invariance, as well as Planck scale physics and quantum gravity.

Observations of the multi-TeV spectra of the nearby BL objects Mkn 421 and Mkn 501 exhibit the high energy cutoffs predicted to be the result of intergalactic annihilation interactions, primarily with infrared photons having a flux level as determined by various astronomical observations. After correction for this absorption effect, the derived intrinsic spectra of these multi-TeV sources can be explained within the framework of simple synchrotron self-Compton emission models. Stecker and Glashow have shown that the existence of such annihilations via electron-positron pair production interactions up to an energy of 20 TeV puts strong constraints on Lorentz invariance violation. Such constraints have important implications for Lorentz invariance violating (LIV) quantum gravity models as well as LIV models involving large extra dimensions. We also discuss the implications of observations of high energy gamma-rays from the Crab Nebula on constraining quantum gravity models.

Observations of the multi-TeV spectra of the nearby BL Lac objects Mkn 421 and Mkn 501 exhibit the high energy cutoffs predicted to be the result of intergalactic annihilation interactions, primarily with IR photons having a flux level as determined by various astronomical observations. After correcting for such intergalactic absorption, these spectra can be explained within the framework of synchrotron self-Compton emission models. Stecker & Glashow have shown that the existence of this annihilation via electron-positron pair production puts strong constraints on Lorentz violation. Such constraints have important implications for quantum gravity models and large extra dimension models. We also discuss the implications of observations of high energy gamma-rays from the Crab Nebula on constraining quantum gravity models.

We discuss the results of modelling the TeV gamma-ray spectra of two AGNs, Mkn 501 and Mkn 421 that have almost the same redshifts: z=0.031 and z=0.034, respectively. The effect of intergalactic gamma-ray absorption is treated as an uncertainty in the measurement of the intrinsic spectrum. Although the objects differ, we obtain satisfactory fits for both of them in a synchrotron self-Compton scenario. Compared to previous models, our fits are characterised by higher values of the Doppler factor (>= 50) and an electron injection spectrum extending to higher energies (Gmax = 1.5x10^5). In the case of Mkn 421, the observed difference in spectral slope in X-rays and TeV gamma-rays between the high and low states can be explained as a variation of a single parameter - the maximum energy Gmax mc^2 at which electrons are injected.

I discuss recent advances being made in the physics and astrophysics of cosmic rays and cosmic gamma-rays at the highest observed energies as well as the related physics and astrophysics of very high energy cosmic neutrinos. I also discuss the connections between these topics.

The primary ultrahigh energy particles which produce giant extensive air showers in the Earth's atmosphere present an intriguing mystery from two points of view: (1) How are these particles produced with such astounding energies, eight orders of magnitude higher than those produced by the best man-made terrestrial accelerators? (2) Since they are most likely extragalactic in origin, how do they reach us from extragalactic distances without suffering the severe losses expected from interactions with the 2.7 K thermal cosmic background photons -- the so-called GZK effect? The answers to these questions may involve new physics: violations of special relativity, grand unification theories, and quantum gravity theories involving large extra dimensions. They may involve new astrophysical sources, "zevatrons". Or some heretofore totally unknown physics or astrophysics may hold the answer. I will discuss here the mysteries involving the production and extragalactic propagation of ultrahigh energy cosmic rays and some suggested possible solutions.

I will give here an overview of the present observational and theoretical situation regarding the question of the matter-antimatter asymmetry of the universe and the related question of the existence of antimatter on a cosmological scale. I will also give a simple discussion of the role of CP violation in this subject.

Using the recent models of Malkan & Stecker (2001) for the infrared background radiation and extrapolating them into the optical and UV range using recent galaxy count data, we rederive the optical depth of the Universe to high energy gamma-rays as a function of energy and redshift for energies between 50 GeV and 100 TeV and redshifts between 0.03 and 0.3. We then use these results to derive the intrinsic gamma-ray spectrum of Mkn 501 during its 1997 high state. We find that the time averaged spectral energy distribution of Mkn 501 while flaring had a broad, flat peak in the 5 to 10 TeV range which corresponds to a broad, flat time averaged X-ray peak in the 50 to 100 keV range observed during the flare. The spectral index of our derived intrinsic differential photon spectrum for Mkn 501 at energies below about 2 TeV was found to be apx. 1.6 to 1.7. This corresponds to a time averaged spectral index of 1.76 found in soft X-rays at energies below the X-ray (synchrotron) peak. These results appear to favor a synchrotron-self Compton origin for the TeV emission together wuth jet parameters which are consistent with time variability constraints within the context of a simple SSC model.

The COMPTEL and EGRET detectors aboard the Compton Gamma-Ray Observatory measured an extragalactic gamma-ray background extending from MeV energies up to about 100 GeV. Calculations performed making reasonable assumptions indicate that blazars can account for the background between about 10 MeV and at least 10 GeV. Below 30 MeV, the background flux and spectrum are not very well determined and a dedicated satellite detector will be required to remedy this situation. Below 10 MeV, supernovae and possibly AGN may contribute to the extragalactic background flux. Above 10 GeV, the role of blazars in contributing to the background is unclear because we do not have data on their spectra at these energies and because theoretical models predict that many of them will have spectra which should cut off in this energy range. At these higher energies, a new component, perhaps from topological defects, may contribute to the background, as well as X-ray selected BL Lac objects. The future GLAST detector should provide important data on the emission of extragalactic sources above 10 GeV and help resolve this issue. GLAST may also be able to detect the signature of intergalactic absorption by pair production interactions of background gamma-rays of energy above 20 GeV with starlight photons, this signature being a steepening of the background spectrum.

It has been suggested that cosmological gamma-ray bursts (GRBs) can produce the observed flux of cosmic rays at the highest energies. However, recent studies of GRBs indicate that their redshift distribution likely follows that of the average star formation rate and that GRBs were more numerous at high redshifts. As a consequence, we show that photomeson production energy losses suffered by ultrahigh energy cosmic rays coming from GRBs would produce too sharp a spectral high energy cutoff to be consistent with the air shower data. Furthermore, we show that cosmological GRBs fail to supply the energy input required to account for the cosmic ray flux above 10 EeV by a factor of 100-1000.

We use the recent reanalysis of multi-TeV gamma-ray observations of Mrk 501 to constrain the Lorentz invariance breaking parameter involving the maximum electron velocity. Our limit is two orders of magnitude better than that obtained from the maximum observed cosmic-ray electron energy.

We predict luminosity functions and number counts for extragalactic infrared sources at various wavelengths using the framework of our empirically based model. Comparisons of our galaxy count results with existing data indicate that either galaxy luminosity evolution is not much stronger than Q=3.1, where L is proportional to (1+z)^{Q}, or that this evolution does not continue beyond a redshift of 2. However, a derivation of the far infrared background from COBE-DIRBE data suggests a stronger evolution for far infrared emission with Q > 4 in the redshift range beteen 0 and 1. We discuss several interpretations of these results and also discuss how future observations can reconcile this apparent conflict. We also make predictions of the redshift distributions of extragalactic infrared sources at selected flux levels which can be tested by planned detectors. Finally, we predict the fluxes at which various future surveys will become confusion limited.

In this paper, I will take a synoptic approach to determining the intergalactic infrared radiation field (IIRF). This approach draws on both the multi-TeV gamma-ray observations and the infrared background observations and relates them via the semi-empirical modelling method of Malkan and Stecker. I discuss the evidence for an intergalactic infrared background obtained by an analysis of the HEGRA observations of the high energy gamma-ray spectrum of Mrk 501 and from constraints from Mrk 421 deduced from the Whipple air Cherenkov telescope results. I will show that this evidence is in accord with the predictions made by Malkan and Stecker for the intergalactic infrared spectral energy distribution produced by galaxies. The Malkan-Stecker predictions are also in excellent agreement with mid- and far infrared galaxy counts. However, there may be potential problems relating these predictions with the results of the analysis of COBE-DIRBE far infrared data. The gamma-ray and COBE-DIRBE observations may also need to be reconciled. I will discuss possible ways to resolve this situation including a partial nullification of the $\gamma$-ray absorption process which can hypothetically occur if Lorentz invariance is broken.

These lectures discuss the mysteries involving the production and extragalactic propagation of ultrahigh energy cosmic rays and suggested possible solutions.

It has been suggested that cosmological gamma-ray bursts (GRBs) can produce the observed flux and spectrum of cosmic rays at the highest energies. However, recent observations indicate that the redshift distribution of GRBs most likely follows that of the star formation rate in the universe, a rate which was much higher at redshifts 1.5-2 than it is today. Thus, most GRBs are at high redshifts. As a consequence, any cosmic rays emitted by these GRBs at energies above 2-3 X 10^{19} eV would be strongly attenuated by interactions with the 3K background radiation. If one assumes rough equality between the energy released in 10^{-2} to 1 MeV photons and that released in 10^{20} eV cosmic rays, then less than 10 per cent of the cosmic rays observed above 10^{20} eV can be accounted for by GRBs.

Astronomy at the highest energies observed must be performed by studying neutrinos rather than photons because the universe is opaque to photons of these energies. By making observations of neutrinos with energies above 10 EeV one can deduce information about the distribution and cosmological history of cosmic rays which may have been accelerated to energies above 100 EeV. Very large fluxes of neutrinos at these energies which exhibit a hard spectrum would be evidence of a GUT phase transition or superheavy dark matter particles. Neutrinos at energies around a PeV may be produced in observable quantities in active galaxies from the acceleration and interactions of cosmic rays in the vicinity of massive black holes or in relativistic jets produced by these black holes. Cosmological gamma-ray bursts may also produce PeV neutrinos. Neutrinos at these energies may oscillate into tau neutrinos which can travel through the Earth by the process of regeneration and produce upward moving air-showers. Neutrinos of energies above of 1000 EeV can produce ultrahigh energy cosmic rays thorugh interactions with big-bang neutrinos, producing Z bosons. This "Z-burst" phenomenon can be used to study the mass and distribution of big-bang neutrinos.

In this paper, we suggest a new hypothesis for explaining the spectrum of the extragalactic MeV gamma-ray background as observed by COMPTEL and SMM. We propose that both the flux level and spectrum can be accounted for as a superposition of non-thermal MeV tails in the spectra of Seyfert galaxies and other AGN. Although present detectors are not sensitive enough to obtain MeV data from individual extragalactic sources, indirect evidence in support of our hypothesis is found in OSSE and COMPTEL observations of the galactic black hole candidate Cygnus X-1.

We show that GLAST should be able to support or rule out the unresolved-blazar hypothesis for the origin of the diffuse, extragalactic gamma-ray background.

The pair-production absorption of high-energy gamma-rays by intergalactic low-energy photons is expected to produce a high-energy cutoff in the spectra of extragalactic sources which is a sensitive function of redshift. We first discuss the expected absorption coefficient as a function of energy and redshift derived by Stecker and De Jager by making use of a new empirically based calculation of the spectral energy distribution of the intergalactic infrared radiation field as given by Malkan and Stecker. We then discuss the fact that new data on the high energy gamma-ray source Mrk 501 appear to show the amount of intergalactic absorption predicted. The implications of this new HEGRA data, should they be confirmed, are significant for the astrophysics of this source, implying that (1) there is no significant intrinsic absorption inside the source, and (2) the physics of the emission produces a power-law spectrum to energies above 20 TeV. As a further test for intergalactic absorption, we give a predicted spectrum, with absorption included, for PKS 2155-304. This XBL lies at a redshift of 0.12, the highest redshift source yet observed at an energy above 0.3 TeV. We also discuss the determination of the gamma-ray opacity of the universe at higher redshifts (out to $z=3$), following the treatment of Salamon and Stecker.

The recent HEGRA observations of the blazar Mkn 501 show strong curvature in the very high energy gamma-ray spectrum. Applying the gamma-ray opacity derived from an empirically based model of the intergalactic infrared background radiation field (IIRF), to these observations, we find that the intrinsic spectrum of this source is consistent with a power-law: dN/dE~ E^-alpha with alpha=2.00 +/- 0.03 over the range 500 GeV - 20 TeV. Within current synchrotron self-Compton scenarios, the fact that the TeV spectral energy distribution of Mkn 501 does not vary with luminosity, combined with the correlated, spectrally variable emission in X-rays, as observed by the BeppoSAX and RXTE instruments, also independently implies that the intrinsic spectrum must be close to alpha=2. Thus, the observed curvature in the spectrum is most easily understood as resulting from intergalactic absorption.

We discuss the determination of the intergalactic pair-production absorption coefficient (Stecker & De Jager 1998) making use of an empirically based calculation of the IR background spectrum (Malkan & Stecker 1998) which agrees well with recent data and constraints on the IR background. While the Whipple observations of Mrk421 hint at extragalactic absorption, the new HEGRA observations of Mrk501 appear to strongly indicate such absorption. We discuss the determination of absorption at higher redshifts (Salamon & Stecker 1998). We also give a predicted spectrum, with absorption included, for PKS2155-304, which, at a redshift of 0.12, is the highest redshift TeV source yet observed.

We present the results of a new calculation of the photodisintegration of ultrahigh energy cosmic ray (UHCR) nuclei in intergalactic space. The critical interactions for energy loss and photodisintegration of UHCR nuclei occur with photons of the 2.73K cosmic background radiation (CBR) and with photons of the infrared background radiation (IBR). We have reexamined this problem making use of a new determination of the IBR based on empirical data, primarily from IRAS galaxies, consistent with direct measurements and upper limits from TeV gamma- ray observations. We have also improved the calculation by including the specific threshold energies for the various photodisintegration interactions in our Monte Carlo calculation. With the new smaller IBR flux, the steepness of the Wien side of the now relatively more important CBR makes their inclusion essential for more accurate results. Our results indicate a significant increase in the propagation time of UHCR nuclei of a given energy over previous results. We discuss the possible significance of this for UHCR origin theory.

We present a new calculation of the intergalactic gamma-ray absorption coefficient as a function of both energy and redshift. In reexamining this problem, we make use of a new, empirically based calculation (as opposed to previous model calculations) of the intergalactic infrared radiation field. We find smaller opacities than those given previously (Stecker & De Jager 1997). We apply our results to the new observations of the flaring gamma-ray spectra of Mrk421 and Mrk501, both at a redshift of apx. 0.03. Our new calculations indicate that there should be no significant curvature in the spectra of these sources for energies below 10 TeV, as indicated by recent observations. However, the intrinsic spectra of these sources should be harder by apx. 0.2 to 0.45 in the spectral index in the 1 to 10 TeV range with an intergalactic absorption cutoff above apx. 20 TeV.

Cosmic topological defects in a wide class of supersymmetric theories can simultaneously be sources of higgs particles of mass comparable to the supersymmetry breaking scale $\sim$ TeV, as well as superheavy gauge bosons of mass $\sim\eta$, where $\eta$ (>> 1 TeV) denotes the associated gauge symmetry breaking scale. For cosmic strings with $\eta\sim 10^{14}$ GeV, the higgs decay can account for the extragalactic diffuse gamma ray background above $\sim$ 10 GeV, while the gauge boson decay may explain the highest energy cosmic ray flux above $\sim 10^{11}$ GeV, provided that particle production is the dominant energy loss mechanism for cosmic strings, as recent simulations suggest. By the same token, cosmic strings with $\eta$ much above $\sim 10^{14}$ GeV will be ruled out.

We present the results of a new estimation of the photodisintegration and propagation of ultrahigh energy cosmic ray (UHCR) nuclei in intergalactic space. The critical interactions for photodisintegration and energy loss of UHCR nuclei occur with photons of the infrared background radiation (IBR). We have reexamined this problem making use of a new determination of the IBR based on empirical data, primarily from IRAS galaxies, and also collateral information from TeV gamma-ray observations of two nearby BL Lac objects. Our results indicate that a 200 EeV Fe nucleus can propagate apx. 100 Mpc through the IBR. We argue that it is possible that the highest energy cosmic rays observed may be heavy nuclei.

We discuss the problem of the absorption of very high-energy gamma-rays by pair production interactions with extragalactic photons which originate from stellar emission in the near IR-UV and reradiation of starlight in the mid- and far-IR. The absorption of gamma-rays above 1 TeV is dominated by interactions with infrared photons. We make a new determination of the optical depth of the universe to multi-TeV photons as a function of energy and redshift and use the results to compare with recent spectral data of Mrk 421 and Mrk 501, sources that have been observed in the flaring state up to apx. 10 TeV energy. For the optical depth calculations, we have made use of a new, empirically based calculation of the intergalactic radiation field by Malkan & Stecker which we consider to be more accurate than that based on previous theoretical modeling. We also discuss the absorption of sub-TeV gamma-rays by starlight photons at high redshifts.

Using the excellent observed correlations among various infrared wavebands with 12 and 60 micron luminosities, we calculate the 2-300 micron spectra of galaxies as a function of luminosity. We then use 12 micron and 60 micron galaxy luminosity functions derived from IRAS data, together with recent data on the redshift evolution of galaxy emissivity, to derive a new, empirically based IR background spectrum from stellar and dust emission in galaxies. Our best estimate for the IR background is of order 2-3 nW/m^2/sr with a peak around 200 microns reaching 6-8 nW/m^2/sr. Our empirically derived background spectrum is fairly flat in the mid-IR, as opposed to spectra based on modeling with discrete temperatures which exhibit a "valley" in the mid-IR. We also derive a conservative lower limit to the IR background which is more than a factor of 2 lower than our derived flux.

The all-sky surveys in gamma-rays by the imaging Compton telescope (COMPTEL) and the Energetic Gamma Ray Experiment Telescope (EGRET) on board the Compton Gamma Ray Observatory for the first time allows detailed studies of the extragalactic diffuse emission at gamma-ray energies greater 1 MeV. A preliminary analysis of COMPTEL data indicates a significant decrease in the level of the derived cosmic diffuse emission from previous estimates in the 1-30 MeV range, with no evidence for an MeV-excess, at least not at the levels claimed previously. The 1-30 MeV flux measurements are compatible with power-law extrapolation from lower and higher energies. These new results indicate that the possible contributions to the extragalactic emission from processes that could explain the MeV-excess, such as matter-antimatter annihilation, is significantly reduced. At high energies ($>$ 30 MeV), the extragalactic emission is well described by a power law photon spectrum with an index of -(2.10+-0.03) in the 30 MeV to 100 GeV energy range. No large scale spatial anisotropy or changes in the energy spectrum are observed in the deduced extragalactic emission. The most likely explanation for the origin of this extragalactic gamma-ray emission above 10 MeV, is that it arises primarily from unresolved gamma-ray-emitting blazars. The extension of the power law spectrum to 100 GeV implies the average emission from gamma-ray blazars may extend to 100 GeV.

We extend earlier calculations of the attenuation suffered by gamma rays during their propagation from extragalactic sources, obtaining new extinction curves for gamma rays down to 10 GeV in energy, from sources up to a redshift of z=3.

We present results of a calculation of absorption of 10-500 GeV gamma-rays at high redshifts. This calculation requires the determination of the high- redshift evolution of the full spectral energy distribution of the intergalactic photon field. For this, we have primarily followed the recent analysis of Fall, Charlot and Pei. We give our results for the gamma-ray opacity as a function of redshift out to a redshift of 3. We then give predicted gamma-ray spectra for selected blazars and also extend our results on the background from unresolved blazars to an energy of 500 GeV. Absorption effects are predicted to significantly steepen the background spectrum above 20 GeV. Our absorption calculations can be used to place limits on the redshifts of gamma-ray bursts. Our background calculations can be used to determine the observability of multi-GeV lines from dark matter neutralino particles.

We calculate the absorption of 10-500 GeV gamma-rays at high redshifts. This calculation requires the determination of the high-redshift evolution of the intergalactic starlight photon field, including its IR-UV spectral energy distribution. To estimate this evolution, we have followed a recent analysis of Fall, Charlot and Pei which gives results consistent with recent data. We give our results for the gamma-ray opacity as a function of redshift out to a redshift of 3. We also give predicted gamma-ray spectra for selected blazars and give an extragalactic unresolved blazar background spectrum up to 500 GeV. Our results indicate that this background should steepen significantly above 20 GeV owing to intergalactic absorption. Future observations of this steepening would provide a test for the blazar background origin hypothesis. We have used our results to discuss upper limits on the redshifts of gamma-ray bursts. We note that the 17 Feb. 1994 burst observed by EGRET must have originated at a redshift less than 2. We also use our estimates of the background to determine the observability of multi-GeV gamma-ray lines from the annihilation of supersymmetric dark matter particles in the galactic halo.

We suggest that low-redshift XBLs (X-ray selected BL Lacertae objects) may be the only extragalactic gamma-ray sources observable at TeV energies. We use simple physical considerations involving synchrotron and Compton component spectra for blazars to suggest why the observed TeV sources are XBLs, whereas mostly RBLs and FSRQs are seen at GeV energies. These considerations indicate that the differences between XBLs and RBLs cannot be explained purely as relativistic jet orientation effects. We note that the only extragalactic TeV sources which have been observed are XBLs and that a nearby RBL with a very hard spectrum in the GeV range has not been seen at TeV energies. We also note that of the 14 BL Lacs observed by EGRET, 12 are RBLs, whereas only 2 are XBLs. We give a list of nearby XBLs which we consider to be good candidate TeV sources and predict estimated TeV fluxes for these objects.

We present a new calculation of the intergalactic $\gamma$-ray pair-production absorption coefficient as a function of both energy and redshift up to the redshift of 3C279, z = 0.54. In reexamining this problem, we make use of new observational data on the intergalactic infrared radiation field (IIRF), together with recent theoretical models of the galactic spectral energy distributions of the IIRF from stars and dust reradiation and estimates of the IIRF from galaxy counts and {\it COBE} results. We present our results for two fairly well defined IIRF spectral energy distributions, one of which is within $1 \sigma $ of our previous estimate of the IIRF at $ \sim 20$ $\mu$m. We then apply our results to the $\gamma$-ray spectrum of Mrk 421, and obtain good agreement with the observational data, including the recent results of the {\it HEGRA} group.

High energy gamma-rays and neutrinos can be produced both by the annihilation and by the possible slow decay of dark matter particles. We discuss the fluxes and spectra of such secondaries produced by dark matter particles in the universe and their observability in competition with other astrophysical gamma-ray signals and with atmospheric neutrinos. To do this, we work within the assumption that the dark matter particles are neutralinos which are the lightest supersymmetric particles (LSPs) predicted by supersymmetry theory.

We present the results of a new model calculation of the gamma-ray background produced by unresolved blazars, using the second EGRET catalogue and taking account of flaring. These results are compared to the preliminary gamma-ray background spectrum reported recently by the EGRET team. We find that blazars can account for the entire extragalactic gamma-ray background observed by EGRET. In addition, the EGRET spectrum shows some evidence of a signature for the blazar background, in the form of a concavity in the spectrum.

We review and clarify the assumptions of our basic model for neutrino production in the cores of quasars, as well as those modifications to the model subsequently made by other workers. We also present a revised estimate of the neutrino background flux and spectrum obtained using more recent empirical studies of quasars and their evolution. We compare our results with other thoeretical calculations and experimental upper limits on the AGN neutrino background flux. We also estimate possible neutrino fluxes from the jets of blazars detected recently by the EGRET experiment on the Compton Gamma Ray Observatory. We discuss the theoretical implications of these estimates.

Following our previously proposed technique, we have used the recent gamma-ray observations of Mrk 421 to place theoretically significant constraints on and possible estimates of the intergalactic infrared radiation field (IIRF) which are consistent with normal galactic IR production by stars and dust and rule out exotic mechanisms proposed to produce a larger IIRF. Using models for the low energy intergalactic photon spectrum from microwave to UV energies, we calculate the opacity of inter- galactic space to gamma-rays as a function of energy and redshift. These calculations indicate that the GeV gamma-ray burst recently observed by the EGRET experiment on CGRO originates at a redshift less than approximately 1.5.