Parallel beam dynamics calculations on high performance computers
Faced with a backlog of nuclear waste and weapons plutonium, as well as an ever-increasing public concern about safety and environmental issues associated with conventional nuclear reactors, many countries are studying new, accelerator-driven technologies that hold the promise of providing safe and effective solutions to these problems. Proposed projects include accelerator transmutation of waste (ATW), accelerator-based conversion of plutonium (ABC), accelerator-driven energy production (ADEP), and accelerator production of tritium (APT). Also, next-generation spallation neutron sources based on similar technology will play a major role in materials science and biological science research. The design of accelerators for these projects will require a major advance in numerical modeling capability. For example, beam dynamics simulations with approximately 100 million particles will be needed to ensure that extremely stringent beam loss requirements (less than a nanoampere per meter) can be met. Compared wi...

title: COMPUTATIONAL CHALLENGES IN HIGH INTENSITY ION BEAM PHYSICS

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: 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: 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: 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: Measurement of isolated photon production in deep-inelastic scattering at HERA
abstract: The production of isolated photons in deep-inelastic scattering ep→eγX is measured with the H1 detector at HERA. The measurement is performed in the kinematic range of negative four-momentum transfer squared 4<Q2<150 GeV2 and a mass of the hadronic system WX>50 GeV. The analysis is based on a total integrated luminosity of 227 pb-1. The production cross section of isolated photons with a transverse energy in the range 3<ETγ<10 GeV and pseudorapidity range -1.2<ηγ<1.8 is measured as a function of ETγ, ηγ and Q2. Isolated photon cross sections are also measured for events with no jets or at least one hadronic jet. The measurements are compared with predictions from Monte Carlo generators modelling the photon radiation from the quark and the electron lines, as well as with calculations at leading and next to leading order in the strong coupling. The predictions significantly underestimate the measured cross sections.

title: Tests of QCD factorisation in the diffractive production of dijets in deep-inelastic scattering and photoproduction at HERA
abstract: Measurements are presented of differential dijet cross sections in diffractive photoproduction (Q2<0.01 GeV2) and deep-inelastic scattering processes (DIS, 4<Q2<80 GeV2). The event topology is given by ep→eXY, in which the system X, containing at least two jets, is separated from a leading low-mass baryonic system Y by a large rapidity gap. The dijet cross sections are compared with NLO QCD predictions based on diffractive parton densities previously obtained from a QCD analysis of inclusive diffractive DIS cross sections by H1. In DIS, the dijet data are well described, supporting the validity of QCD factorisation. The diffractive DIS dijet data are more sensitive to the diffractive gluon density at high fractional parton momentum than the measurements of inclusive diffractive DIS. In photoproduction, the predicted dijet cross section has to be multiplied by a factor of approximately 0.5 for both direct and resolved photon interactions to describe the measurements. The ratio of measured dijet cross section to NLO prediction in photoproduction is a factor 0.5±0.1 smaller than the same ratio in DIS. This suppression is the first clear observation of QCD hard scattering factorisation breaking at HERA. The measurements are also compared to the two soft colour neutralisation models SCI and GAL. The SCI model describes diffractive dijet production in DIS but not in photoproduction. The GAL model fails in both kinematic regions.

title: Measurement of inclusive jet production in deep-inelastic scattering at high Q2 and determination of the strong coupling
abstract: Inclusive jet production is studied in neutral current deep-inelastic positron-proton scattering at large four momentum transfer squared Q^2>150 GeV^2 with the H1 detector at HERA. Single and double differential inclusive jet cross sections are measured as a function of Q^2 and of the transverse energy E_T of the jets in the Breit frame. The measurements are found to be well described by calculations at next-to-leading order in perturbative QCD. The running of the strong coupling is demonstrated and the value of alpha_s(M_Z) is determined. The ratio of the inclusive jet cross section to the inclusive neutral current cross section is also measured and used to extract a precise value for alpha_s(M_Z)=0.1193+/-0.0014(exp.)^{+0.0047}_{-0.0030}(th.)+/-0.0016(pdf).

title: Search for lepton flavour violation in ep collisions at HERA
abstract: A search for the lepton flavour violating processes ep→μX and ep→τX is performed with the H1 experiment at HERA. Final states with a muon or tau and a hadronic jet are searched for in a data sample corresponding to an integrated luminosity of 66.5 pb-1 for e+p collisions and 13.7 pb-1 for e-p collisions at a centre-of-mass energy of 319 GeV. No evidence for lepton flavour violation is found. Limits are derived on the mass and the couplings of leptoquarks inducing lepton flavour violation in an extension of the Buchmüller–Rückl–Wyler effective model. Leptoquarks produced in ep collisions with a coupling strength of λ=0.3 and decaying with the same coupling strength to a muon–quark pair or a tau–quark pair are excluded at 95% confidence level up to masses of 459 GeV and 379 GeV, respectively.