#### Abstract

Exploring and investigating new chaotic systems is a popular topic in nonlinear science. Although numerous chaotic systems have been introduced in the literature, few of them focus on torus-chaotic system. The aim of our short work is to widen the current knowledge of torus chaos. In this paper, a new torus-chaotic system is proposed, which has one positive Lyapunov exponent, two zero Lyapunov exponents, and two negative Lyapunov exponents. The dynamic behavior is investigated by Lyapunov exponents, bifurcations, and stability. The analysis shows that this system has an interesting route leading to chaos. Furthermore, the pseudorandom properties of output sequence are well studied and a random number generator algorithm is proposed, which has the potential of being used in several cyber security systems such as the verification code, secure QR code, and some secure communication protocols.

#### 1. Introduction

Over the last four decades, chaos theory and its applications have been a central topic of nonlinear science [1]. As it is well known, chaos is associated with stochasticity, complex, irregular motion, and so on. It has some peculiar properties such as ergodicity, highly initial value sensitivity, nonperiodicity, and long-term unpredictability. These pseudorandom features lead chaotic systems to enormous cryptographic applications like random number generator, image encryption, and secure communication, to name a few [2, 3].

There are two main paths for studying chaotic systems. The first path involves studying novel dynamic behaviors of chaotic systems. For example, Faghani et al. proposed some three-dimensional chaotic systems based on the general form of quadratic jerk systems. They analyzed the equilibrium, eigenvalue, and Lyapunov exponents of the proposed systems and proved that there are 12 simple chaotic systems with identical eigenvalues. This is the first time that different chaotic systems with identical eigenvalues are described [4]. Qi and Liang investigated the angular momentum, torques, and Casimir energy of Qi system. They discovered different types of torques of the chaotic attractors [5]. Zhu et al. analyzed a sine-cosine compound function system. They proved that the proposed system satisfies Devaney’s chaos definition and it has better unpredictability and more complex chaotic behaviors than classical sine map and Chebyshev map [6].

The second path proposes engineering applications of chaotic systems. For example, Liu et al. applied a simple chaotic system to differential chaos shift keying (DCSK) scheme. They demonstrated that the proposed spread spectrum communication scheme can reduce the bit error rate in the presence of noise [7]. Liu et al. proposed a chaos-based color image encryption scheme, which uses one-time keys from environmental noise [8]. Zhu et al. proposed a 2D compound homogeneous hyperchaotic system (CHHCS). They demonstrated that the proposed scheme provides sufficient security against known attacks [9]. Hua et al. proposed a cosine-transform-based chaotic system (CTBCS). They proved that CTBCS can produce chaotic maps with complex dynamical behaviors. Therefore, the further proposed image encryption scheme provides a higher level of security than several advanced image encryption schemes [10].

In fact, the secure communication technology and data encryption technology based on chaos are essentially implemented using the pseudorandom nature of the chaotic systems. Therefore, the performance of the pseudorandom sequence generated by the chaotic system greatly affects the scores of these schemes. In recent years, a large number of studies have been conducted on random number generators (RNG). Singh et al. proposed a new chaotic oscillator by using a single op amp, two capacitors, one resistor, one inductor, and memristive diode bridge cascaded with an inductor. Based on this system, they designed a pseudorandom number generator and implemented it using Raspberry Pi 3 [11]. Alcin et al. proposed a novel type of high‐speed true random number generator based on chaos and ANN and implemented in a Xilinx field-programmable gate array (FPGA) chip. They demonstrated that the implemented TRNG offers throughput up to 115.794 Mbps. Besides, the generated random numbers have been tested with the FIPS 140-1 and NIST 800.22 test suites [12]. Wang and Cheng proposed a pseudorandom number generator based on a modified logistic chaotic system. They proved that the pseudorandom sequence generated by this method has perfect randomness and cryptographic properties and can pass the statistical tests [13].

As it is well known, sensitive dependence on initial conditions, which can be quantified by Lyapunov exponent (LE), is the most important feature of chaos. For example, a three-dimensional bounded dynamic system can exhibit chaotic behavior if and only if its LEs are . When all the LEs are negative, the orbit starting from any point will eventually converge to a stable point. From another perspective, LEs can be regarded as the notion of divergence of orbits or degradation of information, the presence of positive LE suggests that the original information is lost, while zero LEs are equivalent to some form of invariant of two orbits or preserving the information of dynamic systems. Therefore, LEs can be a practical way to characterize and classify chaos [14]. According to Witold Kinser, a continuous-time system may have the following distinct dynamic behaviors [15]:(1)Limit cycle, whose LEs are (2)Torus, whose LEs are (3)Chaos, whose LEs are (4)Hyperchaos, whose LEs are (5)Torus chaos, whose LEs are

Compared with chaotic systems, hyperchaotic systems have more complex behavior. With more than one positive Lyapunov exponents, the dynamics of hyperchaotic systems are expanded in more than one direction giving rise to a more complex attractor. As Li demonstrated, hyperchaotic signals have great instability in two dimensions and relatively more complicated dynamics, and because of this, hyperchaotic circuits can improve security which are useful in secure communications and radar engineering [16].

To the best of our knowledge, torus chaos has seldom been discussed among numerous proposed literatures. In 2009, Yang et al. proposed a 4D hyperchaotic system by adding a linear controller to the second equation of the 3D Lorenz system. They have shown that the proposed system can exhibit torus-chaos behavior when the parameters satisfied some special values [17]. In 2011, his team reported a 5D hyperchaotic system based on modified generalized Lorenz system. This system can also exhibit torus-chaos behavior [18]. In 2017, Singh and Roy proposed a simplest 4D dissipative autonomous chaotic system which is derived from Rossler type-IV 3D chaotic system. They proved that this system has chaotic 2-torus nature of Lyapunov exponents when bifurcation parameter satisfies some critical value [19]. However, the dynamic behavior of torus-chaotic system has not been well studied. This motivated us to further study torus-chaotic systems [20].

In this paper, we proposed a new torus-chaotic system. The dynamic behavior is investigated by Lyapunov exponents, bifurcations, and stability. The analysis shows that this system has an interesting route leading to chaos. A period-doubling behavior is observed in one subspace, but in other subspaces, there is no period-doubling behavior. It is like the saddle point: the system is stable in one direction but unstable in the other direction. Furthermore, the pseudorandom properties of output sequences are well studied. It is proved that the torus-chaotic system has worse pseudorandom properties than hyperchaotic system but has better pseudorandom properties than chaotic system. Based on its pseudorandom properties, a random number generator algorithm is proposed, which has the potential of being used in several cyber security systems such as the verification code, secure QR code, and some secure communication protocols.

The rest of this paper is organized as follows. In Section 2, the novel torus-chaotic system is introduced. The analysis of dissipativity, equilibrium, stability, and Lyapunov vectors was presented in Section 3. Route to chaos is discussed in Section 4. Section 5 is the pseudorandom properties of output sequence and the random number generator algorithm. Conclusions are finally drawn in Section 6.

#### 2. A Novel Torus-Chaotic System

In 2018, Liu et al. introduced a general method to construct nth-order chaotic systems with hyperbolic sine nonlinearity [21]. When *n* = 5, there is a fifth-order chaotic system which is described by

By introducing a nonlinear feedback controller to system equation (1), the following system which can exhibit torus-chaos behavior is obtained:where , and are the state variables. , , and , which have been chosen to facilitate circuit implementation using diodes in future. *c* is considered as the bifurcation parameter.

To numerically calculate the phase space plots, Lyapunov exponents, and Poincaré sections, a fourth-order Runge–Kutta integrator with fixed step of 0.001 is utilized. To avoid transient response, the initial 1e6 iterated data are omitted, and the following 2e6 data are used in plotting phase space plots and Poincaré sections, while the following 1e9 data are used in calculating some critical Lyapunov exponents.

When *c* *=* 1.55 and the initial conditions are set to be and , the system has two coexisting attractors as shown in Figures 1(a) and 1(b). When , the Lyapunov exponents converged to , which suggests that the system exhibits two-torus-chaotic behavior [22].

**(a)**

**(b)**

In order to catch the trajectory of the attractor, the Poincaré sections of and plane are plotted in Figure 2, of which *c* = 1 and initial conditions is set to be . From Figure 2, the system is exhibiting chaotic behavior.

**(a)**

**(b)**

#### 3. Dynamical Behaviors of the Novel System

##### 3.1. Dissipativity of the Novel System

For system equation (2), it is noticed that

When *c* > −1, the system is dissipative with an exponential contraction rate: . In this situation, this system is an open system which has dissipative structures [23], all system orbits are ultimately confined to some subset of zero volumes as , and the asymptotic motion settles on some attractors [24].

##### 3.2. Equilibrium and Stability

Suppose that the system , possesses an equilibrium point . That is, . Then is a solution for all of the time. Therefore, the equilibrium point can be solved by .

For this system, the equilibrium points are given by

Therefore, it has three equilibrium points, which are , , and .

The Jacobian matrix at the equilibrium is

The characteristic equation of (5) iswhere . When *c* = 1, (6) exhibits two-torus-chaos behavior with LEs of . The eigenvalues of corresponding Jacobian matrix are

Since cosh(*x*) is an even function, the eigenvalues of equilibrium and are the same. Noticing that and of , , and are all pairs of complex conjugate eigenvalues with positive real parts, is all real eigenvalue and less than zero and and are all pairs of complex conjugate eigenvalues with negative real parts. Therefore, the three equilibrium points are all saddle points which have 1D knot manifold, 2D stable focus manifold, and 2D unstable focus manifold [25].

#### 4. Visualization of Lyapunov Vectors

Lyapunov vectors, which have generally received less attention in the literature than Lyapunov exponents, are the time-dependent spatial structures associated with the corresponding Lyapunov exponents, which are in turn the asymptotic exponential growth rates of linear disturbances in general time-dependent flows. With the Lyapunov vectors, expanding/shrinking subspaces can be identified. In this section, we present the visualization of Lyapunov vectors in this torus-chaotic system, so that expanding/shrinking subspaces can be observed.

We utilize quiver function of Matlab to visualize Lyapunov vectors. As it is known, this system has five LEs, and each LE has five vector components. In this experiment, all initial conditions are set to be (0.1, 0.1, 0.1, 0.1, 0.1). The positive Lyapunov vector, zero Lyapunov vectors, and negative Lyapunov vector are shown in Figures 3–5.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

From Figures 3–8, the following can be observed:(1)All Lyapunov vectors rotate greatly along the orbit, which indicate that the expanding/shrinking subspaces changed along with the orbit.(2)Figures 4 and 5 describe two zero local Lyapunov vectors but none of them follows the tangent direction of the trajectory, which indicates that there is an expanding/shrinking mechanism along with the direction of the trajectory. When time approaches infinity, the average effect of expanding/shrinking of these two directions is zero.(3)Two close adjacent trajectories have completely different or even opposite Lyapunov vectors, which may suggest that chaotic system has multiple stretching and folding mechanisms in infinitesimal neighborhoods.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

#### 5. Routes to Chaos

There are five important routes leading to chaos that have been identified and confirmed by numerical and physical experiments. They are period-doubling cascade, Ruelle–Takens–Newhouse, intermittency, crisis, and quasiperiodic routes.

The period-doubling cascade, which is also known as Feigenbaum’s route to chaos, is the most common behavior which exists in numerous chaotic systems. In this route to chaos, a fixed point loses stability to an attracting period-2 limit cycle as some parameters passed critical value. At later parameter value, this limit cycle also becomes unstable and turns into period-4 limit cycle. The new limit cycle then goes through a cascade of period doublings until chaos [26].

Ruelle–Takens–Newhouse has proved that if the attractor of a system has three independent frequencies, then a small perturbation of this system has a hyperbolic strange attractor, a Plykin attractor. This scenario is called Ruelle–Takens–Newhouse route to chaos [27].

Intermittency denotes a type of behavior where the dynamics vary chaotically between two different phases of motion. One of these phases is regular (close to stationary, periodic, or quasiperiodic motion) and is called laminar phase. The laminar phase is interrupted by turbulent bursts, which corresponds to some irregular phases of motion. In chaotic systems, there exist different types of intermittencies. Pomeau and Manneville described three different types of intermittency. Each type of intermittency is related to a different kind of bifurcation, which are saddle-node bifurcation, subcritical Hopf bifurcation, and inverse period-doubling bifurcation [28].

In the crisis route, a chaotic attractor is suddenly created to replace a nonattracting chaotic saddle as the parameter passes through the crisis value. Common types of crisis are boundary crisis, interior crisis, and symmetry restoring (breaking) crisis. These routes are due to the collision of the chaotic attractor with an unstable invariant set, typically an unstable periodic orbit (equivalently, a collision between a chaotic attractor and a coexisting unstable fixed point or periodic orbits) [29].

Quasiperiodic routes are also known as torus breakdown. In this mechanism, the torus (quasiperiodic orbit) loses its smoothness and gives birth to chaos as some parameter passes some critical value. This scenario has been observed in many natural and man-made systems including the Rayleigh–Benard convection [30], the Taylor–Couette flow [31], and electrical circuits [32].

In 2017, Singh and Roy proposed a 4D dissipative autonomous chaotic system which can exhibit 2-torus, 3-torus, and torus-chaos behavior. The numerical calculation shows that there are a period-doubling route and an intermittency route to chaos [19].

To study the route from torus to torus-chaos in this system, the LEs spectrum, Kaplan–Yorke dimension spectrum, and bifurcations [33] of (2) are plotted as the coefficient *c* is varied over the range *c* ∈ [0.3*,* 2]. Note that Figure 9(c) shows the value of every time reaches a local maximum or equivalently whenever crosses zero in the downward direction [34]. Figure 9(d) shows the value of every time reaches a local maximum or equivalently whenever crosses zero in the downward direction. Figure 9(e) shows the value of every time reaches a local maximum or equivalently whenever crosses zero in the downward direction. Figure 9(f) shows the value of every time reaches a local maximum or equivalently whenever crosses zero in the downward direction. Figure 9(g) shows the value of every time . Figure 9(h) shows the value of every time . According to the LEs spectrum, one can locate the region of torus chaos, and then the dimension of a torus-chaotic attractor determined the Kaplan–Yorke dimension spectrum. Since the system cannot exhibit hyperchaotic behavior within the parameters, the dimension of a torus-chaotic attractor is between 3.1 and 3.5, the dimension of 2-torus attractor is 2, and 3-torus attractor is 3.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

**(g)**

**(h)**

Figure 9 suggests that this system has an interesting route leading to chaos:(1)When *c* ∈ [0.3*,* 0.4639], there exists a period-doubling behavior along with and subspace. However, the system shows torus behavior along with and subspace. It is like saddle point; the system is stable in one direction but unstable in the other direction.(2)When *c* ∈ [0.4640*,* 0.5574], the system exhibits two-torus-chaos behavior except for some 2-torus windows. When the parameter passed *c* = 0.4639 to *c* = 0.4640, two-torus chaos is born by replacing the 2-torus behavior. The LEs at these two critical values are (LE1, LE2, LE3, LE4, LE5) = (0, 0, −0.01, −0.57, −0.88) for *c* = 0.4639 and (LE1, LE2, LE3, LE4, LE5) = (0.02, 0, 0, −0.60, −0.88) for *c* = 0.4640. This may be caused by the period-doubling route along with and subspace**.**(3)When *c* ∈ [0.5575, 0.5901], the system exhibits 2-torus behavior.(4)When *c* ∈ [0.5902, 1.5575], the system exhibits two-torus-chaos behavior except for 2-torus windows. The route leading to chaos is similar to point 3.(5)When *c* ∈ [1.5575*,* 2], the system exhibits 2-torus behavior, except for some 3-torus windows like *c* = 1.6157.

#### 6. Pseudorandom Properties

Since chaos is sensitive to initial conditions, it is very suitable to apply them to the purpose of random number generator. But before that, one should notice that the probability density distributions (PDDs) of chaotic systems are not uniform distribution [35]. Figure 10 is the PDDs of , which indicates that it has physical characteristic. Thus, it cannot resist side channel attack.

In fact, all raw data generated by chaotic or hyperchaotic system cannot be used for cryptographic purposes. We compared pseudorandom properties of chaotic system, torus-chaotic system, and hyperchaotic system described in [36] and [18]. The results are shown in Tables 1∼3.

The testing results indicate that hyperchaotic system has the best pseudorandom properties; it passed frequency, linear complexity, and cumulative sums tests. Chaotic system has worst pseudorandom properties which passed none of the tests. Torus-chaotic system passed overlapping template which means that its pseudorandom properties are better than chaotic system. However, all of these raw data cannot be used for cryptographic purposes because the physical characteristics have not been removed.

There are many ways to remove the physical characteristics of the output sequence; for example, in [11], the authors introduced an LSB-16 method to generate random number. In this manuscript, we utilize the following decorrelation operation:

We tested the data of torus-chaotic system with decorrelation operation. The result is shown in Table 4.

From the testing results shown in Table 4, all output sequences passed fifteen kinds of random test of NIST 800-22, which indicated that the proposed method can be used as a part of some cyber security systems such as the verification code, secure QR code, and some secure communication protocols.

#### 7. Conclusion

In this paper, we proposed a novel torus-chaotic system. The dynamic behavior is investigated by Lyapunov exponents, bifurcations, and stability. The analysis shows that this system has an interesting route leading to chaos. A period-doubling behavior is observed in one subspace, but in other subspaces, there is no period-doubling behavior. It is like saddle point: the system is stable in one direction but unstable in the other direction. Furthermore, the pseudorandom properties of output sequence are well studied. It is proved that the torus-chaotic system has worse pseudorandom properties than hyperchaotic system but has better pseudorandom properties than chaotic system. Based on its pseudorandom properties, a random number generator algorithm is proposed, which has the potential of being used in several cyber security systems such as the verification code, secure QR code, and some secure communication protocols.

#### Data Availability

All data are generated by the authors. One can repeat their results by the description provided in this manuscript.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

The authors would like to thank Professor Julien Clinton Sprott for helpful discussion. Jizhao Liu had received research grants from Sun Yat-sen University. This study was supported by the Fundamental Research Funds for the Central Universities (no. 19lgpy230).