319 research outputs found
Non-Gaussianity from extragalactic point-sources
The population of compact extragalactic sources contribute to the
non-Gaussianity at Cosmic Microwave Background frequencies. We study their
non-Gaussianity using publicly available full-sky simulations. We introduce a
parametrisation to visualise efficiently the bispectrum and we describe the
scale and frequency dependences of the bispectrum of radio and IR
point-sources. We show that the bispectrum is well fitted by an analytical
prescription. We find that the clustering of IR sources enhances their
non-Gaussianity by several orders of magnitude, and that their bispectrum peaks
in the squeezed triangles. Examining the impact of these sources on primordial
non-Gaussianity estimation, we find that radio sources yield an important
positive bias to local fNL at low frequencies but this bias is efficiently
reduced by masking detectable sources. IR sources produce a negative bias at
high frequencies, which is not dimmed by the masking, as their clustering is
dominated by faint sources.Comment: 4pages, 2 figures, 2 tables. Contribution to the proceedings of the
International Conference on Gravitation and Cosmology, Goa, December 201
Time series irreversibility: a visibility graph approach
We propose a method to measure real-valued time series irreversibility which
combines two differ- ent tools: the horizontal visibility algorithm and the
Kullback-Leibler divergence. This method maps a time series to a directed
network according to a geometric criterion. The degree of irreversibility of
the series is then estimated by the Kullback-Leibler divergence (i.e. the
distinguishability) between the in and out degree distributions of the
associated graph. The method is computationally effi- cient, does not require
any ad hoc symbolization process, and naturally takes into account multiple
scales. We find that the method correctly distinguishes between reversible and
irreversible station- ary time series, including analytical and numerical
studies of its performance for: (i) reversible stochastic processes
(uncorrelated and Gaussian linearly correlated), (ii) irreversible stochastic
pro- cesses (a discrete flashing ratchet in an asymmetric potential), (iii)
reversible (conservative) and irreversible (dissipative) chaotic maps, and (iv)
dissipative chaotic maps in the presence of noise. Two alternative graph
functionals, the degree and the degree-degree distributions, can be used as the
Kullback-Leibler divergence argument. The former is simpler and more intuitive
and can be used as a benchmark, but in the case of an irreversible process with
null net current, the degree-degree distribution has to be considered to
identifiy the irreversible nature of the series.Comment: submitted for publicatio
A flexible method for optimising sharing of healthcare resources and demand in the context of the COVID-19 pandemic.
As the number of cases of COVID-19 continues to grow, local health services are at risk of being overwhelmed with patients requiring intensive care. We develop and implement an algorithm to provide optimal re-routing strategies to either transfer patients requiring Intensive Care Units (ICU) or ventilators, constrained by feasibility of transfer. We validate our approach with realistic data from the United Kingdom and Spain. In the UK, we consider the National Health Service at the level of trusts and define a 4-regular geometric graph which indicates the four nearest neighbours of any given trust. In Spain we coarse-grain the healthcare system at the level of autonomous communities, and extract similar contact networks. Through random search optimisation we identify the best load sharing strategy, where the cost function to minimise is based on the total number of ICU units above capacity. Our framework is general and flexible allowing for additional criteria, alternative cost functions, and can be extended to other resources beyond ICU units or ventilators. Assuming a uniform ICU demand, we show that it is possible to enable access to ICU for up to 1000 additional cases in the UK in a single step of the algorithm. Under a more realistic and heterogeneous demand, our method is able to balance about 600 beds per step in the Spanish system only using local sharing, and over 1300 using countrywide sharing, potentially saving a large percentage of these lives that would otherwise not have access to ICU
The shape of memory in temporal networks
Temporal networks are widely used models for describing the architecture of
complex systems. Network memory -- that is the dependence of a temporal
network's structure on its past -- has been shown to play a prominent role in
diffusion, epidemics and other processes occurring over the network, and even
to alter its community structure. Recent works have proposed to estimate the
length of memory in a temporal network by using high-order Markov models. Here
we show that network memory is inherently multidimensional and cannot be
meaningfully reduced to a single scalar quantity. Accordingly, we introduce a
mathematical framework for defining and efficiently estimating the microscopic
shape of memory, which fully characterises how the activity of each link
intertwines with the activities of all other links. We validate our methodology
on a wide range of synthetic models of temporal networks with tuneable memory,
and subsequently study the heterogeneous shapes of memory emerging in various
real-world networks.Comment: 35 pages (5 main, 30 supplementary), 14 figures (3 main, 11
supplementary), 3 tables (all supplementary), uses tikz-network.sty and
tikz_network.p
Emergence of collective intonation in the musical performance of crowds
To be published in EPLTo be published in EP
A combinatorial framework to quantify peak/pit asymmetries in complex dynamics
LL’s acknowledges funding from an EPSRC Early Career Fellowship EP/P01660X/1
The Visibility Graph: a new method for estimating the Hurst exponent of fractional Brownian motion
Fractional Brownian motion (fBm) has been used as a theoretical framework to
study real time series appearing in diverse scientific fields. Because its
intrinsic non-stationarity and long range dependence, its characterization via
the Hurst parameter H requires sophisticated techniques that often yield
ambiguous results. In this work we show that fBm series map into a scale free
visibility graph whose degree distribution is a function of H. Concretely, it
is shown that the exponent of the power law degree distribution depends
linearly on H. This also applies to fractional Gaussian noises (fGn) and
generic f^(-b) noises. Taking advantage of these facts, we propose a brand new
methodology to quantify long range dependence in these series. Its reliability
is confirmed with extensive numerical simulations and analytical developments.
Finally, we illustrate this method quantifying the persistent behavior of human
gait dynamics.Comment: 5 pages, submitted for publicatio
Multiplex Decomposition of Non-Markovian Dynamics and the Hidden Layer Reconstruction Problem
Elements composing complex systems usually interact in several different ways
and as such the interaction architecture is well modelled by a multiplex
network. However often this architecture is hidden, as one usually only has
experimental access to an aggregated projection. A fundamental challenge is
thus to determine whether the hidden underlying architecture of complex systems
is better modelled as a single interaction layer or results from the
aggregation and interplay of multiple layers. Here we show that using local
information provided by a random walker navigating the aggregated network one
can decide in a robust way if the underlying structure is a multiplex or not
and, in the former case, to determine the most probable number of hidden
layers. As a byproduct, we show that the mathematical formalism also provides a
principled solution for the optimal decomposition and projection of complex,
non-Markovian dynamics into a Markov switching combination of diffusive modes.
We validate the proposed methodology with numerical simulations of both (i)
random walks navigating hidden multiplex networks (thereby reconstructing the
true hidden architecture) and (ii) Markovian and non-Markovian continuous
stochastic processes (thereby reconstructing an effective multiplex
decomposition where each layer accounts for a different diffusive mode). We
also state and prove two existence theorems guaranteeing that an exact
reconstruction of the dynamics in terms of these hidden jump-Markov models is
always possible for arbitrary finite-order Markovian and fully non-Markovian
processes. Finally, we showcase the applicability of the method to experimental
recordings from (i) the mobility dynamics of human players in an online
multiplayer game and (ii) the dynamics of RNA polymerases at the
single-molecule level.Comment: 40 pages, 24 figure
Development of a tomato spotted wilt virus (TSWV) risk evaluation methology for a processing tomato region
A risk map for the Tomato spotted wilt virus (TSWV) was elaborated for the main Portuguese processing tomato
producing region, the “Ribatejo e Península de Setúbal” region, where periodically this virus causes severe losses.
Forty nine tomato fields were monitored. Risk factors for TSWV infection were identified and quantified according to
their relative importance in TSWV incidence. The risk factors considered for each field were: (1) presence of TSWV
in tomato plants; (2) presence of TSWV in weeds which are hosts of TSWV vectors; (3) presence of TSWV vector
thrips; (4) presence of TSWV host crops previously (in the two years before), namely, tomato, potato and sweet pepper;
and (5) presence of greenhouses, urban areas or TSWV host crops next to the field (up to about 100m from its borders).
A risk estimator was calculated for each field. Among the thrips (Thysanoptera) identified, belonging to 11 genera, four
vector thrips species were detected: Frankliniella occidentalis (Pergande) and Thrips tabaci Lindman, the two most
abundant ones, and F. intonsa (Trybom) and F. schultzei (Trybom). Blue sticky traps placed up to about 75 cm above
the crop canopy caught F. occidentalis and T. tabaci more efficiently than the beating technique. The weeds Datura
stramonium L., Arctotheca calendula (L.), and Conyza bonariensis (L.) were identified as TSWV winter repositories.
This study proposes a methodology to be used by field technicians for the annual evaluation of TSWV risk at a regional
level, for an improved planning of processing tomato crop in the following season
Critical behavior of a Ginzburg-Landau model with additive quenched noise
We address a mean-field zero-temperature Ginzburg-Landau, or \phi^4, model
subjected to quenched additive noise, which has been used recently as a
framework for analyzing collective effects induced by diversity. We first make
use of a self-consistent theory to calculate the phase diagram of the system,
predicting the onset of an order-disorder critical transition at a critical
value {\sigma}c of the quenched noise intensity \sigma, with critical exponents
that follow Landau theory of thermal phase transitions. We subsequently perform
a numerical integration of the system's dynamical variables in order to compare
the analytical results (valid in the thermodynamic limit and associated to the
ground state of the global Lyapunov potential) with the stationary state of the
(finite size) system. In the region of the parameter space where metastability
is absent (and therefore the stationary state coincide with the ground state of
the Lyapunov potential), a finite-size scaling analysis of the order parameter
fluctuations suggests that the magnetic susceptibility diverges quadratically
in the vicinity of the transition, what constitutes a violation of the
fluctuation-dissipation relation. We derive an effective Hamiltonian and
accordingly argue that its functional form does not allow to straightforwardly
relate the order parameter fluctuations to the linear response of the system,
at odds with equilibrium theory. In the region of the parameter space where the
system is susceptible to have a large number of metastable states (and
therefore the stationary state does not necessarily correspond to the ground
state of the global Lyapunov potential), we numerically find a phase diagram
that strongly depends on the initial conditions of the dynamical variables.Comment: 8 figure
- …