# An Exploration of the Range of Noise Intensity That Affects the Membrane Potential of Neurons

Authors: Rubin Wang; Guanzheng Wang; Jinchao Zheng
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

Neuronal activity in the human brain occurs in a complex physiologic environment, and noise from all aspects in this physiologic environment affects all aspects of nervous-system function. An essential issue of neural information processing is whether the environmental noise in a neural system can be estimated and quantified in a proper way. In this paper, we calculated the neural energy to estimate the range of critical values of thermal noise intensity that markedly affect the membrane potential and the energy waveform, in order to define such a noisy environment which neuronal activity relies on.


---

## Body

## 1. Introduction

Thermal noise in neural system is critically vital for information processing, because it has a great influence on a variety of aspects of the central nervous system [1–6]. How to evaluate the intensity of such noise has long been a focus for scientists attempting to understand network behaviors and perception [7–11]. So far, estimations of thermal noise intensity in neurons or neural networks in published papers with regard to neural models are generally subjective, lacking a definite basis that is applicable to the range of thermal noise intensity in real systems. Sometimes, authors define thermal noise themselves in order to obtain an ideal result; however, scientific definition of thermal noise and its intensity range is extremely difficult [12–18]. From the viewpoint of quantification, there have been a lot of studies on information coding and the dynamic behaviors of networks in a noisy environment, but these studies have not attempted to answer one of the most basic questions—what principles should be used to choose the thermal noise values in models? Theoretically speaking, there will be deviations between the calculated results and experimental data, and a further discussion on the effectiveness of calculated results is also needed. Taking into account all the considerations above and based on the literature [19–21], in this paper, a range of values of the different membrane potential under the thermal noise condition studied by using the energy method and their corresponding energy waveforms.We wanted to estimate the range of critical values of noise intensity that is capable of markedly changing the energy waveform. The principal idea of our research was inspired by the fact that since it is impossible to measure the thermal noise intensity at a neuronal level which is great enough to affect neuronal activity in an experiment, and hence, according to the rule of the only corresponding relationship between the membrane potential and its energy function, the range of thermal noise intensity obtained by our neuronal energy function can be estimated to be what it should be in a real neural system. Any kind of membrane potential can be obtained by adjusting the noise intensity, but there is no intrinsic relationship between such a membrane potential and the real neuronal energy. Therefore, we studied the membrane potential starting from the point of view of neuronal energy and observing the type of order of magnitude in noise intensity that can greatly affect the energy function of the membrane potential. As a result, we obtained a range of thermal noise intensity that neurons might have in an actual thermal noisy environment.It follows that a further discussion of thermal noise intensity range is possible for levels of networks. This part of the research is not only significant for the application of thermal noise intensity when modeling a neuron, but is also able to provide an adequate scientific basis for estimating the range of thermal noise intensity in networks and helping to establish neural network models.Finally, it should be emphasized that we did not consider signal to noise ratio (SNR) that is beyond the scope of this paper.


## 2. Methods

Compared with the traditional simple single neuronal model, a voltage source, a current source, and an inductor are innovatively proposed in the biophysical model presented in this paper, which is shown in Figure1. The difference in concentration between inner and external ions of neurons creates the electric voltage source that drives the ions. The current source is created by, on one hand, the chemical gradient of ions and, on the other hand, the stimuli from neighboring neurons. Moreover, in the ion channels, the flow of charged particles, such as sodium ions, potassium ions, and calcium ions, forms a loop current which is equivalent to an inductor Lm, producing self-induction. Cm is used to denote the membrane capacitance. Im denotes the total electrical current formed by the inputs of all the neighboring neurons to the Mth neuron. U denotes the voltage source. rm and r0m denote the resistance of the current source Im and voltage source U, respectively, showing the wastage produced by the actual sources. Due to the actions of the electric voltage source and the electric current source at different points in the neuron, the membrane resistance is divided into three parts r1m, r2m, and r3m.Figure 1
Physical model of theMth neuron with coupling to other neurons.As indicated in Figure1, the power Pm of the Mth neuron produced by the voltage source Um and the current source Im is obtained using the formula
(1)Pm=UmI0m+UimIm,
where
(2)Um=r0mI0m+r1mI1m+LmI˙1m,I0m=I1m-Im+Uimrm+CmI˙1m,Uim=Cmr3mU˙0m+U0m,(3)LmI˙1m+r1mIm=K1mU˙0m+K2mU0m-r2mIm.The current source is calculated from the formula(4)Im=i1m+∑j=1n[i0m(j-1)sinωm(j-1)(tj-tj-1)]+i0m(n)sinωm(n)(t-tn)tn<t<tn+1,n=1,2,…,t0=0,
where i1m is the current used to maintain resting potential, i0m denotes the electric current being produced after the single neuron receives stimuli from all other neighboring neurons connected at the subthreshold level, and ωm denotes the spiking frequency of the action potential, where
(5)K1m=Cm(r2m+r3m+r2mr3mrm),K2m=1+r2mrm.The solution of (3) is as follows:
(6)I1m=e-at(K+1Lm(∫t0tImeatdtK1mK1mU0meat+(K2m-aK1m)K1mU0meatK1mU0meat×∫t0tU0meatdt-r2m∫t0tImeatdt)),
where
(7)a=r1mLm,K=I1m(t0)-rm+r2mrmr1mU0m(t0).
When (6) satisfies the following conditions:
(8)K2m-aK1m=0,Lm=K1mK2mr1m=Cmr1m(r3m+rmr2mrm+r2m),
we can obtain the following equation:
(9)I1m=e-at(K+rmr2mrm+r2mU0meat-r2mLm∫t0tImeatdt).Inserting (9) and (2) into (1), the power of consumption Pm of mth neuron is obtained in the following:
(10)Pm=d1mU˙0m2+d2mU˙0m+d3mU˙0mU0m+d4mU0m2+d5mU0m+d6m,
where d1m, d2m, d3m, d4m, d5m, and d6m can be obtained from [20, 21].According to the above equations, we can obtain the solution of the action potentialUom:
(11)U0m=-g^1λm2-g^2e-a(t-tn)λm2-a2-1λm2+ωm2×(g^2e-a(t-tn)λ2m-a2g^3sin⁡(ωm(n)(t-tn))+g^4cos⁡⁡(ωm(n)(t-tn))×(U0m(tn)+g^1λm2+g^2λm2-a2+g^4λm2+ωm2(n)))e-λm(t-tn).To clarify our point of view, we present in a straightforward manner the action potential and its neuronal energy function represented by the corresponding power obtained by our proposed method, as shown in Figure2.Action potential and its corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, and i0m=70.7μA.
(a)
(b)As indicated in Figure2, the peak voltage of the action potential is around 25 mV and the peak power of the energy function is around 12 nW. The energy waveforms of an excitatory postsynaptic potential (EPSP) and inhibitory postsynaptic potential (IPSP) are shown in Figures 3 and 4, respectively.EPSP and its corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, and i0m=7.155μA.
(a)
(b)IPSP and its corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, and i0m=-0.7μA.
(a)
(b)In the cerebral cortex, the excitatory neurons comprise 85% of the neurons, and the remaining neurons are inhibitory [22]. At the subthreshold state, the sum of the power of the EPSP and IPSP of a single neuron is negative (PP = 0.0367 + (−0.0532) = −0.0165 nW). From the overall observation, the total consumed energy is positive (PP′ = 0.0367 × 85% −0.0532 × 15% = 0.0232 nW). As seen from the results, the neural system consumes the energy supplied from blood flow, even if the neurons act at a subthreshold state; also as seen from the calculated results, neurons consume more energy at suprathreshold activities than at subthreshold activities. It is reported that a massive number of neurons in subthreshold activities consume only 20% of the total energy, whereas a small amount of neurons consume more than 80% due to action potential firing [22]. Our results conform to this conclusion.A single neuron acts under the condition of a neural network. In other words, the interaction among the neurons makes their functional effectiveness emerge. It is in this sense that neuronal activity is performing the process of metabolism in the thermal noisy physiologic environment. To obtain the size of thermal noise intensity of the neuronal activity in the actual environment and to further obtain the noise circumstance under the condition of networks in the brain, we need to first understand the neuronal membrane potential and its corresponding energy function under the condition of no noise interference [19]. Only by obtaining the neuronal energy under an ideal condition can the physiologic influence of thermal noise on the membrane potential and its neuronal energy be better understood. Once that is obtained, it is possible to comprehend the noisy environment in a neural system and evaluate the behaviors of the network in such an environment.


## 3. Results

### 3.1. The Suprathreshold Membrane Potential and Its Neuronal Energy

Because the signal intensities of AP, EPSP, and IPSP are7.07×10-5μA, 7.155×10-6μA, and 7×10-7μA, so if we superpose thermal noise on the current i0m, when the thermal noise intensity Q is more than 10-4μA, the influence of thermal noise on membrane potential and its energy function is far greater than that of the current i0m, and when the thermal noise intensity Q is less than 10-7μA, the influence of thermal noise on membrane potential and its energy function is far weaker than that of i0m. Therefore, the values of thermal noise intensity should be set as 10-4μA, 10-5μA, 10-6μA, and 10-7μA, respectively.After adding thermal noise, the current takes the form(12)Im=i1m+∑j=1n[(i0m(j-1)+Q)sinωm(j-1)(tj-tj-1)]+[i0m(n)+Q]sinωm(n)(t-tn),
where Q denotes the synaptic noise as Gaussian white noise with the unit μA.The action potential together with its corresponding power plotted by the energy function calculated under the condition of different noise intensity is shown as follows.In all of Figures5, 6, 7, 8, 9, and 10, the dotted line represents the case of no noise and the solid line represents the case of noise.The action potential and its corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=1×10-4μA.
(a)
(b)The action potential and its corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=1×10-5μA.
(a)
(b)The action potential and its corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=1×10-6μA.
(a)
(b)The action potential and its corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=6×10-6μA.
(a)
(b)The action potential and its energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=6×10-6μA.
(a)
(b)The action potential and its corresponding function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=6×10-6μA.
(a)
(b)When the noise intensityQ reaches 10-4μA, the peak of the neuronal membrane potential is around 140 mV and far higher than 25 mV of the situation without noise, while the peak of the energy function is around 40 nW and far higher than the original 12 nW. Additionally, on the horizontal axis, between 2.5 ms and 3 ms, the amplitude of the energy waveform is much larger than the one in the situation without noise. Therefore, it fundamentally changes the action potential and its corresponding energy function under the ideal conditions when Q reaches 10-4μA.When the noise intensityQ reaches 10-5μA, the influence of noise on the action potential deceases, but it is still significant. Here, the membrane potential is around 40 mV, whereas the peak of the energy function is around 15 nW. As can be seen, such noise intensity still has a great influence on the membrane potential and its corresponding energy waveform.When the noise intensityQ reaches 10-6μA, there is hardly any influence on the action potential and its energy function. It is in such a noisy environment that the peaks of the membrane potential and energy waveform are exactly the same with those under the ideal conditions where the peaks of action potential and the energy function are 25 mV and 12 nW, respectively. In the presence of noise, the waveforms of the membrane potential and the energy function are nearly the same as without noise.To summarize, we cannot estimate an accurate value of noise intensity for both the suprathreshold membrane potential and the corresponding neuronal energy. In estimating the range of the critical values of noise intensity, we found that when the noise intensityQ reaches about 6×10-6μA (error ±0.5×10-6μA), the noise has great influence on the action potential and its corresponding energy function.As shown in Figures8–10, when the noise intensity Q is 6×10-6μA, the fluctuation of the membrane potential and the energy function turns out to be very large. The peak of the action potential fluctuates between 10 mV and 40 mV, while that of the energy function fluctuates from 10 nW to 15 nW. Such noise intensity has great influence on the membrane potential and the energy function, and as a result, the critical value of the noise intensity Q can be assessed at 6×10-6μA (error ±0.5×10-6μA).


### 3.2. Subthreshold Membrane Potential and Its Neuronal Energy

According to the calculated results of Wang et al. [21], we obtained the EPSP and the corresponding energy waveform as well as the IPSP and the corresponding energy waveform by adding noise (Figures 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, and 25).The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=1×10-4μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=1×10-5μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=1×10-6μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=2×10-6μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=3×10-6μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=3×10-6μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=3×10-6μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=1×10-7μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=1×10-4μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=1×10-5μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=1×10-6μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=1×10-7μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=2×10-7μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=3×10-7μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=4×10-7μA.
(a)
(b)When the level of noise intensityQ reached 1×10-4μA, it had a great effect on the EPSP and the corresponding energy waveform. Without the effect of noise, we found that the EPSP curve increased first and then decreases with a maximum of −68 mV and a minimum of −70 mV. The energy function fell first and then rose towards stability, fluctuating in the range of −0.1 to 0.1 nW. After adding noise, we found that the EPSP curve fell first and then rose towards a minimum of −95 mV, while the energy function increased first and then decreased towards stability, fluctuating in the range of −10 to −80 nW.When the level of noise intensityQ reached 1×10-5μA, it also had a great effect on the EPSP and the corresponding energy waveform. Under an ideal condition without noise, we know that the maximum of the EPSP is about −68 mV and the peak of the energy function is about 0.1 nW. After adding noise, the peak of the EPSP hardly changed, but the oscillatory curve changed sharply in the period from 2.5 ms to 3 ms. The oscillatory amplitude of the energy function with noise was more intense than that without noise, varying between 1.5 ms and 2.5 ms.When the noise intensityQ reached 1×10-6μA, it had little effect on the EPSP and the corresponding energy waveform. This figure showing the EPSP and energy function with noise is almost the same with one without noise. We could not estimate the range of the noise intensity, which affects the EPSP’s corresponding energy function according to the results shown in Figures 11–13. As a solution, we addressed the problem by adjusting the coefficient of the noise level.When the level of noise intensityQ reached 2×10-6μA, it did not have an obvious effect on the EPSP and the corresponding energy waveform. This figure shows the EPSP and energy function with noise is almost the same as one without noise. But the energy function with noise is slightly different from that without noise ranging between 1.5 ms and 2.5 ms.When the noise intensityQ reached 3×10-6μA, it had a great effect on the EPSP and the corresponding energy waveform. Taking the noise into account, we found that the EPSP fluctuated in a small interval from −67 mV to −69 mV, while in the period from 1.5 ms to 2.5 ms, the energy function with noise had greater fluctuations than that without noise.When the level of noiseQ reached 1×10-7μA, it had a little effect on the EPSP and the corresponding energy waveform. The figures of the EPSP and energy function with noise are almost the same as those without noise.In summary, when the noise intensityQ was about 3×10-6μA (error ±0.5×10-6μA), the noise begins to have a significant influence on EPSP and its corresponding energy waveform. As a result, the critical value of the noise was assessed at 3×10-6μA (error ±0.5×10-6μA).When the noise intensityQ reached 1×10-4μA, it had a great effect on the IPSP and the corresponding energy waveform. Without the effect of noise, we found that the IPSP curve fell first and then rose and the energy function curve was almost in the stable state of −0.04 nW. However, after adding noise, we found that the IPSP curve rose first and then fell, while most of the energy function curve was above 0 and there were huge fluctuations of the energy function between 1.5 ms and 2.5 ms.When the noise intensityQ reached 1×10-5μA, it had a critical effect on the IPSP and the corresponding energy waveform. Without the effect of noise, we observed that the IPSP curve fell initially and then rose and the energy function curve was almost in a stable state of −0.04 nW. However, after adding noise, the IPSP curve rose first and then fell and there were huge fluctuations of the energy function between 2.5 ms and 3 ms.When the noise intensityQ reached 1×10-6μA, it had a relatively great effect on the IPSP and the corresponding energy waveform. Without the effect of noise, the IPSP curve fell initially and then rose, while the energy function curve was almost in a stable state of −0.04 nW. After adding noise, the IPSP curve rose first and then fell and there were also huge fluctuations of the energy function from 2.5 ms to 3 ms.When the noise intensityQ reached 1×10-7μA, it had a little effect on the IPSP and the corresponding energy waveform. The figures of the IPSP and energy function with noise are almost the same as those without noise.We still could not estimate the range of the noise that affects the energy function of the IPSP according to the results shown in Figures19–22. So as a solution, we addressed it by adjusting the coefficient of the level of noise intensity.When the noise intensityQ reached 2×10-7μA, it had a little effect on the IPSP and the corresponding energy waveform. The figures of the IPSP and energy function with noise are almost the same as those without noise.When the level of noise intensityQ reached 3×10-7μA, it had a little effect on the IPSP and the corresponding energy waveform. The figures of the IPSP and energy function with noise are almost the same as those without noise.When the level of noise intensityQ reached about 4×10-7μA (error ±0.5×10-7μA), the noise began to significantly influence the IPSP and its corresponding energy waveform. Without the effect of noise, the IPSP curve initially fell and then rose; however, with the effect of noise, the IPSP curve did not monotonically decrease in the time interval from 0 ms to 1.5 ms, while the energy function began to appear to fluctuate within the interval from 2.5 ms to 3 ms.In summary, when the level of noise intensityQ reached about 4×10-7μA (error ±0.5×10-7μA), the noise began to have significant influence on the IPSP and its corresponding energy waveform. As a consequence, its critical value was assessed at 4×10-7μA (error ±0.5×10-7μA).


## 3.1. The Suprathreshold Membrane Potential and Its Neuronal Energy

Because the signal intensities of AP, EPSP, and IPSP are7.07×10-5μA, 7.155×10-6μA, and 7×10-7μA, so if we superpose thermal noise on the current i0m, when the thermal noise intensity Q is more than 10-4μA, the influence of thermal noise on membrane potential and its energy function is far greater than that of the current i0m, and when the thermal noise intensity Q is less than 10-7μA, the influence of thermal noise on membrane potential and its energy function is far weaker than that of i0m. Therefore, the values of thermal noise intensity should be set as 10-4μA, 10-5μA, 10-6μA, and 10-7μA, respectively.After adding thermal noise, the current takes the form(12)Im=i1m+∑j=1n[(i0m(j-1)+Q)sinωm(j-1)(tj-tj-1)]+[i0m(n)+Q]sinωm(n)(t-tn),
where Q denotes the synaptic noise as Gaussian white noise with the unit μA.The action potential together with its corresponding power plotted by the energy function calculated under the condition of different noise intensity is shown as follows.In all of Figures5, 6, 7, 8, 9, and 10, the dotted line represents the case of no noise and the solid line represents the case of noise.The action potential and its corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=1×10-4μA.
(a)
(b)The action potential and its corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=1×10-5μA.
(a)
(b)The action potential and its corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=1×10-6μA.
(a)
(b)The action potential and its corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=6×10-6μA.
(a)
(b)The action potential and its energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=6×10-6μA.
(a)
(b)The action potential and its corresponding function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=70.7μA, and Q=6×10-6μA.
(a)
(b)When the noise intensityQ reaches 10-4μA, the peak of the neuronal membrane potential is around 140 mV and far higher than 25 mV of the situation without noise, while the peak of the energy function is around 40 nW and far higher than the original 12 nW. Additionally, on the horizontal axis, between 2.5 ms and 3 ms, the amplitude of the energy waveform is much larger than the one in the situation without noise. Therefore, it fundamentally changes the action potential and its corresponding energy function under the ideal conditions when Q reaches 10-4μA.When the noise intensityQ reaches 10-5μA, the influence of noise on the action potential deceases, but it is still significant. Here, the membrane potential is around 40 mV, whereas the peak of the energy function is around 15 nW. As can be seen, such noise intensity still has a great influence on the membrane potential and its corresponding energy waveform.When the noise intensityQ reaches 10-6μA, there is hardly any influence on the action potential and its energy function. It is in such a noisy environment that the peaks of the membrane potential and energy waveform are exactly the same with those under the ideal conditions where the peaks of action potential and the energy function are 25 mV and 12 nW, respectively. In the presence of noise, the waveforms of the membrane potential and the energy function are nearly the same as without noise.To summarize, we cannot estimate an accurate value of noise intensity for both the suprathreshold membrane potential and the corresponding neuronal energy. In estimating the range of the critical values of noise intensity, we found that when the noise intensityQ reaches about 6×10-6μA (error ±0.5×10-6μA), the noise has great influence on the action potential and its corresponding energy function.As shown in Figures8–10, when the noise intensity Q is 6×10-6μA, the fluctuation of the membrane potential and the energy function turns out to be very large. The peak of the action potential fluctuates between 10 mV and 40 mV, while that of the energy function fluctuates from 10 nW to 15 nW. Such noise intensity has great influence on the membrane potential and the energy function, and as a result, the critical value of the noise intensity Q can be assessed at 6×10-6μA (error ±0.5×10-6μA).


## 3.2. Subthreshold Membrane Potential and Its Neuronal Energy

According to the calculated results of Wang et al. [21], we obtained the EPSP and the corresponding energy waveform as well as the IPSP and the corresponding energy waveform by adding noise (Figures 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, and 25).The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=1×10-4μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=1×10-5μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=1×10-6μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=2×10-6μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=3×10-6μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=3×10-6μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=3×10-6μA.
(a)
(b)The EPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=7.155μA, and Q=1×10-7μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=1×10-4μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=1×10-5μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=1×10-6μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=1×10-7μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=2×10-7μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=3×10-7μA.
(a)
(b)The IPSP and the corresponding energy function.r0m=0.0001Ω, r1m=0.1Ω, r2m=1000Ω, r3m=0.1Ω, rm=1000Ω, Cm=1μF, Lm=50μH, i0m=-0.7μA, and Q=4×10-7μA.
(a)
(b)When the level of noise intensityQ reached 1×10-4μA, it had a great effect on the EPSP and the corresponding energy waveform. Without the effect of noise, we found that the EPSP curve increased first and then decreases with a maximum of −68 mV and a minimum of −70 mV. The energy function fell first and then rose towards stability, fluctuating in the range of −0.1 to 0.1 nW. After adding noise, we found that the EPSP curve fell first and then rose towards a minimum of −95 mV, while the energy function increased first and then decreased towards stability, fluctuating in the range of −10 to −80 nW.When the level of noise intensityQ reached 1×10-5μA, it also had a great effect on the EPSP and the corresponding energy waveform. Under an ideal condition without noise, we know that the maximum of the EPSP is about −68 mV and the peak of the energy function is about 0.1 nW. After adding noise, the peak of the EPSP hardly changed, but the oscillatory curve changed sharply in the period from 2.5 ms to 3 ms. The oscillatory amplitude of the energy function with noise was more intense than that without noise, varying between 1.5 ms and 2.5 ms.When the noise intensityQ reached 1×10-6μA, it had little effect on the EPSP and the corresponding energy waveform. This figure showing the EPSP and energy function with noise is almost the same with one without noise. We could not estimate the range of the noise intensity, which affects the EPSP’s corresponding energy function according to the results shown in Figures 11–13. As a solution, we addressed the problem by adjusting the coefficient of the noise level.When the level of noise intensityQ reached 2×10-6μA, it did not have an obvious effect on the EPSP and the corresponding energy waveform. This figure shows the EPSP and energy function with noise is almost the same as one without noise. But the energy function with noise is slightly different from that without noise ranging between 1.5 ms and 2.5 ms.When the noise intensityQ reached 3×10-6μA, it had a great effect on the EPSP and the corresponding energy waveform. Taking the noise into account, we found that the EPSP fluctuated in a small interval from −67 mV to −69 mV, while in the period from 1.5 ms to 2.5 ms, the energy function with noise had greater fluctuations than that without noise.When the level of noiseQ reached 1×10-7μA, it had a little effect on the EPSP and the corresponding energy waveform. The figures of the EPSP and energy function with noise are almost the same as those without noise.In summary, when the noise intensityQ was about 3×10-6μA (error ±0.5×10-6μA), the noise begins to have a significant influence on EPSP and its corresponding energy waveform. As a result, the critical value of the noise was assessed at 3×10-6μA (error ±0.5×10-6μA).When the noise intensityQ reached 1×10-4μA, it had a great effect on the IPSP and the corresponding energy waveform. Without the effect of noise, we found that the IPSP curve fell first and then rose and the energy function curve was almost in the stable state of −0.04 nW. However, after adding noise, we found that the IPSP curve rose first and then fell, while most of the energy function curve was above 0 and there were huge fluctuations of the energy function between 1.5 ms and 2.5 ms.When the noise intensityQ reached 1×10-5μA, it had a critical effect on the IPSP and the corresponding energy waveform. Without the effect of noise, we observed that the IPSP curve fell initially and then rose and the energy function curve was almost in a stable state of −0.04 nW. However, after adding noise, the IPSP curve rose first and then fell and there were huge fluctuations of the energy function between 2.5 ms and 3 ms.When the noise intensityQ reached 1×10-6μA, it had a relatively great effect on the IPSP and the corresponding energy waveform. Without the effect of noise, the IPSP curve fell initially and then rose, while the energy function curve was almost in a stable state of −0.04 nW. After adding noise, the IPSP curve rose first and then fell and there were also huge fluctuations of the energy function from 2.5 ms to 3 ms.When the noise intensityQ reached 1×10-7μA, it had a little effect on the IPSP and the corresponding energy waveform. The figures of the IPSP and energy function with noise are almost the same as those without noise.We still could not estimate the range of the noise that affects the energy function of the IPSP according to the results shown in Figures19–22. So as a solution, we addressed it by adjusting the coefficient of the level of noise intensity.When the noise intensityQ reached 2×10-7μA, it had a little effect on the IPSP and the corresponding energy waveform. The figures of the IPSP and energy function with noise are almost the same as those without noise.When the level of noise intensityQ reached 3×10-7μA, it had a little effect on the IPSP and the corresponding energy waveform. The figures of the IPSP and energy function with noise are almost the same as those without noise.When the level of noise intensityQ reached about 4×10-7μA (error ±0.5×10-7μA), the noise began to significantly influence the IPSP and its corresponding energy waveform. Without the effect of noise, the IPSP curve initially fell and then rose; however, with the effect of noise, the IPSP curve did not monotonically decrease in the time interval from 0 ms to 1.5 ms, while the energy function began to appear to fluctuate within the interval from 2.5 ms to 3 ms.In summary, when the level of noise intensityQ reached about 4×10-7μA (error ±0.5×10-7μA), the noise began to have significant influence on the IPSP and its corresponding energy waveform. As a consequence, its critical value was assessed at 4×10-7μA (error ±0.5×10-7μA).


## 4. Discussion

In this study, we obtained the action potential, EPSP, IPSP, and their corresponding energy waveforms with thermal noise added to the current. By changing the thermal noise intensity, we found the estimated range of the critical value of thermal noise that can significantly influence the neuronal membrane potential and the corresponding energy waveform. When the thermal noise intensityQ reached about 6×10-6μA (error ±0.5×10-6μA), it had a great influence on the action potential and its corresponding energy waveform. When the thermal noise intensity Q reached about 3×10-6μA (error ±0.5×10-6μA), it had a great influence on the EPSP and its corresponding energy waveform. When the thermal noise intensityQ reached about 4×10-7μA (error ±0.5×10-7μA), it had a great influence on the IPSP and its corresponding energy waveform. Based on these results, we conclude that noise has less effect on the action potential and its energy waveform than on the EPSP, IPSP, and their energy waveforms, and the thermal noise has less effect on the EPSP and its energy waveform than the IPSP and its waveform. By numerical calculation, we have also obtained the estimated range of the critical value of the noise that has a great effect on the action potential, EPSP, IPSP, and their corresponding energy waveforms. As a result, when we again establish a neural dynamic model in the future, such an estimation of noise intensity will enable us to avoid being subjective in considering the influence of noise on neural firing, so that the calculated results and experimental data will be more consistent. Therefore, the estimation of the thermal noise intensity provides a scientific basis for research on neural information coding in a noisy environment.The physical model which we used in this paper was proposed by Wang et al. [20, 21] from the viewpoint of energy. Compared with other models, the biggest advantage of this model is that it simplifies the calculation of activity of neuronal firing. Due to the energy coding model’s ability to describe mechanisms of brain information processing in biophysical terms, this idea and method can be applied to future research of information encoding in neural ensembles as well as cognitive functioning.



---

*Source: 801642-2014-02-12.xml*