<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Comput. Neurosci.</journal-id>
<journal-title>Frontiers in Computational Neuroscience</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Comput. Neurosci.</abbrev-journal-title>
<issn pub-type="epub">1662-5188</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fncom.2023.1105505</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Computational Neuroscience</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Improving focality and consistency in micromagnetic stimulation</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Ye</surname> <given-names>Hui</given-names></name>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1793778/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Hall</surname> <given-names>Vincent</given-names></name>
</contrib>
<contrib contrib-type="author">
<name><surname>Hendee</surname> <given-names>Jenna</given-names></name>
</contrib>
</contrib-group>
<aff><institution>Department of Biology, Loyola University Chicago</institution>, <addr-line>Chicago, IL</addr-line>, <country>United States</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Ursula van Rienen, University of Rostock, Germany</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Thomas R. Kn&#x00F6;sche, Max Planck Institute for Human Cognitive and Brain Sciences, Germany; Hisayuki Osanai, University of Texas Southwestern Medical Center, United States; Mattia Stefano, Campus Bio-Medico University, Italy</p></fn>
<corresp id="c001">&#x002A;Correspondence: Hui Ye, <email>hye1@luc.edu</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>02</day>
<month>02</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>17</volume>
<elocation-id>1105505</elocation-id>
<history>
<date date-type="received">
<day>22</day>
<month>11</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>17</day>
<month>01</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2023 Ye, Hall and Hendee.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Ye, Hall and Hendee</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>The novel micromagnetic stimulation (&#x03BC;MS) technology aims to provide high resolution on neuronal targets. However, consistency of neural activation could be compromised by a lack of surgical accuracy, biological variation, and human errors in operation. We have recently modeled the activation of an unmyelinated axon by a circular micro-coil. Although the coil could activate the axon, its performance sometimes lacked focality and consistency. The site of axonal activation could shift by several experimental factors, including the reversal of the coil current, displacement of the coil, and changes in the intensity of the stimulation. Current clinical practice with transcranial magnetic stimulation (TMS) has suggested that figure-eight coils could provide better performance in magnetic stimulation than circular coils. Here, we estimate the performance of &#x03BC;MS by a figure-eight micro-coil, by exploring the impact of the same experimental factors on its focality and consistency in axonal activation. We derived the analytical expression of the electric field and activating function generated by the figure-eight micro-coil, and estimated the location of axonal activation. Using NEURON modeling of an unmyelinated axon, we found two different types (A and B) of axon activation by the figure-eight micro-coil, mediated by coil currents of reversed direction. Type A activation is triggered by membrane hyperpolarization followed by depolarization; Type B activation is triggered by direct membrane depolarization. Consequently, the two types of stimulation are governed by distinct ion channel mechanisms. In comparison to the circular micro-coil, the figure-eight micro-coil requires significantly less current for axonal activation. Under figure-eight micro-coil stimulation, the site of axonal activation does not change with the reversal of the coil current, displacement of the coil, or changes in the intensity of the stimulation. Ultimately, the figure-eight micro-coil provides a more efficient and consistent site of activation than the circular micro-coil in &#x03BC;MS.</p>
</abstract>
<kwd-group>
<kwd>micromagnetic stimulation (&#x03BC;MS)</kwd>
<kwd>figure-eight micro-coil</kwd>
<kwd>activating function</kwd>
<kwd>NEURON modeling</kwd>
<kwd>ion channels</kwd>
<kwd>consistency</kwd>
</kwd-group>
<contract-sponsor id="cn001">National Institutes of Health<named-content content-type="fundref-id">10.13039/100000002</named-content></contract-sponsor>
<counts>
<fig-count count="12"/>
<table-count count="1"/>
<equation-count count="27"/>
<ref-count count="66"/>
<page-count count="18"/>
<word-count count="13077"/>
</counts>
</article-meta>
</front>
<body>
<sec id="S1" sec-type="intro">
<title>1. Introduction</title>
<p>Transcranial magnetic stimulation (TMS) is a noninvasive brain stimulation technique used in the clinic to treat several neurological disorders and psychiatric diseases. One of TMS&#x2019;s major limitations is its low spatial resolution, which often results in a mismatch between the target area in the brain and the stimulation site on the scalp. Early TMS practices using a single circular coil were soon found to experience difficulties in locally stimulating the targeted areas in the brain (<xref ref-type="bibr" rid="B3">Barker et al., 1985</xref>). It took only a few years before the figure-eight coil was developed for more localized brain stimulation (<xref ref-type="bibr" rid="B55">Ueno et al., 1988</xref>). Using idealized coils consisting of one or two wire loops, it was soon demonstrated that figure-eight coils allow for more focused stimulation than simple round coils (<xref ref-type="bibr" rid="B44">Ravazzani et al., 1996</xref>). In recent years, different variants of the figure-eight coils have been developed for basic research and clinical purposes (<xref ref-type="bibr" rid="B54">Ueno and Sekino, 2021</xref>), such as the quadruple butterfly coil (<xref ref-type="bibr" rid="B40">Rastogi et al., 2017</xref>), eccentric figure-eight coil (<xref ref-type="bibr" rid="B49">Sekino et al., 2015</xref>), and figure-eight coil with an iron core (<xref ref-type="bibr" rid="B58">Yamamoto et al., 2016</xref>), etc. Applying these novel designs to neural stimulation has led to an in-depth understanding of the functional organization of the human brain, dynamic neuronal connectivity, and neuronal plasticity in the cortex.</p>
<p>An effective way to increase the spatial resolution of magnetic stimulation is by reducing the dimensionality of the stimulating coil (<xref ref-type="bibr" rid="B51">Talebinejad and Musallam, 2010</xref>; <xref ref-type="bibr" rid="B53">Tischler et al., 2011</xref>). This led to the recent development of micromagnetic stimulation (&#x03BC;MS) technology, which significantly increases the focality of magnetic stimulation. This technology uses millimeter or submillimeter coils to target individual neurons (<xref ref-type="bibr" rid="B5">Bonmassar et al., 2012</xref>; <xref ref-type="bibr" rid="B61">Ye and Barrett, 2021</xref>) or specific neuronal populations (<xref ref-type="bibr" rid="B27">Lee and Fried, 2017</xref>). Due to their microscopic size, these coils are meant to be implanted close to the targeted neurons or axons under the cover of biocompatible materials (<xref ref-type="bibr" rid="B34">Park et al., 2013</xref>). This prevents direct contact between the coil and the neural tissue, mitigating numerous adverse effects that arise at the brain-electrode interface when metal electrodes are implanted for neural stimulation (<xref ref-type="bibr" rid="B37">Polikov et al., 2005</xref>; <xref ref-type="bibr" rid="B7">Cogan, 2008</xref>; <xref ref-type="bibr" rid="B24">Koivuniemi et al., 2011</xref>). Several works have used &#x03BC;MS in axonal activation, including the axons of the hippocampus (<xref ref-type="bibr" rid="B45">Saha et al., 2022a</xref>), apical dendrites of layer V pyramidal neurons in the cortex (<xref ref-type="bibr" rid="B27">Lee and Fried, 2017</xref>), and axons in the dorsal cochlear nucleus (<xref ref-type="bibr" rid="B11">Golestanirad et al., 2018</xref>).</p>
<p>The ultimate goal of the novel &#x03BC;MS technology is to achieve highly specific neural stimulation with improved consistency, defined by a micro-coil&#x2019;s ability to repetitively generate the same neural response. However, there is a considerable chance that &#x03BC;MS could generate inconsistent results, introduced by a lack of surgical accuracy, multiple biological variations, and human errors in operation.</p>
<p>The permanent and irreversible implantation surgery could cause inconsistency in &#x03BC;MS. Because of its small dimension, it is difficult to implant the micro-coil at the same anatomical location across different subjects. Recent technological developments in surgical robotics can mitigate the effects of human hand tremors within the micrometer scale (<xref ref-type="bibr" rid="B1">Ahmed et al., 2018</xref>). However, placement of implants is still at the millimeter (mm) level of precision (<xref ref-type="bibr" rid="B9">Fernandes de Oliveira Santos et al., 2021</xref>; <xref ref-type="bibr" rid="B66">Zhang et al., 2021</xref>). Therefore, even with assistance from surgical robots and software, it is difficult to achieve implantation that is precise to the micrometer (&#x03BC;m) scale. It is highly important to ensure the stimulation apparatus (i.e., a micro-coil) is as focal and consistent as possible during &#x03BC;MS.</p>
<p>Variance in the biological system could cause inconsistency in &#x03BC;MS. First, neural activation is determined by the direction, intensity, and duration of the electric field and its spatial derivative along the neural tissue (<xref ref-type="bibr" rid="B36">Pashut et al., 2014</xref>). Consequently, the orientation of the targeted neurons to the micro-coil is crucial for neural activation (<xref ref-type="bibr" rid="B11">Golestanirad et al., 2018</xref>). Yet, ensuring the implanted coil has the ideal orientation to the targeted tissue is challenging. Second, most clinical purposes require a high degree of accuracy and precision when small neuronal volumes are targeted for treatment. However, numerous factors could cause migration of the implanted micro-device (<xref ref-type="bibr" rid="B56">van den Munckhof et al., 2010</xref>; <xref ref-type="bibr" rid="B35">Park et al., 2011</xref>), such as pulsing of blood in the vessels, breathing, or simple locomotion. Device migration would alter the neuronal population being activated.</p>
<p>Human errors in operating the stimulation protocols could cause inconsistency in &#x03BC;MS. Accidentally choosing the wrong intensity and polarity for stimulation could lead to inconsistent or unwanted outcomes. In TMS practice, the stimulation outcome is routinely monitored with fMRI or EEG recording (<xref ref-type="bibr" rid="B4">Bohning et al., 1999</xref>; <xref ref-type="bibr" rid="B13">Habibollahi Saatlou et al., 2018</xref>). This is not an option in &#x03BC;MS, for which the stimulation site is microscopic and only a small population of neurons are affected by the micro-coil. After the micro-coil is chronically implanted, it is difficult to detect the consequence of these human errors by acquiring feedback at the single cell level. Understanding these variabilities and their respective impacts on stimulation outcomes could significantly improve the development of &#x03BC;MS technology and maximize the device&#x2019;s success.</p>
<p>To further increase the focality of &#x03BC;MS and the directionality of the neural response, researchers in the &#x03BC;MS field have started to explore the possibility of combining multiple coils to increase stimulation effectiveness. <xref ref-type="bibr" rid="B20">Jeong et al. (2021)</xref> designed a figure-eight spiral micro-coil to generate high magnetic flux for ultra-focal stimulation. Saito constructed a figure-eight micro-coil using commercially available chip inductors, which effectively suppressed synchronized bursting activity in a cultured neural network (<xref ref-type="bibr" rid="B47">Saito, 2021</xref>). <xref ref-type="bibr" rid="B8">Colella et al. (2021)</xref> simulated axonal stimulation of peripheral nerves with the figure-eight micro-coil, and found that the threshold for axonal activation was lower for the figure-eight coil than for a circular micro-coil.</p>
<p>Computational works have provided invaluable insights to estimate the focality of micro-coil stimulation. Researchers have estimated the induced electric field distribution around the single coil (<xref ref-type="bibr" rid="B5">Bonmassar et al., 2012</xref>; <xref ref-type="bibr" rid="B33">Osanai et al., 2018</xref>; <xref ref-type="bibr" rid="B63">Ye et al., 2020</xref>) and multi-coil systems (<xref ref-type="bibr" rid="B31">Minusa et al., 2019</xref>). However, when the strength of the field induced by the micro-coil is small, the coil can still produce a spatial gradient that can activate the neuron (<xref ref-type="bibr" rid="B25">Lee et al., 2022</xref>). This field gradient, or activating function, is directly related to the site of neural activation (<xref ref-type="bibr" rid="B41">Rattay, 1986</xref>; <xref ref-type="bibr" rid="B27">Lee and Fried, 2017</xref>). Therefore, several research endeavors have computed the activating function for the single coil system to identify the site of neural activation (<xref ref-type="bibr" rid="B28">Lee et al., 2016</xref>; <xref ref-type="bibr" rid="B60">Ye, 2022</xref>). We recently simulated axonal activation by a circular micro-coil using activating function analysis and NEURON simulation (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>). This work suggests that &#x03BC;MS stimulation with a circular micro-coil could be inconsistent, and that the location of axonal activation could shift by several experimental factors, including the reversal of the coil current, increases in the coil-axon distance, and changes in the intensity of stimulation. When a micro-coil with a dimension of 1 mm was tested, the activation site on the axon could shift by several hundred micrometers (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>). So far, activating function analysis has not been performed to estimate the focality and consistency of stimulation in a figure-eight coil system for the activation of a straight axon. Furthermore, although it was believed that the figure-eight micro-coil could provide a better &#x201C;coupling between nerve fibers and the induced electric field&#x201D; (<xref ref-type="bibr" rid="B8">Colella et al., 2021</xref>), there is a lack of understanding of the neuronal mechanisms underlying axonal activation by such a coil.</p>
<p>This paper will examine the consistency of axonal activation by a figure-eight micro-coil, and investigate the underlying ion channel mechanisms using combined activating function analysis and NEURON modeling. We will compare these results from our recent analysis on a single, circular micro-coil (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>). This work provides valuable insights to the further development of &#x03BC;MS technology for focal and consistent neural stimulation.</p>
</sec>
<sec id="S2" sec-type="materials|methods">
<title>2. Materials and methods</title>
<sec id="S2.SS1">
<title>2.1. Axon model under magnetic stimulation in a figure-eight coil system</title>
<p>Previously, we (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>) have modeled a submillimeter, circular micro-coil with radius <italic>Rc</italic> (<xref ref-type="fig" rid="F1">Figure 1A</xref>). To model a figure-eight micro-coil, we positioned two identical micro-coils next to each other, but with reversed winding direction (<xref ref-type="fig" rid="F1">Figure 1B</xref>). The center of the left coil was at point O<sub>1</sub> (<italic>&#x2212;Rc</italic>, 0, 0), and the center of the right coil was at point O<sub>2</sub> (<italic>Rc</italic>, 0, 0). Electric currents in the two coils were synchronized, but with reversed directions. This figure-eight micro-coil was positioned so that its induced electric field was in parallel with the axon, to generate effective stimulation (<xref ref-type="bibr" rid="B18">Jefferys, 1981</xref>; <xref ref-type="bibr" rid="B10">Gluckman et al., 1996</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p>Coordination system that describes the orientation of a circular micro-coil versus a figure-eight micro-coil for the activation of an unmyelinated axon. The coil current (<italic>I</italic>) generated a time-varying magnetic field (<inline-formula><mml:math id="INEQ1"><mml:mover accent='true'><mml:mi>B</mml:mi><mml:mo>&#x2192;</mml:mo></mml:mover></mml:math></inline-formula>), which, in turn, induced an electric field (<inline-formula><mml:math id="INEQ2"><mml:mover accent='true'><mml:mi>E</mml:mi><mml:mo>&#x2192;</mml:mo></mml:mover></mml:math></inline-formula>) to stimulate the axon. The multi-compartment model of the unmyelinated axon was 20,000 &#x03BC;m in length and was divided into 200 segments. Each segment was a cylinder of length 100 &#x03BC;m and diameter 15 &#x03BC;m. Hodgkin-Huxley type ion channels were inserted into each segment. <bold>(A)</bold> The circular micro-coil had a radius <italic>Rc</italic>. The center of the coil was at O (0,0,0) in the Cartesian coordinate system. The axon was in the x-y plan and was parallel to the x-axis. Increase in the coil current (<italic>I, dI/dt</italic> &#x003E; 0) in the clockwise direction generated a counterclockwise electric field and a field gradient along the axon. <bold>(B)</bold> The figure-eight micro-coil was made with two identical, circular micro-coils (L and R) located next to each other. The center of the left coil was O<sub>1</sub> (<italic>&#x2013;Rc</italic>, 0, 0) and the center for the right coil was O<sub>2</sub> (<italic>Rc</italic>, 0, 0). Coil current (I) was clockwise in the L-coil, and counterclockwise in the R-coil. The current was increasing with time (<italic>dI/dt</italic> &#x003E; 0).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g001.tif"/>
</fig>
<p>The unmyelinated axon was represented with a multi-compartment model in NEURON (v.7.8) simulation environment package (<xref ref-type="bibr" rid="B14">Hines and Carnevale, 1997</xref>). The model simulated the axon as a cylinder 20,000 &#x03BC;m in length and 15 &#x03BC;m in diameter. The axon was divided evenly into 200 node segments [Table 1 in <xref ref-type="bibr" rid="B60">Ye (2022)</xref>]. Each segment contained the Hodgkin-Huxley (H-H) types of ion channels, including the fast sodium, slow potassium, and leakage channels in the membrane nodes (<xref ref-type="bibr" rid="B15">Hodgkin and Huxley, 1952</xref>). The ionic current (I) at the n-th segment of the neuron was described as,</p>
<disp-formula id="S2.E1">
<label>(1)</label>
<mml:math id="M1">
<mml:mrow>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
</mml:mpadded>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>g</mml:mi>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mi>m</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>g</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:msup>
<mml:mi>n</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>g</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where V<sub><italic>Na</italic></sub>, V<sub><italic>K</italic></sub>, and V<sub><italic>L</italic></sub> were the equilibrium membrane potentials for sodium, potassium, and leakage channels, respectively. g<sub><italic>Na</italic></sub>, g<sub><italic>k</italic></sub>, and g<sub><italic>L</italic></sub> were the maximal conductance of sodium, potassium, and leakage channels, respectively. <italic>m</italic> and <italic>h</italic> represented the activation and inactivation of the sodium channels, respectively, whereas <italic>n</italic> represented the activation of potassium channels. The state variables <italic>m, h</italic>, and <italic>n</italic> were defined in,</p>
<disp-formula id="S2.E2">
<label>(2)</label>
<mml:math id="M2">
<mml:mrow>
<mml:mpadded width="+3.3pt">
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mpadded>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mn>&#x03B1;</mml:mn>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mn>&#x03B2;</mml:mn>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E3">
<label>(3)</label>
<mml:math id="M3">
<mml:mrow>
<mml:mpadded width="+3.3pt">
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mpadded>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mn>&#x03B1;</mml:mn>
<mml:mi>h</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mn>&#x03B2;</mml:mn>
<mml:mi>h</mml:mi>
</mml:msub>
<mml:mpadded width="+5pt">
<mml:mi>h</mml:mi>
</mml:mpadded>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E4">
<label>(4)</label>
<mml:math id="M4">
<mml:mrow>
<mml:mpadded width="+3.3pt">
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mpadded>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mn>&#x03B1;</mml:mn>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mn>&#x03B2;</mml:mn>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where &#x03B1;<sub><italic>m</italic></sub>, &#x03B2;<sub><italic>m</italic></sub>, &#x03B1;<sub><italic>h</italic></sub>, &#x03B2;<sub><italic>h</italic></sub>, &#x03B1;<sub><italic>n</italic></sub>, and &#x03B2;<sub><italic>n</italic></sub> were rate constants defined previously [Table 2 in <xref ref-type="bibr" rid="B60">Ye (2022)</xref>].</p>
</sec>
<sec id="S2.SS2">
<title>2.2. Electric field induced by a figure-eight micro-coil</title>
<p>Previously, we (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>) have derived the analytical expression of the induced electric field generated by a circular micro-coil (<xref ref-type="fig" rid="F1">Figure 1A</xref>). The electric field generated by the figure-eight micro-coil is, therefore, the summation of the two individual circular micro-coils (Left coil and Right coil, <xref ref-type="fig" rid="F1">Figure 1B</xref>). We applied a single voltage pulse (positive or negative) to the figure-eight coil, since this pulse was sufficient to activate various axons (<xref ref-type="bibr" rid="B5">Bonmassar et al., 2012</xref>; <xref ref-type="bibr" rid="B28">Lee et al., 2016</xref>). Electric fields generated by the left coil (EL) and the right coil (ER) were calculated as follows:</p>
<p>For the onset of the positive stimulation pulse, the electric field generated by the left coil was calculated with Eqs. (13, 14) in <xref ref-type="bibr" rid="B60">Ye (2022)</xref>,</p>
<disp-formula id="S2.E5">
<label>(5-1)</label>
<mml:math id="M5">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mpadded>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E6">
<label>(5-2)</label>
<mml:math id="M6">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mpadded>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Electric field generated by the right coil was,</p>
<disp-formula id="S2.E7">
<label>(5-3)</label>
<mml:math id="M7">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mpadded>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E8">
<label>(5-4)</label>
<mml:math id="M8">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mpadded>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>N</italic> is the number of coil loops, <italic>l</italic> the coil length, <italic>R</italic> the coil resistance, and L the coil inductance. &#x03BC;<sub>0</sub> is the vacuum permeability and <italic>V</italic> is the voltage across the coil. <italic>x</italic> and <italic>y</italic> defined the location considered in the Cartesian coordinate system (<xref ref-type="fig" rid="F1">Figure 1B</xref>), and <italic>t</italic> defined time.</p>
<p>The summarized electric field generated by the figure-eight coil at the onset of a positive stimulation pulse was:</p>
<disp-formula id="S2.Ex9">
<label>(6-1)</label>
<mml:math id="M9">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:mi>x</mml:mi>
</mml:mpadded>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:mi>x</mml:mi>
</mml:mpadded>
</mml:mrow>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E10">
<mml:math id="M10">
<mml:mrow>
<mml:mrow>
<mml:mo maxsize="260%" minsize="260%">[</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo maxsize="260%" minsize="260%">]</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.Ex11">
<label>(6-2)</label>
<mml:math id="M11">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:mi>y</mml:mi>
</mml:mpadded>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mpadded>
</mml:mrow>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E12">
<mml:math id="M12">
<mml:mrow>
<mml:mrow>
<mml:mo maxsize="260%" minsize="260%">[</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo maxsize="260%" minsize="260%">]</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>The intensity of the electric field was calculated with,</p>
<disp-formula id="S2.E13">
<label>(7)</label>
<mml:math id="M13">
<mml:mrow>
<mml:mpadded width="+3.3pt">
<mml:mi>E</mml:mi>
</mml:mpadded>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mpadded>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">+</mml:mo>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>For the offset of the stimulation pulse, the electric field generated by the left coil was calculated by Eqs. (15, 16) from <xref ref-type="bibr" rid="B60">Ye (2022)</xref>.</p>
<disp-formula id="S2.E14">
<label>(8-1)</label>
<mml:math id="M14">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mpadded>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E15">
<label>(8-2)</label>
<mml:math id="M15">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mpadded>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Electric field generated by the right coil was,</p>
<disp-formula id="S2.E16">
<label>(8-3)</label>
<mml:math id="M16">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mpadded>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E17">
<label>(8-4)</label>
<mml:math id="M17">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mpadded>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>The summarized electric field generated by the figure-eight micro-coil at the offset of the stimulation pulse was:</p>
<disp-formula id="S2.Ex18">
<label>(9-1)</label>
<mml:math id="M18">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:mi>x</mml:mi>
</mml:mpadded>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mpadded>
</mml:mrow>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E19">
<mml:math id="M19">
<mml:mrow>
<mml:mrow>
<mml:mo maxsize="260%" minsize="260%">[</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo maxsize="260%" minsize="260%">]</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.Ex20">
<label>(9-2)</label>
<mml:math id="M20">
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:mi>y</mml:mi>
</mml:mpadded>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mpadded>
</mml:mrow>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E21">
<mml:math id="M21">
<mml:mrow>
<mml:mrow>
<mml:mo maxsize="260%" minsize="260%">[</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo maxsize="260%" minsize="260%">]</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Here, <italic>L/R</italic> defined the time constant. Therefore, the single pulse in the figure-eight micro-coil induced a biphasic electric field around the axon. Axon stimulation should occur during the onset phase and/or the offset phase of the stimulus pulse that effectively generates the electric field.</p>
</sec>
<sec id="S2.SS3">
<title>2.3. Measuring the induced electric field generated by the figure-eight micro-coil</title>
<p>To validate the mathematic derivation that the magnetically induced electric field by the figure-eight coil is indeed biphasic, we assembled a figure-eight micro-coil (<xref ref-type="supplementary-material" rid="FS1">Supplementary Figure 1</xref>) using two commercial submillimeter, multilayer surface mount inductors (100 nH, 2 &#x03A9;, MLG1005SR10JTD25, TDK USA. Corporation, Uniondale, NY, USA). The internal structure of the coil was visualized by removing the ceramic core and epoxy coating using 40% liquid hydrofluoric acid and 10 N HCL, based on a published protocol (<xref ref-type="bibr" rid="B63">Ye et al., 2020</xref>). The coil contained 20 wire loops (<italic>N</italic> = 20).</p>
<p>Two copper wires (magnetic wire 32-AWG, GC electronics, IL, L3-616) were inserted into the shaft of a 20 AWG needle and the barrel of a syringe. They were soldered to the metal leads of the two inductors for electric current delivery. The two inductors were positioned side by side but with reversed orientation to ensure opposite winding of the coil. To ensure electrical insulation and water impermeability of the exposed coil terminals during electrophysiological experiments, the coil was then coated with acrylate copolymer enamel (<xref ref-type="bibr" rid="B34">Park et al., 2013</xref>) and allowed 24 h to dry before the experiment. The other end of the copper wires were connected to the pin-connectors. We delivered electric pulses to the figure-eight micro-coil using an arbitrary function generator (AFG1022, Tektronix) and a power amplifier (Pyramid PB 717X 2 channel, Pyramid Car Audio, Brooklyn, NY, 11204, USA). To measure the induced electric field around the figure-eight micro-coil, we filled a petri dish with conductive saline and submerged the figure-eight micro-coil under the saline. We positioned a glass electrode next to the figure-eight micro-coil. The recorded waveform was amplified by a model 1,700 differential AC amplifier (A-M Systems) and stored on a computer with Spike 2 software (v. 7.2 Cambridge Electronic Design Limited).</p>
<p>The assembled figure-eight coils were tested to ensure there was no leakage or current. When the coil was immersed in saline solution, we measured the impedance between one end of the coil to the ground, and found it was greater than 10 M&#x03A9;, suggesting a good insulation by the cover material [5 M&#x03A9; in <xref ref-type="bibr" rid="B5">Bonmassar et al. (2012)</xref>]. In addition, we applied a long pulse (150 ms) to the coil and measured the electric voltage in the petri dish. If the coil was perfectly covered and there was no leaking current, we recorded a quick and large voltage change at the onset and offset of the pulse (<xref ref-type="fig" rid="F2">Figure 2</xref>). If leakage occurred and the terminal end was exposed to the saline, we recorded a much greater voltage at the onset and offset of the pulse. We also recorded a large, non-zero voltage during the whole stimulus pulse period, due to the leakage current. This method provides an easy way to check the encapsulation and isolation of the coil by the cover material during an electrophysiology experiment.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p>Waveform of the induced electric field generated by a figure-eight micro-coil. Positive voltage pulses (0.2 V) of various durations (1.0, 3.5, and 10 ms) were used to drive the figure-eight micro-coil. The biphasic electric field was generated at the onset and offset phases of the pulse due to the quick change in the coil current.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g002.tif"/>
</fig>
</sec>
<sec id="S2.SS4">
<title>2.4. Activating function for estimating the site of activation by the figure-eight micro-coil</title>
<p>The component of the electric field gradient along the axon, or the activating function (<xref ref-type="bibr" rid="B42">Rattay, 1989</xref>), represents the driving force for activation of the axon. It defines the location and speed of depolarization or hyperpolarization by the extracellular stimulation (<xref ref-type="bibr" rid="B41">Rattay, 1986</xref>; <xref ref-type="bibr" rid="B27">Lee and Fried, 2017</xref>). In &#x03BC;MS technology, the activating function has been a powerful tool in estimating the site of axonal activation (<xref ref-type="bibr" rid="B28">Lee et al., 2016</xref>). When deriving the activating function for the figure-eight micro-coil, we considered the biphasic waveform of the induced electric field, since both the onset and offset phases of the stimulation pulse could play essential roles in axonal activation.</p>
<p>For the onset of a positive pulse, the activating function (gradient of electric field) in the direction of the axon was,</p>
<disp-formula id="S2.Ex22">
<label>(10-1)</label>
<mml:math id="M22">
<mml:mrow>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:mi>F</mml:mi>
</mml:mpadded>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mpadded width="+3.3pt">
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mpadded>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mi/>
</mml:mrow>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mo>&#x007B;</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mfrac>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mfrac>
</mml:mrow>
<mml:mo>&#x007D;</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>For the offset of the pulse, the activating function (gradient of electric field) in the direction of the axon was,</p>
<disp-formula id="S2.Ex23">
<label>(10-2)</label>
<mml:math id="M23">
<mml:mrow>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:mi>F</mml:mi>
</mml:mpadded>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mpadded width="+3.3pt">
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mpadded>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mi/>
</mml:mrow>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mo>&#x007B;</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mfrac>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mfrac>
</mml:mrow>
<mml:mo>&#x007D;</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Equations were derived and result figures were illustrated with Mathematica (version 12.3, Wolfram).</p>
</sec>
<sec id="S2.SS5">
<title>2.5. Combining the figure-eight micro-coil model with the axonal model in NEURON simulation</title>
<p>For NEURON simulation, we need to provide NEURON with the electric potential distribution along the modeled axon. This was calculated by integrating the scalar component of the electric field (Eqs. 6-1, 9-1) along the path of the axon. For the onset of a positive stimulation pulse, the electric potential distribution along the axon was,</p>
<disp-formula id="S2.Ex24">
<label>(11-1)</label>
<mml:math id="M24">
<mml:mrow>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo rspace="5.8pt" stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mo largeop="true" symmetric="true">&#x222B;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo mathvariant="italic" rspace="0pt">d</mml:mo>
<mml:mpadded width="+3.3pt">
<mml:mi>x</mml:mi>
</mml:mpadded>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E25">
<mml:math id="M25">
<mml:mrow>
<mml:mn>&#x2003;&#x2003;&#x2003;&#x2003;&#x2003;&#x2003;&#x2003;&#x2003;&#x2003;</mml:mn>
<mml:mtext>[atan</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mi>y</mml:mi>
</mml:mfrac>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>-</mml:mo>
<mml:mtext>atan</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mi>y</mml:mi>
</mml:mfrac>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo stretchy="false">]</mml:mo>
<mml:mi>e</mml:mi>
<mml:msup>
<mml:mi/>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>For the offset of the pulse, the electric potential distribution along the axon was,</p>
<disp-formula id="S2.Ex26">
<label>(11-2)</label>
<mml:math id="M26">
<mml:mrow>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo rspace="5.8pt" stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mo largeop="true" symmetric="true">&#x222B;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo mathvariant="italic" rspace="0pt">d</mml:mo>
<mml:mpadded width="+3.3pt">
<mml:mi>x</mml:mi>
</mml:mpadded>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo lspace="0pt" rspace="2.8pt">=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mn>&#x03BC;</mml:mn>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>N</mml:mi>
<mml:mmultiscripts>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
<mml:none/>
<mml:none/>
<mml:mn>2</mml:mn>
</mml:mmultiscripts>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>L</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E27">
<mml:math id="M27">
<mml:mrow>
<mml:mn>&#x2003;&#x2003;&#x2003;&#x2003;&#x2003;&#x2003;&#x2003;&#x2003;&#x2003;</mml:mn>
<mml:mtext>[atan</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mi>y</mml:mi>
</mml:mfrac>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mtext>atan</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mi>y</mml:mi>
</mml:mfrac>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo stretchy="false">]</mml:mo>
<mml:mi>e</mml:mi>
<mml:msup>
<mml:mi/>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>During NEURON simulation, the figure-eight micro-coil was positioned at the middle point of the axon (Left coil center: <italic>x</italic> = 10,000 &#x03BC;m; Right coil center: <italic>x</italic> = 10,500 &#x03BC;m), with the center of the coil 300 &#x03BC;m away from the axon (<xref ref-type="fig" rid="F1">Figure 1B</xref>). Resting membrane potential was set to be &#x2212;65 mV at the beginning of the simulation. The model was ran at 20 degrees Celsius. 250 ms after the initiation of the simulation, a single pulse (2.5 ms in duration) was sent to the coil for axonal activation. The electric voltages induced by the figure-eight micro-coil were used to create the extracellular stimuli. The waveform was defined as biphasic short pulses (1 ms in duration) with alternating directions, as measured experimentally (<xref ref-type="fig" rid="F2">Figure 2</xref>). We defined a vector to store the waveform for each time step during stimulation. In each step of the simulation, the value of the extracellular potential was updated using &#x201C;e_extracellular&#x201D; at each compartment using the vector class&#x2019; &#x201C;play&#x201D; method (<xref ref-type="bibr" rid="B21">Joucla et al., 2014</xref>). <xref ref-type="supplementary-material" rid="FS2">Supplementary Figure 2</xref> compared the difference between the extracellular potential for the figure-eight micro-coil and the circular micro-coil (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>). Threshold of axonal activation was defined as the least coil voltage (V) that could initiate an action potential in the axon.</p>
<p>The parameters of the inductor provided by the manufacturer were used in this model, including the length of the coil (<italic>l</italic> = 0.5 mm), inductance of the coil (<italic>L</italic> = 100 nH), and resistance of the coil (<italic>R</italic> = 2 &#x03A9;). Each circular coil has a radius of <italic>R<sub>c</sub></italic> = 0.25 mm. &#x03BC;<sub>0</sub> = 4&#x03C0; &#x00D7; 10<sup>&#x2013;7</sup> H/m.</p>
</sec>
</sec>
<sec id="S3" sec-type="results">
<title>3. Results</title>
<sec id="S3.SS1">
<title>3.1. Temporal and spatial characterization of the electric field generated by a figure-eight micro-coil</title>
<p>Under time-varying magnetic stimulation with a magnetic coil, neural tissue is activated by the magnetically induced electric field <italic>via</italic> electromagnetic induction (<xref ref-type="bibr" rid="B57">Walsh and Pascual-Leone, 2003</xref>; <xref ref-type="bibr" rid="B62">Ye and Steiger, 2015</xref>). The temporal and spatial properties of the induced electric field are directly responsible for the activation of neural tissue with the micro-coil technology (<xref ref-type="bibr" rid="B5">Bonmassar et al., 2012</xref>; <xref ref-type="bibr" rid="B26">Lee and Fried, 2014</xref>). In most experimental settings, a single pulse or trains of electric pulses were delivered to the coil for neural activation. Due to electromagnetic induction theory, the electric field will be induced at the onset and offset of the stimulus pulse, when coil current is dramatically changed in a short period of time. Therefore, the electric field induced by the micro-coils are normally biphasic in shape, as has been measured by us (<xref ref-type="bibr" rid="B61">Ye and Barrett, 2021</xref>; <xref ref-type="bibr" rid="B60">Ye, 2022</xref>) and others (<xref ref-type="bibr" rid="B27">Lee and Fried, 2017</xref>; <xref ref-type="bibr" rid="B11">Golestanirad et al., 2018</xref>; <xref ref-type="bibr" rid="B32">Minusa et al., 2018</xref>).</p>
<p>For the figure-eight micro-coil, we have derived the analytical expression of the induced electric field (Eqs. 9-1, 9-2). Similar to a circular micro-coil, if a voltage pulse were to be delivered to the figure-eight micro-coil, the induced electric field is also a biphasic signal, triggered at the onset and offset of the single stimulus pulse. To validate this model prediction, we designed and implemented a figure-eight micro-coil using two commercial inductors that have been previously used for micromagnetic stimulation of various neural tissues (<xref ref-type="bibr" rid="B5">Bonmassar et al., 2012</xref>; <xref ref-type="bibr" rid="B34">Park et al., 2013</xref>; <xref ref-type="bibr" rid="B63">Ye et al., 2020</xref>; <xref ref-type="bibr" rid="B61">Ye and Barrett, 2021</xref>). The internal structure of the coil (<xref ref-type="supplementary-material" rid="FS1">Supplementary Figure 1</xref>) was revealed by chemically dissolving the cover of the inductor (<xref ref-type="bibr" rid="B63">Ye et al., 2020</xref>). We measured the induced electric field around the micro-coil using an electrophysiology setup. When a pulse signal of varying duration was delivered to the micro-coil, the shape of the induced electric field was biphasic in shape (<xref ref-type="fig" rid="F2">Figure 2</xref>).</p>
<p>Due to its submillimeter scale, it was difficult to measure the spatial distribution of the induced electric field generated a micro-coil with conventional electrodes. We, therefore, computed and visualized the electric field generated by the figure-eight micro-coil (<xref ref-type="fig" rid="F3">Figure 3</xref>). <xref ref-type="fig" rid="F3">Figure 3A</xref> illustrates the electric field generated by the onset and offset of a positive stimulus pulse, respectively. <xref ref-type="fig" rid="F3">Figure 3B</xref> illustrates the field generated by the onset and offset of a negative stimulus pulse, respectively. The coil currents in each of the two circular micro-coils have opposite directions. The intensity of the electric field was the summation of the two individual, circular micro-coils. This generated a large local field around the middle line of the figure-eight micro-coil. This field decayed quickly with increased distance from the coil.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p>Spatial distribution of the induced electric field generated by a figure-eight micro-coil at the onset and offset of the stimulus pulse, respectively. Vector plot and contour plot were illustrated. The dashed black line represented the axon, located 300 &#x03BC;m away from the center of the coil. <bold>(A)</bold> Electric field generated by a positive voltage pulse in the coil. Coil current (I) was counterclockwise in the L-coil, and clockwise in the R-coil. The coil current was dramatically increased at the onset and decreased at the offset of the positive pulse. <bold>(B)</bold> Electric field generated by a negative voltage pulse in the coil. Coil current (I) was clockwise in the L-coil, and counterclockwise in the R-coil. The coil current dramatically increased at the onset and decreased at the offset of the negative pulse.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g003.tif"/>
</fig>
</sec>
<sec id="S3.SS2">
<title>3.2. Activating function predicts the location of membrane polarization and its dependency on coil parameters</title>
<p>The activating function is defined as the gradients of the electric field along the axon (<xref ref-type="bibr" rid="B41">Rattay, 1986</xref>). It predicts the location and speed of depolarization or hyperpolarization by the extracellular stimulation (<xref ref-type="bibr" rid="B28">Lee et al., 2016</xref>, <xref ref-type="bibr" rid="B29">2019</xref>; <xref ref-type="bibr" rid="B27">Lee and Fried, 2017</xref>). The activating function generated by the figure-eight micro-coil was derived in (Eqs. 10-1, 10-2). <xref ref-type="fig" rid="F4">Figure 4A</xref> illustrates the activating function (<italic>dEx/dx</italic>) generated by the onset and offset of a positive pulse, respectively. <xref ref-type="fig" rid="F4">Figure 4B</xref> illustrates the activating function (<italic>dEx/dx</italic>) generated by the onset and offset of a negative pulse, respectively. The gradient of the induced electric field is at its maximum around the midline between the two coils and at the flank of the individual coils with opposite sign.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption><p>Gradient of the induced electric field (activating function) generated by the figure-eight micro-coil in the axon direction (<italic>dEx/dx</italic>), and its dependency on the coil parameters. Vector plot and contour plot were illustrated. The dashed black line represented the axon, located 300 &#x03BC;m away from the center from the coil. <bold>(A)</bold> Activating function generated by a positive voltage pulse in the coil. Coil current (I) was counterclockwise in the L-coil, and clockwise in the R-coil. The coil current was dramatically increased at the onset and decreased at the offset of the positive pulse. <bold>(B)</bold> Activating function generated by a negative voltage pulse in the coil. Coil current (I) was clockwise in the L-coil, and counterclockwise in the R-coil. The coil current was dramatically increased at the onset and decreased at the offset of the negative pulse.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g004.tif"/>
</fig>
<p>Previously, we have calculated the activating function of a circular micro-coil (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>). This activating function allowed us to analyze the key factors that could cause shifting of activating location under circular micro-coil stimulation (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>). Here, we took a similar approach to understand the location of axonal activation by a figure-eight micro-coil.</p>
<p>For a positive pulse, at its onset, the figure-eight micro-coil generated an activating function that contained a large &#x201C;hyperpolarization peak&#x201D; (virtual cathode) between the two coils (<italic>x</italic> = 0), and two smaller &#x201C;depolarization peaks&#x201D; (virtual anode, <xref ref-type="fig" rid="F5">Figure 5A1</xref>). At the offset, the direction of the induced electric field was reversed, leading to the switch of the depolarization/hyperpolarization sites on the axon (<xref ref-type="fig" rid="F5">Figure 5A2</xref>). For a negative pulse, at its onset, the figure-eight micro-coil generated an activating function that contained a large &#x201C;depolarization peak&#x201D; (virtual anode) between the two coils (<italic>x</italic> = 0), and two smaller &#x201C;hyperpolarization peaks&#x201D; (virtual cathode, <xref ref-type="fig" rid="F5">Figure 5B1</xref>). At the offset, the direction of the induced electric field was reversed, leading to the switch of the depolarization/hyperpolarization sites on the axon (<xref ref-type="fig" rid="F5">Figure 5B2</xref>). To estimate the impact of changes in the coil-axon distance on the site of activation, we increased the coil-axon distance from 300 to 800 &#x03BC;m. This led to a much smaller measurement of the activating function (<xref ref-type="fig" rid="F5">Figures 5A3</xref>, <xref ref-type="fig" rid="F5">5B3</xref>). However, the large peaks in the activating functions did not change their locations.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption><p>Activating function along the axon and its dependency on the coil parameters. <bold>(A)</bold> Activating function along the axon when a positive pulse was delivered into the figure-eight micro-coil. (1) Activating function at the onset of the positive pulse, when the coil-axon distance was 300 &#x03BC;m. (2) Activating function at the offset of the positive pulse, when the coil-axon distance was 300 &#x03BC;m. (3) Activating function at the onset of the positive pulse, when the coil-axon distance was 800 &#x03BC;m. <bold>(B)</bold> Activating function along the axon when a negative pulse was delivered into the figure-eight micro-coil. (1) Activating function at the onset of the negative pulse, when the coil-axon distance was 300 &#x03BC;m. (2) Activating function at the offset of the negative pulse, when the coil-axon distance was 300 &#x03BC;m. (3) Activating function at the onset of the negative pulse, when the coil-axon distance was 800 &#x03BC;m.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g005.tif"/>
</fig>
<p>These analyses allow us to make the following four predictions about axonal activation with a figure-eight micro-coil: (1) Since the sharp, larger depolarization peaks on the axon under figure-eight micro-coil stimulation are generated by the summation effects of two individual, circular micro-coils, the figure-eight micro-coil will need less coil voltage to initiate action potentials than a circular micro-coil. (2) Reversal of the stimulation pulse in the figure-eight micro-coil will reverse the depolarization/hyperpolarization sites, which may produce different mechanisms of axonal activation. (3) Increasing the distance between the figure-eight micro-coil and the axon will not change the activation site. (4) Increasing stimulus intensity will not cause shifting of the activation site. In the following sections, we will use a multi-compartment model of an unmyelinated axon to test these predictions in NEURON.</p>
</sec>
<sec id="S3.SS3">
<title>3.3. NEURON modeling of axonal activation with a figure-eight micro-coil demonstrates two types (A and B) of axonal activation, caused by the reversal of the coil current</title>
<p>A multi-compartment NEURON model of an unmyelinated axon (<xref ref-type="fig" rid="F1">Figure 1B</xref>) was built to test the accuracy of these predictions. The figure-eight micro-coil is located in the middle segment of the axon. The distance between the center of the coil (O<sub>1</sub> or O<sub>2</sub>) and the axon is <italic>y</italic> = 300 &#x03BC;m. In NEURON simulation, a single pulse (2.5 ms) was sent to the coil model. The model computed the voltage distribution along the axon (Eqs. 11-1, 11-2) and applied this information to the NEURON simulation environment (<xref ref-type="supplementary-material" rid="FS2">Supplementary Figure 2</xref>). This externally applied voltage triggered an action potential.</p>
<p>Previously, we have investigated the activation of the unmyelinated axon by a circular micro-coil (<xref ref-type="fig" rid="F1">Figure 1A</xref>). We found two distinct ways in which the circular micro-coil could trigger an action potential in the axon (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>). In Type I activation, the onset phase preconditioned the membrane for activation, and the offset phase of the stimulus pulse could trigger an action potential [Figure 6 in <xref ref-type="bibr" rid="B60">Ye (2022)</xref>]. This type of activation needs less stimulus intensity, with an activating function of 250,000 V/m<sup>2</sup>. For Type II activation, the onset phase of the stimulus pulse could directly trigger an action potential [Figure 8 in <xref ref-type="bibr" rid="B60">Ye (2022)</xref>]. This type of axonal activation requires a relatively larger stimulus intensity.</p>
<p>For the figure-eight micro-coil, depending on the direction of the coil current, we also observed two different types of axonal activation (Type A and Type B activations).</p>
<sec id="S3.SS3.SSS1">
<title>3.3.1. Type A activation by a figure-eight micro-coil</title>
<p><xref ref-type="fig" rid="F6">Figure 6A</xref> demonstrated a positive pulse applied to the figure-eight micro-coil, the induced electric field, and the location of axonal activation. <xref ref-type="fig" rid="F6">Figure 6B</xref> simulated the initiation and traveling of an action potential with a sequence of frames (1 ms intervals). The onset of the pulse generated two depolarization peaks and one larger hyperpolarization peak (<italic>t</italic> = 1 ms in <xref ref-type="fig" rid="F6">Figure 6B</xref>), as predicted with activating function analysis (<xref ref-type="fig" rid="F5">Figure 5A1</xref>). The small depolarization peaks failed to trigger an action potential. Instead, an action potential was triggered during the subsequent offset phase of the coil current (<italic>t</italic> = 5 ms, <xref ref-type="fig" rid="F6">Figure 6B</xref>), by depolarizing the membrane patch where the hyperpolarization peak presented during the onset phase. We refer to this activation pattern as Type A activation. In Type A activation, axonal activation happens at the membrane patch between the two circular micro-coils, due to membrane hyperpolarization followed by depolarization. The onset of the positive pulse conditions the membrane, and the axonal activation is eventually triggered by the offset phase. The action potential traveled bi-directionally to each end of the axon. The threshold for Type A activation by the figure-eight micro-coil is 83.5% of that required for Type I activation by the circular micro-coil (confirms prediction 1).</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption><p>Activation of an unmyelinated axon by the figure-eight micro-coil with a positive pulse (Type A activation). The action potential was triggered by the offset phase of a positive stimulation pulse. <bold>(A)</bold> Coil current, direction of the induced electric field, and location of activation (+). The coil center was located 300 &#x03BC;m away from the axon. <bold>(B)</bold> A sequence of movie clips demonstrated the initiation and propagation of an action potential. Onset of the pulse generated a local hyperpolarization of the membrane, and the offset phase of the pulse generated a local depolarization, which triggered an action potential. The action potential was initiated at the midline between the two coils, and propagated in both directions. The green box, which shows the location of activation, will be used for comparison in the following <xref ref-type="fig" rid="F7">Figures 7</xref>, <xref ref-type="fig" rid="F9">9</xref>, <xref ref-type="fig" rid="F1">11</xref>.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g006.tif"/>
</fig>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption><p>Activation of an unmyelinated axon by the figure-eight micro-coil (Type B activation) with a negative pulse. The action potential was triggered by the onset phase of a negative stimulation pulse. <bold>(A)</bold> Coil current, direction of the induced electric field, and location of activation (+). The coil center was located 300 &#x03BC;m away from the axon. <bold>(B)</bold> A sequence of movie clips demonstrated the initiation and propagation of an action potential. The action potential was initiated by the onset phase of the coil current, at the midline between two circular micro-coils, and propagated in both directions. The red box, which shows the location of activation, will be used for comparison in the following <xref ref-type="fig" rid="F10">Figures 10</xref>, <xref ref-type="fig" rid="F12">12</xref>. <bold>(C)</bold> Comparison of locations of activation in Type A activation (<xref ref-type="fig" rid="F6">Figure 6</xref>, green box) and Type B activation (red box) revealed no shifting of the activation sites.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g007.tif"/>
</fig>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption><p>Comparison of ion channel dynamics for action potential initiation by Type A and Type B activations with a figure-eight micro-coil in an unmyelinated axon. Membrane potential (1), Na<sup>+</sup> current (2), K<sup>+</sup> current (3), Na<sup>+</sup> channel activation m (4), Na<sup>+</sup> channel inactivation h (5), and K<sup>+</sup> channel activation n (6) were plotted at the locations where action potentials were initiated [axon (0.5125)]. <bold>(A)</bold> For Type A activation, the induced electric field during the onset phase hyperpolarized the cell membrane (thin arrow in A1) and de-inactivated the sodium channels (arrow in A5). Action potential was then initiated in the offset phase, which depolarized the membrane (thick arrow in A1). <bold>(B)</bold> For Type B activation, the induced electric field during the onset phase initiated the action potential simply by depolarizing the axonal membrane (thin arrow in B1) and activating the sodium channels. The offset phase did slightly hyperpolarize the membrane (thick arrow in B1), but it did not significantly affect the initiation and propagation of the action potential.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g008.tif"/>
</fig>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption><p>Impact of axon-coil distance on the location of Type A activation with a figure-eight micro-coil. <bold>(A)</bold> Coil current, direction of the induced electric field, and location of activation (+). The coil-axon distance was increased from 300 &#x03BC;m (in <xref ref-type="fig" rid="F6">Figure 6</xref>) to 800 &#x03BC;m. Stimulation intensity was adjusted to trigger the action potential. <bold>(B)</bold> A sequence of movie clips demonstrated the initiation and propagation of an action potential. The blue box showed the location of activation. <bold>(C)</bold> Comparison of locations of Type A activation when the axon-coil distance was 300 &#x03BC;m (<xref ref-type="fig" rid="F6">Figure 6</xref>, green box) and 800 &#x03BC;m (blue box). Increase in coil-axon distance did not significantly change the location of activation in Type A activation.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g009.tif"/>
</fig>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption><p>Impact of axon-coil distance on the location of activation in Type B activation with a figure-eight micro-coil. <bold>(A)</bold> Coil current, direction of the induced electric field, and location of activation (+). The coil-axon distance was increased from 300 &#x03BC;m (<xref ref-type="fig" rid="F7">Figure 7</xref>) to 800 &#x03BC;m. Stimulation intensity was adjusted to trigger action potentials. <bold>(B)</bold> A sequence of movie clips demonstrated the initiation and propagation of an action potential. The purple box showed the location of activation. <bold>(C)</bold> Comparison of locations of Type B activation when the axon-coil distance was 300 &#x03BC;m (<xref ref-type="fig" rid="F7">Figure 7</xref>, red box) and the 800 &#x03BC;m (purple box). Increase in coil-axon distance did not significantly change the location of activation in Type B activation.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g010.tif"/>
</fig>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption><p>Impact of stimulation intensity on the location of Type A activation. <bold>(A)</bold> A strong Type A activation with 200% threshold intensity was applied to the axon. A sequence of movie clips demonstrated the initiation and propagation of the action potential. The black box showed the location of activation. <bold>(B)</bold> Comparison of the locations of activation with threshold stimulation (<xref ref-type="fig" rid="F6">Figure 6</xref>, green box) and the black box when the stimulation intensity was doubled. Increasing stimulation intensity caused a quicker axonal activation, but it did not change the location of activation in Type A activation.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g011.tif"/>
</fig>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption><p>Impact of stimulation intensity on the location of Type B activation. <bold>(A)</bold> A strong Type B activation with 200% threshold intensity was applied to the axon. A sequence of movie clips demonstrated the initiation and propagation of the action potential. The yellow box showed the location of activation. <bold>(B)</bold> Comparison of the locations of activation with threshold stimulation (<xref ref-type="fig" rid="F7">Figure 7</xref>, red box) and the yellow box when the stimulation intensity was doubled. Increasing stimulation intensity caused a faster activation, but it did not change location of activation in Type B activation.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-17-1105505-g012.tif"/>
</fig>
</sec>
<sec id="S3.SS3.SSS2">
<title>3.3.2. Type B activation by a figure-eight micro-coil</title>
<p><xref ref-type="fig" rid="F7">Figure 7A</xref> demonstrated a negative pulse applied to the figure-eight micro-coil, the induced electric field, and the location of axonal activation. <xref ref-type="fig" rid="F7">Figure 7B</xref> simulated the initiation and traveling of an action potential with a sequence of frames (1 ms intervals). The onset of the pulse generated two small hyperpolarization peaks and one larger depolarization peak (<italic>t</italic> = 1 ms in <xref ref-type="fig" rid="F7">Figure 7B</xref>), as predicated with activating function analysis (<xref ref-type="fig" rid="F5">Figure 5B1</xref>). The large depolarization peak triggered an action potential during the onset phase of the coil current (<italic>t</italic> = 3 ms, <xref ref-type="fig" rid="F7">Figure 7B</xref>). The offset phase of the pulse generated a local hyperpolarization (<italic>t</italic> = 4 ms, <xref ref-type="fig" rid="F7">Figure 7B</xref>), but it did not affect the action potential. In Type B activation, axonal activation occurred at the membrane patch between the two circular micro-coils due to a strong depolarization. The action potential traveled bi-directionally to each end of the axon. We refer to this activation pattern as Type B activation, in which the onset phase of a negative pulse directly activates the axon by membrane depolarization. Since Type B activation is triggered by the onset of the stimulation pulse, axonal activation was faster in Type B activation than in Type A activation. The threshold for Type B activation is slightly (3.7%) higher than Type A activation. The threshold for Type B activation by the figure-eight micro-coil is 46.8% of that required for Type II activation by the circular micro-coil (confirms prediction 1).</p>
<p>In conclusion, reversal of the stimulation pulse in the figure-eight micro-coil will reverse the depolarization/hyperpolarization sites, which produces two distinct types of axonal activation (confirms prediction 2). However, reversal of pulse polarity did not cause the activation site to shift. The axon was activated at the same location in both Type A and Type B activations (<xref ref-type="fig" rid="F7">Figure 7C</xref>), which ensures that the consistency of stimulation is not affected by accidental reversal of coil current.</p>
</sec>
</sec>
<sec id="S3.SS4">
<title>3.4. Distinct ion channel mechanisms underlying type A vs. type B axonal activation by a figure-eight micro-coil</title>
<p>To further investigate the ionic mechanism underlying action potential initiation in Type A and Type B activations with the figure-eight micro-coil, we studied the membrane dynamics at the locations where action potentials were initiated. We compared the inward sodium current (INa+), outward potassium current (IK+), sodium channel activation (m) and inactivation (h) variables, and potassium channel activation (n) variables. The center of the two coils were at (10,000 &#x03BC;m, &#x2212;300 &#x03BC;m) and (10,500, &#x2212;300 &#x03BC;m), respectively. We found that the largest depolarization happened at <italic>x</italic> = 10, 250 &#x03BC;m, which was the midpoint between the two micro-coils. In NEURON, this point is Axon (0.5125). We, therefore, plotted the channel dynamics at this point in both Type A and B activations. Type A and Type B activations with the figure-eight micro-coil demonstrated different ionic mechanisms for triggering action potentials.</p>
<p><xref ref-type="fig" rid="F8">Figure 8A</xref> demonstrated Type A activation, in which the membrane patch experienced a brief hyperpolarization during the onset phase of the positive stimulation pulse, followed by a depolarization induced by the offset phase of the stimulation pulse, which eventually triggered an action potential (<xref ref-type="fig" rid="F8">Figure 8A1</xref>). The initial hyperpolarization did not cause significant changes in the state variable <italic>m</italic> (around 0, <xref ref-type="fig" rid="F8">Figure 8A4</xref>) and <italic>n</italic> (around 0.3, <xref ref-type="fig" rid="F8">Figure 8A6</xref>) values. However, the hyperpolarization removed sodium channel inactivation and caused a clear increase in the <italic>h</italic> value (from 0.7 to 0.8, <xref ref-type="fig" rid="F8">Figure 8A5</xref>). This de-inactivation of the sodium channels led to an increase in the sodium conductance, defined as m3h (<xref ref-type="bibr" rid="B15">Hodgkin and Huxley, 1952</xref>). The large <italic>h</italic> value was maintained after the termination of the onset phase until the offset phase, which caused a membrane depolarization, leading to a large inward sodium current (<xref ref-type="fig" rid="F8">Figure 8A2</xref>) and action potential initiation (<xref ref-type="fig" rid="F8">Figure 8A1</xref>). Due to the hyperpolarization during the onset phase, a larger amount of sodium channels were preconditioned for activation during the offset phase, which is why less stimulation intensity was needed in Type A activation than in Type B activation.</p>
<p><xref ref-type="fig" rid="F8">Figure 8B</xref> demonstrated Type B activation, in which the membrane patch experienced a large depolarization during the onset phase of the stimulation pulse (<xref ref-type="fig" rid="F8">Figure 8B1</xref>). The depolarization mainly activated the sodium channels by increasing the <italic>m</italic> value (from 0 to 0.3, <xref ref-type="fig" rid="F8">Figure 8B4</xref>). This increase in sodium channel activation led to the initiation of an action potential. In comparison with Type A activation (<xref ref-type="fig" rid="F8">Figure 8A5</xref>), this depolarization did not alter the <italic>h</italic> value (<xref ref-type="fig" rid="F8">Figure 8B5</xref>) to de-inactivate more sodium channels. The following offset phase introduced a brief membrane hyperpolarization; however, it did not affect the action potential.</p>
<p>In summary, under Type A activation, the local membrane was hyperpolarized during the onset phase of the coil current, which de-inactivated the sodium channels. This same membrane patch was then activated by the offset phase of the coil current, and an action potential was initiated. In contrast, during Type B activation (which requires greater coil current than Type A activation), the action potential was generated at the same location, where the membrane was directly depolarized by the onset phase of the coil current.</p>
</sec>
<sec id="S3.SS5">
<title>3.5. Activation location does not change with varying coil-axon distance</title>
<p>Previously, when modeling the circular micro-coil, we found that changing the coil-axon distance could cause shifting of the activation site [Figure 7 in <xref ref-type="bibr" rid="B60">Ye (2022)</xref>]. Activating function analysis for the figure-eight micro-coil, in contrast, suggests that the site of membrane depolarization/hyperpolarization is not affected by the increase in coil-axon distance in Type A (<xref ref-type="fig" rid="F5">Figure 5A3</xref>) or Type B (<xref ref-type="fig" rid="F5">Figure 5B3</xref>) activations.</p>
<p>In NEURON simulation, we increased the coil-axon distance from 300 to 800 &#x03BC;m and pinpointed the site where the action potential was initiated. For both Type A (<xref ref-type="fig" rid="F9">Figures 9A, B</xref>) and Type B (<xref ref-type="fig" rid="F10">Figures 10A, B</xref>) activations, a greater stimulation intensity was required to trigger an action potential. However, increasing the coil-axon distance caused no shifting of the activation site in Type A activation (<xref ref-type="fig" rid="F9">Figure 9C</xref>). Similarly, increasing the coil-axon distance caused no shifting of the activation site in Type B activation (<xref ref-type="fig" rid="F10">Figure 10C</xref>) (confirms prediction 3).</p>
</sec>
<sec id="S3.SS6">
<title>3.6. Location of activation does not change with varying stimulation intensities</title>
<p>Previously, when modeling the single, circular micro-coil, we found that the activation site was dependent on the intensity of coil stimulation [Figure 8 in <xref ref-type="bibr" rid="B60">Ye (2022)</xref>]. Activating function analysis for the figure-eight micro-coil, in contrast, suggests that the site of membrane depolarization/hyperpolarization was not affected by the intensity of the stimulus (Eqs. 10-1, 10-2). We hypothesize that increasing the stimulus intensity will not cause any shifting in the activation sites for both Type A and Type B activation.</p>
<p>In <xref ref-type="fig" rid="F11">Figure 11</xref>, for Type A activation, we doubled the stimulation intensity in the figure-eight micro-coil. Onset of the pulse generated a local hyperpolarization, followed by a depolarization and the initiation of the action potential, which traveled in both directions along the axon (<xref ref-type="fig" rid="F11">Figure 11A</xref>). In comparison to <xref ref-type="fig" rid="F6">Figure 6</xref>, local depolarization and hyperpolarization were more prominent when the stimulus intensity doubled, leading to a quick initiation of the action potential. However, the location of activation was not shifted by the intense stimulation in <xref ref-type="fig" rid="F11">Figure 11B</xref> (confirms prediction 4).</p>
<p>For Type B activation, we doubled the stimulation intensity in the figure-eight micro-coil. Onset of the pulse generated a large local depolarization and the initiation of the action potential, which traveled in both directions along the axon (<xref ref-type="fig" rid="F12">Figure 12A</xref>). In comparison to <xref ref-type="fig" rid="F7">Figure 7</xref>, this local depolarization was more prominent when the stimulus intensity was doubled, leading to a quicker initiation of the action potential. However, the location of activation was not shifted by the intense stimulation in <xref ref-type="fig" rid="F12">Figure 12B</xref> (confirms prediction 4).</p>
<p>In summary, NEURON simulation confirms the four predictions made in the activating function analysis. With less coil current, the figure-eight micro-coil provides a more focal activation than the circular micro-coil. This activation site does not shift by the reversal of the coil current, although it will result in different mechanisms of axonal activation. The location of axonal activation by the figure-eight micro-coil is not affected by the intensity of coil current or the translational movement of the coil away from the axon. These results suggest that the figure-eight micro-coil could provide significantly improved consistency over the circular micro-coil in axonal activation.</p>
</sec>
</sec>
<sec id="S4" sec-type="discussion">
<title>4. Discussion</title>
<p>One of the major goals of magnetic stimulation with a micro-coil is to provide consistent stimulation of individual neurons. Previously, we have investigated the consistency of neural simulation with a circular micro- coil, using combined biophysics modeling and NEURON simulation (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>). This paper further examines the consistency of micro-coil stimulation with a figure-eight micro-coil, and compares it with the circular micro-coil. The figure-eight micro-coil provides significantly improved consistency in neural stimulation. <xref ref-type="table" rid="T1">Table 1</xref> compares axonal activation by a figure-eight micro-coil and a circular micro-coil [detailed in <xref ref-type="bibr" rid="B60">Ye (2022)</xref>], respectively.</p>
<table-wrap position="float" id="T1">
<label>TABLE 1</label>
<caption><p>Comparison between the circular micro-coil and figure-eight micro-coil in axon activation.</p></caption>
<table cellspacing="5" cellpadding="5" frame="box" rules="all">
<thead>
<tr>
<td valign="top" align="left" style="color:#ffffff;background-color: #7f8080;">Coil type</td>
<td valign="top" align="left" style="color:#ffffff;background-color: #7f8080;">Circular</td>
<td valign="top" align="left" style="color:#ffffff;background-color: #7f8080;">Circular</td>
<td valign="top" align="left" style="color:#ffffff;background-color: #7f8080;">Figure&#x2013;eight</td>
<td valign="top" align="left" style="color:#ffffff;background-color: #7f8080;">Figure&#x2013;eight</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left"><bold>Type of activation</bold></td>
<td valign="top" align="left">Type I</td>
<td valign="top" align="left">Type II</td>
<td valign="top" align="left">Type A</td>
<td valign="top" align="left">Type B</td>
</tr>
<tr>
<td valign="top" align="left"><bold>Stimulus signal</bold></td>
<td valign="top" align="left">Single pulse (positive or negative)</td>
<td valign="top" align="left">Single pulse (positive or negative)</td>
<td valign="top" align="left">Single pulse (positive)</td>
<td valign="top" align="left">Single pulse (negative)</td>
</tr>
<tr>
<td valign="top" align="left"><bold>Waveform of induced electric field</bold></td>
<td valign="top" align="left">Biphasic (onset and offset phases)</td>
<td valign="top" align="left">Biphasic (onset and offset phases)</td>
<td valign="top" align="left">Biphasic (onset and offset phases)</td>
<td valign="top" align="left">Biphasic (onset and offset phases)</td>
</tr>
<tr>
<td valign="top" align="left"><bold>Action potential triggering time</bold></td>
<td valign="top" align="left">Offset phase</td>
<td valign="top" align="left">Onset phase</td>
<td valign="top" align="left">Offset phase</td>
<td valign="top" align="left">Onset phase</td>
</tr>
<tr>
<td valign="top" align="left"><bold>Membrane dynamic responsible for action potential initiation</bold></td>
<td valign="top" align="left">Hyperpolarization&#x2013;depolarization</td>
<td valign="top" align="left">Depolarization</td>
<td valign="top" align="left">Hyperpolarization&#x2013;depolarization</td>
<td valign="top" align="left">Depolarization</td>
</tr>
<tr>
<td valign="top" align="left"><bold>Ion channel dynamics underlying action potential initiation</bold></td>
<td valign="top" align="left">Na channel de-inactivation followed by activation</td>
<td valign="top" align="left">Na channel activation</td>
<td valign="top" align="left">Na channel de-inactivation followed by activation</td>
<td valign="top" align="left">Na channel activation</td>
</tr>
<tr>
<td valign="top" align="left"><bold>Activation site:</bold><break/> <bold>Coil current direction dependent</bold></td>
<td valign="top" align="left">Yes</td>
<td valign="top" align="left">Yes</td>
<td valign="top" align="left">No</td>
<td valign="top" align="left">No</td>
</tr>
<tr>
<td valign="top" align="left"><bold>Activation site:</bold><break/> <bold>Coil axon distance dependent</bold></td>
<td valign="top" align="left">Yes</td>
<td valign="top" align="left">Yes</td>
<td valign="top" align="left">No</td>
<td valign="top" align="left">No</td>
</tr>
<tr>
<td valign="top" align="left"><bold>Activation site:</bold><break/> <bold>Coil current intensity dependent</bold></td>
<td valign="top" align="left">Yes</td>
<td valign="top" align="left">Yes</td>
<td valign="top" align="left">No</td>
<td valign="top" align="left">No</td>
</tr>
<tr>
<td valign="top" align="left"><bold>Threshold of activation</bold><break/> <bold>a. Same coil comparison</bold><break/> <bold>b. Different coil comparison</bold></td>
<td valign="top" align="left">250,000 V/m<sup>2</sup> (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>).<break/> &#x2212;<break/>&#x2212;</td>
<td valign="top" align="left">184.6% of Type I<break/> &#x2212;</td>
<td valign="top" align="left">&#x2212;<break/> 83.5% of Type I</td>
<td valign="top" align="left">103.7% of Type A<break/> 46.8% of Type II</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn><p>For detailed analysis of the circular micro-coil, see reference <xref ref-type="bibr" rid="B60">Ye (2022)</xref>.</p></fn>
</table-wrap-foot>
</table-wrap>
<p>First, a single stimulus pulse delivered to either of these coils generates a biphasic electric field, due to the quick coil current changes at the onset and offset phases of the stimulus pulse. The biphasic waveform is a favored waveform when an electrode is used to deliver currents for neural blockage in peripheral nerves (<xref ref-type="bibr" rid="B6">Cattell and Gerard, 1935</xref>; <xref ref-type="bibr" rid="B52">Tanner, 1962</xref>). When compared to monophasic stimulation, biphasic stimulation causes less tissue damage due to the neutralization properties of electrochemical reactions (<xref ref-type="bibr" rid="B50">Tai et al., 2005</xref>). The intrinsic biphasic nature of the induced electric field suggests that micro-coil stimulation can lead to a null charge integrated over time, and provide biocompatible stimulation (<xref ref-type="bibr" rid="B19">Jeong et al., 2022</xref>).</p>
<p>Second, both the figure-eight micro-coil and circular micro-coil could activate the axon <italic>via</italic> two distinct mechanisms. These activation patterns share many commonalities. In Type I (circular micro-coil) and Type A (figure-eight micro-coil) activations, the membrane is hyperpolarized at the onset phase, and the sodium channels are sufficiently de-inactivated. Action potential initiation occurs during the offset phase, which depolarizes the membrane. In Type II (circular) and Type B (figure-eight) activations, the membrane is depolarized at the onset phase, with sufficiently inactivates the sodium channels and elicits the action potential. As the consequence of ion channel preconditioning, the threshold for Type I (circular) and Type A (figure-eight) activations are significantly lower than those in Type II (circular) and Type B (figure-eight) activations, respectively.</p>
<p>Third, the figure-eight micro-coil provides a much more focal and consistent axonal activation mechanism than the circular micro-coil. In circular micro-coil stimulation, the activation site shifts with the reversal of coil current direction, changes in the coil-axon distance, and increases in the coil current intensity. In figure-eight micro-coil stimulation, however, changes in these experimental factors have a minimal impact on the site of activation.</p>
<p>Finally, stimulation with the figure-eight micro-coil, overall, requires less coil current than the circular micro-coil for axonal activation. Furthermore, activation due to the preconditioning of the sodium channel (Type I and Type A activations) normally require less coil current than activation without such preconditioning (Type II and Type B activations). These results could potentially provide a solution to the heating effects observed in some &#x03BC;MS practices (<xref ref-type="bibr" rid="B61">Ye and Barrett, 2021</xref>).</p>
<sec id="S4.SS1">
<title>4.1. Implication to the development of &#x03BC;MS technology</title>
<p>Results from this work provide several insights to the novel design of micro-coils for efficient neural stimulation with improved consistency.</p>
<p>First, figure-eight micro-coil design could provide a simple design strategy to improve stimulating efficacy. In circular micro-coil stimulation, to increase the intensity of the induced electric field for stimulation, one could include more winding or encapsulate a core with high magnetic permeability. This could potentially increase the size of the coil and provide challenges in implantation and increase the risk of brain damage. This problem could be avoided, however, if a figure-eight micro-coil were to be used for neural stimulation. A recent work simulated the micro-coil stimulation of the radial nerve fibers in arms with a circular micro-coil and a figure-eight micro-coil. The authors found that the threshold current in the figure-eight micro-coil was significantly less than the circular micro-coil for neural activation (<xref ref-type="bibr" rid="B8">Colella et al., 2021</xref>). Our results are generally in agreeance with this report. The figure-eight coil could generate a large local field gradient without the need to significantly increase the winding and coil current.</p>
<p>Second, it is essential to have a detailed understanding of both spatial and temporal properties of the electric field generated by the figure-eight micro-coil. Our results demonstrate that these properties play significant roles in the pattern of membrane depolarization/hyperpolarization (<xref ref-type="fig" rid="F6">Figures 6</xref>, <xref ref-type="fig" rid="F7">7</xref>), which define the site and speed of axonal activation. Various technologies are available to measure the electric field generated by the micro-coil. This includes electrophysiological recording electrodes to measure the field potential (<xref ref-type="bibr" rid="B27">Lee and Fried, 2017</xref>; <xref ref-type="bibr" rid="B11">Golestanirad et al., 2018</xref>), custom-made pickup coils to measure the electromotive force (<xref ref-type="bibr" rid="B32">Minusa et al., 2018</xref>), and nitrogen vacancy in diamond magnetometers to directly measure the magnetic flux density (<xref ref-type="bibr" rid="B22">Khalifa et al., 2021</xref>). The spatial resolution of these measurements could be limited by the size of the sensors. Alternatively, finite element modeling provides an easy estimation of the induced electric field with high resolution (<xref ref-type="bibr" rid="B29">Lee et al., 2019</xref>; <xref ref-type="bibr" rid="B45">Saha et al., 2022a</xref>).</p>
<p>Third, it is essential to design the stimulus waveforms and understand their significance in neural stimulation using &#x03BC;MS technology. Previously, using a circular micro-coil as a model, we have demonstrated the importance of the stimulus waveform in affecting the location of axonal activation (<xref ref-type="bibr" rid="B60">Ye, 2022</xref>). Here, using a figure-eight micro-coil, we demonstrated that, depending on the direction of coil current, axonal activation could be triggered by either the offset of the pulse (Type A, <xref ref-type="fig" rid="F6">Figure 6</xref>) or onset of the pulse (Type B, <xref ref-type="fig" rid="F7">Figure 7</xref>). Therefore, deliberate control of coil current could speed up or delay the onset time of action potentials.</p>
<p>Fourth, it is essential to understand the biophysics mechanisms of neural activation when applying a stimulus. For example, the stimulus threshold for Type B activation is lower than Type A activation. This is because the membrane patch in Type A activation experienced a priming period during the onset phase of the stimulus pulse, leading to the de-inactivation of the sodium channels (<xref ref-type="fig" rid="F8">Figure 8A</xref>). Therefore, when choosing a specific stimulation intensity needed for axonal activation, one must consider the properties of the neurons, and the dynamic interaction between the neuron and the field (<xref ref-type="bibr" rid="B62">Ye and Steiger, 2015</xref>). On the other hand, deliberate control of the stimulus could lead to reduced coil current, and prevent unnecessary heating effects associated with micro-coil technology (<xref ref-type="bibr" rid="B46">Saha et al., 2022b</xref>).</p>
<p>Finally, it is essential to further develop multi-coil technology. Determining the optimal configuration for efficient stimulation has been a long quest in neuroscience, both in fundamental research and for therapeutic development. One method is to combine the source of stimuli (<xref ref-type="bibr" rid="B65">Zaaimi et al., 2013</xref>; <xref ref-type="bibr" rid="B2">Alekseichuk et al., 2019</xref>). For example, in electrode stimulation, experimental strategies aiming to improve the stimulating focality around the activation sites generally combines two (<xref ref-type="bibr" rid="B43">Rattay and Resatz, 2007</xref>), three (<xref ref-type="bibr" rid="B48">Sankarasubramanian et al., 2011</xref>), or multiple electrodes (<xref ref-type="bibr" rid="B65">Zaaimi et al., 2013</xref>). Such multipolar electrodes are usually associated with higher activation thresholds than the monopolar, single electrode configuration. Since tissue damage and electrode deterioration are directly related to the injected electric current <italic>via</italic> the multiple electrodes, electric current intensity must be constrained for safety reasons (<xref ref-type="bibr" rid="B12">Gunter et al., 2019</xref>). In this paper, when two circular micro-coils were combined to construct the figure-eight micro-coil, we could obtain a significant increase in the focality and field gradient with a decreased current in the coil. Future research in &#x03BC;MS technology must explore the possibility of combining multiple micro-coils for highly specific stimulation.</p>
</sec>
<sec id="S4.SS2">
<title>4.2. Implication to neural activation with highly consistent &#x03BC;MS</title>
<p>Recent &#x03BC;MS research focuses on achieving high resolution stimulation at the cellular level. There have been several preparations in which &#x03BC;MS has been successful in neural activation. The goal of these works is to constrain the stimulation site within certain geometrical dimensions of the neural tissue. In order to achieve consistent stimulation results, it is essential to know the dimensions of the target tissue and estimate the allowable range of shifting in the activation sites.</p>
<p>One of the early &#x03BC;MS successes was stimulating retinal ganglion cells (RGCs) (<xref ref-type="bibr" rid="B5">Bonmassar et al., 2012</xref>). These cells can vary in soma size and dendritic field size, stratification, shape, and compactness (<xref ref-type="bibr" rid="B23">Kim et al., 2021</xref>). For example, midget RGCs, which functionally serve to distinguish red and green, have a dendritic field size range of 10&#x2013;100 &#x03BC;m and, overall, have smaller cell bodies and receptive fields. Meanwhile, parasol RGCs, which have the function of movement, have a dendritic field size range of 30&#x2013;300 &#x03BC;m and have larger cell bodies and receptive fields. RGC axons can drastically vary in length from 0.5 &#x03BC;m to 50 mm (<xref ref-type="bibr" rid="B64">Yu et al., 2013</xref>), and are typically myelinated in the optic nerve and unmyelinated in the retinal nerve fiber layer (<xref ref-type="bibr" rid="B59">Yang et al., 2013</xref>). Therefore, the size of the soma, dendritic field, and axon of RGCs as well as the function and myelination of the cell necessitates highly specific stimulation, most typically in the micrometer range. This ultimately requires delicate control with high specificity and consistency, which could potentially be achieved by the figure-eight micro-coil or another combined coil strategy.</p>
<p>&#x03BC;MS was also used to stimulate presynaptic axons in the Schaffer collateral for synaptic transmission in the hippocampus (<xref ref-type="bibr" rid="B45">Saha et al., 2022a</xref>). The synaptic cleft is 20&#x2013;30 nm at synapses between neurons in the central nervous system and about 50 nm at the neuromuscular junction. A shifting of stimulation site could potentially miss the presynaptic site, rather activating the postsynaptic neurons. A figure-eight micro-coil could possibly provide a more reliable stimulation by ensuring the presynaptic axon is the sole target for stimulation. Alternatively, if a circular micro-coil were to be used, to avoid side effects of shifting in the activation site, it should be positioned far from the CA1 area where synaptic transmission occurs.</p>
<p>&#x03BC;MS was also used to stimulate apical dendrites of the cortex&#x2019;s pyramidal neurons (PNs) without activating horizontally oriented passing axons in mice (<xref ref-type="bibr" rid="B28">Lee et al., 2016</xref>). In these experiments, the length of a single dendrite is usually several hundred micrometers. This is a dimension that is unlikely to be missed by the &#x03BC;MS technology.</p>
</sec>
</sec>
<sec id="S5" sec-type="conclusion">
<title>5. Conclusion and future work</title>
<p>This paper demonstrates that figure-eight &#x03BC;MS could generate consistent axonal activation with high spatial resolution. Because the site of activation was unaffected by the changes in the direction and intensity of the coil current, and was likewise unaffected by the increase in coil-axon distance, neural stimulation using the figure-eight micro-coil could have a higher tolerance for human errors in clinical applications. The coil needs less intensive coil current for axonal activation, which could be a tremendous advantage in avoiding unwanted thermal effects associated with the &#x03BC;MS technology. Although simulation results coming from this work must be fully validated with future experiments, the obvious advantages provided by the figure-eight micro-coil advocate for further efforts in the design, manufacturing, and development of the figure-eight coil in &#x03BC;MS technology.</p>
<p>Previous studies have suggested that axonal bending (<xref ref-type="bibr" rid="B30">Maccabee et al., 1993</xref>) and the presence of dendritic/axonal terminals (<xref ref-type="bibr" rid="B38">Radman et al., 2009</xref>; <xref ref-type="bibr" rid="B39">Rahman et al., 2013</xref>) could be the target of electric field activation, which are not integrated in the current model. It would be interesting and clinically relevant to see if the figure-eight coil could provide focal and consistent neural activation when these morphological details are represented in the model. Finite element modeling would be a fitting tool to represent such complicated neural morphology. The excellent agreement between activating function analysis and the site of neural activation (provided by NEURON modeling) allows us to speculate that, if such complicated neural geometry should be considered, activating function analysis could still serve as a powerful tool to predict the site of neural activation.</p>
<p>Increased resolution of stimulation with the figure-eight micro-coil could have a significant impact on basic research. For example, microelectric stimulation on a single neuron has been used to alter the activity of large networks and animal behavior (<xref ref-type="bibr" rid="B16">Houweling and Brecht, 2008</xref>). Optogenetically driving a fraction of neurons in the mouse primary somatosensory cortex could control behavior in freely moving mice (<xref ref-type="bibr" rid="B17">Huber et al., 2008</xref>). These studies establish causal links between the activity of single neurons and perceptual and cognitive functions, as well as animal behavior. It is expected that future &#x03BC;MS technology could achieve single cell resolution, therefore providing a more biocompatible tool to selectively activate or inhibit neurons, and could provide causal evidence of the functional roles of specific neurons. Future work will test the specificity and consistency of single neuron stimulation by the figure-eight micro-coil configuration using <italic>in vitro</italic> or <italic>in vivo</italic> animal models.</p>
</sec>
<sec id="S6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in this study are included in the article/<xref ref-type="supplementary-material" rid="FS1">Supplementary material</xref>, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="S7" sec-type="author-contributions">
<title>Author contributions</title>
<p>HY and VH: biophysics modeling. HY: NEURON simulation. All authors manuscript drafting and read and approved the final manuscript.</p>
</sec>
</body>
<back>
<sec id="S8" sec-type="funding-information">
<title>Funding</title>
<p>This work was supported by the National Institute of Neurological Disorders and Stroke of the National Institutes of Health (R03NS130387) and a research support grant (1100) from Loyola University Chicago. VH was supported by a Biology Summer Fellowship from Loyola University Chicago.</p>
</sec>
<ack>
<p>Dr. Vincent Chen provided valuable suggestions in the development of the biophysics model of the coil.</p>
</ack>
<sec id="S10" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="S11" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="S12" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fncom.2023.1105505/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fncom.2023.1105505/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Image_1.tif" id="FS1" mimetype="image/tiff" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Figure 1</label>
<caption><p>Figure-eight micro-coil used in neural stimulation experiments. The coil was assembled from two commercially available inductors. <bold>(A)</bold> Coil internal structure revealed by the chemical removal of the encapsulation. <bold>(B)</bold> Assembly of the figure-eight coil with the two inductors stacked on top of each other, with opposite wire winding.</p></caption>
</supplementary-material>
<supplementary-material xlink:href="Image_2.tif" id="FS2" mimetype="image/tiff" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Figure 2</label>
<caption><p>Extracellular potential along the modeled axon in NEURON simulation. The potential was generated by a positive pulse delivered through a circular micro-coil <bold>(A)</bold> versus a figure-eight micro-coil <bold>(B)</bold>.</p></caption>
</supplementary-material>
</sec>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ahmed</surname> <given-names>S. I.</given-names></name> <name><surname>Javed</surname> <given-names>G.</given-names></name> <name><surname>Mubeen</surname> <given-names>B.</given-names></name> <name><surname>Bareeqa</surname> <given-names>S. B.</given-names></name> <name><surname>Rasheed</surname> <given-names>H.</given-names></name> <name><surname>Rehman</surname> <given-names>A.</given-names></name><etal/></person-group> (<year>2018</year>). <article-title>Robotics in neurosurgery: A literature review.</article-title> <source><italic>J. Pak Med. Assoc.</italic></source> <volume>68</volume> <fpage>258</fpage>&#x2013;<lpage>263</lpage>.</citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Alekseichuk</surname> <given-names>I.</given-names></name> <name><surname>Falchier</surname> <given-names>A. Y.</given-names></name> <name><surname>Linn</surname> <given-names>G.</given-names></name> <name><surname>Xu</surname> <given-names>T.</given-names></name> <name><surname>Milham</surname> <given-names>M. P.</given-names></name> <name><surname>Schroeder</surname> <given-names>C. E.</given-names></name><etal/></person-group> (<year>2019</year>). <article-title>Electric field dynamics in the brain during multi-electrode transcranial electric stimulation.</article-title> <source><italic>Nat. Commun.</italic></source> <volume>10</volume>:<fpage>2573</fpage>. <pub-id pub-id-type="doi">10.1038/s41467-019-10581-7</pub-id> <pub-id pub-id-type="pmid">31189931</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Barker</surname> <given-names>A. T.</given-names></name> <name><surname>Jalinous</surname> <given-names>R.</given-names></name> <name><surname>Freeston</surname> <given-names>I. L.</given-names></name></person-group> (<year>1985</year>). <article-title>Non-invasive magnetic stimulation of human motor cortex.</article-title> <source><italic>Lancet</italic></source> <volume>1</volume> <fpage>1106</fpage>&#x2013;<lpage>1107</lpage>. <pub-id pub-id-type="doi">10.1016/s0140-6736(85)92413-4</pub-id> <pub-id pub-id-type="pmid">34813082</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bohning</surname> <given-names>D. E.</given-names></name> <name><surname>Shastri</surname> <given-names>A.</given-names></name> <name><surname>McConnell</surname> <given-names>K. A.</given-names></name> <name><surname>Nahas</surname> <given-names>Z.</given-names></name> <name><surname>Lorberbaum</surname> <given-names>J. P.</given-names></name> <name><surname>Roberts</surname> <given-names>D. R.</given-names></name><etal/></person-group> (<year>1999</year>). <article-title>A combined TMS/fMRI study of intensity-dependent TMS over motor cortex.</article-title> <source><italic>Biol. Psychiatry</italic></source> <volume>45</volume> <fpage>385</fpage>&#x2013;<lpage>394</lpage>. <pub-id pub-id-type="doi">10.1016/s0006-3223(98)00368-0</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bonmassar</surname> <given-names>G.</given-names></name> <name><surname>Lee</surname> <given-names>S. W.</given-names></name> <name><surname>Freeman</surname> <given-names>D. K.</given-names></name> <name><surname>Polasek</surname> <given-names>M.</given-names></name> <name><surname>Fried</surname> <given-names>S. I.</given-names></name> <name><surname>Gale</surname> <given-names>J. T.</given-names></name></person-group> (<year>2012</year>). <article-title>Microscopic magnetic stimulation of neural tissue.</article-title> <source><italic>Nat. Commun.</italic></source> <volume>3</volume>:<fpage>921</fpage>. <pub-id pub-id-type="doi">10.1038/ncomms1914</pub-id> <pub-id pub-id-type="pmid">22735449</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cattell</surname> <given-names>M.</given-names></name> <name><surname>Gerard</surname> <given-names>R. W.</given-names></name></person-group> (<year>1935</year>). <article-title>The &#x201C;inhibitory&#x201D; effect of high-frequency stimulation and the excitation state of nerve.</article-title> <source><italic>J. Physiol.</italic></source> <volume>83</volume> <fpage>407</fpage>&#x2013;<lpage>415</lpage>. <pub-id pub-id-type="doi">10.1113/jphysiol.1935.sp003238</pub-id> <pub-id pub-id-type="pmid">16994640</pub-id></citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cogan</surname> <given-names>S. F.</given-names></name></person-group> (<year>2008</year>). <article-title>Neural stimulation and recording electrodes.</article-title> <source><italic>Annu. Rev. Biomed. Eng.</italic></source> <volume>10</volume> <fpage>275</fpage>&#x2013;<lpage>309</lpage>. <pub-id pub-id-type="doi">10.1146/annurev.bioeng.10.061807.160518</pub-id> <pub-id pub-id-type="pmid">18429704</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Colella</surname> <given-names>M.</given-names></name> <name><surname>Liberti</surname> <given-names>M.</given-names></name> <name><surname>Apollonio</surname> <given-names>F.</given-names></name> <name><surname>Bonmassar</surname> <given-names>G.</given-names></name></person-group> (<year>2021</year>). &#x201C;<article-title>A miniaturized ultra-focal magnetic stimulator and its preliminary application to the peripheral nervous system</article-title>,&#x201D; in <source><italic>Brain and human body modeling 2020: Computational human models presented at EMBC 2019 and the BRAIN initiative(R) 2019 meeting</italic></source>, <role>eds</role> <person-group person-group-type="editor"><name><surname>Makarov</surname> <given-names>S. N.</given-names></name> <name><surname>Noetscher</surname> <given-names>G. M.</given-names></name> <name><surname>Nummenmaa</surname> <given-names>A.</given-names></name></person-group> (<publisher-loc>Cham</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>167</fpage>&#x2013;<lpage>176</lpage>. <pub-id pub-id-type="doi">10.1007/978-3-030-45623-8_9</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fernandes de Oliveira Santos</surname> <given-names>B.</given-names></name> <name><surname>de Araujo Paz</surname> <given-names>D.</given-names></name> <name><surname>Fernandes</surname> <given-names>V. M.</given-names></name> <name><surname>Dos Santos</surname> <given-names>J. C.</given-names></name> <name><surname>Chaddad-Neto</surname> <given-names>F. E. A.</given-names></name> <name><surname>Sousa</surname> <given-names>A. C. S.</given-names></name><etal/></person-group> (<year>2021</year>). <article-title>Minimally invasive supratentorial neurosurgical approaches guided by Smartphone app and compass.</article-title> <source><italic>Sci. Rep.</italic></source> <volume>11</volume>:<fpage>6778</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-021-85472-3</pub-id> <pub-id pub-id-type="pmid">33762597</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gluckman</surname> <given-names>B. J.</given-names></name> <name><surname>Neel</surname> <given-names>E. J.</given-names></name> <name><surname>Netoff</surname> <given-names>T. I.</given-names></name> <name><surname>Ditto</surname> <given-names>W. L.</given-names></name> <name><surname>Spano</surname> <given-names>M. L.</given-names></name> <name><surname>Schiff</surname> <given-names>S. J.</given-names></name></person-group> (<year>1996</year>). <article-title>Electric field suppression of epileptiform activity in hippocampal slices.</article-title> <source><italic>J. Neurophysiol.</italic></source> <volume>76</volume> <fpage>4202</fpage>&#x2013;<lpage>4205</lpage>. <pub-id pub-id-type="doi">10.1152/jn.1996.76.6.4202</pub-id> <pub-id pub-id-type="pmid">8985916</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Golestanirad</surname> <given-names>L.</given-names></name> <name><surname>Gale</surname> <given-names>J. T.</given-names></name> <name><surname>Manzoor</surname> <given-names>N. F.</given-names></name> <name><surname>Park</surname> <given-names>H. J.</given-names></name> <name><surname>Glait</surname> <given-names>L.</given-names></name> <name><surname>Haer</surname> <given-names>F.</given-names></name><etal/></person-group> (<year>2018</year>). <article-title>Solenoidal micromagnetic stimulation enables activation of axons with specific orientation.</article-title> <source><italic>Front. Physiol.</italic></source> <volume>9</volume>:<fpage>724</fpage>. <pub-id pub-id-type="doi">10.3389/fphys.2018.00724</pub-id> <pub-id pub-id-type="pmid">30140230</pub-id></citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gunter</surname> <given-names>C.</given-names></name> <name><surname>Delbeke</surname> <given-names>J.</given-names></name> <name><surname>Ortiz-Catalan</surname> <given-names>M.</given-names></name></person-group> (<year>2019</year>). <article-title>Safety of long-term electrical peripheral nerve stimulation: review of the state of the art.</article-title> <source><italic>J. Neuroeng. Rehabil.</italic></source> <volume>16</volume>:<fpage>13</fpage>. <pub-id pub-id-type="doi">10.1186/s12984-018-0474-8</pub-id> <pub-id pub-id-type="pmid">30658656</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Habibollahi Saatlou</surname> <given-names>F.</given-names></name> <name><surname>Rogasch</surname> <given-names>N. C.</given-names></name> <name><surname>McNair</surname> <given-names>N. A.</given-names></name> <name><surname>Biabani</surname> <given-names>M.</given-names></name> <name><surname>Pillen</surname> <given-names>S. D.</given-names></name> <name><surname>Marshall</surname> <given-names>T. R.</given-names></name><etal/></person-group> (<year>2018</year>). <article-title>MAGIC: An open-source MATLAB toolbox for external control of transcranial magnetic stimulation devices.</article-title> <source><italic>Brain Stimul.</italic></source> <volume>11</volume> <fpage>1189</fpage>&#x2013;<lpage>1191</lpage>. <pub-id pub-id-type="doi">10.1016/j.brs.2018.05.015</pub-id> <pub-id pub-id-type="pmid">29885859</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hines</surname> <given-names>M. L.</given-names></name> <name><surname>Carnevale</surname> <given-names>N. T.</given-names></name></person-group> (<year>1997</year>). <article-title>The NEURON simulation environment.</article-title> <source><italic>Neural Comput.</italic></source> <volume>9</volume> <fpage>1179</fpage>&#x2013;<lpage>1209</lpage>. <pub-id pub-id-type="doi">10.1162/neco.1997.9.6.1179</pub-id> <pub-id pub-id-type="pmid">9248061</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hodgkin</surname> <given-names>A. L.</given-names></name> <name><surname>Huxley</surname> <given-names>A. F.</given-names></name></person-group> (<year>1952</year>). <article-title>A quantitative description of membrane current and its application to conduction and excitation in nerve.</article-title> <source><italic>J. Physiol.</italic></source> <volume>117</volume> <fpage>500</fpage>&#x2013;<lpage>544</lpage>. <pub-id pub-id-type="doi">10.1113/jphysiol.1952.sp004764</pub-id> <pub-id pub-id-type="pmid">12991237</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Houweling</surname> <given-names>A. R.</given-names></name> <name><surname>Brecht</surname> <given-names>M.</given-names></name></person-group> (<year>2008</year>). <article-title>Behavioural report of single neuron stimulation in somatosensory cortex.</article-title> <source><italic>Nature</italic></source> <volume>451</volume> <fpage>65</fpage>&#x2013;<lpage>68</lpage>. <pub-id pub-id-type="doi">10.1038/nature06447</pub-id> <pub-id pub-id-type="pmid">18094684</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huber</surname> <given-names>D.</given-names></name> <name><surname>Petreanu</surname> <given-names>L.</given-names></name> <name><surname>Ghitani</surname> <given-names>N.</given-names></name> <name><surname>Ranade</surname> <given-names>S.</given-names></name> <name><surname>Hromadka</surname> <given-names>T.</given-names></name> <name><surname>Mainen</surname> <given-names>Z.</given-names></name><etal/></person-group> (<year>2008</year>). <article-title>Sparse optical microstimulation in barrel cortex drives learned behaviour in freely moving mice.</article-title> <source><italic>Nature</italic></source> <volume>451</volume> <fpage>61</fpage>&#x2013;<lpage>64</lpage>. <pub-id pub-id-type="doi">10.1038/nature06445</pub-id> <pub-id pub-id-type="pmid">18094685</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jefferys</surname> <given-names>J. G.</given-names></name></person-group> (<year>1981</year>). <article-title>Influence of electric fields on the excitability of granule cells in guinea-pig hippocampal slices.</article-title> <source><italic>J. Physiol.</italic></source> <volume>319</volume> <fpage>143</fpage>&#x2013;<lpage>152</lpage>. <pub-id pub-id-type="doi">10.1113/jphysiol.1981.sp013897</pub-id> <pub-id pub-id-type="pmid">7320909</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jeong</surname> <given-names>H.</given-names></name> <name><surname>Cho</surname> <given-names>A.</given-names></name> <name><surname>Ay</surname> <given-names>I.</given-names></name> <name><surname>Bonmassar</surname> <given-names>G.</given-names></name></person-group> (<year>2022</year>). <article-title>Short-pulsed micro-magnetic stimulation of the vagus nerve.</article-title> <source><italic>Front. Physiol.</italic></source> <volume>13</volume>:<fpage>938101</fpage>. <pub-id pub-id-type="doi">10.3389/fphys.2022.938101</pub-id> <pub-id pub-id-type="pmid">36277182</pub-id></citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jeong</surname> <given-names>H.</given-names></name> <name><surname>Deng</surname> <given-names>J.</given-names></name> <name><surname>Bonmassar</surname> <given-names>G.</given-names></name></person-group> (<year>2021</year>). <article-title>Planar figure-8 coils for ultra-focal and directional micromagnetic brain stimulation.</article-title> <source><italic>J. Vac. Sci. Technol. B Nanotechnol. Microelectron.</italic></source> <volume>39</volume>:<fpage>063202</fpage>. <pub-id pub-id-type="doi">10.1116/6.0001281</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Joucla</surname> <given-names>S.</given-names></name> <name><surname>Gliere</surname> <given-names>A.</given-names></name> <name><surname>Yvert</surname> <given-names>B.</given-names></name></person-group> (<year>2014</year>). <article-title>Current approaches to model extracellular electrical neural microstimulation.</article-title> <source><italic>Front. Comput. Neurosci.</italic></source> <volume>8</volume>:<fpage>13</fpage>. <pub-id pub-id-type="doi">10.3389/fncom.2014.00013</pub-id> <pub-id pub-id-type="pmid">24600381</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Khalifa</surname> <given-names>A.</given-names></name> <name><surname>Zaeimbashi</surname> <given-names>M.</given-names></name> <name><surname>Zhou</surname> <given-names>T. X.</given-names></name> <name><surname>Abrishami</surname> <given-names>S. M.</given-names></name> <name><surname>Sun</surname> <given-names>N.</given-names></name> <name><surname>Park</surname> <given-names>S.</given-names></name><etal/></person-group> (<year>2021</year>). <article-title>The development of microfabricated solenoids with magnetic cores for micromagnetic neural stimulation.</article-title> <source><italic>Microsyst. Nanoeng.</italic></source> <volume>7</volume>:<fpage>91</fpage>. <pub-id pub-id-type="doi">10.1038/s41378-021-00320-8</pub-id> <pub-id pub-id-type="pmid">34786205</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname> <given-names>U. S.</given-names></name> <name><surname>Mahroo</surname> <given-names>O. A.</given-names></name> <name><surname>Mollon</surname> <given-names>J. D.</given-names></name> <name><surname>Yu-Wai-Man</surname> <given-names>P.</given-names></name></person-group> (<year>2021</year>). <article-title>Retinal ganglion cells-diversity of cell types and clinical relevance.</article-title> <source><italic>Front. Neurol.</italic></source> <volume>12</volume>:<fpage>661938</fpage>. <pub-id pub-id-type="doi">10.3389/fneur.2021.661938</pub-id> <pub-id pub-id-type="pmid">34093409</pub-id></citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Koivuniemi</surname> <given-names>A.</given-names></name> <name><surname>Wilks</surname> <given-names>S. J.</given-names></name> <name><surname>Woolley</surname> <given-names>A. J.</given-names></name> <name><surname>Otto</surname> <given-names>K. J.</given-names></name></person-group> (<year>2011</year>). <article-title>Multimodal, longitudinal assessment of intracortical microstimulation.</article-title> <source><italic>Prog. Brain Res.</italic></source> <volume>194</volume> <fpage>131</fpage>&#x2013;<lpage>144</lpage>. <pub-id pub-id-type="doi">10.1016/B978-0-444-53815-4.00011-X</pub-id> <pub-id pub-id-type="pmid">21867800</pub-id></citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>J. I.</given-names></name> <name><surname>Seist</surname> <given-names>R.</given-names></name> <name><surname>McInturff</surname> <given-names>S.</given-names></name> <name><surname>Lee</surname> <given-names>D. J.</given-names></name> <name><surname>Brown</surname> <given-names>M. C.</given-names></name> <name><surname>Stankovic</surname> <given-names>K. M.</given-names></name><etal/></person-group> (<year>2022</year>). <article-title>Magnetic stimulation allows focal activation of the mouse cochlea.</article-title> <source><italic>Elife</italic></source> <volume>11</volume>:<fpage>e76682</fpage>. <pub-id pub-id-type="doi">10.7554/eLife.76682</pub-id> <pub-id pub-id-type="pmid">35608242</pub-id></citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>S. W.</given-names></name> <name><surname>Fried</surname> <given-names>S. I.</given-names></name></person-group> (<year>2014</year>). <article-title>The response of L5 pyramidal neurons of the PFC to magnetic stimulation from a micro-coil.</article-title> <source><italic>Conf. Proc. IEEE Eng. Med. Biol. Soc.</italic></source> <volume>2014</volume> <fpage>6125</fpage>&#x2013;<lpage>6128</lpage>. <pub-id pub-id-type="doi">10.1109/EMBC.2014.6945027</pub-id> <pub-id pub-id-type="pmid">25571395</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>S. W.</given-names></name> <name><surname>Fried</surname> <given-names>S. I.</given-names></name></person-group> (<year>2017</year>). <article-title>Enhanced control of cortical pyramidal neurons with micromagnetic stimulation.</article-title> <source><italic>IEEE Trans. Neural Syst. Rehabil. Eng.</italic></source> <volume>25</volume> <fpage>1375</fpage>&#x2013;<lpage>1386</lpage>. <pub-id pub-id-type="doi">10.1109/TNSRE.2016.2631446</pub-id> <pub-id pub-id-type="pmid">27893396</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>S. W.</given-names></name> <name><surname>Fallegger</surname> <given-names>F.</given-names></name> <name><surname>Casse</surname> <given-names>B. D.</given-names></name> <name><surname>Fried</surname> <given-names>S. I.</given-names></name></person-group> (<year>2016</year>). <article-title>Implantable microcoils for intracortical magnetic stimulation.</article-title> <source><italic>Sci. Adv.</italic></source> <volume>2</volume>:<fpage>e1600889</fpage>. <pub-id pub-id-type="doi">10.1126/sciadv.1600889</pub-id> <pub-id pub-id-type="pmid">27957537</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>S. W.</given-names></name> <name><surname>Thyagarajan</surname> <given-names>K.</given-names></name> <name><surname>Fried</surname> <given-names>S. I.</given-names></name></person-group> (<year>2019</year>). <article-title>Micro-coil design influences the spatial extent of responses to intracortical magnetic stimulation.</article-title> <source><italic>IEEE Trans. Biomed. Eng.</italic></source> <volume>66</volume> <fpage>1680</fpage>&#x2013;<lpage>1694</lpage>. <pub-id pub-id-type="doi">10.1109/Tbme.2018.2877713</pub-id> <pub-id pub-id-type="pmid">30369434</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Maccabee</surname> <given-names>P. J.</given-names></name> <name><surname>Amassian</surname> <given-names>V. E.</given-names></name> <name><surname>Eberle</surname> <given-names>L. P.</given-names></name> <name><surname>Cracco</surname> <given-names>R. Q.</given-names></name></person-group> (<year>1993</year>). <article-title>Magnetic coil stimulation of straight and bent amphibian and mammalian peripheral nerve in vitro: locus of excitation.</article-title> <source><italic>J. Physiol.</italic></source> <volume>460</volume> <fpage>201</fpage>&#x2013;<lpage>219</lpage>. <pub-id pub-id-type="doi">10.1113/jphysiol.1993.sp019467</pub-id> <pub-id pub-id-type="pmid">8487192</pub-id></citation></ref>
<ref id="B31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Minusa</surname> <given-names>S.</given-names></name> <name><surname>Muramatsu</surname> <given-names>S.</given-names></name> <name><surname>Osanai</surname> <given-names>H.</given-names></name> <name><surname>Tateno</surname> <given-names>T.</given-names></name></person-group> (<year>2019</year>). <article-title>A multichannel magnetic stimulation system using submillimeter-sized coils: system development and experimental application to rodent brain in vivo.</article-title> <source><italic>J. Neural. Eng.</italic></source> <volume>16</volume>:<fpage>066014</fpage>. <pub-id pub-id-type="doi">10.1088/1741-2552/ab3187</pub-id> <pub-id pub-id-type="pmid">31642445</pub-id></citation></ref>
<ref id="B32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Minusa</surname> <given-names>S.</given-names></name> <name><surname>Osanai</surname> <given-names>H.</given-names></name> <name><surname>Tateno</surname> <given-names>T.</given-names></name></person-group> (<year>2018</year>). <article-title>Micromagnetic stimulation of the mouse auditory cortex in vivo using an implantable solenoid system.</article-title> <source><italic>IEEE Trans. Biomed. Eng.</italic></source> <volume>65</volume> <fpage>1301</fpage>&#x2013;<lpage>1310</lpage>. <pub-id pub-id-type="doi">10.1109/TBME.2017.2748136</pub-id> <pub-id pub-id-type="pmid">28880154</pub-id></citation></ref>
<ref id="B33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Osanai</surname> <given-names>H.</given-names></name> <name><surname>Minusa</surname> <given-names>S.</given-names></name> <name><surname>Tateno</surname> <given-names>T.</given-names></name></person-group> (<year>2018</year>). <article-title>Micro-coil-induced Inhomogeneous electric field produces sound-driven-like neural responses in microcircuits of the mouse auditory cortex in vivo.</article-title> <source><italic>Neuroscience</italic></source> <volume>371</volume> <fpage>346</fpage>&#x2013;<lpage>370</lpage>. <pub-id pub-id-type="doi">10.1016/j.neuroscience.2017.12.008</pub-id> <pub-id pub-id-type="pmid">29246784</pub-id></citation></ref>
<ref id="B34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Park</surname> <given-names>H. J.</given-names></name> <name><surname>Bonmassar</surname> <given-names>G.</given-names></name> <name><surname>Kaltenbach</surname> <given-names>J. A.</given-names></name> <name><surname>Machado</surname> <given-names>A. G.</given-names></name> <name><surname>Manzoor</surname> <given-names>N. F.</given-names></name> <name><surname>Gale</surname> <given-names>J. T.</given-names></name></person-group> (<year>2013</year>). <article-title>Activation of the central nervous system induced by micro-magnetic stimulation.</article-title> <source><italic>Nat. Commun.</italic></source> <volume>4</volume>:<fpage>2463</fpage>. <pub-id pub-id-type="doi">10.1038/Ncomms3463</pub-id> <pub-id pub-id-type="pmid">24030203</pub-id></citation></ref>
<ref id="B35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Park</surname> <given-names>Y. S.</given-names></name> <name><surname>Kim</surname> <given-names>J.</given-names></name> <name><surname>Chang</surname> <given-names>W. S.</given-names></name> <name><surname>Chang</surname> <given-names>J. W.</given-names></name></person-group> (<year>2011</year>). <article-title>Management of a DBS system in patients with traumatic brain injury: Case report.</article-title> <source><italic>Neuromodulation</italic></source> <volume>14</volume> <fpage>214</fpage>&#x2013;<lpage>218; discussion218</lpage>. <pub-id pub-id-type="doi">10.1111/j.1525-1403.2011.00348.x</pub-id> <pub-id pub-id-type="pmid">21992242</pub-id></citation></ref>
<ref id="B36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pashut</surname> <given-names>T.</given-names></name> <name><surname>Magidov</surname> <given-names>D.</given-names></name> <name><surname>Ben-Porat</surname> <given-names>H.</given-names></name> <name><surname>Wolfus</surname> <given-names>S.</given-names></name> <name><surname>Friedman</surname> <given-names>A.</given-names></name> <name><surname>Perel</surname> <given-names>E.</given-names></name><etal/></person-group> (<year>2014</year>). <article-title>Patch-clamp recordings of rat neurons from acute brain slices of the somatosensory cortex during magnetic stimulation.</article-title> <source><italic>Front. Cell Neurosci.</italic></source> <volume>8</volume>:<fpage>145</fpage>. <pub-id pub-id-type="doi">10.3389/fncel.2014.00145</pub-id> <pub-id pub-id-type="pmid">24917788</pub-id></citation></ref>
<ref id="B37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Polikov</surname> <given-names>V. S.</given-names></name> <name><surname>Tresco</surname> <given-names>P. A.</given-names></name> <name><surname>Reichert</surname> <given-names>W. M.</given-names></name></person-group> (<year>2005</year>). <article-title>Response of brain tissue to chronically implanted neural electrodes.</article-title> <source><italic>J. Neurosci. Methods</italic></source> <volume>148</volume> <fpage>1</fpage>&#x2013;<lpage>18</lpage>. <pub-id pub-id-type="doi">10.1016/j.jneumeth.2005.08.015</pub-id> <pub-id pub-id-type="pmid">16198003</pub-id></citation></ref>
<ref id="B38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Radman</surname> <given-names>T.</given-names></name> <name><surname>Ramos</surname> <given-names>R. L.</given-names></name> <name><surname>Brumberg</surname> <given-names>J. C.</given-names></name> <name><surname>Bikson</surname> <given-names>M.</given-names></name></person-group> (<year>2009</year>). <article-title>Role of cortical cell type and morphology in subthreshold and suprathreshold uniform electric field stimulation in vitro.</article-title> <source><italic>Brain Stimul.</italic></source> <volume>2</volume> <fpage>e211</fpage>&#x2013;<lpage>e213</lpage>. <pub-id pub-id-type="doi">10.1016/j.brs.2009.03.007</pub-id> <pub-id pub-id-type="pmid">20161507</pub-id></citation></ref>
<ref id="B39"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rahman</surname> <given-names>A.</given-names></name> <name><surname>Reato</surname> <given-names>D.</given-names></name> <name><surname>Arlotti</surname> <given-names>M.</given-names></name> <name><surname>Gasca</surname> <given-names>F.</given-names></name> <name><surname>Datta</surname> <given-names>A.</given-names></name> <name><surname>Parra</surname> <given-names>L. C.</given-names></name><etal/></person-group> (<year>2013</year>). <article-title>Cellular effects of acute direct current stimulation: Somatic and synaptic terminal effects.</article-title> <source><italic>J. Physiol.</italic></source> <volume>591</volume> <fpage>2563</fpage>&#x2013;<lpage>2578</lpage>. <pub-id pub-id-type="doi">10.1113/jphysiol.2012.247171</pub-id> <pub-id pub-id-type="pmid">23478132</pub-id></citation></ref>
<ref id="B40"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rastogi</surname> <given-names>P.</given-names></name> <name><surname>Lee</surname> <given-names>E. G.</given-names></name> <name><surname>Hadimani</surname> <given-names>R. L.</given-names></name> <name><surname>Jiles</surname> <given-names>D. C.</given-names></name></person-group> (<year>2017</year>). <article-title>Transcranial magnetic stimulation-coil design with improved focality.</article-title> <source><italic>AIP Adv.</italic></source> <volume>7</volume>:<fpage>056705</fpage>. <pub-id pub-id-type="doi">10.1063/1.4973604</pub-id></citation></ref>
<ref id="B41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rattay</surname> <given-names>F.</given-names></name></person-group> (<year>1986</year>). <article-title>Analysis of models for external stimulation of axons.</article-title> <source><italic>IEEE Trans. Biomed. Eng.</italic></source> <volume>33</volume> <fpage>974</fpage>&#x2013;<lpage>977</lpage>. <pub-id pub-id-type="doi">10.1109/TBME.1986.325670</pub-id> <pub-id pub-id-type="pmid">3770787</pub-id></citation></ref>
<ref id="B42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rattay</surname> <given-names>F.</given-names></name></person-group> (<year>1989</year>). <article-title>Analysis of models for extracellular fiber stimulation.</article-title> <source><italic>IEEE Trans. Biomed. Eng.</italic></source> <volume>36</volume> <fpage>676</fpage>&#x2013;<lpage>682</lpage>. <pub-id pub-id-type="doi">10.1109/10.32099</pub-id></citation></ref>
<ref id="B43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rattay</surname> <given-names>F.</given-names></name> <name><surname>Resatz</surname> <given-names>S.</given-names></name></person-group> (<year>2007</year>). <article-title>Dipole distance for minimum threshold current to stimulate unmyelinated axons with microelectrodes.</article-title> <source><italic>IEEE Trans. Biomed. Eng.</italic></source> <volume>54</volume> <fpage>158</fpage>&#x2013;<lpage>162</lpage>. <pub-id pub-id-type="doi">10.1109/TBME.2006.883730</pub-id> <pub-id pub-id-type="pmid">17260868</pub-id></citation></ref>
<ref id="B44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ravazzani</surname> <given-names>P.</given-names></name> <name><surname>Ruohonen</surname> <given-names>J.</given-names></name> <name><surname>Grandori</surname> <given-names>F.</given-names></name> <name><surname>Tognola</surname> <given-names>G.</given-names></name></person-group> (<year>1996</year>). <article-title>Magnetic stimulation of the nervous system: induced electric field in unbounded, semi-infinite, spherical, and cylindrical media.</article-title> <source><italic>Ann. Biomed. Eng.</italic></source> <volume>24</volume> <fpage>606</fpage>&#x2013;<lpage>616</lpage>. <pub-id pub-id-type="doi">10.1007/BF02684229</pub-id> <pub-id pub-id-type="pmid">8886241</pub-id></citation></ref>
<ref id="B45"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Saha</surname> <given-names>R.</given-names></name> <name><surname>Faramarzi</surname> <given-names>S.</given-names></name> <name><surname>Bloom</surname> <given-names>R. P.</given-names></name> <name><surname>Benally</surname> <given-names>O. J.</given-names></name> <name><surname>Wu</surname> <given-names>K.</given-names></name> <name><surname>di Girolamo</surname> <given-names>A.</given-names></name><etal/></person-group> (<year>2022a</year>). <article-title>Strength-frequency curve for micromagnetic neurostimulation through excitatory postsynaptic potentials (EPSPs) on rat hippocampal neurons and numerical modeling of magnetic microcoil (mu coil).</article-title> <source><italic>J. Neural Eng.</italic></source> <volume>19</volume>. <pub-id pub-id-type="doi">10.1088/1741-2552/ac4baf</pub-id> <pub-id pub-id-type="pmid">35030549</pub-id></citation></ref>
<ref id="B46"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Saha</surname> <given-names>R.</given-names></name> <name><surname>Wu</surname> <given-names>K.</given-names></name> <name><surname>Bloom</surname> <given-names>R.</given-names></name> <name><surname>Liang</surname> <given-names>S.</given-names></name> <name><surname>Tonini</surname> <given-names>D.</given-names></name> <name><surname>Wang</surname> <given-names>J. P.</given-names></name></person-group> (<year>2022b</year>). <article-title>A review on magnetic and spintronic neurostimulation: challenges and prospects.</article-title> <source><italic>Nanotechnology</italic></source> <volume>33</volume>. <pub-id pub-id-type="doi">10.1088/1361-6528/ac49be</pub-id> <pub-id pub-id-type="pmid">35013010</pub-id></citation></ref>
<ref id="B47"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Saito</surname> <given-names>A.</given-names></name></person-group> (<year>2021</year>). <article-title>Fabrication of a miniature figure-of-eight coil for micromagnetic stimulation on neuronal tissue.</article-title> <source><italic>URSI Radio Sci. Lett.</italic></source> <volume>3</volume>. <pub-id pub-id-type="doi">10.46620/21-0008</pub-id></citation></ref>
<ref id="B48"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sankarasubramanian</surname> <given-names>V.</given-names></name> <name><surname>Buitenweg</surname> <given-names>J. R.</given-names></name> <name><surname>Holsheimer</surname> <given-names>J.</given-names></name> <name><surname>Veltink</surname> <given-names>P.</given-names></name></person-group> (<year>2011</year>). <article-title>Electrode alignment of transverse tripoles using a percutaneous triple-lead approach in spinal cord stimulation.</article-title> <source><italic>J. Neural. Eng.</italic></source> <volume>8</volume>:<fpage>016010</fpage>. <pub-id pub-id-type="doi">10.1088/1741-2560/8/1/016010</pub-id></citation></ref>
<ref id="B49"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sekino</surname> <given-names>M.</given-names></name> <name><surname>Ohsaki</surname> <given-names>H.</given-names></name> <name><surname>Takiyama</surname> <given-names>Y.</given-names></name> <name><surname>Yamamoto</surname> <given-names>K.</given-names></name> <name><surname>Matsuzaki</surname> <given-names>T.</given-names></name> <name><surname>Yasumuro</surname> <given-names>Y.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>Eccentric figure-eight coils for transcranial magnetic stimulation.</article-title> <source><italic>Bioelectromagnetics</italic></source> <volume>36</volume> <fpage>55</fpage>&#x2013;<lpage>65</lpage>. <pub-id pub-id-type="doi">10.1002/bem.21886</pub-id> <pub-id pub-id-type="pmid">25399864</pub-id></citation></ref>
<ref id="B50"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tai</surname> <given-names>C.</given-names></name> <name><surname>de Groat</surname> <given-names>W. C.</given-names></name> <name><surname>Roppolo</surname> <given-names>J. R.</given-names></name></person-group> (<year>2005</year>). <article-title>Simulation analysis of conduction block in unmyelinated axons induced by high-frequency biphasic electrical currents.</article-title> <source><italic>IEEE Trans. Biomed. Eng.</italic></source> <volume>52</volume> <fpage>1323</fpage>&#x2013;<lpage>1332</lpage>. <pub-id pub-id-type="doi">10.1109/tbme.2005.847561</pub-id> <pub-id pub-id-type="pmid">16041996</pub-id></citation></ref>
<ref id="B51"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Talebinejad</surname> <given-names>M.</given-names></name> <name><surname>Musallam</surname> <given-names>S.</given-names></name></person-group> (<year>2010</year>). <article-title>Effects of TMS coil geometry on stimulation specificity.</article-title> <source><italic>Annu. Int. Conf. IEEE Eng. Med. Biol. Soc.</italic></source> <volume>2010</volume> <fpage>1507</fpage>&#x2013;<lpage>1510</lpage>. <pub-id pub-id-type="doi">10.1109/IEMBS.2010.5626840</pub-id> <pub-id pub-id-type="pmid">21096368</pub-id></citation></ref>
<ref id="B52"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tanner</surname> <given-names>J. A.</given-names></name></person-group> (<year>1962</year>). <article-title>Reversible blocking of nerve conduction by alternating-current excitation.</article-title> <source><italic>Nature</italic></source> <volume>195</volume> <fpage>712</fpage>&#x2013;<lpage>713</lpage>. <pub-id pub-id-type="doi">10.1038/195712b0</pub-id> <pub-id pub-id-type="pmid">13919574</pub-id></citation></ref>
<ref id="B53"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tischler</surname> <given-names>H.</given-names></name> <name><surname>Wolfus</surname> <given-names>S.</given-names></name> <name><surname>Friedman</surname> <given-names>A.</given-names></name> <name><surname>Perel</surname> <given-names>E.</given-names></name> <name><surname>Pashut</surname> <given-names>T.</given-names></name> <name><surname>Lavidor</surname> <given-names>M.</given-names></name><etal/></person-group> (<year>2011</year>). <article-title>Mini-coil for magnetic stimulation in the behaving primate.</article-title> <source><italic>J. Neurosci. Methods</italic></source> <volume>194</volume> <fpage>242</fpage>&#x2013;<lpage>251</lpage>. <pub-id pub-id-type="doi">10.1016/j.jneumeth.2010.10.015</pub-id> <pub-id pub-id-type="pmid">20974177</pub-id></citation></ref>
<ref id="B54"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ueno</surname> <given-names>S.</given-names></name> <name><surname>Sekino</surname> <given-names>M.</given-names></name></person-group> (<year>2021</year>). <article-title>Figure-eight coils for magnetic stimulation: from focal stimulation to deep stimulation.</article-title> <source><italic>Front. Hum. Neurosci.</italic></source> <volume>15</volume>:<fpage>805971</fpage>. <pub-id pub-id-type="doi">10.3389/fnhum.2021.805971</pub-id> <pub-id pub-id-type="pmid">34975440</pub-id></citation></ref>
<ref id="B55"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ueno</surname> <given-names>S.</given-names></name> <name><surname>Tashiro</surname> <given-names>T.</given-names></name> <name><surname>Harada</surname> <given-names>K.</given-names></name></person-group> (<year>1988</year>). <article-title>Localized stimulation of neural tissues in the brain by means of a paired configuration of time-varying magnetic fields.</article-title> <source><italic>J. Appl. Phys.</italic></source> <volume>64</volume> <fpage>5862</fpage>&#x2013;<lpage>5864</lpage>. <pub-id pub-id-type="doi">10.1063/1.342181</pub-id></citation></ref>
<ref id="B56"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>van den Munckhof</surname> <given-names>P.</given-names></name> <name><surname>Contarino</surname> <given-names>M. F.</given-names></name> <name><surname>Bour</surname> <given-names>L. J.</given-names></name> <name><surname>Speelman</surname> <given-names>J. D.</given-names></name> <name><surname>de Bie</surname> <given-names>R. M.</given-names></name> <name><surname>Schuurman</surname> <given-names>P. R.</given-names></name></person-group> (<year>2010</year>). <article-title>Postoperative curving and upward displacement of deep brain stimulation electrodes caused by brain shift.</article-title> <source><italic>Neurosurgery</italic></source> <volume>67</volume> <fpage>49</fpage>&#x2013;<lpage>53; discussion53&#x2013;44</lpage>. <pub-id pub-id-type="doi">10.1227/01.NEU.0000370597.44524.6D</pub-id></citation></ref>
<ref id="B57"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Walsh</surname> <given-names>V.</given-names></name> <name><surname>Pascual-Leone</surname> <given-names>A.</given-names></name></person-group> (<year>2003</year>). <source><italic>Transcranial magnetic stimulation: A neurochronometrics of mind.</italic></source> <publisher-loc>Cambridge, MA</publisher-loc>: <publisher-name>The MIT Press</publisher-name>.</citation></ref>
<ref id="B58"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yamamoto</surname> <given-names>K.</given-names></name> <name><surname>Miyawaki</surname> <given-names>Y.</given-names></name> <name><surname>Saitoh</surname> <given-names>Y.</given-names></name> <name><surname>Sekino</surname> <given-names>M.</given-names></name></person-group> (<year>2016</year>). <article-title>Improvement in efficiency of transcranial magnetic stimulator coil by combination of iron core plates laminated in different directions.</article-title> <source><italic>IEEE Trans. Magn.</italic></source> <volume>52</volume>:<fpage>5100504</fpage>. <pub-id pub-id-type="doi">10.1109/TMAG.2016.2514321</pub-id></citation></ref>
<ref id="B59"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname> <given-names>X.</given-names></name> <name><surname>Zou</surname> <given-names>H.</given-names></name> <name><surname>Jung</surname> <given-names>G.</given-names></name> <name><surname>Richard</surname> <given-names>G.</given-names></name> <name><surname>Linke</surname> <given-names>S. J.</given-names></name> <name><surname>Ader</surname> <given-names>M.</given-names></name><etal/></person-group> (<year>2013</year>). <article-title>Nonneuronal control of the differential distribution of myelin along retinal ganglion cell axons in the mouse.</article-title> <source><italic>Invest. Ophthalmol. Vis. Sci.</italic></source> <volume>54</volume> <fpage>7819</fpage>&#x2013;<lpage>7827</lpage>. <pub-id pub-id-type="doi">10.1167/iovs.13-12596</pub-id> <pub-id pub-id-type="pmid">24222305</pub-id></citation></ref>
<ref id="B60"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ye</surname> <given-names>H.</given-names></name></person-group> (<year>2022</year>). <article-title>Finding the location of axonal activation by a miniature magnetic coil.</article-title> <source><italic>Front. Comput. Neurosci.</italic></source> <volume>16</volume>:<fpage>932615</fpage>. <pub-id pub-id-type="doi">10.3389/fncom.2022.932615</pub-id> <pub-id pub-id-type="pmid">35847967</pub-id></citation></ref>
<ref id="B61"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ye</surname> <given-names>H.</given-names></name> <name><surname>Barrett</surname> <given-names>L.</given-names></name></person-group> (<year>2021</year>). <article-title>Somatic inhibition by microscopic magnetic stimulation.</article-title> <source><italic>Sci. Rep.</italic></source> <volume>11</volume>:<fpage>13591</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-021-93114-x</pub-id> <pub-id pub-id-type="pmid">34193906</pub-id></citation></ref>
<ref id="B62"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ye</surname> <given-names>H.</given-names></name> <name><surname>Steiger</surname> <given-names>A.</given-names></name></person-group> (<year>2015</year>). <article-title>Neuron matters: Electric activation of neuronal tissue is dependent on the interaction between the neuron and the electric field.</article-title> <source><italic>J. Neuroeng. Rehabil.</italic></source> <volume>12</volume>:<fpage>65</fpage>. <pub-id pub-id-type="doi">10.1186/s12984-015-0061-1</pub-id> <pub-id pub-id-type="pmid">26265444</pub-id></citation></ref>
<ref id="B63"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ye</surname> <given-names>H.</given-names></name> <name><surname>Chen</surname> <given-names>V. C.</given-names></name> <name><surname>Helon</surname> <given-names>J.</given-names></name> <name><surname>Apostolopoulos</surname> <given-names>N.</given-names></name></person-group> (<year>2020</year>). <article-title>Focal suppression of epileptiform activity in the hippocampus by a high-frequency magnetic field.</article-title> <source><italic>Neuroscience</italic></source> <volume>432</volume> <fpage>1</fpage>&#x2013;<lpage>14</lpage>. <pub-id pub-id-type="doi">10.1016/j.neuroscience.2020.02.018</pub-id> <pub-id pub-id-type="pmid">32105740</pub-id></citation></ref>
<ref id="B64"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname> <given-names>D. Y.</given-names></name> <name><surname>Cringle</surname> <given-names>S. J.</given-names></name> <name><surname>Balaratnasingam</surname> <given-names>C.</given-names></name> <name><surname>Morgan</surname> <given-names>W. H.</given-names></name> <name><surname>Yu</surname> <given-names>P. K.</given-names></name> <name><surname>Su</surname> <given-names>E. N.</given-names></name></person-group> (<year>2013</year>). <article-title>Retinal ganglion cells: Energetics, compartmentation, axonal transport, cytoskeletons and vulnerability.</article-title> <source><italic>Prog. Retin. Eye Res.</italic></source> <volume>36</volume> <fpage>217</fpage>&#x2013;<lpage>246</lpage>. <pub-id pub-id-type="doi">10.1016/j.preteyeres.2013.07.001</pub-id> <pub-id pub-id-type="pmid">23891817</pub-id></citation></ref>
<ref id="B65"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zaaimi</surname> <given-names>B.</given-names></name> <name><surname>Ruiz-Torres</surname> <given-names>R.</given-names></name> <name><surname>Solla</surname> <given-names>S. A.</given-names></name> <name><surname>Miller</surname> <given-names>L. E.</given-names></name></person-group> (<year>2013</year>). <article-title>Multi-electrode stimulation in somatosensory cortex increases probability of detection.</article-title> <source><italic>J. Neural Eng.</italic></source> <volume>10</volume>:<fpage>056013</fpage>. <pub-id pub-id-type="doi">10.1088/1741-2560/10/5/056013</pub-id></citation></ref>
<ref id="B66"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>D.</given-names></name> <name><surname>Cui</surname> <given-names>X.</given-names></name> <name><surname>Zheng</surname> <given-names>J.</given-names></name> <name><surname>Zhang</surname> <given-names>S.</given-names></name> <name><surname>Wang</surname> <given-names>M.</given-names></name> <name><surname>Lu</surname> <given-names>W.</given-names></name><etal/></person-group> (<year>2021</year>). <article-title>Neurosurgical robot-assistant stereoelectroencephalography system: Operability and accuracy.</article-title> <source><italic>Brain Behav.</italic></source> <volume>11</volume>:<fpage>e2347</fpage>. <pub-id pub-id-type="doi">10.1002/brb3.2347</pub-id> <pub-id pub-id-type="pmid">34520631</pub-id></citation></ref>
</ref-list>
</back>
</article>
