SciDAC advances and applications in computational beam dynamics
SciDAC has had a major impact on computational beam dynamics and the design of particle accelerators. Particle accelerators -- which account for half of the facilities in the DOE Office of Science Facilities for the Future of Science 20 Year Outlook -- are crucial for US scientific, industrial, and economic competitiveness. Thanks to SciDAC, accelerator design calculations that were once thought impossible are now carried routinely, and new challenging and important calculations are within reach. SciDAC accelerator modeling codes are being used to get the most science out of existing facilities, to produce optimal designs for future facilities, and to explore advanced accelerator concepts that may hold the key to qualitatively new ways of accelerating charged particle beams. In this poster we present highlights from the SciDAC Accelerator Science and Technology (AST) project Beam Dynamics focus area in regard to algorithm development, software development, and applications.

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: 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: IMPACT: An Integrated-Map & Particle Accelerator Tracking Code for Linac Simulation

title: Wigner Functions and Density Matrices in Curved Spaces as Computational Tools
abstract: This paper contrasts two alternative approaches to statistical quantum field theory in curved spacetimes, namely (1) a canonical Hamiltonian approach, in which the basic object is a density matrix ϱ characterizing the noncovariant, but globally defined, modes of the field; and (2) a Wigner function approach, in which the basic object is a Wigner function f defined quasilocally from the Hadamard, or correlation, function G 1 ( x 1 , x 2 ). The key object is to isolate on the conceptual biases underlying each of these approaches and then to assess their utility and limitations in effecting concrete calculations. The following questions are therefore addressed and largely answered. What sorts of spacetimes (e.g., de Sitter or Friedmann-Robertson-Walker) are comparatively easy to consider? What sorts of objects (e.g., average fields or renormalized stress energies) are easy to compute approximately? What, if anything, can be computed exactly? What approximations are intrinsic to each approach or convenient as computational tools? What sorts of “field entropies” are natural to define?

title: Quantum Fields in Curved Spaces: Kinetic Theory.

title: Search for squarks in R-parity violating supersymmetry in ep collisions at HERA
abstract: A search for squarks in R-parity violating supersymmetry is performed in e±p collisions at HERA using the H1 detector. The full data sample taken at a centre-of-mass energy $\sqrt{s}=319$ GeV is used for the analysis, corresponding to an integrated luminosity of 255 pb−1 of e+p and 183 pb−1 of  collision data. The resonant production of squarks via a Yukawa coupling λ′ is considered, taking into account direct and indirect R-parity violating decay modes. Final states with jets and leptons are investigated. No evidence for squark production is found and mass dependent limits on λ′ are obtained in the framework of the Minimal Supersymmetric Standard Model and in the Minimal Supergravity Model. In the considered part of the parameter space, for a Yukawa coupling of electromagnetic strength λ′=0.3, squarks of all flavours are excluded up to masses of 275 GeV at 95% confidence level, with down-type squarks further excluded up to masses of 290 GeV.

title: Neutral and charged current cross sections and extraction of structure functions and parton distributions

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 baryonic resonances decaying to Ξπ in deep-inelastic scattering at HERA
abstract: A search for narrow baryonic resonances decaying into Ξ-π- or Ξ-π+ and their antiparticles is carried out with the H1 detector using deep inelastic scattering events at HERA in the range of negative photon four-momentum transfer squared 2 < Q2 < 100 GeV2. No signal is observed for a new baryonic state in the mass range 1600–2300 MeV in either the doubly charged or the neutral decay channels. The known Ξ(1530)0 is observed through its decay mode into Ξ-π+. Upper limits are given on the ratio of the production rates of new baryonic states, such as the hypothetical pentaquark states Ξ–5q or Ξ05q, relative to the Ξ(1530)0 baryon state.

title: Measurement of dijet production in diffractive deep-inelastic scattering with a leading proton at HERA
abstract: The cross section of diffractive deep-inelastic scattering ep→eXp is measured, where the system X contains at least two jets and the leading final state proton is detected in the H1 Forward Proton Spectrometer. The measurement is performed for fractional proton longitudinal momentum loss xℙ<0.1 and covers the range 0.1<|t|<0.7 GeV2 in squared four-momentum transfer at the proton vertex and 4<Q2<110 GeV2 in photon virtuality. The differential cross sections extrapolated to |t|<1 GeV2 are in agreement with next-to-leading order QCD predictions based on diffractive parton distribution functions extracted from measurements of inclusive and dijet cross sections in diffractive deep-inelastic scattering. The data are also compared with leading order Monte Carlo models.