Toward an accurate mass function for precision cosmology

Cosmological surveys aim to use the evolution of the abundance of galaxy clusters to accurately constrain the cosmological model. In the context of LCDM, we show that it is possible to achieve the required percent level accuracy in the halo mass function with gravity-only cosmological simulations, and we provide simulation start and run parameter guidelines for doing so. Some previous works have had sufficient statistical precision, but lacked robust verification of absolute accuracy. Convergence tests of the mass function with, for example, simulation start redshift can exhibit false convergence of the mass function due to counteracting errors, potentially misleading one to infer overly optimistic estimations of simulation accuracy. Percent level accuracy is possible if initial condition particle mapping uses second order Lagrangian Perturbation Theory, and if the start epoch is between 10 and 50 expansion factors before the epoch of halo formation of interest. The mass function for halos with fewer than ~1000 particles is highly sensitive to simulation parameters and start redshift, implying a practical minimum mass resolution limit due to mass discreteness. The narrow range in converged start redshift suggests that it is not presently possible for a single simulation to capture accurately the cluster mass function while also starting early enough to model accurately the numbers of reionisation era galaxies, whose baryon feedback processes may affect later cluster properties. Ultimately, to fully exploit current and future cosmological surveys will require accurate modeling of baryon physics and observable properties, a formidable challenge for which accurate gravity-only simulations are just an initial step.Comment: revised in response to referee suggestions, MNRAS accepte

Quantifying the heart of darkness with GHALO - a multi-billion particle simulation of our galactic halo

We perform a series of simulations of a Galactic mass dark matter halo at different resolutions, our largest uses over three billion particles and has a mass resolution of 1000 M_sun. We quantify the structural properties of the inner dark matter distribution and study how they depend on numerical resolution. We can measure the density profile to a distance of 120 pc (0.05% of R_vir) where the logarithmic slope is -0.8 and -1.4 at (0.5% of R_vir). We propose a new two parameter fitting function that has a linearly varying logarithmic density gradient which fits the GHALO and VL2 density profiles extremely well. Convergence in the density profile and the halo shape scales as N^(-1/3), but the shape converges at a radius three times larger at which point the halo becomes more spherical due to numerical resolution. The six dimensional phase-space profile is dominated by the presence of the substructures and does not follow a power law, except in the smooth under-resolved inner few kpc.Comment: 6 pages, 4 figures, submitted to MNRAS Letters, for full sized images, see http://www.itp.uzh.ch/news.htm

Relationship between ecosystem productivity and photosynthetically-active radiation for northern peatlands

We analyzed the relationship between net ecosystem exchange of carbon dioxide (NEE) and irradiance (as photosynthetic photon flux density or PPFD), using published and unpublished data that have been collected during midgrowing season for carbon balance studies at seven peatlands in North America and Europe. NEE measurements included both eddy-correlation tower and clear, static chamber methods, which gave very similar results. Data were analyzed by site, as aggregated data sets by peatland type (bog, poor fen, rich fen, and all fens) and as a single aggregated data set for all peatlands. In all cases, a fit with a rectangular hyperbola (NEE = α PPFD Pmax/(α PPFD + Pmax) + R) better described the NEE-PPFD relationship than did a linear fit (NEE = β PPFD + R). Poor and rich fens generally had similar NEE-PPFD relationships, while bogs had lower respiration rates (R = −2.0μmol m−2s−1 for bogs and −2.7 μmol m−2s−1 for fens) and lower NEE at moderate and high light levels (Pmax = 5.2 μmol m−2s−1 for bogs and 10.8 μmol m−2s−1 for fens). As a single class, northern peatlands had much smaller ecosystem respiration (R = −2.4 μmol m−2s−1) and NEE rates (α = 0.020 and Pmax = 9.2μmol m−2s−1) than the upland ecosystems (closed canopy forest, grassland, and cropland) summarized by Ruimy et al. [1995]. Despite this low productivity, northern peatland soil carbon pools are generally 5–50 times larger than upland ecosystems because of slow rates of decomposition caused by litter quality and anaerobic, cold soils

Fossil remnants of reionization in the halo of the Milky Way

Our recently completed one billion particle Via Lactea II simulation of a Milky Way-sized dark matter halo resolves over 50,000 gravitationally bound clumps orbiting today within the virialized region of the main host. About 2,300 of these subhalos have one or more "progenitors" above 1e6 Msun at redshift 11, i.e. massive enough for their gas to have cooled via excitation of H2 and fragmented prior to the epoch of cosmic reionization. We count 4,500 such progenitors: if these were able to convert a fraction of their gas content into very metal-poor stars with a Salpeter initial mass function (IMF), they would be shining today with a visual magnitude M_V=6.7 per solar mass in stars. Assuming a universal baryon fraction, we show that mean star formation efficiencies as low as 0.1% in progenitors below 1e8 Msun would overproduce the abundance of the faint Galatic dwarf spheroidals observed by the Sloan Digital Sky Survey. Star formation at first light must have occurred either with an IMF lacking stars below 0.9 Msun, or was intrinsically very inefficient in small dark matter halos. If the latter, our results may be viewed as another hint of a minimum scale in galaxy formation.Comment: 5 pages, 4 figures, ApJL, in pres

Investigation of initiation of gigantic jets connecting thunderclouds to the ionosphere

The initiation of giant electrical discharges called as "gigantic jets" connecting thunderclouds to the ionosphere is investigated by numerical simulation method in this paper. Using similarity relations, the triggering conditions of streamer formation in laboratory situations are extended to form a criterion of initiation of gigantic jets. The energy source causing a gigantic jet is considered due to the quasi-electrostatic field generated by thunderclouds. The electron dynamics from ionization threshold to streamer initiation are simulated by the Monte Carlo technique. It is found that gigantic jets are initiated at a height of ~18-24 km. This is in agreement with the observations. The method presented in this paper could be also applied to the analysis of the initiation of other discharges such as blue jets and red sprites.Comment: 12th International Congress on Plasma Physics, 25-29 October 2004, Nice (France

Glycan shifting on hepatitis C virus (HCV) E2 glycoprotein is a mechanism for escape from broadly neutralizing antibodies

Hepatitis C virus (HCV) infection is a major cause of liver disease and hepatocellular carcinoma. Glycan shielding has been proposed to be a mechanism by which HCV masks broadly neutralizing epitopes on its viral glycoproteins. However, the role of altered glycosylation in HCV resistance to broadly neutralizing antibodies is not fully understood. Here, we have generated potent HCV neutralizing antibodies hu5B3.v3 and MRCT10.v362 that, similar to the previously described AP33 and HCV1, bind to a highly conserved linear epitope on E2. We utilize a combination of in vitro resistance selections using the cell culture infectious HCV and structural analyses to identify mechanisms of HCV resistance to hu5B3.v3 and MRCT10.v362. Ultra deep sequencing from in vitro HCV resistance selection studies identified resistance mutations at asparagine N417 (N417S, N417T and N417G) as early as 5 days post treatment. Comparison of the glycosylation status of soluble versions of the E2 glycoprotein containing the respective resistance mutations revealed a glycosylation shift from N417 to N415 in the N417S and N417T E2 proteins. The N417G E2 variant was glycosylated neither at residue 415 nor at residue 417 and remained sensitive to MRCT10.v362. Structural analyses of the E2 epitope bound to hu5B3.v3 Fab and MRCT10.v362 Fab using X-ray crystallography confirmed that residue N415 is buried within the antibody–peptide interface. Thus, in addition to previously described mutations at N415 that abrogate the β-hairpin structure of this E2 linear epitope, we identify a second escape mechanism, termed glycan shifting, that decreases the efficacy of broadly neutralizing HCV antibodies

A Study of Time-Dependent CP-Violating Asymmetries and Flavor Oscillations in Neutral B Decays at the Upsilon(4S)

We present a measurement of time-dependent CP-violating asymmetries in neutral B meson decays collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at the Stanford Linear Accelerator Center. The data sample consists of 29.7 ${\rm fb}^{-1}$ recorded at the $\Upsilon(4S)$ resonance and 3.9 ${\rm fb}^{-1}$ off-resonance. One of the neutral B mesons, which are produced in pairs at the $\Upsilon(4S)$, is fully reconstructed in the CP decay modes $J/\psi K^0_S$, $\psi(2S) K^0_S$, $\chi_{c1} K^0_S$, $J/\psi K^{*0}$ ($K^{*0}\to K^0_S\pi^0$) and $J/\psi K^0_L$, or in flavor-eigenstate modes involving $D^{(*)}\pi/\rho/a_1$ and $J/\psi K^{*0}$ ($K^{*0}\to K^+\pi^-$). The flavor of the other neutral B meson is tagged at the time of its decay, mainly with the charge of identified leptons and kaons. The proper time elapsed between the decays is determined by measuring the distance between the decay vertices. A maximum-likelihood fit to this flavor eigenstate sample finds $\Delta m_d = 0.516\pm 0.016 {\rm (stat)} \pm 0.010 {\rm (syst)} {\rm ps}^{-1}$. The value of the asymmetry amplitude $\sin2\beta$ is determined from a simultaneous maximum-likelihood fit to the time-difference distribution of the flavor-eigenstate sample and about 642 tagged $B^0$ decays in the CP-eigenstate modes. We find $\sin2\beta=0.59\pm 0.14 {\rm (stat)} \pm 0.05 {\rm (syst)}$, demonstrating that CP violation exists in the neutral B meson system. (abridged)Comment: 58 pages, 35 figures, submitted to Physical Review

Biotransformation of Silver Released from Nanoparticle Coated Titanium Implants Revealed in Regenerating Bone

Antimicrobial silver nanoparticle coatings have attracted interest for reducing prosthetic joint infection. However, few studies report in vivo investigations of the biotransformation of silver nanoparticles within the regenerating tissue and its impact on bone formation. We present a longitudinal investigation of the osseointegration of silver nanoparticle-coated additive manufactured titanium implants in rat tibial defects. Correlative imaging at different time points using nanoscale secondary ion mass spectrometry, transmission electron microscopy (TEM), histomorphometry, and 3D X-ray microcomputed tomography provided quantitative insight from the nano- to macroscales. The quality and quantity of newly formed bone is comparable between the uncoated and silver coated implants. The newly formed bone demonstrates a trabecular morphology with bone being located at the implant surface, and at a distance, at two weeks. Nanoscale elemental mapping of the bone−implant interface showed that silver was present primarily in the osseous tissue and colocalized with sulfur. TEM revealed silver sulfide nanoparticles in the newly regenerated bone, presenting strong evidence that the previously in vitro observed biotransformation of silver to silver sulfide occurs in vivo

Measurement of the quasi-elastic axial vector mass in neutrino-oxygen interactions

The weak nucleon axial-vector form factor for quasi-elastic interactions is determined using neutrino interaction data from the K2K Scintillating Fiber detector in the neutrino beam at KEK. More than 12,000 events are analyzed, of which half are charged-current quasi-elastic interactions nu-mu n to mu- p occurring primarily in oxygen nuclei. We use a relativistic Fermi gas model for oxygen and assume the form factor is approximately a dipole with one parameter, the axial vector mass M_A, and fit to the shape of the distribution of the square of the momentum transfer from the nucleon to the nucleus. Our best fit result for M_A = 1.20 \pm 0.12 GeV. Furthermore, this analysis includes updated vector form factors from recent electron scattering experiments and a discussion of the effects of the nucleon momentum on the shape of the fitted distributions.Comment: 14 pages, 10 figures, 6 table