Theoretical Analysis of Transcranial Magneto-Acoustical Stimulation with Hodgkin-Huxley Neuron Model

Transcranial magneto-acoustical stimulation (TMAS) is a novel stimulation technology in which an ultrasonic wave within a magnetostatic field generates an electric current in an area of interest in the brain to modulate neuronal activities. As a key part of the neural network, neurons transmit information in the nervous system. However, the effect of TMAS on the neuronal firing pattern remains unknown. To address this problem, we investigated the stimulatory mechanism of TMAS on neurons, by using a Hodgkin-Huxley neuron model. The simulation results indicated that the magnetostatic field intensity and ultrasonic power affect the amplitude and interspike interval of neuronal action potential under a continuous wave ultrasound. The simulation results also showed that the ultrasonic power, duty cycle and repetition frequency can alter the firing pattern of neural action potential under pulsed wave ultrasound. This study may help to reveal and explain the biological mechanism of TMAS and to provide a theoretical basis for TMAS in the treatment or rehabilitation of neuropsychiatric disorders.

Introduction

Transcranial magnetic stimulation (TMS), a noninvasive brain stimulation tool, has been used for treating and rehabilitating neurological and psychiatric disorders (Hallett, 2000; Golestanirad et al., 2012; Laakso and Hirata, 2012). However, the spatial resolution of TMS is more than several centimeters because the magnetic field cannot be effectively concentrated (Bystritsky et al., 2011). Additionally, the stimulation depth of TMS is limited because the alternating magnetic fields obey Laplace's Equation ∇²φ = 0, where φ is the electric potential (Norton, 2003).

Compared with TMS, transcranial focused ultrasound stimulation (tFUS) is a noninvasive method for brain stimulation that can complete a stimulation at a high millimeter-resolution level (< 2 mm; Tufail et al., 2010). More than 40 years ago, Fry et al and Ballantine et al have demonstrated that tFUS had neuromodulatory potential (Fry, 1958; Ballantine et al., 1960). Its effectiveness in neurostimulation has been confirmed in mice, rats, monkeys and humans (Tufail et al., 2010; Deffieux et al., 2013; King et al., 2014; Legon et al., 2014). tFUS can stimulate intact brain circuits of mice and can modulate visuomotor behavior in monkeys as well as the activity of the primary somatosensory cortex in humans.

In this paper, we study the effect of transcranial magneto-acoustical stimulation (TMAS) on the brain. TMAS is a novel brain stimulation technology that was proposed by Norton in 2003 (Norton, 2003). In Norton's paper, this method was analyzed to stimulate nerve tissue, and the calculation results showed that it could be used to locally stimulate active tissue, according to Maxwell's theory. TMAS can generate an electric current in a static magnetic field, using ultrasonic waves to stimulate brain tissues. The spatial resolution of TMAS is determined by the size of the ultrasonic spot, which is approximately 2 mm in diameter. Because the ultrasound and the magnetostatic field strength have good penetration depths, TMAS possesses a high spatial resolution and high penetration for brain stimulation.

The neuron is the basic unit of information transmission in the nervous system. A key neuroscience research issue has been determining how information from external stimuli are encoded and transmitted by neurons (Prescott et al., 2006). Because different neuronal firing patterns carry different stimulation signals, the understanding the firing pattern is very important in research on neural information coding. The mechanism of the neuronal firing pattern and the change of different firing patterns are important to the nervous system. Different external stimuli induce different neuronal firing rhythms and physiological effects. The exploration of the neuronal firing pattern and its underlying mechanisms are of great significance to the study of TMAS and its clinical application.

The relationship between the changes in neuronal firing patterns and TMAS has yet to be determined. Here, we investigated the effect of TMAS on the amplitude (AMP), interspike interval (ISI: the interval of adjacent action potentials per burst) and firing rate (FR: the firing number of the action potential per stimulation cycle) of neuronal action potentials, on the basis of a Hodgkin-Huxley (H-H) neuron model.

Methods

Principle of TMAS

The principle of TMAS is that the action of a focused ultrasonic wave moves charged ions in the nerve tissue. Because a magnetostatic field is perpendicular to the movement direction of the charged ions, the Lorentz force on the ions can be induced in the tissue (Ammari et al., 2015; Grasland-Mongrain et al., 2015). The Lorentz force separates the positive and negative ions in opposite directions, thereby forming the electric current I_ext to stimulate neurons.

In our study, standard Cartesian coordinate axes were used. We assumed that the pressure waves were longitudinal and propagated along the z axis and that the magnetostatic field was along the x axis, thereby placing the current density along the y axis (Figures 1A,B).

FIGURE 1

Figure 1. (A) The usual Cartensian co-ordinate axes, the pressure waves longitudinal and propagating along z axis, and the magetostatic field along x axis, the created electric current density along y axis. (B) The schematic of TMAS principle.

In the Cartesian coordinate axes, a longitudinal pressure wave propagating along the z axis obeys the classical wave Equation

$\begin{array}{ll}\frac{{\partial }^{2}u}{\partial z^2}=\frac{1}{c_0^2}\frac{{\partial }^{2}u}{\partial t^2}& (1)\end{array}$

where u is the distance of the ion from its equilibrium position. In the case of a progressive sine wave, the instantaneous speed v_z of the ion can be expressed as

$\begin{array}{ll}{v}_z=w\text{ sin }\left(\omega \text{t}-\phi \right)& (2)\end{array}$

where w is the magnitude of the speed of the ion. ω = 2πf is the angular frequency, where f is the ultrasound frequency.

The relationship between the magnitude of the speed of a fluid element and the magnitude of the instantaneous pressure P can be expressed as

$\begin{array}{ll}P=\rho {\text{c}}_{0}w& (3)\end{array}$

where ρ is the tissue density and c₀ is the ultrasound speed.

According to Montalibet's theory (Montalibet et al., 2001), the current density J_y along the y axis generated by ultrasound and magnetostatic fields in a biological medium can be expressed as

$\begin{array}{ll}{J}_{\text{y}}=\sigma \frac{w{B}_x}{1+{\text{tan}}^2\psi }\text{sin}\left(\omega t-\psi \right)& (4)\end{array}$

where σ is the conductivity of the tissue. A typical value of the conductivity of tissue is 0.5 Siemens/m (Norton, 2003). B_x is the intensity of the magnetostatic field. Ψ is the phase angle and obeys the following Equation

$\begin{array}{ll}\text{tan}\psi =\omega \tau & (5)\end{array}$

where τ is the time constant and is on the order of femtoseconds for typical electrolytes. At the ultrasound frequency generally used (200–700 kHz) in ultrasound stimulation (Tufail et al., 2011), the quantities of tanΨ and Ψ are small and negligible with the respect to unity. Therefore, the Equation can be expressed as

$\begin{array}{ll}{J}_{\text{y}}\approx \sigma w{B}_x\text{ sin }\omega \text{t}& (6)\end{array}$

The relationship between the intensity of the ultrasonic power and ultrasound pressure satisfies the following Equation (Hendee and Ritenour, 2002)

$\begin{array}{ll}\Gamma =\frac{1}{2}\frac{P^2}{\rho c_0}& (7)\end{array}$

where Γ is the intensity of the ultrasonic power.

Combining Equations (6) and (7), we obtained the following Equation

$\begin{array}{ll}{J}_{\text{y}}\approx \sigma B_x\sqrt{\frac{2\Gamma }{\rho c_0}}\text{sin}(2\pi ft)& (8)\end{array}$

If the ultrasonic wave is a cosine wave, the current density can be expressed as

$\begin{array}{ll}{J}_{\text{y}}\approx \sigma B_x\sqrt{\frac{2\Gamma }{\rho c_0}}\text{cos}(2\pi ft)& (9)\end{array}$

The fixed parameters for Equation (8) and (9) are listed in Table 1. The value of electric current density J_y, which corresponds to the electric current I_ext, can be used to stimulate the neuron and was used for simulation in the H-H neuron model.

TABLE 1

Table 1. Fixed parameters for Equation (8) and (9).

H-H Neuron Model

The H-H neuron model includes the following differential Equations (Hodgkin and Huxley, 1952)

$\begin{array}{l}{C}_m\frac{dV}{dt}=I_{ext}-[{ḡ}_{Na}{m}^{3}h(V-V_{Na})+{ḡ}_K{n}^4(V-V_K)\\ +g_L(V-V_L)]\\ \frac{dm}{dt}=\phi \left[{\alpha }_m(V)(1-m)-{\beta }_m(V)m\right]\\ \frac{dh}{dt}=\phi \left[{\alpha }_h(V)(1-h)-{\beta }_h(V)h\right]\\ \frac{dn}{dt}=\phi \left[{\alpha }_n(V)(1-n)-{\beta }_n(V)n\right]& (10)\end{array}$

where C_m is the membrane capacitance, and I_ext is the external current generated by the ultrasound and magnetic field in nerve tissue. V is the membrane potential and can be expressed as V = V_intra−V_extra. V_intra and V_extra are the intracellular and extracellular membrane potentials, respectively. φ = 3^{(T−6.3)∕10} modifies the time constants of gating variables depending on temperature T. T = 6.3°C was chosen in the simulation because we did not consider the effect of temperature on the electrical activity of neurons. m and h are the gating variables representing the activation and inactivation of the Na⁺ current, respectively. n is the gating variable representing the activation of the K⁺ current. V_Na, V_K, and V_L are the equilibrium potentials of sodium, potassium and the leak current, respectively. g_Na, g_K, and g_L are the maximum conductances of the corresponding ionic currents and are nonlinear functions of V, given by the following Equations:

$\begin{array}{l}{\alpha }_m(V)=0.1(25-V)∕\left[\text{exp}((25-V)∕10)-1\right]\\ {\beta }_m(V)=4\text{ exp}(-V∕18)\\ {\alpha }_h(V)=0.07\text{ exp}(-V∕20)\\ {\beta }_h(V)=1∕\left[\text{exp}((-V+30)∕10)+1\right]\\ {\alpha }_n(V)=0.01(10-V)∕\left[\text{exp}((10-V)∕10)-1\right]\\ {\beta }_n(V)=0.125\text{ exp}(-V∕80)& (11)\end{array}$

The fixed parameters used for the simulation are listed in Table 2. To make the resting potential equal to zero in the H-H neuron model, the value of the membrane potential was shifted by 65 mV. The initial values of V, m, h, n are 0, 0.053, 0.596, and 0.317, respectively. The simulation was performed with Matlab Simulink software (MathWorks, USA).

TABLE 2

Table 2. Fixed parameters for Hodgkin-Huxley neuron model.

The Fundamental Wave of Ultrasound

We assumed that the fundamental wave of ultrasound is a cosine wave. The simulation results are shown in Figure 2. When ultrasound and the magnetic field are under the same intensity, we found that the stimulation cannot make the neurons generate sodium current and action potentials if the fundamental wave of ultrasound is a sine wave (sin(2πft)), (Figure 2A). If the fundamental wave of ultrasound is a sine wave with offset (sin(2πft) + 1), the neurons are able to generate sodium current and action potentials with stimulation (Figure 2B). According to Equation (9), if the fundamental wave of ultrasound is a cosine wave, the electric current changes with the cosine, the sodium channel opens and then closes in a very short period of time ($\frac{{10}^{-5}}{7}~\frac{{10}^{-5}}{2}s$), and there is not enough time to generate an action potential. If the fundamental wave of ultrasound is a sine wave with offset, there is no negative phase in the electric current, the sodium channel can be opened for a long period of time, and the action potential can be generated.

FIGURE 2

Figure 2. (A) The fundamental wave of ultrasound with sine wave (top), the sodium current (center), the neuronal action potentials (bottom). (B) The fundamental wave of ultrasound with sine wave with offset (top), the sodium current (center), the neuronal action potentials (bottom).

Ultrasonic Types and Parameters

If the intensity of the magnetostatic field is constant, the current density of the nerve tissue is a function of the ultrasonic power, according to Equation (9). In this study, two types of ultrasound were used in the simulation. The first type, shown in Figure 3A, is continuous wave ultrasound, and the fundamental wave of ultrasound is a sine wave with offset (sin(2πft) + 1). The second type is pulsed ultrasound, with the modulation of a square wave and continuous wave (Figure 3B). It obeys the following Equation:

$\begin{array}{ll}x\left(t\right)=\{\begin{array}{l}1\left(n-1\right)\frac{1}{RF}<t\le \left[\left(n-1\right)+DC\right]\frac{1}{RF},n=1,2,3\dots \\ 0Others\\ \end{array}& (12)\end{array}$

where RF is the repetition frequency and DC is the duty cycle.

FIGURE 3

Figure 3. Two types of ultrasound and the related electric current. (A) Continuous wave ultrasound. (B) Pulsed ultrasound with the modulation of square wave and continuous wave.

Results

Effect of TMAS under Continuous Wave Ultrasound on the Firing Pattern of Neuronal Action Potential

Magnetostatic Field Intensity

First, to evaluate the effect of TMAS on the firing pattern of neuronal action potential, the action potential was simulated on the basis of the H-H neuron model under continuous wave ultrasound with various magnetostatic field intensities, from 0.5 to 7 Teslas (T). The values of ultrasonic power and ultrasound frequency were 3 Watt/cm² (W/cm²) and 500 kHz, respectively. Figures 4A–D shows the waveforms of action potentials corresponding to magnetostatic field intensities of 0.5, 1, 4, and 7 T, respectively. Additionally, the AMP and ISI of the action potential were calculated to quantitatively analyze the effect of the magnetostatic field potential on the neuronal firing pattern. The AMP of the action potential in relation to the magnetostatic field intensity is shown in Figure 4E. The AMP decreased with the increase in magnetostatic field intensity. Figure 4F shows the ISI of the action potential in relation to ultrasonic power. The results showed that the ISI significantly decreased with increasing magnetostatic field intensity from 0.5 to 7 T.

FIGURE 4

Figure 4. (A–D) Waveforms of neuronal action potentials generated by TMAS under continuous wave ultrasound with different magnetostatic field intensity, (A) 0.5 T, (B) 1 T, (C) 4 T, (D) 7 T. (E–F). The AMP and ISI of action potentials vs. magnetostatic field intensity, (E) AMP, (F) ISI.

Ultrasonic Power

Next, we evaluated the effect of TMAS on the firing pattern of neuronal action potential under continuous wave ultrasound with various ultrasonic powers. According to the Hopf bifurcation theorem (Wang et al., 2004), if the current density is greater than 9.78 μA/cm² (corresponding to an ultrasonic power of 0.73 W/cm²), the neuron can generate periodic action potential. Ultrasonic powers from 1 to 100 W/cm² were used in the simulation to generate action potential. These values were less than 190 W/cm², which is the maximum recommended limit for diagnostic imaging applications (Nyborg, 2001). The values of the magnetostatic field potential and the ultrasound frequency were 3 T and 500 kHz, respectively. Figures 5A–E shows the waveforms of action potentials with ultrasonic powers of 1, 10, 60, and 100 W/cm², respectively. The results showed that there was a decrease in the AMP of the action potential, and the ISI decreased as the ultrasonic power increased. The quantitatively calculated results showed that the AMP significantly decreased with the increase in ultrasonic power (Figure 5E). Figure 5F shows the ISI of the action potentials in relation to ultrasonic power. The results showed a dramatically shortened ISI with increasing ultrasonic powers from 1 to 100 W/cm², after which the narrowing of the ISI range was very limited.

FIGURE 5

Figure 5. (A–D) Waveforms of neuronal action potentials generated by TMAS under continuous wave ultrasound with different ultrasonic powers, (A) 1 W/cm², (B) 10 W/cm², (C) 60 W/cm², (D) 100 W/cm². (E,F). The AMP and ISI of action potentials vs. ultrasonic powers, (E) AMP, (F) ISI.

Ultrasound Frequency

Finally, we chose ultrasound frequencies from 200 to 700 kHz to evaluate the effect of TMAS on the neuronal firing pattern under continuous wave ultrasound. The magnetostatic field intensity and ultrasonic power were 3 T and 3 W/cm², respectively. The action potentials with ultrasound frequencies of 200, 500, and 700 kHz are shown in Figures 6A–C, respectively. The AMP and ISI of the action potential in relation to ultrasound frequency are shown in Figures 6D,E. The results showed that there were no significant changes in the AMP and ISI of the action potential.

FIGURE 6

Figure 6. (A–C) Waveforms of neuronal action potentials generated by TMAS under continuous wave ultrasound with ultrasound frequency, (A) 200 kHz, (B) 500 kHz, (C) 700 kHz. (D,E) The AMP and ISI of action potentials vs. ultrasound frequency, (D) AMP, (E) ISI.

Effect of TMAS under Pulsed Ultrasound on the Firing Pattern of Neuronal Action Potential

Ultrasonic Power

First, to evaluate the effect of ultrasonic power on the firing pattern of neuronal action potential, the action potential was simulated under pulsed wave ultrasound with various ultrasonic powers. The values of the magnetostatic field strength, ultrasound frequency, duty cycle (DC), and repetition frequency (RF) were 3 T, 500 kHz, 50%, and 10 Hz, respectively. Figures 7A–E shows the waveforms of action potentials with ultrasonic powers of 1, 10, 30, 60, and 100 W/cm², respectively. The neuronal firing form was a cluster discharge with the characteristics of a periodic variation and equal time intervals between the clusters. The number of pulses per cluster increased with the increase in ultrasonic power. The AMP, ISI, and FR of the action potential were quantitatively calculated, and the results are shown in Figures 7F–H. The results showed that the AMP significantly decreased with the increase in ultrasonic power. We also found that the ISI decreased with the increase in ultrasonic power from 1 to 100 W/cm². The FR gradually increased in a multilevel ladder shape. The above results indicate that the firing pattern of neuronal action potential is selective for ultrasonic powers.

FIGURE 7

Figure 7. (A–C) Waveforms of neuronal action potentials generated by TMAS under pulsed ultrasound with different ultrasound power, (A) 1 W/cm², (B) 10 W/cm², (C) 30 W/cm², (D) 60 W/cm², (E) 100 W/cm². (F–H). The AMP, ISI and FR of action potentials vs. ultrasonic powers, (F) AMP, (G) ISI, (H) FR.

DC of Ultrasound

Next, DCs from 5 to 95% in intervals of 5% were selected to evaluate the effect of TMAS on the neuronal firing pattern under pulsed ultrasound. The values of the magnetostatic field strength, ultrasound frequency, ultrasonic power, and RF were 3 T, 500 kHz, 3 W/cm², and 10 Hz, respectively. The action potentials with DCs of 5, 50, and 95% are shown in Figures 8A–C, respectively. The firing form of the neuron was also a cluster discharge. The time interval between the clusters decreased with the increase in repetition frequency. The AMP, ISI, and FR of the action potential in relation to the duty cycle are shown in Figures 8D–F, respectively. The results showed a weak decrease in the AMP and ISI with the increase in DC. We also found that the FR gradually decreased with the increase in DC.

FIGURE 8

Figure 8. (A–C) Waveforms of neuronal action potentials generated by TMAS under pulsed ultrasound with different DC, (A) 5 %, (B) 50 %, (C) 95 %. (D–F). The AMP, ISI, and FR of action potentials vs. DC, (D) AMP, (E) ISI, (F) FR.

RF of Ultrasound

Finally, we used RFs from 1 to 100 Hz to evaluate the effect of TMAS on the neuronal firing pattern under pulsed ultrasound. The magnetostatic field strength, ultrasound frequency, ultrasonic power and DC were 3 T, 500 kHz, 3 W/cm², and 50%, respectively. The action potentials with repetition frequencies of 5, 20, and 100 Hz are shown in Figures 9A–C, respectively. The AMP, ISI, and FR of the action potential compared to repetition frequency are shown in Figures 9D–F, respectively. The results showed that both the AMP and ISI gradually decreased in a multilevel ladder shape route with an increase in repetition frequency. We also clearly observed that the FR significantly decreased as the RF increased from 1 to 100 Hz, after which the narrowing of the FR range was very limited.

FIGURE 9

Figure 9. (A–C) Waveforms of neuronal action potentials generated by TMAS under pulsed ultrasound with different RF, (A) 5 Hz, (B) 20 Hz, (C) 100 Hz. (D–F). The AMP, ISI, and FR of action potentials vs. RF, (D) AMP, (E) ISI, (F) FR.

Discussion and Conclusion

In this study, we investigated the effect of TMAS on the neuronal firing pattern. Our results showed that neurons can generate a rich firing pattern with TMAS. The simulation results showed that neurons have a high sensitivity to changes in ultrasound parameters in TMAS. Small changes in ultrasound parameters can lead to changes in the neuronal firing pattern, thus affecting the neural information coding. These findings lay the foundation for further studies on the effects of TMAS on the neural network and reveal the underlying mechanisms of TMAS in the biological neural system.

Regarding the mechanisms of tFUS and TMAS, a previous study has shown that the mechanism of tFUS is that ultrasound-induced cavitation of these nanometric bilayer sonophores can induce a complex mechanoelectrical interplay that leads to excitation, primarily through the effect of currents induced by membrane capacitance changes (Plaksin et al., 2014). The mechanism of TMAS is that an electric current is generated by the static magnetic field and ultrasonic waves in nerve tissues to modulate neuronal activity. Therefore, the mechanisms of ultrasound stimulation and magneto-acoustical stimulation are different.

Regarding the effect of the parameters of tFUS and TMAS on stimulation results, King et al have studied the effect of different types of stimuli (continuous-wave stimuli and pulsed stimuli) and ultrasonic parameters (acoustic intensity, sonication duration, pulse repetition frequency, etc.) on the stimulation success of ultrasound-induced in vivo neurostimulation (King et al., 2014). In our study, we analyzed the effects of magnetic intensity and types of ultrasound (continuous and pulsed ultrasound) with various ultrasonic parameters on the neuronal firing pattern. The simulation results showed that the type of ultrasound, the magnetic intensity, and the ultrasonic parameters induce different neuronal firing patterns. Our results may provide important theoretical guidance for the observation of neuron firing activity in future experiments.

The mechanical index (MI) is the value of the related biological mechanical effects and can be used to estimate the potential biomechanical effects. The thermal index (TI) is used to indicate the temperature rise of the tissue caused by the absorption effect of the ultrasonic energy. The MI and TI formulas presented in the body of the “Standard for real-time display of thermal and mechanical acoustic output indices on diagnostic ultrasound equipment” (Abbott, 1999) are very important for evaluating the MI and TI for TMAS. According to IEC 60601-2-37 standard regulation, the MI can be expressed as $MI=\frac{p_{r.3}}{\sqrt{f}{C}_{MI}}$, where C_MI equals 1 MPa·MHz^1∕2, p_r.3 is the attenuated peak-rarefactional pressure in MPa, and f is the acoustic working frequency in MHz. The TI for unscanned modes can be expressed as $TI=\frac{Wf}{C_{TIS}}$, where C_TIS is 210 mW·MHz. In a previous study, an ultrasonic power of 3 W/cm² has been used for ultrasound brain modulation in humans. In our study, the simulation results showed that neurons can also generate an action potential with an ultrasonic power of 3 W/cm², an ultrasound frequency of 0.5 MHz and a magnetic intensity of 3 T. If the ultrasonic power is 3 W/cm², the ultrasound frequency is 0.5 MHz, and the diameter of the ultrasonic spot is 3 mm. With a p_r.3 of 0.2275 MPa, an ultrasound spot area of 7.065 × 10⁻² cm² and a W of 212 mW, we obtained an MI of 0.32 and a TI of 0.505, which were less than the MI and TI threshold (equal 1) that can induce any tissue injury.

Previous studies have shown that abnormal neuronal firing patterns cause some diseases of the nervous system (Valdez et al., 2013). The pathogenesis of epilepsy is excessive synchronization of brain neuronal discharge, resulting in a brain function imbalance (Ding et al., 2006). Our study shows that TMAS can alter the neural firing pattern, thus indicating that it has great potential to intervene in an epileptic seizure. We also found that the ultrasonic parameters can affect the neuronal firing pattern; hence, we can change these parameters during the treatment of epileptic seizures to achieve an optimal therapeutic effect.

In summary, the current investigation demonstrates that TMAS can alter the neuronal firing pattern in two ways: (i) the magnetostatic field intensity and ultrasonic power are able to affect the AMP and ISI of neuronal action potential under continuous wave ultrasound, and (ii) the ultrasonic power, DC and RF can alter the firing pattern of neural action potential under pulsed ultrasound. Our results suggest that TMAS has potential as a powerful noninvasive method to interfere with brain rhythms for diagnostic and therapeutic purposes.

Author Contributions

YY and XL designed and coordinated the study, YY and YC carried out numerical implementation of the TMAS, done the simulation and drafted the manuscript. All authors gave final approval for publication.

Conflict of Interest Statement

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.

Acknowledgments

This research was supported by National Natural Science Foundation of China (61503321, 61273063), Natural science foundation of Hebei Province (F2014203161).

Abbreviations

TMAS, Transcranial magneto-acoustical stimulation; H-H, Hodgkin-Huxley; TMS, Transcranial magnetic stimulation; tFUS, Transcranial focused ultrasound stimulation; AMP, amplitude; ISI, interspike internal; FR, firing rate; DC, duty cycle; RF, repetition frequency; MI, mechanical index; TI, thermal index.

References

Abbott, J. G. (1999). Rationale and derivation of MI and TI-a review. Ultrasound Med. Biol. 25, 431–441. doi: 10.1016/S0301-5629(98)00172-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Ammari, H., Grasland-Mongrain, P., Millien, P., Seppecher, L., and Seo, J. K. (2015). A mathematical and numerical framework for ultrasonically-induced Lorentz force electrical impedance tomography. J. de Mathematiques Pures et Appliquées. 103, 1390–1409. doi: 10.1016/j.matpur.2014.11.003

CrossRef Full Text | Google Scholar

Ballantine, H. T. Jr., Bell, E., and Manlapaz, J. (1960). Progress and problems in the neurological applications of focused ultrasound. J. Neurosurgery. 17, 858–876. doi: 10.3171/jns.1960.17.5.0858

PubMed Abstract | CrossRef Full Text | Google Scholar

Bystritsky, A., Korb, A. S., Douglas, P. K., Cohen, M. S., Melega, W. P., Mulgaonkar, A. P., et al. (2011). A review of low-intensity focused ultrasound pulsation. Brain Stimul. 4, 125–136. doi: 10.1016/j.brs.2011.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Deffieux, T., Younan, Y., Wattiez, N., Tanter, M., Pouget, P., and Aubry, J.-F. (2013). Low-intensity focused ultrasound modulates monkey visuomotor behavior. Curr. Biol. 23, 2430–2433. doi: 10.1016/j.cub.2013.10.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, L., Worrell, G. A., Lagerlund, T. D., and He, B. (2006). 3D source localization of interictal spikes in epilepsy patients with MRI lesions. Phys. Med. Biol. 51, 4047–4062. doi: 10.1088/0031-9155/51/16/011

PubMed Abstract | CrossRef Full Text | Google Scholar

Fry, W. J. (1958). Use of intense ultrasound in neurological surgery research. Am. J. Phys. Med. 37, 143–147. doi: 10.1097/00002060-195806000-00009

CrossRef Full Text | Google Scholar

Golestanirad, L., Rouhani, H., Elahi, B., Shahim, K., Chen, R., Mosig, J. R., et al. (2012). Combined use of transcranial magnetic stimulation and metal electrode implants: a theoretical assessment of safety considerations. Phys. Med. Biol. 57, 7813–7827. doi: 10.1088/0031-9155/57/23/7813

PubMed Abstract | CrossRef Full Text | Google Scholar

Grasland-Mongrain, P., Destrempes, F., Mari, J. M., Souchon, R., Catheline, S., Chapelon, J. Y., et al. (2015). Acousto-electrical speckle pattern in Lorentz force electrical impedance tomography. Phys. Med. Biol. 60:3747. doi: 10.1088/0031-9155/60/9/3747

PubMed Abstract | CrossRef Full Text | Google Scholar

Hallett, M. (2000). Transcranial magnetic stimulation and the human brain. Nature 406, 147–150. doi: 10.1038/35018000

PubMed Abstract | CrossRef Full Text | Google Scholar

Hendee, W. R., and Ritenour, E.R.. (2002). Medical Imaging Physics, 4th Edn. New York, NY: John Wiley & Sons.

Google Scholar

Hodgkin, A. L., and Huxley, A. F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. 117, 500–544. doi: 10.1113/jphysiol.1952.sp004764

PubMed Abstract | CrossRef Full Text | Google Scholar

King, R. L., Brown, J. R., and Pauly, K. B. (2014). Localization of ultrasound-induced in vivo neurostimulation in the mouse model. Ultrasound Med. Biol. 40, 1512–1522. doi: 10.1016/j.ultrasmedbio.2014.01.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Laakso, I., and Hirata, A. (2012). Fast multigrid-based computation of the induced electric field for transcranial magnetic stimulation. Phys. Med. Biol. 57, 7753–7765. doi: 10.1088/0031-9155/57/23/7753

PubMed Abstract | CrossRef Full Text | Google Scholar

Legon, W., Sato, T. F., Opitz, A., Mueller, J., Barbour, A., Williams, A., et al. (2014). Transcranial focused ultrasound modulates the activity of primary somatosensory cortex in humans. Nat. Neurosci. 3620, 1–8. doi: 10.1038/nn.3620

CrossRef Full Text | Google Scholar

Montalibet, A., Jossinet, J., Matias, A., and Cathignol, D. (2001). Electric current generated by ultrasonically induced Lorentz force in biological media. Med. Biol. Eng. Comput. 39, 15–20. doi: 10.1007/BF02345261

PubMed Abstract | CrossRef Full Text | Google Scholar

Norton, S. J. (2003). Can ultrasound be used to stimulate nerve tissue. Biomed. Eng. Online 2, 1–9. doi: 10.1186/1475-925X-2-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Nyborg, W. L. (2001). Biological effects of ultrasound: development of safety guidelines. Part II: general review. Ultrasound Med. Biol. 27, 301–333. doi: 10.1016/S0301-5629(00)00333-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Plaksin, M., Shoham, S., and Kimmel, E. (2014). Intramembrane cavitation as a predictive bio-piezoelectric mechanism for ultrasonic brain stimulation. Phys. Rev. X. 4:011004. doi: 10.1103/PhysRevX.4.011004

CrossRef Full Text | Google Scholar

Prescott, S. A., Ratté, S., De Koninck, Y., and Sejnowski, T. J. (2006). Nonlinear interaction between shunting and adaptation controls a switch between integration and coincidence detection in pyramidal neurons. J. Neurosci. 26, 9084–9097. doi: 10.1523/JNEUROSCI.1388-06.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Tufail, Y., Matyushov, A., Baldwin, N., Tauchmann, M. L., Georges, J., Yoshihiro, A., et al. (2010). Transcranial pulsed ultrasound stimulates intact brain circuits. Neuron 66, 681–694. doi: 10.1016/j.neuron.2010.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Tufail, Y., Yoshihiro, A., Pati, S., Li, M. M., and Tyler, W. J. (2011). Ultrasonic neuromodulation by brain stimulation with transcranial ultrasound. Nat. Protocols. 6, 1453–1470. doi: 10.1038/nprot.2011.371

PubMed Abstract | CrossRef Full Text | Google Scholar

Valdez, A. B., Hickman, E. N., Treiman, D. M., Smith, K. A., and Steinmetz, P. N. (2013). A statistical method for predicting seizure onset zones from human single-neuron recordings. J. Neural Eng. 10:016001. doi: 10.1088/1741-2560/10/1/016001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Tsang, K. M., and Zhang, H. (2004). Hopf bifurcation in the Hodgkin-Huxley model exposed to ELF electric field. Chaos, Solitons Fractals 20, 759–764. doi: 10.1016/j.chaos.2003.08.011

CrossRef Full Text | Google Scholar