Abstract
OBJECTIVE: We present a PNS oracle, which solves these computation time and linearity problems and is, therefore, well-suited for fast optimization of voltage distributions in contact electrode arrays and current drive patterns in non-contact magnetic coil arrays. APPROACH: The PNS oracle metric for a nerve fiber is computed from an electric field map using only linear operations (projection, differentiation, convolution, scaling). Due to its linearity, this PNS metric can be precomputed for a set of coil or electrode segments, allowing rapid PNS prediction and comparison of any possible coil or electrode stimulation configuration constructed from this set. The PNS oracle is closely related to the classical activating function and modified driving functions but is adjusted to better correlate with full neurodynamic modeling of myelinated mammalian nerves. MAIN RESULTS: We validated the PNS oracle in three MRI gradient coils and two body models and found good correlation between the PNS oracle and the full neurodynamic modeling approach (R (2) > 0.995). Finally, we demonstrated its potential utility by optimizing the driving currents and voltages of arrays of 108 magnetic coils or 108 contact electrodes to selectively stimulate target nerves in the lower leg. SIGNIFICANCE: Peripheral nerve stimulation (PNS) by electromagnetic fields can be accurately simulated using coupled electromagnetic and neurodynamic modeling. Such simulations are slow and non-linear in the electric field, which makes it difficult to iteratively optimize coil and electrode configurations or drive patterns aiming to avoid PNS or to initiate it for therapeutic purposes.