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

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 \phi, temperature T, and aspect ratio \alpha. We are able to collapse the \phi and T-dependent structural relaxation times for disks, and dimers and ellipses over a wide range of \alpha, onto a universal scaling function {\cal F}_{\pm}(|\phi-\phi_0|,T,\alpha), 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 \alpha-dependent scaling exponents \mu and \delta and packing fraction \phi_0(\alpha) 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 \phi_0, m(\phi_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.Comment: 10 pages, 17 figures, and 1 tabl

Attributing Tropical Cyclogenesis to Equatorial Waves in the Western North Pacific

The direct influences of equatorial waves on the genesis of tropical cyclones are evaluated. Tropical cyclogenesis is attributed to an equatorial wave when the filtered rainfall anomaly exceeds a threshold value at the genesis location. For an attribution threshold of 3 mm/day, 51% of warm season western North Pacific tropical cyclones are attributed to tropical depression (TD)-type disturbances, 29% to equatorial Rossby waves, 26% to mixed Rossby-Gravity waves, 23% to Kelvin waves, 13% to the Madden-Julian oscillation (MJO), and 19% are not attributed to any equatorial wave. The fraction of tropical cyclones attributed to TD-type disturbances is consistent with previous findings. Past studies have also demonstrated that the MJO significantly modulates tropical cyclogenesis, but fewer storms are attributed to the MJO than any other wave type. This disparity arises from the difference between attribution and modulation. The MJO produces broad regions of favorable conditions for cyclogenesis, but the MJO alone might not determine when and where a storm will develop within these regions. Tropical cyclones contribute less than 17% of the power in any portion of the equatorial wave spectrum because tropical cyclones are relatively uncommon equatorward of 15deg latitude. In regions where they are active, however, tropical cyclones can contribute more than 20% of the warm season rainfall and up to 50% of the total variance. Tropical cyclone-related anomalies can significantly contaminate wave-filtered precipitation at the location of genesis. To mitigate this effect, the tropical cyclone-related rainfall anomalies were removed before filtering in this study

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 $\Delta \phi$. While amorphous packings of spherical particles at jamming onset ($\Delta \phi=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 $\Delta \phi$, 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 $\delta^4$ for perturbations by amplitude $\delta$ along these `quartic' modes, ${\hat e}_0$. These findings illustrate the significant differences between static packings of spherical and ellipsoidal particles.Comment: 18 pages, 21 figure