Source author record

Ofer Biham

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
13topics
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

20 published item(s)

preprint2022arXiv

The mean and variance of the distribution of shortest path lengths of random regular graphs

The distribution of shortest path lengths (DSPL) of random networks provides useful information on their large scale structure. In the special case of random regular graphs (RRGs), which consist of $N$ nodes of degree $c \ge 3$, the DSPL, denoted by $P(L=\ell)$, follows a discrete Gompertz distribution. Using the discrete Laplace transform we derive a closed-form expression for the moment generating function of the DSPL of RRGs. From the moment generating function we obtain closed-form expressions for the mean and variance of the DSPL. More specifically, we find that the mean distance between pairs of distinct nodes is given by $\langle L \rangle = \frac{\ln N}{\ln (c-1)} + \frac{1}{2} - \frac{ \ln c - \ln (c-2) +γ}{\ln (c-1)} + \mathcal{O} \left( \frac{\ln N}{N} \right)$, where $γ$ is the Euler-Mascheroni constant. While the leading term is known, this result includes a novel correction term, which yields very good agreement with the results obtained from direct numerical evaluation of $\langle L \rangle$ via the tail-sum formula and with the results obtained from computer simulations. However, it does not account for an oscillatory behavior of $\langle L \rangle$ as a function of $c$ or $N$. These oscillations are negligible in sparse networks but detectable in dense networks. We also derive an expression for the variance ${\rm Var}(L)$ of the DSPL, which captures the overall dependence of the variance on $c$ but does not account for the oscillations. The oscillations are due to the discrete nature of the shell structure around a random node. They reflect the profile of the filling of new shells as $N$ is increased. The results for the mean and variance are compared to the corresponding results obtained in other types of random networks. The relation between the mean distance and the diameter is discussed.

preprint2020arXiv

Convergence towards an Erd{\H o}s-Rényi graph structure in network contraction processes

In a highly influential paper twenty years ago, Barabási and Albert [Science 286, 509 (1999)] showed that networks undergoing generic growth processes with preferential attachment evolve towards scale-free structures. In any finite system, the growth eventually stalls and is likely to be followed by a phase of network contraction due to node failures, attacks or epidemics. Using the master equation formulation and computer simulations we analyze the structural evolution of networks subjected to contraction processes via random, preferential and propagating node deletions. We show that the contracting networks converge towards an Erd{\H o}s-Rényi network structure whose mean degree continues to decrease as the contraction proceeds. This is manifested by the convergence of the degree distribution towards a Poisson distribution and the loss of degree-degree correlations.

preprint2020arXiv

Statistical analysis of edges and bredges in configuration model networks

A bredge (bridge-edge) is an edge whose deletion would split the network component on which it resides into two components. Bredges are vulnerable links that play an important role in network collapse processes, which may result from node or link failures, attacks or epidemics. Therefore, the abundance and properties of bredges affect the resilience of the network. We present analytical results for the statistical properties of bredges in configuration model networks. Using a generating function approach based on the cavity method, we calculate the probability $\hat P(e\in{\rm B})$ that a random edge e in a configuration model network with degree distribution P(k) is a bredge (B). We also calculate the joint degree distribution $\hat P(k,k'|{\rm B})$ of the end-nodes of a random bredge. We examine the distinct properties of bredges on the giant component (GC) and on the finite tree components (FC) of the network. On the finite components all the edges are bredges and there are no degree-degree correlations. We calculate the probability $\hat P(e\in{\rm B}|{\rm GC})$ that a random edge on the giant component is a bredge. We also calculate the joint degree distribution $\hat P(k,k'|{\rm B},{\rm GC})$ of the end-nodes of bredges and the joint degree distribution $\hat P(k,k'|{\rm NB},{\rm GC})$ of the end-nodes of non-bredge (NB) edges on the giant component. Surprisingly, it is found that the degrees k and k' of the end-nodes of bredges are correlated, while the degrees of the end-nodes of NB edges are uncorrelated. We thus conclude that all the degree-degree correlations on the giant component are concentrated on the bredges. We calculate the covariance of end-nodes of bredges and show it is negative, namely bredges tend to connect high degree nodes to low degree nodes. The implications of the results are discussed in the context of common attack scenarios and dismantling processes.

preprint2016arXiv

Distance distribution in configuration model networks

We present analytical results for the distribution of shortest path lengths between random pairs of nodes in configuration model networks. The results, which are based on recursion equations, are shown to be in good agreement with numerical simulations for networks with degenerate, binomial and power-law degree distributions. The mean, mode and variance of the distribution of shortest path lengths are also evaluated. These results provide expressions for central measures and dispersion measures of the distribution of shortest path lengths in terms of moments of the degree distribution, illuminating the connection between the two distributions.

preprint2016arXiv

The distribution of path lengths of self avoiding walks on Erdős-Rényi networks

We present an analytical and numerical study of the paths of self avoiding walks (SAWs) on random networks. Since these walks do not retrace their paths, they effectively delete the nodes they visit, together with their links, thus pruning the network. The walkers hop between neighboring nodes, until they reach a dead-end node from which they cannot proceed. Focusing on Erdős-Rényi networks we show that the pruned networks maintain a Poisson degree distribution, $p_t(k)$, with an average degree, $\langle k \rangle_t$, that decreases linearly in time. We enumerate the SAW paths of any given length and find that the number of paths, $n_T(\ell)$, increases dramatically as a function of $\ell$. We also obtain analytical results for the path-length distribution, $P(\ell)$, of the SAW paths which are actually pursued, starting from a random initial node. It turns out that $P(\ell)$ follows the Gompertz distribution, which means that the termination probability of an SAW path increases with its length.

preprint2015arXiv

Analytical results for the distribution of shortest path lengths in random networks

We present two complementary analytical approaches for calculating the distribution of shortest path lengths in Erdos-Rényi networks, based on recursion equations for the shells around a reference node and for the paths originating from it. The results are in agreement with numerical simulations for a broad range of network sizes and connectivities. The average and standard deviation of the distribution are also obtained. In the case that the mean degree scales as $N^α$ with the network size, the distribution becomes extremely narrow in the asymptotic limit, namely almost all pairs of nodes are equidistant, at distance $d=\lfloor 1/α\rfloor$ from each other. The distribution of shortest path lengths between nodes of degree $m$ and the rest of the network is calculated. Its average is shown to be a monotonically decreasing function of $m$, providing an interesting relation between a local property and a global property of the network. The methodology presented here can be applied to more general classes of networks.

preprint2014arXiv

Stochastic analysis of bistability in coherent mixed feedback loops combining transcriptional and post-transcriptional regulations

Mixed feedback loops combining transcriptional and post-transcriptional regulations are common in cellular regulatory networks. They consist of two genes, encoding a transcription factor and a small non-coding RNA (sRNA), which mutually regulate each other's expression. We present a theoretical and numerical study of coherent mixed feedback loops of this type, in which both regulations are negative. Under suitable conditions, these feedback loops are expected to exhibit bistability, namely two stable states, one dominated by the transcriptional repressor and the other dominated by the sRNA. We use deterministic methods based on rate equation models, in order to identify the range of parameters in which bistability takes place. However, the deterministic models do not account for the finite lifetimes of the bistable states and the spontaneous, fluctuation-driven transitions between them. Therefore, we use stochastic methods to calculate the average lifetimes of the two states. It is found that these lifetimes strongly depend on rate coefficients such as the transcription rates of the transcriptional repressor and the sRNA. In particular, we show that the fraction of time the system spends in the sRNA dominated state follows a monotonically decreasing sigmoid function of the transcriptional repressor transcription rate. The biological relevance of these results is discussed in the context of such mixed feedback loops in {\it Escherichia coli}.

preprint2010arXiv

Binomial moment equations for stochastic reaction systems

A highly efficient formulation of moment equations for stochastic reaction networks is introduced. It is based on a set of binomial moments that capture the combinatorics of the reaction processes. The resulting set of equations can be easily runcated to include moments up to any desired order. The number of equations is dramatically reduced compared to the master equation. This formulation enables the simulation of complex reaction networks, involving a large number of reactive species much beyond the feasibility limit of any existing method. It provides an equation-based paradigm to the analysis of stochastic networks, complementing the commonly used Monte Carlo simulations.

preprint2010arXiv

Entanglement of Periodic States, the Quantum Fourier Transform and Shor's Factoring Algorithm

The preprocessing stage of Shor's algorithm generates a class of quantum states referred to as periodic states, on which the quantum Fourier transform is applied. Such states also play an important role in other quantum algorithms that rely on the quantum Fourier transform. Since entanglement is believed to be a necessary resource for quantum computational speedup, we analyze the entanglement of periodic states and the way it is affected by the quantum Fourier transform. To this end, we derive a formula that evaluates the Groverian entanglement measure for periodic states. Using this formula, we explain the surprising result that the Groverian entanglement of the periodic states built up during the preprocessing stage is only slightly affected by the quantum Fourier transform.

preprint2010arXiv

Interaction of Atomic and Molecular Hydrogen with Tholin Surfaces at Low Temperatures

We study the interaction of atomic and molecular hydrogen with a surface of tholin, a man-made polymer considered to be an analogue of aerosol particles present in Titan's atmosphere, using thermal programmed desorption at low temperatures below 30 K. The results are fitted and analyzed using a fine-grained rate equation model that describes the diffusion, reaction and desorption processes. We obtain the energy barriers for diffusion and desorption of atomic and molecular hydrogen. These barriers are found to be in the range of 30 to 60 meV, indicating that atom/molecule-surface interactions in this temperature range are dominated by weak adsorption forces. The implications of these results for the understanding of the atmospheric chemistry of Titan are discussed.

preprint2009arXiv

Accurate rate coefficients for models of interstellar gas-grain chemistry

The methodology for modeling grain-surface chemistry has been greatly improved by taking into account the grain size and fluctuation effects. However, the reaction rate coefficients currently used in all practical models of gas-grain chemistry are inaccurate by a significant amount. We provide expressions for these crucial rate coefficients that are both accurate and easy to incorporate into gas-grain models. We use exact results obtained in earlier work, where the reaction rate coefficient was defined by a first-passage problem, which was solved using random walk theory. The approximate reaction rate coefficient presented here is easy to include in all models of interstellar gas-grain chemistry. In contrast to the commonly used expression, the results that it provides are in perfect agreement with detailed kinetic Monte Carlo simulations. We also show the rate coefficient for reactions involving multiple species.

preprint2007arXiv

Molecular Hydrogen Formation on Amorphous Silicates Under Interstellar Conditions

Experimental results on the formation of molecular hydrogen on amorphous silicate surfaces are presented for the first time and analyzed using a rate equation model. The energy barriers for the relevant diffusion and desorption processes are obtained. They turn out to be significantly higher than those obtained earlier for polycrystalline silicates, demonstrating the importance of grain morphology. Using these barriers we evaluate the efficiency of molecular hydrogen formation on amorphous silicate grains under interstellar conditions. It is found that unlike polycrystalline silicates, amorphous silicate grains are efficient catalysts of H$_{2}$ formation within a temperature range which is relevant to diffuse interstellar clouds. The results also indicate that the hydrogen molecules are thermalized with the surface and desorb with low kinetic energy. Thus, they are unlikely to occupy highly excited states.

preprint2005arXiv

Formation of molecular hydrogen on analogues of interstellar dust grains: experiments and modelling

Molecular hydrogen has an important role in the early stages of star formation as well as in the production of many other molecules that have been detected in the interstellar medium. In this review we show that it is now possible to study the formation of molecular hydrogen in simulated astrophysical environments. Since the formation of molecular hydrogen is believed to take place on dust grains, we show that surface science techniques such as thermal desorption and time-of-flight can be used to measure the recombination efficiency, the kinetics of reaction and the dynamics of desorption. The analysis of the experimental results using rate equations gives useful insight on the mechanisms of reaction and yields values of parameters that are used in theoretical models of interstellar cloud chemistry.

preprint2000arXiv

Analysis of Generalized Grover's Quantum Search Algorithms Using Recursion Equations

The recursion equation analysis of Grover's quantum search algorithm presented by Biham et al. [PRA 60, 2742 (1999)] is generalized. It is applied to the large class of Grover's type algorithms in which the Hadamard transform is replaced by any other unitary transformation and the phase inversion is replaced by a rotation by an arbitrary angle. The time evolution of the amplitudes of the marked and unmarked states, for any initial complex amplitude distribution is expressed using first order linear difference equations. These equations are solved exactly. The solution provides the number of iterations T after which the probability of finding a marked state upon measurement is the highest, as well as the value of this probability, P_max. Both T and P_max are found to depend on the averages and variances of the initial amplitude distributions of the marked and unmarked states, but not on higher moments.

preprint1999arXiv

Electromigration-Induced Flow of Islands and Voids on the Cu(001) Surface

Electromigration-induced flow of islands and voids on the Cu(001) surface is studied at the atomic scale. The basic drift mechanisms are identified using a complete set of energy barriers for adatom hopping on the Cu(001) surface, combined with kinetic Monte Carlo simulations. The energy barriers are calculated by the embedded atom method, and parameterized using a simple model. The dependence of the flow on the temperature, the size of the clusters, and the strength of the applied field is obtained. For both islands and voids it is found that edge diffusion is the dominant mass-transport mechanism. The rate limiting steps are identified. For both islands and voids they involve detachment of atoms from corners into the adjacent edge. The energy barriers for these moves are found to be in good agreement with the activation energy for island/void drift obtained from Arrhenius analysis of the simulation results. The relevance of the results to other FCC(001) metal surfaces and their experimental implications are discussed.

preprint1999arXiv

Grover's Quantum Search Algorithm for an Arbitrary Initial Amplitude Distribution

Grover's algorithm for quantum searching is generalized to deal with arbitrary initial complex amplitude distributions. First order linear difference equations are found for the time evolution of the amplitudes of the marked and unmarked states. These equations are solved exactly. New expressions are derived for the optimal time of measurement and the maximal probability of success. They are found to depend on the averages and variances of the initial amplitude distributions of the marked and unmarked states, but not on higher moments. Our results imply that Grover's algorithm is robust against modest noise in the amplitude initialization procedure.

preprint1998arXiv

Generalized Grover Search Algorithm for Arbitrary Initial Amplitude Distribution

Grover's algorithm for quantum searching of a database is generalized to deal with arbitrary initial amplitude distributions. First order linear difference equations are found for the time evolution of the amplitudes of the r marked and N-r unmarked states. These equations are solved exactly. An expression for the optimal measurement time T \sim O(\sqrt{N/r}) is derived which is shown to depend only on the initial average amplitudes of the marked and unmarked states. A bound on the probability of measuring a marked state is derived, which depends only on the standard deviation of the initial amplitude distributions of the marked or unmarked states.

preprint1998arXiv

Scaling Range and Cutoffs in Empirical Fractals

Fractal structures appear in a vast range of physical systems. A literature survey including all experimental papers on fractals which appeared in the six Physical Review journals (A-E and Letters) during the 1990's shows that experimental reports of fractal behavior are typically based on a scaling range $Δ$ which spans only 0.5 - 2 decades. This range is limited by upper and lower cutoffs either because further data is not accessible or due to crossover bends. Focusing on spatial fractals, a classification is proposed into (a) aggregation; (b) porous media; (c) surfaces and fronts; (d) fracture and (e) critical phenomena. Most of these systems, [except for class (e)] involve processes far from thermal equilibrium. The fact that for self similar fractals [in contrast to the self affine fractals of class (c)] there are hardly any exceptions to the finding of $Δ\le 2$ decades, raises the possibility that the cutoffs are due to intrinsic properties of the measured systems rather than the specific experimental conditions and apparatus. To examine the origin of the limited range we focus on a class of aggregation systems. In these systems a molecular beam is deposited on a surface, giving rise to nucleation and growth of diffusion-limited-aggregation-like clusters. Scaling arguments are used to show that the required duration of the deposition experiment increases exponentially with $Δ$. Furthermore, using realistic parameters for surfaces such as Al(111) it is shown that these considerations limit the range of fractal behavior to less than two decades in agreement with the experimental findings. It is conjectured that related kinetic mechanisms that limit the scaling range are common in other nonequilibrium processes which generate spatial fractals.

preprint1997arXiv

Limited Range Fractality of Randomly Adsorbed Rods

Multiple resolution analysis of two dimensional structures composed of randomly adsorbed penetrable rods, for densities below the percolation threshold, has been carried out using box-counting functions. It is found that at relevant resolutions, for box-sizes, $r$, between cutoffs given by the average rod length $<\ell>$ and the average inter-rod distance $r_1$, these systems exhibit apparent fractal behavior. It is shown that unlike the case of randomly distributed isotropic objects, the upper cutoff $r_1$ is not only a function of the coverage but also depends on the excluded volume, averaged over the orientational distribution. Moreover, the apparent fractal dimension also depends on the orientational distributions of the rods and decreases as it becomes more anisotropic. For box sizes smaller than $<\ell>$ the box counting function is determined by the internal structure of the rods, whether simple or itself fractal. Two examples are considered - one of regular rods of one dimensional structure and rods which are trimmed into a Cantor set structure which are fractals themselves. The models examined are relevant to adsorption of linear molecules and fibers, liquid crystals, stress induced fractures and edge imperfections in metal catalysts. We thus obtain a distinction between two ranges of length scales: $r < <\ell>$ where the internal structure of the adsorbed objects is probed, and $<\ell> < r < r_1$ where their distribution is probed, both of which may exhibit fractal behavior. This distinction is relevant to the large class of systems which exhibit aggregation of a finite density of fractal-like clusters, which includes surface growth in molecular beam epitaxy and diffusion-limited-cluster-cluster-aggregation models.