ROBUSTNESS OF COSMOLOGICAL SIMULATIONS. I. LARGE-SCALE STRUCTURE
The gravitationally driven evolution of cold dark matter dominates the formation of structure in the universe over a wide range of length scales. While the longest scales can be treated by perturbation theory, a fully quantitative understanding of nonlinear effects requires the application of large-scale particle simulation methods. Additionally, precision predictions for next-generation observatio-ns, such as weak gravitational lensing, can only be obtained from numerical simulations. In this paper, we compare results from several N-body codes using test problems and a diverse set of diagnostics, focusing on a medium-resolution regime appropriate for studying many observationally relevant aspects of structure formation. Our conclusions are that—despite the use of different algorithms and error-control methodologies—overall, the codes yield consistent results. The agreement over a wide range of scales for the cosmological tests is test-dependent. In the best cases, it is at the 5% level or better, however, for other cases it can be significantly larger than 10%. These include the halo mass function at low masses and the mass power spectrum at small scales. While there exist explanations for most of the discrepancies, our results point to the need for significant improvement in N-body errors and their understanding to match the precision of near-future observations. The simulation results, including halo catalogs, and initial conditions used, are publicly available.

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: COMPUTATIONAL CHALLENGES IN HIGH INTENSITY ION BEAM PHYSICS

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: 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: Independent measurement of the total active 8B solar neutrino flux using an array of 3He proportional counters at the Sudbury Neutrino Observatory.
abstract: The Sudbury Neutrino Observatory (SNO) used an array of 3He proportional counters to measure the rate of neutral-current interactions in heavy water and precisely determined the total active (nu_x) 8B solar neutrino flux. This technique is independent of previous methods employed by SNO. The total flux is found to be 5.54_-0.31;+0.33(stat)-0.34+0.36(syst)x10(6) cm(-2) s(-1), in agreement with previous measurements and standard solar models. A global analysis of solar and reactor neutrino results yields Deltam2=7.59_-0.21;+0.19x10(-5) eV2 and theta=34.4_-1.2;+1.3 degrees. The uncertainty on the mixing angle has been reduced from SNO's previous results.

title: Search for dark matter and large extra dimensions in monojet events in pp collisions at $ \sqrt {s} = {7} $ TeV
abstract: A bstractA search has been made for events containing an energetic jet and an imbalance in transverse momentum using a data sample of pp collisions at a center-of-mass energy of 7 TeV. This signature is common to both dark matter and extra dimensions models. The data were collected by the CMS detector at the LHC and correspond to an integrated luminosity of 5.0 fb−1. The number of observed events is consistent with the standard model expectation. Constraints on the dark matter-nucleon scattering cross sections are determined for both spin-independent and spin-dependent interaction models. For the spin-independent model, these are the most constraining limits for a dark matter particle with mass below 3.5 GeV/c2, a region unexplored by direct detection experiments. For the spin-dependent model, these are the most stringent constraints over the 0.1–200 GeV/c2 mass range. The constraints on the Arkani-Hamed, Dimopoulos, and Dvali model parameter MD determined as a function of the number of extra dimensions are also an improvement over the previous results.

title: Measurement of inclusive and dijet D∗ meson cross sections in photoproduction at HERA
abstract: The inclusive photoproduction of D∗ mesons and of D∗-tagged dijets is investigated with the H1 detector at the ep collider HERA. The kinematic region covers small photon virtualities Q2<2 GeV2 and photon–proton centre-of-mass energies of 100<Wγp<285 GeV. Inclusive D∗ meson differential cross sections are measured for central rapidities |η(D∗)|<1.5 and transverse momenta pT(D∗)>1.8 GeV. The heavy quark production process is further investigated in events with at least two jets with transverse momentum pT(jet)>3.5 GeV each, one containing the D∗ meson. Differential cross sections for D∗-tagged dijet production and for correlations between the jets are measured in the range |η(D∗)|<1.5 and pT(D∗)>2.1 GeV. The results are compared with predictions from Monte Carlo simulations and next-to-leading order perturbative QCD calculations.

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: Inclusive deep inelastic scattering at high Q2 with longitudinally polarised lepton beams at HERA
abstract: A BSTRACTInclusive e±p single and double differential cross sections for neutral and charged current deep inelastic scattering processes are measured with the H1 detector at HERA. The data were taken at a centre-of-mass energy of $ \sqrt {s} = {319} $ GeV with a total integrated luminosity of 333.7 pb−1 shared between two lepton beam charges and two longitudinal lepton polarisation modes. The differential cross sections are measured in the range of negative four-momentum transfer squared, Q2, between 60 and 50 000 GeV2, and Bjorken x between 0.0008 and 0.65. The measurements are combined with earlier published unpolarised H1 data to improve statistical precision and used to determine the structure function $ xF_{3}^{{\gamma Z}} $. Ameasurementoftheneutralcurrentparityviolating structure function $ F_{2}^{{\gamma Z}} $ is presented for the first time. The polarisation dependence of the charged current total cross section is also measured. The new measurements are well described by a next-to-leading order QCD fit based on all published H1 inclusive cross section data which are used to extract the parton distribution functions of the proton.