The Late Devonian was a protracted period of low speciation resulting in biodiversity decline, culminating in extinction events near the Devonian-Carboniferous boundary. Recent evidence indicates that the final extinction event may have coincided with a dramatic drop in stratospheric ozone, possibly due to a global temperature rise. Here we study an alternative possible cause for the postulated ozone drop: a nearby supernova explosion that could inflict damage by accelerating cosmic rays that can deliver ionizing radiation for up to $\sim 100$ kyr. We therefore propose that end-Devonian extinction was triggered by one or more supernova explosions at $\sim 20 \ \rm pc$, somewhat beyond the ``kill distance'' that would have precipitated a full mass extinction. Nearby supernovae are likely due to core-collapses of massive stars in clusters in the thin Galactic disk in which the Sun resides. Detecting any of the long-lived radioisotopes \sm146, \u235 or \pu244 in one or more end-Devonian extinction strata would confirm a supernova origin, point to the core-collapse explosion of a massive star, and probe supernova nucleosythesis. Other possible tests of the supernova hypothesis are discussed.

Atom interferometers (AIs) on earth and in space offer good capabilities for measuring gravitational waves (GWs) in the mid-frequency deciHz band, complementing the sensitivities of the LIGO/Virgo and LISA experiments and enabling probes of possible modifications of the general relativity predictions for GW propagation. We illustrate these capabilities using the projected sensitivities of the AION (terrestrial) and AEDGE (space-based) AI projects, showing that AION could improve the present LIGO/Virgo direct limit on the graviton mass by a factor $\sim 40$ to $\simeq 10^{-24}\,$eV, and AEDGE could improve the limit by another order of magnitude. AION and AEDGE will also have greater sensitivity than LIGO to some scenarios for Lorentz violation.

We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al.; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision multi-band measurements in several pulsars indicate frequency-dependent DMs. The TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses.

We present time-of-arrival measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of time-of-arrival (TOA) measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for fourteen pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications.

We survey systematically the general parametrisations of particle-physics models for a first-order phase transition in the early universe, including models with polynomial potentials both with and without barriers at zero temperature, and Coleman-Weinberg-like models with potentials that are classically scale-invariant. We distinguish three possibilities for the transition - detonations, deflagrations and hybrids - and consider sound waves and turbulent mechanisms for generating gravitational waves during the transitions in these models, checking in each case the requirement for successful percolation. We argue that in models without a zero-temperature barrier and in scale-invariant models the period during which sound waves generate gravitational waves lasts only for a fraction of a Hubble time after a generic first-order cosmological phase transition, whereas it may last longer in some models with a zero-temperature barrier that feature severe supercooling. We illustrate the implications of these results for future gravitational-wave experiments.

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

We discuss the observability of circular polarisation of the stochastic gravitational-wave background (SGWB) generated by helical turbulence following a first-order cosmological phase transition, using a model that incorporates the effects of both direct and inverse energy cascades. We explore the strength of the gravitational-wave signal and the dependence of its polarisation on the helicity fraction, $\zeta_*$, the strength of the transition, $\alpha$, the bubble size, $R_*$, and the temperature, $T_*$, at which the transition finishes. We calculate the prospective signal-to-noise ratios of the SGWB strength and polarisation signals in the LISA experiment, exploring the parameter space in a way that is minimally sensitive to the underlying particle physics model. We find that discovery of SGWB polarisation is generally more challenging than measuring the total SGWB signal, but would be possible for appropriately strong transitions with large bubble sizes and a substantial polarisation fraction.

We outline the experimental concept and key scientific capabilities of AION (Atom Interferometer Observatory and Network), a proposed UK-based experimental programme using cold strontium atoms to search for ultra-light dark matter, to explore gravitational waves in the mid-frequency range between the peak sensitivities of the LISA and LIGO/Virgo/ KAGRA/INDIGO/Einstein Telescope/Cosmic Explorer experiments, and to probe other frontiers in fundamental physics. AION would complement other planned searches for dark matter, as well as probe mergers involving intermediate mass black holes and explore early universe cosmology. AION would share many technical features with the MAGIS experimental programme in the US, and synergies would flow from operating AION in a network with this experiment, as well as with other atom interferometer experiments such as MIGA, ZAIGA and ELGAR. Operating AION in a network with other gravitational wave detectors such as LIGO, Virgo and LISA would also offer many synergies.

Motivated by recent measurements of deposits of $^{60}$Fe on the ocean floor and the lunar surface, we model the transport of dust grains containing $^{60}$Fe from a near-Earth (i.e., within 100 pc) supernova (SN). We inject dust grains into the environment of an SN remnant (SNR) and trace their trajectories by applying a 1D hydrodynamic description assuming spherical symmetry to describe the plasma dynamics, and include a rudimentary, 3D magnetic field description to examine its influence on charged dust grains. We assume the interstellar medium (ISM) magnetic fields are turbulent, and are amplified by the SNR shock, while the SN wind and ejecta fields are negligible. We examine the various influences on the dust grains within the SNR to determine when/if the dust decouples from the plasma, how much it is sputtered, and where within the SNR the dust grains are located. We find that Rayleigh-Taylor instabilities are important for dust survival, as they influence the location of the SN's reverse shock. We find that the presence of a magnetic field within the shocked ISM material limits the passage of SN dust grains, with the field either reflecting or trapping the grains within the heart of the SNR. These results have important implications for {\it in situ} $^{60}$Fe measurements, and for dust evolution in SNRs generally.

We have recently proposed attractor models for modulus fixing, inflation, supersymmetry breaking and dark energy based on no-scale supergravity. In this paper we develop phenomenological and cosmological aspects of these no-scale attractor models that underpin their physical applications. We consider models in which inflation is driven by a modulus field ($T$-type) with supersymmetry broken by a Polonyi field, or a matter field ($\phi$-type) with supersymmetry broken by the modulus field. We derive the possible patterns of soft supersymmetry-breaking terms, which depend in $T$-type models whether the Polonyi and/or matter fields are twisted or not, and in $\phi$-type models on whether the inflaton and/or other matter fields are twisted or not. In $\phi$-type models, we are able to directly relate the scale of supersymmetry breaking to the inflaton mass. We also discuss cosmological constraints from entropy considerations and the density of dark matter on the mechanism for stabilizing the modulus field via higher-order terms in the no-scale K\"ahler potential.

We search for extrasolar planets around millisecond pulsars using pulsar timing data and seek to determine the minimum detectable planetary masses as a function of orbital period. Using the 11-year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), we look for variations from our models of pulse arrival times due to the presence of exoplanets. No planets are detected around the millisecond pulsars in the NANOGrav 11-year data set, but taking into consideration the noise levels of each pulsar and the sampling rate of our observations, we develop limits that show we are sensitive to planetary masses as low as that of the moon. We analyzed potential planet periods, P, in the range 7 days < P < 2000 days, with somewhat smaller ranges for some binary pulsars. The planetary mass limit for our median-sensitivity pulsar within this period range is 1 M_moon (P / 100 days)^(-2/3).

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $\eta M \sim 2\times10^{10}M_\odot$ and inclination of $\iota=\pi/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.

We present a minimal SU(2,1)/SU(2) x U(1) no-scale model that incorporates in an economical way modulus fixing, Starobinsky-like inflation, an adjustable scale for supersymmetry breaking and the possibility of a small cosmological constant, a.k.a. dark energy.

We develop a string-inspired model for particle cosmology, based on a flipped SU(5)$\times$U(1) gauge group formulated in a no-scale supergravity framework. The model realizes Starobinsky-like inflation, which we assume to be followed by strong reheating, with the GUT symmetry being broken subsequently by a light `flaton' field whose decay generates a second stage of reheating. We discuss the production of gravitinos and the non-thermal contribution made by their decays to the density of cold dark matter, which is assumed to be provided by the lightest neutralino. We also discuss the masses of light and heavy neutrinos and leptogenesis. As discussed previously [1], a key r\^ole is played by a superpotential coupling between the inflaton, matter and GUT Higgs fields, called $\lambda_6$. We scan over possible values of $\lambda_6$, exploring the correlations between the possible values of observables. We emphasize that the release of entropy during the GUT transition allows large regions of supersymmetry-breaking parameter space that would otherwise lead to severe overdensity of dark matter. Furthermore, we find that the Big Bang nucleosynthesis lower limit on the reheating temperature of $\sim 1$ MeV restricts the supersymmetry-breaking scale to a range ${\cal O}(10)$ TeV that is consistent with the absence of supersymmetric particles at the LHC.

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

We study the generation of intergalactic magnetic fields in two models for first-order phase transitions in the early Universe that have been studied previously in connection with the generation of gravitational waves (GWs): the Standard Model supplemented by an $|H|^6$ operator (SM+$H^6$) and a classically scale-invariant model with an extra gauged U(1) $B - L$ symmetry (SM$_{B-L}$). We consider contributions to magnetic fields generated by bubble collisions and by turbulence in the primordial plasma, and we consider the hypotheses that helicity is seeded in the gauge field or kinetically. We study the conditions under which the intergalactic magnetic fields generated may be larger than the lower bounds from blazar observations, and correlate them with the observability of GWs and possible collider signatures. In the SM+$H^6$ model bubble collisions alone cannot yield large enough magnetic fields, whereas turbulence may do so. In the SM$_{B-L}$ model bubble collisions and turbulence may both yield magnetic fields above the blazar bound unless the B$-$L gauge boson is very heavy. In both models there may be observable GW and collider signatures if sufficiently large magnetic fields are generated.

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.

Despite its importance to our understanding of physics at supranuclear densities, the equation of state (EoS) of matter deep within neutron stars remains poorly understood. Millisecond pulsars (MSPs) are among the most useful astrophysical objects in the Universe for testing fundamental physics, and place some of the most stringent constraints on this high-density EoS. Pulsar timing - the process of accounting for every rotation of a pulsar over long time periods - can precisely measure a wide variety of physical phenomena, including those that allow the measurement of the masses of the components of a pulsar binary system (Lorimer & Kramer 2005). One of these, called relativistic Shapiro delay (Shapiro 1964), can yield precise masses for both an MSP and its companion; however, it is only easily observed in a small subset of high-precision, highly inclined (nearly edge-on) binary pulsar systems. By combining data from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 12.5-year data set with recent orbital-phase-specific observations using the Green Bank Telescope, we have measured the mass of the MSP J0740+6620 to be $\mathbf{2.14^{+0.10}_{-0.09}}$ solar masses (68.3% credibility interval; 95.4% credibility interval is $\mathbf{2.14^{+0.20}_{-0.18}}$ solar masses). It is highly likely to be the most massive neutron star yet observed, and serves as a strong constraint on the neutron star interior EoS.

Measurements of the cosmic microwave background (CMB) favour models of inflation with a small tensor-to-scalar ratio r, as predicted by the Starobinsky R + R^2 model. It has been shown previously that various models based on no-scale supergravity with different forms of superpotential make predictions similar to those of the Starobinsky model. In this paper we present a unified and general treatment of Starobinsky avatars of no-scale supergravity, using the underlying non-compact SU(2,1)/SU(2) x U(1) symmetry to demonstrate equivalences between different models, exhibiting 6 specific equivalence classes.

We have presented previously a general treatment of Starobinsky-like inflation in no-scale supergravity where the tensor-to-scalar ratio $r = 3(1 - n_s)^2$, and $n_s$ is the tilt of the scalar perturbations. In particular, we have shown how this scenario can be unified with modulus fixing, supersymmetry breaking and a small cosmological constant. In this paper we extend these constructions to inflationary models based on generalized no-scale structures. In particular, we consider alternative values of the curvature parameter, $\alpha < 1$, as may occur if not all the complex K\"ahler moduli contribute to driving inflation, as well as $\alpha > 1$, as may occur if complex structure moduli also contribute to driving inflation. In all cases, we combine these $\alpha$-Starobinsky inflation models with supersymmetry breaking and a present-day cosmological constant, allowing for additional contributions to the vacuum energy from stages of gauge symmetry breaking.

We propose a complete cosmological scenario based on a flipped SU(5) $\times$ U(1) GUT model that incorporates Starobinsky-like inflation, taking the subsequent cosmological evolution carefully into account. A single master coupling, $\lambda_6$, connects the singlet, GUT Higgs and matter fields, controlling 1) inflaton decays and reheating, 2) the gravitino production rate and therefore the non-thermal abundance of the supersymmetric cold dark matter particle, 3) neutrino masses and 4) the baryon asymmetry of the Universe.

We derive efficiency factors for the production of gravitational waves through bubble collisions and plasma-related sources in strong phase transitions, and find the conditions under which the bubble collisions can contribute significantly to the signal. We use lattice simulations to clarify the dependence of the colliding bubbles on their initial state. We illustrate our findings in two examples, the Standard Model with an extra $|H|^6$ interaction and a classically scale-invariant $U(1)_{\rm B-L}$ extension of the Standard Model. The contribution to the GW spectrum from bubble collisions is found to be negligible in the $|H|^6$ model, whereas it can play an important role in parts of the parameter space in the scale-invariant $U(1)_{\rm B-L}$ model. In both cases the sound-wave period is much shorter than a Hubble time, suggesting a significant amplification of the turbulence-sourced signal. We find, however, that the peak of the plasma-sourced spectrum is still produced by sound waves with the slower-falling turbulence contribution becoming important off-peak.

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.

We report the results of a global analysis of dark matter simplified models (DMSMs) with leptophobic mediator particles of spin one, considering the cases of both vector and axial-vector interactions with dark matter (DM) particles and quarks. We require the DMSMs to provide all the cosmological DM density indicated by Planck and other observations, and we impose the upper limits on spin-independent and -dependent scattering from direct DM search experiments. We also impose all relevant LHC constraints from searches for monojet events and measurements of the dijet mass spectrum. We model the likelihood functions for all the constraints and combine them within the MasterCode framework, and probe the full DMSM parameter spaces by scanning over the mediator and DM masses and couplings, not fixing any of the model parameters. We find, in general, two allowed regions of the parameter spaces: one in which the mediator couplings to Standard Model (SM) and DM particles may be comparable to those in the SM and the cosmological DM density is reached via resonant annihilation, and one in which the mediator couplings to quarks are $\lesssim 10^{-3}$ and DM annihilation is non-resonant. We find that the DM and mediator masses may well lie within the ranges accessible to LHC experiments. We also present predictions for spin-independent and -dependent DM scattering, and present specific results for ranges of the DM couplings that may be favoured in ultraviolet completions of the DMSMs.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

Models of quantum gravity suggest that the vacuum should be regarded as a medium with quantum structure that may have non-trivial effects on photon propagation, including the violation of Lorentz invariance. Fermi Large Area Telescope (LAT) observations of gamma-ray bursts (GRBs) are sensitive probes of Lorentz invariance, via studies of energy-dependent timing shifts in their rapidly-varying photon emissions. In this paper we analyze the Fermi-LAT measurements of high-energy gamma rays from GRBs with known redshifts, allowing for the possibility of energy-dependent variations in emission times at the sources as well as a possible non-trivial refractive index in vacuo for photons. We use statistical estimators based on the irregularity, kurtosis and skewness of bursts that are relatively bright in the 100 MeV to multi-GeV energy band to constrain possible dispersion effects during propagation. We find that the energy scale characterizing a linear energy dependence of the refractive index should exceed a few $\times 10^{17}$ GeV, and we estimate the sensitivity attainable with additional future sources to be detected by Fermi-LAT.

What is the maximum possible strength of a first-order electroweak phase transition and the resulting gravitational wave (GW) signal? While naively one might expect that supercooling could increase the strength of the transition to very high values, for strong supercooling the Universe is no longer radiation-dominated and the vacuum energy of the unstable minimum of the potential dominates the expansion, which can jeopardize the successful completion of the phase transition. After providing a general treatment for the nucleation, growth and percolation of broken phase bubbles during a first-order phase transition that encompasses the case of significant supercooling, we study the conditions for successful bubble percolation and completion of the electroweak phase transition in theories beyond the Standard Model featuring polynominal potentials. For such theories, these conditions set a lower bound on the temperature of the transition. Since the plasma cannot be significantly diluted, the resulting GW signal originates mostly from sound waves and turbulence in the plasma, rather than bubble collisions. We find the peak frequency of the GW signal from the phase transition to be generically $f \gtrsim 10^{-4}$ Hz. We also study the condition for GW production by sound waves to be long-lasting (GW source active for approximately a Hubble time), showing it is generally not fulfilled in concrete scenarios. Because of this the sound wave GW signal could be weakened, with turbulence setting in earlier, resulting in a smaller overall GW signal as compared to current literature predictions.

There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact on the Earth, offering new probes of stellar evolution, nuclear astrophysics, the astrophysics of the solar neighborhood, cosmic-ray sources and acceleration, multi-messenger astronomy, and astrobiology. Interdisciplinary connections reach broadly to include heliophysics, geology, and evolutionary biology. Objectives for the future include pinning down the nature and location of the established near-Earth supernova explosions, seeking evidence for others, and searching for other short-lived isotopes such as 26Al and 244Pu. The unique information provided by geological and lunar detections of radioactive 60Fe to assess nearby supernova explosions make now a compelling time for the astronomy community to advocate for supporting multi-disciplinary, cross-cutting research programs.

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the millisecond pulsars PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821-24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the $\sigma_z$ measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities $\sigma_z \approx 3 \times 10^{-14}$ for PSR B1937+21 and on the order of $10^{-12}$ for PSRs B1821$-$24 and J0218+4232. Within the span of our \textit{NICER} data we do not yet see the characteristic break point in the slope of $\sigma_z$; detection of such a break would indicate that further improvement in the cumulative root-mean-square (RMS) timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821-24 and PSR B1937+21 on time scales of $> 2$ years.

No-scale supergravity and the flipped SU(5)$\times$U(1) gauge group provide an ambitious prototype string-inspired scenario for physics below the string scale, which can accommodate the Starobinsky-like inflation favoured by observation when the inflaton is associated with one of the singlet fields associated with neutrino mass generation. During inflation, the vacuum remains in the unbroken GUT phase, and GUT symmetry breaking occurs later when a field with a flat direction (the flaton) acquires a vacuum expectation value. Inflaton decay and the reheating process depend crucially on GUT symmetry breaking, as decay channels open and close, depending on the value of the flaton vacuum expectation value. Here, we consider the simultaneous cosmological evolution of both the inflaton and flaton fields after inflation. We distinguish weak, moderate and strong reheating regimes, and calculate in each case the entropy produced as all fields settle to their global minima. These three reheating scenarios differ in the value of a Yukawa coupling that introduces mass mixing between the singlets and the ${\bf 10}$s of SU(5). The dynamics of the GUT transition has an important impact on the production of gravitinos, and we also discuss the pattern of neutrino masses we expect in each of the three cases. Finally, we use recent CMB limits on neutrino masses to constrain the reheating models, finding that neutrino masses and the cosmological baryon asymmetry can both be explained if the reheating is strong.

Big-Bang nucleosynthesis (BBN) represents one of the earliest phenomena that can lead to observational constraints on the early Universe properties. It is well-known that many important mechanisms and phase transitions occurred before BBN. We discuss the possibility of gaining insight into the primordial Universe through studies of dark matter in cosmology, astroparticle physics and colliders. For this purpose, we assume that dark matter is a thermal relic, and show that combining collider searches with dark matter observables can lead to strong constraints on the period of freeze-out before BBN.

We revisit the uncertainties in the calculation of spin-independent scattering matrix elements for the scattering of WIMP dark matter particles on nuclear matter. In addition to discussing the uncertainties due to limitations in our knowledge of the nucleonic matrix elements of the light quark scalar densities < N |{\bar u} u, {\bar d} d, {\bar s} s| N>, we also discuss the importances of heavy quark scalar densities < N |{\bar c} c, {\bar b} b, {\bar t} t| N >, and comment on uncertainties in quark mass ratios. We analyze estimates of the light-quark densities made over the past decade using lattice calculations and/or phenomenological inputs. We find an uncertainty in the combination < N |{\bar u} u + {\bar d} d | N > that is larger than has been assumed in some phenomenological analyses, and a range of < N |{\bar s} s| N > that is smaller but compatible with earlier estimates. We also analyze the importance of the {\cal O}(\alpha_s^3) calculations of the heavy-quark matrix elements that are now available, which provide an important refinement of the calculation of the spin-independent scattering cross section. We use for illustration a benchmark CMSSM point in the focus-point region that is compatible with the limits from LHC and other searches.

We report timing results for PSR J2234+0611, a 3.6-ms pulsar in a 32-day, eccentric (e = 0.13) orbit with a helium white dwarf companion discovered as part of the Arecibo Observatory 327 MHz drift scan survey. The precise timing and the eccentric nature of the orbit allow precise measurements of an unusual number of parameters: a) a precise proper motion of 27.10(3) mas/yr and a parallax of 1.05(4) mas resulting in a pulsar distance of 0.95(4) kpc; this allows a precise estimate of the transverse velocity, 123(5) km/s. Together with previously published spectroscopic measurements of the systemic radial velocity, this allows a full 3-D determination of the system's velocity; b) precise measurements of the rate of advance of periastron, which after subtraction of the contribution of the proper motion yields a total system mass of $1.6518^{+0.0033}_{-0.0035}$ solar masses; c) a Shapiro delay measurement, h_3 = $82 \pm 14$ ns despite the orbital inclination not being near 90 deg; combined with the measurement of the total mass, this yields a pulsar mass of $1.353^{+0.014}_{-0.017}$ solar masses and a companion mass of $0.298^{+0.015}_{-0.012}$ solar masses; d) we measure precisely the secular variation of the projected semi-major axis and detect significant annual orbital parallax; together these allow a determination of the full 3-D orbital geometry, including an unambiguous orbital inclination (i = $138.7^{+2.5}_{-2.2}$ deg) and a position angle for the line of nodes (Omega = $44^{+5}_{-4}$ deg). We discuss the component masses to investigate hypotheses previously advanced to explain the origin of eccentric MSPs. The unprecedented determination of the full 3-D position, motion and orbital orientation of the system, plus the precisely measured pulsar and WD mass and the latter's optical detection make this system an unique test of our understanding of white dwarfs and their atmospheres.

The observation by the IceCube Collaboration of a high-energy ($E \gtrsim 200$ TeV) neutrino from the direction of the blazar TXS 0506+056 and the coincident observations of enhanced $\gamma$-ray emissions from the same object by MAGIC and other experiments can be used to set stringent constraints on Lorentz violation in the propagation of neutrinos that is linear in the neutrino energy: $\Delta v = - E/M_1$, where $\Delta v$ is the deviation from the velocity of light, and $M_1$ is an unknown high energy scale to be constrained by experiment. Allowing for a difference in neutrino and photon propagation times of $\sim 10$ days, we find that $M_1 \gtrsim 3 \times 10^{16}$ GeV. This improves on previous limits on linear Lorentz violation in neutrino propagation by many orders of magnitude, and the same is true for quadratic Lorentz violation.

Access to 50 years of data has led to the discovery of pulsar emission and rotation variability on timescales of months and years. Most of this long-term variability has been seen in long-period pulsars, with relatively little focus on recycled millisecond pulsars. We have analyzed a 38-pulsar sub-set of the 45 millisecond pulsars in the NANOGrav 11-year data set, in order to review their pulse profile stability. The most variability, on any timescale, is seen in PSRs J1713+0747, B1937+21 and J2145-0750. The strongest evidence for long-timescale pulse profile changes is seen in PSRs B1937+21 and J1643-1224. We have focused our analyses on these four pulsars in an attempt to elucidate the causes of their profile variability. Effects of scintillation seem to be responsible for the profile modifications of PSR J2145-0750. We see evidence that imperfect polarization calibration contributes to the profile variability of PSRs J1713+0747 and B1937+21, along with radio frequency interference around 2 GHz, but find that propagation effects also have an influence. The changes seen in PSR J1643-1224 have been reported previously, yet elude explanation beyond their astrophysical nature. Regardless of cause, unmodeled pulse profile changes are detrimental to the accuracy of pulsar timing and must be incorporated into the timing models where possible.

We show how knowledge of the cold dark matter (CDM) density can be used, in conjunction with measurements of the parameters of a scenario for beyond the Standard Model (BSM) physics, to provide information about the evolution of the Universe before Big Bang Nucleosynthesis (BBN). As examples of non-standard evolution, we consider models with a scalar field that may decay into BSM particles, and quintessence models. We illustrate our calculations using various supersymmetric models as representatives of classes of BSM scenarios in which the CDM density is either larger or smaller than the observed density when the early Universe is assumed to be radiation-dominated. In the case of a decaying scalar field, we show how the CDM density can constrain the initial scalar density and the reheating temperature after it decays in BSM scenarios that would yield overdense dark matter in standard radiation-dominated cosmology, and how the decays of the scalar field into BSM particles can be constrained in scenarios that would otherwise yield underdense CDM. We also show how the early evolution of the quintessence field can be constrained in BSM scenarios.

We have suggested earlier that D-particles, which are stringy space-time defects predicted in brane-inspired models of the Universe, might constitute a component of dark matter, and that they might contribute to the masses of singlet fermions that could provide another component. Interactions of the quantum-fluctuating D-particles with matter induce vector forces that are mediated by a massless effective U(1) gauge field, the "D-photon", which is distinct from the ordinary photon and has different properties from dark photons. We discuss the form of interactions of D-matter with conventional matter induced by D-photon exchange and calculate their strength, which depends on the density of D-particles. Observations of the hydrogen 21~cm line at redshifts >= 15 can constrain these interactions and the density of D-matter in the early Universe.

Low-frequency gravitational-wave experiments require the highest timing precision from an array of the most stable millisecond pulsars. Several known sources of noise on short timescales in single radio-pulsar observations are well described by a simple model of three components: template-fitting from a finite signal-to-noise ratio, pulse phase/amplitude jitter from single-pulse stochasticity, and scintillation errors from short-timescale interstellar scattering variations. Currently template-fitting errors dominate, but as radio telescopes push towards higher signal-to-noise ratios, jitter becomes the next dominant term for most millisecond pulsars. Understanding the statistics of jitter becomes crucial for properly characterizing arrival-time uncertainties. We characterize the radio-frequency dependence of jitter using data on 48 pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing program. We detect significant jitter in 43 of the pulsars and test several functional forms for its frequency dependence; we find significant frequency dependence for 30 pulsars. We find moderate correlations of rms jitter with pulse width (R = 0.62) and number of profile components (R = 0.40); the single-pulse rms jitter is typically ~1% of pulse phase. The average frequency dependence for all pulsars using a power-law model has index -0.42. We investigate the jitter variations for the interpulse of PSR~B1937+21 and find no significant deviations from the main pulse rms jitter. We also test the time-variation of jitter in two pulsars and find that systematics likely bias the results for high-precision pulsars. Pulsar timing array analyses must properly model jitter as a significant component of the noise within the detector.

The well-known Bayes theorem assumes that a posterior distribution is a probability distribution. However, the posterior distribution may no longer be a probability distribution if an improper prior distribution (non-probability measure) such as an unbounded uniform prior is used. Improper priors are often used in the astronomical literature to reflect a lack of prior knowledge, but checking whether the resulting posterior is a probability distribution is sometimes neglected. It turns out that 23 articles out of 75 articles (30.7%) published online in two renowned astronomy journals (ApJ and MNRAS) between Jan 1, 2017 and Oct 15, 2017 make use of Bayesian analyses without rigorously establishing posterior propriety. A disturbing aspect is that a Gibbs-type Markov chain Monte Carlo (MCMC) method can produce a seemingly reasonable posterior sample even when the posterior is not a probability distribution (Hobert and Casella, 1996). In such cases, researchers may erroneously make probabilistic inferences without noticing that the MCMC sample is from a non-existing probability distribution. We review why checking posterior propriety is fundamental in Bayesian analyses, and discuss how to set up scientifically motivated proper priors.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has observed dozens of millisecond pulsars for over a decade. We have accrued a large collection of dispersion measure (DM) measurements sensitive to the total electron content between Earth and the pulsars at each observation. All lines of sight cross through the solar wind which produces correlated DM fluctuations in all pulsars. We develop and apply techniques for extracting the imprint of the solar wind from the full collection of DM measurements in the recently released NANOGrav 11-yr data set. We filter out long time scale DM fluctuations attributable to structure in the interstellar medium and carry out a simultaneous analysis of all pulsars in our sample that can differentiate the correlated signature of the wind from signals unique to individual lines of sight. When treating the solar wind as spherically symmetric and constant in time, we find the electron number density at 1 A.U. to be $7.9\pm0.2$ cm$^{-3}$. Our data shows little evidence of long-term variation in the density of the wind. We argue that our techniques paired with a high cadence, low radio frequency observing campaign of near-ecliptic pulsars would be capable of mapping out large-scale latitudinal structure in the wind.

Observations have revealed that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, these SMBHs form SMBH binaries (SMBHBs) that emit low-frequency gravitational waves (GWs). The incoherent superposition of these sources produce a stochastic GW background (GWB) that can be observed by pulsar timing arrays (PTAs). The optimal statistic is a frequentist estimator of the amplitude of the GWB that specifically looks for the spatial correlations between pulsars induced by the GWB. In this paper, we introduce an improved method for computing the optimal statistic that marginalizes over the red noise in individual pulsars. We use simulations to demonstrate that this method more accurately determines the strength of the GWB, and we use the noise-marginalized optimal statistic to compare the significance of monopole, dipole, and Hellings-Downs (HD) spatial correlations and perform sky scrambles.

The frequency dependence of radio pulse arrival times provides a probe of structures in the intervening media. Demorest et al. 2013 was the first to show a short-term (~100-200 days) reduction in the electron content along the line of sight to PSR J1713+0747 in data from 2008 (approximately MJD 54750) based on an apparent dip in the dispersion measure of the pulsar. We report on a similar event in 2016 (approximately MJD 57510), with average residual pulse-arrival times of approximately 3.0,-1.3, and -0.7 microseconds at 820, 1400, and 2300 MHz, respectively. Timing analyses indicate possible departures from the standard nu^-2 dispersive-delay dependence. We discuss and rule out a wide variety of potential interpretations. We find the likeliest scenario to be lensing of the radio emission by some structure in the interstellar medium, which causes multiple frequency-dependent pulse arrival-time delays.

We present the polarization pulse profiles for 28 pulsars observed with the Arecibo Observatory by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing project at 2.1 GHz, 1.4 GHz, and 430 MHz. These profiles represent some of the most sensitive polarimetric millisecond pulsar profiles to date, revealing the existence of microcomponents (that is, pulse components with peak intensities much lower than the total pulse peak intensity). Although microcomponents have been detected in some pulsars previously, we present microcomponents for PSRs B1937+21, J1713+0747, and J2234+0944 for the first time. These microcomponents can have an impact on pulsar timing, geometry, and flux density determination. We present rotation measures for all 28 pulsars, determined independently at different observation frequencies and epochs, and find the Galactic magnetic fields derived from these rotation measures to be consistent with current models. These polarization profiles were made using measurement equation template matching, which allows us to generate the polarimetric response of the Arecibo Observatory on an epoch-by-epoch basis. We use this method to describe its time variability, and find that the polarimetric responses of the Arecibo Observatory's 1.4 and 2.1 GHz receivers vary significantly with time.

We study possible effects of a dark matter (DM) core on the maximum mass of a neutron star (NS), on the mass-radius relation and on the NS tidal deformability parameter $\Lambda$. We show that all these quantities would in general be reduced in the presence of a DM core. In particular, our calculations indicate that the presence of a DM core with a mass fraction $\sim 5\%$ could affect significantly the interpretation of these NS data as constraints on the nuclear equation of state (EOS), potentially excluding some EOS models on the basis of the measured mass of PSR J0348+0432, while allowing other EOS models to become consistent with the LIGO/Virgo upper limit on $\Lambda$. Specific scenarios for generating such DM cores are explored in an Appendix.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $\Omega_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.

Motivated by the recent detection of the gravitational wave signal emitted by a binary neutron star merger, we analyse the possible impact of dark matter on such signals. We show that dark matter cores in merging neutron stars may yield an observable supplementary peak in the gravitational wave power spectral density following the merger, which could be distinguished from the features produced by the neutron components.

We use MasterCode to perform a frequentist analysis of the constraints on a phenomenological MSSM model with 11 parameters, the pMSSM11, including constraints from ~ 36/fb of LHC data at 13 TeV and PICO, XENON1T and PandaX-II searches for dark matter scattering, as well as previous accelerator and astrophysical measurements, presenting fits both with and without the $(g-2)_{\mu}$ constraint. The pMSSM11 is specified by the following parameters: 3 gaugino masses $M_{1,2,3}$, a common mass for the first-and second-generation squarks $m_{\tilde{q}}$ and a distinct third-generation squark mass $m_{\tilde{q}_3}$, a common mass for the first-and second-generation sleptons $m_{\tilde l}$ and a distinct third-generation slepton mass $m_{\tilde \tau}$, a common trilinear mixing parameter $A$, the Higgs mixing parameter $\mu$, the pseudoscalar Higgs mass $M_A$ and $\tan\beta$. In the fit including $(g-2)_{\mu}$, a Bino-like $\tilde\chi^0_1$ is preferred, whereas a Higgsino-like $\tilde \chi^0_1$ is favoured when the $(g-2)_{\mu}$ constraint is dropped. We identify the mechanisms that operate in different regions of the pMSSM11 parameter space to bring the relic density of the lightest neutralino, $\tilde\chi^0_1$, into the range indicated by cosmological data. In the fit including $(g-2)_{\mu}$, coannihilations with $\tilde \chi^0_2$ and the Wino-like $\tilde\chi^{\pm}_1$ or with nearly-degenerate first- and second-generation sleptons are favoured, whereas coannihilations with the $\tilde \chi^0_2$ and the Higgsino-like $\tilde\chi^{\pm}_1$ or with first- and second-generation squarks may be important when the $(g-2)_{\mu}$ constraint is dropped. Prospects remain for discovering strongly-interacting sparticles at the LHC as well as for discovering electroweakly-interacting sparticles at a future linear $e^+ e^-$ collider such as the ILC or CLIC.

We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, with six high--timing-precision pulsars observed weekly, and all were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and, if necessary, binary parameters, in addition to time-variable dispersion delays and parameters that quantify pulse-profile evolution with frequency. The new timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of large orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. Future papers will use these data to constrain or detect the signatures of gravitational-wave signals.

We describe a likelihood analysis using MasterCode of variants of the MSSM in which the soft supersymmetry-breaking parameters are assumed to have universal values at some scale $M_{in}$ below the supersymmetric grand unification scale $M_{GUT}$, as can occur in mirage mediation and other models. In addition to $M_{in}$, such `sub-GUT' models have the 4 parameters of the CMSSM, namely a common gaugino mass $m_{1/2}$, a common soft supersymmetry-breaking scalar mass $m_0$, a common trilinear mixing parameter $A$ and the ratio of MSSM Higgs vevs $\tan\beta$, assuming that the Higgs mixing parameter $\mu > 0$. We take into account constraints on strongly- and electroweakly-interacting sparticles from $\sim 36$/fb of LHC data at 13 TeV and the LUX and 2017 PICO, XENON1T and PandaX-II searches for dark matter scattering, in addition to the previous LHC and dark matter constraints as well as full sets of flavour and electroweak constraints. We find a preference for $M_{in} \sim 10^5$ to $10^9$ GeV, with $M_{in} \sim M_{GUT}$ disfavoured by $\Delta \chi^2 \sim 3$ due to the ${\rm BR}(B_{s, d} \to \mu^+\mu^-)$ constraint. The lower limits on strongly-interacting sparticles are largely determined by LHC searches, and similar to those in the CMSSM. We find a preference for the LSP to be a Bino or Higgsino with $\tilde{\chi^0_1} \sim 1$ TeV, with annihilation via heavy Higgs bosons $H/A$ and stop coannihilation, or chargino coannihilation, bringing the cold dark matter density into the cosmological range. We find that spin-independent dark matter scattering is likely to be within reach of the planned LUX-Zeplin and XENONnT experiments. We probe the impact of the $(g-2)_\mu$ constraint, finding similar results whether or not it is included.

We construct a no-scale model of inflation with a single modulus whose real and imaginary parts are fixed by simple power-law corrections to the no-scale K{\" a}hler potential. Assuming an uplift of the minimum of the effective potential, the model yields a suitable number of e-folds of expansion and values of the tilt in the scalar cosmological density perturbations and of the ratio of tensor and scalar perturbations that are compatible with measurements of the cosmic microwave background radiation.

Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|\Delta|< 0.002$ for the UFF test and $-3\times10^{-20} < \hat{\alpha}_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hat{\alpha}_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests.

We show that $R^2$ gravity coupled conformally to scalar fields is equivalent to the real bosonic sector of SU(N,1)/SU(N)$\times$U(1) no-scale supergravity, where the conformal factor can be identified with the K\"ahler potential, and we review the construction of Starobinsky-like models of inflation within this framework.

In-spiraling supermassive black holes should emit gravitational waves, which would produce characteristic distortions in the time of arrival residuals from millisecond pulsars. Multiple national and regional consortia have constructed pulsar timing arrays by precise timing of different sets of millisecond pulsars. An essential aspect of precision timing is the transfer of the times of arrival to a (quasi-)inertial frame, conventionally the solar system barycenter. The barycenter is determined from the knowledge of the planetary masses and orbits, which has been refined over the past 50 years by multiple spacecraft. Within the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), uncertainties on the solar system barycenter are emerging as an important element of the NANOGrav noise budget. We describe what is known about the solar system barycenter, touch upon how uncertainties in it affect gravitational wave studies with pulsar timing arrays, and consider future trends in spacecraft navigation.

Gravitational wave (GW) detection with pulsar timing arrays (PTAs) requires accurate noise characterization. The noise of our Galactic-scale GW detector has been systematically evaluated by the Noise Budget and Interstellar Medium Mitigation working groups within the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. Intrinsically, individual radio millisecond pulsars (MSPs) used by NANOGrav can have some degree of achromatic red spin noise, as well as white noise due to pulse phase jitter. Along any given line-of-sight, the ionized interstellar medium contributes chromatic noise through dispersion measure (DM) variations, interstellar scintillation, and scattering. These effects contain both red and white components. In the future, with wideband receivers, the effects of frequency-dependent DM will become important. Having anticipated and measured these diverse sources of detector noise, the NANOGrav PTA remains well-poised to detect low-frequency GWs.

Supermassive black hole binaries (SMBHBs) in the 10 million to 10 billion $M_\odot$ range form in galaxy mergers, and live in galactic nuclei with large and poorly constrained concentrations of gas and stars. There are currently no observations of merging SMBHBs--- it is in fact possible that they stall at their final parsec of separation and never merge. While LIGO has detected high frequency GWs, SMBHBs emit GWs in the nanohertz to millihertz band. This is inaccessible to ground-based interferometers, but possible with Pulsar Timing Arrays (PTAs). Using data from local galaxies in the 2 Micron All-Sky Survey, together with galaxy merger rates from Illustris, we find that there are on average $91\pm7$ sources emitting GWs in the PTA band, and $7\pm2$ binaries which will never merge. Local unresolved SMBHBs can contribute to GW background anisotropy at a level of $\sim20\%$, and if the GW background can be successfully isolated, GWs from at least one local SMBHB can be detected in 10 years.

We embed a flipped ${\rm SU}(5) \times {\rm U}(1)$ GUT model in a no-scale supergravity framework, and discuss its predictions for cosmic microwave background observables, which are similar to those of the Starobinsky model of inflation. Measurements of the tilt in the spectrum of scalar perturbations in the cosmic microwave background, $n_s$, constrain significantly the model parameters. We also discuss the model's predictions for neutrino masses, and pay particular attention to the behaviours of scalar fields during and after inflation, reheating and the GUT phase transition. We argue in favor of strong reheating in order to avoid excessive entropy production which could dilute the generated baryon asymmetry.

We perform a likelihood analysis of the constraints from accelerator experiments and astrophysical observations on supersymmetric (SUSY) models with SU(5) boundary conditions on soft SUSY-breaking parameters at the GUT scale. The parameter space of the models studied has 7 parameters: a universal gaugino mass $m_{1/2}$, distinct masses for the scalar partners of matter fermions in five- and ten-dimensional representations of SU(5), $m_5$ and $m_{10}$, and for the $\mathbf{5}$ and $\mathbf{\bar 5}$ Higgs representations $m_{H_u}$ and $m_{H_d}$, a universal trilinear soft SUSY-breaking parameter $A_0$, and the ratio of Higgs vevs $\tan \beta$. In addition to previous constraints from direct sparticle searches, low-energy and flavour observables, we incorporate constraints based on preliminary results from 13 TeV LHC searches for jets + MET events and long-lived particles, as well as the latest PandaX-II and LUX searches for direct Dark Matter detection. In addition to previously-identified mechanisms for bringing the supersymmetric relic density into the range allowed by cosmology, we identify a novel ${\tilde u_R}/{\tilde c_R} - \tilde{\chi}^0_1$ coannihilation mechanism that appears in the supersymmetric SU(5) GUT model and discuss the role of ${\tilde \nu_\tau}$ coannihilation. We find complementarity between the prospects for direct Dark Matter detection and SUSY searches at the LHC.

The photon mass, $m_\gamma$, can in principle be constrained using measurements of the dispersion measures (DMs) of fast radio bursts (FRBs), once the FRB redshifts are known. The DM of the repeating FRB 121102 is known to $< 1$\%, a host galaxy has now been identified with high confidence,and its redshift, $z$, has now been determined with high accuracy: $z = 0.19273(8)$. Taking into account the plasma contributions to the DM from the Intergalactic medium (IGM) and the Milky Way, we use the data on FRB 121102 to derive the constraint $m_\gamma \lesssim 2.2 \times 10^{-14}$ eV c$^{-2}$ ($3.9 \times 10^{-50}$ kg). Since the plasma and photon mass contributions to DMs have different redshift dependences, they could in principle be distinguished by measurements of more FRB redshifts, enabling the sensitivity to $m_\gamma$ to be improved.

We reconsider the minimal SU(5) Grand Unified Theory (GUT) in the context of no-scale supergravity, assuming that the soft supersymmetry-breaking parameters satisfy universality conditions at some input scale M_in above the GUT scale M_GUT. When setting up such a no-scale super-GUT model, special attention must be paid to avoiding the Scylla of rapid proton decay and the Charybdis of an excessive density of cold dark matter, while also having an acceptable mass for the Higgs boson. We do not find consistent solutions if none of the matter and Higgs fields are assigned to twisted chiral supermultiplets, even in the presence of Giudice-Masiero terms. However, consistent solutions may be found if at least one fiveplet of GUT Higgs fields is assigned to a twisted chiral supermultiplet, with a suitable choice of modular weights. Spin-independent dark matter scattering may be detectable in some of these consistent solutions.

We perform a likelihood analysis of the minimal Anomaly-Mediated Supersymmetry Breaking (mAMSB) model using constraints from cosmology and accelerator experiments. We find that a wino-like or a Higgsino-like neutralino LSP, $m_{\tilde \chi^0_{1}}$, may provide the cold dark matter (DM) with similar likelihood. The upper limit on the DM density from Planck and other experiments enforces $m_{\tilde \chi^0_{1}} \lesssim 3~TeV$ after the inclusion of Sommerfeld enhancement in its annihilations. If most of the cold DM density is provided by the $\tilde \chi_0^1$, the measured value of the Higgs mass favours a limited range of $\tan \beta \sim 5$ (or for $\mu > 0$, $\tan \beta \sim 45$) but the scalar mass $m_0$ is poorly constrained. In the wino-LSP case, $m_{3/2}$ is constrained to about $900~TeV$ and ${m_{\tilde \chi^0_{1}}}$ to $2.9\pm0.1~TeV$, whereas in the Higgsino-LSP case $m_{3/2}$ has just a lower limit $\gtrsim 650TeV$ ($\gtrsim 480TeV$) and $m_{\tilde \chi^0_{1}}$ is constrained to $1.12 ~(1.13) \pm0.02~TeV$ in the $\mu>0$ ($\mu<0$) scenario. In neither case can the anomalous magnetic moment of the muon, ${(g-2)_\mu}$, be improved significantly relative to its Standard Model (SM) value, nor do flavour measurements constrain the model significantly, and there are poor prospects for discovering supersymmetric particles at the LHC, {though there} are some prospects for direct DM detection. On the other hand, if the ${m_{\tilde \chi^0_{1}}}$ contributes only a fraction of the cold DM density, {future LHC $E_T$-based searches for gluinos, squarks and heavier chargino and neutralino states as well as disappearing track searches in the wino-like LSP region will be relevant}, and interference effects enable ${\rm BR}(B_{s, d} \to \mu^+\mu^-)$ to agree with the data better than in the SM in the case of wino-like DM with $\mu > 0$.

We analyze dispersion measure (DM) variations of 37 millisecond pulsars in the 9-year NANOGrav data release and constrain the sources of these variations. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 13 pulsars show both effects. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized ISM features and, in one case, an upper limit to the size of the dispersing region of 13.2 AU. Finally, five pulsars show very nearly quadratic structure functions, which could be indicative of an underlying Kolmogorov medium. Four pulsars show roughly Kolmogorov structure functions and another four show structure functions less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and show that the presence of other trends in the data is the most likely cause.

Gravitational wave astronomy using a pulsar timing array requires high-quality millisecond pulsars, correctable interstellar propagation delays, and high-precision measurements of pulse times of arrival. Here we identify noise in timing residuals that exceeds that predicted for arrival time estimation for millisecond pulsars observed by the North American Nanohertz Observatory for Gravitational Waves. We characterize the excess noise using variance and structure function analyses. We find that 26 out of 37 pulsars show inconsistencies with a white-noise-only model based on the short timescale analysis of each pulsar and we demonstrate that the excess noise has a red power spectrum for 15 pulsars. We also decompose the excess noise into chromatic (radio-frequency-dependent) and achromatic components. Associating the achromatic red-noise component with spin noise and including additional power-spectrum-based estimates from the literature, we estimate a scaling law in terms of spin parameters (frequency and frequency derivative) and data-span length and compare it to the scaling law of Shannon \& Cordes (2010). We briefly discuss our results in terms of detection of gravitational waves at nanohertz frequencies.

In pulsar timing, timing residuals are the differences between the observed times of arrival and the predictions from the timing model. A comprehensive timing model will produce featureless residuals, which are presumably composed of dominating noise and weak physical effects excluded from the timing model (e.g. gravitational waves). In order to apply the optimal statistical methods for detecting the weak gravitational wave signals, we need to know the statistical properties of the noise components in the residuals. In this paper we utilize a variety of non-parametric statistical tests to analyze the whiteness and Gaussianity of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 5-year timing data which are obtained from the Arecibo Observatory and the Green Bank Telescope from 2005 to 2010 (Demorest et al. 2013). We find that most of the data are consistent with white noise; Many data deviate from Gaussianity at different levels, nevertheless, removing outliers in some pulsars will mitigate the deviations.

We analyze 24 binary radio pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) nine-year data set. We make fourteen significant measurements of Shapiro delay, including new detections in four pulsar-binary systems (PSRs J0613$-$0200, J2017+0603, J2302+4442, and J2317+1439), and derive estimates of the binary-component masses and orbital inclination for these MSP-binary systems. We find a wide range of binary pulsar masses, with values as low as $m_{\rm p} = 1.18^{+0.10}_{-0.09}\text{M}_{\odot}$ for PSR J1918$-$0642 and as high as $m_{\rm p} = 1.928^{+0.017}_{-0.017}\text{M}_{\odot}$ for PSR J1614$-$2230 (both 68.3\% credibility). We make an improved measurement of the Shapiro timing delay in the PSR J1918$-$0642 and J2043+1711 systems, measuring the pulsar mass in the latter system to be $m_{\rm p} = 1.41^{+0.21}_{-0.18}\text{M}_{\odot}$ (68.3\% credibility) for the first time. We measure secular variations of one or more orbital elements in many systems, and use these measurements to further constrain our estimates of the pulsar and companion masses whenever possible. In particular, we used the observed Shapiro delay and periastron advance due to relativistic gravity in the PSR J1903+0327 system to derive a pulsar mass of $m_{\rm p} = 1.65^{+0.02}_{-0.02}\text{M}_{\odot}$ (68.3\% credibility). We discuss the implications that our mass measurements have on the overall neutron-star mass distribution, and on the "mass/orbital-period" correlation due to extended mass transfer.

The use of pulsars as astrophysical clocks for gravitational wave experiments demands the highest possible timing precision. Pulse times of arrival (TOAs) are limited by stochastic processes that occur in the pulsar itself, along the line of sight through the interstellar medium, and in the measurement process. On timescales of seconds to hours, the TOA variance exceeds that from template-fitting errors due to additive noise. We assess contributions to the total variance from two additional effects: amplitude and phase jitter intrinsic to single pulses and changes in the interstellar impulse response from scattering. The three effects have different dependencies on time, frequency, and pulse signal-to-noise ratio. We use data on 37 pulsars from the North American Nanohertz Observatory for Gravitational Waves to assess the individual contributions to the overall intraday noise budget for each pulsar. We detect jitter in 22 pulsars and estimate the average value of RMS jitter in our pulsars to be $\sim 1\%$ of pulse phase. We examine how jitter evolves as a function of frequency and find evidence for evolution. Finally, we compare our measurements with previous noise parameter estimates and discuss methods to improve gravitational wave detection pipelines.

Using a no-scale supergravity framework, we construct an SO(10) model that makes predictions for cosmic microwave background observables similar to those of the Starobinsky model of inflation, and incorporates a double-seesaw model for neutrino masses consistent with oscillation experiments and late-time cosmology. We pay particular attention to the behaviour of the scalar fields during inflation and the subsequent reheating.

Dense, continuous pulsar timing observations over a 24-hr period provide a method for probing intermediate gravitational wave (GW) frequencies from 10 microhertz to 20 millihertz. The European Pulsar Timing Array (EPTA), the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the Parkes Pulsar Timing Array (PPTA), and the combined International Pulsar Timing Array (IPTA) all use millisecond pulsar observations to detect or constrain GWs typically at nanohertz frequencies. In the case of the IPTA's nine-telescope 24-Hour Global Campaign on millisecond pulsar J1713+0747, GW limits in the intermediate frequency regime can be produced. The negligible change in dispersion measure during the observation minimizes red noise in the timing residuals, constraining any contributions from GWs due to individual sources. At 10$^{-5}$Hz, the 95% upper limit on strain is 10$^{-11}$ for GW sources in the pulsar's direction.

We extend our previous study of supersymmetric Higgs inflation in the context of no-scale supergravity and grand unification, to include models based on the flipped SU(5) and the Pati-Salam group. Like the previous SU(5) GUT model, these yield a class of inflation models whose inflation predictions interpolate between those of the quadratic chaotic inflation and Starobinsky-like inflation, while avoiding tension with proton decay limits. We further analyse the reheating process in these models, and derive the number of e-folds, which is independent of the reheating temperature. We derive the corresponding predictions for the scalar tilt and the tensor-to-scalar ratio in cosmic microwave background perturbations, as well as discussing the gravitino production following inflation.

PSR J1024$-$0719 is a millisecond pulsar that was long thought to be isolated. However, puzzling results concerning its velocity, distance, and low rotational period derivative have led to reexamination of its properties. We present updated radio timing observations along with new and archival optical data that show PSR J1024$-$0719 is most likely in a long period (2$-$20 kyr) binary system with a low-mass ($\approx 0.4\,M_\odot$) low-metallicity ($Z \approx -0.9\,$ dex) main sequence star. Such a system can explain most of the anomalous properties of this pulsar. We suggest that this system formed through a dynamical exchange in a globular cluster that ejected it into a halo orbit, consistent with the low observed metallicity for the stellar companion. Further astrometric and radio timing observations such as measurement of the third period derivative could strongly constrain the range of orbital parameters.

The observation of gravitational waves from the Laser Interferometer Gravitational-Wave Observatory (LIGO) event GW150914 may be used to constrain the possibility of Lorentz violation in graviton propagation, and the observation by the Fermi Gamma-Ray Burst Monitor of a transient source in apparent coincidence may be used to constrain the difference between the velocities of light and gravitational waves: $c_g - c_\gamma < 10^{-17}$.

Several searches have found evidence of $^{60}$Fe deposition, presumably from a near-Earth supernova (SN), with concentrations that vary in different locations on Earth. This paper examines various influences on the path of interstellar dust carrying $^{60}$Fe from a SN through the heliosphere, with the aim of estimating the final global distribution on the ocean floor. We study the influences of magnetic fields, angle of arrival, wind and ocean cycling of SN material on the concentrations at different locations. We find that the passage of SN material through the mesosphere/lower thermosphere (MLT) is the greatest influence on the final global distribution, with ocean cycling causing lesser alteration as the SN material sinks to the ocean floor. SN distance estimates in previous works that assumed a uniform distribution are a good approximation. Including the effects on surface distributions, we estimate a distance of $46^{+10}_{-6}$ pc for a $8-10 \ M_{\odot}$ SN progenitor. This is consistent with a SN occurring within the Tuc-Hor stellar group $\sim$2.8 Myr ago with SN material arriving on Earth $\sim$2.2 Myr ago. We note that the SN dust retains directional information to within $1^{\circ}$ through its arrival in the inner Solar System, so that SN debris deposition on inert bodies such as the Moon will be anisotropic, and thus could in principle be used to infer directional information. In particular, we predict that existing lunar samples should show measurable $^{60}$Fe differences.

The frequency-dependent time delays in fast radio bursts (FRBs) can be used to constrain the photon mass, if the FRB redshifts are known, but the similarity between the frequency dependences of dispersion due to plasma effects and a photon mass complicates the derivation of a limit on $m_\gamma$. The dispersion measure (DM) of FRB 150418 is known to $\sim 0.1$%, and there is a claim to have measured its redshift with an accuracy of $\sim 2$%, but the strength of the constraint on $m_\gamma$ is limited by uncertainties in the modelling of the host galaxy and the Milky Way, as well as possible inhomogeneities in the intergalactic medium (IGM). Allowing for these uncertainties, the recent data on FRB 150418 indicate that $m_\gamma \lesssim 1.8 \times 10^{-14}$eV c$^{-2}$ ($3.2 \times 10^{-50}$kg), if FRB 150418 indeed has a redshift $z = 0.492$ as initially reported. In the future, the different redshift dependences of the plasma and photon mass contributions to DM can be used to improve the sensitivity to $m_\gamma$ if more FRB redshifts are measured. For a fixed fractional uncertainty in the extra-galactic contribution to the DM of an FRB, one with a lower redshift would provide greater sensitivity to $m_\gamma$.

Decade-long timing observations of arrays of millisecond pulsars have placed highly constraining upper limits on the amplitude of the nanohertz gravitational-wave stochastic signal from the mergers of supermassive black-hole binaries ($\sim 10^{-15}$ strain at $f = 1/\mathrm{yr}$). These limits suggest that binary merger rates have been overestimated, or that environmental influences from nuclear gas or stars accelerate orbital decay, reducing the gravitational-wave signal at the lowest, most sensitive frequencies. This prompts the question whether nanohertz gravitational waves are likely to be detected in the near future. In this letter, we answer this question quantitatively using simple statistical estimates, deriving the range of true signal amplitudes that are compatible with current upper limits, and computing expected detection probabilities as a function of observation time. We conclude that small arrays consisting of the pulsars with the least timing noise, which yield the tightest upper limits, have discouraging prospects of making a detection in the next two decades. By contrast, we find large arrays are crucial to detection because the quadrupolar spatial correlations induced by gravitational waves can be well sampled by many pulsar pairs. Indeed, timing programs which monitor a large and expanding set of pulsars have an $\sim 80\%$ probability of detecting gravitational waves within the next ten years, under assumptions on merger rates and environmental influences ranging from optimistic to conservative. Even in the extreme case where $90\%$ of binaries stall before merger and environmental coupling effects diminish low-frequency gravitational-wave power, detection is delayed by at most a few years.

Cosmology and particle physics have long been dominated by theoretical paradigms: Einstein's general theory of relativity in cosmology and the Standard Model of particle physics. The time may have come for paradigm shifts. Does cosmological inflation require a modification of Einstein's gravity? Have experiments at the LHC discovered a new particle beyond the Standard Model? It is premature to answer these questions, but we theorists can dream about the possibilities.

We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.

The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.

We report on an effort to extract and monitor interstellar scintillation parameters in regular timing observations collected for the NANOGrav pulsar timing array. Scattering delays are measured by creating dynamic spectra for each pulsar and observing epoch of wide-band observations centered near 1500 MHz and carried out at the Green Bank Telescope and the Arecibo Observatory. The ~800-MHz wide frequency bands imply dramatic changes in scintillation bandwidth across the bandpass, and a stretching routine has been included to account for this scaling. For most of the 10 pulsars for which the scaling has been measured, the bandwidths scale with frequency less steeply than expected for a Kolmogorov medium. We find estimated scattering delay values that vary with time by up to an order of magnitude. The mean measured scattering delays are similar to previously published values and slightly higher than predicted by interstellar medium models. We investigate the possibility of increasing the timing precision by mitigating timing errors introduced by the scattering delays. For most of the pulsars, the uncertainty in the time of arrival of a single timing point is much larger than the maximum variation of the scattering delay, suggesting that diffractive scintillation remains only a negligible part of their noise budget.

We present high-precision timing observations spanning up to nine years for 37 millisecond pulsars monitored with the Green Bank and Arecibo radio telescopes as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We describe the observational and instrumental setups used to collect the data, and methodology applied for calculating pulse times of arrival; these include novel methods for measuring instrumental offsets and characterizing low signal-to-noise ratio timing results. The time of arrival data are fit to a physical timing model for each source, including terms that characterize time-variable dispersion measure and frequency-dependent pulse shape evolution. In conjunction with the timing model fit, we have performed a Bayesian analysis of a parameterized timing noise model for each source, and detect evidence for excess low-frequency, or "red," timing noise in 10 of the pulsars. For 5 of these cases this is likely due to interstellar medium propagation effects rather than intrisic spin variations. Subsequent papers in this series will present further analysis of this data set aimed at detecting or limiting the presence of nanohertz-frequency gravitational wave signals.

The modeling of intrinsic noise in pulsar timing residual data is of crucial importance for Gravitational Wave (GW) detection and pulsar timing (astro)physics in general. The noise budget in pulsars is a collection of several well studied effects including radiometer noise, pulse-phase jitter noise, dispersion measure (DM) variations, and low frequency spin noise. However, as pulsar timing data continues to improve, non-stationary and non-powerlaw noise terms are beginning to manifest which are not well modeled by current noise analysis techniques. In this work we use a trans-dimensional approach to model these non-stationary and non-powerlaw effects through the use of a wavelet basis and an interpolation based adaptive spectral modeling. In both cases, the number of wavelets and the number of control points in the interpolated spectrum are free parameters that are constrained by the data and then marginalized over in the final inferences, thus fully incorporating our ignorance of the noise model. We show that these new methods outperform standard techniques when non-stationary and non-powerlaw noise is present. We also show that these methods return results consistent with the standard analyses when no such signals are present.

Using the nine-year radio-pulsar timing data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), collected at Arecibo Observatory and the Green Bank Telescope, we have measured the positions, proper motions, and parallaxes for 37 millisecond pulsars. We report twelve significant parallax measurements and distance measurements, and eighteen lower limits on distance. We compare these measurements to distances predicted by the NE2001 interstellar electron density model and find them to be in general agreement. We use measured orbital-decay rates and spin-down rates to confirm two of the parallax distances and to place distance upper limits on other sources; these distance limits agree with the parallax distances with one exception, PSR J1024-0719, which we discuss at length. Using the proper motions of the 37 NANOGrav pulsars in combination with other published measurements, we calculate the velocity dispersion of the millisecond pulsar population in Galactocentric coordinates. We find the radial, azimuthal, and perpendicular dispersions to be 46, 40, and 24 km s-1, respectively, in a model that allows for high-velocity outliers; or 81, 58, and 62 km s-1 for the full population. These velocity dispersions are far smaller than those of the canonical pulsar population, and are similar to older Galactic disk populations. This suggests that millisecond pulsar velocities are largely attributable to their being an old population rather than being artifacts of their birth and evolution as neutron star binary systems. The components of these velocity dispersions follow similar proportions to other Galactic populations, suggesting that our results are not biased by selection effects.

We revisit gravitino production following inflation. As a first step, we review the standard calculation of gravitino production in the thermal plasma formed at the end of post-inflationary reheating when the inflaton has completely decayed. Next we consider gravitino production prior to the completion of reheating, assuming that the inflaton decay products thermalize instantaneously while they are still dilute. We then argue that instantaneous thermalization is in general a good approximation, and also show that the contribution of non-thermal gravitino production via the collisions of inflaton decay products prior to thermalization is relatively small. Our final estimate of the gravitino-to-entropy ratio is approximated well by a standard calculation of gravitino production in the post-inflationary thermal plasma assuming total instantaneous decay and thermalization at a time $t \simeq 1.2/\Gamma_\phi$. Finally, in light of our calculations, we consider potential implications of upper limits on the gravitino abundance for models of inflation, with particular attention to scenarios for inflaton decays in supersymmetric Starobinsky-like models.

We present the SuperNova Explosion Code (SNEC), an open-source Lagrangian code for the hydrodynamics and equilibrium-diffusion radiation transport in the expanding envelopes of supernovae. Given a model of a progenitor star, an explosion energy, and an amount and distribution of radioactive nickel, SNEC generates the bolometric light curve, as well as the light curves in different broad bands assuming black body emission. As a first application of SNEC, we consider the explosions of a grid of 15 Msun (at zero-age main sequence) stars whose hydrogen envelopes are stripped to different extents and at different points in their evolution. The resulting light curves exhibit plateaus with durations of ~20-100 days if >~1.5-2 Msun of hydrogen-rich material is left and no plateau if less hydrogen-rich material is left. If these shorter plateau lengths are not seen for Type IIP supernovae in nature, it suggests that, at least for zero-age main sequence masses <~ 20 Msun, hydrogen mass loss occurs as an all or nothing process. This perhaps points to the important role binary interactions play in generating the observed mass-stripped supernovae (i.e., Type Ib/c events). These light curves are also unlike what is typically seen for Type IIL supernovae, arguing that simply varying the amount of mass loss cannot explain these events. The most stripped models begin to show double-peaked light curves similar to what is often seen for Type IIb supernovae, confirming previous work that these supernovae can come from progenitors that have a small amount of hydrogen and a radius of ~500 Rsun.

We consider two potential non-accelerator signatures of generalizations of the well-studied constrained minimal supersymmetric standard model (CMSSM). In one generalization, the universality constraints on soft supersymmetry-breaking parameters are applied at some input scale $M_{in}$ below the grand unification (GUT) scale $M_{GUT}$, a scenario referred to as `sub-GUT'. The other generalization we consider is to retain GUT-scale universality for the squark and slepton masses, but to relax universality for the soft supersymmetry-breaking contributions to the masses of the Higgs doublets. As with other CMSSM-like models, the measured Higgs mass requires supersymmetric particle masses near or beyond the TeV scale. Because of these rather heavy sparticle masses, the embedding of these CMSSM-like models in a minimal SU(5) model of grand unification can yield a proton lifetime consistent with current experimental limits, and may be accessible in existing and future proton decay experiments. Another possible signature of these CMSSM-like models is direct detection of supersymmetric dark matter. The direct dark matter scattering rate is typically below the reach of the LUX-ZEPLIN (LZ) experiment if $M_{in}$ is close to $M_{GUT}$, but may lie within its reach if $M_{in} \lesssim 10^{11}$ GeV. Likewise, generalizing the CMSSM to allow non-universal supersymmetry-breaking contributions to the Higgs offers extensive possibilities for models within reach of the LZ experiment that have long proton lifetimes.

We report 21-yr timing of one of the most precise pulsars: PSR J1713+0747. Its pulse times of arrival are well modeled by a comprehensive pulsar binary model including its three-dimensional orbit and a noise model that incorporates correlated noise such as jitter and red noise. Its timing residuals have weighted root mean square $\sim 92$ ns. The new dataset allows us to update and improve previous measurements of the system properties, including the masses of the neutron star ($1.31\pm0.11$ $M_{\odot}$) and the companion white dwarf ($0.286\pm0.012$ $M_{\odot}$) and the parallax distance $1.15\pm0.03$ kpc. We measured the intrinsic change in orbital period, $\dot{P}^{\rm Int}_{\rm b}$, is $-0.20\pm0.17$ ps s$^{-1}$, which is not distinguishable from zero. This result, combined with the measured $\dot{P}^{\rm Int}_{\rm b}$ of other pulsars, can place a generic limit on potential changes in the gravitational constant $G$. We found that $\dot{G}/G$ is consistent with zero [$(-0.6\pm1.1)\times10^{-12}$ yr$^{-1}$, 95\% confidence] and changes at least a factor of $31$ (99.7\% confidence) more slowly than the average expansion rate of the Universe. This is the best $\dot{G}/G$ limit from pulsar binary systems. The $\dot{P}^{\rm Int}_{\rm b}$ of pulsar binaries can also place limits on the putative coupling constant for dipole gravitational radiation $\kappa_D=(-0.9\pm3.3)\times10^{-4}$ (95\% confidence). Finally, the nearly circular orbit of this pulsar binary allows us to constrain statistically the strong-field post-Newtonian parameters $\Delta$, which describes the violation of strong equivalence principle, and $\hat{\alpha}_3$, which describes a breaking of both Lorentz invariance in gravitation and conservation of momentum. We found, at 95\% confidence, $\Delta<0.01$ and $\hat{\alpha}_3<2\times10^{-20}$ based on PSR J1713+0747.

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. We set upper limits for a GWB from supermassive black hole binaries under power law, broken power law, and free spectral coefficient GW spectrum models. We place a 95\% upper limit on the strain amplitude (at a frequency of yr$^{-1}$) in the power law model of $A_{\rm gw} < 1.5\times 10^{-15}$. For a broken power law model, we place priors on the strain amplitude derived from simulations of Sesana (2013) and McWilliams et al. (2014). We find that the data favor a broken power law to a pure power law with odds ratios of 22 and 2.2 to one for the McWilliams and Sesana prior models, respectively. The McWilliams model is essentially ruled out by the data, and the Sesana model is in tension with the data under the assumption of a pure power law. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. We then place the most stringent limits so far on the energy density of relic GWs, $\Omega_\mathrm{gw}(f)\,h^2 < 4.2 \times 10^{-10}$, yielding a limit on the Hubble parameter during inflation of $H_*=1.6\times10^{-2}~m_{Pl}$, where $m_{Pl}$ is the Planck mass. Our limit on the cosmic string GWB, $\Omega_\mathrm{gw}(f)\, h^2 < 2.2 \times 10^{-10}$, translates to a conservative limit of $G\mu<3.3\times 10^{-8}$ - a factor of 4 better than the joint Planck and high-$l$ CMB data from other experiments.

Different mechanisms operate in various regions of the MSSM parameter space to bring the relic density of the lightest neutralino, neutralino_1, assumed here to be the LSP and thus the Dark Matter (DM) particle, into the range allowed by astrophysics and cosmology. These mechanisms include coannihilation with some nearly-degenerate next-to-lightest supersymmetric particle (NLSP) such as the lighter stau (stau_1), stop (stop_1) or chargino (chargino_1), resonant annihilation via direct-channel heavy Higgs bosons H/A, the light Higgs boson h or the Z boson, and enhanced annihilation via a larger Higgsino component of the LSP in the focus-point region. These mechanisms typically select lower-dimensional subspaces in MSSM scenarios such as the CMSSM, NUHM1, NUHM2 and pMSSM10. We analyze how future LHC and direct DM searches can complement each other in the exploration of the different DM mechanisms within these scenarios. We find that the stau_1 coannihilation regions of the CMSSM, NUHM1, NUHM2 can largely be explored at the LHC via searches for missing E_T events and long-lived charged particles, whereas their H/A funnel, focus-point and chargino_1 coannihilation regions can largely be explored by the LZ and Darwin DM direct detection experiments. We find that the dominant DM mechanism in our pMSSM10 analysis is chargino_1 coannihilation: {parts of its parameter space can be explored by the LHC, and a larger portion by future direct DM searches.

We analyse the prospects for constructing hybrid models of inflation that provide a dynamical realisation of the apparent closeness between the supersymmetric GUT scale and the possible scale of cosmological inflation. In the first place, we consider models based on the flipped SU(5)$\times$U(1) gauge group, which has no magnetic monopoles. In one model, the inflaton is identified with a sneutrino field, and in the other model it is a gauge singlet. In both cases we find regions of the model parameter spaces that are compatible with the experimental magnitudes of the scalar perturbations, $A_s$, and the tilt in the scalar perturbation spectrum, $n_s$, as well as with an indicative upper limit on the tensor-to-scalar perturbation ratio, $r$. We also discuss embeddings of these models into SO(10), which is broken at a higher scale so that its monopoles are inflated away.

Supersymmetry is the most natural framework for physics above the TeV scale, and the corresponding framework for early-Universe cosmology, including inflation, is supergravity. No-scale supergravity emerges from generic string compactifications and yields a non-negative potential, and is therefore a plausible framework for constructing models of inflation. No-scale inflation yields naturally predictions similar to those of the Starobinsky model based on $R + R^2$ gravity, with a tilted spectrum of scalar perturbations: $n_s \sim 0.96$, and small values of the tensor-to-scalar perturbation ratio $r < 0.1$, as favoured by Planck and other data on the cosmic microwave background (CMB). Detailed measurements of the CMB may provide insights into the embedding of inflation within string theory as well as its links to collider physics.

We discuss inflaton decays and reheating in no-scale Starobinsky-like models of inflation, calculating the effective equation-of-state parameter, $w$, during the epoch of inflaton decay, the reheating temperature, $T_{\rm reh}$, and the number of inflationary e-folds, $N_*$, comparing analytical approximations with numerical calculations. We then illustrate these results with applications to models based on no-scale supergravity and motivated by generic string compactifications, including scenarios where the inflaton is identified as an untwisted-sector matter field with direct Yukawa couplings to MSSM fields, and where the inflaton decays via gravitational-strength interactions. Finally, we use our results to discuss the constraints on these models imposed by present measurements of the scalar spectral index $n_s$ and the tensor-to-scalar perturbation ratio $r$, converting them into constraints on $N_*$, the inflaton decay rate and other parameters of specific no-scale inflationary models.

What is the Higgs boson telling us? What else is there, maybe supersymmetry and/or dark matter? How do we find it? These are now the big questions in collider physics that I discuss in this talk, from a personal point of view.

We present a frequentist analysis of the parameter space of the pMSSM10, in which the following 10 soft SUSY-breaking parameters are specified independently at the mean scalar top mass scale Msusy = Sqrt[M_stop1 M_stop2]: the gaugino masses M_{1,2,3}, the 1st-and 2nd-generation squark masses M_squ1 = M_squ2, the third-generation squark mass M_squ3, a common slepton mass M_slep and a common trilinear mixing parameter A, the Higgs mixing parameter mu, the pseudoscalar Higgs mass M_A and tan beta. We use the MultiNest sampling algorithm with 1.2 x 10^9 points to sample the pMSSM10 parameter space. A dedicated study shows that the sensitivities to strongly-interacting SUSY masses of ATLAS and CMS searches for jets, leptons + MET signals depend only weakly on many of the other pMSSM10 parameters. With the aid of the Atom and Scorpion codes, we also implement the LHC searches for EW-interacting sparticles and light stops, so as to confront the pMSSM10 parameter space with all relevant SUSY searches. In addition, our analysis includes Higgs mass and rate measurements using the HiggsSignals code, SUSY Higgs exclusion bounds, the measurements B-physics observables, EW precision observables, the CDM density and searches for spin-independent DM scattering. We show that the pMSSM10 is able to provide a SUSY interpretation of (g-2)_mu, unlike the CMSSM, NUHM1 and NUHM2. As a result, we find (omitting Higgs rates) that the minimum chi^2/dof = 20.5/18 in the pMSSM10, corresponding to a chi^2 probability of 30.8 %, to be compared with chi^2/dof = 32.8/24 (31.1/23) (30.3/22) in the CMSSM (NUHM1) (NUHM2). We display 1-dimensional likelihood functions for SUSY masses, and show that they may be significantly lighter in the pMSSM10 than in the CMSSM, NUHM1 and NUHM2. We discuss the discovery potential of future LHC runs, e+e- colliders and direct detection experiments.

We discuss phenomenological aspects of no-scale supergravity inflationary models motivated by compactified string models, in which the inflaton may be identified either as a K\"ahler modulus or an untwisted matter field, focusing on models that make predictions for the scalar spectral index $n_s$ and the tensor-to-scalar ratio $r$ that are similar to the Starobinsky model. We discuss possible patterns of soft supersymmetry breaking, exhibiting examples of the pure no-scale type $m_0 = B_0 = A_0 = 0$, of the CMSSM type with universal $A_0$ and $m_0 \ne 0$ at a high scale, and of the mSUGRA type with $A_0 = B_0 + m_0$ boundary conditions at the high input scale. These may be combined with a non-trivial gauge kinetic function that generates gaugino masses $m_{1/2} \ne 0$, or one may have a pure gravity mediation scenario where trilinear terms and gaugino masses are generated through anomalies. We also discuss inflaton decays and reheating, showing possible decay channels for the inflaton when it is either an untwisted matter field or a K\"ahler modulus. Reheating is very efficient if a matter field inflaton is directly coupled to MSSM fields, and both candidates lead to sufficient reheating in the presence of a non-trivial gauge kinetic function.

We consider the production and deposition on Earth of isotopes with half-lives in the range 10$^{5}$ to 10$^{8}$ years that might provide signatures of nearby stellar explosions, extending previous analyses of Core-Collapse Supernovae (CCSNe) to include Electron-Capture Supernovae (ECSNe), Super-Asymptotic Giant Branch (SAGBs) stars, Thermonuclear/Type Ia Supernovae (TNSNe), and Kilonovae/Neutron Star Mergers (KNe). We revisit previous estimates of the $^{60}$Fe and $^{26}$Al signatures, and extend these estimates to include $^{244}$Pu and $^{53}$Mn. We discuss interpretations of the $^{60}$Fe signals in terrestrial and lunar reservoirs in terms of a nearby stellar ejection ~2.2 Myr ago, showing that (i) the $^{60}$Fe yield rules out the TNSN and KN interpretations, (ii) the $^{60}$Fe signals highly constrain a SAGB interpretation but do not completely them rule out, (iii) are consistent with a CCSN origin, and (iv) are highly compatible with an ECSN interpretation. Future measurements could resolve the radioisotope deposition over time, and we use the Sedov blast wave solution to illustrate possible time-resolved profiles. Measuring such profiles would independently probe the blast properties including distance, and would provide additional constraints the nature of the explosion.

Among efforts to detect gravitational radiation, pulsar timing arrays are uniquely poised to detect "memory" signatures, permanent perturbations in spacetime from highly energetic astrophysical events such as mergers of supermassive black hole binaries. The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) observes dozens of the most stable millisecond pulsars using the Arecibo and Green Bank radio telescopes in an effort to study, among other things, gravitational wave memory. We herein present the results of a search for gravitational wave bursts with memory (BWMs) using the first five years of NANOGrav observations. We develop original methods for dramatically speeding up searches for BWM signals. In the directions of the sky where our sensitivity to BWMs is best, we would detect mergers of binaries with reduced masses of $10^9$ $M_\odot$ out to distances of 30 Mpc; such massive mergers in the Virgo cluster would be marginally detectable. We find no evidence for BWMs. However, with our non-detection, we set upper limits on the rate at which BWMs of various amplitudes could have occurred during the time spanned by our data--e.g., BWMs with amplitudes greater than $10^{-13}$ must occur at a rate less than 1.5 yr$^{-1}$.

Higgs inflation is attractive because it identifies the inflaton with the electroweak Higgs boson. In this work, we construct a new class of supersymmetric Higgs inflationary models in the no-scale supergravity with an SU(5) GUT group. Extending the no-scale Kahler potential and SU(5) GUT superpotential, we derive a generic potential for Higgs inflation that includes the quadratic monomial potential and a Starobinsky-type potential as special limits. This type of models can accommodate a wide range of the tensor-to-scalar ratio $r = O(10^{-3}-10^{-1})$, as well as a scalar spectral index $n_s \sim 0.96$.

We describe several new techniques which accelerate Bayesian searches for continuous gravitational-wave emission from supermassive black-hole binaries using pulsar timing arrays. These techniques mitigate the problematic increase of search-dimensionality with the size of the pulsar array which arises from having to include an extra parameter per pulsar as the array is expanded. This extra parameter corresponds to searching over the phase of the gravitational-wave as it propagates past each pulsar so that we can coherently include the pulsar-term in our search strategies. Our techniques make the analysis tractable with powerful evidence-evaluation packages like MultiNest. We find good agreement of our techniques with the parameter-estimation and Bayes factor evaluation performed with full signal templates, and conclude that these techniques make excellent first-cut tools for detection and characterisation of continuous gravitational-wave signals with pulsar timing arrays. Crucially, at low to moderate signal-to-noise ratios the factor by which the analysis is sped up can be > 100, permitting rigorous programs of systematic injection and recovery of signals to establish robust detection criteria within a Bayesian formalism.

Supermassive black hole binaries, cosmic strings, relic gravitational waves from inflation, and first order phase transitions in the early universe are expected to contribute to a stochastic background of gravitational waves in the 10^(-9) Hz-10^(-7) Hz frequency band. Pulsar timing arrays (PTAs) exploit the high precision timing of radio pulsars to detect signals at such frequencies. Here we present a time-domain implementation of the optimal cross-correlation statistic for stochastic background searches in PTA data. Due to the irregular sampling typical of PTA data as well as the use of a timing model to predict the times-of-arrival of radio pulses, time-domain methods are better suited for gravitational wave data analysis of such data. We present a derivation of the optimal cross-correlation statistic starting from the likelihood function, a method to produce simulated stochastic background signals, and a rigorous derivation of the scaling laws for the signal-to-noise ratio of the cross-correlation statistic in the two relevant PTA regimes: the weak signal limit where instrumental noise dominates over the gravitational wave signal at all frequencies, and a second regime where the gravitational wave signal dominates at the lowest frequencies.

Since the building-blocks of supersymmetric models include chiral superfields containing pairs of effective scalar fields, a two-field approach is particularly appropriate for models of inflation based on supergravity. In this paper, we generalize the two-field analysis of the inflationary power spectrum to supergravity models with arbitrary K\"ahler potential. We show how two-field effects in the context of no-scale supergravity can alter the model predictions for the scalar spectral index $n_s$ and the tensor-to-scalar ratio $r$, yielding results that interpolate between the Planck-friendly Starobinsky model and BICEP2-friendly predictions. In particular, we show that two-field effects in a chaotic no-scale inflation model with a quadratic potential are capable of reducing $r$ to very small values $\ll 0.1$. We also calculate the non-Gaussianity measure $f_{\rm NL}$, finding that is well below the current experimental sensitivity.

The magnitude of B-mode polarization in the cosmic microwave background as measured by BICEP2 favours models of chaotic inflation with a quadratic $m^2 \phi^2/2$ potential, whereas data from the Planck satellite favour a small value of the tensor-to-scalar perturbation ratio $r$ that is highly consistent with the Starobinsky $R + R^2$ model. Reality may lie somewhere between these two scenarios. In this paper we propose a minimal two-field no-scale supergravity model that interpolates between quadratic and Starobinsky-like inflation as limiting cases, while retaining the successful prediction $n_s \simeq 0.96$.

The arrival of TeV-energy photons from distant galaxies is expected to be affected by their QED interaction with intergalactic radiation fields through electron-positron pair production. In theories where high-energy photons violate Lorentz symmetry, the kinematics of the process $\gamma + \gamma\rightarrow e^+ + e^-$ is altered and the cross-section suppressed. Consequently, one would expect more of the highest-energy photons to arrive if QED is modified by Lorentz violation than if it is not. We estimate the sensitivity of Cherenkov Telescope Array (CTA) to changes in the $\gamma$-ray horizon of the Universe due to Lorentz violation, and find that it should be competitive with other leading constraints.

The Pulsar Arecibo L-band Feed Array (PALFA) Survey uses the ALFA 7-beam receiver to search both inner and outer Galactic sectors visible from Arecibo ($32^{\circ}\lesssim \ell \lesssim 77^{\circ}$ and $168^{\circ}\lesssim \ell \lesssim 214^{\circ}$) close to the Galactic plane ($|b|\lesssim5^{\circ}$) for pulsars. In this paper we detail a precursor survey of this region with PALFA, which observed a subset of the full region (slightly more restrictive in $\ell$ and $|b|\lesssim1^{\circ}$) and detected 45 pulsars. For both Galactic millisecond and normal pulsar populations, we compare the survey's detections with simulations to model these populations and, in particular, to estimate the number of observable pulsars in the Galaxy. We place 95\% confidence intervals of 82,000 to 143,000 on the number of detectable normal pulsars and 9,000 to 100,000 on the number of detectable millisecond pulsars in the Galactic disk. These are consistent with previous estimates. Given the most likely population size in each case (107,000 and 15,000 for normal and millisecond pulsars, respectively) we extend survey detection simulations to predict that, when complete, the full PALFA survey should have detected $1,000\substack{+330 \\ -230}$ normal pulsars and $30\substack{+200 \\ -20}$ millisecond pulsars. Identical estimation techniques predict that $490\substack{+160 \\ -115}$ normal pulsars and $12\substack{+70 \\ -5}$ millisecond pulsars would be detected by the beginning of 2014; at the time, the PALFA survey had detected 283 normal pulsars and 31 millisecond pulsars, respectively.

In a recent paper we showed how Starobinsky-like inflation could emerge from dilaton dynamics in brane cosmology scenarios based on string theory, in which our universe is represented as a three-brane and the effective potential acquires a constant term from a density of effectively point-like non-perturbative defects on the brane: "D-particles". Here we explore how quantum fluctuations of the ensemble of D-particles during the inflationary period may modify the effective inflationary potential due to the dilaton. We then discuss two specific ways in which an enhanced ratio of tensor to scalar perturbations may arise: via either a condensate of vector fields with a Born-Infeld action that may be due to such recoil fluctuations, or graviton production in the D-particle vacuum.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of $h_0\lesssim 3.8\times 10^{-14}$ at a frequency of 10 nHz. Furthermore, we place 95% \emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the $d_L \gtrsim 425$ Mpc for sources at a frequency of 10 nHz and chirp mass $10^{10}{\rm M}_{\odot}$. We find that for gravitational wave sources near our best timed pulsars in the sky, the sensitivity of the pulsar timing array is increased by a factor of $\sim$4 over the sky-averaged sensitivity. Finally we place limits on the coalescence rate of the most massive supermassive black hole binaries.

We discuss how Starobinsky-like inflation may emerge from dilaton dynamics in brane cosmology scenarios based on string theory, in which our universe is represented as a three-brane. The effective potential may acquire a constant term from a density of effectively point-like non-pertubative defects on the brane. Higher-genus corrections generate corrections to the effective potential that are exponentially damped at large field values, as in the Starobinsky model, but at a faster rate, leading to a smaller prediction for the tensor-to scalar perturbation ratio r. This may be compensated partially by logarithmic deformations on the world-sheet due to recoil of the defects due to scattering by string matter on the brane, which tend to enhance the tensor-to-scalar ratio.

The magnitude of primordial tensor perturbations reported by the BICEP2 experiment is consistent with simple models of chaotic inflation driven by a single scalar field with a power-law potential \propto \phi^n: n \simeq 2, in contrast to the WMAP and Planck results, which favored models resembling the Starobinsky R + R^2 model if running of the scalar spectral index could be neglected. While models of inflation with a quadratic potential may be constructed in simple N=1 supergravity, these constructions are more challenging in no-scale supergravity. We discuss here how quadratic inflation can be accommodated within supergravity, focussing primarily on the no-scale case. We also argue that the quadratic inflaton may be identified with the supersymmetric partner of a singlet (right-handed) neutrino, whose subsequent decay could have generated the baryon asymmetry via leptogenesis.

We propose a minimal model framework for physics below the Planck scale with the following features: (i) it is based on no-scale supergravity, as favoured in many string compactifications, (ii) it incorporates Starobinsky-like inflation, and hence is compatible with constraints from the Planck satellite, (iii) the inflaton may be identified with a singlet field in a see-saw model for neutrino masses, providing an efficient scenario for reheating and leptogenesis, (iv) supersymmetry breaking occurs with an arbitrary scale and a cosmological constant that vanishes before radiative corrections, (v) regions of the model parameter space are compatible with all LHC, Higgs and dark matter constraints.

We explore inflation via the effective potential of the minimal Wess-Zumino model, considering both the real and imaginary components of the complex field. Using transport techniques, we calculate the full allowed range of $n_s$, $r$ and $f_{\rm NL}$ for different choices of the single free parameter, $v$, and present the probability distribution of these signatures given a simple choice for the prior distribution of initial conditions. Our work provides a case study of multi-field inflation in a simple but realistic setting, with important lessons that are likely to apply more generally. For example, we find that there are initial conditions consistent with observations of $n_s$ and $r$ for values of $v$ that would be excluded if only evolutions in the real field direction were to be considered, and that these may yield enhanced values of $f_{\rm NL}$. Moreover, we find that initial conditions fixed at high energy density, where the potential is close to quartic in form, can still lead to evolutions in a concave region of the potential during the observable number of e-folds, as preferred by present data. The Wess-Zumino model therefore provides an illustration that multi-field dynamics must be taken into account when seeking to understand fully the phenomenology of such models of inflation.

Inflationary models based on a single scalar field $\phi$ with a quadratic potential $V = \frac{1}{2} m^2 \phi^2$ are disfavoured by the recent Planck constraints on the scalar index, $n_s$, and the tensor-to-scalar ratio for cosmological density perturbations, $r_T$. In this paper we study how such a quadratic inflationary model can be rescued by postulating additional fields with quadratic potentials, such as might occur in sneutrino models, which might serve as either curvatons or supplementary inflatons. Introducing a second scalar field reduces but does not remove the pressure on quadratic inflation, but we find a sample of three-field models that are highly compatible with the Planck data on $n_s$ and $r_T$. We exhibit a specific three-sneutrino example that is also compatible with the data on neutrino mass difference and mixing angles.

We discuss the emergence of an inflationary phase in supergravity with the super-Higgs effect due to dynamical spontaneous breaking of supersymmetry, in which the role of the inflaton is played by the gravitino condensate. Realistic models compatible with the Planck satellite CMB data are found in conformal supergravity scenarios with dynamical gravitino masses that are small compared to the Planck mass, as could be induced by a non-trivial vacuum expectation value of the dilaton superfield of appropriate magnitude.

We have assembled a tiled array (220 cm2) of fine pixel (0.6 mm) imaging CZT detectors for a balloon borne wide-field hard X-ray telescope, ProtoEXIST2. ProtoEXIST2 is a prototype experiment for a next generation hard X-ray imager MIRAX-HXI on board Lattes, a spacecraft from the Agencia Espacial Brasilieira. MIRAX will survey the 5 to 200 keV sky of Galactic bulge, adjoining southern Galactic plane and the extragalactic sky with 6' angular resolution. This survey will open a vast discovery space in timing studies of accretion neutron stars and black holes. The ProtoEXIST2 CZT detector plane consists of 64 of 5 mm thick 2 cm x 2 cm CZT crystals tiled with a minimal gap. MIRAX will consist of 4 such detector planes, each of which will be imaged with its own coded-aperture mask. We present the packaging architecture and assembly procedure of the ProtoEXIST2 detector. On 2012, Oct 10, we conducted a successful high altitude balloon experiment of the ProtoEXIST1 and 2 telescopes, which demonstrates their technology readiness for space application. During the flight both telescopes performed as well as on the ground. We report the results of ground calibration and the initial results for the detector performance in the balloon flight.

We present a model for cosmological inflation based on a no-scale supergravity sector with an SU(2,1)/U(1) Kahler potential, a single modulus T and an inflaton superfield Phi described by a Wess-Zumino model with superpotential parameters (mu, lambda). This model yields a scalar spectral index n_s and a tensor-to-scalar ratio r that are compatible with the Planck measurements for values of lambda simeq mu/3M_P. For the specific choice lambda = mu/3M_P, the model is a no-scale supergravity realization of the R+R^2 Starobinsky model.

Models of cosmological inflation resembling the Starobinsky R + R^2 model emerge naturally among the effective potentials derived from no-scale SU(N,1)/SU(N) x U(1) supergravity when N > 1. We display several examples in the SU(2,1)/SU(2) x U(1) case, in which the inflaton may be identified with either a modulus field or a matter field. We discuss how the modulus field may be stabilized in models in which a matter field plays the role of the inflaton. We also discuss models that generalize the Starobinsky model but display different relations between the tilt in the spectrum of scalar density perturbations, n_s, the tensor-to-scalar ratio, r, and the number of e-folds, N_*. Finally, we discuss how such models can be probed by present and future CMB experiments.

We studied previously the impact on light-element abundances of gravitinos decaying during or after Big-Bang nucleosynthesis (BBN). We found regions of the gravitino mass m_{3/2} and abundance zeta_{3/2} plane where its decays could reconcile the calculated abundance of Li7 with observation without perturbing the other light-element abundances unacceptably. Here we revisit this issue in light of LHC measurements of the Higgs mass and constraints on supersymmetric model parameters, as well as updates in the astrophysical measurements of light-element abundances. In addition to the constrained minimal supersymmetric extension of the Standard Model with universal soft supersymmetry-breaking masses at the GUT scale (the CMSSM) studied previously, we also study models with universality imposed below the GUT scale and models with non-universal Higgs masses (NUHM1). We calculate the total likelihood function for the light-element abundances, taking into account the observational uncertainties. We find that gravitino decays provide a robust solution to the cosmological Li7 problem along strips in the (m_{3/2}, zeta_{3/2}) plane along which the abundances of deuterium, He4 and Li7 may be fit with chi^2_min < 3, compared with chi^2 ~ 34 if the effects of gravitino decays are unimportant. The minimum of the likelihood function is reduced to chi^2 < 2 when the uncertainty on D/H is relaxed and < 1 when the lithium abundance is taken from globular cluster data.

We derive scaling laws for the signal-to-noise ratio of the optimal cross-correlation statistic, and show that the large power-law increase of the signal-to-noise ratio as a function of the the observation time $T$ that is usually assumed holds only at early times. After enough time has elapsed, pulsar timing arrays enter a new regime where the signal to noise only scales as $\sqrt{T}$. In addition, in this regime the quality of the pulsar timing data and the cadence become relatively un-important. This occurs because the lowest frequencies of the pulsar timing residuals become gravitational-wave dominated. Pulsar timing arrays enter this regime more quickly than one might naively suspect. For T=10 yr observations and typical stochastic background amplitudes, pulsars with residual RMSs of less than about $1\,\mu$s are already in that regime. The best strategy to increase the detectability of the background in this regime is to increase the number of pulsars in the array. We also perform realistic simulations of the NANOGrav pulsar timing array, which through an aggressive pulsar survey campaign adds new millisecond pulsars regularly to its array, and show that a detection is possible within a decade, and could occur as early as 2016.

The direct detection of Gravitational Waves (GWs) by Pulsar Timing Arrays (PTAs) is very likely within the next decade. While the stochastic GW background is a promising candidate for detection it is also possible that single resolvable sources may be detectable as well. In this work we will focus on the detection and characterization of single GW sources from supermassive black hole binaries (SMBHBs). We introduce a fully Bayesian data analysis pipeline that is meant to carry out a search, characterization, and evaluation phase. This will allow us to rapidly locate the global maxima in parameter space, map out the posterior, and finally weigh the evidence of a GW detection through a Bayes Factor. Here we will make use of an adaptive metropolis (AM) algorithm and parallel tempering. We test this algorithm on realistic simulated data that are representative of modern PTAs.

We discuss cosmological inflation in the minimal Wess-Zumino model with a single massive chiral supermultiplet. With suitable parameters and assuming a plausible initial condition at the start of the inflationary epoch, the model can yield scalar perturbations in the Cosmic Microwave Background (CMB) of the correct strength with a spectral index n_s ~ 0.96 and a tensor-to-scalar perturbation ratio r < 0.1, consistent with the Planck CMB data. We also discuss the possibility of topological inflation within the Wess-Zumino model, and the possibility of combining it with a seesaw model for neutrino masses. This would violate R-parity, but at such a low rate that the lightest supersymmetric particle would have a lifetime long enough to constitute the astrophysical cold dark matter.

We discuss the stability and masses of topological solitons in QCD and strongly-interacting models of electroweak symmetry breaking with arbitrary combinations of two inequivalent Lagrangian terms of fourth order in the field spatial derivatives. We find stable solitons for only a restricted range of the ratio of these combinations, in agreement with previous results, and we calculate the corresponding soliton masses. In QCD, the experimental constraints on the fourth-order terms force the soliton to resemble the original Skyrmion solution. However, this is not necessarily the case in strongly-interacting models of electroweak symmetry breaking, in which a non-Skyrmion-like soliton is also possible. This possibility will be constrained by future LHC measurements and dark matter experiments. Current upper bounds on the electroweak soliton mass range between 18 and 59 TeV, which would be reduced to 4.6 to 8.1 TeV with the likely sensitivity of LHC data to the fourth-order electroweak Lagrangian parameters.

Direct detection of gravitational waves by pulsar timing arrays will become feasible over the next few years. In the low frequency regime ($10^{-7}$ Hz -- $10^{-9}$ Hz), we expect that a superposition of gravitational waves from many sources will manifest itself as an isotropic stochastic gravitational wave background. Currently, a number of techniques exist to detect such a signal; however, many detection methods are computationally challenging. Here we introduce an approximation to the full likelihood function for a pulsar timing array that results in computational savings proportional to the square of the number of pulsars in the array. Through a series of simulations we show that the approximate likelihood function reproduces results obtained from the full likelihood function. We further show, both analytically and through simulations, that, on average, this approximate likelihood function gives unbiased parameter estimates for astrophysically realistic stochastic background amplitudes.

The 2012 International Pulsar Timing Array (IPTA) Mock Data Challenge (MDC) is designed to test current Gravitational Wave (GW) detection algorithms. Here we will briefly outline two detection algorithms for a stochastic background of gravitational waves, namely, a first-order likelihood method and an optimal statistic method and present our results from the closed MDC data sets.

We present an experimental platform for measuring hydrogen Balmer emission and absorption line profiles for plasmas with white dwarf (WD) photospheric conditions (T_e ~ 1 eV, n_e ~ 10^17 cm^-3). These profiles will be used to benchmark WD atmosphere models, which, used with the spectroscopic method, are responsible for determining fundamental parameters (e.g., effective temperature, mass) for tens of thousands of WDs. Our experiment, performed at the Z Pulsed Power Facility at Sandia National Laboratories, uses the large amount of x-rays generated from a z-pinch dynamic hohlraum to drive plasma formation in a gas cell. The platform is unique compared to past hydrogen line profile experiments in that the plasma is radiation-driven. This decouples the heating source from the plasma to be studied in the sense that the radiation temperature causing the photoionization is independent of the initial conditions of the gas. For the first time we measure hydrogen Balmer lines in absorption at these conditions in the laboratory for the purpose of benchmarking Stark-broadened line shapes. The platform can be used to study other plasma species and to explore non-LTE, time-dependent collisional-radiative atomic kinetics.

Astrophysics experiments by Falcon et al. to create white dwarf photospheres in the laboratory are currently underway. The experimental platform measures Balmer line profiles of a radiation-driven, pure hydrogen plasma in emission and in absorption for conditions at T_e ~ 1 eV, n_e ~ 10^17 cm^-3. These will be used to compare and test line broadening theories used in white dwarf atmosphere models. The flexibility of the platform allows us to expand the direction of our experiments using other compositions. We discuss future prospects such as exploring helium plasmas and carbon/oxygen plasmas relevant to the photospheres of DBs and hot DQs, respectively.

We consider the effects of metastable charged sparticles on Big-Bang Nucleosynthesis (BBN), including bound-state reaction rates and chemical effects. We make a new analysis of the bound states of negatively-charged massive particles with the light nuclei most prominent in BBN, and present a new code to track their abundances, paying particular attention to that of Li7. Assuming, as an example, that the gravitino is the lightest supersymmetric particle (LSP), and that the lighter stau slepton, stau_1, is the metastable next-to-lightest sparticle within the constrained minimal supersymmetric extension of the Standard Model (CMSSM), we analyze the possible effects on the standard BBN abundances of stau_1 bound states and decays for representative values of the gravitino mass. Taking into account the constraint on the CMSSM parameter space imposed by the discovery of the Higgs boson at the LHC, we delineate regions in which the fit to the measured light-element abundances is as good as in standard BBN. We also identify regions of the CMSSM parameter space in which the bound state properties, chemistry and decays of metastable charged sparticles can solve the cosmological Li7 problem.

Global frequentist fits to the CMSSM and NUHM1 using the MasterCode framework predicted m_h \simeq 119 GeV in fits incorporating the g_mu-2 constraint and \simeq 126 GeV without it. Recent results by ATLAS and CMS could be compatible with a Standard Model-like Higgs boson around m_h \simeq 125 GeV. We use the previous MasterCode analysis to calculate the likelihood for a measurement of any nominal Higgs mass within the range of 115 to 130 GeV. Assuming a Higgs mass measurement at m_h \simeq 125 GeV, we display updated global likelihood contours in the (m_0, m_{1/2}) and other parameter planes of the CMSSM and NUHM1, and present updated likelihood functions for m_gluino, m_squark, B to mu mu, and the spin-independent dark matter cross section \sigma^si. The implications of dropping g_mu-2 from the fits are also discussed. We furthermore comment on a hypothetical measurement of m_h \simeq 119 GeV.

Supermassive black hole binaries (SMBHBs) are expected to emit continuous gravitational waves in the pulsar timing array (PTA) frequency band ($10^{-9}$--$10^{-7}$ Hz). The development of data analysis techniques aimed at efficient detection and characterization of these signals is critical to the gravitational wave detection effort. In this paper we leverage methods developed for LIGO continuous wave gravitational searches, and explore the use of the $\mathcal{F}$-statistic for such searches in pulsar timing data. Babak & Sesana 2012 have already used this approach in the context of PTAs to show that one can resolve multiple SMBHB sources in the sky. Our work improves on several aspects of prior continuous wave search methods developed for PTA data analysis. The algorithm is implemented fully in the time domain, which naturally deals with the irregular sampling typical of PTA data and avoids spectral leakage problems associated with frequency domain methods. We take into account the fitting of the timing model, and have generalized our approach to deal with both correlated and uncorrelated colored noise sources. We also develop an incoherent detection statistic that maximizes over all pulsar dependent contributions to the likelihood. To test the effectiveness and sensitivity of our detection statistics, we perform a number of monte-carlo simulations. We produce sensitivity curves for PTAs of various configurations, and outline an implementation of a fully functional data analysis pipeline. Finally, we present a derivation of the likelihood maximized over the gravitational wave phases at the pulsar locations, which results in a vast reduction of the search parameter space.

We discuss the prospects for improved upper limits on neutrino masses that may be provided by a core-collapse supernova explosion in our galaxy, if it exhibits time variations in the neutrino emissions on the scale of a few milliseconds as suggested by recent two-dimensional simulations. Analyzing simulations of such neutrino emissions using the wavelet technique adopted in [1], we find that an upper limit m_nu ~ 0.14 eV could be established at the 95% confidence level if the time variations in emissions were to be preserved during neutrino propagation to the Earth.

We investigate the possibility that the anomalous acceleration of the Pioneer 10 and 11 spacecraft is due to the recoil force associated with an anisotropic emission of thermal radiation off the vehicles. To this end, relying on the project and spacecraft design documentation, we constructed a comprehensive finite-element thermal model of the two spacecraft. Then, we numerically solve thermal conduction and radiation equations using the actual flight telemetry as boundary conditions. We use the results of this model to evaluate the effect of the thermal recoil force on the Pioneer 10 spacecraft at various heliocentric distances. We found that the magnitude, temporal behavior, and direction of the resulting thermal acceleration are all similar to the properties of the observed anomaly. As a novel element of our investigation, we develop a parameterized model for the thermal recoil force and estimate the coefficients of this model independently from navigational Doppler data. We find no statistically significant difference between the two estimates and conclude that once the thermal recoil force is properly accounted for, no anomalous acceleration remains.

Taking into account the available accelerator and astrophysical constraints, the mass of the lightest neutral Higgs boson h in the minimal supersymmetric extension of the Standard Model with universal soft supersymmetry-breaking masses (CMSSM) has been estimated to lie between 114 and ~ 130 GeV. Recent data from ATLAS and CMS hint that m_h ~ 125 GeV, though m_h ~ 119 GeV may still be a possibility. Here we study the consequences for the parameters of the CMSSM and direct dark matter detection if the Higgs hint is confirmed, focusing on the strips in the (m_1/2, m_0) planes for different tan beta and A_0 where the relic density of the lightest neutralino chi falls within the range of the cosmological cold dark matter density allowed by WMAP and other experiments. We find that if m_h ~ 125 GeV focus-point strips would be disfavoured, as would the low-tan beta stau-chi and stop -chi coannihilation strips, whereas the stau-chi coannihilation strip at large tan beta and A_0 > 0 would be favoured, together with its extension to a funnel where rapid annihilation via direct-channel H/A poles dominates. On the other hand, if m_h ~ 119 GeV more options would be open. We give parametrizations of WMAP strips with large tan beta and fixed A_0/m_0 > 0 that include portions compatible with m_h = 125 GeV, and present predictions for spin-independent elastic dark matter scattering along these strips. These are generally low for models compatible with m_h = 125 GeV, whereas the XENON100 experiment already excludes some portions of strips where m_h is smaller.

Gravitational Waves (GWs) are tiny ripples in the fabric of space-time predicted by Einstein's General Relativity. Pulsar timing arrays (PTAs) are well poised to detect low frequency ($10^{-9}$ -- $10^{-7}$ Hz) GWs in the near future. There has been a significant amount of research into the detection of a stochastic background of GWs from supermassive black hole binaries (SMBHBs). Recent work has shown that single continuous sources standing out above the background may be detectable by PTAs operating at a sensitivity sufficient to detect the stochastic background. The most likely sources of continuous GWs in the pulsar timing frequency band are extremely massive and/or nearby SMBHBs. In this paper we present detection strategies including various forms of matched filtering and power spectral summing. We determine the efficacy and computational cost of such strategies. It is shown that it is computationally infeasible to use an optimal matched filter including the poorly constrained pulsar distances with a grid based method. We show that an Earth-term-matched filter constructed using only the correlated signal terms is both computationally viable and highly sensitive to GW signals. This technique is only a factor of two less sensitive than the computationally unrealizable optimal matched filter and a factor of two more sensitive than a power spectral summing technique. We further show that a pairwise matched filter, taking the pulsar distances into account is comparable to the optimal matched filter for the single template case and comparable to the Earth-term-matched filter for many search templates. Finally, using simulated data optimal quality, we place a theoretical minimum detectable strain amplitude of $h>2\times 10^{-15}$ from continuous GWs at frequencies on the order $\sim1/T_{\rm obs}$.

Supernova explosions provide the most sensitive probes of neutrino propagation, such as the possibility that neutrino velocities might be affected by the foamy structure of space-time thought to be generated by quantum-gravitational (QG) effects. Recent two-dimensional simulations of the neutrino emissions from core-collapse supernovae suggest that they might exhibit variations in time on the scale of a few milliseconds. We analyze simulations of such neutrino emissions using a wavelet technique, and consider the limits that might be set on a linear or quadratic violation of Lorentz invariance in the group velocities of neutrinos of different energies, v/c = [1 \pm (E/M_{nuLV1})] or [1 \pm (E/M_{\nuLV2})^2], if variations on such short time scales were to be observed, where the mass scales M_{nuLVi} might appear in models of quantum gravity. We find prospective sensitivities to M_{nuLV1} ~ 2 X 10^{13} GeV and M_{nuLV2} ~ 10^6 GeV at the 95% confidence level, up to two orders of magnitude beyond estimates made using previous one-dimensional simulations of core-collapse supernovae. We also analyze the prospective sensitivities to scenarios in which the propagation times of neutrinos of fixed energies are subject to stochastic fluctuations.

We present an analysis of high-precision pulsar timing data taken as part of the North American Nanohertz Observatory for Gravitational waves (NANOGrav) project. We have observed 17 pulsars for a span of roughly five years using the Green Bank and Arecibo radio telescopes. We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Sub-microsecond timing residuals are obtained in nearly all cases, and the best root-mean-square timing residuals in this set are ~30-50 ns. We present methods for analyzing post-fit timing residuals for the presence of a gravitational wave signal with a specified spectral shape. These optimally take into account the timing fluctuation power removed by the model fit, and can be applied to either data from a single pulsar, or to a set of pulsars to detect a correlated signal. We apply these methods to our dataset to set an upper limit on the strength of the nHz-frequency stochastic supermassive black hole gravitational wave background of h_c (1 yr^-1) < 7x10^-15 (95%). This result is dominated by the timing of the two best pulsars in the set, PSRs J1713+0747 and J1909-3744.

In many strongly-interacting models of electroweak symmetry breaking the lowest-lying observable particle is a pseudo-Goldstone boson of approximate scale symmetry, the pseudo-dilaton. Its interactions with Standard Model particles can be described using a low-energy effective nonlinear chiral Lagrangian supplemented by terms that restore approximate scale symmetry, yielding couplings of the pseudo-dilaton that differ from those of a Standard Model Higgs boson by fixed factors. We review the experimental constraints on such a pseudo-dilaton in light of new data from the LHC and elsewhere. The effective nonlinear chiral Lagrangian has Skyrmion solutions that may be identified with the `electroweak baryons' of the underlying strongly-interacting theory, whose nature may be revealed by the properties of the Skyrmions. We discuss the finite-temperature electroweak phase transition in the low-energy effective theory, finding that the possibility of a first-order electroweak phase transition is resurrected. We discuss the evolution of the Universe during this transition and derive an order-of-magnitude lower limit on the abundance of electroweak baryons in the absence of a cosmological asymmetry, which suggests that such an asymmetry would be necessary if the electroweak baryons are to provide the cosmological density of dark matter. We revisit estimates of the corresponding spin-independent dark matter scattering cross section, with a view to direct detection experiments.

Lorentz invariance is such an important principle of fundamental physics that it should constantly be subjected to experimental scrutiny as well as theoretical questioning. Distant astrophysical sources of energetic photons with rapid time variations, such as active galactic nuclei (AGNs) and gamma-ray bursters (GRBs), provide ideal experimental opportunities for testing Lorentz invariance. The Cherenkov Telescope Array (CTA) is an excellent experimental tool for making such tests with sensitivities exceeding those possible using other detectors.

We study the production of gamma rays via LSP annihilations in the core of the Galaxy as a possible experimental signature of the constrained minimal supersymmetric extension of the Standard Model (CMSSM), in which supersymmetry-breaking parameters are assumed to be universal at the GUT scale, assuming also that the LSP is the lightest neutralino chi. The part of the CMSSM parameter space that is compatible with the measured astrophysical density of cold dark matter is known to include a stau_1 - chi coannihilation strip, a focus-point strip where chi has an enhanced Higgsino component, and a funnel at large tanb where the annihilation rate is enhanced by the poles of nearby heavy MSSM Higgs bosons, A/H. We calculate the total annihilation rates, the fractions of annihilations into different Standard Model final states and the resulting fluxes of gamma rays for CMSSM scenarios along these strips. We observe that typical annihilation rates are much smaller in the coannihilation strip for tanb = 10 than along the focus-point strip or for tanb = 55, and that the annihilation branching ratios differ greatly between the different dark matter strips. Whereas the current Fermi-LAT data are not sensitive to any of the CMSSM scenarios studied, and the calculated gamma-ray fluxes are probably unobservably low along the coannihilation strip for tanb = 10, we find that substantial portions of the focus-point strips and rapid-annihilation funnel regions could be pressured by several more years of Fermi-LAT data, if understanding of the astrophysical background and/or systematic uncertainties can be improved in parallel.

Residual late-time dark matter particle annihilations during and after Big Bang Nucleosynthesis (BBN) may alter the predicted cosmological abundances of the light elements. Within the constrained minimal supersymmetric extension of the Standard Model (the CMSSM) with a neutralino LSP, we find negligible effects on the abundances of Deuterium, He3, He4 and Li7 predicted by homogeneous BBN, but potentially a large enhancement in the predicted abundance of Li6. This enhancement may be as much as two orders of magnitude in the focus-point WMAP strip and in the coannihilation and funnel regions for large tan beta for small m_{1/2}, and the effect is still significant at large m_{1/2}. However, the potential Li6 enhancement is negligible in the part of the coannihilation strip for tan beta = 10 that survives the latest LHC constraints. A similar enhancement of the \li6 abundance may also be found in a model with common, non-universal Higgs masses (the NUHM1).

Gravitational waves are predicted by Einstein's theory of general relativity as well as other theories of gravity. The rotational stability of the fastest pulsars means that timing of an array of these objects can be used to detect and investigate gravitational waves. Simultaneously, however, pulsar timing is used to estimate spin period, period derivative, astrometric, and binary parameters. Here we calculate the effects that a stochastic background of gravitational waves has on pulsar timing parameters through the use of simulations and data from the millisecond pulsars PSR J0437--4715 and PSR J1713+0747. We show that the reported timing uncertainties become underestimated with increasing background amplitude by up to a factor of $\sim10$ for a stochastic gravitational-wave background amplitude of $A=5\times 10^{-15}$, where $A$ is the amplitude of the characteristic strain spectrum at one-year gravitational wave periods. We find evidence for prominent low-frequency spectral leakage in simulated data sets including a stochastic gravitational-wave background. We use these simulations along with independent Very Long Baseline Interferometry (VLBI) measurements of parallax to set a 2--sigma upper limit of $A\le9.1\times 10^{-14}$. We find that different supermassive black hole assembly scenarios do not have a significant effect on the calculated upper limits. We also test the effects that ultralow--frequency (10$^{-12}$--10$^{-9}$ Hz) gravitational waves have on binary pulsar parameter measurements and find that the corruption of these parameters is less than those due to $10^{-9}$--$10^{-7}$ Hz gravitational waves.

We explore the phenomenology of the minimal supersymmetric flipped SU(5) GUT model (CFSU(5)), whose soft supersymmetry-breaking (SSB) mass parameters are constrained to be universal at some input scale, $M_{in}$, above the GUT scale, $M_{GUT}$. We analyze the parameter space of CFSU(5) assuming that the lightest supersymmetric particle (LSP) provides the cosmological cold dark matter, paying careful attention to the matching of parameters at the GUT scale. We first display some specific examples of the evolutions of the SSB parameters that exhibit some generic features. Specifically, we note that the relationship between the masses of the lightest neutralino and the lighter stau is sensitive to $M_{in}$, as is the relationship between the neutralino mass and the masses of the heavier Higgs bosons. For these reasons, prominent features in generic $(m_{1/2}, m_0)$ planes such as coannihilation strips and rapid-annihilation funnels are also sensitive to $M_{in}$, as we illustrate for several cases with tan(beta)=10 and 55. However, these features do not necessarily disappear at large $M_{in}$, unlike the case in the minimal conventional SU(5) GUT. Our results are relatively insensitive to neutrino masses.

The prospects for detecting a candidate supersymmetric dark matter particle at the LHC are reviewed, and compared with the prospects for direct and indirect searches for astrophysical dark matter. The discussion is based on a frequentist analysis of the preferred regions of the Minimal supersymmetric extension of the Standard Model with universal soft supersymmetry breaking (the CMSSM). LHC searches may have good chances to observe supersymmetry in the near future - and so may direct searches for astrophysical dark matter particles, whereas indirect searches may require greater sensitivity, at least within the CMSSM.

We make frequentist analyses of the CMSSM, NUHM1, VCMSSM and mSUGRA parameter spaces taking into account all the public results of searches for supersymmetry using data from the 2010 LHC run and the Xenon100 direct search for dark matter scattering. The LHC data set includes ATLAS and CMS searches for jets + ETslash events (with or without leptons) and for the heavier MSSM Higgs bosons, and the upper limit on bs to mu mu including data from LHCb as well as CDF and D0. The absences of signals in the LHC data favour somewhat heavier mass spectra than in our previous analyses of the CMSSM, NUHM1 and VCMSSM, and somewhat smaller dark matter scattering cross sections, all close to or within the pre-LHC 68% CL ranges, but do not impact significantly the favoured regions of the mSUGRA parameter space. We also discuss the impact of the Xenon100 constraint on spin-independent dark matter scattering, stressing the importance of taking into account the uncertainty in the pi-nucleon sigma term, that affects the spin-independent scattering matrix element, and we make predictions for spin-dependent dark matter scattering. Finally, we discuss briefly the potential impact of the updated predictions for sparticle masses in the CMSSM, NUHM1, VCMSSM and mSUGRA on future e+ e- colliders.

We extend our previous studies of the neutrino fluxes expected from neutralino LSP annihilations inside the Sun to include variants of the minimal supersymmetric extension of the Standard Model (MSSM) with squark, slepton and gaugino masses constrained to be universal at the GUT scale, but allowing one or two non-universal supersymmetry-breaking parameters contributing to the Higgs masses (NUHM1,2). As in the constrained MSSM (CMSSM) with universal Higgs masses, there are large regions of the NUHM parameter space where the LSP density inside the Sun is not in equilibrium, so that the annihilation rate may be far below the capture rate, and there are also large regions where the capture rate is not dominated by spin-dependent LSP-proton scattering. The spectra possible in the NUHM are qualitatively similar to those in the CMSSM. We calculate neutrino-induced muon fluxes above a threshold energy of 10 GeV, appropriate for the IceCube/DeepCore detector, for points where the NUHM yields the correct cosmological relic density for representative choices of the NUHM parameters. We find that the IceCube/DeepCore detector can probe regions of the NUHM parameter space in addition to analogues of the focus-point strip and the tip of the coannihilation strip familiar from the CMSSM. These include regions with enhanced Higgsino-gaugino mixing in the LSP composition, that occurs where neutralino mass eigenstates cross over. On the other hand, rapid-annihilation funnel regions in general yield neutrino fluxes that are unobservably small.

We make a frequentist analysis of the parameter space of minimal supergravity (mSUGRA), in which, as well as the gaugino and scalar soft supersymmetry-breaking parameters being universal, there is a specific relation between the trilinear, bilinear and scalar supersymmetry-breaking parameters, A_0 = B_0 + m_0, and the gravitino mass is fixed by m_{3/2} = m_0. We also consider a more general model, in which the gravitino mass constraint is relaxed (the VCMSSM). We combine in the global likelihood function the experimental constraints from low-energy electroweak precision data, the anomalous magnetic moment of the muon, the lightest Higgs boson mass M_h, B physics and the astrophysical cold dark matter density, assuming that the lightest supersymmetric particle (LSP) is a neutralino. In the VCMSSM, we find a preference for values of m_{1/2} and m_0 similar to those found previously in frequentist analyses of the constrained MSSM (CMSSM) and a model with common non-universal Higgs masses (NUHM1). On the other hand, in mSUGRA we find two preferred regions: one with larger values of both m_{1/2} and m_0 than in the VCMSSM, and one with large m_0 but small m_{1/2}. We compare the probabilities of the frequentist fits in mSUGRA, the VCMSSM, the CMSSM and the NUHM1: the probability that mSUGRA is consistent with the present data is significantly less than in the other models. We also discuss the mSUGRA and VCMSSM predictions for sparticle masses and other observables, identifying potential signatures at the LHC and elsewhere.

Proceedings of the 13th International Conference on Elastic and Diffractive Scattering (Blois Workshop) - Moving Forward into the LHC Era

The prospects for detecting a candidate supersymmetric dark matter particle at the LHC are reviewed, and compared with the prospects for direct and indirect searches for astrophysical dark matter, on the basis of a frequentist analysis of the preferred regions of the Minimal supersymmetric extension of the Standard Model with universal soft supersymmetry breaking (the CMSSM) and a model with equal but non-universal supersymmetry-breaking contributions to the Higgs masses (the NUHM1). LHC searches may have good chances to observe supersymmetry in the near future - and so may direct searches for astrophysical dark matter particles.

We consider the effects of uncertainties in nuclear reaction rates on the cosmological constraints on the decays of unstable particles during or after Big-Bang nucleosynthesis (BBN). We identify the nuclear reactions due to non-thermal hadrons that are the most important in perturbing standard BBN, then quantify the uncertainties in these reactions and in the resulting light-element abundances. These results also indicate the key nuclear processes for which improved cross section data would allow different light-element abundances to be determined more accurately, thereby making possible more precise probes of BBN and evaluations of the cosmological constraints on unstable particles. Applying this analysis to models with unstable gravitinos decaying into neutralinos, we calculate the likelihood function for the light-element abundances measured currently, taking into account the current experimental errors in the determinations of the relevant nuclear reaction rates. We find a region of the gravitino mass and abundance in which the abundances of deuterium, He4 and Li7 may be fit with chi^2 = 5.5, compared with chi^2 = 31.7 if the effects of gravitino decays are unimportant. The best-fit solution is improved to chi^2 ~ 2.0 when the lithium abundance is taken from globular cluster data. Some such re-evaluation of the observed light-element abundances and/or nuclear reaction rates would be needed if this region of gravitino parameters is to provide a complete solution to the cosmological Li7 problem.

We study models in which soft supersymmetry-breaking parameters of the MSSM become universal at some unification scale, $M_{in}$, above the GUT scale, $\mgut$. We assume that the scalar masses and gaugino masses have common values, $m_0$ and $m_{1/2}$ respectively, at $M_{in}$. We use the renormalization-group equations of the minimal supersymmetric SU(5) GUT to evaluate their evolutions down to $\mgut$, studying their dependences on the unknown parameters of the SU(5) superpotential. After displaying some generic examples of the evolutions of the soft supersymmetry-breaking parameters, we discuss the effects on physical sparticle masses in some specific examples. We note, for example, that near-degeneracy between the lightest neutralino and the lighter stau is progressively disfavoured as $M_{in}$ increases. This has the consequence, as we show in $(m_{1/2}, m_0)$ planes for several different values of $\tan \beta$, that the stau coannihilation region shrinks as $M_{in}$ increases, and we delineate the regions of the $(M_{in}, \tan \beta)$ plane where it is absent altogether. Moreover, as $M_{in}$ increases, the focus-point region recedes to larger values of $m_0$ for any fixed $\tan \beta$ and $m_{1/2}$. We conclude that the regions of the $(m_{1/2}, m_0)$ plane that are commonly favoured in phenomenological analyses tend to disappear at large $M_{in}$.

We revisit the constraints that the non-observation of ultra-high-energy photons due to the GZK cutoff can impose on models of Lorentz violation in photon propagation, following recent work by Maccione, Liberati and Sigl [arXiv:1003.5468] that carries further an earlier analysis by the present authors (Phys. Rev. D 63,12402 (2001), [hep-th/0012216]). We argue that the GZK cutoff constraint is naturally evaded in the D-brane model of space-time foam presented recently by the present authors (Phys. Lett. B 665, 412 (2008) [arXiv:0804.3566]), in which Lorentz-violating effects on photon propagation are independent of possible effects during interactions. We also note a novel absorption mechanism that could provide a GZK-like cutoff for photons in low-scale string models.

We evaluate the neutrino fluxes to be expected from neutralino LSP annihilations inside the Sun, within the minimal supersymmetric extension of the Standard Model with supersymmetry-breaking scalar and gaugino masses constrained to be universal at the GUT scale (the CMSSM). We find that there are large regions of typical CMSSM $(m_{1/2}, m_0)$ planes where the LSP density inside the Sun is not in equilibrium, so that the annihilation rate may be far below the capture rate. We show that neutrino fluxes are dependent on the solar model at the 20% level, and adopt the AGSS09 model of Serenelli et al. for our detailed studies. We find that there are large regions of the CMSSM $(m_{1/2}, m_0)$ planes where the capture rate is not dominated by spin-dependent LSP-proton scattering, e.g., at large $m_{1/2}$ along the CMSSM coannihilation strip. We calculate neutrino fluxes above various threshold energies for points along the coannihilation/rapid-annihilation and focus-point strips where the CMSSM yields the correct cosmological relic density for tan(beta) = 10 and 55 for $\mu$ > 0, exploring their sensitivities to uncertainties in the spin-dependent and -independent scattering matrix elements. We also present detailed neutrino spectra for four benchmark models that illustrate generic possibilities within the CMSSM. Scanning the cosmologically-favored parts of the parameter space of the CMSSM, we find that the IceCube/DeepCore detector can probe at best only parts of this parameter space, notably the focus-point region and possibly also at the low-mass tip of the coannihilation strip.

After reviewing the theoretical, phenomenological and experimental motivations for supersymmetric extensions of the Standard Model, we recall that supersymmetric relics from the Big Bang are expected in models that conserve R parity. We then discuss possible supersymmetric dark matter candidates, focusing on the lightest neutralino and the gravitino. In the latter case, the next-to-lightest supersymmetric particle is expected to be long-lived, and possible candidates include spartners of the tau lepton, top quark and neutrino. We then discuss the roles of the renormalization-group equations and electroweak symmetry breaking in delimiting the supersymmetric parameter space. We discuss in particular the constrained minimal extension of the Standard Model (CMSSM), in which the supersymmetry-breaking parameters are assumed to be universal at the grand unification scale, presenting predictions from a frequentist analysis of its parameter space. We also discuss astrophysical and cosmological constraints on gravitino dark matter models, as well as the parameter space of minimal supergravity (mSUGRA) models in which there are extra relations between the trilinear and bilinear supersymmetry-breaking parameters, and between the gravitino and scalar masses. Finally, we discuss models with non-universal supersymmetry-breaking contributions to Higgs masses, and models in which the supersymmetry-breaking parameters are universal at some scale below that of grand unification. http://cambridge.org/us/catalogue/catalogue.asp?isbn=9780521763684

In a D-brane model of space-time foam, there are contributions to the dark energy that depend on the D-brane velocities and on the density of D-particle defects. The latter may also reduce the speeds of photons linearly with their energies, establishing a phenomenological connection with astrophysical probes of the universality of the velocity of light. Specifically, the cosmological dark energy density measured at the present epoch may be linked to the apparent retardation of energetic photons propagating from nearby AGNs. However, this nascent field of `D-foam phenomenology' may be complicated by a dependence of the D-particle density on the cosmological epoch. A reduced density of D-particles at redshifts z ~ 1 - a `D-void' - would increase the dark energy while suppressing the vacuum refractive index, and thereby might reconcile the AGN measurements with the relatively small retardation seen for the energetic photons propagating from GRB 090510, as measured by the Fermi satellite.

The decays of massive gravitinos into neutralino dark matter particles and Standard Model secondaries during or after Big-Bang nucleosynthesis (BBN) may alter the primordial light-element abundances. We present here details of a new suite of codes for evaluating such effects, including a new treatment based on PYTHIA of the evolution of showers induced by hadronic decays of massive, unstable particles such as a gravitino. We also develop an analytical treatment of non-thermal hadron propagation in the early universe, and use this to derive analytical estimates for light-element production and in turn on decaying particle lifetimes and abundances. We then consider specifically the case of an unstable massive gravitino within the constrained minimal supersymmetric extension of the Standard Model (CMSSM). We present upper limits on its possible primordial abundance before decay for different possible gravitino masses, with CMSSM parameters along strips where the lightest neutralino provides all the astrophysical cold dark matter density. We do not find any CMSSM solution to the cosmological Li7 problem for small m_{3/2}. Discounting this, for m_{1/2} ~ 500 GeV and tan beta = 10 the other light-element abundances impose an upper limit m_{3/2} n_{3/2}/n_\gamma < 3 \times 10^{-12} GeV to < 2 \times 10^{-13} GeV for m_{3/2} = 250 GeV to 1 TeV, which is similar in both the coannihilation and focus-point strips and somewhat weaker for tan beta = 50, particularly for larger m_{1/2}. The constraints also weaken in general for larger m_{3/2}, and for m_{3/2} > 3 TeV we find a narrow range of m_{3/2} n_{3/2}/n_\gamma, at values which increase with m_{3/2}, where the Li7 abundance is marginally compatible with the other light-element abundances.

We have used a stringy model of quantum space-time foam to suggest that the vacuum may exhibit a non-trivial refractive index depending linearly on gamma-ray energy: eta -1 ~ E_gamma/M_QG1, where M_QG1 is some mass scale typical of quantum gravity that may be ~ 10^18 GeV: see Phys. Lett. B 665, 412 (2008) and references therein. The MAGIC, HESS and Fermi gamma-ray telescopes have recently probed the possible existence of such an energy-dependent vacuum refractive index. All find indications of time-lags for higher-energy photons, but cannot exclude the possibility that they are due to intrinsic delays at the sources. However, the MAGIC and HESS observations of time-lags in emissions from AGNs Mkn 501 and PKS 2155-304 are compatible with each other and a refractive index depending linearly on the gamma-ray energy, with M_QG1 ~ 10^18 GeV. We combine their results to estimate the time-lag Delta t to be expected for the highest-energy photon from GRB 080916c measured by the Fermi telescope, which has an energy ~ 13.2 GeV, assuming the redshift z = 4.2 \pm 0.3 measured by GROND. In the case of a refractive index depending linearly on the gamma-ray energy we predict Delta t = 25 \pm 11 s. This is compatible with the time-lag Delta t <= 16.5 s reported by the Fermi Collaboration, whereas the time-lag would be negligible in the case of a refractive index depending quadratically on the gamma-ray energy. We suggest a strategy for future observations that could distinguish between a quantum-gravitational effect and other interpretations of the time-lags observed by the MAGIC, HESS and Fermi gamma-ray telescopes.

We analyze the timing of photons observed by the MAGIC telescope during a flare of the active galactic nucleus Mkn 501 for a possible correlation with energy, as suggested by some models of quantum gravity (QG), which predict a vacuum refractive index \simeq 1 + (E/M_{QGn})^n, n = 1,2. Parametrizing the delay between gamma-rays of different energies as \Delta t =\pm\tau_l E or \Delta t =\pm\tau_q E^2, we find \tau_l=(0.030\pm0.012) s/GeV at the 2.5-sigma level, and \tau_q=(3.71\pm2.57)x10^{-6} s/GeV^2, respectively. We use these results to establish lower limits M_{QG1} > 0.21x10^{18} GeV and M_{QG2} > 0.26x10^{11} GeV at the 95% C.L. Monte Carlo studies confirm the MAGIC sensitivity to propagation effects at these levels. Thermal plasma effects in the source are negligible, but we cannot exclude the importance of some other source effect.

In view of the imminent start of the LHC experimental programme, we use the available indirect experimental and cosmological information to estimate the likely range of parameters of the constrained minimal supersymmetric extension of the Standard Model (CMSSM), using a Markov-chain Monte Carlo (MCMC) technique to sample the parameter space. The 95% confidence-level area in the (m_0, m_1/2) plane of the CMSSM lies largely within the region that could be explored with 1/fb of integrated luminosity at 14 TeV, and much of the 68% confidence-level area lies within the region that could be explored with 50/pb of integrated luminosity at 10 TeV. A same-sign dilepton signal could well be visible in most of the 68% confidence-level area with 1/fb of integrated luminosity at 14 TeV. We discuss the sensitivities of the preferred ranges to variations in the most relevant indirect experimental and cosmological constraints and also to deviations from the universality of the supersymmetry-breaking contributions to the masses of the Higgs bosons.

We analyze scenarios in which some flavour of sneutrino is the next-to-lightest supersymmetric particle (NLSP), assuming that the gravitino is the lightest supersymmetric particle (LSP) and provides the cold dark matter. Such scenarios do not arise in the constrained supersymmetric extension of the Standard Model (CMSSM) with universal gaugino and scalar masses input at the GUT scale. However, models with non-universal Higgs masses (NUHM) do allow scenarios with a sneutrino NLSP, which are quite generic. We illustrate how such scenarios may arise, analyze the possible metastable sneutrino lifetime, and explore the theoretical, phenomenological, experimental and cosmological constraints on such scenarios. We also discuss the collider signatures of such scenarios, how they may be distinguished from neutralino LSP scenarios, and how different flavours of sneutrino NLSP may be distinguished.

The LHC will probe the nature of the vacuum that determines the properties of particles and the forces between them. Of particular importance is the fact that our current theories allow the Universe to be trapped in a metastable vacuum, which may decay in the distant future, changing the nature of matter. This could be the case in the Standard Model if the LHC finds the Higgs boson to be light. Supersymmetry is one favoured extension of the Standard Model which one might invoke to try to avoid such instability. However, many supersymmetric models are also condemned to vacuum decay for different reasons. The LHC will be able to distinguish between different supersymmetric models, thereby testing the stability of the vacuum, and foretelling the fate of the Universe.

It has been suggested that energetic photons propagating in vacuo should experience a non-trivial refractive index due to the foamy structure of space-time induced by quantum-gravitational fluctuations. The sensitivity of recent astrophysical observations, particularly of AGN Mk501 by the MAGIC Collaboration, approaches the Planck scale for a refractive index depending linearly on the photon energy. We present here a new derivation of this quantum-gravitational vacuum refraction index, based on a stringy analogue of the interaction of a photon with internal degrees of freedom in a conventional medium. We model the space-time foam as a gas of D-particles in the bulk space-time of a higher-dimensional cosmology where the observable Universe is a D3-brane. The interaction of an open string representing a photon with a D-particle stretches and excites the string, which subsequently decays and re-emits the photon with a time delay that increases linearly with the photon energy and is related to stringy uncertainty principles. We relate this derivation to other descriptions of the quantum-gravitational refractive index in vacuo.

We review the uncertainties in the spin-independent and -dependent elastic scattering cross sections of supersymmetric dark matter particles on protons and neutrons. We propagate the uncertainties in quark masses and hadronic matrix elements that are related to the $\pi$-nucleon $\sigma$ term and the spin content of the nucleon. By far the largest single uncertainty is that in spin-independent scattering induced by our ignorance of the $<N | {\bar q} q | N>$ matrix elements linked to the $\pi$-nucleon $\sigma$ term, which affects the ratio of cross sections on proton and neutron targets as well as their absolute values. This uncertainty is already impacting the interpretations of experimental searches for cold dark matter. {\it We plead for an experimental campaign to determine better the $\pi$-nucleon $\sigma$ term.} Uncertainties in the spin content of the proton affect significantly, but less strongly, the calculation of rates used in indirect searches.

We consider the potentials of the LHC and a linear e^+e^- collider (LC) for discovering supersymmetric particles in variants of the MSSM with soft supersymmetry-breaking mass parameters constrained to be universal at the GUT scale (CMSSM) or at some lower scale M_{in} (GUT-less models), as may occur in some scenarios with mirage unification. Whereas the LHC should be able to discover squarks and/or gluinos along all the CMSSM coannihilation strip where the relic neutralino LSP density lies within the range favoured for cold dark matter, many GUT-less models could escape LHC detection. In particular, if M_{in} < 10^{11} GeV, the LHC would not detect sparticles if the relic density lies within the favoured range. For any given discovery of supersymmetry at the LHC, in such GUT-less models the lightest neutralino mass and hence the threshold for sparticle pair production at a LC increases as M_{in} decreases, and the CMSSM offers the best prospects for measuring sparticles at a LC. For example, if the LHC discovers sparticles with 1 fb^{-1} of data, within the CMSSM a centre-of-mass energy of 600 GeV would suffice for a LC to to produce pairs of neutralinos, if they provide the cold dark matter, whereas over 1 TeV might be required in a general GUT-less model. These required energies increase to 800 GeV in the CMSSM and 1.4 TeV in GUT-less models if the LHC requires 10 fb^{-1} to discover supersymmetry.

We correct the fitting formula used in refs. [1,2] to obtain a robust limit on a violation of Lorentz invariance that depends linearly on the photon energy. The correction leads to a slight increase of the limit on the scale of the violation, to M > 1.4 x 10^{16} GeV.

We study models in which supersymmetry breaking appears at an intermediate scale, M_{in}, below the GUT scale. We assume that the soft supersymmetry-breaking parameters of the MSSM are universal at M_{in}, and analyze the morphology of the constraints from cosmology and collider experiments on the allowed regions of parameter space as M_{in} is reduced from the GUT scale. We present separate analyses of the (m_{1/2},m_0) planes for tan(beta)=10 and tan(beta)=50, as well as a discussion of non-zero trilinear couplings, A_0. Specific scenarios where the gaugino and scalar masses appear to be universal below the GUT scale have been found in mirage-mediation models, which we also address here. We demand that the lightest neutralino be the LSP, and that the relic neutralino density not conflict with measurements by WMAP and other observations. At moderate values of M_{in}, we find that the allowed regions of the (m_{1/2},m_0) plane are squeezed by the requirements of electroweak symmetry breaking and that the lightest neutralino be the LSP, whereas the constraint on the relic density is less severe. At very low M_{in}, the electroweak vacuum conditions become the dominant constraint, and a secondary source of astrophysical cold dark matter would be necessary to explain the measured relic density for nearly all values of the soft SUSY-breaking parameters and tan(beta). We calculate the neutralino-nucleon cross sections for viable scenarios and compare them with the present and projected limits from direct dark matter searches.

We analyze the possibility that the lighter stop {\tilde t_1} could be the next-to-lightest supersymmetric particle (NLSP) in models where the gravitino is the lightest supersymmetric particle (LSP). We do not find any possibility for a stop NLSP in the constrained MSSM with universal input soft supersymmetry-breaking masses at the GUT scale (CMSSM), but do find small allowed regions in models with non-universal Higgs masses (NUHM). We discuss the cosmological evolution of stop hadrons. Most {\tilde t_1}qq `sbaryons' and the corresponding `antisbaryons' annihilate with conventional antibaryons and baryons into {\tilde t_1}{\bar q} `mesinos' and the corresponding `antimesinos', respectively, shortly after the quark-hadron transition in the early Universe, and most mesinos and antimesinos subsequently annihilate. As a result, insufficient metastable charged stop hadrons survive to alter Big Bang nucleosynthesis.

We calculate the spectra of ultra-high-energy cosmic rays (UHECRs) in an explicit top-down model based on the decays of metastable neutral `crypton' states in a flipped SU(5) string model. For each of the eight specific 10th-order superpotential operators that might dominate crypton decays, we calculate the spectra of both protons and photons, using a code incorporating supersymmetric evolution of the injected spectra. For all the decay operators, the total UHECR spectra are compatible with the available data. Also, the fractions of photons are compatible with all the published upper limits, but may be detectable in future experiments.

ESO and ESA agreed to establish a number of Working Groups to explore possible synergies between these two major European astronomical institutions. This Working Group's mandate was to concentrate on fundamental questions in cosmology, and the scope for tackling these in Europe over the next ~15 years. One major resulting recommendation concerns the provision of new generations of imaging survey, where the image quality and near-IR sensitivity that can be attained only in space are naturally matched by ground-based imaging and spectroscopy to yield massive datasets with well-understood photometric redshifts (photo-z's). Such information is essential for a range of new cosmological tests using gravitational lensing, large-scale structure, clusters of galaxies, and supernovae. Great scope in future cosmology also exists for ELT studies of the intergalactic medium and space-based studies of the CMB and gravitational waves; here the synergy is less direct, but these areas will remain of the highest mutual interest to the agencies. All these recommended facilities will produce vast datasets of general applicability, which will have a tremendous impact on broad areas of astronomy.

We discuss fits of cosmological dark energy models to the available data on high-redshift supernovae. We consider a conventional model with Cold Dark Matter and a cosmological constant (LambdaCDM), a model invoking super-horizon perturbations (SHCDM) and models based on Liouville strings in which dark energy is provided by a rolling dilaton field (Q-cosmology). We show that a complete treatment of Q-cosmology requires a careful discussion of non-equilibrium situations (off-shell effects). The two main high-redshift supernova data sets give compatible constraints on LambdaCDM and the other models. We recover the well-known result that LambdaCDM fits very well the combined supernova data sets, as does the super-horizon model. We discuss the model-dependent off-shell corrections to the Q-cosmology model that are relevant to the supernova data, and show that this model fits the data equally well. This analysis could be extended to other aspects of cosmological phenomenology, in particular to the CMB and Baryon Acoustic Oscillations, which have so far been treated using on-shell models.

If the gravitino is the lightest supersymmetric particle and the long-lived next-to-lightest sparticle (NSP) is the stau, the charged partner of the tau lepton, it may be metastable and form bound states with several nuclei. These bound states may affect the cosmological abundances of Li6 and Li7 by enhancing nuclear rates that would otherwise be strongly suppressed. We consider the effects of these enhanced rates on the final abundances produced in Big-Bang nucleosynthesis (BBN), including injections of both electromagnetic and hadronic energy during and after BBN. We calculate the dominant two- and three-body decays of both neutralino and stau NSPs, and model the electromagnetic and hadronic decay products using the PYTHIA event generator and a cascade equation. Generically, the introduction of bound states drives light element abundances further from their observed values; however, for small regions of parameter space bound state effects can bring lithium abundances in particular in better accord with observations. We show that in regions where the stau is the NSP with a lifetime longer than 10^3-10^4 s, the abundances of Li6 and Li7 are far in excess of those allowed by observations. For shorter lifetimes of order 1000 s, we comment on the possibility in minimal supersymmetric and supergravity models that stau decays could reduce the Li7 abundance from standard BBN values while at the same time enhancing the Li6 abundance.

The three-year data from WMAP are in stunning agreement with the simplest possible quadratic potential for chaotic inflation, as well as with new or symmetry-breaking inflation. We investigate the possibilities for incorporating these potentials within supergravity, particularly of the no-scale type that is motivated by string theory. Models with inflation driven by the matter sector may be constructed in no-scale supergravity, if the moduli are assumed to be stabilised by some higher-scale dynamics and at the expense of some fine-tuning. We discuss specific scenarios for stabilising the moduli via either D- or F-terms in the effective potential, and survey possible inflationary models in the presence of D-term stabilisation.

We constrain the possibility of a non-trivial refractive index in free space corresponding to an energy-dependent velocity of light: c(E) \simeq c_0 (1 - E/M), where M is a mass scale that might represent effect of quantum-gravitational space-time foam, using the arrival times of sharp features observed in the intensities of radiation with different energies from a large sample of gamma-ray bursters (GRBs) with known redshifts. We use wavelet techniques to identify genuine features, which we confirm in simulations with artificial added noise. Using the weighted averages of the time-lags calculated using correlated features in all the GRB light curves, we find a systematic tendency for more energetic photons to arrive earlier. However, there is a very strong correlation between the parameters characterizing an intrinsic time-lag at the source and a distance-dependent propagation effect. Moreover, the significance of the earlier arrival times is less evident for a subsample of more robust spectral structures. Allowing for intrinsic stochastic time-lags in these features, we establish a statistically robust lower limit: M > 0.9x10^{16} GeV on the scale of violation of Lorentz invariance.

The `Little Bangs' made in particle collider experiments reproduce the conditions in the Big Bang when the age of the Universe was a fraction of a second. It is thought that matter was generated, the structures in the Universe were formed and cold dark matter froze out during this very early epoch when the equation of state of the Universe was dominated by the quark-gluon plasma (QGP). Future Little Bangs may reveal the mechanism of matter generation and the nature of cold dark matter. Knowledge of the QGP will be an essential ingredient in quantitative understanding of the very early Universe.

We reconsider the effects of unstable particles on the production and destruction of the primordial light elements, with a view to reconciling the high primordial Li7 abundance deduced from Big Bang Nucleosynthesis (BBN), as implied by the baryon-to-photon ratio now inferred from the anisotropies of the Cosmic Microwave Background (CMB), with the lower abundance of Li7 observed in halo stars. The potential destruction of Li7 is strongly constrained by observations of Deuterium (D), He3 and Li6. We identify ranges for the unstable particle abundance and lifetime which would deplete Li7 while remaining consistent with the abundance of Li6. However, in these regions either the D abundance is unacceptably low or the ratio He3/D is unacceptably large. We conclude that late particle decay is unable to explain both the discrepancy of the calculated Li7 abundance and the observed Li7 plateau. In the context of supersymmetric theories with neutralino or gravitino dark matter, we display the corresponding light-element constraints on the model parameters.

Liouville string theory is a natural framework for discussing the non-equilibrium evolution of the Universe. It enables non-critical strings to be treated in mathematically consistent manner, in which target time is identified with a world-sheet renormalization-group scale parameter, preserving target-space general coordinate invariance and the existence of an S-matrix. We review our proposals for a unified treatment of inflation and the current acceleration of the Universe. We link the current acceleration of the Universe with the value of the string coupling. In such a scenario, the dilaton plays an essential background role, driving the acceleration of the Universe during the present era after decoupling as a constant during inflation.

We compare updated predictions for the elastic scattering of supersymmetric neutralino dark matter with the improved experimental upper limit recently published by CDMS II. We take into account the possibility that the \pi-nucleon \Sigma term may be somewhat larger than was previously considered plausible, as may be supported by the masses of exotic baryons reported recently. We also incorporate the new central value of m_t, which affects indirectly constraints on the supersymmetric parameter space, for example via calculations of the relic density. Even if a large value of \Sigma is assumed, the CDMS II data currently exclude only small parts of the parameter space in the constrained MSSM (CMSSM) with universal soft supersymmetry-breaking Higgs, squark and slepton masses. None of the previously-proposed CMSSM benchmark scenarios is excluded for any value of \Sigma, and the CDMS II data do not impinge on the domains of the CMSSM parameter space favoured at the 90 % confidence level in a recent likelihood analysis. However, some models with non-universal Higgs, squark and slepton masses and neutralino masses \lappeq 700 GeV are excluded by the CDMS II data.

In the context of a cosmological string model describing the propagation of strings in a time-dependent Robertson-Walker background space-time, we show that the asymptotic acceleration of the Universe can be identified with the square of the string coupling. This allows for a direct measurement of the ten-dimensional string coupling using cosmological data. We conjecture that this is a generic feature of a class of non-critical string models that approach asymptotically a conformal (critical) sigma model whose target space is a four-dimensional space-time with a dilaton background that is linear in sigma-model time. The relation between the cosmic acceleration and the string coupling does not apply in critical strings with constant dilaton fields in four dimensions.

We present a specific model for cosmological inflation driven by the Liouville field in a non-critical supersymmetric string framework, in which the departure from criticality is due to open strings stretched between the two moving Type-II 5-branes. We use WMAP and other data on fluctuations in the cosmic microwave background to fix parameters of the model, such as the relative separation and velocity of the 5-branes, respecting also the constraints imposed by data on light propagation from distant gamma-ray bursters. The model also suggests a small, relaxing component in the present vacuum energy that may accommodate the breaking of supersymmetry.

Recent cosmological observations of unprecented accuracy, by WMAP in particular, have established a `Standard Model' of cosmology, just as LEP established the Standard Model of particle physics. Both Standard Models raise open questions whose answers are likely to be linked. The most fundamental problems in both particle physics and cosmology will be resolved only within a framework for Quantum Gravity, for which the only game in town is string theory. We discuss novel ways to model cosmological inflation and late acceleration in a non-critical string approach, and discuss possible astrophysical tests.

Live 60Fe has recently been detected in a deep-ocean ferromanganese crust, isolated in layers dating from about 3 Myr ago. Since 60Fe has a mean life of 2.2 Myr, a near-Earth supernova is the only likely source for such a signal, and we explore here the consequences of a supernova origin. We combine the 60Fe data with several supernova nucleosynthesis models to calculate the supernova distance as a function of progenitor mass, finding an allowed range of 15-120 pc. We also predict the signals expected for several other radioisotopes, which are independent of the supernova distance. Species likely to be present near or above background levels are 10Be, 26Al, 53Mn, 182Hf and 244Pu. Of these, 182Hf and 244Pu are nearly background-free, presenting the best opportunities to provide strong confirmation of the supernova origin of the 60Fe signal, and to demonstrate that at least some supernovae are the source for the r-process. The accuracies of our predictions are hampered by large uncertainties in the predicted 60Fe yields for supernovae of different masses, so the new crust data motivate a redoubled theoretical attack on this problem.

We discuss the prospects for detecting supersymmetric particles in variants of the minimal supersymmetric extension of the Standard Model (MSSM), in light of laboratory and cosmological constraints. We first assume that the lightest supersymmetric particle (LSP) is the lightest neutralino chi, and present scatter plots of the masses of the two lightest visible supersymmetric particles when the input scalar and gaugino masses are constrained to be universal (CMSSM), when the input Higgs scalar masses are non-universal (NUHM), and when the squark and slepton masses are also non-universal and the MSSM is regarded as a low-energy effective field theory valid up to the GUT scale (LEEST) or just up to 10 TeV (LEEST10). We then present similar plots in various scenarios when the LSP is the gravitino. We compare the prospects for detecting supersymmetry at linear colliders (LCs) of various energies, at the LHC, and as astrophysical dark matter. We find that, whilst a LC with a centre-of-mass energy E_{CM} <= 1000 GeV has some chance of discovering the lightest and next-to-lightest visible supersymmetric particles, E_{CM} >= 3000 GeV would be required to `guarantee' finding supersymmetry in the neutralino LSP scenarios studied, and an even higher E_{CM} might be required in certain gravitino dark matter scenarios. Direct dark matter experiments could explore part of the low-mass neutralino LSP region, but would not reveal all the models accessible to a low-energy LC.

Cryptons are metastable bound states of fractionally-charged particles that arise generically in the hidden sectors of models derived from heterotic string. We study their properties and decay modes in a specific flipped SU(5) model with long-lived four-particle spin-zero bound states called {\it tetrons}. We show that the neutral tetrons are metastable, and exhibit the tenth-order non-renormalizable superpotential operators responsible for their dominant decays. By analogy with QCD, we expect charged tetrons to be somewhat heavier, and to decay relatively rapidly via lower-order interactions that we also exhibit. The expected masses and lifetimes of the neutral tetrons make them good candidates for cold dark matter (CDM), and a potential source of the ultra-high energy cosmic rays (UHECRs) which have been observed, whereas the charged tetrons would have decayed in the early Universe.

The interactions of different particle species with the foamy space-time fluctuations expected in quantum gravity theories may not be universal, in which case different types of energetic particles may violate Lorentz invariance by varying amounts, violating the equivalence principle. We illustrate this possibility in two different models of space-time foam based on D-particle fluctuations in either flat Minkowski space or a stack of intersecting D-branes. Both models suggest that Lorentz invariance could be violated for energetic particles that do not carry conserved charges, such as photons, whereas charged particles such electrons would propagate in a Lorentz-inavariant way. The D-brane model further suggests that gluon propagation might violate Lorentz invariance, but not neutrinos. We argue that these conclusions hold at both the tree (lowest-genus) and loop (higher-genus) levels, and discuss their implications for the phenomenology of quantum gravity.

Photons may evade a synchrotron radiation constraint on quantum gravity by violating the equivalence principle.

The seesaw model of neutrinos might explain the size, age, flatness and near-homogeneity of the Universe via sneutrino inflation, as well as explaining the origin of matter via leptogenesis. The sneutrino inflation hypothesis makes specific, testable predictions for cosmic microwave background observables, which are compatible with the first release of data from WMAP, and for flavour-violating charged-lepton decays. In particular, mu -> e gamma should occur with a branching ratio very close to the present experimental upper limit, whilst tau -> mu gamma and e gamma should occur further below the present limits.

It has been argued by Jacobson, Liberati and Mattingly that synchrotron radiation from the Crab Nebula imposes a stringent constraint on any modification of the dispersion relations of the electron that might be induced by quantum gravity. We supplement their analysis by deriving the spectrum of synchrotron radiation from the coupling of an electrically-charged particle to an external magnetic fields in the presence of quantum-gravity effects of the general form $(E/M_{QG})^\alpha$. We find that the synchrotron constraint from the Crab Nebula practically excludes $\alpha \lsim 1.74$ for $M_{QG} \sim m_P = 1.2 \times 10^{19}$ GeV. On the other hand, this analysis does not constrain any modification of the dispersion relation of the photon that might be induced by quantum gravity. We point out that such quantum-gravity effects need not obey the equivalence principle, a point exemplified by the Liouville-string D-particle model of space-time foam. This model suggests a linear modification of the dispersion relation for the photon, but not for the electron, and hence is compatible with known constraints from the Crab Nebula and elsewhere.

Big-Bang cosmology and ideas for possible physics beyond the Standard Model of particle physics are introduced. The density budget of the Universe is audited, and the issues involved in calculating the baryon density from microphysics are mentioned, as is the role of cold dark matter in the formation of cosmological structures. Candidates for cold dark matter are introduced, with particular attention to the lightest supersymmetric particle and metastable superheavy relics. Prospects for detecting supersymmetric dark matter in non-accelerator experiments are assessed, and the possible role of decays in generating ultra-high-energy cosmic rays is discussed. More details of these and other astroparticle topics are presented during the rest of this Summer Institute.

A National Research Council study on connecting quarks with the cosmos has recently posed a number of the more important open questions at the interface between particle physics and cosmology. These questions include the nature of dark matter and dark energy, how the Universe began, modifications to gravity, the effects of neutrinos on the Universe, how cosmic accelerators work, and whether there are new states of matter at high density and pressure. These questions are discussed in the context of the talks presented at this Summer Institute.

Specific models of supersymmetry breaking predict relations between the trilinear and bilinear soft supersymmetry breaking parameters A_0 and B_0 at the input scale. In such models, the value of tan beta can be calculated as a function of the scalar masses m_0 and the gaugino masses m_{1/2}, which we assume to be universal. The experimental constraints on sparticle and Higgs masses, b to s gamma decay and the cold dark matter density Omega_{CDM} h^2 can then be used to constrain tan beta in such specific models of supersymmetry breaking. In the simplest Polonyi model with A_0 = (3 - sqrt{3})m_0 = B_0 + m_0, we find 11 < tan beta < 20 (tan beta ~ 4.15) for mu > 0 (mu < 0). We also discuss other models with A_0 = B_0 + m_0, finding that only the range -1.9 < A_0/m_0 < 2.5 is allowed for mu > 0, and the range 1.25 < A_0/m_0 < 4.8 for mu < 0. In these models, we find no solutions in the rapid-annihilation `funnels' or in the `focus-point' region. We also discuss the allowed range of tan beta in the no-scale model with A_0 = B_0 = 0. In all these models, most of the allowed regions are in the chi - stau_1 coannihilation `tail'.

In the first Lecture, the Big Bang and the Standard Model of particle physics are introduced, as well as the structure of the latter and open issues beyond it. Neutrino physics is discussed in the second Lecture, with emphasis on models for neutrino masses and oscillations. The third Lecture is devoted to supersymmetry, including the prospects for discovering it at accelerators or as cold dark matter. Inflation is reviewed from the viewpoint of particle physics in the fourth Lecture, including simple models with a single scalar inflaton field: the possibility that this might be a sneutrino is proposed. Finally, the fifth Lecture is devoted to topics further beyond the Standard Model, such as grand unification, baryo- and leptogenesis - that might be due to sneutrino inflaton decays - and ultra-high-energy cosmic rays - that might be due to the decays of metastable superheavy dark matter particles.

This paper reviews the progress reported at this Royal Society Discussion Meeting and advertizes some possible future directions in our drive to understand dark matter and dark energy. Additionally, a first attempt is made to place in context the exciting new results from the WMAP satellite, which were published shortly after this Meeting. In the first part of this review, pieces of observational evidence shown here that bear on the amounts of dark matter and dark energy are reviewed. Subsequently, particle candidates for dark matter are mentioned, and detection strategies are discussed. Finally, ideas are presented for calculating the amounts of dark matter and dark energy, and possibly relating them to laboratory data.

We re-examine the parameter space of the constrained minimal supersymmetric extension of the Standard Model (CMSSM), taking account of the restricted range of \Omega_{CDM} h^2 consistent with the WMAP data. This provides a significantly reduced upper limit on the mass of the lightest supersymmetric particle LSP: m_\chi < 500 GeV for \tan \beta < 45 and \mu > 0, or \tan \beta < 30 and \mu < 0, thereby improving the prospects for measuring supersymmetry at the LHC, and increasing the likelihood that a 1-TeV linear e^+ e^- collider would be able to measure the properties of some supersymmetric particles.

In some models of quantum gravity, space-time is thought to have a foamy structure with non-trivial optical properties. We probe the possibility that photons propagating in vacuum may exhibit a non-trivial refractive index, by analyzing the times of flight of radiation from gamma-ray bursters (GRBs) with known redshifts. We use a wavelet shrinkage procedure for noise removal and a wavelet `zoom' technique to define with high accuracy the timings of sharp transitions in GRB light curves, thereby optimizing the sensitivity of experimental probes of any energy dependence of the velocity of light. We apply these wavelet techniques to 64 ms and TTE data from BATSE, and also to OSSE data. A search for time lags between sharp transients in GRB light curves in different energy bands yields the lower limit $M \ge 6.9 \cdot 10^{15}$ GeV on the quantum-gravity scale in any model with a linear dependence of the velocity of light $~ E/M$. We also present a limit on any quadratic dependence.

We calculate dark matter scattering rates in the minimal supersymmetric extension of the Standard Model (MSSM), allowing the soft supersymmetry-breaking masses of the Higgs multiplets, m_{1,2}, to be non-universal (NUHM). Compared with the constrained MSSM (CMSSM) in which m_{1,2} are required to be equal to the soft supersymmetry-breaking masses m_0 of the squark and slepton masses, we find that the elastic scattering cross sections may be up to two orders of magnitude larger than values in the CMSSM for similar LSP masses. We find the following preferred ranges for the spin-independent cross section: 10^{-6} pb \ga \sigma_{SI} \ga 10^{-10} pb, and for the spin-dependent cross section: 10^{-3} pb \ga \sigma_{SD}, with the lower bound on \sigma_{SI} dependent on using the putative constraint from the muon anomalous magnetic moment. We stress the importance of incorporating accelerator and dark matter constraints in restricting the NUHM parameter space, and also of requiring that no undesirable vacuum appear below the GUT scale. In particular, values of the spin-independent cross section another order of magnitude larger would appear to be allowed, for small \tan \beta, if the GUT vacuum stability requirement were relaxed, and much lower cross-section values would be permitted if the muon anomalous magnetic moment constraint were dropped.

We extend previous analyses of the violation of Lorentz invariance induced in a non-critical string model of quantum space-time foam, discussing the propagation of low-energy particles through a distribution of non-relativistic D-particles.We argue that nuclear and atomic physics experiments do not constitute sensitive probes of this approach to quantum gravity due to a difference in the dispersion relations for massive probes as compared to those for massless ones, predicted by the model.

We explore the parameter space of the minimal supersymmetric extension of the Standard Model (MSSM), allowing the soft supersymmetry-breaking masses of the Higgs multiplets, m_{1,2}, to be non-universal (NUHM). Compared with the constrained MSSM (CMSSM) in which m_{1,2} are required to be equal to the soft supersymmetry-breaking masses m_0 of the squark and slepton masses, the Higgs mixing parameter mu and the pseudoscalar Higgs mass m_A, which are calculated in the CMSSM, are free in the NUHM model. We incorporate accelerator and dark matter constraints in determining allowed regions of the (mu, m_A), (mu, M_2) and (m_{1/2}, m_0) planes for selected choices of the other NUHM parameters. In the examples studied, we find that the LSP mass cannot be reduced far below its limit in the CMSSM, whereas m_A may be as small as allowed by LEP for large tan \beta. We present in Appendices details of the calculations of neutralino-slepton, chargino-slepton and neutralino-sneutrino coannihilation needed in our exploration of the NUHM.

We revisit the upper limits on the abundance of unstable massive relic particles provided by the success of Big-Bang Nucleosynthesis calculations. We use the cosmic microwave background data to constrain the baryon-to-photon ratio, and incorporate an extensively updated compilation of cross sections into a new calculation of the network of reactions induced by electromagnetic showers that create and destroy the light elements deuterium, he3, he4, li6 and li7. We derive analytic approximations that complement and check the full numerical calculations. Considerations of the abundances of he4 and li6 exclude exceptional regions of parameter space that would otherwise have been permitted by deuterium alone. We illustrate our results by applying them to massive gravitinos. If they weigh ~100 GeV, their primordial abundance should have been below about 10^{-13} of the total entropy. This would imply an upper limit on the reheating temperature of a few times 10^7 GeV, which could be a potential difficulty for some models of inflation. We discuss possible ways of evading this problem.

Interfaces between high-energy physics, astrophysics and cosmology are reviewed, with particular emphasis on the important roles played by high-energy cosmic-ray physics. These include the understanding of atmospheric neutrinos, the search for massive cold dark matter particles and possible tests of models of quantum gravity. In return, experiments at the LHC may be useful for refining models of ultra-high-energy cosmic rays, and thereby contributing indirectly to understanding their origin. Only future experiments will be able to tell whether these are due to some bottom-up astrophysical mechanism or some top-down cosmological mechanism.

Current and future experiments at CERN are reviewed,with emphasis on those relevant to astrophysics and cosmology. These include experiments related to nuclear astrophysics, matter-antimatter asymmetry, dark matter, axions, gravitational waves, cosmic rays, neutrino oscillations, inflation, neutron stars and the quark-gluon plasma. The centrepiece of CERN's future programme is the LHC, but some ideas for perspectives after the LHC are also presented.

Without assuming that Higgs masses have the same values as other scalar masses at the input GUT scale, we combine constraints on the minimal supersymmetric extension of the Standard Model (MSSM) coming from the cold dark matter density with the limits from direct searches at accelerators such as LEP, indirect measurements such as b to s gamma decay and the anomalous magnetic moment of the muon. The requirement that Higgs masses-squared be positive at the GUT scale imposes important restrictions on the MSSM parameter space, as does the requirement that the LSP be neutral. We analyze the interplay of these constraints in the (mu, m_A), (mu, m_{1/2}), (m_{1/2}, m_0) and (m_A, tan beta) planes. These exhibit new features not seen in the corresponding planes in the constrained MSSM in which universality is extended to Higgs masses.

Salient aspects of the meeting are summarized, including our knowledge of dark matter at different cosmological and astrophysical distance scales, ranging from large-scale structure to the cores of galaxies, and our speculations on particle candidates for dark matter, including neutrinos, neutralinos, axinos, gravitinos and cryptons. Comments are also made on prospects for detecting dark matter particles and on the dark energy problem.

We present detailed calculations of the neutralino-stop coannihilation channels that have the largest impact on the relic neutralino density in the constrained minimal supersymmetric extension of the Standard Model (CMSSM), in which scalar masses m_0, gaugino masses m_1/2 and the trilinear soft supersymmetry-breaking parameters A_0 are each assumed to be universal at some input grand unification scale. The most important stop-stop* and stop-stop annihilation channels are also calculated, as well as stop-slepton coannihilation channels. We illustrate the importance of these new coannihilation calculations when A_0 is relatively large. While they do not increase the range of m_1/2 and hence neutralino mass allowed by cosmology, these coannihilation channels do open up new `tails' of parameter space extending to larger values of m_0.

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.

We explore the allowed ranges of cross sections for the elastic scattering of neutralinos \chi on nucleons in the constrained minimal supersymmetric extension of the Standard Model (CMSSM), in which scalar and gaugino masses are each assumed to be universal at some input grand unification scale. We extend previous calculations to larger \tan \beta and investigate the limits imposed by the recent LEP lower limit on the mass of the Higgs boson and by b \to s \gamma, and those suggested by g_\mu - 2. The Higgs limit and b \to s \gamma provide upper limits on the cross section, particularly at small and large \tan \beta, respectively, and the value of g_\mu - 2 suggests a lower limit on the cross section for \mu > 0. The spin-independent nucleon cross section is restricted to the range 6 \times 10^{-8} pb > \sigma_{SI} > 2 \times 10^{-10} pb for \mu > 0, and the spin-dependent nucleon cross section to the range 10^{-5} pb > \sigma_{SD} > 2 \times 10^{-7} pb. Lower values are allowed if \mu <0.

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.

It has recently been pointed out that generic models of quantum gravity must contend with severe phenomenological constraints imposed by gravitational Cerenkov radiation, neutrino oscillations and the cosmic microwave background radiation. We show how the non-critical Liouville-string model of quantum gravity we have proposed escapes these constraints. It gives energetic particles subluminal velocities, obviating the danger of gravitational Cerenkov radiation. The effect on neutrino propagation is naturally flavour-independent, obviating any impact on oscillation phenomenology. Deviations from the expected black-body spectrum and the effects of time delays and stochastic fluctuations in the propagation of cosmic microwave background photons are negligible, as are their effects on observable spectral lines from high-redshift astrophysical objects.

We propose a new set of supersymmetric benchmark scenarios, taking into account the constraints from LEP, b to s gamma, g_mu - 2 and cosmology. We work in the context of the constrained MSSM (CMSSM) with universal soft supersymetry-breaking masses and assume that R parity is conserved. We propose benchmark points that exemplify the different generic possibilities, including focus-point models, points where coannihilation effects on the relic density are important, and points with rapid relic annihilation via direct-channel Higgs poles. We discuss the principal decays and signatures of the different classes of benchmark scenarios, and make initial estimates of the physics reaches of different accelerators, including the Tevatron collider, the LHC, and e+ e- colliders in the sub- and multi-TeV ranges. We stress the complementarity of hadron and lepton colliders, with the latter favoured for non-strongly-interacting particles and precision measurements. We mention features that could usefully be included in future versions of supersymmetric event generators.

An accelerating Universe can be accommodated naturally within non-critical string theory, in which scattering is described by a superscattering matrix \$ that does not factorize as a product of $S$- and $S^\dagger$-matrix elements and time evolution is described by a modified Liouville equation characteristic of open quantum-mechanical systems. We describe briefly alternative representations in terms of the stochastic Ito and Fokker-Planck equations. The link between the vacuum energy and the departure from criticality is stressed. We give an explicit example in which non-marginal \$tring couplings cause a departure from criticality, and the corresponding cosmological vacuum energy relaxes to zero \`a la quintessence.

After recalling the motivations for expecting supersymmetry to appear at energies below 1 TeV, the reasons why the lightest supersymmetric particle is an ideal candidate for cold dark matter are reviewed from a historical perspective. Recent calculations of the relic density including coannihilations and rapid annihilations through direct-channel Higgs boson poles are presented. The experimental constraints from LEP and elsewhere on supersymmetric dark matter are reviewed, and the prospects for its indirect or direct detection are mentioned. The potential implications of a Higgs boson weighing about 115 GeV and the recent measurement of the anomalous magnetic moment of the muon are summarized.

We extend previous combinations of LEP and cosmological relic density constraints on the parameter space of the constrained MSSM, with universal input supersymmetry-breaking parameters, to large tan beta. We take account of the possibility that the lightest Higgs boson might weigh about 115 GeV, but also retain the possibility that it might be heavier. We include the most recent implementation of the b to s gamma constraint at large tan beta. We refine previous relic density calculations at large tan beta by combining a careful treatment of the direct-channel Higgs poles in annihilation of pairs of neutralinos chi with a complete treatment of chi - stau coannihilation, and discuss carefully uncertainties associated with the mass of the b quark. We find that coannihilation and pole annihilations allow the CMSSM to yield an acceptable relic density at large tan beta, but it is consistent with all the constraints only if m_chi > 140 (180) GeV for mu > 0 (mu < 0) for our default choices m_b(m_b) = 4.25 GeV, m_t = 175 GeV, and A_0 = 0.

We combine the constraint suggested by the recent BNL E821 measurement of the anomalous magnetic moment of the muon on the parameter space of the constrained MSSM (CMSSM) with those provided previously by LEP, the measured rate of b to s gamma decay and the cosmological relic density Omega-hsquared. Our treatment of Omega-hsquared includes carefully the direct-channel Higgs poles in annihilation of pairs of neutralinos chi and a complete analysis of chi - slepton coannihilation. We find excellent consistency between all the constraints for tan beta > 10 and mu > 0, for restricted ranges of the CMSSM parameters m_0 and m_1/2. All the preferred CMSSM parameter space is within reach of the LHC, but may not be accessible to the Tevatron collider, or to a first-generation e^+ e^- linear collider with centre-of-mass energy below 1.2 TeV.

Modelling space-time foam using a non-critical Liouville-string model for the quantum fluctuations of D branes with recoil, we discuss the issues of momentum and energy conservation in particle propagation and interactions. We argue that momentum should be conserved exactly during propagation and on the average during interactions, but that energy is conserved only on the average during propagation and is in general not conserved during particle interactions, because of changes in the background metric. We discuss the possible modification of the GZK cutoff on high-energy cosmic rays, in the light of this energy non-conservation as well as the possible modification of the usual relativistic momentum-energy relation.

Cosmic rays may provide opportunities for probing fundamental physics. For example, ultra-high-energy cosmic rays might originate from the decays of metastable heavy particles, and astrophysical gamma rays can be used to test models of quantum gravity. Both scenarios offer ways to avoid the GZK cut-off.

Assuming that the lightest supersymmetric particle (LSP) is the lightest neutralino, we present a detailed exploration of neutralino-stau coannihilation channels, including analytical expressions and numerical results. We also include neutralino coannihilations with the selectron and smuon. We evaluate the implications of coannihilations for the cosmological relic density of the LSP, which is assumed to be stable, in the constrained minimal supersymmetric extension of the Standard Model (CMSSM), in which the soft supersymmetry-breaking parameters are universal at the supergravity GUT scale. We evaluate the changes due to coannihilations in the region of the MSSM parameter space that is consistent with the cosmological upper limit on the relic LSP density. In particular, we find that the upper limit on the mass of the neutralino is increased from about 200 GeV to about 600 GeV in the CMSSM, and estimate a qualitatively similar increase for gauginos in the general MSSM.

We explore the possible cross sections for the elastic scattering of neutralinos chi on nucleons p,n in the minimal supersymmetric extension of the standard model (MSSM). Universality of the soft supersymmetry-breaking scalar masses for the Higgs multiplets is not assumed, but the MSSM parameters are nevertheless required to lead consistently to an electroweak vacuum. We explore systematically the region of MSSM parameter space where LEP and other accelerator constraints are respected, and the relic neutralino density lies in the range 0.1 < Omega_chi h^2 < 0.3 preferred by cosmology. We also discuss models with Omega_chi h^2 < 0.1, in which case we scale the density of supersymmetric dark matter in our galactic halo by Omega_chi h^2 / 0.1, allowing for the possible existence of some complementary form of cold dark matter. We find values of the cross sections that are considerably lower than the present experimental sensitivities. At low neutralino masses, m_chi < 100 GeV, the cross sections may be somewhat higher than in the constrained MSSM with universal soft Higgs masses, though they are generally lower. In the case of large m_chi, the cross sections we find may be considerably larger than in the constrained model, but still well below the present experimental sensitivity.

Classical superstring vacua have zero vacuum energy and are supersymmetric and Lorentz-invariant. We argue that all these properties may be destroyed when quantum aspects of the interactions between particles and non-perturbative vacuum fluctuations are considered. A toy calculation of string/D-brane interactions using a world-sheet approach indicates that quantum recoil effects - reflecting the gravitational back-reaction on space-time foam due to the propagation of energetic particles - induce non-zero vacuum energy that is linked to supersymmetry breaking and breaks Lorentz invariance. This model of space-time foam also suggests the appearance of microscopic event horizons.

We analyze the accelerator constraints on the parameter space of the Minimal Supersymmetric extension of the Standard Model, comparing those now available from LEP II and anticipating the likely sensitivity of Tevatron Run II. The most important limits are those from searches for charginos, neutralinos and Higgs bosons at LEP, and searches for stop squarks, charginos and neutralinos at the Tevatron Collider. We also incorporate the constraints derived from b --> s + gamma decay, and discuss the relevance of charge- and colour-breaking minima in the effective potential. We combine and compare the different constraints on the Higgs-mixing parameter mu, the gaugino-mass parameter m_{1/2} and the scalar-mass parameter m0, incorporating radiative corrections to the physical particle masses. We focus on the resulting limitations on supersymmetric dark matter, assumed to be the lightest neutralino, incorporating coannihilation effects in the calculation of the relic abundance. We find that m_chi > 51 GeV and tan(beta) > 2.2 if all soft supersymmetry-breaking scalar masses are universal, including those of the Higgs bosons, and that these limits weaken to m_chi > 46 GeV and tan(beta) > 1.9 if non-universal scalar masses are allowed. Light neutralino dark matter cannot be primarily Higgsino in composition.

We examine the cross sections for the elastic scattering of neutralinos $\chi$ on nucleons $p,n$, as functions of $m_\chi$ in the constrained minimal supersymmetric standard model. We find narrow bands of possible values of the cross section, that are considerably lower than some previous estimates. The constrained model is based on the minimal supergravity-inspired framework for the MSSM, with universal scalar and gaugino masses $m_0, m_{1/2}$, and $\mu$ and the MSSM Higgs masses treated as dependent parameters. We explore systematically the region of the $(m_{1/2}, m_0)$ plane where LEP and other accelerator constraints are respected, and the relic neutralino density lies in the range $0.1 \le \Omega_{\chi} h^2 \le 0.3$ preferred by cosmology. We update previous discussions of both the spin-independent and -dependent scattering matrix elements on protons and neutrons, using recent analyses of low-energy hadron experiments.

We suggest that supersymmetric dark matter be used to set the energy scale of a linear $e^+ e^-$ collider. Assuming that the lightest supersymmetric particle (LSP) is a stable neutralino $\chi$, as in many incarnations of the MSSM with conserved R parity, previous calculations that include coannihilation effects have delineated the region of the $(m_{1/2}, m_0)$ plane where the LSP cosmological relic density lies in the preferred range $0.1 \la \Omega_{\chi} h^2 \la 0.3$. We evaluate here the total cross section for $e^+ e^- \to$ visible pairs of supersymmetric particles, for different values of $m_{1/2}$ and $m_0$, and investigate how much of the dark matter region can be explored by $e^+ e^-$ colliders with different centre-of-mass energies $E_{CM}$. We find that a collider with $E_{CM} = 500$ GeV or 1 TeV can only explore part of the cosmological region, and that a collider with $E_{CM} = 1.5$ TeV with sufficient luminosity can explore all of the supersymmetric dark matter region.

The density budget of the Universe is reviewed, and then specific particle candidates for non-bayonic dark matter are introduced, with emphasis on the relevance of cosmic-ray physics. The sizes of the neutrino masses indicated by recent atmospheric and solar neutrino experiments may be too small to contribute much hot dark matter. My favoured candidate for the dominant cold dark matter is the lightest supersymmetric particle, which probably weighs between about 50 GeV and about 600 GeV. Strategies to search for it via cosmic rays due to annihilations in the halo, Sun and Earth, or via direct scattering experiments, are mentioned. Possible superheavy relic particles are also discussed, in particular metastable string- or M-theory cryptons, that may be responsible for the ultra-high-energy cosmic rays. Finally, it is speculated that a non-zero contribution to the cosmological vacuum energy might result from incomplete relaxation of the quantum-gravitational vacuum.

Motivated by the possible experimental opportunities to test quantum gravity via its effects on high-energy neutrinos propagating through space-time foam, we discuss how to incorporate spin structures in our D-brane description of gravitational recoil effects in vacuo. We also point to an interesting analogous condensed-matter system. We use a suitable supersymmetrization of the Born-Infeld action for excited D-brane gravitational backgrounds to argue that energetic fermions may travel slower than the low-energy velocity of light: \delta c / c \sim -E/M. It has been suggested that Gamma-Ray Bursters may emit pulses of neutrinos at energies approaching 10^{19} eV: these would be observable only if M \gsim 10^{27} GeV.

We review the possibility that quantum fluctuations in the structure of space-time at the Planck scale might be subject to experimental probes. We study the effects of space-time foam in an approach inspired by string theory, in which solitonic D-brane excitations are taken into account when considering the ground state. We model the properties of this medium by analyzing the recoil of a D particle which is induced by the scattering of a closed-string state. We find that this recoil causes an energy-dependent perturbation of the background metric, which in turn induces an energy-dependent refractive index in vacuo, and stochastic fluctuations of the light cone. We show how distant astrophysical sources such as Gamma-Ray Bursters (GRBs) may be used to test this possibility, making an illustrative analysis of GRBs whose redshifts have been measured. Within this framework, we also discuss the propagation of massive particles and the possible appearance of cosmological vacuum energy that relaxes towards zero. We also discuss D-brane recoil in models with `large' extra dimensions.

We discuss possible tests of the constancy of the velocity of light using distant astrophysical sources such as gamma-ray bursters (GRBs), Active Galactic Nuclei (AGNs) and pulsars. This speculative quest may be motivated by some models of quantum fluctuations in the space-time background, and we discuss explicitly how an energy-dependent variation in photon velocity \delta c/ c \sim - E / M arises in one particular quantum-gravitational model. We then discuss how data on GRBs may be used to set limits on variations in the velocity of light, which we illustrate using BATSE and OSSE observations of the GRBs that have recently been identified optically and for which precise redshifts are available. We show how a regression analysis can be performed to look for an energy-dependent effect that should correlate with redshift. The present data yield a limit M \gsim 10^{15} GeV for the quantum gravity scale. We discuss the prospects for improving this analysis using future data, and how one might hope to distinguish any positive signal from astrophysical effects associated with the sources.

We examine the effect of light-cone broadening induced by quantum-gravity foam in the context of theories with ``large'' extra dimensions stretching between two parallel brane worlds. We consider the propagation of photon probes on one of the branes, including the response to graviton fluctuations, from both field- and string-theoretical viewpoints. In the latter approach, the dominant source of light-cone broadening may be the recoil of the D-brane, which scales linearly with the string coupling. Astrophysical constraints then place strong restrictions on consistent string models of macroscopic extra dimensions. The broadening we find in the field-theoretical picture seems to be close to the current sensitivity of gravity-wave interferometers, and therefore could perhaps be tested experimentally in the foreseeable future.

We present a microscopic model for light-cone fluctuations ``in vacuo'', which incorporates a treatment of quantum-gravitational recoil effects induced by energetic particles. Treating defects in space-time as solitons in string theory, we derive an energy-dependent refractive index and a stochastic spread in the arrival times of mono-energetic photons due to quantum diffusion through space-time foam, as found previously using an effective Born-Infeld action. Distant astrophysical sources provide sensitive tests of these possible quantum-gravitational phenomena.

A satisfactory theory of quantum gravity may necessitate a drastic modification of our perception of space-time, by giving it a foamy structure at distances comparable to the Planck length. It is argued in this essay that the experimental detection of such structures may be a realistic possibility in the foreseeable future. After a brief review of different theoretical approaches to quantum gravity and the relationships between them, we discuss various possible experimental tests of the quantum nature of space-time. Observations of photons from distant astrophysical sources such as Gamma-Ray Bursters and laboratory experiments on neutral kaon decays may be sensitive to quantum-gravitational effects if they are only minimally suppressed. Experimental limits from the Whipple Observatory and the CPLEAR Collaboration are already probing close to the Planck scale, and significant increases in sensitivity are feasible.

We argue that quantum-gravitational fluctuations in the space-time background give the vacuum non-trivial optical properties that include diffusion and consequent uncertainties in the arrival times of photons, causing stochastic fluctuations in the velocity of light ``in vacuo''. Our proposal is motivated within a Liouville string formulation of quantum gravity that also suggests a frequency-dependent refractive index of the particle vacuum. We construct an explicit realization by treating photon propagation through quantum excitations of $D$-brane fluctuations in the space-time foam. These are described by higher-genus string effects, that lead to stochastic fluctuations in couplings, and hence in the velocity of light. We discuss the possibilities of constraining or measuring photon diffusion ``in vacuo'' via $\gamma$-ray observations of distant astrophysical sources.

We reconsider the possibility that the masses of the three light neutrinos of the Standard Model might be almost degenerate and close to the present upper limits from Tritium beta decay and cosmology. In such a scenario, the cancellations required by the latest upper limit on neutrinoless double-beta decay enforce near-maximal mixing that may be compatible only with the vacuum-oscillation scenario for solar neutrinos. We argue that the mixing angles yielded by degenerate neutrino mass-matrix textures are not in general stable under small perturbations. We evaluate within the MSSM the generation-dependent one-loop renormalization of neutrino mass-matrix textures that yielded degenerate masses and large mixing at the tree level. We find that m_{nu_e} > m_{nu_mu} > m_{nu_tau} after renormalization, excluding MSW effects on solar neutrinos. We verify that bimaximal mixing is not stable, and show that the renormalized masses and mixing angles are not compatible with all the experimental constraints, even for tanbeta as low as unity. These results hold whether the neutrino masses are generated by a see-saw mechanism with heavy neutrinos weighing approx. 10^{13} GeV or by non-renormalizable interactions at a scale approx. 10^5 GeV. We also comment on the corresponding renormalization effects in the minimal Standard Model, in which m_{nu_e} < m_{nu_mu} < m_{nu_tau}. Although a solar MSW effect is now possible, the perturbed neutrino masses and mixings are still not compatible with atmospheric- and solar-neutrino data.

Arguments are given that the lightest supersymmetric particle should be a neutralino $\chi$. Minimizing the fine tuning of the gauge hierarchy favours $\Omega_{\chi} h^2 \sim 0.1$. There are important constraints on the parameter space of the MSSM from the stability of the electroweak vacuum. Co-annihilation with the next-to-lightest supersymmetric particle is potentially significant. Incorporating the latest accelerator constraints from LEP and elsewhere, we find that 50 GeV $\lappeq m_\chi \lappeq$ 600 GeV and $\tan\beta \gappeq$ 2.5, if soft supersymmetry breaking parameters are assumed to be universal.

The Cosmic Microwave Background (CMB) provides a precious window on fundamental physics at very high energy scales, possibly including quantum gravity, GUTs and supersymmetry. The CMB has already enabled defect-based rivals to inflation to be discarded, and will be able to falsify many inflationary models. In combination with other cosmological observations, including those of high-redshift supernovae and large-scale structure, the CMB is on the way to providing a detailed budget for the density of the Universe, to be compared with particle-physics calculations for neutrinos and cold dark matter. Thus CMB measurements complement experiments with the LHC and long-baseline neutrino beams.

We consider the possibility that higher-curvature corrections could drive inflation after the compactification to four dimensions. Assuming that the low-energy limit of the fundamental theory is eleven-dimensional supergravity to the lowest order, including curvature corrections and taking the descent from eleven dimensions to four via an intermediate five-dimensional theory, as favored by recent considerations of unification at some scale around $\sim 10^{16}$ GeV, we may obtain a simple model of inflation in four dimensions. The effective degrees of freedom are two scalar fields and the metric. The scalars arise as the large five-dimensional modulus and the self-interacting conformal mode of the metric. The effective potential has a local maximum in addition to the more usual minimum. However, the potential is quite flat at the top, and admits topological inflation. We show that the model can resolve cosmological problems and provide a mechanism for structure formation with very little fine tuning.

Some particle candidates for dark matter are reviewed in the light of recent experimental and theoretical developments. Models for massive neutrinos are discussed in the light of the recent atmospheric-neutrino data, and used to motivate comments on the plausibility of different solutions to the solar neutrino problem. Arguments are given that the lightest supersymmetric particle should be a neutralino \chi, and accelerator and astrophysical constraints used to suggest that 50 GeV \lappeq m_\chi \lappeq 600 GeV. Minimizing the fine tuning of the gauge hierarchy favours \Omega_\chi h^2 \sim 0.1. The possibility of superheavy relic particles is mentioned, and candidates from string and $M$ theory are reviewed. Finally, the possibility of non-zero vacuum energy is discussed: its calculation is a great opportunity for a quantum theory of gravity, and the possibility that it is time dependent should not be forgotten.

Live $^{60}$Fe has recently been reported in a deep-ocean ferromanganese crust. Analysis of the isotopic ratios in the sample suggests that the measured $^{60}$Fe abundance exceeds the levels generated by terrestrial and cosmogenic sources, and it has been proposed that the excess of $^{60}$Fe is a signature of a supernova that exploded near the earth several Myr ago. In this paper, we consider the possible background sources, and confirm that the measured $^{60}$Fe is significantly higher than all known backgrounds, in contrast with the reported abundance of live $^{53}$Mn. We discuss scenarios in which the data are consistent with a supernova event at a distance $D \sim 30$ pc and an epoch $t_{\rm SN} \sim 5$ Myr ago. We propose tests that could confirm or refute the interpretation of the $^{60}$Fe discovery, including searches for $^{10}$Be, $^{129}$I and $^{146}$Sm. Such a nearby supernova event might have had some impact on the earth's biosphere, principally by enhancing the cosmic-ray flux. This might have damaged the earth's ozone layer, enhancing the penetration of solar ultraviolet radiation. In this connection, we comment on the Middle Miocene and Pliocene mini-extinction events. We also speculate on the possibility of a supernova-induced "cosmic-ray winter," if cosmic rays play a significant role in seeding cloud formation.

We discuss possible signatures of quantum gravity for the propagation of light, including an energy-dependent velocity (refractive index), dispersion in velocity at a given energy, and birefringence. We also compare the sensitivities of different astrophysical observations, including BATSE data on GRB 920229, BeppoSAX data on GRB 980425, the possible HEGRA observation of GRB 920925c, and Whipple observations of the active galaxy Mrk 421. Finally, we discuss the prospective sensitivities of AMS and GLAST.

We consider cosmology in the framework of a `material reference system' of D particles, including the effects of quantum recoil induced by closed-string probe particles. We find a time-dependent contribution to the cosmological vacuum energy, which relaxes to zero as $\sim 1/ t^2$ for large times $t$. If this energy density is dominant, the Universe expands with a scale factor $R(t) \sim t^2$. We show that this possibility is compatible with recent observational constraints from high-redshift supernovae, and may also respect other phenomenological bounds on time variation in the vacuum energy imposed by early cosmology.

We consider the effects of neutralino-stau coannihilations on the cosmological relic density of the lightest supersymmetric particle (LSP) $\chi$ in the minimal supersymmetric extension of the Standard Model (MSSM), particularly in the constrained MSSM in which universal supergravity inputs at the GUT scale are assumed. For much of the parameter space in these models, $\chi$ is approximately a U(1) gaugino ${\tilde B}$, and constraints on the cosmological relic density $\Omega_{\tilde B} h^2$ yield an upper bound on the bino mass. We show that in regions of parameter space for which the cosmological bound is nearly saturated, coannihilations of the $\tilde B$ with the stau, the next lightest sparticle, are important and may reduce significantly the bino relic density. Including also bino coannihilations with the selectrons and smuons, we find that the upper limit on $m_{\chi}$ is increased from about 200 GeV to about 600 GeV in the constrained MSSM, with a similar new upper limit expected in the MSSM.

We analyze the LEP constraints from searches for charginos $\chi^{\pm}$ and neutralinos $\chi_i$, taking into account radiative corrections to the relations between their masses and the underlying Higgs-mixing and gaugino-mass parameters $\mu, m_{1/2}$ and the trilinear mass parameter $A_t$. Whilst radiative corrections do not alter the excluded domain in $m_{\chi^{\pm}}$ as a function of $m_{\chi^{\pm}} - m_{\chi}$, its mapping into the $\mu, m_{1/2}$ plane is altered. We update our previous lower limits on the mass of gaugino dark matter and on tan$\beta$, the ratio of Higgs vacuum expectation values, in the light of the latest LEP data and these radiative corrections. We also discuss the viability of Higgsino dark matter, incorporating co-annihilation effects into the calculation of the Higgsino relic abundance. We find that Higgsino dark matter is viable for only a very limited range of $\mu$ and $m_{1/2}$, which will be explored completely by upcoming LEP runs.

The recent confirmation that at least some gamma-ray bursters (GRBs) are indeed at cosmological distances raises the possibility that observations of these could provide interesting constraints on the fundamental laws of physics. Here we demonstrate that the fine-scale time structure and hard spectra of GRB emissions are very sensitive to the possible dispersion of electromagnetic waves in vacuo with velocity differences $\delta v \sim E/E_{\QG}$, as suggested in some approaches to quantum gravity. A simple estimate shows that GRB measurements might be sensitive to a dispersion scale $E_{QG}$ comparable to the Planck energy scale $E_{P} \sim 10^{19}$ GeV, sufficient to test some of these theories, and we outline aspects of an observational programme that could address this goal.

We reconsider superheavy dark matter candidates in string and M theory, in view of the possibility that inflation might generate superheavy particles with an abundance close to that required for a near-critical Universe. We argue that cryptons - stable or metastable bound states of matter in the hidden sector - are favoured over other possible candidates in string or $M$ theory, such as the Kaluza-Klein states associated with extra dimensions. We exhibit a specific string model that predicts cryptons as hidden-sector bound states weighing $\sim 10^{12}$ GeV, and discuss their astrophysical observability.

A significant lower limit on the mass of the lightest neutralino \chi can be obtained by combining the results from sparticle searches at LEP at centre-of-mass energies up to 172 GeV with cosmological considerations, if it is assumed that the \chi is stable. Exclusion domains from slepton searches close m_\chi \sim 0 loopholes that were left open by previous lower-energy LEP searches for charginos and neutralinos, leading to the lower limit $m_\chi \ga 17$ GeV. The constraints on supersymmetric parameter space are strengthened significantly if LEP constraints on supersymmetric Higgs bosons are taken into account, and further if the relic neutralino density is required to fall within the range favoured by astrophysics and cosmology. These bounds are considerably strengthened if universality at the GUT scale is assumed for soft supersymmetry-breaking scalar masses, including those of the Higgs bosons. In this case, the Higgs searches play a dramatic role, and we find that m_\chi \ga 40 GeV. Furthermore, we find that if tan\beta \la 1.7 for \mu<0, or tan\beta \la 1.4 for \mu >0, the cosmological relic density is too large for all values of m_\chi.

We re-examine the very stringent limits on the axion mass based on the strength and duration of the neutrino signal from SN 1987A, in the light of new measurements of the axial-vector coupling strength of nucleons, possible suppression of axion emission due to many-body effects, and additional emission processes involving pions. The suppression of axion emission due to nucleon spin fluctuations induced by many-body effects degrades previous limits by a factor of about 2. Emission processes involving thermal pions can strengthen the limits by a factor of 3-4 within a perturbative treatment that neglects saturation of nucleon spin fluctuations. Inclusion of saturation effects, however, tends to make the limits less dependent on pion abundances. The resulting axion mass limit also depends on the precise couplings of the axion and ranges from 0.5x10**(-3) eV to 6x10**(-3) eV.

We discuss the lower limit on the mass of the neutralino $\chi$ that can be obtained by combining data from $e^+e^-$ annihilation at LEP and elsewhere with astrophysical and theoretical considerations. Loopholes in the purely experimental analysis of ALEPH data from the Z peak and LEP 1.5, which appear when $\mu<0$ for certain values of the sneutrino mass $m_{\tilde\nu}$ and the ratio $\tan\beta$ of supersymmetric Higgs vacuum expectation values, may be largely or totally excluded by data from lower-energy $e^+e^-$ data, the hypothesis that most of the cosmological dark matter consists of $\chi$ particles, and the assumption that electroweak symmetry breaking is triggered by radiative corrections due to a heavy top quark. The combination of these inputs imposes $m_{\chi} \ge 21.4~\gev$, if soft supersymmetry-breaking masses are assumed to be universal at the grand-unification scale.

Nearby supernova explosions may cause geological isotope anomalies via the direct deposition of debris or by cosmic-ray spallation in the earth's atmosphere. We estimate the mass of material deposited terrestrially by these two mechanisms, showing the dependence on the supernova distance. A number of radioactive isotopes are identified as possible diagnostic tools, such as Be-10, Al-26, Cl-36, Mn-53, Fe-60, and Ni-59, as well as the longer-lived I-129, Sm-146, and Pu-244. We discuss whether the 35 and 60 kyr-old Be-10 anomalies observed in the Vostok antarctic ice cores could be due to supernova explosions. Combining our estimates for matter deposition with results of recent nucleosynthesis yields, we calculate the expected signal from nearby supernovae using ice cores back to $\sim 300$ kyr ago, and we discuss using deep ocean sediments back to several hundred Myr. In particular, we examine the prospects for identifying isotope anomalies due to the Geminga supernova explosion, and signatures of the possibility that supernovae might have caused one or more biological mass extinctions.

Neutrino telescopes of large area offer the possibility of searching for indirect signals of relic neutralinos in the galactic halo, due to annihilations in the Sun or the Earth. Here we investigate the sensitivity, using a supergravity scheme where the soft scalar mass terms are not constrained by universality conditions at the grand unification scale. We first discuss in which regions of the supersymmetric parameter space the neutralino may be considered as a good candidate for cold dark matter. The discovery potential of the search using neutrino telescopes is then compared to that of the direct search for relic neutralinos.

At the request of the organizers, this talk surveys some of the hot topics discussed at this meeting, giving my {\it subjective views} on them. Subjects covered include the present age and Hubble expansion rate of the Universe - {\it inflation theorists need not yet abandon $\Omega = 1$}, theories of structure formation in the light of COBE and other data - {\it my favourite is a flat spectrum of initial perturbations subsequently amplified by mixed hot and cold dark matter}, neutrino masses and oscillations - {\it the only experimental indication I take seriously at the moment is the persistent solar neutrino deficit}, the lightest supersymmetric particle - {\it which may behave differently if conventional assumptions are relaxed}, and the axion - {\it much of the window between limits from SN 1987a and cosmology will be explored in an ongoing experiment}. Finally, I present a chronology of some possible interesting future experiments.

Neutralino dark matter is studied in the context of a supergravity scheme where the scalar mass terms are not constrained by universality conditions at the grand unification scale. We analyse in detail the consequences of the relaxation of this universality assumption on the supersymmetric parameter space, on the neutralino relic abundance and on the event rate for the direct detection of relic neutralinos.

We investigate the light gravitino regeneration rate in the early Universe in models based on $N=1$ supergravity. Motivated by a recent claim by Fischler, we evaluate finite-temperature effects on the gravitino regeneration rate due to the hot primordial plasma for a wide range of the supersymmetry-breaking scale $F$. We find that thermal corrections to the gravitino pole mass and to the Goldstino coupling are negligible for a wide range of temperatures, thereby justifying the extension of the equivalence theorem for the helicity-1/2 gravitino and Goldstino to a hot primordial plasma background. Utilizing the Braaten-Pisarski resummation method, the helicity-1/2 gravitino regeneration rate is found to be $0.25 \alpha_s(T) \log(1/\alpha_s(T))|{m_{\rm soft}/F}|^2 T^3(1 + \alpha_s(T) \log(1/\alpha_s(T)) + T^2 / |F|)$ up to a calculable, model-dependent ${\cal O}(1)$ numerical factor. We review the implications of this regeneration rate for supergravity cosmology, focusing in particular on scenaria for baryogenesis.

We describe a scenario for inflation in the framework of non-critical string theory, which does not employ an inflaton field. There is an exponential expansion of the volume of the Universe, induced by enormous entropy production in the early stages of cosmological evolution. This is associated with the loss of information carried by global string modes that cross the particle horizon. It is the same loss of information that induces irreversible time flow when target time is identified with the world-sheet Liouville mode. The resulting scenario for inflation is described by a string analogue of the Fokker-Planck equation that incorporates diffusion and dissipative effects. Cosmological density perturbations are naturally small.

We review our approach to time and quantum dynamics based on non-critical string theory, developing its relationship to previous work on non-equilibrium quantum statistical mechanics and the microscopic arrow of time. We exhibit specific non-factorizing contributions to the ${\nd S}$ matrix associated with topological defects on the world sheet, explaining the r\^ole that the leakage of ${W_{\infty}}$ charges plays in the loss of quantum coherence. We stress the analogy with the quantum Hall effect, discuss the violation of $CPT$, and also apply our approach to cosmology.

In the non-critical string framework that we have proposed recently, the time $t$ is identified with a dynamical local renormalization group scale, the Liouville mode, and behaves as a statistical evolution parameter, flowing irreversibly from an infrared fixed point - which we conjecture to be a topological string phase - to an ultraviolet one - which corresponds to a static critical string vacuum. When applied to a toy two-dimensional model of space-time singularities, this formalism yields an apparent renormalization of the velocity of light, and a $t$-dependent form of the uncertainty relation for position and momentum of a test string. We speculate within this framework on a stringy alternative to conventional field-theoretical inflation, and the decay towards zero of the cosmological constant in a maximally-symmetric space.

We discuss baryogenesis using the flipped $SU(5)$ model for lepton mass matrices. We show that the generalized see-saw mechanism in this model can not only provide MSW neutrino mixing suitable for solving the solar neutrino problem, and supply a hot dark matter candidate ($\nu_\tau$) with mass $0(10)eV$ as indicated by recent COBE results, but can also naturally account for the baryon asymmetry of the universe. Heavy singlet neutrino decay generates a net lepton asymmetry which is subsequently reprocessed by nonperturbative electroweak interactions. We evaluate the baryon asymmetry so produced in light of the constraints that the COBE observations put on inflationary cosmologies, finding it comfortably consistent with observation.