Regularized Einstein-Gauss-Bonnet (EGB) theory of gravity in four dimensions is a new attempt to include nontrivial contributions of Gauss-Bonnet term. In this paper, we make a detailed analysis on possible constraints of the model parameters of the theory from recent cosmological observations, and some theoretical constraints as well. Our results show that the theory with vanishing bare cosmological constant, $\Lambda_0$, is ruled out by the current observational value of $w_{de}$, and the observations of GW170817 and GRB 170817A as well. For nonvanishing bare cosmological constant, instead, our results show that the current observation of the speed of GWs measured by GW170817 and GRB 170817A would place a constraints on $\tilde{\alpha}$, a dimensionless parameter of the theory, as $-7.78\times10^{-16}\le \tilde{\alpha} \leq 3.33\times10^{-15}$.

In the third catalog of active galactic nuclei detected by the $Fermi$ Large Area Telescope Clean (3LAC) sample, there are 402 blazars candidates of uncertain type (BCU). The proposed analysis will help to evaluate the potential optical classification flat spectrum radio quasars (FSRQs) versus BL Lacertae (BL Lacs) objects of BCUs, which can help to understand which is the most elusive class of blazar hidden in the Fermi sample. By studying the 3LAC sample, we found some critical values of $\gamma$-ray photon spectral index ($\Gamma_{\rm ph}$), variability index (VI) and radio flux (${\rm F_R}$) of the sources separate known FSRQs and BL Lac objects. We further utilize those values to defined an empirical "high-confidence" candidate zone that can be used to classify the BCUs. Within such a zone ($\Gamma_{\rm ph}<2.187$, log${\rm F_R}<2.258$ and ${ \rm logVI <1.702}$), we found that 120 BCUs can be classified BL Lac candidates with a higher degree of confidence (with a misjudged rate $<1\%$). Our results suggest that an empirical "high confidence" diagnosis is possible to distinguish the BL Lacs from the Fermi observations based on only on the direct observational data of $\Gamma_{\rm ph}$, VI and ${\rm F_R}$.

We propose a new class of gravity theories which are characterized by a nontrivial coupling between the gravitational metric and matter mediated by an auxiliary rank-2 tensor. The actions generating the field equations are constructed so that these theories are equivalent to general relativity in a vacuum, and only differ from general relativity theory within a matter distribution. We analyze in detail one of the simplest realizations of these generalized coupling theories. We show that in this case the propagation speed of gravitational radiation in matter is different from its value in vacuum and that this can be used to weakly constrain the (single) additional parameter of the theory. An analysis of the evolution of homogeneous and isotropic spacetimes in the same framework shows that there exist cosmic histories with both an inflationary phase and a dark era characterized by a different expansion rate.

The thickness of a slim disk is determined by the balance between the radiation force and the vertical component of the gravity of the black hole (BH). It was found that vertical gravity increases with the disk height, and it will decrease with the disk height if the disk thickness is above a critical value, which implies that gas at the disk surface may be driven into outflows by radiation force when the disk thickness surpasses the critical value. In this work, we derive a global solution to a slim disk with radiation-driven outflows. We find that the outflows are driven from the disk surface if the mass accretion rate $\dot{m} \gtrsim 1.78-1.91$ ($\dot{m}=\dot{M}/\dot{M}_{\rm Edd}$, and $\dot{M}_{\rm Edd}=L_{\rm Edd}/0.1c^2$) depending on BH mass, while the outflows are suppressed when the mass accretion rate is lower than this critical value. The mass accretion rate decreases with decreasing radius in the disk with outflows, and the rate of the gas swallowed by the BH is always limited to $\dot{m}_{\rm in}\sim 1.78-1.91$ even if the mass accretion rate at the outer edge of the disk is very high. This may set constraints on the massive BH growth through accretion in the early Universe. Due to the presence of outflows, there is an upper limit on the radiation flux of the disk, which leads to saturation of the continuum spectra of the disk with outflows at the high energy band. This may be tested by the observations of quasars or/and BH X-ray binaries.

Precision measurements of charged cosmic rays have recently been carried out by space-born (e.g. AMS-02), or ground experiments (e.g. HESS). These measured data are important for the studies of astro-physical phenomena, including supernova remnants, cosmic ray propagation, solar physics and dark matter. Those scenarios usually contain a number of free parameters that need to be adjusted by observed data. Some techniques, such as Markov Chain Monte Carlo and MultiNest, are developed in order to solve the above problem. However, it is usually required a computing farm to apply those tools. In this paper, a genetic algorithm for finding the optimum parameters for cosmic ray injection and propagation is presented. We find that this algorithm gives us the same best fit results as the Markov Chain Monte Carlo but consuming less computing power by nearly 2 orders of magnitudes.

In the third catalog of active galactic nuclei detected by the Fermi-LAT (3LAC) Clean Sample, there are 402 blazars candidates of uncertain type (BCU). Due to the limitations of astronomical observation or intrinsic properties, it is difficult to classify blazars using optical spectroscopy. The potential classification of BCUs using machine learning algorithms is essential. Based on the 3LAC Clean Sample, we collect 1420 Fermi blazars with 8 parameters of {\gamma}-ray photon spectral index, radio flux, flux density, curve significance, the integral photon flux in 100 to 300 MeV, 0.3 to 1 GeV, 10 to 100 GeV and variability index. Here, we apply 4 different supervised machine learning (SML) algorithms (\emph{Decision trees, Random forests, support vector machines and Mclust Gaussian finite mixture models}) to evaluate the classification of BCUs based on the direct observational properties. All the 4 methods can perform exceedingly well with a more accuracy and can effective forecast the classification of Fermi BCUs. The evaluating results show the results of these methods (SML) are valid and robust, where, about 1/4 sources are FSRQs and 3/4 are BL Lacs in 400 BCUs, which are consistent with some other recent results. Although a number of factors influence the accuracy of SML, the results are stable at a fixed ratio 1:3 between FSRQs and BL Lacs, which suggests that the SML can provides an effective method to evaluate the potential classification of BCUs. Among the 4 methods, Mclust Gaussian Mixture Modelling has the highest accuracy for our training sample (4/5, seed=123).

The AMS-02 experiment has reported a new measurement of the antiproton/proton ratio in Galactic cosmic rays (CRs). In the energy range $E\sim\,$60-450 GeV, this ratio is found to be remarkably constant. Using recent data on CR proton, helium, and carbon fluxes, 10Be/9Be and B/C ratios, we have performed a global Bayesian analysis based on a Markov-Chain Monte-Carlo sampling algorithm under a "two halo model" of CR propagation. In this model, CRs are allowed to experience a different type of diffusion when they propagate in the region close of the Galactic disk. We found that the vertical extent of this region is about 900 pc above and below the disk, and the corresponding diffusion coefficient scales with energy as $D\sim\,E^{0.15}$, describing well the observations on primary CR spectra, secondary/primary ratios and anisotropy. Under this model we have carried out improved calculations of antiparticle spectra arising from secondary CR production and their corresponding uncertainties. We made use of Monte-Carlo generators and accelerator data to assess the antiproton production cross-sections and their uncertainties. While the positron excess requires the contribution of additional unknown sources, we found that the new AMS-02 antiproton data are consistent, within the estimated uncertainties, with our calculations based on secondary production.

The antiparticle energy spectra of Galactic cosmic rays (CRs) have several exciting features such as the unexpected positron excess at $E\sim$10-200\,GeV and the remarkably hard antiproton flux at $E\sim$\,60--450\,GeV recently measured by the \AMS{} experiment. In this paper, we report calculations of antiparticle CR spectra arising from secondary production and their corresponding uncertainties. Using the most recent data on CR protons, helium, carbon, and nuclear ratios $^{10}$Be/$^{9}$Be and B/C, we have performed a global Bayesian analysis, based on a Markov Chain Monte Carlo sampling algorithm, under a scenario of spatial-dependent CR diffusion in the Galaxy which reproduces well the observed spectral hardening in the CR hadron fluxes. While the high-energy positron excess requires the contribution of additional unknown sources, we found that the antiproton data are consistent within the estimated uncertainties, with our predictions based on secondary production.

A new analysis of cosmic ray (CR) data collected by the SOKOL experiment in space found that the deuteron-to-helium ratio at energies between 500 and 2000 GeV/nucleon takes the value d/He ~ 1.5. As we will show, this result cannot be explained by standard models of secondary CR production in the interstellar medium and points to the existence of a high-energy source of CR deuterons. To account for the deuteron excess in CRs, we argue that the only viable solution is hadronic interaction processes of accelerated particles inside old supernova remnants (SNRs). From this mechanism, however, the B/C ratio is also expected to increase at energy above ~50 of GeV/nucleon, in conflict with new precision data just released by the AMS-02 experiment. Hence, if this phenomenon is a real physical effect, hadronic production of CR deuterons must occur in SNRs characterized by low metal abundance. In such a scenario, the sources accelerating C-N-O nuclei are not the same as those accelerating helium or protons, so that the connection between d/He ratio and B/C ratio is broken, and the latter cannot be used to place constraints on the production of light isotopes or antiparticles.

Antimatter nuclei in cosmic rays (CRs) are a promising tool for the indirect detection of dark-matter annihilation signatures. However, the search of new-physics signals in CRs relies on our knowledge of the astrophysical antimatter background which, in turns, depends critically on the several fragmentation cross-sections that regulate production and destruction of antiparticles in the interstellar medium. In this work, we have re-evaluated the astrophysical background of CR antiproton, antineutron, and antihelium nuclei in Galactic CRs using improved calculations. The production cross-sections of individual antinucleons are constrained using updated calculations that make use of recent accelerator data. The production of antideuteron and antihelium nuclei is calculated using an improved model of nuclear coalescence that accounts for the asymmetry in antineutron and antiproton production. We discuss the cross-section induced uncertainties and show that they are dominating in comparison with other uncertainties of astrophysical origin.

The Alpha Magnetic Spectrometer (AMS) experiment has released new measurements of primary and secondary nuclei in cosmic rays at the TeV energy scale. Using these data, we present new results based on a global Bayesian analysis of our two halo model of CR transport.

Dark matter search in space has been carried out for many years. Measurements of cosmic ray photons, charged antiparticles and neutrinos are useful tools for dark matter indirect search. The antiparticle energy spectra of cosmic rays have several exciting features such as the unexpected positron excess at energy about 10 -- 500 GeV and the remarkably flattening antiproton/proton at energy about 60--450 GeV precisely measured by the AMS-02 experiment, which can not be explained simultaneously by secondary production in interstellar medium. In this work, we report a combined analysis of cosmic ray antiproton and positron spectra arising from dark matter on the top of a secondary production in a spatial-dependent propagation model. We discuss the systematics from antiproton production cross section using the two latest Monte Carlo generators, i.e. EPOS LHC and QGSJET-II-04m, respectively. We compare their results. In the case of EPOS LHC, we find that the dark matter pair annihilating into tau leptons channel with 100% branching ratio is the only possible one channel scenario to explain data. On the other hand, there is not a single possible channel in the case of QGSJET-II-04m. We also propose possible two-channel scenarios based on these two Monte Carlo generators.

The curvature of the $\gamma$-ray spectrum in blazars may reflect the intrinsic distribution of the emitting electron distribution, which will further give some information on the possible acceleration and cooling processes in the emitting region. The $\gamma$-ray spectra of Fermi blazars are normally fitted either by a single power-law (PL) or a log-normal (call Logarithmic Parabola, LP) form. The possible reason for this differnece is not unclear. We statistically explore this issue based on the different observational properties of 1419 Fermi blazars in the 3LAC Clean sample. We find that the $\gamma$-ray flux (100 MeV-100 GeV) and variability index follow bimodal distributions for PL and LP blazars, where $\gamma$-ray flux and variability index show {a positive correlation}. However, the distributions of the $\gamma$-ray luminosity and redshift follow a unimodal distribution. Our results suggest that the bimodal distribution of $\gamma$-ray flux for LP and PL blazars may be not intrinsic and all blazars may have an intrinsic curved $\gamma$-ray spectrum and the PL spectrum is just caused by the fitting effect due to the less photons.

NGC 4993, as a host galaxy of the electromagnetic counterpart of the first gravitational-wave detection of a binary neutron-star merger, was observed by many powerful telescopes from radio to $\gamma$-ray waveband. The weak nuclear activities of NGC 4993 suggest that it is a low-luminosity active galactic nuclear (LLAGN). We build the multi-waveband spectral energy distributions (SEDs) of NGC 4993 from literatures. We find that the radio spectrum at $\sim 100-300$ GHz is much steeper than that of low-frequency waveband (e.g., 6-100 GHz), where this break was also found in the supermassive black holes in our galaxy center (Sgr A*), and in some other nearby AGNs. The radio emission above and below this break may has different physical origins, which provide an opportunity to probe the accretion and jet properties. We model the multi-waveband SEDs of NGC 4993 with an advection-dominated accretion flow (ADAF)-jet model. We find that the high-frequency steep radio emission at millimeter waveband is consistent with the prediction of the ADAF, while the low-frequency flat radio spectrum is better fitted by the jet. Further more, the X-ray emission can be also explained by the ADAF model simultaneously. From the model fits, we estimate important parameters of the central engine (e.g., accretion rate near the horizon of the black hole and mass-loss rate in jet) for NGC 4993. This result strengthen that the millimiter, submillimer and deep X-ray observations are crucial to understand the weak or quiescent activities in supermassive black-hole systems. Further simultaneous millimeter and X-ray monitoring of this kind of LLAGNs will help us to better understand the physical origin of multiwaveband emission.

Dark matter may be charged under dark electromagnetism with a dark photon that kinetically mixes with the Standard Model photon. In this framework, dark matter will collect at the center of the Earth and annihilate into dark photons, which may reach the surface of the Earth and decay into observable particles. We determine the resulting signal rates, including Sommerfeld enhancements, which play an important role in bringing the Earth's dark matter population to their maximal, equilibrium value. For dark matter masses $m_X \sim$ 100 GeV - 10 TeV, dark photon masses $m_{A'} \sim$ MeV - GeV, and kinetic mixing parameters $\varepsilon \sim 10^{-10} - 10^{-8}$, the resulting electrons, muons, photons, and hadrons that point back to the center of the Earth are a smoking-gun signal of dark matter that may be detected by a variety of experiments, including neutrino telescopes, such as IceCube, and space-based cosmic ray detectors, such as Fermi-LAT and AMS. We determine the signal rates and characteristics, and show that large and striking signals---such as parallel muon tracks---are possible in regions of the $(m_{A'}, \varepsilon)$ plane that are not probed by direct detection, accelerator experiments, or astrophysical observations.

Dark matter may interact with the Standard Model through the kinetic mixing of dark photons, $A'$, with Standard Model photons. Such dark matter will accumulate in the Sun and annihilate into dark photons. The dark photons may then leave the Sun and decay into pairs of charged Standard Model particles that can be detected by the Alpha Magnetic Spectrometer. The directionality of this "dark sunshine" is distinct from all astrophysical backgrounds, providing an opportunity for unambiguous dark matter discovery by AMS. We perform a complete analysis of this scenario including Sommerfeld enhancements of dark matter annihilation and the effect of the Sun's magnetic field on the signal, and we define a set of cuts to optimize the signal probability. With the three years of data already collected, AMS may discover dark matter with mass 1 TeV $\lesssim m_X \lesssim$ 10 TeV, dark photon masses $m_{A'} \sim \mathcal O(100)$ MeV, and kinetic mixing parameters $10^{-10} \lesssim \varepsilon \lesssim 10^{-8}$. The proposed search extends beyond existing beam dump and supernova bounds, and it is complementary to direct detection, probing the same region of parameter space for elastic dark matter, but potentially far more in the case of inelastic dark matter.

Dark photons in the MeV to GeV mass range are important targets for experimental searches. We consider the case where dark photons $A'$ decay invisibly to hidden dark matter $X$ through $A' \to XX$. For generic masses, proposed accelerator searches are projected to probe the thermal target region of parameter space, where the $X$ particles annihilate through $XX \to A' \to \text{SM}$ in the early universe and freeze out with the correct relic density. However, if $m_{A'} \approx 2m_X$, dark matter annihilation is resonantly enhanced, shifting the thermal target region to weaker couplings. For $\sim 10\%$ degeneracies, we find that the annihilation cross section is generically enhanced by four (two) orders of magnitude for scalar (pseudo-Dirac) dark matter. For such moderate degeneracies, the thermal target region drops to weak couplings beyond the reach of all proposed accelerator experiments in the scalar case and becomes extremely challenging in the pseudo-Dirac case. Proposed direct detection experiments can probe moderate degeneracies in the scalar case. For greater degeneracies, the effect of the resonance can be even more significant, and both scalar and pseudo-Dirac cases are beyond the reach of all proposed accelerator and direct detection experiments. For scalar dark matter, we find an absolute minimum that sets the ultimate experimental sensitivity required to probe the entire thermal target parameter space, but for pseudo-Dirac fermions, we find no such thermal target floor.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

The millimeter bump, as found in high-resolution multi-waveband observations of M87, most possibly comes from the synchrotron emission of thermal electrons in advection dominated accretion flow(ADAF). It is possible to constrain the accretion rate near the horizon if both the nuclear millimeter emission and its polarization are produced by the hot plasma in the accretion flow. The jet power of M87 has been extensively explored, which is around $8_{\rm -3}^{+7}\times10^{42} {\rm erg/s}$ based on the analysis of the X-ray cavity. The black hole(BH) spin can be estimated if the jet power and the accretion rate near the horizon are known. We model the multi-wavelength spectral energy distribution (SED) of M87 with a coupled ADAF-jet model surrounding a Kerr BH, where the full set of relativistic hydrodynamical equations of the ADAF are solved. The hybrid jet formation model, as a variant of Blandford-Znajek model, is used to model the jet power. We find that the SMBH should be fast rotating with a dimensionless spin parameter $a_{*}\simeq0.98_{\rm -0.02}^{+0.012}$.

The cosmic ray energy spectra of protons and helium nuclei, which are the most abundant components of cosmic radiation, exhibit a remarkable hardening at energies above one hundred GeV/nucleon. Recent data from AMS-02 confirms this feature with a higher significance. This data challenges the current models of cosmic ray acceleration in Galactic sources and propagation in the Galaxy. Here, we explain the observed break in the spectra of proton and helium nuclei in light of recent advances of cosmic ray diffusion theories in turbulent astrophysical sources as a result of a transition between different cosmic ray diffusion regimes. We reconstruct the observed cosmic ray spectra using the fact that transition from normal diffusion to superdiffusion changes the efficiency of particle acceleration and causes the change in the spectral index. We find that calculated proton and helium spectra match very well with the data.

Recently the AMS-02 collaboration has published the measurement of the cosmic antiproton to proton ratio $\bar{p}/p$ and the $\bar{p}$ flux with a high precision up to $\sim 450\,\mathrm{GeV}$. In this work, we perform a systematic analysis of the secondary antiproton flux generated by the cosmic ray interaction with the interstellar gas. The uncertainty of the prediction originates from the cosmic ray propagation process and the hadronic interaction models. Although the cosmic ray propagation parameters have been well controlled by the AMS-02 $B/C$ ratio data for a specified model, different propagation models can not be discriminated by the $B/C$ data. The $\bar{p}$ flux is also calculated for several hadronic interaction models, which are generally adopted by the cosmic ray community. However, the results for different hadronic models do not converge. We find the EPOS LHC model, which seems to fit the collider data very well, predicts a slightly lower $\bar{p}/p$ ratio than the AMS-02 data at the high energy end. Finally we derive the constraints on the dark matter annihilation cross section from the AMS-02 $\bar{p}/p$ ratio for different propagation and hadronic interaction models.

MSSM4G models, in which the minimal supersymmetric standard model is extended to include vector-like copies of standard model particles, are promising possibilities for weak-scale supersymmetry. In particular, two models, called QUE and QDEE, realize the major virtues of supersymmetry (naturalness consistent with the 125 GeV Higgs boson, gauge coupling unification, and thermal relic neutralino dark matter) without the need for fine-tuned relations between particle masses. We determine the implications of these models for dark matter and collider searches. The QUE and QDEE models revive the possibility of heavy Bino dark matter with mass in the range 300-700 GeV, which is not usually considered. Dark matter direct detection cross sections are typically below current limits, but are naturally expected above the neutrino floor and may be seen at next-generation experiments. Indirect detection prospects are bright at the Cherenkov Telescope Array, provided the 4th-generation leptons have mass above 350 GeV or decay to taus. In a completely complementary way, discovery prospects at the LHC are dim if the 4th-generation leptons are heavy or decay to taus, but are bright for 4th-generation leptons with masses below 350 GeV that decay either to electrons or to muons. We conclude that the combined set of direct detection, CTA, and LHC experiments will discover or exclude these MSSM4G models in the coming few years, assuming the Milky Way has an Einasto dark matter profile.

M87 is arguably the best supermassive black hole (BH) to explore the jet and/or accretion physics due to its proximity and fruitful high-resolution multi-waveband observations. We model the multi-wavelength spectral energy distribution (SED) of M87 core that observed at a scale of 0.4 arcsec ($\sim 10^5R_{\rm g}$, $R_{\rm g}$ is gravitational radius) as recently presented by Prieto et al. Similar to Sgr A*, we find that the millimeter bump as observed by Atacama Large Millimeter/submillimeter Array (ALMA) can be modeled by the synchrotron emission of the thermal electrons in advection dominated accretion flow (ADAF), while the low-frequency radio emission and X-ray emission may dominantly come from the jet. The millimeter radiation from ADAF dominantly come from the region within $10R_{\rm g}$, which is roughly consistent with the recent very long baseline interferometry observations at 230\,GHz. We further calculate the Faraday rotation measure (RM) from both ADAF and jet models, and find that the RM predicted from the ADAF is roughly consistent with the measured value while the RM predicted from the jet is much higher if jet velocity close to the BH is low or moderate (e.g., $v_{\rm jet}\lesssim0.6\,c$). With the constraints from the SED modeling and RM, we find that the accretion rate close to the BH horizon is $\sim (0.2-1)\times10^{-3}{M}_{\odot} \rm yr^{-1}\ll\dot{\it M}_{\rm B}\sim 0.2\it {M}_{\odot} \rm yr^{-1}$ ($\dot{M}_{\rm B}$ is Bondi accretion rate), where the electron density profile, $n_{\rm e}\propto r^{\sim -1}$, in the accretion flow is consistent with that determined from X-ray observation inside the Bondi radius and recent numerical simulations.

AMS-02 recently published its lepton spectra measurement. The results show that the positron fraction no longer increases above $\sim$200 GeV. The aim of this work is to investigate the possibility that the excess of positron fraction is due to pulsars. Nearby known pulsars from ATNF catalogue are considered as a possible primary positron source of the high energy positrons. We find that the pulsars with age $T\simeq (0.45\sim4.5)\times10^{5}$ yr and distance $d<0.5$ kpc can explain the behavior of positron fraction of AMS-02 in the range of high energy. We show that each of the four pulsars --- Geminga, J1741-2054, Monogem and J0942-5552 --- is able to be a single source satisfying all considered physical requirements. We also discuss the possibility that these high energy $e^{\pm}$ are from multiple pulsars. The multiple pulsars contribution predicts a positron fraction with some structures at higher energies.

We supplement the minimal supersymmetric standard model (MSSM) with vector-like copies of standard model particles. Such 4th generation particles can raise the Higgs boson mass to the observed value without requiring very heavy superpartners, improving naturalness and the prospects for discovering supersymmetry at the LHC. Here we show that these new particles are also motivated cosmologically: in the MSSM, pure Bino dark matter typically overcloses the Universe, but 4th generation particles open up new annihilation channels, allowing Binos to have the correct thermal relic density without resonances or co-annihilation. We show that this can be done in a sizable region of parameter space while preserving gauge coupling unification and satisfying constraints from collider, Higgs, precision electroweak, and flavor physics.

There is strong evidence in favor of the idea that dark matter is self interacting, with the cross section-to-mass ratio $\sigma / m \sim 1\,\mathrm{cm^2/g} \sim 1\,\mathrm{barn/GeV}$. We show that viable models of dark matter with this large cross section are straightforwardly realized with non-Abelian hidden sectors. In the simplest of such models, the hidden sector is a pure gauge theory, and the dark matter is composed of hidden glueballs with a mass around $100\,\mathrm{MeV}$. Alternatively, the hidden sector may be a supersymmetric pure gauge theory with a $\sim 10\,\mathrm{TeV}$ gluino thermal relic. In this case, the dark matter is largely composed of glueballinos that strongly self interact through the exchange of light glueballs. We present a unified framework that realizes both of these possibilities in anomaly-mediated supersymmetry breaking, where, depending on a few model parameters, the dark matter may be composed of hidden glueballinos, hidden glueballs, or a mixture of the two. These models provide simple examples of multicomponent dark matter, have interesting implications for particle physics and cosmology, and include cases where a subdominant component of dark matter may be extremely strongly self interacting, with interesting astrophysical consequences.

We consider a simple supersymmetric hidden sector: pure SU(N) gauge theory. Dark matter is made up of hidden glueballinos with mass $m_X$ and hidden glueballs with mass near the confinement scale $\Lambda$. For $m_X \sim 1\,\text{TeV}$ and $\Lambda \sim 100\,\text{MeV}$, the glueballinos freeze out with the correct relic density and self-interact through glueball exchange to resolve small-scale structure puzzles. An immediate consequence is that the glueballino spectrum has a hyperfine splitting of order $\Lambda^2 / m_X \sim 10\,\text{keV}$. We show that the radiative decays of the excited state can explain the observed 3.5 keV X-ray line signal from clusters of galaxies, Andromeda, and the Milky Way.

In this report we summarize the many dark matter searches currently being pursued through four complementary approaches: direct detection, indirect detection, collider experiments, and astrophysical probes. The essential features of broad classes of experiments are described, each with their own strengths and weaknesses. The complementarity of the different dark matter searches is discussed qualitatively and illustrated quantitatively in two simple theoretical frameworks. Our primary conclusion is that the diversity of possible dark matter candidates requires a balanced program drawing from all four approaches.

Astronomers have been using the measured luminosity to estimate the {\em luminous mass} of stars, based on empirically established mass-to-light ratio which seems to be only applicable to a special class of stars---the main-sequence stars---with still considerable uncertainties. Another basic tool to determine the mass of a system of stars or galaxies comes from the study of their motion, as Newton demonstrated with his law of gravitation, which yields the {\em gravitational mass}. Because the luminous mass can at best only represent a portion of the gravitational mass, finding the luminous mass to be different or less than the gravitational mass should not be surprising. Using such an apparent discrepancy as a compelling evidence for the so-called dark matter, which has been believed to possess mysterious nonbaryonic properties and present a dominant amount in galaxies and the universe, seems to be too far a stretch when seriously examining the facts and uncertainties in the measurement techniques. In our opinion, a galaxy with star type distribution varying from its center to edge may have a mass-to-light ratio varying accordingly. With the thin-disk model computations based on measured rotation curves, we found that most galaxies have a typical mass density profile that peaks at the galactic center and decreases rapidly within $\sim 5%$ of the cut-off radius, and then declines nearly exponentially toward the edge. The predicted mass density in the Galactic disk is reasonably within the reported range of that observed in interstellar medium. This leads us to believe that ordinary baryonic matter can be sufficient for supporting the observed galactic rotation curves; speculation of large amount of non-baryonic matter may be based on an ill-conceived discrepancy between gravitational mass and luminous mass which appears to be unjustified.

An accurate computational method is presented to determine the mass distribution in a rotating thin-disk galaxy from given rotation curve by applying Newtonian dynamics for an axisymmetrically rotating thin disk of finite size with or without a central spherical bulge. The governing integral equation for mass distribution, resulting from the balance between the Newtonian gravitational force and centrifugal force due to rotation at every point on the disk, is transformed via a boundary-element method into a linear algebra matrix equation that can be solved numerically for rotation curves with a wide range of shapes. To illustrate the effectiveness of this computational method, mass distributions in several mature spiral galaxies are determined from their measured rotation curves. All the surface mass density profiles predicted by our model exhibit approximately a common exponential law of decay, qualitatively consistent with the observed surface brightness distributions. When a central spherical bulge is present, the total galactic mass increases only slightly but the mass distribution in the galaxy is altered in such a way that the periphery mass density is reduced while more mass appears toward the galactic center. By extending the computational domain beyond the galactic edge, we can determine rotation velocity outside the cut-off radius which appears to continuously decrease and gradually approach the Keplerian rotation velocity out over twice the cut-off radius. Am examination of the circular orbit stability suggests that galaxies with flat or increasing rotation velocities with radius are more stable than those with decreasing rotation velocities especially in the region near the galactic edge. Our results demonstrate the fact that Newtonian dynamics can be adequate for describing the observed rotation behavior of mature spiral galaxies.

A passive and multi-channel microwave sounder onboard Chang'E-2 orbiter has successfully performed microwave observation of the lunar surface and subsurface structure. Compared with Chang'E-1 orbiter, Chang'E-2 orbiter obtained more accurate and comprehensive microwave brightness temperature data which is helpful for further research. Since there is a close relationship between microwave brightness temperature data and some related properties of the lunar regolith, such as the thickness, temperature and dielectric constant, so precise and high resolution brightness temperature is necessary for such research. However, through the detection mechanism of the microwave sounder, the brightness temperature data acquired from the microwave sounder is weighted by the antenna radiation pattern, so the data is the convolution of the antenna radiation pattern and the lunar brightness temperature. In order to obtain the real lunar brightness temperature, a deconvolution method is needed. The aim of this paper is to solve the problem in performing deconvolution of the lunar brightness temperature. In this study, we introduce the maximum entropy method(MEM) to process the brightness temperature data and achieve excellent results. The paper mainly includes the following aspects: firstly, we introduce the principle of the MEM, secondly, through a series of simulations, the MEM has been verified an efficient deconvolution method, thirdly, the MEM is used to process the Chang'E-2 microwave data and the results are significant.

Fluctuations in the intensity and polarization of the cosmic microwave background (CMB) and the large-scale distribution of matter in the universe each contain clues about the nature of the earliest moments of time. The next generation of CMB and large-scale structure (LSS) experiments are poised to test the leading paradigm for these earliest moments---the theory of cosmic inflation---and to detect the imprints of the inflationary epoch, thereby dramatically increasing our understanding of fundamental physics and the early universe. A future CMB experiment with sufficient angular resolution and frequency coverage that surveys at least 1% of the sky to a depth of 1 uK-arcmin can deliver a constraint on the tensor-to-scalar ratio that will either result in a 5-sigma measurement of the energy scale of inflation or rule out all large-field inflation models, even in the presence of foregrounds and the gravitational lensing B-mode signal. LSS experiments, particularly spectroscopic surveys such as the Dark Energy Spectroscopic Instrument, will complement the CMB effort by improving current constraints on running of the spectral index by up to a factor of four, improving constraints on curvature by a factor of ten, and providing non-Gaussianity constraints that are competitive with the current CMB bounds.

These reports present the results of the 2013 Community Summer Study of the APS Division of Particles and Fields ("Snowmass 2013") on the future program of particle physics in the U.S. Chapter 4, on the Cosmic Frontier, discusses the program of research relevant to cosmology and the early universe. This area includes the study of dark matter and the search for its particle nature, the study of dark energy and inflation, and cosmic probes of fundamental symmetries.

These reports present the results of the 2013 Community Summer Study of the APS Division of Particles and Fields ("Snowmass 2013") on the future program of particle physics in the U.S. Chapter 1 contains the Executive Summary and the summaries of the reports of the nine working groups.

In this Report we discuss the four complementary searches for the identity of dark matter: direct detection experiments that look for dark matter interacting in the lab, indirect detection experiments that connect lab signals to dark matter in our own and other galaxies, collider experiments that elucidate the particle properties of dark matter, and astrophysical probes sensitive to non-gravitational interactions of dark matter. The complementarity among the different dark matter searches is discussed qualitatively and illustrated quantitatively in several theoretical scenarios. Our primary conclusion is that the diversity of possible dark matter candidates requires a balanced program based on all four of those approaches.

Isospin-violating dark matter (IVDM) generalizes the standard spin-independent scattering parameter space by introducing one additional parameter, the neutron-to-proton coupling ratio f_n/f_p. In IVDM the implications of direct detection experiments can be altered significantly. We review the motivations for considering IVDM and present benchmark models that illustrate some of the qualitatively different possibilities. IVDM strongly motivates the use of a variety of target nuclei in direct detection experiments.

We determine the prospects for direct and indirect detection of thermal relic neutralinos in supersymmetric theories with multi-TeV squarks and sleptons. We consider the concrete example of the focus point region of minimal supergravity, but our results are generically valid for all models with decoupled scalars and mixed Bino-Higgsino or Higgsino-like dark matter. We determine the parameter space consistent with a 125 GeV Higgs boson including 3-loop corrections in the calculation of the Higgs mass. These corrections increase m_h by 1-3 GeV, lowering the preferred scalar mass scale and decreasing the fine-tuning measure in these scenarios. We then systematically examine prospects for dark matter direct and indirect detection. Direct detection constraints do not exclude these models, especially for \mu < 0. At the same time, the scenario generically predicts spin-independent signals just beyond current bounds. We also consider indirect detection with neutrinos, gamma rays, anti-protons, and anti-deuterons. Current IceCube neutrino constraints are competitive with direct detection, implying bright prospects for complementary searches with both direct and indirect detection.

We consider models of xenophobic dark matter, in which isospin-violating dark matter-nucleon interactions significantly degrade the response of xenon direct detection experiments. For models of near-maximal xenophobia, with neutron-to-proton coupling ratio $f_n / f_p \approx -0.64$, and dark matter mass near 8 GeV, the regions of interest for CoGeNT and CDMS-Si and the region of interest identified by Collar and Fields in CDMS-Ge data can be brought into agreement. This model may be tested in future direct, indirect, and collider searches. Interestingly, because the natural isotope abundance of xenon implies that xenophobia has its limits, we find that this xenophobic model may be probed in the near future by xenon experiments. Near-future data from the LHC and Fermi-LAT may also provide interesting alternative probes of xenophobic dark matter.

For decades, the unnaturalness of the weak scale has been the dominant problem motivating new particle physics, and weak-scale supersymmetry has been the dominant proposed solution. This paradigm is now being challenged by a wealth of experimental data. In this review, we begin by recalling the theoretical motivations for weak-scale supersymmetry, including the gauge hierarchy problem, grand unification, and WIMP dark matter, and their implications for superpartner masses. These are set against the leading constraints on supersymmetry from collider searches, the Higgs boson mass, and low-energy constraints on flavor and CP violation. We then critically examine attempts to quantify naturalness in supersymmetry, stressing the many subjective choices that impact the results both quantitatively and qualitatively. Finally, we survey various proposals for natural supersymmetric models, including effective supersymmetry, focus point supersymmetry, compressed supersymmetry, and R-parity-violating supersymmetry, and summarize their key features, current status, and implications for future experiments.

In the early years, cosmic rays contributed essentially to particle physics through the discovery of new particles. Will history repeat itself? As with the discovery of the charged pion, the recent discovery of a Higgs-like boson may portend a rich new set of particles within reach of current and near future experiments. These may be discovered and studied by cosmic rays through the indirect detection of dark matter.

Recent indications of a 125 GeV Higgs boson are challenging for gauge-mediated supersymmetry breaking (GMSB), since radiative contributions to the Higgs boson mass are not enhanced by significant stop mixing. This challenge should not be considered in isolation, however, as GMSB also generically suffers from two other problems: unsuppressed electric dipole moments and the absence of an attractive dark matter candidate. We show that all of these problems may be simultaneously solved by considering heavy superpartners, without extra fields or modified cosmology. Multi-TeV sfermions suppress the EDMs and raise the Higgs mass, and the dark matter problem is solved by Goldilocks cosmology, in which TeV neutralinos decay to GeV gravitinos that are simultaneously light enough to solve the flavor problem and heavy enough to be all of dark matter. The implications for collider searches and direct and indirect dark matter detection are sobering, but EDMs are expected near their current bounds, and the resulting non-thermal gravitino dark matter is necessarily warm, with testable cosmological implications.

We present a model with dark matter in an anomaly-mediated supersymmetry breaking hidden sector with a U(1)xU(1) gauge symmetry. The symmetries of the model stabilize the dark matter and forbid the introduction of new mass parameters. As a result, the thermal relic density is completely determined by the gravitino mass and dimensionless couplings. Assuming non-hierarchical couplings, the thermal relic density is ~ 0.1, independent of the dark matter's mass and interaction strength, realizing the WIMPless miracle. The model has several striking features. For particle physics, stability of the dark matter is completely consistent with R-parity violation in the visible sector, with implications for superpartner collider signatures; also the thermal relic's mass may be ~ 10 GeV or lighter, which is of interest given recent direct detection results. Interesting astrophysical signatures are dark matter self-interactions through a long-range force, and massless hidden photons and fermions that contribute to the number of relativistic degrees of freedom at BBN and CMB. The latter are particularly interesting, given current indications for extra degrees of freedom and near future results from the Planck observatory.

Searches for dark matter scattering off nuclei are typically compared assuming that the dark matter's spin-independent couplings are identical for protons and neutrons. This assumption is neither innocuous nor well motivated. We consider isospin-violating dark matter (IVDM) with one extra parameter, the ratio of neutron to proton couplings, and include the isotope distribution for each detector. For a single choice of the coupling ratio, the DAMA and CoGeNT signals are consistent with each other and with current XENON constraints, and they unambiguously predict near future signals at XENON and CRESST. We provide a quark-level realization of IVDM as WIMPless dark matter that is consistent with all collider and low-energy bounds.

In anomaly-mediated supersymmetry breaking, superpartners in a hidden sector have masses that are proportional to couplings squared, and so naturally freeze out with the desired dark matter relic density for a large range of masses. We present an extremely simple realization of this possibility, with WIMPless dark matter arising from a hidden sector that is supersymmetric QED with N_F flavors. Dark matter is multi-component, composed of hidden leptons and sleptons with masses anywhere from 10 GeV to 10 TeV, and hidden photons provide the thermal bath. The dark matter self-interacts through hidden sector Coulomb scatterings that are potentially observable. In addition, the hidden photon contribution to the number of relativistic degrees of freedom is in the range \Delta N_eff ~ 0 - 2, and, if the hidden and visible sectors were initially in thermal contact, the model predicts \Delta N_eff ~ 0.2 - 0.4. Data already taken by Planck may provide evidence of such deviations.

By numerically solving the mass distribution in a rotating disk based on Newton's laws of motion and gravitation, we demonstrate that the observed flat rotation curves for most spiral galaxies correspond to exponentially decreasing mass density from galactic center for the most of the part except within the central core and near periphery edge. Hence, we believe the galaxies described with our model are consistent with that seen through the eyes of Newton. Although Newton's laws and Kepler's laws seem to yield the same results when they are applied to the planets in the solar system, they are shown to lead to quite different results when describing the stellar dynamics in disk galaxies. This is because that Keplerian dynamics may be equivalent to Newtonian dynamics for only special circumstances, but not generally for all the cases. Thus, the conclusions drawn from calculations based on Keplerian dynamics are often likely to be erroneous when used to describe rotating disk galaxies.

We present an efficient, robust computational method for modeling the Newtonian dynamics for rotation curve analysis of thin-disk galaxies. For a disk galaxy with a typical flat rotation curve, our modeling results show that the surface mass density monotonically decreases from the galactic center toward periphery, according to Newtonian dynamics. In a large portion of the galaxy, the surface mass density follows an approximately exponential law of decay with respect to the galactic radial coordinate. Yet the radial scale length for the surface mass density seems to be generally larger than that of the measured brightness distribution, suggesting an increasing mass-to-light ratio with the radial distance in a disk galaxy. In a nondimensionalized form, our mathematical system contains a dimensionless parameter which we call the "galactic rotation number" that represents the gross ratio of centrifugal force and gravitational force. The value of this galactic rotation number is determined as part of the numerical solution. Through a systematic computational analysis, we have illustrated that this galactic rotation number remains within $\pm 10%$ of 1.70 for a wide variety of rotation curves. This implies that the total mass in a disk galaxy is proportional to $V_0^2\,R_g$, with $V_0$ denoting the characteristic rotation velocity and $R_g$ the radius of the galactic disk. The predicted total galactic mass of the Milky Way is in good agreement with the star-count data.

WIMPless dark matter provides a framework in which dark matter particles with a wide range of masses naturally have the correct thermal relic density. We show that WIMPless dark matter with mass around 2-10 GeV can explain the annual modulation observed by the DAMA experiment without violating the constraints of other dark matter searches. This explanation implies distinctive and promising signals for other direct detection experiments, GLAST, and the LHC.

The direct detection of dark matter through its elastic scattering off nucleons is among the most promising methods for establishing the particle identity of dark matter. The current bound on the spin-independent scattering cross section is sigma^SI < 10 zb for dark matter masses m_chi ~ 100 GeV, with improved sensitivities expected soon. We examine the implications of this progress for neutralino dark matter. We work in a supersymmetric framework well-suited to dark matter studies that is simple and transparent, with models defined in terms of four weak-scale parameters. We first show that robust constraints on electric dipole moments motivate large sfermion masses mtilde > 1 TeV, effectively decoupling squarks and sleptons from neutralino dark matter phenomenology. In this case, we find characteristic cross sections in the narrow range 1 zb < sigma^SI < 40 zb for m_chi > 70 GeV. As sfermion masses are lowered to near their experimental limit mtilde ~ 400 GeV, the upper and lower limits of this range are extended, but only by factors of around two, and the lower limit is not significantly altered by relaxing many particle physics assumptions, varying the strange quark content of the nucleon, including the effects of galactic small-scale structure, or assuming other components of dark matter. Experiments are therefore rapidly entering the heart of dark matter-favored supersymmetry parameter space. If no signal is seen, supersymmetric models must contain some level of fine-tuning, and we identify and analyze several possibilities. Barring large cancellations, however, in a large and generic class of models, if thermal relic neutralinos are a significant component of dark matter, experiments will discover them as they probe down to the zeptobarn scale.

In anomaly-mediated supersymmetry breaking (AMSB) models, superpartner masses are proportional to couplings squared. Their hidden sectors therefore naturally contain WIMPless dark matter, particles whose thermal relic abundance is guaranteed to be of the correct size, even though they are not weakly-interacting massive particles (WIMPs). We study viable dark matter candidates in WIMPless AMSB models with non-Abelian hidden sectors and highlight unusual possibilities that emerge in even the simplest models. In one example with a pure SU(N) hidden sector, stable hidden gluinos freeze out with the correct relic density, but have an extremely low, but natural, confinement scale, providing a framework for self-interacting dark matter. In another simple scenario, hidden gluinos freeze out and decay to visible Winos with the correct relic density, and hidden glueballs may either be stable, providing a natural framework for mixed cold-hot dark matter, or may decay, yielding astrophysical signals. Last, we present a model with light hidden pions that may be tested with improved constraints on the number of non-relativistic degrees of freedom. All of these scenarios are defined by a small number of parameters, are consistent with gauge coupling unification, preserve the beautiful connection between the weak scale and the observed dark matter relic density, and are natural, with relatively light visible superpartners. We conclude with comments on interesting future directions.

The annihilation cross section of thermal relic dark matter determines both its relic density and indirect detection signals. We determine how large indirect signals may be in scenarios with Sommerfeld-enhanced annihilation, subject to the constraint that the dark matter has the correct relic density. This work refines our previous analysis through detailed treatments of resonant Sommerfeld enhancement and the effect of Sommerfeld enhancement on freeze out. Sommerfeld enhancements raise many interesting issues in the freeze out calculation, and we find that the cutoff of resonant enhancement, the equilibration of force carriers, the temperature of kinetic decoupling, and the efficiency of self-interactions for preserving thermal velocity distributions all play a role. These effects may have striking consequences; for example, for resonantly-enhanced Sommerfeld annihilation, dark matter freezes out but may then chemically recouple, implying highly suppressed indirect signals, in contrast to naive expectations. In the minimal scenario with standard astrophysical assumptions, and tuning all parameters to maximize the signal, we find that, for force-carrier mass m_phi = 250 MeV and dark matter masses m_X = 0.1, 0.3, and 1 TeV, the maximal Sommerfeld enhancement factors are S_eff = 7, 30, and 90, respectively. Such boosts are too small to explain both the PAMELA and Fermi excesses. Non-minimal models may require smaller boosts, but the bounds on S_eff could also be more stringent, and dedicated freeze out analyses are required. For concreteness, we focus on 4 mu final states, but we also discuss 4 e and other modes, deviations from standard astrophysical assumptions and non-minimal particle physics models, and we outline the steps required to determine if such considerations may lead to a self-consistent explanation of the PAMELA or Fermi excesses.

Light gravitinos, with mass in the eV to MeV range, are well-motivated in particle physics, but their status as dark-matter candidates is muddled by early-Universe uncertainties. We investigate how upcoming data from colliders may clarify this picture. Light gravitinos are produced primarily in the decays of the next-to-lightest supersymmetric particle, resulting in spectacular signals, including di-photons, delayed and non-pointing photons, kinked charged tracks, and heavy metastable charged particles. We find that the Tevatron with 20/fb and the 7 TeV LHC with 1/fb may both see evidence for hundreds of light-gravitino events. Remarkably, this collider data is also well suited to distinguish between currently viable light-gravitino scenarios, with striking implications for structure formation, inflation, and other early-Universe cosmology.

The identity of dark matter is a question of central importance in both astrophysics and particle physics. In the past, the leading particle candidates were cold and collisionless, and typically predicted missing energy signals at particle colliders. However, recent progress has greatly expanded the list of well-motivated candidates and the possible signatures of dark matter. This review begins with a brief summary of the standard model of particle physics and its outstanding problems. We then discuss several dark matter candidates motivated by these problems, including WIMPs, superWIMPs, light gravitinos, hidden dark matter, sterile neutrinos, and axions. For each of these, we critically examine the particle physics motivations and present their expected production mechanisms, basic properties, and implications for direct and indirect detection, particle colliders, and astrophysical observations. Upcoming experiments will discover or exclude many of these candidates, and progress may open up an era of unprecedented synergy between studies of the largest and smallest observable length scales.

Dark matter with Sommerfeld-enhanced annihilation has been proposed to explain observed cosmic ray positron excesses in the 10 GeV to TeV energy range. We show that the required enhancement implies thermal relic densities that are too small to be all of dark matter. We also show that the dark matter is sufficiently self-interacting that observations of elliptical galactic dark matter halos exclude large Sommerfeld enhancement for light force carriers. Resonant Sommerfeld enhancement does not modify these conclusions, and the astrophysical boosts required to resolve these discrepancies are disfavored, especially when significant self-interactions suppress halo substructure.

We determine the prospects for finding dark matter at the Tevatron and LHC through the production of exotic 4th generation quarks T' that decay through T' \to t X, where X is dark matter. The resulting signal of t \bar{t} + \met has not previously been considered in searches for 4th generation quarks, but there are both general and specific dark matter motivations for this signal, and with slight modifications, this analysis applies to any scenario where invisible particles are produced in association with top quarks. Current direct and indirect bounds on such exotic quarks restrict their masses to be between 300 and 600 GeV, and the dark matter's mass may be anywhere below m_T'. We simulate the signal and main backgrounds with MadGraph/MadEvent-Pythia-PGS4. For the Tevatron, we find that an integrated luminosity of 20 fb^-1 will allow 3\sigma discovery up to m_T' = 400 GeV and 95% exclusion up to m_T' = 455 GeV. For the 10 TeV LHC with 300 pb^-1, the discovery and exclusion sensitivities rise to 490 GeV and 600 GeV. These scenarios are therefore among the most promising for dark matter at colliders. Perhaps most interestingly, we find that dark matter models that can explain results from the DAMA, CDMS and CoGeNT Collaborations can be tested with high statistical significance using data already collected at the Tevatron and have extraordinarily promising implications for early runs of the LHC.

Non-WIMP dark matter candidates include particles motivated by minimality, candidates motivated by experimental anomalies, and exotic possibilities motivated primarily by the desire of clever iconoclasts to highlight how truly ignorant we are about the nature of dark matter. In this review, I discuss candidates that are not WIMPs, but nevertheless share the same theoretical motivations as WIMPs and also naturally have the correct relic density. There are two classes: superWIMP dark matter, where the desired relic density is inherited through decays, and WIMPless dark matter, where the dark matter's mass and couplings scale together to maintain the desired thermal relic density.

I review recent developments in the direct and indirect detection of dark matter and new candidates beyond the WIMP paradigm.

Can dark matter be stabilized by charge conservation, just as the electron is in the standard model? We examine the possibility that dark matter is hidden, that is, neutral under all standard model gauge interactions, but charged under an exact U(1) gauge symmetry of the hidden sector. Such candidates are predicted in WIMPless models, supersymmetric models in which hidden dark matter has the desired thermal relic density for a wide range of masses. Hidden charged dark matter has many novel properties not shared by neutral dark matter: (1) bound state formation and Sommerfeld-enhanced annihilation after chemical freeze out may reduce its relic density, (2) similar effects greatly enhance dark matter annihilation in protohalos at redshifts of z ~ 30, (3) Compton scattering off hidden photons delays kinetic decoupling, suppressing small scale structure, and (4) Rutherford scattering makes such dark matter self-interacting and collisional, potentially impacting properties of the Bullet Cluster and the observed morphology of galactic halos. We analyze all of these effects in a WIMPless model in which the hidden sector is a simplified version of the minimal supersymmetric standard model and the dark matter is a hidden sector stau. We find that charged hidden dark matter is viable and consistent with the correct relic density for reasonable model parameters and dark matter masses in the range 1 GeV < m_X < 10 TeV. At the same time, in the preferred range of parameters, this model predicts cores in the dark matter halos of small galaxies and other halo properties that may be within the reach of future observations. These models therefore provide a viable and well-motivated framework for collisional dark matter with Sommerfeld enhancement, with novel implications for astrophysics and dark matter searches.

We consider the prospects for testing the dark matter interpretation of the DAMA/LIBRA signal with the Super-Kamiokande experiment. The DAMA/LIBRA signal favors dark matter with low mass and high scattering cross section. We show that these characteristics imply that the scattering cross section that enters the DAMA/LIBRA event rate determines the annihilation rate probed by Super-Kamiokande. Current limits from Super-Kamiokande through-going events do not test the DAMA/LIBRA favored region. We show, however, that upcoming analyses including fully-contained events with sensitivity to dark matter masses from 5 to 10 GeV may corroborate the DAMA/LIBRA signal. We conclude by considering three specific dark matter candidates, neutralinos, WIMPless dark matter, and mirror dark matter, which illustrate the various model-dependent assumptions entering our analysis.

We propose that dark matter is composed of particles that naturally have the correct thermal relic density, but have neither weak-scale masses nor weak interactions. These WIMPless models emerge naturally from gauge-mediated supersymmetry breaking, where they elegantly solve the dark matter problem. The framework accommodates single or multiple component dark matter, dark matter masses from 10 MeV to 10 TeV, and interaction strengths from gravitational to strong. These candidates enhance many direct and indirect signals relative to WIMPs and have qualitatively new implications for dark matter searches and cosmological implications for colliders.

Dark matter may be hidden, with no standard model gauge interactions. At the same time, in WIMPless models with hidden matter masses proportional to hidden gauge couplings squared, the hidden dark matter's thermal relic density may naturally be in the right range, preserving the key quantitative virtue of WIMPs. We consider this possibility in detail. We first determine model-independent constraints on hidden sectors from Big Bang nucleosynthesis and the cosmic microwave background. Contrary to conventional wisdom, large hidden sectors are easily accommodated. A flavour-free version of the standard model is allowed if the hidden sector is just 30% colder than the observable sector after reheating. Alternatively, if the hidden sector contains a 1-generation version of the standard model with characteristic mass scale below 1 MeV, even identical reheating temperatures are allowed. We then analyze hidden sector freezeout in detail for a concrete model, solving the Boltzmann equation numerically and understanding the results from both observable and hidden sector points of view. We find that WIMPless dark matter indeed obtains the correct relic density for masses in the range keV < m_X < TeV. The upper bound results from the requirement of perturbativity, and the lower bound assumes that the observable and hidden sectors reheat to the same temperature and is raised to the MeV scale if the hidden sector is 10 times colder. WIMPless dark matter therefore generalizes the WIMP paradigm to the largest mass range possible for viable thermal relics and provides a unified framework for exploring dark matter signals across nine orders of magnitude in dark matter mass.

In the coming year, the Large Hadron Collider will begin colliding protons at energies nearly an order of magnitude beyond the current frontier. The LHC will, of course, provide unprecedented opportunities to discover new particle physics. Less well-known, however, is that the LHC may also provide insights about gravity and the early universe. I review some of these connections, focusing on the topics of dark matter and dark energy, and highlight outstanding prospects for breakthroughs at the interface of particle physics and cosmology.

Observations reveal that mature spiral galaxies consist of stars, gases and plasma approximately distributed in a thin disk of circular shape, usually with a central bulge. The rotation velocities quickly increase from the galactic center and then achieve a constant velocity from the core to the periphery. The basic dynamic behavior of a mature spiral galaxy, such as the Milky Way, is well described by simple models balancing Newtonian gravitational forces against the centrifugal forces associated with a rotating thin axisymmetric disk. In this research, we investigate the effects of adding central bulges to thin disk gravitational models. Even with the addition of substantial central bulges, all the critical essential features of our thin disk gravitational models are preserved. (1) Balancing Newtonian gravitational and centrifugal forces at every point within the disk yields computed radial mass distributions that describe the measured rotation velocity profiles of mature spiral galaxies successfully. (2) There is no need for gravity deviations or ``massive peripheral spherical halos of mysterious Dark Matter''. (3) The calculated total galactic masses are in good agreement with star count data. (4) The addition of central bulges increases the calculated total galactic masses, possibly more consistent with the presence of galactic gases, dust, grains, lumps, planets and plasma in addition to stars. (5) Compared with the light distribution, our mass distributions within the disk are larger out toward the galactic periphery which is cooler with lower opactiy/emissivity (and thus darker). This is apparent from edge-on views of galaxies which display a dark disk-line against a much brighter galactic halo.

For mature spiral galaxies, the rotation velocities quickly increase from the galactic center and achieve a constant velocity from the core to the periphery. This dynamic behavior is described by models balancing Newtonian gravitational and centrifugal forces in rotating thin axisymmetric disks. Freeman's disk assumes a mass density decreasing exponentially with radius which correctly produces rotational velocities which increase from the galactic center to a maximum near the outer core, but then decreases out to the periphery contrary to measurements. Mestel's disk assumes a mass distribution decreasing more slowly (inversely with radius) that yields a constant rotational velocity across the entire disk, but has an unrealistic central mass singularity and does not describe the core rotation properly. Thus combine the Freeman and Mestel disks to utilize their strengths and eliminate their deficiencies. Utilize the Freeman formula for the central core, and the Mestel formula beyond the core to the galactic rim. This combined model produces rotation curves comparable to the measurements of mature spiral galaxies. For a more general thin-disk model, we develop an alternative computational method to solve for the mass density distributions for various measured rotation curves We compute radial mass densities that balance the Newtonian gravitational and centrifugal forces at every point. The computational solutions show mass densities which decrease approximately exponentially in the central core (similar to Freeman), and transition to a slower inverse radial decrease (similar to Mestel) to the periphery. Thus these diverse approaches yield similar results and are mutaully self-consistent with Newtonian gravity/dynamics.

The measured rotation velocity profiles of mature spiral galaxies are successfully described with a gravitational model consisting of a thin axisymmetric disk of finte radius. The disk is assumed uniformly thin but with variable radial mass density. The governing integral equation is based on mechanical balance between Newtonian gravitational and centrifugal forces (due to galaxy rotation) at each and every point in a finite set of concentric rings. The nondimensionalized mathematical system contains a dimensionless parameter we call ``galactic rotation parameter'' which concisely crystallizes perspective. Computational solutions are obtained for the radial mass distributions that satisfy the measured rotational velocity profiles. Together with a constraint equation for mass conservation, the galactic rotational parameter is also determined from which the total galactic mass is calculated from measured galactic radii and maximum rotation velocities. These calculated total galactic masses are in good agreement with data. Our deduced exponentially decreasing mass distributions in the central galactic core are in agreement with almost all others. However our mass distributions differ toward the galactic periphery with more ordinary baryonic mass in these outer disk regions which are cooler with lower opactiy/emissivity (and thus darker).

Neutralino dark matter is well motivated, but also suffers from two shortcomings: it requires gravity-mediated supersymmetry breaking, which generically violates flavor constraints, and its thermal relic density \Omega is typically too large. We propose a simple solution to both problems: neutralinos freezeout with \Omega ~10-100, but then decay to ~1 GeV gravitinos, which are simultaneously light enough to satisfy flavor constraints and heavy enough to be all of dark matter. This scenario is naturally realized in high-scale gauge-mediation models, ameliorates small scale structure problems, and implies that ``cosmologically excluded'' models may, in fact, be cosmologically preferred.

Dark matter particles need not be completely stable, and in fact they may be decaying now. We consider this possibility in the frameworks of universal extra dimensions and supersymmetry with very late decays of WIMPs to Kaluza-Klein gravitons and gravitinos. The diffuse photon background is a sensitive probe, even for lifetimes far greater than the age of the Universe. Remarkably, both the energy spectrum and flux of the observed MeV gamma ray excess may be naturally explained by decaying dark matter with MeV mass splittings. Future observations of continuum and line photon fluxes will test this explanation and may provide novel constraints on cosmological parameters.

The instability of dark matter may produce visible signals in the spectrum of cosmic gamma-rays. We consider this possibility in frameworks with additional spatial dimensions and supersymmetry. Examples of particles include superweakly-interacting massive particles such as gravitinos in supersymmetry models, the lightest Kaluza-Klein (KK) state in models with universal extra dimensions, and weakly-interacting massive particles such as branons in flexible brane-worlds.

Minimal universal extra dimensions (mUED) is often thought to predict that the lightest Kaluza-Klein particle (LKP) is the Kaluza-Klein gauge boson B^1, leading to conventional missing energy signals at colliders and WIMP dark matter. In fact, the implications of mUED are far richer: the B^1, charged Higgs boson H^{\pm 1}, and graviton G^1 are all possible LKPs, leading to many different ``phases'' with distinct signatures. Considering the complete phase diagram, we find predictions for charged or neutral particles with decay lengths of microns to tens of meters; WIMP, superWIMP, or charged relic particles; metastable particles with lifetimes of the order of or in excess of the age of the universe; and scenarios combining two or more of these phenomena. In the cosmologically preferred region, the Higgs boson mass is between 180 and 245 GeV, the LKP mass is between 810 and 1400 GeV, and the maximal splitting between first Kaluza-Klein modes is less than 320 GeV. This region predicts a variety of exotic collider signals, such as slow charged particles, displaced vertices, tracks with non-vanishing impact parameters, track kinks, and even vanishing charged tracks, all of which provide early discovery possibilities at the Large Hadron Collider.

We summarize the main features of several dark matter candidates in extra-dimensional theories. In particular, we review Kaluza-Klein (KK) gravitons in universal extra dimensions and branons in brane-world models. KK gravitons are superWIMP (superweakly-interacting massive particle) dark matter, and branons are WIMP (weakly-interacting massive particle) dark matter. Both dark matter candidates are naturally produced in the correct amount to form much or all of dark matter.

Ultra-high energy cosmic rays and neutrinos probe energies far above the weak scale. Their usefulness might appear to be limited by astrophysical uncertainties; however, by simultaneously considering up- and down-going events, one may disentangle particle physics from astrophysics. We show that present data from the AMANDA experiment in the South Pole ice already imply an upper bound on neutrino cross sections at energy scales that will likely never be probed at man-made accelerators. The existing data also place an upper limit on the neutrino flux valid for any neutrino cross section. In the future, similar analyses of IceCube data will constrain neutrino properties and fluxes at the O(10%) level.

Recent breakthroughs in cosmology reveal that a quarter of the Universe is composed of dark matter, but the microscopic identity of dark matter remains a deep mystery. I review recent progress in resolving this puzzle, focusing on two well-motivated classes of dark matter candidates: WIMPs and superWIMPs. These possibilities have similar motivations: they exist in the same well-motivated particle physics models, the observed dark matter relic density emerges naturally, and dark matter particles have mass around 100 GeV, the energy scale identified as interesting over 70 years ago by Fermi. At the same time, they have widely varying implications for direct and indirect dark matter searches, particle colliders, Big Bang nucleosynthesis, the cosmic microwave background, and halo profiles and structure formation. If WIMPs or superWIMPs are a significant component of dark matter, we will soon be entering a golden era in which dark matter will be studied through diverse probes at the interface of particle physics, astroparticle physics, and cosmology. I outline a program of dark matter studies for each of these scenarios, and discuss the prospects for identifying dark matter in the coming years.

Collisionless, cold dark matter in the form of weakly-interacting massive particles (WIMPs) is well-motivated in particle physics, naturally yields the observed relic density, and successfully explains structure formation on large scales. On small scales, however, it predicts too much power, leading to cuspy halos, dense cores, and large numbers of subhalos, in apparent conflict with observations. We consider superWIMP dark matter, produced with large velocity in late decays at times 10^5 s - 10^8 s. As analyzed by Kaplinghat in a more general setting, we find that superWIMPs have sufficiently large free-streaming lengths and low phase space densities to help resolve small scale structure problems while preserving all of the above-mentioned WIMP virtues.

Recent breakthroughs in cosmology pose questions that require particle physics answers. I review the problems of dark matter, baryogenesis, and dark energy and discuss how particle colliders, particularly the International Linear Collider, may advance our understanding of the contents and evolution of the Universe.

We investigate supergravity models in which the lightest supersymmetric particle (LSP) is a stable gravitino. We assume that the next-lightest supersymmetric particle (NLSP) freezes out with its thermal relic density before decaying to the gravitino at time t ~ 10^4 s - 10^8 s. In contrast to studies that assume a fixed gravitino relic density, the thermal relic density assumption implies upper, not lower, bounds on superpartner masses, with important implications for particle colliders. We consider slepton, sneutrino, and neutralino NLSPs, and determine what superpartner masses are viable in all of these cases, applying CMB and electromagnetic and hadronic BBN constraints to the leading two- and three-body NLSP decays. Hadronic constraints have been neglected previously, but we find that they provide the most stringent constraints in much of the natural parameter space. We then discuss the collider phenomenology of supergravity with a gravitino LSP. We find that colliders may provide important insights to clarify BBN and the thermal history of the Universe below temperatures around 10 GeV and may even provide precise measurements of the gravitino's mass and couplings.

In 2002, Ray Davis and Masatoshi Koshiba were awarded the Nobel Prize in Physics ``for pioneering contributions to astrophysics, in particular for the detection of cosmic neutrinos.'' However, while astronomy has undergone a revolution in understanding by synthesizing data taken at many wavelengths, the universe has only barely been glimpsed in neutrinos, just the Sun and the nearby SN 1987A. An entire universe awaits, and since neutrinos can probe astrophysical objects at densities, energies, and distances that are otherwise inaccessible, the results are expected to be particularly exciting. Similarly, the revolution in quantitative cosmology has heightened the need for very precise tests that depend on the effects of neutrinos, and prominent among them is the search for the effects of neutrino mass, since neutrinos are a small but known component of the dark matter. In this report, we highlight some of the key opportunties for progress in neutrino astrophysics and cosmology, and the implications for other areas of physics.

Dark matter may be composed of superWIMPs, superweakly-interacting massive particles produced in the late decays of other particles. We focus here on the well-motivated supersymmetric example of gravitino LSPs. Gravitino superWIMPs share several virtues with the well-known case of neutralino dark matter: they are present in the same supersymmetric frameworks (supergravity with R-parity conservation) and naturally have the desired relic density. In contrast to neutralinos, however, gravitino superWIMPs are impossible to detect by conventional dark matter searches, may explain an existing discrepancy in Big Bang nucleosynthesis, predict observable distortions in the cosmic microwave background, and imply spectacular signals at future particle colliders.

Cosmology now provides unambiguous, quantitative evidence for new particle physics. I discuss the implications of cosmology for supersymmetry and vice versa. Topics include: motivations for supersymmetry; supersymmetry breaking; dark energy; freeze out and WIMPs; neutralino dark matter; cosmologically preferred regions of minimal supergravity; direct and indirect detection of neutralinos; the DAMA and HEAT signals; inflation and reheating; gravitino dark matter; Big Bang nucleosynthesis; and the cosmic microwave background. I conclude with speculations about the prospects for a microscopic description of the dark universe, stressing the necessity of diverse experiments on both sides of the particle physics/cosmology interface.

Dark matter may be composed of superWIMPs, superweakly-interacting massive particles produced in the late decays of other particles. We focus on the case of gravitinos produced in the late decays of sleptons or sneutrinos and assume they are produced in sufficient numbers to constitute all of non-baryonic dark matter. At leading order, these late decays are two-body and the accompanying energy is electromagnetic. For natural weak-scale parameters, these decays have been shown to satisfy bounds from Big Bang nucleosynthesis and the cosmic microwave background. However, sleptons and sneutrinos may also decay to three-body final states, producing hadronic energy, which is subject to even more stringent nucleosynthesis bounds. We determine the three-body branching fractions and the resulting hadronic energy release. We find that superWIMP gravitino dark matter is viable and determine the gravitino and slepton/sneutrino masses preferred by this solution to the dark matter problem. In passing, we note that hadronic constraints disfavor the possibility of superWIMPs produced by neutralino decays unless the neutralino is photino-like.

The gravitational interactions of elementary particles are suppressed by the Planck scale M_P ~ 10^18 GeV and are typically expected to be far too weak to be probed by experiments. We show that, contrary to conventional wisdom, such interactions may be studied by particle physics experiments in the next few years. As an example, we consider conventional supergravity with a stable gravitino as the lightest supersymmetric particle. The next-lightest supersymmetric particle (NLSP) decays to the gravitino through gravitational interactions after about a year. This lifetime can be measured by stopping NLSPs at colliders and observing their decays. Such studies will yield a measurement of Newton's gravitational constant on unprecedentedly small scales, shed light on dark matter, and provide a window on the early universe.

I review prospects for detecting dark matter in space-based experiments, with an emphasis on recent developments. I propose the ``Martha Stewart criterion'' for identifying dark matter candidates that are particularly worth investigation and focus on three that satisfy it: neutralino dark matter, Kaluza-Klein dark matter, and superWIMP gravitino dark matter.

Black hole production in elementary particle collisions is among the most promising probes of large extra spacetime dimensions. Studies of black holes at particle colliders have assumed that all of the incoming energy is captured in the resulting black hole. We incorporate the inelasticity inherent in such processes and determine the prospects for discovering black holes in colliders and cosmic ray experiments, employing a dynamical model of Hawking evolution. At the Large Hadron Collider, inelasticity reduces rates by factors of 10^3 to 10^6 in the accessible parameter space, moderating, but not eliminating, hopes for black hole discovery. At the Pierre Auger Observatory, rates are suppressed by a factor of 10. We evaluate the impact of cosmic ray observations on collider prospects.

In this talk we discuss aspects of TeV-scale gravitational collapse to black holes and string balls and their subsequent evaporation. Special emphasis is placed on the interplay of the string/black hole correspondence principle. These ideas are then explored in the context of cosmic ray physics. First, the potential for observing showers mediated by black holes or superstring excitations is examined. Next, existing data from neutrino telescopes are used to constrain the parameter space for the unseen dimensions of the universe. Finally, we close with a discussion of future prospects.

In supergravity theories, a natural possibility is that neutralinos or sleptons freeze out at their thermal relic density, but then decay to gravitinos after about a year. The resulting gravitinos are then superWIMPs -- superweakly-interacting massive particles that naturally inherit the desired relic density from late decays of conventional WIMPs. SuperWIMP dark matter escapes all conventional searches. However, the late decays that produce superWIMPs provide new and promising early universe signatures for cold dark matter.

Cold dark matter may be made of superweakly-interacting massive particles, superWIMPs, that naturally inherit the desired relic density from late decays of metastable WIMPs. Well-motivated examples are weak-scale gravitinos in supergravity and Kaluza-Klein gravitons from extra dimensions. These particles are impossible to detect in all dark matter experiments. We find, however, that superWIMP dark matter may be discovered through cosmological signatures from the early universe. In particular, superWIMP dark matter has observable consequences for Big Bang nucleosynthesis and the cosmic microwave background (CMB), and may explain the observed underabundance of 7Li without upsetting the concordance between deuterium and CMB baryometers. We discuss implications for future probes of CMB black body distortions and collider searches for new particles. In the course of this study, we also present a model-independent analysis of entropy production from late-decaying particles in light of WMAP data.

In models of universal extra dimensions, gravity and all standard model fields propagate in the extra dimensions. Previous studies of such models have concentrated on the Kaluza-Klein (KK) partners of standard model particles. Here we determine the properties of the KK gravitons and explore their cosmological implications. We find the lifetimes of decays to KK gravitons, of relevance for the viability of KK gravitons as dark matter. We then discuss the primordial production of KK gravitons after reheating. The existence of a tower of KK graviton states makes such production extremely efficient: for reheat temperature T_RH and d extra dimensions, the energy density stored in gravitons scales as T_RH^{2+3d/2}. Overclosure and Big Bang nucleosynthesis therefore stringently constrain T_RH in all universal extra dimension scenarios. At the same time, there is a window of reheat temperatures low enough to avoid these constraints and high enough to generate the desired thermal relic density for KK WIMP and superWIMP dark matter.

We revise existing limits on the D-dimensional Planck scale M_D from the nonobservation of microscopic black holes produced by high energy cosmic neutrinos in scenarios with D=4+n large extra dimensions. Previous studies have neglected the energy radiated in gravitational waves by the multipole moments of the incoming shock waves. We include the effects of energy loss, as well as form factors for black hole production and recent null results from cosmic ray detectors. For n>4, we obtain M_D > 1.0 - 1.4 TeV. These bounds are among the most stringent and conservative to date.

We investigate a new class of dark matter: superweakly-interacting massive particles (superWIMPs). As with conventional WIMPs, superWIMPs appear in well-motivated particle theories with naturally the correct relic density. In contrast to WIMPs, however, superWIMPs are impossible to detect in all conventional dark matter searches. We consider the concrete examples of gravitino and graviton cold dark matter in models with supersymmetry and universal extra dimensions, respectively, and show that superWIMP dark matter satisfies stringent constraints from Big Bang nucleosynthesis and the cosmic microwave background.

In previous works, a generic dynamical model has been suggested by Huang et al., which is shown to be correct for both adiabatic and radiative blastwaves, and in both ultra-relativistic and non-relativistic phases. In deriving their equations, Huang et al. have assumed that the radiative efficiency of the fireball is constant. They then applied their model directly to realistic cases where the radiative efficiency evolves with time. In this paper, we abandon the above assumption and re-derive a more accurate dynamical equation for gamma-ray burst remnants. Numerical results show that Huang et al.'s model is accurate enough in general cases.

Ultra-high energy cosmic neutrinos are incisive probes of both astrophysical sources and new TeV-scale physics. Such neutrinos would create extensive air showers deep in the atmosphere. The absence of such showers implies upper limits on incoming neutrino fluxes and cross sections. Combining the exposures of AGASA, the largest existing ground array, with the exposure of the Fly's Eye fluorescence detector integrated over all its operating epochs, we derive 95% CL bounds that substantially improve existing limits. We begin with model-independent bounds on astrophysical fluxes, assuming standard model cross sections, and model-independent bounds on new physics cross sections, assuming a conservative cosmogenic flux. We then derive model-dependent constraints on new components of neutrino flux for several assumed power spectra, and we update bounds on the fundamental Planck scale M_D in extra dimension scenarios from black hole production. For large numbers of extra dimensions, we find M_D > 2.0 (1.1) TeV for \mbhmin = M_D (5 M_D), comparable to or exceeding the most stringent constraints to date.

We propose that cold dark matter is made of Kaluza-Klein particles and explore avenues for its detection. The lightest Kaluza-Klein state is an excellent dark matter candidate if standard model particles propagate in extra dimensions and Kaluza-Klein parity is conserved. We consider Kaluza-Klein gauge bosons. In sharp contrast to the case of supersymmetric dark matter, these annihilate to hard positrons, neutrinos and photons with unsuppressed rates. Direct detection signals are also promising. These conclusions are generic to bosonic dark matter candidates.

If extra spacetime dimensions and low-scale gravity exist, black holes will be produced in observable collisions of elementary particles. For the next several years, ultra-high energy cosmic rays provide the most promising window on this phenomenon. In particular, cosmic neutrinos can produce black holes deep in the Earth's atmosphere, leading to quasi-horizontal giant air showers. We determine the sensitivity of cosmic ray detectors to black hole production and compare the results to other probes of extra dimensions. With n \ge 4 extra dimensions, current bounds on deeply penetrating showers from AGASA already provide the most stringent bound on low-scale gravity, requiring a fundamental Planck scale M_D > 1.3 - 1.8 TeV. The Auger Observatory will probe M_D as large as 4 TeV and may observe on the order of a hundred black holes in 5 years. We also consider the implications of angular momentum and possible exponentially suppressed parton cross sections; including these effects, large black hole rates are still possible. Finally, we demonstrate that even if only a few black hole events are observed, a standard model interpretation may be excluded by comparison with Earth-skimming neutrino rates.

Neutrinos with energies above 10^8 GeV are expected from cosmic ray interactions with the microwave background and are predicted in many speculative models. Such energetic neutrinos are difficult to detect, as they are shadowed by the Earth, but rarely interact in the atmosphere. Here we propose a novel detection strategy: Earth-skimming neutrinos convert to charged leptons that escape the Earth, and these leptons are detected in ground level fluorescence detectors. With the existing HiRes detector, neutrinos from some proposed sources are marginally detectable, and improvements of two orders of magnitude are possible at the proposed Telescope Array.

In spacetimes with asymmetric extra dimensions, cosmic neutrino interactions may be extraordinarily enhanced by p-brane production. Brane formation and decay may then initiate showers deep in the Earth's atmosphere at rates far above the standard model rate. We explore the p-brane discovery potential of cosmic ray experiments. The absence of deeply penetrating showers at AGASA already provides multi-TeV bounds on the fundamental Planck scale that significantly exceed those obtained from black hole production in symmetric compactification scenarios. This sensitivity will be further enhanced at the Auger Observatory. We also examine the possibility that p-brane formation resolves the GZK paradox. For flat compactifications, astrophysical bounds exclude this explanation. For warped scenarios, a solution could be consistent with the absence of deep showers only for extra dimensions with fine-tuned sizes well below the fundamental Planck length. In addition, it requires moderately penetrating showers, so far not reported, and ~100% modifications to standard model phenomenology at 100 GeV energies.

If spacetime has more than four dimensions, ultra-high energy cosmic rays may create microscopic black holes. Black holes created by cosmic neutrinos in the Earth will evaporate, and the resulting hadronic showers, muons, and taus may be detected in neutrino telescopes below the Earth's surface. We simulate such events in detail and consider black hole cross sections with and without an exponential suppression factor. We find observable rates in both cases: for conservative cosmogenic neutrino fluxes, several black hole events per year are observable at the IceCube detector; for fluxes at the Waxman-Bahcall bound, tens of events per year are possible. We also present zenith angle and energy distributions for all three channels. The ability of neutrino telescopes to differentiate hadrons, muons, and possibly taus, and to measure these distributions provides a unique opportunity to identify black holes, to experimentally constrain the form of black hole production cross sections, and to study Hawking evaporation.

Ultra-high energy cosmic rays create black holes in scenarios with extra dimensions and TeV-scale gravity. In particular, cosmic neutrinos will produce black holes deep in the atmosphere, initiating quasi-horizontal showers far above the standard model rate. At the Auger Observatory, hundreds of black hole events may be observed, providing evidence for extra dimensions and the first opportunity for experimental study of microscopic black holes. If no black holes are found, the fundamental Planck scale must be above 2 TeV for any number of extra dimensions.

Supersymmetry may be discovered at high energy colliders, through low energy precision measurements, and by dark matter searches. We present a comprehensive analysis of all available probes in minimal supergravity. This work extends previous analyses by including the focus point branch of parameter space and the full array of promising indirect dark matter searches. We find that particle and astrophysical searches underway are highly complementary: each separately provides only partial coverage of the available parameter space, but together they probe almost all models. Cosmology does {\em not} provide upper bounds on superpartner masses useful for future colliders. At the same time, in the cosmologically preferred region, if supersymmetry is to be observable at a 500 GeV linear collider, some signature of supersymmetry must appear {\em before} the LHC.

We review the prospects for discovering supersymmetric dark matter in a recently proposed set of post-LEP supersymmetric benchmark scenarios. We consider direct detection through spin-independent nuclear scattering, as well as indirect detection through relic annihilations to neutrinos, photons, and positrons. We find that several of the benchmark scenarios offer good prospects for direct detection through spin-independent nuclear scattering, as well as indirect detection through muons produced by neutrinos from relic annihilations in the Sun, and photons from annihilations in the galactic center.

A new set of supersymmetric benchmark scenarios has recently been proposed in the context of the constrained MSSM (CMSSM) with universal soft supersymmetry-breaking masses, taking into account the constraints from LEP, $b \to s \gamma$ and $g_\mu - 2$. These points have previously been used to discuss the physics reaches of different accelerators. In this paper, we discuss the prospects for discovering supersymmetric dark matter in these scenarios. We consider direct detection through spin-independent and spin-dependent nuclear scattering, as well as indirect detection through relic annihilations to neutrinos, photons, and positrons. We find that several of the benchmark scenarios offer good prospects for direct detection via spin-independent nuclear scattering and indirect detection via muons produced by neutrinos from relic annihilations inside the Sun, and some models offer good prospects for detecting photons from relic annihilations in the galactic centre.

Dark matter candidates arising in models of particle physics incorporating weak scale supersymmetry may produce detectable signals through their annihilation into neutrinos, photons, or positrons. A large number of relevant experiments are planned or underway. The `logically possible' parameter space is unwieldy. By working in the framework of minimal supergravity, we can survey the implications of the experiments for each other, as well as for direct searches, collider searches, low-energy experiments, and naturalness in a transparent fashion. We find that a wide variety of experiments provide interesting probes. Particularly promising signals arise in the mixed gaugino-Higgsino region. This region is favored by low-energy particle physics constraints and arises naturally from minimal supergravity due to the focus point mechanism. Indirect dark matter searches and traditional particle searches are highly complementary. In cosmologically preferred models, if there are charged superpartners with masses below 250 GeV, then some signature of supersymmetry must appear before the LHC begins operation.

In recent work, it has been argued that multi-TeV masses for scalar superpartners are not unnatural. Indeed, they appear to have significant phenomenological virtues. Here we explore the implications of such `focus point' supersymmetry for the dark matter problem. We find that constraints on relic densities do not place upper bounds on neutralino or scalar masses. We demonstrate that, in the specific context of minimal supergravity, a cosmologically stable mixed gaugino-Higgsino state emerges as an excellent, robust dark matter candidate. We estimate that, over a wide range of the unknown parameters, the spin-independent proton-neutralino cross sections fall in the range accessible to planned search experiments.

We study the physics of spontaneously broken family symmetries acting on the third generation. Massless familons (or Majorons) $f$ associated with such broken symmetries are motivated especially by cosmological scenarios with decaying tau neutrinos. We first note that, in marked contrast with the case for the first two generations, constraints on third generation familon couplings are poor, and are, in fact, non-existent at present in the hadronic sector. We derive new bounds from $B^0$--$\bar{B}^0$ mixing, $B^0 \to l^+ l'^-$, $b\to s\nu\bar{\nu}$, and astrophysics. The resulting constraints on familon decay constants are still much weaker than those for the first and second generation. We then discuss the promising prospects for significant improvements from searches for $\tau\to l f$, $B\to (\pi, K) f$, and $b\to (d,s) f$ with the current CLEO, ARGUS, and LEP data. Finally, we note that future constraints from CLEO III and the $B$ factories will probe decay constants beyond 10^8 GeV, well within regions of parameter space favored by proposed scenarios in neutrino cosmology.