Speeding up the EM algorithm for mixture model-based segmentation of magnetic resonance images

Mixture models implemented via the expectation-maximization (EM) algorithm are being increasingly used in a wide range of problems in pattern recognition such as image segmentation. However, the EM algorithm requires considerable computational time in its application to huge data sets such as a three-dimensional magnetic resonance (MR) image of over 10 million voxels. Recently, it was shown that a sparse, incremental version of the EM algorithm could improve its rate of convergence. In this paper, we show how this modified EM algorithm can be speeded up further by adopting a multiresolution kd-tree structure in performing the E-step. The proposed algorithm outperforms some other variants of the EM algorithm for segmenting MR images of the human brain. (C) 2004 Pattern Recognition Society. Published by Elsevier Ltd. All rights reserved.

Studies on the membrane interactions of the cyclotides kalata B1 and kalata B6 on model membrane systems by surface plasmon resonance

In this study we have demonstrated the interactions of kalata B1 and its naturally occurring analogue kalata B6 with five model lipid membranes and have analyzed the binding kinetics using surface plasmon resonance. Two kalata peptides showed a higher affinity for the phosphatidylethanolamine-containing membranes, indicating that the peptides would bind selectively to bacterial membranes. Also we have optimized the procedure for the immobilization of five liposome mixtures and have shown that the procedure provides reproducible levels of immobilized liposomes and could be used to screen the selective binding of putative antimicrobial peptides to model mammalian or microbial phospholipid membranes. (C) 2004 Elsevier Inc. All rights reserved.

Quantum transport and integrability of the Anderson model for a quantum dot with multiple leads

We show that an Anderson Hamiltonian describing a quantum dot connected to multiple leads is integrable. A general expression for the nonlinear conductance is obtained by combining the Bethe ansatz exact solution with Landauer-Buttiker theory. In the Kondo regime, a closed form expression is given for the matrix conductance at zero temperature and when all the leads are close to the symmetric point. A bias-induced splitting of the Kondo resonance is possible for three or more leads. Specifically, for N leads, with each at a different chemical potential, there can be N-1 Kondo peaks in the conductance.

Numerical modelling of thermal effects in rats due to high-field magnetic resonance imaging (0.5-1 GHz)

A finite-difference time-domain (FDTD) thermal model has been developed to compute the temperature elevation in the Sprague Dawley rat due to electromagnetic energy deposition in high-field magnetic resonance imaging (MRI). The field strengths examined ranged from 11.75-23.5 T (corresponding to H-1 resonances of 0.5-1 GHz) and an N-stub birdcage resonator was used to both transmit radio-frequency energy and receive the MRI signals. With an in-plane resolution of 1.95 mm, the inhomogeneous rat phantom forms a segmented model of 12 different tissue types, each having its electrical and thermal parameters assigned. The steady-state temperature distribution was calculated using a Pennes 'bioheat' approach. The numerical algorithm used to calculate the induced temperature distribution has been successfully validated against analytical solutions in the form of simplified spherical models with electrical and thermal properties of rat muscle. As well as assisting with the design of MRI experiments and apparatus, the numerical procedures developed in this study could help in future research and design of tumour-treating hyperthermia applicators to be used on rats in vivo.

An FDTD model for calculation of gradient-induced eddy currents in MRI system

In most magnetic resonance imaging (MRI) systems, pulsed magnetic gradient fields induce eddy currents in the conducting structures of the superconducting magnet. The eddy currents induced in structures within the cryostat are particularly problematic as they are characterized by long time constants by virtue of the low resistivity of the conductors. This paper presents a three-dimensional (3-D) finite-difference time-domain (FDTD) scheme in cylindrical coordinates for eddy-current calculation in conductors. This model is intended to be part of a complete FDTD model of an MRI system including all RF and low-frequency field generating units and electrical models of the patient. The singularity apparent in the governing equations is removed by using a series expansion method and the conductor-air boundary condition is handled using a variant of the surface impedance concept. The numerical difficulty due to the asymmetry of Maxwell equations for low-frequency eddy-current problems is circumvented by taking advantage of the known penetration behavior of the eddy-current fields. A perfectly matched layer absorbing boundary condition in 3-D cylindrical coordinates is also incorporated. The numerical method has been verified against analytical solutions for simple cases. Finally, the algorithm is illustrated by modeling a pulsed field gradient coil system within an MRI magnet system. The results demonstrate that the proposed FDTD scheme can be used to calculate large-scale eddy-current problems in materials with high conductivity at low frequencies.

Coherent molecular bound states of bosons and fermions near a Feshback resonance

We analyze molecular bound states of atomic quantum gases near a Feshbach resonance. A simple, renormalizable field theoretic model is shown to have exact solutions in the two-body sector, whose binding energy agrees well with observed experimental results in both Bosonic and Fermionic cases. These solutions, which interpolate between BEC and BCS theories, also provide a more general variational ansatz for resonant superfluidity and related problems.

Model for growth hormone receptor activation based on subunit rotation within a receptor dimer

Growth hormone is believed to activate the growth hormone receptor (GHR) by dimerizing two identical receptor subunits, leading to activation of JAK2 kinase associated with the cytoplasmic domain. However, we have reported previously that dimerization alone is insufficient to activate full-length GHR. By comparing the crystal structure of the liganded and unliganded human GHR extracellular domain, we show here that there is no substantial change in its conformation on ligand binding. However, the receptor can be activated by rotation without ligand by inserting a defined number of alanine residues within the transmembrane domain. Fluorescence resonance energy transfer ( FRET), bioluminescence resonance energy transfer (BRET) and coimmunoprecipitation studies suggest that receptor subunits undergo specific transmembrane interactions independent of hormone binding. We propose an activation mechanism involving a relative rotation of subunits within a dimeric receptor as a result of asymmetric placement of the receptor-binding sites on the ligand.

A neuro-mechanical model for interpersonal coordination

The present study investigates the coordination between two people oscillating handheld pendulums, with a special emphasis on the influence of the mechanical properties of the effector systems involved. The first part of the study is an experiment in which eight pairs of participants are asked to coordinate the oscillation of their pendulum with the other participant's in an in-phase or antiphase fashion. Two types of pendulums, A and B, having different resonance frequencies (Freq A=0.98 Hz and Freq B=0.64 Hz), were used in different experimental combinations. Results confirm that the preferred frequencies produced by participants while manipulating each pendulum individually were close to the resonance frequencies of the pendulums. In their attempt to synchronize with one another, participants met at common frequencies that were influenced by the mechanical properties of the two pendulums involved. In agreement with previous studies, both the variability of the behavior and the shift in the intended relative phase were found to depend on the task-effector asymmetry, i.e., the difference between the mechanical properties of the effector systems involved. In the second part of the study, we propose a model to account for these results. The model consists of two cross-coupled neuro-mechanical units, each composed of a neural oscillator driving a wrist-pendulum system. Taken individually, each unit reproduced the natural tendency of the participants to freely oscillate a pendulum close to its resonance frequency. When cross-coupled through the vision of the pendulum of the other unit, the two units entrain each other and meet at a common frequency influenced by the mechanical properties of the two pendulums involved. The ability of the proposed model to address the other effects observed as a function of the different conditions of the pendulum and intended mode of coordination is discussed.

High-field magnetic resonance imaging with reduced field/tissue RF artefacts - A modeling study using hybrid MoM/FEM and FDTD technique

Due to complex field/tissue interactions, high-field magnetic resonance (MR) images suffer significant image distortions that result in compromised diagnostic quality. A new method that attempts to remove these distortions is proposed in this paper and is based on the use of transceiver-phased arrays. The proposed system uses, in the examples presented herein, a shielded four-element transceive-phased array head coil and involves performing two separate scans of the same slice with each scan using different excitations during transmission. By optimizing the amplitudes and phases for each scan, antipodal signal profiles can be obtained, and by combining both the images together, the image distortion can be reduced several fold. A combined hybrid method of moments (MoM)/finite element method (FEM) and finite-difference time-domain (FDTD) technique is proposed and used to elucidate the concept of the new method and to accurately evaluate the electromagnetic field (EMF) in a human head model. In addition, the proposed method is used in conjunction with the generalized auto-calibrating partially parallel acquisitions (GRAPPA) reconstruction technique to enable rapid imaging of the two scans. Simulation results reported herein for 11-T (470-MHz) brain imaging applications show that the new method with GRAPPA reconstruction theoretically results in improved image quality and that the proposed combined hybrid MoM/FEM and FDTD technique is. suitable for high-field magnetic resonance imaging (MRI) numerical analysis.

Determinants of visual 2-D echo assessment of left ventricular wall motion: Comparison with myocardial thickening on cardiac magnetic resonance imaging

Background: left ventricular wall motion on 2d echo (2de) is usually scored visually. we sought to examine the determinants of visually assessed wall motion scoring on 2de by comparison with myocardial thickening quantified on MRI. Methods: using a 16 segment model, we studied 287 segments in 30 patients aged 61+/ -11 years (6 female), with ischaemic LV dysfunction (defined by at least 2 segments dysfunctional on 2de). 2de was performed in 5 views and wall motion scores (WMS) assigned: 1 (normal) 103 segments, 2 (hypokinetic) 93 segments, 3 (akinetic) 87 segments. MRI was used to measure end systolic wall thickness (ESWT), end diastolic wall thickness (EDWT) and percentage systolic wall thickening (SWT%) in the plane of the 2de and to assess WMS in the same planes visually. No patient had a clinical ischemic event between the tests. Results: visual assessment of wall motion by 2de and MRI showed moderate agreement (kappa = 0.425). Resting 2de wall motion correlated significantly (p

XSophe, a computer simulation software suite for the analysis of electron paramagnetic resonance spectra

The XSophe computer simulation software suite consisting of a daemon, the XSophe interface and the computational program Sophe is a state of the art package for the simulation of electron paramagnetic resonance spectra. The Sophe program performs the computer simulation and includes a number of new technologies including; the SOPHE partition and interpolation schemes, a field segmentation algorithm, homotopy, parallelisation and spectral optimisation. The SOPHE partition and interpolation scheme along with a field segmentation algorithm greatly increases the speed of simulations for most systems. Multidimensional homotopy provides an efficient method for accurately tracing energy levels and hence tracing transitions in the presence of energy level anticrossings and looping transitions and allowing computer simulations in frequency space. Recent enhancements to Sophe include the generalised treatment of distributions of orientational parameters, termed the mosaic misorientation linewidth model and a faster more efficient algorithm for the calculation of resonant field positions and transition probabilities. For complex systems the parallelisation enables the simulation of these systems on a parallel computer and the optimisation algorithms in the suite provide the experimentalist with the possibility of finding the spin Hamiltonian parameters in a systematic manner rather than a trial-and-error process. The XSophe software suite has been used to simulate multifrequency EPR spectra (200 MHz to 6 00 GHz) from isolated spin systems (S > ~½) and coupled centres (Si, Sj _> I/2). Griffin, M.; Muys, A.; Noble, C.; Wang, D.; Eldershaw, C.; Gates, K.E.; Burrage, K.; Hanson, G.R."XSophe, a Computer Simulation Software Suite for the Analysis of Electron Paramagnetic Resonance Spectra", 1999, Mol. Phys. Rep., 26, 60-84.

Brain activity during automatic semantic priming revealed by event-related functional magnetic resonance imaging

Semantic priming occurs when a subject is faster in recognising a target word when it is preceded by a related word compared to an unrelated word. The effect is attributed to automatic or controlled processing mechanisms elicited by short or long interstimulus intervals (ISIs) between primes and targets. We employed event-related functional magnetic resonance imaging (fMRI) to investigate blood oxygen level dependent (BOLD) responses associated with automatic semantic priming using an experimental design identical to that used in standard behavioural priming tasks. Prime-target semantic strength was manipulated by using lexical ambiguity primes (e.g., bank) and target words related to dominant or subordinate meaning of the ambiguity. Subjects made speeded lexical decisions (word/nonword) on dominant related, subordinate related, and unrelated word pairs presented randomly with a short ISI. The major finding was a pattern of reduced activity in middle temporal and inferior prefrontal regions for dominant versus unrelated and subordinate versus unrelated comparisons, respectively. These findings are consistent with both a dual process model of semantic priming and recent repetition priming data that suggest that reductions in BOLD responses represent neural priming associated with automatic semantic activation and implicate the left middle temporal cortex and inferior prefrontal cortex in more automatic aspects of semantic processing.