Source author record

Pedro G. Lind

Pedro G. Lind appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
21topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2026arXiv

Explainability of Complex AI Models with Correlation Impact Ratio

Complex AI systems make better predictions but often lack transparency, limiting trustworthiness, interpretability, and safe deployment. Common post hoc AI explainers, such as LIME, SHAP, HSIC, and SAGE, are model agnostic but are too restricted in one significant regard: they tend to misrank correlated features and require costly perturbations, which do not scale to high dimensional data. We introduce ExCIR (Explainability through Correlation Impact Ratio), a theoretically grounded, simple, and reliable metric for explaining the contribution of input features to model outputs, which remains stable and consistent under noise and sampling variations. We demonstrate that ExCIR captures dependencies arising from correlated features through a lightweight single pass formulation. Experimental evaluations on diverse datasets, including EEG, synthetic vehicular data, Digits, and Cats-Dogs, validate the effectiveness and stability of ExCIR across domains, achieving more interpretable feature explanations than existing methods while remaining computationally efficient. To this end, we further extend ExCIR with an information theoretic foundation that unifies the correlation ratio with Canonical Correlation Analysis under mutual information bounds, enabling multi output and class conditioned explainability at scale.

preprint2018arXiv

Bridging between Load-Flow and Kuramoto-like Power Grid Models: A Flexible Approach to Integrating Electrical Storage Units

In future power systems, electrical storage will be the key technology for balancing feed-in fluctuations. With increasing share of renewables and reduction of system inertia, the focus of research expands towards short-term grid dynamics and collective phenomena. Against this backdrop, Kuramoto-like power grids have been established as a sound mathematical modeling framework bridging between the simplified models from nonlinear dynamics and the more detailed models used in electrical engineering. However, they have a blind spot concerning grid components, which cannot be modeled by oscillator equations, and hence do not allow to investigate storage-related issues from scratch. We remove this shortcoming by bringing together Kuramoto-like and algebraic load-flow equations. This is a substantial extension of the current Kuramoto framework with arbitrary grid components. Based on this concept, we provide a solid starting point for the integration of flexible storage units enabling to address current problems like smart storage control, optimal siting and rough cost estimations. For demonstration purpose, we here consider a wind power application with realistic feed-in conditions. We show how to implement basic control strategies from electrical engineering, give insights into their potential with respect to frequency quality improvement and point out their limitations by maximum capacity and finite-time response.

preprint2016arXiv

The Langevin Approach: An R Package for Modeling Markov Processes

We describe an R package developed by the research group Turbulence, Wind energy and Stochastics (TWiSt) at the Carl von Ossietzky University of Oldenburg, which extracts the (stochastic) evolution equation underlying a set of data or measurements. The method can be directly applied to data sets with one or two stochastic variables. Examples for the one-dimensional and two-dimensional cases are provided. This framework is valid under a small set of conditions which are explicitly presented and which imply simple preliminary test procedures to the data. For Markovian processes involving Gaussian white noise, a stochastic differential equation is derived straightforwardly from the time series and captures the full dynamical properties of the underlying process. Still, even in the case such conditions are not fulfilled, there are alternative versions of this method which we discuss briefly and provide the user with the necessary bibliography.

preprint2015arXiv

From empirical data to continuous Markov processes: a systematic approach

We present an approach for testing for the existence of continuous generators of discrete stochastic transition matrices. Typically, the known approaches to ascertain the existence of continuous Markov processes are based in the assumption that only time-homogeneous generators exist. Here, a systematic extension to time-inhomogeneity is presented, based in new mathematical propositions incorporating necessary and sufficient conditions, which are then implemented computationally and applied to numerical data. A discussion concerning the bridging between rigorous mathematical results on the existence of generators to its computational implementation. Our detection algorithm shows to be effective in more than $80\%$ of tested matrices, typically $90\%$ to $95\%$, and for those an estimate of the (non-homogeneous) generator matrix follows. We also solve the embedding problem analytically for the particular case of three-dimensional circulant matrices. Finally, a discussion of possible applications of our framework to problems in different fields is briefly addressed.

preprint2015arXiv

Parameter-free resolution of the superposition of stochastic signals

This paper presents a direct method to obtain the deterministic and stochastic contribution of the sum of two independent sets of stochastic processes, one of which is composed by Ornstein-Uhlenbeck processes and the other being a general (non-linear) Langevin process. The method is able to distinguish between all stochastic process, retrieving their corresponding stochastic evolution equations. This framework is based on a recent approach for the analysis of multidimensional Langevin-type stochastic processes in the presence of strong measurement (or observational) noise, which is here extended to impose neither constraints nor parameters and extract all coefficients directly from the empirical data sets. Using synthetic data, it is shown that the method yields satisfactory results.

preprint2015arXiv

Uncovering the evolution of non-stationary stochastic variables: the example of asset volume-price fluctuations

We present a framework for describing the evolution of stochastic observables having a non-stationary distribution of values. The framework is applied to empirical volume-prices from assets traded at the New York stock exchange. Using Kullback-Leibler divergence we evaluate the best model out from four biparametric models standardly used in the context of financial data analysis. In our present data sets we conclude that the inverse $Γ$-distribution is a good model, particularly for the distribution tail of the largest volume-price fluctuations. Extracting the time-series of the corresponding parameter values we show that they evolve in time as stochastic variables themselves. For the particular case of the parameter controlling the volume-price distribution tail we are able to extract an Ornstein-Uhlenbeck equation which describes the fluctuations of the largest volume-prices observed in the data. Finally, we discuss how to bridge from the stochastic evolution of the distribution parameters to the stochastic evolution of the (non-stationary) observable and put our conclusions into perspective for other applications in geophysics and biology.

preprint2014arXiv

A thermostatistical approach to scale-free networks

We describe an ensemble of growing scale-free networks in an equilibrium framework, providing insight into why the exponent of empirical scale-free networks in nature is typically robust. In an analogy to thermostatistics, to describe the canonical and microcanonical ensembles, we introduce a functional, whose maximum corresponds to a scale-free configuration. We then identify the equivalents to energy, Zeroth-law, entropy and heat capacity for scale-free networks. Discussing the merging of scale-free networks, we also establish an exact relation to predict their final "equilibrium" degree exponent. All analytic results are complemented with Monte Carlo simulations. Our approach illustrates the possibility to apply the tools of equilibrium statistical physics to study the properties of growing networks, and it also supports the recent arguments on the complementarity between equilibrium and nonequilibrium systems.

preprint2014arXiv

Are credit ratings time-homogeneous and Markov?

We introduce a simple approach for testing the reliability of homogeneous generators and the Markov property of the stochastic processes underlying empirical time series of credit ratings. We analyze open access data provided by Moody's and show that the validity of these assumptions - existence of a homogeneous generator and Markovianity - is not always guaranteed. Our analysis is based on a comparison between empirical transition matrices aggregated over fixed time windows and candidate transition matrices generated from measurements taken over shorter periods. Ratings are widely used in credit risk, and are a key element in risk assessment; our results provide a tool for quantifying confidence in predictions extrapolated from these time series.

preprint2014arXiv

Daily pollution forecast using optimal meteorological data at synoptic and local scales

We present a simple framework to easily pre-select the most essential data for accurately forecasting the concentration of the pollutant PM$_{10}$, based on pollutants observations for the years 2002 until 2006 in the metropolitan region of Lisbon, Portugal. Starting from a broad panoply of different data sets collected at several meteorological stations, we apply a forward stepwise regression procedure that enables us not only to identify the most important variables for forecasting the pollutant but also to rank them in order of importance. We argue the importance of this variable ranking, showing that the ranking is very sensitive to the urban spot where measurements are taken. Having this pre-selection, we then present the potential of linear and non-linear neural network models when applied to the concentration of pollutant PM$_{10}$. Similarly to previous studies for other pollutants, our validation results show that non-linear models in average perform as well or worse as linear models for PM$_{10}$. Finally, we also address the influence of Circulation Weather Types, characterizing synoptic scale circulation patterns and the concentration of pollutants.

preprint2014arXiv

Fatigue Loads Estimation Through a Simple Stochastic Model

We propose a procedure to estimate the fatigue loads on wind turbines, based in a recent framework used for reconstructing data series of stochastic properties measured at wind turbines. Through a standard fatigue analysis, we show that it is possible to accurately estimate fatigue loads in any wind turbine within one wind park, using only the load measurements at one single turbine and the set of wind speed measurements. Our framework consists of deriving a stochastic differential equation that describes the evolution of the torque at one wind turbine driven by the wind speed. The stochastic equation is derived directly from the measurements and is afterwards used for predicting the fatigue loads at neighboring turbines. Such a framework could be used to mitigate the financial efforts usually necessary for placing measurement devices in all wind turbines within one wind farm. Finally, we also discuss the limitations and possible improvements of the proposed procedure.

preprint2014arXiv

From human mobility to renewable energies: Big data analysis to approach worldwide multiscale phenomena

We address and discuss recent trends in the analysis of big data sets, with the emphasis on studying multiscale phenomena. Applications of big data analysis in different scientific fields are described and two particular examples of multiscale phenomena are explored in more detail. The first one deals with wind power production at the scale of single wind turbines, the scale of entire wind farms and also at the scale of a whole country. Using open source data we show that the wind power production has an intermittent character at all those three scales, with implications for defining adequate strategies for stable energy production. The second example concerns the dynamics underlying human mobility, which presents different features at different scales. For that end, we analyze $12$-month data of the Eduroam database within Portuguese universities, and find that, at the smallest scales, typically within a set of a few adjacent buildings, the characteristic exponents of average displacements are different from the ones found at the scale of one country or one continent.

preprint2014arXiv

Modeling the functional network of primary intercellular Ca$^{2+}$ wave propagation in astrocytes and its application to study drug effects

We introduce a simple procedure of multivariate signal analysis to uncover the functional connectivity among cells composing a living tissue and describe how to apply it for extracting insight on the effect of drugs in the tissue. The procedure is based on the covariance matrix of time resolved activity signals. By determining the time-lag that maximizes covariance, one derives the weight of the corresponding connection between cells. Introducing simple constraints, it is possible to conclude whether pairs of cells are functionally connected and in which direction. After testing the method against synthetic data we apply it to study intercellular propagation of Ca$^{2+}$ waves in astrocytes following an external stimulus, with the aim of uncovering the functional cellular connectivity network. Our method proves to be particularly suited for this type of networking signal propagation where signals are pulse-like and have short time-delays, and is shown to be superior to standard methods, namely a multivariate Granger algorithm. Finally, based the statistical analysis of the connection weight distribution, we propose simple measures for assessing the impact of drugs on the functional connectivity between cells.

preprint2014arXiv

Optimal models of extreme volume-prices are time-dependent

We present evidence that the best model for empirical volume-price distributions is not always the same and it strongly depends in (i) the region of the volume-price spectrum that one wants to model and (ii) the period in time that is being modelled. To show these two features we analyze stocks of the New York stock market with four different models: Gamma, inverse-gamma, log-normal, and Weibull distributions. To evaluate the accuracy of each model we use standard relative deviations as well as the Kullback-Leibler distance and introduce an additional distance particularly suited to evaluate how accurate are the models for the distribution tails (large volume-price). Finally we put our findings in perspective and discuss how they can be extended to other situations in finance engineering.

preprint2014arXiv

Principal wind turbines for a conditional portfolio approach to wind farms

We introduce a measure for estimating the best risk-return relation of power production in wind farms within a given time-lag, conditioned to the velocity field. The velocity field is represented by a scalar that weighs the influence of the velocity at each wind turbine at present and previous time-steps for the present "state" of the wind field. The scalar measure introduced is a linear combination of the few turbines, that most influence the overall power production. This quantity is then used as the condition for computing a conditional expected return and corresponding risk associated to the future total power output.

preprint2014arXiv

Reconstructing the intermittent dynamics of the torque in wind turbines

We apply a framework introduced in the late nineties to analyze load measurements in off-shore wind energy converters (WEC). The framework is borrowed from statistical physics and properly adapted to the analysis of multivariate data comprising wind velocity, power production and torque measurements, taken at one single WEC. In particular, we assume that wind statistics drives the fluctuations of the torque produced in the wind turbine and show how to extract an evolution equation of the Langevin type for the torque driven by the wind velocity. It is known that the intermittent nature of the atmosphere, i.e. of the wind field, is transferred to the power production of a wind energy converter and consequently to the shaft torque. We show that the derived stochastic differential equation quantifies the dynamical coupling of the measured fluctuating properties as well as it reproduces the intermittency observed in the data. Finally, we discuss our approach in the light of turbine monitoring, a particular important issue in off-shore wind farms.

preprint2014arXiv

Stochastic Evolution of Stock Market Volume-Price Distributions

Using available data from the New York stock market (NYSM) we test four different bi-parametric models to fit the correspondent volume-price distributions at each $10$-minute lag: the Gamma distribution, the inverse Gamma distribution, the Weibull distribution and the log-normal distribution. The volume-price data, which measures market capitalization, appears to follow a specific statistical pattern, other than the evolution of prices measured in similar studies. We find that the inverse Gamma model gives a superior fit to the volume-price evolution than the other models. We then focus on the inverse Gamma distribution as a model for the NYSM data and analyze the evolution of the pair of distribution parameters as a stochastic process. Assuming that the evolution of these parameters is governed by coupled Langevin equations, we derive the corresponding drift and diffusion coefficients, which then provide insight for understanding the mechanisms underlying the evolution of the stock market.

preprint2013arXiv

Air quality prediction using optimal neural networks with stochastic variables

We apply recent methods in stochastic data analysis for discovering a set of few stochastic variables that represent the relevant information on a multivariate stochastic system, used as input for artificial neural networks models for air quality forecast. We show that using these derived variables as input variables for training the neural networks it is possible to significantly reduce the amount of input variables necessary for the neural network model, without considerably changing the predictive power of the model. The reduced set of variables including these derived variables is therefore proposed as optimal variable set for training neural networks models in forecasting geophysical and weather properties. Finally, we briefly discuss other possible applications of such optimized neural network models.

preprint2013arXiv

The bounds of heavy-tailed return distributions in evolving complex networks

We consider the evolution of scale-free networks according to preferential attachment schemes and show the conditions for which the exponent characterizing the degree distribution is bounded by upper and lower values. Our framework is an agent model, presented in the context of economic networks of trades, which shows the emergence of critical behavior. Starting from a brief discussion about the main features of the evolving network of trades, we show that the logarithmic return distributions have bounded heavy-tails, and the corresponding bounding exponent values can be derived. Finally, we discuss these findings in the context of model risk.

preprint2013arXiv

The dynamics of financial stability in complex networks

We address the problem of banking system resilience by applying off-equilibrium statistical physics to a system of particles, representing the economic agents, modelled according to the theoretical foundation of the current banking regulation, the so called Merton-Vasicek model. Economic agents are attracted to each other to exchange `economic energy', forming a network of trades. When the capital level of one economic agent drops below a minimum, the economic agent becomes insolvent. The insolvency of one single economic agent affects the economic energy of all its neighbours which thus become susceptible to insolvency, being able to trigger a chain of insolvencies (avalanche). We show that the distribution of avalanche sizes follows a power-law whose exponent depends on the minimum capital level. Furthermore, we present evidence that under an increase in the minimum capital level, large crashes will be avoided only if one assumes that agents will accept a drop in business levels, while keeping their trading attitudes and policies unchanged. The alternative assumption, that agents will try to restore their business levels, may lead to the unexpected consequence that large crises occur with higher probability.

preprint2013arXiv

Uncovering wind turbine properties through two-dimensional stochastic modeling of wind dynamics

Using a method for stochastic data analysis, borrowed from statistical physics, we analyze synthetic data from a Markov chain model that reproduces measurements of wind speed and power production in a wind park in Portugal. We first show that our analysis retrieves indeed the power performance curve, which yields the relationship between wind speed and power production and we discuss how this procedure can be extended for extracting unknown functional relationships between pairs of physical variables in general. Second, we show how specific features, such as the rated speed of the wind turbine or the descriptive wind speed statistics, can be related with the equations describing the evolution of power production and wind speed at single wind turbines.

preprint2012arXiv

Extracting strong measurement noise from stochastic series: applications to empirical data

It is a big challenge in the analysis of experimental data to disentangle the unavoidable measurement noise from the intrinsic dynamical noise. Here we present a general operational method to extract measurement noise from stochastic time series, even in the case when the amplitudes of measurement noise and uncontaminated signal are of the same order of magnitude. Our approach is based on a recently developed method for a nonparametric reconstruction of Langevin processes. Minimizing a proper non-negative function the procedure is able to correctly extract strong measurement noise and to estimate drift and diffusion coefficients in the Langevin equation describing the evolution of the original uncorrupted signal. As input, the algorithm uses only the two first conditional moments extracted directly from the stochastic series and is therefore suitable for a broad panoply of different signals. To demonstrate the power of the method we apply the algorithm to synthetic as well as climatological measurement data, namely the daily North Atlantic Oscillation index, shedding new light on the discussion of the nature of its underlying physical processes.

preprint2012arXiv

Heavy-tails in economic data: fundamental assumptions, modelling and analysis

The study of heavy-tailed distributions in economic and financial systems has been widely addressed since financial time series has become a research subject.After the eighties, several "highly improbable" market drops were observed (e.g. the 1987 stock market drop known as "Black Monday" and on even more recent ones, already in the 21st century) that produce heavy losses that were unexplainable in a GN environment. The losses incurred in these large market drop events did not change significantly the market practices or the way regulation is done but drove some attention back to the study of heavy-tails and their underlying mechanisms. Some recent findings in these context is the scope of this manuscript.

preprint2012arXiv

Searching for optimal variables in real multivariate stochastic data

By implementing a recent technique for the determination of stochastic eigendirections of two coupled stochastic variables, we investigate the evolution of fluctuations of NO2 concentrations at two monitoring stations in the city of Lisbon, Portugal. We analyze the stochastic part of the measurements recorded at the monitoring stations by means of a method where the two concentrations are considered as stochastic variables evolving according to a system of coupled stochastic differential equations. Analysis of their structure allows for transforming the set of measured variables to a set of derived variables, one of them with reduced stochasticity. For the specific case of NO2 concentration measures, the set of derived variables are well approximated by a global rotation of the original set of measured variables. We conclude that the stochastic sources at each station are independent from each other and typically have amplitudes of the order of the deterministic contributions. Such findings show significant limitations when predicting such quantities. Still, we briefly discuss how predictive power can be increased in general in the light of our methods.

preprint2012arXiv

Self-organized criticality in a network of economic agents with finite consumption

We introduce a simple model for addressing the controversy in the study of financial systems, sometimes taken as brownian-like processes and other as critical systems with fluctuations of arbitrary magnitude. The model considers a collection of economical agents which establish trade connections among them according to basic economical principles properly translated into physical properties and interaction. With our model we are able to reproduce the evolution of macroscopic quantities (indices) and to correctly retrieve the common exponent value characterizing several indices in financial markets, relating it to the underlying topology of connections.

preprint2012arXiv

The spread of gossip in American schools

We study different mechanisms of gossip propagation on several network topologies and introduce a new network property, the "spread factor", describing the fraction of neighbors that get to know the gossip. We postulate that for scale-free networks the spreading time grows logarithmically with the degree of the victim and prove this statement for the case of the Apollonian network. Applying our concepts to real data from an American school survey, we confirm the logarithmic law and disclose that there exists an ideal number of acquaintances minimizing the fraction attained by the gossip. The similarity between the school survey and scale-free networks remains even for cases when gossip propagation only occurs with some probability $q<1$. The spreading times follow an exponential distribution that can also be calculated analytically for the Apollonian network. When gossip also spreads through strangers the situation changes substantially: the spreading time becomes a constant and there exists no ideal degree of connectivity anymore.

preprint2011arXiv

Deposition of general ellipsoidal particles

We present a systematic overview of granular deposits composed of ellipsoidal particles with different particle shapes and size polydispersities. We study the density and anisotropy of such deposits as functions of size polydispersity and two shape parameters that fully describe the shape of a general ellipsoid. Our results show that, while shape influences significantly the macroscopic properties of the deposits, polydispersity plays apparently a secondary role. The density attains a maximum for a particular family of non-symmetrical ellipsoids, larger than the density observed for prolate or oblate ellipsoids. As for anisotropy measures, the contact forces show are increasingly preferred along the vertical direction as the shape of the particles deviates for a sphere. The deposits are constructed by means of an efficient molecular dynamics method, where the contact forces are efficiently and accurately computed. The main results are discussed in the light of applications for porous media models and sedimentation processes.

preprint2011arXiv

Size distribution and structure of Barchan dune fields

Barchans are isolated mobile dunes often organized in large dune fields. Dune fields seem to present a characteristic dune size and spacing, which suggests a cooperative behavior based on dune interaction. In Duran et al. (2009), we propose that the redistribution of sand by collisions between dunes is a key element for the stability and size selection of barchan dune fields. This approach was based on a mean-field model ignoring the spatial distribution of dune fields. Here, we present a simplified dune field model that includes the spatial evolution of individual dunes as well as their interaction through sand exchange and binary collisions. As a result, the dune field evolves towards a steady state that depends on the boundary conditions. Comparing our results with measurements of Moroccan dune fields, we find that the simulated fields have the same dune size distribution as in real fields but fail to reproduce their homogeneity along the wind direction.

preprint2010arXiv

Superdiffusion of massive particles induced by multi-scale velocity fields

We study drag-induced diffusion of massive particles in scale-free velocity fields, where superdiffusive behavior emerges due to the scale-free size distribution of the vortices of the underlying velocity field. The results show qualitative resemblance to what is observed in fluid systems, namely the diffusive exponent for the mean square separation of pairs of particles and the preferential concentration of the particles, both as a function of the response time.

preprint2008arXiv

Avalanches in anisotropic sheared granular media

We study the influence of particle shape anisotropy on the occurrence of avalanches in sheared granular media. We use molecular dynamic simulations to calculate the relative movement of two tectonic plates. % with transform boundaries. Our model considers irregular polygonal particles constituting the material within the shear zone. We find that the magnitude of the avalanches is approximately independent on particle shape and in good agreement with the Gutenberg-Richter law, but the aftershock sequences are strongly influenced by the particle anisotropy yielding variations on the exponent characterizing the empirical Omori's law. Our findings enable one to identify the presence of anisotropic particles at the macro-mechanical level only by observing the avalanche sequences of real faults. In addition, we calculate the probability of occurrence of an avalanche for given values of stiffness or frictional strength and observe also a significant influence of the particle anisotropy.

preprint2008arXiv

The dune size distribution and scaling relations of barchan dune fields

Barchan dunes emerge as a collective phenomena involving the generation of thousands of them in so called barchan dune fields. By measuring the size and position of dunes in Moroccan barchan dune fields, we find that these dunes tend to distribute uniformly in space and follow an unique size distribution function. We introduce an analyticalmean-field approach to show that this empirical size distribution emerges from the interplay of dune collisions and sand flux balance, the two simplest mechanisms for size selection. The analytical model also predicts a scaling relation between the fundamental macroscopic properties characterizing a dune field, namely the inter-dune spacing and the first and second moments of the dune size distribution.

preprint2007arXiv

Modeling slow deformation of polygonal particles using DEM

We introduce two improvements in the numerical scheme to simulate collision and slow shearing of irregular particles. First, we propose an alternative approach based on simple relations to compute the frictional contact forces. The approach improves efficiency and accuracy of the Discrete Element Method (DEM) when modeling the dynamics of the granular packing. We determine the proper upper limit for the integration step in the standard numerical scheme using a wide range of material parameters. To this end, we study the kinetic energy decay in a stress controlled test between two particles. Second, we show that the usual way of defining the contact plane between two polygonal particles is, in general, not unique which leads to discontinuities in the direction of the contact plane while particles move. To solve this drawback, we introduce an accurate definition for the contact plane based on the shape of the overlap area between touching particles, which evolves continuously in time.

preprint2007arXiv

New approaches to model and study social networks

We describe and develop three recent novelties in network research which are particularly useful for studying social systems. The first one concerns the discovery of some basic dynamical laws that enable the emergence of the fundamental features observed in social networks, namely the nontrivial clustering properties, the existence of positive degree correlations and the subdivision into communities. To reproduce all these features we describe a simple model of mobile colliding agents, whose collisions define the connections between the agents which are the nodes in the underlying network, and develop some analytical considerations. The second point addresses the particular feature of clustering and its relationship with global network measures, namely with the distribution of the size of cycles in the network. Since in social bipartite networks it is not possible to measure the clustering from standard procedures, we propose an alternative clustering coefficient that can be used to extract an improved normalized cycle distribution in any network. Finally, the third point addresses dynamical processes occurring on networks, namely when studying the propagation of information in them. In particular, we focus on the particular features of gossip propagation which impose some restrictions in the propagation rules. To this end we introduce a quantity, the spread factor, which measures the average maximal fraction of nearest neighbors which get in contact with the gossip, and find the striking result that there is an optimal non-trivial number of friends for which the spread factor is minimized, decreasing the danger of being gossiped.

preprint2007arXiv

Numerical Improvement of the Discrete Element Method applied to Shear of Granular Media

We present a detailed analysis of the bounds on the integration step in Discrete Element Method (DEM) for simulating collisions and shearing of granular assemblies. We show that, in the numerical scheme, the upper limit for the integration step, usually taken from the average time $t_c$ of one contact, is in fact not sufficiently small to guarantee numerical convergence of the system during relaxation. In particular, we study in detail how the kinetic energy decays during the relaxation stage and compute the correct upper limits for the integration step, which are significantly smaller than the ones commonly used. In addition, we introduce an alternative approach, based on simple relations to compute the frictional forces, that converges even for integration steps above the upper limit.

preprint2007arXiv

Obtaining the size distribution of fault gouges with polydisperse bearings

We generalize the recent study of random space-filling bearings to a more realistic situation, where the spacing offset varies randomly during the space-filling procedure, and show that it reproduces well the size-distributions observed in recent studies of real fault gouges. In particular, we show that the fractal dimensions of random polydisperse bearings sweep predominantly the low range of values in the spectrum of fractal dimensions observed along real faults, which strengthen the evidence that polydisperse bearings may explain the occurrence of seismic gaps in nature. In addition, the influence of different distributions for the offset is studied and we find that the uniform distribution is the best choice for reproducing the size-distribution of fault gouges.

preprint2007arXiv

Spreading gossip in social networks

We study a simple model of information propagation in social networks, where two quantities are introduced: the spread factor, which measures the average maximal fraction of neighbors of a given node that interchange information among each other, and the spreading time needed for the information to reach such fraction of nodes. When the information refers to a particular node at which both quantities are measured, the model can be taken as a model for gossip propagation. In this context, we apply the model to real empirical networks of social acquaintances and compare the underlying spreading dynamics with different types of scale-free and small-world networks. We find that the number of friendship connections strongly influences the probability of being gossiped. Finally, we discuss how the spread factor is able to be applied to other situations.

preprint2006arXiv

Networks based on collisions among mobile agents

We investigate in detail a recent model of colliding mobile agents [Phys. Rev. Lett.~96, 088702], used as an alternative approach to construct evolving networks of interactions formed by the collisions governed by suitable dynamical rules. The system of mobile agents evolves towards a quasi-stationary state which is, apart small fluctuations, well characterized by the density of the system and the residence time of the agents. The residence time defines a collision rate and by varying the collision rate, the system percolates at a critical value, with the emergence of a giant cluster whose critical exponents are the ones of two-dimensional percolation. Further, the degree and clustering coefficient distributions and the average path length show that the network associated with such a system presents non-trivial features which, depending on the collision rule, enables one not only to recover the main properties of standard networks, such as exponential, random and scale-free networks, but also to obtain other topological structures. Namely, we show a specific example where the obtained structure has topological features which characterize accurately the structure and evolution of social networks in different contexts, ranging from networks of acquaintances to networks of sexual contacts.

preprint2005arXiv

Coherence in complex networks of oscillators

We study fully synchronized (coherent) states in complex networks of chaotic oscillators, reviewing the analytical approach of determining the stability conditions for synchronizability and comparing them with numerical criteria. As an example, we present detailed results for networks of chaotic logistic maps having three different scale-free topologies: random scale-free topology, deterministic pseudo-fractal scale-free network and Apollonian network. For random scale-free topology we find that the lower boundary of the synchronizability region scales approximately as $k^{-μ}$, where $k$ is the outgoing connectivity and $μ$ depends on the local nonlinearity. For deterministic scale-free networks coherence is observed only when the coupling is heterogeneous, namely when it is proportional to some power of the neighbor connectivity. In all cases, stability conditions are determined from the eigenvalue spectrum of the Laplacian matrix and agree well with numerical results based on histograms of coherent states in parameter space. Additionally, we show that almost everywhere in the synchronizability region the basin of attraction of coherent states fills the entire phase space, and that the transition to coherence is of first-order.

preprint2005arXiv

Cycles and clustering in bipartite networks

We investigate the clustering ability in bipartite networks where cycles of size three are absent and therefore the standard definition of clustering coefficient cannot be used. Instead, we use another coefficient given by the fraction of cycles with size four, showing that both coefficients yield the same clustering properties. The new coefficient is computed for two networks of sexual contacts, one monopartite and another bipartite. In both cases the clustering ability is similar. Furthermore, combining both clustering coefficients we deduce an expression for estimating cycles of larger size, which improves previous estimations and is suitable for either monopartite and multipartite networks.

preprint2005arXiv

Impact of bistability in the synchronization of chaotic maps with delayed coupling

We investigate the impact of bistability in the emergence of synchronization in networks of chaotic maps with delayed coupling. The existence of a single finite attractor of the uncoupled map is found to be responsible for the emergence of synchronization. No synchronization is observed when the local dynamics has two competing chaotic attractors whose orbits are dense on the same interval. This result is robust for regular networks with variable ranges of interaction and for more complex topologies.