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

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 \beta_\nu^{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 $\beta_\nu^{a,b}\approx 0.3$ for vapor phases to $\beta_\nu^{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 $\beta_\nu^{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

Minkowski tensors and local structure metrics: Amorphous and crystalline sphere packings

Robust and sensitive tools to characterise local structure are essential for investigations of granular or particulate matter. Often local structure metrics derived from the bond network are used for this purpose, in particular Steinhardt's bond-orientational order parameters ql . Here we discuss an alternative method, based on the robust characterisation of the shape of the particles' Voronoi cells, by Minkowski tensors and derived anisotropy measures. We have successfully applied these metrics to quantify structural changes and the onset of crystallisation in random sphere packs. Here we specifically discuss the expectation values of these metrics for simple crystalline unimodal packings of spheres, consisting of single spheres on the points of a Bravais lattice. These data provide an important reference for the discussion of anisotropy values of disordered structures that are typically of relevance in granular systems. This analysis demonstrates that, at least for sufficiently high packing fractions above φ > 0.61, crystalline sphere packs exist whose Voronoi cells are more anisotropic with respect to a volumetric moment tensor than the average value of Voronoi cell anisotropy in random sphere packs

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.Comment: Chapter of the forthcoming book "Tensor Valuations and their Applications in Stochastic Geometry and Imaging" in Lecture Notes in Mathematics edited by Markus Kiderlen and Eva B. Vedel Jense

Millisecond-Timescale Local Network Coding in the Rat Primary Somatosensory Cortex

Correlation among neocortical neurons is thought to play an indispensable role in mediating sensory processing of external stimuli. The role of temporal precision in this correlation has been hypothesized to enhance information flow along sensory pathways. Its role in mediating the integration of information at the output of these pathways, however, remains poorly understood. Here, we examined spike timing correlation between simultaneously recorded layer V neurons within and across columns of the primary somatosensory cortex of anesthetized rats during unilateral whisker stimulation. We used Bayesian statistics and information theory to quantify the causal influence between the recorded cells with millisecond precision. For each stimulated whisker, we inferred stable, whisker-specific, dynamic Bayesian networks over many repeated trials, with network similarity of 83.3±6% within whisker, compared to only 50.3±18% across whiskers. These networks further provided information about whisker identity that was approximately 6 times higher than what was provided by the latency to first spike and 13 times higher than what was provided by the spike count of individual neurons examined separately. Furthermore, prediction of individual neurons' precise firing conditioned on knowledge of putative pre-synaptic cell firing was 3 times higher than predictions conditioned on stimulus onset alone. Taken together, these results suggest the presence of a temporally precise network coding mechanism that integrates information across neighboring columns within layer V about vibrissa position and whisking kinetics to mediate whisker movement by motor areas innervated by layer V

Membrane Potential-Dependent Modulation of Recurrent Inhibition in Rat Neocortex

Dynamic balance of excitation and inhibition is crucial for network stability and cortical processing, but it is unclear how this balance is achieved at different membrane potentials (Vm) of cortical neurons, as found during persistent activity or slow Vm oscillation. Here we report that a Vm-dependent modulation of recurrent inhibition between pyramidal cells (PCs) contributes to the excitation-inhibition balance. Whole-cell recording from paired layer-5 PCs in rat somatosensory cortical slices revealed that both the slow and the fast disynaptic IPSPs, presumably mediated by low-threshold spiking and fast spiking interneurons, respectively, were modulated by changes in presynaptic Vm. Somatic depolarization (>5 mV) of the presynaptic PC substantially increased the amplitude and shortened the onset latency of the slow disynaptic IPSPs in neighboring PCs, leading to a narrowed time window for EPSP integration. A similar increase in the amplitude of the fast disynaptic IPSPs in response to presynaptic depolarization was also observed. Further paired recording from PCs and interneurons revealed that PC depolarization increases EPSP amplitude and thus elevates interneuronal firing and inhibition of neighboring PCs, a reflection of the analog mode of excitatory synaptic transmission between PCs and interneurons. Together, these results revealed an immediate Vm-dependent modulation of cortical inhibition, a key strategy through which the cortex dynamically maintains the balance of excitation and inhibition at different states of cortical activity

Desynchronization of Neocortical Networks by Asynchronous Release of GABA at Autaptic and Synaptic Contacts from Fast-Spiking Interneurons

An activity-dependent long-lasting asynchronous release of GABA from identified fast-spiking inhibitory neurons in the neocortex can impair the reliability and temporal precision of activity in a cortical network

A biophysical model of dynamic balancing of excitation and inhibition in fast oscillatory large-scale networks

Over long timescales, neuronal dynamics can be robust to quite large perturbations, such as changes in white matter connectivity and grey matter structure through processes including learning, aging, development and certain disease processes. One possible explanation is that robust dynamics are facilitated by homeostatic mechanisms that can dynamically rebalance brain networks. In this study, we simulate a cortical brain network using the Wilson-Cowan neural mass model with conduction delays and noise, and use inhibitory synaptic plasticity (ISP) to dynamically achieve a spatially local balance between excitation and inhibition. Using MEG data from 55 subjects we find that ISP enables us to simultaneously achieve high correlation with multiple measures of functional connectivity, including amplitude envelope correlation and phase locking. Further, we find that ISP successfully achieves local E/I balance, and can consistently predict the functional connectivity computed from real MEG data, for a much wider range of model parameters than is possible with a model without ISP

Brief Bursts Self-Inhibit and Correlate the Pyramidal Network

A multi-cell patch clamp study reveals the summation properties of frequency-dependent disynaptic inhibition between neocortical pyramidal cells and shows how brief bursts of activity in a few cells can synchronize the entire microcircuit