Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.

We study tensor modes in pure natural inflation (arXiv:1706.08522), a recently-proposed inflationary model in which an axionic inflaton couples to pure Yang-Mills gauge fields. We find that the tensor-to-scalar ratio r is naturally bounded from below. This bound originates from the finiteness of the number of metastable branches of vacua in pure Yang-Mills theories. Details of the model can be probed by future cosmic microwave background experiments and improved lattice gauge theory calculations of the theta-angle dependence of the vacuum energy.

We point out that a simple inflationary model in which the axionic inflaton couples to a pure Yang-Mills theory may give the scalar spectral index (n_s) and tensor-to-scalar ratio (r) in complete agreement with the current observational data.

We present a simple and natural dark sector model in which dark matter particles arise as composite states of hidden strong dynamics and their stability is ensured by accidental symmetries. The model has only a few free parameters. In particular, the gauge symmetry of the model forbids the masses of dark quarks, and the confinement scale of the dynamics provides the unique mass scale of the model. The gauge group contains an Abelian symmetry $U(1)_D$, which couples the dark and standard model sectors through kinetic mixing. This model, despite its simple structure, has rich and distinctive phenomenology. In the case where the dark pion becomes massive due to $U(1)_D$ quantum corrections, direct and indirect detection experiments can probe thermal relic dark matter which is generically a mixture of the dark pion and the dark baryon, and the Large Hadron Collider can discover the $U(1)_D$ gauge boson. Alternatively, if the dark pion stays light due to a specific $U(1)_D$ charge assignment of the dark quarks, then the dark pion constitutes dark radiation. The signal of this radiation is highly correlated with that of dark baryons in dark matter direct detection.

An interesting possibility for dark matter is a scalar particle of mass of order 10 MeV-1 GeV, interacting with a U(1) gauge boson (dark photon) which mixes with the photon. We present a simple and natural model realizing this possibility. The dark matter arises as a composite pseudo Nambu-Goldstone boson (dark pion) in a non-Abelian gauge sector, which also gives a mass to the dark photon. For a fixed non-Abelian gauge group, SU(N), and a U(1) charge of the constituent dark quarks, the model has only three free parameters: the dynamical scale of the non-Abelian gauge theory, the gauge coupling of the dark photon, and the mixing parameter between the dark and standard model photons. In particular, the gauge symmetry of the model does not allow any mass term for the dark quarks, and stability of the dark pion is understood as a result of an accidental global symmetry. The model has a significant parameter space in which thermal relic dark pions comprise all of the dark matter, consistently with all experimental and cosmological constraints. In a corner of the parameter space, the discrepancy of the muon g-2 between experiments and the standard model prediction can also be ameliorated due to a loop contribution of the dark photon. Smoking-gun signatures of the model include a monophoton signal from the e+ e- collision into a photon and a "dark rho meson." Observation of two processes in e+ e- collision, the mode into the dark photon and that into the dark rho meson, would provide strong evidence for the model.

We propose a simple mechanism to suppress axion isocurvature fluctuations using hidden sector magnetic monopoles. This allows for the Peccei-Quinn scale to be of order the unification scale consistently with high scale inflation.

A class of supersymmetric grand unified theories is introduced that has a single scale below the cutoff, that of the supersymmetry breaking masses $\tilde{m}$. For a wide range of the dimensionless parameters, agreement with the observed mass of the Higgs boson determines $\tilde{m} \sim 10^9-10^{13} {\rm GeV}$, yielding Intermediate Scale Supersymmetry. We show that within this framework it is possible for seesaw neutrino masses, axions, and inflation to be described by the scale $\tilde{m}$, offering the possibility of a unified origin of disparate phenomena. Neutrino masses allowing for thermal leptogenesis can be obtained, and the axion decay constant lies naturally in the range $f_a \sim 10^9-10^{11} {\rm GeV}$, consistent with a recent observational suggestion of high scale inflation. A minimal $SU(5)$ model is presented that illustrates these features. In this model, the only states at the grand unified scale are those of the heavy gauge supermultiplet. The grand unified partners of the Higgs doublets have a mass of order $\tilde{m}$, leading to the dominant proton decay mode $p \rightarrow \bar{\nu} K^+$, which may be probed in upcoming experiments. Dark matter may be winos, with mass environmentally selected to the TeV scale, and/or axions. Gauge coupling unification is found to be successful, especially if the wino is at the TeV scale.

Models of cosmic inflation posit an early phase of accelerated expansion of the universe, driven by the dynamics of one or more scalar fields in curved spacetime. Though detailed assumptions about fields and couplings vary across models, inflation makes specific, quantitative predictions for several observable quantities, such as the flatness parameter ($\Omega_k = 1 - \Omega$) and the spectral tilt of primordial curvature perturbations ($n_s - 1 = d \ln {\cal P}_{\cal R} / d \ln k$), among others---predictions that match the latest observations from the {\it Planck} satellite to very good precision. In the light of data from {\it Planck} as well as recent theoretical developments in the study of eternal inflation and the multiverse, we address recent criticisms of inflation by Ijjas, Steinhardt, and Loeb. We argue that their conclusions rest on several problematic assumptions, and we conclude that cosmic inflation is on a stronger footing than ever before.

The discovery of accelerating expansion of the universe has led us to take the dramatic view that our universe may be one of the many universes in which low energy physical laws take different forms: the multiverse. I explain why/how this view is supported both observationally and theoretically, especially by string theory and eternal inflation. I then describe how quantum mechanics plays a crucial role in understanding the multiverse, even at the largest distance scales. The resulting picture leads to a revolutionary change of our view of spacetime and gravity, and completely unifies the paradigm of the eternally inflating multiverse with the many worlds interpretation of quantum mechanics. The picture also provides a solution to a long-standing problem in eternal inflation, called the measure problem, which I briefly describe.

The eternally inflating multiverse provides a consistent framework to understand coincidences and fine-tuning in the universe. As such, it provides the possibility of finding another coincidence: if the amount of slow-roll inflation was only slightly more than the anthropic threshold, then spatial curvature might be measurable. We study this issue in detail, particularly focusing on the question: "If future observations reveal nonzero curvature, what can we conclude?" We find that whether an observable signal arises or not depends crucially on three issues: the cosmic history just before the observable inflation, the measure adopted to define probabilities, and the nature of the correlation between the tunneling and slow-roll parts of the potential. We find that if future measurements find positive curvature at \Omega_k < -10^-4, then the framework of the eternally inflating multiverse is excluded with high significance. If the measurements instead reveal negative curvature at \Omega_k > 10^-4, then we can conclude (1) diffusive (new or chaotic) eternal inflation did not occur in our immediate past; (2) our universe was born by a bubble nucleation; (3) the probability measure does not reward volume increase; and (4) the origin of the observed slow-roll inflation is an accidental feature of the potential, not due to a theoretical mechanism. Discovery of \Omega_k > 10^-4 would also give us nontrivial information about the correlation between tunneling and slow-roll; e.g. a strong correlation favoring large N would be excluded in certain measures. We also ask whether the current constraint on \Omega_k is consistent with multiverse expectations, finding that the answer is yes, except for certain cases. In the course of this work we were led to consider vacuum decay branching ratios, and found that it is more likely than one might guess that the decays are dominated by a single channel.

We consider the multiverse in the intrinsically quantum mechanical framework recently proposed in Refs. [1,2]. By requiring that the principles of quantum mechanics are universally valid and that physical predictions do not depend on the reference frame one chooses to describe the multiverse, we find that the multiverse state must be static---in particular, the multiverse does not have a beginning or end. We argue that, despite its naive appearance, this does not contradict observation, including the fact that we observe that time flows in a definite direction. Selecting the multiverse state is ultimately boiled down to finding normalizable solutions to certain zero-eigenvalue equations, analogous to the case of the hydrogen atom. Unambiguous physical predictions would then follow, according to the rules of quantum mechanics.

Quantum mechanics introduces the concept of probability at the fundamental level, yielding the measurement problem. On the other hand, recent progress in cosmology has led to the "multiverse" picture, in which our observed universe is only one of the many, bringing an apparent arbitrariness in defining probabilities, called the measure problem. In this paper, we discuss how these two problems are intimately related with each other, developing a complete picture for quantum measurement and cosmological histories in the quantum mechanical universe. On one hand, quantum mechanics eliminates the arbitrariness of defining probabilities in the multiverse, as discussed in arXiv:1104.2324. On the other hand, the multiverse allows for understanding why we observe an ordered world obeying consistent laws of physics, by providing an infinite-dimensional Hilbert space. This results in the irreversibility of quantum measurement, despite the fact that the evolution of the multiverse state is unitary. In order to describe the cosmological dynamics correctly, we need to identify the structure of the Hilbert space for a system with gravity. We argue that in order to keep spacetime locality, the Hilbert space for dynamical spacetime must be defined only in restricted spacetime regions: in and on the (stretched) apparent horizon as viewed from a fixed reference frame. This requirement arises from eliminating all the redundancies and overcountings in a general relativistic, global spacetime description of nature. It is responsible for horizon complementarity as well as the "observer dependence" of horizons/spacetime---these phenomena arise to represent changes of the reference frame in the relevant Hilbert space. This can be viewed as an extension of the Poincare transformation in the quantum gravitational context, as the Lorentz transformation is viewed as an extension of the Galilean transformation.

Recently, a new framework for describing the multiverse has been proposed which is based on the principles of quantum mechanics. The framework allows for well-defined predictions, both regarding global properties of the universe and outcomes of particular experiments, according to a single probability formula. This provides complete unification of the eternally inflating multiverse and many worlds in quantum mechanics. In this paper we elucidate how cosmological parameters can be calculated in this framework, and study the probability distribution for the value of the cosmological constant. We consider both positive and negative values, and find that the observed value is consistent with the calculated distribution at an order of magnitude level. In particular, in contrast to the case of earlier measure proposals, our framework prefers a positive cosmological constant over a negative one. These results depend only moderately on how we model galaxy formation and life evolution therein.

We present a framework in which well-defined predictions are obtained in an eternally inflating multiverse, based on the principles of quantum mechanics. We show that the entire multiverse is described purely from the viewpoint of a single "observer," who describes the world as a quantum state defined on his/her past light cones bounded by the (stretched) apparent horizons. We find that quantum mechanics plays an essential role in regulating infinities. The framework is "gauge invariant," i.e. predictions do not depend on how spacetime is parametrized, as it should be in a theory of quantum gravity. Our framework provides a fully unified treatment of quantum measurement processes and the multiverse. We conclude that the eternally inflating multiverse and many worlds in quantum mechanics are the same. Other important implications include: global spacetime can be viewed as a derived concept; the multiverse is a transient phenomenon during the world relaxing into a supersymmetric Minkowski state. We also present a theory of "initial conditions" for the multiverse. By extrapolating our framework to the extreme, we arrive at a picture that the entire multiverse is a fluctuation in the stationary, fractal "mega-multiverse," in which an infinite sequence of multiverse productions occurs. The framework discussed here does not suffer from problems/paradoxes plaguing other measures proposed earlier, such as the youngness paradox, the Boltzmann brain problem, and a peculiar "end" of time.

We explore the possibility that both the weak scale and the thermal relic dark matter abundance are environmentally selected in a multiverse. An underlying supersymmetric theory containing the states of the MSSM and singlets, with supersymmetry and R symmetry broken at unified scales, has just two realistic low energy effective theories. One theory, (SM + \tilde{w}), is the Standard Model augmented only by the wino, having a mass near 3 TeV, and has a Higgs boson mass in the range of (127 - 142) GeV. The other theory, (SM + \tilde{h}/\tilde{s}), has Higgsinos and a singlino added to the Standard Model. The Higgs boson mass depends on the single new Yukawa coupling of the theory, y, and is near 141 GeV for small y but grows to be as large as 210 GeV as this new coupling approaches strong coupling at high energies. Much of the parameter space of this theory will be probed by direct detection searches for dark matter that push two orders of magnitude below the present bounds; furthermore, the dark matter mass and cross section on nucleons are correlated with the Higgs boson mass. The indirect detection signal of monochromatic photons from the galactic center is computed, and the range of parameters that may be accessible to LHC searches for trilepton events is explored. Taking a broader view, allowing the possibility of R symmetry protection to the TeV scale or axion dark matter, we find four more theories: (SM + axion), two versions of Split Supersymmetry, and the E-MSSM, where a little supersymmetric hierarchy is predicted. The special Higgs mass value of (141 \pm 2) GeV appears in symmetry limits of three of the six theories, (SM + axion), (SM + \tilde{w}) and (SM + \tilde{h}/\tilde{s}), motivating a comparison of other signals of these three theories.

There is a deep cosmological mystery: although dependent on very different underlying physics, the timescales of structure formation, of galaxy cooling (both radiatively and against the CMB), and of vacuum domination do not differ by many orders of magnitude, but are all comparable to the present age of the universe. By scanning four landscape parameters simultaneously, we show that this quadruple coincidence is resolved. We assume only that the statistical distribution of parameter values in the multiverse grows towards certain catastrophic boundaries we identify, across which there are drastic regime changes. We find order-of-magnitude predictions for the cosmological constant, the primordial density contrast, the temperature at matter-radiation equality, the typical galaxy mass, and the age of the universe, in terms of the fine structure constant and the electron, proton and Planck masses. Our approach permits a systematic evaluation of measure proposals; with the causal patch measure, we find no runaway of the primordial density contrast and the cosmological constant to large values.

Cosmologically long-lived, composite states arise as natural dark matter candidates in theories with a strongly interacting hidden sector at a scale of 10 - 100 TeV. Light axion-like states, with masses in the 1 MeV - 10 GeV range, are also generic, and can decay via Higgs couplings to light standard model particles. Such a scenario is well motivated in the context of very low energy supersymmetry breaking, where ubiquitous cosmological problems associated with the gravitino are avoided. We investigate the astrophysical and collider signatures of this scenario, assuming that dark matter decays into the axion-like states via dimension six operators, and we present an illustrative model exhibiting these features. We conclude that the recent data from PAMELA, FERMI, and H.E.S.S. points to this setup as a compelling paradigm for dark matter. This has important implications for future diffuse gamma ray measurements and collider searches.

A leading interpretation of the electron/positron excesses seen by PAMELA and ATIC is dark matter annihilation in the galactic halo. Depending on the annihilation channel, the electron/positron signal could be accompanied by a galactic gamma ray or neutrino flux, and the non-detection of such fluxes constrains the couplings and halo properties of dark matter. In this paper, we study the interplay of electron data with gamma ray and neutrino constraints in the context of cascade annihilation models, where dark matter annihilates into light degrees of freedom which in turn decay into leptons in one or more steps. Electron and muon cascades give a reasonable fit to the PAMELA and ATIC data. Compared to direct annihilation, cascade annihilations can soften gamma ray constraints from final state radiation by an order of magnitude. However, if dark matter annihilates primarily into muons, the neutrino constraints are robust regardless of the number of cascade decay steps. We also examine the electron data and gamma ray/neutrino constraints on the recently proposed "axion portal" scenario.

Motivated by the galactic positron excess seen by PAMELA and ATIC/PPB-BETS, we propose that dark matter is a TeV-scale particle that annihilates into a pseudoscalar "axion." The positron excess and the absence of an anti-proton or gamma ray excess constrain the axion mass and branching ratios. In the simplest realization, the axion is associated with a Peccei-Quinn symmetry, in which case it has a mass around 360-800 MeV and decays into muons. We present a simple and predictive supersymmetric model implementing this scenario, where both the Higgsino and dark matter obtain masses from the same source of TeV-scale spontaneous symmetry breaking.

The NSF has chosen the site for the Deep Underground Science and Engineering Laboratory (DUSEL) to be in Lead, South Dakota. In fact, the state of South Dakota has already stepped up to the plate and contributed its own funding for the proposed lab, see http://www.sanfordlaboratoryathomestake.org/index.html. The final decision by NSF for funding the Initial Suite of Experiments for DUSEL will be made early in 2009. At that time the NSF Science Board must make a decision. Of order 200 experimentalists have already expressed an interest in performing experiments at DUSEL. In order to assess the interest of the theoretical community, the Center for Cosmology and Astro-Particle Physics (CCAPP) at The Ohio State University (OSU) organized a 3-day DUSEL Theory Workshop in Columbus, Ohio from April 4 - 6, 2008. The workshop focused on the scientific case for six proposed experiments for DUSEL: long baseline neutrino oscillations, proton decay, dark matter, astrophysical neutrinos, neutrinoless double beta decay and N-Nbar oscillations. The outcome of this workshop is the DUSEL Theory White paper addressing the scientific case at a level which may be useful in the decision making process for policy makers at the NSF and in the U.S. Congress. In order to assess the physics interest in the DUSEL project we have posted the DUSEL Theory White paper on the following CCAPP link http://ccapp.osu.edu/whitepaper.html . Please read the white paper and, if you are interested, use the link to show your support by co-signing the white paper.

In any theory it is unnatural if the observed parameters lie very close to special values that determine the existence of complex structures necessary for observers. A naturalness probability, P, is introduced to numerically evaluate the unnaturalness. If P is small in all known theories, there is an observer naturalness problem. In addition to the well-known case of the cosmological constant, we argue that nuclear stability and electroweak symmetry breaking (EWSB) represent significant observer naturalness problems. The naturalness probability associated with nuclear stability is conservatively estimated as P_nuc < 10^{-(3-2)}, and for simple EWSB theories P_EWSB < 10^{-(2-1)}. This pattern of unnaturalness in three different arenas, cosmology, nuclear physics, and EWSB, provides evidence for the multiverse. In the nuclear case the problem is largely solved even with a flat multiverse distribution, and with nontrivial distributions it is possible to understand both the proximity to neutron stability and the values of m_e and m_d - m_u in terms of the electromagnetic contribution to the proton mass. It is reasonable that multiverse distributions are strong functions of Lagrangian parameters due to their dependence on various factors. In any EWSB theory, strongly varying distributions typically lead to a little or large hierarchy, and in certain multiverses the size of the little hierarchy is enhanced by a loop factor. Since the correct theory of EWSB is unknown, our estimate for P_EWSB is theoretical. The LHC will determine P_EWSB more robustly, which may remove or strengthen the observer naturalness problem of EWSB. For each of the three arenas, the discovery of a natural theory would eliminate the evidence for the multiverse; but in the absence of such a theory, the multiverse provides a provisional understanding of the data.

We show that the solid lower bound of about 10^{-44} cm^2 is obtained for the cross section between the supersymmetric dark matter and nucleon in a theory in which the supersymmetric fine-tuning problem is solved without extending the Higgs sector at the weak scale. This bound arises because of relatively small superparticle masses and a fortunate correlation that the two dominant diagrams for the dark matter detection always interfere constructively if the constraint from the b -> s \gamma measurements is obeyed. It is, therefore, quite promising in the present scenario that the supersymmetric dark matter is discovered before the LHC, assuming that the dark matter is the lightest supersymmetric particle.

Theories of evolving quintessence are constructed that generically lead to deviations from the w = -1 prediction of non-evolving dark energy. The small mass scale that governs evolution, m_\phi \approx 10^{-33} eV, is radiatively stable, and the ``Why Now?'' problem is solved. These results rest crucially on seesaw cosmology: in broad outline, fundamental physics and cosmology can be understood from only two mass scales, the weak scale, v, and the Planck scale, M. Requiring a scale of dark energy \rho_{DE}^{1/4} governed by v^2/M, and a radiatively stable evolution rate m_\phi given by v^4/M^3, leads to a distinctive form for the equation of state w(z) that follows from a cosine quintessence potential. An explicit hidden axion model is constructed. Dark energy resides in the potential of the axion field which is generated by a new QCD-like force that gets strong at the scale \Lambda \approx v^2/M \approx \rho_{DE}^{1/4}. The evolution rate is given by a second seesaw that leads to the axion mass, m_\phi \approx \Lambda^2/f, with f \approx M.

Simple models are constructed for "acceleressence" dark energy: the latent heat of a phase transition occurring in a hidden sector governed by the seesaw mass scale v^2/M_Pl, where v is the electroweak scale and M_Pl the gravitational mass scale. In our models, the seesaw scale is stabilized by supersymmetry, implying that the LHC must discover superpartners with a spectrum that reflects a low scale of fundamental supersymmetry breaking. Newtonian gravity may be modified by effects arising from the exchange of fields in the acceleressence sector whose Compton wavelengths are typically of order the millimeter scale. There are two classes of models. In the first class the universe is presently in a metastable vacuum and will continue to inflate until tunneling processes eventually induce a first order transition. In the simplest such model, the range of the new force is bounded to be larger than 25 microns in the absence of fine-tuning of parameters, and for couplings of order unity it is expected to be \approx 100 microns. In the second class of models thermal effects maintain the present vacuum energy of the universe, but on further cooling, the universe will "soon" smoothly relax to a matter dominated era. In this case, the range of the new force is also expected to be of order the millimeter scale or larger, although its strength is uncertain. A firm prediction of this class of models is the existence of additional energy density in radiation at the eV era, which can potentially be probed in precision measurements of the cosmic microwave background. An interesting possibility is that the transition towards a matter dominated era has occurred in the very recent past, with the consequence that the universe is currently decelerating.

We study gauge mediation models of supersymmetry breaking with neutralino LSP. These models are naturally realized by embedding the usual four-dimensional gauge mediation models into a higher-dimensional spacetime such as $M^4 \times S^1 / Z_{2}$. We calculate the relic abundance of the neutralino LSP in these models and show that there exist wide parameter regions where the neutralino LSP constitutes the dominant component of the cold dark matter. These regions evade constraints from collider experiments such as Higgs mass bounds and $b \ra s \gamma$, and also provide the value for the muon anomalous magnetic moment which is consistent with the SUSY explanation of the deviation from the standard model prediction.

Recent cosmological observations suggest the presence of small but nonzero cosmological constant $\Lambda_{cos}^4 \simeq (2 \times 10^{-3} eV)^4$. It is an intriguing possibility that such a small cosmological constant is supplied by the potential energy density of an ultralight axion-like field (called as quintessence axion). If this axion couples to the electroweak gauge fields, its potential may be generated by the electroweak instanton effects. We calculate the axion potential assuming the supersymmetric standard model and obtain a surprising result that the induced energy density of the quintessence axion field is very close to the value suggested from the observations.

If a global chiral symmetry is explicitly broken by anomalies in nonabelian gauge theories, a pseudo Nambu-Goldstone boson (axion) associated with a spontaneous breakdown of such a global symmetry acquires a mass through nonperturbative instanton effects. We calculate the axion mass assuming a supersymmetric SU(2) gauge theory and show that the axion obtains an extremely small mass when the SU(2) gauge symmetry is broken down at very high energy, say at the Planck scale. We identify the axion with a hypothetical ultralight boson field proposed to account for a small but nonzero cosmological constant suggested from recent cosmological observations.

Superheavy particles of masses $\simeq 10^{13}-10^{14} GeV$ with lifetimes $\simeq 10^{10}-10^{22} years$ are very interesting, since their decays may account for the ultra-high energy (UHE) cosmic rays discovered beyond the Greisen-Zatsepin-Kuzmin cut-off energy $E \sim 5 \times 10^{10} GeV$. We show that the presence of such long-lived superheavy particles is a generic prediction of QCD-like SU(N_c) gauge theories with N_f flavors of quarks and antiquarks and the large number of colors N_c. We construct explicit models based on supersymmetric SU(N_c) gauge theories and show that if the dynamical scale $\Lambda \simeq 10^{13}-10^{14} GeV$ and N_c = 6-10 the lightest composite baryons have the desired masses and lifetimes to explain the UHE cosmic rays. Interesting is that in these models the gaugino condensation necessarily occurs and hence these models may play a role of so-called hidden sector for supersymmetry breaking in supergravity.

The recently observed ultra-high energy (UHE) cosmic rays beyond the Greisen-Zatsepin-Kuzmin bound can be explained by the decays of some superheavy $X$ particles forming a part of dark matter in our universe. We consider various discrete gauge symmetries ${\bf Z}_N$ to ensure the required long lifetime ($\tau_X \simeq 10^{10}-10^{22} years$) of the $X$ particle to explain the UHE cosmic rays in the minimal supersymmetric standard model (MSSM) with massive Majorana neutrinos. We show that there is no anomaly-free discrete gauge symmetry to make the lifetime of the $X$ particle sufficiently long in the MSSM with the $X$ particle. We find, however, possible solutions to this problem especially by enlarging the particle contents in the MSSM. We show a number of solutions introducing an extra pair of singlets $Y$ and $\bar{Y}$ which have fractional ${\bf Z}_N$ (N=2,3) charges. The present experimental constraints on the $X$ particle are briefly discussed.

We show that there are discrete gauge symmetries protect naturally heavy X particles from decaying into the ordinary light particles in the supersymmetric standard model. This makes the proposal very attractive that the superheavy X particles constitute a part of the dark matter in the present universe. It is more interesting that there are a class of discrete gauge symmetries which naturally accommodate a long-lived unstable X particle. We find that in some discrete Z_{10} models, for example, a superheavy X particle has lifetime \tau_X \simeq 10^{11}-10^{26} years for its mass M_X \simeq 10^{13}-10^{14} GeV. This long lifetime is guaranteed by the absence of lower dimensional operators (of light particles) couple to the X. We briefly discuss a possible explanation for the recently observed ultra-high-energy cosmic ray events by the decay of this unstable X particle.