Two-dimensional flux-corrected transport solver for convectively dominated flows

A numerical technique designed to solve a wide class of convectively dominated flow problems is presented. An attractive feature of the technique is its ability to resolve the behavior of field quantities possessing large gradients and/or shocks. The method is a finite-difference technique known as flux-corrected transport (FCT) that maintains four important numerical considerations - stability, accuracy, monotonicity, and conservation. The theory and methodology of two-dimensional FCT is presented. The method is applied in demonstrative example calculations of a 2-D Riemann problem with known exact solutions and to the Euler equations in a study of classical Rayleigh-Taylor and Kelvin-Helmholtz instability problems. The FCT solver has been vectorized for execution on the Cray 1S - a typical call with a 50 by 50 mesh requires about 0.00428 cpu seconds of execution time per call to the routine. Additionally, we have maintained a modular structure for the solver that eases its implementation. Fortran listings of two versions of the 2-D FCT solvers are appended with a driver main program illustrating the call sequence for the modules. 59 refs., 49 figs

Probing mSUGRA via the Extreme Universe Space Observatory

An analysis is carried out within mSUGRA of the estimated number of events originating from upward moving ultra-high energy neutralinos that could be detected by the Extreme Universe Space Observatory (EUSO). The analysis exploits a recently proposed technique that differentiates ultra-high energy neutralinos from ultra-high energy neutrinos using their different absorption lengths in the Earth's crust. It is shown that for a significant part of the parameter space, where the neutralino is mostly a Bino and with squark mass $\sim 1$ TeV, EUSO could see ultra-high energy neutralino events with essentially no background. In the energy range 10^9 GeV < E < 10^11 GeV, the unprecedented aperture of EUSO makes the telescope sensitive to neutralino fluxes as low as 1.1 \times 10^{-6} (E/GeV)^{-1.3} GeV^{-1} cm^{-2} yr^{-1} sr^{-1}, at the 95% CL. Such a hard spectrum is characteristic of supermassive particles' $N$-body hadronic decay. The case in which the flux of ultra-high energy neutralinos is produced via decay of metastable heavy particles with uniform distribution throughout the universe is analyzed in detail. The normalization of the ratio of the relics' density to their lifetime has been fixed so that the baryon flux produced in the supermassive particle decays contributes to about 1/3 of the events reported by the AGASA Collaboration below 10^{11} GeV, and hence the associated GeV gamma-ray flux is in complete agreement with EGRET data. For this particular case, EUSO will collect between 4 and 5 neutralino events (with 0.3 of background) in ~ 3 yr of running. NASA's planned mission, the Orbiting Wide-angle Light-collectors (OWL), is also briefly discussed in this context.Comment: Some discussion added, final version to be published in Physical Review

Magnetic Field Effects on Neutron Diffraction in the Antiferromagnetic Phase of UPt3UPt_3

We discuss possible magnetic structures in UPt$_3$ based on our analysis of elastic neutron-scattering experiments in high magnetic fields at temperatures $T<T_N$. The existing experimental data can be explained by a single-{\bf q} antiferromagnetic structure with three independent domains. For modest in-plane spin-orbit interactions, the Zeeman coupling between the antiferromagnetic order parameter and the magnetic field induces a rotation of the magnetic moments, but not an adjustment of the propagation vector of the magnetic order. A triple-{\bf q} magnetic structure is also consistent with neutron experiments, but in general leads to a non-uniform magnetization in the crystal. New experiments could decide between these structures.Comment: 5 figures included in the tex

Spectral method for the time-dependent Gross-Pitaevskii equation with a harmonic trap

We study the numerical resolution of the time-dependent Gross-Pitaevskii equation, a non-linear Schroedinger equation used to simulate the dynamics of Bose-Einstein condensates. Considering condensates trapped in harmonic potentials, we present an efficient algorithm by making use of a spectral Galerkin method, using a basis set of harmonic oscillator functions, and the Gauss-Hermite quadrature. We apply this algorithm to the simulation of condensate breathing and scissors modes.Comment: 23 pages, 5 figure

Top-squark searches at the Tevatron in models of low-energy supersymmetry breaking

We study the production and decays of top squarks (stops) at the Tevatron collider in models of low-energy supersymmetry breaking. We consider the case where the lightest Standard Model (SM) superpartner is a light neutralino that predominantly decays into a photon and a light gravitino. Considering the lighter stop to be the next-to-lightest Standard Model superpartner, we analyze stop signatures associated with jets, photons and missing energy, which lead to signals naturally larger than the associated SM backgrounds. We consider both 2-body and 3-body decays of the top squarks and show that the reach of the Tevatron can be significantly larger than that expected within either the standard supergravity models or models of low-energy supersymmetry breaking in which the stop is the lightest SM superpartner. For a modest projection of the final Tevatron luminosity, L = 4 fb-1, stop masses of order 300 GeV are accessible at the Tevatron collider in both 2-body and 3-body decay modes. We also consider the production and decay of ten degenerate squarks that are the supersymmetric partners of the five light quarks. In this case we find that common squark masses up to 360 GeV are easily accessible at the Tevatron collider, and that the reach increases further if the gluino is light.Comment: 32 pages, 9 figures; references adde

CDMS, Supersymmetry and Extra Dimensions

The CDMS experiment aims to directly detect massive, cold dark matter particles originating from the Milky Way halo. Charge and lattice excitations are detected after a particle scatters in a Ge or Si crystal kept at ~30 mK, allowing to separate nuclear recoils from the dominating electromagnetic background. The operation of 12 detectors in the Soudan mine for 75 live days in 2004 delivered no evidence for a signal, yielding stringent limits on dark matter candidates from supersymmetry and universal extra dimensions. Thirty Ge and Si detectors are presently installed in the Soudan cryostat, and operating at base temperature. The run scheduled to start in 2006 is expected to yield a one order of magnitude increase in dark matter sensitivity.Comment: To be published in the proceedings of the 7th UCLA symposium on sources and detection of dark matter and dark energy in the universe, Marina del Rey, Feb 22-24, 200

The fully differential single-top-quark cross section in next-to-leading order QCD

We present a new next-to-leading order calculation for fully differential single-top-quark final states. The calculation is performed using phase space slicing and dipole subtraction methods. The results of the methods are found to be in agreement. The dipole subtraction method calculation retains the full spin dependence of the final state particles. We show a few numerical results to illustrate the utility and consistency of the resulting computer implementations.Comment: 37 pages, latex, 2 ps figure

A Search for Selectrons and Squarks at HERA

Data from electron-proton collisions at a center-of-mass energy of 300 GeV are used for a search for selectrons and squarks within the framework of the minimal supersymmetric model. The decays of selectrons and squarks into the lightest supersymmetric particle lead to final states with an electron and hadrons accompanied by large missing energy and transverse momentum. No signal is found and new bounds on the existence of these particles are derived. At 95% confidence level the excluded region extends to 65 GeV for selectron and squark masses, and to 40 GeV for the mass of the lightest supersymmetric particle.Comment: 13 pages, latex, 6 Figure

The New Look pMSSM with Neutralino and Gravitino LSPs

The pMSSM provides a broad perspective on SUSY phenomenology. In this paper we generate two new, very large, sets of pMSSM models with sparticle masses extending up to 4 TeV, where the lightest supersymmetric particle (LSP) is either a neutralino or gravitino. The existence of a gravitino LSP necessitates a detailed study of its cosmological effects and we find that Big Bang Nucleosynthesis places strong constraints on this scenario. Both sets are subjected to a global set of theoretical, observational and experimental constraints resulting in a sample of \sim 225k viable models for each LSP type. The characteristics of these two model sets are briefly compared. We confront the neutralino LSP model set with searches for SUSY at the 7 TeV LHC using both the missing (MET) and non-missing ET ATLAS analyses. In the MET case, we employ Monte Carlo estimates of the ratios of the SM backgrounds at 7 and 8 TeV to rescale the 7 TeV data-driven ATLAS backgrounds to 8 TeV. This allows us to determine the pMSSM parameter space coverage for this collision energy. We find that an integrated luminosity of \sim 5-20 fb^{-1} at 8 TeV would yield a substantial increase in this coverage compared to that at 7 TeV and can probe roughly half of the model set. If the pMSSM is not discovered during the 8 TeV run, then our model set will be essentially void of gluinos and lightest first and second generation squarks that are \lesssim 700-800 GeV, which is much less than the analogous mSUGRA bound. Finally, we demonstrate that non-MET SUSY searches continue to play an important role in exploring the pMSSM parameter space. These two pMSSM model sets can be used as the basis for investigations for years to come.Comment: 54 pages, 22 figures; typos fixed, references adde