Regression analysis with missing data and unknown colored noise: application to the MICROSCOPE space mission

The analysis of physical measurements often copes with highly correlated noises and interruptions caused by outliers, saturation events or transmission losses. We assess the impact of missing data on the performance of linear regression analysis involving the fit of modeled or measured time series. We show that data gaps can significantly alter the precision of the regression parameter estimation in the presence of colored noise, due to the frequency leakage of the noise power. We present a regression method which cancels this effect and estimates the parameters of interest with a precision comparable to the complete data case, even if the noise power spectral density (PSD) is not known a priori. The method is based on an autoregressive (AR) fit of the noise, which allows us to build an approximate generalized least squares estimator approaching the minimal variance bound. The method, which can be applied to any similar data processing, is tested on simulated measurements of the MICROSCOPE space mission, whose goal is to test the Weak Equivalence Principle (WEP) with a precision of $10^{-15}$. In this particular context the signal of interest is the WEP violation signal expected to be found around a well defined frequency. We test our method with different gap patterns and noise of known PSD and find that the results agree with the mission requirements, decreasing the uncertainty by a factor 60 with respect to ordinary least squares methods. We show that it also provides a test of significance to assess the uncertainty of the measurement.Comment: 12 pages, 4 figures, to be published in Phys. Rev.

Finite-size and correlation-induced effects in Mean-field Dynamics

The brain's activity is characterized by the interaction of a very large number of neurons that are strongly affected by noise. However, signals often arise at macroscopic scales integrating the effect of many neurons into a reliable pattern of activity. In order to study such large neuronal assemblies, one is often led to derive mean-field limits summarizing the effect of the interaction of a large number of neurons into an effective signal. Classical mean-field approaches consider the evolution of a deterministic variable, the mean activity, thus neglecting the stochastic nature of neural behavior. In this article, we build upon two recent approaches that include correlations and higher order moments in mean-field equations, and study how these stochastic effects influence the solutions of the mean-field equations, both in the limit of an infinite number of neurons and for large yet finite networks. We introduce a new model, the infinite model, which arises from both equations by a rescaling of the variables and, which is invertible for finite-size networks, and hence, provides equivalent equations to those previously derived models. The study of this model allows us to understand qualitative behavior of such large-scale networks. We show that, though the solutions of the deterministic mean-field equation constitute uncorrelated solutions of the new mean-field equations, the stability properties of limit cycles are modified by the presence of correlations, and additional non-trivial behaviors including periodic orbits appear when there were none in the mean field. The origin of all these behaviors is then explored in finite-size networks where interesting mesoscopic scale effects appear. This study leads us to show that the infinite-size system appears as a singular limit of the network equations, and for any finite network, the system will differ from the infinite system

Matter-gravity couplings and Lorentz violation

The gravitational couplings of matter are studied in the presence of Lorentz and CPT violation. At leading order in the coefficients for Lorentz violation, the relativistic quantum hamiltonian is derived from the gravitationally coupled minimal Standard-Model Extension. For spin-independent effects, the nonrelativistic quantum hamiltonian and the classical dynamics for test and source bodies are obtained. A systematic perturbative method is developed to treat small metric and coefficient fluctuations about a Lorentz-violating and Minkowski background. The post-newtonian metric and the trajectory of a test body freely falling under gravity in the presence of Lorentz violation are established. An illustrative example is presented for a bumblebee model. The general methodology is used to identify observable signals of Lorentz and CPT violation in a variety of gravitational experiments and observations, including gravimeter measurements, laboratory and satellite tests of the weak equivalence principle, antimatter studies, solar-system observations, and investigations of the gravitational properties of light. Numerous sensitivities to coefficients for Lorentz violation can be achieved in existing or near-future experiments at the level of parts in 10^3 down to parts in 10^{15}. Certain coefficients are uniquely detectable in gravitational searches and remain unmeasured to date.Comment: 59 pages two-column REVTe

A Markovian event-based framework for stochastic spiking neural networks

In spiking neural networks, the information is conveyed by the spike times, that depend on the intrinsic dynamics of each neuron, the input they receive and on the connections between neurons. In this article we study the Markovian nature of the sequence of spike times in stochastic neural networks, and in particular the ability to deduce from a spike train the next spike time, and therefore produce a description of the network activity only based on the spike times regardless of the membrane potential process. To study this question in a rigorous manner, we introduce and study an event-based description of networks of noisy integrate-and-fire neurons, i.e. that is based on the computation of the spike times. We show that the firing times of the neurons in the networks constitute a Markov chain, whose transition probability is related to the probability distribution of the interspike interval of the neurons in the network. In the cases where the Markovian model can be developed, the transition probability is explicitly derived in such classical cases of neural networks as the linear integrate-and-fire neuron models with excitatory and inhibitory interactions, for different types of synapses, possibly featuring noisy synaptic integration, transmission delays and absolute and relative refractory period. This covers most of the cases that have been investigated in the event-based description of spiking deterministic neural networks

Testing the equivalence principle: why and how?

Part of the theoretical motivation for improving the present level of testing of the equivalence principle is reviewed. The general rationale for optimizing the choice of pairs of materials to be tested is presented. One introduces a simplified rationale based on a trichotomy of competing classes of theoretical models.Comment: 11 pages, Latex, uses ioplppt.sty, submitted to Class. Quantum Gra

Runaway dilaton and equivalence principle violations

In a recently proposed scenario, where the dilaton decouples while cosmologically attracted towards infinite bare string coupling, its residual interactions can be related to the amplitude of density fluctuations generated during inflation, and are large enough to be detectable through a modest improvement on present tests of free-fall universality. Provided it has significant couplings to either dark matter or dark energy, a runaway dilaton can also induce time-variations of the natural "constants" within the reach of near-future experiments.Comment: 4 pages, minor change

On the equivalence principle and gravitational and inertial mass relation of classical charged particles

We show that the locally constant force necessary to get a stable hyperbolic motion regime for classical charged point particles, actually, is a combination of an applied external force and of the electromagnetic radiation reaction force. It implies, as the strong Equivalence Principle is valid, that the passive gravitational mass of a charged point particle should be slight greater than its inertial mass. An interesting new feature that emerges from the unexpected behavior of the gravitational and inertial mass relation, for classical charged particles, at very strong gravitational field, is the existence of a critical, particle dependent, gravitational field value that signs the validity domain of the strong Equivalence Principle. For electron and proton, these critical field values are $g_{c}\simeq 4.8\times 10^{31}m/s^{2}$ and $g_{c}\simeq 8.8\times 10^{34}m/s^{2}$, respectively

Quantum fluctuations for drag free geodesic motion

The drag free technique is used to force a proof mass to follow a geodesic motion. The mass is protected from perturbations by a cage, and the motion of the latter is actively controlled to follow the motion of the proof mass. We present a theoretical analysis of the effects of quantum fluctuations for this technique. We show that a perfect drag free operation is in principle possible at the quantum level, in spite of the back action exerted on the mass by the position sensor.Comment: 4 pages, 1 figure, RevTeX, minor change

Time Variation of the Fine Structure Constant Driven by Quintessence

There are indications from the study of quasar absorption spectra that the fine structure constant $\alpha$ may have been measurably smaller for redshifts $z>2.$ Analyses of other data ($^{149}$Sm fission rate for the Oklo natural reactor, variation of $^{187}$Re $\beta$-decay rate in meteorite studies, atomic clock measurements) which probe variations of $\alpha$ in the more recent past imply much smaller deviations from its present value. In this work we tie the variation of $\alpha$ to the evolution of the quintessence field proposed by Albrecht and Skordis, and show that agreement with all these data, as well as consistency with WMAP observations, can be achieved for a range of parameters. Some definite predictions follow for upcoming space missions searching for violations of the equivalence principle.Comment: Final version, to be published in Phys Rev

Acceleration disturbances and requirements for ASTROD I

ASTRODynamical Space Test of Relativity using Optical Devices I (ASTROD I) mainly aims at testing relativistic gravity and measuring the solar-system parameters with high precision, by carrying out laser ranging between a spacecraft in a solar orbit and ground stations. In order to achieve these goals, the magnitude of the total acceleration disturbance of the proof mass has to be less than 10−13 m s−2 Hz−1/2 at 0.1 m Hz. In this paper, we give a preliminary overview of the sources and magnitude of acceleration disturbances that could arise in the ASTROD I proof mass. Based on the estimates of the acceleration disturbances and by assuming a simple controlloop model, we infer requirements for ASTROD I. Our estimates show that most of the requirements for ASTROD I can be relaxed in comparison with Laser Interferometer Space Antenna (LISA).Comment: 19 pages, two figures, accepted for publication by Class. Quantum Grav. (at press