Algorithms for 3D rigidity analysis and a first order percolation transition

A fast computer algorithm, the pebble game, has been used successfully to study rigidity percolation on 2D elastic networks, as well as on a special class of 3D networks, the bond-bending networks. Application of the pebble game approach to general 3D networks has been hindered by the fact that the underlying mathematical theory is, strictly speaking, invalid in this case. We construct an approximate pebble game algorithm for general 3D networks, as well as a slower but exact algorithm, the relaxation algorithm, that we use for testing the new pebble game. Based on the results of these tests and additional considerations, we argue that in the particular case of randomly diluted central-force networks on BCC and FCC lattices, the pebble game is essentially exact. Using the pebble game, we observe an extremely sharp jump in the largest rigid cluster size in bond-diluted central-force networks in 3D, with the percolating cluster appearing and taking up most of the network after a single bond addition. This strongly suggests a first order rigidity percolation transition, which is in contrast to the second order transitions found previously for the 2D central-force and 3D bond-bending networks. While a first order rigidity transition has been observed for Bethe lattices and networks with ``chemical order'', this is the first time it has been seen for a regular randomly diluted network. In the case of site dilution, the transition is also first order for BCC, but results for FCC suggest a second order transition. Even in bond-diluted lattices, while the transition appears massively first order in the order parameter (the percolating cluster size), it is continuous in the elastic moduli. This, and the apparent non-universality, make this phase transition highly unusual.Comment: 28 pages, 19 figure

Hierarchical models of rigidity percolation

We introduce models of generic rigidity percolation in two dimensions on hierarchical networks, and solve them exactly by means of a renormalization transformation. We then study how the possibility for the network to self organize in order to avoid stressed bonds may change the phase diagram. In contrast to what happens on random graphs and in some recent numerical studies at zero temperature, we do not find a true intermediate phase separating the usual rigid and floppy ones.Comment: 20 pages, 8 figures. Figures improved, references added, small modifications. Accepted in Phys. Rev.

Self-organization with equilibration: a model for the intermediate phase in rigidity percolation

Recent experimental results for covalent glasses suggest the existence of an intermediate phase attributed to the self-organization of the glass network resulting from the tendency to minimize its internal stress. However, the exact nature of this experimentally measured phase remains unclear. We modify a previously proposed model of self-organization by generating a uniform sampling of stress-free networks. In our model, studied on a diluted triangular lattice, an unusual intermediate phase appears, in which both rigid and floppy networks have a chance to occur, a result also observed in a related model on a Bethe lattice by Barre et al. [Phys. Rev. Lett. 94, 208701 (2005)]. Our results for the bond-configurational entropy of self-organized networks, which turns out to be only about 2% lower than that of random networks, suggest that a self-organized intermediate phase could be common in systems near the rigidity percolation threshold.Comment: 9 pages, 6 figure

Variation of the glass transition temperature with rigidity and chemical composition

The effects of flexibility and chemical composition in the variation of the glass transition temperature are obtained by using the Lindemann criteria, that relates melting temperature with atomic vibrations. Using this criteria and that floppy modes at low frequencies enhance in a considerable way the average cuadratic displacement, we show that the consequence is a modified glass transition temperature. This approach allows to obtain in a simple way the empirically modified Gibbs-DiMarzio law, which has been widely used in chalcogenide glasses to fit the changes in the glass transition temperature with the chemical composition . The method predicts that the constant that appears in the law depends upon the ratio of two characteristic frequencies (or temperatures). Then, the constant for the Se-Ge-As glass is estimated by using the experimental density of vibrational states, and the result shows a very good agreement with the experimental fit from glass transition temperature variation

Self-organized criticality in the intermediate phase of rigidity percolation

Experimental results for covalent glasses have highlighted the existence of a new self-organized phase due to the tendency of glass networks to minimize internal stress. Recently, we have shown that an equilibrated self-organized two-dimensional lattice-based model also possesses an intermediate phase in which a percolating rigid cluster exists with a probability between zero and one, depending on the average coordination of the network. In this paper, we study the properties of this intermediate phase in more detail. We find that microscopic perturbations, such as the addition or removal of a single bond, can affect the rigidity of macroscopic regions of the network, in particular, creating or destroying percolation. This, together with a power-law distribution of rigid cluster sizes, suggests that the system is maintained in a critical state on the rigid/floppy boundary throughout the intermediate phase, a behavior similar to self-organized criticality, but, remarkably, in a thermodynamically equilibrated state. The distinction between percolating and non-percolating networks appears physically meaningless, even though the percolating cluster, when it exists, takes up a finite fraction of the network. We point out both similarities and differences between the intermediate phase and the critical point of ordinary percolation models without self-organization. Our results are consistent with an interpretation of recent experiments on the pressure dependence of Raman frequencies in chalcogenide glasses in terms of network homogeneity.Comment: 20 pages, 18 figure

Evaluation and management implications of uncertainty in a multispecies size-structured model of population and community responses to fishing

1. Implementation of an ecosystem approach to fisheries requires advice on trade-offs among fished species and between fisheries yields and biodiversity or food web properties. However, the lack of explicit representation, analysis and consideration of uncertainty in most multispecies models has limited their application in analyses that could support management advice. 2. We assessed the consequences of parameter uncertainty by developing 78 125 multispecies size-structured fish community models, with all combinations of parameters drawn from ranges that spanned parameter values estimated from data and literature. This unfiltered ensemble was reduced to 188 plausible models, the filtered ensemble (FE), by screening outputs against fish abundance data and ecological principles such as requiring species' persistence. 3. Effects of parameter uncertainty on estimates of single-species management reference points for fishing mortality (FMSY, fishing mortality rate providing MSY, the maximum sustainable yield) and biomass (BMSY, biomass at MSY) were evaluated by calculating probability distributions of estimated reference points with the FE. There was a 50% probability that multispecies FMSY could be estimated to within ±25% of its actual value, and a 50% probability that BMSY could be estimated to within ±40% of its actual value. 4. Signal-to-noise ratio was assessed for four community indicators when mortality rates were reduced from current rates to FMSY. The slope of the community size spectrum showed the greatest signal-to-noise ratio, indicating that it would be the most responsive indicator to the change in fishing mortality F. Further, the power of an ongoing international monitoring survey to detect predicted responses of size spectrum slope was higher than for other size-based metrics. 5. Synthesis and applications: Application of the ensemble model approach allows explicit representation of parameter uncertainty and supports advice and management by (i) providing uncertainty intervals for management reference points, (ii) estimating working values of reference points that achieve a defined reduction in risk of not breaching the true reference point, (iii) estimating the responsiveness of population, community, food web and biodiversity indicators to changes in F, (iv) assessing the performance of indicators and monitoring programmes and (v) identifying priorities for data collection and changes to model structure to reduce uncertainty

Floppy modes and the free energy: Rigidity and connectivity percolation on Bethe Lattices

We show that negative of the number of floppy modes behaves as a free energy for both connectivity and rigidity percolation, and we illustrate this result using Bethe lattices. The rigidity transition on Bethe lattices is found to be first order at a bond concentration close to that predicted by Maxwell constraint counting. We calculate the probability of a bond being on the infinite cluster and also on the overconstrained part of the infinite cluster, and show how a specific heat can be defined as the second derivative of the free energy. We demonstrate that the Bethe lattice solution is equivalent to that of the random bond model, where points are joined randomly (with equal probability at all length scales) to have a given coordination, and then subsequently bonds are randomly removed.Comment: RevTeX 11 pages + epsfig embedded figures. Submitted to Phys. Rev.

Reinforced structural plastics

Reinforced polyimide structures are described. Reinforcing materials are impregnated with a suspension of polyimide prepolymer and bonded together by heat and pressure to form a cured, hard-reinforced, polyimide structure

Rigidity transitions and constraint counting in amorphous networks: beyond the mean-field approach

Subj-class: Disordered Systems and Neural NetworksComment: 12 pages, revtex, 3 figure

Critical scaling and heterogeneous superdiffusion across the jamming/rigidity transition of a granular glass

The dynamical properties of a dense horizontally vibrated bidisperse granular monolayer are experimentally investigated. The quench protocol produces states with a frozen structure of the assembly, but the remaining degrees of freedom associated with contact dynamics control the appearance of macroscopic rigidity. We provide decisive experimental evidence that this transition is a critical phenomenon, with increasingly collective and heterogeneous rearrangements occurring at length scales much smaller than the grains' diameter, presumably reflecting the contact force network fluctuations. Dynamical correlation time and length scales soar on both sides of the transition, as the volume fraction varies over a remarkably tiny range ($\delta \phi/\phi \sim 10^{-3}$). We characterize the motion of individual grains, which becomes super-diffusive at the jamming transition $\phi_J$, signaling long-ranged temporal correlations. Correspondingly, the system exhibits long-ranged four-point dynamical correlations in space that obey critical scaling at the transition density.Comment: 4 pages, 8 figure