Subject Areas: Biophysics
The concepts of stability, conservation law, symmetry are important in order to understand natural phenomena in physical, biological and engineering systems  -  . In our previous study, we proposed a conservative nonlinear dynamical model in 2n-variables which is primarily for nonlinear competitive interactions. The model provided a method for studying the relation between a conservation law and stability of a 2n-dimensional competitive system and restoration phenomena   .
Though symmetry and conservation laws have important information on nonlinear systems   , it is not possible to know at the outset what kind of conservation laws and nonlinearity should exist in complicated interacting systems, such as ecosystems of mammals   , microbes  , systems of cells and organs   . Nonlinear dynamical systems are characterized by self-interactions, dissipative structure, nonlinear coope- rative phenomena. They are highly evolved and interacting systems.
Complex systems in natural sciences from microbiology to ecology, economy and environmental sciences are regarded as realizations of nonlinear dynamical interactions. Nonlinear phenomena are difficult to handle because of their complex interactions and structures, self-organizations, spontaneous emergence of order, which exhibit no simple laws or orders. However, complex systems constitute cells, organs, and living animals which make stable systems. There may be some conservation laws corresponding to the stability observed in micro- biological systems. This is a motive that we applied nonlinear analyses based on Noether’s theorem to dynamics, and interesting results for complex phenomena are discussed   .
Our approach to nonlinear phenomena based on conservation laws is briefly reviewed for 2-variable model (), and the solution of 2-variable nonlinear equation is used to construct a solution of the 4-variable () conservative nonlinear equation. The approximate periodic amplitude and time of 4-variable nonlinear solutions are employed to explain stability of the net 4-variable system. The analysis of nonlinear systems with con- servation laws would help understand complex systems on a fundamental level.
The conventional nonlinear approach is reviewed in Section 1.1, and the atto-fox problem of Lotka-Volterra (LV) type nonlinear equation is explained in Section 1.2. Then, a new approach with conservation law and stability is explained in Section 1.3. In Section 2, the time-dependent perturbations and restoration from the perturbation are discussed with applications to the 10-year cycle of Canadian lynx and snowshoe hair and food-web of microbes in Okanagan Lake. In Section 3, the conservative nonlinear 2-variable model is extended to the 4-variable conservative model and properties of the 4-variable calculations are discussed. Conclusion and remarks are in Section 4.
1.1. Conventional Nonlinear Differential Equations for Interacting Systems
It is often explained that a model of interacting populations is written in general as,
where and are two autonomous differential equations which indicate that the time variable t does not appear explicitly in the functions and with the condition and . The meaning of functions depends on problems, where the nonlinear equations are applied, and this model is known as a general differential equation often called Kolmogorov’s predator-prey model  .
The simple and important applications in the field of ecology are those of Malthus (1959) for a population analysis, A. Lotka (1925)  and V. Volterra (1926)  for predator-prey differential equations known as Lotka-Volterra (LV) equation. The nonlinear equations have been applied to vast fields, such as biological pattern formations  -  , restoration phenomena from injuries  -  , epidemiology  , branching pro- cesses of growth and extinction of populations  , the 10-year cycle of Canadian lynx and snowshoe hare  -  , the food chain of microbes in a lake   , economy and Lanchester strategic management   , and so forth.
The mathematical expressions of and in differential Equation (1) have been proposed and physical meanings of, and stability as an equilibrium are investigated in terms of limit cycles and atractors of Lyapunov functions   . Lyapunov functions are scalar functions denoted as and, instead of and in (1), and nonlinear functions are generalized to include t with the condition and. There are two types of Lyapunov functions which are strict Lyapunov (having negative definite time derivative along trajectories) and non-strict Lyapunov functions (negative semi-definite time derivatives along trajectories)  . The notion indicates how strongly a nonlinear system approaches to its equilibrium state when. In other words, a nonlinear interacting system in the Lyapunov type corresponds to a dissipative system or a nonconservative system leading to an equilibrium after a long time. This property is different from 2n-variable conservative nonlinear differential equations   .
1.2. Note on the Atto-Fox (10−18-Fox) Problem on the Lotka-Volterra Equation
The typical Lotka-Volterra (LV) equation used in applications is given by the following simple nonlinear equation:
where x and y are the population number of prey and predator respectively; and represent variations of populations with respect to time. The constants, , are arbitrary nonlinear coefficients whose interpretations depend on problems at hand, and different nonlinear terms are readily included in the equation. It is not generally possible to obtain an analytical solution to the nonlinear differential equation, and a solution to (1.2) must be obtained by numerical computations with given coefficients as shown in the Figure 1 and Figure 2.
However, the nonlinear LV equation of the type (1.2) has an intrinsic problem known as atto-fox (10−18-fox) problem   which is readily explained in the Figure 2. As we know in a common sense, living animals
Figure 1. The values of nonlinear coefficients are, , , , with initial values, and.
Figure 2. The phase space solution.
interact with other animals through reasonably large integer numbers, however, the lower left of the Figure 2 shows that the small numbers such as at certain time t can interact with each other and continue to live, which is absurd and meaningless.
More specifically, the LV equation (2) is invariant with the scale change:
where a and b are arbitrary numbers. The equation for and with and is given by:
Hence, the solution in Figure 2 can be equivalently constructed in the phase space. Because of the arbitrariness of a and b, if we choose, for example, and, the solution can be expressed as and. The solution is obtained, and the atto-fox (10−18-fox) problem is proved. The LV equation of the type (1.2) has a serious potential problem when it is applied to real biological and ecological problems.
The problem of LV equations may be summarized as:
(a) the LV Equation (1.2) is too simple to apply to ecological and interactive systems.
(b) there are no restrictions on the nonlinear coefficients, , whose values would be only determined by experimental data through numerical simulations. However, the arbitrariness of nonlinear coefficients produces a serious problem (the atto-fox problem).
(c) density-independent external perturbations should be included in numerical simulations in order to include changes of climate, temperature, seasons and landscape and so forth. The external perturbations are independent of and.
Hence, it is necessary to find a new approach to resolve (a)-(c), which restricts values of nonlinear coef- ficients and numerically simulates real data with external perturbations.
1.3. A New Approach with Conservation Laws and Stability
It may seem natural to assume that a nonlinear system eventually dissipates energy (or an activation of the system) and arrives at an equilibrium state in the long run. However, the equilibrium state evolved as a dissipative process should be a maximum entropy state, and if the state needs to be activated, it needs some external inputs of energy. Hence, a nonlinear system demands certain external density-independent inputs or external perturbations which are not included in Lotka-Volterra and Kolmogorov type nonlinear Equations (1.1) and (1.2).
There are numerous examples of stable dynamical systems from microbiological systems such as DNA and RNA, amino acids and proteins, cells and organs, to macroscopic systems such as the 10-year cycle of Canadian lynx and snowshoe hair, wolves and caribous, creatures in marine coral reefs. The stability suggests that some conservation laws and symmetries are inherent from complex microbiological systems to macroscopic ecological systems. Hence, it could be possible to investigate conservation laws corresponding to a biological stability, which leads to a fundamental assumption for the model of conservative binary-coupled interactions   based on Noether’s theorem in dynamical systems  -  .
The meaning of a conservation law in biological complex systems may be very different from conservation laws in physics described by way of a Lagrangian or a Hamiltonian. It is difficult to use direct mechanical concepts and analogies such as potential and kinetic energy, but it should be useful to employ the concept of conservation law and symmetry in terms of Noether’s theorem in order to investigate conservation laws corresponding to biodiversity. Conservation laws and stability in biological systems could be employed to protect, understand and sustain the biosphere.
By employing the 2-variable conservative nonlinear model  , it is possible to numerically simulate biological data such as the 10-year cycle of Canadian lynx and snowshoe hair, microbe food-web in the Okanagan lake  . We review briefly the 2-variable conservative nonlinear interactions, because solutions of the 2-variable model is used to construct the solution of the 4-variable conservative nonlinear model.
The Lagrangian of 2-variable conservative nonlinear system is described with the following Lagrangian,
From (1.5), we get the following nonlinear differential equation,
The parameter is given by.
The conserved function, Y, of this system is given by,
The Psi-function, Y, may correspond to the Hamiltonian or the Lagrangian in a mechanical system, however, one should note that the Y-function does not have any velocity (time-derivative) terms, showing no individual mechanical energy terms. The terms of the Y-function look like potential energy terms in order to conserve the Y-function as a constant value, whereas the interpretation of the function depends on nonlinear systems.
The solutions and are rapidly changing in time as shown in Figure 3 (ND1) and Figure 4 (ND2), which are computed as examples with parameters and initial starting values explained in each figure captions. One should be careful that solutions are sensitive to the value of chosen nonlinear parameters and initial starting values. The nonlinear parameters and initial values cannot be changed arbitrary and should be chosen carefully by numerical simulations with external perturbations.
The solutions in the phase space, and, are shown in Figure 5, which are closed lines denoted as ND1 and ND2. The closed line in the phase space of solutions show that solutions respectively have definite periodic time. Though solutions () are rapidly changing in time, the conserved function, Y, of nonlinear solutions is constant in time as shown in Figure 6. When the solution of conservative nonlinear equation cannot exist or be found with given parameters, the Y-function diverges or is simply not constant. This is helpful to examine the accuracy of numerical solutions to conservative nonlinear differential equations.
In widely used nonlinear differential equations, the eight interaction terms on the right-hand sides of (1.6) and (1.7) require eight independent parameters (and are common factors and irrelevant). However, in case of the nonlinear interactions, regarding () as one parameter, we have only five independent parameters at
Figure 3. Numerical simulation of 2-variable ND1. The parameters of 2-variable ND1:, , , , , , initial values are
Figure 4. Numerical simulation of 2-variable ND1. The parameters of 2-variable ND2:, , , , , , , , initial values are,.
Figure 5. Numerical solutions of 2-variable ND1 and ND2.
the beginning in order to solve (1.6) and (1.7). This is helpful for numerical simulations of real data. The property of the 2n-variable conservative nonlinear equation is essentially different from properties of Lotka- Volterra and Kolmogorov equations.
We would like to answer some questions concerning difficulties on nonlinear differential equations in order to elucidate applicability of the 2n-variable conservative nonlinear differential equations.
(Question 1) There may be infinitely many periodic solutions in the nonlinear model that depend on the values of nonlinear coefficients, , and initial starting values,.
(Answer) This is restated as that solutions should be transformed in any form of solutions in the phase space, , resulting in the atto-fox problem or a physically meaningless problem. It may be correct for simple LV type nonlinear equations. However, one should be careful with the 2-variable conservative interaction model. The nonlinear coefficients and initial values in the 2-variable model are restricted as (1.6) and (1.7) by con- servation law, and positivity can be added as a constraint if necessary. Nonlinear coefficients and initial starting values are restricted in certain values and not entirely free adjustable parameters. The phase space solutions in Figure 5 of the 2-variable nonlinear model cannot be transformed continuously to
Figure 6. Numerical solutions of conserved quantity Y of 2-variable ND1 and ND2. Note that Y is constant in both cases.
those of the atto-fox (10−18-fox) solutions.
When initial starting values and nonlinear coefficients are inconsistent, one obtains no solutions and the Y-function diverges  . At present, we do not know of any theoretical relations among values of nonlinear coefficients and initial starting values in conservative nonlinear equations. The consistent values of initial starting values and nonlinear coefficients should be determined from numerical simulations by including external perturbations. The Y-function can be used to check accuracy of numerical solutions.
(Question 2) Selecting initial conditions of scarce prey and plentiful predators, unphysical solutions to nonlinear equations could occur. For example, may decrease even if it is zero when is positive and increasing, and may decrease or increase even if is negative. This means that non-existing prey can be eaten by predators, and predators can increase or decrease by eating or interacting with negative prey.
(Answer) This is a typical surmise derived from LV type nonlinear equations, but the assumption is not correct at all with the conservative nonlinear model. Because nonlinear coefficients are not entirely free parameters, it is not possible to find solutions by selecting any nonlinear coefficients and initial conditions of scarce prey and plentiful predators, and atto-fox type problems are restricted. A specific numerical example is shown in the paper  . The solutions to 2n-variable conservative nonlinear equations are different from those of conventional, nonconservative nonlinear equations. The (Question 2) is correct if one is solving simple non- linear equations with no constraints on nonlinear coefficients.
2. Time-Dependent Perturbations, Restoration from Perturbations
Lotka-Volterra and Kolmogorov equations are not sufficient to study nonlinear systems because they do not include density-independent external perturbations, such as climate and seasonal changes, pesticides and herbicides, hunting, forest fires and storms, etc. These external factors are density-independent and they are not expressible as (1.1). Hence, it is necessary to simulate density-independent and time-dependent external factors, and they are introduced as piecewise continuous, constant perturbations in the nonlinear equations discussed in detail in the papers   .
2.1. Piecewise Continuous External Perturbations
In order to investigate numerical responses of a system to external perturbations, piecewise continuous constants are introduced by using q-function such that
where represents a step function:
and coefficients are positive or negative constants to express strength of external perturbations.
An execution of an external perturbation to (), is shown in Figure 7 and Figure 8. The phase space solution denoted as St 1 (state 1) jumps to the solution Sp 1 (starting-point 1)-Ep 1 (end-point 1) when the external perturbation is switched on at the moment Sp 1, and goes back to the solution St 2 when the perturbation is switched off at Ep 1.
The instant that an external perturbation is added to conservative nonlinear equations, the solution jumps to possible another conserved system. When the perturbation was switched off (see, Ep 1 in Figure 8), the system recovered the original state or a possible stable state.
We gradually changed the strength and duration time of a negative external perturbation. The system advanced to unstable states and disintegrated, , corresponding to the extinction of species. The system cannot find solutions, , which is shown by a non-constancy of Y. Hence, there is a threshold for the negative external perturbation to each component, which makes the component amplitude negative and diverge, resulting in extinction of species  .
The strength and duration time of negative external perturbations which make the system unstable were evaluated numerically. The threshold of negative perturbation is related to the amplitude and periodic time of each component, , and if the strength and duration time of an external perturbation are not large compared to the amplitude and periodic time of species, the nonlinear interacting system converges to a possible stable solution. It should be possible to recover a complex interacting system with the assistance of external perturbations if we know the standard rhythm   or approximate amplitude and periodic time of component variables.
The threshold of strength and duration time of negative perturbations is important information to regulate nonlinear complex systems. The 2-variable conservative nonlinear model indicates that if we know the relation
Figure 7. The () phase-space transition with the negative perturbation as in (a). Solid line, St 1, (State 1) is the initial state, and St 2 (State 2) is the recovered state after the end of perturbation. Dashed line is phase-space during Sp 1 (starting-point 1)-Ep 1 (endpoint 1).
Figure 8. Conservation law Y of 2-variable ND with one perturbation. Y changed to by introducing perturbation. It recovers after Ep 1 (end-point 1).
between stability and rhythm of interacting systems and exert appropriate external perturbations at a correct time, it could be possible to save species from extinction.
2.2. The 10-Year Population Cycle of Canadian Lynx and Snowshoe Hair
As a specific example, the 2-variable nonlinear model with perturbations is applied to analyses of the 10-year cycles of Canadian Lynx and snowshoe hair  -  . The 2-variable conservative nonlinear model with external perturbations can successfully simulate the reported data of the Lynx-hair, and we found that there is a stable cycle, the standard rhythm, inferred from the analysis of conservative nonlinear interactions. It is possible to understand the stability and standard rhythm of the 10-year cycles of Canadian Lynx and snowshoe hair in terms of a conservation law of a prey-predator system  .
The standard rhythm is discussed in the paper  , which explains stable population cycles, such as the 10- year cycle of Canadian lynx and snowshoe hair and microbes’ food-web in Okanagan Lake. The population cycles simulated by the current conservative nonlinear interaction are stable and show recovery and restoration from external perturbations as long as the environmental systems for living and interacting animals are not seriously damaged. Figure 7 is a numerical example of recovery and restoration phenomena.
The current conservative nonlinear model is applied to the data of food-web in Okanagan Lake   ; mysis introduction to the Lake is known as an effective method to enhance ecological interactions among microbes and other creatures so as to increase fisheries’ production. The maximum and minimum density- amplitude of microbes are also well simulated with the 2-variable conservative nonlinear model, and the results suggest that stability and periodic restoration phenomena in ecosystems may be a manifestation of a conser- vation law with external perturbations.
The important property of the 2-variable nonlinear model with conservation law is that coupled systems can have certain recovering strength to external perturbations. As a predator needs a prey for its food, a prey needs a predator for the control and conservation of their own species. The conservation law and the standard rhythm of species are considered to be naturally constructed by species to live for a long time in nature. Therefore, the rhythm and stability of population cycle as a whole in an ecosystem would be interpreted as a manifestation of the survival of the fittest adapted to the balance in an ecosystem.
The stability and conservation law are constructed by species in mutual dependency or cooperation to survive for long-time periods in severe environmental conditions. The standard rhythm should be regarded as the result of strategy for species to live in nature. Whatever roles they have to play, the species that can fit and balance with other creatures can survive in nature. A strong predator cannot even survive if it ignores the law of the standard rhythm and conservation law of nature constructed by other members and the environment. The conservative BCF nonlinear model with perturbations enables one to examine a fundamental law of nature in a form of differential equations. We hope that this study will help understand both activities of animals and humans in natural life.
3. The 4-Variable Dynamical Model
Based on the reviews and discussions, we extend the 2-variable conservative model to the 4-variable model for applications to more realistic interacting systems. Several examples can be found in microbiology and medical fields   . Statistical methods are normally used to investigate correlations among physical data (mor- phogens, responses of cells, concentrations, populations, etc.). Dynamical analysis can be possible by applying the 4-variable conservative nonlinear differential equation.
A morphogen is considered a substance governing the pattern of tissue development in the process of mor- phogenesis, or a signaling molecule that acts directly on cells to produce specific cellular responses depending on its concentration. When morphogens are identified, a dynamical method based on the 2-variable or 4-variable conservative nonlinear dynamical analyses could be possible in microbiology.
However, the 4-variable nonlinear differential equations have too many nonlinear parameters to be set at the beginning, and they are difficult to solve and apply to numerical simulations. Therefore, we discussed a strategy to construct solutions to 4-variable, 6-variable, … , conservative nonlinear differential equations. The construc- tion of 2n-variable solution is based on the solution of 2-variable solutions.
The 4-variable conservative nonlinear model for complex systems is described as:
where coupling constants, , and, and, represent interaction strength among species (). The numbering of nonlinear coefficients, , is chosen for our convenience by starting from a lagrangian to generate (3.1)-(3.4). It is ex- plained in the Appendix.
The external perturbations are expressed as arbitrary, piecewise continuous constants, , which represent strengths of external effects (such as, temperature changes, reductions of species by huntings or fishings, etc., as cases may be). They are pulse-like constants expressed by step-functions.
The 4-variable conserved nonlinear equation is derived from Euler-Lagrange equation in Lagrangian dynamics. The Noether’s theorem in dynamics is employed to obtain the conservation law in complicated non- linear systems. Hence, coefficients of the right-hand side of (3.1)-(3.4) are restricted such that the system has conserved quantities.
In conventional 4-variable nonlinear differential equations of the type (3.1)-(3.4), nine terms in the right-hand side are supposed to be independent and adjusted freely. Hence, the number of nonlinear coefficients basically adds up to 40 parameters (and 4 parameters,). In the current nonlinear equation, coef- ficients of the type, (), , can be handled as one parameter, and the number of independent para- meters adds up to 22. There are still many free parameters, however, admissible solutions would be confined further when initial starting values and positivity, , are included. These conditions and the conservation law, , are helpful for numerical simulations.
In order to make the 4-variable nonlinear equation useful for numerical simulations, we propose a strategy to construct a 4-variable solution by employing 2-variable solutions.
1) One should construct two sets of 2-variable numerical solutions independently by dividing interactions of 4 variables () as 2-variables () + 2-variables (). The two sets of the 2-variable conservative nonlinear solutions, () and (), should be prepared from respective data. External perturbations could be included to create reasonable numerical simulations in the 2-variable form   .
2) The 2-variable numerical simulations, (), are regarded as solutions for (3.1) and (3.2) by assuming coefficients related to and as zero. Similarly, the 2-variable solutions, (), are also regarded as solutions for (3.3) and (3.4) by assuming coefficients related to and as zero. Hence, the value of coefficients in the 4-variable solution (3.1)-(3.4) are almost fixed at the beginning.
3) The next step is to weakly couple the set of 2-variable solutions by using the 4-variable solution (3.1)-(3.4). The weak coupling of () and () can be performed by changing cross-coupling coefficients (assumed zero) slowly from 0 to (). We have four coupling coefficients, , and four coefficients from linear terms, , which can be regarded as four para- meters. Hence, in the 4-variable conservative nonlinear equation, we are supposed to change eight parameters to numerically simulate data in question. By referring to combinations of data: (), (), () and (), reasonable eight parameters should be searched in numerical simulations by trial and error.
The coupling constants, , and four coefficients of linear terms characteristically change amplitude and wave-shape. Even if coupling constants are small, one would encounter divergences when nonlinear parameters and initial conditions are not compatible, which can be specifically perceived by checking the value of Y function. The Y-function is very useful to check the existence of solutions.
Following theoretical discussions of invariant variational method in dynamics, the Y-function for the 4- variable conserved nonlinear dynamical system is expressed as,
The physical quantities () change in time, but the sum of the combination given by (3.5) is constant in time. The strategy to find a solution, (1)-(3), and experiences with nonlinear differential equations help obtain a numerical solution to the system of nonlinear differential equations. A set of solution to (3.1)-(3.4) and the value of Y-function are shown in Figures 9(a)-(d).
The 4-variables () are time-changing as in Figure 9(a), but Y-function is constant in time as shown in Figure 9(b), which confirms numerical solutions are accurate. The phase spaces of solutions for and are shown in Figure 9(c) and Figure 9(d), respectively. One can see that the phase space of solutions is confined in each volume of phase space. The other phase spaces of solutions are also shown in Figures 10(a)-(d).
The Results of the 4-Variable Conservative Nonlinear Model
The 2-variable nonlinear solution shows a recurrence time, T, defined by, and the simple closed line as Figure 5 indicates that the recurrence time of solution (periodicity) is finite. However, trajectories of solutions in the 4-variable conservative nonlinear model in Figure 9(c), Figure 9(d) and Figures 10(a)-(d) are not closed lines, but they enclose a finite area in each phase space. The solutions, , oscillate within a finite, closed surface of phase space and exhibit approximate periodic time, () produced by nonlinear interactions.
It is well known in linear differential equations that the net response (amplitude) y at a given place and time is
the sum of independent responses ():, which is known as the superposition principle.
The superposition principle does not hold in nonlinear systems. However, we examined the net amplitude
defined by, because the net amplitude may have important information for the
In our conservative nonlinear interactions, we found that component solutions approximately produced a periodic shape of the net amplitude. The shape of the net amplitude is very complicated but oscillates at an approximate pitch like a resonance produced by a tuning fork on resonance box.
The approximate periodic times of each variable in Figure 9(a) are respectively within (). The net periodic time is found from the resulting periodic shape of the am- plitude, y, shown in Figure 11. The approximate net periodic time of the system is in the current example as observed in Figure 11. One can examine that the almost similar shape of wave pattern within is repeated, which is also checked by changing time-span longer in the numerical calculations.
Figure 9. The solutions of 4-ND system and phase-space with respect to time. Figure 9(a) represents solutions and Figure 9(b) is the conservation law Y with respect to time. Figure 9(c) and Figure 9(d) are phase-space of solutions of and , and, respectively.
Figure 10. The phase-space of 4-variable ND model for another combination of solutions. Figure 10(a) represents phase-space of and, Figure 10(b) is a phase space of and. Figure 10(c) and Figure 10(d) are phase-space of solutions of and, and, respectively.
Figure 11. The sum of amplitude:. The approximate periodicity can be observed from the data.
Hence, conservative nonlinear interactions of component variables lead to a longer periodic time and a net stable system against external perturbations. In ecosystems, it could be interpreted so that living animals tend to have interactions with other animals as much as possible for their species to live long. The approximate periodic time of components and the net periodic time in nonlinear interactions may be one of important properties to understand stability of nonlinear interacting systems.
Our strategy of constructing solutions of nonlinear differential equations is to find a class of solutions having approximate periodic times. This may be true, and we investigated the class of solutions because the amplitude, y, and net periodic time, , are important physical observables to understand nonlinear dynamical systems. Therefore, it could be understood that the stable periodicity, the standard rhythm, discussed in papers   is the consequence of an approximate periodic phenomena of conservative nonlinear interactions.
The ten-year cycles of Canadian lynx and hair, interactions of microbes and other stable cycles are considered as examples of the approximate periodicity of nonlinear interacting systems. Hence, a possible condition to sustain stability of complex symbiotic systems can be examined by applying a conservative nonlinear dif- ferential equation.
The remarkable properties derived from conservative nonlinear interactions are summarized as follows:
1) The net periodic oscillation and time can be approximately defined in the conservative nonlinear interactions, and the net periodic time is longer than the periodic time of component variables:
2) The 4-variable system becomes more stable than the 2-variable system against negative external per- turbations. Nonlinear interactions of diverse component variables seem to maintain stability of the net system against negative external perturbations.
These results are different from those discussed in dissipative, nonconservative nonlinear interactions in many literatures which discuss limit cycles and attractors. The conservation law and durability of a conservative nonlinear system against external perturbations could be a key to understand stability of complex systems. The pattern simulations combined with external perturbations are possible for conserved nonlinear interactions, and it may be helpful to understand characteristic periodicities and responses of complex interacting systems.
In physics and engineering applications, it is standard to analyze physical quantities caused by heat conduction, electromagnetic and mechanical problems by way of the superposition principle   . The superposition principle does not hold in nonlinear systems, and nonlinear interactions generate chaos in dynamical systems  . However, by applying the 2-variable solutions, we extracted periodic solutions to complex systems des- cribed with conservative nonlinear interactions. The superposition of component solutions exhibited an approximate net periodic amplitude and time, which is important to examine complex systems.
We reviewed characteristic properties for the 2-variable conservative nonlinear differential equation which generalizes Lotka-Volterra type prey-predator, competitive interactions, and explained that the atto-fox problem found in a simple Lotka-Volterra equation is not intrinsic to nonlinear differential equations. The 2-variable and 4-variable conservative nonlinear equations with external perturbations are useful to simulate real data num- erically.
The conservative nonlinear interaction with external perturbations numerically suggested the existence of the standard rhythm (the 10-year population cycle) of Canadian lynx and snowshoe hare. The analysis of the 2- variable conserved nonlinear equation indicates that nonlinear interactions could be a manifestation of strategy to live and thrive in nature   . It is also confirmed in the 4-variable nonlinear model that the approximate periodic time of the net system is longer than each periodic time of component variables. A system in diverse interactions could be more stable than a simple system of 2-body interaction.
It is naturally observed that many complex systems from microscopic to macroscopic world seem to maintain stable structure, recoveries and self-organizations  -  . The conservation laws and symmetries could be a key to understand stability of complex systems. We employed the concept of conservation law or symmetry of Noether’s theorem in dynamics. A conservative nonlinear interaction is defined with nonlinear interactions generated by a function Y which may be considered as Hamiltonian or Lagrangian in dynamical systems. The Y-function can be used to check the accuracy of a numerical solution to the corresponding nonlinear differential equation   .
The solution of the 2-variable model is applied to construct solutions of the 4-variable conservative nonlinear equations. It is difficult to find an appropriate numerical solution for the 4-variable nonlinear differential equation, which tends to make numerical analyses hard to carry out. However, we showed a strategy to find a solution by employing the 2-variable solutions, which helps construct solutions and accomplish numerical analyses ri- gorously by employing Y-function. This approach can be extended to the 6 and 8-variable conservative non- linear equations.
The conservative nonlinear model has not been actively applied yet, though it revealed certain properties of the 10-year cycles of lynx and hair, microbe interactions. We are searching for ecological data or medical applications as microscopic interacting systems, such as diabetes or other chronic diseases exerted by multiple interacting factors. Collaborations to analyze big data with the current model are anticipated.
The nonlinear differential equation approach is useful to understand physical phenomena. We hope that the 4-variable conservative nonlinear equation may find useful and practical applications from microscopic complex systems to macroscopic ecosystems.
Nonlinear Parameters in 4-Variable Lagrangian
The numbering of nonlinear coefficients, in the 4-variable Lagrangian in the Section 3, is shown in order to avoid complications. The nonlinear coefficients are numbered as:
and the nonlinear Equations (3.1)-(3.4), are derived from the lagrangian (5.2) by way of Euler-Lagrange equation:
The numbering of nonlinear coefficients of parentheses, , can be respectively handled as one independent parameter. The coefficients, , are related to overall factors and become trivial in numerical calculations for (3.1)-(3.4). Hence, the number of nonlinear coefficients for the 4-variable con- servative nonlinear equation is 22. It should be noticed that if one directly starts from 4-variable nonlinear nonconservative differential equations in the same type, one is naturally directed to introduce 40 nonlinear coefficients.
The number of 4-variable nonlinear parameters is still large, but if 2-variable numerical solutions are used to construct the 4-variable solution, one is supposed to control only 8 parameters. The method of solution is ex- plained in the Section 3.1.
Submit or recommend next manuscript to OALib Journal and we will provide best service for you:
Publication frequency: Monthly
9 subject areas of science, technology and medicine
Fair and rigorous peer-review system
Fast publication process
Article promotion in various social networking sites (LinkedIn, Facebook, Twitter, etc.)
Maximum dissemination of your research work
Submit Your Paper Online: Click Here to Submit
Contact Us: email@example.com
 Ramos, A.F., Innocentini, G.C.P., Forger, F.M. and Hornos, J.E.M. (2010) Symmetry in Biology: From Genetic Code to Stochastic Gene Regulation. IET Systems Biology, 4, 311-329.
 Allesina, S. and Tang, S. (2012) Stability Criteria for Complex Ecosystems. Nature, 483, 205-208.
 Thébault, E. and Fontaine, C. (2010) Stability of Ecological Communities and the Architecture of Mutualistic and Trophic Networks. Science, 329, 853-856.
 Daafouz, J., Riedinger, P. and Iung, C. (2002) Stability Analysis and Control Synthesis for Switched Systems: A Switched Lyapunov Function Approach. IEEE Transactions on Automatic Control, 47, 1883-1887.
 Uechi, L. and Akutsu, T. (2012) Conservation Laws and Symmetries in Competitive Systems. Progress of Theoretical Physics Supplement, 194, 210-222.
 Naz, R (2010) Conservation Laws for a Complexly Coupled KDV System, Coupled Burgers’ System and Drinfeld-Sokolov-Wilson System via Multiplier Approach. Communications in Nonlinear Science and Numerical Simulation, 15, 1177-1182.
 Zhang, Y. and Zhou, Y. (2013) Symmetries and Conserved Quantities for Fractional Action-Like Pfaffian Variational Problems. Nonlinear Dynamics, 73, 783-793.
 Stenseth, N.C., Falck, W., Chan, K.S., Bjornstad, O.N., O’Donoghue, M., Tong, H., Boonstra, R., Boutin, S., Krebs, C., J. and Yoccoz, N.G. (1998) From Patterns to Processes: Phase and Density Dependencies in the Canadian Lynx Cycle. Proceedings of the National Academy of Sciences, 95, 15430-15435.
 Schindler, D.E., Carter, J.L., Francis, T.B., Lisi, P.J., Askey, P.J. and Sebastian, D.C. (2012) Mysis in the Okanagan Lake Food Web: A Time-Series Analysis of Interaction Strengths in an Invaded Plankton Community. Aquatic Ecology, 46, 1-13.
 Chyu, J., Fonarow, G.C., Tseng, C.H. and Horwich, T.B. (2014) Four-Variable Risk Model in Men and Women with Heart Failure. Circulation: Heart Failure, 7, 88-95.
 Nair, S., Mishra, V., Hayden, K., Lisboa, P.J.G., Pandya, B., Vinjamuri, S., Hardy, K.J. and Wilding, J.P.H. (2011) The Four-Variable Modification of Diet in Renal Disease Formula Underestimates Glomerular Filtration Rate in Obese Type 2 Diabetic Individuals with Chronic Kidney Disease. Diabetologia, 54, 1304-1307.
 Gierer, A. and Meinhardt, H. (1972) A Theory of Biological Pattern Formation. Biological Cybernetics, 12, 30-39.
 Turing, A.M. (1952) The Chemical Basis of Morphogenesis. Philosophical Transactions of the Royal Society B, 237, 37-72.
 Kondo, S. and Asai, R. (1995) A Reaction—Diffusion Wave on the Skin of the Marine Angelfish Pomacanthus. Nature, 376, 765-768.
 Lengyel, I. and Epstein, I.R. (1991) Modeling of Turing Structures in the Chlorite-Iodide-Malonic Acid-Starch Reaction System. Science, 251, 650-652.
 Kumar, N. and Horsthemke, W. (2011) Effects of Cross Diffusion on Turing Bifurcations in Two-Species Reaction-Transport Systems. Physical Review E, 83, Article ID: 036105.
 Haccou, P., Jagers, P. and Vatutin, V.A. (2005) Branching Processes. Cambridge University Press, Cambridge.
 Stenseth, N.C., Falck, W., Chan, K.S., Bjornstad, O.N., M. Donoghue, O., Tong, H., Boonstra, R., Boutin, S., Krebs, C.J. and Yoccoz. N.G. (1998) From Patterns to Processes: Phase and Density Dependencies in the Canadian Lynx Cycle. Proceedings of the National Academy of Sciences of the United States of America, 95, 15430-15435.
 Elton, C. and Nicholson, M. (1942) The Ten-Year Cycle in Numbers of the Lynx in Canada. Journal of Animal Ecology, 11, 215-244.
 Stenseth, N.C., Falck, W., Bjørnstad, O.N. and Krebs, C.J. (1997) Population Regulation in Snowshoe Hare and Canadian Lynx: Asymmetric Food Web Configurations between Hare and Lynx. Proceedings of the National Academy of Sciences of the United States of America, 94, 5147-5152.
 Krebs, C.J., Boonstra, R., Boutin, S. and Sinclair, A.R.E. (2001) What Drives the 10-Year Cycle of Snowshoe Hares? BioScience, 51, 25-35.
 Basille, M., Fortin, D., Dussault, C., Ouellet, J.P. and Courtois, R. (2013) Ecologically Based Definition of Seasons Clarifies Predator—Prey Interactions. Ecography, 36, 220-229.
 Blasius, B., Huppert, A. and Stone, L. (1999) Complex Dynamics and Phase Synchronization in Spatially Extended Ecological Systems. Nature, 399, 354-359.
 Maquet, J., Letellier, C. and Aguirre, L.A. (2007) Global Models from the Canadian Lynx Cycles as a Direct Evidence for Chaos in Real Ecosystems. Journal of Mathematical Biology, 55, 21-39.
 Schindler, D.E., Carter, J.L., Francis, T.B., Lisi, P.J., Askey, P.J. and Sebastian, D.C. (2012) Mysis in the Okanagan Lake Food Web: A Time-Series Analysis of Interaction Strengths in an Invaded Plankton Community. Aquatic Ecology, 46, 1-13.
 Gazi, N.H. (2012) Dynamics of a Marine Plankton System: Diffusive Instability and Pattern Formation. Applied Mathematics and Computation, 218, 8895-8905.
 Taniguchi, M., Bando, M. and Nakayama, A. (2008) Business Cycle and Conserved Quantity in Economics. Journal of the Physical Society of Japan, 77, Article ID: 114001.
 González, E. and Villena. M. (2011) Spatial Lanchester Models. European Journal of Operational Research, 210, 706-715.
 Mazenc, F. and Malisoff, M. (2010) Strict Lyapunov Function Constructions under LaSalle Conditions with an Application to Lotka-Volterra Systems. IEEE Transactions on Automatic Control, 55, 841-854.
 Smith, H.L. and Waltman, P. (1995) The Theory of the Chemostat. Cambridge University Press, Cambridge.
 Mollison, D. (1991) Dependence of Epidemic and Population Velocities on Basic Parameters. Mathematical Biosciences, 107, 255-287.
 Kloeden, P.E. and Pötzsche, C. (2010) Dynamics of Modified Predator-Prey Models. International Journal of Bifurcation and Chaos, 20, 2657-2669.
 Sieniutycz, S. (1994) Conservation Laws in Variational Thermo-Hydrodynamics. Kluwer Academic Pub., Norwell.
 Ott, E. (2002) Chaos in Dynamical Systems. Cambridge University Press, Cambridge.