Sebastian C. Kapfer

Minimal surface scaffold designs for tissue engineering

Triply-periodic minimal surfaces are shown to be a more versatile source of biomorphic scaffold designs than currently reported in the tissue engineering literature. A scaffold architecture with sheetlike morphology based on minimal surfaces is discussed, with significant structural and mechanical advantages over conventional designs. These sheet solids are porous solids obtained by inflation of cubic minimal surfaces to sheets of finite thickness, as opposed to the conventional network solids where the minimal surface forms the solid/void interface. Using a finite-element approach, the mechanical stiffness of sheet solids is shown to exceed that of conventional network solids for a wide range of volume fractions and material parameters. We further discuss structure-property relationships for mechanical properties useful for custom-designed fabrication by rapid prototyping. Transport properties of the scaffolds are analyzed using Lattice-Boltzmann computations of the fluid permeability. The large number of different minimal surfaces, each of which can be realized as sheet or network solids and at different volume fractions, provides design flexibility essential for the optimization of competing design targets.

Minkowski tensor shape analysis of cellular, granular and porous structures

Predicting physical properties of materials with spatially complex structures is one of the most challenging problems in material science. One key to a better understanding of such materials is the geometric characterization of their spatial structure. Minkowski tensors are tensorial shape indices that allow quantitative characterization of the anisotropy of complex materials and are particularly well suited for developing structure-property relationships for tensor-valued or orientation-dependent physical properties. They are fundamental shape indices, in some sense being the simplest generalization of the concepts of volume, surface and integral curvatures to tensor-valued quantities. Minkowski tensors are based on a solid mathematical foundation provided by integral and stochastic geometry, and are endowed with strong robustness and completeness theorems. The versatile definition of Minkowski tensors applies widely to different types of morphologies, including ordered and disordered structures. Fast linear-time algorithms are available for their computation. This article provides a practical overview of the different uses of Minkowski tensors to extract quantitative physically-relevant spatial structure information from experimental and simulated data, both in 2D and 3D. Applications are presented that quantify (a) alignment of co-polymer films by an electric field imaged by surface force microscopy; (b) local cell anisotropy of spherical bead pack models for granular matter and of closed-cell liquid foam models; (c) surface orientation in open-cell solid foams studied by X-ray tomography; and (d) defect densities and locations in molecular dynamics simulations of crystalline copper.

Two-dimensional melting: from liquid-hexatic coexistence to continuous transitions.

The phase diagram of two-dimensional continuous particle systems is studied using the event-chain Monte Carlo algorithm. For soft disks with repulsive power-law interactions ∝r^{-n} with n≳6, the recently established hard-disk melting scenario (n→∞) holds: a first-order liquid-hexatic and a continuous hexatic-solid transition are identified. Close to n=6, the coexisting liquid exhibits very long orientational correlations, and positional correlations in the hexatic are extremely short. For n≲6, the liquid-hexatic transition is continuous, with correlations consistent with the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario. To illustrate the generality of these results, we demonstrate that Yukawa particles likewise may follow either the KTHNY or the hard-disk melting scenario, depending on the Debye-Hückel screening length as well as on the temperature.

Minkowski tensors of anisotropic spatial structure

This paper 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 generalizations 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 paper 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 formalism more readily accessible for future application in the physical sciences.

Tensorial Minkowski functionals and anisotropy measures for planar patterns

Quantitative measures for anisotropic characteristics of spatial structure are needed when relating the morphology of microstructured heterogeneous materials to tensorial physical properties such as elasticity, permeability and conductance. Tensor‐valued Minkowski functionals, defined in the framework of integral geometry, provide a concise set of descriptors of anisotropic morphology. In this article, we describe the robust computation of these measures for microscopy images and polygonal shapes. We demonstrate their relevance for shape description, their versatility and their robustness by applying them to experimental data sets, specifically microscopy data sets of non‐equilibrium stationary Turing patterns and the shapes of ice grains from Antarctic cores.

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

Polycontinuous geometries for inverse lipid phases with more than two aqueous network domains.

Inverse bicontinuous cubic phases with two aqueous network domains separated by a smooth bilayer are firmly established as equilibrium phases in lipid/water systems. The purpose of this article is to highlight the generalisations of these bicontinuous geometries to polycontinuous geometries, which could be realised as lipid mesophases with three or more network-like aqueous domains separated by a branched bilayer. An analysis of structural homogeneity in terms of bilayer width variations reveals that ordered polycontinuous geometries are likely candidates for lipid mesophase structures, with similar chain packing characteristics to the inverse micellar phases (that once were believed not to exist due to high packing frustration). The average molecular shape required by global geometry to form these multi-network phases is quantified by the surfactant shape parameter, v/(al); we find that it adopts values close to those of the known lipid phases. We specifically analyse the 3etc(187 193) structure of hexagonal symmetry P6(3) /mcm with three aqueous domains, the 3dia(24 220) structure of cubic symmetry I43d composed of three distorted diamond networks, the cubic chiral 4srs(24 208) with cubic symmetry P4232 and the achiral 4srs(5 133) structure of symmetry P42/nbc, each consisting of four intergrown undistorted copies of the srs net (the same net as in the QII(G) gyroid phase). Structural homogeneity is analysed by a medial surface approach assuming that the headgroup interfaces are constant mean curvature surfaces. To facilitate future experimental identification, we provide simulated SAXS scattering patterns that, for the 4srs(24 208) and 3dia(24 220) structures, bear remarkable similarity to those of bicontinuous QII(G)-gyroid and QII(D)-diamond phases, with comparable lattice parameters and only a single peak that cannot be indexed to the well-established structures. While polycontinuous lipid phases have, to date, not been reported, the likelihood of their formation is further indicated by the reported observation of a solid tricontinuous mesoporous silicate structure, termed IBN-9, which formed in the presence of surfactants [Han et al., Nat. Chem., 2009, 1, 123].

Morphology and linear-elastic moduli of random network solids

The effective linear-elastic moduli of disordered network solids are analyzed by voxel-based finite element calculations. We analyze network solids given by Poisson-Voronoi processes and by the structure of collagen fiber networks imaged by confocal microscopy. The solid volume fraction ϕ is varied by adjusting the fiber radius, while keeping the structural mesh or pore size of the underlying network fixed. For intermediate ϕ, the bulk and shear modulus are approximated by empirical power-laws K(phi)proptophin and G(phi)proptophim with n≈1.4 and m≈1.7. The exponents for the collagen and the Poisson-Voronoi network solids are similar, and are close to the values n=1.22 and m=2.11 found in a previous voxel-based finite element study of Poisson-Voronoi systems with different boundary conditions. However, the exponents of these empirical power-laws are at odds with the analytic values of n=1 and m=2, valid for low-density cellular structures in the limit of thin beams. We propose a functional form for K(ϕ) that models the cross-over from a power-law at low densities to a porous solid at high densities; a fit of the data to this functional form yields the asymptotic exponent n≈1.00, as expected. Further, both the intensity of the Poisson-Voronoi process and the collagen concentration in the samples, both of which alter the typical pore or mesh size, affect the effective moduli only by the resulting change of the solid volume fraction. These findings suggest that a network solid with the structure of the collagen networks can be modeled in quantitative agreement by a Poisson-Voronoi process.

Non-universal Voronoi cell shapes in amorphous ellipsoid packs

In particulate systems with short-range interactions, such as granular matter or simple fluids, local structure determines the macroscopic physical properties. We analyse local structure metrics derived from the Voronoi diagram of oblate ellipsoids, for various aspect ratios and global packing fractions φg. We focus on jammed static configurations of frictional ellipsoids, obtained by tomographic imaging and by discrete element method simulations. The rescaled distribution of local packing fractions φl, defined as the ratio of particle volume and its Voronoi cell volume, is found to be independent of the particle aspect ratio, and coincide with results for sphere packs. By contrast, the typical Voronoi cell shape, quantified by the Minkowski tensor anisotropy index β = β02,0, points towards a difference between random packings of spheres and those of oblate ellipsoids. While the average cell shape β of all cells with a given value of is 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.

Set Voronoi diagrams of 3D assemblies of aspherical particles

Abstract Several approaches to quantitative local structure characterization for particulate assemblies, such as structural glasses or jammed packings, use the partition of space provided by the Voronoi diagram. The conventional construction for spherical mono-disperse particles, by which the Voronoi cell of a particle is that of its centre point, cannot be applied to configurations of aspherical or polydisperse particles. Here, we discuss the construction of a Set Voronoi diagram for configurations of aspherical particles in three-dimensional space. The Set Voronoi cell of a given particle is composed of all points in space that are closer to the surface (as opposed to the centre) of the given particle than to the surface of any other; this definition reduces to the conventional Voronoi diagram for the case of mono-disperse spheres. An algorithm for the computation of the Set Voronoi diagram for convex particles is described, as a special case of a Voronoi-based medial axis algorithm, based on a triangulation of the particles’ bounding surfaces. This algorithm is further improved by a pre-processing step based on morphological erosion, which improves the quality of the approximation and circumvents the problems associated with small degrees of particle–particle overlap that may be caused by experimental noise or soft potentials. As an application, preliminary data for the volume distribution of disordered packings of mono-disperse oblate ellipsoids, obtained from tomographic imaging, is computed.