Catalog footprint

What is connected

54works
36topics
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

54 published item(s)

preprint2022arXiv

Discrete-to-Continuous Extensions: Lovász extension and Morse theory

This is the first of a series of papers that develop a systematic bridge between constructions in discrete mathematics and the corresponding continuous analogs. In this paper, we establish an equivalence between Forman's discrete Morse theory on a simplicial complex and the continuous Morse theory (in the sense of any known non-smooth Morse theory) on the associated order complex via the Lovász extension. Furthermore, we propose a new version of the Lusternik-Schnirelman category on abstract simplicial complexes to bridge the classical Lusternik-Schnirelman theorem and its discrete analog on finite complexes. More generally, we can suggest a discrete Morse theory on hypergraphs by employing piecewise-linear (PL) Morse theory and Lovász extension, hoping to provide new tools for exploring the structure of hypergraphs.

preprint2022arXiv

Geometry of Data

Topological data analysis asks when balls in a metric space $(X,d)$ intersect. Geometric data analysis asks how much balls have to be enlarged to intersect. We connect this principle to the traditional core geometric concept of curvature. This enables us, on one hand, to reconceptualize curvature and link it to the geometric notion of hyperconvexity. On the other hand, we can then also understand methods of topological data analysis from a geometric perspective.

preprint2022arXiv

Local Detour Centrality: A Novel Local Centrality Measure for Weighted Networks

Centrality, in some sense, captures the extent to which a vertex controls the flow of information in a network. Here, we propose Local Detour Centrality as a novel centrality-based betweenness measure that captures the extent to which a vertex shortens paths between neighboring vertices as compared to alternative paths. After presenting our measure, we demonstrate empirically that it differs from other leading central measures, such as betweenness, degree, closeness, and the number of triangles. Through an empirical case study, we provide a possible interpretation for Local Detour Centrality as a measure that captures the extent to which a word is characterized by contextual diversity within a semantic network. We then examine the relationship between our measure and the accessibility to knowledge stored in memory. To do so, we show that words that occur in several different and distinct contexts are significantly more effective in facilitating the retrieval of subsequent words than are words that lack this contextual diversity.

preprint2022arXiv

Scale Free Avalanches in Excitatory-Inhibitory Populations of Spiking Neurons with Conductance Based Synaptic Currents

We investigate spontaneous critical dynamics of excitatory and inhibitory (EI) sparsely connected populations of spiking leaky integrate-and-fire neurons with conductance-based synapses. We use a bottom-up approach to derive a single neuron gain function and a linear Poisson neuron approximation which we use to study mean-field dynamics of the EI population and its bifurcations. In the low firing rate regime, the quiescent state loses stability due to saddle-node or Hopf bifurcations. In particular, at the Bogdanov-Takens (BT) bifurcation point which is the intersection of the Hopf bifurcation and the saddle-node bifurcation lines of the 2D dynamical system, the network shows avalanche dynamics with power-law avalanche size and duration distributions. This matches the characteristics of low firing spontaneous activity in the cortex. By linearizing gain functions and excitatory and inhibitory nullclines, we can approximate the location of the BT bifurcation point. This point in the control parameter phase space corresponds to the internal balance of excitation and inhibition and a slight excess of external excitatory input to the excitatory population. Due to the tight balance of average excitation and inhibition currents, the firing of the individual cells is fluctuation-driven. Around the BT point, the spiking of neurons is a Poisson process and the population average membrane potential of neurons is approximately at the middle of the operating interval $[V_{Rest}, V_{th}]$. Moreover, the EI network is close to both oscillatory and active-inactive phase transition regimes.

preprint2022arXiv

Self Organized Criticality in a Mesoscopic Model of Excitatory-Inhibitory Neuronal Populations by Short-term and Long-term Synaptic Plasticity

In [1], we have shown that the dynamics of an interconnected population of excitatory and inhibitory spiking neurons wandering around a Bogdanov-Takens (BT)bifurcation point can generate the observed scale-free avalanches at the population level and the highly variable spike patterns of individual neurons. These characteristics match experimental findings for spontaneous intrinsic activity in the brain. In this paper, we address the mechanisms causing the system to get and remain near this BT point. We propose an effective stochastic neural field model which captures the dynamics of the mean-field model. We show how the network tunes itself through local long-term synaptic plasticity by STDP and short-term synaptic depression to be close to this bifurcation point. The mesoscopic model that we derive matches the directed percolation model at the absorbing state phase transition.

preprint2021arXiv

Geometric Sampling of Networks

Motivated by the methods and results of manifold sampling based on Ricci curvature, we propose a similar approach for networks. To this end we make appeal to three types of discrete curvature, namely the graph Forman-, full Forman- and Haantjes-Ricci curvatures for edge-based and node-based sampling. We present the results of experiments on real life networks, as well as for square grids arising in Image Processing. Moreover, we consider fitting Ricci flows and we employ them for the detection of networks' backbone. We also develop embedding kernels related to the Forman-Ricci curvatures and employ them for the detection of the coarse structure of networks, as well as for network visualization with applications to SVM. The relation between the Ricci curvature of the original manifold and that of a Ricci curvature driven discretization is also studied.

preprint2021arXiv

It Means More if It Sounds Good: Yet Another Hypothesis Concerning the Evolution of Polysemous Words

This position paper looks into the formation of language and shows ties between structural properties of the words in the English language and their polysemy. Using Ollivier-Ricci curvature over a large graph of synonyms to estimate polysemy it shows empirically that the words that arguably are easier to pronounce also tend to have multiple meanings.

preprint2021arXiv

Network geometry and market instability

The complexity of financial markets arise from the strategic interactions among agents trading stocks, which manifest in the form of vibrant correlation patterns among stock prices. Over the past few decades, complex financial markets have often been represented as networks whose interacting pairs of nodes are stocks, connected by edges that signify the correlation strengths. However, we often have interactions that occur in groups of three or more nodes, and these cannot be described simply by pairwise interactions but we also need to take the relations between these interactions into account. Only recently, researchers have started devoting attention to the higher-order architecture of complex financial systems, that can significantly enhance our ability to estimate systemic risk as well as measure the robustness of financial systems in terms of market efficiency. Geometry-inspired network measures, such as the Ollivier-Ricci curvature and Forman-Ricci curvature, can be used to capture the network fragility and continuously monitor financial dynamics. Here, we explore the utility of such discrete Ricci curvatures in characterizing the structure of financial systems, and further, evaluate them as generic indicators of the market instability. For this purpose, we examine the daily returns from a set of stocks comprising the USA S&P-500 and the Japanese Nikkei-225 over a 32-year period, and monitor the changes in the edge-centric network curvatures. We find that the different geometric measures capture well the system-level features of the market and hence we can distinguish between the normal or `business-as-usual' periods and all the major market crashes. This can be very useful in strategic designing of financial systems and regulating the markets in order to tackle financial instabilities.

preprint2020arXiv

A Simple Differential Geometry for Complex Networks

We introduce new definitions of sectional, Ricci and scalar curvature for networks and their higher dimensional counterparts, derived from two classical notions of curvature for curves in general metric spaces, namely, the Menger curvature and the Haantjes curvature. These curvatures are applicable to unweighted or weighted and undirected or directed networks, and are more intuitive and easier to compute than other network curvatures. In particular, the proposed curvatures based on the interpretation of Haantjes definition as geodesic curvature allow us to give a network analogue of the classical local Gauss-Bonnet theorem. Furthermore, we propose even simpler and more intuitive proxies for the Haantjes curvature that allow for even faster and easier computations in large-scale networks. In addition, we also investigate the embedding properties of the proposed Ricci curvatures. Lastly, we also investigate the behaviour, both on model and real-world networks, of the curvatures introduced herein with more established notions of Ricci curvature and other widely-used network measures.

preprint2020arXiv

Coupled Dynamics on Hypergraphs: Master Stability of Steady States and Synchronization

In the study of dynamical systems on networks/graphs, a key theme is how the network topology influences stability for steady states or synchronized states. Ideally, one would like to derive conditions for stability or instability that instead of microscopic details of the individual nodes/vertices rather make the influence of the network coupling topology visible. The master stability function is an important such tool to achieve this goal. Here we generalize the master stability approach to hypergraphs. A hypergraph coupling structure is important as it allows us to take into account arbitrary higher-order interactions between nodes. As for instance in the theory of coupled map lattices, we study Laplace type interaction structures in detail. Since the spectral theory of Laplacians on hypergraphs is richer than on graphs, we see the possibility of new dynamical phenomena. More generally, our arguments provide a blueprint for how to generalize dynamical structures and results from graphs to hypergraphs.

preprint2020arXiv

Deriving pairwise transfer entropy from network structure and motifs

Transfer entropy is an established method for quantifying directed statistical dependencies in neuroimaging and complex systems datasets. The pairwise (or bivariate) transfer entropy from a source to a target node in a network does not depend solely on the local source-target link weight, but on the wider network structure that the link is embedded in. This relationship is studied using a discrete-time linearly-coupled Gaussian model, which allows us to derive the transfer entropy for each link from the network topology. It is shown analytically that the dependence on the directed link weight is only a first approximation, valid for weak coupling. More generally, the transfer entropy increases with the in-degree of the source and decreases with the in-degree of the target, indicating an asymmetry of information transfer between hubs and low-degree nodes. In addition, the transfer entropy is directly proportional to weighted motif counts involving common parents or multiple walks from the source to the target, which are more abundant in networks with a high clustering coefficient than in random networks. Our findings also apply to Granger causality, which is equivalent to transfer entropy for Gaussian variables. Moreover, similar empirical results on random Boolean networks suggest that the dependence of the transfer entropy on the in-degree extends to nonlinear dynamics.

preprint2020arXiv

Differentiation of measures on complete Riemannian manifolds

In this note we give a new proof of a version of the Besicovitch covering theorem, given in \cite{EG1992}, \cite{Bogachev2007} and extended in \cite{Federer1969}, for locally finite Borel measures on finite dimensional complete Riemannian manifolds $(M,g)$. As a consequence, we prove a differentiation theorem for Borel measures on $(M,g)$, which gives a formula for the Radon-Nikodym density of two nonnegative locally finite Borel measures $ν_1, ν_2$ on $(M, g)$ such that $ν_1 \ll ν_2$, extending the known case when $(M, g)$ is a standard Euclidean space.

preprint2020arXiv

From the Jordan product to Riemannian geometries on classical and quantum states

The Jordan product on the self-adjoint part of a finite-dimensional $C^{*}$-algebra $\mathscr{A}$ is shown to give rise to Riemannian metric tensors on suitable manifolds of states on $\mathscr{A}$, and the covariant derivative, the geodesics, the Riemann tensor, and the sectional curvature of all these metric tensors are explicitly computed. In particular, it is proved that the Fisher--Rao metric tensor is recovered in the Abelian case, that the Fubini--Study metric tensor is recovered when we consider pure states on the algebra $\mathcal{B}(\mathcal{H})$ of linear operators on a finite-dimensional Hilbert space $\mathcal{H}$, and that the Bures--Helstrom metric tensors is recovered when we consider faithful states on $\mathcal{B}(\mathcal{H})$. Moreover, an alternative derivation of these Riemannian metric tensors in terms of the GNS construction associated to a state is presented. In the case of pure and faithful states on $\mathcal{B}(\mathcal{H})$, this alternative geometrical description clarifies the analogy between the Fubini--Study and the Bures--Helstrom metric tensor.

preprint2020arXiv

Spherical Bernstein theorems for codimension 1 and 2

A result of B.Solomon (On the Gauss map of an area-minimizing hypersurface. 1984. Journal of Differential Geometry, 19(1), 221-232.) says that a compact minimal hypersurface $M^k$ of the sphere $S^{k+1}$ with $H^1(M)=0$, whose Gauss map omits a neighborhood of an $S^{k-1}$ equator, is totally geodesic in $S^{k+1}$. We develop a new proof strategy which can also obtain an analogous result for codimension 2 compact minimal submanifolds of $S^{k+1}$.

preprint2019arXiv

Manifolds of classical probability distributions and quantum density operators in infinite dimensions

The manifold structure of subsets of classical probability distributions and quantum density operators in infinite dimensions is investigated in the context of $C^{*}$-algebras and actions of Banach-Lie groups. Specificaly, classical probability distributions and quantum density operators may be both described as states (in the functional analytic sense) on a given $C^{*}$-algebra $\mathscr{A}$ which is Abelian for Classical states, and non-Abelian for Quantum states. In this contribution, the space of states $\mathscr{S}$ of a possibly infinite-dimensional, unital $C^{*}$-algebra $\mathscr{A}$ is partitioned into the disjoint union of the orbits of an action of the group $\mathscr{G}$ of invertible elements of $\mathscr{A}$. Then, we prove that the orbits through density operators on an infinite-dimensional, separable Hilbert space $\mathcal{H}$ are smooth, homogeneous Banach manifolds of $\mathscr{G}=\mathcal{GL}(\mathcal{H})$, and, when $\mathscr{A}$ admits a faithful tracial state $τ$ like it happens in the Classical case when we consider probability distributions with full support, we prove that the orbit through $τ$ is a smooth, homogeneous Banach manifold for $\mathscr{G}$.

preprint2019arXiv

Topology and curvature of metric spaces

We develop a new concept of non-positive curvature for metric spaces, based on intersection patterns of closed balls. In contrast to the synthetic approaches of Alexandrov and Buesemann, our concept also applies to metric spaces that might be discrete. The natural comparison spaces that emerge from our discussion are no longer Euclidean spaces, but rather tripod spaces. These tripod spaces include the hyperconvex spaces which have trivial Cech homology. This suggests a link of our geometrical method to the topological method of persistent homology in topological data analysis. We also investigate the geometry of general tripod spaces.

preprint2016arXiv

Can one see the shape of a network?

Traditionally, network analysis is based on local properties of vertices, like their degree or clustering, and their statistical behavior across the network in question. This paper develops an approach which is different in two respects. We investigate edge-based properties, and we define global characteristics of networks directly. The latter will provide our affirmative answer to the question raised in the title. More concretely, we start with Forman's notion of the Ricci curvature of a graph, or more generally, a polyhedral complex. This will allow us to pass from a graph as representing a network to a polyhedral complex for instance by filling in triangles into connected triples of edges and to investigate the resulting effect on the curvature. This is insightful for two reasons: First, we can define a curvature flow in order to asymptotically simplify a network and reduce it to its essentials. Second, using a construction of Bloch, which yields a discrete Gauss-Bonnet theorem, we have the Euler characteristic of a network as a global characteristic. These two aspects beautifully merge in the sense that the asymptotic properties of the curvature flow are indicated by that Euler characteristic.

preprint2016arXiv

Characterizing Complex Networks with Forman-Ricci Curvature and Associated Geometric Flows

We introduce Forman-Ricci curvature and its corresponding flow as characteristics for complex networks attempting to extend the common approach of node-based network analysis by edge-based characteristics. Following a theoretical introduction and mathematical motivation, we apply the proposed network-analytic methods to static and dynamic complex networks and compare the results with established node-based characteristics. Our work suggests a number of applications for data mining, including denoising and clustering of experimental data, as well as extrapolation of network evolution.

preprint2016arXiv

Ergodicity of scalar stochastic differential equations with Hölder continuous coefficients

It is well-known that for a one dimensional stochastic differential equation driven by Brownian noise, with coefficient functions satisfying the assumptions of the Yamada-Watanabe theorem \cite{yamada1,yamada2} and the Feller test for explosions \cite{feller51,feller54}, there exists a unique stationary distribution with respect to the Markov semigroup of transition probabilities. We consider systems on a restricted domain $D$ of the phase space $\mathbb{R}$ and study the rate of convergence to the stationary distribution. Using a geometrical approach that uses the so called {\it free energy function} on the density function space, we prove that the density functions, which are solutions of the Fokker-Planck equation, converge to the stationary density function exponentially under the Kullback-Leibler {divergence}, thus also in the total variation norm. The results show that there is a relation between the Bakry-Emery curvature dimension condition and the dissipativity condition of the transformed system under the Fisher-Lamperti transformation. Several applications are discussed, including the Cox-Ingersoll-Ross model and the Ait-Sahalia model in finance and the Wright-Fisher model in population genetics.

preprint2016arXiv

Forman curvature for complex networks

We adapt Forman's discretization of Ricci curvature to the case of undirected networks, both weighted and unweighted, and investigate the measure in a variety of model and real-world networks. We find that most nodes and edges in model and real networks have a negative curvature. Furthermore, the distribution of Forman curvature of nodes and edges is narrow in random and small-world networks, while the distribution is broad in scale-free and real-world networks. In most networks, Forman curvature is found to display significant negative correlation with degree and centrality measures. However, Forman curvature is uncorrelated with clustering coefficient in most networks. Importantly, we find that both model and real networks are vulnerable to targeted deletion of nodes with highly negative Forman curvature. Our results suggest that Forman curvature can be employed to gain novel insights on the organization of complex networks.

preprint2016arXiv

Forman-Ricci flow for change detection in large dynamic data sets

We present a viable solution to the challenging question of change detection in complex networks inferred from large dynamic data sets. Building on Forman's discretization of the classical notion of Ricci curvature, we introduce a novel geometric method to characterize different types of real-world networks with an emphasis on peer-to-peer networks. Furthermore we adapt the classical Ricci flow that already proved to be a powerful tool in image processing and graphics, to the case of undirected and weighted networks. The application of the proposed method on peer-to-peer networks yields insights into topological properties and the structure of their underlying data.

preprint2015arXiv

A hierarchical extension scheme for solutions of the Wright-Fisher model

We develop a global and hierarchical scheme for the forward Kolmogorov (Fokker-Planck) equation of the diffusion approximation of the Wright-Fisher model of population genetics. That model describes the random genetic drift of several alleles at the same locus in a population. The key of our scheme is to connect the solutions before and after the loss of an allele. Whereas in an approach via stochastic processes or partial differential equations, such a loss of an allele leads to a boundary singularity, from a biological or geometric perspective, this is a natural process that can be analyzed in detail. Our method depends on evolution equations for the moments of the process and a careful analysis of the boundary flux.

preprint2014arXiv

A hierarchical extension scheme for backward solutions of the Wright-Fisher model

We develop an iterative global solution scheme for the backward Kolmogorov equation of the diffusion approximation of the Wright-Fisher model of population genetics. That model describes the random genetic drift of several alleles at the same locus in a population from a backward perspective. The key of our scheme is to connect the solutions before and after the loss of an allele. Whereas in an approach via stochastic processes or partial differential equations, such a loss of an allele leads to a boundary singularity, from a biological or geometric perspective, this is a natural process that can be analyzed in detail. A clarification of the role of the boundary resolves certain uniqueness issues and enlucidates the construction of hierarchical solutions.

preprint2014arXiv

A notion of nonpositive curvature for general metric spaces

We introduce a new definition of nonpositive curvature in metric spaces and study its relationship to the existing notions of nonpositive curvature in comparison geometry. The main feature of our definition is that it applies to all metric spaces and does not rely on geodesics. Moreover, a scaled and a relaxed version of our definition are appropriate in discrete metric spaces, and are believed to be of interest in geometric data analysis.

preprint2014arXiv

Natural statistics of binaural sounds

Binaural sound localization is usually considered a discrimination task, where interaural time (ITD) and level (ILD) disparities at pure frequency channels are utilized to identify a position of a sound source. In natural conditions binaural circuits are exposed to a stimulation by sound waves originating from multiple, often moving and overlapping sources. Therefore statistics of binaural cues depend on acoustic properties and the spatial configuration of the environment. In order to process binaural sounds efficiently, the auditory system should be adapted to naturally encountered cue distributions. Statistics of cues encountered naturally and their dependence on the physical properties of an auditory scene have not been studied before. Here, we performed binaural recordings of three auditory scenes with varying spatial properties. We have analyzed empirical cue distributions from each scene by fitting them with parametric probability density functions which allowed for an easy comparison of different scenes. Higher order statistics of binaural waveforms were analyzed by performing Independent Component Analysis (ICA) and studying properties of learned basis functions. Obtained results can be related to known neuronal mechanisms and suggest how binaural hearing can be understood in terms of adaptation to the natural signal statistics.

preprint2014arXiv

Quantifying unique information

We propose new measures of shared information, unique information and synergistic information that can be used to decompose the multi-information of a pair of random variables $(Y,Z)$ with a third random variable $X$. Our measures are motivated by an operational idea of unique information which suggests that shared information and unique information should depend only on the pair marginal distributions of $(X,Y)$ and $(X,Z)$. Although this invariance property has not been studied before, it is satisfied by other proposed measures of shared information. The invariance property does not uniquely determine our new measures, but it implies that the functions that we define are bounds to any other measures satisfying the same invariance property. We study properties of our measures and compare them to other candidate measures.

preprint2014arXiv

Reconsidering unique information: Towards a multivariate information decomposition

The information that two random variables $Y$, $Z$ contain about a third random variable $X$ can have aspects of shared information (contained in both $Y$ and $Z$), of complementary information (only available from $(Y,Z)$ together) and of unique information (contained exclusively in either $Y$ or $Z$). Here, we study measures $\widetilde{SI}$ of shared, $\widetilde{UI}$ unique and $\widetilde{CI}$ complementary information introduced by Bertschinger et al., which are motivated from a decision theoretic perspective. We find that in most cases the intuitive rule that more variables contain more information applies, with the exception that $\widetilde{SI}$ and $\widetilde{CI}$ information are not monotone in the target variable $X$. Additionally, we show that it is not possible to extend the bivariate information decomposition into $\widetilde{SI}$, $\widetilde{UI}$ and $\widetilde{CI}$ to a non-negative decomposition on the partial information lattice of Williams and Beer. Nevertheless, the quantities $\widetilde{UI}$, $\widetilde{SI}$ and $\widetilde{CI}$ have a well-defined interpretation, even in the multivariate setting.

preprint2014arXiv

The uniqueness of hierarchically extended backward solutions of the Wright-Fisher model

The diffusion approximation of the Wright-Fisher model of population genetics leads to partial differentiable equations, the so-called Kolmogorov equations, with an operator that degenerates at the boundary. Standard tools do not apply, and in fact, solutions lack regularity properties. In this paper, we develop a regularising blow-up scheme for a certain class of solutions of the backward Kolmogorov equation, the iteratively extended global solutions presented in \cite{THJ5}, and establish their uniqueness. As the model describes the random genetic drift of several alleles at the same locus from a backward perspective, the singularities result from the loss of an allele. While in an analytical approach, this causes substantial difficulties, from a biological or geometric perspective, this is a natural process that can be analyzed in detail. The presented scheme regularises the solution via a tailored successive transformation of the domain.

preprint2013arXiv

Geometric analysis aspects of infinite semiplanar graphs with nonnegative curvature

In the present paper, we apply Alexandrov geometry methods to study geometric analysis aspects of infinite semiplanar graphs with nonnegative combinatorial curvature in the sense of Higuchi. We obtain the metric classification of these graphs and construct the graphs embedded in the projective plane minus one point. Moreover, we show the volume doubling property and the Poincaré inequality on such graphs. The quadratic volume growth of these graphs implies the parabolicity. In addition, we prove the polynomial growth harmonic function theorem analogous to the case of Riemannian manifolds.

preprint2013arXiv

Information geometry and sufficient statistics

Information geometry provides a geometric approach to families of statistical models. The key geometric structures are the Fisher quadratic form and the Amari-Chentsov tensor. In statistics, the notion of sufficient statistic expresses the criterion for passing from one model to another without loss of information. This leads to the question how the geometric structures behave under such sufficient statistics. While this is well studied in the finite sample size case, in the infinite case, we encounter technical problems concerning the appropriate topologies. Here, we introduce notions of parametrized measure models and tensor fields on them that exhibit the right behavior under statistical transformations. Within this framework, we can then handle the topological issues and show that the Fisher metric and the Amari-Chentsov tensor on statistical models in the class of symmetric 2-tensor fields and 3-tensor fields can be uniquely (up to a constant) characterized by their invariance under sufficient statistics, thereby achieving a full generalization of the original result of Chentsov to infinite sample sizes. More generally, we decompose Markov morphisms between statistical models in terms of statistics. In particular, a monotonicity result for the Fisher information naturally follows.

preprint2013arXiv

Ollivier-Ricci curvature and the spectrum of the normalized graph Laplace operator

We prove the following estimate for the spectrum of the normalized Laplace operator $Δ$ on a finite graph $G$, \begin{equation*}1- (1- k[t])^{\frac{1}{t}}\leq λ_1 \leq \cdots \leq λ_{N-1}\leq 1+ (1- k[t])^{\frac{1}{t}}, \,\forall \,\,\text{integers}\,\, t\geq 1. \end{equation*} Here $k[t]$ is a lower bound for the Ollivier-Ricci curvature on the neighborhood graph $G[t]$, which was introduced by Bauer-Jost. In particular, when $t=1$ this is Ollivier's estimates $k\leq λ_1\leq \ldots \leq λ_{N-1}\leq 2-k$. For sufficiently large $t$ we show that, unless $G$ is bipartite, our estimates for $λ_1$ and $λ_{N-1}$ are always nontrivial and improve Ollivier's estimates for all graphs with $k\leq 0$. By definition neighborhood graphs are weighted graphs which may have loops. To understand the Ollivier-Ricci curvature on neighborhood graphs, we generalize a sharp estimate of the Ricci curvature given by Jost-Liu to weighted graphs with loops and relate it to the relative local frequency of triangles and loops.

preprint2013arXiv

Ollivier's Ricci curvature, local clustering and curvature dimension inequalities on graphs

In this paper, we explore the relationship between one of the most elementary and important properties of graphs, the presence and relative frequency of triangles, and a combinatorial notion of Ricci curvature. We employ a definition of generalized Ricci curvature proposed by Ollivier in a general framework of Markov processes and metric spaces and applied in graph theory by Lin-Yau. In analogy with curvature notions in Riemannian geometry, we interpret this Ricci curvature as a control on the amount of overlap between neighborhoods of two neighboring vertices. It is therefore naturally related to the presence of triangles containing those vertices, or more precisely, the local clustering coefficient, that is, the relative proportion of connected neighbors among all the neighbors of a vertex. This suggests to derive lower Ricci curvature bounds on graphs in terms of such local clustering coefficients. We also study curvature dimension inequalities on graphs, building upon previous work of several authors.

preprint2012arXiv

Analysis of inverse stochastic resonance and the long-term firing of Hodgkin-Huxley neurons with Gaussian white noise

In previous articles we have investigated the firing properties of the standard Hodgkin-Huxley (HH) systems of ordinary and partial differential equations in response to input currents composed of a drift (mean) and additive Gaussian white noise. For certain values of the mean current, as the noise amplitude increased from zero, the firing rate exhibited a minimum and this phenomenon was called inverse stochastic resonance (ISR). Here we analyse the underlying transitions from a stable equilibrium point to the limit cycle and vice-versa. Focusing on the case of a mean input current density $μ=6.8$ at which repetitive firing occurs and ISR had been found to be pronounced, some of the properties of the corresponding stable equilibrium point are found. A linearized approximation around this point has oscillatory solutions from whose maxima spikes tend to occur. A one dimensional diffusion is also constructed for small noise based on the correlations between the pairs of HH variables and the small magnitudes of the fluctuations in two of them. Properties of the basin of attraction of the limit cycle (spike) are investigated heuristically and also the nature of distribution of spikes at very small noise corresponding to trajectories which do not ever enter the basin of attraction of the equilibrium point. Long term trials of duration 500000 ms are carried out for values of the noise parameter $σ$ from 0 to 2.0, with results appearing in Section 3. The graph of mean spike count versus $σ$ is divided into 4 regions $R_1,...,R_4,$ where $R_3$ contains the minimum associated with ISR.

preprint2012arXiv

Bipartite and neighborhood graphs and the spectrum of the normalized graph Laplacian

We study the spectrum of the normalized Laplace operator of a connected graph $Γ$. As is well known, the smallest nontrivial eigenvalue measures how difficult it is to decompose $Γ$ into two large pieces, whereas the largest eigenvalue controls how close $Γ$ is to being bipartite. The smallest eigenvalue can be controlled by the Cheeger constant, and we establish a dual construction that controls the largest eigenvalue. Moreover, we find that the neighborhood graphs $Γ[l]$ of order $l\geq2$ encode important spectral information about $Γ$ itself which we systematically explore. In particular, the neighborhood graph method leads to new estimates for the smallest nontrivial eigenvalue that can improve the Cheeger inequality, as well as an explicit estimate for the largest eigenvalue from above and below. As applications of such spectral estimates, we provide a criterion for the synchronizability of coupled map lattices, and an estimate for the convergence rate of random walks on graphs.

preprint2012arXiv

Minimum vertex covers and the spectrum of the normalized Laplacian on trees

We show that, in the graph spectrum of the normalized graph Laplacian on trees, the eigenvalue 1 and eigenvalues near 1 are strongly related to minimum vertex covers. In particular, for the eigenvalue 1, its multiplicity is related to the size of a minimum vertex cover, and zero entries of its eigenvectors correspond to vertices in minimum vertex covers; while for eigenvalues near 1, their distance to 1 can be estimated from minimum vertex covers; and for the largest eigenvalue smaller than 1, the sign graphs of its eigenvectors take vertices in a minimum vertex cover as representatives.

preprint2012arXiv

Relations Between Graphs

Given two graphs G and H, we ask under which conditions there is a relation R that generates the edges of H given the structure of graph G. This construction can be seen as a form of multihomomorphism. It generalizes surjective homomorphisms of graphs and naturally leads to notions of R-retractions, R-cores, and R-cocores of graphs. Both R-cores and R-cocores of graphs are unique up to isomorphism and can be computed in polynomial time.

preprint2012arXiv

Shared Information -- New Insights and Problems in Decomposing Information in Complex Systems

How can the information that a set ${X_{1},...,X_{n}}$ of random variables contains about another random variable $S$ be decomposed? To what extent do different subgroups provide the same, i.e. shared or redundant, information, carry unique information or interact for the emergence of synergistic information? Recently Williams and Beer proposed such a decomposition based on natural properties for shared information. While these properties fix the structure of the decomposition, they do not uniquely specify the values of the different terms. Therefore, we investigate additional properties such as strong symmetry and left monotonicity. We find that strong symmetry is incompatible with the properties proposed by Williams and Beer. Although left monotonicity is a very natural property for an information measure it is not fulfilled by any of the proposed measures. We also study a geometric framework for information decompositions and ask whether it is possible to represent shared information by a family of posterior distributions. Finally, we draw connections to the notions of shared knowledge and common knowledge in game theory. While many people believe that independent variables cannot share information, we show that in game theory independent agents can have shared knowledge, but not common knowledge. We conclude that intuition and heuristic arguments do not suffice when arguing about information.

preprint2011arXiv

Interlacing inequalities for eigenvalues of discrete Laplace operators

The term interlacing refers to systematic inequalities between the sequences of eigenvalues of two operators defined on objects related by a specific oper- ation. In particular, knowledge of the spectrum of one of the objects then implies eigenvalue bounds for the other one. In this paper, we therefore develop topological arguments in order to de- rive such analytical inequalities. We investigate, in a general and systematic manner, interlacing of spectra for weighted simplicial complexes with arbi- trary weights. This enables us to control the spectral effects of operations like deletion of a subcomplex, collapsing and contraction of a simplex, cover- ings and simplicial maps, for absolute and relative Laplacians. It turns out that many well-known results from graph theory become special cases of our general results and consequently admit improvements and generalizations. In particular, we derive a number of effective eigenvalue bounds.

preprint2011arXiv

Spectra of combinatorial Laplace operators on simplicial complexes

We first develop a general framework for Laplace operators defined in terms of the combinatorial structure of a simplicial complex. This includes, among others, the graph Laplacian, the combinatorial Laplacian on simplicial complexes, the weighted Laplacian, and the normalized graph Laplacian. This framework then allows us to define the normalized Laplace operator $Δ_{i}^{up}$ on simplicial complexes which we then systematically investigate. We study the effects of a wedge sum, a join and a duplication of a motif on the spectrum of the normalized Laplace operator, and identify some of the combinatorial features of a simplicial complex that are encoded in its spectrum.

preprint2010arXiv

Genotype networks in metabolic reaction spaces

Background: A metabolic genotype comprises all chemical reactions an organism can catalyze via enzymes encoded in its genome. A genotype is viable in a given environment if it is capable of producing all biomass components the organism needs to survive and reproduce. Previous work has focused on the properties of individual genotypes while little is known about how genome-scale metabolic networks with a given function can vary in their reaction content. Results: We here characterize spaces of such genotypes. Specifically, we study metabolic genotypes whose phenotype is viability in minimal chemical environments that differ in their sole carbon sources. We show that regardless of the number of reactions in a metabolic genotype, the genotypes of a given phenotype typically form vast, connected, and unstructured sets -- genotype networks -- that nearly span the whole of genotype space. The robustness of metabolic phenotypes to random reaction removal in such spaces has a narrow distribution with a high mean. Different carbon sources differ in the number of metabolic genotypes in their genotype network; this number decreases as a genotype is required to be viable on increasing numbers of carbon sources, but much less than if metabolic reactions were used independently across different chemical environments. Conclusions: Our work shows that phenotype-preserving genotype networks have generic organizational properties and that these properties are insensitive to the number of reactions in metabolic genotypes.

preprint2010arXiv

Non-divergence harmonic maps

We describe work on solutions of certain non-divergence type and therefore non-variational elliptic and parabolic systems on manifolds. These systems include Hermitian and affine harmonics which should become useful tools for studying Hermitian and affine manifolds, resp. A key point is that in addition to the standard condition of nonpositive image curvature that is well known and understood in the theory of ordinary harmonic maps (which arise from a variational problem), here we also need in addition a global topological condition to guarantee the existence of solutions.

preprint2010arXiv

STDP-driven networks and the \emph{C. elegans} neuronal network

We study the dynamics of the structure of a formal neural network wherein the strengths of the synapses are governed by spike-timing-dependent plasticity (STDP). For properly chosen input signals, there exists a steady state with a residual network. We compare the motif profile of such a network with that of a real neural network of \emph{C. elegans} and identify robust qualitative similarities. In particular, our extensive numerical simulations show that this STDP-driven resulting network is robust under variations of the model parameters.

preprint2010arXiv

Symbolic dynamics and synchronization of coupled map networks with multiple delays

We use symbolic dynamics to study discrete-time dynamical systems with multiple time delays. We exploit the concept of avoiding sets, which arise from specific non-generating partitions of the phase space and restrict the occurrence of certain symbol sequences related to the characteristics of the dynamics. In particular, we show that the resulting forbidden sequences are closely related to the time delays in the system. We present two applications to coupled map lattices, namely (1) detecting synchronization and (2) determining unknown values of the transmission delays in networks with possibly directed and weighted connections and measurement noise. The method is applicable to multi-dimensional as well as set-valued maps, and to networks with time-varying delays and connection structure.

preprint2009arXiv

Affine Harmonic Maps

We introduce a class of maps from an affine flat into a Riemannian manifold that solve an elliptic system defined by the natural second order elliptic operator of the affine structure and the nonlinear Riemann geometry of the target. These maps are called affine harmonic. We show an existence result for affine harmonic maps in a given homotopy class when the target has non positive sectional curvature and some global non triviality condition is met. An example shows that such a condition is necessary. The analytical part is made difficult by the absence of a variational structure underlying affine harmonic maps. We therefore need to combine estimation techniques from geometric analysis and PDE theory with global geometric considerations.

preprint2007arXiv

Graph spectra as a systematic tool in computational biology

We present the spectrum of the (normalized) graph Laplacian as a systematic tool for the investigation of networks, and we describe basic properties of eigenvalues and eigenfunctions. Processes of graph formation like motif joining or duplication leave characteristic traces in the spectrum. This can suggest hypotheses about the evolution of a graph representing biological data. To this data, we analyze several biological networks in terms of rough qualitative data of their spectra.

preprint2007arXiv

On the spectrum of the normalized graph Laplacian

The spectrum of the normalized graph Laplacian yields a very comprehensive set of invariants of a graph. In order to understand the information contained in those invariants better, we systematically investigate the behavior of this spectrum under local and global operations like motif doubling, graph joining or splitting. The eigenvalue 1 plays a particular role, and we therefore emphasize those constructions that change its multiplicity in a controlled manner, like the iterated duplication of nodes.

preprint2007arXiv

Spectral plots and the representation and interpretation of biological data

It is basic question in biology and other fields to identify the char- acteristic properties that on one hand are shared by structures from a particular realm, like gene regulation, protein-protein interaction or neu- ral networks or foodwebs, and that on the other hand distinguish them from other structures. We introduce and apply a general method, based on the spectrum of the normalized graph Laplacian, that yields repre- sentations, the spectral plots, that allow us to find and visualize such properties systematically. We present such visualizations for a wide range of biological networks and compare them with those for networks derived from theoretical schemes. The differences that we find are quite striking and suggest that the search for universal properties of biological networks should be complemented by an understanding of more specific features of biological organization principles at different scales.