Source author record

Sebastian C. Kapfer

Sebastian C. Kapfer 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

14works
7topics
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

14 published item(s)

preprint2019arXiv

Phasonic Diffusion and Self-confinement of Decagonal Quasicrystals in Hyperspace

We introduce a novel simulation method that is designed to explore fluctuations of the phasonic degrees of freedom in decagonal colloidal quasicrystals. Specifically, we attain and characterise thermal equilibrium of the phason ensemble via Monte Carlo simulations with particle motions restricted to elementary phasonic flips. We find that, at any temperature, the random tiling ensemble is strongly preferred over the minimum phason-strain quasicrystal. Phasonic flips are the dominant carriers of diffusive mass transport in physical space. Sub-diffusive transients suggest cooperative flip behaviour on short time scales. In complementary space, particle mobility is geometrically restricted to a thin ring around the acceptance domain, resulting in self-confinement and persistent phasonic order.

preprint2016arXiv

Cell-veto Monte Carlo algorithm for long-range systems

We present a rigorous efficient event-chain Monte Carlo algorithm for long-range interacting particle systems. Using a cell-veto scheme within the factorized Metropolis algorithm, we compute each single-particle move with a fixed number of operations. For slowly decaying potentials such as Coulomb interactions, screening line charges allow us to take into account periodic boundary conditions. We discuss the performance of the cell-veto Monte Carlo algorithm for general inverse-power-law potentials, and illustrate how it provides a new outlook on one of the prominent bottlenecks in large-scale atomistic Monte Carlo simulations.

preprint2016arXiv

Cuddling Ellipsoids: Densest local structures of uniaxial ellipsoids

Connecting the collective behavior of disordered systems with local structure on the particle scale is an important challenge, for example in granular and glassy systems. Compounding complexity, in many scientific and industrial applications, particles are polydisperse, aspherical or even of varying shape. Here, we investigate a generalization of the classical kissing problem in order to understand the local building blocks of packings of aspherical grains. We numerically determine the densest local structures of uniaxial ellipsoids by minimizing the Set Voronoi cell volume around a given particle. Depending on the particle aspect ratio, different local structures are observed and classified by symmetry and Voronoi coordination number. In extended disordered packings of frictionless particles, knowledge of the densest structures allows to rescale the Voronoi volume distributions onto the single-parameter family of $k$-Gamma distributions. Moreover, we find that approximate icosahedral clusters are found in random packings, while the optimal local structures for more aspherical particles are not formed.

preprint2016arXiv

Splitting of the universality class of anomalous transport in crowded media

We investigate the emergence of subdiffusive transport by obstruction in continuum models for molecular crowding. While the underlying percolation transition for the accessible space displays universal behavior, the dynamic properties depend in a subtle non-universal way on the transport through narrow channels. At the same time, the different universality classes are robust with respect to introducing correlations in the obstacle matrix as we demonstrate for quenched hard-sphere liquids as underlying structures. Our results confirm that the microscopic dynamics can dominate the relaxational behavior even at long times, in striking contrast to glassy dynamics.

preprint2015arXiv

2D Melting: From Liquid-Hexatic Coexistence to Continuous Transitions

The phase diagram of two-dimensional continuous particle systems is studied using Event-Chain Monte Carlo. For soft disks with repulsive power-law interactions $\propto r^{-n}$ with $n \gtrsim 6$, the recently established hard-disk melting scenario ($n \to \infty$) 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\lesssim 6$, the liquid-hexatic transition is continuous, with correlations consistent with the Kosterlitz-Thouless-Halperin-Nelson-Yong (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.

preprint2015arXiv

Liquid-solid transitions in the three-body hard-core model

We determine the phase diagram for a generalisation of two-and three-dimensional hard spheres: a classical system with three-body interactions realised as a hard cut-off on the mean-square distance for each triplet of particles. Quantum versions of this model are important in the context of the unitary Bose gas, which is currently under close theoretical and experimental scrutiny. In two dimensions, the three-body hard-core model possesses a conventional atomic liquid phase and a peculiar solid phase formed by dimers. These dimers interact effectively as hard disks. In three dimensions, the solid phase consists of isolated atoms that arrange in a simple-hexagonal lattice.

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.

preprint2013arXiv

Generalized event-chain Monte Carlo: Constructing rejection-free global-balance algorithms from infinitesimal steps

In this article, we present an event-driven algorithm that generalizes the recent hard-sphere event-chain Monte Carlo method without introducing discretizations in time or in space. A factorization of the Metropolis filter and the concept of infinitesimal Monte Carlo moves are used to design a rejection-free Markov-chain Monte Carlo algorithm for particle systems with arbitrary pairwise interactions. The algorithm breaks detailed balance, but satisfies maximal global balance and performs better than the classic, local Metropolis algorithm in large systems. The new algorithm generates a continuum of samples of the stationary probability density. This allows us to compute the pressure and stress tensor as a byproduct of the simulation without any additional computations.

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

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

Sampling from a polytope and hard-disk Monte Carlo

The hard-disk problem, the statics and the dynamics of equal two-dimensional hard spheres in a periodic box, has had a profound influence on statistical and computational physics. Markov-chain Monte Carlo and molecular dynamics were first discussed for this model. Here we reformulate hard-disk Monte Carlo algorithms in terms of another classic problem, namely the sampling from a polytope. Local Markov-chain Monte Carlo, as proposed by Metropolis et al. in 1953, appears as a sequence of random walks in high-dimensional polytopes, while the moves of the more powerful event-chain algorithm correspond to molecular dynamics evolution. We determine the convergence properties of Monte Carlo methods in a special invariant polytope associated with hard-disk configurations, and the implications for convergence of hard-disk sampling. Finally, we discuss parallelization strategies for event-chain Monte Carlo and present results for a multicore implementation.

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

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

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.