Workshop report on large-scale matrix diagonalization methods in chemistry theory institute

The Large-Scale Matrix Diagonalization Methods in Chemistry theory institute brought together 41 computational chemists and numerical analysts. The goal was to understand the needs of the computational chemistry community in problems that utilize matrix diagonalization techniques. This was accomplished by reviewing the current state of the art and looking toward future directions in matrix diagonalization techniques. This institute occurred about 20 years after a related meeting of similar size. During those 20 years the Davidson method continued to dominate the problem of finding a few extremal eigenvalues for many computational chemistry problems. Work on non-diagonally dominant and non-Hermitian problems as well as parallel computing has also brought new methods to bear. The changes and similarities in problems and methods over the past two decades offered an interesting viewpoint for the success in this area. One important area covered by the talks was overviews of the source and nature of the chemistry problems. The numerical analysts were uniformly grateful for the efforts to convey a better understanding of the problems and issues faced in computational chemistry. An important outcome was an understanding of the wide range of eigenproblems encountered in computational chemistry. The workshop covered problems involving self- consistent-field (SCF), configuration interaction (CI), intramolecular vibrational relaxation (IVR), and scattering problems. In atomic structure calculations using the Hartree-Fock method (SCF), the symmetric matrices can range from order hundreds to thousands. These matrices often include large clusters of eigenvalues which can be as much as 25% of the spectrum. However, if Cl methods are also used, the matrix size can be between 10{sup 4} and 10{sup 9} where only one or a few extremal eigenvalues and eigenvectors are needed. Working with very large matrices has lead to the development o

Testing Logselfsimilarity of Soil Particle Size Distribution: Simulation with Minimum Inputs

Particle size distribution (PSD) greatly influences other soil physical properties. A detailed textural analysis is time-consuming and expensive. Soil texture is commonly reported in terms of mass percentages of a small number of size fractions (typically, clay, silt and sand). A method to simulate the PSD from such a poor description or even from the poorest description, consisting in the mass percentages of only two soil size fractions, would be extremly useful for prediction purposes. The goal of this paper is to simulate soil PSDs from the minimum number of inputs, i.e., two and three textural fraction contents, by using a logselfsimilar model and an iterated function system constructed with these data. High quality data on 171 soils are used. Additionally, the characterization of soil texture by entropy-based parameters provided by the model is tested. Results indicate that the logselfsimilar model may be a useful tool to simulate PSD for the construction of pedotransfer functions related to other soil properties when textural information is limited to moderate textural data

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

W boson polarization measurement in the ttbar dilepton channel using the CDF II Detector

We present a measurement of $W$ boson polarization in top-quark decays in $t\bar{t}$ events with decays to dilepton final states using $5.1 {\rm fb^{-1}}$ of integrated luminosity in $p\bar{p}$ collisions collected by the CDF II detector at the Tevatron. A simultaneous measurement of the fractions of longitudinal ($f_0$) and right-handed ($f_+$) $W$ bosons yields the results $f_0 = 0.71 ^{+0.18}_{-0.17} {\rm (stat)} \pm 0.06 {\rm (syst)}$ and $f_+ = -0.07 \pm 0.09 {\rm (stat)} \pm 0.03 {\rm (syst)}$. Combining this measurement with our previous result based on single lepton final states, we obtain $f_0 = 0.84 \pm 0.09 {\rm (stat)} \pm 0.05 {\rm (syst)}$ and $f_{+} = -0.16 \pm 0.05 {\rm (stat)} \pm 0.04 {\rm (syst)}$. The results are consistent with standard model expectation.Comment: Published in Phys. Lett.

Measurement of the Bs Lifetime in Fully and Partially Reconstructed Bs -> Ds- (phi pi-)X Decays in pbar-p Collisions at sqrt(s) = 1.96 TeV

We present a measurement of the Bs lifetime in fully and partially reconstructed Bs -> Ds(phi pi)X decays in 1.3 fb-1 of pbar-p collisions at sqrt(s) = 1.96 TeV collected by the CDF II detector at the Fermilab Tevatron. We measure tau(Bs) = 1.518 +/- 0.041 (stat.) +/- 0.027 (syst.) ps. The ratio of this result and the world average B0 lifetime yields tau(Bs)/tau(B0) = 0.99 +/-0.03, which is in agreement with recent theoretical predictions.Comment: submitted to Phys. Rev. Let