Information Transfer Among the Components in Multi-Dimensional Complex Dynamical Systems

In this paper, a rigorous formalism of information transfer within a multi-dimensional deterministic dynamic system is established for both continuous flows and discrete mappings. The underlying mechanism is derived from entropy change and transfer during the evolutions of multiple components. While this work is mainly focused on three-dimensional systems, the analysis of information transfer among state variables can be generalized to high-dimensional systems. Explicit formulas are given and verified in the classical Lorenz and Chua’s systems. The uncertainty of information transfer is quantified for all variables, with which a dynamic sensitivity analysis could be performed statistically as an additional benefit. The generalized formalisms can be applied to study dynamical behaviors as well as asymptotic dynamics of the system. The simulation results can help to reveal some underlying information for understanding the system better, which can be used for prediction and control in many diverse fields.


Introduction
Uncertainty quantification in complex dynamical systems is an important topic in prediction models. By integrating information-theoretic methods to investigate potential physics and measure indices, the uncertainty can be quantified better in ensemble practical predictions of complex dynamical systems. For instance, one of the important motivations is the couplings among variables of dynamical systems generating information at a nonzero rate [1], which produces information exchange [2]. Entropy can be used to quantitatively describe production, gathering, exchange and transfer of information [3]. Information transfer analysis can be used to detect asymmetry in the interactions of subsystems [1,4]. The emergent phenomena cannot be simply derived or solely predicted from the knowledge of the structure or from the interactions among individual elements in complex systems [5]. The dynamics of information transportation plays a critical role in complex systems, resulting in the system prediction [6,7], controls of a system [8,9] and causal analysis [10,11]. It emphasizes further understanding and investigating information transportation in complex dynamical systems. It has been applied to quantify nonlinear interactions based on the information transfer by several underlying efficient estimation strategies in complex dynamical systems [12][13][14]. Simple examples are used to illustrate various complex phenomena. The formalisms about information transfer are mostly based on two time series [1,[15][16][17].
Recently, a new approach on information flow between the components of two-dimensional (2D) systems was adapted by Liang and Kleeman [6], which can be used to deal with the change of the uncertainty of one component given by the other component. This idea is based on specific interactions between two components in complex dynamical systems. For a system with dynamics given, a measure of information transfer can be rigorously formulated (referred as LK2005 formalism henceforth in [6]).
In the forms of continuous flows and discrete mappings, the information flow has been analyzed using the Liouville equations [18] and the Frobenius-Perron operators [18]. These two equations are the evolution equations of the joint probability distributions, respectively. The present formalism is consistent with the transfer entropy of Schreiber [1] in both transfer asymmetry and quantification. A variety of generalizations and applications of the work in Reference [4] are developed in [19][20][21][22][23][24][25]. Majda and Harlim [26] applied the strategy to study subspaces of complex dynamical systems. For 2D systems, Liang and Kleeman discovered a concise law on the entropy evolution of deterministic autonomous systems and obtained the time rate of information flow from one component to the other [6]. Until now, the 2D formalism has been extended to some dynamical systems in different forms and scales with successful applications between two variables [23,25]. In the light of these applications, by thoroughly describing the statistical behavior of a system, this rigorous LK2005 formalism has yielded remarkable results [3].
However, the uncertainty of many real-world systems needs to be quantified among the variables for revealing the nonlinear relationships, so as to better understand the intrinsic mechanism and predict the forthcoming states of the systems [27]. Besides, many physical systems are affected by the interactions between multiple components in diverse fields [28]. For example, sensitivity analysis of an aircraft system with respect to design variables, parameters and uncertainty factors can be used to estimate the effects on the objective function or constraint function. The uncertainty analysis and sensitivity analysis (UASA) process is one of the key steps for determining the optimal search direction and guiding the design and decision-making, which aims at predicting complex computer models by quantifying the sensitivity information of the coupling variables. It can be offered to quick guide of determining design parameters which lead to high performance aircraft designs. Some preceding tools [29,30] related to sensitivity analysis are applicable for low-dimensional static problems and an urgent problem of high dimensionality arises when outputting variables of numerical models with spatially and temporally need to be solved [31]. The rigorous formalism of information flow has the potential to revolutionize the ability to analyze and measure uncertainty and sensitivity information in dynamical systems.
Hence, considering realistic applications, we generalize the LK2005 formalism to several variables of multi-dimensional dynamical systems in this paper. More precisely, we extend the results in [6,25] to the information flow between groups of components, rather than individual components. We aim to demonstrate that the formalism is feasible among several variables in arbitrary multi-dimensional dynamical systems when dynamics is fully known. In addition, the generalized formalisms can be reduced to two-dimensional formalisms as a special case. We also highlight the relationship between the LK2005 formalism and our generalized formalisms. Two applications are proposed with the classical Lorenz system and Chua's system as validations of our formalisms. Compared with the LK2005 formalism and the transfer mutual information method [32], the generalized formalisms are beneficial for revealing more information among variables. It can better explore the complexity of evolution and intrinsic regularity of multi-dimensional dynamical systems. Meanwhile, it can provide a simple and versatile method to analyze sensitivity in dynamical models. These generalized formulas enable one to understand the relationship between information transfer and the behavior of a system. It can be used to perform sensitivity analysis as a measure in multi-dimensional complex dynamical systems. Therefore, the generalized formalisms have much wider applications and are significant to investigate real-world problems.
The structure of this paper is as follows: Section 2 recalls a systematic introduction of the theories and the formalisms about information flow in 2D systems; In Section 3, the formalisms are generalized to adapt to multi-dimensional complex dynamical system components based on the LK2005 formalism. Details on the derivations of the formalisms and the related properties are demonstrated; Section 4 gives a description about the formalisms with multi-dimensional applications; the summary of this paper is given in Section 5.

Continuous Flows
For 2D continuous and deterministic autonomous systems with fully known dynamics, where x 2 ) for any i = 1, 2 is known as the flow vector and x = (x 1 , x 2 ) ∈ Ω = Ω 1 × Ω 2 . A stochastic process X = (X 1 , X 2 ) ∈ Ω with joint probability density ρ(x 1 , x 2 , t) at time t is the random variables corresponding to the sample values (x 1 , x 2 ). For convenience, we will use the notation ρ or ρ(x 1 , x 2 ) instead of the notation ρ(x 1 , x 2 , t) throughout Section 2, including the same expression at multi-dimensional cases in Section 3. In addition, the integral domain is the whole sample space Ω, except where noted. The probability density ρ associated with Equation (1) satisfies the Liouville equation [18]: ∂ρ ∂t The rate of change of joint entropy of X 1 and X 2 , where E means the mathematical expectation with respect to ρ and E(∇ · F) = Ω ρ(∇ · F)dx 1 dx 2 . That is to say, when a system evolves with time, the change of its joint entropy is totally controlled by the contraction or expansion of the phase space [6]. Later on, Liang and Kleeman showed that this property holds for deterministic systems of arbitrary dimensionality [20]. Liang and Kleeman [6] provided a very efficient heuristic argument to describe the decomposition of the various evolutionary mechanisms of information transfer in terms of the individual and joint time rates of entropy changes of X 1 , X 2 and (X 1 , X 2 ). Firstly, they computed dH 1 dt and dH 2 dt , where H i is the entropy of X i defined according to the marginal density, ρ i . Secondly, they employed the novel idea of frozen variables to analyze the individual time rates of entropy changes. When X i is fixed and X j evolves on its own in 2D systems, they found its temporal rate of change of entropy depends only on E( In the presence of interactions between X i and X j , they observed that Therefore, Liang and Kleeman [6] concluded that the difference between dH j dt and E ∂F j ∂x j should equal to the rate of entropy transfer from X i to X j . In the meantime, they denoted the rate of flow from X i to X j by T i→j (T stands for "transfer") and defined information flow/transfer as ρ(x j ,t) and i, j = 1, 2 with different i, j at the same time.

Discrete Mappings
Similarly, Liang and Kleeman [6] also gave the formalism about a system in the discrete mapping form. Considering a 2D transformation where x = (x 1 , x 2 ) ∈ Ω and Ω := Ω 1 × Ω 2 . The evolution of the density of Φ is driven by the Frobenius-Perron operator (F − P operator) P : L 1 (Ω) → L 1 (Ω) [18]. The entropy increases as where J −1 is the Jacobian matrix of Φ. When Φ j is invertible in 2D transformations, The entropy of X j increases as where ρ j is the marginal density of X j . When Φ j is noninvertible in 2D transformations, where P j is the F − P operator when x i is frozen as a parameter in P j . The entropy transferring from X i to X j is where i, j = 1, 2 with different i, j at the same time.

Continuous Flows
Firstly, we consider a three-dimensional (3D) continuous autonomous system, where F = (F 1 , F 2 , F 3 ) is a known flow vector. Similarly, the probability density ρ associated with Equation (8) satisfies the Liouville equation [18]: Analogous to the derivation in [6], firstly, multiplying by (1 + log ρ) for Equation (9), after some algebraic manipulations: Then, integrating for Equation (10), Assuming that ρ vanishes at the boundaries (the compact support assumption for ρ and the assumption is reasonable in real-world problems [6]), it is found that the time rate of the joint entropy change of X 1 , X 2 and X 3 , where E(∇ · F) = Ω ρ(∇ · F)dx 1 dx 2 dx 3 . As mentioned above, the time rate of change of H equals to the mathematical expectation of the divergence of the flow vector F. When we are interested in the entropy evolution of a component, x k in 3D systems, the marginal density is The evolution equation of ρ k is derived by taking the integral of Equation (9) with respect to x i and x j over the subspace Ω i × Ω j : ∂ρ k ∂t The third and fourth terms in Equation (9) have been integrated out with the compact support assumption for ρ. So the entropy for the component The Equation (11) states how H k evolves with time. The evolutionary mechanism of H k derives from two parts: One is from the evolution itself, dH * k dt ; another from the transfers of X i and X j according to the coupling in the joint density distribution ρ. From Section 2, we know that when X k evolves on its own, then Therefore, the rate of information flow/transfer from X i , X j to X k is where and i, j, k = 1, 2, 3 with different i, j, k at the same time. In particular, if F 1 = F 1 (x 1 ) has no dependence on x 2 , then T 2→1 = 0. There is no information transfer from random variable component X 2 to X 1 . This holds true with the transfers defined in LK2005 formalism. Obviously, in system (8) , when F 1 has no dependence on x 2 , x 3 , there should be no information transfer from X 2 , X 3 to X 1 , but there is possibility that the transfers in other directions may be nonzero when F 2 depends on x 1 , x 3 or F 3 depends on x 1 , x 2 . This is consistent with the information transfer defined in Equation (12). As a matter of fact, an important property of the transfer is given below.
Proof of Theorem 1. According to the formalism of information transfer for system (8), with the notation of F k = F k (x k ), It is worth noting that, while X k gains information from X i or X j , or X i and X j , X i or X j might have no dependence on X k in 3D systems. An important property about information transfer is its asymmetry among the components [1]. In addition, it is interesting to note that the formalism of 3D systems can be reduced to 2D cases under the condition that one variable does not depend on another variable. For example, If the evolution of X k is independent of X i , then In particular, when X k is independent of X i and X j , According to Theorem 1, the results are apparent. Furthermore, when X k depends on X i and X j , From the above derivations, we can see that our formalisms are further intensified by emphasizing the inherent relation with the formalisms in 2D systems. The information flows from two variables and the high order interactions between them to another variable are quantified by formula (12). These are generalized forms of the LK2005 formalism. In Section 4, we will validate the conclusions by the applications of all formulas in the Lorenz and Chua's systems. Moreover, when several variables are involved, the formalisms are capable to tackle information transfers of a multi-dimensional system. Combining the Liouville equation with Equation (3), dH dt = E(∇ · F) in n-dimensional situations, we can generalize the formalism to n-dimensional continuous and deterministic autonomous systems in the same way. For example, the transfer of information from components X 2 , X 3 , . . . , X n to X 1 is Hence, Theorem 1 can be generalized to multi-dimensional cases.

Discrete Mappings
For a 3D transformation Φ : , the evolution of its density is driven by the Frobenius-Perron operator (F − P operator) P : L 1 (Ω) → L 1 (Ω) [18]. Similar to the 2D case, after some efficient computations, the entropy transfer from X i , X j to X k in three-dimensional mappings has the following form: We also give a theorem for the discrete mappings and highlight the relationship between two-dimensional formalisms and generalized formalisms. The formalisms can be extended to high-dimensional situations as well. The detailed processes are demonstrated in Appendix A.

The Lorenz System
In this section, we propose an application to study the information flows about the Lorenz system [33]: where σ, r and b are parameters, x 1 , x 2 and x 3 are the system state variables, and t is time. A chaotic attractor of Lorenz system with σ = 10, r = 28, b = 8 3 is shown in Figure 1. Firstly, we need to obtain the joint probability density function ρ(x 1 , x 2 , x 3 ) of X to calculate information flows among the variables. For a deterministic system with known dynamics, the underlying evolution of the joint density ρ(x 1 , x 2 , x 3 ) can be obtained by solving the Liouville equation. Taking into account of the computational load, we estimate the joint density ρ(x 1 , x 2 , x 3 ) via numerical simulations. The steps are summarized as follows: • Initialize the joint density ρ(x 1 , x 2 , x 3 ) with a preset distribution ρ 0 , then generate an ensemble through drawing samples randomly according to the initial distribution ρ 0 . • Partition the sample space Ω into "bins".

•
Obtain an ensemble prediction for the Lorenz system at every time step.

•
Estimate the three-variable joint probability density function ρ via bin counting at every time step.
The Lorenz system is solved by applying a fourth order Runge-Kutta method with a time step ∆t = 0.01. According to Figure 1, the computation domain is restricted to Ω ≡ [− 30,30] × [− 30,30] × [0, 60], which includes the attractor of the Lorenz system. We discretize the sample space into 60 × 60 × 60 = 21, 600 bins to ensure covering the whole attractor and one draw per bin on average via making 21,600 random draws. Initially, we assume X is distributed as a Gaussian process N(u(t), Σ(t)), with a mean u and a covariance matrix Σ: Although we have used different parameters u and σ 2 d (d = 1, 2, 3) to compute information flows for the Lorenz system, the final results are the same and the trends stay invariant. The parameters u and σ 2 d can be adjusted for different experiments. Here we only show the results of one experiment with u d = 4 and σ 2 d = 4. The ensemble is developed by drawing sample randomly in the light of a pre-established distribution ρ 0 (x). We obtain an ensemble of X and estimate the three-variable joint probability density function ρ(x 1 , x 2 , x 3 , t) by the way of counting the bins, at every time step. As the equations are integrated forward in the Lorenz system, ρ can be estimated as a function of time and describe the statistics of the system. A detailed discussion on probability estimation through bin counting are referred to [20,25]. The sample data with initial value (1,1,1) and an estimated marginal density of x 1 , x 2 and x 3 are displayed in Figure 2. Through formula (12), the information transfer within three variables can be computed. There are nine transfer series in the Lorenz system, but here we mainly focus on the couple effect from two components to another component, that is, T i,j→k , i, j, k = 1, 2, 3 with different i, j, k at the same time. A nonzero T i,j→k means that X i and X j are causal to X k , and the value means how much uncertainty that X i and X j bring to X k . Among all the transfers, it is clearly shown that any two variables drive the other variable in the dynamics except the evolution of X 1 which only depends on X 2 . For the sake of revealing some underlying information in the chaotic dynamical system better, we also give information transfer among the components over space with a Gaussian distribution initialization and the averaged density over time via using the following formula: ∂x k dx i dx j , which characterizes the strength of information transfer at different planes of x = x i . That is to say, it demonstrates the information transfer of x j and x k to x i plane, whose relative values represent the magnitudes of information transfer. The calculation results are plotted in the left panel and the right one of Figure 3, respectively. According to the magnitude of parameters in the Lorenz system and the definition of rigorous 3D formalisms, the information transfer from X 1 , X 2 to X 3 is the smallest. The results are just as we expected, |T 1,2→3 | < |T 1,3→2 | < |T 2,3→1 |, as shown in the left panel of Figure 3. Meanwhile, we can get much information through numerical simulations. For example, the information transfer from X 2 and X 3 on X 1 is larger than that of X 1 and X 3 on X 2 in the Lorenz system, which is helpful for us to better analyze the system and the fields of interest. Only the absolute value of T measures the information transfer among the variables [23]. As the ensemble evolution is carried forth, any two variables aim to reduce the uncertainty of the other variable [24]; in other words, any two variables tend to stabilize the other variable. All information flows go to constants, which means that the system tends to be stable simultaneously. Comparing the left panel with the right one in Figure 3, we can find that not only the information flow from X 2 and X 3 to X 1 is the largest at different times, but also the total information transfer is the largest at x 1 plane, and the strength of information transfer obeys a distribution in each direction of x. Repeated experiments are found to be in line with the results no matter whatever the initialization is given. In particular, we compute the transfer, T 2→1 , then compare T 2→1 with the transfer, T 2,3→1 in Figure 4, as well as plot the transfers T 1→2 , T 3→2 and T 1,3→2 in Figure 5.  Since the evolution of X 1 is independent of X 3 and the evolution of X 2 depends on X 1 and X 3 in the Lorenz system, the transfer T 2→1 should be equal to T 2,3→1 and neither the transfer T 1→2 nor T 3→2 should not be equal to T 1,3→2 according to the derivations in Section 3.1. As expected, there is almost no difference between the two flows in Figure 4. The interpretation of the results is that X 3 is not causal to X 1 in the Lorenz system. The result agrees well with theoretical analysis, which also validates our formalisms. But the graphs T 1→2 and T 3→2 are quite different from the graph T 1,3→2 in Figure 5, as that both X 1 and X 3 are causal to X 2 in the Lorenz system. From Figures 4 and 5, we can find that the information flow T 2→1 is different from T 1→2 , as a property of asymmetry of the information transfer. There exists hidden sensitivity information in information transfer processes of high-dimensional dynamical systems: whether or not one variable brings more uncertainty to another variable. Comparing the magnitudes of three flows in Figure 5, we can say that X 3 is more sensitive to X 2 than X 1 to X 2 from the sensitivity analysis point of view. All the above differences are exactly the embodiment of the differences between the information flows in multi-dimensional dynamical systems and the LK2005 formalism. The proposed formalisms can be used to measure information transfers among the variables in dynamical systems and the numerical results can show how the measurement behavior with time, compared with the qualification of information transfer between two variables [4] and the transfer mutual information method [32]. For example, it can be quantified the influence that x 3 on the relationship between x 1 and x 2 using the transfer mutual information method in the Lorenz system. With our generalized formalisms, we can give quantitatively the influence from x 3 to the relationship between x 1 and x 2 as a dynamical process and other relationships (such as the asymmetrical influence between two variables) among the variables for analyzing the system better. To test the influence of error propagation on the measurement of information transfers, we use a different natural interval extension to compute information transfers according to the striking method [34]. In other words, we compute information transfers using formula (12) in the Lorenz system with the rewritten second equation, that is, rx 1 − x 1 x 3 is used to replace x 1 (r − x 3 ). For the Lorenz system, the results show that, the algorithm performs well (the relative error < 2%). All simulations are performed on a 64-bit Matlab R2016a environment. The physical consistency of the proposed approach in this paper can be explained as that a direction of the phase space is frozen in order to extract information transfers from the other two directions [3]. In addition, nonlinearity may lead a deterministic system to chaos, which causes the "spikes" in the right panel of Figure 3 and corresponds to intermittent switching in the chaotic dynamics. As the remarkable theory stated in [35], it indicates when the dynamics are about to switch lobes of the attractor in the Lorenz system.
Since Liouville equations and Frobenius-Perron analysis describe an ensemble of trajectories, we can use the generated formalisms of information flow as a sensitivity analysis index to perform dynamic sensitivity information analysis instead of the preceding widely used methods such as repeated calculation of principal component coefficients [36,37], construction of functional metamodels [31,38], calculation of moving average of the sensitivity index [39] and direct perturbation analysis of a dynamical system [40]. Using information flow to identify sensitive variables is directly based on the statistical perspective, which can improve numerical accuracy and efficiency while reduce the calculation load, compared with conventional dynamic sensitivity analysis methods. We cannot only quantify how much the uncertainty among variables of a system, but also understand how they influence system behavior, so it may be measured and used for prediction and control in realistic applications.
Furthermore, we use Equation (15) to compute information transfers, T yzw→x , T xzw→y , T xyw→z , and T xyz→w with the same strategy in the four-dimensional (4D) dynamical system: , whose results are shown in Figure 6.  The generalized formalisms are useful to deal with universal problems, which is not difficult to be applied to higher-dimensional cases.

The Chua's System
Since it is the first analog circuit to realize chaos in experiments, the initial Chua's system is a well-known dynamical model [41]. The Chua's system is described in reference [42] and there are many researches on its dynamical behavior [43,44]. Here we present an investigation of the information flows within the smooth Chua's system [45]: where p, q are parameters, x, y and z are state variables in R and t ∈ R + . When p = 11 and q = 14.87, a chaotic attractor of the Chua's system is shown in Figure 7.  Figure 7. The attractor of Chua's system with x(0) = −3, y(0) = 2, z(0) = 1. The former three trajectories are x, z-plane,x, y-plane and y, z-plane, respectively. The last trajectory is a 3D plot of x, y and z.
As mentioned before, using the same estimation procedures, we can obtain the density ρ(x, y, z) of R by counting the bins at each step. From Figure 7, the appropriate computation domain Ω ≡ [−10, 10] × [−10, 10] × [−10, 10] which includes an attractor of the Chua's system can be selected to estimate the three-variable joint probability density function. The following computation is demonstrated by applying a fourth order Runge-Kutta method. Similarly, we only show the results of one experiment after computing information flows multiple times by using different parameters. Suppose that R is distributed as a Gaussian process N(u(t), Σ(t)), with a mean u and a covariance matrix Σ in the initial state: Due to the additional fact that the smooth Chua's circuit has a highly non-coherent dynamics [46], we discretize the sample space into 200 × 200 × 200 = 8, 000, 000 bins to adequately understand the information transfer and the behavior of the system over time. A sample data and an estimation result of three marginal densities are shown in Figure 8, and we can find that the dynamical behaviors of the system are consistent with the results, such as symmetry. Using formula (12) to compute the information transfers within three variables of Chua's system. Firstly, we discuss the coupling effect from two components to the other component, the calculation results are demonstrated in Figure 9. Secondly, we compute the transfers, T y→z and T z→y , then compare T y→z with the transfer, T x,y→z and T z→y with T x,z→y in Figures 10 and 11, respectively. We also show the corresponding results of the strength of information transfer among the components with a Gauss distribution initialization and the averaged density over time in Figure 9. Since X causes Y but does not cause Z in the Chua's system, the numerical results of Figures 10 and 11 conform with the derivations of Equations (13) and (14) in Section 3.1. More specifically, there is almost no difference between the two flows in Figure 10, however, there exists large disparity between the two flows in Figure 11. The results also verify our formalisms. In addition, as shown in Figures 10 and 11, we can see that the information flow T y→z is different from T z→y due to the asymmetry of information transfer. All simulations are performed on a 64-bit Matlab R2016a environment. We are able to estimate that one variable makes another variable more uncertain or more predictable via the generalized formalisms. Besides, we can identify sensitive variables by computing information transfers among the variables in dynamical systems.
Compared with the Lorenz system, the Chua's system embodies in engineering systems besides that their discoveries were extraordinary and changed scientific thinking [46]. It can be used as another means to research, experiment and think about humanity, identity and art, etc. [47,48]. In studying visualization of the dynamics of Chua's circuit through computational models, the quantitative transformations of behavior are being taken into account [46]. The multi-dimensional formalisms of information flow enable us to improve our ability to estimate, predict, and control complex systems in many diverse fields. Furthermore, most existing approaches in control and synchronization of chaotic systems require adjusting the parameters of the model and estimating system parameters, which become an active area of research [49], and an additional benefit provided by the multi-dimensional formalisms of information flow is parameter estimation. We can compute information flows of the simulation model with different sets of parameters and do the same procedure for obtaining a group of feedback, then determine the optimal parameters that cater for the actual needs in order to put insight into complex behavior of models by comparing the change rates.

Conclusions
Based on the LK2005 formalism, we propose a rigorous and general formalism of the information transfer among multi-dimensional complex dynamical system components, for continuous flows and discrete mappings, respectively. Information transfers are quantified through entropy transfers from some components to another component, enabling us to better understand the physical mechanism underlying the superficial behavior and explore deeply hidden information in the evolution of multi-dimensional dynamical systems. When the generalized formalisms are reduced to 2D cases, the results are consistent with the LK2005 formalism. We mainly focus on the study of 3D systems and apply the formalisms to investigate information transfers for the Lorenz system and the Chua's system. In the above-mentioned two cases, we show that information flows of the whole evolution and the strength of information transfer at different planes, which implies that how uncertainty propagates and how dynamic essential information in the system transports. The results of experiments on the generalized formalisms conform with observations and empirical analysis in the literature, whose application may benefit many diverse fields. Compared with the qualification of information transfer between two variables [4] and the transfer mutual information method [32], the generalized formalisms are helpful for analyzing the relationships among the variables in dynamical systems and the research of complex systems. Moreover, since the formalism is built on the statistical nature of information, it has the potential to perform sensitivity analysis in multi-dimensional complex dynamical systems and advance our ability to estimate, predict and control these systems. In practice, for complex high-dimensional dynamical systems, it is not easy to give the dynamics analytically. Considering many critical data-driven problems are primed to take advantage of progress in the data-driven discovery of dynamics [35], we are developing a dynamic-free formulation to analyze information flows of multi-dimensional dynamical systems.
In the future, the formalism will be further generated to high-dimensional stochastic dynamical systems and time-delay systems. Meanwhile, future research should investigate how the information flow as a new indicator can be deployed in the frame of dynamic sensitivity analysis.
Author Contributions: Y.Y. proposed the original idea, implemented the experiments, analyzed the data and wrote the paper. X.D. contributed to the theoretical analysis and simulation designs and revised the manuscript. All authors read and approved the manuscript. and the Frobenius-Perron operator (F − P operator) P : L 1 (Ω) → L 1 (Ω) [18] which steers the evolution of its density. Loosely, given a density ρ = (x 1 , x 2 , x 3 ), P is defined as that where w represents any subset of Ω. When Φ is invertible, P can be expressed clearly as is the determinant of the Jacobian matrix of Φ. Similar to the two-dimensional case, the entropy increases concisely rewritten as Meantime, in the case when Φ k is invertible of 3D transformations, The entropy of X k increases as where ρ k is the marginal density of X k . When Φ k is noninvertible, where P k is the F − P operator when x i , x j is frozen as parameters in P k . It is easy to find that Equation (A5) reduces to Equation (A3) when Φ k is invertible. Therefore, the entropy transfers from X i , X j to X k can be unified into a form where i, j, k = 1, 2, 3 with different i, j, k at the same time.
Just as the former case with continuous variables, the information flow obtained by Equation (A6) has the following property: Theorem A1. If Φ k is independent of x i , x j in system (A1) with different i, j, k, then T i,j→k = 0.
The detailed proof of Theorem A1 is presented in Appendix A.1. Moreover, the formalism of 3D system can be reduced to the formalism in 2D cases with the previously mentioned conditions being satisfied. For example, when Φ k has no dependence on x i , In particular, when Φ k has no dependence on x i and x j , Furthermore, when Φ k has dependence on x i and x j , The above formalisms can also be generalized to n-dimensional systems by efficient processing of the relationship between the F − P operator w Pρ(x 1 , x 2 , . . . x n )dx 1 dx 2 . . . dx n = Φ −1 (w) ρ(x 1 , x 2 , . . . x n )dx 1 dx 2 . . . dx n and the entropy evolution at different time steps. For example, the transfer of entropy from X 2 , X 3 , . . . , X n to X 1 is  (Φ 1 (x 1 , x 2 , . . . , x n )) log P 1 ρ 1 (Φ 1 (x 1 , x 2 , . . . , x n )) · ρ(x 2 , x 3 , . . . , x n |x 1 )|J 1 |dx 1 dx 2 . . . dx n .
Here Ω 2...n is the simplified script of Ω 2 × Ω 3 × . . . × Ω n . Similar to the continuous cases, the generalized version of the property of Theorem A1 is also suitable for multi-dimensional discrete mappings.
Proof of Theorem A1. We only need to show that when Φ k is independent of x i , x j in 3D system, According to Equation (A4) and Equation (A5), we only need to prove According to the definition of the F − P operator and the condition that Φ k is independent of x i , x j at the same time,