Beam dynamics calculations and particle tracking using massively parallel processors
During the past decade massively parallel processors (MPPs) have slowly gained acceptance within the scientific community. At present these machines typically contain a few hundred to one thousand off-the-shelf microprocessors and a total memory of up to 32 GBytes. The potential performance of these machines is illustrated by the fact that a month long job on a high end workstation might require only a few hours on an MPP. The acceptance of MPPs has been slow for a variety of reasons. For example, some algorithms are not easily parallelizable. Also, in the past these machines were difficult to program. But in recent years the development of Fortran-like languages such as CM Fortran and High Performance Fortran have made MPPs much easier to use. In the following we will describe how MPPs can be used for beam dynamics calculations and long term particle tracking.

title: CHAOS IN TIME-DEPENDENT VARIATIONAL APPROXIMATIONS TO QUANTUM DYNAMICS
abstract: Dynamical chaos has recently been shown to exist in the Gaussian approximation in quan-tum mechanics and in the self-consistent mean ﬁeld approach to studying the dynamics ofquantum ﬁelds. In this study, we ﬁrst show that any variationalapproximationto the dynam-ics of a quantum system based on the Dirac action principle leads to a classical Hamiltoniandynamics for the variational parameters. Since this Hamiltonian is generically nonlinearand nonintegrable, the dynamics thus generated can be chaotic, in distinction to the exactquantum evolution. We then restrict attention to a system of two biquadratically coupledquantum oscillators and study two variational schemes, the leading order large N (fourcanonical variables) and Hartree (six canonical variables) approximations. The chaos seenin the approximate dynamics is an artifact of the approximations: this is demonstrated bythe fact that its onset occurs on the same characteristic time scale as the breakdown of theapproximations when compared to numerical solutions of the time-dependent Schr¨odingerequation.PACS numbers: 05.45. +b, 03.65. Sq, 2.30 Wd,

title: Joint Efficient Dark-energy Investigation (JEDI): a Candidate Implementation of the NASA-DOE Joint Dark Energy Mission (JDEM)
abstract: JEDI will probe dark energy in three independent ways by measuring the expansion history of the universe: (1) using type Ia supernovae as cosmological standard candles over a range of distances, (2) using baryon oscillations as a cosmological standard ruler over a range of cosmic epochs, (3) mapping the weak gravitational lensing distortion by foreground galaxies of the images of background galaxies at different distances. JEDI will unravel the nature of dark energy with accuracy and precision. 
JEDI is a 2m-class space telescope with the unique ability of simultaneous wide-field imaging (0.8-4.2 micron in five bands) and multi-slit spectroscopy (0.8-3.2 micron) with a field of view of 1 square degree. What makes JEDI efficient is its ability to simultaneously obtain high signal-to-noise ratio, moderate resolution slit spectra for all supernovae and ~ 5000 galaxies in its wide field of view, and to combine imaging and spectroscopy so that the appropriate balance of time is devoted to each. 
JEDI will measure the cosmic expansion history H(z) as a free function to 10^9 galaxies to H \~ 25.

title: COSMIC EMULATION: THE CONCENTRATION-MASS RELATION FOR wCDM UNIVERSES
abstract: The concentration-mass relation for dark matter-dominated halos is one of the essential results expected from a theory of structure formation. We present a simple prediction scheme, a cosmic emulator, for the concentration-mass (c-M) relation as a function of cosmological parameters for wCDM models. The emulator is constructed from 37 individual models, with three nested N-body gravity-only simulations carried out for each model. The mass range covered by the emulator is 2 × 1012 M ☉ < M < 1015 M ☉ with a corresponding redshift range of z = 0-1. Over this range of mass and redshift, as well as the variation of cosmological parameters studied, the mean halo concentration varies from c ~ 2 to c ~ 8. The distribution of the concentration at fixed mass is Gaussian with a standard deviation of one-third of the mean value, almost independent of cosmology, mass, and redshift over the ranges probed by the simulations. We compare results from the emulator with previously derived heuristic analytic fits for the c-M relation, finding that they underestimate the halo concentration at high masses. Using the emulator to investigate the cosmology dependence of the c-M relation over the currently allowable range of values, we find—not surprisingly—that σ8 and ω m influence it considerably, but also that the dark energy equation-of-state parameter w has a substantial effect. In general, the concentration of lower-mass halos is more sensitive to changes in cosmological parameters as compared to cluster mass halos. The c-M emulator is publicly available from http://www.hep.anl.gov/cosmology/CosmicEmu.

title: COMPUTATIONAL CHALLENGES IN HIGH INTENSITY ION BEAM PHYSICS

title: DARK MATTER HALO PROFILES OF MASSIVE CLUSTERS: THEORY VERSUS OBSERVATIONS
abstract: Dark-matter-dominated cluster-scale halos act as an important cosmological probe and provide a key testing ground for structure formation theory. Focusing on their mass profiles, we have carried out (gravity-only) simulations of the concordance ΛCDM cosmology, covering a mass range of 2 × 1012 to 2 × 1015 h –1 M ☉ and a redshift range of z = 0-2, while satisfying the associated requirements of resolution and statistical control. When fitting to the Navarro-Frenk-White profile, our concentration-mass (c-M) relation differs in normalization and shape in comparison to previous studies that have limited statistics in the upper end of the mass range. We show that the flattening of the c-M relation with redshift is naturally expressed if c is viewed as a function of the peak height parameter, ν. Unlike the c-M relation, the slope of the c-ν relation is effectively constant over the redshift range z = 0-2, while the amplitude varies by ~30% for massive clusters. This relation is, however, not universal: using a simulation suite covering the allowed wCDM parameter space, we show that the c-ν relation varies by about ±20% as cosmological parameters are varied. At fixed mass, the c(M) distribution is well fit by a Gaussian with σ c /c 1/3, independent of the radius at which the concentration is defined, the halo dynamical state, and the underlying cosmology. We compare the ΛCDM predictions with observations of halo concentrations from strong lensing, weak lensing, galaxy kinematics, and X-ray data, finding good agreement for massive clusters (M vir > 4 × 1014 h –1 M ☉), but with some disagreements at lower masses. Because of uncertainty in observational systematics and modeling of baryonic physics, the significance of these discrepancies remains unclear.

title: Measurement of beauty photoproduction near threshold using di-electron events with the H1 detector at HERA
abstract: The cross section for $ep \rightarrow e\, b\bar{b} X$ in photoproduction is measured with the H1 detector at the ep-collider HERA. The decay channel $b\bar{b} \rightarrow ee X'$ is selected by identifying the semi-electronic decays of the b-quarks. The total production cross section is measured in the kinematic range given by the photon virtuality Q2≤1 GeV2, the inelasticity 0.05≤y≤0.65 and the pseudorapidity of the b-quarks $|\eta(b)|, |\eta(\bar{b})|\leq2$. The differential production cross section is measured as a function of the average transverse momentum of the beauty quarks 〈PT(b)〉 down to the threshold. The results are compared to next-to-leading-order QCD predictions.

title: Charged particle production in high Q 2 deep-inelastic scattering at HERA
abstract: The average charged track multiplicity and the normalised distribution of the scaled momentum, $\xp$, of charged final state hadrons are measured in deep-inelastic $\ep$ scattering at high $Q^2$ in the Breit frame of reference. The analysis covers the range of photon virtuality $100 < Q^2 < 20 000 \GeV^{2}$. Compared with previous results presented by HERA experiments this analysis has a significantly higher statistical precision and extends the phase space to higher $Q^{2}$ and to the full range of $\xp$. The results are compared with $e^+e^-$ annihilation data and with various calculations based on perturbative QCD using different models of the hadronisation process.

title: Search for heavy long-lived charged particles in $pp$ collisions at $s=7$ TeV
abstract: The result of a search for heavy long-lived charged particles produced in pp collisions at sqrt(s) = 7 TeV at the LHC is described. The data sample has been collected using the CMS detector and corresponds to an integrated luminosity of 5.0 inverse femtobarns. The inner tracking detectors are used to define a sample of events containing tracks with high momentum and high ionization energy loss. A second sample of events with high-momentum and high-ionization tracks satisfying muon identification and long time-of-flight criteria is analyzed independently. In both cases, the results are consistent with the expected background estimated from data. The results are used to establish cross section limits as a function of mass within the context of models with long-lived gluinos, scalar top quarks and scalar taus. Cross section limits on hyper-meson particles, containing new elementary long-lived hyper-quarks predicted by a vector-like confinement model, are also presented. Lower limits at 95% confidence level on the mass of gluinos (scalar top quarks) are found to be 1098 (737) GeV. A limit of 928 (626) GeV is set for a gluino (scalar top quark) that hadronizes into a neutral bound state before reaching the muon detectors. The lower mass limit for a pair produced scalar tau is found to be 223 GeV. Mass limits for a hyper-kaon are placed at 484, 602, and 747 GeV for hyper-rho masses of 800, 1200, and 1600 GeV, respectively.

title: Diffractive open charm production in deep-inelastic scattering and photoproduction at HERA
abstract: Measurements are presented of diffractive open charm production at HERA. The event topology is given by ep→eXY where the system X contains at least one charmed hadron and is well separated by a large rapidity gap from a leading low-mass proton remnant system Y. Two analysis techniques are used for the cross section measurements. In the first, the charm quark is tagged by the reconstruction of a D*±(2010) meson. This technique is used in deep-inelastic scattering (DIS) and photoproduction (γp). In the second, a method based on the displacement of tracks from the primary vertex is used to measure the open charm contribution to the inclusive diffractive cross section in DIS. The measurements are compared with next-to-leading order QCD predictions based on diffractive parton density functions previously obtained from a QCD analysis of the inclusive diffractive cross section at H1. A good agreement is observed in the full kinematic regime, which supports the validity of QCD factorization for open charm production in diffractive DIS and γp.

title: Search in leptonic channels for heavy resonances decaying to long-lived neutral particles
abstract: A bstractA search is performed for heavy resonances decaying to two long-lived massive neutral particles, each decaying to leptons. The experimental signature is a distinctive topology consisting of a pair of oppositely charged leptons originating at a separated secondary vertex. Events were collected by the CMS detector at the LHC during pp collisions at $ \sqrt{s}=7 $ TeV, and selected from data samples corresponding to 4.1 (5.1) fb−1 of integrated luminosity in the electron (muon) channel. No significant excess is observed above standard model expectations, and an upper limit is set with 95% confidence level on the production cross section times the branching fraction to leptons, as a function of the long-lived massive neutral particle lifetime.