Reply to "Comment on `Jamming at zero temperature and zero applied stress: The epitome of disorder' "

We answer the questions raised by Donev, Torquato, Stillinger, and Connelly in their "Comment on "Jamming at zero temperature and zero applied stress: The epitome of disorder.' " We emphasize that we follow a fundamentally different approach than they have done to reinterpret random close packing in terms of the "maximally random jammed" framework. We define the "maximally random jammed packing fraction" to be where the largest number of initial states, chosen completely randomly, have relaxed final states at the jamming threshold in the thermodynamic limit. Thus, we focus on an ensemble of states at the jamming threshold, while DTSC are interested in determining the amount of order and degree of jamming for a particular configuration. We also argue that soft-particle systems are as "clean" as those using hard spheres for studying jammed packings and point out the benefits of using soft potentials

Stabilization of nonlinear velocity profiles in athermal systems undergoing planar shear flow

We perform molecular dynamics simulations of model granular systems undergoing boundary-driven planar shear flow in two spatial dimensions with the goal of developing a more complete understanding of how dense particulate systems respond to applied shear. In particular, we are interested in determining when these systems will possess linear velocity profiles and when they will develop highly localized velocity profiles in response to shear. In previous work on similar systems we showed that nonlinear velocity profiles form when the speed of the shearing boundary exceeds the speed of shear waves in the material. However, we find that nonlinear velocity profiles in these systems are unstable at very long times. The degree of nonlinearity slowly decreases in time; the velocity profiles become linear when the granular temperature and density profiles are uniform across the system at long times. We measure the time $t_l$ required for the velocity profiles to become linear and find that $t_l$ increases as a power-law with the speed of the shearing boundary and increases rapidly as the packing fraction approaches random close packing. We also performed simulations in which differences in the granular temperature across the system were maintained by vertically vibrating one of the boundaries during shear flow. We find that nonlinear velocity profiles form and are stable at long times if the difference in the granular temperature across the system exceeds a threshold value that is comparable to the glass transition temperature in an equilibrium system at the same average density. Finally, the sheared and vibrated systems form stable shear bands, or highly localized velocity profiles, when the applied shear stress is lowered below the yield stress of the static part of the system.Comment: 11 pages, 14 figure

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 $\phi_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 $\phi_J$. We find that amorphous, isostatic packings exist over a finite range of packing fractions from $\phi_{\rm min} \le \phi_J \le \phi_{\rm max}$ in the large-system limit, with $\phi_{\rm max} \approx 0.853$. In agreement with previous calculations, we obtain $\phi_{\rm min} \approx 0.84$ for $r > r^*$, where $r^*$ is the rate above which $\phi_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.Comment: 11 pages, 15 figure

Determining the Onset of Hydrodynamic Erosion in Turbulent Flow

We revisit the longstanding question of the onset of sediment transport driven by a turbulent fluid flow via laboratory measurements. We use particle tracking velocimetry to quantify the fluid flow as well as the motion of individual grains. As we increase the flow speed above the transition to sediment transport, we observe that an increasing fraction of grains are transported downstream, although the average downstream velocity of the transported grains remains roughly constant. However, we find that the fraction of mobilized grains does not vanish sharply at a critical flow rate. Additionally, the distribution of the fluctuating velocities of non-transported grains becomes broader with heavier tails, meaning that unambiguously separating mobile and static grains is not possible. As an alternative approach, we quantify the statistics of grain velocities by using a mixture model consisting of two forms for the grain velocities: a decaying-exponential tail, which represents grains transported downstream, and a peaked distribution centered at zero velocity, which represents grains that fluctuate due to the turbulent flow but remain in place. Our results suggest that more sophisticated statistical measures may be required to quantify grain motion near the onset of sediment transport, particularly in the presence of turbulence.Comment: 9 pages, 7 figure

Reconciling Semiclassical and Bohmian Mechanics: II. Scattering states for discontinuous potentials

In a previous paper [J. Chem. Phys. 121 4501 (2004)] a unique bipolar decomposition, Psi = Psi1 + Psi2 was presented for stationary bound states Psi of the one-dimensional Schroedinger equation, such that the components Psi1 and Psi2 approach their semiclassical WKB analogs in the large action limit. Moreover, by applying the Madelung-Bohm ansatz to the components rather than to Psi itself, the resultant bipolar Bohmian mechanical formulation satisfies the correspondence principle. As a result, the bipolar quantum trajectories are classical-like and well-behaved, even when Psi has many nodes, or is wildly oscillatory. In this paper, the previous decomposition scheme is modified in order to achieve the same desirable properties for stationary scattering states. Discontinuous potential systems are considered (hard wall, step, square barrier/well), for which the bipolar quantum potential is found to be zero everywhere, except at the discontinuities. This approach leads to an exact numerical method for computing stationary scattering states of any desired boundary conditions, and reflection and transmission probabilities. The continuous potential case will be considered in a future publication.Comment: 18 pages, 8 figure

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.Comment: 18 pages, 23 figures, 5 table

Discovery of high-entropy ceramics via machine learning

AbstractAlthough high-entropy materials are attracting considerable interest due to a combination of useful properties and promising applications, predicting their formation remains a hindrance for rational discovery of new systems. Experimental approaches are based on physical intuition and/or expensive trial and error strategies. Most computational methods rely on the availability of sufficient experimental data and computational power. Machine learning (ML) applied to materials science can accelerate development and reduce costs. In this study, we propose an ML method, leveraging thermodynamic and compositional attributes of a given material for predicting the synthesizability (i.e., entropy-forming ability) of disordered metal carbides. The relative importance of the thermodynamic and compositional features for the predictions are then explored. The approach’s suitability is demonstrated by comparing values calculated with density functional theory to ML predictions. Finally, the model is employed to predict the entropy-forming ability of 70 new compositions; several predictions are validated by additional density functional theory calculations and experimental synthesis, corroborating the effectiveness in exploring vast compositional spaces in a high-throughput manner. Importantly, seven compositions are selected specifically, because they contain all three of the Group VI elements (Cr, Mo, and W), which do not form room temperature-stable rock-salt monocarbides. Incorporating the Group VI elements into the rock-salt structure provides further opportunity for tuning the electronic structure and potentially material performance

Viral MicroRNA Effects on Pathogenesis of Polyomavirus SV40 Infections in Syrian Golden Hamsters

Shaojie Zhang, Vojtech Sroller, Preeti Zanwar, Steven J. Halvorson, Nadim J. Ajami, Corey W. Hecksel, Jody L. Swain, Connie Wong, Janet S. Butel, Department of Molecular Virology and Microbiology, Baylor College of Medicine, Houston, Texas, United States of AmericaChun Jung Chen, Christopher S. Sullivan, Department of Molecular Biosciences, The University of Texas at Austin, Austin, Texas, United States of AmericaJody L. Swain, Center for Comparative Medicine, Baylor College of Medicine, Houston, Texas, United States of AmericaEffects of polyomavirus SV40 microRNA on pathogenesis of viral infections in vivo are not known. Syrian golden hamsters are the small animal model for studies of SV40. We report here effects of SV40 microRNA and influence of the structure of the regulatory region on dynamics of SV40 DNA levels in vivo. Outbred young adult hamsters were inoculated by the intracardiac route with 1×107 plaque-forming units of four different variants of SV40. Infected animals were sacrificed from 3 to 270 days postinfection and viral DNA loads in different tissues determined by quantitative real-time polymerase chain reaction assays. All SV40 strains displayed frequent establishment of persistent infections and slow viral clearance. SV40 had a broad tissue tropism, with infected tissues including liver, kidney, spleen, lung, and brain. Liver and kidney contained higher viral DNA loads than other tissues; kidneys were the preferred site for long-term persistent infection although detectable virus was also retained in livers. Expression of SV40 microRNA was demonstrated in wild-type SV40-infected tissues. MicroRNA-negative mutant viruses consistently produced higher viral DNA loads than wild-type SV40 in both liver and kidney. Viruses with complex regulatory regions displayed modestly higher viral DNA loads in the kidney than those with simple regulatory regions. Early viral transcripts were detected at higher levels than late transcripts in liver and kidney. Infectious virus was detected infrequently. There was limited evidence of increased clearance of microRNA-deficient viruses. Wild-type and microRNA-negative mutants of SV40 showed similar rates of transformation of mouse cells in vitro and tumor induction in weanling hamsters in vivo. This report identified broad tissue tropism for SV40 in vivo in hamsters and provides the first evidence of expression and function of SV40 microRNA in vivo. Viral microRNA dampened viral DNA levels in tissues infected by SV40 strains with simple or complex regulatory regions.This work was supported in part by research grants R01 CA134524 (JSB) and R01 AI077746 (CSS) from the National Institutes of Health. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Molecular BiosciencesEmail: [email protected]