BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Mark D. Shattuck

Mark D. Shattuck contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.softcond-mat.mtrl-sciBiomoleculescond-mat.dis-nnEmerging TechnologiesNeural and Evolutionary Computing

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

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.

0 citations0 reviewsNo code linkedOpen 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.

0 citations0 reviewsNo code linkedOpen 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.

0 citations0 reviewsNo code linkedOpen 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 $R_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 $\bar ε_{AB}$ between particles A and B. We find several important results. First, the $\log_{10} R_c$ contours in the $\bar ε_{AB}$-$ε_\_$ plane are ellipsoidal for all diameter ratios, and thus $R_c$ is determined by the Mahalanobis distance $d_M$ from a given point in the $\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, $d_M(ε_\_,\bar ε_{AB})$ is determined by the relative Voronoi volume difference between particles and local chemical order $S_{AB}$, which gives the average fraction of nearest-neighbor B particles surrounding an A particle and vice-versa. In particular, the shifted Mahalanobis distance $d_M - d^0_M$ versus the shifted chemical order $S_{AB}-S_{AB}^0$ collapses onto a hyperbolic master curve for all diameter ratios.

0 citations0 reviewsNo code linkedOpen 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 $\langle G \rangle$ increases with pressure $P$ as $\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)}/G_0 = (P/P_0)^{(α-2)/(α-1)} - P/P_0$, where $P_0 \sim N^{-2(α-1)}$ is a characteristic pressure that separates the two power-law scaling regions and $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 $\langle G (P)\rangle$ is not simply a sum of two power-laws, but $\langle G \rangle \sim (P/P_c)^a$, where $a \approx (α-2)/(α-1)$ in the $P \rightarrow 0$ limit and $\langle G \rangle \sim (P/P_c)^b$, where $b \gtrsim (α-3/2)/(α-1)$ above a characteristic pressure $P_c$. In addition, the magnitudes of both contributions to $\langle G\rangle$ from geometrical families and changes in the contact network remain comparable in the large-system limit for $P >P_c$.

0 citations0 reviewsNo code linkedOpen 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 $Q$, 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 $Q$-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 $Q$ versus the kinetic energy per atom $K$ of the excitation. In the linear and nonlinear response regimes where there are no atomic rearrangements, we find that $Q \rightarrow \infty$ (since we do not consider coupling to the environment). We identify a characteristic $K_r$ above which atomic rearrangements occur, and there is significant energy leakage from the fundamental mode to higher frequencies, causing finite $Q$. Thus, $K_r$ is a critical parameter determining resonator performance. We show that $K_r$ decreases as a power-law, $K_r\sim N^{-k},$ with increasing system size $N$, where $k \approx 1.3$. We estimate the critical strain $\langle γ_r \rangle \sim 10^{-8}$ for micron-sized resonators below which atomic rearrangements do not occur, and thus large $Q$-values can be obtained when they are operated below $γ_r$. We find that $K_r$ for amorphous resonators is comparable to that for resonators with crystalline order.

0 citations0 reviewsNo code linkedOpen 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 $p$, the ensemble-averaged static shear modulus $\langle G-G_0 \rangle$ scales with $p^α$, where $α\approx 1$, but above a characteristic pressure $p^{**}$, $\langle G-G_0 \rangle \sim p^β$, where $β\approx 0.5$. However, we find that the shear modulus $G^i$ for an individual packing typically decreases linearly with $p$ 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, $\langle G \rangle$ also depends on a contribution from discontinuous jumps in $\langle G \rangle$ that occur at the transitions between geometrical families. For $p > p^{**}$, geometrical-family and rearrangement contributions to $\langle G \rangle$ are of opposite signs and remain comparable for all system sizes. $\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.

0 citations0 reviewsNo code linkedOpen access