This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Emergent noise-aided logic through synchronization

Manaoj Aravind Department of Physics, Indian Institute of Technology Bombay, Powai, Mumbai 400 076, India    Sudeshna Sinha Indian Institute of Science Education and Research Mohali, Knowledge City, SAS Nagar, Sector 81, Manauli, Punjab, PO 140 306, India    P. Parmananda Department of Physics, Indian Institute of Technology Bombay, Powai, Mumbai 400 076, India
Abstract

In this article, we present a dynamical scheme to obtain a reconfigurable noise-aided logic gate, that yields all six fundamental 2-input logic operations, including the XOR operation. The setup consists of two coupled bistable subsystems that are each driven by one subthreshold logic input signal, in the presence of a noise floor. The synchronization state of their outputs, robustly maps to 2-input logic operations of the driving signals, in an optimal window of noise and coupling strengths. Thus the interplay of noise, nonlinearity and coupling, leads to the emergence of logic operations embedded within the collective state of the coupled system. This idea is manifested using both numerical simulations and proof-of-principle circuit experiments. The regions in parameter space that yield reliable logic operations were characterized through a stringent measure of reliability, using both numerical and experimental data.

I Introduction

Exploiting the richness in the behaviour of nonlinear dynamical systems for computational tasks has attracted extensive research interest. Various theoretical schemes have been proposed to realize reconfigurable devices using the rich patterns in chaotic systems and varied experimental implementations of these schemes have been been realized[1, 2, 3].

An interesting recent line of enquiry is the effect of noise in such schemes. This approach is crucial as this leads to the possibility of noise-aided computational devices that can utilize noise to facilitate computing. Logical Stochastic Resonance (LSR)[4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15], is one such scheme in which a bistable system driven by a subthreshold stream of inputs can generate responses consistent with a desired logical operation for an optimal window of noise. This again has been implemented in a wide variety of experimental systems[16, 17, 18, 19] that range from synthetic gene networks[20, 21, 22, 23, 24, 25], optical systems[26, 27] to coulomb coupled quantum dots[28]. More recent efforts have focused on using chaotic attractor hopping [29, 30] and strange non-chaos [31, 32] to yield logic gates.

The effect of coupling on such noisy bistable systems has garnered new interest [33, 34]. Coupling induced LSR [35] demonstrated that logic operations maybe obtained from the state variable of coupled bistable systems in an optimal window of noise.

In this work, we show a novel possibility that the logical output can be embedded in collective state of the coupled system, rather than in the state variable of one of the subsystems. This change leads to a rich construct, which can yield all 2-input logic operations for an optimal window of noise strength and coupling strength.

II Scheme

Refer to caption
Figure 1: A schematic representation of the concept.

Consider coupled bistable systems with two kinds of coupling schemes defined as follows (refer Fig. 1),

x1˙=F(x1)+c(x2x1)+I1(t)+b+Dη1(t)\displaystyle\dot{x_{1}}=F(x_{1})+c(x_{2}-x_{1})+I_{1}(t)+b+D\ \eta_{1}(t)
x2˙=F(x2)+c(x1x2)I2(t)b+Dη2(t)\displaystyle\dot{x_{2}}=F(x_{2})+c(x_{1}-x_{2})-I_{2}(t)-b+D\ \eta_{2}(t) (1)

and

x1˙=F(x1)+c(x2x1)+I1(t)+b+Dη1(t)\displaystyle\dot{x_{1}}=F(x_{1})+c(-x_{2}-x_{1})+I_{1}(t)+b+D\ \eta_{1}(t)
x2˙=F(x2)+c(x1x2)+I2(t)+b+Dη2(t)\displaystyle\dot{x_{2}}=F(x_{2})+c(-x_{1}-x_{2})+I_{2}(t)+b+D\ \eta_{2}(t) (2)

Here, F can be any nonlinear function that yields a bistable potential. The terms η1\eta_{1} and η2\eta_{2} are two uncorrelated, zero mean, univariate Gaussian noises of noise strength DD. The drive signals I1I_{1} and I2I_{2} are two low amplitude (subthreshold) input streams that encode the two binary inputs to the system. A constant asymmetrizing bias bb acts as the tether that morphs the bistable potential, leading to reconfigurability in our scheme. In the above equations, two kinds of coupling terms have been used. In Eq. 1, the coupling between the two subsystems is bidirectional and attractive in nature with coupling strength cc. This typical diffusive kind of coupling interaction [36] tends to synchronize the subsystems i.e. in the context of bistable systems, the two subsystems are pulled to the same potential well. In Eq. 2, the coupling interaction is bidirectional and repulsive in nature. This form of coupling repels the states of the subsystems, thus tending to anti-synchronize the subsystems [34, 37, 38], i.e. for bistable systems, the subsystems are pushed to different potential wells. Thus, Equation 1 describes attractive coupling and Equation 2 describes repulsive coupling. In conjunction with the bias bb, we demonstrate that changing between these two coupling forms offers another degree of control, that allows us to obtain all six fundamental logic operations (c.f. Table 2).

Table 1: Relationship between the two inputs and the output of the fundamental OR, AND, NOR, NAND, XOR and XNOR operations, for the four distinct possible input sets (0,0)(0,0), (0,1)(0,1), (1,0)(1,0) and (1,1)(1,1).
Input Set (I1I_{1},I2I_{2}) OR AND NOR NAND XOR XNOR
(0,0) 0 0 1 1 0 1
(0,1)/(1,0) 1 0 0 1 1 0
(1,1) 1 1 0 0 0 1

The signals I1I_{1} and I2I_{2} encode the stream of binary inputs to be processed by the logic gate, where without loss of generality I1I_{1}(or I2I_{2}) <0<0 corresponds to a 0 (OFF state) and I1I_{1}(or I2I_{2}) >0>0 corresponds to a 11 (ON state). It is important to note that, these input signals are subthreshold i.e., they cannot cause a transition between the potential wells on their own and the transitions are actively facilitated by the noise floor, reminiscent of logical stochastic resonance (LSR). The crucial differentiating factor of this scheme from all previous attempts to achieve noise-aided logic operations is that the output from the system is embedded in the collective state of the coupled system. Specifically, the synchronization state of the two subsystems embeds outputs corresponding to logical operations on the input streams. As a convention, the synchronized state (both systems in same potential well) is taken to encode 0 and the anti-synchronized state (both systems in different potential wells) is taken to encode 11. The binary input output relations represented by the six fundamental 2-input logic gates are detailed in Table 1. With these conventions in place, we now show that this scheme is capable of producing all the six fundamental logic operations in a robust, reliable manner over a large region of parameter space.

III Implementation

Refer to caption
Refer to caption
Figure 2: The top two panels show the input streams I1I_{1} and I2I_{2} that take the value 0.5-0.5 when the logic input is 0 and +0.5+0.5 when the logic input is 11. The bottom three panels depict the timeseries of the state variables x1x_{1} (Blue) and x2x_{2} (Orange) obtained from the numerical simulation of (a) the attractive coupling scheme (c.f. Eq. 1) and (b) the repulsive coupling scheme (c.f. Eq. 2). The expected logic output (green) corresponding to each logic operation are plotted over the timeseries as a visual aid. When x1x_{1} and x2x_{2} are synchronized, the output is interpreted as 0, when x1x_{1} and x2x_{2} are anti-synchronized the outputs are interpreted as 11. Panel 3 shows AND logic operation in the attractive scheme and NAND logic operation in the repulsive scheme for noise strength D=0.2D=0.2, coupling strength c=1.2c=1.2 and bias b=0.5b=0.5. Panel 4 shows NOR logic in the attractive scheme and OR logic in the repulsive scheme for D=0.2D=0.2, c=1.2c=1.2 and b=0.5b=-0.5. Panel 5 shows XNOR operation in the attractive scheme and XOR operation in the repulsive scheme for D=0.25D=0.25, c=0c=0 and b=0b=0.
Refer to caption
Refer to caption
Figure 3: The schematic circuit diagram of (a) the attractively coupled system represented by Eqns. 1 and 3. (b) the repulsively coupled system represented by Eqns. 2 and 3. All component values are indicated in the diagram. The diodes used in the circuit are 1N4148 diodes. Both systems are studied for two values of coupling resistances RC=300ΩR_{C}=300~{}\Omega and R=10KΩR=10~{}K\Omega. The system variables x1x_{1} and x2x_{2} in Eq. 1 and Eq. 2 are proportional to the voltages V1V_{1} and V2V_{2} across the capacitors C1C_{1} and C2C_{2}.
Refer to caption
Refer to caption
Figure 4: Top panels in on both sides are the oscilloscope trace of the logic input signals used to drive the circuits. Left: Oscilloscope traces of the voltages V1V_{1} (yellow) and V2V_{2} (blue) across capacitors C1C_{1} and C2C_{2} of the attractive circuit (c.f. Fig. 3(a)). The expected logical output for each logic gate is presented (in red) as a visual aid. Panel 2 shows AND logic operation obtained for bias b=400mVb=400~{}mV, noise strength D=0.55VD=0.55~{}V and coupling resistance RC=300ΩR_{C}=300~{}\Omega. Panel 3 shows NOR logic for b=400mVb=-400~{}mV, D=0.55VD=0.55~{}V and RC=300ΩR_{C}=300~{}\Omega. In panel 4, XNOR logic is obtained for b=0mVb=0~{}mV, D=0.45VD=0.45~{}V and RC=10KΩR_{C}=10~{}K\Omega. Right: Oscilloscope traces of V1V_{1} (yellow) and V2V_{2} (blue) from the repulsive circuit (c.f. Fig. 3(b)). Panel 2 shows NAND logic operation is obtained for bias b=400mVb=400~{}mV, noise strength D=0.55VD=0.55~{}V and coupling resistance RC=300ΩR_{C}=300~{}\Omega. Panel 3 shows OR logic is obtained for b=400mVb=-400~{}mV, D=0.55VD=0.55~{}V and RC=300ΩR_{C}=300~{}\Omega. In panel 4, XOR logic is obtained for b=0mVb=0~{}mV, D=0.45VD=0.45~{}V and RC=10KΩR_{C}=10~{}K\Omega.

The scheme is first implemented in silico, by numerically simulating Eqns. 1 and 2 using the Euler-Maruyama method. For the bistable potential, we use the simple cubic function F(xi)=4(xi5xi3)F(x_{i})=4(x_{i}-5x_{i}^{3}). The input signals I1I_{1} and I2I_{2} are taken to assume the values 0.5-0.5 for the binary input 0 and +0.5+0.5 for binary input 11. The timetrials of the system variables x1x_{1} and x2x_{2} thus obtained, are depicted in Fig. 2 for various values of noise strength DD, coupling strength cc and bias bb. In this figure, the top two panels, show the input streams I1I_{1}(blue) and I2I_{2}(orange) being fed into the two subsystems and the expected logical output(green) of these inputs for each of the six logic operations are overlaid on the timetrails as a visual aid to perceive the correct logical operation. For the six cases depicted, the timetrails of x1x_{1} and x2x_{2} constantly alternate between the synchronized and anti-synchronized states modulated by the input streams. As defined earlier, the logical output is considered 0 when the two state variables are in the same potential wells (synchronized) and 11 when the state variables are in the opposite potential wells (anti-synchronized). From Fig. 2 it is apparent that robust logic operations are obtained at specific parameter values and coupling forms, for all six types of logic gates.

Next we construct a proof-of-principle electronic implementation of the scheme. Two piecewise-linear bistable circuits were built, using simple passive components and two operational amplifiers(op-amps) as depicted in Fig. 3. The detailed description and characterization of the bistable circuit used can be found in reference [39]. The two bistable units were coupled attractively via a resistor as shown in Fig. 3a or repulsively via inverting amplifiers as shown in Fig. 3b. The inputs and bias to the bistable system are fed through the inverting input of the op-amps, hence the signals and biases were inverted to stay consistent with the scheme description and the numerical exploration. The non-dimensionalized form of the coupled equation governing the circuits in Fig. 3a (attractively coupled circuit) and Fig. 3b (repulsively coupled circuit) assumes the form described in Eq.1 and Eq.2 respectively. Where, F of bistable unit is given by the piece-wise linear function,

F(xi)={(xi+1)xi<0.5xi0.5xi0.5(xi1)xi>0.5F(x_{i})=\left\{\begin{array}[]{ll}-(x_{i}+1)&\quad x_{i}<-0.5\\ x_{i}&\quad-0.5\leq x_{i}\leq 0.5\\ -(x_{i}-1)&\quad x_{i}>0.5\end{array}\right. (3)

where the state variables x1x_{1} and x2x_{2} are proportional to the voltages V1V_{1} and V2V_{2} across the capacitors C1C_{1} and C2C_{2}. The timetrails of the experimental systems were obtained using a Tektronics 2104B Digital Storage Oscilloscope. A high speed data acquisition device (Measurement Computing USB-1616HS) was used to both generate signals (I1I_{1}, I2I_{2}, η1\eta_{1}, η2\eta_{2}) and collect high throughput voltage data (V1V_{1}, V2V_{2}) for further analysis. All signal generation and data collection were done at a rate of 2×1042\times 10^{4} samples per second. All specific component values used in the construction of the circuit are indicated in the circuit schematic (c.f. Fig. 3).

In Fig. 4 oscilloscope trails of V1V_{1}(yellow) and V2V_{2}(blue) obtained from the attractive circuit (left) and repulsive circuit (right) are presented. The top panel shows the input streams (I1I_{1} and I2I_{2}) driving the coupled system while the remaining six panels clearly show robust logic response akin to the behaviour observed in Fig. 2. The expected logic output of the input streams are again overlaid(in red) as a visual aid. The synchronized segments faithfully map to the output 0 while the anti-synchronized segments map to 11 for all six logical operations.

The attractive and repulsive coupling schemes yield three logic gates each. The attractive scheme yields AND, NOR and XNOR gates, while the repulsive scheme yields NAND, OR and XOR gate. Notice that the two schemes yield complementary logic operations for the same set of parameter values, suggesting that these two coupling schemes are symmetric counterparts of one another. Within a specific coupling form, the value of the constant bias bb (positive, negative or zero) determines the logical operation obtained from the system. The bias ranges and the coupling forms where each of the six logic operations occur, are detailed in Table 2.

Table 2: Logic operations obtained for the various coupling schemes and bias bb values.
Positive bias Negative bias Zero bias
Attractive coupling AND NOR XNOR
Repulsive coupling NAND OR XOR
Refer to caption
Figure 5: Probability of obtaining reliable logic operations P(logic)P(logic) for all six fundamental logic operations (P(AND), P(NOR), etc.) is plotted as a function of noise strength DD and coupling strength cc. The plots are made for both attractive and repulsive coupling schemes at three specific values of bias bb, depicting broad regions in parameter space where all six logic operations are consistently obtained.

IV Measure of reliability

Robust operation of the scheme, has been demonstrated both using simulation and experiment for a stream of inputs, with specific values of system parameters. The performance of this system is now quantified with a large number of input sets over a significant section of the parameter space. To do this, the system is subject to a large number of (I1I2I_{1}-I_{2}) sets, divided into runs where each run consists of a permutation of the four input sets (0,0)(0,0), (0,1)(0,1), (1,0)(1,0), (1,1)(1,1) and the response of the systems to these inputs are recorded. A run is considered successful, only if the system produced the correct logical output (corresponding to each truth table c.f. 1) throughout the run-time of each signal pulse for all four input sets. The probability of obtaining a specific logic operation P(logic)P(logic) is then defined as the ratio of successful runs to the total number of runs sampled. A small transient time amounting to one tenth the duration of each input pulse is allowed for the system to respond to each new input set. The P(logic)P(logic) corresponding to specific logic operations are denoted as P(AND), P(XOR), etc. This measure is then obtained for a large range of parameter values to ascertain the prevalence of reliable logic operations.

In Fig. 5, numerically obtained probabilities of obtaining different logic gates are plotted for a range of noise strengths DD and coupling strengths cc. In these simulations the P(logic)P(logic) was determined by subjecting the system to 100 runs as described earlier, for each combination of parameter values. Reliable logic operations were seen to occur in large sections of parameter space where the P(logic)P(logic) value tends to 11 (bright yellow regions in Fig. 5). The probability plots were made for each of the six logic operations at the corresponding coupling scheme and bias bb values mentioned in Table 2. Thus, we see that for a window of optimal noise strengths and coupling strengths, all logic operations can be obtained. The coupling form and bias bb act as the control to morph the coupled system from one logic operation to the other. Note that contrary to coupling induced LSR[35], where robust logic occurs for all values above a critical coupling strength, we find an optimal range of coupling strengths yield reliable logic operations.

Refer to caption
Refer to caption
Figure 6: Probability of logic P(logic) calculated from experimentally obtained voltage data, plotted as a function of noise strength DD. Top: For coupling resistance RC=300ΩR_{C}=300~{}\Omega, P(AND) and P(NAND) logic plotted for bias b=400mVb=400~{}mV. P(NOR) and P(OR) are plotted for b=400mVb=-400~{}mV. Bottom: For coupling resistance RC=10KΩR_{C}=10~{}K\Omega, P(XNOR) and P(XOR) logic plotted for bias b=0mVb=0~{}mV.

To further strengthen this assertion, we also obtain the same stringent measure of reliable operations P(logic)P(logic), from the experimental implementation of the scheme. This was made possible by interfacing and voltage data analysis through the high speed DAQ. The experimental circuits were driven by an input stream containing 100 runs with a permutation of the four input sets (0,0)(0,0), (0,1)(0,1), (1,0)(1,0), (1,1)(1,1) in every run. Again, P(logic)P(logic) is obtained from the ratio of successful runs to total number of runs sampled. The inputs were fed into the system at 10Hz (10 binary inputs per second) and again one tenth of that pulse width was allowed as transient and was not included in the computation of P(logic)P(logic). Experimentally, two values of coupling resistances were studied, R=300ΩR=300~{}\Omega and R=10KΩR=10~{}K\Omega. The distribution of P(logic)P(logic) thus obtained for all the six logic gates are shown in Fig. 6. A clear maximization of reliable logic operations, occurs for a broad window of noise strength DD for all the logic operations. Surprisingly, contrary to the numerical exploration experimentally we find a more robust and broad region of XOR and XNOR operations. Note that the XOR/XNOR operations were the hardest to implement in the past, and their realization necessitated the use of more complicated triple well potentials [40], where specific output definitions were assigned for each logic operation. Here, the multiple gates emerge from the collective dynamics of the coupled system and the binary outputs are inferred from their synchronization state.

V Conclusion

A novel scheme to make reconfigurable noise aided logic gates, based on synchronization of coupled nonlinear systems, was introduced. This scheme was implemented both through numerical simulations and electronic experiments. The robustness of the logical operations and the reconfigurability of the scheme was elucidated for both attractive and repulsive coupling. A quantitative measure of performance was used to characterize the region in parameter space of coupling strength and noise strength where reliable logic operations can be obtained. This was done both in simulations and in circuit experiments through live collection and processing of high throughput voltage data. Importantly, all six fundamental logic operations were reliably obtained, using a bias and the coupling form to morph between the logic functionalities. So our results suggest the potential of exploiting synchronization, arising from the interplay of noise and the nature of coupling, to implement flexible logic gates in the presence of a noise floor.

Acknowledgements.
The authors gratefully acknowledge the discussions with Dr. K. Murali, Anna University, Chennai, which helped shape the experimental implementation of the repulsive coupling scheme used in this work.

References

  • Sinha and Ditto [1998] S. Sinha and W. L. Ditto, Physical Review Letters 81, 2156 (1998).
  • Crutchfield et al. [2010] J. Crutchfield, W. Ditto, and S. Sinha, Chaos: An Interdisciplinary Journal of Nonlinear Science 20, 037101 (2010).
  • Murali and Sinha [2007] K. Murali and S. Sinha, Physical Review E (Rapid Communications) 75, 025201(R) (2007).
  • Murali et al. [2009a] K. Murali, S. Sinha, W. L. Ditto, and A. R. Bulsara, Physical Review Letters 102, 104101 (2009a).
  • Sinha et al. [2009] S. Sinha, J. Cruz, T. Buhse, and P. Parmananda, Europhysics Letters 86, 60003 (2009).
  • Bulsara et al. [2010] A. R. Bulsara, A. Dari, W. L. Ditto, K. Murali, and S. Sinha, Chemical Physics 375, 424 (2010).
  • Gupta et al. [2011] A. Gupta, A. Sohane, V. Kohar, K. Murali, and S. Sinha, Physical Review E 84, 055201(R) (2011).
  • Zhang et al. [2010] L. Zhang, A. Song, and J. He, Physical Review E 82, 051106 (2010).
  • Zhang et al. [2012] H. Zhang, Y. Xu, W. Xu, and X. Li, Chaos: An Interdisciplinary Journal of Nonlinear Science 22, 043130 (2012).
  • Kohar and Sinha [2012] V. Kohar and S. Sinha, Physics Letters A 376, 957 (2012).
  • Gui et al. [2020] R. Gui, H. Zhang, G. Cheng, and Y. Yao, Chaos: An Interdisciplinary Journal of Nonlinear Science 30, 023119 (2020).
  • Kohar et al. [2014] V. Kohar, K. Murali, and S. Sinha, Communications in Nonlinear Science and Numerical Simulation 19, 2866 (2014).
  • Das and Kantz [2019] M. Das and H. Kantz, Phys. Rev. E 100, 032108 (2019).
  • Mingjie et al [2020] H. Mingjie et al, The European Physical Journal Plus 135, 747 (2020).
  • Cheng et al. [2020] G. Cheng, W. Liu, R. Gui, and Y. Yao, Chaos, Solitons and Fractals 131, 109514 (2020).
  • Murali et al. [2009b] K. Murali, I. Rajamohamed, S. Sinha, W. L. Ditto, and A. R. Bulsara, Applied Physics Letters 95, 194102 (2009b).
  • Guerra et al. [2010] D. N. Guerra, A. R. Bulsara, W. L. Ditto, S. Sinha, K. Murali, and P. Mohanty, Nano Letters 10, 1168 (2010).
  • Jaimes-Reátegui et al [2014] R. Jaimes-Reátegui et al, Eur. Phys. J. Spec. Top. 223, 2837 (2014).
  • Hellen et al. [2013] E. H. Hellen, S. K. Dana, J. Kurths, E. Kehler, and S. Sinha, Plos one 8, e76032 (2013).
  • Ando et al. [2011] H. Ando, S. Sinha, R. Storni, and K. Aihara, Europhysics Letters 93, 50001 (2011).
  • Dari et al. [2011] A. Dari, B. Kia, A. R. Bulsara, and W. L. Ditto, Chaos: An Interdisciplinary Journal of Nonlinear Science 21, 047521 (2011).
  • Sharma et al. [2014] A. Sharma, V. Kohar, M. D. Shrimali, and S. Sinha, Nonlinear Dynamics 76, 431 (2014).
  • Lei et al. [2018] Z. Lei, W. Zheng, and A. Song, Chaos: An Interdisciplinary Journal of Nonlinear Science 28, 043117 (2018).
  • Wu et al. [2017] J. Wu, Y. Xu, H. Wang, and J. Kurths, Chaos: An Interdisciplinary Journal of Nonlinear Science 27, 063105 (2017).
  • Wang and Song [2016] N. Wang and A. Song, IEEE transactions on neural networks and learning systems 27, 2736 (2016).
  • Singh and Sinha [2011] K. P. Singh and S. Sinha, Physical Review E 83, 046219 (2011).
  • Zamora-Munt and Masoller [2010] J. Zamora-Munt and C. Masoller, Optics Express 18, 16418 (2010).
  • Pfeffer et al. [2015] P. Pfeffer, F. Hartmann, S. Höfling, M. Kamp, and L. Worschech, Physical Review Applied 4, 014011 (2015).
  • Murali et al. [2018] K. Murali, S. Sinha, V. Kohar, B. Kia, and W. L. Ditto, Plos One 13, e0209037 (2018).
  • Murali et al. [2021] K. Murali, S. Sinha, V. Kohar, and W. L. Ditto, The European Physical Journal Special Topics , 1 (2021).
  • Aravindh et al. [2018] M. S. Aravindh, A. Venkatesan, and M. Lakshmanan, Physical Review E 97, 052212 (2018).
  • Sathish Aravindh et al. [2020] M. Sathish Aravindh, A. Venkatesan, and M. Lakshmanan, Chaos: An Interdisciplinary Journal of Nonlinear Science 30, 093137 (2020).
  • Aravind et al. [2020] V. M. Aravind, K. Murali, and S. Sinha, in Nonlinear dynamics and control (Springer, 2020) pp. 325–334.
  • Aravind et al. [2021] M. Aravind, S. Sinha, and P. Parmananda, Chaos: An Interdisciplinary Journal of Nonlinear Science 31, 061106 (2021).
  • Aravind et al. [2018] M. Aravind, K. Murali, and S. Sinha, Physics Letters A 382, 1581 (2018).
  • Stankovski et al. [2017] T. Stankovski, T. Pereira, P. V. McClintock, and A. Stefanovska, Reviews of Modern Physics 89, 045001 (2017).
  • Hens et al. [2013] C. Hens, O. I. Olusola, P. Pal, and S. K. Dana, Physical Review E 88, 034902 (2013).
  • Dixit and Shrimali [2020] S. Dixit and M. D. Shrimali, Chaos: An Interdisciplinary Journal of Nonlinear Science 30, 033114 (2020).
  • Adomaitienė et al. [2008] E. Adomaitienė, A. V. Tamaševičius, G. Mykolaitis, S. Bumelienė, and E. Lindberg, Nonlinear Analysis: Modelling and Control 13, 241 (2008).
  • Storni et al. [2012] R. Storni, H. Ando, K. Aihara, K. Murali, and S. Sinha, Physics Letters A 376, 930 (2012).