Partial wave treatment of Supersymmetric Dark Matter in the presence of CP - violation

We present an improved partial wave analysis of the dominant LSP annihilation channel to a fermion-antifermion pair which avoids the non-relativistic expansion being therefore applicable near thresholds and poles. The method we develop allows of contributions of any partial wave in the total angular momentum J in contrast to partial wave analyses in terms of the orbital angular momentum L of the initial state, which is usually truncated to p-waves, and yields very accurate results. The method is formulated in such a way as to allow easy handling of CP-violating phases residing in supersymmetric parameters. We apply this refined partial wave technique in order to calculate the neutralino relic density in the constrained MSSM (CMSSM) in the presence of CP-violating terms occurring in the Higgs - mixing parameter \mu and trilinear A coupling for large tanb. The inclusion of CP-violating phases in mu and A does not upset significantly the picture and the annihilation of the LSP's to a b b_bar, through Higgs exchange, is still the dominant mechanism in obtaining cosmologically acceptable neutralino relic densities in regions far from the stau-coannihilation and the `focus point'. Significant changes can occur if we allow for phases in the gaugino masses and in particular the gluino mass.Comment: 23 pages LaTeX, 10 eps figures, version to appear in PR

Exact Cross Sections for the Neutralino-Slepton Coannihilation

Coannihilation processes provide an important additional mechanism for reducing the density of stable relics in the Universe. In the case of the stable lightest neutralino of the MSSM, and in particular the Constrained MSSM (CMSSM), the coannihilation with sleptons plays a major role in opening up otherwise cosmologically excluded ranges of supersymmetric parameters. In this paper, we derive a full set of exact, analytic expressions for the coannihilation of the lightest neutralino with the sleptons into all two--body tree--level final states in the framework of minimal supersymmetry. We make no simplifying assumptions about the neutralino nor about sfermion masses and mixings other than the absence of explicit CP--violating terms and inter--family mixings. The expressions should be particularly useful in computing the neutralino WIMP relic abundance without the approximation of partial wave expansion. We illustrate the effect of our analytic results with numerical examples and demonstrate a sizeable difference with approximate expressions available in the literature.Comment: LaTeX, 46 pages, 8 eps figure

Gravitino Dark Matter in the CMSSM and Implications for Leptogenesis and the LHC

In the framework of the CMSSM we study the gravitino as the lightest supersymmetric particle and the dominant component of cold dark matter in the Universe. We include both a thermal contribution to its relic abundance from scatterings in the plasma and a non--thermal one from neutralino or stau decays after freeze--out. In general both contributions can be important, although in different regions of the parameter space. We further include constraints from BBN on electromagnetic and hadronic showers, from the CMB blackbody spectrum and from collider and non--collider SUSY searches. The region where the neutralino is the next--to--lightest superpartner is severely constrained by a conservative bound from excessive electromagnetic showers and probably basically excluded by the bound from hadronic showers, while the stau case remains mostly allowed. In both regions the constraint from CMB is often important or even dominant. In the stau case, for the assumed reasonable ranges of soft SUSY breaking parameters, we find regions where the gravitino abundance is in agreement with the range inferred from CMB studies, provided that, in many cases, a reheating temperature \treh is large, \treh\sim10^{9}\gev. On the other side, we find an upper bound \treh\lsim 5\times 10^{9}\gev. Less conservative bounds from BBN or an improvement in measuring the CMB spectrum would provide a dramatic squeeze on the whole scenario, in particular it would strongly disfavor the largest values of \treh\sim 10^{9}\gev. The regions favored by the gravitino dark matter scenario are very different from standard regions corresponding to the neutralino dark matter, and will be partly probed at the LHC.Comment: JHEP version, several improvements and update

Upper and Lower Limits on Neutralino WIMP Mass and Spin--Independent Scattering Cross Section, and Impact of New (g-2)_{mu} Measurement

We derive the allowed ranges of the spin--independent interaction cross section \sigsip for the elastic scattering of neutralinos on proton for wide ranges of parameters of the general Minimal Supersymmetric Standard Model. We investigate the effects of the lower limits on Higgs and superpartner masses from colliders, as well as the impact of constraints from \bsgamma and the new measurement of \gmtwo on the upper and lower limits on \sigsip. We further explore the impact of the neutralino relic density, including coannihilation, and of theoretical assumptions about the largest allowed values of the supersymmetric parameters. For $\mu>0$, requiring the latter to lie below 1\tev leads to \sigsip\gsim 10^{-11}\pb at \mchi\sim100\gev and \sigsip\gsim 10^{-8}\pb at \mchi\sim1\tev. When the supersymmetric parameters are allowed above 1\tev, for 440\gev \lsim \mchi\lsim 1020 \gev we derive a {\em parameter--independent lower limit} of \sigsip \gsim 2\times 10^{-12}\pb. (No similar lower limits can be set for $\mu<0$ nor for 1020\gev\lsim\mchi\lsim2.6\tev.) Requiring \abundchi<0.3 implies a {\em parameter--independent upper limit} \mchi\lsim2.6\tev. The new \epem--based measurement of $(g-2)_{\mu}$ restricts \mchi\lsim 350\gev at $1 \sigma$ CL and \mchi\lsim515\gev at $2 \sigma$ CL, and implies $\mu>0$. The largest allowed values of \sigsip have already become accessible to recent experimental searches.Comment: LaTeX, 17 pages, 9 eps figures. Version to appear in JHE

Effect of River Levee with Geosynthetic-Reinforced Soil against Overflow Erosion and Infiltration

Overflows from huge floods have caused levee breaches in a great number of places, including Japan. To prevent such destruction and thereby increase the resistance of armored levees to overflow erosion, in this study, we examined the performances of Geosynthetic-Reinforced Soil (GRS) levees against overflow erosion under various conditions, such as reinforcement, back slopes, and geo-grid layers. In addition, we investigated the effect of geo-grid layers on the infiltration of levees. The model tests revealed that 1) with scour protection in front of the toe of the back slope, the GRS levee exhibits much higher resistance against overflow erosion than the armored levee; 2) the armored levee with a steep back slope (= 1:0.5) collapsed faster than that with a normal slope (= 1:2). However, the GRS levee with a steep back slope of 1:0.5 maintained high resistance against overflow erosion after the target time. 3) The GRS levee with partial and full reinforcements had a comparably high resistance against overflow erosion. 4) The GRS levee using a small-sized geo-grid maintained a high residual ratio of the cross-sectional area over a long period. 5) The infiltration discharge of the GRS levee was lesser than that of the levee with no reinforcement due to the reduction in infiltration erosion in the GRS levee. These facts suggest that the GRS levee with partial reinforcement can be applied to the reinforcement of existing levees, and appropriately sized geo-grid layers should be selected

Exact Cross Sections for the Neutralino WIMP Pair-Annihilation

We derive a full set of exact, analytic expressions for the annihilation of the lightest neutralino pairs into all two-body tree-level final states in the framework of minimal supersymmetry. We make no simplifying assumptions about the neutralino nor about sfermion masses and mixings other than the absence of explicit CP--violating terms. The expressions should be particularly useful in computing the neutralino WIMP relic abundance without the usual approximation of partial wave expansion.Comment: LaTeX, 46 pages, no figures. Several minor typographical errors correcte

Computation of 3D Frequency-Domain Waveform Kernals for c(x,y,z) Media

Seismic tomography, as typically practiced on both the exploration, crustal, and global scales, considers only the arrival times of selected sets of phases and relies primarily on WKBJ theory during inversion. Since the mid 1980’s, researchers have explored, largely on a theoretical level, the possibility of inverting the entire seismic record. Due to the ongoing advances in CPU performance, full waveform inversion is finally becoming feasible on select problems with promising results emerging from frequency-domain methods. However, frequency-domain techniques using sparse direct solvers are currently constrained by memory limitations in 3D where they exhibit a O(n4) worst-case bound on memory usage. We sidestep this limitation by using a hybrid approach, calculating frequency domain Green’s functions for the scalar wave equation by driving a high-order, time-domain, finite-difference (FDTD) code to steady state using a periodic source. The frequency-domain response is extracted using the phase sensitive detection (PSD) method recently developed by Nihei and Li (2006). The resulting algorithm has an O(n3) memory footprint and is amenable to parallelization in the space, shot, or frequency domains. We demonstrate this approach by generating waveform inversion kernels for fully c(x,y,z) models. Our test examples include a realistic VSP experiment using the geometry and velocity models obtained from a site in Western Wyoming, and a deep crustal reflection/refraction profile based on the LARSE II geometry and the SCEC community velocity model. We believe that our 3D solutions to the scalar Helmholtz equation, for models with upwards of 100 million degrees of freedom, are the largest examples documented in the open geophysical literature. Such results suggest that iterative 3D waveform inversion is an achievable goal in the near future.Shell GameChangerMassachusetts Institute of Technology. Earth Resources Laborator

Non-minimal Split Supersymmetry

We present an extension of the minimal split supersymmetry model, which is capable of explaining the baryon asymmetry of the Universe. Instead of MSSM we start from NMSSM and split its spectrum in such a way that the low energy theory contains neutral particles, in addition to the content of minimal split supersymmetry. They trigger the strongly first order electroweak phase transition (EWPT) and provide an additional source of CP-violation. In this model, we estimate the amount of the baryon asymmetry produced during EWPT, using WKB approximation for CP-violating sources in diffusion equations. We also examine the contribution of CP-violating interactions to the electron and neutron electric dipole moments and estimate the production of the neutralino dark matter. We find that both phenomenological and cosmological requirements can be fulfilled in this model.Comment: 31 pages, 9 figures, typos correcte

Inflation and Gauge Hierarchy in Randall-Sundrum Compactification

We obtain the general inflationary solutions for the slab of five-dimensional AdS spacetime where the fifth dimension is an orbifold $S^1/Z_2$ and two three-branes reside at its boundaries, of which the Randall-Sundrum model corresponds to the static limit. The investigation of the general solutions and their static limit reveals that the RS model recasts both the cosmological constant problem and the gauge hierarchy problem into the balancing problem of the bulk and the brane cosmological constants.Comment: 9 pages, revtex, minor changes and more references adde

Lattice QCD calculation of the proton decay matrix element in the continuum limit

We present a quenched lattice QCD calculation of the \alpha and \beta parameters of the proton decay matrix element. The simulation is carried out using the Wilson quark action at three values of the lattice spacing in the range a\approx 0.1-0.064 fm to study the scaling violation effect. We find only mild scaling violation when the lattice scale is determined by the nucleon mass. We obtain in the continuum limit, |\alpha(NDR,2GeV)|=0.0090(09)(^{+5}_{-19})GeV^3 and |\beta(NDR,2GeV)|=0.0096(09)(^{+6}_{-20})GeV^3 with \alpha and \beta in a relatively opposite sign, where the first error is statistical and the second is due to the uncertainty in the determination of the physical scale.Comment: 4 pages, 3 figure