Electromagnetic fields inside a lossy, multilayered spherical head phantom excited by MRI coils: models and methods

The precise evaluation of electromagnetic field (EMF) distributions inside biological samples is becoming an increasingly important design requirement for high field MRI systems. In evaluating the induced fields caused by magnetic field gradients and RF transmitter coils, a multilayered dielectric spherical head model is proposed to provide a better understanding of electromagnetic interactions when compared to a traditional homogeneous head phantom. This paper presents Debye potential (DP) and Dyadic Green's function (DGF)-based solutions of the EMFs inside a head-sized, stratified sphere with similar radial conductivity and permittivity profiles as a human head. The DP approach is formulated for the symmetric case in which the source is a circular loop carrying a harmonic-formed current over a wide frequency range. The DGF method is developed for generic cases in which the source may be any kind of RF coil whose current distribution can be evaluated using the method of moments. The calculated EMFs can then be used to deduce MRI imaging parameters. The proposed methods, while not representing the full complexity of a head model, offer advantages in rapid prototyping as the computation times are much lower than a full finite difference time domain calculation using a complex head model. Test examples demonstrate the capability of the proposed models/methods. It is anticipated that this model will be of particular value for high field MRI applications, especially the rapid evaluation of RF resonator (surface and volume coils) and high performance gradient set designs.

Neuro-electromagnetic imaging of the human visual cortex

Methods of solving the neuro-electromagnetic inverse problem are examined and developed, with specific reference to the human visual cortex. The anatomy, physiology and function of the human visual system are first reviewed. Mechanisms by which the visual cortex gives rise to external electric and magnetic fields are then discussed, and the forward problem is described mathematically for the case of an isotropic, piecewise homogeneous volume conductor, and then for an anisotropic, concentric, spherical volume conductor. Methods of solving the inverse problem are reviewed, before a new technique is presented. This technique combines prior anatomical information gained from stereotaxic studies, with a probabilistic distributed-source algorithm to yield accurate, realistic inverse solutions. The solution accuracy is enhanced by using both visual evoked electric and magnetic responses simultaneously. The numerical algorithm is then modified to perform equivalent current dipole fitting and minimum norm estimation, and these three techniques are implemented on a transputer array for fast computation. Due to the linear nature of the techniques, they can be executed on up to 22 transputers with close to linear speedup. The latter part of the thesis describes the application of the inverse methods to the analysis of visual evoked electric and magnetic responses. The CIIm peak of the pattern onset evoked magnetic response is deduced to be a product of current flowing away from the surface areas 17, 18 and 19, while the pattern reversal P100m response originates in the same areas, but from oppositely directed current. Cortical retinotopy is examined using sectorial stimuli, the CI and CIm ;peaks of the pattern onset electric and magnetic responses are found to originate from areas V1 and V2 simultaneously, and they therefore do not conform to a simple cruciform model of primary visual cortex.

Quantitative Analysis of Kilohertz-Frequency Neurostimulation

Mainstream electrical stimulation therapies, e.g., spinal cord stimulation (SCS) and deep brain stimulation, use pulse trains that are delivered at rates no higher than 200 Hz. In recent years, stimulation of nerve fibers using kilohertz-frequency (KHF) signals has received increased attention due to the potential to penetrate deeper in the tissue and to the ability to block conduction of action potentials. As well, there are a growing number of clinical applications that use KHF waveforms, including transcutaneous electrical stimulation (TES) for overactive bladder and SCS for chronic pain. However, there is a lack of fundamental understanding of the mechanisms of action of KHF stimulation. The goal of this research was to analyze quantitatively KHF neurostimulation.

We implemented a multilayer volume conductor model of TES including dispersion and capacitive effects, and we validated the model with in vitro measurements in a phantom constructed from dispersive materials. We quantified the effects of frequency on the distribution of potentials and fiber excitation. We also quantified the effects of a novel transdermal amplitude modulated signal (TAMS) consisting of a non-zero offset sinusoidal carrier modulated by a square-pulse train. The model revealed that high-frequency signals generated larger potentials at depth than did low frequencies, but this did not translate into lower stimulation thresholds. Both TAMS and conventional rectangular pulses activated more superficial fibers in addition to the deeper, target fibers, and at no frequency did we observe an inversion of the strength-distance relationship. In addition, we performed in vivo experiments and applied direct stimulation to the sciatic nerve of cats and rats. We measured electromyogram and compound action potential activity evoked by pulses, TAMS and modified versions of TAMS in which we varied the amplitude of the carrier. Nerve fiber activation using TAMS showed no difference with respect to activation with conventional pulse for carrier frequencies of 20 kHz and higher, regardless the size of the carrier. Therefore, TAMS with carrier frequencies >20 kHz does not offer any advantage over conventional pulses, even with larger amplitudes of the carrier, and this has implications for design of waveforms for efficient and effective TES.

We developed a double cable model of a dorsal column (DC) fiber to quantify the responses of DC fibers to a novel KHF-SCS signal. We validated the model using in vivo recordings of the strength-duration relationship and the recovery cycle of single DC fibers. We coupled the fiber model to a model of SCS in human and applied the KHF-SCS signal to quantify thresholds for activation and conduction block for different fiber diameters at different locations in the DCs. Activation and block thresholds increased sharply as the fibers were placed deeper in the DCs, and decreased for larger diameter fibers. Activation thresholds were > 5 mA in all cases and up to five times higher than for conventional (~ 50 Hz) SCS. For fibers exhibiting persistent activation, the degree of synchronization of the firing activity to the KHF-SCS signal, as quantified using the vector strength, was low for a broad amplitude range, and the dissimilarity between the activities in pairs of fibers, as quantified using the spike time distance, was high and decreased for more closely positioned fibers. Conduction block thresholds were higher than 30 mA for all fiber diameters at any depth and well above the amplitudes used clinically (0.5 – 5 mA). KHF-SCS appears to activate few, large, superficial fibers, and the activated fibers fire asynchronously to the stimulation signal and to other activated fibers.

The outcomes of this work contribute to the understanding of KHF neurostimulation by establishing the importance of the tissue filtering properties on the distribution of potentials, assessing quantitatively the impact of KHF stimulation on nerve fiber excitation, and developing and validating a detailed model of a DC fiber to characterize the effects of KHF stimulation on DC axons. The results have implications for design of waveforms for efficient and effective nerve fiber stimulation in the peripheral and central nervous system.

High-Pressure Nmr And Compressibility Evidence For A Phase-Transition In The Protonic Conductor (NH4)4Fe(CN)6.1.5H2o

1H NMR at high hydrostatic pressures and compressibility studies show that the protonic conductor (NH4)4Fe(CN)6·1.5H2O undergoes a phase transition around 0.45 GPa. The transition is characterized by a large hysteresis. From the NMR studies, an activation volume of 6% is obtained below the phase transition, indicating the dominance of Frenkel defects.

Studies on dielectric behavior of insulator-conductor composites

This paper describes the dielectric behavior of an insulator-conductor composite, namely, the wax-graphite composite. The variation of specific capacitance of these composites with parameters such as volume fraction and grain size of the conducting particles and temperature has been studied. These observed variations have been explained using the same model [C. Rajagopal and M. Satyam, J. Appl. Phys. 49, 5536 (1978)] which explains electrical conduction in composites. The specific capacitance of these materials appears to be governed by the contact capacitance between the conducting particles and the number of contacts each particle has with its neighbors. The variation of specific capacitance with temperature is attributed to the change in contact area.

Studies on electrical conductivity of insulator-conductor composites

The variation of electrical resistivity of an insulator-conductor composite, namely, wax-graphite composite, with parameters such as volume fraction, grain size, and temperature has been studied. A model is proposed to explain the observed variations, which assumes that the texture of the composite consists of insulator granules coated with conducting particles. The resistivity of these materials is controlled mainly by the contact resistance between the conducting particles and the number of contacts each particle has with its neighbors. The variation of resistivity with temperature has also been explained with the help of this model and it is attributed to the change in contact area. Journal of Applied Physics is copyrighted by The American Institute of Physics.

High-pressure NMR investigations of the protonic conductor (NH4)4Fe(CN)6.1.5H2O

The effect of high hydrostatic pressure up to 1.5 GPa on ionic motion in (NH4)4Fe(CN)6.1.5H2O has been studied by wide-line 1H NMR experiments performed in the temperature range from room temperature to 77 K. The experiments at room temperature have shown a large increase in the second moment at 0.45 GPa as a result of a pressure-induced phase transition. The temperature dependence study up to 0.425 GPa has shown a gradual increase in the values of activation energy and attempt frequency with increase in pressure. The activation volume for motion at 300 K has been estimated to be 6% of molar volume. Vacancy-assisted ionic jumps are concluded to be the mode of charge transport. Second moments estimated at 77 K show evidence for tunnelling reorientation of at least one of the two NH4+ groups in the compound.

Studies on dielectric behavior of insulator‐conductor composites

This paper describes the dielectric behavior of an insulator‐conductor composite, namely, the wax‐graphite composite. The variation of specific capacitance of these composites with parameters such as volume fraction and grain size of the conducting particles and temperature has been studied. These observed variations have been explained using the same model [C. Rajagopal and M. Satyam, J. Appl. Phys. 49, 5536 (1978)] which explains electrical conduction in composites. The specific capacitance of these materials appears to be governed by the contact capacitance between the conducting particles and the number of contacts each particle has with its neighbors. The variation of specific capacitance with temperature is attributed to the change in contact area.

A second order control-volume finite-element least-squares strategy for simulating diffusion in strongly anisotropic media

An unstructured mesh �nite volume discretisation method for simulating di�usion in anisotropic media in two-dimensional space is discussed. This technique is considered as an extension of the fully implicit hybrid control-volume �nite-element method and it retains the local continuity of the ux at the control volume faces. A least squares function recon- struction technique together with a new ux decomposition strategy is used to obtain an accurate ux approximation at the control volume face, ensuring that the overall accuracy of the spatial discretisation maintains second order. This paper highlights that the new technique coincides with the traditional shape function technique when the correction term is neglected and that it signi�cantly increases the accuracy of the previous linear scheme on coarse meshes when applied to media that exhibit very strong to extreme anisotropy ratios. It is concluded that the method can be used on both regular and irregular meshes, and appears independent of the mesh quality.