Search

The homogenization procedure developed here is conducted on a laminate with periodic space-time modulation on the fine scale: at leading order, this modulation creates convection in the low-wavelength regime if both parameters are modulated. However, if only one parameter is modulated, which is more realistic, this convective term disappears and one recovers a standard diffusion equation with effective homogeneous parameters; this does not describe the non-reciprocity and the propagation of the field observed from exact dispersion diagrams. This inconsistency is corrected here by considering second-order homogenization which results in a non-reciprocal propagation term that is proved to be non-zero for any laminate and verified via numerical simulation. The same methodology is also applied to the case when the density is modulated in the heat equation, leading therefore to a corrective advective term which cancels out non-reciprocity at the leading order but not at the second order.

Laminated media with material properties modulated in space and time in the form of travelling waves have long been known to exhibit non-reciprocity. However, when using the method of low frequency homogenisation, it was so far only possible to obtain non-reciprocal effective media when both material properties are modulated in time, in the form of a Willis-coupling (or bi-anisotropy in electromagnetism) model. If only one of the two properties is modulated in time, while the other is kept constant, it was thought impossible for the method of homogenisation to recover the expected non-reciprocity since this Willis-coupling coefficient then vanishes. Contrary to this belief, we show that effective media with a single time-modulated parameter are non-reciprocal, provided homogenization is pushed to the second order. This is illustrated by numerical experiments (dispersion diagrams and time-domain simulations) for a bilayered modulated medium.

We analyze a spatially extended version of a well-known model of forest-savanna dynamics, which presents as a system of nonlinear partial integro-differential equations, and study necessary conditions for pattern-forming bifurcations. Analytically, we show that homogeneous solutions dominate the dynamics of the standard forest-savanna model, regardless of the length scales of the various spatial processes considered. However, several different pattern-forming scenarios are possible upon including spatial resource limitation, such as competition for water, soil nutrients, or herbivory effects. Using numerical simulations and continuation, we study the nature of the resulting patterns as a function of system parameters and length scales, uncovering subcritical pattern-forming bifurcations and observing significant regions of multistability for realistic parameter regimes. Finally, we discuss our results in the context of extant savanna-forest modeling efforts and highlight ongoing challenges in building a unifying mathematical model for savannas across different rainfall levels., Comment: 27 pages, 7 figures

High-frequency homogenization is used to study dispersive media, containing inclusions placed periodically, for which the properties of the material depend on the frequency (Lorentz or Drude model with damping, for example). Effective properties are obtained near a given point of the dispersion diagram in frequency-wavenumber space. The asymptotic approximations of the dispersion diagrams, and the wavefields, so obtained are then cross-validated via detailed comparison with finite element method simulations in both one and two dimensions.

Spatial systems with heterogeneities are ubiquitous in nature, from precipitation, temperature and soil gradients controlling vegetation growth to morphogen gradients controlling gene expression in embryos. Such systems, generally described by nonlinear dynamical systems, often display complex parameter dependence and exhibit bifurcations. The dynamics of heterogeneous spatially extended systems passing through bifurcations are still relatively poorly understood, yet recent theoretical studies and experimental data highlight the resulting complex behaviors and their relevance to real-world applications. We explore the consequences of spatial heterogeneities passing through bifurcations via two examples strongly motivated by applications. These model systems illustrate that studying heterogeneity-induced behaviors in spatial systems is crucial for a better understanding of ecological transitions and functional organization in brain development., Comment: 28 pages

The space mission MICROSCOPE dedicated to the test of the Equivalence Principle (EP) operated from April 25, 2016 until the deactivation of the satellite on October 16, 2018. In this analysis we compare the free-fall accelerations ($a_{\rm A}$ and $a_{\rm B}$) of two test masses in terms of the E\"otv\"os parameter $\eta({\rm{A, B}}) = 2 \frac{a_{\rm A}- a_{\rm B}}{a_{\rm A}+ a_{\rm B}}$. No EP violation has been detected for two test masses, made from platinum and titanium alloys, in a sequence of 19 segments lasting from 13 to 198 hours down to the limit of the statistical error which is smaller than $10^{-14}$ for $ \eta({\rm{Ti, Pt}})$. Accumulating data from all segments leads to $\eta({\rm{Ti, Pt}}) =[-1.5\pm{}2.3{\rm (stat)}\pm{}1.5{\rm (syst)}] \times{}10^{-15}$ showing no EP violation at the level of $2.7\times{}10^{-15}$ if we combine stochastic and systematic errors quadratically. This represents an improvement of almost two orders of magnitude with respect to the previous best such test performed by the E\"ot-Wash group. The reliability of this limit has been verified by comparing the free falls of two test masses of the same composition (platinum) leading to a null E\"otv\"os parameter with a statistical uncertainty of $1.1\times{}10^{-15}$., Comment: Class. Quantum Grav. 39 204009

The MICROSCOPE mission was designed to test the Weak Equivalence Principle (WEP), stating the equality between the inertial and the gravitational masses, with a precision of $10^{-15}$ in terms of the E\"otv\"os ratio $\eta$. Its experimental test consisted of comparing the accelerations undergone by two collocated test masses of different compositions as they orbited the Earth, by measuring the electrostatic forces required to keep them in equilibrium. This was done with ultra-sensitive differential electrostatic accelerometers onboard a drag-free satellite. The mission lasted two and a half years, cumulating five-months-worth of science free-fall data, two thirds with a pair of test masses of different compositions -- Titanium and Platinum alloys -- and the last third with a reference pair of test masses of the same composition -- Platinum. We summarize the data analysis, with an emphasis on the characterization of the systematic uncertainties due to thermal instabilities and on the correction of short-lived events which could mimic a WEP violation signal. We found no violation of the WEP, with the E\"otv\"os parameter of the Titanium and Platinum pair constrained to $\eta({\rm Ti, Pt})~=~ [-1.5 \pm 2.3~{\rm (stat)} \pm 1.5~{\rm (syst)}]~\times 10^{-15}$ at $1\sigma$ in statistical errors., Comment: Phys. Rev. Lett. 129, 121102

A variety of nonlinear models of biological systems generate complex chaotic behaviors that contrast with biological homeostasis, the observation that many biological systems prove remarkably robust in the face of changing external or internal conditions. Motivated by the subtle dynamics of cell activity in a crustacean central pattern generator (CPG), this paper proposes a refinement of the notion of chaos that reconciles homeostasis and chaos in systems with multiple timescales. We show that systems displaying relaxation cycles while going through chaotic attractors generate chaotic dynamics that are regular at macroscopic timescales and are thus consistent with physiological function. We further show that this relative regularity may break down through global bifurcations of chaotic attractors such as crises, beyond which the system may also generate erratic activity at slow timescales. We analyze these phenomena in detail in the chaotic Rulkov map, a classical neuron model known to exhibit a variety of chaotic spike patterns. This leads us to propose that the passage of slow relaxation cycles through a chaotic attractor crisis is a robust, general mechanism for the transition between such dynamics. We validate this numerically in three other models: a simple model of the crustacean CPG neural network, a discrete cubic map, and a continuous flow.

Testing the Weak Equivalence Principle (WEP) to a precision of $10^{-15}$ requires a quantity of data that give enough confidence on the final result: ideally, the longer the measurement the better the rejection of thestatistical noise. The science sessions had a duration of 120 orbits maximum and were regularly repeated and spaced out to accommodate operational constraints but also in order to repeat the experiment in different conditions and to allow time to calibrate the instrument. Several science sessions were performed over the 2.5 year duration of the experiment. This paper aims to describe how the data have been produced on the basis of a mission scenario and a data flow process, driven by a tradeoff between the science objectives and the operational constraints. The mission was led by the Centre National d'Etudes Spatiales (CNES) which provided the satellite, the launch and the ground operations. The ground segment was distributed between CNES and Office National d'Etudes et de Recherches A\'erospatiales (ONERA). CNES provided the raw data through the Centre d'Expertise de Compensation de Tra\^{i}n\'{e}e (CECT: Drag-free expertise centre). The science was led by the Observatoire de la C\^ote d{'}Azur (OCA) and ONERA was in charge of the data process. The latter also provided the instrument and the Science Mission Centre of MICROSCOPE (CMSM)., Comment: To be released in special MICROSCOPE edition of CQG

Since the MICROSCOPE instrument aims to measure accelerations as low as a few 10$^{-15}$\,m\,s$^{-2}$ and cannot operate on ground, it was obvious to have a large time dedicated to its characterization in flight. After its release and first operation, the characterization experiments covered all the aspects of the instrument design in order to consolidate the scientific measurements and the subsequent conclusions drawn from them. Over the course of the mission we validated the servo-control and even updated the PID control laws for each inertial sensor. Thanks to several dedicated experiments and the analysis of the instrument sensitivities, we have been able to identify a number of instrument characteristics such as biases, gold wire and electrostatic stiffnesses, non linearities, couplings and free motion ranges of the test-masses, which may first impact the scientific objective and secondly the analysis of the instrument good operation., Comment: This paper is part of set of articles dedicated to the final results of MICROSCOPE space mission and submitted to CQG

The MICROSCOPE experiment was designed to test the weak equivalence principle in space, by comparing the low-frequency dynamics of cylindrical "free-falling" test masses controlled by electrostatic forces. We use data taken during technical sessions aimed at estimating the electrostatic stiffness of MICROSCOPE's sensors to constrain a short-range Yukawa deviation from Newtonian gravity. We take advantage of the fact that in the limit of small displacements, the gravitational interaction (both Newtonian and Yukawa-like) between nested cylinders is linear, and thus simply characterised by a stiffness. By measuring the total stiffness of the forces acting on a test mass as it moves, and comparing it with the theoretical electrostatic stiffness (expected to dominate), it is a priori possible to infer constraints on the Yukawa potential parameters. However, we find that measurement uncertainties are dominated by the gold wires used to control the electric charge of the test masses, though their related stiffness is indeed smaller than the expected electrostatic stiffness. Moreover, we find a non-zero unaccounted for stiffness that depends on the instrument's electric configuration, hinting at the presence of patch-field effects. Added to significant uncertainties on the electrostatic model, they only allow for poor constraints on the Yukawa potential. This is not surprising, as MICROSCOPE was not designed for this measurement, but this analysis is the first step to new experimental searches for non-Newtonian gravity in space., Comment: Accepted for publication in CQG's MICROSCOPE special issue

This article is dedicated to the use the MICROSCOPE mission's data to test chameleon theory of gravity. We take advantage of the technical sessions aimed to characterize the electrostatic stiffness of MICROSCOPE's instrument intrinsic to its capacitive measurement system. Any discrepancy between the expected and measured stiffness may result from unaccounted-for contributors, i.e. extra-forces. This work considers the case of chameleon gravity as a possible contributor. It was previously shown that in situations similar to these measurement sessions, a chameleon fifth force appears and acts as a stiffness for small displacements. The magnitude of this new component of the stiffness is computed over the chameleon's parameter space. It allows us to derive constraints by excluding any force inconsistent with the MICROSCOPE data. As expected --since MICROSCOPE was not designed for the purpose of such an analysis--, these new bounds are not competitive with state-of-the-art constraints, but they could be improved by a better estimation of all effects at play in these sessions. Hence our work illustrates this novel technique as a new way of constraining fifth forces., Comment: Submitted to PRD

In networks of nonlinear oscillators, symmetries place hard constraints on the system that can be exploited to predict universal dynamical features and steady-states, providing a rare generic organizing principle for far-from-equilibrium systems. However, the robustness of this class of theories to symmetry-disrupting imperfections is untested. Here, we develop a model experimental reaction-diffusion network of chemical oscillators to test applications of this theory in the context of self-organizing systems relevant to biology and soft robotics. The network is a ring of 4 identical microreactors containing the oscillatory Belousov-Zhabotinsky reaction coupled to nearest neighbors via diffusion. Assuming perfect symmetry, theory predicts 4 categories of stable spatiotemporal phase-locked periodic states and 4 categories of invariant manifolds that guide and structure transitions between phase-locked states. In our experiments, we observed the predicted symmetry-derived synchronous clustered transients that occur when the dynamical trajectories coincide with invariant manifolds. However, we observe only 3 of the 4 phase-locked states that are predicted for the idealized homogeneous system. Quantitative agreement between experiment and numerical simulations is found by accounting for the small amount of experimentally determined heterogeneity. This work demonstrates that a surprising degree of the network's dynamics are constrained by symmetry in spite of the breakdown of the assumption of homogeneity and raises the question of why heterogeneity destabilizes some symmetry predicted states, but not others., Comment: *Ian and Mike contributed equally to this work. This work is under review by Physical Review X

Dedicated accelerometers have been developed for the MICROSCOPE mission taking into account the specific range of acceleration to be measured on board the satellite. Considering one micro-g and even less as the full range of the instrument, leads to a customized concept and a high performance electronics for the sensing and servo-actuations of the accelerometer test-masses. In addition to a very accurate geometrical sensor core, a high performance electronics architecture provides the measurement of the weak electrostatic forces and torques applied to the test-masses. A set of capacitive sensors delivers the position and the attitude of the test-mass with respect to a very steady gold coated cage made in silica. The voltages applied on the electrodes surrounding each test-mass are finely controlled to generate the adequate electrical field and so the electrostatic pressures on the test-mass. This field maintains the test-mass motionless with respect to the instrument structure. Digital control laws are implemented in order to enable instrument operation flexibility and a weak position sensor noise. These electronics provide both the scientific data for MICROSCOPE's test of General Relativity and the data for the satellite drag-free and attitude control system (DFACS)., Comment: To be submitted to CQG's MICROSCOPE special issue

MICROSCOPE's space test of the weak equivalence principle (WEP) is based on the minute measurement of the difference of accelerations experienced by two test masses as they orbit the Earth. A detection of a violation of the WEP would appear at a well-known frequency $f_{\rm EP}$ depending on the satellite's orbital and spinning frequencies. Consequently, the experiment was optimised to miminise systematic errors at $f_{\rm EP}$. Glitches are short-lived events visible in the test masses' measured acceleration, most likely originating in cracks of the satellite's coating. In this paper, we characterise their shape and time distribution. Although intrinsically random, their time of arrival distribution is modulated by the orbital and spinning periods. They have an impact on the WEP test that must be quantified. However, the data available prevents us from unequivocally tackling this task. We show that glitches affect the test of the WEP, up to an a priori unknown level. Discarding the perturbed data is thus the best way to reduce their effect., Comment: References updated

After performing highly sensitive acceleration measurements during two years of drag-free flight around the Earth, MICROSCOPE provided the best constraint on the Weak Equivalence Principle (WEP) to date. Beside being a technological challenge, this experiment required a specialised data analysis pipeline to look for a potential small signal buried in the noise, possibly plagued by instrumental defects, missing data and glitches. This paper describes the frequency-domain iterative least-square technique that we developed for MICROSCOPE. In particular, using numerical simulations, we prove that our estimator is unbiased and provides correct error bars. This paper therefore justifies the robustness of the WEP measurements given by MICROSCOPE., Comment: References updated

This paper focus on the description of the design and performance of the MICROSCOPE satellite and its Drag-Free and Attitude Control System (DFACS). The satellite is derived from CNES' Myriade platform family, albeit with significant upgrades dictated by the unprecedented MICROSCOPE's mission requirements. The 300kg drag-free microsatellite has completed its 2-year flight with higher-than-expected performances. Its passive thermal concept allowed for variations smaller than 1 $\mu$K at the measurement frequency $f_{\rm{EP}}$. The propulsion system provided a 6 axis continuous and very low noise thrust from zero to some hundreds of micronewtons. Finally, the performance of its DFACS (aimed at compensating the disturbing forces and torques applied to the satellite) is the finest ever achieved in low Earth orbit, with residual accelerations along the three axes are lower than $10^{-12} {\rm m/s}^2$ at $f_{\rm{EP}}$ over 8 days., Comment: References updated

The MICROSCOPE mission aimed to test the Weak Equivalence Principle (WEP) to a precision of $10^{-15}$. The WEP states that two bodies fall at the same rate on a gravitational field independently of their mass or composition. In MICROSCOPE, two masses of different compositions (titanium and platinum alloys) are placed on a quasi-circular trajectory around the Earth. They are the test-masses of a double accelerometer. The measurement of their accelerations is used to extract a potential WEP violation that would occur at a frequency defined by the motion and attitude of the satellite around the Earth. This paper details the major drivers of the mission leading to the specification of the major subsystems (satellite, ground segment, instrument, orbit...). Building upon the measurement equation, we derive the objective of the test in statistical and systematic error allocation and provide the mission's expected error budget., Comment: References updated

Considering two random variables with different laws to which we only have access through finite size iid samples, we address how to reweight the first sample so that its empirical distribution converges towards the true law of the second sample as the size of both samples goes to infinity. We study an optimal reweighting that minimizes the Wasserstein distance between the empirical measures of the two samples, and leads to an expression of the weights in terms of Nearest Neighbors. The consistency and some asymptotic convergence rates in terms of expected Wasserstein distance are derived, and do not need the assumption of absolute continuity of one random variable with respect to the other. These results have some application in Uncertainty Quantification for decoupled estimation and in the bound of the generalization error for the Nearest Neighbor Regression under covariate shift.

Experimental characterization of neuronal dynamics involves recording both of spontaneous activity patterns and of responses to transient and sustained inputs. While much theoretical attention has been devoted to the spontaneous activity of neurons, less is known about the dynamic mechanisms shaping their responses to transient inputs, although these bear significant physiological relevance. Here, we study responses to transient inputs in a widely used class of neuron models (nonlinear adaptive hybrid models) well-known to reproduce a number of biologically realistic behaviors. We focus on responses to transient inputs that have been previously associated with Type III neurons, arguably the least studied category in Hodgkin's classification, which are those neurons that never exhibit continuous firing in response to sustained excitatory currents. The two phenomena that we study are post-inhibitory facilitation, in which an otherwise subthreshold excitatory input can induce a spike if it is applied with proper timing after an inhibitory pulse, and slope detection, in which a neuron spikes to a transient input only when the input's rate of change is in a specific, bounded range. We analyze the origin of these phenomena in nonlinear hybrid models and provide a geometric characterization of dynamical structures associated with PIF in the system and an analytical study of slope detection for tent inputs. While the necessary and sufficient conditions for these behaviors are easily satisfied in neurons with Type III excitability, our proofs are quite general and valid for neurons that do not exhibit Type III excitability as well. This study therefore provides a framework for the mathematical analysis of these responses to transient inputs associated with Type III neurons in other systems and for advancing our understanding of these systems' computational properties.

Many studies have found evidence that the brain operates at a critical point, a processus known as self-organized criticality. A recent paper found remarkable scalings suggestive of criticality in systems as different as neural cultures, anesthetized or awake brains. We point out here that the diversity of these states would question any claimed role of criticality in information processing. Furthermore, we show that two non-critical systems pass all the tests for criticality, a control that was not provided in the original article. We conclude that such false positives demonstrate that the presence of criticality in the brain is still not proven and that we need better methods that scaling analyses.

We use data from the T-SAGE instrument on board the MICROSCOPE space mission to search for Lorentz violation in matter-gravity couplings as described by the Lorentz violating Standard-Model Extension (SME) coefficients $(\bar{a}_\text{eff})_\mu^w$, where ($\mu = T,X,Y,Z$) and ($w = e,p,n$) for the electron, proton and neutron. One of the phenomenological consequences of a non-zero value of those coefficients is that test bodies of different composition fall differently in an external gravitational field. This is similar to "standard" tests of the universality of free fall, but with a specific signature that depends on the orbital velocity and rotation of the Earth. We analyze data from five measurement sessions of MICROSCOPE spread over a year finding no evidence for such a signature, but setting constraints on linear combinations of the SME coefficients that improve on best previous results by one to two orders of magnitude. Additionally, our independent linear combinations are different from previous ones, which increases the diversity of available constraints, paving the way towards a full decorrelation of the individual coefficients.

Deterministic models of vegetation often summarize, at a macroscopic scale, a multitude of intrinsically random events occurring at a microscopic scale. We bridge the gap between these scales by demonstrating convergence to a mean-field limit for a general class of stochastic models representing each individual ecological event in the limit of large system size. The proof relies on classical stochastic coupling techniques that we generalize to cover spatially extended interactions. The mean-field limit is a spatially extended non-Markovian process characterized by nonlocal integro-differential equations describing the evolution of the probability for a patch of land to be in a given state (the generalized Kolmogorov equations of the process, GKEs). We thus provide an accessible general framework for spatially extending many classical finite-state models from ecology and population dynamics. We demonstrate the practical effectiveness of our approach through a detailed comparison of our limiting spatial model and the finite-size version of a specific savanna-forest model, the so-called Staver-Levin model. There is remarkable dynamic consistency between the GKEs and the finite-size system, in spite of almost sure forest extinction in the finite-size system. To resolve this apparent paradox, we show that the extinction rate drops sharply when nontrivial equilibria emerge in the GKEs, and that the finite-size system's quasi-stationary distribution (stationary distribution conditional on non-extinction) closely matches the bifurcation diagram of the GKEs. Furthermore, the limit process can support periodic oscillations of the probability distribution, thus providing an elementary example of a jump process that does not converge to a stationary distribution. In spatially extended settings, environmental heterogeneity can lead to waves of invasion and front-pinning phenomena., Comment: 38 pages

The Weak Equivalence Principle (WEP), stating that two bodies of different compositions and/or mass fall at the same rate in a gravitational field (universality of free fall), is at the very foundation of General Relativity. The MICROSCOPE mission aims to test its validity to a precision of $10^{-15}$, two orders of magnitude better than current on-ground tests, by using two masses of different compositions (titanium and platinum alloys) on a quasi-circular trajectory around the Earth. This is realised by measuring the accelerations inferred from the forces required to maintain the two masses exactly in the same orbit. Any significant difference between the measured accelerations, occurring at a defined frequency, would correspond to the detection of a violation of the WEP, or to the discovery of a tiny new type of force added to gravity. MICROSCOPE's first results show no hint for such a difference, expressed in terms of E\"otv\"os parameter $\delta(Ti,Pt)=[-1\pm{}9{\rm (stat)}\pm{}9{\rm (syst)}] \times{}10^{-15}$ (both 1$\sigma$ uncertainties) for a titanium and platinum pair of materials. This result was obtained on a session with 120 orbital revolutions representing 7\% of the current available data acquired during the whole mission. The quadratic combination of 1$\sigma$ uncertainties leads to a current limit on $\delta$ of about $1.3\times{}10^{-14}$., Comment: To appear in CQG

We study the statistical physics of a surprising phenomenon arising in large networks of excitable elements in response to noise: while at low noise, solutions remain in the vicinity of the resting state and large-noise solutions show asynchronous activity, the network displays orderly, perfectly synchronized periodic responses at intermediate level of noise. We show that this phenomenon is fundamentally stochastic and collective in nature. Indeed, for noise and coupling within specific ranges, an asymmetry in the transition rates between a resting and an excited regime progressively builds up, leading to an increase in the fraction of excited neurons eventually triggering a chain reaction associated with a macroscopic synchronized excursion and a collective return to rest where this process starts afresh, thus yielding the observed periodic synchronized oscillations. We further uncover a novel anti-resonance phenomenon: noise-induced synchronized oscillations disappear when the system is driven by periodic stimulation with frequency within a specific range. In that anti-resonance regime, the system is optimal for measures of information capacity. This observation provides a new hypothesis accounting for the efficiency of Deep Brain Stimulation therapies in Parkinson's disease, a neurodegenerative disease characterized by an increased synchronization of brain motor circuits. We further discuss the universality of these phenomena in the class of stochastic networks of excitable elements with confining coupling, and illustrate this universality by analyzing various classical models of neuronal networks. Altogether, these results uncover some universal mechanisms supporting a regularizing impact of noise in excitable systems, reveal a novel anti-resonance phenomenon in these systems, and propose a new hypothesis for the efficiency of high-frequency stimulation in Parkinson's disease.

We introduce one- and two-dimensional `exponential shapelets': orthonormal basis functions that efficiently model isolated features in data. They are built from eigenfunctions of the quantum mechanical hydrogen atom, and inherit mathematics with elegant properties under Fourier transform, and hence (de)convolution. For a wide variety of data, exponential shapelets compress information better than Gauss-Hermite/Gauss-Laguerre (`shapelet') decomposition, and generalise previous attempts that were limited to 1D or circularly symmetric basis functions. We discuss example applications in astronomy, fundamental physics and space geodesy., Comment: Accepted for publication in MNRAS. 13+4 pages