Source author record

Carl F. Schreck

Carl F. Schreck appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

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.

Published work

12 published item(s)

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.

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 $T_0$ and packing fraction deviation above and below jamming onset $Δϕ=ϕ- ϕ_J$ using numerical simulations of purely repulsive frictionless disks. Iso-coordinated solids (ICS) only occur above jamming for $Δϕ> Δϕ_c(T_0)$; they possess average coordination number equal to the isostatic value ($< 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 $ω_k^d$ of the dynamical matrix evaluated at $T_0=0$. Hypo-coordinated solids (HCS) occur both above and below jamming onset within the region defined by $Δϕ> Δϕ^*_-(T_0)$, $Δϕ< Δϕ^*_+(T_0)$, and $Δϕ> Δϕ_{cb}(T_0)$. In this region, the network of interparticle contacts fluctuates with $< 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(ω)$ is not a collection of sharp peaks at $ω_k^d$, and its precise form depends on the measurement method. For $Δϕ> Δϕ_{cb}(T_0)$ and $Δϕ< Δϕ^*_{-}(T_0)$, the system behaves as a hard-particle liquid.

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.

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$) are isostatic and possess the minimum contact number $z_{\rm iso}$ required for them to be collectively jammed, amorphous packings of ellipsoidal particles generally possess fewer contacts than expected for collective jamming ($z < 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-S$ for static packings of ellipsoidal particles into two important components: the stiffness $H$ and stress $S$ matrices. We find that the stiffness matrix possesses $N(z_{\rm iso} - z)$ eigenmodes ${\hat e}_0$ with zero eigenvalues even at finite compression, where $N$ is the number of particles. In addition, these modes ${\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$ for perturbations by amplitude $δ$ along these `quartic' modes, ${\hat e}_0$. These findings illustrate the significant differences between static packings of spherical and ellipsoidal particles.

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.

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.

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.

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\rightarrow \infty$) for all pressures $p$ (or compressions $Δϕ$), and for all $N$ near jamming onset $p\rightarrow 0$. We perform simulations on 2D frictionless bidisperse mechanically stable disk packings over a range of packing fractions $Δϕ= ϕ-ϕ_J$ above jamming onset $ϕ_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$, 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 >\sim Δϕ/N^λ$, where $1 > λ> 0.5$, and thus jammed particulate systems are inherently nonharmonic with no linear vibrational response regime as $N\rightarrow \infty$ over the full range of $Δϕ$, and as $Δϕ\rightarrow 0$ at any $N$.

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.

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 $r$ and initial packing fractions followed by compression and decompression in small steps to reach packing fractions $ϕ_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$. We find that amorphous, isostatic packings exist over a finite range of packing fractions from $ϕ_{\rm min} \le ϕ_J \le ϕ_{\rm max}$ in the large-system limit, with $ϕ_{\rm max} \approx 0.853$. In agreement with previous calculations, we obtain $ϕ_{\rm min} \approx 0.84$ for $r > r^*$, where $r^*$ is the rate above which $ϕ_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.

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.