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 perform a detailed study of an effective field theory which includes the Standard Model particle content extended by a pair of Weyl fermionic SU(2)-doublets with opposite hypercharges. A discrete symmetry guarantees that a linear combination of the doublet components is stable and can act as a candidate particle for Dark Matter. The dark sector fermions interact with the Higgs and gauge bosons through renormalizable $d=4$ operators, and non-renormalizable $d=5$ operators that appear after integrating out extra degrees of freedom above the TeV scale. We study collider, cosmological and astrophysical probes for this effective theory of Dark Matter. We find that a WIMP with a mass nearby to the electroweak scale, and thus observable at LHC, is consistent with collider and astrophysical data only when fairly large magnetic dipole moment transition operators with the gauge bosons exist, together with moderate Yukawa interactions.

In this talk, we present the package GravitinoPack that calculates decays of unstable supersymmetric particles, involving gravitinos in the final or initial state. If the gravitino is the dark matter particle and therefore stable, the package calculates the decays of the lightest neutralino, and the lighter stau or stop NLSP into the gravitino LSP and one or two Standard Model particles. On the other hand, assuming that the gravitino is unstable, GravitinoPack calculates all its two-body and the three-body decay widths to the neutralino LSP and Standard Model particles. Since all these decays, involving the gravitino, are of gravitational nature, the lifetime of the decaying particle can be of the order of seconds are more, hence called "late decays". The precise knowledge of all these partial decay widths enables the user to apply accurately the relevant cosmological constraints.

We discuss challenges to a dark matter annihilation interpretation of the excess positron fraction in the cosmic rays observed by the PAMELA, Fermi-LAT and AMS-02 collaborations. The spectra of positrons from annihilations into the leptonic two-body final states like $\tau^+ \tau^-$ or $\mu^+ \mu^-$, fit well the AMS-02 data for the positron ratio and the electron flux. Furthermore, we discuss the hadronic annihilation channels $b \bar{b}$ and $W^+ W^-$. However, this interpretation requires a very large annihilation cross section especially for the hadronic channels, conflicting with the unitarity upper limit, unless the positron flux due to annihilations is boosted by a large factor due to inhomogeneities in the galactic halo. In addition, we present predictions within this interpretation for the positron fraction at higher energies for the antiproton flux, and discuss constraints from $\gamma$-ray measurements.

We present results from a new calculation of two- and three-body decays of the gravitino including all possible Feynman graphs. The work is done in the R-parity conserving Minimal Supersymmetric Standard Model, assuming that the gravitino is unstable. We calculate the full width for the three-body decay process Gravitino -> \chi_1 X Y, where \chi_1 is the lightest neutralino that plays the role of the lightest supersymmetric particle, and X, Y are standard model particles. For calculating the squared amplitudes we use the packages FeynArts and FormCalc, with appropriate extensions to accommodate the gravitino with spin 3/2. Our code treats automatically the intermediate exchanged particles that are on their mass shell, using the narrow width approximation. This method considerably simplifies the complexity of the evaluation of the total gravitino width and stabilizes it numerically. At various points of the supersymmetric parameter space, we compare the full two-body and the three-body results with approximations used in the past. In addition, we discuss the relative sizes of the full two-body and the three-body gravitino decay widths. Our results are obtained using the package GravitinoPack, that will be publicly available soon, including various handy features.

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 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.

It has been pointed out recently that the presence of dilaton field in the early Universe can dilute the neutralino dark matter (DM) abundance, if Universe is not radiation dominated at DM decoupling, due to its dissipative-like coupling to DM. In this scenario two basic mechanisms compete, the modified Hubble expansion rate tending to increase the relic density and a dissipative force that tends to decrease it. The net effect can lead to an overall dramatic decrease of the predicted relic abundance, sometimes by amounts of the order of O(10^2) or so. This feature is rather generic, independent of any particular assumption on the underlying string dynamics, provided dilaton dominates at early eras after the end of inflation but before Big Bang Nucleosynthesis (BBN). The latter ensures that BBN is not upset by the presence of the dilaton. In this paper, within the context of such a scenario, we study the phenomenology of the constrained minimal supersymmetric model (CMSSM) by taking into account all recent experimental constraints, including those from the LHC searches. We find that the allowed parameter space is greatly enlarged and includes regions that are beyond the reach of LHC. The allowed regions are compatible with Direct Dark Matter searches since the small neutralino annihilation rates, that are now in accord with the cosmological data on the relic density, imply small neutralino-nucleon cross sections below the sensitivities of the Direct Dark Matter experiments. It is also important that the new cosmologically accepted regions are compatible with Higgs boson masses larger than 120 GeV, as it is indicated from the LHC experimental data. The smaller annihilation cross sections needed to explain WMAP data require that the detector performances of current and planned indirect DM search experiments through gamma rays should be greatly improved in order to probe the CMSSM regions.

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).

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 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 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.

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.

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.

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.

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.

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'.

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.

The recent data from E821 Brookhaven experiment in conjunction with a new determination of the hadronic vacuum polarization contribution to the anomalous magnetic moment of muon, put new bounds on the parameters of the Constrained Minimal Supersymmetric Standard Model. We study the impact this experimental information, along with the $b \goes s \gamma$ branching ratio and light Higgs boson mass bound from LEP, to constrain regions of the model which are consistent with the cosmological data. The effect of these to Dark Matter direct searches is also discussed.

In the framework of the Constrained Minimal Supersymmetric Standard Model we discuss the impact of the recent experimental information, especially from the E821 Brookhaven experiment on $g_{\mu}-2$ along with the light Higgs boson mass bound from LEP, in delineating regions of the parameters which are consistent with cosmological data. The effect of these to the Dark Matter direct searches is also discussed.

In SUGRA inspired supersymmetric models with universal boundary conditions for the soft masses, the scalar cross section $\sigma_{scalar}$ for the elastic neutralino--nucleon scattering is in general several orders of magnitude below the sensitivity of current experiments. For large $\tan \beta$ and low $M_{1/2}, m_0$ values, the theoretically predicted $\sigma_{scalar}$ can approach the sensitivity of these experiments ($\approx 10^{-6} pb$) being at the same time in agreement with recent cosmological data, which impose severe restrictions on the CDM relic density, and with accelerator experiments which put lower bounds on sparticle and Higgs boson masses. Further improvement of the sensitivity of DAMA and CDMS experiments will probe the large $\tan \beta$ region of the parameter space in the vicinity of the boundaries of the parameter space allowed by chargino and Higgs searches.

We discuss the relic density of the lightest of the supersymmetric particles in view of new cosmological data, which favour the concept of an accelerating Universe with a non-vanishing cosmological constant. Recent astrophysical observations provide us with very precise values of the relevant cosmological parameters. Certain of these parameters have direct implications on particle physics, e.g., the value of matter density, which in conjunction with electroweak precision data put severe constraints on the supersymmetry breaking scale. In the context of the Constrained Minimal Supersymmetric Standard Model (CMSSM) such limits read as: $M_{1/2} \simeq 300 \GeV - 340 \GeV$, $m_0 \simeq 80 \GeV - 130 \GeV$. Within the context of the CMSSM a way to avoid these constraints is either to go to the large $\tan \beta$ and $\mu > 0$ region, or make ${\tilde \tau}_R$, the next to lightest supersymmetric particle (LSP), be almost degenerate in mass with LSP.

We discuss the relic density of the lightest of the supersymmetric particles ({\small LSP}) in view of new cosmological data, which favour the concept of an accelerating Universe with a non-vanishing cosmological constant. The new bound on the Cold Dark Matter density, $\Omega_{\mathrm{CDM}} h_0^2 \lesssim 0.22$, puts stringent constraints on supersymmetry preferring low supersymmetry breaking scales, in sharp contrast to electroweak precision measurements favouring large supersymmetry breaking scales. Supersymmetric predictions are in agreement with cosmological data and electroweak precision data in the window of the parameter space: $m_0<200\GeV$, $300\GeV<M_{1/2}<400\GeV$, putting bounds on sparticle masses, which may be evaded if $m_{LSP} <m_{{\tilde{\tau}}_R} \lesssim 1.2 m_{LSP}$.