Abstract
Synapse coupling can benefit signal exchange between neurons and information encoding for neurons, and the collective behaviors such as synchronization and pattern selection in neuronal network are often discussed under chemical or electric synapse coupling. Electromagnetic induction is considered at molecular level when ion currents flow across the membrane and the ion concentration is fluctuated. Magnetic flux describes the effect of timevarying electromagnetic field, and memristor bridges the membrane potential and magnetic flux according to the dimensionalization requirement. Indeed, field coupling can contribute to the signal exchange between neurons by triggering superposition of electric field when synapse coupling is not available. A chain network is designed to investigate the modulation of field coupling on the collective behaviors in neuronal network connected by electric synapse between adjacent neurons. In the chain network, the contribution of field coupling from each neuron is described by introducing appropriate weight dependent on the position distance between two neurons. Statistical factor of synchronization is calculated by changing the external stimulus and weight of field coupling. It is found that the synchronization degree is dependent on the coupling intensity and weight, the synchronization, pattern selection of network connected with gap junction can be modulated by field coupling.
Introduction
Neurons are the basic structural and functional units of nervous system. The information processing and encoding are approached by exchange of neurotransmitter via chemical and electrical action, and the signal receiving is dependent on synapses between neurons^{1,2,3,4,5,6,7}. For example, Qin et al.^{5} discussed the biological function of autapse connection to neurons and confirmed that autapse driving can enhance selfadaption of neurons. The biological nervous system is usually very complex and hierarchical, and it contains very complex inner structures, such as multiple subsystems. Therefore, the research in internal competition and cooperation behavior of complex system has to resort to the sampling and system analysis. The related studies have very great significance on understanding the evolution of collective behaviors of complex systems, and it is also helpful to predict the dynamical behavior of control systems. For an isolated neuron, geometry and anatomical structure often account for biological function, for example, autapse formation could be associated with injury in axon thus auxiliary loops are developed to propagate the blocked signals^{8}. Moreover, the impact of different topology networks in various neuron models has been widely investigated by using bifurcation analysis. Indeed, coupling strength^{9,10}, time delay^{11}, spatial in homogeneity stimulus^{12}, and noise^{13,14} can affect the synchronization and stability of network. The network is composed of many nodes, the statistical analysis based on mean field theory and nonlinear analysis from sampled time series of observable variables are often used to detect the synchronization stability and pattern selection^{15,16}. However, it could be difficult for generic sampling and statistical analysis when network is composed of many subnetworks, but pattern selection could be effective to understand the transition in collective behaviors. On the other hand, the use of speckle graph dynamics is much helpful to study the group discharge behaviors of neuronal networks, thus the spiral wave characteristics of spatial distribution of electrical activity in the cerebral cortex^{17} can be understood. The spiral wave is selfsustained, and the developed spiral waves in cortex can change the collective behaviors in electrical activities as a stable pacemaker. The potential mechanism for emergence of spiral wave and its dynamics are discussed in ref.^{18}, which local and diffusive poisoning in ion channels^{19,20,21} can generate defects to trigger continuous waves and further occurrence of spiral waves after symmetrical breaking^{22} in the profile. For example, Hu et al.^{22} confirmed that symmetry breaking of target wave in the media accounts for emergence of multiarmed spiral wave, which can be developed from a group of spiral waves with single arm under appropriate condition, thus the potential formation mechanism of multiarmed spiral wave in the media is explained. Gong et al^{23}. discussed the influence of time delay and channel blocking on multiple coherence resonance in network composed of HodgkinHuxley neurons. Qin et al.^{24} explained the defect mechanism for autapse driving in network by using negative feedback and the wave propagation was also discussed. Xiao et al.^{25} reported the dynamics in spatiotemporal patterns in network by changing the excitability. Xu et al.^{26} explained the new mechanism of spatial coherence resonance by imposing noise and periodical pacing in different local areas in the network, and it is found that local forcinginduced wave can be broken in the area driven by noise and spiral waves are formed under appropriate noise intensity in the network. Li et al.^{27} discussed the selection and breakup of spiral wave in a coupled network driven by Gaussian colored noise. Qian et al.^{28} analyzed all oscillatory complex networks consisting of nonoscillatory nodes by using the dominant phaseadvanced driving method, the oscillation sources and wave propagation paths of the selfsustained target waves in excitable smallworld networks were explored explicitly.
In fact, spatial patterns such as spiral wave, target wave and Turing patterns can also be found in Reactiondiffusion systems and neuronal networks as well. Pattern formation in nonlinear complex systems is one of the central problems of the natural, social, and technological sciences. Since Turing firstly proposed reactiondiffusion (RD) theory to model the range of spatial patterns observed in the developing embryo^{29}, RD models have been studied extensively to explain patterns in both epidemiology and ecosystems^{30}. It is found that cross diffusion^{31}, time delay^{32,33} and functional response^{34} can induce the appearance of stationary patterns. It is confirmed that the standard multiplescale analysis and exact Turing regions can be obtained^{35} to develop kinds of spatial patterns. In fact, spiral wave, target wave or patched invasion can be observed^{36,37} when the diffusion coefficients are the same (Turing theory does not hold). Based on the data observations, it is revealed that multiple scale spatial patterns^{38}, isolation degree^{39} and spatial scaling laws^{40} play important roles on the resilience of the biological systems and thus they have implications on population protection or diseases control. In cardiac tissue and also cortex, the emergence and survival of spiral waves, continuous pulses could disturb signal propagation and information processing, thus normal electrical activities in cortex and wave propagation in cardiac tissue can be perturbed for inducing possible diseases. Therefore, it is important to apply effective schemes^{40,41} to remove and suppress spiral waves. Based on the drivensynchronization, the rotating electric field (REF) can be utilized to “smooth” a heterogeneous medium and suppress turbulence^{42}. Chen et al.^{43,44} found that the REF can induce wave emission from heterogeneity in excitable media, which may provide a potential application to utilized the existed defects on cardiac tissue as the second sources to remove pinned spiral waves^{45}.
That is, the effect of electric field on cardiac tissue can often be described by adding gradient terms in the membrane potential of media described by the ReactionDiffusion system. In the cortex and neuronal network, it is a challenge to deal with this problem that the effect of electromagnetic induction can be understood from physical view. The influence of the biological effects of electromagnetic fields on human health has been investigated, and electromagnetic radiationinduced disease and abnormality in nervous systems (such as memory loss and other symptoms) were also discussed^{46,47,48,49,50,51,52,53,54,55}. On the other hand, the dynamics of neuronal activities in presence of electromagnetic field also became attractive^{56,57,58}. In fact, by further improving the computational models for neuron^{59,60} and cardiac tissue^{61,62}, the physical effect can be described by using magnetic flux associated with timevarying electromagnetic field. For example, Lv et al.^{63} proposed a fourvariable neuron model developed from the previous HindmarshRose neuron model^{59,60} by adding the variable for magnetic flow, which memristor is also used to bridge the coupling between membrane potential and magnetic flux^{63}. Inspired by this model, the transition of electrical activities in neurons induced by electromagnetic radiation has been investigated extensively^{64,65,66,67,68}, particularly, the same scheme was carried out on cardiac tissue, and two death mechanisms^{61,62} in heart induced by electromagnetic radiation are explained. For a further survey and guidance, readers can find clues in refs^{69,70} and references therein. In most of the previous works, synapse connection and gap junction are regarded as the main bridge to benefit signal exchange and wave propagation between neurons. In this paper, we argued that field coupling could be another effective way for signal propagation in network because field coupling can realize phase synchronization between neurons^{71}. The potential mechanism is explained that each neuron can contribute to the distribution of electric field, and the electrical activity of each neuron will be affected by other neurons by triggering different electric fields in the network.
Model and scheme
The Hindmarshrose neuron model^{59,60} mainly highlights the main nonlinear dynamic characteristics of mollusk neurons, the dynamical kinetics can be approached by the ordinary differential equations (ODEs) composed of three variables shown in ref.^{59}. As proposed and discussed in refs^{63,64}, magnetic flux is used to describe the effect of electromagnetic induction^{72}, and memristor is considered to realize coupling between membrane potential and magnetic flux thus the induction field and action potential can be bridged in physical view. The dynamical neuron model is described by
where the variable x, y, z represents the membrane potential, slow current for recovery variable, and adaption current, respectively. The variable I defines the external forcing current, and the fourth variable φ describes the magnetic flux across membrane. The ρ(φ) is the memoryconductance of a magnetic fluxcontrolled memristor and it is approached as ρ(φ) = α + 3βφ^{2}, kρ(φ)x denotes induction current and k is induction coefficient. Further detailed description for parameters can be found in refs^{63,64}. The induction current is calculated as follows
That is, the effect of electromagnetic induction and change of electric field is described by induction current. The memductance of memristor is dependent on the magnetic flux and timevarying induction current is imposed on the neuron. The memristor^{73} is often used as nonlinear electric device in setting nonlinear circuits and synchronization can be approached by using memristor coupling^{74,75}. It is also effective to describe the effect of electromagnetic induction in neuronal activities. From physical view, each neuron and cell can be thought as an nonuniform charged body because the distribution of charged ions are not distributed in uniform way, and thus complex field distribution can be triggered in intracellular and extracellular space. The effect of electric field could be distinct and the timevarying electromagnetic field plays important role in changing the exchange of charged ions. In biological and nervous systems, contribution to the field distribution from each neuron could be different and shows diversity, therefore, appropriate weight is considered in the neuronal network. For simplicity, in a chainlike network, the position of neuron is marked with subscript “i”, another neuron is marked with subscript “j”, it is supposed that the contribution to field intensity is inversely proportional to the number error as Γ_{ ij } = W/│i − j│, and W is the weight value contributed to the field intensity. As a result, the network can be approached by gap junction coupling and field coupling; the dynamical equations are described by
The subscript i indicates the position of the neuron in the network, D is the coupling coefficient between adjacent neurons via electric synapse coupling, D_{0} describes the field interaction between neurons. W represents the intensity of the field effect associated with distance between neurons. According to Eq. (3), increasing the value for intensity of field coupling, positive feedback will enhance the electromagnetic induction effect, and stronger induction current will be induced to decrease the excitability of neurons. The parameters are selected as a = 1.0, b = 3.0, c = 1.0, d = 5.0, s = 4.0, α = 0.1, β = 0.02, k = 0.9, k_{1} = 0.4, k_{2} = 0.5. The field coupling is illustrated in Fig. 1.
To discern the collective behavior of neurons, a statistical factor of synchronization^{24,26,69} is calculated to find synchronization approach based on mean field theory, and it reads as follows
where N is the node number of neuronal network, < * > describes the average of variable over time, for simplicity, a transient period T = 2000 time units will be used for numerical studies. Perfect synchronization is approached when R is close to 1, and nonperfect synchronization is detected together with emergence of appropriate pattern distribution when R is close to zero.
Numerical results and discussion
In the numerical section, the fourth order RungeKutta algorithm is used to approach solutions for the dynamical equations. The chain network is composed of 100 neurons, time step is set as h = 0.01. The initial values are selected random values as (0.3 + ξ, 0.1 + ξ, 5 + ξ, 0, 0.2 + ξ), where ξ is a random number ranged from 0 to 1. Transient period for calculation is about 2000 time units, coupling coefficient D = 0.5. In order to illustrate the effect of field superposition on neuronal discharges, different external stimuli I are applied when the chain neuron W is set to 1. The spatiotemporal evolution of membrane potentials are calculated for dynamical analysis when the node position driven by external stimulus current is changed. Furthermore, the distribution for factors of synchronization on the network is calculated by selecting different external stimuli I and D_{0.} At first, the weight of electric field contribution is fixed at W = 1, the intensity of field coupling is adjusted to observe the development of spatial patterns(spatial distribution of membrane potential), the results are plotted in Fig. 2 and sampled time series for membrane potential of neuron connected to node (i = 80) are shown in Fig. 3.
The excitability of neuron mainly depends on the external forcing current, and larger external stimuli are much helpful to excite neurons and the collective behaviors of network will present oscillating state. When the subthreshold forcing is applied on the network, diversity in initial settings can trigger short transient oscillation and then the network become quiescent state occupied the network completely. By further increasing the intensity of external stimuli, the excitability of neuron is also enhanced, the quiescent neurons are waken up to present spiking and bursting states. Indeed, the field coupling contributes magnetic flux and induction current, as a result, the mode in electrical activities are controlled. As shown in Figs 2d, 3d, further increasing the intensity of field coupling, diversity in magnetic flux and induction currents are induced, as a result, each neuron shows different excitability to trigger appropriate oscillating behaviors in electrical activities.
The results in Fig. 4 confirmed that the field coupling plays important role in regulating the collective behaviors of the network, while slight difference in external stimuli can also enhance the developed state when the intensity of field coupling is highly increased. The sampled time series for some nodes found that bursting behaviors are enhanced under field coupling even stimuli is much weak, the results are shown in Fig. 5_{.}
The change in the discharge period of the neuron is caused by the change of the field interaction between each neuron and the other neurons, that is, D_{0} is positively correlated with the neuronal discharge period T. It is found that the bursting states show distinct difference and the response intervals are much different. The field coupling can drive all neurons to give appropriate response in time. The distribution for factors of synchronization is calculated in the twoparameter space, and the contribution from field coupling and external forcing is estimated, the results are shown in Fig. 6.
It is confirmed that the synchronization degree is high when external forcing current is weak that excitability is low, as a result, the effect of diversity in initials setting can be suppressed by field coupling. To discern the effect of field coupling, the intensity of electric coupling is removed as D = 0, then external forcing current and field coupling are considered, the results are shown in Figs 7, 8, 9, 10.
It is confirmed that weak excitability can’t support continuous oscillating, while field coupling can enhance the oscillating behaviors in electrical activities. The external stimuli mainly contribute the electrical activities when the intensity of field coupling is weak. Furthermore, the field coupling is enhanced to observe the evolution of collective behaviors, the results are shown in Figs 9, 10.
In the case of high excitability by applying stronger external stimuli, the oscillating behaviors and bursting states can be increased greatly under field coupling. As shown in Fig. 10, multiple modes in electrical activities are induced by further increasing the intensity of field coupling, and all the neurons can give rapid response under field coupling. Extensive numerical results confirmed synchronization can be approached when electric coupling is set as D = 0.5 when diversity in initials setting is removed. As a result, the developed states show some difference when field coupling is triggered in different transient periods. Furthermore, the effect of weight contributed to the supimposed field is discussed, the transition of synchronization factors is estimated by changing the intensity of field coupling, the results are plotted in Fig. 11.
It is found that the weight in field coupling also contributes to the synchronization approach of the network. When the weight for field coupling is much small, the synchronization factors show slight difference and synchronization mainly depends on the electric coupling. The role of weight becomes most important when field coupling is increased in intensity because the contribution from each neuron becomes distinct different than other neurons. When the external stimulus is decreased to I = 1.2, the field coupling can modulate the synchronization behavior greatly (decreases the synchronization) with increasing the intensity of field coupling. However, in case of high excitability and stronger stimulus, the oscillating behaviors and evolution of network mainly depend on the stimuli and the effect of weight shows slight difference on synchronization approach.
In a summary, the evolution of collective behaviors and synchronization degree are dependent on the electric coupling via gap junction, intensity of field coupling, weight for superimposed field, diversity in initials setting. Due to the memory effect of memristor associated with magnetic flux and effect of initials diversity becomes important and distinct. The external forcing current changes the excitability of neuron and the electrical activities of neuronal network can be changed directly. Field coupling can modulate the collective behaviors of network by inducing induction currents with diversity and diversity in excitability of network is triggered. The development and evolution of collective electrical activities of neuronal network are contributed by field coupling, electric coupling when external forcing is fixed. On the other hand, the field coupling between neurons can be considered as external electromagnetic radiation on isolated neuron. Readers can extend this study by setting other networks with different connection types when field coupling is considered.
Conclusions
An improved neuron model, which the membrane potential and magnetic flux are bridged by using memristor, is used to describe the local kinetics of chain network of neurons. The effect of electromagnetic field is described by magnetic flux. Each neuron is regarded as an charged body and is controlled by the electric field triggered by other neurons. A weight value is introduced to discern the contribution to superimposed field from each neuron, field coupling between neurons is described by exchange of magnetic flux. It is found that field coupling between neurons can change the magnetic flux and induction current, and then the excitability of neurons are changed to modulate the collective behaviors of electrical activities in neuronal network. It could give instructive clues to understand the signal encoding and exchange when synapse coupling is suppressed.
References
 1.
Pfrieger, F. Roles of glial cells in synapse development, Cell Molec. Life Sci 66, 2037 (2009).
 2.
Liu, Q., Xu, Q. & Kang, J. et al. Astrocyte activation of presynaptic metabotropic glutamate receptors modulates hippocampal inhibitory synaptic transmission. Neuron Glia Biol 1, 307 (2004).
 3.
Lichtman, J. & Colman, H. Synapse elimination and indelible memory. Neuron 25, 269–78 (2000).
 4.
Bi, G. & Poo, M. Synaptic modifications in cultured hippocampal neurons: dependence on spike timing, synaptic strength, and postsynaptic cell type. J Neurosci 18, 10464–72 (1998).
 5.
Qin, H. X., Ma, J. & Jin, W. Y. et al. Dynamics of electric activities in neuron and neurons of network induced by autapses. Sci China Technol Sci 57, 936–946 (2014).
 6.
Wang, H., Wang, L. & Yu, L. et al. Response of MorrisLecar neurons to various stimuli. Phys Rev E 83, 021915 (2011).
 7.
Yilmaz, E., Uzuntarla, M. & Ozer, M. et al. Stochastic resonance in hybrid scalefree neuronal networks. Physica A 392, 5735–5741 (2013).
 8.
Wang, C. N., Guo, S. L. & Xu, Y. et al. Formation of autapse connected to neuron and its biological function. Complexity 2017, 543673 (2017).
 9.
Shi, X. & Xi, W. Synchronization and rhythm dynamics of a neuronal network consisting of mixed bursting neurons with hybrid synapses. Int J Mod Phys B 30(16), 1650091 (2016).
 10.
Xu, Y., Jin, W. & Ma, J. Emergence and robustness of target waves in a neuronal network. Int J Mod Phys B 29(23), 1550164 (2015).
 11.
Wang, Q., Perc, M. & Duan, Z. et al. Spatial coherence resonance in delayed HodgkinHuxley neuronal neuronal networks. Int J Mod Phys B 24(9), 1201–1213 (2010).
 12.
Wang, Q., Perc, M. & Duan, Z. et al. Delayenhanced coherence of spiral waves in noisy Hodgkin–Huxley neuronal networks. Phys Lett A 372(35), 5681–5687 (2008).
 13.
Jin, W., Lin, Q., Wang, A. & Wang, C. Computer simulation of noise effects of the neighborhood of stimulus threshold for a mathematical model of homeostatic regulation of sleepwake cycles. Complexity 2017, 4797545 (2017).
 14.
Yilmaz, E., Baysal, V. & Perc, M. et al. Enhancement of pacemaker induced stochastic resonance by an autapse in a scalefree neuronal network. Sci China Technol Sci 59, 364–370 (2016).
 15.
Ma, J., Xu, Y. & Ren, G. D. et al. Prediction for breakup of spiral wave in a regular neuronal network. Nonlinear Dyn 84, 497–509 (2017).
 16.
Song, X. L., Wang, C. N. & Ma, J. et al. Collapse of ordered spatial pattern in neuronal network. Physica A 451, 95–112 (2016).
 17.
Huang, X., Xu, W. & Liang, J. et al. Spiral wave dynamics in neocortex. Neuron 68, 978–990 (2010).
 18.
Wu, X. Y. & Ma, J. The Formation Mechanism of Defects, Spiral Wave in the Network of Neurons. PLoS One 8, e55403 (2013).
 19.
Wu, X., Ma, J. & Li, F. et al. Development of spiral wave in a regular network of excitatory neurons due to stochastic poisoning of ion channels. Commun Nonlinear Sci Numer Simulat 18, 3350–3364 (2013).
 20.
Ma, J., Huang, L. & Tang, J. et al. Spiral wave death, breakup induced by ion channel poisoning on regular HodgkinHuxley neuronal networks. Commun Nonlinear Sci Numer Simulat 17, 4281–4293 (2012).
 21.
Huang, L., Ma, J. & Tang, J. et al. Transition of ordered waves in neuronal network induced by diffusive poisoning of ion channels. J Biol Syst 21, 1350002 (2013).
 22.
Hu, B., Ma, J. & Tang, J. Selection of Multiarmed Spiral Waves in a Regular Network of Neurons. Plos One 8(7), e69251 (2013).
 23.
Gong, Y., Hao, Y. & Lin, X. et al. Influence of time delay and channel blocking on multiple coherence resonance in HodgkinHuxley neuron networks, Biosyst 106(2–3), 76–81 (2011).
 24.
Qin, H. X., Wu, Y. & Wang, C. N. et al. Emitting waves from defects in network with autapses. Commun Nonlinear Sci Numer Simulat 23, 164–174 (2015).
 25.
Xiao, W. W., Gu, H. & Liu, M. Spatiotemporal dynamics in a network composed of neurons with different excitabilities and excitatory coupling. Sci China Technol Sci 59, 1943–1952 (2016).
 26.
Xu, Y., Wang, C. N. & Lv, M. et al. Local pacing, noise induced ordered wave in a 2D lattice of neurons. Neurocomputing 207, 398–407 (2016).
 27.
Li, F. & Ma, J. Selection of spiral wave in the coupled network under Gaussion colored noise. Int J Mod Phys B 27(21), 1350115 (2013).
 28.
Qian, Y., Huang, X. D. & Hu, G. et al. Structure and control of selfsustained target waves in excitable smallworld networks. Phys Rev E 81, 036101 (2010).
 29.
Turing, A. M. On the chemical basis of morphogenesis, Philos. Trans R Soc Lond B 237, 37–72 (1952).
 30.
Klausmeier, C. A. Regular and irregular patterns in semiarid vegetation. Science 284, 1826–1828 (1999).
 31.
Sun, G. Q., Jin, Z. & Li, L. et al. Spatial patterns of a predatorprey model with cross diffusion. Nonlinear Dyn. 69, 1631–1638 (2012).
 32.
Sun, G. Q., Wang, S. L. & Ren, Q. et al. Effects of time delay and space on herbivore dynamics: linking inducible defenses of plants to herbivore outbreak. Sci Rep 5, 11246 (2015).
 33.
Li, L., Jin, Z. & Li, J. Periodic solutions in a herbivoreplant system with time delay and spatial diffusion, Appl Math. Mod 40, 4765–4777 (2016).
 34.
Alonso, D., Bartumeus, F. & Catalan, J. Mutual interference between predators can give rise to Turing spatial patterns. Ecology 83, 28–34 (2002).
 35.
Sun, G. Q., Wang, C. H. & Wu, Z. Y. Pattern dynamics of a GiererMeinhardt model with spatial effects. Nonlinear Dyn 88, 1385–1396 (2017).
 36.
Sun, G. Q., Jin, Z. & Li, L. et al. Selforganized wave pattern in a predatorprey model. Nonlinear Dyn 60, 265–275 (2010).
 37.
Sun, G. Q. Mathematical modeling of population dynamics with Allee effect. Nonlinear Dyn 85, 1–12 (2016).
 38.
Liu, Q. X., Herman, P. M. J. & Mooij, W. M. et al. Pattern formation at multiple spatial scales drives the resilience of mussel bed ecosystems. Nature Commun 5, 5234 (2014).
 39.
Sun, G. Q., Wu, Z. Y. & Jin, Z. et al. Influence of isolation degree of spatial patterns on persistence of populations. Nonlinear Dyn 83, 811–819 (2016).
 40.
Ritchie, M. E. & Olff, H. Spatial scaling laws yield a synthetic theory of biodiversity. Nature 400, 557–560 (1999).
 41.
Yuan, G., Bao, X. & Yang, S. et al. Dynamics of pulses and spiral waves in excitable media with an anomalous diffusion. Int J Mod Phys B 30(20), 1650127 (2016).
 42.
Ma, J., Ying, H. P. & Li, Y. L. Suppression of spiral waves using intermittent local electric shock. Chinese Phys 16, 955–961 (2007).
 43.
Chen, J. X., Mao, J. W. & Hu, B. et al. Suppression of spirals and turbulence in inhomogeneous excitable media. Phys Rev E 79, 066209 (2009).
 44.
Chen, J. X., Peng, L. & Ma, J. et al. Liberation of a pinned spiral wave by a rotating electric pulse. EPL 107, 38001 (2014).
 45.
Chen, J. X., Zhang, H. & Qiao, L. Y. et al. Interaction of excitable waves emitted from two defects by pulsed electric fields. Commun Nonlinear Sci Numer Simulat 54, 202–209 (2018).
 46.
LacyHulbert, A., Metcalfe, J. C. & Hesketh, R. Biological response to electromagnetic fields. Faseb J 12, 395–420 (1998).
 47.
Carpenter, D. O. Human disease resulting from exposure to electromagnetic fields. Rev Environ Heal 28, 159–172 (2017).
 48.
MejiaArangure, J. M., FajardoGutierrez, A. & PerezSaldivar, M. L. et al. Magnetic fields and acute leukemia in children with Down Syndrome. Epidemiol 18, 158–161 (2007).
 49.
Foliart, D. E., Pollock, B. H. & Mezei, G. et al. Magnetic field exposure and longterm survival among children with leukaemia. Br J Cancer 94, 161–164 (2006).
 50.
Garcia, A. M., Sisternas, A. & Hoyos, S. P. Occupational exposure to extremely low frequency electric and magnetic fields and Alzheimer disease: a metaanalysis. Int J Epidemiol 37, 329–340 (2008).
 51.
Feychting, M., Jonsson, F. & Pedersen, N. L. et al. Occupational magnetic field exposure and neurodegenerative disease. Epidemiol 14, 413–419 (2003).
 52.
BurdakRothkamm, S., Rothkamm, K. & Folkard, M. et al. DNA and chromosomal damage in response to intermittent extremely lowfrequency magnetic fields. Mut Res 672, 82–89 (2009).
 53.
Salehi, I., Sani, K. G. & Zamani, A. Exposure of rats to extremely lowfrequency electromagnetic fields (ELFEMF) alters cytokines production, Electromag. Biol Med 32, 1–8 (2012).
 54.
Gultekin, D. H. & Moeller, L. NMR imaging of cell phone radiation absorption in brain tissue. PNAS 110, 58–63 (2013).
 55.
Clark, V. P. & Parasuraman, R. Neuroenhancement: Enhancing brain and mind in health and in disease. NeuroImage 85, 889–894 (2014).
 56.
Ren, G., Xu, Y. & Wang, C. Synchronization behavior of coupled neuron circuits composed of memristors. Nonlinear Dyn 88, 893–901 (2017).
 57.
Ma, J., Wu, F. & Wang, C. Synchronization behaviors of coupled neurons under electromagnetic radiation. Int J Mod Phys B 31(2), 1650251 (2017).
 58.
Yu, K., Wang, J. & Deng, B. et al. Synchronization of neuron population subject to steady DC electric field induced by magnetic stimulation. Cogn Neurodyn 7(3), 237–252 (2013).
 59.
Hindmarsh, J. L. & Rose, R. M. A Model of Neuronal Bursting Using Three Coupled First Order. Differential Equations, Proc Roy Soc London B 221(1222), 87–102 (1984).
 60.
Izhikevich, E. M. Which model to use for cortical spiking neurons? IEEE Tr Neural Netw15 5, 1063–1070 (2004).
 61.
Wu, F. Q., Wang, C. N. & Xu, Y. et al. Model of electrical activity in cardiac tissue under electromagnetic induction. Sci Rep 6, 28 (2016).
 62.
Ma, J., Wu, F. Q. & Hayat, T. et al. Electromagnetic induction and radiationinduced abnormality of wave propagation in excitable media. Physica A 486, 508–516 (2017).
 63.
Lv, M., Wang, C. & Ren, G. et al. Model of electrical activity in a neuron under magnetic flow effect. Nonlinear Dyn 85, 1479–1490 (2016).
 64.
Lv, M. & Ma, J. Multiple modes of electrical activities in a new neuron model under electromagnetic radiation. Neurocomputing 205, 375–381 (2016).
 65.
Xu, Y., Ying, H. P. & Jia, Y. et al. Autaptic regulation of electrical activities in neuron under electromagnetic induction. Sci Rep 7, 43452 (2017).
 66.
Wang, Y., Ma, J. & Xu, Y. et al. The electrical activity of neurons subject to electromagnetic induction and Gaussian White noise. Int J Bifurcat Chaos 27, 1750030 (2017).
 67.
Mvogo, A., Takembo, C. N. & Fouda, H. P. E. et al. Pattern formation in diffusive excitable systems under magnetic flow effects. Phys Lett A 381, 2264–2271 (2017).
 68.
Yi, G. S., Wang, J. & Deng, B. et al. Morphology controls how hippocampal CA1 pyramidal neuron responds to uniform electric fields: a biophysical modeling study. Sci Rep 7, 3210 (2017).
 69.
Ma, J. & Tang, J. A review for dynamics in neuron and neuronal network. Nonlinear Dyn 89, 1569–1578 (2017).
 70.
Ma, J. & Tang, J. A review for dynamics of collective behaviors of network of neurons. Sci China Technol Sci 58, 1007–1014 (2015).
 71.
Ma, J., Mi, L. & Zhou, P. et al. Phase synchronization between two neurons induced by coupling of electromagnetic field. Appl Math Comput 307, 321–328 (2017).
 72.
Wu, F. Q., Wang, C. N. & Jin, W. Y. et al. Dynamical responses in a new neuron model subjected to electromagnetic induction and phase noise. Physica A 469, 81–88 (2017).
 73.
Ren, G. D., Zhou, P. & Ma, J. et al. Dynamical response of electrical activities in digital neuron circuit driven by autapse. Int. J. Bifurcat. Chaos 27, 1750187 (2017).
 74.
Xu, Y., Jia, Y. & Ma, J. et al. Synchronization between neurons coupled by memristor. Chaos, Solitons & Fractals 104, 435–442 (2017).
 75.
Guo, S. L., Xu, Y. & Wang, C. et al. Collective response, synapse coupling and field coupling in neuronal network. Chaos, Solitons & Fractals 105, 120–127 (2017).
Acknowledgements
This project is supported by National Natural Science Foundation of China under Grants. 11672122, 11365014,11765011. The authors thank Ms. Lulu Lu(Central China Normal University) for her help in plotting Fig. 1.
Author information
Affiliations
Contributions
Y.X., YJ, J.M. implemented project and analyzed data; T.H. and A.A. extensively shared and discussed results; J.M. wrote the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Xu, Y., Jia, Y., Ma, J. et al. Collective responses in electrical activities of neurons under field coupling. Sci Rep 8, 1349 (2018). https://doi.org/10.1038/s41598018198581
Received:
Accepted:
Published:
Further reading

Synchronization and chimeras in a network of photosensitive FitzHugh–Nagumo neurons
Nonlinear Dynamics (2021)

Study on propagation efficiency and fidelity of subthreshold signal in feedforward hybrid neural network under electromagnetic radiation
Nonlinear Dynamics (2021)

Mathematical modeling and mechanisms of pattern formation in ecological systems: a review
Nonlinear Dynamics (2021)

Propagation characteristics of weak signal in feedforward Izhikevich neural networks
Nonlinear Dynamics (2020)

Pattern control of external electromagnetic stimulation to neuronal networks
Nonlinear Dynamics (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.