
Citation: Kanishka Jha, Ravinder Kataria, Jagesvar Verma, Swastik Pradhan. Potential biodegradable matrices and fiber treatment for green composites: A review[J]. AIMS Materials Science, 2019, 6(1): 119-138. doi: 10.3934/matersci.2019.1.119
[1] | V. Volpert, B. Xu, A. Tchechmedjiev, S. Harispe, A. Aksenov, Q. Mesnildrey, A. Beuter . Characterization of spatiotemporal dynamics in EEG data during picture naming with optical flow patterns. Mathematical Biosciences and Engineering, 2023, 20(6): 11429-11463. doi: 10.3934/mbe.2023507 |
[2] | Sridevi Sriram, Hayder Natiq, Karthikeyan Rajagopal, Ondrej Krejcar, Hamidreza Namazi . Dynamics of a two-layer neuronal network with asymmetry in coupling. Mathematical Biosciences and Engineering, 2023, 20(2): 2908-2919. doi: 10.3934/mbe.2023137 |
[3] | Diego Fasoli, Stefano Panzeri . Mathematical studies of the dynamics of finite-size binary neural networks: A review of recent progress. Mathematical Biosciences and Engineering, 2019, 16(6): 8025-8059. doi: 10.3934/mbe.2019404 |
[4] | Anna Cattani . FitzHugh-Nagumo equations with generalized diffusive coupling. Mathematical Biosciences and Engineering, 2014, 11(2): 203-215. doi: 10.3934/mbe.2014.11.203 |
[5] | Zhenglong Tang, Chao Chen . Spatio-temporal information enhance graph convolutional networks: A deep learning framework for ride-hailing demand prediction. Mathematical Biosciences and Engineering, 2024, 21(2): 2542-2567. doi: 10.3934/mbe.2024112 |
[6] | Xiaowen Jia, Jingxia Chen, Kexin Liu, Qian Wang, Jialing He . Multimodal depression detection based on an attention graph convolution and transformer. Mathematical Biosciences and Engineering, 2025, 22(3): 652-676. doi: 10.3934/mbe.2025024 |
[7] | Mahtab Mehrabbeik, Fatemeh Parastesh, Janarthanan Ramadoss, Karthikeyan Rajagopal, Hamidreza Namazi, Sajad Jafari . Synchronization and chimera states in the network of electrochemically coupled memristive Rulkov neuron maps. Mathematical Biosciences and Engineering, 2021, 18(6): 9394-9409. doi: 10.3934/mbe.2021462 |
[8] | Felicia Maria G. Magpantay, Xingfu Zou . Wave fronts in neuronal fields with nonlocal post-synaptic axonal connections and delayed nonlocal feedback connections. Mathematical Biosciences and Engineering, 2010, 7(2): 421-442. doi: 10.3934/mbe.2010.7.421 |
[9] | P. E. Greenwood, L. M. Ward . Rapidly forming, slowly evolving, spatial patterns from quasi-cycle Mexican Hat coupling. Mathematical Biosciences and Engineering, 2019, 16(6): 6769-6793. doi: 10.3934/mbe.2019338 |
[10] | Xin-You Meng, Tao Zhang . The impact of media on the spatiotemporal pattern dynamics of a reaction-diffusion epidemic model. Mathematical Biosciences and Engineering, 2020, 17(4): 4034-4047. doi: 10.3934/mbe.2020223 |
As an important branch in the field of nonlinear dynamics, the research scope of pattern dynamics involves biology, physics, mathematics, chemistry, medicine, astronomy and other disciplines, and the research in various interdiscipline has become a hot spot for researchers [1]. Pattern formation is a kind of non-uniform macro structure with certain regularity in time or space, which extensively exists in nature, such as pattern formed by self-organization in organic polymers, crystal structure in inorganic chemistry, stripes on the surface of animal fur, streak clouds in the sky, etc [2]. From the thermodynamic point of view, the first two types of patterns are those existing in the thermodynamic equilibrium state, and the last two examples are those generated when leaving the thermodynamic equilibrium state [3]. Pattern dynamics is an important branch in the field of nonlinear science [4]. Its research purpose is to explore the basic laws of pattern formation and evolution, which are common among various systems in the objective world and have universal guiding significance [5]. In recent years, as a group of interconnected basic units, complex dynamic networks have attracted more and more attention in many fields such as social science, biology, mathematics and engineering science [6]. In complex dynamic networks, one of the interesting and remarkable phenomena is to show collective behavior [7]. Collective behavior is an important attribute of complex systems, which means that the whole is greater than the sum of parts [8]. Because complex networks are natural spatiotemporal systems, patterns can be observed in these systems [9]. Spatial patterns, such as target patterns, spiral waves and vortex waves, are examples of patterns that occur in complex systems, such as reaction-diffusion systems and neural networks [10].
Spiral wave is a kind of important non-equilibrium pattern, which reflects the macroscopic structure of nonlinear system with some special laws in time or space, and widely exists in excitable, oscillatory and bistable systems [11]. It's mentioned in many literatures that spiral wave is a traveling wave rotating outward from the spiral center in a two-dimensional excitable medium [12]. Generally speaking, the formation of spiral wave is either the free end generated in the propagation wave or the wave collision generated in the inhomogeneous medium [13]. Many interesting works have been carried out, most of which have been proved to be effective in removing spiral waves and preventing spiral waves from breaking [14]. It's believed that the spiral wave in the heart tissue is harmful, so many schemes have been proposed to inhibit the spiral wave and prevent the spiral wave from breaking [15]. Neural network is a complex dynamic network, which can show many active states of spatial structures, such as spatiotemporal chaos, stochastic resonance, synchronization, chimeric state and spiral wave [16]. Network structure, neural model, even considering disturbance and noise, which are inherent properties of neural network, will affect the formation of patterns [17]. Studying these phenomena can bring new insights into the function of neurons [18,19]. Among the tools for studying neural networks, wave propagation mechanism is one of the most effective mechanisms [20]. Spiral wave is a kind of collective behavior, which exists widely in nature [21]. In fact, spiral wave is a special propagation of nonlinear wave [22]. It rotates around a center (called seed), which determines the wave dynamics [23]. It is very important to study the dynamics of spiral waves, because they have been observed in neocortex and arrhythmia of mammals [24]. It has been proved that both atrial fibrillation and ventricular fibrillation are caused by spiral waves [25]. The spiral seed in the heart fiber rotates more frequently than the natural frequency of the heart, making the heartbeat irregular [26]. Therefore, it can cause fibrous fibrillation [27]. Therefore, modeling and identifying the structures formed by spiral waves can help design methods to eliminate spiral waves and control fibrillation [28].
The combination of chemical, physical, electrical and structural characteristics of neurons makes them highly complex dynamic units [10]. Zhao et al. introduced some recent developments in the dynamic behavior of complex networks (CNs) and complex networks with multi-weights (CNMWs) under various control methods, and several output synchronization criteria for multiple output coupled complex networks (MOCCNs) were formulated by using the Lyapunov functional method and inequality techniques. [29,30]. One of the best tools to deal with this advanced complexity is to study neural spatiotemporal patterns [31]. The electrophysiological and structural characteristics of neurons in neural networks lead to complex characterization [32,33]. These performances are the driving force behind our biological behavior as human beings [34]. Technically, the firing rhythm of neurons determines the function of the brain [35]. These rhythms mainly include spikes and bursts, or a combination of these two modes [36]. In other words, they help us better understand their functions, although there is little detailed information about their chemical, electrical, morphological and structural properties [37]. In this field, one of the most important models is the spiral model. This pattern can be observed in many chemical, biological, physical and ecological systems [38]. Spirals are unique because they are self-organizing and self-sustaining [10,39]. They can play a regulatory role in a system and change its dynamics [40]. For example, there is experimental evidence that spiral waves are crucial to some ongoing cortical activities, because they act as rhythmic regulators in neural populations [41]. Sleep disorders, seizures and attention deficit hyperactivity disorder are just a few examples [42]. Spiral waves can also cause arrhythmias. They are considered to be the main cause of reentry wave front [43]. Reentry is one of the most prominent types of arrhythmia, which can lead to sudden death [44]. From the perspective of experimental research, spiral waves widely exist in nature, such as the oxidation of carbon monoxide on platinum [45]. Biological experiments found that spiral waves also exist in the cerebral cortex [46]. The appearance of spiral wave is also observed in human cardiac muscle [47]. Using a simple active medium model, Kuklik et al. studied the influence of spatial spreading inhomogeneity of transverse element coupling on spiral wave trajectory [48]. They also used the FitzHugh-Nagumo model of an excitable medium to investigate the effect of the random perturbation of cell coupling on the stability of a spiral wave in 2010 [49]. They believed that electrical cardioversion could lead to one of three outcomes, such as immediate termination of arrhythmic activity, delayed termination or unsuccessful termination, and they propose a model of atrial fibrillation as a coexistence of several spiral waves fixed in an active medium with inhomogeneity [50]. In Kumar and Amita Das's work, molecular dynamics simulation was used to prove the excitation of two-dimensional dusty plasma at the particle level [51]. Kwon et al. clearly demonstrated that the new mechanism can create a period-2 helix by computer simulation of a simple mathematical model describing the dynamics of the spiral wavefront [52]. Lacitignola et al. tested some findings by numerical approximations of the complete model and found interesting scenarios that led to spiral fracture to obtain appropriate changes in system parameters [53]. Studies have shown that researchers studying the firing activity of neurons in the cerebral cortex have also observed patches with spiral waves that are closely related to the information transmission between neurons in the brain's neural network [54]. In experimental studies, researchers have found that cardiac patients can also discharge myocardial tissue cells in spiral waves, and more seriously, spiral wave rupture may cause heart fibrillation, resulting in sudden cardiac death and causing serious consequences [55]. In general, the study of spiral waves has very important practical significance [8].
Little problem has been reported about the spatiotemporal dynamics in neural networks constructed by Izhikevich neurons derived from the modeling of cortical neurons at present [56]. At the same time, considering that Izhikevich neurons have many different types of discharge patterns, people basically only consider them as excitatory or inhibitory neurons to build neural networks when they are used for analysis, but rarely consider the formation and rupture of spiral waves in the neural networks [57]. According to the existing research, the neural network of the cerebral cortex has a 5-layer structure, so the bi-layer network is the most basic component. Neurons will inevitably connect with other neurons around them, which is actually the channel for information exchange and connection. Under the influence of these complex factors, neural networks will show different synchronization properties and spatiotemporal patterns. Moreover, it is unclear whether the neural networks constructed by different types of Izhikevich neurons can induce spiral waves [58]. In this paper, the matrix neural networks constructed by several different types of Izhikevich neurons under the random boundary conditions are discussed, and the influence of the coupling strength between neurons on spatiotemporal dynamics are investigated. A square neural network driven by random boundary is constructed firstly, and then multiple connection regions are set on the network to connect with the second layer neural network. The first layer generates spiral waves, which are transmitted from multiple regions to the second layer, and then the formation and rupture of spiral waves are observed from the second layer. In order to have a clearer understanding of cortical neural networks, the synchronization of the second layer network is also studied by changing the coupling strength between neurons as well as the inter-layer connection strength, and the change of synchronization factor with the coupling strength between neurons in the second layer also shows very interesting results. The arrangement of the paper is as follows: In Section 2, the Izhikevich neuronal model is introduced and the matrix network is constructed; In Section 3, the numerical simulation results are analyzed; Section 4 summarizes the important research conclusions of the research.
To understand how the brain works, people needs to combine experimental studies of the animal and human nervous systems with numerical simulations of large-scale brain models [59]. Eugene M. Izhikevich proposed a neuronal model in 2003 which is computationally simple, but capable of producing the rich firing patterns exhibited by real biological neurons [60]. The Izhikevich neuronal model is biologically as plausible as the Hodgkin-Huxley neural model, and is as computationally efficient as the integrate-and-fire neural model [61]. The reason why we choose Izhikevich neuronal model from various neuronal models is that the main modeling object of this model is cortical and thalamic neurons, which can reproduce all their known neuronal firing behaviors, and has simple structure, rich physiological significance, and high computing efficiency. And more importantly, it is applicable to network simulation. The Izhikevich neuronal model driven by external stimulation currents can be represented as follows
{dvdt=0.04v2+5v+140−u+I,dudt=a(bv−u). | (1) |
where v represents the membrane potential of the cerebral cortical neuron, and u represents the recovery variable, both of them are dimensionless variables [62]. Constants a and b are used to control different types of neurons, I is considered as external stimulus current [63].
When the value of neuronal membrane potential is greater than the peak value, that is, if v > 30 mV, the membrane potential is reset in the following way
{v→c,u→u+d. | (2) |
where c and d are constants.
Different setting of the constants a, b, c and d yields several typical models of the Izhikevich neuron with different firing types, for example, regular spiking (RS), fast spiking (FS), Chattering (CH) and intrinsically bursting (IB) [64]. The specific values and corresponding discharge types are given in the table as below.
In order to study the collective properties and spatiotemporal patterns of neural network, a matrix neural network which contains 200 × 200 nodes is constructed in the first step, and all neurons are evenly placed in each node [65]. A certain neuron is connected to other neurons at four locations, upper, lower, left and right, with connection strength D, and the schematic diagram of the connections between neurons is plotted in Figure 1. For the boundary of the matrix neural network, the no-flow boundary condition is considered [66]. The no-flow boundary condition considers that the current value inside the boundary is equal to that outside the boundary, that is, the current value outside the boundary is the same as the value setting inside the boundary, thus the current flowing into the network is zero, so it is called "no-flow boundary". Next, the second layer of the matrix neural network is constructed again in this manner and it is considered connecting each other with channels from multiple regions [67]. Izhikevich neural model can be used to represent the dynamical equations of two-layer network, each of which is connected to the nearest neighbor type in a two-dimensional matrix [68]. The collective behaviors of two-layer network can be represented by
{dv1ijdt=0.04v21ij+5v1ij+140−u1ij+I1ext+D1(v1i−1,j+v1i+1,j+v1i,j−1+v1i,j+1−4v1i,j),du1ijdt=a(bv1ij−u1ij),dv2ijdt=0.04v22ij+5v2ij+140−u2ij+I2ext+D2(v2i−1,j+v2i+1,j+v2i,j−1+v2i,j+1−4v2i,j) +k(v1ij−v2ij)δiαδjβ,du2ijdt=a(bv2ij−u2ij). | (3) |
where subscripts 1 and 2 represent the first and second layer network, and the subscript (ij) denotes the location of the node in the same layer. It can be seen from the structure of the network that each neuron is connected to four neurons, which means that the degree of the network is 4, and the total number of neurons in the network is 40,000, which indicates that the scale of the neural network is defined as 40,000.
In addition, D1 and D2 in Eq (3) are used to represent the coupling intensity of nearest adjacent nodes in the bi-layer network, and each layer is placed in a two-dimensional matrix, as shown in Figure 2. The intensity of the channel between the two layers is expressed by k. δiα = 1 for α = i and δjβ = 1 for β = j [69]. Otherwise, δiα = 0 and δjβ = 0. i, j, α and β are integers. For bi-layer neural networks, multiple connection regions are opened at specified locations for generating information exchange between the bi-layer networks, respectively [70]. For example, in Figure 2(a) the connection region is opened between nodes 99 to 102. Figure 2(b) expresses the case that the two layers connect at two local areas (99 ≤ α ≤ 102, 65 ≤ β ≤ 68), (99 ≤ α ≤ 102,131 ≤ β ≤ 134). For the three local coupling areas (65 ≤ α, β ≤ 68), (65 ≤ α ≤ 68,131 ≤ β ≤ 134), (131 ≤ α ≤ 134, 99 ≤ β ≤ 102) are investigated in Figure 2(c) and the four local coupling areas (65 ≤ α, β ≤ 68), (65 ≤ α ≤ 68,131 ≤ β ≤ 134), (131 ≤ α ≤ 134, 65 ≤ β ≤ 68), (131 ≤ α, β ≤ 134) are displayed in Figure 2(d).
To investigate the statistical features of the collective dynamics in the neuronal network more systematically, the synchronization factor R of the neural network is calculated by using the mean-field theory [71]. Synchronization factor R can be calculated as follows
{F=1N2N∑j=1N∑i=1vij;R=⟨F2⟩−⟨F⟩21N2N∑j=1N∑i=1(⟨v2ij⟩−⟨vij⟩2). | (4) |
where vij denotes the membrane potential of each layer nodes (i, j) neuron and it could be calculated from Eq (1). N represents the location of the neuronal node, N2 represents the number of nodes in the network, and symbols < > means that the variables are averaged over time [72]. In particular, if the value of R approaches to 1, it means that the firing behavior of all neurons exhibits a fully synchronized state; when the value of R approaches to 0, it indicates that the neuronal system is in a incompletely synchronized state [73] Previous studies have shown that smaller values of synchronization will support ordered spatial patterns, while larger values of synchronization can develop homogeneous states [74]. Appropriate values of the synchronization factor can generate graceful spatial waves in the network [75].
The collective dynamics are calculated by using Euler algorithm with a time step of 0.02 when exploring the spatiotemporal properties of neural networks [76]. In each layer, the 200 × 200 (Ni, Nj) neuron nodes are uniformly embedded into a two-dimensional square array, with a near-neighbor coupling action between the neurons and considering the no-flow boundary conditions [77]. If there is no special instruction, each neuron in the first layer will be applied an direct current as external force, and the direct current signal applied to each Izhikevich neuron is Iext = 10 [78]. On the boundaries of the first layer of the matrix neural network, all random initial values will be generated by random functions as v0 = 0.8ξln(i) – 0.2ξln(j) – 3, u0 = –0.8ξln(i) + 0.2ξln(j) – 5, where ξ represents a random number between 0 and 1 [79,80,81,82]. And initial values of all other neuron will be selected with the same values as (v0, u0) = (0, 0).
The time-series of the neuronal membrane potential are calculated using Euler algorithm method for a certain time period according to the parameters given in Table 1, and the results are plotted in Figure 3. The results in Figure 3(a) illustrate that the membrane potential of Izhikevich neuron in the CH firing mode exhibits a series of cluster firing forms with periodic properties, and each cluster spike exhibits a modal morphology. For the FS firing pattern, as shown in Figure 3(b), the modelling data is derived from inhibitory cortical neurons, where the membrane potential exhibits a periodic firing sequence with an extremely high frequency. The IB discharge pattern is first showing patterned discharge clusters followed by a series of repeated discharge spikes, as shown in Figure 3(c). The neurons in RS states are the most typical neurons in the cortex, and their firing frequency is not too fast because the limitation of their parameters determines that the neurons in RS states have the characteristics of spike frequency adaptation, as shown in Figure 3(d). Unlike neurons in FS states, neurons in RS, CH, and IB states are all used to imitate excitatory neurons, and previous studies have shown that the ratio of excitatory and inhibitory neurons is generally 4:1 in cortical neural network.
Types | a | b | c | d |
RS | 0.02 | 0.2 | –65 | 8 |
FS | 0.1 | 0.2 | –65 | 2 |
CH | 0.02 | 0.2 | –50 | 2 |
IB | 0.02 | 0.2 | –55 | 4 |
In order to understand the specific properties of the four different discharge modes of Izhikevich neuron more accurately, the bifurcation of the membrane potential is drawn in Figure 4. It is known that bifurcation often appears in the mathematical research of dynamic systems, which refers to small and continuous changes in system parameters, resulting in sudden changes in the nature of the system. Here, how the change of external stimulation current affect the spike interval of membrane potential is analyzed, and the analysis results obtained are represented by bifurcation of ISI (inter-spike interval). Figure 4(a), (b) make it clear that the peak interval of RS and FS neuronal membrane potential gradually decreases with increasing stimulus current intensity I, and the ISI curve shows a tendency to decrease monotonically with the increasing of stimulation current intensity from an overall perspective. These results mean that the RS and FS neurons are in spiking state, and the time interval between discharge spikes keeps getting smaller. Above conclusions can also be drawn from the time series of membrane potential in Figure 3. However, a obvious difference between the two case is that the peak interval of neuron in FS discharge state is smaller overall than the ISI of neuron in RS discharge state. Figure 4(c), (d) demonstrate that neurons in CH and IB discharge states exhibit more abundant firing properties when the external stimulation current intensity is relatively small. For example, for neuron in CH discharge state, the discharge mode of Izhikevich neuron undergoes the transition from period-2 bursting state to bursting state when the external stimulus current intensity is greater than 3.5. If the external stimulation current intensity is less than 3.5, the discharge mode includes chaotic state, period-3 bursting state and period-2 bursting state. For neuron in IB discharge state, when the external stimulation current intensity increases from 0 to 30, the discharge mode of neurons experiences period-2 bursting state and chaotic discharge, then evolves to period-2 bursting state and finally to spiking state.
In the first layer of the network, the influence of the random boundary value of the matrix network on the spatiotemporal pattern at different times can be observed due to the influence of the external current and the random boundary, as shown in Figure 5. It can be seen that the first layer generates spiral wave induced by random values of boundary under appropriate coupling intensity and external force, and the second layer is in the different states. When the time is relatively small, the advance of traveling wave induced by random boundary can be observed, and some spiral seeds and broken spiral waves gradually appear. As time increases, some single armed and double armed spiral seeds can also be observed.
Coupling channels between the two layers are set in multiple areas and spiral wave of first layer affect second layer via the coupling channels. When only one connection channel is set in the middle of the first layer network, the spatiotemporal pattern of the second layer network at different time units are shown in Figure 6. An obvious characteristic of Figure 6 is that all the images at different times in the figure are similar to target waves, except for breakage at 500 time units. Target-like waves tend to move from the boundary to the center of the matrix network, thus different spatiotemporal patterns can be seen at different times. Since the change of the spatiotemporal pattern is not of analytical significance when changing the observation time, we chose to fix the time at 5000 time units and study the impact of changing the neuron coupling strength on the spatiotemporal pattern in the later research.
As shown in Figure 7, the developed pattern of second layer is calculated under different coupling intensity D2 at t = 5000 time units when the two layers connect at one local areas. It can be observed that the random boundary value propagates from the boundary of the matrix network to the middle area and after the collision traveling wave it starts to form a spiral seed and a target-like wave in the case of D2 = 0.9. Continue to increase the coupling strength between neurons in the second layer, it can be found that the target-like wave breaks and gradually disappears, and then a target-like wave with high potential is observed in the central area of the network. By further increasing the coupling strength between neurons, it can be observed that the spatiotemporal patterns of the neural network experience the process of convergence → diffusion → re-convergence → re-diffusion. Moreover, when the coupling strength D2 = 3.0, only the neurons in the central area of the neural network are in the discharge state, and the neurons in other areas are inhibited, which may be caused by too large coupling strength.
According to Eq (4), the synchronization factor is calculated when the coupling strength between neurons in the second layer neural network is changed, and the results are shown in Figure 8. Once you see the image of synchronization factor, you can easily analyze that it is an inverted bell-like shaped curve. In other words, there is an optimal D2 value, which can minimize the synchronization factor R. Further analysis shows that when the value of D2 is between 0.4 and 1.3, the synchronization factors are close to the minimum, which means that the synchronization of the neural network is low, but an unexpected phenomenon is that spiral seeds or target-like waves appear in the network. In essence, it can be understood that the synchronization factor is small in this case, and the spatiotemporal patterns of the neural system are orderly arranged, so we can see the above situation. When the coupling strength is small or large, the neural system tends to be stay in homogeneous state, so the synchronization factor is also large. These results are consistent with the previous statement.
In the following research, the influence of the inter-layer connection strength on the spatiotemporal patterns in the second layer neural network is studied. As shown in Figure 9, the developed pattern of second layer is calculated under different inter-layer connection strength k at t = 5000 time units when the two layers connect at one local areas. It can be observed that when the inter-layer connection strength is small (for k = 0.1 and k = 0.1), all neurons are basically in a resting state, and only a few neurons oscillate below the threshold. With the increase of inter-layer coupling strength, several relatively complete target-like wave patterns can be observed, and then these ring patterns begin to break, showing a random change rule. This means that the random boundary plays a greater role than the inter-layer coupling strength in the process of generating spiral wave patterns.
The change trend of synchronization factor when changing the inter-layer coupling strength is shown in Figure 10. It is easy to find that the curve is a monotone decreasing function of the inter-layer coupling strength, which indicates that the synchronization of neurons in the neural network will rapidly decrease to a value close to 0 with the increase of the inter-layer coupling strength, and then basically keep fluctuating in a small range. By comparing Figures 9 and 10, it can be found that the spatiotemporal pattern shows that there is neither spiral wave nor spiral seed when the synchronization factor is relatively large. With the further increase of the inter-layer coupling strength, the synchronization factor decreases rapidly, and the target-like wave appears in the spatiotemporal pattern. When the inter-layer coupling strength is further increased and the synchronization factor is reduced to close to 0, the occurred target-like waves begin to break. As mentioned earlier, this is also related to the randomness of random boundary values.
In order to understand the impact of the number of connecting channels between layers on the spatiotemporal pattern of the bi-layer neural network, the spatiotemporal pattern in the case of two area connections are studied as well, as shown in Figure 11. From Figure 11(a)–(c), it can be clearly seen that the specific positions of the two connecting regions are obviously different from their surrounding regions, which indicates that the spiral waves generated by the first layer of network have a great impact on the spatiotemporal pattern of the second layer. When the connection strength D2 between neurons in the second layer increases to 0.4, obvious spiral seeds can be observed, including single arm and double arm. If we continue to increase the connection strength between neurons, it is found that the spatiotemporal pattern starts to change into the form of target-like waves, but the appearance of complete target-like waves can hardly be observed. Another obvious feature is that the area where the two connecting channels are located has a great impact on the neural network of the second layer, and no matter what the coupling strength is, this impact basically exists.
When the two layers are connected at two local areas, the synchronization factor changing with the coupling intensity of the second layer is shown in Figure 12, and the curve in the form of anti-resonance is observed again in this image. It can be found that when the variation range of D2 is from 0.5 to 1.2, the value of synchronization factor R is the relative minimum, however the synchronization factor is relatively large when the coupling strength is small. This anti-resonance curve shows that there is a certain range of coupling strength that can inhibit the synchronization of neurons in the neural network.
Considering that the inter-layer connection strength also affects the spatiotemporal dynamics of the second layer neural network, the developed pattern of second layer under different inter-layer connection strength is shown in Figure 13. At the same time, if you compare Figures 13 and 14, one can find the same rule as described above, that is, the synchronization factor is relatively large when the inter-layer coupling strength is small, and no spiral seeds or spiral waves are displayed in the network. If the inter-layer coupling strength k is selected as 0.2, the synchronization factor shows a maximum value that is not very obvious, showing the form of "resonance". If the inter-layer connection strength is continued to increase, it will again get the result that the synchronization factor decreases rapidly and remains near zero. However, no matter whether the value of synchronization factor is 0 or not, broken spiral seeds and incomplete target like waves can be observed in the spatiotemporal model. Among these broken spiral seeds, the traces of the connection channels in the two areas can also be clearly seen, indicating that the location of the connection channels really plays a crucial role in the spatiotemporal pattern in the second layer network.
In the following research, the bisection and trisection points of the spatial location in the second layer matrix networks are found out and they are determines as three connection areas to connect the two layers respectively. The developed pattern of second layer is calculated under different coupling intensity D2 at t = 5000 time units when the two layers connect at three local areas, which is plotted in Figure 15. It can be observed that the three connection regions have a great impact on the spatiotemporal pattern of the second layer neural network. The diversity of the spatiotemporal pattern is derived from the changes in the signals of the three connection channels. When the signals come from the three channels collide and interact, rich and varied spatiotemporal pattern can be generated in the entire network. With the generation of the collision and the further propagation of the signal, the broken spiral wave and double armed spiral seed appear. If the coupling strength between neurons is further increased, the appearance of target-like waves can also be observed, and a large number of neurons may appear in some certain regular area and discharge intensively at the same time.
Synchronization factor varies with coupling intensity of second layer D2 is plotted in Figure 16 when the two layers connect at three local areas. As in the previous R-D2 curves, the relationship between synchronization factor and coupling strength between neurons in the second layer presents the form of "anti-resonance", indicating that there is a certain range of coupling strength which can make the synchronization of neural network the lowest. When the connection channel becomes four regions, the same conclusion can be obtained, as shown in Figure 20. In general, the synchronization factor changes with D2 in the form of anti-resonance when the coupling strength between neurons in the second layer is changed. This conclusion is valid regardless of how many connection channels exist between the two layers of networks.
If the inter-layer connection strength is changed, whether the connection channels between the two layers of neural networks are selected as 3 or 4, the overall change trend of the spatiotemporal pattern is basically consistent, as shown in Figures 17 and 21. No matter how many connection areas are between the two layers, one can clearly observe the effect of the connection channel on the second layer network in the spatiotemporal pattern. In the above two cases, the synchronization factor changing with inter-layer connection strength are plotted in Figures 18 and 22. Although these two figures are basically consistent with the previous R-k curve on the whole, they show resonance peaks similar to resonance when the inter-layer coupling strength is small. Under other inter-layer coupling strengths, there is still a random transition from the broken spiral seed to the target-like wave.
Figure 19 shows the spatiotemporal pattern when there are four connection regions between the two layers of networks. It can be seen that when the coupling strength of neurons in the second layer is 0.1, the spatiotemporal pattern in the network presents a basically symmetrical pattern. Then, this symmetry is destroyed with the increase of coupling strength between neurons. The reason may be that the increase of coupling strength also leads to the increase of randomness, which leads to the loss of symmetry. In the subsequent images, spiral seeds and broken target-like waves can still be observed. The analysis of these images tells us that when the coupling strength between neurons in the second layer is changed, the spatiotemporal patterns in the second layer of neural network may appear in a variety of forms, and the corresponding synchronization factor may appear in the form of anti- resonance with the change of coupling strength. On the other hand, a curve that is a monotone decreasing function as a whole can be observed when the coupling strength between layers is changed. Figures 20–22 have been discussed previously and will not be repeated here.
Based on four kinds of firing patterns of Izhikevich neuronal model, a bi-layer neural network with multi-channel connection is constructed using numerical simulation. Each layer of network is composed of 200 × 200 Izhikevich neurons to model a two-dimensional matrix network, and random functions are applied at the boundary of first layer of network. The spiral wave in the first layer is transmitted to the second layer through multiple channels, and the spatiotemporal mode and network synchronization in the second layer are studied. By means of simulation, four different firing modes of neurons are discussed respectively, and the bifurcation of the membrane potential are studied. Then the formation and breaking mechanism of spiral wave in the two-layer network are studied as the coupling strength between neurons in the second layer and the inter-layer connection strength increase. With the development of the research, the synchronization properties of neural networks are explored by changing the coupling strength between neurons in the second layer as well as the inter-layer connection strength.
Obtained results indicate that only when the firing mode of the Izhikevich neuron constituting the matrix neural network is in RS state, the emergence and disappearance of spiral waves can be observed in the network, while the formation of spiral wave seeds cannot be observed in the network if it is composed of other firing modes such as FS, IB and CH. Further research shows that the variation of synchronization factor with coupling strength between adjacent neurons in the second layer shows an inverse bell-like curve in the form of "inverse resonance", but the variation of synchronization factor with inter-layer connection strength is a curve which is approximately monotonically decreasing. No matter how many connection areas are between the two layers of neural networks, spiral waves can be observed in the second layer of matrix neural networks, even the appearance of target-like waves can be observed under some certain conditions. Furthermore, we find that the more important phenomenon is that lower synchronicity is helpful to develop spatiotemporal patterns.
The study of the coupling strength and the influence of random boundaries in the neural networks constructed by Izhikevich neurons on the spatiotemporal behavior can be instructive to further explore the signal propagation in the cerebral cortical neural networks.
This work is supported by Science and Technology Project of Jiangxi Provincial Department of Education under Grants No. GJJ203111, Science and Technology Project of Yuzhang Normal University under Grants No. YZYB-21-17, and Talent Introduction Project (No. NGRCZX-22-07), Special Fund for Doctor of Science and Technology Program of Nanchang Institute of Science and Technology under Grants No. NGKJ-21-03.
The authors declare there is no conflict of interest.
Guowei Wang: Conceptualization, Methodology, Software, Visualization, Investigation, Writing- Reviewing and Editing; Yan Fu: Data curation, Writing-Original draft preparation, Supervision, Software, Validation.
The authors confirm that the data supporting the findings of this study are available within the article.
[1] | Bogoeva-Gaceva G, Dimeski D, Srebrenkoska V (2013) Biocomposites based on poly(lactic acid) and kenaf fibers: Effect of micro-fibrillated cellulose. Maced J Chem Chem En 32: 331–335. |
[2] | Dundar T, Ayrilmis N, Büyüksari U (2010) Utilization of waste pine cone in manufacture of wood/plastic composite. Second International Conference on Sustainable Construction Materials and Technologies, Ancona, Italy. |
[3] | Kabir MM, Wang H, Aravinthan T, et al. (2011) Effects of natural fibre surface on composite properties: A review. Proceedings of the 1st international postgraduate conference on engineering, designing and developing the built environment for sustainable wellbeing (eddBE2011), Queensland University of Technology, 94–99. |
[4] |
Mittal G, Dhand V, Rhee KY, et al. (2015) A review on carbon nanotubes and graphene as fillers in reinforced polymer nanocomposites. J Ind Eng Chem 21: 11–25. doi: 10.1016/j.jiec.2014.03.022
![]() |
[5] |
Reddy MM, Vivekanandhan S, Misra M, et al. (2013) Biobased plastics and bionanocomposites: Current status and future opportunities. Prog Polym Sci 38: 1653–1689. doi: 10.1016/j.progpolymsci.2013.05.006
![]() |
[6] |
Rahman A, Ali I, Al Zahrani SM, et al. (2011) A review of the applications of nanocarbon polymer composites. Nano 6: 185–203. doi: 10.1142/S179329201100255X
![]() |
[7] |
Ogunsona EO, Misra M, Mohanty AK (2017) Impact of interfacial adhesion on the microstructure and property variations of biocarbons reinforced nylon 6 biocomposites. Compos Part A-Appl S 98: 32–44. doi: 10.1016/j.compositesa.2017.03.011
![]() |
[8] |
Bledzki AK, Gassan J (1999) Composites reinforced with cellulose based fibres. Prog Polym Sci 24: 221–274. doi: 10.1016/S0079-6700(98)00018-5
![]() |
[9] | Cyras VP, Iannace S, Kenny JM, et al. (2001) Relationship between processing and properties of biodegradable composites based on PCL/starch matrix and sisal fibers. Polym Composite 22: 104–110. |
[10] |
Chiellini E, Cinelli P, Chiellini F, et al. (2004) Environmentally degradable bio-based polymeric blends and composites. Macromol Biosci 4: 218–231. doi: 10.1002/mabi.200300126
![]() |
[11] | Lee SG, Choi SS, Park WH, et al. (2003) Characterization of surface modified flax fibers and their biocomposites with PHB. Macromolecular symposia, Weinheim: WILEY-VCH Verlag, 197: 89–100. |
[12] |
Barari B, Omrani E, Moghadam AD, et al. (2016) Mechanical, physical and tribological characterization of nano-cellulose fibers reinforced bio-epoxy composites: an attempt to fabricate and scale the 'Green' composite. Carbohyd Polym 147: 282–293. doi: 10.1016/j.carbpol.2016.03.097
![]() |
[13] | Omrani E, Menezes PL, Rohatgi PK (2016) State of the art on tribological behavior of polymer matrix composites reinforced with natural fibers in the green materials world. JESTECH 19: 717–736. |
[14] |
Varma IK, Krishnan SRA, Krishnamoorthy S (1989) Composites of glass/modified jute fabric and unsaturated polyester resin. Composites 20: 383–388. doi: 10.1016/0010-4361(89)90664-2
![]() |
[15] |
Valadez-Gonzalez A, Cervantes-Uc JM, Olayo R, et al. (1999) Effect of fiber surface treatment on the fiber–matrix bond strength of natural fiber reinforced composites. Compos Part B-Eng 30: 309–320. doi: 10.1016/S1359-8368(98)00054-7
![]() |
[16] | Mohanty AK, Misra MA, Hinrichsen G (2000) Biofibres, biodegradable polymers and biocomposites: An overview. Macromol Mater Eng 276: 1–24. |
[17] | Wallenberger FT, Weston N (2003) Natural fibers, plastics and composites, Springer Science & Business Media. |
[18] | Cicala G, Cristaldi G, Recca G, et al. (2010) Composites based on natural fibre fabrics. In: Woven fabric engineering, InTech. |
[19] |
Abdelmouleh M, Boufi S, Belgacem MN, et al. (2007) Short natural-fibre reinforced polyethylene and natural rubber composites: effect of silane coupling agents and fibres loading. Compos Sci Technol 67: 1627–1639. doi: 10.1016/j.compscitech.2006.07.003
![]() |
[20] | Joseph PV (2001) Studies on short sisal fibre reinforced isotactic polypropylene composites. |
[21] | Taj S, Munawar MA, Khan S (2007) Natural fiber-reinforced polymer composites. Proc Pakistan Acad Sci 44: 129–144. |
[22] | Sutton A, Black D, Walker P (2011) Natural Fiber Insulation: An Introduction to Low Impact Building Material, IHS BRE Press. |
[23] |
Dhakal HN, Zhang ZY, Richardson MOW (2007) Effect of water absorption on the mechanical properties of hemp fibre reinforced unsaturated polyester composites. Compos Sci Technol 67: 1674–1683. doi: 10.1016/j.compscitech.2006.06.019
![]() |
[24] |
Razak NIA, Ibrahim NA, Zainuddin N, et al. (2014) The influence of chemical surface modification of kenaf fiber using hydrogen peroxide on the mechanical properties of biodegradable kenaf fiber/poly(lactic acid) composites. Molecules 19: 2957–2968. doi: 10.3390/molecules19032957
![]() |
[25] |
Misra S, Misra M, Tripathy SS, et al. (2002) The influence of chemical surface modification on the performance of sisal-polyester biocomposites. Polym Composite 23: 164–170. doi: 10.1002/pc.10422
![]() |
[26] |
Krishnaiah P, Ratnam CT, Manickam S (2017) Enhancements in crystallinity, thermal stability, tensile modulus and strength of sisal fibres and their PP composites induced by the synergistic effects of alkali and high intensity ultrasound (HIU) treatments. Ultrason Sonochem 34: 729–742. doi: 10.1016/j.ultsonch.2016.07.008
![]() |
[27] |
Asim M, Jawaid M, Abdan K, et al. (2016) Effect of alkali and silane treatments on mechanical and fibre-matrix bond strength of kenaf and pineapple leaf fibres. J Bionic Eng 13: 426–435. doi: 10.1016/S1672-6529(16)60315-3
![]() |
[28] | Mathew L, Joseph KU, Rani J (2004) Isora fibres and their composites with natural rubber. Prog Rubber Plast Re 20: 337. |
[29] |
Mohanty AK, Misra M, Drzal LT (2001) Surface modifications of natural fibers and performance of the resulting biocomposites: an overview. Compos Interface 8: 313–343. doi: 10.1163/156855401753255422
![]() |
[30] | Fengel D, Wegener G (1983) Wood: chemistry, ultrastructure, reactions, Walter de Gruyter. |
[31] |
Behera AK, Avancha S, Basak RK, et al. (2012) Fabrication and characterizations of biodegradable jute reinforced soy based green composites. Carbohyd Polym 88: 329–335. doi: 10.1016/j.carbpol.2011.12.023
![]() |
[32] |
Aziz SH, Ansell MP (2004) The effect of alkalization and fibre alignment on the mechanical and thermal properties of kenaf and hemp bast fibre composites: Part 1-polyester resin matrix. Compos Sci Technol 64: 1219–1230. doi: 10.1016/j.compscitech.2003.10.001
![]() |
[33] |
Xie Y, Hill CA, Xiao Z, et al. (2010) Silane coupling agents used for natural fiber/polymer composites: A review. Compos Part A-Appl S 41: 806–819. doi: 10.1016/j.compositesa.2010.03.005
![]() |
[34] |
Devi LU, Bhagawan SS, Thomas S (1997) Mechanical properties of pineapple leaf fiber-reinforced polyester composites. J Appl Polym Sci 64: 1739–1748. doi: 10.1002/(SICI)1097-4628(19970531)64:9<1739::AID-APP10>3.0.CO;2-T
![]() |
[35] |
George J, Sreekala MS, Thomas S (2001) A review on interface modification and characterization of natural fiber reinforced plastic composites. Polym Eng Sci 41: 1471–1485. doi: 10.1002/pen.10846
![]() |
[36] |
Wang B, Panigrahi S, Tabil L, et al. (2007) Pre-treatment of flax fibers for use in rotationally molded biocomposites. J Reinf Plast Comp 26: 447–463. doi: 10.1177/0731684406072526
![]() |
[37] |
Valadez-Gonzalez A, Cervantes-Uc JM, Olayo R, et al. (1999) Chemical modification of henequen fibers with an organosilane coupling agent. Compos Part B-Eng 30: 321–331. doi: 10.1016/S1359-8368(98)00055-9
![]() |
[38] |
Rong MZ, Zhang MQ, Liu Y, et al. (2001) The effect of fiber treatment on the mechanical properties of unidirectional sisal-reinforced epoxy composites. Compos Sci Technol 61: 1437–1447. doi: 10.1016/S0266-3538(01)00046-X
![]() |
[39] | Tserki V, Zafeiropoulos NE, Simon F, et al. (2005) A study of the effect of acetylation and propionylation surface treatments on natural fibres. Compos Part A-Appl S 36: 1110–1118. |
[40] |
Hill CA, Khalil HA, Hale MD (1998) A study of the potential of acetylation to improve the properties of plant fibres. Ind Crop Prod 8: 53–63. doi: 10.1016/S0926-6690(97)10012-7
![]() |
[41] |
Bledzki AK, Mamun AA, Lucka-Gabor M, et al. (2008) The effects of acetylation on properties of flax fibre and its polypropylene composites. Express Polym Lett 2: 413–422. doi: 10.3144/expresspolymlett.2008.50
![]() |
[42] |
Paul S, Nanda P, Gupta R (2003) PhCOCl-Py/basic alumina as a versatile reagent for benzoylation in solvent-free conditions. Molecules 8: 374–380. doi: 10.3390/80400374
![]() |
[43] |
Nair KM, Thomas S, Groeninckx G (2001) Thermal and dynamic mechanical analysis of polystyrene composites reinforced with short sisal fibres. Compos Sci Technol 61: 2519–2529. doi: 10.1016/S0266-3538(01)00170-1
![]() |
[44] | Li X, Tabil LG, Panigrahi S (2007) Chemical treatments of natural fiber for use in natural fiber-reinforced composites: a review. J Polym Environ 15: 25–33. |
[45] |
Joseph K, Thomas S, Pavithran C (1996) Effect of chemical treatment on the tensile properties of short sisal fiber reinforced polyethylene composites. Polymer 37: 5139–5149. doi: 10.1016/0032-3861(96)00144-9
![]() |
[46] | Rao CHC, Madhusudan S, Raghavendra G, et al. (2012) Investigation in to wear behavior of coir fiber reinforced epoxy composites with the Taguchi method. Int J Eng Res Appl 2: 2248–9622. |
[47] |
Paul A, Joseph K, Thomas S (1997) Effect of surface treatments on the electrical properties of low-density polyethylene composites reinforced with short sisal fibers. Compos Sci Technol 57: 67–79. doi: 10.1016/S0266-3538(96)00109-1
![]() |
[48] |
Khan MA, Hassan MM, Drzal LT (2005) Effect of 2-hydroxyethyl methacrylate (HEMA) on the mechanical and thermal properties of jute-polycarbonate composite. Compos Part A-Appl S 36: 71–81. doi: 10.1016/S1359-835X(04)00178-2
![]() |
[49] | Aji IS, Sapuan SM, Zainudin ES, et al. (2009) Kenaf fibres as reinforcement for polymeric composites: a review. Int J Mech Mater Eng 4: 239–248. |
[50] |
Wang B, Panigrahi S, Tabil L, et al. (2007) Pre-treatment of flax fibers for use in rotationally molded biocomposites. J Reinf Plast Comp 26: 447–463. doi: 10.1177/0731684406072526
![]() |
[51] | Rivera-Armenta JL, Flores-Hernández CG, Del Angel-Aldana RZ, et al. (2012) Evaluation of graft copolymerization of acrylic monomers onto natural polymers by means infrared spectroscopy, In: Infrared Spectroscopy-Materials Science, Engineering and Technology, InTech. |
[52] | Singha AS, Rana AK (2012) A comparative study on functionalization of cellulosic biofiber by graft copolymerization of acrylic acid in air and under microwave radiation. BioResources 7: 2019–2037. |
[53] | Wang B, Panigrahi S, Crerar W, et al. (2003) Application of pre-treated flax fibers in composites. CSAE/SCGR Paper No. 03-367. |
[54] |
George J, Sreekala MS, Thomas S (2001) A review on interface modification and characterization of natural fiber reinforced plastic composites. Polym Eng Sci 41: 1471–1485. doi: 10.1002/pen.10846
![]() |
[55] |
Kalia S, Kaith BS, Kaur I (2009) Pretreatments of natural fibers and their application as reinforcing material in polymer composites-a review. Polym Eng Sci 49: 1253–1272. doi: 10.1002/pen.21328
![]() |
[56] | Wallenberger FT, Weston N (2004) Natural fibers plastics and composites, Springer Science & Business Media. |
[57] |
Rahman MM, Mallik AK, Khan MA (2007) Influences of various surface pre-treatments on the mechanical and degradable properties of photo grafted oil palm fibres. J Appl Polym Sci 105: 3077–3086. doi: 10.1002/app.26481
![]() |
[58] |
Paul SA, Joseph K, Mathew GG, et al. (2010) Influence of polarity parameters on the mechanical properties of composites from polypropylene fiber and short banana fiber. Compos Part A-Appl S 41: 1380–1387. doi: 10.1016/j.compositesa.2010.04.015
![]() |
[59] | Kiattipanich N, Kreua-Ongarjnukool N, Pongpayoon T, et al. (2007) Properties of polypropylene composites reinforced with stearic acid treated sugarcane fiber. J Polym Eng 27: 411–428. |
[60] |
Kalaprasad G, Francis B, Thomas S, et al. (2004) Effect of fibre length and chemical modifications on the tensile properties of intimately mixed short sisal/glass hybrid fibre reinforced low density polyethylene composites. Polym Int 53: 1624–1638. doi: 10.1002/pi.1453
![]() |
[61] |
Torres FG, Cubillas ML (2005) Study of the interfacial properties of natural fibre reinforced polyethylene. Polym Test 24: 694–698. doi: 10.1016/j.polymertesting.2005.05.004
![]() |
[62] | Pickering KL, Li Y, Farrell RL, et al. (2007) Interfacial modification of hemp fiber reinforced composites using fungal and alkali treatment. J Biobased Mater Bio 1: 109–117. |
[63] | Jafari MA, Nikkhah A, Sadeghi AA, et al. (2007) The effect of Pleurotus spp. fungi on chemical composition and in vitro digestibility of rice straw. Pak J Biol Sci 10: 2460–2464. |
[64] | Kolybaba M, Tabil LG, Panigrahi S, et al. (2006) Biodegradable polymers: past, present, and future. ASABE/CSBE North Central Intersectional Meeting, American Society of Agricultural and Biological Engineers. |
[65] | Mohanty AK, Misra MA, Hinrichsen G (2000) Biofibres, biodegradable polymers and biocomposites: An overview. Macromol Mater Eng 276: 1–24. |
[66] |
Riedel U, Nickel J (1999) Natural fibre-reinforced biopolymers as construction materials-new discoveries. Die Angewandte Makromolekulare Chemie 272: 34–40. doi: 10.1002/(SICI)1522-9505(19991201)272:1<34::AID-APMC34>3.0.CO;2-H
![]() |
[67] | Mohanty AK, Wibowo A, Misra M, et al. (2004) Effect of process engineering on the performance of natural fiber reinforced cellulose acetate biocomposites. Compos Part A-Appl S 35: 363–370. |
[68] | Ghanbarzadeh B, Almasi H (2013) Biodegradable Polymer, In: Chamy R, Rosenkranz F, Biodegradation-Life of Science, InTech. |
[69] |
Albertsson AC, Karlsson S (1995) Degradable polymers for the future. Acta Polym 46: 114-123. doi: 10.1002/actp.1995.010460203
![]() |
[70] |
Chandra R, Rustgi R (1998) Biodegradable polymers. Prog Polym Sci 23: 1273–1335. doi: 10.1016/S0079-6700(97)00039-7
![]() |
[71] |
Liu D, Tian H, Zhang L, et al. (2008) Structure and properties of blend films prepared from castor oil-based polyurethane/soy protein derivative. Ind Eng Chem Res 47: 9330–9336. doi: 10.1021/ie8009632
![]() |
[72] | Kumar R, Liu D, Zhang L (2008) Advances in proteinous biomaterials. J Biobased Mater Bio 2: 1–24. |
[73] |
Cao X, Chen Y, Chang PR, et al. (2008) Green composites reinforced with hemp nanocrystals in plasticized starch. J Appl Polym Sci 109: 3804–3810. doi: 10.1002/app.28418
![]() |
[74] |
Vazquez A, Dominguez VA, Kenny JM (1999) Bagasse fiber-polypropylene based composites. J Thermoplast Compos 12: 477–497. doi: 10.1177/089270579901200604
![]() |
[75] |
Chen X, Guo Q, Mi Y (1998) Bamboo fiber-reinforced polypropylene composites: A study of the mechanical properties. J Appl Polym Sci 69: 1891–1899. doi: 10.1002/(SICI)1097-4628(19980906)69:10<1891::AID-APP1>3.0.CO;2-9
![]() |
[76] |
Thakore IM, Desai S, Sarawade BD, et al. (2001) Studies on biodegradability, morphology and thermo-mechanical properties of LDPE/modified starch blends. Eur Polym J 37: 151–160. doi: 10.1016/S0014-3057(00)00086-0
![]() |
[77] |
Zhang JF, Sun X (2004) Mechanical and thermal properties of poly(lactic acid)/starch blends with dioctyl maleate. J Appl Polym Sci 94: 1697–1704. doi: 10.1002/app.21078
![]() |
[78] | Wattanakornsiri A, Pachana K, Kaewpirom S, et al. (2011) Green composites of thermoplastic corn starch and recycled paper cellulose fibers. Songklanakarin J Sci Technol 33: 461–467. |
[79] | Avérous L (2008) Polylactic acid: synthesis, properties and applications, In: Belgacem MN, Gandini A, Monomers, polymers and composites from renewable resources, Elsevier Ltd., 433–450. |
[80] | Auras R, Harte B, Selke S (2005) Polylactides. A new era of biodegradable polymers for packaging application. In Ann Tech Conf–ANTEC Conf Proc, 8: 320–324. |
[81] |
Oksman K, Skrifvars M, Selin JF (2003) Natural fibres as reinforcement in polylactic acid (PLA) composites. Compos Sci Technol 63: 1317–1324. doi: 10.1016/S0266-3538(03)00103-9
![]() |
[82] |
Fuqua MA, Huo S, Ulven CA (2012) Natural fiber reinforced composites. Polym Rev 52: 259–320. doi: 10.1080/15583724.2012.705409
![]() |
[83] |
Simoes CL, Viana JC, Cunha AM (2009) Mechanical properties of poly(ɛ-caprolactone) and poly(lactic acid) blends. J Appl Polym Sci 112: 345–352. doi: 10.1002/app.29425
![]() |
[84] |
Jain S, Reddy MM, Mohanty AK, et al. (2010) A new biodegradable flexible composite sheet from poly(lactic acid)/poly(ɛ-caprolactone) blends and Micro-Talc. Macromol Mater Eng 295: 750–762. doi: 10.1002/mame.201000063
![]() |
[85] | Tuil R, Fowler P, Lawther M, et al. (2000) Properties of biobased packaging materials. In: Production of Biobased Packaging Materials for the Food Industry, Center for Skov, Landskab og Planlægning/Københavns Universitet. |
[86] |
Lehermeir HJ, Dorgan JR, Way JD (2001) Gas permeation properties of poly(lactic acid). J Membrane Sci 190: 243–251. doi: 10.1016/S0376-7388(01)00446-X
![]() |
[87] | Rasal RM, Hirt DE (2009) Toughness decrease of PLA-PHBHHx blend films upon surface-confined photopolymerization. J Biomed Mater Res A 88: 1079–1086. |
[88] |
Hiljanen-Vainio M, Varpomaa P, Seppälä J, et al. (1996) Modification of poly(L-lactides) by blending: mechanical and hydrolytic behavior. Macromol Chem Phys 197: 1503–1523. doi: 10.1002/macp.1996.021970427
![]() |
[89] |
Sinclair RG (1996) The case for polylactic acid as a commodity packaging plastic. J Macromol Sci A 33: 585–597. doi: 10.1080/10601329608010880
![]() |
[90] |
Martin O, Averous L (2001) Poly(lactic acid): plasticization and properties of biodegradable multiphase systems. Polymer 42: 6209–6219. doi: 10.1016/S0032-3861(01)00086-6
![]() |
[91] |
Ljungberg N, Wesslen B (2002) The effects of plasticizers on the dynamic mechanical and thermal properties of poly(lactic acid). J Appl Polym Sci 86: 1227–1234. doi: 10.1002/app.11077
![]() |
[92] | Pongtanayut K, Thongpin C, Santawitee O (2013) The effect of rubber on morphology, thermal properties and mechanical properties of PLA/NR and PLA/ENR blends. Energy Procedia 34: 888–897. |
[93] | Senawi R, Alauddin SM, Saleh RM, et al. (2013) Polylactic acid/empty fruit bunch fiber biocomposite: Influence of alkaline and silane treatment on the mechanical properties. Int J Biosci Biochem Bioin 3: 59. |
[94] |
Huda MS, Drzal LT, Misra M, et al. (2006) Wood-fiber-reinforced poly(lactic acid) composites: evaluation of the physico-mechanical and morphological properties. J Appl Polym Sci 102: 4856–4869. doi: 10.1002/app.24829
![]() |
[95] |
Oksman K, Skrifvars M, Selin JF (2003) Natural fibres as reinforcement in polylactic acid (PLA) composites. Compos Sci Technol 63: 1317–1324. doi: 10.1016/S0266-3538(03)00103-9
![]() |
[96] |
Tang G, Zhang R, Wang X, et al. (2013) Enhancement of flame retardant performance of bio-based polylactic acid composites with the incorporation of aluminum hypophosphite and expanded graphite. J Macromol Sci A 50: 255–269. doi: 10.1080/10601325.2013.742835
![]() |
[97] | Huang SJ, Edelman PG (1995) An overview of biodegradable polymers and biodegradation of polymers, In: Scott G, Gilead D, Degradable Polymers, Dordrecht: Springer. |
[98] |
Wu CS (2003) Physical properties and biodegradability of maleated-polycaprolactone/starch composite. Polym Degrad Stabil 80: 127–134. doi: 10.1016/S0141-3910(02)00393-2
![]() |
[99] |
Nair LS, Laurencin CT (2007) Biodegradable polymers as biomaterials. Prog Polym Sci 32: 762–798. doi: 10.1016/j.progpolymsci.2007.05.017
![]() |
[100] |
Honma T, Zhao L, Asakawa N, et al. (2006) Poly(ɛ-caprolactone)/chitin and poly(ɛ-caprolactone)/chitosan blend films with compositional gradients: fabrication and their biodegradability. Macromol Biosci 6: 241–249. doi: 10.1002/mabi.200500216
![]() |
[101] |
Wu CS (2005) A comparison of the structure, thermal properties, and biodegradability of polycaprolactone/chitosan and acrylic acid grafted polycaprolactone/chitosan. Polymer 46: 147–155. doi: 10.1016/j.polymer.2004.11.013
![]() |
[102] |
Wu KJ, Wu CS, Chang JS (2007) Biodegradability and mechanical properties of polycaprolactone composites encapsulating phosphate-solubilizing bacterium Bacillus sp. PG01. Process Biochem 42: 669–675. doi: 10.1016/j.procbio.2006.12.009
![]() |
[103] | Jha K, Chamoli S, Tyagi YK, et al. (2018) Characterization of biodegradable composites and application of preference selection index for deciding optimum phase combination. Mater Today Proc 5: 3353–3360. |
[104] |
Steinbüchel A, Valentin HE (1995) Diversity of bacterial polyhydroxyalkanoic acids. FEMS Microbiol Lett 128: 219–228. doi: 10.1111/j.1574-6968.1995.tb07528.x
![]() |
[105] |
Reis KC, Pereira J, Smith AC, et al. (2008) Characterization of polyhydroxybutyrate-hydroxyvalerate (PHB-HV)/maize starch blend films. J Food Eng 89: 361–369. doi: 10.1016/j.jfoodeng.2008.04.022
![]() |
[106] |
Barham PJ, Organ SJ (1994) Mechanical properties of polyhydroxybutyrate-hydroxybutyrate-hydroxyvalerate copolymer blends. J Mater Sci 29: 1676–1679. doi: 10.1007/BF00368945
![]() |
[107] | Mohammadi M, Ghaffari-Moghaddam M (2014) Recovery and Extraction of Polyhydroxyalkanoates (PHAs), In: Polyhydroxyalkanoate (PHA) based Blends, Composites and Nanocomposites, 30: 47. |
[108] |
Macedo JS, Costa MF, Tavares MIB, et al. (2010) Preparation and characterization of composites based on polyhydroxybutyrate and waste powder from coconut fibers processing. Polym Eng Sci 50: 1466–1475. doi: 10.1002/pen.21669
![]() |
[109] |
Fragassa C, de Camargo FV, Pavlovic A, et al. (2018) Experimental evaluation of static and dynamic properties of low styrene emission vinylester laminates reinforced by natural fibres. Polym Test 69: 437–449. doi: 10.1016/j.polymertesting.2018.05.050
![]() |
[110] |
Fragassa C, Pavlovic A, Živković I (2018) The accelerated aging effect of salt water on lignocellulosic fibre reinforced composites. Tribol Ind 40: 1–9. doi: 10.24874/ti.2018.40.01.01
![]() |
[111] | Fragassa C, Pavlovic A, Santulli C (2018) Mechanical and impact characterisation of flax and basalt fibre vinylester composites and their hybrids. Compos Part B-Eng 137: 247–259. |
[112] | Mundera F (2003) Advanced Technology for Processing of NFP for Industrial Applications. 7th International Conference on Wood Plastic Composites, Madison, WI, May 19–20. |
[113] |
Shakoor A, Muhammad R, Thomas NL, et al. (2013) Mechanical and thermal characterisation of poly(l-lactide) composites reinforced with hemp fibers. J Phys Conf Ser 451: 012010. doi: 10.1088/1742-6596/451/1/012010
![]() |
[114] |
Dey K, Ganguly S, Khan RA, et al. (2013) Surface treatment of areca-nut fiber using silane and gamma irradiation: fabrication of polycaprolactone based composite. J Compos Biodegrad Polym 1: 1–7. doi: 10.12974/2311-8717.2013.01.01.1
![]() |
[115] | Khandanlou R, Ahmad MB, Shameli K, et al. (2014) Mechanical and thermal stability properties of modified rice straw fiber blend with polycaprolactone composite. J Nanomater 2014: 93. |
[116] | Sandeep Laxmeshwar S, Viveka S, Madhu Kumar DJ, et al. (2012) Preparation and properties of composite films from modified cellulose fiber-reinforced with PLA. Pharma Chemica 4: 159–168. |
[117] |
Wu CS (2010) Preparation and characterizations of polycaprolactone/green coconut fiber composites. J Appl Polym Sci 115: 948–956. doi: 10.1002/app.30955
![]() |
[118] | Phua YJ, Chow WS, Mohd Ishak ZA (2013) Mechanical properties and structure development in poly(butylene succinate)/organo-montmorillonite nanocomposites under uniaxial cold rolling. Express Polym Lett 5: 93–103. |
[119] |
Barkoula NM, Garkhail SK, Peijs T (2010) Biodegradable composites based on flax/polyhydroxybutyrate and its copolymer with hydroxyvalerate. Ind Crop Prod 31: 34–42. doi: 10.1016/j.indcrop.2009.08.005
![]() |
[120] | Jha K, Tyagi YK, Yadav AS (2018) Mechanical and thermal behaviour of biodegradable composites based on polycaprolactone with pine cone particle. Sādhanā 43: 135. |
1. | G.-W. Wang, Y. Fu, Modes transition and network synchronization in extended Hindmarsh–Rose model driven by mutation of adaptation current under effects of electric field, 2023, 0973-1458, 10.1007/s12648-023-02613-2 | |
2. | Yan Fu, Guowei Wang, Firing patterns transition and network dynamics of an extended Hindmarsh-Rose neuronal system, 2024, 0973-1458, 10.1007/s12648-024-03228-x | |
3. | Ranjit Kumar Upadhyay, Debasish Pradhan, Sanjeev Kumar Sharma, Arnab Mondal, Emergence of spiral and antispiral patterns and its CGLE analysis in leech-heart interneuron model with electromagnetic induction, 2024, 128, 0307904X, 154, 10.1016/j.apm.2024.01.013 | |
4. | Yan Fu, Tian Lu, Meng Zhou, Dongwei Liu, Qihang Gan, Guowei Wang, Effect of color cross-correlated noise on the growth characteristics of tumor cells under immune surveillance, 2023, 20, 1551-0018, 21626, 10.3934/mbe.2023957 | |
5. | Yan Fu, Guowei Wang, Firing patterns transitions and resonance effects of the extended Hindmarsh-Rose neural model with Gaussian noise and transcranial magneto-acousto-electrical stimulation, 2024, 99, 0031-8949, 115266, 10.1088/1402-4896/ad85a0 |
Types | a | b | c | d |
RS | 0.02 | 0.2 | –65 | 8 |
FS | 0.1 | 0.2 | –65 | 2 |
CH | 0.02 | 0.2 | –50 | 2 |
IB | 0.02 | 0.2 | –55 | 4 |