Uncertainty Quantification in SAR Induced by Ultra-High-Field MRI RF Coil via High-Dimensional Model Representation.

As magnetic field strength in Magnetic Resonance Imaging (MRI) technology increases, maintaining the specific absorption rate (SAR) within safe limits across human head tissues becomes challenging due to the formation of standing waves at a shortened wavelength. Compounding this challenge is the uncertainty in the dielectric properties of head tissues, which notably affects the SAR induced by the radiofrequency (RF) coils in an ultra-high-field (UHF) MRI system. To this end, this study introduces a computational framework to quantify the impacts of uncertainties in head tissues' dielectric properties on the induced SAR. The framework employs a surrogate model-assisted Monte Carlo (MC) technique, efficiently generating surrogate models of MRI observables (electric fields and SAR) and utilizing them to compute SAR statistics. Particularly, the framework leverages a high-dimensional model representation technique, which constructs the surrogate models of the MRI observables via univariate and bivariate component functions, approximated through generalized polynomial chaos expansions. The numerical results demonstrate the efficiency of the proposed technique, requiring significantly fewer deterministic simulations compared with traditional MC methods and other surrogate model-assisted MC techniques utilizing machine learning algorithms, all while maintaining high accuracy in SAR statistics. Specifically, the proposed framework constructs surrogate models of a local SAR with an average relative error of 0.28% using 289 simulations, outperforming the machine learning-based surrogate modeling techniques considered in this study. Furthermore, the SAR statistics obtained by the proposed framework reveal fluctuations of up to 30% in SAR values within specific head regions. These findings highlight the critical importance of considering dielectric property uncertainties to ensure MRI safety, particularly in 7 T MRI systems.

A Compact Dual-Polarized Vivaldi Antenna with High Gain for Tree Radar Applications.

A dual-polarized compact Vivaldi antenna with high gain performance is proposed for tree radar applications. The proposed design introduces an array configuration consisting of two pairs of two Vivaldi elements to optimize the operating bandwidth and gain while providing dual-polarization capability. To enhance the gain of the proposed antenna over a certain frequency range of interest, directors and edge slots are incorporated into each Vivaldi element. To further enhance the overall antenna gain, a metal back reflector is used. The measurement results of the fabricated antenna show that the proposed antenna achieves a high gain of 5.5 to 14.8 dBi over broadband from 0.5 GHz to 3 GHz. Moreover, it achieves cross-polarization discrimination larger than 20 dB, ensuring high polarization purity. The fabricated antenna is used to detect and image the defects inside tree trunks. The results show that the proposed antenna yields a better-migrated image with a clear defect region compared to that obtained by a commercial Horn antenna.

Fast computational E-field dosimetry for transcranial magnetic stimulation using adaptive cross approximation and auxiliary dipole method (ACA-ADM).

Transcranial Magnetic Stimulation (TMS) is a non-invasive brain stimulation technique that uses a coil to induce an electric field (E-field) in the brain and modulate its activity. Many applications of TMS call for the repeated execution of E-field solvers to determine the E-field induced in the brain for different coil placements. However, the usage of solvers for these applications remains impractical because each coil placement requires the solution of a large linear system of equations. We develop a fast E-field solver that enables the rapid evaluation of the E-field distribution for a brain region of interest (ROI) for a large number of coil placements, which is achieved in two stages. First, during the pre-processing stage, the mapping between coil placement and brain ROI E-field distribution is approximated from E-field results for a few coil placements. Specifically, we discretize the mapping into a matrix with each column having the ROI E-field samples for a fixed coil placement. This matrix is approximated from a few of its rows and columns using adaptive cross approximation (ACA). The accuracy, efficiency, and applicability of the new ACA approach are determined by comparing its E-field predictions with analytical and standard solvers in spherical and MRI-derived head models. During the second stage, the E-field distribution in the brain ROI from a specific coil placement is determined by the obtained rows and columns in milliseconds. For many applications, only the E-field distribution for a comparatively small ROI is required. For example, the solver can complete the pre-processing stage in approximately 4 hours and determine the ROI E-field in approximately 40 ms for a 100 mm diameter ROI with less than 2% error enabling its use for neuro-navigation and other applications. Highlight: We developed a fast solver for TMS computational E-field dosimetry, which can determine the ROI E-field in approximately 40 ms for a 100 mm diameter ROI with less than 2% error.

DeeptDCS: Deep Learning-Based Estimation of Currents Induced During Transcranial Direct Current Stimulation.

OBJECTIVE: Transcranial direct current stimulation (tDCS) is a non-invasive brain stimulation technique used to generate conduction currents in the head and disrupt brain functions. To rapidly evaluate the tDCS-induced current density in near real-time, this paper proposes a deep learning-based emulator, named DeeptDCS. METHODS: The emulator leverages Attention U-net taking the volume conductor models (VCMs) of head tissues as inputs and outputting the three-dimensional current density distribution across the entire head. The electrode configurations are also incorporated into VCMs without increasing the number of input channels; this enables the straightforward incorporation of the non-parametric features of electrodes (e.g., thickness, shape, size, and position) in the training and testing of the proposed emulator. RESULTS: Attention U-net outperforms standard U-net and its other three variants (Residual U-net, Attention Residual U-net, and Multi-scale Residual U-net) in terms of accuracy. The generalization ability of DeeptDCS to non-trained electrode configurations can be greatly enhanced through fine-tuning the model. The computational time required by one emulation via DeeptDCS is a fraction of a second. CONCLUSION: DeeptDCS is at least two orders of magnitudes faster than a physics-based open-source simulator, while providing satisfactorily accurate results. SIGNIFICANCE: The high computational efficiency permits the use of DeeptDCS in applications requiring its repetitive execution, such as uncertainty quantification and optimization studies of tDCS.

A Butterfly-Accelerated Volume Integral Equation Solver for Broad Permittivity and Large-Scale Electromagnetic Analysis.

A butterfly-accelerated volume integral equation (VIE) solver is proposed for fast and accurate electromagnetic (EM) analysis of scattering from heterogeneous objects. The proposed solver leverages the hierarchical off-diagonal butterfly (HOD-BF) scheme to construct the system matrix and obtain its approximate inverse, used as a preconditioner. Complexity analysis and numerical experiments validate the O(N log 2 N) construction cost of the HOD-BF-compressed system matrix and O(N 1.5 log N) inversion cost for the preconditioner, where N is the number of unknowns in the high-frequency EM scattering problem. For many practical scenarios, the proposed VIE solver requires less memory and computational time to construct the system matrix and obtain its approximate inverse compared to a H matrix-accelerated VIE solver. The accuracy and efficiency of the proposed solver have been demonstrated via its application to the EM analysis of large-scale canonical and real-world structures comprising of broad permittivity values and involving millions of unknowns.

Compact Dual-Polarized Vivaldi Antenna with High Gain and High Polarization Purity for GPR Applications.

A compact ultra-wideband dual-polarized Vivaldi antenna is proposed for full polarimetric ground-penetrating radar (GPR) applications. A shared-aperture configuration comprising four Vivaldi elements for orthogonal polarizations is designed to reduce the low-end operating frequency and improve the port isolation with a compact antenna size. The directivity of the antenna is enhanced by the oblique position of the radiators and the implementation of a square loop reflector. Experimental results demonstrate that the antenna has very good impedance matching, port isolation, and dual-polarized radiation performance, with low dispersion characteristics across band of interest from 0.4 GHz to 3.0 GHz. GPR measurements with the designed antenna show that the antenna maintains good detection capability even for objects buried in a highly conductive soil.

An FMM-FFT Accelerated SIE Simulator for Analyzing EM Wave Propagation in Mine Environments Loaded With Conductors.

A fast and memory efficient three-dimensional full-wave simulator for analyzing electromagnetic (EM) wave propagation in electrically large and realistic mine tunnels/galleries loaded with conductors is proposed. The simulator relies on Muller and combined field surface integral equations (SIEs) to account for scattering from mine walls and conductors, respectively. During the iterative solution of the system of SIEs, the simulator uses a fast multipole method-fast Fourier transform (FMM-FFT) scheme to reduce CPU and memory requirements. The memory requirement is further reduced by compressing large data structures via singular value and Tucker decompositions. The efficiency, accuracy, and real-world applicability of the simulator are demonstrated through characterization of EM wave propagation in electrically large mine tunnels/galleries loaded with conducting cables and mine carts.

The ICVSIE: A General Purpose Integral Equation Method for Bio-Electromagnetic Analysis.

OBJECTIVE: An internally combined volume surface integral equation (ICVSIE) for analyzing electromagnetic (EM) interactions with biological tissue and wide ranging diagnostic, therapeutic, and research applications, is proposed. METHOD: The ICVSIE is a system of integral equations in terms of volume and surface equivalent currents in biological tissue subject to fields produced by externally or internally positioned devices. The system is created by using equivalence principles and solved numerically; the resulting current values are used to evaluate scattered and total electric fields, specific absorption rates, and related quantities. RESULTS: The validity, applicability, and efficiency of the ICVSIE are demonstrated by EM analysis of transcranial magnetic stimulation, magnetic resonance imaging, and neuromuscular electrical stimulation. CONCLUSION: Unlike previous integral equations, the ICVSIE is stable regardless of the electric permittivities of the tissue or frequency of operation, providing an application-agnostic computational framework for EM-biomedical analysis. SIGNIFICANCE: Use of the general purpose and robust ICVSIE permits streamlining the development, deployment, and safety analysis of EM-biomedical technologies.

Uncertainty quantification in transcranial magnetic stimulation via high-dimensional model representation.

A computational framework for uncertainty quantification in transcranial magnetic stimulation (TMS) is presented. The framework leverages high-dimensional model representations (HDMRs), which approximate observables (i.e., quantities of interest such as electric (E) fields induced inside targeted cortical regions) via series of iteratively constructed component functions involving only the most significant random variables (i.e., parameters that characterize the uncertainty in a TMS setup such as the position and orientation of TMS coils, as well as the size, shape, and conductivity of the head tissue). The component functions of HDMR expansions are approximated via a multielement probabilistic collocation (ME-PC) method. While approximating each component function, a quasi-static finite-difference simulator is used to compute observables at integration/collocation points dictated by the ME-PC method. The proposed framework requires far fewer simulations than traditional Monte Carlo methods for providing highly accurate statistical information (e.g., the mean and standard deviation) about the observables. The efficiency and accuracy of the proposed framework are demonstrated via its application to the statistical characterization of E-fields generated by TMS inside cortical regions of an MRI-derived realistic head model. Numerical results show that while uncertainties in tissue conductivities have negligible effects on TMS operation, variations in coil position/orientation and brain size significantly affect the induced E-fields. Our numerical results have several implications for the use of TMS during depression therapy: 1) uncertainty in the coil position and orientation may reduce the response rates of patients; 2) practitioners should favor targets on the crest of a gyrus to obtain maximal stimulation; and 3) an increasing scalp-to-cortex distance reduces the magnitude of E-fields on the surface and inside the cortex.