Multivariate variance-components analysis in DTI

Twin studies are a major research direction in imaging genetics, a new field, which combines algorithms from quantitative genetics and neuroimaging to assess genetic effects on the brain. In twin imaging studies, it is common to estimate the intraclass correlation (ICC), which measures the resemblance between twin pairs for a given phenotype. In this paper, we extend the commonly used Pearson correlation to a more appropriate definition, which uses restricted maximum likelihood methods (REML). We computed proportion of phenotypic variance due to additive (A) genetic factors, common (C) and unique (E) environmental factors using a new definition of the variance components in the diffusion tensor-valued signals. We applied our analysis to a dataset of Diffusion Tensor Images (DTI) from 25 identical and 25 fraternal twin pairs. Differences between the REML and Pearson estimators were plotted for different sample sizes, showing that the REML approach avoids severe biases when samples are smaller. Measures of genetic effects were computed for scalar and multivariate diffusion tensor derived measures including the geodesic anisotropy (tGA) and the full diffusion tensors (DT), revealing voxel-wise genetic contributions to brain fiber microstructure.

Probabilistic multi-tensor estimation using the tensor distribution function

Diffusion weighted magnetic resonance (MR) 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 6 directions, second-order tensors can be computed to model dominant diffusion processes. However, conventional DTI is not sufficient to resolve crossing fiber tracts. Recently, a number of high-angular resolution schemes with greater than 6 gradient directions have been employed to address this issue. In this paper, we introduce the Tensor Distribution Function (TDF), a probability function defined on the space of symmetric positive definite matrices. 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 diffusion orientation distribution function (ODF) can easily be computed by analytic integration of the resulting displacement probability function.

White matter integrity measured by fractional anisotropy correlates poorly with actual individual fiber anisotropy

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. High-angular resolution diffusion imaging (HARDI) can resolve more complex diffusion geometries than standard DTI, including fibers crossing or mixing. The tensor distribution function (TDF) can be used to reconstruct multiple underlying fibers per voxel, representing the diffusion profile as a probabilistic mixture of tensors. Here we found that DTIderived mean diffusivity (MD) correlates well with actual individual fiber MD, but DTI-derived FA correlates poorly with actual individual fiber anisotropy, and may be suboptimal when used to detect disease processes that affect myelination. Analysis of the TDFs revealed that almost 40% of voxels in the white matter had more than one dominant fiber present. To more accurately assess fiber integrity in these cases, we here propose the differential diffusivity (DD), which measures the average anisotropy based on all dominant directions in each voxel.

How does angular resolution affect diffusion imaging measures?

A key question in diffusion imaging is how many diffusion-weighted images suffice to provide adequate signal-to-noise ratio (SNR) for studies of fiber integrity. Motion, physiological effects, and scan duration all affect the achievable SNR in real brain images, making theoretical studies and simulations only partially useful. We therefore scanned 50 healthy adults with 105-gradient high-angular resolution diffusion imaging (HARDI) at 4T. From gradient image subsets of varying size (6 ≤ N ≤ 94) that optimized a spherical angular distribution energy, we created SNR plots (versus gradient numbers) for seven common diffusion anisotropy indices: fractional and relative anisotropy (FA, RA), mean diffusivity (MD), volume ratio (VR), geodesic anisotropy (GA), its hyperbolic tangent (tGA), and generalized fractional anisotropy (GFA). SNR, defined in a region of interest in the corpus callosum, was near-maximal with 58, 66, and 62 gradients for MD, FA, and RA, respectively, and with about 55 gradients for GA and tGA. For VR and GFA, SNR increased rapidly with more gradients. SNR was optimized when the ratio of diffusion-sensitized to non-sensitized images was 9.13 for GA and tGA, 10.57 for FA, 9.17 for RA, and 26 for MD and VR. In orientation density functions modeling the HARDI signal as a continuous mixture of tensors, the diffusion profile reconstruction accuracy rose rapidly with additional gradients. These plots may help in making trade-off decisions when designing diffusion imaging protocols.

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.

Anisotropic Gaussian Schell-model beams: Passage through optical systems and associated invariants

Anisotropic Gaussian Schell-model (AGSM) fields and their transformation by first-order optical systems (FOS’s) forming Sp(4,R) are studied using the generalized pencils of rays. The fact that Sp(4,R), rather than the larger group SL(4,R), is the relevant group is emphasized. A convenient geometrical picture wherein AGSM fields and FOS’s are represented, respectively, by antisymmetric second-rank tensors and de Sitter transformations in a (3+2)-dimensional space is developed. These fields are shown to separate into two qualitatively different families of orbits and the invariants over each orbit, two in number, are worked out. We also develop another geometrical picture in a (2+1)-dimensional Minkowski space suitable for the description of the action of axially symmetric FOS’s on AGSM fields, and the invariants, now seven in number, are derived. Interesting limiting cases forming coherent and quasihomogeneous fields are analyzed.

ESR studies of phase transitions in double propionates

Detailed ESR investigations of Mn2+ substituting for Ca2+ in Ca2Sr(C2H5COO)6 (DSP), Ca2Pb(C2H5COO)6 (DLP) and Ca2Ba(C2H5COO)6 (DBP), in single crystals and powders, over the temperature range from 200°C to -180°C have been carried out to study the successive phase transitions in these compounds. (DSP: [Tetragonal] ← 8.5°C → [tetragonal, ferroelectric] [tetragonal] ← -169°C → [monoclinic, ferroelectric]; DLP : [tetragonal] ← 60°C → [tetragonal, ferroelectric] ← -71.5°C → [monoclinic, ferroelectric]; [Cubic] ← -6°C → [orthorhombic] ← -75°C → [?]). Spectra have been analysed in terms of axial spin Hamiltonians and the temperature dependences of the parameters studied. In DSP and DLP across the I ↔ II transition, new physically and chemically inequivalent sites appear indicating the disappearance of the diad axes on which the propionate groups are located, bringing out the connection between the motional states of the propionate groups and the occurence of ferroelectricity. The II ↔ III transition also causes chemically inequivalent sites to develop, indicating that the transitions may not be isomorphous as believed previously. In DBP, the -6°C transition leads to (i) a doubling of both physically and chemically inequivalent sites (ii) a small (150 G at -6°C to 170 G at -8°C), but abrupt change in the magnitude of the zero-field splitting tensor D, and (iii) displacements of the orientations of the D tensors. Results are interpreted in terms of alternate rotations of the oxygen octahedra, showing participation of the carboxyl oxygens in the transition. No drastic changes in the parameters occur across the -75°C transition consistent with its second order nature. Similarities and dissimilarities of the ESR spectra of the three compounds in relation to the phase transitions, are discussed.

ESR study of copper diethyldithiophosphate

ESR investigations are reported in single crystals of copper diethyldithiophosphate, magnetically diluted with the corresponding diamagnetic nickel complex. The spectrum at normal gain shows hyperfine components from 63Cu, 65Cu, and 31P nuclei. At much higher gain, hyperfine interaction from 33S nuclei in the ligand is detected. The spin Hamiltonian parameters relating to copper show tetragonal symmetry. The measured parameters are g = 2.085, g =2.025, A63Cu = 149.6 × 10−4 cm−1, A65Cu = 160.8 × 10−4 cm−1, BCu = 32.5 × 10−4 cm−1 and QCu 5.5 × 10−4cm−1. The 31P interaction is isotropic with a coupling constant AP = 9.6 × 10−4 cm−1. Angular variation of the 33S lines shows two different hyperfine tensors indicating the presence of two chemically inequivalent Cu S bonds. The experimentally determined hyperfine constants are A =34.9×10−4 cm−1, B =26.1×10−4 cm−1, A =60.4×10−4 cm−1, B =55.5×10−4 cm−1. The hyperfine parameters show that the hybridization of the ligand orbitals is very sensitive to the symmetry around the ligand. The g values and Cu hyperfine parameters are not much affected by the distortions occurring in the ligand. The energies of the d-d transitions are determined by optical absorption measurements on Cu diethyldithiophosphate in solution. Using the spin Hamiltonian parameters together with optical absorption results, the MO parameters for the complex are calculated. It is found that in addition to the bond, the bonds are also strongly covalent. ©1973 The American Institute of Physics

ESR study of copper diethyldithiophosphate

ESR investigations are reported in single crystals of copper diethyldithiophosphate, magnetically diluted with the corresponding diamagnetic nickel complex. The spectrum at normal gain shows hyperfine components from 63Cu, 65Cu, and 31P nuclei. At much higher gain, hyperfine interaction from 33S nuclei in the ligand is detected. The spin Hamiltonian parameters relating to copper show tetragonal symmetry. The measured parameters are g|| = 2.085, g[perpendicular]=2.025, A63Cu = 149.6 × 10−4 cm−1, A65Cu = 160.8 × 10−4 cm−1, BCu = 32.5 × 10−4 cm−1 and QCu [infinity] 5.5 × 10−4cm−1. The 31P interaction is isotropic with a coupling constant AP = 9.6 × 10−4 cm−1. Angular variation of the 33S lines shows two different hyperfine tensors indicating the presence of two chemically inequivalent Cu[Single Bond]S bonds. The experimentally determined hyperfine constants are A ₁^s=34.9×10−4 cm−1, B ₁^s=26.1×10−4 cm−1, A ₂^s=60.4×10−4 cm−1, B₂^s=55.5×10−4 cm−1. The hyperfine parameters show that the hybridization of the ligand orbitals is very sensitive to the symmetry around the ligand. The g values and Cu hyperfine parameters are not much affected by the distortions occurring in the ligand. The energies of the d-d transitions are determined by optical absorption measurements on Cu diethyldithiophosphate in solution. Using the spin Hamiltonian parameters together with optical absorption results, the MO parameters for the complex are calculated. It is found that in addition to the sigma bond, the pi bonds are also strongly covalent. ©1973 The American Institute of Physics.

Multicarrier conduction and Boltzmann transport analysis of heavy hole mobility in HgCdTe near room temperature

Magnetotransport measurements in pulsed fields up to 15 T have been performed on mercury cadmium telluride (Hg1-xCdxTe, x similar to 0.2) bulk as well as liquid phase epitaxially grown samples to obtain the resistivity and conductivity tensors in the temperature range 220-300 K. Mobilities and densities of various carriers participating in conduction have been extracted using both conventional multicarrier fitting (MCF) and mobility spectrum analysis. The fits to experimental data, particularly at the highest magnetic fields, were substantially improved when MCF is applied to minimize errors simultaneously on both resistivity and conductivity tensors. The semiclassical Boltzmann transport equation has been solved without using adjustable parameters by incorporating the following scattering mechanisms to fit the mobility: ionized impurity, polar and nonpolar optical phonons, acoustic deformation potential, and alloy disorder. Compared to previous estimates based on the relaxation time approximation with outscattering only, polar optical scattering and ionized impurity scattering limited mobilities are shown to be larger due to the correct incorporation of the inscattering term taking into account the overlap integrals in the valence band.

Dynamics and control of buckling type devices using SMA wire integrated beam

An analysis and design study using Shape Memory Alloy (SMA) wire integrated beam and its buckling shape control are reported. The dynamical system performance is analyzed with a mathematical set-up involving nonlocal and rate sensitive kinetics of phase transformation in the SMA wire. A standard phenomenological constitutive model reported by Brinson (1993) is modified by considering certain consistency conditions in the material property tensors and by eliminating spurious singularity. Considering the inhomogeneity effects, a finite element model of the SMA wire is developed. Simulations are carried out to study the buckling shape control of a beam integrated with SMA wire.

Improved procedures for static and dynamic analyses of wrinkled membranes

An analysis of large deformations of flexible membrane structures within the tension field theory is considered. A modification-of the finite element procedure by Roddeman et al. (Roddeman, D. G., Drukker J., Oomens, C. W J., Janssen, J. D., 1987, ASME J. Appl. Mech. 54, pp. 884-892) is proposed to study the wrinkling behavior of a membrane element. The state of stress in the element is determined through a modified deformation gradient corresponding to a fictive nonwrinkled surface. The new model uses a continuously modified deformation gradient to capture the location orientation of wrinkles more precisely. It is argued that the fictive nonwrinkled surface may be looked upon as an everywhere-taut surface in the limit as the minor (tensile) principal stresses over the wrinkled portions go to zero. Accordingly, the modified deformation gradient is thought of as the limit of a sequence of everywhere-differentiable tensors. Under dynamic excitations, the governing equations are weakly projected to arrive at a system of nonlinear ordinary differential equations that is solved using different integration schemes. It is concluded that, implicit integrators work much better than explicit ones in the present context.

Cubic and ultimate groups of complete symmetry

It is shown that the systems of definite actions described by polar and axial tensors of the second rank and their combinations during the superposition of their elements of complete symmetry with the elements of complete symmetry of the "grey" cube, result in 11 cubic crystallographical groups of complete symmetry. There are 35 ultimate groups (i.e., the groups having the axes of symmetry of infinite order) in complete symmetry of finite figures. 14 out of these groups are ultimate groups of symmetry of polar and axial tensors of the second rank and 24 are new groups. All these 24 ultimate groups are conventional groups since they cannot be presented by certain finite figures possessing the axes of symmetry {Mathematical expression}. Geometrical interpretation for some of the groups of complete symmetry is given. The connection between complete symmetry and physical properties of the crystals (electrical, magnetic and optical) is shown.