Novel Numerical Basis Sets for Electromagnetic Field Expansion in Arbitrary Inhomogeneous Objects.

We investigated how to construct low-order subspace basis sets to accurately represent electromagnetic fields generated within inhomogeneous arbitrary objects by radio-frequency sources external to a Huygen's surface. The basis generation relies on the singular value decomposition of Green's functions integro-differential operators which makes it feasible to derive a reduced-order yet stable model. We present a detailed study of the theoretical and numerical requisites for generating such basis, and show how it can be used to calculate performance limits in magnetic resonance imaging applications. Finally, we propose a novel numerical framework for the computation of characteristic modes of arbitrary inhomogeneous objects. We validated accuracy and convergence properties of the numerical basis against a complete analytical basis in the case of a uniform spherical object. We showed that the discretization of the Huygens's surface has a minimal effect on the accuracy of the calculations, which mainly depended on the electromagnetic solver resolution and order of approximation.

Magnetic-Resonance-Based Electrical Property Mapping Using Global Maxwell Tomography With an 8-Channel Head Coil at 7 Tesla: A Simulation Study.

OBJECTIVE: Global Maxwell Tomography (GMT) is a recently introduced volumetric technique for noninvasive estimation of electrical properties (EP) from magnetic resonance measurements. Previous work evaluated GMT using ideal radiofrequency (RF) excitations. The aim of this simulation study was to assess GMT performance with a realistic RF coil. METHODS: We designed a transmit-receive RF coil with 8 decoupled channels for 7T head imaging. We calculated the RF transmit field ( B1+) inside heterogeneous head models for different RF shimming approaches, and used them as input for GMT to reconstruct EP for all voxels. RESULTS: Coil tuning/decoupling remained relatively stable when the coil was loaded with different head models. Mean error in EP estimation changed from [Formula: see text] to [Formula: see text] and from [Formula: see text] to [Formula: see text] for relative permittivity and conductivity, respectively, when changing head model without re-tuning the coil. Results slightly improved when an SVD-based RF shimming algorithm was applied, in place of excitation with one coil at a time. Despite errors in EP, RF transmit field ( B1+) and absorbed power could be predicted with less than [Formula: see text] error over the entire head. GMT could accurately detect a numerically inserted tumor. CONCLUSION: This work demonstrates that GMT can reliably reconstruct EP in realistic simulated scenarios using a tailored 8-channel RF coil design at 7T. Future work will focus on construction of the coil and optimization of GMT's robustness to noise, to enable in-vivo GMT experiments. SIGNIFICANCE: GMT could provide accurate estimations of tissue EP, which could be used as biomarkers and could enable patient-specific estimation of RF power deposition, which is an unsolved problem for ultra-high-field magnetic resonance imaging.

A formalism to investigate the optimal transmit efficiency in radiofrequency shimming.

Transmit efficiency specifies the amplitude of the magnetic resonance excitation field produced over a region of interest with respect to the radiofrequency (RF) power deposited in the sample. This metric is highly important at ultra-high field magnetic resonance imaging (≥7 T), where excitation inhomogeneities and electric field interference effects could prevent achieving the desired flip angle distribution while satisfying the power safety limits. The aim of this work was to introduce an approach to calculate a theoretical upper bound on the transmit efficiency (OPTXE) for RF shimming, independent from any particular coil design. We computed the OPTXE for head-mimicking uniform spherical samples and a realistic heterogeneous head model by maximizing the square of the net transmit field per unit power deposition. The corresponding RF shimming weights were used to combine the analytical surface current modes into ideal current patterns. OPTXE grew monotonically as the target excitation voxel approached the surface of the object, and overall decreased at higher field strengths, presenting similar trends in both the uniform sphere and heterogeneous head model. Arrays with an increasing number of loops could closely approach OPTXE in the central region of the object, but performance decreased closer to the surface and at higher magnetic field strengths. The performance of 32 loops for a two-dimensional excitation region at 7 T increased from 34% to 93% when they were arranged based on the shape of the ideal current patterns. OPTXE provides an absolute reference to evaluate coil designs and RF shimming algorithms, whereas ideal current patterns could serve as guidelines for novel coil designs at ultra-high field. The uniform sphere model enables rapid analytic simulations and provides a good approximation of the OPTXE distribution in a realistic heterogeneous head model with comparable dimensions.

Noninvasive Estimation of Electrical Properties From Magnetic Resonance Measurements via Global Maxwell Tomography and Match Regularization.

OBJECTIVE: In this paper, we introduce global Maxwell tomography (GMT), a novel volumetric technique that estimates electric conductivity and permittivity by solving an inverse scattering problem based on magnetic resonance measurements. METHODS: GMT relies on a fast volume integral equation solver, MARIE, for the forward path, and a novel regularization method, match regularization, designed specifically for electrical property estimation from noisy measurements. We performed simulations with three different tissue-mimicking numerical phantoms of different complexity, using synthetic transmit sensitivity maps with realistic noise levels as the measurements. We performed an experiment at 7 T using an eight-channel coil and a uniform phantom. RESULTS: We showed that GMT could estimate relative permittivity and conductivity from noisy magnetic resonance measurements with an average error as low as 0.3% and 0.2%, respectively, over the entire volume of the numerical phantom. Voxel resolution did not affect GMT performance and is currently limited only by the memory of the graphics processing unit. In the experiment, GMT could estimate electrical properties within 5% of the values measured with a dielectric probe. CONCLUSION: This work demonstrated the feasibility of GMT with match regularization, suggesting that it could be effective for accurate in vivo electrical property estimation. GMT does not rely on any symmetry assumption for the electromagnetic field, and can be generalized to estimate also the spin magnetization, at the expense of increased computational complexity. SIGNIFICANCE: GMT could provide insight into the distribution of electromagnetic fields inside the body, which represents one of the key ongoing challenges for various diagnostic and therapeutic applications.

Circulant preconditioning in the volume integral equation method for silicon photonics.

Recently, the volume integral equation (VIE) approach has been proposed as an efficient simulation tool for silicon photonics applications [J. Lightw. Technol.36, 3765 (2018)JLTEDG0733-872410.1109/JLT.2018.2842054]. However, for the high-frequency and strong-contrast problems arising in photonics, the convergence of iterative solvers for the solution of the linear system can be extremely slow. The uniform discretization of the volume integral operator leads to a three-level Toeplitz matrix, which is well suited to preconditioning via its circulant approximation. In this paper, we describe an effective circulant preconditioning strategy based on the multilevel circulant preconditioner of Chan and Olkin [Numer. Algorithms6, 89 (1994)NUALEG1017-139810.1007/BF02149764]. We show that this approach proves ideal in the canonical photonics problem of propagation within a uniform waveguide, in which the flow is unidirectional. For more complex photonics structures, such as Bragg gratings, directional couplers, and disk resonators, we generalize our preconditioning strategy via geometrical partitioning (leading to a block-diagonal circulant preconditioner) and homogenization (for inhomogeneous structures). Finally, we introduce a novel memory reduction technique enabling the preconditioner's memory footprint to remain manageable, even for extremely long structures. The range of numerical results we present demonstrates that the preconditioned VIE is fast and has great utility for the numerical exploration of prototype photonics devices.

Computation of ultimate SAR amplification factors for radiofrequency hyperthermia in non-uniform body models: impact of frequency and tumour location.

PURPOSE: We introduce a method for calculation of the ultimate specific absorption rate (SAR) amplification factors (uSAF) in non-uniform body models. The uSAF is the greatest possible SAF achievable by any hyperthermia (HT) phased array for a given frequency, body model and target heating volume. METHODS: First, we generate a basis-set of solutions to Maxwell's equations inside the body model. We place a large number of electric and magnetic dipoles around the body model and excite them with random amplitudes and phases. We then compute the electric fields created in the body model by these excitations using an ultra-fast volume integral solver called MARIE. We express the field pattern that maximises the SAF in the target tumour as a linear combination of these basis fields and optimise the combination weights so as to maximise SAF (concave problem). We compute the uSAFs in the Duke body models at 10 frequencies in the 20-900 MHz range and for twelve 3 cm-diameter tumours located at various depths in the head and neck. RESULTS: For both shallow and deep tumours, the frequency yielding the greatest uSAF was ∼900 MHz. Since this is the greatest frequency that we simulated, we hypothesise that the globally optimal frequency is actually greater. CONCLUSIONS: The uSAFs computed in this work are very large (40-100 for shallow tumours and 4-17 for deep tumours), indicating that there is a large room for improvement of the current state-of-the-art head and neck HT devices.

The ultimate signal-to-noise ratio in realistic body models.

PURPOSE: We compute the ultimate signal-to-noise ratio (uSNR) and G-factor (uGF) in a realistic head model from 0.5 to 21 Tesla. METHODS: We excite the head model and a uniform sphere with a large number of electric and magnetic dipoles placed at 3 cm from the object. The resulting electromagnetic fields are computed using an ultrafast volume integral solver, which are used as basis functions for the uSNR and uGF computations. RESULTS: Our generalized uSNR calculation shows good convergence in the sphere and the head and is in close agreement with the dyadic Green's function approach in the uniform sphere. In both models, the uSNR versus B0 trend was linear at shallow depths and supralinear at deeper locations. At equivalent positions, the rate of increase of the uSNR with B0 was greater in the sphere than in the head model. The uGFs were lower in the realistic head than in the sphere for acceleration in the anterior-posterior direction, but similar for the left-right direction. CONCLUSION: The uSNR and uGFs are computable in nonuniform body models and provide fundamental performance limits for human imaging with close-fitting MRI array coils. Magn Reson Med 78:1969-1980, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Fundamental limits to optical response in absorptive systems.

At visible and infrared frequencies, metals show tantalizing promise for strong subwavelength resonances, but material loss typically dampens the response. We derive fundamental limits to the optical response of absorptive systems, bounding the largest enhancements possible given intrinsic material losses. Through basic conservation-of-energy principles, we derive geometry-independent limits to per-volume absorption and scattering rates, and to local-density-of-states enhancements that represent the power radiated or expended by a dipole near a material body. We provide examples of structures that approach our absorption and scattering limits at any frequency; by contrast, we find that common "antenna" structures fall far short of our radiative LDOS bounds, suggesting the possibility for significant further improvement. Underlying the limits is a simple metric, |χ|2/Im χ for a material with susceptibility χ, that enables broad technological evaluation of lossy materials across optical frequencies.

Fast Electromagnetic Analysis of MRI Transmit RF Coils Based on Accelerated Integral Equation Methods.

A fast frequency domain full-wave electromagnetic simulation method is introduced for the analysis of MRI coils loaded with the realistic human body models. The approach is based on integral equation methods decomposed into two domains: 1) the RF coil array and shield, and 2) the human body region where the load is placed. The analysis of multiple coil designs is accelerated by introducing the precomputed magnetic resonance Green functions (MRGFs), which describe how the particular body model used responds to the incident fields from external sources. These MRGFs, which are precomputed once for a given body model, can be combined with any integral equation solver and reused for the analysis of many coil designs. This approach provides a fast, yet comprehensive, analysis of coil designs, including the port S-parameters and the electromagnetic field distribution within the inhomogeneous body. The method solves the full-wave electromagnetic problem for a head array in few minutes, achieving a speed up of over 150 folds with root mean square errors in the electromagnetic field maps smaller than 0.4% when compared to the unaccelerated integral equation-based solver. This enables the characterization of a large number of RF coil designs in a reasonable time, which is a first step toward an automatic optimization of multiple parameters in the design of transmit arrays, as illustrated in this paper, but also receive arrays.