Best individual template selection from deformation tensor minimization

We study the influence of the choice of template in tensor-based morphometry. Using 3D brain MR images from 10 monozygotic twin pairs, we defined a tensor-based distance in the log-Euclidean framework [1] between each image pair in the study. Relative to this metric, twin pairs were found to be closer to each other on average than random pairings, consistent with evidence that brain structure is under strong genetic control. We also computed the intraclass correlation and associated permutation p-value at each voxel for the determinant of the Jacobian matrix of the transformation. The cumulative distribution function (cdf) of the p-values was found at each voxel for each of the templates and compared to the null distribution. Surprisingly, there was very little difference between CDFs of statistics computed from analyses using different templates. As the brain with least log-Euclidean deformation cost, the mean template defined here avoids the blurring caused by creating a synthetic image from a population, and when selected from a large population, avoids bias by being geometrically centered, in a metric that is sensitive enough to anatomical similarity that it can even detect genetic affinity among anatomies.

A novel measure of fractional anisotropy based on the tensor distribution function

Fractional anisotropy (FA), a very widely used measure of fiber integrity based on diffusion tensor imaging (DTI), is a problematic concept as it is influenced by several quantities including the number of dominant fiber directions within each voxel, each fiber's anisotropy, and partial volume effects from neighboring gray matter. With High-angular resolution diffusion imaging (HARDI) and the tensor distribution function (TDF), one can reconstruct multiple underlying fibers per voxel and their individual anisotropy measures by representing the diffusion profile as a probabilistic mixture of tensors. We found that FA, when compared with TDF-derived anisotropy measures, correlates poorly with individual fiber anisotropy, and may sub-optimally detect disease processes that affect myelination. By contrast, mean diffusivity (MD) as defined in standard DTI appears to be more accurate. Overall, we argue that novel measures derived from the TDF approach may yield more sensitive and accurate information than DTI-derived measures.

Analyzing multi-fiber reconstruction in high angular resolution diffusion imaging using the tensor distribution function

High-angular resolution diffusion imaging (HARDI) can reconstruct fiber pathways in the brain with extraordinary detail, identifying anatomical features and connections not seen with conventional MRI. HARDI overcomes several limitations of standard diffusion tensor imaging, which fails to model diffusion correctly in regions where fibers cross or mix. As HARDI can accurately resolve sharp signal peaks in angular space where fibers cross, we studied how many gradients are required in practice to compute accurate orientation density functions, to better understand the tradeoff between longer scanning times and more angular precision. We computed orientation density functions analytically from tensor distribution functions (TDFs) which model the HARDI signal at each point as a unit-mass probability density on the 6D manifold of symmetric positive definite tensors. In simulated two-fiber systems with varying Rician noise, we assessed how many diffusionsensitized gradients were sufficient to (1) accurately resolve the diffusion profile, and (2) measure the exponential isotropy (EI), a TDF-derived measure of fiber integrity that exploits the full multidirectional HARDI signal. At lower SNR, the reconstruction accuracy, measured using the Kullback-Leibler divergence, rapidly increased with additional gradients, and EI estimation accuracy plateaued at around 70 gradients.

The tensor distribution function

Diffusion weighted magnetic resonance imaging is a powerful tool that can be employed to study white matter microstructure by examining the 3D displacement profile of water molecules in brain tissue. By applying diffusion-sensitized gradients along a minimum of six directions, second-order tensors (represented by three-by-three positive definite matrices) can be computed to model dominant diffusion processes. However, conventional DTI is not sufficient to resolve more complicated white matter configurations, e.g., crossing fiber tracts. Recently, a number of high-angular resolution schemes with more than six gradient directions have been employed to address this issue. In this article, we introduce the tensor distribution function (TDF), a probability function defined on the space of symmetric positive definite matrices. Using the calculus of variations, we solve the TDF that optimally describes the observed data. Here, fiber crossing is modeled as an ensemble of Gaussian diffusion processes with weights specified by the TDF. Once this optimal TDF is determined, the orientation distribution function (ODF) can easily be computed by analytic integration of the resulting displacement probability function. Moreover, a tensor orientation distribution function (TOD) may also be derived from the TDF, allowing for the estimation of principal fiber directions and their corresponding eigenvalues.

Conformal invariance and scalar-tensor theories

A new generalisation of Einstein’s theory is proposed which is invariant under conformal mappings. Two scalar fields are introduced in addition to the metric tensor field, so that two special choices of gauge are available for physical interpretation, the ‘Einstein gauge’ and the ‘atomic gauge’. The theory is not unique but contains two adjustable parameters ζ anda. Witha=1 the theory viewed from the atomic gauge is Brans-Dicke theory (ω=−3/2+ζ/4). Any other choice ofa leads to a creation-field theory. In particular the theory given by the choicea=−3 possesses a cosmological solution satisfying Dirac’s ‘large numbers’ hypothesis.

Unitary Processes with Independent Increments and Representations of Hilbert Tensor Algebras

The aim of this article is to characterize unitary increment process by a quantum stochastic integral representation on symmetric Fock space. Under certain assumptions we have proved its unitary equivalence to a Hudson-Parthasarathy flow.

R1STM: One-class support tensor machine with randomised kernel

Identifying unusual or anomalous patterns in an underlying dataset is an important but challenging task in many applications. The focus of the unsupervised anomaly detection literature has mostly been on vectorised data. However, many applications are more naturally described using higher-order tensor representations. Approaches that vectorise tensorial data can destroy the structural information encoded in the high-dimensional space, and lead to the problem of the curse of dimensionality. In this paper we present the first unsupervised tensorial anomaly detection method, along with a randomised version of our method. Our anomaly detection method, the One-class Support Tensor Machine (1STM), is a generalisation of conventional one-class Support Vector Machines to higher-order spaces. 1STM preserves the multiway structure of tensor data, while achieving significant improvement in accuracy and efficiency over conventional vectorised methods. We then leverage the theory of nonlinear random projections to propose the Randomised 1STM (R1STM). Our empirical analysis on several real and synthetic datasets shows that our R1STM algorithm delivers comparable or better accuracy to a state-of-the-art deep learning method and traditional kernelised approaches for anomaly detection, while being approximately 100 times faster in training and testing.

Top-spin analysis of new scalar and tensor interactions in e(+)e(-) collisions with beam polarization

We utilize top polarization in the process e(+)e(-) -> t (t) over bar at the International Linear Collider ( ILC) with transverse beam polarization to probe interactions of the scalar and tensor type beyond the standard model and to disentangle their individual contributions. Ninety percent confidence level limits on the interactions with realistic integrated luminosity are presented and are found to improve by an order of magnitude compared to the case when the spin of the top quark is not measured. Sensitivities of the order of a few times 10(-3) TeV-2 for real and imaginary parts of both scalar and tensor couplings at root s = 500 and 800 GeV with an integrated luminosity of 500 fb(-1) and completely polarized beams are shown to be possible. A powerful model-independent framework for inclusive measurements is employed to describe the spin-momentum correlations, and their C, P, and T properties are presented in a technical appendix.

Top-spin analysis of new scalar and tensor interactions in e(+) e (-) collisions with transverse beam polarization

The top polarization at the International Linear Collider (ILC) with transverse beam polarization is utilized in the process to probe interactions of the scalar and tensor type beyond the Standard Model and to disentangle their individual contributions. Confidence level limits of 90% are presented on the interactions with realistic integrated luminosity and are found to improve by an order of magnitude compared to the case when the spin of the top quark is not measured. Sensitivities of the order of a few times 10 (-aEuro parts per thousand 3) TeV (-aEuro parts per thousand 2) for real and imaginary parts of both scalar and tensor couplings at and 800 GeV with an integrated luminosity of 500 fb (-aEuro parts per thousand 1) and completely polarized beams are shown to be possible.

Effect of compressibility on the velocity gradient tensor at free shear flow boundaries

Analyses of the invariants of the velocity gradient ten- sor were performed on flow fields obtained by DNS of compressible plane mixing layers at convective Mach num- bers Mc=0:15 and 1.1. Joint pdfs of the 2nd and 3rd invariants were examined at turbulent/nonturbulent (T/NT) boundaries—defined as surfaces where the local vorticity first exceeds a threshold fraction of the maximum of the mean vorticity. By increasing the threshold from very small lev-els, the boundary points were moved closer into the turbulent region, and the effects on the pdfs of the invariants were ob-served. Generally, T/NT boundaries are in sheet-like regions at both Mach numbers. At the higher Mach number a distinct lobe appears in the joint pdf isolines which has not been ob-served/reported before. A connection to the delayed entrain-ment and reduced growth rate of the higher Mach number flow is proposed.

Entanglement entropy from the holographic stress tensor

We consider entanglement entropy in the context of gauge/gravity duality for conformal field theories in even dimensions. The holographic prescription due to Ryu and Takayanagi (RT) leads to an equation describing how the entangling surface extends into the bulk geometry. We show that setting to zero, the timetime component of the Brown-York stress tensor evaluated on the co-dimension 1 entangling surface, leads to the same equation. By considering a spherical entangling surface as an example, we observe that the Euclidean actionmethods in AdS/CFT will lead to the RT area functional arising as a counterterm needed to regularize the stress tensor. We present arguments leading to a justification for the minimal area prescription.

Holographic stress tensor at finite coupling

We calculate one, two and three point functions of the holographic stress tensor for any bulk Lagrangian of the form L (g(ab), R-abcd, del(e) R-abcd). Using the first law of entanglement, a simple method has recently been proposed to compute the holographic stress tensor arising from a higher derivative gravity dual. The stress tensor is proportional to a dimension dependent factor which depends on the higher derivative couplings. In this paper, we identify this proportionality constant with a B-type trace anomaly in even dimensions for any bulk Lagrangian of the above form. This in turn relates to C-T, the coefficient appearing in the two point function of stress tensors. We use a background field method to compute the two and three point function of stress tensors for any bulk Lagrangian of the above form in arbitrary dimensions. As an application we consider general situations where eta/s for holographic plasmas is less than the KSS bound.

Approximation of Spatio-Temporal Random Processes Using Tensor Decomposition

A new representation of spatio-temporal random processes is proposed in this work. In practical applications, such processes are used to model velocity fields, temperature distributions, response of vibrating systems, to name a few. Finding an efficient representation for any random process leads to encapsulation of information which makes it more convenient for a practical implementations, for instance, in a computational mechanics problem. For a single-parameter process such as spatial or temporal process, the eigenvalue decomposition of the covariance matrix leads to the well-known Karhunen-Loeve (KL) decomposition. However, for multiparameter processes such as a spatio-temporal process, the covariance function itself can be defined in multiple ways. Here the process is assumed to be measured at a finite set of spatial locations and a finite number of time instants. Then the spatial covariance matrix at different time instants are considered to define the covariance of the process. This set of square, symmetric, positive semi-definite matrices is then represented as a third-order tensor. A suitable decomposition of this tensor can identify the dominant components of the process, and these components are then used to define a closed-form representation of the process. The procedure is analogous to the KL decomposition for a single-parameter process, however, the decompositions and interpretations vary significantly. The tensor decompositions are successfully applied on (i) a heat conduction problem, (ii) a vibration problem, and (iii) a covariance function taken from the literature that was fitted to model a measured wind velocity data. It is observed that the proposed representation provides an efficient approximation to some processes. Furthermore, a comparison with KL decomposition showed that the proposed method is computationally cheaper than the KL, both in terms of computer memory and execution time.

On the stability of the L-p-norm of the Riemannian curvature tensor

We consider the Riemannian functional defined on the space of Riemannian metrics with unit volume on a closed smooth manifold M where R(g) and dv (g) denote the corresponding Riemannian curvature tensor and volume form and p a (0, a). First we prove that the Riemannian metrics with non-zero constant sectional curvature are strictly stable for for certain values of p. Then we conclude that they are strict local minimizers for for those values of p. Finally generalizing this result we prove that product of space forms of same type and dimension are strict local minimizer for for certain values of p.