Scale-invariant moving finite elements for nonlinear partial differential equations in two dimensions

A scale-invariant moving finite element method is proposed for the adaptive solution of nonlinear partial differential equations. The mesh movement is based on a finite element discretisation of a scale-invariant conservation principle incorporating a monitor function, while the time discretisation of the resulting system of ordinary differential equations is carried out using a scale-invariant time-stepping which yields uniform local accuracy in time. The accuracy and reliability of the algorithm are successfully tested against exact self-similar solutions where available, and otherwise against a state-of-the-art h-refinement scheme for solutions of a two-dimensional porous medium equation problem with a moving boundary. The monitor functions used are the dependent variable and a monitor related to the surface area of the solution manifold

Projective and Coarse Projective Integration for Problems with Continuous Symmetries

Temporal integration of equations possessing continuous symmetries (e.g. systems with translational invariance associated with traveling solutions and scale invariance associated with self-similar solutions) in a ``co-evolving'' frame (i.e. a frame which is co-traveling, co-collapsing or co-exploding with the evolving solution) leads to improved accuracy because of the smaller time derivative in the new spatial frame. The slower time behavior permits the use of {\it projective} and {\it coarse projective} integration with longer projective steps in the computation of the time evolution of partial differential equations and multiscale systems, respectively. These methods are also demonstrated to be effective for systems which only approximately or asymptotically possess continuous symmetries. The ideas of projective integration in a co-evolving frame are illustrated on the one-dimensional, translationally invariant Nagumo partial differential equation (PDE). A corresponding kinetic Monte Carlo model, motivated from the Nagumo kinetics, is used to illustrate the coarse-grained method. A simple, one-dimensional diffusion problem is used to illustrate the scale invariant case. The efficiency of projective integration in the co-evolving frame for both the macroscopic diffusion PDE and for a random-walker particle based model is again demonstrated

A moving mesh finite element algorithm for the adaptive solution of time-dependent partial differential equations with moving boundaries

A moving mesh finite element algorithm is proposed for the adaptive solution of nonlinear diffusion equations with moving boundaries in one and two dimensions. The moving mesh equations are based upon conserving a local proportion, within each patch of finite elements, of the total “mass” that is present in the projected initial data. The accuracy of the algorithm is carefully assessed through quantitative comparison with known similarity solutions, and its robustness is tested on more general problems. Applications are shown to a variety of problems involving time-dependent partial differential equations with moving boundaries. Problems which conserve mass, such as the porous medium equation and a fourth order nonlinear diffusion problem, can be treated by a simplified form of the method, while problems which do not conserve mass require the full theory

Neutrino Quasielastic Scattering on Nuclear Targets: Parametrizing Transverse Enhancement (Meson Exchange Currents)

We present a parametrization of the observed enhancement in the transverse electron quasielastic (QE) response function for nucleons bound in carbon as a function of the square of the four momentum transfer ($Q^2$) in terms of a correction to the magnetic form factors of bound nucleons. The parametrization should also be applicable to the transverse cross section in neutrino scattering. If the transverse enhancement originates from meson exchange currents (MEC), then it is theoretically expected that any enhancement in the longitudinal or axial contributions is small. We present the predictions of the "Transverse Enhancement" model (which is based on electron scattering data only) for the $\nu_\mu, \bar{\nu}_\mu$ differential and total QE cross sections for nucleons bound in carbon. The $Q^2$ dependence of the transverse enhancement is observed to resolve much of the long standing discrepancy in the QE total cross sections and differential distributions between low energy and high energy neutrino experiments on nuclear targets.Comment: Revised Version- July 21, 2011: 17 pages, 20 Figures. To be published in Eur. Phys. J.

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

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

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

Search for the standard model Higgs boson decaying to a bbˉb\bar{b} pair in events with no charged leptons and large missing transverse energy using the full CDF data set

We report on a search for the standard model Higgs boson produced in association with a vector boson in the full data set of proton-antiproton collisions at $\sqrt{s} = 1.96$ TeV recorded by the CDF II detector at the Tevatron, corresponding to an integrated luminosity of 9.45 fb$^{-1}$. We consider events having no identified charged lepton, a transverse energy imbalance, and two or three jets, of which at least one is consistent with originating from the decay of a $b$ quark. We place 95% credibility level upper limits on the production cross section times standard model branching fraction for several mass hypotheses between 90 and $150 \mathrm{GeV}/c^2$. For a Higgs boson mass of $125 \mathrm{GeV}/c^2$, the observed (expected) limit is 6.7 (3.6) times the standard model prediction.Comment: Accepted by Phys. Rev. Let