Catalog footprint

What is connected

63works
30topics
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

63 published item(s)

preprint2022arXiv

Deep learning estimation of complex reverberant wave fields by a programmable metasurface

Electromagnetic environments are becoming increasingly complex and congested, creating a growing challenge for systems that rely on electromagnetic waves for communication, sensing, or imaging, particularly in reverberating environments. The use of programmable metasurfaces provides a potential means of directing waves to optimize wireless channels on-demand, ensuring reliable operation and protecting sensitive electronic components. Here we introduce a technique that combines a deep learning network with a binary programmable metasurface to shape waves in complex reverberant electromagnetic environments, in particular ones where there is no direct line of sight. We applied this technique for wavefront reconstruction and control, and accurately determined metasurface configurations based on measured system scattering responses in a chaotic microwave cavity. The state of the metasurface that realizes desired electromagnetic wave field distribution properties was successfully determined even in cases previously unseen by the deep learning algorithm. Our technique is enabled by the reverberant nature of the cavity, and is effective with a metasurface that covers only $\sim$1.5\% of the total cavity surface area.

preprint2022arXiv

Short-wavelength Reverberant Wave Systems for Physical Realization of Reservoir Computing

Machine learning (ML) has found widespread application over a broad range of important tasks. To enhance ML performance, researchers have investigated computational architectures whose physical implementations promise compactness, high-speed execution, physical robustness, and low energy cost. Here, we experimentally demonstrate an approach that uses the high sensitivity of reverberant short wavelength waves for physical realization and enhancement of computational power of a type of ML known as reservoir computing (RC). The potential computation power of RC systems increases with their effective size. We here exploit the intrinsic property of short wavelength reverberant wave sensitivity to perturbations to expand the effective size of the RC system by means of spatial and spectral perturbations. Working in the microwave regime, this scheme is tested experimentally on different ML tasks. Our results indicate the general applicability of reverberant wave-based implementations of RC and of our effective reservoir size expansion technique

preprint2021arXiv

Parallel Machine Learning for Forecasting the Dynamics of Complex Networks

Forecasting the dynamics of large complex networks from previous time-series data is important in a wide range of contexts. Here we present a machine learning scheme for this task using a parallel architecture that mimics the topology of the network of interest. We demonstrate the utility and scalability of our method implemented using reservoir computing on a chaotic network of oscillators. Two levels of prior knowledge are considered: (i) the network links are known; and (ii) the network links are unknown and inferred via a data-driven approach to approximately optimize prediction.

preprint2020arXiv

Backpropagation Algorithms and Reservoir Computing in Recurrent Neural Networks for the Forecasting of Complex Spatiotemporal Dynamics

We examine the efficiency of Recurrent Neural Networks in forecasting the spatiotemporal dynamics of high dimensional and reduced order complex systems using Reservoir Computing (RC) and Backpropagation through time (BPTT) for gated network architectures. We highlight advantages and limitations of each method and discuss their implementation for parallel computing architectures. We quantify the relative prediction accuracy of these algorithms for the longterm forecasting of chaotic systems using as benchmarks the Lorenz-96 and the Kuramoto-Sivashinsky (KS) equations. We find that, when the full state dynamics are available for training, RC outperforms BPTT approaches in terms of predictive performance and in capturing of the long-term statistics, while at the same time requiring much less training time. However, in the case of reduced order data, large scale RC models can be unstable and more likely than the BPTT algorithms to diverge. In contrast, RNNs trained via BPTT show superior forecasting abilities and capture well the dynamics of reduced order systems. Furthermore, the present study quantifies for the first time the Lyapunov Spectrum of the KS equation with BPTT, achieving similar accuracy as RC. This study establishes that RNNs are a potent computational framework for the learning and forecasting of complex spatiotemporal systems.

preprint2020arXiv

Combining Machine Learning with Knowledge-Based Modeling for Scalable Forecasting and Subgrid-Scale Closure of Large, Complex, Spatiotemporal Systems

We consider the commonly encountered situation (e.g., in weather forecasting) where the goal is to predict the time evolution of a large, spatiotemporally chaotic dynamical system when we have access to both time series data of previous system states and an imperfect model of the full system dynamics. Specifically, we attempt to utilize machine learning as the essential tool for integrating the use of past data into predictions. In order to facilitate scalability to the common scenario of interest where the spatiotemporally chaotic system is very large and complex, we propose combining two approaches:(i) a parallel machine learning prediction scheme; and (ii) a hybrid technique, for a composite prediction system composed of a knowledge-based component and a machine-learning-based component. We demonstrate that not only can this method combining (i) and (ii) be scaled to give excellent performance for very large systems, but also that the length of time series data needed to train our multiple, parallel machine learning components is dramatically less than that necessary without parallelization. Furthermore, considering cases where computational realization of the knowledge-based component does not resolve subgrid-scale processes, our scheme is able to use training data to incorporate the effect of the unresolved short-scale dynamics upon the resolved longer-scale dynamics ("subgrid-scale closure").

preprint2020arXiv

Critical Network Cascades with Re-excitable nodes: Why tree-like approximations usually work, when they breakdown, and how to correct them

Network science is a rapidly expanding field, with a large and growing body of work on network-based dynamical processes. Most theoretical results in this area rely on the so-called \emph{locally tree-like approximation}. This is, however, usually an `uncontrolled' approximation, in the sense that the magnitudes of the error are typically unknown, although numerical results show that this error is often surprisingly small. In this paper, we place this approximation on more rigorous footing by calculating the magnitude of deviations away from tree-based theories in the context of discrete-time critical network cascades with re-excitable nodes. We discuss the conditions under which tree-like approximations give good results for calculating network criticality, and also explain the reasons for deviation from this approximation, in terms of the density of certain kinds of network motifs. Using this understanding, we derive results for network criticality that apply to general networks that explicitly do not satisfy the locally tree-like approximation. In particular, we focus on the bi-parallel motif, the smallest motif relevant to the failure of a tree-based theory in this context, and we derive the corrections due to such motifs on the conditions for criticality. We verify our claims on computer-generated networks, and we confirm that our theory accurately predicts the observed deviations from criticality. Using our theory, we explain why numerical simulations often show that deviations from a tree-based theory are surprisingly small. More specifically, we show that these deviations are negligible for networks whose average degree is even modestly large compared to one, justifying why tree-based theories appear to work well for most real-world networks.

preprint2020arXiv

Efficient Statistical Model for Predicting Electromagnetic Wave Distribution in Coupled Enclosures

The Random Coupling Model (RCM) has been successfully applied to predicting the statistics of currents and voltages at ports in complex electromagnetic (EM) enclosures operating in the short wavelength limit. Recent studies have extended the RCM to systems of multi-mode aperture-coupled enclosures. However, as the size (as measured in wavelengths) of a coupling aperture grows, the coupling matrix used in the RCM increases as well, and the computation becomes more complex and time-consuming. A simple Power Balance Model (PWB) can provide fast predictions for the \textit{averaged} power density of waves inside electrically-large systems for a wide range of cavity and coupling scenarios. However, the important interference induced fluctuations of the wavefield retained in the RCM is absent in PWB. Here we aim to combine the best aspects of each model to create a hybrid treatment and study the EM fields in coupled enclosure systems. The proposed hybrid approach provides both mean and fluctuation information of the EM fields without the full computational complexity of coupled-cavity RCM. We compare the hybrid model predictions with experiments on linear cascades of over-moded cavities. We find good agreement over a set of different loss parameters and for different coupling strengths between cavities. The range of validity and applicability of the hybrid method are tested and discussed.

preprint2020arXiv

Wavefront Shaping with a Tunable Metasurface: Creating Coldspots and Coherent Perfect Absorption at Arbitrary Frequencies

Modern electronic systems operate in complex electromagnetic environments and must handle noise and unwanted coupling. The capability to isolate or reject unwanted signals for mitigating vulnerabilities is critical in any practical application. In this work, we describe the use of a binary programmable metasurface to (i) control the spatial degrees of freedom for waves propagating inside an electromagnetic cavity and demonstrate the ability to create nulls in the transmission coefficient between selected ports; and (ii) create the conditions for coherent perfect absorption. Both objectives are performed at arbitrary frequencies. In the first case a novel and effective optimization algorithm is presented that selectively generates coldspots over a single frequency band or simultaneously over multiple frequency bands. We show that this algorithm is successful with multiple input port configurations and varying optimization bandwidths. In the second case we establish how this technique can be used to establish a multi-port coherent perfect absorption state for the cavity.

preprint2019arXiv

Dynamic regulation of resource transport induces criticality in interdependent networks of excitable units

Various functions of a network of excitable units can be enhanced if the network is in the `critical regime', where excitations are, on average, neither damped nor amplified. An important question is how can such networks self-organize to operate in the critical regime. Previously it was shown that regulation via resource transport on a secondary network can robustly maintain the primary network dynamics in a balanced state where activity doesn't grow or decay. Here we show that this inter-network regulation process robustly produces a power-law distribution of activity avalanches, as observed in experiments, over ranges of model parameters spanning orders of magnitude. We also show that the resource transport over the secondary network protects the system against the destabilizing effect of local variations in parameters and heterogeneity in network structure. For homogeneous networks, we derive a reduced 3-dimensional map which reproduces the behavior of the full system.

preprint2019arXiv

Separation of Chaotic Signals by Reservoir Computing

We demonstrate the utility of machine learning in the separation of superimposed chaotic signals using a technique called Reservoir Computing. We assume no knowledge of the dynamical equations that produce the signals, and require only training data consisting of finite time samples of the component signals. We test our method on signals that are formed as linear combinations of signals from two Lorenz systems with different parameters. Comparing our nonlinear method with the optimal linear solution to the separation problem, the Wiener filter, we find that our method significantly outperforms the Wiener filter in all the scenarios we study. Furthermore, this difference is particularly striking when the component signals have similar frequency spectra. Indeed, our method works well when the component frequency spectra are indistinguishable - a case where a Wiener filter performs essentially no separation.

preprint2019arXiv

Wave scattering properties of multiple weakly-coupled complex systems

The statistics of scattering of waves inside single ray-chaotic enclosures have been successfully described by the Random Coupling Model (RCM). We expand the RCM to systems consisting of multiple complex ray-chaotic enclosures with variable coupling scenarios. The statistical properties of the model-generated quantities are tested against measured data of electrically large multi-cavity systems of various designs. The statistics of model-generated trans-impedance and induced voltages on a load impedance agree well with the experimental results. The RCM coupled chaotic enclosure model is general and can be applied to other physical systems including coupled quantum dots, disordered nanowires, and short-wavelength electromagnetic propagation through rooms in buildings, aircraft and ships.

preprint2016arXiv

Focusing Waves at Arbitrary Locations in a Ray-Chaotic Enclosure Using Time-Reversed Synthetic Sonas

Time reversal methods are widely used to achieve wave focusing in acoustics and electromagnetics. A typical time reversal experiment requires that a transmitter be initially present at the target focusing point, which limits the application of this technique. In this paper, we propose a method to focus waves at an arbitary location inside a complex enclosure using a numerically calculated wave excitation signal. We use a semi-classical ray algorithm to calculate the signal that would be received at a transceiver port resulting from the injection of a short pulse at the desired target location. The time-reversed version of this signal is then injected into the transceiver port and an approximate reconstruction of the short pulse is created at the target. The quaility of the pulse reconstruction is quantified in three different ways and the values of these metrics are predicted by the statistics of the scattering-parameter $|S_{21}|^2$ between the transceiver and target points in the enclosure. We experimentally demonstrate the method using a flat microwave billiard and quantify the reconstruction quality as a function of enclosure loss, port coupling and other considerations.

preprint2016arXiv

Intermodulation in Nonlinear SQUID Metamaterials: Experiment and Theory

The response of nonlinear metamaterials and superconducting electronics to two-tone excitation is critical for understanding their use as low-noise amplifiers and tunable filters. A new setting for such studies is that of metamaterials made of radio frequency superconducting quantum interference devices (rf-SQUIDs). The two-tone response of self-resonant rf-SQUID meta-atoms and metamaterials is studied here via intermodulation (IM) measurement over a broad range of tone frequencies and tone powers. A sharp onset followed by a surprising strongly suppressed IM region near the resonance is observed. Using a two time scale analysis technique, we present an analytical theory that successfully explains our experimental observations. The theory predicts that the IM can be manipulated with tone power, center frequency, frequency difference between the two tones, and temperature. This quantitative understanding potentially allows for the design of rf-SQUID metamaterials with either very low or very high IM response.

preprint2016arXiv

Metabolite transport through glial networks stabilizes the dynamics of learning

Learning and memory are acquired through long-lasting changes in synapses. In the simplest models, such synaptic potentiation typically leads to runaway excitation, but in reality there must exist processes that robustly preserve overall stability of the neural system dynamics. How is this accomplished? Various approaches to this basic question have been considered. Here we propose a particularly compelling and natural mechanism for preserving stability of learning neural systems. This mechanism is based on the global processes by which metabolic resources are distributed to the neurons by glial cells. Specifically, we introduce and study a model comprised of two interacting networks: a model neural network interconnected by synapses which undergo spike-timing dependent plasticity (STDP); and a model glial network interconnected by gap junctions which diffusively transport metabolic resources among the glia and, ultimately, to neural synapses where they are consumed. Our main result is that the biophysical constraints imposed by diffusive transport of metabolic resources through the glial network can prevent runaway growth of synaptic strength, both during ongoing activity and during learning. Our findings suggest a previously unappreciated role for glial transport of metabolites in the feedback control stabilization of neural network dynamics during learning.

preprint2015arXiv

A statistical model for the excitation of cavities through apertures

In this paper, a statistical model for the coupling of electromagnetic radiation into enclosures through apertures is presented. The model gives a unified picture bridging deterministic theories of aperture radiation, and statistical models necessary for capturing the properties of irregular shaped enclosures. A Monte Carlo technique based on random matrix theory is used to predict and study the power transmitted through the aperture into the enclosure. Universal behavior of the net power entering the aperture is found. Results are of interest for predicting the coupling of external radiation through openings in irregular enclosures and reverberation chambers.

preprint2015arXiv

Defining Chaos

In this paper we propose, discuss and illustrate a computationally feasible definition of chaos which can be applied very generally to situations that are commonly encountered, including attractors, repellers and non-periodically forced systems. This definition is based on an entropy-like quantity, which we call "expansion entropy", and we define chaos as occurring when this quantity is positive. We relate and compare expansion entropy to the well-known concept of topological entropy, to which it is equivalent under appropriate conditions. We also present example illustrations, discuss computational implementations, and point out issues arising from attempts at giving definitions of chaos that are not entropy-based.

preprint2015arXiv

Frequency assortativity can induce chaos in oscillator networks

We investigate the effect of preferentially connecting oscillators with similar frequency to each other in networks of coupled phase oscillators (i.e., frequency assortativity). Using the network Kuramoto model as an example, we find that frequency assortativity can induce chaos in the macroscopic dynamics. By applying a mean-field approximation in combination with the dimension reduction method of Ott and Antonsen, we show that the dynamics can be described by a low dimensional system of equations. We use the reduced system to characterize the macroscopic chaos using Lyapunov exponents, bifurcation diagrams, and time-delay embeddings. Finally, we show that the emergence of chaos stems from the formation of multiple groups of synchronized oscillators, i.e., meta-oscillators.

preprint2014arXiv

Controlling systems that drift through a tipping point

Slow parameter drift is common in many systems (e.g., the amount of greenhouse gases in the terrestrial atmosphere is increasing). In such situations, the attractor on which the system trajectory lies can be destroyed, and the trajectory will then go to another attractor of the system. We consider the case where there are more than one of these possible final attractors, and we ask whether we can control the outcome (i.e., the attractor that ultimately captures the trajectory) using only small controlling perturbations. Specifically, we consider the problem of controlling a noisy system whose parameter slowly drifts through a saddle-node bifurcation taking place on a fractal boundary between the basins of multiple attractors. We show that, when the noise level is low, a small perturbation of size comparable to the noise amplitude applied at a single point in time can ensure that the system will evolve toward a target attracting state with high probability. For a range of noise levels, we find that the minimum size of perturbation required for control is much smaller within a time period that starts some time after the bifurcation, providing a "window of opportunity" for driving the system toward a desirable state. We refer to this procedure as tipping point control.

preprint2014arXiv

Dynamical Transitions in Large Systems of Mean Field-Coupled Landau-Stuart Oscillators: Extensive Chaos and Clumped States

In this paper, we study dynamical systems in which a large number $N$ of identical Landau-Stuart oscillators are globally coupled via a mean-field. Previously, it has been observed that this type of system can exhibit a variety of different dynamical behaviors including clumped states in which each oscillator is in one of a small number of groups for which all oscillators in each group have the same state which is different from group to group, as well as situations in which all oscillators have different states and the macroscopic dynamics of the mean field is chaotic. We argue that this second type of behavior is $^{\backprime}$extensive$^{\prime}$ in the sense that the chaotic attractor in the full phase space of the system has a fractal dimension that scales linearly with $N$ and that the number of positive Lyapunov exponents of the attractor also scales with linearly $N$. An important focus of this paper is the transition between clumped states and extensive chaos as the system is subjected to slow adiabatic parameter change. We observe explosive (i.e., discontinuous) transitions between the clumped states (which correspond to low dimensional dynamics) and the extensively chaotic states. Furthermore, examining the clumped state, as the system approaches the explosive transition to extensive chaos, we find that the oscillator population distribution between the clumps continually evolves so that the clumped state is always marginally stable. This behavior is used to reveal the mechanism of the explosive transition. We also apply the Kaplan-Yorke formula to study the fractal structure of the extensively chaotic attractors.

preprint2014arXiv

Inhibition causes ceaseless dynamics in networks of excitable nodes

The collective dynamics of a network of excitable nodes changes dramatically when inhibitory nodes are introduced. We consider inhibitory nodes which may be activated just like excitatory nodes but, upon activating, decrease the probability of activation of network neighbors. We show that, although the direct effect of inhibitory nodes is to decrease activity, the collective dynamics becomes self-sustaining. We explain this counterintuitive result by defining and analyzing a "branching function" which may be thought of as an activity-dependent branching ratio. The shape of the branching function implies that for a range of global coupling parameters dynamics are self-sustaining. Within the self-sustaining region of parameter space lies a critical line along which dynamics take the form of avalanches with universal scaling of size and duration, embedded in ceaseless timeseries of activity. Our analyses, confirmed by numerical simulation, suggest that inhibition may play a counterintuitive role in excitable networks.

preprint2014arXiv

Mean field theory of assortative networks of phase oscillators

Employing the Kuramoto model as an illustrative example, we show how the use of the mean field approximation can be applied to large networks of phase oscillators with assortativity. We then use the ansatz of Ott and Antonsen [Chaos 19, 037113 (2008)] to reduce the mean field kinetic equations to a system of ordinary differential equations. The resulting formulation is illustrated by application to a network Kuramoto problem with degree assortativity and correlation between the node degrees and the natural oscillation frequencies. Good agreement is found between the solutions of the reduced set of ordinary differential equations obtained from our theory and full simulations of the system. These results highlight the ability of our method to capture all the phase transitions (bifurcations) and system attractors. One interesting result is that degree assortativity can induce transitions from a steady macroscopic state to a temporally oscillating macroscopic state through both (presumed) Hopf and SNIPER (saddle-node, infinite period) bifurcations. Possible use of these techniques to a broad class of phase oscillator network problems is discussed.

preprint2014arXiv

Mitigating the effect of non-uniform loss on time reversal mirrors

Time reversal mirrors work perfectly only for lossless wave propagation. Here, the performance of time-reversal mirrors is quantitatively defined, and the adverse effect of dissipation on their performance is investigated. An application of the technique of exponential amplification is proposed to overcome the effect of dissipation in the case of uniform loss distributions, and, to some extent, in the case of non-uniform loss distributions. A numerical model of a star graph was employed to test the applicability of this technique on realizations with various random spatial distributions of loss. A subset of the numerical results are also verified by experimental results from an electromagnetic time-reversal mirror. The exponential amplification technique should improve the performance of emerging technologies based on time-reversed wave propagation such as directed communication and wireless power transfer.

preprint2014arXiv

The Impact of Imperfect Information on Network Attack

This paper explores the effectiveness of network attack when the attacker has imperfect information about the network. For Erdős-Rényi networks, we observe that dynamical importance and betweenness centrality-based attacks are surprisingly robust to the presence of a moderate amount of imperfect information and are more effective compared with simpler degree-based attacks even at moderate levels of network information error. In contrast, for scale-free networks the effectiveness of attack is much less degraded by a moderate level of information error. Furthermore, in the Erdőos-Rényi case the effectiveness of network attack is much more degraded by missing links as compared with the same number of false links.

preprint2014arXiv

The Scaling of Chaos vs Periodicity: How Certain is it that an Attractor is Chaotic?

The character of the time-asymptotic evolution of physical systems can have complex, singular behavior with variation of a system parameter, particularly when chaos is involved. A perturbation of the parameter by a small amount $ε$ can convert an attractor from chaotic to non-chaotic or vice-versa. We call a parameter value where this can happen $ε$-uncertain. The probability that a random choice of the parameter is $ε$-uncertain commonly scales like a power law in $ε$. Surprisingly, two seemingly similar ways of defining this scaling, both of physical interest, yield different numerical values for the scaling exponent. We show why this happens and present a quantitative analysis of this phenomenon.

preprint2013arXiv

Impedance and Scattering Variance Ratios of Complicated Wave Scattering Systems in the Low Loss Regime

Random matrix theory (RMT) successfully predicts universal statistical properties of complicated wave scattering systems in the semiclassical limit, while the random coupling model offers a complete statistical model with a simple additive formula in terms of impedance to combine the predictions of RMT and nonuniversal system-specific features. The statistics of measured wave properties generally have nonuniversal features. However, ratios of the variances of elements of the impedance matrix are predicted to be independent of such nonuniversal features and thus should be universal functions of the overall system loss. In contrast with impedance variance ratios, scattering variance ratios depends on nonuniversal features unless the system is in the high loss regime. In this paper, we present numerical tests of the predicted universal impedance variance ratios and show that an insufficient sample size can lead to apparent deviation from the theory, particularly in the low loss regime. Experimental tests are carried out in three two-port microwave cavities with varied loss parameters, including a novel experimental system with a superconducting microwave billiard, to test the variance-ratio predictions in the low loss time-reversal-invariant regime. It is found that the experimental results agree with the theoretical predictions to the extent permitted by the finite sample size.

preprint2013arXiv

Modeling the dynamics of bivalent histone modifications

Epigenetic modifications to histones may promote either activation or repression of the transcription of nearby genes. Recent experimental studies show that the promoters of many lineage-control genes in stem cells have "bivalent domains" in which the nucleosomes contain both active (H3K4me3) and repressive (H3K27me3) marks. It is generally agreed that bivalent domains play an important role in stem cell differentiation, but the underlying mechanisms remain unclear. Here we formulate a mathematical model to investigate the dynamic properties of histone modification patterns. We then illustrate that our modeling framework can be used to capture key features of experimentally observed combinatorial chromatin states.

preprint2013arXiv

Nonlinear Time-Reversal of Classical Waves: Experiment and Model

We consider time-reversal of electromagnetic waves in a closed, wave-chaotic system containing a discrete, passive, harmonic-generating nonlinearity. An experimental system is constructed as a time-reversal mirror, in which excitations generated by the nonlinearity are gathered, time-reversed, transmitted, and directed exclusively to the location of the nonlinearity. Here we show that such nonlinear objects can be purely passive (as opposed to the active nonlinearities used in previous work), and develop a high-rate secure communication system based on nonlinear time-reversal. A model of the experimental system is developed, using a star-graph network of transmission lines, with one of the lines terminated by a model diode. The model simulates time-reversal of linear and nonlinear signals, demonstrates features seen in the experimental system, and supports our interpretation of the experimental results.

preprint2013arXiv

Phase and Amplitude Dynamics in Large Systems of Coupled Oscillators: Growth Heterogeneity, Nonlinear Frequency Shifts and Cluster States

This paper addresses the behavior of large systems of heterogeneous, globally coupled oscillators each of which is described by the generic Landau-Stuart equation, which incorporates both phase and amplitude dynamics of individual oscillators. One goal of our paper is to investigate the effect of a spread in the amplitude growth parameter of the oscillators and of the effect of a homogeneous nonlinear frequency shift. Both of these effects are of potential relevance to recently reported experiments. Our second goal is to gain further understanding of the macroscopic system dynamics at large coupling strength, and its dependence on the nonlinear frequency shift parameter. It is proven that at large coupling strength, if the nonlinear frequency shift parameter is below a certain value, then there is a unique attractor for which the oscillators all clump at a single amplitude and uniformly rotating phase (we call this a single-cluster "locked state"). Using a combination of analytical and numerical methods, we show that at higher values of the nonlinear frequency shift parameter, the single-cluster locked state attractor continues to exist, but other types of coexisting attractors emerge. These include two-cluster locked states, periodic orbits, chaotic orbits, and quasiperiodic orbits.

preprint2013arXiv

Predictability and suppression of extreme events in complex systems

In many complex systems, large events are believed to follow power-law, scale-free probability distributions, so that the extreme, catastrophic events are unpredictable. Here, we study coupled chaotic oscillators that display extreme events. The mechanism responsible for the rare, largest events makes them distinct and their distribution deviates from a power-law. Based on this mechanism identification, we show that it is possible to forecast in real time an impending extreme event. Once forecasted, we also show that extreme events can be suppressed by applying tiny perturbations to the system.

preprint2013arXiv

Predicting the statistics of wave transport through chaotic cavities by the Random Coupling Model: a review and recent progress

In this review, a model (the Random Coupling Model) that gives a statistical description of the coupling of radiation into and out of large enclosures through localized and/or distributed channels is presented. The Random Coupling Model combines both deterministic and statistical phenomena. The model makes use of wave chaos theory to extend the classical modal description of the cavity fields in the presence of boundaries that lead to chaotic ray trajectories. The model is based on a clear separation between the universal statistical behavior of the isolated chaotic system, and the deterministic coupling channel characteristics. Moreover, the ability of the random coupling model to describe interconnected cavities, aperture coupling, and the effects of short ray trajectories is discussed. A relation between the random coupling model and other formulations adopted in acoustics, optics, and statistical electromagnetics, is examined. In particular, a rigorous analogy of the random coupling model with the Statistical Energy Analysis used in acoustics is presented.

preprint2013arXiv

Quantifying Volume Changing Perturbations in a Wave Chaotic System

A sensor was developed to quantitatively measure perturbations which change the volume of a wave chaotic cavity while leaving its shape intact. The sensors work in the time domain by using either scattering fidelity of the transmitted signals or time reversal mirrors. The sensors were tested experimentally by inducing volume changing perturbations to a one cubic meter mixed chaotic and regular billiard system. Perturbations which caused a volume change that is as small as 54 parts in a million were quantitatively measured. These results were obtained by using electromagnetic waves with a wavelength of about 5cm, therefore, the sensor is sensitive to extreme sub-wavelength changes of the boundaries of a cavity. The experimental results were compared with Finite Difference Time Domain (FDTD) simulation results, and good agreement was found. Furthermore, the sensor was tested using a frequency domain approach on a numerical model of the star graph, which is a representative wave chaotic system. These results open up interesting applications such as: monitoring the spatial uniformity of the temperature of a homogeneous cavity during heating up / cooling down procedures, verifying the uniform displacement of a fluid inside a wave chaotic cavity by another fluid, etc.

preprint2013arXiv

Spatially embedded growing small-world networks

Networks in nature are often formed within a spatial domain in a dynamical manner, gaining links and nodes as they develop over time. We propose a class of spatially-based growing network models and investigate the relationship between the resulting statistical network properties and the dimension and topology of the space in which the networks are embedded. In particular, we consider models in which nodes are placed one by one in random locations in space, with each such placement followed by configuration relaxation toward uniform node density, and connection of the new node with spatially nearby nodes. We find that such growth processes naturally result in networks with small-world features, including a short characteristic path length and nonzero clustering. These properties do not appear to depend strongly on the topology of the embedding space, but do depend strongly on its dimension; higher-dimensional spaces result in shorter path lengths but less clustering.

preprint2013arXiv

Stability of Boolean networks: The joint effects of topology and update rules

We study the stability of orbits in large Boolean networks with given complex topology. We impose no restrictions on the form of the update rules, which may be correlated with local topological properties of the network. While recent past work has addressed the separate effects of nontrivial network topology and certain special classes of update rules on stability, only crude results exist about how these effects interact. We present a widely applicable solution to this problem. Numerical experiments confirm our theory and show that local correlations between topology and update rules can have profound effects on the qualitative behavior of these systems.

preprint2013arXiv

Weakly Explosive Percolation in Directed Networks

Percolation, the formation of a macroscopic connected component, is a key feature in the description of complex networks. The dynamical properties of a variety of systems can be understood in terms of percolation, including the robustness of power grids and information networks, the spreading of epidemics and forest fires, and the stability of gene regulatory networks. Recent studies have shown that if network edges are added "competitively" in undirected networks, the onset of percolation is abrupt or "explosive." The unusual qualitative features of this phase transition have been the subject of much recent attention. Here we generalize this previously studied network growth process from undirected networks to directed networks and use finite-size scaling theory to find several scaling exponents. We find that this process is also characterized by a very rapid growth in the giant component, but that this growth is not as sudden as in undirected networks.

preprint2012arXiv

Continuum modeling of the equilibrium and stability of animal flocks

Groups of animals often tend to arrange themselves in flocks that have characteristic spatial attributes and temporal dynamics. Using a dynamic continuum model for a flock of individuals, we find equilibria of finite spatial extent where the density goes continuously to zero at a well-defined flock edge, and we discuss conditions on the model that allow for such solutions. We also demonstrate conditions under which, as the flock size increases, the interior density in our equilibria tends to an approximately uniform value. Motivated by observations of starling flocks that are relatively thin in a direction transverse to the direction of flight, we investigate the stability of infinite, planar-sheet flock equilibria. We find that long- wavelength perturbations along the sheet are unstable for the class of models that we investigate. This has the conjectured consequence that sheet-like flocks of arbitrarily large transverse extent relative to their thickness do not occur. However, we also show that our model admits approximately sheet-like, 'pancake-shaped', three-dimensional ellipsoidal equilibria with definite aspect ratios (transverse length- scale to flock thickness) determined by anisotropic perceptual/response characteristics of the flocking individuals, and we argue that these pancake-like equilibria are stable to the previously mentioned sheet instability.

preprint2012arXiv

Detecting Functional Communities in Complex Networks

We consider an alternate definition of community structure that is functionally motivated. We define network community structure-based on the function the network system is intended to perform. In particular, as a specific example of this approach, we consider communities whose function is enhanced by the ability to synchronize and/or by resilience to node failures. Previous work has shown that, in many cases, the largest eigenvalue of the network's adjacency matrix controls the onset of both synchronization and percolation processes. Thus, for networks whose functional performance is dependent on these processes, we propose a method that divides a given network into communities based on maximizing a function of the largest eigenvalues of the adjacency matrices of the resulting communities. We also explore the differences between the partitions obtained by our method and the modularity approach (which is based solely on consideration of network structure). We do this for several different classes of networks. We find that, in many cases, modularity-based partitions do almost as well as our function-based method in finding functional communities, even though modularity does not specifically incorporate consideration of function.

preprint2012arXiv

Dynamical Instability in Boolean Networks as a Percolation Problem

Boolean networks, widely used to model gene regulation, exhibit a phase transition between regimes in which small perturbations either die out or grow exponentially. We show and numerically verify that this phase transition in the dynamics can be mapped onto a static percolation problem which predicts the long-time average Hamming distance between perturbed and unperturbed orbits.

preprint2012arXiv

Echoes and revival echoes in systems of anharmonically confined atoms

We study echoes and what we call 'revival echoes' for a collection of atoms that are described by a single quantum wavefunction and are confined in a weakly anharmonic trap. The echoes and revival echoes are induced by applying two, successive temporally localized potential perturbations to the confining potential, one at time $t=0$, and a smaller one at time $t=τ$. Pulse-like responses in the expectation value of position $<x(t)>$ are predicted at $t \approx nτ$ ($n=2,3,...$) and are particularly evident at $t \approx 2τ$. A novel result of our study is the finding of 'revival echoes'. Revivals (but not echoes) occur even if the second perturbation is absent. In particular, in the absence of the second perturbation, the response to the first perturbation dies away, but then reassembles, producing a response at revival times $mT_x$ ($m=1,2,...$). Including the second perturbation at $t=τ$, we find temporally localized responses, revival echoes, both before and after $t\approx mT_x$, e.g., at $t\approx m T_x-n τ$ (pre-revival echoes) and at $t\approx mT_x+nτ$, (post-revival echoes), where $m$ and $n$ are $1,2,...$ . Depending on the form of the perturbations, the 'principal' revival echoes at $t \approx T_x \pm τ$ can be much larger than the echo at $t \approx 2τ$. We develop a perturbative model for these phenomena, and compare its predictions to the numerical solutions of the time-dependent Schrödinger Equation. The scaling of the size of the various echoes and revival echoes as a function of the symmetry and size of the perturbations applied at $t=0$ and $t=τ$ is investigated. We also study the presence of revivals and revival echoes in higher moments of position, $<x^p (t)>$, $p>1$, and the effect of atom-atom interactions on these phenomena.

preprint2012arXiv

Statistical Properties of Avalanches in Networks

We characterize the distributions of size and duration of avalanches propagating in complex networks. By an avalanche we mean the sequence of events initiated by the externally stimulated `excitation' of a network node, which may, with some probability, then stimulate subsequent firings of the nodes to which it is connected, resulting in a cascade of firings. This type of process is relevant to a wide variety of situations, including neuroscience, cascading failures on electrical power grids, and epidemology. We find that the statistics of avalanches can be characterized in terms of the largest eigenvalue and corresponding eigenvector of an appropriate adjacency matrix which encodes the structure of the network. By using mean-field analyses, previous studies of avalanches in networks have not considered the effect of network structure on the distribution of size and duration of avalanches. Our results apply to individual networks (rather than network ensembles) and provide expressions for the distributions of size and duration of avalanches starting at particular nodes in the network. These findings might find application in the analysis of branching processes in networks, such as cascading power grid failures and critical brain dynamics. In particular, our results show that some experimental signatures of critical brain dynamics (i.e., power-law distributions of size and duration of neuronal avalanches), are robust to complex underlying network topologies.

preprint2011arXiv

A first-principles model of time-dependent variations in transmission through a fluctuating scattering environment

Fading is the time-dependent variation in transmitted signal strength through a complex medium, due to interference or temporally evolving multipath scattering. In this paper we use random matrix theory (RMT) to establish a first-principles model for fading, including both universal and non-universal effects. This model provides a more general understanding of the most common statistical models (Rayleigh fading and Rice fading) and provides a detailed physical basis for their parameters. We also report experimental tests on two ray-chaotic microwave cavities. The results show that our RMT model agrees with the Rayleigh/Rice models in the high loss regime, but there are strong deviations in low-loss systems where the RMT approach describes the data well.

preprint2011arXiv

Approximating the Largest Eigenvalue of the Modified Adjacency Matrix of Networks with Heterogeneous Node Biases

Motivated by its relevance to various types of dynamical behavior of network systems, the maximum eigenvalue $λ_Q$ of the adjacency matrix $A$ of a network has been considered, and mean-field-type approximations to $λ_Q$ have been developed for different kinds of networks. Here $A$ is defined by $A_{ij} = 1$ ($A_{ij} = 0$) if there is (is not) a directed network link to $i$ from $j$. However, in at least two recent problems involving networks with heterogeneous node properties (percolation on a directed network and the stability of Boolean models of gene networks), an analogous but different eigenvalue problem arises, namely, that of finding the largest eigenvalue $λ_Q$ of the matrix $Q$, where $Q_{ij} = q_i A_{ij}$ and the `bias' $q_i$ may be different at each node $i$. (In the previously mentioned percolation and gene network contexts, $q_i$ is a probability and so lies in the range $0 \le q_i \le 1$.) The purposes of this paper are to extend the previous considerations of the maximum eigenvalue $λ_A$ of $A$ to $λ_Q$, to develop suitable analytic approximations to $λ_Q$, and to test these approximations with numerical experiments. In particular, three issues considered are (i) the effect of the correlation (or anticorrelation) between the value of $q_i$ and the number of links to and from node $i$; (ii) the effect of correlation between the properties of two nodes at either end of a network link (`assortativity'); and (iii) the effect of community structure allowing for a situation in which different $q$-values are associated with different communities.

preprint2011arXiv

Cluster Synchrony in Systems of Coupled Phase Oscillators with Higher-Order Coupling

We study the phenomenon of cluster synchrony that occurs in ensembles of coupled phase oscillators when higher-order modes dominate the coupling between oscillators. For the first time, we develop a complete analytic description of the dynamics in the limit of a large number of oscillators and use it to quantify the degree of cluster synchrony, cluster asymmetry, and switching. We use a variation of the recent dimensionality-reduction technique of Ott and Antonsen [Chaos {\bf 18}, 037113 (2008)] and find an analytic description of the degree of cluster synchrony valid on a globally attracting manifold. Shaped by this manifold, there is an infinite family of steady-state distributions of oscillators, resulting in a high degree of multi-stability in the cluster asymmetry. We also show how through external forcing the degree of asymmetry can be controlled, and suggest that systems displaying cluster synchrony can be used to encode and store data.

preprint2011arXiv

Dynamics and Pattern Formation in Large Systems of Spatially-Coupled Oscillators with Finite Response Times

We consider systems of many spatially distributed phase oscillators that interact with their neighbors. Each oscillator is allowed to have a different natural frequency, as well as a different response time to the signals it receives from other oscillators in its neighborhood. Using the ansatz of Ott and Antonsen (Ref. \cite{OA1}) and adopting a strategy similar to that employed in the recent work of Laing (Ref. \cite{Laing2}), we reduce the microscopic dynamics of these systems to a macroscopic partial-differential-equation description. Using this macroscopic formulation, we numerically find that finite oscillator response time leads to interesting spatio-temporal dynamical behaviors including propagating fronts, spots, target patterns, chimerae, spiral waves, etc., and we study interactions and evolutionary behaviors of these spatio-temporal patterns.

preprint2011arXiv

Effects of network topology, transmission delays, and refractoriness on the response of coupled excitable systems to a stochastic stimulus

We study the effects of network topology on the response of networks of coupled discrete excitable systems to an external stochastic stimulus. We extend recent results that characterize the response in terms of spectral properties of the adjacency matrix by allowing distributions in the transmission delays and in the number of refractory states, and by developing a nonperturbative approximation to the steady state network response. We confirm our theoretical results with numerical simulations. We find that the steady state response amplitude is inversely proportional to the duration of refractoriness, which reduces the maximum attainable dynamic range. We also find that transmission delays alter the time required to reach steady state. Importantly, neither delays nor refractoriness impact the general prediction that criticality and maximum dynamic range occur when the largest eigenvalue of the adjacency matrix is unity.

preprint2011arXiv

Ensemble regional data assimilation using joint states

We propose a data assimilation scheme that produces the analyses for a global and an embedded limited area model simultaneously, considering forecast information from both models. The purpose of the proposed approach is twofold. First, we expect that the global analysis will benefit from incorporation of information from the higher resolution limited area model. Second, our method is expected to produce a limited area analysis that is more strongly constrained by the large scale flow than a conventional limited area analysis. The proposed scheme minimizes a cost function in which the control variable is the joint state of the global and the limited area models. In addition, the cost function includes a constraint term that penalizes large differences between the global and the limited area state estimates. The proposed approach is tested by idealized experiments, using `toy' models introduced by Lorenz in 2005. The results of these experiments suggest that the proposed approach improves the global analysis within and near the limited area domain and the regional analysis near the lateral boundaries. These analysis improvements lead to forecast improvements in both the global and the limited area models.

preprint2011arXiv

Multiscale Dynamics in Communities of Phase Oscillators

We investigate the dynamics of systems of many coupled phase oscillators with het- erogeneous frequencies. We suppose that the oscillators occur in M groups. Each oscillator is connected to other oscillators in its group with "attractive" coupling, such that the coupling promotes synchronization within the group. The coupling between oscillators in different groups is "repulsive"; i.e., their oscillation phases repel. To address this problem, we reduce the governing equations to a lower-dimensional form via the ansatz of Ott and Antonsen . We first consider the symmetric case where all group parameters are the same, and the attractive and repulsive coupling are also the same for each of the M groups. We find a manifold L of neutrally stable equilibria, and we show that all other equilibria are unstable. For M \geq 3, L has dimension M - 2, and for M = 2 it has dimension 1. To address the general asymmetric case, we then introduce small deviations from symmetry in the group and coupling param- eters. Doing a slow/fast timescale analysis, we obtain slow time evolution equations for the motion of the M groups on the manifold L. We use these equations to study the dynamics of the groups and compare the results with numerical simulations.

preprint2011arXiv

On the propagation of information and the use of localization in ensemble Kalman filtering

Several localized versions of the ensemble Kalman filter have been proposed. Although tests applying such schemes have proven them to be extremely promising, a full basic understanding of the rationale and limitations of localization is currently lacking. It is one of the goals of this paper to contribute toward addressing this issue. The second goal is to elucidate the role played by chaotic wave dynamics in the propagation of information and the resulting impact on forecasts. To accomplish these goals, the principal tool used here will be analysis and interpretation of numerical experiments on a toy atmospheric model introduced by Lorenz in 2005. Propagation of the wave packets of this model is shown. It is found that, when an ensemble Kalman filter scheme is employed, the spatial correlation function obtained at each forecast cycle by averaging over the background ensemble members is short ranged, and this is in strong contrast to the much longer range correlation function obtained by averaging over states from free evolution of the model. Propagation of the effects of observations made in one region on forecasts in other regions is studied. The error covariance matrices from the analyses with localization and without localization are compared. From this study, major characteristics of the localization process and information propagation are extracted and summarized.

preprint2010arXiv

Comment on "Long Time Evolution of Phase Oscillator Systems" [Chaos 19,023117 (2009), arXiv:0902.2773]

A previous paper (arXiv:0902.2773, henceforth referred to as I) considered a general class of problems involving the evolution of large systems of globally coupled phase oscillators. It was shown there that, in an appropriate sense, the solutions to these problems are time asymptotically attracted toward a reduced manifold of system states (denoted M). This result has considerable utility in the analysis of these systems, as has been amply demonstrated in recent papers. In this note, we show that the analysis of I can be modified in a simple way that establishes significant extensions of the range of validity of our previous result. In particular, we generalize I in the following ways: (1) attraction to M is now shown for a very general class of oscillator frequency distribution functions g(ω), and (2) a previous restriction on the allowed class of initial conditions is now substantially relaxed.

preprint2010arXiv

Dynamic synchronization of a time-evolving optical network of chaotic oscillators

We present and experimentally demonstrate a technique for achieving and maintaining a global state of identical synchrony of an arbitrary network of chaotic oscillators even when the coupling strengths are unknown and time-varying. At each node an adaptive synchronization algorithm dynamically estimates the current strength of the net coupling signal to that node. We experimentally demonstrate this scheme in a network of three bidirectionally coupled chaotic optoelectronic feedback loops and we present numerical simulations showing its application in larger networks. The stability of the synchronous state for arbitrary coupling topologies is analyzed via a master stability function approach.

preprint2010arXiv

Experimental Examination of the Effect of Short Ray Trajectories in Two-port Wave-Chaotic Scattering Systems

Predicting the statistics of realistic wave-chaotic scattering systems requires, in addition to random matrix theory, introduction of system-specific information. This paper investigates experimentally one aspect of system-specific behavior, namely the effects of short ray trajectories in wave-chaotic systems open to outside scattering channels. In particular, we consider ray trajectories of limited length that enter a scattering region through a channel (port) and subsequently exit through a channel (port). We show that a suitably averaged value of the impedance can be computed from these trajectories and that this can improve the ability to describe the statistical properties of the scattering systems. We illustrate and test these points through experiments on a realistic two-port microwave scattering billiard.

preprint2010arXiv

Quantum chaos of a mixed, open system of kicked cold atoms

The quantum and classical dynamics of particles kicked by a gaussian attractive potential are studied. Classically, it is an open mixed system (the motion in some parts of the phase space is chaotic, and in some parts it is regular). The fidelity (Lochshmidt echo) is found to exhibit oscillations that can be determined from classical considerations but are sensitive to phase space structures that are smaller than Planck's constant. Families of quasi-energies are determined from classical phase space structures. Substantial differences between the classical and quantum dynamics are found for time dependent scattering. It is argued that the system can be experimentally realized by cold atoms kicked by a gaussian light beam.

preprint2010arXiv

Sensing Small Changes in a Wave Chaotic Scattering System

Classical analogs of the quantum mechanical concepts of the Loschmidt Echo and quantum fidelity are developed with the goal of detecting small perturbations in a closed wave chaotic region. Sensing techniques that employ a one-recording-channel time-reversal-mirror, which in turn relies on time reversal invariance and spatial reciprocity of the classical wave equation, are introduced. In analogy with quantum fidelity, we employ Scattering Fidelity techniques which work by comparing response signals of the scattering region, by means of cross correlation and mutual information of signals. The performance of the sensing techniques is compared for various perturbations induced experimentally in an acoustic resonant cavity. The acoustic signals are parametrically processed to mitigate the effect of dissipation and to vary the spatial diversity of the sensing schemes. In addition to static boundary condition perturbations at specified locations, perturbations to the medium of wave propagation are shown to be detectable, opening up various real world sensing applications in which a false negative cannot be tolerated.

preprint2010arXiv

Spontaneous synchronization of coupled oscillator systems with frequency adaptation

We study the synchronization of Kuramoto oscillators with all-to-all coupling in the presence of slow, noisy frequency adaptation. In this paper we develop a new model for oscillators which adapt both their phases and frequencies. It is found that this model naturally reproduces some observed phenomena that are not qualitatively produced by the standard Kuramoto model, such as long waiting times before the synchronization of clapping audiences. By assuming a self-consistent steady state solution, we find three stability regimes for the coupling constant k, separated by critical points k1 and k2: (i) for k<k1, only the stable incoherent state exists; (ii) for k>k2, the incoherent state becomes unstable and only the synchronized state exists; (iii) for k1<k<k2, both the incoherent and synchronized states are stable. In the bistable regime spontaneous transitions between the incoherent and synchronized states are observed for finite ensembles. These transitions are well described as a stochastic process on the order parameter r undergoing fluctuations due to the system's finite size, leading to the following conclusions: (a) in the bistable regime, the average waiting time of an incoherent-to-coherent transition can be predicted by using Kramer's escape time formula and grows exponentially with the number of oscillators; (b) when the incoherent state is unstable (k>k2), the average waiting time grows logarithmically with the number of oscillators.

preprint2009arXiv

Adaptive synchronization of coupled chaotic oscillators

We experimentally demonstrate and numerically simulate a new adaptive method to maintain synchronization between coupled nonlinear chaotic oscillators, when the coupling between the systems is unknown and time-varying (e.g., due to environmental parameter drift). The technique is applied to optoelectronic feedback loops exhibiting high dimensional chaotic dynamics. In addition to keeping the two systems isochronally synchronized in the presence of a priori unknown time-varying coupling strength, the technique provides an estimate of the time-varying coupling.

preprint2009arXiv

Complex Dynamics and Synchronization of Delayed-Feedback Nonlinear Oscillators

We describe a flexible and modular delayed-feedback nonlinear oscillator that is capable of generating a wide range of dynamical behaviours, from periodic oscillations to high-dimensional chaos. The oscillator uses electrooptic modulation and fibre-optic transmission, with feedback and filtering implemented through real-time digital-signal processing. We consider two such oscillators that are coupled to one another, and we identify the conditions under which they will synchronize. By examining the rates of divergence or convergence between two coupled oscillators, we quantify the maximum Lyapunov exponents or transverse Lyapunov exponents of the system, and we present an experimental method to determine these rates that does not require a mathematical model of the system. Finally, we demonstrate a new adaptive control method that keeps two oscillators synchronized even when the coupling between them is changing unpredictably.

preprint2009arXiv

Long Time Evolution of Phase Oscillator Systems

It is shown, under weak conditions, that the dynamical evolution of an important class of large systems of globally coupled, heterogeneous frequency, phase oscillators is, in an appropriate physical sense, time-asymptotically attracted toward a reduced manifold of system states. This manifold, which is invariant under the system evolution, was previously known and used to facilitate the discovery of attractors and bifurcations of such systems. The result of this paper establishes that attractors for the order parameter dynamics obtained by restriction to this reduced manifold are, in fact, the only such attractors of the full system. Thus all long time dynamical behavior of the order parameters of these systems can be obtained by restriction to the reduced manifold.

preprint2009arXiv

Obstacle and predator avoidance by a flock

The modeling and investigation of the dynamics and configurations of animal groups is a subject of growing attention. In this paper, we present a continuum model of flocking and use it to investigate the reaction of a flock to an obstacle or an attacking predator. We show that the flock response is in the form of density disturbances that resemble Mach cones whose configuration is determined by the anisotropic propagation of waves through the flock. We analytically and numerically test relations that predict the Mach wedge angles, disturbance heights, and wake widths. We find that these expressions are insensitive to many of the parameters of the model.

preprint2009arXiv

The effect of network topology on the stability of discrete state models of genetic control

Boolean networks have been proposed as potentially useful models for genetic control. An important aspect of these networks is the stability of their dynamics in response to small perturbations. Previous approaches to stability have assumed uncorrelated random network structure. Real gene networks typically have nontrivial topology significantly different from the random network paradigm. In order to address such situations, we present a general method for determining the stability of large Boolean networks of any specified network topology and predicting their steady-state behavior in response to small perturbations. Additionally, we generalize to the case where individual genes have a distribution of `expression biases,' and we consider non-synchronous update, as well as extension of our method to non-Boolean models in which there are more than two possible gene states. We find that stability is governed by the maximum eigenvalue of a modified adjacency matrix, and we test this result by comparison with numerical simulations. We also discuss the possible application of our work to experimentally inferred gene networks.

preprint2009arXiv

The effect of short ray trajectories on the scattering statistics of wave chaotic systems

In many situations, the statistical properties of wave systems with chaotic classical limits are well-described by random matrix theory. However, applications of random matrix theory to scattering problems require introduction of system specific information into the statistical model, such as the introduction of the average scattering matrix in the Poisson kernel. Here it is shown that the average impedance matrix, which also characterizes the system-specific properties, can be expressed in terms of classical trajectories that travel between ports and thus can be calculated semiclassically. Theoretical results are compared with numerical solutions for a model wave-chaotic system.

preprint2009arXiv

The stability of adaptive synchronization of chaotic systems

In past works, various schemes for adaptive synchronization of chaotic systems have been proposed. The stability of such schemes is central to their utilization. As an example addressing this issue, we consider a recently proposed adaptive scheme for maintaining the synchronized state of identical coupled chaotic systems in the presence of a priori unknown slow temporal drift in the couplings. For this illustrative example, we develop an extension of the master stability function technique to study synchronization stability with adaptive coupling. Using this formulation, we examine local stability of synchronization for typical chaotic orbits and for unstable periodic orbits within the synchronized chaotic attractor (bubbling). Numerical experiments illustrating the results are presented. We observe that the stable range of synchronism can be sensitively dependent on the adaption parameters, and we discuss the strong implication of bubbling for practically achievable adaptive synchronization.

preprint2009arXiv

Universal and non-universal properties of wave chaotic scattering systems

The application of random matrix theory to scattering requires introduction of system-specific information. This paper shows that the average impedance matrix, which characterizes such system-specific properties, can be semiclassically calculated in terms of ray trajectories between ports. Theoretical predictions are compared with experimental results for a microwave billiard, demonstrating that the theory successfully uncovered universal statistics of wave-chaotic scattering systems.

preprint2005arXiv

Formation of Multifractal Population Patterns from Reproductive Growth and Local Resettlement

We consider the general character of the spatial distribution of a population that grows through reproduction and subsequent local resettlement of new population members. We present several simple one and two-dimensional point placement models to illustrate possible generic behavior of these distributions. We show, numerically and analytically, that these models all lead to multifractal spatial distributions of population. Additionally, we make qualitative links between our models and the example of the Earth at Night image, showing the Earth's nighttime man-made light as seen from space. The Earth at Night data suffer from saturation of the sensing photodetectors at high brightness (`clipping'), and we account for how this influences the determined dimension spectrum of the light intensity distribution.