Shortcomings of the Bond Orientational Order Parameters for the Analysis of Disordered Particulate Matter

Local structure characterization with the bond-orientational order parameters q4, q6, ... introduced by Steinhardt et al. has become a standard tool in condensed matter physics, with applications including glass, jamming, melting or crystallization transitions and cluster formation. Here we discuss two fundamental flaws in the definition of these parameters that significantly affect their interpretation for studies of disordered systems, and offer a remedy. First, the definition of the bond-orientational order parameters considers the geometrical arrangement of a set of neighboring spheres NN(p) around a given central particle p; we show that procedure to select the spheres constituting the neighborhood NN(p) can have greater influence on both the numerical values and qualitative trend of ql than a change of the physical parameters, such as packing fraction. Second, the discrete nature of neighborhood implies that NN(p) is not a continuous function of the particle coordinates; this discontinuity, inherited by ql, leads to a lack of robustness of the ql as structure metrics. Both issues can be avoided by a morphometric approach leading to the robust Minkowski structure metrics ql'. These ql' are of a similar mathematical form as the conventional bond-orientational order parameters and are mathematically equivalent to the recently introduced Minkowski tensors [Europhys. Lett. 90, 34001 (2010); Phys. Rev. E. 85, 030301 (2012)]

Space-resolved dynamics of a tracer in a disordered solid

The dynamics of a tracer particle in a glassy matrix of obstacles displays slow complex transport as the free volume approaches a critical value and the void space falls apart. We investigate the emerging subdiffusive motion of the test particle by extensive molecular dynamics simulations and characterize the spatio-temporal transport in terms of two-time correlation functions, including the time-dependent diffusion coefficient as well as the wavenumber-dependent intermediate scattering function. We rationalize our findings within the framework of critical phenomena and compare our data to a dynamic scaling theory.Comment: 10 pages, 7 figures, submitted to Journal of Non-Crystalline Solid

Group Theory of Chiral Photonic Crystals with 4-fold Symmetry: Band Structure and S-Parameters of Eight-Fold Intergrown Gyroid Nets

The Single Gyroid, or srs, nanostructure has attracted interest as a circular-polarisation sensitive photonic material. We develop a group theoretical and scattering matrix method, applicable to any photonic crystal with symmetry I432, to demonstrate the remarkable chiral-optical properties of a generalised structure called 8-srs, obtained by intergrowth of eight equal-handed srs nets. Exploiting the presence of four-fold rotations, Bloch modes corresponding to the irreducible representations E- and E+ are identified as the sole and non-interacting transmission channels for right- and left-circularly polarised light, respectively. For plane waves incident on a finite slab of the 8-srs, the reflection rates for both circular polarisations are identical for all frequencies and transmission rates are identical up to a critical frequency below which scattering in the far field is restricted to zero grating order. Simulations show the optical activity of the lossless dielectric 8-srs to be large, comparable to metallic metamaterials, demonstrating its potential as a nanofabricated photonic material

Group Theory of Circular-Polarization Effects in Chiral Photonic Crystals with Four-Fold Rotation Axes, Applied to the Eight-Fold Intergrowth of Gyroid Nets

We use group or representation theory and scattering matrix calculations to derive analytical results for the band structure topology and the scattering parameters, applicable to any chiral photonic crystal with body-centered cubic symmetry I432 for circularly-polarised incident light. We demonstrate in particular that all bands along the cubic [100] direction can be identified with the irreducible representations E+/-,A and B of the C4 point group. E+ and E- modes represent the only transmission channels for plane waves with wave vector along the ? line, and can be identified as non-interacting transmission channels for right- (E-) and left-circularly polarised light (E+), respectively. Scattering matrix calculations provide explicit relationships for the transmission and reflectance amplitudes through a finite slab which guarantee equal transmission rates for both polarisations and vanishing ellipticity below a critical frequency, yet allowing for finite rotation of the polarisation plane. All results are verified numerically for the so-called 8-srs geometry, consisting of eight interwoven equal-handed dielectric Gyroid networks embedded in air. The combination of vanishing losses, vanishing ellipticity, near-perfect transmission and optical activity comparable to that of metallic meta-materials makes this geometry an attractive design for nanofabricated photonic materials

Pomelo, a tool for computing Generic Set Voronoi Diagrams of Aspherical Particles of Arbitrary Shape

We describe the development of a new software tool, called "Pomelo", for the calculation of Set Voronoi diagrams. Voronoi diagrams are a spatial partition of the space around the particles into separate Voronoi cells, e.g. applicable to granular materials. A generalization of the conventional Voronoi diagram for points or monodisperse spheres is the Set Voronoi diagram, also known as navigational map or tessellation by zone of influence. In this construction, a Set Voronoi cell contains the volume that is closer to the surface of one particle than to the surface of any other particle. This is required for aspherical or polydisperse systems. Pomelo is designed to be easy to use and as generic as possible. It directly supports common particle shapes and offers a generic mode, which allows to deal with any type of particles that can be described mathematically. Pomelo can create output in different standard formats, which allows direct visualization and further processing. Finally, we describe three applications of the Set Voronoi code in granular and soft matter physics, namely the problem of packings of ellipsoidal particles with varying degrees of particle-particle friction, mechanical stable packings of tetrahedra and a model for liquid crystal systems of particles with shapes reminiscent of pearsComment: 4 pages, 9 figures, Submitted to Powders and Grains 201

Low-temperature statistical mechanics of the QuanTizer problem: fast quenching and equilibrium cooling of the three-dimensional Voronoi Liquid

The Quantizer problem is a tessellation optimisation problem where point configurations are identified such that the Voronoi cells minimise the second moment of the volume distribution. While the ground state (optimal state) in 3D is almost certainly the body-centered cubic lattice, disordered and effectively hyperuniform states with energies very close to the ground state exist that result as stable states in an evolution through the geometric Lloyd's algorithm [Klatt et al. Nat. Commun., 10, 811 (2019)]. When considered as a statistical mechanics problem at finite temperature, the same system has been termed the 'Voronoi Liquid' by [Ruscher et al. EPL 112, 66003 (2015)]. Here we investigate the cooling behaviour of the Voronoi liquid with a particular view to the stability of the effectively hyperuniform disordered state. As a confirmation of the results by Ruscher et al., we observe, by both molecular dynamics and Monte Carlo simulations, that upon slow quasi-static equilibrium cooling, the Voronoi liquid crystallises from a disordered configuration into the body-centered cubic configuration. By contrast, upon sufficiently fast non-equilibrium cooling (and not just in the limit of a maximally fast quench) the Voronoi liquid adopts similar states as the effectively hyperuniform inherent structures identified by Klatt et al. and prevents the ordering transition into a BCC ordered structure. This result is in line with the geometric intuition that the geometric Lloyd's algorithm corresponds to a type of fast quench.Comment: 11 pages, 6 figure

Permeability of porous materials determined from the Euler characteristic

We study the permeability of quasi two-dimensional porous structures of randomly placed overlapping monodisperse circular and elliptical grains. Measurements in microfluidic devices and lattice Boltzmann simulations demonstrate that the permeability is determined by the Euler characteristic of the conducting phase. We obtain an expression for the permeability that is independent of the percolation threshold and shows agreement with experimental and simulated data over a wide range of porosities. Our approach suggests that the permeability explicitly depends on the overlapping probability of grains rather than their shape

Jammed Spheres: Minkowski Tensors Reveal Onset of Local Crystallinity

The local structure of disordered jammed packings of monodisperse spheres without friction, generated by the Lubachevsky-Stillinger algorithm, is studied for packing fractions above and below 64%. The structural similarity of the particle environments to fcc or hcp crystalline packings (local crystallinity) is quantified by order metrics based on rank-four Minkowski tensors. We find a critical packing fraction \phi_c \approx 0.649, distinctly higher than previously reported values for the contested random close packing limit. At \phi_c, the probability of finding local crystalline configurations first becomes finite and, for larger packing fractions, increases by several orders of magnitude. This provides quantitative evidence of an abrupt onset of local crystallinity at \phi_c. We demonstrate that the identification of local crystallinity by the frequently used local bond-orientational order metric q_6 produces false positives, and thus conceals the abrupt onset of local crystallinity. Since the critical packing fraction is significantly above results from mean-field analysis of the mechanical contacts for frictionless spheres, it is suggested that dynamic arrest due to isostaticity and the alleged geometric phase transition in the Edwards framework may be disconnected phenomena