Source author record

Cristian F. Moukarzel

Cristian F. Moukarzel 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

10works
10topics
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

10 published item(s)

preprint2019arXiv

Sustained rotation in a vibrated disk with asymmetric supports

A single frictional elastic disk, supported against gravity by two others, rotates steadily when the supports are vibrated and the system is tilted with respect to gravity. Rotation is here studied using Molecular Dynamics Simulations, and a detailed analysis of the dynamics of the system is made. The origin of the observed rotational ratcheting is discussed by considering simplified situations analytically. This shows that the sense of rotation is not fixed by the tilt but depends on the details of the excitation as well.

preprint2015arXiv

Comment on "Penrose Tilings as Jammed Solids"

In a recent letter, Stenull and Lubensky claim that periodic approximants of Penrose tilings, which are generically isostatic, have a nonzero bulk modulus B when disordered, and, therefore, Penrose tilings are good models of jammed packings. The claim of a nonzero B, which is made on the basis of a normal mode analysis of periodic Penrose approximants for a single value of the disorder epsilon, is the central point of their letter: other properties of Penrose tilings, such as the vanishing of the shear modulus, and a flat density of vibrational states, are already shared by most geometrically disordered isostatic networks studied so far. In this comment, Conjugate Gradient is used to solve the elastic equations on approximants with up to 8x10^4 sites for several values of epsilon, to show beyond reasonable doubt that Stenull and Lubensky's claim is incorrect. The bulk modulus of generic Penrose tilings is zero asymptotically. According to our results, B grows as (epsilon^2 L^3) when (epsilon^2 L^3) << 10^2, then saturates, and finally decays as (epsilon^2 L^3)^{-2/3} ~ 1/L^2 for epsilon^2 L^3 >> 10^2. Stenull and Lubensky seem to have only analyzed one value of epsilon for which saturation is reached at the largest size studied. This led them to a wrong conclusion. We support our results by also considering generic Penrose approximants with fixed boundaries, whose bulk modulus constitutes a strict upper bound for that of periodic systems, finding that these have a vanishing B as well for large L. We conclude that the main point in Stenull and Lubensky letter is unjustified. Penrose tilings are no better models of jammed packings than any of the previously studied isostatic networks with geometric disorder.

preprint2013arXiv

Two rigidity percolation transitions on binary Bethe networks and the intermediate phase in glass

Rigidity Percolation is studied analytically on randomly bonded networks with two types of nodes, respectively with coordination numbers $z_1$ and $z_2$, and with $g_1$ and $g_2$ degrees of freedom each. For certain cases that model chalcogenide glass networks, two transitions, both of first order, are found, with the first transition usually rather weak. The ensuing intermediate pase, although is not isostatic in its entirety, has very low self-stress. Our results suggest a possible mechanism for the appearance of intermediate phases in glass, that does not depend on a self-organization principle.

preprint2012arXiv

Elastic collapse in disordered isostatic networks

Isostatic networks are minimally rigid and therefore have, generically, nonzero elastic moduli. Regular isostatic networks have finite moduli in the limit of large sizes. However, numerical simulations show that all elastic moduli of geometrically disordered isostatic networks go to zero with system size. This holds true for positional as well as for topological disorder. In most cases, elastic moduli decrease as inverse power-laws of system size. On directed isostatic networks, however, of which the square and cubic lattices are particular cases, the decrease of the moduli is exponential with size. For these, the observed elastic weakening can be quantitatively described in terms of the multiplicative growth of stresses with system size, giving rise to bulk and shear moduli of order exp{-bL}. The case of sphere packings, which only accept compressive contact forces, is considered separately. It is argued that these have a finite bulk modulus because of specific correlations in contact disorder, introduced by the constraint of compressivity. We discuss why their shear modulus, nevertheless, is again zero for large sizes. A quantitative model is proposed that describes the numerically measured shear modulus, both as a function of the loading angle and system size. In all cases, if a density p>0 of overconstraints is present, as when a packing is deformed by compression, or when a glass is outside its isostatic composition window, all asymptotic moduli become finite. For square networks with periodic boundary conditions, these are of order sqrt{p}. For directed networks, elastic moduli are of order exp{-c/p}, indicating the existence of an "isostatic length scale" of order 1/p.

preprint2012arXiv

Yard-Sale exchange on networks: Wealth sharing and wealth appropriation

Yard-Sale (YS) is a stochastic multiplicative wealth-exchange model with two phases: a stable one where wealth is shared, and an unstable one where wealth condenses onto one agent. YS is here studied numerically on 1d rings, 2d square lattices, and random graphs with variable average coordination, comparing its properties with those in mean field (MF). Equilibrium properties in the stable phase are almost unaffected by the introduction of a network. Measurement of decorrelation times in the stable phase allow us to determine the critical interface with very good precision, and it turns out to be the same, for all networks analyzed, as the one that can be analytically derived in MF. In the unstable phase, on the other hand, dynamical as well as asymptotic properties are strongly network-dependent. Wealth no longer condenses on a single agent, as in MF, but onto an extensive set of agents, the properties of which depend on the network. Connections with previous studies of coalescence of immobile reactants are discussed, and their analytic predictions are successfully compared with our numerical results.

preprint2011arXiv

Multiplicative Asset Exchange with Arbitrary Return Distributions

The conservative wealth-exchange process derived from trade interactions is modeled as a multiplicative stochastic transference of value, where each interaction multiplies the wealth of the poorest of the two intervening agents by a random gain eta=(1+kappa), with kappa a random return. Analyzing the kinetic equation for the wealth distribution P(w,t), general properties are derived for arbitrary return distributions pi(kappa). If the geometrical average of the gain is larger than one, i.e. if <ln eta> >0, in the long time limit a nontrivial equilibrium wealth distribution P(w) is attained. Whenever <ln eta> <0, on the other hand, Wealth Condensation occurs, meaning that a single agent gets the whole wealth in the long run. This concentration phenomenon happens even if the average return <kappa> of the poor agent is positive. In the stable phase, P(w) behaves as w^{(T-1)} for w -> 0, and we find T exactly. This exponent is nonzero in the stable phase but goes to zero on approach to the condensation interface. The exact wealth distribution can be obtained analytically for the particular case of Kelly betting, and it turns out to be exponential. We show, however, that our model is never reversible, no matter what pi(kappa) is. In the condensing phase, the wealth of an agent with relative rank x is found to be w(x,t) \sim e^{x t <ln eta>} for finite times t. The wealth distribution is consequently P(w) \sim 1/w for finite times, while all wealth ends up in the hands of a single agent for large times. Numerical simulations are carried out, and found to satisfactorily compare with the above mentioned analytic results.

preprint2002arXiv

Shortest paths on systems with power-law distributed long-range connections

We discuss shortest-path lengths $\ell(r)$ on periodic rings of size L supplemented with an average of pL randomly located long-range links whose lengths are distributed according to $P_l \sim l^{-\xpn}$. Using rescaling arguments and numerical simulation on systems of up to $10^7$ sites, we show that a characteristic length $ξ$ exists such that $\ell(r) \sim r$ for $r<ξ$ but $\ell(r) \sim r^{θ_s(\xpn)}$ for $r>>ξ$. For small p we find that the shortest-path length satisfies the scaling relation $\ell(r,\xpn,p)/ξ= f(\xpn,r/ξ)$. Three regions with different asymptotic behaviors are found, respectively: a) $\xpn>2$ where $θ_s=1$, b) $1<\xpn<2$ where $0<θ_s(\xpn)<1/2$ and, c) $\xpn<1$ where $\ell(r)$ behaves logarithmically, i.e. $θ_s=0$. The characteristic length $ξ$ is of the form $ξ\sim p^{-ν}$ with $ν=1/(2-\xpn)$ in region b), but depends on L as well in region c). A directed model of shortest-paths is solved and compared with numerical results.

preprint1996arXiv

An efficient algorithm for two-dimensional central-force rigidity

Given a structure made up of n sites connected by b bars, the problem of recognizing which subsets of sites form rigid units is not a trivial one, because of the non-local character of rigidity in central-force systems. Even though this is a very old problem of statics, no simple algorithms are available for it so the most usual approach has been to solve the elastic equations, which is very time-consuming for large systems. Recently an integer algorithm was proposed for this problem in two dimensions, using matching methods from graph theory and Laman's theorem for two-dimensional graphs. The method is relatively simple, but its time complexity grows as $n^2$ in the worst case, and almost as fast on practical cases, so that an improvement is highly desirable. I describe here a further elaboration of that procedure, which relies upon the description of the system as a collection of rigid bodies connected by bars, instead of sites connected by bars. Sets of rigidly connected objects are replaced by a unique body, and this is done recursively as more rigid connections between bodies are discovered at larger scales. As a consequence of this ``rescaling transformation'', our algorithm has much improved average behavior, even when its worst-case complexity remains $n^2$. The time complexity of the body-bar algorithm is found to scale as $n^{1.12}$ for the randomly diluted triangular lattice, while the original site-bar version scales as $n^{1.9}$ on the same problem.

preprint1996arXiv

Infinite-cluster geometry in central-force networks

We show that the infinite percolating cluster (with density P_inf) of central-force networks is composed of: a fractal stress-bearing backbone (Pb) and; rigid but unstressed ``dangling ends'' which occupy a finite volume-fraction of the lattice (Pd). Near the rigidity threshold pc, there is then a first-order transition in P_inf = Pd + Pb, while Pb is second-order with exponent Beta'. A new mean field theory shows Beta'(mf)=1/2, while simulations of triangular lattices give Beta'_tr = 0.255 +/- 0.03.