An advanced approach for the electrical responses of discrete fractional-order biophysical neural network models and their dynamical responses

The multiple activities of neurons frequently generate several spiking-bursting variations observed within the neurological mechanism. We show that a discrete fractional-order activated nerve cell framework incorporating a Caputo-type fractional difference operator can be used to investigate the impacts of complex interactions on the surge-empowering capabilities noticed within our findings. The relevance of this expansion is based on the model’s structure as well as the commensurate and incommensurate fractional-orders, which take kernel and inherited characteristics into account. We begin by providing data regarding the fluctuations in electronic operations using the fractional exponent. We investigate two-dimensional Morris–Lecar neuronal cell frameworks via spiked and saturated attributes, as well as mixed-mode oscillations and mixed-mode bursting oscillations of a decoupled fractional-order neuronal cell. The investigation proceeds by using a three-dimensional slow-fast Morris–Lecar simulation within the fractional context. The proposed method determines a method for describing multiple parallels within fractional and integer-order behaviour. We examine distinctive attribute environments where inactive status develops in detached neural networks using stability and bifurcation assessment. We demonstrate features that are in accordance with the analysis’s findings. The Erdös–Rényi connection of asynchronization transformed neural networks (periodic and actionable) is subsequently assembled and paired via membranes that are under pressure. It is capable of generating multifaceted launching processes in which dormant neural networks begin to come under scrutiny. Additionally, we demonstrated that boosting connections can cause classification synchronization, allowing network devices to activate in conjunction in the future. We construct a reduced-order simulation constructed around clustering synchronisation that may represent the operations that comprise the whole system. Our findings indicate the influence of fractional-order is dependent on connections between neurons and the system’s stored evidence. Moreover, the processes capture the consequences of fractional derivatives on surge regularity modification and enhance delays that happen across numerous time frames in neural processing.

The aforementioned models' evaluation revealed multiple broader characteristics, including bifurcations, constraint cycles, and chaotic behaviour 4,5 .Despite the many research investigations that have been conducted to examine the evolving behaviours of continuous-time structures, discrete-time methods and mechanisms have attracted much fewer resources.Discrete-time frameworks have distinguished evolving characteristics, and they can be employed to symbolize a wide range of beneficial complications within the reality of life 6,7 .Because of such features, research into independent neural network mechanisms is critical and has culminated in major discoveries in physical science, mathematical concepts, and numerous other disciplines 8,9 .
Discrete fractional calculus encounters have sparked the interest of many scientists over the past 20 years, owing to their significance for disciplines as diverse as computer application, photo encoding, and encrypted interaction.A plethora of publications on this popular subject have recently been published, with the authors offering an assortment of discrete-time fractional operators, durability assessments, and numerous mathematical findings [10][11][12] .Alsharidi et al. 13 investigated a short-memory discrete fractional difference equation wind turbine model using a Caputo-type fractional difference operator.Al-Qurashi et al. 14 expounded the complexity analysis and discrete fractional difference implementation of the Hindmarsh-Rose neuron system.Chu et al. 15 derived the new configurations of the fuzzy fractional differential Boussinesq model with application in ocean engineering.Wu and Baleanu 16 established their inaugural research, which concerned the modelling of the fractional chaotic map employing the Caputo difference operator and researching its turbulent features.These have led to the proposal of more commensurate fractional discrete chaotic systems, such as 17 , and more incommensurate fractional discrete chaotic systems, such as 18 , along with a variety of control strategies and synchronization schemes between different fractional chaotic maps, such as 19 .
Numerous researchers have been investigating the oscillatory patterns of fractional continuous-time neural network frameworks due to their prospective functions in various disciplines such as pattern identification, computational optimization, connected retention, and data processing 20,21 .However, when implementing continuoustime structures for computerized processing, exploration, or modelling, it is crucial to isolate continuous-time structures regardless of whether they ought to be implemented for technological device estimation, exploration, or modelling.Discrete-time neural networks have been developed and are used in an extensive spectrum of domains 22,23 .Mohamad and Gopalsamy 19 discussed the exponential stability of continuous and discrete-time cellular neural networks with delays.As a result, an examination of the functioning of discrete-time neural network is required.A while ago, a couple of studies were published that investigated the evolving dynamics of fractional discrete-time modelling in neural network [23][24] .Weinberg 25 addressed how membrane conductive storage affects spikes in nerve cells characterized by the fractional-order Hodgkin-Huxley neural network framework.Lundstrom et al. 26 contemplated the fractional operators by neurons in pyramids of the neocortex.Huang et al. 27 expounded the new fractional discrete-time neural network based on the variable order.Allehiany 28 presented the chaotic regulatory and fractional-order dynamics of neural network under electromagnetic radiation.A significant amount of the previously stated discrete-time neural network studies focus on frameworks with corresponding or incommensurate fractional-orders; however, currently, as far as we comprehend, only a handful of improvements have been published that examine the exciting practises of fractional-order discretetime neural network frameworks [29][30][31] .As a consequence, studying the unpredictable behaviour of neural network frameworks driven by fractional differences with commensurate and incommensurate values, in addition to their synchronization and oversight, is an appealing topic.
Influenced by the previous debate, we aim to investigate and analyze the fluctuating patterns of the intriguing two-and three-dimensional discrete fractional-order Morris-Lecar neural network paradigm using electrically powered operations of independent Morris-Lecar firing neural network regarding Groups i and ii activities and slow-fast Morris-Lecar neural network 32,33 concerning its presence design employing commensurate and incommensurate fractional-orders in the present article.The key elements of this discrete fractional-order system will be investigated by implementing hypothetical and computational assessments.Take into account a conductancedependent simulation that investigates the inherent behaviour underpinning fractional-order differences when investigating these terminating features and transformations.Prior studies examined different periods of spiked remarks based on varying setting regimes; despite this, discrete fractional-order can also evaluate broad-sack reactions 34 .The Morris-Lecar frameworks are used because of their respective positions in various reactions, which range from spikes to ruptures.In the parametric forum of the Morris-Lecar approach, it includes two regimes: dominant spiked and swift spiked behaviours.Furthermore, the representation has been transformed into a three-dimensional equivalent in which the generated electricity, ℑ , remains steady but fluctuates with energy.We take the discrete fractional-order technique to be the most significant method in the framework, considering that it modifies gradually, spike adjustments happen, and mixed mode oscillations and mutated monarch butterfly optimizations are observed.It constitutes one of the particularly fascinating neurological variations that arise from electrical activity 35 .Mixed mode oscillations indicate the reversed pathways of tiny and massive intensity resonances 36 .This makes the structure captivating, and its outcome supplies exciting and potentially advantageous uses in an adaptive structure.The resulting appearance of mutated monarch butterfly optimizations results in the creation of an elevated incorporation process.Previously, it was discovered that mixed mode oscillations examined the evolving and neurological behaviour of manoeuvring or breathing 37 .Calcium signalling and electrocardiac applications were both affected 38 .Krupa et al. 39 investigated the process of mixed mode oscillation fluctuations in the mammal cerebral cortex using a two-dimensional representation of dopamine signalling.We additionally glance at how electronic connections affect an assortment of cells that are, in two ways, inactive or periodic.In this case, the synapses of neural networks are considered to be associated via the pathways of the Erdös-Rényi structure.To identify all components in the structure, the interaction causes intricate terminating operations that involve repeated overflowing or surge regularity adaptation mechanisms.Depending on the discovered asymptotic occurrences, a diminished-order emulate is established that can generate the whole system's actions.
Employing the fractional exponent, we discovered commensurate and incommensurate differences in the features of neurologically effective behaviour.Fractional-order power interactions may drastically impact the spikes in the characteristics of various unattached neuron models 40,41 .Practical characteristics are capable of rendering the mathematical framework more responsive to neuronal activity interactions, especially in the network's prospective collaborative behaviour, in which previous dynamic behaviour could impact current conditions.
The remaining content of this article is structured as follows: The qualitative description of the discrete fractional-order neural network system in terms of the electrical nature of independent Morris-Lecar spike neurons via category i and ii activities is presented in a detailed manner.Furthermore, the slow-fast Morris-Lecar cells employing the Caputo difference expression are presented in Section 2. In Section 3, we investigate the framework's fundamentally evolving characteristics using computation and analytical techniques of the twodimensional discrete fractional-order Morris-Lecar neural network model.Section 4 demonstrates the detailed analysis of discrete fractional-order Morris-Lecar neural networks with commensurate and incommensurate fractional-orders.Moreover, the various oscillatory responses are presented with detailed simulation techniques.To stabilize and synchronize the unstable pathways of the suggested fractional discrete-time neural network approach, we put forward responsive processors.Ultimately, Section 6 presents the overall paper's ending.

Mathematical model
In this part, we indicate the fractional-order conductance-based simulation and discuss the presence of multiple neural characteristics 42 .We determine a specific setting structure that enables terminating characteristics via fractional exponent modifications.We investigate the two-and three-dimensional Morris-Lecar frameworks via specific settings and transmit interactions for producing broad surges employing fractional-order behaviour.We chose one of these frameworks to distinguish the influence of fractional formulations on the system's dynamic practises.Morris and Lecar 32 suggested a simple mathematical design that described fluctuations in the barnacle's enormous muscle fibre, which included the cellular electricity formula alongside immediately apparent calcium up-to-date stimulation and a supplementary restoration formula explaining slow potassium energy initialization.In the fractional-order sense, the two-dimensional Morris-Lecar system is explained as follows: A voltage-gated Ca 2+ energy, an induced converter K + energy, and a leaky energy are all part of the biophysi- cally stimulated activated approach.The membrane power interactions are measured by u 1 , and the stimulation attribute of K + charge activity is u 2 .The highest conductivity operates to Ca + , K + and release flows are indicated by the factors W Ca , W K and W L , respectively.The transformations of prospective to various ionised energy performs are F Ca , F K and F L .The protective layer capaciousness is measured by C symbolizes the ambient temperature developing rate for revealing the K + pathway.The settings for F 1 , F 2 , F 3 and F 4 have constant positive values.The letter ℑ denotes the substance that was used stimulus.We want to take consideration of the influence of different inserted energy stimuli on the fractional-order framework based on ϑ ∈ (0, 1].

Diverse dynamical aspects
Initially, suppose the neuron is just starting to fire and is generating spikes as the regulation variable advances cautiously.Slow-fast interactions 42 can be computationally represented as where ẏ(ξ ) = �(y, z) indicates fast spiking and ż(ξ ) = ζ h(y, z) represents the slow modulation.The rapidly changing parameters are represented by y ∈ R r and the slow factors are represented by z ∈ R s , with 0 < ζ ≪ 1 assessing the duration of the split parameter estimation.
The slow-fast three-dimensional Morris-Lecar model is represented by a corresponding structure of ordinary differential equations, where (u 1 , u 2 ) indicate the fast component and u 3 the slow variable.The fractional-order improved three-dimensional Morris-Lecar model 32,42 is shown as follows: The framework has a few additional features.Externally administered energy, which implies index-law interactions in the fractional framework and illustrates the memory impact of the cellular within the voltage 32,42 , is the mechanism factor u 3 .The hyperbolic mappings' settings F 1 , F 2 , F 3 and F 4 are effectively chosen to demonstrate how they are capable of reaching their fixed points instantly.The specified worth value θ is smaller than one, i.e., θ ∈ (0, 1) , and it determines the length of time scale ratio between resonances and modification.To identify this kind of feedback, Lundstrom et al 26 discovered that neuronal pyramids are able to function as fractional differentiation of the challenge intensity envelope.To determine the initial differentiation of a mapping, differentiate twice by implementing the fractional-order derivative of ϑ = 1/2 which leads to the first derivative 43,44 .It adjust- ments the reaction using a decomposing time delay that is proportional to ϑ.
To identify all simulated outcomes, the attribute distinguishes are (for (1), see 33,42 ): 40, ϕ = 0.0670 (for the excited cell simulation of category i ) with various ℑ.Group ii: ℑ = 40.Group ii: ℑ = 45 and for the group ii membrane framework, the specifications are actually identical mentioned previously, aside for Group iii: W Ca = 4.4, F 3 = 2.0, F 4 = 30, ϕ = 0.04 and ℑ = 100.To investigate the functioning of systems, we initially examine neutral states and then bifurcation.afterwards that, we employ the subsequent variable distinguishes for interacting with structures (1) and its improved forms, alongside ℑ(u 3 ) = 0.08 − 0.03u 3 and C = 1 for attribute initiates i, ii and iii, respectively.
Group i: W Ca = 0.9, Group iii: W Ca = 0.9, To investigate the fractional-order Mor- ris-Lecar framework, we employ the well-known description of the fractional difference, namely the Caputo-type fractional difference operator.The commensurate fractional-order framework with the fractional exponent ϑ ∈ (0, 1) is defined as follows: where q 1 , q 2 and q 3 are the fractional-order such that qι ∈ (0, 1).c D q 1 is the Caputo fractional derivative which can be evaluated by the following formula: It is emphasized that when q 1 = q 2 = q 3 , a commensurate fractional-order model is developed; otherwise, an incommensurate fractional-order model is formed.

Preliminaries on discrete fractional calculus
Prior to defining our fractional discrete framework, we will go over certain crucial terms and mathematical principles of discrete fractional calculus.Definition 3.1 ( 10 ) Assume that there is a time scale N a = {a, a + 1, a + 2, . ..}, a ∈ R. The fractional sum of order ϑ for a mapping F can be defined as Definition 3.2 ( 11 ) The ϑ-Caputo fractional difference operator is described as follows where ξ ∈ N a+q−ϑ , ϑ / ∈ N and q = ⌈ϑ⌉ + 1. � q F (ϒ) and (ξ − 1 − ϒ) (q−ϑ−1) represents the m th integer dif- ference operator and falling factorial mapping, respectively, indicated as and Currently, we require the resulting hypothesis 45 to figure out the stability manipulations to stay on the critical highlights of a fractional discrete mechanism with commensurate fractional-order characteristics: Vol.:(0123456789) www.nature.com/scientificreports/Theorem 3.1 ( 44 ) Assume that there is a fractional-order (ξ ) = 1 (ξ ), . . ., q (ξ ) T , ϑ ∈ (0, 1) and W ∈ R q×q , the zero steady states of the commensurate discrete fractional-order framework where ℘ indicates the eigenvalues of the matrix W.
However, the reliability and robustness concept of the dynamic fractional-order incommensurate framework is formulated as follows: If every root of the equation If x 0 = x(0) occurs inside the collection C/K ξ , then the straightforward solution of the structure (9) is locally asymp- totically stable (LAS), where ξ = 1/H, J is the Jacobian matrix of (9) and Remark 3.1 The classical Hartman-Grobman linearization essentially contributes a vital part to the analysis of the locally stable characteristics of a constantly changing system's equilibrium, which contends that the specific behaviour of an unstable system in the neighbourhood of a hyperbolic neutral state is subjectively similar to the behaviour of its linearization within the equilibrium.It is worth noting that a fractional-order equivalent for this linearization theory was discovered in 46 .If X 1 * is a state of equilibrium of (9), i.e. �(X 1 ) , the linearized framework at X 1 * is as follows: where J � (X * 1 ) is the Jacobian matrix of the mapping determined at X * 1 .As a result, the dynamic technique's equilibrium X * 1 is asymptotically stable iff the insignificant algorithm for the dynamical framework ( 9) is asymptotically stable 47 .Also, in view of the well-noted Matignon's principal 47 , the dynamical fractional-order model ( 13) is asymptotically stable iff arg( ) > ϑπ 2 , for some eigenvalue of the Jacobian matrix J � (X * 1 ).

Remark 3.2
According to Morris-Lecar system (3), a steady state X 1 * is known to be a saddle of rank one if an eigenvalue of J � (X 1 * ) is unsteady (that is arg( 1 ) < ϑπ 2 ) and remaining two eigenvalues are stable arg( 2,3 ) > ϑπ 2 .However, the two eigenvalues connected to the equilibrium X 1 * are unstable, when merely one eigenvalue stays constant, the saddle point X 1 * is referred to as the saddle of rank two 48 .

Two-dimensional model
The framework (1) is an example of a two-dimensional fractional-order conductance-dependent activated structure.The discrete version of the fractional-order Morris-Lecar system can be constructed by substituting the Caputo fractional-order formulation c D q ι a 1 in framework (1) with the fractional-order difference scheme defined in (6) c � � ι a 1 , which is demonstrated below: for ξ ∈ N a+1−� .� ι , ι = 1, .., 3 are the fractional-order values such that � ι ∈ (0, 1], ι = 1, 2, 3. Also, u 1 and u 2 are the neuron's cell power and restricting parameter, ℑ is an implemented electricity, Vol:.( 1234567890) www.nature.com/scientificreports/energy, ℓ(u 2 ) is the frequency steady for releasing ionized pathways, and u 2∞ is the amount of accessible channel openings at the state of equilibrium.The numerical procedure for the fractional discrete system can be constructed using the subsequent hypothesis: Theorem 3.3 ( 49 ) Suppose there is a fractional difference equation the unique solution of this initial value problem (15) is presented as where Therefore, the mathematical representation of the discrete fractional-order Morris-Lecar neural networks framework ( 14) is developed using the aforesaid approach: where u 1 (0) and u 2 (0) indicate initial conditions.This is an entirely novel kind of Morris-Lecar system that has "memory impacts." As shown in (18), the configurations of u 1 (n) and u 2 (n) are interdependent concerning every previous factor u 1 (0), u 1 (1), . . ., u 1 (n − 1) and u 2 (0), u 2 (1), . . ., u 2 (n − 1).

Dynamic analysis of Morris-Lecar model
In this section, we will investigate whether the formerly suggested discrete fractional-order Morris-Lecar neural network model ( 19) is steady or in chaos in both instances: commensurate and non-commensurate fractionalorders, respectively.The present research will make use of an assortment of computational techniques, among them bifurcation illustrations and the visual representation of phase profiles in multifaceted estimations.In addition, we use the reduced order evaluation to determine whether or not chaos is present.

Commensurate fractional-order for two-dimensional Morris-Lecar model
In what follows, the robustness of the stable state points within the discrete fractional-order Morris-Lecar neural networks model (18) with commensurate fractional-order is investigated.For the purpose of determining the system's steady state, we consider the model (18): . The mathematical method's solutions are the fixed points of model ( 18): equivalently, The mapping ℑ(u 1 ) has a handful of fundamental characteristics: www.nature.com/scientificreports/ It is critical to emphasize that these features are fulfilled in the specific scenario of the Morris-Lecar neural networks that includes the mapping ℑ(u 1 , u 2 ) supplied by (18).

Stability of steady states
For an undetermined equilibrium state such as u 1r = u 1max or u 1r = u 1min , and considering that is sufficiently small, the variant χ(u 1r ) > 0 is capable of being noticed.As a result, the mapping χ(u 1 ) may have two roots, u 1 ′ ∈ (F K , u 1max ) and u 1 ′′ ∈ (u 1min , 1), respectively.In view of the numerical information, we might additionally suppose that when such roots occur, they constitute something particular within the preceding time frames.
Finally, the stability of equilibria conditions is only affected by fractional-order in the subsequent three scenarios: (c 1 ) the steady state is associated with the initial region and u 1 * ∈ (u 1 ′ , u 1max ); (c 2 ) the steady state is associated with the third region and u 1 * ∈ (u 1min , u 1 ′′ ); (c 3 ) the steady state is associated with the second region is unstable regardless of fractional-order.Now, the threshold ϑ * provided by ( 21) is equivalent to a bifurcation in the present instance (that is., the Jacobian matrix possesses two sets of imaginary type eigenvalues that result in arg( ) = ϑπ 2 ).In a nutshell, contingent on the fractional-order evaluated in framework (18), the states of the Hopf bifurcation are indicated in the (ℑ, u 1 )-plane, located on the primary and/or third branches, accordingly.Evidently, such will have an extremely significant influence on the form of spiked and brimming behaviour in the two-dimensional framework (18) and the three-dimensional slow-fast framework, as will be revealed in the subsection that follows.
As an illustration, take into account ℑ as a bifurcation variable while looking at the bifurcation setting of the conventional two-dimensional Morris-Lecar neural networks framework that includes Hopf-bifurcating highlights (Fig. 1a).Consequently, we examined the influence of fractional-order on the behaviour of the framework (18) and demonstrated the way that it stabilizes the entire structure as the quantity decreases (Fig. 1b).The period depictions of the two-dimensional Morris-Lecar structure using its settings for Group ii are then shown in Fig. 2a-j for various fractional-order outcomes.In the present scenario, the framework has exclusively one unsteady equilibrium point, (u 1 * , u 2 * ) = (6.23101,1/24, 232) (at the point of convergence of the nullclines), which is positioned at the subsequent location.According to the procedure (21), the significant worth of the fractional-order that describes the Hopf-bifurcation at the steady state is ϑ * = 0.89342.A large- amplitude constraint process attractor, similar to spiked behaviour, can be observed in the classical situation, ϑ = 1 .The large-amplitude appealing quasiperiodic restrict process addresses the unsteady state of equilib- rium as the fractional-order reduces, and as it tackles the important threshold for Hopf splitting, an additional quasiperiodic pass over arises, featuring smaller-amplitude fluctuations near the point of stability in addition to large-amplitude rises.When ϑ < ϑ * , the balance of power turns into asynchronously steady.Furthermore, we demonstrate associated data sets to additionally validate the mathematical findings, mentioning that the fractional-order determined thresholds are ϑ * = 0.83241 for group i and ϑ * = 0.96720 for group ii presented in Fig. 3a-d for ℑ = 43, Fig. 3e-h for ℑ = 50 and Fig. 3i-l for ℑ = 90 with the settings of Group iii for various fractional-order outcomes, respectively.

Commensurate fractional-order for three-dimensional Morris-Lecar model
In accordance with the already mentioned features, the three-dimensional slow-fast fractional-order designs (19) are capable of being composed in (20) with the parameter settings Ī(u 1 , u 2 ) as follows: -0.4 -0.www.nature.com/scientificreports/intensity for every rise lowers in tandem bursts (see the bifurcation plots of Fig. 4a,b,d).The distinctive point of reference for the aforementioned parameter setting is a saddle corresponding to measure two.Spike's condition rate modification has been noticed alongside decreasing fractional-order values.Subsequently, following a greater reduction of ϑ = 0.72 , it induces an additional surge in the rate of adjustment (see the bifurcation plots on the left single panel in Fig. 4g).Furthermore, the standard independent neuronal cell simulation exhibits overflow via setting ii.Through reducing ϑ = 0.90 and 0.80, it exhibits different periods of overflowing and spikes in  www.nature.com/scientificreports/behaviour via rise regularity modification and initially surge delays at ϑ = 0.70 (see the bifurcation plots on the left single panel in Fig. 4c,e,h).Ultimately, within group iii, an individual neuronal cell switches from overflowing to swiftly spiked while changing via one to ϑ = 0.95 and 0.85.With ϑ = 0.72 , it enters an uninterrupted state of equilibrium, that is, it settles at the technique's locally asynchronous centre of gravity (see Fig. 4g).

Incommensurate discrete fractional-order two-dimensional Morris-Lecar neural networks model
Here, we examine the evolution of a discrete fractional-order Morris-Lecar neural network simulation using an unpredictable network configuration, whereby every neural network communicates unintentionally and is assigned an interaction likelihood p .To identify computational modelling, we build an Erdös-Rényi network 51 of N = 100 Morris-Lecar oscillating components with an average node-index of �η� ≈ 7 .A comprehensive review of the structure design will be provided in the sections that follow: where W e > 0 denotes the system's electrically powered interaction.The system's connection matrix is expressed by M = (c ι ) N×N .Additionally, we separate a group of size N into two particular groups based on fractional- order expressed as  , where r components possess similar fractional-order, indicating fluctuating behaviour, while the rest of the s points have fractional-order, indicating activated behaviour.Con- sequently, the general population capacity N is given by N = r + s.We begin by investigating the behaviour of independently associated group i activated Morris-Lecar neural networks with two fractional-order factors, ϑ = 1 and = 0.80 .The entire amount of components in the system is N = 100, using r = 65 and s = 35 indi- cating that we are considering a system of neuronal populations with 65% oscillated neural networks and 35% activated nerve cells.In a condition without any association ( W ǫ = 0.001 ), every fluctuating nerve cell in the structure exhibits dominant exploding, while the additional nerve cells stay dormant.In addition to a negligible improvement in electrical power attachment W ǫ = 0.0005 , the fluctuating subgroup continues to have dominant spikes in setting, while an alternate subdivision goes into quiescence.Red and blue lines denote the cumulative information of two independently interconnected components representing two particular groups (see Fig. 5a  and b).The red line connection is picked at random based on the inactive components, and the blue line connection is compared to the increasing nodes in the network.The space-time plot shows that the spikes in components (1-65) are asymptotic (see Fig. 5e,f).When the relationship is increased 100 times ( W ǫ = 0.05 ), the framework's behaviour alters.The previously dormant subgroup is presently exhibiting overflowing behaviour.In particular, the time time-frame between each wave of activity is not constant.An additional subgroup exhibits asymptotic unsteady active spiked (see Fig. 5i,j).It is obvious that as the degree of association in the transformed group increases, the sporadic and additionally dormant structure of the structure disappears, and inconsistent overflowing or jumping comes up.The neural network as a whole exhibits overflowing interactions via a limited number of surges during every explosion and a modest rise of W ǫ = 0.007 (see Fig. 5g,h).Two layers possess distinct intensities; however, identical stages develop in this particular instance.Subsequently, at W ǫ = 1, the interconnected structure displays nearly synchronization practises by switching to active jumping (see Fig. 5k,l).When we improve the momentary stimulation to ℑ = 45, the fluctuating subgroup exhibits swift spikes in while another subsection keeps quiescent.At minimal attachment ( W ǫ = 0.0005 ), neither subgroup's behaviour changes.These include specific groups that begin terminating and exhibit overflowing interactions such that the dependence W ǫ = 0.005 and 0.06 increases.Ultimately, at W ǫ = 1, the sporadic network operation transforms into synchronized regular spikes (Fig. 5c,d).The group ii activated Morris-Lecar neural networks are then studied in an identical structure configuration, featuring a subgroup dormant ( = 0.72 ) and an additional showing mixed mode oscillations ( ϑ = 0.84 ) in their dearth of interaction.Nerve cells are asynchronized when there is an inadequate connection ( W ǫ = 0.0005, 0.05 ).As integrating increases, every link in the structure exhibits surge regularity modification ( W ǫ = 0.05 and W ǫ = 1 ).It is obvious that every one of the two separate arrays of oscillating elements demonstrates nearly perfect synchronization, via analogous steps and intensities (see Fig. 5).

Incommensurate discrete fractional-order three-dimensional Morris-Lecar neural networks model
Here, a modified three-dimensional discrete fractional-order Morris-Lecar neural networks system is presented to demonstrate that the reduction in fractional-order analyses the identical characteristics exhibited by the executed arbitrary network techniques.We discovered an additional accumulation of synchronization status in the framework during the intermediary interaction.Inspired by this information, we are able to compose u 11 = u 12 = u 13 = . . .= u 160 = u 1ϑ and u 161 = u 162 = . . .= u 1100 = u 1� .Given our knowledge of the Erdös-Rényi plot, one can estimate the level of every component or neural network by the standard extent of the system under consideration 52 .As a result, we are able to conclude that the pitch of the  cluster is η  = �η�.The total quantity of spiked oscillating elements in the vicinity for every oscillating device is anticipated to be (1 − ℘ ǫ )η = ℘ 0 η, while the result associated with inactive oscillating elements will be ℘ ǫ η .As a result, we can create a reduced-order simulation involving two oscillating elements, as shown below:  We selected these two structures from two samples to display the time indications.The interval for the assessment of the particular components identified via a blue region is selected from a subclass of components in inactive states (when ).The red indicator is drawn from points that were maintained in jumping assertions in the dearth of interaction.
where ℘ ǫ = s N , ℘ 0 = r n and ℘ 1 = x N are the likelihoods of activated and fluctuating neurons in general.For group i and group ii Morris-Lecar scenarios, we operated computational experiments utilizing lower-order twoconnected frameworks whereby every group is represented through an alike fractional-order exponent demonstrating concurrent behaviour.The computational findings indicate that whenever the two specific populations adhere to collection synchronization, the fluctuations of the reduced-order simulation obey the same structure as the entirety chart (see Fig. 6a-d).For example, overflowing involving two surges might occur for temporary connections in a group i activated framework detached by two layers in the entire system (see Fig. 6e-h).The minimized control system has an analogous bursting structure (see Fig. 6i-l).In the field of neuroscience, Hopf bifurcation occurs when neuron behaviours fluctuate from halting to spiking (the stable, steady options match the state of rest, and the spiking state demonstrates the presence of oscillating solutions).
A further typical form of dynamic activity for a neuron cell happens when, as the controlling value is increased, a saddle point and a limit cycle come together, resulting in a saddle-homoclinic bifurcation.The period of the periodic orbit that seems to be at the point of bifurcation approaches infinity, and as the controlling value is increased more, the periodic orbit disintegrates.It is shown that the formation and demise of saddlehomoclinic bifurcation in model ( 26) depend on implemented energy ℑ ∈ [23, 70] (see Fig. 6), similar to the original model ( 3), but for fractional-order models of orders ϑ = 1.0, 0.95 and 0.90, the neuron requires greater applied current ℑ to bifurcate.When saddle homoclinic bifurcation occurs, neurological research predicts the occurrence or elimination of spiking practices (see Figs. 2 and 4).

Conclusion
We started with the behaviour of just one neuron, which can evolve its reactions to different feedback data analyses contingent on the discrete fractional-order commensurate and incommensurate cases.We addressed the activities of groups i and ii by electrically conducting two-dimensional Morris-Lecar in pragmatic neural networks of commensurate and incommensurate order and reduced orders that are capable of being captured by experimentation.We demonstrated that variations in the discrete fractional-order and system settings influence the stipulates of the suggested approach, and we acquired an assortment of interactions involving steady pathways, regular actions and chaotic behaviours.We investigated the classical-order model's associated bifurcation assessment.The discrete fractional-order may generate dismissing distinctions that are not visible in integerorder fluctuations.With the aid of a reduction in the discrete fractional-order scheme, the potential of the twodimensional Morris-Lecar simulation for varied surge reactions involving mixed mode oscillations and mutated monarch butterfly optimizations via fractional difference is investigated.The findings show that the slow-fast three-dimensional Morris-Lecar in the fractional context produces a variety of packed sequences.Furthermore, it transforms its behaviour from overflowing into a sequence of rises in addition to swift spikes into inconsistent abundance depending on an appropriate collection of factors, which the classical-order simulation is unable to measure for a specific set of factors.The discrete fractional-order commensurate depicts a new illustration of the spiking and exploding reactions with fractional exponent transforms in the unpredictable representations.The method summarized the numerous responses of an individual-activated simulation to an assigned trigger variability using memory-contingent procedures.We examined the significance of electromagnetic connections in an indeterminate system with a subset of nodes in a dormant condition.The whole population would demonstrate rising regularity in adjusting if the fluctuating components were in mutated monarch butterfly optimizations.When the disconnected oscillated components persist in the rapidly active jumping region, all of the population splits into a pair of layers, demonstrating continuous overflowing in temporary association and active shooting in more substantial interaction.We were additionally capable of simplifying the network's complexity into separate, complemented behaviours that effectively represented the changing patterns of the whole system throughout accumulation synchronization, thanks to the clustered synchronization phenomenon.Using the discrete fractional-order incommensurate, we discovered unique impacts on various membrane electrical attributes.Moreover, variations in neurological attributes depend on the memory consequences.The discrete fractional-order difference has potential interactions and can be used to investigate pragmatic occurrences.These findings show that the mathematical framework and networks present a comprehensible strategy for studying neurological behaviour.( 26) https://doi.org/10.1038/s41598-023-45227-8

Figure 1 .
Figure 1.According to group i and ii, bifurcation cases of the two-dimensional Morris-Lecar neural networks.(a) ℑ as a bifurcation factor: ( ℑ = 90, ϑ = 0.97 ) and ( ℑ = 40, ϑ = 0.1 ) indicate the presence of Hopf and saddle component type plots in the framework(18), respectively.The blue and red lines represent the structure's balance and imbalance branches, respectively.(b) reflects the appearance of an unsteady limit process.