Finite difference time domain (FDTD) method for modeling the effect of switched gradients on the human body in MRI

Numerical modeling of the eddy currents induced in the human body by the pulsed field gradients in MRI presents a difficult computational problem. It requires an efficient and accurate computational method for high spatial resolution analyses with a relatively low input frequency. In this article, a new technique is described which allows the finite difference time domain (FDTD) method to be efficiently applied over a very large frequency range, including low frequencies. This is not the case in conventional FDTD-based methods. A method of implementing streamline gradients in FDTD is presented, as well as comparative analyses which show that the correct source injection in the FDTD simulation plays a crucial rule in obtaining accurate solutions. In particular, making use of the derivative of the input source waveform is shown to provide distinct benefits in accuracy over direct source injection. In the method, no alterations to the properties of either the source or the transmission media are required. The method is essentially frequency independent and the source injection method has been verified against examples with analytical solutions. Results are presented showing the spatial distribution of gradient-induced electric fields and eddy currents in a complete body model.

Analysis of transient eddy currents in MRI using a cylindrical FDTD method

Most magnetic resonance imaging (MRI) spatial encoding techniques employ low-frequency pulsed magnetic field gradients that undesirably induce multiexponentially decaying eddy currents in nearby conducting structures of the MRI system. The eddy currents degrade the switching performance of the gradient system, distort the MRI image, and introduce thermal loads in the cryostat vessel and superconducting MRI components. Heating of superconducting magnets due to induced eddy currents is particularly problematic as it offsets the superconducting operating point, which can cause a system quench. A numerical characterization of transient eddy current effects is vital for their compensation/control and further advancement of the MRI technology as a whole. However, transient eddy current calculations are particularly computationally intensive. In large-scale problems, such as gradient switching in MRI, conventional finite-element method (FEM)-based routines impose very large computational loads during generation/solving of the system equations. Therefore, other computational alternatives need to be explored. This paper outlines a three-dimensional finite-difference time-domain (FDTD) method in cylindrical coordinates for the modeling of low-frequency transient eddy currents in MRI, as an extension to the recently proposed time-harmonic scheme. The weakly coupled Maxwell's equations are adapted to the low-frequency regime by downscaling the speed of light constant, which permits the use of larger FDTD time steps while maintaining the validity of the Courant-Friedrich-Levy stability condition. The principal hypothesis of this work is that the modified FDTD routine can be employed to analyze pulsed-gradient-induced, transient eddy currents in superconducting MRI system models. The hypothesis is supported through a verification of the numerical scheme on a canonical problem and by analyzing undesired temporal eddy current effects such as the B-0-shift caused by actively shielded symmetric/asymmetric transverse x-gradient head and unshielded z-gradient whole-body coils operating in proximity to a superconducting MRI magnet.

Theoretical investigation of convective instability in inclined and fluid-saturated three-dimensional fault zones.

The convective instability of pore-fluid flow in inclined and fluid-saturated three-dimensional fault zones has been theoretically investigated in this paper. Due to the consideration of the inclined three-dimensional fault zone with any values of the inclined angle, it is impossible to use the conventional linear stability analysis method for deriving the critical condition (i.e., the critical Rayleigh number) which can be used to investigate the convective instability of the pore-fluid flow in an inclined three-dimensional fault zone system. To overcome this mathematical difficulty, a combination of the variable separation method and the integration elimination method has been used to derive the characteristic equation, which depends on the Rayleigh number and the inclined angle of the inclined three-dimensional fault zone. Using this characteristic equation, the critical Rayleigh number of the system can be numerically found as a function of the inclined angle of the three-dimensional fault zone. For a vertically oriented three-dimensional fault zone system, the critical Rayleigh number of the system can be explicitly derived from the characteristic equation. Comparison of the resulting critical Rayleigh number of the system with that previously derived in a vertically oriented three-dimensional fault zone has demonstrated that the characteristic equation of the Rayleigh number is correct and useful for investigating the convective instability of pore-fluid flow in the inclined three-dimensional fault zone system. The related numerical results from this investigation have indicated that: (1) the convective pore-fluid flow may take place in the inclined three-dimensional fault zone; (2) if the height of the fault zone is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone stabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; (3) if the thickness of the stratum is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone destabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; and that (4) the shape of the inclined three-dimensional fault zone may affect the convective instability of pore-fluid flow in the system. (C) 2004 Published by Elsevier B.V.

Tailoring magnetic field gradient design to magnet cryostat geometry

Eddy currents induced within a magnetic resonance imaging (MRI) cryostat bore during pulsing of gradient coils can be applied constructively together with the gradient currents that generate them, to obtain good quality gradient uniformities within a specified imaging volume over time. This can be achieved by simultaneously optimizing the spatial distribution and temporal pre-emphasis of the gradient coil current, to account for the spatial and temporal variation of the secondary magnetic fields due to the induced eddy currents. This method allows the tailored design of gradient coil/magnet configurations and consequent engineering trade-offs. To compute the transient eddy currents within a realistic cryostat vessel, a low-frequency finite-difference time-domain (FDTD) method using total-field scattered-field (TFSF) scheme has been performed and validated

FDTD study of an UWB radar technique for breast tumor detection

A finite difference time domain (FDTD) method is applied to investigate capabilities of an ultra-wide band (UWB) radar system to detect a breast tumor. The first part of the investigations concerns FDTD simulations of a phantom formed by a plastic container with liquid and a small reflecting target. The second part focuses on a three-dimensional numerical breast model with a small tumor. FDTD simulations are carried out assuming a planar incident wave. Various time snap shots of the electromagnetic field are recorded to learn about the physical phenomenon of reflection and scattering in different layers of the phantom.

A high definition, finite difference time domain method

A high definition, finite difference time domain (HD-FDTD) method is presented in this paper. This new method allows the FDTD method to be efficiently applied over a very large frequency range including low frequencies, which are problematic for conventional FDTD methods. In the method, no alterations to the properties of either the source or the transmission media are required. The method is essentially frequency independent and has been verified against analytical solutions within the frequency range 50 Hz-1 GHz. As an example of the lower frequency range, the method has been applied to the problem of induced eddy currents in the human body resulting from the pulsed magnetic field gradients of an MRI system. The new method only requires approximately 0.3% of the source period to obtain an accurate solution. (C) 2003 Elsevier Science Inc. All rights reserved.

Full-field Fourier-domain optical coherence tomography

Full-field Fourier-domain optical coherence tomography (3F-OCT) is a full-field version of spectraldomain/swept-source optical coherence tomography. A set of two-dimensional Fourier holograms is recorded at discrete wavenumbers spanning the swept-source tuning range. The resultant three-dimensional data cube contains comprehensive information on the three-dimensional morphological layout of the sample that can be reconstructed in software via three-dimensional discrete Fourier-transform. This method of recording of the OCT signal confers signal-to-noise ratio improvement in comparison with "flying-spot" time-domain OCT. The spatial resolution of the 3F-OCT reconstructed image, however, is degraded due to the presence of a phase cross-term, whose origin and effects are addressed in this paper. We present theoretical and experimental study of imaging performance of 3F-OCT, with particular emphasis on elimination of the deleterious effects of the phase cross-term.

A method for the design of MRI radiofrequency coils based on triangular and pulse basis functions

This paper presents a numerical technique for the design of an RF coil for asymmetric magnetic resonance imaging (MRI) systems. The formulation is based on an inverse approach where the cylindrical surface currents are expressed in terms of a combination of sub-domain basis functions: triangular and pulse functions. With the homogeneous transverse magnetic field specified in a spherical region, a functional method is applied to obtain the unknown current coefficients. The current distribution is then transformed to a conductor pattern by use of a stream function technique. Preliminary MR images acquired using a prototype RF coil are presented and validate the design method. (C) 2002 Elsevier Science B.V. All rights reserved.

Finite-difference time-domain-based studies of MRI pulsed field gradient-induced eddy currents inside the human body

In modern magnetic resonance imaging (MRI), patients are exposed to strong, rapidly switching magnetic gradient fields that, in extreme cases, may be able to elicit nerve stimulation. This paper presents theoretical investigations into the spatial distribution of induced current inside human tissues caused by pulsed z-gradient fields. A variety of gradient waveforms have been studied. The simulations are based on a new, high-definition, finite-difference time-domain method and a realistic inhomogeneous 10-mm resolution human body model with appropriate tissue parameters. it was found that the eddy current densities are affected not only by the pulse sequences but by many parameters such as the position of the body inside the gradient set, the local biological material properties and the geometry of the body. The discussion contains a comparison of these results with previous results found in the literature. This study and the new methods presented herein will help to further investigate the biological effects caused by the switched gradient fields in a MRI scan. (C) 2002 Wiley Periodicals, Inc.

Expression, purification, crystallization, and NMR studies of the helicase interaction domain of Escherichia coli DnaG primase

in Escherichia coli, the DnaG primase is the RNA polymerase that synthesizes RNA primers at replication forks. It is composed of three domains, a small N-terminal zinc-binding domain, a larger central domain responsible for RNA synthesis, and a C-terminal domain comprising residues 434-581 [DnaG(434-581)] that interact with the hexameric DnaB helicase. Presumably because of this interaction, it had not been possible previously to express the C-terminal domain in a stably transformed E coli strain. This problem was overcome by expression of DnaG(434-581) under control of tandem bacteriophage gimel-promoters, and the protein was purified in yields of 4-6 mg/L of culture and studied by NMR. A TOCSY spectrum of a 2 mM solution of the protein at pH 7.0, indicated that its structured core comprises residues 444-579. This was consistent with sequence conservation among most-closely related primases. Linewidths in a NOESY spectrum of a 0.5 mM sample in 10 mM phosphate, pH 6.05, 0.1 M NaCl, recorded at 36 degreesC, indicated the protein to be monomeric. Crystals of selenomethionine-substituted DnaG(434-581) obtained by the hanging-drop vapor-diffusion method were body-centered tetragonal, space group I4(1)22, with unit cell parameters a = b 142.2 Angstrom, c = 192.1 Angstrom, and diffracted beyond 2.7 Angstrom resolution with synchrotron radiation. (C) 2003 Elsevier Inc. All rights reserved.

A comparison of low-storage strategies for spectral analysis in dissipative systems: filter diagonalisation in the Lanczos representation and harmonic inversion of the Chebychev-order-domain autocorrelation function

We compare the performance of two different low-storage filter diagonalisation (LSFD) strategies in the calculation of complex resonance energies of the HO2, radical. The first is carried out within a complex-symmetric Lanczos subspace representation [H. Zhang, S.C. Smith, Phys. Chem. Chem. Phys. 3 (2001) 2281]. The second involves harmonic inversion of a real autocorrelation function obtained via a damped Chebychev recursion [V.A. Mandelshtam, H.S. Taylor, J. Chem. Phys. 107 (1997) 6756]. We find that while the Chebychev approach has the advantage of utilizing real algebra in the time-consuming process of generating the vector recursion, the Lanczos, method (using complex vectors) requires fewer iterations, especially for low-energy part of the spectrum. The overall efficiency in calculating resonances for these two methods is comparable for this challenging system. (C) 2001 Elsevier Science B.V. All rights reserved.

Characterization of a glycine receptor domain that controls the binding and gating mechanisms of the beta-amino acid agonist, taurine

The beta -amino acid, taurine, is a full agonist of the human glycine receptor al subunit when recombinantly expressed in a mammalian (HEK293) cell line, but a partial agonist of the same receptor when expressed in Xenopus oocytes. Several residues in the Ala101-Thr112 domain have previously been identified as determinants of beta -amino acid binding and gating mechanisms in Xenopus oocyte-expressed receptors. The present study used the substituted cysteine accessibility method to investigate the role of this domain in controlling taurine-specific binding and gating mechanisms of glycine receptors recombinantly expressed in mammalian cells. Asn102 and Glu103 are identified as taurine and glycine binding sites, whereas Ala101 is eliminated as a possible binding site. The N102C mutation also abolished the antagonistic actions of taurine, indicating that this site does not discriminate between the putative agonist- and antagonist-bound conformations of beta -amino acids. The effects of mutations from Lys104-Thr112 indicate that the mechanism by which this domain controls beta -amino acid-specific binding and gating processes differs substantially depending on whether the receptor is expressed in mammalian cells or Xenopus oocytes. Thr112 is the only domain element in mammalian cell-expressed GlyRs which was demonstrated to discriminate between glycine and taurine.

A single beta subunit M2 domain residue controls the picrotoxin sensitivity of alpha beta heteromeric glycine receptor chloride channels

This study investigated the residues responsible for the reduced picrotoxin sensitivity of the alpha beta heteromeric glycine receptor relative to the alpha homomeric receptor. By analogy with structurally related receptors, the beta subunit M2 domain residues P278 and F282 were considered the most likely candidates for mediating this effect. These residues align with G254 and T258 of the alpha subunit. The T258A, T258C and T258F mutations dramatically reduced the picrotoxin sensitivity of the alpha homomeric receptor. Furthermore, the converse F282T mutation in the beta subunit increased the picrotoxin sensitivity of the alpha beta heteromeric receptor. The P278G mutation in the beta subunit did not affect the picrotoxin sensitivity of the alpha beta heteromer. Thus, a ring of five threonines at the M2 domain depth corresponding to alpha subunit T258 is specifically required for picrotoxin sensitivity. Mutations to alpha subunit T258 also profoundly influenced the apparent glycine affinity. A substituted cysteine accessibility analysis revealed that the T258C sidechain increases its pore exposure in the channel open state. This provides further evidence for an allosteric mechanism of picrotoxin inhibition, but renders it unlikely that picrotoxin las an allosterically acting 'competitive' antagonist) binds to this residue.