Observation of muon intensity variations by season with the MINOS far detector

The temperature of the upper atmosphere affects the height of primary cosmic ray interactions and the production of high-energy cosmic ray muons which can be detected deep underground. The MINOS far detector at Soudan, MN, has collected over 67 X 10(6) cosmic ray induced muons. The underground muon rate measured over a period of five years exhibits a 4% peak-to-peak seasonal variation which is highly correlated with the temperature in the upper atmosphere. The coefficient, alpha(T), relating changes in the muon rate to changes in atmospheric temperature was found to be alpha(T) 0: 873 +/- 0: 009(stat) +/- 0.010(syst). Pions and kaons in the primary hadronic interactions of cosmic rays in the atmosphere contribute differently to alpha(T) due to the different masses and lifetimes. This allows the measured value of alpha(T) to be interpreted as a measurement of the K/pi ratio for E(p) greater than or similar to 7 TeV of 0.12(-0.05)(+0.07), consistent with the expectation from collider experiments.

Search for Muon-Neutrino to Electron-Neutrino Transitions in MINOS

This Letter reports on a search for nu(mu)->nu(e) transitions by the MINOS experiment based on a 3.14x10(20) protons-on-target exposure in the Fermilab NuMI beam. We observe 35 events in the Far Detector with a background of 27 +/- 5(stat)+/- 2(syst) events predicted by the measurements in the Near Detector. If interpreted in terms of nu(mu)->nu(e) oscillations, this 1.5 sigma excess of events is consistent with sin(2)(2 theta(13)) comparable to the CHOOZ limit when |Delta m(2)|=2.43x10(-3) eV(2) and sin(2)(2 theta(23))=1.0 are assumed.

First Direct Observation of Muon Antineutrino Disappearance

This Letter reports the first direct observation of muon antineutrino disappearance. The MINOS experiment has taken data with an accelerator beam optimized for (nu) over bar (mu) production, accumulating an exposure of 1.71 x 10(20) protons on target. In the Far Detector, 97 charged current (nu) over bar (mu) events are observed. The no-oscillation hypothesis predicts 156 events and is excluded at 6.3 sigma. The best fit to oscillation yields vertical bar Delta(m) over bar (2)vertical bar = [3.36(-0.40)(+0.46)(stat) +/- 0.06(sys)] x 10(-3) eV(2), sin(2)(2 (theta) over bar) = 0.86(-0.12)(+0.11)(stat) +/- 0.01(syst). The MINOS nu(mu) and (nu) over bar (mu) measurements are consistent at the 2.0% confidence level, assuming identical underlying oscillation parameters.

Search for Active Neutrino Disappearance Using Neutral-Current Interactions in the MINOS Long-Baseline Experiment

We report the first detailed comparisons of the rates and spectra of neutral-current neutrino interactions at two widely separated locations. A depletion in the rate at the far site would indicate mixing between nu(mu) and a sterile particle. No anomalous depletion in the reconstructed energy spectrum is observed. Assuming oscillations occur at a single mass-squared splitting, a fit to the neutral- and charged-current energy spectra limits the fraction of nu(mu) oscillating to a sterile neutrino to be below 0.68 at 90% confidence level. A less stringent limit due to a possible contribution to the measured neutral-current event rate at the far site from nu(e) appearance at the current experimental limit is also presented.

Testing Lorentz Invariance and CPT Conservation with NuMI Neutrinos in the MINOS Near Detector

A search for a sidereal modulation in the MINOS near detector neutrino data was performed. If present, this signature could be a consequence of Lorentz and CPT violation as predicted by the effective field theory called the standard-model extension. No evidence for a sidereal signal in the data set was found, implying that there is no significant change in neutrino propagation that depends on the direction of the neutrino beam in a sun-centered inertial frame. Upper limits on the magnitudes of the Lorentz and CPT violating terms in the standard-model extension lie between 10(-4) and 10(-2) of the maximum expected, assuming a suppression of these signatures by a factor of 10(-17).

Measurement of neutrino oscillations with the MINOS detectors in the NuMI beam

This Letter reports new results from the MINOS experiment based on a two-year exposure to muon neutrinos from the Fermilab NuMI beam. Our data are consistent with quantum-mechanical oscillations of neutrino flavor with mass splitting vertical bar Delta m(2)vertical bar = (2.43 +/- 0.13) x 10(-3) eV(2) (68% C.L.) and mixing angle sin(2)(2 theta) > 0.90 (90% C.L.). Our data disfavor two alternative explanations for the disappearance of neutrinos in flight: namely, neutrino decays into lighter particles and quantum decoherence of neutrinos, at the 3.7 and 5.7 standard-deviation levels, respectively.

Active to Sterile Neutrino Mixing Limits from Neutral-Current Interactions in MINOS

Results are reported from a search for active to sterile neutrino oscillations in the MINOS long-baseline experiment, based on the observation of neutral-current neutrino interactions, from an exposure to the NuMI neutrino beam of 7.07 x 10(20) protons on target. A total of 802 neutral-current event candidates is observed in the Far Detector, compared to an expected number of 754 +/- 28(stat) +/- 37(syst) for oscillations among three active flavors. The fraction f(s) of disappearing nu(mu) that may transition to nu(s) is found to be less than 22% at the 90% C.L.

Measurement of the Neutrino Mass Splitting and Flavor Mixing by MINOS

Measurements of neutrino oscillations using the disappearance of muon neutrinos from the Fermilab NuMI neutrino beam as observed by the two MINOS detectors are reported. New analysis methods have been applied to an enlarged data sample from an exposure of 7.25 x 10(20) protons on target. A fit to neutrino oscillations yields values of vertical bar Delta m(2)vertical bar = (2.32(-0.08)(+0.12) x 10(-3) eV(2) for the atmospheric mass splitting and sin(2)(2 theta) > 0.90 (90% C.L.) for the mixing angle. Pure neutrino decay and quantum decoherence hypotheses are excluded at 7 and 9 standard deviations, respectively.

TESTING STATISTICAL HYPOTHESIS ON RANDOM TREES AND APPLICATIONS TO THE PROTEIN CLASSIFICATION PROBLEM

Efficient automatic protein classification is of central importance in genomic annotation. As an independent way to check the reliability of the classification, we propose a statistical approach to test if two sets of protein domain sequences coming from two families of the Pfam database are significantly different. We model protein sequences as realizations of Variable Length Markov Chains (VLMC) and we use the context trees as a signature of each protein family. Our approach is based on a Kolmogorov-Smirnov-type goodness-of-fit test proposed by Balding et at. [Limit theorems for sequences of random trees (2008), DOI: 10.1007/s11749-008-0092-z]. The test statistic is a supremum over the space of trees of a function of the two samples; its computation grows, in principle, exponentially fast with the maximal number of nodes of the potential trees. We show how to transform this problem into a max-flow over a related graph which can be solved using a Ford-Fulkerson algorithm in polynomial time on that number. We apply the test to 10 randomly chosen protein domain families from the seed of Pfam-A database (high quality, manually curated families). The test shows that the distributions of context trees coming from different families are significantly different. We emphasize that this is a novel mathematical approach to validate the automatic clustering of sequences in any context. We also study the performance of the test via simulations on Galton-Watson related processes.

GEVREY SOLVABILITY AND GEVREY REGULARITY IN DIFFERENTIAL COMPLEXES ASSOCIATED TO LOCALLY INTEGRABLE STRUCTURES

In this work we study some properties of the differential complex associated to a locally integrable (involutive) structure acting on forms with Gevrey coefficients. Among other results we prove that, for such complexes, Gevrey solvability follows from smooth solvability under the sole assumption of a regularity condition. As a consequence we obtain the proof of the Gevrey solvability for a first order linear PDE with real-analytic coefficients satisfying the Nirenberg-Treves condition (P).

An oracle approach for interaction neighborhood estimation in random fields

We consider the problem of interaction neighborhood estimation from the partial observation of a finite number of realizations of a random field. We introduce a model selection rule to choose estimators of conditional probabilities among natural candidates. Our main result is an oracle inequality satisfied by the resulting estimator. We use then this selection rule in a two-step procedure to evaluate the interacting neighborhoods. The selection rule selects a small prior set of possible interacting points and a cutting step remove from this prior set the irrelevant points. We also prove that the Ising models satisfy the assumptions of the main theorems, without restrictions on the temperature, on the structure of the interacting graph or on the range of the interactions. It provides therefore a large class of applications for our results. We give a computationally efficient procedure in these models. We finally show the practical efficiency of our approach in a simulation study.

Chaos-Galerkin solution of stochastic Timoshenko bending problems

This paper presents an accurate and efficient solution for the random transverse and angular displacement fields of uncertain Timoshenko beams. Approximate, numerical solutions are obtained using the Galerkin method and chaos polynomials. The Chaos-Galerkin scheme is constructed by respecting the theoretical conditions for existence and uniqueness of the solution. Numerical results show fast convergence to the exact solution, at excellent accuracies. The developed Chaos-Galerkin scheme accurately approximates the complete cumulative distribution function of the displacement responses. The Chaos-Galerkin scheme developed herein is a theoretically sound and efficient method for the solution of stochastic problems in engineering. (C) 2011 Elsevier Ltd. All rights reserved.

Galerkin Solution of Stochastic Beam Bending on Winkler Foundations

In this paper, the Askey-Wiener scheme and the Galerkin method are used to obtain approximate solutions to stochastic beam bending on Winkler foundation. The study addresses Euler-Bernoulli beams with uncertainty in the bending stiffness modulus and in the stiffness of the foundation. Uncertainties are represented by parameterized stochastic processes. The random behavior of beam response is modeled using the Askey-Wiener scheme. One contribution of the paper is a sketch of proof of existence and uniqueness of the solution to problems involving fourth order operators applied to random fields. From the approximate Galerkin solution, expected value and variance of beam displacement responses are derived, and compared with corresponding estimates obtained via Monte Carlo simulation. Results show very fast convergence and excellent accuracies in comparison to Monte Carlo simulation. The Askey-Wiener Galerkin scheme presented herein is shown to be a theoretically solid and numerically efficient method for the solution of stochastic problems in engineering.

Post-failure behaviour of impulsively loaded circular plates

The dynamic plastic response of a simply supported circular plate is analysed. Emphasis is given to the plate behaviour after it has broken free from the supports due to a local material failure. The theoretical rigid plastic analysis predicts various features of the response such as the time to failure, residual kinetic energy and the critical velocity at failure. The residual kinetic energy of the plate could be significant enough to cause secondary impact damage. It is shown that the shape of the plate changes after breaking free from the supports, which is important for forensic investigations. The solution for various cases were proven to be exact in the context of the upper and lower bounds theorems of the theory of plasticity. (C) 2009 Elsevier Ltd. All rights reserved.

A Separation Principle for the Continuous-Time LQ-Problem With Markovian Jump Parameters

In this paper, we devise a separation principle for the finite horizon quadratic optimal control problem of continuous-time Markovian jump linear systems driven by a Wiener process and with partial observations. We assume that the output variable and the jump parameters are available to the controller. It is desired to design a dynamic Markovian jump controller such that the closed loop system minimizes the quadratic functional cost of the system over a finite horizon period of time. As in the case with no jumps, we show that an optimal controller can be obtained from two coupled Riccati differential equations, one associated to the optimal control problem when the state variable is available, and the other one associated to the optimal filtering problem. This is a separation principle for the finite horizon quadratic optimal control problem for continuous-time Markovian jump linear systems. For the case in which the matrices are all time-invariant we analyze the asymptotic behavior of the solution of the derived interconnected Riccati differential equations to the solution of the associated set of coupled algebraic Riccati equations as well as the mean square stabilizing property of this limiting solution. When there is only one mode of operation our results coincide with the traditional ones for the LQG control of continuous-time linear systems.