Local Stereology of Tensors of Convex Bodies

In this paper, we present local stereological estimators of Minkowski tensors defined on convex bodies in ℝ d . Special cases cover a number of well-known local stereological estimators of volume and surface area in ℝ3, but the general set-up also provides new local stereological estimators of various types of centres of gravity and tensors of rank two. Rank two tensors can be represented as ellipsoids and contain information about shape and orientation. The performance of some of the estimators of centres of gravity and volume tensors of rank two is investigated by simulation.

Estimating Particle Shape and Orientation Using Volume Tensors

We develop statistical procedures for estimating shape and orientation of arbitrary three-dimensional particles. We focus on the case where particles cannot be observed directly, but only via sections. Volume tensors are used for describing particle shape and orientation, and we derive stereological estimators of the tensors. These estimators are combined to provide consistent estimators of the moments of the so-called particle cover density. The covariance structure associated with the particle cover density depends on the orientation and shape of the particles. For instance, if the distribution of the typical particle is invariant under rotations, then the covariance matrix is proportional to the identity matrix. We develop a non-parametric test for such isotropy. A flexible Lévy-based particle model is proposed, which may be analysed using a generalized method of moments in which the volume tensors enter. The developed methods are used to study the cell organization in the human brain cortex.

Stereological estimation of particle shape and orientation from volume tensors

In the present paper, we describe new robust methods of estimating cell shape and orientation in 3D from sections. The descriptors of 3D cell shape and orientation are based on volume tensors which are used to construct an ellipsoid, the Miles ellipsoid, approximating the average cell shape and orientation in 3D. The estimators of volume tensors are based on observations in several optical planes through sampled cells. This type of geometric sampling design is known as the optical rotator. The statistical behaviour of the estimator of the Miles ellipsoid is studied under a flexible model for 3D cell shape and orientation. In a simulation study, the lengths of the axes of the Miles ellipsoid can be estimated with CVs of about 2% if 100 cells are sampled. Finally, we illustrate the use of the developed methods in an example, involving neurons in the medial prefrontal cortex of rat.

Rotational integral geometry of tensor valuations

We derive a new rotational Crofton formula for Minkowski tensors. In special cases, this formula gives (1) the rotational average of Minkowski tensors defined on linear subspaces and (2) the functional defined on linear subspaces with rotational average equal to a Minkowski tensor. Earlier results obtained for intrinsic volumes appear now as special cases.

Proton diffusion tensor spectroscopy of metabolites in human muscle in vivo

PURPOSE To study the apparent diffusivity and its directionality for metabolites of skeletal muscle in humans in vivo by (1) H magnetic resonance spectroscopy. METHODS The diffusion tensors were determined on a 3 Tesla MR system using optimized acquisition and processing methods including an adapted STEAM sequence with orientation-dependent diffusion weighting, pulse-triggering with individually adapted delays, eddy-current correction schemes, median filtering, and simultaneous prior-knowledge fitting of all related spectra. RESULTS The average apparent diffusivities, as well as the fractional anisotropies of taurine (ADCav  = 0.74 × 10(-3) s/mm(2) , FA = 0.46), creatine (ADCav  = 0.41 × 10(-3)  s/mm(2) , FA = 0.33), trimethylammonium compounds (ADCav  = 0.48 × 10(-3)  s/mm(2) , FA = 0.34), carnosine (ADCav  = 0.46 × 10(-3)  s/mm(2) , FA = 0.47), and water (ADCav  = 1.5 × 10(-3)  s/mm(2) , FA = 0.36) were estimated. The diffusivities of most metabolites and water were significantly different from each other. Diffusion was found to be anisotropic and the diffusion tensors showed tensor correlation coefficients close to 1 and were hence found to be essentially coaligned. The magnitudes of apparent metabolite diffusivities were largely ordered according to molecular weight, with taurine as the smallest molecule diffusing fastest, both along and across the fiber direction. CONCLUSION Diffusivities, directional dependence of diffusion and fractional anisotropies of (1) H MRS-visible muscle metabolites were presented. It was shown that metabolites share diffusion directionality with water and have similar fractional anisotropies, hinting at similar diffusion barriers. Magn Reson Med, 2014. © 2014 Wiley Periodicals, Inc.

Measurement of structural anisotropy in femoral trabecular bone using clinical-resolution CT images

Discrepancies in finite-element model predictions of bone strength may be attributed to the simplified modeling of bone as an isotropic structure due to the resolution limitations of clinical-level Computed Tomography (CT) data. The aim of this study is to calculate the preferential orientations of bone (the principal directions) and the extent to which bone is deposited more in one direction compared to another (degree of anisotropy). Using 100 femoral trabecular samples, the principal directions and degree of anisotropy were calculated with a Gradient Structure Tensor (GST) and a Sobel Structure Tensor (SST) using clinical-level CT. The results were compared against those calculated with the gold standard Mean-Intercept-Length (MIL) fabric tensor using micro-CT. There was no significant difference between the GST and SST in the calculation of the main principal direction (median error=28°), and the error was inversely correlated to the degree of transverse isotropy (r=−0.34, p<0.01). The degree of anisotropy measured using the structure tensors was weakly correlated with the MIL-based measurements (r=0.2, p<0.001). Combining the principal directions with the degree of anisotropy resulted in a significant increase in the correlation of the tensor distributions (r=0.79, p<0.001). Both structure tensors were robust against simulated noise, kernel sizes, and bone volume fraction. We recommend the use of the GST because of its computational efficiency and ease of implementation. This methodology has the promise to predict the structural anisotropy of bone in areas with a high degree of anisotropy, and may improve the in vivo characterization of bone.

Calculation of crystal optical properties from molecular electron density

We have recently developed a method to obtain distributed atomic polarizabilities adopting a partitioning of the molecular electron density (for example, the Quantum Theory of Atoms in Molecules, [1]), calculated with or without an applied electric field. The procedure [2] allows to obtained atomic polarizability tensors, which are perfectly exportable, because quite representative of an atom in a given functional group. Among the many applications of this idea, the calculation of crystal susceptibility is easily available, either from a rough estimation (the polarizability of the isolated molecule is used) or from a more precise estimation (the polarizability of a molecule embedded in a cluster representing the first coordination sphere is used). Lorentz factor is applied to include the long range effect of packing, which is enhancing the molecular polarizability. Simple properties like linear refractive index or the gyration tensor can be calculated at relatively low costs and with good precision. This approach is particularly useful within the field of crystal engineering of organic/organometallic materials, because it would allow a relatively easy prediction of a property as a function of the packing, thus allowing "reverse crystal engineering". Examples of some amino acid crystals and salts of amino acids [3] will be illustrated, together with other crystallographic or non-crystallographic applications. For example, the induction and dispersion energies of intermolecular interactions could be calculated with superior precision (allowing anisotropic van der Waals interactions). This could allow revision of some commonly misunderstood intermolecular interactions, like the halogen bonding (see for example the recent remarks by Stone or Gilli [4]). Moreover, the chemical reactivity of coordination complexes could be reinvestigated, by coupling the conventional analysis of the electrostatic potential (useful only in the circumstances of hard nucleophilic/electrophilic interaction) with the distributed atomic polarizability. The enhanced reactivity of coordinated organic ligands would be better appreciated. [1] R. F. W. Bader, Atoms in Molecules: A Quantum Theory. Oxford Univ. Press, 1990. [2] A. Krawczuk-Pantula, D. Pérez, K. Stadnicka, P. Macchi, Trans. Amer. Cryst. Ass. 2011, 1-25 [3] A. S. Chimpri1, M. Gryl, L. H.R. Dos Santos1, A. Krawczuk, P. Macchi Crystal Growth & Design, in the press. [4] a) A. J. Stone, J. Am. Chem. Soc. 2013, 135, 7005−7009; b) V. Bertolasi, P. Gilli, G. Gilli Crystal Growth & Design, 2013, 12, 4758-4770.

Petrov classification and holographic reconstruction of spacetime

Using the asymptotic form of the bulk Weyl tensor, we present an explicit approach that allows us to reconstruct exact four-dimensional Einstein spacetimes which are algebraically special with respect to Petrov’s classification. If the boundary metric supports a traceless, symmetric and conserved complex rank-two tensor, which is related to the boundary Cotton and energy-momentum tensors, and if the hydrodynamic congruence is shearless, then the bulk metric is exactly resummed and captures modes that stand beyond the hydrodynamic derivative expansion. We illustrate the method when the congruence has zero vorticity, leading to the Robinson-Trautman spacetimes of arbitrary Petrov class, and quote the case of non-vanishing vorticity, which captures the Plebański-Demiański Petrov D family.