Fast computation of Bernoulli, Tangent and Secant numbers

We consider the computation of Bernoulli, Tangent (zag), and Secant (zig or Euler) numbers. In particular, we give asymptotically fast algorithms for computing the first n such numbers in O(n^2.(log n)^(2+o(1))) bit-operations. We also give very short in-place algorithms for computing the first n Tangent or Secant numbers in O(n^2) integer operations. These algorithms are extremely simple, and fast for moderate values of n. They are faster and use less space than the algorithms of Atkinson (for Tangent and Secant numbers) and Akiyama and Tanigawa (for Bernoulli numbers).Comment: 16 pages. To appear in Computational and Analytical Mathematics (associated with the May 2011 workshop in honour of Jonathan Borwein's 60th birthday). For further information, see http://maths.anu.edu.au/~brent/pub/pub242.htm

Metallic liquid hydrogen and likely Al2O3 metallic glass

Dynamic compression has been used to synthesize liquid metallic hydrogen at 140 GPa (1.4 million bar) and experimental data and theory predict Al2O3 might be a metallic glass at ~300 GPa. The mechanism of metallization in both cases is probably a Mott-like transition. The strength of sapphire causes shock dissipation to be split differently in the strong solid and soft fluid. Once the 4.5-eV H-H and Al-O bonds are broken at sufficiently high pressures in liquid H2 and in sapphire (single-crystal Al2O3), electrons are delocalized, which leads to formation of energy bands in fluid H and probably in amorphous Al2O3. The high strength of sapphire causes shock dissipation to be absorbed primarily in entropy up to ~400 GPa, which also causes the 300-K isotherm and Hugoniot to be virtually coincident in this pressure range. Above ~400 GPa shock dissipation must go primarily into temperature, which is observed experimentally as a rapid increase in shock pressure above ~400 GPa. The metallization of glassy Al2O3, if verified, is expected to be general in strong oxide insulators. Implications for Super Earths are discussed.Comment: 8 pages, 5 figures, 14th Liquid and Amorphous Metals Conference, Rome 201

Shrinking a large dataset to identify variables associated with increased risk of Plasmodium falciparum infection in Western Kenya

Large datasets are often not amenable to analysis using traditional single-step approaches. Here, our general objective was to apply imputation techniques, principal component analysis (PCA), elastic net and generalized linear models to a large dataset in a systematic approach to extract the most meaningful predictors for a health outcome. We extracted predictors for Plasmodium falciparum infection, from a large covariate dataset while facing limited numbers of observations, using data from the People, Animals, and their Zoonoses (PAZ) project to demonstrate these techniques: data collected from 415 homesteads in western Kenya, contained over 1500 variables that describe the health, environment, and social factors of the humans, livestock, and the homesteads in which they reside. The wide, sparse dataset was simplified to 42 predictors of P. falciparum malaria infection and wealth rankings were produced for all homesteads. The 42 predictors make biological sense and are supported by previous studies. This systematic data-mining approach we used would make many large datasets more manageable and informative for decision-making processes and health policy prioritization

Global Search for New Physics with 2.0/fb at CDF

Data collected in Run II of the Fermilab Tevatron are searched for indications of new electroweak-scale physics. Rather than focusing on particular new physics scenarios, CDF data are analyzed for discrepancies with the standard model prediction. A model-independent approach (Vista) considers gross features of the data, and is sensitive to new large cross-section physics. Further sensitivity to new physics is provided by two additional algorithms: a Bump Hunter searches invariant mass distributions for "bumps" that could indicate resonant production of new particles; and the Sleuth procedure scans for data excesses at large summed transverse momentum. This combined global search for new physics in 2.0/fb of ppbar collisions at sqrt(s)=1.96 TeV reveals no indication of physics beyond the standard model.Comment: 8 pages, 7 figures. Final version which appeared in Physical Review D Rapid Communication

Observation of Orbitally Excited B_s Mesons

We report the first observation of two narrow resonances consistent with states of orbitally excited (L=1) B_s mesons using 1 fb^{-1} of ppbar collisions at sqrt{s} = 1.96 TeV collected with the CDF II detector at the Fermilab Tevatron. We use two-body decays into K^- and B^+ mesons reconstructed as B^+ \to J/\psi K^+, J/\psi \to \mu^+ \mu^- or B^+ \to \bar{D}^0 \pi^+, \bar{D}^0 \to K^+ \pi^-. We deduce the masses of the two states to be m(B_{s1}) = 5829.4 +- 0.7 MeV/c^2 and m(B_{s2}^*) = 5839.7 +- 0.7 MeV/c^2.Comment: Version accepted and published by Phys. Rev. Let

Measurement of the Bottom-Strange Meson Mixing Phase in the Full CDF Data Set

We report a measurement of the bottom-strange meson mixing phase \beta_s using the time evolution of B0_s -> J/\psi (->\mu+\mu-) \phi (-> K+ K-) decays in which the quark-flavor content of the bottom-strange meson is identified at production. This measurement uses the full data set of proton-antiproton collisions at sqrt(s)= 1.96 TeV collected by the Collider Detector experiment at the Fermilab Tevatron, corresponding to 9.6 fb-1 of integrated luminosity. We report confidence regions in the two-dimensional space of \beta_s and the B0_s decay-width difference \Delta\Gamma_s, and measure \beta_s in [-\pi/2, -1.51] U [-0.06, 0.30] U [1.26, \pi/2] at the 68% confidence level, in agreement with the standard model expectation. Assuming the standard model value of \beta_s, we also determine \Delta\Gamma_s = 0.068 +- 0.026 (stat) +- 0.009 (syst) ps-1 and the mean B0_s lifetime, \tau_s = 1.528 +- 0.019 (stat) +- 0.009 (syst) ps, which are consistent and competitive with determinations by other experiments.Comment: 8 pages, 2 figures, Phys. Rev. Lett 109, 171802 (2012

Reparative cell therapy for the heart: critical internal appraisal of the field in response to recent controversies

The concept that cell-based repair of myocardial injury might be possible was introduced almost two decades ago; however, the field of cardiovascular reparative medicine has been criticized as translation to clinically effective approaches has been slow. The recent retraction of a series of papers has further impacted perception of this area of research. As researchers, clinicians, and teachers in this field, we felt it incumbent to critically appraise the current state of cardiac cell repair, determine what can be learned from past mistakes, and formulate best practices for future work. This special communication summarizes an introspective assessment of what has fallen short, how to prevent similar issues, and how the field might best move forward in the service of science and patients.Cardiolog

Search for direct production of charginos and neutralinos in events with three leptons and missing transverse momentum in √s = 7 TeV pp collisions with the ATLAS detector

A search for the direct production of charginos and neutralinos in final states with three electrons or muons and missing transverse momentum is presented. The analysis is based on 4.7 fb−1 of proton–proton collision data delivered by the Large Hadron Collider and recorded with the ATLAS detector. Observations are consistent with Standard Model expectations in three signal regions that are either depleted or enriched in Z-boson decays. Upper limits at 95% confidence level are set in R-parity conserving phenomenological minimal supersymmetric models and in simplified models, significantly extending previous results

Jet size dependence of single jet suppression in lead-lead collisions at sqrt(s(NN)) = 2.76 TeV with the ATLAS detector at the LHC

Measurements of inclusive jet suppression in heavy ion collisions at the LHC provide direct sensitivity to the physics of jet quenching. In a sample of lead-lead collisions at sqrt(s) = 2.76 TeV corresponding to an integrated luminosity of approximately 7 inverse microbarns, ATLAS has measured jets with a calorimeter over the pseudorapidity interval |eta| < 2.1 and over the transverse momentum range 38 < pT < 210 GeV. Jets were reconstructed using the anti-kt algorithm with values for the distance parameter that determines the nominal jet radius of R = 0.2, 0.3, 0.4 and 0.5. The centrality dependence of the jet yield is characterized by the jet "central-to-peripheral ratio," Rcp. Jet production is found to be suppressed by approximately a factor of two in the 10% most central collisions relative to peripheral collisions. Rcp varies smoothly with centrality as characterized by the number of participating nucleons. The observed suppression is only weakly dependent on jet radius and transverse momentum. These results provide the first direct measurement of inclusive jet suppression in heavy ion collisions and complement previous measurements of dijet transverse energy imbalance at the LHC.Comment: 15 pages plus author list (30 pages total), 8 figures, 2 tables, submitted to Physics Letters B. All figures including auxiliary figures are available at http://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/HION-2011-02

CASA: An Efficient Automated Assignment of Protein Mainchain NMR Data Using an Ordered Tree Search Algorithm

Rapid analysis of protein structure, interaction, and dynamics requires fast and automated assignments of 3D protein backbone triple-resonance NMR spectra. We introduce a new depth-first ordered tree search method of automated assignment, CASA, which uses hand-edited peak-pick lists of a flexible number of triple resonance experiments. The computer program was tested on 13 artificially simulated peak lists for proteins up to 723 residues, as well as on the experimental data for four proteins. Under reasonable tolerances, it generated assignments that correspond to the ones reported in the literature within a few minutes of CPU time. The program was also tested on the proteins analyzed by other methods, with both simulated and experimental peaklists, and it could generate good assignments in all relevant cases. The robustness was further tested under various situations.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/43050/1/10858_2005_Article_4079.pd