This study describes a whole-body, non-contact electromagnetic stimulation device based on the concept of a conventional MRI Radio Frequency (RF) resonating coil, but at a much lower resonant frequency (100-150 kHz), with a field modulation option (0.5-100 Hz) and with an input power of up to 3 kW. Its unique features include a high electric field level within the biological tissue due to the resonance effect and a low power dissipation level, or a low Specific Absorption Rate (SAR), in the body itself. Because of its large resonator volume together with non-contact coupling, the subject may be located anywhere within the coil over a longer period at moderate and safe electric field levels. The electric field effect does not depend on body position within the resonator. However, field penetration is deep anywhere within the body, including the extremities where muscles, bones, and peripheral tissues are mostly affected. A potential clinical application of this device is treatment of chronic pain. Substantial attention is paid to device safety; this includes both AC power safety and exposure of human subjects to electromagnetic fields. In the former case, we employ inductive coupling which eliminates a direct current path from AC power to the coil. Our design enhances overall device safety at any power level, even when operated under higher-power conditions. Human exposure to electromagnetic fields within the coil is evaluated by performing modeling with two independent numerical methods and with an anatomically realistic multi-tissue human phantom. We show that SAR levels within the body correspond to International Electrotechnical Commission (IEC) safety standards when the input power level of the amplifier driver does not exceed 3 kW. We also show that electric field levels generally comply with International Commission on Non-Ionizing Radiation Protection safety standards if the input power level does not exceed 1.5 kW.

We present a general numerical approach for solving the forward problem in high-resolution. This approach can be employed in the analysis of noninvasive electroencephalography (EEG) and magnetoencephalography (MEG) as well as invasive electrocorticography (ECoG), stereoencephalography (sEEG), and local field potential (LFP) recordings. The underlying algorithm is our recently developed boundary element fast multipole method (BEM-FMM) that simulates anatomically realistic head models with unprecedented numerical accuracy and speed. This is achieved by utilizing the adjoint double layer formulation and zeroth-order basis functions in conjunction with the FMM acceleration. We present the mathematical formalism in detail and validate the method by applying it to the canonical multilayer sphere problem. The numerical error of BEM-FMM is 2-10 times lower while the computational speed is 1.5-20 times faster than those of the standard first-order FEM. We present four practical case studies: (i) evaluation of the effect of a detailed head model on the accuracy of EEG/MEG forward solution; (ii) demonstration of the ability to accurately calculate the electric potential and the magnetic field in the immediate vicinity of the sources and conductivity boundaries; (iii) computation of the field of a spatially extended cortical equivalent dipole layer; and (iv) taking into account the effect a fontanel for infant EEG source modeling and comparison of the results with a commercially available FEM. In all cases, BEM-FMM provided versatile, fast, and accurate high-resolution modeling of the electromagnetic field and has the potential of becoming a standard tool for modeling both extracranial and intracranial electrophysiological signals.