Source author record

Klaus Mecke

Klaus Mecke 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

25works
16topics
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

25 published item(s)

preprint2021arXiv

Characterization of Anisotropic Gaussian Random Fields by Minkowski Tensors

Gaussian random fields are among the most important models of amorphous spatial structures and appear across length scales in a variety of physical, biological, and geological applications, from composite materials to geospatial data. Anisotropy in such systems can sensitively and comprehensively be characterized by the so-called Minkowski tensors from integral geometry. Here, we analytically calculate the expected Minkowski tensors of arbitrary rank for the level sets of Gaussian random fields. The explicit expressions for interfacial Minkowski tensors are confirmed in detailed simulations. We demonstrate how the Minkowski tensors detect and characterize the anisotropy of the level sets, and we clarify which shape information is contained in the Minkowski tensors of different rank. Using an irreducible representation of the Minkowski tensors in the Euclidean plane, we show that higher-rank tensors indeed contain additional anisotropy information compared to a rank two tensor. Surprisingly, we can nevertheless predict this information from the second-rank tensor if we assume that the random field is Gaussian. This relation between tensors of different rank is independent of the details of the model. It is, therefore, useful for a null hypothesis test that detects non-Gaussianities in anisotropic random fields.

preprint2016arXiv

Cell shape analysis of random tessellations based on Minkowski tensors

To which degree are shape indices of individual cells of a tessellation characteristic for the stochastic process that generates them? Within the context of stochastic geometry and the physics of disordered materials, this corresponds to the question of relationships between different stochastic models. In the context of image analysis of synthetic and biological materials, this question is central to the problem of inferring information about formation processes from spatial measurements of resulting random structures. We address this question by a theory-based simulation study of shape indices derived from Minkowski tensors for a variety of tessellation models. We focus on the relationship between two indices: an isoperimetric ratio of the empirical averages of cell volume and area and the cell elongation quantified by eigenvalue ratios of interfacial Minkowski tensors. Simulation data for these quantities, as well as for distributions thereof and for correlations of cell shape and volume, are presented for Voronoi mosaics of the Poisson point process, determinantal and permanental point processes, and Gibbs hard-core and random sequential absorption processes as well as for Laguerre tessellations of polydisperse spheres and STIT- and Poisson hyperplane tessellations. These data are complemented by mechanically stable crystalline sphere and disordered ellipsoid packings and area-minimising foam models. We find that shape indices of individual cells are not sufficient to unambiguously identify the generating process even amongst this limited set of processes. However, we identify significant differences of the shape indices between many of these tessellation models. Given a realization of a tessellation, these shape indices can narrow the choice of possible generating processes, providing a powerful tool which can be further strengthened by density-resolved volume-shape correlations.

preprint2016arXiv

Genus dependence of the number of (non-)orientable surface triangulations

Topological triangulations of orientable and non-orientable surfaces with arbitrary genus have important applications in quantum geometry, graph theory and statistical physics. However, until now only the asymptotics for 2-spheres are known analytically, and exact counts of triangulations are only available for both small genus and small triangulations. We apply the Wang-Landau algorithm to calculate the number $N(m,h)$ of triangulations for several order of magnitudes in system size $m$ and genus $h$. We verify that the limit of the entropy density of triangulations is independent of genus and orientability and are able to determine the next-to-leading and the next-to-next-to-leading order terms. We conjecture for the number of surface triangulations the asymptotic behavior \begin{equation*} N(m,h) \rightarrow (170.4 \pm 15.1)^h m^{-2(h - 1)/5} \left( \frac{256}{27} \right)^{m / 2}\;, \end{equation*} what might guide a mathematicians proof for the exact asymptotics.

preprint2016arXiv

Spectral properties of unimodular lattice triangulations

Random unimodular lattice triangulations have been recently used as an embedded random graph model, which exhibit a crossover behaviour between an ordered, large-world and a disordered, small-world behaviour. Using the ergodic Pachner flips that transform such triangulations into another and an energy functional that corresponds to the degree distribution variance, Markov chain Monte-Carlo simulations can be applied to study these graphs. Here, we consider the spectra of the adja cency and the Laplacian matrix as well as the algebraic connectivity and the spectral radius. Power law dependencies on the system size can clearly be identified and compared to analytical solutions for periodic ground states. For random triangulations we find a qualitative agreement of the spectral properties with well-known random graph models. In the microcanonical ensemble analytical approximations agree with numerical simulations. In the canonical ensemble a crossover behaviour can be found for the algebraic connectivity and the spectral radius, thus combining large-world and small-world behavior in one model. The considered spectral properties can be applied to transport problems on triangulation graphs and the crossover behaviour allows a tuning of important transport quantities.

preprint2015arXiv

Non-universal Voronoi cell shapes in amorphous ellipsoid packings

In particulate systems with short-range interactions, such as granular matter or simple fluids, local structure plays a pivotal role in determining the macroscopic physical properties. Here, we analyse local structure metrics derived from the Voronoi diagram of configurations of oblate ellipsoids, for various aspect ratios $α$ and global volume fractions $ϕ_g$. We focus on jammed static configurations of frictional ellipsoids, obtained by tomographic imaging and by discrete element method simulations. In particular, we consider the local packing fraction $ϕ_l$, defined as the particle's volume divided by its Voronoi cell volume. We find that the probability $P(ϕ_l)$ for a Voronoi cell to have a given local packing fraction shows the same scaling behaviour as function of $ϕ_g$ as observed for random sphere packs. Surprisingly, this scaling behaviour is further found to be independent of the particle aspect ratio. By contrast, the typical Voronoi cell shape, quantified by the Minkowski tensor anisotropy index $β=β_0^{2,0}$, points towards a significant difference between random packings of spheres and those of oblate ellipsoids. While the average cell shape $β$ of all cells with a given value of $ϕ_l$ is very similar in dense and loose jammed sphere packings, the structure of dense and loose ellipsoid packings differs substantially such that this does not hold true. This non-universality has implications for our understanding of jamming of aspherical particles.

preprint2015arXiv

Unimodular lattice triangulations as small-world and scale-free random graphs

Real-world networks, e.g. the social relations or world-wide-web graphs, exhibit both small-world and scale-free behaviour. We interpret lattice triangulations as planar graphs by identifying triangulation vertices with graph nodes and one-dimensional simplices with edges. Since these triangulations are ergodic with respect to a certain Pachner flip, applying different Monte-Carlo simulations enables us to calculate average properties of random triangulations, as well as canonical ensemble averages using an energy functional that is approximately the variance of the degree distribution. All considered triangulations have clustering coefficients comparable with real world graphs, for the canonical ensemble there are inverse temperatures with small shortest path length independent of system size. Tuning the inverse temperature to a quasi-critical value leads to an indication of scale-free behaviour for degrees $k \geq 5$. Using triangulations as a random graph model can improve the understanding of real-world networks, especially if the actual distance of the embedded nodes becomes important.

preprint2014arXiv

Deriving fundamental measure theory from the virial series: Consistency with the zero-dimensional limit

Fundamental measure theory (FMT) for hard particles has great potential for predicting the phase behavior of colloidal and nanometric shapes. The modern versions of FMT are usually derived from the zero-dimensional limit, a system of at most one particle confined in a collection of cavities in the limit that all cavities shrink to the size of the particle. In [Phys. Rev. E 85, 041150 (2012)], a derivation from an approximated and resummed virial expansion was presented, whose result was not fully consistent with the FMT from the zero-dimensional limit. Here we improve upon this derivation and obtaining exactly the same FMT functional as was obtained earlier from the zero-dimensional limit. As a result, further improvements of FMT based on the virial expansion can now be formulated, some of which we suggest in the outlook.

preprint2014arXiv

Dynamical States in Driven Colloidal Liquid Crystals

We study a model colloidal liquid crystal consisting of hard spherocylinders under the influence of an external aligning potential by Langevin dynamics simulation. The external field that rotates in a plane acts on the orientation of the individual particles and induces a variety of collective nonequilibrium states. We characterize these states by the time-resolved orientational distribution of the particles and explain their origin using the single particle behavior. By varying the external driving frequency and the packing fraction of the spherocylinders we construct the dynamical state diagram.

preprint2014arXiv

Entropy of unimodular Lattice Triangulations

Triangulations are important objects of study in combinatorics, finite element simulations and quantum gravity, where its entropy is crucial for many physical properties. Due to their inherent complex topological structure even the number of possible triangulations is unknown for large systems. We present a novel algorithm for an approximate enumeration which is based on calculations of the density of states using the Wang-Landau flat histogram sampling. For triangulations on two-dimensional integer lattices we achive excellent agreement with known exact numbers of small triangulations as well as an improvement of analytical calculated asymptotics. The entropy density is $C=2.196(3)$ consistent with rigorous upper and lower bounds. The presented numerical scheme can easily be applied to other counting and optimization problems.

preprint2014arXiv

Simulation of a Hard-Spherocylinder Liquid Crystal with the pe

The pe physics engine is validated through the simulation of a liquid crystal model system consisting of hard spherocylinders. For this purpose we evaluate several characteristic parameters of this system, namely the nematic order parameter, the pressure, and the Frank elastic constants. We compare these to the values reported in literature and find a very good agreement, which demonstrates that the pe physics engine can accurately treat such densely packed particle systems. Simultaneously we are able to examine the influence of finite size effects, especially on the evaluation of the Frank elastic constants, as we are far less restricted in system size than earlier simulations.

preprint2014arXiv

Surface tension of isotropic-nematic interfaces: Fundamental Measure Theory for hard spherocylinders

A fluid constituted of hard spherocylinders is studied using a density functional theory for non-spherical hard particles, which can be written as a function of weighted densities. This is based on an extended deconvolution of the Mayer $f$-function for arbitrarily shaped convex hard bodies in tensorial weight functions, which depend each only on the shape and orientation of a single particle. In the course of an examination of the isotropic- nematic interface at coexistence the functional is applied to anisotropic and inhomogeneous problems for the first time. We find good qualitative agreement with other theoretical predictions and also with Monte-Carlo simulations.

preprint2013arXiv

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.

preprint2013arXiv

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.

preprint2013arXiv

Minkowski tensor density formulas for Boolean models

A stationary Boolean model is the union set of random compact particles which are attached to the points of a stationary Poisson point process. For a stationary Boolean model with convex grains we consider a recently developed collection of shape descriptors, the so called Minkowski tensors. By combining spatial and probabilistic averaging we define Minkowski tensor densities of a Boolean model. These densities are global characteristics of the union set which can be estimated from observations. In contrast local characteristics like the mean Minkowski tensor of a single random particle cannot be observed directly, since the particles overlap. We relate the global to the local properties by density formulas for the Minkowski tensors. These density formulas generalize the well known formulas for intrinsic volume densities and are obtained by applying results from translative integral geometry. For an isotropic Boolean model we observe that the Minkowski tensor densities are proportional to the intrinsic volume densities, whereas for a non-isotropic Boolean model this is usually not the case. Our results support the idea that the degree of anisotropy of a Boolean model may be expressed in terms of the Minkowski tensor densities. Furthermore we observe that for smooth grains the mean curvature radius function of a particle can be reconstructed from the Minkowski tensor densities. In a simulation study we determine numerically Minkowski tensor densities for non-isotropic Boolean models based on ellipses and on rectangles in two dimensions and find excellent agreement with the derived analytic density formulas. The tensor densities can be used to characterize the orientational distribution of the grains and to estimate model parameters for non-isotropic distributions.

preprint2013arXiv

Minkowski Tensors of Anisotropic Spatial Structure

This article describes the theoretical foundation of and explicit algorithms for a novel approach to morphology and anisotropy analysis of complex spatial structure using tensor-valued Minkowski functionals, the so-called Minkowski tensors. Minkowski tensors are generalisations of the well-known scalar Minkowski functionals and are explicitly sensitive to anisotropic aspects of morphology, relevant for example for elastic moduli or permeability of microstructured materials. Here we derive explicit linear-time algorithms to compute these tensorial measures for three-dimensional shapes. These apply to representations of any object that can be represented by a triangulation of its bounding surface; their application is illustrated for the polyhedral Voronoi cellular complexes of jammed sphere configurations, and for triangulations of a biopolymer fibre network obtained by confocal microscopy. The article further bridges the substantial notational and conceptual gap between the different but equivalent approaches to scalar or tensorial Minkowski functionals in mathematics and in physics, hence making the mathematical measure theoretic method more readily accessible for future application in the physical sciences.

preprint2013arXiv

Network-like propagation of cell-level stress in random foams

Quasistatic simple shearing flow of random monodisperse soap froth is investigated by analyzing Surface Evolver simulations of spatially periodic foams. Elastic-plastic behavior is caused by irreversible topological rearrangements (T1s) that occur when Plateau's laws are violated; the first T1s occur at the elastic limit and at large strains frequent cascades of T1s, composed of one or more individual T1s, sustain the yield-stress plateau. The stress and shape anisotropy of individual cells is quantified by $Q$, a scalar measure derived from the interface tensor that gauges each cell's contribution to the global stress. During each T1 cascade, the connected set of cells with decreasing $Q$, called the \textit{stress release domain}, is network-like and highly non-local. Geometrically, the network-like nature of the stress release domains is corroborated through morphological analysis using the Euler characteristic. The stress release domain is distinctly different from the set of cells that change topology during a T1 cascade. Our results highlight the unique rheological behavior of foams, where complex large-scale cooperative rearrangements of foam cells are observed as a consequence of distinctly local events.

preprint2013arXiv

Percolation transition by random vertex splitting of diamond networks?

We propose a statistical model defined on the three-dimensional diamond network where the splitting of randomly selected nodes leads to a spatially disordered network, with decreasing degree of connectivity. The terminal state, that is reached when all nodes have been split, is a dense configuration of self-avoiding walks on the diamond network. Starting from the crystallographic diamond network, each of the four-coordinated nodes is replaced with probability p by a pair of two edges, each connecting a pair of the adjacent vertices. For all values of 0 <= p <= 1 the network percolates, yet the fraction fp of the system that belongs to a percolating cluster drops sharply at pc=1 to a finite value. This transition has the signature of a phase transition with scaling exponents for p close to pc that are different from the critical exponents of the second-order phase transition of standard percolation models. As is the case for percolation transitions, this transition significantly affects the mechanical properties of linear-elastic realisations, obtained by replacing edges with solid circular struts to give an effective density phi. Finite element methods demonstrate that, as a low-density cellular structure, the bulk modulus K shows a cross-over from a compression-dominated behaviour (with K proportional to phi) at p=0 to a bending-dominated behaviour (with quadratic dependence of K on phi) at p=1.

preprint2013arXiv

Time-Dependent Hartree-Fock Approach to Nuclear Pasta at Finite Temperature

We present simulations of neutron-rich matter at subnuclear densities, like supernova matter, with the time-dependent Hartree-Fock approximation at temperatures of several MeV. The initial state consists of $α$ particles randomly distributed in space that have a Maxwell-Boltzmann distribution in momentum space. Adding a neutron background initialized with Fermi distributed plane waves the calculations reflect a reasonable approximation of astrophysical matter. This matter evolves into spherical, rod-like, and slab-like shapes and mixtures thereof. The simulations employ a full Skyrme interaction in a periodic three-dimensional grid. By an improved morphological analysis based on Minkowski functionals, all eight pasta shapes can be uniquely identified by the sign of only two valuations, namely the Euler characteristic and the integral mean curvature. In addition, we propose the variance in the cell density distribution as a measure to distinguish pasta matter from uniform matter.

preprint2012arXiv

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 ϕ_c \approx 0.649, distinctly higher than previously reported values for the contested random close packing limit. At ϕ_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 ϕ_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.

preprint2012arXiv

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.

preprint2012arXiv

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)].

preprint2011arXiv

All good things come in threes - Three beads learn to swim with lattice Boltzmann and a rigid body solver

We simulate the self-propulsion of devices in a fluid in the regime of low Reynolds numbers. Each device consists of three bodies (spheres or capsules) connected with two damped harmonic springs. Sinusoidal driving forces compress the springs which are resolved within a rigid body physics engine. The latter is consistently coupled to a 3D lattice Boltzmann framework for the fluid dynamics. In simulations of three-sphere devices, we find that the propulsion velocity agrees well with theoretical predictions. In simulations where some or all spheres are replaced by capsules, we find that the asymmetry of the design strongly affects the propelling efficiency.

preprint2010arXiv

Anomalous transport of a tracer on percolating clusters

We investigate the dynamics of a single tracer exploring a course of fixed obstacles in the vicinity of the percolation transition for particles confined to the infinite cluster. The mean-square displacement displays anomalous transport, which extends to infinite times precisely at the critical obstacle density. The slowing down of the diffusion coefficient exhibits power-law behavior for densities close to the critical point and we show that the mean-square displacement fulfills a scaling hypothesis. Furthermore, we calculate the dynamic conductivity as response to an alternating electric field. Last, we discuss the non-gaussian parameter as an indicator for heterogeneous dynamics.

preprint2010arXiv

Local Anisotropy of Fluids using Minkowski Tensors

Statistics of the free volume available to individual particles have previously been studied for simple and complex fluids, granular matter, amorphous solids, and structural glasses. Minkowski tensors provide a set of shape measures that are based on strong mathematical theorems and easily computed for polygonal and polyhedral bodies such as free volume cells (Voronoi cells). They characterize the local structure beyond the two-point correlation function and are suitable to define indices $0\leq β_ν^{a,b}\leq 1$ of local anisotropy. Here, we analyze the statistics of Minkowski tensors for configurations of simple liquid models, including the ideal gas (Poisson point process), the hard disks and hard spheres ensemble, and the Lennard-Jones fluid. We show that Minkowski tensors provide a robust characterization of local anisotropy, which ranges from $β_ν^{a,b}\approx 0.3$ for vapor phases to $β_ν^{a,b}\to 1$ for ordered solids. We find that for fluids, local anisotropy decreases monotonously with increasing free volume and randomness of particle positions. Furthermore, the local anisotropy indices $β_ν^{a,b}$ are sensitive to structural transitions in these simple fluids, as has been previously shown in granular systems for the transition from loose to jammed bead packs.

preprint2006arXiv

Radial Distribution Function for Semiflexible Polymers Confined in Microchannels

An analytic expression is derived for the distribution $G(\vec{R})$ of the end-to-end distance $\vec{R}$ of semiflexible polymers in external potentials to elucidate the effect of confinement on the mechanical and statistical properties of biomolecules. For parabolic confinement the result is exact whereas for realistic potentials a self-consistent ansatz is developed, so that $G(\vec{R})$ is given explicitly even for hard wall confinement. The theoretical result is in excellent quantitative agreement with fluorescence microscopy data for actin filaments confined in rectangularly shaped microchannels. This allows an unambiguous determination of persistence length $L_P$ and the dependence of statistical properties such as Odijk's deflection length $λ$ on the channel width $D$. It is shown that neglecting the effect of confinement leads to a significant overestimation of bending rigidities for filaments.