Directed force chain networks and stress response in static granular materials

A theory of stress fields in two-dimensional granular materials based on directed force chain networks is presented. A general equation for the densities of force chains in different directions is proposed and a complete solution is obtained for a special case in which chains lie along a discrete set of directions. The analysis and results demonstrate the necessity of including nonlinear terms in the equation. A line of nontrivial fixed point solutions is shown to govern the properties of large systems. In the vicinity of a generic fixed point, the response to a localized load shows a crossover from a single, centered peak at intermediate depths to two propagating peaks at large depths that broaden diffusively.Comment: 18 pages, 12 figures. Minor corrections to one figur

Random trees between two walls: Exact partition function

We derive the exact partition function for a discrete model of random trees embedded in a one-dimensional space. These trees have vertices labeled by integers representing their position in the target space, with the SOS constraint that adjacent vertices have labels differing by +1 or -1. A non-trivial partition function is obtained whenever the target space is bounded by walls. We concentrate on the two cases where the target space is (i) the half-line bounded by a wall at the origin or (ii) a segment bounded by two walls at a finite distance. The general solution has a soliton-like structure involving elliptic functions. We derive the corresponding continuum scaling limit which takes the remarkable form of the Weierstrass p-function with constrained periods. These results are used to analyze the probability for an evolving population spreading in one dimension to attain the boundary of a given domain with the geometry of the target (i) or (ii). They also translate, via suitable bijections, into generating functions for bounded planar graphs.Comment: 25 pages, 7 figures, tex, harvmac, epsf; accepted version; main modifications in Sect. 5-6 and conclusio

Random tensor models in the large N limit: Uncoloring the colored tensor models

Tensor models generalize random matrix models in yielding a theory of dynamical triangulations in arbitrary dimensions. Colored tensor models have been shown to admit a 1/N expansion and a continuum limit accessible analytically. In this paper we prove that these results extend to the most general tensor model for a single generic, i.e. non-symmetric, complex tensor. Colors appear in this setting as a canonical book-keeping device and not as a fundamental feature. In the large N limit, we exhibit a set of Virasoro constraints satisfied by the free energy and an infinite family of multicritical behaviors with entropy exponents \gamma_m=1-1/m.Comment: 15 page

Multicritical continuous random trees

We introduce generalizations of Aldous' Brownian Continuous Random Tree as scaling limits for multicritical models of discrete trees. These discrete models involve trees with fine-tuned vertex-dependent weights ensuring a k-th root singularity in their generating function. The scaling limit involves continuous trees with branching points of order up to k+1. We derive explicit integral representations for the average profile of this k-th order multicritical continuous random tree, as well as for its history distributions measuring multi-point correlations. The latter distributions involve non-positive universal weights at the branching points together with fractional derivative couplings. We prove universality by rederiving the same results within a purely continuous axiomatic approach based on the resolution of a set of consistency relations for the multi-point correlations. The average profile is shown to obey a fractional differential equation whose solution involves hypergeometric functions and matches the integral formula of the discrete approach.Comment: 34 pages, 12 figures, uses lanlmac, hyperbasics, eps

On static shells and the Buchdahl inequality for the spherically symmetric Einstein-Vlasov system

In a previous work \cite{An1} matter models such that the energy density $\rho\geq 0,$ and the radial- and tangential pressures $p\geq 0$ and $q,$ satisfy $p+q\leq\Omega\rho, \Omega\geq 1,$ were considered in the context of Buchdahl's inequality. It was proved that static shell solutions of the spherically symmetric Einstein equations obey a Buchdahl type inequality whenever the support of the shell, $[R_0,R_1], R_0>0,$ satisfies $R_1/R_0<1/4.$ Moreover, given a sequence of solutions such that $R_1/R_0\to 1,$ then the limit supremum of $2M/R_1$ was shown to be bounded by $((2\Omega+1)^2-1)/(2\Omega+1)^2.$ In this paper we show that the hypothesis that $R_1/R_0\to 1,$ can be realized for Vlasov matter, by constructing a sequence of static shells of the spherically symmetric Einstein-Vlasov system with this property. We also prove that for this sequence not only the limit supremum of $2M/R_1$ is bounded, but that the limit is $((2\Omega+1)^2-1)/(2\Omega+1)^2=8/9,$ since $\Omega=1$ for Vlasov matter. Thus, static shells of Vlasov matter can have $2M/R_1$ arbitrary close to $8/9,$ which is interesting in view of \cite{AR2}, where numerical evidence is presented that 8/9 is an upper bound of $2M/R_1$ of any static solution of the spherically symmetric Einstein-Vlasov system.Comment: 20 pages, Late

Rate-dependent propagation of cardiac action potentials in a one-dimensional fiber

Action potential duration (APD) restitution, which relates APD to the preceding diastolic interval (DI), is a useful tool for predicting the onset of abnormal cardiac rhythms. However, it is known that different pacing protocols lead to different APD restitution curves (RCs). This phenomenon, known as APD rate-dependence, is a consequence of memory in the tissue. In addition to APD restitution, conduction velocity restitution also plays an important role in the spatiotemporal dynamics of cardiac tissue. We present new results concerning rate-dependent restitution in the velocity of propagating action potentials in a one-dimensional fiber. Our numerical simulations show that, independent of the amount of memory in the tissue, waveback velocity exhibits pronounced rate-dependence and the wavefront velocity does not. Moreover, the discrepancy between waveback velocity RCs is most significant for small DI. We provide an analytical explanation of these results, using a system of coupled maps to relate the wavefront and waveback velocities. Our calculations show that waveback velocity rate-dependence is due to APD restitution, not memory.Comment: 17 pages, 7 figure

Critical collapse of collisionless matter - a numerical investigation

In recent years the threshold of black hole formation in spherically symmetric gravitational collapse has been studied for a variety of matter models. In this paper the corresponding issue is investigated for a matter model significantly different from those considered so far in this context. We study the transition from dispersion to black hole formation in the collapse of collisionless matter when the initial data is scaled. This is done by means of a numerical code similar to those commonly used in plasma physics. The result is that for the initial data for which the solutions were computed, most of the matter falls into the black hole whenever a black hole is formed. This results in a discontinuity in the mass of the black hole at the onset of black hole formation.Comment: 22 pages, LaTeX, 7 figures (ps-files, automatically included using psfig