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

The orbit rigidity matrix of a symmetric framework

A number of recent papers have studied when symmetry causes frameworks on a graph to become infinitesimally flexible, or stressed, and when it has no impact. A number of other recent papers have studied special classes of frameworks on generically rigid graphs which are finite mechanisms. Here we introduce a new tool, the orbit matrix, which connects these two areas and provides a matrix representation for fully symmetric infinitesimal flexes, and fully symmetric stresses of symmetric frameworks. The orbit matrix is a true analog of the standard rigidity matrix for general frameworks, and its analysis gives important insights into questions about the flexibility and rigidity of classes of symmetric frameworks, in all dimensions. With this narrower focus on fully symmetric infinitesimal motions, comes the power to predict symmetry-preserving finite mechanisms - giving a simplified analysis which covers a wide range of the known mechanisms, and generalizes the classes of known mechanisms. This initial exploration of the properties of the orbit matrix also opens up a number of new questions and possible extensions of the previous results, including transfer of symmetry based results from Euclidean space to spherical, hyperbolic, and some other metrics with shared symmetry groups and underlying projective geometry.Comment: 41 pages, 12 figure

Internal states of model isotropic granular packings. I. Assembling process, geometry and contact networks

This is the first paper of a series of three, reporting on numerical simulation studies of geometric and mechanical properties of static assemblies of spherical beads under an isotropic pressure. Frictionless systems assemble in the unique random close packing (RCP) state in the low pressure limit if the compression process is fast enough, slower processes inducing traces of crystallization, and exhibit specific properties directly related to isostaticity of the force-carrying structure. The different structures of frictional packings assembled by various methods cannot be classified by the sole density. While lubricated systems approach RCP densities and coordination number z^*~=6 on the backbone in the rigid limit, an idealized "vibration" procedure results in equally dense configurations with z^*~=4.5. Near neighbor correlations on various scales are computed and compared to available laboratory data, although z^* values remain experimentally inaccessible. Low coordination packings have many rattlers (more than 10% of the grains carry no force), which should be accounted for on studying position correlations, and a small proportion of harmless "floppy modes" associated with divalent grains. Frictional packings, however slowly assembled under low pressure, retain a finite level of force indeterminacy, except in the limit of infinite friction.Comment: 29 pages. Published in Physical Review

Dissecting the Shared Genetic Architecture of Suicide Attempt, Psychiatric Disorders, and Known Risk Factors

Background Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders. Methods We conducted a genome-wide association study (GWAS) of 29,782 suicide attempt (SA) cases and 519,961 controls in the International Suicide Genetics Consortium (ISGC). The GWAS of SA was conditioned on psychiatric disorders using GWAS summary statistics via multitrait-based conditional and joint analysis, to remove genetic effects on SA mediated by psychiatric disorders. We investigated the shared and divergent genetic architectures of SA, psychiatric disorders, and other known risk factors. Results Two loci reached genome-wide significance for SA: the major histocompatibility complex and an intergenic locus on chromosome 7, the latter of which remained associated with SA after conditioning on psychiatric disorders and replicated in an independent cohort from the Million Veteran Program. This locus has been implicated in risk-taking behavior, smoking, and insomnia. SA showed strong genetic correlation with psychiatric disorders, particularly major depression, and also with smoking, pain, risk-taking behavior, sleep disturbances, lower educational attainment, reproductive traits, lower socioeconomic status, and poorer general health. After conditioning on psychiatric disorders, the genetic correlations between SA and psychiatric disorders decreased, whereas those with nonpsychiatric traits remained largely unchanged. Conclusions Our results identify a risk locus that contributes more strongly to SA than other phenotypes and suggest a shared underlying biology between SA and known risk factors that is not mediated by psychiatric disorders.Peer reviewe

Avoidance Expression in Rats as a Function of Signal-Shock Interval: Strain and Sex Differences

Inbred Wistar Kyoto (WKY) rats express inhibited temperament, increased sensitivity to stress, and exaggerated expressions of avoidance. A long-standing observation for lever press escape/avoidance learning in rats is the duration of the warning signal (WS) determines whether avoidance is expressed over escape. Outbred female Sprague-Dawley (SD) rats trained with a 10-s WS efficiently escaped, but failed to exhibit avoidance; avoidance was exhibited to a high degree with WSs longer than 20-s. We examined this longstanding WS duration function and extended it to male SD and male and female WKY rats. A cross-over design with two WS durations (10 s or 60 s) was employed. Rats were trained (20 trials/session) in four phases: acquisition (10 sessions), extinction (10 sessions), re-acquisition (8 sessions) and re-extinction (8 sessions). Consistent with the literature, female and male SD rats failed to express avoidance to an appreciable degree with a 10-s WS. When these rats were switched to a 60-s WS, performance levels in the initial session of training resembled the peak performance of rats trained with a 60-s WS. Therefore, the avoidance relationship was acquired, but not expressed at 10-s WS. Further, poor avoidance at 10-s does not adversely affect expression at 60-s. Failure to express avoidance with a 10-s WS likely reflects contrasting reinforcement value of avoidance, not a reduction in the amount of time available to respond or competing responses. In contrast, WKY rats exhibited robust avoidance with a 10-s WS, which was most apparent in female WKY rats. Exaggerated expression of avoidances by WKY rats, especially female rats, further confirms this inbred strain as a model of anxiety vulnerability