BZPEER
More

Discover

About BZPEERSearchGraph

Research tools

HomeCollectionsResearchCollaborationPublish

Network

ExpertsCommunities

Opportunities

Listings

Account

Log inCreate account
Log inCreate account

Source author record

Corey S. O'Hern

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechcond-mat.mtrl-sciBiomoleculescond-mat.dis-nnBiological Physicsphysics.comp-phchao-dyncomp-gascond-matcond-mat.mes-hallEmerging TechnologiesNeural and Evolutionary Computingnlin.CDnlin.CGnlin.PSpatt-solphysics.data-anphysics.opticsQuantitative Methods

Catalog footprint

Research overview

48works
20topics
4close collaborators

Actions

Is this your work?

Log in to claim
Open graphBrowse papers

Published work

48 publications

preprint2026arXiv

Assessment of scoring functions for computational models of protein-protein interfaces

A goal of computational studies of protein-protein interfaces (PPIs) is to predict the binding site between two monomers that form a heterodimer. The simplest version of this problem is to rigidly re-dock the bound forms of the monomers, which involves generating computational models of the heterodimer and then scoring them to determine the most native-like models. Scoring functions have been assessed previously using rank- and classification-based metrics, however, these methods are sensitive to the number and quality of models in the scoring function training set. We assess the accuracy of seven PPI scoring functions by comparing their scores to a measure of structural similarity to the x-ray crystal structure (i.e. the DockQ score) for a non-redundant set of heterodimers from the Protein Data Bank. For each heterodimer, we generate re-docked models uniformly sampled over DockQ and calculate the Spearman correlation between the PPI scores and DockQ. For some targets, the scores and DockQ are highly correlated; however, for many targets, there are weak correlations. Several physical features can explain the difference between difficult- and easy-to-score targets. For example, strong correlations exist between the score and DockQ for targets with highly intertwined monomers and many interface contacts. We also develop a new score based on only three physical features that matches or exceeds the performance of current PPI scoring functions. These results emphasize that PPI prediction can be improved by focusing on correlations between the PPI score and DockQ and incorporating more discriminating physical features into PPI scoring functions.

Open access
preprint2025arXiv

Particle-scale origin of quadrupolar non-affine displacement fields in granular solids

In this work, we identify the local structural defects that control the non-affine displacement fields in jammed disk packings subjected to athermal, quasistatic (AQS) simple shear. While complex non-affine displacement fields typically occur during simple shear, isolated effective quadrupoles are also observed and their probability increases with increasing pressure. We show that the emergence of an isolated effective quadrupole requires the breaking of an interparticle contact that is aligned with low-frequency, spatially extended vibrational modes. Since the Eshelby inhomogeneity problem gives rise to quadrupolar displacement fields in continuum materials, we reformulate and implement Eshelby's equivalent inclusion method (EIM) for jammed disk packings. Using EIM, we show that we can reconstruct the non-affine displacement fields for jammed disk packings in response to applied shear as a sum of discrete Eshelby-like defects that are caused by mismatches in the local stiffnesses of triangles formed from Delaunay triangulation of the disk centers.

Open access
preprint2022arXiv

Evolving Programmable Computational Metamaterials

Granular metamaterials are a promising choice for the realization of mechanical computing devices. As preliminary evidence of this, we demonstrate here how to embed Boolean logic gates (AND and XOR) into a granular metamaterial by evolving where particular grains are placed in the material. Our results confirm the existence of gradients of increasing "AND-ness" and "XOR-ness" within the space of possible materials that can be followed by evolutionary search. We measure the computational functionality of a material by probing how it transforms bits encoded as vibrations with zero or non-zero amplitude. We compared the evolution of materials built from mass-contrasting particles and materials built from stiffness-contrasting particles, and found that the latter were more evolvable. We believe this work may pave the way toward evolutionary design of increasingly sophisticated, programmable, and computationally dense metamaterials with certain advantages over more traditional computational substrates.

Open access
preprint2022arXiv

The glass-forming ability of binary Lennard-Jones systems

The glass-forming ability (GFA) of alloys, colloidal dispersions, and other particulate materials, as measured by the critical cooling rate RcR_c, can span more than ten orders of magnitude. Even after numerous previous studies, the physical features that control the GFA are still not well understood. For example, it is well-known that mixtures are better glass-formers than monodisperse systems and that particle size and cohesive energy differences among constituents improve the GFA, but it is not currently known how particle size differences couple to cohesive energy differences to determine the GFA. We perform molecular dynamics simulations to determine the GFA of equimolar, binary Lennard-Jones (LJ) mixtures versus the normalized cohesive energy difference ε_ε_\_ and mixing energy εˉAB\bar ε_{AB} between particles A and B. We find several important results. First, the log⁡10Rc\log_{10} R_c contours in the εˉAB\bar ε_{AB}-ε_ε_\_ plane are ellipsoidal for all diameter ratios, and thus RcR_c is determined by the Mahalanobis distance dMd_M from a given point in the εˉAB\bar ε_{AB}-ε_ε_\_ plane to the center of the ellipsoidal contours. Second, LJ systems for which the larger particles have larger cohesive energy are generally better glass formers than those for which the larger particles have smaller cohesive energy. Third, dM(ε_,εˉAB)d_M(ε_\_,\bar ε_{AB}) is determined by the relative Voronoi volume difference between particles and local chemical order SABS_{AB}, which gives the average fraction of nearest-neighbor B particles surrounding an A particle and vice-versa. In particular, the shifted Mahalanobis distance dM−dM0d_M - d^0_M versus the shifted chemical order SAB−SAB0S_{AB}-S_{AB}^0 collapses onto a hyperbolic master curve for all diameter ratios.

Open access
preprint2022arXiv

The Jamming Transition and the Marginally Stable Solid

We review the physics of jamming from the theoretical, experimental and numerical perspectives. We summarize the mean-field theory of jamming and the marginally stable solid phase, with particular emphasis on the connection with the Replica Symmetry Breaking theory of glasses. We report validations of the mean-field theory of jamming from experimental and numerical studies of critical behaviors near the transition. In particular, we describe the physics of jamming of frictionless, spherical particles, as well as more recent work on jamming of frictionless, non-spherical particles and frictional, nearly spherical particles. We also present current efforts in expanding the mean-field theory to systems that more closely resemble externally driven granular media, cell aggregates, and active colloidal suspensions.

Open access
preprint2020arXiv

Shear response of granular packings compressed above jamming onset

We investigate the mechanical response of jammed packings of repulsive, frictionless spherical particles undergoing isotropic compression. Prior simulations of the soft-particle model, where the repulsive interactions scale as a power-law in the interparticle overlap with exponent αα, have found that the ensemble-averaged shear modulus ⟨G⟩\langle G \rangle increases with pressure PP as ∼P(α−3/2)/(α−1)\sim P^{(α-3/2)/(α-1)} at large pressures. However, a deep theoretical understanding of this scaling behavior is lacking. We show that the shear modulus of jammed packings of frictionless, spherical particles has two key contributions: 1) continuous variations as a function of pressure along geometrical families, for which the interparticle contact network does not change, and 2) discontinuous jumps during compression that arise from changes in the contact network. We show that the shear modulus of the first geometrical family for jammed packings can be collapsed onto a master curve: G(1)/G0=(P/P0)(α−2)/(α−1)−P/P0G^{(1)}/G_0 = (P/P_0)^{(α-2)/(α-1)} - P/P_0, where P0∼N−2(α−1)P_0 \sim N^{-2(α-1)} is a characteristic pressure that separates the two power-law scaling regions and G0∼N−2(α−3/2)G_0 \sim N^{-2(α-3/2)}. Deviations from this form can occur when there is significant non-affine particle motion near changes in the contact network. We further show that ⟨G(P)⟩\langle G (P)\rangle is not simply a sum of two power-laws, but ⟨G⟩∼(P/Pc)a\langle G \rangle \sim (P/P_c)^a, where a≈(α−2)/(α−1)a \approx (α-2)/(α-1) in the P→0P \rightarrow 0 limit and ⟨G⟩∼(P/Pc)b\langle G \rangle \sim (P/P_c)^b, where b≳(α−3/2)/(α−1)b \gtrsim (α-3/2)/(α-1) above a characteristic pressure PcP_c. In addition, the magnitudes of both contributions to ⟨G⟩\langle G\rangle from geometrical families and changes in the contact network remain comparable in the large-system limit for P>PcP >P_c.

Open access
preprint2020arXiv

Using physical features of protein core packing to distinguish real proteins from decoys

The ability to consistently distinguish real protein structures from computationally generated model decoys is not yet a solved problem. One route to distinguish real protein structures from decoys is to delineate the important physical features that specify a real protein. For example, it has long been appreciated that the hydrophobic cores of proteins contribute significantly to their stability. As a dataset of decoys to compare with real protein structures, we studied submissions to the bi-annual CASP competition (specifically CASP11, 12, and 13), in which researchers attempt to predict the structure of a protein only knowing its amino acid sequence. Our analysis reveals that many of the submissions possess cores that do not recapitulate the features that define real proteins. In particular, the model structures appear more densely packed (because of energetically unfavorable atomic overlaps), contain too few residues in the core, and have improper distributions of hydrophobic residues throughout the structure. Based on these observations, we developed a deep learning method, which incorporates key physical features of protein cores, to predict how well a computational model recapitulates the real protein structure without knowledge of the structure of the target sequence. By identifying the important features of protein structure, our method is able to rank decoys from the CASP competitions equally well, if not better than, state-of-the-art methods that incorporate many additional features.

Open access
preprint2019arXiv

Analyses of protein cores reveal fundamental differences between solution and crystal structures

There have been several studies suggesting that protein structures solved by NMR spectroscopy and x-ray crystallography show significant differences. To understand the origin of these differences, we assembled a database of high-quality protein structures solved by both methods. We also find significant differences between NMR and crystal structures---in the root-mean-square deviations of the Cα_α atomic positions, identities of core amino acids, backbone and sidechain dihedral angles, and packing fraction of core residues. In contrast to prior studies, we identify the physical basis for these differences by modelling protein cores as jammed packings of amino-acid-shaped particles. We find that we can tune the jammed packing fraction by varying the degree of thermalization used to generate the packings. For an athermal protocol, we find that the average jammed packing fraction is identical to that observed in the cores of protein structures solved by x-ray crystallography. In contrast, highly thermalized packing-generation protocols yield jammed packing fractions that are even higher than those observed in NMR structures. These results indicate that thermalized systems can pack more densely than athermal systems, which suggests a physical basis for the structural differences between protein structures solved by NMR and x-ray crystallography.

Open access
preprint2019arXiv

Intrinsic dissipation mechanisms in metallic glass resonators

Micro- and nano-resonators have important applications including sensing, navigation, and biochemical detection. Their performance is quantified using the quality factor QQ, which gives the ratio of the energy stored to the energy dissipated per cycle. Metallic glasses are a promising materials class for micro- and nano-scale resonators since they are amorphous and can be fabricated precisely into complex shapes on these lengthscales. To understand the intrinsic dissipation mechanisms that ultimately limit large QQ-values in metallic glasses, we perform molecular dynamics simulations to model metallic glass resonators subjected to bending vibrations. We calculate the vibrational density of states, redistribution of energy from the fundamental mode of vibration, and QQ versus the kinetic energy per atom KK of the excitation. In the linear and nonlinear response regimes where there are no atomic rearrangements, we find that Q→∞Q \rightarrow \infty (since we do not consider coupling to the environment). We identify a characteristic KrK_r above which atomic rearrangements occur, and there is significant energy leakage from the fundamental mode to higher frequencies, causing finite QQ. Thus, KrK_r is a critical parameter determining resonator performance. We show that KrK_r decreases as a power-law, Kr∼N−k,K_r\sim N^{-k}, with increasing system size NN, where k≈1.3k \approx 1.3. We estimate the critical strain ⟨γr⟩∼10−8\langle γ_r \rangle \sim 10^{-8} for micron-sized resonators below which atomic rearrangements do not occur, and thus large QQ-values can be obtained when they are operated below γrγ_r. We find that KrK_r for amorphous resonators is comparable to that for resonators with crystalline order.

Open access
preprint2019arXiv

Pressure-dependent shear response of jammed packings of spherical particles

The mechanical response of packings of purely repulsive, spherical particles to athermal, quasistatic simple shear near jamming onset is highly nonlinear. Previous studies have shown that, at small pressure pp, the ensemble-averaged static shear modulus ⟨G−G0⟩\langle G-G_0 \rangle scales with pαp^α, where α≈1α\approx 1, but above a characteristic pressure p∗∗p^{**}, ⟨G−G0⟩∼pβ\langle G-G_0 \rangle \sim p^β, where β≈0.5β\approx 0.5. However, we find that the shear modulus GiG^i for an individual packing typically decreases linearly with pp along a geometrical family where the contact network does not change. We resolve this discrepancy by showing that, while the shear modulus does decrease linearly within geometrical families, ⟨G⟩\langle G \rangle also depends on a contribution from discontinuous jumps in ⟨G⟩\langle G \rangle that occur at the transitions between geometrical families. For p>p∗∗p > p^{**}, geometrical-family and rearrangement contributions to ⟨G⟩\langle G \rangle are of opposite signs and remain comparable for all system sizes. ⟨G⟩\langle G \rangle can be described by a scaling function that smoothly transitions between the two power-law exponents αα and ββ. We also demonstrate the phenomenon of {\it compression unjamming}, where a jammed packing can unjam via isotropic compression.

Open access
preprint2016arXiv

Reconstruction of Ordinary Differential Equations From Time Series Data

We develop a numerical method to reconstruct systems of ordinary differential equations (ODEs) from time series data without {\it a priori} knowledge of the underlying ODEs using sparse basis learning and sparse function reconstruction. We show that employing sparse representations provides more accurate ODE reconstruction compared to least-squares reconstruction techniques for a given amount of time series data. We test and validate the ODE reconstruction method on known 1D, 2D, and 3D systems of ODEs. The 1D system possesses two stable fixed points; the 2D system possesses an oscillatory fixed point with closed orbits; and the 3D system displays chaotic dynamics on a strange attractor. We determine the amount of data required to achieve an error in the reconstructed functions to less than 0.1%0.1\%. For the reconstructed 1D and 2D systems, we are able to match the trajectories from the original ODEs even at long times. For the 3D system with chaotic dynamics, as expected, the trajectories from the original and reconstructed systems do not match at long times, but the reconstructed and original models possess similar Lyapunov exponents. Now that we have validated this ODE reconstruction method on known models, it can be employed in future studies to identify new systems of ODEs using time series data from deterministic systems for which there is no currently known ODE model.

Open access
preprint2015arXiv

Asymmetric crystallization during cooling and heating in model glass-forming systems

We perform molecular dynamics (MD) simulations of the crystallization process in binary Lennard-Jones systems during heating and cooling to investigate atomic-scale crystallization kinetics in glass-forming materials. For the cooling protocol, we prepared equilibrated liquids above the liquidus temperature TlT_l and cooled each sample to zero temperature at rate RcR_c. For the heating protocol, we first cooled equilibrated liquids to zero temperature at rate RpR_p and then heated the samples to temperature T>TlT > T_l at rate RhR_h. We measured the critical heating and cooling rates Rh∗R_h^* and Rc∗R_c^*, below which the systems begin to form a substantial fraction of crystalline clusters during the heating and cooling protocols. We show that Rh∗>Rc∗R_h^* > R_c^*, and that the asymmetry ratio Rh∗/Rc∗R_h^*/R_c^* includes an intrinsic contribution that increases with the glass-forming ability (GFA) of the system and a preparation-rate dependent contribution that increases strongly as Rp→Rc∗R_p \rightarrow R_c^* from above. We also show that the predictions from classical nucleation theory (CNT) can qualitatively describe the dependence of the asymmetry ratio on the GFA and preparation rate RpR_p from the MD simulations and results for the asymmetry ratio measured in Zr- and Au-based bulk metallic glasses (BMG). This work emphasizes the need for and benefits of an improved understanding of crystallization processes in BMGs and other glass-forming systems.

Open access
preprint2015arXiv

Beyond packing of hard spheres: The effects of core softness, non-additivity, intermediate-range repulsion, and many-body interactions on the glass-forming ability of bulk metallic glasses

When a liquid is cooled well below its melting temperature at a rate that exceeds the critical cooling rate RcR_c, the crystalline state is bypassed and an amorphous glassy state forms instead. RcR_c (or the corresponding critical casting thickness dcd_c) characterizes the glass-forming ability (GFA) of each material. While silica is an excellent glass-former with small Rc<10−2R_c<10^{-2} K/s, pure metals and most alloys are poor glass-formers with large Rc>1010R_c>10^{10} K/s. Only in the past thirty years have bulk metallic glasses (BMGs) been identified with RcR_c approaching that for silica. Recent simulations have shown that hard-sphere models are able to identify the atomic size ratio and number fraction regime where BMGs exist with critical cooling rates more than 13 orders of magnitude smaller than those for pure metals. However, there are many other features of interatomic potentials beyond hard-core interactions. How do these other features affect the glass-forming ability of BMGs? We perform molecular dynamics simulations to determine how variations in the softness and non-additivity of the repulsive core and form of the interatomic pair potential at intermediate distances affect the GFA of binary alloys. These variations in the interatomic pair potential allow us to introduce geometric frustration and change the crystal phases that compete with glass formation. We also investigate the effect of tuning the strength of the many-body interactions from zero to the full embedded atom model on the GFA for pure metals. We then employ the full embedded atom model for binary BMGs and show that hard-core interactions play the dominant role in setting the GFA of alloys, while other features of the interatomic potential only change the GFA by one to two orders of magnitude.

Open access
preprint2015arXiv

Onset and cessation of motion in hydrodynamically sheared granular beds

We performed molecular dynamics simulations of granular beds driven by a model hydrodynamic shear flow to elucidate general grain-scale mechanisms that determine the onset and cessation of sediment transport. By varying the Shields number (the nondimensional shear stress at the top of the bed) and particle Reynolds number (the ratio of particle inertia to viscous damping), we explore how variations of the fluid flow rate, particle inertia, and fluid viscosity affect the onset and cessation of bed motion. For low to moderate particle Reynolds numbers, a critical boundary separates mobile and static states. Transition times between these states diverge as this boundary is approached both from above and below. At high particle Reynolds number, inertial effects become dominant, and particle motion can be sustained well below flow rates at which mobilization of a static bed occurs. We also find that the onset of bed motion (for both low and high particle Reynolds numbers) is described by Weibullian weakest-link statistics, and thus is crucially dependent on the packing structure of the granular bed, even deep beneath the surface.

Open access
preprint2015arXiv

Outcome prediction in mathematical models of immune response to infection

Clinicians need to predict patient outcomes with high accuracy as early as possible after disease inception. In this manuscript, we show that patient-to-patient variability sets a fundamental limit on outcome prediction accuracy for a general class of mathematical models for the immune response to infection. However, accuracy can be increased at the expense of delayed prognosis. We investigate several systems of ordinary differential equations (ODEs) that model the host immune response to a pathogen load. Advantages of systems of ODEs for investigating the immune response to infection include the ability to collect data on large numbers of `virtual patients', each with a given set of model parameters, and obtain many time points during the course of the infection. We implement patient-to-patient variability vv in the ODE models by randomly selecting the model parameters from Gaussian distributions with variance vv that are centered on physiological values. We use logistic regression with one-versus-all classification to predict the discrete steady-state outcomes of the system. We find that the prediction algorithm achieves near 100%100\% accuracy for v=0v=0, and the accuracy decreases with increasing vv for all ODE models studied. The fact that multiple steady-state outcomes can be obtained for a given initial condition, i.e. the basins of attraction overlap in the space of initial conditions, limits the prediction accuracy for v>0v>0. Increasing the elapsed time of the variables used to train and test the classifier, increases the prediction accuracy, while adding explicit external noise to the ODE models decreases the prediction accuracy. Our results quantify the competition between early prognosis and high prediction accuracy that is frequently encountered by clinicians.

Open access
preprint2015arXiv

Random close packing in protein cores

Shortly after the determination of the first protein x-ray crystal structures, researchers analyzed their cores and reported packing fractions ϕ≈0.75ϕ\approx 0.75, a value that is similar to close packing equal-sized spheres. A limitation of these analyses was the use of `extended atom' models, rather than the more physically accurate `explicit hydrogen' model. The validity of using the explicit hydrogen model is proved by its ability to predict the side chain dihedral angle distributions observed in proteins. We employ the explicit hydrogen model to calculate the packing fraction of the cores of over 200200 high resolution protein structures. We find that these protein cores have ϕ≈0.55ϕ\approx 0.55, which is comparable to random close-packing of non-spherical particles. This result provides a deeper understanding of the physical basis of protein structure that will enable predictions of the effects of amino acid mutations and design of new functional proteins.

Open access
preprint2015arXiv

Shear accelerated crystallization in a supercooled atomic liquid

A bulk metallic glass forming alloy is subjected to shear flow in its supercooled state by compression of a short rod to produce a flat disc. The resulting material exhibits enhanced crystallization kinetics during isothermal annealing as reflected in the decrease of the crystallization time relative to the non-deformed case. The transition from quiescent to shear-accelerated crystallization is linked to strain accumulated during shear flow above a critical shear rate γ˙c≈0.3\dotγ_c\approx 0.3 s−1^{-1} which corresponds to Péclet number, Pe∼O(1)Pe\sim\mathcal{O}(1). The observation of shear accelerated crystallization in an atomic system at modest shear rates is uncommon. It is made possible here by the substantial viscosity of the supercooled liquid which increases strongly with temperature in the approach to the glass transition. We may therefore anticipate the encounter of non-trivial shear-related effects during thermoplastic deformation of similar systems.

Open access
preprint2014arXiv

Calibrated Langevin dynamics simulations of intrinsically disordered proteins

We perform extensive coarse-grained (CG) Langevin dynamics simulations of intrinsically disordered proteins (IDPs), which possess fluctuating conformational statistics between that for excluded volume random walks and collapsed globules. Our CG model includes repulsive steric, attractive hydrophobic, and electrostatic interactions between residues and is calibrated to a large collection of single-molecule fluorescence resonance energy transfer data on the inter-residue separations for 36 pairs of residues in five IDPs: αα-, ββ-, and γγ-synuclein, the microtubule-associated protein ττ, and prothymosin αα. We find that our CG model is able to recapitulate the average inter-residue separations regardless of the choice of the hydrophobicity scale, which shows that our calibrated model can robustly capture the conformational dynamics of IDPs. We then employ our model to study the scaling of the radius of gyration with chemical distance in 11 known IDPs. We identify a strong correlation between the distance to the dividing line between folded proteins and IDPs in the mean charge and hydrophobicity space and the scaling exponent of the radius of gyration with chemical distance along the protein.

Open access
preprint2014arXiv

Can the packing efficiency of binary hard spheres explain the glass-forming ability of bulk metallic glasses?

We perform molecular dynamics simulations to compress binary hard spheres into jammed packings as a function of the compression rate RR, size ratio αα, and number fraction xSx_S of small particles to determine the connection between the glass-forming ability (GFA) and packing efficiency in bulk metallic glasses (BMGs). We define the GFA by measuring the critical compression rate RcR_c, below which jammed hard-sphere packings begin to form "random crystal" structures with defects. We find that for systems with α≳0.8α\gtrsim 0.8 that do not de-mix, RcR_c decreases strongly with ΔϕJΔϕ_J, as Rc∼exp⁡(−1/ΔϕJ2)R_c \sim \exp(-1/Δϕ_J^2), where ΔϕJΔϕ_J is the difference between the average packing fraction of the amorphous packings and random crystal structures at RcR_c. Systems with α≲0.8α\lesssim 0.8 partially de-mix, which promotes crystallization, but we still find a strong correlation between RcR_c and ΔϕJΔϕ_J. We show that known metal-metal BMGs occur in the regions of the αα and xSx_S parameter space with the lowest values of RcR_c for binary hard spheres. Our results emphasize that maximizing GFA in binary systems involves two competing effects: minimizing αα to increase packing efficiency, while maximizing αα to prevent de-mixing.

Open access
preprint2014arXiv

Quasi-One Dimensional Models for Glassy Dynamics

We describe numerical simulations and analyses of a quasi-one-dimensional (Q1D) model of glassy dynamics. In this model, hard rods undergo Brownian dynamics through a series of narrow channels connected by JJ intersections. We do not allow the rods to turn at the intersections, and thus there is a single, continuous route through the system. This Q1D model displays caging behavior, collective particle rearrangements, and rapid growth of the structural relaxation time, which are also found in supercooled liquids and glasses. The mean-square displacement Σ(t)Σ(t) for this Q1D model displays several dynamical regimes: 1) short-time diffusion Σ(t)∼tΣ(t) \sim t, 2) a plateau in the mean-square displacement caused by caging behavior, 3) single-file diffusion characterized by anomalous scaling Σ(t)∼t0.5Σ(t) \sim t^{0.5} at intermediate times, and 4) a crossover to long-time diffusion Σ(t)∼tΣ(t) \sim t for times tt that grow with the complexity of the circuit. We develop a general procedure for determining the structural relaxation time tDt_D, beyond which the system undergoes long-time diffusion, as a function of the packing fraction ϕϕ and system topology. This procedure involves several steps: 1) define a set of distinct microstates in configuration space of the system, 2) construct a directed network of microstates and transitions between them, 3) identify minimal, closed loops in the network that give rise to structural relaxation, 4) determine the frequencies of `bottleneck' microstates that control the slow dynamics and time required to transition out of them, and 5) use the microstate frequencies and lifetimes to deduce tD(ϕ)t_D(ϕ). We find that tDt_D obeys power-law scaling, tD∼(ϕ∗−ϕ)−αt_D\sim (ϕ^* - ϕ)^{-α}, where both ϕ∗ϕ^* (signaling complete kinetic arrest) and α>0α>0 depend on the system topology.

Open access
preprint2014arXiv

The glass-forming ability of model metal-metalloid alloys

Bulk metallic glasses (BMGs) are amorphous alloys with desirable mechanical properties and processing capabilities. To date, the design of new BMGs has largely employed empirical rules and trial-and-error experimental approaches. Ab initio computational methods are currently prohibitively slow to be practically used in searching the vast space of possible atomic combinations for bulk glass formers. Here, we perform molecular dynamics simulations of a coarse-grained, anisotropic potential, which mimics interatomic covalent bonding, to measure the critical cooling rates for metal-metalloid alloys as a function of the atomic size ratio σS/σLσ_S/σ_L and number fraction xSx_S of the metalloid species. We show that the regime in the space of σS/σLσ_S/σ_L and xSx_S where well-mixed, optimal glass formers occur for patchy and LJ particle mixtures coincides with that for experimentally observed metal-metalloid glass formers. Our simple computational model provides the capability to perform combinatorial searches to identify novel glass-forming alloys.

Open access
preprint2014arXiv

The statistics of frictional families

We develop a theoretical description for mechanically stable frictional packings in terms of the difference between the total number of contacts required for isostatic packings of frictionless disks and the number of contacts in frictional packings, m=Nc0−Ncm=N_c^0-N_c. The saddle order mm represents the number of unconstrained degrees of freedom that a static packing would possess if friction were removed. Using a novel numerical method that allows us to enumerate disk packings for each mm, we show that the probability to obtain a packing with saddle order mm at a given static friction coefficient μμ, Pm(μ)P_m(μ), can be expressed as a power-series in μμ. Using this form for Pm(μ)P_m(μ), we quantitatively describe the dependence of the average contact number on friction coefficient for static disk packings obtained from direct simulations of the Cundall-Strack model for all μμ and NN.

Open access
preprint2013arXiv

Computational studies of the glass-forming ability of model bulk metallic glasses

Bulk metallic glasses (BMGs) are produced by rapidly thermally quenching supercooled liquid metal alloys below the glass transition temperature at rates much faster than the critical cooling rate R_c below which crystallization occurs. The glass-forming ability of BMGs increases with decreasing R_c, and thus good glass-formers possess small values of R_c. We perform molecular dynamics simulations of binary Lennard-Jones (LJ) mixtures to quantify how key parameters, such as the stoichiometry, particle size difference, attraction strength, and heat of mixing, influence the glass-formability of model BMGs. For binary LJ mixtures, we find that the best glass-forming mixtures possess atomic size ratios (small to large) less than 0.92 and stoichiometries near 50:50 by number. In addition, weaker attractive interactions between the smaller atoms facilitate glass formation, whereas negative heats of mixing (in the experimentally relevant regime) do not change R_c significantly. These studies represent a first step in the development of computational methods for quantitatively predicting glass-formability.

Open access
preprint2013arXiv

Isostaticity at Frictional Jamming

Amorphous packings of frictionless, spherical particles are isostatic at jamming onset, with the number of constraints (contacts) equal to the number of degrees of freedom. Their structural and mechanical properties are controlled by the interparticle contact network. In contrast, amorphous packings of frictional particles are typically hyperstatic at jamming onset. We perform extensive numerical simulations in two dimensions of the geometrical asperity (GA) model for static friction, to further investigate the role of isostaticity. In the GA model, interparticle forces are obtained by summing up purely repulsive central forces between periodically spaced circular asperities on contacting grains. We compare the packing fraction, contact number, mobilization distribution, and vibrational density of states using the GA model to those generated using the Cundall-Strack (CS) approach. We find that static packings of frictional disks obtained from the GA model are mechanically stable and isostatic when we consider interactions between asperities on contacting particles. The crossover in the structural and mechanical properties of static packings from frictionless to frictional behavior as a function of the static friction coefficient coincides with a change in the type of interparticle contacts and the disappearance of a peak in the density of vibrational modes for the GA model. These results emphasize that mesoscale features of the model for static friction play an important role in determining the properties of granular packings.

Open access
preprint2013arXiv

Minimal energy packings of nearly flexible polymers

We extend recent studies of the minimal energy packings of short flexible polymers with hard-core-like repulsions and short-range attractions to include bond-angle interactions with the aim of describing the collapsed conformations of `colloidal' polymers. We find that flexible tangent sticky-hard-sphere (t-SHS) packings provide a useful perturbative basis for analyzing polymer packings with nonzero bending stiffness only for {\it small} ratios of the stiffnesses for the bond-angle (kbk_b) and pair (kck_c) interactions, i.e. kbcrit/kc≲0.01k_b^{\rm crit}/k_c \lesssim 0.01 for N<10N<10 monomers, and the critical ratio decreases with NN. Below kbcritk_b^{crit}, angular interactions give rise to an exponential (in NN) increase in the number of distinct angular energies arising from the diversity of covalent backbone paths through t-SHS packings. As kbk_b increases above kbcritk_b^{crit}, the low-lying energy landscape changes dramatically as finite bending stiffness alters the structure of the polymer packings. This study lays the groundwork for exact-enumeration studies of the collapsed states of t-SHS-like models with larger bending stiffness.

Open access
preprint2013arXiv

Particle-scale reversibility in athermal particulate media below jamming

We perform numerical simulations of athermal repulsive frictionless disks and spheres in two and three spatial dimensions undergoing cyclic quasi-static simple shear to investigate particle-scale reversible motion. We identify three classes of steady-state dynamics as a function of packing fraction ϕand maximum strain amplitude per cycle γ_{\rm max}. Point-reversible states, where particles do not collide and exactly retrace their intra-cycle trajectories, occur at low ϕand γ_{\rm max}. Particles in loop-reversible states undergo numerous collisions and execute complex trajectories, but return to their initial positions at the end of each cycle. Loop-reversible dynamics represents a novel form of self-organization that enables reliable preparation of configurations with specified structural and mechanical properties over a broad range of ϕfrom contact percolation to jamming onset at ϕ_J. For sufficiently large ϕand γ_{\rm max}, systems display irreversible dynamics with nonzero self-diffusion.

Open access
preprint2013arXiv

Vibrations in jammed solids: Beyond linear response

We propose a `phase diagram' for particulate systems that interact via purely repulsive contact forces, such as granular media and colloidal suspensions. We identify and characterize two distinct classes of behavior as a function of the input kinetic energy per degree of freedom T0T_0 and packing fraction deviation above and below jamming onset Δϕ=ϕ−ϕJΔϕ=ϕ- ϕ_J using numerical simulations of purely repulsive frictionless disks. Iso-coordinated solids (ICS) only occur above jamming for Δϕ>Δϕc(T0)Δϕ> Δϕ_c(T_0); they possess average coordination number equal to the isostatic value (<z>=ziso< z> = z_{\rm iso}) required for mechanically stable packings. ICS display harmonic vibrational response, where the density of vibrational modes from the Fourier transform of the velocity autocorrelation function is a set of sharp peaks at eigenfrequencies ωkdω_k^d of the dynamical matrix evaluated at T0=0T_0=0. Hypo-coordinated solids (HCS) occur both above and below jamming onset within the region defined by Δϕ>Δϕ−∗(T0)Δϕ> Δϕ^*_-(T_0), Δϕ<Δϕ+∗(T0)Δϕ< Δϕ^*_+(T_0), and Δϕ>Δϕcb(T0)Δϕ> Δϕ_{cb}(T_0). In this region, the network of interparticle contacts fluctuates with <z>≈ziso/2< z> \approx z_{\rm iso}/2, but cage-breaking particle rearrangements do not occur. The HCS vibrational response is nonharmonic, {\it i.e} the density of vibrational modes D(ω)D(ω) is not a collection of sharp peaks at ωkdω_k^d, and its precise form depends on the measurement method. For Δϕ>Δϕcb(T0)Δϕ> Δϕ_{cb}(T_0) and Δϕ<Δϕ−∗(T0)Δϕ< Δϕ^*_{-}(T_0), the system behaves as a hard-particle liquid.

Open access
preprint2013arXiv

Vibrations of Jammed Disk Packings with Hertzian Interactions

Contact breaking and Hertzian interactions between grains can both give rise to nonlinear vibrational response of static granular packings. We perform molecular dynamics simulations at constant energy in 2D of frictionless bidisperse disks that interact via Hertzian spring potentials as a function of energy and measure directly the vibrational response from the Fourier transform of the velocity autocorrelation function. We compare the measured vibrational response of static packings near jamming onset to that obtained from the eigenvalues of the dynamical matrix to determine the temperature above which the linear response breaks down. We compare packings that interact via single-sided (purely repulsive) and double-sided Hertzian spring interactions to disentangle the effects of the shape of the potential from contact breaking. Our studies show that while Hertzian interactions lead to weak nonlinearities in the vibrational behavior (e.g. the generation of harmonics of the eigenfrequencies of the dynamical matrix), the vibrational response of static packings with Hertzian contact interactions is dominated by contact breaking as found for systems with repulsive linear spring interactions.

Open access
preprint2012arXiv

Constraints and vibrations in static packings of ellipsoidal particles

We numerically investigate the mechanical properties of static packings of ellipsoidal particles in 2D and 3D over a range of aspect ratio and compression ΔϕΔϕ. While amorphous packings of spherical particles at jamming onset (Δϕ=0Δϕ=0) are isostatic and possess the minimum contact number zisoz_{\rm iso} required for them to be collectively jammed, amorphous packings of ellipsoidal particles generally possess fewer contacts than expected for collective jamming (z<zisoz < z_{\rm iso}) from naive counting arguments, which assume that all contacts give rise to linearly independent constraints on interparticle separations. To understand this behavior, we decompose the dynamical matrix M=H−SM=H-S for static packings of ellipsoidal particles into two important components: the stiffness HH and stress SS matrices. We find that the stiffness matrix possesses N(ziso−z)N(z_{\rm iso} - z) eigenmodes e^0{\hat e}_0 with zero eigenvalues even at finite compression, where NN is the number of particles. In addition, these modes e^0{\hat e}_0 are nearly eigenvectors of the dynamical matrix with eigenvalues that scale as ΔϕΔϕ, and thus finite compression stabilizes packings of ellipsoidal particles. At jamming onset, the harmonic response of static packings of ellipsoidal particles vanishes, and the total potential energy scales as δ4δ^4 for perturbations by amplitude δδ along these `quartic' modes, e^0{\hat e}_0. These findings illustrate the significant differences between static packings of spherical and ellipsoidal particles.

Open access
preprint2012arXiv

Molecular Simulations of the Fluctuating Conformational Dynamics of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) do not possess well-defined three-dimensional structures in solution under physiological conditions. We develop all-atom, united-atom, and coarse-grained Langevin dynamics simulations for the IDP alpha-synuclein that include geometric, attractive hydrophobic, and screened electrostatic interactions and are calibrated to the inter-residue separations measured in recent smFRET experiments. We find that alpha-synuclein is disordered with conformational statistics that are intermediate between random walk and collapsed globule behavior. An advantage of calibrated molecular simulations over constraint methods is that physical forces act on all residues, not only on residue pairs that are monitored experimentally, and these simulations can be used to study oligomerization and aggregation of multiple alpha-synuclein proteins that may precede amyloid formation.

Open access
preprint2012arXiv

Rods are less fragile than spheres: Structural relaxation in dense liquids composed of anisotropic particles

We perform extensive molecular dynamics simulations of dense liquids composed of bidisperse dimer- and ellipse-shaped particles in 2D that interact via repulsive contact forces. We measure the structural relaxation times obtained from the long-time decay of the self-part of the intermediate scattering function for the translational and rotational degrees of freedom (DOF) as a function of packing fraction ϕ, temperature T, and aspect ratio α. We are able to collapse the ϕand T-dependent structural relaxation times for disks, and dimers and ellipses over a wide range of α, onto a universal scaling function {\cal F}_{\pm}(|ϕ-ϕ_0|,T,α), which is similar to that employed in previous studies of dense liquids composed of purely repulsive spherical particles in 3D. {\cal F_{\pm}} for both the translational and rotational DOF are characterized by the α-dependent scaling exponents μand δand packing fraction ϕ_0(α) that signals the crossover in the scaling form {\cal F}_{\pm} from hard-particle dynamics to super-Arrhenius behavior for each aspect ratio. We find that the fragility at ϕ_0, m(ϕ_0), decreases monotonically with increasing aspect ratio for both ellipses and dimers. Moreover, the results for the slow dynamics of dense liquids composed of dimer- and ellipse-shaped particles are qualitatively the same, despite the fact that zero-temperature static packings of dimers are isostatic, while static packings of ellipses are hypostatic.

Open access
preprint2012arXiv

Structure of finite sphere packings via exact enumeration: Implications for colloidal crystal nucleation

We analyze the geometric structure and mechanical stability of a complete set of isostatic and hyperstatic sphere packings obtained via exact enumeration. The number of nonisomorphic isostatic packings grows exponentially with the number of spheres NN, and their diversity of structure and symmetry increases with increasing NN and decreases with increasing hyperstaticity H≡Nc−NISOH \equiv N_c - N_{ISO}, where NcN_c is the number of pair contacts and NISO=3N−6N_{ISO} = 3N-6. Maximally contacting packings are in general neither the densest nor the most symmetric. Analyses of local structure show that the fraction ff of nuclei with order compatible with the bulk (RHCP) crystal decreases sharply with increasing NN due to a high propensity for stacking faults, 5- and near-5-fold symmetric structures, and other motifs that preclude RHCP order. While ff increases with increasing HH, a significant fraction of hyperstatic nuclei for NN as small as 11 retain non-RHCP structure. Classical theories of nucleation that consider only spherical nuclei, or only nuclei with the same ordering as the bulk crystal, cannot capture such effects. Our results provide an explanation for the failure of classical nucleation theory for hard-sphere systems of N≲10N\lesssim 10 particles; we argue that in this size regime, it is essential to consider nuclei of unconstrained geometry. Our results are also applicable to understanding kinetic arrest and jamming in systems that interact via hard-core-like repulsive and short-ranged attractive interactions.

Open access
preprint2011arXiv

Glassy dynamics of crystallite formation: The role of covalent bonds

We examine nonequilibrium features of collapse behavior in model polymers with competing crystallization and glass transitions using extensive molecular dynamics simulations. By comparing to "colloidal" systems with no covalent bonds but the same non-bonded interactions, we find three principal results: (i) Tangent-sphere polymers and colloids, in the equilibrium-crystallite phase, have nearly identical static properties when the temperature T is scaled by the crystallization temperature T_{cryst}; (ii) Qualitative features of nonequilibrium relaxation below T_{cryst}, measured by the evolution of local structural properties (such as the number of contacts) toward equilibrium crystallites, are the same for polymers and colloids; and (iii) Significant quantitative differences in rearrangements in polymeric and colloidal crystallites, in both far-from equilibrium and near-equilibrium systems, can be understood in terms of chain connectivity. These results have important implications for understanding slow relaxation processes in collapsed polymers, partially folded, misfolded, and intrinsically disordered proteins.

Open access
preprint2011arXiv

Photonic Band Gaps in 3D Network Structures with Short-range Order

We present a systematic study of photonic band gaps (PBGs) in three-dimensional (3D) photonic amorphous structures (PAS) with short-range order. From calculations of the density of optical states (DOS) for PAS with different topologies, we find that tetrahedrally connected dielectric networks produce the largest isotropic PBGs. Local uniformity and tetrahedral order are essential to the formation of PBGs in PAS, in addition to short-range geometric order. This work demonstrates that it is possible to create broad, isotropic PBGs for vector light fields in 3D PAS without long-range order.

Open access
preprint2011arXiv

The contact percolation transition

Typical quasistatic compression algorithms for generating jammed packings of athermal, purely repulsive particles begin with dilute configurations and then apply successive compressions with relaxation of the elastic energy allowed between each compression step. It is well-known that during isotropic compression athermal systems with purely repulsive interactions undergo a jamming transition at packing fraction ϕJϕ_J from an unjammed state with zero pressure to a jammed, rigid state with nonzero pressure. Using extensive computer simulations, we show that a novel second-order-like transition, the contact percolation transition, which signals the formation of a system-spanning cluster of mutually contacting particles, occurs at ϕP<ϕJϕ_P < ϕ_J, preceding the jamming transition. By measuring the number of non-floppy modes of the dynamical matrix, and the displacement field and time-dependent pressure following compression, we find that the contact percolation transition also heralds the onset of complex spatiotemporal response to applied stress. Thus, highly heterogeneous, cooperative, and non-affine particle motion occurs in unjammed systems significantly below the jamming transition for ϕP<ϕ<ϕJϕ_P < ϕ< ϕ_J, not only for jammed systems with ϕ>ϕJϕ> ϕ_J.

Open access
preprint2010arXiv

A comparison of jamming behavior in systems composed of dimer- and ellipse-shaped particles

We compare the structural and mechanical properties of static packings composed of frictionless convex (ellipses) and concave (rigid dimers) particles in two dimensions. We employ numerical simulations to generate static packings and measure the shear stress in response to applied simple shear strain as a function of the aspect ratio and amount of compression. We find that the behavior near jamming is significantly different for ellipses and dimers even though both shapes are roughly characterized by the aspect ratio and possess the same number of translational and rotational degrees of freedom per particle. For example, we find that ellipse packings are hypostatic (not isostatic as found for dimers), display novel power-law scaling of the static linear shear modulus and contact number with the amount of compression, and possess stress-strain relations that are qualitatively different from that for dimers. Thus, we observe that important macroscopic properties of static packings of anisotropic particles can depend on the microscale geometrical features of individual particles.

Open access
preprint2010arXiv

Jammed particulate systems are inherently nonharmonic

Jammed particulate systems, such as granular media, colloids, and foams, interact via one-sided forces that are nonzero only when particles overlap. We find that systems with one-sided repulsive interactions possess no linear response regime in the large system limit (N→∞N\rightarrow \infty) for all pressures pp (or compressions ΔϕΔϕ), and for all NN near jamming onset p→0p\rightarrow 0. We perform simulations on 2D frictionless bidisperse mechanically stable disk packings over a range of packing fractions Δϕ=ϕ−ϕJΔϕ= ϕ-ϕ_J above jamming onset ϕJϕ_J. We apply perturbations with amplitude δδ to the packings along each eigen-direction from the dynamical matrix and determine whether the response of the system evolving at constant energy remains in the original eigenmode of the perturbation. For δ>δcδ> δ_c, which we calculate analytically, a single contact breaks and fluctuations abruptly spread to all harmonic modes. As δδ increases further all discrete harmonic modes disappear into a continuous frequency band. We find that <δc>∼Δϕ/Nλ<δ_c >\sim Δϕ/N^λ, where 1>λ>0.51 > λ> 0.5, and thus jammed particulate systems are inherently nonharmonic with no linear vibrational response regime as N→∞N\rightarrow \infty over the full range of ΔϕΔϕ, and as Δϕ→0Δϕ\rightarrow 0 at any NN.

Open access
preprint2010arXiv

Minimal energy packings and collapse of sticky tangent hard-sphere polymers

We enumerate all minimal energy packings (MEPs) for small single linear and ring polymers composed of spherical monomers with contact attractions and hard-core repulsions, and compare them to corresponding results for monomer packings. We define and identify ``dividing surfaces" in polymer packings, which reduce the number of arrangements that satisfy hard-sphere and covalent bond constraints. Compared to monomer MEPs, polymer MEPs favor intermediate structural symmetry over high and low symmetries. We also examine the packing-preparation dependence for longer single chains using molecular dynamics simulations. For slow temperature quenches, chains form crystallites with close-packed cores. As quench rate increases, the core size decreases and the exterior becomes more disordered. By examining the contact number, we connect suppression of crystallization to the onset of isostaticity in disordered packings. These studies represent a significant step forward in our ability to predict how the structural and mechanical properties of compact polymers depend on collapse dynamics.

Open access
preprint2010arXiv

Photonic Band Gap Effects in Two-dimensional Polycrystalline and Amorphous Structures

We study numerically the density of optical states (DOS) in two-dimensional photonic structures with short-range positional order, and observe a clear transition from polycrystalline to amorphous photonic systems. In polycrystals, photonic band gaps (PBGs) are formed within individual do- mains, which leads to a depletion of the DOS similar to that in periodic structures. In amorphous photonic media, the domain sizes are too small to form PBGs, thus the depletion of the DOS is weakened significantly. The critical domain size that separates the polycrystalline and amorphous regimes is determined by the attenuation length of Bragg scattering, which depends not only on the degree of positional order but also the refractive index contrast of the photonic material. Even with relatively low refractive index contrast, we find that modest short-range positional order in photonic structures enhances light confinement via collective scattering and interference.

Open access
preprint2010arXiv

Tuning Jammed Frictionless Disk Packings from Isostatic to Hyperstatic

We perform extensive computational studies of two-dimensional static bidisperse disk packings using two distinct packing-generation protocols. The first involves thermally quenching equilibrated liquid configurations to zero temperature over a range of thermal quench rates rr and initial packing fractions followed by compression and decompression in small steps to reach packing fractions ϕJϕ_J at jamming onset. For the second, we seed the system with initial configurations that promote micro- and macrophase-separated packings followed by compression and decompression to ϕJϕ_J. We find that amorphous, isostatic packings exist over a finite range of packing fractions from ϕmin≤ϕJ≤ϕmaxϕ_{\rm min} \le ϕ_J \le ϕ_{\rm max} in the large-system limit, with ϕmax≈0.853ϕ_{\rm max} \approx 0.853. In agreement with previous calculations, we obtain ϕmin≈0.84ϕ_{\rm min} \approx 0.84 for r>r∗r > r^*, where r∗r^* is the rate above which ϕJϕ_J is insensitive to rate. We further compare the structural and mechanical properties of isostatic versus hyperstatic packings. The structural characterizations include the contact number, bond orientational order, and mixing ratios of the large and small particles. We find that the isostatic packings are positionally and compositionally disordered, whereas bond-orientational and compositional order increase with contact number for hyperstatic packings. In addition, we calculate the static shear modulus and normal mode frequencies of the static packings to understand the extent to which the mechanical properties of amorphous, isostatic packings are different from partially ordered packings. We find that the mechanical properties of the packings change continuously as the contact number increases from isostatic to hyperstatic.

Open access
preprint2010arXiv

Viscoplasticity and large-scale chain relaxation in glassy-polymeric strain hardening

A simple theory for glassy polymeric mechanical response which accounts for large scale chain relaxation is presented. It captures the crossover from perfect-plastic response to strong strain hardening as the degree of polymerization NN increases, without invoking entanglements. By relating hardening to interactions on the scale of monomers and chain segments, we correctly predict its magnitude. Strain activated relaxation arising from the need to maintain constant chain contour length reduces the NN dependence of the characteristic relaxation time by a factor ∼ε˙N\sim \dotεN during active deformation at strain rate ε˙\dotε. This prediction is consistent with results from recent experiments and simulations, and we suggest how it may be further tested experimentally.

Open access
preprint2009arXiv

Entropic origin of stress correlations in granular materials

We study the response of granular materials to external stress using experiment, simulation, and theory. We derive an entropic, Ginzburg-Landau functional that enforces mechanical stability and positivity of contact forces. In this framework, the elastic moduli depend only on the applied stress. A combination of this feature and the positivity constraint leads to stress correlations whose shape and magnitude are extremely sensitive to the applied stress. The predictions from the theory describe the stress correlations for both simulations and experiments semiquantitatively.

Open access
preprint2009arXiv

Experimental Demonstration of Nonuniform Frequency Distributions of Granular Packings

We developed a novel experimental technique to generate mechanically stable (MS) packings of frictionless granular disks. We performed a series of coordinated experiments and numerical simulations to enumerate the MS packings in small 2D systems composed of bidisperse disks. We find that frictionless MS packings occur as discrete, well-separated points in configuration space and obtain excellent quantitative agreement between MS packings generated in experiments and simulations. In addition, we observe that MS packing probabilities can vary by many orders of magnitude and are robust with respect to the packing-generation procedure. These results suggest that the most frequent MS packings may dominate the structural and mechanical properties of granular systems. We argue that these results for small systems represent a crucial first-step in constructing a statistical description for large granular systems from the `bottom-up'.

Open access
preprint2009arXiv

Geometrical families of mechanically stable granular packings

We enumerate and classify nearly all of the possible mechanically stable (MS) packings of bidipserse mixtures of frictionless disks in small sheared systems. We find that MS packings form continuous geometrical families, where each family is defined by its particular network of particle contacts. We also monitor the dynamics of MS packings along geometrical families by applying quasistatic simple shear strain at zero pressure. For small numbers of particles (N < 16), we find that the dynamics is deterministic and highly contracting. That is, if the system is initialized in a MS packing at a given shear strain, it will quickly lock into a periodic orbit at subsequent shear strain, and therefore sample only a very small fraction of the possible MS packings in steady state. In studies with N>16, we observe an increase in the period and random splittings of the trajectories caused by bifurcations in configuration space. We argue that the ratio of the splitting and contraction rates in large systems will determine the distribution of MS-packing geometrical families visited in steady-state. This work is part of our long-term research program to develop a master-equation formalism to describe macroscopic slowly driven granular systems in terms of collections of small subsystems.

Open access
preprint2009arXiv

Protein folding on rugged energy landscapes: Conformational diffusion on fractal networks

We employ simulations of model proteins to study folding on rugged energy landscapes. We construct ``first-passage'' networks as the system transitions from unfolded to native states. The nodes and bonds in these networks correspond to basins and transitions between them in the energy landscape. We find power-laws between the folding time and number of nodes and bonds. We show that these scalings are determined by the fractal properties of first-passage networks. Reliable folding is possible in systems with rugged energy landscapes because first passage networks have small fractal dimension.

Open access
preprint2008arXiv

Jamming in Systems Composed of Frictionless Ellipse-Shaped Particles

We study the structural and mechanical properties of jammed ellipse packings, and find that the nature of the jamming transition in these systems is fundamentally different from that for spherical particles. Ellipse packings are generically hypostatic with more degrees of freedom than constraints. The spectra of low energy excitations possess two gaps and three distinct branches over a range of aspect ratios. In the zero compression limit, the energy of the modes in the lowest branch increases {\it quartically} with deformation amplitude, and the density of states possesses a δδ-function at zero frequency. We identify scaling relations that collapse the low-frequency part of the spectra for different aspect ratios. Finally, we find that the degree of hypostaticity is determined by the number of quartic modes of the packing.

Open access
preprint1995arXiv

Lyapunov spectral analysis of a nonequilibrium Ising-like transition

By simulating a nonequilibrium coupled map lattice that undergoes an Ising-like phase transition, we show that the Lyapunov spectrum and related dynamical quantities such as the dimension correlation length~ξδξ_δ are insensitive to the onset of long-range ferromagnetic order. As a function of lattice coupling constant~gg and for certain lattice maps, the Lyapunov dimension density and other dynamical order parameters go through a minimum. The occurrence of this minimum as a function of~gg depends on the number of nearest neighbors of a lattice point but not on the lattice symmetry, on the lattice dimensionality or on the position of the Ising-like transition. In one-space dimension, the spatial correlation length associated with magnitude fluctuations and the length~ξδξ_δ are approximately equal, with both varying linearly with the radius of the lattice coupling.

Open access

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.